A two-dimensional model for crack closure effect in plates under bending

A two-dimensional model is proposed for taking into account the establishment of contact on the compression side of crack faces in plates under bending. An approximate but simple method is developed for evaluating reduction of stress intensity factor due to such ‘crack closure’. Analysis is first carried out permitting interference of the crack faces. Contact forces are then introduced on the crack faces and their magnitudes determined from the consideration that the interference is just eliminated. The method is based partly on finite element analysis and partly on a continuum analysis using Irwin's solution for point loads on the crack line.

Study of the random vibration of nonlinear systems by the Gaussian closure technique

A technique is developed to study random vibration of nonlinear systems. The method is based on the assumption that the joint probability density function of the response variables and input variables is Gaussian. It is shown that this method is more general than the statistical linearization technique in that it can handle non-Gaussian excitations and amplitude-limited responses. As an example a bilinear hysteretic system under white noise excitation is analyzed. The prediction of various response statistics by this technique is in good agreement with other available results.

Observation of crack closure using a crack mouth opening displacement gauge

A crack mouth opening displacement (CMOD) gauge has been designed to estimate the crack length as well as the crack opening stress level in an automated fatigue crack propagation test programme. The CMOD gauge accurately predicts crack tip opening levels in K-controlled tests with constant K, K increasing or K decreasing. In all three K-controlled tests with a single 100% overload cycle, the CMOD gauge does not measure the crack tip opening stress level over a large range of crack lengths after the application of the overload. The CMOD gauge measures the stress level at which the overload plastic zone site opens. Caution should thus be exercised in using the crack opening stress level, estimated by the CMOD gauge, to explain fatigue crack propagation under arbitrary load sequences from the viewpoint of crack closure phenomena.

Effects of temper level on the dependence of fatigue crack growth threshold and crack closure on the prior austenitic grain size

An attempt has been made to systematically investigate the effects of microstructural parameters, such as the prior austenite grain size (PAGS), in influencing the resistance to fatigue crack growth (FCG) in the near-threshold region under three different temper levels in a quenched and tempered high-strength steel. By austenitizing at various temperatures, the PAGS was varied from about 0.7 to 96 μm. The microstructures with these grain sizes were tempered at 200 °C, 400 °C, and 530 °C and tested for fatigue thresholds and crack closure. It has been found that, in general, three different trends in the dependence of both the total threshold stress intensity range, ΔK th , and the intrinsic threshold stress intensity range, ΔK eff, th , on the PAGS are observable. By considering in detail the factors such as cyclic stress-strain behavior, environmental effects on FCG, and embrittlement during tempering, the present observations could be rationalized. The strong dependence of ΔK th and ΔK eff, th on PAGS in microstructures tempered at 530 °C has been primarily attributed to cyclic softening and thereby the strong interaction of the crack tip deformation field with the grain boundary. On the other hand, a less strong dependence of ΔK th and ΔK eff, th on PAGS is suggested to be caused by the cyclic hardening behavior of lightly tempered microstructures occurring in 200 °C temper. In both microstructures, crack closure influenced near-threshold FCG (NTFCG) to a significant extent, and its magnitude was large at large grain sizes. Microstructures tempered at the intermediate temperatures failed to show a systematic variation of ΔKth and ΔKeff, th with PAGS. The mechanisms of intergranular fracture vary between grain sizes in this temper. A transition from “microstructure-sensitive” to “microstructure-insensitive” crack growth has been found to occur when the zone of cyclic deformation at the crack tip becomes more or less equal to PAGS. Detailed observations on fracture morphology and crack paths corroborate the grain size effects on fatigue thresholds and crack closure.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Modified Crack Closure Integral (MCCI) For 3-D Problems Using 20-Noded Brick Elements

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Modified crack closure integral (MCCI) for 3-d problems using 20 noded brick elememnts

The Modified Crack Closure Integral (MCCI) technique based on Irwin's crack closure integral concept is very effective for estimation of strain energy release rates G in individual as well as mixed-mode configurations in linear elastic fracture mechanics problems. In a finite element approach, MCCI can be evaluated in the post-processing stage in terms of nodal forces and displacements near the crack tip. The MCCI expressions are however, element dependent and require a systematic derivation using stress and displacement distributions in the crack tip elements. Earlier a general procedure was proposed by the present authors for the derivation of MCCI expressions for 3-dimensional (3-d) crack problems modelled with 8-noded brick elements. A concept of sub-area integration was proposed to estimate strain energy release rates at a large number of points along the crack front. In the present paper a similar procedure is adopted for the derivation of MCCI expressions for 3-d cracks modelled with 20-noded brick elements. Numerical results are presented for centre crack tension and edge crack shear specimens in thick slabs, showing a comparison between present results and those available in the literature.

Assessment of conditional moment closure models of turbulent autoignition using DNS data

A systematic assessment of the submodels of conditional moment closure (CMC) formalism for the autoignition problem is carried out using direct numerical simulation (DNS) data. An initially non-premixed, n-heptane/air system, subjected to a three-dimensional, homogeneous, isotropic, and decaying turbulence, is considered. Two kinetic schemes, (1) a one-step and (2) a reduced four-step reaction mechanism, are considered for chemistry An alternative formulation is developed for closure of the mean chemical source term , based on the condition that the instantaneous fluctuation of excess temperature is small. With this model, it is shown that the CMC equations describe the autoignition process all the way up to near the equilibrium limit. The effect of second-order terms (namely, conditional variance of temperature excess sigma(2) and conditional correlations of species q(ij)) in modeling is examined. Comparison with DNS data shows that sigma(2) has little effect on the predicted conditional mean temperature evolution, if the average conditional scalar dissipation rate is properly modeled. Using DNS data, a correction factor is introduced in the modeling of nonlinear terms to include the effect of species fluctuations. Computations including such a correction factor show that the species conditional correlations q(ij) have little effect on model predictions with a one-step reaction, but those q(ij) involving intermediate species are found to be crucial when four-step reduced kinetics is considered. The "most reactive mixture fraction" is found to vary with time when a four-step kinetics is considered. First-order CMC results are found to be qualitatively wrong if the conditional mean scalar dissipation rate is not modeled properly. The autoignition delay time predicted by the CMC model compares excellently with DNS results and shows a trend similar to experimental data over a range of initial temperatures.