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.     

An insight into mode II fracture toughness testing using SCB specimen

The effect of friction forces between the test specimen and its bottom supports on the mode II fracture toughness values obtained using the semicircular bend (SCB) specimen is investigated. First, a number of experiments were conducted on SCB specimen in order to determine the mode II fracture toughness of polymethyl methacrylate (PMMA) according to the conventional approaches available in the literature. Three different types of supports that have been frequently employed by researchers in recent years were used to evaluate the effect of support type on the fracture loads. It was found that the friction forces between the supports and the SCB specimen have a significant effect on the value of mode II fracture toughness measured using the SCB samples. Then, the specimen was simulated using finite element method for more detailed investigation on the near crack tip stress field evolution when friction forces increase between the supports and the SCB specimen. The finite element results confirmed that the type of support affects not only the stress intensity factors K_I and K_II but also the T‐stress. The experimental and numerical results showed that the use of the crack tip parameters available in literature for frictionless contact between the supports and the SCB specimen can result in significant errors when the mode II experiments are performed by using the fixed or roller‐in‐grove types of supports.     

An improved semi-circular bend specimen for investigating mixed mode brittle fracture

An edge cracked semi-circular specimen subjected to asymmetric three-point bend loading was suggested for investigating mixed mode fracture in brittle materials. Using finite element analysis, the crack parameters were obtained for various crack lengths and different locations of loading points. It was shown that by selecting appropriate positions for the loading points, full mode mixities from pure mode I to pure mode II could be achieved. Then, a series of fracture tests were conducted on PMMA using the proposed specimen. Very good agreement was found between the experimental results and those predicted from the generalized maximum tangential stress criterion.     

Analysis of a new specimen for mixed mode fracture tests on brittle materials

Numerical and experimental studies were performed on a new fracture test configuration called the diagonally loaded square plate (DLSP) specimen. The mode I and mode II stress intensity factors were computed for different crack lengths and crack orientation angles using finite element analysis. The numerical results show that the DLSP specimen is able to provide pure mode I, pure mode II and any mixed mode loading conditions in between. Fracture experiments were also conducted on Plexiglas using the DLSP specimen. It is shown that the results obtained from the fracture tests are consistent very well with mixed mode fracture theories.     

On the use of Brazilian disc specimen for calculating mixed mode I-II fracture toughness of rock materials

The centrally cracked Brazilian disc specimen has been used frequently in the past for investigating mixed mode I-II fracture toughness in rock materials. However, a review of the available test results reveals that the conventional fracture criteria like the maximum tangential stress criterion always underestimate the mixed mode I-II fracture toughness data obtained from the Brazilian disc specimen. In this paper, a generalized maximum tangential stress criterion which takes into account the effects of the three fracture parameters K_I, K_II and T-stress is used for predicting the mixed mode fracture toughness data available in the literature for several types of rock materials tested with the Brazilian disc specimen. It is shown that the generalized maximum tangential stress criterion provides significantly improved predictions for the experimental results.     

A material model to reproduce mixed‐mode fracture in concrete

This paper presents a material model to reproduce crack propagation in cement‐based material specimens under mixed‐mode loading. Its numerical formulation is based on the cohesive crack model, proposed by Hillerborg, and extended for the mixed‐mode case. This model is inspired by former works by Gálvez et al but implemented for its use in a finite element code at a material level, that is to say, at an integration point level. Among its main features, the model is able to predict the crack orientation and can reproduce the fracture behaviour under mixed‐mode fracture loading. In addition, several experimental results found in the literature are properly reproduced by the model.     

Determining the double‐K fracture parameters for three‐point bending notched concrete beams using weight function

Parameters of universal form of weight functions having four terms and five terms are derived for edge cracks in finite width of plate. The standard Tada Green's function is taken as the basis for the derivation. The shape of universal form of weight functions considered enables closed form expressions for cohesive toughness of three‐point bending test geometry of notched concrete beams due to linear cohesive stress distribution in the fictitious fracture zone. This solution provides a viable method to determine the double‐K fracture parameters: the initiation toughness, and the unstable toughness for mode I fracture of concrete beam. A comparison with existing analytical method shows that the weight function method for determination of the double‐K fracture parameters yields results without any appreciable error. The use of weight function will not only simplify the calculation to obtain the double‐K fracture parameters, and but also it will avoid the need of skilled numerical integration technique due to singularity problem at the integral boundary.     

The role of first non‐singular stress terms in mixed mode brittle fracture of V‐notched components: an experimental study

This paper investigates the effects of the first non‐singular stress terms on the fracture assessment of sharp V‐notches under mixed mode loading. First, numerical studies have been performed on a fracture test configuration called single V‐notched ring (SVR) specimen. Then, the notch stress intensity factors as well as the coefficients of the first non‐singular stress terms, which are vital parameters in brittle fracture of V‐notched components, were calculated via a finite element over‐deterministic algorithm for a wide range of loading and geometry conditions. The obtained results demonstrate that the SVR specimen is able to provide a complete range of mode mixities from pure mode I to pure mode II loading conditions. The numerical results, next, have been converted to dimensionless parameters and are illustrated in several graphs. Indeed, these graphs can be easily employed by the engineers for rapid calculation of the corresponding notch stress intensity factors and the coefficients of the first non‐singular stress terms in the SVR specimen. The obtained fracture parameters are then submitted to the maximum tangential stress criterion to assess the effects of the first non‐singular terms on fracture behaviour of the specimen. Finally, an experimental study has been performed on the SVR specimen made of Nayriz Marble rock for two notch angles with a complete range of mode mixities. The obtained experimental data confirm the significant role of the first non‐singular stress terms. In fact, these results show that considering only the singular stress terms may induce an average error of 38% in the predicted fracture loads, which can be decreased to about 12% just by adding the contribution of the first non‐singular terms to the maximum tangential stress criterion.     

Coupled meshfree and fractal finite element method for mixed mode two‐dimensional crack problems

This paper presents a coupling technique for integrating the element‐free Galerkin method (EFGM) with the fractal finite element method (FFEM) for analyzing homogeneous, isotropic, and two‐dimensional linear‐elastic cracked structures subjected to mixed‐mode (modes I and II) loading conditions. FFEM is adopted for discretization of the domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprise both the EFG and the finite element (FE) shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no special enriched basis functions or no structured mesh with special FEs are necessary and no post‐processing (employing any path independent integrals) is needed to determine fracture parameters, such as stress‐intensity factors (SIFs) and T‐stress. The numerical results show that SIFs and T‐stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also, a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions. A numerical example on mixed‐mode condition is presented to simulate crack propagation. As in the proposed coupled EFGM–FFEM at each increment during the crack propagation, the FFEM mesh (around the crack tip) is shifted as it is to the new updated position of the crack tip (such that FFEM mesh center coincides with the crack tip) and few meshless nodes are sprinkled in the location where the FFEM mesh was lying previously, crack‐propagation analysis can be dramatically simplified. Copyright © 2010 John Wiley & Sons, Ltd.     

On the analysis of a mixed mode bending sandwich specimen for debond fracture characterization

The mixed mode bending specimen originally developed for mixed mode delamination fracture characterization of unidirectional composites has been extended to the study of debond propagation in foam cored sandwich specimens. The compliance and strain energy release rate expressions for the mixed mode bending sandwich specimen are derived based on a superposition analysis of solutions for the double cantilever beam and cracked sandwich beam specimens by applying a proper kinematic relationship for the specimen deformation combined with the loading provided by the test rig. This analysis provides also expressions for the global mode mixities. An extensive parametric analysis to improve the understanding of the influence of loading conditions, specimen geometry and mechanical properties of the face and core materials has been performed using the derived expressions and finite element analysis. The mixed mode bending compliance and energy release rate predictions were in good agreement with finite element results. Furthermore, the numerical crack surface displacement extrapolation method implemented in finite element analysis was applied to determine the local mode mixity at the tip of the debond.     

Effect of compressive transverse normal stress on mode II fracture toughness in polymeric composites

Effect of transverse normal stress on mode II fracture toughness of unidirectional fiber reinforced composites was studied experimentally in conjunction with finite element analyses. Mode II fracture tests were conducted on the S2/8552 glass/epoxy composite using off-axis specimens with a through thickness crack. The finite element method was employed to perform stress analyses from which mode II fracture toughness was extracted. In the analysis, crack surface contact friction effect was considered. It was found that the transverse normal compressive stress has significant effect on mode II fracture toughness of the composite. Moreover, the fracture toughness measured using the off-axis specimen was found to be quite different from that evaluated using the conventional end notched flexural (ENF) specimen in three-point bending. It was found that mode II fracture toughness cannot be characterized by the crack tip singular shear stress alone; nonsingular stresses ahead of the crack tip appear to have substantial influence on the apparent mode II fracture toughness of the composite.     

A hybrid approach to determine fracture resistance for mode I and mixed‐mode I and II fracture specimens

This paper proposes a hybrid approach to determine the fracture resistance for mode I and mixed‐mode I and II fracture specimens, combining both numerically computed and experimentally measured load (P) versus load‐line displacement (LLD or Δ) relationships for metallic fracture specimens. The hybrid approach predicates on the same principle as the conventional, multiple‐specimen experimental method in determining the energy release rate. The hybrid method computes the P–Δ curves from multiple finite element (FE) models, each with a different crack depth. The experimental procedure measures the P–Δ curve from a standard fracture specimen with a growing crack. The intersections between the experimental P–Δ curve and the numerical P–Δ curves from multiple FE models dictate the LLD levels to compute the strain energy (U) using the area under the numerical P–Δ curves. This method provides accurate estimates of the J resistance data for both SE(B) specimen under mode I loading and single‐edge notched specimens under mixed‐mode I and II loading.     

A computational strategy for damage‐tolerant design of hollow shafts under mixed‐mode loading condition

Three‐dimensional numerical analyses, using the finite element method (FEM), have been adopted to simulate fatigue crack propagation in a hollow cylindrical specimen, under pure axial or combined axial‐torsion loading conditions. Specimens, made of Al alloys B95AT and D16T, have been experimentally tested under pure axial load and combined in‐phase constant amplitude axial and torsional loadings. The stress intensity factors (SIFs) have been calculated, according to the J‐integral approach, along the front of a part through crack, initiated in correspondence of the outer surface of a hollow cylindrical specimen. The crack path is evaluated by using the maximum energy release rate (MERR) criterion, whereas the Paris law is used to calculate crack growth rates. A numerical and experimental comparison of the results is presented, showing a good agreement in terms of crack growth rates and paths.     

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%.     

Error estimation of stress intensity factors for mixed‐mode cracks

This paper presents an a posteriori error estimator for mixed‐mode stress intensity factors in plane linear elasticity. A surface integral over an arbitrary crown is used for the separate calculation of the combined mode's stress intensity factors. The error in the quantity of interest is based on goal‐oriented error measures and estimated through an error in the constitutive relation. Copyright © 2006 John Wiley & Sons, Ltd.     

On the measurement and physical meaning of the cleavage fracture stress in steel

Elastic-plastic two-dimensional (2D) and three-dimensional (3D) finite element models (FEM) are used to analyze the stress distributions ahead of notches of four-point bending (4PB) and three-point bending (3PB) specimens with various sizes of a C-Mn steel. By accurately measuring the location of the cleavage initiation sites, the local cleavage fracture stress _f and the macroscopic cleavage fracture stress _F is accurately measured. The _f and _F measured by 2D FEM are higher than that by 3D FEM. _f values are lower than the _F, and the _f values could be predicted by _f=(0.8––1.0)_F. With increasing specimen sizes (W,B and a) and specimen widths (B) and changing loading methods (4PB and 3PB), the fracture load P _f changes considerably, but the _F and _f remain nearly constant. The stable lower boundary _F and _f values could be obtained by using notched specimens with sizes larger than the Griffiths–Owen specimen. The local cleavage fracture stress _f could be accurately used in the analysis of fracture micromechanism, and to characterize intrinsic toughness of steel. The macroscopic cleavage fracture stress _F is suggested to be a potential engineering parameter which can be used to assess fracture toughness of steel and to design engineering structure.     

A mode II weight function for subsurface cracks in a two-dimensional half-space

The general properties of a mode II Weight Function for a subsurface crack in a two‐dimensional half‐space are discussed. A general form for the WF is proposed, and its analytical expression is deduced from the asymptotic properties of the displacements field near the crack tips and from some reference cases obtained by finite elements models. Although the WF has general validity, the main interest is on its application to the study of rolling contact fatigue: its properties are explored for a crack depth range within which the most common failure phenomena in rolling contact are experimentally observed, and for a crack length range within the field of short cracks. The accuracy is estimated by comparison with several results obtained by FEM models, and its validity in the crack depth range explored is shown.     

Tangential strain‐based criteria for mixed‐mode I/II fracture toughness of cement concrete

Experimental and theoretical works are performed on the mixed‐mode I/II brittle fracture of cement concrete tested by edge cracked semicircular bend specimens. Theoretical background of the traditional fracture criteria including strain energy density, maximum tangential stress, and maximum tangential strain (MTSN) are introduced. The ability of each fracture criterion in prediction of the fracture test data is investigated. The comparison between the evaluations by the traditional criteria and the experimental data shows that none of them are capable of successfully estimating the fracture resistance of cement concrete. An enhanced version of the MTSN criterion is then employed to predict the test data. It is demonstrated that the extended MTSN criterion can successfully predict the test data in a higher accuracy than traditional criteria.     

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.     

Determination of T-stress from finite element analysis for mode I and mixed mode I/II loading

The elastic T-stress has been recognised as a measure of constraint around the tip of a crack in contained yielding problems. A review of the literature indicates that most methods for obtaining T are confined to simple geometry and loading configurations. This paper explores direct use of finite element analysis for calculating T. It is shown that for mode I more reliable results with less mesh refinement can be achieved if crack flank nodal displacements are used. Methods are also suggested for calculating T for any mixed mode I/II loading without having to calculate stress intensity factors. There is good agreement between the results from the proposed methods and analytical results. T-stress is determined for a test configuration designed to investigate brittle and ductile fracture in mixed mode loading. It is shown that in shear loading of a cracked specimen T vanishes only when a truly antisymmetric field of deformation is provided. However this rarely happens in practice and the presence of T in shear is often inevitable. It is shown that for some cases the magnitude of T in shear is much more than that for tension. The effect of crack length is also investigated. This revised version was published online in August 2006 with corrections to the Cover Date.