In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors (SIFs) of cracks emanating from a surface semi-spherical cavity in a finite body is derived us...In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors (SIFs) of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.展开更多
Because of the wicked service environment of the high speed train, it is possible that the hollow axle of the train may encounter the foreign object damage and form a sharp notch. Under the fatigue loading a crack can...Because of the wicked service environment of the high speed train, it is possible that the hollow axle of the train may encounter the foreign object damage and form a sharp notch. Under the fatigue loading a crack can initiate from the notch and propagate to failure. It is noted that the stress intensity factor is the control parameter of the crack propagating, for the purpose of getting the more exact propagation characteristics, the stress intensity factor is studied mainly. The service loads of hollow axles are defined, and the stress distribution of hollow axles is obtained according to the load spectrum. The semi-ellipse crack configuration is defined with three parameters: the aspect ratio, the relative depth and the relative location along the crack front. Quarter point 20-node isoparametric degenerate singular elements are used for the region near the crack tip. The finite element model of crack extension of hollow axle is created, and the crack front is dispersed which can realize orthogonal extension. Based on this the stress intensity factors of crack front were calculated, and the distribution rules of the stress intensity factors of different initial crack shapes are obtained. The conclusions are compared with that of the analytic method and they agree with each other very well, and the calculating results show that there is a close relationship between the stress intensity factor and the initial crack shape. For a round crack the stress intensity factor at the surface point increases faster than the one at the center point with the crock propagation. However, for a narrow crack, the results are in contrast with that of a round one. So, all the cracks with different shapes propagate toward to a similar shape, and they grow at this shape to end. The study may contribute to the crack propagate characteristics research.展开更多
The hatches for inspecting are usually designed with elliptical holes in airplane structures, so computation of the stress intensity factor of three dimensional crack at elliptical holes is pivotal for damage toleranc...The hatches for inspecting are usually designed with elliptical holes in airplane structures, so computation of the stress intensity factor of three dimensional crack at elliptical holes is pivotal for damage tolerance analysis of these structures. In this paper, weight function is derived for a two dimensional through cracks at elliptical holes by applying a compounding method. Stress intensity factor formulas for an internal surface semi-elliptical crack in elliptical holes are obtained wing the three dimensional weight function method. Stress intensity factors for an internal surface semi-elliptical crack in elliptical holes under remote tension are computed. At the same time, research on how radius of curvature for elliptical holes affect stress intensity factors was conducted. Stress intensity factors decrease when radius of curvature increases. Some results and conclusions which are of practical value are given.展开更多
The method of complex function and the method of Green's function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an...The method of complex function and the method of Green's function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an elliptical hole, and the solution of dynamic stress intensity factor at the crack tip was given. A Green's function was constructed for the problem, which is a basic solution of displacement field for an elastic half space containing a half elliptical gap impacted by anti-plane harmonic linear source force at any point of its horizontal boundary. With division of a crack technique, a series of integral equations can be established on the conditions of continuity and the solution of dynamic stress intensity factor can be obtained. The influence of an elliptical hole on the dynamic stress intensity factor at the crack tip was discussed.展开更多
In this work,the stress intensity factors(SIFs)for a wide range of semi-elliptical surface cracks with different inclination angles in a cylindrical bar subjected to pure tension were investigated numerically.The main...In this work,the stress intensity factors(SIFs)for a wide range of semi-elliptical surface cracks with different inclination angles in a cylindrical bar subjected to pure tension were investigated numerically.The main parameters considered to evaluate the SIFs were the crack size ratio(a/d),the crack aspect ratio(a/c),and the crack inclination angle(0).The dual-boundary element met hod implemented in software BEASY was used to compute the SIF values for cracks.Furt her more,the general closed-form solutions were proposed to evaluate the corresponding SIFs for Mode I,Mode II and Mode III types of fracture through curve fitting approaches.These solutions can provide accurate and reliable values of stress intensity factors for a crack on a cylindrical bar under pure tension in a rapid way compared to those obtained using computational models.In addition,these results can be used as inputs for failure studies and life evaluations of cracked cylinder under working conditions.展开更多
In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a ...In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid's boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.展开更多
The authors have developed a new line-spring boundary element method in the present paper, which combines the advantage of the line-spring model with that of the boundary element method. This method reduces the three-...The authors have developed a new line-spring boundary element method in the present paper, which combines the advantage of the line-spring model with that of the boundary element method. This method reduces the three-dimension problem of the surface cracks into a quasi-one-dimension problem and can be used to analyze the surface cracked plate under various loading conditions. In this paper theoretical analyses and numerical verifications are carried out. The calculated results are reported, which indicate that the present method is efficient and can be used to analyze the surface crack problem on a personal computer.展开更多
Calculating interacting stress intensity factors(SIFs)of multiple ellipticalholes and cracks is very important for safety assessment,stop-hole optimization design and resource exploitation production in underground ro...Calculating interacting stress intensity factors(SIFs)of multiple ellipticalholes and cracks is very important for safety assessment,stop-hole optimization design and resource exploitation production in underground rock engineering,e.g.,buried tunnels,deep mining,geothermal and shale oil/gas exploitation by hydraulic fracturing technology,where both geo-stresses and surface stresses are applied on buried tunnels,horizontal wells and natural cracks.However,current literatures are focused mainly on study of interacting SIFs of multiple elliptical-holes(or circularholes)and cracks only under far-field stresses without consideration of arbitrary surface stresses.Recently,our group has proposed a new integral method to calculate interacting SIFs of multiple circular-holes and cracks subjected to far-filed and surface stresses.This new method will be developed to study the problem of multiple elliptical-hole and cracks subjected to both far-field and surface stresses.In this study,based on Cauchy integral theorem,the exact fundamental stress solutions of single elliptical-hole under arbitrarily concentrated surface normal and shear forces are derived to establish new integral equation formulations for calculating interacting SIFs of multiple elliptical-holes and cracks under both far-field and arbitrary surface stresses.The new method is proved to be valid by comparing our results of interacting SIFs with those obtained by Green’s function method,displacement discontinuity method,singular integral equation method,pseudo-dislocations method and finite element method.Computational examples of one elliptical-hole and one crack in an infinite elastic body are given to analyze influence of loads and geometries on interacting SIFs.Research results show that whenσ_(xx)^(∞)≥σ^(yy′)^(∞),there appears a neutral crack orientation angle b0(without elliptical-hole’s effect).Increasing s¥xx/s¥yy and b/a(close to circularhole)usually decreases b0 of KI and benefits to the layout of stop-holes.The surface compressive stresses applied onto elliptical-hole(n)and crack(p)have significant influence on interacting SIFs but almost no on b0.Increasing n and p usually results in increase of interacting SIFs and facilitates crack propagation and fracture networks.The elliptical-hole orientation angle(a)and holed-cracked distance(t)have great influence on the interacting SIFs while have little effect on b0.The present method is not only simple(without any singular parts),high-accurate(due to exact fundamental stress solutions)and wider applicable(under far-field stresses and arbitrarily distributed surface stress)than the common methods,but also has the potential for the anisotropic problem involving multiple holes and cracks.展开更多
In the present work the fatigue crack growth in AISI304 specimens is investigated experimentally. In 3D finite element analysis the virtual crack closure technique is applied to calculate distributions and variations ...In the present work the fatigue crack growth in AISI304 specimens is investigated experimentally. In 3D finite element analysis the virtual crack closure technique is applied to calculate distributions and variations of the stress intensity factor along the surface crack front. It is confirmed that the stress intensity factor along the surface crack front varies non-uniformly with crack growth. Crack growth rate is proportional to the stress intensity factor distribution in the 3D cracked specimen. The fatigue crack growth in surface cracked specimens can be described by the Forman model identified in conventional compact tension specimens. For crack growth in the free specimen surface the arc length seems more suitable to quantify crack progress. Geometry and loading configuration of the surface cracked specimen seem to not affect the fatigue crack growth substantially.展开更多
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.展开更多
In order to perfectly reflect the dynamic corrosion of reinforced concrete (RC) cover in practical engineering,an analytic model of non-uniform corrosion induced cracking was presented based on the elastic-plastic fra...In order to perfectly reflect the dynamic corrosion of reinforced concrete (RC) cover in practical engineering,an analytic model of non-uniform corrosion induced cracking was presented based on the elastic-plastic fracture mechanics theory.Comparisons with the published experimental data show that the predictions given by the present model are in good agreement with the results both for natural exposed experiments and short-time indoor tests (the best difference is about 2.7%).Also it obviously provides much better precision than those models under the assumption of uniform corrosion (the maximal improved precision is about 48%).Therefore,it is pointed out that the so-called uniform corrosion models to describe the cover cracking of RC should be adopted cautiously.Finally,the influences of thickness of local rusty layer around the reinforcing steel bar on the critical corrosion-induced crack indexes were investigated.It is found that the thickness of local rusty layer has great effect on the critical mass loss of reinforcing steel,threshold expansion pressure,and time to cover cracking.For local rusty layer thickness with a size of a=0.5 mm,the time to cover cracking will increase by about one times when a/b (a,semi-minor axis;b,semi-major axis) changes from 0.1 to 1 mm.展开更多
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.展开更多
This paper discusses stress intensity factor (SIF) calculations for surface cracks in round bars subjected to combined torsion and bending loadings. Different crack aspect ratios, a/b, ranging from 0.0 to 1.2 and re...This paper discusses stress intensity factor (SIF) calculations for surface cracks in round bars subjected to combined torsion and bending loadings. Different crack aspect ratios, a/b, ranging from 0.0 to 1.2 and relative crack depths, a/D, ranging from 0.1 to 0.6 were considered. Since the loading was non-symmetrical for torsion loadings, a whole finite element model was constructed. Then, the individual and combined bending and torsion loadings were remotely applied to the model. The equivalent SIF method, F*EQ, was then used explicitly to combine the individual SIFs from the bending and torsion loadings. A comparison was then carried out with the combined SIF, F*FE, obtained using the finite element analysis (FEA) under similar loadings. It was found that the equivalent SIF method successfully predicted the combined SIF for Mode I. However, discrepancies between the results determined from the different approaches occurred when Fm was involved. It was also noted that the predicted F* EE using FEA was higher than the F*EQ predicted through the equivalent SIF method due to the difference in crack face interactions.展开更多
文摘In this paper, a new semi-analytical and semi-engineering method of the closed form solution of stress intensity factors (SIFs) of cracks emanating from a surface semi-spherical cavity in a finite body is derived using the energy release rate theory. A mode of crack opening displacements of a normal slice is established, and the normal slice relevant functions are introduced. The proposed method is both effective and accurate for the problem of three-dimensional cracks emanating from a surface cavity. A series of useful results of SIFs are obtained.
基金supported by National Basic Research and Development Program of China (973 Program, Grant No. 2007CB714705)
文摘Because of the wicked service environment of the high speed train, it is possible that the hollow axle of the train may encounter the foreign object damage and form a sharp notch. Under the fatigue loading a crack can initiate from the notch and propagate to failure. It is noted that the stress intensity factor is the control parameter of the crack propagating, for the purpose of getting the more exact propagation characteristics, the stress intensity factor is studied mainly. The service loads of hollow axles are defined, and the stress distribution of hollow axles is obtained according to the load spectrum. The semi-ellipse crack configuration is defined with three parameters: the aspect ratio, the relative depth and the relative location along the crack front. Quarter point 20-node isoparametric degenerate singular elements are used for the region near the crack tip. The finite element model of crack extension of hollow axle is created, and the crack front is dispersed which can realize orthogonal extension. Based on this the stress intensity factors of crack front were calculated, and the distribution rules of the stress intensity factors of different initial crack shapes are obtained. The conclusions are compared with that of the analytic method and they agree with each other very well, and the calculating results show that there is a close relationship between the stress intensity factor and the initial crack shape. For a round crack the stress intensity factor at the surface point increases faster than the one at the center point with the crock propagation. However, for a narrow crack, the results are in contrast with that of a round one. So, all the cracks with different shapes propagate toward to a similar shape, and they grow at this shape to end. The study may contribute to the crack propagate characteristics research.
文摘The hatches for inspecting are usually designed with elliptical holes in airplane structures, so computation of the stress intensity factor of three dimensional crack at elliptical holes is pivotal for damage tolerance analysis of these structures. In this paper, weight function is derived for a two dimensional through cracks at elliptical holes by applying a compounding method. Stress intensity factor formulas for an internal surface semi-elliptical crack in elliptical holes are obtained wing the three dimensional weight function method. Stress intensity factors for an internal surface semi-elliptical crack in elliptical holes under remote tension are computed. At the same time, research on how radius of curvature for elliptical holes affect stress intensity factors was conducted. Stress intensity factors decrease when radius of curvature increases. Some results and conclusions which are of practical value are given.
文摘The method of complex function and the method of Green's function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an elliptical hole, and the solution of dynamic stress intensity factor at the crack tip was given. A Green's function was constructed for the problem, which is a basic solution of displacement field for an elastic half space containing a half elliptical gap impacted by anti-plane harmonic linear source force at any point of its horizontal boundary. With division of a crack technique, a series of integral equations can be established on the conditions of continuity and the solution of dynamic stress intensity factor can be obtained. The influence of an elliptical hole on the dynamic stress intensity factor at the crack tip was discussed.
文摘In this work,the stress intensity factors(SIFs)for a wide range of semi-elliptical surface cracks with different inclination angles in a cylindrical bar subjected to pure tension were investigated numerically.The main parameters considered to evaluate the SIFs were the crack size ratio(a/d),the crack aspect ratio(a/c),and the crack inclination angle(0).The dual-boundary element met hod implemented in software BEASY was used to compute the SIF values for cracks.Furt her more,the general closed-form solutions were proposed to evaluate the corresponding SIFs for Mode I,Mode II and Mode III types of fracture through curve fitting approaches.These solutions can provide accurate and reliable values of stress intensity factors for a crack on a cylindrical bar under pure tension in a rapid way compared to those obtained using computational models.In addition,these results can be used as inputs for failure studies and life evaluations of cracked cylinder under working conditions.
基金Project supported by the National Natural Science Foundation of China (No.10176003).
文摘In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid's boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.
文摘The authors have developed a new line-spring boundary element method in the present paper, which combines the advantage of the line-spring model with that of the boundary element method. This method reduces the three-dimension problem of the surface cracks into a quasi-one-dimension problem and can be used to analyze the surface cracked plate under various loading conditions. In this paper theoretical analyses and numerical verifications are carried out. The calculated results are reported, which indicate that the present method is efficient and can be used to analyze the surface crack problem on a personal computer.
基金supports by National Natural Science Foundation of China(Nos.51874351,51474251 and 12072309)Excellent Postdoctoral Innovative Talents Project of Hunan Province(No.2020RC2001).
文摘Calculating interacting stress intensity factors(SIFs)of multiple ellipticalholes and cracks is very important for safety assessment,stop-hole optimization design and resource exploitation production in underground rock engineering,e.g.,buried tunnels,deep mining,geothermal and shale oil/gas exploitation by hydraulic fracturing technology,where both geo-stresses and surface stresses are applied on buried tunnels,horizontal wells and natural cracks.However,current literatures are focused mainly on study of interacting SIFs of multiple elliptical-holes(or circularholes)and cracks only under far-field stresses without consideration of arbitrary surface stresses.Recently,our group has proposed a new integral method to calculate interacting SIFs of multiple circular-holes and cracks subjected to far-filed and surface stresses.This new method will be developed to study the problem of multiple elliptical-hole and cracks subjected to both far-field and surface stresses.In this study,based on Cauchy integral theorem,the exact fundamental stress solutions of single elliptical-hole under arbitrarily concentrated surface normal and shear forces are derived to establish new integral equation formulations for calculating interacting SIFs of multiple elliptical-holes and cracks under both far-field and arbitrary surface stresses.The new method is proved to be valid by comparing our results of interacting SIFs with those obtained by Green’s function method,displacement discontinuity method,singular integral equation method,pseudo-dislocations method and finite element method.Computational examples of one elliptical-hole and one crack in an infinite elastic body are given to analyze influence of loads and geometries on interacting SIFs.Research results show that whenσ_(xx)^(∞)≥σ^(yy′)^(∞),there appears a neutral crack orientation angle b0(without elliptical-hole’s effect).Increasing s¥xx/s¥yy and b/a(close to circularhole)usually decreases b0 of KI and benefits to the layout of stop-holes.The surface compressive stresses applied onto elliptical-hole(n)and crack(p)have significant influence on interacting SIFs but almost no on b0.Increasing n and p usually results in increase of interacting SIFs and facilitates crack propagation and fracture networks.The elliptical-hole orientation angle(a)and holed-cracked distance(t)have great influence on the interacting SIFs while have little effect on b0.The present method is not only simple(without any singular parts),high-accurate(due to exact fundamental stress solutions)and wider applicable(under far-field stresses and arbitrarily distributed surface stress)than the common methods,but also has the potential for the anisotropic problem involving multiple holes and cracks.
基金supported by the German Science Foundation(DFG, YU119/5-2)
文摘In the present work the fatigue crack growth in AISI304 specimens is investigated experimentally. In 3D finite element analysis the virtual crack closure technique is applied to calculate distributions and variations of the stress intensity factor along the surface crack front. It is confirmed that the stress intensity factor along the surface crack front varies non-uniformly with crack growth. Crack growth rate is proportional to the stress intensity factor distribution in the 3D cracked specimen. The fatigue crack growth in surface cracked specimens can be described by the Forman model identified in conventional compact tension specimens. For crack growth in the free specimen surface the arc length seems more suitable to quantify crack progress. Geometry and loading configuration of the surface cracked specimen seem to not affect the fatigue crack growth substantially.
文摘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.
基金Project(50925829) supported by the National Science Fund for Distinguished Young Scholars of ChinaProject(50908148) supported by the National Natural Science Foundation of ChinaProjects(2009-K4-23,2010-11-33) supported by the Research of Ministry of Housing and Urban Rural Development of China
文摘In order to perfectly reflect the dynamic corrosion of reinforced concrete (RC) cover in practical engineering,an analytic model of non-uniform corrosion induced cracking was presented based on the elastic-plastic fracture mechanics theory.Comparisons with the published experimental data show that the predictions given by the present model are in good agreement with the results both for natural exposed experiments and short-time indoor tests (the best difference is about 2.7%).Also it obviously provides much better precision than those models under the assumption of uniform corrosion (the maximal improved precision is about 48%).Therefore,it is pointed out that the so-called uniform corrosion models to describe the cover cracking of RC should be adopted cautiously.Finally,the influences of thickness of local rusty layer around the reinforcing steel bar on the critical corrosion-induced crack indexes were investigated.It is found that the thickness of local rusty layer has great effect on the critical mass loss of reinforcing steel,threshold expansion pressure,and time to cover cracking.For local rusty layer thickness with a size of a=0.5 mm,the time to cover cracking will increase by about one times when a/b (a,semi-minor axis;b,semi-major axis) changes from 0.1 to 1 mm.
基金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.
文摘This paper discusses stress intensity factor (SIF) calculations for surface cracks in round bars subjected to combined torsion and bending loadings. Different crack aspect ratios, a/b, ranging from 0.0 to 1.2 and relative crack depths, a/D, ranging from 0.1 to 0.6 were considered. Since the loading was non-symmetrical for torsion loadings, a whole finite element model was constructed. Then, the individual and combined bending and torsion loadings were remotely applied to the model. The equivalent SIF method, F*EQ, was then used explicitly to combine the individual SIFs from the bending and torsion loadings. A comparison was then carried out with the combined SIF, F*FE, obtained using the finite element analysis (FEA) under similar loadings. It was found that the equivalent SIF method successfully predicted the combined SIF for Mode I. However, discrepancies between the results determined from the different approaches occurred when Fm was involved. It was also noted that the predicted F* EE using FEA was higher than the F*EQ predicted through the equivalent SIF method due to the difference in crack face interactions.