四节点矩形单元有限元分析.pptVIP

四节点矩形单元有限元分析演示文稿现在是1页\一共有32页\编辑于星期六

本节内容提要1、分析提高有限元法求解精度的途径2、简要回顾三节点三角形单元有限元分析过程3、全面介绍四节点矩形单元有限元分析过程4、总结TianjinUniversity现在是2页\一共有32页\编辑于星期六

分析提高有限元求解精度的途径TianjinUniversity一、三节点三角形单元的缺点三节点三角形单元精度低，收敛慢，由于单元内应力和应变均为常量，故在单元内不能很好地反映应力和应变的变化。该单元只有三个节点，单元自由度少，单元位移插值函数（位移模式）只能是线性函数，描述单元内位移变化的能力差。现在是3页\一共有32页\编辑于星期六

分析提高有限元求解精度的途径TianjinUniversity二、提高有限元求解精度的途径第一个途径是对某一种特定类型的单元采用网格划分加密，依靠单元的收敛性提高求解精度。第二个途径是对一定的单元网格和单元尺寸，采用高精度单元来提高求解精度。现在是4页\一共有32页\编辑于星期六

分析提高有限元求解精度的途径TianjinUniversity三、建立高精度单元的原理和途径原理：提高单元位移插值函数多项式的阶次，从而提单元拟合局部区域位移、应力变化的能力。途径：增加单元的节点数目。对于平面有限元问题，除三节点三角形单元外，还可以考虑六节点三角形单元和四节点矩形单元。现在是5页\一共有32页\编辑于星期六

三节点三角形单元有限元分析过程设位移函数求位移函数中的未知量代入函数中整理可得形函数（性质？）几何方程求解应变（几何矩阵）物理方程求解应力（弹性矩阵应力矩阵）运用虚功原理求解由合成(方法？)建立节点荷载列阵（方法？处理位移约束条件（方法？）组成？)TianjinUniversity现在是6页\一共有32页\编辑于星期六

TianjinUniversity四节点矩形单元有限元分析过程一、四节点矩形单元位移函数单元节点编号为k，l，m，n（逆时针）单元节点位移列阵为：设位移函数为：或写为：现在是7页\一共有32页\编辑于星期六

TianjinUniversity四节点矩形单元有限元分析过程二、求解位移函数中的未知系数将节点坐标代入函数中，并写成矩阵形式：解上述方程组可得：现在是8页\一共有32页\编辑于星期六

TianjinUniversity四节点矩形单元有限元分析过程三、将所求值代入位移函数中现在是9页\一共有32页\编辑于星期六

TianjinUniversity四节点矩形单元有限元分析过程四、整理位移函数可得形函数展开上式可得：现在是10页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程其中，形函数为：TianjinUniversity现在是11页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程单元位移插值函数可以由单元形状函数与节点位移值的乘积表示：即可以表示为：TianjinUniversity现在是12页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程由此可见，位移插值函数完全由形函数决定；因此抛开节点位移，只讨论形函数的性质，就可以了解单元的变形性质。例如：四节点矩形单元，若则由可得：因此可以看出，单元变形完全由形函数决定。TianjinUniversity现在是13页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程另外，可以验证形函数另外两个性质：（1）同理对于其余三个形函数TianjinUniversity现在是14页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程TianjinUniversity（2）即在单元内任意一点处的形函数之和等于1。现在是15页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程TianjinUniversity五、几何方程求解应变将位移插值函数代入几何方程中：形函数矩阵经过微分算子矩阵作用后得到3×8几何矩阵：现在是16页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程TianjinUniversity六、物理方程求解应力由平面问题物理方程可得：其中：因此，应力矩阵为：现在是17页\一共有32页\编辑于星期六

四节点矩形单元有限元分析过程TianjinUniversity结论：对于平面四节点矩形单元，其单元上的应力、应变不再是常数，而是在一定程度上呈线性变化，即：方向的正应力和正应变随坐标线性变化；方向的正应力和正应变随坐标