Abstract
Considering the uncertainties of seismic excitation and stochastic parameters of the Interstory Isolation Structure, the implicit structural response function is transformed into an explicit one by the Stochastic Response Surface Method (SRSM), combined with a probabilistic collocation method based on the linear independence principle. By using the reliability analysis method based on the firstpassage failure, the reliability of isolation structure is analyzed and the sensitivity of the reliability of each substructure is further analyzed. The results show that conducting sensitivity research based on reliability analysis can provide an important theoretical basis for an optimized and robust design of a random vibration system. The changes of the stiffness and damping ratio of the Isolation Layer have a great influence on the reliability of the system and have different effects on the reliability of different parts of the system. The analysis results of reliability sensitivity could provide system reliability indexes to rank stochastic parameters by their importance, thus improving the efficiency of the structure reliability analysis and optimization design.
Highlights
 Reliability sensitivity analysis of stochastic structures is introduced into the Interstory Isolation System.
 Considering the uncertainties of seismic excitation and stochastic parameters of the Interstory Isolation Structure.
 The random response surface method and the reliability analysis method based on the first over failure are used to study the influence of random variables on reliability sensitivity.
1. Introduction
Reliability sensitivity reflects the impact of the change of design parameters on failure probability, which is often used in the reliabilitybased optimization design of structures to determine a search direction for an optimal solution and plays an important role in reliabilitybased optimization design of structures. Reliability sensitivity is defined as the ratio of the variation of the failure probability caused by change in the distribution parameters of basic random variables [13]. Mathematically, it is the partial derivative of the failure probability to the distribution parameters of basic random variables. In multivariate reliability analysis, Hohenbichler and Rackwitz [4] first put forward the traditional concepts of sensitivity and importance. Karamchandani and Cornell [5] analyzed the reliability sensitivity of structures of distributed parameters based on firstorder and secondorder reliability analysis. Wu [6] put forward an effective adaptive importance sampling method to analyze the reliability and reliability sensitivity of structures. Robert et al. [7] put forward the MPPSS (Most Probable Point System Simulation), which is suitable for calculating series systems, parallel systems, and seriesparallel hybrid systems and can be used to determine the system sensitivity factor, that is, the importance of each random variable to the system reliability. Lu et al. [8] used the point estimation method and highorder moment standardized technology to analyze the sensitivity of the failure probability to the change of randomly distributed parameters. Due to the low efficiency of the Monte Carlo method, Wu et al. [9] put forward the adaptive importance sampling method for calculating structure reliability and reliability sensitivity, in which sampling density in proportion to the joint probability density of random variables was used to improve the sampling efficiency. Lu et al. [10] put forward the reliability sensitivity analysis method for subset sampling and line sampling.
For the description of uncertain factors, Guo et al. [11] analyzed the sensitivity of structure parameters under the variables of mixed types, i.e., random variables and interval variables. Song et al. [12] developed a method based on a matrix to transform the solution of system reliability into matrix operation and analyzed the sensitivity of each distributed parameter under system reliability on this basis. Zhang et al. [1315] dealt with the reliability sensitivity problem of single degreeoffreedom nonlinear vibration systems with random parameters and put forward a numerical method for reliability sensitivity analysis of mechanical parts with multiple failure modes. Yuan et al. [16] developed two reliability sensitivity analysis methods under the nonnormal random variables for nonlinear limit state equations. Valdebenito et al. [17] analyzed linear systems under random excitations in depth and put forward a new method for analyzing reliability and reliability sensitivity of linear systems. Song et al. [18] developed fuzzy reliability sensitivity by using the moment method on the basis of the fuzzy theory. Qiao et al. [19, 20] established the reliability sensitivity analysis method based on the stochastic response surface method and applied it to the reliabilitybased optimization design of structures. In order to quantitatively evaluate the contribution of various uncertainties to the uncertainty of windinduced response of structures. Wu et al. [21] conducted parametric sensitivity studies on windinduced response of typical longspan roof structures. Sun et al. [22] took the stiffness ratio, the postyield stiffness coefficient, the yield displacement of the isolation layer, the friction coefficient, the mass ratio and the natural period of the superstructure as the influencing factors, sensitivity analysis of the total earthquake input energy is carried out for rural building structures supported on sliding baseisolation systems. Dong et al. [23] deduced the recursion formula of honeycombtype cable dome with two struts under prestress distribution state, and studied the distribution characteristic and rule of prestress distribution state through parameter analysis, and deeply analyzed the sensitivity of each parameter. Khaneghahi et al. [24] used global sensitivity analysis based on improved Latin hypercube sampling to indicate the influence of each random variable and their interaction on variance of the responses, which could show how each random variable and its interactions affect the response variance.
As a new type of structural system, the Interstory Isolation System shows the different working mechanisms and main parameters affecting structural reliability due to the different locations of the Isolation Layer. Therefore, it is of great significance to study the main parameters affecting the structural reliability due to the change of the Isolation Layer location for the economy and safety of the Interstory Isolation System. In this paper, reliability sensitivity analysis of stochastic structures is introduced into the Interstory Isolation System. Considering the uncertainty of seismic excitation and random parameters, the random response surface method and the reliability analysis method based on the first over failure are used to study the influence of random variables on reliability sensitivity. Through the sensitivity analysis of structural reliability, it can be determined which design parameter is the most effective to adjust, so as to select as few optimization variables as possible, to help engineers and technicians reasonably establish the safety tolerance of the structural system and control the impact of random parameters on the safety of the structural system. This analysis has guiding significance for the reliability design, modification, and optimization of the Interstory Isolation System.
2. Modeling of the Interstory isolation system
2.1. Equations of motion
Fig 1. illustrates the calculation model of a Interstory isolation system consisting of two substructures and isolation layer of ${n}_{s}$, ${n}_{u}$ and one story, respectively. For the sheartype buildings, the equation of the structural system can be expressed as:
where $\mathbf{M}$, $\mathbf{C}$ and $\mathbf{K}$ are the $N\times N$dimensional ($N={n}_{s}+{n}_{u}+1$) mass, damping, and stiffness matrices of the system, respectively. $\mathbf{u}$ is the $N$dimensional relative displacement vector with respect to the ground and a dot over the symbol indicates differentiation with respect to time. $\mathbf{I}$ is the load index vector with all its elements equal to 1. ${\ddot{u}}_{g}$ is the ground motion acceleration.
Fig. 1Computational model of the Interstory isolation structure
In order to derive the covariance matrix equation, the equation of motion can be transformed to the statespace form:
where:
This equation can be solved by Complex mode method, which is much easier. It is unnecessary to go into details. Given the transfer function matrix such that [25]:
where ${\lambda}_{r}$ and ${\psi}_{r}$ are the complex eigenvalue and complex exponential, respectively, which is performed by Complex modal analysis method, $j$ is the imaginary unit $\sqrt{1}$. The sign “*” represents the complex conjugate and the superscript “$T$” denotes the transpose of a matrix.
2.2. Seismic response under filtered white noise excitation
The CloughPenzien filtered white noise ground motion model, which has been widely used in earthquake engineering, can be represented by the following spectrum [26]:
where ${\omega}_{g}$, ${\omega}_{f}$, ${\xi}_{g}$, ${\xi}_{f}$ and ${S}_{0}$ may be regarded as the characteristics parameters of the soil surrounding.
The autopower crosspower spectral density function of the displacement response of the system under the ground acceleration can be determined by [27]:
The displacement response standard deviation of each floor can be evaluated through the numerical integration of the autopower crosspower spectral density function:
And its standard deviation of acceleration response is given by:
3. Sensitivity analysis of structural reliability
3.1. Dynamic reliability based on the Firstpassage Failure
Dynamic reliability analysis of structures subjected to earthquake action is a typical analysis based on the Firstpassage Failure [27]. Assume that both seismic excitation and structural response are stationary random processes with zero mean. Recording ${y}_{e}$ and ${\sigma}_{y}$ are extremes (absolute maximum) and standard deviations of structural response $y\left(t\right)$, respectively, and the dimensionless parameters are defined here:
According to the assumption that the number of horizontal crossings is a Poisson process, the expected value of the extreme value of $\eta $ is approximately equal to:
The standard deviation of the extreme value of $\eta $ is:
The expectation of the extreme value of the stochastic process $y\left(t\right)$ can be expressed as:
The state equation of the Interstory Isolation System expressed by the limit state of the interstory displacement angle is as follows:
where ${y}_{max}$ is the bounded interlayer displacement angle, $\stackrel{}{y}\left(x\right)$ is the extreme mean of the interlayer displacement angle, and $x$ is the vector of random variables.
Limit state function: $G\left(x\right)\le 0$ is failure state. $G\left(x\right)>0$ is a safety state.
The reliability index is defined as:
If $y\left(x\right)$ obeys a normal distribution, reliability can be obtained:
3.2. Stochastic response surface model
The Stochastic Response Surface Method (SRSM) [28] is based on the assumption that random variables with square integrable probability density function can be expressed as functional forms of independent standard distribution random variables. SRSM converts input random variables into functions of independent and identically distributed standard normal random variables, and uses chaotic polynomials of independent standard normal random variables to describe random responses.
The input random parameters of the structure affect the output response of the structure. If the input random parameters are represented by independent standard normal random variables ${\left\{{\xi}_{i}\right\}}_{i=1}^{n}$, the output model of the structure can be expressed as a chaotic polynomial model with only standard normal distribution random variables [29]:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{{i}_{1}=1}^{n}\sum _{{i}_{2}=1}^{{i}_{1}}\sum _{{i}_{3}=1}^{{i}_{2}}{b}_{{i}_{1}{i}_{2}{i}_{3}}{\mathrm{\Gamma}}_{3}\left({\xi}_{{i}_{1}},{\xi}_{{i}_{2}},{\xi}_{{i}_{3}}\right)+\cdots ,$
where $y$ represents the structural response variable, $n$ is the number of random variables in the input model, ${\xi}_{{i}_{j}}$ is the random variable of $j$ standard normal distribution, and ${\mathrm{\Gamma}}_{p}\left({\xi}_{{i}_{1}}\cdots {\xi}_{{i}_{p}}\right)$ is a $p$order $n$dimensional Hermite polynomial:
The results of the stochastic response surface method depend on the probabilistic collocation selection of the parameter ${b}_{{i}_{1}}$, ${b}_{{i}_{1}{i}_{2}}$,…, ${b}_{{i}_{1}{i}_{2}\cdots {i}_{n}}$ to be estimated. Generally speaking, the expansion collocation of porder random polynomials is determined by the root of ($p+1$)order Hermite polynomials. According to this method, the number of collocation points available is:
For $n$ random variables, the number of expansion terms of $p$order Hermite polynomials (undetermined coefficient) ${N}_{b}$ is:
Generally speaking, $N\gg {N}_{b}$, so choosing collocation ${N}_{b}$ is only part of $N$. In this paper, the probability collocation points are selected based on the principle of linear independence, and a large number of linear correlation collocations are eliminated to ensure that the coefficient matrix of chaotic polynomials is linearindependent in rows. Moreover, the coefficient matrix is a full rank matrix, and the undetermined coefficient has a unique solution. This not only ensures the stability of the calculation results, but also significantly improves the calculation efficiency.
The above sample points are all in the standard normal space and are converted to the original space according to need. Variables selected in this paper are assumed to obey normal distribution, and the transformation form is as follows:
where ${\mu}_{i}$ and ${\sigma}_{i}$ are the mean and standard deviation of the variable ${x}_{ip}$, respectively. The original space sample points are transformed to form test samples, and the real response of the structure is obtained by a dynamic analysis method. The undetermined coefficients can be obtained by solving the linear equations with the least squares method, i.e., the explicit expression function of the stochastic response surface model.
To ensure the adaptability of the stochastic response surface model, the multiple fitting coefficient ${R}^{2}$ and the modified multiple fitting coefficient ${R}_{A}^{2}$ are used to test adaptability, and the coefficients are defined as follows:
where $\mathbf{B}$ is the coefficient matrix of the parameter to be estimated, $\mathbf{X}$ is the Hermite polynomial combination matrix, $\mathbf{Y}$ is the response matrix, $k$ is the number of samples for model validation, and $p$ is the number of random variables.
3.3. Sensitivity analysis of reliability
The reliability sensitivity reflects the influence of changing the design parameters on the failure probability. It helps to determine the search direction of the optimal solution of structural reliability optimization design, which is an important part of structural reliability optimization design [30]. In our paper, a reliability sensitivity analysis method based on Monte Carlo simulation was established on the explicit function generated by SRSM. The sensitivity of the mean matrix $\mu $ and the mean square deviation matrix $\mathbf{D}$ of the basic random parameter vector is expressed as follows:
where:
$\frac{\partial {\mu}_{y}}{\partial {\mu}^{T}}=\left(\frac{\partial {\mu}_{y}}{\partial {\mu}_{1}},\frac{\partial {\mu}_{y}}{\partial {\mu}_{2}},...,\frac{\partial {\mu}_{y}}{\partial {\mu}_{n}}\right),\frac{\partial {\mu}_{y}}{\partial {\mathbf{D}}^{T}}=\left(\frac{\partial {\mu}_{y}}{\partial {D}_{1}},\frac{\partial {\mu}_{y}}{\partial {D}_{2}},...,\frac{\partial {\mu}_{y}}{\partial {D}_{n}}\right),$
$\frac{\partial {D}_{y}}{\partial {\mu}^{T}}=\left(\frac{\partial {D}_{y}}{\partial {\mu}_{1}},\frac{\partial {D}_{y}}{\partial {\mu}_{2}},...,\frac{\partial {D}_{y}}{\partial {\mu}_{n}}\right),\frac{\partial {D}_{y}}{\partial {\mathbf{D}}^{T}}=\left(\frac{\partial {D}_{y}}{\partial {D}_{1}},\frac{\partial {D}_{y}}{\partial {D}_{2}},...,\frac{\partial {D}_{y}}{\partial {D}_{n}}\right).$
4. Analytical results of Interstory isolation structure
4.1. Structural model parameters
Taking the minimum base shear force as the optimization goal, superstructure and substructure optimal frequency ratio ${R}_{\omega}=$0.165 and Isolation layer optimal damping ratio $\xi =$0.25 are obtained by the NSGAII optimization method [31], in mass ratios ${R}_{m}=$2.0. The average value of mass, the shear stiffness and damping ratios of the 12story to 6story are uniform for each story. The mass, the shear stiffness and damping ratios is 2.0×10^{6 }kg, 4.0×10^{9} N/m and 0.05. The Isolation layer average value of mass and the shear stiffness is 2.0×10^{6 }kg and 7.93×10^{7} N/m. Assume that the random variables of the structure are the Isolation Layer stiffness, ${k}_{iso}$, the Isolation Layer damping ratio, ${\xi}_{iso}$, the superstructure and substructure stiffness, ${k}_{sup}$ and ${k}_{sub}$, and the superstructure and substructure mass, ${m}_{sup}$ and ${m}_{sub}$, and the coefficient of variation are 0.10, 0.12, 0.10, 0.10, 0.08, 0.08, respectively. It is also assumed that all the above random variables obey normal distribution. The heights of the superstructure and substructure are 3.0 m and 2.9 m, respectively.
The CloughPenzien filtered whitenoise spectrum is used as the ground acceleration spectrum in the computation. The characteristic ground frequency, the characteristic damping ration, and the intensity of an earthquake in the ground acceleration spectrum are selected as ${S}_{0}=$ 9.537 cm^{2}/s, ${\omega}_{f}=$ 17.95 rad/s, ${\xi}_{f}=$ 0.72, ${\omega}_{c}=$ 0.15${\omega}_{f}$, ${\xi}_{c}={\xi}_{f}$.
4.2. Analysis of structural model examples
Using SRSM, the input variables ${k}_{iso}$, ${\xi}_{iso}$, ${k}_{sup}$, ${k}_{sub}$, ${m}_{sup}$, ${m}_{sub}$ are converted into functions of independent and identically distributed standard normal random variables. The interlayer displacement angle, the output variable of the structure, is expressed as a secondorder SRSM. The corresponding collocation uses the roots of the thirdorder onedimensional Hermite polynomial.
The onedimensional thirdorder Hermite polynomial is ${\xi}^{3}3\xi $, and the corresponding roots are: –$\sqrt{\text{3}}$, 0, $\sqrt{\text{3}}$. The output variable is a thirdorder SRSM:
According to the example parameters, there are 6 random variables. Using the thirdorder Hermite polynomial, we can see that the number of collocations $N$ to be selected is 729 and the number of undetermined coefficients ${N}_{b}$ in the SRSM expansion formula is 28, according to Eq. (18). Selecting the probability collocation based on the linear independent principle can ensure the stability of the calculation results and improve the calculation efficiency.
The random variables are $\xi =\{{k}_{iso},{\xi}_{iso},{k}_{sub},{k}_{sup},{m}_{sub},{m}_{sup}\}$, and according to the substructure responses, the Stochastic Response Surface Function coefficients are expressed as follows:
(1) The Stochastic Response Surface Function coefficients of the substructure response:
${b}_{0}=${–0.073};
${b}_{ii}=${7.164×10^{11}, 0.053, 2.728×10^{13}, –1.414×10^{12}, 1.129×10^{9}, 2.091×10^{9}};
${b}_{ij}=${1.017×10^{20}, –0.066, –1.453×10^{22}, 1.084×10^{22}, 1.664×10^{17}, 4.303×10^{17}, –5.044×10^{10}, 6.712×10^{21}, 1.004×10^{22}, 3.699×10^{18}, –1.236×10^{18}, 1.481×10^{12}, 1.315×10^{12}, –4.229×10^{11}, –3.353×10^{10}, 1.244×10^{22}, –1.235×10^{19}, –2.819×10^{20}, –4.101×10^{21}, –2.253×10^{19}, –3.947×10^{16}}.
(2) The Stochastic Response Surface Function coefficients of the Isolation Layer response:
${b}_{0}=${–13.183};
${b}_{ii}=${5.715×10^{9}, 9.427, –1.956×10^{11}, –1.624×10^{10}, 3.889×10^{7}, –1.080×10^{8}};
${b}_{ij}=${3.374×10^{17}, –12.412, 1.353×10^{20}, 7.319×10^{21}, 2.097×10^{15}, 1.575×10^{14}, –8.252×10^{8}, 7.411×10^{19}, 8.088×10^{19}, –1.745×10^{15}, –2.776×10^{17}, 1.664×10^{10}, 2.585×10^{10}, –1.518×10^{7}, –2.188×10^{9}, 1.190×10^{20}, –1.039×10^{17}, 1.024×10^{17}, –6.397×10^{18}, –1.141×10^{17}, –2.418×10^{14}}.
(3) The Stochastic Response Surface Function coefficients of the superstructure response:
${b}_{0}=${0.077};
${b}_{ii}=${–5.640×10^{11}, –0.054, –2.318×10^{13}, 8.207×10^{13}, 2.047×10^{10}, 3.762×10^{10}};
${b}_{ij}=${–1.919×10^{20}, 0.070, 1.078×10^{22}, –3.230×10^{23}, –7.282×10^{17}, –5.954×10^{17}, 4.279×10^{10}, –4.245×10^{21}, –4.491×10^{21}, 3.995×10^{18}, 4.750×10^{19}, –5.280×10^{13}, –1.499×10^{12}, 3.978×10^{10}, –4.060×10^{11}, –5.620×10^{23}, –7.592×10^{20}, –4.855×10^{20}, 9.520×10^{20}, 7.051×10^{21}, 1.330×10^{17}}.
According to the Eqs. (20) and (21), the adaptability test of the function fitted by the SRSM was conducted. The multifit coefficient ${R}^{2}$ and the modified multiple fitting coefficient ${R}_{A}^{2}$ can be gotten by testing the interlayer displacement angle, the mean and the mean square deviation of the displacement of the Isolation Layer. The results are as shown in Fig. 2. Using the Monte Carlo method to sample and simulate 1.0×10^{6} times, the probability and statistical parameters of the state function of the Interstory Isolation System can be calculated, as shown in Table 1. Combined with the fitting coefficients and probability and statistical parameters, it can be seen that the established chaotic polynomial response surface function has a good fitting effect, and the SRSM has high precision for simulating random structures.
Fig. 2Multifit coefficients of each substructure response extreme value
Table 1Comparison of probability and statistical information
Number  Structuring form  Mean of limit state  Mean variance of limit state  
SRSM  Monte Carlo Method  SRSM  Monte Carlo Method  
1  Superstructure  6.3146×10^{4}  6.3139×10^{4}  3.9315×10^{4}  3.9216×10^{4} 
2  Isolation layer  0.0463  0.0463  0.0758  0.0872 
3  Substructure  6.1340×10^{4}  6.1341×10^{4}  3.5214×10^{4}  3.5045×10^{4} 
4.3. Reliability analysis
The dynamic reliability was obtained by using the reliability analysis method based on the firstpassage failure in Section 3.1, as shown in Fig. 3, wherein the interlayer displacement angle limit of the superstructure was taken as ${d}_{max}=$ 1/550 rad, and the interlayer displacement angle limit of the substructure was taken as ${d}_{max}=$ 1/450 rad. The horizontal displacement of the Isolation Layer takes the maximum value in the value 0.55 times the effective diameter of the isolation bearing and the value 3.0 times the total thickness of the bearing rubber.
Fig. 3Reliability coefficient of each layer
It can be seen from Fig. 3 that the reliability of the first layer of the substructure is the lowest, and its value is 44.35 %. As the number of layers increases, the reliability increases, with the reliability of the Isolation Layer and the superstructure 99.34 % and 100 %, respectively.
4.4. Sensitivity analysis of reliability
Assuming that the random parameters are independent of each other, let the random variable mean matrix and the mean square deviation matrix be:
To investigate the influence of different random parameters on the reliability of the system, the sensitivity analysis method of system reliability was established, using the chaotic polynomial of independent standard normal random variables to describe the random response, and the influence degree of different random parameters on the reliability of the system was discussed.
It can be seen from Fig. 6. that in the Isolation Structure System, due to the interaction between the substructure and the superstructure, the random earthquake action causes mutual coupling of the substructures, so that the same parameters have different effects on the reliability of each substructure. From the overall analysis results, the damping ratio and stiffness of the Isolation layer have a greater impact on the reliability of each substructure. An increase in these parameters will increase the reliability of the Isolation layer and, in turn, a decrease will reduce the reliability of the substructure and the superstructure.
For the superstructure, the degree of influence of different random parameter changes on its reliability is ordered as ${\xi}_{iso}()>{k}_{iso}()>{k}_{sup}()>{m}_{sub}()>{m}_{sup}()>{k}_{sub}()$, which is shown in Fig. 6(a). For the Isolation Layer, the degree of influence of different random parameter changes on its reliability is ordered as ${k}_{iso}\left(+\right)>{m}_{sup}\left(\right)>{k}_{sub}\left(\right)>{\xi}_{iso}\left(+\right)>{k}_{sup}\left(\right)>{m}_{sub}\left(+\right)$, which is shown in Fig. 6(b). For the substructure, the degree of influence of different random parameter changes on its reliability is ordered as ${\xi}_{iso}\left(\right)>{k}_{iso}\left(\right)>{m}_{sup}\left(+\right)>{k}_{sup}\left(+\right)>{k}_{sub}\left(\right)>{m}_{sub}\left(\right)$, which is shown in Fig. 6(c). Among them, the positive sign indicates that the system reliability increases with an increase of the mean value of the random variable, and the system tends to be safe; the negative sign indicates that the system reliability decreases with an increase of the mean value of the random variable, and the system tends to fail.
Fig. 6Sensitivity of extreme reliability and relative sensitivity of normalization
a) Substructure
b) Isolation layer
c) Superstructure
4.5. Results in time domain
4.5.1. Analytical model
To gain a better understanding of the dynamic behavior of the building, and further verify the sensitivity of random parameters on structural reliability under seismic excitation, threedimensional finite element analysis software (ETABS) was used to analyse the Interstory isolation structure. Threedimensional beam and beam elements were selected for the beams and columns, floor elements and Isolation bearing were used for the floor slabs and Rubber Isolator. Fig. 7 shows the side view of the analytical model, in which the longitudinal direction is defined as axis $X$ and the lateral direction as axis $Y$. The heights of the superstructure and substructure are 3.0 m and 2.9 m, respectively. Assume that the coefficient of variation are 0.10, 0.12, 0.10, 0.10, 0.08, 0.08 of the Isolation Layer stiffness, the Isolation Layer damping ratio, the superstructure and substructure stiffness, and the superstructure and substructure mass, respectively.
Fig. 7Analytical model of Interstory isolation structure
4.5.2. Analytical dynamic characteristics
The first four period are 2.586, 2.572, 2.205 and 0.588 s, respectively and the corresponding modes are translation in $X$ of superstructure, translation in $Y$ of superstructure, torsion of superstructure and translation in $X$ of structure, as shown in Fig. 8.
Fig. 8First four vibration modes
a) First vibration mode
b) Second vibration mode
c) Third vibration mode
d) Fourth vibration mode
4.5.3. Time history analysis
The Interstory isolation structure was analysed using the El centro (1940), Kobe (1995), and Northridge (1994) accelogram. The time history curve and timefrequency diagram actual earthquake are displayed in Fig. 9. It can be seen from the spectrogram that the spectral components of the seismic wave are selected mainly in the middle and low frequency, and the energy distribution is wider. The ground motion input was onedimensional and peak acceleration was set to be 70 gal corresponding to minor earthquake.
Fig. 9The historically recorded curve and timefrequency diagram
a) El centro
b) Kobe
c) Northridge
The time history analysis is conducted on the structural system excited by seismic acceleration records, with the peak accelerations being scaled to 70 gal. Each random parameter are increased separately by one time mean square error. Fig. 10(a) reports the average values of isolation layer displacement with the change of each parameter under the three earthquakes, each graph reports the parameters in term of the Isolation Layer stiffness (SI, □), the Isolation Layer damping ratio (DRI, ■), the superstructure stiffness (SSP, △), the substructure stiffness (SSB, ○), the superstructure mass (MSP, ▲), the substructure mass (MSB, ●). The virtual line reports the average values of isolation layer displacement under the three earthquakes. It can be seen from that the change of each parameter has different effects on the isolation layer displacement. With the increase of the stiffness and the damping ratio of the isolation layer as well as the mass of the substructure, the isolation layer displacement will be reduced, conversely, the other parameters will increase the displacement of the isolation layer and reduce the reliability of the isolation layer. These results also coincide with the law revealed in Fig. 6, but the degree of sensitivity has difference of parameters, which is mainly caused by the small number of selected seismic wave samples. Fig. 10(b) and Fig. 10(c) shows the interlayer displacement angle of the superstructure and substructure. For the substructure, the stiffness and the damping ratio of the isolation layer have larger influence on the interlayer displacement angle. With the increase of the stiffness and the damping ratio, the interlayer displacement angle will be increased and reduce the reliability of the substructure, and others have the same effect. In the same way, For the superstructure, with the increase of the stiffness and the damping ratio of the isolation layer, the interlayer displacement angle will be increased and the system tends to be dangerous, substructure and superstructure stiffness increased that could reduce interlayer displacement angle, and the system tends to be safe, other parameters have relatively little effect. These results are consistent with the sensitivity analysis of system in section 4.4.
Fig. 10Interlayer displacement angle of each part of the structure
a) The isolation layer
b) Substructure
c) Supstructure
5. Conclusions
1) By conducting sensitivity research based on the reliability analysis of the Interstory Isolation System, a quantitative relationship between the reliability (or failure probability) and the fluctuation of the distribution parameter of the basic random variable can be obtained, which can provide an important theoretical basis for an optimized and robust Interstory Isolation System design.
2) In the Interstory Isolation System, the change of the stiffness and damping ratio of the Isolation Layer has a great influence on the reliability of the system. The change of this parameter will have different effects on the reliability of different parts of the system. In the design of this system, the stiffness and damping ratio values of the Isolation Layer should be strictly optimized to reduce their impact on the reliability of each part of the system.
3) Reliability sensitivity analysis results of the Interstory Isolation Structure can provide an order of importance of structural reliability indicators for random parameter changes. For the less influential parameters, it can be used as a quantitative treatment to improve the reliability of largescale complex structural reliability analysis and optimization design.
Acknowledgements
The support of the National Natural Science Foundation of China (Grant No. 51608283) and The Project of Young Scholars in the West (Grant No. XAB2017AW09) and Natural Science Foundation of Ningxia (Grant No. 2019AAC03041) are much appreciated.
References

Helton J. C., Johnson J. D., Sallaberry C., et al. Survey of samplingbased methods for uncertainty and sensitivity analysis. Reliability Engineering and System Safety, Vol. 91, Issue 10, 2006, p. 11751209.

Au S. Reliabilitybased design sensitivity by efficient simulation. Computers and Structures, Vol. 83, Issue 14, 2005, p. 10481061.

Song S., Lu Z., Qiao H. Subset simulation for structural reliability sensitivity analysis. Reliability Engineering and System Safety, Vol. 94, Issue 2, 2009, p. 658665.

Hohenbichler M., Rackwitz R. Sensitivity and importance measures in structural reliability. Civil Engineering System, Vol. 3, Issue 4, 1986, p. 203209.

Karamchandani A., Cornell C. A. Sensitivity estimation within first and second order reliability method. Structural Safety, Vol. 11, Issue 2, 1992, p. 95107.

Wu Y.T. Computational methods for efficient structural reliability and reliability sensitivity analysis. AIAA Journal, Vol. 32, Issue 8, 1994, p. 17171723.

Robert H. S., Mark A. C. System reliability and sensitivity factors via the MPPSS method. Probabilistic Engineering Mechanics, Vol. 20, Issue 2, 2005, p. 148157.

Lu Z., Song J., Song S., et al. Reliability sensitivity by method of moments. Applied Mathematical Modelling, Vol. 34, Issue 10, 2010, p. 28602871.

Wu Y. T. Computational methods for efficient structural reliability sensitivity analysis. AIAA Journal, Vol. 32, Issue 8, 1994, p. 17171723.

Lu Z., Song S., Yue Z., et al. Reliability sensitivity method by line sampling. Structural Safety, Vol. 30, Issue 6, 2008, p. 517532.

Guo J., Du X. Reliability sensitivity analysis with random and interval variables. International Journal for Numerical Methods in Engineering, Vol. 78, Issue 13, 2009, p. 15851617.

Song J., Kang W. H. System reliability and sensitivity under statistical dependence by matrixbased system reliability method. Structural Safety, Vol. 31, Issue 2, 2009, p. 148156.

Zhang Y. M., Liu Q. L., Wen C. B. Sensitivity of reliability in nonlinear random systems with independent failure modes. Acta Mechanica Sinica, Vol. 35, Issue 1, 2003, p. 117120.

Wang X. G., Zhang Y. M., Wang B. Y. Dynamic reliability sensitivity analysis of mechanical components. Journal of Mechanical Engineering, Vol. 46, Issue 10, 2010, p. 188192.

Lu H., Zhang Y. M., Zhan C. L., et al. Reliability sensitivity estimation of mechanical components with multiple failure modes. Journal of Mechanical Engineering, Vol. 48, Issue 2, 2012, p. 6367.

Yuan X. K., Lv Z. Z. Reliability sensitivity algorithm and its application in fatigue failure. Chinese Journal of Computational Mechanics, Vol. 24, Issue 1, 2007, p. 6973.

Valdebenito M. A., Jensen H. A., Schuëller G. I., et al. Reliability sensitivity estimation of linear systems under stochastic excitation. Computers and Structures, Vol. 92, Issue 93, 2012, p. 257268.

Jun S., Zhenzhou L. Moment method based on fuzzy reliability sensitivity analysis for a degradable structural system. Chinese Journal of Aeronautics, Vol. 21, Issue 6, 2008, p. 518525.

Qiao H. W., Lv Z. Z., Zhan X. P. The topology optimization design of the continuum structure based on the meshless numerical technique. Chinese Journal of Computational Mechanics, Vol. 27, Issue 2, 2010, p. 207212.

Qiao H. W., Lv Z. Z. Dynamic reliability sensitivity analysis for stochastic structures under random excitations. Journal of Vibration Engineering, Vol. 21, Issue 4, 2008, p. 404408.

Wu D., Wu Y., Yang Q. S., et al. Parameter sensitivity analysis of windinduced response of longspan roofs. Engineering Mechanics, Vol. 32, Issue 2, 2015, p. 171177.

Sun M., Tong L. P., Zhu Y. Z. Parameter sensitivity analysis of the earthquake input energy for rural building structures supported on sliding baseisolation systems. Engineering Mechanics, Vol. 33, Issue 5, 2016, p. 200210.

Dong S. L., Zhu X. L., Tu Y., et al. Configuration, simplified calculation method and parameter sensitivity analysis on prestress distribution of honeycombtype cable dome with two struts. Journal of Building Structures, Vol. 40, Issue 2, 2019, p. 128135.

Khaneghahi M. H., Alembagheri M., Soltani N. Reliability and variancebased sensitivity analysis of arch dams during construction and reservoir impoundment. Frontiers of Structural and Civil Engineering, Vol. 13, Issue 3, 2019, p. 526541.

Fu Z. F., Hua H. X. Theory and application of Modal Analysis. Shanghai Jiao Tong University Press, Shanghai, 2000.

Chopra Anil K. Dynamics of Structures. Fourth Edition, Prentice Hall Inc., Upper Saddle River, 2011.

Zhu W. Q. Random Vibration. Science Press, Beijing, 1992.

Isulapalli S. S., Roy A., Georgopoulos P. G. Stochastic response surface methods for uncertainty propagation: application to environmental and biological systems. Risk Analysis, Vol. 18, Issue 3, 1998, p. 351363.

Hu R., Li D. Q., Zhou C. B., et al. Structural reliability analysis using stochastic response surface method. Journal of Mechanical Engineering, Vol. 27, Issue 9, 2010, p. 192200.

Isukapalli S. S. An uncertainty analysis of transport transformation models. The State University of New Jersey, New Brunswick, 1999.

Zhang S. R., Tan P., Du Y. F., et al. Parameters optimization and energy analysis of interstory system based on Nondominated Sortig Genetic Algorithm with Elitism. Journal of Hefei University of Technology, Vol. 38, Issue 7, 2015, p. 949954.