• Journal of Korean Tunnelling and Underground Space Association
• /
• v.12 no.1
• /
• pp.31-41
• /
• 2010

This paper presents both numerical and physical modeling approaches for the dilation and spalling of rock recognized as typical process of rock around an underground opening at depth. For physical approach, laboratory testing of rectangular beams using a synthetic rock was used to investigate the onset of dilation and spalling. The beams are axially compressed and subjected to 4-point bending to provide non-uniform compressive stresses which are similar to the maximum tangential stress distribution around circular openings. Discrete element numerical analyses using commercial code $PFC^{2D}$ (Particle Flow Code) were performed to evaluate the stress path at various locations in the beams. The findings from these approaches suggest that the onset of dilation in laboratory tests appears to be a good indicator for assessing the stress magnitudes required to initiate spalling.

• Explosives and Blasting
• /
• v.22 no.1
• /
• pp.33-43
• /
• 2004

In this study, a comparison was made between the resulting behaviors of scaled model test and particle flow analysis for blasting demolition of a reinforced concrete structure. For the test and analysis, a progressive failure of a five-story structure was considered. The dimension analysis was carried out to properly scale down the real structure into the laboratory size. The test model was made of the mixture of gypsum, sand and water along with soldering lead to analogy reinforcing steel bars. The ratio of mixing components was chosen to best represent the scaled down strength and deformation modulus. The columns and girders of the structure were precasted in the laboratory and assembled right before the blasting test. The numerical analysis of the blasting demolition was carried out using PFC2D (Particle Flow Analysis 2-Dimension by Itasca). The results of the blasting of concrete lahmen structure showed roughly identical demolition behavior between scaled model test and numerical test. For the blasting of the reinforced concrete structure, the results were more identical and closer to the real demolition behavior, since the demolition behavior was better represented in this case due to the increased tensile strength of the component.

• Computers and Concrete
• /
• v.31 no.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.

• Structural Engineering and Mechanics
• /
• v.69 no.1
• /
• pp.11-20
• /
• 2019

The biaxial failure mechanism of transversally bedding concrete layers was numerically simulated using a sophisticated two-dimensional discrete element method (DEM) implemented in the particle flow code (PFC2D). This numerical modelling code was first calibrated by uniaxial compression and Brazilian testing results to ensure the conformity of the simulated numerical model's response. Secondly, 21 rectangular models with dimension of $54mm{\times}108mm$ 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 were 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 noted 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.

• Structural Engineering and Mechanics
• /
• v.86 no.5
• /
• pp.607-619
• /
• 2023

This study, it was tried to evaluate the asphalt behavior under tensile loading conditions through indirect Brazilian and direct tensile tests, experimentally and numerically. This paper is important from two points of view. The first one, a new test method was developed for the determination of the direct tensile strength of asphalt and its difference was obtained from the indirect test method. The second one, the effects of particle size and loading rate have been cleared on the tensile fracture mechanism. The experimental direct tensile strength of the asphalt specimens was measured in the laboratory using the compression-to-tensile load converting (CTLC) device. Some special types of asphalt specimens were prepared in the form of slabs with a central hole. The CTLC device is then equipped with this specimen and placed in the universal testing machine. Then, the direct tensile strength of asphalt specimens with different sizes of ingredients can be measured at different loading rates in the laboratory. The particle flow code (PFC) was used to numerically simulate the direct tensile strength test of asphalt samples. This numerical modeling technique is based on the versatile discrete element method (DEM). Three different particle diameters were chosen and were tested under three different loading rates. The results show that when the loading rate was 0.016 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis till coalescence to the model boundary. When the loading rate was 0.032 mm/sec, two tensile cracks were initiated from the left and right of the hole and propagated perpendicular to the loading axis. The branching occurs in these cracks. This shows that the crack propagation is under quasi-static conditions. When the loading rate was 0.064 mm/sec, mixed tensile and shear cracks were initiated below the loading walls and branching occurred in these cracks. This shows that the crack propagation is under dynamic conditions. The loading rate increases and the tensile strength increases. Because all defects mobilized under a low loading rate and this led to decreasing the tensile strength. The experimental results for the direct tensile strengths of asphalt specimens of different ingredients were in good accordance with their corresponding results approximated by DEM software.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2010.03a
• /
• pp.314-323
• /
• 2010

Recently, granular soils having a large particle size are frequently used as a filling material in the construction of foundation, harbor, dam, and so on. The shear behavior of this granular soil plays a key role in the stability of structures. For example, soil particle crushing occurring at the interface between structure and soil and/or within soil mass can cause the disturbance of ground characteristics and consequently induce an issues in respect of stability of structures. In order to investigate the shear behavior according to an existence and nonexistence of particle crushing, numerical analyses were conducted by using the DEM(Discrete Element Method)-based software program PFC(Particle Flow Code). Using the crushing model and non-crushing model which were created in this study, numerical analyses of ring shear test were conducted and their results were analyzed and compared. In general, landslide and slope stability are accompanied by a large displacement and consequently not only a peak strength but also a residual strength are very important in the analysis of landslide and slope stability. However the direct shear test which has been commonly used in the determination of shear strength parameters has a limitation on displacement therefore the residual strength parameters can not be obtained. The characteristics of residual shear behavior were investigated through the numerical analyses in this study.

• Computers and Concrete
• /
• v.20 no.1
• /
• pp.39-47
• /
• 2017

In this paper, mode I fracture toughness of rock was determined by direct and indirect methods using Particle Flow Code simulation. Direct methods are compaction tension (CT) test and hollow centre cracked quadratic sample (HCCQS). Indirect methods are notched Brazilian disk (NBD) specimen, the semi-circular bend (SCB) specimen, hollow centre cracked disc (HCCD), the single edge-notched round bar in bending (SENRBB) specimen and edge notched disk (END). It was determined that which one of indirect fracture toughness values is close to direct one. For this purpose, initially calibration of PFC was undertaken with respect to data obtained from Brazilian laboratory tests to ensure the conformity of the simulated numerical models response. Furthermore, the simulated models in five introduced indirect tests were cross checked with the results from direct tests. By using numerical testing, the failure process was visually observed. Discrete element simulations demonstrated that the macro fractures in models are caused by microscopic tensile breakages on large numbers of bonded discs. Mode I fracture toughness of rock in direct test was less than other tests results. Fracture toughness resulted from semi-circular bend specimen test was close to direct test results. Therefore semi-circular bend specimen can be a proper test for determination of Mode I fracture toughness of rock in absence of direct test.

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

• Journal of the Korea Concrete Institute
• /
• v.20 no.2
• /
• pp.147-156
• /
• 2008

In an uniaxial compressive test of a concrete standard specimen (150$\times$300 mm) the crack initiation and extension with the stress increase are the major reason of the failure, which is similar to the breakage of the particle bonding in the simulation by using particle bonded model, especially particle flow code in 3 dimensions (PFC3D) developed by Itasca Consulting Group Inc. That is the main motive to study the possibility of an uniaxial compressive strength test simulation. It is important to investigate the relationship between the micro-parameters and the macro-properties because the 3-dimensional particle bonded model uses the spherical particles to analyze the physical phenomena. Contact bonded model used herein has eight micro-parameters and there are five macro-properties; Young's modulus, Poisson's ratio, uniaxial compressive strength and the crack initiation stress and the ratio concerning the crack propagation with the stress. To simulate the compressive test we made quantitative relationships between the micro-parameters and the macro-properties by using the fractional factorial design and various sensitivity analyses including regression analysis, which result in the good agreement with the previous studies. Also, the stress-stain curve and the crack distribution over the specimen given by PFC3D showed the mechanical behavior of the concrete standard specimen under the uniaxial compression. It is concluded that the particle bonded model can be a good tool for the analyzing the mechanical behavior of concrete under the uniaxial compressive load.

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