网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于近场动力学的裂纹扩展及轨面剥离掉块成因机理  PDF

  • 程中宁 1,2
  • 周宇 1,2
  • 李骏鹏 1,2
  • 黄旭炜 3
1. 同济大学 道路与交通工程教育部重点实验室,上海 201804; 2. 同济大学 上海市轨道交通结构耐久与系统安全重点实验室, 上海 201804; 3. 上海市隧道工程轨道交通设计研究院,上海 200235

中图分类号: U213.4+2

最近更新:2023-06-13

DOI:10.11908/j.issn.0253-374x.22045

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

建立了态基近场动力学钢轨表面滚动接触疲劳裂纹扩展模型,研究钢轨疲劳裂纹的扩展以及由此引起的轨面剥离掉块形成过程。结果表明,轨面单裂纹或裂纹间距>1.5 mm情况下,随着通过总重的累积,裂纹扩展的路径逐渐转向钢轨内部,且扩展速率随着通过总重累积和裂纹长度的增加逐渐降低;轨面多裂纹情况下,通过总重增量低于8.8 Mt时,裂纹在荷载作用下主要发生沿长度方向扩展。通过总重增量超过8.8 Mt后,裂纹之间的次生损伤形成;而且,这种情况在轨面多条裂纹间距≤1.5 mm时较为显著,裂纹间距>1.5 mm时,裂纹间的相互影响减小;轨面多条裂纹同时存在且裂纹间距≤1.5 mm时,通过总重增量达到一定程度(>10.0 Mt),裂纹间的次生损伤形成剥离掉块,剥离掉块深度约为裂纹深度的1/2以上。

钢轨表面材料在车轮荷载的反复作用下,容易产生滚动接触疲劳,形成疲劳裂

1,并进一步引发剥离掉2。这类剥离掉块影响超声波探伤,引起轮轨振动和噪声,严重时有导致断轨的风险。

既有钢轨疲劳裂纹扩展与剥离掉块的研

13大多基于断裂力学Paris理论和有限元方法。裂纹和钢轨模型受有限元网格划分影响大,且仅能考虑裂纹尖端处的参数,由此推算一定荷载次数下的扩展速率和角度,其本质是一个静态计算过程。研究认为动态解和静态解有较大差4。此外,有限元方法也较难以处理多裂纹之间互相交叉连通的问5。在裂纹动态扩展研究中,近场动力学理6将物质离散成物质点,通过物质点之间的“键”将模型无网格化;利用积分项代替连续介质力学本构方程中的微分项,计算中对每一根键的损伤程度的更新,符合裂纹的动态变化过程。该理论已经在多晶体断7、结构稳定性和破坏分8、路径扩9等领域得到应用,且键基(bond)近场动力学方法已经用于钢轨疲劳裂纹萌生10,但钢轨材料在荷载作用下发生材料断裂、形成裂纹并扩展和交互,还需要考虑材料点之间的相互影响和作用,引入“态”(state)基近场动力学方法。

本文根据态基近场动力学理论,建立钢轨表面多裂纹疲劳扩展模型,分析钢轨表面疲劳裂纹的动态扩展和剥离掉块的形成,研究不同裂纹间距、裂纹数量在扩展中引起钢轨表面剥离掉块形成过程,为研究钢轨表面疲劳裂纹扩展机理和剥离掉块形成原因提供参考。

1 基于近场动力学的钢轨滚动接触疲劳裂纹扩展

1.1 态基近场动力学方法

近场动力学将连续物体离散为物质点xi,每个物质点xi包含了位置坐标PX,PY,PZ、密度ρ、体积V等信息。物质点xi与其周围一定范围Hxi内的其他任意一物质点xj相互之间存在作用。这样的范围Hxi称作近场作用域(peridynamic horizon)。

上述物质点之间的互相作用,通过“键”(bond)来表达,当键满足一定准则断裂后,即可判定键发生断裂,而穿过这些断裂键所在的位置即是裂纹,如图 1所示。一般以φ(xi)表征物质点xi处发生损伤的程度,即一个物质点xi的断裂键数量占其近场作用域内总的键数量的比,公式如下:

图 1  近场动力学的物质点、作用域与键的断裂

Fig. 1  Material point, horizon and bond breakage in peridynamics

φ(xi)=1-Hμ dVxjHdVxj (1)

式中:μ为一个具有时间历程的标量函数;V为体积;H为作用范围(本文域的取值主要参考Ramulu

11提出的特征距离r0,在允许伤损区定义域大小为0.15 mm,在非裂纹扩展区域不存在尺度效12,大小为1.5 mm)。

考虑到钢轨材料在出现裂纹及其扩展过程中,材料点间的作用力大小不一定相等、方向也不一定一致,进一步引入态基动力

6的“态”(state)来表示此种情况,其中“态”包含了材料点相互作用信息的无限维数组。这样,任意向量ξH在态A̲的运算下的映射是一个m阶张量,记为A̲ξ,可表示为函数

A̲ξ=HLm (2)

式中:Lm为所有m阶张量的集合,其余参数含义同式(1)

态基近场动力学中键之间的关系分变形态和力态,态基近场动力学的物质点运动方程为

ρxiu¨xi,t=Hfxj,xi,tdVxj+bxi,t (3)

式中:ρ为物质点的密度;u¨为加速度矢量;f为点对互相作用的力密度;xj为物质点xi的近场作用域内其他物质点;xi为物质点的位置矢量;b为物质点xi受到的体力密度。

可见,式(3)中将时间、物质点及其之间的关系考虑进来,可以描述物质点相互之间的动态关系,其中力密度函数的具体形式如下:

fxj,xi,t=txj,xi,t-txi,xj,t (4)
txj,xi,t=T̲xi,txj-xi (5)
txi,xj,t=T̲xj,txi-xj (6)

式中:T̲为力态;中为力态的映射对象,表示两物质点之间的位置矢量,也就是联系两个物质点的键ξ,又可以记作T̲ξ,即力态T̲是键ξ和时间t的函数。因此,T̲完整表达写作T̲x,tx'-x

1.2 基于近场动力学的疲劳裂纹扩展方法

采用各向同性线弹性材料疲劳断裂准

13来判断近场动力学模型中键的断裂。假定一个物体受到循环外载作用,而每个键ξ都有剩余寿命λ,初始值设定为1,并且随着循环荷载加载周期N增加而减少。

dλdN(N)=-Aεmλ(0)=1 (7)

式中:Am为疲劳损伤待定参数;ε=s+-s-为键的“循环应变”,其中 s+s-分别为在一次载荷循环中某个键的伸长率的最大值和最小值。

假设在有限次循环周期中,裂纹每一周期恒定向前扩展da/dN长度,a表示裂纹深度,N表示应力循环次数。裂纹尖端位置表达式如下:

z=x-dadNN (8)

式中:x为沿裂尖轴线方向的空间坐标;z为裂纹尖端的空间坐标,当z=0时表示键即将断裂。

记该键的循环应变ε¯(0)=εcore,则推导出裂纹扩展速率同裂纹尖端即将断裂的键循环应变之间的关系式为

dadN=Cεcorem (9)

式中: C为材料常数。

根据线弹性断裂力学,可知εcore与应力强度因子范围ΔK成正比。同时,式(9)与描述疲劳裂纹扩展速度变化规律的Paris公

14相似。因此可直接利用Paris公式获取参数C、m。为获得参数A必须进行一次疲劳试算,然后再根据试算得到的疲劳损伤发展数据进行校准。首先设定任意A值,记作A',施加循环荷载,得到该模型中的裂纹扩展速率(da/dN)',据此进行材料实际疲劳参数A的计算。

A=A'da/dN(da/dN)' (10)

根据近场动力学理论,参数A同近场作用域大小δ的值相关。因此,一般转换为独立于作用域的参数A^,再根据模型情况进行放缩。

A(δ)=A^δm-22 (11)

根据上述理论和方法,建立如下材料疲劳伤损计算方法:

(1) 根据轮轨接触状态,采用式(3)计算钢轨受轮轨接触影响的每个物质点位移。

(2) 选取任一物质点,确定式(7)中的循环应变。

(3) 更新式(7)中键的剩余寿命。

(4) 若剩余寿命小于等于0,则认为键断开。

(5) 遍历该物质点所连接的键,并更新式(1)的损伤度。

(6) 遍历所有物质点,循环(2)~(5)步骤。

(7) 直到循环荷载达到上限值。

以上步骤采用Peridigm计算平台进行二次开发,编写疲劳损伤键断裂和裂纹扩展算法。

1.3 疲劳裂纹扩展建模

1.3.1 二维模型

首先,根据裂纹在钢轨上的空间位

15,以过裂纹在轨面上开口的中心点、沿钢轨纵向的铅垂面为纵剖面,建立含裂纹的钢轨二维模型,如图 2所示。其中,二维平面模型尺寸为长40 mm×高10 mm×厚0.1 mm。整个平面从上至下划分为4个部分:荷载施加区、允许损伤区、常规(材料)区和固定边界区。以二维平面模型的左下角为坐标原点。裂纹按其扩展时与轨面的夹角设置在荷载施加区和允许损伤区。

图 2  钢轨二维平面模型与区域划分

Fig. 2  Rail model and area division in 2 dimensional

根据图 2的二维平面模型区域划分,做如下设定:

(1)对固定边界区内的物质点施加位移约束。

(2)对荷载施加区的物质点施加移动轮轨接触荷载。

(3)允许损伤区物质点之间键的损伤和断裂。

(4)设定组成裂纹的两平面物质点间存在法向接触作用,引入“罚函数”思想,即两裂纹面之间存在有限的穿透量和较大的接触刚度,从而约束裂纹面两侧物质点之间的法向渗透,使得两者渗透较小,且整体刚度矩阵不出现病态,保证计算结果收敛。

(5)采用库伦摩擦模

16描述裂纹面的切向接触属性,以摩擦系数及裂纹面受轮轨接触荷载作用下的法向压力表示裂纹面两侧物质点之间的切向摩擦力。

1.3.2 裂纹数量

裂纹扩展研究发

17,钢轨表面疲劳裂纹常以多裂纹形态出现,本文为做对比,设定钢轨中存在单裂纹和三裂纹工况。每种工况再根据裂纹间距和裂纹长度进行分类,其中裂纹间距分别为1.0、1.5以及2.0 mm,裂纹初始长度为1.5和2.0 mm。

1.3.3 施加荷载

采用车辆‒轨道动力学模型计算车轮在轨面上运行时产生的轮轨接触位置和面积、接触斑黏着‒滑动区分布、接触法向应力和切向应力分布等。根据接触斑中心与轨面裂纹的位置关系在裂纹开口后方施加接触斑法向和切向应力,并使上述接触应力沿钢轨纵向通过裂纹开口。

1.3.4 裂纹的不同位置

对真实裂纹进行观测发

1518,裂纹尖端各点受到的轮轨接触作用不同,因而其扩展情况也相应不同。为了评价裂纹不同尖端位置的扩展情况,使用尖端为半椭圆或样条曲线的平面来模拟裂19,如图 3所示。裂纹开口位于轨面,裂纹向钢轨内部扩展方向与车轮运行方向一致,且裂纹与轨面纵向水平线呈一定角度。这里将裂纹开口长度为1.5 mm、开口与车轮运行方向呈90°、裂纹向钢轨内部扩展角度为30°且具有一定长度和间距的平面状裂纹加入到钢轨三维模型中。再对裂纹沿车轮运行方向纵剖,形成不同裂纹位置的二维断面,分别按1.3节的方法建立含裂纹的二维模型。同时,每个二维模型上对应的轮轨接触应力为该铅垂面所切位置的轮轨接触应力。

图 3  裂纹扩展三维模型

Fig. 3  3-dimensional model for crack propagation

1.4 参数确定

通过车辆‒轨道动力学模型计算轮轨接触荷载。其中,线路为半径500 m曲线;车辆为重载货运C70型敞车,轴重23 t,ZK6转向架,列车速度60 km·h-1,LM磨耗型车轮踏面;钢轨为75 kg⋅m标准型面U75V热处理钢轨,1:40轨底坡,轮轨摩擦系数0.3。接触斑相应参数如表 1所示。

表 1  接触斑参数
Tab. 1  Contact patch parameters
轮轨接触状态参数数值
纵向蠕滑率 6.8×10-4
横向蠕滑率 3.9×10-4
自旋蠕滑率 0.19
纵向蠕滑力(分布力)/MPa 87.5
横向蠕滑力(分布力)/MPa -171.1
法相接触力(分布力)/MPa 1 853.43
接触斑长半轴长/mm 8.3
接触斑短半轴长/mm 3.4
接触斑中心到轨顶中心的距离/mm 6.5

根据轮轨接触荷载,并参考文献[

20-21]的疲劳试验结果,确定的U75V钢轨的Paris公式为

         dadN=2.165 67×10-4+             4.234×10-15(ΔK-650)3.656                ΔK650 MPam1/2 (12)

根据式(12),确定参数m2=3.656。参数A2(式中下标1代表萌生阶段,2代表扩展阶段,)不能从Paris曲线数据直接获得,需要进行一次疲劳试算然后再根据数据校准。该试算根据规

22和文献[23]用单边缺口疲劳试样近场动力学模型来进行。试样长50 mm、宽8 mm、厚3 mm,模型左侧预设单条长5 mm的裂纹缺口,物质点尺寸为0.5 mm×0.5 mm×0.5 mm,上下面物质点层施加均布循环拉荷载σmax=50 MPa,应力循环比R=0.1,疲劳损伤试算参数取A2'=100 000,m2=3.656。计算发现,在荷载循环次数405 017次后,数值求解不再收敛。计算上述模型中疲劳裂纹扩展速率与裂尖应力强度因子,并进行数据拟合,如图 4所示。

图 4  试算模型疲劳裂纹扩展速率与U75V材料Paris曲线

Fig. 4  Fatigue crack propagation rate in test model with U75V rail material Paris curve

根据图 4中的裂纹扩展速率数据与拟合结果,结合1.2节中的式(9),计算U75V钢轨材质对应的独立于物质点尺寸的疲劳损伤参数A^为2 558.224,m为3.656。

1.5 模型验证

根据1.1~1.4节方法,建立如图 2 所示的基于近场动力学的钢轨位移模型,同时建立相同尺寸的有限元模型,材料采用 U75V 钢轨线弹性模型,在上表面中心施加最大接触压力 1 800 MPa,接触斑长半轴 8 mm 的标准 Hertz 理论法向应力,纵向切应力设置为垂向应力数值大小的 0.3倍,模型底部施加固定边界条件,分别计算在相同车辆荷载和轨道条件下两模型的位移结果,如表2所示。

表2  有限元与近场动力学计算的最大、最小位移结果
Tab. 2  Maximum and minimum displacement by the finite element and peridynamics models
最值X方向位移/mmY方向位移/mm
有限元近场动力学相差/%有限元近场动力学相差/%
最大 5.70×10-4 5.03×10-4 -11.75 6.86×10-3 6.60×10-3 -3.79
最小 -4.36×10-2 -3.99×10-2 -8.49 -7.25×10-2 -7.02×10-2 -3.17

图5表2的结果对比来看,近场动力学模型与有限元模型计算结果之间具有较好的一致性。两模型在X方向(钢轨纵向)的位移相差8.49%~11.75%,在Y方向(钢轨横向)的位移相差3.17%~3.79%。

图5  位移云图对比

Fig. 5  Comparison of displacement

2 预测结果与分析

2.1 单条裂纹的动态扩展

以向钢轨内部扩展角度为30°的裂纹为例,不同车轮荷载通过次数下的钢轨接触疲劳裂纹扩展过程如图6所示。每一计算工况都如1.3节中所述包括3个钢轨二维切面,切面坐标以左下角为原点。每一切面上设有初始裂纹,右侧条带表示物质点损伤度。行车方向为从左至右。此外,本文的分析是基于已经有一定长度的预设裂纹,即已经有一定的通过总重,因此,下文所述通过总重均是指在该通过总重下继续积累形成的增量。一般认为当损伤程度>0.385,则物质点之间存在裂纹,为了使结果更加明确,文中的色图中已经将数据经过前处理,将损伤程度低于0.385的部分转化为与损伤程度为0相同的颜色。下图相同。

图6  单裂纹扩展情况

Fig. 6  Single crack propagation

根据图 3的3个切面,确定预设裂纹被切面1、3所切得长度为1.5 mm;被切面2所切得长度为2.0 mm,初始裂纹如图6a所示。

图6b可以看出,当通过总重增量达到约8.8 Mt时(车轮通过次数约N=3.55×105次),疲劳裂纹尖端出现向钢轨内部扩展的现象。图6c反映的是通过总重增量约26.4 Mt时裂纹进一步向钢轨内部扩展,且路径转向与轨面垂直近90°方向,也即在通过总重增量8.8~26.4 Mt时,裂纹扩展的路径开始发生变化。当通过总重增量达到约44.0 Mt时(图6d),裂纹扩展路径进一步朝向钢轨内部并有轻微转向与车轮滚动方向相反方向的趋势(图6d切面2)。进一步,将图6所示的裂纹尖端3处位置的扩展情况结合,得到裂纹尖端三维形态示意图,如图7所示。

图7  裂纹三维扩展及路径变化示意

Fig. 7  3D crack propagation and its route change

尽管图7描述的裂纹尖端点只有3处,但综合图6图7仍然可以看出,单条裂纹情况下,随着通过总重的累积,浅角度裂纹扩展会发生路径变化,转向钢轨内部;由于裂纹尖端各个点受到的接触斑荷载条件不同,尖端各点的发展速率不同,从三维形态来看,裂纹扩展呈现“凸”的形态。

3个切面上疲劳裂纹长度随着通过总重增加的变化,以及裂纹扩展速率分别如图8图9所示。

图8  疲劳裂纹长度变化

Fig. 8  Variation in fatigue crack length

图9  裂纹各切面的阶段扩展速率

Fig. 9  Propagation rate at each section of the crack during different wear phases

图8可以看出,裂纹长度随通过总重累积而逐渐增加,且切面2(裂纹尖端的中点)的裂纹长度在各通过总重时期均大于切面1、3(裂纹尖端中点两侧的点),这与尖端中点上方的轨面承受接触斑中心的荷载有关,即接触荷载大,裂纹的长度长,裂纹沿其尖端呈现“凸”形发展(同图6图7),这与文献[

24]所测情况吻合。

图9可以看出,裂纹扩展速率随着通过总重累积和裂纹长度的增加呈降低趋势,最大扩展速率约0.7 nm·cycle-1;在通过总重在26.4~35.2 Mt阶段时,扩展速率降低到约0.2 nm·cycle-1,并呈上下波动趋势,说明随着裂纹增长,接触斑荷载对裂纹扩展的作用降低,扩展减缓。

在本节预测的通过总重增量情况下,裂纹平均扩展速率约0.33 nm·cycle-1,扩展了约0.5~0.7 mm,转向角度为79°~110°。

2.2 多裂纹的动态扩展

实际情况下,钢轨表面同时会存在多条裂纹,且研究认为多条裂纹在扩展中相互影

17。因此,在模型中考虑多条裂纹的情况,裂纹截取的切面方法同1.3节。裂纹间距分别为1.0、1.5和2.0 mm,裂纹数量为3条。裂纹扩展情况如图 10~ 11所示。

图 10  车轮通过次数N=1.77×105(4.5 Mt)

Fig. 10  Wheel cycles N=1.77×105 (4.5 Mt)

图 11  车轮通过次数N=5.32×105(13.2 Mt)

Fig. 11  Wheel cycles N=5.32×105 (13.2 Mt)

从图 10~ 11可以看出,随着通过总重的累积,裂纹沿其长度方向逐渐向钢轨内部扩展,同时轨面纵向的相邻裂纹之间发生损伤,特别是裂纹间距小于1.5 mm时(多发于硬质钢轨)。在通过总重增量为4.5 Mt时,各裂纹相互之间的影响不明显,但对于间距小的裂纹,裂纹间已经出现损伤迹象(图 10a),特别是裂纹尖端中部位置(图 10a切面2);通过总重增量为13.2 Mt,裂纹间距为1.0~1.5 mm时,裂纹在尖端尚未出现明显的路径转向时,就已经发生裂纹间的损伤,形成裂纹之间的次生裂纹(图 11a、11b),这种次生裂纹与水平面的夹角约45~73°,将会引起剥离掉块。即多裂纹情况下,裂纹间距小于1.5 mm,通过总重累积使得裂纹间的材料先发生破坏。而裂纹间距为2.0 mm时,裂纹间的材料损伤不明显,而裂纹仍然随着通过总重累积向钢轨内部扩展,并在裂纹尖端中点发生路径转变(图 11c切面2)。

可以看出,轨面存在多裂纹情况是实际钢轨使用中出现疲劳裂纹的基本特征,而裂纹间距与钢轨材质有关。随着通过总重的累积,裂纹在荷载作用下扩展,是持续向内部扩展还是出现剥离掉块,则受到裂纹间距(钢轨材质)的影响。

2.3 结果验证

2.3.1 现场取样观测数据

图 11a切面2放大后如图 12a所示,当钢轨轨面存在多裂纹,且裂纹间距≤1.5 mm时,通过总重增量达到一定程度(>10.0 Mt),裂纹间的材料发生水平方向的损伤,损伤深度约为裂纹深度的1/2以上。这样,裂纹间出现损伤区域并发展成次生裂纹,与相邻裂纹串联引起掉块。类似情况也在重载铁路现场取样钢轨的纵断面金相观测中发现,如图 12b、 12c所示。其中图 12b、 12c均为沿列车运行方向(从左向右)的钢轨纵向断面观测图。

图 12  多裂纹情况

Fig. 12  Multi-crack

12b、 12c显示的是同一条曲线的内轨顶面裂纹和掉块情况,对应通过总重分别为60.0 Mt和80.0 Mt

18,可以发现,通过总重增量约20.0 Mt时,裂纹之间已经形成深度约0.45~0.60 mm的剥离掉块。考虑到现场钢轨磨耗速度约0.005 5 mm·Mt-125,出现剥离掉块时磨耗深度约0.22 mm,则加上磨耗量,剥离掉块实际深度为0.67~0.82 mm。而本文预设的裂纹发展到剥离掉块的深度为0.70~1.10 mm,对应的通过总重增量约13.2~17.6 Mt。

进一步对通过总重60.0~150.0 Mt范围的重载铁路曲线内轨表面试块进行金相观测,可以观测到从主裂纹扩展出分枝裂纹或是在裂纹间出现次生裂纹的现象,如图 13 所示。从图 13a可以看出,相邻裂纹间的次生裂纹,呈与裂纹扩展角度更小、近似平行钢轨表面的角度发展,扩展至一定长度后与相邻裂纹串联,导致裂纹上方的材料与钢轨整体的连接减少,最终脱离钢轨表面母材。当裂纹范围内大面积区域发生材料破坏,就形成宏观的钢轨表面剥离掉块,多条裂纹的存在促进了这一过程。对现场取样的钢轨表面疲劳损伤发展观测中(图 13b)也能够观察到分枝裂纹与次生裂纹的出现和发展,与仿真形态相似。

图 13  次生裂纹情况

Fig. 13  Secondary crack

2.3.2 与现有文献的对比

文献[

26]的重载铁路钢轨裂纹扩展,线路条件与本文仿真的线路条件相近。上述文献和本文的预测结果如表 3所示。

表 3  预测结果对比
Tab. 3  Comparison of forecast results
对比钢轨材质扩展速率/(nm·cycle-1扩展角度/(°)
文献[26 U75V 0.325 145~150
本文 U75V 0.270~0.340 79~110

表 3可以看出,本文提出的基于态基近场动力学的疲劳裂纹扩展方法,预测的裂纹扩展速率和裂纹扩展角度与观测结果相近。

3 结论

(1)轨面单裂纹(或裂纹间距>1.5 mm)情况下,通过总重增量8.8~26.4Mt时,裂纹扩展的路径开始发生变化,裂纹有向钢轨内部扩展的趋势,且裂纹扩展速率随着通过总重累积和裂纹长度的增加逐渐降低。

(2)轨面多裂纹情况下,通过总重增量低于8.8 Mt时,轮轨接触荷载为裂纹扩展的主要推动因素,裂纹沿长度方向扩展,但未发生剥离掉块现象。通过总重增量超过8.8 Mt后,裂纹在沿长度方向扩展的同时,裂纹之间的次生损伤形成,同时,裂纹间距会直接影响裂纹是否通过次生损伤发生交叉联通形成剥离掉块。

(3)轨面多条裂纹同时存在且裂纹间距≤1.5 mm时,在裂纹扩展过程中裂纹间出现次生裂纹和损伤。裂纹间距大于>1.5 mm,裂纹间的相互影响减小。

(4)轨面多条裂纹同时存在且裂纹间距≤1.5 mm时,通过总重增量达到一定程度(>10.0 Mt),形成剥离掉块,剥离掉块深度(0.7~1.10 mm)约为裂纹深度的1/2以上,且剥离掉块深度随裂纹间距增大(裂纹间距不超过1.5 mm)呈增大趋势。

作者贡献声明

程中宁:论文撰写,仿真方法设计,结果可视化。 周 宇:论文指导,文本内容校准。 李骏鹏:电镜图片收集整理。 黄旭炜:仿真数据整理与管理。

参考文献

1

CANADINC DSEHITOGLU HVERZAL K. Analysis of surface crack growth under rolling contact fatigue[J]. International Journal of Fatigue2008309): 1678. [百度学术] 

2

张斌付秀琴. 铁路车轮、轮箍踏面剥离的类型及形成机理[J]. 中国铁道科学20012): 76. [百度学术] 

ZHANG BinFU Xiuqin. Type and formation mechanism of railway wheel and tire tread spall[J]. China Railway Science20012): 76. [百度学术] 

3

MAI S HGRAVOUIL ANGUYEN-TAJAN M Let al. Numerical simulation of rolling contact fatigue crack growth in rails with the rail bending and the frictional contact[J]. Engineering Fracture Mechanics2017174196. [百度学术] 

4

ZHAO XZHAO XLIU Cet al. A study on dynamic stress intensity factors of rail cracks at high speeds by a 3D explicit finite element model of rolling contact[J]. Wear2016366/36760. [百度学术] 

5

黄丹章青乔丕忠. 近场动力学方法及其应用[J]. 力学进展2010404): 448. [百度学术] 

HUANG DanZHANG QingQIAN Pizhonget al. A review on peridynamics(PD) method and its applications[J]. Advances in Mechanics2010404): 448. [百度学术] 

6

SILLING S AEPTON MWECKNER Oet al. Peridynamic states and constitutive modeling[J]. Journal of Elasticity2007882): 151. [百度学术] 

7

ASKARI EBOBARU FLEHOUCQ R Bet al. Peridynamics for multiscale materials modeling[J]. Journal of Physics Conference200812512078. [百度学术] 

8

KILIC BMADENCI E. Structural stability and failure analysis using peridynamic theory[J]. International Journal of Non-Linear Mechanics2009448): 845. [百度学术] 

9

SILLING S AWECKNER OASKARI Eet al. Crack nucleation in a peridynamic solid[J]. International Journal of Fracture20101621): 219. [百度学术] 

10

马晓川刘林芽张鹏飞. 近场动力学框架下钢轨疲劳裂纹萌生预测的数值方法研究[J]. 摩擦学学报2020405): 608. [百度学术] 

MA XiaochuanLIU LinyaZHANG Pengfeiet al. Numerical method for predicting rail fatigue crack initiation with peridynamic theory[J]. Tribology2020405): 608. [百度学术] 

11

RAMULU MKOBAYASHI A SKANG Bet al. Further studies on dynamic crack branching[J]. Experimental Mechanics1983234): 431. [百度学术] 

12

BOBARU FHU W. The meaning, selection, and use of the peridynamic horizon and its relation to crack branching in brittle materials[J]. International Journal of Fracture20121762): 215. [百度学术] 

13

SILLING S AASKARI A. Peridynamic model for fatigue cracking[R]. AlbuquerqueSandia National Laboratories2014. [百度学术] 

14

PARIS P. A critical analysis of crack propagation laws[J]. J Basic Eng1963854): 528. [百度学术] 

15

ZHOU YZHENG XJIANG Jet al. Modeling of rail head checks by x-ray computed tomography scan technology[J]. International Journal of Fatigue201710021. [百度学术] 

16

KALKER J J. A fast algorithm for the simplified theory of rolling contact[J]. Vehicle System Dynamics1982111): 1. [百度学术] 

17

周宇孙鼎人王树国. 钢轨滚动接触疲劳多裂纹相互影响[J]. 华东交通大学学报2020372): 7. [百度学术] 

ZHOU YuSUN DingrenWANG Shuguoet al. Interactive influence on multiple head checks propagation in rail[J]. Journal of East China Jiaotong University2020372): 7. [百度学术] 

18

ZHOU YWANG SWANG Tet al. Field and laboratory investigation of the relationship between rail head check and wear in a heavy-haul railway[J]. Wear20143151): 68. [百度学术] 

19

周宇邝迪峰郑晓峰. 基于三维重构的钢轨滚动接触疲劳裂纹扩展预测[J]. 机械工程学报20184): 158. [百度学术] 

ZHOU YuKUANG DifengZHENG Xiaofenget al. Prediction of the rail head checks propagation based on three dimensional reconstruction[J]. Journal of Mechanical Engineering20184): 158. [百度学术] 

20

吕斌. U75V钢轨疲劳裂纹扩展行为研究[D]. 杭州浙江工业大学2013. [百度学术] 

LÜ Bin. Research on fatigue crack propagation behavior of U75V rail steel[D]. HangzhouZhejiang University of Technology2013. [百度学术] 

21

王文健刘启跃. 轮轨滚动接触疲劳与磨损耦合关系及预防措施研究[J]. 中国铁道科学2009304): 137. [百度学术] 

WANG WenjianLIU Qiyue. Study on the coupling relationship between the rolling contact fatigue and the wear of wheel-rail and prevention measures[J]. China Railway Science2009304): 137. [百度学术] 

22

钢铁研究总院北京航空材料研究院冶金工业信息标准研究院. 金属材料 疲劳试验 疲劳裂纹扩展方法GB/T 6398—2017[S]. 北京中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会2017. [百度学术] 

Institute Central Iron Steel Research, Materials Bejing Institute of Aeronautical, Institute China Metallurgical Information and Standardization Institute. Metallic materials-fatigue testing-fatigue crack growth methodGB/T 6398—2017[S]. BeijingGeneral Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration of the People’s Republic of China2017. [百度学术] 

23

FREIMANIS AKAEWUNRUEN S. Peridynamic analysis of rail squats[J]. Applied Sciences2018811): 2299. [百度学术] 

24

陈朝阳张银花刘丰收. 朔黄铁路曲线下股热处理钢轨剥离伤损成因分析[J]. 中国铁道科学20084): 28. [百度学术] 

CHEN ZhaoyangZHANG YinhuaLIU Fengshouet al. Analysis on the formation cause of spalling and damage of the heat-treated low rail on Shuohuang railway curve[J]. China Railway Science20084): 28. [百度学术] 

25

周宇姜俊楠张杰. 重载铁路曲线内轨疲劳裂纹与磨耗关系现场试验分析[R]. 北京中国铁道学会2014. [百度学术] 

ZHOU YuJIANG JunnanZHANG Jieet al. Field test and analysis on the relationship between fatigue crack and wear of curved inner rail on heavy haul railway[R]. BeijingChina Railway Society2014. [百度学术] 

26

木东升. 钢轨滚动接触疲劳弯钩形裂纹扩展路径变化研究[D]. 上海同济大学2019. [百度学术] 

MU Dongsheng. Research on the propagation route change of the hook-shaped rolling contact fatigue crack in the rail[D]. ShanghaiTongji University2019. [百度学术]