Brittle fracture dynamics with arbitrary paths I. Kinking of a dynamic crack in general antiplane loading

We present a new method for determining the elasto-dynamic stress fields associated with the propagation of anti-plane kinked or branched cracks. Our approach allows the exact calculation of the corresponding dynamic stress intensity factors. The latter are very important quantities in dynamic brittle fracture mechanics, since they determine the crack path and eventual branching instabilities. As a first illustration, we consider a semi-infinite anti-plane straight crack, initially propagating at a given time-dependent velocity, that changes instantaneously both its direction and its speed of propagation. We will give the explicit dependence of the stress intensity factor just after kinking as a function of the stress intensity factor just before kinking, the kinking angle and the instantaneous velocity of the crack tip.     

Brittle fracture dynamics with arbitrary paths. II. Dynamic crack branching under general antiplane loading

The dynamic propagation of a bifurcated crack under antiplane loading is considered. The dependence of the stress intensity factor just after branching is given as a function of the stress intensity factor just before branching, the branching angle and the instantaneous velocity of the crack tip. The jump in the dynamic energy release rate due to the branching process is also computed. Similar to the single crack case, a growth criterion for a branched crack is applied. It is based on the equality between the energy flux into each propagating tip and the surface energy which is added as a result of this propagation. It is shown that the minimum speed of the initial single crack which allows branching is equal to 0.39c, where c is the shear wave speed. At the branching threshold, the corresponding bifurcated cracks start their propagation at a vanishing speed with a branching angle of approximately 40°.     

Asymptotic analysis of growing plane strain tensile cracks in elastic-ideally plastic solids

An exact asymptotic analysis is presented of the stress and deformation fields near the tip of a quasistatically advancing plane strain tensile crack in an elastic-ideally plastic solid. In contrast to previous approximate analyses, no assumptions which reduce the yield condition, a priori, to the form of constant in-plane principal shear stress near the crack tip are made, and the analysis is valid for general Poisson ratio ν. Specific results are given for ν = 0.3 and 0.5, the latter duplicating solutions in previous work by L.I. Slepyan, Y.-C. Gao and the present authors. The crack tip field is shown to divide into five angular sectors of four different types ; in the order in which these sweep across a point in the vicinity of the advancing crack, they are : two plastic sectors which can be described asymptotically (i.e., as r → 0, where r is distance from the crack tip) in slip-line terminology as ‘constant stress’ and ‘centered fan’ sectors, respectively ; a plastic sector of non-constant stress which cannot be described asymptotically in terms of slip lines; an elastic unloading sector; and a trailing plastic sector of the same type as that directly preceding the elastic sector. Further, these four different sector types constitute the full set of asymptotically possible solutions at the crack tip. As is known from prior work, the plastic strain accumulated by a material point passing through such a moving ‘centered fan’ sector is O(ln r) as r → 0 ; it is proved in the present work that the plastic strain accumulated by a material point passing through the ‘constant stress’ sector ahead of a growing crack must be less singular than In r as r → 0. As suggested also in earlier studies, the rate of increase of opening gap δ at a point currently at a distance r behind, but very near, the crack tip is given for crack advance under contained yielding by

$\text{δ}\text{?}\text{=}\text{α}\text{J}\text{?}\text{σ}{}_0\text{+β(}\text{σ}{}_0\text{E}\text{)}\text{a}\text{?}\text{ln(}\text{R}\text{r}\text{)}$

where a is crack length, σ0 is tensile yield strength, E is Young's modulus, J is the value of the J-integral taken in surrounding elastic material, and the parameters α and R are undetermined by the asymptotic analysis. The exact solution for ν = 0.3 gives β = 5.462, which agrees very closely with estimates obtained from finite element solutions. An approximate analysis based on use of slip line representations in all plastic sectors is outlined in the Appendix.     

Characterization of crack tip stresses in plane-strain fracture specimens having weld center crack

Fracture toughness of metals depends strongly on the state of stress near the crack tip. The existing standards (like R-6, SINTAP) are being modified to account for the influence of stress triaxiality in the flaw assessment procedures. These modifications are based on the ability of so-called ‘constraint parameters’ to describe near tip stresses. Crack tip stresses in homogeneous fracture specimens are successfully described in terms of two parameters like J–Q or J–T. For fracture specimens having a weld center crack, strength mismatch ratio between base and weld material and weld width are the additional variables, along with the magnitude of applied loading, type of loading, and geometry of specimen that affect the crack tip stresses. In this work, a novel three-parameter scheme was proposed to estimate the crack tip opening stress accounting for the above-mentioned variables. The first and second parameters represent the crack tip opening stress in a homogeneous fracture specimen under small-scale yielding and are well known. The third parameter accounts for the effect of constraint developed due to weld strength mismatch. It comprises of weld strength mismatch ratio (M, i.e. ratio of yield strength of weld material to that of base material), and a plastic interaction factor (I_p) that scales the size of the plastic zone with the width of the weld material. The plastic interaction factor represents the degree of influence of weld strength mismatch on crack tip constraint for a given mismatch ratio. The proposed scheme was validated with detailed FE analysis using the Modified Boundary Layer formulation.     

静止裂纹顶端塑性区内应力场

本文对不可压缩的理想塑性材料裂纹顶端塑性区内的应力场进行了数学分析,证明了当塑性区包围着裂纹顶端而应力函数可用分离变匱型的级数展开且该级数展开的首项与第一类渐近解相同时,第一类渐近解即是塑性区内应力场的精确解。本文又提出了第二类渐近解,说明应力场的渐近解不是唯一的。     

Stability of an advancing crack to small perturbation of its path

In this work we investigate the stability of a nominally straight two-dimensional quasistatically growing crack to a small perturbation of its path. Formulae for perturbations of stress intensity factors induced by slight deviation of the crack trajectory were developed by Movchan et al. (Int. J. Solids Struct. 35, 3419) Their solution is exploited to derive an equation for the perturbation of the crack path on the assumption that the crack advances in pure “opening” mode (i.e. local K_II=0). Various types of loading conditions are considered, including a cracked body loaded by a pair of point body forces and a crack whose faces are subjected to given tractions acting in the direction normal to the crack boundary. The body is also subjected to a remotely maintained uniaxial stress, aligned with the direction of the unperturbed crack. The loading is assumed to advance as the crack advances, to maintain the critical value of Mode I stress intensity factor. Numerical computations of possible crack paths have been performed, extending results on crack stability obtained by Cotterell and Rice (Int. J. Fract. 16, 155). The results show that in the case of loading by point body forces the stability of the crack path depends on the positions of the points of application of the applied forces and the magnitude of the applied stress acting parallel to the crack. There exists a critical value of this stress such that the crack path is stable for values less than critical and unstable otherwise. It is shown that the crack is always unstable in the case of point force tractions applied normal to the crack faces.     

Role of elasticity in elastic–plastic fracture: Analytical bi-linear crack-tip fields and finite element analysis

A solution for Model-I plane strain crack tip fields in a bi-linear elastic–plastic material is presented. The elastic–plastic Poisson's ratio is introduced to characterize the influence of elastic deformation on the near tip constraint. Attention is focused on the distribution of elastic/plastic strain energy in the sensitive region of the forward sector ahead of a crack tip. The present study shows that the elastic strain energy can be higher than the plastic strain energy in this sensitive sector while large amount of the plastic strain energy develops outside this sector around the crack tip. The effect of elastic deformation in this sensitive region on the structure of crack-tip fields is considerable and the assumption in some important solutions for crack-tip fields reported in literature that the elastic deformation is small and can be ignored is therefore not physically reasonable. Besides, finite element analysis is carried out to validate the analytical solution and good agreement between them is found. It is seen that the present solution with T-stress can properly describe the crack-tip fields under various constraints for different specimens and an analytical relation is established between the critical value of J-integral, J_c, and T-stress for elastic–plastic fracture.     

On the fracture toughness of ferroelastic materials

The toughness enhancement due to domain switching near a steadily growing crack in a ferroelastic material is analyzed. The constitutive response of the material is taken to be characteristic of a polycrystalline sample assembled from randomly oriented tetragonal single crystal grains. The constitutive law accounts for the strain saturation, asymmetry in tension versus compression, Bauschinger effects, reverse switching, and strain reorientation that can occur in these materials due to the non-proportional loading that arises near a propagating crack. Crack growth is assumed to proceed at a critical level of the crack tip energy release rate. Detailed finite element calculations are carried out to determine the stress and strain fields near the growing tip, and the ratio of the far field applied energy release rate to the crack tip energy release rate. The results of the finite element calculations are then compared to analytical models that assume the linear isotropic K-field solution holds for either the near tip stress or strain field. Ultimately, the model is able to account for the experimentally observed toughness enhancement in ferroelastic ceramics.     

Mixed-mode I/II fields around a crack with a cohesive zone ahead of the crack tip

Considering a cohesive zone ahead of the crack tip, the linear elastic crack problem under Mode I, Mode II or mixed-mode conditions is formulated in an elliptic coordinate system, so that the cohesive surfaces are conveniently represented by straight line segments. It is shown that the displacement and stress fields around the crack tip and the cohesive zone, expressed in terms of elliptic coordinates, have a simple mathematical form, which does not contain a stress singularity at the crack tip due to the existence of the cohesive zone.     

Energy dissipation mechanisms in ductile fracture

The objective is to investigate energy dissipation mechanisms that operate at different length scales during fracture in ductile materials. A dimensional analysis is performed to identify the sets of dimensionless parameters which contribute to energy dissipation via dislocation-mediated plastic deformation at a crack tip. However, rather than using phenomenological variables such as yield stress and hardening modulus in the analysis, physical variables such as dislocation density, Burgers vector and Peierls stress are used. It is then shown via elementary arguments that the resulting dimensionless parameters can be interpreted in terms of competitions between various energy dissipation mechanisms at different length scales from the crack tip; the energy dissipations mechanisms are cleavage, crack tip dislocation nucleation and also dislocation nucleation from a Frank-Read source. Therefore, the material behavior is classified into three groups. The first two groups are the well-known intrinsic brittle and intrinsic ductile behavior. The third group is designated to be extrinsic ductile behavior for which Frank-Read dislocation nucleation is the initial energy dissipation mechanism. It is shown that a material is predicted to exhibit extrinsic ductility if the dimensionless parameter bρ_disl^1/2 (b is Burgers vector, ρ_disl is dislocation density) is within a certain range defined by other dimensionless parameters, irrespective of the competition between cleavage and crack tip dislocation nucleation. The predictions compare favorably to the documented behavior of a number of different classes of materials.     

Dynamic stability of a propagating crack

In this work we investigate the stability of a straight two-dimensional dynamically propagating crack to small perturbation of its path. Willis and Movchan (J. Mech. Phys. Solids 43 (1995) 319; J. Mech. Phys. Solids 45 (1997) 591) constructed formulae for the perturbations of the stress intensity factors induced by a small three-dimensional dynamic perturbation of a nominally plane crack. Their solution is exploited here to derive equations for the in-plane and out-of-plane perturbations of the crack path making use of the Griffith fracture criterion and the principle of “local symmetry” (i.e the crack propagates so that local K_II=0). We consider a crack propagating in a body loaded by a pair of point body forces and subjected to a remote uniaxial stress, aligned with the direction of the unperturbed crack. We assume that the loading follows the crack as the crack advances and is such that the unperturbed crack is subjected to Mode I loading. We perform an analysis of the stability of the dynamic crack in a similar way as in earlier work (Obrezanova et al., J. Mech. Phys. Solids 50 (2002) 57) on the quasistatically advancing crack. We present numerical results illustrating the influence of the crack velocity on the crack stability. Numerical computations of the possible crack paths have been performed which show that at velocities of crack propagation exceeding about one-third of the speed of Rayleigh waves the crack may admit one or more oscillatory modes of instability.     

Wedge indentation into elastic-plastic single crystals, 1: Asymptotic fields for nearly-flat wedge

Asymptotic stress and deformation fields under the contact point singularities of a nearly-flat wedge indenter and of a flat punch are derived for elastic ideally-plastic single crystals with three effective in-plane slip systems that admit a plane strain deformation state. Face-centered cubic (FCC), body-centered cubic (BCC), and hexagonal-close packed (HCP) crystals are considered. The asymptotic fields for the flat punch are analogous to those at the tip of a stationary crack, so a potential solution is that the deformation field consists entirely of angular constant stress plastic sectors separated by rays of plastic deformation across which stresses change discontinuously. The asymptotic fields for a nearly-flat wedge indenter are analogous to those of a quasistatically growing crack tip fields in that stress discontinuities can not exist across sector boundaries. Hence, the asymptotic fields under the contact point singularities of a nearly-flat wedge indenter are significantly different than those under a flat punch. A family of solutions is derived that consists entirely of elastically deforming angular sectors separated by rays of plastic deformation across which the stress state is continuous. Such a solution can be found for FCC and BCC crystals, but it is shown that the asymptotic fields for HCP crystals must include at least one angular constant stress plastic sector. The structure of such fields is important because they play a significant role in the establishment of the overall fields under a wedge indenter in a single crystal. Numerical simulations—discussed in detail in a companion paper—of the stress and deformation fields under the contact point singularity of a wedge indenter for a FCC crystal possess the salient features of the analytical solution.     

Elastic-plastic fields near crack tip growing step-by-step in a perfectly plastic solid

In this paper, based on the three-dimensional flow theory of plasticity, the fundametal equations for plane strain problem of elastic-perfectly plastic solids are presented. By using these equations the elastic-plastic fields near the crack tip growing step-by-step in an elastic incompressible-perfectly plastic solid are analysed.The first order asymptotic solutions for the stress field and velocity fields near the crack tip are obtained. The solutions show the evolution process of elastic unloading domain and the development process of central fan domain and reveal the possibility of the presence of the secondary plastic domain. The second order asymptotic solution for stress field is also presented.     

A numerical analysis of perpendicular cracks under general in-plane loading with a hybrid displacement discontinuity method

In this paper, a numerical analysis of perpendicular cracks under general in-plane loading is performed by using a hybrid displacement discontinuity method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results show that the numerical approach is simple, yet very accurate for calculating numerically stress intensity factors for perpendicular cracks under general in-plane loading.     

Finite deformation analysis of crack tip fields in plastically compressible hardening–softening–hardening solids

Crack tip fields are calculated under plane strain small scale yielding conditions. The material is characterized by a finite strain elastic–viscoplastic constitutive relation with various hardening–softening–hardening hardness functions. Both plastically compressible and plastically incompressible solids are considered. Displacements corresponding to the isotropic linear elastic mode I crack field are prescribed on a remote boundary. The initial crack is taken to be a semi-circular notch and symmetry about the crack plane is imposed. Plastic compressibility is found to give an increased crack opening displacement for a given value of the applied loading. The plastic zone size and shape are found to depend on the plastic compressibility, but not much on whether material softening occurs near the crack tip.On the other hand, the near crack tip stress and deformation fields depend sensitively on whether or not material softening occurs. The combination of plastic compressibility and softening(or softening–hardening) has a particularly strong effect on the near crack tip stress and deformation fields.     

可压缩理想塑性体中逐步扩展裂纹顶端的弹塑性场

本文利用理想塑性固体平面应变问题的基本方程,分析了可压缩理想塑性体中逐步扩展裂纹顶端的弹塑性场,得到了关于应力的渐近场,分析了弹性卸载区的演变过程和修正的中心扇形区的发展过程,预示了出现二次塑性区的可能性,弹性可压缩性的影响明显表现在经典的中心扇形区必需加以修正,垂直于板面方向的应力偏量不再为零,而且随着新裂纹面的形成,裂纹前方的均匀应力场和紧连着的修正的中心扇形区的应力偏量将发生变化,这种变化是由于垂直于板面方向的应力偏量发生变化造成的。     

Thin film cracking modulated by underlayer creep

In devices that integrate dissimilar materials in small dimensions, crack extension in one material often accompanies inelastic deformation in another. In this paper we analyze a channel crack advancing in an elastic film, while an underlayer creeps. The film is subject to a tensile stress. As the underlayer creeps, the stress field in the film relaxes in the crack wake, and intensifies around the crack tip. In a blanket film, the crack can attain a steady velocity, set by two rate processes: subcritical decohesion at the crack tip, and creep in the underlayer. In a thin-film microbridge over a viscous stripe, the crack cannot grow when the bridge is short, and can grow at a steady velocity when the bridge is long. We use a two-dimensional shear lag model to approximate the three-dimensional fracture process, and an extended finite element method to simulate the moving crack with an invariant, relatively coarse mesh. On the basis of the theoretical findings, we propose new experiments to measure fracture toughness and creep laws in small structures. As a byproduct, an analytical formula is found for the growth rate per temperature cycle of a channel crack in a brittle film, induced by ratcheting plastic deformation in a metal underlayer.     

A new criterion for mixed mode fracture initiation based on the crack tip plastic core region

In this work, we propose a new criterion for mixed mode I-II crack initiation angles based on the characteristics of the plastic core region surrounding the crack tip. The shape and size of the plastic core region are thoroughly analyzed under different loading conditions and a new formulation for the non-dimensional variable radius of the core region is presented for mixed mode (K_I−K_II) fracture. The proposed criterion states that the crack extends in the direction of the local or global minimum of the plastic core region boundary depending on the resultant stress state at the crack tip. The results show a well-defined correlation between the plastic core region characteristics and crack extension angles predicted by other criteria. The proposed criterion is formulated for various loading conditions and is compared with other available criteria against the limited available experimental data. It is shown that the proposed criterion provides a better agreement with the experimental data.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

考虑应力松驰复合型裂纹顶端塑性区分析

Ｓｈｉｈ［１］应用奇异单元，获得了不考虑应力松驰小范围屈服条件下复合型裂纹尖端塑性区形状。Ｚ．Ｚ．Ｚｕ等［２］采用Ｒｉｃｅ［５］给出的裂纹尖端应力关系式，利用有限元分析获得了不考虑应力松驰下复合型裂纹尖端塑性区，本文基于静力学中内力与外力平衡条件，用线弹性的全场解代替局部解，给出了考虑应力松驰下复合型裂纹尖端塑性区边界方程，获得了考虑应力松驰下的任意方向的塑性区尺寸及塑性区形状