• Smart Structures and Systems
• /
• v.23 no.5
• /
• pp.479-493
• /
• 2019

In this paper, the interaction between notch and micro pore under uniaxial compression has been performed experimentally and numerically. Firstly calibration of PFC2D was performed using Brazilian tensile strength, uniaxial tensile strength and biaxial tensile strength. Secondly uniaxial compression test consisting internal notch and micro pore was performed experimentally and numerically. 9 models consisting notch and micro pore were built, experimentally and numerically. Dimension of these models are 10 cm*1 cm*5 cm. the length of joint is 2 cm. the angularities of joint are $30^{\circ}$, $45^{\circ}$ and $60^{\circ}$. For each joint angularity, micro pore was situated 2 cm above the lower tip of the joint, 2 cm above the middle of the joint and 2 cm above the upper of the joint, separately. Dimension of numerical models are 5.4 cm*10.8 cm. The size of the cracks was 2 cm and its orientation was $30^{\circ}$, $45^{\circ}$ and $60^{\circ}$. Diameter of pore was 1cm which situated at the upper of the notch i.e., 2 cm above the upper notch tip, 2 cm above the middle of the notch and 2 cm above the lower of the notch tip. The results show that failure pattern was affected by notch orientation and pore position while uniaxial compressive strength is affected by failure pattern.

• Geomechanics and Engineering
• /
• v.33 no.2
• /
• pp.133-140
• /
• 2023

Soil erosion can cause scouring and failures of underwater structures, therefore, various soil improvement techniques are used to increase the soil erosion resistance. The microbially induced calcium carbonate precipitation (MICP) method is proposed to increase the erosion resistance, however, there are only limited experimental and numerical studies on the use of MICP treatment for improvement of surface erosion resistance. Therefore, this study investigates the improvement in surface erosion resistance of sands by MICP through laboratory experiments and numerical modeling. The surface erosion behaviors of coarse sands with various calcium carbonate contents were first investigated via the erosion function apparatus (EFA). The test results showed that MICP treatment increased the overall erosion resistance, and the contribution of the precipitated calcium carbonate to the erosion resistance and critical shear stress was quantified in relation to the calcium carbonate contents. Further, these surface erosion processes occurring in the EFA test were simulated through the coupled computational fluid dynamics (CFD) and discrete element method (DEM) with the cohesion bonding model to reflect the mineral precipitation effect. The simulation results were compared with the experimental results, and the developed CFD-DEM model with the cohesion bonding model well predicted the critical shear stress of MICP-treated sand. This work demonstrates that the MICP treatment is effective in improving soil erosion resistance, and the coupled CFD-DEM with a bonding model is a useful and promising tool to analyze the soil erosion behavior for MICP-treated sand at a particle scale.

• Structural Engineering and Mechanics
• /
• v.81 no.1
• /
• pp.81-91
• /
• 2022

The effects of interaction between concrete-gypsum interface and edge crack on the failure behavior of the specimens in senicircular bend (SCB) test were studied in the laboratory and also simulated numerically using the discrete element method. Some quarter circular specimens of gypsum and concrete with 5 cm radii and hieghts were separately prepared. Then the semicircular testing specimens were made by attaching one gypsum and one concrete sample to one another using a special glue and one edge crack is produced (in the interface) by do not using the glue in that part of the interface. The tensile strengths of concrete and gypsum samples were separately measured as 2.2 MPa and 1.3 MPa, respectively. during all testing performances a constant loading rate of 0.005 mm/s were stablished. The proposed testing method showed that the mechanism of failure and fracture in the brittle materials were mostly governed by the dimensions and number of discontinuities. The fracture toughnesses of the SCB samples were related to the fracture patterns during the failure processes of these specimens. The tensile behaviour of edge notch was related to the number of induced tensile cracks which were increased by decreasing the joint length. The fracture toughness of samples was constant by increasing the joint length. The failure process and fracture pattern in the notched semi-circular bending specimens were similar for both methods used in this study (i.e., the laboratory tests and the simulation procedure using the particle flow code (PFC2D)).

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.36 no.10
• /
• pp.1139-1146
• /
• 2012

An alternative approach for topology optimization of steady incompressible Navier-Stokes flow problems is presented by using P1 nonconforming finite elements. This study is the extended research of the earlier application of P1 nonconforming elements to topology optimization of Stokes problems. The advantages of the P1 nonconforming elements for topology optimization of incompressible materials based on locking-free property and linear shape functions are investigated if they are also valid in fluid equations with the inertia term. Compared with a mixed finite element formulation, the number of degrees of freedom of P1 nonconforming elements is reduced by using the discrete divergence-free property; the continuity equation of incompressible flow can be imposed by using the penalty method into the momentum equation. The effect of penalty parameters on the solution accuracy and proper bounds will be investigated. While nodes of most quadrilateral nonconforming elements are located at the midpoints of element edges and higher order shape functions are used, the present P1 nonconforming elements have P1, {1, x, y}, shape functions and vertex-wisely defined degrees of freedom. So its implentation is as simple as in the standard bilinear conforming elements. The effectiveness of the proposed formulation is verified by showing examples with various Reynolds numbers.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5A
• /
• pp.861-872
• /
• 2006

The element free Galerkin method is extended by the local partition of unity method to model the cohesive cracks in two dimensional continuum. The shape function of a particle whose domain of influence is completely cut by a crack is enriched by the step enrichment function. If the domain of influence contains a crack tip inside, it is enriched by a branch enrichment function which does not have the LEFM stress singularity. The discrete equations are obtained directly from the standard Galerkin method since the enrichment is only for the displacement field, which satisfies the local partition of unity. Because only particles whose domains of influence are influenced by a crack are enriched, the system matrix is still sparse so that the increase of the computational cost is minimized. The condition for crack growth in dynamic problems is obtained from the material instability; when the acoustic tensor loses the positive definiteness, a cohesive crack is inserted to the point so as to change the continuum to a discontiuum. The crack speed is naturally obtained from the criterion. It is found that this method is more accurate and converges faster than the classical meshless methods which are based on the visibility concept. In this paper, several well-known static and dynamic problems were solved to verify the method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.2
• /
• pp.505-513
• /
• 2014

This paper investigated the effect of joint on the earth pressure against an excavation wall in rockmass with the consideration of various rock and joint conditions. For this purpose, this study briefly reviewed of the previous earth pressure studies, and then numerical parametric studies were conducted based on the Discrete Element Method (DEM) to overcome the limitations of the previous studies. The numerical tests were carried out with the controlled parameters including rock types and joint conditions (joint shear strength, joint inclination angle, and joint set), and the magnitude and distribution characteristics of the induced earth pressure were investigated considering the interactions between the ground and the excavation wall. In addition, the earth pressures induced in rock stratum were compared with Peck's earth pressure for soil ground. The results showed that the earth pressure against an excavation wall in jointed rockmass were highly affected by different rock and joint conditions and thus different from Peck's empirical earth pressure for soil ground.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.2
• /
• pp.43-52
• /
• 2014

This paper examined the effect of step-wise excavation depth on the earth pressure against an excavation wall in rock mass. Numerical parametric studies were conducted based on the Discrete Element Method (DEM) to carry out the problems in rock mass. Controlled parameters included step-wise excavation depth, rock types, and joint conditions (joint shear strength and joint inclination angle). The magnitude and distribution characteristics of the induced earth pressure in a jointed rock mass were investigated and compared with Peck's earth pressure for soil ground. The results showed that the earth pressure against an excavation wall in rock mass were highly affected by different rock and joint conditions, and the effect of step-wise excavation depth increased as a rock type is deteriorated more. In addition, it was found that the earth pressure against an excavation wall in rock mass might be considerably different from Peck's empirical earth pressure for soil ground.

• Proceedings of the Korean Society for Agricultural Machinery Conference
• /
• 2017.04a
• /
• pp.9-9
• /
• 2017

입자거동해석소프트웨어(EDEM)은 DEM(Discrete Element Method)기법을 이용한 입자 거동 전용 해석툴로 입자 유입량, 위치 등을 조절하여 입자거동과 관련된 제품 개발, 프로세스 최적화를 위한 비용 및 시간 절감에 활용도가 뛰어난 소프트웨어이다. EDEM을 활용하기 위해선 적용대상에 대한 물성치를 적용하여야 한다. 따라서 본 연구에서는 EDEM를 이용하여 현재 연구개발 중인 카드클리너 방식의 고추 선별기의 성능을 분석을 수행하기 위해 고추 물성측정 실험을 수행 하였다. EDEM을 이용한 입자거동해석에 필요한 개인 물성치에는 포아송비, 전단탄성계수, 밀도가 있다. 또한 입자-입자, 입자-Geometry 간의 상호관계를 위한 물성치인 반발계수, 정지마찰계수, 구름마찰계수가 필요하다. 공시 시료인 고추는 광주광역시 남구 승촌동 소재의 개인농가 Plastic 온실로 재배된 '천상'품종을 사용하였다. 푸아송 비와, 전단 탄성계수를 측정하기 위한 인장시험기기로는 만능인장시험기(TA-XT2, Stable Micro, 영국)를 이용하였으며, 인장에 의한 고추의 변형량 축정은 초고속카메라(NX4-SI, IDT, 미국)을 이용하였다. 밀도는 비중병법에 기초하여 질량과 부피를 측정하여 밀도를 계산하였다. 반발계수는 고추의 충돌 실험을 통해 변화한 높이를 이용하여 계산하였고, 충동 실험을 통해 변화한 높이는 초고속카메라를 이용하여 측정하였다. 정지마찰계수와 구름마찰계수는 고추의 미끄러짐이 시작하는 각도와 등속도 운동으로 구르는 각도를 초고속카메라를 이용하여 측정 후 계산하였다. 모든 실험은 3번 반복을 통해 평균값을 시험 결과 값으로 이용하였다. 고추의 대한 물성치 실험결과 고추의 푸아송 비는 0.294(std : 0.2), 전단탄성계수 4.624E+06 Pa, 밀도 $600kg/m^3$로 나타났다. 또한 입자-입자 간의 물성치인 반발계수는 0.383, 정지마찰계수는 0.455, 구름마찰 계수는 0.043로 나타났다. 추후 본 연구에서 측정한 고추의 물성치를 적용한 EDEM 입자거동해석 시뮬레이션을 통해 카드클리너 방식의 고추 선별기의 성능에 대한 분석을 하고자 한다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.1
• /
• pp.175-183
• /
• 2014

This study investigates the response of structures to tunnelling-induced ground movements in clay ground, varying tunnel excavation condition (tunnel depth and diameter), tunnel construction condition (ground loss), and tunnel ground condition (soft clay and stiff clay). Four-story block-bearing structures have been used because the structures can easily be characterized of the extent of damages with crack size and distribution. Numerical parametric studies have been used to investigate of the response of structures to varying tunnelling conditions. Numerical analysis has been conducted using Discrete Element Method (DEM) to have real cracks when the shear and tensile stress exceed the maximum shear and tensile strength. The results of structure responses from various parametric studies have been integrated to consider tunnel excavation condition, tunnel construction condition, and tunnel ground condition and provide a relationship chart among them. Using the chart, the response of structures to tunnelling can easily be evaluated in practice in clay ground.

• Korean Journal of Materials Research
• /
• v.28 no.3
• /
• pp.148-158
• /
• 2018

This study investigated the effect of the grinding media of a ball mill under various conditions on the raw material of copper powder during the milling process with a simulation of the discrete element method. Using the simulation of the three-dimensional motion of the grinding media in the stirred ball mill, we researched the grinding mechanism to calculate the force, kinetic energy, and medium velocity of the grinding media. The grinding behavior of the copper powder was investigated by scanning electron microscopy. We found that the particle size increased with an increasing rotation speed and milling time, and the particle morphology of the copper powder became more of a plate type. Nevertheless, the particle morphology slightly depended on the different grinding media of the ball mill. Moreover, the simulation results showed that rotation speed and ball size increased with the force and energy.