Crack growth under plane stress conditions in an elastic perfectly-plastic material

Mode-I crack growth under conditions of generalized plane stress has been investigated. It has been assumed that near the plane of the crack in the loading zone, the simple stress components corresponding to a central fan field maintain validity up to the elastic-plastic boundary. By the use of expansions of the particle velocities in the coordinate y, and by matching of the relevant stress components and particle velocities to the dominant terms of appropriate elastic fields at the elastic-plastic boundary, a complete solution has been obtained for ε_y in the plane of the crack. The solution applies from the propagating crack tip up to the moving elastic-plastic boundary. The strain fields for a self-similar crack nucleating at a point and for steady-state propagation of a crack have been considered as special cases.     

A crack of arbitrary shape inside an infinite transversely isotropic cylinder under arbitrary load

Summary  In the first part of the article an infinite circular cylinder is considered, made of transversely isotropic elastic material and weakened by a plane crack perpendicular to its axis O z. The crack is opened by an arbitrary normal stress. The second part is devoted to the same crack loaded by an arbitrary tangential stress. The complete solution in both cases is presented as a sum of the solution of a similar problem of a crack in an infinite space and an integral transform term, the parameters of which are determined from a set of linear algebraic equations derived from the boundary conditions. Governing integral equations with respect to the yet unknown crack displacement discontinuities are obtained. In the case of a circular crack, these equations can be inverted and solved by the method of consecutive interations. Received 30 November 2000; accepted for publication 3 May 2001     

The stress field near the front of an arbitrarily shaped crack in a three-dimensional elastic body

The problem stated in the title is investigated with special emphasis on the first three terms of the stress expansion, proportional to r ^-1/2, r ⁰=1 and r ^1/2 respectively, where r denotes the distance to the crack front. The particular case of a plane crack with a straight front and of stresses independent of the distance along the latter is studied first. It is shown that the classical plane strain and antiplane solutions must be supplemented by a few additional particular solutions to obtain the full stress expansion. The general case is then considered. The stress expansion is studied by writing the field equations (equilibrium, strain compatibility and boundary conditions) in a system of suitable curvilinear coordinates. It is shown that the number of independent constants in the stress expansion is the same as in the particular case considered previously but that the curvatures of the crack and its front and the non-uniformity of the stresses along the latter induce the appearance of corrective terms in this expansion.     

A generalized force measure of conditions at crack tips

A tentative measure of the forces tending to cause crack growth—the apparent crack extension force—is defined within the framework of continuum mechanics. By an associated fracture criterion initiation of growth may be predicted as well as the direction of preferred growth. The theory is specialized to elastic, viscoelastic and elastic-plastic materials. Under specified conditions the apparent crack extension force may be expressed by surface integrals over the boundary of an arbitrary part of the body for quasi-static deformation and for steady-state propagation of the crack. For plane deformation and for infinitesimal deformation under plane stress conditions these surface integrals reduce to path independent line integrals which include the J integral by Rice[1] and the G integral by Sih[2] as special cases.     

A dynamic model of bridging fiber pull-out of composite materials

An elastic analysis of an internal central crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane is carried out. In this paper a dynamic model of bridging fiber pull-out is presented for analyzing the distributions stress and displacement of composite materials with the internal central crack under the loading conditions of an applied non-uniform stress and the traction forces on crack faces yielded by the fiber pull-out model. Thus the fiber failure is determined by maximum tensile stress, the fiber breaks and hence the crack propagation should occur in self-similar fashion. By reducing the dynamic model to the Keldysh–Sedov mixed boundary value problem, a straightforward and easy analytical solution can be attained. When the crack extends, its fibers continue to break. Analytical study on the crack extension under the action of an inhomogeneous point force Px/t, Pt is obtained for orthotropic anisotropic body, respectively; and it can be utilized to attain the concrete solutions of the model by the ways of superposition.     

A Mode- II Griffith Crack in Decagonal Quasicrystals

By using the method of stress functions, the problem of mode-II Griffith crack in decagonal quasicrystals was solved. First, the crack problem of two-dimensional quasicrystals was decomposed into a plane strain state problem superposed on anti-plane state problem and secondly, by introducing stress functions, the18 basic elasticity equations on coupling phonon-phason field of decagonal quasicrystals were reduced to a single higher-order partial differential equations. The solution of this equation under mixed boundary conditions of mode-II Griffith crack was obtained in terms of Fourier transform and dual integral equations methods. All components of stresses and displacements can be expressed by elemental functions and the stress intensity factor and the strain energy release rate were determined. Biography: GUO Yu-cui (1962-), Associate professor, Doctor     

Crack tip field andJ-integral with strain gradient effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventionalJ ₂ deformation theory. No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory. Two typical crack problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-IK-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory. The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, theJ-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases. The project supported by the National Natural Science Foundation of China (19704100 and 10202023) and the Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20)     

Stress intensity factors for asymmetric branched cracks in plane extension by using crack-tip displacement discontinuity elements

Stress intensity factors are important in the analysis of cracked materials. They are directly related to the fracture propagation and fatigue crack growth criteria. Based on the analytical solution (Crouch, S.L., 1976. Solution of plane elasticity problems by displacement discontinuity method, Int. J. Numer. Methods Eng. 10, pp. 301–343; Crouch, S.L., Starfield, A.M., 1983. Boundary Element Method in Solid Mechanics, with Application in Rock Mechanics and Geological Mechanics, London, Geore Allon and Unwin, Bonton, Sydney) to the problem of a constant discontinuity in displacement over a finite line segment in the x, y plane of an infinite elastic solid, recently, the crack-tip displacement discontinuity element which can be classified as the left and right crack-tip displacement discontinuity elements are developed by the author Yan, X., (in press. A special crack-tip displacement discontinuity element, Mechanics Research Communications) to model the crack-tip fields to more accurately compute the stress intensity factors of cracks in general plane elasticity. In the boundary element implementation the left or the right crack-tip displacement discontinuity element is placed locally at the corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundaries. To prove further the efficiency of the suggested approach and provide more results of the stress intensity factors, in this study, analysis of an asymmetric branched crack bifurcated from a main crack in plane extension is carried out.     

THE PROBLEM OF I-II COMBINED PLANE CRACK SOLVED WITH WESTERGAARD STRESS FUNCTION

直接获得I-II复合型平面裂纹问题裂纹尖端区域的应力场是一个比较复杂的问题,在此应用Westergaard应力函数求解I-II复合型平面裂纹问题,导出了裂纹尖端区域应力分量的表达式。该方法推导过程简单,物理概念清晰,其结果与一般断裂力学教材和文献中的结果一致。同时,应用叠加原理将裂纹面上的作用力转化为裂纹外边界的受力,给出了解决裂纹面上有作用力的I-II复合型平面裂纹问题的解题方法。     

考虑摩擦作用闭合裂纹应力强度因子计算的边界元分区算法

本文采用边界元分区算法研究考虑摩擦的闭合裂纹问题,简单裂纹系问题及非均匀介质中裂纹问题.在裂纹尖端采用了1/4面力奇异单元,并对相应的奇异积分给出了数值处理.对于复杂载荷下裂纹面计算给出了增量迭代算法,并采用方程减缩技术使迭代仅在裂纹上进行.计算实例表明方法是可行的.     

Elastoplastic near field solution for mode II crack loaded by remote shear stress in an infinite plate

A suitable elastic stress field near the crack line which satisfied the far field boundary conditions and the boundary conditions of the crack surfaces has been obtained and successful analysis has been made of a near crack line field for an infinite elastic-perfectly plastic medium containing a quasi-statically propagating plane stress crack subjected to far field shear stress. It is shown that the solutions of the problem of mode II crack loaded by remote shear stress from the Westergaard method in some previous papers is used as the elastic stress field near the crack line, are inappropriate.     

Influence of non-singular stress terms and specimen geometry on small-scale yielding at crack tips in elastic-plastic materials

The plane strain elastic-plastic state at a crack tip is determined for compact tension, bend, double edge-cracked and centre-cracked specimens using a finite element method with triangular constant-strain elements. The solutions are found to differ by 10 to 30 per cent at the ASTM-limit as regards fracture surface displacement, normal stress and plastic zone size. In order to bring the boundary layer solution for the crack problem into agreement with the solution for a specific specimen one has to modify this solution. The modification consists of an addition to the boundary tractions for the boundary layer problem of tractions corresponding to the non-singular, constant second term in a series expansion of the normal stress parallel to the crack plane.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Anti-plane dynamic hole–crack interaction in a functionally graded piezoelectric media

The anti-plane dynamic problem of a functionally graded piezoelectric plane containing a hole–crack system is treated by a non-hypersingular traction-based boundary integral equation method. The material parameters vary exponentially in the same manner in an arbitrary direction. The system is loaded by an incident SH-type wave, and impermeable boundary conditions are assumed. Using a frequency-dependent fundamental solution of the wave equation, the boundary value problem is transformed into a system of integro-differential equations along the boundary of the hole and on the crack line. Its numerical solution yields the dynamic stress intensity factors and stress concentration factors. A parametric study reveals their dependence on the hole–crack scenario and its geometry, characteristics of the dynamic load and magnitude and direction of material inhomogeneity.     

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

弹性半平面中边界裂纹研究的新方法

讨论了载荷作用在裂纹面上的弹性半平面边界裂纹问题.研究以线弹性断裂力学为基础,采用复变函数方法以及Riemann-Hilbert(R-H)边值问题的一般理论,将问题分拆为含有限裂纹的全平面问题与无裂纹的半平面问题的叠加,计算得到裂纹尖端的应力强度因子.与文献结果比较,该方法具有精度高的优点.     

Application of complex path-independent integrals to the solution of the problem of a straight crack in a finite plane isotropic elastic medium

The classical problem of a straight crack in a finite, plane, isotropic, elastic medium of arbitrary shape is reconsidered by the well-known method of Muskhelishvili for such a crack (but in an infinite medium). Both the crack and the boundary of the medium are assumed loaded in an arbitrary way. It is shown that this problem can be completely solved if the numerical values of the first complex potential (z) of Muskhelishvili are known along a closed contour surrounding the crack, probably along the boundary of the medium. To this end, complex path-independent integrals associated with (z) and Chebyshev polynomials have been used. Numerical results for the stress intensity factors are displayed in an application. Generalizations of the method are also proposed and the second fundamental crack problem, the problem of a crack in an anisotropic medium and the problem of an interface crack between two isotropic media are considered in some detail.     

Axisymmetric displacement boundary value problem for a penny-shaped crack

Summary A solution is derived from equations of equilibrium in an infinite isotropic elastic solid containing a penny-shaped crack where displacements are given. Abel transforms of the second kind stress and displacement components at an arbitrary point of the solid are known in the literature in terms of jumps of stress and displacement components at a crack plane. Limiting values of these expressions at the crack plane together with the boundary conditions lead to Abel-type integral equations, which admit a closed form solution. Explicit expressions for stress and displacement components on the crack plane are obtained in terms of prescribed face displacements of crack surfaces. Some special cases of the crack surface shape functions have been given in the paper.     

Surface contact of elliptical crack under normally incident tension–compression wave

Surface contact interaction of a plane elliptical crack under normally incident tension–compression wave is solved by the method of boundary integral equations. The contact forces and the displacement discontinuity of the crack edges are examined. The dependence of the mode I stress intensity factor on the wave number is studied. The solution is compared with the results obtained for an elliptical crack without allowance for crack edges contact interaction.