• Interaction and multiscale mechanics
• /
• 제1권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
• /
• 제55권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년도 추계 학술대회논문집
• /
• 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년 춘계학술대회논문집
• /
• 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

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 2011년도 춘계학술대회 논문집
• /
• pp.539-540
• /
• 2011

• 한국분말야금학회:학술대회논문집
• /
• 한국분말야금학회 2006년도 Extended Abstracts of 2006 POWDER METALLURGY World Congress Part 1
• /
• 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.

• 한국철도학회논문집
• /
• 제6권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.

• 한국지반공학회논문집
• /
• 제23권8호
• /
• pp.59-68
• /
• 2007

사질토의 직교 이방 탄성계수를 실험적 방법과 개별 요소법을 이용하여 측정하고 결과를 정량적으로 비교-분석하였다. 실내 시험에서는 미소 응력 제어 반복재하시험과 벤더 엘리먼트의 측정 결과를 조합하여 탄성 특성을 구하였으며, 개별 요소법에서는 실내 시험 방법과 상응하는 반복재하시험을 모사하였다. 축 방향 및 전단방향에 관련된 탄성 특성의 비교 결과는 제한된 조건에서도 개별 요소법이 복잡한 실내 시험 과정을 모사할 수 있음을 보여주었다. 방사 방향 재하시험 결과에서 나타난 두 방법 간 차이의 비교 분석 결과는 개별 요소법이 실내 시험으로 구한 개별 측정치의 정량적 타당성을 확인하는데 충분히 활용될 수 있음을 보여준다.

• 한국윤활학회:학술대회논문집
• /
• 한국윤활학회 2001년도 제34회 추계학술대회 개최
• /
• 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
• /
• 제47권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.