7轴对称旋转单元.pptVIP

南京航空航天大学能源与动力学院机械振动冲击仿真研究室(PC:210016) Tel:(025)4892220 Fax:(025)4895966 E-mail:nuaacw@263.net 7.1 应力与应变关系 第七章 轴对称旋转单元 7.2 单元刚度矩阵 7.3 ANSYS轴对称旋转单元计算示例 轴对称结构在工程中应用比较广泛。轴对称结构是由任意平面图形绕着某直线旋转一周而形成的回转体，该直线称为对称轴；旋转平面称为子午面。如果轴对称结构的约束条件以及作用的载荷都对称于对称轴的话，则在载荷作用下产生的位移、应变和应力也对称于此轴，这种问题称为轴对称问题。 图7-1 轴对称旋转体 第七章 轴对称旋转单元 第七章 轴对称旋转单元 当几何形状、约束条件以及作用的载荷全部沿着对称轴周向保持不变时 ，结构产生轴对称变形，且其位移、应变、应力都与夹角θ无关，而仅是径向坐标r和轴向坐标x的函数。也就是说，在任何一个过 x 轴的子午面上的位移、应变和应力的分布规律都相同。这类轴对称问题可称之为完全轴对称问题，结构分析时可简化为平面问题。否则按空间问题处理。对完全轴对称问题进行计算时，只需取任一个子午面离散化后，进行分析即可。 7.1 应力与应变关系 图7-2 轴对称问题中应变与应力状态 第七章 轴对称旋转单元 右图给出了轴对称旋转体问题中的应力和应变分量。令旋转轴为x、径向轴为r、环向坐标为θ。沿轴 x和轴 r的位移分量为u和v，它们都是坐标r、x的函数。在轴对称情况下，任一径向位移都会引起周向应变，所以必须考虑周向应变和应力分量。 第七章 轴对称旋转单元 式中 ——弹性矩阵； ——初应变； ——初应力。 应变位移关系为 应力应变关系矩阵表达式为 对于各向同性材料有 E ——材料弹性模量；μ ——料泊松比。 由于温度变化而引起的初 应变，对于各向同性材料有 T ——单元的平均温升； ——热膨胀系数。 第七章 轴对称旋转单元 如上所述完全轴对称问题的分析，可以转化为对其任意子午面的分析，可将此截面离散为许多三结点三角形单元。与平面问题相似，每个单元的位移向量记为 假定单元内位移为坐标的线性函数，与平面问题三角形单元类似，位移函数可表示为 用形函数表达的单元内任意点位移 第七章 轴对称旋转单元 7.2 单元刚度矩阵 或写为 图7-3 四边形旋转体单元 四边形旋转体单元，设旋转轴为x轴，经向轴为y轴，结点按逆时针方向排序。 单元坐标变换为 ——形函数 单元位移列向量用单元结点位移插值来表示 ——插值函数，它与形函数相同。 第七章 轴对称旋转单元 写成矩阵形式 i=1， 2， 3， 4 第七章 轴对称旋转单元 单元刚度矩阵的积分是在环体上进行的，因此有 采用高斯数值积分方法，得到四边形旋转体单元的刚度矩阵 第七章 轴对称旋转单元 圆筒直径0.4m，高度0.6m，壁厚0.005m；材料Q235，弹性模量E=2.1e11Pa，泊松比μ=0.3；约束：圆筒的下部在轴线方向固定，其它方向自由；载荷：顶部环线上承受轴向线压力P＝-200000N/m。 图7-4 圆筒示意图 图7-5单元类型对话框 第七章 轴对称旋转单元 7.3.1 问题描述 7.3 ANSYS轴对称旋转单元计算示例 第七章 轴对称旋转单元 （1）选择单元类型 运行PreprocessorElement TypeAdd/Edit/Delete，弹出Element Types对话框单击Add，弹出Library of Element Types对话框，如图7-6所示，选择SHELL51单元。 7.3.2 ANSYS求解操作过程 图7-6 单元类型库对话框 图7-7 选择材料属性对话框 （2）设置材料属性 运行PreprocessorMaterial PropsMaterial Models，弹出Define Material Model Behavior对话框，如图7-7所示。双击Isotropic选项，弹出Linear Isotropic Properties for Material Number1对话框，如图7-8所示。 图7-8 设置材料属性对话框 第七章 轴对称旋转单元 （3）定义单元实常数 选择Main MenuPreprocessor Real Constants Ad