Full Text:   <401>

Summary:  <109>

CLC number: U213.7

On-line Access: 2019-06-04

Received: 2019-03-29

Revision Accepted: 2019-08-17

Crosschecked: 2019-08-27

Cited: 0

Clicked: 465

Citations:  Bibtex RefMan EndNote GB/T7714


Jie-ling Xiao


Gan-zhong Liu


-   Go to

Article info.
Open peer comments

Journal of Zhejiang University SCIENCE A 2019 Vol.20 No.9 P.685-700


Parameters of a discrete element ballasted bed model based on a response surface method

Author(s):  Jie-ling Xiao, Gan-zhong Liu, Jian-xing Liu, Jia-cheng Dai, Hao Liu, Ping Wang

Affiliation(s):  MOE Key Laboratory of High-speed Railway Engineering, Southwest Jiaotong University, Chengdu 610031, China; more

Corresponding email(s):   kevinlau@my.swjtu.edu.cn

Key Words:  Ballasted track, Ballast, Discrete element method, Parameter, Calibration, Response surface method

Jie-ling Xiao, Gan-zhong Liu, Jian-xing Liu, Jia-cheng Dai, Hao Liu, Ping Wang. Parameters of a discrete element ballasted bed model based on a response surface method[J]. Journal of Zhejiang University Science A, 2019, 20(9): 685-700.

@article{title="Parameters of a discrete element ballasted bed model based on a response surface method",
author="Jie-ling Xiao, Gan-zhong Liu, Jian-xing Liu, Jia-cheng Dai, Hao Liu, Ping Wang",
journal="Journal of Zhejiang University Science A",
publisher="Zhejiang University Press & Springer",

%0 Journal Article
%T Parameters of a discrete element ballasted bed model based on a response surface method
%A Jie-ling Xiao
%A Gan-zhong Liu
%A Jian-xing Liu
%A Jia-cheng Dai
%A Hao Liu
%A Ping Wang
%J Journal of Zhejiang University SCIENCE A
%V 20
%N 9
%P 685-700
%@ 1673-565X
%D 2019
%I Zhejiang University Press & Springer
%DOI 10.1631/jzus.A1900133

T1 - Parameters of a discrete element ballasted bed model based on a response surface method
A1 - Jie-ling Xiao
A1 - Gan-zhong Liu
A1 - Jian-xing Liu
A1 - Jia-cheng Dai
A1 - Hao Liu
A1 - Ping Wang
J0 - Journal of Zhejiang University Science A
VL - 20
IS - 9
SP - 685
EP - 700
%@ 1673-565X
Y1 - 2019
PB - Zhejiang University Press & Springer
ER -
DOI - 10.1631/jzus.A1900133

Discrete element simulation on ballasted beds is an important method to study the service characteristics of ballasted track%29&ck%5B%5D=abstract&ck%5B%5D=keyword'>ballasted tracks; an effective simulation should be based on proper ballast parameters. ballast contact parameter, which exhibits a high discreteness affected by factors such as material, shape, and gradation, can effectively be calibrated by an angle of repose test. Based on the testing principles of a multi-parameter response surface method, the Box–Behnken method is adopted to design the angle of repose test under the influence of restitution, static friction, and rolling friction coefficients; laboratory-measured results are combined with the simulation; regression analyzed angle of repose is considered as the goal; parameters optimization and ballasted bed resistance simulations are verified for multiple parameters. The results demonstrate that Chinese special-grade ballasts exhibit an average laboratory-measured angle of repose of (39.78±1.27)°, and the optimal combination of parameters in this discrete element simulation based on the response surface method are as follows: the restitution coefficient is 0.72, the static friction coefficient is 0.56, and the rolling friction coefficient is 0.27. The results of the lateral resistance simulation are in accordance with the laboratory test, indicating that the optimal parameters are usable. The multi-parameter response surface method effectively helps calibrate the parameters of the discrete element simulation on ballasted beds.

The manuscript presents a study on the use of discrete element method (DEM) to study the angle of repose of the special-grade ballast. A discrete element model is developed to examine the lateral resistance forces of the ballast aggregates and sleepers. Regression analysis approach is used to analyse results and the predicted DEM results are also comparable with field measurements. I find this study is interesting and will be useful for readership.


目的:采用响应面方法研究特级道砟休止角的离散元参数仿真试验,建立多次回归模型并对其进行优化,以及对道砟接触参数进行优选. 通过多参数的响应面法为碎石道床离散元参数的快速标定提供有效途径.
创新点:1. 基于响应面法对离散元道砟接触参数进行统计分析,提出道砟最优参数的回归方程及回归曲面. 2. 以道砟堆积体在多个正交平面的均值为道砟休止角,构建基于一致线性描述方式的道砟休止角实测试验与离散元仿真试验. 3. 提出道砟参数的动态标定思想.
方法:1. 通过道砟休止角的室内试验,测量出道砟休止角度的有效均值. 2. 以实测结果为目标,采用响应面法对道砟离散元仿真参数进行优选. 3. 验证所提方法的可行性和有效性.
结论:1. 通过采用休止角的4个正交面取均值的方法对中国特级道砟休止角进行实测所得到的平均值为(39.78±1.27)°. 2. 以实测结果为目标,采用响应面方法对特级道砟离散元仿真参数进行优选,由方差分析可得2个显著的一次项参数(静摩擦系数和滚动摩擦系数)以及多个显著的多次项参数组合; 最优参数组合为:泊松比为0.24,密度为2600 kg/m3,杨氏模量为5.45×1010 Pa,碰撞恢复系数为0.72,静摩擦系数为0.56,滚动摩擦系数为0.27. 3. 利用特级道砟离散元最优参数建立了离散元轨排模型,并进行了轨枕横向阻力试验; 仿真结果与室内试验实测阻力曲线趋势基本一致,表明响应面法获得的接触参数取值可用于碎石道床相关的离散元仿真.

关键词:有砟轨道; 道砟; 离散元; 参数; 标定; 响应面法

Darkslateblue:Affiliate; Royal Blue:Author; Turquoise:Article


[1]Aursudkij B, McDowell GR, Collop AC, 2009. Cyclic loading of railway ballast under triaxial conditions and in a railway test facility. Granular Matter, 11(6):391-401.

[2]Bar-Gera H, 2017. The target parameter of adjusted R-squared in fixed-design experiments. The American Statistician, 71(2):112-119.

[3]Biabani MM, Indraratna B, Ngo NT, 2016. Modelling of geocell-reinforced subballast subjected to cyclic loading. Geotextiles and Geomembranes, 44(4):489-503.

[4]Box GEP, Wilson KB, 1951. On the experimental attainment of optimum conditions. Journal of the Royal Statistical Society: Series B (Methodological), 13(1):1-38.

[5]Chaubey YP, 1993. Resampling-based multiple testing: examples and methods for p-value adjustment. Technometrics, 35(4):450-451.

[6]Chehreghani S, Noaparast M, Rezai B, et al., 2017. Bonded-particle model calibration using response surface methodology. Particuology, 32:141-152.

[7]Chen HH, Bian XC, 2018. Discrete element simulation study of contact pressure distribution between sleeper and ballasts. In: Bian XC, Chen YM, Ye XW (Eds.), Environmental Vibrations and Transportation Geodynamics. Springer, Singapore, p.189-195.

[8]Chen R, Chen JY, Wang P, et al., 2017. Numerical investigation on wheel-turnout rail dynamic interaction excited by wheel diameter difference in high-speed railway. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 18(8):660-676.

[9]Cheng YP, Bolton MD, Nakata Y, 2005. Grain crushing and critical states observed in DEM simulations. In: García-Rojo R, Herrmann HJ, McNamara S (Eds.), Powders and Grains. Taylor & Francis Group, London, UK, p.1393-1397.

[10]Cundall PA, 1971a. A computer model for simulating progressive large scale movements in blocky rock systems. Proceedings of the Symposium of the International Society for Rock Mechanics.

[11]Cundall PA, 1971b. The Measurement and Analysis of Accelerations in Rock Slopes. PhD Thesis, University of London, London, UK.

[12]Cundall PA, Strack ODL, 1980. Discussion: a discrete numerical model for granular assemblies. Géotechnique, 30(3):331-336.

[13]Diyaljee V, 2013. Discussion of “Stress-strain degradation response of railway ballast stabilized with geosynthetics” by Buddhima Indraratna and Sanjay Nimbalkar. Journal of Geotechnical and Geoenvironmental Engineering, 139(12):684-700.

[14]Dong YX, Song ZP, Cui SJ, 2008. Perspectives on the measurement of angle of repose. Journal of China Pharmaceutical University, 39(4):317-320 (in Chinese).

[15]Esveld C, 2001. Modern Railway Track, 2nd Edition. Delft University of Technology, Delft, The Netherlands.

[16]Ferreira SLC, Bruns RE, Ferreira HS, et al., 2007. Box-Behnken design: an alternative for the optimization of analytical methods. Analytica Chimica Acta, 597(2):179-186.

[17]Guo ZG, Chen XL, Liu HF, et al., 2014. Theoretical and experimental investigation on angle of repose of biomass– coal blends. Fuel, 116:131-139.

[18]Hertz H, 1881. On the contact of elastic solids. Journal für die reine und angewandte Mathematik, 92:156-171.

[19]Hochberg Y, 1988. A sharper Bonferroni procedure for multiple tests of significance. Biometrika, 75(4):800-802.

[20]Huang H, Tutumluer E, 2014. Image-aided element shape generation method in discrete-element modeling for railroad ballast. Journal of Materials in Civil Engineering, 26(3):527-535.

[21]Huang H, Tutumluer E, Dombrow W, 2009. Laboratory characterization of fouled railroad ballast behavior. Transportation Research Record: Journal of the Transportation Research Board, 2117(1):93-101.

[22]Ibragimov R, Müller UK, 2010. t-statistic based correlation and heterogeneity robust inference. Journal of Business & Economic Statistics, 28(4):453-468.

[23]Indraratna B, Salim W, 2005. Mechanics of Ballasted Rail Tracks: a Geotechnical Perspective. Taylor & Francis, London, UK.

[24]Indraratna B, Salim W, Rujikiatkamjorn C, 2011. Advanced Rail Geotechnology–Ballasted Track. CRC Press/Balkema, The Netherlands.

[25]Indraratna B, Ngo NT, Rujikiatkamjorn C, et al., 2014. Behavior of fresh and fouled railway ballast subjected to direct shear testing: discrete element simulation. International Journal of Geomechanics, 14(1):34-44.

[26]Khuri AI, Mukhopadhyay S, 2010. Response surface methodology. Wiley Interdisciplinary Reviews: Computational Statistics, 2(2):128-149.

[27]Komar PD, 1978. Angle of repose. In: Sedimentology. Springer, Heidelberg, Germany, p.25-26.

[28]Li YJ, 2005. The Experimental-simulative Study on the Granular Piling Using the Discrete Element Method. MS Thesis, China Agricultural University, Beijing, China (in Chinese).

[29]Liu FY, Zhang J, Li B, et al., 2016. Calibration of parameters of wheat required in discrete element method simulation based on repose angle of particle heap. Transactions of the Chinese Society of Agricultural Engineering, 32(12):247-253 (in Chinese).

[30]Liu GZ, 2018. Study on Stability of Ballast Beds on Bridges Based on the Discrete Element Method. MS Thesis, Southwest Jiaotong University, Chengdu, China (in Chinese).

[31]Liu JX, Xiao JL, Liu H, et al., 2019. Random generation method of ballast 2D topology based on particle characteristics. Construction and Building Materials, 221:762-771.

[32]Lu M, McDowell GR, 2006. Discrete element modelling of ballast abrasion. Géotechnique, 56(9):651-655.

[33]Matuttis HG, Luding S, Herrmann HJ, 2000. Discrete element simulations of dense packings and heaps made of spherical and non-spherical particles. Powder Technology, 109(1-3):278-292.

[34]McDowell GR, Harireche O, Konietzky H, et al., 2006. Discrete element modelling of geogrid-reinforced aggregates. Proceedings of the Institution of Civil Engineers-Geotechnical Engineering, 159(1):35-48.

[35]Mindlin RD, 1949. Compliance of elastic bodies in contact. Journal of Applied Mechanics, 16:259-268.

[36]Mindlin RD, Deresiewicz H, 1953. Elastic spheres in contact under varying oblique forces. Journal of Applied Mechanics, 20:327-344.

[37]MOR (Ministry of Railways of the People’s Republic of China), 2008. Railway ballast, TB/T 2140-2008. National Standards of People’s Republic of China (in Chinese).

[38]Morris AS, Langari R, 2016. Statistical analysis of measurements subject to random errors. In: Morris AS, Langari R (Eds.), Measurement and Instrumentation, 2nd Edition. Academic Press, San Diego, USA, p.75-130.

[39]Ngo NT, Indraratna B, Rujikiatkamjorn C, 2014. DEM simulation of the behaviour of geogrid stabilised ballast fouled with coal. Computers and Geotechnics, 55:224-231.

[40]Ngo NT, Indraratna B, Rujikiatkamjorn C, 2017. A study of the geogrid–subballast interface via experimental evaluation and discrete element modelling. Granular Matter, 19(3):54.

[41]Ngo NT, Indraratna B, Rujikiatkamjorn C, 2017. Simulation ballasted track behavior: numerical treatment and field application. International Journal of Geomechanics, 17(6):04016130.

[42]Qian Y, Mishra D, Tutumluer E, et al., 2015. Characterization of geogrid reinforced ballast behavior at different levels of degradation through triaxial shear strength test and discrete element modeling. Geotextiles and Geomembranes, 43(5):393-402.

[43]Rackl M, Grötsch FE, 2018. 3D scans, angles of repose and bulk densities of 108 bulk material heaps. Scientific Data, 5:180102.

[44]Roessler T, Katterfeld A, 2018. Scaling of the angle of repose test and its influence on the calibration of DEM parameters using upscaled particles. Powder Technology, 330: 58-66.

[45]Sakaguchi H, Ozaki E, Igarashi T, 1993. Plugging of the flow of granular materials during the discharge from a silo. International Journal of Modern Physics B, 07(09n10):1949-1963.

[46]Shorts DC, Feitosa K, 2018. Experimental measurement of the angle of repose of a pile of soft frictionless grains. Granular Matter, 20(1):2.

[47]Tsuji Y, Tanaka T, Ishida T, 1992. Lagrangian numerical simulation of plug flow of cohesionless particles in a horizontal pipe. Powder Technology, 71(3):239-250.

[48]Tutumluer E, Huang H, Bian XC, 2012. Geogrid-aggregate interlock mechanism investigated through aggregate imaging-based discrete element modeling approach. International Journal of Geomechanics, 12(4):391-398.

[49]Wang Y, Yang JH, 2001. Study on rock mass structure and mechanic property of basalt. Chinese Journal of Rock Mechanics and Engineering, 21(9):1307-1310 (in Chinese).

[50]Wu AX, Sun YZ, Liu XP, 2002. Granular Dynamic Theory and Its Applications. Metallurgical Industry Press, Beijing, China (in Chinese).

[51]Xiao JL, Liu H, Xu JM, et al., 2017. Longitudinal resistance performance of granular ballast beds under cyclic symmetric displacement loading. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 18(8):648-659.

[52]Xiao JL, Liu H, Wang P, et al., 2018. Evolution of longitudinal resistance performance of granular ballast track with durable dynamic reciprocated changes. Advances in Materials Science and Engineering, 2018:3189434.

[53]Xu Y, Gao L, Wang H, et al., 2016. Study of fractal method and influence of ballast gradation on ballast bed shear behavior. Journal of the China Railway Society, 38(12):94-101 (in Chinese).

Open peer comments: Debate/Discuss/Question/Opinion


Please provide your name, email address and a comment

Journal of Zhejiang University-SCIENCE, 38 Zheda Road, Hangzhou 310027, China
Tel: +86-571-87952783; E-mail: cjzhang@zju.edu.cn
Copyright © 2000 - Journal of Zhejiang University-SCIENCE