Plastic yielding on inclined slip-planes at a crack tip

Plastic yield at crack tips on singular slip-planes, inclined to the crack plane, has been studied under plane-strain conditions for combined tension, hydrostatic stress, and in-plane shear. The singular integral equation, which represents the equilibrium condition of edge dislocations on the slip-planes, is transformed into a Fredholm integral equation in order to avoid difficulties that occur with its numerical solution. Results are presented for the slip-band length, the plastic crack-tip opening displacement, stress fields, and crack-opening contours. A series expansion of the results obtained numerically confirms approximate analytical expressions given by J.R. Rice (1974), up to the third-order in the applied stresses. The results of finite element methods agree with values of the crack-tip opening displacement obtained here to within 10 per cent. Ahead of the crack tip, the principal tensile stresses exceed the principal shear stresses by a factor of 10, approximately.     

A Peierls criterion for the onset of deformation twinning at a crack tip

A criterion for the onset of deformation twinning (DT) is derived within the Peierls framework for dislocation emission from a crack tip due to Rice (J. Mech. Phys. Solids 40(2) (1992) 239). The critical stress intensity factor (SIF) is obtained for nucleation of a two-layer microtwin, which is taken to be a precursor to DT. The nucleation of the microtwin is controlled by the unstable twinning energyγ_ut, a new material parameter identified in the analysis. γ_ut plays the same role for DT as γ_us, the unstable stacking energy introduced by Rice, plays for dislocation emission. The competition between dislocation emission and DT at the crack tip is quantified by the twinning tendencyT defined as the ratio of the critical SIFs for dislocation nucleation and microtwin formation. DT is predicted when T>1 and dislocation emission when T<1. For the case where the external loading is proportional to a single load parameter, T is proportional to . The predictions of the criterion are compared with atomistic simulations for aluminum of Hai and Tadmor (Acta Mater. 51 (2003) 117) for a number of different crack configurations and loading modes. The criterion is found to be qualitatively exact for all cases, predicting the correct deformation mode and activated slip system. Quantitatively, the accuracy of the predicted nucleation loads varies from 5% to 56%. The sources of error are known and may be reduced by appropriate extensions to the model.     

Stress distribution at the tip of a growing crack

High-speed motion-picture photography and an optical method of stress analysis have been used to study the distribution of elastic stress fields at the tip of a crack growing at a fast rate. The existence of several specific properties of the field characteristic of fast crack propagation rates has been established, and the results obtained are used to explain the branching of cracks.The authors convey their thanks to G. I. Barenblatt for sponsoring this work.The authors convey their thanks to G. I. Barenblatt for sponsoring this work.     

Heat generation at the tip of a moving crack

The high energy concentration at the tip of a moving crack causes irreversible deformations and produces heat as a consequence. The resulting temperatures were calculated by consideration of the crack tip as a moving heat-source of rectangular shape. In brittle materials with very small plastic zones and high crack velocities, these temperatures are predicted to be higher than 1000 K. For the experimental verification of these calculations, a very sensitive radiation thermometer was developed. It registers the intensity of the radiation at four wavelengths. By comparison of these intensities with that of black body radiation, the temperature was determined as 3200 K for glass and 4700 K for quartz.     

Singular stress at the ends of a brittle crack for zero value of the J integral

Novosibirsk State Technical University, 630092 Novosibirsk. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 36, No. 6, pp. 108–112, November–December, 1995     

A note on the stress singularity at a non-uniformly moving crack tip

The stress-singularity at a crack tip moving arbitrarily in an elastic plate under plane strain conditions is investigated. By formulating the wave-equations in a polar coordinate system attached to the crack-tip, it is found by an asymptotic analysis that the angle-dependence of the singularity is only dependent on the instantaneous cracktip velocity. This result is used to derive a relation between the dynamic stress-intensity factor and the energyrelease rate.     

Reverse shear on inclined planes at the tip of a sharp crack

Plane strain plastic yielding at a crack tip has been represented by edge dislocations with Burgers vectors parallel to symmetrical planes inclined at 70° and 45° to the plane of the crack. The plastic displacement and the stresses near the crack tip were calculated by a numerical method and the effect of a reduction in applied stress was determined. Removal of the whole or a part of the initial load produces reverse shear in regions of the slip band nearest the crack tip. The amount of reverse shear depends only on the reduction in the load and not on its initial value. The reverse shear is associated with the presence of negative dislocations and the stresses near the crack tip may become compressive even though the applied (remote) stress is still tensile. The degree and extent of compression depends on the reduction in applied stress and on its original value. It is argued that the residual compressive stresses produced under fluctuating loads may produce crack closure and crack arrest. The effect of residual plasticity in a slip band left behind a growing crack has been estimated. It is shown that after an overload the excess residual plasticity opposing crack opening rises to a maximum value when the crack tip has advanced some distance from the point where the overload was applied.     

Energy domain integral applied to solve center and double-edge crack problems in three-dimensions

A three-dimensional Boundary Element Method (BEM) implementation of the energy domain integral for the numerical computation of the crack energy release rate is presented in this paper. The domain expression of the energy release rate is naturally compatible with the BEM, since stresses, strains and derivatives of displacements at internal points can be evaluated using the appropriate boundary integral equations. The pointwise crack energy release rate is evaluated along the three-dimensional crack front over a cylindrical domains that surround a segment of the crack front. The accuracy of the implementation is demonstrated by solving several problems, which include geometries containing straight as well as curved crack fronts.     

On the temperature rise at the tip of a fast running crack2

The high energy concentration at the tip of a running crack leads to irreversible deformations, and a great amount of the deformation energy is set free as heat. Assuming that this moving heat source is of circular shape, the temperature distribution around the crack tip has been calculated. The temperatures are dependent on the radius of the heat source and the crack velocity. Some examples for the material glass are given. The very high temperatures computed lead to the supposition that the observed light emission during fast fracture is of thermal origin.     

Energy condition at the end of a nonstationary crack

The plane problem of the theory of elasticity is considered. It is assumed that in the neighborhood of the tip of an arbitrarily moving crack the stresses have a singularity of order r^–1/2. On this assumption a general expression is obtained for the distribution of the stresstensor components in the given neighborhood. This distribution is determined by the two parameters N and P. In the case of stresses symmetrical about the line of the crack (P=0) the angular distribution does not depend on the intensity coefficient N and is determined only by the velocity of the crack at the given instant and the transverse and longitudinal wave velocities. On the same assumptions it is shown that the energy condition obtained by Craggs for the particular case of steady-state motion is a necessary condition for the arbitrarily moving crack. Irwin [1] and Cherepanov [2] have studied these questions in the quasi-static approximation.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 10, No. 3, pp. 175–178, May–June, 1969.     

Strength reduction of metals upon hydrogen chemisorption at the tip of a crack

When hydrogen interacts with the freshly formed surface, the interatomic bonds at the tip of a crack in the single crystal of a metal are known to break. The decrease in the brittle strength of the single crystals of metals in the presence of hydrogen in the crack compared with the strength of the same metal in the absence of hydrogen is estimated quantitatively using comparative criteria of brittle strength. Lavrent'ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 173–178, May–June, 1998.     

Deformation and stability loss of a part of the atomic chain at the tip of a crack

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 3, pp. 160–172, May–June, 1993.     

Calculation of the evolution of plastic yielding at the tip of a crack and related phenomena

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 3, pp. 154–160, May–June, 1993.     

Yielding on inclined planes at the tip of a crack loaded in uniform tension

Plane-strain yielding from a crack in an infinite elastic body is represented here by a distribution of edge dislocations on two planes inclined at angles ±ga to the crack plane, and the equilibrium condition is solved numerically. Approximate analytical expressions are obtained for the plastic-zone length, the crack opening displacement, and the J-integral, as functions of the applied stress and α. A comparison with a co-planar model of the plastic zone gives very similar results for α ≈ 65°. It is shown that fracture criteria based either on a critical crack opening displacement (COD) or on a critical value of J are always different, and the use of the former may lead to critical defect-sizes which are twice as large as those given by the latter. Furthermore, COD appears not to be a well-defined material property. The critical J criterion gives a fracture stress which is α-dependent : this may be responsible for deviations towards results of linear elastic fracture mechanics when post-yield fracture mechanics is used to analyse extensive yielding. The changes in the stress field of the crack due to the existence of the plastic zone are also discussed.     

Finite strains at the tip of a crack in a sheet of hyperelastic material: I. Homogeneous case

This paper describes an asymptotic analysis of the strain and stress fields at the tip of a crack in a sheet of incompressible hyperelastic material. The investigations are carried out within the framework of finite elastostatics and for the class of Generalized Neo-Hookean materials. Both the symmetric (mode I) and non-symmetric (mixed-mode) cases are considered. It is shown that the latter situation corresponds locally to a rigid body rotation of the symmetric fields. The effect of the hardening parameter on crack tip blunting is investigated analytically and numerically.     

Asymptotic mechanical fields at the tip of a mode I crack in rubber-like solids

Asymptotic analyses of the mechanical fields in front of stationary and propagating cracks facilitate the understanding of the mechanical and physical state in front of crack tips, and they enable prediction of crack growth and failure. Furthermore, efficient modelling of arbitrary crack growth by use of XFEM (extended finite element method) requires accurate knowledge of the asymptotic crack tip fields. In the present work, we perform an asymptotic analysis of the mechanical fields in the vicinity of a propagating mode I crack in rubber. Plane deformation is assumed, and the material model is based on the Langevin function, which accounts for the finite extensibility of polymer chains. The Langevin function is approximated by a polynomial, and only the term of the highest order contributes to the asymptotic solution. The crack is predicted to adopt a wedge-like shape, i.e. the crack faces will be straight lines. The angle of the wedge and the order of the stress singularity depend on the hardening of the strain energy function. The present analysis shows that in materials with a significant hardening, the inertia term in the equations of motion becomes negligible in the asymptotic analysis. Hence, there is no upper theoretical limit to the crack speed.     

Anisotropic plastic stress field at a mixed-mode crack tip under plane and anti-plane strain

On condition that any perfectly plastic stress component at a crack tip is nothingbut the function ofθ.by making use of equilibrium equations,anisotropic plastic stress-strain-rate relations,compatibility equations and Hill anisotropic plastic yieldcondition,in the present paper,we derive the generally analytical expressions of theanisotropic plastic stress field at a mixed-mode crack tip under plane and anti-planestrain.Applying these generally analytical expressions to the mixed-mode cracks,wecan obtain the analytical expressions of anisotropic plastic stress fields at the tips ofmixed-modeⅠ-Ⅲ,Ⅱ-ⅢandⅠ-Ⅱ-Ⅲcracks.     

On perfectly stress field at a mixed-mode crack tip under plane and anti-plane strain

In [1], under the condition that all the perfectly plastic stress components at a crack tip are functions of ϕ only, making use of equilibrium equations, stress-strain rate relations, compatibility equations and yield condition. Lin derived the general analytical expressions of the perfectly plastic stress field at a mixed-mode crack tip under plane and anti-plane strain. But in [1] there were several restrictions on the proportionality factor γ in the stress-strain rate relations, such as supposing that γ is independent of ϕ and supposing that γ=c or cr⁻¹. In this paper, we abolish these restrictions. The cases in [1], γ=cr^d (n=0 or-1) are the special cases of this paper.