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

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 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

Some important infrastructures like roads, railway tracks or dams were constructed in places threatened by natural hazards. With the purpose of preserving these infrastructures from landslides and rock-falls, different containment systems are installed, and one of the most popular are the flexible metallic fences. The development of full-scale laboratory tests to evaluate the behaviour of flexible metallic fences is unfunctional, accounting to the huge magnitude of the event. On the other hand, small-scale testing may lead to inaccurate results, due to the distortion in the contours (e.g. anchors of the metallic fences). These problems in laboratory testing have led to the popularization of the use of numerical methods. In this study, the bonded Discrete Element Method (DEM) is used for the analysis of the behaviour of flexible metallic fences for rockfall protection. The bonded DEM is a modification of the classical DEM which assumes that bonds exist between particles, resisting their separation. In this case, the net cables are represented using rigid spheres joined by bond elements that are deformed according to an elasto-plastic law. Calculations were carried out using the DEMPack program, a specific software developed in CIMNE for modelling with the bonded DEM. This software allows considering the inter-action between discrete and finite elements, which can be useful to represent the boundaries of the domain, such as the surface of the slope. The code is firstly validated reproducing benchmark tests available in the literature. Finally, full-scale tests are computed in order to evaluate the energy dissipation capacity of the fence during a rockfall event.

Abstract

We present some developments in the Particle Finite Element Method (PFEM) for analysis of complex coupled problems in mechanics involving fluid-soil-structure interaction (FSSI). The PFEM uses an updated Lagrangian description to model the motion of nodes (particles) in both the fluid and the solid domains (the later including soil/rock and structures). A mesh connects the particles (nodes) defining the discretized domain where the governing equations for each of the constituent materials are solved as in the standard FEM. The stabilization for dealing with an incompressibility continuum is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the non linear transient coupled FSSI problem is described. The procedure to model frictional contact conditions and material erosion at fluid-solid and solid-solid interfaces is described. We present several examples of application of the PFEM to solve FSSI problems such as the motion of rocks by water streams, the erosion of a river bed adjacent to a bridge foundation, the stability of breakwaters and constructions sea waves and the study of landslides.

Abstract

Wedge shaped blocks spillways are an innovative solution that allows spilling over the downstream shoulder of earth and rock-fill dams in a safe way. However, they have been barely used as main spillway due to the lack of practical design criteria. An innovative procedure for computer-aided design of wedge shaped blocks spillways is presented in this paper. It includes the design of the drainage and supporting layer considering its seepage capacity and stability. The leakage flow through the joints between blocks is estimated by means of a numerical model calibrated and validated from experimental results. Stability against sliding of the downstream shoulder or the drainage layer is ensured, considering the properties of the granular material selected by the designer and non-linear resistance laws. The procedure will also suggest the shape of a toe protection.

Abstract

The static and seismic analysis of Janneh arch-gravity dam (157 m) is carried out by considering a combination of self-weight, hydrostatic, uplift and seismic loads. Linear and nonlinear analyses are performed for both static and seismic cases. Nonlinear behavior is studied by means of joint elements in the contact between the rock foundation and the dam following a bilinear cohesive law. The hydrodynamic effects derived to seismic loads are also considered. Hydrodynamic Westergaard’s approach is applied in the pseudo-static analysis, and the Westergaard’s generalized added mass is used for time-history case.

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

Landslides in reservoirs constitute one of the main threats to dams. Although they are not frequent, they may generate severe damages in the population and goods placed near the shores. It is a complex phenomenon, because of the interaction between the landslide, the still water in the reservoir, and the dam. The Particle Finite Element Method (PFEM) is an innovative numerical scheme which combines a lagrangian approach with the resolution of FEM equations via mesh generation. In previous works, it had been applied to the simulation of landslides in reservoirs, considering the moving body either as a rigid solid or as a Newtonian fluid. Despite the promising results of the simulation of some reference cases, these constitutive models showed poor performance when applied to other relevant events, such as Vaiont case. Current work introduces the implementation of a non Newtonian Bingham model in the particle finite element method scheme, and the results of its application to several benchmark cases. Some of them have been used to calibrate the parameters of the Bingham model, whereas others served as test cases. The phenomenon presents several intrinsic uncertainties, such as the lack of precise information about the physical properties of the landslide, its geometry, its granulometry, the drag coefficient along the sliding surface, among others. In spite of that, results show that the Bingham model implemented can be useful for estimating the potential consequences of landslides in reservoirs. This implementation makes the PFEM a versatile numerical tool for the analysis of landslides in reservoirs, provided that from now on both rigid solids, Newtonian and non Newtonian fluids can be accounted for. The authors are currently working in the optimization of the numerical code, so that more complex problems could be afforded in more detail. In addition, the erosion module is being validated, in order to reproduce effects such as the dragging of the underlying material due to the landslide, as well as the effect of overtopping on dam stability, particularly for rockfill dams.

Abstract

Double curvature arch dams feature geometrical complexity with a significant amount of parameters involved. Different criteria exist to assist in the design task, from simplified geometrical approaches to optimization procedures. However, most of them present a lack of flexibility and are not integrated in computer-aided design tools. In this contribution, an interactive and flexible software tool is presented to support the complete design process: geometrical definition, FEM model generation (including the mesh, the loads and the boundary conditions) and thermo-mechanical analysis. The design can be performed with different levels of detail to adapt to the information available in each stage of the project. The tool allows defining the shape of the reference cylinder, the excavation depth and slope along the foundation, the crown cantilever thickness and curvature, the shape and location of horizontal arcs; all these steps were described in former contributions. Here, special attention is paid to the introduction of additional features such as joints, spillways, abutments of varying shape and outlet works. All steps have been defined with a high degree of flexibility in the design process. The tool is integrated with the pre and post process software GiD, which allows taking advantage of its functionalities, such as mesh generation and results analysis. It is also coupled with a specific application for thermo-mechanical analysis of dams, developed in Kratos Multiphysics – a framework for building parallel multi-disciplinary simulation software. The whole design process can be followed in a unique environment, because the structural response of preliminary designs can be computed and the results considered to refine the dam geometry.