有限单元法 第4章 空间轴对称问题的有限元分析.pptVIP

第4章空间轴对称问题的有限元分析 4.1　概述 4.2　三角形单元 4.2.1　位移函数的选取 4.2.2　单元刚度矩阵 4.2.3　等效结点荷载 　　从薄板内取出一个平行六面体，它在x和y方向上具有单位宽度，在z方向的高度为t（图6.3）。在x为常量的横截面上，作用有　和　。由于　和　都和z成正比，所以它们在薄板全厚度上的代数和分别等于零。只可能分别合成弯矩和扭矩。用　表示由　所合成的单位宽度上的弯矩，则得： 将式（6-8）中的　的表达式代入上式，并对z积分，得到 与此相似，应力分量　将合成扭矩 将式（6-8）中的　表达式代入上式，并对z积分，得到 （6-11） （6-10） 同样，在y为常量的横截面上，每单位宽度上的　和　也分别合成如下的弯矩和扭矩 　　在这里可以看到，由剪应力　和　的互等关系，得到扭矩　和　的互等关系。 　　将式（6-10）、（6-11）和（6-12）中的第一式合并起来，用矩阵表示，则有 （6-12） 其中 是薄板弯曲问题的弹性矩阵，它等于平面应力问题中的弹性矩阵乘以　　　。式（6-13）表示了薄板的内力与应变两者之间的关系，因而是薄板弯曲问题中的物理方程。 （6-13） （6-14） 　　弯矩和扭矩　、　、　、　的方向及其作用面的位置示于图6.4a中。按右手螺旋法则用双箭头矢量来表示力偶，如图6.4b所示（图中所示各力偶的方向均为正）。 　　从式（6-13）中解出　　　　、　　　　、　　　　，再代入式（6-8），就得到各应力分量与内力之间的关系式。 图6.4 （a） （b） 　　从式（6-13）中解出　　　　、　　　　、　　　　，再代入式（6-8），就得到各应力分量与内力之间的关系式。 （6-15） 6.2　三角形单元 6.2.1　位移模式问题 　为保证挠度 w 为坐标的全三次多项式，从帕斯卡三角形可知必须要有10项，但三角形3个结点只能有9个自由度，若舍去三次项中的任一项，显然都无法保证对坐标的不变性，为此Tocher提出了一种解决方案如下： 但是当三角形的二边分别平行坐标x轴、y轴时，从式（6-16）无法通过结点的位移条件来确定广义坐标参数　，为此必须在离散化时设法避免边界（内部）单元不出现上述现象，分析结果表明此位移模式对应的单元能得到很好的结果。 （6-16） 　　另一解决方案是Adini提出的 此方案由于舍去了二次项xy，故常扭率　　　无法得到保证，从而位移偏小，精度很差。这也可以从挠曲函数的Taylor级数展开来说明，由于不包含完全二次项，故只能有一阶精确度。 　　还有一种Bell提出的解决方案，他对三角形单元除了3个角点作为结点外，取形心点挠度作为位移参数，从而取位移模式为全三次多项式 利用10个位移参数条件确定广义坐标参数，通过分析建立起单元刚度、等效结点荷载矩阵，在整体分析之前采取如下措施消去内点自由度。 　　单元刚度方程 　　　　 （6-19） 式中 从式（6-19）的第二个方程可解出： 　　　 （a） 将式（a）代回式（6-19）第一个方程并进行整理，可得： 　　　　 　　设 单元刚度矩阵： 等效结点荷载： 　 （6-20） 则单元刚度方程为 （6-21） 　　通过静力凝聚最后获得了9自由度的三角形单元。但是Zienkiowicz曾指出这样所得到的单元不能保证收敛性，因此，通常采用另外一种位移模式，即下面介绍的面积坐标下的位移模式。 6.2.2　面积坐标下的位移模式 设面积坐标变量为　　　　。其定义如下：如图6.6所示，三角形单元中任意一点P的位置可以由如下三个比值来确定： 　 A为三角形单元的面积， 为小三角形面积 根据三角形面积公式，可以得到面积坐标与直角坐标的关系为： 　　　 　 式中 　　Zienkiowicz等采用面积坐标解决了直角坐标下所产生的困难，提出了如下的位移模式 　　利用结点的位移参数条件可以确定广义坐标参数α，从而建立形函数（当然也可象前面介绍的那样由形函数性质来确定）如下 （a） 式中脚标轮换规则为 1→2→3→1， 如果注意到： 则式（a）可改写成： * 有限单元法 　　教学提示：空间轴对称结构是工程中一类特殊的结构。对于这类结构，由于其几何形状、约束情况、所受荷载都对称于某一根轴（如图4.1所示），因此，通常是在柱坐标系中进行分析。本章以三角形单元为例介绍了这一类结构的有限元分析过程。 　　教学