施青公路典型边坡治理措施分析评价.doc

6 施青公路典型边坡治理措施分析评价 6.1 初始地应力的模拟 在ANSYS中，有两种方法可用来模拟初始地应力，一种是只考虑岩体的自重应力，忽略其构造应力，在分析的第一步，首先计算岩体的自重应力场。这种方法计算简单方便，只需要给出岩体的各项参数即可计算。不足指出在与计算出的应力场与实际应力场有偏差，而且岩体在自重作用下还产生了初始位移，在继续分析后续施工工序时，得到的位移结果是累加了初始位移的结果，而现实中初始位移早就结束，对边坡的开挖不产生影响，因此在以后的每个施工阶段分析位移场时，都需要减去初始位移场（LCDEF,1,1 LCOPER,SUB,1）。 在进行结构分析时，ANSYS中可以使用输入文件来把初应力指定为一种荷载，因此当具有实测的初始地应力资料时，可将初始地应力写成初应力荷载文件，然后读入ANSYS作为荷载条件（ISFILE,READ,文件名,ist），就可直接进行第一步的开挖计算。 6.2 开挖与支护过程的模拟 在ANSYS中，可以用杀死和激活单元来模拟材料的消去和添加。利用ANSYS这种单元的生死功能，可以简单有效的模拟边坡的开挖和支护过程。 杀死单元时，程序将通过用一个非常小的数乘以单元的刚度，并从总的质量矩阵消去单元的质量来实现“杀死”单元。边坡开挖时，可直接选择将被开挖掉的单元，然后ekill,all，即可实现开挖的模拟；施加支护时，可首先将相应支护部分在开挖时被杀死的单元激活（EALIVE），然后改变其材料性质（MPCHG）。当单元被激活时，他们具有零应变状态，且若大变形选项打开的话（NLGEOM,ON），他们的几何性质被修改来与他们现在的偏移位置相适应。 在一些情况下，单元的生死状态可以根据ANSYS的计算结果来决定，如应力、应变等。利用这一功能，在模拟过程中根据计算结果，可以将超过许用应力（线弹性分析时）或许用应变（弹塑性分析时）的单元杀死，以此来模拟围岩或结构的破坏。 6.3 连续施工的模拟 ANSYS程序中的荷载步（LOAD STEP）功能可以实现不同工况间的连续计算，用来模拟边坡的连续施工过程方便有效。开始建立整体有限元模型，包括要杀死或者激活的部分，模拟的过程中不需要重新划分网格。在一个荷载步计算结束后，可直接进行下道工序的施工：杀死或激活单元，改变虚拟支撑力等，然后求解计算。如此继续一直到施工结束。 值得注意的是，整个的连续计算过程应在求解器中完成。若想在每一荷载步求解完成后进入通用后处理器，以便及时查看每一工序的计算结果，应在开始计算之前设置重启动选项（RESTART），且生死单元也恢复到初始状态，需重新设置单元生死。 6.4 有限元模型 在模拟施青公路路基边坡时，本次的分析范围以地面线为上边界范围，左边界距地表54米，右边界距地表16米，模型共划分1053个单元，930个节点，其中锚杆135个单元，梁单元50个。围岩采用四节点二维实体单元（plane42平面应变）模拟，锚杆采用一维杆单元（link1）模拟；喷射混凝土采用二维梁（beam3）单元，边坡围岩构材料按均质弹塑性考虑，采用Drucker-Prager屈服准则，参数可由地质勘测资料得出； 为保证计算的准确性，模型尺寸为：边坡取m，边坡取m，竖直向上取至地表，。Tab.3.1 Surrounding rock tunnel structural materials basic physical and mechanical parameters 图3.3 有限元分析模型 Fig.3.3 Finite element analysis model 表3.1 边坡围岩及结构材料基本物理力学参数 Tab.3.1 Surrounding rock tunnel structural materials basic physical and mechanical parameters E /MPa /KN/m3 C /MPa /° Q4al+pl 18.27 160 0.30 14.22° 41.0 Q4al+pl 20.29 220 0.40 17.0° 47.0 ∈1p 26.85 400 0.40 24.0° 47.0 初期支护 28000 23 2 60 0.3 二次衬砌 31000 25 3 60 0.3 锚杆 200000 78 / / 0.3 6.5 施工过程模拟 计算过程为了尽可能真实地模拟边坡实际开挖支护作业，本次模拟顺序如下： 加自重进行初始应力计算； 边坡开挖数值模拟； 边坡支护数值模拟。 6.6计算结果分析 6.6.1初始状态分析 计算模型图 第一主应力S1分布图 第一主应力S3分布图 围岩Y方向位移UY 围岩初始主应力数值计算过程中,外载荷考