1 Introduction
Hybrid simulation is an advanced structural testing method in which the critical component of the structure expected to experience instability or significant nonlinearity is tested physically in the laboratory while the rest of the structure is numerically analyzed using a parallel finite element model [12]. In earthquake engineering, there are generally two approaches to perform hybrid simulation: 1) pseudodynamic hybrid simulation (PsDHS), where the simulation is conducted at a slower pace compared to the time of ground motion acceleration, and 2) realtime hybrid simulation (RTHS) where the simulation is performed in realtime, offering a favorable option for structures sensitive to strain rate effects.
Hybrid simulation offers a viable and economical response evaluation method over the conventional quasistatic and shaking table testing methods. However, the results of hybrid simulation may become biased when only one or a limited number of critical elements (e.g., seismic fuses) in the structure is physically tested due to limited experimental testing resources, and the remaining critical elements are numerically modelled along with the rest of the structure using the mathematical model assumptions that may not be suitable to predict the response of such critical elements [5, 7]. For instance, Imanpour et al. [10] physically tested one of two critical columns of a steel multitiered concentrically braced frame, which may affect the overall response of the frame or even instability of the column tested in the laboratory. There are several other challenges associated with hybrid simulation, such as errors arising from test controller and hydraulic actuators, signal delays, and noise generated in instrumentation.
Machinelearning algorithms such as Prandtl neural network
[15], and leastsquares supportvector machine
[6]offer an efficient solution to some of the challenges described above. Such algorithms have been used in the past to develop metamodels for experimental substructures. Additionally, metamodels were created for finite element models in RTHS using the linear regression algorithm, and the recurrent neural network
[2].This paper proposes a new substructuring technique that employs a novel databased model for the seismic response evaluation of steel braced frames, and presents preliminary results used for verification purposes. The proposed databased model can learn the nonlinear hysteretic behaviour of structural elements (i.e., steel braces) using the sparse identification algorithm. This model is first trained using a pushover analysis performed on a numerical model of a steel buckling restrained brace (BRB), replacing actual experimental test data. A hybrid model of a prototype bucklingrestrained braced frame (BRBF) is then developed with a trained model simulating the hysteretic response of the BRB. The adequacy of datadriven hybrid simulation technique is finally evaluated by comparing the results of a dynamic analysis obtained using the proposed technique with those predicted by a pure numerical model of the same BRBF.
2 Datadriven Hybrid Simulation (DDHS) Technique
A new datadriven hybrid simulation technique is proposed here by incorporating two substructures, including a datadriven substructure that typically represents the structure’s critical component and a finite element model that simulates the remaining components of the structure. The architecture of the proposed technique is shown in Fig. 1. The datadriven substructure is powered by a machine learning algorithm trained using experimental test data capable of predicting the nonlinear dynamic response of critical components (e.g., seismic fuses in steel seismic forceresisting frames). The numerical substructure acting parallel to the datadriven substructure consists of validated elements that simulate the structural elements that essentially remain elastic (e.g., capacityprotected members in steel seismic forceresisting frames) under earthquake excitations.
2.1 Sparse Identification of Nonlinear Dynamics (SINDy) Algorithm
A new model discovery algorithm referred to as Sparse Identification of Nonlinear Dynamics (SINDy) was recently proposed by Brunton et al. [4] to extract the governing equation of dynamical systems from data using the sparse identification algorithm. SINDy will allow us to define the mathematical relationship, , between the component’s force and the component’s displacement using the experimental data. For this purpose, various potential functions are stacked in a library matrix, :
(1) 
in which,
is the displacement vector obtained from the experimental data used for training,
are the candidate functions (i.e. basis functions), and are the number of measured data points and candidate functions selected to form , respectively. The basis functions should be selected based on the nature of the phenomenon [4]. The higher number of basis functions involved in the library matrix can help improve discovering the relationship by the algorithm, but with more computational effort. By defining the library matrix, , approximated restoring force vector of the component can be represented as:(2) 
in which, is the coefficient vector that determines which combination of the basis functions should be used to best approximate the nonlinearity involved in the problem. Sparse regression algorithms such as regularized regression (LASSO) can be employed to find the active entries of the
vector. LASSO regression is chosen here due to its ability in shrinking irrelevant terms to zero, which helps promote sparsity and avoid overfitting
[16]. LASSO regression is formulated as the following minimization problem:(3) 
in which, is the restoring force vector obtained from training experimental data, and represent norm and norm, respectively. denotes the regularization parameter that controls the sparsity of the solution. By solving the LASSO minimization problem, is obtained, which is then used to discover the the function, , as follows:
(4) 
where, is the symbolic function of (displacement), in contrast to that is the data library matrix of displacement defined in Eq. 1.
2.2 PrandtlIshlinskii (PI) Model
Given that the efficacy of SINDy highly depends on the basis functions selected in the library matrix, it is crucial to appropriately select such functions to well reproduce the hysteretic response of nonlinear systems. One method to complement SINDy with a hysteretic memory is achieved by utilizing hysteretic operators such as the stop operators in PrandtlIshlinskii (PI) model [3].
PI model is a phenomenological operatorbased model defined as the weighted superposition of multiple stop (elasticperfectly plastic) operators. The stop operator, which was first introduced in the field of continuum mechanics to describe the elastoplastic behaviour of materials [3], is recursively defined as follows:
(5) 
(6) 
(7) 
in which, is the output of the stop operator for the given input signal , which is defined using a threshold . Within each stop operator, the time domain needs to be discretized into subintervals of such that the input signal, , becomes monotonic within each time step . By linearly combining multiple stop operators with different thresholds (r values), the discrete PI model can be defined:
(8) 
where, is the total output signal when the hysteretic system is subjected to an exciting input signal (). The PI model parameters, the weights and thresholds , should be optimized to fit the training data.
2.3 PISINDy Model for the Simulation of Hysteretic Response
The proposed model, which is referred to as PISINDy, combines the PI model and the SINDy algorithm to simulate the nonlinear hysteretic response of structural components. In the PISINDy model, the stop operators play the role of basis functions of the library matrix (Eq. 1), and the vector contains optimized weights () as given in Eq. 8 after the LASSO regression is carried out on the experimental data used for training. The PISINDy model is then defined as follows:
(9) 
In the proposed model, several potential thresholds need to be assumed for stop operators. Such thresholds should then be included in the library matrix so that SINDy can determine which combination of thresholds can best describe the hysteretic response of the experimental data used for training. The initial set of thresholds are defined based on the maximum deformation () as , in which, can be obtained using a sensitivity analysis. It is worth noting that the larger value of would result in a finer discretization but a computationally lesseffective optimization.
3 Numerical Verification
The efficiency of DDHS was evaluated using the nonlinear response history analysis (NLRHA) of a prototype twodimensional (2D) bucklingrestrained braced frame (BRBF). The properties of the frame are shown in Fig. 2. The pure numerical model of the BRBF is developed in the OpenSees program [13]. The beam and columns were modelled using elastic beamcolumn elements, and the GiuffréMenegottoPinto material model assigned to a truss element was used to construct the BRB [8]. For simplicity, the influence of frictional forces developed between the BRB core and restraint was neglected in the model, resulting in a symmetric hysteretic response under cyclic loading in tension and compression. The fundamental period of this structure is 0.492 seconds.
3.1 Training of the BRB using the PISINDy model
In lieu of actual experimental test data, fictitious test data is used in training, which was created using a cyclic nonlinear static (pushover) analysis performed on an identical isolated BRB under an increasing cyclic displacement protocol shown in Fig. 3a. PISINDy model was then trained with this input as described in Algorithm 1. 50 stop operators were selected to be used in the data library matrix, and was taken as 0.1 in LASSO regression. Normalized rootmeansquareerror (NRMSE) as given in Eq. 10 was used to compute the error between the approximated hysteretic response by the PISINDy model and that obtained from the pure numerical model.
(10) 
Fig. 3b shows the cyclic response of the isolated BRB obtained from the cyclic pushover analysis, compared with the BRB response simulated by PISINDy. As shown, the PISINDy successfully captured the nonlinear cyclic behaviour of the BRB with NRMSE = .
3.2 Verification of the DDHS Technique
The BRB in the BRBF of Fig. 2 was replaced with a PISINDypowered datadriven substructure to verify the capability of the proposed DDHS technique in predicting the restoring force of the BRB when the frame is subjected to dynamic loading. The training process of the PISINDy model is described in Section 3.1. PISINDy developed in MATLAB [11] was linked to the numerical substructure in OpenSees. The UTSIM framework [9, 14] was used to establish communication between two substructures. An NLRHA [1] was performed under the 2011 Miyagi MYG016 record applied in the horizontal direction of the frame in the OpenSees program. At the analysis time step, the displacement was sent to MATLAB to predict the restoring force. The predicted force vector in interfacial DOFs was then fed back to OpenSees to complete the numerical integration of the equation of motion at that time step.
Fig. 4 portrays the BRB hysteretic response and the history of storey drift ratio obtained from the NLRHA of the DDHS technique and the reference model. As shown, BRB responded in the inelastic range and showed a symmetrical response as expected. The frame underwent a maximum drift ratio of 1.43%, while the BRB experienced a maximum axial strain equal to 0.66%. The NRMSE of the storey drift ratio is 0.7702%, which implies the accuracy of the DDHS technique using the trained symmetrical hysteresis response of the BRB.
4 Conclusions
This paper proposes a new datadriven substructuring technique for the seismic evaluation of steel braced frame structures that is accomplished through incorporating a novel datadriven model (PISINDy) predicting the nonlinear cyclic response of the frame’s critical elements combined with a fibrebased numerical substructure simulating the wellunderstood components of the frame. PISINDy was developed using the sparse identification algorithm and inspired by the PrandtlIshlinskii model to predict the hysteretic response of BRB using available experimental test data. The preliminary results confirm that once PISINDy is trained, it can properly emulate the nonlinear cyclic response of a physical substructure. Moreover, for a symmetric hysteresis response, cyclic test data obtained from quasistatic testing would be enough to train PISINDy and reproduce the brace hysteretic response under random signals such as earthquake accelerations. The numerical verification results on the prototype single storey frame show the feasibility and suitability of the PISINDy model. Future studies should be devoted to further improve the proposed datadriven model for more complex hysteretic responses typically observed in steel braced frames. Furthermore, verification tests should be carried out in the future to verify the application of the proposed DDHS technique.
5 Acknowledgments
Funding for this project has been provided by the National Sciences and Engineering Research Council (NSERC). The authors gratefully acknowledge the CISC Centre for Steel Structures Education and Research at the University of Alberta for their continuous support, and the UTSIM framework team at the University of Toronto for their assistance in the development of the hybrid simulation platform. The first author acknowledges scholarships from the government of Alberta, the Norman and Tess Reid Family, and the Parya foundation.
References
 [1] (2016) Multidirectional structural component hybrid testing system for the assessment of the seismic response of steel ishaped columns. In Proceedings of the 11th Pacific Structural Steel Conference, Cited by: §3.2.
 [2] (2020) Using machine learning approach for computational substructure in realtime hybrid simulation. arXiv preprint arXiv:2004.02037. Cited by: §1.

[3]
(1996)
Hysteresis and phase transitions
. Vol. 121, Springer Science & Business Media. Cited by: §2.2, §2.2.  [4] (2016) Discovering governing equations from data by sparse identification of nonlinear dynamical systems. Proceedings of the national academy of sciences 113 (15), pp. 3932–3937. Cited by: §2.1.
 [5] (2010) Largescale cyclic and hybrid simulation testing and development of a controlledrocking steel building system with replaceable fuses. University of Illinois at UrbanaChampaign. Cited by: §1.
 [6] (2018) Hysteresis simulation using leastsquares support vector machine. Journal of Engineering Mechanics 144 (9), pp. 04018084. Cited by: §1.
 [7] (2020) Experimental evaluation of inelastic highermode effects on the seismic behavior of rc structural walls. Journal of Structural Engineering 146 (4), pp. 04020016. Cited by: §1.
 [8] (1983) Effects of bond deterioration on hysteretic behavior of reinforced concrete joints. Cited by: §3.
 [9] (2020) A generalized numerical/experimental distributed simulation framework. Journal of Earthquake Engineering 24 (4), pp. 682–703. Cited by: §3.2.
 [10] (2018) Development of a hybrid simulation computational model for steel braced frames. In Key Engineering Materials, Vol. 763, pp. 609–618. Cited by: §1.
 [11] (2019) MATLAB version 9.7.0.1216025 (R2019b) Update 1. The Mathworks, Inc., Natick, Massachusetts. Cited by: §3.2.
 [12] (2016) An overview of seismic hybrid testing of engineering structures. Engineering Structures 118, pp. 240–261. Cited by: §1.
 [13] (1999) Objectoriented finite element programming: frameworks for analysis, algorithms and parallel computing.. Cited by: §3.
 [14] Cited by: §3.2.
 [15] (2019) Neurohybrid simulation of nonlinear frames using prandtl neural networks. Proceedings of the Institution of Civil EngineersStructures and Buildings, pp. 1–18. Cited by: §1.
 [16] (1996) Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological) 58 (1), pp. 267–288. Cited by: §2.1.
Comments
There are no comments yet.