ch18 杆与杆系分析的计算机方法(3rd).docVIP

第18章 杆与杆系分析的计算机方法 18-1 试建立桁架单元在局部坐标系中的刚度矩阵。 解：在线弹性范围内，设在局部坐标系中杆单元的刚度方程为 f = k d (a) 其中，d,f和k分别为杆单元的端点位移矢量（图18-1a）、端点力矢量（图18-1b）及杆单元刚度矩阵，具体为 d = (ui vi uj vj)T (b) f = (Fxi Fyi Fxj Fyj)T (c) 由于桁架杆单元为二力杆，故有Fyi = Fyj = 0,式（c）便成为 f = (Fxi 0 Fxj 0)T (c′) 图18-1 下面推导杆单元刚度矩阵k。 假设单元的j端固定，使i端产生位移ui（图18-1c），从而有 由平衡条件可得 假设单元的i端固定，使j端产生位移uj（图18-1d），从而有 由平衡条件可得 将以上两组端点力叠加，得 (d) 将式（d）写成矩阵形式并扩大为包括所有节点位移的形式 f = (e) 对比式（a）与式（e），得到桁架单元在局部坐标系中的刚度矩阵为 k 18-2 试建立桁架单元在整体坐标系中的刚度矩阵。 解：在整体坐标系中，设杆单元的方位角为(（见图18-2a），那么在局部坐标系中的端点位移与整体坐标系中的端点位移的关系为 (a) (b) 图18-2 整体坐标系中的端点位移矢量为 (c) 将局部坐标系中的端点位移矢量 d (d) 写成列阵形式，并注意到式（a）和（b），可得 d 令 T= (e) 则两种端点位移矢量之间的关系可简记为 d = T (f) 其中，T称为转换矩阵。 由图14-2（b）与图14-2（a）的相似性不难找到，两种端点力矢量之间的关系为 (g) 其中， 式（g）中的T由式（e）给出。 将式（g）和(f)代入式f =kd中，有 (h) 两边同乘以逆矩阵T-1后，得 (i) 注意到本题的逆矩阵与转置矩阵是相等的，即 (j) 式（i）还可改写成 (k) 这表明，在整体坐标系中，桁架的单元刚度矩阵为 (l) 其中，T仍由式（e）给出。 18-3 图示连续梁，同时承受集中载荷FP、分布载荷q与矩为Me的集中力偶作用。已知FP = 6kN，q = 4 N/mm，Me= 4×106N·mm，l = 1m，试写出节点载荷矢量。 题18-3图 解：本题求解可以分三步来完成。 首先，求出图18-3（a）所示简支梁两端的转角(1和(3。 图18-3 为此，可先求出图18-3（b）所示简支梁两端的转角，它们分别为 (a) (b) 对于分布载荷情况（见图18-3c），两端的转角分别为 (c) (d) 注意到图18-3（a）中两段的分布载荷依次为 q (0x≤l) (e) (l≤x≤2l) 于是，将式（e）分别代入式（c）和(d)，即可求得该梁两端的转角分别为 (() (f) (() (g) 其次，求出图18-3（d）所示梁两固支端的约束力偶矩。 由变形协调条件(1=0及(3=0，得 (1) (2) 解此联立方程，得到 最后，求题图所示连续梁的节点载荷矢量。 根据上列结果、教材（Ⅱ）表18-1及书中规定的符号，本题有 于是，得该连续梁的节点载荷矢量为 M = 5 q(x)=