Anti-plane shear fields with discontinuous deformation gradients near the tip of a crack in finite elastostatics

Summary This paper reconsiders the problem of determining the elastostatic field near the tip of a crack in an all-round infinite body deformed by a Mode III loading at infinity to a state of anti-plane shear. The problem is treated for a class of incompressible, homogeneous, isotropic elastic materials whose constitutive laws permit a loss of ellipticity in the governing displacement equation of equilibrium at sufficiently severe shearing strains. The analysis represents a generalization of that reported in an earlier study and, as before, is carried out for the small-scale nonlinear crack problem, in which a crack of finite length is replaced by a semi-infinite one, and the nonlinear field far from the crack-tip is matched to the near field predicted by the linearized theory. The methods employed in the present paper are necessarily largely qualitative, since they apply to all materials in the class considered. The principal feature of the resulting elastic field is the presence of two symmetrically located curves issuing from the crack-tip and bearing discontinuities in displacement gradient and stress.The results communicated in this paper were obtained in the course of an investigation supported in part by Contract N00014-75-C-0196 with the Office of Naval Research in Washington, D.C.     

Anti-plane shear cracks in ideally plastic crystals

Cracks in ductile single crystals are analyzed here for geometries and orientations such that two-dimensional states of anti-plane shear constitute possible deformation fields. The crystals are modelled as ideally plastic and yield according to critical resolved shear stresses on their slip systems. Restrictions on the asymptotic forms of stress and deformation fields at crack tips are established for anti-plane loading of stationary and quasistatically growing cracks, and solutions are presented for several specific orientations in f.c.c. and b.c.c. crystals. The asymptotic solutions are complemented by complete elastic-plastic solutions for stationary and growing cracks under small scale yielding, based on previous work by Rice (1967, 1984) and Freund (1979). Remarkably, the plastic zone at a stationary crack tip collapses into discrete planes of displacement and stress discontinuity emanating from the tip; plastic flow consists of concentrated shear on the displacement discontinuities. For the growing crack these same planes, if not coincident with the crack plane, constitute collapsed plastic zones in which velocity and plastic strain discontinuities occur, but across which the stresses and anti-plane displacement are fully continuous. The planes of discontinuity are in several cases coincident with crystal slip planes but it is shown that this need not be the case, e.g., for orientations in which anti-plane yielding occurs by multi-slip, or for special orientations in which the crack tip and the discontinuity planes are perpendicular to the activated slip plane.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. 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 crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

Singular fields of a mode III interface crack in a power-law hardening bimaterial

A high order of asymptotic solution of the singular fields near the tip of a mode III interface crack for pure power-law hardening bimaterials is obtained by using the hodograph transformation. It is found that the zero order of the asymptotic solution corresponds to the assumption of a rigid substrate at the interface, and the first order of it is deduced in order to satisfy completely two continuity conditions of the stress and displacement across the interface in the asymptotic sense. The singularities of stress and strain of the zeroth order asymptotic solutions are −1/(n ₁+1) and −n/(n ₁+1) respectively. (n=n ₁,n ₂ is the hardening exponent of the bimaterials.) The applicability conditions of the asymptotic solutions are determined for both zeroth and first orders. It is proved that the Guo-Keer solution^[10] is limited in some conditions. The angular functions of the singular fields for this interface crack problem are first expressed by closed form. The project supported by National Natural Science Foundation of China     

Analysis of deformation-induced crack tip toughening in ductile single crystals by nano-indentation

This paper reports on the experimental examination of the deformation characteristics near a crack tip in a cyclically work-hardened copper single crystal using a 2D surface scans with nano-indentation. The experimental methodology enables the characterization of the primary deformation field near a crack tip via the modulation of the imposed secondary deformation field by a nano-indentation. In a heavily deformed 4-point bend specimen, the measurements showed an existence of an asymptotic field around the crack tip at a distance of R ⩾ 2.5J/σ₀. The measurements also showed the qualitative details of toughness evolution within the high-gradient deformation field around the crack tip. The nature of dislocation distribution (i.e. statistically distributed vs. distributions necessitated by geometry) around the crack tip is quantified. The measurements indicate the dominance of the geometrically necessary dislocation within the finite deformation zone ahead of the tip up to a distance of R ≈ 3J/σ₀. Thereafter, it is confined in radial rays coinciding with the sector boundaries around the crack tip. These measurements elucidate the origin of the inhomogeneous hardening and the size dependent macroscopic response close to crack tip.     

The near tip stress and strain fields in cracked single crystals

The numerical analyses of stationary mathematically sharp Mode I crack in FCC and BCC crystals with elastic-ideally plastic (EIP) and fast hardening saturation (FHS) law are carried out in the present paper. From the calculated results, it is shown that: for the cases of small strain, EIP crystal cracks, the features of concentrated deformation patterns and the stress state in near-crack tip deformation fields are identical to the earlier analytical solutions, but along the angular sector boundaries, there exist narrow complex stress zones. The overall characteristics of deformation patterns for the cases of EIP and FHS are similar. The behaviours of crack tip opening can be characterized by crack-tip-opening-displacement (CTOD). For the case of FHS, finite deformation BCC crystal crack, our calculations are qualitatively in agreement with recent experimental observations. The project supported by National Natural Science Foundation of China     

Higher-order crack tip fields for functionally graded material plate with transverse shear deformation

The crack tip fields are investigated for a cracked functionally graded material (FGM) plate by Reissner’s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson’s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner’s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner’s effect when the in-homogeneity parameter approaches zero.     

A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals

A phase-field theory of dislocation dynamics, strain hardening and hysteresis in ductile single crystals is developed. The theory accounts for: an arbitrary number and arrangement of dislocation lines over a slip plane; the long-range elastic interactions between dislocation lines; the core structure of the dislocations resulting from a piecewise quadratic Peierls potential; the interaction between the dislocations and an applied resolved shear stress field; and the irreversible interactions with short-range obstacles and lattice friction, resulting in hardening, path dependency and hysteresis. A chief advantage of the present theory is that it is analytically tractable, in the sense that the complexity of the calculations may be reduced, with the aid of closed form analytical solutions, to the determination of the value of the phase field at point-obstacle sites. In particular, no numerical grid is required in calculations. The phase-field representation enables complex geometrical and topological transitions in the dislocation ensemble, including dislocation loop nucleation, bow-out, pinching, and the formation of Orowan loops. The theory also permits the consideration of obstacles of varying strengths and dislocation line-energy anisotropy. The theory predicts a range of behaviors which are in qualitative agreement with observation, including: hardening and dislocation multiplication in single slip under monotonic loading; the Bauschinger effect under reverse loading; the fading memory effect, whereby reverse yielding gradually eliminates the influence of previous loading; the evolution of the dislocation density under cycling loading, leading to characteristic ‘butterfly’ curves; and others.     

Higher order fields at crack and notch tips in power-law materials under longitudinal shear

Summary Singular fields and higher order fields near a sharp notch in a power-law material under longitudinal shear are investigated. Using the perturbation method the whole set of eigenvalues is determined. The higher order eigensolutions are constructed by use of the dominant singular solution. Some examples, including the special case of a crack are discussed for different boundary conditions.

---

Felder höherer Ordnung an Riß- und Kerbspitzen unter nichtebener Schubbelastung in Materialien mit Potenzverfestigungsgesetz

Übersicht Untersucht werden singuläre Felder und Felder höherer, Ordnung an Spitzkerben unter nichtebener Schubbelastung in Materialien mit Potenzgesetz. Die Eigenwerte werden mittels der Störungsrechnung ermittelt. Die Bestimmung der Eigenlösungen höherer Ordnung erfolgt unter Verwendung der dominierenden singulären Lösung. Einige Beispiele, die auch den Spezialfall des Risses einschließen, werden für unterschiedliche Randbedingungen diskutiert.

---

---

Dedicated to Prof. Dr.-Ing. F. G. Kollmann on the occasion of his sixtieth Birthday     

Theoretical latent hardening of crystals in double slip—I. F.C.C. crystals with a common slip plane

The hardening of all slip systems in f.c.c. crystals, deforming in finite double-slip on two systems with a common slip plane, is determined according to the “simple theory” of rotation-dependent plastic anisotropy. Both tensile and compressive axial-loading are considered. Of particular interest are predictions of crystal response after the loading axis has rotated into a higher symmetry position (6-fold in tension and 4-fold in compression). In contrast with classical Taylor hardening, the simple theory predicts that the axis will “overshoot” the higher symmetry position. A postulate of minimum plastic work plays a significant role in the theoretical analyses of multiple-slip positions.     

Theoretical latent hardening in crystals—I. General equations for tension and compression with application to F.C.C. crystals in tension

Equations for latent strengths in single slip, based upon the simple theory of finite distortional crystal hardening introduced by K.S. Havner and A.H. Shalaby (1977), are derived for both tensile and compression tests without restriction as to crystal class. Detailed comparisons between theoretical results and the experiments of P.J. Jackson and Z.S. Basinski (1967) on copper crystals in tension are presented. There is good qualitative agreement between theory and experiment regarding the diversity of anisotropic hardening among slip systems. Moreover, there is satisfactory quantitative agreement between the theory and the extrapolated experimental data in the stage III, large-strain range. It is suggested that further experimental investigation of latent hardening at large prestrains would be desirable.The simple theory predicts anisotropic hardening and the perpetuation of single slip in axial loading of cubic crystals initially oriented for single slip, but predicts symmetric, isotropic hardening of specimens initially oriented in positions of 4, 6 or 8-fold multiple-slip. These predictions are in general accord with experimental observations from tests of f.c.c. and b.c.c. crystals.     

The singular field near a sharp notch with mixed boundary conditions for hardening materials under longitudinal shear

The field behavior near a sharp notch tip with mixed homogeneous stress and displacement boundary conditions is examined for a power law hardening material. Using the hodograph transformation, the singularity and the angular distribution of the fields are determined. Special cases as those for linear elastic and perfect plastic materials are discussed.     

Anti-plane shear waves in a fibre-reinforced composite with a non-linear imperfect interface

The propagation of non-linear elastic anti-plane shear waves in a unidirectional fibre-reinforced composite material is studied. A model of structural non-linearity is considered, for which the non-linear behaviour of the composite solid is caused by imperfect bonding at the “fibre–matrix” interface. A macroscopic wave equation accounting for the effects of non-linearity and dispersion is derived using the higher-order asymptotic homogenisation method. Explicit analytical solutions for stationary non-linear strain waves are obtained. This type of non-linearity has a crucial influence on the wave propagation mode: for soft non-linearity, localised shock (kink) waves are developed, while for hard non-linearity localised bell-shaped waves appear. Numerical results are presented and the areas of practical applicability of linear and non-linear, long- and short-wave approaches are discussed.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Theoretical latent hardening in crystals—II. BCC crystals in tension and compression

The general latent hardening law of single slip derived in the first paper of this series (Havner, Baker and Vause, 1979) is applied to an analysis of “overshooting” phenomena in bcc crystals in tension and compression. This new law, which predicts anisotropic hardening of latent slip systems, is based upon the simple theory of finite distortional crystal hardening introduced by Havner and Shalaby (1977).Because of historical ambiguities regarding identification of the slip plane in bcc metals, parallel analyses are presented corresponding to two separate criteria: (i) slip on {110}, {112} and {123} crystallographic planes only; and (ii) slip on the plane of maximum resolved shear stress containing a 〈111〉 direction. It is established that the new hardening law is a theory of “overshooting” in bcc crystals according to either identification of the slip plane.A qualitative comparison between theoretical results and two experimental papers on Fe crystals is included. The general difficulties in making comparisons with the experimental literature on finite distortional latent hardening are briefly discussed.     

The near tip fields and temperature distribution around a crack in a body of hardening material containing small damage

In this paper, the following conclusions are reached: The influence of damage on the stress and strain feilds can be neglected in an asymptotic sense for the solutions of damage field in a plastic solid containing small damage. The formulation of the problem is simplified with an uncoupled approach. Based on experimental results of plastic damage, most of the damage in the material are considered as small damage with the critacal damage variable ω _c ≪1. Using this approach, closed form expressions of the near tip damage fields for mode III, mode I and the temperature distribution induced by plastic dissipation in a hardening material containing damage are deduced. We point out that the temperature distribution in the process zone is strongly dependent on the damage of materials even for the small damage case. The results of the predicted value of the temperature rise near the tip region ignoring the damage effect is appreciably higher than the observed data. The main reason of this discrepancy is the presence of damage dissipation and the fact that its influence on the calculation of plastic dissipation have not been appropriately taken account of. The calculation is improved by taking into account the damage effect on the temperature rise, then theT _max value is in better accord with the experimental value. The project supported by the National Natural Science Foundation of China.     

Crack growth resistance in non-perfect plasticity: Isotropic versus kinematic hardening

The Strain Energy Density Theory is applied for analyzing energy dissipation and crack growth in the three-point bending specimen when the material behavior follows a multilinear strain-hardening stress-strain relationship. The problem is solved through the application of incremental theory of plasticity and finite element method.The rate of change of the strain energy density factor S with crack length a is verified to be governed by the relation . Results are obtained for isotropic and kinematic hardening. Moreover, the effects of loading step and specimen size are pointed out.