• Structural Engineering and Mechanics
• /
• v.55 no.2
• /
• pp.315-334
• /
• 2015

This study's primary aim is to check the existence of a representative volume element for granular materials and determine the link between the properties (responses) of macro structures and the size of the discrete particle assembly used to represent a constitutive relation in a two-scale model. In our two-scale method the boundary value problem on the macro level was solved using finite element method, based on the Cosserat continuum; the macro stresses and modulus were obtained using a solution of discrete particle assemblies at certain element integration points. Meanwhile, discrete particle assemblies were solved using discrete element method under boundary conditions provided by the macro deformation. Our investigations focused largely on the size effects of the discrete particle assembly and the radius of the particle on macro properties, such as deformation stiffness, bearing capacity and the residual strength of the granular structure. According to the numerical results, we suggest fitting formulas linking the values of different macro properties (responses) and size of discrete particle assemblies. In addition, this study also concerns the configuration and displacement fluctuation of discrete particle assemblies on the micro level, accompanied with the evolution of bearing capacity and deformation on the macro level.

• Computers and Concrete
• /
• v.6 no.2
• /
• pp.133-153
• /
• 2009

Discrete element modeling (DEM) in concrete technology is concerned with design and use of models that constitute a schematization of reality with operational potentials. This paper discusses the material science principles governing the design of DEM systems and evaluates the consequences for their operational potentials. It surveys the two families in physical discrete element modeling in concrete technology, only touching upon probabilistic DEM concepts as alternatives. Many common DEM systems are based on random sequential addition (RSA) procedures; their operational potentials are limited to low configuration-sensitivity features of material structure, underlying material performance characteristics of low structure-sensitivity. The second family of DEM systems employs concurrent algorithms, involving particle interaction mechanisms. Static and dynamic solutions are realized to solve particle overlap. This second family offers a far more realistic schematization of reality as to particle configuration. The operational potentials of this family involve valid approaches to structure-sensitive mechanical or durability properties. Illustrative 2D examples of fresh cement particle packing and pore formation during maturation are elaborated to demonstrate this. Mainstream fields of present day and expected application of DEM are sketched. Violation of the scientific knowledge of to day underlying these operational potentials will give rise to unreliable solutions.

• Journal of the Korea Concrete Institute
• /
• v.11 no.3
• /
• pp.181-192
• /
• 1999

The progressive failure following strain localization in concrete can be analyzed effectively using finite element modeling of fracture process zone of concrete with a finite element embedded discontinuity. In this study, a finite element with embedded discontinuous line is utilized for the analysis of progressive failure in concrete. The finite element with embedded discontinuity is a kind of discrete crack element, but the difficulties in discrete crack approach such as remeshing or adding new nodes along with crack growth can be avoided. Using a discontinuous shape function for this element, the displacement discontinuity is embedded within an element and its constitutive equation is modeled from the modeling of fracture process zone. The element stiffness matrix is derived and its dual mapping technique for numerical integration is employed. Then, a finite element analysis program with employed algorithms is developed and failure analysis results using developed finite element program are verified through the comparison with experimental data and other analysis results.

• Journal of the Korean Society for Railway
• /
• v.6 no.3
• /
• pp.174-178
• /
• 2003

This paper presents the development of 2-dimensional discrete element algorithm to generate circle and line elements for the simulation of the ballast and sleeper in railway. An example of randomly distributed circle elements show a good applicability of this algorithm for the modeling of the behaviors of ballast. The output about unbalaned force, particle velocity, and total energy conservation from the code is evaluated to check if the calculation is conducted properly.

• Geomechanics and Engineering
• /
• v.33 no.4
• /
• pp.353-367
• /
• 2023

The construction of a protective embankment is a suitable strategy to stop and control high-energy rock blocks' impacts during the rockfall phenomenon. In this paper, based on the discrete element numerical method, by modeling an existing embankment reinforced with geogrid, its stability status under the impact of a rock block with two types of low and high kinetic energy, namely 2402 and 4180 kJ, respectively, has been investigated. The modeling results show that the use of geogrid has caused the displacement in the front and back of the embankment to decrease by more than 30%. In this case, the reinforced embankment has stopped the rock block earlier. The displacements obtained from the DEM modeling are compared with the displacements measured from an actual practical experiment to evaluate the results' validity. Comparison between the results shows that the displacement values are close together, while the maximum percentage error in previous studies by an analytical method and the finite element method was 76.4% and 36.6%, respectively. Therefore, the obtained results indicate the discrete numerical method's high ability compared to other numerical and analytical methods to simulate and design the geogrid-reinforced soil embankment under natural disasters such as rockfall with a minor error.

• Earthquakes and Structures
• /
• v.11 no.4
• /
• pp.563-582
• /
• 2016

The failure mechanism and maximum collapse load of masonry structures may change significantly under static and dynamic excitations depending on their internal arrangement and material properties. Hence, it is important to understand correctly the nonlinear behavior of masonry structures in order to adequately assess their safety and propose efficient strengthening measures, especially for historical constructions. The discrete element method (DEM) can play an important role in these studies. This paper discusses possible collapse mechanisms and provides a set of parametric analyses by considering the influence of material properties and cross section morphologies on the out of plane strength of masonry walls. Detailed modeling of masonry structures may affect their mechanical strength and displacement capacity. In particular, the structural behavior of stacked and rubble masonry walls, portal frames, simple combinations of masonry piers and arches, and a real structure is discussed using DEM. It is further demonstrated that this structural analysis tool allows obtaining excellent results in the description of the nonlinear behavior of masonry structures.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1839-1844
• /
• 2003

In order to provide a tool for designing more efficient methods of mixing fabric, a simplified discrete element computational model was developed for modeling fabric dynamics in a rotating horizontal drum. Because modeling the interactions between actual pieces of fabric is quite complex, a simplified model was developed where individual pieces of bundled fabric are represented by spherical particles. The simulations are used to investigate fabric bundle kinematics, the power required to drive the rotating drum, and the power dissipated through normal and tangential contacts. Parametric studies indicate only fill percentage, drum rotation speed, and friction coefficient play significant roles in the fabric bundle dynamics.

• Computers and Concrete
• /
• v.27 no.4
• /
• pp.297-304
• /
• 2021

The hybrid Finite-Discrete Element (FDEM) approach combines aspects of both finite elements and discrete elements with fracture mechanics principles, and therefore it is well suited for realistic simulation of quasi-brittle materials. Notwithstanding, in the literature its application for the analysis of concrete is rather limited. In this paper, the proprietary FDEM code ELFEN is used to model concrete specimens under uniaxial compression and indirect tension (Brazilian tests) of different sizes. The results show that phenomena such as size effect and influence of strain-rate are captured using this modeling technique. In addition, a preliminary model of a slab subjected to dynamic shear punching due to progressive collapse is presented. The resulting fracturing pattern of the impacted slab is similar to observations from actual collapse.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.343-358
• /
• 2008

This article focuses on concrete structures submitted to impact loading and is aimed at predicting local damage in the vicinity of an impact zone as well as the global response of the structure. The Discrete Element Method (DEM) seems particularly well suited in this context for modeling fractures. An identification process of DEM material parameters from macroscopic data (Young's modulus, compressive and tensile strength, fracture energy, etc.) will first be presented for the purpose of enhancing reproducibility and reliability of the simulation results with DE samples of various sizes. The modeling of a large structure by means of DEM may lead to prohibitive computation times. A refined discretization becomes required in the vicinity of the impact, while the structure may be modeled using a coarse FE mesh further from the impact area, where the material behaves elastically. A coupled discrete-finite element approach is thus proposed: the impact zone is modeled by means of DE and elastic FE are used on the rest of the structure. The proposed approach is then applied to a rock impact on a concrete slab in order to validate the coupled method and compare computation times.

• Structural Engineering and Mechanics
• /
• v.85 no.3
• /
• pp.369-379
• /
• 2023

Reliability-Based Design Optimization (RBDO) is an appropriate framework for obtaining optimal designs by taking uncertainties into account. Large-scale problems with implicit limit state functions and problems with discrete design variables are two significant challenges to traditional RBDO methods. To overcome these challenges, this paper proposes a hybrid method to perform RBDO of structures that links Firefly Algorithm (FA) as an optimization tool to advanced (finite element) reliability methods. Furthermore, the Genetic Algorithm (GA) and the FA are compared based on the design cost (objective function) they achieve. In the proposed method, Weighted Simulation Method (WSM) is utilized to assess reliability constraints in the RBDO problems with explicit limit state functions. WSM is selected to reduce computational costs. To performing RBDO of structures with finite element modeling and implicit limit state functions, a First-Order Reliability Method (FORM) based on the Direct Differentiation Method (DDM) is utilized. Four numerical examples are considered to assess the effectiveness of the proposed method. The findings illustrate that the proposed RBDO method is applicable and efficient for RBDO problems with discrete and continuous design variables and finite element modeling.