城市道路轻型汽车行驶工况构建

张宏; 姚延钢; 杨晓勤

doi:10.3969/j.issn.0258-2724.20180224

城市道路轻型汽车行驶工况构建

doi: 10.3969/j.issn.0258-2724.20180224

基金项目: 内蒙古自治区自然科学基金资助项目（2019MS07021）；内蒙古自治区交通运输厅科技项目（NJ-2017-8）；中央高校基本科研业务费资助资金项目（310822171134）；中国国家留学基金资助项目（201806815002）

详细信息

作者简介:
张宏（1981—），女，副教授，博士，研究方向为城市交通数据科学及应用，E-mail：zhang2001hong@126.com

中图分类号: U461
计量
- 文章访问数: 681
- HTML全文浏览量: 339
- PDF下载量: 34
- 被引次数: 0
出版历程
- 收稿日期: 2018-03-27
- 修回日期: 2018-09-07
- 网络出版日期: 2018-10-24
- 刊出日期: 2019-12-01

Light-Duty Vehicles Driving Cycle Construction Based on Urban Roads

摘要

摘要: 为了制定标准车辆能耗规范和标定车辆排放，优化汽车使用性能，以呼和浩特市区道路上轻型汽车为研究对象，对轻型汽车行驶工况进行了分析. 首先通过中国汽车检测工况研究和开发（China automotive test cycle，CATC）专用数据采集设备，收集了74台车辆的行驶工况样本数据，数据采集覆盖所有时段类型、道路类型、轻型车类型和驾驶员类型；其次通过加权二次构建控制不同车辆类型的比例，采用主成分分析和聚类分析预处理数据，制定短行程规则；最后将运动学片段进行裁剪和特征值分类，构建了城市道路轻型汽车行驶工况. 研究结果表明：CATC平均速度为25.87 km/h，运行平均速度为33.92 km/h、匀速比例为20.59%，怠速比例为23.72%，加速比例为28.56%，减速比例为27.13%. 与欧盟提出的轻型车循环测试工况相对比，平均速度、运行平均速度、匀速比例低于欧洲工况，加速比例、减速比例、怠速比例高于欧洲工况.
- 轻型汽车 /
- 行驶工况 /
- 短片段 /
- 城市交通
Abstract: In order to regulate vehicle energy consumption, calibrate vehicle emission and optimize vehicle performance, the driving cycles of light-duty vehicles on the urban roads of Hohhot is analyzed. Firstly, the 74 sample data of driving cycles were collected through special data acquisition equipments of CATC (China automotive test cycle). These samples cover tall types of roads, light-duty vehicles and drivers. Secondly, weights were assigned to different vehicles to set vehicle type proportion. The sample data are processed by principal component analysis and clustering analysis and the short segment rules are formulated. Finally the kinematic segments are clipped and their characteristics values are classified. Thus, the operating conditions of light-duty vehicles on urban roads are constructed. The results show that the average speed of CATC is 25.87 km/h, average running speed is 33.92 km/h, uniform speed ratio is 20.59%, idling ratio is 23.72%, acceleration ratio is 28.56% and deceleration ratio is 27.13%. Compared to the European Union (EU) cycle test condition, the average speed, average running speed and uniform speed ratio are lower than those of the EU standard, while the acceleration ratio, deceleration ratio and idling ratio are higher.
- light-duty vehicles /
- driving cycle /
- short segment /
- urban transportation

HTML全文

随着我国基础设施建设的持续推进，多座隧道在施工过程中遇到了高地应力软弱围岩大变形^[1-3]，针对这一复杂工程难题，众多学者提出了包括改进型钢拱架、高强预应力锚杆（索）、让压支护、多层支护等多种变形控制方式^[4-7]，其中，超前应力释放（超前导洞）是高地应力软岩大变形控制常用的措施，相关工程案例证明，超前应力释放可使围岩变形量下降30%～40%^[8-10].

周宗青等^[11-13]采用数值模拟对超前导洞的围岩变形控制效果进行了验证，并对导洞断面形状、布置方式以及其他影响因素进行了分析；崔光耀等^[14]在兰渝铁路某隧道通过现场试验对比了双台阶法、超前钻孔、超前导洞3种开挖方式对围岩的变形控制；王道远等^[15]在新莲隧道大变形段采用“超前导洞 + 扩挖”的方式对大变形进行控制后发现，围岩变形仍然有近30 cm；尤显明等^[16-17]以木寨岭隧道为依托验证了超前导洞与多层支护结构联合作用下对软岩变形的控制效果，但该方案围岩变形仍然较大，且施工工序复杂. 对文献分析后可知，尽管已有不少工程案例证明超前应力释放可以有效控制软岩隧道围岩大变形，但单独使用该方法时围岩变形仍然较大，且稳定性较差，需要尝试联合其他措施共同用于软岩大变形控制.

根据软岩隧道大变形特征，可以将大变形分为挤压变形和松动变形. 超前应力释放是释放挤压变形的重要方式，但是，仅释放挤压变形不对松动变形进行控制，对整个大变形的控制效果并不十分理想. 因此，本文首先提出了“超前应力释放 + 注浆加固 + 加长锚杆”的高地应力软岩大变形控制措施；然后，借助适用于围岩大变形计算的软岩峰后刚度与强度统一劣化模型，基于统一强度准则得到了超前导洞围岩及正洞围岩的弹塑性解；最后，借助Visual studio 2010软件及FLAC3D相关插件实现峰后刚度与强度统一劣化模型的本构开发后，对超前导洞围岩及正洞注浆围岩进行了数值模拟研究.

1. 超前导洞应用存在的问题分析及大变形控制措施

超前导洞常用于高地应力软岩隧道围岩大变形控制，其目的是通过开挖超前导洞使围岩中存储的高地应力得到释放，待开挖隧道正断面时围岩变形量有所下降，达到释放围岩挤压变形的目的. 但目前，在部分工程案例^[9,14-15]中发现，采用超前导洞应力释放后，待扩挖释放层时正洞围岩仍会产生大变形，而且扩挖过程中释放层也容易产生塌方，造成围岩不稳定且影响施工效率. 超前导洞掌子面与扩挖掌子面的间距以及超前导洞的临时支护强度也没有明确的确定方法. 相关工程案例表明^[18-20]，可以将高地应力软岩隧道大变形分为挤压变形和松动变形，挤压变形需要释放，以此来降低作用在支护结构上的形变压力，但是在释放挤压变形的过程中靠近洞壁的部分围岩会进入松动状态，并以松动压力的形式作用在支护结构上，当松动压力过大时支护结构仍然会产生大变形，造成侵限、失效等问题. 因此，可以将超前导洞应力释放后正洞围岩产生的大变形归类为松动变形，释放层不稳定且容易产生塌方也可以证明这一点.

基于上述分析，本文提出一种基于超前应力释放及环向注浆的高地应力软岩大变形支护方式，具体为“超前应力释放 + 环向注浆 + 加长锚杆”，如图1所示. 其施工顺序为：超前导洞开挖与临时支护→环向注浆→释放层扩挖→施作加长锚杆与永久初期支护→浇筑二次衬砌. 采用上述开挖与支护方式时需要注意以下3点：1）导洞开挖后施作临时支护，待其围岩变形结束且稳定后进行注浆加固；2）合理控制超前导洞断面面积，使其既能释放围岩大部分变形，又能保证围岩在注浆加固前保持稳定，且临时支护结构处于安全稳定状态；3）由于高地应力软岩在原始状态下注浆施工较为困难且注浆效果不理想，因此，需要合理控制超前导洞断面面积，使其开挖造成的围岩塑性软化区范围大于正洞断面范围，且大于未注浆条件下正断面开挖完毕后造成的塑性残余区范围.

图 1 “超前导洞 + 环向注浆 + 加长锚杆”支护示意

Figure 1. Support of “advance guide tunnel + circumferential grouting + lengthened anchor rod”

下载: 全尺寸图片幻灯片

2. 理论分析模型

2.1 超前导洞及正洞围岩力学模型

超前导洞及正洞开挖时的假定条件：1）围岩为各向同性均质围岩，超前导洞及正洞埋深足够深，原岩应力场为P₀；2）超前导洞及正洞开挖断面均为圆形，超前导洞开挖半径为r₀，正洞开挖半径为R₀；3）力学模型如图2所示. 图中：r_p、r_b分别为超前导洞开挖造成的围岩塑性软化区半径和塑性残余区半径；p_i为临时支护提供的支护抗力； R₀−r₀为释放层厚度；R_p为由于扩挖造成的围岩塑性软化区半径； R_g为环向注浆半径；R_L为系统锚杆锚杆半径，存在R_L>R_p；p_j 为永久初期支护抗力.

图 2 围岩分区示意

Figure 2. Zoning scheme of surrounding rock

下载: 全尺寸图片幻灯片

2.2 统一强度准则与围岩劣化模型

统一强度准则^[18]充分考虑了中间主应力对材料强度的影响，如式（1a）、（1b）所示.

当σ₂≤0.5(σ₁ + σ₃) + 0.5sin φ(σ₁−σ₃)时：

${\sigma _1} - \frac{{\left( {1 - \sin \;\varphi } \right)\left( {b {\sigma _2} + {\sigma _3}} \right)}}{{\left( {1 + b} \right)\left( {1 + \sin \;\varphi } \right)}} = \frac{{2 c \cos \;\varphi }}{{1 + \sin \;\varphi }}；$

(1a)

当σ₂>0.5(σ₁ + σ₃) + 0.5sin φ(σ₁−σ₃)时：

$\frac{{\left( {{\sigma _1} + b{\sigma _2}} \right)}}{{1 + b}} - \frac{{1 - \sin \;\varphi }}{{1 + \sin \;\varphi }}{\sigma _2} = \frac{{\left( {2c\cos\; \varphi } \right)}}{{1 + \sin \;\varphi }},$

(1b)

式中：σ₁、σ₃分别为最大主应力与最小主应力，σ₂为中间主应力，c、φ分别为材料黏聚力与内摩擦角，b为反映中间主应力σ₂影响的权系数.

令m=(σ₂−σ₃)/(σ₁−σ₃)，将m代入式（1）可得

${\sigma _1} = \frac{N}{M}{\sigma _3} + \frac{k}{M},$

(2)

式中：M、N、k均为计算参数.

当b≤（1 + sin φ）/2时，有$\dfrac{N}{M} = \dfrac{{\left[ {1 + \left( {1 - b} \right)m} \right]\left( {1 + \sin \;\varphi } \right)}}{{1 + \left( {1 - b} \right)m - \left[ {1 + \left( {1 + b} \right)m} \right]\sin \;\varphi }} $，$\dfrac{k}{M} = \dfrac{{2\left( {1 + m} \right)c\cos \;\varphi }}{{1 + \left( {1 - b} \right)m - \left[ {1 + \left( {1 + b} \right)m} \right]\sin \;\varphi }} $；

当b>（1+sin φ）/2时，有$\dfrac{N}{M} = \dfrac{{1 + bm + \left[ {1 + \left( {2 - b} \right)m} \right]\sin \;\varphi }}{{\left( {1 + bm} \right)\left( {1 - \sin \;\varphi } \right)}} $，$ \dfrac{k}{M} = \dfrac{{2\left( {1 + m} \right)c\cos \;\varphi }}{{\left( {1 + bm} \right)\left( {1 - \sin \;\varphi } \right)}} $.

软弱围岩进入塑性阶段后，如前所述的围岩大变形特征以松动变形为主，岩体完整性较差，自身承载力明显下降，松散破碎，处于塑性流动（残余）状态. 传统应变软化模型仅考虑了岩体进入塑性阶段后强度参数的劣化^[21-23]，忽略了塑性区内围岩松散破碎、刚度下降的特点. 因此，本文采用软岩峰后刚度与强度统一劣化模型来反映软岩峰后刚度与强度的统一劣化状态. 如图3所示：当围岩处于弹性阶段时，岩石内部的等效应变累积值ε小于进入塑性软化阶段的临界值ε_p，此时黏聚力初始值c₀、内摩擦角初始值φ₀、弹性模量初始值E₀均可以取峰值，将c₀、φ₀代入式（2）可得到弹性阶段的屈服准则表达式；当岩石内部的等效应变累积值ε达到或大于进入塑性软化阶段的临界值ε_p后，黏聚力、内摩擦角、弹性模量均取软化值c₁、φ₁、E₁，如式（3）所示，将c₁、φ₁代入式（2）便可得到塑性软化阶段的屈服准则表达式（即可得到相应的M₁、N₁、k₁）；当岩石内部的等效应变累积值ε达到或大于进入塑性残余阶段的临界值ε_b后，黏聚力、内摩擦角、弹性模量均取残余值c₂、φ₂、E₂，将c₂、φ₂代入式（2）可得到塑性残余阶段的屈服准则表达式（即可得到相应的M₂、N₂、k₂）.

图 3 岩石参数随塑性应变的劣化规律

Figure 3. Degradation of rock parameters with plastic strain

下载: 全尺寸图片幻灯片

锚杆锚固后的围岩通常被视为等效复合岩体，借助参数等效原则通过引入锚杆密度因子^[24]得到复合岩体刚度与强度参数，当锚固围岩内部应变累积达到一定程度后会进入塑性软化（残余）阶段，将塑性锚固围岩参数代入式（3）后可得到不同阶段的屈服准则表达式.

当ε_p≤ε<ε_b时，有

$\left\{ {\begin{array}{*{20}{l}} {{c_1} = {c_0} - {M_c}\left( {\varepsilon - {\varepsilon _{\text{p}}}} \right)}, \\ {{\varphi _1} = {\varphi _0} - {M_\varphi }\left( {\varepsilon - {\varepsilon _{\text{p}}}} \right)}, \\ {{E_1} = {E_0} - {M_E}\left( {\varepsilon - {\varepsilon _{\text{p}}}} \right)} , \end{array}} \right.$

(3)

式中： M_c、M_φ、M_E分别为黏聚力、内摩擦角、弹性模量的软化模量.

2.3 塑性围岩扩容规律

当塑性区围岩内的等效应变小于ε_b时，流动法则如式（4）所示. 假设：岩石塑性变形服从非关联流动法则，且令塑性势函数g等于屈服函数f，最后将f中的内摩擦角替换成剪胀角ψ即可得到g的表达式^[25].

${\varepsilon _{{\text{rp}}}} + {\eta _1}{\varepsilon _{{\text{θp}}}} = 0 \text{，}$

(4)

式中：η₁为扩容系数，ε_rp与ε_θp分别为塑性软化区内径向累积应变、环向累积应变.

由塑性位势理论可知：

${\mathrm{d}}{\varepsilon _{{\mathrm{ij}} {\text{-}} {\mathrm{p}}}} = {\mathrm{d}}\lambda \frac{{\partial g}}{{\partial {\sigma _{{\mathrm{ij}} {\text{-}} {\mathrm{p}}}}}},$

(5)

式中：dε_ij-p为塑性应变增量；σ_ij-p为应力张量；dλ为与塑性势函数相关联的比例系数，dλ≥0.

根据式（2）及式（5）有

$\left\{ {\begin{array}{*{20}{c}} {{\mathrm{d}}{\varepsilon _{{\text{θp}}}} = {\mathrm{d}}\lambda }, \\ {{\mathrm{d}}{\varepsilon _{{\text{rp}}}} = - {\mathrm{d}}\lambda \dfrac{{{N_{\psi }}}}{{{M_{{\psi }}}}}} , \end{array}} \right.$

(6)

式中：当b≤(1 + sin φ)/2时，N_ψ=[1 + (1−b)m](1−sin ψ)，M_ψ=1 + (1−b)m−[1 + (1−b)m]sin ψ；当b≥(1 + sin φ)/2时，N_ψ=1 + bm + [1 + (2−b)m]sin ψ，M_ψ=(1 + bm)(1−sin ψ).

定义扩容系数η₁为最小塑性主应变与最大塑性主应变之比，根据以上推导可以得到η₁=N_ψ/M_ψ.

当塑性区内的部分锚固围岩等效应变大于ε_b时，这部分锚固围岩便进入塑性残余状态^[24]，此时考虑扩容的流动法则为

${\varepsilon _{{\text{rb}}}} + {\eta _2}{\varepsilon _{{\text{θb}}}} = 0,$

(7)

式中：ε_rb与ε_θb分别为塑性残余区内径向累积应变、环向累积应变；η₂为塑性残余区的扩容系数，可令η₂=1 + μ，μ多介于0.3～0.5，因此，η₂取1.3～1.5.

3. 围岩弹塑性分析

3.1 基本方程

平衡微分方程（不计体力）为

$\frac{{{\text{d}}{\sigma _{\text{r}}}}}{{{\text{d}}r}} + \frac{{{\sigma _{\text{r}}} - {\sigma _{\text{θ }}}}}{r} = 0；$

(8a)

几何方程为

$\left\{ {\begin{array}{*{20}{c}} {{\varepsilon _{\text{r}}} = \dfrac{{{\text{d}}u}}{{{\text{d}}r}}}, \\ {{\varepsilon _{\text{θ }}} = \dfrac{u}{r}} ; \end{array}} \right.$

(8b)

本构方程（平面应变）为

$\left\{ {\begin{array}{*{20}{c}} {{\varepsilon _{\text{r}}} = \dfrac{{1 - {v^2}}}{E}\left( {{\sigma _{\text{r}}} - \dfrac{v}{{1 - v}}{\sigma _{\text{θ }}}} \right)}, \\ {{\varepsilon _{\text{θ }}} = \dfrac{{1 - {v^2}}}{E}\left( {{\sigma _{\text{θ }}} - \dfrac{v}{{1 - v}}{\sigma _{\text{r}}}} \right)} , \end{array}} \right.$

(8c)

式中：σ_r为围岩径向应力，σ_θ为围岩环向应力；ɛ_r为围岩径向总应变，ɛ_θ为围岩环向总应变，u为围岩径向位移；r为开挖隧道围岩极径；v为泊松比，E为弹性模量.

3.2 超前导洞开挖后围岩弹塑性分析

3.2.1 弹性区

根据拉梅应力计算式以及弹性区边界条件^[25]：r→∞时，径向应力σ_re=P₀，超前导洞围岩弹塑性交界面处，即r=r_p时，径向应力σ_re=σ_r=(2P₀M−k)/(M + N)，可得弹性区围岩径向应力σ_re、环向应力σ_θe为

$\left\{ {\begin{array}{*{20}{l}} {{\sigma _{{\text{re}}}} = {P_0} - \left( {{P_0} - {\sigma _{\text{r}}}} \right) {{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^2}} ，\\ {{\sigma _{{\text{θe}}}} = {P_0} + \left( {{P_0} - {\sigma _{\text{r}}}} \right) {{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^2}}. \end{array}} \right.$

(9)

根据式（8）、（9）可得弹性区围岩径向应变ε_re、环向应变ε_θe分别为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{{\text{re}}}} = \dfrac{{{\sigma _{\text{r}}} - {P_0}}}{{2G}}{{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^2}}, \\ {{\varepsilon _{{\text{θe}}}} = \dfrac{{{P_0} - {\sigma _{\text{r}}}}}{{2G}}{{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^2}}, \end{array}} \right.$

(10)

式中：G为剪切模量，G=E/[2(1 + v)].

将围岩残余强度值c₂、φ₂代入式（2）求得M₂、N₂、k₂，求得超前导洞塑性软化区半径r_p为

${r_{\text{p}}} = {r_0} {\left\{ {\frac{{\left( {{M_2} + {N_2}} \right)\left[ {{p_{\text{i}}}\left( {{M_2} - {N_2}} \right) - {k_2}} \right]}}{{\left( {2{P_0} - {M_2}} \right)\left( {{M_2} - {N_2}} \right) - {k_2}\left( {{M_2} + {N_2}} \right)}}} \right\}^{\frac{{{M_2}}}{{{M_2} - {N_2}}}}},$

(11)

式中：p_i=p_f + p_ts，p_f为超前导洞掌子面提供的虚拟支护力；p_ts为喷射混凝土提供的支护抗力.

3.2.2 塑性软化区

塑性软化区内围岩应力满足式（2）及式（8a），联立两式及边界条件：r=r_p时，σ_rp=σ_re，可得软化区环向、径向应力为

$\left\{ {\begin{array}{*{20}{l}} {\sigma _{{\text{rp}}}} = \left( {\dfrac{{2{M_1}{P_0} - {k_1}}}{{{M_1} + {N_1}}} - \dfrac{{{k_1}}}{{{M_1} - {N_1}}}} \right){{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{\frac{{{M_1} - {N_1}}}{{{M_1}}}}} +\\ \quad \dfrac{{{k_1}}}{{{M_1} - {N_1}}}， \\ {\sigma _{{\text{θp}}}} = \dfrac{{{N_1}}}{{{M_1}}}\left( {\dfrac{{2{M_1}{P_0} - {k_1}}}{{{M_1} + {N_1}}} - \dfrac{{{k_1}}}{{{M_1} - {N_1}}}} \right){{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{\frac{{{M_1} - {N_1}}}{{{M_1}}}}} +\\ \quad \dfrac{{{k_1}}}{{{M_1} - {N_1}}} . \end{array}} \right.$

(12)

塑性软化区内围岩径向总应变ɛ_r和环向总应变ɛ_θ为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{\text{r}}} = {\varepsilon _{{\text{re}}}}\left| {_{r = {r_{\text{p}}}}} \right. + {\varepsilon _{{\text{rp}}}}}, \\ {{\varepsilon _{\text{θ }}} = {\varepsilon _{{\text{θe}}}}\left| {_{r = {r_{\text{p}}}}} \right. + {\varepsilon _{{\text{θp}}}}} . \end{array}} \right.$

(13)

根据式（4）、（8）及边界条件r=r_p时，u=r_p(P₀−σ_R)/(2G)，可得到塑性软化区围岩位移为

${u_{{\text{rp}}}} = r\frac{{{P_0} - {\sigma _{\text{r}}}}}{{\left( {{\eta _1} + 1} \right){G_1}}}{\left( {\frac{{{r_{\text{p}}}}}{r}} \right)^{{\eta _1} + 1}} + r\frac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _r}} \right)}}{{2{G_1}\left( {{\eta _1} + 1} \right)}}.$

(14)

根据式（8b）、（14）可以得到

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{{\text{rp}}}} = - {\eta _1}\dfrac{{{P_0} - {\sigma _{\text{r}}}}}{{\left( {{\eta _1} + 1} \right){G_1}}}{{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{{\eta _1} + 1}} + \dfrac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _{\text{r}}}} \right)}}{{2{G_1}\left( {{\eta _1} + 1} \right)}}}, \\ {{\varepsilon _{{\text{θp}}}} = \dfrac{{{P_0} - {\sigma _{\text{r}}}}}{{\left( {{\eta _1} + 1} \right)}}{{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{{\eta _1} + 1}} + \dfrac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _{\text{r}}}} \right)}}{{2{G_1}\left( {{\eta _1} + 1} \right)}}} . \end{array}} \right.$

(15)

3.2.3 塑性残余区

塑性残余区内围岩应力满足式（2）、（8a），联立两式及边界条件r=r₀时，σ_rb=p_i（t），可得塑性残余区围岩径向应力σ_rb、环向应力σ_θb为

$\left\{ {\begin{array}{*{20}{l}} {{\sigma _{{\text{rb}}}} = \left( {{p_{\text{i}}} - \dfrac{{{k_2}}}{{{M_2} - {N_2}}}} \right) {{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{\frac{{{M_2} - {N_2}}}{{{M_2}}}}} + \dfrac{{{k_2}}}{{{M_2} - {N_2}}}}, \\ {{\sigma _{{\text{θb}}}} = \dfrac{{{N_2}}}{{{M_2}}}\left( {{p_{\text{i}}} - \dfrac{{{k_2}}}{{{M_2} - {N_2}}}} \right) {{\left( {\dfrac{{{r_{\text{p}}}}}{r}} \right)}^{\frac{{{M_2} - {N_2}}}{{{M_2}}}}} + \dfrac{{{k_2}}}{{{M_2} - {N_2}}}} . \end{array}} \right.$

(16)

塑性残余区围岩径向总应变 ɛ_r和环向总应变 ɛ_θ为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{\text{r}}} = {\varepsilon _{{\text{rp}}}}\left| {_{r = {r_{\text{b}}}}} \right. + {\varepsilon _{{\text{rb}}}}}, \\ {{\varepsilon _{\text{θ}}} = {\varepsilon _{{\text{θp}}}}\left| {_{r = {r_{\text{b}}}}} \right. + {\varepsilon _{{\text{θb}}}}} . \end{array}} \right.$

(17)

根据式（7）、（8b）及r=r_b时的边界条件可以得到软化区位移为

$\begin{split} & {u_{{\text{rb}}}} = r\frac{{{P_0} - {\sigma _{\text{r}}}}}{{\left( {{\eta _2} + 1} \right){G_2}}}{\left( {\frac{{{r_{\text{p}}}}}{{{r_{\text{b}}}}}} \right)^{{\eta _1} + 1}}{\left( {\frac{{{r_{\text{b}}}}}{r}} \right)^{{\eta _2} + 1}} + r {\left( {\frac{{{r_{\text{p}}}}}{{{r_{\text{b}}}}}} \right)^{{\eta _1} + 1}} \times \\ &\quad \frac{{\left( {{\eta _2} - {\eta _1}} \right)\left( {{P_0} - {\sigma _{\text{r}}}} \right)}}{{\left( {{\eta _2} + 1} \right)\left( {{\eta _1} + 1} \right){G_2}}} + r\frac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _{\text{r}}}} \right)}}{{2{G_2}\left( {{\eta _1} + 1} \right)}}, \end{split}$

(18)

式中：${r_{\text{b}}} = {{{r_{\text{p}}}} / {{{\left[ {\dfrac{{\left( {{c_0} - {c_2}} \right)\left( {{\eta _1} + 1} \right) {G_2}}}{{{M_{\text{c}}}\left( {{P_0} - {\sigma _{\text{r}}}} \right)}} + 1} \right]}^{\frac{1}{{{\eta _1} + 1}}}}}} $.

3.3 扩挖后围岩弹塑性分析

扩挖后的围岩分区（由外到内）：普通弹塑性区→普通塑性软化区→注浆加固区（注浆区半径为R_g），其中普通弹性区和普通塑性软化区应力及位移如式（9）、（12）及式（14）所示，计算时只需把式中r₀替换为R₀，p_i替换为p_f.

注浆加固区内围岩应力满足式（2）、（8a），联立两式及边界条件可得注浆区围岩径向应力σ_rg、环向应力σ_θg为

$\left\{ {\begin{array}{*{20}{l}} {{\sigma _{{\text{rg}}}} = \left( {{p_{\text{f}}} - \dfrac{{{k_{\text{g}}}}}{{{M_{\text{g}}} - {N_{\text{g}}}}}} \right) {{\left( {\dfrac{{{R_0}}}{r}} \right)}^{\frac{{{M_{\text{g}}} - {N_{\text{g}}}}}{{{M_{\text{g}}}}}}} + \dfrac{{{k_{\text{g}}}}}{{{M_{\text{g}}} - {N_{\text{g}}}}}}, \\ {{\sigma _{{\text{θg}}}} = \dfrac{{{N_{\text{g}}}}}{{{M_{\text{g}}}}}\left( {{p_{\text{f}}}\left( t \right) - \dfrac{{{k_{\text{g}}}}}{{{M_{\text{g}}} - {N_{\text{g}}}}}} \right) {{\left( {\dfrac{{{R_0}}}{r}} \right)}^{\frac{{{M_{\text{g}}} - {N_{\text{g}}}}}{{{M_{\text{g}}}}}}} + \dfrac{{{k_{\text{g}}}}}{{{M_{\text{g}}} - {N_{\text{g}}}}}} , \end{array}} \right.$

(19)

式中：M_g、N_g、k_g分别为围岩注浆后的计算参数，可将围岩注浆黏聚力c_g、内摩擦角φ_g代入式（2）求得.

注浆区围岩径向总应变 ɛ_r和环向总应变 ɛ_θ为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{\text{r}}} = {\varepsilon _{{\text{rp}}}}\left| {_{r = {r_{\text{g}}}}} \right. + {\varepsilon _{{\text{rg}}}}}, \\ {{\varepsilon _{\text{θ}}} = {\varepsilon _{{\text{θp}}}}\left| {_{r = {r_{\text{g}}}}} \right. + {\varepsilon _{{\text{θg}}}}} . \end{array}} \right.$

(20)

根据式（7）、（8b）及边界条件可得

$\begin{split} & {u_{{\text{rb}}}} = r\frac{{{P_0} - {\sigma _{\text{R}}}}}{{\left( {{\eta _2} + 1} \right){G_{\text{g}}}}}{\left( {\frac{{{R_{\text{p}}}}}{{{R_{\text{g}}}}}} \right)^{{\eta _1} + 1}}{\left( {\frac{{{R_{\text{g}}}}}{r}} \right)^{{\eta _2} + 1}} + r{\left( {\frac{{{R_{\text{p}}}}}{{{R_{\text{g}}}}}} \right)^{{\eta _1} + 1}}\times \\ &\quad \frac{{\left( {{\eta _2} - {\eta _1}} \right)\left( {{P_0} - {\sigma _{\text{R}}}} \right)}}{{\left( {{\eta _2} + 1} \right)\left( {{\eta _1} + 1} \right){G_{\text{g}}}}} + r\frac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _{\text{R}}}} \right)}}{{2{G_{\text{g}}}\left( {{\eta _1} + 1} \right)}}, \end{split}$

(21)

式中：R_p、σ_R为直接开挖半径为R₀的圆形隧道时产生的塑性软化区半径与弹塑性交界面处径向应力，R_p如式（22）.

${R_{\text{p}}} = {R_0}{\left\{ {\frac{{\left( {{M_2} + {N_2}} \right)\left( {{p_{\text{f}}}\left( {{M_2} - {N_2}} \right) - {k_2}} \right)}}{{\left( {2{P_0} - {M_2}} \right)\left( {{M_2} - {N_2}} \right) - {k_2}\left( {{M_2} + {N_2}} \right)}}} \right\}^{\frac{{{M_2}}}{{{M_2} - {N_2}}}}}.$

(22)

3.4 永久初期支护施作后围岩弹塑性分析

3.4.1 锚固弹性区

锚固弹性区围岩可以看成内外受径向应力的圆环，p₂为锚固区与非锚固区交界处（r=R_L）的径向应力，p₁为锚固围岩弹塑性交界面处（r=R_sp）的径向应力，因此锚固围岩应力表达式为

$\left\{ \begin{array}{*{20}{c}}\sigma_{\text{re}}=\dfrac{\sigma_{\text{R}_{\text{L}}}R_{\text{L}}^{\text{2}}-\sigma_{\text{R}_{\text{p}}}R_{\text{sp}}^{\text{2}}}{2\left(R_{\text{L}}^{\text{2}}-R_{\text{sp}}^{\text{2}}\right)}-\dfrac{\left(\sigma_{\text{R}_{\text{L}}}-\sigma_{\text{R}_{\text{p}}}\right)R_{\text{L}}^{\text{2}}R_{\text{sp}}^{\text{2}}}{r^2\left(R_{\text{L}}^{\text{2}}-R_{\text{sp}}^{\text{2}}\right)}， \\ \sigma_{\text{θe}}=\dfrac{\sigma_{\text{R}_{\text{L}}}R_{\text{L}}^{\text{2}}-\sigma_{\text{R}_{\text{p}}}R_{\text{sp}}^{\text{2}}}{2\left(R_{\text{L}}^{\text{2}}-R_{\text{sp}}^{\text{2}}\right)}+\dfrac{\left(\sigma_{\text{R}_{\text{L}}}-\sigma_{\text{R}_{\text{p}}}\right)R_{\text{L}}^{\text{2}}R_{\text{sp}}^{\text{2}}}{r^2\left(R_{\text{L}}^{\text{2}}-R_{\text{sp}}^{\text{2}}\right)}，\end{array}\right.$

(23)

式中：R_sp、${\sigma _{{{\text{R}}_{\text{L}}}}} $、${\sigma _{{{\text{R}}_{\text{sp}}}}} $、${\sigma _{{{\text{R}}_{\text{p}}}}} $分别为直接开挖半径为R₀的圆形隧道，且仅有系统锚杆支护时产生的围岩塑性软化区半径、r=R_L、r=R_sp、r=R_p处的围岩径向应力，其计算公式可借用式（22）.

根据本构方程，减去应力分量中的原岩应力成分后，由式（8c）、（19）可得锚固围岩弹性区应变分布为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{{\text{re}}}} = \dfrac{1}{{2{G_{\text{s}}}}}\left[ {\left( {2C - {P_0}} \right)\left( {1 - 2v} \right) - \dfrac{A}{{{r^2}}}} \right]}, \\ {{\varepsilon _{{\text{θe}}}} = \dfrac{1}{{2{G_{\text{s}}}}}\left[ {\left( {2C - {P_0}} \right)\left( {1 - 2v} \right) + \dfrac{A}{{{r^2}}}} \right]} , \end{array}} \right.$

(24)

式中：$ A = \dfrac{{\left( {{\sigma _{{{\text{R}}_{\text{L}}}}} - {\sigma _{{{\text{R}}_{\text{p}}}}}} \right)R_{\text{L}}^{\text{2}}R_{{\text{sp}}}^{\text{2}}}}{{R_{\text{L}}^{\text{2}} - R_{{\text{sp}}}^{\text{2}}}} $，$ C = \dfrac{{{\sigma _{{{\text{R}}_{\text{L}}}}}R_{\text{L}}^{\text{2}} - {\sigma _{{{\text{R}}_{\text{p}}}}}R_{{\text{sp}}}^{\text{2}}}}{{2\left( {R_{\text{L}}^{\text{2}} - R_{{\text{sp}}}^{\text{2}}} \right)}} $.

3.4.2 锚固塑性软化区

锚固塑性软化区内围岩应力满足式（2）、（8a），联立两式及边界条件可得锚固塑性软化区径向应力σ_r-sp、环向应力σ_θ-sp为

$\left\{ {\begin{array}{*{20}{l}} {\sigma _{{\text{r-sp}}}}\left( {\dfrac{{2{M_{{\text{1s}}}}{P_0} - {k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} + {N_{{\text{1s}}}}}} - \dfrac{{{k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}}} \right) {{\left( {\dfrac{{{R_{{\text{sp}}}}}}{r}} \right)}^{\frac{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}}{{{M_{{\text{1s}}}}}}}} + \\ \quad \dfrac{{{k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}} ,\\ {\sigma _{{\text{θ-sp}}}} = \dfrac{{{N_{{\text{1s}}}}}}{{{M_{{\text{1s}}}}}}\left( {\dfrac{{2{M_{{\text{1s}}}}{P_0} - {k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} + {N_{{\text{1s}}}}}} - \dfrac{{{k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}}} \right) \times \\ \quad {{\left( {\dfrac{{{R_{{\text{sp}}}}}}{r}} \right)}^{\frac{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}}{{{M_{{\text{1s}}}}}}}} + \dfrac{{{k_{{\text{1s}}}}}}{{{M_{{\text{1s}}}} - {N_{{\text{1s}}}}}} . \end{array}} \right.$

(25)

式中：M_1s、N_1s、k_1s为系统锚杆加固后的围岩计算参数.

根据3.2.2小节所述塑性软化区内的应变组成，根据式（4）、（8）及r=R_sp时的边界条件可得锚固塑性软化区围岩的位移u_r-sp为

$\begin{split} & {u_{{\text{r-sp}}}} = r\frac{A}{{\left( {{\eta _1} + 1} \right){G_{{\text{1s}}}}R_{{\text{sp}}}^{\text{2}}}}{\left( {\frac{{{R_{{\text{sp}}}}}}{r}} \right)^{{\eta _1} + 1}} + \\ &\quad \frac{r}{{2{G_{{\text{1s}}}}}}\left( {\frac{{{\eta _1} - 1}}{{{\eta _1} + 1}} \frac{A}{{R_{{\text{sp}}}^{\text{2}}}} + \left( {2C - {P_0}} \right)\left( {1 - 2v} \right)} \right) . \end{split}$

(26)

3.4.3 锚固注浆区

锚固塑性注浆区内围岩应力满足式（2）、（8a），联立两式及边界条件r=R₀时，σ_r-gs=p_j，可得锚固注浆区径向应力σ_r-gs、环向应力σ_θ-gs为

$\left\{ {\begin{array}{*{20}{l}} {\sigma _{{\mathrm{r}} {\text{-}} {\mathrm{gs}}}} = \left( {\dfrac{{2{M_{{\mathrm{gs}}}}{P_0} - {k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} + {N_{{\mathrm{gs}}}}}} - \dfrac{{{k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}}} \right) \times\\ \quad {{\left( {\dfrac{{{R_0}}}{r}} \right)}^{\frac{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}}}}}} + \dfrac{{{k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}}， \\ {\sigma _{ {\text{θ-}}{\mathrm{ gs}}}} = \dfrac{{{N_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}}}}\left( {\dfrac{{2{M_{{\mathrm{gs}}}}{P_0} - {k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} + {N_{{\mathrm{gs}}}}}} - \dfrac{{{k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}}} \right)\times \\ \quad {{\left( {\dfrac{{{R_0}}}{r}} \right)}^{\frac{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}}}}}} + \dfrac{{{k_{{\mathrm{gs}}}}}}{{{M_{{\mathrm{gs}}}} - {N_{{\mathrm{gs}}}}}} , \end{array}} \right.$

(27)

式中：M_gs、N_gs、k_gs为同时考虑系统锚杆与注浆加固后的围岩计算参数.

锚固注浆区围岩径向总应变ɛ_r和环向总应变ɛ_θ为

$\left\{ {\begin{array}{*{20}{l}} {{\varepsilon _{\text{r}}} = {\varepsilon _{{\text{rp}}}}\left| {_{r = {r_{{\text{sp}}}}}} \right. + {\varepsilon _{{\text{r-gs}}}}}, \\ {{\varepsilon _{\text{θ}}} = {\varepsilon _{{\text{θp}}}}\left| {_{r = {r_{{\text{sp}}}}}} \right. + {\varepsilon _{{\text{θ-gs}}}}} . \end{array}} \right.$

(28)

根据式（4）、（8b）及边界条件可得锚固注浆区围岩的位移u_r-gs为

$\begin{split} & {u_{{\text{r-gs}}}} = r\frac{{{P_0} - {\sigma _{\text{R}}}}}{{{\eta _2} + 1}}{\left( {\frac{{{R_{{\text{sp}}}}}}{{{R_{\text{g}}}}}} \right)^{{\eta _1} + 1}}{\left( {\frac{{{R_{\text{g}}}}}{r}} \right)^{{\eta _2} + 1}} +{\left( {\frac{{{R_{{\text{sp}}}}}}{{{R_{\text{g}}}}}} \right)^{{\eta _1} + 1}}\times \\ &\quad r\frac{{\left( {{\eta _2} - {\eta _1}} \right)\left( {{P_0} - {\sigma _{\text{R}}}} \right)}}{{\left( {{\eta _2} + 1} \right)\left( {{\eta _1} + 1} \right)}} + r\frac{{\left( {{\eta _1} - 1} \right)\left( {{P_0} - {\sigma _{\text{R}}}} \right)}}{{2{G_{{\mathrm{gs}}}}\left( {{\eta _1} + 1} \right)}} , \end{split}$

(29)

式中：G_gs为同时考虑系统锚杆与注浆加固后的围岩计算参数，R_sp如式（30）.

$\begin{split} & {R_{{\text{sp}}}} = \\ &\quad {R_0}{\left\{ {\frac{{\left[ {{p_{\text{j}}}\left( {{M_{{\text{2s}}}} - {N_{{\text{2s}}}}} \right) - {k_{{\text{2s}}}}} \right]\left( {{M_{{\text{2s}}}} + {N_{{\text{2s}}}}} \right)}}{{\left( {2{P_0} - {M_{{\text{2s}}}}} \right)\left( {{M_{{\text{2s}}}} - {N_{{\text{2s}}}}} \right) - {k_{{\text{2s}}}}\left( {{M_{{\text{2s}}}} + {N_{{\text{2s}}}}} \right)}}} \right\}^{\frac{{{M_{{\text{2s}}}}}}{{{M_{{\text{2s}}}} - {N_{{\text{2s}}}}}}}}, \end{split}$

(30)

式中：p_j=p_f + p_ts，p_ts为喷射混凝土提供的支护抗力，M_2s、N_2s、k_2s为考虑系统锚杆加固后的围岩残余计算参数.

4. 算例验证及数值模拟

某隧道工程所穿越地层主要为志留系下统页岩，洞身围岩大多为微、中风化页岩，本文选取该隧道深埋段K34 + 500～K34 + 575作为研究对象. 采用面积代换的方式测得隧道开挖等价半径为5.25 m，计算过程中超前导洞开挖半径为3.5 m，超前导洞临时支护采用C25喷射混凝土（厚度为20 cm） + HW200钢拱架（间距为1.0 m），正洞永久初期支护包括全黏结式锚杆（长度为6.5 m） + C25喷射混凝土（厚度为20 cm） + HW200钢拱架（间距为0.8 m），系统锚杆看作是对围岩参数的提高^[21]. 借助FLAC3D有限差分软件中的二次开发功能，以Visual Studio 2010为开发平台，编写软岩峰后劣化黏弹塑性本构模型，对该隧道进行数值模拟分析，模型尺寸为100 m × 100 m × 60 m，计算参数如表1所示.

表 1 围岩计算参数

Table 1. Calculation parameters of surrounding rock

条件	E/MPa	c/MPa	φ/（º）	v	容重 γ/kN·m⁻³	M_c/MPa	M_φ/ MPa	P₀/MPa
初始值	2 800	3.2	32	0.2	20.5	400	1 800	12.2
残余值	1 000	1.2	20
注浆后	2 500	2.6	28

下载: 导出CSV

| 显示表格

4.1 算例验证

由图4可知：在高地应力软岩中开挖隧道（无超前导洞、无注浆加固）时，围岩最大位移为19.45 cm，此计算结果接近施工现场实际监测的20.82 cm；开挖后的围岩位移主要集中在前30 d，前30 d的围岩变形占总位移的89.6%；超前导洞开挖及支护（临时支护）后的围岩最大位移为11.78 cm，大约需要30 d的时间达到稳定状态；无注浆时正洞围岩最终位移为12.38 cm，有注浆时正洞围岩最终位移为3.21 cm，二者相比后者下降了74%.

图 4 围岩变形控制效果对比

Figure 4. Comparison of deformation control effects of surrounding rock

下载: 全尺寸图片幻灯片

4.2 数值模拟分析

图5为高地应力软岩中直接开挖隧道时围岩的位移、最小主应力、塑性区分布云图. 由图可知：1）隧道围岩最大位移为19.53 cm，塑性软化区半径R_p为13.6 m，数值模拟结果与理论计算结果均相近，说明理论计算结果可靠，围岩最小主应力的最大值位于弹塑性交界面处，约为18.98 MPa；2）当r>R_p时，随着r的增大，围岩最小主应力逐渐下降，最终最小主应力接近P₀；当r<R_p时，随着r的减小，围岩最小主应力同样下降，当r=R₀时，围岩最小主应力仅有4.46 MPa.

图 5 无导洞无注浆时正洞围岩数值模拟计算结果

Figure 5. Numerical simulation results of main tunnel surrounding rock without guide tunnel and grouting

下载: 全尺寸图片幻灯片

图6为超前导洞围岩位移、最小主应力、塑性区分布云图，由图可知：超前导洞围岩最大位移为12.56 cm，r_p为12 m，与图5中的数据相比，围岩位移有大幅下降，塑性区范围也有明显下降，可见超前导洞开挖半径是影响围岩位移的重要因素；超前导洞围岩最小主应力为18.67 MPa，与图5中的数据相比有小幅下降，但下降幅值较小.

图 6 超前导洞围岩数值模拟计算结果

Figure 6. Numerical simulation results of surrounding rock of advance guide tunnel

下载: 全尺寸图片幻灯片

图7为注浆加固后对缓释层开挖时的正洞围岩位移、最小主应力、塑性区分布云图，由图可知：围岩最大位移为1.92 cm，R_p为6.25 cm，围岩最小主应力最大值同样位于弹塑性交界面处，但塑性区范围内围岩的承载能力有所增加，当r=R₀时，围岩最小主应力为7.83 MPa.

图 7 正洞围岩（注浆后）数值模拟计算结果

Figure 7. Numerical simulation results of main tunnel surrounding rock （after grouting）

下载: 全尺寸图片幻灯片

5. 影响因素分析

5.1 软化模量（M_c、M_φ、M_E）的影响

图8（a）为软化模量（M_c、M_φ）取不同值时的超前导洞围岩位移. 由图可知：1）理论计算与数值模拟结果均显示软化模量M_c、M_φ取值越小时（超前导洞）围岩位移越小；2） M_c、M_φ取值越大时，M_c、M_φ每增加100 MPa/200 MPa时得到的围岩位移增长幅值更大，以M_c为例分析可知，M_c从400 MPa增加到500 MPa时围岩位移增加了4.57 cm，M_c从500 MPa增加到600 MPa时围岩位移增加了10.01 cm.

图 8 软化模量M_c、M_φ、M_E的影响

Figure 8. Influence of softening modulus M_c、M_φ and M_E

下载: 全尺寸图片幻灯片

图8（b）为为软化模量（M_E）取不同值时的超前导洞围岩位移. 由图可知：1）超前导洞围岩位移大小与弹性模量（刚度参数）劣化密不可分，考虑弹性模量劣化后计算得到的（超前导洞）围岩位移值迅速增大；2） M_E取值越大时围岩位移值越大；3）当M_E越大时，M_E每增加250 MPa时得到的围岩位移增长幅值更大. 如当M_E为2 000 MPa时围岩位移为6.83 cm，当M_E增加到2 250、2 500、2 750、3 000 MPa时，围岩位移分别增加了30.7%、89.8%、158.2%、251.7%.

5.2 注浆参数的影响

图9为注浆参数取不同值时正洞围岩位移大小. 由图可知：1）注浆参数c_g、φ_g、扩挖后围岩注浆的弹性模量E_g取值越大时正洞围岩的位移越小；2）注浆参数c_g、φ_g取值越小时，c_g（φ_g）每增加0.2 MPa（2°）时正洞围岩位移下降幅值越大；3）与c_g、φ_g相比，注浆弹性模量E_g的取值对正洞围岩位移的影响更明显；4）当E_g取值越小时，每增加200 MPa时得到的正洞围岩位移下降幅值越大，例如E_g从1 600 MPa增加至1 800 MPa时围岩位移下降了约1.6 cm，E_g从2 200 MPa增加至2 400 MPa时围岩位移下降了约0.4 cm.

图 9 注浆参数的影响

Figure 9. Influence of grouting parameters

下载: 全尺寸图片幻灯片

5.3 超前导洞半径及两掌子面间距的影响分析

图10、11分别为超前导洞开挖半径、两掌子面（超前导洞掌子面及正洞缓释层掌子面）间距取不同值时的围岩位移. 由图中可知：1）超前导洞开挖半径r₀越大时初始地应力释放越充分，正洞围岩位移越小；2）当r₀越小时，r₀每增加0.5 m时得到的正洞围岩位移下降幅值越大，例如：当r₀取2.0、3.0、4.0 m时正洞围岩位移分别为28.91、16.41、12.38 m；3）当两掌子面间距越大时，超前应力释放越充分，正洞围岩位移越小；4）两掌子面间距对正洞围岩位移的影响明显小于超前导洞开挖半径.

图 10 掌子面间距的影响

Figure 10. Influence of distance between tunnel faces

下载: 全尺寸图片幻灯片

图 11 超前导洞半径r₀的影响

Figure 11. Influence of radius r₀ of advance guide tunnel

下载: 全尺寸图片幻灯片

6. 结　论

本文根据软岩大变形特征对超前导洞在用于软岩大变形控制时存在的问题进行了分析，提出了基于“超前应力释放 + 注浆加固 + 加长锚杆”的软岩大变形控制措施. 借助软岩峰后刚度与强度劣化模型对超前导洞及正洞围岩进行弹塑性解析，并对超前导洞及正洞围岩进行数值模拟研究，主要得出以下结论：

1）超前导洞释放了围岩挤压变形，岩体松动破碎是造成正洞围岩变形过大及释放层岩体失稳的主要原因，注浆加固能够有效控制松动变形.

2）与未注浆相比，注浆后正洞围岩位移下降了66.8%，注浆加固控制了围岩松动变形，并改善了围岩应力分布，减小了塑性区围岩分布范围，提高了围岩承载力.

3）软化模量M_c、M_E、M_φ取值越大时，超前导洞围岩位移越大，反之超前导洞围岩位移越小；注浆参数c_g、φ_g、E_g取值越大时，正洞围岩位移越小，反之正洞围岩位移越大.

4）超前导洞开挖半径r₀及两掌子面间距越大时，超前应力释放越充分，但增大r₀对超前应力释放的效果更理想.

图 1 实时数据查询

Figure 1. Real-time data query

下载: 全尺寸图片幻灯片

图 2 历史数据查询

Figure 2. Historical data query

下载: 全尺寸图片幻灯片

图 3 城市道路轻型汽车行驶工况构建基本过程

Figure 3. Basic process of constructing the automobile driving cycle on urban roads

下载: 全尺寸图片幻灯片

图 4 私家车、出租车、公务车出行特征分析

Figure 4. Travel feature analysis of private，taxi and financing vehicles

下载: 全尺寸图片幻灯片

图 5 城市道路轻型汽车行驶工况

Figure 5. Light-duty vehicles condition on urban roads

下载: 全尺寸图片幻灯片

表 1 CATC信息化系统菜单功能

Table 1. Menu function of CATC information system

一级菜单	二级菜单
车辆监控	全屏地图检测历史轨迹查询与下载历史轨迹对比实时数据查询
数据管理	历史数据查询与下载文件下载
车辆与终端	车辆信息管理终端管理车辆与终端绑定车辆与终端快速注册 SIM卡管理 SIM卡流量流水查询
城市监控	全国城市监控各城市监控详情单车数据监控

下载: 导出CSV

表 2 部分数据格式要求

Table 2. Data format requirements

数据表示内容	数据类型	描述及要求
VIN	STRING	识别码是识别的唯一标识，由 17 位字码构成，不足 17 位后补“\0”，例如车架号、产品条码、SIM卡号等
经度	DWORD	以度为单位的经度值 × 10⁶，精确到百万分之一度
纬度	DWORD	以度为单位的纬度值 × 10⁶，精确到百万分之一度
方向	WORD	0～360^o，正北为 0，顺时针
速度	WORD	0～220 km/h，最小计量单元：0.1 km/h
CAN 车速	WORD	0～220 km/h，最小计量单元：0.1 km/h，“0xFF，0xFE”表示异常，“0xFF，0xFF”表示无效
里程	DWORD	GPS 车速积分累计里程，有效值范围 0～429 496 729.5 km，最小计量单元： 0.1 km
海拔	WORD	单位：m
电机转速	WORD	0～65 532 r/min（数值偏移量 20 000，有效转速 −20 000～45 532 r/min），最小计量单元：1 r/min，“0xFF，0xFE”表示默认值，“0xFF，0xFF”表示无效
档位	WORD	0 表示空档，1～6 分别表示 1～6 档，254 表示自动档，255 表示倒退档
电机负荷百分比	BYTE	0～100%，最小计量单元 1%，0xFE 表示默认值，0xFF 表示无效
发动机转速	WORD	0～65 532 r/min（数值偏移量 20 000，有效转速 −20 000～45 532 r/min），最小计量单元：1 r/min，“0xFF，0xFE” 表示异常，“0xFF，0xFF” 表示无效
发动机扭矩百分比	BYTE	0～100%，最小计量单元 1%， 0xFE 表示默认值，0xFF 表示无效
发动机负荷百分比	BYTE	0～100%，最小计量单元 1%， 0xFE 表示默认值，0xFF 表示无效效
电池电流	WORD	最小计量单元：0.05 A，偏移 −1 600 A，有效范围：−80～1 000 A
电池电压	WORD	最小计量单元：0.05 A，“0xFF，0xFE” 表示默认值，“0xFF，0xFF” 表示无效
空调 on/off	BYTE	1 开，0 关，0xFE 表示默认值，0xFF 表示无效
环境温度	BYTE	0～250 ℃，最小计量单元 1 ℃（数值偏移量 50 表示 −50～200 ℃），0xFE 表示默认值，0xFF 表示无效
环境湿度	BYTE	0～100%，最小计量单元 1%，0xFE 表示默认值，0xFF 表示无效
大气压力	DWORD	单位 Pa
油门踏板开度	BYTE	有效范围：0～200 （表示 0～100%），最小计量单元 0.5%，“0xFE” 表示异常，“0xFF” 表示无效
空燃比	WORD	0～100，最小计量单元 0.01，“0xFF，0xFE” 表示默认值，“0xFF，0xFF” 表示无效
瞬时油耗	WORD	0～65 532 L，最小计量单元 0.01 L，“0xFF，0xFE” 表示默认值，“0xFF，0xFF” 表示无效

下载: 导出CSV

表 3 参数筛选前后的对比

Table 3. The parameters are compared before and after the screening

时间	筛选前怠速比例/%	筛选后怠速比例/%	筛选前平均加速度/ （m•s⁻²）	筛选后平均加速度/ （m•s⁻²）	筛选前平均速度/ （km•h ⁻¹）	筛选后平均速度/ （km•h ⁻¹）	筛选前平均减速度/ （m•s⁻²）	筛选后平均减速度/ （m•s⁻²）
2016 年 11 月	23.30	24.40	0.50	0.45	26.3	27.5	−0.55	−0.54
2016 年 12 月	24.50	25.20	0.40	0.36	24.6	23.6	−0.50	−0.53
2017 年 01 月	27.30	28.30	0.45	0.46	23.3	23.5	−0.56	−0.55
2017 年 02 月	25.60	24.50	0.48	0.47	25.5	24.3	−0.52	−0.53
2017 年 03 月	26.40	27.10	0.47	0.46	24.4	24.6	−0.51	−0.51
2017 年 04 月	24.10	23.80	0.42	0.45	25.9	24.5	−0.52	−0.53
2017 年 05 月	24.80	24.50	0.46	0.43	25.5	24.8	−0.54	−0.53
2017 年 06 月	25.00	24.90	0.49	0.45	25.3	25.0	−0.52	−0.52
2017 年 07 月	23.80	23.20	0.48	0.44	26.3	24.8	−0.53	−0.54
2017 年 08 月	23.20	23.50	0.45	0.46	25.4	24.8	−0.51	−0.53
2017 年 09 月	24.70	24.80	0.46	0.45	24.8	24.9	−0.52	−0.53
2017 年 10 月	23.40	23.50	0.48	0.46	26.5	26.3	−0.53	−0.52

下载: 导出CSV

表 4 各主成分贡献率及累计贡献率

Table 4. Contribution rate and cumulative contribution rate of each principal component %

主成分	贡献率	累计贡献率
F₁	45.900 000	45.900 000
F₂	20.500 000	66.400 000
F₃	13.800 000	80.200 000
F₄	8.600 000	88.800 000
F₅	5.300 000	94.100 000
F₆	2.800 000	96.900 000
F₇	1.400 000	98.300 000
F₈	0.800 000	99.100 000
F₉	0.600 000	99.700 000
F₁₀	0.200 000	99.900 000
F₁₁	0.090 000	99.980 000
F₁₂	0.006 000	99.980 000
F₁₃	0.000 900	99.989 000
F₁₄	0.000 010	99.989 100
F₁₅	0.000 008	99.989 170
F₁₆	0.000 000	100.000 000

下载: 导出CSV

表 5 聚类分析迭代表

Table 5. Cluster analysis

聚类类别	迭代次数
聚类类别	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15
类1	0.369	0.298	0.220	0.126	0.071	0.072	0.044	0.032	0.021	0.010	0
类2	4.235	0.296	0.196	0.135	0.056	0.055	0.033	0.024	0
类3	5.768	0.996	0.632	0.438	0.239	0.152	0.073	0.062	0.057	0.032	0.021	0.015	0.005	0.003	0

下载: 导出CSV

表 6 实际采集数据与典型法规工况比较

Table 6. Comparison of actual data and typical statutory conditions

工况	平均速度/（km•h⁻¹）	运行平均速度/（km•h⁻¹）	匀速比例%	怠速比例%	加速比例%	减速比例%
实际行驶工况	24.50	32.60	22.90	24.84	26.60	25.66
CATC	25.87	33.92	20.59	23.72	28.56	27.13
FTP75	33.89	40.93	24.65	17.20	31.10	27.05
WLTP	46.42	53.15	27.83	12.67	30.94	28.56
NEDC	33.34	43.48	37.54	22.63	23.22	16.61

下载: 导出CSV

表 7 不同城市轻型汽车行驶工况对比

Table 7. Comparison of driving conditions in different cities of China

参数	北京	上海	广州	深圳	呼和浩特
匀速比例/%	27.34	22.28	25.95	18.40	20.59
怠速比例/%	16.52	31.61	17.77	20.15	23.72
加速比例/%	25.29	22.83	29.11	32.30	28.56
减速比例/%	30.85	23.28	27.16	29.15	27.13
平均速度/（km•h⁻¹）	19.98	14.96	14.14	32.39	25.87

下载: 导出CSV

参考文献(23)

孙强,白书战,韩尔樑,等. 基于试验测量的瞬时行驶工况构建[J]. 吉林大学学报(工学版),2015,45(2): 364-370.

SUN Qiang, BAI Shuzhan, HAN Erliang, et al. Instantaneous driving cycle construction based no experimental measurement[J]. Journal of Jilin University (Engineering and Technology Edition), 2015, 45(2): 364-370.

杨正军. WLTP用于中国道路轻型车测试循环可行性研究[D]. 北京: 清华大学, 2014.

BIELACZYC P, WOODBURN J, SZCZOTKA A. A comparison of carbon dioxide exhaust emissions and fuel consumption for vehicles tested over the NEDC, FTP-75 and WLTC chassis dynamometer test cycles[J]. SAE Technical Papers, 2015(1): 1065-1076.

GALGAMUWA U, PERERA L, BANDARA S. Developing a general methodology for driving cycle construction: comparison of various established driving cycles in the world to propose a general approach[J]. Journal of Transportation Technologies, 2015, 5(4): 191-203. doi: 10.4236/jtts.2015.54018

MINGYUE M, BENEDIKT W, FERIT K, et al. A statistical method for driving cycle construction based on path geometry[C]//Proceedings of 2013 the International Conference on Remote Sensing, Environment and Transportation Engineering. Nanjing: Atlantis Press, 2013: 890-893.

GAO X, ZHANG B, XIONG X, et al. Construction and analysis of the dalian driving cycle[J]. International Journal of Control & Automation, 2015, 8(6): 363-368.

KNEZ M, MUNEER T, JEREB B, et al. The estimation of a driving cycle for Celje and a comparison to other European cities[J]. Sustainable Cities & Society, 2014, 11(2/3): 56-60.

DIMARATOS A, TSOKOLIS D, FONTARAS G, et al. Comparative evaluation of the effect of various technologies on light-duty vehicle CO2, emissions over NEDC and WLTP[J]. Transportation Research Procedia, 2016, 14: 3169-3178. doi: 10.1016/j.trpro.2016.05.257

PREVEDOUROS P, MITROPOULOS L. Life cycle emissions and cost study of light duty vehicles[J]. Transportation Research Procedia, 2016, 15: 749-760. doi: 10.1016/j.trpro.2016.06.062

DEMUYNCK J, BOSTEELS D, PAEPE M D, et al. Recommendations for the new WLTP cycle based on an analysis of vehicle emission measurements on NEDC and CADC[J]. Energy Policy, 2012, 49(1): 234-242.

张宏,姚延钢,张露. 呼和浩特市轻型汽车运行工况构建研究[J]. 公路交通科技,2019,36(11): 104-110.

ZHANG Hong, YAO Yangang, ZHANG Lu. Study on construction of light-duty vehicle driving cycle in Hohhot city[J]. Journal of Highway and Transportation Research and Development, 2019, 36(11): 104-110.

沈姝,田冬莲,杜常清,等. 国Ⅴ在用车实际道路下IUPR分布特征研究[J]. 汽车技术,2016(6): 45-48. doi: 10.3969/j.issn.1000-3703.2016.06.011

SHEN Shu, TIAN Donglian, DU Changqing, et al. The research on IUPR distribution characteristic of China V vehicles based on the real road conditions[J]. Automobile Technology, 2016(6): 45-48. doi: 10.3969/j.issn.1000-3703.2016.06.011

路尧. 中国轻型汽车行驶工况开发[D]. 北京: 北京理工大学, 2016.

石琴,郑与波. 因子分析在汽车行驶工况构建过程中的应用[J]. 中国机械工程,2011(9): 1123-1127.

SHI Qin, ZHENG Yubo. Application of factor analysis in driving cycle development[J]. China Mechanical Engineering, 2011(9): 1123-1127.

王中奇,关伟,毕军,等. 基于短时实际行驶工况的机动车微观排放模型[J]. 环境工程学报,2016(10): 5803-5807. doi: 10.12030/j.cjee.201505184

WANG Zhongqi, GUAN Wei, BI Jun, et al. Microscopic emission model of motor vehicle based on short-time real driving cycle[J]. Chinese Journal of Environmental Engineering, 2016(10): 5803-5807. doi: 10.12030/j.cjee.201505184

杜常清,汪雨航,张佩. 基于GPS/GIS数据的短行程工况构建方法研究[J]. 武汉理工大学学报(交通科学与工程版),2016,40(5): 803-808.

DU Changqing, WANG Yuhang, ZHANG Pei. Research on short trip driving cycle development based on GPS/GIS data[J]. Journal of Wuhan University of Technology (Transportation Science & Engineering), 2016, 40(5): 803-808.

秦大同,詹森,漆正刚,等. 基于K-均值聚类算法的行驶工况构建方法[J]. 吉林大学学报(工学版),2016,46(2): 383-389.

QIN Datong, ZHAN Sen, QI Zhenggang, et al. Driving cycle construction using K-means clustering method[J]. Journal of Jilin University (Engineering and Technology Edition), 2016, 46(2): 383-389.

石琴,仇多洋,周洁瑜. 基于组合聚类法的行驶工况构建与精度分析[J]. 汽车工程,2012,34(4): 164-169.

SHI Qin, CHOU Duoyang, ZHOU Jieyu. Driving cycle construction and accuracy analysis based on combined clustering technique[J]. Automotive Engineering, 2012, 34(4): 164-169.

胡志远,秦艳,谭丕强,等. 基于大样本的上海市乘用车行驶工况构建[J]. 同济大学学报(自然科学版),2015,43(10): 1523-1527. doi: 10.11908/j.issn.0253-374x.2015.10.011

HU Zhiyuan, QIN Yan, TAN Piqiang, et al. Large sample based car driving cycle in Shanghai city[J]. Journal of Tongji University (Natural Science), 2015, 43(10): 1523-1527. doi: 10.11908/j.issn.0253-374x.2015.10.011

潘登. 混合动力汽车城市循环工况构建及运行工况多尺度预测[D]. 北京: 北京理工大学, 2015.

钟森鸣,黄坚,吴嘉健,等. 新能源汽车动态工况EMI测试系统构架设计与关键技术分析[J]. 中国测试,2017,43(8): 76-79.

ZHONG Senming, HUANG Jian, WU Jiajian, et al. Frame design and key technical analysis of EMI test system for new energy vehicle dynamic condition[J]. China Measurement & Test, 2017, 43(8): 76-79.

徐小俊,李君,刘宇,等. 电动汽车城市行驶工况构建[J]. 科学技术与工程,2017,17(35): 330-336. doi: 10.3969/j.issn.1671-1815.2017.35.055

XU Xiaojun, LI Jun, LIU Yu, et al. The development of electric vehicles urban driving cycle[J]. Science Technology and Engineering, 2017, 17(35): 330-336. doi: 10.3969/j.issn.1671-1815.2017.35.055

姚延钢. 城市道路轻型汽车运行工况构建[D]. 呼和浩特: 内蒙古大学, 2017.

施引文献

附加材料(0)

访问统计

点击查看大图

图(5) / 表(7)

计量

文章访问数: 681
HTML全文浏览量: 339
PDF下载量: 34
被引次数: 0

1. 超前导洞应用存在的问题分析及大变形控制措施
2. 理论分析模型
2.1 超前导洞及正洞围岩力学模型
2.2 统一强度准则与围岩劣化模型
2.3 塑性围岩扩容规律
3. 围岩弹塑性分析
3.1 基本方程
3.2 超前导洞开挖后围岩弹塑性分析
3.3 扩挖后围岩弹塑性分析
3.4 永久初期支护施作后围岩弹塑性分析
4. 算例验证及数值模拟
4.1 算例验证
4.2 数值模拟分析
5. 影响因素分析
5.1 软化模量（Mc、Mφ、ME）的影响
5.2 注浆参数的影响
5.3 超前导洞半径及两掌子面间距的影响分析
6. 结　论

1. 超前导洞应用存在的问题分析及大变形控制措施
2. 理论分析模型
2.1 超前导洞及正洞围岩力学模型
2.2 统一强度准则与围岩劣化模型
2.3 塑性围岩扩容规律
3. 围岩弹塑性分析
3.1 基本方程
3.2 超前导洞开挖后围岩弹塑性分析
3.3 扩挖后围岩弹塑性分析
3.4 永久初期支护施作后围岩弹塑性分析
4. 算例验证及数值模拟
4.1 算例验证
4.2 数值模拟分析
5. 影响因素分析
5.1 软化模量（M_c、M_φ、M_E）的影响
5.2 注浆参数的影响
5.3 超前导洞半径及两掌子面间距的影响分析
6. 结　论

参考文献(23)

施引文献

附加材料(0)

访问统计

城市道路轻型汽车行驶工况构建

doi: 10.3969/j.issn.0258-2724.20180224

作者简介: 张宏（1981—），女，副教授，博士，研究方向为城市交通数据科学及应用，E-mail：zhang2001hong@126.com

计量

出版历程

Light-Duty Vehicles Driving Cycle Construction Based on Urban Roads

1. 超前导洞应用存在的问题分析及大变形控制措施

2. 理论分析模型

2.1 超前导洞及正洞围岩力学模型

2.2 统一强度准则与围岩劣化模型

2.3 塑性围岩扩容规律

3. 围岩弹塑性分析

3.1 基本方程

3.2 超前导洞开挖后围岩弹塑性分析

3.2.1 弹性区

3.2.2 塑性软化区

3.2.3 塑性残余区

3.3 扩挖后围岩弹塑性分析

3.4 永久初期支护施作后围岩弹塑性分析

3.4.1 锚固弹性区

3.4.2 锚固塑性软化区

3.4.3 锚固注浆区

4. 算例验证及数值模拟

4.1 算例验证

4.2 数值模拟分析

5. 影响因素分析

5.1 软化模量（Mc、Mφ、ME）的影响

5.2 注浆参数的影响

5.3 超前导洞半径及两掌子面间距的影响分析

6. 结 论

计量

出版历程

目录

1. 超前导洞应用存在的问题分析及大变形控制措施

2. 理论分析模型

2.1 超前导洞及正洞围岩力学模型

2.2 统一强度准则与围岩劣化模型

2.3 塑性围岩扩容规律

3. 围岩弹塑性分析

3.1 基本方程

3.2 超前导洞开挖后围岩弹塑性分析

3.3 扩挖后围岩弹塑性分析

3.4 永久初期支护施作后围岩弹塑性分析

4. 算例验证及数值模拟

4.1 算例验证

4.2 数值模拟分析

5. 影响因素分析

5.1 软化模量（Mc、Mφ、ME）的影响

5.2 注浆参数的影响

5.3 超前导洞半径及两掌子面间距的影响分析

6. 结 论

作者简介:
张宏（1981—），女，副教授，博士，研究方向为城市交通数据科学及应用，E-mail：zhang2001hong@126.com

5.1 软化模量（M_c、M_φ、M_E）的影响

6. 结　论

5.1 软化模量（M_c、M_φ、M_E）的影响

6. 结　论