黏土脏污对道砟集料的应力-剪胀关系影响

陈静; 高睿; 刘洋泽鹏; 石知政; 张荣隆

doi:10.3969/j.issn.0258-2724.20200627

黏土脏污对道砟集料的应力-剪胀关系影响

doi: 10.3969/j.issn.0258-2724.20200627

陈静^{1, 2,},
高睿^{1, 2, ,},
刘洋泽鹏^{1, 2},
石知政^{1, 2},
张荣隆^{1, 2}

1.
武汉大学土木建筑工程学院，湖北武汉 430072
2.
武汉大学岩土与结构工程安全湖北省重点实验室，湖北武汉 430072

基金项目: 国家自然科学基金（51178358，51878521）

详细信息

作者简介:
陈静（1994—），女，博士研究生，研究方向为道路与铁道工程，E-mail：chenjing_whu@hotmail.com

通讯作者:
高睿（1975—），男，教授，博士，研究方向为道路与铁道工程，E-mail：gaorui@whu.edu.cn

中图分类号: U213.7
计量
- 文章访问数: 389
- HTML全文浏览量: 144
- PDF下载量: 18
- 被引次数: 5
出版历程
- 收稿日期: 2020-09-14
- 修回日期: 2021-04-21
- 网络出版日期: 2022-08-03
- 刊出日期: 2021-04-29

Effect of Clay Contamination on Stress-Dilatancy Relationships of Ballast Aggregate

CHEN Jing^{1, 2
,},
GAO Rui^{1, 2
, ,},
LIU Yangzepeng^{1, 2},
SHI Zhizheng^{1, 2},
ZHANG Ronglong^{1, 2}

1.
School of Civil Engineering, Wuhan University, Wuhan 430072, China
2.
Key Laboratory of Safety for Geotechnical and Structural Engineering of Hubei Province, Wuhan University, Wuhan 430072, China

摘要

摘要:
列车循环荷载下，铁路地基层中的黏土颗粒会逐渐侵入道砟层，降低道床承载性能. 通过开展一系列大型直剪试验，对黏土脏污质污染下的受格栅加筋道砟集料和未受格栅加筋道砟集料的强度及变形发展、应力-剪胀关系进行了研究. 试验结果表明：随着脏污程度的增加，道砟集料在剪切过程中强度及法向变形均降低；干净道砟集料的应力比与剪胀比呈一阶线性关系，而黏土污染的道砟集料塑性增强，在剪应力峰值状态下剪胀比更高而抗剪强度发展缓慢，导致应力比与剪胀比呈二阶多项式关系；高法向压力下，道砟集料具有更低的剪胀比；黏土脏污质的存在会减小道砟集料的变形破坏速率，但降低了集料的抗剪强度，而格栅的加筋作用可弥补由于黏土脏污引起的强度降低.
- 格栅加固 /
- 道床脏污 /
- 直剪试验 /
- 应力-剪胀关系
Abstract:
Clay fines from subgrade would gradually intrude into the ballast layer under cyclic loadings of passing trains, which would reduce the bearing capacity and impede the free drainage of track beds. A series of large-scale direct shear tests were carried out to investigate the strength and deformation characteristics and stress-dilatancy relationship of the geogrid-reinforced and unreinforced ballast contaminated by clay fines. The results showed that the strength and normal displacement of ballast aggregate decrease with an increase in the contamination level. The stress ratio of clean ballast is linear with the dilatancy ratio, while the addition of clay fines would increase the plasticity of the ballast aggregate. For fouled ballast in the peak state of shear stress, the dilatancy ratio of aggregate increases while the shear strength decreases, and a second-order polynomial relationship between stress ratio and dilatancy ratio is observed. Under a higher normal pressure, the aggregates have a lower dilatancy ratio. The reduction in the dilatancy rate and the shear strength of clay-contaminated ballast can be remedied by the inclusion of geogrid in the aggregate.
- geogrid reinforcement /
- ballast fouling /
- direct shear test /
- stress-dilatancy relationship

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. Apparatus of large-scale direct shear tests

下载: 全尺寸图片幻灯片

图 2 道砟及黏土脏污质级配图

Figure 2. Size distribution of ballast and clay fine

下载: 全尺寸图片幻灯片

图 3 试验所用材料

Figure 3. Materials for the large-scale direct shear test

下载: 全尺寸图片幻灯片

图 4 受格栅加筋与未受格栅加筋道砟τ_s及δ_n与δ_s关系

Figure 4. Relationships of shear stress τ_s vs. shear displacement δ_s and relationships of normal displacement δ_n vs. shear displacement δ_s for geogrid-reinforced and unreinforced ballast

下载: 全尺寸图片幻灯片

图 5 τ_s/σ_n-δ_n/δ_s关系

Figure 5. Relationships between τ_s/σ_n and δ_n/δ_s

下载: 全尺寸图片幻灯片

图 6 τ_s/σ_n-δ_n/δ_s关系曲线

Figure 6. Relationships between τ_s/σ_n and δ_n/δ_s

下载: 全尺寸图片幻灯片

表 1 道砟及黏土脏污质工程性能指标表

Table 1. Engineering properties of clean ballast and clay fines

材料	比重	堆积密度/ （kg·m⁻³）	e	ω_L/ %	ω_P/ %	D₅₀/ mm	C_u	C_c
道砟	2.66	1432	0.858			39.0	1.15	1.75
黏土	2.70	1178	1.207	42.10	22.40	1.3	24.29	1.34

下载: 导出CSV

表 2 格栅参数表

Table 2. Mechanical properties of geogrid

参数	规格
孔径形式	双向孔径 55 mm × 55 mm
材料	聚丙烯
2% 应变时抗拉强度/（kN·m⁻¹）	11
5% 应变时抗拉强度/（kN·m⁻¹）	15
峰值抗拉强度/（kN·m⁻¹）	30
屈服点伸长率/%	13

下载: 导出CSV

表 3 试样设置情况表

Table 3. Details of experimental specimens

VCI/%	试样编号	格栅布置	黏土质量/kg
0	S1-1	格栅加筋
0	S1-2	格栅加筋
20	S2-1	格栅加筋	21.70
20	S2-2	格栅加筋	21.70
40	S3-1	格栅加筋	43.40
40	S3-2	格栅加筋	43.40

下载: 导出CSV

参考文献(27)

[1]	INDRARATNA B, SINGH M, NGUYEN T T. The mechanism and effects of subgrade fluidisation under ballasted railway tracks[J]. Railway Engineering Science, 2020, 28(2): 113-128. doi: 10.1007/s40534-020-00210-1
[2]	TENNAKOON N, INDRARATNA B. Behaviour of clay-fouled ballast under cyclic loading[J]. Géotechnique, 2014, 64(6): 502-506.
[3]	INDRARATNA B, TENNAKOON N, NIMBALKAR S, et al. Behaviour of clay-fouled ballast under drained triaxial testing[J]. Géotechnique, 2013, 63(5): 410-419.
[4]	NGO T, INDRARATNA B. Analysis of deformation and degradation of fouled ballast: experimental testing and DEM modeling[J]. International Journal of Geomechanics, 2020, 20(9): 06020020.1-06020020.8.
[5]	HUANG H, TUTUMLUER E, DOMBROW W. Laboratory characterization of fouled railroad ballast behavior[J]. Transportation Research Record: Journal of the Transportation Research Board, 2009, 2117(1): 93-101. doi: 10.3141/2117-12
[6]	DANESH A, PALASSI M, MIRGHASEMI A A. Effect of sand and clay fouling on the shear strength of railway ballast for different ballast gradations[J]. Granular Matter, 2018, 20(3): 1-14.
[7]	TRINH V N, TANG A M, CUI Y J, et al. Mechanical characterisation of the fouled ballast in ancient railway track substructure by large-scale triaxial tests[J]. Soils and Foundations, 2012, 52(3): 511-523. doi: 10.1016/j.sandf.2012.05.009
[8]	ROWE P W. The stress-dilatancy relation for static equilibrium of an assembly of particles in contact[J]. Proceedings of the Royal Society of London Series A Mathematical and Physical Sciences, 1962, 269(1339): 500-527.
[9]	MITCHELL J K, SOGA K. Fundamentals of Soil Behavior[M]. Manhattan: John Wiley & Sons, 1976.
[10]	XIAO Y, MENG M Q, CHEN Q S, et al. Friction and dilatancy angles of granular soils incorporating effects of shearing modes[J]. International Journal of Geomechanics, 2018, 18(11): 06018027.1-06018027.9.
[11]	XIAO Y, LIU H L, LIU H, et al. Strength and dilatancy behaviors of dense modeled rockfill material in general stress space[J]. International Journal of Geomechanics, 2016, 16(5): 04016015.1-04016015.15.
[12]	DOŁŻYK-SZYPCIO K. Stress-dilatancy relationship for railway ballast[J]. Studia Geotechnica et Mechanica, 2018, 40(2): 79-85. doi: 10.2478/sgem-2018-0018
[13]	SZYPCIO Z. Stress-dilatancy of gravel for triaxial compression tests[J]. Annals of Warsaw University of Life Sciences-SGGW. Land Reclamation, 2018, 50(2): 119-128. doi: 10.2478/sggw-2018-0010
[14]	DOŁŻYK-SZYPCIO K. Direct shear test for coarse granular soil[J]. International Journal of Civil Engineering, 2019, 17(12): 1871-1878. doi: 10.1007/s40999-019-00417-2
[15]	SAROJINIAMMA B K, INDRARATNA B, VINOD J S. A semi-empirical dilatancy model for ballast fouled with plastic fines[J]. Geomechanics and Geoengineering, 2019, 14(1): 12-17. doi: 10.1080/17486025.2018.1476737
[16]	陈静,高睿,刘洋泽鹏,等. 不同脏污质对格栅加筋道砟性能的影响[J]. 西南交通大学学报,2022,57(1): 200-206. CHEN Jing, GAO Rui, LIU Yangzepeng, et al. Influence of various fouling materials on geogrid-reinforced ballast performance[J]. Journal of Southwest Jiaotong University, 2022, 57(1): 200-206.
[17]	铁道科学研究院. 铁路碎石道砟: TBT 2140—2008. 北京: 中国铁道出版社, 2008.
[18]	MCDOWELL G R, HARIRECHE O, KONIETZKY H, et al. Discrete element modelling of geogrid-reinforced aggregates[J]. Proceedings of the Institution of Civil Engineers−Geotechnical Engineering, 2006, 159(1): 35-48. doi: 10.1680/geng.2006.159.1.35
[19]	INDRARATNA B, SU L J, RUJIKIATKAMJORN C. A new parameter for classification and evaluation of railway ballast fouling[J]. Canadian Geotechnical Journal, 2011, 48(2): 322-326. doi: 10.1139/T10-066
[20]	SELIG E T, WATERS J M. Track geotechnology and substructure management[M]. London: Homas Telford Publishing, 1994
[21]	FELDMAN F, NISSEN D. Alternative testing method for the measurement of ballast fouling: percentage void contamination[C]//CORE 2002, Cost Efficient Railways Through Engineering, Conference on Railway Engineering. Wollongong: [s.n.], 2002: 10-13.
[22]	BIABANI M M, INDRARATNA B, NIMBALKAR S. Assessment of interface shear behaviour of sub-ballast with geosynthetics by large-scale direct shear test[J]. Procedia Engineering, 2016, 143: 1007-1015. doi: 10.1016/j.proeng.2016.06.094
[23]	高睿,石知政,刘洋泽鹏,等. 土工格栅对受污道砟直剪特性影响的试验研究[J]. 西南交通大学学报,2021,56(6): 1185-1191. GAO Rui, SHI Zhizheng, LIU Yangzepeng, et al. Experimental study on effect of geogrid on direct shear behavior of contaminated ballast[J]. Journal of Southwest Jiaotong University, 2021, 56(6): 1185-1191.
[24]	INDRARATNA B, NGO N T, RUJIKIATKAMJORN C. Behavior of geogrid-reinforced ballast under various levels of fouling[J]. Geotextiles and Geomembranes, 2011, 29(3): 313-322. doi: 10.1016/j.geotexmem.2011.01.015
[25]	SUN Y F, GAO Y F, CHEN C. Critical-state fractional model and its numerical scheme for isotropic granular soil considering state dependence[J]. International Journal of Geomechanics, 2019, 19(3): 04019001.1-04019001.14.
[26]	BOLTON M D. The strength and dilatancy of sands[J]. Géotechnique, 1986, 36(1): 65-78.
[27]	INDRARATNA B, NGO N T, RUJIKIATKAMJORN C, et al. Behavior of fresh and fouled railway ballast subjected to direct shear testing: discrete element simulation[J]. International Journal of Geomechanics, 2014, 14(1): 34-44. doi: 10.1061/(ASCE)GM.1943-5622.0000264

施引文献

期刊类型引用(2)

1.	马少坤，田发派，黄海均，张加兵，段智博，龚健. 地铁循环荷载下黏土污染圆砾的大型动三轴试验及动力特性研究. 振动与冲击. 2024(07): 245-254 . 百度学术
2.	王萌，肖源杰，张冲冲，杨涛，谭攀，卢明皎. 不同煤灰脏污程度下道砟宏观强度及细观颗粒转动特性分析. 岩石力学与工程学报. 2024(08): 2027-2041 . 百度学术

其他类型引用(3)

附加材料(0)

访问统计

点击查看大图

图(6) / 表(3)

计量

文章访问数: 389
HTML全文浏览量: 144
PDF下载量: 18
被引次数: 5

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

黏土脏污对道砟集料的应力-剪胀关系影响

doi: 10.3969/j.issn.0258-2724.20200627

作者简介: 陈静（1994—），女，博士研究生，研究方向为道路与铁道工程，E-mail：chenjing_whu@hotmail.com

通讯作者: 高睿（1975—），男，教授，博士，研究方向为道路与铁道工程，E-mail：gaorui@whu.edu.cn

计量

出版历程