船舶结构力学能量法、矩阵法.pptVIP

解:(1)根据结构的受力特点，将它离散为3个单元，4个节点，并建立杆元的局部坐标及结构的总坐标如上图所示。 (2)计算杆元的刚度矩阵。 杆元①: 杆元①需进行坐标转换，因 ,故坐标转换矩阵为 则杆元①在总坐标系中的刚度矩阵为 杆元②与③的局部坐标与总坐标一致，注意到此二杆元的长度为L/2，故有 (3)根据各杆元刚度矩阵的分割子矩阵，组成结构刚度矩阵: (4)求节点外载荷矩阵从而写出节点平衡方程式。 杆元①因三角形分布荷重引起固端弯矩及固端剪力，在单元坐标系中，固端力矩阵为 所以有 由此可列出节点平衡方程式形式如下： (5)约束处理。本例中，因不计杆件的轴向变形又有 ,因此在上式中需划去与 对应的七个行与列。 在节点4处y方向有弹性支座，在总刚相应行的主对角线元素上加上弹性支座的刚度后可不计其外载荷项。故经约束处理后的方程式为 解此方程式组，即可得节点位移。 例6 用矩阵法建立下图中梁的总刚度矩阵，并对节点平衡方程式进行约束处理， 计算时节点号码自左向右取为1、2、3…，杆元号码自左向右取为① ②… 将各子块代入得： 结构力学习题 ——能量法、矩阵法 一、能量法 例1 结构体系如下图所示，试利用应力能原理计算结构右端点的挠度位移。不考虑剪力影响。 解：由叠加法原理，可将上述杆系拆分为： 由力的平衡条件可以知F=3P/2; 图（a）中 其应变能为 图（b）中 故杆系的总变形能 对于线性体系而言，由卡式第二定理知右端挠度为 例2 用李兹法求图6-b中变断面梁的中点挠度。已知 , 。计算时试取挠曲线函数 。 解：取挠曲线函数为 时, 满足梁两端的位移边界条件，即 x=0时 x=L时 说明此挠曲线函数满足李兹法的要求，下面进行计算。 (1)计算应变能。此梁的应变能包括两部分，一是梁本身的弯曲应变能 ，二是弹性支座的应变能 。注意到梁是变断面的，故有 将 及 带入可算出 总的应变能为 (2)计算力函数。此梁的力函数包括集中力P引起 及分布荷重 引起的两部分。 计算时，先写出分布荷重的表达式。对图示坐标，有 故 则总的力函数为 (3)计算总位能 由 得 由此解得 故梁的挠曲线方程为 梁中点的挠度为 例3 用最小功原理求解下图所示梁的静不定性。已知弹性支座的柔性系数 。 解：设弹性支座处的支反力为F；两端支 反力为 。 该双跨梁的弯矩表达式为 由于是线性结构，该杆系的余能等于应变能，即 ；且该杆系的应变能有两部分组成，杆本身的应变能 和弹性支座所具备的应变能 。 故总应变能为 由最小功定理 知 弹性支座处的挠度为 二、矩阵法 例4 用矩阵位移法求解该梁的剪力和弯矩。 y x ① ② 解：把该连续梁离散为3个节点1、2、3，2个单元（1）、（2）。如上图所示。 计算杆元刚度矩阵： 总刚度矩阵为 列平衡方程 由题意得初始条件 划去1、2、5行列得 : 解得 例5 用矩阵法计算下图中的平面刚架，写出结构刚度矩阵及经约束处理后的平衡方程式组。已知 ,计算时杆件的轴向变形不计。