平面问题有限元分析四节点矩形单元.ppt

四节点矩形单元有限元分析过程 六、物理方程求解应力 由平面问题物理方程可得： 其中： 因此，应力矩阵 为： 四节点矩形单元有限元分析过程 结论： 对于平面四节点矩形单元，其单元上的应力、应变不再是常数，而是在一定程度上呈线性变化，即： 方向的正应力和正应变随 坐标线性变化； 方向的正应力和正应变随 坐标线性变化；剪应力沿 坐标和 坐标均成线性变化。 因此，若在弹性体中采用相同数目的节点时，矩形单元的精度要比常应变三角形单元的精度高。 四节点矩形单元有限元分析过程 七、运用虚功原理求解 由虚功原理，节点力在节点的虚位移上所做的虚功应等于单元内部应力在虚应变上所做的虚功，即内力虚功=外力虚功，也即： 将 和 代入上式，可得： 由此，可得： 四节点矩形单元有限元分析过程 四节点矩形单元有限元分析过程 引入无量纲坐标： 四节点矩形单元有限元分析过程 由几何方程可得单元应变场表达式： 可记为： 四节点矩形单元有限元分析过程 几何矩阵可表示成分块形式： 其中： 四节点矩形单元有限元分析过程 由应力与应变关系，可得单元应力场表达式： 应力矩阵可表示成分块形式： 其中： 对于平面应变问题： 四节点矩形单元有限元分析过程 其中： 即： 单元刚度矩阵： 四节点矩形单元有限元分析过程 其中： 四节点矩形单元有限元分析过程 八、由 合成 刚度集成法：首先求出各单元的贡献矩阵，然后将它们叠加形成整体刚度矩阵。但是由于编程时需先将各单元的贡献矩阵储存起来，而贡献矩阵的阶数与整体刚度矩阵阶数相同，因此占用非常大空间，不利于节约空间资源。 单元定位数组法：将单元的节点位移编码按照节点顺序排成一行形成一维数组，利用各单元的定位数组，采用“边定位，边累加”的方法。 四节点矩形单元有限元分析过程 九、建立节点荷载列阵 节点荷载列阵的组成： 其中， 为节点荷载， 为等效节点荷载。 可按照虚功等效原则求解，即将单元内的荷载移置到节点上后，应当与原荷载所作虚功等效。 集中力、分布体力（均质等厚单元自重）、分布面力（均布侧压、X方向均布荷载、X方向三角形荷载） 四节点矩形单元有限元分析过程 十、处理位移约束条件 （1）降阶法：若第r个自由度方向位移分量为0，则将整体刚度矩阵第r行，第r列划掉，后一行上移，右一列左移，这样总刚减少一阶，未知数减少一个。 例： 四节点矩形单元有限元分析过程 （2）对角元素置1法： 例：已知位移边界条件 （可为零）  平面问题有限元分析 四节点矩形单元 本节内容提要 1、分析提高有限元法求解精度的途径 2、简要回顾三节点三角形单元有限元分析过程 3、全面介绍四节点矩形单元有限元分析过程 4、总结 分析提高有限元求解精度的途径 一、三节点三角形单元的缺点 三节点三角形单元精度低，收敛慢，由于单元内应力和应变均为常量，故在单元内不能很好地反映应力和应变的变化。 该单元只有三个节点，单元自由度少，单元位移插值函数（位移模式）只能是线性函数，描述单元内位移变化的能力差。 分析提高有限元求解精度的途径 二、提高有限元求解精度的途径 第一个途径是对某一种特定类型的单元采用网格划分加密，依靠单元的收敛性提高求解精度。 第二个途径是对一定的单元网格和单元尺寸，采用高精度单元来提高求解精度。 分析提高有限元求解精度的途径 三、建立高精度单元的原理和途径 原理：提高单元位移插值函数多项式的阶次，从而提单元拟合局部区域位移、应力变化的能力。 途径：增加单元的节点数目。 对于平面有限元问题，除三节点三角形单元外，还可以考虑六节点三角形单元和四节点矩形单元。 三节点三角形单元有限元分析过程 设位移函数 求位移函数中的未知量 代入函数中 整理可得形函数 （性质？） 几何方程求解应变（几何矩阵 ） 物理方程求解应力 （弹性矩阵 应力矩阵 ） 运用虚功原理求解 由 合成