• 제목/요약/키워드: discrete element method

검색결과 499건 처리시간 0.031초

Experimental and numerical studies of the pre-existing cracks and pores interaction in concrete specimens under compression

  • Haeri, Hadi;Sarfarazi, Vahab;Zhu, Zheming;Marji, Mohammad Fatehi
    • Smart Structures and Systems
    • /
    • 제23권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.

Surface erosion of MICP-treated sands: Erosion function apparatus tests and CFD-DEM bonding model

  • Soo-Min Ham;Min-Kyung Jeon;Tae-Hyuk Kwon
    • Geomechanics and Engineering
    • /
    • 제33권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.

Mechanism of failure in the Semi-Circular Bend (SCB) specimen of gypsum-concrete with an edge notch

  • Fu, Jinwei;Sarfarazi, Vahab;Haeri, Hadi;Marji, Mohammad Fatehi;Guo, Mengdi
    • Structural Engineering and Mechanics
    • /
    • 제81권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)).

P1 비순응 요소를 이용한 비압축성 유동 문제의 위상최적화 (Topology Optimization of Incompressible Flow Using P1 Nonconforming Finite Elements)

  • 장강원;장세명
    • 대한기계학회논문집A
    • /
    • 제36권10호
    • /
    • pp.1139-1146
    • /
    • 2012
  • P1 비순응 요소를 이용하여 정상 비압축성 Navier-Stokes 유동의 위상최적화 문제를 푸는 방법을 제시한다. 본 연구는 Stokes 유동의 위상최적화 문제에 P1 비순응 요소를 적용하여 그 수치적 효용성을 보인바 있는 이전 연구에 대한 후속 연구이다. 비압축성 물질 해석에서 잠김현상이 발생하지 않으며 선형형상함수를 가지는 P1 비순응 요소의 장점이 관성항을 가지는 유체 문제의 해석과 설계에도 유효한 지를 파악하고자 한다. 일반적으로 사용되는 혼합정식화법과 비교하여 P1 비순응 요소의 사용은 벌칙 함수를 이용하여 연속 방정식을 따로 사용하지 않고 운동방정식에 부과할 수 있기 때문에 자유도의 개수를 감소시킬 수 있다. 벌칙 파라미터가 해의 정확도에 주는 영향과 적정 범위는 수치적으로 검토하도록 한다. 또한 보통의 사각 비순응 요소들이 요소면의 중앙에 절점을 가지고 고차의 형상함수를 지니는데 비하여, 본 연구에서 제시하는 P1 비순응 요소는 요소의 꼭지점에 절점을 가지고 {1, x, y}의 P1 형상함수로 구성됨으로써 수치적인 구현의 용이함이 일반 선형 사절점 요소와 동일하다. 제안한 방법의 효용성을 다양한 레이놀즈수에 따른 유동최적화 문제들을 살펴봄으로써 검증하도록 한다.

점성균열 모델을 위한 국부단위분할이 적용된 무요소법 (A Meshless Method Using the Local Partition of Unity for Modeling of Cohesive Cracks)

  • 지광습;정진규;김병민
    • 대한토목학회논문집
    • /
    • 제26권5A호
    • /
    • pp.861-872
    • /
    • 2006
  • 본 연구에서는 이차원 연속체에 존재하는 점성균열을 무요소법에서 국부 단위분할 원리에 근거하여 정식화하였다. 균열이 한 절점의 영향영역(domain of influence)을 완전히 통과하는 경우 그 절점의 형상함수는 계단함수로 확장되고, 균열 끝이 영향영역 내에 위치하는 경우 특이성이 제거된 가지함수(branch function)로 확장된다. 이러한 해의 영역의 확장은 국부 단위분할 원리를 만족하는 변위계에서만 이루어지므로, 약형 정식화는 표준 Galerkin방법에 의해서 얻어진다. 균열과 상호작용하는 영향영역만 확장되기 때문에, 성긴 형태의 시스템의 행렬을 유지하게 된다. 그러므로 확장에 의해 발생하는 계산비용의 증가는 최소화된다. 동적인 문제에서 균열성장에 관한 조건은 재료안정론으로부터 얻어졌다. 즉, 재료 한 점에서 어느 방향으로든 변형열화가 집중하게 되면, 그 방향에 점성균열을 삽입하여 연속체가 비연속체로 되도록 하였다. 균열의 성장속도도 같은 조건으로부터 자연스럽게 얻어졌다. 전통적인 무요소법보다 더 나은 정확도와 빠른 수렴성을 보이는 것이 확인되었으며, 이 기법의 적용성을 보이기 위해 잘 알려진, 정적 및 동적문제에 적용하였다.

암반지층 굴착벽체에 작용하는 토압에 대한 절리의 영향 (Effect of Joint on the Earth Pressure Against an Excavation Wall in Rockmass)

  • 손무락;솔로몬 아데도쿤
    • 대한토목학회논문집
    • /
    • 제34권2호
    • /
    • pp.505-513
    • /
    • 2014
  • 본 논문은 암반지층에서 다양한 암반종류 및 절리조건들을 고려하면서 굴착벽체에 작용하는 토압에 대한 절리의 영향을 조사한 것이다. 이를 위해 본 연구는 먼저 기존 선행연구에 대해서 살펴보고 그 다음으로 기존 선행연구의 한계를 극복하기 위하여 개별요소법에 근거한 수치해석적 매개변수해석을 수행하였다. 수치해석은 매개변수로서 암반종류 및 절리조건(절리면의 전단강도, 절리경사각 및 절리군의 수)을 고려하였으며, 지반과 굴착벽체의 상호작용하에 발생된 토압의 크기 및 분포특성이 조사되었다. 이와 더불어, 암반지층에서 발생된 토압과 토사지반에서의 경험토압인 Peck토압이 서로 비교되었다. 비교결과 절리가 형성된 암반지층에서 발생된 토압은 암반의 종류 및 절리조건에 따라서 크게 영향을 받았으며 토사지반에 대한 Peck의 경험토압과 서로 다른 결과를 나타내었다.

암반지층 굴착벽체 발생토압에 대한 단계별 굴착깊이의 영향 (Effect of Step-Wise Excavation Depth on the Earth Pressure against an Excavation Wall in Rock Mass)

  • 손무락;솔로몬 아데도쿤
    • 한국지반공학회논문집
    • /
    • 제30권2호
    • /
    • pp.43-52
    • /
    • 2014
  • 본 논문은 암반지층에 설치되는 굴착벽체에 발생하는 토압에 대한 단계별 굴착깊이의 영향을 조사한 것이다. 암반지층에서 발생하는 문제들을 다루기 위해 개별요소법에 근거한 수치해석적 매개변수해석을 수행하였다. 수치해석은 매개변수로서 단계별 굴착깊이, 암반종류 및 절리조건(절리면의 전단강도 및 절리경사각)을 고려하였다. 암반지층에서 발생된 토압과 토사지반에서의 경험토압인 Peck토압이 서로 비교되었다. 비교결과 암반지층에서 발생된 토압은 암반의 종류 및 절리조건에 따라서 크게 영향을 받았으며, 단계별 굴착깊이의 영향은 암반상태가 나빠질수록 증가하는 것으로 나타났다. 더불어, 암반지층에 설치된 굴착벽체에 발생하는 토압은 토사지반에 대한 Peck의 경험토압과는 크게 다를 수 있다는 결과를 얻었다.

고추수확기용 선별장치의 입자 거동 해석을 위한 고추 물성측정 (Measurement of Physical Properties of Pepper for Particle Behavior analysis of sorting system for Pepper Harvester)

  • 변준희;김수빈;김명호;김대철
    • 한국농업기계학회:학술대회논문집
    • /
    • 한국농업기계학회 2017년도 춘계공동학술대회
    • /
    • 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 입자거동해석 시뮬레이션을 통해 카드클리너 방식의 고추 선별기의 성능에 대한 분석을 하고자 한다.

  • PDF

점토지반에서 터널굴착에 따른 상부 블록구조물의 거동분석 (Response Analysis of Block-Bearing Structure due to Tunnel Excavation in Clay Ground)

  • 손무락
    • 대한토목학회논문집
    • /
    • 제34권1호
    • /
    • pp.175-183
    • /
    • 2014
  • 본 연구는 점토지반에서 터널의 굴착조건(터널깊이 및 터널직경)과 시공조건(지반손실량), 지반조건(연약점토, 단단한 점토)을 변화시키면서 터널상부에 위치한 구조물의 거동변화를 조사 및 분석한 것이다. 분석에 사용된 구조물은 4층 블록식구조물로서 변형 등에 의한 균열발생과 균열폭의 차이로 인해 구조물의 손상정도를 쉽게 파악할 수 있는 특징이 있다. 다양한 터널 굴착조건 및 시공조건, 지반조건에 대해서 발생할 수 있는 터널상부 블록구조물의 거동상태를 파악하기 위해 수치해석적 매개변수 해석을 수행하였으며, 수치해석은 구조물의 실제크랙 발생을 묘사할 수 있도록 개별요소법(DEM)에 근거하여 수행하였다. 다양한 매개변수 해석으로부터 얻어진 구조물의 거동상태에 대한 결과는 터널 굴착조건 및 시공조건, 점토 지반조건과 상호 연관하여 함께 반영될 수 있도록 도표화 하였으며, 이를 이용하여 향후 점토지반에서 다양한 터널굴착 및 시공조건, 지반조건으로 인해 유발되는 터널 상부구조물의 손상정도를 보다 용이하게 파악할 수 있을 것으로 기대된다.

교반볼밀을 이용한 밀링공정에서 각종실험조건에 따른 구리분말의 입자형상 변화 및 DEM 시뮬레이션에 의한 정량적 에너지 변화 (Particle Morphology Change and Quantitative Input Energy Variation during Stirred Ball Milling Process by DEM Simulation on Various Experimental Conditions)

  • 보르 암갈란;오양가;자갈사이항 바체첵;이재현;최희규
    • 한국재료학회지
    • /
    • 제28권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.