Next Article in Journal
Study on the Evolution and Coupling Coordinated Development of Passenger and Freight Transport Network of New Western Land-Sea Corridor
Next Article in Special Issue
Experimental Study on the Seismic Performance of New Energy Dissipation Connectors in an Autoclaved Aerated Concrete Panel with Assembled Steel Frame
Previous Article in Journal
Intelligent Fault Diagnosis of Broken Wires for Steel Wire Ropes Based on Generative Adversarial Nets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Study of Meta-Modeling for Response Estimation of Stochastic Nonlinear MDOF Systems Using MIMO-NARX Models

1
Key Laboratory of Concrete and Prestressed Concrete Structures, Ministry of Education, Southeast University, Nanjing 210096, China
2
School of Civil Engineering, Shandong University, Jinan 250061, China
3
School of Engineering, San Francisco State University, San Francisco, CA 94132, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(22), 11553; https://doi.org/10.3390/app122211553
Submission received: 20 September 2022 / Revised: 8 November 2022 / Accepted: 9 November 2022 / Published: 14 November 2022

Abstract

:
Complex dynamic behavior of nonlinear structures makes it challenging for uncertainty analysis through Monte Carlo simulations (MCS). Surrogate modeling presents an efficient and accurate computational alternative for a large number of MCS. The previous study has demonstrated that the multi-input multi-output nonlinear autoregressive with exogenous input (MIMO-NARX) model provides good discrete-time representations of deterministic nonlinear multi-degree-of-freedom (MDOF) structural dynamic systems. Model order reduction (MOR) is executed to eliminate insignificant modes to reduce the computational burden due to too many degrees of freedom. In this study, the MIMO-NARX strategy is integrated with different meta-modeling techniques for uncertainty analysis. Different meta-models including Kriging, polynomial chaos expansion (PCE), and arbitrary polynomial chaos (APC) are used to surrogate the NARX coefficients for system uncertainties. A nine-DOF structure is used as an MDOF dynamic system to evaluate different meta-models for the MIMO-NARX. Good fitness of statistical responses is observed between the MCS results of the original system and all surrogated MIMO-NARX predictions. It is demonstrated that the APC-NARX model with the advantage of being data-driven is the most efficient and accurate tool for uncertainty quantification of nonlinear structural dynamics.

1. Introduction

Uncertainty analysis in structural engineering has attracted considerable research efforts in recent years [1,2,3,4]. It accounts for the influence of various uncertain factors including material property, external load, ground motion excitation, and experimental measurement, therefore is important for the design and optimization of structures and experiments to meet the performance target of robustness or reliability and to reduce the adverse impact of uncertainty propagation on system performance. Monte Carlo simulation (MCS), as the most direct and effective method for uncertainty analysis [5,6], however, requires significant computational efforts to obtain the statistical information of the uncertain outputs. Meta modeling also known as surrogate, has recently emerged as an efficient method for uncertainty analysis. The meta-models are trained by a number of samples of the original model to approximately predict the response quantities of interests in the entire sample space of uncertain inputs. The most commonly used meta-models include Kriging [7], polynomial chaos expansion (PCE) [8], arbitrary polynomial chaos expansion (APC) [9], etc.
Engineering structures generally involve nonlinear complex mechanical behavior related to the time history-dependent states of responses. Required computational time significantly increases for accurate dynamic time-history simulations. The complex behavior exhibited by nonlinear systems brings great challenges to uncertainty analysis due to the requirement of a large number of samples for MCS. Dynamic analysis of a deterministic nonlinear structure can be considered as solving differential equations in the discrete-time domain. Nonlinear autoregressive with exogenous input (NARX) models provide good discrete-time representations of system dynamics as a function of previous response values as well as current and/or previous input signals [10,11]. The NARX model thus has been successfully applied for describing nonlinear dynamic behavior for single-degree-of-freedom (SDOF) systems [12]. In the case of multi-degree-of-freedom (MDOF) systems, the single-input single-output (SISO) strategy establishes a separate NARX model for each output of the MDOF system. Multi-input multi-output (MIMO) formulation considers the inevitable response coupling between the multiple outputs of the NARX model [13,14]. Nevertheless, a large number of degrees of freedom imposes significant computational demands for building the NARX model through the MIMO strategy. The model order reduction (MOR) method is usually utilized to reduce the degrees of freedom for the MIMO-NARX model, which has been demonstrated to significantly improve computational efficiency with little effect on its accuracy [15,16,17].
The NARX model can capture the nonlinear dynamic behavior of the deterministic system. A stochastic NARX model is characterized by a set of specified NARX model terms and associated random coefficients [18]. Once the NARX model is built from a series of nonlinear dynamic responses, the NARX coefficients need to be calibrated based on the specific system conditions including structural property and external excitation. This means that the NARX coefficients vary for different nonlinear dynamic behavior due to system uncertainties such as ground motions and structural properties. In uncertainty analysis, the nonlinear dynamic behavior of an MDOF system from a large number of MCS is difficult to obtain due to the unbearable computational cost of dynamic analysis. It means that the NARX coefficients cannot be calibrated directly by the target dynamic behaviors. Meta models have been used to govern the uncertainties in the nonlinear system represented by SISO-NARX in previous research [18,19]. Different from the Li et al. [14] method of adopting the constant mean value of the MIMO-NARX coefficients, this study builds different meta-models to surrogate the system uncertainties into the MIMO-NARX model coefficients for nonlinear system response prediction. More specifically, Kriging, PCE, and APC are utilized as meta-models to surrogate the NARX coefficients based on uncertain system parameters. The MIMO-NARX model is built in MOR-reduced coordinates using a series of numerical simulation samples to replicate the nonlinear dynamic behavior of the MDOF system with uncertainties. A nine-DOF structure is selected as a numerical case to evaluate the MIMO-NARX meta-models for nonlinear dynamic history prediction and to compare the effectiveness of different meta-models.

2. Brief Review of MIMO-NARX Modeling

2.1. MIMO-NARX Model for MDOF System

In structural dynamics, a nonlinear MDOF dynamic system can be expressed as:
M x ¨ t + C x ˙ t + f r x t = f t
where M and C are the mass and damping matrix, respectively; t denotes time; x ¨ t , x ˙ t and x t are response vectors of acceleration, velocity, and displacement, respectively; f r is the vector of nonlinear hysteresis forces; f t is the external excitation force, which is often expressed in earthquake engineering as f t = M ι x ¨ g t with an all-ones vector ι and the seismic acceleration x ¨ g t .
The nonlinear dynamic behavior of the MDOF system in Equation (1) can be captured in a NARX model by considering several steps of the current and past responses of input and output. A MIMO-NARX technique is adopted in this study to capture coupled nonlinear dynamics since a SISO-NARX model is unable to represent the response coupling between multiple outputs for general nonlinearity [13]. The current output x t of nonlinear MDOF dynamic system depend on its past output values x t Δ t , , x t n o Δ t with maximum time lag n o and time interval Δ t , and current and past input values f t , f t Δ t , , f t n i Δ t with maximum time lag n i and can be expressed as:
x t = G z t + ϵ t
where z t = x T t Δ t , , x T t n o Δ t , f T t , , f T t n i Δ t is the regression vector of current and past values of input and output; G · represents the MIMO-NARX model to be identified based on the input and output; and ϵ t is vector-valued Gaussian error process [20]. This study adopts the linear-in-the-parameter form for the MIMO-NARX model which can be expressed as:
G z t = Θ T g z t
where g · = g 1 T · ; g 2 T · , , g n T · T contains the NARX terms for the n DOFs system; and Θ = diag Θ 1 , , Θ n contains the NARX coefficients of the n NARX models. For each degree of freedom, the output of a MIMO-NARX model can be formulated as:
x j t = Θ j T g j z t + ϵ j t
where j represents the jth degree of freedom; g j z t i is the most relevant NARX terms selected from all potential NARX terms, that is, predetermined basis functions. While only the jth degree of freedom is considered in the SISO-NARX strategy, the MIMO-NARX strategy includes the input and output histories from all degrees of freedom.
Identification of the MIMO-NARX model requires the selection of NARX terms and calibration of their coefficient. An efficient approach [18] has been proposed for identifying the NARX model for SISO systems by implementing the least-angle regression (LARs) algorithm for structure determination and the ordinary least square (OLS) method for coefficient calibration. The LARs algorithm [21] is employed to select candidate NARX model terms by calculating the correlation between each potential model term and the system output [13]. The corresponding NARX coefficients can then be calculated based on the OLS method:
Θ j , k = Z j , k T Z j , k 1 Z j , k T X j T
where Z j , k is the terms matrix of the kth candidate NARX model; and X j is the response time series of the jth degree of freedom. The most appropriate MIMO-NARX model is selected after the terms and coefficients of the candidate NARX models are determined. An appropriate error measure is conducted by decoupling the identification of the NARX models for each degree of freedom. For the jth MIMO-NARX model, only the outputs of the jth degree of freedom are estimated by recursively running the model, while the outputs of the other degrees of freedom are directly obtained from the simulation results. The error measure is defined as:
e ˜ SE , j , k = || X j X ˜ j , k || 2 || X j T ι E t X j || 2
where E t X j is the mean value of X j ; and X ˜ j , k represents the response time series of the jth degree of freedom estimated from the kth candidate NARX model. The optimal model with model terms g j · and coefficients Θ j is selected as the candidate, which reaches a sufficiently small total error by satisfying a predefined threshold value. It is worth noting that as few as possible numbers of NARX terms have been found to avoid deleterious effects on the accuracy of the model, for example, overfitting and spurious dynamics [18,22,23].

2.2. Model Order Reduction

Although MIMO-NARX is suitable for multiple outputs, it is computationally expensive to establish the NARX model for each degree of freedom. The fact that the first few modes dominate the dynamic response of a multi-degree-of-freedom structure can be utilized for eliminating the burden due to a large number of degrees of freedom. The MOR method [15,16,17] thus can be utilized to lower the degrees of freedom of the MDOF nonlinear dynamic system in this study. Using a small number of lower structural modes corresponding to the number of DOF, an accurate approximation of the response history can be acquired. The response vector in the reduced space can be obtained through the coordinate transformation matrix and the response vector in the original space:
q t = Φ n r T x t
where q t is the response in the reduced n r -dimensional space; x t is the response vector in the original n-dimensional space; Φ n r T is the n r × n coordinate transformation matrix with n r n ;. The transformation matrix can be obtained by singular value decomposition (SVD) [24,25]. For the n s samples of the simulation results for the MDOF dynamic system under uncertainties, a snapshot matrix can be formulated as follows:
X = x 1 t 1 , , x 1 t n t , , x n s t 1 , , x n s t n t
where X n × n t n s contains n t time step responses of the n-DOF dynamic system. The SVD can then be used to decompose the snapshot matrix X as:
X = U Λ V T
where U is a n × n orthonormal matrix that contains the left singular vectors of X ; V is the n t n s × n t n s orthonormal matrix that contains the right singular vectors; Λ is a n × n t n s pseudo-diagonal matrix with the singular values Λ j , j = λ j , j = 1 , , n . A computationally effective economy-size SVD is used in view of the extremely large size of V with n t n s n [14]. The left singular vectors U are the proper orthogonal decomposition (POD) modes.
Then the coordinate transformation matrix is defined as the first n r POD modes as:
Φ n r = U 1 , U 2 , , U n r
where U 1 , U 2 , , U n r is the first n r columns of U ; n r can be determined by the energy criterion which ensures the energy captured in the truncated representation of x reaches η of the total energy [17]:
j = 1 n r λ j η j = 1 n λ j
where η is typically assumed to be close to 1, for example, 0.99. A properly selected η can bring significant dimensional reduction to the system, that is, nrn. The transformation error of MOR can be measured by the comparison between the original response x t from Equation (7) and the reverse-transformation response x _ t , which can be calculated as:
x _ t = Φ n r q t

2.3. NARX Model Pruning

The optimal NARX terms might still contain spurious terms due to the LARs approach selecting candidate terms through correlation analysis, which is not always a reflection of the contribution of a term to the model [22]. An iterative error-based pruning procedure [14] is adopted in this study to progressively delete the most deleterious NARX term at each iteration until a user-defined error threshold is reached. Each iteration intends to remove one unique term from the current MIMO-NARX model, and the coefficients of each trial model are recalculated through the OLS method. To evaluate the performance of trial models, an error indicator is defined to describe the accuracy of the MIMO-NARX model based on the responses of the mth DOF in the original coordinate:
e ^ SE , M = || X m Φ n r m Q ^ || 2 || X m T ι E t X m || 2
where m is the DOF of interest; X m is the response time series of the mth DOF; Φ n r m is the mth row of the transformation matrix Φ n r ; Q ^ is the response matrix predicted by the MIMO-NARX model. The accuracy of the current MIMO-NARX model in each iteration can be evaluated by the mean value of error indicators overall training samples e ^ ¯ SE , m . While the error indicator for each trial model can be denoted as e ^ ¯ SE , m ' . The current MIMO-NARX model is replaced by the optimal trial model with min e ^ ¯ SE , m ' e ^ ¯ SE , m in each iteration. The next iteration of the pruning process then proceeds with the new MIMO-NARX model serving as the current model until the predefined error threshold is satisfied, that is, e ^ ¯ SE , m ' e ^ ¯ SE , m > E ^ .

3. Meta-Modeling Techniques for Surrogating NARX Model Coefficients

In this study, meta-models are used to surrogate the NARX model coefficients to account for system uncertainties, which can be formulated as the d-dimensional input random variables ξ = [ ξ 1 , ξ 2 , , ξ d ] T d . System outputs can be generalized as Y = y ( ξ 1 , y ( ξ 2 ) , , y ( ξ N ) ] T for N groups of training samples. Three different meta-models are evaluated and compared in this study including Kriging, PCE, and APC.

3.1. Kriging Meta-Model

Kriging uses a weighted linear combination of all observed output values to estimate the system output of a random function or random process as a Gaussian process. These weights are generally described by the specified correlation function and are based on the distances between the location to be predicted and the locations already observed. The Kriging meta-model can be mathematically described as [26]:
y ξ = β T f ξ + ε ξ ,   β T f ξ = i = 1 D β i f i ξ
where β T f ξ is the mean value of the Gaussian process; f i ξ is the ith polynomial basis function with corresponding coefficients β i ; D represents the total number of basis function terms which depends on the degree of the polynomials. The residual term ε ξ is assumed to have zero mean and the following covariance:
Cov ε ξ i ,   ε ξ j = σ Z 2 𝓡 θ ,   ξ i ,   ξ j
where ξ i and ξ j represent the two different sample inputs; σ Z 2 is the process variance; and 𝓡 θ ,   ξ i ,   ξ j represents the correlation function with parameter θ. Several auto-correlation functions have been used in the literature [26,27,28,29]. In this study, the Gaussian correlation function is adopted, which can be expressed as:
𝓡 ξ i , ξ j = k = 1 d exp [ θ k ( ξ i , k ξ j , k ) 2 ]
where the parameters θ k are hyper-parameters of auto-correlation functions.
The Kriging model predicts the response for untried input ξ u through the best linear unbiased predictor:
y ^ ξ u = β ^ T f ξ u + 𝕽 T 𝓡 1 β ^ T f ξ u
σ M ^ 2 ξ u = σ ^ Z 2 1 f T ξ u   𝕽 T ξ u 0 𝓛 𝓛 𝓡 1 f ξ u 𝕽 ξ u
where 𝕽 ξ u = 𝓡 ξ u ,   ξ 1 ,   𝓡 ξ u ,   ξ 2 ,   ,   𝓡 ξ u ,   ξ N T defines the correlation matrix between untried sample point ξ u and training samples; 𝓛 is the matrix of basis functions at training points of dimension N × D; and 𝓡 is the correlation matrix for training samples.
The model parameters β ^ and σ ^ Z 2 , and the hyper-parameters θ can be determined using the maximum-likelihood estimation (MLE) [27,28]:
β ^ = ( 𝓛 T 𝓡 1 𝓛 ) 1 𝓛 T 𝓡 1 Y
σ ^ Z 2 = 1 N Y 𝓛 β ^ T 𝓡 1 Y 𝓛 β ^
More details of the Kriging meta-model can be referred to in the work of Santner et al. [28]. The MATLAB-based toolbox DACE [30] is used in this study for the Kriging model.

3.2. PCE Meta-Model

PCE is a weighted linear combination of orthogonal polynomial basis functions of input random variables [31] and has been widely used in many fields due to its solid mathematical foundation and good performance. For the d-dimensional input random variables ξ = [ ξ 1 , ξ 2 , , ξ d ] T d , the single system output y can be expanded through polynomial chaos expansion as:
y ξ = i = 0 D 1 C i Φ i ξ + ε
where Φ i ξ is the ith orthogonal multivariate basis function composed of Winer–Askey polynomials; the terms number D = d + r ! / d ! r ! with the PCE order r; and ε is the truncation error. The coefficients c = c 0 ,   c 1 , , c D 1 T can be determined through least square regression as:
c = ( A T A ) 1 A T Y
A = Φ 0 ξ 1 Φ 1 ξ 1 Φ D 1 ξ 1 Φ 0 ξ 2 Φ 1 ξ 2 Φ D 1 ξ 2 Φ 0 ξ N Φ 1 ξ N Φ D 1 ξ N ,   Y = Y ξ 1 Y ξ 2 Y ξ N
where A is the matrix of orthonormal multivariate polynomials; Y is the real value vector of the sample output; and N is the number of samples, which is suggested to be no less than twice the number of expansion terms D. In this study, the MATLAB-based UQLab toolbox [32] is utilized to construct the PCE model for the NARX model coefficients obtained from the least square method. The orthogonal polynomials are chosen based on the distribution of input variables. When the system has multiple outputs under uncertain inputs, one PCE surrogate model needs to be constructed for each system output.

3.3. APC Meta-Model

APC uses a data-driven technique to build the optimal orthogonal polynomial basis function without assuming the distribution types of inputs, which extends PCE to apply to arbitrary distributions of random input variables [9]. The system output y can be expressed as:
y ξ = i = 0 D 1 c i P i ξ + ε
where D is the number of APC expansion terms; P i ξ is the optimal orthogonal multivariate polynomial obtained through data-driven technique; ci is the corresponding expansion coefficient obtained from the least square method; The ith term orthogonal multivariate polynomial P i ξ is multiplied by univariate orthogonal polynomials as:
P i ξ = j = 1 d P j l i , j ξ j
j = 1 d l i , j = 0 , 1 , , r
where P j l i , j ξ j represents the l i , j th order univariate orthogonal polynomial of the jth dimensional input variable; the sum of the orders l i , j of the d univariate orthogonal polynomials ranges from 0 to r, and r represents the highest order of the multivariate APC models. The univariate orthogonal polynomial P j l ξ j is obtained based on statistical moments as in the following expression:
P j l ξ j = k = 0 l p k l X k
μ 0 μ 1 μ l μ 1 μ 2 μ l + 1 μ l 1 μ l μ 2 l 1 0 0 1 p 0 l p 1 l p l 1 l p l l = 0 0 0 1  
where the interior coefficients p k l of the univariate orthogonal polynomial are calculated from μ 0 to μ 2 l 1 , that is, the statistical moment with order 0 to (2l − 1) of the sample data of the univariate input variable ξ j . APC meta-model is directly obtained from the sample data and has demonstrated faster error convergence with the increasing order and better accuracy than PCE of the same order [33].

4. Training of MIMO-NARX Meta-Models

The MIMO-NARX model is trained offline using numerical simulation data of the structural system in this study. The procedure includes generating sample data from the nonlinear dynamic history simulations, training MIMO-NARX terms using reduced coordinate from model order reduction, and surrogating NARX coefficients through meta-models using the uncertain parameters of the nonlinear dynamic system.

4.1. Nine-DOF Shear Structure

A nine-DOF shear structure is used as the MDOF nonlinear dynamic system in this study, which is simplified from a nine-story benchmark frame [34] as shown in Figure 1a. The first three natural frequencies and mode shapes of this nine-DOF structure are shown in Figure 1b.
The equation of motion of the nine-DOF under seismic excitation can be expressed by Equation (1) and the nonlinear restoring vector f r is emulated using the Bouc–Wen model [35,36], which can be formulated as:
f r = k b w α b w x + 1 α b w z
z ˙ = A b w x ˙ γ b w x ˙ z n b w β b w x ˙ z n b w 1 z
where f r is the hysteresis force for each story; x is the inter-story displacement between two neighboring floors; and the Bouc–Wen model parameters ( α b w , Abw, γbw, βbw, and nbw) are listed in Table 1. Under a typically selected ground motion excitation, moderate nonlinearity can be observed in Figure 2 for the Bouc–Wen models of each story.

4.2. Uncertainties in Ground Motion and Structural Properties

The dynamic responses are affected by the system uncertainties in ground motion and structural properties. By varying ground motions and structural parameters, structural responses can be obtained and used as training samples. The stochastic ground motion model [37] is adopted in this study to generate ground motion uncertainties. Each ground motion is determined by six parameters I ¯ a , D 5 95 , t m i d , ω m i d , ω , ζ f and white noise. In this study, ω m i d and ζ f , which control the filter frequency and damping ratio, respectively, are selected as uncertain parameters of ground motion while the other parameters ( I ¯ a , D 5 95 , t m i d , ω ) are set as (0.0314 s·g, 11.23 s, 7.85 s, −0.04 Hz/s). It is worth noting that the white noise remains unchanged to generate seismic acceleration time histories. A typical acceleration time history is shown in Figure 3.
A total of n s = 200 parameter sets of ( ω m i d , ζ f ) are generated following the distribution in Figure 4a,b, which leads to 200 stochastic ground motions for MIMO-NARX training. Figure 4c and Figure 5 present the distribution of peak ground accelerations and the mean response spectra of these 200 acceleration time histories, respectively.
The concentrated masses of each story are considered structural uncertainties, which reflects the uncertainty of live load or temporary load changes in the building. The masses of floor 1, floor 2~8, and floor 9 are assumed to follow the Gaussian distribution as shown in Figure 6.
The two ground motion parameters along with the three structural parameters are considered uncertain input variables in this study and are listed in Table 2. The uncertain variables of the dynamic system can be generated randomly from the distribution in Table 2 for dynamic time history analysis. The sample data used for training the MIMO-NARX model are extracted from the current and past state values of structural responses.

4.3. MIMO-NARX Model Training

Firstly, the MOR strategy is applied to transform the nine-DOF to the reduced coordinates for the convenience of MIMO-NARX modeling. The snapshot matrix X is extracted from the structural responses of the two hundred dynamic history simulations. The transformation matrix Φ n r is obtained through SVD as the first n r POD modes according to Equations (9) and (10). The coordinate is transformed from the original order n = 9 to the reduced order n r = 4, which is determined by the energy criteria in Equation (11). The fitness indicator of the nine-DOF response in original coordinates and reverse-transformation response is shown in Figure 7a, in which the red line indicates the median, the blue box indicates the 25th and 75th percentiles, and the red plus symbol indicates outliers. It can be observed from Figure 7a that the maximum error of MOR occurs in the first DOF because the value of absolute response of the first floor is the smallest compared with the other floors. For a randomly selected sample, the comparison in Figure 7b shows a good fitness of 92% between original and reverse-transformation velocity responses x ˙ t and x ˙ _ t for the first floor. This implies that the error introduced by MOR is negligible in this study.
To capture the nonlinear dynamic behavior of the MDOF system, the MIMO-NARX model is trained from the n s = 200 dynamic history simulation samples. Multiple inputs and outputs in the reduced coordinates are then used to build the MIMO-NARX model to avoid significant computational efforts associated with the original number of DOFs. The simulation samples are first examined by a predefined threshold value to screen responses of interest that exhibit suitable nonlinearity [18]. NARX candidates are identified by the LARs algorithm and the corresponding term coefficients are calculated by the OLS approach as shown in Figure 8. The optimal NARX model is selected based on the minimum mean errors of these candidates and terms pruning is then applied to improve precision. If the accuracy requirement is not satisfied for the optimal NARX model, some measures to improve NARX accuracy need to be taken such as redesigning the formulation of the NARX full model.
In this study, the MIMO-NARX full model is designed as shown in Figure 9. The ground motion accelerations are considered as the input, while responses of the n r = 4 MOR reduced coordinates are used as the outputs. The maximum time lag of n i = n o = 4 is used for the input and output in the NARX model to capture the memories. The current value for each output is replicated by one MIMO-NARX model considering the coupling of the past signals of all four outputs.
The velocity responses in the reduced coordinates are selected as NARX outputs because previous studies showed that it is more effective to describe nonlinear systems with respect to velocity. The full NARX model consists of potential NARX terms including seismic acceleration from 0 to l steps forward x ¨ g t l Δ t , velocity responses of the reduced coordinates from 1 to l steps forward q ˙ j t l Δ t , one step ahead absolute value of velocity responses of the reduced coordinates q ˙ j t Δ t , q ˙ j t Δ t q ˙ j t l Δ t , q ˙ j t Δ t x ¨ g t l Δ t and 1, with l = 1, 2, 3, 4, j = 1, 2, 3, 4, which leads to 5 + 4 × 4 + 4 + 4 × 4 × 4 + 4 × 5 + 1 = 110 terms. It should be noticed that the absolute value terms directly use velocity responses of the reduced coordinates rather than velocity responses in the original coordinates Φ n r m q ˙ (t Δ t ) used in previous studies [14]. This results in a smaller number of terms to avoid redundant terms, which brings the benefits of faster computation and almost similar accuracy. The final MIMO-NARX model can be expressed as follows:
q ˙ 1 t = y 1 , 1 x ¨ g t + y 1 , 2 x ¨ g t Δ t + y 1 , 3 x ¨ g t 2 Δ t + y 1 , 4 x ¨ g t 3 Δ t + y 1 , 5 x ¨ g t 4 Δ t + + y 1 , 55 q ˙ 2 t Δ t x ¨ g t 4 Δ t + y 1 , 56 q ˙ 3 t Δ t x ¨ g t 4 Δ t + y 1 , 57 q ˙ 4 t Δ t x ¨ g t 4 Δ t
q ˙ 2 t = y 2 , 1 x ¨ g t + y 2 , 2 x ¨ g t Δ t + y 2 , 3 x ¨ g t 2 Δ t + y 2 , 4 x ¨ g t 3 Δ t + y 2 , 5 x ¨ g t 4 Δ t + + y 2 , 60 q ˙ 2 t Δ t x ¨ g t 4 Δ t + y 2 , 61 q ˙ 3 t Δ t x ¨ g t 4 Δ t + y 2 , 62 q ˙ 4 t Δ t x ¨ g t 4 Δ t
q ˙ 3 t = y 3 , 1 x ¨ g t + y 3 , 2 x ¨ g t 4 Δ t + y 3 , 3 q ˙ 2 t 3 Δ t + y 3 , 4 q ˙ 2 t 4 Δ t + y 3 , 5 q ˙ 3 t Δ t + + y 3 , 27 q ˙ 2 t Δ t x ¨ g t + y 3 , 28 q ˙ 3 t Δ t x ¨ g t + y 3 , 29 q ˙ 4 t Δ t x ¨ g t
q ˙ 4 t = y 4 , 1 x ¨ g t + y 4 , 2 x ¨ g t Δ t + y 4 , 3 x ¨ g t 2 Δ t + y 4 , 4 x ¨ g t 3 Δ t + y 4 , 5 x ¨ g t 4 Δ t + + y 4 , 47 q ˙ 2 t Δ t x ¨ g t 4 Δ t + y 4 , 48 q ˙ 3 t Δ t x ¨ g t 4 Δ t + y 4 , 49 q ˙ 4 t Δ t x ¨ g t 4 Δ t
where y i , j represents the NARX coefficient for the jth NARX term of the ith output. It can be observed that there are 57, 62, 29, and 49 terms for the four outputs q ˙ 1 t , q ˙ 2 t , q ˙ 3 t , and q ˙ 4 t , respectively. Higher order of DOF does not necessarily lead to more terms for the MIMO-NARX model.

4.4. Surrogating NARX Coefficients by Meta-Models

The MIMO-NARX model trained by the simulation samples is aimed to determine the optimal NARX candidate terms from Figure 9. The corresponding NARX coefficients need to be calibrated based on the system responses in the specific dynamic simulation case. For the n s = 200 simulation samples, the five parameters of ground motion and structural masses in Table 2 are considered uncertain factors affecting the dynamic responses. Therefore, the coefficients of the MIMO-NARX model are taken as output variables varying with the five input variables of system uncertainties. The meta-models are then built to surrogate the NARX coefficients using the five uncertain parameters under the condition that the responses of a large number of MCS are unknown as shown in Figure 10. This study explores three meta-models including Kriging, PCE, and APC. The order of 3 is used for PCE and APC as the recommended number of samples is twice the number of polynomial terms.

5. Validation of Surrogated MIMO-NARX Model

To explore the uncertainty propagation in a nonlinear dynamic system with MDOF, a total of 300 validation Monte Carlo simulations with different ground motions and structural properties are performed in this section by traditional dynamic history analysis using the CR algorithm [38]. The MIMO-NARX models obtained from the previous section are adopted to capture the nonlinear dynamic behavior while the corresponding coefficients are calibrated by the meta-models for each set of uncertainty variable values. Finally, structural responses are compared between the MIMO-NARX model prediction and MCS to evaluate its effectiveness.
An error evaluation criterion is introduced to measure the fitness between the predicted and original value and expressed as:
Fitness = 1 x i x s i 2 x i x ¯ 2 × 100 %
where x i represent the original response from MCS with mean x ¯ ; and x s i represent the predicted response by the surrogated MIMO-NARX model.

5.1. Accuracy Assessment of MIMO-NARX Outputs in Reduced Coordinates

The MIMO-NARX model replicates four outputs of the order-reduced model in reduced coordinates as shown in Figure 11 for a typically selected ground motion with the four reduced coordinates presented from low to high. It can be observed in Figure 11 that the outputs of the MIMO-NARX model with the coefficients calibrated by Kriging, PCE, and APC meta-models are basically consistent with the MCS responses transformed to reduced coordinates. The matching accuracy of the four reduced coordinates, however, decreases with the increase of frequency due to the difficulty in capturing high-frequency nonlinearity as shown in Figure 12. For this specific selected ground motion, the APC meta-model has better accuracy for the fourth output prediction than that of Kriging and APC meta-models, while the three meta-models have nearly similar accuracy for the first three outputs.

5.2. Accuracy Assessment of Absolute Structural Responses

Since the absolute velocity responses are used in NARX training to better capture nonlinearity, Figure 13a first presents the comparison of absolute velocity responses for the nine-DOF structure, which are obtained from the four outputs of the reduced coordinates from the MIMO-NARX model. Good fitness can be observed between Kriging/PCE/APC-NARX model prediction and MCS results. The absolute displacement responses are then acquired through the integration of corresponding velocity responses and presented in Figure 13b. Slight errors can be observed at the end of time history analysis which can be attributed to integral error and error accumulation. Figure 14 compares the fitness index of absolute responses between three meta-models. Among the three meta-models, almost exactly the same fitness results can be observed for absolute velocity responses. Figure 14a shows that the first floor still has the lowest accuracy compared with the other floors, which is consistent with MOR errors in Figure 7. It can also be observed from Figure 14b that the Kriging meta-model has worse accuracy in absolute displacement prediction than the PCE and APC meta-models due to the drift errors at the end of the time history response. The drift errors are obvious in the lower and upper floors but less significant in the middle floors for the Kriging-NARX results of the selected ground motion.

5.3. Accuracy Assessment of Inter-Story Structural Responses

Inter-story velocity and displacement responses are further calculated from the absolute velocity and displacement responses of the nine-DOF structure. Figure 15a,b presents the comparison of inter-story velocity and displacement responses for the nine-DOF structure, respectively. As can be seen from the figures, the inter-story velocity curves still have good fitting accuracy, while the inter-story displacement curves show more obvious drift error from integral accumulation error, especially in the upper floors. Figure 16 compares the fitness of inter-story responses between three meta-models. It can be observed that the accuracy of both inter-story velocity and displacement responses slightly decreases compared with those of absolute responses. The Kriging-NARX model shows a significant reduction in the accuracy of upper floors compared with PCE/APC-NARX models due to the error accumulation at the end of time history analysis. Good fitness of inter-story velocity can still be observed between Kriging/PCE/APC-NARX model prediction and MCS results.

5.4. Comparison of Peak Story Responses

The safety and serviceability requirements for structural design mainly depend on the maximum peak responses. Figure 17a compares the peak inter-story velocity responses between Kriging/PCE/APC-NARX model prediction and MCS. It can be found that peak values from Kriging/PCE/APC-NARX models align well with actual values from MCS validation cases. The peak values of inter-story displacement responses are further compared in Figure 17b, where good agreement again can be observed. It demonstrates that the maximum peak displacements are generally close to each other and are not significantly affected by the accumulation errors from integrating the time history of velocity responses. Correlation coefficients between the actual and predicted maximum peak values of velocity and displacement responses are further compared in Figure 18. It should be noticed that the correlation coefficients are relatively low due to some extreme data points commonly referred to as outliers. Figure 18 implies that the correlation coefficients for peak displacement are lower than those of peak velocity due to the existence of more outliers. For peak responses, it can be further observed that PCE and APC meta-models have more stable performance than the Kriging meta-models.

5.5. Statistics for All Validation Cases

The statistical information of the total 300 validation samples is further summarized for comparison of the Kriging/PCE/APC-NARX models and traditional MCS in Figure 19. Figure 19a shows the box plot for the mean fitness of the four outputs in reduced coordinates, in which the red line indicates the median, the blue box indicates the 25th and 75th percentiles, and the red plus symbol indicates outliers. It demonstrates that all three meta-model-surrogated MIMO-NARX models have nearly similar accuracy for most MCS validation cases except for some outliers. Figure 19b shows the histogram with a distribution fitting curve of the mean fitness of the four outputs in reduced coordinates, which illustrates that APC-NARX has more counts for higher fitness than the corresponding Kriging- and PCE-NARX models. Figure 19c–f shows the statistical information for the mean fitness of the nine-DOF absolute velocity and displacement responses. It can be observed from Figure 19c,e that the Kriging/PCE/APC-NARX strategy is quite accurate in velocity prediction while the accuracy of displacement obtained by integration decreases to a certain extent. Figure 19d,f indicates that PCE- and APC-NARX are superior to Kriging-NARX in terms of the counts of good fitness. Figure 19g–j shows the statistical information for the mean fitness of the nine inter-story responses of velocity and displacement. The accuracy is observed to decrease when compared with those in Figure 19c–f. Lower accuracy brings more outliers, which interferes with the comparison of the three meta-model surrogated MIMO-NARX models.
Figure 20 shows the variation in the accuracy of Kriging/PCE/APC-NARX model prediction with respect to the mean ductility ratio of the nine inter-story Bouc–Wen hysteresis. The fitness and ductility in Figure 20 represent the mean value of the fitness and ductility ratio of the four outputs or nine responses. As shown in Figure 20a–e, the accuracy of the Kriging/PCE/APC-NARX strategy significantly decrease with the increase of inter-story ductility. Figure 20c,e illustrates this phenomenon again that the predicted accuracy of displacement is generally inferior to that of velocity due to the accumulation error. It can be demonstrated that the Kriging/PCE/APC-NARX strategy is basically accurate for most MCS validation cases within a certain nonlinearity range.

6. Conclusions

This research presents a computational study of different surrogate strategies to account for uncertainty propagation in nonlinear MDOF structural systems. A nine-DOF structure is selected for the numerical study for generating the training and validation samples. MIMO-NARX modeling is established to capture the nonlinear dynamic behavior in MOR-reduced coordinates using the results from dynamic time history simulations. The NARX coefficients are surrogated using meta-models including Kriging, PCE, and APC with uncertain variables of ground motion and structural mass.
The MCS results show good fitness of responses between MIMO-NARX models and the traditional dynamic simulation using the integration algorithm. It is worth noting that the absolute velocity responses are used in NARX training for improving the nonlinear capture effectiveness. The accuracy of the four reduced coordinates decreases with the increase of frequency due to difficulty in capturing high-frequency nonlinearity. Good performances are observed for all three meta-models in absolute velocity response prediction when compared with MCS. Generally, a slight decrease in accuracy is observed in displacement response prediction obtained by integrating the velocity responses due to the drift errors caused by error accumulation at the end of time history analysis. The accuracy of both inter-story velocity and displacement responses show different degrees of decrease compared with absolute velocity and displacement responses. The integral and accumulation drift errors do not necessarily affect the meta-model predicted peak displacement responses, which are more concerned with structural safety and serviceability. From a statistical point of view, all three meta-models surrogated MIMO-NARX models have nearly similar accuracy for most of the three hundred MCS validation cases except for some outliers. PCE- and APC- are shown to be superior to Kriging-NARX in terms of the counts of good fitness according to some statistical information. Considering the advantages of APC over PCE in the scope of application, APC-NARX is best suitable to replace the original model in a large number of MCS. It is illustrated that the MIMO-NARX modeling with surrogated coefficients by meta-models provides an effective and efficient tool for uncertainty quantification of nonlinear structures in earthquake engineering.
This study takes a nine-DOF structure as a proof of concept to verify the effectiveness of the proposed method. The nonlinear performance of the Bouc–Wen model in the MDOF system changes with different stochastic seismic excitations and structural properties. The choice of the basic functions in MIMO-NARX, which depends on different nonlinear types of structural models, remains to be studied in the future.
It should be noticed that the accuracy of the MIMO-NARX model is affected and limited by the intensity of nonlinearity. Although the accuracy of the three meta-model-surrogated MIMO-NARX models is significantly decreased with the intensity of structural ductility increase, it still can be demonstrated that the Kriging/PCE/APC-NARX strategy is basically accurate for most of the MCS validation cases within a certain nonlinearity range. Future research may focus on finding more suitable NARX basis functions to eliminate the adverse effect of high nonlinearity on accuracy.

Author Contributions

Conceptualization, C.C.; methodology, C.C., X.G. and M.C.; software, X.G. and M.C.; validation, M.C.; formal analysis, M.C.; investigation, X.G. and M.C.; resources, X.G. and M.C.; data curation, M.C.; writing—original draft preparation, M.C.; writing—review and editing, C.C.; visualization, M.C.; supervision, T.G.; project administration, T.G.; funding acquisition, T.G. and W.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of the People’s Republic of China (Grant No. 2018YFE0206100), the National Natural Science Foundation of China (Grant No. 52178114 and No. 51808111), Young scientific and technological talents promotion project of Jiangsu Association for science and technology (No. 2021-79).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kundu, A.; Adhikari, S. Transient Response of Structural Dynamic Systems with Parametric Uncertainty. J. Eng. Mech. 2014, 140, 315–331. [Google Scholar] [CrossRef] [Green Version]
  2. Silva, F.M.A.; Goncalves, P.B. The influence of uncertainties and random noise on the dynamic integrity analysis of a system liable to unstable buckling. Nonlinear Dyn. 2015, 81, 707–724. [Google Scholar] [CrossRef]
  3. Chakraborty, S.; Chowdhury, R. Polynomial Correlated Function Expansion for Nonlinear Stochastic Dynamic Analysis. J. Eng. Mech. 2015, 141, 04014132. [Google Scholar] [CrossRef]
  4. Mai, C.V.; Sudret, B. Surrogate Models for Oscillatory Systems Using Sparse Polynomial Chaos Expansions and Stochastic Time Warping. SIAM-ASA J. Uncertain. Quantif. 2017, 5, 540–571. [Google Scholar] [CrossRef] [Green Version]
  5. Echard, B.; Gayton, N.; Lemaire, M. AK-MCS: An active learning reliability method combining Kriging and Monte Carlo Simulation. Struct. Saf. 2011, 33, 145–154. [Google Scholar] [CrossRef]
  6. Zhang, Z.; Ma, X.; Hua, H.; Liang, X. Nonlinear stochastic dynamics of a rub-impact rotor system with probabilistic uncertainties. Nonlinear Dyn. 2020, 102, 2229–2246. [Google Scholar] [CrossRef]
  7. Gidaris, I.; Taflanidis, A.A.; Mavroeidis, G.P. Kriging metamodeling in seismic risk assessment based on stochastic ground motion models. Earthq. Eng. Struct. Dyn. 2015, 44, 2377–2399. [Google Scholar] [CrossRef]
  8. Blatman, G.; Sudret, B. An adaptive algorithm to build up sparse polynomial chaos expansions for stochastic finite element analysis. Probabilistic Eng. Eng. Mech. 2010, 25, 183–197. [Google Scholar] [CrossRef]
  9. Oladyshkin, S.; Nowak, W. Data-driven uncertainty quantification using the arbitrary polynomial chaos expansion. Reliab. Eng. Syst. Saf. 2012, 106, 179–190. [Google Scholar] [CrossRef]
  10. Leontaritis, I.J.; Billings, S.A. Input-output parametric models for non-linear systems Part I: Deterministic non-linear systems. Int. J. Control 1985, 41, 303–328. [Google Scholar] [CrossRef]
  11. Leontaritis, I.J.; Billings, S.A. Input-output parametric models for non-linear systems Part II: Stochastic non-linear systems. Int. J. Control 1985, 41, 329–344. [Google Scholar] [CrossRef]
  12. Gao, X.S.; Hou, H.T.; Huang, L.; Yu, G.Q.; Chen, C. Evaluation of Kriging-NARX modeling for uncertainty quantification of nonlinear SDOF systems with degradation. Int. J. Struct. Stab. Dyn. 2021, 21, 2150060. [Google Scholar] [CrossRef]
  13. Billings, S.A.; Chen, S.; Korenberg, M.J. Identification of MIMO non-linear systems using a forward-regression orthogonal estimator. Int. J. Control 1989, 49, 2157–2189. [Google Scholar] [CrossRef]
  14. Li, B.W.; Chuang, W.C.; Spence, S.M.J. Response Estimation of Multi-Degree-of-Freedom Nonlinear Stochastic Structural Systems through Metamodeling. J. Eng. Mech. 2021, 147, 04021082. [Google Scholar] [CrossRef]
  15. Grigoriu, M. Reduced order models for random functions. Application to stochastic problems. Appl. Math. Model. 2009, 33, 161–175. [Google Scholar] [CrossRef]
  16. Jensen, H.A.; Munoz, A.; Papadimitriou, C.; Millas, E. Model reduction techniques for reliability-based design problems of complex structural systems. Reliab. Eng. Syst. Saf. 2016, 149, 204–217. [Google Scholar] [CrossRef]
  17. Bamer, F.; Amiri, A.K.; Bucher, C. A new model order reduction strategy adapted to nonlinear problems in earthquake engineering. Earthq. Eng. Struct. Dyn. 2017, 46, 537–559. [Google Scholar] [CrossRef]
  18. Mai, C.V.; Spiridonakos, M.D.; Chatzi, E.N.; Sudret, B. Surrogate modeling for stochastic dynamical systems by combining nonlinear autoregressive with exogenous input models and polynomial chaos expansions. Int. J. Uncertain. Quantif. 2016, 6, 313–339. [Google Scholar] [CrossRef]
  19. Bhattacharyya, B.; Jacquelin, E.; Brizard, D. A Kriging-NARX model for uncertainty quantification of nonlinear stochastic dynamical systems in time domain. J. Eng. Mech. 2020, 146, 1–21. [Google Scholar] [CrossRef]
  20. Leontaritis, I.J.; Billings, S.A. Experimental design and identifiability for non-linear systems. Int. J. Syst. Sci. 1987, 18, 189–202. [Google Scholar] [CrossRef]
  21. Efron, B.; Hastie, T.; Johnstone, I.; Tibshirani, R. Least angle regression. Ann. Stat. 2004, 32, 407–451. [Google Scholar] [CrossRef] [Green Version]
  22. Piroddi, L.; Spinelli, W. An identification algorithm for polynomial NARX models based on simulation error minimization. Int. J. Control 2003, 76, 1767–1781. [Google Scholar] [CrossRef]
  23. Billings, S.A. Nonlinear System Identification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains; Wiley: Chichester, UK, 2013. [Google Scholar]
  24. Kerschen, G.; Golinval, J.C. Physical interpretation of the proper orthogonal modes using the singular value decomposition. J. Sound Vib. 2002, 249, 849–865. [Google Scholar] [CrossRef] [Green Version]
  25. Phalippou, P.; Bouabdallah, S.; Breitkopf, P.; Villon, P.; Zarroug, M. Sparse POD modal subsets for reduced-order nonlinear explicit dynamics. Int. J. Numer. Methods Eng. 2019, 121, 763–777. [Google Scholar] [CrossRef]
  26. Kaymaz, I. Application of kriging method to structural reliability problems. Struct. Saf. 2005, 27, 133–151. [Google Scholar] [CrossRef]
  27. Sacks, J.; Welch, W.J.; Mitchell, T.J.; Wynn, H.P. Design and analysis of computer experiments. Stat. Sci. 1989, 4, 409–423. [Google Scholar] [CrossRef]
  28. Santner, T.J.; Williams, B.J.; Notz, W.I.; Williams, B.J. The Design and Analysis of Computer Experiments; Springer: New York, NY, USA, 2003. [Google Scholar]
  29. Schobi, R.; Sudret, B.; Wiart, J. Polynomial-chaos-based Kriging. Int. J. Uncertain. Quantif. 2015, 5, 171–193. [Google Scholar] [CrossRef] [Green Version]
  30. Nielsen, H.B.; Lophaven, S.N.; Søndergaard, J. DACE, a MATLAB Kriging Toolbox, Informatics and Mathematical Modelling; Technical University of Denmark: Lyngby, Denmark, 2002. [Google Scholar]
  31. Xiu, D.; Karniadakis, G.E. The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations. SIAM J. Sci. Comput. 2002, 24, 619–644. [Google Scholar] [CrossRef]
  32. Marelli, S.; Sudret, B. UQLab: A framework for uncertainty quantification in Matlab. In Proceedings of the 2nd International Conference on Vulnerability, Risk Analysis and Management (ICVRAM2014), Liverpool, UK, 13–16 July 2014. [Google Scholar]
  33. Chen, M.H.; Guo, T.; Chen, C.; Xu, W.J. Data-driven Arbitrary Polynomial Chaos Expansion on Uncertainty Quantification for Real-time Hybrid Simulation Under Stochastic Ground Motions. Exp. Tech. 2020, 44, 751–762. [Google Scholar] [CrossRef]
  34. Ohtori, Y.; Christenson, R.E.; Spencer, B.F.; Dyke, S.J. Benchmark control problems for seismically excited nonlinear buildings. J. Eng. Mech. 2004, 130, 366–385. [Google Scholar] [CrossRef]
  35. Wen, Y.K. Method for random vibration of hysteretic systems. J. Eng. Mech. ASCE 1976, 102, 249–263. [Google Scholar] [CrossRef]
  36. Ikhouane, F.; Hurtado, J.E.; Rodellar, J. Variation of the hysteresis loop with the Bouc-Wen model parameters. Nonlinear Dyn. 2007, 48, 361–380. [Google Scholar] [CrossRef]
  37. Rezaeian, S.; Kiureghian, A.D. Simulation of synthetic ground motions for specified earthquake and site characteristics. Earthq. Eng. Struct. Dyn. 2010, 39, 1155–1180. [Google Scholar] [CrossRef]
  38. Chen, C.; Ricles, J.M. Development of Direct Integration Algorithms for Structural Dynamics Using Discrete Control Theory. J. Eng. Mech. 2008, 134, 676–683. [Google Scholar] [CrossRef]
Figure 1. Nine-DOF shear structure simplified from nine-story frame. (a) Target structure; (b) Natural frequencies and modes.
Figure 1. Nine-DOF shear structure simplified from nine-story frame. (a) Target structure; (b) Natural frequencies and modes.
Applsci 12 11553 g001
Figure 2. Hysteresis behavior of Bouc–Wen models (a typical sample).
Figure 2. Hysteresis behavior of Bouc–Wen models (a typical sample).
Applsci 12 11553 g002
Figure 3. Time history of stochastic ground motion (a typical sample).
Figure 3. Time history of stochastic ground motion (a typical sample).
Applsci 12 11553 g003
Figure 4. Uncertainties of ground motions (200 samples). (a) Filter frequency ω m i d ; (b) Filter damping ratio ζ f ; (c) Peak ground acceleration.
Figure 4. Uncertainties of ground motions (200 samples). (a) Filter frequency ω m i d ; (b) Filter damping ratio ζ f ; (c) Peak ground acceleration.
Applsci 12 11553 g004
Figure 5. Response spectra of ground motions with damping ratio 0.05 (200 samples).
Figure 5. Response spectra of ground motions with damping ratio 0.05 (200 samples).
Applsci 12 11553 g005
Figure 6. Variable mass for structural uncertainties (200 samples). (a) Floor 1; (b) Floor 2~8; (c) Floor 9.
Figure 6. Variable mass for structural uncertainties (200 samples). (a) Floor 1; (b) Floor 2~8; (c) Floor 9.
Applsci 12 11553 g006
Figure 7. MOR error between the original and reverse-transformation responses. (a) Fitness for nine DOFs (200 samples); (b) 1st DOF response (a random sample).
Figure 7. MOR error between the original and reverse-transformation responses. (a) Fitness for nine DOFs (200 samples); (b) 1st DOF response (a random sample).
Applsci 12 11553 g007
Figure 8. Procedure for MIMO-NARX training.
Figure 8. Procedure for MIMO-NARX training.
Applsci 12 11553 g008
Figure 9. Structure of MIMO-NARX full model.
Figure 9. Structure of MIMO-NARX full model.
Applsci 12 11553 g009
Figure 10. Meta-models to surrogate MIMO-NARX coefficients.
Figure 10. Meta-models to surrogate MIMO-NARX coefficients.
Applsci 12 11553 g010
Figure 11. Time history of the four outputs in reduced coordinates (a typical sample).
Figure 11. Time history of the four outputs in reduced coordinates (a typical sample).
Applsci 12 11553 g011
Figure 12. Fitness for the four outputs (a typical sample).
Figure 12. Fitness for the four outputs (a typical sample).
Applsci 12 11553 g012
Figure 13. Time history of the absolute responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Figure 13. Time history of the absolute responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Applsci 12 11553 g013aApplsci 12 11553 g013b
Figure 14. Fitness for the absolute responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Figure 14. Fitness for the absolute responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Applsci 12 11553 g014
Figure 15. Time history of the inter-story responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Figure 15. Time history of the inter-story responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Applsci 12 11553 g015aApplsci 12 11553 g015b
Figure 16. Fitness for the inter-story responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Figure 16. Fitness for the inter-story responses (a typical sample). (a) Velocity responses; (b) Displacement responses.
Applsci 12 11553 g016
Figure 17. Peak values of velocity responses (300 samples). (a) Peak velocity; (b) Peak displacement.
Figure 17. Peak values of velocity responses (300 samples). (a) Peak velocity; (b) Peak displacement.
Applsci 12 11553 g017aApplsci 12 11553 g017b
Figure 18. Correlation coefficients for peak responses (300 samples). (a) Peak velocity; (b) Peak displacement.
Figure 18. Correlation coefficients for peak responses (300 samples). (a) Peak velocity; (b) Peak displacement.
Applsci 12 11553 g018
Figure 19. Boxplot and histogram of the mean fitness of different time history curves (300 samples). (a) The mean fitness of the four outputs in reduced coordinates: Boxplot; (b) Histogram; (c) The mean fitness of the nine-DOF absolute velocity responses: Boxplot; (d) Histogram; (e) The mean fitness of the nine-DOF absolute displacement responses: Boxplot; (f) Histogram; (g) The mean fitness of the nine inter-story velocity responses: Boxplot; (h) Histogram; (i) The mean fitness of the nine inter-story displacement responses: Boxplot; (j) Histogram.
Figure 19. Boxplot and histogram of the mean fitness of different time history curves (300 samples). (a) The mean fitness of the four outputs in reduced coordinates: Boxplot; (b) Histogram; (c) The mean fitness of the nine-DOF absolute velocity responses: Boxplot; (d) Histogram; (e) The mean fitness of the nine-DOF absolute displacement responses: Boxplot; (f) Histogram; (g) The mean fitness of the nine inter-story velocity responses: Boxplot; (h) Histogram; (i) The mean fitness of the nine inter-story displacement responses: Boxplot; (j) Histogram.
Applsci 12 11553 g019
Figure 20. Relationship between the mean ductility ratio and the mean fitness (300 samples). (a) Outputs in reduced coordinates; (b) Absolute velocity responses; (c) Absolute displacement responses; (d) Inter-story velocity responses; (e) Inter-story displacement responses.
Figure 20. Relationship between the mean ductility ratio and the mean fitness (300 samples). (a) Outputs in reduced coordinates; (b) Absolute velocity responses; (c) Absolute displacement responses; (d) Inter-story velocity responses; (e) Inter-story displacement responses.
Applsci 12 11553 g020
Table 1. Bouc–Wen model parameters.
Table 1. Bouc–Wen model parameters.
Floor k b w   ( N / mm ) α b w A b w γ b w β b w n b w
12.9793 × 1050.0110.020.020.01
27.3787 × 1050.0110.10.10.01
36.9252 × 1050.0110.10.10.01
46.0177 × 1050.0110.10.10.01
55.2397 × 1050.0110.080.080.01
64.4108 × 1050.0110.060.060.01
73.6986 × 1050.0110.060.060.01
83.4706 × 1050.0110.080.080.01
93.2748 × 1050.0110.10.10.01
Table 2. Marginal distributions of uncertainties.
Table 2. Marginal distributions of uncertainties.
VariableDistribution TypeBoundsMeanStandard Deviation
ω m i d / 2 π Gamma[0, +∞]5.873.11
ζ f Beta[0.02, 1]0.2130.143
Mass 1Gaussian[392.77, 729.43]561.156.11
Mass 2Gaussian[384.58, 714.22]549.454.94
Mass 3Gaussian[416.08, 772.72]594.459.44
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, M.; Gao, X.; Chen, C.; Guo, T.; Xu, W. A Comparative Study of Meta-Modeling for Response Estimation of Stochastic Nonlinear MDOF Systems Using MIMO-NARX Models. Appl. Sci. 2022, 12, 11553. https://doi.org/10.3390/app122211553

AMA Style

Chen M, Gao X, Chen C, Guo T, Xu W. A Comparative Study of Meta-Modeling for Response Estimation of Stochastic Nonlinear MDOF Systems Using MIMO-NARX Models. Applied Sciences. 2022; 12(22):11553. https://doi.org/10.3390/app122211553

Chicago/Turabian Style

Chen, Menghui, Xiaoshu Gao, Cheng Chen, Tong Guo, and Weijie Xu. 2022. "A Comparative Study of Meta-Modeling for Response Estimation of Stochastic Nonlinear MDOF Systems Using MIMO-NARX Models" Applied Sciences 12, no. 22: 11553. https://doi.org/10.3390/app122211553

APA Style

Chen, M., Gao, X., Chen, C., Guo, T., & Xu, W. (2022). A Comparative Study of Meta-Modeling for Response Estimation of Stochastic Nonlinear MDOF Systems Using MIMO-NARX Models. Applied Sciences, 12(22), 11553. https://doi.org/10.3390/app122211553

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop