Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classic...Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classical plane elastic crack model, only the known conditions were revised in the new formulation, which are greatly convenient to solve the problem, and no other new condition was given. Results and Conclusion The general exact analytic solution is given here based on the formulation though the problem is very complicated. Furthermore, the stress intensity factors K Ⅰ, K Ⅱ of the problem are also given.展开更多
The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed...The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.展开更多
Industry design standards such as BS 7910 deployed some empirical formulas for the prediction of stress intensity factor(SIF) based on simulation results from traditional finite element method(FEM).However,such FEM si...Industry design standards such as BS 7910 deployed some empirical formulas for the prediction of stress intensity factor(SIF) based on simulation results from traditional finite element method(FEM).However,such FEM simulation occasionally failed to convince people due to the large discrepancies compared with engineering practice.As a consequence,inaccuracy predictions via such formulas in engineering standards inevitably occur,which will compromise the safety of structures.In our previous research work,an abnormal phenomenon of SIF in a cracked T-butt joint accounting for welding effect has been observed.Compared with BS 7910,the calculation results of SIF at the surface points of welded specimens cannot be well predicted,with a large discrepancy appearing.In order to explore such problem with an abnormal increase at the surface points of cracked welded specimens,a numerical investigation in terms of SIF among BS 7910,XFEM,and FEM is performed in this paper.Numerical models on both a simple cracked plate without welding effect and a cracked T-butt joint with welding effect are developed through ABAQUS.Parametric studies in terms of the effects of varied crack depth to thickness ratio(a/T) and the effects of crack depth to crack half-length ratio(a/c) are carried out.Empirical solutions from BS 7910 are used for comparison.It is found that the XFEM can provide predictions of SIF at both the crack deepest point and crack surface point of a simple cracked plate as accurate as FEM.For a T-butt joint with a transverse stiffener,a large discrepancy in terms of the weld magnification factors(Mk) occurs at the crack surface point compared with empirical predictions.An exceptional increase of von Mises stress gradient in regions close to the weld-toe is found through the simulation of FEM,whereas a constant stress gradient is obtained through XFEM.The comparison results indicate an inappropriate prediction of SIF by the utilization of the empirical formulas in BS 7910.A more reasonable prediction of the SIF at the surface point of a crack is obtained by the XFEM.Therefore,further updating of the empirical solutions in BS7910 for SIF accounting for welding effect is recommended.展开更多
Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equation...Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.展开更多
Based on the theory of Muskhelishvili, the general solutions for stress and strain of conjugate cracks in cubic quasicrystal are obtained, with which the stress intensity factors of cubic quasicrystal at crack tips an...Based on the theory of Muskhelishvili, the general solutions for stress and strain of conjugate cracks in cubic quasicrystal are obtained, with which the stress intensity factors of cubic quasicrystal at crack tips and the stress distribution functions of phonon and phason fields are given. The results show that though phason field is coupled with phonon field by constitutive equations, the stress intensity factors are not coupled with any other factors.展开更多
Within the framework of nonlinear eleetroelasticity, the stress field near to the crack tip in an infinite piezoelectric media subject to a far field uniform loading is studied by using an electrical strip saturation ...Within the framework of nonlinear eleetroelasticity, the stress field near to the crack tip in an infinite piezoelectric media subject to a far field uniform loading is studied by using an electrical strip saturation model and the complex variable method. And the emphasis is placed on the stress field near to the crack tip. The obtained solutions show that the normalized stress components at an arbitrary point near to the crack tip are determined by the angle of the point. Moreover, the stress components are independent of the distance from the point to the ori- gin of the coordinate. The distributions of in-plane stress components near to the crack tip are analyzed based on numerical results for PZT-SH. Compared with some related solutions, results show that the solutions are valid.展开更多
The fracture theory of cubic quasicrystal was developed. The exact analytic solution of a Mode Ⅲ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, and so t...The fracture theory of cubic quasicrystal was developed. The exact analytic solution of a Mode Ⅲ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, and so the displacement and stress fields, the stress intensity factor and strain energy release rate were determined. The results show that the stress intensity factor is independent of material constants, and the strain energy release rate is dependent on all material constants. These provide important information for studying the deformation and fracture of the new solid material.展开更多
This study investigates the interaction and influence of surface cracks on the spherical pressure hull of a deep-sea manned submersible.The finite element model of the spherical hull is established,and a semi-elliptic...This study investigates the interaction and influence of surface cracks on the spherical pressure hull of a deep-sea manned submersible.The finite element model of the spherical hull is established,and a semi-elliptical surface crack is inserted in the welding toe of the spherical hull as the main crack.Considering the combined effect of external uniform pressure and welding residual stress at the weld toe,the stress intensity factor(SIF)is obtained based on the M-integral method.Inserting disturbing cracks at different positions on the spherical hull surface,the interaction and influence between multi-cracks are revealed by numerical calculation.The results show that the existence of the disturbing crack has a great influence on the stress intensity factor of the main crack,and the influence is different with the different location of disturbing crack.The study of the interaction of multiple cracks under different interference factors and the influence of disturbing cracks on the main crack can provide some reference for future engineering applications.展开更多
Based on the finite element method, a numerical investigation into the bonded repair efficiency of cracked plates under in plane biaxial loadings is presented. The main considerations are: reduction in stress intensi...Based on the finite element method, a numerical investigation into the bonded repair efficiency of cracked plates under in plane biaxial loadings is presented. The main considerations are: reduction in stress intensity factor (SIF) at the crack tip, the maximum tensile stress in the composite patch and the maximum shear stress in the adhesive bond between the patch and the plate. Without the patch, a tensile or compressive stress parallel to the crack has no effect on the SIF at the crack tip. While with a composite patch, there exists coupling effect between the normal stress parallel to the crack and the SIF, and the coupling effect depends significantly on ply orientation of the patch and the biaxial stress ratio of the plate.展开更多
Using the complex variable function method and the technique of conformal mapping, the anti-plane shear problem of an elliptic hole with asymmetric colfinear cracks in a one-dimensional hexagonal quasi-crystal is solv...Using the complex variable function method and the technique of conformal mapping, the anti-plane shear problem of an elliptic hole with asymmetric colfinear cracks in a one-dimensional hexagonal quasi-crystal is solved, and the exact analytic solutions of the stress intensity factors （SIFs） for mode Ⅲ problem are obtained. Under the limiting conditions, the present results reduce to the Griffith crack and many new results obtained as well, such as the circular hole with asymmetric collinear cracks, the elliptic hole with a straight crack, the mode T crack, the cross crack and so on. As far as the phonon field is concerned, these results, which play an important role in many practical and theoretical applications, are shown to be in good agreement with the classical results.展开更多
The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack grow...The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack growth is proposed in this study.By adopting the proposed method,the potential collapse location of strata is derivable in accordance with a static model,the durability of roof strata can be estimated,a dynamic time step control strategy is achieved to balance the accuracy and speed of computing,and the initial crack size of rock can be estimated.In addition to the above,a mechanical model of underground excavation with non-uniformly distributed loads and partially yielded foundation is presented as the prototypical case.A set of case studies is carried out that showcase a power correlation between applied stress and roof durability.The allowable applied tensile stress for a 100-year life cycle is about 76%of the tensile strength.By using the proposed subcritical crack growth computation scheme,the roof stability in an underground excavation can be identified not only from the spatial view but also from the temporal perspective.展开更多
In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under ...In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.展开更多
A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to c...A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.展开更多
Wave propagation in an infinite elastic piezoelectric medium with a circular cavity and an impermeable crack subjected to steady-state anti-plane shearing was studied based on Green's function and the crack-divisi...Wave propagation in an infinite elastic piezoelectric medium with a circular cavity and an impermeable crack subjected to steady-state anti-plane shearing was studied based on Green's function and the crack-division technique.Theoretical solutions were derived for the whole elastic displacement and electric potential field in the interaction between the circular cavity and the impermeable crack.Expressions were obtained on the dynamic stress concentration factor(DSCF) at the cavity's edge,the dynamic stress intensity factor(DSIF) and the dynamic electric displacement intensity factor(DEDIF) at the crack tip.Numerical solutions were performed and plotted with different incident wave numbers,parameters of piezoelectric materials and geometries of the structure.Finally,some of the calculation results were compared with the case of dynamic anti-plane interaction of a permeable crack and a circular cavity in an infinite piezoelectric medium.This paper can provide a valuable reference for the design of piezoelectric actuators and sensors widely used in marine structures.展开更多
A research on the stable fatigue crack propagation of 16MnR steel is investigated systematically in this paper.First,control experiments of 16MnR with compact tension specimen is conducted to study the effect of R-rat...A research on the stable fatigue crack propagation of 16MnR steel is investigated systematically in this paper.First,control experiments of 16MnR with compact tension specimen is conducted to study the effect of R-ratios,specimen thickness and notch sizes.The experiments show that the fatigue crack growth(FCG)rate in stable propagation was insensitive to these factors.Then,the stress intensity factor(SIF)is computed and compared by displacement interpolation method,J integral and interaction integral method respectively.The simulation shows that optimization on the mesh density and the angle of singular element improved the computational efficiency and accuracy of SIF and the interaction integral method has an obvious advantage on stability.Finally,the FCG rate is modeled by the Jiang fatigue damage criterion and the extended finite element method(XFEM)respectively.The simulation results of FCG rate are in line with experiments data and indicate that XFEM method is more accurate than Jiang fatigue damage method.展开更多
Using the boundary integral equation method, the problem of an external circular crack in a three_dimensional infinite elastic body under asymmetric loadings is investigated. The two_dimensional singular boundary inte...Using the boundary integral equation method, the problem of an external circular crack in a three_dimensional infinite elastic body under asymmetric loadings is investigated. The two_dimensional singular boundary integral equations of the problem were reduced to a system of Abel integral equations by means of Fourier series and hypergeometric functions. The exact solutions of stress intensity factors are obtained for the problem of an external circular crack under asymmetric loadings, which are even more universal than the results obtained by the use of Hankel transform method. The results demonstrate that the boundary integral equation method has great potential as a new analytic method.展开更多
文摘Aim The general arbitrary cracked problem in an elastic plane was discussed. Methods For the purpose of acquiring the solution of the problem, a new formulation on the problem was proposed. Compared with the classical plane elastic crack model, only the known conditions were revised in the new formulation, which are greatly convenient to solve the problem, and no other new condition was given. Results and Conclusion The general exact analytic solution is given here based on the formulation though the problem is very complicated. Furthermore, the stress intensity factors K Ⅰ, K Ⅱ of the problem are also given.
文摘The dynamic stress intensity factor (DSIF) and the scattering of SH wave by circle canyon and crack are studied with Green's function. In order to solve the problem, a suitable Green's function is constructed first, which is the solution of displacement fields for elastic half space with circle canyon under output plane harmonic line loading at horizontal surface. Then the integral equation for determining the unknown forces in the problem can be changed into the algebraic one and solved numerically so that crack DSIF can be determined. Last when the medium parameters are altered, the influence on the crack DSIF is discussed partially with the displacement between circle canyon and crack.
文摘Industry design standards such as BS 7910 deployed some empirical formulas for the prediction of stress intensity factor(SIF) based on simulation results from traditional finite element method(FEM).However,such FEM simulation occasionally failed to convince people due to the large discrepancies compared with engineering practice.As a consequence,inaccuracy predictions via such formulas in engineering standards inevitably occur,which will compromise the safety of structures.In our previous research work,an abnormal phenomenon of SIF in a cracked T-butt joint accounting for welding effect has been observed.Compared with BS 7910,the calculation results of SIF at the surface points of welded specimens cannot be well predicted,with a large discrepancy appearing.In order to explore such problem with an abnormal increase at the surface points of cracked welded specimens,a numerical investigation in terms of SIF among BS 7910,XFEM,and FEM is performed in this paper.Numerical models on both a simple cracked plate without welding effect and a cracked T-butt joint with welding effect are developed through ABAQUS.Parametric studies in terms of the effects of varied crack depth to thickness ratio(a/T) and the effects of crack depth to crack half-length ratio(a/c) are carried out.Empirical solutions from BS 7910 are used for comparison.It is found that the XFEM can provide predictions of SIF at both the crack deepest point and crack surface point of a simple cracked plate as accurate as FEM.For a T-butt joint with a transverse stiffener,a large discrepancy in terms of the weld magnification factors(Mk) occurs at the crack surface point compared with empirical predictions.An exceptional increase of von Mises stress gradient in regions close to the weld-toe is found through the simulation of FEM,whereas a constant stress gradient is obtained through XFEM.The comparison results indicate an inappropriate prediction of SIF by the utilization of the empirical formulas in BS 7910.A more reasonable prediction of the SIF at the surface point of a crack is obtained by the XFEM.Therefore,further updating of the empirical solutions in BS7910 for SIF accounting for welding effect is recommended.
文摘Using the single crack solution and the regular solution of plane harmonic function, the problem of Saint_Venant bending of a cracked cylinder by a transverse force was reduced to solving two sets of integral equations and its general solution was then obtained. Based on the obtained solution, a method to calculate the bending center and the stress intensity factors of the cracked cylinger whose cross_section is not thin_walled, but of small torsion rigidity is proposed. Some numerical examples are given.
基金Science Research Foundation of Shang-hai in China (No.2000SG31& 2004096)Shanghai Leading Academic Discipline Project (No.T0601)
文摘Based on the theory of Muskhelishvili, the general solutions for stress and strain of conjugate cracks in cubic quasicrystal are obtained, with which the stress intensity factors of cubic quasicrystal at crack tips and the stress distribution functions of phonon and phason fields are given. The results show that though phason field is coupled with phonon field by constitutive equations, the stress intensity factors are not coupled with any other factors.
基金the Innovation Project for Graduates in Jiangsu Province~~
文摘Within the framework of nonlinear eleetroelasticity, the stress field near to the crack tip in an infinite piezoelectric media subject to a far field uniform loading is studied by using an electrical strip saturation model and the complex variable method. And the emphasis is placed on the stress field near to the crack tip. The obtained solutions show that the normalized stress components at an arbitrary point near to the crack tip are determined by the angle of the point. Moreover, the stress components are independent of the distance from the point to the ori- gin of the coordinate. The distributions of in-plane stress components near to the crack tip are analyzed based on numerical results for PZT-SH. Compared with some related solutions, results show that the solutions are valid.
文摘The fracture theory of cubic quasicrystal was developed. The exact analytic solution of a Mode Ⅲ Griffith crack in the material was obtained by using the Fourier transform and dual integral equations theory, and so the displacement and stress fields, the stress intensity factor and strain energy release rate were determined. The results show that the stress intensity factor is independent of material constants, and the strain energy release rate is dependent on all material constants. These provide important information for studying the deformation and fracture of the new solid material.
基金This work was supported by the State Key Program of National Natural Science of China(Grant No.51439004)the Natural Science Foundation of the Jiangsu Higher Education Institutions of China(Grant No.19KJA530002)the Open Project Foundation of State Key Laboratory of Ocean Engineering(No.2006).
文摘This study investigates the interaction and influence of surface cracks on the spherical pressure hull of a deep-sea manned submersible.The finite element model of the spherical hull is established,and a semi-elliptical surface crack is inserted in the welding toe of the spherical hull as the main crack.Considering the combined effect of external uniform pressure and welding residual stress at the weld toe,the stress intensity factor(SIF)is obtained based on the M-integral method.Inserting disturbing cracks at different positions on the spherical hull surface,the interaction and influence between multi-cracks are revealed by numerical calculation.The results show that the existence of the disturbing crack has a great influence on the stress intensity factor of the main crack,and the influence is different with the different location of disturbing crack.The study of the interaction of multiple cracks under different interference factors and the influence of disturbing cracks on the main crack can provide some reference for future engineering applications.
文摘Based on the finite element method, a numerical investigation into the bonded repair efficiency of cracked plates under in plane biaxial loadings is presented. The main considerations are: reduction in stress intensity factor (SIF) at the crack tip, the maximum tensile stress in the composite patch and the maximum shear stress in the adhesive bond between the patch and the plate. Without the patch, a tensile or compressive stress parallel to the crack has no effect on the SIF at the crack tip. While with a composite patch, there exists coupling effect between the normal stress parallel to the crack and the SIF, and the coupling effect depends significantly on ply orientation of the patch and the biaxial stress ratio of the plate.
基金supported by the National Natural Science Foundation of China (Grant No 10761005)the Inner Mongolia Natural Science Foundation of China (Grant No 200607010104)
文摘Using the complex variable function method and the technique of conformal mapping, the anti-plane shear problem of an elliptic hole with asymmetric colfinear cracks in a one-dimensional hexagonal quasi-crystal is solved, and the exact analytic solutions of the stress intensity factors （SIFs） for mode Ⅲ problem are obtained. Under the limiting conditions, the present results reduce to the Griffith crack and many new results obtained as well, such as the circular hole with asymmetric collinear cracks, the elliptic hole with a straight crack, the mode T crack, the cross crack and so on. As far as the phonon field is concerned, these results, which play an important role in many practical and theoretical applications, are shown to be in good agreement with the classical results.
基金China Scholarship Council(CSC)The University of Queensland for a Ph D fellowship。
文摘The long-term stability of the roof is particularly important in designing underground rock structures.To estimate the durability of roof strata in underground excavation,a computation scheme of subcritical crack growth is proposed in this study.By adopting the proposed method,the potential collapse location of strata is derivable in accordance with a static model,the durability of roof strata can be estimated,a dynamic time step control strategy is achieved to balance the accuracy and speed of computing,and the initial crack size of rock can be estimated.In addition to the above,a mechanical model of underground excavation with non-uniformly distributed loads and partially yielded foundation is presented as the prototypical case.A set of case studies is carried out that showcase a power correlation between applied stress and roof durability.The allowable applied tensile stress for a 100-year life cycle is about 76%of the tensile strength.By using the proposed subcritical crack growth computation scheme,the roof stability in an underground excavation can be identified not only from the spatial view but also from the temporal perspective.
基金Supported by the Opening Fund of the State Key Laboratory of Explosion Science and Technology of Beijing Institute of Technology(KFJJ17-01M)
文摘In this paper,dynamic caustic method is applied to analyze the blast-induced crack propagation and distribution of the dynamic stress field around an empty circular hole in polymethyl methacrylate(PMMA)material under delayed ignition blasting loads.The following experimental results are obtained.(1)In directional-fracture-controlled blasting,the dynamic stress intensity factors(DSIFs)and the propagation paths of the blast-induced cracks are obviously influenced by the delayed ignition.(2) The circular hole situated between the two boreholes poses a strong guiding effect on the coelesence of the cracks,causing them to propagate towards each other when cracks are reaching the circular hole area.(3)Blast-induced cracks are not initiated preferentially because of the superimposed effect from the explosive stress waves on the cracking area.(4) By using the scanning electron microscopy(SEM)method,it is verified that the roughness of crack surfaces changes along the crack propagation paths.
基金supported by the National Basic Research Program of China(Grant No.2012CB026200)the National Natural Science Foundation of China(Grant No.50878048)
文摘A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.
基金Supported by the Natural Science Foundation of Heilongjiang Province of China (A00-10)the Basis Research Foundation of Harbin Engineering University (HEUF04008)
文摘Wave propagation in an infinite elastic piezoelectric medium with a circular cavity and an impermeable crack subjected to steady-state anti-plane shearing was studied based on Green's function and the crack-division technique.Theoretical solutions were derived for the whole elastic displacement and electric potential field in the interaction between the circular cavity and the impermeable crack.Expressions were obtained on the dynamic stress concentration factor(DSCF) at the cavity's edge,the dynamic stress intensity factor(DSIF) and the dynamic electric displacement intensity factor(DEDIF) at the crack tip.Numerical solutions were performed and plotted with different incident wave numbers,parameters of piezoelectric materials and geometries of the structure.Finally,some of the calculation results were compared with the case of dynamic anti-plane interaction of a permeable crack and a circular cavity in an infinite piezoelectric medium.This paper can provide a valuable reference for the design of piezoelectric actuators and sensors widely used in marine structures.
文摘A research on the stable fatigue crack propagation of 16MnR steel is investigated systematically in this paper.First,control experiments of 16MnR with compact tension specimen is conducted to study the effect of R-ratios,specimen thickness and notch sizes.The experiments show that the fatigue crack growth(FCG)rate in stable propagation was insensitive to these factors.Then,the stress intensity factor(SIF)is computed and compared by displacement interpolation method,J integral and interaction integral method respectively.The simulation shows that optimization on the mesh density and the angle of singular element improved the computational efficiency and accuracy of SIF and the interaction integral method has an obvious advantage on stability.Finally,the FCG rate is modeled by the Jiang fatigue damage criterion and the extended finite element method(XFEM)respectively.The simulation results of FCG rate are in line with experiments data and indicate that XFEM method is more accurate than Jiang fatigue damage method.
文摘Using the boundary integral equation method, the problem of an external circular crack in a three_dimensional infinite elastic body under asymmetric loadings is investigated. The two_dimensional singular boundary integral equations of the problem were reduced to a system of Abel integral equations by means of Fourier series and hypergeometric functions. The exact solutions of stress intensity factors are obtained for the problem of an external circular crack under asymmetric loadings, which are even more universal than the results obtained by the use of Hankel transform method. The results demonstrate that the boundary integral equation method has great potential as a new analytic method.
