A novel predictor–corrector algorithm for sub‐structure pseudo‐dynamic testing

A new predictor–corrector (P–C) method for multi‐site sub‐structure pseudo‐dynamic (PSD) test is proposed. This method is a mixed time integration method in which computational components separable from experimental components are solved by implicit time integration method (Newmark β method). The experiments are performed quasi‐statically based on explicit prediction of displacement. The proposed P–C method has an important advantage as it does not require the determination of the initial stiffness values of experimental components and is thus suitable for representing elastic and inelastic systems. A parameter relating to quality of displacement prediction at boundaries nodes is introduced. This parameter is determined such that P–C method can be applicable to many practical problems. Error‐propagation characteristics of P–C method are also presented. A series of examples including linear and non‐linear soil–foundation–structure interaction problem demonstrate the performance of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

Full operator algorithm for hybrid simulation

One of the weaknesses of the operator splitting method (OSM) is that its corrector step employs the approximation that incremental forces are linearly related to the tested structure's initial stiffness matrix. This paper presents a new predictor–corrector technique in which the assumptions about the tested structure's response are shifted to the predictor step, which results in an enhancement in overall simulation accuracy, especially for nonlinear structures. Unlike OSM, which splits the displacement and velocity operators into explicit and implicit terms, the new method uses predicted accelerations to compute fully explicit displacement and velocity values in the predictor step. Another advantage of the proposed technique, termed the full operator method (FOM) is that its formulation makes it suitable for both quasi‐static and real‐time hybrid simulation. The effectiveness of FOM is first evaluated by investigating error propagation in an undamped single degree‐of‐freedom model. It is shown that the corrector step in FOM is able to significantly suppress aberrant simulation results caused by incorrect estimation of the structure's stiffness matrix. The performance of FOM is demonstrated by exercising two additional models, which exhibit significant inelastic behavior under the prescribed excitation. The simulation results show that the proposed FOM algorithm is capable of producing accurate solutions and that the corrector step is influential in effectively reducing simulation errors. It is also shown that FOM suppresses actuator displacement control errors because of its reliance on measured quantities in the corrector step. Copyright © 2009 John Wiley & Sons, Ltd.     

Pseudo‐dynamic testing with substructuring at the ELSA Laboratory

The paper describes the distinctive features of the pseudo‐dynamic test method as implemented at the ELSA reaction‐wall facility. Both hardware and software aspects are considered. Particular attention is devoted to the digital control system and to a coupled numerical–experimental substructuring technique allowing realistic earthquake testing of very large structures. Mathematical and implementation details corresponding to this testing technique are given for both synchronous and asynchronous input motion. Selected test results illustrate the advantages of the presented features. Copyright © 2000 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.     

Bond‐graph based substructuring of dynamical systems

A bond graph approach to hybrid simulation of dynamical systems using numerical–experimental real‐time substructuring is presented. The bond graph concepts of a virtual junction and a virtual actuator, hitherto used in the context of physical‐model based control, are used to perform the substructuring in an intuitively appealing way. The approach is illustrated by the reworking of a previously‐published example. The approach is verified experimentally using a bench‐top multiple mass–spring system for the physical substructure and automatically generated real‐time code is used to implement the numerical substructure. Copyright © 2005 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.     

Seismic behavior of bridges with pier and foundation strengthening: PsD tests and analytical study

By means of a simplified three degrees of freedom model, seismic behavior of reinforced concrete bridge piers and foundations were evaluated based on pseudo‐dynamic (PsD) tests for cases where pier strengthening and foundation strengthening are implemented. In addition, analysis based on PsD test results was conducted to investigate the influence of pier strengthening on seismic damage to the foundation. The PsD tests and the analysis show that the foundation suffers increased hysteretic response when pier strengthening is applied. The results also show that the foundation strengthening can prevent foundation damage. Copyright © 2011 John Wiley & Sons, Ltd.     

Online hybrid test by internet linkage of distributed test‐analysis domains

Online hybrid tests (called the online tests), particularly when combined with substructuring techniques, are able to conduct large‐scale tests. An extension of this technique is to combine multiple loading tests conducted in remote locations and to integrate the tests with large numerical analysis codes. In this study, a new Internet online test system is developed in which a physical test is conducted in one place, the associated numerical analysis is performed in a remote location, and the two locations communicate over the Internet. To implement the system, a technique that links test and analysis domains located at different places is proposed, and an Internet data exchange interface is devised to allow data communication across Internet. A practical method that utilizes standard protocols implemented by operating systems for sharing files and folders is adopted to ensure stable and robust communication between remotely located servers that commonly protect themselves by strict firewalls. To combine the online test with a finite element program formulated in an incremental form and adopting an implicit integration scheme, a tangent stiffness prediction procedure is proposed. In this procedure, a tangent stiffness is estimated based on a few previous steps of experimental data. Using the system devised, tests on a base‐isolated structure were carried out. Here, the base‐isolation layer was taken as the tested part and tested in Kyoto University, Japan, and the superstructure was modelled by means of a finite element program and analysed in a computer located in Osaka University. A series of physical Internet online tests were carried out, with the integration time interval and the method of tangent stiffness prediction as the major parameters. The tests demonstrated that the Internet communication was very stable and robust, without malfunctions. The proposed method of stiffness prediction was effective even when the experimental hysteresis curves exhibited complex behaviour, thereby ensuring accurate simulation for the earthquake response of the entire structure. Copyright © 2005 John Wiley & Sons, Ltd.     

Non‐iterative time integration scheme for non‐linear dynamic FEM analysis

This paper proposes a non‐iterative time integration (NITI) scheme for non‐linear dynamic FEM analysis. The NITI scheme is constructed by combining explicit and implicit schemes, taking advantage of their merits, and enables stable computation without an iteration process for convergence even when used for non‐linear dynamic problems. Formulation of the NITI scheme is presented and its stability is studied. Although the NITI scheme is not unconditionally stable when applied to non‐linear problems, it is stable in most cases unless stiffness hardening occurs or the problem has a large velocity‐dependent term. The NITI scheme is applied to dynamic analysis of the non‐linear soil–structure system and computation results are compared with those by the central difference method (CDM). Comparison shows that the stability of the NITI scheme is superior to that of the CDM. Accuracy of the NITI scheme is verified because its results are identical with those by the CDM in which the time step is set as 1/10 of that for the NITI scheme. The application of the NITI scheme to the mesh‐partitioned FEM is also proposed. It is applied to dynamic analysis of the linear soil–structure system. It yields the same results as a conventional single‐domain FEM analysis using the Newmark β method. This result verifies the usability of mesh‐partitioned FEM analysis using the NITI scheme. Copyright © 2003 John Wiley& Sons, Ltd.     

Implementation and application of the unconditionally stable explicit parametrically dissipative KR‐α method for real‐time hybrid simulation

In real‐time hybrid simulations (RTHS) that utilize explicit integration algorithms, the inherent damping in the analytical substructure is generally defined using mass and initial stiffness proportional damping. This type of damping model is known to produce inaccurate results when the structure undergoes significant inelastic deformations. To alleviate the problem, a form of a nonproportional damping model often used in numerical simulations involving implicit integration algorithms can be considered. This type of damping model, however, when used with explicit integration algorithms can require a small time step to achieve the desired accuracy in an RTHS involving a structure with a large number of degrees of freedom. Restrictions on the minimum time step exist in an RTHS that are associated with the computational demand. Integrating the equations of motion for an RTHS with too large of a time step can result in spurious high‐frequency oscillations in the member forces for elements of the structural model that undergo inelastic deformations. The problem is circumvented by introducing the parametrically controllable numerical energy dissipation available in the recently developed unconditionally stable explicit KR‐α method. This paper reviews the formulation of the KR‐α method and presents an efficient implementation for RTHS. Using the method, RTHS of a three‐story 0.6‐scale prototype steel building with nonlinear elastomeric dampers are conducted with a ground motion scaled to the design basis and maximum considered earthquake hazard levels. The results show that controllable numerical energy dissipation can significantly eliminate spurious participation of higher modes and produce exceptional RTHS results. Copyright © 2014 John Wiley & Sons, Ltd.     

Verification of system identification utilizing shaking table tests of a full‐scale 4‐story steel building

This paper verifies the feasibility of the proposed system identification methods by utilizing shaking table tests of a full‐scale four‐story steel building at E‐Defense in Japan. The natural frequencies, damping ratios and modal shapes are evaluated by single‐input‐four‐output ARX models. These modal parameters are prepared to identify the mass, damping and stiffness matrices when the objective structure is modelled as a four degrees of freedom (4DOF) linear shear building in each horizontal direction. The nonlinearity in stiffness is expressed as a Bouc–Wen hysteretic system when it is modelled as a 4DOF nonlinear shear building. The identified hysteretic curves of all stories are compared to the corresponding experimental results. The simple damage detection is implemented using single‐input‐single‐output ARX models, which require only two measurements in each horizontal direction. The modal parameters are equivalent‐linearly evaluated by the recursive Least Squares Method with a forgetting factor. When the structure is damaged, its natural frequencies decrease, and the corresponding damping ratios increase. The fluctuation of the identified modal properties is the indirect information for damage detection of the structure. Copyright © 2015 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.     

Substructuring of dynamical systems via the adaptive minimal control synthesis algorithm

In this paper we consider the concept of modelling dynamical systems using numerical–experimental substructuring. This type of modelling is applicable to large or complex systems, where some part of the system is difficult to model numerically. The substructured model is formed via the adaptive minimal control synthesis (MCS) algorithm. The aim of this paper is to demonstrate that substructuring can be carried out in real time, using the MCS algorithm. Thus, we reformulate the MCS algorithm into a substructuring form. We introduce the concepts of a transfer system, and carry out numerical simulations of the substructuring process using a coupled three mass example. These simulations are compared with direct simulations of a three mass system. In addition we consider the stability of the substructuring algorithm, which we discuss in detail for a class of second‐order transfer systems. A numerical–experimental system is considered, using a small‐scale experimental system, for which the substructuring algorithm is implemented in real time. Finally we discuss these results, with particular reference to the future application of this method to modelling large‐scale structures subject to earthquake excitation. Copyright © 2001 John Wiley & Sons, Ltd.     

Control issues relating to real‐time substructuring experiments using a shaking table

Real‐time substructuring is a method of dynamically testing a structure without experimentally testing a physical model of the entire system. Instead the structure can be split into two linked parts, the region of particular interest, which is tested experimentally, and the remainder which is tested numerically. A transfer system, such as a hydraulic actuator or a shaking table, is used to impose the displacements at the interface between the two parts on the experimental substructure. The corresponding force imposed by the substructure on the transfer system is fed back to the numerical model. Control of the transfer system is critical to the accuracy of the substructuring process. A study of two controllers used in conjunction with the University of Bristol shaking table is presented here. A proof‐of‐concept one degree‐of‐freedom mass–spring–damper system is substructured such that a portion of the mass forms the experimental substructure and the remainder of the mass plus the spring and the damper is modelled numerically. Firstly a linear controller is designed and tested. Following this an adaptive substructuring strategy is considered, based on the minimal control synthesis algorithm. The deleterious effect of oil‐column resonance common to shaking tables is examined and reduced through the use of filters. The controlled response of the experimental specimen is compared for the two control strategies. Copyright © 2005 John Wiley & Sons, Ltd.     

Comparison between real‐time dynamic substructuring and shake table testing techniques for nonlinear seismic applications

Results from real‐time dynamic substructuring (RTDS) tests are compared with results from shake table tests performed on a two‐storey steel building structure model. At each storey, the structural system consists of a cantilevered steel column resisting lateral loads in bending. In two tests, a slender diagonal tension‐only steel bracing member was added at the first floor to obtain an unsymmetrical system with highly variable stiffness. Only the first‐storey structural components were included in the RTDS test program and a Rosenbrock‐W linearly implicit integration scheme was adopted for the numerical solution. The tests were performed under seismic ground motions exhibiting various amplitude levels and frequency contents to develop first and second mode‐dominated responses as well as elastic and inelastic responses. A chirp signal was also used. Coherent results were obtained between the shake table and the RTDS testing techniques, indicating that RTDS testing methods can be used to successfully reproduce both the linear and nonlinear seismic responses of ductile structural steel seismic force resisting systems. The time delay introduced by actuator‐control systems was also studied and a novel adaptive compensation scheme is proposed. Copyright © 2010 John Wiley & Sons, Ltd.     

Development of a family of unconditionally stable explicit direct integration algorithms with controllable numerical energy dissipation

The implicit dissipative generalized‐ α method is analyzed using discrete control theory. Based on this analysis, a one‐parameter family of explicit direct integration algorithms with controllable numerical energy dissipation, referred to as the explicit KR‐α method, is developed for linear and nonlinear structural dynamic numerical analysis applications. Stability, numerical dispersion, and energy dissipation characteristics of the proposed algorithms are studied. It is shown that the algorithms are unconditionally stable for linear elastic and stiffness softening‐type nonlinear systems, where the latter indicates a reduction in post yield stiffness in the force–deformation response. The amount of numerical damping is controlled by a single parameter, which provides a measure of the numerical energy dissipation at higher frequencies. Thus, for a specific value of this parameter, the resulting algorithm is shown to produce no numerical energy dissipation. Furthermore, it is shown that the influence of the numerical damping on the lower mode response is negligible. It is further shown that the numerical dispersion and energy dissipation characteristics of the proposed explicit algorithms are the same as that of the implicit generalized‐ α method. A numerical example is presented to demonstrate the potential of the proposed algorithms in reducing participation of undesired higher modes by using numerical energy dissipation to damp out these modes. Copyright © 2014 John Wiley & Sons, Ltd.     

A component approach for the hysteretic behaviour of connections in cross‐laminated wooden structures

In this paper, a numerical model to estimate the dissipative capacity and describe the cyclic response of cross‐laminated (X‐lam) timber buildings is presented. The connections between panels and to the foundation (metal hold‐downs and angle brackets, and screwed connectors) are modelled with nonlinear hysteretic multispring elements taking into account the strength interaction between different degrees of freedom according to a predefined domain. The timber components (solid X‐lam floors and wall panels) are modelled using elastic shell elements. By calibration on experimental cyclic tests carried out on each degree of freedom, important features of timber connection behaviour such as post‐peak strength, pinching and stiffness degradation can all be considered. In addition, the effect of friction at the interface between panels and with foundation can be taken into account. These springs have been implemented as external subroutines in a widespread software package such as Abaqus. By comparison with the experimental results of cyclic tests carried out on single X‐lam walls, coupled X‐lam walls and a single‐storey X‐lam building, the accuracy of the proposed model is demonstrated. Copyright © 2013 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.     

A technique for the reduction of dynamic degrees of freedom

This paper presents a variation on the component mode technique for the dynamic substructuring of large‐scale structural analysis of building and bridge frames. The principal innovation of the proposed method of dynamic reduction is that the resulting mass matrix of the reduced substructures remains diagonal. As in the component mode technique, the reduction is accomplished by transforming the degrees of freedom in the substructure using boundary shapes and internal shapes. The diagonal mass matrix is achieved by orthogonalization of the boundary shapes to the internal shapes, and a selective row‐by‐row summation of the mass matrix into the diagonal entry (where off‐diagonal terms are unavoidable). To aid in recovering the accuracy of the rigid‐body inertias that is lost in the diagonalization process, additional pseudo‐rigid‐body‐mode shapes are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

Amplification system for concentrated and distributed energy dissipation devices

Recent analytic, experimental, and practical studies are developing energy dissipation devices combined with amplifying mechanisms (AM) to enhance structural behavior. This research presents the theoretical and experimental development of the eccentric lever‐arm system (ELAS), a new system generically called amplified added damping (AAD), which is a combination of an AM with one or more dampers capable of supporting large deformations. The proposed AM device is a variant of the well‐known lever‐arm system. This work is divided in four parts: (1) kinematics of the ELAS and definition of an equivalent AAD; (2) parametric analysis of a linear single‐story structure with ELAS; (3) numerical analysis of a multi‐degree‐of‐freedom structure with frictional damping with and without AM; and (4) pseudo‐dynamic tests of a full scale asymmetric one story steel structure with and without frictional AAD. Parametric analyses demonstrate that using high‐amplification ratios and low supplemental damping could be a good practice. On the other hand, similar to systems without AMs, dissipation efficiency increases conformably with the stiffness of the secondary structure. As expected, it was observed that deformation was highly concentrated in the flexible edge of the asymmetric test model without damper. Conversely, the structure with frictional AAD clearly showed uniform plane deformation. The implemented AM, which has a large amplifying ratio of α≈11, performed with close accordance with numerical simulations and a high mechanical efficiency (≈95%) using a frictional damper with a very low force capacity. Copyright © 2016 John Wiley & Sons, Ltd.