基于有限元法的車架仿真分析與試驗(yàn)研究

黃麗芳，尹曉春，姜亮，肖欣

(南京理工大學(xué) 理學(xué)院，江蘇 南京 210094)

0 引言

煉鋼廠的鋼包車用于將鋼水包運(yùn)輸?shù)骄珶捾囬g或連鑄生產(chǎn)線。由于鋼水包內(nèi)裝載著高溫液態(tài)鋼水，一旦在運(yùn)輸過程出現(xiàn)轉(zhuǎn)向失靈或者車體斷裂，極可能導(dǎo)致盛裝鋼水的鋼包傾翻，造成嚴(yán)重的安全事故[1]。因此，提高鋼包車的結(jié)構(gòu)可靠性，降低鋼包車的損壞，是鋼包車設(shè)計的主要任務(wù)，具有重要的工程意義。

某鋼廠5號250 t鋼包車于2012年投入生產(chǎn)。使用中，車體兩側(cè)出現(xiàn)了長約1 m的裂縫(圖1和圖2黑線)，開裂位置位于傳動側(cè)距軌面標(biāo)高1.51 m、與貫穿橫板焊接的東、西兩側(cè)立板焊縫處。在原始設(shè)計中，鋼包車車體的最大計算應(yīng)力為245MPa，位于臺車中部。實(shí)際斷裂處的設(shè)計應(yīng)力約為50MPa，斷裂處計算的應(yīng)力約為屈服強(qiáng)度的1/7，安全系數(shù)足夠。而在實(shí)際工作中，斷裂位置與設(shè)計計算不符。因此，需要進(jìn)行精確的三維有限元分析，重新校核強(qiáng)度，找到開裂原因，進(jìn)而提出結(jié)構(gòu)優(yōu)化的設(shè)計方案，修復(fù)結(jié)構(gòu)，并使修復(fù)后的鋼包車滿足安全生產(chǎn)要求。

圖1 現(xiàn)場車體西側(cè)開裂處

圖2 現(xiàn)場車體東側(cè)開裂處

1 鋼包車結(jié)構(gòu)與材料

鋼包車與鋼包整體結(jié)構(gòu)是由車架整體與鋼包、鋼水組成，其中鋼包與鋼水作為輔助結(jié)構(gòu)，用于驗(yàn)證鋼包車結(jié)構(gòu)的可靠性與穩(wěn)定性。鋼包車主要由彈簧支座(用于放置鋼包)、車架體、蓋板、車輪組(包括從動車輪和主動車輪)、溜槽等結(jié)構(gòu)組成。該鋼包車為全鋼焊接的空間結(jié)構(gòu)，體積龐大，噸位重，車體承受鋼水包質(zhì)量和自質(zhì)量、運(yùn)輸過程的沖擊與振動以及可能的溫度載荷。

車體長12.57 m、寬6.15 m、高(含鋼水罐)5.75 m，體積68.197 m3，總質(zhì)量535.35 t。鋼包車車架材料為Q345-B，密度ρ=7 850kg/m3，彈性模量E=206GPa，泊松比μ=0.3。各種規(guī)格材料主要設(shè)計參數(shù)見表1。

表1 Q345-B許用應(yīng)力一覽表

2 有限元理論基礎(chǔ)

鋼包車車架的靜態(tài)特征主要體現(xiàn)在車架結(jié)構(gòu)的靜剛度上。靜剛度是結(jié)構(gòu)在特定的動態(tài)激擾下抵抗變形的能力，一般用結(jié)構(gòu)在靜載荷作用下的變形情況來衡量。分析車架結(jié)構(gòu)靜剛度的主要目的是：

1)避免車架上其他功能部件工作時由于過大的變形產(chǎn)生干涉現(xiàn)象，導(dǎo)致鋼包傾覆；

2)保證鋼包車不會因剛度不足等原因明顯地影響澆鑄設(shè)備之間的連接。

在已知外載和各節(jié)點(diǎn)的約束條件下，就可以得到各節(jié)點(diǎn)的位移δ、應(yīng)變ε及應(yīng)力σ。對鋼包車結(jié)構(gòu)采用8節(jié)點(diǎn)實(shí)體單元進(jìn)行三維建模，能夠準(zhǔn)確體現(xiàn)車架的應(yīng)力變化情況和位移變化情況[2]。

8節(jié)點(diǎn)實(shí)體單元的幾何矩陣是[3]

B=[B1B2B3…B8]

(1)

式中

根據(jù)形函數(shù)的定義可以發(fā)現(xiàn)，形函數(shù)的計算精度與單元體的尺寸有關(guān)。豎向隔板和橫向墊板之間的垂直焊接處存在應(yīng)力集中現(xiàn)象，該處應(yīng)力場變化劇烈。因此為了提高計算結(jié)果準(zhǔn)確性，在網(wǎng)格劃分的處理上，縮小該部位的網(wǎng)格尺寸，即該部位網(wǎng)格加密，使載荷準(zhǔn)確傳遞，并提高計算結(jié)果準(zhǔn)確性。

8節(jié)點(diǎn)實(shí)體單元的形函數(shù)表示為

其中i=1,2，…，8。

應(yīng)力矩陣為

S=D×[B1B2B3…B8]

(2)

式中：

式中：μ為泊松比；E為彈性模量。

其中：Be{σ}e={ε},Se{σ}e={δ}，按照虛功原理方程為

?{ε}T{σ}dxdydz={ε}T{R}e

(3)

有單元剛度矩陣為

(4)

通過單元剛度矩陣可得整體剛度矩陣K，整體剛度矩陣方程為

K{σ}={R}

(5)

3 有限元仿真

3.1 建立幾何模型

模型的建立，一方面要盡量保留原實(shí)體結(jié)構(gòu)的細(xì)節(jié)，反映局部危險部位對車架應(yīng)力的影響，另一方面，為減少有限元計算量，亦需要對危險部位進(jìn)行結(jié)構(gòu)簡化[4]。采用SolidWorks建立幾何模型，在滿足分析要求的前提下，對現(xiàn)有結(jié)構(gòu)的一些零件進(jìn)行簡化，忽略小零件的作用。根據(jù)鋼包車整體圖樣的設(shè)計，主要簡化了蓋板和車輪組，去掉了溜槽部位。簡化后的整體主要由鋼包、彈簧座、車輪組和車架組成。圖3為經(jīng)過簡化后的鋼包車模型。

圖3 車體幾何模型

3.2 劃分網(wǎng)格

在車架結(jié)構(gòu)的有限元模型建立所選擇的單元類型中，早期時主要采用殼單元進(jìn)行結(jié)構(gòu)離散，但殼單元的有限元模型計算的結(jié)果誤差較大，并對復(fù)雜結(jié)構(gòu)的模擬能力有限[5]。現(xiàn)在主要選擇8節(jié)點(diǎn)實(shí)體單元進(jìn)行有限元網(wǎng)格劃分，8節(jié)點(diǎn)實(shí)體單元(Solid 185)對應(yīng)力傳遞和結(jié)構(gòu)變形都能夠準(zhǔn)確表述，并且加上荷載后，能夠模擬實(shí)際工況中可能出現(xiàn)的變形情況[6]。根據(jù)網(wǎng)格劃分的規(guī)則[7]，保證計算的準(zhǔn)確性：1)在應(yīng)力分布均勻處進(jìn)行規(guī)則離散；2)在與車輪組連接的不規(guī)則處進(jìn)行手動過渡，保證網(wǎng)格疏密有致；3)在實(shí)際工況中斷裂部位加密網(wǎng)格劃分；4)有效地控制總體網(wǎng)格數(shù)目和質(zhì)量。圖4顯示了開裂部位周圍的網(wǎng)格劃分。

圖4 開裂周邊網(wǎng)格分布

3.3 邊界條件

由于鋼包車在運(yùn)行過程中，主要承擔(dān)的是鋼包與鋼水及其自身的靜力作用。主要有3個約束條件：1)實(shí)際情況中鋼包是可以脫離鋼包車的，因此彈簧支座與鋼包之間建立接觸約束，約束鋼包在Z方向(即寬度方向)的位移；2)在車輪組處，軸輪與大車輪外殼兩組結(jié)構(gòu)建立接觸約束，車輪組與車架之間建立剛性連接；3)為了準(zhǔn)確地模擬鋼包車在運(yùn)行時的受力與變形，將車輪組放置在軌道上，建立接觸約束，并對軌道底面施加完全約束。主要約束位置如圖5所示。

圖5 邊界約束布置

3.4 施加荷載

鋼包車荷載參數(shù)見表2。

表2 鋼包車荷載參數(shù) 單位：t

在分析鋼包車的結(jié)構(gòu)強(qiáng)度時，考慮了在各種工況中的受力情況。本文展示了對其中最危險的一種工況的分析，即考慮作用在鋼包車上的所有荷載。荷載的傳力路線主要是：鋼包和鋼水的質(zhì)量最初施加在車架的彈簧支座處，隨后再由彈簧支座通過剛性連接傳遞給車架其他部分，最后通過車架整體與車輪組的剛性連接傳遞至車輪組及軌道上。

根據(jù)《YB/T4224—2010冶金用鋼水罐車和鐵水罐車技術(shù)規(guī)范》 行業(yè)標(biāo)準(zhǔn)，初步設(shè)計計算只選擇最惡劣的工況，即鋼包滿載時垂直沖擊工況，并且考慮鋼包車在運(yùn)行過程中的橫向沖擊荷載[8]：F(豎向模型施加載荷)=G(滿載鋼包自質(zhì)量載荷)×h(動載沖擊系數(shù))=250×1.5=375t。將375t的豎向模型施加載荷作為面載荷均勻施加在罐座載荷面上，并且在X方向施加1.5m/s2的加速度。將有限元模型導(dǎo)入ANSYS結(jié)構(gòu)靜力學(xué)分析模塊，進(jìn)行計算。

4 有限元計算結(jié)果與分析

4.1 鋼包車有限元云圖分析

圖6顯示了車架體的整體應(yīng)力云圖。位于車架體腹板與邊部相焊接處，傳動側(cè)東側(cè)和西側(cè)各構(gòu)件交接處，最大應(yīng)力達(dá)到413MPa，超過Q345-B的屈服強(qiáng)度。因此，可以推測傳動側(cè)東側(cè)和西側(cè)裂紋是由此處開裂并延展的。在車架傳動西側(cè)板，最大應(yīng)力為557 MPa，在車架傳動?xùn)|側(cè)板，最大應(yīng)力為489 MPa，兩側(cè)最大應(yīng)力遠(yuǎn)超Q345-B的屈服強(qiáng)度。由于底部應(yīng)力傳遞到豎向隔板，導(dǎo)致豎向隔板和橫向墊板之間的垂直焊接處出現(xiàn)了裂縫。圖7顯示了車架體的整體位移分布。由圖可見，中間腹板部分位移最大，車體與車輪組連接處位移最小，車架體的位移分布符合受力情況。根據(jù)有限元結(jié)果，沒有出現(xiàn)車架體發(fā)生位移突變現(xiàn)象，雖然車架體表面有局部破裂，但是不影響整體變形。

圖6 車架應(yīng)力分布圖

圖7 車架位移分布圖

鋼包車的車架體結(jié)構(gòu)復(fù)雜，為了清晰地呈現(xiàn)有限元計算結(jié)果，繪制應(yīng)力分布曲線與位移分布曲線。選取沿X方向具有應(yīng)力分布代表性的某一直線提取數(shù)據(jù)，該直線通過東、西側(cè)破裂點(diǎn)。如圖8和圖9所示(本刊黑白印刷，相關(guān)疑問咨詢作者)，車架體東、西側(cè)的應(yīng)力曲線和位移曲線基本吻合，由于車架體的南北方向(X方向)不是對稱結(jié)構(gòu)，東西方向(Z方向)為對稱結(jié)構(gòu)，在理論計算下，對稱結(jié)構(gòu)的應(yīng)力應(yīng)變結(jié)果也基本對稱，證明有限元模型契合實(shí)際且可靠。圖8顯示，車架體應(yīng)力在車體破損處急劇升高，超過了鋼的屈服強(qiáng)度，且明顯高于兩側(cè)應(yīng)力，產(chǎn)生了顯著的應(yīng)力梯度。

圖8 應(yīng)力分布曲線圖

圖9 位移分布曲線圖

4.2 鋼包車開裂原因

通過現(xiàn)場照片和現(xiàn)場察看可以明顯看出：焊縫裂紋沿焊縫方向延伸，焊肉完整但與母材脫離，母材表面完好無撕裂。據(jù)此分析，鋼包車車體存在未焊透、未熔合等焊接缺陷。未焊透會降低焊縫的疲勞強(qiáng)度，可能使裂紋源造成焊縫破壞。通過有限元模型的計算，發(fā)現(xiàn)鋼包車在實(shí)際使用過程中的開裂部分，與有限元模型計算所得的結(jié)果相吻合，即出現(xiàn)裂縫處達(dá)到了Q345-B的屈服強(qiáng)度，力從折角連接處傳遞到了上部焊接處，導(dǎo)致焊接處出現(xiàn)開裂。

4.3 車架結(jié)構(gòu)優(yōu)化方案

由于現(xiàn)場條件的限制，無法快速地重新設(shè)計鋼包車，即在無法改變鋼包車的長、寬、高的情況下，初步按照以下步驟對現(xiàn)有損壞的鋼包車進(jìn)行調(diào)整[9]：

1)將損壞部位去除，包括部分貫穿橫板和裂紋處的上部部分鋼板；

2)上部鋼板(立板厚度為50mm，中間增加抗彎人字筋)按圖樣要求預(yù)制好，將與原貫穿橫板及中部連接梁等結(jié)合部位處打剖口。添加構(gòu)件具體部位如圖10所示，鋼板結(jié)構(gòu)尺寸如圖11所示；

圖10 優(yōu)化添加的構(gòu)件

圖11 加固構(gòu)件結(jié)構(gòu)圖

3)與原車體焊接，剖口處必須焊透，焊后超聲波探傷檢查；

4)新焊接部位去應(yīng)力處理。

結(jié)構(gòu)優(yōu)化后，相同約束條件和沖擊荷載下，車架東西側(cè)板部位的最大應(yīng)力從413 MPa降低至180 MPa左右(圖12-圖13)。通過添加構(gòu)件的優(yōu)化方式，極大地提高了車架的承載能力。

圖12 優(yōu)化后西側(cè)局部應(yīng)力圖

圖13 優(yōu)化后東側(cè)局部應(yīng)力圖

5 結(jié)語

1)結(jié)合某鋼廠的實(shí)際工況，有限元分析結(jié)果與現(xiàn)場實(shí)際工況一致，表明有限元模型準(zhǔn)確可靠，計算結(jié)果真實(shí)可信。有限元分析表明，在實(shí)際工況下，鋼包車整體位移分布無異常，局部破裂處位移沒有突變。局部應(yīng)力突然變大，應(yīng)力梯度明顯，符合現(xiàn)場破損情況。

2)基于鋼包車的應(yīng)力和應(yīng)變分析，對車體結(jié)構(gòu)進(jìn)行優(yōu)化，使局部應(yīng)力從400MPa下降至180MPa左右，證明該優(yōu)化方案符合結(jié)構(gòu)設(shè)計要求，結(jié)構(gòu)優(yōu)化過程操作簡單，經(jīng)濟(jì)性良好。

3)車體的有限元建模方法適用于各種大型結(jié)構(gòu)件，能夠準(zhǔn)確分析車體的強(qiáng)度和剛度性能，為車體設(shè)計、校核與優(yōu)化提供可靠的理論依據(jù)。