有限元第四章-杆系结构的整体分析.pptVIP

第4章　杆系结构的整体分析 坐标变换 由局部坐标下的单元刚度矩阵和单元等效节点荷载矩阵，通过坐标转换形成整体坐标下的单元刚度矩阵和等效荷载矩阵。 集成规则 用最小势能原理（或虚位移原理）建立结构的的整体刚度方程，并有次导出直接刚度法的集成规则。 约束处理 采用先处理法或后处理法进行边界约束处理，使结构的边界已知位移得以自动满足。 解方程 选取与整体（结构）刚度矩阵存贮方法相应的线性方程组，求得结构的结点位移。 求内力 从结构全部结点位移中取出单元的结点位移，经坐标转换变成局部坐标下杆端位移，在由单元刚度方程求得杆端内力，进一步还可求出反力和截面内力，为结构设计提供数据。 4-1坐标转换 在讲结构离散时，需要建立二套坐标杆系：杆轴为 另两轴为惯性主轴的局部坐标。整个结构统一的笛卡尔坐标。上一章单元分析是在局部坐标下进行的，而实际结构的中的每个杆件（也即单元）方位除了连续梁之外各不相同，要考虑结点位移协调、受力平衡，很自然应该用统一的结构整体坐标系。显然，二套坐标下对应的物理量必然存在相互转换关系，在进行具体整体分析之前应该将局部的量转换成整体的量，或将整体的量转换成局部的量。这项工作称为坐标转换。 附：杆系结构分析算例 4-5-1 “划零置一”法 设结点位移向量中第 个位移 已知，对应整体刚度方程有如下修改。 （1）将 改为 （2）整体刚度矩阵 中的第 行第 列中，将主对角线元素 改为1，其他元素改为0。 经修改后第 个方程就改为 ，满足了边界条件，同时也考虑了他对其他方程的影响 4-5-2 乘大数法 仍设结点位移向量中第 个位移 已知，应对第 个方程进行如下修改。 在 和 中将对角线元素 改为 GKrr ，将 改为 ，其中 G为一大数通常为106~108 （双精度是也可更大），经过修改后，第 个方程变为 在上式等式两边同除以 GKrr，由于G是很大的数，故可得 即 对于后处理法建立的刚度方程，经过上述修改后在整个刚度方程完全满足了结构在支座处的位移边界条件，如果边界位移条件足以消除结构的刚性位移，则修改后的刚度矩阵改称为整体刚度矩阵，此时他已是非奇异的对称矩阵，通过求解结构刚度方程，即可求得未知数的结点位移。 4-5-3 斜支撑处理 （ 略） 附;　边界条件的处理 　 　1　铰结点 　　在杆系结构中，除了刚性结点外，通常会遇到一些杆件通过铰结点与其它杆件联结，通过铰支座与其它杆件铰接和刚性接触。对于这样的铰结点，具有如下的性质： 　　①　铰结点上各杆具有相同的线位移，但截面的转角位移不相同； 　　②　结点上具有铰接杆端不承受弯矩作用。 　　对于这样的结点，我们在对其进行单元划分时，通常考虑在D处设置2个结点。按照先处理法，对图示结构进行位移编码，如图4-10所示。 图4-10　铰结点的处理示意图 　　2.4.2.2　弹性支承点 　　在实际工程中，有时会遇到弹性支承的情况如下图所示，这时一般将弹性支座看作是在结构约束点沿约束方向的一个弹簧，弹簧的刚度系数为 ， 在数值上等于使弹簧支座沿约束方向产生单位位移时所需施加的力。 （a）结构离散示意 ③ ② ① ① ② ③ （b）单元结点受力 图 4-6 结点外力势计算图 显然式（a）结论适合于任何情况。由此式（4-26）可具体写为 例如图4-6所示结构结点的外力势为 式（a）中: （a） （b） （4-28） 若引入如下矩阵符号 （4-31） （4-30） （4-29） （4-32） 式（4-31）中，diag表示“元素”对角线排列，也即为块对角矩阵。 根据矩阵乘法规则，不难验证式（4-28）可改写为： （4-33） 2、位移转换 由单元杆端位移和结点位移之间的协调条件，可建立 ， 之间的对应关系 式中: 为反映位移协调的位移变换矩阵。 例如图4-6所示结构，其 ， 为 由此可得位移变换矩阵 为 式中 为单位矩阵。 （4-34） 3、整体刚度方程 将式（4-3