基于離散元法的散體道床與基床表層耦合變形機理研究

高 亮， 殷 浩， 徐 旸， 楊國濤,2

(1.北京交通大學 土木建筑工程學院，北京 100044；2.中國鐵路總公司，北京 100844)

有砟軌道具有價格低廉、便于養(yǎng)護維修、振動噪聲小的優(yōu)點。但由于散體道床主要由碎石道砟組成，其上部承受了經(jīng)軌枕傳遞而來的列車荷載，而下部又與軟質(zhì)路基相接觸，在循環(huán)荷載作用下會發(fā)生散體道床的累積沉降及道砟顆粒的破碎、粉化等現(xiàn)象。如養(yǎng)護維修不善，就會在水的作用下形成道床翻漿冒泥、道砟囊等病害，使有砟軌道彈性逐漸喪失，導致線路幾何不平順的產(chǎn)生，嚴重時會威脅高鐵列車的安全運營。國內(nèi)外學者對這一問題開展了大量研究工作。

在試驗研究方面，佐藤吉彥等[1]提出線路變形由道床下沉和路基下沉導致，并結(jié)合現(xiàn)場試驗、運營經(jīng)驗分別總結(jié)了道床和路基的下沉公式，但并未對兩者下沉關(guān)系進行深入研究；平野雅之等[2]對散體道床下沉的室內(nèi)試驗和現(xiàn)場數(shù)據(jù)進行統(tǒng)計，估算得到了行車速度為70 km/h條件下道床與路基下沉量比值為0.12～3.13；美國進行FAST軌道試驗[3]研究重復荷載作用下軌面總下沉中道砟層、墊砟層和路基的下沉比例，指出軌道總下沉中的一半以上來自道床；胡仁偉等[4-5]通過室內(nèi)道砟動三軸試驗，對循環(huán)荷載作用下散體道床的累積變形率以及道床彈性隨荷載作用次數(shù)的變化規(guī)律進行了研究，但并未考慮下部路基基礎(chǔ)對道床沉降的影響。

僅采用試驗的手段無法獲得散體道床與路基結(jié)構(gòu)相互作用的細觀機理，國內(nèi)外學者還采用數(shù)值仿真和機理分析的方法。Mohammad等[6]結(jié)合有限元和邊界元法研究道床、路基、下臥土的力學性能；Ferreira等[7]建立車輛、軌道、路基系統(tǒng)有限元模型，分析軌道沉降規(guī)律；高建敏等[8]利用車輛-軌道耦合動力學模型，分析運營條件和軌道結(jié)構(gòu)參數(shù)對道床下沉變形的影響；徐鵬等[9]采用連續(xù)體建模方法，建立列車-有砟軌道-路基空間耦合模型，并采用Galerkin離散法研究道床、路基的變形特性；徐旸等[10]基于離散元法對道砟進行精細化建模，通過建立道砟箱模型研究道床累積沉降規(guī)律；付龍龍等[11]采用二維顆粒流方法建立軌道離散元模型，分析路基不均勻沉降及列車荷載變化對軌面沉降的影響；王平等[12]利用散粒體介質(zhì)理論中的機理分析法，建立理想碎石道床的平均增量彈性本構(gòu)，分析道床彈性性能。

綜上所述，已有研究的不足之處主要包括如下幾個方面：① 已有試驗研究多是針對道砟集料所開展的模型試驗，難以從結(jié)構(gòu)角度對散體道床與路基的耦合作用進行研究，且受實驗條件的限制，難以考慮高速運營的行車條件；② 數(shù)值仿真方面多是采用單一的離散元法僅對散粒體碎石道床進行仿真，或是采用有限元法將散體道床視為連續(xù)介質(zhì)層，對有砟軌道整體結(jié)構(gòu)力學特性進行分析，均無法真實反映出散體道床與多層體系間的點-點、點-面接觸特性以及道砟顆粒之間的相互作用，與真實情況存在較大的差別；③ 軌道沉降主要由道床沉降與路基沉降兩部分組成，而已有研究多是通過模型試驗對散體道床的局部沉降力學特性進行研究，缺乏從細觀角度對考慮道床與軌排之間、道床與路基之間耦合作用下有砟軌道垂向變形機理的研究。

本文針對已有研究的不足，考慮真實道砟顆粒外形，采用離散元法建立了“軌枕-道床-基床表層耦合模型”，從軌道結(jié)構(gòu)受力、彈性變形、塑性變形等方面研究路基基礎(chǔ)上散體道床垂向變形機理。并分析了高速運營條件下，行車參數(shù)對散體道床及下部基礎(chǔ)結(jié)構(gòu)變形特性的影響規(guī)律，以期為高速鐵路散體道床的養(yǎng)護維修及狀態(tài)評估提供科學依據(jù)。

1 模型建立及驗證

1.1 模型建立

道砟顆粒的形狀會對散體道床的力學特性產(chǎn)生顯著的影響，本文基于真實道砟顆粒外形，應用PFC2D離散元軟件，采用clump單元對道砟顆粒和軌枕進行模擬。我國高速鐵路有砟軌道結(jié)構(gòu)的基床表層與道床直接接觸，且路基結(jié)構(gòu)中基床表層內(nèi)部動應力衰減最為急劇，其下結(jié)構(gòu)所承受列車荷載明顯減小，故路基結(jié)構(gòu)只模擬到基床表層部分。同時，基床表層級配碎石粒徑相對較小，結(jié)構(gòu)密實程度較高，其外形參數(shù)對結(jié)構(gòu)垂向變形特性影響較小，故根據(jù)《高速鐵路設(shè)計規(guī)范》[13]中基床表層級配碎石粒徑級配曲線要求，采用ball單元進行模擬。有砟道床厚度取0.35 m，基床表層厚度取0.7 m。

本文選用典型道砟顆粒，結(jié)合圖像識別手段，考慮計算精度和效率的最優(yōu)化，結(jié)合已有研究進行計算后確定每個clump單元由9～15個ball單元組成，再根據(jù)道砟粒徑反算ball單元半徑。基于MATLAB編寫了道砟顆粒單元球坐標矩陣的生成程序，建立了高速鐵路道砟顆粒外形參數(shù)數(shù)據(jù)庫。同時采用FISH語言編寫了模型生成子程序，實現(xiàn)了在模型生成過程中從典型顆粒形狀中對道砟顆粒進行隨機抽取。本文所選用的典型道砟顆粒及其離散元模型，如圖1所示。圖1中clump單元是不可變形的剛體結(jié)構(gòu)，只有最外層ball單元參與力學作用，結(jié)構(gòu)內(nèi)部間隙對顆粒間的相互作用沒有影響。

為減少邊界效應的影響，采用遠置邊界條件的方法選取5跨軌枕進行建模，并將邊界wall單元同與之接觸的顆粒設(shè)置為具有相同的摩擦系數(shù)。在加載前對各層結(jié)構(gòu)充分壓實，使其達到《高速鐵路設(shè)計規(guī)范》中高速列車運行所需要求。經(jīng)充分壓實，并施加列車循環(huán)荷載至沉降穩(wěn)定后的軌枕-道床-基床表層耦合模型，如圖2所示。

圖2 軌枕-道床-基床表層耦合模型

模型中各顆粒單元和wall單元均滿足剛性假設(shè)，顆粒間相互作用基于Mohr-Coulomb滑動摩擦準則。模型中clump單元和ball單元均以接觸重疊量作為接觸判定以及力-位移參數(shù)計算的媒介物理量，計算時在模型結(jié)構(gòu)之間以時步為單位進行迭代。在參數(shù)的選取上參考Lim[14]、Lu等[15]以及課題組已有相關(guān)研究，并通過改變模型各細觀參數(shù)進行大量試驗分析，最終確定模型各細觀參數(shù)，如表1所示。

表1 模型參數(shù)

1.2 荷載施加

參考道砟箱試驗和RTF(Railway Test Facility)列車荷載模擬試驗[16]，在軌枕正上方施加垂向余弦荷載，荷載公式為

F(t)=F0(cos 2πft-1)

(1)

式中：F(t)為枕上壓力；t為加載時間；荷載峰值F0取決于列車軸重，并受行車速度影響。USACE[17]通過現(xiàn)場測試和理論分析指出，每組車輪荷載主要通過5跨軌枕向下傳遞，車輪正下方軌枕承受荷載作用最大，為0.4倍車輪荷載。故列車通過時枕下最大壓強P可以通過Jeffs等[18]的彈性理論計算得到

β

(2)

式中：m為列車軸重；g為重力加速度；α為動載系數(shù)；S為軌枕底面積；β為軌道結(jié)構(gòu)影響系數(shù)，取決于軌枕類型、道床狀態(tài)等因素(本文β=1)。通過調(diào)整F0大小，保證枕下壓強的模型計算值與理論值P相一致。

加載頻率f取決于列車通過速度，由f=v/L計算得到，其中v為行車速度(m/s)，L選取時考慮最不利加載條件[19]，取轉(zhuǎn)向架軸距2.5 m。

沿線路方向模型長度為3 m，各軌枕承受的列車荷載存在時間差，故不可將列車荷載同步施加至各軌枕。圖3為列車以某速度通過時模擬列車荷載曲線示意圖，此時一組車輪位于軌枕3正上方，該軌枕已完成一次加載并處于待卸載階段，軌枕1、軌枕4處于加載階段，軌枕2、軌枕5處于卸載階段。此時軌枕3承受余弦荷載曲線峰值荷載(C點)，軌枕2、軌枕4只承受余弦荷載曲線峰值約1/5大小荷載(B點、D點)，應進一步考慮對多跨軌枕同時加載時存在的荷載作用時間差現(xiàn)象，故采用改進后的余弦荷載公式進行加載，見式(3)和式(4)

圖3 荷載作用示意圖(軌枕-道床結(jié)構(gòu))

F(t)=F0[cos 2πf(t+xTv)-1]

(3)

Tv=d/v

(4)

式中：Tv為施加在軌枕上的列車余弦荷載時間差；d為軌枕間距；v為列車通過速度，m/s；x取值由軌枕位置決定：以軌枕3為原點，水平向右為x軸正方向建立直線坐標系，相鄰軌枕間距為一個單位長度，則軌枕1～軌枕5分別取x=-2、x=-1、x=0、x=1、x=2。

1.3 模型驗證

文獻[20]中對秦沈客運專線路基基礎(chǔ)上有砟軌道結(jié)構(gòu)進行綜合測試，得到了四種機車車輛通過時軌道結(jié)構(gòu)受力與變形數(shù)據(jù)。本文應用所建模型，對秦沈線上兩列試驗車型通過時，基床表層的彈性變形量和動應力進行計算。模型加載參數(shù)與現(xiàn)場條件一致：通過式(2)計算得到兩種車型軸重對應的枕下壓強；行車速度為250 km/h和300 km/h時，對應加載頻率分別為22.78 Hz、33.33 Hz。計算對比結(jié)果如表2所示。

表2 基床表層彈性變形、動應力實測與計算結(jié)果對比

表2中彈性位移和動應力的計算值均取自模型經(jīng)加載達沉降穩(wěn)定后基床表層頂面位置，由表2可知，相同運營條件下，本文軌枕-道床-基床表層耦合模型計算結(jié)果與實測數(shù)據(jù)較為吻合，但存在微小誤差。分析原因為雖然仿真模型中采用的clump單元能較精細地模擬道砟顆粒外形，但考慮到計算效率問題，在保證計算精度的前提下對模型簡化，模型與真實道砟存在微小差異，使仿真模型中道砟顆粒間咬合力略小于真實情況，故仿真結(jié)果稍大于實測數(shù)據(jù)；模型通過設(shè)置measure圓監(jiān)測基床表層頂面應力，當對軌道結(jié)構(gòu)加載時，底層道砟會產(chǎn)生向下的位移而進入到measure圓中，使應力監(jiān)測結(jié)果偏大，故表中動應力計算誤差稍大于彈性變形量誤差。但模型誤差在可接受范圍內(nèi)，表明模型準確可靠。

2 軌道結(jié)構(gòu)垂向變形機理

本節(jié)針對行車速度250 km/h、列車軸重14 t條件下，軌道結(jié)構(gòu)受力特性和變形機理進行討論。

2.1 細觀接觸力特性

道砟顆粒之間接觸力過大是導致道砟劣化、道床沉降的直接因素，為研究軌道結(jié)構(gòu)垂向變形機理，首先需分析結(jié)構(gòu)受力特性。在列車荷載作用下，道床、基床表層結(jié)構(gòu)內(nèi)部各離散單元之間的接觸力形成力鏈，力鏈在列車荷載加載、卸載過程中逐漸演變。圖4顯示了列車荷載作用于中間軌枕3正上方時，軌道結(jié)構(gòu)接觸力分布情況，力鏈線條粗細與接觸力大小成正比。

圖4 軌道結(jié)構(gòu)接觸力分布圖

由圖4可知，道砟顆粒與軌枕底面有效接觸面較小，造成該位置處接觸力最大且分布不均勻，會加大枕下道砟劣化，加快道床沉降；接觸力在傳遞過程中逐漸減小，作用范圍逐漸擴大，相鄰軌枕下接觸力力鏈交織在一起構(gòu)成力鏈網(wǎng)絡(luò)；道床與基床表層交界面接觸力為點-點接觸狀態(tài)下的非均勻分布；由于軌道結(jié)構(gòu)的擴散作用，以及基床表層碎石粒徑相對較小，使接觸力在基床表層中的分布相比于道床中更為均勻、密集，但其大小已明顯降低。

枕下接觸力經(jīng)道床、基床表層呈錐臺狀向下擴散，但由于同一時刻不同軌枕承受的列車荷載大小不同，且接觸力擴散具有時效性，導致不同枕下接觸力發(fā)生不均勻疊加，使接觸力擴散的錐臺形狀存在“斷層”(如圖5中虛線部分)，該區(qū)域內(nèi)接觸力相對較小。軌枕2、軌枕4下接觸力會受到軌枕3上車輪荷載影響，模擬列車自右向左運行，此時軌枕4下承受的上一組車輪荷載并未擴散完全，這部分荷載與受軌枕3影響的荷載相疊加，導致圖中自左向右各軌枕下接觸力分布范圍逐漸增大。因此，在接觸力錐臺形狀斷層處形成一組斜面，向結(jié)構(gòu)右下方延伸。

圖5 接觸力傳遞規(guī)律示意圖

在道床與基床表層交界面處存在道砟與基床表層碎石相互入侵現(xiàn)象，如圖6所示。入侵量與道床厚度、道砟顆粒外形參數(shù)、基床表層壓實度等因素有關(guān)，并受線路運營條件影響。這種現(xiàn)象的存在會增加道床臟污，降低道床彈性，加速道床下沉，同時引起基床病害，應在設(shè)計工作中予以重視。

圖6 道砟與基床表層碎石相互入侵

為研究軌道結(jié)構(gòu)在列車移動荷載作用下的變形機理，對列車自右向左運行過程中軌道結(jié)構(gòu)細觀力學特性的演變規(guī)律進行分析。當車輪分別位于軌枕4、軌枕2正上方時，接觸力分布情況，如圖7(a)、圖7(b)所示。在模擬列車運營條件不變時，枕下接觸力的錐臺形狀擴散角度、錐臺狀斷層部分形成斜面的斜率均保持不變。由于枕下各結(jié)構(gòu)單元之間接觸均為點-點接觸，故接觸力存在位置已基本固定，其傳遞路徑不會隨荷載作用位置的改變發(fā)生明顯的變化；而最大接觸力存在位置會隨著荷載作用位置不同而發(fā)生變化，但始終存在于車輪荷載正下方。在接觸力分布圖中表現(xiàn)為，當荷載作用于不同軌枕時，不同枕下位置處接觸力力鏈線條粗細發(fā)生明顯變化，而力鏈位置保持不變。

2.2 彈性、塑性變形機理

有砟軌道結(jié)構(gòu)彈性性能主要由道床彈性決定，來源于相接觸道砟顆粒間的彈性壓縮和彈性滑移，與軌道結(jié)構(gòu)基本狀態(tài)有關(guān)，受列車運營條件影響。彈性變形量的大小能直接反應出軌道結(jié)構(gòu)的彈性性能和健康狀態(tài)：若彈性變形量過小，會造成鋼軌傷損，振動、噪聲增大，旅客舒適性降低，導致軌道結(jié)構(gòu)在線路運營初期發(fā)生破壞；若彈性變形量過大，會加劇道砟磨耗，導致道床幾何形狀改變，降低道床服役性能。

圖8為列車荷載作用下軌枕、基床表層頂面碎石垂向位移曲線，圖中對局部曲線進行放大。由圖8可知，軌道結(jié)構(gòu)垂向位移由彈性變形和塑性變形組成，單個加載周期內(nèi)軌道結(jié)構(gòu)的垂向位移曲線可以反應出結(jié)構(gòu)的彈、塑性變形特性：在卸載過程中可恢復的變形為結(jié)構(gòu)彈性變形，不可恢復的微小殘余變形體現(xiàn)為結(jié)構(gòu)塑性變形，在列車重復荷載作用下，這種微小塑性變形的逐漸累積就導致了軌道結(jié)構(gòu)沉降。由于單個加載周期內(nèi)結(jié)構(gòu)塑性變形量非常小，故可近似將軌道結(jié)構(gòu)單元在單個加載周期內(nèi)的垂向位移曲線視為其彈性變形曲線，反應出軌道結(jié)構(gòu)在高速列車荷載作用下垂向彈性變形特性。

(a)荷載作用于4號軌枕正上方

(b)荷載作用于2號軌枕正上方

圖8 軌枕、基床表層頂面碎石垂向位移曲線

碎石顆粒的相互錯位和重新排列是產(chǎn)生軌道結(jié)構(gòu)塑性變形的原因之一，塑性變形主要來自道床變形和路基變形，會影響列車運行安全性和舒適性，是評價道床服役狀態(tài)的重要指標。為觀察軌道結(jié)構(gòu)的塑性變形規(guī)律，將監(jiān)測到的軌枕和基床表層頂面碎石垂向位移數(shù)據(jù)進行處理：只保留監(jiān)測小球在每個加載周期內(nèi)的最大位移值，得到軌道結(jié)構(gòu)累積塑性變形曲線(即沉降曲線)，如圖9所示。由于道床沉降時，位于下部的基床表層也發(fā)生相應沉降，即道床下邊界為可動界面，故道床的沉降曲線是通過將軌枕沉降與基床表層沉降相減后處理得到。

圖9 軌枕、道床、基床表層累積塑性變形曲線

3 列車運營條件對軌道垂向變形的影響

軌道結(jié)構(gòu)垂向變形主要由列車荷載引起，與列車運營條件密切相關(guān)，為分析列車運營條件對軌道垂向變形的影響，本節(jié)研究不同行車速度、列車軸重下軌道結(jié)構(gòu)垂向變形機理。模擬列車軸重為14 t時，行車速度分別為250 km/h、300 km/h、350 km/h、400 km/h，以及行車速度為250 km/h時，列車軸重分別為14 t、16 t、18 t、20 t。

3.1 對彈性變形的影響

軌枕的彈性變形直接影響鋼軌的振動和輪軌受力，是軌道結(jié)構(gòu)彈性性能的直接體現(xiàn)，故本節(jié)著重分析軌枕的彈性變形特性。由軌道結(jié)構(gòu)垂向位移曲線(見圖8)可知，不同加載階段軌枕彈性變形量也不相同，故在加載至100次、200次、300次、400次、500次時分別統(tǒng)計附近50個彈性變形量值，分別計算其平均值作為加載至該階段時軌道結(jié)構(gòu)彈性變形量，得到不同行車速度、列車軸重條件下，位于不同加載階段時軌枕彈性變形量變化曲線，如圖10所示。

由圖10可知，軌枕彈性變形量隨著加載次數(shù)的增加逐漸減小并趨于穩(wěn)定，說明在列車荷載長期作用下，軌道結(jié)構(gòu)剛度逐漸增加，彈性性能逐漸降低。此外，軌枕彈性變形量隨著行車速度的提高、列車軸重的增大非線性增大，且車速越高、軸重越大，對其影響越大：當行車速度由350 km/h提升至400 km/h時，軌枕彈性變形量增大18.2%，是行車速度由200 km/h提升至250 km/h時軌枕彈性變形增量的3倍多；當軸重由18 t提升至20 t時，軌枕彈性變形量增大16.8%，約為軸重由14 t提升至16 t時軌枕彈性變形增量的2倍。

(a)車速影響(b)軸重影響

圖10 不同運營條件下軌枕彈性變形曲線

Fig.10 Elastic deformation curves in different operating conditions

3.2 對塑性變形的影響

圖11、圖12為不同行車速度、列車軸重條件下道床和基床表層累積塑性變形曲線，圖中不同運營條件下軌道結(jié)構(gòu)沉降表現(xiàn)出相同規(guī)律：軌道結(jié)構(gòu)沉降速率在荷載作用前期較快，后期逐漸減慢并趨于穩(wěn)定；軌枕完成前期、后期沉降所需時間均明顯長于基床表層，基床表層最先達到沉降穩(wěn)定階段。此外，道床和基床表層沉降速率、總沉降量均隨著行車速度、列車軸重的提高非線性增大。

圖11 行車速度對軌道結(jié)構(gòu)累積塑性變形影響曲線

圖12 列車軸重對軌道結(jié)構(gòu)累積塑性變形影響曲線

對加載至穩(wěn)定階段后軌道結(jié)構(gòu)沉降量進行統(tǒng)計得到表3、表4所示數(shù)據(jù)，表中“軌枕沉降”表示軌道結(jié)構(gòu)整體沉降，由道床沉降和基床表層沉降組成。由表可知，軌道結(jié)構(gòu)沉降主要發(fā)生在道床，與國外相關(guān)試驗研究結(jié)果相吻合。結(jié)合“2.1”節(jié)中軌道結(jié)構(gòu)受力特性，分析造成這種現(xiàn)象的原因為：基床表層所受荷載較道床已大幅衰減且分布均勻，使其產(chǎn)生較小沉降；基床表層結(jié)構(gòu)孔隙率相對較小，在列車荷載的長期反復作用下不易發(fā)生顆粒的錯動和重排，更容易達到沉降穩(wěn)定；底層道砟向基床表層的陷入也是道床下沉所占比例較大的主要原因，且隨著高鐵列車行車速度的提高、發(fā)車間隔的縮短，道床向基床表層的入侵量越來越大，使該因素對道床下沉的影響作用更為顯著。

表3 不同行車速度軌道結(jié)構(gòu)沉降數(shù)據(jù)

表4 不同列車軸重軌道結(jié)構(gòu)沉降數(shù)據(jù)

隨著行車速度的提高、列車軸重的增大，軌道結(jié)構(gòu)沉降中道床沉降所占比例也在逐漸增大，即行車速度、列車軸重的增大對道床沉降的影響更為顯著。其原因在于，荷載的增大會導致由道床下陷和由道砟錯動引起的道床沉降變大；同時，行車速度的提高會導致荷載頻率增大，使軌道結(jié)構(gòu)彈性變形恢復時間縮短，累積塑性變形量明顯增大。

此外，可根據(jù)軌道垂向變形情況對線路養(yǎng)護維修工作量進行估算：行車速度由250 km/h提升至400 km/h，軌枕彈性變形量增大39.4%，軌道沉降中道床沉降所占比例由83.9%提升至86.8%，道床沉降增大42.1%，基床表層沉降增大13.0%；列車軸重由14 t提升至20 t，軌枕彈性變形量增大44.0%，軌道沉降中道床沉降所占比例由83.9%提升至85.3%，道床沉降增大33.9%，基床表層沉降增大19.9%。兩種情況下養(yǎng)護維修工作量均會成倍增加，需在高速鐵路有砟軌道的設(shè)計、建設(shè)和運營過程中重點關(guān)注。

4 結(jié) 論

本文考慮散體道床與基床表層相互作用，采用離散元法建立軌枕-道床-基床表層耦合模型，計算結(jié)果與秦沈客運專線現(xiàn)場實測結(jié)果較一致，驗證了模型的可靠性。利用該模型從軌道結(jié)構(gòu)受力、彈性變形、塑性變形等角度研究路基基礎(chǔ)上散體道床垂向變形機理，并分析了列車運營條件對結(jié)構(gòu)變形特性的影響，主要結(jié)論如下：

(1) 軌道結(jié)構(gòu)內(nèi)部接觸力呈錐臺狀向下擴散，在錐臺形狀斷層處以斜面形式向結(jié)構(gòu)右下方(與行車方法相反)延伸。加載位置會影響最大接觸力存在位置，但不會影響接觸力傳遞路徑。

(2) 存在道砟與基床表層碎石相互入侵現(xiàn)象。軌道結(jié)構(gòu)沉降主要發(fā)生在道床，但基床表層沉降量也不可忽略，在本文模擬列車運營條件下，道床沉降占軌道整體沉降的比例為83.9%～86.8%。

(3) 軌道結(jié)構(gòu)垂向彈性和塑性變形量、道床沉降所占比例均隨著行車速度、列車軸重的增大非線性增加，且增加速率逐漸變大，應在高速鐵路有砟軌道的設(shè)計、建設(shè)和運營過程中重點關(guān)注。

[1] 佐藤吉彥,周柔琴.新軌道破壞理論[J].鐵道建筑, 1980(1):25-27.

SATO Yoshihiko, ZHOU Rouqin. The new failure theory of track structure[J].Railway Engineering, 1980(1):25-27.

[2] 平野雅之,周柔琴.采用新軌道破壞理論試算軌道下沉量[J].鐵道建筑,1981(6):31-36.

HIRANO Masayuki, ZHOU Rouqin. Calculation of track settlement based on the new failure theory[J]. Railway Engineering, 1981(6):31-36.

[3] 曾樹谷.鐵路散粒體道床[M].北京:中國鐵道出版社,1997.

[4] 胡仁偉,王紅,趙國堂.道碴動三軸試驗研究[J].中國鐵道科學,2001,22(2):101-106.

HU Renwei, WANG Hong, ZHAO Guotang. Study on dynamic triaxial test and accumulated deformation of ballast[J]. China Railway Science, 2001,22(2):101-106.

[5] 王紅.道碴彈性與累積變形的試驗研究[J].中國鐵道科學,2001,22(6):106-110.

WANG Hong. Test and study on ballast elasticity and accumulated deformation[J]. China Railway Science, 2001, 22(6):106-110.

[6] MOHAMMAD F, ESTORFFLI O V. Simulation of the dynamic behaviour of bedding-foundation soil in The time domain[J]. System Dynamics and Long-term Behaviour of Railway Vehicles, Track and Subgrade, 2003(6):357-374.

[7] FERREIRA P A, LPEZ-PITA A. Numerical modeling of high-speed train/track system to assess track vibrations and settlement prediction[J]. Journal of Transportation Engineering, 2013,139(3):330-337.

[8] 高建敏,翟婉明,徐涌.有碴軌道下沉變形參數(shù)影響分析[J].交通運輸工程學報,2007,7(4):15-20.

GAO Jianmin, ZHAI Wanming, XU Yong. Analysis of parameters’ effects on settlement of ballasted track[J]. Journal of Traffic and Transportation Engineering, 2007,7(4): 15-20.

[9] 徐鵬,蔡成標.列車-有砟軌道-路基空間耦合動力學模型[J].工程力學,2011,28(3):191-197.

XU Peng, CAI Chengbiao. Spatial dynamic model of train-ballast track-subgrade coupled system[J]. Engineering Mechanics, 2011,28(3):191-197.

[10] 徐旸,高亮,楊國濤，等.道砟顆粒的精細化建模及對道床力學性能影響[J].鐵道學報,2014,36(9):73-78.

XU Yang, GAO Liang, YANG Guotao, et al. Fine modelling of ballast grains and influence on mechanical properties of track bed[J]. Journal of the China Railway Society, 2014,36 (9):73-78.

[11] 付龍龍,宮全美,周順華，等.列車荷載作用下有砟軌道軌面沉降與路基不均勻沉降間的相關(guān)關(guān)系[J].振動與沖擊, 2013,32(14):23-28.

FU Longlong, GONG Quanmei, ZHOU Shunhua, et al. Relationships between ballasted track settlement and differential subgrade settlement under train load[J]. Journal of Vibration and Shock, 2013,32(14):23-28.

[12] 王平,萬復光.鐵路碎石道床彈性特性研究初探[J].鐵道學報,1997,19(4):108-114.

WANG Ping, WAN Fuguang. Research on the elasticity of ballast[J]. Journal of the China Railway Society, 1997,19(4): 108-114.

[13] 中華人民共和國鐵道部. 高速鐵路設(shè)計規(guī)范：TB 10621—2014[S].北京:中國鐵道出版社,2014.

[14] LIM W L. Mechanics of railway ballast behavior[D]. Nottingham: University of Nottingham.2004.

[15] LU M, MC-DOWELL G R. The importance of modelling ballast particle shape in the discrete element method[J]. Granular Matter,2007,9(1):69-80.

[16] AURSUDKIJ B, MC-DOWELL G R, COLLOP A C. Cyclic loading of railway ballast under triaxial conditions and in a railway test facility[J]．Granular Matter, 2009(11): 391-401.

[17] USACE (U.S. Army Corps of Engineers) Railroad design and construction[M]. Washington, DC: Track Design Manual, 2000.

[18] TEW J. A review of track design procedures: Sleepers and ballast[M]. Melbourne, Australia: Railways of Australia, 1991.

[19] SUN Q D, INDRARATNA B, NIMBALKAR S. Deformation and degradation mechanisms of railway ballast under high frequency cyclic loading[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016,142(1)：04015056.

[20] 董亮,趙成剛,蔡德鉤，等. 高速鐵路路基的動力響應分析方法[J].工程力學,2008,25(11):231-240.

DONG Liang, ZHAO Chenggang, CAI Degou, et al. Method for dynamic response of subgrade subjected to high-speed moving load[J]. Engineering Mechanics, 2008,25(11):231-240.