• 터널과지하공간
• /
• 제16권6호
• /
• pp.476-485
• /
• 2006

발파에 의한 암반의 손상이나 파쇄는 폭약의 폭굉 과정에서 발생하는 충격파와 가스팽창의 영향에 의해 야기된다. 발파에 의한 파괴 메커니즘을 완전히 이해하기 위해서는 두 메커니즘을 같이 연구해야한다. 본 연구에서는 개별 요소법에 기초한 수치해석 프로그램인 PFC2D를 이용하여 발파공 벽면에 작용하는 폭굉압과 가스압을 동시에 모델링 할 수 있고 이에 따른 암반 내 균열 발생을 확인할 수 있는 알고리즘을 개발하였다. 또한 시멘트-모르타르 블록에서의 모형 발파시험을 수치해석을 수행함으로써, 개발된 알고리즘을 검증하였다.

• Smart Structures and Systems
• /
• 제22권6호
• /
• pp.643-662
• /
• 2018

One of the methods for investigation of mechanical behavior of materials is numerical simulation. For simulation, its need to model behavior is close to real condition. PFC is one of the rock mechanics software that needs calibration for models simulation. The calibration was performed based on simulation of unconfined compression test and Brazilian test. Indeed the micro parameter of models change so that the UCS and Brazilian test results in numerical simulation be close to experimental one. In this paper, the effect of four micro parameters has been investigated on the uniaxial compression test and Brazilian test. These micro parameters are friction angle, Accumulation factor, expansion coefficient and disc distance. The results show that these micro parameters affect the failure pattern in UCS and Brazilian test. Also compressive strength and tensile strength are controlled by failure pattern.

• Structural Engineering and Mechanics
• /
• 제69권1호
• /
• pp.43-50
• /
• 2019

In this research, the effect of bedding layer on the tensile failure mechanism of rocks has been investigated using PFC2D. For this purpose, firstly calibration of PFC2d was performed using Brazilian tensile strength. Secondly Brazilian test was performed on the bedding layer. Thickness of layers were 5 mm, 10 mm and 20 mm. in each thickness layer, layer angles changes from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Totally, 21 model were simulated and tested by loading rate of 0.016 mm/s. The results show that when layer angle is less than 15, tensile cracks initiates between the layers and propagate till coalesce with model boundary. Its trace is too high. With increasing the layer angle, less layer mobilizes in failure process. Also, the failure trace is very short. It's to be noted that number of cracks decrease with increasing the layer thickness. Also, Brazilian tensile strength is minimum when bedding layer angle is between $45^{\circ}$ and $75^{\circ}$. The maximum one is related to layer angle of $90^{\circ}$.

• 터널과지하공간
• /
• 제14권1호
• /
• pp.54-68
• /
• 2004

본 연구에서는 경제성장기에 고층 구조물 시공에서 널리 사용된 철근콘크리트(RC, Reinforced Concrete) 구조물을 대상으로 발파해체 축소모형실험을 수행하고 전산실험결과와 이를 비교하였다. 발파해체 공법으로는 파괴거동을 비교적 쉽게 확인할 수 있는 점진붕괴공법을 적용하였으며, 축소모형실험은Hobbs(1969)에 의한 축소율의 개념에 따라 차원해석을 실시하여 축소된 강도특성을 계산하였다. 사용재료로는 석고, 모래, 물의 혼합액을 콘크리트 대용으로 사용하였다. 모래와 석고의 중량 비를 다양하게 변화시키면서 이에 따른 강도의 변화를 측정하고 최적의 강도 값을 갖는 배합 비를 결정하였다. 또한 연성을 가지고 있으며 축소강도로 비교할 때 철근과 유사한 특성을 지니는 땜용 납선을 철근대용 재료로 사용하였다. 수치해석 프로그램으로는 요소의 파괴거동을 육안으로 확인할 수 있는 개별 요소법(DEM Distinct Element Method)에 의해 수행되는 상용코드인 PFC2D(Particle Flow Code 2-Dimension)을 사용하였다. 모형의 제작은 실내에서 미리 양생된 부재들을 현장으로 옮겨 연결부만을 타설하여 일체화시키는 방법으로 이루어졌다. 먼저 3차원 무근 콘크리트 라멘 구조의 모형을 설계하고 그 축소 모형을 발파해체하여 그 거동을 촬영하였다. 이를 수치적인 해석과 비교하는 과정을 통해 2차원 해석이라는 한계성은 존재하지만 대체로 유사한 형태의 거동을 보임을 알 수 있었다. 사전해석의 경험과 RC 보의 실내 굴곡 실험결과를 근거로 하여 RC 구조모형의 발파해체 사전해석을 실시하였다. 시차는 200㎳로 하여 점진적으로 붕괴되도록 설계하였다. 모형실험과는 달리 2차원 해석이라는 한계에도 불구하고 900㎳까지 매우 유사한 거동을 보이며 붕괴되었다.

• Geomechanics and Engineering
• /
• 제23권5호
• /
• pp.455-466
• /
• 2020

This paper presents analysis of the interaction between tunnel and Qanat with a particular interest for the optimization of Qanat shape using the discrete element code, PFC2D, and the results will be compared with the FEM results of PLAXIS2D. For these concerns, using software PFC2D based on Discrete Element Method (DEM), a model with dimension of 100m ^⁎ 100 m was prepared. A circular tunnel with dimension of 9 m was situated 20 m below the ground surface. Also one Qanat was situated perpendicularly above the tunnel roof. Distance between Qanat center and ground surface was 8 m. Five different shapes for Qanat were selected i.e., square, semi-circular, vertical ellipse, circular and horizontal ellipse. Confining pressure of 5 MPa was applied to the model. The vertical displacement of balls situated in ground surface was picked up to measure the ground subsidence. Also two measuring circles were situated at the tunnel roof and at the Qanat roof to check the vertical displacements. The properties of the alluvial soil of Tehran city are: γ_dry=19 (KN/㎥), E= 750 (kg/㎠), ν=0.35, c=0.3(kg/㎠), φ=34°. In order to validate the DEM results, a comparison between the numerical results (obtained in this study) and analytical and field monitoring have been done. The PFC2D results are compared with the FEM results. The results shows that when Qanat has rectangular shape, the tensile stress concentration at the Qanat corners has maximum value while it has minimum value for vertical ellipse shape. The ground subsidence for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. The vertical displacements at the tunnel roof for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. Historical shape of Qante approved the finding of this research.

• Geomechanics and Engineering
• /
• 제13권6호
• /
• pp.1045-1055
• /
• 2017

Tensile strength is considered key properties for characterizing rock material in engineering project. It is determined by direct and indirect methods. Point load test is a useful testing method to estimate the tensile strengths of rocks. In this paper, the effects of rock shape on the point load index of gypsum are investigated by PFC2D simulation. For PFC simulating, initially calibration of PFC was performed with respect to the Brazilian experimental data to ensure the conformity of the simulated numerical models response. In second step, nineteen models with different shape were prepared and tested under point load test. According to the obtained results, as the size of the models increases, the point load strength index increases. It is also found that the shape of particles has no major effect on its tensile strength. Our findings show that the dominant failure pattern for numerical models is breaking the model into two pieces. Also a criterion was rendered numerically for determination of tensile strength of gypsum. The proposed criteria were cross checked with the results of experimental point load test.

• Structural Engineering and Mechanics
• /
• 제85권3호
• /
• pp.337-351
• /
• 2023

In this paper, the shear behavior of soft filling in rectangular-hollow concrete specimens was simulated using the 2D particle flow code (PFC2D). The laboratory-measured properties were used to calibrate some PFC2D micro-properties for modeling the behavior of geo-materials. The dimensions of prepared and modeled samples were 100 mm×100 mm. Some disc type narrow bands were removed from the central part of the model and different lengths of bridge areas (i.e., the distance between internal tips of two joints) with lengths of 30 mm, 50 mm, and 70 mm were produced. Then, the middle of the rectangular hollow was filled with cement material. Three filling sizes with dimensions of 5 mm×5 mm, 10 mm×5 mm, and 15 mm×5 mm were provided for different modeled samples. The parallel bond model was used to calibrate and re-produce these modeled specimens. Therefore, totally, 9 different types of samples were designed for the shear tests in PFC2D. The shear load was gradually applied to the model under a constant loading condition of 3 MPa (σc/3). The loading was continued till shear failure occur in the modeled concrete specimens. It has been shown that both tensile and shear cracks may occur in the fillings. The shear cracks mainly initiated from the crack (joint) tips and coalesced with another one. The shear displacements and shear strengths were both increased as the filling dimensions increased (for the case of a bridge area with a particular fixed length).

• 한국터널지하공간학회 논문집
• /
• 제14권3호
• /
• pp.261-275
• /
• 2012

본 논문은 이수식 쉴드 TBM 굴착시 굴착면 안정에 큰 영향을 미치는 굴착면의 이수 침투에 따른 이막형성여부에 대하여 해석적 고찰을 수행하였다. 이러한 해석적 접근은 이수식 쉴드 TBM의 계획 및 설계시 최적 장비 형식 선정 및 시공계획 수립에 참고가 될 수 있으나, 이제까지 이에 대한 효과적인 해석적 방법에 대하여 제안된 것이 없는 상황이다. 따라서 본 논문에서는 사질토 입자의 모델링에 효과적인 $PFC^{2D}$의 Fluid Coupling simulation을 이용하여 해석을 수행하였으며, 해석결과, 이수침투 특성은 굴착 대상지반의 투수성과 이수특성(비중, 점성계수등)에 따라 변화됨을 파악할 수 있었다. 본 해석적 접근방법을 통하여 이수식 쉴드 TBM의 계획 설계단계에서 대상지반 특성을 고려한 이수 점도관리 방안 수립과 함께 시공중의 초기굴진단계에서의 최적 이수점도의 재검증을 통하여 안전한 터널시공에 도움이 될 것이라 판단된다.

• 터널과지하공간
• /
• 제15권2호
• /
• pp.119-128
• /
• 2005

절리 암반의 역학적 물성 및 거동 평가가 터널 및 지하구조물의 설계에 매우 중요하다 할지라도, 그것은 항상 매우 어려운 문제로 간주되어 왔다. 암반 거동을 모사하는데 있어서 어려움중의 하나는 적절한 구성 모델을 선정하는 것이다. 이러한 한계점은 PFC와 같이 사용자로 하여금 암반의 구성 모델을 요구하지 않는 개별요소 프로그램의 개발과 함께 극복되어질 것이다. 본 연구에서는 도로터널 현장의 30\;m\;\times\;30\;m\;\times\;30\;m 절리 암반블록을 대상으로, 시추 및 지표 지질조사를 통해 얻어진 절리의 기하학적 형태자료를 근거로 개별균열망이 작성되었다. 개별균열망 모델의 절리 형상을 근거로 절리가 없는 상태에서 점차적으로 절리군을 추가해가면서 2차원 PFC모델이 만들어졌다. 또한 각각의 PFC모델에 대한 수치모사를 통하여 각 모델의 응력-변형율 곡선이 얻어졌다. 응력-변형율 곡선으로부터 절리 암반의 역학적 물성이 결정되었다. 절리의 존재는 암반의 역학적 물성에 상당한 영향을 미쳤으며, 더욱 중요한 것은 PFC모델의 역학적 거동은 기존의 수치모델에서 요구되는 구성 모델에 의하여 결정되지 않는다는 것이다.

• Computers and Concrete
• /
• 제31권3호
• /
• pp.207-221
• /
• 2023

In this investigation, the interaction between opening space and neighboring joint has been examined by experimental test and Particle flow code in two dimension (PFC2D) simulation. Since, firs of all PFC was calibrated using Brazilian experimental test and uniaxial compression test. Secondly, diverse configurations of opening and neighboring joint were provided and tested by uniaxial test. 12 rectangular sample with dimension of 10 cm*10 cm was prepared from gypsum mixture. One quarter of tunnel and one and or two joint were drilled into the sample. Tunnel diameter was 5.5 cm. The angularities of joint in physical test were 0°, 45° and 90°. The angularities of joint in numerical simulation were 0°, 30°, 60°, -30°, -45°, -60° and its length were 2cm and 4cm. Loading rate was 0.016 m/s. Tensile strength of material was 4.5 MPa. Results shows that dominant type of crack which took place in the model was tensile cracks and or several shear bands develop within the model. The Final stress is minimum in the cases where oriented angle is negative. The failure stress decrease by decreasing the joint angle from 30° to 60°. In addition, the failure stress decrease by incrementing the joint angle from -30° to -60°. The failure stress was incremented by decreasing the number of notches. The failure stress was incremented by decreasing the joint length. The failure stress was incremented by decreasing the number of notches. Comparing experimental results and numerical one, showed that the failure stress is approximately identical in both conditions.