Bearing Characteristics of X-Shaped Sleeper Based on Scaled Test and Discrete Element Simulation
-
摘要:
为研究新型X形轨枕道床承载力和横向阻力特性,开展X形轨枕与Ⅲ型轨枕道床刚度和横向阻力的缩尺模型试验,通过离散元法建立有砟轨道三维模型,并从细观上对比分析2种轨枕道床的竖向荷载传递机制及横向阻力. 研究结果表明:在最大竖向荷载时,X形轨枕道床的竖向位移(刚度)相较于Ⅲ型轨枕降低了约26.3%(提高了约46.6%);与Ⅲ型轨枕相比,X形轨枕的横向极限阻力提高了约22.4%,有效提高了轨道横向稳定性; X形轨枕与枕间道砟的接触面积和应力均明显增加,轨枕接触力沿X形4个夹角范围扩散,使得枕间道砟充分参与分担荷载;由于X形轨枕的叉形结构能够提升枕间道砟的参与,使得道床刚度和横向阻力分别提升约29.2%和31.6%,与试验结论较接近.
Abstract:To investigate the trackbed bearing capacity and lateral resistance characteristics of a new X-shaped sleeper, scaled tests comparing the stiffness and lateral resistance between X-shaped and Type Ⅲ sleepers were conducted. A 3D model of a ballast track was established by the discrete element method to analyze the vertical load transmission mechanism and lateral resistance of these two types of sleepers at a micro level. The results indicate that at the maximum vertical load, the X-shaped sleepers significantly reduce vertical displacement (stiffness) by approximately 26.3% compared to Type Ⅲ sleepers (an increase of about 46.6%). Furthermore, the ultimate lateral resistance of the X-shaped sleeper is increased by 22.4%, which effectively improves the lateral stability of the track. The X-shaped sleepers exhibit a substantial increase in the contact area and stress with ballast between the sleepers. The contact forces on the X-shaped sleeper are distributed over four angular segments, making the ballast between the sleepers fully participate in the load sharing. Because the structure of the X-shaped sleeper can increase the participation of the ballast between the sleepers, the stiffness and transverse resistance of the trackbed are increased by about 29.2% and 31.6%, respectively, which is close to the experimental conclusion.
-
Key words:
- railway sleeper /
- discrete element method /
- trackbed stiffness /
- lateral resistance /
- scaled test
-
高层建筑作为典型的风敏结构,建筑上的风荷载与其气动外形密切相关[1-2]. 随着建造技术的发展,越来越多的辅助构件,如阳台、外伸板、粗糙条[3-4]等被安装在高层建筑外立面. 这些附属物的存在不仅影响建筑的气动外形,也会改变建筑表面风压分布、风荷载及周围风环境. 因此,弄清附属物对建筑风效应的影响十分必要.
部分学者通过风洞试验探索了附属物对风压和风荷载的影响. 阳台可以显著改变迎风面上部的平均和峰值风压分布[5-6],建筑物表面风压降低约30%[7]. 不同布置形式的肋板对风荷载的影响不同[8],连续水平板会略微增加顺风向平均基底弯矩,而间隔水平板能降低横风向脉动基底弯矩[9]. 迎风面边缘附近的竖向肋板会增加迎风角区域的吸力[10],同时,还会显著降低表面脉动风压,最大降幅约为57%[11].
虽然试验证明了附属物能减弱结构上的风荷载,但其抗风工作机理尚不清楚. 因此,为阐明附属物的流场修正机理,一些学者尝试利用流动可视化 技术(PIV)和LES (large eddy simulation)研究附属物对建筑周围流场特征的影响. Ai等[12]发现带栏杆的水平阳台能显著抑制迎风面和背风面附近的竖向流,但同时增强了近壁湍流脉动. Zheng等[13]认为相邻阳台间的局部旋涡是影响风速和风压分布的重要因素. 建筑边缘附近的阳台能有效抑制流动分离,并破坏涡脱过程[14].
除水平附属物外,Liu等[15]发现竖向肋板能显著降低侧风面附近湍流脉动,并将高湍流区推离侧风面;Liu等[16-17]认为迎风角竖向肋板附近形成的局部涡旋可以明显改变模型周围的流场结构;Hui等[18]发现不同位置的肋板对流场结构的影响明显不同,迎风面肋板抗风作用最明显,而侧风面和背风面肋板对流场结构的影响较弱.
尽管部分学者探究了外伸肋板的流场修正机理,但由于肋板布置的复杂性,水平肋板的抗风工作机理仍需进行系统的研究.
本文采用LES方法研究大气边界层来流(ABL)下水平肋板密度和布置形式对高层建筑流场和气动特性的影响,并详细分析肋板对风场特征的影响,如近壁流、涡结构和竖向流等,借助流场信息阐明肋板的抗风工作机理,为实际高层建筑外辅助构件的抗风设计和布置提供参考.
1. 数值方法
1.1 模型信息
以刚性模型为目标高层建筑,模型长、宽、高(D、B、H)的比例为1∶1∶4. 模型缩尺比为1∶200,模型高度和对应实际建筑高度分别为0.4 m和80 m. 外伸肋板安装在模型的外立面,如图1所示. 图中:d为外伸长度,b、h分别为2个相邻板间的水平和垂直距离, α为风向角. 模型1~3研究肋板密度的影响,模型2~5研究肋板布局的影响. 外伸肋板主要有3个参数, d固定为0.06D (实际尺寸为1.2 m),该尺寸与阳台的尺寸相似[6]. 肋板的具体尺寸信息如表1所示.
表 1 肋板尺寸信息Table 1. Dimension information of ribs模型编号 d/D/% b/D//% h/D//% 1 6 50.0 2 6 25.0 3 6 12.5 4 6 20 25.0 5 6 60 25.0 1.2 控制方程及求解策略
本次模拟选择LES湍流模型,亚格子尺度模型为标准SGS (subgrid scale)模型,Smagorinsky常数Cs取值为0.13,Cs在0.10~0.14范围内模拟结果较好[19]. 不可压缩流滤波后的N-S方程为
∂¯ui∂t+∂(¯ui¯uj)∂xj=−1ρ∂¯p∂xi+ν∂2¯ui∂xj∂xj−∂τij∂xj, (1) τij = −2νt¯Sij+13τkkδij, (2) ¯Sij = 12(∂¯ui∂xj+∂¯uj∂xi), (3) ∂¯ui∂xi=0, (4) 式中:xi (xj)为笛卡尔坐标系中x、y、z 3个方向的坐标值,¯ui(uj)为x、y、z 3个方向的速度,¯p为滤波后的压力,ρ、ν分别为空气密度和运动黏度,τij为SGS亚格子应力张量,νt为SGS涡黏性系数,¯Sij为滤波后的应变率张量,t为时间.
利用开源软件OpenFOAM进行瞬态模拟,求解器选择PimpleFoam. 在控制方程的求解中,数值模型离散格式选择如下:二阶后向差分方案用于时间离散,控制方程中的扩散项和对流项分别采用二阶中心差分格式和线性迎风稳定输运格式(LUST). PIMPLE算法用于解决压力-速度耦合问题. 几何代数多重网格(GAMG)和预处理双共轭梯度(PBiCG)求解器用于计算迭代过程中的压力和速度[20]. 本次模拟的时间步长设置为5 × 10−5,模拟中的平均库郎数(CFL)为0.008,最大库郎数小于1.200. 每个工况的总模拟时长为9 s.
1.3 网格信息和边界条件
计算域和网格信息如图2所示. 图中:①~④ 为加密区. 计算域入口处施加湍流来流(ABL),出口为自由流条件,侧壁和计算域顶部为滑移边界条件,建筑物表面和地面为无滑移边界条件,阻塞率为0.6%.
利用网格生成工具snappyHexMesh生成网格系统. 本次模拟采用4个加密区实现网格的逐级加密(图2(a)、(b)),加密区 ①~④ 尺寸分别为9H × 3.5H × 2H、1.75H × 1.75H × 1.5H、0.75H × 0.75H × 1.1H、0.04H × 0.04H × 1.03H. 计算域内主要采用结构化网格,非结构化网格只出现在加密区的交界处. 3种不同的网格布局用于检查网格无关性,如表2所示. 网格2的近壁首层网格在水平方向上高度为1/125D. 建筑物附近竖直方向的网格尺寸远小于水平方向的网格尺寸,近地面的网格尺寸为1/600D,网格2的数量约为7百万.
表 2 网格信息Table 2. Mesh information网格 最小网格 最大网格 数量/百万 x/y z x/y z 网格 1 1/62.5D 1/300D 0.6D 1D 3.7 网格 2 1/125.0D 1/600D 0.6D 1D 7.2 网格 3 1/250.0D 1/ 1200 D0.6D 1D 13.6 使用一致离散随机流生成(CDRFG)方法[21]生成入口的脉动风速时程. 该方法基于速度功率谱的离散化和合成技术,可以很好地保持湍流能量和速度相干性,也可以准确评估作用在建筑物上的风荷载[21]. 顺流向平均风速U、湍流强度Iu剖面[22-23]及顺风向脉动风速功率谱如图3所示:速度剖面指数为1/6,建筑物屋顶处的顺风向平均风速UH为8 m/s,对应的雷诺数Re=5.5 × 104,建筑物高度处的湍流积分尺度LuH为0.24 m. 目标建筑物处的风场特性决定结构上的风荷载. 由图3可知,平均风速和湍流强度剖面与已有文献[22-23]结果吻合较好,但由于LES自身的数值耗散,模型位置处的湍流强度相比入口处有一定程度的减小;顺风向脉动风速谱也与von Karma谱[20]一致,表明CDRFG方法生成的来流湍流是可靠的.
1.4 气动力定义和结果验证
无量纲风压系数Cp、阻力系数Cd和升力系数Cl定义如下:
Cp=Δp1/2ρU2H, (5) Cd=Fd1/2ρU2HA, (6) Cl=Fl1/122ρU2HA, (7) 式中:Δp为压差,Fd、Fl分别为总阻力和总升力,A为投影面积.
大涡模拟的结果验证如图4[22-23]和表3所示. 图中:a、b、c、d表示模型边缘. 表中:AR为高宽比,¯Cd、Cd1为阻力的时均值和均方根值,Cl1为升力的均方根值. 可以看出,本次LES得到的平均风压系数¯CP的分布与已有研究结果接近[22-23],但LES得到的脉动风压均方根Cp1与已有结果略微差异,可能因为湍流生成方法、数值方案和网格分辨率等方面的差异造成的. 风压的差异会导致模型上风荷载不同(如表3所示). 与已有研究结果相比,本次模拟结果在合理范围内,表明本次模拟结果是可靠的. 网格1会高估脉动风压均方根,网格2、3的结果更接近已有研究结果. 因此,本研究采用网格2方案来平衡计算成本和结果精度.
表 3 参考方柱LES结果及验证Table 3. LES results of square column and verification方法 文献 Re/104 AR ¯Cd Cd1 Cl1 St LES 网格 1 5.5 4 1.220 0.09 0.330 0.100 网格 2 5.5 4 1.147 0.10 0.300 0.096 网格 3 5.5 4 1.156 0.11 0.310 0.094 试验 文献 [22] 5.5 4 1.200 0.410 0.100 文献 [23] 5.0 4 1.280 0.13 0.320 0.085 文献 [24] 6.0 4 1.270 0.08 0.296 0.097 文献 [25] 7.3 3 1.290 0.102 文献 [11] 5.5 5 0.930 0.350 0.100 文献 [10] 5.5 5 1.080 2. 外伸板对流场结构的影响
2.1 平均流场结构
大气边界层来流下的建筑物周围会产生明显的三维流动结构,不同高度处的平均流场如图5所示. 可以看出:不同高度处的旋涡大小差异明显;相比模型1~3,在z=0.52H高度处,参考模型的尾涡相对较短,尺寸较小,表明肋板能拉长尾涡;模型顶部自由端效应导致建筑物屋顶附近尾涡缺失(z=0.95 H).
图6所示为竖直中心面(y=0D)上的顺风向平均风速和平均流线,以阐明水平肋板对竖向流动的影响. 对于参考模型,来流在迎风面0.73H高度附近停滞. 当模型上安装水平肋板时,迎风面上的驻点相比参考模型无明显变化,肋板附近形成小尺寸局部旋涡.
肋板密度会明显改变尾流区鞍点高度,鞍点高度随密度增加而减小,模型2、3的尾部鞍点高度分别为1.17D和1.23D(图6(c)、(d)),远低于参考模型的高度(1.50D). 对比图6(c)、(d)的结果可以发现,沿建筑物高度的肋板密度较大时,肋板对尾部鞍点影响较小. 图6(c)中的点F表示流动再附点,局部旋涡再附在模型迎风面上.
另外,本节还研究了肋板水平间距(模型4和模型5)对平均流场的影响. 对比模型2、4、5的结果可以看出,随着水平间隙的增加,尾涡逐渐变小,尾部鞍点的高度明显变高(图6(c)、(e)~(f)),说明外伸板随着水平距离的增大,其对尾涡的影响逐渐变弱.
3. 外伸板对风压和风荷载的影响
3.1 外伸板密度影响
流场结构的变化会影响模型壁面风压分布(图7). 由图7(a)可知,参考模型中心线上的¯CP沿竖向平滑变化,¯CP最大值出现在0.8H左右. 当模型上安装外伸板时,¯CP分布变化明显,模型1~3的¯CP呈“之”字形分布,肋板附近出现¯CP的局部最大值和最小值. Hui等[9]也发现了迎风面和背风面¯CP的“之”字形模式. 这种结果可用图6所示的近壁流场来解释,因为肋板上、下形成的局部旋涡会导致局部压力的出现. 靠近模型自由端(0.8H~1.0H)或模型底部(0~0.3H)时,压力变化更显著(图7(a)),因为上、下部肋板附近形成的旋涡尺寸相对较大,大尺寸旋涡对¯CP的影响更大. ¯Cp的变化幅度随着肋板密度的增加而减小,因为密集的肋板明显影响近壁流并减小肋板附近局部旋涡的尺寸(图6(b)~(d)中的区域L). 较大的局部旋涡将导致模型1上出现较大Cp1. 对于迎风面,Cp1的变化也随着肋板密度的增加而减小,如图7(b)所示.
在z=0.52H处,侧壁面上的风压分布如图8所示. 图中:A、B为侧风面边缘. 4个模型在0.52H处的风压系数在前缘和后缘处基本保持不变,(图8(a)),但模型边缘附近风压系数因流动分离和次级涡的影响而发生突变. 外伸板能明显降低侧壁面风压,相比参考模型,模型1~3在0.52H处的¯CP和Cp1小于参考模型,¯CP和Cp1的最大降幅分别为20%和17%.
带板和无板光面模型的平均阻力系数¯Cd、Cd1和Cl1如表4所示,脉动力是荷载系数的标准差. 由表可知,外伸板虽然影响迎风面的压力分布,但其对阻力无明显影响,然而,带板模型1~3的Cl1明显较小,减幅约为27%,是因为模型1~3侧壁面风压较小(图8).
表 4 参考方柱和模型1~3的风荷载系数Table 4. Wind load coefficients for square column and models 1–3名称 ¯Cd Cd1 Cl1 方柱 1.147 0.100 0.300 模型 1 0.146 0.102 0.232 模型 2 0.149 0.102 0.219 模型 3 0.145 0.101 0.218 3.2 外伸板布置形式的影响
图9为不同布置形式肋板模型的平均压力分布. 由图可知:连续和间断布置形式明显影响压力的分布;对于迎风面,肋板的存在使带板模型上的¯CP间断分布,模型4、5的¯CP表面压力分布特征逐渐消失;模型2中的连续肋板明显增加了模型迎风面顶部的¯CP;肋板显著降低整个侧壁面上的¯CP,尤其侧壁面底部,相比参考模型,¯CP最大减小22%.
模型2、4、5的风荷载如表5所示. 由表可知,3个模型的阻力相似,但脉动升力从模型2到模型5逐渐增加,模型5的升力最大,这说明仅在模型角部布置水平板对升力影响较小.
表 5 模型2、4、5的风荷载Table 5. Wind loads of models 2, 4, 5名称 ¯Cd Cd1 Cl1 模型 2 1.149 0.102 0.219 模型 4 1.156 0.106 0.229 模型 5 1.153 0.103 0.261 4. 结 论
本文采用LES方法评估大气边界层流下外伸肋板对高层建筑气动特性的影响. 详细阐述了肋板密度和布置形式对模型流场特征的影响,如近壁流、涡结构和竖向流等,并借助流场信息阐明了壁面风压和整体风荷载的变化,以评估不同肋板的抗风性能,阐明湍流来流下外伸板的流场机理. 得到以下结论:
1) 水平板明显影响尾涡尺寸和尾部鞍点,水平板能明显拉长尾涡. 尾部鞍点高度随着肋板密度的增加而减小,沿建筑物高度的肋板密度过大时,尾部鞍点无明显变化. 外伸板随着水平间距的增大,其对尾涡的影响逐渐变弱.
2) 水平板可以有效抑制建筑立面附近的竖向流动,导致近壁流动形式明显变化. 连续竖直板明显影响近壁流,而间隔肋板对近壁流的影响相对较弱.
3) 流场的变化会改变壁面风压分布. 水平板使迎风面的平均风压¯CP呈“之”字形分布. 水平板可以明显降低侧壁面上平均和脉动风压,最大降幅分别约为20%和17%.
4) 水平板不改变模型的平均阻力¯Cd,但会明显降低模型上的Cl1,最大减幅为27%. 连续肋板相比间隔肋板能更明显的降低模型上的气动力.
5) 连续布置且密度适中的水平肋板对高层建筑周围的流场结构和建筑上的风荷载影响更显著,因此,实际高层建筑外辅助构件的设计和布置时应采用连续水平肋板,肋板布置密度适中即可.
致谢:中国电建集团成都院科技项目(P49321,PA2202)
-
表 1 3D打印材料(PLA)参数
Table 1. 3D printing material (PLA) parameters
参数 线径/mm 密度/(g·cm−3) 熔点/(℃) 拉伸强度/MPa 断裂伸长率/% 弹性模量/MPa 取值 1.75 1.26 176 48 7 3500 表 2 缩尺试验相似数值
Table 2. Similarity parameters of scaled test
物理量 相似关系 相似常数 长度/位移/m CL 8.7 道砟密度/(kg·m−3) Cρ 1.0 内摩擦角 Cφ 1.0 轨枕弹性模量/MPa CE 10.3 轨枕集中荷载/N CF=CECL2 779.6 表 3 模型物理参数与接触参数
Table 3. Physical and contact parameters of model
模型部件 泊松比 剪切模量/
(×1010 Pa)密度/ (kg·m−3) 摩擦系数 道砟 0.30 2.0 2700 0.55 轨枕 0.23 1.5 2800 0.85 钢轨 0.23 1.5 7800 0.85 -
[1] 段海滨. 再生复合轨枕有砟轨道垂向动力特性研究[D]. 成都:西南交通大学,2017. [2] 闫雪,沈毓婷,耿浩,等. 弹性轨枕 + 道砟垫有砟轨道疲劳性能试验研究[J]. 铁道标准设计,2018,62(12): 40-45.YAN Xue, SHEN Yuting, GENG Hao, et al. Experimental research on fatigue behavior of elastic sleeper ballast track with under-ballast mat[J]. Railway Standard Design, 2018, 62(12): 40-45. [3] 刘晓春,陈佳琦,罗延亮,等. 钢管混凝土双块式轨枕无砟轨道结构疲劳性能试验研究[J]. 铁道科学与工程学报,2024,21(1): 106-115.LIU Xiaochun, CHEN Jiaqi, LUO Yanliang, et al. Experimental study on fatigue performance of ballastless track with concrete-filled tube connected double-block sleepers[J]. Journal of Railway Science and Engineering, 2024, 21(1): 106-115. [4] HAYANO K, KOIKE Y, NAKAMURA T, et al. Effects of sleeper shape on lateral resistance of railway ballasted tracks[C]//Advances in Soil Dynamics and Foundation Engineering. Shanghai: American Society of Civil Engineers, 2014: 491-499. [5] 杜香刚. 重载铁路新型高稳定性预应力混凝土轨枕的研发[J]. 铁道建筑,2016,56(3): 113-117. doi: 10.3969/j.issn.1003-1995.2016.03.28DU Xianggang. Research and development of new-type prestressed concrete sleeper with high stability for heavy haul railway[J]. Railway Engineering, 2016, 56(3): 113-117. doi: 10.3969/j.issn.1003-1995.2016.03.28 [6] 井国庆,付豪,贾文利,等. 框架优化型Ⅲc轨枕道床横向阻力试验研究[J]. 西南交通大学学报,2018,53(4): 727-732. doi: 10.3969/j.issn.0258-2724.2018.04.009JING Guoqing, FU Hao, JIA Wenli, et al. Experimental study on lateral resistance of optimized ⅲc sleeper with different frame types[J]. Journal of Southwest Jiaotong University, 2018, 53(4): 727-732. doi: 10.3969/j.issn.0258-2724.2018.04.009 [7] 郭云龙,王新雨,廉栋,等. 摩擦型轨枕道床的横向阻力研究[J]. 西南交通大学学报,2022,57(2): 301-305,368. doi: 10.3969/j.issn.0258-2724.20200464GUO Yunlong, WANG Xinyu, LIAN Dong, et al. Lateral resistance of frictional sleeper ballast bed[J]. Journal of Southwest Jiaotong University, 2022, 57(2): 301-305,368. doi: 10.3969/j.issn.0258-2724.20200464 [8] JING G Q, FU H, AELA P. Lateral displacement of different types of steel sleepers on ballasted track[J]. Construction and Building Materials, 2018, 186: 1268-1275. doi: 10.1016/j.conbuildmat.2018.07.095 [9] ZAKERI J A, TALEBI R. Experimental investigation into the effect of steel sleeper vertical stiffeners on railway track lateral resistance[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2017, 231(1): 104-110. doi: 10.1177/0954409715622500 [10] 井国庆,贾文利,强伟乐,等. 有砟道床梯形轨枕横向阻力试验与构成分析[J]. 西南交通大学学报,2019,54(1): 9-13.JING Guoqing, JIA Wenli, QIANG Weile, et al. Ladder sleeper lateral resistance test and contribution analysis of ballasted track[J]. Journal of Southwest Jiaotong University, 2019, 54(1): 9-13. [11] JING G Q, AELA P, FU H, et al. Numerical and experimental analysis of lateral resistance of biblock sleeper on ballasted tracks[J]. International Journal of Geomechanics, 2020, 20(6): 04020051.1-04020051.17. [12] LIEGNER N. Investigation of the internal forces of the first track constructed with y-shape steel sleepers under operation in Hungary summary of results of research[J]. Periodica Polytechnica Civil Engineering, 2004, 48(1/2): 115-130. [13] 马春生. 高速铁路有砟道床质量评价指标优化方法研究[J]. 铁道标准设计,2016,60(5): 20-24.MA Chunsheng. Research on evaluation index optimization method of high-speed railway ballasted bed[J]. Railway Standard Design, 2016, 60(5): 20-24. [14] 陈成,李建峰,王磊,等. 一种X形抗滑动轨枕:CN209686156U[P]. 2019-11-26. [15] 陈成,王威,芮瑞,等. 装配式可调间距的X形抗滑动轨枕研究[J]. 铁道标准设计,2023,67(2): 49-54.CHEN Cheng, WANG Wei, RUI Rui, et al. Study on X-shaped anti-sliding sleeper with adjustable spacing[J]. Railway Standard Design, 2023, 67(2): 49-54. [16] 陈成,唐子桉,芮瑞,等. 传统轨枕与新型X轨枕下道床劣化及动力特性研究[J]. 北京交通大学学报,2021,45(6): 117-124. doi: 10.11860/j.issn.1673-0291.20210047CHEN Cheng, TANG Zian, RUI Rui, et al. Study on dynamic behavior and ballast degradation using the traditional sleeper and new X-shape sleeper[J]. Journal of Beijing Jiaotong University, 2021, 45(6): 117-124. doi: 10.11860/j.issn.1673-0291.20210047 [17] 万轶,林红松,刘学毅,等. 有碴轨道道床纵向阻力缩尺模型试验研究[J]. 铁道建筑,2006,46(5): 80-81. doi: 10.3969/j.issn.1003-1995.2006.05.030WAN Yi, LIN Hongsong, LIU Xueyi, et al. Experimental study on scale model of longitudinal resistance of ballast track bed[J]. Railway Engineering, 2006, 46(5): 80-81. doi: 10.3969/j.issn.1003-1995.2006.05.030 [18] SHAER A, DUHAMEL D, SAB K, et al. Experimental settlement and dynamic behavior of a portion of ballasted railway track under high speed trains[J]. Journal of Sound and Vibration, 2008, 316(1/2/3/4/5): 211-233. [19] 邹锦华,王荣辉,魏德敏. 车辆—轨道系统缩尺模型振动试验研究[J]. 铁道建筑,2011,51(8): 112-116. doi: 10.3969/j.issn.1003-1995.2011.08.036ZOU Jinhua, WANG Ronghui, WEI Demin. Experimental study on vibration of scale model of vehicle-track system[J]. Railway Engineering, 2011, 51(8): 112-116. doi: 10.3969/j.issn.1003-1995.2011.08.036 [20] 朱建华. 基于相似理论的列车-轨道力学特性分析试验平台研究[D]. 北京:北京交通大学,2022. [21] 中华人民共和国铁道部. 铁路碎石道砟:TB/T 2140—2008[S]. 北京:中国铁道出版社,2008. [22] LIM W L. Mechanics of railway ballast behavior[D]. Nottingham: University of Nottingham, 2004. [23] 李朋. 铁路碎石道床车致垂向振动特性分析[D]. 成都:西南交通大学,2020. [24] ZAKERI J A, HASSANREZAEI H. Experimental investigation on effect of winged sleeper on lateral resistance of ballasted track[J]. Scientia Iranica, 2021, 28(2): 656-665. [25] KOIKE Y, NAKAMURA T, HAYANO K, et al. Numerical method for evaluating the lateral resistance of sleepers in ballasted tracks[J]. Soils and Foundations, 2014, 54(3): 502-514. doi: 10.1016/j.sandf.2014.04.014 [26] LE PEN L M, POWRIE W. Contribution of base, crib, and shoulder ballast to the lateral sliding resistance of railway track: a geotechnical perspective[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2011, 225(2): 113-128. doi: 10.1177/0954409710397094 -