0403空间等参数单元的力学分析.docVIP

§4-3 空间等参数单元的力学分析 现在来对n个节点的空间等参数单元进行力学分析，从而建立这种单元的载荷列阵，应力矩阵和刚度矩阵。 一、等效节点载荷计算 空间等参数单元的载荷移置，仍遵循静力等效原则，因而前面几章介绍的普遍公式仍然适用。 1. 集中力 当单元在任意点受有集中力时，等效节点载何为 (4-32) 但在这里 (4-33) 而矩阵［N］可用三阶单位矩阵I及各个形函数表示成为 (4-34) 由于其中Ni是由受力点的局部坐标ξ、η、ζ表示的，为计算方便，单元划分时应将集中力作用点取为节点。否则 ，由整体坐标求局部坐标很不方便。 2. 体力 单元受有分布体力，其强度为，此时有 上式中［N］是局部坐标的显函数，且单元的真实形状及边界描述起来非常困难，一般将上式化成局部坐标积分比较方便。利用(4-24)式，有 (4-35) 对于变量的分布体力，须将体力分量X、Y、Z表示成局部坐标ξ，η，ζ的函数，再进行积分。 3. 面力 当单元在某一边界上(例如在ξ=±1的ηζ面上)，受有分布的面力时，可利用(4-34)式的积分求得等效节点载荷 应用(4-29)式，可以将上式表示成为 (4-36) 对于变量的分布面力，须将面力分量、、表示成局部坐标的函数，再进行积分。 如果上述分布面力为静水压力q，则 将(4-31)式代入，然后再代入(4-36)式，得 (4-37) 当面力作用在η＝±１的ζξ面上或ζ＝±1的ξη面上时，只须将(4-37)中的ξ，η，ζ轮换，即可得出各种状况下的载荷列阵表达式。 二、应变、应力矩阵 1．应变 将位移函数表达式(4-17)代入空间问题的几何方程，可得单元内的应变表达式 (4-38) 其中是单元的节点位移列阵，即 (4-39) 而 (i=1,2,…，n) (4-40) 2．应力 和以前一样，单元内的应力可以表示成 (4-41) 将其中应力矩阵［S］写成分块形式 (4-42) 则有 (i=1,2,…，n) (4-43) 三、刚度矩阵 将单元上的节点力表示成为 则由虚功方程可以得出和以前一样的公式 其中单元刚度矩阵为 (4-44) 单元刚度矩阵是一个3n×3n的矩阵，可以写成分块形式： (4-45) 其中每一个子矩阵是3×3方阵（3个自由度时），为 (i=1,2,…,n;j=1,2,…,n) (4-46) 有了单元的载荷列阵和刚度矩阵之后，即可利用以前多次用过的方法，建立那些求解节点位移所需的节点平衡方程，在求出节点位移之后，再按(4-41)式求单元的应力。 四、几点说明 1．坐标变换 对于等参单元中的任一点，由它的局部坐标求出它的整体坐标时，只须直接应用(4-18)式； 要由整体坐标算出局部坐标，则须利用(4-18)式，求解三元非线性联立方程组，这是非常繁琐的，甚至不可能。 2．参数计算方法 在计算单元中的应力等参数时，一般是先设定一组局部坐标（一般是节点或积分点），代入应力矩阵以求出应力；并根据这组局部坐标算出相应的整体坐标，从而可知所求应力是单元(构件)那一点处的应力。对于单元中以整体坐标标定的一点（除非是节点或积分点等特殊点），要计算应力则是非常复杂和不现实的。 计算节点处的应力，则比较方便。因为其局部坐标是已知的，且非常简单，不外乎0，±1以及简单的分数±，±等等。 因此，在整理应力等计算结果时，一般宜采用绕点平均法。 3．精度问题 对于等参数单元，单元内点的应力精度优于边界点的精度，角节点的应力精度最差。 如果单元的形态比较差，则计算求得的节点处的应力的表征性也就较差，通过平均后其表征性可能仍然不够好。对于这样的单元，就应整理单元内点处的应力，由内点处的应力来外插求节点的应力。 3