黄土地区某铁路专用线路基动力响应规律
Dynamic Responses Law of Subgrade of a Special Railway Line in Loess Region
通讯作者: 胡志平,男,教授,博士生导师,电话(Tel.):029-82337356;E-mail:huzhping@chd.edu.cn.
责任编辑: 陈晓燕
收稿日期: 2021-04-22
基金资助: |
|
Received: 2021-04-22
作者简介 About authors
王瑞(1992-),男,河北省邯郸市人,讲师,现主要从事路基工程相关研究.
为了研究移动列车荷载下黄土地区某铁路专用线动力响应规律,结合现场测试数据及2.5D有限元计算结果评估了新建货运线路和临近既有线列车运行情况下路基动应力分布规律及影响深度,讨论了沿线办公区环境振动规律.结果表明,轨道系统的过度抽象会导致路基内部动应力幅值和影响范围偏大,道砟-枕木-紧固件系统的实体化建模使得路基内部动应力分布更为合理.路基处于拟静力状态时动力蠕变引发的能量消耗对路基动力响应的影响可以忽略.货运列车低速运行时,车速对路堤内部动应力分布影响微弱,动应力影响深度约为4.2 m.此外,车速越低,轨道不平顺对场地环境振动的影响越明显.改善轨道短波不平顺是控制环境振动的有效措施.
关键词:
In order to investigate the dynamic response law of subgrade of a special railway line in loess region under moving train loading, the field test data and 2.5D finite element results were adopted to assess the dynamic stress distribution and its influence depth, and the environmental vibration law in the office area along the line was also analyzed. The results show that the excessive abstraction of the track system might enlarge the amplitude and influence the range of dynamic stress. The entity modeling of the ‘ballast-sleeper-fastener’ system makes the distribution of the dynamic stress in the subgrade more reasonable. The influence of energy consumption caused by the dynamic creep of soil have a negligible influence when the subgrade is in a quasi-static state. The speed has a negligible effect on the distribution of the dynamic stress of the subgrade when a freight train is running at a low speed, and the influence depth of the dynamic stress is about 4.2 m. As the speed decreases, the influence of track irregularity on the environmental vibration increases. The environmental vibration can be effectively controlled by improving the short wave length irregularity.
Keywords:
本文引用格式
王瑞, 胡志平, 殷珂, 马甲宽, 任翔.
WANG Rui, HU Zhiping, YIN Ke, MA Jiakuan, REN Xiang.
交通强国战略的提出再次肯定了交通运输在国民经济发展中的重要驱动作用.铁路运输以其运量大、速度较快等特点在公共交通领域备受重视.长期以来,从业人员一直致力于改善路基服役性能,提升铁路运输的社会经济效益.列车运行荷载下路基动力响应研究涉及附加动应力分布及地表环境振动两方面主要内容,是路基工程设计、施工和减隔振的重要参考.
目前主流的两种动应力影响深度界定标准分别是“动应力幅值衰减至最大幅值的10%时对应的深度”和“动应力幅值沿深度分布曲线与1/10或1/5自重应力线交点的深度”[7].由于围压是影响路基土体力学和强度参数的重要因素,一般推荐采用动应力幅值曲线与自重应力线交点确定动应力影响深度.总体而言,附加动应力分布受路基土质条件、道床结构、车速和轴重等多重因素影响,衰减曲线大同小异,动应力影响深度普遍在3~6 m范围内.
目前的大量研究习惯针对具体工程背景评估列车运行诱发的环境振动对周边敏感建筑或区域的影响.孟庆成等[8]基于现场测试和数值模拟评估了某火车站高架候车厅结构车致振动特性及振动对舒适度的影响.Yang等[9]采用2.5D有限元方法研究了高速列车运行情况下空沟的隔振效果.Thompson等[10]通过2.5D有限元方法研究了路基内部横向波阻块的减振效果及机制.Coulier等[11]以西班牙El Realengo地区某铁路沿线场地为背景研究了混凝土灌注桩挡墙对地表环境振动的隔振效果.此外,Dijckmans等[12]以瑞典Furet地区某铁路沿线场地为背景研究了波纹板桩板墙的隔振效果.Yarmohammadi等[13]结合列车荷载激发环境振动的特点发展了一种耦合方法,研究指出双隔振沟的减振效果比单沟更好.马骙骙等[14]通过开展宝兰高铁路堑段地面垂向振动现场测试研究了高速列车荷载下路堑的地面垂向振动随距离传播规律,指出场地速度特性(地基覆盖层与下卧层模量比、覆盖层厚度)是影响地面振动剧烈程度的重要因素.相比于动应力分布规律,环境振动研究涉及振动理论知识、情况各异的路基断面甚至精细化的地基-建筑物基础相互作用等问题,因此获得的量化成果相对较少.
在高速铁路大量修建的政策背景下,铁路路基动力响应研究大多集中在我国东南沿海软土地区,黄土地区相关研究的系统性和深度均有所欠缺.由于黄土地区沟壑纵横,修筑铁路线路时大量采用填土路堤结构.阐述列车运行荷载引发的路堤-地基系统响应规律是细化设计参数、评估服役性能的重要举措.
本文依托某铁路专用线工程,旨在通过现场测试和2.5D有限元模拟揭示新建货运线路和临近既有客运线列车运行情况下路基动力响应规律,评估货运列车运行情况下路基动应力分布及影响深度,讨论线路附近办公区域的环境振动规律.以2.5D有限元为主要研究手段,辅以现场实测数据,在验证数值模拟精度的前提下开展参数分析.
1 现场测试概况
依托陕北黄土梁峁沟壑区某铁路专用线工程开展,测试断面位于洛河一级阶地,原始地基未经强夯处理.地层由上向下依次为黏质黄土(Ⅰ级非自重湿陷性黄土)及下覆基岩(页岩夹砂岩).黄土覆盖层厚度约10.0 m,地层分布较均匀,起伏变化较小.路堤中心填高10.8 m,边坡坡率1∶1.5.地下水位标高在路堤底面下部5 m处,年变化幅度在2 m左右.路堤主体采用压实黄土填筑,与既有线采用台阶搭接.
根据路堤填筑高程分层多次埋设土压力传感器,待路堤填筑完成且变形稳定后开展现场测试工作.路基断面轮廓、传感器分布及现场测试情况如图1所示.根据动应力分布规律及影响深度,在路堤中线处埋设4层土压力盒,最大埋深8 m.
图1
测试采用BY-1型电阻式双油腔动土压力盒,盘面尺寸20 cm,量程0.4 MPa.土压力盒盘面较薄,遇尖锐物体易损坏,因此并未在道床底面及基床表层范围内布设.数采仪器型号为DH5922N,采样频率 1000 Hz.测试重载列车由1节车头及10节满载棚车组成,如图2所示.车头为东风4B型货运内燃机车,轴重23 t;棚车为P64K型,满载轴重为 22.5 t,测试车速30 km/h.
图2
图2
车型及轮轴荷载示意(m)
Fig.2
Vehicle type and geometric profile of train wheel loads (m)
2 2.5D有限元模型验证及优化
2.1 2.5D有限元简介
2.5 D有限元计算方法是基于路基断面(位于x-y平面)在列车运行方向z保持不变的基本假定建立的.计算时首先采用双重傅里叶变换得到断面各节点在频率-波数域内的响应,然后由傅里叶反变换获得系统的实际响应.沿列车运行方向z和时间t的傅里叶变换可以表示为
式中:u为位移;ξz为波数;ω为圆频率;i为虚数单位;上标“-”和“~”分别代表波数和频率变换.
2.1.1 列车荷载
三维空间中沿z轴移动的荷载可以表示为
式中:函数ψ(x, y)表示荷载在与z轴垂直断面上的分布特征;函数ϕ(z-vt)表示荷载的移动特征,v为荷载的移动速度;函数q(t)表示荷载的时变特征.本文列车荷载在变换域中的表达式参照文献[20].
2.1.2 地基系统动力方程
系统在变换域中的动力方程可表达为
式中:
2.1.3 轨道系统动力方程
将轨道系统模拟为欧拉梁[20],轮轴荷载作用下变换域内轨道的振动变形可以描述为
式中:K为欧拉梁的弯曲刚度,E为弹性模量,I为截面惯性矩;m为欧拉梁的单位长度质量;
2.2 程序验证
图3
图3
X2000列车运行下路基动力响应数值重现
Fig.3
Numerical recurrence of subgrade dynamic response when Train X2000 is running
2.3 现场测试规律的数值重现
图4
表1 模型土体参数
Tab.1
土层 | 厚度/ m | 密度/ (kg·m-3) | 模量/ MPa | 阻尼比 | 泊松比 |
---|---|---|---|---|---|
基床表层(新建) | 0.7 | 1900 | 86.5 | 0.04 | 0.25 |
基床底层(新建) | 2.3 | 1850 | 75.0 | 0.04 | 0.25 |
路堤本体(新建) | 1.75 | 1750 | 50.0 | 0.04 | 0.25 |
场地填方体 | 7.0 | 1750 | 50.0 | 0.04 | 0.25 |
基床表层(既有) | 0.7 | 1900 | 86.5 | 0.04 | 0.25 |
基床底层(既有) | 2.3 | 1850 | 75.0 | 0.04 | 0.25 |
路堤本体(既有) | 4.7 | 1750 | 50.0 | 0.04 | 0.25 |
地基土层 | 10 | 1680 | 26.3 | 0.05 | 0.25 |
基岩 | 10.0 | 2000 | 126.0 | 0.02 | 0.32 |
图5
图5
现场测试数据与有限元结果对比
Fig.5
Comparison of field test data and finite element results
即便如此,对比响应幅值随深度的变化规律可以看出,有限元计算获得的幅值衰减速度比现场测试数据慢.在2.5D数值模拟方面,目前针对其精度影响因素的争论主要集中在轨道系统抽象的合理性、路基土体的非线性及轮-轨相互作用等方面.在重现路基实际响应规律时,土体动力参数的代表性也会影响预测效果.相关研究指出,交通荷载作用下土体动力蠕变会带来额外的能量消耗[22].本文有意依据本次现场实测数据开展进一步研究,分析轨道系统建模方式和土体动力蠕变对路基动力响应的影响.相关研究指出,交通荷载下压实黄土的动应力-应变关系在一定起始应力范围内的非线性表现不明显[22].结合列车运行情况下路基内部动应力幅值分布情况,暂不考虑土体非线性的影响.
2.4 轨道系统各向异性实体建模对模拟结果的影响
2.4.1 道砟-枕木-紧固件各向异性单元建模信息
图6
图6
道砟-枕木-紧固件系统实体单元示意(m)
Fig.6
Entity finite element of ballast-sleeper-fastener system (m)
参照已有文献,将轨道等效为欧拉梁,单根轨道抗弯刚度K=6.4 MN·m2,单位长度质量m=60 kg/m.紧固件单元只考虑其竖向等效模量,其他方向的模量和单元质量设为一极小值.等效模量计算式为
式中:krp为紧固件竖向刚度(取值为1.75×
枕木单元的密度
式中:Esl为枕木混凝土弹性模量(为 12500 MPa);bsl为枕木平均宽度(为0.25 m);ρsl为材料密度(为 2500 kg/m3);υsl为材料泊松比(为0.2)[25].
文献[26]基于大型三轴试验设备开展了道砟材料各向异性参数的试验.依据循环加载试验结果,道砟材料的9个各向异性参数取值分别为弹性模量Ey,b=900 MPa,Ex,b=Ez,b=600 MPa;泊松比 υxy,b=0.31,υxz,b=0.32,υyz,b=0.28;切变模量Gxy,b=(Ex,b+Ey,b)/(2+2υxy,b)=572 MPa,Gxz,b=(Ex,b+Ez,b)/(2+2υxz,b)=455 MPa,Gyz,b=(Ey,b+Ez,b)/(2+2υyz,b)=586 MPa.密度ρb=1900 kg/m3,阻尼比η=0.2.
2.4.2 路基响应规律
图7
图7
不同轨道系统建模情况下动应力
Fig.7
Dynamic stress in different track system modeling strategies
图8
图8
不同轨道系统建模情况下动应力幅值横向分布
Fig.8
Lateral distributions of dynamic stress amplitude in different track system modeling strategies
图9
图9
不同轨道系统建模情况下路堤中心线动应力衰减规律
Fig.9
Dynamic stress attenuation along the centerline of embankment in different track system modeling strategies
由图8可知,由于考虑了枕木的横向刚度,轨道系统实体化后浅层动应力幅值横向分布呈马鞍形,更符合路基动应力分布情况.采用欧拉梁模拟轨道系统时基床表层顶面动应力均匀分布且幅值及范围均偏大.由此导致动应力向下传递时影响范围和幅值偏大.可见在考虑了真实的荷载传递路径之后,道砟-枕木-紧固件的实体化建模使得路基内部动应力分布更为合理.轨道系统的过度抽象会导致路基内部动应力幅值和影响范围偏大.
对比图9所示路堤中心线处动应力衰减规律可以看出,轨道系统实体化建模后路基内部动应力幅值与现场实测数据更加吻合.另外,本文实体化建模后单元数仅增加44个(共计 1672 个),可不考虑其对计算效率的影响.
2.5 土体动力蠕变能耗对路基动力响应的影响
图10
动应力在路基中迅速衰减,因此对不同深度的土层采用差异化的附加阻尼比系数,为了充分体现动力蠕变的影响,设置基床表层、基床底层和路堤本体的附加系数分别为1.0、0.5及0.2.提取附加阻尼前后路堤内部动应力衰减规律如图11所示.
图11
图11
附加阻尼前后路堤中心线动应力衰减规律
Fig.11
Dynamic stress attenuation along the centerline of embankment before and after appending damping
低车速情况下附加阻尼对路基内部动力衰减无影响;高车速时阻尼比增加会轻微减小路基内部动应力幅值.已知40 m/s的车速远小于基床土体瑞雷波速(约130 m/s),此时路基属于拟静力响应范畴,土体中并无波动传播,可能导致阻尼比对系统响应影响微弱.
当荷载移动速度大于地表瑞雷波速时会产生波动传播现象,路基由拟静力状态转换为动力响应状态.阻尼比的增大会强化路基材料对波动能量的吸收效果,从而减小动应力幅值.本文所涉工程背景不存在高速运行工况,因此在后续分析中不考虑附加阻尼的影响.
3 填方路基动力响应规律
基于上述分析,本文采用各向异性实体单元模拟轨道系统,在不考虑动力蠕变能耗的情况下讨论新建铁路和既有线列车运行时图4中的路基动力响应规律.左侧填方场地布置有设备间及办公、休息区域,因此重点关注新建线路路堤内部动应力分布及填方场地地表环境振动规律.
新建线路为货运线路,列车运行速度偏低,因此采用图2所示车辆模型,讨论车速较低情况下(10~40 m/s),路基内部动应力分布及地表环境振动规律.由于缺乏车体动力学参数,暂不考虑轨道不平顺影响.既有线为客货共线铁路,采用常见的CRH3型车体参数,讨论列车在较高速度运行情况下引发的地表环境振动规律,同时评估轨道不平顺的影响.
3.1 新建线路列车运行情况下路基响应规律
提取不同车速下路基内部动应力分布规律,如图12所示.
图12
图12
新建线路列车运行情况下路基动应力衰减规律
Fig.12
Dynamic stress attenuation of subgrade under a moving train on a new line
可以看出,货运列车低速运行情况下,车速对路堤内部动应力分布的影响十分微弱.定义1/10自重应力线与动应力幅值分布曲线的交点所在深度为列车荷载的影响深度.参照《铁路路基设计规范》TB 10001—2016[5]将道床等上部结构等效为14 kPa的均匀荷载,路基土重度统一取18.5 kN/m3.本文工程背景下23 t轴重货运列车荷载对路基的影响深度约为4.2 m.
采用1/3倍频程中心频率f1/3的分频z振级VLz表征列车荷载引发的环境振动,其表达式如下:
选取左侧地表与新建路堤中线水平距离约10、16、24及35 m的4个测点,获得40 m/s运行速度情况下各测点加速度频谱及分频z振级随1/3倍频中心频率变化,如图13所示.图中:aamp为加速度幅值;f为频率.可以看出,测点加速度响应集中在 0~4 Hz低频范围内.各测点分频z振级随中心频率变化规律类似.
图13
图13
测点加速度频谱曲线及分频z振级随1/3倍频中心频率变化规律(v=40 m/s)
Fig.13
Acceleration spectrum curve and z vibration law of variation at 1/3 octave band center frequency (v=40 m/s)
为了分析列车运行情况下场地环境振动沿地表横向衰减规律,进一步提取地表节点最大分频z振级VLz,max分布规律,如图14所示,图中Hsur为与线路垂直的水平横向地表位置.随着列车运行速度的提高,地表环境振动强度逐渐增大但衰减规律类似.行业标准《城市轨道交通引起建筑物振动与二次辐射噪声限值及其测量方法标准》JGJ/T 170—2009[27]规定4类区域(交通干线两侧区域)建筑物室内振动限值为昼间75 dB、夜间72 dB.当前工况下距离路基中心线15 m以外设备间及办公区域的最大分频振级基本不超过72 dB.货运线路车速一般较低,以往研究重点关注轴重对路基响应的影响,以至于本文所涉货运内燃机车及棚车车体动力学参数十分鲜见,因此本文暂不考虑新建线路轨道不平顺对环境振动的影响.
图14
图14
路基表面最大分频z振级分布规律(新建线路加载)
Fig.14
Maximum z vibration law on surface of subgrade (new line loading)
3.2 既有线列车运行引发的地表环境振动规律
图15
图15
路基表面最大分频z振级分布规律(既有线加载)
Fig.15
Maximum z vibration law on surface of subgrade (current line loading)
列车在竖向不平顺的轨道上运行时,车体会发生振动从而导致轮轨接触力的改变.一般车速越高不平顺对轮轨力的影响越大.在研究客运列车高速运行诱发的铁路路基动力响应规律时不宜忽略轨道不平顺的影响.本文参考德国高速铁路高低不平顺谱,根据其界定的波长-波高对应关系设定3种不平顺参数(参数I~III 的波长分别为2、10、50 m,波高分别为0.15、0.6、2.4 mm)初步评估轨道不平顺对环境振动的影响.10自由度车体参数及轨道不平顺在2.5D有限元中的实现详见文献[28].
考虑轨道不平顺时,既有线列车运行情况下路基环境振动规律如图16所示.对比平顺情况下的响应可以看出,当客运列车以常规车速(40 m/s)行驶时,不平顺波长越大,轨道不平顺对环境振动的影响越小.在空间上,随着与加载点距离的增大,不平顺对环境振动的影响逐渐显现,并且存在分叉现象.波长越大,最大分频z振级曲线的分叉点距离加载点越远.其原因可能是不平顺荷载会激发地基内部波动传播,使得不平顺荷载产生的附加响应在路基内的衰减规律与平顺荷载产生的响应衰减规律产生差异,加之式(3)对数取值会放大远场微小的响应差异,由此产生了明显的分叉现象.上述现象表明改善轨道短波不平顺可以有效控制列车运行诱发的环境振动.
图16
图16
考虑轨道不平顺时既有线列车运行情况下路基环境振动规律
Fig.16
Environmental vibration of subgrade under current line loading considering track irregularity
观察车速较高时(80 m/s)的环境振动规律可以看出,轨道不平顺对环境振动的影响减弱.由于本文既有线为非高速铁路线路,所以对此工况不做过多讨论.
4 结论
本文以黄土地区某铁路专用线工程为背景,通过现场测试得到了填方路堤内部的动应力响应规律,据此验证了2.5D有限元数值模拟精度并探讨了轨道系统各向异性实体化建模及土体动力蠕变对模型动力响应的影响.在此基础上,分别讨论了新建线路运行货运列车及既有线运行客运列车情况下路基动力响应规律,得出以下几点结论.
(1) 轨道系统的过度抽象会导致路基内部动应力幅值和影响范围偏大.由于考虑了真实的荷载传递路径,道砟-枕木-紧固件系统的各向异性实体化建模使得路基内部动应力分布更为合理.
(2) 路基处于拟静力响应状态时动力蠕变引发的能量消耗对路基动应力衰减规律的影响可以忽略,路基处于动力响应状态时动力蠕变引发的能量消耗会减缓路基内部动应力的衰减.
(3) 23 t轴重货运列车低速运行时,车速对黄土填方路堤内部动应力分布影响微弱,动应力影响深度约为4.2 m.
(4) 随着地表测点与加载位置间距的增大,轨道不平顺对场地环境振动的影响逐渐显现.车速越低轨道不平顺对场地环境振动的影响越明显.改善轨道短波不平顺是控制环境振动的有效措施.
参考文献
无碴轨道路基基床动力特性的研究
[J]. ,
Study of dynamic characteristics of soil subgrade bed for ballastless track
[J]. ,
高铁膨胀土新型路堑基床动力特性与参数敏感性
[J]. ,
Dynamic characteristics and parameter sensitivities of new cutting subgrade for high-speed railway in expansive soil area
[J]. ,
高速铁路路基基床结构分析及设计方法
[J]. ,
Structural analysis and design method for subgrade bed of high-speed railway
[J]. ,
重载铁路路基基床结构分析及设计方法
[J]. ,
Structural analysis and design method for subgrade bed of heavy haul railway
[J]. ,
A new method for modelling and simulation of the dynamic behaviour of the wheel-rail contact
[J]. ,DOI:10.1007/s11633-012-0640-6 URL [本文引用: 1]
高架候车厅车致振动特性及减振控制研究
[J]. ,
Research of vibration characteristics and vibration control of elevated waiting halls
[J]. ,
2.5D vibration of railway-side buildings mitigated by open or infilled trenches considering rail irregularity
[J]. ,DOI:10.1016/j.soildyn.2017.12.027 URL [本文引用: 2]
Mitigation of railway-induced vibration by using subgrade stiffening
[J]. ,DOI:10.1016/j.soildyn.2015.09.005 URL [本文引用: 2]
Experimental and numerical evaluation of the effectiveness of a stiff wave barrier in the soil
[J]. ,DOI:10.1016/j.soildyn.2015.04.007 URL [本文引用: 1]
Efficacy of a sheet pile wall as a wave barrier for railway induced ground vibration
[J]. ,DOI:10.1016/j.soildyn.2016.02.001 URL [本文引用: 1]
Design of wave barriers for mitigation of train-induced vibrations using a coupled genetic-algorithm/finite-element methodology
[J]. ,DOI:10.1016/j.soildyn.2019.03.007 URL [本文引用: 1]
高速铁路路堑段地面振动试验研究及数值分析
[J]. ,
Experimental study and numerical analysis of ground vibration in cutting section of high-speed railway
[J]. ,
A 2.5D finite/infinite element approach for modelling visco-elastic bodies subjected to moving loads
[J]. ,DOI:10.1002/nme.208 URL [本文引用: 1]
Simulation of track-ground vibrations due to a high-speed train: The case of X-2000 at Ledsgard
[J]. ,DOI:10.1016/S0022-460X(02)01007-6 URL [本文引用: 2]
Critical speed of railway tracks. Detailed and simplified approaches
[J]. ,DOI:10.1016/j.trgeo.2014.09.003 URL [本文引用: 1]
Investigating ground vibration induced by moving train loads on unsaturated ground using 2.5D FEM
[J]. ,DOI:10.1016/j.soildyn.2019.05.034 URL [本文引用: 2]
Pore pressure generation in a poro-elastic soil under moving train loads
[J]. ,DOI:10.1016/j.soildyn.2019.105711 URL [本文引用: 1]
基于2.5维有限元方法模拟高速列车产生的地基振动
[J]. ,
Numerical simulation of high-speed train induced ground vibrations using 2.5D finite element approach
[J]. ,
2.5D有限元建模关键问题: 边界条件、网格划分及计算域选取
[J]. ,
Key issues in modeling process of 2.5D finite element method—Boundary conditions, meshing and computing range selection
[J]. ,
交通荷载引起的静偏应力对压实黄土动力特性的影响
[J]. ,
Study on dynamic characteristics of compacted loess subjected to static deviatoric stress induced by traffic loading
[J]. ,
Long-term dynamic stability of improved loess subgrade for high-speed railways
[J]. ,DOI:10.1680/jgeen.19.00088 URL [本文引用: 1]
2.5维有限元分析饱和地基列车运行引起的地面振动
[J]. ,
Ground vibration induced by trains moving on saturated ground using 2.5D FEM
[J]. ,
Using the 2.5D FE and transfer matrix methods to study ground vibration generated by two identical trains passing each other
[J]. ,DOI:10.1016/j.soildyn.2018.06.025 URL [本文引用: 2]
Railway ballast anisotropy testing via true triaxial apparatus
[J]. ,DOI:10.1016/j.trgeo.2020.100355 URL [本文引用: 1]
/
〈 | 〉 |