第四章空间轴对称问题精选PPT.ppt

《第四章空间轴对称问题精选PPT.ppt》由会员分享，可在线阅读，更多相关《第四章空间轴对称问题精选PPT.ppt（60页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章空间轴对称问题第1页，此课件共60页哦 由于空间轴对称问题的几何形状，约束情况及轴对称弹性体所受的外载荷都对称于z轴，如图，故这种弹性体内各点的各项应力分量、应变分量和位移分量也都对称于z轴，而与环向坐标 无关，所谓各项应力分量、应变分量和位移分量都与 坐标无关，其含义是，在任何一个过z轴的子午面上的位移、应变和应力的分布规律都相同。4-1 4-1 轴对称问题的弹性力学基本方程轴对称问题的弹性力学基本方程第2页，此课件共60页哦 由于各分量都对称于z轴，与 无关，因此弹性体内各点只可能存在着径向位移u和轴向位移w（此时，u,w只是r,z的函数），而环向位移v=0。则必有剪应变 ，否则这种

2、轴对称问题的弹性体将不能保持轴对称状态而发生歪扭，这在实际中是不可能出现的。根据轴对称问题弹性体的上述特点，即v=0，。其平衡微分方程为：(4-1)第3页，此课件共60页哦几何方程为物理方程式中轴对称问题弹性体的弹性矩阵；轴对称问题弹性体的应变列阵。(4-2)(4-3)第4页，此课件共60页哦 根据轴对称问题的特点，我们只须考查轴对称结构的任意子午面上各单元的各应力分量和各结点的位移分量。在轴对称结构的任意子午面上（即rz平面）任取一个三角形单元i,j,m，基本未知量仍然取结点位移，单元的结点位移可用列阵表示为(4-4)4-2 4-2 三角形截面环单元三角形截面环单元图5-1 轴对称结构第5页

3、，此课件共60页哦 仿照平面问题，取线性位移模式 类似于平面三角形单元的推导，可得 (b)(a)其中形函数(c)第6页，此课件共60页哦而 (4-5)(4-6)第7页，此课件共60页哦(b)式也可写成矩阵形式 (d)将(b)式代入几何方程(4-2)式，得到单元体内的应变，即 (e)第8页，此课件共60页哦其中(e)式仍然还可以简化成其中(f)(4-7)第9页，此课件共60页哦 由此可见，单元中的应变分量 都是常量，但是环向正应变 不是常量，它与 中的r有关。单元的应力分量仍可表示为 (h)第10页，此课件共60页哦其中而(i)(4-8)第11页，此课件共60页哦 显然，只有应力分量 在单元中为

4、常量外，其余三个正应力在单元中都不是常量。在实用上，为了简化计算和消除对称轴上由于r=0所引起的麻烦，常把各个单元中的r及z近似地当作常量，并且分别等于各单元形心的坐标，即于是(f)式成为 这样就可把各单元近似地当作常应变单元。将(j)、(k)式代入(4-7)和(4-8)式求得的是单元形心处应变和应力的近似值。(j)第12页，此课件共60页哦 现在，再运用虚功原理求导轴对称结构上任意单元的刚度矩阵 。由虚功原理知：三角形断面的环形单元体积所吸收的虚变形能应等于单元结点力所做的虚功：假设单元的虚位移为 则单元的虚应变为4-3 4-3 单元刚度矩阵单元刚度矩阵(a)(b)第13页，此课件共60页哦

5、将上式代入(a)式，并注意到，得由于虚位移是任意的，所以有(d)第14页，此课件共60页哦上式右边与单元结点位移列阵 相乘的矩阵便是单元刚度矩阵它也可以写成下列分块形式其中的子矩阵为(4-11)(4-9)(4-10)第15页，此课件共60页哦 由于在轴对称问题的矩阵 中出现坐标r、z，所以(4-11)式的积分运算比平面问题要复杂得多。现在仍取单元形心的坐标替代 矩阵中的坐标r、z作为一次近似，得到一个近似的单元刚度矩阵。此时，(4-11)式成为 (4-12)第16页，此课件共60页哦 对于整体刚度矩阵，如果弹性体被划分为 个单元和 n 个结点，于是就可得到 个型如(d)式的方程组。与平面问题的

6、情况完全相类似的处理，把各单元的 、等都加以扩大到整个结构的自由度的维数，然后叠加得到(f)引进记号：载荷列阵 (4-13)第17页，此课件共60页哦 整体刚度矩阵于是(f)式便可以写成与平面问题相同的标准形式这就是求解结点位移的平衡方程组。(4-14)(g)整体刚度矩阵也可以写成分块形式 第18页，此课件共60页哦第19页，此课件共60页哦其中子矩阵为 和平面问题一样，整体刚度矩阵 是对称的带状稀疏阵，在消除刚体位移后，它是正定的。第20页，此课件共60页哦 (g)式右边的载荷列阵展开的形式为其中(b)4-4 4-4 等效结点力的计算，载荷列阵等效结点力的计算，载荷列阵(a)与平面问题一样，

7、等效结点力也是由作用在环形单元上的集中力、表面力和体积力分别移置到结点上而得到的。移置的原则也是根据这些力和等效结点力在任意虚位移上所作的虚功相等，即第21页，此课件共60页哦式中 为集中载荷 作用点的径向坐标。将上节(b)式代入上式并考虑到 ，上式可以化成式中右边第一项是环形单元上的集中力 移置到结点的等效结点力，第二项是环形单元上表面力 的等效结点力，第三项是环形单元体积力 的等效结点力。(c)第22页，此课件共60页哦 采用平面问题中相同的符号：集中力的等效结点力 表面力的等效结点力 体积力的等效结点力(4-19)(4-17)(4-18)第23页，此课件共60页哦于是(c)式可以改写成再

8、将上式代入(a)式，等效载荷列阵可写成 将(4-18)、(4-19)式和(2-39)、(2-40)式比较可见，在轴对称情况下积分号后的被积函数比平面问题的多一个变量r，所以虽然也是采用线性位移模式，但是不能象平面问题那样利用刚体的静力等效原则求得结点等效力。当体积力或表面力可表示为坐标r和z的多项式时，不难利用(2-25)或(2-26)式精确积分得到等效结点力。(d)(e)第24页，此课件共60页哦 1.体积力 (1)自重。在此情况下 ；其中 为密度。于是单元的自重移置到结点i,j,m上的等效结点力为 由类似等参元的坐标变换式，将r写成(f)(f)第25页，此课件共60页哦这样就得到代入(f)

9、式即得 如果单元离开对称轴较远 ，可以认为将的自重移置到每个结点上。(4-20)第26页，此课件共60页哦 (2)离心力。在此情况下 ，其中 为角速度。于是单元的离心力移置到结点i,j,m上的等效结点力注意到(f)式，得代入(g)式即得(4-21)(g)第27页，此课件共60页哦 2.表面力 设rz平面上单元ijm的ij边上受有线性分布的径向表面力。在此情况下有 ，于是结点i的等效结点力为(h)第28页，此课件共60页哦注意到在ij边上 ，于是积分代入(h)式得到 (4-22)第29页，此课件共60页哦 经过类似推导可得移置到结点j和m上的等效结点力(4-23)(4-24)第30页，此课件共6

10、0页哦 两种特殊情况 (1)只有当单元离开对称轴较远时，才可以认为 与 大致相等，此时可由(4-25)式得出简单的结果；即将面力的 移置到结点i，移置到结点j。(4-25)第31页，此课件共60页哦 (2)显然，只有当单元离开对称轴较远时，才可以认为 、大致相等，则由上式得出简单的结果；即将面力的 移置到结点i，移置到结点j。需要注意，在轴对称问题中的结点力实际上是整个结圆上的力，这是与平面问题不同的。(4-26)第32页，此课件共60页哦 3考虑温度改变的影响，由于应力应变关系包含温度应变，于是公式(e)还应加上热载荷单元e中结点i上的温度改变引起的结点力为式中 为温度改变引起的应变，将(4

11、-7)、(4-3)式中的 矩阵和上式代入(6-27)式得到(4-28)(4-27)第33页，此课件共60页哦上式 中的两个积分通常可用数值积分法求积，或简单地取近似表达式 将以上两式代入(4-28)式，并注意到4-2中的(i)式，使得温度改变引起的等效结点力的近似表达式(4-29)第34页，此课件共60页哦 在4-3中我们导得单元刚度矩阵 ，它可以划分为3个22的子矩阵，即 由于子矩阵 、与坐标r、z有关，上式的积分不可能象平面问题那样简单地进行。为了避免复杂的积分运算，在4-3中用单元形心的坐标值 代替其中的坐标r和z，导得近似刚度矩阵(4-12)式。现在我们来推导精确的刚度矩阵。为此，把子

12、矩阵分成两部分(a)4-5 4-5 精确刚度矩阵的计算精确刚度矩阵的计算第35页，此课件共60页哦其中 是象4-3中那样，用单元形心坐标 代入 后得到的；第36页，此课件共60页哦而 是它的变化部分。由(4-7)式可知(b)第37页，此课件共60页哦把(a)式代入(4-11)式，有注意到矩阵 、的元素都是常量，它们可以提到积分号的外面，而且(c)第38页，此课件共60页哦故有于是(c)式可以化成其中第一项即(4-12)式给出的近似刚度矩阵，第二项是它的修正部分(4-30)第39页，此课件共60页哦 (4-31)第40页，此课件共60页哦上式使用了缩写符号(4-32)按照(4-3)式计算修正项

13、必须先求出In，也就是需要求出上述三个积分，它的积分区域在三角形ijm上（左图）。当具体积分时，可以在三个梯形i1imm1、m1mjj1和i1ijj1上进行，由前面二个梯形面积分的和中减去第三个梯形面积分。经过运算并加以整理后，得到第41页，此课件共60页哦其中在(4-33)式中的和号 表示大括号中各项对i,j,m轮换后再相加。(4-33)(4-34)第42页，此课件共60页哦 利用(4-33)式进行计算时，有两种特殊情况会出现奇异性。一种情况是当单元的某个结点i位于对称轴上（即ri=0），于是(4-33)式中包含对数的项变成0。这项的极限可由罗彼塔(LHopital)法则来确定，事实上这个极

14、限总是零。因此，如果单元的某个结点的r为零时，只要将其对应的对数项除去即可。第二种情况是当两个结点的径向坐标相等时，例如rj=rm即ci=0；此时，相应的Amj、Bmj都成为无限大，利用(4-33)式计算又会引起困难。然而，在此情况下，梯形j1jmm1的面积等于零，因此，对于jm无需进行积分（上页右图），我们只需令系数Amj、Bmj等于零。第43页，此课件共60页哦 结合以上两种情况，可以将(4-34)式改写成如下的形式以消除计算中的奇异性，即(4-35)第44页，此课件共60页哦 对于有两个结点在对称轴上的单元，例如结点i和j，则有ri=rj=0，此时除了子矩阵 外，其它子矩阵的积分均发散。

15、但是，由于此时ui=uj=0，正好可将刚度矩阵中的第2i-1行和2i-1列以及2j-1行和2j-1列划去。因此，在计算程序中，可以统一应用公式(4-35)使(4-33)式的右端获得一个有限值。这些有限值出现在将要划去的行和列上，因此，并不影响计算的正确性。第45页，此课件共60页哦 实际上，一些例子表明，精确刚度矩阵与近似刚度矩阵，对于主对角元，修正影响甚微，修正项在精确刚度矩阵的相应元素中所占的比例均在5%以下；对于非主对角元，一些元素的修正项所占的比例可达33%。尽管如此，对于位移与应力的计算，其影响并不大。如果采取较大的单元进行同样的计算，结论基本是相同的。因此，近似刚度矩阵（单元不很大

16、的情况下）就位移和应力而言，在工程误差(5%)要求的范围内完全可以采用的。但使用近似等效结点力通常对计算结果影响较大。所以使用精确等效结点力，在实际计算中是有意义的。第46页，此课件共60页哦 实际上，原来三角形环单元现在将用八结点等参元代替，它的位移公式与平面八结点等参元相同，即坐标变换式和位移模式分别具有下列形式(4-36)4-6 4-6 轴对称等参单元轴对称等参单元第47页，此课件共60页哦应变的计算公式(4-37)第48页，此课件共60页哦其中记号 和 分别表示Ni对r和z的偏导数(4-39)(4-38)第49页，此课件共60页哦其中 和 就是平面八结点等参元中的(3-29)式，矩阵是

17、雅可比矩阵 的逆矩阵雅可比行列式是(4-41)(4-40)第50页，此课件共60页哦其中 (4-42)应力公式形式上与前面轴对称三角形单元相应公式相同，即 (4-43)第51页，此课件共60页哦但是 (4-44)其中常数(4-45)第52页，此课件共60页哦 当r=0时，即对称轴上有 ，因此，可以用 来代替 单元刚度矩阵具有以下的形式其中每一个子矩阵是(4-47)以消除(4-44)式中的奇异项。(4-46)第53页，此课件共60页哦而(4-48)式中常数A1、A2、A3按(4-45)式计算。第54页，此课件共60页哦 等效结点力可以按照以下公式计算 (1)集中力 设在单元径向坐标为rc的某点作

18、用有集中力移置到结点上的等效结点力，可按下式计算式中(Ni)c是形函数Ni在载荷作用点的值。(4-49)第55页，此课件共60页哦 (2)体积力 设单元体积力是 ，则作用在单元各结点上的等效结点力分别按下式计算式中r用(4-36)式的第一式代入即可。(4-50)第56页，此课件共60页哦 (3)表面力 设单元的某边界上作用有表面力 ，其中 和分别是单元表面力在作用边外法线方向和切线方向的投影。于是，此边上各结点的等效结点力分别按下式计算：式中r可用(4-36)式的第一式代入。比如，是对应 的那条边，那么(4-51)第57页，此课件共60页哦 (4)温度应力 如果考虑温度变化产生的初应变则任意结点上的等效结点力是(4-52)第58页，此课件共60页哦 计算应力的公式(5-52)(浙大有限元P65)式要改写为 (4-53)第59页，此课件共60页哦10.什么叫轴对称问题？其位移有何特点？11.试定义具有线性位移模式的三角形环单元。12.三角形环单元的应力和应变的特点。第四次作第四次作业业：第60页，此课件共60页哦