Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass...Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.展开更多
In this paper,Brazilian test was performed on disk samples of analogue materials with defined structural planes.The surface strain evolution process of the disk samples during loading was analyzed via digital image co...In this paper,Brazilian test was performed on disk samples of analogue materials with defined structural planes.The surface strain evolution process of the disk samples during loading was analyzed via digital image correlation.The damage evolution process was explored from a microscopic perspective by combining discrete element numerical simulation technology.The criterion of the failure mode of the disc specimen in the split state was theoretically deduced.The influence of structural surface roughness and loading inclination angle on the stress state at the center of the specimen was explored.The results showed that the failure modes of the samples could be divided into three typical modes as matrix failure,structural plane failure and combination failure.The rough structural plane improves the failure strength of the specimen by limiting its lateral deformation,and the degree of improvement weakens continuously with the increase of the inclination angle of the structural plane.As the inclination angle of the structural plane increases,the main type of microcracks in the structural plane changes from shear microcracks to tensile microcracks.This study contributes to a better understanding of macro-and meso-failure characteristics of rock masses with structural planes under a splitting state.展开更多
Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode o...Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode of the vertical well hydraulic fracturing,the true triaxial hydraulic fracturing test and numerical simulation are carried out,and the influence of the principal stress difference,water injection displacement,perforation angle and natural fracture on fracture propagation is analyzed.The results show that the fracture propagation mode of limestone is mainly divided into two types:the single vertical fracture and the transverse-longitudinal crossed complex fracture.Under high displacement,the fracturing pressure is larger,and the secondary fracture is more likely to occur,while variable displacement loading is more likely to induce fracture network.Meanwhile,the amplitude of acoustic emission(AE)waveform of limestone during fracturing is between 0.01 and 0.02 mV,and the main frequency is maintained in the range of 230−300 kHz.When perforation angleθ=45°,it is easy to produce the T-type fracture that connects with the natural fracture,while X-type cracks are generated whenθ=30°.The results can be used as a reference for further study on the mechanism of limestone hydraulic fracturing.展开更多
It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the mi...It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.展开更多
In this paper,a new micro-creep model of salt rock is proposed based on a linear parallel bonded model(LPBM)using the two-dimensional particle flow code(PFC2D).The power function weakening form is assumed to describe ...In this paper,a new micro-creep model of salt rock is proposed based on a linear parallel bonded model(LPBM)using the two-dimensional particle flow code(PFC2D).The power function weakening form is assumed to describe the variation of the parallel bonded diameter(PBD)over time.By comparing with the parallel-bonded stress corrosion(PSC)model,a smaller stress fluctuation and smoother creep strain−time curves can be obtained by this power function model at the same stress level.The validity and adaptability of the model to simulate creep deformation of salt rock are verified through comparing the laboratory creep test curves and the Burgers model fitting result.The numerical results reveal that this model can be capable of capturing the creep deformation and damage behavior from the laboratory observations.展开更多
The mechanical behaviors of the interface between coarse-grained soil and concrete were investigated by simple shear tests under condition of mixed soil slurry (bentonite mixed with cement grout).For comparison,the in...The mechanical behaviors of the interface between coarse-grained soil and concrete were investigated by simple shear tests under condition of mixed soil slurry (bentonite mixed with cement grout).For comparison,the interfaces both without slurry and with bentonite slurry were analyzed.The experimental results show that different slurries exert much influence on the strength and deformation of soil/structure interface.Under mixed soil slurry,strain softening and shear dilatation are observed,while shear dilatation appears under the small normal stress of the interface without slurry,and shear contraction is significant under the condition of the bentonite slurry.The thickness of the interface was determined by analyzing the disturbed height of the sample with both simple shear test and particle flow code (PFC).An elasto-plastic constitutive model incorporating strain softening and dilatancy for thin layer element of interface was formulated in the framework of generalized potential theory.The relation curves of shear stress and shear strain,as well as the relation curves of normal strain and shear strain,were fitted by a piecewise function composed by hyperbolic functions and resembling normal functions.The entire model parameters can be identified by tests.The new model is verified by comparing the measured data of indoor cut-off wall model tests with the predictions from finite element method (FEM).The FEM results indicate that the stress of wall calculated by using Goodman element is too large,and the maximum deviation between the test data and prediction is about 45%.While the prediction from the proposed model is close to the measured data,and the error is generally less than 10%.展开更多
The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and tw...The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and twodimensional particle flow models were established.The principal stress and principal strain distributions surrounding the four-arc-shaped and inverted U-shaped tunnels were investigated,respectively.Numerical results indicated that the dip angle combination of preexisting fractures directly affects the principal stress,principal strain distribution and the failure characteristics around the tunnel.The larger the absolute value of the preexisting fracture inclination angle,the higher the crushing degree of compression splitting near the hance and the larger the V-shaped failure zone.With a decrease in the absolute value of the preexisting fracture inclination angle,the compressive stress concentration of the sidewall with preexisting fractures gradually increases.The types of cracks initiated around the four-arc-shaped tunnel and the inverted U-shape tunnel are different.When the fractures are almost vertical,they have a significant influence on the stress of the sidewall force of the four-arc-shaped tunnel.When the fractures are almost horizontal,they have a significant influence on the stress of the sidewall of the inverted U-shaped tunnel.The findings provide a theoretical support for the local strengthening design of the tunnel supporting structure.展开更多
Failure of rock under impact loadings involves complex micro-fracturing and progressive damage. Strength increase and splitting failure have been observed during dynamic tests of rock materials. However, the failure m...Failure of rock under impact loadings involves complex micro-fracturing and progressive damage. Strength increase and splitting failure have been observed during dynamic tests of rock materials. However, the failure mechanism still remains unclear. In this work, based on laboratory tests, numerical simulations with the particle flow code(PFC) were carried out to reproduce the micro-fracturing process of granite specimens. Shear and tensile cracks were both recorded to investigate the failure mode of rocks under different loading conditions. At the same time, a dynamic damage model based on the Weibull distribution was established to predict the deformation and degradation behavior of specimens. It is found that micro-cracks play important roles in controlling the dynamic deformation and failure process of rock under impact loadings. The sharp increase in the number of cracks may be the reason for the strength increase of rock under high strain rates. Tensile cracks tend to be the key reason for splitting failure of specimens. Numerical simulation of crack propagation by PFC can give vivid description of the failure process. However, it is not enough for evaluation of material degradation. The dynamic damage model is able to predict the stress-strain relationship of specimens reasonably well, and can be used to explain the degradation of specimens under impact loadings at macro-scale. Crack and damage can describe material degradation at different scales and can be used together to reveal the failure mechanism of rocks.展开更多
基金Project（BK20150005） supported by the Natural Science Foundation of Jiangsu Province for Distinguished Young Scholars, China Project（2015XKZD05） supported by the Fundamental Research Funds for the Central Universities, China
文摘Permeability is a vital property of rock mass, which is highly affected by tectonic stress and human engineering activities. A comprehensive monitoring of pore pressure and flow rate distributions inside the rock mass is very important to elucidate the permeability evolution mechanisms, which is difficult to realize in laboratory, but easy to be achieved in numerical simulations. Therefore, the particle flow code （PFC）, a discrete element method, is used to simulate permeability behaviors of rock materials in this study. Owe to the limitation of the existed solid-fluid coupling algorithm in PFC, an improved flow-coupling algorithm is presented to better reflect the preferential flow in rock fractures. The comparative analysis is conducted between original and improved algorithm when simulating rock permeability evolution during triaxial compression, showing that the improved algorithm can better describe the experimental phenomenon. Furthermore, the evolution of pore pressure and flow rate distribution during the flow process are analyzed by using the improved algorithm. It is concluded that during the steady flow process in the fractured specimen, the pore pressure and flow rate both prefer transmitting through the fractures rather than rock matrix. Based on the results, fractures are divided into the following three types： I） fractures link to both the inlet and outlet, II） fractures only link to the inlet, and III） fractures only link to the outlet. The type I fracture is always the preferential propagating path for both the pore pressure and flow rate. For type II fractures, the pore pressure increases and then becomes steady. However, the flow rate increases first and begins to decrease after the flow reaches the stop end of the fracture and finally vanishes. There is no obvious pore pressure or flow rate concentration within type III fractures.
基金Project(52274148)supported by the National Natural Science Foundation of ChinaProject(2022XJLJ01)supported by the Fundamental Research Funds for the Central Universities,China。
文摘In this paper,Brazilian test was performed on disk samples of analogue materials with defined structural planes.The surface strain evolution process of the disk samples during loading was analyzed via digital image correlation.The damage evolution process was explored from a microscopic perspective by combining discrete element numerical simulation technology.The criterion of the failure mode of the disc specimen in the split state was theoretically deduced.The influence of structural surface roughness and loading inclination angle on the stress state at the center of the specimen was explored.The results showed that the failure modes of the samples could be divided into three typical modes as matrix failure,structural plane failure and combination failure.The rough structural plane improves the failure strength of the specimen by limiting its lateral deformation,and the degree of improvement weakens continuously with the increase of the inclination angle of the structural plane.As the inclination angle of the structural plane increases,the main type of microcracks in the structural plane changes from shear microcracks to tensile microcracks.This study contributes to a better understanding of macro-and meso-failure characteristics of rock masses with structural planes under a splitting state.
基金Projects(51879148,51709159,51911530214)supported by the National Natural Science Foundation of ChinaProject(2019GSF111030)supported by Shandong Provincial Key R&D Program of ChinaProject(KT201804)supported by the Project of Special Fund for Science and Technology of Water Resources Department of Guizhou Province,China。
文摘Hydraulic fracturing,as a key technology of deep energy exploitation,accelerates the rapid development of the modern petroleum industry.To study the mechanisms of hydraulic fracture propagation and rock failure mode of the vertical well hydraulic fracturing,the true triaxial hydraulic fracturing test and numerical simulation are carried out,and the influence of the principal stress difference,water injection displacement,perforation angle and natural fracture on fracture propagation is analyzed.The results show that the fracture propagation mode of limestone is mainly divided into two types:the single vertical fracture and the transverse-longitudinal crossed complex fracture.Under high displacement,the fracturing pressure is larger,and the secondary fracture is more likely to occur,while variable displacement loading is more likely to induce fracture network.Meanwhile,the amplitude of acoustic emission(AE)waveform of limestone during fracturing is between 0.01 and 0.02 mV,and the main frequency is maintained in the range of 230−300 kHz.When perforation angleθ=45°,it is easy to produce the T-type fracture that connects with the natural fracture,while X-type cracks are generated whenθ=30°.The results can be used as a reference for further study on the mechanism of limestone hydraulic fracturing.
基金Projects(51474251,51874351)supported by the National Natural Science Foundation,China。
文摘It is important to calibrate micro-parameters for applying partied flow code(PFC)to study mechanical characteristics and failure mechanism of rock materials.Uniform design method is firstly adopted to determine the microscopic parameters of parallel-bonded particle model for three-dimensional discrete element particle flow code(PFC3D).Variation ranges of microscopic of the microscopic parameters are created by analyzing the effects of microscopic parameters on macroscopic parameters(elastic modulus E,Poisson ratio v,uniaxial compressive strengthσc,and ratio of crack initial stress to uniaxial compressive strengthσci/σc)in order to obtain the actual uniform design talbe.The calculation equations of the microscopic and macroscopic parameters of rock materials can be established by the actual uniform design table and the regression analysis and thus the PFC3D microscopic parameters can be quantitatively determined.The PFC3D simulated results of the intact and pre-cracked rock specimens under uniaxial and triaxial compressions(including the macroscopic mechanical parameters,stress−strain curves and failure process)are in good agreement with experimental results,which can prove the validity of the calculation equations of microscopic and macroscopic parameters.
基金Projects(51621006,51874274)supported by the National Natural Science Foundation of ChinaProject(2018YFC0808401)supported by the National Key Research and Development Program of China
文摘In this paper,a new micro-creep model of salt rock is proposed based on a linear parallel bonded model(LPBM)using the two-dimensional particle flow code(PFC2D).The power function weakening form is assumed to describe the variation of the parallel bonded diameter(PBD)over time.By comparing with the parallel-bonded stress corrosion(PSC)model,a smaller stress fluctuation and smoother creep strain−time curves can be obtained by this power function model at the same stress level.The validity and adaptability of the model to simulate creep deformation of salt rock are verified through comparing the laboratory creep test curves and the Burgers model fitting result.The numerical results reveal that this model can be capable of capturing the creep deformation and damage behavior from the laboratory observations.
基金Project(20110094110002) supported by the Specialized Research Fund for the Doctoral Program of Higher Education of ChinaProject(200801014) supported by the Ministry of Water Resources of ChinaProject(50825901) supported by the National Natural Science Foundation of China
文摘The mechanical behaviors of the interface between coarse-grained soil and concrete were investigated by simple shear tests under condition of mixed soil slurry (bentonite mixed with cement grout).For comparison,the interfaces both without slurry and with bentonite slurry were analyzed.The experimental results show that different slurries exert much influence on the strength and deformation of soil/structure interface.Under mixed soil slurry,strain softening and shear dilatation are observed,while shear dilatation appears under the small normal stress of the interface without slurry,and shear contraction is significant under the condition of the bentonite slurry.The thickness of the interface was determined by analyzing the disturbed height of the sample with both simple shear test and particle flow code (PFC).An elasto-plastic constitutive model incorporating strain softening and dilatancy for thin layer element of interface was formulated in the framework of generalized potential theory.The relation curves of shear stress and shear strain,as well as the relation curves of normal strain and shear strain,were fitted by a piecewise function composed by hyperbolic functions and resembling normal functions.The entire model parameters can be identified by tests.The new model is verified by comparing the measured data of indoor cut-off wall model tests with the predictions from finite element method (FEM).The FEM results indicate that the stress of wall calculated by using Goodman element is too large,and the maximum deviation between the test data and prediction is about 45%.While the prediction from the proposed model is close to the measured data,and the error is generally less than 10%.
基金Project(41807241) supported by the National Natural Science Foundation of ChinaProject(2021M693544) supported by China Postdoctoral Science FoundationProject(2022JM-160) supported by the Natural Science Basic Research Program of Shaanxi,China。
文摘The complexity of a rock masses structure can lead to high uncertainties and risk during underground engineering construction.Laboratory tests on fractured rock-like materials containing a tunnel were conducted,and twodimensional particle flow models were established.The principal stress and principal strain distributions surrounding the four-arc-shaped and inverted U-shaped tunnels were investigated,respectively.Numerical results indicated that the dip angle combination of preexisting fractures directly affects the principal stress,principal strain distribution and the failure characteristics around the tunnel.The larger the absolute value of the preexisting fracture inclination angle,the higher the crushing degree of compression splitting near the hance and the larger the V-shaped failure zone.With a decrease in the absolute value of the preexisting fracture inclination angle,the compressive stress concentration of the sidewall with preexisting fractures gradually increases.The types of cracks initiated around the four-arc-shaped tunnel and the inverted U-shape tunnel are different.When the fractures are almost vertical,they have a significant influence on the stress of the sidewall force of the four-arc-shaped tunnel.When the fractures are almost horizontal,they have a significant influence on the stress of the sidewall of the inverted U-shaped tunnel.The findings provide a theoretical support for the local strengthening design of the tunnel supporting structure.
基金Projects(51274254,51322403)supported by the National Natural Science Foundation of ChinaProject(NCET-11-0528)supported by theProgram for New Century Excellent Talents in University,ChinaProject(2013SK2011)supported by Hunan Province Science andTechnology Plan,China
文摘Failure of rock under impact loadings involves complex micro-fracturing and progressive damage. Strength increase and splitting failure have been observed during dynamic tests of rock materials. However, the failure mechanism still remains unclear. In this work, based on laboratory tests, numerical simulations with the particle flow code(PFC) were carried out to reproduce the micro-fracturing process of granite specimens. Shear and tensile cracks were both recorded to investigate the failure mode of rocks under different loading conditions. At the same time, a dynamic damage model based on the Weibull distribution was established to predict the deformation and degradation behavior of specimens. It is found that micro-cracks play important roles in controlling the dynamic deformation and failure process of rock under impact loadings. The sharp increase in the number of cracks may be the reason for the strength increase of rock under high strain rates. Tensile cracks tend to be the key reason for splitting failure of specimens. Numerical simulation of crack propagation by PFC can give vivid description of the failure process. However, it is not enough for evaluation of material degradation. The dynamic damage model is able to predict the stress-strain relationship of specimens reasonably well, and can be used to explain the degradation of specimens under impact loadings at macro-scale. Crack and damage can describe material degradation at different scales and can be used together to reveal the failure mechanism of rocks.
