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.展开更多
The stability of rock slopes is frequently controlled by the initiation and propagation of inherent dominant cracks.This study systematically investigated these processes in valley slopes by combining fracture-mechani...The stability of rock slopes is frequently controlled by the initiation and propagation of inherent dominant cracks.This study systematically investigated these processes in valley slopes by combining fracture-mechanics analysis with transparent soil model tests.An analytical expression for the stress field at the dominant crack tip was derived from the slope stress distribution by superposing the corresponding stress intensity factors(SIFs).The theoretical predictions were then validated against observations from transparent soil model tests.The influences of slope angle(β),crack inclination angle(α),crack position parameter(b),and crack length parameter(h)on crack initiation and propagation were quantified.The results indicated that:(1)cracks at the slope crest tended to propagate in shear mode,and the shear crack initiation angle(θ_(s))was approximately 8°.Cracks at the slope toe might propagate in either tensile or shear mode.(2)θ_(s) at the slope crest increased withβ,b,and l,and decreased withα.The maximum change inθ_(s) induced by the considered parameters was approximately 30°.(3)The tensile crack initiation angle(θ_(t))at the slop toe decreased withβ,α,and l,while the influence of b was comparatively minor.The maximum change inθ_(t) caused by individual parameters ranged approximately from 25°to 60°.Predicted crack propagation modes and directions showed good agreement with experimental results.These findings provide theoretical guidance for stability assessments of valley slopes controlled by dominant crack propagation.展开更多
This research centers on structural health monitoring of bridges,a critical transportation infrastructure.Owing to the cumulative action of heavy vehicle loads,environmental variations,and material aging,bridge compon...This research centers on structural health monitoring of bridges,a critical transportation infrastructure.Owing to the cumulative action of heavy vehicle loads,environmental variations,and material aging,bridge components are prone to cracks and other defects,severely compromising structural safety and service life.Traditional inspection methods relying on manual visual assessment or vehicle-mounted sensors suffer from low efficiency,strong subjectivity,and high costs,while conventional image processing techniques and early deep learning models(e.g.,UNet,Faster R-CNN)still performinadequately in complex environments(e.g.,varying illumination,noise,false cracks)due to poor perception of fine cracks andmulti-scale features,limiting practical application.To address these challenges,this paper proposes CACNN-Net(CBAM-Augmented CNN),a novel dual-encoder architecture that innovatively couples a CNN for local detail extraction with a CBAM-Transformer for global context modeling.A key contribution is the dedicated Feature FusionModule(FFM),which strategically integratesmulti-scale features and focuses attention on crack regions while suppressing irrelevant noise.Experiments on bridge crack datasets demonstrate that CACNNNet achieves a precision of 77.6%,a recall of 79.4%,and an mIoU of 62.7%.These results significantly outperform several typical models(e.g.,UNet-ResNet34,Deeplabv3),confirming their superior accuracy and robust generalization,providing a high-precision automated solution for bridge crack detection and a novel network design paradigm for structural surface defect identification in complex scenarios,while future research may integrate physical features like depth information to advance intelligent infrastructure maintenance and digital twin management.展开更多
According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing ...According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.展开更多
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.展开更多
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 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.展开更多
基金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.
基金financially supported by the National Nature Science Foundation of China(Nos.52379110 and 42207222)the Key Technologies for Accurate Diagnosis and Intelligent Prevention and Control of Slope Hazards in Open Pit Mines,181 Major R&D projects of Metallurgical Corporation of China Ltd。
文摘The stability of rock slopes is frequently controlled by the initiation and propagation of inherent dominant cracks.This study systematically investigated these processes in valley slopes by combining fracture-mechanics analysis with transparent soil model tests.An analytical expression for the stress field at the dominant crack tip was derived from the slope stress distribution by superposing the corresponding stress intensity factors(SIFs).The theoretical predictions were then validated against observations from transparent soil model tests.The influences of slope angle(β),crack inclination angle(α),crack position parameter(b),and crack length parameter(h)on crack initiation and propagation were quantified.The results indicated that:(1)cracks at the slope crest tended to propagate in shear mode,and the shear crack initiation angle(θ_(s))was approximately 8°.Cracks at the slope toe might propagate in either tensile or shear mode.(2)θ_(s) at the slope crest increased withβ,b,and l,and decreased withα.The maximum change inθ_(s) induced by the considered parameters was approximately 30°.(3)The tensile crack initiation angle(θ_(t))at the slop toe decreased withβ,α,and l,while the influence of b was comparatively minor.The maximum change inθ_(t) caused by individual parameters ranged approximately from 25°to 60°.Predicted crack propagation modes and directions showed good agreement with experimental results.These findings provide theoretical guidance for stability assessments of valley slopes controlled by dominant crack propagation.
基金supported by the National Natural Science Foundation of China(No.52308332)the General Scientific Research Project of the Education Department of Zhejiang Province(No.Y202455824).
文摘This research centers on structural health monitoring of bridges,a critical transportation infrastructure.Owing to the cumulative action of heavy vehicle loads,environmental variations,and material aging,bridge components are prone to cracks and other defects,severely compromising structural safety and service life.Traditional inspection methods relying on manual visual assessment or vehicle-mounted sensors suffer from low efficiency,strong subjectivity,and high costs,while conventional image processing techniques and early deep learning models(e.g.,UNet,Faster R-CNN)still performinadequately in complex environments(e.g.,varying illumination,noise,false cracks)due to poor perception of fine cracks andmulti-scale features,limiting practical application.To address these challenges,this paper proposes CACNN-Net(CBAM-Augmented CNN),a novel dual-encoder architecture that innovatively couples a CNN for local detail extraction with a CBAM-Transformer for global context modeling.A key contribution is the dedicated Feature FusionModule(FFM),which strategically integratesmulti-scale features and focuses attention on crack regions while suppressing irrelevant noise.Experiments on bridge crack datasets demonstrate that CACNNNet achieves a precision of 77.6%,a recall of 79.4%,and an mIoU of 62.7%.These results significantly outperform several typical models(e.g.,UNet-ResNet34,Deeplabv3),confirming their superior accuracy and robust generalization,providing a high-precision automated solution for bridge crack detection and a novel network design paradigm for structural surface defect identification in complex scenarios,while future research may integrate physical features like depth information to advance intelligent infrastructure maintenance and digital twin management.
基金support of the National Natural Science Foundation of China(No.12172023)。
文摘According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.
基金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.
基金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(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.
