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.     

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.     

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.     

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.     

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.     

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.     

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.     

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

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

The ρ-family of algorithms for time-step integration with improved numerical dissipation

The dynamic analysis of complex non-linear structural systems by the finite element approach requires the use of time-step algorithms for solving the equations of motion in the time domain. Both an implicit and an explicit version of such a time-step algorithm, called the ρ-method, the parameter ρ being used for controlling numerical damping in the higher modes, are presented in this paper. For the implicit family of algorithms unconditional stability, consistency, convergence, accuracy and overshoot properties are first discussed and proved. On the basis of the algorithmic damping ratio (dissipation) and period elongation (dispersion) the ρ-method is then compared with the well-known implicit algorithms of Hilber, Newmark, Wilson, Park and Houbolt. An explicit version of the algorithm is also derived and briefly discussed. This shows numerical properties similar to the central difference method. Both versions of the algorithm have been implemented in a general purpose computer program which has been often used for both numerical tests and practical applications.     

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.     

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.     

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.     

一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用

在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.     

基于一种新显式时间积分算法的场地非线性地震反应分析

针对线弹性结构动力学方程,作者已提出一种具有良好稳定性的二阶精度单步显式时间积分算法。本文将该方法推广到求解材料非线性结构动力学方程中,采用带误差控制的修正欧拉算法计算单元应力,提高显式时间积分算法的精度。将求解非线性问题的显式算法应用于地震波垂直入射时非线性地震反应分析中,使用黏性边界模拟场地土层底部半空间基岩的辐射阻尼,并考虑地震动输入。与中心差分法计算结果进行对比,以表明新显式算法的有效性。     

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.     

An error analysis of truncated starting conditions in step-by-step time integration: Consequences for structural dynamics

Most of the step-by-step time integration algorithms for structural dynamics require an initial acceleration vector to be specified, in addition to displacement and velocity vectors. A consistent initial acceleration vector may be calculated by solving the equations of motion at the initial time, while a truncated initial acceleration vector is obtained by setting the acceleration values to zero. Although the truncated starting procedure decreases computational effort, it is shown to affect accuracy adversely. For the structural dynamics algorithms considered herein, the rate at which the numerical solution converges to the exact solution is O(Δt) when the truncated starting procedure is used, compared to O(Δt²) when consistent initial acceleration values are used.     

Velocity plus displacement equivalent force control for real-time substructure testing

This paper employs a velocity plus displacement(V+D)-based equivalent force control(EFC) method to solve the velocity/displacement difference equation in a real-time substructure test. This method uses type 2 feedback control loops to replace mathematical iteration to solve the nonlinear dynamic equation. A spectral radius analysis of the amplification matrix shows that the type 2 EFC-explicit, Newmark-β method has beneficial numerical characteristics for this method. Its stability limit of Ω = 2 remains unchanged regardless of the system damping because the velocity is achieved with very high accuracy during simulation. In contrast, the stability limits of the central difference method using direct velocity prediction and the EFC-average acceleration method with linear interpolation are shown to decrease with an increase in system damping. In fact, the EFC-average acceleration method is shown to change from unconditionally stable to conditionally stable. We also show that if an over-damped system with a damping ratio of 1.05 is considered, the stability limit is reduced to Ω =1.45. Finally, the results from an experiment with a single-degree-of-freedom structure installed with a magneto-rheological(MR) damper are presented. The results demonstrate that the proposed method is able to follow both displacement and velocity commands with moderate accuracy, resulting in improved test performance and accuracy for structures that are sensitive to both velocity and displacement inputs. Although the findings of the study are promising, additional test data and several further improvements will be required to draw general conclusions.     

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.     

Numerical failure analysis of a continuous reinforced concrete bridge under strong earthquakes using multi-scale models

Previous failure analyses of bridges typically focus on substructure failure or superstructure failure separately. However, in an actual bridge, the seismic induced substructure failure and superstructure failure may influence each other. Moreover, previous studies typically use simplified models to analyze the bridge failure; however, there are inherent defects in the calculation accuracy compared with using a detailed three-dimensional (3D) finite element (FE) model. Conversely, a detailed 3D FE model requires more computational costs, and a proper erosion criterion of the 3D elements is necessary. In this paper, a multi-scale FE model, including a corresponding erosion criterion, is proposed and validated that can significantly reduce computational costs with high precision by modelling a pseudo-dynamic test of an reinforced concrete (RC) pier. Numerical simulations of the seismic failures of a continuous RC bridge based on the multi-scale FE modeling method using LS-DYNA are performed. The nonlinear properties of the bridge, various connection strengths and bidirectional excitations are considered. The numerical results demonstrate that the failure of the connections will induce large pounding responses of the girders. The nonlinear deformation of the piers will aggravate the pounding damages. Furthermore, bidirectional earthquakes will induce eccentric poundings to the girders and different failure modes to the adjacent piers.     

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.