有限元分析和应用案例分析.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 为各个矩阵的维数（即行和列）。 4.4.2 局部坐标系中的平面梁单元 为推导平面问题中的梁单元的坐标变换公式，我们在纯弯梁的基础上叠加轴向位移（由于为线弹性问题，满足叠加原理），如下图所示 上图所示平面梁单元的结点位移 和结点力 为 相应的刚度方程为 将杆单元刚度矩阵与纯弯梁单元刚度矩阵进行组合，可得到刚度矩阵 4.4.3 平面问题中梁单元的坐标变换 局部坐标系下的结点位移 整体坐标系中的结点位移 注意：转角 和 在两个坐标系中是相同的。 按照两个坐标系中的位移向量相等效的原则，可推导出以下变换关系。 写成矩阵形式有 其中T为坐标变换矩阵，即 则整体坐标系中的刚度方程为 其中 空间