‘Pressurized’ Mode III Crack in an Elastic Perfectly Plastic Solid

This paper presents the solution for the elastic plastic boundary of a fully ‘pressurized’ mode III crack in an elastic perfectly plastic solid. By fully ‘pressurized’ (in analogy to an internally pressurized mode I crack) is meant that the crack faces are subjected to an externally applied (shear) stress equal to the (shear) yield stress. The dislocations in stress space method is used to find the solution. At the crack plane a plastic zone for the fully pressurized crack extends out a distance 0.273 a from a crack tip, where a is the half width of the crack. (The comparable distance found from a Bilby–Cottrell–Swinden crack is 0.414 a.) The plastic zones originating at the two crack tips merge, in a point, at the center of the crack. This revised version was published online in August 2006 with corrections to the Cover Date.     

Interface crack in periodically layered bimaterial composite

A directional crack growth prediction in a compressed homogenous elastic isotropic material under plane strain conditions is considered. The conditions at the parent crack tip are evaluated for a straight stationary crack. Remote load is a combined biaxial compressive normal stress and pure shear. Crack surfaces are assumed to be frictionless and to remain closed during the kink formation wherefore the mode I stress intensity factor K _I is vanishing. Hence the mode II stress intensity factor K _II remains as the single stress intensity variable for the kinked crack. An expression for the local mode II stress intensity factor k ₂ at the tip of a straight kink has been calculated numerically with an integral equation using the solution scheme proposed by Lo (1978) and refined by He and Hutchinson (1989). The confidence of the solution is strengthened by verifications with a boundary element method and by particular analytical solutions. The expression has been found as a function of the mode II stress intensity factor K _II of the parent crack, the direction and length of the kink, and the difference between the remote compressive normal stresses perpendicular to, and parallel with, the plane of the parent crack. Based on the expression, initial crack growth directions have been suggested. At a sufficiently high non-isotropic compressive normal stress, so that the crack remains closed, the crack is predicted to extend along a curved path that maximizes the mode II stress intensity factor k ₂. Only at an isotropic remote compressive normal stress the crack will continue straight ahead without change of the direction. Further, an analysis of the shape of the crack path has revealed that the propagation path is, according the model, required to be described by a function y=cx , where the exponent is equal to 3/2. In that case, when =3/2, predicts the analytical model a propagation path that is self-similar (i.e. the curvature c is independent of any length of a crack extension), and which can be described by a function of only the mode II stress intensity factor K _II at the parent crack tip and the difference between the remote compressive normal stress perpendicular to, and parallel with, the parent crack plane. Comparisons with curved shear cracks in brittle materials reported in literature provide limited support for the model discussed.     

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.     

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.     

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 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.     

Stresses near Crack Tips on the Boundary of Two Media in the Presence of Plastic Strips

Under the conditions of plane deformation, we study stresses in a piecewise-homogeneous isotropic body near the tip of a mode I crack appearing at the angular point on the boundary of two media. It is assumed that plastic strips (modeled by plastic slip lines) are formed on the boundary of the media. To determine the stresses, we use the Mellin integral transformation and the Wiener–Hopf method. The angular point is a stress concentrator with power singularity and, therefore, immediately after the appearance of lateral plastic strips (zones), a new plastic zone begins to develop from this point. We study the dependence of the power of singularity of stresses on the angle made by the boundary and the elastic characteristics of the media.     

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.     

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic–perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near‐tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T ‐stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T ‐stress becomes more negative irrespective of lattice orientation. Also, the near‐tip field for a range of constraint levels can be characterized by two‐parameters such as K–T or J–Q as in isotropic plastic solids.     

Quasi-static crack tip fields in rate-sensitive FCC single crystals

In this work, the effects of loading rate, material rate sensitivity and constraint level on quasi-static crack tip fields in a FCC single crystal are studied. Finite element simulations are performed within a mode I, plane strain modified boundary layer framework by prescribing the two term (K − T) elastic crack tip field as remote boundary conditions. The material is assumed to obey a rate-dependent crystal plasticity theory. The orientation of the single crystal is chosen so that the crack surface coincides with the crystallographic (010) plane and the crack front lies along [10[`1]][10\overline 1] direction. Solutions corresponding to different stress intensity rates [(K)\dot]\dot{{K}}, T-stress values and strain rate exponents m are obtained. The results show that the stress levels ahead of the crack tip increase with [(K)\dot]\dot{{K}} which is accompanied by gradual shrinking of the plastic zone size. However, the nature of the shear band patterns around the crack tip is not affected by the loading rate. Further, it is found that while positive T-stress enhances the opening and hydrostatic stress levels ahead of crack tip, they are considerably reduced with imposition of negative T-stress. Also, negative T-stress promotes formation of shear bands in the forward sector ahead of the crack tip and suppresses them behind the tip.     

Tensile crack tip fields in elastic-ideally plastic crystals

Crack tip stress and deformation fields are analyzed for tensile-loaded ideally plastic crystals. The specific cases of (0 1 0) cracks growing in the [1 0 1] direction, and (1 0 1) cracks in the [0, 1, 0] direction, are considered for both fcc and bcc crystals which flow according to the critical resolved shear stress criterion. Stationary and quasistatically growing crack fields are considered. The analysis is asymptotic in character; complete elastic-plastic solutions have not been determined. The near-tip stress state is shown to be locally constant within angular sectors that are stressed to yield levels at a stationary crack tip, and to change discontinuously from sector to sector. Near tip deformations are not uniquely determined but fields involving shear displacement discontinuities at sector boundaries are required by the derived stress state. For the growing crack both stress and displacement must be fully continuous near the tip. An asymptotic solution is given that involves angular sectors at the tip that elastically unload from, and then reload to, a plastic state. The associated near-tip velocity field then has discontinuities of slip type at borders of the elastic sectors. The rays, emanating from the crack tip, on which discontinuities occur in the two types of solutions are found to lie either parallel or perpendicular to the family of slip plane traces that are stressed to yield levels by the local stresses. In the latter case the mode of concentrated shear along a ray of discontinuity is of kink type. Some consequences of this are discussed in terms of the dislocation generation and motion necessary to allow the flow predicted macroscopically.     

A MODEL TO PREDICT THE FATIGUE LIFE OF FIBRE-REINFORCED TITANIUM MATRIX COMPOSITES UNDER CONSTANT AMPLITUDE LOADING

A model based on micro-mechanical concepts has been developed for predicting fatigue crack growth in titanium alloy matrix composites. In terms of the model, the crack system is composed of three zones: the crack, the plastic zone and the fibre. Crack tip plasticity is constrained by the fibres and remains so until certain conditions are met. The condition for crack propagation is that fibre constraint is overcome when the stress at the location of the fibre ahead of the crack tip attains a critical level required for debonding. Crack tip plasticity then increases and the crack is able to propagate round the fibre. The debonding stress is calculated using the shear lag model from values of interfacial shear strength and embedded fibre length published in the literature. If the fibres in the crack wake remain unbroken, friction stresses on the crack flanks are generated, as a result of the matrix sliding along the fibres. The friction stresses (known as the bridging effect) shield the crack tip from the remote stress, reducing the crack growth relative to that of the matrix alone. The bridging stress is calculated by adding together the friction stresses, at each fibre row bridging the crack, which are assumed to be a function of crack opening displacement and sliding distance at each row. The friction stresses at each fibre row will increase as the crack propagates further until a critical level for fibre failure is reached. Fibre failure is modelled through Weibull statistics and published experimental results. Fibre failure will reduce the bridging effect and increase the crack propagation rate. Calculated fatigue lives and crack propagation rates are compared with experimental results for three different materials (32% SCS6/Ti-15-3, 32% and 38% SCS6/Ti-6-4) subjected to mode I fatigue loading. The good agreement shown by these comparisons demonstrates the applicability of the model to predict the fatigue damage in Ti-based MMCs.     

Hybrid stress finite element analysis of bending of a plate with a through flaw

A hybrid stress finite element procedure for the solution of bending stress intensity factors of a plate with a through-the-thickness crack is presented. Reissner's sixth-order plate theory including the effects of transverse shear deformation is used. The dominant singular crack tip stress field is embedded in the crack tip singular elements and only regular polynomial functions are assumed in the far field elements. The stress intensity factors can be calculated directly from the crack tip singular stress solution functions. The effects of the plate thickness, the ratio between the crack size and the inplane dimension of the plate, and the singular element size on the stress intensity factor solution are investigated. The effects of the explicit enforcement of traction-free conditions along crack surfaces, which are the natural boundary conditions in the present hybrid stress finite element model, are also investigated. The numerical results of bending of a plate with a straight central crack compare favourably with analytical solutions. It is also found that the explicit enforcement of traction-free conditions along crack surfaces is mandatory to obtain meaningful results for the Mode I type of bending stress intensity factor.     

A two-dimensional elastic–plastic finite element analysis of friction effects on sliding crack surfaces in full or partial contact

Sliding crack surfaces are analysed, that are completely or partially in contact, using a two-dimensional plane-stress elastic–plastic finite element technique. Our in-house program was modified to account for the friction which acts between two rough mating surfaces. The analysis is applied to a cantilever beam cracked along its span through its centroidal plane. Twelve cracks with length-to-span ratios ranging from zero to 0.5 were analysed. The effect of friction was investigated by considering 0, 0.4, 0.8, 1.2 and 1.6 as values for the coefficient of friction with each crack length. The results show the influence of friction on the beam stiffness, strain energy release rate, modes of crack tip and surface displacements, and the development of plastic deformation. The present finite element outputs assist in the explanation of experimental events associated with mode II crack tip displacement data found in the literature.     

Crack tip fields in ductile crystals

Results on the asymptotic analysis of crack tip fields in elastic-plastic single crystals are presented and some preliminary results of finite element solutions for cracked solids of this type are summarized. In the cases studied, involving plane strain tensile and anti-plane shear cracks in ideally plastic f c c and b c c crystals, analyzed within conventional small displacement gradient assumptions, the asymptotic analyses reveal striking discontinuous fields at the crack tip.For the stationary crack the stress state is found to be locally uniform in each of a family of angular sectors at the crack tip, but to jump discontinuously at sector boundaries, which are also the surfaces of shear discontinuities in the displacement field. For the quasi-statically growing crack the stress state is fully continuous from one near-tip angular sector to the next, but now some of the sectors involve elastic unloading from, and reloading to, a yielded state, and shear discontinuities of the velocity field develop at sector boundaries. In an anti-plane case studied, inclusion of inertial terms for (dynamically) growing cracks restores a discontinuous stress field at the tip which moves through the material as an elastic-plastic shock wave. For high symmetry crack orientations relative to the crystal, the discontinuity surfaces are sometimes coincident with the active crystal slip planes, but as often lie perpendicular to the family of active slip planes so that the discontinuities correspond to a kinking mode of shear.The finite element studies so far attempted, simulating the ideally plastic material model in a small displacement gradient type program, appear to be consistent with the asymptotic analyses. Small scale yielding solutions confirm the expected discontinuities, within limits of mesh resolution, of displacement for a stationary crack and of velocity for quasi-static growth. Further, the discontinuities apparently extend well into the near-tip plastic zone. A finite element formulation suitable for arbitrary deformation has been used to solve for the plane strain tension of a Taylor-hardening crystal panel containing, a center crack with an initially rounded tip. This shows effects due to lattice rotation, which distinguishes the regular versus kinking shear modes of crack tip relaxation. and holds promise for exploring the mechanics of crack opening at the tip.     

An elastic-plastic finite element analysis of crack tip fields under biaxial loading conditions

A square plate containing a central crack and subjected to biaxial stresses has been studied by a finite element analysis. An elastic analysis shows that the crack opening displacement and stress of separation ahead of the crack tip are not affected by the mode of biaxial loading and therefore the stress intensity factor adequately describes the crack tip states in an elastic continuum.

An elastic-plastic analysis involving more than localized yielding at the crack tip provides different solutions of crack tip stress fields and crack face displacements for the different modes of biaxial loading.

The equi-biaxial loading mode causes the greatest separation stress but the smallest plastic shear ear and crack displacement. The shear loading system induces the maximum size of shear ear and crack displacement but the smallest value of crack tip separation stress.

     

Application of finite element analysis techniques for predicting crack propagation in lugs

The plane elastic problem corresponding to a through radial crack emanating from the internal boundary of a symmetrical lug is considered. A pin bearing pressure distribution was developed by utilizing photoelastic test data and differs considerably from the usually assumed uniform or cosine pressure distributions. The stress intensity factors at the crack tip were obtained by using recently derived quadratic isoparametric finite elements which embody the inverse square root singularity. Fatique crack growth tests of 17 aluminium, titanium and steel lugs were utilized to verify stress intensity factor solutions.     

A crack in a linear-elastic material under mode II loading, revisited

A sharp crack in a two-dimensional infinite linear-elastic material, under pure shear (mode II) loading is re-examined. Several criteria have been proposed for the prediction of the onset and direction of crack extension along a path emanating from the tip of the initial crack. These criteria date back some three decades and are well documented in the literature. All the predictions from the different criteria are close and indicate that the crack extension takes a direction at an angle of ≈ −70° measured counterclockwise from the positive x -axis, in the case of a remotely applied positive shear stress. However, the possibility seems to have been overlooked that the crack extension may initiate not from the crack tip itself, but instead may initiate on the free surface at an infinitesimal distance behind the crack tip. The effect of crack tip plasticity on the relevant stresses in the region of the crack tip is investigated by the application of an elastic–plastic finite element program.     

Mode II stress intensity factors for central slant cracks with frictional surfaces in uniaxially compressed plates

The effect of crack surface friction on mode II stress intensity factor (SIF) of a central slant crack in a plate uniformly loaded in uniaxial compression is quantified. A previously developed two-dimensional finite element analysis was utilised after its modification to accommodate the friction between the crack surfaces. The plane strain state was assumed. A new numerical technique was devised to avoid the iteration procedures, which had to be employed due to the existence of frictional forces.

The crack inclination angle varied between zero and 75° measured from the horizontal direction. The coefficient of friction of the crack surfaces changed from zero to 1. In case of relatively sliding crack surfaces, mode II SIF existed. As is well known, the resulting mode II SIF decreased with increasing the coefficient of friction of the crack surfaces. Further, mode II SIF increased with increasing crack line inclination angle and then decreased after reaching a maximum value. The angle corresponding to that maximum SIF increased as the coefficient of friction of the crack surfaces increased.