Dugdale plastic zone of a penny-shaped crack in a piezoelectric material under axisymmetric loading

The Dugdale plastic zone ahead of a penny-shaped crack in a piezoelectric material, subjected to electric and axisymmetric mechanical loadings, is evaluated analytically. Hankel transform is employed to reduce the mixed boundary-value problem of the penny-shaped crack to dual integral equations, which are solved exactly under the assumption of electrically permeable crack face conditions. A closed-form solution to the mixed boundary-value problem is obtained to predict the relationship between the length of the plastic zone and the applied loading. The stress distribution in and outside of the yield zone has been derived analytically, and the crack opening displacement has been investigated. The electric displacement has a constant value in the strip yield zone. The current Dugdale crack model leads to non-singular stress and electric fields near the crack front, and it is observed that the material properties affect the crack opening displacement.     

Plastic deformation in a thick transversely isotropic layer containing a penny-shaped crack

Summary The width of a thin plastic annular zone formed during the deformation of a pennyshaped crack in a transversely isotropic layer of an ideal elasto-plastic material is determined. Considered are the cases where the penny-shaped crack is extended by normal stresses and by torsional stresses. The faces of the layer are shear free and deformation of the plastic zone around the penny-shaped crack occurs according to the Dugdale hypothesis. For each case, the solution of the problem is reduced to a Fredholm integral equation of the second kind. Iterative solutions are obtained for small values of the parameters and numerical results for the width of the plastic zone are determined. Graphical results showing the effect of transverse isotropy upon the width of the plastic zone are also presented.With 6 Figures     

Dugdale model solutions for a penny-shaped crack in three-dimensional transversely isotropic piezoelectric media by boundary-integral equation method

Making use of the Displacement Discontinuity Boundary Integral Equation Method (DDBIEM), the dimension of the plastic zone at the tip of a penny-shaped crack in a three-dimensional elastic medium is determined by the application of the Dugdale model; Furthermore, the solutions for a penny-shaped crack in three-dimensional piezoelectric media are obtained by the use of the Dugdale-like model proposed by Gao et al.[Gao H, Zhang T, Tong P. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 1997;45:491–510], in which the electrical polarization is assumed to reach a saturation limit in a thin annular region in front of a crack while the mechanical stresses have the ordinary singularity.     

Ductile penny-shaped crack in a transversely isotropic cylinder

The problem of a penny-shaped crack contained in a transversely isotropic cylinder of elastic perfectly-plastic material is considered for the case when the crack is extended by an axial load. The problem is reduced to solving numerically a Fredholm integral equation of the second kind for the width of the plastic zone. Graphical results are presented showing the effect of transverse isotropy upon the width of the plastic zone and these are compared with the results for isotropic materials.     

Effects of notch radius and anisotropy on the crack tip plastic zone

Factors which influence the shape and size of the plastic zone in the immediate vicinity of a crack tip in isotropic materials at small loads are investigated. The plastic zone dimensions for the opening mode (Mode I) have been calculated over a range of values for the crack tip radius. An increase in tip radius results in an increase in the plastic zone dimension. In anisotropic materials, the orientation of crack slit and the anisotropic yield constants are other factors that affect the plastic zone size and shape. In this paper, typical curves for the shape and size of plastic zone are given to illustrate the influence of normal or shear anisotropic yield constants. For sheet metals the effects of anisotropy on the plastic zone dimensions can be evaluated in terms of R values. Suggested values of constant b for isotropic materials are given if the “radius” approximation is employed for small applied stresses.     

Elastic–plastic analysis of interaction between an interfacial crack and a subinterfacial microcrack in bi-materials

In the present work the finite element method is used to analyze the effect of interaction between an interfacial crack and a microcrack in ceramic/aluminum bi-materials. The behaviour is analyzed by the determination of the J integral, the plastic zone at the tips of the interfacial crack and the microcrack. The effects of longitudinal and transversal distance between the tips of the two cracks and the rotation of the microcrack are examined. The obtained results allow us to deduce a correlation between the position of the microcrack and the J integral and the plastic zone.The obtained results shows that the J integral at the interfacial crack tip reaches a maximum value when the microcrack is moved in the vicinity of the interfacial crack. With this distance the effect of interaction is marked more; the stress field at the microcrack tip and that of the interfacial crack generates only one plastic zone at the interfacial crack tip. The maximum size of the plastic zone is localised at the interfacial crack tip. Those of the two tips of the microcrack are very weak and even negligible in front of the zone plasticized at the interfacial crack tip.     

The interface penny-shaped crack reconsidered

The penny-shaped crack at the interface between two bonded dissimilar media is reconsidered on the basis of recent developments on the elimination of oscillatory singularities. This is accomplished by assuming an annular frictionless contact zone at the crack circumference and reducing the problem to a Fredholm integral equation. Expressions for the strain energy, crack opening force and bond stresses are obtained and numerical results given for specific material combinations.     

Fatigue Crack Propagation Behavior of TiNi50.6Shape Memory Alloy

The fatigue crack propagation behavior of TiNi50.6 shape memory alloy was studied. The experiment results showed that the crack propagation properties of this alloy display difference and similarity in comparison with common metallic materials. Because of the stress concentration there was stress induced martensite transformation (SIMT) near the crack tip though thenominal stress was lower than the threshold stress of SIMT. The position and the amount of SIMT was in situ observed by a quester remote measurement system (QRMS). The observationresults showed that the position of SIMT was beside the crack tip and was not in the plastic zone of common     

Plastic zone size estimation under cyclic loadings using in situ optical microscopy fatigue testing

An in situ optical microscopy fatigue testing is proposed in this paper to investigate the forward and reversed plastic zone size under cyclic loadings for Al‐7075‐T6. This experimental study is used to verify the hypotheses in a recently developed small time scale formulation of fatigue crack growth. During the testing, the entire cyclic loading cycle is divided into a certain number of steps. Images of the crack tip are taken at each step. The full strain field around the crack tip is determined using the digital image correlation (DIC) technique. The plastic zone size is obtained by combining the DIC results and the material constitutive relationship. Experimental measurements from the proposed study are compared with theoretical predictions. It is observed that the crack closure has a large effect on the reversed plastic zone size. The plastic zone size remains almost constant when the unloading path is below a certain stress level, which is one of the hypotheses used in a previous crack growth model. Discussions are given for the modelling of plastic zone size variation under cyclic loadings and several conclusions are drawn based on the current investigation.     

Plastic energy dissipation due to a penny-shaped crack

International Journal of Fracture - A penny-shaped crack in a material which is ideally elastic-plastic has been envisaged with the assumption that the plastic zone forms a very thin layer...     

CRACK CLOSURE AND PLASTIC ZONE SIZES IN FATIGUE

Abstract— An elastic-plastic finite element simulation of growing fatigue cracks which accounts for plasticity-induced crack closure is used to study the size of the forward and reversed plastic zones at the crack tip. Forward plastic zone widths for fatigue cracks and stationary, monotonically loaded cracks are compared and found to be similar. The width of the forward plastic zone at the tip of a fatigue crack is not significantly influenced by closure. The traditional Irwin-Rice estimate for crack tip plastic zone size in plane stress is found to be generally consistent with the finite element results. The width of the reversed plastic zone at the tip of a growing fatigue crack in plane stress is found to be considerably less than one-fourth the size of the forward plastic zone, the traditional Rice estimate. This decrease appears to be due to fatigue crack closure. A simple model is developed which permits estimation of the reversed plastic zone size for any stress ratio from information about maximum and minimum stresses and the closure stress. The predictions of this model agree closely with plastic zone sizes calculated by the finite element analysis. These observations appear to be consistent with experimental measurements of forward and reversed plastic zones sizes reported in the literature.     

On the full set of elastic T-stress terms of internal elliptical cracks under mixed-mode loading condition

Analytical expressions for the elastic constant stress terms of the asymptotic field, the so called T-stresses, for internal mixed-mode elliptical cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed-mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses for mixed-mode elliptical cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all singular and constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

On the full set of elastic T-stress terms of internal circular cracks under mixed-mode loading conditions

Analytical expressions for all non-singular stress terms of the asymptotic crack tip field, the so-called T-stresses of internal mixed mode circular (penny shaped) cracks in infinite homogeneous and isotropic elastic solids are addressed. To solve the problem the mixed mode crack problem is divided into sub-problems using the superposition method, and each sub-problem is then solved for the asymptotic stress field. Considering the expansion of the local stress field at the crack front, the elastic T-stress terms are derived for each sub-problem. The results are superimposed to give the analytical expressions of the so far missing elastic T-stresses of internal mixed mode penny shaped cracks.The effect of the T-stresses on the size and shape of the plastic zone at the crack tip is discussed, and analytical results are compared to the ones from finite element analyses, both for the T-stress components and the size of the plastic zone. For an accurate prediction of the plastic zone all the singular terms and the constant terms (T-stresses) in the stress expansion formulae should be considered. It is observed that negative T-stresses increase the size of the plastic zone, while positive ones reduce it.     

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.     

The effects of inter-surface cohesive tractions on linear and penny-shaped cracks

A mesoscopic fracture model of equilibrium slit cracks in brittle solids, including inter-surface cohesive tractions acting near the crack tip, is analyzed and the effects of the cohesive tractions on the in-plane stress fields, crack-opening displacement profiles, and crack driving forces examined quantitatively for linear and penny-shaped cracks. The (numerical) analysis method is described in detail, along with results for four different cohesive forces. The assumed distribution of cohesive tractions were found to suppress the in-plane stress field adjacent to cracks in a homogeneous, isotropic medium when uniformly loaded in mode-I, and the suppression was a function of crack length. The crack-opening displacement profile was also perturbed and a new regime identified between the near-field Barenblatt zone and the far-field continuum zone. The extent of this `cohesive zone' was quantified by use of an interpolating function fit to the calculated profiles and found to be independent of crack size for a given cohesive tractions distribution. Furthermore, the crack-opening displacement at the edge of the cohesive zone was found to be independent of crack size, implying that despite significant perturbations to the stress field, the crack driving force at unstable equilibrium remains unchanged with crack size.     

Interface fracture in adhesively bonded shell structures

Two methods for the prediction of crack propagation through the interface of adhesively bonded shells are discussed. One is based on a fracture mechanics approach; the other is based on a cohesive zone approach. Attention is focussed on predicting the shape of the crack front and the critical stress required to propagate the crack under quasi-static conditions. The fracture mechanical model is theoretically sound and it is accurate and numerically stable. The cohesive zone model has some advantages over the fracture mechanics based model. It is easier to generalise the cohesive zone model to take into account effects such as plastic deformation in the adhering shells, and to take into account effects of large local curvatures of the interface crack front. The comparison shows a convergence of the results based on the cohesive zone model towards the results based on a fracture mechanics approach in the limit where the size of the cohesive zone becomes smaller than other relevant geometrical lengths for the problem. However, convergence issues and numerical stability must be addressed.     

Stress intensity factors for penny-shaped cracks perpendicular to graded interfacial zone of bonded bi-materials

This paper examines the stress intensity factors that are associated with a penny-shaped crack perpendicular to the interface of a bi-material bonded with a graded interfacial zone. Elastic modulus of the graded interfacial zone is assumed to be an exponential function of the depth. The stress intensity factors are calculated numerically using a so-called generalized Kelvin solution based boundary element method. Three cases of normal or shear tractions acting on the crack surfaces are examined. Values of the stress intensity factors are examined by taking into account the effects of the following four parameters: (a) the crack front position; (b) the non-homogeneity parameter of the graded interfacial zone; (c) the crack distance to the graded interfacial zone; and (d) the graded interfacial zone thickness. The numerical results are compared well with existing solutions under some degenerated conditions. These results are useful to furthering our knowledge on fracture behavior of bi-material systems with or without a graded interfacial zone.     

Delamination vs matrix cracking in layered composites subjected to transverse loading

An algorithm is presented to compute the distribution of the strain energy release rate along the crack front of a penny-shaped delamination in a layered orthotropic body. The method applies a finite element recently proposed for three-dimensional analysis of layered orthotropic circular plates. The algorithm is economical even though it treats a full three-dimensional state of stress. The method requires only a single virtual crack extension to accurately compute the strain energy release rate at a point along the crack front. The method is applied to the study of delamination crack growth in a nine layer cross-ply laminate. The variation of strain energy release rate, G, along the crack front, is determined. The significance of the plate aspect ratio, as well as length scale, on the fracture process is studied. The establishment of a loading case where a distributed transverse compressive loading causes delamination growth is given.     

The indentation of a precompressed penny-shaped crack

The paper examines the axisymmetric problem related to the indentation of the plane surface of a penny-shaped crack by a smooth rigid disc inclusion. The crack is also subjected to a far-field compressive stress field which induces closure over a part of the crack. The paper presents the Hankel integral transform development of the governing mixed boundary value problem and its reduction to a single Fredholm integral equation of the second kind and an appropriate consistency condition which considers the stress state at the boundary of the crack closure zone. A numerical solution of this integral equation is used to develop results for the axial stiffness of the inclusion and for the stress intensity factors at the tip of the penny-shaped crack.     

Steady crack growth in a thin, ductile plate under small-scale yielding conditions: Three-dimensional modeling

This work employs high resolution, finite element computations to investigate key features of the elastic–plastic fields near a steadily advancing crack at quasi-static rates under three-dimensional, small-scale yielding conditions. The model represents a structurally thin component constructed of a material (e.g., Al and Ti alloys) with flow stress and fracture toughness properties that together limit the size of the in-plane plastic zone during steady-growth to no more than several multiples of the plate thickness. The computational approach generalizes the streamline integration procedure used previously for two-dimensional studies into three dimensions to represent steady-state growth on a fixed mesh in a boundary-layer framework. The plate thickness provides the only geometrical length scale. Crack extension occurs at the remotely applied, fixed loading without the need for a local growth criterion. In the first computations of this type, the present work considers a straight crack front advancing under local and global mode I loading with zero T-stress in a moderately hardening material. Applied remote loads at steady growth generate plastic zone sizes ahead of the advancing crack front ranging from 0.25 to 6.4 times the thickness. Key results include: (1) the crack-front fields exhibit a self-similar scaling characterized by a non-dimensional loading parameter; (2) three-dimensional effects extend to distances of approximately 1.5–2.5 times the thickness ahead of the advancing crack front for key values of this loading parameter, beyond which the fields (elastic–plastic then linear-elastic at greater distances) become uniform over the thickness; and (3) crack opening profiles on the outside surface reveal a “wedge-like”, opening shape which simplifies the definition of a crack-tip opening angle.