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.     

Modified state–space procedures for pseudodynamic testing

The existing on‐line numerical integration algorithms are derived from the Newmark method, which is based on an approximation of derivatives in the differential equation. The state–space procedure (SSP), based on an interpolation of the discrete excitation signals for piecewise convolution integral, has been confirmed as more reliable than the Newmark method in terms of numerical accuracy and stability. In an attempt to enhance the pseudodynamic test, this study presents an on‐line integration algorithm (referred to as the OS–SSP method) via an integration of the state–space procedure with Nakashima's operator‐splitting concept. Numerical stability and accuracy assessment of the proposed algorithm in addition to the explicit Newmark method and the OS method were investigated via an eigenvalue, frequency‐domain and time‐domain analysis. Of the on‐line integration algorithms investigated, the OS–SSP method is demonstrated as the most accurate method with an acceptable stability (although not unconditionally stable) characteristic. Therefore, the OS–SSP method is the most desirable method for pseudodynamic testing if the numerical stability criterion (Δt/T⩽0.5) is ensured for every vibration mode involved. Copyright © 2001 John Wiley & Sons, Ltd.     

Phase and amplitude error indices for error quantification in pseudodynamic testing

Real-time pseudodynamic (PSD) and hybrid PSD testing methods are displacement controlled experimental techniques that are used to investigate the dynamic behaviour of complex and load rate-dependent structures. Because the imposed command displacements are not predefined but generated during the test based on measured feedback, these methods are inherently prone to error propagation, which can affect the accuracy and even the stability of the entire experiment. As a result, to have these experimental methods as reliable tools, the accuracy of the test results needs to be assessed by carefully monitoring, and if possible, quantifying the errors involved. In this paper, phase and amplitude error indices (PAEI) are introduced to identify the experimental errors through uncoupled closed-form equations. Unlike the indicators that have been previously introduced in the literature for error identification purposes, PAEI do not use test setup specific parameters in their formulation, and can quantify the errors independent of the amplitude of the command displacements. As such, PAEI can be used as standard tools for assessing the quality of the experiments performed in different laboratories or under different conditions. Additionally, because they can quantify the error, when implemented online, PAEI have the potential to be incorporated in the control law and thereby improve the actuator control during the tests. The formulation and implementation of PAEI are provided in this paper. The enhanced performance of the proposed indices is demonstrated by processing several different measured and command signals using PAEI and comparing the results with those revealed by the previous indicators. Copyright © 2011 John Wiley & Sons, Ltd.     

Analytical exploration of γ-function explicit method for pseudodynamic testing of nonlinear systems

It has been well studied that the γ-function explicit method can be effective in providing favorable numerical dissipation for linear elastic systems. However, its performance for nonlinear systems is unclear due to a lack of analytical evaluation techniques. Thus, a novel technique is proposed herein to evaluate its efficiency for application to nonlinear systems by introducing two parameters to describe the stiffness change. As a result, the numerical properties and error propagation characteristics of the γ-function explicit method for the pseudodynamic testing of a nonlinear system are analytically assessed. It is found that the upper stability limit decreases as the step degree of nonlinearity increases; and it increases as the current degree of nonlinearity increases. It is also shown that this integration method provides favorable numerical dissipation not only for linear elastic systems but also for nonlinear systems. Furthermore, error propagation analysis reveals that the numerical dissipation can effectively suppress the severe error propagation of high frequency modes while the low frequency responses are almost unaffected for both linear elastic and nonlinear systems.     

Stability and accuracy analysis of outer loop dynamics in real‐time pseudodynamic testing of SDOF systems

Real‐time pseudodynamic (PSD) testing is an experimental technique for evaluating the dynamic behaviour of a complex structure. During the test, when the targeted command displacements are not achieved by the test structure, or a delay in the measured restoring forces from the test structure exists, the reliability of the testing method is impaired. The stability and accuracy of real‐time PSD testing in the presence of amplitude error and a time delay in the restoring force is presented. Systems consisting of an elastic single degree of freedom (SDOF) structure with load‐rate independent and dependent restoring forces are considered. Bode plots are used to assess the effects of amplitude error and a time delay on the steady‐state accuracy of the system. A method called the pseudodelay technique is used to derive the exact solution to the delay differential equation for the critical time delay that causes instability of the system. The solution is expressed in terms of the test structure parameters (mass, damping, stiffness). An error in the restoring force amplitude is shown to degrade the accuracy of a real‐time PSD test but not destabilize the system, while a time delay can lead to instability. Example calculations are performed for determining the critical time delay, and numerical simulations with both a constant delay and variable delay in the restoring force are shown to agree well with the stability limit for the system based on the critical time delay solution. The simulation models are also used to investigate the effects of a time delay in the PSD test of an inelastic SDOF system. The effect of energy dissipation in an inelastic structure increases the limit for the critical time delay, due to the energy removed from the system by the energy dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Selection of time step for pseudodynamic testing

Although the step degree of nonlinearity has been introduced to conduct basic analysis and error propagation analysis for the pseudodynamic testing of nonlinear systems, it cannot be reliably used to select an appropriate time step before performing a pseudodynamic test. Therefore, a novel parameter of instantaneous degree of nonlinearity is introduced to monitor the stiffness change at the end of a time step, and can thus be used to evaluate numerical and error propagation properties for nonlinear systems. Based on these properties, it is possible to select an appropriate time step to conduct a pseudodynamic test in advance. This possibility is very important in pseudodynamic testing, since the use of an arbitrary time step might lead to unreliable results or even destroy the test specimen. In this paper, guidelines are proposed for choosing an appropriate time step for accurate integration of nonlinear systems. These guidelines require estimation of the maximum instantaneous degree of nonlinearity and the solution of the initial natural frequency. The Newmark explicit method is chosen for this study. All the analytical results and the guidelines proposed are thoroughly confirmed with numerical examples.     

Improved time integration for pseudodynamic tests

Converting the second-order differential equation to a first-order equation by integrating it with respect to time once as the governing equation of motion for a structural system can be very promising in the pseudodynamic testing. This was originally found and developed by Chang. The application of this time-integration technique to the Newmark explicit method is implimented and investigated in this paper. The main advantages of using the integral form of Newmark explicit method instead of the commonly used Newmark explicit method in a pseudodynamic test are: a less-error propagation effect, a better capability in capturing the rapid changes of dynamic loading and in eliminating the adverse linearization errors. All these improvements have been verified by theoretical studies and experimental tests. Consequently, for a same time step this time-integration technique may result in less-error propagation and achieve more accurate test results than applying the original form of Newmark explicit method in a pseudodynamic test due to these significant improvements. Thus, the incorporation of this proposed time-integration technique into the direct integration method for pseudodynamic testings is strongly recommended. © 1998 John Wiley & Sons, Ltd.     

A family of explicit algorithms for general pseudodynamic testing

A new family of explicit pseudodynamic algorithms is proposed for general pseudodynamic testing. One particular subfamily seems very promising for use in general pseudodynamic testing since the stability problem for a structure does not need to be considered. This is because this subfamily is unconditionally stable for any instantaneous stiffness softening system, linear elastic system and instantaneous stiffness hardening system that might occur in the pseudodynamic testing of a real structure. In addition, it also offers good accuracy when compared to a general second-order accurate method for both linear elastic and nonlinear systems.     

Improved numerical dissipation for explicit methods in pseudodynamic tests

There is no second-order accurate, dissipative, explicit method in the currently available step-by-step integration algorithms. Two new families of second-order accurate, dissipative, explicit methods have been successfully developed for the direct integration of equations of motion in structural dynamics. These two families of methods are numerically equivalent and possess the desired numerical dissipation which can be continuously controlled. These two families of algorithms are very useful for pseudodynamic tests since the favourable numerical damping can be used to suppress the spurious growth of high-frequency modes due to the presence of numerical and/or experimental errors in performing a pseudodynamic test. © 1997 by John Wiley & Sons, Ltd.     

On the use of multi‐algorithm,genetically adaptive multi‐objective method for multi‐site calibration of the SWAT model

With the availability of spatially distributed data, distributed hydrologic models are increasingly used for simulation of spatially varied hydrologic processes to understand and manage natural and human activities that affect watershed systems. Multi‐objective optimization methods have been applied to calibrate distributed hydrologic models using observed data from multiple sites. As the time consumed by running these complex models is increasing substantially, selecting efficient and effective multi‐objective optimization algorithms is becoming a nontrivial issue. In this study, we evaluated a multi‐algorithm, genetically adaptive multi‐objective method (AMALGAM) for multi‐site calibration of a distributed hydrologic model—Soil and Water Assessment Tool (SWAT), and compared its performance with two widely used evolutionary multi‐objective optimization (EMO) algorithms (i.e. Strength Pareto Evolutionary Algorithm 2 (SPEA2) and Non‐dominated Sorted Genetic Algorithm II (NSGA‐II)). In order to provide insights into each method's overall performance, these three methods were tested in four watersheds with various characteristics. The test results indicate that the AMALGAM can consistently provide competitive or superior results compared with the other two methods. The multi‐method search framework of AMALGAM, which can flexibly and adaptively utilize multiple optimization algorithms, makes it a promising tool for multi‐site calibration of the distributed SWAT. For practical use of AMALGAM, it is suggested to implement this method in multiple trials with relatively small number of model runs rather than run it once with long iterations. In addition, incorporating different multi‐objective optimization algorithms and multi‐mode search operators into AMALGAM deserves further research. Copyright © 2009 John Wiley & Sons, Ltd.     

Multi‐directional response of unreinforced masonry walls: experimental and computational investigations

This paper describes the results of an experimental and numerical study that focused on multi‐directional behavior of unreinforced masonry walls and established the requisite of the related proposed design equations. The tests were conducted following several sets of multi‐directional loading combinations imposed on the top plane of the wall along with considering monotonic and cyclic quasi‐static loading protocols. Various boundary conditions, representing possible wall–roof connections, were also considered for different walls to investigate the influence of rotation of the top plane of the wall on the failure modes. The results of the tests were recorded with a host of high precision data acquisition systems, showing three‐dimensional displacements of a grid on the surface of the wall. Finite element models of the walls are developed using the commercial software package ABAQUS/Explicit compiled with a FORTRAN subroutine (VUMAT) written by the authors. The experimental results were then used to validate the finite element models and the developed user‐defined material models. With the utility of validated models, a parametric study was performed on a set of parameters with dominant influence on the behavior of the wall system under in‐plane and out‐of‐plane loading combinations. The experimental and numerical results are finally used to investigate the adequacy of ASCE 41 empirical equations, and some insights and recommendations are made. Copyright © 2016 John Wiley & Sons, Ltd.     

Extensions of nonlinear error propagation analysis for explicit pseudodynamic testing

Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.     

Considerations for the development of real‐time dynamic testing using servo‐hydraulic actuation

This paper presents a study of the use of servo‐hydraulic systems in the implementation of real‐time large‐scale structural testing methods in force control such as effective force testing (EFT) and in displacement control such as real‐time pseudodynamic testing (RPsD). Mathematical models for both types of control systems are presented and used to investigate the influences of servo‐systems on the overall system performance. Parameters investigated include the overall system dynamics, nonlinearities of servo‐systems, actuator damping, system mass including piston mass, and system response delay. Results of both numerical simulations and experiments showed that many of the influences of the servo‐hydraulic system that significantly affect the real‐time dynamic tests can be properly compensated through control schemes identified in this paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Model‐based multi‐metric control of uniaxial shake tables

Shake tables provide a direct means by which to evaluate structural performance under earthquake excitation. Because the entire structure is mounted on the base plate and subjected to the ground motion in real time, dynamic effects and rate‐dependent behavior can be accurately represented. Shake table control is not straightforward as the desired signal is an acceleration record, while most actuators operate in displacement feedback for stability. At the same time, the payload is typically large relative to the capacity of the actuator, leading to pronounced control‐structure interaction. Through this interaction, the dynamics of the specimen influence the dynamics of the shake table, which can be problematic when specimens change behavior because of damage or other nonlinearities. Moreover, shake tables are themselves inherently nonlinear, making it difficult to accurately recreate a desired acceleration record over a broad range of amplitudes and frequencies. A model‐based multi‐metric shake table control strategy is proposed to improve tracking of the desired acceleration of a uniaxial shake table, remaining robust to nonlinearities including changes in specimen condition. The proposed strategy is verified for the shake table testing of both linear and nonlinear structures. Copyright © 2013 John Wiley & Sons, Ltd.     

Seismic response of three‐dimensional r/c multi‐storey frame building under uni‐ and bi‐directional input ground motion

This paper deals with seismic analysis of plan‐asymmetric r/c frame multi‐storey buildings. Non‐linear numerical analyses are carried out by using a lumped plasticity model for beams and a multi‐spring model for columns, the latter one introduced to account for axial force–biaxial bending moment interaction. A comparison between numerical analyses and experimental test results is reported in order to calibrate the numerical model, showing that the adopted model is very suitable. In order to study the effects of the earthquake orthogonal component, the seismic response of the modelled structure under uni‐directional excitation is compared to the one under bi‐directional excitation. Such comparison shows that the maximum base shear and the top displacement are not very sensitive to the presence of the orthogonal component, which, conversely, leads to large increase in the column plastic excursions. 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.     

Applying a real‐coded multi‐population genetic algorithm to multi‐reservoir operation

The primary objective of this study is to propose a real‐coded hypercubic distributed genetic algorithm (HDGA) for optimizing reservoir operation system. A conventional genetic algorithm (GA) is often trapped into local optimums during the optimization procedure. To prevent premature convergence and to obtain near‐global optimal solutions, the HDGA is designed to have various subpopulations that are processed using separate and parallel GAs. The hypercubic topology with a small diameter spreads good solutions rapidly throughout all of the subpopulations, and a migration mechanism, which exchanges chromosomes among the subpopulations, exchanges information during the joint optimization to maintain diversity and thus avoid a systematic premature convergence toward a single local optimum. Three genetic operators, i.e. linear ranking selection, blend‐α crossover and Gaussian mutation, are applied to search for the optimal reservoir releases. First, a benchmark problem, the four‐reservoir operation system, is considered to investigate the applicability and effectiveness of the proposed approach. The results show that the known global optimal solution can be effectively and stably achieved by the HDGA. The HDGA is then applied in the planning of a multi‐reservoir system in northern Taiwan, considering a water reservoir development scenario to the year 2021. The results searched by an HDGA minimize the water deficit of this reservoir system and provide much better performance than the conventional GA in terms of obtaining lower values of the objective function and avoiding local optimal solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

Identification of the hysteretic behaviour of a partial‐strength steel–concrete moment‐resisting frame structure subject to pseudodynamic tests

In low‐rise steel‐concrete composite structures, moment‐resisting frames can be designed to develop a ductile response in beam‐to‐column joints and column bases by activating flexural yielding of beams and end plates, shear yielding of column web panel zones and yielding of anchors. To evaluate the performance of these components under differing earthquake intensities, a series of pseudodynamic, quasistatic cyclic and vibration tests were carried out on a two‐storey two‐bay moment resisting structure. The performance‐based seismic design and control of these structures requires that stiffness degradation, strength deterioration and slip are properly modelled. In this context, compact hysteretic models can play a key role and must therefore be striven for. Nonetheless, relevant techniques, like nonlinear system identification, are far from representing standard and reliable tools for the dynamic characterization of full‐scale structural systems. With this objective in mind, we present a restoring force surface‐based technique applied to pseudodynamic test data, in view of the nonlinear identification of multistorey frames. The technique is developed by means of a parametric approach, where a time‐variant stiffness operator is coupled to a modified Bouc–Wen model that allows both for slip and for degradation in stiffness. Strength deterioration is indirectly taken into account too. We also show how model‐based parameters can be correlated to the damage process progressively observed both in the structure and in its components. Finally, the predictive capabilities of the identified model are highlighted. Copyright © 2012 John Wiley & Sons, Ltd.     

Performance of a real‐time pseudodynamic test system considering nonlinear structural response

This paper presents the implementation details of a real‐time pseudodynamic test system that adopts an implicit time integration scheme. The basic configuration of the system is presented. Physical tests were conducted to evaluate the performance of the system and validate a theoretical system model that incorporates the dynamics and nonlinearity of a test structure and servo‐hydraulic actuators, control algorithm, actuator delay compensation methods, and the flexibility of an actuator reaction system. The robustness and accuracy of the computational scheme under displacement control errors and severe structural softening are examined with numerical simulations using the model. Different delay compensation schemes have been implemented and compared. One of the schemes also compensates for the deformation of an actuator reaction system. It has been shown that the test method is able to attain a good performance in terms of numerical stability and accuracy. However, it has been shown that test results obtained with this method can underestimate the inelastic displacement drift when severe strain softening develops in a test structure. This can be attributed to the fact that the numerical damping effect introduced by convergence errors becomes more significant as a structure softens. In a real‐time test, a significant portion of the convergence errors is caused by the time delay in actuator response. Hence, a softening structure demands higher precision in displacement control. Copyright © 2007 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.