结构力学第10矩阵位移法d.pptVIP

结构力学第10章 矩阵位移法 主要内容 * * 1 基本概念 2 局部坐标系下单元刚度矩阵 3 整体坐标系下单元刚度矩阵 4 直接刚度法 5 引入边界条件的方法 6 等效结点荷载 7 直接刚度法的另一种形式——先处理法 8 用先处理法计算矩形刚架 §10.6等效结点荷载 在前面的讨论中，结构仅承受结点荷载，然而，在实际结构中，经常会遇见荷载不作用于结点的情况，即非结点荷载问题。在进行结构分析时，必须将非结点荷载转化为等效结点荷载。 等效的原则：结构在结点荷载作用下，其结构的结点位移与原荷载作用下的结点位移相同。这种经过变换后所得到的结点荷载称为等效结点荷载。 2 图10-10平面刚架 ① ② ③ 1 3 4 q Fp l/2 l/2 l （a） q Fp M1=-Fpl/8 M2=Fpl/8-ql2/12 M3=ql2/12 Fp/2 Fp/2 M1 M3 -ql/2 -ql/2 M2 （b） -ql/2 -ql/2 M3 M2 Fp/2 Fp/2 M1 + = 若结构同时承受结点和非结点荷载作用，此时荷载分为两类：一类为原直接作用于结点上的荷载——结点荷载{F}D；另一类为由非结点荷载转换来的等效结点荷载{F}E 。将这两部分荷载叠加即得综合结点荷载{F}= {F}D+ {F}E 。用综合结点荷载进行直接刚度法进行计算。 应当注意，在综合结点荷载作用下的杆端力，并不是实际荷载作用下的杆端力。只有综合结点荷载作用下的杆端力与相应的两端固定梁的杆端力相叠加后，才能得到结构在原荷载作用下的杆端力。 在先处理法中，首先将结构的结点位移分量统一编码（?1，?2，?3……?n）即总码，确定单元刚度矩阵中的元素在结构刚度矩阵中的位置是通过单元定位向量来实现的。 §10.7直接刚度法的另一种形式——先处理法 前面介绍了直接刚度法的后处理法。它的优点是简明易懂，易为初学者掌握，且适用于各种边界条件。但由于在结构的原始刚度方程中包含了所有结点（包括了支座结点）的平衡方程，因此原始刚度矩阵的阶数较高，占机内存较多。下面介绍直接刚度法的先处理法可以解决这个问题。 单元定位向量是指由结点位移分量编码（总码）组成的表示单元杆端位移编码（局部码）的向量，由{?}e表示。 i1 i2 i3 M1 M2 M3 图10-11 三跨连续梁 ?2 ?3 ?1 其位移法典型方程为 其中 i1 i2 i3 M1 M2 M3 图10-11 三跨连续梁 ?2 ?3 ?1 采用相同的结点位移分量，用直接刚度法的先处理法来组装结构的刚度矩阵看看是否相同。 (1)结构离散化：2个单元，4个结点 ① ② 1 2 3 4 ③ 结点位移分量统一编码，?1，?2，?3，结点位移分量已知的均变为0码。 ?1 ?2 ?3 (2)结点位移列向量： 显然结构刚度矩阵为3?3方阵。 (3)单元分析，并确定单元刚度矩阵中的元素在结构刚度矩阵中的位置 单元① i=1，j=2，{?}1={i j}T={1 2}T 单元② i=2，j=3，{?}1={i j}T={2 3}T ① ② 1 2 3 4 ③ ?1 ?2 ?3 单元③ i=3，j=0，{?}1={i j}T={3 0}T (4)组装结构刚度矩阵 赋0 加单元① 加单元② 加单元③ 注意： (1)用总码表示单元刚度矩阵中的元素在结构刚度矩阵中的位置时，若其下标有一个为0，则在结构刚度矩阵中不出现。即结构刚度矩阵不受此类元素大小的影响。 (2)求出结构刚度矩阵后，即得结构刚度方程 [K]{?}={F} 由于前处理法中，在形成结构刚度矩阵中已考虑了结构的边 界条件，因此矩阵是非奇异的，因此 例4 如图示刚架，已知EA=1.5?106kN，EI=1.25?105kN.m2，作M图。 5m 5m 60? 图10-12 例4题图 5kN.m q=4.8kN/m 6kN 解： 图（a） x y ① ② 1 2 3 ?1 ?3 ?2 (1) 结构离散化，2个单元，3个结点，整体坐标系如图示。 (2)结构的位移列向量 {?} ={?1 ?2 ?3 }T 显然结构的刚度矩阵为3?3方阵 结点荷载列向量为3?1列阵 {F} ={F1 F2 F3 }T。 (3)单元分析 {?}1={1 2 3 4 5 6}T={1 2 3 0 0 0}T 单元① i=1，j=3，?=4?/3, 局部码 整体码 等效结点荷载为 单元② i=1，j=2，?=0 , l=5m {?}2={1 2 3 4 5 6}T={