路基冻胀地区CRTS I型板式无砟轨道结构变形与离缝特征分析
蔡小培, 梁延科, 谭诗宇, 沈宇鹏
北京交通大学 土木建筑工程学院,北京 100044

第一作者:蔡小培(1982—),男,江苏徐州人,副教授,博士,博士生导师.研究方向为铁路轨道设计理论与试验技术.email: caixiaopei@126.com

摘要

为揭示严寒地区高速铁路路基冻胀变形对无砟轨道平顺性的影响,建立CRTS I型板式无砟轨道-路基空间耦合有限元模型,分析了不同路基冻胀条件下轨道结构的变形特征,探讨了层间离缝的发展演变过程,以及层间粘结强度和底座板刚度对离缝发展的影响规律.结果表明:路基冻胀位置对轨道结构变形影响较大,冻胀变形基本能反映到轨面,当冻胀作用在轨道板中间位置时,底座板与基床表层之间的离缝值最大;最大离缝值随冻胀量增加呈线性增长,随冻胀波长的增大而减小;离缝在轨道结构横向位置是从中间向两边逐渐扩展的;随着底座板与基床表层之间粘结强度的增大,层间离缝值和离缝长度逐渐减小;离缝值随着底座板刚度的减小而减小,当底座板刚度减小为原来的60%时,离缝值减小了近20%.

关键词: 无砟轨道; 路基冻胀; 有限元模型; 层间离缝; 粘结强度
中图分类号:U213.2 文献标志码:A
Deformation and seam characteristics analysis of CRTSⅠslab ballastless track in subgrade frost heaving zone
CAI Xiaopei, LIANG Yanke, TAN Shiyu, SHEN Yupeng
School of Civil Engineering, Beijing Jiaotong University, Beijing 100044,China
Abstract

In order to reveal the effects of subgrade frost heaving deformation on riding comfort of high-speed railway ballastless track in severe cold area, a spatial coupling finite element model for a CRTS I ballastless track-subgrade system was established. Based on the model, the deformation characteristics of track structure under different frost heaving conditions were analyzed. Besides, the developing process of interlayer seam was studied, the influences of coherence strength between layers, and bed plate stiffness on the development of layer seam were explored. The results show that frost heaving position has a big effect on track structure deformation, and the frost heaving deformation can reflect to rail surface. The maximum value of seam between bed plate and bedding surface occurs when frost heaving acts on the middle of track plate. The maximumvalue of the seam increases linearly with the increase of frost heaving amount, and it decreases with the increase of frost heaving wavelength. The seam extends from the middle to both sides in the transverse direction of trade structure. Both the value and the length of seam diminish as the coherence strength between bed plate and bedding surface increases. The value of seam decreases with the decrease of bed plate stiffness. Specifically, the value of seam decreases by 20% when bed plate stiffness is 60% of the initial.

Keyword: ballastless track; subgrade frost heaving; finite element model; inter layer seam; coherence strength

随着高速铁路的快速发展, 我国严寒地区的高速铁路里程已达3 000 km[1].在严寒的西北和东北地区, 路基冻胀的问题普遍存在[2, 3], 导致无砟轨道的高低出现异常.在列车荷载长时间作用下, 高低不平路基上的无砟轨道结构很容易在底座板和基床表层之间产生层间离缝, 即底座板与基床表层之间粘结强度失效, 影响轨道结构受力, 造成轨道不平顺, 给列车的安全平稳运行及线路设备耐久性带来不利影响.因此, 有必要研究路基冻胀区无砟轨道结构变形及离缝特征.

在对轨道结构层间离缝及路基冻胀的研究中, 杨俊斌等[4]分析了I型板端不同离缝长度和高度对轨道结构及列车力学性能的影响; 邵丕彦等[5]分析了I型板在不同温度下轨道板的翘曲变形和CA砂浆层离缝变化情况; Ma Fuxun等[2]利用GPS分析了在温度变化下路基冻胀的变化规律; 赵国堂[1]针对严寒地区的路基冻胀进行特征分析, 提出了路基冻胀管理标准的确定方法.

然而, 现有研究成果大都集中在离缝对轨道结构和列车的影响及温度作用下离缝和路基冻胀的机理研究.离缝发生位置主要考虑在CA砂浆层, 针对路基冻胀区底座板与基床表层之间离缝的研究成果较少, 模型很少考虑底座板伸缩缝及重力的作用, 且常将层间接触关系进行简化, 以上均会影响层间离缝破坏机理的研究.但底座板与基床表层之间的离缝对轨道结构及列车有着不容忽视的影响, 所以对不同路基冻胀条件下轨道结构变形及离缝的研究具有理论意义和实际的工程意义.

本文作者运用ABAQUS有限元分析软件, 以我国严寒地区应用较多的CRTS I型板无砟轨道为研究对象, 建立无砟轨道-路基空间耦合有限元模型, 分析了局部冻胀地区冻胀位置、冻胀量、冻胀波长、层间粘结强度及底座板刚度对轨道结构变形的影响, 揭示了不同路基冻胀条件下无砟轨道结构变形及离缝的发展规律.

1 有限元模型的建立
1.1 有限元模型

CRTSⅠ 型板式无砟轨道主要由钢轨、扣件系统、轨道板、CA砂浆层、凸型挡台、环形树脂、混凝土底座板和基床组成[6], 其横断面见图1.为了更好地反映工程实际情况, 除扣件外, 其他部件均选用实体单元.所建立的无砟轨道-路基空间耦合有限元模型的长度为60 m.建立的有限元模型见图2.

图1 CRTSⅠ 型板式无砟轨道横断面图(单位:mm)Fig.1 Cross section of CRTSⅠ slab ballastless track(unit:mm)

图2 CRTSⅠ 型板式无砟轨道-路基空间耦合有限元模型Fig.2 Spatial coupling finite element model for CRTSⅠ slab ballastless track-subgrade system

有限元模型中, 钢轨为我国高速铁路采用的60 kg/m钢轨, 扣件为J-7B型扣件, 扣件采用连接钢轨与轨道板对应单元节点的线性弹簧模拟, 并考虑了纵、横、垂3个方向的刚度, 约束弹簧下部端点的转动, 以模拟扣件垫板的作用.扣件间距为0.629 m, 纵、横、垂3个方向的刚度分别为30 kN/mm、30 kN/mm、50 kN/mm.

对于严寒地区选用的CRTSⅠ 型板式无砟轨道, 底座板伸缩缝的存在会对离缝的产生与发展产生不容忽视的影响[1], 所以本文所建立的模型将重点考虑伸缩缝.工程实际中, 严寒地区底座板每隔2块或3块轨道板设置1条伸缩缝, 本文混凝土底座板每隔2块轨道板设置1条伸缩缝, 伸缩缝宽度为20 mm.

无砟轨道结构各部件几何尺寸按工程实际考虑, 轨道各结构的材料参数见表1.

表1 轨道结构材料特性 Tab.1 Material properties of track components
1.2 接触、边界条件及荷载定义

研究底座板与基床表层之间的离缝时, 将轨道结构本身视为整体, 故在定义接触时, 不考虑轨道板与CA砂浆层、CA砂浆层与混凝土底座板、环形树脂与凸型挡台之间的相对位移, 采用绑定约束连接[7]; 底座板与基床表层之间的接触考虑为摩擦系数为0.5的硬性接触, 即两接触面在压紧状态下会传递法向力, 并会产生相对的滑移[8, 9, 10], 这样有利于模拟离缝的产生和发展, 且符合工程实际情况.

边界设置如下:钢轨两端约束纵向位移, 以模拟无缝线路; 底座板两端约束纵向位移; 基床底部则定义为全约束.

实际工程中, 在轨道结构重力作用下, 路基会产生累积变形, 最终达到平衡状态.为了避免基床表层在计算过程中产生过大变形, 模拟路基的真实状态, 进行了地应力平衡, 同时考虑了轨道各结构和基床的重力作用.基床冻胀变形考虑为余弦曲线函数[1], 运用ABAQUS中的displacement对基床表层施加变形荷载, 模拟基床的冻胀作用.

2 路基冻胀的影响
2.1 冻胀位置的影响

路基冻胀作用位置考虑3种情况, 如图3所示.a位置为底座板伸缩缝处, b位置为轨道板中间, c位置为两伸缩缝中间.计算中, 冻胀考虑为路基局部冻胀, 冻胀曲线选取波长5 m, 冻胀量20 mm的半余弦曲线, 不考虑基床表层的粘结强度, 且底座板刚度保持不变.

图3 冻胀作用位置Fig.3 Frost heaving positions

图4是作用在不同位置时底座板与基床表层离缝值变化情况.由图4可知, 由于伸缩缝的存在, 冻胀作用位置对底座板变形影响较大, 当冻胀分别作用在a位置、b位置和c位置时, 最大离缝值分别为11.26 mm、21.52 mm、14.46 mm, 可见冻胀发生在轨道板中间位置时, 离缝值最大, 约为发生在伸缩缝处时的2倍.

图4 冻胀位置对离缝值的影响Fig.4 Influence of frost heaving positions on seam value

图5是不同冻胀作用位置轨道结构应力云图.由图5可知, 冻胀作用在a位置时, 轨道结构基本不会出现损伤; 冻胀作用在b位置时, 由于轨道结构的层间粘结作用, 使冻胀波峰处单块轨道板产生较大拉应力; 冻胀作用在c位置时, 冻胀作用对轨道结构拉应力影响范围最广, 在冻胀波峰位置处的相邻2块轨道板均产生了较大的拉应力, 超过混凝土的抗拉强度1.89 MPa, 轨道结构最易出现损伤.从离缝大小的角度考虑, 冻胀作用最不利位置为轨道板中间处; 从轨道结构应力考虑, 冻胀作用最不利位置为两伸缩缝中间处.选取冻胀量20 mm, 波长10 m, 位置为两伸缩缝中间处, 进行变形传递规律的研究, 得出冻胀作用下各层位移的变化如图6所示.由图6可知, 路基冻胀变形基本能传递到轨面, 轨道不平顺波形与路基冻胀波形基本一致.考虑到实际工程中轨道结构的破坏损伤, 选取冻胀位置为两伸缩缝中间处, 进行冻胀波长、冻胀量、粘结强度及底座板刚度的影响分析.

图5 冻胀作用位置对轨道结构应力的影响Fig.5 Influence of frost heaving positions on track structure stress

图6 冻胀作用下各层位移图Fig.6 Displacements diagram of each layer under the effect of frost heaving

2.2 冻胀波长的影响

根据工程实际情况, 冻胀量取20 mm, 不考虑基床表层的粘结强度, 且底座板刚度保持不变, 分析冻胀波长对离缝的影响, 冻胀波长分别取5 m、10 m、15 m、20 m、30 m.

冻胀波长对最大离缝值的影响见图7.由图7可知, 随着冻胀波长的增加, 最大离缝值呈不断减小的趋势, 且变化速率逐渐减小, 冻胀波长大于30m时, 最大离缝值小于1.0 mm, 当波长为40 m时, 基本没有离缝.可见, 波长40 m以上的冻胀不会导致底座板和基床表层间产生离缝.除冻胀波长为5 m和10 m的最大离缝值出现在冻胀波峰附近外, 其他波长冻胀产生的最大离缝值都出现在伸缩缝处.

图7 冻胀波长对最大离缝值的影响Fig.7 Influence of frost heaving wavelengths on the maximum seam value

2.3 冻胀量的影响

根据实际情况, 冻胀波长取20 m, 不考虑基床表层的粘结强度, 且底座板刚度保持不变, 分析冻胀量对离缝的影响, 冻胀量d分别取5 mm、10 mm、15 mm、20 mm、30 mm、40 mm, 冻胀量d对底座板变形的影响见图8.

图8 冻胀量对底座板变形的影响Fig.8 Influence of frost heaving amounts on bed plate deformation

从整体来看, 轨道结构的跟随性良好, 在不同冻胀量下轨道结构冻胀曲线呈现相同的规律, 只是底座板变形的大小不同, 见图8.现以冻胀量为20 mm为例, 离缝值在不同位置的变化情况见图9.由图9可知, 离缝值主要出现在底座板伸缩缝处和冻胀起始位置, 离缝最大值位于伸缩缝处, 达到2.2 mm.图10所示为最大离缝值随冻胀量的变化情况.

图9 离缝值纵向分布规律Fig.9 Longitudinal distribution law of seam value

图10 冻胀量对最大离缝值的影响Fig.10 Influence of frost heaving amounts on the maximum seam value

由图10可见, 最大离缝值随着冻胀量的增大而不断增大, 基本呈线性变化, 在冻胀量为11.3 mm时, 最大离缝值达到1.0 mm.根据工程实际情况, 认为当离缝值大于1.0 mm时层间粘结力失效, 轨道结构出现脱空现象, 会对行车产生不利影响.可见, 在路基冻胀波长为20 m, 冻胀量达到11.3 mm时, 应进行冻胀防治或治理.

3 基床表层粘结强度的影响

内聚力模型最早是由Xu等提出的一种描述断裂的典型模型[11].内聚力模型在模拟不同层与层之间的粘结失效而产生的病害问题上被广泛使用, 它在一定程度上还可以模拟裂纹的扩展, 但具有局限性, 只能模拟粘结层内部的裂纹扩展.在大多数的研究中[12, 13, 14], 内聚力模型主要是用来模拟胶结层间的脱离情况, 在轨道工程领域, 诸多学者利用内聚力模型对CA砂浆层的粘结进行模拟, 并研究了砂浆层破坏发展机理, 但针对底座板与基床表层之间的粘结模拟较少.

本文采用内聚力模型来模拟基床表层与混凝土底座板之间的粘弹性, 在应用内聚力模型模拟裂纹的产生和扩展时, 需要在预计产生裂纹的区域加入cohesive层.本文在原来模型的基础上设置cohesive层厚度为0, 与底座板与基床表层之间为tie连接, 如图11所示.不同基床表层粘结强度影响下轨道结构变形与离缝发展情况如图12~图15所示.

图11 内聚力模型Fig.11 Cohesive model

图12 粘结层破坏发展历程Fig.12 Developing process of adhesive layers’ damage

图13 离缝横向发展示意图Fig.13 Schematic diagram of the transverse development of seam

图14 粘结强度对离缝值的影响Fig.14 Influence of coherence strengths on seam value

图15 冻胀量对离缝长度的影响Fig.15 Influence of frost heaving amounts on seam length

选取冻胀量20 mm, 冻胀波长20 m, 冻胀作用在伸缩缝处, 底座板刚度保持不变, 研究粘结强度为0.1 MPa时, 底座板与基床表层之间的离缝发展情况.图12为一块轨道板范围内粘结层的发展破坏过程, 冻胀波峰位于图中粘结层底部.

图12中运用损伤因子来反映粘结层的破坏程度, 当损伤因子大于1时说明完全破坏.由图12可知, 在冻胀量d为3.8 mm时, 在冻胀波峰位置处粘结层部分单元开始失效, 离缝开始产生, 之后扩展比较缓慢; 在冻胀量d为3.9 mm时, 冻胀侧边的粘结层开始失效, 并呈放射状向周围方向扩展, 发展比较迅速; 在冻胀量d为4.5 mm时, 轨道结构横向位置形成了贯通离缝.

离缝在轨道结构横向位置从板中向两边发展, 见图13.这是由于轨道结构的组成从整体来看是中间刚度大, 两侧刚度小.在应变一定的条件下, 刚度越大所产生的应力越大, 所以离缝是在板中间产生并向两边逐渐扩展, 直至离缝形成横向贯通.

选择路基冻胀量20 mm, 波长20 m, 冻胀作用在两伸缩缝的中间位置, 且底座板刚度保持不变.研究粘结强度P为0.1 MPa、0.4 MPa、1.0 MPa时离缝在无砟轨道结构纵向的发展情况, 主要分析在半个冻胀波长下离缝值和离缝长度随粘结强度的变化情况, 见图14.由图14知, 随着底座板与基床表层之间粘结强度的增加, 离缝值不断减小, 当粘结强度达到1.0 MPa时, 底座板与基床表层无离缝产生.

图15是粘结强度P为0.1 MPa、0.4 MPa和1.0 Mpa时, 离缝长度在不同冻胀量下的变化情况.由图15可知, 当粘结强度为0.1 MPa时, 冻胀量达到4.5 mm时开始产生离缝, 离缝长度开始增长迅速, 之后逐渐趋于平稳, 最终离缝长度为13.37 m.当粘结强度为0.4 MPa时, 冻胀量达到13.5 mm时开始产生离缝, 离缝发展迅速, 当冻胀量为20 mm时离缝长度达到9.84 m.当粘结强度为1.0 MPa时, 基本没有离缝产生.

4 底座板刚度折减的影响

由于在浇筑过程中养护条件不佳、预应力设置不足、外界自然条件恶劣及列车荷载的冲击作用, 在混凝土底座板内部会产生拉应力, 从而产生细小的内部裂纹, 造成混凝土刚度折减.在严寒地区, 巨大的温度力甚至使底座板开裂, 混凝土刚度折减的问题更为突出.在以往的研究中, 很少考虑混凝土刚度折减, 对计算结果产生很大影响[15].因此, 本文考虑底座板的刚度折减具有实际的工程指导意义.

本文分析底座板刚度折减对轨道结构变形的影响时, 刚度折减系数分别取0.6和0.8.计算中, 冻胀曲线选取波长10 m, 冻胀量20 mm的半余弦曲线, 冻胀作用于底座板两伸缩缝的中间位置, 且不考虑基床表层的粘结强度.底座板刚度对离缝值的影响见图16.

图16 底座板刚度对离缝值的影响Fig.16 Influence of bed plate stiffnesses on seam value

由图16可知, 最大离缝值出现在底座板伸缩缝处, 随着底座板刚度的降低, 最大离缝值呈线性减小的趋势.当不考虑刚度折减时, 即刚度为33 GPa时, 离缝最大值为7.75 mm; 当底座板刚度降为原来的80%时, 离缝最大值为7.02 mm, 减幅为9.4%; 当底座板的刚度降为原来的60%时, 离缝最大值减小到6.22 mm, 减幅达到19.7%.可见, 底座板刚度折减对底座板与基床表层之间的离缝值影响较大.

5 结论

1)冻胀由路基传到轨面, 形成的轨道不平顺与冻胀波形基本一致.冻胀作用位置对轨道结构变形影响较大, 当冻胀作用于轨道板中间位置时, 离缝值最大.

2)随着冻胀波长的增加, 最大离缝值逐渐减小; 当冻胀量为20 mm, 冻胀波长大于30 m时, 最大离缝值小于1.0 mm, 可认为不产生脱空区域.

3)随着冻胀量的增加, 最大离缝值呈线性增大; 当冻胀波长为20 m, 冻胀量大于11.3 mm时, 出现了大于1.0 mm的最大离缝值, 应及时进行冻胀防治或治理.

4)底座板与基床表层之间的离缝在横向位置是从中间向两边逐渐扩展的, 而且发展速度较快, 很快形成横向贯通, 纵向位置离缝长度增长较快.

5)底座板刚度折减对最大离缝值有较大的影响, 当刚度减小为原来的60%时, 离缝值减小了近20%, 说明在冻胀地区可适当减小底座板刚度来控制离缝值.

The authors have declared that no competing interests exist.

参考文献
[1] 赵国堂. 严寒地区高速铁路无砟轨道路基冻胀管理标准的研究[J]. 铁道学报, 2016, 38(3): 1-8.
ZHAO Guotang. Study on management stand ard of frost heaving of ballastless track subgrade on high-speed railway in severe cold regions[J]. Journal of the China Railway Society, 2016, 38(3): 1-8. (in Chinese) [本文引用:4]
[2] MA Fuxun, XI Ruijie, XU Nan. Analysis of railway subgrade frost heave deformation based on GPS[J]. Geodesy and Geodynamics, 2016, 7(2): 143-147. [本文引用:2]
[3] 赵国堂, 刘秀波, 高亮, . 哈大高速铁路路基冻胀区轨道不平顺特征分析[J]. 铁道学报, 2016, 38(7): 105-109.
ZHAO Guotang, LIU Xiubo, GAO Liang, et al. Characteristic analysis of track irregularity in subgrade frost heave area of Harbin-Dalian high-speed railway[J]. Journal of the China Railway Society, 2016, 38(7): 105-109. (in Chinese) [本文引用:1]
[4] 杨俊斌, 刘学毅, 刘永孝, . Ⅰ型轨道板端离缝对轨道结构及车辆动力特性的影响[J]. 西南交通大学学报, 2014, 49(3): 432-437.
YANG Junbin, LIU Xueyi, LIU Yongxiao, et al. Influence of seam in type Ⅰ slab ends on dynamic characteristics of track and vehicle[J]. Journal of Southwest Jiaotong University, 2014, 49(3): 432-437. (in Chinese) [本文引用:1]
[5] 邵丕彦, 李海燕, 吴韶亮, . CRTS I型轨道板温度变形及与砂浆垫层间离缝的测试研究[J]. 中国铁道科学, 2013, 34(2): 18-22.
SHAO Piyan, LI Haiyan, WU Shaoliang, et al. Measurement and research on temperature warping of CRTS I track slab and crack between track slab and cement asphalt mortar cushion[J]. China Railway Science, 2013, 34(2): 18-22. (in Chinese) [本文引用:1]
[6] ZHU Shengyang, CAI Chengbiao. Fatigue life prediction of CRTS I ballastless slab track[C] //Third International Conference on Transportation Engineering. Chengdu, China, 2011: 1714-1719. [本文引用:1]
[7] ZENG X H, XIE Y J, DENG D H, et al. A study of the dynamic mechanical properties of CRTS I type CA mortar[J]. Construction and Building Material, 2016, 112(7): 93-99. [本文引用:1]
[8] KÜBLER L, EBERHARD P, GEISLER J. Flexible multibody systems with large deformations and nonlinear structural damping using absolute nodal coordinates[J]. Nonlinear Dynamics, 2003, 34(1): 31-52. [本文引用:1]
[9] SCHENK O, GÄRTNER K, FICHTNER W, et al. PARDISO: a high-performance serial and parallel sparse linear solver in semiconductor device simulation[J]. Future Generation Computer Systems, 2001, 18(1): 69-78. [本文引用:1]
[10] GONZÁLEZ F, LUACES A, LUGRIS U, et al. Non-intrusive parallelization of multibody system dynamic simulations[J]. Computational Mechanics, 2009, 44(4): 493-504. [本文引用:1]
[11] XU X P, NEEDLEMAN A. Numerical simulations of fast crack growth in brittle solids[J]. Journal of the Mechanics and Physics of Solids, 1994, 42(9): 1397-1434. [本文引用:1]
[12] SILVA C M, FLORES-COLEN I, GASPAR S. Numerical analysis of renders' adhesion using an interface model[J]. Construction & Building Materials, 2013, 38(1): 292-305. [本文引用:1]
[13] EVANGELISTA J F, ROESLER J R, PROENCA S P. Three-dimensional cohesive zone model for fracture of cementitious materials based on the thermodynamics of irreversible processes[J]. Engineering Fracture Mechanics, 2013, 97(1): 261-280. [本文引用:1]
[14] GHOVANLOU M K, JAHED H, KHAJEPOUR A. Cohesive zone modeling of fatigue crack growth in brazed joints[J]. Engineering Fracture Mechanics, 2014, 120(4): 43-59. [本文引用:1]
[15] 李思云, 杨荣山. CRTS-I型板式无砟轨道疲劳寿命研究[J]. 铁道标准设计, 2016, 60(3): 34-37.
LI Siyun, YANG Rongshan. Research on fatigue life of CRTS-I slab ballastless track[J]. Railway Stand ard Design, 2016, 60(3): 34-37. (in Chinese) [本文引用:1]