摘要:針對傳統(tǒng)優(yōu)化算法全局性較弱,容易陷入局部解,以及遺傳算法的局部速率較慢和局部搜索能力不足等問題,采用將遺傳算法和復(fù)合形算法相結(jié)合的遺傳復(fù)合形算法。在此基礎(chǔ)上對拖拉機彈性懸架結(jié)構(gòu)進行了優(yōu)化設(shè)計研究,并且與傳統(tǒng)優(yōu)化方法進行了比較。結(jié)果表明,對同樣初始條件的拖拉機彈性懸架進行結(jié)構(gòu)設(shè)計時,采用遺傳復(fù)合形算法比采用傳統(tǒng)算法可以使得拖拉機彈性懸架鋼板彈簧的體積減?。玻叮ィ虼诉z傳復(fù)合形算法在求解彈性懸架多變量、多約束優(yōu)化問題的時候是有效的和正確的。
關(guān)鍵詞:優(yōu)化設(shè)計;彈性懸架;遺傳復(fù)合形算法;復(fù)合形法;拖拉機
中圖分類號:S219;TH135 文獻標(biāo)識碼:A 文章編號:0439-8114(2012)17-3859-03
Optimized Design of Tractor Elastic Suspension Structure Based on Genetic
Compound Method
HUANG Wei-hua,LIANG Ai-qin
(School of Mechanical and Electronical Engineering,Henan Agricultural University,Zhengzhou 450002,China)
Abstract: According to the global property of traditional optimization algorithm was weak, easily fall into local solutions, as well as local rate slow and lack of local search of genetic algorithm. Tractor elastic suspension was optimizated designed by using genetic compound algorithm method which combined the genetic algorithm and compound method, and it was compared with the traditional optimization methods. The results showed that the size of the tractor elastic suspension leaf spring reduced 2.6% when using the genetic complex algorithm rather than traditional algorithm under the same initial conditions. The algorithm is correct and effective in solving multivariable and constraint optimization problem of elastic suspension.
Key words: optimizated design; elastic suspension; genetic complex algorithm method; complex method; tractor
拖拉機具有良好的平穩(wěn)性,國外一些主要拖拉機公司生產(chǎn)的拖拉機產(chǎn)品,采用彈性懸架的越來越多,裝用彈性懸架極大地提高了駕駛操縱的舒適性和在公路上進行運輸作業(yè)的安全性。理論研究和實踐證明,拖拉機前后彈性懸架是影響車輛行駛平順性的重要指標(biāo)之一,而且對于衰減車身的自由振動和抑制車身、車輪的共振,減少車身的垂直振動加速度和車輪的振幅有重要的影響。彈性懸架對于拖拉機行駛穩(wěn)定性有著重要的影響,可以使車輛具有某種程度的轉(zhuǎn)向不足性,減少加速制動時候的仰頭和點頭,轉(zhuǎn)彎時側(cè)傾小。
遺傳復(fù)合形算法是結(jié)合了傳統(tǒng)算法的局部較快,現(xiàn)代算法[1]的全局性較好的方法。遺傳算法具有全域搜索的性能,可以求得全域的最優(yōu)值,但是其局部搜索時效率較低,速度不高。復(fù)合形算法作為一種傳統(tǒng)的優(yōu)化方法,在局部搜索的時候具有比較高的效率,但是不具備全局性的特征,結(jié)合這兩種算法,取長補短,使用遺傳復(fù)合形算法[2]較好地克服了復(fù)合形算法易陷入局部解和遺傳優(yōu)化算法收斂于最優(yōu)解較慢的特點,既有復(fù)合形算法局部搜索快的優(yōu)點,又有遺傳算法的全局搜索好的特點[3],具有較好的全局性和收斂速度。以此為基礎(chǔ)對懸架系統(tǒng)的結(jié)構(gòu)參數(shù)進行優(yōu)化設(shè)計,及以懸架鋼板彈簧體積最小為目標(biāo)函數(shù),并取得了較好的結(jié)果。
1 遺傳復(fù)合形算法[4]
算法的基本思想:首先由遺傳算法產(chǎn)生初始群體,并由概率選擇出父代種群,再按照復(fù)合形算法對父代種群進行改造,找出種群中的最好點、次好點以及最差點,通過反射、擴張、壓縮或向最好點壓縮操作以產(chǎn)生一個較好的頂點替代最差的頂點,并傳回種群替代最差點[5],再對新產(chǎn)生的群體進行交配,變異生成下一代的群體,反復(fù)進行,直至取得最優(yōu)值。
算法的具體步驟:①隨機產(chǎn)生N個染色體的初始群體CH(k),令k=0;②若算法滿足終止準(zhǔn)則,則輸出優(yōu)化結(jié)果,否則執(zhí)行下一步;③對群體中每個染色體計算Chi(k),i=1,2,…,N,計算適應(yīng)函數(shù)fi=fitness(Chi(k));④進行概率計算確定進入交配池的父代;用復(fù)合形算法對父代中的最差點進行改造,產(chǎn)生新的父代群體NGH(k+1);⑤對新群體進行交配,變異形成BCH(k+1);⑥令k=k+1,Chi(k)=BCH(k+1),返回第②步。
2 拖拉機懸架結(jié)構(gòu)優(yōu)化設(shè)計[6,7]
2.1 目標(biāo)函數(shù)和設(shè)計變量
從節(jié)省材料和降低成本的角度考慮,以鋼板彈簧的體積最小為設(shè)計目標(biāo),目標(biāo)函數(shù)為
V(X)min=V1+V2+…+Vi+…+Vn
=■Vi=■bhiLi (1)
式中,■Vi為鋼板彈簧總成的體積;n為鋼板彈簧葉片數(shù)目;b為鋼板彈簧葉片寬度;hi為鋼板彈簧第i片葉片厚度;Li為鋼板彈簧第i片葉片伸直的長度。
從目標(biāo)函數(shù)可以看出,懸架優(yōu)化設(shè)計變量有4個,即
X=[b,n,Li,hi]T i=1,2,…,n
2.2 約束條件
1)厚度限制條件[8]。為了保證鋼板彈簧良好的性能和足夠的壽命,應(yīng)使葉片厚度具有可淬透性,葉片厚度應(yīng)限制在最大淬透范圍之內(nèi),即
[h]-h≥0 (2)
式中,[h]為簧片最大淬透厚度,h為葉片厚度,對60Si2MnAl最大淬透厚[h]=10 mm。
2)鋼板彈簧伸直長度限制條件。從鋼板彈簧縱向角剛度來看,希望加長彈簧長度來提高縱向角剛度,從而減少車輛的縱向角振動。另一面彈簧伸直長度又受到總體布置的限制,因此要求其長度限制在某一允許范圍內(nèi),即
[l]-l≥0
l-0.95[l]≥0 (3)
式中,[l]為根據(jù)總體布置要求所允許的最大伸直長度之半;l為板簧最大伸直長度之半。
3)卷耳應(yīng)力的約束條件。其工作應(yīng)力應(yīng)小于許用應(yīng)力,即
[σ]-σ≥0 (4)
式中,[σ]為卷耳許用應(yīng)力,一般取值為350 N/mm2。
4)固定端應(yīng)力和端部接觸應(yīng)力的約束條件。
[σy]-(σi)0≥0 (5)
[σy]-(σi)c≥0 (6)
式中,[σy]為許用應(yīng)力,取值600 N/mm2;(σi)0為鋼板彈簧第i片固定段應(yīng)力;(σi)c為端部接觸應(yīng)力。
5)板簧剛度約束條件。彈簧實際剛度C2對設(shè)計要求剛度C2′相對誤差應(yīng)在3%~5%,才能保證車輛平順性要求。
設(shè)計剛度C2′=G2/fc (7)
式中,G2為滿載時懸架承載(一副板簧);fc為鋼板彈簧的靜撓度。
板簧實際剛度G2=■ (8)
式中,L為板簧最大伸直長度;E為材料彈性模量,?。玻薄粒保埃?N/mm2;k為非工作長度系數(shù),?。埃担唬訛轵T馬螺栓中心距;δ為擾度系數(shù),取1.8;J0為板簧總慣性矩。
約束條件■≤ε (9)
式中,ε為剛度收斂精度,?。常ァ担?。
6)板簧葉片厚度與寬度比值約束條件。為了保證板簧具有足夠側(cè)向剛度,能承受一定的側(cè)向力,故得
10-b/h≥0 (10)
b/h-6≥0 (11)
7)板簧葉片數(shù)限制條件。為了減少板簧在變形過程中摩擦損失,故得
n-4≥0 (12)
10-h≥0 (13)
8)板簧比應(yīng)力限制條件[9]。板簧使用壽命很大程度上決定于板簧工作應(yīng)力幅度大小,故板簧單位變形的應(yīng)力,即比應(yīng)力應(yīng)控制在允許范圍內(nèi)。比應(yīng)力
σ=■≤[σ] (14)
式中,[σ]為許用比應(yīng)力,取6 N/mm;σ為板簧比應(yīng)力;δ為撓度系數(shù),?。保?;k為非工作長度系數(shù),?。埃?。
9)板簧變性最大時最大應(yīng)力約束條件[10]。鋼板彈簧在極限撓度(沖擊載荷下彈簧最大變形)即鋼板彈簧變形時,最大應(yīng)力應(yīng)小于最大許用應(yīng)力。
σmax=σ(fc+fd)≤[σmax] (15)
式中,fc為板簧靜撓度;fd為板簧動撓度;[σmax]為最大許用應(yīng)力,取1 200 N/mm2。
10)制動或驅(qū)動時板簧最大應(yīng)力限制條件。制動時前鋼板彈簧在它的后半段出現(xiàn)最大應(yīng)力
σ1max=■≤[σ1max] (16)
式中,m1為重量重新分配系數(shù),?。保础保?;[σ1max]為制動或驅(qū)動時最大許用應(yīng)力,?。?200 N/mm2;G3為一個車輪上滿載負(fù)荷;l1為彈簧最大伸直長度之半;φ為地面附著系數(shù),?。埃?;Y為卷耳中心距地面高度;W0為板簧的總截面矩。
驅(qū)動時后鋼板彈簧在前半段出現(xiàn)最大應(yīng)力
σ2max=■+■≤[σ1max] (17)
式中,m2為重量重新分配系數(shù),?。保薄保?。
11)板簧中部應(yīng)力限制條件。車輛通過不平路面,垂直力達最大值時
σ=■≤[σ] (18)
式中,k′為動載荷系數(shù),k′=1+fd/fc。
由上述可知,拖拉機鋼板彈簧的懸架系統(tǒng)體積的最優(yōu)化設(shè)計,是一個由18個不等式約束組成的多維非線性規(guī)劃問題。
3 優(yōu)化應(yīng)用實例
已知一副板簧的簧下質(zhì)量m=103 kg,滿載時簧上質(zhì)量M=850 kg,空載時簧上載荷G0=3 000 N,騎馬螺栓中心距S=130 mm,滿載時卷耳中心至地面高度Y=375 mm,板簧最大伸直長度之半的取值范圍650~680 mm,板簧寬度取值范圍60~80 mm,板簧葉片數(shù)取值范圍4~10,滿載時板簧的動撓度fd=60 mm。其優(yōu)化結(jié)果為:n=6,b=80 mm,h=6.5 mm,l=1 059 mm,其體積為1 653 158.5 mm3。采用常規(guī)的隨機方向法優(yōu)化,其優(yōu)化結(jié)果為:n=6,b=80 mm,h=6.7 mm,l=1 064 mm,其體積為1 696 582.9 mm3。這兩種優(yōu)化方法進行比較,前一種方法鋼板彈簧體積減小了2.6%。所以用遺傳復(fù)合形算法的效果是非常明顯的。
4 小結(jié)
采用的遺傳復(fù)合形算法結(jié)合了遺傳算法的全局搜索和復(fù)合形算法的局部搜索的特點,運用遺傳復(fù)合形算法進行拖拉機懸架鋼板彈簧體積的優(yōu)化,比普通優(yōu)化有著更佳的結(jié)果,是一種全局性較好、效率較高、求解有效的優(yōu)化算法,其具有重要的現(xiàn)實意義,可以進一步推廣應(yīng)用。
參考文獻:
[1] GEN M, CHENG R W. Gentic Algorithms & Engineering Optimization[M]. New York: John Wiley & Sons,2000.
[2] 黃偉華. 基于混合遺傳模擬退火算法的離合器蝶形壓緊彈簧優(yōu)化設(shè)計[J].機械強度,2009,31(4):675-677.
[3] 張智文.基于遺傳算法的機械傳動方案可拓優(yōu)化設(shè)計[J].機械制造與自動化,2011,40(2)22-24.
[4] 孫靖民.機械優(yōu)化設(shè)計[M].北京:機械工業(yè)出版社,2005.
[5] 張寶吉.基于自適應(yīng)小生境遺傳算法的船型優(yōu)化[J].計算機工程,2011,37(8):207-209.
[6] 程賢福,吳志強,張 鑫.車輛轉(zhuǎn)向機構(gòu)運動精度的可靠性穩(wěn)健優(yōu)化設(shè)計[J].機械設(shè)計與制造,2011(3):8-10.
[7] 高象平.拖拉機零部件優(yōu)化設(shè)計[M].廣州:廣東科技出版社,1995.
[8] 李 軍.基于MATLAB優(yōu)化工具箱的農(nóng)用運輸車轉(zhuǎn)向梯形優(yōu)化設(shè)計[J].農(nóng)機化研究,2011(5):147-150.
[9] 章一鳴,張錫清.車輛懸掛系統(tǒng)優(yōu)化設(shè)計[M].北京:北京理工大學(xué)出版社,2002.
[10] 商躍進,王 紅.鐵路貨車變剛度彈簧組疲勞強度設(shè)計方法研究[J].機械強度,2006,28(1):141-145.