弹性力学与有限元分析第三章平面刚架有限元分析及程序设计.ppt

§3.6 程序设计 §3.6.1 程序框图 输入数据 单元局 部荷载 坐标转换矩阵 单元整体荷载 集成整体刚 度矩阵元素 约束条件处理、解方程 计算单元轴力、约束反力 单元循环 包括单元、结点、材料、荷载、约束数据 单元局 部刚度 单元整体刚度 集成整体等效结点荷载向量元素 形成整体刚度方程 §3.6 程序设计 §3.6.2 程序说明 * * 注意这里的弯矩是作为外力的，符号应与外力符号一致。但作为内力时（双重作用,即作为杆端力时是外力，又作为杆件的内力）。 本质上是建立杆端位移引起的荷载与外荷载引起的杆端力之间的结点力的平衡关系，所以关键是找到杆端位移引起的荷载与外荷载引起的杆端力。 讲清楚位移法的基本思路（以位移表示的平衡方程）。 系数对应着位移编号。 单位位移。 注意这里的正负号（杆端力-外力-与坐标是否一致） 对于连续梁单元，方程左端为左右端弯矩，右端的位移为左右端的转角。 4 解整体刚度方程，求结点位移 §3.5 计算步骤和算例 5 求各单元的杆端内力 (1) 确定整体坐标系下的杆端位移 ① ② ③ ① ② ③ ① ② 2 1 A B C D ③ 4 3 5 §3.5 计算步骤和算例 ① ① ① ② ② ③ ③ ③ (2) 确定局部坐标系下的杆端位移 §3.5 计算步骤和算例 5 求各单元的杆端内力 ① ① ① ① §3.5 计算步骤和算例 ② ② ② ② ③ ③ ③ ③ §3.5 计算步骤和算例 ③ ② ① 1kN/m ① 1.24kN 2.09kNm 4.76kN 8.49kNm 0.43kN 0.43kN 0.43kN 0.43kN 2.09kNm 3.04kNm ② 1.24kN 1.24kN ③ 1.24kN 3.04kNm 1.24kN 4.38kNm 0.43kN 0.43kN §3.5 计算步骤和算例 6 画出内力图 2.09 8.49 4.38 3.04 M图(kNm) 1.24 4.76 1.24 0.43 FQ图(kN) 0.43 0.43 1.24 FN图(kN) §3.5 计算步骤和算例 C1 A B D C2 ① ② ③ 1 4 1 2 3 ② ③ ① 例题 求整体刚度矩阵. 已知整体坐标系的单元刚度矩阵. §3.5.3 忽略轴向变形刚架的分析 §3.5 计算步骤和算例 集成单元① 解 ① ① ② ③ 1 C1 A B D 4 C2 1 2 3 §3.5 计算步骤和算例 集成单元② ② ① ② ③ 1 C1 A B D 4 C2 1 2 3 §3.5 计算步骤和算例 集成单元③ ① ② ③ 1 C1 A B D 4 C2 1 2 3 ② §3.5 计算步骤和算例 整体刚度矩阵 §3.5 计算步骤和算例 例题 做图示刚架的内力图 12m 1kN/m 6m A B C D ① ② 2 1 A B C D ③ 1 3 解 1 确定整体和局部坐标系、单元和结点位移编码 §3.5 计算步骤和算例 2 形成刚度矩阵 (1)形成局部坐标系下的单元刚度矩阵。 ① ③ ② §3.5 计算步骤和算例 ① ③ (2)形成整体坐标系下的单元刚度矩阵。 ① ③ ① ① ① ② ② §3.5 计算步骤和算例 (3) 定位累加，形成整体坐标系下的整体刚度矩阵 集成单元① ① ② 2 1 A B C D ③ 1 3 ① §3.5 计算步骤和算例 集成单元② ① ② 2 1 A B C D ③ 1 3 ② §3.5 计算步骤和算例 集成单元③ ① ② 2 1 A B C D ③ 1 3 ③ §3.5 计算步骤和算例 整体刚度矩阵 3 形成等效结点荷载 (1)形成局部坐标系下的单元等效结点荷载向量 (3)定位累加形成整体坐标系下的等效结点荷载列向量 (2)形成整体坐标系下的单元单元等效结点荷载向量 ① §3.5 计算步骤和算例 ① ① ① ① 4 解整体刚度方程，求结点位移 §3.5 计算步骤和算例 5 求各单元的杆端内力 (1)由结点位移，确定各单元的局部坐标系下的杆端位移 ① ② ③ ① ② ③ §3.5 计算步骤和算例 ① ① ① ② ② ③ ③ ③ §3.5 计算步骤和算例 5 求各单元的杆端内力 ① ① ① ① §3.5 计算步骤和算例 ② ② ② ② ③ ③ ③ ③ §3.5 计算步骤和算例 ① ② ③ 1kN/m ① 1.25kN 2.09kNm 4.75kN 8.41kNm ③ 1.25kN 3.09kNm 1.75kN 4.47kNm 0.43kN 0.43kN 2.09kNm 3.09kNm ② 6 画出内力图 §3.5 计算步骤和算例 2.09 8.41 4.47 3.09 M图(kNm) 1.25 4.75 1.26 0.43 FQ图(kN) 0.43 0.43 1.25 FN图(kN)