Finite element analysis of axial splits in composite Iosipescu specimens

In this paper, a finite element analysis of skew-symmetric splits along the fiber direction in unidirectional composite Iosipescu specimens is performed. The energy release rates G _I, G _II, and G _total associated with axial splits in cracked Iosipescu specimens under external biaxial loading conditions are computed by four different numerical schemes: displacement correlation, displacement extrapolation, J-integral, and the modified crack closure integral. Using beam theory analysis, an analytical solution for the energy release rates is also proposed. Axial splits in Iosipescu specimen propagate under mixed mode conditions, with G _I and G _II varying with the crack length a. For short and medium crack lengths G _I>G _II, while for long cracks, G _II is dominant. The energy release rates G _I, G _II, and G _total are strongly dependent on the biaxial type of loading. The G-estimates obtained by the modified crack closure integral schemes are found to be the most accurate among all the numerical schemes chosen in this study. In the analyses of axial splits in composite Iosipescu specimens, the displacement correlation and extrapolation techniques yielded poor results. For long crack lengths, the analytical results from the beam theory analysis are in fair agreement with those from the modified crack closure integral schemes; however, for short and medium crack lengths, there is a significant difference between the analytical and numerical results. In composite Iosipescu specimens, stable crack propagation (mode I dominant) can be achieved by increasing the tension/shear ratio in the external loading boundary conditions.     

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.     

J integral and COD fracture criteria in 40CrNiMo steel under mixed mode I+II loading

Abstract

By defining the J integral and crack opening displacement (COD) under mixed mode I + II loading, a mixed J integral fracture criterion is proposed and the relationship between the J integral and COD in 40CrNiMo steel is discussed. The mixed J integral J _M and its mode I and II components J _I and J _II were calculated by the finite element method, while the mixed COD and its mode I and II components CTOD and CTSD were measured using a duplicated grid. The critical values J _Mc and COD_c for mixed crack initiation were determined by a resistance curve. The results show that mode II loading lowers both J _Mc and COD_c for 40CrNiMo steel. The variation of J _Mcfrommode I tomode II loading is found to be in accordance with the linked equations J _Mc=J _Ii+J _IIi;(J _Ii/J _Ic)+(J _IIi/J _IIc)=1, where J _Icand J _IIcare the critical J integrals of pure mode I and II cracks,and J _Iiand J _IIiare the mode I and II components of J _Mc at arbitrary mixed K_I/K_II ratio respectively; the J _Mc value for a given K_I/K_II ratio can be obtained if J _Ic and J _IIc are known. Finally, under valid loads, J_M and the mixed COD satisfy the relation J _M=J _I+J _II=dσ₀CTOD+d _sτ₀CTSD. When unified by yield stress σ_y the relation becomes J _M=dσ_yCOD, where d _n, d _s and d are coefficients, and σ₀ and τ₀ are the tensile and shear stress at the crack tip strip respectively. While d _n and d _s vary with K_I/K_II ratio and materials, d was found to have a constant value of about 0.98.     

A fracture mechanics analysis on the fatigue behaviour of cruciform joints of duplex stainless steel

The paper presents the results of an experimental and numerical study on the fatigue behaviour of cruciform load carrying joints made from the duplex stainless steel and failing from the weld root through the weld metal. Fatigue crack growth (FCG) data, obtained in specimens of the weld metal, are presented, as well as threshold data, both obtained for R= 0 and 0.5. The influence of stress ratio is discussed, and the FCGR results are compared with data for low carbon structural steels. S–N data were obtained in the joints, both for R= 0.05 and 0.5, and the fatigue cracking mechanisms were analysed in detail with the SEM. It was found that the cracks propagated very early in the lifetime of the joints, under mixed mode conditions (I + II), but the mode I component was found to be predominant over mode II. The geometries of the cracks were defined in detail from measurements taken in the fracture surfaces. A 2D FE analysis was carried out for the mixed mode inclined cracks obtained at the weld root, and the J‐integral formulations were obtained as a function of crack length and crack propagation angle. The values of the crack propagation angle, θ_i, were obtained for the J_max conditions, and it was found that, in the fatigue tests, the cracks propagated in directions very close to the predicted directions of maximum J. K_I and K_II formulations were obtained, and the K_I data were compared with the formulations given in the PD6493 (BS7910) document, and some differences were found. A more general formulation for K under mixed mode conditions was derived. The derived K solutions were applied to predict the fatigue lives of the joints under crack propagation, and an extremely good agreement was found with the experimental results obtained in the fatigue tests.     

Effects of microstructure on mixed-mode, high-cycle fatigue crack-growth thresholds in Ti-6Al-4V alloy

Effect of microstructure on mixed‐mode (mode I + II), high‐cycle fatigue thresholds in a Ti‐6Al‐4V alloy is reported over a range of crack sizes from tens of micrometers to in excess of several millimeters. Specifically, two microstructural conditions were examined—a fine‐grained equiaxed bimodal structure (grain size ～20 µm) and a coarser lamellar structure (colony size ～500 µm). Studies were conducted over a range of mode‐mixities, from pure mode I (ΔK_II/ΔK_I = 0) to nearly pure mode II (ΔK_II/ΔK_I ～ 7.1), at load ratios (minimum load/maximum load) between 0.1 and 0.8, with thresholds characterized in terms of the strain‐energy release rate (ΔG) incorporating both tensile and shear‐loading components. In the presence of through‐thickness cracks—large (> 4 mm) compared to microstructural dimensions—significant effects of mode‐mixity and load ratio were observed for both microstructures, with the lamellar alloy generally displaying the better resistance. However, these effects were substantially reduced if allowance was made for crack‐tip shielding. Additionally, when thresholds were measured in the presence of cracks comparable to microstructural dimensions, specifically short (～200 µm) through‐thickness cracks and microstructurally small (< 50 µm) surface cracks, where the influence of crack‐tip shielding would be minimal, such effects were similarly markedly reduced. Moreover, small‐crack ΔG_TH thresholds were some 50–90 times smaller than corresponding large crack values. Such effects are discussed in terms of the dominant role of mode I behaviour and the effects of microstructure (in relation to crack size) in promoting crack‐tip shielding that arises from significant changes in the crack path in the two structures.     

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.     

Continuum shape sensitivity analysis of mixed-mode fracture using fractal finite element method

This paper presents a new fractal finite element based method for continuum-based shape sensitivity analysis for a crack in a homogeneous, isotropic, and two-dimensional linear-elastic body subject to mixed-mode (modes I and II) loading conditions. The method is based on the material derivative concept of continuum mechanics, and direct differentiation. Unlike virtual crack extension techniques, no mesh perturbation is needed in the proposed method to calculate the sensitivity of stress-intensity factors. Since the governing variational equation is differentiated prior to the process of discretization, the resulting sensitivity equations predicts the first-order sensitivity of J-integral or mode-I and mode-II stress-intensity factors, K_I and K_II, more efficiently and accurately than the finite-difference methods. Unlike the integral based methods such as J-integral or M-integral no special finite elements and post-processing are needed to determine the first-order sensitivity of J-integral or K_I and K_II. Also a parametric study is carried out to examine the effects of the similarity ratio, the number of transformation terms, and the integration order on the quality of the numerical solutions. Four numerical examples which include both mode-I and mixed-mode problems, are presented to calculate the first-order derivative of the J-integral or stress-intensity factors. The results show that first-order sensitivities of J-integral or stress-intensity factors obtained using the proposed method are in excellent agreement with the reference solutions obtained using the finite-difference method for the structural and crack geometries considered in this study.     

A crack propagation criterion based on ΔCTOD measured with 2D‐digital image correlation technique

The fatigue cracks growth rate of a forged HSLA steel (AISI 4130) was investigated using thin single edge notch tensile specimen to simulate the crack development on a diesel train crankshafts. The effect of load ratio, R, was investigated at room temperature. Fatigue fracture surfaces were examined by scanning electron microscopy. An approach based on the crack tip opening displacement range (ΔCTOD) was proposed as fatigue crack propagation criterion. ΔCTOD measurements were carried out using 2D‐digital image correlation techniques. J‐integral values were estimated using ΔCTOD. Under test conditions investigated, it was found that the use of ΔCTOD as a fatigue crack growth driving force parameter is relevant and could describe the crack propagation behaviour, under different load ratio R.     

Effect of residual stresses on crack behaviour in single edge bending specimens

Residual stresses due to manufacturing processes, such as welding, change the load bearing capacity of cracked components. The effects of residual stresses on crack behaviour in single edge bending specimens were investigated using Finite element analyses. Three parameters (J, Q and R) were used to study the crack behaviour. The J‐integral predicts the size scale over which large stresses and strains exist, the constraint parameter Q describes the crack‐tip constraint as a result of geometry, loading mode and crack depth and the constraint parameter R is used to describe the constraint resulting from residual stresses. To carry out a systematic investigation on the effect of residual stresses on the J‐integral and crack‐tip constraints, models under different combinations of residual stresses and external loads with different crack depths were analysed. It has been shown that the crack‐tip constraint R increased by tensile residual stresses around the crack‐tip. On the other hand, the constraint parameter R decreased and tended to zero at high external load levels.     

Part-circular surface cracks in round bars under tension,bending and twisting

Circular-fronted cracks in round bars subject to tension, bending and twisting are considered. Numerical expressions are given allowing the calculation of stress intensity factors K _I, K _II, K _III at every point on the crack front for a wide range of crack geometries. Comparisons are made with analytical, experimental and numerical results abailable in the literature. Crack shapes satisfying the iso-K _I criterion are also determined, making it possible to investigate the problem of crack growth behaviour under tensile or bending fatigue loads.     

Three-dimensional fracture analysis of thin-film debonding

A simplified mixed-mode fracture analysis combining nonlinear thin-plate stress solutions with crack-tip elasticity results has been developed to account for local variations of G _I, G _II and G _III in thin-film debond problems associated with large film deformations. Membrane and bending stresses from the plate analysis are matched with the crack-tip singularity solution over a small boundary region at the crack tip where the effect of geometric nonlinearity is small. Local variations in each of the individual components of the energy release rate are directly related to the jump in these stresses across the crack border.Specific results are presented for 1-D and elliptical planeform cracks. Deformations were induced either by a pressure acting normal to the film surface or biaxial compression or tension stresses applied to the substrate in which the loading axes and debond axes coincide. The latter type of loading involves buckling of the delaminated film. The model predictions compare well with more rigorous solutions provided the film thickness is small compared to the debond dimensions. In all cases analyzed, G _III was negligible. The ratio G _I/G _II typically decreases with increasing load or film deformation, the rate was moderate for pressure loading while generally sharp for compression loading. Film-substrate overlap may occur for certain debond geometry and loading conditions. Prevention of this by the substrate may critically increase the energy available for crack propagation.     

J II fracture testing of a plastically deforming material

A compact model II fracture specimen was previously analyzed and employed to determine the mode II fracture toughness K _IIc , of perspex. In employing this specimen for a more ductile material such as aluminium, it was observed that the load vs. crack sliding displacement record becomes nonlinear for small loads. Thus, concepts of linear elastic fracture mechanics cannot be employed. To this end, the specimen was calibrated for J-integral testing, so that J _IIc mesurements can be performed.In this study, mode I and II tests are carried out on an aircraft aluminium alloy, AI 7075-T7351. First, standard K _Ic tests are performed leading to a value of 27.9 15-1 which would be equivalent to a J _Ic of 10.7 kN/m. Then standard J _Ic tests are carried out on this material with specimen thicknesses, of 5, 7.5 and 9.9 mm, leading to an average J _Ic value of 10.5 kN/m. Methods for J _II testing are proposed; a series of specimens of six thicknesses between 5 and 16 mm are employed for testing. An average J _IIc value was found to be 40.2 kN/m which yields a K _IIc value of 54.1 15-2. Thus, K _IIc is seen to be approximately twice that of K _Ic for this material.     

Mixed mode interface crack in a pure power-hardening bimaterial

Analytical and numerical analysis of the dominant singularity solutions of the stress and strain field near an interface crack in a pure power-hardening bimaterial indicates that the crack stress singularity is –1/(n _II+1) for hardening power of n _I and n _II(n _I<n _II). This result is obtained by solving the non-linear eigenvalue equations of the stress field near a plane crack while observing the conservation properties of the J-integral and the continuity conditions of the interface. Numerical results are presented for the distribution of stress, strain and displacement field of various mixed mode interface cracks when n _I=1 and n _II=5. Possible further destruction of the bimaterial with interface cracks is also discussed.

---

Résumé On procède à una analyse théorique et numérique des solutions de singularité dominantes dans les champs de contraintes et de dilatation voisins d'une fissure d'interface dans un complexe bimétallique répondant à une loi simple de durcissement parabolique. On montre que la singularité de la contrainte vaut –1/n _II+1 pour des modules d'écrouissage n _I et n _II(n _I<n _II).On obtient ce résultat en résolvant les équations non linéaires d'eigenvalues du champ de contraintes au voisinage d'une fissure plane, tout en respectant les propriétés de conservation de l'intégrale J, et les conditions de continuité numérique pour la distribution des champs de contraintes, de dilatations et de déplacements correspondant à diverses fissures d'interfaces sollicitées selon des modes mixtes, avec n _I=1 et n _II=5.On discute également d'une destruction ultérieure possible du bimatériau en présence de fissures d'interface.

     

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.     

Analysis of dynamic characteristics of through‐wall cracks between 2 boreholes in the directed fracture controlled blasting

Through‐wall cracks between 2 boreholes in the directed fracture controlled blasting have been always concerned by researchers. Dynamic characteristics of through‐wall cracks between 2 boreholes and lateral crack propagation of boreholes in the double‐borehole slot mode and the synchronous blasting of boreholes were examined using the explosive loading digital dynamic caustics experiment system. And the effects of borehole loading mode and borehole clearance on through‐wall cracks between boreholes were examined using distinct lattice spring model numerical analysis, based on the experiment model. Findings show that the tips of through‐wall cracks between boreholes did not meet directly but staggered, continuously propagated after meeting, and moved closer to the existing anisotropic crack direction. The velocity and acceleration of crack propagation fluctuated. K_I rapidly decreased from the maximum value, then gradually increased after a repeated volatility, and began to decrease after it reached the second peak. During the process of crack propagation, K_II was basically smaller than K_I. The dynamic energy release rate rapidly decreased from the maximum value, reached the second peak after the volatility, and gradually decreased again. The borehole loading mode and borehole clearance had significant effects on through‐wall cracks between boreholes.     

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.     

Comparison of mode I and mode II elastic-plastic fracture toughness for two low alloyed high strength steels

In the present work, mode I and mode II tests were carried out on two low alloyed high strength steels. An asymmetrical four point bend specimen and J _II-integral vs. crack growth resistance curve technique were used for determining the mode II elastic-plastic fracture toughness, J _IIc · J _II-integral expression of the specimen was calibrated by finite element method. The results indicate that the present procedure for determining the J _IIc values is easy to use. Moreover, the mode I fracture toughness J _Ic is very sensitive to the rolling direction of the test steels, but the mode II fracture toughness J _IIc is completely insensitive to the rolling direction of the steels, and the J _IIc /J _Ic ratio is not a constant for the two steels, including the same steel with different orientations. Finally, the difference of the fracture toughness between the mode I and mode II is discussed with consideration of the different fracture mechanisms.     

Effect of biaxial loads on elastic-plastic <Emphasis Type="Italic">J</Emphasis> and crack tip constraint for cracked plates: finite element study

Based on detailed two-dimensional (2-D) and three-dimensional (3-D) finite element (FE) analyses, this paper attempts to quantify in-plane and out-of-plane constraint effects on elastic-plastic J and crack tip stresses for a plate with a through-thickness crack and semi-elliptical surface crack under positive biaxial loading. For the plate with a through-thickness crack, plate thickness and relative crack length are systematically varied, whereas for the plate with a semi-elliptical surface crack, the relative crack depth and aspect ratio of the semi-elliptical crack are systematically varied. It is found that the reference stress based approach for uniaxial loading can be applied to estimate J under biaxial loading, provided that the limit load specific to biaxial loading is used, implying that quantification of the biaxiality effect on the limit load is important. Investigation on the effect of biaxiality on the limit load suggests that for relatively thin plates with small cracks, in particular with semi-elliptical surface cracks, the effect of biaxiality on the limit load can be neglected for positive biaxial loading, and thus elastic-plastic J for a biaxially loaded plate could be estimated, assuming that such plate is subject to uniaxial load. Regarding the effect of biaxiality on crack tip stress triaxiality, it is found that such effect is more pronounced for a thicker plate. For plates with semi-elliptical surface cracks, the crack aspect ratio is found to be more important than the relative crack depth, and the effect of biaxiality on crack tip stress triaxiality is found to be more pronounced near the surface points along the crack front.     

Fatigue cracks emanating from sharp notches in high-strength aluminium alloys: The effect of loading direction, kinking, notch geometry and microstructure

The effect of notch geometry on the propagation of fatigue cracks emanating from sharp V-shaped notches was investigated. To this purpose, an experimental campaign has been conducted on Al-7075–T651 specimens carrying notches with aperture angles of 45°, 90°, and 135°. In order to investigate the role of microstructure texture, specimens were extracted from the plates with the main axis either in the longitudinal rolling direction (L-samples) or in the transversal direction (T-samples), or 45° inclined with respect to both directions (LT-samples). The effect of stress amplitude was investigated by performing tests at two load levels. Three loading directions θ = 0°, 45° and 90° were considered. Some specimens experienced pure Mode I loading condition, whereas the remaining ones were subjected to combined Mode I and Mode II loading condition. The crack deflection induced by the variation in loading direction was determined by measuring the kinking angle. A linear elastic fracture mechanics approach was adopted for the analysis of experimental results. Stress intensity factors (SIF) of straight cracks were calculated using an appropriate weight function set up for studying inclined edge cracks emanating from sharp V-notches. On the contrary, a finite element model has been built up to derive the SIFs at the tip of the kinked cracks. The influence of K_II on the crack propagation was discussed on the basis of theoretical and semi-empirical models. It has been found that (i) the crack initiation at the notch root occurred in mixed mode conditions, (ii) a decreasing Mode II component with growing crack length was observed under initial loading direction θ₀ = 45° and θ₀ = 90°, (iii) a crack deflection was observed after 45° rotation of the initial loading direction; a good prediction of the kinking angle was obtained using the maximum tangential stress criterion, and (iv) a fairly good rationalization of all the collected crack growth rate data is obtained if the driving force for crack propagation is expressed in terms of K_I.     

Contour integral method for stress intensity factors of interface crack

A general Betti's reciprocal work theorem with interface cracks of a bimaterial is established in this paper, and a path independent contour integral method for the stress intensity factor (SIF) of the interface crack was obtained. When the stress and displacement fields in a specimen are calculated by the finite element method, the SIF K _I and K _II of interface cracks can be obtained immediately by a contour integral. Some solutions of interesting examples, such as two collinear interface cracks, are also given.Presented at the Far East Fracture Group (FEFG) International Symposium on Fracture and Strength of Solids, 4–7 July 1994 in Xi'an China.