• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

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

• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.197-203
• /
• 1998

The discrete element method is a numerical model capable of describing the mechanical behaviour of assemblies of discs and spheres. The method is based on the use of an explicit numerical scheme in which the interaction of the particles is monitored contact by contact and the motion of the particles modelled particle by particle. In this paper, A two-dimensional model for computing contacts and motions of granular particles of unform, inelasticity is presented. And, code is developed. The primary aim of this paper is to approv computational result of continuum alaysis which is on processing. The end of this paper, that code is tested with several examples.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.412-415
• /
• 2010

Research about pedestrian flow in the view of sociology and psychology has been studied for more than a few decade. Due to the advance of computational facility, computational study for pedestrian flow extended to the field of architecture and traffic engineering. However, there is few study for the extremely high dense condition where pedestrian flow is driven by contact force among pedestrian. In this research, we analyze highly dense pedestrian flow using discrete element method

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2011.06a
• /
• pp.539-540
• /
• 2011

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 2006.09a
• /
• pp.187-188
• /
• 2006

In contrast with the Finite Element Method, the Discrete Element Method (DEM) takes explicitly into account the particulate nature of powders. DEM exhibits some drawbacks and many advantages. Simulations can be computationally expensive and they are only able to represent a volume element. However, these simulations have the great advantage of providing a wealth of information at the microstructural level. Here we demonstrate that the method is well suited for modelling, in coordination with FEM, the compaction of ceramic $UO_2$ particles that have been aggregated. Aggregates of individual ceramic crystallites that are strongly bonded together are represented by porous spheres.

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

• Journal of the Korean Geotechnical Society
• /
• v.23 no.8
• /
• pp.59-68
• /
• 2007

The comparative study using the lab-scale experiment and the discrete element analysis is attempted to analyze the cross-anisotropic elasticity of granular soils. The lab-scale experiment consists of the small stress-controlled triaxial cyclic tests and the bender element tests. In the discrete element analysis the simulations of lab-scale cyclic tests are conducted in the various directions. Good agreement between the experimental data and the simulation on the elastic properties in the axial and shear directions confirms the usefulness of the discrete element method. The comparative analysis of the difference in the experimental data and the simulation of radial cyclic tests shows that the discrete element method can successfully be used to check the reasonable magnitude of each measurement in the experiments.

• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2001.11a
• /
• pp.166-171
• /
• 2001

In order to improve the efficiency of the design process and the quality of the resulting design for the application-based exclusive rolling element bearings, this study propose design methodologies by using a genetic-based combinatorial optimization. By the presence of discrete variables such as the number of rolling element (standard component) and by the engineering point of views, the design problem of the rolling element bearing can be characterized by the combinatorial optimization problem as a fully discrete optimization. A genetic algorithm is used to efficiently find a set of the optimum discrete design values from the pre-defined variable sets. To effectively deal with the design constraints and the multi-objective problem, a ranking penalty method is suggested for constructing a fitness function in the genetic-based combinatorial optimization. To evaluate the proposed design method, a robust performance analyzer of ball bearing based on quasi-static analysis is developed and the computer program is applied to some design problems, 1) maximize fatigue life, 2) maximize stiffness, 3) maximize fatigue life and stiffness, of a angular contact ball bearing. Optimum design results are demonstrate the effectiveness of the design method suggested in this study. It believed that the proposed methodologies can be effectively applied to other multi-objective discrete optimization problems.

• Steel and Composite Structures
• /
• v.47 no.3
• /
• pp.365-374
• /
• 2023

This paper presents a new scheme for constructing locking-free finite elements in thick and thin plates, called Discrete Shear Gap element (DSG), using multiphase material topology optimization for triangular elements of Reissner-Mindlin plates. Besides, common methods are also presented in this article, such as quadrilateral element (Q4) and reduced integration method. Moreover, when the plate gets too thin, the transverse shear-locking problem arises. To avoid that phenomenon, the stabilized discrete shear gap technique is utilized in the DSG3 system stiffness matrix formulation. The accuracy and efficiency of DSG are demonstrated by the numerical examples, and many superior properties are presented, such as being a strong competitor to the common kind of Q4 elements in the static topology optimization and its computed results are confirmed against those derived from the three-node triangular element, and other existing solutions.