02-04 单元刚度矩阵

《02-04 单元刚度矩阵》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、§2-4单元刚度矩阵第四步：利用平衡方程，建立节点力和节点位移之间的关系，即用单元节点位移表示节点力。上节己给出了用节点位移表示单元应力和应变。本节来推导单元节点力和节点位移之间的关系。一、节点力和节点位移间的关系节点力是指弹性体离散化之后，外载、约束和其他单元通过节点作用在某一单元上的力。图2-11三角形单元的节点位移和节点力对于己从整体结构中取出来的单元来说，作用在其上的节点力就是外力。这些节点力在单元内部会引起相应的应力。当整体处于平衡状态时，单元在节点力作用下也处于平衡状态。在平面问题中节点力有二个分量，分别用U和V加节点号下标表示该节点水平和垂直节点力分量(有时还再加

2、单元号上标表示该单元上的节点力)。节点力的方向以节点对单元的力沿坐标正方向为正，反之为负。对三节点三角形单元来讲，共有六个节点力分量(如图2-11所示)。用列阵表示为：；(2-24)1．虚位移原理为了推导单元的节点力与单元节点位移之间的关系，要用到虚位移原理。2．节点力和节点位移间的关系虚位移原理在一处于平衡状态的单元上的数学描述为：单元上节点力(外力)在某一虚位移上所作的虚功应等于单元应力(内力)在相应虚应变上所作的虚功。设单元节点处的虚位移为；=(i,j,m)(2-25)采用和真实位移相同的位移模式，则单元内各点的虚位移为(a)相应虚应变为(b)7于是虚功方程可写成(2

3、-26)将(b)式及(2-18)式代入上式，得根据矩阵乘法逆序法则，上式可以写成由于列阵中的元素是常量，即与单元内点的位置坐标x，y无关，上式右边的可以提到积分号前面去。又由于虚位移可以是任意的，所以矩阵也是任意的，因此等式两边与它相乘的矩阵应当相等，于是得(2-27)令(2-28)则(2-27)式得(2-29)这就是节点力与节点位移之间的关系式。由于〔D〕中的元素是常量，而且在线性位移模式的情况下，〔B〕中的元素也是常量，同时，(单元面积)，则公式(2-27)与(2-28)就简化为(2-30)(2-31)矩阵反映了单元节点力与节点位移之间的关系，所以称其为

4、单元的刚度矩阵。矩阵表示该单元各节点沿坐标方向发生单位位移时引起的节点力，它决定于单元的形状，大小，方位和弹性常数，而与单元的位置无关，即不随单元或坐标轴的平行移动而改变。7一、单元刚度矩阵的计算为了计算单元刚度矩阵的值，将(2-31)式展开可得(2-32)对于6自由度的三角形三节点单元，其单元刚度矩阵为6×6阶方阵。其中子矩阵(r,s=i,j,m)是九个2×2阶方阵，它们写成一个统一的计算格式，即(r,s=i,j,m)(2-33)式中(r，s=i，j，m)对于平面应力问题(r,s=i,j,m)(2-34)对于平面应变问题(r,s=i,j,m)(2-35)7一、单元刚度矩阵的物理

5、意义单元刚度矩阵表示三角形单元抵抗变形的能力，即表示了节点位移与节点力之间的关系。将(2-29)式写成下面形式(2-36)根据矩阵乘法可得(2-37)从(2-37)式的第一式可以看出，、、分别表示使节点i、节点j与节点m产生单位位移时，在节点i上所需要施加的节点力的大小，即在单元节点i处产生的作用力。u则表示节点j产生位移时，在节点i处所产生的弹性力，其余类推。由此可知：u子矩阵表示使节点s(s=i,j,m〕产生单位位移时，在节点r(r=i,j,m)处所产生弹性力的大小。现将单元刚度矩阵子矩阵展开成如下形式并将(2-37)式中节点力列阵与节点位移列阵均展开成6×1阶列矩阵，则（2

6、-29）式所列平衡方程为(2-38)7由此可见，单元刚度矩阵中的每一个元素都有其明显的物理意义，即表示后一下标所示的节点的某一方向上产生单位位移时，在前一下标所示的节点的某一方向上所产生的作用力。其中元素(r、s=i,j,m)的下标上有“一”的表示水平方向，没有“一”的表示垂直方向。例如：u表示节点i在水平方向产生单位位移时，在节点i所产生的水平作用力。u表示节点j在垂直方向产生单位位移时，在节点i所产生的水平作用力。u表示节点j在水平方向产生单位位移时，在节点i所产生的垂直作用力。u表示节点j在垂直方向产生单位位移时，在节点i所产生的垂直作用力。其余类推。同时可以得出：单元一个

7、节点上的节点力分量是本单元上三个节点六个位移分量对其提供贡献的总和，即r=i,j,m(2-39)这其中包括该节点对自己做出的贡献和其它节点对该节点提供的贡献，即节点之间互相做出贡献。但应注意分清所求节点力之点(r)与提供贡献之点(s)，即s节点对r节点提供贡献，其单元刚度子矩阵为。一、单元刚度矩阵的特性1.对称性，就是说单元刚度矩阵是一个对称矩阵。根据(2-34)式，将单元刚度矩阵中各元素相互对照便得余类推。由此可见，在单元刚度矩阵这个6×6阶方阵中，对称于主对角线的元素都两两相