Effect of Plate Thickness on Crack-Tip Plasticity

This paper presents an analytical method for determining the three-dimensional stress fields in plates with a through-the-thickness crack, especially under elastic-plastic conditions. Using the generalised plane strain theory in conjunction with the deformation theories of plasticity, exact solutions are obtained for the effects of plate thickness on the crack-tip plastic zone size and a plastic constraint factor, which is shown to correlate well with published finite element solutions.     

A strain-based Dugdale model for cracks under generally yielding conditions

To develop an analytical method for quantifying the growth behaviour of short cracks embedded in notch plastic zones, a critical assessment of the Dugdale model is first made by comparison against finite element analysis for an edge-cracked plate subjected to an applied strain varying linearly along the crack path. It is shown that the conventional stress-based Dugdale model provides accurate estimates for the crack-tip opening displacement and the plastic zone size provided that the applied strain does not exceed one third of the yield strain. These estimates become significantly inaccurate at higher strain levels. To overcome this limitation of the conventional model, a strain-based implementation of the Dugdale model is proposed in which the conventional equilibrium equation is replaced by strain compatibility. Comparison with finite element results shows that this strain-based model provides accurate values for both the crack-tip-opening displacement and the plastic zone size for applied strains up to four times the yield strain and with no evidence of decreasing accuracy with increasing strain. Furthermore, it is shown that the relevant plastic constraint factor to be used for plane strain is that appropriate for the notch plastic zone in the absence of a crack, rather than the more usual choice which is appropriate only for small-scale yielding conditions. This provides a practical and physically plausible approach for extending the scope of current predictive software for fatigue crack growth based on the Dugdale model to include conditions of large-scale yielding.     

Three-dimensional numerical simulations of dynamic fracture in silicon carbide reinforced aluminum

The three-point bending test by Kolsky-bar apparatus is a convenient technique to test the dynamic fracture properties of materials. This paper presents detailed three-dimensional finite element simulations of a silicon particle reinforced aluminum (SiC_p/Al) experiment (Li et al., [Proceedings of the US Army Symposium on Solid Mechanics]. In the simulations, the interaction between the input bar and the specimen is modeled by coupled boundary conditions. The material model includes large plastic deformations, strain-hardening and strain-rate hardening mechanisms. Furthermore, crack initiation and propagation processes are simulated by a cohesive element model. The simulation results quantitatively agree with the experimental measurements on three fronts: (1) the structural response of the specimen, (2) the time of unstable crack propagation, and (3) the local deformations at the crack-tip zone. The simulations reveal crack propagation characteristics, including crack-tip plastic deformation, crack front curving, and crack velocity profile. The effectiveness of Kolsky-bar type fracture tests is verified. It is shown that a rate-independent cohesive model can describe the complicated dynamic elastic–plastic fracture process in the SiC_p/Al material.     

Modelling of mode-I stable crack growth under hydrogen assisted stress corrosion cracking

The paper presents a new strategy based on combined analytical and finite element (FE) solution to hydrogen assisted stress corrosion crack growth. The diffusion process is solved analytically through both one-and two-dimensional modelling. These solutions are adopted with two-dimensional FE based cohesive zone model of crack extension study. The results fit well with published experimental data and show improvement over the predictions by full FE approach. The new solution approach helps to reduce time required for simulation/computation. The study has produced a relationship between concentration dependent reduction in cohesive strength and plastic strain rate.     

A finite element investigation of quasi-static and dynamic asymptotic crack-tip fields in hardening elastic-plastic solids under plane stress

The asymptotic structures of crack-tip stress and deformation fields are investigated numerically for quasi-static and dynamic crack growth in isotropic linear hardening elastic-plastic solids under mode I, plane stress, and small-scale yielding conditions. An Eulerian type finite element scheme is employed. The materials are assumed to obey the von Mises yield criterion and the associated flow rule. The ratio of the crack-tip plastic zone size to that of the element nearest to the crack tip is of the order of 1.6 × 10⁴. The results of this study strongly suggest the existence of crack-tip stress and strain singularities of the type r ^s (s < 0) at r=0, where r is the distance to the crack tip, which confirms the asymptotic solutions of variable-separable type by Amazigo and Hutchinson [1] and Ponte Castañeda [2] for quasi-static crack growth, and by Achenbach, Kanninen and Popelar [3] for dynamic crack propagation. Both the values of the parameter s and the angular stress and velocity field variations from the present full-field finite element analysis agree very well with those from the analytical solutions. It is found that the dominance zone of the r ^s-singularity is quite large compared to the size of the crack-tip active plastic zone. The effects of hardening and inertia on the crack-tip fields as well as on the shape and size of the crack-tip active plastic zone are also studied in detail. It is discovered that as the level of hardening decreases and the crack propagation speed increases, a secondary yield zone emerges along the crack flank, and kinks in stress and velocity angular variations begin to develop. This dynamic phenomenon observable only for rapid crack growth and for low hardening materials may explain the numerical difficulties, in obtaining solutions for such cases, encountered by Achenbach et al. who, in their asymptotic analysis, neglected the existence of reverse yielding zones along the crack surfaces.     

A perturbation solution for a crack in a power-law material under gross yielding

Abstract In order to develop an analytical method for quantifying the plastic-blunting behaviour of a short crack embedded in a notch plastic zone, the perturbation solution of He and Hutchinson is extended to include the effect of strain gradient. An edge-cracked plate subjected to a linearly varying remote strain is considered in this work to simulate the plastic deformation associated with a small crack at a notch root. The strain hardening of the material is assumed to obey a power-law. Comparison with finite-element (FE) computations demonstrates that this perturbation solution provides accurate values for the crack-tip opening displacement (CTOD) under gross-yielding conditions for a range of hardening parameters.     

A moving cohesive interface model for fracture in creeping materials

We present a new cohesive interface model for quasi-static creep crack growth that is implemented within a moving-grid finite element model. A pseudo crack tip separates the cohesive process zone from the free surfaces of the crack. The moving-grid formulation models continuous crack advance by describing relative motion between the pseudo crack tip and the material. This eliminates the need for extensive mesh refinement away from the current crack-tip location and supports both transient and direct steady-state solutions. A traction-separation law determines the energetic properties of the decohesion process and generates a simple criterion for crack advance. The new formulation remedies a problem in earlier models which permit a crack to heal on unloading. Numerical examples demonstrate the moving cohesive interface model in studies of steady-state crack growth. Adaptive grid refinement is used to control the accuracy of the solution.     

A cohesive finite element formulation for modelling fracture and delamination in solids

In recent years, cohesive zone models have been employed to simulate fracture and delamination in solids. This paper presents in detail the formulation for incorporating cohesive zone models within the framework of a large deformation finite element procedure. A special Ritz-finite element technique is employed to control nodal instabilities that may arise when the cohesive elements experience material softening and lose their stress carrying capacity. A few simple problems are presented to validate the implementation of the cohesive element formulation and to demonstrate the robustness of the Ritz solution method. Finally, quasi-static crack growth along the interface in an adhesively bonded system is simulated employing the cohesive zone model. The crack growth resistance curves obtained from the simulations show trends similar to those observed in experimental studies     

Plastic Deformation at the Tip of a Tensile Crack in a Non-linear Kinematic Hardening Material

The plastic deformation at the tip of a tensile crack in a non-linear kinematic hardening material under small-scale yielding conditions is investigated, with a view to quantifying the functional dependence of crack-tip plastic blunting size on material's strain hardening parameters. It is shown by dimensional analysis that, for materials being characterised by the Armstrong-Frederick non-linear kinematic hardening rule, the crack-tip blunting parameter depends parametrically on only two non-dimensional parameters; the functional dependence is determined using a parametric finite element analysis.     

Constraint effects on plastic crack-tip fields for plane strain mode-I, II and III cracks in non-hardening materials

To explore constraint effects on fully plastic crakc-tip fields, analytical solutions are examined for mode-I, II and III loading in non-hardening materials under plane strain conditions. The results reveal that under mode-II and III loading the crack-tip stress fields are unique, and thus can be characterized by a `single parameter'. Under mode-I loading, however, the crack-tip stress field is non-unique but can be characterized by two sets of solutions or `two parameters'. One set of the solutions is the well-known Prandtl field and the other is a plastic T-stress field. This conclusion corroborates the observation of McClintock (1971) that the slip-line field is non-unique for plane strain tensile cracks. A two-term plastic solution which combines the Prandtl field and the plastic T-stress field with two parameters B ₁ and B ₂ can then characterize the crack-tip stress field of plane strain mode-I crack over the plastic region and quantify the magnitude of crack-tip constraints. These characters are similar to those for hardening materials. Analyses and examples show that the two-term plastic solution can match well with the slip-line field or finite element results over plastic region. Thus the parameters B ₁ and B ₂ can be used to characterize the constraint level for mode-I finite-sized crack specimens in non-hardening materials under plane strain conditions.     

Finite deformation finite element analyses of tensile growing crack fields in elastic-plastic material

Finite deformation finite element analyses of plane strain stationary and quasi-statically growing crack fields in fully incompressible elastic-ideally plastic material are reported for small-scale yielding conditions. A principal goal is to determine the differences between solutions of rigorous finite deformation formulation and those of the usual small-displacement-gradient formulation, and thereby assess the validity of the (nearly all) extant studies of ductile crack growth that are based on a small-displacement-gradient formulation. The stationary crack case with a significantly blunted tip is studied first; excellent agreement in stress characteristics at all angles about the crack tip and up to a radius of about three times the crack tip opening displacement is shown between Rice and Johnson's [1] approximate analytical solution and our numerical solution. Outside this radius, the numerical results agree very well with Drugan and Chen's [2] small-displacement-gradient analytical characteristics solution in the region of principal plastic deformation. Thus we identify accurate analytical representations for the stress field throughout the plastic zone of a blunted stationary crack. For the growing crack case, the macroscopic difference in crack tip opening profiles between previous small-displacement-gradient solutions and the present results is shown to be negligible, as is the difference in the stress fields in plastic regions. The stress characteristics again agree very well with analytical results of [2]. The numerical results suggest—in agreement with a recent analytical finite deformation study by Reid and Drugan [3]—that it is the finite geometry changes rather than the additional spin terms in the objective constitutive equation that cause any differences between the small-displacement-gradient and the finite deformation solutions, and that such differences are nearly indistinguishable for growing cracks.     

Analytical and numerical study of cohesive zones for a crack in an adhesive layer between identical isotropic materials

A crack in a thin adhesive elastic-perfectly plastic layer between two identical isotropic elastic half-spaces is considered. Uniformly distributed normal stress is applied to the substrates at infinity. First, stress distribution in the cohesive zones and the J-integral values are defined numerically by the finite element method (FEM). Further, a mathematical formulation of the problem is given and its analytical solution is proposed. It is assumed that, at the crack continuations, there exist cohesive zones. The interlayer thickness is neglected since it is much smaller than the crack length. The distribution of the normal stress, which was obtained by means of the FEM, is now approximated by a piecewise-constant function and assumed to be applied at the faces of the cohesive zones. The formulated problem is solved analytically and an equation for determination of the cohesive zone lengths is derived. Also, closed expressions for the crack tip opening displacement and for the J-integral are obtained in an analytical form. These parameters are found with respect to the values of the normal stress applied at infinity. Finally, a universal approximating function, which describes the stress distribution in the cohesive zones, is constructed. This function depends on the ratio between the interlayer thickness and the crack length and on the ratio between the normal stress applied at infinity and the yield limit of the interlayer’s material. Once again, the problem is solved analytically, but this time for the stress distribution prescribed by the universal approximating function. The cohesive zone lengths, the values of the crack tip opening displacement and of the J-integral are calculated. A comparative analysis of the obtained results is carried out. A good agreement of the J-integral values calculated by means of the developed analytical models and by the associated finite element analysis is demonstrated.     

A finite element analysis of fracture initiation in ductile/brittle periodically layered composites

A near-tip plane strain finite element analysis of a crack terminating at and normal to the interface in a laminate consisting of alternate brittle and ductile layers is conducted under mode-I loading. The studies are carried out for a system representing steel/alumina composite laminate. The Gurson constitutive model, which accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence, is employed within the framework of small deformation plasticity theory. Evolution of plastic zone and damage in the ductile layer is monitored with increasing load. High plastic strain localization and microvoid damage accumulation are found to occur along the brittle/ductile interface at the crack-tip. Fracture initiation in the ductile phase is predicted and the conditions for crack renucleation in the brittle layer ahead of the crack are established for the system under consideration. Ductile fracture initiation has been found to occur before plasticity spreads in multiple ductile layers. Effects of material mismatch and yield strength on the plastic zone evolution are briefly discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

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.     

An incremental analysis of plane strain fully plastic crack growth in strain-hardening materials under extension

Analyses of inelastic fracture have mainly followed two directions. One is crack-tip-field analysis in strain-hardening materials (e.g., the HRR solution). The other is whole field analysis in non-hardening materials (e.g., McClintock's slip-line approach). In this paper, a theoretical approach that combines the two directions is presented to account for large crack growth. As an example, plane strain mode I fully plastic crack growth in a single-edge-cracked-specimen under extension is analyzed. The incremental analysis based on the deformed configuration is developed for large crack growth in strain-hardening materials. A kinematically admissible displacement increment field with crack-tip singularity is first constructed in a presumed symmetrical triangular deformation zone extending from the crack tip to the back flank of the ligament (whole field). Then the size of the deformation zone is determined by minimizing the total force in each incremental step. The strain histories of all material points in the ligament are traced and a fracture criterion based on the hole growth theory is applied. The comparisons between the present study and the experiments existing in the literature show the validity of the present approach.     

An extended finite element method with analytical enrichment for cohesive crack modeling

A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.     

Fully plastic crack-tip fields for plane strain shallow-cracked specimens under pure bending

Detailed finite element (PE) analyses are performed to study the effect of crack depth on crack-tip constraint at full yielding for pure bending of plane strain single-edge-cracked specimens. Analyses are based on small-strain formulations and perfect plasticity. The crack depth a/W ranges from 0.1 to 0.7, and the deformation is applied up to the limiting state of full plasticity where crack-tip stresses reach steady-state limiting values.At load levels smaller than the limit load (contained yielding), the crack-tip constraint (stress triaxiality) gradually decreases as a/W decreases, but, at load levels close to the limit load (or at the limit load), it decreases very sharply. In terms of a/W, tractable closed-form approximations for fully plastic crack-tip stress and strain fields are proposed, and fully plastic values of crack-tip stresses are re-phrased in terms of the Q-parameter [1, 2]. The role of crack-tip strains on fracture of shallow-cracked bending specimens is briefly discussed.     

Thermal effect of plastic dissipation at the crack tip on the stress intensity factor under cyclic loading

Plastic dissipation at the crack tip under cyclic loading is responsible for the creation of an heterogeneous temperature field around the crack tip. A thermomechanical model is proposed in this paper for the theoretical problem of an infinite plate with a semi-infinite through crack under mode I cyclic loading both in plane stress or in plane strain condition. It is assumed that the heat source is located in the reverse cyclic plastic zone. The proposed analytical solution of the thermo-mechanical problem shows that the crack tip is under compression due to thermal stresses coming from the heterogeneous stress field around the crack tip. The effect of this stress field on the stress intensity factor (its maximum and its range) is calculated analytically for the infinite plate and by finite element analysis. The heat flux within the reverse cyclic plastic zone is the key parameter to quantify the effect of dissipation at the crack tip on the stress intensity factor.     

A FINITE ELEMENT ANALYSIS OF A CRACK GROWING UNDER CYCLIC LOADING

Abstract— Crack growth under cyclic loading has been studied by the finite element method. The calculation was made for plane stress conditions. The crack tip zone was modelled as a cohesive zone. The displacement of the free crack surface during unloading was found to be governed by the surrounding continuum and was independent of the details in the fracture zone. This means that crack closure upon unloading is directly related to the ultimate separation, of the cohesive zone, which in turn controls the residual plastic deformation left in the wake of the growing crack. If the distance over which closure takes place is rather small, closure may be very difficult to detect by the compliance technique.