第十章矩阵位移法 (2)精选文档.ppt
《第十章矩阵位移法 (2)精选文档.ppt》由会员分享,可在线阅读,更多相关《第十章矩阵位移法 (2)精选文档.ppt(84页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第十章矩阵位移法本讲稿第一页,共八十四页矩阵知识矩阵知识1 1、矩阵、矩阵:mnmn个数个数a aij ij 排成排成mm行、行、n n列的矩形阵列列的矩形阵列 记作记作方阵:方阵:行矢量:行矢量:本讲稿第二页,共八十四页列矢量:列矢量:2、行列式:、行列式:n阶方阵阶方阵A相应的行列式相应的行列式D,记作,记作若若D D0 0,A A为奇异矩阵为奇异矩阵本讲稿第三页,共八十四页3 3、矩阵运算、矩阵运算相等:相等:A Amnmn=B=Bmnmn,则,则a aijij=b=bijij加减:加减:C CmnmnA Amnmn+B+Bmnmn,则,则c cijij=a=aijij+b+bijij数
2、乘:数乘:C Cmnmn=k*A=k*Amnmn,则,则c cijij=k*a=k*aijij乘法:乘法:C Cmnmn=A=Amlml*B*Blnln,则,则转秩:转秩:B Bmnmn=A=AT Tmnmn,则,则b bijij=a=ajiji(A+B)(A+B)T T=A=AT T+B+BT T(kA)(kA)T T=kA=kAT T(AB)(AB)T T=B=BT T*A*AT T(反序定律)(反序定律)本讲稿第四页,共八十四页4、特殊矩阵、特殊矩阵单位矩阵单位矩阵对角矩阵对角矩阵对称矩阵:对称矩阵:A An*nn*n,a aijij=a=ajiji正定矩阵:正定矩阵:特征值都大于零的实
3、对称矩阵特征值都大于零的实对称矩阵充要条件:所有的主子式都大于零充要条件:所有的主子式都大于零即即|A|Ai i|0|0本讲稿第五页,共八十四页分块矩阵:分块矩阵:正交矩阵:正交矩阵:A AT T=A=A-1-15、逆矩阵、逆矩阵 B=AB=A-1-1AB=BA=IAB=BA=I(单位矩阵)(单位矩阵)A A、B B互为逆矩阵互为逆矩阵矩阵可逆,矩阵可逆,detAdetA0 0反序定律:反序定律:(AB)(AB)-1-1=B=B-1-1A A-1-1本讲稿第六页,共八十四页102 单元刚度矩阵(局部坐标系)单元刚度矩阵(局部坐标系)单元分析:杆端力杆端位移单元分析:杆端力杆端位移刚度关系刚度关
4、系(10101 1)转角位移方程)转角位移方程矩阵形式(轴向变形)矩阵形式(轴向变形)1 1、编码、编码 单元单元 e e(e=1,2(e=1,2,n)n)等截面直杆:单元等截面直杆:单元 I I、E E、A A、l l。2 2、局部坐标系、局部坐标系 ijij杆轴正方向杆轴正方向 (局部编码,箭头)(局部编码,箭头)ijexxyM,坐标系:(右手系)坐标系:(右手系)本讲稿第七页,共八十四页3、杆端力与杆端位移、杆端力与杆端位移 “”局部坐标标志局部坐标标志 正负号规则正负号规则 与坐标系对应与坐标系对应(列向量)(列向量)(10 3、4)F FSiSieF FNjNjM Mj jM Mi
5、iF FNiNiF FSjSjvieujjiuivj本讲稿第八页,共八十四页4、刚度方程、刚度方程杆端力与杆端位移的刚度关系杆端力与杆端位移的刚度关系(102)本讲稿第九页,共八十四页(表(表81)*本讲稿第十页,共八十四页(10 1)刚度方程刚度方程本讲稿第十一页,共八十四页刚度矩阵:刚度矩阵:行行数数杆端力杆端力列向量分量数列向量分量数 列列数数杆端位移杆端位移列向量分量数列向量分量数本讲稿第十二页,共八十四页记忆:记忆:小子块小子块 12-6-6-4 (主)(主)12-6-6-2 (副)(副)4、5 行、列,除主元素外,均为负值行、列,除主元素外,均为负值行行杆端力(杆端力(X、Y、M)
6、列列杆端位移(杆端位移(u、v、)排列顺序:排列顺序:ij刚度系数度系数 kij(物理意(物理意义)对应杆端位移杆端位移j1时,引起引起对应i方向的杆端力方向的杆端力Fi。本讲稿第十三页,共八十四页5、单元刚度矩阵的性质、单元刚度矩阵的性质(1)物理意义)物理意义)物理意义)物理意义单刚单刚单刚单刚表示单元杆端力与杆端位移的刚度关系表示单元杆端力与杆端位移的刚度关系表示单元杆端力与杆端位移的刚度关系表示单元杆端力与杆端位移的刚度关系单元刚度系数单元刚度系数单元刚度系数单元刚度系数当其所在列对应的杆端位移等于当其所在列对应的杆端位移等于当其所在列对应的杆端位移等于当其所在列对应的杆端位移等于1
7、1(其余杆端(其余杆端(其余杆端(其余杆端位移等于位移等于位移等于位移等于0 0)时,所引起的所在行对应的杆端力)时,所引起的所在行对应的杆端力)时,所引起的所在行对应的杆端力)时,所引起的所在行对应的杆端力(2 2)重要性质)重要性质)重要性质)重要性质对称性对称性kij ij=k=kji ji 奇异性奇异性奇异性奇异性|K|=0第第第第1 1和和和和2 2 行行行行(列列列列)与第与第与第与第4 4 或或或或5 5行行行行(列列列列)相加,所得一行相加,所得一行相加,所得一行相加,所得一行(列列列列)元素全为零元素全为零元素全为零元素全为零物理概念:已知杆端位移物理概念:已知杆端位移物理概
8、念:已知杆端位移物理概念:已知杆端位移杆端力,反之不成。因杆端力,反之不成。因杆端力,反之不成。因杆端力,反之不成。因为讨论为讨论的是自由式的是自由式的是自由式的是自由式单单元,元,元,元,存在任意的存在任意的存在任意的存在任意的刚刚性位移。性位移。性位移。性位移。分块性质分块性质本讲稿第十四页,共八十四页6、特殊单元、特殊单元 简支单元简支单元M1M212拉压杆单元拉压杆单元(108)便于坐标转换便于坐标转换(106)本讲稿第十五页,共八十四页弯曲杆单元弯曲杆单元(忽略轴向变形)(忽略轴向变形)本讲稿第十六页,共八十四页103 单元刚度矩阵的坐标转换单元刚度矩阵的坐标转换 整体分析整体分析统
9、一坐标系统一坐标系 变形条件、平衡条件变形条件、平衡条件整体坐标系整体坐标系 整体坐标系的单元刚度方程整体坐标系的单元刚度方程本讲稿第十七页,共八十四页 坐坐标变换标变换方法:方法:1、单单元坐元坐标转换标转换矩矩阵阵设设为为从从 x 轴轴到到 x 轴的夹角,逆时针为正y yx xy yx x转角与平面坐标系的变换无关转角与平面坐标系的变换无关本讲稿第十八页,共八十四页 力的转换:用力的转换:用X、Y、M表示表示 本讲稿第十九页,共八十四页 矩矩阵阵形式形式本讲稿第二十页,共八十四页 简简写:写:单元坐标转换矩阵单元坐标转换矩阵 T T正交矩阵正交矩阵 T T-1-1=T=TT T 即即TTT
10、TT T=T=TT TT=IT=I本讲稿第二十一页,共八十四页逆转换式逆转换式 同理可得杆端位移的转换关系同理可得杆端位移的转换关系本讲稿第二十二页,共八十四页2、整体坐标系中的单元刚度矩阵、整体坐标系中的单元刚度矩阵 由局部坐标系中的刚度方程代入整体坐标系中的杆端力、杆端位移的转换关系代入整体坐标系中的杆端力、杆端位移的转换关系等式的两边同乘等式的两边同乘TTT T写为:写为:本讲稿第二十三页,共八十四页整体坐标系中的刚度方程:整体坐标系中的刚度方程:同阶,相同的性质:同阶,相同的性质:k kijij整体坐标系中整体坐标系中j j1 1时在时在i i方向引起的杆端力方向引起的杆端力 对称性,
11、对称性,k kijijk kjiji 奇异性,自由式单元奇异性,自由式单元|k|k|0 0 分块性质分块性质整体坐标系的刚度矩阵:整体坐标系的刚度矩阵:写成分块形式写成分块形式:本讲稿第二十四页,共八十四页刚度矩阵计算公式刚度矩阵计算公式其中其中c=cosc=cos,s=sin,s=sin对称、奇异对称、奇异(10 2310 23)本讲稿第二十五页,共八十四页平面桁架杆件平面桁架杆件轴力轴力杆端力、杆端位移向量杆端力、杆端位移向量坐标转换矩阵坐标转换矩阵本讲稿第二十六页,共八十四页拉压杆单元单刚(整体坐标系)拉压杆单元单刚(整体坐标系)简支单元单刚简支单元单刚(整体坐标系与局部坐标系相同)(整
12、体坐标系与局部坐标系相同)本讲稿第二十七页,共八十四页【例例】连续梁连续梁连续连续梁的整体梁的整体刚刚度矩度矩阵阵整体整体刚刚度方程:度方程:FK 其中:其中:F结结点力向量点力向量 结结点位移向量点位移向量则则:K整体整体刚刚度矩度矩阵阵结结构的构的结结点力点力与与结结点位移点位移的的刚刚度关系度关系二跨二跨二跨二跨连续连续连续连续梁:梁:梁:梁:结结结结点位移向量点位移向量点位移向量点位移向量 =1 1 1 1 2 2 2 2 3 3 3 3 T T结结结结点力向量点力向量点力向量点力向量 F=F1 1 F F2 2 F F3 3 T T龙:图龙:图108F1F2123F3104 结构(原
13、始)刚度矩阵结构(原始)刚度矩阵本讲稿第二十八页,共八十四页位移法:加位移法:加约约束束按主按主动动位移考位移考虑虑:每个每个结结点位移点位移对对Fi的的贡贡献献迭加迭加刚刚度方程度方程整体整体刚刚度矩度矩阵阵位移法建立方程:位移法建立方程:连续连续梁梁 简简支支单单元元单刚单刚龙:图龙:图109F1F2123F3F1F21F3F1F2F32F1F23F3本讲稿第二十九页,共八十四页1.单单元集成法(直接元集成法(直接刚刚度法)度法)每个每个单单元元对对F的的贡贡献献单单元杆端位移所需要的力元杆端位移所需要的力令令i10叠加叠加单单元集成:元集成:KK贡献矩阵贡献矩阵令令i20F1F2123F
14、3F1F2123F3KK贡献矩阵贡献矩阵本讲稿第三十页,共八十四页2 2、定位向量法、定位向量法、定位向量法、定位向量法 定位向量定位向量定位向量定位向量 e确定位移局部编码与整体编码的对应关系确定位移局部编码与整体编码的对应关系确定位移局部编码与整体编码的对应关系确定位移局部编码与整体编码的对应关系局部码:单元杆端局部码:单元杆端局部码:单元杆端局部码:单元杆端 ijij整体码:对应结点编码整体码:对应结点编码整体码:对应结点编码整体码:对应结点编码 mnmne e单元杆端位移对应结点位移总码组成的向量单元杆端位移对应结点位移总码组成的向量 【例例例例】对号入座对号入座对号入座对号入座根据单
15、元定位向量将单刚元素根据单元定位向量将单刚元素kij ij对应放入总刚对应放入总刚对应放入总刚对应放入总刚 i i行、行、行、行、j j列位置列位置列位置列位置【例例例例】定位向量法定位向量法定位向量法定位向量法引入支座条件引入支座条件引入支座条件引入支座条件位移为零的结点位移分量,编码为零。位移为零的结点位移分量,编码为零。位移为零的结点位移分量,编码为零。位移为零的结点位移分量,编码为零。对应总码为零的单刚元素不置入整体刚度矩阵对应总码为零的单刚元素不置入整体刚度矩阵对应总码为零的单刚元素不置入整体刚度矩阵对应总码为零的单刚元素不置入整体刚度矩阵单刚元素,定位向量,对号入座,同号相加单刚元
16、素,定位向量,对号入座,同号相加单刚元素,定位向量,对号入座,同号相加单刚元素,定位向量,对号入座,同号相加【例例】本讲稿第三十一页,共八十四页连续梁整体刚度矩阵的性质连续梁整体刚度矩阵的性质整体刚度系数整体刚度系数Kij 当当当当j j j j1 1时(其它位移为零)所产生的对应时(其它位移为零)所产生的对应时(其它位移为零)所产生的对应时(其它位移为零)所产生的对应i i方向方向方向方向的结点力的结点力的结点力的结点力对称性:对称性:K Kij ijKji非奇异:非奇异:定位向量法,考虑支座约束定位向量法,考虑支座约束定位向量法,考虑支座约束定位向量法,考虑支座约束刚度矩阵可逆刚度矩阵可逆
17、刚度矩阵可逆刚度矩阵可逆稀疏矩阵与带状矩阵:稀疏矩阵与带状矩阵:连续梁连续梁简支单元简支单元本讲稿第三十二页,共八十四页整体分析整体分析建立结构的刚度方程建立结构的刚度方程即结构即结构结点力结点位移结点力结点位移的刚度关系的刚度关系 F F K K 其中:其中:F F 结构结构结点荷载结点荷载 结构结构结点位移结点位移 K K 结构结构刚度矩阵刚度矩阵本讲稿第三十三页,共八十四页一般结构(图一般结构(图10106 6)FF结构结点荷载列向量结构结点荷载列向量R1R4本讲稿第三十四页,共八十四页结构结点位移列向量结构结点位移列向量u2u3v2v323几何条件几何条件变形协调:变形协调:杆端位移相
18、应的结点位移杆端位移相应的结点位移本讲稿第三十五页,共八十四页刚度方程刚度方程分块形式分块形式分块展开:分块展开:引入支座条件引入支座条件的结构刚度方的结构刚度方程:程:消除了刚体位消除了刚体位移,移,总刚总刚K非非奇异奇异,可解未,可解未知位移。知位移。可求支座反力可求支座反力(1032)(1033)本讲稿第三十六页,共八十四页位移法位移法典型方程典型方程可解可解位移位移各杆各杆杆端位移杆端位移e单元刚度方程单元刚度方程求求杆端力杆端力(整体整体)坐标转换坐标转换局部坐标局部坐标系的系的杆端力杆端力(1029、34)(1032)关键关键总刚总刚本讲稿第三十七页,共八十四页整体分析整体分析集成
19、总刚(集成总刚(10105 5支承条件引入)支承条件引入)定位向量法定位向量法单元定位向量单元定位向量e其中:其中:Ni i表示杆端表示杆端i i对应的对应的整体整体结点编码结点编码;n ni i 表示杆端位移表示杆端位移i i所对所对应的应的整体结点位移编码整体结点位移编码支座结点相应的支座结点相应的约束位移约束位移方向的编码为方向的编码为 0 0 。本讲稿第三十八页,共八十四页1、编码、编码结点号:结点号:0、1、2、0(支座结点编码为(支座结点编码为0 0)单元号:单元号:、坐标系:坐标系:整体坐标系整体坐标系XOYXOY局部坐标系局部坐标系xoyxoy1212、2323、3434TTe
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第十章矩阵位移法 2精选文档 第十 矩阵 位移 精选 文档
限制150内