#### Abstract

This paper considers structural response analysis when structural uncertainty parameters distribution cannot be specified precisely due to lack of information and there are complex dependencies in the variables. Uncertainties in parameter are quantified by probability boxes (-boxes) and dependence among uncertain parameters is modeled by copula. To calculate uncertainty structural response, a sampling-based method is proposed. In this method, a sampling strategy is used to sample random intervals from dependent -boxes according to the copula theory and the metamodel-based optimization method is applied to solve a range of structural interval response problems. Two types of errors are presented to evaluate the error of different -boxes. Four numerical examples are investigated to demonstrate the effectiveness of the present method.

#### 1. Introduction

Traditional quantification of the uncertainties existing in a system and the corresponding calculation of uncertainty propagation are generally based on the probability model, in which random distributions are used to describe the uncertainty. Unfortunately, for practical engineering problems, sufficient experimental samples are not always available or sometimes very expensive to obtain. Thus alternative imprecise probabilities have been proposed, including the interval theory [1], the evidence theory [2], and probability boxes (-boxes) [3]. The focus of this paper is on -boxes, which are more expressive generalization of both traditional probability distributions and interval representations. The -box incorporates facilities from probability theory for modeling correlations and dependencies and from interval analysis expressing ignorance by expressing interval bounds on the cumulative distribution function (CDF) for a random variable [4]. On the other hand, the -box has a clear behavioral interpretation and therefore is easy acceptable for practicing engineers [5].

In the field of uncertainty analysis based on -boxes, a lot of exploratory work has been reported in theoretical and application aspects. For example, Williamson and Downs [6] proposed an algorithm to compute binary arithmetic operations (addition, subtraction, multiplication, and division) on pairs of independent discrete -boxes for risk assessment. Karanki et al. [7] used the two-phase Monte Carlo simulation to calculate the multiplication of two -boxes for the safety assessment. Bruns [8] compared three uncertainty propagation methods to find the appropriate computational methods for propagating -boxes through the black box engineering models. Zhang et al. [9–11] combined the interval finite element analysis (FEA) and the interval sampling method and proposed an interval Monte Carlo (IMC) method for structural analysis with -boxes which are defined by distributions with interval parameters. Bai et al. [12] developed a numerical method to compute the linear elastic static and dynamic response of structures with epistemic uncertainty represented by evidence variables, which integrated the moment concept and finite element method and also suited to -boxes variables. Ghosh and Olewnik [13] improved the performance of the two-phase Monte Carlo simulation by replacing the outer loop by optimization algorithms and replacing the inner loop by the sparse grid numerical integration.

Among the above-mentioned methods, they focused on discussing uncertainty propagation of independent -boxes. However, the complex dependencies are common in physical systems and can have profound impacts on the numerical results of calculations [14–16]. Ferson et al. [17] illustrated several methods which are used to model the dependence among variables that are cannot be generalized easily for use with -boxes and proposed that copulas could represent easily the dependence in -boxes. Based on the analysis, they presented an approach known as the dependency bounds convolution (DBC) method. The DBC method could calculate the propagation of dependent -boxes whose dependence is expressed by given parametric copula. Despite the fact that DBC is useful for a sequence of basic arithmetic operations or elementary function with a small number of uncertain variables, it must overcome some obstacles for application in engineering design. Among these problems, the low efficiency difficulty seems to be the severest one. The Cartesian product method is applied in the DBC, which can impose significant computational burden.

In this paper, a computational method is proposed to propagate dependent -boxes through structural model, where dependence in -boxes is described by copula. By the copula sampling technique, the input epistemic uncertainties represented by -boxes are transformed into a range of intervals. The metamodel-based optimization is employed as the solver to calculate the structural response with input intervals. A mathematical function has been tested to show the procedure and the efficiency of the proposed method and then three structural numerical examples are investigated by the present method.

#### 2. Probability Boxes Theory under Dependence

##### 2.1. Probability Boxes

A -box is a class of distribution functions delimited by an upper and a lower bound which collectively represent the epistemic uncertainty about the distribution function of a random variable. Let denote the space of distribution functions on the real numbers . A -box will be defined aswhere , , and are the lower and upper bounds of a -box, which is also written as .

There are several ways to construct -boxes, depending on the type of information available. In this paper, -boxes are constructed based on 95% confidence intervals on the parameters of a known distribution type. The distribution types are determined on the basis of the theoretical knowledge or previous experience.

##### 2.2. Describing Dependence in Probability Boxes with Copula

Now, consider the dependence in -boxes which can be described by copulas. A copula is a multivariate probability distribution for which the marginal probability of each variable is uniformly distributed. By Sklar’s theorem [18], copulas are simply the dependence functions that knit together marginal distributions to form their joint distribution. Consider a random vector with CDFs which are represented by -boxes in this paper, the multivariate cumulative distribution function can be written aswhere is a multivariate copula.

##### 2.3. Dependency Bounds Convolution Method

Ferson et al. [17] firstly used dependency bounds convolution method to calculate binary arithmetic operations on pairs of dependent discrete Dempster-Shafter structures. This method mainly contains discretization of the -boxes and Cartesian product. Consider that , are two dependent -boxes connected by copula function , and is the response quantity. The -boxes and are discretized into focal elements, which are represented by and , respectively. By Yager’s Cartesian product, the response quantity in the th focal element can be defined by an interval , where . The probability mass associated with the th focal element can be calculated bywhere and are the cumulative masses, and . The result is written as focal elements which can be transferred into a -box bywhere and are upper bound and lower bound, respectively.

#### 3. Structural Response Analysis

Using -boxes to quantify the structural parameters, the structural response will be also a -box variable. In this section, we will develop an efficient response analysis method for structures by combining -box theory and the optimization method.

##### 3.1. Sampling Random Intervals from the Dependent -Boxes

By virtue of Sklar’s theorem [18], pseudorandom samples can be generated from (2). That is, given a procedure to generate a sample from the copula , the required sample can be constructed aswhere are quasi-inverse function of CDF and are a sample from a random vector . So the random intervals in accordance with the dependent probability boxes can be generated by following two steps.

The first is to generate a sample from copula. There are lots of algorithms which can generate a sample from a specified copula function [19]. As conditional method may be applied for every chosen copula, it is used as the sample method and illustrated in this section. Just to explain this method in a simple way, let us assume a bivariate copula in which all of its parameters are known and pairs need to be generated from bivariate copula function . Two independent uniform (0, 1) variates and are generated, and then , where denotes a quasi-inverse of which equals the partial derivative of the copula .

The other step is to compute dependent random intervals from -boxes. Suppose that are a sample from the copula and , are their CDFs. For each and corresponding to , a random interval is generated:where is random interval from . The superscript represents the interval distribution parameters and the superscripts and , respectively, represent the lower and upper bounds of interval. and denote the quasi-inverse function of upper and lower bounds of a -box. Thus using will simulate a range of dependent intervals .

Zhang et al. [9] presented an IMC method for sampling intervals from the independent -boxes. When random variables contain both dependent and independent -boxes, the present sampling method and the IMC method can be used to perform random sampling for dependent and independent -boxes, respectively.

##### 3.2. Metamodel-Based Interval Analysis

Through above section sampled treatments, a range of interval vectors is acquired from -boxes. Consider the vector of uncertain parameters defined to be contained within an interval vector (or hypercube) . Interval vectors sampling times from -boxes are represented as (). For a practical engineering problem, the input-output relationship between these parameters and the output quantity of interest is represented by the function applied on these parameters. The interval analysis for the function is numerically equivalent to solving the following equation:where contains all which are obtained from applying the function on all possible vectors within the interval vector . In many cases, an exact description of is extremely difficult to find. Therefore, an interval is used to approximate the exact solution . For the black box model in the practical engineer problem, solution of will make the computational cost extremely expensive [20]. Besides, a great number of repeating interval analyses will be conducted, which will generate unaffordable computational burden. So metamodeling techniques are used to approximate for improving the computational cost.

To improve the accuracy of metamodel, design of experiments (DOE) should be applied. An experimental design domain can be acquired by truncating the -box at a probability of ( is level of significance), which are represented as . The vector form is . In this paper, the Latin hypercube design (LHD) method is adopted, which is a space-filling design with constrainedly stratified sampling method. The sample of this method does not increase exponentially with the number of variables. Radial basis function (RBF) [21] is used to approximate the function , which can be expressed as where is the sampling point selected through LHS, is the number of sampling points, is mass coefficient, and is the basis function. The common basis functions include thin-plate spline, Gaussian, and multiquadrics. The metamodel-based solution of can be transferred into the following optimization problems: The gradient-based optimization techniques such as sequential quadratic programming can be used to solve this problem, and the initial point can be selected on the bound point of intervals.

##### 3.3. Computation of Output Quantity

Repeating times interval analysis, we will acquire intervals , which can be transformed into a -box by [9] with and , where is the indicator function. Summarizing the above procedure, the present method can be summarized as follows.

*Step 1. *Collect all -boxes and sample interval vectors from by the method in Section 3.1.

*Step 2. *Conduct the design of experiments with LHS in the experimental design domain .

*Step 3. *Construct a RBF response surface for approximating (7).

*Step 4. *Conduct optimization analysis over each and acquire interval vector .

*Step 5. *Generate the -box of the output quantity by (10).

#### 4. Error Analysis

The absolute error and the relative error are popular to measure the accuracy between the experimental data and the true value, which is extended to analysis errors of different -boxes. In this paper, the absolute errors contain the area error and the boundary error. Suppose that the two -boxes and ; then the left area error and the right area error between two -boxes are, respectively, calculated by

The area error between two -boxes is defined as the sum of the left area error and the right area error, which is written as . As the area error cannot reflect the discrepancy of boundary between two -boxes, the discrepancy of the boundary is used as the boundary error. Plotting a line with the cumulative probability will produce two intersection points with the boundary of a -box. For , two intersection points are and . Two intersection points of are and . Thus the boundary error iswhere and . The illustration of the area error and the boundary area is shown in Figure 1.

The relative error is used for comparing the precision of different -boxes, which have the relative area error and the relative boundary error. Just to explain the definition, two -boxes are named, respectively, the relative -box and the contrastive -box. The definition of relative area error is the absolute error of two -boxes dividing the area between two bounds of the relative -box. Suppose that and are taken as the relative -box and the contrastive -box, respectively. The relative area error is defined by where is the area error between two -boxes. The relative boundary error in the probability value is computed bywhere is the boundary error between two -boxes. Two types of error are used to quantify the errors between -boxes, when the discrepancy of different -boxes can be distinguished from the figure.

#### 5. Numerical Examples and Discussions

##### 5.1. Sum of and

The test functionwhere and are random variables represented by dependent -boxes with the Gaussian copula. The Pearson correlation of the two random variables is 0.3. is bounded by normal distribution whose means are between 2 and 8 and whose standard deviation is 0.5. also is specified by the class of normal distributions with means [15, 20] and standard deviation 1.5. According to the statistics theory, the theoretical cumulative probability of can be calculated by the following equations:where , , , and . The numerical boundaries of are computed by the Simpson quadrature [22].

As a comparison, the DBC method with a range of the discretization numbers is adopted to compute this problem here. According to the Cartesian product, analyses would be required if two -boxes are discretized into focal elements, respectively. The absolute error values with the theoretical -boxes under the different discretization number are shown in Figures 2 and 3 by the red line with asterisk markers. Now, use the present method with the sample sizes to solve the same problem, and their absolute error values are plotted in Figures 2 and 3 with the blue line with filled circular markers.

**(a)**Boundary error lines using two methods with

**(b)**Boundary error lines using two methods with

**(c)**Boundary error lines using two methods withFrom the area error lines in Figure 2, two lines have the obvious convergence procedure with the increasing iterations, while the area error of the present method has the faster convergence procedure. Under the same calculated amount, the circular-marked values are always under the asterisk-marked values. In other words, the present method has more accuracy than the DBC method. Figure 3 shows the boundary error line under three boundary points , and 0.975 with two different methods. These boundary error lines have the trend of approximate 0 with increasing the calculation time. In three boundary points, the filled circular-marked error lines are always below the asterisk-marked error lines, except little points when that the asterisk-marked values are slightly bigger than the circular-marked lines. Using the present method will produce less error than using the DBC method under the same computation time; thus the present method has more calculation efficiency. For the multiple random variables, this advantage is more obvious. The calculated amount using the DBC method is exponential growth with increasing number of random variables and discretization number. This numerical test demonstrated the validity of the proposed present method. The calculation of structural responses is performed in the following section using the present method.

##### 5.2. A 10-Bar Aluminum Truss

A linear elastic 10-bar aluminum truss [23] is shown in Figure 4. The certain parameters for the 10-bar aluminum truss are the density kg/m^{3}, the length of the horizontal and vertical bars m, and the bars areas mm^{2} . The truss is subjected to two vertical forces and , and a horizontal force , which are uncertain variables of the structures. The Young modulus is treated as the random variables. The interesting system response is the vertical displacement of the joint 2. The uncertain parameters presented by -boxes are shown in Table 1.

, , , and are truncated with as the design domain of input data. LHS method is applied to generate sampling points and finite element method is used to solve . RBF approximation model is constructed, which can be written aswhere and are dependent and the Young modulus and the force are independent variables. Dependence in and is represented by the Gaussian copula with Pearson coefficient 0.9. The probability bounds of are computed by the present method involving by 10^{5} replications. For comparison purpose, we assume that , , , and are mutually independent due to lacking enough dependent information. The result of is also calculated by the present method with 10^{5} sample size. Two -boxes of considering dependent or independent input data are plotted in Figure 5 which are called and , respectively. is regarded as the theoretical value and is taken as proximate value due to lacking enough dependent information. The relative errors between and are shown in Table 2. It can be found that the relative area error is 3.43%. When or , the boundary error are 8.24% and 7.85%, respectively. It can be noticed that the relative boundary error is close to zero when . From this example, the dependence in the input -boxes has great impact on probability bounds of structural responses, which is particularly obvious in the tail of probability bounds.

##### 5.3. A 25-Bar Steel Truss

This case is a 25-bar steel truss [24] which is shown in Figure 6. The length of each horizontal or vertical bar is 15.24 m and the cross-sectional area of all bars is 5000 mm^{2}. The Poisson ratio is 0.3. The joint 12 is hinge-supported, and the joints 6, 8, and 10 are roller-supported. The joints 7, 9, and 11 are subjected to the vertical loads , , and , respectively, and the joint 1 is subjected to a horizontal load . The forces , , , and and the Young modulus are treated as -boxes, which are shown in Table 3. The interesting value is the vertical displacement of the joint 9. RBF is used to replace the finite element model, which can be written as where , , , and and are truncated at the probability of 0.05 and LHS is used to sample points.

Four cases are designed to analyse the influence from dependence in input variables and shown in Table 4. For cases 1–3, each of two forces is mutually dependent and their copulas are Frank copula. is regarded as the independent variable. The relatively large correlation coefficient 0.9 for random variables is used. In case 4, five input variables are assumed to be mutually independent. -boxes of in four cases are obtained by the present method with 10^{4} samples that are shown in Figure 7, where -boxes of with the condition in cases 1–3 are taken as the theoretical value and -box of with case 4 is used as approximate value. The relative area error between the theoretical value and approximate value is listed in Table 5. By comparing the relative error under different cases, the relative area error or the relative bounds area under case 2 is more than the error under case 1 and case 3. For each case 1–case 3, the relative bounds error for is less than the relative bounds error for or , which is very small and nearly 0. Overall, dependence under case 2 has more influence for the structural response of in this example. This example demonstrates that neglecting the dependence in input variables will produce unpredictable error for bounds of structural response.

**(a)**-boxes for with case 1 and case 4

**(b)**-boxes for with case 2 and case 4

**(c)**-boxes for with case 3 and case 4##### 5.4. Application to a Taylor Impact Test

Taylor impact tests are often performed to confirm the plastic behavior of a material under severe strain conditions. Extremely high plastic strains and strain rates occur at the crushed end of the rod, resulting in severe local deformation. The experimental data usually consist of measurements of its final deformed profile. As shown in Figure 8, the cylindrical sample total length cm and the diameter cm. Testing material is the high-strength low-alloy (HSLA-100) steel [25]. The modulus of elasticity is taken to be GPa and the material density kg/m^{3}. The initial velocity of the sample is 246 m/s. For the Taylor test, the final deformed profile is an important index to evaluate the material performance. Thus in the application the uncertainty response analysis is conducted for the final total length and the final diameter . To compute the two values, the finite element method is employed, and Zerilli-Armstrong (Z-A) modelis used for rate-dependent plasticity, where is the equivalent plastic strain, is resulting stress, and is the temperature at the start of the experiment. The parameters are material-specific. The Belytschko-Tsay shell element is used to create finite element model with 2000 elements. The initial temperature is 298 K. The time duration of the impact process is 20 ms.

Due to the manufacturing and measuring errors, Z-A model parameters are thus treated as -boxes variables, which are given in Table 6. There is dependence in these parameters [26] and their correlation coefficients are given in Table 7.

According to the present method, the random variables are truncated at firstly. For generating the sample points, LHS method is used. The approximation models of and are constructed by RBF, which are expressed as The analysis results of dependence are plotted in Figures 9 and 10, respectively. From this example, the present method is also effective for the nonlinear model including dependent uncertain parameters.

#### 6. Conclusion

A computational method has been developed for the calculation of the structural responses under uncertain parameters represented by dependent -boxes. The -boxes are discretized into a range of intervals based on the copula sampling technique. The calculation of structural response interval is carried out by the metamodel-based optimization method. The resulting intervals are then assembled into a -box. The present method takes both the advantages of efficiency of sampling method and the advantage of handling dependent -boxes of DBC method.

By investigating dependence in the structural uncertainty parameters based on the present method, dependence in the variables cannot be neglected. When the dependent variables are regarded as the independent variables, the probability bound of structural responses will produce error, which are unpredictable and may lead to the failure decision-making. We have demonstrated, through the use of the present method for solving a benchmark problem and three numerical examples, that this method is effective and practical for the structural response analysis.

#### Conflict of Interests

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

#### Acknowledgment

This work is supported by the grant from the Fundamental Research Funds for the Central Universities (Grant no. 227201401203).