An investigation of T‐stresses in a plate with an inclined crack under biaxial loadings

For plates with an inclined crack of wide‐range aspect ratios under biaxial loadings, T‐stress values are calculated with three‐dimensional finite element method. The results show that the normalized T‐stress is crack length and orientation dependent. A linear equation for the relationship between normalized T‐stresses and biaxility factors is proposed to describe the normalized T‐stresses for different crack lengths and crack angles under different biaxial loadings, which is more convenient and involves wider biaxility ratios compared with the existing solutions. The plate thickness effect and the trend of normalized T‐stresses along the crack front thickness are also studied for mode I and I–II mixed‐mode cracks. Based on the analyses and comparisons, it is necessary to take the thickness effect into consideration when the crack length is long enough (a/W = 7/10). When the component of mode II is significant (β > 45°), and the biaxility ratios are negative, T‐stresses near the free surface are lower than those at other positions, which are the opposite of mode I crack and most of I–II mixed‐mode crack.     

Buckling analysis of orthotropic rectangular plate resting on Pasternak elastic foundation under biaxial in-plane loading

In this study, buckling of rectangular orthotropic plates resting on a Pasternak elastic foundation under biaxial in-plane loading by the power series method (the method of Frobenius) was analyzed. Similar to many studies, two opposite edges of loading are simply supported and two other edges are assumed clamped. In order to extract the characteristic equations of orthotropic rectangular plate under in-plane loading resting on a Pasternak elastic foundation, the classical plate theory, by considering the interaction between plate and foundation, is used. The results showed that in the aspect ratio of less than 2, the existing Pasternak foundation caused the buckling load to increase severely, but by increasing the aspect ratio, the effect of the foundation is negligible. Applying the in-plane load in the y-direction caused the buckling load to decrease, but by increasing the aspect ratios the effect of the load is negligible.     

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

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

The effect of T‐stress on crack‐tip plastic zones under mixed‐mode loading conditions

In this paper, the influence of T‐stress on crack‐tip plastic zones under mixed‐mode I and II loading conditions is examined. The crack‐tip stress field is defined in terms of the mixed‐mode stress intensity factors and the T‐stress using William's series expansion. The crack‐tip stress field is incorporated into the Von Mises yield criteria to develop an expression that determines the crack‐tip plastic zone. Using the resultant expression, the plastic zone is plotted for various combinations of mode II to mode I stress intensity factor ratios and levels of T‐stress. The properties of the plastic zone affected by T‐stress and mixed‐mode phase angle are discussed. The observations obtained on plastic zones variations are important for further fatigue and fracture analyses for defects in engineering structures under mixed‐mode loading conditions.     

Role of plastic zone in crack growth direction criterion under mixed mode loading

An approach to determine the crack growth direction under mixed-mode loading conditions is presented. The plastic zone shape around the crack tip is applied for evaluating angle of crack propagation. It is proposed that a mixed-mode crack will extend along the plastic zone radius with a minimum value. The prediction of the proposed criterion is compared with the experimental data and other models. The agreement is fairly good.     

A method for the assessment of hyperstatic cracked structures in the elastic–plastic regime

This paper presents a methodology for the assessment of hyperstatic cracked structures constituted of high toughness materials which therefore break under Elastic–Plastic conditions. The method presented combines calculations based on Strength of Materials and Elastic–Plastic fracture mechanics to calculate the crack driving force applied to the cracked section by making its movements and those of the rest of the structure compatible. The application of the proposed assessment procedure is illustrated through the resolution of a practical example. The results of an experimental test, with the same geometrical configuration as the practical example, are also provided.     

Three‐dimensional mixed‐mode (I and II) crack‐front fields in ductile thin plates — effects of T‐stress

It has been well‐established that the non‐singular T‐stress provides a first‐order estimate of geometry and loading mode (e.g. tension versus bending) effects on elastic–plastic crack‐front field under mode I loading conditions. The objective of this paper is to exam the T‐stress effect on three‐dimensional (3D) crack‐front fields under mixed‐mode (modes I and II) loading. To this end, detailed 3D small strain, elastic–plastic simulations are carried out using a 3D boundary layer (small‐scale yielding) formulation. Characteristics of near crack‐front fields are investigated for a wide range of T‐stresses (T/σ₀ = ?0.8, ?0.4, 0.0, 0.4, 0.8). The plastic zones and thickness and angular and radial variations of the stresses are studied, corresponding to two values of the remote elastic mixity parameters M_e = 0.3 and 0.7, under both low and high levels of applied loads. It is found that different T‐stresses have a significant effect on the plastic zones size and shapes, regardless of the mode mixity and load level. The thickness, angular and radial distributions of stresses are also affected markedly by T‐stress. It is important to include these effects when investigating the mixed‐mode ductile fracture failure process in thin‐walled structural components.     

Quantification of constraint effects in fracture mechanism transition for cracked structures under mixed mode loading

The objective of the study is to evaluate the effects of plastic constraint on transition between tensile-type and shear-type fracture. The T -stress is employed as the quantifying parameter for constraint and is incorporated into the existing theoretical criteria for modelling this transition. It is found that different constraint levels can dramatically alter the transition point. To verify this finding, two sets of mixed mode tests with different constraint levels are carried out. Alongside the theoretical and experimental study, finite element simulation is performed to verify and support these findings. Substantially improved agreement is observed with experimental data if the effect of plastic constraint on transition is included.     

Fatigue crack propagation under in‐phase and out‐of‐phase biaxial loading

Fatigue damage characteristics of aluminium alloy under complex biaxial loads such as in‐phase and out‐of‐phase loading conditions and different biaxiality ratios have been investigated. The effects of microscale phenomena on macroscale crack growth were studied to develop an in‐depth understanding of crack nucleation and growth. Material characterization was conducted to study the microstructure variability. Scanning electron microscopy was used to identify the second phase particles, and energy dispersive X‐ray spectroscopy was performed to analyse their phases and elements. Extensive quasi‐static and fatigue tests were conducted on Al7075‐T651 cruciform specimens over a wide range of load ratios and phases. Detailed fractography analysis was conducted to understand the crack growth behaviour observed during the fatigue tests. Significant differences in crack initiation and propagation behaviour were observed when a phase difference was applied. Primarily, crack retardation and splitting were observed because of the constantly varying mode mixity caused by phase difference. The crack growth behaviour and fatigue lives under out‐of‐phase loading were compared with those under in‐phase loading to understand the effect of mixed‐mode fracture.     

Investigation of ductile rupture in U‐notched Al 6061‐T6 plates under mixed mode loading

The aim of the present research is to evaluate ductile failure of U‐notched components under mixed mode I/II loading conditions. For this purpose, first, several rectangular plates made of the aluminium alloy Al 6061‐T6 and weakened by central bean‐shaped slit with two U‐shaped ends are tested under mixed mode I/II loading conditions, and the load‐carrying capacity of the specimens are experimentally measured. Then, using the equivalent material concept, Al 6061‐T6, which is a highly ductile material, is equated with a virtual brittle material, and the load‐carrying capacity of the same U‐notched specimens virtually made of the equivalent material is theoretically predicted by using two well‐known stress‐based brittle fracture criteria. Finally, the theoretical failure loads of the virtual specimens are compared with the experimental ones of the real Al 6061‐T6 specimens. It is revealed that the experimental results could very well be predicted by means of both brittle fracture criteria without conducting time‐consuming elastic–plastic analyses.     

Formulation and numerical treatment of incompressibility constraints in large strain elastic–plastic analysis

The present paper is concerned with an efficient framework for a nonlinear finite element procedure for the rate‐independent finite strain analysis of solids undergoing large elastic‐isochoric plastic deformations. The formulation relies on the introduction of a mixed‐variant metric deformation tensor which will be multiplicatively decomposed into a plastic and an elastic part. This leads to the definition of an appropriate logarithmic strain measure which can be additively decomposed into the exact isochoric (deviatoric) and volumetric (spheric) strain measures. This fact may be seen as the basic idea in the formulation of appropriate mixed finite elements which guarantee the accurate computation of isochoric strains. The mixed‐variant logarithmic elastic strain tensor provides a basis for the definition of a local isotropic hyperelastic stress response whereas the plastic material behavior is assumed to be governed by a generalized J₂ yield criterion and rate‐independent isochoric plastic strain rates are computed using an associated flow rule. On the numerical side, the computation of the logarithmic strain tensors is based on higher‐order Padé approximations. To be able to take into account the plastic incompressibility constraint a modified mixed variational principle is considered which leads to a quasi‐displacement finite element procedure. Finally, the numerical solution of finite strain elastic‐plastic problems is presented to demonstrate the efficiency and the accuracy of the algorithm. Copyright © 1999 John Wiley & Sons, Ltd.     

Experimental fracture investigation of CNT/epoxy nanocomposite under mixed mode II/III loading conditions

For the first time, the brittle fracture of epoxy‐based nanocomposite reinforced with MWCNTs (multi‐walled carbon nanotubes) and subjected to mixed mode II/III loading conditions is investigated. This experimental investigation is carried out using a newly developed test configuration. Araldite LY 5052 epoxy, which is a resin frequently used in aerospace industry, is utilized to fabricate pure epoxy and nanocomposite test specimens with two different MWCNTs contents of 0.1 and 0.5 wt%. The obtained experimental results reveal that adding MWCNTs to epoxy resin up to 0.5 wt% improves the fracture toughness under pure mode II and pure mode III loading with an increasing trend. This is while the improvement under mixed mode II/III loading is reduced by adding nanotubes more than 0.1 wt%. To justify the variations of fracture toughness in terms of nanoparticles content, SEM (scanning electron microscopy) photographs of the fracture surfaces of the specimens in the vicinity of the initial crack front are prepared. Additional fracture mechanisms caused by adding carbon nanotubes are discussed in detail based on the provided SEM images.     

Digital image correlation measurement of near‐tip fatigue crack displacement fields: constant amplitude loading and load history effects

Recent work by de Matos and colleagues employed digital image correlation to measure near tip displacement fields for fatigue cracks in 6082 T6 aluminium alloy. The main focus of this work was to directly measure fatigue crack closure, but the measurements can also be used to examine conditions at and ahead of the crack tip. In this paper, the results are re‐analysed and compared to two crack‐tip deformation models. The first assumes simple elastic deformation (according the Westergaard solution). This allows the history of crack‐tip stress intensity to be examined. Reasonable agreement with the elastic model is obtained, although there is a residual stress intensity caused by the plastic wake, which gives rise to crack closure. The second model examined is a simple elastic–plastic assumption, proposed by Pommier and colleagues. This can be applied to constant amplitude loading, although the results obtained here are very similar to the elastic case. A slightly more complex load case (a single overload in an otherwise constant amplitude variation of load) gives a much more complicated crack‐tip history. Here, the importance of crack‐tip plastic displacement, represented by the second term in Pommier's model becomes much clearer. Load history effects are captured by the residual value of this term and its associated displacement fields as well as by stress intensity factor. The implications for further modelling and experimental work are discussed.     

In‐plane and out‐of‐plane constraint for single edge notched bending specimen and cruciform specimen under uniaxial and biaxial loading

Considering fracture constraint is an efficient way to describe stress–strain field and fracture toughness more accurately, so it is necessary to realise the relationship with in‐plane and out‐of‐plane constraint for different standard specimens. In this paper, three‐dimensional finite element method is applied to study the in‐plane and out‐of‐plane constraint for both cruciform specimen and single edge notched bending specimen made from commercial pure titanium. Crack length and in‐plane loading as the factors affecting in‐plane constraint, and thickness as the factor affecting the out‐of‐plane constraint are used to study the effect on both in‐plane and out‐of‐plane constraint in this paper. From the results, in‐plane and out‐of‐plane constraint are both related to specimen geometries and loading styles. And there exist relationships with in‐plane and out‐of‐plane constraint because of factors for different specimens. Depending on crack length, out‐of‐plane constraint increases with in‐plane constraint. While depending on transverse loading, out‐of‐plane constraint decreases with in‐plane constraint. In addition, when the in‐plane constraint of a specimen is higher, in‐plane constraint increases with out‐of‐plane constraint (thickness). When the in‐plane constraint is lower, in‐plane constraint almost remains unchanged with out‐of‐plane constraint.     

Advances in J‐integral estimation of circumferentially surface cracked pipes

This paper describes enhanced J‐integral estimation schemes for pipes with circumferential semi‐elliptical cracks subjected to tensile loading, global bending and internal pressure. These schemes are given in two different forms to cover the wide ranges of geometries and material parameters; the modified GE/EPRI method and the modified reference stress method. In the former method, new plastic influence functions for fully plastic J‐integral estimation are developed based on extensive three‐dimensional finite element calculations. In the latter method, new optimized reference loads are suggested and utilized to predict the J values. To verify the feasibility of these two schemes, J‐integral values obtained from further detailed FE analyses are compared to those from the proposed schemes. Because the estimated J‐integrals agree fairly well with the detailed FE analysis results, the new solutions can be applied for accurate structural integrity assessment of different size pipes with a circumferential surface crack.     

On improved crack tip plastic zone estimates based on T‐stress and on complete stress fields

Cracked ductile structures yield locally to form a plastic zone (pz) around their crack tips, which size and shape controls their structural behaviour. Classical pz estimates are based solely on stress intensity factors (SIF), but their precision is limited to very low σ_n/S_Y nominal stress to yield strength ratios. T‐stresses are frequently used to correct SIF‐based pz estimates, but both SIF and SIF plus T‐stress pz estimates are based on truncated linear elastic (LE) stress fields that do not satisfy boundary conditions. Using Griffith's plate complete LE stress field to avoid such truncated pz estimates, the influence of its Williams’ series terms on pz estimation is evaluated, showing that T‐stress improvements are limited to medium σ_n/S_Y values. Then, corrections are proposed to introduce equilibrium requirements into LE pz estimates. Finally, these improved estimates are compared with pz calculated numerically by an elastic–plastic finite element analysis.     

Evaluation of a phenomenological elastic‐plastic approach for magnesium alloys under multiaxial loading conditions

Magnesium alloys are greatly appreciated due to their high strength to weight ratio, stiffness, and low density; however, they can exhibit complex types of cyclic plasticity like twinning, de‐twinning, or Bauschinger effect. Recent studies indicate that these types of cyclic plastic deformations cannot be fully characterized using the typical tools used in cyclic characterization of steels and aluminium alloys; thus, it is required new approaches to fully capture their cyclic deformation and plasticity. This study aims to propose and evaluate a phenomenological cyclic elastic‐plastic approach designed to capture the cyclic deformation of magnesium alloys under multiaxial loading conditions. Series of experimental tests were performed to characterize the cyclic mechanical behaviour of the magnesium alloy AZ31BF considering proportional loadings with different strain amplitude ratios and a nonproportional loading with a 45° phase shift. The experimental results were modulated using polynomial functions in order to implement a cyclic plasticity model for the AZ311BF based on the phenomenological approach proposed. Results show good correlations between experiments and estimates.     

Mode I–III Decomposition of the J‐integral from DIC Displacement Data

This paper presents the implementation of the decomposition method on digital image correlation (DIC) obtained displacement fields to obtain J‐integral results (J) and respective stress intensity factors (SIFs). DIC is increasingly used with the J‐integral approach in experimental mechanics to obtain J estimates from complex fracture processes. In this approach, the decomposition method is applied to DIC displacement fields for the first time. Here, displacement fields are separated before stresses and strains are computed, so that subsequent computation of separate J or SIF components may follow the classical full‐field J‐integral approach. The sensitivity of the decomposition method to experimental errors is investigated using synthetically generated errors imposed on crack tip displacement fields (Williams' series), from which improvements to the procedure are proposed. The method is experimentally tested on PMMA Arcan specimens under mode I, II, and III, and mixed‐mode I–III loading. Test results were compared to fracture toughness values obtained from ASTM tests and literature with close agreement.