2、空间问题有限元法.pptVIP

空间轴对称问题，采用圆柱坐标系 ，根据轴对称问题的特点，任一对称面为 面，没有 方向的位移。物体上任一点的位移、应变、应力均与 角无关，轴对称问题可以把三维问题转化为以 为自变量的二维问题。 图3-4代表三角形环形单元，用在 平面上的三角形 表示。轴向位移为 ，沿半径方向的位移为 。 则，三角形单元的节点位移为： 与平面问题三角形单元相仿，取位移模式： 带入节点位移后，可解出 ，用形函数表达： 矩阵中含有变量r,z,因此，不是常数矩阵，三角形环形单元也不是常应变单元。 根据弹性力学理论，空间轴对称问题的弹性方程： 其中，弹性矩阵为： 轴对称问题的单元刚度矩阵 三角形环形单元的积分可简化为三角形截面上的积分： 其体积力，面力的表达式： 车辆工程技教研室 机电工程学院 * 第三章 空间问题有限单元法 实际工程中，对于那些形体复杂，三个方向尺寸同量级的结构，必须按空间（三维）问题求解。 空间问题的有限单元法中的位移仍然只有平动位移，所以仍属于C0连续问题，因此构造单元并不难。将平面问题有限元法“稍加变动”并“加以推广”便可用于空间问题。 由平面问题转为空间问题，给有限元分析带来两个主要难题： 1、空间离散化不太直观，人工离散很容易出错。 2、未知量的数量剧增，对计算机的存储和计算时间要求较高。 解决问题： 1、编程建模 2、采用高精度单元 由于通用软件有很好的前后处理功能，因此空间问题基本上都靠软件来解决。 一、空间问题常用单元 二、常应变四面体单元 一、空间问题常用单元 2. 按位移函数阶次分 1. 按形状分： 四面体单元（三棱锥） 五面体单元（三棱柱） 六面体单元（立方体） 线性单元：四结点四面体，六结点五面体、八结点六面体等 二阶单元：十结点四面体，二十结点六面体等 三阶单元：二十结点四面体，三十二结点六面体等 一、空间问题常用单元 4. 构造曲面单元 3. 形函数构造方法： 1）广义坐标法：仅用在常应变单元 等参元：利用规则单元作母元，通过等参变换构造曲面单元 2）试凑法：在自然坐标下直接写出形函数 四面体单元的自然坐标是体积坐标 二、常应变四面体单元 1. 基本变量 单元内任一点位移: 单元内任一点应变: 单元内任一点应力: 二、常应变四面体单元 1. 基本变量 单元结点位移: 结点位移: 2. 单元位移插值函数： 设单元内任一点的位移为坐标的线性函数: 即为广义坐标 二、常应变四面体单元 将结点坐标代入u(x,y,z)，得结点x方向位移: （四个方程、四个未知量） 2. 单元位移插值函数： 二、常应变四面体单元 解方程组得 后，可将u的表达整理成： 式中: 按1、2、3、4的顺序变换下标，可得其它系数 令: 则: 同样的过程可得到: 形函数 则单元位移模式可写成： （由结点位移表示的单元内位移） 或： 形函数矩阵 3. 单元几何方程： 二、常应变四面体单元 由结点位移求单元内应变 将位移表达式代入，得： 其中： [B]中各元素为常数，则{?}也为常量 — 常应变单元 4. 单元物理方程： 二、常应变四面体单元 由结点位移求单元内应力 — 应力矩阵 令: 其中： 其中： 5. 单元基本方程： 二、常应变四面体单元 利用变分原理建立单元平衡方程： 其中： 单元刚度矩阵 等效结点荷载 三、空间轴对称问题的有限元法 三、空间轴对称问题的有限元法 三、空间轴对称问题的有限元法 该式与平面问题的有限元问题中的形函数形式一样，系数相同，只是将平面问题的坐标轴x,y换成 其中， 车辆工程技教研室 机电工程学院