張和平,徐文濤,唐運(yùn)軍,巫紹寧,趙木青
(1.武漢理工大學(xué)機(jī)電工程學(xué)院,湖北 武漢 430070;2.上汽通用五菱汽車股份有限公司技術(shù)中心,廣西 柳州 545007)
汽車驅(qū)動(dòng)橋殼在傳統(tǒng)設(shè)計(jì)中往往只考慮典型工況下的靜態(tài)結(jié)構(gòu)強(qiáng)度,然而汽車在各種道路上行駛時(shí)會(huì)受到復(fù)雜的交變載荷作用,當(dāng)設(shè)計(jì)不當(dāng)或者制造工藝有問(wèn)題時(shí),會(huì)引起橋殼疲勞,發(fā)生斷裂[1]。同時(shí)微型汽車橋殼為沖壓焊接件,焊縫的疲勞強(qiáng)度一般低于母材的強(qiáng)度,因此焊縫疲勞強(qiáng)度對(duì)橋殼的疲勞壽命有很大影響。有限元法作為一種虛擬仿真技術(shù)在汽車領(lǐng)域應(yīng)用越來(lái)越廣泛,特別是模擬分析對(duì)象在各種工況下的受力情況,可以在分析結(jié)果的基礎(chǔ)上進(jìn)行疲勞壽命計(jì)算,避免通過(guò)實(shí)車道路試驗(yàn)來(lái)獲取汽車的各項(xiàng)性能參數(shù)。將該方法與試驗(yàn)相結(jié)合,減少試驗(yàn)次數(shù),可大大縮短產(chǎn)品開發(fā)周期,減少研發(fā)成本。
筆者針對(duì)橋殼進(jìn)行有限元分析,在分析結(jié)果的基礎(chǔ)上進(jìn)行疲勞壽命計(jì)算,計(jì)算流程如圖1所示。首先對(duì)橋殼模擬實(shí)際工況進(jìn)行靜力學(xué)分析,得到其應(yīng)力分布,在此基礎(chǔ)上對(duì)載荷-時(shí)間歷程計(jì)數(shù),提取疲勞載荷譜,由于載荷是非對(duì)稱變幅載荷,應(yīng)用平均應(yīng)力修正法對(duì)其進(jìn)行修正,結(jié)合材料的S-N曲線計(jì)算單個(gè)循環(huán)的損傷。最后運(yùn)用Miner累積理論對(duì)循環(huán)損傷進(jìn)行計(jì)算,即可求得橋殼的疲勞壽命。
一般情況下,機(jī)械結(jié)構(gòu)疲勞失效的研究主要有兩個(gè)方向:①根據(jù)載荷譜確定加載程序,在試驗(yàn)臺(tái)架上進(jìn)行疲勞試驗(yàn)得到實(shí)際壽命;②根據(jù)零件結(jié)構(gòu)、材料特性,以及載荷譜用Miner累積理論估算疲勞壽命。以上兩種方法的關(guān)鍵部分是獲取準(zhǔn)確的載荷譜,結(jié)構(gòu)正常工況下測(cè)得的載荷-時(shí)間歷程稱為工作譜,由于隨機(jī)載荷的不確定性,需要對(duì)其進(jìn)行統(tǒng)計(jì)計(jì)數(shù)。目前針對(duì)隨機(jī)的載荷-時(shí)間歷程進(jìn)行計(jì)數(shù)的主要方法有雨流計(jì)數(shù)法,它的主要功能是將實(shí)測(cè)的載荷-時(shí)間歷程以離散載荷循環(huán)的形式表示出來(lái)[2]。
圖1 疲勞壽命計(jì)算流程圖
疲勞壽命預(yù)測(cè)的準(zhǔn)確程度既取決于載荷譜的精確統(tǒng)計(jì),更取決于損傷模型的合理選取。結(jié)構(gòu)零件經(jīng)常承受復(fù)雜的多軸比例與非比例交互循環(huán)載荷作用。早期經(jīng)常將比例載荷作用下的多軸疲勞問(wèn)題用單軸理論處理。隨著疲勞理論的發(fā)展,預(yù)測(cè)多軸疲勞壽命的臨界平面法得到應(yīng)用,該方法基于斷裂模型及裂紋萌生機(jī)理,認(rèn)為裂紋發(fā)生在某一特定平面上,疲勞損傷的累積和壽命預(yù)測(cè)都在該平面上進(jìn)行[3]。當(dāng)應(yīng)用臨界平面法找到最大損傷平面時(shí),得到其剪切和法向應(yīng)力,由于它們很復(fù)雜,故將其轉(zhuǎn)換為等效的應(yīng)力和應(yīng)變,然后應(yīng)用單軸疲勞理論來(lái)計(jì)算其疲勞壽命。
單軸疲勞壽命的出發(fā)點(diǎn)是應(yīng)變-壽命關(guān)系式:
或分開寫成:
式中:σ'f為疲勞強(qiáng)度系數(shù);ε'f為疲勞延續(xù)系數(shù);b為疲勞強(qiáng)度指數(shù);c為疲勞延續(xù)指數(shù);E為彈性模量;Δεe/2為彈性應(yīng)變分量;Δεp/2為塑性應(yīng)變分量;Δε/2為總應(yīng)變量。
由于該關(guān)系式是在對(duì)稱載荷作用下得到的,對(duì)于變幅載荷作用下的疲勞問(wèn)題,需要進(jìn)行平均應(yīng)力修正。
當(dāng)材料處于彈性范圍內(nèi)時(shí),平均應(yīng)力對(duì)疲勞壽命影響比較大,因此需要考慮在內(nèi)。而當(dāng)處于塑性范圍內(nèi)時(shí),由于平均應(yīng)力的松弛效應(yīng),影響可以忽略不計(jì),對(duì)式(2)進(jìn)行修正的修正公式為:
式中:σa為應(yīng)力幅;σm為應(yīng)力平均值;σr為等效應(yīng)力幅。
修正后的應(yīng)變-壽命關(guān)系式為:
應(yīng)用該損傷模型時(shí),首先將求得的載荷-時(shí)間歷程轉(zhuǎn)化為名義應(yīng)力-時(shí)間歷程,然后利用雨流計(jì)數(shù)法進(jìn)行循環(huán)計(jì)數(shù),計(jì)算每個(gè)循環(huán)的疲勞損傷,最后應(yīng)用Miner累積理論對(duì)損傷進(jìn)行累積計(jì)算,所得的結(jié)果即為某一點(diǎn)的疲勞壽命。
該橋殼為整體式橋殼,是沖壓焊接件,焊縫部位的材料與母材的力學(xué)性能不匹配,常常伴有應(yīng)力集中、焊接缺陷等,對(duì)橋殼整體的疲勞壽命影響很大,因此要特別考慮焊縫的疲勞強(qiáng)度。焊縫疲勞強(qiáng)度主要取決于焊縫的結(jié)構(gòu)、焊接工藝和焊接質(zhì)量。焊縫結(jié)構(gòu)的疲勞破壞是大多數(shù)事故發(fā)生的主要原因,為了對(duì)焊縫進(jìn)行相關(guān)的疲勞分析,國(guó)際焊接協(xié)會(huì)(IIW)和歐洲相關(guān)機(jī)構(gòu)分別制定了一些標(biāo)準(zhǔn),主要有美國(guó)AAR貨車結(jié)構(gòu)疲勞設(shè)計(jì)標(biāo)準(zhǔn)、英國(guó)BS7608標(biāo)準(zhǔn)和IIW標(biāo)準(zhǔn)等。筆者主要采用英國(guó)標(biāo)準(zhǔn)BS7608-1993,它根據(jù)焊接的接頭類型、承載方向等將焊縫的疲勞特性分成不同的等級(jí),處于相同等級(jí)的焊接結(jié)構(gòu)具有相同的疲勞特性[4]。每個(gè)等級(jí)都有相應(yīng)的 S-N曲線與之對(duì)應(yīng),如圖2所示。
圖2 焊縫S-N曲線
通過(guò)靜力學(xué)分析得到橋殼的應(yīng)力分布,進(jìn)而獲得疲勞載荷譜,因此靜力學(xué)分析的結(jié)果要盡量準(zhǔn)確。筆者運(yùn)用CAD軟件UG對(duì)橋殼進(jìn)行建模,在保證橋殼力學(xué)性能的前提下,同時(shí)又使建模簡(jiǎn)單,對(duì)橋殼進(jìn)行了簡(jiǎn)化,省略了一些對(duì)橋殼受力沒有影響的油孔和螺栓孔,同時(shí)用環(huán)形圈模擬橋包與半軸套管、半軸套管與半軸法蘭間的焊縫。在ANSYS Workbench中建立橋殼的有限元模型,對(duì)其采用自動(dòng)網(wǎng)格劃分,最終單元數(shù)為71718,節(jié)點(diǎn)數(shù)為147608,后橋殼整體材料為20號(hào)鋼,彈性模量為213 GPa,泊松比為0.282,許用應(yīng)力為245 MPa。有限元模型如圖3所示。
圖3 橋殼有限元模型
橋殼滿載后軸荷為1040 kg,根據(jù)QC/T533-1999《汽車驅(qū)動(dòng)橋臺(tái)架試驗(yàn)方法》中的規(guī)定,對(duì)橋殼兩端板簧處施加2.5倍的滿載軸荷,同時(shí)為了模擬橋殼實(shí)際工況又消除橋殼的剛性位移,在一端約束Y軸、Z軸平動(dòng)自由度和繞Y軸、Z軸的轉(zhuǎn)動(dòng)自由度,另一端約束X軸、Y軸、Z軸的平動(dòng)自由度和繞Y軸、Z軸的轉(zhuǎn)動(dòng)自由度[5]。
橋殼應(yīng)力分析界面圖如圖4所示,半軸套管與橋包焊縫處和板簧座處應(yīng)力較大,其中半軸法蘭變截面處的最大應(yīng)力為220.23 MPa,小于材料許用應(yīng)力245 MPa。同時(shí)最大變形量出現(xiàn)在橋包處,為 0.636 mm,橋殼每米輪距變形為 0.463 mm,遠(yuǎn)小于QT/C 534-1999《汽車驅(qū)動(dòng)橋臺(tái)架試驗(yàn)評(píng)價(jià)指標(biāo)》所規(guī)定的1.5 mm,橋殼滿足強(qiáng)度、剛度要求。
圖4 橋殼應(yīng)力分布界面圖
一般情況下,根據(jù)疲勞壽命分析結(jié)構(gòu)應(yīng)力分布,首先得到單位載荷響應(yīng),根據(jù)實(shí)測(cè)載荷譜來(lái)按比例獲得整個(gè)結(jié)構(gòu)的應(yīng)力和應(yīng)變-時(shí)間歷程[6]。由于橋殼是恒幅加載,可以直接提取其工作載荷的結(jié)果響應(yīng),因此直接在橋殼2.5倍滿載軸荷下應(yīng)力分布的基礎(chǔ)上計(jì)算疲勞壽命。同時(shí),由于平均應(yīng)力對(duì)疲勞極限影響較大,對(duì)應(yīng)力-應(yīng)變時(shí)間歷程進(jìn)行平均應(yīng)力修正,獲得名義應(yīng)力-應(yīng)變-時(shí)間歷程。
橋殼非焊接部位材料為20號(hào)鋼,根據(jù)材料手冊(cè),得到其不同應(yīng)力所對(duì)應(yīng)的循環(huán)次數(shù),擬合成的S-N曲線如圖5所示[7],焊縫的S-N曲線如圖2所示的G等級(jí),將兩種S-N曲線賦予有限元模型中。
圖5 橋殼材料的S-N曲線
當(dāng)施加2.5倍滿載軸荷時(shí),得到的疲勞壽命分布界面圖如圖6所示,從圖6中可以看出,橋包上與半軸套管結(jié)合處是薄弱位置,其最低疲勞壽命為829300次,達(dá)到80萬(wàn)次的標(biāo)準(zhǔn)要求。
圖6 橋殼疲勞壽命分布界面圖
為了驗(yàn)證橋殼疲勞壽命計(jì)算的準(zhǔn)確性,在QK-WQ-6000NST-I型橋殼彎曲疲勞試驗(yàn)臺(tái)上進(jìn)行了橋殼疲勞試驗(yàn),對(duì)3個(gè)樣件進(jìn)行正弦波載荷的加載,最大載荷為2.5倍的滿載載荷,最小載荷為0.1倍的滿載載荷,直到其發(fā)生疲勞破壞為止。最終結(jié)果如表1所示,橋殼破壞位置為橋包上與半軸套管結(jié)合處,橋殼破壞時(shí)的壽命與計(jì)算結(jié)果一致,驗(yàn)證了計(jì)算結(jié)果的準(zhǔn)確性。圖7所示為橋殼最終破壞處的圖像。
表1 橋殼疲勞試驗(yàn)結(jié)果
圖7 橋殼破壞處的圖像
(1)機(jī)械結(jié)構(gòu)疲勞壽命預(yù)測(cè)主要是建立在有限元分析的載荷-時(shí)間歷程結(jié)果上,由于載荷復(fù)雜多變,需要對(duì)其進(jìn)行平均應(yīng)力修正,針對(duì)修正的載荷譜進(jìn)行循環(huán)計(jì)數(shù),計(jì)算單個(gè)循環(huán)的損傷并運(yùn)用Miner累積理論計(jì)算疲勞壽命。
(2)針對(duì)微型車橋殼進(jìn)行了疲勞壽命與臺(tái)架試驗(yàn)驗(yàn)證,結(jié)果顯示仿真計(jì)算的疲勞位置與臺(tái)架試驗(yàn)吻合,最低疲勞壽命偏差在合理的范圍內(nèi),證實(shí)仿真計(jì)算法預(yù)測(cè)橋殼的疲勞壽命能縮短開發(fā)周期,降低產(chǎn)品成本,為新產(chǎn)品的試驗(yàn)設(shè)計(jì)提供了指導(dǎo)。
[1]劉惟信.汽車車橋設(shè)計(jì)[M].北京:清華大學(xué)出版社,2004:76-132.
[2]SHAO Y.Drive axle housing failure analysis of a mining dump truck based on the load spectrum[J].Engineering Failure Analysis,2011(18):1049 -1057.
[3]秦大同,謝里陽(yáng).疲勞強(qiáng)度與可靠性設(shè)計(jì)[M].北京:化學(xué)工業(yè)出版社,2013:38-98.
[4]BS7608 -1993,F(xiàn)atigue design and assessment of steel structures[S].
[5]高翔.商用車驅(qū)動(dòng)橋橋殼強(qiáng)度及疲勞壽命分析[J].機(jī)械工程與自動(dòng)化,2010(6):50-52.
[6]高晶,宋健,朱濤.隨機(jī)載荷作用下汽車驅(qū)動(dòng)橋殼疲勞壽命預(yù)估[J].機(jī)械強(qiáng)度,2008,30(6):982 -987.
[7]朱崢濤.江鈴汽車驅(qū)動(dòng)橋橋殼有限元分析[J].汽車工程,2007,29(10):896 -899.