有限元及工程软件第三章.ppt

* * 第三章 轴对称问题的有限元法 第一节 轴对称问题 几何形状，所受载荷对称于中心轴，则其变形 也对称于此轴，这种问题称为轴对称问题。 坐标系：柱坐标系，以中心轴为z轴 通过中心轴的平截面－－子午面。 位移、应力和应变与角坐标 无关，只是径向r和轴向 z的函数。 研究子午面的变形 二维问题 子午面内各点位移为 当圆周上各点都存在径向位移时， 圆周被拉伸，产生环向应变 几何方程 应力向量为 物理方程 写成矩阵形式 弹性矩阵 第二节 轴对称问题的简单三角形单元 3个节点的坐标分别为 三角形单元的节点位移向量为 三角形单元的节点力 单元刚度矩阵 注意： 单元为圆环状 、单元位移函数和形函数 选取如下位移函数 将3个节点的坐标代入写成矩阵表达形式 ――三角形面积 （i＝1，2，3）为 求逆过程中的代数余子式 如 其中 （i=1，2，3） 插值函数，是坐标的函数，反映单元的位移形态， 在有限元法中称为形函数。 单元的形函数矩阵 二、应变矩阵 单元位移函数为 几何方程 应变矩阵 应变矩阵是由节点位移求单元内任一点应变的转化矩阵 提问：轴对称问题的三角形单元是否为常应变单元？ 应力矩阵 三、应力矩阵 应力矩阵是由节点位移求单元内任一点应力的转化矩阵。 四、单元刚度矩阵 变形体的虚功原理：要使变形体在某一形变位置处 于平衡，其充要条件是，在这一变形位置，所有内 力和外力在任何虚位移上所做的虚功之和为零。 --变形体虚功方程 设单元产生虚位移，单元节点虚位移为 单元内部的虚应变为 外力虚功 内力虚功的负值 注意：三角形单元代表 环状单元体 子午面 内三角 形单元 的面积 域 任意性 其中的元素为常量 单元刚度矩阵 单元节点位移求节点力的转换矩阵 单元刚阵的近似积分 为了简化计算和消除在对称轴上r＝0所引起的 麻烦，把单元中随点而变化的r，z用单元截面形心 处的 和 来近似，即 应力矩阵和应变矩阵为常量阵 三角形环状 单元截面积 代入D，B 其中 单元刚度矩阵的精确积分 将应变矩阵分成两部分 其中： 用单元形心代入获得的应变矩阵不变部分 变化部分 由 常量阵 第二项和第三项积分为零 其中： 近似积分简单精度足够 一般采用近似计算 的积分包含下面的关系 五、 单元等效节点载荷向量 单元所受的非节点载荷一般包括：集中载荷、 分布体力、边界面力。 移置到节点上 依据虚功原理－－原载荷与等效节点载荷在任 意虚位移上的虚功相等。 1、集中载荷的等效节点载荷 作用集中载荷 等效节点载荷向量为 假设该单元产生虚位移 单元内各节点的虚位移 由 M点虚位移 根据虚功原理