鐵路無縫線路溫度應(yīng)力的有限元分析

潘文彬，葉 淵，韓洪江，楊 韜，宋家旺*

（1.浙江工業(yè)大學(xué) 機械工程學(xué)院，浙江 杭州 310000；2.吉林大學(xué) 超塑性與塑性研究所，吉林 長春 130022）

0 引 言

無縫線路是將標(biāo)準(zhǔn)長度的鋼軌焊接成長軌條并鋪設(shè)到線路上，當(dāng)環(huán)境溫度發(fā)生變化時，由于軌枕等附屬設(shè)施的存在，使得焊接長軌條不能進(jìn)行自由伸縮，鋼軌內(nèi)部會產(chǎn)生巨大的溫度應(yīng)力，同時會破壞軌道結(jié)構(gòu)。

無縫線路穩(wěn)定性的研究始于德國。1902年德國科學(xué)家哈爾曼（A.Harrmann）首次討論了無縫線路臌曲的可能性。我國一直十分重視無縫線路穩(wěn)定性的理論研究工作。1977年以鐵道科學(xué)研究院和長沙鐵道學(xué)院為主的科研小組，在總結(jié)以往研究成果的基礎(chǔ)上，提出了“統(tǒng)一無縫線路穩(wěn)定性計算公式”[1-3]。羅雁云等[4-6]通過建立無縫線路軌道脹臌曲理論模型，分析無縫線路脹軌時的位移變化規(guī)律，研究溫度應(yīng)力作用下無縫線路軌道臌曲的變化特征以及軌道參數(shù)對其的影響。石現(xiàn)峰等[7]利用傳熱學(xué)的基本理論，采用有限元分析軟件ABAQUS對板式無砟軌道結(jié)構(gòu)在溫度作用下的影響進(jìn)行仿真計算，分析不同支撐形式及不同軌道板寬度和厚度對無砟軌道結(jié)構(gòu)溫度效應(yīng)的影響。此外，國內(nèi)外學(xué)者在分析無縫線路穩(wěn)定性方面做了大量的工作，并取得了一定的成果[8-12]。

本研究以無縫線路臌曲理論為依據(jù)，結(jié)合有限元方法，利用大型有限元分析軟件ABAQUS對無縫線路進(jìn)行溫度應(yīng)力的分析研究。

1 有限元分析

1.1 鋼軌模型的建立

本研究參考的無砟軌道模型如圖1所示[13]。

圖1 I型板式無砟軌道結(jié)構(gòu)示意圖

由于ABAQUS的建模功能并不強大，筆者使用三維建模軟件CATIA建立模型，再將模型導(dǎo)入ABAQUS。鋼軌的模型選擇60 kg/m軌型。

無縫線路模型的參數(shù)取值如表1[14-15]所示。本研究主要針對鋼軌內(nèi)部溫度應(yīng)力的變化，軌道板及其底座的物理參數(shù)數(shù)非常接近，故在模型中將軌道板和底座簡化在一個模型中。由于每組扣件都對鋼軌的垂向、縱向以及橫向運動限制，筆者將每個扣件簡化為3個彈簧單元。鋼軌與軌道板的連接主要通過3個組模擬扣件的彈簧連接。

表1 無縫線路模型力學(xué)性能和材料參數(shù)

無縫線路模型如圖2所示。

圖2 初始無縫線路模型

1.2 有限元分析過程

結(jié)合文獻(xiàn)[1]可得到鋼軌溫度變化數(shù)值，分析過程主要分為3步：

（1）建立鋼軌模型的初始邊界條件。初始邊界條件分為兩個部分，其一是初始溫度場，即鋼軌的鎖定軌溫，這里取22℃。

（2）模擬氣溫上升到最高氣溫。最高軌溫取值62.7℃，模擬鋼軌軌溫上升到相應(yīng)的溫度，并輸出最高軌溫時鋼軌的溫度應(yīng)力截圖。

（3）模擬氣溫從最高點下降至最低氣溫-27℃，鋼軌軌溫也降至相應(yīng)的最低點，并輸出鋼軌模型的溫度應(yīng)力。

1.3 有限元分析結(jié)果

本研究采用有限方法對得到的結(jié)果進(jìn)行分析[16-17]，取模型中間的中心的一截鋼軌，鋼軌與扣件接觸面位于底面正中間，在軌溫最高點的溫度應(yīng)力模型如圖3所示。從圖3中可以看出，溫度應(yīng)力最大處主要集中在扣件所在位置，應(yīng)力集中區(qū)域沿著扣件逐漸減小。

圖3 最高軌溫下無縫線路的mises應(yīng)力圖

由于mises應(yīng)力適用于第4強度理論，而無縫線路中的鋼軌材料為脆性材料，多用第1強度理論，即只考察最大主應(yīng)力，本研究繼續(xù)對各向應(yīng)力分量進(jìn)行討論。

由于扣件、軌道板等基礎(chǔ)設(shè)施的約束，溫度變化時鋼軌將不能自由伸縮。因而無縫線路中最危險的就是鋼軌縱向溫度應(yīng)力，縱向溫度應(yīng)力分布如圖4、圖5所示。鋼軌縱向溫度應(yīng)力最大點位于扣件與鋼軌接觸面積邊緣，以最高軌溫為例進(jìn)行分析（如圖4所示）：由于扣件的直接約束，使得鋼軌與扣件接觸面不能自由膨脹，在接觸面受到壓應(yīng)力時，最大值為164.9 MPa，同時導(dǎo)致接觸面左右兩端的鋼軌底面受到過大的縱向拉應(yīng)力，最大值為131.7 MPa。同理可得圖5中縱向溫度應(yīng)力最大壓應(yīng)力為158.5 MPa，拉應(yīng)力為198.5 MPa。由于縱向鋼軌溫度應(yīng)力對無縫線路影響最大。查文獻(xiàn)[1]得無縫線路鋼軌強度為457 MPa。經(jīng)比較得，縱向溫度應(yīng)力滿足鋼軌強度。

無縫線路鋼軌切向溫度應(yīng)力分布圖如圖6、圖7所示。以最高軌溫為例（如圖6所示），軌溫升高時，由于扣件的約束，鋼軌與扣件接觸面積處為應(yīng)力最大區(qū)域，最大值壓應(yīng)力為165.1 MPa。圖7中最大拉應(yīng)力為198.7 MPa。

無縫線路中鋼軌的垂向溫度應(yīng)力分布如圖8、圖9所示。以圖8為例，軌溫升高時，應(yīng)力集中區(qū)域還是位于扣件與鋼軌接觸面上，沿著接觸面積的輪廓分布，最大溫度應(yīng)力值為219.4 MPa。在圖9中，軌溫降低，溫度應(yīng)力最大值為264.2 MPa。

2 結(jié)束語

圖4 最高軌溫時鋼軌縱向溫度應(yīng)力分布圖

圖5 最低軌溫時鋼軌縱向溫度應(yīng)力分布圖

圖6 最高軌溫時鋼軌切向溫度應(yīng)力分布圖

圖7 最低軌溫時鋼軌切向溫度應(yīng)力分布圖

圖8 最高軌溫時鋼軌垂向溫度應(yīng)力分布圖

圖9 最低軌溫時鋼軌垂向溫度應(yīng)力分布圖

本研究以無砟軌道為理論模型，考慮了扣件、軌道板等附屬設(shè)施的約束，并按照對應(yīng)的材料參數(shù)，完成了無縫線路的有限元建模，在軌溫升降的基礎(chǔ)上，得到了相應(yīng)的溫度應(yīng)力分布模型圖。通過進(jìn)一步分析mises應(yīng)力、縱向、切向和垂向的溫度應(yīng)力分布圖，得到了鋼軌在軌溫變化時鋼軌的受壓以及受拉應(yīng)力部位，并在圖像中顯示出來。根據(jù)第一強度理論得出：溫度變化時鋼軌內(nèi)部最大應(yīng)力為縱向溫度應(yīng)力，并且經(jīng)過比較，該應(yīng)力滿足鋼軌強度。

筆者的研究工作為進(jìn)一步研究無縫線路溫度應(yīng)力打下了基礎(chǔ)，提供了一定的思路。但在模型上還具有一定的局限性，以后的研究工作將在此基礎(chǔ)上繼續(xù)完善軌道模型，不斷深入研究。

（References）：

[1] 廣鐘巖，高慧安.鐵路無縫線路[M].4版.北京：中國鐵道出版社，2005.

[2] KISH A，SAMAVEDAM G，JEONG D.Analysis of Thermal Bucking Tests on U.S.Railroads[R].Technque Report，NTISPB83-203554，1982：21-102.

[3] JEONG D，SAMAVEDAM C，KISH A.Determination of Track Lateral Resistance from Lateral Pull Tests[R].Final Report，NTIS-PB86-225612，1986：1-49.

[4] 金壽延，朱劍月，羅雁云，等.無縫線路軌道的溫度力與振動特性關(guān)系的研究[J].機械強度，2002，24（1）：144-147.

[5] 羅雁云，施董燕，譚曉春.縱向力作用下無縫線路動態(tài)特性有限元分析[J].力學(xué)季刊，2008，29（2）：284-290.

[6] 羅雁云，李振廷.軌道參數(shù)變化對無縫線路穩(wěn)定性影響[J].中國鐵道科學(xué)，2008，29（2）：34-36.

[7] 石現(xiàn)峰，李建斌.溫度對板式無砟軌道結(jié)構(gòu)的影響研究[J].鐵道工程學(xué)報，2008（5）：30-32.

[8] DAXUN T.Distribution of Temperature Stresses along the Continuouslv Welded Track[R].American Railroad Engi?neering AssociaLion，Bulletin 692，1983：35-51.

[9] 張 斌，雷曉燕.基于車輛-軌道單元的無砟軌道動力特性有限元分析[J].鐵道學(xué)報，2011，33（7）：78-85.

[10] 雷曉燕.軌道力學(xué)與工程新方法[M].北京：中國鐵道出版社，2002.

[11] 童大塤.鐵路軌道[M].3版.北京：中國鐵道出版社，1996.

[12] 鐘方千，周小林，褚衛(wèi)松.無砟軌道無縫道岔溫度縱向力研究[J].山西建筑，2008，34（18）：24-25.

[13] 李中華.CRTSⅠ型與CRTSⅡ型板式無砟軌道結(jié)構(gòu)特點分析[J].華東交通大學(xué)學(xué)報，2010，27（1）：22-28.

[14] 李 偉，溫澤峰，金學(xué)松，等.滾滑接觸下鋼軌熱力耦合分析[J].工程力學(xué)，2010，27（8）：199-205.

[15] 李春霞，殷明旻，車曉娟.路基上無砟軌道機床反力系數(shù)取值的探討[J].鐵道建筑，2009，27（8）：199-205.

[16] 陳樂平，單 磊，陳 罡.基于Moldflow和ANSYS的防眩板有限元分析[J].輕工機械,2012,30(1):5-9.

[17] 石亦平，周玉蓉.ABAQUS有限元分析實例講解[M].北京：機械工業(yè)出版社，2006.