【3.1.17】基于非支配排序遗传算法的变速器多目标优化设计
发布时间:2024-10-18
发布时间:2024-10-18
第25卷第12期2008年12月
机械设计
JOURNALOFMACHINEDESIGN
V01.25Dec.
No.122008
基于非支配排序遗传算法的变速器多目标优化设计’
范钦满1’2,吴永海1’2,徐诚1
(1.南京理工大学机械工程学院。江苏南京210094;2.淮阴工学院,江苏淮安223003)
摘要:以某新型货车变速器为研究对象,基于多目标遗传算法和有限元法,建立其多目标优化模型,包括优化数学模型和参数化的有限元模型,以多学科优化软件iSIGHT为平台,集成有限元软件ANSYS,使用基于Pareto最优概念的非支配排序遗传算法,寻求变速器在最小体积、最大传动总重合度和最小中心距目标下的最优结构设计方案,并获得了问题的Pareto最优解集。文中提出的分析方法能够有效地解决复杂结构的多参数多目标优化问题。
关键词:变速器;优化设计;NSGA—II;ANSYS中图分类号:THll2
文献标识码:A
文章编号:1001—2354(2008)12—0062—04
解,文中引入一个重要的概念:Pareto最优解。
多目标优化问题一直是结构优化领域的难点和热点问题。由于多目标优化问题的各个目标之间大多相互联系、制约,甚至相互对立,每个目标具有不同的物理意义和量纲。这种竞争性和复杂性使得对其优化变得十分困难。自1896年Pareto
V
定义(Pareto最优解)坤J:设u∈舻是多目标优化模型的约束集以菇)∈矿是向量值目标函数。若戈’∈U,并且不存在其
他可行点菇,满足:
Z(x)≤五(*‘)
Vi=1,2,…,p
提出多目标优化问题以来,相关的理论和应用研究一直非常活跃。Koopmm强TC提出了Pareto最优解的概念。Arron提出了凸集有效点的概念。Zadeh
L
并且Z(菇)<Z(戈。)解。
了i=l,2,…,P
A从控制论的角度提出多目标控
则称名‘是多目标极小化模型的Pareto最优解,或称为非劣
制理论。JohnsenZ系统地提出了多目标优化问题的理论¨二J。由此发展出的算法有加权和法、8约束法、目标规划法等,这些算法的共同特点是将多目标问题转化为单目标处理,往往只能得到一个最优解,且不能保证Pareto最优。直到Holland在1975年提出遗传算法,由于其具有全局最优性、多方向性和隐含并行性等特点,使得带有潜在解的种群能够一代一代地维持下来。目前,遗传算法已成为解决多目标优化问题的得力工具,国内外学者对此进行了大量的研究p。1。基于Pareto最优概念的非支配排序遗传算法(Non-dominated
SortingGeneticAI-
多目标遗传算法和单目标遗传算法的区别主要体现在适应值分配、多样性保持和精英策略3方面。目前,多目标遗传算法已发展到第3代。NSGA—II算法是该领域最具代表性的算法之一。NSGA算法最早由SrinivasN和DebK提出,他们将非支配排序思想引入遗传算法,把对具有多个目标函数的计算转化为适应值计算,用于求解多目标问题。NSGA—II算法是以NSGA为基础发展起来的,通过采用快速非劣分类方法、拥挤距离选择算子和精英保留策略,克服了传统NSGA的计算复杂度高、非精英保存策略和共享参数难以确定的缺点,提高了算法的收敛性、鲁棒性和计算效率。NSGA—II的算法流程见图1。图中,只为第n代种群,p。为只产生的子代种群,两个种群联合在一起形成规模为2N的混合种群R。。使用非劣分类方法对种群R。进行分类。与只对子代种群Q。进行非劣分类相比,尽管计算量较大,但是它可以在整个子代和亲代进行全局
gorithm,NSGA一Ⅱ)是在基本遗传算法基础上发展起来的针对多目标问题的遗传算法,对于求解多目标问题十分有效。多档变速器是汽车传动系统的关键部件,结构复杂,设计目标较多,有关其优化设计领域一直为国内学者所关注并进行了相关研究,但研究方法大都局限于单目标或传统的多目标优化算法№.7J。文中以某新型货车变速器为研究对象,基于有限元分析方法和多目标遗传算法,进行变速器的多目标优化设计。建立了变速器的多目标优化模型,包括参数化的有限元模型和优化数学模型,以多学科优化软件iSIGHT为平台,集成著名有限元软件ANSYS,采用NSGA一Ⅱ算法,寻求变速器在最小体积、最大传动总重合度和最小中心距目标下的最优结构设计方案并获得了问题的Pareto最优解集。
性的非劣检验。新种群P川由不同非劣等级的个体填充。填
充过程从最优非劣等级开始,接着是次优非劣等级,依此类推。由于整个种群R。的规模为2N,新种群的Ⅳ个位置不能容纳所有的非劣解。当决定最劣等级的解时,该非劣等级可能存在比新种群剩下位置更多的个体,使用密度估计方法选择稀疏区域的个体填充该种群,保证了解的多样性。
1
多目标遗传算法与NSGA一Ⅱ算法简介
对于多目标优化问题,很难得到使所有目标都达到最优的
文中优化目标为变速器具有最小体积、最大传动齿轮重合度和最小中心距,属多目标优化问题,优化框架如图2所示。在优化框架中,NSGA一Ⅱ算法为编制优化程序的依据。包括优化模型中设计变鼍的编码、遗传操作的实现、种群分类、集成用于约束计算的有限元计算程序等。
收稿日期:2007—01—16:修订日期:2007一05一08
作者简介:范钦满(1966一),男.江苏靖江人.副教授,博士研究生,研究方向:现代设计理论与方法。CAD/CAE/CAM/PDM等,发表论文20余
篇。
万方数据
2008年12月
范钦满,等:基于非支配排序遗传算法的变速器多目标优化设计
63
墨学
初始化种群P_(*o)
计算目标函数值
二二1=
适值分配
l
遗传操作.产生子I
二二二[二代种群Q。l
二二二[二合并操作,产生混l合种群t=只uQl二二二[二
对瞄行非劣分类I十算拥挤距离并排序,l
P肿。=u:。.P,l
图1算法流程图
图2优化框架
2
变速器的优化数学模型
2.1目标函数的建立
文中的研究对象为某新型货车的变速器,3轴4档,使用斜
齿轮传动,如图3所示。减少变速器的体积有利于提高汽车动力性能、降低成本及节省材料。因此,文中以变速器体积最小作为第l目标函数。此外,选取传动齿轮传动总重合度、变速器中心距为第2、第3目标函数,优化框架如图2所示。图2中,体积、重合度和中心距计算使用数学模型,强度和刚度计算使用有限元软件ANSYS,建立的参数化有限元模型如图4所示。建立的优化目标函数如下所示:
图34档变速箱结构简图
目标函数:
minmize:F(x)=∽以石)7
式中斫——变速器体积函数。由于变速器的整体尺寸主要取决于其传
动轮系,故以变速器齿轮体积之和作为变速器体积,
石=}[至(蠢)2(钿确’岛】;
万
方数据五——轴中心距函数以=mnl(毛+龟)/2cos历;
五——传动齿轮重合度之和函数.由于优化目标是函数最小化的
形式.故取总重合度的相反数为目标函数,
五=一妻【(1.88-3.2(士+丢))c∞晟+警1;
毛——齿轮齿数。毛(i;1,2。…。8)5“_——模数,"k(f=1.2,3,4)5反——螺旋角,晟(i=l。2,3,4);bi——齿宽,虬(i=1,2,3,4)。
图4
变速箱齿轮有限元模型
2.2设计变量的选择
影响变速器设计的参数很多,文中取常啮合齿轮齿数、模数、螺旋角、齿宽、传动比屯(i=1,2,3)等主要参数为设计变量。
设计变量:
工=[it,如,岛,句。龟,m血,“吐,4d,““,卢l,尼,岛,风,
bl,62,b3,b4]’
式中五=争×争;
如墨ZI×詈;
i3=詈x詈;
%——齿轮齿数。t(‘=I,2.….8)。其中7.3一旄为中间变量,其值分别为:
m“c08岛刁i3(气+≈)
白2
ii西而i了万
知5瓦鬲丽瓦两
mdc08尾92(毛+z2)
≈。ii而丽i了万
m。lc08卢3刁i2(毛+龟)
毛2磊面丽瓦而
in甜cos岛龟(毛+z2)
mmIcos角卸iI(4+z2)
勿2
%2瓦面百瓦而
ii硒忑i7万
“nlCOS反Z2(q+免)
2,3约束的建立
2.3.1
传动比约束
汽车变速器一挡传动比‘受到汽车最大爬坡度、附着力及汽车最低稳定车速的约束。根据驱动轮与路面的附着条件及最大爬坡度要求,得到以下约束:
堂芒掣乱《拣
式中:a一——汽车最大爬坡度;
i0——主减速器传动比;
64
机械,.道路滚动阻力系数;
旷一地面附着系数;
,,T——传动系效率;G.——汽车总质量;只——驱动轮法向反作用力;r。——车轮滚动半径;t。——发动机最大转距。
变速器各挡传动比的比值对变速器的使用性能有很大的影响,比值过大会造成换挡困难,一般认为比值不宜大于1.7—
1.8
o如果比值过小,则挡位数增多,使得变速器结构复杂化。因
此,得到以下约束:
1.5《fl/屯≤1.8
1.5《i21i3≤1.8
2.3.2轴向力平衡约束
由于斜齿轮螺旋角较大,传动时会产生很大的轴向力。为尽量减小变速器箱体的受力,设计时应力求使变速器中间轴上的轴向力大致平衡。于是得到下列约束:
sin卢,m。224—8in历mnlZ2—20≤08in岛rttnIz2一sin卢lm正孔一20《0
sin
J93mnI乏一sin角m曲%一20≤0
sin卢l“D3气一sin岛rtznlz2—20≤0sin反mnIz2一sinplml..钿一20≤08in芦lm耐硇一sin风mnIz2—20≤0
2.3.3齿轮几何参数约束
模数约束,对于轻型货车,通常取2.25≤m。≤2.75。于是得到下列约束:
2.25《m“《2.75
(f=1,2,3,4)
齿宽约束,对于轻型货车,通常取6m。≤b≤8.5m。。于是得到下列约束:
6m碰《bf≤8.5m面
(i=1,2,3,4)
螺旋角约束,对于轻型货车,通常取18。≤p≤34。。于是得到下列约束:
18≤鹿≤34
(i=1,2,3。4)
最小齿数约束,中间轴一挡小齿轮齿数一般要求:
12≤缸≤17
2.3.4中心距约束
中心距约束是保证两啮合齿轮间正常传动的空间位置约束。选择的原则是在保证传递发动机最大转矩t一,变速器具有最大传动比和齿轮具有足够强度的条件下,应尽量减小中心距。对于轻型货车来说,中心距的选择范围受t。限制,于是得到下列约束:
14瓦≤掣≤17K
2.3.5
斜齿重迭系数约束
为了保证斜齿轮传动的平稳性,要求斜齿轮传动的轴向重
迭系数不小于l,于是得到约束条件:
wm耐一blsin觑≤0
(i=1,2,3,4)
2.3.6齿轮强度约束
根据国标GB3480w1997,齿面接触疲劳强度约束为:
zHzE%x√、/面b__Ft
11+咖i,一粕‰足鼬≤[盯]肿
(i=1,2,3,4)
万
方数据设计第25卷第12期
式中:,t——分度圆上的圆周力;
d一小齿轮分度圆直径;
ZH——节点区域系数;磊——材料弹性系数;
z胡——接触强度计算的重合度和螺旋角系数5K^——使用系数;粕——动载系数;
x啦——齿向载荷分配系数;x№——齿问载荷分配系数;[矿]帅——齿面接触疲劳强度极限。
根据国标GB3480--1997,齿根弯曲疲劳强度约束为:
矿Fi≤[盯]FP
(i=1.2,…,8)
式中:口Fi——齿轮i齿根应力,由内嵌的标准斜齿轮强度模块求解出,
如图4.该模块在有限元软件ANSYS中开发,在处理应力
约束时。考虑了齿轮的安全系数(取2.0);
[矿]阡——齿根疲劳强度许用应力。
假设齿面接触疲劳强度、齿根弯曲疲劳强度、齿面接触应力、齿根弯曲应力均服从正态分布,取可靠度不小于0.997,通过标准正态分布函数表查到对应的可靠性系数为2.75,得到可靠性约束条件为:
2.75一—In—d='H=PJ=-=I=n=O'HSl≤0
2.75一半坠墼≤0
 ̄/《HP+C2m ̄/谚FJ+c:聊
式中:,删,巳咿——分别为齿面接触疲劳强度极限应力均值和变差系
数;矿啦,乞髓——分别为齿面接触应力均值和变差系数;
矿H,C,订——分别为齿根弯曲疲劳强度极限应力均值和变差系
数;
盯脚,c羽——分别为齿根弯曲应力均值和变羞系数。3
变速器的多目标优化设计
采用文献[6]的算例进行计算:某3轴4挡机械式汽车变
速器,已知q=1
835
kg,t一=147
N m,Fz=10100N,rr=
313
mm,站=4.11,设计最高时速‰一120
km/h。齿轮材料
为渗碳合金钢20CrMnTi。
以多学科优化软件iSIGHT为平台进行多目标协同优化。iSIGHT软件主要具有自动化、集成化与最优化等方面的功能。通过iSIGHT软件的过程集成接口,可以方便地实现与有限元软件ANSYS的集成。需要在ANSYS中完成有限元模型的建立、加载和施加约束等工作,创建后缀为.txt的APDL命令流文件input.txt,将其作为iSIGHT调用ANSYS软件的输入文件。该文件记录了变速器齿轮的有限元模型和载荷信息,在iSlGHT中解析该输入文件,将齿轮的结构参数如模数、齿宽等作为解析参数,输入iSIGHT中的对应值进行参数替代。通过执行AN—SYS软件的批处理文件木.bat来生成输出文件output.txt,解析输出文件从中得到齿轮的最大应力和最大位移,作为校核齿轮强度的依据。每次ANSYS计算的结果,将与约束条件进行对比,判断解是否违反约束。
使用改进型非支配排序遗传算法(NSGA—II),对变速器优化模型进行多目标优化计算。设置NSGA一Ⅱ算法的控制参
2008年12月范钦满,等:基于非支配排序遗传算法的变速器多目标优化设计
数:种群规模为20,交叉概率为0.8,变异概率为0.05,最大进化代数为100代。经过100代运行,获得了该变速器的多目标优化Pareto解集。图5为优化运算结束时的种群分布图,为便于理解,图中使用相对坐标系,以初始方案的重合度、中心距和体积为基准值,具体数值见表l。从图5中可以看出,点最密集的区域围成一个立体,该立体的表面即为Pareto最优解曲面,立体区域中的点为Pareto支配解。Pareto最优解曲面具有以下特点:在整个可行解集中没有比Pareto最优解所对应的个体性能更好的解,各Pareto最优解之间也没有优劣之分。设计者可以.根据经验和对各目标的重视程度,从中选出优化设计最满意的解。可见,Pareto曲面对最优解的选取具有通用性,即在一种约束条件下选择某种方案,而当对优化目标的约束条件变化后仍能在该Pareto曲面围成的其他区域寻找符合要求的解,无须修改后重新进行运算。从图5可以看出,与初始设计方案相比。在整个Pareto最优解集上,变速器体积与中心距值均小于1,传动齿轮的重合度均大于1.1,变速器的体积目标、重合度目标和轴中心距目标均有较好的改善。经均衡,选择优化方案如表1所示。为了便于评价曲线上所对应各点的优化效果,表1同时列出了变速器原设计方案、文献[6]的优化方案和基于Pareto多目标优化的方案。表1为3种方案的优化结果比较。从表1可看出,与文献[6]相比,文中优化后的变速器体积有所增加,增加约0.25%,但齿轮传动总重合度大幅增加,增长约15.3%,同时中心距有所减小,减小约6.9%。另外,变速器第1档传动比略有增加,增加了变速器的变速范围。从理论上来说,变速器优化的所有Pareto最优解都集中于图5中的Pareto曲面上。文献[6]采用线性加权法进行多目标优化计算,它的优化结果对应于图5中的一点,选用不同的权系数对应于图中不同的点,且不能保证Pareto最优。
图5第100代种群分布表1优化性能比较变速器体积
总重合度
轴中心距
/mm3
/mm
第1档传动比
初始方案6198509.712770.70
3.52
文献[6]的优化方案562989lO.319l70.57
3.54文中的优化方案
564373
11.8980
65.70
3.79
结论
基于多目标遗传算法NSGA一Ⅱ和有限元分析法,以
iSIGHT为平台,获得了某变速器在最小体积、最大传动总重合
万
方数据度和最小中心距目标下的Pareto最优解集。优化结果形象直观,信息量大,可操作性强。文中所使用的研究方法在解决复杂结构多参数多目标优化问题上可行性较强。
参考文献
Santom。Esamuele.Globalmethodsin
multi—objective
optimization
andtheir印plication
to
a
mechanicaldesignproblem[J].Computers
in
Industry,1991,18(2):169—175.
【2]Tamaki.1-lisashi.Multi—objectiveoptimization
by
genetic出orithms:
aⅢ,iew[C].In:Proceedings
oftheIEEEConference
on
Evolution—
ary,ICEC96:517—522.
[3]
DebK,Pratap
A,AgnrwalS,eta1.A
fastandelitistmulti-objective
学eneticalgorithm:NSGA—n
EJ].IEEE
Transactions
on
Evolution-
ary
Computation,2002.6(2):182—197.
[4]KnowlesJD.ComeDW.Approximatingthenon-dominatedfront
U—
sing
theParetoarehivedevolution
stmte盯【J].EvolutionaryCompu—
tation,2000,8(2):149—172.
[5]
Zitzhr
E,LannumnsM,Thiele
L.Evolutionary
methodsfordesign
[J].OptimizationandControl,2002,28(1):25—33.
[6]蒋春明.阮米庆.汽车机械式变速器多目标可靠性优化设计[J].汽车工程,2007,29(12):l090一l093.
[7]蔡炳炎.机械式汽车变速器速比优化设计及扭转震动分析[D].
武汉:武汉理工大学,2005.
【8]
周明,孙树栋.遗传算法原理及应用[M].北京:国防工业出版
社,2002.
Multi-objectiveoptimization
design
of
speed
variator
based
on
genetic
algorithmofnon-dominatedordering
FAN
Qin.manl”,WU
Yo喈hall一,XU
Chen91
(I.School
ofMechanical
Engineering,Nanjing
Universityof
ScienceandTechnology,Nanjing210094,China;2.HuaiyinPoly
technicInstitute,Huai’821
223003,China)
Abstract:Takingthespeedvariatorofcertain
new
typedtruck
as
the
object
of
study,basing
on
themulti-objectivegeneticalso-
rithmandfiniteelementmethod
the
multi-objectiveoptimization
model
including
theoptimal
mathematicalmodelandtheparameter-
izedfiniteelement
model
was
established.Usingthemultidiseipli-
naryoptimizationsoftwareiSIGHTas
theplatform
to
integratethefa-
mousfinite
elementsoftwareANSYS.Utilizingthenon—dominated
orderinggenetic
algorithmbased
on
theParetooptimalconception
to
seek
out
theoptimalstructuraldesigningschemeofthespeedvaria-
totunderthetargetoftheleastvolume,thelargestoveralltrsD_smis-
sion
superposition
andtheminimal
center
distanceandthusoh-
mined
theParetooptimalsolution
set
oftheproblem.The
analytical
methodbeingputforwardinthis
papercouldresolveeffectivelythe
problemsofmulti parameterandmulti-objectiveoptimizationsfor
the
complex
structures.
Keywords:speedvariator;optimization
design;NSGA—II;
ANSYS
Fig5
Tab1Ref8
“JixieSheji”8294
4
基于非支配排序遗传算法的变速器多目标优化设计
作者:作者单位:
范钦满, 吴永海, 徐诚, FAN Qin-man, WU Yong-hai, XU Cheng
范钦满,吴永海,FAN Qin-man,WU Yong-hai(南京理工大学,机械工程学院,江苏,南京
,210094;淮阴工学院,江苏,淮安,223003), 徐诚,XU Cheng(南京理工大学,机械工程学院,江苏,南京,210094)机械设计
JOURNAL OF MACHINE DESIGN2008,25(12)2次
刊名:英文刊名:年,卷(期):被引用次数:
参考文献(8条)
1.Santoro;Esamuele Global methods in multi-objective optimization and their application to amechanical design problem 1991(02)
2.Tamaki,Hisashi Multi-objective optimization by genetic algorithms:a review
3.DebK;Pratap A;Agarwal S A fast and elitist multi-objective genetic algorithm:NSGA-Ⅱ 2002(02)4.Knowles J D;Come D W Approximating the non-dominated front u-sing the Pareto archived evolutionstrategy 2000(02)
5.Zitzhr E;Laumanns M;Thiele L Evolutionary methods for design 2002(01)
6.蒋春明;阮米庆 汽车机械式变速器多目标可靠性优化设计[期刊论文]-汽车工程 2007(12)7.蔡炳炎 机械式汽车变速器速比优化设计及扭转震动分析[学位论文] 20058.周明;孙树栋 遗传算法原理及应用 2002
本文读者也读过(4条)
1. 关志华 非支配排序遗传算法(NSGA)算子分析[期刊论文]-管理工程学报2004,18(1)2. 范钦满 强化创新意识 实践创新教育[期刊论文]-教育探索2003(5)
3. 颜伏伍.王洪建.田韶鹏.袁智军.Yan Fuwu.Wang Hongjian.Tian Shaopeng.Yuan Zhijun 基于第二代非支配排序遗传算法的变速器齿轮系多目标可靠性优化[期刊论文]-汽车工程2010,32(3)
4. 高俣.顾克秋.GAO Yu.GU Ke-qiu 陶瓷内衬多层复合身管结构多目标优化设计[期刊论文]-弹道学报2006,18(4)
引证文献(2条)
1.王晓云.原思聪.罗丹 灰色集成优化算法在齿轮变速箱设计中的应用[期刊论文]-广西大学学报(自然科学版)2011(2)
2.王晓云.原思聪.罗丹 基于混沌遗传算法的齿轮变速箱优化设计[期刊论文]-机械传动 2011(9)
本文链接:http://www.77cn.com.cn/Periodical_jxsj200812021.aspx