考慮氣動載荷的動車組鋁合金司機(jī)室焊接接頭應(yīng)力因子分析

謝素明，齊碩 ，程亞軍，范文杰

(1.大連交通大學(xué) 機(jī)車車輛工程學(xué)院，遼寧 大連 116028；2.中車長春軌道客車股份有限公司 工程實驗室，吉林 長春 130062)*

列車高速會車和穿越隧道時，車體頭部和尾部承受的氣動壓力瞬間發(fā)生劇烈變化，這種頻繁交變的氣動載荷已成為高速動車組車體設(shè)計過程中重點關(guān)注的外部載荷[1].

當(dāng)前，在高速列車空氣動力學(xué)試驗和氣動載荷下車體疲勞問題方面開展了比較深入地研究.熊小慧[2]等闡述了膠濟(jì)線CRH2型動車組列車交會空氣壓力波實車測試情況，總結(jié)出CHR2型動車組上的交會壓力波幅值近似與同型號等速交會動車組運行速度的平方成正比；盧耀輝等[3]將氣動載荷與其他疲勞載荷相組合，基于Goodman疲勞曲線對車體疲勞強度進(jìn)行評定，指出氣動載荷容易導(dǎo)致側(cè)墻疲勞失效；余思均[4]開展了高速列車車體受氣動載荷作用的瞬態(tài)過程研究，分析了氣動載荷的瞬態(tài)作用對車體側(cè)墻結(jié)構(gòu)強度的影響.

目前高速動車組制造部門進(jìn)行車體焊接接頭設(shè)計時執(zhí)行EN15085-3:2007標(biāo)準(zhǔn)[5]，要求車體焊接接頭即要滿足疲勞壽命設(shè)計要求，還需確定接頭的應(yīng)力因子和應(yīng)力狀態(tài)等級.為此，謝素明等[6]利用DVS1608-2011及BS EN15085-3標(biāo)準(zhǔn)，在BS EN 12663標(biāo)準(zhǔn)提供的加速度疲勞載荷作用下，研究動車組鋁合金車體焊接接頭的應(yīng)力因子和應(yīng)力狀態(tài)等級.然而，由于DVS 1608-2011標(biāo)準(zhǔn)僅適用于ENAW-6005A和EN AW-6082的擠壓型材、板材電弧焊接接頭承受載荷循環(huán)次數(shù)為107的情況[7]，實際上高速動車組車體的氣動疲勞載荷的循環(huán)次數(shù)達(dá)不到這樣高，所以，當(dāng)考慮氣動疲勞載荷時DVS 1608-2011標(biāo)準(zhǔn)就無法用于動車組車體焊接接頭的應(yīng)力因子分析.

在研究EN1999-1-3:2007[8]和IIW-2008[9]標(biāo)準(zhǔn)的基礎(chǔ)上，本文以某高速動車組鋁合金車體為分析對象，基于疲勞損傷等效原則，研究氣動載荷和加速度疲勞載荷共同作用下的動車組車體司機(jī)室焊接接頭疲勞壽命，進(jìn)而根據(jù)BS EN15085-3標(biāo)準(zhǔn)計算出焊接接頭的應(yīng)力因子，確定焊接接頭的應(yīng)力狀態(tài)等級.

1 車體焊接接頭應(yīng)力因子的分析方法

由于國際焊接學(xué)會IIW-2008文件提供了較多的鋁合金焊接接頭的S-N曲線數(shù)據(jù)，所以鋁合金車體設(shè)計部門在鋁合金結(jié)構(gòu)焊接接頭疲勞設(shè)計時經(jīng)常采用該文件的名義應(yīng)力法的疲勞性能參數(shù).隨著鋁合金車體研制質(zhì)量的升級，現(xiàn)已開展將鋁合金結(jié)構(gòu)疲勞設(shè)計的英國標(biāo)準(zhǔn)BS EN 1999-1-3: 2007用于高速鋁合金車體焊接結(jié)構(gòu)的設(shè)計工作.

BS EN 1999-1-3: 2007適用于鋁合金車體常用的EN AW-6005A和EN AW-6082的鋁合金擠壓型材、板材以及鍛造件母材和電弧焊接結(jié)構(gòu)的疲勞強度評價.其疲勞強度數(shù)據(jù)曲線由試驗數(shù)據(jù)的均值減去2個標(biāo)準(zhǔn)差表示，應(yīng)力與壽命次數(shù)的關(guān)系曲線以對數(shù)表達(dá)，見圖1.選擇主應(yīng)力范圍作為名義應(yīng)力法評估焊接接頭疲勞強度的依據(jù).由于車體結(jié)構(gòu)有限元模型中包括了部件幾何結(jié)構(gòu)不連續(xù)的主要細(xì)部結(jié)構(gòu)，所以，從有限元分析結(jié)果中提取出焊接接頭的主應(yīng)力可視為修正標(biāo)稱應(yīng)力.焊接接頭的應(yīng)力選取點位置應(yīng)距離焊縫一定距離，建議10 mm.

圖1 BS EN 1999-1-3中l(wèi)ogΔσ-logN曲線

圖1中a為疲勞強度曲線；b為參考疲勞強度(2×106次循環(huán)Nc對應(yīng)應(yīng)力范圍Δσc)；c為常幅疲勞極限(5×106次循環(huán)ND對應(yīng)應(yīng)力范圍ΔσD)；d為截止極限(108次循環(huán)NL對應(yīng)應(yīng)力范圍ΔσL)；m1和m2為曲線Δσ-N的反斜率.該曲線的第一個拐點為5×106，第二個拐點為108.曲線的分段關(guān)系式為

(1)

(5×106～108)

(2)

式中:γFf為荷載譜和響應(yīng)分析中考慮到不確定性的分項系數(shù)(建議取值1.0)；γMf為材料和施工不確定性的分項系數(shù)(建議取值1.0)；對于焊接接頭，m2=m1+2.除全熔透雙面磨平的對接接頭的logΔσ-logN曲線的m1取值為7外，其他焊接接頭的m1取值為3.4或4.3.

國際焊接學(xué)會IIW-2008文件“焊接接頭與部件的疲勞設(shè)計”提供了材料屈服強度高達(dá)700MPa的鐵素體/珠光體或貝氏體結(jié)構(gòu)鋼鍛造或擠壓產(chǎn)品制成的焊接結(jié)構(gòu)和常用于制造焊接結(jié)構(gòu)的鋁合金的至少有95%存活率的抗疲勞數(shù)據(jù).文件中給出焊接接頭雙對數(shù)坐標(biāo)系下的應(yīng)力范圍Δσ與達(dá)到疲勞的循環(huán)次數(shù)N雙斜率的曲線，第一個曲線拐點對應(yīng)的循環(huán)次數(shù)Nc= 107，第二個拐點即為曲線截止循環(huán)次數(shù)Ncut= 1×109.與BS EN 1999-1-3: 2007相同，將N= 2×106循環(huán)次數(shù)對應(yīng)的疲勞強度值記為FAT等級.Δσ-N曲線的兩階段的函數(shù)關(guān)系式均可寫成

N=C/Δσm

(3)

當(dāng)Δσ>Δσ1時，上式中C為C1，m為m1，值為3；當(dāng)Δσ2≤Δσ≤Δσ1時，C為C2，m為m2，常幅載荷時值為22；變幅載荷時值為5.Δσ1和Δσ2為Δσ-N曲線的拐點；C1和C2為常系數(shù)；m1和m2為Δσ-N曲線斜率.

兩標(biāo)準(zhǔn)中焊接接頭細(xì)部、承載方向、質(zhì)量要求一致，疲勞強度等級相同或接近的四個典型接頭如圖2所示，圖2(a)中類型I對接接頭的承載方向與焊縫平行，連接板對接部位沒有圓角半徑；類型II對接接頭的承載方向與焊縫垂直，焊縫全溶透且雙面磨平；類型III搭接接頭的承載方向與焊縫垂直，焊縫潛在裂紋初始位置為焊喉；類型IV接頭的承載方向與焊縫平行，縱向角焊縫連接板長度大于100 mm且雙面焊接.四個接頭的logΔσ-logN曲線典型參數(shù)數(shù)據(jù)見表1，其中類型I和II的logΔσ-logN曲線如圖3所示.

(a) 類型I (b) 類型II

(c) 類型III (d) 類型IV圖2 焊接接頭細(xì)部與承載方向

結(jié)合表1和圖3可以看出：接頭疲勞強度等級相同的接頭I和II，以及接頭疲勞強度等級相近的接頭III和IV，它們的循環(huán)次數(shù)107對應(yīng)的應(yīng)力變化范圍IIW中的值要小于BS EN 1999-1-3中的值.

表1 接頭疲勞強度等級接近的疲勞參數(shù)

(a) 類型I

(b) 類型II圖3 接頭疲勞強度等級相同的logΔσ-logN曲線

由于軌道車輛鋁合金車體加速度疲勞設(shè)計載荷來自EN12663，加載次數(shù)均為107.所以，采用IIW文件評估客車車體焊接接頭疲勞強度偏于保守.

目前高速動車組鋁合金車體的疲勞載荷工況可主要歸納為兩類：一類是由多個單載荷工況組成的疲勞工況，加載次數(shù)為107；另一類是由多個疲勞工況組成且各疲勞載荷的加載次數(shù)不同，見表2.

表2 車體兩種疲勞載荷工況

兩標(biāo)準(zhǔn)的名義應(yīng)力方法評定焊接接頭疲勞強度的過程相同.承受第一類疲勞載荷工況的焊接接頭疲勞強度主要步驟如下：

(1)分析每一單工況載荷作用下焊接接頭評估點主應(yīng)力大小和方向，并將所有載荷工況下評估點主應(yīng)力值最大的確定為最大主應(yīng)力σmax，其對應(yīng)的方向余弦(nx、ny、nz)作為基本應(yīng)力矢量方向；

圖4 主應(yīng)力變化范圍的計算過程示意圖

(4)

(4)判斷評估點的主應(yīng)力變化范圍Δσ與ΔσR的大小.同時，若Δσ小于ΔσR，且大于標(biāo)準(zhǔn)中疲勞曲線的截止應(yīng)力變化范圍時，需利用式(1)～(3)計算損傷；

(5)當(dāng)Δσ小于ΔσR時，接頭疲勞評估結(jié)束.否則，需要修改局部結(jié)構(gòu)或者提高焊接接頭的疲勞等級，返回(1)重新進(jìn)行計算，直到滿足設(shè)計要求.

當(dāng)焊接接頭承受第二類疲勞載荷時，需對單疲勞工況執(zhí)行上述步驟(1)～(4)，獲得接頭的每一單疲勞工況的損傷ni/Ni，然后再進(jìn)行損傷累積

(5)

若D小于1.0時，接頭疲勞評估結(jié)束.否則，需要修改局部結(jié)構(gòu)或者提高焊接接頭的疲勞等級重新計算，直到滿足設(shè)計要求.

執(zhí)行EN15085-3:2007標(biāo)準(zhǔn)進(jìn)行高速動車組車體焊接接頭設(shè)計時，車體焊接接頭首先要滿足疲勞壽命設(shè)計要求，然后還需確定接頭的應(yīng)力因子和應(yīng)力狀態(tài)等級.對于承受第一類疲勞載荷工況的焊接接頭的應(yīng)力因子計算比較簡單，Δσ與ΔσR的比值即為接頭應(yīng)力因子；當(dāng)接頭承受第二類疲勞載荷工況時情況相對復(fù)雜，需先依據(jù)累積損傷D求出等效循環(huán)次數(shù)Neq

(6)

式中，NT為每一單疲勞工況的加載次數(shù)之和.然后利用式(1)～(3)獲得等效應(yīng)力σeq和NT對應(yīng)的應(yīng)力變化范圍ΔσRT.接頭的應(yīng)力因子為等效應(yīng)力σeq與ΔσRT的比值.

2 車體焊接接頭應(yīng)力因子分析

某動車組鋁合金車體長度為25 498 mm，寬度為3 360 mm，高度為4 050 mm(距軌面)，車輛定距為17 800 mm，車輛定員質(zhì)量為44 t.考慮氣動載荷的車體疲勞載荷工況屬于第二類，參見表2.氣動載荷作用于車體側(cè)墻和車頂?shù)耐獗砻?，此時NT應(yīng)為3 200萬次.該車體有限元分析模型以任意四節(jié)點薄殼單元為主，單元總數(shù)為2 467 282，節(jié)點總數(shù)為2 240 556，模型重10.5t.車體有限元分析模型如圖5所示.

圖5 車體有限元模型及典型接頭示意

整車車體由不同形狀的鋁合金型材及板材組焊而成，存在大量形式不同和方向各異的焊縫，車體型材典型接頭參見圖5.基于IIW-2008和BS EN 1999-1-3: 2007標(biāo)準(zhǔn)開展如此眾多焊縫的疲勞分析及接頭應(yīng)力狀態(tài)等級的評定工作，評估周期長.為提高計算效率，應(yīng)考慮編制程序來實現(xiàn).

由于車體結(jié)構(gòu)強度分析采用ANSYS軟件完成的，同時ANSYS批處理模式執(zhí)行分析過程不需要顯示模型就可以在后臺執(zhí)行APDL編寫的分析程序流，所以為便于讀取ANSYS的分析結(jié)果數(shù)據(jù)，編制接頭疲勞分析及應(yīng)力因子計算程序的基本思路確定為：利用C#開發(fā)可視化程序，以批處理模式啟動ANSYS軟件，然后讀取編寫的APDL程序流進(jìn)行計算[10].其中，C#程序和ANSYS軟件之間的通信為重要核心問題.

基于APDL參數(shù)化語言編寫的程序包括：①根據(jù)C#編寫的可視化頁面上用戶的選擇確定鋁合金車體需要分析的關(guān)切焊縫、疲勞載荷工況以及焊縫在不同工況下所對應(yīng)的疲勞強度等級參數(shù)；②利用*GET語句訪問ANSYS數(shù)據(jù)庫，從ANSYS計算結(jié)果rst文件中提取整體坐標(biāo)系下所有疲勞計算工況的焊縫節(jié)點的第一主應(yīng)力；③程序包含3層循環(huán)，每層循環(huán)通過循環(huán)語句*DO-*ENDDO編制而成，同時循環(huán)過程中穿插著多個條件語句*IF對整個循環(huán)體進(jìn)行完善，外側(cè)循環(huán)以焊縫組為循環(huán)控制變量，中間層循環(huán)以載荷工況編號為循環(huán)控制變量，內(nèi)層循環(huán)以焊縫分組內(nèi)的節(jié)點編號為循環(huán)控制變量.當(dāng)完成焊縫組別內(nèi)的所有節(jié)點的計算后，程序?qū)⒆詣咏Y(jié)束；④通過創(chuàng)建宏文件語句*CFWRITE將評估的每條焊縫在每個疲勞載荷下的分析結(jié)果以txt文件的形式輸出保存.

基于C#編寫交互界面及ANSYS事件處理程序：交互界面要考慮到使用的易用性，同時還必須提供分析焊接接頭必要的參數(shù)輸入窗口.軟件會自動讀取車體有限元模型中的焊縫名稱以及所加載的疲勞載荷工況供用戶選擇.交互界面編寫完成后，在計算按鈕的Click事件中新建進(jìn)程啟動ANSYS.交互界面如圖6所示.軟件會在后臺使用Stream Write和Stream Reader類來讀取和寫入由APDL腳本輸出的txt結(jié)果文檔，并對其中的數(shù)據(jù)進(jìn)行處理，將關(guān)心的數(shù)據(jù)顯示在交互界面上[11].關(guān)切的應(yīng)力云圖以及位移云圖會通過Image. From File方法顯示在交互界面上，參見圖6.

圖6 接頭應(yīng)力因子計算程序的交互界面

在疲勞載荷工況作用下，利用編制的IW-2008標(biāo)準(zhǔn)和BS EN 1999-1-3: 2007標(biāo)準(zhǔn)評估程序，確定了車體司機(jī)室疲勞薄弱的四條焊縫分別為：司機(jī)室立柱與邊梁焊縫(焊縫1)、司機(jī)室立板加強板與邊梁焊縫(焊縫2)、司機(jī)車車窗立柱與邊梁焊縫(焊縫3)、司機(jī)室立板與邊梁焊縫(焊縫4).

由于各疲勞工況的載荷方向和作用位置不同，因此，同一條焊縫在不同疲勞載荷工況下的主應(yīng)力矢量方向不同(例如焊縫3在不同疲勞載荷下的主應(yīng)力矢量圖如圖7所示).

圖7 不同載荷下焊縫3的主應(yīng)力方向

根據(jù)疲勞載荷工況下的主應(yīng)力矢量方向與焊接接頭的類型確定出焊縫在每個工況下的FAT等級.表3給出了這四條焊縫在每工況下的FAT等級、主應(yīng)力大小及損傷.由表3可以看出：這四條焊縫的累積損傷值均是IIW的計算結(jié)果大于EN 1999-1-3的計算結(jié)果；焊縫1的累積損傷最大，值為0.753.

表3 車體關(guān)鍵焊縫的累積損傷

表4給出了這四條焊縫的應(yīng)力因子以及應(yīng)力狀態(tài).BS EN15085-3標(biāo)準(zhǔn)規(guī)定：接頭應(yīng)力因子大于等于0.9時，應(yīng)力狀態(tài)等級為高；應(yīng)力因子大于等于0.750且小于0.9時，應(yīng)力狀態(tài)等級為中.由表4可以看出：四條焊縫IIW的應(yīng)力因子計算結(jié)果均大于EN 1999-1-3的計算結(jié)果；使用IIW-2008標(biāo)準(zhǔn)評定的四條焊縫的應(yīng)力狀態(tài)等級均為高；而用EN 1999-1-3評定時焊縫1的應(yīng)力因子為0.94，應(yīng)力狀態(tài)等級為高；其他三條焊縫的應(yīng)力因子均小于0.9，應(yīng)力狀態(tài)等級均為中.

表4 車體關(guān)鍵焊縫的應(yīng)力狀態(tài)等級

3 結(jié)論

(1)兩標(biāo)準(zhǔn)中疲勞強度等級相同及相近的接頭， IIW中的循環(huán)次數(shù)107對應(yīng)的應(yīng)力變化范圍的值小于BS EN 1999-1-3中的相應(yīng)值.由于來自EN12663的軌道車輛車體加速度疲勞設(shè)計載荷的加載次數(shù)為107，所以，采用IIW文件對焊接接頭進(jìn)行疲勞設(shè)計比采用BS EN 1999-1-3偏于保守；

(2)結(jié)合C#和APDL并以批處理模式啟動ANSYS，開發(fā)了適用于多疲勞載荷共同作用下的基于名義應(yīng)力法進(jìn)行焊接接頭疲勞分析及應(yīng)力因子計算程序，可提高疲勞評定工作效率，降低分析工作強度；

(3)氣動載荷和加速度疲勞載荷同時作用下，兩標(biāo)準(zhǔn)分析的司機(jī)室焊接接頭的累積損傷均大于0.7，且是由氣動疲勞載荷引起的；焊縫1(司機(jī)室立柱與邊梁焊縫)的累積損傷和應(yīng)力因子最大，分別為0.753和0.987.所以，高速動車組司機(jī)室焊接結(jié)構(gòu)抗疲勞設(shè)計時應(yīng)重點關(guān)注氣動疲勞載荷.