Abstract

This paper presents a novel procedure based on first-order reliability method (FORM) for structural reliability analysis with hybrid variables, that is, random and interval variables. This method can significantly improve the computational efficiency for the abovementioned hybrid reliability analysis (HRA), while generally providing sufficient precision. In the proposed procedure, the hybrid problem is reduced to standard reliability problem with the polar coordinates, where an n-dimensional limit-state function is defined only in terms of two random variables. Firstly, the linear Taylor series is used to approximate the limit-state function around the design point. Subsequently, with the approximation of the n-dimensional limit-state function, the new bidimensional limit state is established by the polar coordinate transformation. And the probability density functions (PDFs) of the two variables can be obtained by the PDFs of random variables and bounds of interval variables. Then, the interval of failure probability is efficiently calculated by the integral method. At last, one simple problem with explicit expressions and one engineering application of spacecraft docking lock are employed to demonstrate the effectiveness of the proposed methods.

1. Introduction

In recent years, reliability analysis has been a research hotspot of structural design and analysis. The influence of uncertainty arising in loads, material properties, dimensions, and geometry becomes more and more profound [1]. Probability theory is one of the most universal tools to estimate structural reliability and safety by calculating the probability of failure Pf and the reliability index β [24]. Hence, FORM is the traditional and most commonly used method in the field of structural reliability, and it has been widely employed in a wide range of industrial fields [58].

In FORM for structural reliability analysis, by the first-order Taylor series, the limit-state function is approximated at the design point in the standard Gauss space of transformed independent random variables [9, 10]. In this space, reliability analysis is usually converted into an optimization problem to solve the reliability index [8]. A review of FORM can be found in Breitung [11]. Its application and accuracy in most structure engineering are well recognized, for instance by Rackwitz [12]. Derived from FORM, many probability approximation methods have further developments and more refined alternatives [1115]. However, some problems may arise in the application of these methods, such as accuracy, convergence, and computational cost [1113]. Specially, there is a defect that both these methods require precise probability distributions of the random variables based on a great amount of experimental samples as in many engineering applications and the samples are often limited. It causes the distributions of some variables not being precisely known. Generally, we can only obtain the variation ranges of some parameters based on the limited data and descript the stochastic characteristic of them with subjective assumptions [16]. It is likely to bring about a serious error in structural reliability analysis [1720].

In real-world engineering, it is obvious that random and interval variables coexist frequently [21, 22]. Quantifying both types of the abovementioned uncertainty is called hybrid reliability analysis (HRA) [1]. HRA has been an important issue in the research field in structural reliability [22]. In recent years, many research studies have been carried out on the HRA methods for structure [2326]. The following numerical methods have been proposed: the function approximation technique [27], the probabilistic transformation technique [2830], the iterative rescaling method [31], the probability bounds approach [32], the mixed perturbation Monte Carlo (MC) method [33], the optimization algorithm with single-layer nesting [34, 35], two-layer nesting [36, 37], and the complex nesting [22], among others [3840]. Nevertheless, the HRA with FORM is still in its infancy, as there are some calculating problems in the practical applications. Algorithms with high efficiency and robust convergence must be developed [21, 22].

Hence, a new procedure, based on first-order reliability method (FORM), is proposed to deal with random and interval hybrid variables. This procedure can significantly improve the computational efficiency for the abovementioned HRA problem, preserving sufficient precision. In the proposed method, we involve the distance coordinate and the cosine of the angle coordinate in a polar space. In this space, the hybrid problem is reduced to standard reliability problem, where an n-dimensional limit-state function is defined only in terms of the two random variables. And the probability density functions (PDFs) of the two variables can be obtained by the PDFs of random variables and boundaries of interval variables. Finally, the interval of failure probability is efficiently calculated with the integral method, for the new bidimensional function.

The main contributions of this paper lie twofold in the following: (1)In the framework of the proposed method, the novel procedure is more useful and efficient in analyzing the hybrid reliability problem with randomness and interval.(2)The new PDFs of two nonlinear features have been deduced, using the probability distribution of randomness and the range of intervals. It is distinct from the assumption of uniform density function model for interval variables. The new PDFs of the two features contribute to the evaluation of the lower and upper probability of failure, respectively.

The remainder of this paper is organized as follows. A brief summary of the classical structural reliability analysis with FORM is presented in Section 2. The whole efficient reliability analysis method has been proposed in Section 3. In Section 3.1, the polar transformation with random and interval uncertainty has been discussed, and the distribution of two new variables in the polar space has been derived. The new calculated models and algorithm procedures for the range of failure probability have been proposed in Section 3.2. Finally, two numerical applications to demonstrate its efficiency and practicability are presented. Firstly, the hybrid reliability problem with one simple linear explicit expression is used to prove that this method is effective and accurate, due to low analytical complexity of the linear limit-state function. Secondarily, the method is applied to deal with the hybrid reliability problem in a real structural system, in this specific case a spacecraft docking lock, to demonstrate the effectiveness of the proposed method.

2. Traditional Structural Reliability Analysis

In the reliability design and analysis of structures, the limit-state function is always constructed according to specified requirements. Here, the n-dimensional variable vector denotes the independent random variables that affect the function, such as load, material properties, geometry, dimensions, and environmental factors. can be used to determine if the structural system works normally by checking whether or not. The probability of occurrence of failure event Pf is obtained as follows [2, 3]: where is the probabilistic measure. According to the dual event , the probability of reliability can be formulated as follows: which can be computed by the following integral: where denotes the joint PDF of the random variables X. Generally, the above integral cannot be solved analytically because of its complex integration boundary and the high dimensions of X. Therefore, some approximation techniques have been well established to compute the integral efficiently, which include the well-known first-order reliability method (FORM) [16]. In FORM, by a following transformation from the original space of X to an independent standard normal space of u, the limit-state function becomes as follows: where T is the transformation that converts u into X. Additionally, this equation will in the sequel be represented simply as with the understanding that there is an underlying transformation to the Gaussian space u from the given space X.

Thus, the integration in (3) can be rewritten as

is the n-dimensional standard Gaussian density function. By using the reliability index approach, the above equation can be further changed to the following optimization problem [23]: where denotes the norm of a vector. The optimum u is called the design point where the joint PDF reaches a maximum value among all of the points on the limit-state function. is a reliability index which geometrically represents a minimum distance from the origin to the limit-state surface in u space.

Despite that this problem and its solution are well known, for the ensuing developments, it is necessary to recall its main steps. The first-order Taylor expansion of the limit-state function about a general point u+ can be written as follows: where stands for evaluated at u+. In FORM, the design point is the solution of the optimization problem (7). Its value is as follows: where is the design point unit vector, which is normal to the tangent hyper plane at the design point.

Moreover, we can suggest that a source of difficulties for solving the structural reliability problem by any method is its dimensionality, determined by the number of random variables n [41]. In essence, this is due to the fact that the estimation of the normally very small probability failure requires a large number of samples in high dimensional spaces, either for establishing a surrogate of the limit-state function or for MC solution. This also makes the visualization of the reliability problem impossible, which is a desirable goal for the design practice [29]. Therefore, some methods are proposed for reducing the dimensionality of the problem to only two dimensions, by means of a polar representation for random samples used in MC simulation, based on the vector of the center of mass of the failure probability and the design point [23]. There are two nonlinear features: (a) their distance to the origin and (b) the cosine of the angle they make with the design point unit vector α. Moreover, these methods are developed to estimate the reliability of a structure system in which the input variables are modeled using any representation provided by possibility distributions, intervals, probability boxes, CDFs, or Dempster-Shafer structures based on random set theory [29]. The two new variables are given as follows:

For the large number n, v1 obeys a chi distribution [41]:

and v2 obeys the following distribution:

It was derived by Hurtado that the two new variables are independent [23, 41]. In order to avoid the serious error of the reliability analysis brought by subjective assumptions on description of the uncertainty characteristics, such as the uniform distribution or truncated Gaussian distribution representation for interval uncertainty, the new polar transformation with random and interval uncertainty is presented in next the section.

3. Efficient Reliability Analysis Method with FORM

3.1. Polar Transformation with Random and Interval Uncertainty

In structural reliability problem, let us assume that the limit-state function is , where is a n-dimensional vector of random variables and is a m-dimensional vector of interval variables. In terms of the standard independent Gaussian vector and the normalized independent interval vector , the limit-state function becomes as follows: where T is the standard Gaussian transformation that converts u into x and is the normalized interval transformation that converts δ into y. In order to keep the notation uncluttered, this equation will in the sequel be represented simply as

With , the structural reliability problem is defined in the n + m dimensional space w of n independent standard Gaussian variables and m normalized interval variables. The main purpose of this method is to transform it to a problem of only two dimensions. By (10),, and, we can use the polar transformation from space ω to the polar space by the two nonlinear features: where is the solution of the new optimization problem

Hence, in the polar space, the distribution of two new variables can be derived as follows.

For (N is the nonnegative integer domain), n is the number of the random variables and m is the number of the interval variables. (1)Like v2 in Section 2, obeys a distribution: (2)For let , with ; it is derived that . As is the independent standard Gaussian variable, obeys the distribution shown as follows: where n is the number of the random variables. Hence, for , obeys the distribution written as follows: where m is the number of the interval variables. And the cumulative distribution function (CDF) of with ∆ denotes . It can be obtained that is nonincreasing for the variable ∆.

Proof. Let , the partial derivative of for ∆ can be written as follows:

Considering , , , and , thus it is obvious that . Hence, is nonincreasing for the variable ∆.

Moreover, because of the two new variables independent in the polar space, with the limit-state function , the probability of failure in the structure can be written as follows: where obeys (19) and obeys (17).

Considering the nonincreasing of CDF of for ∆, the interval of is ; the lower boundary and the upper boundary can be calculated by where obeys (20) and obeys (18).

3.2. The Efficient Method with FORM

Let us consider the first-order approximation at the design point of the limit-state function given by (8) with the last term disregarded, where obtained by the optimization problem (16). can be written as follows:

Let where

Equation (26) can be written as follows by taking (24): where

Distinctly, ; in the polar space, the new safe domain becomes , the failure domain becomes , and the limit-state surface is .

We can calculate the interval of the failure probability in the structure with (29) and (31) as follows. (1)The lower boundary of failure probability :(2)The upper boundary of failure probability : where obeys (19) and obeys (17).

The method proceeds as follows: (1)For the random variables and interval variables in the limit-state function of a structure, convert into the standard independent Gaussian vector and into the normalized independent interval vector . The limit-state function becomes with .(2)Let the initial design point be , obtained by the optimization problem (16).(3)Calculate the gradient of the function at the design point .(4)Transform the function into in the polar space and calculate the parameters(5)In the polar space, use (32) and (33) to calculate the initial interval of failure probability of the structure written as follows:(6)Calculate the new gradient by the following:(7)Calculate the new design point with by the following:

In the polar space, the parameters and by (35) and (36), respectively, are written as follows:

And the new limit-state function is . (8)Calculate the interval of failure probability of the structure by the following:(9)For the robust convergence of this procedure, let the iterative error of the calculation be . If and , stop calculating.Otherwise, go to step 2. Set the new design point to be the initial design point.

4. Numerical Examples

4.1. Numerical Example with Linear Limit-State Function

(1) In the structural reliability problems, let the limit-state function of a structure be , with the random variable obeying the standard normal distribution N(0,1) and the interval variable in the range of [0,1]. They are independent of each other.

With the polar transformation, let

And can be transformed to in polar space, where the PDF of the random variable is (18), and the PDF of is (20).

We can see that uncertain variables are transformed to two uncertain variables in the polar space. Subsequently, we will discuss the calculated interval of failure probability with different n and m as shown in the following tables. For comparison purposes, a full-interval MC simulation involving the generation of 10,000 random realizations of input parameters u and the computation of the failure probability for each of them with 10,000 samples of vector was performed.

The relevant information about the two methods applied is summarized in Tables 1 and 2. Table 1 indicates that the estimates calculated by the proposed method in the paper are exactly in the interval Monte Carlo sense. Also, by the relative errors of the bounds in Table 2, it can be seen that the method succeeds in giving accurate reliability analyses with the maximum value 11.6%. This example shows the efficiency of the method, inasmuch as only 191 iterations (maximum value) were required in comparison to the 100,000,000 sample evaluations required by interval MC in the evaluation of and , achieving the same precision of the latter (interval MC).

On the other hand, from the same total number of uncertain variables with different n random variables and m interval variables, it is suggested that the number of interval variables influences the results stronger than randomness. From the different total numbers of uncertain variables, it is suggested that the number influences the efficiency (iteration).

4.2. Structural Latch of Space Docking Mechanism

Consider the structural latch of space docking mechanism shown in Figure 1. It is an important part of the docking mechanism, which is the key structure to realize the rigid connection and separation between two space vehicles. The connection force is generated by the interaction between the series of latches in two different docking mechanisms.

The working procedure of the workflow of structural latch is shown in Figure 2.

During the procedure of locking every structural latch, the load F is an interval variable with the range of [30, 40] kN. The modulus of elasticity of the material E is a normal variable with the mean 117,000 MPa and standard deviation 1170 MPa, and the cross section depth d is a normal variable with the mean 5 mm and standard deviation 0.05 mm. Also, there are two random variables Poisson ratio and the material density . The reliability problem is defined by the limit-state function where is the yield strength with the constant and is the contact stress on the connected surface of the two latches.

It is supposed that the four random variables are independent. A simplified finite element analysis (FEA) model is created to compute the contact stress σ (Figure 3), with 25,642 elements.

The proposed algorithm in this paper was employed to perform the new FORM analysis. After 127 iterations, the lower and upper bound of the failure probability were calculated as [4.3580e−06, 4.3930e−05]. It coincides with the result calculated using Monte Carlo simulation, inasmuch as all of the failure samples of the simulation were correctly identified. However, the relative errors are 8.12% and 11.4%, which indicate that this method employs the linear Taylor series to preserve the efficiency, sacrificing the accuracy requirements. The calculated result for the structural latch with comparison is summarized in Table 3.

The difference between the results of MC and the proposed method is caused by the nonlinearity of the limit-state function. However, the error of the upper bound is only 8.12% and that of the lower bound is 11.4%, which indicates that the nonlinearity of this problem is not too strong, but still directly influences the result precision. Compared with the relative error, the results of the presented method can be accepted. On the other hand, the evaluation number of FEA for the proposed method is only 223. Since each evaluation of the FEA model takes almost 10 min, the increased computational cost for MC would be more than 15 days. Therefore, if we pay more attention to the computational cost in this problem, then this proposed method in this paper is preferred, although losing better accuracy.

5. Conclusion

In this paper, an efficient method for the hybrid reliability analysis of structures for random and interval uncertainty is presented. It has been shown that the design point concept of probability-based reliability analysis can still be applied in this article. Specifically, the unit vector of the point is the direction of the rapid changes of the output variable implied in the limit-state function. In order to avail of this property and reduce the computation cost, a bidimensional mapping of two nonlinear features is established in the polar space: (a) the distance to the origin and (b) the cosine of the angle with the design point unit vector. In this case, the hybrid problem is reduced to standard reliability problem, where n-dimensional limit-state function with the linear Taylor series is approximated around the design point defined only in terms of two random variables, obeying new PDFs with the number of random and interval variables. They are distinct from the assumption of uniform density function model for interval variables and contribute to the evaluation of the lower and upper probability of failure, respectively. The numerical experiments confirm the solid theoretical foundation, the agreeable accuracy, and the low computational cost of the proposed methodology.

Conflicts of Interest

The authors declare that they have no competing interests.

Acknowledgments

This study was supported by the National Basic Research Program of China (no. 2013CB733002) and the National Natural Science Foundation of China (no. 51675026).