矩阵位移法习题课件

广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE1例例1 1：已知某单元的单元坐标下的杆端力：已知某单元的单元坐标下的杆端力：作出单元的内力图。作出单元的内力图。1.241.52.61.241.53.4jiyx-1.241.52.63.4N图图Q图图M图图单位（单位（KN）单位（单位（KN.m）广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE23.044.381.240.43M(KN.m)Q(KN)N(KN)例例2 2：已知单元内力图，列出单元的杆端力向量：已知单元内力图，列出单元的杆端力向量xy1.240.434.381.240.433.04单元坐标表示：单元坐标表示：整体坐标表示：整体坐标表示：广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE32020701234(1)(2)(3)xy作业：已知单元和结点的离散如图，给定荷载作用作业：已知单元和结点的离散如图，给定荷载作用下各结点整体坐标下的位移：下各结点整体坐标下的位移：写出单元坐标及整体坐标表示的单元杆端写出单元坐标及整体坐标表示的单元杆端位移向量，整体坐标表示的单元杆端力向位移向量，整体坐标表示的单元杆端力向量，并作出单元的内力图。量，并作出单元的内力图。单元坐标下：各单元的杆端力向量如下：单元坐标下：各单元的杆端力向量如下：广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE4习题习题1 1：列出下列结构各单元的刚度矩阵。：列出下列结构各单元的刚度矩阵。习题广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE5习题n习题习题2 2：求图示各单元用单元坐标和整体坐标描述求图示各单元用单元坐标和整体坐标描述的单元刚度矩阵。已知：各单元坐标系如图所示，的单元刚度矩阵。已知：各单元坐标系如图所示，各杆各杆EIEI常量。如果忽略轴向变形，有何不同常量。如果忽略轴向变形，有何不同?yl lx2143广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE6作业作业1 1：按概念列出两端刚结点平面梁单元在单元：按概念列出两端刚结点平面梁单元在单元及整体坐标系下的单元刚度方程，画出计算单元杆及整体坐标系下的单元刚度方程，画出计算单元杆端力和变形的示意图。端力和变形的示意图。l作业作业作业2 2：写出课本：写出课本P429P429页题页题9-69-6、9-79-7中的各单元刚中的各单元刚度矩阵。度矩阵。7例例3.列出先处理法各单元整体坐标下的刚度列出先处理法各单元整体坐标下的刚度矩阵。已知各杆矩阵。已知各杆EI、EA均为常量。均为常量。1 2 3 41 2 3 42 3 4 234例例4.平面桁架如图所示，各杆截面平面桁架如图所示，各杆截面EA均均为常数。已知为常数。已知P1=15kN,P2=20kN，试桁架各杆轴力。试桁架各杆轴力。1.结构离散如图所示结构离散如图所示2.单元坐标描述的各单元刚度矩阵单元坐标描述的各单元刚度矩阵后处理法后处理法先处理法先处理法3.整体坐标描述的各单元刚度矩阵整体坐标描述的各单元刚度矩阵单元（单元（1 1）的单元坐标和整体坐标一致）的单元坐标和整体坐标一致后处理法后处理法先处理法先处理法单元（单元（2 2）需要坐标变换，）需要坐标变换，a=30=300 0先处理法先处理法后处理法后处理法4.装配结构总刚度矩阵，建立结构的刚装配结构总刚度矩阵，建立结构的刚度方程度方程先处理法先处理法u 总刚度矩阵：总刚度矩阵：u 荷载向量：荷载向量：u 结构刚度方程：结构刚度方程：u 位移向量：位移向量：5.求解结构的刚度方程可得：求解结构的刚度方程可得：6.单元坐标下各单元的杆端位移向量：单元坐标下各单元的杆端位移向量：整体坐标下：整体坐标下：单元坐标下：单元坐标下：7.单元坐标下各单元的杆端力向量：单元坐标下各单元的杆端力向量：单元内力：单元内力：例例 4.平面刚架如图所示，各杆截面相同。平面刚架如图所示，各杆截面相同。E=1107kN/m2,A=0.24m2,I=0.0072m4，求各杆端力，并画出内力图。求各杆端力，并画出内力图。解：解：1.结构离散如图所示结构离散如图所示;单元单元 a=0单元单元a=45 2.列出单元参数表列出单元参数表;3.单元坐标表示的单单元坐标表示的单元刚度矩阵元刚度矩阵 先处理法先处理法例例 4.平面刚架如图所示，各杆截面相同。平面刚架如图所示，各杆截面相同。E=1107kN/m2,A=0.24m2,I=0.0072m4，求各杆端力，并画出内力图。求各杆端力，并画出内力图。解：解：1.结构离散如图所示结构离散如图所示;单元单元 a=0单元单元a=45 2.列出单元参数表列出单元参数表;3.单元坐标表示的单单元坐标表示的单元刚度矩阵元刚度矩阵 先处理法先处理法3.单元坐标表示的单元刚度矩阵单元坐标表示的单元刚度矩阵 将各单元的参数代入梁单元的统一刚度矩阵公式：后处理法后处理法4.整体坐标表示的单元刚度矩阵整体坐标表示的单元刚度矩阵 单元（1）(3)的单元坐标和整体坐标一致，所以先处理法先处理法单元（2）的单元坐标和整体坐标成450角，要坐标变换4.整体坐标表示的单元刚度矩阵整体坐标表示的单元刚度矩阵 后处理法后处理法广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE19结构的结点位移和结点力向量结构的结点位移和结点力向量v例例5。写出图示结构的结点位移向量写出图示结构的结点位移向量和结点力向量和结点力向量P。先处理法：先处理法：后处理法：后处理法：支座反力支座反力支座反力支座反力自由结点外力自由结点外力20v例例6。写出图示结构的结点位移向量和结点力向量。写出图示结构的结点位移向量和结点力向量。先处理法先处理法后处理法后处理法如果忽略轴向变形，则：如果忽略轴向变形，则：21v写出图示结构矩阵位移法的基本写出图示结构矩阵位移法的基本未知量。未知量。先处理法先处理法后处理法后处理法v例例7。写出图示结构的结点位移向量。写出图示结构的结点位移向量。先处理法：先处理法：后处理法：后处理法：17*3=51考虑轴向变形：考虑轴向变形：14*3=4214*3=42忽略轴向变形：忽略轴向变形：14+5=1914+5=1922例例8.列出先处理法的总刚度方程。列出先处理法的总刚度方程。EI、EA均为常量。均为常量。2 3 4 234贡献矩阵贡献矩阵整体坐标描述的单元刚度矩阵整体坐标描述的单元刚度矩阵1 2 3 41 2 3 4定位向量定位向量贡献矩阵直接代数相加贡献矩阵直接代数相加23常用作法，不写出贡献矩阵，根据定位向量直接装配总刚。常用作法，不写出贡献矩阵，根据定位向量直接装配总刚。1 2 3 41 2 3 42 3 4 2341 2 3 41 2 3 4装配总刚矩阵：装配总刚矩阵：荷载向量：荷载向量：结构的刚度方程结构的刚度方程:24例例9 9：写出图示连续梁的整体刚度矩阵。：写出图示连续梁的整体刚度矩阵。xy(0,1,0)(0,0,2)(0,0,3)(0,4,0)1 2 12定位向量定位向量2 3 233 4 341 2 3 41 2 3 4装配总刚矩阵：装配总刚矩阵：单元刚度矩阵：单元刚度矩阵：如：图示平面刚架的总刚度矩阵的集成如：图示平面刚架的总刚度矩阵的集成n各结点编号如图。用矩阵可记为：n如果求出各单元整如果求出各单元整体坐标下的单元刚体坐标下的单元刚度矩阵：度矩阵：1 2 3 4 5 6123456n根据结点编号把各根据结点编号把各单元的子块搬入总单元的子块搬入总刚刚K中的对应位置。中的对应位置。同一位置各子块的同一位置各子块的对应元素相加；空对应元素相加；空位补位补0。n如果修改各结点编如果修改各结点编号如图。则用矩阵：号如图。则用矩阵：1 2 3 4 5 6123456n把各单元的子块搬把各单元的子块搬入总刚度矩阵入总刚度矩阵K中中的相应位置。的相应位置。v单元的子块搬入单元的子块搬入总刚度矩阵中的总刚度矩阵中的位置，完全取决位置，完全取决于结构结点编号。于结构结点编号。对同一结构，如对同一结构，如果改变了结点的果改变了结点的编号，则总刚度编号，则总刚度矩阵完全不同。矩阵完全不同。u判断图示结构矩阵位移法（后处理法）总刚判断图示结构矩阵位移法（后处理法）总刚矩阵的最大半带宽。矩阵的最大半带宽。d=(3+1)*3=1227d=(6+1)*3=21n总刚矩阵中元素的排列与结点的顺序直接相关。总刚矩阵中元素的排列与结点的顺序直接相关。d=(7+1)*2=1628例例10.列出后处理法的总刚度方程。列出后处理法的总刚度方程。EI、EA均为常量。均为常量。1 2 122.2.整体坐标描述的单元刚度矩阵整体坐标描述的单元刚度矩阵2 32 31.结构离散如图结构离散如图123 1 2 33.3.装配总刚矩阵：装配总刚矩阵：1 2 31 2 12 1 2 3装配总刚矩阵：装配总刚矩阵：1 2 32 3 2331结构刚度方程：结构刚度方程：荷载向量：荷载向量：位移向量：位移向量：主元素置主元素置1 1 法修正后的结构刚度方程：法修正后的结构刚度方程：1234等效移置的方法等效移置的方法:首先求出首先求出基本结构基本结构在非在非结点荷载作用下引起的结点荷载作用下引起的固端力固端力。最后将。最后将各固端力反向作用到结构的结点上去，各固端力反向作用到结构的结点上去，即为该结点的等效结点荷载。即为该结点的等效结点荷载。1234等效结点荷载向量等效结点荷载向量:荷载向量荷载向量:2312荷载等效图荷载等效图基本结构基本结构:两端固定的梁。两端固定的梁。例例1111：计算图示结构的荷载向量。：计算图示结构的荷载向量。荷载等效图荷载等效图 原始荷载向量：原始荷载向量：2.先处理法先处理法1.后处理法后处理法考虑边界后荷载向量：考虑边界后荷载向量：考虑轴向变形：考虑轴向变形：忽略轴向变形：忽略轴向变形：例例1212：计算图示结构的荷载向量。：计算图示结构的荷载向量。荷载等效图荷载等效图 原始荷载向量：原始荷载向量：2.先处理法先处理法1.后处理法后处理法考虑边界后荷载向量：考虑边界后荷载向量：考虑轴向变形：考虑轴向变形：忽略轴向变形：忽略轴向变形：例例1313：写出图示结构的荷载向量。写出图示结构的荷载向量。原始荷载向量：原始荷载向量：1.后处理法后处理法考虑边界后荷载向量不变：考虑边界后荷载向量不变：2.先处理法先处理法例例14：写出图示结构矩阵位移法的刚度方程。写出图示结构矩阵位移法的刚度方程。返回目录返回目录111 212232 3（0,0,1）（0,0,2）（0,0,3）1231 2 3结构的刚度方程结构的刚度方程:例例1515：计算图示结构的内力。：计算图示结构的内力。忽略轴向变形，已知求得结忽略轴向变形，已知求得结构的结点位移向量：构的结点位移向量：2.单元坐标下各单元的位移向量单元坐标下各单元的位移向量1.单元坐标下各单元的刚度矩阵单元坐标下各单元的刚度矩阵为了同时计算弯矩和剪力，取四项：为了同时计算弯矩和剪力，取四项：单单元元信信息息例例1515：计算图示结构的内力。：计算图示结构的内力。忽略轴向变形，已知求得结忽略轴向变形，已知求得结构的结点位移向量：构的结点位移向量：4.单元杆端力向量单元杆端力向量3.单元固端力向量单元固端力向量394041弯矩图弯矩图剪力图剪力图5.根据杆端力向量作内力图根据杆端力向量作内力图广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE1.1.计算图示结构计算图示结构先处理法先处理法的刚度方程。的刚度方程。习题训练广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE43习题训练2.2.求图示各单元用单元坐标和整体坐标描述的单元求图示各单元用单元坐标和整体坐标描述的单元刚度矩阵，形成结构的总刚度矩阵。已知：各单刚度矩阵，形成结构的总刚度矩阵。已知：各单元坐标系如图所示，各杆元坐标系如图所示，各杆EIEI、EAEA均为常量。均为常量。ll l1yx243广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE3.写出图示结构矩阵位移法结点荷载列阵写出图示结构矩阵位移法结点荷载列阵P。（包括。（包括后处理和先处理法）后处理和先处理法）习题训练yx广西大学土木建筑工程学院广西大学土木建筑工程学院COLLEGEOFCIVILENGINEERINGANDARCHITECTURECOLLEGEOFCIVILENGINEERINGANDARCHITECTURE4.求图示结构各元的杆端力，并画出弯矩图。已知：求图示结构各元的杆端力，并画出弯矩图。已知：结点结点1的位移：的位移：习题训练