第4节 杆系结构单元.ppt

4.6 坐标变换矩阵计算式 坐标变换矩阵的计算式因单元类型不同而异。 1、平面铰接杆单元 设OXY为结构坐标，oxy为单元坐标。?为从单元 i 端出发的任一矢量。它在结构坐标系中的分量为?X、?Y；在单元坐标系中的分量为?x、?y。结构坐标系中的分量?X、?Y 在单元坐标x轴上投影的代数和给出?x 。同理， ?X、?Y 在单元坐标y轴上投影的代数和给出?y 。 由图4-12得： X Y x y ?X ?Y ?x ?y 图4-12 （4-61） ? 写成矩阵形式， i ? a b c ab+bc 取 ——i节点在单元坐标系中的位移向量 ——i节点在结构坐标系中的位移向量 x对X、Y的方向余弦 y对X、Y的方向余弦 同理可得单元j节点在单元坐标系和结构坐标系中的位移向量： 有 组合上述结果，得平面铰接杆单元的单元坐标单元位移和结构坐标单元位移之间关系： i、j两节点间的位移变换关系互不耦合。 上式可写成 坐标变换矩阵[R]的计算式： 用单元节点在结构坐标系中的坐标值（Xi、Yi和Xj、Yj ）描述方向余弦： （4-62） ○ ○ ○ ○ ○ X Y i j (e) x （4-62a） ? 式中，（Xi，Yi）和（Xj，Yj）分别为节点i和节点j在结构坐标系中的坐标值。 2、两端承受剪力、弯矩的平面梁单元 如果在连续梁中使用这类单元，通常可将单元坐标和结构坐标取得一致。此时，无须进行坐标变换。 ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ ○ i j (e) X Y i j (e) x y 3、两端承受轴力、剪力、弯矩的平面梁单元 此时，节点自由度为3，见图4-13。 i j x y ?2 ?3 ?5 ?6 l ?1 ?4 图4-13 i j x y ?2 ?3 ?5 ?6 l ?1 ?4 于是得到： X Y 注意?1 、 ?2 、?4、 ?5的性质和平面铰接杆相同，因而有相同的[R]矩阵。又，单元坐标系xy平面和结构坐标系XY平面在同一平面上， 因而单元坐标系z轴和结构坐标系的Z轴总有相同指向，所以恒有： （4-63） i j x y z 4、两端承受扭矩和面外剪力、弯矩的平面梁单元 此时，xy平面和结构坐标系XY平面仍在同一平面上，因而z轴和结构坐标系的Z轴指向相同，只须取 l 图4-14 ?xi ?xj ?yi ?yj Wi Wj X Y Z 在单元坐标系中，单元每个节点（如i）有3个位移分量： ?xi、?yi和wi (5) 等价节点力 对于梁上作用的集中力或集中力矩，在划分单元时可将其作用点取为结点，按结构的节点载荷处理。 这里考虑的是把单元上的横向分布载荷转化为等价节点力问题。当梁单元上作用有横向分布荷载qy(x)时（图4-6）， x y i j l 图4-6 qy(x) 图4-7 l x y i j qy(x) i j x dx v(x) qy (x)dx 横向分布荷载qy(x)的势能Vq为： （4-36） 形函数矩阵由式（4-30）和（4-31）给出。对于具体问题，只要将qy(x)代入上式进行积分即可。表1给出了几种特殊情况的等价节点力。 - 5ql2/96 ql/4 5ql2/96 ql/4 - ql2/20 7ql/20 ql2/30 3ql/20 - ql2/12 ql/2 ql2/12 ql/2 Mj Qj Mi Qi 荷载分布 i j q q i j q i j 几种横向分布荷载等价节点力 表 1 (6) 单元坐标单元刚度矩阵 梁单元刚度矩阵公式为 将式（4-34）代入上式进行积分，并注意到 Iz——梁截面对Z轴（主轴）的惯性矩 得单元坐标单元刚度矩阵[k]e: （4-37） 单元刚度矩阵式(4-38)适合于连续梁分析。 （4-38） 2、两端承受轴力、剪力、弯矩的平面梁单元 i j x y i j x y ?2 ?3 ?5 ?6 l ?1 ?4 F2 F3 F5 F6 l F1 F4 图4-8 （1）单元坐标单元位移和单元力 ① 单元位移 （4-39） 其中， u——x方向（轴向）位移。 v——y方向位移，即挠度。 ?——角位移。 ② 单元力 （4-40） 其中， N——轴向力