Some problems in the J-integral

For a blunt crack the j-integral is path dependent on contours which are very close to the crack tip even for elastic material. Using the incremental J-integral theory we introduce a new parameter J_t, characterizing the behavior of a crack tip and prove that the J-integral is almost path independent on contours whose radii are greater than several COD if σ_ij,1Δε_ij — ε_ijΔσ_ij = 0 in plastic regions for elasto-plastic material.     

Some remarks on the Det.-criterion and caustics

A comprehensive study is made of the Det.(σ_ij)-criterion that involves the weighted difference of the dilatational and distortional component of the strain energy density function. Certain inconsistencies and contradictions when applied to the Mode I crack extension problem are pointed out. Clear distinctions are also made on the meaning of the core region in contrast to the loci of yield initiation and the initial curve of the caustics. Examples are provided to illustrate how overspecification in the initial assumption could lead to obscure results.     

Strain energy density fracture criterion in elastodynamic mixed mode crack propagation

Sih's fracture criterion based on strain energy density, S, for mixed mode crack extension under static loading is extended to dynamic mixed mode, K_I and K_II, crack propagation. Influence of the second order term, σ_ox, which represents the non-singular constant stress acting parallel to the direction of crack propagation, on the S distribution surrounding the crack tip, is demonstrated. Numerical studies show that positive σ_ox enhances the fracture angle and negative σ_oxreduces the fracture angle irrespective of the sign of K_II/K_I, when S is measured at a critical distance r_c from the crack tip. This fracture criterion is verified by the crack curving results of dynamic photoelastic fracture specimens. Omission of σ_ox term leads to predicted fracture angles which are at variance with experimental data.     

A higher-order approximation for the T-criterion of fracture in biaxial fields

The T-criterion of fracture is based on the principle that crack propagates when the maximum value of the distribution of the dilatational component of strain energy density T_v, evaluated along contour lines of constant distortional energy density T_D around the crack tip, attains a limiting value T_vo The angle of this maximum defines also the direction of initiation of crack propagation. Then, the study of the distribution of T_v around the crack tip presents a special interest for understanding mechanisms of fracture.

In this investigation an exhaustive theoretical analysis of the distribution of t_v-component around the tip of crack under in-plane modes of loading was undertaken. The T_v-distribution was evaluated along the elastic-plastic boundary, developed around the crack tip for impending plasticity, according to the Mises yield condition (T_D = T_D0 = const.). The mode of loading of the cracked plate was assumed biaxial with different biaxiality ratios k and a two-term approximation for the respective complex stress function was considered, according to the studies of Liebowitz et al.[1], instead of only the singular term considered up-to-now.

It was found that the T_v-distribution along the Mises initial elastic-plastic boundary presents always a maximum in front of the crack tip, whose position and magnitude depend on the biaxiality factor k and the angle of loading β. The position and the magnitude of this maximum for the two-term approximation of φ(z) showed differences in some regions with the respective values for the singular solution.     

A simple procedure for the accurate determination of stress intensity factors by finite element method

A simple procedure for the accurate determination of stress intensity factors K_I, K_II by the conventional finite element method is proposed. The first step of the method is to calculate the stress σ₂ of the plate without a crack. The second step is to calculate the stress σ_tip, of the plate with the crack. The value of (σ_tip−σ_g) at the crack tip element is regarded to have the intimate relation with K_I, K_II K_I, and K_II are determined from the value of (σ_tip−σ_g) and a standard solution. It is shown that the results obtained for many problems by the proposed method are in excellent coincidence with the analytical solutions. The error is below 1–3% for the most cases.     

Weight function calculations for mixedmode fracture problems with the virtual crack extension technique

An efficient finite element method is presented for calculating the stress intensity factors (K_I and K_II) and the weight functions for mixed-mode cracks with one virtual crack extension. The computational efficiency is enhanced through the use of singular elements and the application of colinear virtual crack extension (VCE) technique to symmetric mesh in cracktip neighborhood. This symmetric mesh in crack-tip vicinity permits the analytical separation of strain energy release rate into G_I for Mode I and G_II for Mode II for the mixed fracture problems with the colinear virtual crack extension.

Rice's displacement derivative representation of weight function vector for symmetric crack has been extended to the mixed fracture mode at nodal location (x_i,y_i) with crack length (a) and inclination angle (β) as h_I(II)(x_i, y_i, a, β) = (H/2K_I(II)(∂U_I(II)(x_i, y_i, a, β/∂a).

This equation permits explicit determination of weight functions for the entire structure of a given asymmetric crack geometry with colinear VCE technique. The explicit weight functions for mixed fracture mode depend strongly on the constraint conditions. The method of obtaining the required stress intensity factors of a given asymmetric crack geometry, from the weight function concept under the selected constraint conditions, which are different from constraint conditions used in the available weight functions for the same crack geometry, is also presented in this paper. This is accomplished by combining the predetermined explicit weight functions with the self-equilibrium forces at their application locations. These self-equilibrium forces include both the applied surface tractions and the reaction forces induced from the constraint conditions.     

Experimental micromechanical evaluation of strength of granular materials: Effects of particle rolling

Biaxial compression tests have been performed on assemblies of oval cross-sectional rods, in an effort to evaluate the effects of interparticle friction, particle shape, and initial fabric on the overall strength of granular materials. The variation in the spatial arrangement of the particles (fabric) and particle rolling and sliding are monitored by taking photoelastic pictures at various stages during the course of deformation. Based on this, the following conclusions are obtained. (1) Particle rolling appears to be a major microscopic deformation mechanism, especially when interparticle friction is large. (2) There are relatively few contacts at which relative sliding is dominant, and this seems to be true even when the assembly reaches the overall failure state; this observation is in contradiction to the common assumption that particle sliding is the major microscopic deformation mode. (3) During the course of deformation and up to the peak stress, new contacts are continually formed in such a manner that the contact unit normals tend to concentrate more in a direction parallel to the maximum principal compression. This concentration of unit normals seems to be closely related to the formation of new column-like load paths which carry the increasing axial stress under constant lateral force. After the peak stress, such a column-like microstructure disappears and considerable rearrangement of the load paths takes place, leading to a more diffused (homogeneous) microstructure in the critical state. (4) If a fabric tensor F_ij, i, J = 1,2,3, is defined to be proportional to the volume average of the quantity m_im_j, where m_i are the rectangular Cartesian components of a unit vector along a vector that connects the centroids of two typical contacting granules, then it appears that the overall stress with components σ_ij tends to become coaxial with the fabric tensor F_ij, as the overall deformation continues. For two-dimensional granules the result σ_ij = _OF_ij + β_OF_jkF_kj (k summed) obtained by Mebrabadi, Nemat-Nasser and Oda (1980) by microchemical modeling is confirmed experimentally; _O and β_O are material parameters.     

Fatigue crack initiation and growth under mixed mode loading in aluminum alloys 2017-T3 and 7075-T6

Fatigue crack initiation and growth characteristics under mixed mode loading have been investigated on aluminum alloys 2017-T3 and 7075-T6, using a newly developed apparatus for mixed mode loading tests. In 2017-T3, the fatigue crack initiation and growth characteristics from a precrack under mixed mode loading are divided into three regions—shear mode growth, tensile mode growth and no growth—on the ΔK_I-ΔK_II plane. The shear mode growth is observed in the region expressed approximately by ΔK_II > 3MPa√m and ΔK_II/ΔK_I > 1.6. In 7075-T6, the condition of shear mode crack initiation is expressed by ΔK_II > 8 MPa√m and ΔK_II/ΔK_I > 1.6, and continuous crack growth in shear mode is observed only in the case of ΔK_I/ΔK_II, 0. The threshold condition of fatigue crack growth in tensile mode is described by the maximum tensile stress criterion, which is given by Δσ_θmax √2πr 1.6MPa√m, in both aluminum alloys. The direction of shear mode crack growth approaches the plane in which K_I decreases and K_II increases towards the maximum with crack growth. da/dN-ΔK_II relations of the curved cracks growing in shear mode under mixed mode loading agree well with the da/dN-ΔK_II relation of a straight crack under pure mode II loading.     

Fracture angle and strain-energy-density-factor of a crack at hole at an arbitrary angle

For both the maximum stress criterion and strain-energy-density-factor (S) theory, fracture angle (the initial angle of crack growth) − θ_o is predicted by using opening and sliding mate stress intensity factors, k₁ and k₂. These theoretical predictions are consistent with experimental fracture angles.

For the S theory, the crack spreads in the negative θ_o-direction in a plane for which S is a minimum, S_min. This quantity was obtained analytically. The experimental data of the critical S (S_cr) on plexiglass fracture specimens remains essentially constant.     

The theory of non-steady state fracture kinetics analysis—1: General theory of crack propagation

The kinetics theory of thermally activated time dependent crack propagation is extended to describe the crack size distribution in non-steady state. The distribution is represented by a series of n differential equations, each expressing the rate of crack tip concentration change over the system of n consecutive energy barriers. The general solution for the set of homogeneous linear first order differential equations developed in this report is of the form ρ_i = Σ_jc_i,je^{−λj_t + C_i}, where ρ_i is the crack tip concentration in the ith valley. The theory takes into consideration the discrete character of solids in contrast with the usual continuum models of fracture studies. The analysis is readily applicable to Regions I and II of stress corrosion cracking.     

Characteristics of brittle fracture under general combined modes including those under bi-axial tensile loads

In this paper, the brittle fracture initiation characteristics under general combination of the opening mode (Mode I), sliding mode (Mode II) and tearing mode (Mode III) were investigated both theoretically and experimentally.

First, the perfectly brittle fracture tests were conducted on specimens of PMMA (Polymethylmethacrylate) for all possible combinations of the fracture modes including respective pure modes. The experimental fracture strengths were compared with those predicted by the fracture criteria which are represented in terms of: (1) maximum tangential stress, [σ_gq]_max, extended to general combined modes, (2) maximum energy release rate at the propagation of a small kinked crack, [G^k(γ)]_max, and (3) newly derived maximum energy release rate at the initiation of a small kinked crack, [G(γ)]_max. It was found that the [G^k(γ)]_max or [G(γ)]_max criterion was very effective to predict both the direction of initial crack propagation and the fracture strength. These energy release rates are expressed in closed forms, and the interaction curves of the brittle fracture strength under arbitrary combinations of Modes I, II and III were derived.

Next, for fracture accompanied by plastic deformation, tests were carried out on specimens of mild steel (SM 41) imposing bi-axial tensile loads at various low temperatures. Then, brittle fracture with plastic deformation occurs under a combination of Modes I and II. In the case of brittle fracture with small scale yielding, the [G(γ)]_max criterion predicts well the direction of initial crack propagation but estimates only lower fracture strength than the experimental one. In the cases of brittle fracture with large scale yielding and under general yielding, it was found from the fracture tests that the direction of initial crack propagation was nearly normal to the resultant vector of the crack opening displacements in the opening and sliding modes at the notch tip. To this type of fracture, the modified COD criterion predicts well the direction of initial crack propagation, but lower fracture strength.

When brittle fracture occurs under the influence of plastic deformation, in such cases as the last three mentioned above, the actual fracture strength is higher than what the most reliable criterion predicts and it increases as deformation in Mode II becomes larger.     

The influence of the stress state on the plastic zone size

“Linear elastic” fracture mechanics is based on the “plane strain” condition which means dε_zz = 0. Really stress states can occur in structures or components where dε_zz > 0. Two examples are discussed. Using the v. Mises criterion and the Sneddon equations the influence of an external stress σ^e_zz acting parallel to the crack front on the plastic zone size is calculated and demonstrated by examples.     

Analysis of the temperature dependence of the fracture stress of sharp notched bending specimens of a structural steel with respect to intermediate transitions

The flow and fracture stresses, σ_y^g and σ_f^g (δ = crack tip displacement), of sharply notched bending specimens of a structural steel U St 37-1 are measured in the temperature range from full scale to small scale yielding. The best adaption of the experimental results for σ_f^g is obtained by a curve which exhibits an intermediate transition, i.e. which follows in a temperature range between an upper, T_tM¹ and a lower, T_tl¹, transition temperature to the curve σ_y^g(T) for the flow stress with a constant δ = δ₁. This transition corresponds to that of the slip to the twin nucleated fracture. Two analyses [3,5] according to the local fracture stress, σ_f^*, concept show that the amount and the temperature dependence of σ_f^* are somewhat different for both methods, but that both exhibit an increase of σ_f^* in the transition range. It is concluded that each transition in the nucleation mode of the fracture is connected with such a transition in the fracture stress. It may, however, become indistinct or even be covered by the scatter of the experimental points.     

Influence of stress ratio at baseline loading on delayed retardation phenomena of fatigue crack growth in A553 steel and A5083 aluminium alloy

The delayed retardation phenomena of fatigue crack growth following a single application of tensile overload were investigated under the baseline loading with the stress ratio, R = σ_min/σ_max, ranging from −1 to 0.5 for A553 steel and A5083 aluminium alloy. Two different overload cycles were applied; the one is the case that the ratio of peak stress range to baseline stress range, r = Δσ₂/Δσ₁, is equal to two and the other is the case that the ratio of maximum peak stress to maximum baseline stress, σ_2max/σ_1max, is equal to two. The retardation took place stronger in aluminium than in steel. Under the condition of r = 2 the normalized number of cycles, N_D/N_C, (N_D: the number of cycles during retardation, N_C: the number of cycles required for propagation through the overload-affected-zone size) decreased slightly as the R ratio increased from −1 to 0.5, while under the condition of σ_2max/σ_1max = 2 the N_D/N_C-values increased drastically as the R ratio increased from −1 to 0 (or the overload ratio, r, increased from 1.5 to 2) in both the materials. These retardation behaviors were expressed theoretically according to the model proposed by Matsuoka and Tanaka [1, 3] by using four parameters: the overload ratio, r, the exponent in Paris equation, m, the overload-affected-zone size, ω_D, and the distance at the inflection point, ω_B.     

The mixed mode I/II fracture behaviour of lightly tempered HY130 steel at room temperature

The mixed mode I/II fracture behaviour at room temperature of HY130 steel tempered at 350°C has been investigated using edge-cracked bend bar specimens loaded in anti-symmetric and symmetric four point bend configurations. In all cases fracture occurred by a localized shear decohesion mechanism that could not be characterized by the stress intensity factors, K_I and K_II, but for which the crack tip displacements, δ_I, and δ_II, appear to provide a first level of characterization. The results suggest that fracture is described by a maximum shear criterion, and this is consistent with the present understanding of fibrous fracture micro-mechanisms in the material.     

A fracture criterion for three-dimensional crack problems

A criterion for predicting the growth of three-dimensional cracks is developed on the basis of the strain energy density concept which has been used successfully for treating two-dimensional crack problems. Fracture is assumed to initiate from the nearest neighbor element located by a set of spherical coordinates (r, θ, φ) attached to the crack border. The new fracture surface is described by a locus of these elements whose locations correspond to the strain energy function, dW/dV, being a minimum. The function dW/dV is found to be singular of the type 1/r and is of quadratic form in the three stress intensity factors k₁, k₂ and k₃ expressed through the strain energy density factor S. It is postulated that unstable crack propagation initiates from a region where S reaches a critical value S_cr = r₀(dW/dV)_cr. The locations of failure lying on the fracture surface is determined by holding (dW/dV)_cr = S_min/r₀ constant. The quantity S_min stands for the value of S minimized with respect to θ and φ and r₀ is a radial distance measured from the crack border.

An example of failure prediction for an embedded elliptical crack subjected to both normal and shear loads is presented. According to the S-criterion, fracture initiation takes place at the ends of the minor axis. An unexpected result is that for a narrow elliptical crack and Poisson's ratio of 1/3 the lowest failure load occurs when the uniaxial tensile load makes an angle of approximately 60° with the crack surface and is in the plane of the major axis. This is in contrast to the expectation that the lowest critical load occurs when the uniaxial tension is perpendicular to the crack surface. In the limit as the elliptical crack becomes increasingly narrower, the result reduces to the two dimensional line crack case of Mode I and III loading. The S-criterion is also applied to the failure prediction of three dimensional cracks under compressive loads.     

A model for high cycle fatigue

Failure due to fatigue consists of such macroscopic events as crack initiation and propagation. Microscopic events including microcrack nucleation, microcrack growth and coalescence of some of the microcracks are also important in that such crack interactions can be considered to contribute to the development of a critical defect, i.e. a defect which can self propagate and lead to failure. Since crack initiation is important at high cycles, this paper considers a microcrack and computes its growth to advance a high cycle stress-life (S-N_f) formulation for metals and alloys based on crack initiation. In addition to a material's Burgers' vector and grain size, an indirect effort is also made to include the role of its propensity to cracking through a ratio S(σ_I)/S(σ_f), where S(σ_I) and S(σ_f) simply represent steady state crack spacings at stress amplitudes σ_I and σ_f (endurance limit), respectively. Conceptually, a decrease of this ratio suggests an increased tendency for cracking and vice versa. As shown in the text, the model predicts that the crack initiation period varies increasingly with this ratio, and the value of the steady state crack spacing ratio is in and of itself quite sufficient to model the experimental stress-life data in instances where life is controlled by crack initiation period and not the stage II crack propagation. Because of this limitation, extreme care must thus be taken with regard to its application.     

The micropolar thermoelastic problem of uniform heat-flow interrupted by an insulated penny-shaped crack

The classical problem of uniform heat-flow disturbed by an insulated penny-shaped crack is solved in the context of micropolar elasticity. The mode II stress intensity factor, K_II is found to depend on two new non-dimensional parameters N and τ − N is a measure of the coupling of the displacement field with the microstructure or the medium (0 N √2) and τ is the ratio of a material characteristic length to the crack radius. K_II remains higher than its classical value when N > 0, τ > 0 and attains the classical value as N and τ vanish. A closed-form expression to K_II is obtained in the physically important limiting case of τ → 0 with N fixed. In this limit the relative increment in K_II, over its classical value, is found to be (1 − v')N² where v' is the micropolar Poisson's ratio.     

Fracture mechanisms of ceramic coating during wear

A typical buckling phenomenon of the coating on the wear groove caused by the residual compressed stress was analyzed by the interface fracture mechanics and the buckling theory. It has been found that there is a critical thickness of coating on the wear groove for the buckling. The critical thickness can be calculated by t_b/c_d = [12(1 − v²_f)σ_R/π²E_f]1/2 (here t_b is the coating thickness, c_d the length of the interfacial crack, v_f the Poisson's ratio of the coating, σ_R the residual compressed stress in the coating, and E_f the elastic modulus of the coating).     

Growth mechanism of stress corrosion cracking in high strength steel

Stress corrosion crack growth rates are measured at sveral stress intensity levels for low-tempered 4340 steel in 0.1N H₂SO₄ solution. The characteristics of the growth rates are divided into three regions of stress intensity factors: Region I near K_1SCC; Region III near unstable fracture toughness, K_1SC; and Region II, which lies between the two. K_1SCC is the value of K at which no crack growth can be detected after 240 hr.

In order to explain these experimental results, the crack initiation analysis reported in a previous paper is extended to the growth rates. A detached crack initiates and grows at the tip of an already existing crack. When the detached crack reaches the tip of the main crack, the process repeats as a new existing crack.

A relationship between crack growth rate, v, and stress intensity factor, K, is obtained as a function of b/a and a = b + d, where b is the distance from the tip of the main crack to the detached crack, and d is the ydrogen atom saturated domain.

The experimental data are in good agreement with the theoretical values in Region II when a = 0.02 mm, b/a = 0.8, c₁/c₀ = 2.8 for 200°C tempered specimens and a = 0.015 mm, b/a = 0.7, c₁/c₀ = 3.0, ρ_b = 0.055 mm for 400°C tempered specimens, where ρ_b is a fictitious notch radius. The plateau part in Region II for 400°C tempered specimens is also successfully explained by the present theory. For Region III, the value of b/a will be almost equal to 1 because v → ∞ for b/a → 1. On the other hand, for Region I, b/a will be zero, since the value of v becomes negligibly small and no crack growth is observable.