基于虛擬樣機(jī)聯(lián)合仿真的礦用汽車車架強(qiáng)度計算方法

金 純，孟慶勇，高小凡

(北京科技大學(xué)機(jī)械工程學(xué)院，北京 100083)

礦業(yè)是我國的基礎(chǔ)產(chǎn)業(yè)，礦用自卸汽車是礦石運(yùn)輸?shù)闹饕ぞ摺Ｓ捎诟咝浴⒔?jīng)濟(jì)性的需求，礦用汽車正向大型化和高速化發(fā)展。礦用汽車的發(fā)展趨勢對整車各個系統(tǒng)，尤其是車架的可靠性提出了更高的要求。礦用車輛行駛工況復(fù)雜，載荷譜采集困難，目前通常采用靜載估算的方法作為設(shè)計輸入，由于缺乏數(shù)據(jù)支撐，導(dǎo)致設(shè)計對象笨重且應(yīng)力集中現(xiàn)象嚴(yán)重，大大降低了礦用汽車的行駛經(jīng)濟(jì)性及可靠性。因此，針對礦用汽車的行駛工況，提出可行的強(qiáng)度計算方法，具有非常重要的意義。

載荷譜的測試及載荷譜的編制是可靠性計算的關(guān)鍵技術(shù)。載荷譜的獲取首先要對原始信號進(jìn)行降噪分析，通過雨流計數(shù)方法，得到研究對象壽命里程載荷的統(tǒng)計特性，在全壽命里程內(nèi)擴(kuò)展與重構(gòu)工作壽命內(nèi)的載荷譜，然后確定載荷幅值和均值的最大值。結(jié)合幅值和均值分級計算壽命里程內(nèi)各級的循環(huán)次數(shù)，對載荷譜進(jìn)行構(gòu)建。由于載荷譜和設(shè)計譜間匹配度較低，且費(fèi)用高、周期長，因此對于小批量生產(chǎn)的礦用汽車應(yīng)用較少。虛擬樣機(jī)技術(shù)融合了先進(jìn)建模方法和仿真技術(shù)，與傳統(tǒng)產(chǎn)品的設(shè)計方式相比，虛擬樣機(jī)技術(shù)主要強(qiáng)調(diào)系統(tǒng)整機(jī)性能及多領(lǐng)域協(xié)同仿真設(shè)計[1-2]。機(jī)械系統(tǒng)運(yùn)動學(xué)、動力學(xué)、有限元分析和控制理論是虛擬樣機(jī)技術(shù)的核心。通過虛擬樣機(jī)技術(shù)，可以避免產(chǎn)品設(shè)計缺陷，提前預(yù)測產(chǎn)品的運(yùn)行性能，從而獲得最佳的優(yōu)化設(shè)計方案[3-4]。

針對上述問題，本文通過建立虛擬樣機(jī)及其行駛環(huán)境的方法獲得礦用汽車的載荷譜，在此基礎(chǔ)上通過有限元的方法對車架進(jìn)行分析，以期獲得準(zhǔn)確的計算結(jié)果，滿足設(shè)計需求。

1 建立整車虛擬樣機(jī)

以載重量80 t的礦用汽車為對象，建立整車的三維模型，在此基礎(chǔ)上確定整車各部件的約束關(guān)系，構(gòu)建輪胎與路面模型，形成整車及運(yùn)行環(huán)境模型。

1.1 約束條件

在整車約束中，包括車架與貨箱、貨箱與舉升缸、貨物與貨箱、車架與前懸后懸等14個約束點(diǎn)，共有轉(zhuǎn)動、移動、球絞和固定四種運(yùn)動副，具體見表1。

表1 各部件的約束關(guān)系Table 1 Constraints of components

1.2 輪胎模型

通過輪胎模型，可由已知車輪運(yùn)動參數(shù)，得到輪胎沿X軸、Y軸和Z軸所受的力，描述輪胎在特定工作條件下的輸入輸出的關(guān)系。Magic Formula輪胎模型利用三角函數(shù)的組合公式擬合試驗輪胎數(shù)據(jù)。該模型通用性好，用一套形式相同的公式就能完整表達(dá)輪胎所受縱向力、側(cè)向力、回正力矩以及同時有側(cè)向和縱向運(yùn)動的聯(lián)合工況下的受力，形式簡單，計算方便，精度較高，且能最大限度地反映出車輛的實(shí)際運(yùn)行狀態(tài)，模型的表達(dá)式見式(1)。

Y(x)=

Dsin{Carctan[Bx-E(Bx-arctan(Bx))]}+Sv

(1)

式中：Y為輸出變量；x為輸入變量，一般包括縱向滑移率、輪胎側(cè)偏角等；系數(shù)B、C、D、E分別為剛度因子、形狀因子、峰值因子和曲率因子；Sv為曲線垂直方向漂移。

“魔術(shù)公式”輪胎模型的輸入輸出如圖1所示，輪胎模型的輸入為車輛的垂直載荷、滑轉(zhuǎn)率、側(cè)偏角及路面附著系數(shù)，輸出為輪胎的縱向力及側(cè)向力。

1.3 路面模型

路面模型采用3D組合路面，模擬礦用汽車的制動與越障工況，Adams環(huán)境中建立的路面模型如圖2所示。整車多體動力學(xué)虛擬樣機(jī)如圖3所示，考慮計算量的問題，對該虛擬樣機(jī)進(jìn)行簡化，共有70個運(yùn)動副。

圖1 輪胎輸入輸出模型Fig.1 Tire input and output model

圖2 Adams環(huán)境下的路面模型Fig.2 Pavement model in Adams environment

圖3 整車多體動力學(xué)虛擬樣機(jī)Fig.3 Virtual prototype of vehicle multi-body dynamics

2 確定安全系數(shù)

由于材料強(qiáng)度和應(yīng)力載荷不是確定值，而是服從一定分布規(guī)律的隨機(jī)變量，因此，目前采用材料強(qiáng)度和應(yīng)力載荷的均值對強(qiáng)度進(jìn)行計算來確定安全系數(shù)方法誤差較大。材料的極限應(yīng)力δlim是服從概率密度函數(shù)閥f(c)的隨機(jī)變量，而作用在零件的工作應(yīng)力s是服從于概率密度函數(shù)g(s)的隨機(jī)變量。通常情況下，材料極限應(yīng)力的平均值遠(yuǎn)大于工作應(yīng)力，但仍有可能工作應(yīng)力大于材料的極限應(yīng)力，即安全系數(shù)小于1的概率不為零。尤其是隨著材料疲勞強(qiáng)度降低，零件的承載能力將隨時間不斷衰減。依據(jù)可靠性設(shè)計理論，計算見式(2)和式(3)。

(2)

(3)

式中：R為可靠度；ZR為可靠性系數(shù)，計算見式(4)。

(4)

考慮礦用汽車的行駛工況及材料統(tǒng)計特性，通過上述計算，靜載下工況下安全系數(shù)n=1.32，動載工況下安全系數(shù)n=1.1。

3 仿真計算

基于分析研究的目的，本文選取了礦用汽車具有典型性的滿載制動與滿載過障兩個工況，計算車架的應(yīng)力水平，確定不同工況下的可靠性。

3.1 滿載制動工況

滿載制動工況下，由于整車質(zhì)量較大，為避免急加速導(dǎo)致慣性負(fù)載過大，整車從啟動到加速至額定最高車速47 km/h的時間設(shè)置為30 s，保持額定車速5 s后開始制動。歐洲經(jīng)濟(jì)委員會(ECE)和歐洲經(jīng)濟(jì)共同體(EEC)法規(guī)中規(guī)定總質(zhì)量大于12 t的貨車制動距離為ST=0.1v+v2/115，計算得到該工況下車輛的制動距離為23.9 m。

滿載制動時，車架主要承受來自于各縱向拉桿對車架的縱向載荷，多體動力學(xué)仿真結(jié)果見圖4。

圖4 滿載制動工況下作用于車架的縱向力Fig.4 Longitudinal force acting on the frame

車架的等效應(yīng)力計算結(jié)果如圖5和圖6所示。由圖5可知，滿載制動時，縱拉桿給予車架較大的作用力，最大應(yīng)力出現(xiàn)在車架拉桿支架處(圖6(c))，最大等效應(yīng)力為396.59 MPa。

圖5 滿載制動車架應(yīng)力圖Fig.5 Stress diagram of full load brake frame

圖6 滿載制動車架等效應(yīng)力局部放大圖Fig.6 Partial enlarged drawing of equal effect force of frame

3.2 滿載前輪越障

滿載前輪越障工況下，整車從啟動到加速至車速20 km/h的時間設(shè)置為5 s?？紤]到金屬礦山凹坑深度一般不大于0.1 m，因此該工況下設(shè)置凹坑深度為0.1 m。車架縱向載荷主要包括各縱向拉桿對車架的縱向力以及前懸處縱向力，多體動力學(xué)分析結(jié)果如圖7所示。

車架的等效應(yīng)力計算結(jié)果如圖8和圖9所示。由圖8可知，滿載制動時，縱拉桿給予車架較大的作用力，最大應(yīng)力出現(xiàn)在車架拉桿支架處(圖9(c))，最大等效應(yīng)力為309.85 MPa。

滿載越障與滿載制動工況是車輛行駛時最惡劣的兩個工況，由表2可知，在此兩種工況下，采用Q460C材料的車架結(jié)構(gòu)能夠滿足強(qiáng)度的要求。

圖7 滿載越障工況下作用于車架的縱向力Fig.7 Longitudinal force acting on the frame under full load obstacle crossing condition

圖8 滿載前輪過坑車架應(yīng)力圖Fig.8 Stress diagram of front wheel passing through pit with full load

圖9 滿載前輪越障車架等效應(yīng)力局部放大圖Fig.9 Partial enlarged drawing of full load front wheel obstacle crossing frame

表2 滿載上坡、越障工況下車架強(qiáng)度分析表Table 2 Frame strength analysis under full load uphill and obstacle crossing conditions

4 結(jié) 論

1) 在滿載越障與滿載制動等惡劣工況下，礦用汽車車架的危險截面發(fā)生在懸架與車架的連接處，因此，在車架的設(shè)計中，此處結(jié)構(gòu)尤其要避免應(yīng)力集中。

2) 在滿載制動工況和前輪滿載越障工況下，作用于車架的縱向力的變化趨勢和幅值符合其在實(shí)際工況下的變化，說明采用虛擬樣機(jī)技術(shù)能夠反映車輛的實(shí)際運(yùn)行工況。

3) 相比于載荷估算方法，基于虛擬樣機(jī)的車架強(qiáng)度計算方法可以準(zhǔn)確計算出車架在不同工況下所受的縱向力，從而有助于獲得可信度更高的強(qiáng)度計算結(jié)果。