斜齿轮的优化设计与有限元分析

《斜齿轮的优化设计与有限元分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、现代设计方法三级项目报告斜齿轮的优化设计与有限元分析姓名：课程名称：现代设计方法指导教师：-------------------------------------------------------来自燕大2013年5月摘要：目录1任务分工12问题描述13基于matlab的斜齿轮参数优化13.1目标函数的建立13.2约束条件的建立23.3建立数学模型33.4斜齿轮参数34基于ansys的斜齿轮有限元分析44.1网格划分本44.2加载54.3受力分析64.4分析结果65总结76参考文献7斜齿轮的优化设计与有限元分析徐航，赵航，骆华玥（燕山大学机械工程学院）摘要：本文利用m

2、atlab和ansys对二级同轴斜齿轮减速器进行了优化设计。通过对中心距的优化得到了最理想的齿轮参数，即在满足使用强度的前提下，最大限度的降低了成本。1任务分工徐航负责Matlab与Ansys的模拟仿真赵航负责模型的建立及数值的分析计算骆华玥负责演示文稿与说明书的制作。2问题描述齿轮减速器广泛应用在煤炭、机械等行业,传统设计全由设计人员手工完成，但在性能更好、使用更可靠方便、成本更低、体积或质量更小的指标要求下，希望能从一系列可行的设计方案中精选最优，传统的设计方法做不到，因而有必要采用matlab优化方法来确定其设计参数。再运用Ansys软件来对其进行受力模拟，通过An

3、sys就可以辨别优化方案的可靠性，对其进行筛选，通过Matlab与Ansys软件的共同使用就可以对方案进行提前鉴别，避免了不必要的损失，更有利于资源的优化使用和效益的产生。3基于matlab的斜齿轮参数优化3.1目标函数的建立据优化目标的不同，齿轮减速器设计可以有多种最优化方案，文中讨论的是在满足齿轮传动强度、刚度和寿命条件下，使减速器体积最小或质量最小。显然，若减速器结构紧凑，则其重量和体积为最小，而结构的紧凑与否，关键在于减速器的总中心距，因此以总中心距最小为优化目标，建立优化设计数学模型。二级斜齿圆柱齿轮减速器总中心距A的数学表达式为式中mn12，i12和mn34,

4、i34———高速级和低速级齿轮的法向模数和传动比 Z1，Z3———高速级和低速级小齿轮的齿数 β———斜齿轮螺旋角因为总传动比i已知，则i12=i34=√2。又因为是同轴减速器mn12=mn34，Z1=Z3，β12=β34。所以目标函数有3个独立的设计变量： 令f（x）=A，所以目标函数的表达式是： 6要求解的是目标函数f（x）的最小值，即minf（x）。1.1约束条件的建立（1）齿根弯曲疲劳强度计算式中K———齿轮传动的载荷系数； T1———主动轮的转矩； YFa，YSa，Yε，Yβ———齿形修正系数，应力修正系数，重合度系数，螺旋角系数 σF———齿轮许用弯曲应力。

5、（2）齿面接触疲劳强度式中ZE，ZH，Zε，Zβ———材料弹性系数，节点区域系数重合度 系数，螺旋角系数； d1———小齿轮分度圆直径； [σF]———齿轮许用接触应力。（3）将齿根弯曲疲劳强度条件、齿面接触疲劳强度条件分别按高速级大小齿轮、低速级大小齿轮进行参数化，得到以下6个强度约束条件，均为非线性约束。①高速级齿轮齿根弯曲疲劳强度约束条件②高速级齿轮齿面接触疲劳强度约束条件（4）其他约束条件考虑到齿轮传动平稳、斜齿轮轴向力不可太大、满足短期过载条件、高速级与低速级大齿轮浸油深度大致相近、齿轮分度圆直径不能太小等因素，可建立以下隐性约束，均为线性的。 a1≤mn12≤

6、b1,a2≤Z1≤b2,a3≤β≤b3 6可根据经验确定设计变量的上下边界bj和ai： 1≤x1≤6,20≤x2≤40,8≤x3≤301.1建立数学模型把上述数据分别代入目标函数及约束条件，可得，目标函数标准式约束条件标准式：1≤x1≤620≤x2≤408≤x3≤301.2斜齿轮参数1.工件材料：小齿轮选用45钢，调质，HB1=240HBS;大齿轮选用45钢，正火，HB2=200HBS；齿轮精度8级；传动效率取i=12.6，i12=i34=3.5；载荷系K=2.2；齿宽系数ψd=b/d1=0.8。转矩：T1=22.13N•mT2=73.60N•mT3=245.00N•m。

7、2.齿轮齿面接触许用应力3.齿轮齿根弯曲许用应力4.修正系数5.matlab编程计算采用MATLAB软件优化工具箱求解最优结果，进行非线性有约束多元函数最小值计算[4]，命令函数为fmincon。 函数m文件function[c,ceq]=m(x) c(1)=259.34*cos(x(3)/180*pi)/(0.8*x(1)^3*x(2)^2+6*x(1)^2*x(2))-450 c(2)=862.59*cos(x(3)/180*pi)/(0.8*x(1)^3*x(2)^2)-390 c(3)=354.35*sqrt(125.