0 引 言

航空齿轮减速器通常工作于高转速高负荷工况，其结构复杂，是一个故障多发、可靠性差的部件,对其结构进行研究和优化具有重要的工程实际意义。传统的齿轮传动系统优化设计多未考虑齿轮所受载荷和转速等工况参数的随机性、工艺水平的差异性，以及材料强度等力学性能参数的波动对优化目标和约束条件的影响，因此，优化结果不可避免地与工程实际具有较大差异，即传动系统的性能仍可能随着噪声因素的变化而产生较大范围内的波动。工程实际中通常通过严格控制加工制造精度以尽可能消除影响因素的波动性，从而提高优化目标和约束条件的稳健性，但是实际中却难以实现，也会付出很大代价[1]。稳健优化设计理论[2]的蓬勃发展为齿轮减速器的结构优化设计提供了新思路。

1998年，Isaias C.Regalado[3]首次通过正交试验和田口法，以圆柱齿轮传动中大齿轮和小齿轮的弯曲疲劳寿命相差最小、接触疲劳寿命相差最小、噪声最小、效率最高和体积最小为优化目标，以齿数、压力角、螺旋角、齿宽、刀具变位系数和法面模数等为设计变量，同时加入噪声因素、压力角偏差、螺旋角偏差和扭矩偏差，进行稳健优化设计；R.Dong等[4]发现齿轮的安全系数对噪声因素比较敏感，采用插值分析方法对风电齿轮箱齿轮传动系统的强度进行稳健优化设计，并通过随机抽样对比了常规优化设计和稳健优化设计的结果，证明稳健优化设计的结果能够满足设计要求；Shaul Salomon等[5]采用主动稳健优化设计方法，考虑载荷等工况参数的变化，以减速器的能源效率最高和产品花费最低为目标，对减速器的级数和各级传动比进行优化；董恩国等[6]应用双响应面方法，建立了行星齿轮机构体积均值和方差的双响应面模型，对其进行稳健设计研究；周娜[7]将可靠性敏感性理论与稳健设计方法相结合，基于可靠性优化模型研究了齿轮传动的可靠性稳健设计问题，并提出了可靠性稳健设计的数值方法[7]。综上所述，齿轮传动系统稳健优化设计的研究大多只考虑优化目标的稳健性，而忽略约束条件的稳健性，并且多为单目标优化。

针对航空齿轮减速器优化设计中出现的诸多问题，本文以航空齿轮减速器中常用的斜齿圆柱齿轮传动系统为研究对象，综合考虑噪声因素的波动性，基于6σ稳健优化设计理论，以重量、效率、可靠度为优化目标，建立齿轮传动系统的多目标稳健优化设计数学模型，并采用NSGA-Ⅱ多目标优化算法进行优化求解，得到Pareto前沿解集。

1 6σ稳健优化设计基本理论

6σ稳健优化设计是从统计学的角度出发，用概率模型分析不确定因素对产品性能产生的影响，再依据概率分析的方法来控制随机变量对产品性能的影响，获得满足产品性能、可靠性和成本等各方面要求的最优解集。它是将6σ质量管理理念、可靠性设计、优化设计以及稳健设计相结合的一种现代设计方法，其中，对不确定性因素的处理全面综合了蒙特卡洛分析法、可靠性分析和基于可靠性的优化设计以及田口法等方法的优势，是一套完整的评估、改进可靠性和稳健性的优化设计方法[8]。在6σ稳健优化设计过程中，通过改进设计变量的均值和方差，有效地降低产品违反性能约束的概率，提高产品的可靠度。

6σ稳健优化设计的通用数学模型为

④保证必要的齿顶厚度

s.t. μgj(x,p)+mσgj(x,p)≤0 (j=1,2,…,n)

xL+mσx≤x≤xU-mσx

(1)

式中：μf(x)和σf(x)分别为优化目标的均值和方差；μgj和σgj分别为第j个约束条件的均值和方差；m为σ水平，本文m取值为6；x为设计变量；xL和xU分别为设计变量x的上限值和下限值；σx为设计变量的方差。

从式(1)可以看出：6σ稳健优化设计的目标函数中包含了设计目标的方差，从而使优化目标在接近最优解的同时，其波动量也达到了最小化，即保证了设计目标的稳健性；在约束中考虑了约束函数可能存在的容差，采用使约束函数偏离约束边界的方法，以此提高设计的可靠性，可靠性的度量是6σ水平。

2 航空齿轮减速器6σ稳健优化设计

以单级斜齿圆柱齿轮减速器为研究对象，通过分析并考虑可控因素和不可控因素，建立6σ稳健优化数学模型，并求取稳健解的。

2.1 可控因素和不可控因素

参考ISO 6336齿轮标准[9-11]，总结出影响齿轮强度、重量、效率等性能的相关参数，如图1所示，结合可控因素和噪声因素的概念，可以看出：可控因素包括结构参数和加工参数，其中，结构参数包含齿数、模数、压力角、螺旋角、变位系数和齿宽系数，加工参数包含刀具齿顶圆角半径和精度等级；不可控因素包括工况参数和材料参数，其中，工况参数包含功率和转速，材料参数包含接触疲劳极限应力和弯曲疲劳极限应力。

  

图1 齿轮性能相关参数

 

Fig.1 Parameters related to gear performance

可控因素和不可控因素中各随机变量的均值主要通过实验或者参数手册计算获得，而其标准差通过变异系数C获得[12]。

2011年版《全日制义务教育语文课程标准》指出：“语文课程是实践性课程，应着重培养学生的语文实践能力，而培养这种能力的主要途径也应是语文实践。”“应该让学生多读多写，日积月累，在大量的语文实践中体会、把握运用语文的规律。”可见，培养学生写作能力应通过大量实践去达到目标。

2.2 确定性优化设计模型

按照航空齿轮传动系统轻量化、高效率和高安全性的要求，在满足承载能力的条件下，把传动系统的重量和效率作为优化目标。由于齿轮弯曲疲劳失效和齿面接触疲劳破坏是航空齿轮最普遍的失效形式，而在航空发动机的维修中，失效的齿轮中有85%是由于齿面的接触疲劳破坏而失效的[13]。因此，将齿面接触强度也作为优化目标之一。

在多元参与态势中，基于分工合作基础上的政府自上而下的主动作为和企业、民间社会力量自下而上的多领域积极参与应成为其主要特点。[9] 一方面，政府必须在整个机制中主动作为，在针对公共消防安全服务方面的政策制定、对消防中介组织的引导以及相关法律的制度保障方面给予大力支持。 政府制定政策也要因势而变，根据每个地区的不同情况，在提供公共消防安全服务的过程中发现政策的不足之处，推动政策改革从而建立更完善的公共消防安全供给体系。 同时，政府要帮助社区引入资源，搭建起消防中介组织和社区的对接桥梁，在政府的推动下，消防中介组织的资源进入也会变得更加快速和有效，能够为以后的合作打下扎实的基础。

(1) 目标函数

设计目标有三个，分别为传动系统的重量最小、齿面接触强度最高、传递效率最高。

①传动系统的重量最小表示为

 

(2)

式中：ρ为齿轮的材料密度；V为传动系统的体积；b为齿宽；d1为主动轮的分度圆直径；i为传动比。

公办养老机构由于已经形成相当的规模和较为完善的管理制度，康复院、图书馆、理发室、健身房和室外活动场地一应俱全，社会信誉度高，入住率高，甚至申请已经排到了几年以后才能入住。

政府在社会主义市场经济条件下发挥的作用，主要包括经济调节、市场监管、社会管理、公共服务、环境保护等职能。水资源由于流动性、外部性、不确定性、自然垄断等原因，水资源配置过程中存在着广泛的“市场失灵”，单纯依靠自由市场无法导致有效的社会产出。因此，水管理中需要广泛的政府介入和政府干预。

②齿面接触强度最高表示为

Y2=-lnSH

(3)

③传递效率最高表示为

Y3=1-η

(4)

式中：η为传递效率。

可是，谁说“可”与“不可”都没有用，在词的演进过程中，整着整着就变化了。先是苏东坡，后是辛弃疾，把诗的言志的部分搞入了词的写作，打破了传统文人关于“诗庄词媚”的游戏规则。对此，有些诗词专家赞同，说是“一洗绮罗香泽之态”，不那样“娘”了，而有的诗词专家反对，认为词的表述就是“娘”的表述，所谓“诗化”是对词的离经叛道。清人朱彝尊、汪森编了一本《词综》，苏东坡的词只选了15首，周邦彦却选了37首，人们耳熟能详的岳飞的《满江红》根本没能入选，理由是他们的词太诗了。

(2) 设计变量

设计变量有六个，分别为主动轮齿数z1、法面模数mn、压力角α、螺旋角β、齿宽系数Фd以及精度等级Ag。

(3) 约束条件

约束条件有6个，分别为齿轮不发生根切，保证必要的重合度，不发生过渡曲线干涉，保证必要的齿顶厚度，强度约束以及设计变量的边界条件。

式中：gμ为函数g(x)在xi=μi的值。

①齿轮不发生根切

加工变位齿轮时，主、从动齿轮不产生根切的最小变位系数分别为

从图3可知，电压补偿器作为系统中用于接收和反馈信号的核心部件，位于辐射大厅，一端通过长约10 m的信号电缆与辐照晶体管连接，另一端则通过长约100 m屏蔽双绞电缆与测量间的测试系统相连，将VCC，VBB，VRB和VRC反馈至测试转接盒。为避免电阻直接受中子辐照引起性能下降，产生测量误差，将辐照板上的电阻RB，RC与晶体管分开，放置于电压补偿器中。

 

(5)

 

(6)

主动轮齿根与从动轮齿顶不发生干涉的条件为

②保证必要的重合度

重合度与齿轮传动的平稳性直接相关，一般对重合度的要求为

 

(7)

式中：εα为端面重合度；αat1和αat2分别为主动轮和从动轮的齿顶圆端面压力角；为端面啮合角。

③不发生过渡曲线干涉

式中：xn1、xn2分别为主、从动齿轮的法向变位系数；为法向齿顶高系数；z1、z2分别为主、从动齿轮齿数；αt为端面压力角；β为螺旋角。

 

(8)

从动轮齿根与主动轮齿顶不发生干涉的条件为

 

(9)

min μf(x)+mσf(x)

正变位系数增大会导致齿轮齿顶厚度减小，因此，变位齿轮的齿顶厚度应满足一定的约束条件。

主动轮齿顶厚度的约束为

 

(10)

从动轮齿顶厚度的约束为

 

(11)

⑥设计变量的边界条件

会计作为一门应用型学科，其内容从来不是一成不变的。全球经济形势的不断发展，要求会计、审计规则和方法与之相适应，其提供和反映的信息必须更可靠、更合理。人区别于机器和其他生物的最大特征之一在于人类具有优秀的学习实践和创新能力。而会计领域十分强调学习能力，继续教育种类和有关认证考试繁多。只有不断学习掌握高精尖的会计领域专业知识，具备探索和创新精神，才能在人工智能的浪潮中站稳脚跟、不被淘汰。

⑤强度约束

SF≥SFmin

(12)

SH≥SHmin

(13)

式中：SF和SH分别为齿轮的弯曲疲劳强度安全系数和接触疲劳强度安全系数；SFmin和SHmin分别为齿轮的最小弯曲疲劳强度安全系数和最小接触疲劳强度安全系数。

式中：da1和da2分别为主动轮和从动轮的齿顶圆直径；αn为分度圆的法向压力角；mt为端面模数。

 

(14)

2.3 6σ稳健优化设计数学模型

在上述齿轮传动系统确定性优化设计模型的基础上，基于6σ稳健优化设计理论，齿轮传动系统的稳健优化设计模型可表示为

采用人体能量检测仪对我校170名普通大学生进行测试发现，处于健康状况的有120人，占总人数的70.6%，处于亚健康的有50人，亚健康现患率为29.4%.其中：女生亚健康人数为34人，其亚健康现患率为35.8%；男生亚健康人数为16人，其亚健康现患率为21.3%，调查结果见表1.表1说明绍兴文理学院学生的亚健康状况普遍存在，女生亚健康状况要比男生严重.

min fk=μYk(x)+6σYk(x) (k=1,2,3)

数据经统计学软件(SPSS 18.0)分析，计量资料用“”表示，计数资料用“%”表示，经t/χ2检验。P<0.05提示差异有统计学意义。

s.t. μgj(x)+6σgj(x)≤0 (j=1,2,…,10)

xL+6σx≤xμ≤xU-6σx

(15)

6σ稳健优化设计方法必须计算输出的均值和方差，常用方法有蒙特卡洛模拟法和一次二阶矩法[14]，且二者实现起来并不十分复杂，但蒙特卡洛方法为随机抽样，其所需采样次数巨大，因此计算效率较低。本文采用一次二阶矩方法实现对期望和方差的求解。将函数g(x1,x2,…,xn)在其均值点处进行泰勒展开：

 

 

(16)

从图4可以看出，对于该数据集，使用CLIQUE算法聚类时，外部圆环部分缺失大量边界，圆心部分被聚类成两个不同的类，虽然GPCLIQUE算法的聚类质量有所提升，但效果也不理想，使用改进算法聚类时，成功找回外部圆环的边界和圆心，聚类效果明显。因此，使用本文提出的方法改进CLIQUE算法，能够有效的提高聚类结果的准确率。

设各xi相互独立，对式(16)取一次近似，可得：

对梨花女子大学的课程设置进行分析，其中通识必修课程36学分，专业基础课程22学分，专业深化课程60学分，专业必修课程9学分。在梨花女子大学教育工程学又分为不同的方向，每个方向学习的具体课程又有所不同。

 

=gμ=g(μi,μ2,…,μn)

(17)

 

 

(18)

或者写为

μg=gμ

(19)

 

(20)

把上述常规优化设计数学模型中的约束条件转化为以下形式：

g≤0

(21)

通过一次二阶矩法可以求得相应约束变量的标准差。

对于约束条件中σSF和σSH的求解，由于

结合σF、σFG、σH和σHG的变异系数公式：

 

(22)

 

(23)

式中：σHG和σFG分别为齿轮的许用接触应力和许用弯曲应力；σH和σF分别为齿轮的接触应力和弯曲应力。

对式(22)～式(23)等式两边取对数，可得：

lnSH=lnσHG-lnσH

(24)

lnSF=lnσFG-lnσF

(25)

由于σHG、σFG、σH和σF分别服从对数正态分布，结合一次二阶矩法，得：

 

(26)

 

(27)

（2）对于施工材料做不到严格的查验。施工材料是保证整个工程质量的基础前提，任何一个行业都不例外，尤其是建筑行业，对材料的要求标准更高。但在实际的施工过程中，对于施工材料的要求并没有严格的把控，材料监管员起不到严格检测的作用，基本上都是处于应付的状态，使建筑工程质量从最根本上失去了保障。

 

(28)

 

(29)

 

(30)

 

(31)

式中：C为变异系数，下标代表齿轮许用应力和应力计算过程的中间变量，其中，KA为使用系数；KV为动载系数；KHβ和KFβ分别为计算接触应力和弯曲应力的齿向载荷分布系数；KHα和KFα分别为计算接触应力和弯曲应力的齿间载荷分配系数；YF为齿形系数；YS为应力修正系数；KC为计算系数；ZE为弹性系数；YST为试验齿轮的应力修正系数；YNT为计算弯曲应力的寿命系数；YδrelT为相对齿根圆角敏感系数；YRrelT为相对齿根表面状况系数；YX和ZX为尺寸系数；ZN为计算接触应力的寿命系数；ZL为润滑系数；ZV为速度系数；ZR为粗糙度系数；ZW为齿面工作硬化系数。

就拿言情小说和武侠小说来说，一般读书是为了追逐故事情节而读，在阅读的过程很少思考，以浏览了解为主，可以说是浅阅读。但也有特别爱好者，会钻研里面的情节、语言及人物形象，甚至于还有人会自己琢磨它写法，自己也学习写小说。这些都是需要深入研读的。同样，如果以敷衍或者简单了解的态度，来阅读经典名著，可以称之为深阅读吗？在一个以故事情节为主线的阅读者心中，红楼梦和言情小说并没有什么太大区别。所以，依据阅读的客体来判断深浅阅读也是非常不合理的。

计算许用应力和应力所需要的中间变量的变异系数选取参见文献[15]。

由于服从对数正态分布的变量x的变异系数Cx和lnx的标准差之间存在如下关系：

σlnx=ln(+1)

(32)

根据Zhang Y M等[12]给出的中间变量的变异系数，可求得相应的σF、σFG、σH和σHG的变异系数，然后根据式(32)求出相应的对数标准差。

齿轮优化数学模型中的优化目标为重量、接触强度和效率。重量的标准差可由一次二阶矩方法求得；接触强度安全系数和效率的标准差则可以采用上述计算方法。

2.4 稳健优化设计模型的求解

齿轮传动系统稳健优化设计模型为多目标优化问题。处理多目标优化问题的传统方法为将各目标函数乘以权重系数并进行线性相加，从而将多目标优化问题转化为单目标优化问题。显然，在不同的权重系数下，模型求解所得到的结果也有所不同。随着优化理论的发展，基于Pareto最优解机制的多目标优化算法被提出，其基本思想为获取在无法改进任何目标函数的同时又不削弱至少一个其他目标函数的一组解集。本文采用K.Deb[16]于2002年提出的NSGA-Ⅱ多目标优化算法进行求解，该方法是目前非常优秀的进化多目标算法，已在工程优化中得到了广泛应用。

NSGA-Ⅱ优化算法的基本流程为：首先随机生成包含一定数量个体的初始种群，然后对初始种群进行非支配排序，再依次对种群进行选择、交叉和变异三种遗传算法的基本操作，从而得到下一代种群；其次，将上述操作得到的子代种群和父代种群合并再进行快速非支配排序，并计算各非支配层中个体的拥挤度；最后根据个体的拥挤度和非支配关系选择同等规模的个体生成父代种群并继续采用选择、交叉和变异三种遗传算法的基本操作产生子代种群。依此类推，直到达到所要计算的遗传代数为止。

3 算例分析

涡轮轴航空发动机某减速器由单级斜齿圆柱齿轮传动系统组成，减速比为3.54。输入功率为113.4 kW，输入转速为20 900 r/min，输出转速为5 915 r/min。

要建立齿轮系统的稳健优化设计数学模型，必须首先建立各个变量的方差或变异系数。根据随机变量的均值和方差的选取原则[12]，各变量的变异系数如表1所示。

 

表1 各变量的变异系数

 

Table 1 Coefficient of variation of variables

  

变量名称及符号变异系数C齿数z0法向模数mn0.002 5法向压力角αn0.002 5法向变位系数xn0.002 5刀具齿顶圆角半径系数ρ*fp0.1螺旋角β0.002 5齿宽b0.002 5接触疲劳极限应力σHlim0.08弯曲疲劳极限应力σFlim0.08功率P0.08转速n0.08

结合齿轮传动系统的6σ稳健优化设计数学模型，运用NSGA-Ⅱ多目标优化算法对稳健优化设计数学模型进行求解。采用NSGA-Ⅱ多目标优化算法求解齿轮传动系统的多目标稳健优化模型得到的Pareto最优解集如图2所示。

  

图2 NSGA-Ⅱ优化得到的Pareto最优解

 

Fig.2 Pareto optimal solution obtained by NSGA-Ⅱ

从优化得到的Pareto最优解集中找出一组和原始方案进行对比。齿轮传动系统优化设计前后设计参数的对比如表2所示。

 

表2 齿轮传动系统优化前后设计参数对比

 

Table 2 Comparison of parameters before and after optimization

  

设计参数原始方案优化方案主动轮齿数4141模数/mm21.5齿宽系数0.550.54精度等级65压力角/(°)2024螺旋角/(°)2119

齿轮传动系统的稳健设计结果和初始设计结果的对比如表3所示，包含齿轮的设计变量、重量、接触疲劳强度安全系数、效率和弯曲疲劳强度安全系数。

 

表3 齿轮传动系统优化前后性能参数对比

 

Table 3 Comparison of gear transmission system performance before and after optimization

  

性能参数原始方案优化方案质量/kg30.4+0.54-0.5412.3+0.15-0.15安全系数2.3+0.05-0.052.02+0.04-0.04传动效率99.37+0.18-0.1899.42+0.1-0.1

从表3可以看出：传动系统的整体重量与优化之前相比减少了47.4%；在工况波动和齿轮参数误差的影响下，效率相对提高了0.08%，效率的波动降低了44.4%；在保证弯曲疲劳强度安全系数和各个约束条件可靠性的前提下，接触疲劳强度安全系数降低了5.6%，波动量降低了20%。总体来说，稳健优化设计得到的设计结果，在稍微降低接触疲劳强度安全系数的情况下，极大地减少了重量并提高了效率，降低了在可控因素和不可控因素作用下的系统性能的波动，保证了系统的稳健性。

4 结 论

(1) 通过将6σ稳健设计基本理论和NSGA-Ⅱ多目标优化算法相结合，本文建立的齿轮传动系统多目标稳健设计优化模型，可以充分考虑航空齿轮减速器所特有的多种不确定性因素，满足了工程实际的需求。

(2) 以某直升机的涡轮轴发动机内部的单级斜齿圆柱齿轮传动减速器为例，发现基于6σ稳健优化设计方法在保证不违背约束条件的前提下，不仅优化了齿轮传动系统的性能，而且降低了在可控因素和不可控因素作用下传动系统性能的波动量。

参考文献

[1] R Dong, W Sun, H Xu. Robust design optimization of gear parameters for a wind turbine using interval analysis[J]. Journal of Mechanical Engineering Science, 2010(1): 2235-2245.

[2] Hans-Georg Beyer, Bernhard Sendhoff. Robust optimization-a comprehensive survey[J]. Computer Methods in Mechanics and Engineering, 2007, 196: 3190-3218.

[3] Isaias C Regalado. Application of robustness in the multi-objective optimization of gears[D]. Columbus: The Ohio State University, 1998.

[4] R Dong, W Sun, H Xu. Robust design optimization of gear parameters for a wind turbine using interval analysis[J]. Journal of Mechanical Engineering Science, 2010, 1: 2235-2245.

[5] Shaul Salomon, Gideon Avigad. Gearbox design for uncertain load requirements using active robust optimization[J]. Engineering Optimization, 2016, 48(4): 652-671.

[6] 董恩国, 张蕾, 孙奇涵. 基于双响应面法的行星齿轮机构稳健设计研究[J]. 机械传动, 2009, 33(2): 35-37.

Dong Enguo, Zhang Lei, Sun Qihan. Robust design of planet gear based on dual response surface methodology[J]. Journal of Mechanical Transmission, 2009, 33(2): 35-37.(in Chinese)

[7] 周娜. 多失效模式齿轮传动的可靠性稳健设计理论研究[D]. 沈阳: 东北大学, 2012.

Zhou Na. Theory research on reliability-based robust design of gear transmission with multiple failure modes[D]. Shenyang: Northeastern University, 2012.(in Chinese)

[8] 王丰超, 潘毅, 宋袁曾. 基于6σ的民机制造工艺稳健优化设计研究[J]. 工业工程与管理, 2014, 19(3): 127-134.

Wang Fengchao, Pan Yi, Song Yuanceng. Study on robust design and optimization technology of civil aircraft manufacturing process based on six-sigma[J]. Industrial Engineering and Management, 2014, 19(3): 127-134.(in Chinese)

[9] International Orgranization for Standards. ISO 6636-1-2006 Calculation of load capacity of spur and helical gears, Part 1: Basic principles, introduction and general influence factors[S]. Geneva, Switzerland, 2006.

[10] International Orgranization for Standards. ISO 6636-2-2006 Calculation of load capacity of spur and helical gears, Part 2: Calculation of surface durability(pitting)[S]. Geneva, Switzerland, 2006.

[11] International Orgranization for Standards. ISO 6636-3-2006 Calculation of load capacity of spur and helical gears, Part 3: Calculation of bending strength[S]. Geneva, Switzerland, 2006.

[12] Y M Zhang, Q L Liu, B C Wen. Practical reliability-based design of gear pairs[J]. Mechanism and Machine Theory, 2003, 38(12): 1363-1370.

[13] 刘滨春. 航空齿轮疲劳失效机理探究[J]. 中国科技信息, 2011(12): 108-109.

Liu Binchun. Research of three basic fatigue failure modes for the aero-engine transmission gear[J]. China Science and Technology Information, 2011(12): 108-109.(in Chinese)

[14] 李玉强, 崔振山, 陈军. 基于响应面模型的6σ稳健设计方法[J]. 上海交通大学学报, 2006, 40(2): 201-205.

Li Yuqiang, Cui Zhenshan, Chen Jun. Six sigma robust design methodology based on response surface model[J]. Journal of Shanghai Jiaotong University, 2006, 40(2): 201-205.(in Chinese)

[15] 黄洪钟. 机械传动可靠性理论与应用[M]. 北京: 中国科学技术出版社, 1995: 56-86.

Huang Hongzhong. Theory and application of mechanical transmission reliability[M]. Beijing: China Science and Technology Press, 1995: 56-86.(in Chinese)

[16] Deb K. An efficient constraint handling method for genetic algorithms[J]. Computer Methods in Mechanics and Engineering, 2000, 186(2): 311-338.