考虑枕下胶垫超弹性本构的弹性长枕轨道动力仿真

第51卷第7期2020年7月
中南大学学报(自然科学版)
Journal of Central South University (Science and Technology)
V ol.51No.7
Jul.2020
考虑枕下胶垫超弹性本构的弹性长枕轨道动力仿真
王启好，蔡小培，常文浩，赵闻强
(北京交通大学土木建筑工程学院，北京，100044)
摘要：为揭示枕下胶垫超弹性本构特征对弹性长枕轨道动力响应的影响，基于Mooney-Rivlin 超弹性本构理论，建立车辆−轨道−隧道−土体空间耦合动力学模型，分析车辆荷载作用下轨道−隧道系统的动力响应，与等效刚度法进行对比，并探讨超弹性本构中枕下胶垫的合理硬度。

研究结果表明：基于超弹性本构仿真得到的轨道各位置位移、加速度数据更真实准确；随着枕下胶垫硬度减小，隧道壁加速度有效值减小，但会放大轨枕垂向位移及下部结构75Hz 频段的振动；综合考虑轨道动力响应、振动传递特性及减振效果，建议弹性长枕轨道枕下胶垫的绍尔硬度取52。

关键词：弹性长枕轨道；枕下胶垫；超弹性本构；动力响应；减振效果；合理硬度中图分类号：U213
文献标志码：A
开放科学(资源服务)标识码(OSID)
文章编号：1672-7207（2020）07-2021-07
Dynamic simulation of long elastic sleeper track based on hyperelastic constitutive models of ruber cushion under sleeper
WANG Qihao,CAI Xiaopei,CHANG Wenhao,ZHAO Wenqiang
(School of Civil Engineering,Beijing Jiaotong University,Beijing 100044,China)
Abstract:In order to reveal the influence of rubber elastic cushion under the sleeper on the dynamic behaviors of elastic long sleeper track,a 3D coupling dynamic model of vehicle −track −tunnel −soil was established based on Mooney-Rivlin hyperelastic constitutive theory.The dynamic responses of track-tunnel system were analyzed and compared with equivalent stiffness method,and the reasonable Shore hardness of elastic cushion was discussed.The results show that the displacement and acceleration data obtained by the hyperelastic constitutive method are closer to the measured results at different locations of metro track.The effective value of tunnel wall acceleration decreases with the decrease of the hardness of elastic cushion.The maximum of sleeper displacement and low-frequency vibration near 75Hz of the structure below the cushion are amplified.Considering the dynamic behaviors,vibration transmission and reduction,the ruber with about 52Shore hardness is recommended as elastic cushion.
Key words:long elastic sleeper track;rubber cushion under sleeper;hyperelastic constitutive law;dynamic responses;vibration reduction;reasonable hardness
DOI:10.11817/j.issn.1672-7207.2020.07.027
收稿日期：2019−09−15；修回日期：2019−12−25
基金项目(Foundation item)：中央高校基本科研业务费专项资金资助项目(2018JBZ003)；国家自然科学基金资助项目(51778050)
(Project(2018JBZ003)supported by Fundamental Research Funds for the Central Universities;Project(51778050)supported by the National Natural Science Foundation of China)
通信作者：蔡小培，博士，教授，从事高速铁路、轨道结构及轨道动力学研究；E-mail:
**************.cn
第51卷中南大学学报(自然科学版)
随着城市轨道交通线网密度的增大，大量线路不可避免地穿过住宅、文教等振动敏感区域。

为解决环境振动问题，减振轨道的应用必不可少。

弹性长枕轨道作为一种中等减振轨道形式，可提供良好的减振效果[1]。

枕下胶垫是影响弹性长枕轨道动力行为和减振效果的关键。

对于橡胶本构和轨道动力仿真，国内外学者已进行了大量研究。

自20世纪50年代MOONEY等[2−3]研究橡胶的弹性大变形特性开始，YEOH[4]比较分析了不同橡胶应变能本构模型的大变形行为；KIM等[5−6]通过试验和数值仿真比较了不同本构模型的仿真效果，探讨了橡胶有限元分析中的误差；杨晓翔[7]对橡胶材料的有限元计算进行了深入研究。

橡胶超弹性本构的应变能方法已十分成熟。

基于多体动力学和有限元方法的轨道动力仿真方法也较为完善。

侯博文等[8]建立了精细化车辆−道岔动力分析模型，研究了车岔系统动力特性；谭诗宇等[9]通过车辆−有砟轨道−桥梁空间耦合动力学模型，分析了桥上道砟垫的减振效果；马晓川等[10]建立了道岔区车辆−轨道系统耦合振动模型，研究了钢轨轧制不平顺
对道岔区垂向轮轨力和轮重减载率的影响；对于橡胶材质的减振垫结构，THOMPSON等[11−12]考虑了刚度的变化，采用随频率变化的刚度进行模拟。

此外，赵才友等[13]对桥上无砟轨道的橡胶减振垫进行了试验研究，定量描述了橡胶垫的减振作用。

可见，既有轨道动力分析中对扣件胶垫、枕下胶垫、道砟垫等橡胶材料的仿真多采用等效刚度法，未考虑橡胶材料的超弹性。

本文作者考虑枕下胶垫橡胶材料的超弹性，建立车辆−弹性长枕轨道−隧道−土体空间耦合动力学模型，分析列车荷载下轨道及隧道的位移、振动加速度等数据，并与等效刚度仿真结果进行对比，最后从轨道动力响应、振动传递特性及减振效果等方面综合考虑，探讨超弹性本构模型中枕下胶垫的合理绍尔硬度。

1橡胶超弹性本构模型及其仿真实现
1.1橡胶超弹性本构模型
橡胶是一种超弹性材料，应力应变关系具有明显的非线性特征[14]，力学性能与金属、混凝土等材料有很大差别。

橡胶本构模型的选取是准确模拟橡胶材料力学行为的关键。

本文在模拟橡胶超弹性时采用基于应变能函数的唯象模型，以应变不变量表征材料特性。

连续介质力学理论中，超弹性材料的应力可通过对应变能函数求偏导得到。

假设橡胶材料为各向同性体，且为不可压缩材料，将应变能密度函数表示成变形张量不变量的函数。

应变能密度函数W由应变偏量能函数f和体积应变能函数g 组成：
W=f(I
1
-3,I
2
-3)+g(J-1)(1)
式中：I
1
和I
2
为Cauchy-Green变形张量不变量；J 为压缩材料系数。

将式(1)泰勒展开可得：
W=
∑i+j=1N C ij(I1-3)i(I2-3)j+∑i=1N1D i(J-1)2i(2)
式中：C
ij
为展开式系数，i，j分别取0或1；D
i
为材料系数；N为多项式阶数，当N=1时，只有线性部分的应变能保留下来，即Mooney-Rivlin模型。

模拟时橡胶材料视为不可压缩，J取1，忽略体积应变能部分。

最终可得不可压缩橡胶材料的Mooney-Rivlin本构模型为
W=C
10
(I
1
-3)+C
01
(I
2
-3)(3) Mooney-Rivlin模型对不同橡胶材料力学行为的模拟适应性较强，适用于中小变形，可用于弹性长枕轨道中枕下胶垫的模拟。

1.2橡胶硬度与仿真实现
Mooney-Rivlin本构模型通过系数C
10
和C
01
确
定。

当橡胶应变在150%以内时，C
01
和C
10
与剪切
模量G、弹性模量E
有以下关系：
G=2(C
01
+C
10
)=E
/3(4)
橡胶弹性模量E
与绍尔硬度H
r
相关，二者关系可通过试验拟合得到。

目前常用的拟合公式有对数和分式[15]形式2种：
lg E
=0.0198H
r
-0.5432(5)
E
=
15.75+2.15H
r
100-H
r
(6)
图1所示为2种拟合公式下橡胶弹性模量随绍尔硬度的变化曲线，仿真模型中取二者平均值。

在有限元模型中将超弹性参数赋予橡胶材料以实现枕下胶垫的超弹性仿真。

2022
第7期
王启好，等：考虑枕下胶垫超弹性本构的弹性长枕轨道动力仿真
2车辆−弹性长枕轨道−隧道−土体
模型
车辆−轨道−隧道−土体空间耦合动力学模型由车辆、轮轨接触、弹性长枕轨道、隧道及土体模型组成，如图2所示。

2.1车辆及轮轨接触模型
车辆模型基于多体动力学建立，由车体、前
后转向架和四轮对组成。

车体、转向架、轮对均视为刚体，考虑车体和转向架的沉浮、横移、点头、摇头、侧滚以及车轮的沉浮、横移、摇头、侧滚共计31个自由度，通过一系二系悬挂连接。

车轮半径为0.42m ，轴距为2.5m ，定距为12.6m ，一系悬挂刚度为1.38kN/mm ，阻尼为45kN·s/m ，二系悬挂刚度为0.35kN/mm ，阻尼为50kN·s/m [16]。

车辆和轨道通过轮轨相互作用耦合，轮轨法向相互作用采用赫兹非线性接触模型，切向相互作用采用非线性摩擦因数模拟[17]。

轨道不平顺采用
北京地铁15号线实测数据。

2.2
弹性长枕轨道、隧道及土体模型
弹性长枕轨道模型基于有限元建立，如图3所示。

钢轨、轨枕、枕下胶垫、道床均采用实体单元模拟，扣件简化为弹簧阻尼单元，考虑纵横垂向刚度和阻尼。

扣件间距为0.6m ，垂向刚度为30kN/mm ，纵横向刚度为40kN/mm ，纵横垂向阻尼为50kN·s/m ；轨枕长为2300mm ，宽为300mm ，高为180mm ；道床下表面为圆弧形，半
径为2.45m ，中间排水沟留空。

橡胶套靴和枕下胶垫密贴，视为一体的枕下胶垫，厚度取15mm 。

道床与隧道、隧道与土体之间视为密贴，不考虑其相对位移。

隧道及土体采用实体单元模拟，土体采用Mohr-Coulomb 本构，考虑其摩擦角和黏聚力。

模型中隧道埋深为12m ，外径为5.80m ，衬砌厚为0.45m 。

土体宽度和深度均取80m ，满足隧道直径的15倍左右和隧道埋深的7倍左右[12]要求。

轨道、隧道及土体模型中材料参数见表1。

2.3
枕下胶垫参数取值
采用超弹性本构模拟枕下胶垫，并与等效刚度法进行对比。

等效刚度仿真时，枕下胶垫刚度取120kN/mm [18]。

由式(4)，(5)和(6)可知，只需确定橡胶硬度及C 01/C 10即可确定唯一的Mooney-Rivlin
本构模型。

超弹性仿真时，为保证橡胶各向
图2车辆−轨道−隧道−土体耦合模型Fig.2
Vehicle −track −tunnel −soil
model
(a)轨道不平顺实测数据；(b)轨道有限元模型
1—垂向不平顺；2—横向不平顺。

图3轨道不平顺实测数据及弹性长枕轨道有限元模型Fig.3
Track irregularity measured data and long elastic sleeper track finite element
model
1—对数拟合；2—分式拟合；3—平均值。

图1
绍尔硬度与弹性模量的拟合关系
Fig.1
Fitting relationship between Shore hardness and
elastic module
2023
第51卷中南大学学报(自然科学版)
刚度与实测值偏差最小，模型中C
01/C
10
取0.05[19]，
枕下胶垫硬度参照现场情况取56，同时选取硬度48，52，60和64进行参数影响分析。

不同绍尔硬度下，枕下胶垫超弹性参数见表2。

3超弹性本构模型对比验证
3.1计算时间分析
在动力学有限元模型仿真中，计算精确度与计算效率相矛盾。

若想得到更精确的结果，需划分更细的网格，会导致计算时间倍增；若想提高计算精度同时不影响计算效率，可在本构算法等方面改进。

超弹性本构模型使用应变能势的响应曲线来评估材料的力学行为，将单位参考体积储存的应变能定义为该点应变的函数。

因此，超弹性本构仿真模型在提高计算精度的同时，仍能保证较高的计算效率。

在相同网格密度、相同分析步长、相同服务器计算能力等条件下，等效刚度法与超弹性本构法耗时基本相同，相差仅2.25%。

3.2仿真结果对比验证
北京地铁15号线现场实测与有限元仿真所得轨道各位置位移、加速度对比分别如图4和图5所示。

从图4可以看出，在2种仿真工况下，轨枕垂向位移均与现场实测相差不大。

超弹性仿真比等效刚度仿真的轨枕垂向位移更接近实测位移。

这是因为超弹性本构可以模拟枕下胶垫刚度的非线性变化，仿真过程中变形更接近实际情况。

从图5可见：超弹性本构仿真比等效刚度仿真有较大优势，各位置加速度及振动衰减量更接近实测数据。

从轨枕到道床，道床到隧道壁，现场实测振动加速度衰减分别为36.980m/s2和1.620m/s2，超弹性仿真时分别为39.632m/s2和1.815m/s2，等效刚度仿真时分别为29.617m/s2和1.004m/s2。

无论是轨枕到道床的衰减，还是道床到隧道壁的衰减，超弹性本构仿真均更贴近实测值，说明超弹性本构工况能更真实地模拟枕下胶垫对振动的衰减效果。

表2橡胶超弹性参数取值Table2Parameters of rubber hyperelasticity
绍尔硬度
48
52
56
60
64C
01
/C
10
0.05
0.05
0.05
0.05
0.05
C
01
/MPa
0.0192
0.0227
0.0269
0.0320
0.0379
C
10
/MPa
0.384
0.454
0.538
0.640
0.758
图5仿真与实测不同位置加速度对比
Fig.5Comparison of accelerations between simulation
and measurement at different
locations
图4仿真与实测轨枕垂向位移对比
Fig.4Comparison of vertical sleeper displacement
between simulation and measurement
表1有限元模型中轨道及隧道材料参数Table1Material parameters of track and tunnel in
finite element model
钢轨类型/ (kg·m−1)
60
衬砌弹性模量/GPa
31.5钢轨弹性
模量/GPa
210.0
土体弹性
模量/GPa
0.3
轨枕弹性
模量/GPa
36.0
土体摩擦
角/(°)
45
道床弹性
模量/GPa
36.0
土体黏聚力/
GPa
1.2
2024
第7期
王启好，等：考虑枕下胶垫超弹性本构的弹性长枕轨道动力仿真
4枕下胶垫合理硬度
弹性长枕轨道通过减小枕下垂向刚度，降低结构的固有频率，使其远离激振频率，减小传递到下部结构的能量。

但是枕下垂向刚度过小会降低轨道承载能力，导致轨道垂向位移过大，进而影响轨道寿命和列车运行安全舒适，因此，有必要对枕下胶垫的合理硬度进行探讨。

4.1
轨道及隧道动力响应有效值分析
当枕下胶垫取不同硬度时，轨枕垂向位移、振动加速度及道床、隧道壁振动加速度有效值如图6所示。

由图6可知：当枕下胶垫硬度为48时，轨枕垂向位移有效值达到最大，为0.205mm ；在枕下胶垫硬度取64时，道床和隧道壁振动加速度达到最大，分别为0.653m/s 2和0.173m/s 2；随着枕下胶垫硬度的增加，轨枕垂向位移有效值逐渐减小，道床和隧道壁加速度有效值逐渐增加，且在硬度48~64范围内呈线性变化；轨枕振动加速度有效值先增加后保持平稳。

4.2
轨道及隧道加速度功率谱密度分析
为揭示振动在轨枕、道床、隧道壁间的传递规律，进一步分析了轨道及隧道的加速度功率谱密度，为在同一坐标系中比较，将道床功率谱放大100倍，隧道壁功率谱放大1000倍。

图7所示
为枕下胶垫硬度56时轨枕、道床、隧道壁的加速
度功率谱密度，图8~10所示为不同枕下胶垫硬度下轨枕、隧道壁的加速度功率谱密度及振级。

由图7可以看出，从轨枕到道床再到隧道壁，440Hz 频段加速度高频分量衰减至0；75Hz 以及145Hz 频段分量衰减相对较慢；从轨枕到道床再到隧道壁，75Hz 频段振动占比逐渐增大；215Hz 频段振动一直占主导地位。

75Hz 处振动衰减最
慢，可以预见振动继续在土层中传播时，该频率振动将逐渐占据主导地位。

图8所示为不同枕下胶垫硬度下轨枕加速度功率谱密度。

从图8可以看出，峰值大多出现在215Hz 和440Hz 频段；在枕下胶垫硬度为64时，还会在138Hz 处出现较大峰值；在215Hz 频段，随着枕下胶垫硬度的增大，轨枕加速度功率谱密度峰值先增大后减小；在440Hz 频段，随着枕下胶垫硬度增大，轨枕加速度功率谱密度峰值减小，表明轨枕振动与轨道结构的模态振型有关[20]；当枕下胶垫硬度为64时，易引起轨枕138Hz 频段的共振；当枕下胶垫硬度为56~60时，易引起轨枕215Hz 频段的共振。

图9所示为不同枕下胶垫硬度下隧道壁加速度功率谱密度。

从图9可知：频率为75Hz 和225Hz 处均出现明显峰值。

频率225Hz 峰值随着枕下胶垫硬度的增大而增大，说明较硬的枕下胶垫会放大隧道壁的高频振动；在75Hz 频段，峰值随着枕
下胶垫硬度的减小而增大，说明当枕下胶垫硬度
1—轨枕振动加速度；2—轨枕垂向位移；3—道床振动加速度；4—隧道壁1.5m 加速度。

图6
动力响应与绍尔硬度的关系
Fig.6
Relationship between dynamic response and Shore
hardness
1—轨枕；2—道床扩大100倍的结果；3—隧道壁扩大1000倍的结果。

图7
不同位置加速度功率谱密度对比Fig.7
Comparison of acceleration power spectral density
(PSD)for different locations
2025
第51卷
中南大学学报(自然科学版)较小时，低频振动会放大；随着枕下胶垫硬度的增大，在145Hz 频段，加速度功率谱密度先增大后减小。

图10所示为不同枕下胶垫硬度下轨枕及隧道壁的加速度级。

从图10可知：当枕下胶垫硬度从48增大到56时，轨枕加速度级增大3.01dB ，之后保持平稳；当枕下胶垫硬度大于56时，易引起轨枕共振，因此，设计中应避开共振频率，枕下胶垫硬度小于56。

隧道壁加速度级与枕下胶垫硬度呈线性关系，枕下胶垫硬度每增加4，隧道壁加速度级增大约1dB 。

综上所述，为保证较小的轨枕位移和加速度、
较好的减振效果，同时考虑振动的传播特性，建
议枕下胶垫绍尔硬度取52。

5结论
1)与现有的等效刚度仿真相比，超弹性本构仿真能更准确地模拟枕下胶垫的力学行为，得到的位移、加速度等更贴近轨道各位置实测结果。

2)随枕下胶垫硬度增加，轨枕垂向位移有效值线性减小，道床和隧道壁加速度有效值线性增加；在枕下胶垫硬度为48~56时，轨枕加速度增加迅速，硬度达到56以后保持平稳。

3)轨道及隧道加速度功率谱密度主频均在
225Hz 频段，75Hz 频段振动衰减最慢；枕下胶垫绍尔硬度为48时会放大隧道壁75Hz 频段振动，硬度为64时会放大隧道壁225Hz 频段振动。

4)为达到较好的减振效果，同时保证轨道动力响应和隧道壁低频振动分量不过大，建议弹性长枕轨道枕下胶垫绍尔硬度取为52。

参考文献：
[1]
CAI Xiaopei,LI Dacheng,ZHANG Yanrong,et al.Experimental study on the vibration control effect of long elastic sleeper track in subways[J].Shock and Vibration,2018,526:1−13.[2]MOONEY M.A theory of large elastic deformation[J].Journal of Applied Physics,1940,11(9):582−592.
[3]
RIVLIN R S,SAUNDERS D rge elastic
deformations
1—轨枕；2—隧道壁。

图10不同枕下胶垫硬度下轨枕及隧道壁加速度级Fig.10Acceleration level of sleeper and tunnel wall
under different hardness of elastic
cushion
绍尔硬度：1—48；2—52；3—56；4—60；5—64。

图8
不同枕下胶垫硬度下轨枕加速度功率谱密度对比Fig.8
Comparison of sleeper acceleration PSD under different hardness of elastic
cushion
绍尔硬度：1—48；2—52；3—56；4—60；5—64。

图9
不同枕下胶垫硬度下隧道壁加速度功率谱
密度对比
Fig.9Comparison of tunnel wall acceleration PSD under
different hardness of elastic cushion
2026
第7期王启好，等：考虑枕下胶垫超弹性本构的弹性长枕轨道动力仿真
of isotropic materials VII:experiments on the deformation of rubber[J].Philosophical Transactions of the Royal Society of London Series A:Mathematical and Physical Sciences, 1951,243:251−288.
[4]YEOH O H.Some forms of the strain energy function for
rubber[J].Rubber Chemistry and Technology,1993,66(5): 754−771.
[5]KIM B,LEE S B,LEE J,et al.A comparison among Neo-
Hookean model,Mooney-Rivlin model,and Ogden model for chloroprene rubber[J].International Journal of Precision Engineering and Manufacturing,2012,13(5):759−764. [6]CHARLTON D J,YANG J,TEH K K.A review of methods
to characterize rubber elastic behavior for use in finite element analysis[J].Rubber Chemistry&Technology,1994, 67(3):481−503.
[7]杨晓翔.非线性橡胶材料的有限单元法[M].北京:石油工
业出版社,1999:70−82.
YANG Xiaoxiang.Finite element method for nonlinear rubber materials[M].Beijing:Petroleum Industry Press, 1999:70−82.
[8]侯博文,高亮,刘启宾.重载车辆−道岔耦合动力特性及岔
区加强研究[J].西南交通大学学报,2015,50(4):604−609.
HOU Bowen,GAO Liang,LIU Qibin.Vehicle-turnout coupled dynamics and improvement measures for heavy haul lines[J].Journal of Southwest Jiaotong University,2015,50
(4):604−609.
[9]谭诗宇,蔡小培,崔日新,等.环境敏感区桥上有砟轨道铺
设道砟垫的减振效果[J].振动与冲击,2017,36(10):38−44.
TAN Shiyu,CAI Xiaopei,CUI Rixin,et al.Vibration reduction effect of laying ballast mats under ballast tracks on bridge in environmental sensitive areas[J].Journal of Vibration and Shock,2017,36(10):38−44.
[10]马晓川,王平,徐金辉,等.钢轨轧制不平顺对车岔耦合系
统垂向动力特性的影响[J].中南大学学报(自然科学版), 2017,48(7):1942−1950.
MA Xiaochuan,WANG Ping,XU Jinhui,et al.Effect of rail straightening irregularity on vertical dynamic characteristics of vehicle-turnout coupling system[J].Journal of Central South University(Science and Technology),2017,48(7): 1942−1950.
[11]THOMPSON D J,V AN VLIET W J,VERHEIJ J W.
Developments of the indirect method for measuring the high frequency dynamic stiffness of resilient elements[J].Journal of Sound and Vibration,1998,213(1):169−188.
[12]韦凯,杨帆,王平.扣件胶垫刚度的频变性对地铁隧道环境
振动的影响[J].铁道学报,2015,37(4):80−86.
WEI Kai,YANG Fan,WANG Ping.Influence of frequency-
dependent stiffness of rail pads on environment vibration induced by subway in tunnel[J].Journal of the China Railway Society,2015,37(4):80−86.
[13]赵才友,王平.桥上无砟轨道橡胶减振垫减振性能试验研
究[J].中国铁道科学,2013,34(4):8−13.
ZHAO Caiyou,WANG Ping.Experimental study on the vibration damping performance of rubber absorbers for ballastless tracks on viaduct[J].China Railway Science, 2013,34(4):8−13.
[14]刘艳,罗雁云.钢轨扣件减振橡胶动态特性分析[J].中南大
学学报(自然科学版),2011,42(9):2875−2881.
LIU Yan,LUO Yanyun.Dynamic characteristic analysis of rubber absorber in rail fastening[J].Journal of Central South University(Science and Technology),2011,42(9):2875−2881.
[15]FREAKLEY P K,PAYNE A R.Theory and practice of
engineering with rubber[M].London:Applied Science Publishing Ltd.,1978:443−484.
[16]CAI Xiaopei,WANG Qihao,LI Dacheng,et al.Theoretical
study of long elastic sleeper track to reduce vibrations induced by subway trains[J].Proceedings of the Institution of Mechanical Engineers,Part F:Journal of Rail and Rapid Transit,2020,234(5):538−549.
[17]ZHAI Wanming,SUN Xiang.A detailed model for
investigating vertical interaction between railway vehicle and track[J].Vehicle System Dynamics,1994,23(S1):603−615.
[18]赵信洋.弹性长枕式无砟轨道枕长与支撑长度比例分析
[J].路基工程,2010(4):202−204.
ZHAO Xinyang.Analysis on sleeper length and supporting length ratio in long elastic sleeper ballastless track[J].
Subgrade Engineering,2010(4):202−204.
[19]郑明军,王文静,陈政南,等.橡胶Mooney-Rivlin模型力学
性能常数的确定[J].橡胶工业,2003,50(8):462−465.
ZHENG Mingjun,WANG Wenjing,CHEN Zhengnan,et al.
Determination for mechanical constants of rubber Mooney-Rivlin model[J].China Rubber Industry,2003,50(8):462−465.
[20]付娜,李成辉,杨荣山,等.高速铁路减振双块式无砟轨道
的减振性能[J].中南大学学报(自然科学版),2018,49(2): 502−509.
FU Na,LI Chenghui,YANG Rongshan,et al.Damping performance of double-block ballastless damping track in high-speed railway[J].Journal of Central South University (Science and Technology),2018,49(2):502−509.
(编辑赵俊)
2027。