Path prediction of I + III mixed mode fatigue crack propagation

In this paper, compact tension specimens with tilted cracks under monotonic fatigue loading were tested to investigate I + III mixed mode fatigue crack propagation in the material of No. 45 steel with the emphasis on the propagation rate expression and the path prediction. It is found that during the mode transformation process, the crack propagation rate is still controlled by the mode I stress intensity factor; and Paris equation also holds for the relationship between and ΔK_I . Crack propagation path can be predicted only when both the crack mode transformation rate and propagation rate are available.     

Experimental study of I + III mixed mode fatigue crack transformation propagation

In this paper, compact tension specimens with tilted cracks under monotonic fatigue loading were tested to investigate I + III mixed mode fatigue crack propagation in the material of No. 45 steel with the emphasis on the mode transformation process. It is found that with the crack growth, I + III mixed mode changes to Mode I. Crack mode transformation is governed by the Mode III component and the transformation rate is a function of the relative magnitude of the Mode III stress intensity factor. However, even in the process of the crack mode transformation the fatigue crack propagation is controlled by the Mode I deformation.     

Evaluation of fatigue crack propagation by mode I and mixed mode in 1050 aluminium

The mechanism of mixed‐mode fatigue crack propagation was investigated in pure aluminum. Push‐pull fatigue tests were performed using two types of specimens. One was a round bar specimen having a blind hole, one was a plate specimen having a slit. The slit direction cut in the specimen was perpendicular or inclined 45 degrees relative to the centre of the specimen axis. In both cases, cracks propagated by mode I or by the mixed mode combining mode I and shear mode, depending on the testing conditions. In these cases the crack propagation rate was evaluated with a modified effective stress intensity factor range. Crack propagation retardation was observed in some specimens. However, it was found that the crack propagation rate could also be evaluated by the effective stress intensity factor range independent of the crack propagation mode.     

Can pure mode III fatigue loading contribute to crack propagation in metallic materials?

The possibility of pure mode III crack growth is analysed on the background of theoretical and experimental results obtained in the last 20 years. Unlike for modes I and II, there is no plausible micromechanistic model explaining a pure mode III crack growth in ductile metals. In order to realize 'plain' mode III fracture surface, we propose the propagation of a series of pure mode II cracks along the crack front. Fractographical observations on crack initiation and propagation in a low alloy steel under cyclic torsion support such a model. The authors have not seen any clear indication of a pure mode III crack growth micromechanism in metals until now.     

An embedded cohesive crack model for finite element analysis of mixed mode fracture of concrete*

An embedded cohesive crack model is proposed for the analysis of the mixed mode fracture of concrete in the framework of the Finite Element Method. Different models, based on the strong discontinuity approach, have been proposed in the last decade to simulate the fracture of concrete and other quasi‐brittle materials. This paper presents a simple embedded crack model based on the cohesive crack approach. The predominant local mode I crack growth of the cohesive materials is utilized and the cohesive softening curve (stress vs. crack opening) is implemented by means of a central force traction vector. The model only requires the elastic constants and the mode I softening curve. The need for a tracking algorithm is avoided using a consistent procedure for the selection of the separated nodes. Numerical simulations of well‐known experiments are presented to show the ability of the proposed model to simulate the mixed mode fracture of concrete.     

Mixed mode fatigue crack growth: A literature survey

The applications of fracture mechanics have traditionally concentrated on crack growth problems under an opening or mode I mechanism. However, many service failures occur from growth of cracks subjected to mixed mode loadings. This paper reviews the various criteria and parameters proposed in the literature for predictions of mixed mode crack growth directions and rates. The physical basis and limitations for each criterion are briefly reviewed, and the corresponding experimental supports are discussed. Results from experimental studies using different specimen geometries and loading conditions are presented and discussed. The loading conditions discussed consist of crack growth under mode II, mode III, mixed mode I and II, and mixed mode I and III loads. The effects of important variables such as load magnitudes, material strength, initial crack tip condition, mean stress, load non-proportionality, overloads and crack closure on mixed mode crack growth directions and/or rates are also discussed.     

Mixed mode cracks in Reissner plates

Based on the sixth order Reissner plate theory, the generalized displacement functions for a cracked plate are derived by eigenfunction expansion method. The fractal two-level finite element method is employed to obtain the stress (moment and shear) intensity factors for the center cracked plate subjected to out-of-plane bending and twisting loads. The numerical results from the present method are checked with those available in literature. Highly accurate stress intensity factors are predicted for a wide range of thickness to crack length ratio and a full range of PoissonÆs ratio provided that the radius of fractal mesh to thickness ratio is not less than .     

Influence of shear parameters on mixed–mode fracture of concrete

The influence of the mode II fracture parameters on the mixed mode fracture experimental tests of quasibrittle materials is studied. The study is based on experimental results and numerical analyses. For the numerical study, a procedure for mixed mode fracture of quasibrittle materials is presented. The numerical procedure is based on the cohesive crack approach, and extends it to mixed mode fracture. Four experimental sets of mixed mode fracture were modelled, one from Arrea and Ingraffea and another from a nonproportional loading by the authors, both with bending concrete beams. Two other sets of experimental fracture were modelled, based on double-edge notched testing; in these tests an important mode II is beforehand expected. The numerical results agree quite well with experimental records. The influence of the main parameters for mode II fracture on the mixed mode fracture is studied for the four experimental set of tests and compared with these results. In all them, large changes in the mode II fracture energy hardly modify the numerical results. The tangential and normal stresses along the crack path during the loading proccess are obtained, also with different values of the mode II fracture energy. For the studied experimental tests it is concluded that the crack is initiated under mixed mode but propagated under predominant mode I. This allows a development of mixed mode fracture models, mainly based on standard properties of the material measured by standard methods, avoiding the problems associated with the measurement of mode II fracture parameters, such as mode II fracture energy and cohesion.     

Mixed mode crack in a periodically-layered composite

The problem on a crack in a bimaterial periodically-layered composite is considered. The single finite length crack parallel to the interfaces is loaded by normal opening tractions but the fracture mode is the mixed one as a result of non–symmetric crack location within the layer. The crack is presented as distributed dislocations with unknown density and the problem is reduced to a system of singular integral equations of the first kind. The coefficients of the system are derived from the application of the Green function for a single dislocation which is obtained in a closed form with the help of the representative cell approach. The dependence of the stress intensity factors K_I and K_II upon the geometric and elastic mismatch parameters is examined. The numerical study allowed to point out the cases in which the simplified sandwich model can be employed for the analysis. On the other hand, for the case of very thin and stiff non–cracked layers essentially dissimilar behavior of the stress intensity factors was revealed. In particular, we discovered that K_II may vanish not only for the symmetric crack position in the midplane of the layer but also in several additional ones. For some limiting cases the solution is seen to coincide with known results.     

Numerical analysis of 2-D crack propagation problems using the numerical manifold method

The numerical manifold method is a cover-based method using mathematical covers that are independent of the physical domain. As the unknowns are defined on individual physical covers, the numerical manifold method is very suitable for modeling discontinuities. This paper focuses on modeling complex crack propagation problems containing multiple or branched cracks. The displacement discontinuity across crack surface is modeled by independent cover functions over different physical covers, while additional functions, extracted from the asymptotic near tip field, are incorporated into cover functions of singular physical covers to reflect the stress singularity around the crack tips. In evaluating the element matrices, Gaussian quadrature is used over the sub-triangles of the element, replacing the simplex integration over the whole element. First, the method is validated by evaluating the fracture parameters in two examples involving stationary cracks. The results show good agreement with the reference solutions available. Next, three crack propagation problems involving multiple and branched cracks are simulated. It is found that when the crack growth increment is taken to be 0.5h≤da≤0.75h, the crack growth paths converge consistently and are satisfactory.     

A modified specimen for evaluating the mixed mode fracture toughness of adhesives

This article introduces a specimen geometry that allows the separation of fracture energy release rates G_I and G_II in adhesively joined beams made of disparate materials. The analysis is based upon a Green's functions formulation for shear deformable beams and circumvents the need to employ finite element computations. The current method results in a system of non-singular integral equations, that can be discretized and reduced to a system of algebraic equations which may be solved by common numerical techniques. The analysis accounts for the dimensions and properties of the adhesive and provides results for a wide range of G_I, G_II and their ratio. Those results agree with finite element computational values to within less than 4%.     

Face/core interface fracture characterization of mixed mode bending sandwich specimens

Debonding of the core from the face sheets is a critical failure mode in sandwich structures. This paper presents an experimental study on face/core debond fracture of foam core sandwich specimens under a wide range of mixed mode loading conditions. Sandwich beams with E‐glass fibre face sheets and PVC H45, H100 and H250 foam core materials were evaluated. A methodology to perform precracking on fracture specimens in order to achieve a sharp and representative crack front is outlined. The mixed mode loading was controlled in the mixed mode bending (MMB) test rig by changing the loading application point (lever arm distance). Finite element analysis was performed to determine the mode‐mixity at the crack tip. The results showed that the face/core interface fracture toughness increased with increased mode II loading. Post failure analysis of the fractured specimens revealed that the crack path depends on the mode‐mixity at the crack tip, face sheet properties and core density.     

Influence of effective stress intensity factor range on mixed mode fatigue crack propagation

ABSTRACT The behaviour of fatigue crack propagation of rectangular spheroidal graphite cast iron plates, each consisting of an inclined semi‐elliptical crack, subjected to axial loading was investigated both experimentally and theoretically. The inclined angle of the crack with respect to the axis of loading varied between 0° and 90°. In the present investigation, the growth of the fatigue crack was monitored using the AC potential drop technique, and a series of modification factors, which allow accurate sizing of such defects, is recommended. The rate of fatigue crack propagation db/dN is postulated to be a function of the effective strain energy density factor range, ΔS_eff. Subsequently, this concept is applied to predict crack growth due to fatigue loads. The mixed mode crack growth criterion is discussed by comparing the experimental results with those obtained using the maximum stress and minimum strain energy density criteria. The threshold condition for nongrowth of the initial crack is established based on the experimental data.     

Automated dynamic fracture procedure for modelling mixed-mode crack propagation using explicit time integration brick finite elements

Several fracture codes have been developed in recent years to perform analyses of dynamic crack propagation in arbitrary directions. However, general-purpose, commercial finite-element software which have capabilities to do fracture analyses are still limited in their use to stationary cracks and crack propagation along trajectories known a priori . In this paper, we present an automated fracture procedure implemented in the large-scale, nonlinear, explicit, finite-element code DYNA3D which can be used to simulate dynamic crack propagation in arbitrary directions. The model can be used to perform both generation- and application-phase simulations of self-similar as well as non-self-similar dynamic crack propagation in linear elastic structures without user intervention. It is developed based on dynamic fracture mechanics concepts and implemented for three-dimensional solid elements. Energy approach is used in the model to check for crack initiation/propagation. Dynamic energy release rate and stress intensity factors are determined from far-field finite-element field solutions using finite-domain integrals. Fracture toughness is input as a function of crack-tip velocity, and when the criterion for crack growth is satisfied, an element deletion-and-replacement re-meshing procedure is used along with a gradual nodal release technique to update the crack geometry and model the crack propagation. Direction of crack propagation is determined using the maximum circumferential stress criterion. Numerical simulations of experiments involving non-self-similar crack propagation are performed, and results are presented as verification examples.     

Two parameter approach for elastic-plastic fracture of short cracked specimens under mixed mode loading

For mode-I loading, in order to describe the near-tip stress field in a specimen under large scaled yielding, two parameter approaches such as J-T, J-Q and J-A₂ theories have been developed and proved well for their validity and limit. In this work elastic-plastic finite element analysis were performed to investigate the effects of mode mixity and T-stress upon near-tip stress distribution for a small-scale-yield model with the modified boundary layer and CTS (Compact Tension-Shear) configuration under large-scale-yield state. As the results, some peculiar characteristics were found as follows; As the mode mixity increases, normal stresses _rr and near the crack tip in the small-scale-yield model get significantly affected by the positive T-stress as well as the negative T-stress, while the shear stress _r is little affected by T-stress. Also, the near-tip stress distribution of short cracked CTS specimens under the large-scale-yield state agree fairly well with that of the small-scale-yield model with an appropriate positive T-stress. The two parameters approach with J-integral and T-stress seems to be a good tool for describing the near-tip stress field under a mixed mode loading and large-scale-yield state.     

3D fracture mechanics investigation on surface fatigue crack propagation

Surface fatigue crack propagation is the typical failure mode of engineering structures. In this study, the experiment on surface fatigue crack propagation in 15MnVN steel plate is carried out, and the crack shape and propagation life are obtained. With the concept of ‘equivalent thickness’ brought into the latest three‐dimensional (3D) fracture mechanics theory, one closure model applicable to 3D fatigue crack is put forward. By using the above 3D crack‐closure model, the shape and propagation life of surface fatigue crack in 15MnVN plates are predicted. The simulative results show that the 3D fracture mechanics‐based closure model for 3D fatigue crack is effective.     

Mixed-mode I/II fracture behaviour of a high strength steel

Mixed mode fracture in a certain high strength steel has been investigated through physical experiments and numerical calculations. The main objective has been to investigate the implications of local crack tip processes on the macroscopic mixed mode fracture behaviour. A scrutiny of the fractured specimens revealed evidence of crack branching in all cases where β_{_eq} < 40° ( β_{_eq =} atan ( K_{_I} /K_{_II} ) ). The appearance of branching was found to be accompanied to a rather abrupt increase in the macroscopic mixed mode fracture toughness. In the numerical calculations the effective plastic strain criterion suggested by Hallbäck and Nilsson (1994) was applied to the present material. Crack tip branching was from the analysis predicted to occur at β_{_eq =} 60°. Besides the presence of branching, it was however questionable whether the analysis corroborated the experimentally observed behaviour.     

Modelling the fracture of concrete under mixed loading

A simple and efficient numerical procedure for mixed mode fracture of quasibrittle materials is shown: This technique predicts crack trajectories as well as load-displacement or load-CMOD responses. The model is based on the cohesive crack concept and uses the local mode I approach. Numerical results agree quite well with three experimental sets of mixed mode fracture of concrete beams; one from Arrea and Ingraffea, another from García, Gettu and Carol and from a nonproportional loading by the authors. In constrast to more sophisticated models, this method offers two major advantages: it requires only material properties measured by standardized methods and it can easily be implemented with general multipurpose finite element codes.     

Three-dimensional mixed mode analysis of a cracked body by fractal finite element method

A semi-analytical method namely fractal finite element method is presented for the determination of stress intensity factor for the straight three-dimenisonal plane crack. Using the concept of fractal geometry, infinite many of finite elements are generated virtually around the crack border. Based on the analytical global displacement function, numerous DOFs are transformed to a small set of generalised coordinates in an expeditious way. No post-processing and special finite elements are required to develop for extracting the stress intensity factors. Examples are given to illustrate the accuracy and efficiency of the present method. Very good accuracy (with less than 3% errors) is obtained for the maximum value of SIFs for different modes.     

An analytical method for mixed-mode propagation of pressurized fractures in remotely compressed rocks

An analytical method for mixed-mode (mode I and mode II) propagation of pressurized fractures in remotely compressed rocks is presented in this paper. Stress intensity factors for such fractured rocks subjected to two-dimensional stress system are formulated approximately. A sequential crack tip propagation algorithm is developed in conjunction with the maximum tensile stress criterion for crack extension. For updating stress intensity factors during crack tip propagation, a dynamic fictitious fracture plane is used. Based on the displacement correlation technique, which is usually used in boundary element/finite element analyses, for computing stress intensity factors in terms of nodal displacements, further simplification in the estimation of crack opening and sliding displacements is suggested. The proposed method is verified comparing results (stress intensity factors, propagation paths and crack opening and sliding displacements) with that obtained from a boundary element based program and available in literatures. Results are found in good agreements for all the verification examples, while the proposed method requires a trivial computing time.