ccg注浆过程的数值模拟

CCG注浆过程的数值模拟研究

本文基于有限变形理论，采用点面接触模型模拟注浆过程中浆泡与土体之间的相互作用关系，通过拉格朗日方法引入这种相互作用关系，用约束变分原理建立适合模拟CCG注浆过程的控制方程。然后根据在试验中观察到现象建立的注浆过程模型编制的有限元程序可以动态模拟CCG注浆过程，为CCG注浆施工的优化设计奠定了基础。

【关键词】：CCG注浆有限变形点面接触模型数值模拟Abstract:Onthebasisofthefinitedeformationtheory，thetheoreticalstudyandnumericalsimulationaremadeoncompactiongrouthing.ByusingmodifiedLagrangianformula，thesimulationistakenonviscouscontactbetweenmortarbubbleandsoilgroundbypointansurfacecontactmodel.Acontrolformulaissetuptoadapttocompactiongroutingbyfiniteelementmethod，andthecompactiongroutingmodelisestablishedonthebasisoftestingobservation.Finally，thefiniteelementprogramisdevelopedusingFortranlanguageintendedtodynamicallysimulatetheinjectionprocess.

Keywords:compactiongrouting，finitedeformation，pointandsurfacecontactmodel，numericalanalysis.

1引言

在加固地基土层的各种注浆方法中，无论是劈裂注浆、渗透注浆还是喷射注浆，注入的浆液均是易于流动的，而CCG注浆注入的浆材是坍落度很小、甚至近似为零的不流动的浆材。这使得CCG注浆方法具有比其他注浆方式更好的可控制性，注浆施工时浆液可以准确地注入到土层中的预定位置，浆液注入到地下土层后不进入土体的孔隙，也不劈裂土层，在均匀土体中形成一个整体向周围扩大并挤密周围土体的浆泡，通过提升或下降注浆管，在注浆管的端部就可以产生一系列局部相“重叠”的浆泡，形成一柱状注浆体（图1），并挤密周围土层，达到加固地基的目的。

因为注浆浆液的特殊性，CCG注浆的加固机理不同于其他的注浆方法，加上诸多的注浆变量如注浆孔的间距、注浆体长度、注浆速率、注浆压力、注浆量等，CCG注浆施工设计涉及到一组复杂的参数。从世界范围来看虽然人们对这种注浆方法的机理进行理论研究并提出了各种不同的注浆设计方法，但理论研究水平仍然滞后于实际应用发展的需要，迄今为止，还没有一套成熟的注浆设计方法能够合理地考虑涉及到具体工程的不同的参数并对这种注浆方式进行优化设计。要解决这个问题，最关键的方面在于有效地动态模拟CCG注浆过程。采用何种方法描述CCG注浆过程中浆泡与土体之间的相互关系也是一个关键问题。K.Komiya、AdelM.EI-Kelesh和K.Soga等人都采用了穴孔扩张（Cavityexpansion）理论模拟浆泡与土体之间的接触，在浆泡与土体之间的分界面上通过人为地施加各向均匀的分布力模拟CCG注浆过程，但用有限元进行分析时,在每个荷载步中施加多大的力是不易确定的。另外，采用这个方法只能模拟单个浆泡的情况，不能模拟连续的注浆过程。对于单个浆泡的注入过程和连续的注浆过程，可以采用点面接触模型，在CCG注浆过程中浆泡与土体之间没有相对滑动，变形前后接触点的局部坐标值相同，它们之间的接触是粘滞接触。研究中将注浆泡与土体作为一个系统中的两个相互接触的物体，那么，注浆泡体和土体离散化后，注浆泡与土体之间的相互接触就可以模拟为接触面上土体的结点与注浆泡表面的任意点的接触。另外，应该注意到当介质在荷载作用下发生的位移（与自身的几何尺度相比）和应变（与1相比）很小时，建立物体和微元体的平衡条件时可以不考虑物体的位置和形状的变化，因此，分析中不必区分变形前和变形后的构形，相应的变形可以用无限小的线性应变进行度量。CCG注浆过程中靠近浆泡的土体不符合小变形和小应变假设，此时必须考虑变形对平衡的影响，即平衡条件应建立在变形后的构形上，所以，应采用有限变形理论分析CCG注浆过程中注浆体周围土体的变形问题。

本文基于有限变形理论，采用点面接触模型模拟注浆过程中浆泡与土体之间的相互作用关系，通过拉格朗日方法引入这种相互作用关系，用约束变分原理建立适合模拟CCG注浆过程的控制方程。然后编制了可以动态模拟CCG注浆过程的有限元程序，用该程序对鞍山路地铁车站基坑加固注浆进行了模拟计算和分析，并将模拟计算结果与现场实测值进行了对比。

2有限变形条件下的CCG注浆过程的控制方程

属

之

元

资

疼

绞赵3种

听己模崭拟牲C瓦C阳G扑注未浆钩过胁程脚的营有乓限啄元棋分由析

寿

隆

划

若

钳

弯异3涉.姥1鸦泽注滴浆赛模盐型恒的捎建截立牧与绝系栽统缎的露离只散桌化

扁

脏

女

鬼

出

验挑通每过震试僻验坏发冰现页，具注蒸入筐土桥体骄中怎的纠浆暑液秤不梅进剑入腥土惰体使孔哨隙德，庸在分均视质旗土卵体教中战形披成辜一拒个父球桥体寿，陡或准通业过驾均泻匀灶提慰升迷注阿浆数管衫可叉形龄成克均鸟匀奋的忘柱抱状船注友浆砌体律，晃注怀浆贼体冠与疼土曾体枕之像间栋有碧明达显扇的裕分最界魂面旨。

根据以上观测到的现象，可以建立图1和图2膊所键示残的谷注含浆扎模订型代，赞其歪中剧，境图饮1诱表翼示仪土纳体贤中萌已织形荷成陡的负一扬个泻长蔽度埋（甩高繁度功）乏为兵H泻，诚直贴径炮约匹为剧d垫的焦柱南桩攻注招浆励体派，魄图拳2际表扯示费柱率状阻注铜浆全体彩的宝形治成段过漏程犁，裤先狭形妨成阿一袭个巡独窝立颗的缘球约体站，月这悟个炒球晓体妨由摩小核到倦大到，脆其学半圆径启由佳R深1砌增吧大枕到板半廊径展R您2汇，抗再斗增施大良到脾半乔径亩R迎3袋；奥第属一司个乌球壮体衫生未成寇好同之会后宽，傍接刊着艳在阔第桌一阁个岩球刑体孙上捕方扎生钞成伏第唱二耽个飞球佳体偿，毫第殖二个个挖球蚁体伞的巩生潜成肚也然是季由染小群到隐大伍，蜂同筑样蛙其厉半痰径么由盈R撕1沃增假大阳到糠半潜径归R顶2工，确再鸭增涌大殿到拌半恢径狸R榆3脊，得但氧是叶第劲二河个酷球搭体之与药第辰一比个晋球锄体薪有的部疼分孕重说迭且，卖注锋浆强的挨实肯际百情散况故就剧是搂这秩样台。屡同舍理置，感可胁生导成以第湾三喝个恨、景第筛四涨个陵、夸第珠五为个取和球更漏多浮个胃球园体哀，哥形蠢成系不牧断呈增偶长妖的筐柱邻状新注眯浆第体使。支本啦文拔模忆拟欢注架浆索过垃程驼的恐计爹算皇程刃序旨就邀是校按键照凡这屋样犯的吼注腾浆印模饶型孕编毛制明的芒，老从男而杨模尿拟座出桐连芬续忍的近注释浆袜过不程卖。

提

蚂

令

裂

哈

冷占然终后端，哀将朴柱丈状尚注影浆针体亡和档土跟体零作桂为膀一瓜个摘整展体房或戴一灵个吹系泳统亲来刻考邪虑姻，茧但形对吓这议二匆者镰进嫌行寨离认散贡化任时群可第以倍分键开虹处敏理惠，护注铸浆恶泡晴的镜离舰散买见舅图昨3汤所也示拜，岂第绢一少个假浆器泡愉表凑面猛离刮散务成陶1血2盒个世线津性歇单阿元树，骨其坐他滥浆剧泡业也旨作言相降应凶类跟似割离残散他。获在宁注辜浆弟压戚力朝的王作申用刮下渔，值注柏浆层过膨程勤中抗浆天泡侵的仍变贤形捞相提对托于肠土帮体货来筒说清较户小形，界因斩此蛮，刻在涝实天际慧处仔理什中歌可涨以区将我其周近创似必地荡看螺成习刚哈体告，果也忙可良以选将村其芝当彩成质变率形副体绣，记这徒可匠以硬通素过寄设膛定跪变匙形状模摊量杆的料大以小干来听实耀现策。孕土贸体该的似离焰散租需达要饲仔网细值考源虑烦，刮应煤该吸适阿当丘进秃行疲土全体嗽的捞网俗格持划烦分叨，糊若抱网控格畜划骗分不得凭太深细准，疏一燃方阶面罩随握着瞎节钩点拍数牵和现单虎元帜数猾的怒增茎加眯，魂计冻算佳工骄作惯量忠会卷大岭量圈增胀加仔；锤另显一史方坝面干，毫若衔划睁分症的赏网礼格模太恢小哈，抽由幕于窜注么浆跌体具变宜面捕附幼近腐的时土佩体蚂会忠发费生倘大鞭变螺形黄，科因薯而凉有拉些乒单败元华形次状薄会滩发禽生棉畸疯变深，什致民使墙该王单昆元袭的业雅核可搏比孕矩沿阵使的购逆意阵粱不难存觉在笋，腿也奏就弟无阁法组进揭行禁计物算促了细。

才

隶

极

枣

透

断凭3低.芹2延铁C罪C场G葱注重浆右的用模叉拟麻计宵算

敌

器

拦

户

芒

而钥本宝文版的死程决序躺是莫按欢空图间转轴忙对盾称忽问冬题架编蚀制涝的碍，屑共习有振3搂5滴个黑子嫌程仍序英和归一浆个乖主插程肃序强。口程寒序馒考欺虑嗓了五几这何漠非梢线院性驼和路材工料书非击线盖性吨以厕及乐接界触颠边此界拨条妻件吗非卫线寿性呆的蚁有煌限酸元棒分击析剑程故序伪。社该绩程爆序箭可岁以鞭模崖拟立连弟续震的赚注听浆挽过勇程洞，没能除够钞计些算要出下注花浆洲过歇程穗中线任腰意泄时浮刻偿注尊浆算体纹周凉围诊土敬体狼的烦位德移势场晕，英还捐能掠计腹算迅出娱注刮浆陈体都周信围数土驶体搭的链塑逐性恨变颂形造范发围现。纺土驼体墓单疲元筋采挑用吧了箩四溉节楚点方线采性棵等芽参智单街元卡，追接晒触扁边巡界袋元丘为挨线览性剪单屠元触。块土笑的重本虚构轰模训型童是贞考饲虑骑M变o晴h投r捡-夫C宰o且u图l各o闻m副b牌屈迁服老准纽则研的刷弹毕塑逝性岗模闹型违。还这削里侵对蓝鞍姜山厌新惜村蝶地惜铁盆站座基杜坑舱的誓中巾间胳标蓝准请段尊的龄C樱C爆G岭注领浆映过马程温进句行巧模壶拟烦计车算翻，叉并啦将芹计纪算锡结招果碧与医在扬标迹准双段穿进脏行尘的吊单垒孔拜注做浆块试闲验京过严程撕中戚的跃实盖测摄结化果霜进由行坐比振较衣。

禁

趁

队

班

动

热傻根召据彻提衣供赛的书地绍质走勘睛察跨报田告屿，属各抢土扶层搬的逼物搭理敞力拉学查性走能元指侵标妖见柱表扫1邀。

计算深度取为21.0m，计算半径取为10.0m。边界条件为：底面和右侧面的径向位移和竖向位移为零；上表面自由；左侧面为对称轴。计算中泊松比取0.49，其中的压缩模量应换算为变形模量。有限元网格划分见图4，采用四节点线性等参单元，结点总数为903，单元总数为840。在车站基坑中间标准段进行的单孔注浆试验中，注浆从地表下16.5m开始，最后形成了一个半径为25～27cm、长度为3.0m的柱状注浆体，模拟计算时也是从地表下16.5m开始注浆，模拟的注浆体长度与半径基本上与试验中情况相同。部分计算结果如图5～图12所示。图5和图6表示注浆体周围土体在竖值方向的位移，其中，图5表示注浆过程中地表面的隆起值与注浆体长度之间的关系。图6表示注入第一个浆泡后不同深度的土层在竖直方向的位移值。图7和图8表示注浆过程中注浆体周围土体在水平方向上的位移。其中，图7表示距离注浆体中心线为2.5m处土体在水平方向的位移值与注浆体长度的变化（即注浆量和注浆深度的变化）之间的关系。图8表示注浆体长度达到3.0m时，距离注浆体中心线不同处的土体的水平方向位移值的变化。可以看到，注浆的开始阶段，位移的增量明显；随着注浆体长度的增加，水平方向的位移量增加不多，水平位移的峰值点随着注浆体长度的增长而逐渐上移。图9～图12中，将试验中的实测值与相应的计算值进行了比较。计算值的大小与实测值的大小接近，但曲线的形状不完全吻合。对于软粘土来说，因为侧斜管本身有一定的刚度，它不可能完全真实地反映出土体的实际变形，模拟的曲线形状比实测的变形曲线更能真实的反映实际变形情况。

4结语

（1）本文基于有限变形理论，采用点面接触模型模拟注浆过程中浆泡与土体之间的粘滞接触形式，根据在试验中观察到现象建立注浆过程的模型是恰当的。

（2）文中通过拉格朗日方法结合变分原理建立适合了模拟CCG注浆过程的控制方程，编制的有限元程序可以动态模拟CCG注浆过程，通过对鞍山路地铁车站基坑加固注浆过程进行了模拟计算和分析，模拟计算结果与现场试验实测值基本吻合。

（3）通过本文对CCG注浆过程的动态模拟方法，为进一步进行CCG注浆施工的优化设计奠定了基础。

参考文献

［1］AdelM.El-Kelesh,IsmailM.Basha,modelofcompactiongrouting,Journalofgeotechnicalandgeoenvironmentalengineering,2001,november,p955～964

［2］PantusoD.,BouzinovPA.,BatheKJ,afiniteelementprocedurefortheanalysisofthermo-mechanicalsolidsincontact,Computers&Structures,2000,Vol.75,No.6

［3］Carter,J.P.,Booker,J.R.&Yeung,S.K.,Cavityexpansionincohesivefrictionalsoils.Geotechnique36,No.3,1986,349-353

［4］Graf,EdwardD.CompactionGroutingTechniqueandObservations,JournaloftheSoilMechanicsandFound