The paper examines the influence of uncertainty parameters on the wave propagation responses at high frequencies for a damaged beam structure in the structural reliability context. The reliability analyses were performed using the perturbation method, First-Order Reliability Method (FORM), and response surface method (RSM) which were compared with Monte Carlo simulation (MCS) under the spectral element method environment. The simulated results were performed to investigate the effects of material property and geometric uncertainties on the response at high frequency modes, such as the computational efficiency of reliability methods. For the first time, the spectral element method is used in the context of reliability analysis at medium and high frequency bands applied to damage detection. It has shown the effects of parameters uncertainty on the dynamic beam response due on an impulsive load and the robustness of each method. Numerical examples in a bending vibrating beam with random parameters are performed to verify the computational efficiency of the present study.

1. Introduction

At medium and high frequency bands approaches like the finite element method (FEM) create large numerical models, which can require very high computational times. Statistical Energy Analysis (SEA) produces small models but without spacial variation in subsystems [1]. An alternative to these problems is the spectral element method (SEM) [25]. It consists in the analytical solution of the displacement wave equation in the frequency domain written in the form of a finite element. Then, a spectral element is equivalent to an infinite number of finite elements. This characteristic and the spectral domain make SEM more suitable to solve high frequency band problems. Its wave propagation formulation is also more adequate to solve damage detection problems. In the last decade, researches about damage detection concentrate on methods that use elastic wave propagation at medium and high frequencies to detect and quantify structural damage [69]. They use the evidence that material discontinuities, such as a crack, generate changes in the elastic waves propagating into the structure [10]. Some particular advantages of elastic wave-based damage detection methods include their capacity to propagate over significant distances and their high sensitivity to discontinuities near the wave propagation path. The presence of a structural damage introduces a local flexibility change that modifies its vibration response [1113]. Therefore, vibration energy can be used to investigate the damaged condition of a structure. The treatment of uncertainties using spectral element method is recent [14, 15], and very few were made related with detection and assessment of the damage. Recently, the authors [1618] and some other researchers [1921] have presented works in damage detection using wave propagation in the context of uncertainty quantification and stochastic SEM model.

Engineering analysis consists in verifying and certifying that the system complies with many performance criteria, safety, and durability under different kinds of solicitation. It is well known that there are many uncertainty sources in external loads and structural parameters. Consequently, the necessity to incorporate this information in the system design is increasing. With the aim to include uncertainty in system design, reliability methods have been developed in the last few decades. First-order second-moment (FOSM), First-Order Reliability Method (FORM), and Second-Order Reliability Method (SORM) [2224] are methods based on Hasofer and Lind transformation [25], which include uncertainty in the reliability analyses. Frequently, in the reliability analysis of complex structures, the limit state function cannot be expressed in a closed form. Typically, it needs to be evaluated implicitly through an approximated solution, such as a finite element method. The reliability analysis can also be carried out using Monte Carlo (MC) simulation [26]. However, MC requires a large number of realizations to converge, which could be an expensive technique and for large structures it may be infeasible. In order to reduce the computation time the response surface method (RSM) has been applied [2732]. Usually, in this case, the actual limit state function is approximated by a polynomial function. Melchers and Choi [33, 34] present the techniques for reliability analysis of engineering structures using probability theory. In the context of damage detection and crack propagation, reliability methods were employed using approximated techniques such as boundary element method (BEM) [35], finite element method [22, 23, 33, 36], and finite element model updating with wave propagation [37]. In this paper a new contribution to the structural damage detection of a nonpropagating crack in the context of reliability methods using the spectral element method is presented.

In this study the structural damage detection problem using the the spectral element method together with reliability algorithms is presented. The damaged beam spectral element [6] is extended to include uncertainty in material property and geometric parameters. A straightforward procedure to estimate the parameter randomness is the Monte Carlo simulation [26]. It consists in obtaining a large number of deterministic analyses with different realizations of random variables. For complex structures with implicit performance function the evaluation of each realization is computationally expensive and the reliability analysis could be impracticable. Therefore, MC simulation becomes computationally infeasible to estimate the failure probability for large structures. Approximated methods like First-Order Reliability Method (FORM) and response surface method (RSM) [28] provide more appropriated tools to estimate the structural failure in these cases.

To consider probabilistic damage detection problems and random crack parameters, it is required to couple the reliability procedures with the spectral element model. This coupling can be performed using either direct coupling or response surface method. In the direct coupling the limit state function derivatives are calculated directly based on the numerical response of the SEM model [36]. In this approach, the limit state function remains implicit and is defined by the SEM model responses. The basic procedure consists in directly coupling the reliability model, FORM, with the numerical beam model based on SEM. In the reliability analysis the RSM is used to approximate the beam response at the vicinity of the most probable failure point in terms of the random variables (material properties, geometries, etc.). The beam response is estimated by a local polynomial approximations in an iterative process, as the search of the design point is evaluated [35]. Numerical tests are presented in order to compare the coupling of the numerical model based on SEM with reliability methods and the computational efficiency of each approach.

2. Spectral Element Method

The spectral element method is similar in style to the finite element method. However, there are two important differences between these methods. The first one is that SEM is a wave propagation formulation written in the frequency domain. The second one is that the element interpolation function is the exact analytical solution of the differential equation. Based on these characteristics the number of elements required for a spectral model will coincide with the number of discontinuities in the structure (Figure 1).

2.1. Beam Spectral Element

For this model it is assumed that the beam is a slender structure subjected to transversal loads and bending moments. It is also considered that transversal and rotational displacement are small. By neglecting shear deformations, the differential equation of movement in its spectral form can be written as follows [3]: with the homogeneous solution given by where is the transversal displacement in the frequency domain and is the beam length. The wave numbers, , , and , are given by where is the circular frequency, is Young’s modulus, is the cross section area, is the density, is the inertia moment, and . By using a complex Young’s modulus, , an internal structural damping is introduced where is the hysteretic structural damping factor. Figure 2 illustrates a two-node undamaged beam spectral element model with two degrees of freedom by nodes ( and ) and nodal loads ( and ).

By applying the displacement and rotation boundary conditions to nodes 1 and 2 (equation (2)) the coefficients, , , , and , can be obtained. By substituting these coefficients in (2), the expression to calculate the displacements and the rotations for any arbitrary point of the beam element can be written as where are interpolation functions, which are omitted here for brevity [3]. By using the relationship between nodal loadings and displacements and applying the boundary conditions to the beam spectral element, the following matrix equation is obtained: where is the undamaged beam spectral element dynamic stiffness matrix, which is symmetrical and usually complex and can be expressed as

The throw-off beam spectral element can be thought of as a particular case of the undamaged beam spectral element where node 2 on right-hand side goes to infinity (Figure 3). In the context of probabilistic analysis this element will be considered as a deterministic element. Following a similar procedure, the throw-off beam spectral element dynamic stiffness matrix can be obtained as

The damaged beam spectral element is a beam including a transverse, open, and nonpropagating crack [6] as illustrated in Figure 4. The crack is modelled by a dimensionless local flexibility, , which is calculated based on Castigliano’s theorem and the laws of fracture mechanics [38].

The solution for (2) must be written in two parts: where and are the vertical displacement on the left- and right-hand side of the crack, respectively. The coefficients , , , , , , , and are determined by the boundary conditions. Displacements, nodal loads, and the element dynamic stiffness matrix are obtained in a similar way as the undamaged beam spectral element. Then, the damaged beam spectral element dynamic stiffness matrix can be written as where , , , , , and .

2.2. Flexibility at the Crack Position

The crack flexibility coefficient () is calculated using Castigliano’s theorem, where the flexibility at the crack location for the one-dimensional beam spectral element can be obtained as where denotes the elastic strain energy due to the crack and is the nodal force on the element.

Figure 5 shows a detail of the cross section area at the crack position. By considering that only crack mode is present in the beam element, the elastic strain energy can be expressed as where is Poisson’s ratio, is the damaged cross section area, and is the stress intensity factor corresponding to the crack mode , which can be written by where is the cross section base, is the cross section height, is the crack depth, is the bending moment at crack position, and is the fracture mechanics correction function given by It can be shown that crack flexibility can be expressed as where is the dimensionless crack depth. The dimensionless local flexibility is given by

A spectral element for undamaged and damaged beam was present in this section. These models will be used to simulate the response of the structural system. Uncertainty parameters are considered in these formulations with the exception of throw-off element. It is used here only to throw off the energy of the system and has no significant influence on the dynamic system response.

3. Coupling Reliability Analysis and SEM Modelling

The aim of reliability analysis is to estimate the failure probability in view of a certain failure scenario [33, 34]. The performance or safety of the structure can be described by limit state functions , where the parameters which are significant sources of uncertainty in the system are grouped in vector . The limit state functions of interest can be defined as , although in real situation a random behaviour can be observed. In this case the limit state function is separate in safe domain, where , and in failure domain, where . However, the limit state function can be calculated at some desired points by running the SEM analysis. The failure probability is calculated by the integral of the limit state function as

Reliability analysis is based on probabilistic algorithms to calculate the consequence of a random input on the system response. Uncertainty quantification is performed by computing approximate response functions based on random variable probability distributions. These methods are more able to compute events with low probability than sampling based approaches, since the sampling size required to solve can be unaffordable. In order to evaluate probabilistic damage detection problems, it is required to couple the reliability procedures with the spectral element model. This coupling will be performed using direct coupling and response surface method. Different treatments of reliability procedures coupling with SEM are used to perform the reliability analysis of structural damage: First-Order Reliability Method (FORM), First-Order Reliability Method with limit state function obtained by response surface first-order approximation (FORM_RSM1), Monte Carlo using response surface first-order (MC_RSM1) and second-order (MC_RSM2) approximations, and Monte Carlo (MC) simulation as reference.

3.1. Direct Coupling

The direct coupling or direct method consists in directly coupling the reliability model, FORM or MC, with the numerical model of SEM. As mentioned, the limit state function defines the interface between safe and failure domains. In a context of structural damaged detection analysis, the limit state function is defined in function of the rms (root mean square) value of deterministic system response, , and random system response, . Both system responses derived from SEM numerical model in time domain.

When the system limit state equation or performance function is known, the FORM method can be used for reliability analysis. The limit state equation can be a linear function of correlated or uncorrelated normal variables or a nonlinear function represented by its first-order approximation. If there are nonnormal variables, the method requires to work with equivalent normal variables [23]. The development of the FORM method can be traced to second-moment methods: first-order second-moment (FOSM) and advanced first-order second-moment (AFOSM). They consider the first and the second moments of random variables, that is, the mean and standard deviation. In the structure field, a failure could be defined as a structural parameter, whose value is higher than a given critical value. In this case we have a simple limit state function defined as where is the critical value of the variable to the random variables . For all cases failure domain is defined as and safe domain as . Then, the failure probability will be defined as where is the joint probability density function of the random variables and is the region of failure . In the FORM method it is approximated by where is the cumulative distribution function of the standard normal variable and is the distance from the origin to the approximated hypersurface called reliability index [25, 39]. For the FORM algorithm the limit state function needs to be solved to find the new design point, although this may be difficult in the case of complex nonlinear or implicit limit state functions. An algorithm suggested by Rackwitz and Fiessler [39] can be applied to overcome this problem, which uses the derivatives to find the next iteration point.

Another direct coupling approach presented in this paper and also used as a reference is the Monte Carlo simulation. A widely used application of Monte Carlo simulation is to determine the value of a multidimensional definite integral. Let us consider that the following integral needs to be evaluated: By assuming as a probability density function we can estimate as where the values have been sampled with density and is the indicator function given by The MC approach uses the so-called Simple Random Sampling, which is quite inefficient in the majority of cases. From (21) it is obvious that the accuracy of the estimation will depend on the number of samples. The failure probability would lead to the true value as the number of samples goes to infinity. Then the sample size to accurately estimate the failure probability is a matter of concern. There are several ways to evaluate the accuracy of (21). One way is to study the error associated with the sample size by approximating the binomial distribution with the normal distribution and estimating the 95% confidence interval of the failure probability [40]. It can be shown that the sampling size with a confidence of 95% can be obtained with In the reliability context, the failure probability is computed based on (22). For a certain number of samples, when the state limit equation is lower then zero () a failure is considered, and the reliability index is calculated by the first and second moment of the number of samples that was considered.

3.2. Response Surface Method (RSM)

The response surface method (RSM) is a set of statistical techniques designed to find the best response value, considering the uncertainty or variations in the values of input variables [41]. This method can be employed to estimate the structural statistic response or the implicit limit state functions, which is represented by the mean and coefficient of variation of the response. In RSM, the true limit state function, , is approximated by a simple and explicit mathematical expression, which is a th order polynomial with undetermined coefficients. The value of the limit state function is evaluated at a number of samples of , to determine the unknown coefficients such that the error of approximation at the samples of is minimized. The selection form of the approximated limit state function should be based on the shape and the nonlinearity of the true limit state function. Since is usually unknown, there has been a tendency to develop a generic form for a response surface which can be applied across a wide range of structural reliability problems. The most common forms are the first- or second-order polynomial. However, to capture the nonlinearity of the limit state function more precisely, cross terms are sometimes included into the polynomial, which can be written in a more general form as where () are the th random variables, is the total number of random variables, , , , and are the unknown coefficients to be determined, and is an approximate representation of the limit state function . The values of the coefficients can be determined via regression analysis using a set of sample points from the limit state function, . The number of sample points must be larger or at least equal to the number of coefficients. Among various sampling methods, a common Design of Experiment (DOE) approach is evaluated for at combinations of and , where and are the mean and standard deviation of and is an arbitrary factor. In this case, the number of sample points is just sufficient for the determination of the coefficients, by this means minimizing the number of sample points and the number of evaluations of the limit state function. The RSM can approximate the limit state function accurately for roughly linear and quadratic limit states. However, when the shape of the limit state function is not close to linear or quadratic, parameter plays an important role in the accuracy of the second-order RSM approximation [31]. In the solution of structural reliability problems RSM is used to approximate the structural response at the surroundings of the most probable failure point, in terms of input variables [42]. From that approximated structural response the failure probability can be evaluated using reliability methods. In this work Monte Carlo simulation and FORM are used to estimate the failure probability.

To compute the failure probability using RSM coupled with FORM the following procedure needs to be done. Different sets of points related with are chosen according to the DOE. The selected points are used as input to calculate the structural response by SEM. Then the structural responses are used to build the response surface approximated by a polynomial expression. The next step is to compute the limit state function in a standard normalized space (24) and to calculate the reliability index (), which is estimated through an optimization procedure. Finally, the failure probability is estimated according to FORM approximation from equations (18) and (19).

As mentioned, the limit state function in RSM approach is constructed using a polynomial approximation fitted according to the structural responses at some desired points. In this application the response surface is used to build the approximate state limit function based on the structural response from SEM. Since the RSM can only provide a local approximation, an iterative process to map the safe and failure domains needs to be done. This procedure is based on the construction of successive surfaces, each of them valid only in the surroundings of failure point. The failure probability is computed based on the amount of sample and when the approximate state limit equation is lower then zero () a failure is considered. The reliability index is calculated for the first and second moment and failure probability obtained. In order to evaluate the efficiency of these methods under the spectral element environment a numerical test will be presented in the next section.

4. Numerical Analysis

The first test is performed with deterministic parameters in order to evaluate the computational code to calculate the dynamic responses for the undamaged and damaged beam spectral element. The second test will verify the efficiency of coupling reliability analysis with spectral element modelling under parameter uncertainty environment.

4.1. Deterministic SEM Dynamic Response

For the numerical tests of SEM dynamic response two structures are analysed: one composed of an undamaged beam spectral element connected with a throw-off beam spectral element and the other composed of a damaged beam spectral element connected with a throw-off beam spectral element. Both structures are illustrated in Figure 6. Assuming a free-free boundary condition the structures are excited at node 2 on the right-hand side by a tone-burst impulsive force. The beam structures are made of steel ( = 210 GPa, = 7850 kg/m3, and = 0.01) with element length = 4 m and rectangular cross section ( = 0.01 m and = 0.03 m). In regard to the crack, its position is fixed at = 2 m and the crack depth will be varying from a small value to a value close to the failure limit ( = 0.02–0.28).

The tone-burst impulsive force is used to improve the wave propagation visualization and to localize the crack position in the beam. Tone-burst pulse is generated using a sine signal windowed over 5 cycles with a triangular window. The sine is generated at the frequency of 15.0 kHz with 100 N of amplitude. Figure 7 shows the tone-burst pulse in time and frequency domain.

For the simulated test the excitation force is applied at node 2 of the beam element, and the acceleration response is obtained at the same node. Figure 8 shows the time acceleration response for the undamaged (a) structure. It presents two pulses, one at the excitation moment (first pulse) and the other at the reflection moment (second pulse). The pulse wave travels through the structure until it reaches node 1 where it is reflected back to node 2 and continues through the throw-off element to the infinity. Throw-off element works as an anechoic termination dissipating the remaining energy into the signal. Due to structural damping () the amplitude of vertical acceleration decreases through the structure length. Figure 8 shows the time acceleration response for the damaged (b) structure. The acceleration response shows two additional pulses as compared to the undamaged case, one between the excitation and reflection moments and the other after the reflection moment. These additional pulses characterize the waves partially reflected and transmitted due to any structural discontinuity, which could be a crack. These results show that the SEM model is able to reproduce the wave propagation behaviour in the undamaged and damaged beam structure and to localize the crack.

4.2. Cracked Beam Reliability Analysis

The reliability analysis will be performed for the damaged beam structure. The limit state function used to verify the structural reliability evaluates the rms value of the acceleration time response vector. It is calculated at node 2 of the beam structure (), which must be greater than a critical value (). By establishing the limit value to the nondimensional crack depth as (30% of beam cross section height) the critical deterministic rms value of the acceleration time response at node 2 of the beam structure can be obtained. From (17) the limit state function can be written as The uncertainity parameters considered are the nondimensional crack depth with mean of [0.02 : 0.08 : 0.28] and coefficient of variation [0.01 : 0.05 : 0.31], Young’s modulus with mean of  GPa and , crack location with  m and coefficient of variation of , and cross section height with  m and . The random parameters are considered uncorrelated random variables with Log-normal distribution. The Log-normal distribution was assumed because the chosen random parameters cannot assume negative values. All other structural parameters are assumed as deterministic.

In order to see the effect of random parameters in the structural beam response as compared with the deterministic response a simulated test was made. Considering the specified values for the means and ’s for the structural random parameters, 100 samples of acceleration responses were calculated and compared with the deterministic response. The results are shown in Figure 9. As demonstrated, the dispersion in the parameters has a great influence in the acceleration responses. This kind of behaviour is normally seen in experimental tests. Based on this idea the random parameter is incorporated in the present study, and methods to treat this kind of randomness are used in the context of reliability analysis. By taking the crack depth with and and keeping all other random and deterministic variables as defined, a Monte Carlo simulation is performed to evaluate the convergence of the failure probability. Figure 10 shows that failure probability convergence can be achieved with a sample size >500,000. By using (23) with the converged failure probability, 0.038, and a 5% error 40,000 samples are obtained. Therefore, from now on a sampling size of 50,000 samples will be taken for all Monte Carlo simulations along the paper.

To evaluate performance and computation time to estimate failure probability using the presented methods, the following simulations are carried out. The reference for the failure probability is calculated by Monte Carlo simulation with SEM structural damaged beam model and is noted as MC. Failure probability computed by RSM first-order approximation uses full factorial DOE. It considers four random variables () with the levels , , and evaluates the performance function for all possible combinations ( design points). By using (24) the coefficients for the first-order polynomial are estimated using regression analysis. For the RSM second-order approximation the Box-Behnken design (BBD) of experiments is used. The BBD is an independent quadratic design in that it does not contain an embedded factorial or fractional factorial design. In this design the treatment combinations are at the edges midpoints of the process space and at the centre. These designs are rotated or near rotated and require three levels of each variable. For this case we consider the four random variables () with the levels and . The limit state function is calculated for the second-order model without cross terms; then coefficients and design points. From (24) the coefficients for the second-order polynomial are estimated using regression analysis. When the failure probability is computed with Monte Carlo simulation with first- and second-order RSM approximation the model is referred to as MC_RSM1 and MC_RSM2, respectively. The FORM with RSM first-order approximation works out the failure probability using the derivative of limit state function computed by finite difference method. The limit state function is approximated by a polynomial and is referred to as FORM_RMS1. The direct coupling which is the FORM coupled directly with SEM is called FORM. In order to evaluate the sensitivity of damaged beam spectral element response related with the crack depth, all simulations are performed, varying the crack depth from a small value to a value close to the allowed failure limit (). Moreover, for each case, changes in crack depth coefficient of variation are made with . Figure 11(a) shows the failure probability calculated by FORM coupled directly with SEM model. The results show that for crack depth values in the range of the probabilities of failure are very small, even for the highest values of . As the crack depth values increase to the failure probability starts to raise as soon as the value starts to increase.

For the failure probability rises in a very fast rate as the value increases. Figure 11(b) shows the failure probability calculated by FORM_RSM1. The results show that for crack depth values in the range of the probabilities of failure are very small as observed in the direct method. As the crack depth values increase to the failure probability starts to raise as soon as the value starts to increase. However, in the limit state function approximated by a surface response function the failure probability can take a smooth increasing with the highest COVs than that obtained with direct FORM.

For the failure probability rises in a very fast rate as the value increases. Figure 11(c) shows the failure probability calculated by MC_RSM1. The failure probability results present a similar behaviour as that of FORM_RSM1.

Figure 11(d) shows the failure probability calculated by MC_RSM2. For this case the general behaviour seems to be similar to the previous cases in that increases as and increase. Nevertheless, the failure probability starts to raise earlier. When the crack depth value and increases starts to raise, but in a slower rate than the previous cases.

Figure 11(e) shows the failure probability calculated by MC, which presents a similar behaviour as MC_RSM2. The general behaviour of these results seems to be in agreement with the physic of the problem, since the probabilities of failure are very small for for all ’s values and become higher for at high values of .

In order to make the result analysis more clear, two additional plots comparing the results of all methods together are presented. In the first the low probabilities of failure (Figure 12(a)) are grouped and in the second the high probabilities of failure (Figure 12(b)) are grouped. Figure 12(a) shows that low probabilities of failure obtained by MC_RSM1, FORM_RSM1, and FORM present a good agreement for all values of crack depth in this range ( = 0.02–0.18). By comparing direct Monte Carlo (MC) simulation with MC_RSM1, FORM_RSM1, and FORM a significant difference appears, which means that the nonlinearities in the damaged beam spectral element do not allow a linear approximation for the state limit function. This can be confirmed by the good agreement among Monte Carlo (MC) simulation, FORM, and Monte Carlo with a second-order response surface (MC_RSM2). Figure 12(b) shows that for high probabilities of failure a similar behaviour is observed.

Monte Carlo with first- and second-order response surface (MC_RSM1 and MC_RSM2) presents very close computational times that are much lower than Monte Carlo simulation (MCS) (Table 1). FORM shows a computational time of less than 10% of MC_RSM1 and MC_RSM2. FORM with linear response surface (FORM_RSM1) presents the lowest computation time. This reduced computation time of FORM_RSM1 as compared with FORM comes from the limit state equation that must be calculated by the last approach. Moreover, stable numerical models (as SEM) lead to smooth and faster convergence of FORM, which makes this approach more efficient than other reliability algorithms as Monte Carlo simulation and response surface method. All runs were performed in a Desktop-PC using the processor Intel Core i7 CPU 930 @ 2.80 GHz, 12.0 Gb RAM memory, and Windows 7, 64 bits as the operating system.

5. Conclusion

The wave propagation responses in a damaged beam structure including uncertainties in material and geometric parameters are analysed. A reliability analysis is made to quantify the failure probability of a damaged structure modelled with beam spectral element at high frequencies. The limit state function evaluates the rms value of the acceleration time response which must be above the limit value. A crack depth of cross section height defines the allowed deterministic rms value of the acceleration time response. By varying crack depth () and its coefficient of variation () the reliability analysis is performed using linear response surface method (RSM1) solved by Monte Carlo simulation (MC_RSM1) and First-Order Reliability Method (FORM_RSM1). Also, a quadratic response surface method (RSM2) is solved by Monte Carlo simulation (MC_RSM2) and all methods are compared with Monte Carlo (MC) simulation. In general, all methods present results in agreement with the physics of the problem; that is, the probabilities of failure are very small for small crack depth and become higher as increases. The agreement between MC_RSM1 and FORM_RSM1 confirms other literature results that a linear RSM is equivalent to computing the FORM solution directly. The disagreement between Monte Carlo (MC) and first-order RSM models (MC_RSM1 and FORM_RSM1) reveals that the degree of nonlinearity in the damaged beam spectral element does not allow a linear approximation for the state limit function. This can be confirmed by the good agreement between Monte Carlo (MC) results and Monte Carlo with a second-order Response Surface (MC_RSM2). As expected computational time results confirm that Response Surface approach reduces meaningfully the computational time as compared with direct Monte Carlo simulation. In this case Monte Carlo with second-order response surface (MC_RSM2) represents a reduction of 12,159.3% in computational time as compared with Monte Carlo (MC).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.