CCG注漿過(guò)程的數(shù)值模擬研究

【提　要】：本文基于有限變形理論，采用點(diǎn)面接觸模型模擬注漿過(guò)程中漿泡與土體之間的相互作用關(guān)系，通過(guò)拉格朗日方法引入這種相互作用關(guān)系，用約束變分原理建立適合模擬CCG注漿過(guò)程的控制方程。然后根據(jù)在試驗(yàn)中觀察到現(xiàn)象建立的注漿過(guò)程模型編制的有限元程序可以動(dòng)態(tài)模擬CCG注漿過(guò)程，為CCG注漿施工的優(yōu)化設(shè)計(jì)奠定了基礎(chǔ)。
【關(guān)鍵詞】：CCG注漿有限變形點(diǎn)面接觸模型數(shù)值模擬

Abstract: On the basis of the finite deformation theory，the theoretical study and numerical simulation are made on compaction grouthing.By using modified Lagrangian formula，the simulation is taken on viscous contact between mortar bubble and soil ground by point an surface contact model.A control formula is set up to adapt to compaction grouting by finite element method，and the compaction grouting model is established on the basis of testing observation.Finally，the finite element program is developed using Fortran language intended to dynamically simulate the injection process.
Keywords: compaction grouting， finite deformation， point and surface contact model， numerical analysis.

1　引言
在加固地基土層的各種注漿方法中，無(wú)論是劈裂注漿、滲透注漿還是噴射注漿，注入的漿液均是易于流動(dòng)的，而CCG注漿注入的漿材是坍落度很小、甚至近似為零的不流動(dòng)的漿材。這使得CCG注漿方法具有比其他注漿方式更好的可控制性，注漿施工時(shí)漿液可以準(zhǔn)確地注入到土層中的預(yù)定位置，漿液注入到地下土層后不進(jìn)入土體的孔隙，也不劈裂土層，在均勻土體中形成一個(gè)整體向周?chē)鷶U(kuò)大并擠密周?chē)馏w的漿泡，通過(guò)提升或下降注漿管，在注漿管的端部就可以產(chǎn)生一系列局部相“重疊”的漿泡，形成一柱狀注漿體（圖1），并擠密周?chē)翆樱_(dá)到加固地基的目的。

因?yàn)樽{漿液的特殊性，CCG注漿的加固機(jī)理不同于其他的注漿方法，加上諸多的注漿變量如注漿孔的間距、注漿體長(zhǎng)度、注漿速率、注漿壓力、注漿量等，CCG注漿施工設(shè)計(jì)涉及到一組復(fù)雜的參數(shù)。從世界范圍來(lái)看雖然人們對(duì)這種注漿方法的機(jī)理進(jìn)行理論研究并提出了各種不同的注漿設(shè)計(jì)方法，但理論研究水平仍然滯后于實(shí)際應(yīng)用發(fā)展的需要，迄今為止，還沒(méi)有一套成熟的注漿設(shè)計(jì)方法能夠合理地考慮涉及到具體工程的不同的參數(shù)并對(duì)這種注漿方式進(jìn)行優(yōu)化設(shè)計(jì)。要解決這個(gè)問(wèn)題，最關(guān)鍵的方面在于有效地動(dòng)態(tài)模擬CCG注漿過(guò)程。采用何種方法描述CCG注漿過(guò)程中漿泡與土體之間的相互關(guān)系也是一個(gè)關(guān)鍵問(wèn)題。K. Komiya、Adel M. EI-Kelesh和K. Soga等人都采用了穴孔擴(kuò)張（Cavity expansion）理論模擬漿泡與土體之間的接觸，在漿泡與土體之間的分界面上通過(guò)人為地施加各向均勻的分布力模擬CCG注漿過(guò)程，但用有限元進(jìn)行分析時(shí),在每個(gè)荷載步中施加多大的力是不易確定的。另外，采用這個(gè)方法只能模擬單個(gè)漿泡的情況，不能模擬連續(xù)的注漿過(guò)程。對(duì)于單個(gè)漿泡的注入過(guò)程和連續(xù)的注漿過(guò)程，可以采用點(diǎn)面接觸模型，在CCG注漿過(guò)程中漿泡與土體之間沒(méi)有相對(duì)滑動(dòng)，變形前后接觸點(diǎn)的局部坐標(biāo)值相同，它們之間的接觸是粘滯接觸。研究中將注漿泡與土體作為一個(gè)系統(tǒng)中的兩個(gè)相互接觸的物體，那么，注漿泡體和土體離散化后，注漿泡與土體之間的相互接觸就可以模擬為接觸面上土體的結(jié)點(diǎn)與注漿泡表面的任意點(diǎn)的接觸。另外，應(yīng)該注意到當(dāng)介質(zhì)在荷載作用下發(fā)生的位移（與自身的幾何尺度相比）和應(yīng)變（與1相比）很小時(shí)，建立物體和微元體的平衡條件時(shí)可以不考慮物體的位置和形狀的變化，因此，分析中不必區(qū)分變形前和變形后的構(gòu)形，相應(yīng)的變形可以用無(wú)限小的線性應(yīng)變進(jìn)行度量。CCG注漿過(guò)程中靠近漿泡的土體不符合小變形和小應(yīng)變假設(shè)，此時(shí)必須考慮變形對(duì)平衡的影響，即平衡條件應(yīng)建立在變形后的構(gòu)形上，所以，應(yīng)采用有限變形理論分析CCG注漿過(guò)程中注漿體周?chē)馏w的變形問(wèn)題。
本文基于有限變形理論，采用點(diǎn)面接觸模型模擬注漿過(guò)程中漿泡與土體之間的相互作用關(guān)系，通過(guò)拉格朗日方法引入這種相互作用關(guān)系，用約束變分原理建立適合模擬CCG注漿過(guò)程的控制方程。然后編制了可以動(dòng)態(tài)模擬CCG注漿過(guò)程的有限元程序，用該程序?qū)Π吧铰返罔F車(chē)站基坑加固注漿進(jìn)行了模擬計(jì)算和分析，并將模擬計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)值進(jìn)行了對(duì)比。

2　有限變形條件下的CCG注漿過(guò)程的控制方程

3 模擬CCG注漿過(guò)程的有限元分析
3.1　注漿模型的建立與系統(tǒng)的離散化
通過(guò)試驗(yàn)發(fā)現(xiàn)，注入土體中的漿液不進(jìn)入土體孔隙，在均質(zhì)土體中形成一個(gè)球體，或通過(guò)均勻提升注漿管可形成均勻的柱狀注漿體，注漿體與土體之間有明顯的分界面。

根據(jù)以上觀測(cè)到的現(xiàn)象，可以建立圖1和圖2所示的注漿模型，其中，圖1表示土體中已形成的一個(gè)長(zhǎng)度（高度）為H，直徑約為d的柱樁注漿體，圖2表示柱狀注漿體的形成過(guò)程，先形成一個(gè)獨(dú)立的球體，這個(gè)球體由小到大，其半徑由R1增大到半徑R2，再增大到半徑R3；第一個(gè)球體生成好之后，接著在第一個(gè)球體上方生成第二個(gè)球體，第二個(gè)球體的生成也是由小到大，同樣其半徑由R1增大到半徑R2，再增大到半徑R3，但是第二個(gè)球體與第一個(gè)球體有部分重迭，注漿的實(shí)際情況就是這樣。同理，可生成第三個(gè)、第四個(gè)、第五個(gè)和更多個(gè)球體，形成不斷增長(zhǎng)的柱狀注漿體。本文模擬注漿過(guò)程的計(jì)算程序就是按照這樣的注漿模型編制的，從而模擬出連續(xù)的注漿過(guò)程。
然后，將柱狀注漿體和土體作為一個(gè)整體或一個(gè)系統(tǒng)來(lái)考慮，但對(duì)這二者進(jìn)行離散化時(shí)可以分開(kāi)處理，注漿泡的離散見(jiàn)圖3所示，第一個(gè)漿泡表面離散成12個(gè)線性單元，其他漿泡也作相應(yīng)類(lèi)似離散。在注漿壓力的作用下，注漿過(guò)程中漿泡的變形相對(duì)于土體來(lái)說(shuō)較小，因此，在實(shí)際處理中可以將其近似地看成剛體，也可以將其當(dāng)成變形體，這可以通過(guò)設(shè)定變形模量的大小來(lái)實(shí)現(xiàn)。土體的離散需要仔細(xì)考慮，應(yīng)該適當(dāng)進(jìn)行土體的網(wǎng)格劃分，若網(wǎng)格劃分得太細(xì)，一方面隨著節(jié)點(diǎn)數(shù)和單元數(shù)的增加，計(jì)算工作量會(huì)大量增加；另一方面，若劃分的網(wǎng)格太小，由于注漿體變面附近的土體會(huì)發(fā)生大變形，因而有些單元形狀會(huì)發(fā)生畸變，致使該單元的雅可比矩陣的逆陣不存在，也就無(wú)法進(jìn)行計(jì)算了。
3.2　CCG注漿的模擬計(jì)算
本文的程序是按空間軸對(duì)稱(chēng)問(wèn)題編制的，共有35個(gè)子程序和一個(gè)主程序。程序考慮了幾何非線性和材料非線性以及接觸邊界條件非線性的有限元分析程序。該程序可以模擬連續(xù)的注漿過(guò)程，能夠計(jì)算出注漿過(guò)程中任意時(shí)刻注漿體周?chē)馏w的位移場(chǎng)，還能計(jì)算出注漿體周?chē)馏w的塑性變形范圍。土體單元采用了四節(jié)點(diǎn)線性等參單元，接觸邊界元為線性單元。土的本構(gòu)模型是考慮Mohr-Coulomb屈服準(zhǔn)則的彈塑性模型。這里對(duì)鞍山新村地鐵站基坑的中間標(biāo)準(zhǔn)段的CCG注漿過(guò)程進(jìn)行模擬計(jì)算，并將計(jì)算結(jié)果與在標(biāo)準(zhǔn)段進(jìn)行的單孔注漿試驗(yàn)過(guò)程中的實(shí)測(cè)結(jié)果進(jìn)行比較。
根據(jù)提供的地質(zhì)勘察報(bào)告，各土層的物理力學(xué)性能指標(biāo)見(jiàn)表1。

計(jì)算深度取為21.0m，計(jì)算半徑取為10.0m。邊界條件為：底面和右側(cè)面的徑向位移和豎向位移為零；上表面自由；左側(cè)面為對(duì)稱(chēng)軸。計(jì)算中泊松比取0.49，其中的壓縮模量應(yīng)換算為變形模量。有限元網(wǎng)格劃分見(jiàn)圖4，采用四節(jié)點(diǎn)線性等參單元，結(jié)點(diǎn)總數(shù)為903，單元總數(shù)為840。在車(chē)站基坑中間標(biāo)準(zhǔn)段進(jìn)行的單孔注漿試驗(yàn)中，注漿從地表下16.5m開(kāi)始，最后形成了一個(gè)半徑為25～27cm、長(zhǎng)度為3.0m的柱狀注漿體，模擬計(jì)算時(shí)也是從地表下16.5m開(kāi)始注漿，模擬的注漿體長(zhǎng)度與半徑基本上與試驗(yàn)中情況相同。部分計(jì)算結(jié)果如圖5～圖12所示。圖5和圖6表示注漿體周?chē)馏w在豎值方向的位移，其中，圖5表示注漿過(guò)程中地表面的隆起值與注漿體長(zhǎng)度之間的關(guān)系。圖6表示注入第一個(gè)漿泡后不同深度的土層在豎直方向的位移值。圖7和圖8表示注漿過(guò)程中注漿體周?chē)馏w在水平方向上的位移。其中，圖7表示距離注漿體中心線為2.5m處土體在水平方向的位移值與注漿體長(zhǎng)度的變化（即注漿量和注漿深度的變化）之間的關(guān)系。圖8表示注漿體長(zhǎng)度達(dá)到3.0m時(shí)，距離注漿體中心線不同處的土體的水平方向位移值的變化。可以看到，注漿的開(kāi)始階段，位移的增量明顯；隨著注漿體長(zhǎng)度的增加，水平方向的位移量增加不多，水平位移的峰值點(diǎn)隨著注漿體長(zhǎng)度的增長(zhǎng)而逐漸上移。圖9～圖12中，將試驗(yàn)中的實(shí)測(cè)值與相應(yīng)的計(jì)算值進(jìn)行了比較。計(jì)算值的大小與實(shí)測(cè)值的大小接近，但曲線的形狀不完全吻合。對(duì)于軟粘土來(lái)說(shuō)，因?yàn)閭?cè)斜管本身有一定的剛度，它不可能完全真實(shí)地反映出土體的實(shí)際變形，模擬的曲線形狀比實(shí)測(cè)的變形曲線更能真實(shí)的反映實(shí)際變形情況。

4　結(jié)語(yǔ)
（1） 本文基于有限變形理論，采用點(diǎn)面接觸模型模擬注漿過(guò)程中漿泡與土體之間的粘滯接觸形式，根據(jù)在試驗(yàn)中觀察到現(xiàn)象建立注漿過(guò)程的模型是恰當(dāng)?shù)摹?BR> （2） 文中通過(guò)拉格朗日方法結(jié)合變分原理建立適合了模擬CCG注漿過(guò)程的控制方程，編制的有限元程序可以動(dòng)態(tài)模擬CCG注漿過(guò)程，通過(guò)對(duì)鞍山路地鐵車(chē)站基坑加固注漿過(guò)程進(jìn)行了模擬計(jì)算和分析，模擬計(jì)算結(jié)果與現(xiàn)場(chǎng)試驗(yàn)實(shí)測(cè)值基本吻合。
（3） 通過(guò)本文對(duì)CCG注漿過(guò)程的動(dòng)態(tài)模擬方法，為進(jìn)一步進(jìn)行CCG注漿施工的優(yōu)化設(shè)計(jì)奠定了基礎(chǔ)。
參考文獻(xiàn)
［1］ Adel M. El-Kelesh, Ismail M. Basha, model of compaction grouting, Journal of geotechnical and geoenvironmental engineering, 2001, november, p955～ 964
［2］ Pantuso D., Bouzinov PA., Bathe KJ, a finite element procedure for the analysis of thermo-mechanical solids in contact , Computers & Structures, 2000, Vol.75, No.6
［3］ Carter, J. P., Booker, J. R. & Yeung, S. K. , Cavity expansion in cohesive frictional soils. Geotechnique 36, No. 3, 1986,349-353
［4］ Graf, Edward D. Compaction Grouting Technique and Observations, Journal of the Soil Mechanics and Foundations Division, ASCE, 1969, p1151～1158
［5］ Ning Hu, A solution method for dynamic contact problems, Computers & Structures, Vol.63, No.6, 1997,1053-1063
［6］ Warner, James, Compaction Grouting Mechanism - What Do We Know, Grouting Proceedings of the Geo-Institute, ASCE, Geotechnical Special Publication No. 66, 1997, p1
［7］ 程曉，張鳳祥. 土建注漿施工與效果檢測(cè). 上海：同濟(jì)大學(xué)出版社，1998
［8］ 李向紅，傅德明，劉建航. CCG注漿及其工程應(yīng)用. 工業(yè)建筑，2003，12
［9］ 李向紅，傅德明. CCG注漿試驗(yàn)研究. 工業(yè)建筑，2003，11
［10］ 蘇世豐，廖洪鈞. 擠壓灌漿后粘土之正規(guī)化不排水剪力強(qiáng)度變化. 第八屆大地工程學(xué)術(shù)研究討論會(huì)論文集，1999，p665～662
［11］ 王勖成，邵敏. 有限元基本原理與數(shù)值方法. 北京：清華大學(xué)出版社，1988
［12］ 殷有泉. 固體力學(xué)非線性有限元引論. 北京：北京大學(xué)出版社、清華大學(xué)出版社，1987

新聞來(lái)源：《城市交通隧道工程最新技術(shù)》