列車循環(huán)荷載下水泥土復(fù)合地基沉降的顆粒離散元計(jì)算

關(guān)振長(zhǎng)，廖重輝，羅志彬，郭光照

(1.福州大學(xué)土木工程學(xué)院，福建福州350116;2.中國(guó)土木工程集團(tuán)福州勘察設(shè)計(jì)研究院有限公司，福建福州350001)

列車循環(huán)荷載下水泥土復(fù)合地基沉降的顆粒離散元計(jì)算

關(guān)振長(zhǎng)1，廖重輝1，羅志彬1，郭光照2

(1.福州大學(xué)土木工程學(xué)院，福建福州350116;2.中國(guó)土木工程集團(tuán)福州勘察設(shè)計(jì)研究院有限公司，福建福州350001)

依托福州江陰港鐵路支線水泥土攪拌樁軟基處理試驗(yàn)工程，探討利用顆粒離散元方法計(jì)算水泥土復(fù)合地基在列車循環(huán)荷載作用下的工后沉降問(wèn)題。首先根據(jù)原場(chǎng)地淤泥質(zhì)土的三軸壓縮試驗(yàn)和水泥土的單軸壓縮試驗(yàn)所得的應(yīng)力—應(yīng)變曲線，利用顆粒離散元建立三維數(shù)值模擬模型，并結(jié)合統(tǒng)計(jì)學(xué)習(xí)方法，對(duì)復(fù)合地基中的淤泥質(zhì)土和水泥土顆粒的細(xì)觀參數(shù)進(jìn)行反演;然后按江陰港水泥土攪拌樁復(fù)合地基的實(shí)際尺寸建立縮尺的離散元數(shù)值模型，進(jìn)行列車荷載下的沉降計(jì)算。結(jié)果表明:在500～1 000趟列車通過(guò)后水泥土攪拌樁復(fù)合地基的工后沉降趨于穩(wěn)定，沉降量約為20 mm。這一結(jié)果與按現(xiàn)有規(guī)范計(jì)算所得沉降量較為接近。

水泥土攪拌樁 列車循環(huán)荷載 工后沉降 顆粒離散元

水泥土攪拌樁是一種利用水泥作為固化劑，就地將原位土和水泥漿液強(qiáng)制攪拌，使原位土硬結(jié)成具有整體性、水穩(wěn)性和一定強(qiáng)度的水泥加固土的軟基處理方法。該法因成本低廉、形式靈活和對(duì)環(huán)境影響小等優(yōu)點(diǎn)，在軟基處理工程中得到廣泛應(yīng)用。國(guó)內(nèi)外許多學(xué)者通過(guò)室內(nèi)試驗(yàn)、現(xiàn)場(chǎng)測(cè)試及數(shù)值模擬等方法，對(duì)水泥土攪拌樁的力學(xué)特性、加固機(jī)理及設(shè)計(jì)方法進(jìn)行了大量的研究［1-3］。但受樁土相互作用的復(fù)雜性以及研究手段的制約，已有研究主要局限于考慮水泥土復(fù)合地基的承載力，而甚少探討它的工后沉降，尤其是未能考慮在列車循環(huán)荷載作用下水泥土攪拌樁復(fù)合地基的工后沉降問(wèn)題［4-6］。

對(duì)于復(fù)合地基的工后沉降分析，常面臨如何模擬樁土接觸這一難題。一般而言，對(duì)于剛性樁復(fù)合地基，在有限元等常規(guī)的數(shù)值方法中可采用接觸面單元加以解決，其樁土接觸面參數(shù)可參考室內(nèi)試驗(yàn)或現(xiàn)場(chǎng)測(cè)試結(jié)果來(lái)進(jìn)行選取［7-8］。但水泥土攪拌樁是原位攪拌固化成型的半剛性樁，受成樁過(guò)程中強(qiáng)力攪拌擾動(dòng)的影響，攪拌樁與原位土之間的接觸面常常并不明顯，若簡(jiǎn)單套用剛性樁的樁土接觸面參數(shù)，將給沉降計(jì)算帶來(lái)誤差甚至錯(cuò)誤。反之，若不考慮樁土接觸問(wèn)題，即假設(shè)水泥土攪拌樁和原位軟土始終共同變形，又與實(shí)際明顯不符。另外，若要深入探討水泥土攪拌樁復(fù)合地基在列車循環(huán)荷載作用下的工后沉降機(jī)理，還須設(shè)定水泥土與原位軟土的動(dòng)力本構(gòu)及相應(yīng)的動(dòng)參數(shù)，而這方面可資利用的研究成果也甚為少見(jiàn)。因此，常規(guī)的數(shù)值方法尚存諸多不足和局限，有必要嘗試新的思路和計(jì)算方法。

自上世紀(jì)七十年代Cundall提出適用于巖土力學(xué)的顆粒離散元方法以來(lái)，在學(xué)者們的推動(dòng)下，該方法已成功應(yīng)用于許多室內(nèi)巖土試驗(yàn)的數(shù)值模擬與機(jī)理研究，并得以在巖土和地下工程領(lǐng)域廣泛應(yīng)用。與采用有限元等常規(guī)方法計(jì)算水泥土復(fù)合地基沉降的不同之處在于，在顆粒離散元方法中可設(shè)定水泥土顆粒、原位軟土顆粒的接觸本構(gòu)及相應(yīng)的接觸參數(shù)，通過(guò)細(xì)觀層面上不同顆粒之間的接觸力與位移變化，來(lái)實(shí)現(xiàn)對(duì)水泥土復(fù)合地基宏觀力學(xué)性能的模擬。此法可巧妙繞開(kāi)半剛性樁復(fù)合地基中樁土接觸模擬和動(dòng)力參數(shù)選取的難題，同時(shí)在計(jì)算中采用顯式時(shí)間差分動(dòng)態(tài)求解技術(shù)，很適用于列車循環(huán)荷載之類的動(dòng)力問(wèn)題研究［9-12］。

福州江陰港鐵路支線為Ⅱ級(jí)貨運(yùn)鐵路，其線路穿越地區(qū)上覆有16～20 m厚的濱海相欠固結(jié)淤泥質(zhì)軟土層。線路下設(shè)計(jì)采用樁長(zhǎng)12 m的水泥土攪拌樁進(jìn)行地基加固，樁徑0.5 m，梅花形布置，水泥摻量12%。本文以福州江陰港鐵路支線水泥土攪拌樁軟基處理工程為背景，嘗試?yán)没赑FC3D數(shù)值平臺(tái)的顆粒離散元方法對(duì)水泥土攪拌樁復(fù)合地基在列車循環(huán)荷載作用下的工后沉降問(wèn)題進(jìn)行深入探討。

1 復(fù)合地基細(xì)觀參數(shù)的反演

1.1 淤泥質(zhì)土與水泥土的基本物性

在攪拌樁加固施工前，對(duì)場(chǎng)地淤泥質(zhì)土進(jìn)行了6個(gè)孔位的原位鉆探取樣，進(jìn)行了常規(guī)土工試驗(yàn)和三軸壓縮試驗(yàn)，試樣直徑D×高H=3.9 cm×8 cm，得到淤泥質(zhì)土的基本物理性能指標(biāo)如表1所示。

表1 加固處理前淤泥質(zhì)土的物理性能指標(biāo)

在水泥土攪拌樁成樁90 d后，對(duì)其中18根水泥土攪拌樁進(jìn)行全長(zhǎng)抽芯取樣，進(jìn)行了常規(guī)土工試驗(yàn)和單軸壓縮試驗(yàn)，試樣尺寸D×H=5 cm×10 cm，得到水泥土的基本物理性能指標(biāo)如表2所示。

表2 加固處理后水泥土的物理性能指標(biāo)

1.2 顆粒接觸本構(gòu)及參數(shù)

在顆粒離散元數(shù)值分析中，顆粒的接觸本構(gòu)及其參數(shù)的選擇至為關(guān)鍵。對(duì)淤泥質(zhì)軟土和水泥土這兩種具有黏聚力的材料選用粘結(jié)型顆粒接觸本構(gòu)是較為合適的。該本構(gòu)模型假設(shè):離散顆粒間的接觸僅發(fā)生在接觸點(diǎn)位置;顆粒間通過(guò)法向彈簧傳遞法向力，通過(guò)剪切彈簧傳遞剪切力;接觸力的大小由顆粒的法向剛度kn、切向剛度ks、法向抗拉強(qiáng)度nb、切向剪切強(qiáng)度sb和顆粒間的摩擦系數(shù)μ這5個(gè)參數(shù)共同決定。當(dāng)任意兩個(gè)顆粒A與B相互接觸時(shí)，總法向剛度Kn與總切向剛度Ks可由下式得到

在計(jì)算中，離散顆粒之間的法向接觸力Fn與切向接觸力Fs可根據(jù)上一計(jì)算步求得的顆粒間相對(duì)法向位移un與切向位移us由式(3)和式(4)求得。需要注意的是，若計(jì)算所得法向接觸力Fn超過(guò)了閾值nb，可認(rèn)為是法向黏結(jié)破壞，F(xiàn)n被重置為0;若切向接觸力Fs超過(guò)了閾值sb+μFn，則認(rèn)為是切向黏結(jié)破壞，F(xiàn)s被重置為μFn。

本文將根據(jù)室內(nèi)試驗(yàn)所得應(yīng)力和應(yīng)變數(shù)據(jù)，通過(guò)顆粒離散元數(shù)值模擬與統(tǒng)計(jì)學(xué)習(xí)相結(jié)合的方法，分別反演得到淤泥質(zhì)土顆粒和水泥土顆粒的kn，ks，nb，sb和μ這5個(gè)細(xì)觀參數(shù)。

1.3 淤泥質(zhì)土三軸壓縮過(guò)程的數(shù)值模擬

本文以PFC3D作為數(shù)值分析平臺(tái)，對(duì)淤泥質(zhì)土UU三軸壓縮試驗(yàn)進(jìn)行全過(guò)程動(dòng)態(tài)仿真模擬，具體思路和步驟如下:

1)根據(jù)淤泥質(zhì)土三軸壓縮試樣的實(shí)際尺寸(即D ×H=3.9 cm×8.0 cm)，構(gòu)建一個(gè)圓柱形的剛性模型箱，如圖1所示。在箱內(nèi)隨機(jī)填充不同尺寸的圓形顆粒，其直徑在0.8～1.2 mm之間均勻分布，以模擬淤泥質(zhì)土試樣的初始狀態(tài)。根據(jù)常規(guī)土工試驗(yàn)測(cè)得的淤泥質(zhì)土試樣孔隙率，計(jì)算填充所需的顆粒數(shù)目，以避免填充過(guò)密或過(guò)疏。

圖1 淤泥質(zhì)土試樣UU三軸壓縮過(guò)程的數(shù)值模擬示意

2)設(shè)定淤泥質(zhì)土顆粒之間的一組接觸本構(gòu)參數(shù)，即顆粒的kn，ks，nb，sb和μ 5個(gè)參數(shù)，在自重條件下計(jì)算至初始平衡。

3)將模型箱周圍的圓柱形墻體設(shè)置為按應(yīng)力伺服控制，即始終保持一定的圍壓。在每一計(jì)算步前，先求得各顆粒與模型箱墻體之間的平均接觸應(yīng)力，再控制本計(jì)算步內(nèi)圓柱形墻體的外擴(kuò)或內(nèi)縮速度，以保證本計(jì)算步后顆粒與模型箱墻體平均接觸應(yīng)力即圍壓保持不變。

4)將墻體模型箱的上下底面墻體設(shè)置為按定常速度進(jìn)行控制，即每一計(jì)算步壓縮變形量相等，總計(jì)算步按試驗(yàn)過(guò)程中試樣的總壓縮量控制;同時(shí)記錄模擬過(guò)程中土樣中部的平均軸向應(yīng)力，繪制主應(yīng)力差與軸向應(yīng)變曲線。此處平均軸向應(yīng)力是指特定體積內(nèi)所有顆粒間接觸力之和在軸向上的集度。

1.4 基于統(tǒng)計(jì)學(xué)習(xí)方法的參數(shù)反演

本文采用數(shù)值模擬與基于神經(jīng)網(wǎng)絡(luò)的統(tǒng)計(jì)學(xué)習(xí)相結(jié)合的方法，對(duì)土顆粒的5個(gè)細(xì)觀參數(shù)進(jìn)行反演，即要找到一組最佳參數(shù)組合，使得由PFC3D數(shù)值模擬得到的軸向應(yīng)力差與軸向應(yīng)變曲線在最大程度上與試驗(yàn)結(jié)果相吻合［13-14］。參數(shù)反演的基本過(guò)程為:

1)利用誤差反向傳播神經(jīng)網(wǎng)絡(luò)構(gòu)建一個(gè)統(tǒng)計(jì)學(xué)習(xí)機(jī)。學(xué)習(xí)機(jī)為一個(gè)5×5×5的神經(jīng)網(wǎng)絡(luò)，即輸入層對(duì)應(yīng)五維的輸入向量(含5個(gè)神經(jīng)元)，輸出層和中間單隱藏層也設(shè)定為5個(gè)神經(jīng)元;

2)根據(jù)計(jì)算經(jīng)驗(yàn)分別確定淤泥質(zhì)土5個(gè)細(xì)觀參數(shù)的取值范圍，并在各參數(shù)的上下限范圍內(nèi)隨機(jī)取值生成一組參數(shù)組合，即輸入向量;

3)將該輸入向量代入PFC3D數(shù)值模型進(jìn)行計(jì)算，記錄軸向應(yīng)變?yōu)?%，3%，5%，7%，9%時(shí)的主應(yīng)力差，即輸出向量;

4)把這樣一組輸入輸出向量對(duì)作為一對(duì)學(xué)習(xí)樣本，重復(fù)以上步驟共生成若干對(duì)學(xué)習(xí)樣本;

5)將以上學(xué)習(xí)樣本輸入到統(tǒng)計(jì)學(xué)習(xí)機(jī)中，通過(guò)不斷訓(xùn)練和調(diào)整自身權(quán)重，最終在統(tǒng)計(jì)意義上用學(xué)習(xí)機(jī)的映射過(guò)程取代PFC3D的數(shù)值計(jì)算;

6)利用訓(xùn)練完畢的學(xué)習(xí)機(jī)，結(jié)合遺傳算法，不斷優(yōu)化輸入向量，最終尋得一個(gè)最佳輸入向量，使得其對(duì)應(yīng)的輸出向量盡可能接近試驗(yàn)結(jié)果;

7)將尋得的最佳輸入向量再次代入PFC3D數(shù)值模型中，驗(yàn)證其對(duì)應(yīng)的輸出向量是否接近試驗(yàn)結(jié)果;

8)若相差較大則舍棄之，并按以上方法重新尋找一個(gè)最佳輸入向量，直到滿足精度要求為止。

1.5 淤泥質(zhì)土細(xì)觀參數(shù)的反演

利用上述方法，以深度5.5 m處的淤泥質(zhì)土試樣為例，在圍壓100 kPa下對(duì)試樣的細(xì)觀參數(shù)進(jìn)行反演，最終得到一組最佳細(xì)觀參數(shù)為:ks=49.1 MPa，kn/ks= 2.01，nb=4.9 kN，sb/nb=9.8，μ=0.05。將其代入PFC3D數(shù)值模擬中進(jìn)行驗(yàn)證，并繪制三軸壓縮試驗(yàn)和數(shù)值模擬所得的應(yīng)力應(yīng)變曲線，如圖2所示。由計(jì)算結(jié)果可知，兩條曲線在軸向應(yīng)變?yōu)?%，3%，5%，7%，9%處的平均誤差僅為3.4%，滿足精度要求。

同樣，對(duì)深度16 m處的淤泥質(zhì)土試樣，在300 kPa圍壓下利用數(shù)值模擬與統(tǒng)計(jì)學(xué)習(xí)相結(jié)合的方法，反演得到深部淤泥質(zhì)土顆粒的一組最佳細(xì)觀參數(shù)為:ks= 48.0 MPa，kn/ks=2.05，nb=5.1 kN，sb/nb=9.4，μ= 0.06。將其代入PFC3D數(shù)值模型中進(jìn)行驗(yàn)證，繪制的應(yīng)力應(yīng)變曲線如圖3所示。與試驗(yàn)曲線相比，數(shù)值模擬平均誤差僅為2.9%，滿足精度要求。

圖2 淤泥質(zhì)土試樣UU三軸壓縮過(guò)程的應(yīng)力—應(yīng)變曲線(100 kPa)

圖3 淤泥質(zhì)土試樣UU三軸壓縮過(guò)程的應(yīng)力—應(yīng)變曲線(300 kPa)

由反演結(jié)果可知，不同深度處淤泥質(zhì)土顆粒的細(xì)觀參數(shù)比較接近，這與原場(chǎng)地內(nèi)淤泥質(zhì)土均勻性較好的實(shí)際情況是一致的。

1.6 水泥土細(xì)觀參數(shù)的反演

對(duì)于成樁90 d后的水泥土，也可進(jìn)行類似的細(xì)觀參數(shù)反演。以深度為2.0 m和8.0 m處的芯樣為例，反演得到淺部水泥土顆粒的一組最佳細(xì)觀參數(shù)為:ks=69.4 MPa，kn/ks=1.45，nb=13.6 kN，sb/nb=3.7，μ=0.75;深部水泥土顆粒的一組最佳細(xì)觀參數(shù)為: ks=53.0 MPa，kn/ks=1.58，nb=8.1 kN，sb/nb=4.2，μ=0.63。

將此2組細(xì)觀參數(shù)分別代入PFC3D數(shù)值模型進(jìn)行驗(yàn)證，得到的應(yīng)力應(yīng)變曲線分別繪制于圖4和圖5中。與試驗(yàn)曲線相比，數(shù)值模擬平均誤差分別為4.8%和5.3%，滿足精度要求。

圖4 淺部水泥土芯樣單軸壓縮過(guò)程的應(yīng)力—應(yīng)變曲線

圖5 深部水泥土芯樣單軸壓縮過(guò)程的應(yīng)力—應(yīng)變曲線

由反演結(jié)果可知，不同深度處水泥土顆粒的細(xì)觀參數(shù)差別較大，這與淺部軟土水泥攪拌樁加固質(zhì)量較好，而深部加固質(zhì)量較差的實(shí)際情況相吻合。

2 列車循環(huán)荷載作用下的工后沉降

2.1 水泥土攪拌樁復(fù)合地基的離散元數(shù)值模型

在上文對(duì)淤泥質(zhì)土和水泥土的細(xì)觀參數(shù)進(jìn)行反演時(shí)，顆粒粒徑設(shè)定為0.8～1.2 mm之間均勻分布。若按這一粒徑來(lái)模擬實(shí)際尺寸的水泥土攪拌樁復(fù)合地基(樁徑0.5 m，樁長(zhǎng)12 m，樁間距1.2 m)，對(duì)其進(jìn)行沉降分析，則需要上億個(gè)離散元顆粒，這大大超出了微型工作站的運(yùn)算能力。因此本文根據(jù)相似理論進(jìn)行縮尺的數(shù)值模擬［15］。由量綱分析法獲得的各物理量的相似比如表3所示。根據(jù)前人研究結(jié)果［7-10］，當(dāng)模型最小幾何尺寸與顆粒最大幾何尺寸的比值≥10時(shí)，可忽略顆粒粒徑對(duì)數(shù)值模擬結(jié)果造成的尺寸效應(yīng)，因此縮尺數(shù)值模擬中n的取值為100。即水泥土攪拌樁復(fù)合地基的縮尺數(shù)值模型，其內(nèi)部水泥土攪拌樁的直徑為5 mm，深度為120 mm;其周邊淤泥質(zhì)土的直徑為12 mm，深度為180 mm。淤泥質(zhì)土和水泥土顆粒的5個(gè)細(xì)觀參數(shù)按上節(jié)選取，但重力加速度根據(jù)相似關(guān)系取為980 m/s2。

表3 縮尺數(shù)值模擬的相似比

2.2 列車循環(huán)荷載的施加

為便于研究，本文根據(jù)文獻(xiàn)［16］將列車豎向荷載表示為

式中:P0為列車恒載，根據(jù)《鐵路路基設(shè)計(jì)規(guī)范》(TB 10001—2005)，江陰港Ⅱ級(jí)貨運(yùn)鐵路支線P0取為60.1 kPa;f為荷載頻率，根據(jù)文獻(xiàn)［16］，當(dāng)列車速度低于250 km/h時(shí)，其振動(dòng)頻率主要集中在4 Hz左右。根據(jù)相似關(guān)系，數(shù)值模擬中所采用的荷載頻率為400 Hz。

由此作用于復(fù)合地基頂面的荷載

圖6 水泥土攪拌樁復(fù)合地基工后沉降的經(jīng)時(shí)變化

2.3 工后沉降計(jì)算

按500 000個(gè)計(jì)算步進(jìn)行模擬，每個(gè)計(jì)算步代表的模擬時(shí)間為0.000 21 s(此為程序根據(jù)顆粒數(shù)量、應(yīng)力水平等因素自行設(shè)定的計(jì)算參數(shù))，則相當(dāng)于進(jìn)行了約105 s的列車循環(huán)加載的動(dòng)態(tài)模擬。由此所得水泥土攪拌樁復(fù)合地基的工后沉降與模擬時(shí)間的關(guān)系曲線如圖6所示。由該圖可以看出，經(jīng)歷模擬時(shí)間100 s后(按相似關(guān)系即實(shí)際時(shí)間為10 000 s后)，水泥土攪拌樁復(fù)合地基的工后沉降趨于穩(wěn)定，其中攪拌樁樁頂沉降為0.202 mm，按相似比其實(shí)際沉降為20.2 mm;攪拌樁中部6 m深度處沉降為0.107 mm，實(shí)際沉降為10.7 mm，攪拌樁底部12 m深度處的沉降僅為0.01 mm，實(shí)際沉降僅為1 mm。

需要說(shuō)明的是，實(shí)際時(shí)間是指列車荷載實(shí)際作用在攪拌樁復(fù)合地基上的時(shí)間，假定一趟列車通過(guò)某處復(fù)合地基的所需時(shí)間為10～20 s，那么實(shí)際時(shí)間10 000 s則相當(dāng)于有500～1 000趟列車通過(guò)。

2.4 應(yīng)力面積法的沉降計(jì)算結(jié)果

為便于比較，根據(jù)《鐵路路基設(shè)計(jì)規(guī)范》(TB 10001—2005)和《建筑地基基礎(chǔ)設(shè)計(jì)規(guī)范》(GB 50007—2011)，計(jì)算了場(chǎng)地水泥土攪拌樁復(fù)合地基在列車荷載下的沉降。將列車視為寬度3.5 m，集度60.1 kPa的條形荷載，則水泥土攪拌樁復(fù)合地基的壓縮模量Esp

式中:m為置換率，根據(jù)設(shè)計(jì)取為0.136;Es為淤泥質(zhì)土的壓縮模量，根據(jù)室內(nèi)土工試驗(yàn)結(jié)果取為3.4 MPa; Ep為水泥土的壓縮模量，根據(jù)芯樣的單軸壓縮試驗(yàn)，其平均值取為50 MPa。

經(jīng)計(jì)算求得復(fù)合地基壓縮模量Esp為9.69 MPa。最后采用應(yīng)力面積法，得到福州江陰港鐵路支線水泥土攪拌樁復(fù)合地基在列車荷載作用下的總沉降量為14.1 mm，與上節(jié)的計(jì)算結(jié)果較為接近，這表明采用顆粒離散元方法進(jìn)行沉降計(jì)算是可行的。

3 結(jié)論

針對(duì)有限元等常規(guī)數(shù)值方法在計(jì)算水泥土復(fù)合地基沉降時(shí)遇到的樁土接觸參數(shù)和動(dòng)力本構(gòu)參數(shù)選取等諸多不易解決的問(wèn)題，本文提出了利用顆粒離散元方法研究水泥土攪拌樁復(fù)合地基在列車循環(huán)荷載作用下的工后沉降計(jì)算問(wèn)題的整體思路和方法:

1)首先根據(jù)室內(nèi)地基土的試驗(yàn)數(shù)據(jù)，對(duì)淤泥質(zhì)土顆粒和水泥土顆粒的細(xì)觀參數(shù)進(jìn)行反演計(jì)算;

2)再按相似原理，對(duì)水泥土攪拌樁復(fù)合地基建立縮尺的數(shù)值模型;

3)最后將反演所得的細(xì)觀參數(shù)應(yīng)用于縮尺模型的數(shù)值模擬，求得水泥土攪拌樁復(fù)合地基工后沉降的經(jīng)時(shí)變化。

針對(duì)福州江陰港Ⅱ級(jí)貨運(yùn)鐵路支線軟基處理的實(shí)際，通過(guò)離散元數(shù)值模擬得到在500～1 000趟列車通過(guò)后，場(chǎng)地水泥土攪拌樁復(fù)合地基的工后沉降約20 mm，并趨于穩(wěn)定。這一結(jié)果與采用傳統(tǒng)的應(yīng)力面積法計(jì)算所得的工后沉降值較為接近，證明了采用顆粒離散元方法研究水泥土攪拌樁復(fù)合地基在列車循環(huán)荷載作用下的工后沉降問(wèn)題的可行性。

［1］HORPIBULSUK S，CHINKULKIJNIWAT A，CHOLPHATSORN A，et al.Consolidation Behavior of Soil-cement Column Improved ground［J］.Computers and Geotechniques，2012，43(1):37-50.

［2］張偉麗，蔡健，林奕禧，等.水泥土攪拌樁復(fù)合地基荷載傳遞機(jī)理的試驗(yàn)研究［J］.土木工程學(xué)報(bào)，2010，43(6):116-121.

［3］向瑋，劉松玉，經(jīng)緋，等.深長(zhǎng)變徑攪拌樁荷載傳遞規(guī)律的試驗(yàn)研究［J］.巖土力學(xué)，2010，31(9):2765-2771.

［4］趙春風(fēng)，李永剛，錢濤.水泥土攪拌樁復(fù)合地基樁土應(yīng)力比的解析算法［J］.中南大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012，43 (6):2390-2395.

［5］MAEDA Y，ICHIKAWA A，TSURUKUBO S，et al.Study on Bearing Capacity Characteristics and Practical Application of Composite Ground Pile Foundation［J］.Japanese Geotechnical Journal，2008，3(1):37-54.

［6］安關(guān)峰，張洪彬，劉添俊.旋噴群樁復(fù)合地基承載特性的數(shù)值分析［J］.巖土力學(xué)，2012，33(3):906-912.

［7］張洪彬，安關(guān)峰，劉添俊，等.接觸面剛度對(duì)雙排樁支護(hù)結(jié)構(gòu)的影響分析［J］.地下空間與工程學(xué)報(bào)，2012，8(2): 301-306.

［8］VOOTTIPRUEX P，SUKSAWAT T，BERGADO D，et al. Numerical Simulations and Parametric Study of SDCM and DCM Piles under Full Scale Axial and Lateral Loads［J］. Computers and Geotechnics，2011，38(3):318-329.

［9］羅勇.土工問(wèn)題的顆粒流數(shù)值模擬及應(yīng)用研究［D］.杭州:浙江大學(xué)，2007.

［10］周健，陳小亮，王冠英，等.開(kāi)口管樁沉樁過(guò)程試驗(yàn)研究與顆粒流模擬［J］.同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012，40 (2):173-178.

［11］ELMEKATIA，SHAMYU.APracticalCo-simulation Approach for Multiscale Analysis of Geotechnical Systems［J］.Computers and Geotechnics，2010，37(4):494-503.

［12］Itasca Consulting Group Inc.PFC3D User Manual Version 3.0［R］.Minnesota:Itasca Consulting Group Inc.，2009.

［13］PICHLER B，LACKNER R，MANG H.Back Analysis of Model Parameters in Geotechnical Engineering by Means of Soft Computing［J］.International Journal of Numerical Methods in Engineering，2003，57(11):1943-1978.

［14］GUAN Z，JIANG Y，TANABASHI Y.Rheological Parameter Estimation for the Prediction of Long-term Deformations in Conventional Tunneling［J］.Tunnelling and Underground Space Technology，2009，24(3):250-259.

［15］杜延齡，韓連兵.土工離心模型試驗(yàn)技術(shù)［M］.北京:中國(guó)水利水電出版社，2010.

［16］王常晶.列車移動(dòng)荷載作用下地基的動(dòng)應(yīng)力及飽和軟粘土特性研究［D］.杭州:浙江大學(xué)，2006.

Particle discrete element calculation for subsidence of cement mixing pile composite foundation under train cyclic load

GUAN Zhenchang1，LIAO Chonghui1，LUO Zhibin1，GUO Guangzhao2
(1.College of Civil Engineering，F(xiàn)uzhou University，F(xiàn)uzhou Fujian 350116，China;2.China Civil Engineering Construction Corporation，F(xiàn)uzhou Survey and Design Institute Co.，Ltd.，F(xiàn)uzhou Fujian 350001，China)

Based on the soft ground treatment with cement mixing piles project in Jiangyin railway sub-line，the particle discrete element method(DEM)is used to explore the post-construction subsidence of cement mixing pile composite foundation with train cyclic loads.According to the stress-strain curve from triaxial compression and uniaxial compression lab tests，the micro-parameters of silty and cement soil particles are back-analyzed through an integration methodology，which consists of statistical analysis method and 3D numerical simulation model based on particle discrete element.T hen built the reduced scale DEM model according to the real scale of cement mixing pile composite foundation in Jiangyin Port to calculate the subsidence with train load.From the result，the postconstruction subsidence of cement-mixing pile composite foundation tends to stable around 20 mm after 500～1 000 times of train cyclic load，and this subsidence value is close to the subsidence calculated by existing standard.

Cement mixing pile;T rain cyclic load;Post-construction subsidence;Particle discrete element method

TU472.6

A

10.3969/j.issn.1003-1995.2015.12.28

1003-1995(2015)12-0104-06

(責(zé)任審編周彥彥)

2015-06-21;

2015-08-11

教育部留學(xué)回國(guó)人員科研啟動(dòng)基金(LXKQ0904);國(guó)家自然科學(xué)基金(51008082)

關(guān)振長(zhǎng)(1980—)，男，副教授，博士。