四节点矩形单元有限元分析演示文稿

四节点矩形单元有限元分析演示文稿当前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结论：

对于平面四节点矩形单元，其单元上的应力、应变不再是常数，而是在一定程度上呈线性变化，即：方向的正应力和正应变随坐标线性变化；方向的正应力和正应变随坐标线性变化；剪应力沿坐标和坐标均成线性变化。因此，若在弹性体中采用相同数目的节点时，矩形单元的精度要比常应变三角形单元的精度高。

当前18页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity七、运用虚功原理求解由虚功原理，节点力在节点的虚位移上所做的虚功应等于单元内部应力在虚应变上所做的虚功，即内力虚功=外力虚功，也即：

将和代入上式，可得：由此，可得：当前19页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity当前20页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity引入无量纲坐标：当前21页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity由几何方程可得单元应变场表达式：可记为：当前22页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity几何矩阵可表示成分块形式：其中：当前23页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity由应力与应变关系，可得单元应力场表达式：应力矩阵可表示成分块形式：其中：对于平面应变问题：当前24页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity其中：即：单元刚度矩阵：当前25页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity其中：当前26页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity八、由合成刚度集成法：首先求出各单元的贡献矩阵，然后将它们叠加形成整体刚度矩阵。但是由于编程时需先将各单元的贡献矩阵储存起来，而贡献矩阵的阶数与整体刚度矩阵阶数相同，因此占用非常大空间，不利于节约空间资源。单元定位数组法：将单元的节点位移编码按照节点顺序排成一行形成一维数组，利用各单元的定位数组，采用“边定位，边累加”的方法。

当前27页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity九、建立节点荷载列阵节点荷载列阵的组成：

其中，为节点荷载，为等效节点荷载。

可按照虚功等效原则求解，即将单元内的荷载移置到节点上后，应当与原荷载所作虚功等效。集中力、分布体力（均质等厚单元自重）、分布面力（均布侧压、X方向均布荷载、X方向三角形荷载）

当前28页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity十、处理位移约束条件（1）降阶法：若第r个自由度方向位移分量为0，则将整体刚度矩阵第r行，第r列划掉，后一行上移，右一列左移，这样总刚减少一阶，未知数减少一个。

例：

当前29页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity（2）对角元素置1法：

例：已知位移边界条件（可为零）

当前30页，总共32页。四节点矩形单元有限元分析过程

TianjinUniversity（3）对角元素乘大数法：

例：已知位移边界条件（可为零）

略去小量有即当前31页，总共32页。总结

TianjinUniversity四节点矩形单元采用双线性位移插值函数，应力和应变沿坐标轴呈线性变化，因而求解精度比三节点三角形单元高。由于位移