有限元刚度矩阵推导.docxVIP

有限元刚度矩阵推导 第 PAGE \* Arabic \* MERGEFORMAT 5 页 改图 改图 一平面桁架及其节点和单元的编号如图所示，各杆的截面积均为A，弹性模量均为E，杆1和杆2的长度均为L。求在图示坐标系下杆系的整体刚度矩阵。 解： (1) 推导任意单一杆件的单元刚度矩阵，如图所示，单一杆件有两个节点1、2，杆件的局部坐标系为o’x’y’，整体坐标系为oxy。 图中θ角是从整体坐标系x轴正向起逆时针转到杆件方向，设u’、v’、Fx’、Fy’对应的是局部坐标系下的位移和作用力，u、v、Fx、Fy对应的是整体坐标系下位移和作用力。 位移和作用力从局部坐标系转换到整体坐标系的表达式为： (1) (2) (3) (4) 由单个弹簧刚度矩阵，可得单一杆件在局部坐标系下力与位移的关系矩阵为： (5) 令，，由式(1)-(4)可得 (6) (7) 综合式(5)-(7)可即可得出该单一杆件在总体坐标系下作用力与位移之间的矩阵关系： (8) 该单一杆件在总体坐标系的刚度矩阵即为： (9) (2) 由题意可知该平面桁架共有三个节点，每个节点有两个自由度，因此力和位移各包含6项，即为 ， 各单元杆件的转角θ，以及λ、μ分别为 杆件 θ sinθ cosθ = 1 \* GB3 ① 0° μ=0 λ=1 = 2 \* GB3 ② 90° μ=1 λ=0 = 3 \* GB3 ③ 135° 改角度，注意：杆的长度可能也改变了 改角度，注意：杆的长度可能也改变了 由(9)式可得3根杆件在整体坐标系下力与位移关系矩阵分别为 杆 = 1 \* GB3 ① ， 杆 = 2 \* GB3 ② ， 杆 = 3 \* GB3 ③ 这三个矩阵要算出来，把3表中的λμ带入公式 这三个矩阵要算出来，把3表中的λμ带入公式8，注意杆长 又知节点1由杆 = 1 \* GB3 ①和杆 = 2 \* GB3 ②共有，节点2由杆 = 1 \* GB3 ①和杆 = 3 \* GB3 ③共有，节点3由杆 = 2 \* GB3 ②和杆 = 3 \* GB3 ③共有，因此，由叠加原理可得： 按照公式算一遍，注意：这里隐藏了AE/L这部分公式 按照公式算一遍，注意：这里隐藏了AE/L这部分公式 将以上六式写成矩阵的形式为 把上面的六个算式写成矩阵形式就可以了 把上面的六个算式写成矩阵形式就可以了 即为该平面桁架在给定坐标系下的整体刚度矩阵。