• Journal of the Korean GEO-environmental Society
• /
• v.12 no.7
• /
• pp.57-65
• /
• 2011

Numerical simulations by three-dimensional Particle Flow Code($PFC^{3D}$, Itasca) considering distinct element method (DEM) were carried out for prediction of triaxial compression test with sand material. The effect of scale conditions for numerical model and distinct material on final prediction results was analyzed by numerical models under various scale conditions, and following observations were made from the numerical experiments. It is very useful to model the initial material condition without any porosity conversion from 2-D to 3-D DEM. Numerical experiments have shown that in all cases considered, 3D distinct element modeling could provide good agreement on stress-strain behavior, volume change and strength properties with laboratory testing results. It was important thing to assess reasonable scale ratio of numerical model and distinct elements for saving calculation time and securing calculation efficiency under condition with accuracy and appropriateness as numerical laboratory. As results of DEM simulations under various scale conditions, most of results show that shear strength properties as cohesion and internal friction angle are similar in condition of $D_{mod}/D_{gmax}$ < 10. It shows that 3-D distinct element method could be used as efficient tool to assess strength properties by numerical laboratory technique.

• Structural Engineering and Mechanics
• /
• v.67 no.5
• /
• pp.493-504
• /
• 2018

In this paper, the tensile failure behaviour of transversally bedding layers was numerically simulated by using particle flow code in two dimensions. Firstly, numerical model was calibrated by uniaxial, Brazilian and triaxial experimental results to ensure the conformity of the simulated numerical model's response. Secondly, 21 circular models with diameter of 54 mm were built. Each model contains two transversely bedding layers. The first bedding layer has low mechanical properties, less than mechanical properties of intact material, and second bedding layer has high mechanical properties, more than mechanical properties of intact material. The angle of first bedding layer, with weak mechanical properties, related to loading direction was $0^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, $75^{\circ}$ and $90^{\circ}$ while the angle of second layer, with high mechanical properties, related to loading direction was $90^{\circ}$, $105^{\circ}$, $120^{\circ}$, $135^{\circ}$, $150^{\circ}$, $160^{\circ}$ and $180^{\circ}$. Is to be note that the angle between bedding layer was $90^{\circ}$ in all bedding configurations. Also, three different pairs of the thickness was chosen in models; i.e., 5 mm/10 mm, 10 mm/10 mm and 20 mm/10 mm. The result shows that In all configurations, shear cracks develop between the weaker bedding layers. Shear cracks angel related to normal load change from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Numbers of shear cracks are constant by increasing the bedding thickness. It's to be note that in some configuration, tensile cracks develop through the intact area of material model. There is not any failure in direction of bedding plane interface with higher strength.

• Geomechanics and Engineering
• /
• v.30 no.1
• /
• pp.93-105
• /
• 2022

With the exploitation of natural resources in China, underground resource extraction and underground space development, as well as other engineering activities are increasing, resulting in the creation of many defective rocks. In this paper, uniaxial compression tests were performed on rocks with double holes and fractures at different angles using particle flow code (PFC2D) numerical simulations and laboratory experiments. The failure behavior and mechanical properties of rock samples with holes and fractures at different angles were analyzed. The failure modes of rock with defects at different angles were identified. The fracture propagation and stress evolution characteristics of rock with fractures at different angles were determined. The results reveal that compared to intact rocks, the peak stress, elastic modulus, peak strain, initiation stress, and damage stress of fractured rocks with different fracture angles around holes are lower. As the fracture angle increases, the gap in mechanical properties between the defective rock and the intact rock gradually decreased. In the force chain diagram, the compressive stress concentration range of the combined defect of cracks and holes starts to decrease, and the model is gradually destroyed as the tensile stress range gradually increases. When the peak stress is reached, the acoustic emission energy is highest and the rock undergoes brittle damage. Through a comparative study using laboratory tests, the results of laboratory real rocks and numerical simulation experiments were verified and the macroscopic failure characteristics of the real and simulated rocks were determined to be similar. This study can help us correctly understand the mechanical properties of rocks with defects and provide theoretical guidance for practical rock engineering.

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

• 한국신재생에너지학회:학술대회논문집
• /
• 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의 열전도도는 접촉된 두 입자의 열전도도와 접촉면의 평균 열전도도를 고려하여 적용하였다.

• Geomechanics and Engineering
• /
• v.8 no.6
• /
• pp.757-767
• /
• 2015

Rock mass is an important engineering material. In hydropower engineering, rock mass of bank slope controlled the stability of an arch dam. However, mechanical characteristics of the rock mass are not only affected by lithology, but also joints. On the basis of field geological survey, this paper built rock mass material containing parallel concentrated joints with different dip angle, different number under different stress conditions by PFC (Particle Flow Code) numerical simulation. Next, we analyzed mechanical property and fracture features of this rock mass. The following achievements have been obtained through this research. (1) When dip angle of joints is $15^{\circ}$ and $30^{\circ}$, with the increase of joints number, peak strength of rock mass has not changed much. But when dip angle increase to $45^{\circ}$, especially increase to $60^{\circ}$ and $75^{\circ}$, peak strength of rock mass decreased obviously with the increase of joints number. (2) With the increase of confining stress, peak strengths of all rock mass have different degree of improvement, especially the rock mass with dip angle of $75^{\circ}$. (3) Under the condition of no confining stress, dip angle of joints is low and joint number is small, existence of joints has little influence on fracture mode of rock mass, but when joints number increase to 5, tensile deformation firstly happened at joints zone and further resulted in tension fracture of the whole rock mass. When dip angle of joints increases to $45^{\circ}$, fracture presented as shear along joints, and with increase of joints number, strength of rock mass is weakened caused by shear-tension fracture zone along joints. When dip angle of joints increases to $60^{\circ}$ and $75^{\circ}$, deformation and fracture model presented as tension fracture zone along concentrated joints. (4) Influence of increase of confining stress on fracture modes is to weaken joints' control function and to reduce the width of fracture zone. Furthermore, increase of confining stress translated deformation mode from tension to shear.

• Geomechanics and Engineering
• /
• v.15 no.4
• /
• pp.1017-1027
• /
• 2018

To explore the influence of coal thickness on the mechanical behavior and the failure characteristics of rock-coal-rock (RCR) mass, the experimental investigation of uniaxial compressive tests was conducted first and then a systematic numerical simulation by particle flow code (PFC2D) was performed to deeply analyze the failure mechanical behavior of RCR specimens with different coal thicknesses in conventional compression tests. The overall elastic modulus and peak stress of RCR specimens lie between the rock and the coal. Inter-particle properties were calibrated to match the physical sample strength and the stiffness response. Numerical simulation results show that the deformation and strength behaviors of RCR specimens depend not only on the coal thickness, but also on the confining pressure. Under low confining pressures, the overall failure mechanism of RCR specimen is the serious damage of coal section when the coal thickness is smaller than 30 mm, but it is shear failure of coal section when the coal thickness is larger than 30 mm. Whereas under high confining pressures, obvious shear bands exist in both the coal section and the rock section when the coal thickness is larger than 30 mm, but when the coal thickness is smaller than 30mm, the failure mechanism is serious damage of coal section and shear failure of rock section.

• Geomechanics and Engineering
• /
• v.22 no.5
• /
• pp.385-396
• /
• 2020

A uniaxial compression test was performed to analyse the mechanical properties and macroscale and mesoscale failure mechanisms of sandstone with pyrite concretions. The effect of the pyrite concretions on the evolution of macroscale cracks in the sandstone was further investigated through numerical simulations with Particle Flow Code in 2D (PFC2D). The results revealed that pyrite concretions substantially influence the mechanical properties and macroscale and mesoscale failure characteristics of sandstone. During the initial loading stage, significant stress concentrations occurred around the edges of the pyrite concretion accompanied by the preferential generation of cracks. Meanwhile, the events and cumulative energy counts of the acoustic emission (AE) signal increased rapidly because of friction sliding between the concretion and sandstone matrix. As the axial stress increased, the degree of the stress concentration remained relatively unchanged around the edges of the concretions. The cracks continued growing rapidly around the edges of the concretions and gradually expanded toward the centre of the sample. During this stage, the AE events and cumulative energy counts increased quite slowly. As the axial stress approached the peak strength of the sandstone, the cracks that developed around the edges of the concretion started to merge with cracks that propagated at the top-left and bottom-right corners of the sample. This crack evolution ultimately resulted in the shear failure of the sandstone sample around the edges of the pyrite concretions.

• Advances in concrete construction
• /
• v.13 no.2
• /
• pp.199-213
• /
• 2022

Experimental and discrete element approaches were used to examine the effects of F shape non-persistent joints on the failure behaviour of concrete under uniaxial compressive test. concrete specimens with dimensions of 200 cm×200 cm×50 cm were provided. Within the specimen, F shape non-persistent joint consisting three joints were provided. The large joint length was 6 cm, and the length of two small joints were 2 cm. Vertical distance between two small joints change from 1.5 cm to 4.5 cm with increment of 1.5 cm. In constant joint lengths, the angle of large joint change from 0° to 90° with increments of 30°. Totally 12 different models were tested under compression test. The axial load rate on the model was 0.05 mm/min. Concurrent with experimental tests, numerical simulation (Particle flow code in two dimension) were performed on the models containing F shape non-persistent joint. Distance between small joints and joint angles were similar to experimental one. the results indicated that the failure process was mostly governed by both of the Distance between small joints and joint angles. The axial loading rate on the model was 0.05 mm/min. The compressive strengths of the samples were related to the fracture pattern and failure mechanism of the discontinuities. Furthermore, it was shown that the compressive behaviour of discontinuities is related to the number of the induced tensile cracks which are increased by increasing the joint angle. In the first, there were only a few acoustic emission (AE) hits in the initial stage of loading, and then AE hits rapidly grow before the applied stress reached its peak. Furthermore, a large number of AE hits accompanied every stress drop. Finally, the failure pattern and failure strength are similar in both approaches i.e., the experimental testing and the numerical simulation approaches.

• Structural Engineering and Mechanics
• /
• v.81 no.2
• /
• pp.243-258
• /
• 2022

In this paper, tensile behavior of joint filling has been investigated under experimental test and numerical simulation (particle flow code). Two concrete slabs containing semi cylinder hole were prepared. These slabs were attached to each other by glue and one cubic specimen with dimension of 19 cm×15 cm×6 cm was prepared. This sample placed in the universal testing machine where the direct tensile stress can be applied to this specimen by implementing a special type of load transferring device which converts the applied compressive load to that of the tensile during the test. In the present work, two different joint filling thickness i.e., 3 mm and 6 mm were prepared and tested in the laboratory to measure their direct tensile strengths. Concurrent with experimental test, numerical simulation was performed to investigate the effect of hole diameter, length of edge notch, filling thickness and filling length on the tensile behavior of joint filling. Model dimension was 19 cm×15 cm. hole diameter was change in four different values of 2.5 cm, 5 cm, 7.5 cm and 10 cm. glue lengths were different based on the hole diameter, i.e., 12.5 cm for hole diameter of 2.5 cm, 10 cm for hole diameter of 5 cm, 7.5 cm for hole diameter of 7.5 cm and 5 cm for hole diameter of 10 cm. length of edge notch were changed in three different value i.e., 10%, 30% and 50% of glue length. Filling thickness were changed in three different value of 3 mm, 6 mm and 9 mm. Tensile strengths of glue and concrete were 2.37 MPa and 6.4 MPa, respectively. The load was applied at a constant rate of 1 kg/s. Results shows that hole diameter, length of edge notch, filling thickness and filling length have important effect on the tensile behavior of joint filling. In fixed glue thinks and fixed joint length, the tensile strength was decreased by increasing the hole diameter. Comparing the results showed that the strength, failure mechanism and fracture patterns obtained numerically and experimentally were similar for both cases.