第3讲 弹性问题有限元方法.ppt

北京航空航天大学 单元形状函数的性质（m个节点） 在单元节点上形状函数的值为1或为0。 在单元中的任意一点上，形状函数之和等于1 。 当i=j时 当i!=j时 (i, j=1, m) 北京航空航天大学 单元刚度矩阵的特点 对称： 2D: 3D: 弹性矩阵D的对称性决定 北京航空航天大学 奇异：经行变换出现零行或列变换出现零列 主元恒正：位移和作用力同向 北京航空航天大学 整体刚度矩阵的特点 对称 由单元刚度矩阵的对称性所决定 奇异 由单元刚度矩阵的奇异性所决定 稀疏 整体刚度矩阵的多数元素为零，非零元素的个数只占较小的部分。 非零元素带状分布 整体刚度矩阵的非零元素分布在以对角线为中心的带形区域内 北京航空航天大学 非零元素带状分布 b: 半带宽 b 北京航空航天大学 北京航空航天大学 节点编号时，尽量使单元节点最大差值最小。 北京航空航天大学 (a)编号半带宽：b=10 (b)编号半带宽：b=14 (c) 编号半带宽：b=18 (a) (b) (c) b 判断：哪种编号方案最好？ 北京航空航天大学 位移边界条件的处理 对角元素乘大数法 北京航空航天大学 位移单元解的下限性 系统表现过刚 F 证明过程参考有关文献 北京航空航天大学 位移解的下限性质可以解释如下： 单元原是连续体的一部分，具有无限多个自由度。在假定了单元的位移函数后，自由度限制为只有以结点位移表示的有限自由度，即位移函数对单元的变形进行了约束和限制，使单元的刚度较实际连续体加强了，因此连续体的整体刚度随之增加，离散后的K较实际的K为大，因此求得的位移近似解总体上（而不是每一点）将小于精确解。 北京航空航天大学 单元等效节点力列阵 单元外力功： 北京航空航天大学 扩充叠加 扩充叠加 Step 4. 单元集成——系统势能 北京航空航天大学 关于单元刚度矩阵的扩充叠加 m i j m i j 单元编号 mij 北京航空航天大学 关于单元等效节点载荷列阵的扩充叠加 m i j 单元编号 mij 北京航空航天大学 Step 5. 变分处理 Step 6: 处理位移边界条件并求解 Step 7: 计算每个单元的应变及应力 北京航空航天大学 关于三角形单元 3节点三角形单元是常应变（常应力）单元。在应变梯度较大的部位（亦即应力梯度较大的部位），单元划分应适当密集，否则不能反映真实的应变变化而导致较大的误差。 提高计算精度的其它措施 采用高精度三角形单元（2次单元、3次单元…） 采用四边形单元（1次单元、2次单元…） 北京航空航天大学 4节点四边形单元 北京航空航天大学 构造位移函数： 对u,v分别利用节点条件： 北京航空航天大学 对于一般四边形，逆矩阵的表达式比较复杂。 北京航空航天大学 北京航空航天大学 N—单元形状函数矩阵qe —单元节点位移矩阵 北京航空航天大学 特例：4节点矩形单元 矩形单元的重心坐标 北京航空航天大学 北京航空航天大学 对于一般的四边形单元，在总体坐标系下构造位移插值函数，则计算形状函数矩阵、单元刚度矩阵及等效节点载荷列阵时十分冗繁；而对于矩形单元，相应的计算要简单的多。 矩形单元明显的缺点是不能很好的符合曲线边界，因此可以采用矩形单元和三角形单元混合使用。更为一般的方法是通过等参变换将局部自然坐标系内的规格化矩形单元变换为总体坐标系内的任意四边形单元（包括高次曲边四边形单元）。 北京航空航天大学 受均布内压作用的长圆筒 轴对称问题的有限元求解过程 北京航空航天大学 研究轴对称问题时通常采用圆柱坐标系(r,θ,z)，以z轴为对称轴 由于对称性: 4个应力分量，4个应变分量，2个位移分量 北京航空航天大学 3节点三角形轴对称单元 360度环形单元 横截面为3节点三角形 r z平面内，单元节点位移和节点力为： 北京航空航天大学 构造位移函数： 利用节点条件： 北京航空航天大学 北京航空航天大学 北京航空航天大学 三维问题的有限元求解过程 离散时采用体单元：四面体或六面体 求解步骤和平面问题完全一样 单元分析的时候将二维扩充到三维 北京航空航天大学 形状函数矩阵 北京航空航天大学 应变矩阵 北京航空航天大学 对称 单元刚度矩阵 北京航空航天大学 单元等效节点力列阵 北京航空航天大学 3.3 有限元分析中的若干问题探讨 单元位移函数的构造 单元形状函数的性质 单元刚度矩阵的特点 整体刚度矩阵的特点 位移边界条件的处理 位移单元解的下限性 北京航空航天大学 单元位移函数的构造满足收敛性要求 收敛——单元尺寸趋于零时，有限元解趋于真解 北京航空航天大学 准则1：完备性—包含常应变项和刚体位移项 如果在势能泛函中所出现的位移函数的最高阶导数是m阶，则选取的位移函数至少是m阶完全多项式。 准则2：协调性—相邻单元公共边界保持位移连续 如果在势能泛函中