数值模拟在粉喷桩复合路基中运用.docVIP

数值模拟在粉喷桩复合路基中运用 　　摘要：本文以在建的大庆至广州高速公路湖北省麻城至浠水段软土路基处理监测与试验研究为背景，选取粉喷桩复合路基试验段，在现场一年多的变形观测的基础上，借助有限差分数值分析方法对粉喷桩复合路基进行了沉降计算和模拟，模拟实际填土加载过程，将模拟计算得到的结果和现场实测结果进行对比分析。结果表明，经过粉喷桩处理后的路基在沉降等方面显著降低，模拟结果和现场试验值比较吻合，可将此方法应用于粉喷桩复合地基沉降计算中。同时对整个复合地基的工作性状进行了综合研究，对路基变形的相关影响因素进行了探讨。 　　关键词：粉喷桩，复合路基，数值模拟，沉降计算 　　 　　1.模型的建立 　　本次模型选取在建的大庆至广州高速公路湖北省麻城至浠水段第四合同段断面K40+500断面处，根据设计资料，填土高度为6m，断面路堤坡度为1：1.5，填土完毕后的路面宽度为26m，路堤两坡脚之间间距为44m，选取5倍路面宽为路基横断面总长，即130m。由于路堤具有几何对称性，可以采用一半模型来进行计算分析。路基模型根据勘察静力触探资料，自上而下分为四层，第一层为亚粘土，厚4m；第二层为淤泥质土，厚2.5m；第三层为亚砂土，厚1.5m；第四层为粗砂，厚2m。粉喷桩复合路基模型以原状地基为基础，在路堤范围内采用粉喷桩加固，等边三角形布置，桩间距1.5m，加固宽度为路堤坡脚外一排桩，桩长7m，穿透软土层，进入粗砂层。 　　1.1 路基模型分析步骤 　　（1）生成原状路基几何模型，对路基各层土进行赋值； 　　（2）初始应力计算，将路堤网格设置为空单元，地基单元赋值为Mohr-Coulomb模型，清除初始应力后填土时采用“group”命令进行激活填土。 　　（3）每步按0.6m分层，由“step”命令来控制填土加载时间。 　　（4）路基竖直沉降、水平位移状态分析。 　　12 模型及结构单元的选择 　　1.2.1 模型的选择 　　（1）本构模型 　　原状路基采用摩尔-库仑（Mohr-Coulomb）为本构模型，路堤填土采用弹性模型进行加载。当路堤未进行填筑时采用空单元模型（null），当进行填土加载时使用“group”命令即可进行路堤的加载填??。 　　（2）计算模型 　　路堤填土施工计算采用分阶段的弹塑性求解方法。在填土加载前，将路堤部分设置成空模型（null），将原状地基赋值为Mohr-Coulomb模型，先将Mohr-Coulomb模型的粘结力c值和抗拉强度赋值为无穷大进行求解，保证在重力作用下单元不至于发生屈服，再将Mohr-Coulomb模型参数赋值为真实值，进行求解计算。在填土加载模拟施工时要讲初始应力计算过程中产生的节点位移和速度进行清零处理，保证与实测值的可比性。 　　1.2.2 结构单元的选择 　　原状路基土和路堤填土都采用六面实体单元（brick）进行模拟。 　　粉喷桩采用桩单元（pile）模拟，pile单元属于二维单元，可在网格单元间传递剪切力、法向力及弯矩，能够承受点荷载和分布荷载，因此用pile单元模拟与地基块体既有切向力及法向力的相互作用的结构物。 　　1.3 模型参数 　　原状路基土采用Mohr-Coulomb模型，所需要进行计算的参数有体积模量K、剪切模量G、粘结力c、摩擦角以及密度p；路堤填土采用弹性模型，所需要进行计算的参数有体积模量K、剪切模量G及密度P。其中粘结力、摩擦角及密度都可以根据实验数据得到，体积模量和剪切模量根据下式得到： 　　体积模量：（2-10） 　　 　　剪切模量：（2-11） 　　 　　式中：E为弹性模量，在此处使用压缩模量进行计算比较，为泊松比。 　　所有模型参数如下表1-1所示： 　　表1-1 土层计算参数 　　 　　复合路基中粉喷桩采用结构桩单元pile，计算所需要的参数有：桩横截面积、截面惯性矩、外周长、桩密度、弹性模量、泊松比、切向耦合弹簧刚度、切向耦合弹簧刚度的粘聚力、切向耦合弹簧的摩擦角、法向耦合弹簧刚度、法向耦合弹簧的粘聚力、法向弹簧的摩擦角。 　　2.粉喷桩复合路基成果计算分析 　　2.1 竖向沉降变形分析 　　 　　图1-1现场实测沉降值与数值模拟值对比图 　　从图1-1看到，数值模拟软件有效地模拟了填土加载过程，在开始和结尾沉降值基本一致，中间过程值稍许偏大，与模拟计算时步和计算沉降结果取值有关，数值计算一直处于差分的平衡计算过程中。总体来说，计算沉降值与实测现场实测值较为吻合，说明将此数值计算方法应用于复合路基沉降计算中，可达到良好的效果，与实际情况吻合。 　　2.2 水平位移变形分析 　　 从图1-2可以看出： 　　（1）随着填土加载的进行，路堤在发生竖直沉降的同时，由于上层加载填土对下部软卧