0903矩阵位移法整体分析-先处理法（谷凡）.PPT

8.4 整体分析Global analysis 8.5 先处理法 后处理法的计算步骤 先处理法和后处理法比较 （1）单元、结点、结点位移分量编码；建立整体坐标系和各单元局部坐标系。 （2）建立结点位移列阵和结点力列阵 （3）建立整体坐标系下单元刚度矩阵 （4）建立结构刚度方程（5）求位移，计算杆端力 先处理法的计算步骤 例：求连续梁的整体刚度矩阵 例：先处理法建立结构刚度方程，杆长均为 l ，不考虑杆件轴向变形。 4. 结构刚度矩阵元素速算法 作业（1） 谢 谢 ! 2010.8 （4）建立结构刚度方程。 Kij 的物理意义： 当且仅当结构的第 j 个位移分量产生单位位移时，与第 i 个位移相应的“附加约束”上的反力大小。 速算法： 利用这一概念和等截面直杆的形常数，不经单元分析和集成而得到结构刚度矩阵中任何指定的元素。 Kij 的速算步骤： 1. 如果第 j 个位移和第 i 个位移分量间没有直接的单元连接（即非相关结点位移），则Kij＝0 ; 2. 如果第 j 个位移和第 i 个位移分量间是相关结点位移，则从结构中取出与?j 相关联的单元部分; 3. 令所取出部分结构仅产生?j=1的位移，作相应内力图； 4. 利用形常数和隔离体平衡条件求与单位位移相应的附加约束上的总反力，其值即为Kij。 例：求结构相应编码的结构刚度矩阵元素K44，K4,13，K66，K6,15. x y 1(0,0,0) 2(0,0,0) 3(0,0,0) 4(1,2,3) 5(4,5,6) 6(7,8,9) 7(10,11,12) 8(13,14,15) ? ? x y 1(0,0,0) 2(0,0,0) 3(0,0,0) 4(1,2,3) 5(4,5,6) 6(7,8,9) 7(10,11,12) 8(13,14,15) 例：求结构相应编码的结构刚度矩阵元素K11，K27，K15，K55 ，K35 。(本题特征是不计轴向变形) x y 1(0,0,0) 2(0,0,0) 3(0,0,0) 4(1,0,3) 5(1,0,4) 6(1,0,5) 7(2,0,6) 8(2,0,7) 6(1,0,5) x y 1(0,0,0) 2(0,0,0) 3(0,0,0) 4(1,0,3) 5(1,0,4) 7(2,0,6) 8(2,0,7) 例：求结构相应编码的结构刚度矩阵元素K11，K27，K15，K55 ，K35 。(本题特征是不计轴向变形) 8-1. 用先处理法计算图示结构刚度矩阵的元素K22、K33、K13、K15 。 * 第八章 矩阵位移法 第八章 矩阵位移法 结点、单元标码，并选择整体坐标系和局部坐标系； 结点位移分量编码，建立整体坐标系下的结点位移列阵和结点力列阵； 建立局部坐标系下单元刚度矩阵，坐标变换，建立整体坐标系下单元刚度矩阵； 利用定位向量形成结构原始刚度矩阵； 引入位移边界条件修改原始刚度矩阵，计算未知结点位移； 计算整体坐标系、局部坐标系下的单元杆端力。 1. 基本未知量 2. 采用的单元 3. 占用计算机存储空间 各有优缺点。 1. 准备工作 （1）结点力列阵和结点位移列阵 FR2y FR1y FR1x F2y F2x F1x F1y F3y F4y F4x F3x 1(0,0) 2(1,0) 3(2,3) 4(4,5) ① ② ③ ④ ⑤ ⑥ 未知的结点位移分量所对应的结点力是已知的； 已知的结点位移分量所对应的结点力是未知的。 FR2y FR1y FR1x F2y F2x F1x F1y F3y F4y F4x F3x 1(0,0) 2(1,0) 3(2,3) 4(4,5) ① ② ③ ④ ⑤ ⑥ FR2y FR1y FR1x F2y F2x F1x F1y F3y F4y F4x F3x 1(1,2) 2(3,4) 3(5,6) 4(7,8) ① ② ③ ④ ⑤ ⑥ 1（1,2） 3（5,6） 4（7,8） 2（3,4） 1（0,0） 3（2,3） 4（4,5） 2（1,0） （2）单元定位向量 由单元结点位移的整体码（总码）按与单元局部码对应的顺序组成的向量。 1(0,0) 2(1,0) 3(2,3) 4(4,5) ① ② ③ ④ ⑤ ⑥ ① ② ③ A B C D FP1 FP2 O y x 1(0,0,0) 2(1,2,3) 3(4,5,6) 4(0,0,7) ① ② ③ FP1 FP2 O y x 1(0,0,0) 2(1,2,3) 3(4,5,6) 4(0,0,7) ① ② ③ 1 2 3 4 5 6 位移局码 1 2 3 4