有限元分析及应用案例分析.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：平面三结点三角形单元的结点位移和坐标变换 由于该单元的结点位移是以整体坐标系中的x方向位移(u1)和y方向位移(v1)来定义的，所以没有坐标变换问题。 讨论2：平面三结点三角形单元的应变矩阵和应力矩阵为常系数矩阵 单元的位移场为线性关系，由于只与(xi，yi)相关，是常系数，因而求出的 和 为常系数矩阵，不随x、y变化，即这种单元在单元内任意一点的应变和应力都相同。因此，三结点三角形单元称为常应变单元，在应变梯度较大（即应力梯度比较大）的部位，单元划分应适当密集，否则将不能反映应变的真实变化而导致较大的误差。 由于有四个位移结点条件，可假设纯弯梁单元的位移场为具有四个待定系数的函数模式，即 (1) 单元位移场的表达 其中 为待定系数。 由该单元的结点位移条件 可求出 中的4个待定系数，即 将上式代入 中，重写位移函数，有 其中， ， 叫做单元的形状函数矩阵，即 (2) 单元应变场的表达 由纯弯梁的几何方程，有梁的应变表达式 其中 为基于中性层的坐标， 叫做单元的几何函数矩阵，即 其中 (3) 单元应力场的表达 其中 弹性模量， 叫做单元的应力函数矩阵 该单元的势能为 (4) 单元势能 的表达 其中应变能 其中 为刚度矩阵，即 其中 为惯性矩，则外力功为 (5) 单元的刚度方程 同样，由最小势能原理，将 对 取一阶极小值，有刚度方程 其中刚度矩阵 和力矩阵 分别在以上的计算中给出。注意上式中的下表(4×4)(4×1)(4×1)为各个矩阵的维数（即行和列）。