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.展开更多
A new exact and universal conformal mapping is proposed. Using Muskhelishvili＇s complex potential method, the plane elasticity problem of power function curved cracks is investigated with an arbitrary power of a natu...A new exact and universal conformal mapping is proposed. Using Muskhelishvili＇s complex potential method, the plane elasticity problem of power function curved cracks is investigated with an arbitrary power of a natural number, and the general solutions of the stress intensity factors （SIFs） for mode I and mode II at the crack tip are obtained under the remotely uniform tensile loads. The present results can be reduced to the well-known solutions when the power of the function takes different natural numbers. Numerical examples are conducted to reveal the effects of the coefficient, the power, and the projected length along the x-axis of the power function curved crack on the SIFs for mode I and mode II.展开更多
An asymptotic algorithm is applied to the problem of a finite, thermo-elastic solid containing a surface breaking crack, when the exterior surface is subjected to oscillatory thermal loading. This algorithm involves t...An asymptotic algorithm is applied to the problem of a finite, thermo-elastic solid containing a surface breaking crack, when the exterior surface is subjected to oscillatory thermal loading. This algorithm involves the study of a model problem. An analytical and numerical study of this model problem of a thermo-elastic half space containing a surface breaking crack and subjected to oscillatory thermal loading is presented. The crack surface is traction free. In particular, the amplitude of the stress intensity factor at the crack vertex is found as a function of the crack depth and the frequency of thermal oscillation.展开更多
Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is ...Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.展开更多
In this paper, a finite crack with constant length （Yoffe type crack） propagating in a functionally graded coating with spatially varying elastic properties bonded to a homogeneous substrate of finite thickness unde...In this paper, a finite crack with constant length （Yoffe type crack） propagating in a functionally graded coating with spatially varying elastic properties bonded to a homogeneous substrate of finite thickness under anti-plane loading was studied. A multi-layered model is employed to model arbitrary variations of material properties based on two linearly-distributed material compliance parameters. The mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. Some numerical examples are given to demonstrate the accuracy, efficiency and versatility of the model. The numerical results show that the graded parameters, the thicknesses of the interfacial layer and the two homogeneous layers, the crack size and speed have significant effects on the dynamic fracture behavior.展开更多
基金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(Nos.10932001,11072015, and 10761005)the Scientific Research Key Program of Beijing Municipal Commission of Education (No.KZ201010005003)+1 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China(No.20101102110016)the Ph.D.Innovation Foundation of Beijing University of Aeronautics and Astronautics(No.300351)
文摘A new exact and universal conformal mapping is proposed. Using Muskhelishvili＇s complex potential method, the plane elasticity problem of power function curved cracks is investigated with an arbitrary power of a natural number, and the general solutions of the stress intensity factors （SIFs） for mode I and mode II at the crack tip are obtained under the remotely uniform tensile loads. The present results can be reduced to the well-known solutions when the power of the function takes different natural numbers. Numerical examples are conducted to reveal the effects of the coefficient, the power, and the projected length along the x-axis of the power function curved crack on the SIFs for mode I and mode II.
文摘An asymptotic algorithm is applied to the problem of a finite, thermo-elastic solid containing a surface breaking crack, when the exterior surface is subjected to oscillatory thermal loading. This algorithm involves the study of a model problem. An analytical and numerical study of this model problem of a thermo-elastic half space containing a surface breaking crack and subjected to oscillatory thermal loading is presented. The crack surface is traction free. In particular, the amplitude of the stress intensity factor at the crack vertex is found as a function of the crack depth and the frequency of thermal oscillation.
基金supported by the China Aviation Industry Corporation I Program (ATPD-1104-02).
文摘Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.
基金Project supported by the National Natural Science Foundation of China (Nos. 10802078 and 10872150)China Postdoctoral Science Foundation (No. 20100471006)
文摘In this paper, a finite crack with constant length （Yoffe type crack） propagating in a functionally graded coating with spatially varying elastic properties bonded to a homogeneous substrate of finite thickness under anti-plane loading was studied. A multi-layered model is employed to model arbitrary variations of material properties based on two linearly-distributed material compliance parameters. The mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. Some numerical examples are given to demonstrate the accuracy, efficiency and versatility of the model. The numerical results show that the graded parameters, the thicknesses of the interfacial layer and the two homogeneous layers, the crack size and speed have significant effects on the dynamic fracture behavior.
