有限元 平面问题的高阶单元.pptVIP

3.1 四结点矩形单元 (PLANE42) (1) 位移模式选择：两个，每个四项， (2) 形函数计算 在式（3-1）中代入结点位移和结点坐标后，可解出待定系数。将这些系数再代入式（3-1），可得形函数： (3) 单元应变 (4) 单元应力： (5) 单元刚度矩阵 把[B]、[D]代入上式，整理后可得： (6)等效结点荷载 单元的体积力和表面力引起的结点力仍可用式（1-36）和（1-37）进行计算。由于位移分量在x为常数及y为常数的直线上是线性变化的，因此，载荷向结点的分配也符合静力等效的原则。 (7) 整体平衡方程 根据各单元的刚度矩阵[k]、等效结点力列阵，按对号入座的方式叠加组装整体刚度矩阵和结点荷载列阵，从而得到整体平衡方程： 3.2 六结点三角形单元(PLANE2) 1、位移模式 在三角形单元i、j、m的各边中点上，各增设一个结点，使每个单元具有6个结点，则得到图示的六结点三角形单元。这种单元具有12个自由度，可以采用完全二次多项式的位移模式： 收敛条件：所取位移模式反映了单元的刚体位移和常应变。另外，在单元边界上位移分量是按抛物线变化的，而每条公共边界上有3个公共结点，正好可以保证相邻两单元位移的连续性。因此，上述位移模式满足收敛的完备协调性条件。 根据假定的位移模式，按照前面的方法和过程，可以确定形函数[N]、应变矩阵[B]、应力