徐中明，楊建國，張志飛，李仕生，李曉靈

(1.重慶大學 機械傳動國家重點實驗室，重慶 400030;2.重慶大學 機械工程學院，重慶 400030;3.重慶建設(shè)摩托車股份有限公司，重慶 400054)

摩托車懸架系統(tǒng)是整車的重要組成部分，對整車的行駛平順性、操縱穩(wěn)定性、騎乘的舒適性等有著很大影響［1］。用戶對摩托車綜合性能要求的不斷提高，國內(nèi)外對摩托車懸架系統(tǒng)的研究越來越受到重視。當前國內(nèi)外對摩托車懸架系統(tǒng)的研究較多的集中在懸架系統(tǒng)零部件性能改善［2－3］、單一工況下的性能提升等方面，從人－車系統(tǒng)和整車綜合性能角度對懸架系統(tǒng)進行設(shè)計研究的還較少。在平順性方面，當前大部分研究都是以隨機路面工況下降低垂向振動為目標［4－5］，并沒有考慮俯仰振動的影響，由于摩托車的懸掛質(zhì)量分配系數(shù)ε遠小于1，前后懸架系統(tǒng)具有高度的耦合性，因此，車身局部垂向振動的改善可能導致俯仰振動的加劇，從而使整車舒適性變差;此外，隨機路面工況下平順性的改善可能導致其它工況下舒適性的降低，由于摩托車前懸架剛度較小，加速制動時俯仰角太大會導致前叉管撞擊限位，從而影響整車舒適性。半主動懸架能夠適應(yīng)復(fù)雜多變的行駛工況，對于改善摩托車的垂向和俯仰振動能夠取得良好的效果［6－8］，并且能夠保證其它工況下的整車舒適性，但由于成本昂貴，在摩托車上的應(yīng)用受到限制，目前在國內(nèi)的摩托車上未見使用。俯仰和垂向振動是影響車輛乘坐舒適性的關(guān)鍵問題［9］，國內(nèi)外從考慮俯仰振動角度對摩托車平順性研究的文獻很少，而從多工況角度對摩托車懸架系統(tǒng)進行匹配研究的文獻還未有報道。

本文針對某125摩托車存在的制動俯仰角過大的問題對懸架系統(tǒng)進行優(yōu)化匹配，通過懸架特性參數(shù)的合理匹配實現(xiàn)制動俯仰角的降低，避免發(fā)生前懸架撞擊限位提高騎乘舒適性，同時保證隨機路面工況下的平順性能。采用摩托車動力學專用軟件BikeSim建立包含駕駛員在內(nèi)的人－車系統(tǒng)仿真模型，以加速/制動工況下的車身俯仰角最大值和B級路面工況下的俯仰、垂向振動加速度均方根值為優(yōu)化目標，在iSIGHT中建立兩種工況的并行優(yōu)化任務(wù)，集成BikeSim、MATLAB進行多目標優(yōu)化，取得了明顯效果。

1 加速/制動與平順性仿真

1.1 BikeSim中人－車系統(tǒng)模型的建立

整車模型的建立。BikeSim是一種基于數(shù)學模型的面向特性的參數(shù)化建模的動力學軟件，該模型共包含10個剛性體29個自由度，將整車分為車體、轉(zhuǎn)向、輪胎、懸架、制動、傳動等子系統(tǒng)，通過輸入各子系統(tǒng)的特性參數(shù)或曲線進行建模，其整車模型如圖1所示。

圖1 摩托車整車模型Fig.1 Vehicle model of motorcycle

由于摩托車后懸架結(jié)構(gòu)形式多樣，BikeSim將后懸架統(tǒng)一簡化為搖臂機構(gòu)，通過搖臂杠桿比來表達不同的懸架結(jié)構(gòu)。對于制動系，BikeSim用杠桿比來表達制動操縱機構(gòu)，將操縱機構(gòu)對制動系輸入的制動力簡化為制動杠桿力，制動杠桿力經(jīng)過增壓器、主缸、比例閥、輪缸等一系列部件的傳遞作用最終形成車輪上的制動力矩，由于本文只是研究相同制動力條件下不同懸架特性參數(shù)的車身俯仰角，只需制動系產(chǎn)生相同的減速度，因此制動系可采用BikeSim自帶的制動系模型。BikeSim提供的輪胎模型包括:“魔術(shù)公式”、查表輪胎模型和第三方輪胎模型，本文采用的是“魔術(shù)公式”。整車模型中制動系、傳動系和輪胎對本文的研究不產(chǎn)生影響，考慮到單獨建模的復(fù)雜性，這里采用了BikeS-im自帶的模型，其他子系統(tǒng)均是根據(jù)對象車的實際結(jié)構(gòu)提取的參數(shù)建立。建立的某125摩托車部分參數(shù)如表1所示。

表1 某125摩托車部分參數(shù)Tab.1 Some parameters of a 125 motorcycle

駕駛員模型的建立。由于摩托車駕駛員的質(zhì)量與摩托車整備質(zhì)量相當，因此在整車建模時必須考慮駕駛員的影響。當前大部分研究中對于駕駛員模型過于簡化，并不能全面的反映整車動力學特性;有的甚至忽略了駕駛員的影響，不符合真實情況。BikeSim中將駕駛員質(zhì)量分為上下半身兩部分，考慮了上下半身的質(zhì)量、轉(zhuǎn)動慣量、尺寸、姿勢等關(guān)鍵因素，上下半身通過轉(zhuǎn)動副連接，上半身可以相對于車架做側(cè)傾和側(cè)向運動。按照NASA的建議駕駛員上半身質(zhì)量占62%，下半身質(zhì)量占38%進行質(zhì)量分配。這里取駕駛員質(zhì)量為70 kg，建立的駕駛員模型如圖2所示。

圖2 駕駛員模型Fig.2 Driver model

1.2 路面模型的建立

加速制動仿真路面的建立。在BikeSim的道路建模部分根據(jù)GB20073的要求建立平直路面，設(shè)置道路附著系數(shù)為0.85，建立的平直路面如圖3所示。

圖3 平直路面Fig.3 Flat straight road model

平順性仿真路面的建立。摩托車是單軌車輛，平順性仿真時只需在BikeSim中建立道路高程隨路長變化的二維路面模型。按照GB/T7031對于道路分級的要求以及道路平度的表示方法，本文采用諧波疊加法［10－11］對 B級路面的平度進行模擬。

式中:θ為［0，2π］上的隨機數(shù);x為道路的x方向。通過離散道路x方向的值，就可以得到隨機路面z方向的值。

根據(jù)以上原理，采用MATLAB編制B級路面程序，得到的B級路面的平度數(shù)據(jù)導入BikeSim中，如圖4。

圖4 B級路面不平度Fig.4 Grade B road stochastic excitation

1.3 加速/制動工況的設(shè)定與仿真計算

圖5 速度變化曲線Fig.5 Target speed curve

在平直路面上進行加速/制動工況的仿真，加速工況設(shè)定為0～5 s速度從0均勻增加到60 km/h，之后進入勻速行駛;從第15 s開始制動，為了獲得更好的制動效果，采用聯(lián)動制動的方式給前后輪同時添加制動力。作用在摩托車上的制動力的大小與眾多因素有關(guān)，包括操縱機構(gòu)的杠桿比、液壓傳動系統(tǒng)、地面附著情況等，但制動力與制動系統(tǒng)并非本文研究重點，這里只是為了在附著良好的情況下獲得滿足要求的制動加速度，因此只要添加合適的制動力即可。在BikeSim中需要添加的是制動杠桿力(即制動踏板力與杠桿比的乘積)，不需要考慮復(fù)雜的制動操縱機構(gòu)，這里給前后輪添加的制動杠桿力分別為120 N和80 N。速度和制動力曲線如圖5、圖6。

設(shè)定仿真時間為25 s，運行仿真，得到的速度、俯仰角和前后懸架變形曲線如下圖所示。雖然制動時俯仰角加速度也會對行駛舒適性產(chǎn)生一定影響，但是由于俯仰角加速度主要受制動力大小的影響，懸架特性參數(shù)的改變對其降低效果并不明顯，因此，這里沒有給出俯仰角加速度的曲線。

從圖7可以發(fā)現(xiàn)，在制動力作用下17.5 s時速度降低為0，前進已經(jīng)停止，計算可知制動減速度大小為6.67 m/s2，完全滿足GB20073對于L3類摩托車制動性能的要求，說明添加的制動杠桿力大小是合適的。

從圖8可以看出加速階段俯仰角較小，俯仰角最大值出現(xiàn)在制動即將停止時，最大值為6.4°，制動階段的車身俯仰運動近似于整車繞后軸的旋轉(zhuǎn)，會導致前叉管發(fā)生壓縮運動，由圖9可以看出前叉管壓縮最大值為199.7 mm，而前叉管的行程僅為105 mm，因此必然會發(fā)生撞擊限位而影響整車的騎乘舒適性。

圖6 前后制動杠桿力曲線Fig.6 The front and rear break lever forces

圖7 加速/制動時速度曲線Fig.7 Speed curve during acceleration-braking

圖8 加速/制動時俯仰角曲線Fig.8 Pitch angle curve during acceleration-braking

圖9 加速/制動時前后懸架變形Fig.9 Suspension deformation during acceleration-braking

圖10 車身俯仰振動角加速度功率譜密度Fig.10 PSD of body pitching vibration acceleration

1.4 平順性工況的設(shè)定與仿真計算

B級路面下的平順性仿真。設(shè)定摩托車以車速60 km/h等速行駛于B級路面上，仿真時間90s，將計算結(jié)果導入MATLAB進行后處理得到俯仰和垂向振動加速度均方根值分別為 5.12 rad/s2、1.69 m/s2，圖 10、圖11分別是車身俯仰、垂向振動加速度功率譜密度。

圖11 車身質(zhì)心垂向振動加速度功率譜密度Fig.11 PSD of sprung mass vertical vibration acceleration

2 多工況并行任務(wù)多目標優(yōu)化的實現(xiàn)

2.1 問題描述

為了實現(xiàn)該車在加速/制動時行駛舒適性的改善，避免前叉管撞擊限位，同時保證隨機路面的平順性，防止俯仰或垂向振動的惡化;以前、后懸架的剛度、阻尼為設(shè)計變量，以加速/制動時俯仰角最大值和隨機路面工況下車身質(zhì)心俯仰、垂向振動加速度均方根值為優(yōu)化目標，采用NSGA－Ⅱ算法進行多目標優(yōu)化。通過前后懸架的剛度、阻尼的優(yōu)化匹配實現(xiàn)該車在這兩種工況下綜合性能的提高。

該優(yōu)化問題的數(shù)學模型可描述為:式中:fyjd是加速/制動工況下俯仰角最大值，rmsp、rmsv分別是俯仰和垂向振動加速度均方根值，F(xiàn)ront_K、Front_C、Rear_K、Rear_C分別是前、后懸架的剛度、阻尼。

2.2 多工況并行任務(wù)優(yōu)化平臺的搭建

BikeSim具有開放的接口，可以與多種軟件連接實現(xiàn)聯(lián)合仿真優(yōu)化。本文的多目標優(yōu)化平臺的搭建采用專業(yè)的優(yōu)化軟件iSIGHT通過接口技術(shù)集成BikeSim來實現(xiàn)，為了處理仿真計算過程中產(chǎn)生的大量數(shù)據(jù)，這里采用MATLAB編寫數(shù)據(jù)處理程序，并集成到iSIGHT中實現(xiàn)數(shù)據(jù)的自動化處理。在這三個軟件的集成中，BikeSim負責不同工況的仿真，MATLAB負責數(shù)據(jù)處理，iSIGHT的優(yōu)化控件提供優(yōu)化算法并驅(qū)動BikeSim和MATLAB實現(xiàn)迭代計算。每一次計算過程中，優(yōu)化控件給出的一組懸架特性參數(shù)將同時傳遞給加速/制動工況和B級路面工況，兩種工況將同時開始仿真，形成并行計算任務(wù)。建立的并行任務(wù)優(yōu)化工作流程如圖12，并行任務(wù)的計算結(jié)果將同時傳遞給優(yōu)化控件。

圖12 并行任務(wù)優(yōu)化工作流程Fig.12 The workflow of parallel optimized tasks

iSIGHT驅(qū)動BikeSim和MATLAB進行聯(lián)合仿真優(yōu)化，能夠使各個軟件發(fā)揮專長使優(yōu)化計算的精度和效率都極大提高。

2.3 參數(shù)設(shè)置與計算

本文以摩托車前后懸架的剛度和阻尼作為設(shè)計變量。由于摩托車的懸掛質(zhì)量分配系數(shù)ε遠小于1，與轎車的差別較大，前后懸架系統(tǒng)耦合緊密，因此不宜采用偏頻法來確定前后懸架剛度的取值范圍。這里是根據(jù)懸架的靜撓度來確定剛度的上下限的，在BikeSim中通過在靜平衡時懸架變形量在一定范圍內(nèi)來確定剛度的上下限，阻尼的上下限按照與剛度初始值擴大相同的比例來確定。這里將表1中懸架原始特性參數(shù)值作為設(shè)計變量初始值，初始值和上下限如表2。

表2 設(shè)計變量參數(shù)表Tab.2 The list of design variables

在平順性工況中，由于懸架動撓度的增大會使撞擊限位的概率增大，使平順性變差，因此在優(yōu)化的同時必須對前后懸架的動撓度進行約束。一般要求為:

式中:fd1，fd2分別為前后懸的動撓度;［fd1］=50 mm，［fd2］=25 mm分別是前后懸的限位行程。

遺傳算法具有魯棒性、全局最優(yōu)性、高效并行性、不要求函數(shù)連續(xù)可導等特點，近年來在車輛懸架方面的應(yīng)用越來越多［12］，然而基本遺傳算法只適用于傳統(tǒng)的單目標優(yōu)化問題。對于本文的多目標優(yōu)化問題，采用iSIGHT提供的NSGA－Ⅱ算法來實現(xiàn)，設(shè)定種群數(shù)量50、遺傳代數(shù)50、交叉率0.9。

經(jīng)過2 500次計算，得到的計算結(jié)果如圖13所示，將計算結(jié)果向坐標平面投影可以得到圖14、圖15，從圖中可以看出計算結(jié)果已收斂，在兩個坐標面上Pareto前沿已經(jīng)形成，由多目標優(yōu)化的基礎(chǔ)理論可知Pareto前沿是各種權(quán)重下最優(yōu)值點的集合，從計算結(jié)果可以發(fā)現(xiàn)兩個Pareto前沿分別近似為直線段。

3 優(yōu)化計算結(jié)果分析

3.1 目標權(quán)重與最優(yōu)值的確定

多目標優(yōu)化的最優(yōu)值是針對某一權(quán)重下的最優(yōu)值，為了找出最優(yōu)值首先要確定目標權(quán)重。但是先決策后尋優(yōu)的模式無法準確把握目標權(quán)重，該模式在應(yīng)用中存在明顯的缺陷，本文采用的是先尋優(yōu)后決策的模式［13－14］。由多目標優(yōu)化的基礎(chǔ)理論可知:多目標優(yōu)化的最優(yōu)值通常存在于Pareto前沿上。因此本文根據(jù)計算得到的Pareto前沿的特點來確定權(quán)重關(guān)系。從圖14、圖15可以看出兩個圖中的Pareto前沿近似為直線，通過數(shù)據(jù)擬合可以得到Pareto前沿的斜率，Pareto前沿的斜率反映了目標之間的權(quán)重關(guān)系。

設(shè)加速/制動工況下俯仰角最大值fyjd的權(quán)重為P;等速B級路面工況下俯仰振動加速度均方根值rmsp的權(quán)重為M，垂向振動加速度均方根值rmsv的權(quán)重為N。在rmsp－fyjd平面上，通過數(shù)據(jù)擬合可知Pareto前沿的斜率為 －0.39，為了計算方便這里取 －0.4。同樣，在rmsv－fyjd平面上得到的 Pareto前沿斜率為－0.83，這里取 －0.8，即權(quán)重關(guān)系為:

設(shè)X、Y、Z分別代表俯仰振動、垂向振動和俯仰角最大值，由于本文的優(yōu)化目標是希望這三個值均能降低，即求函數(shù)T的最小值，其表達式為:

將iSIGHT中2500次優(yōu)化計算的全部結(jié)果導入MATLAB中進行數(shù)據(jù)處理，第2195次計算結(jié)果為最優(yōu)值，從圖14、圖15可以看出2 195次與初始點的對比。

圖13 計算散點圖Fig.13 Optimization 3D scatter plot

圖14 rmsp-fyjd平面投影Fig.14 The rmsp-fyjd plane projection

將圖13的計算散點圖向兩個坐標平面投影形成的Pareto前沿是兩條近似直線，實際上在空間坐標系中Pareto前沿應(yīng)該是一個前沿面，在MATLAB中的數(shù)據(jù)處理即是找出平面T=0.4X+0.8Y+Z與Pareto前沿面的切點，前沿面上的切點即為最優(yōu)值點。

3.2 優(yōu)化前后懸架性能比較

優(yōu)化前后懸架特性參數(shù)如表3所示;優(yōu)化前后的動撓度均方根值對比如表4所示，優(yōu)化后的前后懸動撓度略有增大;優(yōu)化前后的制動俯仰角最大值fyjd、俯仰振動與垂向振動加速度均方根值rmsp、rmsv如表5所示，從表中可以看出，加速制動工況下的俯仰角最大值優(yōu)化效果最為明顯，降低幅度高達57.7%，同時平順性工況下的俯仰和垂向振動均有不同程度下降，完全實現(xiàn)了預(yù)期的優(yōu)化目標。

表3 優(yōu)化前后懸架特性參數(shù)對比Tab.3 Suspension parameters before and after optimization

圖15 rmsv-fyjd平面投影Fig.15 The rmsv-fyjd plane projection

表4優(yōu)化前后動撓度均方根值對比Tab.4 Suspension dynamic deflection RMS values before and after optimization

表5 優(yōu)化前后懸架性能對比Tab.5 Comparison of optimization objectives

圖16 優(yōu)化前后的俯仰角變化Fig.16 Pitch angle before and after optimization

圖17 優(yōu)化前后前懸架變形Fig.17 Front suspension deformation before and after optimization

圖18 優(yōu)化前后的后懸架變形Fig.18 Rear suspension deformation before and after optimization

圖19 優(yōu)化前后俯仰振動對比Fig.19 Comparison of pitching vibration

圖20 優(yōu)化前后車身質(zhì)心垂向振動對比Fig.20 Comparison of vertical vibration

圖16、17是優(yōu)化以后加速/制動工況下的車身俯仰角的變化以及前懸架變形關(guān)系圖，從圖中可以看出俯仰角和懸架變形量明顯降低了，優(yōu)化后前懸架變形最大值降低為88 mm，叉管行程為105 mm，因此前懸架撞擊限位的現(xiàn)象消除，騎乘舒適性得到了提高。圖18是優(yōu)化前后的后懸架變形圖，從圖中可以看出優(yōu)化后的后懸架變形比優(yōu)化前有相應(yīng)的增大，優(yōu)化后的懸架變形在前5 s的加速階段和勻速階段無明顯突變，最大值41 mm出現(xiàn)在制動停止后的回彈階段，小于后懸架的限位行程50 mm，因此在加速階段和停止后的回彈階段都不會引起后懸架的撞擊。圖19、20是優(yōu)化前后俯仰振動和垂向振動的功率譜對比，從圖中可以看出采用匹配后的懸架參數(shù)整車的俯仰和垂向振動均有降低，在整個頻率范圍內(nèi)都沒有出現(xiàn)惡化的現(xiàn)象，并且在低頻范圍內(nèi)垂向振動降低的效果更加明顯，優(yōu)化后的整車性能滿足了預(yù)期的要求。

通過以上優(yōu)化前后兩種工況下的懸架各項性能的對比分析，加速/制動工況的俯仰角峰值明顯降低，騎乘舒適性得到改善，隨機路面工況下的平順性得到了提高，驗證了多工況并行任務(wù)優(yōu)化方法的有效性。

4 結(jié)論

(1)采用摩托車動力學專用軟件BikeSim建立了路面和人－車系統(tǒng)模型，進行了加速制動工況和平順性工況的仿真;以降低加速/制動工況下的俯仰角最大值、等速B級路面工況下的俯仰和垂向振動為目標，采用iSIGHT集成BikeSim和MATLAB的聯(lián)合優(yōu)化仿真的方法建立了兩種工況的并行優(yōu)化任務(wù)，對摩托車懸架系統(tǒng)進行了多目標優(yōu)化，得到了Pareto前沿，實現(xiàn)了預(yù)期的優(yōu)化目標。

(2)采用先優(yōu)化后決策的模式，根據(jù)優(yōu)化計算結(jié)果對Pareto前沿進行了分析，提出了一種確定目標權(quán)重的方法，找出了該權(quán)重下的最優(yōu)值;通過對優(yōu)化前后的俯仰角、前懸架變形、俯仰和垂向振動加速度均方根值、功率譜密度的結(jié)果分析表明，消除了制動時的前懸架撞擊現(xiàn)象，騎乘舒適性得到了提高，同時隨機路面平順性得到了改善，驗證了該方法的有效性，為車輛綜合性能的優(yōu)化提供了重要參考。

［1］李德寬.汽車工程手冊摩托車篇［M］.北京:人民交通出版社，2001.

［2］Cossalter V，Doria A，Pegoraro R，et al.On the non-linear behavior of motorcycle shock absorbers［J］.Journal of Automobile Engineering，2010，224(1):15 －27.

［3］Audenino A L，Belingardi G.Modelling the dynamic behavior ofa motorcycle damper［J］. JournalofAutomobile Engineering，1995，209(4):249 －262.

［4］徐中明，盧少波，文 俐，等.基于MATLAB的摩托車舒適性仿真分析［J］.計算機仿真，2006，23(7):249－252.XU Zhong-ming，LU Shao-bo，WEN Li，et al.Simulation of ride comfort for motorcycle based on MATLAB［J］.Computer Simulation，2006，23(7):249 －252.

［5］陳 松，蹇開林，黃澤好，等.摩托車平順性的仿真研究［J］.機械強度，2006，28(3):465 －469.

CHEN Song， JIAN Kai-Lin， HUANG Ze-hao， etal.Simulation and research of motorcycle's rider comfort［J］.Journal of Mechanical Strength，2006，28(3):465 －469.

［6］Spelta C，Savaresi S M，F(xiàn)abbri L.Experiment analysis of a motorcycle semi-active rear suspension［J］. Control Engineering Pratice，2010，18(11):1239 －1250.

［7］ Wu L，Zhang W J.Hierarchical modeling of semi-active control of a full motorcycle suspension with six degrees of freedoms［J］. International Journal of Automotive Technology，2010，11(1):27 －32.

［8］吳 龍，陳花玲，陳玲莉.摩托車半主動懸架分層預(yù)測控制及仿真［J］.系統(tǒng)仿真學報，2006，18(8):2239 －2243.

WU Long，CHEN Hua-ling，CHEN Ling-li.Hierarchical preview control and simulation of semi-cative motorcycle suspension［J］.Journal of System Simulation，2006，18(8):2239－2243.

［9］于學華，丁康.汽車俯仰與跳動復(fù)合運動的分析［J］.噪聲與振動控制，2005，4:26 －29.

YU Xue-hua，DING Kang.Analysis of the pitching and bouncing motion of the car［J］.Noise and Vibration Control，2005，4:26 －29.

［10］張永林.用諧波疊加法重構(gòu)隨機道路不平順高程的時域模型［J］.農(nóng)業(yè)工程學報，2003，19(6):32 －34.

ZHANG Yong-lin.Time domain model of road irregularities simulated using the harmony superposition method［J］.Transactions of The Chinese Society of Agricultural Engineering，2003，19(6):32 －34.

［11］張志飛，徐中明，賀巖松.基于多體動力學的摩托車平順性分析［J］.中國機械工程，2010，21(8):988 －992.

ZHANG Zhi-fei，XU Zhong-ming，HE Yan-song.Analysis of motorcycle riding comfort based on multi-body dynamics［J］.China Mechanical Engineering，2010，21(8):988 －992.

［12］奉銅明，鐘志華，閆曉磊，等.基于NSGA－Ⅱ算法的多連桿懸架多目標優(yōu)化［J］.汽車工程，2010，32(12):1063－1066.

FENG Tong-ming，ZHONG Zhi-hua，YAN Xiao-lei，et al.Multi-objective optimization for multi-link suspension based on NSGA-Algorithm［J］.Automotive Engineering，2010，32(12):1063－1066.

［13］王 濤.汽車懸架參數(shù)的多目標多標準決策優(yōu)化［J］.農(nóng)業(yè)機械學報，2009，40(4):27 －32.

WANG Tao.Multi-objective and multi-criteriadecision optimization ofautomobile suspension parameters［J］.Transactions of the Chinese Society for Agricultural Machinery，2009，40(4):27 －32.

［14］王 濤，陶 薇.考慮隨機因素的汽車懸架參數(shù)多目標穩(wěn)健優(yōu)化［J］.振動與沖擊，2009，28(11):146 －149.

WANG Tao，TAO Wei.Multi-objective robust optimization of automobile suspension parameters considering random factors［J］.Journal of Vibration and Shock，2009，28(11):146－149.