有限元方法 第二章 矩阵位移法概述.ppt

第二章 矩阵位移法概述 本部分内容结束 集成总刚阵时必须使用整体坐标表示的单元刚度矩阵。 用先处理法集成总刚阵时必须先建立各单元的定位向量λ根据定位向量的指引将单刚阵中的元素逐个搬入总刚度矩阵中。 用后处理法集成总刚阵时必须先集成原始总刚度矩阵。集成原始总刚阵时应根据结点编号情况指示矩阵G 以子块搬家的方式将单刚阵中的子块逐个搬入总刚度矩阵中。 如不考虑轴向变形的单元 由6×6刚度矩阵划去1、4行和列后可得 平面刚架程序的扩大功能: 1. 平面桁架 2. 桁梁组合体系 3. 斜向支座 4. 弹性支座 5. 弹性地基 6. 带铰结点的刚架 矩阵位移法与位移法在理论上并无区别，只是在表达方式上有所不同。 (1)矩阵位移法的理论基础与一般位移法完全相同，只是表达方式不同。用矩阵形式表示具有更强的概括性。 (2)总刚度矩阵是由各单元刚度矩阵装配成的，只要找出了装配的规律，总刚度矩阵不必计算而可直接由单元刚度矩阵装配而成。 (3)矩阵位移法与一般位移法解题步骤的对应关系可以由下表表示： §2-9 总结 任一结点i上的等效结点荷载FEi为： 如果除了非结点荷载的等效结点荷载FEi外，结点i上还作用有直接结点荷载FDi，则i点总的结点荷载为： 结点i的综合结点荷载 整个结构的综合结点荷载 各单元最后的杆端力是固端力和综合结点荷载作用下产生的杆端力之和，即 和 或 表 : 单元固端约束力 （局部坐标系） ? 荷载简图 ? 始 端 1 末 端 2 1 2 12 2 ab F P - 表 : 单元固端约束力 （局部坐标系） 3 4 表 : 单元固端约束力 （局部坐标系） 5 6 表 : 单元固端约束力 （局部坐标系） 7 计算步骤： （1）对结点和单元进行编号，选定整体坐标系和局部坐标系； （2）计算各杆的单元刚度矩阵； （3）形成结构原始刚度矩阵； （4）计算固端力、等效结点荷载和综合结点荷载； （5）引入支承条件，修改结构原始刚度方程，得到缩减总刚； （6）结算结构刚度方程，求出结点位移； （7）计算各单元杆端力。 §2-7 矩阵位移法的计算步骤和示例 [K] 求单元常数 [T] {F} 原始数据、局部码、总码 解方程{F}=[K]{?} 求出结点位移 {?} 开始 单元刚度 矩阵 k e 单元固 端力 e 结束 [K]{?} = {F} {FP} + = 程序设计框图 求杆端力 e e e e 1 2 3 4 30KN/m 100KN 50KN 2m 2m 4m x y ② ① ③ 10-1 求图示刚架的内力。已知各杆材料及截面相同。 （1）将单元、结点编号，确定坐标系， 如图所示。 （2）求出各单元在整体坐标系中的单 元刚度矩阵，见书P247。 （3）将各单刚子块对号入座，形成结构原始刚度矩阵，见书P248。 （4）计算非结点荷载作用下的各单元固端力、等效结点荷载及综 合结点荷载。 对局部坐标和整体坐标不一致的单元，要对刚度、荷载进行坐标转换。 各单元在其局部坐标系下的固端力为： 经过坐标转换，得到各单元在整体坐标系下的固端力为： 结点2、3上的等效结点荷载为： 结点2、3上的综合结点荷载为： 1 2 3 4 30KN/m 100KN 50KN 2m 2m 4m x y ② ① ③ 结构的结点外力列向量为 这里，F1和F4应为综合结点荷载和支座反力的代数和，其中支座反力仍为未知量；引入支承条件时， F1和F4将被划掉，因此不必计算其等效结点荷载和综合结点荷载。 结构原始刚度方程见书P256 （5）引入支承条件，修改原始刚度方程。 结点1、4为固定端，位移已知： 代入原始刚度方程，得到修改后的结构刚度方程为 （6）解方程，求得未知结点位移为： （7）计算各单元杆端力，见书P258。 10-2 平面桁架如图所示，各杆截面EA均为常数。已知F1=15 kN, F2=20 kN, 试求桁架各杆轴力。 解: (1) 对结点和单元编号 cosα sinα α 即Cy 即Cx 单 元 单元坐标x轴方向（即指示矩阵） 单元长度 (m) 1→2 ① 30° 1 0 1.732 ② 3→2 0° 0.866 0.5 2.00 各单元参数表 (2) 列出各单元刚度矩阵 整体坐标表示的单元刚度矩阵 即 (3) 集合总刚度矩阵 原始总刚度