矩阵位移法 结构力学 学习资料.ppt

《矩阵位移法 结构力学 学习资料.ppt》由会员上传分享，免费在线阅读，更多相关内容在PPT专区-天天文库。

1、基本要求：熟练掌握两种坐标系中的单元刚度矩阵、结构的整体刚度矩阵、等效结点荷载的形成，已知结点位移求单元杆端力的计算方法，整体刚度矩阵和结构结点荷载的集成过程。理解单元刚度矩阵中和整体刚度矩阵中的元素的物理意义。了解不计轴向变形时矩形刚架的整体分析.第十三章矩阵位移法MATRIXDISPLACEMENTMETHOD基本概念单元刚度矩阵整体刚度矩阵等效结点荷载计算步骤和算例忽略轴向变形时的整体分析上机作业§13-1概述矩阵位移法以传统的结构力学作为理论基础，以矩阵作为数学表达形式，以电子计算机作为计算手段，三位一

2、体的方法。手算与电算的不同：手算：怕繁，讨厌重复性的大量运算，追求机灵的计算技巧，运算次数较少的方法。电算：怕乱，讨厌头绪太多，零敲碎打的算法，追求计算过程程序化，通用性强的方法。矩阵位移法（有限单元法finiteelementmethod）的基本思路是：先将结构离散成有限个单元，然后再将这些单元按一定条件集合成整体。这样，就使一个复杂结构的计算问题转化为有限个简单单元的分析与集成问题。有限单元法的两个基本环节：1）单元分析：建立单元刚度方程，形成单元刚度矩阵(物理关系)2）整体分析：由单元刚度矩阵形成整体刚度

3、矩阵，建立结构的位移法基本方程(几何关系、平衡条件)单元刚度矩阵是用来表示杆端力与杆端位移之间的物理关系的，不是新东西，但有几点新考虑：重新规定正负规则，以矩阵的形式表示，讨论杆件单元的一般情况。杆端局部编码与局部坐标系eE，A，Il局部坐标系中的杆端位移分量1u2u1v2v2q1q2M2Y2X1X1M1Y局部坐标系中的杆端力分量{}=Feee={}De12yx§13-2单元刚度矩阵（局部坐标系）(elementstiffnessmatrix)单元刚度方程方程1q1u1v2v2q2u1X1Y1M2X2Y2M)(

4、,,1221uulNXNX-=D=-=llEAND=由虎克定律：由转角位移方程，并考虑：2QYBA=1,QYAB-=12,vv-=D()212212212)(6vvlEIlEIY--+-=qq()212212112)(6vvlEIlEIY-++=qq()212212642vvlEIlEIlEIM-++=qq()212211624vvlEIlEIlEIM-++=qqee=ee（13—5）单元刚度方程（13—6）单元刚度矩阵单元刚度矩阵的性质1）单元刚度矩阵是杆端力用杆端位移来表达的联系矩阵。2）其中每个元素称为单

5、元刚度系数，表示由于单位杆端位移引起的杆端力。如第个杆端位移分量=1时引起的第个杆端力2M1q三六ji?反力互等定理3）单元刚度矩阵是对称矩阵。4）第k列元素分别表示当第k个杆端位移=1时引起的六个杆端力分量。5）一般单元的单元刚度矩阵是奇异矩阵。不存在逆矩阵{}De{}Fe正问题力学模型将单元视为“两端有六个人工控制的附加约束的杆件”{}De控制附加约束加以指定。解的性质为任何值时，{}De{}Fe都有唯一的解答。且总是一个平衡力系，不可能是不平衡力系。{}De{}Fe反问题将单元视为“两端自由的杆件”。{}

6、Fe直接加在自由端作为指定的杆端力{}Fe为不平衡力系时没有静力解。{}De{}Fe为平衡力系时有无穷多组解。{}De1X1Y1M2X2Y2M1X1Y1M2X2Y2M2v=10312lEI312lEI-EI26l-EI26l-01v=10312lEI312lEI-EI26lEI26l0e第二列元素变符号即第五列第一列元素变符号即第四列第二行元素变符号即第五行第一行元素变符号即第四行特殊单元单元的某个或某些杆端位移的值已知为零。如梁单元、柱单元。特殊单元的单元刚度矩阵，可由一般单元的单元刚度矩阵删除与零杆端位移对

7、应的行和列得到。为了使计算过程程序化、标准化、自动化，只采用一般单元的刚度矩阵作为标准形式。各种特殊单元的刚度矩阵有计算机程序去自动形成。某些特殊单元的刚度矩阵是可逆的。121q2q122M1M选局部坐标系推导单元刚度矩阵方便且单元刚度矩阵的形式简单。选整体坐标系是为进行整体分析。按一个统一的坐标系来建立各单元的刚度矩阵单元坐标转换矩阵1X1Y1M2X2Y2Myx1X1Y1M2X2Y2Myxαyxyxα§13-3单元刚度矩阵（整体坐标系）局部坐标系中的杆端力整体坐标系中的杆端力单元坐标转换矩阵[T]是一正交矩阵

8、。同理：整体坐标系中的单元刚度矩阵设：将（a）、（b）代入（a）（b）1）表示在整体坐标系第j个杆端位移分量=1时引起的第i个杆端力分量。2）[k]@是对称矩阵。3）一般单元的[k]@是奇异矩阵。例13-1求图示刚架中各单元在整体标系中的单元刚度矩阵。设各杆的几何尺寸相同。l=5m，A=0.5m2,I=1/24m4E=3×107kN/m221与比较[I][k]@，[k]@同阶,性质类似