物理矩阵位移法.pptxVIP

物理矩阵位移法会计学第1页/共84页二、结构离散化 用矩阵位移法求解，首先要将结构离散成单元和结点的组合体系。具体做法就是对结构进行结点编号和单元编号。 离散化的关键就是确定结构的全部结点。因为单元的两端是结点，单元与单元之间、单元与支座均通过结点相连。确定了单元的全部结点，也就确定了单元的划分。 在杆系结构中，杆件的连结点、截面突变点、支座结点以及荷载作用点等均可取作结点。单元编号：①，②，③，……， m1234③ ① ② 5② 2③ ④ 47① 4352⑥ ③ ④ ⑤ 6② ⑤ 17613① ⑥ 第2页/共84页结点编号：1，2，3，……，n曲杆可用多段折线单元代替。 进行结点编号时，要尽量使单元两端结点编号的差值最小。 单元局部坐标系为 ，单元始端指向末端的方向就是 轴的正方向，由 轴正方向顺时针转90。得到 轴的正方向，如下图所示。x1e坐标系2y第3页/共84页三、单元杆端力和杆端位移的坐标变换1．坐标系 结构整体坐标系为xy，坐标轴x、y 遵循右手法则，即从x轴正方向顺时针转 90°得到y轴正方向。12e第4页/共84页2. 单元杆端力和杆端位移1）连续梁单元连续梁单元的杆端没有线位移。xx1122yy单元杆端力单元杆端位移ee第5页/共84页2）平面刚架单元第6页/共84页单元局部坐标系结构整体坐标系 以上杆端力和杆端位移中，下标1表示单元的始端，下标2表示单元的末端。在单元局部坐标系中， ；位移对单元杆端力无贡献。第7页/共84页3）桁架单元单元局部坐标系结构整体坐标系第8页/共84页4）单元杆端力和杆端位移的正负号 以上杆端力和杆端线位移与相应的坐标轴正方向一致为正，相反为负。 以上杆端力矩和杆端转角均以顺时针方向为正，逆时针方向为负。 角：由整体坐标系的x轴顺时针转到单元局部坐标系 轴的正方向所形成的夹角,如下图示。x1e2y第9页/共84页3. 单元坐标转换矩阵x1e2y第10页/共84页 对于平面刚架单元，用整体坐标系中的杆端力来表示单元局部坐标系中的杆端力，得到：第11页/共84页写成矩阵形式即第12页/共84页[T ]称为单元坐标转换矩阵。同理有 。第13页/共84页 对于平面桁架单元，其单元坐标转换矩阵为： 坐标转换矩阵[T]是正交矩阵，即矩阵任意两行或两列元素对应乘积之和等于零。 如果用局部坐标系中的单元杆端力来表示整体坐标系中的杆端力，就可以得出：同理 单元局部坐标系中， 。 ① 结构整体坐标系中， 。 ② 在上面两式中， 分别称为局部坐标系及整体坐标系中的单元刚度矩阵。 通常总是先求出局部坐标系中的单元刚度矩阵 ，然后利用 推得整体坐标系中的单元刚度矩阵 。第14页/共84页四、单元刚度矩阵 单元杆端力与单元杆端位移之间的转换关系式称为单元刚度方程。将①代入 得到： 因为 ，于是得到：第15页/共84页下面讨论 和 之间的转换关系。③④比较式②、④，有：⑤ 下面讨论如何求局部坐标系中的单元刚度矩阵。 为了推导 中各元素，采用单位位移法：即在单元的六个杆端位移中，每次只令一个杆端位移为1，其余杆端位移为0。为此，在单元两端就必须施加一组杆端力，这一组杆端力就构成了 的一列元素。第16页/共84页1．平面刚架单元12EA le12EI le1EI l2e第17页/共84页12EI leEI l12e12EI le第18页/共84页第19页/共84页 由上面六个图中的杆端力就可以写出局部坐标系中单元杆端力和杆端位移之间的刚度方程：第20页/共84页 上式就是平面刚架单元在局部坐标系中的单元刚度矩阵。 中第2 列各元素是 而其余杆端位移等于零时单元的六个杆端力的值。 中第2行各元素是单元的六个杆端位移分别等于1时杆端力 的值。第21页/共84页单元刚度矩阵中元素的物理意义：第22页/共84页2．平面桁架单元EI l21e令 ，则e12EI l第23页/共84页3 ．连续梁单元4． 的性质1）对称性，即 。2）奇异性，即 。 对于平面刚架单元和桁架单元，因为单元两端无约束，在平面内