Novel coupling Rosenbrock‐based algorithms for real‐time dynamic substructure testing

Real‐time testing with dynamic substructuring is a novel experimental technique capable of assessing the behaviour of structures subjected to dynamic loadings including earthquakes. The technique involves recreating the dynamics of the entire structure by combining an experimental test piece consisting of part of the structure with a numerical model simulating the remainder of the structure. These substructures interact in real time to emulate the behaviour of the entire structure. Time integration is the most versatile method for analysing the general case of linear and non‐linear semi‐discretized equations of motion. In this paper we propose for substructure testing, L‐stable real‐time (LSRT) compatible integrators with two and three stages derived from the Rosenbrock methods. These algorithms are unconditionally stable for uncoupled problems and entail a moderate computational cost for real‐time performance. They can also effectively deal with stiff problems, i.e. complex emulated structures for which solutions can change on a time scale that is very short compared with the interval of time integration, but where the solution of interest changes on a much longer time scale. Stability conditions of the coupled substructures are analysed by means of the zero‐stability approach, and the accuracy of the novel algorithms in the coupled case is assessed in both the unforced and forced conditions. LSRT algorithms are shown to be more competitive than popular Runge–Kutta methods in terms of stability, accuracy and ease of implementation. Numerical simulations and real‐time substructure tests are used to demonstrate the favourable properties of the proposed algorithms. Copyright © 2007 John Wiley & Sons, Ltd.     

Real‐time hybrid experiments with Newmark integration,MCSmd outer‐loop control and multi‐tasking strategies

Real‐time hybrid testing is a promising technique for experimental structural dynamics, in which the structure under consideration is split into a physical test of key components and a numerical model of the remainder. The physical test and numerical analysis proceed in parallel, in real time, enabling testing of critical elements at large scale and at the correct loading rate. To date most real‐time hybrid tests have been restricted to simple configurations and have used approximate delay compensation schemes. This paper describes a real‐time hybrid testing approach in which non‐linearity is permitted in both the physical and numerical models, and in which multiple interfaces between physical and numerical substructures can be accommodated, even when this results in very stiff coupling between actuators. This is achieved using a Newmark explicit numerical solver, an advanced adaptive controller known as MCSmd and a multi‐tasking strategy. The approach is evaluated through a series of experiments on discrete mass–spring systems. Copyright © 2006 John Wiley & Sons, Ltd.     

Real‐time hybrid testing using the unconditionally stable explicit CR integration algorithm

Real‐time hybrid testing combines experimental testing and numerical simulation, and provides a viable alternative for the dynamic testing of structural systems. An integration algorithm is used in real‐time hybrid testing to compute the structural response based on feedback restoring forces from experimental and analytical substructures. Explicit integration algorithms are usually preferred over implicit algorithms as they do not require iteration and are therefore computationally efficient. The time step size for explicit integration algorithms, which are typically conditionally stable, can be extremely small in order to avoid numerical stability when the number of degree‐of‐freedom of the structure becomes large. This paper presents the implementation and application of a newly developed unconditionally stable explicit integration algorithm for real‐time hybrid testing. The development of the integration algorithm is briefly reviewed. An extrapolation procedure is introduced in the implementation of the algorithm for real‐time testing to ensure the continuous movement of the servo‐hydraulic actuator. The stability of the implemented integration algorithm is investigated using control theory. Real‐time hybrid test results of single‐degree‐of‐freedom and multi‐degree‐of‐freedom structures with a passive elastomeric damper subjected to earthquake ground motion are presented. The explicit integration algorithm is shown to enable the exceptional real‐time hybrid test results to be achieved. Copyright © 2008 John Wiley & Sons, Ltd.     

Stability analysis of SDOF real‐time hybrid testing systems with explicit integration algorithms and actuator delay

Real‐time hybrid testing is a method that combines experimental substructure(s) representing component(s) of a structure with a numerical model of the remaining part of the structure. These substructures are combined with the integration algorithm for the test and the servo‐hydraulic actuator to form the real‐time hybrid testing system. The inherent dynamics of the servo‐hydraulic actuator used in real‐time hybrid testing will give rise to a time delay, which may result in a degradation of accuracy of the test, and possibly render the system to become unstable. To acquire a better understanding of the stability of a real‐time hybrid test with actuator delay, a stability analysis procedure for single‐degree‐of‐freedom structures is presented that includes both the actuator delay and an explicit integration algorithm. The actuator delay is modeled by a discrete transfer function and combined with a discrete transfer function representing the integration algorithm to form a closed‐loop transfer function for the real‐time hybrid testing system. The stability of the system is investigated by examining the poles of the closed‐loop transfer function. The effect of actuator delay on the stability of a real‐time hybrid test is shown to be dependent on the structural parameters as well as the form of the integration algorithm. The stability analysis results can have a significant difference compared with the solution from the delay differential equation, thereby illustrating the need to include the integration algorithm in the stability analysis of a real‐time hybrid testing system. Copyright © 2007 John Wiley & Sons, Ltd.     

Equivalent force control method for generalized real‐time substructure testing with implicit integration

This paper presents a new method, called the equivalent force control method, for solving the nonlinear equations of motion in a real‐time substructure test using an implicit time integration algorithm. The method replaces the numerical iteration in implicit integration with a force‐feedback control loop, while displacement control is retained to control the motion of an actuator. The method is formulated in such a way that it represents a unified approach that also encompasses the effective force test method. The accuracy and effectiveness of the method have been demonstrated with numerical simulations of real‐time substructure tests with physical substructures represented by spring and damper elements, respectively. The method has also been validated with actual tests in which a Magnetorheological damper was used as the physical substructure. Copyright © 2007 John Wiley & Sons, Ltd.     

Operator‐splitting method for real‐time substructure testing

It has been shown that the operator‐splitting method (OSM) provides explicit and unconditionally stable solutions for quasi‐static pseudo‐dynamic substructure testing. However, the OSM provides only an explicit target displacement but not an explicit target velocity, so that it is essentially an implicit method for real‐time substructure testing (RST) when the velocity‐dependent restoring force is considered. This paper proposes a target velocity formulation based on the forward difference of the predicted displacements so as to render the OSM explicit for RST. The stability and accuracy of the resulting OSM‐RST algorithm are investigated. It is shown that the OSM‐RST is unconditionally stable so long as the non‐linear stiffness and damping are of the softening type (i.e. the tangent stiffness and damping never exceed the initial values). The stability of the OSM‐RST for structures with infinite tangent damping coefficient or stiffness is also proved, and the stability of the method for MDOF structures with a non‐classical damping matrix is demonstrated by an energy criterion. The effects of actuator delay and compensation are analysed based on the bilinear approximation of the actuator step response. Experiments on damped SDOF and MDOF structures verify that the stability of the OSM‐RST is preserved when the experimental substructure generates velocity‐dependent reaction forces, whereas the stability of real‐time substructure tests based on the central difference method is worsened by the damping of the specimen. Copyright © 2005 John Wiley & Sons, Ltd.     

Implementation of online model updating in hybrid simulation

Hybrid simulation combines numerical and experimental methods for cost‐effective, large‐scale testing of structures under simulated earthquake loading. Structural system level response can be obtained by expressing the equation of motion for the combined experimental and numerical substructures, and solved using time‐stepping integration similar to pure numerical simulations. It is often assumed that a reliable model exists for the numerical substructures while the experimental substructures correspond to parts of the structure that are difficult to model. A wealth of data becomes available during the simulation from the measured experiment response that can be used to improve upon the numerical models, particularly if a component with similar structural configuration and material properties is being tested and subjected to a comparable load pattern. To take advantage of experimental measurements, a new hybrid test framework is proposed with an updating scheme to update the initial modeling parameters of the numerical model based on the instantaneously‐measured response of the experimental substructures as the test progresses. Numerical simulations are first conducted to evaluate key algorithms for the selection and calibration of modeling parameters that can be updated. The framework is then expanded to conduct actual hybrid simulations of a structural frame model including a physical substructure in the laboratory and a numerical substructure that is updated during the tests. The effectiveness of the proposed framework is demonstrated for a simple frame structure but is extendable to more complex structural behavior and models. Copyright © 2013 John Wiley & Sons, Ltd.     

Testing of dynamically substructured,base‐isolated systems using adaptive control techniques

Experimental techniques for testing dynamically substructured systems are currently receiving attention in a wide range of structural, aerospace and automotive engineering environments. Dynamic substructuring enables full‐size, critical components to be physically tested within a laboratory (as physical substructures), while the remaining parts are simulated in real‐time (as numerical substructures). High quality control is required to achieve synchronization of variables at the substructuring interfaces and to compensate for additional actuator system(s) dynamics, nonlinearities, uncertainties and time‐varying parameters within the physical substructures. This paper presents the substructuring approach and associated controller designs for performance testing of an aseismic, base‐isolation system, which is comprised of roller‐pendulum isolators and controllable, nonlinear magnetorheological dampers. Roller‐pendulum isolators are typically mounted between the protected structure and its foundation and have a fundamental period of oscillation far‐removed from the predominant periods of any earthquake. Such semi‐active damper systems can ensure safety and performance requirements, whereas the implementation of purely active systems can be problematic in this respect. A linear inverse dynamics compensation and an adaptive controller are tailored for the resulting nonlinear synchronization problem. Implementation results favourably compare the effectiveness of the adaptive substructuring method against a conventional shaking‐table technique. A 1.32% error resulted compared with the shaking‐table response. Ultimately, the accuracy of the substructuring method compared with the response of the shaking‐table is dependent upon the fidelity of the numerical substructure. Copyright © 2009 John Wiley & Sons, Ltd.     

A Rosenbrock‐W method for real‐time dynamic substructuring and pseudo‐dynamic testing

A variant of the Rosenbrock‐W integration method is proposed for real‐time dynamic substructuring and pseudo‐dynamic testing. In this variant, an approximation of the Jacobian matrix that accounts for the properties of both the physical and numerical substructures is used throughout the analysis process. Only an initial estimate of the stiffness and damping properties of the physical components is required. It is demonstrated that the method is unconditionally stable provided that specific conditions are fulfilled and that the order accuracy can be maintained in the nonlinear regime without involving any matrix inversion while testing. The method also features controllable numerical energy dissipation characteristics and explicit expression of the target displacement and velocity vectors. The stability and accuracy of the proposed integration scheme are examined in the paper. The method has also been verified through hybrid testing performed of SDOF and MDOF structures with linear and highly nonlinear physical substructures. The results are compared with those obtained from the operator splitting method. An approach based on the modal decomposition principle is presented to predict the potential effect of experimental errors on the overall response during testing. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluation of numerical time‐integration schemes for real‐time hybrid testing

We present a comparison of methods for the analysis of the numerical substructure in a real‐time hybrid test. A multi‐tasking strategy is described, which satisfies the various control and numerical requirements. Within this strategy a variety of explicit and implicit time‐integration algorithms have been evaluated. Fully implicit schemes can be used in fast hybrid testing via a digital sub‐step feedback technique, but it is shown that this approach requires a large amount of computation at each sub‐step, making real‐time execution difficult for all but the simplest models. In cases where the numerical substructure poses no harsh stability condition, it is shown that the Newmark explicit method offers advantages of speed and accuracy. Where the stability limit of an explicit method cannot be met, one of the several alternatives may be used, such as Chang's modified Newmark scheme or the α‐operator splitting method. Appropriate methods of actuator delay compensation are also discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

Adaptive multi‐rate interface: development and experimental verification for real‐time hybrid simulation

Real‐time hybrid simulation (RTHS) is a powerful cyber‐physical technique that is a relatively cost‐effective method to perform global/local system evaluation of structural systems. A major factor that determines the ability of an RTHS to represent true system‐level behavior is the fidelity of the numerical substructure. While the use of higher‐order models increases fidelity of the simulation, it also increases the demand for computational resources. Because RTHS is executed at real‐time, in a conventional RTHS configuration, this increase in computational resources may limit the achievable sampling frequencies and/or introduce delays that can degrade its stability and performance. In this study, the Adaptive Multi‐rate Interface rate‐transitioning and compensation technique is developed to enable the use of more complex numerical models. Such a multi‐rate RTHS is strictly executed at real‐time, although it employs different time steps in the numerical and the physical substructures while including rate‐transitioning to link the components appropriately. Typically, a higher‐order numerical substructure model is solved at larger time intervals, and is coupled with a physical substructure that is driven at smaller time intervals for actuator control purposes. Through a series of simulations, the performance of the AMRI and several existing approaches for multi‐rate RTHS is compared. It is noted that compared with existing methods, AMRI leads to a smaller error, especially at higher ratios of sampling frequency between the numerical and physical substructures and for input signals with high‐frequency content. Further, it does not induce signal chattering at the coupling frequency. The effectiveness of AMRI is also verified experimentally. Copyright © 2016 John Wiley & Sons, Ltd.     

Rosenbrock‐based algorithms and subcycling strategies for real‐time nonlinear substructure testing

In this paper, Rosenbrock‐based algorithms originally developed for real‐time testing of linear systems with dynamic substructuring are extended for use on nonlinear systems. With this objective in mind and for minimal overhead, both two‐ and three‐stages linearly implicit real‐time compatible algorithms were endowed with the Jacobian matrices requiring only one evaluation at the beginning of each time step. Moreover, these algorithms were improved with subcycling strategies. In detail, the paper briefly introduces Rosenbrock‐based L‐Stable Real‐Time (LSRT) algorithms together with linearly implicit and explicit structural integrators, which are now commonly used to perform real‐time tests. Then, the LSRT algorithms are analysed in terms of linearized stability with reference to an emulated spring pendulum, which was chosen as a nonlinear test problem, because it is able to exhibit a large and relatively slow nonlinear circular motion coupled to an axial motion that can be set to be stiff. The accuracy analysis on this system was performed for all the algorithms described. Following this, a coupled spring‐pendulum example typical of real‐time testing is analysed with respect to both stability and accuracy issues. Finally, the results of representative numerical simulations and real‐time substructure tests, considering nonlinearities both in the numerical and the physical substructure, are explored. These tests were used to demonstrate how the LSRT algorithms can be used for substructuring tests with strongly nonlinear components. Copyright © 2010 John Wiley & Sons, Ltd.     

Stability of central difference method for dynamic real‐time substructure testing

This paper studies the stability of the central difference method (CDM) for real‐time substructure test considering specimen mass. Because the standard CDM is implicit in terms of acceleration, to avoid iteration, an explicit acceleration formulation is assumed for its implementation in real‐time dynamic substructure testing. The analytical work shows that the stability of the algorithm decreases with increasing specimen mass if the experimental substructure is a pure inertia specimen. The algorithm becomes unstable however small the time integration interval is, when the mass of specimen equal or greater than that of its numerical counterpart. For the case of dynamic specimen, the algorithm is unstable when there is no damping in the whole test structure; a damping will make the algorithm stable conditionally. Part of the analytical results is validated through an actual test. Copyright © 2009 John Wiley & Sons, Ltd.     

Large‐scale real‐time hybrid simulation involving multiple experimental substructures and adaptive actuator delay compensation

Real‐time hybrid simulation provides a viable method to experimentally evaluate the performance of structural systems subjected to earthquakes. The structural system is divided into substructures, where part of the system is modeled by experimental substructures, whereas the remaining part is modeled analytically. The displacements in a real‐time hybrid simulation are imposed by servo‐hydraulic actuators to the experimental substructures. Actuator delay compensation has been shown by numerous researchers to vitally achieve reliable real‐time hybrid simulation results. Several studies have been performed on servo‐hydraulic actuator delay compensation involving single experimental substructure with single actuator. Research on real‐time hybrid simulation involving multiple experimental substructures, however, is limited. The effect of actuator delay during a real‐time hybrid simulation with multiple experimental substructures presents challenges. The restoring forces from experimental substructures may be coupled to two or more degrees of freedom (DOF) of the structural system, and the delay in each actuator must be adequately compensated. This paper first presents a stability analysis of actuator delay for real‐time hybrid simulation of a multiple‐DOF linear elastic structure to illustrate the effect of coupled DOFs on the stability of the simulation. An adaptive compensation method then proposed for the stable and accurate control of multiple actuators for a real‐time hybrid simulation. Real‐time hybrid simulation of a two‐story four‐bay steel moment‐resisting frame with large‐scale magneto‐rheological dampers in passive‐on mode subjected to the design basis earthquake is used to experimentally demonstrate the effectiveness of the compensation method in minimizing actuator delay in multiple experimental substructures. Copyright © 2011 John Wiley & Sons, Ltd.     

Modified predictor–corrector numerical scheme for real‐time pseudo dynamic tests using state‐space formulation

This paper deals with an explicit numerical integration method for real‐time pseudo dynamic tests. The proposed method, termed the MPC‐SSP method, is suited to use in real‐time pseudo dynamic tests as no iteration steps are involved in each step of computation. A procedure for implementing the proposed method in real‐time pseudo dynamic tests is described in the paper. A state‐space approach is employed in this study to formulate the equations of motion of the system, which is advantageous in real‐time pseudo dynamic testing of structures with active control devices since most structural control problems are formulated in state space. A stability and accuracy analysis of the proposed method was performed based on linear elastic systems. Owing to an extrapolation scheme employed to predict the system's future response, the MPC‐SSP method is conditionally stable. To demonstrate the effectiveness of the MPC‐SSP method, a series of numerical simulations were performed and the performance of the MPC‐SSP method was compared with other pseudo dynamic testing methods including Explicit Newmark, Central Difference, Operator Splitting, and OS‐SSP methods based on both linear and non‐linear single‐degree‐of‐freedom systems. Copyright © 2004 John Wiley & Sons, Ltd.     

Stability and accuracy analysis of the central difference method for real‐time substructure testing

The central difference method (CDM) that is explicit for pseudo‐dynamic testing is also believed to be explicit for real‐time substructure testing (RST). However, to obtain the correct velocity dependent restoring force of the physical substructure being tested, the target velocity is required to be calculated as well as the displacement. The standard CDM provides only explicit target displacement but not explicit target velocity. This paper investigates the required modification of the standard central difference method when applied to RST and analyzes the stability and accuracy of the modified CDM for RST. Copyright © 2005 John Wiley & Sons, Ltd.     

Efficiency of base isolation systems in structural seismic protection and energetic assessment

This paper concerns the seismic response of structures isolated at the base by means of High Damping Rubber Bearings (HDRB). The analysis is performed by using a stochastic approach, and a Gaussian zero mean filtered non‐stationary stochastic process is used in order to model the seismic acceleration acting at the base of the structure. More precisely, the generalized Kanai–Tajimi model is adopted to describe the non‐stationary amplitude and frequency characteristics of the seismic motion. The hysteretic differential Bouc–Wen model (BWM) is adopted in order to take into account the non‐linear constitutive behaviour both of the base isolation device and of the structure. Moreover, the stochastic linearization method in the time domain is adopted to estimate the statistical moments of the non‐linear system response in the state space. The non‐linear differential equation of the response covariance matrix is then solved by using an iterative procedure which updates the coefficients of the equivalent linear system at each step and searches for the solution of the response covariance matrix equation. After the system response variance is estimated, a sensitivity analysis is carried out. The final aim of the research is to assess the real capacity of base isolation devices in order to protect the structures from seismic actions, by avoiding a non‐linear response, with associated large plastic displacements and, therefore, by limiting related damage phenomena in structural and non‐structural elements. In order to attain this objective the stochastic response of a non‐linear n‐dof shear‐type base‐isolated building is analysed; the constitutive law both of the structure and of the base devices is described, as previously reported, by adopting the BWM and by using appropriate parameters for this model, able to suitably characterize an ordinary building and the base isolators considered in the study. The protection level offered to the structure by the base isolators is then assessed by evaluating the reduction both of the displacement response and the hysteretic dissipated energy. Copyright © 2003 John Wiley & Sons, Ltd.     

Model‐based predicting and correcting algorithms for substructure online hybrid tests

A set of algorithms combined with a substructure technique is proposed for an online hybrid test framework, in which the substructures are encapsulated by a standard interface that implements displacements and forces at the common substructure boundaries. A coordinator equipped with the proposed algorithms is designed to achieve boundary compatibility and equilibrium, thereby endowing the substructures the ability to behave as one piece. A model‐based predictor and corrector, and a noniterative procedure, characterize the set of algorithms. The coordinator solves the dynamics of the entire structure and updates the static boundary state simultaneously by a quasi‐Newton procedure, which gradually formulates the condensed stiffness matrix associated with corresponding degrees of freedom. With the condensed stiffness matrix and dynamic information, a condensed equation of motion is derived and then solved by a typical time integration algorithm. Three strategies for updating the condensed stiffness matrix are incorporated into the proposed algorithms. Each adopts different stiffness matrix during the predicting and correcting stage. These algorithms are validated by two numerical substructure simulations and a hybrid test. The effectiveness and feasibility are fully demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

Hybrid simulation of steel frame structures with sectional model updating

Hybrid simulations that combine numerical computations and physical experiment represent an effective method of evaluating the dynamic response of structures. However, it is sometimes impossible to take all the uncertain or nonlinear parts of the structure as the physical substructure. Thus, the modeling errors of the numerical part can raise concerns. One method of solving this problem is to update the numerical model by estimating its parameters from experimental data online. In this paper, an online model updating method for the hybrid simulation of frame structures is proposed to reduce the errors of nonlinear modeling of numerical substructures. To obtain acceptable accuracy with acceptable extra computation efforts as a result of model parameter estimation, the sectional constitutive model is adopted, therein considering axial‐force and bending‐moment coupling; moreover, the unscented Kalman filter is used for parameter estimation of the sectional model. The effectiveness of the sectional model updating with the unscented Kalman filter is validated via numerical analyses and actual hybrid tests on a full‐scale steel frame structure, with one column as the experimental substructure loaded by three actuators to guarantee the consistency of the boundary conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Stability analysis for real‐time pseudodynamic and hybrid pseudodynamic testing with multiple sources of delay

Real‐time pseudodynamic (PSD) and hybrid PSD test methods are experimental techniques to obtain the response of structures, where restoring force feedback is used by an integration algorithm to generate command displacements. Time delays in the restoring force feedback from the physical test structure and/or the analytical substructure cause inaccuracies and can potentially destabilize the system. In this paper a method for investigating the stability of structural systems involved in real‐time PSD and hybrid PSD tests with multiple sources of delay is presented. The method involves the use of the pseudodelay technique to perform an exact mapping of fixed delay terms to determine the stability boundary. The approach described here is intended to be a practical one that enables the requirements for a real‐time testing system to be established in terms of system parameters when multiple sources of delay exist. Several real‐time testing scenarios with delay that include single degree of freedom (SDOF) and multi‐degree of freedom (MDOF) real‐time PSD/hybrid PSD tests are analyzed to illustrate the method. From the stability analysis of the real‐time hybrid testing of an SDOF test structure, delay‐independent stability with respect to either experimental or analytical substructure delay is shown to exist. The conditions that the structural properties must satisfy in order for delay‐independent stability to exist are derived. Real‐time hybrid PSD testing of an MDOF structure equipped with a passive damper is also investigated, where observations from six different cases related to the stability plane behavior are summarized. Throughout this study, root locus plots are used to provide insight and explanation of the behavior of the stability boundaries. Copyright © 2008 John Wiley & Sons, Ltd.