冰川冻土, 2022, 44(2): 555-565 doi: 10.7522/j.issn.1000-0240.2022.0059

寒区工程与灾害

列车荷载下冻土路基非平稳随机振动分析的时域显式法

陈智锋,1, 李颖雄2, 赖远明,3, 苏成1,4

1.华南理工大学 土木与交通学院,广东 广州 510640

2.西南交通大学 智慧城市与交通学院,四川 成都 611756

3.中国科学院 西北生态环境资源研究院 冻土工程国家重点实验室,甘肃 兰州 730000

4.华南理工大学 亚热带建筑科学国家重点实验室,广东 广州 510640

Non-stationary random vibration analysis of railway embankments in permafrost regions under train loads using the explicit time-domain method

CHEN Zhifeng,1, LI Yingxiong2, LAI Yuanming,3, SU Cheng1,4

1.School of Civil Engineering and Transportation,South China University of Technology,Guangzhou 510640,China

2.Institute of Smart City and Intelligent Transportation,Southwest Jiaotong University,Chengdu 611756,China

3.State Key Laboratory of Frozen Soil Engineering,Northwest Institute of Eco-Environment and Resources,Chinese Academy of Sciences,Lanzhou 730000,China

4.State Key Laboratory of Subtropical Building Science,South China University of Technology,Guangzhou 510640,China

通讯作者: 赖远明,院士,研究员,主要从事寒旱区工程理论和数值分析研究. E-mail: ymlai@lzb.ac.cn

收稿日期: 2021-12-30   修回日期: 2022-04-10  

基金资助: 国家自然科学基金项目.  51678252.  52178479

Received: 2021-12-30   Revised: 2022-04-10  

作者简介 About authors

陈智锋,硕士研究生,主要从事冻土路基随机振动研究.E-mail:317565114@qq.com , E-mail:317565114@qq.com

摘要

随着我国冻土地区铁路运营里程的不断提升,冻土路基在随机列车荷载作用下的动力响应分析成为了急需解决的工程问题。本文以青藏铁路某路基横断面为例,采用时域显式蒙特卡罗模拟法,计算其在随机列车荷载作用下响应的统计特征。首先,提出列车驶过施加在路基顶面荷载的快速计算方法,并引入轨道不平顺及列车行驶速度两类随机参数,以生成随机分析所需的大量荷载样本。然后,采用数值积分方法将运动方程在时域上进行离散,建立任意离散时刻冻土路基动力响应关于列车荷载的显式表达式。基于该显式表达式,可以高效地实施蒙特卡罗模拟,得到在随机列车荷载作用下路基关键响应的均值、标准差和峰值等统计量。采用该方法,分析了路基在夏冬两个季节不同深度处的随机动力响应,发现在夏季路基浅层中的位移和速度响应最为剧烈。数值算例表明,时域显式蒙特卡罗模拟法在冻土路基随机振动分析中具有理想的计算精度和计算效率。

关键词: 冻土路基 ; 随机振动 ; 时域显式法 ; 蒙特卡罗模拟

Abstract

With the rapidly increasing railway operating length in the permafrost regions of China, dynamic analysis of permafrost embankment under random train loads has become an urgent engineering problem to be solved. Taking a cross-section of the embankment of Qinghai-Tibet Railway as an example, this paper adopts the explicit time-domain method to analyze the statistical characteristics of its responses to random train loads. In order to obtain the train load samples, which serve as the inputs in the stochastic analysis, an efficient procedure is proposed incorporating the random parameters of track irregularity and train speed. The numerical integration method is used to discretize the equation of motion in the time domain, and an explicit expression of the dynamic responses of permafrost embankment is established. Based on the established explicit expressions, Monte Carlo simulation can be efficiently conducted, delivering the statistics of the random system outputs, including the mean values, standard deviations and mean peak values of the embankment dynamic responses. Using this method, the random dynamic responses of the permafrost embankment in summer and winter are analyzed. It is found that the displacement and velocity responses in summer are more severe than those in winter. Numerical results indicate that the explicit time-domain method has high accuracy and efficiency in analyzing the random vibration of permafrost embankment.

Keywords: permafrost embankment ; random vibration ; explicit time-domain method ; Monte Carlo simulation

PDF (4569KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈智锋, 李颖雄, 赖远明, 苏成. 列车荷载下冻土路基非平稳随机振动分析的时域显式法[J]. 冰川冻土, 2022, 44(2): 555-565 doi:10.7522/j.issn.1000-0240.2022.0059

CHEN Zhifeng, LI Yingxiong, LAI Yuanming, SU Cheng. Non-stationary random vibration analysis of railway embankments in permafrost regions under train loads using the explicit time-domain method[J]. Journal of Glaciology and Geocryology, 2022, 44(2): 555-565 doi:10.7522/j.issn.1000-0240.2022.0059

0 引言

我国为世界第三冻土大国,多年冻土面积占国土总面积的21.5%,其中高原冻土面积更是居于世界之最1。随着“一带一路”倡议深入推进实施,冻土区的铁路建设与维护得到高度重视。同时,国民经济的发展使得列车行驶频次有所提升,加之受全球变暖2等气候因素的影响,冻土区铁路病害问题越来越突出。为保障冻土区铁路的安全运营,亟须发展高效准确的计算方法,分析列车行驶对冻土路基动力响应的影响。

冻土路基在列车荷载作用下的动力响应是当前的研究热点之一。马巍等3和朱占元等4开展了列车荷载作用下冻土路基动力响应的现场监测研究,为列车通过冻土路段时路基稳定性的预测分析提供了实测资料依据。Stevens5、孔祥兵等6和李双洋等7-9采用数值分析方法,研究了一次行车荷载作用下路基的位移和应力等动力响应问题。此外,还有部分学者10-13进行了列车荷载作用下冻土路基的稳定性、蠕变、轨枕响应和车速对路基响应影响等方面的数值模拟研究。

上述研究中,都采用了确定性的列车荷载,选取某一次列车驶过产生的荷载样本计算路基响应,其结果难以反映铁路实际运行中由轨道不平顺及列车类型和轴重差异导致的荷载随机性14。对于路基的特定横断面,列车驶入和驶出的过程导致其顶部荷载的统计特征具有显著的时变特性。因此,在研究路基动力响应时应将其顶部列车荷载按非平稳随机荷载考虑715。由于冻土独特的物理、热力学性质以及非平稳随机振动分析的复杂性,目前关于随机列车荷载作用下冻土路基动力响应的研究较少。

目前,解决结构非平稳随机振动问题的主要方法有蒙特卡罗模拟法、虚拟激励法16、概率密度演化法17以及时域显式法18等。其中,时域显式法通过构建系统动力响应的时域显式表达式,在随机分析中可以针对任意关键响应进行降维,具有理想的计算效率,已成功应用于高层建筑19、大跨度桥梁20和车桥耦合系统21等大型复杂结构系统的随机振动分析。在随机振动分析中,若已知随机激励的均值和互相关函数等统计特征,则可按统计矩的运算规律根据时域显式表达式直接求得响应的前两阶统计矩,这一方法称为时域显式直接法。若通过数值模拟手段生成大量随机激励样本,则可结合时域显式表达式进行高效蒙特卡罗模拟,这一方法称为时域显式蒙特卡罗模拟法。

为了解决随机列车荷载作用下冻土路基动力响应分析这一工程问题,本文首先提出一种随机列车荷载的简化模拟方法,以获取随机分析所需的荷载样本。然后,通过有限元脉冲响应分析构建路基结构动力响应的时域显式表达式。在此基础上,进一步利用时域显式蒙特卡罗模拟法,对随机列车荷载作用下冻土路基的关键响应开展统计分析,可大幅提高冻土路基结构随机振动的计算效率。最后,通过数值算例展示时域显式蒙特卡罗模拟法在冻土路基非平稳随机振动分析中的计算精度与效率。

1 随机列车荷载模拟方法

在冻土路基随机振动问题中,需要列车驶过引起的随机荷载的统计特性或大量样本作为输入。本文采用常见的二维路基模型7进行计算,此类模型需要列车行驶引起的随机轨枕竖向作用力作为荷载输入15。目前,相关实测数据较少,难以满足冻土路基随机振动分析的需求。为了获取分析所需输入,可以对确定性列车荷载的数值模拟方法22-24加以改进,在其中引入随机参数,进而生成大量随机列车荷载样本。

列车在轨道上行驶的简图如图1所示,其中列车行驶速度为v。为了描述列车在沿轨道方向上的位置,建立如图1所示的直线坐标系Ox。列车前端的初始坐标假定为-L0'

图1

图1   列车行驶简图

Fig. 1   Moving train schematic


考虑如图1所示的列车模型,其车厢的数量为N,包括机车和普通两类车厢。在列车行驶过程中,每个轮对的荷载通过轨道及轨枕传递到路基顶面。对于某一关注轨枕,列车对路基的竖向作用力F(t)的表达式为

F(t)=i=1Nfi(t)   (i=1,2,,N)

式中:fi(t)表示第i节车厢的荷载传递到关注轨枕下方的部分,可进一步表示为该车厢所有轮对产生的竖向轨枕力之和。若该车厢为机车车厢,则其表达式为

fi(t)=j=1wPij(t)K(xij)

式中:w为车厢轮对数量,其中机车和普通两类车厢的轮对数分别为6和4;Pij(t)为第i节车厢第j组轮对在轨道上产生的荷载;K(xij)为关注轨枕对轮对作用力的分担系数,按日本规范的经验假设轮对荷载由附近5根轨枕按0.1∶0.2∶0.4∶0.2∶0.1分担25,可进一步表示为

K(xij)=0.4,xij-5Ls2Ls20.2,Ls2<xij-5Ls23Ls20.1,3Ls2<xij-5Ls25Ls20,xij-5Ls2>5Ls2

式中:Ls为轨枕间距,如图2所示;xij为轮对在Ox坐标系中的坐标。当第i节车厢为机车车厢时,xij具体表示为

图2

图2   道床顶面荷载沿线路方向分布图

Fig. 2   Distribution of loading on top surface of road bed along the road


xi1=-L0'+vt-k=1iLd,k+L0+4ai+bixi2=-L0'+vt-k=1iLd,k+L0+3ai+bixi3=-L0'+vt-k=1iLd,k+L0+2ai+bixi4=-L0'+vt-k=1iLd,k+L0+2aixi5=-L0'+vt-k=1iLd,k+L0+aixi6=-L0'+vt-k=1iLd,k+L0

式中:aii节车厢同一个转向架内相邻两个轮对的轴距;bi为第i节车厢不同转向架相邻两个轮对的轴距;Ld,k为第k节车厢长度;L0为车厢后端与该车厢最后一组轮对水平方向上的距离;-L0'为列车前端在该坐标系内的初始坐标。当第i节车厢为普通车厢时,也可以参考图1采用类似的方式计算轮对坐标。

列车第i节车厢第j组轮对产生的作用力Pij(t)可以采用反映行车平稳性、附加动载和轨面波形磨耗效应三个控制条件的激励力来模拟24,其表达式为

Pij(t)=Pij(0)+Pij(1)sinω1t+Pij(2)sinω2t+Pij(3)sinω3t

式中:Pij(0)为轮对静载;Pij(1)Pij(2)Pij(3)分别为反映前述三个控制条件的某一典型值的振动荷载幅值;ω1ω2ω3分别为反映控制条件的轨道几何不平顺曲线波长下的振动圆频率;t为荷载作用时间。其中Pij(k)ωk的表达式分别为

Pij(k)=Mij(0)hkωk2ωk=2πv/lk    (k=1, 2, 3)

式中:Mij(0)为列车簧下质量;lkhk分别为反映控制条件的轨道几何不平顺曲线波长和正矢;v为列车的运行速度。

轨道几何不平顺曲线波长及相应的正矢和列车行驶速度都具有一定的不确定性14。可以引入相应的随机变量,以量化这些参数的不确定性。对描述相关参数的随机变量进行抽样,将样本代入式(5)计算轮对在轨道上产生的作用力,进而根据式(1)可以计算所关注轨枕的竖向作用力,即可获得后续分析所需的随机列车荷载样本。

2 冻土路基随机振动分析的时域显式蒙特卡罗模拟法

作为一种多孔结构的物质,冻土中包含了土、水、冰、气四种成分,其动态力学行为较为复杂。为了实现对冻土路基这一复杂工程结构的随机振动分析,需要在反映冻土动力响应本质特征的前提下,对问题进行简化以降低建模难度和提升分析计算效率。冻土力学性质对季节性温度变化十分敏感,然而本文关注的列车荷载作用下冻土路基动力响应分析问题最多仅考虑数分钟的时程,因此可将列车荷载作用期间冻土温度及力学性质视为不变26-27。此外,列车荷载引起的冻土应变相对较小,特别是在冻结状态下,一般不会造成路基下方冻土的塑性变形28。因此,为降低后续非平稳随机振动分析的难度,可将冻土路基中各土层假定为各向同性线弹性材料29-30

2.1 动力响应显式表达式

根据上述假定,可将列车荷载作用下冻土路基的动力响应分析问题简化为线性非平稳随机振动问题。时域显式蒙特卡罗模拟法是解决这类问题的高效手段。在该方法中,首先需要建立结构响应的时域显式表达式。将冻土路基离散为具有nd个自由度的有限元模型,在非平稳随机激励作用下该结构的运动微分方程可表示为

MU¨+CU˙+KU=LF(t)

式中:M,C,K分别表示冻土路基结构的质量矩阵、阻尼矩阵和刚度矩阵;U,U˙,U¨分别为位移、速度和加速度向量;L为定位随机激励的nd×1阶常向量;F(t)为具有非平稳性的随机列车荷载。

将时间步长记为Δt,时程积分步数记为n,并把列车荷载F(t)离散为各时刻t1,t2,,tn处的荷载F1,F2,,Fn,定义状态向量为V=[UTU˙T]T。采用任意一种数值积分方法求解式(7),均可得到路基结构动力响应关于各时刻列车荷载的显式表达式。不失一般性,可假设结构初始状态V0=V(0)=0以及荷载F0=F(0)=0,则ti=iΔt时刻结构的状态向量Vi=V(ti)的显式表达式为

Vi=Ai,1F1+Ai,2F2++Ai,iFi    (i=1,2,,n)

式中:Ai,1,Ai,2,,Ai,i为只与路基结构参数有关的系数向量,反映结构参数对结构动力响应的影响。

参考文献[18]和[31]推导了系数向量的闭合公式,如下式所示:

A1,1=Q2   A2,1=TQ2+Q1Ai,1=TAi-1 ,1   (3in)Ai,j=Ai-1 ,  j-1    (2jin)

式中:T,Q1Q2的表达式由求解式(7)时所采用的积分格式决定,当采用Newmark-β积分格式31,可进一步表示为

T=H11  H12H21  H22,Q1=R1R3L,Q2=R2R4LH11=K̂-1(S1-S3M-1K),H12=K̂-1(S2-S3M-1C)H21=a3(H11-I)+a5M-1K,H22=a3H12-a4I+a5M-1CR1=K̂-1S3M-1,R2=K̂-1,R3=a3R1-a5M-1,R4=a3R2K̂=K+a0M+a3CS1=a0M+a3C,S2=a1M+a4C,S3=a2M+a5Ca0=1/(βΔt2),a1=1/(βΔt),a2=1/(2β)-1a3=γ/(βΔt),a4=γ/β-1,a5=(γ/β-2)Δt/2

式中:γ=0.5,β=0.25为Newmark-β算法中的计算参数。

式(9)展示了各系数向量间的内在联系。根据该关系,可将各时刻结构响应显式表达式所需的系数向量排列成如表1所示。由表1可以看出,仅需计算第一列系数向量Ai,1(i=1,2,,n)即可获得所有时刻响应的显式表达式,其计算量与1次确定性时程分析的计算量相当。该列系数向量的物理含义如图3所示,可方便地通过在结构上施加单位三角脉冲激励进行时程分析求得。在存储量方面,由于结构随机响应分析中并不需要关注结构的全部响应,因此Ai,1(i=1,2,,n)中的元素不需要全部储存。假设所关心的结构响应量个数为m,则需要存储的元素个数仅为2mn,其中n为时程分析步数。可见,需储存的系数向量元素个数与结构自由度数无关,即使对于大型复杂路基结构随机振动问题,在矩阵元素的存储方面也不会出现问题。从本质上看,这是由于结构的物理演变机制已经在式(8)中得到全面反映,结构响应已经全部解耦,可以仅对所关心的响应量进行后续随机分析而不涉及其他响应量,这是采用时域显式表达式的主要优势之一。

表1   各时刻结构响应对应的系数向量

Table 1  Coefficient vectors for structural responses at each time instant

时刻系数向量
F1F2F3Fn-2Fn-1Fn
t1A1,1-----
t2A2,1A1,1----
t3A3,1A2,1A1,1---
tn-2An-2,1An-3,1An-4,1A1,1--
tn-1An-1,1An-2 ,1An-3,1A2,1A1,1-
tnAn,1An-1,1An-2 ,1A3,1A2,1A1,1

新窗口打开| 下载CSV


图3

图3   在时刻t1施加全三角单位脉冲激励

Fig. 3   Unit impulse excitation at time t1


2.2 关键响应均值、标准差和平均峰值的计算

在工程结构分析中,通常只关注少量的结构关键响应。对于结构某一关键响应r,由式(8)可以得到关于r的一维显式表达式为

ri=qVi=ai,1F1+ai,2F2++ai,iFi   (i=1,2,,n)

其中

ai,j=qAi,j   (1ji)

式中:q为由状态向量V转换得到关键响应r的转换行向量,当r为某一位移或速度响应时,q由元素0和1组成;当r为某一应力响应时,q依赖于相应的本构关系。

通过本文第1节所述方法生成大量随机列车荷载样本,假设一共得到M个列车荷载F(t)的样本。在第k个向量样本Fk(t) (k=1,2,,M)作用下,由式(11)可以得到关键响应r在各时刻的值为

rik=ai,1F1k+ai,2F2k++ai,iFik
(i=1,2,,n;k=1,2,,M)

则关键响应r的均值、标准差和平均峰值分别为

μri=1Mk=1Mrik(i=1,2,,n)σri=1M-1k=1M(rik-μri)2(i=1,2,,n)rpeak=1Mk=1Mmaxi=1nrik(i=1,2,,n)

可见,利用所建立的时域显式表达式,结合蒙特卡罗模拟,可以高效地计算关键响应的统计值,这种方法称为时域显式蒙特卡罗模拟法。而在传统蒙特卡罗模拟法中,对于每一个列车荷载样本,均需要根据时程分析法求解式(7)所示的运动微分方程。因此,时域显式蒙特卡罗模拟法的计算量远小于传统蒙特卡罗模拟法的计算量,从而可以使蒙特卡罗模拟法直接应用于冻土路基的随机振动分析。

3 数值算例

3.1 随机列车荷载

本文主要考虑轨道几何不平顺曲线波长、相应的正矢以及列车行驶速度的不确定性。引入服从高斯分布的随机变量2432-33,以量化轨道不平顺波长、正矢和列车行驶速度的不确定性,如表2所示。列车荷载模拟过程中采用的时间步长为Δt=0.005 s,考虑列车编组的方式为2节NJ2机车+4节YZ25T客车,即式(1)中N=6,列车模型参数如表3所示。图1中列车前端的初始坐标可以假定为L0'=0图2中轨枕间距为Ls=0.556 m。对表2所列的随机变量进行抽样,将一组样本代入式(5)和式(6)中,就可以得到轮对在轨道上产生的作用力,然后根据式(1)和式(2)即可得到所关注轨枕的竖向作用力,获得一个荷载样本。

表2   波长、正矢和车速随机变量

Table 2  Random variables of wavelength, versine and vehicle speed

参数物理含义均值标准差
l1反映行车平稳性的轨道不平顺波长10 m1 m
h1反映行车平稳性的轨道不平顺正矢3.5 mm0.35 mm
l2反映附加动载的轨道不平顺波长2 m0.2 m
h2反映附加动载的轨道不平顺正矢0.4 mm0.04 mm
l3反映轨面波形磨耗效应的轨道不平顺波长0.5 m0.05 m
h3反映轨面波形磨耗效应的轨道不平顺正矢0.08 mm0.008 mm
v列车行驶速度94 km·h-110 km·h-1

新窗口打开| 下载CSV


表3   青藏铁路列车模型参数34

Table 3  Vehicle model parameters of Qinghai-Tibet Railway34

参数符号含义单位客车机车
YZ25TNJ2
Pij(0)轮对静载kN147215.11
Mij(0)簧下质量kg1 9002 750
Ld,k车厢长度m26.622
ai

同一转向架相邻

两个轮对的轴距

m0.9151.25
bi

不同转向架相邻

两个轮对的轴距

m15.68.4
L0车厢后端与该车厢最后一组 轮对水平方向上的距离m1.20.8

新窗口打开| 下载CSV


为验证本文提出的随机列车荷载模拟方法的准确性,将生成的一个荷载样本与文献[15]中采用车辆-轨道耦合动力学计算所得的轨枕竖向作用力时程进行对比,如图4所示。由图4可知,本文方法得到的荷载样本在幅值和波形上与文献[15]的方法得到的荷载时程结果吻合,说明本文所提的随机列车荷载模拟方法是可行的。

图4

图4   列车荷载时程样本

Fig. 4   Time history sample of train load


3.2 路基结构参数

以青藏铁路某传统道砟路基为研究对象,计算模型如图5所示。计算区域中的土层依次为道碴层、路基填土层、砂土层、粉质黏土层及弱风化岩层,如图5(a)所示。冻土的力学参数不仅依赖于各土层材料,还取决于土层的温度。本文选取青藏铁路路基修建完成后第10年夏冬两个典型季节气温最高时路基内部的温度分布,根据参考文献[7]中的拟合公式计算路基各土层的力学参数,结果见表4~5

图5

图5   二维冻土路基计算模型

Fig. 5   Two-dimensional permafrost embankment model: geometry of embankment7 (a), finite-element model of permafrost embankment (b)


表4   路基各土层的参考温度和力学参数(夏季)

Table 4  Mechanical parameters of the soils in embankment model (Summer)

计算参数密度ρ/(kg⋅m-3温度T/℃弹性模量E/MPa泊松比μ
道碴层2 1006.0200.000.30
路基填土层2 00010.061.000.35
砂土层2 0502.041.000.35
亚黏土层1 920028.000.40
弱风化岩层2 000-0.5211.250.25

新窗口打开| 下载CSV


表5   路基各土层的参考温度和力学参数(冬季)

Table 5  Mechanical parameters of the soils in embankment model (Winter)

计算参数密度ρ/(kg⋅m-3温度T/℃弹性模量E/MPa泊松比μ
道碴层2 100-8.0200.000.30
路基填土层2 000-2.0141.330.34
砂土层2 050-0.573.400.34
亚黏土层1 920028.000.40
弱风化岩层2 000-0.5211.250.25

新窗口打开| 下载CSV


根据表4~5所得的各土层的力学参数,在通用有限元软件ANSYS中建立二维冻土路基有限元计算模型。建模时按平面应变问题考虑34-35,采用PLANE42单元,模型单元数为3 500,节点数为3 631,总自由度数为7 262,如图5(b)所示。考虑在路基两侧及底部布置黏弹性人工边界,在ANSYS中采用COMBIN14单元进行模拟,单元的刚度和阻尼参数根据文献[36]所述的方法确定。分析过程中选用瑞利阻尼7,阻尼参数为α=β=0.03,时程积分步长取Δt=0.005 s。在路基设计和监测过程中,路基填土层顶面和底面的竖向位移、竖向速度和竖向正应力是设计人员较为关注的关键响应,因此选择路基填土层顶面和底面的中点作为关注位置,分别记为1#点和2#点,如图5所示。由于时域显式表达式独特的降维计算优势,在后续计算过程中可以只针对这两个关注位置进行降维计算,从而大幅度提升蒙特卡罗模拟的效率。

3.3 计算结果讨论

为验证所建立的时域显式表达式的计算精度,本文分别采用时域显式表达式和ANSYS直接时程分析法计算冻土路基在一个列车荷载样本作用下的动力响应,结果如图6所示。由图6可知,采用两种计算方法得到的结果一致,从而验证了时域显式表达式在冻土路基动力响应分析中的计算精度。

图6

图6   路基1#点竖向动力响应时程(夏季)

Fig. 6   Time histories of dynamic responses at point 1# (Summer): time history of vertical displacement (a), time history of vertical velocity (b), time history of vertical normal stress (c)


为展示时域显式蒙特卡罗模拟法相比于传统蒙特卡罗模拟法的效率优势,分别采用了不同数量的样本作为输入进行随机分析,两种方法所需计算时间如表6所示。由第2.1节中的分析可知,时域显式蒙特卡罗模拟法的计算时间包括求取系数向量的计算时间和实施蒙特卡罗模拟的计算时间两部分。由于该方法在实施蒙特卡罗模拟时无需反复求解冻土路基结构运动微分方程,第二部分耗时很短。而传统蒙特卡罗模拟法在每次样本分析中均需求解路基结构运动微分方程,计算非常耗时,其效率远低于时域显式蒙特卡罗模拟法的计算效率。由表6还可以看出,随着样本数的增大,传统蒙特卡罗模拟法计算时间呈线性增长,而时域显式蒙特卡罗模拟法的计算时间增长不多,说明当样本数量规模较大时,时域显式蒙特卡罗模拟法的计算优势更加明显。

表6   两种不同方法的计算时间比较

Table 6  Comparison of time costs between two methods

方法样本数计算时间

传统

蒙特卡罗模拟法

1350 s
500350 s × 500 ≈ 49 h
1 000350 s × 1 000 ≈ 97 h
2 000350 s × 2 000 ≈ 194 h

时域显式

蒙特卡罗模拟法

1350 s + 0.01 s = 350.01 s
500350 s + 0.23 s = 350.23 s
1 000350 s + 0.52 s = 350.52 s
2 000350 s + 0.96 s = 350.96 s

新窗口打开| 下载CSV


当随机激励样本数分别取500、1 000和2 000时采用时域显式蒙特卡罗模拟法获得的夏季冻土路基竖向位移、竖向速度和竖向正应力计算结果同时显示于图7中。由图7可以发现,当样本量达到 1 000时,冻土路基各随机响应的均值和标准差结果已经收敛。由表6可知,采用传统蒙特卡罗模拟法处理1 000个样本耗时约4 d,而时域显式蒙特卡罗模拟法耗时约6 min。

图7

图7   不同样本数下路基1#点2#点竖向动力响应统计矩时程(夏季)

Fig. 7   Statistical moment time histories of dynamic responses at point 1# and point 2# with different numbers of samples (Summer): mean value time histories of vertical displacements at point 1# and point 2# (a), mean value time histories of vertical velocities at point 1# and point 2# (b), mean value time histories of vertical normal stresses at point 1# and point 2# (c), stadnard deviation time histories of vertical displacements at point 1# and point 2# (d), standard deviation time histories of vertical velocities at point 1# and point 2# (e), standard deviation time histories of vertical normal stresses at point 1# and point 2# (f)


由时域显式蒙特卡罗模拟法还可以得到夏季和冬季冻土路基1#点和2#点动力响应平均峰值,结果如表7所示。由表7可知,列车荷载作用下冻土路基的动力响应不仅与路基深度有关,还受季节变化影响。总的来说,夏季时路基浅层的位移和速度响应最为剧烈,这是因为冻土是一种对温度敏感的土体。夏季温度相对冬季较高,导致冻土的弹性模量等力学参数降低,从而导致了路基中更为剧烈的位移和速度响应。

表7   夏季和冬季冻土路基响应平均峰值

Table 7  Mean peak value of permafrost embankment in Summer and Winter

响应季节关注位置平均峰值
竖向位移夏季1#1.5×10-3 m
2#0.86×10-3 m
冬季1#0.96×10-3 m
2#0.70×10-3 m
竖向速度夏季1#2.71×10-2 m·s-1
2#1.22×10-2 m·s-1
冬季1#1.54×10-2 m·s-1
2#0.93×10-2 m·s-1
竖向正应力夏季1#2.05×104 Pa
2#7.19×103 Pa
冬季1#2.08×104 Pa
2#6.65×103 Pa
横向正应力夏季1#5.05×103 Pa
2#1.59×103 Pa
冬季1#6.95×103 Pa
2#3.46×103 Pa

新窗口打开| 下载CSV


4 结论

冻土路基随机振动问题十分复杂,目前对冻土路基随机振动的理论研究较少。本文对上述问题展开了探索性的分析与研究,取得了一定的研究成果。基于以上分析研究,得到的主要结论如下:

(1)采用本文提出的随机列车荷载模拟方法,可以快速生成冻土路基随机振动分析所需的大量列车荷载样本用于蒙特卡罗模拟。

(2)本文采用时域显式蒙特卡罗模拟法计算了冻土路基在两个典型季节下的动力响应平均峰值,发现路基的位移和速度响应在夏季最为剧烈。

(3)对比了时域显式蒙特卡罗模拟法和传统蒙特卡罗模拟法在处理不同数量的样本时所需的计算时间,展示了时域显式蒙特卡罗模拟法在处理冻土路基随机振动问题中的计算精度和效率。

(4)本文中的时域显式表达式仅适用于线性系统,因此将冻土理想化为线弹性材料,这一理想化带来的误差尚需进一步分析。

由于目前尚无随机列车荷载的现场实测数据,因此本文仅从计算理论上进行初步探索,以期能为冻土路基的设计、维护及研究提供理论依据和参考。冻土是一种多孔多相材料,特别是在夏季高温时其力学行为十分复杂。为反映这一特性,后续研究将结合等效线性化方法,考虑非线性冻土本构关系,计算冻土路基永久变形和可靠度。

参考文献

Zhou YouwuGuo DongxinQiu Guoqinget al. Geocryology in China[M]. BeijingScience Press20001-9.

[本文引用: 1]

周幼吾郭东信邱国庆. 中国冻土[M]. 北京科学出版社20001-9.

[本文引用: 1]

Zhang Shimin.

Numerical study on the thermal-mechanical stability of Qinghai-Tibet Railway embankment in permafrost regions

[J]. Journal of Civil, Architectural & Environmental Engineering, 201234(): 82-87.

[本文引用: 1]

张世民.

青藏铁路多年冻土路基热-力稳定性数值仿真分析

[J]. 土木建筑与环境工程, 201234(): 82-87.

[本文引用: 1]

Ma WeiLiu DuanWu Qingbai.

Monitoring and analysis of embankment deformation in permafrost regions of Qinghai-Tibet Railway

[J]. Rock and Soil Mechanics, 2008293): 571-579.

[本文引用: 1]

马巍刘端吴青柏.

青藏铁路冻土路基变形监测与分析

[J]. 岩土力学, 2008293): 571-579.

[本文引用: 1]

Zhu ZhanyuanLing XianzhangZhang Fenget al.

Field monitoring on vibration response of railway structure in the seasonally frozen region in summer

[J]. Journal of Harbin Institute of Technology, 20094112): 41-45.

[本文引用: 1]

朱占元凌贤长张锋.

季节冻土区夏季轨道结构振动反应现场监测

[J]. 哈尔滨工业大学学报, 20094112): 41-45.

[本文引用: 1]

Stevens H W.

The response of frozen soils to vibratory loads: CRREL Technical Report 265

[R]. Hanover, NHUS Army CRREL1973.

[本文引用: 1]

Kong XiangbingZhao ShupingMu Yanhuet al.

Research on the calculation of dynamic stress of embankment in permafrost regions under train load

[J]. Journal of Glaciology and Geocryology, 2013356): 1490-1498.

[本文引用: 1]

孔祥兵赵淑萍穆彦虎.

列车荷载作用下冻土路基中的动应力计算研究

[J]. 冰川冻土, 2013356): 1490-1498.

[本文引用: 1]

Li ShuangyangLai YuanmingZhang Shujuanet al.

Dynamic responses of Qinghai-Tibet Railway embankment subjected to train loading in different seasons

[J]. Soil Dynamics and Earthquake Engineering, 2012. 321): 1-14.

[本文引用: 6]

Li ShuangyangZhang MingyiZhang Shujuanet al.

Analysis of the dynamic response of Qinghai-Tibet Railway embankment in permafrost regions under train load

[J]. Journal of Glaciology and Geocryology, 2008305): 860-867.

李双洋张明义张淑娟.

列车荷载下青藏铁路冻土路基动力响应分析

[J]. 冰川冻土, 2008305): 860-867.

Huang ZhijunLai YuanmingLi Shuangyanget al.

Dynamic response of embankment in permafrost regions under traffic load

[J]. Journal of Glaciology and Geocryology, 2012342): 418-426.

[本文引用: 1]

黄志军赖远明李双洋.

交通荷载作用下冻土路基动力响应分析

[J]. 冰川冻土, 2012342): 418-426.

[本文引用: 1]

Liu ChaoChe AilanWu Zhijianet al.

Stability of embankments in permafrost regions considering temperature field coupled with dynamic field

[J]. Chinese Journal of Geotechnical Engineering, 201133(): 466-471.

[本文引用: 1]

刘超车爱兰吴志坚.

温度及机车载荷耦合作用下多年冻土区路基稳定性研究

[J]. 岩土工程学报, 201133(): 466-471.

[本文引用: 1]

Wu ZhijianChen TuoMa Wei.

Creep analysis of plain fill embankment at permafrost regions along Qinghai-Tibet Railway under train dynamic load

[J]. Rock and Soil Mechanics, 201132(): 83-87.

吴志坚陈拓马巍.

机车荷载作用下青藏铁路多年冻土区普通路基的蠕变分析

[J]. 岩土力学, 201132(): 83-87.

Chen TuoWu ZhijianMu Yanhuet al.

Dynamic behavior of the Qinghai-Tibetan Railway embankment in permafrost regions under trained-induced vertical loads

[J]. Sciences in Cold and Arid Regions, 2020124): 189-199.

Wu ZhijianChen TuoZhao Taoet al.

Dynamic response analysis of railway embankments under train loads in permafrost regions of the Qinghai-Tibet Plateau

[J]. Soil Dynamics and Earthquake Engineering, 20181121): 1-7.

[本文引用: 1]

Ma MengLiu WeiningLiu Weifeng.

Research progress of prediction method and uncertainty of train-induced environmental vibration

[J]. Journal of Traffic and Transportation Engineering. 2020203): 1-16.

[本文引用: 2]

马蒙刘维宁刘卫丰.

列车引起环境振动预测方法与不确定性研究进展

[J]. 交通运输工程学报, 2020203): 1-16.

[本文引用: 2]

Zhu ZhanyuanLing XianzhangChen Shijunet al.

Investigation of train -induced vertical applied force on sleeper along Qinghai-Tibet Railway

[J]. Journal of Harbin Institute of Technology, 2011436): 6-10.

[本文引用: 4]

朱占元凌贤长陈士军.

青藏铁路列车行驶引起的轨枕竖向作用力研究

[J]. 哈尔滨工业大学学报, 2011436): 6-10.

[本文引用: 4]

Lin JiahaoZhang Yahui. Pseudo-excitation method of random vibration [M]. BeijingScience Press2004.

[本文引用: 1]

林家浩张亚辉. 随机振动的虚拟激励法[M]. 北京科学出版社2004.

[本文引用: 1]

Li JieChen Jianbing. Stochastic dynamics of structures[M]. SingaporeJohn Wiley and Sons2009.

[本文引用: 1]

Su ChengXu Rui.

Time-domain method for dynamic reliability of structural systems subjected to non-stationary random excitations

[J]. Chinese Journal of Theoretical and Applied Mechanics, 2010433): 512-520.

[本文引用: 2]

苏成徐瑞.

非平稳随机激励下结构体系动力可靠度时域显式解法

[J]. 力学学报, 2010433): 512-520.

[本文引用: 2]

Su ChengHuang ZhijianLiu Xiaolu.

Time-domain explicit random simulations method for seismic analysis of tall buildings

[J]. Journal of Building Structures, 2015361): 13-22.

[本文引用: 1]

苏成黄志坚刘小璐.

高层建筑地震作用计算的时域显式随机模拟法

[J]. 建筑结构学报, 2015361): 13-22.

[本文引用: 1]

Su ChengLi BaomuChen Taiconget al.

Stochastic optimal design of nonlinear viscous dampers for large-scale structures subjected to non-stationary seismic excitations based on dimension-reduced explicit method

[J]. Engineering Structures, 2018175217-230.

[本文引用: 1]

Su ChengWu ZhenzhenXian Jianhua.

Stochastic dynamic analysis of vehicle-bridge coupled systems with nonlinear Hertz contacts by explicit time-domain method

[J]. Vehicle System Dynamics, 2020201-21.

[本文引用: 1]

Liang BoCai Ying.

Dynamic analysis on subgrade of high speed railways in geometric irregular condition

[J]. Journal of the China Railway Society, 1999212): 93-97.

[本文引用: 1]

梁波蔡英.

不平顺条件下高速铁路路基的动力分析

[J]. 铁道学报, 1999212): 93-97.

[本文引用: 1]

Feng Yue.

Dynamic response of warm permafrost subgrade under train-induced dynamic loads

[D]. HarbinHarbin Institute of Technology2012.

冯悦.

列车荷载下高温冻土场地路基的动力反应分析

[D]. 哈尔滨哈尔滨工业大学2012.

Liang BoLuo HongSun Changxin.

Simulated study on vibration load of high speed railway

[J]. Journal of the China Railway Society, 2006284): 89-94.

[本文引用: 3]

梁波罗红孙常新.

高速铁路振动荷载的模拟研究

[J]. 铁道学报, 2006284): 89-94.

[本文引用: 3]

Liu JingleiYe QingzhiSong Xuguoet al.

Research on the load condition and dynamic properties for heavy haul railway subgrade

[J]. Journal of Railway Engineering Society. 2015322): 33-38 53.

[本文引用: 1]

刘晶磊叶庆志宋绪国.

重载铁路路基荷载条件及动力特性研究

[J]. 铁道工程学报, 2015322): 33-38 53.

[本文引用: 1]

Gu HaidongTang AipingYang Min.

Dynamic response analysis of embankment under traffic load in permafrost regions

[J]. Low Temperature Architecture Technology, 2010328): 85-87.

[本文引用: 1]

古海东汤爱平杨敏.

交通荷载作用下冻土路基动力反应分析

[J]. 低温建筑技术, 2010328): 85-87.

[本文引用: 1]

Chang LiwuXu Yanjie.

Dynamic response analysis of high temperature permafrost subgrade

[J]. Subgrade Engineering, 20074): 34-35.

[本文引用: 1]

常利武徐艳杰.

高温冻土路基动力响应的数值模拟研究

[J]. 路基工程, 20074): 34-35.

[本文引用: 1]

Shi Yehui.

Study on stability of the subgrade under train loading and freeze-thaw

[D]. BeijingBeijing Jiaotong University2011.

[本文引用: 1]

施烨辉.

列车荷载和冻融循环作用下冻土路基稳定性研究

[D]. 北京北京交通大学2011.

[本文引用: 1]

Liu JiankunLiu FengxiFang Jianhong.

Modeling and monitoring the train-induced vibration in permafrost in Reshui Coal Mining, Qinghai

[J]. Journal of Glaciology and Geocryology, 2004262): 177-181.

[本文引用: 1]

刘建坤刘奉喜房建宏.

青海热水煤矿多年冻土区列车引起的地面振动检测与模拟

[J]. 冰川冻土, 2004262): 177-181.

[本文引用: 1]

Li Tao.

Dynamic analysis of subgrade structure of Qinghai-Tibet Railway in permafrost region

[J]. Journal of Railway Engineering Society, 20073): 29-32.

[本文引用: 1]

李涛.

青藏铁路多年冻土区路基结构的动力分析

[J]. 铁道工程学报, 20073): 29-32.

[本文引用: 1]

Su ChengHuang HuanMa Haitao.

Fast equivalent linearization method for nonlinear structures under nonstationary random excitations

[J]. Journal of Engineering Mechanics, 20161428): 1-15.

[本文引用: 1]

Li Chenghui. Track engineering[M]. ChongqingChongqing University Press2014.

[本文引用: 1]

李成辉. 轨道工程[M]. 重庆重庆大学出版社2014.

[本文引用: 1]

Chen WenZhou LiangXu Yudeet al.

Analysis of testing data on geometric track irregularity in urban rail transit

[J]. Urban Mass Transit, 2013161): 28-32.

[本文引用: 1]

陈雯周亮许玉德.

城市轨道交通轨道几何不平顺检测数据的应用分析

[J]. 城市轨道交通研究, 2013161): 28-32.

[本文引用: 1]

Zhu Zhanyuan.

Train-induced vibration response and subsidence prediction of permafrost subgrade along Qinghai-Tibet Railway

[D]. HarbinHarbin Institute of Technology2009.

[本文引用: 3]

朱占元.

青藏铁路列车行驶多年冻土场地路基振动反应与振陷预测

[D]. 哈尔滨哈尔滨工业大学2009.

[本文引用: 3]

Wang Ziyu.

Dynamic response characteristics and permanent deformation of subgrade induced by train load in seasonally deep frozen region

[D]. HarbinHarbin Institute of Technology2014.

[本文引用: 1]

王子玉.

深季节冻土区列车荷载下路基振动响应特性与永久变形研究

[D]. 哈尔滨哈尔滨工业大学2014.

[本文引用: 1]

Chen BaokuiWang DongshengCheng Hu.

Research review on the application of viscous-spring artificial boundary in earthquake engineering

[J]. Journal of Seismological Research, 2016391): 137-142 182.

[本文引用: 1]

陈宝魁王东升成虎.

粘弹性人工边界在地震工程中应用研究综述

[J]. 地震研究, 2016391): 137-142 182.

[本文引用: 1]

/