• Structural Engineering and Mechanics
• /
• 제71권1호
• /
• pp.65-75
• /
• 2019

In this article, the dynamic fracture instability characteristics, including dynamic crack propagation and crack branching, in PMMA brittle solids under dynamic loading are investigated using the discrete element method (DEM) simulations. The microscopic parameters in DEM are first calibrated using the comparison with the previous experimental results not only in the field of qualitative analysis, but also in the field of quantitative analysis. The calibrating process illustrates that the selected microscopic parameters in DEM are suitable to effectively and accurately simulate dynamic fracture process in PMMA brittle solids subjected to dynamic loads. The typical dynamic fracture behaviors of solids under dynamic loading are then reproduced by DEM. Compared with the previous experimental and numerical results, the present numerical results are in good agreement with the existing ones not only in the field of qualitative analysis, but also in the field of quantitative analysis. Furthermore, effects of dynamic loading magnitude, offset distance of the initial crack and initial crack length on dynamic fracture behaviors are numerically discussed.

• 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
• /
• 제15권4호
• /
• pp.477-493
• /
• 2003

Masonry is not a simple material, the influence of mortar joints as a plane of weakness is a significant feature and this makes the numerical modelling of masonry very difficult especially when dynamic (seismic) analysis is involved. In order to develop a simple numerical model for masonry under earthquake load, an analytical model based on Distinct Element Method (DEM) is being developed. At the first stage, the model is applied to simulate the in-plane shear behaviour of an unreinforced masonry wall with and without opening where the testing results are available for comparison. In DEM, a solid is represented as an assembly of discrete blocks. Joints are modelled as interface between distinct bodies. It is a dynamic process and specially designed to model the behaviour of discontinuities. The numerical solutions obtained from the distinct element analysis are validated by comparing the results with those obtained from existing experiments and finite element modelling.

• 한국도로학회논문집
• /
• 제17권4호
• /
• pp.25-31
• /
• 2015

PURPOSES: Implementation and verification of the simple linear cohesive viscoelastic contact model that can be used to simulate dynamic behavior of sticky aggregates. METHODS: The differential equations were derived and the initial conditions were determined to simulate a free falling ball with a sticky surface from a ground. To describe this behavior, a combination of linear contact model and a cohesive contact model was used. The general solution for the differential equation was used to verify the implemented linear cohesive viscoelastic API model in the DEM. Sensitivity analysis was also performed using the derived analytical solutions for several combinations of damping coefficients and cohesive coefficients. RESULTS : The numerical solution obtained using the DEM showed good agreement with the analytical solution for two extreme conditions. It was observed that the linear cohesive model can be successfully implemented with a linear spring in the DEM API for dynamic analysis of the aggregates. CONCLUSIONS: It can be concluded that the derived closed form solutions are applicable for the analysis of the rebounding behavior of sticky particles, and for verification of the implemented API model in the DEM. The assumption of underdamped condition for the viscous behavior of the particles seems to be reasonable. Several factors have to be additionally identified in order to develop an enhanced contact model for an asphalt mixture.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2009년도 춘계 학술발표회
• /
• pp.136-145
• /
• 2009

Factors influencing shear behavior of granular material include particle size, shape, distribution, relative density, particle crushing, etc. In this study, these factors are characterized by viewpoint of shear behavior using numerical analysis based on DEM. Geometrical particle shape is represented by a combination of small circular particles and influence of particle shape on crushing is studied through relative comparisons between clump (uncrushable) and cluster (crushable) models which are modeled using DEM. Also, particle shape is quantified by the dimensionless parameters such as circularity and convexity. The results indicate that particle shape indexes have a negative association with internal friction angle. Also, internal friction angle becomes reduced and failure envelop curve becomes nonlinear due to the particle crushing. It is also found that numerical results are quite good agreement with the experimental test conducted in this study.

• 터널과지하공간
• /
• 제11권4호
• /
• pp.344-351
• /
• 2001

발파를 이용한 터널굴착에 있어서, 장약량 산정이나 터널 주변 암반 및 구조물웨 대한 영향 평가 등은 주로 간단한 경험식에 의해 이루어져 왔다. 또한 지금까지 발파의 충격에너지를 고려한 지하구조물의 동해석 연구는 매우 빈약한 실정이라 할 수 있다. 따라서 본 연구에서는 발파하중을 받는 불연속체 지하구조물의 비선형 거동을 평가하기 위하여 2차원 개별요소법의 수치모델기법을 개발하였고, 실제문제에의 적용성을 검토하기 위하여 제방하부에 위치한 터널발파의 수치해석을 통하여 시간별 지하 및 지상 구조물의 변위 및 입자속도의 전파과정을 살펴보았다. 해석결과 본 연구에서 제안한 발파하중의 모형이 실제 문제에 적용될 수 있었고, 이를 개별요소해석에 적용함으로써 발파지역 주변의 구조물의 안정성을 검토할 수 있었다.

• 한국항해항만학회지
• /
• 제31권9호
• /
• pp.785-791
• /
• 2007

In this paper we proposed a model that the deformation of the submerged rubble mound breakwaters composed with materials of various size, induced by wave action, can be computed. The water particle kinematics by waves in porous mound structure are computed by CADMAS-SURF, then the deformation of structure is computed using DEM module. To investigate the interaction of wave and sectional deformation of structures, analysis is accomplished by two steps. Analysis at the first step is executed with incipient mound. And analysis at the second step is executed with deformed mound by wave action. Furthermore, behaviors of materials are influenced by various properties such as the contact stiffness and the friction angle. Therefore, in order to present the behavior of the element caused by various properties, computations are accomplished with random coefficients by using the Monte Carlo simulation.

• 한국터널지하공간학회 논문집
• /
• 제19권6호
• /
• pp.1029-1044
• /
• 2017

개별요소법(Discrete Element Method, DEM)은 다수의 작은 입자들의 운동 및 상호영향을 계산하여 시스템의 거동을 해석하는 수치해석법으로써, 실제 화학공학, 약학, 토목공학, 재료과학, 식품공학 등 다양한 산업현장에서 적용되고 있다. 본 연구에서는 DEM 기법에 근거한 입자 역학 전용 해석 상용 소프트웨어를 사용하여 스포크타입 토압식 쉴드TBM 굴착성능을 평가하기 위한 예비 해석을 수행하였다. TBM에 대한 해석은 커터헤드의 회전속가 다른 2가지 조건에 대해 수행되었다. 해석을 진행하는 동안 커터헤드면에 작용하는 저항 토크, 커터헤드면과 쉴드면에 작용하는 압축력, 스크루 오거를 통해 배출되는 토사의 양을 검토하였다. 해석을 통해 DEM 해석을 이용한 TBM 장비 모델링의 적용성을 검토하였다.

• 터널과지하공간
• /
• 제30권5호
• /
• pp.484-495
• /
• 2020

EPB TBM의 굴진을 수치적으로 해석하기 위해 개별요소법(DEM, discrete element method), 유한요소법(FEM, finite element method), 유한차분법(FDM, finite difference method) 등과 같은 다양한 수치해석 기법이 적용되어 왔다. 본 논문에서는 이중 개별요소법과 유한차분법을 연계하는 방식을 채택하여 EPB TBM 굴진해석 모델링 방법을 제시하였다. 제시한 개별요소법-유한차분법 연계 TBM 굴진해석 모델에서 TBM이 굴착하는 굴착부는 개별요소법을 적용하였으며, 입자 접촉 물성치의 경우 일련의 삼축압축시험을 통해 교정하였다. 굴착부 주변지반은 유한차분법을 연계시켜 정지토압계수를 고려하여 굴착부에 수평지중응력을 구현할 수 있도록 하였다. 또한, 이를 통해 소요 입자 개수를 감소시켜 모델의 해석효율을 증대시켰다. 본 논문에서 제시한 수치해석 모델은 TBM의 굴진율, 커터헤드 및 스크류 컨베이어 회전속도 등을 조절할 수 있으며 TBM 굴진 중 토크, 추력, 챔버압, 배토량을 도출해 낼 수 있다.

• 한국지반신소재학회논문집
• /
• 제17권4호
• /
• pp.129-139
• /
• 2018

개별요소법은 입자형태의 재료를 다루는 분야에서 널리 사용되는 방법이다. 특히, 지반에서와 같이 변형이 크게 발생하는 특성을 지닌 재료의 경우에도 유용하게 사용할 수 있는 해석방법 중 하나이다. 본 연구에서는 개별요소법을 사용하여 지반을 생성하고 지반에 TBM 형상을 관입시킴으로써 굴진을 모사하고자 하였다. 해석에 사용된 TBM은 7.73m 직경의 토압식 쉴드 TBM을 사용하였으며 해석을 통해 굴진 성능을 예측할 수 있는 수치해석 모델을 검토하고자 하였다. 대상 모델의 스포크 개수는 8개이며, 커터헤드의 개구율은 약 21.31%이다. 또한, 커터헤드에는 52개 디스크 커터와 167개 커터비트가 굴착 도구로 장착되어 있으며, 커터헤드는 게이지 커터가 장착된 커터헤드 외주면이 굴곡진 세미돔 타입이다. 해석을 통해 커터헤드와 각각의 절삭 도구에 작용하는 반력과 저항 토크를 검토할 수 있었다. 또한 커터헤드 중심에서부터 거리에 따른 비교를 통해 커터헤드의 절삭 도구 위치별로 발생하는 반력과 토크를 검토하였다.