• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.22-25
• /
• 1999

A dynamic explicit finite element code for simulating sheet forming processes has been developed The code utilises the discrete Kirchhoff shell element and contact force is treated by a conventional penalty method. In order to reduce the computational cost a new and robust contact searching algorithm has been developed and implemented into the code. in the method a hierarchical structure of tool segments called a tree structure is built for each tool at the initial stage of the analysis Tree is built in a way to divide a trunk to 8 sub-trunk 2 in each direction until the lowest level of the tree(leaf) contains exactly one segment of the tool. In order to have a well-balanced tree each box on each sub level contains one eighth of the segments. Then at each time step contact line from a node comes out of the surface of the tool. Simulation of various sheet forming processes were performed to verify the validity of the developed code with main focus on he usefulness of the developed contact searching algorithm.

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

• Transactions of Materials Processing
• /
• v.9 no.3
• /
• pp.274-283
• /
• 2000

A dynamic explicit finite element code for simulating sheet forming processes has been developed. The code utilizes the discrete Kirchhoff shell element and contact force is treated by a conventional penalty method. In order to reduce the computational cost, a new and robust contact searching algorithm has been developed and implemented into the code. In the method, a hierarchical structure of tool segments is built for each tool at the initial stage of the analysis. hierarchical structure is built in a way to divide a box to 8 sub-boxes, 2 in each direction, until the lowest level of the hierarchical structure contains exactly one segment of the tool or empty. Then at each time step, contact is checked from the box to sub-boxes hierarchically for each node. Comparisons of computational results of various examples with the existing ones show the validity of the method.

• Computational Structural Engineering
• /
• v.5 no.2
• /
• pp.105-118
• /
• 1992

Static performance was compared for the triangular plate elements through some numerical experiments. Four Kirchhoff elements and six Mindlin elements were selected for the comparison. Numerical tests were executed for the problems of rectangular plates with regular and distorted meshes, rhombic plates, circular plates and cantilever plates. Among the Kirchhoff 9 DOF elements, the discrete Kirchhoff theory element was the best. Element distortion and the aspect ratio were shown to have negligible effects on the displacement behaviour. The Specht's element resulted in better results than the Bergan's but it was sensitive to the aspect ratio. The element based on the hybrid stress method also resulted in good results but it assumed to be less reliable. Among the linear Mindlin elements, the discrete shear triangle was the best in view of reliability, accuracy and convergence. Since the thin plate behaviour of it was as good as the DKT element, it can be used effectively in the finite element code regardless of the thickness. As a quadratic Mindlin element, the MITC7 element resulted in best results in almost all cases considered. The results were at least as good as those of doubly refined meshes of linear elements.

• Tunnel and Underground Space
• /
• v.7 no.1
• /
• pp.20-30
• /
• 1997

A discrete joint finite element with joint surface degradation was developed to investigate the shear behavior of rough rock joint. Isoparametric formulation was used for facilitating the implementation of the element in existing Finite Element Codes. The elasto-plastic joint deformation model with the discontinuity constitutive law proposed by Plesha was applied to the element. The reliability of the developed finite element code was successfully testified through numerical direct shear tests conducted under both constant normal stress and constant normal displacement conditions. The result of the numerical direct shear test showed that the code can capture characteristic deformation features envisaged in the direct shear test of rough rock joint.

• Computers and Concrete
• /
• v.2 no.5
• /
• pp.345-366
• /
• 2005

This work presents an assessment of the computational performance of a vector-parallel implementation of probabilistic model for concrete cracking in 3D. This paper shows the continuing efforts towards code optimization as reported in earlier works Paz, et al. (2002a,b and 2003). The probabilistic crack approach is based on the direct Monte Carlo method. Cracking is accounted by means of 3D interface elements. This approach considers that all nonlinearities are restricted to interface elements modeling cracks. The heterogeneity governs the overall cracking behavior and related size effects on concrete fracture. Computational kernels in the implementation are the inexact Newton iterative driver to solve the non-linear problem and a preconditioned conjugate gradient (PCG) driver to solve linearized equations, using an element by element (EBE) strategy to compute matrix-vector products. In particular the paper analyzes code behavior using OpenMP directives in parallel vector processors (PVP), such as the CRAY SV1 and CRAY T94. The impact of the memory architecture on code performance, and also some strategies devised to circumvent this issue are addressed by numerical experiment.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.329-342
• /
• 2008

A local behavior law, which includes elasticity, plasticity and damage, is developed in a three dimensional numerical model for concrete. The model is based on the Discrete Element Method (DEM)and the computational implementation has been carried out in the numerical Code YADE. This model was used to study the response of a concrete slab impacted by a rigid missile, and focuses on the extension of the compacted zone. To do so, the model was first used to simulate compression and hydrostatic tests. Once the local constitutive law parameters of the discrete element model were calibrated, the numerical model simulated the impact of a rigid missile used as a reference case to be compared to an experimental data set. From this reference case, simulations were carried out to show the importance of compaction during an impact and how it expands depending on the different impact conditions. Moreover, the numerical results were compared to empirical predictive formulae for penetration and perforation cases, demonstrating the importance of taking into account the local compaction process in the local interaction law between discrete elements.

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

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.123.1-123.1
• /
• 2010

수평형 지중열교환기의 최적설계를 위해서는 되메움재의 광물특성 및 입자크기, 열전도도(thermal conductivity), 열용량(heat capacity)등과 같은 열적 특성을 파악 하는 것은 중요하다. 수평형 지중 열교환기용 되메움재의 열전도도를 파악하기 위해 비정상 열선법을 적용한 QTM-500을 사용하여 포화도에 따른 천연규사-물-공기 혼합물의 열전도도를 측정하였다. 측정된 열전도도를 개별요소법(Discrete Element Mothod)에 근거한 2차원 수치해석 프로그램인 PFC2D(Particle Flow Code in 2 Dimension)를 이용하여 비교 분석하였다. 수치해석에서는 혼합물의 건조밀도를 일정하게 유지한 상태에서 포화도에 따라 가상의 물 입자 개수를 변화시켰다. 개별요소법을 이용한 열전달 수치해석에서는 입자의 접촉을 통해 발생한 thermal pipe에 의해 열전달이 이루어진다. 이러한 thermal pipe의 열전도도는 접촉된 두 입자의 열전도도와 접촉면의 평균 열전도도를 고려하여 적용하였다.