The fracture toughness of rocks is a critical fracturing parameter in geo-energy exploitation playing a significant role in fracture mechanics and hydraulic fracturing.The edge-notched disk bending(ENDB)specimens are ...The fracture toughness of rocks is a critical fracturing parameter in geo-energy exploitation playing a significant role in fracture mechanics and hydraulic fracturing.The edge-notched disk bending(ENDB)specimens are employed to measure the entire range of mixed-modeⅠ/Ⅲfracture toughness of Longmaxi shale.To theoretically interpret the fracture mechanisms,this research first introduces the detailed derivations of three established fracture criteria.By distinguishing the volumetric and distortional strain energy densities,an improved three-dimensional mean strain energy density(MSED)criterion is proposed.As the critical volumetric to distortional MSED ratio decreases,the transition from tensiondominated fracture to shear-dominated fracture is observed.Our results indicate that both peak load and applied energy increase significantly with the transition from pure mode I(i.e.,tension)to pure modeⅢ(i.e.,torsion or tearing)since mode-Ⅲcracking happens in a twisted manner and mode-Ⅰcracking occurs in a coplanar manner.The macroscopic fracture signatures are consistent with those of triaxial hydraulic fracturing.The average ratio of pure mode-Ⅲfracture toughness to pure mode-Ⅰfracture toughness is 0.68,indicating that the obtained mode-Ⅲfracture resistance for a tensionbased loading system is apparent rather than true.Compared to the three mainstream fracture criteria,the present fracture criterion exhibits greater competitiveness and can successfully evaluate and predict mixed-modeⅠ/Ⅲfracture toughness of distinct materials and loading methods.展开更多
Mechanical and electrical damages are introduced to study the fracture mechanics of piezoelectric ceramics in this paper. Two kinds of piezoelectric fracture criteria are established using the method of least squares ...Mechanical and electrical damages are introduced to study the fracture mechanics of piezoelectric ceramics in this paper. Two kinds of piezoelectric fracture criteria are established using the method of least squares combined with a damage analysis of the well-known piezoelectric fracture experiments of Park and Sun’s. One is based on a linear combination of the mechanical and electrical damages and the other on their nonlinear combination. When the combined damage D is up to its critical value Dc, piezoelectric fracture occurs. It is found from the qualitative comparison of the numerical results with the experimental data that the nonlinearly combined damage fracture criterion can give a better prediction of piezoelectric fracture. And it is concluded from the nonlinearly combined damage fracture criterion that a negative electric feld impedes fracture whereas the efect of a positive electric feld on fracture depends on its magnitude.展开更多
A modified MK model combined with ductile fracture criterion（DFC-MK model） is proposed to compute the forming limit diagrams（FLDs） of 5A06-O aluminum alloy sheet at different temperatures.The material constant（C...A modified MK model combined with ductile fracture criterion（DFC-MK model） is proposed to compute the forming limit diagrams（FLDs） of 5A06-O aluminum alloy sheet at different temperatures.The material constant（C） of ductile fracture criterion and initial thickness imperfection parameter（f0） at various temperatures are determined by using a new computing method based on wide sheet bending test.The FLDs at 20 and 200 °C are calculated through the DFC-MK model.The DFC-MK model,which includes the influence of through-thickness normal stress,is written into the subroutine VUMAT embedded in Abaqus/ Explicit.The cylindrical cup hydroforming tests are carried out to verify the model.The results show that compared with experimental observations,the predicted FLDs based on DFC-MK model are more accurate than the conventional MK model;the errors between the simulations and experiments in warm hydroforming are 8.23% at 20 °C and 9.24% at 200 °C,which verify the effectiveness of the proposed model.展开更多
Knowledge of the strength and deformability of fractured rocks is important for design, construction and stability evaluation of slopes, foundations and underground excavations in civil and mining engineering. However...Knowledge of the strength and deformability of fractured rocks is important for design, construction and stability evaluation of slopes, foundations and underground excavations in civil and mining engineering. However, laboratory tests of intact rock samples cannot provide information about the strength and deformation behaviors of fractured rock masses that include many fractures of varying sizes, orientations and locations. On the other hand, large-scale in situ tests of fractured rock masses are economically costly and often not practical in reality at present. Therefore, numerical modeling becomes necessary. Numerical predicting using discrete element methods(DEM) is a suitable approach for such modeling because of their advantages of explicit representations of both fractures system geometry and their constitutive behaviors of fractures, besides that of intact rock matrix. In this study, to generically determine the compressive strength of fractured rock masses, a series of numerical experiments were performed on two-dimensional discrete fracture network models based on the realistic geometrical and mechanical data of fracture systems from feld mapping. We used the UDEC code and a numerical servo-controlled program for controlling the progressive compressive loading process to avoid sudden violent failure of the models. The two loading conditions applied are similar to the standard laboratory testing for intact rock samples in order to check possible differences caused by such loading conditions. Numerical results show that the strength of fractured rocks increases with the increasing confning pressure, and that deformation behavior of fractured rocks follows elasto-plastic model with a trend of strain hardening. The stresses and strains obtained from these numerical experiments were used to ft the well-known Mohr-Coulomb(MC) and Hoek-Brown(H-B) failure criteria, represented by equivalent material properties defning these two criteria. The results show that both criteria can provide fair estimates of the compressive strengths for all tested numerical models. Parameters of the elastic deformability of fractured models during elastic deformation stages were also evaluated, and represented as equivalent Young’s modulus and Poisson’s ratio as functions of lateral confning pressure. It is the frst time that such systematic numerical predicting for strength of fractured rocks was performed considering different loading conditions, with important fndings for different behaviors of fractured rock masses, compared with testing intact rock samples under similar loading conditions.展开更多
With complex fractured-vuggy heterogeneous structures, water has to be injected to facilitate oil pro- duction. However, the effect of different water injection modes on oil recovery varies. The limitation of existing...With complex fractured-vuggy heterogeneous structures, water has to be injected to facilitate oil pro- duction. However, the effect of different water injection modes on oil recovery varies. The limitation of existing numerical simulation methods in representing fractured- vuggy carbonate reservoirs makes numerical simulation difficult to characterize the fluid flow in these reservoirs. In this paper, based on a geological example unit in the Tahe Oilfield, a three-dimensional physical model was designed and constructed to simulate fluid flow in a fractured-vuggy reservoir according to similarity criteria. The model was validated by simulating a bottom water drive reservoir, and then subsequent water injection modes were optimized. These were continuous （constant rate）, intermittent, and pulsed injection of water. Experimental results reveal that due to the unbalanced formation pressure caused by pulsed water injection, the swept volume was expanded and consequently the highest oil recovery increment was achieved. Similar to continuous water injection, intermit- tent injection was influenced by factors including the connectivity of the fractured-vuggy reservoir, well depth, and the injection-production relationship, which led to a relative low oil recovery. This study may provide a constructive guide to field production and for the devel- opment of the commercial numerical models specialized for fractured-vuggy carbonate reservoirs.展开更多
Based on the traditional theory of fracture, time dependence of rock fracture is taken into account and the creep fracture criteria are given. The effect of the microcrack zone on fracture toughness is analyzed. Fina...Based on the traditional theory of fracture, time dependence of rock fracture is taken into account and the creep fracture criteria are given. The effect of the microcrack zone on fracture toughness is analyzed. Finally, a model for time prediction of rheological fracture is derived and its feasibility is proven.展开更多
In recent years,the issue of aircraft icing has gained widespread recognition.The breaking and detachment of dynamic ice can pose a threat to flight safety.However,the shedding and fracture mechanisms of dynamic ice a...In recent years,the issue of aircraft icing has gained widespread recognition.The breaking and detachment of dynamic ice can pose a threat to flight safety.However,the shedding and fracture mechanisms of dynamic ice are unclear and cannot meet the engineering needs of ice-shedding hazard assessment.Therefore,studying the fracture toughness of ice bodies has extremely important practical significance.To address this issue,this article uses a centrally cracked Brazilian disk(CCBD)specimen to measure the pure modeⅠtoughness and pure modeⅡfracture toughness of freshwater ice at different loading rates.The mixed-mode(Ⅰ–Ⅱ)fracture characteristics of ice are discussed,and the experimental results are compared and analyzed with the theoretical values of the generalized maximum tangential stress(GMTS)criterion considering the influence of T-stress.The results indicated that as the loading rate increases,the pure modeⅠtoughness and pure modeⅡfracture toughness of freshwater ice decrease,and the fracture toughness of freshwater ice is more sensitive to the loading rate.Ⅰn terms of fracture criteria,the theoretical value of the ratio of pure modeⅡfracture toughness to pure modeⅠfracture toughness based on the GMTS criterion is in good agreement with the experimental value,while the theoretical value based on the maximum tangential stress(MTS)criterion deviates significantly from the experimental value,indicating that the GMTS criterion considering the influence of T-stress can better predict the experimental results.展开更多
To identify the type of main flow channels of complex porous media in oil and gas reservoirs,the"main flow channel index"is defined as the ratio of comprehensive permeability obtained from well test to matri...To identify the type of main flow channels of complex porous media in oil and gas reservoirs,the"main flow channel index"is defined as the ratio of comprehensive permeability obtained from well test to matrix permeability obtained from core analysis or well logging.Meanwhile,a mathematical model is established based on equivalent flow assumption,the classification method for main flow channels is put forward,and quantitative characterization of main flow channels is realized.The method has been verified by analysis of typical gas reservoirs.The study results show that the"main flow channel index"can quantitatively classify types of flow channels.If the index is less than 3,the matrix pore is the main flow channel;if the index is between 3 and 20,the fracture is the main flow channel and the matrix pore acts as the supplement one;if the index is more than 20,the fracture is the only seepage channel.The dynamic analysis of typical gas reservoirs shows that the"main flow channel index"can be used to identify the type of flow channel in complex porous media,guiding the classified development of gas reservoirs,and avoiding development risk.展开更多
为提高厚煤层采空区定向钻孔的瓦斯抽采效率,针对山西某高瓦斯矿井采煤工作面,采用理论分析和FLUENT数值模拟相结合的方法研究采动裂隙分域演化特性,提出覆岩裂隙场分域准则,确定定向钻孔布置区域与核心抽采布置范围,并在采空区现场开...为提高厚煤层采空区定向钻孔的瓦斯抽采效率,针对山西某高瓦斯矿井采煤工作面,采用理论分析和FLUENT数值模拟相结合的方法研究采动裂隙分域演化特性,提出覆岩裂隙场分域准则,确定定向钻孔布置区域与核心抽采布置范围,并在采空区现场开展定向钻孔分域抽采瓦斯试验。结果表明:破断裂隙密集区内,岩层断裂穿层裂隙发育较明显、瓦斯聚集显著,且钻孔稳定性高,是布置定向钻孔的最佳区域;并将与回风巷中心线水平距离3~13 m,与煤层顶板垂直距离10~18 m的区域设定为核心抽采区域。定向钻孔分域抽采试验中,单孔抽采瓦斯体积分数平均提升22.355%,单孔瓦斯抽采纯量平均提升1.295 m 3/min,该结论验证了厚煤层采空区定向钻孔分域抽采方法的实用性与合理性。展开更多
An FAE (Fuel Air Explosives) device is used to develop a numerical and theoretical analysis of a thin cylindrical shell with inner explosive loading. The dynamic fracture process is simulated numerically in the DYNA...An FAE (Fuel Air Explosives) device is used to develop a numerical and theoretical analysis of a thin cylindrical shell with inner explosive loading. The dynamic fracture process is simulated numerically in the DYNA3D program using the finite element method. The material’s dynamic properties are described by a strain hardening viscoplastic constitution. A damage variable is introduced in the determination of the dynamic fracture criterion. Final rupture of structure is decided by a rupture strain criterion which is deduced in terms of a critical damage variable. The numerical results have been compared with theoretical solutions.展开更多
基金supported by National Natural Science Foundation of China(Grant Nos.52364004,52264006,and 52164001).
文摘The fracture toughness of rocks is a critical fracturing parameter in geo-energy exploitation playing a significant role in fracture mechanics and hydraulic fracturing.The edge-notched disk bending(ENDB)specimens are employed to measure the entire range of mixed-modeⅠ/Ⅲfracture toughness of Longmaxi shale.To theoretically interpret the fracture mechanisms,this research first introduces the detailed derivations of three established fracture criteria.By distinguishing the volumetric and distortional strain energy densities,an improved three-dimensional mean strain energy density(MSED)criterion is proposed.As the critical volumetric to distortional MSED ratio decreases,the transition from tensiondominated fracture to shear-dominated fracture is observed.Our results indicate that both peak load and applied energy increase significantly with the transition from pure mode I(i.e.,tension)to pure modeⅢ(i.e.,torsion or tearing)since mode-Ⅲcracking happens in a twisted manner and mode-Ⅰcracking occurs in a coplanar manner.The macroscopic fracture signatures are consistent with those of triaxial hydraulic fracturing.The average ratio of pure mode-Ⅲfracture toughness to pure mode-Ⅰfracture toughness is 0.68,indicating that the obtained mode-Ⅲfracture resistance for a tensionbased loading system is apparent rather than true.Compared to the three mainstream fracture criteria,the present fracture criterion exhibits greater competitiveness and can successfully evaluate and predict mixed-modeⅠ/Ⅲfracture toughness of distinct materials and loading methods.
基金Project supported by the National Natural Science Foundation of China(No.10172036)and by the Scientific ResearchFoundation for Returned Overseas Chinese Scholars,State Education Ministry.
文摘Mechanical and electrical damages are introduced to study the fracture mechanics of piezoelectric ceramics in this paper. Two kinds of piezoelectric fracture criteria are established using the method of least squares combined with a damage analysis of the well-known piezoelectric fracture experiments of Park and Sun’s. One is based on a linear combination of the mechanical and electrical damages and the other on their nonlinear combination. When the combined damage D is up to its critical value Dc, piezoelectric fracture occurs. It is found from the qualitative comparison of the numerical results with the experimental data that the nonlinearly combined damage fracture criterion can give a better prediction of piezoelectric fracture. And it is concluded from the nonlinearly combined damage fracture criterion that a negative electric feld impedes fracture whereas the efect of a positive electric feld on fracture depends on its magnitude.
基金Project(51175024)supported by the National Natural Science Foundation of China
文摘A modified MK model combined with ductile fracture criterion（DFC-MK model） is proposed to compute the forming limit diagrams（FLDs） of 5A06-O aluminum alloy sheet at different temperatures.The material constant（C） of ductile fracture criterion and initial thickness imperfection parameter（f0） at various temperatures are determined by using a new computing method based on wide sheet bending test.The FLDs at 20 and 200 °C are calculated through the DFC-MK model.The DFC-MK model,which includes the influence of through-thickness normal stress,is written into the subroutine VUMAT embedded in Abaqus/ Explicit.The cylindrical cup hydroforming tests are carried out to verify the model.The results show that compared with experimental observations,the predicted FLDs based on DFC-MK model are more accurate than the conventional MK model;the errors between the simulations and experiments in warm hydroforming are 8.23% at 20 °C and 9.24% at 200 °C,which verify the effectiveness of the proposed model.
文摘Knowledge of the strength and deformability of fractured rocks is important for design, construction and stability evaluation of slopes, foundations and underground excavations in civil and mining engineering. However, laboratory tests of intact rock samples cannot provide information about the strength and deformation behaviors of fractured rock masses that include many fractures of varying sizes, orientations and locations. On the other hand, large-scale in situ tests of fractured rock masses are economically costly and often not practical in reality at present. Therefore, numerical modeling becomes necessary. Numerical predicting using discrete element methods(DEM) is a suitable approach for such modeling because of their advantages of explicit representations of both fractures system geometry and their constitutive behaviors of fractures, besides that of intact rock matrix. In this study, to generically determine the compressive strength of fractured rock masses, a series of numerical experiments were performed on two-dimensional discrete fracture network models based on the realistic geometrical and mechanical data of fracture systems from feld mapping. We used the UDEC code and a numerical servo-controlled program for controlling the progressive compressive loading process to avoid sudden violent failure of the models. The two loading conditions applied are similar to the standard laboratory testing for intact rock samples in order to check possible differences caused by such loading conditions. Numerical results show that the strength of fractured rocks increases with the increasing confning pressure, and that deformation behavior of fractured rocks follows elasto-plastic model with a trend of strain hardening. The stresses and strains obtained from these numerical experiments were used to ft the well-known Mohr-Coulomb(MC) and Hoek-Brown(H-B) failure criteria, represented by equivalent material properties defning these two criteria. The results show that both criteria can provide fair estimates of the compressive strengths for all tested numerical models. Parameters of the elastic deformability of fractured models during elastic deformation stages were also evaluated, and represented as equivalent Young’s modulus and Poisson’s ratio as functions of lateral confning pressure. It is the frst time that such systematic numerical predicting for strength of fractured rocks was performed considering different loading conditions, with important fndings for different behaviors of fractured rock masses, compared with testing intact rock samples under similar loading conditions.
基金supported by China National Science and Technology Major Project(2011ZX05009-004,2011ZX05014-003)National Key Basic Research and Development Program(973 Program),China(2011CB201006)Science Foundation of China University of Petroleum,Beijing(2462014YJRC053)
文摘With complex fractured-vuggy heterogeneous structures, water has to be injected to facilitate oil pro- duction. However, the effect of different water injection modes on oil recovery varies. The limitation of existing numerical simulation methods in representing fractured- vuggy carbonate reservoirs makes numerical simulation difficult to characterize the fluid flow in these reservoirs. In this paper, based on a geological example unit in the Tahe Oilfield, a three-dimensional physical model was designed and constructed to simulate fluid flow in a fractured-vuggy reservoir according to similarity criteria. The model was validated by simulating a bottom water drive reservoir, and then subsequent water injection modes were optimized. These were continuous （constant rate）, intermittent, and pulsed injection of water. Experimental results reveal that due to the unbalanced formation pressure caused by pulsed water injection, the swept volume was expanded and consequently the highest oil recovery increment was achieved. Similar to continuous water injection, intermit- tent injection was influenced by factors including the connectivity of the fractured-vuggy reservoir, well depth, and the injection-production relationship, which led to a relative low oil recovery. This study may provide a constructive guide to field production and for the devel- opment of the commercial numerical models specialized for fractured-vuggy carbonate reservoirs.
文摘Based on the traditional theory of fracture, time dependence of rock fracture is taken into account and the creep fracture criteria are given. The effect of the microcrack zone on fracture toughness is analyzed. Finally, a model for time prediction of rheological fracture is derived and its feasibility is proven.
基金supported by the National Natural Science Foundation of China(Nos.12132019 and 11872042)the Open Fund for Key Laboratory of Deep Underground Science and Engineering of Ministry of Education(No.DESEYU202301)+1 种基金the 2023 Open Project of Failure Mechanics and Engineering Disaster Prevention,Key Lab of Sichuan Province(No.FMEDP202306)the Natural Science Foundation of Sichuan Province(No.2023NSFSC0043).
文摘In recent years,the issue of aircraft icing has gained widespread recognition.The breaking and detachment of dynamic ice can pose a threat to flight safety.However,the shedding and fracture mechanisms of dynamic ice are unclear and cannot meet the engineering needs of ice-shedding hazard assessment.Therefore,studying the fracture toughness of ice bodies has extremely important practical significance.To address this issue,this article uses a centrally cracked Brazilian disk(CCBD)specimen to measure the pure modeⅠtoughness and pure modeⅡfracture toughness of freshwater ice at different loading rates.The mixed-mode(Ⅰ–Ⅱ)fracture characteristics of ice are discussed,and the experimental results are compared and analyzed with the theoretical values of the generalized maximum tangential stress(GMTS)criterion considering the influence of T-stress.The results indicated that as the loading rate increases,the pure modeⅠtoughness and pure modeⅡfracture toughness of freshwater ice decrease,and the fracture toughness of freshwater ice is more sensitive to the loading rate.Ⅰn terms of fracture criteria,the theoretical value of the ratio of pure modeⅡfracture toughness to pure modeⅠfracture toughness based on the GMTS criterion is in good agreement with the experimental value,while the theoretical value based on the maximum tangential stress(MTS)criterion deviates significantly from the experimental value,indicating that the GMTS criterion considering the influence of T-stress can better predict the experimental results.
文摘To identify the type of main flow channels of complex porous media in oil and gas reservoirs,the"main flow channel index"is defined as the ratio of comprehensive permeability obtained from well test to matrix permeability obtained from core analysis or well logging.Meanwhile,a mathematical model is established based on equivalent flow assumption,the classification method for main flow channels is put forward,and quantitative characterization of main flow channels is realized.The method has been verified by analysis of typical gas reservoirs.The study results show that the"main flow channel index"can quantitatively classify types of flow channels.If the index is less than 3,the matrix pore is the main flow channel;if the index is between 3 and 20,the fracture is the main flow channel and the matrix pore acts as the supplement one;if the index is more than 20,the fracture is the only seepage channel.The dynamic analysis of typical gas reservoirs shows that the"main flow channel index"can be used to identify the type of flow channel in complex porous media,guiding the classified development of gas reservoirs,and avoiding development risk.
文摘为提高厚煤层采空区定向钻孔的瓦斯抽采效率,针对山西某高瓦斯矿井采煤工作面,采用理论分析和FLUENT数值模拟相结合的方法研究采动裂隙分域演化特性,提出覆岩裂隙场分域准则,确定定向钻孔布置区域与核心抽采布置范围,并在采空区现场开展定向钻孔分域抽采瓦斯试验。结果表明:破断裂隙密集区内,岩层断裂穿层裂隙发育较明显、瓦斯聚集显著,且钻孔稳定性高,是布置定向钻孔的最佳区域;并将与回风巷中心线水平距离3~13 m,与煤层顶板垂直距离10~18 m的区域设定为核心抽采区域。定向钻孔分域抽采试验中,单孔抽采瓦斯体积分数平均提升22.355%,单孔瓦斯抽采纯量平均提升1.295 m 3/min,该结论验证了厚煤层采空区定向钻孔分域抽采方法的实用性与合理性。
文摘An FAE (Fuel Air Explosives) device is used to develop a numerical and theoretical analysis of a thin cylindrical shell with inner explosive loading. The dynamic fracture process is simulated numerically in the DYNA3D program using the finite element method. The material’s dynamic properties are described by a strain hardening viscoplastic constitution. A damage variable is introduced in the determination of the dynamic fracture criterion. Final rupture of structure is decided by a rupture strain criterion which is deduced in terms of a critical damage variable. The numerical results have been compared with theoretical solutions.
