# Application of a Cohesive Zone Model for Simulating Fatigue Crack Growth from Moderate to High [DELTA]K Levels of Inconel 718.

1. IntroductionMost of the service failures of aerospace components such as landing gears, turbine blades, fastener joints, and so on are caused by fatigue associated with cyclic loading. The fatigue damage of structural metals normally is associated with the initiation of macrocracks (the nucleation and growth of microcracks), such that a few dominated macrocracks grow until the final unstable fracture happens. By taking into account material degradation instead of one single major crack, continuum damage mechanics provides a sophisticated framework for the fatigue crack initiation analysis of metallic materials [1], while the fracture mechanics-based models are dominated approaches in simulating fatigue crack growth [2].

One of the extensively used models for predicting fatigue crack growth is the Paris' law [3] which relates the range of the stress intensity factor [DELTA]K and the crack growth rate da/dN through a power law in double logarithmic coordinates. Generally, the relationship between da/dN and [DELTA]K for most of metallic materials can be divided into three stages [4], as illustrated in Figure 1. Regime I corresponds to the formation of a fatigue crack at the low [DELTA]K, where crack growth cannot be observed if [DELTA]K is less than the threshold [DELTA][K.sub.th]. Regime II usually maintains the Paris' law, which exhibits the moderate fatigue crack growth rate (i.e., [10.sup.-6] mm/cycle < da/dN < [10.sup.-3] mm/cycle). The high crack growth rate (i.e., da/dN > [10.sup.-3] mm/cycle) achieves at Regime III in which [DELTA]K asymptotically approaches the fracture toughness [DELTA][K.sub.c]. Fatigue life predictions in engineering design have much success with the modifications of the Paris' law to determine crack growth rates [5, 6].

Since elastic or the small-scale yielding (SSY) condition at crack-tip is assumed in the K-based models, many attractive models have also been proposed to simulate fatigue crack growth with considerable plastic deformations [6-10]. One representative model proposed by Dowling and Begley [7] is constructed by replacing [DELTA]K with [DELTA]J-integral to extend the Paris' law into the large-scale yielding (LSY) conditions. However, some researchers [11, 12] have pointed out that J-integral is invalid upon unloading. Although other researchers [10] found that [DELTA]J-integral still maintains its path independence under cyclic loadings, nevertheless, it is a considerable issue how to use [DELTA]J-integral correctly in fatigue crack growth simulations.

The cohesive zone model (CZM) is another reliable approach which can describe the gradual degradation of stress-carrying ability at a propagating crack-tip within the framework of the nonlinear fracture mechanics. CZMs were firstly proposed by Dugdale [13] and Barenblatt [14] to consider the presence of process zone at crack-tip, where the local degradation during crack growth is determined by a given constitutive law (cohesive law) to relate the traction with the magnitude of material separation. Except the applications in brittle [15] and composite materials [16], cohesive models have been widely used in simulating the ductile fracture of metallic materials. Tvergaard and Hutchinson [17] introduced a trapeziform cohesive law to investigate crack growth resistance during the fracture process of an elastic-plastic solid. Chen et al. [18] simulated the ductile crack growth of the thick CT specimens made of the pressure vessel steel using a polynomial cohesive law. Uthaisangsuk et al. [19] applied a CZM between martensitic islands and ferritic grains to consider interface debonding in the dual-phase microstructure of DP600 steel. Li et al. [20] studied the ductile fracture of a notched Ti-6Al-4V specimen using a CZM. Several authors [21] summarized classifications and applications of the CZM, which demonstrate the importance of cohesive models in simulating ductile fracture.

For fatigue crack calculations, the material degradation under cyclic loading must be considered; namely, the cyclic cohesive zone model (CCZM) should be applied to simulate fatigue crack growth. So far, there are many existing CCZMs that are able to describe the linear Paris' region (Regime II) behavior. The first CCZM disinterring the Paris-like behavior of fatigue crack growth for metallic materials was proposed by Nguyen et al. [22]. They simulated fatigue crack growth of a central cracked aluminum specimen and demonstrated that the CCZM is capable of reproducing the Paris' law. The unloading-reloading hysteresis was applied to account for dissipative mechanisms, which prevents the "shakedown" effect efficiently in the crack growth process. However, they did not mention the predictive ability of the model in Regime III. Maiti and Geubelle [23] introduced a bilinear cohesive law with an irreversible unloading-loading path to assess capacity of the model for predicting fatigue crack growth in polymers. Although they concluded that the developed model can be used in both Regime II and Regime III for an epoxy, they did not show the applications of the proposed model in simulating ductile materials. Roe and Siegmund [12] considered the effects of the interface fatigue crack growth under different loading modes using an irreversible CZM, and while most of the simulation results were performed under the linear Paris' region, high rates associated with Regime III were not demonstrated. Ural et al. [24] used a damage-based cohesive model to predict the fatigue crack growth behavior in Regime II of an aluminum alloy. However, the application of this model has not been identified in Regime III also. Li and Yuan [25] used a normal stress-dominated cohesive model to simulate fatigue crack growth of cracked and notched specimen and reproduced the linear Paris-like behavior at different loading ratios; nonetheless, the predictive ability in Regime III has not been considered in the CZM. Jimenez and Duddu [26] proposed a strain energy release rate-dependent CZM to model high cycle fatigue delamination of the composite materials under both mode I and mixed mode loadings, and they pointed out that the CZM is suitable for the constant amplitude loading in Regime II and also expected to extended the model to describe all three regimes in fatigue crack simulations. CCZMs can also be applied in analyzing the elastoplastic fatigue cracks by introducing some stress triaxiality-dependent parameters [27, 28]; that is, Jha and Banerjee [28] proposed a triaxiality-dependent cohesive model to simulate the fatigue failure of ductile metals and proved that the CZM is suitable for the geometries with different stress states. Similarly, these models [27, 28] do not mention their applications for fatigue crack growth at the high [DELTA]K levels in Regime III.

Although the Paris-like behavior which corresponds to the moderate [DELTA]K levels in Regime II can be predicted by using the aforementioned CCZMs [11, 22-27], this is not the advantage of CCZMs, since the Paris regime associated with SSY has already been described successfully according to the K-based models. In order to resolve the more challenging task that to simulate the high fatigue crack growth rates in Regime III of the metallic materials, a CCZM was proposed in our previous work [29] in which two damage variables were defined to represent monotonic damage and fatigue damage, respectively. In fact, it is not easy to distinguish these two different damage mechanisms in fatigue crack growth before the final rupture occurs; thus, it is better to describe the material damage evolution uniformly. Meanwhile, it is anticipated that CCZMs have the potential to predict high fatigue crack growth rates in Regime III under the LSY conditions. For the purpose of easing the difficulty level for cohesive modeling of fatigue crack growth at LSY, the predictive ability of CCZMs in simulating high fatigue crack growth rates at SSY should be verified systematically.

Therefore, the objective of this work is to examine a CCZM for predicting high fatigue crack growth rates corresponding to Regime III under the SSY conditions. Meanwhile, the present study attempts to describe the material damage evolution uniformly and identify the effects of the cohesive model parameters on fatigue crack growth under both moderate and high [DELTA]K levels. The contents of this paper include the following aspects: in Section 2, the detailed formulations of the CCZM are introduced. In Section 3, the identification of the model parameters is performed firstly, then the numerical results are presented according to the compact-tension (CT) specimens and the single-edge notched (SEN) specimens made of the IN718 superalloy and validated against the experimental data from literatures. Finally, the major conclusions are summarized in Section 4.

2. Cohesive Zone Model

The cohesive model has the ability to consider the ductile rupture of the metallic materials under various loading conditions. For fatigue crack growth simulations, the CCZM has to be able to describe accumulative damage associated with cyclic loading. Assume that the degradation of the original cohesive envelope occurs with damage accumulation, and the constitutive law under fatigue loading is defined as

[mathematical expression not reproducible], (1)

where the normal traction [T.sub.n] is assumed to be a piecewise function depending upon the normal separation [[delta].sub.n]. [[delta].sub.0] is the characteristic length of the cohesive law, which is denoted as the cohesive length. [[delta].u] denotes the maximum normal separation above which the normal traction vanishes. Following the concept from the continuum damage models, the damage variable D = [[integral].sub.t] [??] dt stands for the material degradation and represents failure if D approaches 1. [[sigma].max] = [[sigma].max],0] (1-D) is the cohesive strength used to identify the peak value of the normal traction and [[sigma].sun.max,0] is the initial cohesive strength. The original cohesive law under monotonic loading can be obtained by substituting D = 0 into (1). The parameter [[delta].sub.c] corresponds to the material separation relating to the damage variable D and [[delta].sub.0] [less than or equal to] [[delta].sub.c] [less than or equal to] [[delta].sub.u]. For both constant and variable amplitude loadings, the damage evolution equation suggested by Roe and Siegmund [12] can be applied to reflect the accumulative damage for the metallic materials, which reads

[mathematical expression not reproducible], (2)

where [d.sub.[summation]] is the accumulated cohesive length used to scale the increment of the normal separation [mathematical expression not reproducible] is the cohesive endurance limit. The Macaulay brackets <> are used to signify nonnegative value of damage evolution, and the Heaviside function H states that damage initiates at [bar.[delta]] > [[delta].sub.0].

The novel cohesive envelope determined by (1) will degrade and change its shape correspondingly associated with the damage accumulation, and this characteristic is the major difference between the current model and other CCZMs [11, 12, 20, 22-28]. For example, in [25], the combination use of (2) and the cohesive envelope of the exponential form are performed to simulate fatigue crack growth; however, this combination results in the contradictory damage evolution at [[delta].sub.n] > [[delta].sub.0] if a slightly larger value of [[sigma].sub.f]/[[sigma].sub.max,0] is provided, and this problem can be solved by using (1). As mentioned above, the similar CCZM has been proposed to investigate the mixed mode fatigue crack growth in our previous simulations [29], the difference between these two models is that the current one does not introduce the monotonic damage, the material degradation under cyclic loadings is described uniformly by (2), and the monotonic damage is determined by the cohesive envelope implicitly. Actually, the combination of (1) and (2) has the capacity to simulate fatigue crack growth at the high [DELTA]K levels, the applications need not be distinguished between monotonic damage and fatigue damage as introduced in [29], the key point in such simulations is the form of (1), and more details will be performed in Section 3.

For the numerical implementation of the model under each unloading-reloading cycle, assume that the residual separation exists after a complete unloading cycle as plastic deformations, and the corresponding unloading-reloading path can be denoted by

[T.sub.n] = [T.sub.n,t] + [k.sub.n] ([[delta].sub.n] - [[delta].sub.n,t]), [[delta].sub.n], [[delta].sub.n,t] [greater than or equal to] 0, (3)

with

[k.sub.n] = [[sigma].sub.max] exp (1)/[[delta].sub.0], (4)

where [T.sub.n] is the traction at time t + [DELTA]t related to the separation [[delta].sub.n] according to the cohesive stiffness [k.sub.n]. [T.sub.n,t] and [[delta].sub.n,t] are the traction and the separation at time t, respectively.

As shown in Figure 2(a), assume that the rising segment (1 [right arrow] 2) of the loading sequence results in the separation of a material point at the cohesive zone. In order to explain the operating mechanism of the model under cyclic loading in detail, two different cases are considered, namely, [[delta].sub.n] [less than or equal to] [[delta].sub.0] (Figure 2(b)) and [[delta].sub.n] > [[delta].sub.0] (Figure 2(c)). Upon unloading, the path (2 [right arrow] 3) follows (3) without damage evolution. The contact of two cohesive surfaces is constrained by introducing a penalized equation as

[T.sub.n] = A [[sigma].sub.max,0] exp (1) [[delta].sub.n]/[[delta].sub.0], [[delta].sub.n] < 0, (5)

where A is the penalty stiffness. Here, we should mention that this contact algorithm is applied to simulate crack closure explicitly. The crack closure effects are related to the plastic deformation of the bulk material around the crack-tip implicitly under cyclic loading.

The reloading path (3 [right arrow] a) specified by (3) deviates from the previous unloading path due to damage evolution. At the subsequent loading segment (3 [right arrow] 4), the value of [T.sub.n] corresponding to reloading may reach the degraded cohesive envelope [T.sup.1.sub.n] at a. Upon further reloading, the value of [T.sub.n] moves along the degraded envelope [T.sup.1.sub.n] at time t; however, the continued damage accumulation through (2) makes it move along the subsequent degraded cohesive envelope [T.sup.2.sub.n] at time t + [DELTA]t. The dot lines in Figures 2(b) and 2(c) are two degraded cohesive envelopes corresponding to [[delta].sup.1.sub.c] and [[delta].sup.2.sub.c], respectively. In both cases, if the loading level makes the reloading path move along the descending part of the degraded cohesive envelop, there is no damage accumulation that will occur and the material degradation is described by (1) implicitly.

The above CCZM has been implemented based on the user interface UMAT of the commercial code ABAQUS via the cohesive element [30]. Since the numerical simulations are performed according to the ABAQUS/Standard solver, both the updated cohesive traction and the material Jacobian matrix should be provided for constructing the internal nodal force vector and the cohesive stiffness matrix in UMAT. The standard Newton method is applied to solve the global nonlinear equations, and the default values for the solution control parameters are adopted. The global convergence rule is checked at the end of each time increment, such that the tolerance limit on the residual force can be fulfilled.

3. Model Validation

In the present study, the influence of model parameters on fatigue crack growth will be verified systematically; meanwhile, the predictive ability of the model for the high fatigue crack rates in Regime III will be validated experimentally. The material under investigation is a nickel-based superalloy IN718 which is widely used in aeroengine components due to its high yield stress, excellent fatigue resistance, and so on in severe conditions [31]. In modeling a CT specimen with the initial crack, the length is a = 25 mm and the thickness is B = 12.5 mm. The material properties are as follows: Young's modulus E = 200 GPa, Poisson's ratio v = 0.3, the initial yield stress [[sigma].sub.y] = 1196 MPa, and the plane strain fracture toughness [K.sub.Ic] = 125 MPa [m.sup.1/2]. The [J.sub.2] plasticity is applied to describe the plastic behavior of the alloy, and the relationship between true stress [sigma] and true plastic strain [[epsilon].sub.p] is described by the Holomon's equation as

[sigma] = K [([epsilon].sub.p]).sup.n], (6)

where the bulk modulus K = 1531 MPa and the strain hardening exponent n = 0.015. We note from the perspective of cohesive modeling for fatigue crack that the isotropic hardening is often assumed to simulate the material response approximately without knowing the exact cyclic deformation behavior of the bulk material [11, 28, 32]. However, the more accurate cyclic plasticity model should be determined experimentally as mentioned in [33]. Due to the symmetry, only the upper half of the specimen is modeled under the plane strain condition, as shown in Figure 3. The whole finite-element model consists of 1725 elements and 1872 nodes. The cohesive zone is predefined along the crack growth path using 240 COH2D4 elements which are uniformly distributed along the ligament of the specimen with the element length of 0.1 mm. The CPE4 elements are used for the rest of the regions.

3.1. Identification of Model Parameters. Based on previous descriptions, the model parameters can be divided into two types. The monotonic model parameters include cohesive strength [[sigma].sub.max,0], cohesive length [[delta].sub.0], and cohesive energy [[GAMMA].sub.0]. The cyclic model parameters are the accumulative length [d.sub.[summation]] and the cohesive endurance limit [[sigma].sub.f]. Omitting shear stress in the cohesive zone, only the normal traction is active under the mode I crack growth condition. The cohesive energy [[GAMMA].sub.0] is determined by integration of the monotonic cohesive envelope as follows:

[mathematical expression not reproducible], (7)

where the form of [T.sub.n] is represented by (1) with D = 0. The monotonic cohesive envelope can be constructed through the initial cohesive strength [[sigma].sub.max,0], the cohesive length [[delta].sub.0], and the cohesive energy [[GAMMA].sub.0]. Under the SSY conditions, the cohesive energy is related to the plain strain fracture toughness [K.sub.Ic] based on the Irwin's relation as follows:

[[GAMMA].sub.0] = (1 - [v.sup.2]) [K.sup.2.sub.Ic]/E. (8)

According to (8), the cohesive energy is [[GAMMA].sub.0] = 7.1 x [10.sup.4] J/[m.sup.2] in the present analysis. The cohesive strength [[sigma].sub.max,0] = 3900 MPa is used, which is nearly three times above the initial yield stress of IN718 according to the suggestion in [34], where the corresponding cohesive length [[sigma].sub.0] = 0.003 m m. Actually, the relationship between the initial cohesive strength and cohesive length can be attributed to the effect of the initial cohesive stiffness [k.sub.n]. The choice related to [k.sub.n] should guarantee that the cohesive zone does not influence the overall compliance before damage initiation and normally a higher value of [k.sub.n] is used. Since the cohesive energy is the work needed to create a unit area of two mating fracture surfaces, the resulting [[sigma].sub.u] = 0.017 mm is obtained according to (1) with D = 0.

Identification of the influence of the model parameters is performed using the CT specimen under constant loading amplitudes at the loading ratio R = 0.1, where R = [F.sub.min]/[F.sub.max] with [F.sub.max] and [F.sub.min] are the maximum and minimum of the uniaxial applied force F, respectively. The nominal stress intensity factor is applied to describe the stress field at crack-tip. The range of the stress intensity factor [DELTA]K for the CT specimens can be calculated through [35]:

[mathematical expression not reproducible], (9)

where [zeta] = a/W. This formulation is accurate within 0.5% error over the range 0.2 [less than or equal to] [zeta] [less than or equal to] 1.

The effect of the accumulative length [d.sub.[summation]] to the crack growth length curve is indicated in Figure 4(a), where the initial stress intensity factor [DELTA][K.sub.i] = 76 MPa [m.sup.1/2]. At a given crack growth length [DELTA]a, the lower value of [d.sub.[summation]] corresponds to the higher crack growth rate for the same value of [[sigma].sub.f] / [[sigma].sub.max,0]. In comparison with the number of loading cycles to the final failure between [d.sub.[summation]] = 20 [[delta].sub.0] and [d.sub.[summation]] = 25 [[delta].sub.0], the critical crack growth length of approx. 5.5 mm is achieved under both [d.sub.[summation]] = 20 [[delta].sub.0] and [d.sub.[summation]] = 25 [[delta].sub.0] at failure. Theoretically, the critical crack length should be reached at [d.sub.[summation]] = 10 [[delta].sub.0] also; however, the FE calculation is aborted since the high damage evolution rate associated with [d.sub.[summation]] = 10 [[delta].sub.0] is so fast that the mesh resolution of 0.1 mm at crack-tip cannot capture the crack growth rate before the final failure. The effect of [[sigma].sub.f]/[[sigma].sub.max,0] on the crack growth curves is shown in Figure 4(b) by comparing two different values of [[sigma].sub.f]/[[sigma].sub.max,0] at [d.sub.[summation]] = 10 [[delta].sub.0], where the higher crack growth rate is accompanied by the lower value of [[sigma].sub.f]/[[sigma].sub.max,0] at a given crack length. Based on Figure 4, we can conclude that the critical crack length is only determined by the monotonic model parameters together with the [DELTA]K level.

The effect of [d.sub.[summation]] on the da/dN versus [DELTA]K curve is shown in Figure 5(a), where the da/dN data are determined based on the corresponding a-N curves (see Figure 4). In order to examine the effects of cyclic model parameters on fatigue crack growth in the Paris' region, the numerical results of the lower loading level [F.sub.max] = 20 kN associated with [DELTA][K.sub.i] = 63 MPa [m.sup.1/2] is also displayed. The accumulative length [d.sub.[summation]] acts like the intercept coefficient C which operates in the Paris' law, that is, da/dN = C[([DELTA]K).sup.m]. It plays a role in increasing or decreasing crack growth rate in all portions at a certain value of [[sigma].sub.f]/[[sigma].sub.max,0]. In contrast, [[sigma].sub.f]/[[sigma].sub.max,0] has a tendency to modify the slope of the da/dN versus [DELTA]K curve as the parameter m of the Paris' law, as shown in Figure 5(b). However, the effects of the cyclic model parameters are also influenced by the monotonic model parameters; namely, the crack growth rate is accelerated inevitably as the [DELTA]K level approaches the critical value [DELTA][K.sub.c] = (1-R) [K.sub.Ic] in Regime III before the final failure. The above results reveal that the high fatigue crack growth rates at high [DELTA]K levels can be captured reasonably based on the current model and the uniform damage evolution equation can be applied in such simulations.

3.2. Comparison with Experiments. After the parameter identification, the next the predictive ability of the model will be examined using the CT specimen as displayed in Figure 3(a). The experimental data in validations are obtained by digitizing the results reported by Yamada and Newman [36] and Newman et al. [37], as shown in Figure 6, where the fatigue crack growth tests were performed by using the CT specimens and the SEN specimens at R = 0.1, respectively. The experimental results are fitted according to the Erdogan-Ratwani model [5]. For the CT specimens, the experimental points can be well approximated by the following equation:

da/dN = 2.554 x [10.sup.-5] [([DELTA]K - 8]).sup.1.887] /[DELTA][K.sub.c] - [DELTA] K, (10)

where da/dN is in mm/cycle and [DELTA]K is in MPa [m.sup.1/2]. For the SEN specimens, all points including the threshold region can be well described by the following equation:

da/dN = 1.782 x [10.sup.-5] [([DELTA]K - 5).sup.1.917]/[DELTA] [K.sub.c] - [DELTA] K. (11)

Three constant loading amplitudes, that is, [F.sub.max] = 17 kN, 18.5 kN, and 20 kN at R = 0.1, are applied to obtain the cyclic cohesive model parameters, where the corresponding initial stress intensity range are [DELTA][K.sub.i] = 54 MPa [m.sup.1/2], [DELTA] [K.sub.i] = 58 MPa [m.sup.1/2], and [DELTA][K.sub.i] = 63 MPa [m.sup.1/2], respectively. The fatigue crack growth simulations are performed to fit the experimental data associated with the moderate growth rates. As depicted in Figure 7, the fitting curve constructed by the aforementioned three [DELTA]K levels correlates the experiments very well in the linear Paris' region. The corresponding crack growth versus the applied loading cycle curves are displayed in Figure 8, where the fastest crack growth rate is achieved under the highest [DELTA]K level at a given crack length. The determined cyclic model parameters of IN718 are listed in Table 1.

For the fatigue crack growth simulation in Regime III, the constant loading amplitude is applied with [F.sub.max] = 24 kN at R = 0.1, which corresponds to the initial stress intensity factor range [DELTA] [K.sub.i] = 76 MPa [m.sup.1/2]. In order to capture the local response of each material point along the crack growth path, four material points 1, 2, 3, and 4 located at the crack growth path are used to trace the material responses associated with the high propagation rates, as illustrated in Figure 3(c). The distance between two adjoining monitory point is 1.2 mm.

The traction-separation responses during crack growth are indicated in Figure 9 for both points 1 and 4, where the material damage is obtained based on (2) for the two material points. The traction of point 1 approaches 87% of the initial cohesive strength [[sigma].sub.max,0] at the first loading cycle while the traction of point 4 is still small; namely, the traction of point 4 is only 10% of [[sigma].sub.max,0] at the first loading cycle, as shown in Figure 10. Except for point 1 loses its stress-carrying ability at the 207th cycle, other points have no damage accumulations at the first 548 loading cycles such that their original cohesive stiffnesses still exist, since there are no damage evolutions for these points (see Figure 11). The damage initiation for point 4 starts after 1500 loading cycles, where the maximal value of the traction increases (up to 80% of [[sigma].sub.max,0]) at first and then decreases with the applied loading cycles until the material point is completely damaged (see Figure 10(b)), and this phenomenon is also true for points 2 and 3.

During the crack growth process, the distributions of the damage zone, the von-Mises stress, and the equivalent plastic strain under both the moderate and the high levels of [DELTA]K are compared in Figure 12. The fatigue crack propagates from length a = 25.5 mm to 29.5 mm by measuring the length of the complete damage zone of the two levels of [DELTA]K. For two cases, the von-Mises stress ahead of the current crack-tip is greater than the initial yield stress 1196 MPa and the obvious plastic zone is created gradually in the crack wake during crack growth. These results show that the material responses of both the cohesive zone and the surrounding bulk material are simulated reasonably.

Figure 13 summarizes the numerical and experimental crack growth results under constant amplitude loading. The simulation associated with [F.sub.max] = 24 kN is the predicting result based on the determined model parameters listed in Table 1. The whole da/dN versus [DELTA]K curves are assembled through these four different [DELTA][K.sub.i] levels. With increment of the crack length, da/dN becomes very high and approaches instability limited by the [DELTA][K.sub.c] at failure (see Figure 13). As indicated in Figure 1, fatigue crack growth in Regime II is largely influenced by the loading ratio and the frequency in the air; therefore, fatigue crack growth with moderate rates is irrelevant to the geometries of specimens. At high propagation rates as [DELTA]K approaches to [DELTA][K.sub.c], there is a good agreement between the prediction and the experiment of the SEN specimens, while a deviation between the simulation and the experiment of the CT specimens exists. Actually, the difference in experiment data also occurs between the CT and SEN specimens in Regime III due to the occurrence of ductile fracture. The ductile fracture process of metallic materials involves nucleation, growth, and coalescence of microvoids, which is greatly influenced by the stress state and may not be described, for example, by (1) simply, and the local constraint effect accompanied by the specimen geometry also affects the crack growth resistance. Nevertheless, based on the cyclic model parameters determined by the experimental data in Regime II, the numerical result captures the accelerated fatigue crack growth phenomenon of the experimental data in Regime III associated with the high [DELTA]K levels.

4. Conclusions

In order to overcome the limitations of fracture mechanics-based models, for example, the Paris' law or the [DELTA]J-integral-based model, for predicting elastic-plastic fatigue crack growth rates, the cyclic cohesive model is verified systematically on the fatigue crack growth of the IN718 alloy from moderate to high [DELTA]K levels. Since ductile fracture occurs as the [DELTA]K level approaches the fracture toughness, the model parameters determined from the experimental data at the moderate [DELTA]K levels cannot predict the fatigue crack growth rates at the high [DELTA]K levels for all the specimens exactly, and the local constraint effect should be considered during the damage evolution in Regime III.

The simulations reveal that the novel cohesive envelope combined with one damage evolution equation can also be used to capture the accelerated fatigue crack growth behavior at the high [DELTA]K levels, and there is no need to distinguish between the monotonic damage and the fatigue damage as introduced in [29]. The crack growth behaviors are largely affected by the accumulative length [d.sub.[SIGMA]] and the cohesive endurance limit [[sigma].sub.f]; meanwhile, the effects of the cyclic model parameters are also influenced by the monotonic model parameters, that is, cohesive energy and the initial cohesive strength as the [DELTA]K level approaches the critical value corresponding to the fracture toughness.

After the model parameter identification and the experimental validation, we confirm that the current model has the capacity to simulate the elastic-plastic fatigue crack growth behavior of IN718 at SSY, especially in Regime III. Furthermore, it remains for the next work to extend the application of the CCZM for simulating fatigue crack growth of other alloys at LSY.

https://doi.org/10.1155/2018/4048386

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was made possible through research funding provided by the National Natural Science Foundation of China (no. 11502204), the National Key Research and Development Program of China (no. 2016YFB0701303), the Natural Science Basic Research Plan in Shaanxi Province of China (no. 2017JM1023), and the China Postdoctoral Science Foundation funded project (no. 2017M623234).

References

[1] J. Lemaitre and R. Desmorat, Engineering Damage Mechanics: Ductile, Creep, Fatigue and Brittle Failures, Springer-Verlag, Berlin Heidelberg, 2005.

[2] Y. Jiang and M. Feng, "Modeling of fatigue crack propagation," Journal of Engineering Materials and Technology, vol. 126, no. 1, pp. 77-86, 2004.

[3] P. Paris, M. Gomez, and W. Anderson, "A rational analytic theory of fatigue," Trend in Engineering, vol. 13, pp. 9-14, 1961.

[4] R. O. Ritchie, "Mechanisms of fatigue-crack propagation in ductile and brittle solids," International Journal of Fracture, vol. 100, no. 1, pp. 55-83, 1999.

[5] F. Erdogan and M. Ratwani, "Fatigue and fracture of cylindrical shells containing a circumferential crack," International Journal of Fracture Mechanics, vol. 6, no. 4, pp. 379-392, 1970.

[6] S. M. Beden, S. Abdullah, and A. Ariffin, "Review of fatigue crack propagation models for metallic components," European Journal of Scientific Research, vol. 28, pp. 364-397, 2009.

[7] N. Dowling and J. Begley, "Fatigue crack growth during gross plasticity and the J-integral," Mechanics of crack growth, vol. 590, pp. 82-103, 1976.

[8] V. Shlyannikov, N. Boychenko, A. Tumanov, and A. Fernandez-Canteli, "The elastic and plastic constraint parameters for three-dimensional problems," Engineering Fracture Mechanics, vol. 127, pp. 83-96, 2014.

[9] V. Shlyannikov, A. Tumanov, A. Zakharov, and A. Gerasimenko, "Surface flaws behavior under tension, bending and biaxial cyclic loading," International Journal of Fatigue, vol. 92, pp. 557-576, 2016.

[10] M. Metzger, T. Seifert, and C. Schweizer, "Does the cyclic J-integral AJ describe the crack-tip opening displacement in the presence of crack closure?," Engineering Fracture Mechanics, vol. 134, pp. 459-473, 2015.

[11] A. de-Andres, J. Perez, and M. Ortiz, "Elastoplastic finite element analysis of three-dimensional fatigue crack growth in aluminum shafts subjected to axial loading," International Journal of Solids and Structures, vol. 36, no. 15, pp. 2231-2258, 1999.

[12] K. L. Roe and T. Siegmund, "An irreversible cohesive zone model for interface fatigue crack growth simulation," Engineering Fracture Mechanics, vol. 70, no. 2, pp. 209-232, 2003.

[13] D. S. Dugdale, "Yielding of steel sheets containing slits," Journal of the Mechanics and Physics of Solids, vol. 8, no. 2, pp. 100-104, 1960.

[14] G. I. Barenblatt, "The mathematical theory of equilibrium cracks in brittle fracture," Advances in Applied Mechanics, vol. 7, pp. 55-129, 1962.

[15] J. Tryding and M. Ristinmaa, "Normalization of cohesive laws for quasi-brittle materials," Engineering Fracture Mechanics, vol. 178, pp. 333-345, 2017.

[16] O. Shor and R. Vaziri, "Application of the local cohesive zone method to numerical simulation of composite structures under impact loading," International Journal of Impact Engineering, vol. 104, pp. 127-149, 2017.

[17] V. Tvergaard and J. Hutchinson, "The relation between crack growth resistance and fracture process parameters in elastic-plastic solids," Journal of the Mechanics and Physics of Solids, vol. 40, no. 6, pp. 1377-1397, 1992.

[18] C. R. Chen, O. Kolednik, J. Heerens, and F. D. Fischer, "Three-dimensional modeling of ductile crack growth: cohesive zone parameters and crack tip triaxiality," Engineering Fracture Mechanics, vol. 72, no. 13, pp. 2072-2094, 2005.

[19] V. Uthaisangsuk, U. Prahl, and W. Bleck, "Modelling of damage and failure in multiphase high strength DP and TRIP steels," Engineering Fracture Mechanics, vol. 78, no. 3, pp. 469-486, 2011.

[20] H. Li, C. Li, and H. Yuan, "Prediction of fatigue crack growth retardation using a cyclic cohesive zone model," Archive of Applied Mechanics, vol. 87, no. 6, pp. 1061-1075, 2017.

[21] K. Park and G. Paulino, "Cohesive zone models: a critical review of traction-separation relationships across fracture surfaces," Applied Mechanics Reviews, vol. 64, no. 6, article 060802, 2013.

[22] O. Nguyen, E. A. Repetto, M. Ortiz, and R. A. Radovitzky, "A cohesive model of fatigue crack growth," International Journal of Fracture, vol. 110, no. 4, pp. 351-369, 2001.

[23] S. Maiti and P. Geubelle, "A cohesive model for fatigue failure of polymers," Engineering Fracture Mechanics, vol. 72, no. 5, pp. 691-708, 2005.

[24] A. Ural, V. Krishnan, and K. Papoulia, "A cohesive zone model for fatigue crack growth allowing for crack retardation," International Journal of Solids and Structures, vol. 46, no. 11-12, pp. 2453-2462, 2009.

[25] H. Li and H. Yuan, "Cohesive zone modelling of low cycle fatigue cracks in cracked and notched specimens," Fatigue & Fracture of Engineering Materials & Structures, vol. 36, no. 12, pp. 1246-1257, 2013.

[26] S. Jimenez and R. Duddu, "On the parametric sensitivity of cohesive zone models for high-cycle fatigue delamination of composites," International Journal of Solids and Structures, vol. 82, pp. 111-124, 2016.

[27] X. Cao and M. Vormwald, "Application of a new cohesive zone model in low cycle fatigue," Solids and Structures, vol. 2, pp. 31-40, 2013.

[28] D. Jha and A. Banerjee, "A cohesive model for fatigue failure in complex stress-states," International Journal of Fatigue, vol. 36, no. 1, pp. 155-162, 2012.

[29] H. Li, H. Yuan, and X. Li, "Assessment of low cycle fatigue crack growth under mixed-mode loading conditions by using a cohesive zone model," International Journal of Fatigue, vol. 75, pp. 39-50, 2015.

[30] ABAQUS, Version 6.11, User Subroutines Reference Manual, Dassault Systemes Simulia Corp., Providence, RI, USA, 2011.

[31] E. Lertora, C. Mandolfino, and C. Gambaro, "Mechanical behaviour of Inconel 718 thin-walled laser welded components for aircraft engines," International Journal of Aerospace Engineering, vol. 2014, Article ID 721680, 9 pages, 2014.

[32] P. Beaurepaire and G. I. Schueller, "Modeling of the variability of fatigue crack growth using cohesive zone elements," Engineering Fracture Mechanics, vol. 78, no. 12, pp. 2399-2413, 2011.

[33] J. L. Chaboche, "A review of some plasticity and viscoplasticity constitutive theories," International Journal of Plasticity, vol. 24, no. 10, pp. 1642-1693, 2008.

[34] A. Needleman, "Some issues in cohesive surface modeling," Procedia IUTAM, vol. 10, pp. 221-246, 2014.

[35] ASTM International, Standard Test Method for Measurement of Fatigue Crack Growth Rates, 2011.

[36] Y. Yamada and J. C. Newman, "Crack closure under high load-ratio conditions for Inconel-718 near threshold behavior," Engineering Fracture Mechanics, vol. 76, no. 2, pp. 209-220, 2009.

[37] J. Newman, S. Smith, S. Willard, and R. Piascik, A Comparison of Weld-Repaired and Base Metal for Inconel 718 and CRES 321 at Cryogenic and Room Temperatures, 2004, NASA/TM-2004-213, 253 ARL-TR-3266.

Huan Li (iD), (1,2) Jinshan Li, (2) and Huang Yuan (3)

(1) School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi'an, China

(2) State Key Laboratory of Solidification Processing, Northwestern Polytechnical University, Xi'an, China

(3) School of Aerospace Engineering, Tsinghua University, Beijing, China

Correspondence should be addressed to Huan Li; lihuan@nwpu.edu.cn

Received 8 September 2017; Accepted 6 December 2017; Published 1 February 2018

Academic Editor: Roberto G. Citarella

Caption: Figure 1: The typical fatigue crack growth rate with regard to the stress intensity factor curve for metallic materials, which is distinguished by three regimes. The corresponding influence factors and the crack propagation velocity at 50 Hz under different levels of the fatigue crack growth rate are also illustrated [4].

Caption: Figure 2: Illustration of the loading-unloading-reloading path together with (a) the cyclic loading spectrum at (b) [[delta].sub.n] < [[delta].sub.0] (first half cycle) and (c) [[delta].sub.n] > [[delta].sub.0] (first half cycle): loading with opening of the cohesive surface (1 [right arrow] 2); unloading (2 [right arrow] 3); reloading with the damage accumulation (3 [right arrow] 4); and subsequent partial unloading (4 [right arrow] 5). The dot lines represent the degraded cohesive envelopes.

Caption: Figure 3: (a) Sketch of a CT specimen (initial crack length a = 25 mm, width W = 50 mm, and thickness B = 12.5 mm) and (b) the FE mesh (a half of the specimen is modeled), and (c) points 1, 2, 3, and 4 correspond to the four integration points located along the crack growth path (the distance between two adjacent points is 1.2 mm), and the arrow indicates the initial crack-tip of the CT specimen.

Caption: Figure 4: The effects of (a) the accumulative length [d.sub.[summation]] and (b) [[sigma].sub.f]/[[sigma].sub.max,0] on the crack growth curve. The constant loading amplitude is carried out with [F.sub.max] = 24 kN at R = 0.1, where the corresponding initial stress intensity factor range [DELTA] [K.sub.i] = 76 MPa [m.sup.1/2] for the CT specimen.

Caption: Figure 5: The effects of (a) the accumulative length [d.sub.[summation]] and (b) [[sigma].sub.f]/[[sigma].sub.max,0] on the fatigue crack growth rate. The constant loading amplitude is carried out with [F.sub.max] = 20 kN and 24 kN at R = 0.1, where the corresponding initial stress intensity factor range [DELTA][K.sub.i] = 63 MPa [m.sup.1/2] and [DELTA][K.sub.i] = 76 MPa [m.sup.1/2] for the CT specimen, respectively.

Caption: Figure 6: Plot of the experimental result and the corresponding fitting curve by using the Erdogan-Ratwani model of the IN718 alloy for (a) the CT specimens and (b) the SEN specimens.

Caption: Figure 7: Plot of the experimental results and the corresponding fitting curves of the CCZM in the Paris' region.

Caption: Figure 8: The crack growth curves with varying [DELTA][K.sub.i] under the constant amplitude loading at R = 01.

Caption: Figure 9: The traction versus separation curves are compared between (a) point 1 and (b) point 4, and the constant loading amplitude is carried out with [F.sub.max] = 24 kN at R = 0.1, which corresponds to the initial stress intensity factor range [DELTA][K.sub.i] = 76 MPa [m.sup.1/2] for the CT specimen.

Caption: Figure 10: The traction against the number of applied cycles is compared between (a) point 1 and (b) point 4, and the constant loading amplitude is carried out with [F.sub.max] = 24 kN at R = 0.1, which corresponds to the initial stress intensity factor range [DELTA][K.sub.i] = 76 MPa [m.sup.1/2] for the CT specimen.

Caption: Figure 11: (a) The damage evolution and (b) the cohesive stiffness with the number of applied cycles are compared from point 1 to point 4, and the constant loading amplitude is carried out with [F.sub.max] = 24 kN at R = 0.1, which corresponds to the initial stress intensity factor range [DELTA][K.sub.i] = 76 MPa [m.sup.1/2] for the CT specimen.

Caption: Figure 12: Contour plots during crack growth: (a, b) damage evolution; (c, d) von-Mises stress; (e, f) equivalent plastic strain (PEEQ). The corresponding crack length a is (a, c, e) 25.5 mm with the moderate [DELTA]K = 77 MPa [m.sup.1/2]; (b, d, f) 29.5 mm with the high [DELTA]K =102 MPa [m.sup.1/2].

Caption: Figure 13: Comparison between experiments and the proposed CCZM under the constant force amplitude with R = 0.1.

Table 1: Model parameters of IN718 used in simulations. [[sigma].sub.max,0] [[delta].sub.0] [[delta].sub.u] (MPa) (mm) (mm) 3900 0.003 0.017 [[sigma].sub.f] [d.sub.[summation]]] (MPa) (MPa) 1847 0.03

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Li, Huan; Li, Jinshan; Yuan, Huang |

Publication: | International Journal of Aerospace Engineering |

Date: | Jan 1, 2018 |

Words: | 7015 |

Previous Article: | New Integrated Guidance and Control of Homing Missiles with an Impact Angle against a Ground Target. |

Next Article: | Efficient Training Data Generation for Reduced-Order Modeling in a Transonic Flight Regime. |

Topics: |