Resolved shear stress intensity coefficient and fatigue crack growth in large crystals

Fatigue crack growth is caused primarily by shear decohesion due to dislocation motion in the crack tip region. The resolved shear stress, which drives dislocation in a crystal, is strongly orientation dependent, and therefore, the cyclic plastic deformation of the shear decohesion process is highly anisotropic.The crack planes are often inclined to the loading axis both in the inplane orientation and in the thickness direction. This inclination induces all three modes of the crack tip stress field, K_I, K_II, and K_III.Fatigue crack growth in large-grain Al 7029 aluminum alloy was studied. The crack tip stress fields of the test specimens are calculated with the finite element method. The values of K_I, K_II, and K_III are evaluated. The orientation of the crystal at a crack tip was determined with the Laue X-ray method. The crystal orientation and the calculated crack tip stress fields are used to compute the resolved shear stress intensity of each of the twelve slip systems of the crystal at the crack tip. The resolved shear stress field of a slip system is linearly proportional to the resolved shear stress intensity coefficient, RSSIC.The values of RSSIC thus evaluated are used to analyze the orientations of the crack plane and to correlate with the shear fatigue crack growth rate.     

Plane-strain crack tip field and a dual-parameter fracture-criterion for non-linear fracture mechanics

Hancock and Cowling measured the critical crack tip opening displacements, δ_f, at fracture initiation in HY-80 steel specimens of six different configurations. δ_f varied from 90 μm in a deeply double-edge-cracked tensile panel to 900 μm in a single-edge-cracked tensile panel.McMeeking and Parks, and Shih and German have shown by their finite element calculations that the characteristics of the plane strain crack tip fields in both large scale yielding and general yielding are strongly dependent on specimen geometry and load level.In this study, the plane strain crack tip fields in the specimens tested by Hancock and Cowling were calculated using the finite element method. The crack tip triaxial tensile stress field is strongly affected by specimen geometric constraint, and the state of the triaxial tensile stress in a crack tip region is monitored by the ratio between the local tensile stress and the effective stress, i.e., ( ), at a distance x=2δ from the crack tip. The values of ( ) vary from 3.1 for the double-edge-cracked tensile panel to 1.7 for the single-edge-cracked tensile panel. The δ_f measured by Hancock and Cowling correlates very well with the ratio ( ). δ_f is a measure of the fracture ductility of the material ahead of the crack tip, and the ductility decreases with an increase in the triaxial tensile stress, i.e., the ratio ( ).     

In Plane Displacement Formulation for Finite Cracked Plates Under Mode I Using Grid Method and Finite Element Analysis

A grid method is used to experimentally determine the in-plane displacement fields around a crack tip in a Single-Edge-Notch (SEN) tensile polyurethane specimen. Horizontal displacement u _x-exp and vertical displacement u _y-exp are expressed as functions of circular coordinates centred on the crack tip. These are compared with the approximate solutions of linear elastic fracture mechanics with a view to studying the applicability to polymers. The results show that this solution is not in agreement with the experiments at the focused on the vicinity of a crack tip. Taking this into account, an FEA program is developed with CAST3M for the purpose of comparing the experimental displacements and the numerical data. New formulations of displacements u _x and u _y are then developed. These formulations are derived from the principle of superposition and based on Arakawa’s formulation. With the displacement gradients obtained from the FEA and the new formulations, the determination of J-integrals is found to be in very good agreement with those derived from numerical calculation. Consequently, the proposed formulations can give displacement fields compatible with the J-integral calculation for the region near the crack tip. An application based on an experimental test is proposed to evaluate the performances of the proposed formulations.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

Stress field of a damaged material near the crack tip

In this paper, a deformation theory of plasticity for damaged materials is proposed. An asymptotic expression forH near a crack tip is obtained. Finally, the stress and strain fields near the crack tip are presented.     

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for theT stresses ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analysed in detail. It is found that normal stresses {ie27-1} and {ie27-2} ahead of the crack tip, are characterised byQ fields if the crack is within a stiff material and the parameters |p _T | and |q _T | are very small, whereQ is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses {ie27-3} and {ie27-4} are dominated by theQ field plusT stress. This work was supported by the Swedish Research Council for Engineering Sciences.     

Anisotropic plastic stress fields at a slowly propagating crack tip

Under the condition that any perfectly plastic stress components at a crack tip are nothing but the functions of 0 only making use of equilibrium equations. Hill anisotropic yield condition and unloading stress-strain relations, in this paper, we derive the general analytical expressions of anisotropic plastic stress fields at the slowly steady propagating tips of plane and anti-plane strain. Applying these general analytical expressions to the concrete cracks, the analytical expressions of anisotropic plastic stress fields at the-slowly steady propagating tips of Mode I and Mode III cracks are obtained. For the isotropic plastic material, the anisotropic plastic stress fields at a slowly propagating crack tip become the perfectly plastic stress fields.     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

The stress-strain fields at crack tip in axially cracked cylindrical shells and the calculation of stress intensity factors

A perturbation solution for stress-strain fields (including modes I, II, III) at crack tip in axially cracked cylindrical shells is given. The analysis, using 10th-order differential equations which take the transverse shear deformations into account, involves perturbation in a curvature parameter λ², (λ²=[12(1-v ²)]^1/2 a ²/Rh). Stress intensity factors for finite size cylindrical shells under bending and internal pressure loading are evaluated. A good accuracy can be obtained without using fine meshes in a region near the crack tip. Besides, the influence of the finite size and the shearing stiffness on bulging factors, which are commonly used in engineering, are analyzed.     

Anti-plane shear of kinked interface crack in bonded dissimilar anisotropic solid

Complex potentials are derived to describe the anti-plane singular shear stress fields around a kinked crack, the main portion of which is embedded along the interface of two dissimilar anisotropic elastic media. This is accomplished by formulating the problem as singular integral equations with generalized Cauchy kernels. The shear stress singularity at the kink differs from the familiar inverse square root of the local distance; it is found to influence the magnitude of the Mode III crack tip stress intensity factor, K₃. Numerical results of K₃ are obtained and displayed in graphical forms for different degree of material anisotropy and crack dimensions.     

Experimental GH singularity field around growing crack-tip

The displacement fieldsu _x ,u _y at growing crack tip of LY12-M specimens with double edge cracks are measured using moire method. The experimental singularity fields are compared with GH theoretical field [12–14]. The size and shape of the experimental GH singularity fields are obtained. The error in both the experimental and theoretical evaluations is controlled within ±10%. The experiments show that there is singularity dominant around a growing crack tip. The shape of this dominant region ranges from butterfly wing to oblate and circular. Inside GH-field, there is a 3-D deformed damage zone where no GH singularity exists. The project suppoted by National Natural Science Foundation of China     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

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     

含裂纹体的单位分解扩展无网格方法

不连续体的数值模拟尤其是动态裂纹的追踪问题一直是工程界研究的热点和难点问题。无网格方法仅仅需要结点信息,非常适合于求解这类问题。基于单位分解思想,在移动最小二乘近似函数(MLS)中根据裂纹面的不连续位移增加一个Heaviside函数,在裂尖则增加四个扩展函数描述渐进裂纹位移场;应用Galerkin方法推导了平衡方程的离散线性方程,并给出了求解裂纹问题应力强度因子的计算公式。与其他类型的扩展无网格相比,在裂尖处近似函数不需要使用可视准则,很容易生成r1/2奇异;另一个优势是影响域并没有因为裂纹的存在而改变,不会降低方程的稀疏性,求解效率较高。数值算例表明,该方法能方便有效地模拟不连续问题,具有十分广阔的应用空间。     

Dynamic fracture analysis by moiré interferometry

Dynamic moiré interferometry was used to measure separately theu- andv-displacement fields surrouding a rapidly propagating crack tip in Homalite-100 and 7075-T6 aluminum-alloy plates. These transient crack-tip-displacement data were then used to compute the dynamic stress-intensity factor and the remote stress component.J-integral values were also estimated using the static approximate procedure of Kanget al. This static analysis provided the correctJ when the contour integral was taken within 3 mm of the crack tip.Paper was presented at the 1991 SEM Spring Conference on Experimental Mechanics held in Milwaukee, WI on June 9–13.     

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.     

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacementD ₂ induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with realH, the normal electric displacementD ₂ is constant along the crack faces and electric fieldE ₂ has the singularity ahead of the crack tip and a jump across the interface. The project is supported by the National Natural Science Foundation of China(No. 19704100) and the Natural Science Foundation of Chinese Academy of Sciences(No. KJ951-1-201).     

High-order asymptotic analysis for the crack in nonlinear material

Accurate high-order asymptotic analyses were carried out for Mode II plane strain crack in power hardening materials. The second-order crack tip fields have been obtained. It is found that the amplitude coefficientk ₂ of the second term of the asymptotic field is correlated to the first order field as the hardening exponentn<n ^* (n ^*≈5), but asn≥n ^*,k ₂ turns to become an independent parameter. Our results also indicated that, the second term of the asymptotic field has little influence on the near-crack-tip field and can be neglected whenn<n ^*. In fact,k ₂ directly reflects the effects of triaxiality near the crack tip, the crack geometry and the loading mode, so that besidesJ-integral it can be used as another characteristic parameter in the two-parameter criterion. The project supported by National Natural Science Foundation of China     

Exponential stress and strain singularity at tip of mode I growing crack in power hardening plastic material

The stress-strain distribution near the tip of a Mode I growing crack in a power hardening plastic material is reconsidered. Two types of asymptotic equations are derived and solved numerically. It is shown that when the crack tip is approached, the stress is singular of the order r^−δ, while the strain is singular of the order r^−nδ, where r is the distance measured from the crack tip. The parameter δ is a constant; it depends on the hardening exponent n being greater than one.