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

1、ccg注浆过程的数值模拟研【提要】：本文基于有限变形理论，采用点面接触模 型模拟注浆过程中浆泡与土体之间的相互作用关系，通过拉 格朗日方法引入这种相互作用关系，用约束变分原理建立 适合模拟ccg注浆过程的控制方程。然后根据在试验中观察 到现象建立的注浆过程模型编制的有限元程序可以动态模拟 ccg注浆过程，为ccg注浆施工的优化设计奠定了基础。【关 键词】:ccg注浆有限变形点面接触模型数值模拟abstract: on the basis of the finite def orma t ion theory, the theoretical study and numerical simula

2、tion are made on compaction grouthing. by using modified lagrangian forrnula, the simulation is taken on viscous contact between mortar bubble and soil ground by point an surface contact modela control formula is set up to adapt to compaction grouting by fini te element method, and the compaction gr

3、outing model is established on the basis of testing observationfinally， the finite element program is developed using fortran language intended to dynamically s imulate the injection process .keywords: compaction grouting, finite deformation, point and surface contact model, numerical analysis.1 弓|

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

5、的注浆变量如注浆孔的间距、注浆 体长度、注浆速率.注浆压力、注浆量等，ccg注浆施工设计 涉及到一组复杂的参数。从世界范围来看虽然人们对这种注 浆方法的机理进行理论研究并提出了各种不同的注浆设计 方法，但理论研究水平仍然滞后于实际应用发展的需要， 迄今为止，还没有一套成熟的注浆设计方法能够合理地考虑 涉及到具体工程的不同的参数并对这种注浆方式进行优化设 计。要解决这个问题，最关键的方面在于有效地动态模拟 ccg注浆过程。采用何种方法描述ccg注浆过程中浆泡与土体 之间的相互关系也是一个关键问题。k. komiya、adel m. ei-kelesh和ksoga等人都采用了穴孔扩张(cavity

6、 expansion )理论模拟浆泡与土体之间的接触，在浆泡与土体之间的分界面上通过人为地施加各向均匀的分布力模拟ccg 注浆过程，但用有限元进行 分析时，在每个荷载步中施加多大的力是不易确定的。另外，采用这个方法只能模拟单个浆 泡的情况，不能模拟连续的注浆过程。对于单个浆泡的注入 浆过程中浆泡与土体之间没有相对滑动，变形前后接触点的 局部坐标值相同，它们之间的接触是粘滞接触。研究中将注 浆泡与土体作为一个系统中的两个相互接触的物体，那么， 注浆泡体和土体离散化后，注浆泡与土体之间的相互接触就过程和连续的注浆过程，可以采用点面接触模型，在ccg注可以模拟为接触面上土体的结点与注浆泡表面的任意点

7、的接 触。另外，应该注意到当介质在荷载作用下发生的位移（与 自身的几何尺度相比）和应变（与1相比）很小时，建立物 体和微元体的平衡条件时可以不考虑物体的位置和形状的变 化，因此，分析中不必区分变形前和变形后的构形，相应的 变形可以用无限小的线性应变进行度量。ccg注浆过程中靠近 浆泡的土体不符合小变形和小应变假设，此时必须考虑变形 对平衡的影响，即平衡条件应建立在变形后的构形上，所 以，应采用有限变形理论分析ccg注浆过程中注浆体周围土 体的变形问题。本文基于有限变形理论，采用点面接触模型模拟注浆过 程中浆泡与土体之间的相互作用关系，通过拉格朗日方法引 入这种相互作用关系，用约束变分原理建立适

8、合模拟ccg注 浆过程的控制方程。然后编制了可以动态模拟ccg注浆过程 的有限元程序，用该程序对鞍山路地铁车站基坑加固注浆进 行了模拟计算和分析，并将模拟计算结果与现场实测值进 行了对比。2有限变形条件下的ccg注浆过程的控制方程3模拟ccg注浆过程的有限元分析3.1注浆模型的建立与系统的离散化通过试验发现，注入土体中的浆液不进入土体孔隙，在 均质土体中形成一个球体，或通过均匀提升注浆管可形成均 匀的柱状注浆体，注浆体与土体之间有明显的分界面。根据 以上观测到的现象，可以建立图1和图2所示的注浆模型， 其中，图1表示土体中已形成的一个长度（高度）为h，直径 约为d的柱桩注浆体，图2表示柱状注浆

9、体的形成过程，先 形成一个独立的球体，这个球体由小到大，其半径由r1增大到半径r2,再增大到半径r3;第一个球体生成好之后，接着 在第一个球体上方生成第二个球体，第二个球体的生成也是 由小到大，同样其半径由r1增大到半径r2，再增大到半径 r3，但是第二个球体与第一个球体有部分重迭，注浆的实际 情况就是这样。同理，可生成第三个、第四个、第五个和更 多个球体，形成不断增长的柱状注浆体。本文模拟注浆过程 的计算程序就是按照这样的注浆模型编制的，从而模拟出连 续的注浆过程。然后，将柱状注浆体和土体作为一个整体或一个系统来 考虑，但对这二者进行离散化时可以分开处理，注浆泡的离 散见图3所示，第一个浆泡

10、表面离散成12个线性单元，其他 浆泡也作相应类似离散。在注浆压力的作用下，注浆过程中 浆泡的变形相对于土体来说较小，因此，在实际处理中可以 将其近似地看成刚体，也可以将其当成变形体，这可以通过 设定变形模量的大小来实现。土体的离散需要仔细考虑，应 该适当进行土体的网格划分，若网格划分得太细，一方面随 着节点数和单元数的增加，计算工作量会大量增加；另一方 面，若划分的网格太小，由于注浆体变面附近的土体会发生 大变形，因而有些单元形状会发生畸变，致使该单元的雅可 比矩阵的逆阵不存在，也就无法进行计算了。3.2 ccg注浆的模拟计算本文的程序是按空间轴对称问题编制的，共有35个子 程序和一个主程序。

11、程序考虑了几何非线性和材料非线性以 及接触边界条件非线性的有限元分析程序。该程序可以模拟连续的注浆过程，能够计算出注浆过程中任意时刻注浆体周围土体的位移场，还能计算出注浆体周土体的塑性变形范围。土体单元采用了四节点线性等参单元，接触边界元为线 性单元。土的本构模型是考虑mohr-coulomb屈服准则的弹塑 性模型。这里对鞍山新村地铁站基坑的中间标准段的ccg注 浆过程进行模拟计算，并将计算结果与在标准段进行的单孔 注浆试验过程中的实测结果进行比较。根据提供的地质勘察报告，各土层的物理力学性能指标 见表1。计算深度取为21. om,计算半径取为10. 0mo边界条件 为：底面和右侧面的径向位移

12、和竖向位移为零；上表面自由； 左侧面为对称轴。计算中泊松比取0.49,其中的压缩模量应 换算为变形模量。有限元网格划分见图4,采用四节点线性等 参单元，结点总数为903，单元总数为840。在车站基坑中间 标准段进行的单孔注浆试验中，注浆从地表下16.5m开始， 最后形成了一个半径为2527cm、长度为3. om的柱状注浆 体，模拟计算时也是从地表下16.5m 始注浆，模拟的注浆 体长度与半径基本上与试验中情况相同。部分计算结果如图 5图12所示。图5和图6表示注浆体周围土体在竖值方向的位移，其中，图5表示注浆过程中地表面的隆起值与注浆 体长度之间的关系。图6表示注入第一个浆泡后不同深度的 土层

13、在竖直方向的位移值。图7和图8表示注浆过程中注浆 体周围土体在水平方向上的位移。其中，图7表示距离注浆 体中心线为2. 5m处土体在水平方向的位移值与注浆体长度的 变化（即注浆量和注浆深度的变化）之间的关系。图8表示 注浆体长度达到3. om时，距离注浆体中心线不同处的土体的 水平方向位移值的变化。可以看到，注浆的开始阶段，位移 的增量明显；随着注浆体长度的增加，水平方向的位移量增 加不多，水平位移的峰值点随着注浆体长度的增长而逐渐上 移。图9图12中，将试验中的实测值与相应的计算值进行 了比较。计算值的大小与实测值的大小接近，但曲线的形状不完全吻合。对于软粘土来说，因为侧斜管本身有一定的刚

14、度，它不可能完全真实地反映出土体的实际变形，模拟的曲 线形状比实测的变形曲线更能真实的反映实际变形情况。4结语(1 )本文基于有限变形 理论，采用点面接触模型模拟注浆过程中浆泡与土体之间的粘滞接触形式，根据在试验 中观察到现象建立注浆过程的模型是恰当的。(2 )文中通过拉格朗日方法 结合变分原理建立适合 了模拟ccg注浆过程的控制方程，编制的有限元程序可以动 态模拟ccg注浆过程，通过对鞍山路地铁车站基坑加固注浆 过程进行了模拟计算和分析，模拟计算结果与现场试验实 测值基本吻合。(3)通过本文对ccg注浆过程的动态模拟方法，为进 一步进行ccg注浆施工的优化设计奠定了基础。参考文献1 adel

15、 m. el-kelesh, ismai 1 m. basha, model of compaction grout ing, journal of geotechnical and geoenvironmental engineering, 2001, november, p955 9642 pantuso dbouzinov pa. , bathe kj, a finite element procedure for the analysis of thermo-mechanica 1 solids in contact , computersstructures, 2000, vol.7

16、5, no. 63 carter, j. p. , booker, j. r. yeung, s. k.,cavity expansion in cohesive f r ict ional soiis geotechnique 36, no. 3, 1986, 349-3534 graf, edward d. compaction grouting technique and observations, journal of the soil mechanics and foundations division, asce, 1969, pll5111585 ning hu, a solution method for dynamic contact problems, computersstructures, vol. 63, no. 6, 1997, 1053-10636 warner, james, compaction grouting mechanism 一 what do we know, grouting proceedings of the geo-institute, asce, geotechnical special publication no. 66, 199