Analysis of stress distribution around tunnels by hybridized FSM and DDM considering the influences of joints parameters

Nikadat, Nooraddin,Marji, Mohammad Fatehi Techno-Press 2016 Geomechanics & engineering Vol.11 No.2

The jointed rock mass behavior often plays a major role in the design of underground excavation, and their failures during excavation and in operation, are usually closely related to joints. This research attempts to evaluate the effects of two basic geometric factors influencing tunnel behavior in a jointed rock mass; joints spacing and joints orientation. A hybridized indirect boundary element code known as TFSDDM (Two-dimensional Fictitious Stress Displacement Discontinuity Method) is used to study the stress distribution around the tunnels excavated in jointed rock masses. This numerical analysis revealed that both the dip angle and spacing of joints have important influences on stress distribution on tunnel walls. For example the tensile and compressive tangential stresses at the boundary of the circular tunnel increase by reduction in the joint spacing, and by increase the dip joint angle the tensile stress in the tunnel roof decreases.

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

Jinwei Fu,Vahab Sarfarazi,Mohammad Fatehi Marji,Mengdi Guo 국제구조공학회 2022 Structural Engineering and Mechanics, An Int'l Jou Vol.81 No.1

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

Numerical simulation of the effect of bedding layer geometrical properties on the shear failure mechanism using PFC3D

Haeri, Hadi,Sarfarazi, Vahab,Zhu, Zheming,Marji, Mohammad Fatehi Techno-Press 2018 Smart Structures and Systems, An International Jou Vol.22 No.5

In this research the effect of bedding layer angle and bedding layer thickness on the shear failure mechanism of concrete has been investigated using PFC3D. For this purpose, firstly calibration of PFC3d was performed using Brazilian tensile strength. Secondly shear 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 $25^{\circ}$. Totally 15 model were simulated and tested by loading rate of 0.016 mm/s. The results shows that when layer angle is less than $50^{\circ}$, 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 mobilize in failure process. Also the failure trace is very short. It's to be note that number of cracks decrease with increasing the layer thickness. The minimum shear test strength was occurred when layer angle is more than $50^{\circ}$. The maximum value occurred in $0^{\circ}$. Also, the shear test tensile strength was increased by increasing the layer thickness.

Numerical simulation of the effect of confining pressure and tunnel depth on the vertical settlement using particle flow code (with direct tensile strength calibration in PFC Modeling)

Hadi Haeri,Vahab Sarfarazi,Mohammad Fatehi Marji 국제구조공학회 2020 Smart Structures and Systems, An International Jou Vol.25 No.4

In this paper the effect of confining pressure and tunnel depth on the ground vertical settlement has been investigated using particle flow code (PFC2D). For this perpuse firstly calibration of PFC2D was performed using both of tensile test and triaxial test. Then a model with dimention of 100 m × 100 m was built. Acircular tunnel with diameter of 20 m was drillled in the middle of the model. Also, a rectangular tunnel with wide of 10 m and length of 20 m was drilled in the model. The center of tunnel was situated 15 m, 20 m, 25 m, 30 m, 35 m, 40 m, 45 m, 50 m, 55 m and 60 m below the ground surface. these models are under confining pressure of 0.001 GPa, 0.005 GPa, 0.01 GPa, 0.03 GPa, 0.05 GPa and 0.07 GPa. The results show that the volume of colapce zone is constant by increasing the distance between ground surface and tunnel position. Also, the volume of colapce zone was increased by decreasing of confining pressure. The maximum of settlement occurs at the top of the tunnel roof. The maximum of settlement occurs when center of tunnel was situated 15 m below the ground surface. The settlement decreases by increasing the distance between tunnel center line and measuring circles in the ground surface. The minimum of settlement occurs when center of circular tunnel was situated 60 m below the surface ground. Its to be note that the settlement increase by decreasing the confining pressure.

Numerical simulation of the effect of bedding layer on the tensile failure mechanism of rock using PFC2D

Vahab Sarfarazi,Hadi Haeri,Mohammad Fatehi Marji 국제구조공학회 2019 Structural Engineering and Mechanics, An Int'l Jou Vol.69 No.1

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° to 90° with increment of 15°. 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° and 75°. The maximum one is related to layer angle of 90°.

Direct tensile strength measurement of granite by the universal tensile testing machine

Hadi Haeri,Vahab Sarfarazi,Mohammad Fdatehi Marji,Mohammad Davood Yavari,Amin Zahedi-khameneh 국제구조공학회 2021 Smart Structures and Systems, An International Jou Vol.27 No.4

The direct tensile strength of a typical hard rock like granite is measured by a novel apparatus known as compression-to-tensile load transfer (CTLT) device. The rock specimen is prepared in form of a slab containing a central hole and 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, some typical hard rock specimens of granite are specially prepared and tested in the laboratory to measure their direct tensile strengths. Then, a new load converting device implemented in the universal tensile testing machine is used to cause the rock specimen to be subjected to a direct tensile loading during the test. The compressive load was applied to the transferring device at the rate of 0.02 MPa/s. Numerical modeling of the tested specimens were accomplished using the discrete element method (DEM) and the higher order displacement discontinuity method (HODDM). The tensile failure of granite rock mainly occurs along the horizontal axis. The experimental results were in a good accordance with DEM results and HODDM outputs.

Simulation of the tensile failure behaviour of transversally bedding layers using PFC2D

Haeri, Hadi,Sarfarazi, Vahab,Zhu, Zheming,Marji, Mohammad Fatehi Techno-Press 2018 Structural Engineering and Mechanics, An Int'l Jou Vol.67 No.5

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.

The effect of particle size on the edge notched disk (END) using particle flow code in three dimension

Hadi Haeri,Vahab Sarfarazi,Zheming Zhu,Mohammad Fatehi Marji 국제구조공학회 2018 Smart Structures and Systems, An International Jou Vol.22 No.6

In this study, the effect of particle size on the cracks propagation and coalescence or cracking pattern of the edge notched disc specimens are investigated. Firstly, calibration of PFC3D was performed using Brazilian experimental test output. Then micro parameters were used to build edge notched disc specimen. The horizontal wall of the assembly is let to move downward with a standard low speed of 0.016 m/s. The numerical results show that the tensile cracks are dominant failure pattern for the modeled discs. These tensile cracks initiate from the pre-existing notch tip and propagate parallel to the loading direction then interact with the upper boundary of the modeled specimen. As the size of the balls (ball diameter) decrease the number of tensile cracks increase. The tensile fracture toughness of the samples also decreases as the particle size increases. Understanding the crack propagation and crack coalescence phenomena in brittle materials such as concretes and rocks is of paramount importance in the stability analyses for engineering structures such as rock slopes, underground structures and tunneling.

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

Hadi Haeri,Vahab Sarfaraz,Zheming Zhu,Mohammad Fatehi Marji 국제구조공학회 2019 Smart Structures and Systems, An International Jou Vol.23 No.5

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 *10 cm* 5 cm. the length of joint is 2 cm. the angularities of joint are 30°, 45° and 60°. For each joint angularity, micro pore was situated 2cm 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°, 45° and 60°. 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.

Numerical simulation of compressive to tensile load conversion for determining the tensile strength of ultra-high performance concrete

Hadi Haeri,Nader Mirshekari,Vahab Sarfarazi,Mohammad Fatehi Marji 국제구조공학회 2020 Smart Structures and Systems, An International Jou Vol.26 No.5

In this study, the experimental tests for the direct tensile strength measurement of Ultra-High Performance Concrete (UHPC) were numerically modeled by using the discrete element method (circle type element) and Finite Element Method (FEM). The experimental tests used for the laboratory tensile strength measurement is the Compressive-to-Tensile Load Conversion (CTLC) device. In this paper, the failure process including the cracks initiation, propagation and coalescence studied and then the direct tensile strength of the UHPC specimens measured by the novel apparatus i.e., CTLC device. For this purpose, the UHPC member (each containing a central hole) prepared, and situated in the CTLC device which in turn placed in the universal testing machine. The direct tensile strength of the member is measured due to the direct tensile stress which is applied to this specimen by the CTLC device. This novel device transferring the applied compressive load to that of the tensile during the testing process. The UHPC beam specimen of size 150 × 60 × 190 mm and internal hole of 75 × 60 mm was used in this study. The rate of the applied compressive load to CTLC device through the universal testing machine was 0.02 MPa/s. The direct tensile strength of UHPC was found using a new formula based on the present analyses. The numerical simulation given in this study gives the tensile strength and failure behavior of the UHPC very close to those obtained experimentally by the CTLC device implemented in the universal testing machine. The percent variation between experimental results and numerical results was found as nearly 2%. PFC2D simulations of the direct tensile strength measuring specimen and ABAQUS simulation of the tested CTLC specimens both demonstrate the validity and capability of the proposed testing procedure for the direct tensile strength measurement of UHPC specimens.