杆梁结构有限元(第四章)描述.pptVIP

4.3 梁单元有限元分析的标准化标准与算例 简支梁的最小势能原理解 仍以如图4-9所示的平面简支梁的弯曲问题为例，为提高计算精度，可以选取多项函数的组合，这里取满足位移边界条件BC(u)的许可位移场： 计算应变能U为 4.3 梁单元有限元分析的标准化标准与算例 相应的外力功W为 则总势能为Π=U-W，为使取极小值，则有 则 解出 4.3 梁单元有限元分析的标准化标准与算例 可以看出，该方法得到的第一项与前面虚功原理求解出来的结果相同，与精确解相比，该结果比前面由虚功原理得到的结果更为精确，这时因为选取两项函数作为试函数，这也是提高计算精度的重要途径。以上求解过程所用的试函数为许可基底函数的线性组合，因此，上述求解方法也是瑞利-里兹方法。 以上的求解，都是基于试函数的能量方法(泛函极值法)，基本要点是不需求解原微分方程，但需要假设一个满足位移边界条件BC(u)的许可位移场。因此，如何寻找或构建满足所需要求的许可位移场是一个关键，并且，还期望这种构建许可位移场的方法还应具有标准化和规范性。下面的重点将讨论通过基于“单元”的位移函数的构建就可以满足这些要求。 4.3 梁单元有限元分析的标准化标准与算例 局部坐标系中的平面梁单元 图4-12所示为一局部坐标系中的纯弯梁单元，其长度为l，弹性模量为E，横截面的惯性矩为Iz。 单元的几何及节点描述 设有两个端节点，节点位移列阵： 节点力列阵： 其 中分别为各节点的挠度和转角 单元位移场的表达 由于有4个位移节点条件，可假设纯弯梁单元的位移场挠度为具有四个待定系数的函数模式，即 ： 4.3 梁单元有限元分析的标准化标准与算例 由该单元的节点位移条件 可求出式4个待定系数，即 重写位移函数，有 其中 叫做单元的形状函数矩阵，即 单元应变场的表达 由纯弯梁的几何方程，有梁的应变表达式 4.3 梁单元有限元分析的标准化标准与算例 其中 是以中性层为起点的y方向的坐标， 叫做单元的几何矩阵： 其中 单元应力场的表达 由梁的物理方程 其中E为弹性模量，S(x)叫做单元的应力矩阵。 4.3 梁单元有限元分析的标准化标准与算例 单元势能的表达 该单元的势能为 其中应变能 刚度矩阵，具体地有 4.3 梁单元有限元分析的标准化标准与算例 外力功为 其中 单元的刚度方程 4.3 梁单元有限元分析的标准化标准与算例 一般平面梁单元的描述 为推导局部坐标系中的一般平面梁单元，在图4-13所示的纯弯梁的基础上叠加进轴向位移(由于为线弹性问题，满足叠加原理)，这时的节点位移自由度(DOF)共有6个，见图4-13。 图4-13所示平面梁单元的节点位移列阵 和节点力列阵 ： 相应的刚度方程： 将杆单元刚度矩阵与纯弯梁单元刚度矩阵进行组合，可得到刚度矩阵： 4.3 梁单元有限元分析的标准化标准与算例 受均布载荷平面梁单元的等效节点载荷 (a) 几种受均布载荷作用的梁构件 (b) 节点等效载荷 4.3 梁单元有限元分析的标准化标准与算例 图4-15 受均布载荷作用的直接静力等效的 节点载荷(每个节点分一半) 对于(a)所示的几种受均布载荷的情况，若需要采用几个或一个梁单元，均可以建立如图(b)所示的单元。 其节点位移列阵为 节点力列阵为 单元的挠度位移场 其中 计算该单元的外力功为W为 4.3 梁单元有限元分析的标准化标准与算例 将形状函数矩阵代入上式，可计算出节点力 讨论1：若凭一种直觉，直接按照静力等效的方式来进行计算，即，每个节点各分一半进行静力等效，见图4-15，则计算出的节点等效力为 将该结果与上面进行比较，显然这样计算出的M1和M2都是错误的。 讨论2：该等效节点载荷是按照外力功进行计算的，是通用的均布载荷的节点等效载荷，与节点的实际约束状态没有关系。也就是说，图(a)中的几种情况的节点等效载荷都用式（本页第一个公式)。 4.3 梁单元有限元分析的标准化标准与算例 平面梁单元的坐标变换 图4-17所示为一整体坐标系中的平面梁单元，它有两个端节点，梁的长度为l，弹性模量为E，横截面的面积为A，惯性矩为Iz。 设局部坐标系(Oxy)下的节点位移列阵为 设整体坐标系(Oxy)下的节点位移列阵为 变换关系 写成矩阵形式有 4.3 梁单元有限元分析的标准化标准与算例 坐标变换矩阵 梁单元在整体坐标系中的刚度方程为 其中 4.3 梁单元有限元分析的标准化标准与算例 空间梁单元及坐标变换