基于apdl语言的齿轮瞬态响应分析.docxVIP

基于apdl语言的齿轮瞬态响应分析 齿裂齿轮的特殊形状是齿轮压力的一种特殊形式。适用于低速和重载机械传动。在泥浆和磨损条件差的工作环境中，它被广泛应用于重量缓慢的大型设备中，如重量、矿山、冶金、化工甚至娱乐行业。 当前,普遍使用的销齿传动齿轮的弯曲强度分析计算方法仍然是基于传统材料力学强度理论和假设的强度计算方法,此方法是否实用有效、精确,有待进一步研究。目前有关基于FEM的销齿副齿轮静应力分析以及模态分析已见报道。事实上,齿轮所承受的载荷是交变载荷,齿根是在脉动循环(单向转动)或对称循环(双向转动)动应力的作用下工作。有关销齿传动齿轮弯曲动强度的精确分析计算迄今为止仍未见报道。为了在现有基础上开发出性能更加完善的销齿传动机构,缩短产品开发周期,提高产品质量,本文抛开传统的基于悬臂梁理论的齿根弯曲应力的计算方法,根据瞬态响应特性分析理论,借助于有限元分析软件ANSYS10.0对某一特定工况下的销齿传动齿轮齿根弯曲动应力进行瞬态响应特性分析,得出齿根等效弯曲动应力的分布区域和变化规律,为其动态性能的进一步优化奠定基础。 1 振动体系的独立参数 瞬态响应特性分析是用于确定承受任意随时间变化载荷的结构的动力学响应的一种方法。使用这种方法确定结构在任意随时间变化载荷作用下的应力、应变和位移等。 齿轮的质量是连续分布的。确定其体系全部质量位置所需要的独立参数(即自由度)将有无穷多个。按照有限元的思想,将无穷多个自由度的振动体系近似地以有限个自由度来代替。像静力学问题一样,将结构离散成有限个单元,按有限元的集成方法,建立结构瞬态响应特性分析,其基本方程为 [M]{ü}+[C]{U˙}+[K]{U}={F(t)}.[Μ]{ü}+[C]{U˙}+[Κ]{U}={F(t)}. 式中:[M]为质量矩阵,[C]为阻尼矩阵,[K]为刚度矩阵,{ü}为节点加速度向量,{U˙U˙}为节点速度向量,{U}为节点位移向量,{F(t)}为载荷向量,是时间任意函数。 2 描述和分析实例 2.1 齿轮齿长度设计 齿轮齿数Z1=13,销轮齿数Z2=360,齿距p=63.16 mm, 齿轮节圆直径d1=261.4 mm,齿轮齿根圆角半径ρf=15.5 mm,齿轮齿根圆角半径中心至节圆圆周距离c=1.5 mm,齿轮齿顶高ha=27.16 mm ,齿轮全齿高h=44.16 mm,齿轮齿廓过渡圆弧半径R=10 mm ,齿轮齿顶圆半径da1=315.7 mm,齿轮齿根圆直径df1=227.4 mm,中心距a=3 749.5 mm,齿轮齿宽b=45 mm,销轮功率P2=2.6 kW,销轮轴转矩T2=5 2461 N·m,重合度ε=1.28,弹性模量E=207 GPa,泊松比μ=0.27,单向传动。 2.2 齿合合度的计算 本例中,齿轮转速n1=10.356 r/min=62.136°/s,齿轮每转1周所用时间t1=5.793 7 s,每齿啮合的时间t2=0.445 673 s。根据啮合原理,由于重合度ε=1.28,因此,两对齿同时啮合的时间t3=0.28t2=0.124 8 s,单齿啮合的时间t4=0.72t2=0.320 9 s。单向传动,齿根弯曲应力循环特性是脉动循环。 3 几何简化和网格细化 根据实际结构,先对实体模型进行适当简化。为了避免有限元网格的尺寸大小相差太悬殊而影响单元质量和计算精度,对一些影响齿轮强度极微的小倒角和圆角进行简化处理。同时,为了节约计算资源,用固定尺寸划分整体网格,在此基础上,将加载齿表面进行网格细化。为了能精确加载,在创建有限元模型时,在加载点处创建硬点。选用SOLID95单元进行离散分网,此单元是分析弹性结构空间问题中应用较广的一种元素。由于采用了20节点单元,因而能利用更复杂的形状函数,对实际变形有更好的表达,计算精度更高,得到16 985个单元,28 332个节点。有限元模型如图1所示。 4 应立即分析噪音 4.1 给安装额外的限制 根据齿轮的实际工作状况,对齿轮内孔表面施加完全约束。 4.2 数据表的加载 由于作用在轮齿上的载荷是随时间移动的,因此,采用数据表的方式加载,用以表达载荷的真实特性。在一个工作循环内取5个加载时刻,共制作5个数据表,如表1所示。 4.3 应立即分析噪音 在求解器中,用完全法进行瞬态响应分析。设定求解结束时间0.436 8 s,加载时间间隔为0.000 03 s,采用斜坡式加载。 5 弯曲应力分布 1)通过在通用后处理器中查看各个时刻等效应力的分布情况,得知当时间t=0.383 3 s时,动应力有最大值,其值为78.65 MPa,应力分布情况如图2所示。当按传统计算方法进行计算时,得轮齿计算弯曲应力为81.61 MPa。与之相比,小3.75%。可见我国目前采用的计算方法可用,但略显保守。 2)在创建几何模型及