有限元法与程序-杆系结构单元1研究.ppt

* 第5章 杆系结构的有限元法 §5.1 概述 杆系结构：梁、拱、框架、桁架等，它们常可离散成杆元和梁元。 ○ ○ ○ ○ ○ ○ ○ ○ ○ 梁 拱 框架 ○ ○ ○ ○ ○ 桁架 坐标系 有限元中的坐标系有整体坐标系和局部坐标系。对于一个结构，整体坐标系一般只有一个；而局部坐标系有很多个，一个单元就有一个局部坐标。并且在局部坐标系中每一个单元的规定都是相同的，这样，同类型单元刚度矩阵相同。 X Y ○ ○ ○ ○ ○ P x y x y 杆系结构单元主要有桁架单元和梁单元两种类型。它们都只有2个节点i、j。 约定：单元坐标系的原点置于节点i；节点i到j的杆轴（形心轴）方向为单元坐标系中x轴的正向。 y轴、z轴都与x轴垂直，并符合右手螺旋法则。 对于梁单元， y轴和z轴分别为横截面上的两个惯性主轴。 x y z i j · · §5.2 杆单元 下图示出了一维杆单元，横截面积为A，长度为l，弹性模量为E，轴向分布载荷为px。单元有2个结点i，j，单元坐标为一维坐标轴x。 · · i j x l LINK px uj ui 1、一维杆单元 单元结点力向量： （1）位移模式和形函数 ① 位移模式 单元结点位移向量 因为只有2个结点，每个结点位移只有1个自由度，因此单元的位移模式可设为： 式中a1、a2为待定常数，可由结点位移条件 x=xi 时， u=ui x=xj 时， u=uj 确定。再将由此确定的a1、a2 其代入上式，得 ② 形函数 将位移模式改写为下列形式 （2）应变矩阵 一维杆单元仅有轴向应变 上式也可写为 式中[B]为应变矩阵 代入位移模式得 由应力应变关系 （3）应力矩阵 将 代入上式，得 式中[S]为应力矩阵 (4) 单元刚度矩阵 单元刚度矩阵由下式推出 对于等截面桁架杆单元（截面积为A ） ，v=Adx，故有： 由 得 (5) 等效节点力 单元上作用分布力px，则等效节点力计算公式仍为以下形式 当分布力集度px为常数时，有 2、平面杆单元（2D LINK1） ?1 ?2 ?3 ?4 i j x y l （1）单元坐标单元位移向量 ?1 ?2 ?3 ?4 i j x y 看成局部坐标下的拉压杆 应力矩阵 应变矩阵 单元刚度矩阵 形函数 等效节点力 i j x y l z 3、空间杆单元（3D LINK8） （1）单元坐标及单元位移向量 ?1 ?2 ?4 ?5 ?3 ?6 （2）形函数 （3）应变矩阵 （4）应力矩阵 (5) 等效节点力 (6) 单元坐标及单元刚度矩阵 对于等截面铰接杆单元， §5.3 梁单元 1、两端承受剪力、弯矩的平面梁单元 i j x y i j x y ?1 ?2 ?3 ?4 l F1 F2 F3 F4 l （1）局部坐标下单元位移和单元力 ① 单元位移 其中， v——y方向位移，即挠度。 ?——角位移。 ② 单元力 其中， Q——剪力 M——弯矩 （2）位移函数和形函数 ① 位移模式 设单元坐标位移模式为 i j x y ?1 ?2 ?3 ?4 l 梁单元内一点有2个位移： v、? 因为， ?=dv/dx；仅一个位移是独立的，取 v 。 式中 ② 形函数 由单元两端点的节点位移条件，解出式中的a1、a2、a3、a4。再代入该式，可将位移模式写为以下形式： 1 x y 1 y ? ? ? （3）应变矩阵 ① 单元弯曲应变?b与节点位移???e的关系。 梁单元上任一点的应变和该点挠度之间关系为： 将式 代入 ，得单元弯曲应变和单元位移之间关系 （4）应力矩阵 [S]=[D][B] (5) 等效节点力 对于梁上作用的集中力或集中力矩，在划分单元时可将其作用点取为结点，按结构的节点载荷处理。 这里仅考虑把单元上的横向分布载荷转化为等价节点力问题。 x y i j l py(x) 将形函数矩阵[N]代入上式，积分可得分布荷载的等效结点力。表5-1给出了几种特殊情况的等价节点力。 - 5ql2/96 ql/4 5ql2/96 ql/4 - ql2/20 7ql/20 ql2/3