Abstract

In the last two decades there has been a great development of high-speed train lines. This advance has led to more demanding loads in railway infrastructures and the appearance of a new problem called ballast flight that happens when some stones rise when the train passes. For these reasons, the development of an application that allows the numerical modelling of the ballast superstructure under different stresses can be very useful, as it will enable design optimization. The DEM is being considered an effective and powerful method for the calculation of engineering problems with granular and discontinuous materials. Railroad ballast layer consists of discrete aggregate particles, so that DEM is one of the most suitable ways to simulate the behavior of particulate ballast material. However, the computational cost of contact calculation between irregular particles is high and limits the calculation capability. From the point of view of micro-scale analysis, it is essential to represent the exact geometry of the particle. On the other hand, if the interest lies in the behavior of the granular material as a whole, the geometry is not a determining factor. Besides that, setting up a simulation of granular material taking care of the exact geometry of each particle will not be efficient. Current work presents the methodology followed to achieve accurate results in the calculation of railway ballast behaviour using DEM and spherical particles. The use of spherical particles reduces the computational cost and makes the simulation set up efficient. Validation results for the calculation of the lateral resistance force against a sleeper moving inside a ballast bed are presented. Regarding ballast flight problem, some high speed ballast collision calculations have also been performed.

Abstract

The development of high-speed train lines has increased during the last twenty years, leading to more demanding loads in railway infrastructures. For these reasons, a discrete element model of granular material was carried out using 3D spherical particles with a rolling resistance, in order to consider non-sphericity of ballast stones. Other discrete element model representing granular materials as sphere aggregates (sphere clusters) has been also implemented. Current work presents the methodology followed to develop the discrete element model able to calculate railway ballast interaction. Interaction between discrete and finite elements is another key point of the calculations that is addressed. This document displays some results evaluating the influence of material parameters and geometric representation in the simulation of railway ballast and presents some laboratory tests calculated with the numerical application developed.

Abstract

The development of high-speed train lines has increased significantly during the last twenty-five years, leading to more demanding loads in railway infrastructures. Most of these infrastructures were constructed using railway ballast, which is a layer of granular material placed under the sleepers whose roles are: resisting to vertical and horizontal loads and facing climate action. Moreover, the Discrete Element Method was found to be an effective numerical method for the calculation of engineering problems involving granular materials. For these reasons, the main objective of the work is the development of a numerical modelling tool based on the Discrete Element Method which allows the users to understand better the mechanical behaviour of railway ballast. The first task was the review of the specifications that ballast material must meet. Then, the features of the available Discrete Elements code, called ''DEMPack'', were analysed. After those revisions, it was found that the code needed some improvement in order to reproduce correctly and efficiently the behaviour of railway ballast. The main deficiencies identified in the numerical code were related to the contact between discrete element particles and planar boundaries and to the geometrical representation of such a irregular material as ballast. Contact interactions between rigid boundaries and Discrete Elements are treated using a new methodology called the Double Hierarchy method. This new algorithm is based on characterising contacts between rigid parts (meshed with a Finite Element-like discretisation) and spherical Discrete Elements. The procedure is described in the course of the work. Moreover, the method validation and the assessment of its limitations are also displayed. The representation of irregular particles using the Discrete Element Method is a very challenging issue, leading to different geometrical approaches. In this work, a deep revision of those approaches was performed. Finally, the most appropriate methods were chosen: spheres with rolling friction and clusters of spheres. The main advantage of the use of spheres is their low computational cost, while clusters of spheres stand out for their geometrical versatility. Some improvements were developed for describing the movement of each kind of particles, specifically, the imposition of the rolling friction and the integration of the rotation of clusters of spheres. In the course of this work the way to fill volumes with particles (spheres or clusters) was also analysed. The aim is to control properly the initial granulometry and compactness of the samples used in the calculations. After checking the correctness of the numerical code with simplified benchmarks, some laboratory tests for evaluating railway ballast behaviour were computed. The aim was to calibrate the ballast material properties and validate the code for the representation of railway ballast. Once the material properties were calibrated, some examples of a real train passing through a railway ballast track were reproduced numerically. This calculations allowed to prove the possibilities of the implemented tool.

This publication is a revised version of the text of the PhD Thesis “Numerical analysis of railway ballast behaviour using the Discrete Element Method” of Joaquín Irazábal, presented at the Technical University of Catalonia on October 6th 2017

Abstract

The development of high-speed train lines has increased during the last twenty years, leading to more demanding loads in railway infrastructures. For these reasons, the implementation of a numerical tool for the calculation of railway ballast behaviour has been found useful, as it will enables design optimization.

Regarding the numerical method, the DEM is considered effective and powerful for the calculation of engineering problems with granular and discontinuous materials. Due to the fact that railroad ballast layer consists of discrete aggregate particles, the DEM is considered suitable for the simulation of particulate ballast material. However, the computational cost of contact calculation between irregular particles is high and limits the calculation capability.

From the point of view of micro-scale analysis, it is essential to represent the exact geometry of the particle. On the other hand, if the interest lies in the behaviour of the granular material as a whole, the geometry is not a determining factor. Besides that, setting up a simulation of granular material taking care of the exact geometry of each particle will not be efficient.

Current work presents different geometrical approaches for the representition of ballast stones: spheric particles with rolling friction, sphere clusters, polyhedrons and superquadrics; showing their advantages and drawbacks.

Finally, some simulation results, using spheric particles and sphere clusters, are displayed in order to evaluate
the convenience or not of using more accurate and computational demanding geometries in each case.

Abstract

Railway ballast is a layer of granular material that resists to vertical and horizontal loads, produced by the passing train over the rail. The calculation of this kind of complex geomechanic problems has been traditionally addressed using refined constitutive models, based in continuum assumptions. Although these models may be suitable in the evaluation of the critical state of soils, or in the calculation of bulk material masses flowing, they are not appropriate to represent the local discontinuities of granular materials, which induce special features such as anisotropy or instabilities.

The Discrete Element Method (DEM) is an alternative approach that considers the discontinuous nature of granular materials, which has proven to be a very useful tool to obtain complete qualitative information on calculations of groups of particles. However, the computational cost of contact evaluation between Discrete Elements (DEs) is high and limits the calculation capability. In this regard, it should be noted that particle shape greatly affects contact calculation computational cost, being spherical DEs the less computational demanding type of particles.

From the point of view of micro-scale analysis, it is essential to represent the exact geometry of the grains. By contrast, if the interest lies in the behaviour of the granular material as a whole, particles geometry is not a determining factor. Therefore, for efficiency purposes, a trade-off between particle shape accuracy and computational cost needs to be achieved.

In this work, different approaches to represent ballast stones are assessed: spheres with rolling friction, sphere clusters, polyhedrons and superquadrics. The first two were chosen for further analysis.

Rolling friction allows avoiding excessive rotation when irregular shaped materials are simulated as spherical particles. This work presents a new insight for its application called the Bounded Rolling Friction model.

Regarding sphere clusters, there is a key point in the friction between elements. As they reproduce irregular particles using clumps of spheres rigidly joined, the cavities between those spheres introduce interlocks that increase friction. To overcome this drawback, a new contact model is proposed.

Finally, results of the application of both approaches are displayed, and conclusions are drawn as regards the convenience of using more accurate and computational demanding geometries.

Abstract

The development of high-speed train lines has increased significantly during the last decades leading to more demanding loads in railway infrastructures. Most of these infrastructures were constructed using railway ballast, whose main roles are resisting to vertical and horizontal loads
and facing climate action. Moreover, new challenges are arising in the railway industry, such as the development of high-speed train lines in locations with extreme weather. For these reasons, the implementation of a numerical code able to represent ballast behaviour, including its interaction with other structures, has become very attractive.

Among a wide range of numerical methods, the Discrete Element Method (DEM) was found to be effective for the calculation of engineering problems with granular materials. This approach considers the discontinuous nature of these materials and has proven to be a very useful tool to obtain complete qualitative information on calculations of groups of particles.

The code used in this work is developed within DEMPack, a specific software tool for modelling physical problems using the DEM. The computer program is adapted to meet the needs for representing the behaviour of railway ballast.

Abstract

Rail transport, both for people and goods, is becoming increasingly significant all over the world, which is reflected in the great growth of conventional and high-speed train lines. Most of these infrastructures are built with railway ballast, a granular material whose main functions are to resist vertical and horizontal loads and to face climatic actions.

The growing popularity of these infrastructures has led to the development of numerical models to evaluate their performance. Among a wide range of numerical methods, the Discrete Element Method (DEM) was found to be effective for evaluating the performance of granular materials. This approach considers their discontinuous nature and has proven to be a useful tool to determine the dynamic behaviour of groups of particles. Moreover, the DEM is also used to compute the behaviour of continuum materials. In this work, rails and bearing plates are characterised in the calculations using this methodology, called the bonded DEM. It is a modification of the classical DEM which assumes that bonds exist between particles, resisting their separation.

The code used is developed within DEMPack, a specific software tool for modelling physical problems using the DEM. Currently, DEMPack allows the use of two different types of geometry: spheres with rolling friction and clusters of spheres. A previous analysis showed that spheres are more effective for studying the macroscopic behaviour of the ballast layer, while clusters are necessary for small-scale tests involving highly compacted particles since their results are greatly influenced by particles and contacts distribution.

After calibrating the code, full-scale tests were performed applying the load of a high-speed train on a railway track section in different situations. Considering the amount of material (about 260,000 particles) and that the aim is to evaluate the deflection of the rails, the calculations are carried out using spheres.

The numerical results correctly capture the effect on the deflection of the rails. It can be concluded that the DEM increases the possibilities for analysing innovative solutions since real case-scenarios can be studied with enough accuracy and feasible time.

Abstract

Conventional and high-speed train lines are being constructed all over the world with the objective of improving the mobility of both people and goods. Most of these infrastructures are built with railway ballast, a granular material whose main functions are resisting train loads and facing climate actions. The growth in popularity of railway infrastructures has led to an increasing interest in the development of numerical models to evaluate their performance and improve their maintenance. To this respect, the Discrete Element Method (DEM) is an approach that considers the discontinuous nature of granular materials, such as railway ballast. Moreover, it can also be used to compute the behaviour of continuum materials, such as rails and bearing plates, applying the so-called bonded DEM. The code used is developed within DEMPack, a specific software tool for modelling physical problems using the DEM. After evaluating the different geometrical alternatives for representing railway ballast and considering the high amount of material involved in the full-scale tests (more than 130,000 particles) the calculations were carried out using the most efficient option, spheres with rolling friction. This geometrical simplification is not suitable for small-scale tests, however, previous analysis showed that it is accurate enough for reproducing the macroscopic behaviour of the ballast layer. The numerical results correctly capture the effect of changing several parameters such as ballast compaction, inter-particle friction or grain size. It can be concluded that the DEM increases the possibilities for analysing innovative solutions aiming to improve ballasted tracks design, maintenance and performance, since real case-scenarios can be studied with enough accuracy and feasible time.

Abstract