On High-Dimensional Time-Variant Reliability Analysis with the Maximum Entropy Principle

Zhou, Fuliang; Hou, Yu; Nie, Hong

doi:https://doi.org/10.1155/2022/6612864

International Journal of Aerospace Engineering

On this page

Abstract Introduction Conclusions Appendix Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 6612864 | https://doi.org/10.1155/2022/6612864

On High-Dimensional Time-Variant Reliability Analysis with the Maximum Entropy Principle

Fuliang Zhou,¹Yu Hou,²and Hong Nie¹

Academic Editor: Xingling Shao

Received21 Oct 2021

Revised07 Jan 2022

Accepted04 Feb 2022

Published19 Mar 2022

Abstract

The structural reliability analysis suffers from the curse of dimensionality if the associated limit state function involves a large number of inputs. This study develops a reliability analysis method that deals with high-dimensional inputs over time. The probability distribution of the structural response is reconstructed by the maximum entropy principle which is achieved by solving an optimization problem derived from the concept of relative entropy. The optimization problem is transformed into a convex one with respect to the orders of fractional moments and the Lagrange multipliers. Additionally, considering the associated computational issues, it is reformulated with side constraints on the parameters of the maximum entropy distribution. Then, a global optimization procedure is performed. The proposed method is successfully applied to the reliability analysis of a linear and a nonlinear structural system, which involves a large number of inputs deriving from the discretization of the input random processes.

1. Introduction

Structural reliability analysis has been intensively studied during the past few decades. Based on the consideration of time factor, structural reliability analysis can be roughly divided into two groups: time-invariant reliability analysis and time-variant reliability analysis. For the first group, as its name implies, the time-invariant reliability methods do not consider the effect of time factor, for example, first-/second-order reliability method (FORM/SORM) [1–3], moment methods [4, 5], response surface methods [6–9], and importance sampling [10, 11]. Therefore, they only take random variables as inputs. However, almost all structural systems involve time-dependent parameters during their period of operation and/or service, e.g., stochastic loading and material property degradation. Thus, time-invariant reliability analysis may provide inaccurate reliability assessment and prediction.

In the past decade, time-variant reliability analysis has drawn much more attention than its counterpart because of its capacity of dealing with time factor. Thus, the corresponding time-variant reliability methods can take random variables, random processes, and even random fields as inputs. Many time-variant reliability methods have been developed in recent years. They can be roughly classified into three categories: (1) the out-crossing-rate-based methods [12–18], (2) the extreme value methods [19–27], and (3) the composite limit state methods [28–32]. A common feature of these methods is that a time-variant reliability analysis problem is usually transformed into a time-invariant one or broken down into a series of time-invariant ones. Then, the aforementioned time-invariant reliability methods are applied to complete the time-variant reliability analysis. For instance, FORM has been integrated into many time-variant reliability methods [14, 17]. Note that the time-variant reliability methods with FORM inherent the drawbacks of FORM.

Many challenges are encountered by time-variant reliability analysis. One main challenge arises from the curse of dimension when a large number of random variables are used to represent the input random processes and/or fields. It should be pointed out that there is also a universal and dimension-free tool for both time-invariant reliability analysis and time-variant ones. It is direct Monte Carlo simulation (MCS). However, MCS suffers from a huge computational effort for estimating a small failure probability when a time-consuming model or experiment is involved, e.g., a large-scale finite element model. Subset simulation [33, 34] is an alternative for MCS for estimating small failure probabilities in both time-invariant reliability analysis and time-variant ones. However, several thousands of samples are still required to obtain a robust and accurate estimation by subset simulation.

In this paper, we are interested in developing an efficient and accurate method for high-dimensional time-variant reliability analysis. The high-dimensional issue is solved through the development of a reliability method using statistical moments.

Zhao and Ono’s moment method is the most attractive one for structural reliability analysis [4, 5] after the year 2000, where the first three or four moments of structural response are used to fit a parametric distribution. The Pearson system [35] was recommended as the reliability index for the moment method. Although the Pearson system of frequency curves is flexible, they still impose spurious information on the unknown distribution. It is also contentious from both theoretical and practical views [36, 37]. The maximum entropy (MaxEnt) principle was pioneered by Jaynes [38, 39], which minimizes the amount of prior information built into the approximate distribution while maximizing entropy, i.e., a measure of uncertainty. Using the MaxEnt principle, Li and Zhang proposed a combined method for structural reliability analysis [37]. The statistical moments of structural response are calculated by the dimension reduction method [40, 41]. Zhang and Pandey analyzed two challenging issues of using integer moments when approximating the unknown distribution by the MaxEnt principle [36]. They are a large number of required moments and the oscillatory phenomena at the tail of MaxEnt distribution. The usage of fractional moments [42] was proposed to overcome these two issues for time-invariant reliability problems since a fractional moment contains the statistical information of a large number of integer moments. A multiplicative dimensional reduction method was also developed to compute the fractional moments of structural response [36]. Considering correlated input random variables, Li et al. suggested using the unscented transformation technique to calculate fractional moments [43]. Dai et al. employed the MaxEnt principle to construct the optimal importance sampling density for the importance sampling method [44]. However, the above investigations are limited to time-invariant reliability problems.

For the time-variant reliability analysis problem, Shi et al. presented the performance of the combination of dimension reduction method and the MaxEnt principle with fractional moments [45]. Xu and his colleagues carried out a series of studies about the MaxEnt principle with fractional moments and their applications for time-variant reliability analysis [46–50]. However, the number of input random variables in the application examples is not beyond 20 [45–50]. Thus, the abovementioned methods cannot be applied to the time-variant reliability analysis problems involving a large number of inputs [47].

The objective of this paper is to propose an efficient method for the high-dimensional time-variant reliability problems by using the MaxEnt principle with fractional moments, where the number of input random variables is larger than 100. The contributions of this study are as follows: (1) a new perspective for estimating a failure probability curve from the combination of fractional moments and the maximum entropy principle and (2) a new tool for high-dimensional and nonlinear time-variant reliability analysis with input random variables and processes. High-dimensional issue is very challenging in many science and engineering fields [51–53]. To address this issue, fractional moment and the maximum entropy distribution are used to approximate the probability density function of structural responses of interested, which has a capacity of dealing with zero-shot events [54] for structural reliability analysis. In other words, the proposed method can estimate structural reliability with zero-failure samples.

The remainder of the paper is organized as follows: Section 2 briefly reviews the problem setup for time-variant reliability analysis. Section 3 describes the procedure and theoretical rationale of the proposed method. Section 4 analyzes the implementation issues and their corresponding possible solutions. Two high-dimensional examples are used to demonstrate the performance of the proposed method in Section 5. Concluding remarks are provided in Section 6.

2. Problem Statement

Limit state functions (LSFs) are employed to capture the intended operating condition of structures in the structural reliability community, which divides the operating space into two domains, i.e., the safe domain and the failure domain. Let be an LSF of the intended operating condition of a structure, where indicates the structure is safe while means there is a failure in the structure. The general form of LSF is defined as where is the response quantity, is the input random variable vector, is the input random process vector, and is the input time, respectively. Based on the different type of inputs in the LSF , the LSF may degrade into simpler types [22]. It is well-known that reliability is a function of time snice the LSF of interest is related to time via the time-dependent random process and the explicit operating time . Suppose that we are interested in the time-variant reliability of the structure over a time interval , where is the final time instant. The time-variant reliability can be expressed as where the symbol means for all time instants. In practical estimation, the corresponding time-variant failure probability is more preferable which is defined as a complementary quantity to the time-variant reliability, i.e., where the symbol has a meaning of “there exists at least one.” To distinguish the failure probability over a time interval from that for a time instant, in Equation (3) is also called as the cumulative failure probability. Obviously, it is a function of the final time .

To handle the input random processes, the time interval of interest is usually divided into a set of discrete time instants with a constant time step . Let represent the time instant at the th step, where the index . Spectral decomposition methods are suggested to deal with the input random processes and transform them into a set of random variables [17, 22, 24, 26]. Let be the transformed random variables from the input random processes. Then, the LSF becomes . Furtherly, at a time instant , an instantaneous LSF is defined as . Using the concept of series-system reliability, the cumulative failure probability in Equation (3) is approximated as where . The advantage of Equation (4) is that the output random process is converted into a random variable .

The magnitude of time step size affects the approximate accuracy of input random processes with a finite number of random variables firstly. Then, this effect may propagate into Equation (4) through the LSF and furtherly influence the estimation accuracy of the cumulative failure probability. In this study, we assume that the magnitude of time step size is sufficiently small and adequate to capture the uncertainty in the input random processes and failure information in the transformed series system. The current study focuses on the development of a new high-dimensional time-variant reliability analysis method based on Equation (4).

3. Proposed Method

3.1. Overview

The new high-dimensional time-variant reliability analysis method has three stages. The first stage is to generate random samples for the input random variables and random processes using MCS sampling strategy and furtherly compute the structural responses. The second stage involves the approximation of the MaxEnt distribution. The third stage is related to the calculation of failure probability using the MaxEnt distribution obtained in the second stage. It is clear that the second stage is the key stage of the proposed method. The implementation procedure of the proposed method is shown in Figure 1.

Stage 1 involves four steps. Step 1: discretize the time interval with a constant time step size and obtain discrete time instant . Step 2: decompose the input random processes with a group of random variables . Step 3: generate random samples for the input random variables and random processes. Step 4: substitute the random samples into the LSF to obtain trajectories of structural response.

Stage 2 is furtherly divided into three steps. Step 1: solve the optimization problem derived from the concept of relative entropy. Step 2: calculate the fractional moments from the trajectories of structural response. Step 3: determine the MaxEnt distribution based on the solution to the optimization problem in Step 1. In the subsequent sections, we will discuss the details of Step 1 and 3 since these steps are theoretical fundamentals of the proposed method.

Stage 3 is straightforward once the probability density function (PDF) of MaxEnt distribution is determined in Stage 2. The estimation of a cumulative failure probability is just a numerical integration over the failure domain.

3.2. Maximum Entropy Principle with Fractional Moments

The advantages of usage of fractional moments in the maximum entropy principle have been proven by several studies in the structural reliability community [36, 43, 45–50]. This subsection will give a brief review of the maximum entropy principle with fractional moments.

Identifying the probability distribution with statistical moments has been investigated extensively in the literature. However, the use of integer statistical moments suffers from accuracy concern. Zhang and Pandey study the need of fractional moments instead of integer moments, which is evoked by the relationship between them [36].

We start with the definition of the fractional binomial coefficient that is given by where is a real number and is an integer number. If applying the Taylor series expansion on a function of , one has where is a reference point.

The -order fractional moment of a positive random variable is defined as

Substituting Equation (6) into Equation (7), one has where is the th central moment of the random variable . If the reference point is put on the mean value of the random variable , the th central moment can be rewritten as

Equations (8) and (9) indicate that a fractional moment includes the information of a series of integer statistical moments, either central moments or raw moments.

If a couple of fractional moments are available at hand, they would provide much more information than the same number of integer moments for identifying an unknown distribution type. In contrast, Gzyl and Tagliani [55] proved that 30 or more integer moments are required by Equations (8) and (9) to guarantee the accuracy of a fractional moment, where the infinite sequence of integer moments is demanded theoretically. This result is extremely important and gives a strong motivation to employ the concept of the fractional moment to reconstruct an unknown distribution type from a group of samples.

The information entropy of a random variable is a measure of the uncertainty containing by the random variable, which is defined as where is the probability density function of the random variable .

The core step of the proposed method in this study is to reconstruct an expression for the unknown PDF of the structural response , given a sample of experimental data. The MaxEnt principle states that the most unbiased estimation of the unknown PDF is the one that maximizes Equation (10) with given constraints. In this study, the given constraints are the fractional moments determined from the experimental data. Therefore, by the MaxEnt principle, the PDF identification problem is converted into an optimization problem [42]. where is the fractional moment estimated from a sample

Note that are the random samples for the random variable and is the number of experimental points in the sample.

Solving the optimization problem with the Lagrange multiplier method, one has the following generic form of the MaxEnt distribution [36, 42]. where denotes the MaxEnt PDF for , are the Lagrange multipliers, are the orders for the fractional moments, and is the number of available fractional moments, respectively. The parameter is a normalizing factor which makes sure that the integration of over its support domain is 1 and satisfies the definition of PDF. Thus, from Equation (13), one has

Based on the concept of relative entropy (or the Kullback-Leibler distance), Zhang and Pandey proved that the following optimization problem is equivalent to Equatio (11) [36].

The minimization of the relative entropy means the minimization of the generalized distance between the true PDF and the estimated PDF , which is quantified by the information entropy. We will further reformulate the optimization problem in Equation (15) and derive an attractive property of it in the next subsection.

3.3. Computation of Maximum Entropy Distribution

Recalling the equivalent optimization problem in Equation (15), its objective function is furtherly rewritten as

It is obvious that minimizing is equivalent to minimize the following integral:

The integral term is defined as a function for the purpose of simplification. Based on the definition of in Equation (17), the MaxEnt distribution can be rewritten as

This study examines the properties of the integral function firstly. The derivative of with respect to is written

Recalling the definition of fractional moment, we have

The first derivative of with respect to is given by

From Equations (20) and (21), we know that the solution point to Equation (17) is a saddle point of at least.

Next, we examine the elements in the Hessian matrix of . There are four types elements in the Hessian matrix, i.e., , , , and which consist of the submatrices , , , and , respectively, as shown in

For the submatrix , the elements are

It is apparent submatrix is a covariance-like and symmetric one. In a similar way, we have

Equation (24) indicates that submatrix is also a covariance-like and symmetric one. For an element in submatrix , it has an expression

Similarly, an element in submatrix has an expression

It can be seen from Equations (25) and (26) that submatrices and are also covariance-like ones. An attracting property of submatrices and is that the columns of are the rows of , i.e., is the transpose of .

According to above discussion, the Hessian matrix of is proven to be a covariance-like and symmetric one. As long as the investigated is a PDF, the Hessian matrix is positive definite and of full rank. Thus, the solution to Equation (16) is the minimum point of and furtherly of . Since the optimization problem in Equation(17), (15) is convex and has no constraint; it seems to provide an opportunity to solve it easily.

4. Computational Issues and Solutions

We have discussed the theoretical rationale of the proposed method in Section 3. However, the technical rationale of the proposed method is more important for implementing this algorithm in practical applications.

4.1. Problem Reformulation

We have proved that the core optimization problem in Equation (17) is a convex one and there is no constraint on it. That is why most of studies have suggested that the optimization procedure can be carried out by the Nelder-Mead algorithm to search the solution [36, 45–50]. However, its technical implementation is not as perfect as the theoretical derivations, based on a lot of numerical implementation experience.

The first reason is that the Nelder-Mead algorithm is only suitable for a small number of optimization variables. As the number of optimization variables increases, the performance of the Nelder-Mead algorithm becomes worse. The current study focuses on the high-dimensional time-variant reliability analysis; then, the number of fractional moments is set to be 3 or larger empirically. Thus, the total number of optimization variables is 6 or larger. We also encountered the phenomena of diversity due to different initial values of the optimization variables. This may be caused by the incapability of the Nelder-Mead algorithm.

The second reason is that the estimate error of fractional moment increases with the order of moment, i.e., [36]. For this reason, the lower-order moments are used in the classical maximum entropy method. To overcome this serious drawback, the values of should have limits during the optimization procedure. For example, was suggested to obtain an unbiased estimation from experimental data [36] and was used by Inverardi and Tagliani [42]. In this study, is used. Additionally, the magnitude of the Lagrange multipliers should also have limits on their values to reduce the searching space and the numerical error associated to the numerical computation. We put the restriction of on the Lagrange multipliers.

The third reason is the numerical error associated with the integral calculation involved in Equation (17). Furthermore, it affects the accuracy of , PDF and CDF of the MaxEnt distribution, and further reliability analysis. The CDF of the MaxEnt distribution may be larger than 1 if the improper fractional exponents and Lagrange multipliers are obtained. A high-order global adaptive quadrature is used to perform the integral calculation in this study, which possesses a relatively high accuracy among the numerical integration algorithms in the literature.

Except the above reasons, there is another aspect that is worth to be discussed since it is directly connected with the interval of integration and further the associated accuracy. That is the support of the MaxEnt distribution. Since the samples of structural responses are available when determining the support of the MaxEnt distribution, the maximum, minimum, mean, and standard deviation can be estimated from the samples easily. To avoid arbitrariness, the well-known principle is suggested to determine the support of the MaxEnt distribution. The upper limit is set to be , where and are the maximum and standard deviation of the sample, respectively. The lower limit is set to be if is larger than 0, where is the minimum observation in the sample. Otherwise, it is set to be 0 as the investigated response quantity has positive support which is required by the definition of fractional moment.

Integrating the above discussion, the reformulated optimization problem has side constraints on the fractional exponents and the Lagrange multipliers. Thus, subset simulation optimization [56, 57] is employed to solve the reformulated optimization problem, instead of the Nelder-Mead algorithm in this study.

4.2. Parameter Estimation Using Subset Simulation Optimization

Subset simulation optimization is a heuristic algorithm for optimization problems, which is extended from a popular reliability method, i.e., subset simulation [56–58]. An attractive feature of subset simulation optimization is its ability for high-dimensional problems.

Based on the discussion in subsection 4.1, we have reformulated the optimization problem in the MaxEnt principle as an unconstrained one with side constraints on optimization variables. In subset simulation optimization, a truncated normal distribution is assigned to each optimization variable according to the associated side constraint. Take a moment order as an example, the PDF of the corresponding truncated normal distribution is given by [56–58] where and are the upper and lower limits of the moment order, respectively; and are the mean value and standard deviation determined by the 3-sigma rule in reliability engineering; and and are the PDF and cumulative distribution function of standard normal distribution, respectively.

Since the side constraints are considered in the truncated normal distributions, there is no need to handle them during the optimization procedure. Along with the spirit of subset simulation optimization, subsets are generated to approach the optimal solution progressively. All subsets have a generic form as where are the thresholds for the objective function in Equation (17). If one arranges the thresholds in a descent order, i.e., , the subsets have a nested structure

Thus, the final subset in the sequence would be or include the minimum solution theoretically.

In subset simulation optimization, the value of the objective function is the driving quantity for optimization, while the optimization procedure is achieved by the evaluation of the level probability that is an optimization parameter for subset simulation optimization. The conditional probability between two contiguous subsets is defined as the level probability and simulated by Markov Chain Monte Carlo simulation, i.e., where is the number of candidate designs (or samples) belonging to the subset , and is an indicator function which is used to check a candidate design falling in or not. In this study, we adopted a strategy in which the value of is reduced as the number of optimization iterations (or simulation level) increases.

The standard deviation of the samples in the th simulation level is employed to check the convergence of the optimal searching. The estimator of standard deviation tends to zero while the searching process is approaching the global minimum. At the same time, we also impose a maximum number of simulation levels as another stopping criterion to avoid the overflow of the affordable computational effort.

More details of subset simulation optimization can be referred to Refs. [56–58].

5. Test Problems

Two dynamic reliability analysis problems are used to demonstrate the performance of the proposed method. The first problem involves a single-degree-of-freedom (SDOF) oscillator which has a linear behavior, while the second problem is a five-DOF structural system with nonlinear behavior.

5.1. A SDOF Oscillator

The investigated SDOF oscillator is shown in Figure 2 with a white noise excitation . The governing equation of this structural system is given by where the natural frequency is and the damping ratio is , respectively. The spectral intensity of the white noise takes 1. The following initial conditions are applied to the system:

The failure event of this system is defined as its maximum displacement exceeding a threshold over a time interval s where is the total number of discrete time instants to represent the Gaussian white noise. Ref. [33] used a constant time step size for discretizing the time interval of interest, which leads to the total number of time instants be . Then, 1501 input standard normal random variables are employed to represent the input white noise, i.e.,

To meet the definition of the fractional moment, the structural response associated to the failure event in Equation (33) is rewritten as where is the normalized structural response. Then, the proposed method was applied to solve this problem with a sample size , 1000, 1500, and 2000. For the purpose of comparison, an MCS with 10⁶ samples was carried out to obtain a reference result.

Figure 3 shows the comparisons of the fractional moments and the relative error for the normalized response in Equation (35). Here, the order of fractional moments varies in the interval of as we discussed in Section 4.1. The fractional moments calculated from different sample sizes almost overlap one another (Figure 3(a)). However, it is clear that the relative error is becoming small with the increase of sample size (Figure 3(b)). Furthermore, the biggest relative error is 1% when the order of fractional moment takes -1.5, which is close to the used values reported by Xu and Wang [48] using a low-discrepancy sequence. Therefore, the relatively accurate fractional moments are employed in the proposed method by limiting the scale of the order exponents.

(a) Fractional moment

(b) Relative error

Inverardi and Tagliani proposed a penalty term , which is proportional to the order of the MaxEnt distribution and inversely proportional to the sample size , to avoid reconstructing a “too elaborate” distribution that cannot be justified by the given experimental data [42]. Table 1 lists the values of the objective function in Equation (17) while Table 2 gives the values of the objective function considering the penalty term. Six digitals must be kept in order to observe the difference among the results for a sample size in Table 1, while four digitals are required in Table 2. The value of the objective function increases as the number of fractional moments used in the MaxEnt distribution and decreases with the number of samples adopted for the estimation of fractional moments. Thus, is preferable for approximating the MaxEnt distribution of the normalized structural response, and 2000 samples provide the most accurate estimation of fractional moment and further reduce the objective function value.

Figures 4 and 5 show the histograms, PDFs, and probability of exceedances (POEs) [36] obtained by the proposed method and MCS. Since all settings of parameters produce very similarly results, we only report the results for with (Figure 4) and 2000 (Figure 5). The term “Exp” refers to the experimental data (sample) used to estimate the fractional moments. The results of the proposed method agree with those of MCS, low to a very small probability (about 10^-6). Thus, the proposed method possesses a good accuracy and efficiency for this investigated linear system.

(a) PDF

(b) .POE

(c) PDF

(d) POE

(a) PDF

(b) POE

(c) PDF

(d) POE

We completed all calculations on a desktop PC with [email protected] GHz and 16 GB RAM. Table 3 reports the average computational time of the SDOF oscillator based on 30 runs of different combinations of and . Although it is hard to draw a common conclusion from these data, the average calculational time that are around 20~ 40 s is attractive for engineering optimization design.

5.2. A Nonlinear Five-DOF Structural System

The second example is a five-DOF shear building adapted from Refs. [34, 51]. The structural system shown in Figure 6 is subjected to an external excitation. There are 25 input random variables to model uncertainties in the system properties, including masses, stiffness, sliding force stiffness ratios, and damping ratios. Table 4 lists the means of these 25 input random variables which are independently and normally distributed. Note that the damping terms in are given by where are the damping ratios, are the masses, and are the stiffness, respectively. Except for the damping variables, the coefficient of variation of the other random property variables is set to be 0.1. The standard deviations of damping variables are set to be 0.01.

Suppose that the external excitation to the structural system is a general Gaussian excitation acting on the horizontal ground. The Karhunen-Loève expansion is used to represent it with 200 independent standard normal variables. Therefore, there are total 225 input random variables which are beyond the capacity of FORM and surrogate model methods.

The failure event is defined by the first-passage of the maximum displacement at the first story with a threshold For the structural dynamic analysis, the time interval of is discretized by a time step size of 0.01 s.

For this problem, 500 experimental points generated from MCS are employed to estimate the fractional moments of the normalized response. The comparisons of fractional moments and the corresponding relative errors are shown in Figure 7. Similar conclusions can be drawn as for the SDOF oscillator, except that the magnitude of relative error grows to 2.75% for .

(a) Fractional moment

(b) Relative error

Additionally, , 4, and 5 are adopted in the proposed method. Table 5 lists the values of the objective function with penalty term. As the value of increases, the values of the objective function also increase, which is consistent with the results of the SDOF oscillator. This indicates that the MaxEnt distribution with has a good balance of accuracy and model complexity. Table 6 lists the orders of fractional moments and Lagrange multipliers obtained by the proposed method.

Based on the parameters in Table 6, we reconstructed the distribution of the normalized response and compared the proposed method and MCS in Figure 8. The PDF curve of the proposed method matches the histogram of the experimental data closely (Figure 8(a)), while we observed a dispersion between it and the histogram of MCS samples in Figure 8(c). Similarly, the POE curve of the proposed method agrees with that of experimental data (Figure 8(b)), whereas it has a deviation from the curve estimated from MCS samples (Figure 8(d)). For the interval of , the POE of the proposed method has a good agreement with that of MCS. The failure probability is estimated as by MCS with 10⁷ samples when the threshold level is 0.039 m (corresponding to the case of the normalized threshold level of 1). The estimation of failure probability yielded by the proposed method is . Therefore, the proposed method overestimated the failure probability for a large threshold level. Compared with the investigated linear system, the overestimations may cause by the nonlinear behavior of the current structural system.

(a) PDF

(b) POE

(c) PDF

(d) POE

Table 7 reports the average computational time of the nonlinear structural system based on 30 runs of different . Similar conclusion can be drawn from this table as that of the SDOF oscillator example.

6. Conclusions

In this study, a maximum entropy distribution method is proposed to perform structural reliability analysis for high-dimensional problems with random process input. With the development of the maximum entropy distribution with fractional moments, we proved that the solution to a function is the global one to the maximum relative entropy between the true PDF and the approximate one. Thus, the corresponding optimization problem is converted into a convex one with a global solution. To estimate the fractional moments efficiently, a set of experimental points are generated for this purpose, which provides the estimator of a fractional moment for any order. Converting the common optimization problem into an equivalent one can avoid the trapping in a locally optimal solution. Furthermore, we reexamined the associated computational issues of the target optimization problem and reformulated it with side constraints on the optimization variables, i.e., the orders of fractional moments and the corresponding Lagrange multipliers. Subset simulation optimization is used to solve the new optimization problem.

Two benchmark problems with high-dimensional inputs show that the proposed method has both good accuracy and efficiency. If the limit state function, however, has a strong nonlinear behavior with respect to inputs, the maximum entropy distribution may not approximate the tail region accurately, and in this case, the proposed method may result in a low efficiency or accuracy, or both.

Future studies will involve developing efficiently sampling methods for the estimation of fractional moments for high-dimensional nonlinear problems.

Appendix

A. Symbols

:limit state function

, :response quantity and its value given an input

, :input random variables and its observation

,:random vector and its observation vector

:input random process vector

,:time and final time

:reliability function

:failure probability function

:time step

:random vector transformed from input random process

:probability operator

:expectation operator

:maximum value variable

, :order of fractional moment

,:probability density function and approximate probability density function

:reference point

:th fractional moment

:Lagrange multiplier

:objective function

:potential function

and submatrix

B.

Recall that the PDF of MaxEnt distribution is given by

Since , one could obtain the following relationship:

Substitute Equation (B.2) into Equation (B.1), the PDF expression of MaxEnt distribution is

Further, the PDF expression of MaxEnt distribution can be rewritten as as long as . With the definition of in Equation (17), one obtain that

Data Availability

Some or all data, models, or code generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

R. Rackwitz and B. Flessler, “Structural reliability under combined random load sequences,” Computers & Structures, vol. 9, no. 5, pp. 489–494, 1978.
View at: Publisher Site | Google Scholar
K. Breitung, “Asymptotic approximations for probability integrals,” Probabilistic Engineering Mechanics, vol. 4, no. 4, pp. 187–190, 1989.
View at: Publisher Site | Google Scholar
X. Du and Z. Hu, “First order reliability method with truncated random variables,” Journal of Mechanical Design, vol. 134, no. 9, article 091005, 2012.
View at: Publisher Site | Google Scholar
Y.-G. Zhao and T. Ono, “Moment methods for structural reliability,” Structural Safety, vol. 23, no. 1, pp. 47–75, 2001.
View at: Publisher Site | Google Scholar
Y.-G. Zhao and T. Ono, “On the problems of the fourth moment method,” Structural Safety, vol. 26, no. 3, pp. 343–347, 2004.
View at: Publisher Site | Google Scholar
C. G. Bucher and U. Bourgund, “A fast and efficient response surface approach for structural reliability problems,” Structural Safety, vol. 7, no. 1, pp. 57–66, 1990.
View at: Publisher Site | Google Scholar
M. R. Rajashekhar and B. R. Ellingwood, “A new look at the response surface approach for reliability analysis,” Structural Safety, vol. 12, no. 3, pp. 205–220, 1993.
View at: Publisher Site | Google Scholar
H. P. Gavin and S. C. Yau, “High-order limit state functions in the response surface method for structural reliability analysis,” Structural Safety, vol. 30, no. 2, pp. 162–179, 2008.
View at: Publisher Site | Google Scholar
H. S. Li, Z. Z. Lu, and H. W. Qiao, “A new high-order response surface method for structural reliability analysis,” Structural Engineering and Mechanics, vol. 34, no. 6, pp. 779–799, 2010.
View at: Publisher Site | Google Scholar
R. E. Melchers, “Importance sampling in structural systems,” Structural Safety, vol. 6, no. 1, pp. 3–10, 1989.
View at: Publisher Site | Google Scholar
S. K. Au and J. L. Beck, “A new adaptive importance sampling scheme for reliability calculations,” Structural Safety, vol. 21, no. 2, pp. 135–158, 1999.
View at: Publisher Site | Google Scholar
P. H. Madsen and S. Krenk, “An integral equation method for the first-passage problem in random vibration,” Journal of Applied Mechanics, vol. 51, no. 3, pp. 674–679, 1984.
View at: Publisher Site | Google Scholar
Ø. Hagen and L. Tvedt, “Vector process out‐crossing as parallel system sensitivity measure,” Journal of Engineering Mechanics, vol. 117, no. 10, pp. 2201–2220, 1991.
View at: Publisher Site | Google Scholar
C. Andrieu-Renaud, B. Sudret, and M. Lemaire, “The PHI2 method: a way to compute time-variant reliability,” Reliability Engineering & System Safety, vol. 84, no. 1, pp. 75–86, 2004.
View at: Publisher Site | Google Scholar
B. Sudret, “Analytical derivation of the outcrossing rate in time-variant reliability problems,” Structure and Infrastructure Engineering, vol. 4, no. 5, pp. 353–362, 2008.
View at: Publisher Site | Google Scholar
J. Zhang and X. Du, “Time-dependent reliability analysis for function generator mechanisms,” Journal of Mechanical Design, vol. 133, no. 3, article 031005, 2011.
View at: Publisher Site | Google Scholar
Z. Hu and X. Du, “Time-dependent reliability analysis with joint upcrossing rates,” Structural and Multidisciplinary Optimization, vol. 48, no. 5, pp. 893–907, 2013.
View at: Publisher Site | Google Scholar
Z. Hu, H. Li, X. du, and K. Chandrashekhara, “Simulation-based time-dependent reliability analysis for composite hydrokinetic turbine blades,” Structural and Multidisciplinary Optimization, vol. 47, no. 5, pp. 765–781, 2013.
View at: Publisher Site | Google Scholar
J.-B. Chen and J. Li, “The extreme value distribution and dynamic reliability analysis of nonlinear structures with uncertain parameters,” Structural Safety, vol. 29, no. 2, pp. 77–93, 2007.
View at: Publisher Site | Google Scholar
J. Li, J.-B. Chen, and W.-L. Fan, “The equivalent extreme-value event and evaluation of the structural system reliability,” Structural Safety, vol. 29, no. 2, pp. 112–131, 2007.
View at: Publisher Site | Google Scholar
Z. Wang and P. Wang, “A nested extreme response surface approach for time-dependent reliability-based design optimization,” Journal of Mechanical Design, vol. 134, no. 12, article 121007, 2012.
View at: Publisher Site | Google Scholar
Z. Hu and X. Du, “A sampling approach to extreme value distribution for time-dependent reliability analysis,” Journal of Mechanical Design, vol. 135, no. 7, article 071003, 2013.
View at: Publisher Site | Google Scholar
Z. Wang and P. Wang, “A double-loop adaptive sampling approach for sensitivity-free dynamic reliability analysis,” Reliability Engineering & System Safety, vol. 142, pp. 346–356, 2015.
View at: Publisher Site | Google Scholar
Z. Hu and S. Mahadevan, “A single-loop kriging surrogate modeling for time-dependent reliability analysis,” Journal of Mechanical Design, vol. 138, no. 6, article 061406, 2016.
View at: Publisher Site | Google Scholar
Z. Wang and W. Chen, “Confidence-based adaptive extreme response surface for time-variant reliability analysis under random excitation,” Structural Safety, vol. 64, pp. 76–86, 2017.
View at: Publisher Site | Google Scholar
Z. Hu and X. Du, “Mixed efficient global optimization for time-dependent reliability analysis,” Journal of Mechanical Design, vol. 137, no. 5, article 051401, 2015.
View at: Publisher Site | Google Scholar
D. Zhang, X. Han, C. Jiang, J. Liu, and Q. Li, “Time-dependent reliability analysis through response surface method,” Journal of Mechanical Design, vol. 139, no. 4, article 041404, 2017.
View at: Publisher Site | Google Scholar
A. Singh, Z. P. Mourelatos, and J. Li, “Design for lifecycle cost using time-dependent reliability,” Journal of Mechanical Design, vol. 132, no. 9, article 091008, 2010.
View at: Publisher Site | Google Scholar
Z. P. Mourelatos, M. Majcher, V. Pandey, and I. Baseski, “Time-dependent reliability analysis using the total probability theorem,” Journal of Mechanical Design, vol. 137, no. 3, article 031405, 2015.
View at: Publisher Site | Google Scholar
G. J. Savage and Y. K. Son, “Dependability-based design optimization of degrading engineering systems,” Journal of Mechanical Design, vol. 131, no. 1, article 011002, 2009.
View at: Publisher Site | Google Scholar
Y. K. Son and G. J. Savage, “Set theoretic formulation of performance reliability of multiple response time-variant systems due to degradations in system components,” Quality and Reliability Engineering International, vol. 23, no. 2, pp. 171–188, 2007.
View at: Publisher Site | Google Scholar
C. Jiang, X. P. Huang, X. Han, and D. Q. Zhang, “A time-variant reliability analysis method based on stochastic process discretization,” Journal of Mechanical Design, vol. 136, no. 9, article 091009, 2014.
View at: Publisher Site | Google Scholar
S.-K. Au and J. L. Beck, “Estimation of small failure probabilities in high dimensions by subset simulation,” Probabilistic Engineering Mechanics, vol. 16, no. 4, pp. 263–277, 2001.
View at: Publisher Site | Google Scholar
S. K. Au, J. Ching, and J. L. Beck, “Application of subset simulation methods to reliability benchmark problems,” Structural Safety, vol. 29, no. 3, pp. 183–193, 2007.
View at: Publisher Site | Google Scholar
M. G. Kendall, A. Stuart, and J. K. Ord, Kendall's Advanced Theory of Statistics, Charles Griffin & Company Ltd, vol. 1, 1987.
X. Zhang and M. D. Pandey, “Structural reliability analysis based on the concepts of entropy, fractional moment and dimensional reduction method,” Structural Safety, vol. 43, pp. 28–40, 2013.
View at: Publisher Site | Google Scholar
G. Li and K. Zhang, “A combined reliability analysis approach with dimension reduction method and maximum entropy method,” Structural and Multidisciplinary Optimization, vol. 43, no. 1, pp. 121–134, 2011.
View at: Publisher Site | Google Scholar
E. T. Jaynes, “Information theory and statistical mechanics,” Physical Review, vol. 106, no. 4, pp. 620–630, 1957.
View at: Publisher Site | Google Scholar
E. T. Jaynes, Probability Theory: The Logic of Science, Cambridge Univeristy Press, Cambridge, UK, 2003.
S. Rahman and H. Xu, “A univariate dimension-reduction method for multi-dimensional integration in stochastic mechanics,” Probabilistic Engineering Mechanics, vol. 19, no. 4, pp. 393–408, 2004.
View at: Publisher Site | Google Scholar
H. Xu and S. Rahman, “A generalized dimension-reduction method for multidimensional integration in stochastic mechanics,” International Journal for Numerical Methods in Engineering, vol. 61, no. 12, pp. 1992–2019, 2004.
View at: Publisher Site | Google Scholar
P. L. N. Inverardi and A. Tagliani, “Maximum entropy density estimation from fractional moments,” Communications in Statistics - Theory and Methods, vol. 32, no. 2, pp. 327–345, 2003.
View at: Publisher Site | Google Scholar
B. Li, L. Zhang, X. Zhu, X. Yu, and X. Ma, “Reliability analysis based on a novel density estimation method for structures with correlations,” Chinese Journal of Aeronautics, vol. 30, no. 3, pp. 1021–1030, 2017.
View at: Publisher Site | Google Scholar
H. Dai, H. Zhang, and W. Wang, “A new maximum entropy-based importance sampling for reliability analysis,” Structural Safety, vol. 63, pp. 71–80, 2016.
View at: Publisher Site | Google Scholar
Y. Shi, Z. Lu, K. Cheng, and Y. Zhou, “Temporal and spatial multi-parameter dynamic reliability and global reliability sensitivity analysis based on the extreme value moments,” Structural and Multidisciplinary Optimization, vol. 56, no. 1, pp. 117–129, 2017.
View at: Publisher Site | Google Scholar
J. Xu, “A new method for reliability assessment of structural dynamic systems with random parameters,” Structural Safety, vol. 60, pp. 130–143, 2016.
View at: Publisher Site | Google Scholar
J. Xu, C. Dang, and F. Kong, “Efficient reliability analysis of structures with the rotational quasi- symmetric point- and the maximum entropy methods,” Mechanical Systems and Signal Processing, vol. 95, pp. 58–76, 2017.
View at: Publisher Site | Google Scholar
J. Xu and D. Wang, “A two-step methodology to apply low-discrepancy sequences in reliability assessment of structural dynamic systems,” Structural and Multidisciplinary Optimization, vol. 57, no. 4, pp. 1643–1662, 2018.
View at: Publisher Site | Google Scholar
J. Xu, D. Wang, and C. Dang, “A marginal fractional moments based strategy for points selection in seismic response analysis of nonlinear structures with uncertain parameters,” Journal of Sound and Vibration, vol. 387, pp. 226–238, 2017.
View at: Publisher Site | Google Scholar
J. Xu, W. Zhang, and R. Sun, “Efficient reliability assessment of structural dynamic systems with unequal weighted quasi-Monte Carlo simulation,” Computers and Structures, vol. 175, pp. 37–51, 2016.
View at: Publisher Site | Google Scholar
G. I. Schuëller and H. J. Pradlwarter, “Benchmark study on reliability estimation in higher dimensions of structural systems - an overview,” Structural Safety, vol. 29, no. 3, pp. 167–182, 2007.
View at: Publisher Site | Google Scholar
Z. Li, F. Nie, X. Chang, L. Nie, H. Zhang, and Y. Yang, “Rank-constrained spectral clustering with flexible embedding,” IEEE Transactions on Neural Networks and Learning Systems, vol. 29, no. 12, pp. 6073–6082, 2018.
View at: Publisher Site | Google Scholar
Z. Li, F. Nie, X. Chang, Y. Yang, C. Zhang, and N. Sebe, “Dynamic affinity graph construction for spectral clustering using multiple features,” IEEE Transactions on Neural Networks and Learning Systems, vol. 29, no. 12, pp. 6323–6332, 2018.
View at: Publisher Site | Google Scholar
Z. Li, L. Yao, X. Chang, K. Zhan, J. Sun, and H. Zhang, “Zero-shot event detection via event-adaptive concept relevance mining,” Pattern Recognition, vol. 88, pp. 595–603, 2019.
View at: Publisher Site | Google Scholar
H. Gzyl and A. Tagliani, “Hausdorff moment problem and fractional moments,” Applied Mathematics and Computation, vol. 216, no. 11, pp. 3319–3328, 2010.
View at: Publisher Site | Google Scholar
H.-S. Li, “Subset simulation for unconstrained global optimization,” Applied Mathematical Modelling, vol. 35, no. 10, pp. 5108–5120, 2011.
View at: Publisher Site | Google Scholar
H.-S. Li and Z.-J. Cao, “Matlab codes of subset simulation for reliability analysis and structural optimization,” Structural and Multidisciplinary Optimization, vol. 54, no. 2, pp. 391–410, 2016.
View at: Publisher Site | Google Scholar
H.-S. Li and S.-K. Au, “Design optimization using subset simulation algorithm,” Structural Safety, vol. 32, no. 6, pp. 384–392, 2010.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Fuliang Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

269

Downloads

614

Citations