考虑浆土应力耦合作用的劈裂注浆机理分析

秦鹏飞; 钟宏伟; 刘坚; 苏丹娜; 孙卓宇

doi:10.3969/j.issn.0258-2724.20220433

考虑浆土应力耦合作用的劈裂注浆机理分析

doi: 10.3969/j.issn.0258-2724.20220433

秦鹏飞^{1, 2, 3,},
钟宏伟³,
刘坚²,
苏丹娜³,
孙卓宇³

1.
兰州大学西部灾害与环境力学教育部重点实验室，甘肃兰州 730000
2.
河海大学岩土力学与堤坝工程教育部重点实验室，江苏南京 370000
3.
郑州铁路职业技术学院铁道工程学院，河南郑州 451000

基金项目: 国家自然科学基金青年基金（62003313）；中央高校基本科研业务费专项资金资助（lzujbky-2021-kb03）；河海大学岩土力学与堤坝工程教育部重点实验室开放基金（2022LP002）

详细信息

作者简介:
秦鹏飞（1985—），男，副教授，研究方向为岩土与地下工程，E-mail：929163723@qq.com

中图分类号: TU43.1
计量
- 文章访问数: 502
- HTML全文浏览量: 116
- PDF下载量: 80
- 被引次数: 20
出版历程
- 收稿日期: 2022-05-20
- 修回日期: 2022-09-02
- 网络出版日期: 2023-04-10
- 刊出日期: 2023-01-05

Analysis of Split Grouting Mechanism Considering Coupling Effect of Slurry and Soil Stress

1.
Key Laboratory of Mechanics on Disaster and Environment in Western China, The Ministry of Education, Lanzhou University, Lanzhou 730000, China
2.
Key Laboratory of Geotechnical Mechanics and Dam Engineering of the Ministry of Education, Hohai University, Nanjing 370000, China
3.
School of Railway Engineering, Zhengzhou University of Railway Vocational Technology, Zhengzhou 451000, China

摘要

摘要:
浆土应力耦合作用对劈裂注浆浆液扩散规律具有显著影响，砂土劈裂注浆设计应充分考虑这种影响作用. 将劈裂注浆视作平面无限域的圆形扩张过程，基于牛顿型本构方程分析了浆液流场变化特征，并将劈裂通道下侧砂层视作半无限空间弹性体，采用弹性力学推导了均布荷载下劈裂通道宽度、浆液压力的分布方程. 通过设置不同的浆液黏度、砂土弹性模量参数，深刻揭示了耦合效应下砂土劈裂注浆基本机理. 结合郑州地铁某在建工程进行了对比验证. 研究结果表明：浆液压力在孔口及远端处急速衰减，而在中间区段呈稳定变化趋势，劈裂通道宽度基本由浆液压力决定，其分布趋势与浆液压力分布趋势相同；浆液黏度、砂土弹性模量是影响劈裂扩散半径的重要因素，黏度和模量均与扩散半径正相关，黏度与劈裂宽度正相关，模量则与劈裂宽度负相关；本文理论计算值与现场开挖实际值偏差12%～15%，基本符合预期要求.
- 浆-土耦合效应 /
- 劈裂注浆 /
- 起劈压力 /
- 扩散机理 /
- 工程应用
Abstract:
The coupling effect of slurry and soil stress has a significant influence on the slurry diffusion law of split grouting, which should be fully considered in the design phase of split grouting. The split grouting is regarded as a circular expansion process in an infinite plane domain. Based on the Newtonian constitutive equation, variation characteristics of the slurry flow field are analyzed. The sand layer under the splitting channel is regarded as a semi-infinite elastic body. The distribution equations of the splitting channel width and the slurry pressure under the uniform load are derived by elastic mechanics. By setting different slurry viscosity and sand elastic modulus parameters, the basic mechanism of split grouting in sandy soil under coupling effect is revealed. The results show that the slurry pressure decreases rapidly at the orifice and the far end, but changes slowly in the middle section. The width of the splitting channel is basically determined by the slurry pressure, and has a distribution trend the same as that of the slurry pressure. The viscosity of slurry and the elastic modulus of sand are important factors affecting the splitting diffusion radius , and both are positively related to the diffusion radius. The slurry viscosity is positively related to the splitting width, while the modulus is negatively related to the splitting width. Finally, a comparison was made through a project under construction, and it was found that the deviation between the theoretical calculation value and the actual excavation value was 12%−15%, which basically met the expected requirements.
- slurry-soil coupling effect /
- split grouting /
- splitting pressure /
- diffusion mechanism /
- engineering application

HTML全文

固定辙叉是道岔重要组成部分，广泛应用于我国重载和提速铁路. 固定辙叉区轮轨动力相互作用异常剧烈，关键钢轨部件磨耗磨损等病害较区间线路更加频发^[1-6]. 我国固定辙叉平均使用寿命仅为美国的40%，欧洲的60%，造成其服役寿命相对较低的主要原因之一是轮载过渡区辙叉钢轨总承载面积相对欧美较小^[7].

辙叉区不规则的钢轨廓形是影响轮载过渡区钢轨承载面积的直接因素. 为此，国内外学者从钢轨廓形参数优化角度进行了大量研究：王树国等^[7]从理论分析、动力仿真和现场试验角度对辙叉心轨顶宽、查照间隔和护轨轮缘槽宽度等参数进行优化；徐井芒等^[8-9]提出基于轮轨廓形净差值的心轨廓形优化方法；曹洋等^[10]针对固定辙叉轮载过渡段轮轨匹配性能对心轨关键断面廓形进行优化；沈钢等^[11]采用粒子群算法求解优化模型对心轨顶面降低值参数进行优化；张鹏飞等^[12]研究翼轨加高值对列车过岔动力特性的影响；Wan等^[13-14]采用数值模拟和现场试验分析辙叉区各轮轨动力响应的权重大小，并结合多目标优化方法提出固定辙叉心轨廓形设计方法.

上述文献研究主要针对辙叉廓形单一几何参数，通过调整局部廓形进行优化设计，并未从固定辙叉几何特征控制参数角度探寻廓形参数化设计方法并分析各参数对辙叉钢轨廓形的影响. 因此，本文基于B样条理论，采用DOE （design of experiments）方法，分析各关键控制参数对心轨廓形变化的影响权重，对实现固定辙叉区钢轨廓形参数化设计具有重要意义.

1. 固定辙叉廓形特征参数

在固定辙叉中，有害空间引起轨距线不连续，且翼轨和心轨廓形呈连续不规则变化. 然而，固定辙叉平面布置、心轨各关键断面廓形以及翼轨断面几何特征具有如下规律：对于任意心轨和翼轨断面，轨顶横坡和轨头侧面纵坡不变；沿辙叉走行方向，心轨和翼轨高度分段线性变化，心轨断面宽度线性增大；心轨和翼轨轮缘槽宽度一定；距轨距测量点标准轨头半宽35.5 mm处的翼轨顶面高度以及翼轨宽度可精确确定. 以60 kg/m钢轨12号固定辙叉^[15]为例，心轨轨顶横坡为1∶20，轨头侧面纵坡为1∶5；翼轨工作边纵坡为1∶5，轨顶横坡为1∶20，非工作边侧面纵坡1∶20，如图1所示，W为翼轨轨头宽度.

图 1 固定辙叉几何特征

Figure 1. Geometry characteristics of fixed-nose crossing

下载: 全尺寸图片幻灯片

由此可见，虽然辙叉钢轨轨头宽度和高度沿辙叉方向不断变化，但仍存在保持不变的设计参数，分别为各断面处心轨和翼轨轨头宽度和高度、轨顶横坡、轨头侧面坡度. 本文将以上参数统称为固定辙叉钢轨全断面廓形特征参数，为后文中不断重构更新辙叉钢轨廓形提供了设计依据.

2. 辙叉钢轨廓形参数化拟合方法

2.1 B样条曲线理论

辙叉区钢轨各断面廓形是由多段不同曲率圆弧组成的复杂曲线，整个廓形需保持连续平滑. B样条曲线在计算机辅助设计中具有广泛应用^[16]，可通过较少的轮廓关键控制顶点实现达到期望逼近精度的廓形拟合，且具有连续平滑性和强凸包性. 因此，通过B样条曲线拟合辙叉钢轨廓形，对利用辙叉廓形特征参数进行廓形参数化设计具有明显优势. 已有相关研究将其应用于钢轨廓形设计：Yang等^[17]采用B样条理论，选取实测廓形33个控制点建立磨耗钢轨廓形曲线，对不同打磨方式钢轨廓形进行磨耗预测；温士明^[18]运用B样条曲线对钢轨廓形进行参数化描述，提出钢轨打磨目标廓形优化方法；刘建桥^[19]结合遗传算法和B样条曲线，以降低钢轨磨耗为目标确定钢轨优化廓形；唐彦玲等^[20]通过B样条曲线对钢轨廓形进行参数化描述，结合遗传算法对曲线地段钢轨廓形进行设计；王亮等^[21]采用B样条曲线和NSGA-Ⅱ算法对控制点权因子寻优计算，对钢轨打磨区域轨头廓形进行设计；林凤涛等^[22]采用B样条理论对高速铁路辙叉钢轨廓形拟合，基于标准辙叉廓形设置16个控制点设计了钢轨经济性打磨廓形；史振帅^[23]以道岔转辙器区域直尖轨为研究对象，采用NURBS理论对道岔直尖轨打磨廓形进行设计重构. 已有文献虽实现了钢轨廓形重构，但仍是通过标准或实测钢轨廓形的稠密离散点数据进行廓形拟合，且大多针对区间线路和可动心辙叉，对固定辙叉钢轨廓形研究较少.

在B样条方法中，任意曲线可用式（1）表示.

$C(t) = \sum\limits_{v = 0}^{n{{ - }}1} {{B_{v,d}}({{t}})Q{}_v},$

(1)

式中：Q_v为序号v的廓形关键控制点；n为控制点Q_v个数，依次连接控制点可以构成一个特征多边形；t为节点，t₀≤t≤t_m-1，T=（t₀，t₁，…，t_v，…t_m-1），T为节点向量，m=n + d + 1，m为节点个数，d为B样条曲线的阶数；B_v,d(t)为由递归方程表示的d−1次B样条基函数，如式（2）所示.

$\left\{\begin{array}{l}{B}_{v,d}({{t}})=\dfrac{{{t}}-{t}_{v}}{t{}_{v + d}-{t}_{v}}{B}_{v,d-1}({{t}}) + \dfrac{t_{v + d + 1}-{{t}}}{t_{v + d + 1}-{t}_{v + 1}}{B}_{v + 1,d-1}({{t}}),\\ {B}_{v,0}=\left\{\begin{array}{l}1\text{，}\quad{{{t}}}_{v}\leqslant {{t}}\leqslant {t}_{v + 1},\\ 0\text{，}\quad {{t}} < {t}_{v}\;或\;{{t}} > {t}_{v + 1}. \end{array}\right. \end{array} \right.$

(2)

根据B样条曲线理论与相关工程应用情况可知，采用3次样条可保证拟合曲线2阶导数连续，满足工程实际需要，因此，本文采用3次非均匀B样条曲线拟合辙叉钢轨廓形.

2.2 关键控制点确定

依据心轨和翼轨廓形关键几何特征参数可确定B样条曲线拟合廓形所需关键控制点，由各控制点组成心轨和翼轨廓形拟合控制多边形，如图2所示. 由于心轨断面廓形轴对称，因此只选取心轨中心轴一侧进行廓形参数化说明.

图 2 辙叉钢轨廓形控制多边形

Figure 2. Fixed-crossing rail profile control polygon

下载: 全尺寸图片幻灯片

在图2（a）中，心轨廓形拟合的关键控制点为{N₁,N₂…,N_i}，N_i=(x_i，y_i)，i =1，2，…，8. 其中：N₁为心轨轨顶高度控制点；N₂、N₃用来调节心轨顶面平直段廓形范围；N₄、N₆、N₈用来调节心轨轨头复合圆弧段廓形；N₅为轨顶横坡和侧面纵坡交点；N₇为心轨宽度控制点；P为△N₄N₅N₆中垂线与N₄N₆的交点，为确定控制点N₈的辅助点，N₈位于N₅P的连线上.

在图2（b）中，翼轨廓形拟合的关键控制点为{W₁,W₂,…,W_j}，j=1，2，…，9. 其中：W₁和W₉为翼轨轨头宽度控制点；W₃为轨顶横坡和翼轨工作边纵坡交点；W₅为翼轨高度控制点；W₇为轨顶横坡和翼轨非工作边纵坡交点；W₂、W₄、W₆、W₈用来调节翼轨轨头两侧复合圆弧段廓形；k_W1、k_W2、k_W3分别为翼轨工作边纵坡、轨顶横坡和非工作边纵坡斜率；h_W为翼轨高度；w_W为标准轨头宽度的一半，即35.5 mm.

由于心轨和翼轨廓形参数化方法一致，且心轨廓形变化较翼轨复杂，因此，以心轨为例详细说明辙叉钢轨廓形参数化过程. 采用向量法确定各控制点间关系，如式（3）.

$\left\{\begin{array}{l} \dfrac{{{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}}{\left|{{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}\right|}=(1\text{，}{k}_{1}),\\ \dfrac{{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}}{\left|{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}\right|}=(1\text{，}{k}_{2}),\\ {{\boldsymbol{N}}}_{4}{{\boldsymbol{N}}}_{5}={\lambda }_{1}{{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{5},\\ {{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{2}={\lambda }_{2}{{\boldsymbol{N}}}_{1}{{\boldsymbol{N}}}_{4},\\ {{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{6}={\lambda }_{3}{{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{7},\\ {{\boldsymbol{N}}}_{2}{{\boldsymbol{N}}}_{3}={\lambda }_{4}{{\boldsymbol{N}}}_{2}{{\boldsymbol{N}}}_{4},\\ {{\boldsymbol{N}}}_{5}{{\boldsymbol{N}}}_{8}={\lambda }_{5}{{\boldsymbol{N}}}_{5}{\boldsymbol{P}}, \end{array}\right.$

(3)

式中：λ₁、λ₂、λ₃、λ₄、λ₅分别为调节控制点N₄、N₂、N₆、N₃、N₈位置的比例系数，为可变参数，记为λ_s∈[0，1]，其中，s=1,2,…，5，λ₅用来调节心轨轨头复合圆弧段廓形，简称心轨复合圆弧段比例系数；k₁、k₂分别为心轨轨顶横坡和侧面纵坡斜率.

当其他控制点一定时，心轨轨头复合圆弧段廓形主要由控制点N₈位置决定，为了保持心轨轨头廓形凸曲线的特点，N₈取值需在N₅P范围内. 依据向量法求三角形垂点P，如式（4）所示.

${\boldsymbol{P}} = {{\boldsymbol{N}}_4} + \frac{{{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_5} {\text{•}} {{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}}}{{\left| {{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}} \right| {\text{•}} \left| {{{\boldsymbol{N}}_4}{{\boldsymbol{N}}_6}} \right|}} {\text{•}} ({{\boldsymbol{N}}_6} - {{\boldsymbol{N}}_4}).$

(4)

结合式（3）、（4），根据定比分点原理确定各关键控制点坐标，如式（5）所示.

$\left\{ \begin{array}{l}\left({x}_{1},{y}_{1}\right)=(0,h),\\ \left({x}_{2},{y}_{2}\right)=(1-{\lambda }_{2})({x}_{1},{y}_{1}) + {\lambda }_{2}({x}_{4},{y}_{4}) + (0\text{，}\varDelta ),\\ \left({x}_{3},{y}_{3}\right)=(1-{\lambda }_{4})({x}_{2},{y}_{2}) + {\lambda }_{4}({x}_{4},{y}_{4}),\\ \left({x}_{4},{y}_{4}\right)={\lambda }_{1}({x}_{1},{y}_{1}) + (1-{\lambda }_{1})({x}_{5},{y}_{5}),\\ \left({x}_{5},{y}_{5}\right)=\left(\dfrac{h + {k}_{2}w}{{k}_{2}-{k}_{1}},h + \dfrac{{k}_{1}(h + {k}_{2}w)}{{k}_{2}-{k}_{1}}\right),\\ \left({x}_{6},{y}_{6}\right)=(1-{\lambda }_{3})({x}_{5},{y}_{5}) + {\lambda }_{3}({x}_{7},{y}_{7}),\\ \left({x}_{7},{y}_{7}\right)=(w,0),\\ \left({x}_{8},{y}_{8}\right)=(1-{\lambda }_{5})({x}_{5},{y}_{5}) + {\lambda }_{5}({x}_{\text{p}},{y}_{\text{p}}), \end{array} \right.$

(5)

式中： (x_p，y_p)为点P坐标；Δ为心轨平直段廓形修正系数，通过调整Δ使控制点N₂与N₁在同一水平线上；h为x处心轨高度，x为心轨任意断面距辙叉理论尖端的距离；w为x处心轨轨头半宽.

由式（3）～（5）可知，心轨任意断面廓形各控制点坐标可表示为各特征参数和比例系数的函数. 翼轨同理，则固定辙叉全断面廓形控制点可表示为

$\left\{ \begin{gathered} {{{\boldsymbol{N}}_i}} = {f_{{\text{cp}}}}(h,w,{k_1},{k_2},\varDelta ,[{\lambda _s}]), \\ {{{\boldsymbol{W}}_j}} = {g_{{\text{wp}}}}({h_{\text{W}}},{w_{\text{W}}},{k_{\text{W}}}_1,{k_{\text{W}}}_2,[{\lambda _t}]), \\ \end{gathered} \right.$

(6)

式中：f_cp(·)、g_wp(·)分别为求解心轨和翼轨任意断面所有控制点的函数；λ_t为翼轨控制点调节比例系数，记为λ_t∈[0，1]，其中，t为其对应的比例系数编号.

综上所述，关键控制点确定所涉及的参数具体包括两类，一类是与廓形相关的有实际物理意义的特征参数，另一类是B样条曲线几何控制点调节比例系数. 为方便后续分析，将上述两类参数统称为辙叉钢轨廓形关键控制参数.

2.3 廓形拟合

以60 kg/m钢轨12号固定辙叉钢轨为例，采用B样条曲线进行廓形构建，心轨和翼轨控制点分布及拟合廓形分别如图3所示.

图 3 辙叉钢轨断面廓形拟合

Figure 3. Section fitting of rail profiles in fixed-nose crossing

下载: 全尺寸图片幻灯片

由图3可知：各控制点位置随心轨轨顶高度和宽度的变化而改变，从而实现心轨廓形连续变化；本文不考虑翼轨高度和轨头宽度的变化，因此翼轨拟合廓形唯一.

在此基础上，结合辙叉区心轨和翼轨走行方向变化，可进一步确定辙叉三维廓形，考虑一侧翼轨的固定辙叉三维廓形如图4所示.

图 4 固定辙叉三维廓形实现

Figure 4. Three-dimensional profile of fixed-nose crossing

下载: 全尺寸图片幻灯片

由图4可知，虽然固定辙叉钢轨廓形复杂多变，心轨和翼轨沿辙叉走行方向廓形不断变化，但结合辙叉结构设计中各关键断面的关键控制参数，基于B样条方法可快速确定沿辙叉走行方向的任意截面廓形，该方法可极大拓展以固定辙叉廓形优化设计为基础的轮轨动力仿真以及轮轨接触关系等的研究.

3. 心轨廓形参数试验设计分析方法（DOE）分析

3.1 评价指标

心轨廓形关键控制参数调整会改变心轨廓形，廓形面积差可以直观量化地表达心轨整体廓形的变化情况；与此同时，不同参数对心轨轨头廓形和心轨侧面廓形变化影响程度不同，而心轨轨头廓形变化量和侧面廓形变化量可以更细致地反映各参数对心轨轨头和侧面廓形的单独影响. 因此，为分析各关键控制参数对心轨廓形的影响，选取ΔS以及D₁、D₂评价各方案廓形与标准廓形的贴合度，如式（7）所示. 本文借鉴《铁路线路修理规则》中对钢轨垂磨和侧磨的定义，选取心轨轨头顶面和侧面某一位置处廓形距离偏差量来定义D₁和D₂.

$\left\{ \begin{gathered} \Delta S = \left| {{S_{{\text{case}}}} - {S_{{\text{std}}}}} \right|, \\ {D_1} = \left| {{Y_{{\text{case}}}} - {Y_{{\text{std}}}}} \right|, \\ {D_2} = \left| {{X_{{\text{case}}}} - {X_{{\text{std}}}}} \right|, \\ \end{gathered} \right.$

(7)

式中：S_case为各方案廓形曲线与X坐标轴围成的面积；S_std为采取表1参数确定的标准廓形与X坐标轴围成的面积；ΔS为各方案廓形与标准廓形面积差，反映心轨整体廓形的变化程度；Y_case和Y_std分别为距心轨中心轴2w/3处，各方案廓形曲线与标准廓形曲线纵坐标值；X_case和X_std分别为距心轨顶面h/2处各方案廓形曲线与标准廓形曲线横坐标值；D₁和D₂分别为各方案廓形与标准廓形在垂向和横向上的廓形偏差.

表 1 DOE全阶乘方案设计

Table 1. DOE full factorial scheme design

方案	k₁	k₂	λ₅
方案 1	1/30	1/6	1/4
方案 2	1/30	1/6	1/2
方案 3	1/30	1/4	1/4
方案 4	1/30	1/4	1/2
方案 5	1/10	1/6	1/4
方案 6	1/10	1/6	1/2
方案 7	1/10	1/4	1/4
方案 8	1/10	1/4	1/2
标准	1/20	1/5	1/3

下载: 导出CSV

| 显示表格

在DOE分析中，为明确各关键控制参数在所选合理取值范围内对心轨廓形变化的影响权重E，除了对D₁、D₂及ΔS各自的均值进行分析，还需对E进行分析，如式（8）所示.

$E = \frac{{\dfrac{{I^+ - I^ - }}{{V_b^ + - V_b^{{ - }}}}}}{{\displaystyle\sum\limits_{b = 1}^3 {\dfrac{{I^ + - I^ - }}{{V_b^ + - V_b^{{ - }}}}} }} \times 100{\text{%}},$

(8)

式中： $I^ -$ 和 $I^ +$ 分别为各控制参数分别取2个水平（−和 + ）时的各贴合度指标均值， $V_b^ +$ 和 $V_b^ -$ 分别为各关键参数的2个水平值（−和 + ）.

3.2 方案设置

在工程实践中，心轨轨头顶面和轨头侧面复合圆弧部分更易发生严重磨耗^[24]，k₁和k₂分别对心轨顶面廓形和轨头侧面廓形起控制作用，比例系数λ₅能够调整控制点N₈的位置，从而实现轨头侧面复合圆弧段廓形的调整. 因此，采用DOE全阶乘设计方法对心轨廓形关键控制参数k₁、k₂及λ₅进行方案设计. 依据《道岔设计手册》60 kg/m钢轨12号道岔^[15]中的心轨标准设计参数，在一定范围内对各参数调整取值进行方案设计，以满足心轨整体廓形符合工程实际. 各方案参数取值如表1.

3.3 重构廓形轮轨匹配特征

轮轨接触点分布规律是轮轨廓形匹配的最直观反映^[25]. 选取LMA型车轮踏面，以心轨20 mm断面和50 mm断面为例，分析轮对不同横移量下各方案轮轨接触点分布规律，进一步验证该参数化设计方法的合理性. 各方案心轨拟合廓形如图5所示，各方案轮轨接触点分布情况如图6所示.

图 5 心轨拟合廓形

Figure 5. Nose rail fitting profiles

下载: 全尺寸图片幻灯片

图 6 心轨关键断面车轮踏面接触点变化规律

Figure 6. Contact point variations on wheel tread for key nose rail sections

下载: 全尺寸图片幻灯片

由图5可知：各方案心轨廓形与标准廓形贴合度良好，说明取得了较好的钢轨廓形重构效果. 所建立的钢轨廓形B样条曲线构造方法为后文参数分析中不断重构更新钢轨廓形奠定了基础.

由图6可知：各方案拟合廓形轮轨接触点分布规律均与现有标准廓形一致；在心轨顶宽20 mm断面，各方案轮轨接触点位置基本一致，仅接触点发生突变对应的横移量有极小差别，说明在心轨顶宽较小时，各参数对轮轨接触点分布影响较小；在心轨顶宽50 mm断面，与标准廓形方案相比，方案1～4中轮轨接触点分布偏向于车轮踏面外侧，变化规律差别较小；而方案5～8则偏向于车轮踏面轮缘侧，变化规律也基本相近；随着心轨顶宽增大，参数k₁对轮轨接触点分布的影响较k₂和λ₅大.

3.4 DOE分析

固定辙叉钢轨廓形影响轮轨匹配性能，因此，明确各控制参数对辙叉钢轨廓形影响权重，对辙叉钢轨廓形参数化设计中的参数调整优化具有重要意义. 以心轨20 mm断面和50 mm断面为例，通过DOE设计方法分析各控制参数对心轨廓形变化的影响效果及影响权重. 图7、8所示为心轨各控制参数对D₁、D₂和ΔS变化范围的影响. 图9为反映各控制参数对心轨廓形绝对影响效应的Pareto图. 图10为各控制参数对心轨廓形变化的影响权重.

图 7 心轨20 mm断面

Figure 7. Nose rail with section of 20 mm

下载: 全尺寸图片幻灯片

图 8 心轨50 mm断面

Figure 8. Nose rail with section of 50 mm

下载: 全尺寸图片幻灯片

图 9 各指标影响效应的Pareto图

Figure 9. Pareto chart of influencing effect of all indexes

下载: 全尺寸图片幻灯片

图 10 各控制参数影响权重

Figure 10. Influencing weights of control parameters

下载: 全尺寸图片幻灯片

由图7、8可知：在心轨20 mm断面，k₂和λ₅对ΔS、D₁和D₂变化范围均影响较大，k₁对三者影响较小；在心轨50 mm断面，k₁对ΔS、D₁影响较大，对D₂影响较小，k₂和λ₅对ΔS、D₁和D₂影响相当. 说明当心轨轨头宽度较小时，心轨整体廓形变化、心轨顶面和侧面廓形均主要受k₂和λ₅影响，k₁影响较小；随着心轨轨头宽度增大，k₁对心轨整体廓形和顶面廓形影响增大，对侧面廓形影响较小；k₂和λ₅对心轨各部位廓形影响效果相当.

由图9可知：在心轨20 mm断面，k₁和λ₅对ΔS影响效应较大，各参数对D₁和D₂影响效应相当且均较小；在心轨50 mm断面，各参数对ΔS影响效应最大，D₁次之，对D₂影响很小，其中，k₁对ΔS影响效应最大，k₂和λ₅对ΔS影响效应相当；各参数对D₁影响效应相当.

由图10可知：1）在心轨20 mm断面，k₁和λ₅对ΔS和D₁影响权重较小，而k₂影响较大. 其中，ΔS对应的E分别为21.08%、56.89%和22.02%，D₁对应的E分别为8.42%、61.95%、29.63%. 说明在心轨20 mm断面，k₁和λ₅对心轨整体廓形和顶面廓形的影响权重远不及k₂，主要原因是在心轨轨头宽度较小时，侧面高度相对较大，心轨侧面廓形对心轨廓形变化起主导作用. 2）在心轨50 mm断面，k₁、k₂和λ₅对ΔS和D₁的影响权重依次减小. 其中，ΔS对应的E分别为55.90%、33.38%、10.72%，D₁对应的E分别为64.81%、26.71%、8.48%，说明随着心轨轨头宽度增大，各关键控制参数在方案所选取值范围内时，k₁对心轨顶面和心轨整体廓形的影响权重逐渐明显，k₂和λ₅的影响逐渐减弱. 3）在心轨20 mm断面，k₁、k₂和λ₅对心轨侧面廓形偏差D₂对应的E分别为9.93%、76.82%、13.25%；在心轨50 mm断面，k₁、k₂和λ₅对心轨侧面廓形偏差D₂对应的E分别为14.09%、66.04%、19.87%；由此可见，由于不同心轨轨头宽度对应的高度变化不大，因此，无论心轨宽度大小，心轨侧面廓形主要受k₂影响较大，受λ₅影响次之，受k₁影响最小.

4. 结　论

1）基于B样条理论，提出了一种考虑固定辙叉全断面钢轨廓形几何特征的参数化设计方法，该方法利用B样条曲线定义工业产品几何形状的优势，通过固定辙叉钢轨全断面廓形关键特征参数确定任意断面廓形控制点，进而实现钢轨廓形参数化重构；提出了钢轨重构廓形贴合度评价指标和关键控制参数对廓形贴合度影响权重指标，可直观量化地反映廓形参数变化对心轨廓形的影响.

2）在心轨20 mm和50 mm断面，各方案心轨重构廓形与标准廓形贴合度良好，不同横移量下轮轨接触点分布规律符合实际. 各参数对轮轨接触点变化规律与心轨顶面宽度有关，在心轨顶宽较小时，各参数对轮轨接触点变化影响较小，在心轨顶宽较大时，k₁对轮轨接触点分布的影响较k₂和 λ₅大.

3）当心轨轨头宽度较小时，心轨整体廓形、心轨顶面和侧面廓形变化范围主要受k₂和λ₅影响，k₁影响较小；k₂对各部分廓形变化影响权重较大，k₁和λ₅影响较小；随着心轨轨头宽度的增大，k₁对心轨整体廓形和顶面廓形变化范围和影响权重均逐渐增大，k₂和λ₅的影响逐渐减弱. 由此可见，各控制参数对心轨各部分廓形影响权重不同，通过考虑各控制参数对钢轨各部分廓形的影响权重，有针对性地重点考虑相应的控制参数进行廓形方案设计，避免了经验—性能计算—修正等反复试凑的传统设计过程，使辙叉廓形设计更具有针对性和可预测性.

图 1 砂土劈裂注浆机理分析

Figure 1. Mechanism analysis of split grouting in sandy soil

下载: 全尺寸图片幻灯片

图 2 劈裂注浆力学机理分析

Figure 2. Mechanical mechanism analysis of split grouting

下载: 全尺寸图片幻灯片

图 3 注浆压力沿程分布

Figure 3. Grouting pressure distribution

下载: 全尺寸图片幻灯片

图 4 劈裂注浆浆-土耦合作用分析

Figure 4. Slurry-soil coupling effect analysis of split grouting

下载: 全尺寸图片幻灯片

图 5 浆液扩散半径与注浆压力差关系

Figure 5. Relationship between slurry diffusion radius and grouting pressure

下载: 全尺寸图片幻灯片

图 6 注浆压力沿程衰减规律

Figure 6. Relationship between slurry diffusion radius and grouting pressure

下载: 全尺寸图片幻灯片

图 7 劈裂宽度沿程衰减规律

Figure 7. Relationship between slurry diffusion radius and splitting gap

下载: 全尺寸图片幻灯片

图 8 隧道开挖揭露浆脉

Figure 8. Slurry veins exposed during tunnel excavation

下载: 全尺寸图片幻灯片

表 1 郑州地铁区间隧道砂土力学参数

Table 1. Mechanical parameters of sandy soil in Zhengzhou Metro Tunnel

参数天然密度/（g·cm⁻³）干密度/（g·cm⁻³）弹性模量/MPa 泊松比含水率/% 含泥量/% 注浆影响范围/m

取值 2.01 1.68 5.94 0.38 23.3 11.8 10

下载: 导出CSV

参考文献(28)

[1]	DONGSOON P, JEHEON O. Permeation grouting for remediation of dam cores[J]. Engineering Geology, 2018, 233: 63-75. doi: 10.1016/j.enggeo.2017.12.011
[2]	张庆松,王洪波,刘人太,等. 考虑浆液扩散路径的多孔介质渗透注浆机理研究[J]. 岩土工程学报,2018,40(5): 918-924. doi: 10.11779/CJGE201805017 ZHANG Qingsong, WANG Hongbo, LIU Rentai, et al. Infiltration grouting mechanism of porous media considering diffusion paths of grout[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(5): 918-924. doi: 10.11779/CJGE201805017
[3]	RASOULI R, HAYASHI K, ZEN K. Controlled permeation grouting method for mitigation of liquefaction[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(11): 1-11.
[4]	刘人太,郑卓,李术才,等. 地下工程裂隙型涌水超前注浆治理方法[J]. 西南交通大学学报,2019,54(1): 48-54. LIU Rentai, ZHENG Zhuo, LI Shucai, et al. Pre-excavation grouting to prevent water inflow of jointed rock in underground construction[J]. Journal of Southwest Jiaotong University, 2019, 54(1): 48-54.
[5]	李召峰,李术才,刘人太,等. 富水破碎岩体注浆加固实验与机制研究[J]. 岩石力学与工程学报,2017,36(1): 198-207. doi: 10.13722/j.cnki.jrme.2016.0119 LI Zhaofeng, LI Shucai, LIU Rentai, et al. Grouting reinforcement experiment for water-rich broken rock mass[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(1): 198-207. doi: 10.13722/j.cnki.jrme.2016.0119
[6]	张庆松,张连震,张霄,等. 基于浆液黏度时空变化的水平裂隙岩体注浆扩散机制[J]. 岩石力学与工程学报,2015,34(6): 1198-1210. doi: 10.13722/j.cnki.jrme.2014.0958 ZHANG Qingsong, ZHANG Lianzhen, ZHANG Xiao, et al. Grouting diffusion in a horizontal crack considering temporal and spatial variation of viscosity[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(6): 1198-1210. doi: 10.13722/j.cnki.jrme.2014.0958
[7]	叶飞,秦楠,梁兴,等. 基于驱替效应的盾构壁后注浆微观模型分析[J]. 西南交通大学学报,2022,57(2): 339-345. doi: 10.3969/j.issn.0258-2724.20200344 YE Fei, QIN Nan, LIANG Xing, et al. Microscopic model analysis of shield tunnel backfill grouting based on displacement effect[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 339-345. doi: 10.3969/j.issn.0258-2724.20200344
[8]	张连震,李志鹏,张庆松,等. 基于土体非线性压密效应的劈裂注浆机制分析[J]. 岩石力学与工程学报,2016,35(7): 1483-1493. doi: 10.13722/j.cnki.jrme.2015.0565 ZHANG Lianzhen, LI Zhipeng, ZHANG Qingsong, et al. Split grouting mechanism based on nonlinear characteristics of compression process of soil[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7): 1483-1493. doi: 10.13722/j.cnki.jrme.2015.0565
[9]	张庆松,张连震,刘人太,等. 基于“浆-土”界面应力耦合效应的劈裂注浆理论研究[J]. 岩土工程学报,2016,38(2): 323-330. ZHANG Qingsong, ZHANG Lianzhen, LIU Rentai, et al. Split grouting theory based on slurry-soil coupling effects[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(2): 323-330.
[10]	秦鹏飞. 不良地质体注浆细观力学模拟研究[J]. 煤炭学报,2020,45(7): 2646-2654. doi: 10.13225/j.cnki.jccs.2019.0634 QIN Pengfei. Study on meso-mechanical simulation of grouting in bad geo-body[J]. Journal of China Coal Society, 2020, 45(7): 2646-2654. doi: 10.13225/j.cnki.jccs.2019.0634
[11]	申时钊,涂小兵,雷进生,等. 不同渗透系数的非均质黏土劈裂注浆数值模拟[J]. 水利水运工程学报,2022(5): 102-112. SHEN Shizhao, TU Xiaobing, LEI Jinsheng, et al. Numerical simulation of split grouting of heterogeneous clay with different permeability coefficient[J]. Hydro-Science and Engineering, 2022(5): 102-112.
[12]	刘海明,卢昊正,南敢,等. 基于黏度时变性的Herschel-Bulkley流体劈裂注浆扩散特性研究[J]. 自然灾害学报,2022,31(3): 213-221. LIU Haiming, LU Haozheng, NAN Gan, et al. Research on diffusion characteristics of Herschel-Bulkley fluid fracture grouting based on time variation of viscosity[J]. Journal of Natural Disasters, 2022, 31(3): 213-221.
[13]	周茗如,张建斌,卢国文,等. 扩孔理论在非饱和黄土劈裂注浆中的应用[J]. 建筑结构学报,2018,39(增1): 368-378. doi: 10.14006/j.jzjgxb.2018.S1.048 ZHOU Mingru, ZHANG Jianbin, LU Guowen, et al. Application of expansion hole theory in fracture grouting of unsaturated loess[J]. Journal of Building Structures, 2018, 39(S1): 368-378. doi: 10.14006/j.jzjgxb.2018.S1.048
[14]	蔡德国,叶飞,曹凯,等. 砂性地层盾构隧道壁后注浆浆液扩散室内试验[J]. 中国公路学报,2018,31(10): 274-283. CAI Deguo, YE Fei, CAO Kai, et al. Test of grout diffusion of shield tunnel backfill grouting in sandy strata[J]. China Journal of Highway and Transport, 2018, 31(10): 274-283.
[15]	WANG S Y, CHAN D H, LAM K C, et al. A new laboratory apparatus for studying dynamic compaction grouting into granular soils[J]. Soils and Foundations, 2013, 53(3): 462-468. doi: 10.1016/j.sandf.2013.04.007
[16]	付艳斌,陈湘生,吕桂阳,等. 基于小孔扩张弹塑性理论的注浆起始劈裂压力研究[J]. 中国公路学报,2020,33(12): 154-163. doi: 10.3969/j.issn.1001-7372.2020.12.012 FU Yanbin, CHEN Xiangsheng, LYU Guiyang, et al. Study on initial fracture pressure of grouting based on the elastic-plastic theory of cavity expansion[J]. China Journal of Highway and Transport, 2020, 33(12): 154-163. doi: 10.3969/j.issn.1001-7372.2020.12.012
[17]	徐芝纶. 弹性力学（上册）[M]. 4版. 北京: 高等教育出版社, 2006.
[18]	姚茂宏,谢长岭,程少振,等. 双孔劈裂条件下地应力及孔间应力的耦合分析[J]. 北京交通大学学报,2021,45(3): 61-68. doi: 10.11860/j.issn.1673-0291.20200140 YAO Maohong, XIE Changling, CHENG Shaozhen, et al. Coupling analysis of in-situ stress and inter hole stress under double-hole splitting[J]. Journal of Beijing Jiaotong University, 2021, 45(3): 61-68. doi: 10.11860/j.issn.1673-0291.20200140
[19]	PANTAZOPOULOS I A, MARKOU I N, CHRISTODOULOU D N, et al. Development of microfine cement grouts by pulverizing ordinary cements[J]. Cement and Concrete Composites, 2012, 34(5): 593-603. doi: 10.1016/j.cemconcomp.2012.01.009
[20]	KANG S P, KWON S J. Effects of red mud and Alkali-activated slag cement on efflorescence in cement mortar[J]. Construction and Building Materials, 2017, 133: 459-467.
[21]	PEPPER R A, COUPERTHWAITE S J, MILLAR G J. Comprehensive examination of acid leaching behaviour of mineral phases from red mud: recovery of Fe, Al, Ti, and Si[J]. Minerals Engineering, 2016, 99: 8-18. doi: 10.1016/j.mineng.2016.09.012
[22]	张淼,邹金锋,陈嘉祺,等. 非对称荷载作用下土体劈裂注浆压力分析[J]. 岩土力学,2013,34(8): 2255-2263. ZHANG Miao, ZOU Jinfeng, CHEN Jiaqi, et al. Analysis of soil fracturing grouting pressure under asymmetric loads[J]. Rock and Soil Mechanics, 2013, 34(8): 2255-2263.
[23]	邹金锋,童无欺,罗恒,等. 基于Hoek-Brown强度准则的裂隙岩体劈裂注浆力学机理[J]. 中南大学学报(自然科学版),2013,44(7): 2889-2896. ZOU Jinfeng, TONG Wuqi, LUO Heng, et al. Mechanism of fracture grouting for fractured rock based on Hoek-Brown failure criterion[J]. Journal of Central South University (Science and Technology), 2013, 44(7): 2889-2896.
[24]	黄明利,管晓明,吕奇峰. 基于弹性力学的诱导劈裂注浆机制分析[J]. 岩土力学,2013,34(7): 2059-2064. HUANG Mingli, GUAN Xiaoming, LÜ Qifeng. Mechanism analysis of induced fracture grouting based on elasticity[J]. Rock and Soil Mechanics, 2013, 34(7): 2059-2064.
[25]	张连震,张庆松,刘人太,等. 基于浆液-岩体耦合效应的微裂隙岩体注浆理论研究[J]. 岩土工程学报,2018,40(11): 2003-2011. ZHANG Lianzhen, ZHANG Qingsong, LIU Rentai, et al. Grouting mechanism in fractured rock considering slurry-rock stress coupling effects[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(11): 2003-2011.
[26]	邹金锋,李亮,杨小礼. 劈裂注浆扩散半径及压力衰减分析[J]. 水利学报,2006,37(3): 314-319. doi: 10.3321/j.issn:0559-9350.2006.03.010 ZOU Jinfeng, LI Liang, YANG Xiaoli. Penetration radius and pressure attenuation law in fracturing grouting[J]. Journal of Hydraulic Engineering, 2006, 37(3): 314-319. doi: 10.3321/j.issn:0559-9350.2006.03.010
[27]	张忠苗,邹健. 桩底劈裂注浆扩散半径和注浆压力研究[J]. 岩土工程学报,2008,30(2): 181-184. ZHANG Zhongmiao, ZOU Jian. Penetration radius and grouting pressure in fracture grouting[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(2): 181-184.
[28]	ZHOU H, LIU H L, ZHA Y H, et al. Influence of stress anisotropy on the cylindrical cavity expansion in undrained elastic-perfectly plastic soil[J]. Science China Technological Sciences, 2018, 61(3): 370-380. doi: 10.1007/s11431-017-9127-7

施引文献

期刊类型引用(12)

1.	邓开来，吴光源，远全闯，谢海清，何庭国. 高铁大跨双肢墩连续刚构桥减震系梁地震易损性分析. 桥梁建设. 2024(03): 61-69 . 百度学术
2.	迟德有. 双幅整体转体跨铁路刚构桥主墩设计方案比选. 桥梁建设. 2024(03): 123-128 . 百度学术
3.	石先林，孟庆利. 某高墩刚构桥抗震性能分析. 西南科技大学学报. 2024(02): 74-82 . 百度学术
4.	吴冠雄. 基于PSO-SVM的大跨PC连续刚构桥施工期可靠度分析. 湖南交通科技. 2024(02): 132-136 . 百度学术
5.	赵灿晖，杨正浩，袁小聪，肖仕周，任小朋，廖文彬. 高桩承台连续刚构桥的调谐隔震设计方法. 中国公路学报. 2024(07): 137-146 . 百度学术
6.	石岩，焦应乾，李军，王文仙，王瑞. 考虑墩-梁碰撞效应的大跨高墩连续刚构桥近断层地震损伤分析. 应用基础与工程科学学报. 2024(04): 1140-1153 . 百度学术
7.	陈俊伟. 减隔震支座对高墩大跨连续刚构桥减震效果分析. 山西交通科技. 2024(04): 62-67 . 百度学术
8.	钟川剑，王戈，蒋建军，何锐. 金阳河特大桥超高桥墩设计关键技术. 世界桥梁. 2024(05): 15-23 . 百度学术
9.	何通，涂金平，邹雨鹤. 高烈度区连续刚构桥弹塑性地震反应及断层走向影响研究. 工程抗震与加固改造. 2024(05): 151-158 . 百度学术
10.	刘仕远，王慧东. 考虑各施工阶段受力影响的连续刚构桥梁地震损伤分析. 工程抗震与加固改造. 2024(05): 159-168+178 . 百度学术
11.	董是，龙志友，毕洁夫，王建伟，邵永军，杨超，左琛，张士远，万昭龙. 基于Savitzky-Golay平滑-小波降噪处理的桥梁结构监测数据分析方法. 测绘通报. 2023(09): 100-106 . 百度学术
12.	刘腾龙，孔君，李震，许梦旗，尹俊红. 桥墩类型对曲线连续刚构桥结构性能影响研究. 山西建筑. 2023(22): 174-176 . 百度学术

其他类型引用(8)

附加材料(0)

访问统计

点击查看大图

图(8) / 表(1)

计量

文章访问数: 502
HTML全文浏览量: 116
PDF下载量: 80
被引次数: 20

1. 固定辙叉廓形特征参数
2. 辙叉钢轨廓形参数化拟合方法
2.1 B样条曲线理论
2.2 关键控制点确定
2.3 廓形拟合
3. 心轨廓形参数试验设计分析方法（DOE）分析
3.1 评价指标
3.2 方案设置
3.3 重构廓形轮轨匹配特征
3.4 DOE分析
4. 结　论

参数	天然密度/（g·cm⁻³）	干密度/（g·cm⁻³）	弹性模量/MPa	泊松比	含水率/%	含泥量/%	注浆影响范围/m
取值	2.01	1.68	5.94	0.38	23.3	11.8	10

考虑浆土应力耦合作用的劈裂注浆机理分析

doi: 10.3969/j.issn.0258-2724.20220433

作者简介: 秦鹏飞（1985—），男，副教授，研究方向为岩土与地下工程，E-mail：929163723@qq.com

计量

出版历程