李新康 ，王蘇秦 ，劉潮濤 ，張 嘉 ，陳逸斐

（1. 中車青島四方機(jī)車車輛股份有限公司, 山東 青島 266111；2. 西南交通大學(xué)牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室, 四川 成都 611756）

車體作為地鐵車輛關(guān)鍵承載零部件，在運(yùn)行過(guò)程中承受來(lái)自轉(zhuǎn)向架、車鉤、車下設(shè)備以及乘客質(zhì)量等一系列復(fù)雜載荷. 新造地鐵車體由于輕量化需求，多采用鋁合金型材整體焊接而成，一節(jié)車廂焊縫多達(dá)數(shù)千條. 以傳統(tǒng)A型地鐵為例，車體側(cè)墻兩邊各開(kāi)有5門4窗，且門窗尺寸遠(yuǎn)大于干線車體，大大增加了車身強(qiáng)度薄弱位置的數(shù)量.基于以上特點(diǎn)，地鐵車體在承受復(fù)雜載荷時(shí)的疲勞問(wèn)題越來(lái)越突出.

傳統(tǒng)A型地鐵鋁合金車體長(zhǎng)約為22 m、寬度約為3 m、高度約為3 m，整備質(zhì)量約為25 t，超員載荷約為50 t. 整車的疲勞試驗(yàn)試驗(yàn)周期長(zhǎng)、經(jīng)費(fèi)高、風(fēng)險(xiǎn)大. 在ANSYS軟件中，子模型方法通過(guò)對(duì)局部關(guān)注區(qū)域的網(wǎng)格進(jìn)行細(xì)化，并在區(qū)域邊界處施加等效的邊界條件，使得計(jì)算結(jié)果更精確，計(jì)算速度更快.故可考慮截取車體受力最惡劣的部分，對(duì)該部分的邊界施加等效邊界條件后進(jìn)行疲勞試驗(yàn). 根據(jù)車體的受力狀態(tài)以及運(yùn)營(yíng)維修經(jīng)驗(yàn)，車體疲勞裂紋往往產(chǎn)生于牽枕緩部位[1]，因此可以截取車體底架端部包含牽枕緩部位的部分進(jìn)行子模型疲勞試驗(yàn).

胡杰鑫等[2]利用Monte Carlo方法，對(duì)鋁合金地鐵車體枕梁進(jìn)行了虛擬疲勞試驗(yàn)，并與實(shí)際疲勞試驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了該方法的有效性；王玉偉[3]基于實(shí)測(cè)載荷編制CRH2型車枕梁關(guān)鍵部位的應(yīng)力譜，研究了運(yùn)營(yíng)里程與損傷之間的關(guān)系；李珊珊[4]利用子模型方法，對(duì)CRH2型車枕梁結(jié)構(gòu)網(wǎng)格進(jìn)行了細(xì)分，提出了枕梁的補(bǔ)強(qiáng)方案；臧偉鋒等[5]針對(duì)飛機(jī)機(jī)身壁板設(shè)計(jì)了一種新型工裝，該工裝使得機(jī)身壁板的應(yīng)力分布與實(shí)際全尺寸機(jī)身的應(yīng)力分布相一致，并進(jìn)行了內(nèi)壓載荷疲勞試驗(yàn)；四方股份聯(lián)合航空工業(yè)飛機(jī)強(qiáng)度研究所對(duì)CRH380中間車車體全長(zhǎng)的1/3結(jié)構(gòu)進(jìn)行了氣密強(qiáng)度疲勞試驗(yàn)[6]；Gutiérrez-Carvajal等[7]設(shè)置了多種載荷工況對(duì)鋁合金車體枕梁進(jìn)行了疲勞試驗(yàn)，利用累積損傷法評(píng)估了枕梁的疲勞壽命.

本文將設(shè)定邊界條件下的靜強(qiáng)度試驗(yàn)結(jié)果與有限元計(jì)算結(jié)果進(jìn)行對(duì)比分析，驗(yàn)證了端部底架子模型（以下簡(jiǎn)稱端部底架）有限元模型的正確性. 研究了端部底架邊界條件模擬方法和試驗(yàn)加載方法，根據(jù)標(biāo)準(zhǔn)EN12663-1：2010 + A1：2014[8]（以下簡(jiǎn)稱EN 12663標(biāo)準(zhǔn)）規(guī)定的載荷進(jìn)行了疲勞損傷計(jì)算，得到了端部底架與全尺寸車體關(guān)鍵位置損傷分布的規(guī)律，驗(yàn)證了以端部底架子模型疲勞試驗(yàn)代替全尺寸車體疲勞試驗(yàn)的可行性. 通過(guò)端部底架疲勞試驗(yàn)結(jié)果對(duì)全尺寸車體的疲勞壽命進(jìn)行了評(píng)估.

1 端部底架邊界條件及載荷

1.1 邊界條件

在全尺寸車體模型中，車體的底架邊梁與車體的側(cè)墻連接，底架的端梁與車體的端墻連接，底架的尾部與車體地板及底架邊梁連接. 截取下來(lái)的端部底架長(zhǎng)5 m，與端墻、側(cè)墻、地板以及邊梁的連接均被截?cái)? 所設(shè)置的邊界條件需保證端部底架關(guān)鍵部位的應(yīng)力分布與應(yīng)力大小與全尺寸車體相一致.

在有限元計(jì)算時(shí)，容易獲取端部底架邊界節(jié)點(diǎn)的力和彎矩，該力和彎矩與施加于端部底架的外部載荷使得底架處于力平衡的狀態(tài). 在疲勞試驗(yàn)過(guò)程中，底架車鉤、空簧和中心銷處的集中力可以通過(guò)作動(dòng)器施加，由于底架邊界在全尺寸車體中變形較小，邊界處的分布力則以約束產(chǎn)生的反力施加.

根據(jù)端部底架與全尺寸車體的連接方式，在端部底架的頭部?jī)蓚?cè)與車體側(cè)墻連接處、枕梁兩側(cè)與車體側(cè)墻連接處約束端部底架的垂向自由度. 在尾部地板處約束底架的縱向和垂向自由度. 橫向約束在運(yùn)營(yíng)過(guò)程中由車體的慣性產(chǎn)生，在臺(tái)架試驗(yàn)中無(wú)法復(fù)現(xiàn)，因此選取與牽枕緩焊縫應(yīng)力無(wú)關(guān)的邊梁側(cè)面進(jìn)行約束，如圖1所示，圖中：X軸所示方向?yàn)榭v向；Y軸所示方向?yàn)闄M向；Z軸所示方向?yàn)榇瓜?

圖1 端部底架邊界條件Fig. 1 End underframe boundary conditions

1.2 試驗(yàn)載荷計(jì)算

根據(jù)EN 12663標(biāo)準(zhǔn)中關(guān)于車體疲勞載荷的規(guī)定：地鐵車體應(yīng)根據(jù)定員載荷（m1+m3）和超員載荷（m1+m4）分別設(shè)計(jì)疲勞載荷工況，其中：m1為車體整備質(zhì)量；m3為定員乘客質(zhì)量；m4為超員乘客質(zhì)量，并按其在服役壽命期間所占的比例進(jìn)行疲勞試驗(yàn).其他類型軌道交通車輛車體根據(jù)其載荷特點(diǎn)，選取恒幅定員或恒幅超員疲勞載荷進(jìn)行疲勞試驗(yàn).

由于地鐵車體定員載荷和超員載荷在服役期間所占比重難以準(zhǔn)確統(tǒng)計(jì)，為得到較為保守的結(jié)果，認(rèn)為車體在服役過(guò)程中一直承受最大載荷，車體的疲勞壽命最短，即以恒幅超員載荷（m1+m4）進(jìn)行數(shù)值計(jì)算和疲勞試驗(yàn)，m1= 28.279 t，m4= 22.440 t.

EN 12663標(biāo)準(zhǔn)中規(guī)定整車縱向和橫向疲勞載荷為 ± 0.15g，垂向載荷為g± 0.15g，但均為單個(gè)載荷分別作用1000萬(wàn)次的疲勞載荷. 對(duì)于組合載荷，EN 12663標(biāo)準(zhǔn)第6.8節(jié)推薦引用VDV 152-2016[9]標(biāo)準(zhǔn)（以下簡(jiǎn)稱VDV 152標(biāo)準(zhǔn)）.

VDV 152標(biāo)準(zhǔn)第5.3.2節(jié)中將X、Y、Z方向的載荷進(jìn)行組合，建立了3個(gè)載荷工況，每個(gè)工況循環(huán)200萬(wàn)次，見(jiàn)表1.

表1 標(biāo)準(zhǔn)之間疲勞載荷工況對(duì)比Tab. 1 Comparison of fatigue load cases between standards

VDV 152標(biāo)準(zhǔn)中載荷工況的組合與EN 12663標(biāo)準(zhǔn)中規(guī)定的3個(gè)方向載荷每個(gè)方向單獨(dú)作用1000萬(wàn)次是等效的，即在疲勞試驗(yàn)中，可用VDV 152標(biāo)準(zhǔn)規(guī)定的疲勞載荷代替EN 12663標(biāo)準(zhǔn)中規(guī)定疲勞載荷[9].

對(duì)于X方向的載荷，VDV 152標(biāo)準(zhǔn)只在一個(gè)工況下循環(huán)了200萬(wàn)次，因此需根據(jù)等損傷原則，將X方向200萬(wàn)次載荷折算成EN 12663標(biāo)準(zhǔn)中與垂向、橫向共同作用的1000萬(wàn)次載荷.

1.2.1 車鉤載荷

車鉤載荷主要源于列車牽引或制動(dòng)時(shí)的不同步，對(duì)于A型地鐵車輛，一般采用4動(dòng)2拖編組，拖車為首尾車，最大車鉤載荷發(fā)生在拖車和動(dòng)車連接的車鉤位置. 按啟動(dòng)加速度a= 1 m/s2計(jì)算，最大車鉤載荷為

根據(jù)VDV 152標(biāo)準(zhǔn)第5.3.2.1節(jié)，最大車鉤力僅在牽引（制動(dòng)）不同步的瞬間產(chǎn)生，對(duì)于動(dòng)力分散列車，車鉤力的疲勞載荷為最大車鉤力的25%. 因此，作用200萬(wàn)次車鉤載荷的動(dòng)態(tài)幅值為FXmax/4=25KN.

根據(jù)疲勞壽命-應(yīng)力曲線冪函數(shù)表達(dá)式：

式中：σ為應(yīng)力；m、C為材料常數(shù)，VDV 152標(biāo)準(zhǔn)第5.3.3.3節(jié)推薦m取4；N為循環(huán)次數(shù).

按VDV 152標(biāo)準(zhǔn)，將0.15g作用200萬(wàn)次，根據(jù)式（2）折算成1000萬(wàn)次，則車鉤載荷動(dòng)態(tài)幅值

1.2.2 中心銷牽引（制動(dòng)）載荷

依據(jù)EN 12663標(biāo)準(zhǔn)按 ± 0.15g施加，中心銷牽引（制動(dòng)）載荷如式（3）.

按式（3）將200萬(wàn)次載荷折算成1000萬(wàn)次，載荷FpivotX=24.94kN.

中心銷牽引（制動(dòng)）載荷的加載高度根據(jù)實(shí)車確定，本文中加載高度距枕梁下表面為586.5 mm.

1.2.3 中心銷橫向載荷

列車在曲線運(yùn)行時(shí)會(huì)承受來(lái)自于轉(zhuǎn)向架橫向止擋的橫向力FY，依據(jù)EN 12663標(biāo)準(zhǔn)按± 0.15g施加，如式（4）.

該載荷由中心銷和空氣彈簧共同承擔(dān)，在試驗(yàn)時(shí)僅施加中心銷承受的橫向載荷. 中心銷橫向載荷加載高度距枕梁下表面為297 mm. 根據(jù)空氣彈簧參數(shù)，每一個(gè)轉(zhuǎn)向架單個(gè)空氣彈簧承擔(dān)的橫向載荷FairsringY=3.78kN.

因此中心銷承擔(dān)的橫向載荷為

1.2.4 空簧垂向載荷

根據(jù)表1超員載荷確定空簧承受的垂向載荷.依據(jù)EN 12663標(biāo)準(zhǔn)，按g± 0.15g施加，由空氣彈簧承擔(dān).

空簧載荷均值為

空簧載荷幅值為

2 端部底架及整車有限元建模分析

2.1 端部底架及整車有限元建模

采用Hypermesh對(duì)端部底架進(jìn)行網(wǎng)格劃分，網(wǎng)格類型為Shell181. 底架所用材料為鋁合金，材料密度為2.7×10?9t/mm3，彈性模量為69 GPa，泊松比為0.3.中心銷所用材料為鋼，材料密度為7.85×10?9t/mm3，彈性模量為210 GPa，泊松比為0.3. 中心銷與枕梁下表面通過(guò)梁?jiǎn)卧˙eam188）模擬的螺栓連接. 端部底架有限元模型共劃分單元213586個(gè)，節(jié)點(diǎn)197352個(gè).底架的約束方式如圖2所示，F(xiàn)Z為Z方向的載荷.通過(guò)剛性元在車鉤座、中心銷和空簧位置施加集中力載荷.

圖2 端部底架有限元模型Fig. 2 FE model of the end underframe

全尺寸車體建模方式與端部底架相同，共劃分單元1956597個(gè)，節(jié)點(diǎn)1779596個(gè)，如圖3所示.

圖3 全尺寸車體有限元模型Fig. 3 FE model of the full-size carbody

車體在4個(gè)空簧位置約束垂向，在中心銷橫向止擋位置約束橫向. 通過(guò)剛性元在車鉤座、中心銷空簧位置施加縱向集中力載荷，垂向載荷及橫向載荷以重力場(chǎng)形式施加. 乘客質(zhì)量通過(guò)柔性元均布在地板上，空調(diào)及車下設(shè)備等考慮其對(duì)整車剛度的加強(qiáng)，通過(guò)剛性元連接到設(shè)備安裝座，車門及車窗通過(guò)剛性元連接到安裝位置.

2.2 模型正確性驗(yàn)證

試驗(yàn)過(guò)程中，利用應(yīng)變片對(duì)端部底架關(guān)鍵部位進(jìn)行應(yīng)力測(cè)量，測(cè)點(diǎn)布置見(jiàn)圖4，S1～S9為應(yīng)變片測(cè)點(diǎn)，R1～R4為應(yīng)變花測(cè)點(diǎn)，應(yīng)變花用相當(dāng)應(yīng)力進(jìn)行評(píng)定.

圖4 測(cè)點(diǎn)布置Fig. 4 Measuring points arrangement

當(dāng)載荷加到最大時(shí)，讀取各個(gè)位置的應(yīng)力值如圖5所示.

由圖5可知：測(cè)點(diǎn)S1、S2、S3、S4、S5、S9、R1、R2、R3、R4的實(shí)測(cè)應(yīng)力與計(jì)算應(yīng)力相對(duì)誤差較小，測(cè)點(diǎn)S3誤差最小為0.64%；測(cè)點(diǎn)S6、S7、S8誤差較大，其中測(cè)點(diǎn)S8誤差最大為122.2%.

圖5 實(shí)測(cè)應(yīng)力與計(jì)算應(yīng)力對(duì)比Fig. 5 Comparison of measured stress and calculated stress

測(cè)點(diǎn)S8位于枕內(nèi)加強(qiáng)筋的焊趾處，該處應(yīng)力梯度較大，且打磨后的幾何輪廓不平整，難以準(zhǔn)確確定焊趾位置，因此應(yīng)變片位置的細(xì)微差別會(huì)導(dǎo)致測(cè)量結(jié)果較大的變化；測(cè)點(diǎn)S6與測(cè)點(diǎn)S3為對(duì)稱點(diǎn)、測(cè)點(diǎn)S7與測(cè)點(diǎn)S1為對(duì)稱點(diǎn)，其仿真應(yīng)力值誤差分別為8.2%和12.2%，且測(cè)點(diǎn)S1和測(cè)點(diǎn)S3的實(shí)測(cè)應(yīng)力值和仿真應(yīng)力值相差較小，因此測(cè)點(diǎn)S6和測(cè)點(diǎn)S7誤差較大的原因可能是試驗(yàn)中的約束條件沒(méi)有仿真中的約束條件理想，約束工裝本身與車體之間存在間隙導(dǎo)致.

本文中端部底架所需評(píng)估的部位主要位于枕梁內(nèi)側(cè)加強(qiáng)框與周邊零部件的焊縫，無(wú)法直接粘貼應(yīng)變片. 測(cè)點(diǎn)S6、S7、S8與所評(píng)估部位無(wú)直接聯(lián)系，且測(cè)點(diǎn)部位在實(shí)際運(yùn)營(yíng)過(guò)程中狀態(tài)良好. 除測(cè)點(diǎn)S6、S7、S8外，其余測(cè)點(diǎn)實(shí)測(cè)應(yīng)力與計(jì)算應(yīng)力誤差較小，因此，認(rèn)為端部底架有限元模型能夠真實(shí)反映實(shí)際模型的應(yīng)力狀態(tài).

3 端部底架及整車疲勞壽命分析

3.1 疲勞理論

基于Miner線性累積損傷理論進(jìn)行疲勞計(jì)算，如式（8）.

式中：D為損傷值，當(dāng)D<1時(shí)，認(rèn)為材料沒(méi)有疲勞失效的風(fēng)險(xiǎn).；ni為第i級(jí)應(yīng)力譜下循環(huán)的次數(shù)；Ni為第i級(jí)應(yīng)力譜下的疲勞壽命.

在ANSYS下計(jì)算單位靜載荷下焊縫單元的名義應(yīng)力，在nCode Designlife中施加恒幅疲勞載荷譜.

端部底架同時(shí)承受來(lái)自車鉤座、中心銷橫向、中心銷縱向和空簧垂向載荷，在試驗(yàn)過(guò)程中其單元的主應(yīng)力方向不斷發(fā)生變化，屬于典型的多軸疲勞. 因此采用帶絕對(duì)值的最大主應(yīng)力進(jìn)行疲勞分析.

在疲勞試驗(yàn)中，端部底架承受空簧載荷均值產(chǎn)生的應(yīng)力，需進(jìn)行平均應(yīng)力修正，本文中選取FKM指導(dǎo)書(shū)[10]中規(guī)定的方法進(jìn)行平均應(yīng)力修正. 列車車體在實(shí)際服役過(guò)程中，結(jié)構(gòu)承載的載荷循環(huán)次數(shù)遠(yuǎn)遠(yuǎn)超過(guò)材料疲勞極限處的循環(huán)次數(shù)，大量小幅值載荷仍能對(duì)車體造成損傷，因此將S-N曲線的第二段延伸至應(yīng)力幅為0處，考慮小載荷循環(huán)的影響[11].

3.2 S-N曲線

本文根據(jù)Eurocode 9[12]標(biāo)準(zhǔn)建立端部底架鋁合金材料不同焊縫類型的S-N曲線. Δ σ代表循環(huán)應(yīng)力范圍. 根據(jù)端部底架焊縫類型選用Eurocode 9中7.2.3號(hào)焊縫、7.4.3號(hào)焊縫、9.4號(hào)焊縫和9.6號(hào)焊縫.焊縫類型如圖6所示.

各焊縫對(duì)應(yīng)的S-N曲線如圖7所示. 由圖7可知：各焊縫S-N曲線的拐點(diǎn)位于500萬(wàn)次循環(huán)處；S-N曲線的第一段斜率k1= ?1/3.4，第二段斜率k2= ?1/5.4.

3.3 疲勞計(jì)算結(jié)果

根據(jù)有限元模型與圖6焊縫類型的對(duì)應(yīng)關(guān)系，在nCode中將圖7各焊縫的S-N曲線分別賦予給對(duì)應(yīng)的焊縫單元，按1.2節(jié)計(jì)算的載荷進(jìn)行疲勞計(jì)算. 得到端部底架和全尺寸車體的損傷云圖如圖8，對(duì)比如表2.

表2 端部底架與全尺寸車體損傷對(duì)比Tab. 2 Damage comparison of the end underframe and the full-scale carbody

圖6 焊縫類型示意Fig. 6 Weld type diagram

圖7 焊縫S-N曲線Fig. 7 S-N curves of welding seams

圖8 損傷云圖Fig. 8 Damage nephogram

從表2可以看出：端部底架損傷最大的3個(gè)部位與全尺寸車體損傷最大的3個(gè)部位位置一致，且同一位置端部底架的損傷均大于全尺寸車體損傷.因此通過(guò)端部底架疲勞試驗(yàn)驗(yàn)證整車的疲勞強(qiáng)度可以獲得較為保守的結(jié)果.

3.4 端部底架疲勞試驗(yàn)

端部底架包括了車體運(yùn)營(yíng)過(guò)程中受力最惡劣的牽枕緩結(jié)構(gòu)，采用接口工裝和MTS液壓作動(dòng)器進(jìn)行疲勞加載，如圖9.

圖9 端部底架疲勞試驗(yàn)Fig. 9 End underframe fatigue testing

端部底架的垂向和橫向運(yùn)動(dòng)由頭部及枕梁兩側(cè)邊梁的“C”字形工裝提供約束，縱向運(yùn)動(dòng)通過(guò)固定尾部地板面提供約束，如圖10所示. 按1.2節(jié)計(jì)算的載荷進(jìn)行1000萬(wàn)次循環(huán)的疲勞試驗(yàn)，對(duì)圖8（a）所示的區(qū)域進(jìn)行滲透探傷，未發(fā)現(xiàn)裂紋，因此由端部底架的疲勞試驗(yàn)結(jié)果認(rèn)為整車的疲勞壽命滿足EN 12663標(biāo)準(zhǔn)的要求.

圖10 端部底架約束示意Fig. 10 Constraint diagram of the end underframe

4 結(jié) 論

1） 全尺寸車體的疲勞試驗(yàn)對(duì)試驗(yàn)設(shè)備、試驗(yàn)成本均有較高要求，用小尺寸子模型代替整車進(jìn)行疲勞試驗(yàn)?zāi)軌蛴行У乜s短研發(fā)周期，降低試驗(yàn)成本. 本文采用端部底架子模型代替全尺寸車體的疲勞損傷評(píng)估結(jié)果偏于保守，從數(shù)值仿真和臺(tái)架試驗(yàn)的角度驗(yàn)證了方法對(duì)車體疲勞壽命進(jìn)行評(píng)估是可行的.

2） 使用子模型進(jìn)行疲勞試驗(yàn)時(shí)，需保證子模型的邊界遠(yuǎn)離可能的裂紋潛在部位，且保證子模型結(jié)構(gòu)潛在危險(xiǎn)點(diǎn)與全尺寸車體結(jié)構(gòu)保持一致. 邊界條件的設(shè)置可根據(jù)邊界處剛度的方向特性，選取剛度最大的方向進(jìn)行約束，釋放他方向的約束. 通過(guò)有限元仿真設(shè)置的邊界條件與試驗(yàn)進(jìn)行對(duì)比，獲得正確的有限元數(shù)值分析模型，為車體的壽命評(píng)估提供依據(jù).

3） 結(jié)合名義應(yīng)力法S-N曲線和線性損傷理論，分別計(jì)算了端部底架與全尺寸車體的疲勞損傷. 發(fā)現(xiàn)該方法依賴于網(wǎng)格質(zhì)量和S-N曲線的選取，針對(duì)車體和其焊縫特點(diǎn)，在計(jì)算中均進(jìn)行了考慮，因此采用的方法適合車體的壽命評(píng)估.

4） 采用端部底架子模型的疲勞試驗(yàn)與數(shù)值計(jì)算，并與全尺寸車體結(jié)構(gòu)的對(duì)比分析，表明該車體能夠滿足EN 12663標(biāo)準(zhǔn)規(guī)定的疲勞壽命要求.

致謝：鋁合金車體疲勞壽命評(píng)估及試驗(yàn)方法研究（2682019CX45）.