有限元分析论文

1机械1003班孙祥和3100301144基于高速旋转齿轮的有限元分析引言：齿轮泵是工程中较为常见的一种泵，在高速运转时齿轮受到多种力的作用，包括齿面受到的压力，啮合时的接触应力以及自身离心力。在此过程中，齿轮将发生形变，为此我们需要对其进行分析，确保其结构的稳定性，这对于齿轮泵安全有效地运行具有很重要的意义。关键词：高速齿轮、平面静力分析、接触应力分析、离心力分析一、分析对象这里我们分析的对象是齿轮泵中高速运转的齿轮，在ANSYS中我们建立了标准齿轮模型，其各项数据如下表所示齿顶直径24mm齿底直径20mm齿数10厚度4mm弹性模量2.06E11pa密度7.8e3kg/m3最大转速62.8rad/s摩擦系数0.1啮合齿轮中心距44mm表1齿轮泵高速齿轮参数二、平面静力分析1、分析问题为了考查齿轮泵在高速运转时，齿轮发生多大的径向位移，从而判断其变形情况，以及齿轮运转过程齿面受到的压力作用。在这里我们将齿轮的空间结构简化为平面模型，并分析其平面应力情况。此处的静力分析为线性静力分析，求解步骤分为建模、施加载荷和边界条件并求解、结果分析和评价三个步骤，下面依序进行。2、建立模型2.1定义单元类型根据齿轮的平面几何对称性和此处分析类型，我们选择四节点矩形单元PLANE42。PLANE42不仅可以用于计算平面应力问题，还可以用于分析平面应变和轴对称问题。每个节点2个自由度：x,y方向。具有塑性，徐变，膨胀，应力强化，大变形，大应变能力。2设定好单元类型后，对选择的PLANE42单元进行设置，在Elementbehavior（单元行为方式）选择Planestressw/wk。2.2定义实常数本处选用带有厚度的平面应力行为方式的PLANE42单元，需要设置器厚度实常数，只需在“Type1PLANE42”中将厚度设为4即可。2.3定义材料属性考虑惯性力的静力分析中必须定义材料的弹性模量和密度。2.4建立齿轮面模型，如下图所示图2建立齿轮面模型2.5对盘面划分网格选择MainMenu：PreprocessorMeshingMeshingTool(网格工具)命令,然后单击Line域选择所有线条(PickAll),之后用线控制单元网格划分,在No.ofelementdivision(划分单元的份数)中输入10,表示所有线条被划分为10份。本处选用PLANE42单元对盘面划分映射网格。3、定义边界条件并求解建立有限元模型后，就需要定义分析类型和施加边界条件及载荷，然后求解。此处齿轮的载荷为62.8rad/s转速形成的离心力，位移边界条件将内孔边缘节点的周向位移固定，具体分为以下几个步骤。3.1施加位移边界由于此处是对圆柱齿轮进行静态受力分析，为了获得较好的弯曲应力特性，3限制各个面上的法相位移，且将内孔边缘节点的周向位移固定。为各节点其施加周向位移，选择在节点上施加位移约束，并选择UY（Y方向位移），Y方向为周向（此处节点坐标系为柱坐标系，X方向为径向），即施加周向位移约束。3.2施加转速惯性载荷及压力载荷并求解1）要施加齿轮高速旋转引起的惯性载荷，从主菜单中选择SolutionDefineLoadApplyInertiaAngularVelocityGlobal,然后施加角速度：打开ApplyAngularVelocity（施加角速度）对话框，在GlobalCartesianZ-comp（总体Z轴角速度分量）中输入62.8rad/s（即齿轮转速），则由转速引起的惯性载荷施加完成；2）定义载荷时打开选择线的对话框中选择，选择两个相邻的齿边，打开ApplyPRESonlines对话框，依据实际工作环境下所测出的受力数据，在LoadPRESvalue文本框中输入5e6N，施加齿轮啮合产生的压力。定义位移约束并施加载荷后的结果如图3所示，然后选择SolutionSolveCurrentLS进行求解。图3施加的位移约束和载荷4、查看结果4.1查看变形齿轮在高速运转时的关键变形为径向变形，径向变形过大，可能大致边缘与齿轮壳发生摩擦，在DOFsolution（自由度解）选项中选择TranslationX（X向位移），X向位移即为径向位移。结果如下图所示。4图4径向变形图可以看出在边缘处的最大径向位移只有0.497mm，整体变形还是很小的，且沿径向方向上形变递减，说明当齿轮高速旋转时其边缘处变形较大易磨损。4.2查看径向应力齿轮高速旋转时的主要应力也是径向应力，因此要查看该方向上的应力。在Stress(应力）选项中选择X-directionSX(X方向应力）选项，结果如下。图5径向应力分布图5可以看出在齿轮高速旋转时，齿顶处所受应力相对较小仅为0.41E7Pa，而啮合处的应力相对较大，齿根处应力变形更大。在实际工作环境中，由于多次重复的弯曲应力和应力集中，可能会造成轮齿的疲劳折断，故需要提高轮齿的抗断折能力，如：增大齿根过渡圆角半径、采用表面强化处理等。三、离心力分析1、分析问题由于齿轮是轴对称结构，在ANSYS中可以利用结构的周期对称性，在建立模型和求解时，只对一个基本扇区建模和分析，这样可以降低分析的规模，节省计算费用。此处我们单独分析一个轮齿，进而进行离心力分析。2、建立模型2.1定义单元类型此处选用八节点六面体单元SOLID45，不需要设定实常数。8个节点，每个节点3个自由度，x,y,z三个方向。该单元有塑性，徐变，膨胀，应力强化，大变形和大应变能力。（各向异性选用solid64。solid45的高次形式使用solid95。）2.2定义材料属性2.3建立齿轮的一个扇形模型3、定义边界条件并求解考虑到在实际过程中一对啮合齿轮的相互作用。选择轮齿的齿顶面和齿面，在选择面上施加对称边界条件：SolutionDefineLoadsApplyStructuralDisplacementSymmetryB.C.OnAreas。此处单独对高速齿轮受离心力作用下的应力及变形，因此为防止刚性位移，约束关键点47（定义的轮齿边缘的一关键点）Z方向位移。同时施加转速惯性载荷及压力载荷，在SolutionDefineLoadsApplyInertiaAngularVelocityGlobal施加角速度62.8rad/s。选择Sparseiterative求解器进行求解。4、查看结果6图6径向变形图由上图可以看出齿轮径向变形沿径向方向递减，这与对齿轮进行静力分析所得出的结论是一致的。图7径向应力分布图7径向应力分布图显示出齿轮所受应力在内孔边缘处较为集中达到了0.86e11Pa，在齿根处的应力也较大，但没有超出材料的受力极限。四、接触应力分析1、分析问题一对啮合的齿轮在工作时产生接触，分析其接触的位置、面积和接触力的大小。这里我们也对分析模型进行简化，分析平面应变。2、建立模型2.1定义单元类型选用四节点四边形板单元PLANE182，此单元不仅可用于计算平面应力问题，还可以用于分析平面应变和轴对称问题。并且要对单元PLANE182进行设置，在Elementbehavior(单元行为方式)中选择PlaneStress（平面应力）选项。2.2定义实常数此处平面应力行为方式的PLANE182单元，需要设置其厚度实常数为4。2.3定义材料属性。考虑惯性力的静力分析中必须定义材料的弹性模量和密度。2.4建立齿轮面模型2.5对齿面划分网格此处选用PLANE182单元对齿面划分映射网格。PLANE182单元由4个节点定义，每个节点2个自由度，x,y方向。可用于平面单元也可用于轴对称单元。具有塑性，超弹性，应力强化，大变形，大应变能力2.6定义接触对1）在应用菜单中选择UtilityMenu：SelectEntities命令先选择线,选择一个齿轮上可能与另一个齿轮相接触的线(即一个轮齿上的两个齿面);2)在实体选择对话框中选择Nodes,采用“Attachedto”选择方式，在单选列表中选择Lines，all；3）从应用菜单中选择UtilityMenu：SelectComp/AssemblyCreateComponent命令,在Componentname文本框中输入“Node1”定义一个接触节点；4）仿照步骤1）,选择另一个齿轮上可能与前一个齿轮相接触的线；5）仿照步骤3），定义“Node2”，这样节点集合就定义好了。6）选中所有实体（UtilityMenu：SelectEverything）,单击工具条中的接触定义向导按钮，ANSYS会打开ContactManger对话框，选择ContactWizard，分别选择Node1和Node2，单击Create则接触对建立完成。8图8啮合齿轮的接触对模型3.定义边界条件并求解建立有限元模型后，就需要定义分析类型和施加边界条件及载荷，然后求解。此处的载荷为第一个齿轮的转角位移，位移边界条件是第一个齿轮内控边缘节点的径向位移固定，另一个齿轮内孔边缘节点的各个方向位移固定。3.1施加位移边界此处的位移边界条件为将第一个齿轮内径边缘节点的径向位移固定，并将节点坐标系旋转到柱坐标系下施加周向位移。选择第一个齿轮内径上的所有节点，在节点上施加UX（X方向）位移约束。3.2施加第一个齿轮位移载荷及第二个齿轮的位移边界条件并求解此处要施加第一个齿轮内径节点的位移载荷，选择第一个齿轮内径上的所有节点，施加周向位移约束，将Displacementvalue设定为-0.2。然后选择第二个齿轮内径上的所有节点，在节点上施加各个方向位移约束（ALLDOF），将Displacementvalue设定为0。在分析类型中将求解控制选项，选择LargeDisplacementStatic，在TimeatendofLoadstep文本框中输入1，在Numberofsubsteps文本框中输入20，然后求解。4、查看结果4.1查看vonMises等效应力如下图所示9图9齿轮啮合vonMises等效应力图由上图可以看出一对接触齿轮啮合时应力集中在内孔边缘和轮齿齿根部分，但整体的受力很小。4.2查看接触应力如下图所示图10啮合齿轮接触应力10如图所示，齿轮啮合的接触应力分布呈对称状态，在接触面的应力较集中，但数值较小。在实际工作环境中，当齿面材料硬度较低，载荷较重时齿轮传动啮合时齿面相对滑动速度方向在节点处发生改变。在主动轮齿面的节线两侧，齿顶和齿根的摩擦力方向相背，在节线附近形成凹槽，在从动轮齿面节线附近形成凸脊，从而发生齿面塑性变形。五、分析综合结果通过以上一系列的分析可以发现齿轮在工作状态下的形变程度较小，接触应力在其极限范围内，能够满足要求，符合工程安全性。本文重点是对高速齿轮进行静力学分析，当然对于一对啮合齿轮，还需要考虑其动载荷。在实际工作环境下由于制造误差和弹性变形等原因，基圆齿距不可能完全相等，这齿轮啮合时的瞬时速比发生变化而产生冲击和振动载荷。齿轮的圆周速度越大、加工精度越低，齿轮动载荷越大。因此针对高速齿轮，为使其正常工作并延长其工作寿命，可以采取以下措施：1）适当提高制造精度、降低齿轮圆周速度、增加轮齿及支承件的刚度，对齿轮齿顶的一小部分齿廓曲线进行适当修削，都能减少内部附加载荷；2）提高齿面硬度，降低表面粗糙值；3）采用粘度较高，抗胶合能力强的润滑油。参考文献[1]张文志.机械结构有限元分析.哈尔滨：哈尔滨工业大学出版社，2006，第1版[2]王庆五，左昉，胡仁喜.ANSYS10.0机械设计高级应用实例，北京：机械工业出版社，2006，第2版[3]SaeedMoaveni.有限元分析-ANSYS理论与应用.北京：电子工业出版社，2013，第3版[4]叶东友.基于ANSYS的渐开线支持圆柱齿轮有限元分析[J].煤矿机械，2004,06:43-45[5]刘海娥，张思婉.基于ANSYS的渐开线啮合齿轮有限元分析[J].机械制造与自动化，2010,01:28-29