二级直齿圆柱齿轮的优化设计.doc

济南大学泉城学院毕业设计

- 24 -
1 前言
1.1 课题来源以及背景
减速器是原动机和工作机之间独立的闭式机械传动装置,用来降低原动机转速或增大转矩,以满足工作机需要，是一种被广泛在机械制造机床，汽车，矿山机械等不同机器的机械部件。国外的减速器，以德国、丹麦和日本处于领先地位，特别在材料和制造工艺方面占据优势，减速器工作可靠性好，使用寿命长。但其传动形式仍以定轴齿轮传动为主，体积和重量问题，也未解决好。最近报导，日本住友重工研制的FA型高精度减速器，美国Alan-Newton公司研制的X-Y式减速器，在传动原理和结构上与本项目类似或相近，都为目前先进的齿轮减速器。当今的减速器是向着大功率、大传动比、小体积、高机械效率以及使用寿命长的方向发展。因此，除了不断改进材料品质、提高工艺水平外，还在传动原理和传动结构上深入探讨和创新，平动齿轮传动原理的出现就是一例。减速器与电动机的连体结构，也是大力开拓的形式，并已生产多种结构形式和多种功率型号的产品。目前，超小型的减速器的研究成果尚不明显。在医疗、生物工程、机器人等领域中，微型发动机已基本研制成功，美国和荷兰近期研制的分子发动机的尺寸在纳米级范围，如能辅以纳米级的减速器，则应用前景远大。
国内的减速器多以齿轮传动、蜗杆传动为主，但普遍存在着功率与重量比小，或者传动比大而机械效率过低的问题。另外，材料品质和工艺水平上还有许多弱点，特别是大型的减速器问题更突出，使用寿命不长。国内使用的大型减速器（500kw以上），多从国外（如丹麦、德国等）进口，花去不少的外汇。60年代开始生产的少齿差传动、摆线针轮传动、谐波传动等减速器具有传动比大，体积小、机械效率高等优点。但受其传动的理论的限制，不能传递过大的功率，功率一般都要小于40kw。由于在传动的理论上、工艺水平和材料品质方面没有突破，因此，没能从根本上解决传递功率大、传动比大、体积小、重量轻、机械效率高等这些基本要求。
减速器设计应在满足需要的前提下追求体积小、尺寸小、重量轻、成本低、润滑方便等优化目标。圆柱齿轮机构是应用最为广泛的一种传动机构。传统圆柱齿轮机构的设计方法大多是依靠分析、试凑或类比的方法来确定复杂的结构参数，这就造成了可行设计方案的纰漏，使设计变得相对被动，往往需要多次重复性的工作，才能得到较满意的结果。随着计算机技术的快速发展和优化理论的日趋成熟，最优化设计在齿轮设计上得到了实际的应用。
1.2 初步设计方案的确定
圆柱齿轮减速器的的设计方法是：通过查阅资料，结合实习所得的资料文献，运用自己平时的理论知识，对以往设计的减速器进行对比，第一，通过对优化设计的理论理解，在理论上对体积进行优化，然后进行设计变量的选取，建立数学模型，最后一步是以带式运输机的二级直齿圆柱齿轮为例，具体的结合数据，建立起约束条件，利用MATLAB软件对将目标函数和约束条件编写进程序，运行MATLAB软件对数据进行优化，得到优化的参数。这样设计方案就初步确定了。
1.3 初步利用数学建模进行优化
由于所有优化的设计变量较多,在优化减速器参数时注意属于整数规划的范畴。但目前对于整数规划的理论研究还不充分,求解方法比较有限,就工程中广泛使用的圆柱二级级圆柱齿轮减速器的优化设计来说，由于齿轮系统优化设计时涉及的参数较多，计算比较复杂，因此，为了减低制造减速器的成本，提高减速器的使用效率，对二级圆柱齿轮进行优化设计是未来设计齿轮的趋势，通过建立数学模型，对设计变量进行约束，应用一般的开发语言，在编程时工作量很大，在数值方面没有明显优势，采用MATLAB 优化工具箱，将上述的模型和公式编入程序利用MATLAB进行求的最优解。随着计算机在工程制造业中的广泛应用，MATLAB作为一款数学软件，在进行数学公式和计算中更容易，更快捷，在工程设计中越来越多的使用MATLAB进行优化设计。
1.4 本章小结
本章主要是阐明课题来源，了解到减速器的广泛用途，但是通过国内外现状的比较，让我们认识到对二级减速器进行优化的可行性和必要性，然后初步的构建起初步方案，利用MATLAB软件进行优化设计。

2 优化设计的基本理论
优化方法在机械设计中有着广泛的应用，构成了机械优化设计体系。机械设计领域设计的面很宽，如机械结构设计，零件设计，部件设计以及整机等产品设计，所需考虑的方面也很多如结构组成与分析、机构运动学以及动力学问题，作为一项好的设计，要使产品具有好的性能，又要满足工艺要求；要达到预期的精度要求，并使用安全可靠，特别是当设计的参数很多，各个参数之间的关系很复杂的情况下，势必需要采取现代设计手段，设计者从大量繁琐的计算中解脱出来。因此，优化设计已经成为当设计中的一个重要方面，并在大量工作的基础上取得了丰硕的成果。
优化设计是从多种方案中选择最佳方案的设计方法。它以数学中的最优化理论为基础，以计算机为手段，根据设计所追求的性能目标，建立目标函数，在满足给定的各中优化条件下，寻求最优的设计方案。然而机械优化设计是最优化技术在机械设计领域的移植和应用,其基本思想是根据机械设计的理论,方法和标准规范等建立一反映工程设计问题和符合数学规划要求的数学模型,然后采用数学规划方法和计算机计算技术自动找出设计问题的最优方案.优化设计是一种格式化的设计方法。对于各式各样的设计问题，都必须先将它按照规定的格式要求建立优化计算的数学模型，然后再在计算机上用优化设计方法解出他的最优解。
2.1优化设计方法解决实际问题的步骤
2.1.1 分析实际问题，建立优化设计的数学模型
分析：①设计的要求、目标和准则。
②设计的限制约束条件。
③设计的参数，确定设计变量。
建立：机械设计优化方法相应的数学模型
2.1.2 数学模型的分析与优化
建立优化数学模型，将工程实际问题以数学模型。将工程实际问题以数学模型表示。为此要恰当学选取设计变量，将设计问题所追求的最佳设计指标与设计变量之间的关系用函数式或其他的方式表示出来，以上三方面构成了优化设计的数学模型，选取优化设计的方法。目前可功优化设计工作使用的优化方法很多，其中包括无约束优化方法和约束优化方法，要根据不同的设计规模选用恰当的优化方法，以使得达到计算的高效率，以期具有预期的设计精度。
2.1.3数学模型的最优解及圆整
工程中的最优设计，是在考虑诸多影响因素的条件下，以得到最佳的设计方案，即最佳设计参数值。用计算机进行求解，优化设计是综合有关各个方面的因素，按优化原理进行求解，是建立在近代数学以及计算机广泛应用基础上的一项新技术，所以以计算机为工具，以人机配合的方式进行自动搜寻最优值。
2.2开始建立优化设计的数学模型
建立优化设计的数学模型必须确定设计变量、目标函数、和约束条件这三个必要的部分。
2.2.1 首先确定设计变量
在构成一项设计方案的全部参数中，可能有一部分参数根据事件情况预先确定他的数值，它们在优化算法过程中始终保持不变，这样的参数成为设计变量，在机械设计中有些参数只能选用规定的离散值，如齿轮的模数、钢材的规格尺寸等这样的参数作为设计变量叫作离散设计变量。
2.2.2 其次确定目标函数
在所列举的设计问题中，如弹簧的体积，剪刀的轨迹误差等都是设计中所追求的目标，而且它是设计变量的函数，成为目标函数。优化设计的目的就是要求合理选择设计变量值使目标函数达到最佳值。由于目标函数是设计变量的函数的函数，故给定一组设计变量就对应的有一个函数值。在解决优化设计问题时，正确选择目标函数是非常重要的，它不仅直接影响优化设计的结果，而且对于整个优化设计的繁简难易也会有一定的影响。
2.2.3最后确定约束条件
优化设计不仅要使选择方案的设计指标达到最佳值，同时还必须满足一些附加的设计条件，这些附加的设计条件是对设计变量却值得相互关系及其大小加以限制，在优化设计中叫约束条件。根据约束性质的不同可以将约束条件分为边界性约束和性能性约束两类。在解决工程问题时，约束条件是优化设计获得工程可接受设计方案的重要条件。不等式约束以及有关概念相当重要的，并且要注意在满足设计要求为原则时避免重复矛盾和线性相关的约束。
2.3序列二次规划法（SQP）函数fmincon
2.3.1序列二次规划法
序列二次规划法又称为逼近法，有约束问题的变尺度法，或拉格朗日牛顿法，为非线性数学规划法，用于求有约束的最优化问题，具有整体收敛性且同时保持局部超一次收敛性，被公认为是当今求解光滑的非曲线性规划问题，最优秀的算法之一，数值计算表明序列二次规划法优于乘子法。序列二次规划法的基本思路是在当前的迭代点处，利用目标函数的二次近似和约束函数的一次近似构成二次规划，通过求解这个二次函数获得下一个迭代点， k-t方程的解形成了许多线性规划算法的基础，通过用类比牛顿法更新过程，给k-t方程积累二节信息，可以保证有约束拟牛顿法的超线性收敛，这样的算法成为序列二次规划法。
对于规划问题

(2.1)

用Kuhn-thcher方程表示为

(2.2)

上述描述了目标函数和约束梯度之间的函数关系。
2.3.2优化函数 fmincon
（1）fmincon函数提供了大型优化算法和中型优化算法。使fmincon函数，将原问题转化为一系列的二次规划问题进行求解，对于每一种函数每一步的求解都是通过选择一种最佳方法来进行。默认时，若在函数中提供了梯度，并且只在上下界存在或只有等式约束时，fmincon函数将使用大型算法，当即又等式约束又有梯度约束时，使用中型算法。
（2）fmincon函数中的中型算法使用的是序列二次规划法。在每一次迭代中求解二次子规划问题，并用BFGS更新拉格朗日矩阵。
（3）fmincon函数可能会给出局部最优解，这与初值的选取有关系。
Fmincon函数的语法与说明
[x, fval]=fmincon(@fun, x0, A, b, Aeq, beq, lb, ub, @nonlcon, options)
fval：返回目标函数在最优解x点的值；
fun：目标函数；
nonlcon：约束值函数；
options：设置优化选项参数。
2.4 本章小结
本章在优化设计的理论上阐述了对二级直齿圆柱齿轮减速器进行优化的可行性，进行优化设计包含的三个内容。第一，建立优化数学模型；第二，选取优化方法；第三，利用计算机进行求解。而且将序列二次规划法和MATLAB软件进行结合，利用MATLBA优化工具箱进行求解。

3 二级直齿圆柱齿轮的优化设计

图3.1 齿轮传动简图

图3.1表示的二级直齿圆柱齿轮的传动系统简图，具体采用的是小齿轮采用齿轮轴的形式，因为大齿轮的直径较大，已经大于160mm,所以采用腹板式结构。采用等变为传动。

图3.2 大齿轮的结构和各部分的尺寸

由于低速轴也采用小齿轮为齿轮轴形式，大齿轮采用腹板式结构，所以都可以用上图表示出来。对于大齿轮的机构我们采用一下计算公式：
(3.1)

其中 m表示齿轮的模数，Z表示齿轮的齿数，b表示齿轮的齿宽，表示轴径。
3.1已知参数模型的建立

图2.5 插图样式

图3.3 带式传送机

图 3.4 传动装置简图
如图3所示，二级圆柱直齿圆柱齿轮在工矿业和农业生产中有着广泛应用，改图就是在农业和工业中运用比较广泛的带式输送机，对于带式输送机大家可能不陌生，它伸缩胶带输送机与普通胶带输送机的工作原理一样，是以胶带作为牵引承载机的运输设备，不过在增加了储带装置和胶带装置等，当游动小车向机尾一端移动时，胶带进入储带装置内，机尾回缩，反之，则机尾延长，因而使输送机具有可伸缩性能。它的机构组成由固定部分和非固定部分两大部分。固定部分由机头传动装置、储带装置、收放胶带装置等组成；非固定部分由无螺栓连接的快速可拆支架、机尾等组成。机头传动装置由传动卷筒、减速器、液力联轴器、机架、卸载滚筒、清扫器组成。所以由此可见减速器是输送机的重要组成部分，对输送机的传动系统进行优化时，首先要对减速器进行优化，我们就以输送机为实例进行优化设计。
已知二级齿轮减速器的输出功率是2.6KW，高速轴转速是1430r/min,总的传动比为15，小齿轮的需用接触应力[]=540MP,许用弯曲应力为[]=303.57MP,[] =238.86MP。齿轮的材料和热处理方式：大齿轮采用45号钢正火HB=187-207,小齿轮采用45号钢调质HB=228-255,总工作时间不少于十年。这样需要优化的参数模型就建立起来了。
3.2 优化目标
在二级齿轮传动系统中，我们所追求的性能指标是比较多的，比如传递功率要求最大，中心距要求最小，要求二级直齿圆柱齿轮传动的质量最轻，也就是要求齿轮的体积最小，所以我们在满足各项几何条件和传动质量指标的条件下，我们以体制最小为优化目标，并且其他的目标也会在体积最小这个目标中，就是确定了体积最小为目标也就让其他的性能指标得到了优化，所以我们在保证承载能力的条件下，以体积最小为目标进行二级直齿圆柱齿轮的优化设计，这样优化的目标也就明确了。
3.3 设计变量的选取
对于优化设计问题，合理地选择设计变量，才能得到最优化结果。第一级齿轮传动的设计参数有模数m1，齿数z1、z2，齿宽b1，为了提高齿轮的承载能力或满足中心距要求，常采用变位齿轮，因此变位系数y1和y2也要作为设计参数。在传动比已知的情况下，由于存在z2= z1·i12关系式，所以z1、z2不是一组独立的变量，取z1作为设计变量。以体积为目标函数，因此大齿轮轴径dzh2也作为设计变量。同样，第二级传动的设计变量为m2，z3，y3，y4，b3，dzh4。由于总的传动比已经确定，由i=i12·i34，所以i12、i34不是一组独立的变量，取i12作为设计变量。将表达目标函数和约束条件所必须的独立参数——设计变量列下：
(3.2)

式中X—设计向量；
——表示x的各个分量；
——表示高，低速齿轮的传动模数；
——表示高低速级小齿轮的齿数；
——表示齿轮传动的变为系数；
——表示两个大齿轮的齿宽；
——表示两个大齿轮的轴径；
——表示高速级齿轮的传动比；
这样设计变量就确定了。
3.4目标函数的建立
3.4.1小齿轮的体积
由于小齿轮采用齿轮轴的形式，因此其公式应为:

小齿轮1的体积： (.3.3)
小齿轮3的体积： (3.4)

3.4.2 大齿轮的体积
由于大齿轮采用腹板式，所以其体积公式由图2可得到:
大齿轮2的体积为：

(3.5)
(3.6)
大齿轮4的体积为：
(3.7)
3.4.3 二级直齿圆柱齿轮传动系统优化设计的目标函数

(3.8)

这样通过求各个齿轮的体积，然后加起来就是总的体积，得到了传动系统的目标函数。
3.5 约束条件的建立
3.5.1 变量上下限的确定
⑴为避免发生根切最小齿数；
⑵根据装配工艺的条件和以往的经验，应该让小齿轮1的直径不超过45mm，小齿轮3的直径不超过75mm；
⑶综合考虑传动的平稳性，轴向力不可太大，为保证齿轮的承载能力同时又避免载荷齿宽分布严重不均，要求齿轮系数应该在0.9到1.4之间，所以就要求

(3.9)

⑷对于传递动力的齿轮，一般模数 ，且应该取标准值；
⑸根据以往的经验，大齿轮的轴径应该在30mm到60mm之间即
(3.10)
⑹变位系数的确定

由里西特公式得到 (3.11)

得到变位系数的取值范围为：
(3.12)

⑺由于减速箱是展开式布置，所以
所以通过上面所叙述的可确定各个变量的上下限为：
(3.13)

3.5.1约束条件的确立
⑴标准齿轮传动具有设计简单，互换性好等一系列优点。但其中也有不足之处，例如要求齿轮齿数最少齿数不少于十七，否则将产生根切现象；标准齿轮不适用于中心距
的场合。因为当时，无法安装；而当时，而又会产生过大的齿侧间隙，影响传动的平稳性，而重合度也随之降低；另外，在一对相互啮合的标准齿轮中，由于小齿轮齿廓渐开线的曲率半径较小，齿根厚度也比较薄，参与啮合的次数又较多，强度较低，影响到整个齿轮传动的承载能力。为了改善标准齿轮的上述不足之处，就必须突破标准齿轮的限制，对齿轮进行修正。所以变为齿轮就出现了。若齿轮采用等变为传动得到
； (3.14)
⑵由工艺装备条件，， (3.15)
⑶由齿宽系数条件

， (3.17)
； (3.18)

⑷所设计的齿轮传动在具体的工作情况下，必须具有足够的、相应的工作能力，以保证在整个工作寿命期间不致失效。因此，针对上述各种工作情况及失效形式，都应分别确立相应的设计准则。但是对于齿面磨损、塑形变形等，由于尚未建立起广为实际工程使用而且行之有效的设计方法及设计数据，所以目前设计一般使用的齿轮时，通常只按保证齿根弯曲疲劳强度及保证齿轮接触疲劳强度两准则进行计算。对于告诉大齿轮的齿轮传动，还要保证齿面抗胶合能力的准则进行计算。至于抵抗其他失效的能力，目前一般不进行计算，但应该采用相应的措施，以增强齿轮抵抗这些失效的能力。
有实践可得，在闭式齿轮穿洞中，通常保证吃面接触疲劳为主。但对于齿面硬度很高齿心强度又低的齿轮或者材质较脆的齿轮，通常则保证齿根弯曲疲劳强度为主。
⑸由齿面接触疲劳强度条件可得到：

, ； (3.19)
由齿面接触疲劳强度计算的基本公式可得到：
(3.20)

上式各参数含义及其取值：K——为载荷系数，适选载荷系数K=1.3；
-为载量的弹性影响系数，由表10-6查（《机械设计》第八版）的材料的弹性影响系数；
-为齿轮的节点区域系数，（标准直齿轮时，）；
-为小齿轮传递的名义转矩；
计算小齿轮1的传递转矩：
(3.21)
同理 通过计算可得到低速轴的转速(3.22)
由于存在功率损失，通过功率损失的计算可得到低速轴的输入功率为2.6KW；
(3.23)
由题意可得到，所以可得到
则 (3.24)
低速轴齿轮的

则 (3.25)
⑹由齿根弯曲强度条件可得到：
；
(3.26)
由齿根弯曲疲劳强度计算的公式可得到：
(3.27)

上式各个参数的含义及其取值：K——为载荷系数，适选载荷系数K=1.3；
-为小齿轮传递的名义转矩；
YF——为齿形系数，齿形系数的计算公式如下：

；
；

(3.28)

具体的计算齿根弯曲应力步骤和建立不等式约束：

(3.29)
②由上面所得到以下不等式：

(3.30)

通过上述不等式我们就建立起了约束条件。
3.6 选用MATLAB优化工具箱对数据模型进行优化
由于模型类型是非线性规划问题，所以采用MATLAB工具箱中的非线性规划函数fmincon.
3.6.1 MATLAB的概述
对于MATLAB, 它是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近Windows的标准界面，人机交互性更强，操作更简单。而且新版本的MATLAB提供了完整的联机查询、帮助系统，极大的方便了用户的使用。简单的编程环境提供了比较完备的调试系统，程序不必经过编译就可以直接运行，而且能够及时地报告出现的错误及进行出错原因分析。
3.6.2 MATLAB语言概述
(1)NEW及其子菜单：允许用户打开一个心的文件，新的图形窗，仿真模型文件，和图形用户界面文件；
(2)OPEN:从指定的相应文件的路径和文件中打开一个已经存在的文件；
(3)CLOSE COMMAND WINDOW:关闭命令窗口；
(4)IMPORT DATA:在MATLAB工作空间中生成一变量，并从指定的文件中和路径中获取数据；
(5)SAVE WORKSPACE AS:将工作空间中的所有变量数据保存在指定的路径下的相应的文件中；
(6)SET PATH:设置MATLAB的搜索路径；
(7)PREFERENCE:允许用户对系统的一些性能参数进行设置，如数据格式、字体大小与颜色等。
3.6.3 MATLAB常用的函数
(1)．在进行二级直齿圆柱齿轮优化设计中涉及到的函数用语；
eps：浮点相对精度
exp：自然对数的底数e
i 或 j：基本虚数单位
inf 或 Inf：无限大， 例如1/0
nan或NaN：非数值（Not a number），例如0/0
pi：圆周率 p（= 3.1415926...）
realmax：系统所能表示的最大数值
realmin：系统所能表示的最小数值
nargin: 函数的输入引数个数
nargout: 函数的输出引数个数
lasterr：存放最新的错误信息
lastwarn：存放最新的警告信息
(2).MATLAB常用基本数学函数
abs(x)：纯量的绝对值或向量的长度
angle(z)：复数z的相角(Phase angle)
sqrt(x)：开平方
real(z)：复数z的实部
imag(z)：复数z的虚部
conj(z)：复数z的共轭复数
round(x)：四舍五入至最近整数
3.6.4 利用MATLAB开始编程求解
目标函数fun定义为：
function f=GearFunct(x)
f=(pi*x(9)/4)*(x(1)*(x(3)+x(5)))^2+(pi*x(10)/4)*(x(2)*(x(4)+x(7)))^2+(pi*x(9)/4)*((x(1)*(x(13)*x(3)+x(6)))^2-(x(1)*x(3)*x(13)-8*x(1))^2)+(pi*x(9)/4)*(0.3*((x(1)*x(3)*x(13)-8*x(1))^2-(1.6*x(11))^2)+1.872*x(11)^2)+(pi*x(10)/4)*((x(2)*(15/x(13)*x(4)+x(8)))^2-(x(2)*15/x(13)*x(4)-8*x(2))^2)+(pi*x(10)/4)*(0.3*((x(2)*15/x(13)*x(4)-8*x(2))^2-(1.6*x(12))^2)+1.872*x(12)^2);函数nonlcon定义为：
function[C Ceq]=GearConstr(x)
C(1)=x(1)*x(3)-30;
C(2)=x(2)*x(4)-30;
C(3)=0.9-x(9)/(x(1)*x(3));
C(4)=x(9)/(x(1)*x(3))-1.4;
C(5)=0.9-x(10)/(x(2)*x(4));
C(6)=x(10)/(x(2)*x(4))-1.4;
C(7)=(3390.2/(x(1)*x(3)))*sqrt((x(13)+1)/(x(9)*x(13)))-540;
C(8)=(12349/(x(2)*x(4)))*sqrt((15/x(13)+1)/(15/x(13)*x(10)))-540;
C(9)=51.0484/(x(1)^2*x(3)*x(9)*(0.169+0.006666*x(3)-0.000845*x(3)^2))-303.57;
C(10)=51.0484/(x(1)^2*x(3)*x(9)*(0.2824+0.0003539*x(13)*x(3)-0.000001576*x(3)^2*x(13)^2))-238.86;
C(11)=677.3/(x(2)^2*x(4)*x(10)*(0.169+0.006666*x(4)-0.000845*x(4)^2))-303.57;
C(12)=677.5/(x(2)^2*x(4)*x(10)*(0.2824+0.005309*x(4)/x(13)-0.0003546*x(4)^2/x(13)^2))-238.86;
Ceq=[ ];
主程序：
clc
x0=[2, 3.5, 21, 20, 0.3, -0.3,0.4,-0.4,42,70,40,52,4.8]; %传统设计值
lb=[1.75,4,20,20,0.3,-0.5,0.3,-0.5,40,65,30,40,4];
ub=[2.5,5,25,25,0.5,-0.3,0.5,-0.3,50,75,50,60,5];
Aeq=[0, 0, 0,0,1,1,1,1,0,0,0,0,0];
beq=[0];
options = optimset('LargeScale','off');
[x, fval] = fmincon (@GearFunct, x0, [ ], [ ], Aeq, beq, lb, ub, @GearConstr, options)执行结果x=Columns 1 through 11 2.0000 2.8043 21.0000 18.8043 0.3000 -0.3000 0.4000 -0.4000 42.0000 70.0000 40.0000 52.0000 4.8000 fval = 1.6718e+006
优化之前的主程序：clc
x=[2, 3.5, 21, 20, 0.3, -0.3,0.4,-0.4,42,70,40,52,4.8];
f=(pi*x(9)/4)*(x(1)*(x(3)+x(5)))^2+(pi*x(10)/4)*(x(2)*(x(4)+x(7)))^2+(pi*x(9)/4)*((x(1)*(x(13)*x(3)+x(6)))^2-(x(1)*x(3)*x(13)-8*x(1))^2)+(pi*x(9)/4)*(0.3*((x(1)*x(3)*x(13)-8*x(1))^2-(1.6*x(11))^2)+1.872*x(11)^2)+(pi*x(10)/4)*((x(2)*(15/x(13)*x(4)+x(8)))^2-(x(2)*15/x(13)*x(4)-8*x(2))^2)+(pi*x(10)/4)*(0.3*((x(2)*15/x(13)*x(4)-8*x(2))^2-(1.6*x(12))^2)+1.872*x(12)^2)
执行的结果为：f =2.2967e+006
因为模数m应该为标准值，齿数z应该为正整数，其他的参数一般也应适当圆整所以上最优解经标准化与圆整，最终解见表3.1。

表3.1 原设计及优化设计结果

V
传统
2
3.5
21
20
0.3
-0.3
0.4
-0.4
42
70
40
52
4.8
2296700
优化
2
3
21
19
0.3
-0.3
0.4
-0.4
42
70
40
52
4.8
1812500
优化之后体积和传统的方法体积相比较：（2296700-1812500）/2296700=21%。
3.7 本章小结
本章通过对带式运输机的二级齿轮减速器进行了具体的优化，在上章的理论基础上我们具体问题具体分析，通过带式运输机的功率转速传动比，以体积最小为目标建立起目标函数，得到优化条件。最后运用MATLAB强大的工具箱对数据进行优化得到了优化之后的数据参数。通过比较优化之后的体积比以前传统优化的结果要减小了21%.

4结 论
4.1 本文总结
通过对二级直齿圆柱齿轮的优化设计，我们可以得到以下结论：
（1）建立直齿圆柱齿轮优化数学模型：以体积最小为优化目标，以齿数，模数变为系数等等为优化变量，通过保证齿根弯曲疲劳强度和保证齿面接触疲劳强度建立约束条件。采用优化方法对二级减速器进行优化设计之后，使得减速器的结构更加紧凑，体积减小，重量减轻，设计优化后，并使设计参数达到最优值，这个对设计水平的提高。由于是对齿轮进行优化的同时，在建立齿轮的约束条件时，已经对轴和齿轮进行了校核，所以省去了再去校核计算。好对我们以后设计新产品提供了较好的设计数据，同时也让我们体会到设计的一步一步优化对产品质量和效率的重要性；由于优化设计一般多在完成初试设计之后进行，进而获得的优化结果，满足了齿轮传动的刚度、强度和使用寿命的要求。
（2）利用MATLAB中优化工具箱命对所建数学模型进行求解，应用MATLAB优化工具想求解优化问题，具有编程工作量小，初始参数输入简单，不需要编写大量算法程序，提高了设计效率，获得了良好的优化结果，提高了优化设计效率等优点，
可以得到最优值，通过对参数的标准值的选取和圆整得到了优化之后的具体体积。
（3）以某直齿圆柱减速器为例，运用所提出的方法对其进行优化，结果表明，在保证强度的条件下，体积减小了21%。从而验证了本文所提出的模型以及优化方法的可行性。
4.2 研究展望
在本文工作的基础上，作者认为还应进行以下方面的研究工作：
（1）在齿轮箱体积计算过程中可去掉其重合部分，这样更加合理。
（2）本文中仅以体积为优化目标，这在现代齿轮设计中是不全面的。因此可以同时以齿面接触应力和齿根弯曲应力最小、中心距最小等多目标进行二级直齿圆柱之论的优化设计。

参 考 文 献
[1] 减速器实用技术手册编委会. 减速器实用技术手册[M]. 北京: 机械工业出版社, 1992.
[2] 濮良贵, 纪名刚. 机械设计( 第七版)[M]. 北京: 高等教育出版社, 2001.
[3] 施阳，李俊.MATLAB语言工具箱—TOOLBOX实用指南. 西安: 西北工业大学出版社, 1998.
[4] 黄洪钟. 机械设计模糊优化原理与应用[M]. 北京: 科学出版社, 1997.
[5] 邱宣怀. 机械设计[M]. 北京: 高等教育出版社, 1997.
[6] 董玉革. 机械模糊可靠性设计[M]. 北京: 机械工业出版社,2001.
[7] 王健, 李露, 等. 模糊理论在机械设计中的应用[J]. 机械设计与制造, 2000( 4) : 5-7.
[8] 孙淑霞, 田芳, 等. 齿轮传动的可靠性优化设计[J]. 机械设计与制造, 2001( 5) : 3-6.
[9] 成经平. 机械产品运动方案创新设计与研究[J]. 机械设计与制造, 2001( 5 ) : 86-88.
[10] 叶元烈.机械优化理论与设计[M]. 北京:中国计量出版社,2001.
[11] 濮良贵, 纪名刚.机械设计[M]. 北京:高等教育出版社,2001.
[12] 李涛, 贺勇军, 刘志俭.M 工具箱应用指南——应用数学篇[M] .北京:电子工业出版社,2000.
[13] 刘惟信. 机械可靠性设计[M]. 清华大学出版社,2002.
[14] 徐锦康, 吕慧瑛, 刘极峰. 机械优化设计[M].北京: 机械工业出版社, 1995.
[15] Spiral and Worm Gearing , by the Machinery Publishing Co.,Sep.1914, pp.30and 31
[16] PATENT ABSTRACTS OF JAPAN vol. 1996, no. 03, 29 March 1996 (1996-03-29) -& JP 07 310807 A (JIDOSHA KIKI CO LTD

1. 基于C8051F单片机直流电动机反馈控制系统的设计与研究
2. 基于单片机的嵌入式Web服务器的研究
3. MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究
4. 基于模糊控制的电阻钎焊单片机温度控制系统的研制
5. 基于MCS-51系列单片机的通用控制模块的研究
6. 基于单片机实现的供暖系统最佳启停自校正（STR）调节器
7. 单片机控制的二级倒立摆系统的研究
8. 基于增强型51系列单片机的TCP/IP协议栈的实现
9. 基于单片机的蓄电池自动监测系统
10. 基于32位嵌入式单片机系统的图像采集与处理技术的研究
11. 基于单片机的作物营养诊断专家系统的研究
12. 基于单片机的交流伺服电机运动控制系统研究与开发
13. 基于单片机的泵管内壁硬度测试仪的研制
14. 基于单片机的自动找平控制系统研究
15. 基于C8051F040单片机的嵌入式系统开发
16. 基于单片机的液压动力系统状态监测仪开发
17. 模糊Smith智能控制方法的研究及其单片机实现
18. 一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制
19. 基于双单片机冲床数控系统的研究
20. 基于CYGNAL单片机的在线间歇式浊度仪的研制
21. 基于单片机的喷油泵试验台控制器的研制
22. 基于单片机的软起动器的研究和设计
23. 基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究
24. 基于单片机的机电产品控制系统开发
25. 基于PIC单片机的智能手机充电器
26. 基于单片