COLLECTIONS

Documents published in Scipedia

• A review on thermo-mechanical modelling of arch dams during construction and operation. Effect of the reference temperature on the stress field

  F. Salazar, D. Vicente, J. Irazábal González, I. de-Pouplana, J. San Mauro

  Archives of Computational Methods in Engineering (2020). 27 (5), 1681-1707

  
  

  Abstract

  Double-curvaturedamsareunique structures for several reasons. Their behaviour changes significantly after joint grouting, when they turn from a set of independent cantilevers into a monolithic structure with arch effect. The construction process has a relevant influence on the stress state, due to the way in which self-weight loads are transmitted, and to the effect on the dissipation of the hydration heat. Temperature variations in the dam body with respect to those existing at joint grouting generate thermal stresses that may be important in the stress state of the structure. It is thus essential to have a realistic estimate of this thermal field, also called reference or

  Abstract
  Double-curvaturedamsareunique structures for several reasons. Their behaviour changes significantly after joint grouting, when they turn from a set of independent cantilevers [...]

  • 10
  •  read

• Evaluation of the performance of ballasted railway structures subjected to high-speed train loads by the Discrete Element Method

  . s:, J. Irazábal, F. Salazar, E. Oñate, J. Irazábal González

  J. Irazábal's Conference proceedings (2019). 1

  
  

  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
  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 [...]

  • 0
  •  read

• Full-scale numerical calculation of ballasted tracks with the Discrete Element Method

  J. Irazábal, F. Salazar, E. Oñate, J. Irazábal González

  J. Irazábal's Conference proceedings (2019). 1

  
  

  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
  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 [...]

  • 16
  •  read

• Effect of the integration scheme on the rotation of non-spherical particles with the discrete element method

  J. Irazábal González, F. Salazar, M. Santasusana, E. Oñate

  Computational Particle Mechanics (2019). Vol. 6 (4), pp 545–559

  
  

  Abstract

  The discrete element method (DEM) is an emerging tool for the calculation of the behaviour of bulk materials. One of the key features of this method is the explicit integration of the motion equations. Explicit methods are rapid, at the cost of a limited time step to achieve numerical stability. First- or second-order integration schemes based on a Taylor series are frequently used in this framework and shown to be accurate for the translational and rotational motion of spherical particles. However, they may lead to relevant inaccuracies when non-spherical particles are used since the orientation implies a modification in the second-order inertia tensor in the inertial reference frame. Specific integration schemes for non-spherical particles have been proposed in the literature, such as the fourth-order Runge–Kutta scheme presented by Munjiza et al. and the predictor–corrector scheme developed by Zhao and van Wachem which applies the direct multiplication algorithm for integrating the orientation. In this work, both methods are adapted to be used together with a velocity Verlet scheme for the translational integration. The performance of the resulting schemes, as well as that of the direct integration method, is assessed, both in benchmark tests with analytical solution and in real-scale problems. The results suggest that the fourth-order Runge–Kutta and the Zhao and van Wachem schemes are clearly more accurate than the direct integration method without increasing the computational time.

  Abstract
  The discrete element method (DEM) is an emerging tool for the calculation of the behaviour of bulk materials. One of the key features of this method is the explicit integration [...]

  • 13
  •  read

• Numerical modelling of landslide‐generated waves with the particle finite element method (PFEM) and a non‐Newtonian flow model

  F. Salazar, J. Irazábal González, A. LARESE, E. Oñate

  Int. J. Numer. Anal. Meth. Geomech. (2016). Vol. 40 (6), pp. 809-826

  
  

  Abstract

  Landslide‐generated impulse waves may have catastrophic consequences. The physical phenomenon is difficult to model because of the uncertainties in the kinematics of the mobilised material and to the intrinsic complexity of the fluid–soil interaction. The particle finite element method (PFEM) is a numerical scheme that has successfully been applied to fluid–structure interaction problems. It uses a Lagrangian description to model the motion of nodes (particles) in both the fluid and the solid domains (the latter including soil/rock and structures). A mesh connecting the particles (nodes) is re‐generated at every time step, where the governing equations are solved. Various constitutive laws are used for the sliding mass, including rigid solid and the Newtonian and non‐Newtonian fluids. Several examples of application are presented, corresponding both to experimental tests and to actual full‐scale case studies. The results show that the PFEM can be a useful tool for analysing the risks associated with landslide phenomena, providing a good estimate to the potential hazards even for full‐scale events.

  Abstract
  Landslide‐generated impulse waves may have catastrophic consequences. The physical phenomenon is difficult to model because of the uncertainties in the kinematics of the [...]

  • 9
  •  read

• Numerical modelling of granular materials with spherical discrete particles and the bounded rolling friction model. Application to railway ballast

  J. Irazábal González, F. Salazar, E. Oñate

  Computers and Geotechnics (2017). Vol. 85, pp. 220-229

  
  

  Abstract

  The Discrete Element Method (DEM) was found to be an effective numerical method for the calculation of engineering problems involving granular materials. However, the representation of irregular particles using the DEM is a very challenging issue, leading to different geometrical approaches. This document presents a new insight in the application of one of those simplifications known as rolling friction, which avoids excessive rotation when irregular shaped materials are simulated as spheric particles. This new approach, called the Bounded Rolling Friction model, was applied to reproduce a ballast resistance test.

  Abstract
  The Discrete Element Method (DEM) was found to be an effective numerical method for the calculation of engineering problems involving granular materials. [...]

  • 7
  •  read

• Advances in the modelling of railway ballast using the Discrete Element Method (DEM)

  J. Irazábal González, F. Salazar, E. Oñate

  (2018). Fourth International Conference on Railway Technology: Research, Development and Maintenance (RAILWAYS 2018))

  
  

  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
  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 [...]

  • 10
  •  read

• Shape characterization of railway ballast stones for discrete element calculations

  J. Irazábal González, F. Salazar, E. Oñate

  (2017). V International Conference on Particle-Based Methods (PARTICLES 2017)

  
  

  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
  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 [...]

  • 10
  •  read

• Geometric representation of railway ballast using the Discrete Element Method (DEM)

  J. Irazábal González, F. Salazar, E. Oñate

  (2016). European Congress on Computational Methods in Applied Sciences and Engineering

  
  

  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
  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 [...]

  • 10
  •  read

• Numerical modelling with discrete elements of rockfall protection systems

  J. Irazábal González, S. Latorre, F. Salazar, M. Celigueta, E. Oñate

  (2017). V International Conference on Particle-Based Methods (PARTICLES 2017)

  
  

  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
  Some important infrastructures like roads, railway tracks or dams were constructed in places threatened by natural hazards. With the purpose of preserving these infrastructures [...]

  • 11
  •  read

Articles (Documents not available in Scipedia)

• A Review on Thermo-mechanical Modelling of Arch Dams During Construction and Operation: Effect of the Reference Temperature on the Stress Field

  Fernando Salazar, David J. Vicente, Ignasi de-Pouplana and Javier San Mauro

  Salazar, F., Vicente, D.J., Irazábal, J. et al. A Review on Thermo-mechanical Modelling of Arch Dams During Construction and Operation: Effect of the Reference Temperature on the Stress Field. Arch Computat Methods Eng 27, 1681–1707 (2020).

• Effect of the integration scheme on the rotation of non-spherical particles with the discrete element method

  Fernando Salazar, Miquel Santasusana, Eugenio Oñate

  Irazábal, J., Salazar, F., Santasusana, M. et al. Effect of the integration scheme on the rotation of non-spherical particles with the discrete element method. Comp. Part. Mech. 6, 545–559 (2019).

• Numerical modelling of granular materials with spherical discrete particles and the bounded rolling friction model. Application to railway ballast

  Fernando Salazar, Eugenio Oñate

  Irazábal, Joaquín, Fernando Salazar, and Eugenio Oñate. Numerical modelling of granular materials with spherical discrete particles and the bounded rolling friction model. Application to railway ballast. Computers and Geotechnics 85 (2017): 220-229.

• The Double Hierarchy Method. A parallel 3D contact method for the interaction of spherical particles with rigid FE boundaries using the DEM

  M Santasusana, J Irazábal, E Oñate, JM Carbonell

  Santasusana, M., Irazábal, J., Oñate, E. et al. The Double Hierarchy Method. A parallel 3D contact method for the interaction of spherical particles with rigid FE boundaries using the DEM. Comp. Part. Mech. 3, 407–428 (2016).

• Análisis experimental y numérico del oleaje producido por la inestabilidad de laderas de embalses

  R Alves, R Morán, MÁ Toledo Municio, J Irazábal González, ...

  Actas de las IV Jornadas de Ingeniería del Agua: la precipitación y los ..., 2015

• Numerical modelling of landslide‐generated waves with the particle finite element method (PFEM) and a non‐Newtonian flow model

  F Salazar, J Irazábal, A Larese, E Oñate

  Salazar, Fernando, et al. "Numerical modelling of landslide‐generated waves with the particle finite element method (PFEM) and a non‐Newtonian flow model." International Journal for Numerical and Analytical Methods in Geomechanics 40.6 (2016): 809-826.

• El papel de los modelos numéricos en la investigación y el diseño de aliviaderos de presas

  F Salazar, J San Mauro, J Irazábal, A Larese, R Rossi, E Oñate Ibáñez, ...

  Red de Laboratorios de Hidráulica de España (RLHE), 2013

• Validation of the Particle Finite Element Method (PFEM) for Simulation of Rock Slides in Lakes and Reservoirs

  J Irazábal, F Salazar, E Oñate

  CIMNE Technical Report, 2012