康維新， 陳 帥， 魏春城， 劉學(xué)毅， 李佳莉， 劉笑凱

(1. 西南交通大學(xué) 高速鐵路線路工程教育部重點(diǎn)實(shí)驗(yàn)室， 四川 成都 610031： 2. 中設(shè)設(shè)計(jì)集團(tuán)股份有限公司， 江蘇 南京 210001)

無砟軌道以其高穩(wěn)定性、高平順性和少維修等特點(diǎn)，在高速鐵路上得到廣泛運(yùn)用[1]。作為長期暴露在自然環(huán)境中的長條狀結(jié)構(gòu)物，氣象條件的年周期波動(dòng)及日周期波動(dòng)都使得軌道結(jié)構(gòu)溫度周期性波動(dòng)[2]?，F(xiàn)有針對無砟軌道溫度場的研究主要集中在正常天氣下的軌道溫度特性，對于持續(xù)高溫天氣作用下的軌道溫度尚未有人研究。已有研究成果可以按照研究手段分為試驗(yàn)和理論兩類，其中試驗(yàn)測試是獲取軌道結(jié)構(gòu)溫度特性的主要手段，楊榮山等[3]對成都地區(qū)CRTSⅠ型雙塊式軌道冬季溫度特性進(jìn)行了實(shí)測，給出了成都地區(qū)冬季無砟軌道溫度荷載的取值。戴公連等[4]對橋上無砟軌道溫度梯度進(jìn)行了實(shí)測研究，發(fā)現(xiàn)軌道板表面平均日溫差可達(dá)19 ℃，軌道板溫度沿深度呈指數(shù)分布。趙坪銳等[5]監(jiān)測了達(dá)坂城地區(qū)雙塊式軌道溫度變化規(guī)律，并依據(jù)試驗(yàn)分析給出了無砟軌道溫度取值方法。劉學(xué)毅等[6]、李佳莉[7]采用頻譜分析方法研究了成都地區(qū)實(shí)測雙塊式無砟軌道溫度與氣象因素間的關(guān)系，給出了無砟軌道溫度簡便計(jì)算方法?？稻S新[8]實(shí)測了上海氣溫高達(dá)42 ℃時(shí)無砟軌道的溫度特性，并針對反射隔熱涂料的使用效果進(jìn)行了研究。已有這些試驗(yàn)研究中都是在正常天氣下或短暫高溫天氣下得到的結(jié)論，而未能證明這些試驗(yàn)得到的溫度荷載取值方法、經(jīng)驗(yàn)公式等在持續(xù)高溫天氣下的適用性。

目前對于正常天氣下無砟軌道溫度場的理論計(jì)算和數(shù)值仿真較少。歐祖敏等[9]推導(dǎo)了無砟軌道溫度場日變化過程的計(jì)算公式，分析了各氣象因素對無砟軌道溫度的影響。閆斌等[10]建立了有限元模型，分析了軌道溫度荷載模式，推薦了部分地區(qū)溫度荷載值。這些研究為無砟軌道溫度荷載取值提供了依據(jù)，但這些研究中建立的模型及其成果并不能直接用來研究持續(xù)高溫天氣的影響。在已有理論推導(dǎo)中，未考慮地溫對軌道溫度的影響，而是將其完全簡化為日周期問題，這顯然無法體現(xiàn)持續(xù)高溫天氣的影響；現(xiàn)有軌道溫度仿真模型中，將軌道下部溫度邊界條件假設(shè)為隔熱條件或某一實(shí)測的恒定溫度，這種模型在計(jì)算短期溫度變化是可行的，但計(jì)算時(shí)間較長時(shí)下部邊界條件會(huì)發(fā)生變化，此時(shí)這種假設(shè)會(huì)帶來較大誤差，且軌道越深處誤差越大，時(shí)間越長誤差也越大。高速鐵路設(shè)計(jì)規(guī)范[11]以現(xiàn)有研究成果為基礎(chǔ)，將無砟軌道溫度荷載分為兩類，即年溫差和溫度梯度，但規(guī)范對于年溫差的計(jì)算方法并未明確給出，對于溫度梯度則簡單地取正溫度梯度90 ℃/m，負(fù)溫度梯度-45 ℃/m，未能考慮持續(xù)高溫天氣的影響。為此本文根據(jù)實(shí)測氣象資料，基于傳熱學(xué)及氣象學(xué)原理，分析自然環(huán)境中無砟軌道的換熱機(jī)理，推導(dǎo)適用于持續(xù)高溫天氣下軌道溫度計(jì)算公式，開展試驗(yàn)研究，分析持續(xù)高溫下無砟軌道溫度特性。

1 無砟軌道結(jié)構(gòu)溫度場模型

1.1 無砟軌道結(jié)構(gòu)溫度場模型的簡化

現(xiàn)有無砟軌道主要以板式軌道和雙塊式軌道為主，板式軌道由預(yù)制軌道板、CA砂漿或自密實(shí)混凝土層及支承層(底座板)組成，雙塊式軌道由現(xiàn)澆道床、預(yù)制軌枕及支承層(底座板)組成[12-14]。作為組合式結(jié)構(gòu)，無砟軌道溫度場實(shí)質(zhì)為復(fù)雜邊界條件下的層狀固體三維熱傳導(dǎo)問題。

為了研究無砟軌道溫度場的主要規(guī)律，需要進(jìn)行適當(dāng)簡化。根據(jù)文獻(xiàn)[15-16]混凝土、砂漿及各類土等材料熱工參數(shù)基本一致的研究結(jié)果，不考慮鋼軌和扣件影響，可認(rèn)為軌道結(jié)構(gòu)是熱傳導(dǎo)各向同性且均一的。無砟軌道作為長條狀結(jié)構(gòu)物，在一定長度內(nèi)可認(rèn)為其溫度分布沿線路方向差別不大。此外，考慮到無砟軌道側(cè)邊與空氣接觸面積小，交換的熱量少，僅對軌道結(jié)構(gòu)的邊角溫度有影響，且無砟軌道材料熱工參數(shù)較為穩(wěn)定，隨溫度變化小，此時(shí)可將復(fù)雜的三維問題簡化為垂向一維線性熱傳導(dǎo)問題。

無砟軌道上表面受氣象變化影響，下表面(支承層或底座板下表面)受基礎(chǔ)溫度變化的影響。氣象因素隨時(shí)間變化規(guī)律復(fù)雜，難以用簡單函數(shù)表示，這是無砟軌道溫度場計(jì)算的一大難點(diǎn)。無砟軌道下表面溫度除實(shí)測外，難以直接獲取。由于無法給出下表面邊界條件，故數(shù)值方法難以根據(jù)實(shí)測氣象資料全面計(jì)算無砟軌道溫度場，這也是軌道結(jié)構(gòu)溫度計(jì)算和預(yù)估的另一大難點(diǎn)，為此本文將無砟軌道簡化見圖1，只有上表面的半無限空間體，以解決下表面溫度未知的困難。

1.2 無砟軌道與環(huán)境的熱量交換

1.2.1 太陽輻射

太陽輻射是無砟軌道的主要熱量來源，到達(dá)地球的太陽輻射強(qiáng)度相對恒定，波動(dòng)很小，但無砟軌道表面能夠接受到的太陽輻射受大氣、地理位置等因素影響很大。進(jìn)入大氣的太陽輻射一部分被云層反射回宇宙空間(見圖1)，一部分被臭氧、水蒸氣等吸收，一部分被空氣分子及塵埃向各個(gè)方向散射，剩余的部分(太陽直接輻射)以一定的角度照射到地表上，即到達(dá)地表的太陽輻射為天空散射的一部分和太陽直接輻射的總和?；炷敛⒉皇墙^對黑體，照射在其上的太陽輻射有些會(huì)被反射掉。混凝土的太陽輻射吸收系數(shù)受其表面亮度、平整度及臟污程度等的影響程度較大，取值為0.5～0.7[17-21]。研究表明[22-23]：混凝土新澆筑時(shí)取0.5；使用時(shí)間較長時(shí)取0.7。太陽輻射帶給無砟軌道的熱流密度為

qs=asQ

( 1 )

式中：qs為太陽輻射折算后的熱流密度，W/m2；as為軌道表面的太陽輻射吸收系數(shù)；Q為單位面積上的太陽輻射，W/m2。

1.2.2 大氣

對流熱交換和紅外電磁輻射換熱是物體與空氣熱交換的主要手段。

(1) 對流換熱及風(fēng)的影響

由冷卻定律知，對流形式的熱流密度為

qc=uc(Ta-Ts)

( 2 )

式中：qc為對流熱流密度，W/m2；Ta為氣溫，℃；Ts為軌道表面溫度，℃；uc為對流系數(shù)，W/(m2·℃)。

對流系數(shù)受風(fēng)速影響很大，算式為[24]

( 3 )

式中：v為日均風(fēng)速，m/s。

(2) 輻射換熱

現(xiàn)代物理學(xué)表明物體內(nèi)粒子的熱運(yùn)動(dòng)運(yùn)動(dòng)導(dǎo)致表面在持續(xù)的輻射電磁波，溫度越高輻射能力越強(qiáng)，波長越短。大氣溫度和無砟軌道表面溫度遠(yuǎn)低于太陽溫度，它們發(fā)出的輻射波長大于太陽輻射的波長，稱為長波輻射。大氣向各個(gè)方向發(fā)射長波輻射，朝向地表的部分稱為大氣逆輻射,見圖1。大氣逆輻射與無砟軌道表面間的差值稱為有效輻射。由Stefan-Boltzmann式[25]，以有效輻射形式進(jìn)入無砟軌道的熱流密度為

ql=εaσ(Ta+273.15)4-εsσ(Ts+273.15)4

( 4 )

式中：ql為有效輻射熱流密度，W/m2；σ為5.67×10-8，W/(m2·℃4)；εs為黑度系數(shù)，本文取0.88；εa為大氣黑度系數(shù)，受濕度、大氣溫度及云量等影響，取值范圍在0.8～0.9，為簡化計(jì)算，本文取為0.88。

利用平方差公式，將式( 4 )展開后寫成對流換熱的形式，即

ql=ul(Ta-Ts)

( 5 )

( 6 )

大氣帶給無砟軌道的總熱流密度為

qz=uz(Ta-Ts)

( 7 )

uz=

( 8 )

式中：qz為由大氣與軌道間的總熱流密度，W/m2；uz為總熱交換系數(shù)，W/(m2·℃)。

綜上所述，自然環(huán)境下無砟軌道表面的熱流密度可表達(dá)為

q=qs+qz=

( 9 )

2 無砟軌道溫度場的求解及修正

2.1 氣溫和太陽輻射的函數(shù)表示

以往研究中多采用半波正弦來表示完全晴天時(shí)到達(dá)地表的太陽輻射日變化過程，對比實(shí)測資料發(fā)現(xiàn)，即便是晴天，太陽輻射日變化與半波正弦也有較大差別，而有云天、陰天則差別更大，顯然采用這種簡單函數(shù)表示的太陽輻射日過程無法準(zhǔn)確分析無砟軌道溫度場。與太陽輻射日變化類似，以往研究中多采用正弦函數(shù)模擬氣溫日變化過程，這種方法模擬白天氣溫變化較為準(zhǔn)確，而夜間誤差較大。事實(shí)上大氣溫度受很多因素影響，其日變化規(guī)律非常復(fù)雜，晴天與陰天不同，不同地區(qū)不同，不同季節(jié)也不相同。由于氣象資料都可通過各地氣象站方便獲取，假設(shè)ci和di(i=0,1,2,3,…,24)分別表示τ=0,1,2,3,…,24時(shí)測得的氣溫和太陽輻射強(qiáng)度，則采用傅里葉級(jí)數(shù)可建立太陽輻射和氣溫表達(dá)式

(10)

(11)

(12)

(13)

式中：Q為太陽輻射強(qiáng)度，W/m2；Ta為氣溫，℃；an、bn分別為太陽輻射、氣溫的傅里葉級(jí)數(shù)系數(shù)；n為無窮級(jí)數(shù)項(xiàng)數(shù)；τ為時(shí)間，h。

2.2 無砟軌道溫度場的求解

若記無砟軌道各深度處溫度為T，則其應(yīng)滿足[25]

(14)

(15)

x=T≠

(16)

求解式(14)，得軌道溫度

(17)

(18)

(19)

(20)

(21)

式中：Sn為軌道表面溫度幅值與等效氣溫的比值；βn為軌下結(jié)構(gòu)深度x處溫度波動(dòng)幅度相對表面溫度波動(dòng)幅度的比值；φxn為x深度處溫度相對表面溫度的滯后相位；φsn為表面溫度相對等效氣溫的滯后相位。當(dāng)n取不同值時(shí)，得到βn與x的關(guān)系，見圖2。

當(dāng)x≥0.85 m，n≥1時(shí)，βn=0，式(17)變?yōu)?/p>

Tx≥0.85=A0

(22)

式中：Tx≥0.85為軌下結(jié)構(gòu)深度0.85 m及以下的溫度，℃；A0為等效氣溫的平均值，℃。

式(22)表明，軌下結(jié)構(gòu)深度0.85 m及以下處溫度與外界因素日變化無關(guān)，稱此處溫度為地溫。實(shí)際上，地溫由氣象因素年變化決定。

2.3 地溫修正

除地球自轉(zhuǎn)導(dǎo)致氣象因素日周期變化外，地球繞太陽公轉(zhuǎn)也會(huì)使得氣象因素產(chǎn)生年周期變化，年變化對無砟軌道表面溫度影響不大，但隨著深度增加，氣象年變化逐漸取代日變化，成為決定無砟軌道溫度的主要因素。

由前述分析知，地溫與氣象因素的日變化無關(guān)，由年變化決定。若已知365 d的逐日平均氣溫和平均太陽輻射，應(yīng)用前述類似方法可求得軌下結(jié)構(gòu)深度0.85 m處溫度為

(23)

(24)

(25)

式中：Syn、φsy與Sn和φsn類似，分別為軌道表面溫度年變化幅值相對等效氣溫年變化幅值的比值和滯后相位；Ayn為年等效氣溫的振幅，℃；αy為混凝土以日為單位的導(dǎo)溫系數(shù)，0.072 m2/d；uzy為以年平均風(fēng)速和年平均氣溫代入式( 8 )算得的總換熱系數(shù)，W/(m2·℃)；τy為時(shí)間，d。

(26)

式中：Tc為某一深度處溫度修正量，℃。

故無砟軌道溫度場最終解為

(27)

而各深度處的日平均溫度可由等效氣溫日平均值與地溫線性插值得到

(28)

3 無砟軌道溫度計(jì)算公式的驗(yàn)證

為研究無砟軌道溫度場的特性，修筑了長10 m的CRTSⅠ型雙塊式軌道，軌道走向?yàn)楸逼?7°, 見圖3(a)。試驗(yàn)地點(diǎn)位于成都(104°1′E，30°5′N)，海拔500 m，常年風(fēng)速較低，氣候濕潤，夏季日照較多，冬季日照較少。按圖3(b)所示位置在道床板上下表面及距離上表面0.16 m深度處布置鉑電阻溫度傳感器監(jiān)測軌道溫度，并在軌道上安裝了氣象站，對氣象參數(shù)進(jìn)行監(jiān)測。

2015年8月21日—2016年8月21日的日平均太陽輻射強(qiáng)度、日平均氣溫見圖4(a)，該段時(shí)間平均風(fēng)速約為1 m/s。2016年8月21日太陽輻射、氣溫見圖4(b)，當(dāng)日平均風(fēng)速0.3 m/s。由于試驗(yàn)是在該段無砟軌道修筑好兩個(gè)月后開始監(jiān)測溫度，此時(shí)混凝土尚新，太陽輻射吸收系數(shù)取0.5。

利用式(23)、式(28)分別算得2015年8月21日—2016年8月21日的地溫和道床板中部溫度，見圖5(a)。式(28)算得試驗(yàn)期間道床板表面和板底面日平均溫度與實(shí)測值的對比見圖5(b)。

圖5(a)表明，地溫變化趨勢與日平均氣溫變化趨勢一致，但變化幅值有所減小，道床板內(nèi)部溫度由于受太陽輻射的影響，日平均溫度高于氣溫，年變化幅度約為40 ℃，故成都地區(qū)無砟軌道年溫差可取40 ℃。圖5(b)表明，道床板向陽處的日平均溫度最高，板中溫度次之，背陽處日平均溫度最小，但變化趨勢一致。無砟軌道上下表面日平均溫度計(jì)算值與實(shí)測值吻合較好，證實(shí)了本文提出的無砟軌道各深度日平均溫度計(jì)算公式的準(zhǔn)確性。

應(yīng)用式( 8 )算得當(dāng)日的總換熱系數(shù)uz為13.1 W/(m2·℃)，算得等效氣溫見圖4(b)虛線，最大值達(dá)60 ℃。8月21日地溫為30 ℃，應(yīng)用式(27)可計(jì)算道床板上下表面及深度0.16 m處的溫度，將其與道床板相應(yīng)測點(diǎn)的實(shí)測溫度對比，見圖6。

圖6表明向陽側(cè)(A處)道床板溫度整體高于背陽側(cè)(C處)溫度約3 ℃，但板邊處溫度變化趨勢與板中心處(B處)溫度變化趨勢一致。由圖6可知，道床板中心上下表面、道床板中部溫度的理論計(jì)算值與實(shí)測值吻合較好，說明本文推導(dǎo)的公式能夠有效計(jì)算無砟軌道溫度場。

根據(jù)式(27)算得當(dāng)日不同時(shí)刻軌道溫度分布見圖7。由圖7可得，軌道表面溫度波動(dòng)18 ℃， 0.2 m處6 ℃，下表面約4 ℃。無砟軌道在0～0.2 m深度內(nèi)日溫度變化較大，容易出現(xiàn)較大溫度梯度。凌晨0:00到8:00，道床板下表面溫度較高，道床板呈負(fù)溫度梯度，支承正溫度梯度。6:00時(shí)出現(xiàn)負(fù)溫度梯度極值-15 ℃/m。8:00時(shí)受太陽影響，軌道溫度快速上升，16:00時(shí)正溫度梯度達(dá)到極值，為70 ℃/m。16:00至次日凌晨0:00，道床板表面溫度逐漸降低，熱量集聚在道床板內(nèi)，形成道床板中心溫度高，上下表面溫度低的非線性分布。

4 持續(xù)高溫天氣下的無砟軌道溫度場

氣象學(xué)根據(jù)人體適應(yīng)性，將高溫天氣分有藍(lán)、黃、橙、紅級(jí)別，最高氣溫連續(xù)三天超過35 ℃時(shí)稱為黃色預(yù)警。2013年夏季長江下游地區(qū)連續(xù)50多天最高氣溫超過35 ℃，部分無砟軌道出現(xiàn)上拱等病害，給行車帶來巨大安全隱患。

杭州2013年日平均氣溫和太陽輻射見圖8(a)，最熱當(dāng)天的逐小時(shí)太陽輻射強(qiáng)度、氣溫和等效氣溫見圖8(b)，杭州當(dāng)年平均風(fēng)速1.5 m/s，持續(xù)高溫期間風(fēng)速0.3 m/s，太陽輻射吸收系數(shù)取0.5。

圖8(a)表明，杭州地區(qū)從2013年第180 d到第230 d，日平均太陽輻射強(qiáng)度在400 W/m2左右，日平均氣溫約為35 ℃，是典型的持續(xù)時(shí)間很長的高溫天氣，圖8(b)表明最熱一天最高氣溫超過40 ℃，正午太陽輻射達(dá)到1 100 W/m2。

利用氣象資料，以雙塊式軌道為例算得杭州地區(qū)2013年地溫及道床板中部溫度見圖9(a)，從第180 d開始，地溫隨日平均氣溫的上升而持續(xù)上升，最熱當(dāng)日地溫達(dá)到42 ℃，道床板中最高日平均溫度更是達(dá)到47.5 ℃，遠(yuǎn)高于日平均氣溫。最熱當(dāng)日各深度處溫度變化見圖9(b)，無砟軌道各深度處溫度遠(yuǎn)高于氣溫，表面最高溫度達(dá)到65 ℃，支承層溫度達(dá)到46 ℃。道床板表面日溫差達(dá)25 ℃，道床板中部日溫差12.5 ℃，支承層日溫度變化較小，溫差較小。

利用式(27)可算得最熱當(dāng)日各深度處溫度分布見圖10。相比圖7，持續(xù)高溫下的無砟軌道溫度分布主要特點(diǎn)是：下部基礎(chǔ)溫度高，無砟軌道整體溫度偏高。在杭州持續(xù)高溫下，軌下結(jié)構(gòu)深度0.2 m范圍內(nèi)，日最大正溫度梯度達(dá)到100 ℃/m，高于TB 10621—2014《高速鐵路設(shè)計(jì)規(guī)范》[11]規(guī)定值，日最大負(fù)溫度梯度達(dá)到-30 ℃/m。

對2013年持續(xù)高溫下無砟軌道溫度場的計(jì)算表明：持續(xù)高溫期間天氣持續(xù)晴朗，太陽輻射持續(xù)較強(qiáng)，平均氣溫持續(xù)上升，無砟軌道整體溫度也隨之上升，由于太陽輻射的影響，無砟軌道整體溫度遠(yuǎn)高于氣溫。道床板中部溫度高于氣溫15 ℃，最高可達(dá)55 ℃，日均溫度為46 ℃。在持續(xù)高溫天氣下，無砟軌道除了整體溫度升高外，其正溫度梯度也高于規(guī)范規(guī)定90 ℃/m。

由上述分析可以看出，氣象學(xué)上對于持續(xù)高溫的定義較難適用于無砟軌道溫度。對于無砟軌道，不僅要關(guān)注最高氣溫，而且要關(guān)注日平均氣溫，當(dāng)出現(xiàn)高溫天氣時(shí)，若日平均氣溫也較高時(shí)就要注意無砟軌道的整體溫度，預(yù)防道床板/軌道板上拱等病害的發(fā)生。

5 結(jié)論

本文在分析無砟軌道與大氣環(huán)境換熱機(jī)理的基礎(chǔ)上，考慮了氣象日變化、年變化及地溫對軌道溫度的影響，推導(dǎo)了基于氣溫、風(fēng)速及太陽輻射的無砟軌道溫度場計(jì)算公式，開展實(shí)測驗(yàn)證，計(jì)算分析了持續(xù)高溫天氣下無砟軌道的溫度特性。

主要結(jié)論為：

(1) 忽略鋼軌扣件等的影響，將無砟軌道簡化為受周期性熱邊界條件影響的導(dǎo)熱均質(zhì)半無限空間體，能夠較為準(zhǔn)確有效地計(jì)算無砟軌道溫度特性。

(2) 無砟軌道溫度不僅受氣象日變化影響，而且受氣象年變化的影響。氣象日變化主要影響深度小于0.85 m的溫度，氣象年變化主要決定深度大于0.85 m的溫度(即地溫)，計(jì)算無砟軌道溫度時(shí)，應(yīng)考慮氣象年變化造成的地溫修正。

(3) 持續(xù)高溫天氣時(shí)，日平均氣溫很高，導(dǎo)致無砟軌道下部基礎(chǔ)溫度(地溫)較高，地溫能夠達(dá)到42 ℃，熱量難以向下傳遞。無砟軌道整體溫度升高且遠(yuǎn)高于氣溫，可比氣溫高約15 ℃，而且出現(xiàn)較大的溫度梯度，最大正溫度梯度能夠達(dá)到100 ℃/m。