有限元分析及应用精编版.ppt

(1) 原连续体（几何上）的逼近离散 对应于连续体的力学分析，有限元分析的一般过程如下： 其中 为单元。 (2) 单元特性的研究 研究单元特性以形成单元刚度矩阵和结点外载矩阵 结点自由度（位移）描述： 位移模式（简单性、完备性、连续性、唯一确定性） 由结点条件确定位移模式中的待定系数，推导出形状函数矩阵： 单元应变场的表达（由几何方程）： ：形状函数矩阵 ：弹性力学中几何方程算子 ：几何方程 单元应力场的表达（由物理方程）： ：弹性力学中的弹性系数矩阵 ：应力矩阵 单元势能的表达 在以上公式中， ：单元刚度矩阵 ：单元结点力矩阵 ：体积力向量 ：面积力向量 其中 几何的集成 对单元势能，应用最小势能原理，可得到单元的平衡关系 (3) 离散单元的装配和集成 结点位移的集成 刚度矩阵的集成 结点外载的集成 形成整体刚度方程 (4) 处理边界条件并且求解结点位移 (5) 求各单元内的应变、应力、支反力 三结点三角形2D单元如下图所示。三个结点为1、2、3，各自的位置坐标为(xi，yi)，i=1，2，3，各自的结点位移（分别沿x方向和y方向）为(ui，vi)，i=1，2，3。 5.2 2D单元（三结点、四结点）的构造 5.2.1 三结点三角形2D单元 上图所示三结点三角形2D单元，结点位移向量 和结点力向量 为 下面，我们需要将所有力学参量用结点位移向量 来表达。 (1) 单元位移场的表达 就三结点三角形2D单元，考虑到简单性、完备性、连续性及待定系数的唯一确定性原则，选取位移模式为 由结点条件，在x=xi，y=yi处，有 （1） （2） 将（1）代入结点条件（2）中，可求解（1）中的待定系数，即 在上述各式中， 上式（1，2，3）表示下标轮换，如1 2，2 3， 3 1。 将各系数代入（1）中，重写位移函数，并以结点位移的形式表示，有 写成矩阵形式， 其中 为形状函数矩阵，即 由弹性力学平面问题的几何方程 (2) 单元应变场的表达 其中几何函数矩阵 为 将形函数代入上式，有 其中 由弹性力学平面问题的物理方程 (3) 单元应力场的表达 其中平面应力问题的弹性系数矩阵为 将几何方程代入物理方程，有 为单元应力矩阵。 其中 t为平面问题的厚度。 (4) 单元的势能的表达 其中 是单元刚度矩阵，即 势能公式中的 为单元结点等效载荷，即 其中 为单元上作用有外载荷的边。 为线积分 (5) 单元的刚度方程 由于有四个位移结点条件，可假设纯弯梁单元的位移场为具有四个待定系数的函数模式，即 (1) 单元位移场的表达 其中 为待定系数。 由该单元的结点位移条件 可求出 中的4个待定系数，即 将上式代入 中，重写位移函数，有 其中， ， 叫做单元的形状函数矩阵，即 (2) 单元应变场的表达 由纯弯梁的几何方程，有梁的应变表达式 其中 为基于中性层的坐标， 叫做单元的几何函数矩阵，即 其中 (3) 单元应力场的表达 其中 弹性模量， 叫做单元的应力函数矩阵 该单元的势能为 (4) 单元势能 的表达 其中应变能 其中 为刚度矩阵，即 其中 为惯性矩，则外力功为 (5) 单元的刚度方程 同样，由最小势能原理，将 对 取一阶极小值，有刚度方程 其中刚度矩阵 和力矩阵 分别在以上的计算中给出。注意上式中的下表(4×4)(4×1)(4×1)为各个矩阵的维数（即行和列）。 4.4.2 局部坐标系中的平面梁单元 为推导平面问题中的梁单元的坐标变换公式，我们在纯弯梁的基础上叠加轴向位移（由于为线弹性问题，满足叠加原理），如下图所示 上图所示平面梁单元的结点位移 和结点力 为 相应的刚度方程为 将杆单元刚度矩阵与纯弯梁单元刚度矩阵进行组合，可得到刚度矩阵 4.4.3 平面问题中梁单元的坐标变换 局部坐标系下的结点位移 整体坐标系中的结点位移 注意：转角 和 在两个坐标系中是相同的。 按照两个坐标系中的位移向量相等效的原则，可推导出以下变换关系。 写成矩阵形式有 其中T为坐标变换矩阵，即 则整体坐标系中的