A finite element analysis of plane strain dynamic crack growth in materials displaying the Bauschinger effect

In this paper dynamic crack growth in an elastic-plastic material is analyzed under mode I plane strain small-scale yielding conditions using a finite element procedure. The main objective of this paper is to investigate the influence of anisotropic strain hardening on the material resistance to rapid crack growth. To this end, materials that obey an incremental plasticity theory with linear isotropic or kinematic hardening are considered. A detailed study of the near-tip stress and deformation fields is conducted for various crack speeds. The results demonstrate that kinematic hardening does not oppose the role of inertia in decreasing the plastic strains and stresses near the crack tip with increase in crack speed to the same extent as isotropic strain hardening. A ductile crack growth criterion based on the attainment of a critical crack opening displacement at a small micro-structural distance behind the tip is used to obtain the dependence of the theoretical dynamic fracture toughness with crack speed. It is found that for any given level of strain hardening, the dynamic fracture toughness displays a much more steep increase with crack speed over the quasi-static toughness for the kinematic hardening material as compared to the isotropic hardening case.     

Mode I plane stress crack-tip constraint for elastic-plastic materials with different strain hardening

Finite element analyses and simulations have been undertaken to investigate the triaxial constraint in the crack-tip regions of a stationary crack and a steady-state growing crack under mode I plane stress for elastic-plastic materials with different strain hardening. The results show that the triaxial constraint in the crack-tip region is independent of specimen geometry, and material strain hardening, both for a stationary and an extending crack quasi-statically. The triaxial constraints for the various configurations examined are in better accordance with those required by the HRR solution for a stationary crack, which defines the low and similar constraints in crack-tip regions for different material strain hardening in the plane stress case. Along the entire ligament ahead of a crack tip, the constraint level transites gradually from that defined by the HRR solution within the near tip zone to that characterized by the stress intensity factor K _I in the far field.     

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.     

Plane strain crack closure and cyclic hardening

This study is focused on the understanding of the mechanical effects of cyclic hardening on crack tip plasticity and on plasticity-induced crack closure. Various finite element analyses were conducted using abaqus. Cyclic hardening is found to affect both crack closure and the shape of the plastic zone at the crack tip. Crack growth modelling in plane strain conditions in a cyclically hardening material is discussed. An empirical formula is provided which allows the calculation of the crack tip plastic zone size under plane strain conditions in a cyclically hardening material. The effects of overloads are also examined.     

A numerical investigation of ductile fracture initiation in a high-strength low-alloy steel

In this work, static and drop-weight impact experiments, which have been conducted using three-point bend fracture specimens of a high-strength low-alloy steel, are analysed by performing finite-element simulations. The Gurson constitutive model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence is employed within the framework of a finite deformation plasticity theory. Two populations of second-phase particles are considered, including large inclusions which initiate voids at an early stage and small particles which require large strains to nucleate voids. The most important objective of the work is to assess quantitatively the effects of material inertia, strain rate sensitivity and local adiabatic temperature rise (due to conversion of plastic work into heat) on dynamic ductile crack initiation. This is accomplished by comparing the evolution histories of void volume fraction near the notch tip in the static analysis with the dynamic analyses. The results indicate that increased strain hardening caused by strain rate sensitivity, which becomes important under dynamic loading, plays a benign role in considerably slowing down the void growth rate near the notch tip. This is partially opposed by thermal softening caused by adiabatic heating near the notch tip.     

Dynamic crack growth along an elastoplastic bimaterial interface

Summary The near-tip stress and deformation rate fields of a crack dynamically propagating along an interface between dissimilar elastic-plastic bimaterials are presented in this paper. The elastic-plastic materials are characterised by theJ ₂-flow theory with linear plastic hardening. The solutions are assumed to be of variable-separable form with a power-law singularity in the radial direction. Two distinct solutions corresponding to the tensile and shear solutions exist with slightly different singularity strengths and very different mixities at the crack tip. The phenomenon of discrete and determinate mixities at the interfacial crack tip is confirmed in dynamic crack growth. This is not an artifact of the variable-separable solution assumption, arising from the linear-hardening material model. The dynamic crack analysis shows that the mixity of the near-tip field is mainly determined by the given material parameters and affected slightly by the crack propagation velocity. A significant variation of the mixity is observed near to the coalescing point of the tensile and shear solutions. The strength of the singularity is almost determined by the smaller strain-hardening alone, and dynamic inertia decreases the stress intensity. The asymptotic solutions reveal that the crack propagation velocity changes only the stress field of the tensile mode significantly. With increasing the crack propagation velocity, the stress singularity of the tensile solutions decreases obviously and the stress triaxiality at the tip (=0) falls considerably at the unity effective stress. These observations imply that the fracture toughness of the interface crack under tensile mode may be significantly higher than that under quasi-static conditions.     

Effeet of Hardening on Elastic-plastic Curved Crack Tip Displacement Under Dynamic Load

In this article,elastic-plastic curved crack tip opening displacement maximum in hardened material under dynamic load has been studied,and curved crack tip opening displacement has been calculated as a practical application of a second order perturbation method and theorem of surname KA,where the effects of dynamic applied stresses and dynamic normal and shear stresses on the boundaries of plastic area are synthetically taken into considerations. Diagrams have been constructed to analyz the transformation relationships between the curved crack tip opening displacement and the work hardening exponents. Curved crack tip opening displacement will decrease with the increasing of hardening exponents in hardened material. The decrease extent of curved crack tip dynamic opening displacement will be more and more severe when hardening exponents increase evenly. Maximum dynamic opening displacement of curved crack tip will decrease when external load decrease with the same hardening exponents.     

Dynamic asymptotic mode III crack tip field in rate dependent materials

The asymptotic field at a dynamically growing crack tip in strain-rate sensitive elastic-plastic materials is investigated under anti-plane shear loading conditions. In the conventional viscoplasticity theory, the rate sensitivity is included only in the flow stress. However, it is often found that the yield strength is also affected by previous strain rates. The strain rate history effects in metallic solids are observed in strain rate change tests in which the flow stress decreases gradually after a rapid drop in strain rate. This material behavior may be explained by introducing the rate sensitivity in the hardening rule in addition to the flow rule. The strain-rate history effect is pronounced near the propagating crack where the change of strain rates take place. Effects of the rate dependency in the flow rule and the hardening rule on the crack propagation are analyzed. The order of the stress singularity in the asymptotic field is determined in terms of material parameters which characterize the rate sensitivity of the material. The results show that an elastic sector is present in the wake zone when the rate-dependency is considered only in the hardening rule. Terminal crack propagation speed is determined by applying the critical stress fracture criterion and the critical strain criterion to the asymptotic fields under the small scale yielding condition.     

Ratchetting strain as a damage parameter in controlling crack growth at elevated temperature

Progressive increase in tensile strains near a crack tip has been observed from finite element studies of stationary and growing cracks (Zhao, 2004, 2008) [1] and [2] under cyclic loading conditions. In this work, the significance of such a phenomenon was further explored. In particular, stress-controlled experiments were carried out to evaluate the uniaxial ratchetting response of a nickel-based superalloy, and the material parameters were re-calibrated using both strain-controlled and stress-controlled experimental data. An additional kinematic hardening term was introduced in the viscoplastic constitutive model and the models were utilised via a user-defined subroutine to study near crack tip ratchetting behaviour of a single edge notch tension (SENT) model geometry at elevated temperature. Loading modes near the crack tip were examined, together with the influence of particular constitutive models on the mechanistic response of the crack tip. The crack tip deformation was found to be predominantly strain-controlled, where the mean ratchetting strain seems to be more relevant to crack growth than the strain range. The former was used as a measure of crack tip damage to correlate crack growth rates at selected loading conditions.     

Studies on the fracture mechanics parameters of weldment with mechanical heterogeneity

The fully plastic solutions of welded centre-cracked strip for plane stress problem were carefully investigated with the fully plastic finite element method. It was introduced for assessing the fracture mechanics parameters of weldment with mechanical heterogeneity that there existed an equivalent yielding stress and equivalent strain hardening exponent in the vicinity of crack tip keeping the assessment of fracture mechanics parameters of weldment in the same way as the homogeneous material. The equivalent yielding stress and equivalent strain hardening exponent of various matched weldment were computed and the effect of weld metal width were calculated and discussed on equivalent yielding stress and equivalent strain hardening exponent near crack tip. The engineering approach was given for estimating the fracture mechanics parameters of weldment with mechanical heterogeneity in elastic-plastic range.     

Element free galerkin study of steady quasi-static crack growth in plane strain tension in elastic-plastic materials

Employing an extension of the convective mesh technique, an Element Free Galerkin (EFG) based formulation for steady quasi-static crack growth in elastic-plastic materials undergoing small scale yielding is established. In accordance with the steady state condition, a parallel path constitutive law integration scheme is implemented into the formulation and a body force type term is introduced for the contribution of crack tip plasticity. A detailed numerical analysis is performed for steady quasi-static growth of a mode I crack under plane strain conditions in elastic-perfectly plastic materials. The numerical solutions to this problem confirm the existence of an elastic unloading wedge. The computed near-tip stress distributions and the crack opening displacements for this problem provide good agreements with the corresponding asymptotic solutions and demonstrate the validity of the EFG based formulation. A rough estimate of the range of validity of the asymptotic solutions is also made based on the numerical solutions. In addition, the effect of strain hardening on the steady quasi-static crack growth is investigated. Both power hardening and linear hardening models are considered.     

Finite deformation cohesive polygonal finite elements for modeling pervasive fracture

In the present study, mode I crack subjected to cyclic loading has been investigated for plastically compressible hardening and hardening–softening–hardening solids using the crack tip blunting model where we assume that the crack tip blunts during the maximum load and re-sharpening of the crack tip takes place under minimum load. Plane strain and small scale yielding conditions have been assumed for analysis. The influence of cyclic stress intensity factor range (\(\Delta \hbox {K})\), load ratio (R), number of cycles (N), plastic compressibility (\({\upalpha })\) and material softening on near tip deformation, stress–strain fields were studied. The present numerical calculations show that the crack tip opening displacement (CTOD), convergence of the cyclic trajectories of CTOD to stable self-similar loops, plastic crack growth, plastic zone shape and size, contours of accumulated plastic strain and hydrostatic stress distribution near the crack tip depend significantly on \(\Delta \hbox {K}\), R, N, \({\upalpha }\) and material softening. For both hardening and hardening–softening–hardening materials, yielding occurs during both loading and unloading phases, and resharpening of the crack tip during the unloading phase of the loading cycle is very significant. The similarities are revealed between computed near tip stress–strain variables and the experimental trends of the fatigue crack growth rate. There was no crack closure during unloading for any of the load cycles considered in the present study.     

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.     

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.     

Energy density approach for calculation of three-dimensional inelastic strain and stress at the crack tip in compact tension specimens of polycarbonate

An energy-based method is utilized for calculating elastic-plastic strains and stresses near fatigue crack tip in specimens of Merlon polycarbonate. The stress redistribution caused by the plastic yielding around the crack tip is taken into account so that theoretical crack tip strain is improved. The estimated values of crack tip strain based on an energy density approach are compared with experimental results obtained from an embedded grid moire technique and embedded strain gages. Large-scale yielding seems to dominate near the crack tip. In fact, the measured strain is in agreement with the elastic solution, which means, in reality, only small-scale yielding takes place near the crack tip. The strain in the mid-plane (plane strain) is found to be higher than in the surface plane (plane stress). The experimental and theoretical results are in good agreement.     

A finite element analysis of plane strain dynamic crack growth at a ductile-brittle interface

In this work, steady, dynamic crack growth under plane strain, small-scale yielding conditions along a ductile-brittle interface is analysed using a finite element procedure. The ductile solid is taken to obey the J ₂ flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behaviour. The objectives of this work are to establish the validity of an asymptotic solution for this problem which has been derived recently [12], and to examine the effect of changing the remote (elastic) mode-mixity on the near-tip fields. Also, the influence of crack speed on the stress fields and crack opening profiles near the propagating interface crack tip is assessed for various bi-material combinations. Finally, theoretical predictions are made for the variation of the dynamic fracture toughness with crack speed for crack growth under a predominantly tensile mode along ductile-brittle interfaces. Attention is focused on the effect of mismatch in stiffness and density of the constituent phases on the above aspects.     

The effect of T-stress on plane strain dynamic crack growth in elastic–plastic materials

In this paper, the effects of T‐stress on steady, dynamic crack growth in an elastic–plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J₂ flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T‐stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T‐stress for both ductile and cleavage mode of crack growth.     

Interfacial crack-tip constraints and <Emphasis Type="Italic">J</Emphasis>-integral for bi-materials with plastic hardening mismatch

This paper investigates interfacial crack tip stress fields and the J-integral for bi-materials with plastic hardening mismatch via detailed elastic-plastic finite element analyses. For small scale yielding, the modified boundary layer formulation with the elastic T-stress is employed. For fully plastic yielding, plane strain single-edge- cracked specimens under pure bending are considered. Interfacial crack tip stress fields are explained by modified Prandtl slip-line fields. It is found that, for bi-materials consisting of two elastic-plastic materials, increasing plastic hardening mismatch increases both crack-tip stress constraint in the lower hardening material and the J-contribution there. The implication of asymmetric J-integral in bi-materials is also discussed.     

The evolution of the stress–strain fields near a fatigue crack tip and plasticity-induced crack closure revisited

The evolution of the stress–strain fields near a stationary crack tip under cyclic loading at selected R‐ratios has been studied in a detailed elastic–plastic finite element analysis. The material behaviour was described by a full constitutive model of cyclic plasticity with both kinematic and isotropic hardening variables. Whilst the stress/strain range remains mostly constant during the cyclic loading and scales with the external load range, progressive accumulation of tensile strain occurs, particularly at high R‐ratios. These results may be of significance for the characterization of crack growth, particularly near the fatigue threshold. Elastic–plastic finite element simulations of advancing fatigue cracks were carried out under plane‐stress, plane‐strain and generalized plane‐strain conditions in a compact tension specimen. Physical contact of the crack flanks was observed in plane stress but not in the plane‐strain and generalized plane‐strain conditions. The lack of crack closure in plane strain was found to be independent of the material studied. Significant crack closure was observed under plane‐stress conditions, where a displacement method was used to obtain the actual stress intensity variation during a loading cycle in the presence of crack closure. The results reveal no direct correlation between the attenuation in the stress intensity factor range estimated by the conventional compliance method and that determined by the displacement method. This finding seems to cast some doubts on the validity of the current practice in crack‐closure measurement, and indeed on the role of plasticity‐induced crack closure in the reduction of the applied stress intensity factor range.     

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.