DEM modelling of surface indentations caused by granular materials: application to wheel-rail sanding.

The presented surface indentation model is one step towards building a DEM model for wheel-rail sanding. In railways, so-called low-adhesion conditions can cause problems in traction and braking, and sanding is used to overcome this problem. Sand grains are blasted towards wheel-rail contact, fracture repeatedly as they enter the nip and are drawn into the contact and then increase adhesion. Research on this topic has mostly been experimental, but focussed on adhesion enhancement measurement. Thus, physical mechanisms increasing the adhesion are not well understood. Previous works involved experiments and DEM modelling of single sand grain crushing tests under realistic wheel-rail contact pressures of 900 MPa, focusing on sand fragment spread and formation of clusters of solidified fragments. In the experiments, indents in the compressing steel plates were also observed, which are also observed on wheel and rail surfaces in railway operation. These are now modelled by adapting an existing surface indentation model from literature to the case of surface indentations caused by granular materials. Two test cases are studied, and experimental spherical indentation tests for model parametrisation are presented. In a proof of concept, the mentioned single sand grain crushing tests under 900 MPa pressure are simulated including the surface indentation model. This work contributes to DEM modelling of wheel-rail sanding, which is believed to be a good approach to deepen the understanding of adhesion increasing mechanisms under sanded conditions.

Efficient DEM simulations of railway ballast using simple particle shapes.

For complex shaped materials, computational efficiency and accuracy of DEM models are usually opposing requirements. In the literature, DEM models of railway ballast often use very complex and computationally demanding particle shapes in combination with very simple contact laws. In contrast, this study suggests efficient DEM models for railway ballast using simple particle shapes together with a contact law including more physical effects. In previous works of the authors, shape descriptors, calculated in a shape analysis of two types of ballast, were used to construct simple particle shapes (clumps of three spheres). Using such a shape in DEM simulations of compression and direct shear tests, accurate results were achieved only when the contact law included additional physical effects e.g. edge breakage. A parametrisation strategy was developed for this contact law comparing DEM simulations with the measurements. Now, all the constructed simple particle shapes are parametrised allowing to study their suitability and relating their shape descriptors to those of railway ballast. The most suitable particle shapes consist of non-overlapping spheres, thus have a high interlocking potential, and have lowest sphericity and highest convexity values. In a micromechanical analysis of the four best performing shapes, three shapes show similar behaviour on the bulk and the micro-scale, while one shape differs clearly on the micro-scale. This analysis shows, which shapes can be expected to produce similar results in DEM simulations of other tests/load cases. The presented approach is a step towards both efficient and accurate DEM modelling of railway ballast.

DEM modelling of railway ballast using the Conical Damage Model: a comprehensive parametrisation strategy.

Despite ongoing research, the parametrisation of a DEM model is a challenging task, as it depends strongly on the particle shape representation used, particle-particle contact law and the simulated applications: for railway ballast e.g. lab tests or track conditions. The authors previously modelled railway ballast with a DEM model using a simple particle shape. The DEM model was parametrised, by trial-and-error, to compression and direct shear test results. A good agreement between DEM model and experimental results was achieved only when the Conical Damage Model (CDM) was used as the contact law. Compared to the well-known linear-spring Cundall-Strack law or the Hertz-Mindlin law, this contact law takes into account additional physical effects (e.g. edge breakage) occurring in the experiment. Little is known on the influence of the CDM model parameters on the simulation results or on possible parameter ambiguities. This lack of knowledge hinders a reliable and efficient parametrisation of DEM models using different particle shapes. Both points are addressed in this work in detail by investigating a DEM model for railway ballast using one simple particle shape. Suggestions for a parametrisation strategy of reduced computational effort are formulated and tested using a second particle shape. In future works, the newly presented parametrisation strategy can help to calibrate different DEM models and to study the influence of particle shape.

Search, reuse and sharing of research data in materials science and engineering-A qualitative interview study.

Open research data practices are a relatively new, thus still evolving part of scientific work, and their usage varies strongly within different scientific domains. In the literature, the investigation of open research data practices covers the whole range of big empirical studies covering multiple scientific domains to smaller, in depth studies analysing a single field of research. Despite the richness of literature on this topic, there is still a lack of knowledge on the (open) research data awareness and practices in materials science and engineering. While most current studies focus only on some aspects of open research data practices, we aim for a comprehensive understanding of all practices with respect to the considered scientific domain. Hence this study aims at 1) drawing the whole picture of search, reuse and sharing of research data 2) while focusing on materials science and engineering. The chosen approach allows to explore the connections between different aspects of open research data practices, e.g. between data sharing and data search. In depth interviews with 13 researchers in this field were conducted, transcribed verbatim, coded and analysed using content analysis. The main findings characterised research data in materials science and engineering as extremely diverse, often generated for a very specific research focus and needing a precise description of the data and the complete generation process for possible reuse. Results on research data search and reuse showed that the interviewees intended to reuse data but were mostly unfamiliar with (yet interested in) modern methods as dataset search engines, data journals or searching public repositories. Current research data sharing is not open, but bilaterally and usually encouraged by supervisors or employers. Project funding does affect data sharing in two ways: some researchers argue to share their data openly due to their funding agency's policy, while others face legal restrictions for sharing as their projects are partly funded by industry. The time needed for a precise description of the data and their generation process is named as biggest obstacle for data sharing. From these findings, a precise set of actions is derived suitable to support Open Data, involving training for researchers and introducing rewards for data sharing on the level of universities and funding bodies.

Simple particle shapes for DEM simulations of railway ballast: influence of shape descriptors on packing behaviour.

ABSTRACT: In any DEM simulation, the chosen particle shape will greatly influence the simulated material behaviour. For a specific material, e.g. railway ballast, it remains an open question how to model the particle shape, such that DEM simulations are computationally efficient and simulation results are in good accordance with measurements. While DEM shape modelling for railway ballast is well addressed in the literature, approaches mainly aim at approximating the stones' actual shape, resulting in rather complex and thus inefficient particle shapes. In contrast, very simple DEM shapes will be constructed, clumps of three spheres, which aim to approximate shape descriptors of the considered ballast material. In DEM simulations of the packing behaviour, a set of clump shapes is identified, which can pack at porosities observed at track sites, as well as in lab tests. The relation between particle shape (descriptors) and obtained packing (characteristic) is investigated in a correlation analysis. The simulated packing's porosity is strongly correlated to four shape descriptors, which are also strongly correlated among each other. Thus, to derive simple shape models of a given particle shape, matching one of these shape descriptors, might be a good first step to bring simulated porosities closer to measured ones. The conducted correlation analysis also shows that packing's coordination number and isotropic fabric are correlated to more shape descriptors, making it more difficult to estimate the effect of particle shape on these quantities.

Shape analysis of railway ballast stones: curvature-based calculation of particle angularity.

Particle shape analysis is conducted, to compare two types of railway ballast: Calcite and Kieselkalk. Focus lies on the characterisation of particle angularity using 3D scanner data. In the literature, angularity is often characterised using 2D data, as these types of data are easier to collect. 3D scanner data contain a vast amount of information (e.g. curvatures) which can be used for shape analysis and angularity characterisation. Literature approaches that use 3D data are often not thoroughly tested, due to a lack of test cases. In this work, two new curvature-based angularity indices are introduced and compared to one from the literature. Analytical test bodies with shapes ranging from spherical towards cubic are used for a first plausibility test. Then, 3D scans of ballast stones are compared to artificially rounded meshes. Only one out of three evaluated angularity indices seem to be suited to characterise angularity correctly in all of the above tests: the newly introduced scaled Willmore energy. A complete shape analysis of the scanned ballast stones is conducted and no difference between the two types of ballast can be seen regarding form, angularity, roughness, sphericity or convexity index. These findings of shape analysis are set in the context of previous works, where experimental results and DEM simulations of uniaxial compression tests and direct shear tests were presented for the same ballast types.

Micro-mechanical investigation of railway ballast behavior under cyclic loading in a box test using DEM: effects of elastic layers and ballast types.

ABSTRACT: Ballasted tracks are the commonly used railway track systems with constant demands for reducing maintenance cost and improved performance. Elastic layers are increasingly used for improving ballasted tracks. In order to better understand the effects of elastic layers, physical understanding at the ballast particle level is crucial. Here, discrete element method (DEM) is used to investigate the effects of elastic layers - under sleeper pad ( USP ) at the sleeper/ballast interface and under ballast mat ( UBM ) at the ballast/bottom interface - on micro-mechanical behavior of railway ballast. In the DEM model, the Conical Damage Model (CDM) is used for contact modelling. This model was calibrated in Suhr et al. (Granul Matter 20(4):70, 2018) for the simulation of two different types of ballast. The CDM model accounts for particle edge breakage, which is an important phenomenon especially at the early stage of a tamping cycle, and thus essential, when investigating the impact of elastic layers in the ballast bed. DEM results confirm that during cyclic loading, USP reduces the edge breakage at the sleeper/ballast interface. On the other hand, UBM shows higher particle movement throughout the ballast bed. Both the edge breakage and particle movement in the ballast bed are found to influence the sleeper settlement. Micro-mechanical investigations show that the force chain in deeper regions of the ballast bed is less affected by USP for the two types of ballast. Conversely, dense lateral forces near to the box bottom were seen with UBM . The findings are in good (qualitative) agreement with the experimental observations. Thus, DEM simulations can aid to better understand the micro-macro phenomena for railway ballast. This can help to improve the track components and track design based on simulation models taking into account the physical behavior of ballast.

Comparison of two different types of railway ballast in compression and direct shear tests: experimental results and DEM model validation.

Railway ballast is an angular and coarse material, which demands careful DEM modelling and validation. Particle shape is often modelled in high accuracy, thus leading to computational expensive DEM models. Whether this effort will increase the DEM model's overall prediction quality will also vitally depend on the used contact law and the validation process. In general, a DEM model validated using different types of principal experiments can be considered more trustworthy in simulating other load cases. Here, two types of railway ballast are compared and DEM model validation is conducted. Calcite and Kieselkalk are investigated under compression and direct shear test. All experimental data will be made openly accessible to promote further research on this topic. In the experiments, the behaviour of Calcite and Kieselkalk is surprisingly similar in the direct shear test, while clear differences can be seen in the stiffnesses in the compression test. In DEM modelling, simple particle shapes are combined with the Conical Damage Model contact law. For each type of ballast, one set of parameters is found, such that simulation and experimental results are in good accordance. A comparison with the simplified Hertz-Mindlin contact law shows several drawbacks of this model. First, the model cannot be calibrated to meet both compression and shear test results. Second, the similar behaviour in shear testing but differences in compression cannot be reproduced using the Hertz-Mindlin model. For these reasons, the CDM model is considered the better choice for the simulation of railway ballast, if simple particle shapes are used.

Parametrisation of a DEM model for railway ballast under different load cases.

The prediction quality of discrete element method (DEM) models for railway ballast can be expected to depend on three points: the geometry representation of the single particles, the used contact models and the parametrisation using principal experiments. This works aims at a balanced approach, where none of the points is addressed with excessive depth. In a first step, a simple geometry representation is chosen and the simplified Hertz-Mindlin contact model is used. When experimental data of cyclic compression tests and monotonic direct shear tests are considered, the model can be parametrised to fit either one of the two tests, but not both with the same set of parameters. Similar problems can be found in literature for monotonic and cyclic triaxial tests of railway ballast. In this work, the comparison between experiment and simulation is conducted using the entire data of the test, e.g. shear force over shear path curve from the direct shear test. In addition to a visual comparison of the results also quantitative errors based on the sum of squares are defined. To improve the fit of the DEM model to both types of experiments, an extension on the Hertz-Mindlin contact law is used, which introduces additional physical effects (e.g. breakage of edges or yielding). This model introduces two extra material parameters and is successfully parametrised. Using only one set of parameters, the results of the DEM simulation are in good accordance with both experimental cyclic compression test and monotonic directs shear test.

Friction phenomena and their impact on the shear behaviour of granular material.

In the discrete element simulation of granular materials, the modelling of contacts is crucial for the prediction of the macroscopic material behaviour. From the tribological point of view, friction at contacts needs to be modelled carefully, as it depends on several factors, e.g. contact normal load or temperature to name only two. In discrete element method (DEM) simulations the usage of Coulomb's law of friction is state of the art in modelling particle-particle contacts. Usually in Coulomb's law, for all contacts only one constant coefficient of friction is used, which needs to reflect all tribological effects. Thus, whenever one of the influence factors of friction varies over a wide range, it can be expected that the usage of only one constant coefficient of friction in Coulomb's law is an oversimplification of reality. For certain materials, e.g. steel, it is known that a dependency of the coefficient of friction on the contact normal load exists. A more tribological tangential contact law is implemented in DEM, where the interparticle friction coefficient depends on the averaged normal stress in the contact. Simulations of direct shear tests are conducted, using steel spheres of different size distributions. The strong influence of interparticle friction on the bulk friction is shown via a variation of the constant interparticle friction coefficient. Simulations with constant and stress-dependent interparticle friction are compared. For the stress-dependent interparticle friction, a normal stress dependency of the bulk friction is seen. In the literature, measurements of different granular materials and small normal loads also show a stress dependency of the bulk friction coefficient. With increasing applied normal stress, the bulk friction coefficient reduces both in the experiments and in the simulations.