Plane-strain mixed-mode near-tip fields in elastic perfectly plastic solids under small-scale yielding conditions

Within the context of the small-strain approach, plane-strain mixed-mode near-tip fields of a stationary crack in an elastic perfectly plastic Mises solid under small-scale yielding conditions are examined by finite element methods. Steady-state stress fields in the immediate vicinity of the crack tip are obtained as the remote loading of the elastic K-field increases. Asymptotic crack-tip solutions consisting of constant stress sectors, centered fan sectors, and an elastic sector are then constructed accordingly. The asymptotic crack-tip stress solutions agree well with the numerical results for a whole spectrum of mixed-mode loadings. Our mixed-mode near-tip solution with an elastic sector differs from that of Saka et al. by one (plastic) constant stress sector situated between the elastic sector and the neighbouring fan sector. The effect of the existence of the elastic sector on the near-tip fields is discussed in the light of the computational results. The plastic mixity factor of the near-tip field is given as a function of the elastic mixity factor of the prescribed K-field. This function is well bounded by that of the perfectly plastic limit of the corresponding solutions for power-law hardening materials given by Shih. Some issues pertaining to the numerical procedures such as the implementation of the small-scale yielding assumption are also addressed.     

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.     

Fracture mechanics parameters for cracks on a slightly undulating interface

Typical bimaterial interfaces are non-planar due to surface facets or roughness. Crack-tip stress fields of an interface crack must be influenced by non-planarity of the interface. Consequently, interface toughness is affected. In this paper, the crack-tip fields of a finite crack on an elastic/rigid interface with periodic undulation are studied. Particular emphasis is given to the fracture mechanics parameters, such as the stress intensity factors, crack-tip energy release rate, and crack-tip mode mixity. When the amplitude of interface undulation is very small relative to the crack length (which is the case for rough interfaces), asymptotic analysis is used to convert the non-planarity effects into distributed dislocations located on the planar interface. Then, the resulting stress fields near the crack tip are obtained by using the Fourier integral transform method. It is found that the stress fields at the crack tip are strongly influenced by non-planarity of the interface. Generally speaking, non-planarity of the interface tends to shield the crack tip by reducing the crack-tip stress concentration.     

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     

Partially stiffened elastic half-plane with an edge crack

A technique, using the Brazilian disk specimen, for measuring the fracture toughness of unidirectional fiber-reinforced composites, over the entire range of crack-tip mode mixities, was developed. The fracture toughness of a graphite/epoxy fiber-reinforced composite was measured, under both mode-I and mode-II loading conditions. We found that for certain material orientations the mode-II fracture toughness is substantially higher than the mode-I toughness. The complete dependence of the fracture toughness on the crack-tip mixity was determined for particular material orientations and the phenomenological fracture toughness curves were constructed. Using the Brazilian disk specimen, together with a hydraulic testing machine, the fracture toughness of the composite under moderate loading rates was measured. We observed that the mode-I fracture toughness was not sensitive to the loading rate at the crack tip, K, while the mode-II ‘dynamic’ fracture toughness increased approximately 50 percent over the quasi-static fracture toughness. A qualitative explanation of the dependency of fracture toughness on crack-tip loading rate is discussed. Finally, a mechanical fracture criterion, at the microscopic level, which governs the crack initiation under mixed-mode loading conditions is presented; these theoretical predictions closely follow the trend of experimental measurements. This revised version was published online in August 2006 with corrections to the Cover Date.     

Strain and stress fields near the blunted tip of a crack under mixed mode loading and the implications for fracture

The strain distribution in the vacinity of a blunted crack-tip is analysed by slip line theory under the conditions of plane-strain, small-scale yielding, and mixed-mode loading of Modes I and II. A generalized crack-tip opening displacement is introduced by which the strain and stress fields near the blunted crack-tip are determined uniquely over a wide range of Mode I and II combinations. Also, coupled experimental and finite-element analyses under the condition of large-scale yielding reveal that the initiation of stable crack growth occurs when the generalized crack-tip opening displacement attains a critical value which is constant for the material tested. The finite-element analysis is based on the finite deformation theory of elastic-plastic materials. The generalized crack-tip opening displacement criterion is found to be superior to the J-integral and the usual COD for the characterization of the initiation of stable crack growth. The plastic work in a small circular region at the crack-tip is found to be equivalent to the generalized crack-tip opening displacement, as a fracture criterion.     

Mixed-mode fracture of ductile thin-sheet materials under combined in-plane and out-of-plane loading

Cracks in thin structures often are subjected to combined in-plane and out-of-plane loading conditions leading to complex mixed mode conditions in the crack tip region. When applied to ductile materials, large out-of-plane displacements make both experimentation and modeling difficult. In this work, the mixed-mode behavior of thin, ductile materials containing cracks undergoing combined in-plane tension (mode I) and out-of-plane shear (mode III) deformation is investigated experimentally. Mixed-mode fracture experiments are performed and full, three-dimensional (3D) surface deformations of thin-sheet specimens from aluminum alloy and steel are acquired using 3D digital image correlation. General characteristics of the fracture process are described and quantitative results are presented, including (a) the fracture surface, (b) crack path, (c) load-displacement response, (d) 3D full-field surface displacement and strain fields prior to crack growth, (e) radial and angular distributions of the crack-tip strain fields prior to crack growth and (f) singularity analysis of the crack-tip strains prior to crack growth. Results indicate that the introduction of a mode III component to the loading process (a) alters the crack tip fields relative to those measured during nominally mode I loading and (b) significantly increases the initial and stable critical crack-opening-displacement. The data on strain fields in both AL6061-T6 aluminum and GM6208 steel consistently show that for a given strain component, the normalized angular and radial strains at all load levels can be reasonably represented by a single functional form over the range of loading considered, confirming that the strain fields in highly ductile, thin-sheet material undergoing combined in-plane tension and out-of-plane shear loading can be expressed in terms of separable angular and radial functions. For both materials, the displacement and strain fields are (a) similar for both mixed-mode loading angles Φ = 30° and Φ = 60° and (b) different from the fields measured for Mode I loading angle Φ = 0°. Relative to the radial distribution, results indicate that the in-plane strain components do not uniformly exhibit the singularity trends implicit in the HRR theory.     

Experimental investigations of three-dimensional effects near a crack tip using computer vision

An experimental study of the near tip deformation fields for a Single Edge-Cracked specimen (SEC) has been completed. The surface deformation fields for a thin SEC plexiglas specimen have been obtained in the region 0.15 <- r/t <- 0.60, where r is radial distance from the crack tip and t is the specimen thickness, by using a novel computer vision method. The results of the study indicate that the value of the J-Integral obtained from the measured surface deformation, and under the assumption of plane stress linearly elastic behavior, is essentially path independent over a region that is considered to be within the crack tip three-dimensional zone by researchers who have performed tests on other materials.However, the near-tip, in-plane displacement field does show a deviation from the traditional linear elastic crack tip singularity over the same region, suggesting the need for further studies to address this inconsistency.     

Dynamic fields generated by rapid crack growth

Stresses and strains near a rapidly propagating crack tip are affected by the mass density of the material. This paper starts with a brief summary of analytical results for near-tip dynamic fields as predicted by linear elastic fracture mechanics. Next, exact expressions are derived for dynamic crack-line strains, for mode-III crack propagation in a nonlinear elastic material and in an elastic perfectly-plastic material. These expressions are valid on the crack line from the moving crack tip to the moving boundary with the region of linearly elastic deformation. For steady-state crack growth, a critical strain criterion is used to compute the relation between external load and crack tip speed. The required external load increases with crack-tip speed.     

An Experimental/Analytical Compression of Three-Dimentional Deformations at the Tip of a Crack in a Plastically deforming Plate III: Material Characterization and finite element analysis

A comparison between the three-dimensional experimental and numerical displacement fields surrounding a notch/crack in a ductile 4340 steel tested in three-point bending is presented. Excellent agreement between computed and measured deformations exists at load levels below 50 to 75 percent of ultimate loads. Experimentally determined crack tunnel profiles are included in the finite element model through nodal release; the evidence of the crack tunnel appears in the displacements at the surface. It is shown that surface measurements of unloading reveal specimen-internal failure initiation in the form of tunneling. Out-of-plane deformations deviate from analytical values earlier than in-plane values; this observation compromises the accuracy with which predictions of in-plane crack tip variables can be made when they are based on measured out-of-plane deformations (caustics, gradient sensing) once significant plasticity arises. Comparison is made between J-integral values calculated from the external boundary conditions and from a domain integral. The tunneling tests provide a method of estimating a critical value of J. The stress intensity factor governs the deformation in the elastic regime, but, because of the finite notch- tip radius underlying the experimental configuration, the HRR field does not describe the deformation well under plastic conditions. Comparison of numerical simulations with and without tunneling provide insight into criteria that could be used to implement an implicit crack propagation scheme into the numerical model.     

Geometrical relationships between side necking and plastic zone size near a crack tip in ductile metals. Part I: Modified boundary layer solutions

The plastic zone size is regarded as the measure of a material's resistance, and it also determines the fracture behaviour. Recently, stereo digital speckle photography (SDSP) has been found to be useful for measuring in-situ the side necking developed on the lateral surfaces of a specimen during a standard fracture test procedure for J_IC . Because plastic deformation occurs without any volume change, the in-plane plastic zone developed around a crack tip should be accompanied by out-of-plane deformation, that is, side necking. With the aid of the new measurement technique, side necking is expected to act as a gauge for indicating the plastic zone size. As a preliminary study, the geometrical relationships between side necking and the plastic zone size near a crack tip in ductile metals are explored by using a finite element model with modified boundary conditions. As parameters representing the geometrical similarity between side necking and the plastic zone, the shapes of each region, the distances from the crack tip to the boundaries of each region, r_p and r_s and the areas of each region, A_p and A_s are examined for their sensitivities to variables such as mode mixity, hardening exponent and so on. Among them, the areas, A_p and A_s seem to be the best for application because an excellent linearity between them is maintained in a wide range of mode mixity and load level regardless of the hardening exponent, specimen thickness and yield stress.     

The plastic stress ahead of a mixed mode I and II plane stress crack

The small scale yielding for mixed mode I and II plane stress crack problems in elastic perfectly-plastic solids is analysed by considering the stress field near the crack line. By expanding the stresses near the crack line and 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 problem is reduced to solving a system of nonlinear algebraic equations. The relationship between the near-field mixity parameter M^p and the far-field mixity parameter M^e is detennined by solving the system of equations numerically. Analogous to Shih's calculation by the finite element method for the small scale yielding of mixed mode plane strain crack problems, the numerical results indicate that the shift from a mixed mode to a pure mode may not be a smooth one.     

Mechanisms and modeling of low cycle fatigue crack propagation in a pressure vessel steel Q345

The fatigue process near crack is governed by highly concentrated strain and stress in the crack tip region. Based on the theory of elastic–plastic fracture mechanics, we explore the cyclic J-integral as breakthrough point, an analytical model is presented in this paper to determine the CTOD for cracked component subjected to cyclic axial in-plane loading. A simple fracture mechanism based model for fatigue crack growth assumes a linear correlation between the cyclic crack tip opening displacement (ΔCTOD) and the crack growth rate (da/dN). In order to validate the model and to calibrate the model parameters, the low cycle fatigue crack propagation experiment was carried out for CT specimen made of Q345 steel. The effects of stress ratio and crack closure on fatigue crack growth were investigated by elastic–plastic finite element stress–strain analysis of a cracked component. A good comparison has been found between predictions and experimental results, which shows that the crack opening displacement is able to characterize the crack tip state at large scale yielding constant amplitude fatigue crack growth.     

Finite-element modelling of crack propagation in elastic–plastic media: Part II Horizontal subsurface cracks

Horizontal subsurface cracks in an elastic–plastic material are analysed using finite-element techniques. The sliding surface is modelled as a rigid cylinder. The effect of such parameters as the friction between the cylinder and the material being indented, the elastic and plastic modulus of the material and the depth of crack location on the J-integral values at the left and right tips of a horizontal subsurface crack is considered. The prospective crack propagation direction is taken as the direction along which the J integral assumes a maximum as the indenter slides along the material surface. The left and right tip cracks were found likely to propagate at about 10° to the horizontal. This propagation direction was found to depend strongly on the location of the crack. Both crack tips are expected to propagate closer to the vertical direction as the depth of crack location is reduced. Also, horizontal cracks closer to the surface are found to have higher J integral values. While friction between the slider and the specimen did not affect the crack propagation direction, the crack-tip plasticity reduced the propagation direction, with respect to the horizontal.     

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.     

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.     

Measurement of mixed-mode fracture parameters near cracks in homogeneous and bimaterial beams

An investigation of deformation fields and evaluation of fracture parameters near mixed-mode cracks in homogeneous and bimaterial specimens under elastostatic conditions is undertaken. A modified edge notched flexural geometry is proposed for testing bimaterial interface fracture toughness. The ability of the specimen in providing a fairly wide range of mode mixities is demonstrated through direct optical measurements and a simple flexural analysis. A full field optical shearing interferometry called Coherent Gradient Sensing (CGS) is used to map crack tip deformations in real time. Experimental measurements and predictions based on beam theory are found to be in good agreement. Also, for a large stiffness mismatch bimaterial system, the interface crack initiation toughness is evaluated as a function of the crack tip mode mixity.     

Connecting the essential work of fracture, stress intensity factor, Hill’s criterion in mixed mode I/II loading

Ductile sheet structures are frequently subjected to mixed mode loading, resulting that the structure is under the influence of a mixed mode stress field. Instances of interest are when stable crack growth occurs and when the crack-tip is propagating in this complex mixed-mode condition, prior to final fracture. Purposely designed apparatus was built to test thin-sheets of steel (Grade: DX51D) under mixed-mode I/II. These tests, under plane stress conditions, also investigated the effect of thickness on the specific essential work of fracture or the fracture toughness of the material under quasi-static cracking conditions. The fracture toughness is evaluated under incremental mixed-mode loading conditions. The direction of the propagating crack path and fracture type were observed and discussed as the loading mixity was varied. Whilst the specific essential work of fracture or fracture toughness was obtained using the energy approach, the theoretical analysis of the fracture type and direction of crack path were based on the crack tip stresses and fracture criterions of maximum hoop stress and maximum shear stress along with the utilisation of Hill’s theory. For mixed-mode I/II loading, the variation in the fracture toughness contributions ratios are evaluated and used predicatively using the established energy criterion approach to the crack tip stress intensity approach. The comparison between the theoretical directions of the crack path, failure mode propagation are in good agreement with those obtained from experimental testing indicating the definite link between both approaches.     

Influence of local fracture energy distribution on maximum fracture load of three-point-bending notched concrete beams

The maximum fracture load of a notched concrete beam has been related to the local fracture energy at the cohesive crack tip region analytically in this paper, and then the correlation between the size effects on the maximum fracture loads and the RILEM specific fracture energy is established. Two extreme conditions have been established, namely zero crack-tip bridging with zero local fracture energy and maximum crack-tip bridging with the maximum size-independent fracture energy. It is concluded that the local fracture energy at the crack tip region indeed varies with the initial crack length and the size of specimen. The tri-linear model for the local fracture energy distribution is confirmed by using the proposed simple analytical solution.     

DCB-specimen loaded with uneven bending moments

A double cantilever beam specimen loaded with uneven bending moments (DCB-UBM) is proposed for mixed mode fracture mechanics characterisation of adhesive joints, laminates and multilayers. A linear elastic fracture mechanics analysis gives the energy release rate and mode mixity analytically for both isotropic and orthotropic materials. By varying the ratio between the two applied moments, the crack tip stress state can be varied from pure mode I to pure mode II for the same specimen geometry. The specimen allows stable crack growth. A special test fixture is developed to create uneven bending moments. As a preliminary example, the DCB-UBM specimen was used for characterising fracture of adhesive joints between two laminates of thermoset glass fibre reinforced plastic.