On the I–II mixed mode fracture of granite using four‐point bend specimen

Four‐point bend experiments on black granite are conducted. The fracture behaviours of granite under pure mode I, pure mode II and I–II mixed mode are investigated, and the corresponding stress intensity factors K_I , K_II and the non‐singular term T‐stress are obtained through numerical–experimental method. The results are compared with the theoretical predictions of generalized maximum tangential stress criterion and other conventional criteria. It shows that generalized maximum tangential stress criterion fits the experimental results better for considering the effect of T‐stress. Contrasting with other loading configurations, the values of T‐stress for asymmetric four‐point bend specimens are much smaller, especially for pure mode II specimens, which provide an asymmetric deformation field where the T‐stress is approaching zero.     

Finite element analysis of a subsurface crack

The problem of a subsurface crack parallel to the surface of a half space was studied by the finite element method. Without using the interface or gap elements over the crack faces, the crack faces would penetrate into each other for the traction-free boundary condition under shear loading, which is physically impossible. Using the gap elements, this problem was avoided, and a contact zone was observed near one crack tip. The size of the contact zone decreases but the maximum contact pressure at the closed crack tip increases as the crack approaches the surface. For tensile and shear loadings, both K _I (mode I stress intensity factor) and K _II (mode II stress intensity factor) increase as the crack approaches the surface. For shear loading there is no K _I at the closed tip and the K _I and K _II at the open tip are comparable as the crack approaches the surface.     

Measurement of mode I and mode II rock dynamic fracture toughness with cracked straight through flattened Brazilian disc impacted by split Hopkinson pressure bar

In order to find an effective and convincing method to measure rock dynamic fracture toughness for mode I and mode II, cracked straight through flattened Brazilian disc specimens of marble, which were geometrically similar for three size, were diametrically impacted by split Hopkinson pressure bar on the flat end of the specimen with three load angle respectively. History of stress intensity factors (K_I(t) for opening mode I, and K_II(t) for sliding mode II), mode mixture ratio (K_I(t)/K_II(t)), as well as mode I and mode II dynamic fracture toughness at crack initiation (K_Id and K_IId) were determined with the experimental–numerical method. It is found that there is a unique size effect for dynamic fracture test with the specimens, the mode mixture ratio is not solely determined by load angle (the angle between load direction and crack line) as in the static loading; the pure mode II load angle is 19° for the ?50 mm specimen, however it is 10° for the ?130 mm and ?200 mm specimens; the mode II load angle decreases with increment of specimen size. Realization of pure mode II is justified by the mode mixture ratio approaching zero, it can be realized under certain load angle and loading rate for the specimen of specified size. K_IId is generally greater than K_Id. Both K_Id and K_IId increase with increment of specimen size, and this trend for K_IId is more remarkable than that for K_Id.     

Mixed-mode fatigue crack growth behaviour in aluminium alloy

Fatigue crack propagation tests in compact mixed-mode specimens were carried out for several stress intensity ratios of mode I and mode II, K_I/K_II, in AlMgSi1-T6 aluminium alloy with 3 mm thickness. The tests were performed in a standard servo-hydraulic machine. A linkage system was developed in order to permit the variation of the K_I/K_II ratio by changing the loading angle. Crack closure loads were obtained through the compliance technique. A finite element analysis was also done in order to obtain the K_I and K_II values for the different loading angles. Crack closure increases under mixed-mode loading conditions in comparison to mode-I loading due the friction between the crack tip surfaces. Moreover, the crack closure level increases with the K_I/K_II ratio decrease. Correlations of the equivalent values of the effective stress intensity factor with the crack growth rates are also performed. Finally, an elastic–plastic finite element analysis was performed to obtain the plastic zones sizes and shapes and model the effect of mixed-mode loading on crack closure.     

Mixed mode crack tip parameters for different wheel positions relative to a vertical crack at the rail foot

Finite element method is used to analyze a rail with a vertical bottom up crack at its foot, under the axle load and surface traction of a wheel. The possibility of crack formation at the foot of the rail in the neighborhood of a welding connection is discussed. A brief review on the importance of T‐stress in brittle fracture is presented. Seven cases with different locations of the crack relative to rail's sleeper contact region are considered. Numerous positions of the wheel are considered, and in each case, 3 crack parameters K_I, K_II, and T‐stress are calculated. Then, the biaxiality ratio and the mixity parameter for each loading and crack condition are calculated. It is shown that the location of crack and wheel can create mixed mode loading in the cracked rail and that the magnitude of crack tip parameters are strongly dependent on these geometric variables. In particular, the magnitudes of T‐stress and biaxiality ratio are significant in some cases. The effect of friction between the crack faces in the presence of compressive mode I loading on the mode II stress intensity factor is studied. Under mixed mode loading, due to the axle load and surface traction, the most critical condition is the formation of vertical cracks near the sleeper contact region.     

Mixed mode fatigue and fracture in planar geometries: Observations on Keq and crack path modelling

The mixed mode I‐II fatigue and fracture is briefly reviewed, addressing experimental and numerical modelling aspects, and focusing on planar specimens. One major challenge concerns the determination of equivalent stress intensity factor (K_eq) in mixed mode situations. Several approaches were compared through the determination of K_eq/K_I over a wide range of values of K_I/K_II or K_II/K_I. Whereas all different approaches converge to the same value as K_I/K_II increases, the same does not happen for large K_II/K_I, where differences between values of K_eq persist. In the regions of 0 < K_I/K_II < 2 and 0 < K_II/K_I < 2, no stable trend of results can be defined. Experimental fatigue crack growth results are presented for Al alloy AA6082‐T6. Compact tension specimens, modified with holes, and four‐point bending specimens under asymmetrical loading promoting mixed mode situations, were subjected to fatigue crack growth tests, where crack path and crack growth rate were measured. The presentation of the fatigue crack growth data was made using a Paris law based upon K_eq. Differences in the Paris law constants were found for the different K_eq criteria. Recent developments in numerical techniques, as the implementation of the extended finite element method (XFEM) in finite element software packages allows to determine accurately crack paths in mixed mode fracture. This article highlights concepts for mixed‐mode fatigue and fracture and supporting data, identifying challenges still to be overcome.     

Mode I and mixed mode fracture of polysilicon for MEMS

An experimental study was carried out to investigate the local and effective fracture behaviour of polycrystalline silicon for microelectromechanical systems (MEMS). The apparent mode I critical stress intensity factor was determined from MEMS‐scale tension specimens containing atomically sharp edge pre‐cracks, while local deformation fields were recorded near the crack tip, with high resolution by the in situ Atomic Force Microscopy (AFM)/Digital Image Correlation (DIC) method previously developed by this group. The effective mode I critical stress intensity factor varied in the range 0.843–1.225 MPa√m. This distribution of values was attributed to local (in grain) cleavage anisotropy and to enhanced grain boundary toughening. The same sources resulted in very different local and macroscopic (apparent) stress intensity factors, which, combined with the small grain size of polysilicon (0.3 μm,) were the reason for subcritical crack growth that was evidenced experimentally by AFM topographic and AFM/DIC displacement measurements. The effect of local in‐grain anisotropy and granular inhomogeneity was stronger under mixed mode loading of edge cracks inclined at angles up to 55° with respect to the applied far‐field load. The K_I–K_II locus was characterized by scatter in the K_Ic values but on average it followed the curves calculated by the maximum tensile stress and the maximum energy release rate criteria calculated assuming isotropy.     

Fatigue crack growth under non-proportional mixed-mode loading in ferritic-pearlitic steel

Mode II fatigue crack growth tests as well as tests in sequential mode I and then mode II were performed on ferritic‐pearlitic steel. For ΔK_II ranging from 7 to , bifurcation occurs after 12–450 μm of coplanar growth at a decreasing speed. By contrast, hundreds of micrometres of constant speed coplanar growth were obtained under sequential mode I and then mode II loading, for and ΔK_I ranging from 0.25 to 1.0 ΔK_II . The crack growth rate is a simple sum of the contributions of each mode for ΔK_I= 0.25 ΔK_II but above this value a synergetic effect is found. The mechanism of this fast‐propagation mode is discussed in the light of strain range maps ahead of the crack tip obtained by digital SEM image correlation and elastic–plastic finite element calculations. The stability of the crack path according to the maximum growth rate criterion is demonstrated.     

The determination of the critical stress intensity factor in mode II loading and the shear fracture strength of pharmaceutical powder specimens

The shear fracture strength and the critical stress intensity factor in mode II loading of lactose monohydrate and acetylsalicylic acid powder compacts has been evaluated. The experimental results of the shear fracture strength and the critical stress intensity factor in mode II loading appeared to be in good agreement with powder behaviour such as lamination and capping during compaction. Values for the critical stress intensity factor in mode II loading depended on the depth of the crack and hence, any reference of such values or their use to calculate a fracture toughness ratio (K _IC ^I/K _IC ^II) must refer to the notch depth applied. The results confirmed that the failure of such powder compacts occurs mainly in tension, but that lactose monohydrate has a tendency also to fail in shear. The latter does not apply to acetylsalicylic acid. Hence, lactose monohydrate should only be used cautiously in layer or press-coated tablets.     

Theoretical analysis of the effects of relative crack length and loading angle on the experimental results for cracked Brazilian disk testing

Based on the closed-form expressions of the stress intensity factors for a central cracked circular disk, four error transfer functions of the stress intensity factors were introduced in order to analyze the effect of the relative crack length and the error of loading angle on the experimental results for Brazilian disk testing. The analyzed results show that the precision of are relevant on the relative crack length and the error of loading angle. Further analysis show that the error of loading angle Δθ has a significant effect on the errors of under mixed-mode loading condition, but Δθ has nearly no effect on the error of K_I under pure mode I loading condition and smaller effect on that of K_II under pure mode II loading condition. Finally, the recommended range of the relative crack length is between 0.4 and 0.6.     

State‐of‐the‐art review on extended stress intensity factor concepts

The stress intensity factor concept for describing the stress field at pointed crack or slit tips is well known from fracture mechanics. It has been substantially extended since Williams' basic contribution (1952) on stress fields at angular corners. One extension refers to pointed V‐notches with stress intensities depending on the notch opening angle. The loading‐mode‐related simple notch stress intensity factors K₁, K₂ and K₃ are introduced. Another extension refers to rounded notches with crack shape or V‐notch shape in two variants: parabolic, elliptic or hyperbolic notches (‘blunt notches’) on the one hand and root hole notches (‘keyholes’ when considering crack shapes) on the other hand. Here, the loading‐mode‐related generalised notch stress intensity factors K_1ρ, K_2ρ and K_3ρ are defined. The concepts of elastic stress intensity factor, notch stress intensity factor and generalised notch stress intensity factor are extended into the range of elastic–plastic (work‐hardening) or perfectly plastic notch tip or notch root behaviour. Here, the plastic notch stress intensity factors K_1p, K_2p and K_3p are of relevance. The elastic notch stress intensity factors are used to describe the fatigue strength of fillet‐welded attachment joints. The fracture toughness of brittle materials may also be evaluated on this basis. The plastic notch stress intensity factors characterise the stress and strain field at pointed V‐notch tips. A new version of the Neuber rule accounting for the influence of the notch opening angle is presented.     

Fracture Toughness and Creep Crack Growth in Polypropylene Under Combined Modes I and II Loading

Fracture toughness and creep crack growth characteristics under combined mode I and II loadings were studied using the compact tension shear (CTS) specimens of polypropylene. The K _I - K _II envelope for crack initiation was obtained under various combined mode loadings. The creep crack growth rates da/dt under combined mode I and mode II loadings can be correlated with a single effective stress intensity factor K _Ieff based on the combined mode fracture envelope.     

Cracked asphalt pavement under traffic loading – A 3D finite element analysis

An asphalt pavement containing a transverse top-down crack is investigated under traffic loading using 3D finite element analysis. The stress intensity factors (SIFs) and the T-stress are calculated for different distances between the crack and the vehicle wheels. It is found that all the three Modes (I, II and III) are present in the crack deformation. The signs and magnitudes of K_I, K_II, K_III and T are significantly dependent on the location of the vehicle wheels with respect to the crack plane. The magnitude of T-stress is considerable, if compared to the stress intensity factors, when one of the wheels is very close to the crack plane.     

Spannungsfaktoren eines Risses in der Umgebung eines Kreisloches und ihr Einfluß auf das Bruchverhalten

Stress Intensity Factors in the Neighbourhood of a Circular Hole and their Influence on Crack Behavior . A photoelastic method was developed to determine the stress intensity factors K_I and K_II for cracks subjected to mixed-mode loading. The constants in the near field expansion about the crack tip were computed using a non-linear optimization program to give a best fit to the observed isochromatics. Copying the latter onto an equal density film increased their sharpness and, thus, the accuracy of the determination. The method was applied to cracks lying perpendicular to the external stress near a circular hole in a plate under uniaxial tension and the results used to describe the paths of cracks in the neighbourhood of a hole.     

Short fatigue crack growth behavior under mixed-mode loading

Mixed-mode loading represents the true loading condition in many practical situations. In addition, most of the fatigue life of many components is often spent in the short crack growth stage. The study of short crack growth behavior under mixed-mode loading has, therefore, much practical significance. This work investigated short crack growth behavior under mixed-mode loading using a common medium carbon steel. The effects of load mixity, crack closure, and load ratio on short crack growth behavior were evaluated by conducting experiments using four-point bending specimens with several initial K _II /K _I mixed-mode ratios and two load ratios. Cracks were observed to grow along the paths with very small K _II /K _I ratios (i.e. mode I). The maximum tangential stress criterion was used to predict the crack growth paths and the predictions were found to be close to the experimental observations. Several parameters including equivalent stress intensity factor range and effective stress intensity factor range were used to correlate short crack growth rates under mixed-mode loading. Threshold values for short cracks were found to be lower than those for long cracks for all the mixed-mode loading conditions. Crack closure was observed for the entire crack length regime with all load mixity conditions at R ≈ 0.05 and for short crack regime under high load mixity condition at R = 0.5. Several models were used to describe mean stress effects and to correlate crack growth rate data.     

Analytical and Numerical Solution of the Stress Field and the Dynamic Stress Intensity Factors in a Cracked Plate under Sinusoidal Loading

A cracked plate subjected to a sinusoidal loading perpendicular to its plane is considered, and the analytical solution of the dynamic vibration behavior of a plate, which allowed the determination of the stress field near the crack tip, is developed. A mixed mode of loading near the crack tip has been established and described with dynamic stress intensity factors K _I (z,t) and K _II (z,t) associated with modes I and II crack openings, respectively. To validate the analytical results, a finite element analysis (FEA) of a 1 × 1 m square plate with a thickness of 1 cm, having a middle crack of 10 cm in length, is made. The results have shown significant agreement between analytical and FEA findings.     

Experimental Evaluation of Mixed Mode Stress Intensity Factor for Prediction of Crack Growth by Phoelastic Method

Determination of stress intensity factors K _I, K _II, and constant stress term, σ _ox is investigated. A theory of determining the stress intensity factors using photo-elastic method is formulated taking three stress terms. Three-parameter method of fracture analysis for determining the mixed mode stress intensity factors under biaxial loading conditions from photo-elastic isochromatic fringe data is used. A special biaxial test rig is designed and fabricated for loading the specimen biaxially. A simplified and accurate method is proposed to collect the data from isochromatic fringes. Taking specimen geometry and boundary conditions into account, regression models are developed for estimation of fracture parameters.     

Mixed-mode Bueckner weight functions using boundary element analysis

A Bueckner singular displacement field is incorporated into the boundary element method in order to calculate weight functions for mixed mode stress intensity factor in a two-dimensional cracked body. Weight functions for modes I and II can be calculated independently and hence the factors K _I and K _II can be obtained for any arbitrary loading on any boundary. Results are obtained for some slant crack problems in finite sheets and compared with known results where available.

---

Résumé On incorpore à la méthode des éléments dans un contour un champ singulier de déplacements de Bueckner en vue de calculer les fonctions pondérales qui caractérisent le facteur d'intensité pour un mode mixte dans un corps bidimensionnel fissuré. Les fonctions pondérales relatives aux Modes I et II peuvent être calculées indépendamment. Dès lors, les facteur K _I et K _II peuvent être obtenus pour toute charge arbitraire agissant sur un contour quelconque. On obtient des résultats pour divers problèmes de fissures et on les compare à des résultats connus, lorsque ceux-ci sont disponibles.

     

Interface crack in periodically layered bimaterial composite

A directional crack growth prediction in a compressed homogenous elastic isotropic material under plane strain conditions is considered. The conditions at the parent crack tip are evaluated for a straight stationary crack. Remote load is a combined biaxial compressive normal stress and pure shear. Crack surfaces are assumed to be frictionless and to remain closed during the kink formation wherefore the mode I stress intensity factor K _I is vanishing. Hence the mode II stress intensity factor K _II remains as the single stress intensity variable for the kinked crack. An expression for the local mode II stress intensity factor k ₂ at the tip of a straight kink has been calculated numerically with an integral equation using the solution scheme proposed by Lo (1978) and refined by He and Hutchinson (1989). The confidence of the solution is strengthened by verifications with a boundary element method and by particular analytical solutions. The expression has been found as a function of the mode II stress intensity factor K _II of the parent crack, the direction and length of the kink, and the difference between the remote compressive normal stresses perpendicular to, and parallel with, the plane of the parent crack. Based on the expression, initial crack growth directions have been suggested. At a sufficiently high non-isotropic compressive normal stress, so that the crack remains closed, the crack is predicted to extend along a curved path that maximizes the mode II stress intensity factor k ₂. Only at an isotropic remote compressive normal stress the crack will continue straight ahead without change of the direction. Further, an analysis of the shape of the crack path has revealed that the propagation path is, according the model, required to be described by a function y=cx , where the exponent is equal to 3/2. In that case, when =3/2, predicts the analytical model a propagation path that is self-similar (i.e. the curvature c is independent of any length of a crack extension), and which can be described by a function of only the mode II stress intensity factor K _II at the parent crack tip and the difference between the remote compressive normal stress perpendicular to, and parallel with, the parent crack plane. Comparisons with curved shear cracks in brittle materials reported in literature provide limited support for the model discussed.     

History effect in fatigue crack growth under mixed-mode loading conditions

Plastic deformation within the crack tip region introduces internal stresses that modify subsequent behaviour of the crack and are at the origin of history effects in fatigue crack growth. Consequently, fatigue crack growth models should include plasticity-induced history effects. A model was developed and validated for mode I fatigue crack growth under variable amplitude loading conditions. The purpose of this study was to extend this model to mixed-mode loading conditions. Finite element analyses are commonly employed to model crack tip plasticity and were shown to give very satisfactory results. However, if millions of cycles need to be modelled to predict the fatigue behaviour of an industrial component, the finite element method becomes computationally too expensive. By employing a multiscale approach, the local results of FE computations can be brought to the global scale. This approach consists of partitioning the velocity field at the crack tip into plastic and elastic parts. Each part is partitioned into mode I and mode II components, and finally each component is the product of a reference spatial field and an intensity factor. The intensity factor of the mode I and mode II plastic parts of the velocity fields, denoted by dρ_I/dt and dρ_II/dt, allow measuring mixed-mode plasticity in the crack tip region at the global scale. Evolutions of dρ_I/dt and dρ_II/dt, generated using the FE method for various loading histories, enable the identification of an empirical cyclic elastic–plastic constitutive model for the crack tip region at the global scale. Once identified, this empirical model can be employed, with no need of additional FE computations, resulting in faster computations. With the additional hypothesis that the fatigue crack growth rate and direction can be determined from mixed-mode crack tip plasticity (dρ_I/dt and dρ_II/dt), it becomes possible to predict fatigue crack growth under I/II mixed-mode and variable amplitude loading conditions. To compare the predictions of this model with experiments, an asymmetric four point bend test system was setup. It allows applying any mixed-mode loading case from a pure mode I condition to a pure mode II. Initial experimental results showed an increase of the mode I fatigue crack growth rate after the application of a set of mode II overload cycles.