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.     

A thermal-medium crack model

In analyzing the fracture behavior of a cracked thermoelastic material, of much importance are the effects of thermal loadings on the crack growth. Under the consideration of a medium in an opening crack, a thermal-medium crack model is proposed in this paper. The heat flux at the crack surfaces is assumed to depend on the jumps of the temperature and the elastic displacement across the crack. The thermally permeable and impermeable crack models are the limiting cases of a thermal-medium one. The proposed crack model is applied to solve the problem of a Griffith crack in a transversely isotropic material under thermal and mechanical loadings. Using two introduced displacement functions and the Fourier transform technique, the thermoelastic field and the elastic T-stress are determined in explicit forms by using elementary functions. Numerical results are presented to show the effects of the thermal conductivity inside a crack and applied mechanical loadings on the heat flux at the crack faces, the jumps of temperature across the crack and mode-II stress intensity factor in graphics respectively. The obtained results reveal that the mode-II stress intensity factor for a thermal-medium crack in a thermoelastic material depends not only on applied thermal loadings but also on applied mechanical ones.     

The extended finite element method in thermoelastic fracture mechanics

The extended finite element method (XFEM) is applied to the simulation of thermally stressed, cracked solids. Both thermal and mechanical fields are enriched in the XFEM way in order to represent discontinuous temperature, heat flux, displacement, and traction across the crack surface, as well as singular heat flux and stress at the crack front. Consequently, the cracked thermomechanical problem may be solved on a mesh that is independent of the crack. Either adiabatic or isothermal condition is considered on the crack surface. In the second case, the temperature field is enriched such that it is continuous across the crack but with a discontinuous derivative and the temperature is enforced to the prescribed value by a penalty method. The stress intensity factors are extracted from the XFEM solution by an interaction integral in domain form with no crack face integration. The method is illustrated on several numerical examples (including a curvilinear crack, a propagating crack, and a three‐dimensional crack) and is compared with existing solutions. Copyright © 2007 John Wiley & Sons, Ltd.     

Statistical investigation of surface fatigue cracks in large-sized turbine rotor shaft steel

Distribution properties of an initiation life N_i and a propagation life N_p of surface cracks, statistical characteristics of a crack growth rate dl/dN, and a relationship between a scatter of the distributions and a gradient a of S-log N curves in rotating bending fatigue tests were investigated for Ni-Cr-Mo-V steel, using for a large-sized turbine rotor shaft. The distributions of N_i and N_p were expressed as Weibull distributions, and the scatter of them for smooth specimens and for lower stress amplitude σ_a tests were larger than those for notched ones and for higher σ_a tests, respectively. The statistical properties of crack propagation rate were almost similar in both smooth and notched specimens. The relationship between the a and a coefficient of variation η for the distributions of N_i, N_p and a final fracture life N_f was expressed as η = c(a)^−b, where c and b are constants.     

Effect of crack face contact and friction on Brazilian disk specimens—A finite difference solution

A full curvilinear transformation is employed to study the effect of contact and friction on Brazilian disk specimens containing a crack and subjected to concentrated loads at angles 0° <  < 90°. Homogeneous and bimaterial disks made of glass and epoxy are considered. The effect of loading angle and friction coefficient on the stress intensity factors, as well as the contact length is studied. The results are compared to available semi-analytical and finite elements results. It is found that when the crack faces are in contact without stick zones, an increase in friction causes a decrease of the normal gap, tangential shift and stress intensity factors. When stick conditions appear in the contact zone, an increase in the coefficient of friction also results in increasing the stick zone within the contact zone.     

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.     

Automatic near-threshold fatigue crack growth rate measurements at liquid helium temperature

The development of a fully automated test apparatus for near-threshold fatigue crack growth rate measurements in a liquid helium environment is described, and some initial results for AISI 300 series stainless steels are presented. The experimental apparatus consists of a servohydraulic test machine and a cryostat, complete with a minicomputer, a programmable arbituary waveform generator, a programmable digital oscilloscope and a fully automatic liquid helium refill system. The technique uses 6.4 mm thick compact specimens subjected to systematically decreasing loads, with 24 h operation at 40 Hz, the crack growth being continuously monitored by specimen compliance measurements. The results presented in this study include da/dN vs ΔK curves and threshold fatigue stress intensity factors, ΔK_th, at 4 K for AISI 304L, 304LN and 316 stainless steels. The near-threshold fatigue behaviours of these materials are similar, and the fatigue crack growth rate trends at intermediate ΔK levels nearly agree with published results.     

A general method for determining mixed-mode stress intensity factors from isochromatic fringe patterns

A general method is presented for determining mixed-mode stress intensity factors K_I and K_II from isochromatic fringes near the crack tip. The method accounts for the effects of the far-field, non-singular stress, σ_ox. A non-linear equation is developed which relates the stress field in terms of K_I, K_II, and σ_ox to the co-ordinates, r and θ, defining the location of a point on an isochromatic fringe of order N.

Four different approaches for the solution of the non-linear equation are given. These include: a selected line approach in which data analysis is limited to the line θ = π and the K---N relation can be linearized and simplified, the classical approach in which two data points at (r_m, θ_m) are selected where r_m/θ = 0; a deterministic method where three arbitrarily located data points are used; and an over-deterministic approach where m (>3) arbitrarily located points are selected from the fringe field.

Except for the selected line approach, the method of solution involves an iteractive numerical procedure based on the Newton-Raphson technique. For the over-deterministic approach, the method of least squares was employed to fit the K-N relation to the field data.

All four methods provide solutions to 0.1% providing that the input parameters r, θ, and N describing the isochromatic field are exact. Convergence of the iterative methods is rapid (3–5 iterations) and computer costs are nominal. When experimental errors in the measurements of r and θ are taken into consideration, the over-deterministic approach which utilizes the method of least squares has a significant advantage. The method is global in nature and the use of multiple-point data available from the full-field fringe patterns permits a significant improvement in accuracy of K_I, K_II, and σ_ox determinations.     

Relaxation of residual stress with fatigue loading

Investigations have been carried out to study the relaxation of the surface residual stress in 0.23% C steel due to the application of fatigue loading. The residual stress was induced in the specimen by pre-straining and was measured by X-ray back reflection method using Cr-K radiation. The surface residual stress induced, depends on the plastic strain and appears to bear a relation of the type σ_R = σ_o(e_p)^0.78. The decay of the residual stress appears to depend on log N, given by the relation σ_R1 = σ_RO−K log N, where N is the number of fatigue cycles. The constant K depends on the initial value of the residual stress.     

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.     

Fatigue crack initiation from notch root (local-strain damage accumulation process on crack initiation)

Local deformation (i.e. local-strain behavior) at the notch root in a crack initiation process in annealed 0.48% carbon steel was investigated by the real-time fine-grid method. The fatigue crack initiation cycle was controlled by local-strain damage accumulation. For a quantitative expression of cumulative fatigue damage, we propose a new parameter, the “average local-strain accumulation value,” , which is defined by the integration of local-strain histories until crack initiation. The relationship between average local-strain accumulation range, , and crack initiation cycles, N_c, showed a line whose slope was nearly −0.5 on a log-log coordinate graph. This line we term the “local-strain damage accumulation curve.” The mean stress effect in cases of R = −1, 0 on this line ( vs N_c) was very small or negligible. From the results of variable-loading tests, the linear cumulative damage law based on the local-strain value was also confirmed.     

Effects of welding residual stresses on fatigue crack growth behaviour in butt welds of a pipeline steel

In the present test the fatigue crack growth rate in the parent plate, weld and cross-bond regions was measured and the results were correlated with the stress intensity range ΔK and the effective stress intensity range ΔK_eff. It is indicated that the welding residual stresses strongly affect the crack growth rate. For the weld metal and cross-bond compact tension specimens in which crack growth is along the weld line the fatigue crack growth rate increases as the crack grows. However, for the T compact tension specimen in which crack growth is perpendicular to the weld line at a constant value of applied ΔK the crack growth rate initially decreases as the crack grows. Particularly, at a low constant value of applied ΔK the crack growth rate obviously decreases and the crack fails to grow after short crack growth. When the crack grows to intersect the welded zone, the fatigue crack growth rate gradually increases as the crack grows further. It is clear that the effect of welding residual stresses on the crack growth rate is related to the position of the crack and its orientation with respect to the weld line. Finally, the models of welding residual stress redistribution in the compact tension specimens with the growing crack and its influence on the fatigue crack closure are discussed. It appears that for a butt-welded joint one of the crack closure mechanisms may be considered by the bend or rotation deformation of crack faces due to the welding residual stress redistribution as the fatigue crack grows in the welded joint.     

A micromechanical analysis of fiber crack propagation in composite materials due to transverse tensile loads

The fiber crack propagation in composites due to transverse tensile loads is studied using a micromechanical model and Linear Elastic Fracture Mechanics. To approach the problem, a three domain cylindrical model is introduced to simulate the fiber cracking. The model problem is then solved by the dislocation and singular integral equation techniques. The stress intensity factors of the fiber crack are calculated for various situations. It is found that fiber anisotropy has hardly any effect on the fiber crack propagation; “reverse composites” (composites in which fiber is less stiff than the matrix, such as Nicalon/SiC ceramic composite) virtually eliminate fiber crack propagation.     

An indirect boundary element method for three-dimensional dynamic fracture mechanics

Indirect boundary element methods (fictitious load and displacement discontinuity) have been developed for the analysis of three-dimensional elastostatic and elastodynamic fracture mechanics problems. A set of boundary integral equations for fictitious loads and displacement discontinuities have been derived. The stress intensity factors were obtained by the stress equivalent method for static loading. For dynamic loading the problem was studied in Laplace transform space where the numerical calculation procedure, for the stress intensity factor K_I(p), is the same: as that for the static problem. The Durbin inversion method for Laplace transforms was used to obtain the stress intensity factors in the time domain K_I(t). Results of this analysis are presented for a square bar, with either a rectangular or a circular crack, under static and dynamic loads.     

Impact fatigue properties of epoxy resin filled with SiO2 particles

Impact fatigue tests were carried out on epoxy resin filled with SiO₂ particles. The effects of the percentage of SiO₂ particles and the impact cyclic loading frequency on the impact fatigue strength was investigated. The micromechanism of impact fatigue failure was examined and correlated with the morphology of the fracture surface. The impact stress amplitude, σ_t, can be estimated by the formula, σ₂(N_f · Te)^mt = D_t where (N_f· Te) is the cumulative duration time, and mt and D_t, are parameters describing impact fatigue characteristics. The impact fatigue strength and the static strength are governed by the percent of SiO₂ particles. Crack initiation under monotonie cyclic impact loading was attributed to decision of the epoxy-SiO₂ interface. Unstable crack propagation occurs when the crack passes through the SiO₂ particles.     

Flexural fatigue characteristics of sandwich structures at different loading frequencies

The effects of frequency on the fatigue behavior of S2 glass fiber–vinylester reinforced sandwich composites with two different PVC cores have been investigated. Flexural fatigue tests were performed on sandwich beams with core densities of 130 and 260 kg/m³ at frequencies of 3 and 15 Hz, at a stress ratio, R=0.1 and at four different load levels viz. 90%, 85%, 80% and 75% of the ultimate load. S–N diagrams were generated and it was observed that the fatigue strength increased with core density, and the number of cycles to failure, N_f, increased with increase in frequency. In all cases failure was dominated by a primary shear crack in the core however, the crack path and crack propagation rates varied with frequency. The fatigue crack growth rate (FCGR) in the core of the H130 sandwich beams was subsequently investigated and the relationship between the crack growth rate, da/dN, and the cyclic stress intensity range, ΔK, was determined. It was found that crack growth rate decreased with increase in loading frequency.     

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 fatigue crack propagation model for strain hardening materials

In this paper a crack propagation model based on Tomkins concept (dl/dN ∝ Δε_p · ω) has been developed using the theoretically developed cyclic plastic zone sizes. The crack propagation rates are found to be functions of stress intensity factor, Elber's effective stress range ratio, cyclic yield strength of material, crack length, specimen width and cyclic strain hardening exponent. Suitably grouped to give the crack growth rate in terms of five constants termed as Loading Constant, Material constant, Crack size constant, specimen Width Constant and Stress Intensity Exponent. The crack growth rates found by theory are compared with the experimental results available in literature and a good agreement is found.     

Transient response of a crack in an anisotropic strip

Summary The problem of an anisotropic elastic strip containing a crack which is opened by stresses suddenly applied on the crack faces is considered here. The problem is reduced to a set of simultaneous Fredholm integral equations of the second kind which may be solved iteratively. Once the solutions of these integral equations are obtained, the dynamic stress intensity factors may be evaluated numerically. Numerical results are obtained for a particular transversely isotropic strip.With 1 Figure     

Thermoelastic analysis for an opening crack in an orthotropic material

The thermoelastic analysis of an opening crack embedded in an orthotropic material is made under applied uniform heat flux and mechanical loadings. To simulate the case of an opening crack filled with a medium, a thermal-medium crack model is proposed. The thermally permeable and impermeable cracks are the limiting ones of the proposed thermal-medium one. The crack-tip thermoelastic fields induced by a crack in an orthotropic material are determined in closed forms. The elastic T-stress can be also obtained explicitly. The effects of applied mechanical loadings and the thermal conductivity of crack interior on the heat flux at the crack surfaces and the mode-II stress intensity factor are investigated through numerical computations. The obtained results reveal that an increase of the thermal conductivity of crack interior decreases the mode-II stress intensity factor. And when an applied mechanical loading is increasing, the mode-II stress intensity factor is rising.