矩阵位移法讲稿.ppt

《矩阵位移法讲稿.ppt》由会员分享，可在线阅读，更多相关《矩阵位移法讲稿.ppt（50页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于矩阵位移法1第一页，讲稿共五十页哦2 边界条件可分为：边界条件可分为：前处理法前处理法后处理法后处理法8-4 结构的整体分析第二页，讲稿共五十页哦3后处理法后处理法单元采用自由式单元刚度矩阵，形成结构的单元采用自由式单元刚度矩阵，形成结构的原原始刚度矩阵始刚度矩阵对原始刚度矩阵进行边界条件处理，形成结对原始刚度矩阵进行边界条件处理，形成结构刚度矩阵构刚度矩阵 K 原始刚度矩阵的阶数由结点总数乘结点的位移原始刚度矩阵的阶数由结点总数乘结点的位移分量来确定，便于编制通用程序。分量来确定，便于编制通用程序。8-4 结构的整体分析第三页，讲稿共五十页哦4前处理法前处理法形成结构刚度方程之前处理边界

2、条件形成结构刚度方程之前处理边界条件结构的结点位移分量只有独立的未知位移分量，由结构的结点位移分量只有独立的未知位移分量，由单元刚度矩阵直接形成结构刚度矩阵。单元刚度矩阵直接形成结构刚度矩阵。在程序编制中，须建立结点位移分量编号数组，在程序编制中，须建立结点位移分量编号数组，来代替后处理法的约束处理数组。来代替后处理法的约束处理数组。8-4 结构的整体分析第四页，讲稿共五十页哦5后处理法后处理法 四个刚结点，共有四个刚结点，共有12个结点位移分量，结构的结个结点位移分量，结构的结点位移列向量为点位移列向量为8-4 结构的整体分析第五页，讲稿共五十页哦6后处理法后处理法设刚架上只有结点荷载作用，

3、结点外力（包括设刚架上只有结点荷载作用，结点外力（包括荷载和反力）列向量为荷载和反力）列向量为形成结构的原始刚度矩阵形成结构的原始刚度矩阵 K结构的原始刚度矩阵为结构的原始刚度矩阵为12阶方阵。阶方阵。8-4 结构的整体分析第六页，讲稿共五十页哦7用后处理法组集图示结构的原始刚度矩阵。用后处理法组集图示结构的原始刚度矩阵。8-4 结构的整体分析第七页，讲稿共五十页哦8结构的原始刚度方程结构的原始刚度方程K8-4 结构的整体分析第八页，讲稿共五十页哦9将单元杆端位移按结点将单元杆端位移按结点i、j划分为两个子阵划分为两个子阵 8-4 结构的整体分析第九页，讲稿共五十页哦10 单元刚度矩阵的四个子

4、块分别为单元刚度矩阵的四个子块分别为 1 3 3 4 2 4结点码结点码各单元的各单元的始、末两端始、末两端i、j 的结的结点号码。点号码。8-4 结构的整体分析第十页，讲稿共五十页哦11结构原始刚度矩阵结构原始刚度矩阵K的组集的组集把每个单元刚度矩阵的四个子块按其两个下把每个单元刚度矩阵的四个子块按其两个下标号码，送到结构原始刚度矩阵中相应的行标号码，送到结构原始刚度矩阵中相应的行和列的位置上去，可得到结构原始刚度矩阵。和列的位置上去，可得到结构原始刚度矩阵。即即“对号入座对号入座”原则。原则。8-4 结构的整体分析第十一页，讲稿共五十页哦12 各单元刚度矩阵的四个子块分别为各单元刚度矩阵的

5、四个子块分别为 1 3 3 4 2 48-4 结构的整体分析第十二页，讲稿共五十页哦1321344321结结点点码码8-4 结构的整体分析第十三页，讲稿共五十页哦14结构的原始刚度方程结构的原始刚度方程“原始原始”表示未进行支承条件处理。表示未进行支承条件处理。8-4 结构的整体分析第十四页，讲稿共五十页哦15 结构原始刚度矩阵结构原始刚度矩阵 K 中中元素元素Kij 的物理意义：的物理意义：当仅发生广义位移当仅发生广义位移j1时，在第时，在第i个广义位移个广义位移对应处所需施加的广义力。对应处所需施加的广义力。K87？K26？8-4 结构的整体分析第十五页，讲稿共五十页哦16K87K878-

6、4 结构的整体分析第十六页，讲稿共五十页哦17原始刚度矩阵的性质：原始刚度矩阵的性质：1.对称性：对称性：KijKji，K 是对称矩阵。是对称矩阵。2.奇异性：奇异性：结构存在刚体位移，在给定的外荷结构存在刚体位移，在给定的外荷载作用下不能确定结构的位移，载作用下不能确定结构的位移，K 是是奇异的。奇异的。8-4 结构的整体分析第十七页，讲稿共五十页哦183.稀疏性稀疏性相关结点和不相关结点。相关结点和不相关结点。i 和和 j 不相关，则不相关，则 K 的子矩阵的子矩阵KijKji=0。结构刚度矩阵结构刚度矩阵 K 只在对角线附近一带状区只在对角线附近一带状区域内有非零元素，域内有非零元素，具

7、有稀疏性。具有稀疏性。8-4 结构的整体分析第十八页，讲稿共五十页哦19 二二.边界条件处理边界条件处理常用的三种边界条件处理方法：常用的三种边界条件处理方法：1划行划列法划行划列法 2乘大数法乘大数法 3置换法置换法 8-4 结构的整体分析第十九页，讲稿共五十页哦20划行划列法划行划列法将方程分块将方程分块 将将 中的行列重新排列，把受中的行列重新排列，把受约束的结点位移分量靠后。约束的结点位移分量靠后。8-4 结构的整体分析第二十页，讲稿共五十页哦21可写成可写成由由第第一一式式求求解解未未知知结结点点位位移移 ：由第二式求解未知结点力由第二式求解未知结点力Fr ，可求得支座反力，可求得支

8、座反力FR8-4 结构的整体分析第二十一页，讲稿共五十页哦22当支座位移为零时，简化为当支座位移为零时，简化为 相当于把原始刚度方程中对应于已知支座位移相当于把原始刚度方程中对应于已知支座位移分量为零的行与列划去，因此称为划行划列法。分量为零的行与列划去，因此称为划行划列法。8-4 结构的整体分析第二十二页，讲稿共五十页哦23乘大数法乘大数法设已知设已知 ，进行替换，进行替换8-4 结构的整体分析第二十三页，讲稿共五十页哦24第第i个方程为个方程为两边同除以两边同除以N，得，得即即 8-4 结构的整体分析第二十四页，讲稿共五十页哦25一般情况下一般情况下此时只需替换此时只需替换kii8-4 结

9、构的整体分析第二十五页，讲稿共五十页哦26此时，第此时，第i个方程为个方程为两边同除以两边同除以N，得，得即即 8-4 结构的整体分析第二十六页，讲稿共五十页哦27置换法置换法设设 (1)将将 Fj 换为换为 (2)将将i行、行、j列非对角线元素换为列非对角线元素换为0；(3)将将kii换为换为1，Fi换为换为Ci 能满足相应的位移条件。能满足相应的位移条件。8-4 结构的整体分析第二十七页，讲稿共五十页哦28前处理法前处理法用前处理法分析结构时，结构的结点位移分用前处理法分析结构时，结构的结点位移分量只引入独立的未知位移分量。量只引入独立的未知位移分量。8-4 结构的整体分析第二十八页，讲稿

10、共五十页哦29与结点位移列向量相对应的结点外力（结点力向与结点位移列向量相对应的结点外力（结点力向量不包括支座反力）列向量为量不包括支座反力）列向量为 8-4 结构的整体分析第二十九页，讲稿共五十页哦30按前处理法对结点位移编码时，将已知为零的按前处理法对结点位移编码时，将已知为零的结点位移分量编号均用零表示。结点位移分量编号均用零表示。8-4 结构的整体分析第三十页，讲稿共五十页哦31结构的结点位移分量只引入未知位移分量；对结构的结点位移分量只引入未知位移分量；对应的结点力向量不包括支座反力应的结点力向量不包括支座反力8-4 结构的整体分析第三十一页，讲稿共五十页哦32由单元刚度矩阵直接形成

11、考虑边界条件由单元刚度矩阵直接形成考虑边界条件的结构刚度方程。的结构刚度方程。K8-4 结构的整体分析第三十二页，讲稿共五十页哦33(1).(1).结点位移分量编码及单元的定位向量结点位移分量编码及单元的定位向量 由单元杆端位移分量对应的结构结点由单元杆端位移分量对应的结构结点位移分量序号所组成的向量，称为单元的位移分量序号所组成的向量，称为单元的定位向量。定位向量。8-4 结构的整体分析第三十三页，讲稿共五十页哦34单元的定位向量为：单元的定位向量为：0 0 0 1 2 30 0 0 1 2 3结点位移码结点位移码8-4 结构的整体分析第三十四页，讲稿共五十页哦35单元的定位向量为：单元的定

12、位向量为：1 2 3 4 5 61 2 3 4 5 6结点位移码结点位移码8-4 结构的整体分析第三十五页，讲稿共五十页哦36单元的定位向量为：单元的定位向量为：0 0 0 4 5 60 0 0 4 5 6结点位移码结点位移码8-4 结构的整体分析第三十六页，讲稿共五十页哦37（2）用单元定位向量确定单元刚度矩阵元素）用单元定位向量确定单元刚度矩阵元素在结构刚度矩阵中位置：在结构刚度矩阵中位置：1.将单元的定位向量分别写在单元刚度将单元的定位向量分别写在单元刚度矩阵的上方和右侧。矩阵的上方和右侧。2.若单元定位向量的某个分量为零，若单元定位向量的某个分量为零，把把 中相应的行和列删去，不送入结

13、构刚度矩中相应的行和列删去，不送入结构刚度矩阵阵 K。8-4 结构的整体分析第三十七页，讲稿共五十页哦38 3.单元定位向量就是单元定位向量就是 中元素中元素 在结构刚度在结构刚度矩阵矩阵 K 中的行码和列码。中的行码和列码。按照单元定位向量中非零分量给出的行码按照单元定位向量中非零分量给出的行码和列码，就能够将单元刚度矩阵的元素正和列码，就能够将单元刚度矩阵的元素正确地累加到结构刚度矩阵确地累加到结构刚度矩阵K 中去。中去。8-4 结构的整体分析第三十八页，讲稿共五十页哦39单元的单元刚度矩阵对单元的单元刚度矩阵对K的贡献的贡献 0 0 0 1 2 30 0 0 1 2 3结点位移码结点位移

14、码遇零不送遇零不送8-4 结构的整体分析第三十九页，讲稿共五十页哦40单元的单元刚度矩阵对单元的单元刚度矩阵对K的贡献的贡献 0 0 0 4 5 60 0 0 4 5 6结点位移码结点位移码8-4 结构的整体分析第四十页，讲稿共五十页哦41 1 2 3 4 5 61 2 3 4 5 68-4 结构的整体分析第四十一页，讲稿共五十页哦42结构刚度矩阵结构刚度矩阵K8-4 结构的整体分析第四十二页，讲稿共五十页哦43主对角线元素是由同一结点相关单元主对角线元素是由同一结点相关单元的刚度矩阵主对角线元素叠加而成，的刚度矩阵主对角线元素叠加而成，一定为正。一定为正。副对角线元素是由定位向量所对应的单元

15、副对角线元素是由定位向量所对应的单元刚度矩阵副对角线元素累加而成，可为正、刚度矩阵副对角线元素累加而成，可为正、负或零值。负或零值。8-4 结构的整体分析第四十三页，讲稿共五十页哦44“对号入座对号入座”法综述法综述把单元杆端位移分量序号换成对应的结点位移分把单元杆端位移分量序号换成对应的结点位移分量序号，通常称为量序号，通常称为换码，换码，搬到搬到K 中相应编号的位中相应编号的位置，通称置，通称对号入座对号入座。K 中同一号码位置上的元素相加，称为中同一号码位置上的元素相加，称为集合集合。“换码换码”的实质是满足变形协调条件。的实质是满足变形协调条件。“集合集合”的实的实质是满足平衡条件。质

16、是满足平衡条件。8-4 结构的整体分析第四十四页，讲稿共五十页哦45 例：试求结构原始刚度矩阵中的子块试求结构原始刚度矩阵中的子块 的的4个元素个元素。8-4 结构的整体分析第四十五页，讲稿共五十页哦46已知各杆件在整体坐标系中的单元刚度矩阵为已知各杆件在整体坐标系中的单元刚度矩阵为：8-4 结构的整体分析第四十六页，讲稿共五十页哦47例：例：求图示结构原始刚度矩阵的元素求图示结构原始刚度矩阵的元素 。8-4 结构的整体分析第四十七页，讲稿共五十页哦48已知已知、单元单元整体坐标系整体坐标系的单元刚度矩阵的单元刚度矩阵 为为 1 3 4 562 12 3 4 568-4 结构的整体分析第四十八页，讲稿共五十页哦49 8-4 结构的整体分析第四十九页，讲稿共五十页哦感感谢谢大大家家观观看看第五十页，讲稿共五十页哦