**Interacting Cracks Analysis Using Finite Element Method**

Ruslizam Daud, Ahmad Kamal Ariffin, Shahrum Abdullah and Al Emran Ismail

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/54358

## **1. Introduction**

358 Applied Fracture Mechanics

39.

Zhang S, Shibata T, Haruna T. (1997) Initiation and propagation of IGSCC for s ensitized type 304 stainless steel in dilute sulfate solutions. *Corrosion Science*; 39(9):1725–

> This chapter aims to introduce the concept of fracture mechanics and numerical approaches to solve interacting cracks problems in solid bodies which involves elastic crack interaction. The elastic crack interaction is a result of changes in stress field distribution as the applied force is given during remote loading. The main emphasis is to address the computational evaluation on mechanistic models based on crack tip displacement, stress fields and energy flows for multiple cracks. This chapter start with a brief discussion on fracture and failure that promoted by interacting cracks from industrial cases to bring the issues of how important the crack interaction behaviour is. The present fracture and failure mechanism is assumed to exhibit the brittle fracture. Thus, the concept of linear elastic fracture mechanics (LEFM) is discussed regarding the crack interaction model formulation. As the elastic crack interaction is concerned, the previous analytical and numerical solution of crack interaction are elaborated comprehensively corresponds to fitness-for-service (FFS) as published by ASME boiler and pressure vessel code (Section XI, Articles IWA-3330), JSME fitness-for-service code and BSI PD6493 and BS7910. A new computational fracture mechanics algorithm is developed by adopting stress singularity approach in finite element (FE) formulation. The result of developed approach is discussed based on the crack interaction limit (CIL) aspects and crack unification limit (CUL) in pertinent to the equality of two cracks to single crack rules in FFS. As a conclusion, the FE formulated approach was found to be at agreeable accuracy with analytical formulation and FFS at certain range of crack interval.

## **2. Fracture and failures by interacting cracks**

This section provides the overview of failure cases in industries and related fitness-forsevice (FFS) codes which used to assess any cracks or flaws that detected in structures. The works on solution models for FFS codes improvement in specific cases of interacting cracks also discussed.

© 2012 Daud et al., licensee InTech. This is an open access chapter distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2012 Daud et al., licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **2.1. Industrial failures**

In this section, the fracture and failure by interacting cracks is explained by examples from industrial failures. Mechanical structures and components are designed with multiple stress concentration features (SCF) such as notches, holes, corners and bends. For example, welding and riveting in joining and fastening process have consumed to the increase of stress concentration factor. In every SCF, there is a critical point that experienced the highest concentrated stress field, named as multiple stress riser points (MRSP). Under multiple mode of loading and environmental effect, the interaction between SCF and MRSP tends to form multiple cracks in various types of cracks (e.g. straight crack, surface crack and curved crack) before the cracks propagate in various path to coalescence, overlapping, overlapping, branching and finally fracture in brittle manner.

Crack interaction that induced from MRSP has caused many catastrophic failures , for example in aircraft fuselage (Hu, Liu, & Barter, 2009), F-18 Hornet bulkhead (Andersson et al., 2009), rotor fault (Sekhtar, 2008), gigantic storage tanks (Chang & Lin, 2006), oil tankers (Garwood, 2001), polypropylene tank (Lewis & Weidmann, 2001) and the most recent is the fail of helicopter longerons (J. A. Newman, Baughman, & Wallace, 2010). The recent lab experimental work on multiple crack initiation, propagation and coalescence by (Park & Bobet, 2010) and metallurgical work by (J. A. Newman et al., 2010) supported the important role of crack interaction in fracture and failure. The above cases proved how crucial and important the research on crack interaction is.

To explain further, failure in aircraft is considered as an example. Al Alloy has been extensively used for the fabrication of fuselage, wing, empennage, supporting structures that involve many fastening and joining points. Under static, cyclic loading and environmental effect, the micro-cracks are initiated from MRSP. To certain extent, out of many factors, brittle failure may happen through catastrophic failure (Hu, Liu, & Barter, 2009). As the distance between MRSP is close, the interaction between cracks is become more critical. The fracture behavior due to interacting cracks as the distance between cracks is closed need more understanding. The conventional fracture mechanics may be insufficient to support.

In this case, the advancement of computational fracture mechanic may contribute a lot in crack interaction research and increase the accuracy of failure prediction (Andersson et al., 2009). Most recent structural failure being reported is dealing with fuselage joints in large aircraft structures (Hu et al., 2009). The aircraft fuselage structure is made by 7050-T7451 aluminum alloy and designed to have multiple shallow notches that purposely used for lap joints. Due to environmental reaction and variable magnitude of operational loading, multisite cracks are formed at MRSP. Therefore, the interaction of multi-site cracks is needed to be quantified for possible coalescence between cracks. The challenge is how to predict the accurate coalescence and fatigue crack growth for the multiple crack problems. Table 1 provides the list of industrial failures that originated from various kinds of multiple cracks. Under loading condition such as mechanical or thermal loading, it is observed that most of the cracks or flaws formation is start on the surface of the body rather than embedded inside the body.


**2.1. Industrial failures** 

branching and finally fracture in brittle manner.

important the research on crack interaction is.

insufficient to support.

the body.

In this section, the fracture and failure by interacting cracks is explained by examples from industrial failures. Mechanical structures and components are designed with multiple stress concentration features (SCF) such as notches, holes, corners and bends. For example, welding and riveting in joining and fastening process have consumed to the increase of stress concentration factor. In every SCF, there is a critical point that experienced the highest concentrated stress field, named as multiple stress riser points (MRSP). Under multiple mode of loading and environmental effect, the interaction between SCF and MRSP tends to form multiple cracks in various types of cracks (e.g. straight crack, surface crack and curved crack) before the cracks propagate in various path to coalescence, overlapping, overlapping,

Crack interaction that induced from MRSP has caused many catastrophic failures , for example in aircraft fuselage (Hu, Liu, & Barter, 2009), F-18 Hornet bulkhead (Andersson et al., 2009), rotor fault (Sekhtar, 2008), gigantic storage tanks (Chang & Lin, 2006), oil tankers (Garwood, 2001), polypropylene tank (Lewis & Weidmann, 2001) and the most recent is the fail of helicopter longerons (J. A. Newman, Baughman, & Wallace, 2010). The recent lab experimental work on multiple crack initiation, propagation and coalescence by (Park & Bobet, 2010) and metallurgical work by (J. A. Newman et al., 2010) supported the important role of crack interaction in fracture and failure. The above cases proved how crucial and

To explain further, failure in aircraft is considered as an example. Al Alloy has been extensively used for the fabrication of fuselage, wing, empennage, supporting structures that involve many fastening and joining points. Under static, cyclic loading and environmental effect, the micro-cracks are initiated from MRSP. To certain extent, out of many factors, brittle failure may happen through catastrophic failure (Hu, Liu, & Barter, 2009). As the distance between MRSP is close, the interaction between cracks is become more critical. The fracture behavior due to interacting cracks as the distance between cracks is closed need more understanding. The conventional fracture mechanics may be

In this case, the advancement of computational fracture mechanic may contribute a lot in crack interaction research and increase the accuracy of failure prediction (Andersson et al., 2009). Most recent structural failure being reported is dealing with fuselage joints in large aircraft structures (Hu et al., 2009). The aircraft fuselage structure is made by 7050-T7451 aluminum alloy and designed to have multiple shallow notches that purposely used for lap joints. Due to environmental reaction and variable magnitude of operational loading, multisite cracks are formed at MRSP. Therefore, the interaction of multi-site cracks is needed to be quantified for possible coalescence between cracks. The challenge is how to predict the accurate coalescence and fatigue crack growth for the multiple crack problems. Table 1 provides the list of industrial failures that originated from various kinds of multiple cracks. Under loading condition such as mechanical or thermal loading, it is observed that most of the cracks or flaws formation is start on the surface of the body rather than embedded inside

**Table 1.** Summary of industrial failures caused by multiple crack interaction

Crack interaction intensity exist in the form of stress shielding and amplification. The failure mechanism by cracks interaction may occurs under brittle, plastic and creep failure. In general, the interacting cracks problems that promotes the fracture and failure of structures are solved using the advancement of fracture mechanics. The fracture mechanics solution can be accomplished through analytical, numerical and experimental work. Based on the above individual approach or combination of them, typically the solution is represented by a model. The model may be defined based on the uncertainty input for the model such as crack geometry, loading and material properties. For crack interaction problems that randomness is relatively small, the deterministic analysis is the best to considered rather than probabilistic analysis. The model is suitable for any deterministic system response. However, when the randomness is relatively high, the model system response required a more robust solution, known as probabilistic approach.

Multiple cracks interaction can be defined into elastic crack interaction and plastic crack interaction that may be referred to theory of elasticity and plasticity, respectively. Under loading condition, the high stresses near crack tips usually accompanied by inelastic deformation and other non-linear effect. If the inelastic deformation and other non-linear effect are relatively small compared to crack sizes and other geometrical body characteristic, the linear theory is most adequate to address the crack interaction behavioral problems. Thus, the role of elastic driving force originated from crack tips can be translated into elastic crack interaction. Then, for every type of interaction, it may classified into interaction without crack propagation (EIWO) where the interaction occurs in the region of SIF *K Kc* and interaction with crack propagation (EIWI) occurs in the region of SIF *K Kc* . In this case of EIWO, the quantification of crack extension is neglected. The EIWO becomes the main issue of interaction in present study since it is inadequate investigation on crack interaction limit and multiple to single crack equivalencies. The study on EIWO is typically measured the fracture parameter and its behavior based on SIF while EIWI more focused to evaluation and prediction of crack path, propagation, coalescence, branching and crack arrestment. Under mechanical or thermal loading condition, the generated interaction will varies depends on type of loading mode (e.g. Mode I and Mode II) being applied. The preceding sections outlined the related fitness-for-service (FFS) codes in pertinent to failures that caused by crack interaction.

## **2.2. Fitness-for-Service (FFS) codes**

This section presents the guidelinea that have been published in fitness-for-service (FFS) codes. The related investigation works also discussed as the FFS codes are evaluated in different case of interacting cracks problems. ASME Boiler and Pressure Vessel Code Section XI (ASME, 1998, 2004) and API 579 (ASME, 2007) defines that the multiple cracks are assumed to be independent until or unless the following conditions are satisfied. In the case of two parallel cracks in solid bodies, if distance between crack planes *d d* 12.7 mm , the cracks are treated as being coplanar. For coplanar cracks, if the distance between cracks *s* 1,2 2 (maximum of { } ) *i i s a* . Indeed, the single enveloping crack is assumed equal to coplanar cracks if the condition of crack depth *a* and surface length *c* satisfy maximum of { } , 1,2 1 2 / 2 *i i a ac ccs* . If the non-coplanar cracks in overlapped condition at *s* 0 , the cracks are assumed to be coplanar, the surface length *c new* 1 2 *c ccs* . Thus, coalescence occur at max(2 ,2 ) 1 1 *s aa* . In United Kingdom, the engineering critical assessment (ECA) of potential or actual defects in engineering structures is codified into two prime standard; British Standard PD6495 (BSI, 1991) and Nuclear Electric CEGB R6 (R6, 2006). BSI PD6495 has replaced by BS7910 (BSI, 1997), and most latest (BSI, 2005). Original BSI PD6495 primarily concerned on assessment of defect welds. The PD 6495 used crack tip open displacement (CTOD) and SIF *K* based analysis while the BS7910 used CTOD, *K* or equivalent *K* that derived from *J-*integral. Both codes define the cracks are assumed to be independent until or unless the following distance between crack planes *d* satisfy 1 2 *d aa* 0.5 and the cracks are treated as being coplanar. For coplanar cracks, if the distance between cracks *s* 1,2 2 (minimum of { } ) *i i s c* , a single enveloping crack is assumed and the coalescence occur at 1 2 *s cc* 0.5(2 2 ) . The R6 is an approach to upgrade the BSI PD6495 by Central Electricity Generating Board (CEGB) that focused on operating equipment at high temperature where the assumption of equal or greater fracture and possibility of plastic collapse together and fracture separately. The concept of failure assessment diagram (FAD) is introduced to occupy the need of fracture parameter of plasticity fracture. R6 provides a special form of *J*-integral analysis to impose the plastic collapse limit. Details of FAD can be found in (R6, 2006). JSME Fitness-for-Service Code provides no prescription for the interference between multiple cracks or flaws. In (JSME, 2000), in example of parallel offset cracks,the multiple cracks are replaced by an equivalent single crack based on the stage of detected cracks with satisfying the condition of *S H* 5mm 10mm and *S HS* 5mm 2 where *S* is relative vertical spacing and *H* is relative horizontal spacing. When the crack tips distance *S* 0 due to overlapped condition, the crack growth evaluation is considered about the coalescence stages. The guideline in JSME Code is based on experimental results. In (JSME, 2008), the interacting cracks are combined in crack growth prediction and the judgment is based on the relative spacing *S* and *H* at the initial condition. If the relative spacing at the beginning of the growth prediction meets the criterion, two cracks are combined when the distance *S* become zero during the crack growth.

#### **2.3. Interacting cracks models for Fitness-for-Service**

362 Applied Fracture Mechanics

Crack interaction intensity exist in the form of stress shielding and amplification. The failure mechanism by cracks interaction may occurs under brittle, plastic and creep failure. In general, the interacting cracks problems that promotes the fracture and failure of structures are solved using the advancement of fracture mechanics. The fracture mechanics solution can be accomplished through analytical, numerical and experimental work. Based on the above individual approach or combination of them, typically the solution is represented by a model. The model may be defined based on the uncertainty input for the model such as crack geometry, loading and material properties. For crack interaction problems that randomness is relatively small, the deterministic analysis is the best to considered rather than probabilistic analysis. The model is suitable for any deterministic system response. However, when the randomness is relatively high, the model system response required a

Multiple cracks interaction can be defined into elastic crack interaction and plastic crack interaction that may be referred to theory of elasticity and plasticity, respectively. Under loading condition, the high stresses near crack tips usually accompanied by inelastic deformation and other non-linear effect. If the inelastic deformation and other non-linear effect are relatively small compared to crack sizes and other geometrical body characteristic, the linear theory is most adequate to address the crack interaction behavioral problems. Thus, the role of elastic driving force originated from crack tips can be translated into elastic crack interaction. Then, for every type of interaction, it may classified into interaction without crack propagation (EIWO) where the interaction occurs in the region of SIF *K Kc* and interaction with crack propagation (EIWI) occurs in the region of SIF *K Kc* . In this case of EIWO, the quantification of crack extension is neglected. The EIWO becomes the main issue of interaction in present study since it is inadequate investigation on crack interaction limit and multiple to single crack equivalencies. The study on EIWO is typically measured the fracture parameter and its behavior based on SIF while EIWI more focused to evaluation and prediction of crack path, propagation, coalescence, branching and crack arrestment. Under mechanical or thermal loading condition, the generated interaction will varies depends on type of loading mode (e.g. Mode I and Mode II) being applied. The preceding sections outlined the related fitness-for-service (FFS) codes in pertinent to failures

This section presents the guidelinea that have been published in fitness-for-service (FFS) codes. The related investigation works also discussed as the FFS codes are evaluated in different case of interacting cracks problems. ASME Boiler and Pressure Vessel Code Section XI (ASME, 1998, 2004) and API 579 (ASME, 2007) defines that the multiple cracks are assumed to be independent until or unless the following conditions are satisfied. In the case of two parallel cracks in solid bodies, if distance between crack planes *d d* 12.7 mm , the cracks are treated as being coplanar. For coplanar cracks, if the distance between cracks *s* 1,2 2 (maximum of { } ) *i i s a* . Indeed, the single enveloping crack is assumed equal to

more robust solution, known as probabilistic approach.

that caused by crack interaction.

**2.2. Fitness-for-Service (FFS) codes** 

The solution for interacting surface or embedded cracks that based on FFS revision is limited in literature especially for interacting parallel edge cracks in finite body. Therefore, a review on available developed technique or models that related to FFS is presented in this section. The need of continues revision on FFS limitations based on ASME Boiler and Pressure Vessel Code Section XI and British Standard BS7910 started by the industrial failures in all major pressure vessels (Burdekin, 1982). The pressure vessels are designed and built to comply with ASME Boiler and Pressure Vessel Code Section XI and British Standard BS7910 codes but the failure occurrences are significantly high. Both codes can be expressed as

$$\begin{array}{ll}\text{ASME} & s \leq 2 \times \text{(maximum of } \{a\_1, a\_2\}\text{)}\\\text{BSIPD6493} & s \leq 2 \times \text{(maximum of } \{c\_1, c\_2\}\text{)} \end{array} \tag{1}$$

To investigate the problem, (Burdekin, 1982) studied the interaction between the collapse and fracture in pressure vessels using the approach that successfully applied to nuclear applications. The approach applied the LEFM and EPFM using COD and *J*-integral based on single crack under bending condition. The study revealed the important of fracture mechanics as a tool for interaction in failures. This work can be considered as among the first work that put concern on the FFS codes.

Similarly, (O'donoghue et al., 1984) investigated the formation of elliptical cracks in aircraft and pressure vessel attachment lugs and identified the formation is due to stress risers and cracks interaction. Two equal coplanar surface cracks under Mode I loading are modeled using the proposed finite element alternating method (FEAM) in finite solid and the FE analysis results are compared to ASME Boiler and Pressure Vessel Code procedure. The interaction effects are defined by proposed magnification factor (normalized SIF) and the magnification factors seem to increase due to the increase interaction of two cracks and the depth of cracks. The Section XI of ASME Boiler and Pressure Vessel Code recommend that two interacting surface flaws in a pressure vessel should be modeled by a single elliptical crack that covers both flaws. It can be seen that SIF for single crack as proposed by FFS code are generally larger than those due to two interacting cracks as proposed by (O'donoghue et al., 1984). This trend of magnification factor shows that the ASME Boiler and Pressure Vessel Code in Section XI procedure will tend to underestimate the design life of multiple flaws structures. The ASME pressure vessel codes (ASME, 1998) and British Standard PD6495 (BSI, 1991) do not quantify the interaction between cracks especially in two close proximity cracks. At sufficiently close distance, the interaction may cause the increase of SIF. The exclusion of crack interaction may result with unrealistic SIF. Therefore, with the concern on the above standard guideline, (Leek & Howard, 1994) presents an empirical method to approximate the interaction factor of two coplanar surface cracks under tension and bending loading. The approximation approach resulted with good agreement with FE analysis using developed BERSAFE program , (Murakami & Nasser, 1982) and (J. C. Newman & Raju, 1981) within ± 5% discrepancy.

$$a = \text{maximum of } \{a\_i\}\_{i=1,2} \text{ } \ c = c\_1 + c\_2 \tag{2}$$

$$(s/\overline{c}) \times (s/\overline{a}) > 3.38 \quad \text{and} \quad s/\overline{a} \ge 2.49 \tag{3}$$

Based on the (ASME, 2004) and (BSI, 1991) design code that expressed by Eq. (1), (Moussa et al., 1999) used FEM to analyze interaction of two identical parallel non-coplanar surface cracks subjected to remote tension and pure bending loads. The interaction factor as a function of stress shielding to cause overlapping in distance is studied using three dimensional linear finite element analyses. The formation of stress relaxation state is introduced near crack front, as a form of shielding effect at sufficient overlapping. *J*-integral is calculated based on models by (Shivakumar & Raju, 1992) and the interaction factor is defined as follows:

$$\gamma = \mathcal{K}\_{in} \;/\; \mathcal{K}\_{is} \tag{4}$$

where *Kin* and *Kis* are the SIFs with and without the influence of interaction, respectively. As conclusion, the interaction effect appears to diminish as the value of *s*/*c* approaches 2.0. The existing rules for re-characterization of interacting cracks as less conservative for high values of *s*/*c* and over-conservative as *s*/*c* is close to 2.0.

364 Applied Fracture Mechanics

first work that put concern on the FFS codes.

Newman & Raju, 1981) within ± 5% discrepancy.

defined as follows:

maximum of { } , *i i* 1,2 1 2 *a ac c c* (2)

( / ) ( / ) 3.38 and / 2.49 *sc sa s a* (3)

is

Based on the (ASME, 2004) and (BSI, 1991) design code that expressed by Eq. (1), (Moussa et al., 1999) used FEM to analyze interaction of two identical parallel non-coplanar surface cracks subjected to remote tension and pure bending loads. The interaction factor as a function of stress shielding to cause overlapping in distance is studied using three dimensional linear finite element analyses. The formation of stress relaxation state is introduced near crack front, as a form of shielding effect at sufficient overlapping. *J*-integral is calculated based on models by (Shivakumar & Raju, 1992) and the interaction factor

(1)

ASME 2 (maximum of { , }) BSIPD6493 2 (maximum of { , })

To investigate the problem, (Burdekin, 1982) studied the interaction between the collapse and fracture in pressure vessels using the approach that successfully applied to nuclear applications. The approach applied the LEFM and EPFM using COD and *J*-integral based on single crack under bending condition. The study revealed the important of fracture mechanics as a tool for interaction in failures. This work can be considered as among the

Similarly, (O'donoghue et al., 1984) investigated the formation of elliptical cracks in aircraft and pressure vessel attachment lugs and identified the formation is due to stress risers and cracks interaction. Two equal coplanar surface cracks under Mode I loading are modeled using the proposed finite element alternating method (FEAM) in finite solid and the FE analysis results are compared to ASME Boiler and Pressure Vessel Code procedure. The interaction effects are defined by proposed magnification factor (normalized SIF) and the magnification factors seem to increase due to the increase interaction of two cracks and the depth of cracks. The Section XI of ASME Boiler and Pressure Vessel Code recommend that two interacting surface flaws in a pressure vessel should be modeled by a single elliptical crack that covers both flaws. It can be seen that SIF for single crack as proposed by FFS code are generally larger than those due to two interacting cracks as proposed by (O'donoghue et al., 1984). This trend of magnification factor shows that the ASME Boiler and Pressure Vessel Code in Section XI procedure will tend to underestimate the design life of multiple flaws structures. The ASME pressure vessel codes (ASME, 1998) and British Standard PD6495 (BSI, 1991) do not quantify the interaction between cracks especially in two close proximity cracks. At sufficiently close distance, the interaction may cause the increase of SIF. The exclusion of crack interaction may result with unrealistic SIF. Therefore, with the concern on the above standard guideline, (Leek & Howard, 1994) presents an empirical method to approximate the interaction factor of two coplanar surface cracks under tension and bending loading. The approximation approach resulted with good agreement with FE analysis using developed BERSAFE program , (Murakami & Nasser, 1982) and (J. C.

*s a a s c c*

> The existing (ASME, 1998; JSME, 2000) FFS combination rules provided no prescription for the interference between multiple cracks for corrosion fatigue. Therefore, (Kamaya, 2003) developed simulation model to extent the condition of coalescence rules in (JSME, 2000) for crack growth process using body force method (BFM). BFM is used to investigate the multiple cracks growth in stress corrosion cracking. Based on JSME code and the SIF value of coalescence behavior from experiments, the new SIF formulation is developed using BFM where focus is given to the interaction between cracks under various relative position and size. The crack propagation direction can be written as

$$\theta\_{\text{max}} = \mp \cos^{-1} \left( \left( 3K\_{II}^2 + K\_I \sqrt{8K\_{II}^2 + K\_I^2} \right) / \left( 9K\_{II}^2 + K\_I^2 \right) \right) \tag{5}$$

where the sign in Eq. (5) positive in the case of / 0 *K K II I* and negative in the case of / 0 *K K II I* . When the crack are close and overlapped, the crack interaction intensity between cracks is almost equivalent to single coalesced crack. The change of inner crack tips direction also found with little influence on the crack growth behavior. The relative crack length and position influenced the crack interaction intensity.

The combination rule in ASME Code is found to provide the relative large overestimation of the actual crack growth since the complex growth phenomena under interaction are summarized in simple combination rules. In order to reduce the conservativeness in existing code, (Kamaya, 2008b) proposed alternative assessment procedures based on the size of area and fatigue crack growth. Experimental analysis and testing is conducted using stainless steel specimens (A-H0S5 and B-H0S5) subjected to cyclic tensile loading. FE analysis is carried out to simulate the crack growth during coalescence. In the simulation, the automatic meshing was generated by command language in PATRAN and the SIF is derived from energy release rate obtained from virtual crack extension integral method using ABAQUS. The normalized SIF of Mode I, *KI* is expresses as

$$F\_I = K\_I / \sigma\_0 \sqrt{\pi a} \tag{6}$$

where *<sup>o</sup>* denotes the applied tensile stress and the *a* is the maximum depth. As a result, the area of the crack face is concluded to be the predominant parameter for the crack growth of interacting cracks under test condition. The cracks of various shapes can be characterized as semi-elliptical cracks of the same area. In extension of parallel semi-elliptical cracks study, (Kamaya, 2008a) investigated the coalescence of adjacent cracks as a result of crack growth with the influence of crack interaction. The magnitude of interaction is represented by driving force of the crack growth (CGF), written as

$$\mathcal{W}\_m = \sum\_{i=1}^{n-1} 0.5 \left( D\_p K\_{I(i)}^{m\_p} + D\_p K\_{I(i+1)}^{m\_p} \right) \Delta \mathbf{g} \tag{7}$$

where *Dp* and *mp* are the material constant, *KI i*( ) denotes the Mode I SIF of the *i* th node from 0 *<sup>o</sup> p* . *g* and *n* are the distance of neighboring mode on the crack front and number of nodes. The CGF formulation proved that the interaction between surface cracks not only dependent on relative spacing but also the position of crack front. In the condition of *S* > 0, as the cracks overlapped, the stress shielding effect influenced the change of CGF. The most important, the study notified the cracks can be replaced with single crack of the same area when the relative spacing is sufficiently close, at crack spacing *H* < *a*. In regular inspection of pressure vessel components, the adjacent defects are found close enough. Under operational loading, the stress field around the crack tips will be magnified and accelerates the crack growth rate. This matter has been referred to current fitness-for-service (FFS) rules such as ASME Boiler and Pressure Vessel Code Section XI (ASME, 1998, 2004), API 579 (ASME, 2007), British Standard PD6495 (BSI, 1991) , BS7910 (BSI, 1997, 2005) and Nuclear Electric CEGB R6 (R6, 2006). The multiple interacting cracks are combined as single crack as the two cracks satisfy the prescribed criterion. As observed, this rule introduced unrealistic discontinuity in the process of crack growth due to the crack interaction is neglected. The evaluation of two interacting coplanar cracks in plates under tension is conducted by (Xuan, Si, & Tu, 2009) and creep interaction factor *creep* is introduced by using *C*\*integral prediction analysis and the FE analysis is executed using ABAQUS to verify the proposed approach. Creep interaction factor *creep* is expressed as

$$\mathcal{N}\_{crep} = \left(\boldsymbol{\mathcal{C}}\_{Double}^{\*} / \boldsymbol{\mathcal{C}}\_{Single}^{\*}\right)^{1/2} \tag{8}$$

In conclusion, the creep crack interaction represented by *C*\*integral is affected by crack configuration (e.g. relative crack distance *c/d*, depth of crack *a/t* and location at crack front 2*ø/*π) and time dependent properties of material such as creep exponent *n*. The increasing crack aspect *a/c* resulted with no significant effect to *C*\* integral.

Most recent, (Kamaya et al., 2010) used S-version finite element method to determine the SIF changes due to the interaction of stress field which caused variation in crack growth rate and cracks shape. The root of interaction problems is referred to (JSME, 2000, 2008) and (ASME, 2004, 2007) for the case of interacting dissimilar crack sizes. However, the effect of difference crack size or relative size effect is not taken into account in the aforementioned code. The results have shown that smaller cracks stopped growing when the difference in size of interaction was large enough. It means, the interaction effect on the fatigue life of the larger cracks was negligibly small. Moreover, the offset distance and the relative size were important parameter for interaction evaluation especially when the *S* 0 and the condition of crack spacing 1 *H c* / and cracks ratio 2 1 *c c* / must be considered most. In present study, the focus is given to determine the stage of crack interaction intensity is equal to single crack in a state of crack interaction limit (CIL) and unification (CUL) using finite element method.

## **3. Finite element analysis**

366 Applied Fracture Mechanics

driving force of the crack growth (CGF), written as

Si, & Tu, 2009) and creep interaction factor *creep*

Creep interaction factor *creep*

with the influence of crack interaction. The magnitude of interaction is represented by

1

*i*

 <sup>1</sup> ( ) ( 1)

(7)

is introduced by using *C*\*integral prediction

*creep Double Single C C* (8)

0.5 *p p <sup>n</sup> m m m p I i p I i*

*W DK DK g*

where *Dp* and *mp* are the material constant, *KI i*( ) denotes the Mode I SIF of the *i* th node from 0 *<sup>o</sup> p* . *g* and *n* are the distance of neighboring mode on the crack front and number of nodes. The CGF formulation proved that the interaction between surface cracks not only dependent on relative spacing but also the position of crack front. In the condition of *S* > 0, as the cracks overlapped, the stress shielding effect influenced the change of CGF. The most important, the study notified the cracks can be replaced with single crack of the same area when the relative spacing is sufficiently close, at crack spacing *H* < *a*. In regular inspection of pressure vessel components, the adjacent defects are found close enough. Under operational loading, the stress field around the crack tips will be magnified and accelerates the crack growth rate. This matter has been referred to current fitness-for-service (FFS) rules such as ASME Boiler and Pressure Vessel Code Section XI (ASME, 1998, 2004), API 579 (ASME, 2007), British Standard PD6495 (BSI, 1991) , BS7910 (BSI, 1997, 2005) and Nuclear Electric CEGB R6 (R6, 2006). The multiple interacting cracks are combined as single crack as the two cracks satisfy the prescribed criterion. As observed, this rule introduced unrealistic discontinuity in the process of crack growth due to the crack interaction is neglected. The evaluation of two interacting coplanar cracks in plates under tension is conducted by (Xuan,

is expressed as

crack aspect *a/c* resulted with no significant effect to *C*\* integral.

analysis and the FE analysis is executed using ABAQUS to verify the proposed approach.

 1/2 \* \* /

In conclusion, the creep crack interaction represented by *C*\*integral is affected by crack configuration (e.g. relative crack distance *c/d*, depth of crack *a/t* and location at crack front 2*ø/*π) and time dependent properties of material such as creep exponent *n*. The increasing

Most recent, (Kamaya et al., 2010) used S-version finite element method to determine the SIF changes due to the interaction of stress field which caused variation in crack growth rate and cracks shape. The root of interaction problems is referred to (JSME, 2000, 2008) and (ASME, 2004, 2007) for the case of interacting dissimilar crack sizes. However, the effect of difference crack size or relative size effect is not taken into account in the aforementioned code. The results have shown that smaller cracks stopped growing when the difference in size of interaction was large enough. It means, the interaction effect on the fatigue life of the larger cracks was negligibly small. Moreover, the offset distance and the relative size were

#### **3.1. Finite element analysis**

The stress in the neighborhood of a crack tip in homogenous isotropic material exist in a form of square-root singular and there have been many special elements or singularity function based approach were described in details in (Banks-Sills, 1991). The square root singular stresses in the neighborhood can be modeled by quarter-point, square and collapsed, triangular elements for two dimensional problems, and by brick and collapsed, prismatic elements in three dimensions. Quarter-points square have been found to produce the most excellent results (Banks-Sills, 2010). The stiffness matrix of the element is evaluated using two-dimensional integral based on Gaussian quadrature approach. The plate is constructed with a consideration of singular element and assigned to both crack tips *Ct1* and *Ct2*. It is because the high gradients of singular stress-strain and deformation fields are concentrated at both crack tips. The SIF calculation is limited to linear elastic problem with a homogeneous, isotropic material near the crack region.

### *3.1.1. Singularity stress field*

The studies are conducted in a pure Mode I loading condition with specified material, Alloy 7475 T7351 solid plate in constant thickness, homogenous isotropic continuum material, linear elastic behavior, small strain and displacements, and crack surface are smooth. According to Westergaard method for single crack, Mode I *KI* and Mode II *KII* SIF can be expressed as:

$$\begin{aligned} \mathbf{K}\_{I} &= \mathbf{F} \;/\, \mathbf{K}\_{o} = \mathbf{F} \{ \mathbf{b} \;/\, a, a \;/\, \mathbf{W} \} \; \sigma\_{o} \sqrt{\pi a} \\ \mathbf{K}\_{II} &= \mathbf{F} \;/\, \mathbf{K}\_{o} = \mathbf{F} \{ \mathbf{b} \;/\, a, a \;/\, \mathbf{W} \} \; \tau\_{o} \sqrt{\pi a} \end{aligned} \tag{9}$$

where *<sup>o</sup>* is nominal stress, *<sup>o</sup>* is shear stress, *W* is width of specimen, *a* and *b* is the length of crack and crack interval, respectively. The work starts by determination of *KI* and *KII* using Eq. (9). The important issue which differs from single crack is the existence of cracks interaction in fracture analysis.

Consider two multiple edge crack of length 1*a* and 2 *a* which occupies the segment of 0.05 ≤ *a/W* ≤ 0.5 and 0.5 ≤ *b/a* ≤ 3.0 in finite plate subjected to uniform equal stress *σ* along the *y* direction, as shown in Fig. 1(a) and (b). The SIF formulation is based on the creation of singular element at the crack tip based on quadratic isoparametric finite element developed in ANSYS evironment based on (Madenci & Guven, 2006), where the element is based on Barsoum (Barsoum, 1974, 1975), as depicted in Fig. 2 (a) and (b). The singularity is obtained by shifting the mid-side node the ¼ point close to the crack tip. To calculate the SIF, the elements are assumed to be in rigid body motion and constant strain modes. The master element mapping in Cartesian space is transformed into curvilinear space using Jacobian transformation which used to interpolate the displacement within the elements (Chandrupatla & Belegundu, 2002). The accuracy of special element has been addressed by (Murakami, 1976) where the crack tip nodal point is enclosed by a number of special element. In analysis, the size, number and compatibility of special elements really affect the accuracy. The special elements also defined as singularity function methods where stress singularity at crack tip is modeled. The condition of continuity between elements is the most important. By using singularity function method, Mode I and Mode II of stress intensity factor may be able to calculate with high accuracy (Shields, Srivatsan, & Padovan, 1992).

**Figure 1.** Barsoum singular element for (a) strong crack interaction and (b) weak crack interaction

According to Eq. (9), the shape correction factor can be converted to new elastic interaction factor

$$\mathcal{N}\_{I,in,D} = \mathcal{K}\_I \mid \mathcal{K}\_o \tag{10}$$

$$\mathcal{Y}\_{\text{II},in,\text{D}} = \mathcal{K}\_{\text{II}} / \mathcal{K}\_o \tag{11}$$

where *<sup>I</sup>* and *II* denotes to elastic interaction factor for Mode I and II fracture, respectively.

The SIF for Mode I and Mode II is determined using Displacement Extrapolation Method (DEM) by written APDL macro code in ANSYS (Madenci & Guven, 2006), and expressed as

$$K\_{l} = \frac{E}{3(1+\nu)(1+\kappa)}\sqrt{\frac{2\pi}{L}}\left(4\left(v\_{2}-v\_{4}\right) - \frac{\left(v\_{3}-v\_{5}\right)}{2}\right) \tag{12}$$

$$K\_{ll} = \frac{E}{3(1+\nu)(1+\kappa)}\sqrt{\frac{2\pi}{L}}\left(4\left(\mu\_2 - \mu\_4\right) - \frac{\left(\mu\_3 - \mu\_5\right)}{2}\right) \tag{13}$$

where, *E=*Young Modulus, 3 4 for plain stress, 3 4 /1 for plain strain, *L* is length of element, *v* and *u* are displacements in a local Cartesian coordinate system and υ is Poisson's ratio.

Source: (Barsoum, 1974, 1976; Henshell & Shaw, 1975)

**Figure 2.** (a) Eight nodes quadratic isoparametric elements (b) Parent element

#### **4. Findings and discussion**

Two parallel edge cracks interaction will mainly referred to shielding effect rather than amplification effect. The crack interaction is proportional to the magnitude of elastic interaction factor *<sup>I</sup>* . The crack interaction will only exist at *b/a* < 3 (Z.D.Jiang, A.Zeghloul, G.Bezine, & J.Petit, 1990; Z.D.Jiang, J.Petit, & G.Bezine, 1990), the analytical formulation can be expressed as

$$K\_I = K\_o F\_n(a / W, b / a) \tag{14}$$

in which

368 Applied Fracture Mechanics

factor

by shifting the mid-side node the ¼ point close to the crack tip. To calculate the SIF, the elements are assumed to be in rigid body motion and constant strain modes. The master element mapping in Cartesian space is transformed into curvilinear space using Jacobian transformation which used to interpolate the displacement within the elements (Chandrupatla & Belegundu, 2002). The accuracy of special element has been addressed by (Murakami, 1976) where the crack tip nodal point is enclosed by a number of special element. In analysis, the size, number and compatibility of special elements really affect the accuracy. The special elements also defined as singularity function methods where stress singularity at crack tip is modeled. The condition of continuity between elements is the most important. By using singularity function method, Mode I and Mode II of stress intensity factor may be able to calculate with high accuracy (Shields, Srivatsan, & Padovan, 1992).

**Figure 1.** Barsoum singular element for (a) strong crack interaction and (b) weak crack interaction

According to Eq. (9), the shape correction factor can be converted to new elastic interaction

, , /

, , /

*I in D I o K K* (10)

*II in D II o K K* (11)

$$\begin{aligned} \mathcal{Y}\_{I,in,li} &= F\_u = F\_2 = \frac{A\_0 + A\_1 \left(a \, \, \dot{w}\right)^{1.5} + A\_2 \left(a \, \, \dot{w}\right)^4}{\sqrt{1 - \left(a \, \, \dot{w}\right)^2}}\\ A\_0 &= 0.79 + 0.07 \left(b \, \, \dot{a}\right) + 0.04 \left(b \, \, \dot{a}\right)^2 - 0.011 \left(b \, \, \dot{a}\right)^3\\ A\_1 &= 1.74 + 2.84 \left(b \, \, \dot{a}\right) - 1.44 \left(b \, \, \dot{a}\right)^2 + 0.206 \left(b \, \, \dot{a}\right)^3\\ A\_2 &= 6.02 + 2.19 \left(b \, \, \dot{a}\right) - 3.26 \left(b \, \, \dot{a}\right)^2 + 0.828 \left(b \, \, \dot{a}\right)^3 \end{aligned} \tag{15}$$

and analytical single edge crack SIF reference by (Brown & Strawley, 1966),

$$\begin{aligned} K\_{I,ref} &= \sigma\_0 \sqrt{\pi a} \left( 1.12 - 0.23 \left( a/\,\,\text{W} \right) + 10.6 \left( a/\,\,\text{W} \right)^2 - 21.7 \left( a/\,\,\text{W} \right)^3 + 30.4 \left( a/\,\,\text{W} \right)^4 \right) \\ K\_{I,ref} &= \sigma\_0 \sqrt{\pi a} \left( f\_{I,ref,BS} \right) \end{aligned} \tag{16}$$

#### **4.1. Crack interaction factor** *I in D* , ,  **comparison with analytical data** *I in Ji* , ,

The mode I fracture of the elastic interaction factors ,, ,,,1 ,,,2 (,) *I in D I in D ct I in D ct* for crack interval ratio *b a* / 1.5 3.0 and *a W*/ 0.05 0.5 are shown in Fig. 3-6. Overall, it can be seen that the interaction factor varies with the different *a W*/ values where the crack interaction factor increases as the *a W*/ increases and vice versa. The point of intersection also observed occurred for all the *b a* / values. For example, the plot of ,, ,,,1 ,,,2 (,) *I in D I in D ct I in D ct* against *a W*/ at *b a* / 3.0 are illustrated in Fig. 3.

The *I in D* , , prediction line is compared with the predicted result of single edge crack *I ref BS* , , *f* . A general good agreement can be observed with a minimum difference 0.6% at *a W*/ 0.1 and maximum difference 5.39 % difference *a W*/ 0.15 . In comparison, the present *I in D* , , has demonstrated more accurate prediction compared with *I in Ji* , , results. For example, in reference to *I ref BS* , , *f* , at *a W*/ 0.5 , the present prediction *I in D* , , is in difference of *I in D ct* ,,,1 0.85 % and *I in D ct* ,,,2 0.4%, while *I in Ji* , , is at 2.15% difference. In terms of the CIL point, the closer the crack interaction to *I ref BS* , , *f* of single crack, a more accurate CIL prediction can be achieved. It is noted that the *I in Ji* , , analytical expression was formulated using the numerical results of *J*-integral analysis. The formulation is unable to calculate the crack interaction factor for both crack tips and become the weakness of *I in Ji* , , . Therefore, the present work of DEM has improved the existing *J*-integral analysis by improving the accuracy of CIL to predict fracture due to crack interaction.

Theoretically, the study of intersection point is most significant in identification of crack interaction limit (CIL) and crack unification limit (CUL). The intersection point of two cracks *I in D* , , with single crack *I ref BS* , , *f* justifies the realization of CIL at higher *a W*/ and CUL at lower *a W*/ . The CIL and CUL also differ with different *b a* / . From Fig. 3, the identified CUL is at *a W*/ 0.1 and CIL approximately at *a W*/ 0.5 .

**Figure 3.** Variation of *I in D* , , against *a W*/ for *b a* / 3.0

370 Applied Fracture Mechanics

0

, ,

*I in Ji n*

*F F*

1

2

*I ref BS*

0.85 % and *I in D ct* ,,,2

improving the accuracy of CIL to predict fracture due to crack interaction.

CIL prediction can be achieved. It is noted that the *I in Ji* , ,

, ,

*a f*

0

,, ,,,1 ,,,2 (,) *I in D I in D ct I in D ct*

 

The *I in D* , , 

present *I in D* , , 

*I in D* , , 

difference of *I in D ct* ,,,1 

 ,

*K*

*K*

*I ref*

,

*I ref*

0

**4.1. Crack interaction factor** *I in D* , ,

/ /

1.5 4

2 3

(15)

(16)

for crack

results.

.

is in

2 3

2 3

 **comparison with analytical data** *I in Ji* , ,

against *a W*/ at *b a* / 3.0 are illustrated in Fig. 3.

prediction line is compared with the predicted result of single edge crack

234

is at 2.15% difference. In terms

analytical expression was

 

*a w*

 

1.12 0.23 / 10.6 / 21.7 / 30.4 /

*a aW aW aW aW*

1 /

*A Aaw A aw*

0 1 2 <sup>2</sup> <sup>2</sup>

0.79 0.07 / 0.04 / 0.011 / 1.74 2.84 / 1.44 / 0.206 / 6.02 2.19 / 3.26 / 0.828 /

*A ba ba ba A ba ba ba A ba ba ba*

 

The mode I fracture of the elastic interaction factors ,, ,,,1 ,,,2 (,) *I in D I in D ct I in D ct*

interval ratio *b a* / 1.5 3.0 and *a W*/ 0.05 0.5 are shown in Fig. 3-6. Overall, it can be seen that the interaction factor varies with the different *a W*/ values where the crack interaction factor increases as the *a W*/ increases and vice versa. The point of intersection also observed occurred for all the *b a* / values. For example, the plot of

*I ref BS* , , *f* . A general good agreement can be observed with a minimum difference 0.6% at *a W*/ 0.1 and maximum difference 5.39 % difference *a W*/ 0.15 . In comparison, the

For example, in reference to *I ref BS* , , *f* , at *a W*/ 0.5 , the present prediction *I in D* , ,

has demonstrated more accurate prediction compared with *I in Ji* , ,

0.4%, while *I in Ji* , ,

of the CIL point, the closer the crack interaction to *I ref BS* , , *f* of single crack, a more accurate

formulated using the numerical results of *J*-integral analysis. The formulation is unable to calculate the crack interaction factor for both crack tips and become the weakness of *I in Ji* , ,

Therefore, the present work of DEM has improved the existing *J*-integral analysis by

Theoretically, the study of intersection point is most significant in identification of crack interaction limit (CIL) and crack unification limit (CUL). The intersection point of two cracks

with single crack *I ref BS* , , *f* justifies the realization of CIL at higher *a W*/ and CUL at

and analytical single edge crack SIF reference by (Brown & Strawley, 1966),

Another significant improvement of present *I in D* , , is the moving intersection point as *a W*/ decreases for every *b a* / , as shown in Fig. 3-6. From Fig. 4, the moving intersection point can be noticed moves from *a W*/ 0.1 to *a W*/ 0.075 . The moving intersection point exhibits similar prediction trend of *I in D* , , . The intersection point also can be denoted as the crack unification limit (CUL) point, which indicates the starting point of strong interaction region start approximately at *a W*/ 0.07 .

The intersection point is also observed to move from *a W*/ 0.075 for *b a* / 2.5 , *a W*/ 0.06 for *b a* / 2.5 , *a W*/ 0.05 for *b a* / 2.0 and *a W*/ 0.05 for *b a* / 1.5 , as shown in Fig. 4-6. It means that the CUL is not in a fix limit, it exist in dynamic condition which depends on crack interval ratio *b a* / . Conversely, the *I in Ji* , , prediction model overruled the FFS codes because it does not lead to a single independent or combined crack because of not having any intersection point. The intersection point could not be defined by *I in Ji* , , prediction model.

**Figure 4.** Variation of *I in D* , , against *a W*/ for *b a* / 2.5

**Figure 5.** Variation of *I in D* , , against *a W*/ for *b a* / 2.0

It also can be seen that *I in Ji* , , prediction is unable to display the unification of crack interaction factor, which defined equivalent to *I ref BS* , , *f* of single crack. In this case, the *I in Ji* , , prediction is expected to encounter some numerical errors since at lower *a W*/ 0.07 . In analysis, the ratio of crack length and width also define the critical stress field, if the path independent radius of *J*-integral for both crack tips are overlapped, the calculation of *J* value might be overestimated and the stress behavior is equal to behavior of single edge crack in finite body. The path integral line should be always apart and controlled in individual condition.

Based on the FFS codes, the multiple cracks are assumed to be independent as single cracks or combined cracks, until or unless certain conditions are satisfied. The intersection point that lies in the present *I in D* , , prediction trend curve, which intersects with the single crack prediction, shows good agreement with the outlined FFS codes. It is seen that at range of 0.05 / 0.15 *a W* , the *I in D* , , value is about to level at 1.072 -1.085. The small changes in these range indicate that shielding effect is very small and promote the unification in interaction at the point of *a W*/ 0.15 . The smaller the *b a* / , the faster unification process starts.

**Figure 6.** Variation of *I in D* , , against *a W*/ for *b a* / 1.5

#### **4.2. Mode II fracture behavior**

372 Applied Fracture Mechanics

**Figure 4.** Variation of *I in D* , ,

0.8

1

1.2

1.4

1.6

1.8

*Crack interaction factor, γI*

2

2.2

2.4

2.6

2.8

3

**Figure 5.** Variation of *I in D* , ,

0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3

*Crack interaction factor, γI*

*b/a* = 2.0

*Intersection point*

> *I in Ji* , ,

against *a W*/ for *b a* / 2.5

Brown & Strawley (1966)

Present (ct1) Present (ct2) Jiang et al. (1990)

Brown & Strawley (1966)

Present (ct1) Present (ct2) Jiang et al. (1990)

*b/a* = 2.5

*Intersectio n point*

> *I in Ji* , ,

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

 *Intersection point*

*Intersection point*

*I ref BS* , , *f*

, , , ,,1 , ,,2 (, ) *I in D I in D ct I in D ct*

 

*I ref BS* , , *f*

, , , ,,1 , ,,2 (, ) *I in D I in D ct I in D ct*

 

*Crack to width ratio , a/W*

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

*Crack to width ratio , a/W*

against *a W*/ for *b a* / 2.0

Fig. 7 depicts the trend of *KI in D* , , and *KII in D* , , against *b a* / for *a W*/ 0.25 0.5 . It can be seen that the increase in *b a* / results in an increase of *KI in D* , , and a decrease of *KII in D* , ,

**Figure 7.** Variation of *KI in D* , , and *KII in D* , , for ( / 0.25 0.5 ) *a W*

values. In strong interaction region 0 / 1.0 *b a* , the values of *KI in D* , , increased rapidly for higher crack-to-width ratio *a W*/ 0.5,0.45 and before grew slowly as the value of *a W*/ decreases to *a W*/ 0.25 . Meanwhile, the values of *KII in D* , , declined significantly at higher crack-to-width ratio *a W*/ 0.5,0.45 and before decreased slightly as the value of *a W*/ decreases to *a W*/ 0.25 . For both crack interaction phase, the value of *KI in D* , , is always much higher than *KII in D* , , . In weak interaction region 1.5 / 3.0 *b a* , it can be seen that the values of *KI in D* , , increased slightly before growing slowly and then maintaining at the same level to steady state at *b a* / 3.0 . At the same region, the value of *KII in D* , , declined moderately before decreased slightly and remain stable at the level of *b a* / 3.0 . It means that mode II SIF is less influenced by damage shielding effect than mode I SIF. It also defined that the crack opening is more affected by damage shielding effect than the crack sliding. This has been clearly indicated in Fig. 7 (a)-(f) and Fig. 8(a)-(d).

374 Applied Fracture Mechanics

**Figure 7.** Variation of *KI in D* , , and *KII in D* , , for ( / 0.25 0.5 ) *a W*

**Figure 8.** Variation of *KI in D* , , and *KII in D* , , for ( / 0.05 0.2) *a W*

From these figures, it also observed that the different between mode II SIF and mode I SIF is reduced as the *a W*/ decreases. In the context of crack interaction limit (CIL) based on *KI in D* , , and *KII in D* , , , by considering the convergence level as the indicator of CIL. It can be seen that the degree and speed of CIL achievement not only depends heavily on the increased of *b a* / , the reduction of *a W*/ from *a W*/ 0.5 0.05 also provides significant impact on CIL determination.

Fig. 8 shows the variation of *KI in D* , , and *KII in D* , , for ( / 0.05 0.2) *a W* . It can be observed that at smallest *a W*/ 0.05 , the value of *KI in D* , , and *KII in D* , , is almost hold at constant value and stabilize.

The identification of CIL in this condition is absent because the value of *KI in D* , , and *KII in D* , , are about the same value for all *b a* / . This equalization condition may be referred to the identification of crack unification limit (CIL). Overall, it can be concluded that the higher ratio of *b a* / and *a W*/ , the more the realization of CIL. Inversely, the lower ratio of *b a* / and *a W*/ , the more indication to CUL can be realized.

## **5. Conclusion**

The numerical solution based on displacement extrapolation method (DEM) has proved to be more consistent in SIF prediction comparedd to for both crack tips. However, the DEM able to predict the SIF for Mode I and Mode II fracture behaviour. The FE results conclude that the interaction of two cracks is directly influence the reduction of SIF magnitude and *<sup>I</sup>* 

at the crack tips. The parallel cracks have experienced decrease shielding effect as the cracks interval *b/a* decrease. The identification of crack interaction limit (CIL) and crack unification limit (CUL) has been accomplished. The SIF of Mode I is found more significant compared to Mode II in higher or lower *b a* / and *a W*/ ratio. Mode II SIF can be neglected due to its small effect to the stress shielding effect. The FFS codes rules that define the combination of two cracks as single crack are well translated as CIL and CUL.

## **Author details**

Ruslizam Daud, Ahmad Kamal Ariffin, Shahrum Abdullah and Al Emran Ismail *Universiti Kebangsaan Malaysia, Malaysia* 

## **6. References**

Andersson, B., Blom, A. F., Falk, U., Wang.G.S, Koski, K., Siljander, A., et al. (2009). *A case study of multiple site fatigue crack growth in the F-18 Hornet bulkhead*. Paper presented at the 25th ICAF Symposium.

ASME. (1998). ASME Boiler and Pressure Vessel Code, Section XI. *New York, USA*.


Barsoum, R. S. (1976). Application of Triangular Quarter-Point Elements as Crack Tip Elements of Power Law Hardening Material. *International Journal of Fracture, 12*, 463-466.

Brown, W. F., & Strawley, J. E. (1966). ASTM STP 410.

376 Applied Fracture Mechanics

impact on CIL determination.

and *a W*/ , the more indication to CUL can be realized.

two cracks as single crack are well translated as CIL and CUL.

Ruslizam Daud, Ahmad Kamal Ariffin, Shahrum Abdullah and Al Emran Ismail

ASME. (1998). ASME Boiler and Pressure Vessel Code, Section XI. *New York, USA*.

value and stabilize.

**5. Conclusion** 

**Author details** 

**6. References** 

*Universiti Kebangsaan Malaysia, Malaysia* 

the 25th ICAF Symposium.

From these figures, it also observed that the different between mode II SIF and mode I SIF is reduced as the *a W*/ decreases. In the context of crack interaction limit (CIL) based on *KI in D* , , and *KII in D* , , , by considering the convergence level as the indicator of CIL. It can be seen that the degree and speed of CIL achievement not only depends heavily on the increased of *b a* / , the reduction of *a W*/ from *a W*/ 0.5 0.05 also provides significant

Fig. 8 shows the variation of *KI in D* , , and *KII in D* , , for ( / 0.05 0.2) *a W* . It can be observed

that at smallest *a W*/ 0.05 , the value of *KI in D* , , and *KII in D* , , is almost hold at constant

The identification of CIL in this condition is absent because the value of *KI in D* , , and *KII in D* , , are about the same value for all *b a* / . This equalization condition may be referred to the identification of crack unification limit (CIL). Overall, it can be concluded that the higher ratio of *b a* / and *a W*/ , the more the realization of CIL. Inversely, the lower ratio of *b a* /

The numerical solution based on displacement extrapolation method (DEM) has proved to be more consistent in SIF prediction comparedd to for both crack tips. However, the DEM able to predict the SIF for Mode I and Mode II fracture behaviour. The FE results conclude that the interaction of two cracks is directly influence the reduction of SIF magnitude and *<sup>I</sup>*

at the crack tips. The parallel cracks have experienced decrease shielding effect as the cracks interval *b/a* decrease. The identification of crack interaction limit (CIL) and crack unification limit (CUL) has been accomplished. The SIF of Mode I is found more significant compared to Mode II in higher or lower *b a* / and *a W*/ ratio. Mode II SIF can be neglected due to its small effect to the stress shielding effect. The FFS codes rules that define the combination of

Andersson, B., Blom, A. F., Falk, U., Wang.G.S, Koski, K., Siljander, A., et al. (2009). *A case study of multiple site fatigue crack growth in the F-18 Hornet bulkhead*. Paper presented at


JSME. (2008). *JSME Fitness-for-Service Code S NA1-2008*.


Moussa, W. A., Bell, R., & Tan, C. L. (1999). The interaction of two parallel non-coplanar identical surface cracks under tension and bending. *International Journal of Pressure Vessels and Piping, 76*, 135-145.

378 Applied Fracture Mechanics

1349.

307-321.

*Mechanics, 51*(2), 233-238.

*International Journal, 46*(1), 15-23.

*of Pressure Vessel and Piping, 47*, 247-257.

*Vessels and Piping, 68*, 181-201.

*Fracture Mechanics, 9*, 17-24.

*ANSYS* (1 ed.): Springer Science-Business.

with short cracks. *J. Biomechanics, 23*(10), 967-975.

factor. *Engineering Fracture Mechanics, 40*(3), 585-592.

JSME. (2000). *JSME Fitness-for-Service Code S NA1-2000*. JSME. (2008). *JSME Fitness-for-Service Code S NA1-2008*.

Jeong, D. Y., & Brewer, J. C. (1995). On the linkup of multiple cracks. *Engineering Fracture* 

Jones, R., Peng, D., & Pitt, S. (2002). Assessment of multiple flat elliptical cracks with

Kamaya, M. (2003). A crack growth evaluation method for interacting multiple cracks. *JSME* 

Kamaya, M. (2008a). Growth evaluation of multiple interacting surface cracks. Part II: Growth evaluation of parallel cracks. *Engineering Fracture Mechanics, 75*, 1350-1366. Kamaya, M. (2008b). Growth evaluation of multiple interacting surface cracks.Part I: Experiments and simulation of coalesced crack. *Engineering Fracture Mechanics, 75*, 1336-

Kamaya, M., Miyokawa, E., & Kikuchi, M. (2010). Growth prediction of two interacting surface cracks of dissimilar sizes. *Engineering Fracture Mechanics, 77*, 3120-3131. Kirkhope, K. J., Bell, R., & Kirkhope, J. (1991). Stress intensity factors for single and multiple semi-elliptical surface cracks in pressurized thick-walled cylinders. *International Journal* 

Kobayashi, H., & Kashima, K. (2000). Overview of JSME flaw evaluation code for nuclear

Lakes, R. S., Nakamura, S., Behiri, J. C., & Bonfield, W. (1990). Fracture mechanics of bone

Lam, K. Y., & Phua, S. P. (1991). Multiple crack interaction and its effect on stress intensity

Leek, T. H., & Howard, I. C. (1994). Estimating the elastic interaction factor of two coplanar surface cracks under Mode I load. *International Journal of Pressure Vessels and Piping, 60*,

Leek, T. H., & Howard, I. C. (1996). An examination of methods of assessing interacting surface cracks by comparison with experimental data. *International Journal of Pressure* 

Lewis, P. R., & Weidmann, G. W. (2001). Catastropic failure of a polypropylene tank Part I:

Madenci, E., & Guven, I. (2006). *The finite element method and application in engineering using* 

Matake, T., & Imai, Y. (1977). Pop-in behaviour induced by interaction of cracks. *Engineering* 

Meizoso, A. M., Esnaola, J. M. M., & Perez, M. F. (1995). Interaction effect of multiple crack

Milwater, H. R. (2010). A simple and accurate method for computing stress intensity factors

Mischinski, S., & Mural, A. (2011). Finite Element Modeling of Microcrack Growth in

growth on fatigue. *Theoretical and Applied Fracture Mechanics, 23*, 219-233.

primary investigation. *Engineering Failure Analysis, 6*(4), 197-214.

of collinear interacting cracks. *Aerospace Science and Technology*.

Cortical Bone. *Journal of Applied Mechanics, 78*, 1-9.

power plants. *International Journal of Pressure Vessels and Piping, 77*, 937-944.

interactions. *Theoretical and Applied Fracture Mechanics, 38*, 281-291ssss.


*Engineering Failure Analysis, 6*(4), 215-232.

*Fracture Mechanics, 37*(5), 1139-1144.

under tension. *Engineering Fracture Mechanics, 76*, 2192-2201.

a finite width sheet. *Engineering Fracture Mechanics, 35*(6), 1073-1079.

Weidmann, G. W., & Lewis, P. R. (2001). Catastropic failure of a polypropylene tank Part II: comparison of the DVS 2205 code of practice and the design of the failed tank.

Xuan, F.-Z., Si, J., & Tu, S. T. (2009). Evaluation of C\* integral for interacting cracks in plates

Yang, Y. H. (2009). Multiple parallel symmetric mode III cracks in a functionally graded material plane. *Journal of Solid Mechanics and Materials Engineering, 3*(5), 819-830. Z.D.Jiang, A.Zeghloul, G.Bezine, & J.Petit. (1990). Stress intensity factor of parallel cracks in

Z.D.Jiang, J.Petit, & G.Bezine. (1990). Fatigue propagation of two parallel cracks. *Engineering* 

## *Edited by Alexander Belov*

The book "Applied Fracture Mechanics" presents a collection of articles on application of fracture mechanics methods to materials science, medicine, and engineering. In thirteen chapters, a wide range of topics is discussed, including strength of biological tissues, safety of nuclear reactor components, fatigue effects in pipelines, environmental effects on fracture among others. In addition, the book presents mathematical and computational methods underlying the fracture mechanics applications, and also developments in statistical modeling of fatigue. The work presented in this book will be useful, effective, and beneficial to mechanical engineers, civil engineers, and material scientists from industry, research, and education.

Photo by VladNikon / iStock

Applied Fracture Mechanics

Applied Fracture Mechanics

*Edited by Alexander Belov*