A finite element study of the near crack tip deformation of a ductile material under mixed mode loading

In ductile fracture, voids near a crack tip play an important role. From this point of view, a large deformation finite element analysis has been made to study the deformation, stress and strain, and void ratio near the crack tip under mixed mode plane strain loading conditions, employing Gurson's constitutive equation which has taken into account the effects of void nucleation and growth. The results show that: (i) one corner of the crack tip sharpens while the other corner blunts, (ii) the stress and strain distributions except for the near crack tip region, can be superimposed by normalizing distance from the crack tip by a crack tip deformation length, i.e., a steady-state solution under a mixed mode condition has been obtained, (iii) the field near a crack tip can be divided into four characteristic fields (K field, HRR field, blunted crack tip field, and damaged region), and (iv) the strain and void volume fraction become concentrated in the sharpened part of a crack tip with increasing Mode II component.     

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     

Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic–viscoplastic constitutive relation with various hardening–softening–hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip.On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening(or softening–hardening) has a particularly strong effect on the near crack tip stress and deformation fields.     

复合型韧性断裂实验和控制参数

通过对不同韧性材料在各种平面复合载荷形式下裂纹启裂阶段裂端变形、启裂位置和扩展方向的系统的宏微观实验验证及计算分析,考察了韧性断裂参数空穴扩张比的分布特征和裂纹启裂及扩展方向的关系．得到：对于不同韧性的材料,在裂端的钝化变形区域,空穴扩张比的极大值区对应干裂纹的启裂位置,裂纹启裂时钝化裂端前缘空穴扩张比的临界值不敏感于复合比的变化．而对于裂纹启裂后的扩展方向,则需根据具体材料在相应的特定区域中比较空穴扩张比参数极大值的分布特征,需经进一步的分析,从而确定裂纹的扩展方向．实验及计算结果表明,尽管复合型断裂时裂纹启裂及扩展的机理极其复杂,用于韧性材料复合型断裂的空穴扩张比参数仍能很好地预测裂端的启裂及扩展方向,可作为复合型韧性断裂过程的控制参数．     

Cracking characteristics of mixed mode dislocations near a lip-like mode crack

This work is concerned with the cracking characteristics of mixed mode dislocations near a lip-like mode crack, stress intensity and strain energy density factor are obtained by using conformal mapping, singularity analysis and Cauchy integrals. Shielding effect generated by screw dislocation near a lip-like mode crack decreases with the increment of the distance between screw dislocation and crack tip. Larger distance between two faces of the crack leads to the shielding effect waning. The strain energy density factor of mode III decreases with the increment of the distance between dislocation and crack tip. Larger distance between two faces of lip-like mode crack also leads to the strain energy density factor waning and encourages crack initiation; the shielding effects generated by edge dislocation near the crack decrease with the increment of the distance between edge dislocation and crack tip.     

ASYMPTOTIC ANALYSIS OF MODE Ⅱ STATIONARY GROWTH CRACK ON ELASTIC-ELASTIC POWER LAW CREEPING BIMATERIAL INTERFACE

A mechanical model was established for mode Ⅱ interfacial crack static growing along an elastic-elastic power law creeping bimaterial interface. For two kinds of boundary conditions on crack faces, traction free and frictional contact, asymptotic solutions of the stress and strain near tip-crack were given. Results derived indicate that the stress and strain have the same singularity, there is not the oscillatory singularity in the field; the creep power-hardening index n and the ratio of Young' s module notably influence the cracktip field in region of elastic power law creeping material and n only influences distribution of stresses and strains in region of elastic material. When n is bigger, the creeping deformation is dominant and stress fields become steady, which does not change with n.Poisson ' s ratio does not affect the distributing of the crack- tip field.     

Crack emanating from an open notch

This paper is concerned with the plane elastostatic problem of a crack emanating from an open notch where it is assumed that the crack is directed parallel to one side of the notch. The problem is reduced to two simultaneous Fredholm integral equations which are then solved numerically for three loading conditions of the crack faces. In particular, it is shown that the crack tip intensity factors are almost independent of crack length if the singular stress field associated with the uncracked notch is used as the loading of the crack faces.     

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005     

A transverse crack tip field in bimaterial

The asymptotic field near an interface crack tip is analyzed with the fully nonlinear theory. By dividing the crack tip field into narrowing sectors and an expanding sector, the asymptotic equations for the crack tip field are derived and solved. The singular characters of stress and strain near the crack tip are revealed.     

A MIXED ELECTRIC BOUNDARY VALUE PROBLEM FOR AN ANTI－PLANE PIEZOELECTRIC CRACK

The analytical continuation method is adopted to solve a mixed electric boundary value problem for a piezoelectric medium under anti-plane deformation. The crack face is partly conductive and partly impermeable. The results show that the stress intensity factor is identical with the mode III stress intensity factor independent of the conducting length. But the electric field and the electric displacement are dependent on the electric boundary conditions on the crack faces and are singular not only at the crack tips but also at the junctures between the impermeable part and conducting portions.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

AN INVESTIGATION ON ASYMPTOTIC NEAR-TIP FIELDS OF DYNAMIC CRACK IN AN ELASTIC-PERFECTLY PLASTIC COMPRESSIBLE MATERIAL

The stress and deformation fields near the tip of a mode-I dynamic crack steadilypropagating in an elastic-perfectly plastic compressible material are considered under plane strain condi-tions. Within the framework of infinitesimal displacement gradient theory, the material is character-ized by the Von Mises yield criterion and the associated J_2 flow theory of plasticity. Through rigorousmathematical analysis, this paper eliminates the possibilities of elastic unloading and continuousasymptotic fields with singular deformation, and then constructs a fully continuous and boundedasymptotic stress and strain field. It is found that in this solution there exists a parameter (?)_0 whichcannot be determined by asymptotic analysis but may characterize the effect of the far field. Lastly thevariations of continuous stresses, velocities and strains around the crack tip are given numerically fordifferent values of (?)_0.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

A mixed electric boundary value problem for a two-dimensional piezoelectric crack

In this paper, a mixed electric boundary value problem for a two-dimensional piezoelectric crack problem is presented, in the sense that the crack face is partly conducting and partly impermeable. By the analytical continuation method, the unknown electric charge distributions on the upper and lower conducting crack faces are reduced to two decoupled singular integral equations and then these two equations are converted into algebraic equations to find the full field solution. Though the results suggest that the stress intensity factors at the crack tip are identical to those of conventional piezoelectric materials, but the electric field and electric displacement are related to the electric boundary conditions on the crack faces. The electric field and electric displacement are singular not only at crack tips but also at the junctures between the impermeable part and conducting parts. Numerical results for the variations of the electric field, electric displacement field and J-integral with respect to the normalized impermeable crack length are shown. Some discussions for the energy release rate and the J-integral are made.     

复合型裂纹应力场的研究及其应力强度因子的全息光弹性测定

本文应用复变函数解法,等出复合型中心裂纹板弹性应力场的精确解及主应力和与主应力差的精确表达式。通过与各自的奇异表达式比较,得到了主应力和与主应力差的远近场关系图谱。利用这些图谱以及全息光弹性试验所获得的远场等和线与等差线条纹,就能确定裂纹尖端的应力强度因子 K_Ⅰ,K_Ⅱ。实例表明:本法概念清晰、演算简便、精度较高。     

Anti-plane shear crack normal to and terminating at the interface of two bonded piezoelectric ceramics

The electroelastic analysis of two bonded dissimilar piezoelectric ceramics with a crack perpendicular to and terminating at the interface is made. By using Fourier integral transform, the associated boundary value problem is reduced to a singular integral equation with generalized Cauchy kernel, the solution of which is given in closed form. Results are presented for a permeable crack under anti-plane shear loading and in-plane electric loading. Obtained results indicate that the electroelastic field near the crack tip in the homogeneous piezoelectric ceramic is dominated by a traditional inverse square-root singularity, while the electroelastic field near the crack tip at the interface exhibits the singularity of power law r^−α, r being distance from the interface crack tip and α depending on the material constants of a bi-piezoceramic. In particular, electric field has no singularity at the crack tip in a homogeneous solid, whereas it is singular around the interface crack tip. Numerical results are given graphically to show the effects of the material properties on the singularity order and field intensity factors.     

Finite element analysis of void growth near a blunting crack tip

Large deformation finite element analysis has been used to study the near crack tip growth of long cylindrical holes aligned parallel to the plane of a mode I plane strain crack. The near crack tip stress and deformation fields are analyzed. The results show that the holes are pulled towards the crack tip and change their shape to approximately elliptical with the major axis radial to the crack. They also grow faster directly ahead of the crack than at an angle to the crack plane. Several crack-hole coalescence criteria are discussed and estimates for the conditions for fracture initiation are given and compared with experimental results. The range of estimates now available from finite element calculations coincides quite well with the range of experimental data for materials containing inclusions which are only loosely bonded to the matrix.     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

新型各向异性核工业石墨断裂力学行为的实验研究和数值分析

本文分别采用激光和白光DSCM(数字散斑相关测量)方法对一种新型各向异性核工业石墨的裂纹尖端位移、应变场进行了实验研究，两种方法都取得了较好的结果。考虑核工业石墨制备过程形成的各向异性特点，本文构建了三维各向异性有限元模型，采用奇异单元，计算模拟石墨裂纹尖端的变形和应力场，通过对实验结果和有限元计算结果的比较可以发现两者具有相近的趋势。     

Perturbed magnetic field of an infinite plate with a centered crack

Deforming a cracked magnetoelastic body in a magnetic field induces a perturbed magnetic field around the crack. The quantitative relationship between this perturbed field and the stress around the crack is crucial in developing a new generation of magnetism-based nondestructive testing technologies. In this paper, an analytical expression of the perturbed magnetic field induced by structural deforma- tion of an infinite ferromagnetic elastic plate containing a centered crack in a weak external magnetic field is obtained by using the linearized magnetoelastic theory and Fourier transform methods. The main finding is that the perturbed magnetic field intensity is proportional to the applied tensile stress, and is dominated by the displacement gradient on the boundary of the magnetoelastic solid. The tangential component of the perturbed magnetic-field intensity near the crack exhibits an antisymmetric distribution along the crack that reverses its direction sharply across its two faces, while the normal component shows a symmetric distribution along the crack with singular points at the crack tips.