Dynamic phasor-based hybrid simulation for multi-inverter grid-connected system

To realize the efficient transient simulation of a grid-connected power generation system based on multiple inverters, this paper proposes a hybrid simulation method integrating the models of electromagnetic transient and dynamic phasors. Based on a demonstration of the concepts and properties of dynamic phasors, the models of single-phase and three-phase inverters described by dynamic phasors are established first. Considering the numerical compatibility problem between dynamic phasors and instantaneous values, an interface scheme between dynamic phasors and instantaneous values is designed, and the efficiency and precision differences of various transformation methods are compared in detail.Finally, by utilizing MATLAB/Simulink, a hybrid simulation platform of a multi-inverter grid-connected system is built, and the efficiency and accuracy of the hybrid simulation are validated via comparison with the full electromagnetic transient simulation.

Evaluation of integration methods for hybrid simulation of complex structural systems through collapse

This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration （iterative） and the operator-splitting （non-iterative） is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom （SDOF） cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1：2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1：2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.

Model-based framework for multi-axial real-time hybrid simulation testing

Real-time hybrid simulation is an efficient and cost-effective dynamic testing technique for performance evaluation of structural systems subjected to earthquake loading with rate-dependent behavior. A loading assembly with multiple actuators is required to impose realistic boundary conditions on physical specimens. However, such a testing system is expected to exhibit significant dynamic coupling of the actuators and suffer from time lags that are associated with the dynamics of the servo-hydraulic system, as well as control-structure interaction （CSI）. One approach to reducing experimental errors considers a multi-input, multi-output （MIMO） controller design, yielding accurate reference tracking and noise rejection. In this paper, a framework for multi-axial real-time hybrid simulation （maRTHS） testing is presented. The methodology employs a real-time feedback-feedforward controller for multiple actuators commanded in Cartesian coordinates. Kinematic transformations between actuator space and Cartesian space are derived for all six-degrees-of- freedom of the moving platform. Then, a frequency domain identification technique is used to develop an accurate MIMO transfer function of the system. Further, a Cartesian-domain model-based feedforward-feedback controller is implemented for time lag compensation and to increase the robustness of the reference tracking for given model uncertainty. The framework is implemented using the 1/5th-scale Load and Boundary Condition Box （LBCB） located at the University of Illinois at Urbana- Champaign. To demonstrate the efficacy of the proposed methodology, a single-story frame subjected to earthquake loading is tested. One of the columns in the fraane is represented physically in the laboratory as a cantilevered steel column. For real- time execution, the numerical substructure, kinematic transformations, and controllers are implemented on a digital signal processor. Results show excellent performance of the maRTHS framework when six-degrees-of-freedom are controUed at the interface between substructures.

Restoring force correction based on online discrete tangent stiffness estimation method for real-time hybrid simulation

In real-time hybrid simulation（RTHS）, it is difficult if not impossible to completely erase the error in restoring force due to actuator response delay using existing displacement-based compensation methods. This paper proposes a new force correction method based on online discrete tangent stiffness estimation（online DTSE） to provide accurate online estimation of the instantaneous stiffness of the physical substructure. Following the discrete curve parameter recognition theory, the online DTSE method estimates the instantaneous stiffness mainly through adaptively building a fuzzy segment with the latest measurements, constructing several strict bounding lines of the segment and calculating the slope of the strict bounding lines, which significantly improves the calculation efficiency and accuracy for the instantaneous stiffness estimation. The results of both computational simulation and real-time hybrid simulation show that：（1） the online DTSE method has high calculation efficiency, of which the relatively short computation time will not interrupt RTHS; and（2） the online DTSE method provides better estimation for the instantaneous stiffness, compared with other existing estimation methods. Due to the quick and accurate estimation of instantaneous stiffness, the online DTSE method therefore provides a promising technique to correct restoring forces in RTHS.

Analysis of actuator delay and its effect on uncertainty quantification for real-time hybrid simulation

Uncertainties in structure properties can result in different responses in hybrid simulations. Quantification of the effect of these tmcertainties would enable researchers to estimate the variances of structural responses observed from experiments. This poses challenges for real-time hybrid simulation （RTHS） due to the existence of actuator delay. Polynomial chaos expansion （PCE） projects the model outputs on a basis of orthogonal stochastic polynomials to account for influences of model uncertainties. In this paper, PCE is utilized to evaluate effect of actuator delay on the maximum displacement from real-time hybrid simulation of a single degree of freedom （SDOF） structure when accounting for uncertainties in structural properties. The PCE is first applied for RTHS without delay to determine the order of PCE, the number of sample points as well as the method for coefficients calculation. The PCE is then applied to RTHS with actuator delay. The mean, variance and Sobol indices are compared and discussed to evaluate the effects of actuator delay on uncertainty quantification for RTHS. Results show that the mean and the variance of the maximum displacement increase linearly and exponentially with respect to actuator delay, respectively. Sensitivity analysis through Sobol indices also indicates the influence of the single random variable decreases while the coupling effect increases with the increase of actuator delay.

Real-time hybrid simulation for structural control performance assessment

Real-time hybrid simulation is an attractive method to evaluate the response of structures under earthquake loads. The method is a variation of the pseudodynamic testing technique in which the experiment is executed in real time, thus allowing investigation of structural systems with rate-dependent components. Real-time hybrid simulation is challenging because it requires performance of all calculations, application of displacements, and acquisition of measured forces, within a very small increment of time. Furthermore, unless appropriate compensation for actuator dynamics is implemented, stability problems are likely to occur during the experiment. This paper presents an approach for real-time hybrid simulation in which compensation for actuator dynamics is implemented using a model-based feedforward compensator. The method is used to evaluate the response of a semi-active control of a structure employing an MR damper. Experimental results show good agreement with the predicted responses, demonstrating the effectiveness of the method for structural control performance assessment.

Real-time hybrid simulation of structures equipped with viscoelastic-plastic dampers using a user-programmable computational platform

A user-programmable computational/control platform was developed at the University of Toronto that offers real-time hybrid simulation （RTHS） capabilities. The platform was verified previously using several linear physical substructures. The study presented in this paper is focused on further validating the RTHS platform using a nonlinear viscoelastic-plastic damper that has displacement, frequency and temperature-dependent properties. The validation study includes damper component characterization tests, as well as RTHS of a series of single-degree-of-freedom （SDOF） systems equipped with viscoelastic-plastic dampers that represent different structural designs. From the component characterization tests, it was found that for a wide range of excitation frequencies and friction slip loads, the tracking errors are comparable to the errors in RTHS of linear spring systems. The hybrid SDOF results are compared to an independently validated thermal- mechanical viscoelastic model to further validate the ability for the platform to test nonlinear systems. After the validation, as an application study, nonlinear SDOF hybrid tests were used to develop performance spectra to predict the response of structures equipped with damping systems that are more challenging to model analytically. The use of the experimental performance spectra is illustrated by comparing the predicted response to the hybrid test response of 2DOF systems equipped with viscoelastic-plastic dampers.

Comparison of delay compensation methods for real-time hybrid simulation using frequency-domain evaluation index

The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.

Efficiency analysis of numerical integrations for finite element substructure in real-time hybrid simulation

Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.

A new integration scheme for application to seismic hybrid simulation

Hybrid simulation is a powerful test method for evaluating the seismic performance of structural systems. This method makes it feasible that only critical components of a structure be experimentally tested. This paper presents a newly proposed integration algorithm for seismic hybrid simulation which is aimed to extend its capabilities to a wide range of systems where existing methods encounter some limitations. In the proposed method, which is termed the variable time step （VTS） integration method, an implicit scheme is employed for hybrid simulation by eliminating the iterative phase on experimental element, the phase which is necessary in regular implicit applications. In order to study the effectiveness of the VTS method, a series of numerical investigations are conducted which show the successfulness of the VTS method in obtaining accurate, stable and converged responses. Then, in a comparative approach, the improved accuracy of the VTS method over commonly used integration methods is demonstrated. The stability of the VTS method is also studied and the results show that it provides conditional stability; however, its stability limit is well beyond the accuracy limit. The effect of time delay on the VTS method results is also investigated and it is shown that the VTS method is quite successful in handling this experimental error.

Numerical characteristics of the full operator hybrid simulation method

The full operator method （FOM） has been proposed to overcome some of the shortcomings of the commonly used operator splitting method （OSM）. In particular, the FOM is improved by increasing the accuracy of both the predictor and corrector using the estimated tangent stiffness of the tested structure. The numerical characteristics of the FOM, including stability and accuracy, are investigated in this study. It is shown that FOM is conditionally stable. The stability and accuracy characteristics are dependent on the accuracy of the estimated tangent stiffness and the parameters associated with the acceleration variation in the time-stepping integration method. Mass-spring systems with different types of nonlinearity, including hardening, stiffening, and softening behavior, are used to evaluate the performance of the FOM. It is found that the FOM can capture these types of nonlinearity with satisfactory accuracy. Using a prototype 12-story composite coupled wall system, the influences of the strong nonlinearity of the system as well as the displacement control errors from hydraulic actuators on the performance of the FOM are explored. The results show that the FOM is capable of generating reasonably accurate results despite the presence of strong structural nonlinearity and displacement control errors.

A phased approach to enable hybrid simulation of complex structures

Hybrid simulation has been shown to be a cost-effective approach for assessing the seismic performance of structures. In hybrid simulation,critical parts of a structure are physically tested,while the remaining portions of the system are concurrently simulated computationally,typically using a finite element model. This combination is realized through a numerical time-integration scheme,which allows for investigation of full system-level responses of a structure in a cost-effective manner. However,conducting hybrid simulation of complex structures within large-scale testing facilities presents significant challenges. For example,the chosen modeling scheme may create numerical inaccuracies or even result in unstable simulations; the displacement and force capacity of the experimental system can be exceeded; and a hybrid test may be terminated due to poor communication between modules(e.g.,loading controllers,data acquisition systems,simulation coordinator). These problems can cause the simulation to stop suddenly,and in some cases can even result in damage to the experimental specimens; the end result can be failure of the entire experiment. This study proposes a phased approach to hybrid simulation that can validate all of the hybrid simulation components and ensure the integrity largescale hybrid simulation. In this approach,a series of hybrid simulations employing numerical components and small-scale experimental components are examined to establish this preparedness for the large-scale experiment. This validation program is incorporated into an existing,mature hybrid simulation framework,which is currently utilized in the Multi-Axial Full-Scale Sub-Structuring Testing and Simulation(MUST-SIM) facility of the George E. Brown Network for Earthquake Engineering Simulation(NEES) equipment site at the University of Illinois at Urbana-Champaign. A hybrid simulation of a four-span curved bridge is presented as an example,in which three piers are experimentally controlled in a total of 18 degrees of freedom(DOFs). This simulation illustrates the effectiveness of the phased approach presented in this paper.

Linear hybrid simulations of low-frequency fishbone instability driven by energetic passing particles in tokamak plasmas

A linear simulation study of energetic passing particle-driven low-frequency fishbone instability in tokamak plasmas has been carried out using the global kinetic-MHD(magnetohydrodynamics)hybrid code M3D-K.This work is focused on the interaction of energetic passing beam ions and n=1 mode with a monotonic safety factor q profile and q_(0)<1.Specifically,the stability and mode frequency as well as mode structure of the n=1mode are calculated for scans of parameter values of beam ion beta,beam ion injection energy,beam ion orbit width,beam ion beta profile,as well as background plasma beta.The excited modes are identified as a low-frequency fishbone with the corresponding resonance of w_(φ)+w_(θ)=w,where w_(φ)is the beam ion toroidal transit frequency and w_(θ)is the beam ion poloidal transit frequency.The simulated mode frequency is approximately proportional to the beam ion injection energy and beam ion orbit width.The mode structure is similar to that of internal kink mode.These simulation results are similar to the analytic theory of Yu et al.

Small-scale multi-axial hybrid simulation of a shear-critical reinforced concrete frame

This study presents a numerical multi-scale simulation framework which is extended to accommodate hybrid simulation （numerical-experimental integration）. The framework is enhanced with a standardized data exchange format and connected to a generalized controller interface program which facilitates communication with various types of laboratory equipment and testing configurations. A small-scale experimental program was conducted using a six degree-of-freedom hydraulic testing equipment to verify the proposed framework and provide additional data for small-scale testing of shear- critical reinforced concrete structures. The specimens were tested in a multi-axial hybrid simulation manner under a reversed cyclic loading condition simulating earthquake forces. The physical models were 1/3.23-scale representations of a beam and two columns. A mixed-type modelling technique was employed to analyze the remainder of the structures. The hybrid simulation results were compared against those obtained from a large-scale test and finite element analyses. The study found that if precautions are taken in preparing model materials and if the shear-related mechanisms are accurately considered in the numerical model, small-scale hybrid simulations can adequately simulate the behaviour of shear-critical structures. Although the findings of the study are promising, to draw general conclusions additional test data are required.

Theoretical and experimental studies on critical time delay of multi-DOF real-time hybrid simulation

This paper aims to investigate the critical stability of a multi-degree-of-freedom(multi-DOF)real-time hybrid simulation(RTHS).First,the critical time-delay analysis models are developed using the continuous-and discrete-time root locus(RL)techniques,respectively.A bilinear transform is introduced into the first-order Padéapproximation while conducting the discrete RL analysis.Based on this technique,the time delay can be explicitly used as the gain factor and thus the instability mechanism of the multi-DOF RTHS system can be analyzed.Subsequently,the critical time delays calculated by the continuous-and discrete-time RL techniques,respectively,are compared for a 2-DOF RTHS system.It is shown that assuming the RTHS system to be a continuous-time system will result in overestimating the critical time delay.Finally,theoretically calculated critical delays are demonstrated and validated by numerical simulation and a set of RTHS experiments.Parametric analysis provides a glimpse of the effects of time step,frequency and damping ratio in a performing partitioning scheme.The constructed analysis model proves to be useful for evaluating the critical time delay to predict stability and performance,therefore facilitating successful RTHS.

Seismic performance evaluation of VCFPB isolated storage tank using real-time hybrid simulation

Variable curvature friction pendulum bearings(VCFPB)effectively reduce the dynamic response of storage tanks induced by earthquakes.Shaking table testing is used to assess the seismic performance of VCFPB isolated storage tanks.However,the vertical pressure and friction coefficient of the scaled VCFPB in the shaking table tests cannot match the equivalent values of these parameters in the prototype.To avoid this drawback,a real-time hybrid simulation(RTHS)test was developed.Using RTHS testing,a 1/8 scaled tank isolated by VCFPB was tested.The experimental results showed that the displacement dynamic magnification factor of VCFPB,peak reduction factors of the acceleration,shear force,and overturning moment at bottom of the tank,were negative exponential functions of the ratio of peak ground acceleration(PGA)and friction coefficient.The peak reduction factors of displacement,acceleration,force and overturning moment,which were obtained from the experimental results,are compared with those calculated by the Housner model.It can be concluded that the Housner model is applicable in estimation of the acceleration,shear force,and overturning moment of liquid storage tank,but not for the sliding displacement of VCFPBs.

Parametric Instabilities of Parallel Propagating Circularly Polarized Alfven Waves:One-Dimensional Hybrid Simulations

By performing one-dimensional （l-D） hybrid simulations, we analyze in detail the parametric instabilities of the Alfv^n waves with a spectrum in a low beta plasma. The parametric instabilities experience two stages. In the first stage, the density modes are excited and immediately couple with the pump Alfv6n waves. In the second stage, each pump Alfv^n wave decays into a density mode and a daughter Alfv6n mode similar to the monochromatic cases. Ftlrthermore, the proton velocity beam will also be formed after the saturation of the parametric instabilities. When the plasma beta is high, the parametric decay in the second stage will be strongly suppressed.

Hybrid Simulation of ±500 kV HVDC Power Transmission Project Based on Advanced Digital Power System Simulator

In order to effectively imitate the dynamic operation characteristics of the HVDC （high voltage direct current） power transmission system at a real ±500kV HVDC transmission project, the electromechanical-electromagnetic transient hybrid simulation was carried out based on advanced digital power system simulator （ADPSS）. In the simulation analysis, the built hybrid model＇s dynamic response outputs under three different fault conditions are considered, and by comparing with the selected fault recording waveforms, the validities of the simulation waveforms are estimated qualitatively. It can be ascertained that the hybrid simulation model has the ability to describe the HVDC system＇s dynamic change trends well under some special fault conditions.

Broadband ground motion simulation using a hybrid approach of the May 21, 2021 M7.4 earthquake in Maduo, Qinghai, China

In this study,the broadband ground motions of the 2021 M7.4 Maduo earthquake were simulated to overcome the scarcity of ground motion recordings and the low resolution of macroseismic intensity map in sparsely populated high-altitude regions.The simulation was conducted with a hybrid methodology,combining a stochastic high-frequency simulation with a low-frequency ground motion simulation,from the regional 1-D velocity structure model and the Wang WM et al.(2022)source rupture model,respectively.We found that the three-component waveforms simulated for specific stations matched the waveforms recorded at those stations,in terms of amplitude,duration,and frequency content.The validation results demonstrate the ability of the hybrid simulation method to reproduce the main characteristics of the observed ground motions for the 2021 Maduo earthquake over a broad frequency range.Our simulations suggest that the official map of macroseismic intensity tends to overestimate shaking by one intensity unit.Comparisons of simulations with empirical ground motion models indicate generally good consistency between the simulated and empirically predicted intensity measures.The high-frequency components of ground motions were found to be more prominent,while the low-frequency components were not,which is unexpected for large earthquakes.Our simulations provide valuable insight into the effects of source complexity on the level and variability of the resulting ground motions.The acceleration and velocity time histories and corresponding response spectra were provided for selected representative sites where no records were available.The simulated results have important implications for evaluating the performance of engineering structures in the epicentral regions of this earthquake and for estimating seismic hazards in the Tibetan regions where no strong ground motion records are available for large earthquakes.

Validation of the current and pressure coupling schemes with nonlinear simulations of TAE and analysis on the linear stability of tearing mode in the presence of energetic particles

Both current and pressure coupling schemes have been adopted in the hybrid kinetic–magnetohydrodynamic code CLT-K recently.Numerical equivalences between these two coupling schemes are strictly verified under different approximations.First,when considering only the perturbed distribution function of energetic particles(EPs),the equivalence can be proved analytically.Second,when both the variations of the magnetic field and the EP distribution function are included,the current and pressure coupling schemes numerically produce the same result in the nonlinear simulations.On this basis,the influences of co-/counter-passing and trapped EPs on the linear stabilities of the m/n=2/1 tearing mode(TM)have been investigated(where m and n represent the poloidal and toroidal mode numbers,respectively).The results of scanningβh of EPs show that the co-passing and trapped EPs are found to stabilize the TM,while the counter-passing EPs tend to destabilize the TM.The behind(de)stabilization mechanisms of the TM by EPs are carefully analyzed.Furthermore,after exceeding critical EP betas,the same branch of the high-frequency mode is excited by co-/counterpassing and trapped EPs,which is identified as the m/n=2/1 energetic particle mode.