Photoelastic study of the stresses near a cleavage crack

Photoelastic analysis has been used to determine the stress distribution near the tip of a cleavage crack. The crack bisects a slender rectangular bar along its length, and forces are applied at one end of it tending to open it out. The region of high stresses is localized in a region of the order of the half-width of the bar in size. Very near the crack tip, the stresses decrease in proportion to the inverse square root of the distance from the crack tip. The maximum gradient of the principal tensile stress lies at an angle of 70 to 80 deg from the plane of the crack. Contour maps of the principal stresses, the 45-deg shear stress and the maximum tensile stress are presented.     

Mode I crack tip field with strain gradient effects

A plane strain mode I crack tip field with strain gradient effects is investigated. A new strain gradient theory is used. An elastic-power law hardening strain gradient material is considered and two hardening laws, i. e. a separation law and an integration law are used respectively. As for the material with the separation law hardening, the angular distributions of stresses are consistent with the HRR field, which differs from the stress results^[19]; the angular distributions of couple stresses are the same as the couple stress results^[19]. For the material with the integration law hardening, the stress field and the couple stress field can not exist simultaneously, which is the same as the conclusion^[19], but for the stress dominated field, the angular distributions of stresses are consistent with the HRR field; for the couple stress dominated field, the angular distributions of couple stresses are consistent with those in Ref. [19]. However, the increase in stresses is not observed in strain gradient plasticity because the present theory is based on the rotation gradient of the deformation only, while the crack tip field of mode I is dominated by the tension gradient, which will be shown in another paper. Supported by the National Science Foundation of China (No. 19704100), Science Foundation of Chinese Academy of Sciences (Project KJ951-1-20), CAS K. C. Wong Post-doctoral Research Award Fund and the Post Doctoral Science Fund of China.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

The effect of couple stresses on the tip of a crack

The effect of couple stresses at a crack tip is investigated by considering two particular problems. A formally exact solution is obtained (for couple-stress and micropolar elasticity) for the case of a semi-infinite crack with a prescribed internal stress. Secondly, the problem of a finite crack in an infinite medium (with couple stresses) under uniform tension at infinity, is solved by matched expansions when the couple stress parameter is small compared with the crack length. In each case it is shown that the energy release rate from a crack tip tends to the classical elastic value as the couple stress (or micropolar) parameter tends to zero.     

Loss of adhesion at the tip of an interface crack

A model is constructed to analyze adhesive bond failure at the tip of an interface crack. The model is based on the assumption that there are zones of bounded cohesive tensile and shear stresses near a crack tip. Within the context of certain broad a-priori assumptions on the distributions of certain stress and displacement components in the cohesive zones, the requirement thatall stresses in the two materials remain bounded provides a method to compute the specific details for these zones. It is assumed that bond failure occurs when the extension of the bond fiber at the crack tip exceeds a critical value. For an interface crack in a uniform tension field computations for two alternate formulations suggest that this failure criterion is independent of the precise distribution of the cohesive stresses, but rather depends only upon their averaged values. Combined loading with a dominant tensile component has also been analyzed. If the critical extension of bond fibers and the maximum value of the cohesive tensile stress are known, the model provides the maximum allowable interface stresses for given crack dimension and material parameters.     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

Numerical study of contact forces for crack closure analysis

Plasticity induced crack closure (PICC) has been widely studied using numerical models. Different numerical parameters can be considered to quantify the opening level, namely one based on the analysis of contact stresses at minimum load. A modified version of this parameter is proposed here, based on nodal contact forces instead of contact stresses. The predictions were found to be similar to those obtained from the contact status of 2nd node behind crack tip. The PICC_contact parameter was also found to be very consistent and adequate for parametric studies of the influence of different physical parameters. The contributions to the opening stress intensity factor of different points along crack flank were found to strongly decrease with distance to crack tip. The cumulative K_open between the crack tip and a distance of 0.1 mm was found to vary from 30% to 100%, increasing with stress ratio, R. Finally, a K solution was developed for punctual forces applied on crack flank and compared with a literature solution for infinite plates. A good agreement was found for plane strain state but significant differences of about 10% were found for plane stress state.     

Photoelastic analysis of the stress field surrounding a fatigue crack

A photelastic analysis was carried out on plane polyester specimens containing a fatigue crack, in order to study the effect of plastic yielding around the crack tip on the elastic stress distribution in the vicinity of the crack. In general, results were in good agreement with values calculated for the case of a sharp-tipped crack. However, very near the crack tip, principal stresses obtained experimentally were slightly lower than calculated stresses, probably due to the bluntness of the fatigue crack. Also lines of constant stress tended to move behind the crack tip, in contrast with the calculated stresses, which occurred further forward over the field of investigation.     

Improved model of the griffith crack

The plane-stress state of a cracked continuous medium in tension is determined using relaxation elements. The stress state is analyzed at the tip of a crack surrounded by a plastically deformed material as a band of localized plastic deformation (LPD) shaped like an elongated ellipse. The plastic deformation considerably decreases the stress concentration at the crack tip. As the localization of the plastic deformation increases, the stresses at the crack sides decrease to zero. The decrease in stresses at the tip is accompanied by an increase in the concentration and gradients of the stresses at the end of theLPD band. Here the region of perturbation of the stress field is comparable with the width of the band. Institute of Physics of Strength and Materials Science, Tomsk 634055. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 6, pp. 132–141, November–December, 1998.     

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.     

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.     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

残余应力下厚壁筒表面裂纹的应力强度因子计算

本文首先介绍了边界元法计算裂纹尖端应力强度因子的基本理论,接着利用边界元法计算了在残余应力下不同厚壁筒内表面椭圆裂纹的应力强度因子,研究了其大不随椭圆裂纹不同而变化的规律,为厚壁筒结构的设计,制造以及疲劳寿命分析提供了许多有价值的参考资料。     

Geometry dependence of crack growth resistance curve for a ligament toughened brittle material

Many brittle materials exhibit a time-independent behaviour, whereby the crack tip stress intensity (K) increases during crack growth (R). This increase is associated with the development of a ligament zone behind the crack tip, the restraining stresses due to the crack-bridging elements within this zone being responsible for the stress intensity increase. Theoretical analyses, based on the assumption of a constant restraining stress within the ligament zone, show that, for a prescribed crack-solid geometry configuration, although the K - R curve shape is independent of the loading conditions up to the full development of a ligament zone, the K and R values associated with this full development are critically dependent on the loading conditions. Particular attention is focused on the difference between bend and tensile loading of a finite width solid containing an edge crack.     

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for theT stresses ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analysed in detail. It is found that normal stresses {ie27-1} and {ie27-2} ahead of the crack tip, are characterised byQ fields if the crack is within a stiff material and the parameters |p _T | and |q _T | are very small, whereQ is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses {ie27-3} and {ie27-4} are dominated by theQ field plusT stress. This work was supported by the Swedish Research Council for Engineering Sciences.     

On the Finite Motions Generated by a Mode I Propagating Crack

The motion field surrounding a rapidly propagating crack, loaded symmetrically about the plane of the crack, is investigated. The problem is formulated within the framework of finite elastodynamics for thin slabs composed of compressible hyperelastic material. Writing the motion equations, the initial and the internal boundary conditions, with respect to a coordinate system that translates with the moving crack tip, we perform an asymptotic local analysis for a traction-free straight crack that suddenly grows at constant velocity. Moreover, the asymptotic Piola–Kirchhoff and Cauchy stress fields are computed, and we discuss the order of singularity of the dynamic stresses. This revised version was published online in August 2006 with corrections to the Cover Date.     

Fatigue crack opening stress based on the strip-yield model

The modified strip-yield model based on the Dugdale model and two-dimensional approximate weight function method were utilized to evaluate the effect of in-plane constraint, transverse stress, on the fatigue crack closure. The plastic zone sizes and the crack opening stresses considering transverse stress were calculated for four specimens: single edge-notched tension (SENT) specimen, single edge-notched bend (SENB) specimen, center-cracked tension (CCT) specimen, double edge-notched tension (DENT) specimen under uniaxial loading. And the crack opening behavior of the center-cracked specimen under biaxial loading was also evaluated. Normalized crack opening stresses σ_op/σ_max for four specimens were successfully described by the normalized plastic zone parameter Δω^′_rev/ω^′ considering transverse stress, where Δω^′_rev and ω^′ are the size of the reversed plastic zone at the moment of first crack tip closure and the size of the forward plastic zone for maximum stress, respectively. The normalized plastic zone parameter with transverse stress also was satisfactorily correlated with the behavior of crack closure for CCT specimen under biaxial loading.     

Crack separation energy-rates for the DBCS model under biaxial modes of loading

A Modified version of the Dugdale-Bilby-Cottrell-Swinden (DBCS) model simulating the effect of plasticity at the tip of a crack in an infinite region was used by kfouri and rice (1978) to calculate the crack separation energy-rate G^Δ corresponding to a finite crack growth step Δa during plane strain mode I crack extension. The loading consisted of a remote uniaxial tension σp applied normally to the plane of the crack. Using Rice's path-independent integral J to characterize the applied load in the crack tip region, and assuming the length R of the crack tip plastic zone to be small compared with the length of the crack, an analytical expression was derived relating the ratios (J/G^Δ) and (2a/R) for small values of (2a/R). The analytical solution was incomplete in itself in that the value assumed in the plastic region of the DBCS model for the normal stress Y acting on the extending crack surfaces was the product of the yield stress in uniaxial tension σ_Y and an unknown parameter C, the value of which depends on the effect of the local hydrostatic stresses in the case of plane strain conditions. The analytical solution was compared with a numerical solution obtained from a plane strain elastic-plastic finite element analysis on a centre-cracked plate (CCP) of material obeying the von Mises yield criterion. The value used for the yield stress was 310 MN/m² and moderate isotropic linear hardening was applied with a tangent modulus of 4830 MN/m². A uniaxial tension σp was applied on the two appropriate sides of the plate. The comparisons showed that the analytical and finite element solutions were mutually consistent and they enabled the value of C to be established at 1.91. In the present paper similar comparisons are made between the analytical solution and the finite element solutions for the CCP of the same material under different biaxial modes of loading. By assuming further that the form of the analytical solution does not depend on the details of the geometry and of the loading at remote boundaries, a comparison has also been made with the results of a finite element analysis on a compact tension specimen (CTS) made of the same material as the CCP. The different values of C obtained in each case are consistent with investigations by other authors on the effect of load biaxiality on crack tip plasticity.     

Dynamic self-similar debonding of interface at very high velocity: Anti-plane case

This paper analyzes the anti-plane problem of dynamic self-similar debonding of interface at very high velocity. The debonding is modeled as an interface crack propagating self-similarly from zero-length. The extending speed is assumed to be transonic or supersonic. We first consider the dynamic debonding under moving concentrated loads. The moving dislocation model of self-similar propagation of an interface crack is used to formulate the problem to a singular integral equation which is solved analytically. The singularity of stresses near the crack tip is discussed and the dynamic stress intensity factors are presented. Finally the solution of dynamic debonding underx ²-type loads is obtained by using the superposition method.     

Stress and strain field near tip of mode III growing crack in materials with creep behavior

Using a proposed constitutive relation for materials with creep behavior, the stress and strain distribution near the tip of a Mode III growing crack is examined. Asymptotic equations of the crack tip field are derived and solved numerically. The stresses remain finite at the crack tip. Obtained qualitatively is the crack tip velocity and the local autonomy of the near tip field solution is discussed.