The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, where...The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples shaw that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks ate in good agreement with other analytical solutions.展开更多
Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both ...Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.展开更多
The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped i...The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.展开更多
Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numeri...Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numerical results is made by other researchers' experimental results, and significant parameters affecting fatigue crack growth are studied.展开更多
In existing studies, most slope stability analyses concentrate on conditions with constant temperature, assuming the slope is intact, and employ the Mohr-Coulomb (M-C) failure criterion for saturated soil to character...In existing studies, most slope stability analyses concentrate on conditions with constant temperature, assuming the slope is intact, and employ the Mohr-Coulomb (M-C) failure criterion for saturated soil to characterize the strength of the backfill. However, the actual working temperature of slopes varies, and natural phenomena such as rainfall and groundwater infiltration commonly result in unsaturated soil conditions, with cracks typically present in cohesive slopes. This study introduces a novel approach for assessing the stability of unsaturated soil stepped slopes under varying temperatures, incorporating the effects of open and vertical cracks. Utilizing the kinematic approach and gravity increase method, we developed a three-dimensional (3D) rotational wedge failure mechanism to simulate slope collapse, enhancing the traditional two-dimensional analyses. We integrated temperature-dependent functions and nonlinear shear strength equations to evaluate the impact of temperature on four typical unsaturated soil types. A particle swarm optimization algorithm was employed to calculate the safety factor, ensuring our method’s accuracy by comparing it with existing studies. The results indicate that considering 3D effects yields a higher safety factor, while cracks reduce slope stability. Each unsaturated soil exhibits a distinctive temperature response curve, highlighting the importance of understanding soil types in the design phase.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
Accurate characterization of three-dimensional burning crack propagation remains pivotal yet challenging for energetic material safety,as conventional diagnostics and models inadequately resolve coupled crack-pressure...Accurate characterization of three-dimensional burning crack propagation remains pivotal yet challenging for energetic material safety,as conventional diagnostics and models inadequately resolve coupled crack-pressure dynamics in confined explosives.This study combines a novel spherical confinement system(with/without sapphire windows)with synchronized high-speed imaging and 3D reconstruction to overcome optical limitations in opaque explosives.Experimental analysis of centrally ignited HMX-based PBX-1 reveals:(1)burning cracks propagate radially with equatorial acceleration and polar deceleration,(2)systematic formation of 3–4 dominant crack branches across geometries,and(3)pressure evolution exhibiting gradual accumulation(subsurface cracking)followed by exponential growth(surface burn-through),with decay governed by cavity expansion.Building on Hill's framework,we develop a model incorporating cavity volume and fracture toughness criteria,validated against PBX explosive(95%HMX-based)experiments.The model demonstrates improved prediction of pressure trends compared to prior approaches,particularly in resolving laminar-phase accumulation and crackinduced surge transitions.Results establish structural cavity volume as a critical modulator of measured pressure and reveal direction-dependent crack kinematics as fundamental features of constrained combustion.This work provides experimentally validated insights into mechanisms of reaction pressure development and burning cracks pathways during constrained PBX explosive combustion.展开更多
The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional （3D） growth and distribution of hy...The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional （3D） growth and distribution of hydrofracturing cracks within heterogeneous rocks is of particular use to the design and implementation of hydrofracturing stimulation of unconventional reservoirs. However, because of the difficulties involved in visually representing and quantitatively characterizing a 3D hydrofracturing crack-network, this issue remains a challenge. In this paper, a novel method is proposed for physically visualizing and quantitatively characterizing the 3D hydrofracturing crack-network distributed through a heterogeneous structure based on a natural glutenite sample. This method incorporates X-ray microfocus computed tomography （μCT）, 3D printing models and hydrofracturing triaxial tests to represent visually the heterogeneous structure, and the 3D crack growth and distribution within a transparent rock model during hydrofracturing. The coupled effects of material heterogeneity and confining geostress on the 3D crack initiation and propagation were analyzed. The results indicate that the breakdown pressure of a heterogeneous rock model is significantly affected by material heterogeneity and confining geostress. The measured breakdown pressures of heterogeneous models are apparently different from those predicted by traditional theories. This study helps to elucidate the quantitative visualization and characterization of the mechanism and influencing factors that determine the hydrofracturing crack initiation and propagation in heterogeneous reservoir rocks.展开更多
In this study, it was assumed that three-dimensional penny-shaped cracks existed in deep rock masses. A new non-Euclidean model was established, in which the effects of penny- shaped cracks and axial in-situ stress on...In this study, it was assumed that three-dimensional penny-shaped cracks existed in deep rock masses. A new non-Euclidean model was established, in which the effects of penny- shaped cracks and axial in-situ stress on zonal disintegration of deep rock masses were taken into account. Based on the non-Euclidean model, the stress intensity factors at tips of the penny- shaped cracks were determined. The strain energy density factor was applied to investigate the occurrence of fractured zones. It was observed from the numerical results that the magnitude and location of fractured zones were sensitive to micro- and macro-mechanical parameters, as well as the value of in-situ stress. The numerical results were in good agreement with the experimental data.展开更多
By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. I...By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
Initiation and propagation of cracks in rotating disks may cause catastrophic failures. Therefore, determination of fracture parameters under different working con- ditions is an essential issue. In this paper, a comp...Initiation and propagation of cracks in rotating disks may cause catastrophic failures. Therefore, determination of fracture parameters under different working con- ditions is an essential issue. In this paper, a comprehensive study of stress intensity factors （SIFs） in rotating disks containing three-dimensional （3D） semi-elliptical cracks subjected to different working conditions is carried out. The effects of mechanical prop- erties, rotational velocity, and orientation of cracks on SIFs in rotating disks under cen- trifugal loading are investigated. Also, the effects of using composite patches to reduce SIFs in rotating disks are studied. The effects of patching design variables such as mechanical properties, thickness, and ply angle are investigated separately. The modeling and analytical procedure are verified in comparison with previously reported results in the literature.展开更多
Three-dimensional edge cracks are analyzed using the Self-SimilarCrack Expansion (SSCE)method with a boundary integral equationtechnique. The boundary integral equations for surface cracks in ahalf space are presented...Three-dimensional edge cracks are analyzed using the Self-SimilarCrack Expansion (SSCE)method with a boundary integral equationtechnique. The boundary integral equations for surface cracks in ahalf space are presented based on a half space Green'sfunction(Mindlin, 1936). By using the SSCe method, the stressintensity factors are determined by crack-opening displacement overthe crack surface. In discrete boundary integral equations, theregular and singular integrals on the crack sur- face elements areevaluated by an analytical method, and the closed form expressions ofthe integrals are given for subsurface cracks and edge cracks.展开更多
By using the Somigliana representation and the concepts offinite-part integrals, a set of hypersingular integral equations ofthe interaction between two parallel planar cracks in a three-dimen-sional finite body subje...By using the Somigliana representation and the concepts offinite-part integrals, a set of hypersingular integral equations ofthe interaction between two parallel planar cracks in a three-dimen-sional finite body subjected to arbitrary loads is derived, and thenits numerical method is proposed by the finite-part integral methodcombined with the boundary element method. According to the analytictheory of hypersingular integral equations, the square root models ofdisplacement discontinuities in the elements near the crack front areapplied, and thus the computational precision is raised. Based onthis, the stress intensity factors can be directly calculated.Finally, the stress intensity factors of sever- al typicalinteraction problems are calculated.展开更多
The relation between the normal displacement on the surface of a dynamical elliptical crack and the normal stress over the crack surface was studied. The three dimensional elastodynamic equations and Fourier Laplace...The relation between the normal displacement on the surface of a dynamical elliptical crack and the normal stress over the crack surface was studied. The three dimensional elastodynamic equations and Fourier Laplace transforms are used. Based on the influence function and the inversion of integral transforms, one can find that if the distribution of normal displacement on the surface of a dynamic elliptical crack is a polynomial of degree n in x 1 and x 2 , then the normal pressure acting over the ellipse is also a polynomial P n(x 1,x 2) of the same degree in x 1 and x 2 .展开更多
The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the rat...The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.展开更多
A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the ...A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.展开更多
Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant s...Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant sources, including force dipole, electric dipole, moment, force dilatation and rotation. Such force and charge sources may model defects like vacancies, foreign particles and dislocations. The locations and orientations of the stress and charge sources with respect to the crack are arbitrary.展开更多
Large solidification ranges and coarse columnar grains in the additively manufacturing of Al-Mg-Si alloys are normally involved in hot cracks during solidification.In this work,we develop novel crack-free Al-Mg_(2) Si...Large solidification ranges and coarse columnar grains in the additively manufacturing of Al-Mg-Si alloys are normally involved in hot cracks during solidification.In this work,we develop novel crack-free Al-Mg_(2) Si alloys fabricated by laser powder-bed fusion(L-PBF).The results indicate that the eutectic Mg_(2) Si phase possesses a strong ability to reduce crack susceptibility.It can enhance the grain growth restriction factor in the initial stage of solidification and promote eutectic filling in the terminal stage of solidifica-tion.The crack-free L-PBFed Al-x Mg_(2) Si alloys(x=6 wt.%,9 wt.%,and 12 wt.%)exhibit the combination of low crack susceptibility index(CSI),superior ability for liquid filling,and grain refinement.Particularly,the L-PBFed Al-9Mg_(2) Si alloy shows improved mechanical properties(e.g.yield strength of 397 MPa and elongation of 7.3%).However,the cracks are more likely to occur in the region near the columnar grain boundaries of the L-PBFed Al-3Mg_(2) Si alloy with a large solidification range and low eutectic content for liquid filling.Correspondingly,the L-PBFed Al-3Mg_(2) Si alloy shows poor bearing capacity of mechanical properties.The precise tuning of Mg_(2) Si eutectic content can offer an innovative strategy for eliminating cracks in additively manufactured Al-Mg-Si alloy.展开更多
The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterio...The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.展开更多
基金the National Institute of Standards and Technologythe Army Office of Research
文摘The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples shaw that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks ate in good agreement with other analytical solutions.
基金supported by the National Natural Science Foundation of China (Grant Nos. 52104125, U1765204 and 51739008)
文摘Similar to hydraulic fracturing(HF), the coalescence and fracture of cracks are induced within a rock under the action of an ultrasonic field, known as ultrasonic fracturing(UF). Investigating UF is important in both hard rock drilling and oil and gas recovery. A three-dimensional internal laser-engraved crack(3D-ILC) method was introduced to prefabricate two parallel internal cracks within the samples without any damage to the surface. The samples were subjected to UF. The mechanism of UF was elucidated by analyzing the characteristics of fracture surfaces. The crack propagation path under different ultrasonic parameters was obtained by numerical simulation based on the Paris fatigue model and compared to the experimental results of UF. The results show that the 3D-ILC method is a powerful tool for UF research.Under the action of an ultrasonic field, the fracture surface shows the characteristics of beach marks and contains powder locally, indicating that the UF mechanism includes high-cycle fatigue fracture, shear and friction, and temperature load. The two internal cracks become close under UF. The numerical result obtained by the Paris fatigue model also shows the attraction of the two cracks, consistent with the test results. The 3D-ILC method provides a new tool for the experimental study of UF. Compared to the conventional numerical methods based on the analysis of stress-strain and plastic zone, numerical simulation can be a good alternative method to obtain the crack path under UF.
文摘The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is devel- oped in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.
文摘Due to high cost of full-scale experimental setup, this study presents a numerical model on fatigue behaviours of offshore pipeline with multiple coplanar cracks under cyclic tensile loadings. The validation on numerical results is made by other researchers' experimental results, and significant parameters affecting fatigue crack growth are studied.
基金Project(51378510) supported by the National Natural Science Foundation of China。
文摘In existing studies, most slope stability analyses concentrate on conditions with constant temperature, assuming the slope is intact, and employ the Mohr-Coulomb (M-C) failure criterion for saturated soil to characterize the strength of the backfill. However, the actual working temperature of slopes varies, and natural phenomena such as rainfall and groundwater infiltration commonly result in unsaturated soil conditions, with cracks typically present in cohesive slopes. This study introduces a novel approach for assessing the stability of unsaturated soil stepped slopes under varying temperatures, incorporating the effects of open and vertical cracks. Utilizing the kinematic approach and gravity increase method, we developed a three-dimensional (3D) rotational wedge failure mechanism to simulate slope collapse, enhancing the traditional two-dimensional analyses. We integrated temperature-dependent functions and nonlinear shear strength equations to evaluate the impact of temperature on four typical unsaturated soil types. A particle swarm optimization algorithm was employed to calculate the safety factor, ensuring our method’s accuracy by comparing it with existing studies. The results indicate that considering 3D effects yields a higher safety factor, while cracks reduce slope stability. Each unsaturated soil exhibits a distinctive temperature response curve, highlighting the importance of understanding soil types in the design phase.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
基金supported by the National Natural Science Foundation of China(Grant No.12402445)the National Defense Foundation Stabilization Support Program(Grant No.JCKYS2024212108)the National Key Laboratory of Shock Wave Physics and Detonation Physics Foundation(Grant No.2024CXPTGFJJ06404)。
文摘Accurate characterization of three-dimensional burning crack propagation remains pivotal yet challenging for energetic material safety,as conventional diagnostics and models inadequately resolve coupled crack-pressure dynamics in confined explosives.This study combines a novel spherical confinement system(with/without sapphire windows)with synchronized high-speed imaging and 3D reconstruction to overcome optical limitations in opaque explosives.Experimental analysis of centrally ignited HMX-based PBX-1 reveals:(1)burning cracks propagate radially with equatorial acceleration and polar deceleration,(2)systematic formation of 3–4 dominant crack branches across geometries,and(3)pressure evolution exhibiting gradual accumulation(subsurface cracking)followed by exponential growth(surface burn-through),with decay governed by cavity expansion.Building on Hill's framework,we develop a model incorporating cavity volume and fracture toughness criteria,validated against PBX explosive(95%HMX-based)experiments.The model demonstrates improved prediction of pressure trends compared to prior approaches,particularly in resolving laminar-phase accumulation and crackinduced surge transitions.Results establish structural cavity volume as a critical modulator of measured pressure and reveal direction-dependent crack kinematics as fundamental features of constrained combustion.This work provides experimentally validated insights into mechanisms of reaction pressure development and burning cracks pathways during constrained PBX explosive combustion.
基金We gratefully acknowledge the financial support of the National Natural Science Foundation of China （Grants 51374213 and 51674251）, National Natural Science Fund for Distinguished Young Scholars of China （Grant 51125017）, Science Fund for Creative Research Groups of the National Natural Science Foundation of China （Grant 51421003）, Fund for Innovative Research and Development Group Program of Jiangsu Province （Grant 2014-27）, and the Priority Academic Program Development of Jiangsu Higher Education Institutions （Grant PAPD 2014）.
文摘The heterogeneity of unconventional reservoir rock tremendously affects its hydrofracturing behavior. A visual representation and accurate characterization of the three-dimensional （3D） growth and distribution of hydrofracturing cracks within heterogeneous rocks is of particular use to the design and implementation of hydrofracturing stimulation of unconventional reservoirs. However, because of the difficulties involved in visually representing and quantitatively characterizing a 3D hydrofracturing crack-network, this issue remains a challenge. In this paper, a novel method is proposed for physically visualizing and quantitatively characterizing the 3D hydrofracturing crack-network distributed through a heterogeneous structure based on a natural glutenite sample. This method incorporates X-ray microfocus computed tomography （μCT）, 3D printing models and hydrofracturing triaxial tests to represent visually the heterogeneous structure, and the 3D crack growth and distribution within a transparent rock model during hydrofracturing. The coupled effects of material heterogeneity and confining geostress on the 3D crack initiation and propagation were analyzed. The results indicate that the breakdown pressure of a heterogeneous rock model is significantly affected by material heterogeneity and confining geostress. The measured breakdown pressures of heterogeneous models are apparently different from those predicted by traditional theories. This study helps to elucidate the quantitative visualization and characterization of the mechanism and influencing factors that determine the hydrofracturing crack initiation and propagation in heterogeneous reservoir rocks.
基金supported by the 973 Project(No.2014CB046903)the National Natural Science Foundation of China(Nos.51325903 and 51279218)the Natural Science Foundation Project of CQ CSTC(Nos.CSTC2013KJRC-1JRCCJ30001 and CSTC2013JCYJYS0005)
文摘In this study, it was assumed that three-dimensional penny-shaped cracks existed in deep rock masses. A new non-Euclidean model was established, in which the effects of penny- shaped cracks and axial in-situ stress on zonal disintegration of deep rock masses were taken into account. Based on the non-Euclidean model, the stress intensity factors at tips of the penny- shaped cracks were determined. The strain energy density factor was applied to investigate the occurrence of fractured zones. It was observed from the numerical results that the magnitude and location of fractured zones were sensitive to micro- and macro-mechanical parameters, as well as the value of in-situ stress. The numerical results were in good agreement with the experimental data.
基金supported by the National Natural Science Foundation of China (No. 10872213)
文摘By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
文摘Initiation and propagation of cracks in rotating disks may cause catastrophic failures. Therefore, determination of fracture parameters under different working con- ditions is an essential issue. In this paper, a comprehensive study of stress intensity factors （SIFs） in rotating disks containing three-dimensional （3D） semi-elliptical cracks subjected to different working conditions is carried out. The effects of mechanical prop- erties, rotational velocity, and orientation of cracks on SIFs in rotating disks under cen- trifugal loading are investigated. Also, the effects of using composite patches to reduce SIFs in rotating disks are studied. The effects of patching design variables such as mechanical properties, thickness, and ply angle are investigated separately. The modeling and analytical procedure are verified in comparison with previously reported results in the literature.
文摘Three-dimensional edge cracks are analyzed using the Self-SimilarCrack Expansion (SSCE)method with a boundary integral equationtechnique. The boundary integral equations for surface cracks in ahalf space are presented based on a half space Green'sfunction(Mindlin, 1936). By using the SSCe method, the stressintensity factors are determined by crack-opening displacement overthe crack surface. In discrete boundary integral equations, theregular and singular integrals on the crack sur- face elements areevaluated by an analytical method, and the closed form expressions ofthe integrals are given for subsurface cracks and edge cracks.
文摘By using the Somigliana representation and the concepts offinite-part integrals, a set of hypersingular integral equations ofthe interaction between two parallel planar cracks in a three-dimen-sional finite body subjected to arbitrary loads is derived, and thenits numerical method is proposed by the finite-part integral methodcombined with the boundary element method. According to the analytictheory of hypersingular integral equations, the square root models ofdisplacement discontinuities in the elements near the crack front areapplied, and thus the computational precision is raised. Based onthis, the stress intensity factors can be directly calculated.Finally, the stress intensity factors of sever- al typicalinteraction problems are calculated.
文摘The relation between the normal displacement on the surface of a dynamical elliptical crack and the normal stress over the crack surface was studied. The three dimensional elastodynamic equations and Fourier Laplace transforms are used. Based on the influence function and the inversion of integral transforms, one can find that if the distribution of normal displacement on the surface of a dynamic elliptical crack is a polynomial of degree n in x 1 and x 2 , then the normal pressure acting over the ellipse is also a polynomial P n(x 1,x 2) of the same degree in x 1 and x 2 .
基金the National Natural Science Foundation of China( No.K19672007)
文摘The dynamic stress intensity factor of a three-dimensionalelliptic crack under impact loading is determined with the finiteelement method. The computation results can take into account theinfluence of time and the ratio of the wave speeds on the stressintensity factor. The present method is suitable not only forthree-dimensional dynamic crack, but also for three-dimensionaldynamic contact.
文摘A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.
基金Project supported by the National Natural Science Foundation of China (No. 10172075)the Yu-Ying Foundation of Hunan University.
文摘Exact solutions in elementary functions are derived for the stress and electric displacement intensity factors of a half-plane crack in a transversely isotropic piezoelectric space interacting with various resultant sources, including force dipole, electric dipole, moment, force dilatation and rotation. Such force and charge sources may model defects like vacancies, foreign particles and dislocations. The locations and orientations of the stress and charge sources with respect to the crack are arbitrary.
基金financially supported by the National Natural Science Foundation of China(Grant No.52071343)the Leading Innovation and Entrepreneurship Team of Zhejiang Province-Automotive Light Alloy Innovation Team(No.2022R01018).
文摘Large solidification ranges and coarse columnar grains in the additively manufacturing of Al-Mg-Si alloys are normally involved in hot cracks during solidification.In this work,we develop novel crack-free Al-Mg_(2) Si alloys fabricated by laser powder-bed fusion(L-PBF).The results indicate that the eutectic Mg_(2) Si phase possesses a strong ability to reduce crack susceptibility.It can enhance the grain growth restriction factor in the initial stage of solidification and promote eutectic filling in the terminal stage of solidifica-tion.The crack-free L-PBFed Al-x Mg_(2) Si alloys(x=6 wt.%,9 wt.%,and 12 wt.%)exhibit the combination of low crack susceptibility index(CSI),superior ability for liquid filling,and grain refinement.Particularly,the L-PBFed Al-9Mg_(2) Si alloy shows improved mechanical properties(e.g.yield strength of 397 MPa and elongation of 7.3%).However,the cracks are more likely to occur in the region near the columnar grain boundaries of the L-PBFed Al-3Mg_(2) Si alloy with a large solidification range and low eutectic content for liquid filling.Correspondingly,the L-PBFed Al-3Mg_(2) Si alloy shows poor bearing capacity of mechanical properties.The precise tuning of Mg_(2) Si eutectic content can offer an innovative strategy for eliminating cracks in additively manufactured Al-Mg-Si alloy.
基金Supported by the National Natural Science Foundation of China (50979052,40872203, 41072234)the Provincial Natural Science Foundation of Shandong (ZR2009FM041,ZR2010EM032,ZR2009AZ001)
文摘The three-dimensional （3D） crack propagation is a hot issue in rock mechanics. To properly simulate 3D crack propagation, a modified maximum tangential tensile stress criterion is proposed. In this modified criterion, it is supposed that cracks propagate only at crack front in the principal normal plane. The tangential tensile stress at crack front in the principal normal plane in local coordinates is employed to determine crack propagation, which is calculated through coordinate transformation from global to local coordinates. New cracks will propagate when the maximum tangential tensile stress at crack front in the principal normal plane reaches the tensile strength of rock-like materials. Compared with the previous crack propagation criteria, the modified crack propagation criterion is helpful in calculating 3D crack stress intensity factor, and can overcome the limitations of propagation step determined by individual experiences in previous studies. Finally, the 3D crack propagation process is traced by element-free Galerkin method. The numerical results agree well with the experimental ones for a frozen resin sample with prefabricated 3D cracks.
