Next Article in Journal
Correction: Speranza, V. et al. Hierarchical Structure of iPP During Injection Molding Process with Fast Mold Temperature Evolution. Materials 2019, 12, 424
Next Article in Special Issue
A Novel Approach to Describe the Time–Temperature Conversion among Relaxation Curves of Viscoelastic Materials
Previous Article in Journal
Role of the La/K Compositional Ratio in the Properties of Waveguides Written by Fs-Laser Induced Element Redistribution in Phosphate-Based Glasses
Previous Article in Special Issue
Influence of the Elastoplastic Strain on Fatigue Durability Determined with the Use of the Spectral Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study on the Variability of Plastic CTOD

by
Pedro André Prates
*,
Armando Eusébio Marques
,
Micael Frias Borges
,
Ricardo Madeira Branco
and
Fernando Ventura Antunes
Centre for Mechanical Engineering, Department of Mechanical Engineering, Materials and Processes (CEMMPRE), University of Coimbra, 3030-788 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Materials 2020, 13(6), 1276; https://doi.org/10.3390/ma13061276
Submission received: 13 January 2020 / Revised: 5 March 2020 / Accepted: 9 March 2020 / Published: 11 March 2020
(This article belongs to the Special Issue Probabilistic Mechanical Fatigue and Fracture of Materials)

Abstract

:
This paper presents a numerical study on the influence of material parameters and loading variability in the plastic crack tip opening displacement (CTOD) results. For this purpose, AA7050-T6 was selected as reference material and a middle-cracked tension specimen geometry was considered. The studied input parameters were the Young’s modulus, Poisson’s ratio, isotropic and kinematic hardening parameters and the maximum and minimum applied loads. The variability of the input parameters follows a Gaussian distribution. First, screening design-of-experiments were performed to identify the most influential parameters. Two types of screening designs were considered: one-factor-at-a-time and fractional factorial designs. Three analysis criteria were adopted, based on: main effect, index of influence and analysis of variance. Afterwards, metamodels were constructed to establish relationships between the most influential parameters and the plastic crack tip opening displacement (CTOD) range, based on two types of designs: Face-Centered Central Composite Design and Box-Behnken design. Finally, the metamodels were validated, enabling the expeditious evaluation of the variability in the plastic CTOD range; in addition, the variability in the fatigue crack growth rate was also evaluated.

1. Introduction

Fatigue is the main failure mechanism in components submitted to cyclic loads. Design against fatigue based on damage tolerance approach assumes the presence of intrinsic defects, such as cracks. In fact, technological processes like casting, machining, welding and additive manufacturing are known to produce small defects. Accurate assessments of fatigue crack growth (FCG) rates are needed to define the time between inspections. Several models have been proposed in the literature to quantify the FCG rate as a function of loading and/or material parameters (e.g., [1,2,3]). These models are deterministic, i.e., assume that there is a well-defined relation between these parameters and FCG rate. However, the variability of the material, loading, geometry, temperature, etc., introduces uncertainty in the results. In this context, there are two main procedures to overcome the influence of sources of uncertainties [4]: (i) the use of safety factors, which is the conservative solution usually adopted to accommodate uncertainties, and (ii) performing stochastic analysis, which is often supported by probabilistic fatigue models (e.g., [5,6]), metamodeling techniques (e.g., [7,8,9,10,11,12,13]) and more recently by digital twin approaches (e.g., [14,15]). In the context of stochastic fatigue analysis, a probabilistic fracture mechanics approach to predict the fatigue life of aircraft wing attachment bulkheads was proposed by White [16]. A sensitivity analysis was performed to quantify the effect of parameters that significantly influence crack growth in metallic airframes. The studied parameters included: equivalent pre-crack size, crack growth rate, fracture toughness of the material and maximum load per flight. The variability in the total life is most sensitive to the variability in the equivalent pre-crack size, and least sensitive to the variability in the fracture toughness. In other works, FCG tests were performed on cracked specimens made of 2024-T3 aluminum alloy, illustrating the scatter in the FCG rate due to material inhomogeneity [17,18]. A predictive procedure was developed by Bahloul et al. [19], using the finite element method coupled with Monte Carlo, to evaluate the residual FCG life of a 7075-T6 aluminum alloy cracked lug component, considering the stochastic behaviour of the material parameters and the crack tip stress field. The proposed procedure showed a good ability in improving the deterministic FCG life evaluation. An uncertainty quantification methodology was proposed by Sankararaman et al. [7], to predict the crack growth as a function of the number of load cycles in a cylindrical structural component, with complex geometry and subjected to multi-axial loading with variable amplitude. The proposed methodology, which is supported by a Gaussian process metamodel trained via finite element analysis, takes into account sources of uncertainty such as loading conditions and material variability. Wentao et al. [8] proposed an approach to assess the probabilistic life of mixed mode FCG in mechanical components, by coupling finite element analysis and a Kriging-based Monte Carlo metamodel. The use of a Kriging-based metamodel for the probabilistic FCG analysis is shown to be more efficient than the Response Surface Methodology. A probabilistic fatigue life prediction model for multiple site damage in structural components was developed by Kim et al. [9], which resort to an extended finite element method coupled with a Gaussian Process metamodel for computational efficiency. Metamodels were developed by Vélez et al. [10], to relate the crack tip opening displacement (CTOD) with material properties in high-thickness offshore steel welded joints. The studied material properties (hardness, chemical composition, toughness and microstructural morphology) were shown to be significantly related with the CTOD, also being simpler and cheaper to measure and more available than the CTOD.
The aim of this paper is to carry out a numerical analysis on the influence of material parameters and loading variability in the variability of the plastic CTOD results obtained from FCG tests on middle-cracked tension (M(T)) specimens. Antunes et al. [20] developed a material law relating FCG rate with the plastic CTOD (δp). The FCG rate was experimentally measured in standard M(T) samples, while δp was numerically predicted in models that simulate the experimental work regarding material properties, sample geometry and imposed loadings. It is required a certain amount of crack propagation in the numerical models for each crack length studied, in order to stabilize the cyclic plastic deformation and the crack closure level. The plastic CTOD is measured only at the end of crack propagation, and is not expected to be affected by the crack growth rate imposed artificially in the numerical model. This approach assumes that plastic deformation at the crack tip is the main mechanism behind crack tip and that the plastic CTOD is the crack driving force. Other authors proposed the use of total CTOD [21,22], but Vasco-Olmo et al. [23] proved that the plastic CTOD is the crack driving force. The da/dN–δp model is supposed to be a material law, which can be used to estimate da/dN in other situations. In this work, the numerical models were changed regarding the material and loading parameters, and the da/dN–δp model was used to estimate the FCG rate.
Following this introduction, the paper is organized as follows: Section 2 describes the numerical model of the FCG tests; Section 3 presents the results and discussion. Firstly, sensitivity analyses are performed to identify the most influential parameters on the plastic CTOD results. Afterwards, metamodels are constructed with the most influential parameters, to expeditiously evaluate the variability of the plastic CTOD. Finally, the proposed metamodels are validated and compared. Section 4 presents the conclusions.

2. Numerical Model

Numerical simulations of FCG tests were performed on a standard middle-cracked tension (M(T)) sample geometry, as shown in Figure 1. A straight crack was defined, with an initial size ao equal to 5 mm (ao/W = 0.083). Only 1/8 of the low cycle fatigue test was simulated, due to sample geometry symmetries. Plane strain state was modeled as illustrated in Figure 1b.
The finite element (FE) model of the M(T) sample is composed of 6639 trilinear solid elements, with 13586 nodes. A mesh refinement was performed near the crack tip, using 8 × 8 μm2 elements. Crack propagation was simulated by successive debonding of nodes at minimum imposed load. Each crack increment corresponds to one FE and two load cycles were applied before each crack propagation. In each cycle, the crack propagated uniformly over the thickness by releasing both current crack front nodes. 200 load cycles (i.e., 100 crack propagations) were performed, which corresponds to a total crack propagation of Δa = (100 − 1) × 8 μm = 792 μm. The first two load cycles were applied without crack increment, i.e., at a = 5 mm. The aim was the stabilization of crack opening values, since a transient behavior is observed as the residual plastic wake is formed. Additionally, the residual plastic wake may produce crack closure, which affects the crack tip fields. Again, the plastic CTOD is measured at the end of crack propagation and is not expected to be affected by the artificial crack growth rate imposed in the numerical model. The applied remote stresses σ were obtained by dividing the imposed loads by the cross-section area A (= 30 × 0.1 mm2).
The numerical simulations were performed with DD3IMP in-house FE solver [24]. The average simulation time for the FCG tests is about 18 h (Intel® Core™i9–7900X 10-Core processor @ 3.3 GHz). The material elastic behavior is considered isotropic and is described by the generalized Hooke’s law, where the Young’s modulus (E) and the Poisson’s ratio (ν) are the elastic parameters. The plastic behavior is described by the von Mises yield criterion, coupled with the Armstrong-Frederick kinematic hardening law. The von Mises yield criterion is:
3 2 ( σ X ) : ( σ X ) Y 0 = 0
where σ is the deviatoric Cauchy stress tensor, X is the deviatoric backstress tensor and Y 0 is the initial yield stress. The Armstrong-Frederick kinematic hardening law is described as follows [25]:
X ˙ = C X [ X Sat σ X σ ¯ X ] ε ¯ ˙ p
where X ˙ is the backstress rate, representing the translational rate of the center of the yield surface due to plastic deformation, σ ¯ is the equivalent stress and ε ¯ ˙ p is the equivalent plastic strain rate; and the material parameters C X and X Sat represent, respectively, the saturation rate and saturation value of the kinematic hardening (exponential) law.
The material used in this study is a 7050-T6 aluminum alloy, whose elastoplastic cyclic behavior was modeled in a previous work [20]; the values of the material parameters E, ν, Y 0 , C X , X Sat that describe the elastoplastic behavior of AA7050-T6 were used as reference. In this context, a Gaussian distribution is used for describing the variability of the material parameters and the maximum and minimum imposed loads per cycle, Fmax and Fmin. The mean (μ) and standard deviation (SD) values of each parameter are detailed in Table 1, where the mean values of the material parameters correspond to the reference ones (see [20]); a coefficient of variation CV = 5% was assumed for each parameter, where CV = SD/μ.
Figure 2a presents the reference numerical simulation results of CTOD vs. applied stress, obtained from the mean values of the input parameters in Table 1 during the 200th load cycle. The CTOD is obtained at the first node at the left of the crack tip (8 μm from the crack tip, as schematized in Figure 2a). At minimum applied stress (point A) the crack is closed, and therefore the CTOD is equal to zero. The increase in the stress value opens the crack at point B. Afterwards, there is a linear region (point B to point C) associated with the material elastic behavior. Plastic deformation initiates after point C, and increases up to the maximum stress at point D. The dashed line shows the elastic CTOD (CTODe); the plastic CTOD (CTODp), is obtained by subtracting the elastic CTOD from the total CTOD. Figure 2b presents the evolution of plastic CTOD with applied stress during the 200th load cycle. Between points C and D, there is a progressive increase of CTODp up to the maximum stress. The range of CTODp, δp, shown in Figure 2b, is assumed to control FCG; accordingly, the variability analysis presented in the subsequent sections will focus on δp.

3. Results and Discussion

3.1. Sensitivity Analysis

The sensitivity analysis enables excluding the least influential input parameters on the δp results, to reduce computational effort in the subsequent metamodeling step. For this purpose, numerical simulations were carried out from two types of screening design of experiments (DOE): One-factor-at-a-time (OFAT) and fractional factorial design (FFD) approaches [26]. Three levels were assumed for the input parameters: the mean value (µ), the 2.5th (P2.5) and the 97.5th (P97.5) percentiles, which correspond to a 95% confidence interval. Table 2 and Table 3 show the δp results obtained from the numerical simulations under the OFAT and FFD screening approaches, respectively.
The influence of each input parameter on the plastic CTOD range was quantified using three analysis criteria [26]: (i) Main Effect, (ii) Index of Influence and (iii) Analysis of Variance (ANOVA). This enables assessing the influence of the analysis criterion and screening DOE on the input parameters sensitivity.
The Main Effect of a given input parameter is given by:
Main   Effect = δ p P 97.5 δ p P 2.5 0.5 × s
where s is the total amount of numerical simulations of the screening DOE under analysis, having input parameter values at levels P2.5 and/or P97.5 (i.e., (s = 14 for OFAT and s = 16 for FFD); δ p P 97.5 and δ p P 2.5 are the sum of plastic CTOD range values obtained at the 97.5th and the 2.5th percentiles, respectively. According to the Main Effect criterion, the input parameters are assumed influential if the absolute value of their main effect is higher than the average main effect of the seven input parameters, for each DOE.
The Index of Influence of a given parameter is given by:
Index   of   Influence = | δ p P 2.5 0.5 × s δ p N | + | δ p P 97.5 0.5 × s δ p N | 2 × δ p N
where δ p N is the nominal value of the plastic CTOD range, and depends on the type of DOE: In case of OFAT, it is obtained from the Simulation 1 in Table 2 (i.e., when considering the mean values of the input parameters shown in Table 1); in case of FFD, it is equal to the average of the δp values from Table 3. According to the Index of Influence criterion, the input parameters are considered influential when the value of their Index of Influence is higher than the average Index of Influence of the seven input parameters, for each DOE.
ANOVA is used to check the significance of each input parameter, by determining the p-value for the F-test at a 95% confidence interval. ANOVA starts with the main effect of each input parameter (Equation (4)), for determining the corresponding sum of squares, SS:
SS = s ( Main   Effect 2 ) 2
In the next step, the mean of squares, MS, is calculated for each input parameter:
MS = SS DF
where DF is the number of degrees of freedom per input parameter, DF = l − 1, where l is the number of levels of variation considered for each input parameter (excluding the mean, µ). Since the number of levels considered is l = 2 (i.e., P2.5 and P97.5), then DF = 1. The sum of squares for error, SSE, is:
SSE = i = 1 s ( δ p i δ ^ p i ) 2
where δ p i is the plastic CTOD range obtained from DOE simulation i and δ ^ p i is the corresponding predicted plastic CTOD range value, obtained by the following linear relationship:
δ ^ p i = α 0 + α 1 E i + α 2 ν i + α 3 Y 0 i + α 4 C X i + α 5 X Sat i + α 6 F max i + α 7 F min i ,
where αj are coefficients, with j = 0, …, 7. The values of the coefficients were obtained by minimization of Equation (7) resorting to the Generalized Reduced Gradient non-linear optimization algorithm [27]. The mean squared error MSE is given by:
MSE = SSE s k 1
where k is the number of input parameters under analysis (k = 7). The F-ratio is given by:
F Ratio = MS MSE
The F-ratio enables obtaining the p-value at a 95% confidence interval. For the 95% confidence interval, the input parameters are assumed influential if the p-value is less than 0.05.
The sensitivity analysis results for OFAT and FFD are respectively indicated in Table 4 and Table 5. Whatever the analysis (Main Effect, Index of Influence and ANOVA) and type of screening DOE (OFAT and FFD), the input parameters E, Y0 and Fmax are consistently shown to be the most influential in the δp results.

3.2. Metamodeling

The variability in the δp results was evaluated resorting to Response Surface Methodology (RSM) metamodels, constructed with the most influential input parameters, E, Y0 and Fmax, to reduce the computational effort. In this context, a quadratic polynomial model is adopted to establish the following relationship between E, Y0, Fmax and δ p :
δ p i RSM = β 0 + β 1 E i + β 2 Y 0 i + β 3 F max i + β 4 E i Y 0 i + β 5 E i F max i + β 6 Y 0 i F max i + β 7 E i 2 + β 8 Y 0 i 2 + β 9 F max i 2
where δ p i RSM is the value of plastic CTOD range predicted by the RSM model for RSM simulation i and βj are fitting RSM coefficients, with j = 0, …, 9. Equation (11) can be also written as a system of linear equations, as follows:
δ p i = H i j β j + ε i
where δ p i (i = 1, …, n) is the vector of plastic CTOD range measurements obtained from n RSM simulations, ε i is an error term and H i j is the linear system matrix,
H i j = [ 1 E 1 Y 0 1 F max 1 E Y 0 1 E F max 1 Y 0 F max 1 E 1 2 Y 0 1 2 F max 1 2 1 E n Y 0 n F max n E Y 0 n E F max n Y 0 F max n E n 2 Y 0 n 2 F max n 2 ]
RSM simulations were performed based on two types of designs: Face-Centered Central Composite Design (FCCCD) and Box-Behnken design (BBD); the aim is to assess the influence of the choice of RSM metamodel on predicting the variability in the δp results. Figure 3 shows the map of points (i.e., sets of input parameters) considered in the FCCCD and BBD design space, each point corresponding to one numerical simulation. Three levels of variation are considered for each input parameter: P2.5, µ and P97.5 (see Table 1); these levels are referred in Figure 3 as −1, 0, 1, respectively. The FCCCD design uses vertex points to assess the plastic CTOD range under extreme combinations between the three input parameters; as such, the resulting response surface is expected to describe extreme (and highly unlikely) behavior more accurately than that obtained using Box-Behnken (which lacks vertex points, having edge points instead). On the other hand, the Box-Behnken design is a computationally cheaper alternative to FCCCD, requiring fewer numerical simulations. The RSM numerical simulations assumed the mean values (µ) concerning the least influential parameters (ν, CX, XSat and Fmin—see Table 1).
Table 6 presents the least-squares solution for the RSM coefficients associated with the FCCCD and BBD designs. Whatever the design, FCCCD or BBD, the number of numerical simulations n is greater than the number of RSM coefficients β (i.e., n > j); in this case, the least-squares solution is the one that minimizes the sum of squared error at the design points, ε i = ( δ p i H i j β j ) 2 , given by [28]:
β j = ( H j i H i j ) 1 H j i δ p i
Table 7 and Table 8 compare the δp responses obtained from the FCCCD or BBD simulations with those predicted by the corresponding RSM metamodels; the metrics R2 and relative root mean squared error (RRMSE) were chosen to quantify the fitting performance of the design points. In general, the obtained metamodels are able to adequately describe the AA7050-T6 plastic CTOD range measurements at the design points; also, the metamodel accuracy is best when considering the Box-Behnken design (R2 = 0.9996, RRMSE = 0.4%) than FCCCD (R2 = 0.9927, RRMSE = 2.0%).

3.3. Metamodel Validation

The R2 analysis shown in Table 7 and Table 8 enabled a first assessment of the predictive ability of the RSM metamodels. At this stage, 60 random numerical simulations are performed to check if the response surfaces still represent a proper approximation of the plastic CTOD range within the range of variation of the main input parameters. Figure 4 presents the correspondence between numerical simulation and response surface δp results, for the FCCCD and BBD based metamodels. In this figure, each point corresponds to one numerical simulation. The scatter of the random points around the straight line is assessed with R2; the straight line indicates the ideal metamodel, in which the predicted response is equal to that obtained by numerical simulation. In general, both FCCCD and BBD metamodeling approaches provide accurate predictions of δp, with RRMSE equal to 1.7% (FCCCD) and 0.9% (BBD), although dispersion is more evident for FCCCD (with R2 = 0.9844—see Figure 4a) than for BBD (with R2 = 0.9963—see Figure 4b). This corroborates the adequacy of the proposed RSM metamodels for establishing the relationship between the most influential inputs and δp.

3.4. Variability Analysis Based on Metamodels

The RSM metamodels are now used to expeditiously predict the variability in the δp results. Accordingly, 100,000 random sets of input parameters E, Y0 and Fmax were generated. The values of each set of input parameters were used as input for the RSM models (i.e., see Equation (11), with the values of the RSM coefficients shown in Table 6), to predict the δp variability. Figure 5 presents the δp variability results predicted by the RSM metamodels based on FCCCD and BBD designs. Both histograms present the predicted mean δp values as well as percentiles P2.5 and P97.5, enabling the definition of a 95% confidence interval. According to the figure, the RSM metamodels predict the variability in the plastic CTOD range in a similar way, both presenting right-skewed distributions. Such skewed distributions contrast with the normal distribution, which was assumed for describing the variability of the input parameters.
Table 9 presents descriptive statistics of the δp histograms shown in Figure 5. The predicted mean δp value (approx. 0.340 µm) is close to the δp value obtained from the reference simulation (=0.334 µm, see Simulation 1 in Table 1), which is assumed to be the exact solution under a deterministic analysis. However, from the perspective of the current stochastic analysis, it is possible to estimate not only the mean δp value, but also the associated variance; accordingly, the predicted variability in the δp results follows a coefficient of variation close to 15% (note that the variability in the input parameters follows σ/µ = 5%). This emphasizes the advantage and the importance of stochastic analysis, in detriment of deterministic analysis, which is key in the robust design and optimization for fatigue life of critical components subjected to time-varying loading.
Additionally, the following linear relationship between plastic CTOD range (δp, in µm) and FCG rate (da/dN, in µm/cycle) from a previous work by Antunes et al. [20] allows predicting the variability in terms of da/dN:
da / dN = 0.5246 × δ p
According to Equation (15), the predicted mean value for da/dN is close to 0.178 µm/cycle, with a coefficient of variation close to 15%. It should be highlighted that the variability analysis results may depend on the assumptions considered in the numerical model, such as the specimen geometry, number of loading cycles between crack propagations, boundary conditions, type of material, among others. This motivates further numerical studies on plastic CTOD range variability. However, linear elastic fracture parameters are commonly adopted to predict da/dN, although fatigue crack propagation is closely linked with plastic deformation occurring at the crack tip. This motivates conducting variability studies focused on the stress intensity factor range ΔK, for comparison purposes with δp.

4. Conclusions

This paper concerns a numerical study on the influence of material parameters and loading variability in the plastic CTOD range of AA7050-T6 on M(T) specimens. First, the most influential input parameters were identified based on the analysis of two types of screening DOE approaches (OFAT and FFD), according to three analysis criteria (Main Effect, Index of Influence and ANOVA). Afterwards, RSM metamodels were established to predict the variability in the δp results in an expeditious way, based on two types of designs (FCCCD and BBD). Additionally, the variability in the FCG rate, da/dN, was also predicted, based on da/dN–δp relationships previously established in a previous work by the authors. The following remarks can be drawn:
  • The type of screening DOE and analysis does not interfere with the identification of the relevant parameters influencing the plastic CTOD range: the parameters E, Y0 and Fmax are consistently shown to be the most influential;
  • Both FCCCD and BBD metamodeling approaches provide similar and accurate predictions of the plastic CTOD range, with RRMSE = 1.7% (FCCCD) and RRMSE = 0.9% (BBD);
  • The predicted variability in the plastic CTOD range results presents right-skewed distributions that follow a coefficient of variation close to 15%.

Author Contributions

Conceptualization, P.A.P.; methodology, P.A.P.; formal analysis, P.A.P.; investigation, P.A.P.; data curation, P.A.P., A.E.M., and M.F.B.; writing—original draft preparation, P.A.P.; writing—review and editing, P.A.P., A.E.M., M.F.B., R.M.B., and F.V.A.; supervision, P.A.P., R.M.B., and F.V.A.; funding acquisition, P.A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research is sponsored by FEDER funds through the program COMPETE – Programa Operacional Factores de Competitividade – and by national funds through FCT – Fundação para a Ciência e a Tecnologia –, under the project UID/EMS/00285/2020. It was also supported by projects: RDFORMING (reference PTDC/EME-EME/31243/2017), co-funded by Portuguese Foundation for Science and Technology, by FEDER, through the programme Portugal-2020 (PT2020), and by POCI, with reference POCI-01-0145-FEDER-031243; EZ-SHEET (reference PTDC/EME-EME/31216/2017), co-funded by Portuguese Foundation for Science and Technology, by FEDER, through the programme Portugal-2020 (PT2020), and by POCI, with reference POCI-01-0145-FEDER-031216. All supports are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paris, P.; Erdogan, F. A critical analysis of crack propagation laws. J. Basic. Eng. 1963, 85, 528–534. [Google Scholar] [CrossRef]
  2. Elber, W. The significance of fatigue crack closure. Damage Toler. Aircr. Struct. 1971, 486, 230–242. [Google Scholar]
  3. Wang, H.; Zhang, W.F.; Sun, F.Q.; Zhang, W. A Comparison Study of Machine Learning Based Algorithms for Fatigue Crack Growth Calculation. Materials 2017, 10, 543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. González, S.B. Phenomenological Approach to Probabilistic Models of Damage Accumulation—Application to the Analysis and Prediction of Fatigue Crack Growth. Ph.D. Thesis, Universidad de Oviedo, Oviedo, Spain, 2019. [Google Scholar]
  5. Castillo, E.; Canteli, A.F.; Siegele, D. Obtaining S–N curves from crack growth curves: An alternative to self-similarity. Int. J. Fract. 2014, 187, 159–172. [Google Scholar] [CrossRef]
  6. Zhu, S.P.; Foletti, S.; Beretta, S. Probabilistic framework for multiaxial LCF assessment under material variability. Int. J. Fatigue 2017, 103, 371–385. [Google Scholar] [CrossRef] [Green Version]
  7. Sankararaman, S.; Ling, Y.; Shantz, C.; Mahadevan, S. Uncertainty quantification in fatigue crack growth prognosis. Int. J. Progn. Health Manag. 2011, 2, 1–14. [Google Scholar]
  8. Wentao, H.; Jingxi, L.; De, X. Probabilistic life assessment on fatigue crack growth in mixed-mode by coupling of Kriging model and finite element analysis. Eng. Fract. Mech. 2015, 139, 56–77. [Google Scholar]
  9. Kim, J.H.; Chau-Dinh, T.; Zi, G.; Lee, W.W.; Kong, J.S. Probabilistic fatigue integrity assessment in multiple crack growth analysis associated with equivalent initial flaw and material variability. Eng. Fract. Mech. 2016, 156, 182–196. [Google Scholar] [CrossRef]
  10. Vélez, A.; Sánchez, A.; Fernández, M.; Muñiz, Z. Multivariate analysis to relate CTOD values with material properties in steel welded joints for the offshore wind power industry. Energies 2019, 12, 4001. [Google Scholar] [CrossRef] [Green Version]
  11. Sankararaman, S.; Ling, L.; Mahadevan, S. Uncertainty quantification and model validation of fatigue crack growth prediction. Eng. Fract. Mech. 2011, 78, 1487–1504. [Google Scholar] [CrossRef]
  12. Ocampo, J.D.; Crosby, N.E.; Millwater, H.R.; Anagnostou, E.L.; Engel, S.J.; Madsen, J.S.; Engel, K.W. Probabilistic damage tolerance for aviation fleets using a kriging surrogate model. In Proceedings of the AIAA SciTech Conference, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar]
  13. Leser, P.; Hochhalter, J.D.; Warner, J.E.; Newman, J.A.; Leser, W.P.; Wawrzynek, P.A.; Yuan, F.-G. Probabilistic fatigue damage prognosis using surrogate models trained via three-dimensional finite element analysis. Struct. Health Monit. 2017, 16, 291–308. [Google Scholar] [CrossRef] [Green Version]
  14. Millwater, H.; Ocampo, J.; Crosby, N. Probabilistic methods for risk assessment of airframe digital twin structures. Eng. Fract. Mech. 2019, 221, 106674. [Google Scholar] [CrossRef]
  15. Karve, P.M.; Guo, Y.; Kapusuzoglu, B.; Mahadevan, S.; Haile, M.A. Digital twin approach for damage-tolerant mission planning under uncertainty. Eng. Fract. Mech. 2020, 225, 106766. [Google Scholar] [CrossRef]
  16. White, P.; Molent, L.; Barter, S. Interpreting fatigue test results using a probabilistic fracture approach. Int. J. Fatigue 2005, 27, 752–767. [Google Scholar] [CrossRef]
  17. Virkler, D.A.; Hillberry, B.M.; Goel, P.K. The Statistical Nature of Fatigue Crack Propagation. J. Eng. Mater. Technol. 1979, 101, 148–153. [Google Scholar] [CrossRef]
  18. Wu, W.F.; Ni, C.C. Statistical aspects of some fatigue crack growth data. Eng. Fract. Mech. 2007, 74, 2952–2963. [Google Scholar] [CrossRef] [Green Version]
  19. Bahloul, A.; Ahmed, A.B.; Bouraoui, C. An engineering predictive approach of fatigue crack growth behavior: The case of the lug-type joint. C R Mec. 2018, 346, 1–12. [Google Scholar] [CrossRef]
  20. Antunes, F.V.; Branco, R.; Prates, P.A.; Borrego, L. Fatigue crack growth modelling based on CTOD for the 7050-T6 alloy. Fatigue Fract. Eng. Mater. Struct. 2017, 40, 1309–1320. [Google Scholar] [CrossRef]
  21. Tvergaard, V. On fatigue crack growth in ductile materials by crack–tip blunting. J. Mech. Phys. Solids 2004, 52, 2149–2166. [Google Scholar] [CrossRef]
  22. Pippan, R.; Grosinger, W. Fatigue crack closure: From LCF to small scale yielding. Int. J. Fatigue 2013, 46, 41–48. [Google Scholar] [CrossRef] [Green Version]
  23. Vasco-Olmo, J.M.; Díaz, F.A.; Antunes, F.A.; James, M.N. Plastic CTOD as fatigue crack growth characterising parameter in 2024-T3 and 7050-T6 aluminium alloys using DIC. Fatigue Fract. Eng. Mater. Struct. 2020, 1–12. [Google Scholar] [CrossRef]
  24. Oliveira, M.C.; Alves, J.L.; Menezes, L.F. Algorithms and strategies for treatment of large deformation frictional contact in the numerical simulation of deep drawing process. Arch. Comput. Meth. Eng. 2008, 15, 113–162. [Google Scholar] [CrossRef] [Green Version]
  25. Chaboche, J.L. A review of some plasticity and viscoplasticity constitutive theories. Int. J. Plast. 2008, 24, 1642–1693. [Google Scholar] [CrossRef]
  26. Prates, P.A.; Adaixo, A.S.; Oliveira, M.C.; Fernandes, J.V. Numerical study on the effect of mechanical properties variability in sheet metal forming processes. Int. J. Adv. Manuf. Technol. 2018, 96, 561–580. [Google Scholar] [CrossRef]
  27. Lasdon, L.S.; Waren, A.D.; Jain, A.; Ratner, M.W. Design and Testing of a Generalized Reduced Gradient Code for Nonlinear Optimization; NTIS National Technical Information Service, U.S. Department of Commerce: Cleveland, OH, USA, 1975. [Google Scholar]
  28. Bonte, M.H.A. Optimisation Strategies for Metal Forming Processes. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 2007. [Google Scholar]
Figure 1. Schematic representation of the M(T) sample: (a) Geometry and dimensions (in mm); (b) Boundary conditions in the frontal view [20].
Figure 1. Schematic representation of the M(T) sample: (a) Geometry and dimensions (in mm); (b) Boundary conditions in the frontal view [20].
Materials 13 01276 g001
Figure 2. Reference numerical simulation results of AA7050-T6: (a) Total crack tip opening displacement (CTOD) vs. applied stress, with schematic representation of CTOD measurement; (b) Plastic CTOD vs. applied stress, with indication of δp.
Figure 2. Reference numerical simulation results of AA7050-T6: (a) Total crack tip opening displacement (CTOD) vs. applied stress, with schematic representation of CTOD measurement; (b) Plastic CTOD vs. applied stress, with indication of δp.
Materials 13 01276 g002
Figure 3. Face-Centered Central Composite Design (FCCCD) and Box-Behnken design (BBD) design points for constructing Response Surface Methodology (RSM) metamodels.
Figure 3. Face-Centered Central Composite Design (FCCCD) and Box-Behnken design (BBD) design points for constructing Response Surface Methodology (RSM) metamodels.
Materials 13 01276 g003
Figure 4. Correlation between the δp values obtained by random numerical simulations and those predicted by the response surfaces constructed with (a) FCCCD and (b) BBD simulations.
Figure 4. Correlation between the δp values obtained by random numerical simulations and those predicted by the response surfaces constructed with (a) FCCCD and (b) BBD simulations.
Materials 13 01276 g004
Figure 5. Plastic CTOD range histograms generated from 100,000 random sets of input parameters using (a) FCCCD and (b) BBD metamodels. The mean δp values and the 2.5th and 97.5th percentiles are also represented.
Figure 5. Plastic CTOD range histograms generated from 100,000 random sets of input parameters using (a) FCCCD and (b) BBD metamodels. The mean δp values and the 2.5th and 97.5th percentiles are also represented.
Materials 13 01276 g005
Table 1. Mean (μ) and standard deviation (SD) values of the studied input parameters; percentiles 2.5th (P2.5) and 97th (P97.5) are also shown.
Table 1. Mean (μ) and standard deviation (SD) values of the studied input parameters; percentiles 2.5th (P2.5) and 97th (P97.5) are also shown.
AA7050-T6E (GPa)νY0 (MPa)CXXSat (MPa)Fmax (N)Fmin (N)
µ71.700.3300420.50228.91198.35385.2919.26
SD3.590.016521.0311.459.9219.260.96
P2.564.670.2977379.29206.48178.91347.5317.37
P97.578.730.3623461.71251.34217.79423.0521.15
Table 2. Numerical simulation results of δp for the one-factor-at-a-time (OFAT) screening design.
Table 2. Numerical simulation results of δp for the one-factor-at-a-time (OFAT) screening design.
SimulationE (GPa)νY0 (MPa)CXXSat (MPa)Fmax (N)Fmin (N)δp (µm)
1µµµµµµµ0.334
2P2.5µµµµµµ0.364
3P97.5µµµµµµ0.304
4µP2.5µµµµµ0.340
5µP97.5µµµµµ0.326
6µµP2.5µµµµ0.379
7µµP97.5µµµµ0.303
8µµµP2.5µµµ0.337
9µµµP97.5µµµ0.326
10µµµµP2.5µµ0.341
11µµµµP97.5µµ0.322
12µµµµµP2.5µ0.264
13µµµµµP97.5µ0.421
14µµµµµµP2.50.335
15µµµµµµP97.50.333
Table 3. Numerical simulation results of δp for the fractional factorial design (FFD) screening design.
Table 3. Numerical simulation results of δp for the fractional factorial design (FFD) screening design.
SimulationE (GPa)νY0 (MPa)CXXSat (MPa)Fmax (N)Fmin (N)δp (µm)
1P2.5P2.5P2.5P2.5P2.5P2.5P2.50.346
2P97.5P2.5P2.5P2.5P97.5P2.5P97.50.273
3P2.5P97.5P2.5P2.5P97.5P97.5P2.50.529
4P97.5P97.5P2.5P2.5P2.5P97.5P97.50.499
5P2.5P2.5P97.5P2.5P97.5P97.5P97.50.418
6P97.5P2.5P97.5P2.5P2.5P97.5P2.50.366
7P2.5P97.5P97.5P2.5P2.5P2.5P97.50.269
8P97.5P97.5P97.5P2.5P97.5P2.5P2.50.220
9P2.5P2.5P2.5P97.5P2.5P97.5P97.50.561
10P97.5P2.5P2.5P97.5P97.5P97.5P2.50.425
11P2.5P97.5P2.5P97.5P97.5P2.5P97.50.492
12P97.5P97.5P2.5P97.5P2.5P2.5P2.50.270
13P2.5P2.5P97.5P97.5P97.5P2.5P2.50.273
14P97.5P2.5P97.5P97.5P2.5P2.5P97.50.232
15P2.5P97.5P97.5P97.5P2.5P97.5P2.50.421
16P97.5P97.5P97.5P97.5P97.5P97.5P97.50.324
Table 4. Main Effect, Index of Influence and ANOVA results obtained from the OFAT simulations (see Table 2). The values in bold indicate that the respective parameters are assumed influential.
Table 4. Main Effect, Index of Influence and ANOVA results obtained from the OFAT simulations (see Table 2). The values in bold indicate that the respective parameters are assumed influential.
OFATE (GPa)νY0 (MPa)CXXSat (MPa)Fmax (N)Fmin (N)
Main Effect0.05980.01330.07520.01090.01920.15610.0016
Index of Influence0.08950.01990.11250.01640.02870.23350.0024
ANOVA p-value0.00100.23090.00030.31440.10200.00000.8778
Table 5. Main Effect, Index of Influence and ANOVA results obtained from the FFD simulations (see Table 3). The values in bold indicate that the respective parameters are assumed influential.
Table 5. Main Effect, Index of Influence and ANOVA results obtained from the FFD simulations (see Table 3). The values in bold indicate that the respective parameters are assumed influential.
FFDE (GPa)νY0 (MPa)CXXSat (MPa)Fmax (N)Fmin (N)
Main Effect0.08740.01620.10870.01010.00140.14600.0270
Index of Influence0.11820.02190.14690.01360.00200.19740.0365
ANOVA p-value0.00870.50390.00310.67510.95160.00070.2815
Table 6. Least-squares solutions of RSM coefficients obtained from the FCCCD and BBD designs.
Table 6. Least-squares solutions of RSM coefficients obtained from the FCCCD and BBD designs.
RSM CoefficientsFace-Centered
Central Composite Design
Box-Behnken
Design
β03.061 × 1001.814 × 10−1
β1−3.124 × 10−2−6.938 × 10−2
β2−6.633 × 10−3−7.472 × 10−4
β3−1.783 × 10−32.630 × 10−3
β43.220 × 10−59.827 × 10−6
β5−7.202 × 10−6−3.067 × 10−5
β6−6.909 × 10−6−1.082 × 10−5
β71.158 × 10−47.009 × 10−5
β87.281 × 10−63.872 × 10−6
β99.765 × 10−68.157 × 10−6
Table 7. Comparison between δp results obtained from the FCCCD simulations with those predicted by the RSM metamodel; the respective R2 value is also shown.
Table 7. Comparison between δp results obtained from the FCCCD simulations with those predicted by the RSM metamodel; the respective R2 value is also shown.
SimulationE (GPa)Y0 (MPa)Fmax (N)δp (µm)δpRSM (µm)
1P2.5P2.5P2.50.3230.334R2 = 0.9927
RRMSE = 2.0%
2P97.5P2.5P2.50.2700.264
3P2.5P97.5P2.50.2710.265
4P97.5P97.5P2.50.2250.233
5P2.5P2.5P97.50.5430.534
6P97.5P2.5P97.50.4520.457
7P2.5P97.5P97.50.4180.423
8P97.5P97.5P97.50.3950.383
9P2.5µµ0.3640.363
10P97.5µµ0.3040.308
11µP2.5µ0.3790.378
12µP97.5µ0.3030.307
13µµP2.50.2640.256
14µµP97.50.4210.431
15µµµ0.3340.330
Table 8. Comparison between δp results obtained from the Box-Behnken simulations with those predicted by the RSM metamodel; the respective R2 value is also shown.
Table 8. Comparison between δp results obtained from the Box-Behnken simulations with those predicted by the RSM metamodel; the respective R2 value is also shown.
SimulationE (GPa)Y0 (MPa)Fmax (N)δp (µm)δpRSM (µm)
1P2.5P2.5µ0.4170.419R2 = 0.9996
RRMSE = 0.4%
2P97.5P2.5µ0.3480.349
3P2.5P97.5µ0.3350.334
4P97.5P97.5µ0.2770.276
5P2.5µP2.50.2930.292
6P97.5µP2.50.2440.244
7P2.5µP97.50.4710.471
8P97.5µP97.50.3890.391
9µP2.5P2.50.2940.293
10µP97.5P2.50.2460.248
11µP2.5P97.50.4930.490
12µP97.5P97.50.3780.378
13µµµ0.3340.334
Table 9. Descriptive statistics of the δp histograms shown in Figure 5.
Table 9. Descriptive statistics of the δp histograms shown in Figure 5.
δpFCCCD-Based Model BBD-Based Model
Mean value, µ0.338 µm0.340 µm
Standard deviation, SD0.050 µm0.049 µm
Coefficient of variation, CV14.9%14.5%
2.5th percentile, P2.50.252 µm0.257 µm
97.5th percentile, P97.50.448 µm0.448 µm

Share and Cite

MDPI and ACS Style

André Prates, P.; Eusébio Marques, A.; Frias Borges, M.; Madeira Branco, R.; Antunes, F.V. Numerical Study on the Variability of Plastic CTOD. Materials 2020, 13, 1276. https://doi.org/10.3390/ma13061276

AMA Style

André Prates P, Eusébio Marques A, Frias Borges M, Madeira Branco R, Antunes FV. Numerical Study on the Variability of Plastic CTOD. Materials. 2020; 13(6):1276. https://doi.org/10.3390/ma13061276

Chicago/Turabian Style

André Prates, Pedro, Armando Eusébio Marques, Micael Frias Borges, Ricardo Madeira Branco, and Fernando Ventura Antunes. 2020. "Numerical Study on the Variability of Plastic CTOD" Materials 13, no. 6: 1276. https://doi.org/10.3390/ma13061276

APA Style

André Prates, P., Eusébio Marques, A., Frias Borges, M., Madeira Branco, R., & Antunes, F. V. (2020). Numerical Study on the Variability of Plastic CTOD. Materials, 13(6), 1276. https://doi.org/10.3390/ma13061276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop