Mixed-mode fracture of ductile thin-sheet materials under combined in-plane and out-of-plane loading

Cracks in thin structures often are subjected to combined in-plane and out-of-plane loading conditions leading to complex mixed mode conditions in the crack tip region. When applied to ductile materials, large out-of-plane displacements make both experimentation and modeling difficult. In this work, the mixed-mode behavior of thin, ductile materials containing cracks undergoing combined in-plane tension (mode I) and out-of-plane shear (mode III) deformation is investigated experimentally. Mixed-mode fracture experiments are performed and full, three-dimensional (3D) surface deformations of thin-sheet specimens from aluminum alloy and steel are acquired using 3D digital image correlation. General characteristics of the fracture process are described and quantitative results are presented, including (a) the fracture surface, (b) crack path, (c) load-displacement response, (d) 3D full-field surface displacement and strain fields prior to crack growth, (e) radial and angular distributions of the crack-tip strain fields prior to crack growth and (f) singularity analysis of the crack-tip strains prior to crack growth. Results indicate that the introduction of a mode III component to the loading process (a) alters the crack tip fields relative to those measured during nominally mode I loading and (b) significantly increases the initial and stable critical crack-opening-displacement. The data on strain fields in both AL6061-T6 aluminum and GM6208 steel consistently show that for a given strain component, the normalized angular and radial strains at all load levels can be reasonably represented by a single functional form over the range of loading considered, confirming that the strain fields in highly ductile, thin-sheet material undergoing combined in-plane tension and out-of-plane shear loading can be expressed in terms of separable angular and radial functions. For both materials, the displacement and strain fields are (a) similar for both mixed-mode loading angles Φ = 30° and Φ = 60° and (b) different from the fields measured for Mode I loading angle Φ = 0°. Relative to the radial distribution, results indicate that the in-plane strain components do not uniformly exhibit the singularity trends implicit in the HRR theory.     

SEMI-ELLIPTICAL FATIGUE CRACK GROWTH UNDER ROTATING OR REVERSED BENDING COMBINED WITH STEADY TORSION

Abstract— An analysis of the influence of steady torsion loading on fatigue crack growth rates under rotating or reversed bending is presented. Mixed-mode (I + III) tests were carried out on cylindrical specimens in DIN Ck45k steel and results are compared for two different testing machines: rotary bending and reversed bending obtained by cyclic Mode I (Δ K ₁) with or without superimposed static Mode III ( K _III) loading, simulating the real conditions on power rotor shafts where many failures occur. The growth and shape evolution of semi-elliptical surface cracks, starting from a chordal notch on the cylindrical specimen surface, was measured for several Mode III/ Mode I ratios. Results have shown that the steady Mode III loading superimposed on the cyclic mode I leads to a significant reduction in the crack growth rates. It is suggested that this retardation is related to an increase of plastic zone size near the cylindrical surface in association with the interlocking of rough fracture surfaces, friction and fretting debris, leading to a decrease of the ΔK effective at the crack tip profile due to the "crack closure effect". This work provides a contribution to a better understanding of crack growth rates under mixed-mode load conditions thereby allowing one to predict remaining lifetimes and to estimate the risks of pre-cracked rotor shafts.     

Modeling of mixed-mode crack growth in ductile thin sheets under combined in-plane and out-of-plane loading

This paper describes a modeling approach for analyzing mixed-mode crack growth events in ductile thin-sheet materials under large deformation and combined in-plane and out-of-plane loading conditions. The remote mixed-mode I/III loading leads to local mixed-mode I/II/III fields near the crack front. Making use of full-field surface deformation measurements, finite element models of mixed-mode I/III stable tearing events in thin-sheet specimens have been developed. Model predictions have been compared with experimental measurements (a) just prior to initial crack growth and (b) during stable tearing crack growth. Analyses of curvilinear crack growth events are carried out using a nodal release option or a local re-meshing option and using a generalized CTOD parameter with experimentally measured critical CTOD values. Results of this study suggest that the modeling approach can be employed to numerically re-construct experimental crack growth events in thin plate specimens. This offers a viable means of analyzing and understanding the mixed-mode crack growth events and provides a tool for further investigations of 3D crack front fields (which are otherwise unavailable experimentally) and for the study of fracture criteria for stable tearing events.     

Experimental and finite element studies on mode I and mixed mode (I and II) stable crack growth—II. finite element analysis

Finite element studies are presented on both mode I and mixed mode stable crack growth under static loadings through an aluminium (D16AT) alloy. A COD based criterion has been used to predict the load-displacement diagram from initiation to instability. The theoretical predictions are compared with experimental results presented in Part I. Results on computed crack profiles, stress-strain distribution ahead of the crack tip, J integrals, J resistance curves, plastic zones, etc., are included. The study indicates that the load-displacement diagram associated with a mixed mode stable crack growth in a compact tension type of specimen geometry can be predicted reasonably accurately using the criterion of a fixed crack opening displacement at a finite distance behind the crack tip provided the crack is allowed to grow in the direction of initial growth in the finite element analysis. The crack assumes a more blunted profile in a mixed mode than in the mode I at all the stages of stable extension. The distributions of normal stress and strain in the direction perpendicular to the crack extension line, ahead of the current crack tip, have similarities between the mode I and mixed mode, irrespective of loading angle. Both the stress and strain levels increase as the crack extension proceeds. In a mixed mode, the J integral at the onset of crack extension is the lowest compared with the values at the later stages of the extension. Further, the tearing modulus associated with initial kinking is very small; it becomes close to the mode I values at the later stages. The tearing modulus remained approximately constant during the whole mode I stable growth and it had a similar trend subsequent to kinking in a mixed mode. The specific work of crack extension is zero as Δa → 0 and it increases gradually with Δa irrespective of the mode of loading; the actual variation depends on the loading angle. The plastic zone size grows as the stable extension progresses; the growth is approximately the maximum along the crack extension line.     

CRACK FACE INTERACTIONS AND NEAR-THRESHOLD FATIGUE CRACK GROWTH

Rough fracture surfaces usually influence substantially the fatigue growth properties of materials in the regime of low growth rates. Friction, abrasion, interlocking of fracture surface asperites and fretting debris reduce the applied load amplitude to a smaller effective value at the crack tip (“sliding crack closure”, or “crack surface interaction” or “crack surface interference”). The influence of these phenomena on the fatigue crack growth properties of structural steel is discussed and compared for the two kinds of mixed mode loading employed in this work. Mixed mode loading was performed by (A): cyclic mode III + superimposed static mode I and (B): cyclic mode I + superimposed static mode III loading. Such loading cases frequently occur in rotating load-transmission devices. Several differences are typical for these two mixed-mode loading cases. A superimposed static mode I load increases the crack propagation rate under cyclic mode III loading whereas cyclic mode I fatigue crack propagation is retarded when a static mode III load is superimposed. Increase of the R -ratio (of the cyclic mode III load) leads to an insignificant increase of fracture surface interaction and subsequently to a small decrease of the crack growth rate for cyclic mode III loading, whereas higher R -values during cyclic mode I+ superimposed static mode III loading lead to a significant reduction of the crack growth rates.     

Failure criteria for mixed mode delamination in glass fibre epoxy composites

Critical strain energy release rate of glass/epoxy laminates using the virtual crack closure technique for mode I, mode II, mixed-mode I + II and mode III were determined. Mode I, mode II, mode III and mixed-mode I + II fracture toughness were obtained using the double cantilever beam test, the end notch flexure test, the edge crack torsion test and the mixed-mode bending test respectively. Results were analysed through the most widely used criteria to predict delamination propagation under mixed-mode loading: the Power Law and the Benzeggagh and Kenane criteria. Mixed-mode fracture toughness results seem to represent the data with reasonable accuracy.     

Connecting the essential work of fracture, stress intensity factor, Hill’s criterion in mixed mode I/II loading

Ductile sheet structures are frequently subjected to mixed mode loading, resulting that the structure is under the influence of a mixed mode stress field. Instances of interest are when stable crack growth occurs and when the crack-tip is propagating in this complex mixed-mode condition, prior to final fracture. Purposely designed apparatus was built to test thin-sheets of steel (Grade: DX51D) under mixed-mode I/II. These tests, under plane stress conditions, also investigated the effect of thickness on the specific essential work of fracture or the fracture toughness of the material under quasi-static cracking conditions. The fracture toughness is evaluated under incremental mixed-mode loading conditions. The direction of the propagating crack path and fracture type were observed and discussed as the loading mixity was varied. Whilst the specific essential work of fracture or fracture toughness was obtained using the energy approach, the theoretical analysis of the fracture type and direction of crack path were based on the crack tip stresses and fracture criterions of maximum hoop stress and maximum shear stress along with the utilisation of Hill’s theory. For mixed-mode I/II loading, the variation in the fracture toughness contributions ratios are evaluated and used predicatively using the established energy criterion approach to the crack tip stress intensity approach. The comparison between the theoretical directions of the crack path, failure mode propagation are in good agreement with those obtained from experimental testing indicating the definite link between both approaches.     

Influence of overloads and block loading sequences on Mode III fatigue crack propagation in A469 rotor steel

A study has been made of the influence of variable amplitude loading on Mode III (anti-plane shear) fatigue crack propagation in circumferentially-notched cylindrical specimens of ASTM A469 rotor steel (yield strength 621 MN/m²), subjected to cyclic torsional loading. Specifically, transient crack growth behavior has been examined following spike and fully-reversed single overloads and for low-high and high-low block loading sequences, and the results compared to equivalent tests for Mode I (tensile opening) fatigue crack growth. It is found that the transient growth rate response following such loading histories is markedly different for the Mode III and Mode I cracks. Whereas Mode I cracks show a pronounced transient retardation following single overloads (in excess of 50% of the baseline stress intensity), Mode III cracks show a corresponding acceleration. Furthermore, following high-low block loading sequences, the transient velocity of Mode I cracks is found to be less than the steady-state velocity corresponding to the lower (current) load level, whereas for Mode III cracks this transient velocity is higher. Such differences are attributed to the fact that during variable amplitude loading histories. Mode III cracks are not subjected to mechanisms such as crack tip blunting/branching and fatigue crack closure, which markedly influence the behavior of Mode I cracks. The effect of arbitrary loading sequences on anti-plane shear crack extension can thus be analyzed simply in terms of the damage accumulated within the reversed plastic zones for each individual load reversal. Based on a micro-mechanical model for cyclic Mode III crack advance, where the crack is considered to propagate via a mechanism of Mode II shear (along the main crack front) of voids initiated at inclusion close to the crack tip, models relying on Coffin-Manson damage accumulation are developed which permit estimation of the cumulative damage, and hence the crack growth rates, for arbitrary loading histories. Such models are found to closely predict the experimental post-overload behavior of Mode III cracks, provided that the damage is confined to the immediate vicinity of the crack tip, a notion which is consistent with fractographic analysis of Mode III fracture surfaces.     

Study of slant fracture in ductile materials

Slant fracture is widely observed during crack growth in thin sheet specimens made of ductile materials, providing a good case for investigating three-dimensional criteria for mixed-mode ductile fracture. To gain an understanding of slant fracture events and to provide insight for establishing a slant fracture criterion, stable tearing fracture experiments on combined tension-torsion (nominal mixed-mode I/III) specimens and nominal Mode I Arcan specimens made of Al 2024-T3 are analyzed using the finite element method under three-dimensional conditions. Two types of finite element models are considered for the study of slant fracture: (a) combined tension-torsion specimens containing stationary, flat and slant cracks subject to loads corresponding to the onset of crack growth, and (b) stable tearing crack growth with slanting in a nominal Mode I Arcan specimen. Analysis results reveal that there exists a strong correlation between certain features of the crack-front effective plastic strain field and the orientation of the slant fracture surface. In particular, it is observed that (a) at the onset of crack growth in the combined tension-torsion experiments, the angular position of the maximum effective plastic strain around the crack front serves as a good indicator for the slant fracture surface orientation during subsequent crack growth; and (b) during stable tearing crack growth in the Mode I Arcan specimen, which experiences a flat-to-slant fracture surface transition, the crack growth path on each section plane through the thickness of the specimen coincides with the angular position of the maximum effective plastic strain around the crack front.     

Simulation of mixed-mode I/III stable tearing crack growth events using the cohesive zone model approach

A cohesive zone model (CZM) approach is applied to simulate mixed-mode I/III stable tearing crack growth events in specimens made of 6061-T6 aluminum alloy and GM 6208 steel. The materials are treated as elastic–plastic following the \(J_{2}\) flow theory of plasticity, and the triangular cohesive law is employed to describe the traction-separation relation in the cohesive zone ahead of crack front. A hybrid numerical/experimental approach is employed in simulations using 3D finite element method. For each material, CZM parameter values are chosen by matching simulation prediction with experimental measurement (Yan et al. in Int J Fract 144:297–321, 2009), of the crack extension-time curve for the \(30^{\circ }\) mixed-mode I/III stable tearing crack growth test. With the same sets of CZM parameter values, simulations are performed for the \(60^{\circ }\) loading cases. Good agreements are reached between simulation predictions of the crack extension-time curve and experimental results. The variations of CTOD with crack extension are calculated from CZM simulations under both \(30^{\circ }\) and \(60^{\circ }\) mixed-mode I/III conditions for the aluminum alloy and steel respectively. The predictions agree well with experimental measurements (Yan et al. in Int J Fract 144:297–321, 2009). The findings of the current study demonstrate the applicability of the CZM approach in mixed-mode I/III stable tearing simulations and reaffirm the connection between CTOD and CZM based simulation approaches shown previously for mixed-mode I/II crack growth events.     

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.     

THE EFFECT OF MEAN STRESS ON MIXED MODE (I + III) FATIGUE THRESHOLDS

Abstract— Mode I, mode III and mixed mode (I + III) fatigue crack growth threshold tests have been performed on a 3.5% NiCrMoV steel at a range of mean stresses corresponding to load ratios R = 0.06 to 0.5. Angled slit three point bend specimens were used for mixed mode (I + III) tests, circumferentially slit round bars for mode III tests and conventional three point bend specimens for mode I tests.
Test results were divided into those specimens which had failed and those in which no extension of the initial slit had occurred. The results were compared with the response predicted by two existing mixed-mode threshold models. It was found that a model based on the magnitude and direction of mode I crack opening gave results which were in better agreement with the experimental results than a model based on branch crack formation.     

On the influence of rubbing fracture surfaces on fatigue crack propagation in mode III

Fatigue crack growth has been studied under fully reversed torsional loading (R = ?1) using AISI 4340 steel, quenched and tempered at 200°, 400° and 650°C. Only at high stress intensity ranges and short crack lengths are all specimens characterized by a microscopically flat Mode III (anti-plane shear) fracture surface. At lower stress intensities and larger crack lengths, fracture surfaces show a local hill-and-valley morphology with Mode I, 45° branch cracks. Since such surfaces are in sliding contact, friction, abrasion and mutual support of parts of the surface can occur readily during Mode III crack advance. Without significant axial loads superimposed on the torsional loading to minimize this interference, Mode III crack growth rates cannot be uniquely characterized by driving force parameters, such as ΔK_III and ΔCTD_III, computed from applied loads and crack length values. However, for short crack lengths (?0.4 mm), where such crack surface interference is minimal in this steel, it is found that the crack growth rate per cycle in Mode III is only a factor of four smaller than equivalent behaviour in Mode I, for the 650°C temper at $\text{ΔK}{}_{\mathrm{III}}\text{= 45 MPa m}\text{1}\text{2}$.     

Dislocation emission criterion for a wedge crack under mixed mode loading

Dislocation emission criterion for a wedge crack under mixed mode loading was investigated using Airy stress function. The order of singularity at the wedge crack tip due to remote loading was found to vary with the loading mode. The plastic zones for plane stress and plane strain were studied based on von Mises' and Tresca criteria. The dislocation emission criterion was examined for both loading modes. The mechanism of crack propagation was believed to be controlled by dislocation emission. Under an action of Mode I loading, the wedge tip movement occurred when a pair of edge dislocations of Burgers vectors be ⁱ and –be ^–i were emitted from the wedge tip where b and were the magnitude of Burgers vector and the angle between the positive x axis and the line connecting from the tip to dislocation. Similarly, under an action of Mode II loading, the wedge crack tip moved as soon as either an edge dislocation of Burgers vector along the x direction was emitted from its tip or a pair of edge dislocations of Burgers vectors be ⁱ and be ^–i were emitted from the wedge tip. The conventional mechanism of crack propagation based on the energy release rate was not expected to occur. The calculated results for a few special cases were presented and compared with those reported in the literature.     

Experimental and Numerical Investigation of Mixed-Mode Interlaminar Fracture of Carbon-Polyester Laminated Woven Composite by Using Arcan Set-up

Composite materials are widely used in marine, aerospace and automobile industries. These materials are often subjected to defects and damages from both in-service and manufacturing process. Delamination is the most important of these defects. This paper reports investigation of mixed-mode fracture toughness in carbon–polyester composite by using numerical and experimental methods. All tests were performed by Arcan set-up. By changing the loading angle, α, from 0° to 90° at 15° intervals, mode-I, mixed-mode and mode-II fracture data were obtained. Correction factors for various conditions were obtained by using ABAQUS software. Effects of the crack length and the loading angle on fracture were also studied. The interaction j-integral method was used to separate the mixed–mode stress intensity factors at the crack tip under different loading conditions. As the result, it can be seen that the shearing mode interlaminar fracture toughness is larger than the opening mode interlaminar fracture toughness. This means that interlaminar cracked specimen is tougher in shear loading condition and weaker in tensile loading condition.     

Fracture behaviour of PC/ABS resin under mixed-mode loading

Fracture behaviour of polycarbonate (PC)/acrylonitrile-butadiene-styrene (ABS) under mixed-mode loading conditions was studied for several weight fractions of PC and ABS. Mode I and mixed-mode fracture tests were carried out by using compact–tension–shear specimens. At a certain value of mixed-mode loading ratio K _II / K _I a crack of the shear type will initiates at the initial crack tip. Fracture toughness increases under mixed-mode loading with an increase in the mode II component, whereas it reduces with the appearance of a shear-type fracture. Fracture toughness and the appearance of a shear-type fracture depends on the blending ratio of PC and ABS. The transition to shear-type fracture occurs at lower value of K _II / K _I for resins with higher fracture toughness.     

FATIGUE CRACK GROWTH BEHAVIOUR OF SM45C STEEL UNDER CYCLIC MODE I WITH SUPERIMPOSED STATIC MODE II LOADINGS

Abstract— The practical applications of studies related to constant amplitude mode I loading are somewhat limited, since mode I crack growth is often influenced by mode II (sliding mode) or mode III (tearing mode) in industrial situations. For these cases, criteria, rules, and laws have to be worked out and verified by experiments. However, it is very difficult to evaluate mixed-mode fatigue cracking due to crack surface interference, crack closure, crack branching, etc. This paper, which defines the length of a branched crack as an effective slant crack with a length equal to the distance between the two crack tips, explains the influences of crack surface interference by introducing concepts of adhesive wear and scrutinizes some related researches on mixed-mode crack growth behaviour. Additionally an effective stress intensity factor range is described which considers crack closure and crack surface interference and is verified with crack growth tests under mode I fatigue loading and cyclic mode I with a superimposed static mode II loading.     

CRACK PROPAGATION UNDER MIXED MODE (I + III) LOADING

Abstract— In this paper are presented the results of fatigue crack propagation tests on angled-slit, three point bend mixed-mode (I + III) specimens manufactured from a low pressure steam turbine rotor forging. The path of crack propagation has been studied for two mixed mode (I + III) loading conditions. It has been observed that crack growth occurs by a mode I mechanism and a model has been developed to correlate crack growth rates in mixed mode (I + III) specimens with data from pure mode I fatigue tests.     

Experimental and numerical investigations on the influence of the loading direction on the fatigue crack growth

During a service loading fatigue cracks can be subjected to a mixed mode loading if, due to the alteration of the loading direction, the basic crack modes (Modes I, II and III) are combined. An alteration of the loading direction, e.g. can occur either occasionally paired with an overload (mixed mode overload) or permanently in terms of a mixed mode block loading as a combination of normal and shear stresses.Within the scope of this paper, experimental investigations on both mixed mode overloads, which are interspersed into a Mode I baseline level loading, and mixed mode block loadings are presented. The experimental investigations show that the retardation effect decreases with an increasing amount of Mode II of the overload. Due to the block loading, the fatigue crack growth rate is retarded as well, and the crack is also deflected. The kinking angle depends on the fraction of shear stresses. Furthermore, a detailed elastic–plastic finite element analysis of the fatigue crack growth after mixed mode overloads is presented in order to understand the mechanism of the load interaction effects. By such numerical simulations, it can be shown that, due to mixed mode overloads, plastic deformations occur, which on the one hand reduce the near-tip closure and on the other hand cause a far-field closure. Also the stress distribution before and after the crack tip changes. A mixed mode overload causes lower closure and the crack tip deformations become asymmetrical, which is a reason for the smaller retardation effect of a mixed mode overload.     

Decomposition of Eshelby’s energy momentum tensor and application to path and domain independent integrals for the crack extension force of a plane circular crack in Mode III loading

Vanishing divergence of Eshelby’s (energy momentum) tensor allows formulation of path or domain independent integral expressions of the crack extension force. In this work, a decomposition scheme of this tensor is presented, which results in zero divergence decomposed parts that allow formulation of expressions yielding the Mode I, II and III crack tip parameters J and K, with particular emphasis on Mode III, at present. By using the Mode III decomposed part of Eshelby’s tensor and the virtual crack extension method, a path and a domain independent integral, both new, for the crack extension force of a plane circular crack in axi-symmetric Mode III loading, are derived as examples of application.