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.     

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.     

实时子结构实验Chang算法的稳定性和精度

与慢速拟动力子结构实验相比,实时子结构实验的优点在于它能真实地反映速度相关型试件的特性。实时子结构实验的逐步积分算法通常借用拟动力算法,但是目前液压伺服作动器很难实现速度反馈控制,因而试件速度不能实现原算法的假定值,这样一来算法的稳定性和计算精度将发生改变。台湾学者S.Y.Chang提出一种无条件稳定的显式拟动力算法,本文分析了这种方法应用于实时子结构实验时的稳定性和计算精度。研究发现在实时子结构实验中该方法由无条件稳定变成了有条件稳定的,精度也发生了改变。     

Stability of average acceleration method for structures with nonlinear damping

The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.     

A method for estimating specimen tangent stiffness for hybrid simulation

Researchers have long recognized the importance and potential benefits of utilizing the tangent stiffness matrix of a test specimen in hybrid simulations employing implicit and mixed‐integration schemes. However, the computation of the tangent stiffness matrix during testing has proved to be challenging, particularly for test specimens with more than one degree of freedom (dof). This paper presents a new methodology that is more straightforward and simpler than existing techniques for computing the tangent stiffness matrix of a multi‐dof test specimen. The proposed method is combined with the operator‐splitting method (OSM), and the capabilities, advantages and limitations of the new formulation are demonstrated through several examples. The accuracy, stability, and error propagation characteristics of the modified OSM are also studied theoretically as well as numerically. The research results show that the proposed algorithm provides results that are better than those produced via the regular OSM alone, especially for damped structures undergoing highly inelastic behavior during testing. Copyright © 2008 John Wiley & Sons, Ltd.     

Improved control algorithm for real‐time substructure testing

Real‐time substructure testing is a novel method of testing structures under dynamic loading. The complete structure is separated into two substructures, one of which is tested physically at large scale and in real time, so that time‐dependent non‐linear behaviour of the substructure is realistically represented. The second substructure represents the surrounding structure, which is modelled numerically. In the current formulation this numerical substructure is assumed to remain linear. The two substructures interact in real‐time so that the response of the complete structure, incorporating the non‐linear behaviour of the physical substructure, is accurately represented. This paper presents several improvements to the linear numerical modelling of substructures for use in explicit time‐stepping routines for real‐time substructure testing. An extrapolation of a first‐order‐hold discretization is used which increases the accuracy of the numerical model over more direct explicit methods. Additionally, an integral form of the equation of motion is used in order to reduce the effects of noise and to take into account variations of the input over a time‐step. In order to take advantage of this integral form, interpolation of the model output is performed in order to smooth the output. The improvements are demonstrated using a series of substructure tests on a simple portal frame. While the testing approach is suitable for cases in which the physical substructure behaves non‐linearly, the results presented here are for fully linear systems. This enables comparisons to be made with analytical solutions, as well as with the results of tests based on the central difference method. Copyright © 2001 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.     

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.     

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.     

一种基于离散控制理论的无条件稳定显式算法

基于离散控制理论,结合CR法和RST法提出一种无条件稳定的动力学显式新算法。以算法精度和稳定性为条件,通过离散传递函数推导参数表达式和极点,使得新算法可满足零振幅衰减率和零周期延长率。算法参数α和γ作为传递格式选择参数,当α和γ分别取1时,新算法对应CR法和RST法的位移速度表达式。对新算法的精度和稳定性理论分析表明:新算法可满足无振幅衰减和周期延长,且对于线性系统和非线性刚度软化系统为无条件稳定,对非线性刚度硬化系统为条件稳定,并给出了非线性刚度硬化系统的稳定性范围。算例分析验证了新算法的精度和稳定性,证明提出的新算法是可靠有效的。     

考虑阻尼影响的修正CD-Newmark组合算法数值特性研究

显-隐式组合数值积分算法结合了显式算法无需迭代和隐式算法无条件稳定的各自优点,是结构抗震拟动力试验顺利运行的关键.在对传统显式中央差分法和隐式Newmark β组合算法进行参数修正的基础上,建立了修正CD-Newmark算法,考虑阻尼的影响分析了组合算法的稳定性条件、周期失真率和数值阻尼比,分别得到了试验子结构的稳定性条件和计算子结构无条件稳定的参数合理取值范围,并对计算精度进行了分析.通过算例分析验证了算法的数值特性,从而初步解决了CD-Newmark算法存在稳定性界限过严的问题,为结构抗震拟动力混合试验提供了研究参考.     

Numerical artifacts associated with Rayleigh and modal damping models of inelastic structures with massless coordinates

The paper presents a detailed reexamination of the effects of three damping models on the inelastic seismic response of structures with massless degrees of freedom. The models considered correspond to (a) Rayleigh damping based on current properties (tangent stiffness), (b) Rayleigh damping based on initial properties, and (c) modal damping. The results suggest that some nonzero damping forces/moments at massless DOFs obtained in multistory frames for the case of Rayleigh damping with tangent stiffness may be numerical artifacts rather than a deficiency of the damping model. The results also indicate that significant artificial numerical oscillations in the velocities of the secondary components of MDOF structures are introduced when modal damping or mass-proportional damping is used.     

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.     

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.     

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.     

Free‐vibration analysis of a three‐dimensional soil–structure system

A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd.     

Real‐time hybrid testing of a smart base isolation system

Real‐time hybrid testing is a very effective technique for evaluating the dynamic responses of rate‐dependent structural systems subjected to earthquake excitation. A smart base isolation system has been proposed by others using conventional low‐damping isolators and controllable damping devices such as magnetorheological (MR) dampers to achieve specified control target performance. In this paper, real‐time hybrid tests of a smart base isolation system are conducted. The simulation is for a base‐isolated two‐degrees‐of‐freedom building model where the superstructure and the low‐damping base isolator are numerically simulated, and the MR damper is physically tested. The target displacement obtained from the step‐by‐step integration of the numerical substructure is imposed on the MR damper, which is driven by three different control algorithms in real‐time. To compensate the actuator delay and improve the accuracy of the test, an adaptive phase‐lead compensator is implemented. The accuracy of each test is investigated by using the root mean square error and the tracking indicator. Experimental results demonstrate that the hybrid testing procedure using the proposed actuator compensation techniques is effective for investigating the control performance of the MR damper in a smart base isolation system. Copyright © 2013 John Wiley & Sons, Ltd.     

Stability analysis of MDOF real‐time dynamic hybrid testing systems using the discrete‐time root locus technique

The time delay resulting from the servo hydraulic systems can potentially destabilize the real‐time dynamic hybrid testing (RTDHT) systems. In this paper, the discrete‐time root locus technique is adopted to investigate the delay‐dependent stability performance of MDOF RTDHT systems. Stability analysis of an idealized two‐story shear frame with two DOFs is first performed to illustrate the proposed method. The delay‐dependent stability condition is presented for various structural properties, time delay, and integration time steps. Effects of delay compensation methods on stability are also investigated. Then, the proposed method is applied to analyze the delay‐dependent stability of a single shaking table RTDHT system with an 18‐DOF finite element numerical substructure, and corresponding RTDHTs are carried out to verify the theoretical results. Furthermore, the stability behavior of a finite element RTDHT system with two physical substructures, loaded by twin shaking tables, is theoretically and experimentally investigated. All experimental results convincingly demonstrate that the delay‐dependent stability analysis on the basis of the discrete‐time root locus technique is feasible. Copyright © 2014 John Wiley & Sons, Ltd.     

Equivalent force control method for substructure pseudo‐dynamic test of a full‐scale masonry structure

The effectiveness of equivalent force control (EFC) method has been experimentally validated through hybrid tests with simple specimens. In this paper, the EFC method is applied for the MDOF pseudo‐dynamic substructure tests in which a three‐storey frame‐supported reinforced concrete masonry shear wall with full scale is chosen as physical substructure. The effects of equivalent force controller parameters on the response performance are studied. Analytical expressions for the controller parameter ranges are derived to avoid response overshooting or oscillation and are verified by numerical simulation. The controller parameters are determined based on analytical and numerical studies and used in the actual full‐scale pseudo‐dynamic test. The test results show good tracking performance of EFC, which indicates a successful test. Copyright © 2013 John Wiley & Sons, Ltd.