A Three-dimensional Numerical Study of Mixed Mode (I and II) Crack Tip Fields in Elastic–plastic Solids

The stress fields near a crack front in a ductile solid are essentially three-dimensional (3D) in nature. The objective of this paper is to investigate the structure of these fields and to establish the validity of two-dimensional (2D) plane stress and plane strain approximations near the crack front under mixed mode (combined modes I and II) loading. To this end, detailed 3D and 2D small strain, elastic–plastic finite element simulations are carried out using a boundary layer (small scale yielding) formulation. The plastic zones and radial, angular and thickness variations of the stresses are studied corresponding to different levels of remote elastic mode mixity and applied load, as measured by the plastic zone size with respect to the plate thickness. The 3D results are compared with those obtained from 2D simulations and asymptotic solutions. It is found that, in general, plane stress conditions prevail at a distance from the crack front exceeding half the plate thickness, although it could be slightly smaller for mode II predominant loading. The implications of the 3D stress distribution on micro-void growth near the crack front are briefly discussed.     

J-Dominance in Mixed Mode Ductile Fracture Specimens

The asymptotic mixed mode crack tip fields in elastic-plastic solids are scaled by the J-integral and parameterized by a near-tip mixity parameter, M _{_p} . In this paper, the validity and range of dominance of these fields are investigated. To this end, small strain elastic-plastic finite element analyses of mixed mode fracture are first performed using a modified boundary layer formulation. Here, a two term expansion of the elastic crack tip field involving the stress intensity factor |K| the elastic mixity parameter M _{_e} as well as the T-stress is prescribed as remote boundary conditions. The analyses are conducted for different values of M _{_e} and the T-stress. Next, several commonly used mixed mode fracture specimens such as Compact Tension Shear (CTS), Four Point Bend (4PB), and modified Compact Tension specimen are considered. Here, the complete range of loading from contained yielding to large scale yielding is analyzed. Further, different crack to width ratios and strain hardening exponents are considered. The results obtained establish that the mixed mode asymptotic fields dominate over physically relevant length scales in the above geometries, except for predominantly mode I loading and under large scale yielding conditions. This revised version was published online in August 2006 with corrections to the Cover Date.     

Three‐dimensional mixed‐mode (I and II) crack‐front fields in ductile thin plates — effects of T‐stress

It has been well‐established that the non‐singular T‐stress provides a first‐order estimate of geometry and loading mode (e.g. tension versus bending) effects on elastic–plastic crack‐front field under mode I loading conditions. The objective of this paper is to exam the T‐stress effect on three‐dimensional (3D) crack‐front fields under mixed‐mode (modes I and II) loading. To this end, detailed 3D small strain, elastic–plastic simulations are carried out using a 3D boundary layer (small‐scale yielding) formulation. Characteristics of near crack‐front fields are investigated for a wide range of T‐stresses (T/σ₀ = ?0.8, ?0.4, 0.0, 0.4, 0.8). The plastic zones and thickness and angular and radial variations of the stresses are studied, corresponding to two values of the remote elastic mixity parameters M_e = 0.3 and 0.7, under both low and high levels of applied loads. It is found that different T‐stresses have a significant effect on the plastic zones size and shapes, regardless of the mode mixity and load level. The thickness, angular and radial distributions of stresses are also affected markedly by T‐stress. It is important to include these effects when investigating the mixed‐mode ductile fracture failure process in thin‐walled structural components.     

Two-parameter characterization of the near-tip stress fields for a bi-material elastic-plastic interface crack

A particular case of interface cracks is considered. The materials at each side of the interface are assumed to have different yield strength and plastic strain hardening exponent, while elastic properties are identical. The problem is considered to be a relevant idealization of a crack at the fusion line in a weldment. A systematic investigation of the mismatch effect in this bi-material plane strain mode I dominating interface crack has been performed by finite strain finite element analyses. Results for loading causing small scale yielding at the crack tip are described. It is concluded that the near-tip stress field in the forward sector can be separated, at least approximately, into two parts. The first part is characterized by the homogeneous small scale yielding field controlled by J for one of the interface materials, the reference material. The second part which influences the absolute value of stresses at the crack tip and measures the deviation of the fields from the first part can be characterized by a mismatch constraint parameter M. Results have indicated that the second part is a very weak function of distance from the crack tip in the forward sector, and the angular distribution of the second part is only a function of the plastic hardening property of the reference material.     

The dynamic near crack-line fields for mode II crack growth in an elastic perfectly-plastic solid

The dynamic near crack-line fields for mode II crack growth in an elastic perfectly-plastic solid are investigated under plane strain and plane stress conditions. In each case, by expanding the plastic fields and the governing equations in the coordinate y, the problem is reduced to solving a system of nonlinear ordinary differential equations which is similar to that of mode III derived by Achenbach and Z.L.Li. An approximate solution for small values of x is obtained and matched with the elastic field of a blunt crack at the elastic-plastic boundary. The crack growth criterion of critical strain is employed to determine the value of K _II of the far-field that would be required for a steadily growing crack.     

Characterization of crack tip stress fields in test specimens using mode mixity parameters

The aim of this study is to represent the combined effect of mode mixity, specimen geometry and relative crack length on the $T$ -stress, elastic–plastic stress fields, integration constant $I_{n}$ , angle of initial crack extension, and the plastic stress intensity factor. The analytical and numerical results are obtained for the complete range of mixed modes of loading between mode I and mode II. For comparison purposes, the reference fields for plane mixed-mode problems governing the asymptotic behavior of the stresses and strains at the crack tip are developed in a power law elastic–plastic material. For the common experimental fracture mechanics specimen geometries considered, the numerical constant of the plastic stress field $I_{n}$ and the $T$ -stress distributions are obtained as a function of the dimensionless crack length and mode mixity. A method is also suggested for calculating the plastic stress intensity factor for any mixed-mode I/II loading based on the $T$ -stress and power law solutions. It is further demonstrated that in both plane stress and the plane strain, the plastic stress intensity factor can be used to characterize the crack tip stress fields for a variety of specimen geometries and different mixed-mode loading. The applicability of the plastic stress intensity factor to analysis of the in-plane and out-of-plane constraint effect is also discussed.     

Precise elastic-plastic analysis of crack line field for mode II plane strain crack

The near crack line analysis method has been used to investigate the exact elastic-plastic solutions of a mode II crack under plane strain condition in an elastic-perfectly plastic solid. The significance of this paper is that the assumptions of the conventional small scale yielding theory have been completely abandoned. The inappropriateness of matching conditions formerly taken at the elastic-plastic boundary ths been corrected as well. By eatching the general solution of the plastic stress (but not the special solution that was adopted) with the exact elastic stresses (but not the crack tip K-dominant field) at the elastic-plastic boundary near the crack line, the plastic stresses, the length of the plastic zone and the unit normal vector of the elastic-plastic boundary, which are sufficiently precise near the crack line region, have been given. The solutions are suitable not only under the condition that the plastic region is sufficiently small but also under the condition that the plastic region is large.     

Evaluation of the elastic-plastic mixity parameters on the base of different crack propagation criteria. Part 1. Preliminary analysis

In order to evaluate the mechanical behavior around small-scale yielding crack tip for both plane strain and plane stress, the asymptotic governing equations and their boundary conditions by considering fracture mechanisms are formulated. A total deformation theory of plasticity with a power-law hardening is used. The analysis of the near-tip fields is carried out for both the maximum tensile and shear stress crack growth direction criteria, as well as for the complete range of mixity parameters and various strain-hardening levels. The new scheme of mixed-mode problem solution is proposed. Realationships between elastic and plastic mixity parameters are given as functions of the crack growth direction criterion and the strain-hardening exponent.Translated from Problemy Prochnosti, No. 1, pp. 60–75, January–February, 2005     

Numerical analysis of blunting of a crack tip in a ductile material under small-scale yielding and mixed mode loading

The blunting of the tip of a crack in a ductile material is analysed under the conditions of plane strain, small-scale yielding, and mixed mode loading of Modes I and II. The material is assumed to be an elastic-perfectly plastic solid with Poisson's ratio being 1/2. The stress and strain fields for a sharp crack under mixed mode loading are first determined by means of elastic-plastic finite element analysis. It is shown that only one elastic sector exists around the crack tip, in contrast with the possibility of existence of two elastic sectors as discussed by Gao. The results obtained for a sharp crack are used as the boundary conditions for the subsequent numerical analysis of crack tip blunting under mixed mode loading, based on slip line theory. The characteristic shapes of the blunted crack tip are obtained for a wide range of Mode I and Mode II combinations, and found to resemble the tip of Japanese sword. Also the stress field around the blunted crack tip is determined.     

Asymptotic deformation and stress fields at the tip of a sharp notch in an elastic-plastic material

The singular elastic-plastic stress, strain and the displacement fields at the tip of a sharp notch for both plane stress and plane strain conditions are investigated analytically. The material is assumed to be governed by the deformation theory of plasticity with linear strain hardening characteristic. Since the elastic strain is retained in the analysis, the final strain and displacement fields can be separated into the elastic and the plastic parts. In the case with zero notch angle, the results reduce to the classical crack problem. The relationship of the amplitude of the near crack tip elastic-plastic field to the elastic far field is obtained. Both mode I and mode II cases are investigated. The mixed mode case is also discussed.     

Analysis of the temperature dependence of the fracture stress of sharp notched bending specimens of a structural steel with respect to intermediate transitions

The flow and fracture stresses, σ_y^g and σ_f^g (δ = crack tip displacement), of sharply notched bending specimens of a structural steel U St 37-1 are measured in the temperature range from full scale to small scale yielding. The best adaption of the experimental results for σ_f^g is obtained by a curve which exhibits an intermediate transition, i.e. which follows in a temperature range between an upper, T_tM¹ and a lower, T_tl¹, transition temperature to the curve σ_y^g(T) for the flow stress with a constant δ = δ₁. This transition corresponds to that of the slip to the twin nucleated fracture. Two analyses [3,5] according to the local fracture stress, σ_f^*, concept show that the amount and the temperature dependence of σ_f^* are somewhat different for both methods, but that both exhibit an increase of σ_f^* in the transition range. It is concluded that each transition in the nucleation mode of the fracture is connected with such a transition in the fracture stress. It may, however, become indistinct or even be covered by the scatter of the experimental points.     

Mixed mode fatigue growth of curved cracks emanating from fastener holes in aircraft lap joints

The finite element alternating method is extended further for analyzing multiple arbitrarily curved cracks in an isotropic plate under plane stress loading. The required analytical solution for an arbitrarily curved crack in an infinite isotropic plate is obtained by solving the integral equations formulated by Cheung and Chen (1987a, b). With the proposed method several example problems are solved in order to check the accuracy and efficiency of the method. Curved cracks emanating from loaded fastener holes, due to mixed mode fatigue crack growth, are also analyzed. Uniform far field plane stress loading on the plate and sinusoidally distributed pin loading on the fastener hole periphery are assumed to be applied. Small cracks emanating from fastener holes are assumed as initial cracks, and the subsequent fatigue crack growth behavior is examined until long arbitrarily curved cracks are formed near the fastener holes under mixed mode loading conditions.     

Plane stress crack-line fields for crack growth in an elastic perfectly-plastic material

Mode-I crack growth in an elastic perfectly-plastic material under conditions of generalized plane stress has been investigated. In the plastic loading zone, near the plane of the crack, the stresses and strains have been expanded in powers of the distance, $\text{y}$, to the crack line. Substitution of the expansions in the equilibrium equations, the yield condition and the constitutive equations yields a system of simple ordinary differential equations for the coefficients of the expansions. This system is solvable if it is assumed that the cleavage stress is uniform on the crack line. By matching 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 depends on crack-line position and time, and applies from the propagating crack tip up to the moving elastic-plastic boundary. Numerical results are presented for the edge crack geometry.     

A perturbation analysis of combined mode I and III dynamic crack propagation

Summary In the present paper the asymptotic stress and deformation fields of dynamic crack extension in materials with linear plastic hardening under combined mode I (plane strain and plane stress) and anti-plane shear loading conditions (mode III) are investigated. The governing equations of the asymptotic crack-tip fields are formulated from two groups of angular functions, one for the in-plane mode and the other for the anti-plane shear mode. It was assumed that all stresses and deformations are of separable functional forms ofr and , which represent the polar coordinates centered at the actual crack tip. Perturbation solutions of the governing equations were obtained. The singularity behavior and the angular functions of the crack-tip in-plane and the anti-plane stresses obtained from the perturbation analysis show that, regardless of the mixity of the crack-tip field and the strain-hardening, the in-plane stresses under the combined mode I and mode III conditions have stronger singularity in the whole mixed mode steady-state crack growth than that of the anti-plane shear stresses. The anti-plane shear stresses perturbed from the plane strain mode I solutions lose their singularity for small strain hardening, whereas the angular stress functions perturbed from the plane stress mode I have a nearly analogous uniform distribution feature compared to pure mode III cases. An obvious deviation from the unperturbed solution is generally to be observed under combined plane strain mode I and anti-plane mode III conditions, especially for a large Mach number in a material with small strain-hardening; but not under plane stress and mode III conditions. The crack propagation velocity decreases the singularities of both pure mode and perturbed crack-tip fields.     

Numerical simulations of dynamic crack growth along an interface

Dynamic crack growth is analyzed numerically for a plane strain bimaterial block with an initial central crack. The material on each side of the bond line is characterized by an isotropic hyperelastic constitutive relation. A cohesive surface constitutive relation is also specified that relates the tractions and displacement jumps across the bond line and that allows for the creation of new free surface. The resistance to crack initiation and the crack speed history are predicted without invoking any ad hoc failure criterion. Full finite strain transient analyses are carried out, with two types of loading considered; tensile loading on one side of the specimen and crack face loading. The crack speed history and the evolution of the crack tip stress state are investigated for parameters characterizing a PMMA/Al bimaterial. Additionally, the separate effects of elastic modulus mismatch and elastic wave speed mismatch on interface crack growth are explored for various PMMA-artificial material combinations. The mode mixity of the near tip fields is found to increase with increasing crack speed and in some cases large scale contact occurs in the vicinity of the crack tip. Crack speeds that exceed the smaller of the two Rayleigh wave speeds are also found.     

An Experimental Study on the Elastic–Plastic Fracture in a Ductile Material under Mixed-Mode Loading

Abstract:  Stereo vision was used to measure the crack-tip parameters, such as J integral, plastic mixity and elastic mixity of mixed-mode fracture specimens, and to study the applicability of the Shih's plane strain solution to the mixed-mode crack-tip fields. The fracture specimen used in this study was a compact tension shear (CTS) specimen made of 2024-O aluminum. The in-plane strain and stress fields near the mixed-mode crack tip of the CTS specimen were determined using the deformation field measured by the stereo vision. It is observed that the J integral values computed along rectangular contours surrounding the mixed-mode crack-tip approach constant values after r / h  > 0.5. The in-plane strains determined experimentally at several points near the crack tip and at several radial lines emerging from the crack tip are compared with the values calculated using Shih's plane-strain solution and the HRR slope, named after the investigations of Hutchinson, Rice and Rosengren respectively. It is found that the measured values follow the trends of the Shih's plane-strain solution. The elastic mixity evaluated using the measured crack-tip stress fields is close to that obtained from analytical solution. However, the evaluated plastic mixity deviates from the analytical solution.     

Constraint effects on ductile fracture processes near a notch tip under mixed-mode loading

In this work, the effect of constraint on ductile fracture process of microvoid growth and coalescence near a notch tip in a ductile material under mode I and mixed mode loading (involving modes I and II) is investigated. To this end, two sets of finite element simulations are carried out under two-dimensional plane strain conditions. In the first set, a modified boundary layer formulation is employed in which the mixed mode elastic K–T field is prescribed as remote boundary conditions. Several analyses are carried out corresponding to different values of T-stress and remote elastic mode-mixity. Next, ductile four point bend specimens subjected to mode I and mixed mode loading are considered. In both sets of simulations, the interaction between a notch tip and a pre-nucleated hole ahead of it is modelled. The background material is represented by the Gurson constitutive model and micro-void nucleation at uniformly distributed small scale particles is also taken into account. The accumulation of matrix plastic strain and porosity in the ligament between the notch tip and the hole as well as the growth of the hole are studied. Finally, the effect of crack tip constraint on the relationship between the fracture toughness and mode mixity is examined.     

The plastic zone ahead of a mode II plane strain crack

The problem of a mode II plane strain crack in an elastic perfectly-plastic solid is analysed. It is shown that the elastic singular field for a crack cannot be matched with the near-tip field. Hence, the usual definition of small scale yielding is inappropriate. By matching the stress field in the plastic zone with the elastic dominant field for a blunt crack near the crack line at the elastic-plastic boundary, the elastic-plastic boundary ahead of the crack is determined.

---

Résumé On analyse le problème d'une fissure de mode II en état de déformation plane dans un solide élastique parfaitement plastique. On montre que le champ singulier élastique relatif à une fissure ne peut être atteint par le champ régnant au voisinage de l'extrémité de la fissure considérée. Dès lors, la définition habituelle de l'écoulement platique à petite échelle n'est pas adéquate. On détermine la limite élastique-plastique en avant de la fissure en faisant se correspondre le champ de contraintes dans la zone plastique avec le champ principal élastique correspondant à une fissure arrondie située près de la ligne de la fissure.

     

Plastic notch stress intensity factors for pointed V-notches under antiplane shear loading

The paper deals with a work-hardening, elastic–plastic, stress analysis of pointed V-notches under antiplane shear deformation loading both under small and large scale yielding. Stress and strain field intensities are expressed in terms of plastic Notch Stress Intensity Factors, which are analytically linked to the corresponding linear elastic ones under small scale yielding. The near tip stress and strain fields are then used to give closed-form expressions for the Strain Energy Density over a circular sector surrounding the notch tip, and for the J-integral parameter, both as a function of the relevant plastic NSIFs, these expressions being valid both under small and large scale yielding.     

Specimen proposals for mixed mode testing of adhesive layer

The experimental methods to determine the fracture properties for adhesives under mixed mode loading is not as well established as compared to such methods for adhesives under pure mode loading. Some controversies exist regarding the decomposition of the mode mixity. For a flexible linear elastic adhesive, the mode mixity of a single-layer adhesive joint is directly related to the deformation of the adhesive layer at the crack tip. The governing equations for linear elastic single-layer adhesive joints show that the mode mixity depends on the external loads, the properties of the adherends and often on the flexibility of the adhesive layer. This implies some fundamental problems that have to be addressed before an experimental method can be established. The purpose of this paper is to investigate different specimen configurations for mixed mode loading. Requirements for the design of a specimen configuration are given. A new specimen configuration is proposed based on the geometry of a semi-infinite symmetric DCB-specimen. According to this study, the proposed test specimen offers exceptional flexibility, variety and stability.