# Multi-Input-Multi-Output Continuous Swept-Sine Vibration Test Realization by Inverse Multistep Prediction Model

Authors: Wei Zheng; Huaihai Chen; Zhengbo Luo
Journal: Shock and Vibration
 (2020)

---

## Abstract

As frequency-varying sine excitations in rotating machines are always emerging during run-ups and shutdowns, the multi-input-multi-output (MIMO) swept-sine test is of utter significance in product validation. At present, swept-sine vibration tests are mostly conducted with frequency-domain methods, where drive spectra are generated and updated by frequency response function (FRF), and drive signals are then generated with sinusoid oscillators. In this paper, a time-domain approach using an inverse system method based on a multistep prediction model is developed to realize the MIMO continuous swept-sine vibration test. First, the multistep prediction model of the original system is estimated in the time domain. Then, the inverse multistep prediction model is derived. After that, this model is truncated to guarantee the robustness of the inverse system and the smoothness of the generated drive signals. At last, the proposed method is validated by a simulation example with a cantilever beam and an actual test by using a three-axis shaker. The results show that the MIMO continuous swept-sine vibration test can be operated effectively by the proposed method.


---

## Body

## 1. Introduction

Rotating machines employed in helicopters, propeller-driven airplanes, washing machines, and so on constantly produce swept-sine excitations to the whole structures during their run-ups and shutdowns. The vibrations induced by swept-sine excitations must be taken into consideration for the safety and the durability of the affected structures [1]. Besides, the swept-sine test is also required in some standard during product development [2] and is often utilized in modal analysis and other experimental practices [3]. It is worth noting that the swept-sine test, as it oftentimes produces cleaner frequency response function (FRF) than the random test, is frequently performed in aircraft ground vibration tests [4, 5]. For better simulation of complex cases such as aircraft with multiple propellers, a number of methods for the multi-input-multi-output (MIMO) test, including the MIMO swept-sine test, are introduced and developed during the last decades [6–8].The MIMO swept-sine test applies frequency-varying sinusoid excitations upon a test article and controls the response spectra at certain points to be within the acceptable limits. As a simplified substitution for the MIMO swept-sine test, the MIMO stepped-sine test is sometimes used in literatures [9, 10] and is often included in commercial vibration controllers. In the MIMO stepped-sine test, the frequency does not vary smoothly but instead steps through a sequence of fixed frequencies. In each step, the controlled responses should be pure sine signals with predetermined cycles. Therefore, discontinuities will occur when the frequency jumps between any two successive steps. This discontinuity can be mitigated but cannot be completely removed by some smoothing techniques [10, 11]. In fact, the stepped-sine test is essentially replacing sweeping process with multiple steady-state sinusoidal tests, which is not really the continuous swept-sine test. Several other frequency-domain methods for the MIMO swept-sine test were given in the literatures. A patented frequency-domain control method with the tracking filters for the sine test was proposed by Underwood [6], where the controlled response is achieved through iterations between learning and control loops. Yin et al. [12] presented an offline swept-sine control method, which generated drive signals by the inverse of frequency response function (FRF), and the drive signals were updated according to the errors in every control loop where the error spectra were obtained by fast Fourier transform (FFT). However, applying FFT to swept-sine signals might introduce some distortions at the beginning and end of the spectra as mentioned in the literature [9]. Besides, Fourier transform is more suitable for processing steady-state signals rather than nonstationary data such as continuous MIMO swept-sine responses.It is noted that the FRF is needed in these frequency-domain methods to connect the input and the output. However, theoretically this relationship is based on the assumption that both the input and output are of the same fixed frequency and in steady state. As this assumption is no longer valid in the swept-sine test, the relationship between the input and output should not be constructed by FRF in the frequency domain.Hence, it seems more reasonable to control the swept-sine test by the time-domain method, even though the frequency-domain methods had been successfully used for years in engineering. There seems to be no examples of the MIMO swept-sine test controlled by the time-domain method in the existing literature, to the best of the authors’ knowledge.Among time-domain approaches, the inverse system method is suitable to solve the problem at hand, as it directly takes the targeted system responses as its input to generate drive signals. This type of method is initially developed and applied in the force estimation area. As a typical inverse system method, the inverse structure filter (ISF) which utilized the inverse of the original system’s state-space model was introduced by Kammer and Steltzner [13, 14]. However, as mentioned in the literature [15], the inverse system created by the ISF is likely to be unstable, which might restrain its practicability. Regarding these disadvantages, a delayed multistep ISF method was proposed by Allen and Carne. They observed that their method could produce more stable inverse systems, though they did not prove it theoretically [15, 16]. Mao et al. [17] presented a time-domain force estimation method based on Markov parameter precise computation, and a regularization technique was utilized to deal with the ill-posed problem.Naturally, using the inverse system for drive signals generation in the vibration test is different from force estimation, where forces are estimated in a relatively short period of time as in most of the studies mentioned previously. Several recent studies on the MIMO random test conducted by Chen et al. [18] and Zheng et al. [19, 20] have taken advantage of inverse system methods. However, drive signals for random vibration are generated and smoothed in a frame-by-frame manner, and hence the initial state at each frame can be ignored, which is not the case in the swept-sine vibration test. So, the inverse system methods developed for the random vibration test cannot be used in the swept-sine vibration test directly. To address this issue, an inverse system method is developed for the MIMO swept-sine vibration test.The main contents of the paper are arranged as follows. In Section2, basic formulae for the swept-sine, multistep prediction model and some other background theories are reviewed. After that, an inverse multistep prediction model is proposed, and its application in the MIMO swept-sine test is studied. In Section 3, a numerical example of the MIMO swept-sine test by a cantilever beam is demonstrated. In Section 4, the proposed method is further validated by a test with a 3-axis shaker.


## 2. MIMO Swept-Sine Test Method

### 2.1. Theoretical Background

#### 2.1.1. Swept Sine

The MIMO swept-sine vibration test imparts frequency-varying sinusoid excitations upon a test article and controls the response spectra at certain points. At these control points, response signals should match the predefined reference swept-sine signals with some tolerable errors within the given margins.These reference swept-sine signals are given in the form of envelopeAt and transient phase φt as shown in equation (1), which are controlled by parameters such as starting frequency F1, ending frequency F2, sweep mode (linear or logarithm), and sweep rate αorβ or sweep duration (T).(1)pt=Atsinφt.The change rate of phaseφt is the instant circular frequency ωt:(2)φ˙t=ωt=2πft,where ft is the instant frequency in Hz. How this frequency varies is determined by the sweep mode.(3)ft=F1+αt,linear sweep,F12βt/60,logarithmic sweep,where α=F2−F1/T is the linear sweep rate in Hz/s and β=60log2F2/F1/T is the logarithmic sweep rate in oct/min. Substituting equation (3) into equation (2) and integrating the latter with zero condition yields(4)φt=2πF1t+12αt2linear sweep,2π60F1βln22βt/60−1logarithmic sweep.For the MIMO swept-sine vibration test, the amplitudeA and the initial phase Φ for a certain frequency are defined in the reference spectra. Since the instant frequency is determined at any given time as shown in equation (3), both the amplitude A and the initial phase Φ are also functions of t. So, the reference signal of the ith control point can be obtained by(5)pit=Aitsinφt+Φit,where Ai and Φii=1,2,…,N are the ith element of A and Φ, respectively. Hence, the reference signals for the MIMO swept-sine vibration test can be written in the vector form:(6)pt=Aitsinφt+Φit.As mentioned in the introduction, stepped sine is often used as a simplified replacement of continuous swept sine, where a number of frequencies are defined prior to the test according to each time steptii=0,…,n−1, with n being the number of steps. At each frequency, a fixed frequency sinusoid signal is determined by substituting t in equation (6) with ti. However, the number of steps has a strong influence on the effectiveness of the test: if too few steps are chosen, the dynamic characteristics of the article under test might not be fully excited; too many steps, on the other hand, may cause the period for each frequency to be very brief, thus preventing the stable responses to emerge.Therefore the stepped-sine test is flawed despite its simplicity. The key problem is that the stepped-sine test cannot achieve continuous sweep frequency and hence cannot reflect the actual continuous change process of rotating machinery. This paper focuses on solving this problem.


#### 2.1.2. Finite Difference Model

The finite difference model is a representation of the vibration system which describes the relationship between the current responses and previous drives and responses. If we letni be the number of drive channels and no be the number of response channels, the finite difference model can be written as(7)yk=−α10yk−1−⋯−αp0yk−p+β00uk+β10uk−1+⋯+βp0uk−p,where the coefficient matrices αi0i=1,…,p of dimension no×no are associated with the no×1 output vectors y at time steps k−1,…,k−p, while no×ni matrices βi0i=0,1,…,p correspond to the ni×1 input vectors u at time steps k,k−1,…,k−p. In addition, parameter p indicates the length of time history used to predict the current responses.


#### 2.1.3. Multistep Prediction Model

From the finite difference model, the multistep prediction model can be derived. Shifting equation (7) by one time step yields(8)yk+1=−α10yk−⋯−αp0yk−p+1+β00uk+1+β10uk+⋯+βp0uk−p+1.Substituting equation (7) into equation (8), one can obtain(9)yk+1=−α11yk−1−⋯−αp−11yk−p+1−αp1yk−p+β00uk+1+β01uk+β11uk−1+⋯+βp−11uk−p+1+βp1uk−p.Keep shifting and substituting equation (7) into the new one to yield(10)yk+j=−α1jyk−1−⋯−αp−1jyk−p+1−αpjyk−p+β00uk+j+β01uk+j−1+⋯+β0juk+β1juk−1+⋯+βp−1juk−p+1+βpjuk−p.The coefficient matrices can be computed by(11)α1j=α2j−1−α1j−1α10,⋮αp−1j=αpj−1−α1j−1αp−10,αpj=−α1j−1αp0,(12)β0j=β1j−1−α1j−1β00,⋮βp−1j=βpj−1−α1j−1βp−10,βpj=−α1j−1βp0,where j=1,2,…,q−1 and q represents the number of time steps of future response to be predicted. From equation (10), one can get(13)yqk=Tuqk+Bupk−p−Aypk−p,where(14)yqk=ykyk+1⋮yk+q−1,uqk=ukuk+1⋮uk+q−1,(15)ypk−p=yk−pyk−p+1⋮yk−1,upk−p=uk−puk−p+1⋮uk−1,(16)T=β000⋯0β01β00⋯0⋮⋮⋱⋮β0q−1β0q−2⋯β00,(17)A=αp0αp−10⋯α10αp1αp−11⋯α11⋮⋮⋱⋮αpq−1αp−1q−1⋯α1q−1,(18)B=βp0βp−10⋯β10βp1βp−11⋯β11⋮⋮⋱⋮βpq−1βp−1q−1⋯β1q−1.Equation (13) is also known as the multistep prediction model [21], which predicts future output signals from the previous input and output signals.The initial coefficientsαi0i=1,…,p and βi0i=0,1,…,p can be identified from test data. After that, αij and βij for any j from 1 to q−1 can be computed using equations (11) and (12), and thus the coefficient matrices T, A, and B can be obtained by equations (16)–(18).


#### 2.1.4. Truncated SVD

T matrix in equation (13) can be decomposed by SVD (singular value decomposition) into three matrices U,Σ, and V as(19)T=UΣVT,where the left and the right matrices U∈Rqno×qno and V∈Rqni×qni are orthogonal and Σ∈Rqno×qni is diagonal.(20)Σ=diagσ1,σ2,…,σk,where k=minqno,qni, and the singular values are arranged in a descending order:(21)σ1≥σ2≥⋯≥σk≥0.For actual vibration systems, some of the singular values might be very close to zero, and the direct inverse of the T matrix might induce numerical instability, so the truncated SVD or TSVD method [22] is often used to obtain a more stable inverse. In the TSVD method, an approximation of the T matrix is produced by(22)Tr=UΣrVT,Σr=diagσ1,…,σr,0,…,0,where r is the number of reserved singularity values which are larger than a given threshold. The pseudoinverse of Tr matrix can be written as(23)Tr+=VΣr+UT,Σr+=diagσ1−1,…,σr−1,0,…,0.Then,Tr+ can be used as a pseudoinverse to preclude the ill-posed problem, and it can be seen that(24)Tr+T=Ir000.When using this method, the parameterr will influence the accuracy of the pseudoinverse matrix, thus affecting the result of drive signal generation. This parameter can be determined by cross validation methods [22].


### 2.2. A New Test Method

From the multistep prediction model of the test system, the inverse multistep prediction model can be derived and used to generate the drive signals. Combining the inverse multistep prediction model with the control strategy based on matrix power, a new time-domain MIMO sweep vibration test method can be constructed.


#### 2.2.1. Inverse Multistep Prediction Model

If the T matrix in equation (13) were always invertible, then the inverse model could be obtained in a quite straight forward way, which is unfortunately not always true for test systems. Therefore, it is necessary to use pseudoinverse derived in equation (23) instead. Left multiply equation (13) by Tr+ on both sides to yield(25)Tr+yqk=Tr+Tuqk+Tr+Bupk−p−Tr+Aypk−p=Ir000uqk+Tr+Bupk−p−Tr+Aypk−p.It can be seen from equation (25) that only the first r elements of uqk are selected. However, as we mentioned in equation (14), uqk comprises drive force data from uk to uk+q−1, so the first r elements might not be meaningful in practice. Therefore, in order to select the drive force data of the length l from uk to uk+l−1, a new T¯ is formed by picking up the first l⋅ni row l⋅ni<r of Tr+. Substituting Tr+ in equation (25) with T¯ yields a new equation:(26)ulk=T¯yqk+B¯ypk−p−A¯upk−p,where yq,yp, and up are the same as in equation (14) and (15) and(27)ulk=ukuk+1⋮uk+l−1,(28)T¯=Il⋅ni0Tr+,B¯=T¯A,A¯=T¯B,in which T¯ is a l⋅ni×q⋅no matrix, A¯ is a l⋅ni×p⋅ni matrix, and B¯ is a l⋅ni×p⋅no matrix. As mentioned before, p is the number of time steps of past drives and responses and q is the number of time steps of future responses. Besides, the new parameter l stands for the time steps of future drives the inverse system generates.Needless to say, equation (26) bears remarkable resemblance to the multistep prediction model in equation (13), except that the T¯ matrix is truncated, which precludes the instability in the inverse system caused by the singularity problem while computing the inverse of T. Hence, the model defined in equation (26) is named as the inverse multistep prediction model hereafter and will be used in the MIMO swept-sine test.


#### 2.2.2. Drive Signal Generation

In the real world, test articles start vibrating from rest conditions which cannot achieve stable response immediately; a segment of transient response will always exist. The reference signals for the swept-sine test actually specify the stable response of the system, so if the reference signals are used as the input for the inverse multistep prediction model directly, the drive signals generated will be unrealistic. Hence, a transition segment should be added to the forepart of the original reference signals in order to generate reasonable drive signals as illustrated in Figure1.Figure 1
Modification of the original reference signal.To compose this transition segment, a segment of fixed frequency sine signal which matches the initial amplitude and phase of the original reference signal is generated, and a 1/4 sine window is applied to ensure the amplitude of this segment varies smoothly from zero to the same initial amplitude of the original reference signal. This transition segment is inserted before the original reference signal to modify the reference signal. Obviously, the response of this transition segment is redundant and should be ignored.Since equation (26) produces ulk instead of uqk, using the reference signals with a length of q can only generate drive signals with a length of l, where l < q. Hence, the reference signals will be used in an overlapped manner in order to generate full length drive signals. In the drive signal generating process, if the reference signals are not enough in the last segment, zeros will be supplemented. The process of drive signal generation is illustrated in Figure 2 (for clarity, the figure shows only one channel, but this process is applicable for MIMO systems).Figure 2
Drive signal is generated by the inverse system method and the reference signal is used in an overlapped manner.From the above discussions, the generation of drive signals for the MIMO swept-sine test can be concluded as follows (also illustrated in Figure3):(a)
Estimate the finite difference model of the testing system(b)
ObtainT, A, and B in equation (13) from the finite difference model(c)
Calculate the pseudoinverse of T by the TSVD method to getTr+ and then to form T¯ from Tr+(d)
Establish the inverse multistep prediction model as the inverse system using equation (28)(e)
Modify the reference signals and use them as the inputs of the inverse system to produce drive signalsFigure 3
Generation of drive signals for the MIMO swept-sine test.


#### 2.2.3. Control Strategy

Due to the noises in the input and output signals and other errors in the modeling of the system, the responses of the first loop may not be satisfactory without control. An offline control strategy is put forward, in which the response spectra are obtained and compared to the reference counterparts, and thus the errors in frequency domain are obtained. Then, the reference spectra for next control loop are updated according to these errors, and the drive signals are updated afterwards. To avoid damaging the test article and equipment, it advisable to control the swept-sine responses within acceptable range at a lower level and then restore the original level test.The reference spectra is updated as(29)Enω=R0iωSyiω,Rn+1ω=EniωεRniω,i=1,2,…,no,n=0,1,2,…where Syω is the response spectra at present loop, Rnω is the reference spectra after n times of modification (R0 is the actual reference spectra), Enω is the spectral error at present loop, Rn+1ω is the modified reference spectrum for the next loop, and ε is a positive constant and ε<1 [23]. All these spectra take the form of vectors, and they are used in equation (29) in an element-by-element manner, namely, Syiω,Eni, andRniωi=1,2,…,no represent the ith elemnt of the corresponding vectors.The control scheme for the swept-sine test is shown in Figure4.Figure 4
Offline swept-sine control scheme where A andΦ are amplitude and initial phase vectors defined in reference spectra.


## 2.1. Theoretical Background

### 2.1.1. Swept Sine

The MIMO swept-sine vibration test imparts frequency-varying sinusoid excitations upon a test article and controls the response spectra at certain points. At these control points, response signals should match the predefined reference swept-sine signals with some tolerable errors within the given margins.These reference swept-sine signals are given in the form of envelopeAt and transient phase φt as shown in equation (1), which are controlled by parameters such as starting frequency F1, ending frequency F2, sweep mode (linear or logarithm), and sweep rate αorβ or sweep duration (T).(1)pt=Atsinφt.The change rate of phaseφt is the instant circular frequency ωt:(2)φ˙t=ωt=2πft,where ft is the instant frequency in Hz. How this frequency varies is determined by the sweep mode.(3)ft=F1+αt,linear sweep,F12βt/60,logarithmic sweep,where α=F2−F1/T is the linear sweep rate in Hz/s and β=60log2F2/F1/T is the logarithmic sweep rate in oct/min. Substituting equation (3) into equation (2) and integrating the latter with zero condition yields(4)φt=2πF1t+12αt2linear sweep,2π60F1βln22βt/60−1logarithmic sweep.For the MIMO swept-sine vibration test, the amplitudeA and the initial phase Φ for a certain frequency are defined in the reference spectra. Since the instant frequency is determined at any given time as shown in equation (3), both the amplitude A and the initial phase Φ are also functions of t. So, the reference signal of the ith control point can be obtained by(5)pit=Aitsinφt+Φit,where Ai and Φii=1,2,…,N are the ith element of A and Φ, respectively. Hence, the reference signals for the MIMO swept-sine vibration test can be written in the vector form:(6)pt=Aitsinφt+Φit.As mentioned in the introduction, stepped sine is often used as a simplified replacement of continuous swept sine, where a number of frequencies are defined prior to the test according to each time steptii=0,…,n−1, with n being the number of steps. At each frequency, a fixed frequency sinusoid signal is determined by substituting t in equation (6) with ti. However, the number of steps has a strong influence on the effectiveness of the test: if too few steps are chosen, the dynamic characteristics of the article under test might not be fully excited; too many steps, on the other hand, may cause the period for each frequency to be very brief, thus preventing the stable responses to emerge.Therefore the stepped-sine test is flawed despite its simplicity. The key problem is that the stepped-sine test cannot achieve continuous sweep frequency and hence cannot reflect the actual continuous change process of rotating machinery. This paper focuses on solving this problem.


### 2.1.2. Finite Difference Model

The finite difference model is a representation of the vibration system which describes the relationship between the current responses and previous drives and responses. If we letni be the number of drive channels and no be the number of response channels, the finite difference model can be written as(7)yk=−α10yk−1−⋯−αp0yk−p+β00uk+β10uk−1+⋯+βp0uk−p,where the coefficient matrices αi0i=1,…,p of dimension no×no are associated with the no×1 output vectors y at time steps k−1,…,k−p, while no×ni matrices βi0i=0,1,…,p correspond to the ni×1 input vectors u at time steps k,k−1,…,k−p. In addition, parameter p indicates the length of time history used to predict the current responses.


### 2.1.3. Multistep Prediction Model

From the finite difference model, the multistep prediction model can be derived. Shifting equation (7) by one time step yields(8)yk+1=−α10yk−⋯−αp0yk−p+1+β00uk+1+β10uk+⋯+βp0uk−p+1.Substituting equation (7) into equation (8), one can obtain(9)yk+1=−α11yk−1−⋯−αp−11yk−p+1−αp1yk−p+β00uk+1+β01uk+β11uk−1+⋯+βp−11uk−p+1+βp1uk−p.Keep shifting and substituting equation (7) into the new one to yield(10)yk+j=−α1jyk−1−⋯−αp−1jyk−p+1−αpjyk−p+β00uk+j+β01uk+j−1+⋯+β0juk+β1juk−1+⋯+βp−1juk−p+1+βpjuk−p.The coefficient matrices can be computed by(11)α1j=α2j−1−α1j−1α10,⋮αp−1j=αpj−1−α1j−1αp−10,αpj=−α1j−1αp0,(12)β0j=β1j−1−α1j−1β00,⋮βp−1j=βpj−1−α1j−1βp−10,βpj=−α1j−1βp0,where j=1,2,…,q−1 and q represents the number of time steps of future response to be predicted. From equation (10), one can get(13)yqk=Tuqk+Bupk−p−Aypk−p,where(14)yqk=ykyk+1⋮yk+q−1,uqk=ukuk+1⋮uk+q−1,(15)ypk−p=yk−pyk−p+1⋮yk−1,upk−p=uk−puk−p+1⋮uk−1,(16)T=β000⋯0β01β00⋯0⋮⋮⋱⋮β0q−1β0q−2⋯β00,(17)A=αp0αp−10⋯α10αp1αp−11⋯α11⋮⋮⋱⋮αpq−1αp−1q−1⋯α1q−1,(18)B=βp0βp−10⋯β10βp1βp−11⋯β11⋮⋮⋱⋮βpq−1βp−1q−1⋯β1q−1.Equation (13) is also known as the multistep prediction model [21], which predicts future output signals from the previous input and output signals.The initial coefficientsαi0i=1,…,p and βi0i=0,1,…,p can be identified from test data. After that, αij and βij for any j from 1 to q−1 can be computed using equations (11) and (12), and thus the coefficient matrices T, A, and B can be obtained by equations (16)–(18).


### 2.1.4. Truncated SVD

T matrix in equation (13) can be decomposed by SVD (singular value decomposition) into three matrices U,Σ, and V as(19)T=UΣVT,where the left and the right matrices U∈Rqno×qno and V∈Rqni×qni are orthogonal and Σ∈Rqno×qni is diagonal.(20)Σ=diagσ1,σ2,…,σk,where k=minqno,qni, and the singular values are arranged in a descending order:(21)σ1≥σ2≥⋯≥σk≥0.For actual vibration systems, some of the singular values might be very close to zero, and the direct inverse of the T matrix might induce numerical instability, so the truncated SVD or TSVD method [22] is often used to obtain a more stable inverse. In the TSVD method, an approximation of the T matrix is produced by(22)Tr=UΣrVT,Σr=diagσ1,…,σr,0,…,0,where r is the number of reserved singularity values which are larger than a given threshold. The pseudoinverse of Tr matrix can be written as(23)Tr+=VΣr+UT,Σr+=diagσ1−1,…,σr−1,0,…,0.Then,Tr+ can be used as a pseudoinverse to preclude the ill-posed problem, and it can be seen that(24)Tr+T=Ir000.When using this method, the parameterr will influence the accuracy of the pseudoinverse matrix, thus affecting the result of drive signal generation. This parameter can be determined by cross validation methods [22].


## 2.1.1. Swept Sine

The MIMO swept-sine vibration test imparts frequency-varying sinusoid excitations upon a test article and controls the response spectra at certain points. At these control points, response signals should match the predefined reference swept-sine signals with some tolerable errors within the given margins.These reference swept-sine signals are given in the form of envelopeAt and transient phase φt as shown in equation (1), which are controlled by parameters such as starting frequency F1, ending frequency F2, sweep mode (linear or logarithm), and sweep rate αorβ or sweep duration (T).(1)pt=Atsinφt.The change rate of phaseφt is the instant circular frequency ωt:(2)φ˙t=ωt=2πft,where ft is the instant frequency in Hz. How this frequency varies is determined by the sweep mode.(3)ft=F1+αt,linear sweep,F12βt/60,logarithmic sweep,where α=F2−F1/T is the linear sweep rate in Hz/s and β=60log2F2/F1/T is the logarithmic sweep rate in oct/min. Substituting equation (3) into equation (2) and integrating the latter with zero condition yields(4)φt=2πF1t+12αt2linear sweep,2π60F1βln22βt/60−1logarithmic sweep.For the MIMO swept-sine vibration test, the amplitudeA and the initial phase Φ for a certain frequency are defined in the reference spectra. Since the instant frequency is determined at any given time as shown in equation (3), both the amplitude A and the initial phase Φ are also functions of t. So, the reference signal of the ith control point can be obtained by(5)pit=Aitsinφt+Φit,where Ai and Φii=1,2,…,N are the ith element of A and Φ, respectively. Hence, the reference signals for the MIMO swept-sine vibration test can be written in the vector form:(6)pt=Aitsinφt+Φit.As mentioned in the introduction, stepped sine is often used as a simplified replacement of continuous swept sine, where a number of frequencies are defined prior to the test according to each time steptii=0,…,n−1, with n being the number of steps. At each frequency, a fixed frequency sinusoid signal is determined by substituting t in equation (6) with ti. However, the number of steps has a strong influence on the effectiveness of the test: if too few steps are chosen, the dynamic characteristics of the article under test might not be fully excited; too many steps, on the other hand, may cause the period for each frequency to be very brief, thus preventing the stable responses to emerge.Therefore the stepped-sine test is flawed despite its simplicity. The key problem is that the stepped-sine test cannot achieve continuous sweep frequency and hence cannot reflect the actual continuous change process of rotating machinery. This paper focuses on solving this problem.


## 2.1.2. Finite Difference Model

The finite difference model is a representation of the vibration system which describes the relationship between the current responses and previous drives and responses. If we letni be the number of drive channels and no be the number of response channels, the finite difference model can be written as(7)yk=−α10yk−1−⋯−αp0yk−p+β00uk+β10uk−1+⋯+βp0uk−p,where the coefficient matrices αi0i=1,…,p of dimension no×no are associated with the no×1 output vectors y at time steps k−1,…,k−p, while no×ni matrices βi0i=0,1,…,p correspond to the ni×1 input vectors u at time steps k,k−1,…,k−p. In addition, parameter p indicates the length of time history used to predict the current responses.


## 2.1.3. Multistep Prediction Model

From the finite difference model, the multistep prediction model can be derived. Shifting equation (7) by one time step yields(8)yk+1=−α10yk−⋯−αp0yk−p+1+β00uk+1+β10uk+⋯+βp0uk−p+1.Substituting equation (7) into equation (8), one can obtain(9)yk+1=−α11yk−1−⋯−αp−11yk−p+1−αp1yk−p+β00uk+1+β01uk+β11uk−1+⋯+βp−11uk−p+1+βp1uk−p.Keep shifting and substituting equation (7) into the new one to yield(10)yk+j=−α1jyk−1−⋯−αp−1jyk−p+1−αpjyk−p+β00uk+j+β01uk+j−1+⋯+β0juk+β1juk−1+⋯+βp−1juk−p+1+βpjuk−p.The coefficient matrices can be computed by(11)α1j=α2j−1−α1j−1α10,⋮αp−1j=αpj−1−α1j−1αp−10,αpj=−α1j−1αp0,(12)β0j=β1j−1−α1j−1β00,⋮βp−1j=βpj−1−α1j−1βp−10,βpj=−α1j−1βp0,where j=1,2,…,q−1 and q represents the number of time steps of future response to be predicted. From equation (10), one can get(13)yqk=Tuqk+Bupk−p−Aypk−p,where(14)yqk=ykyk+1⋮yk+q−1,uqk=ukuk+1⋮uk+q−1,(15)ypk−p=yk−pyk−p+1⋮yk−1,upk−p=uk−puk−p+1⋮uk−1,(16)T=β000⋯0β01β00⋯0⋮⋮⋱⋮β0q−1β0q−2⋯β00,(17)A=αp0αp−10⋯α10αp1αp−11⋯α11⋮⋮⋱⋮αpq−1αp−1q−1⋯α1q−1,(18)B=βp0βp−10⋯β10βp1βp−11⋯β11⋮⋮⋱⋮βpq−1βp−1q−1⋯β1q−1.Equation (13) is also known as the multistep prediction model [21], which predicts future output signals from the previous input and output signals.The initial coefficientsαi0i=1,…,p and βi0i=0,1,…,p can be identified from test data. After that, αij and βij for any j from 1 to q−1 can be computed using equations (11) and (12), and thus the coefficient matrices T, A, and B can be obtained by equations (16)–(18).


## 2.1.4. Truncated SVD

T matrix in equation (13) can be decomposed by SVD (singular value decomposition) into three matrices U,Σ, and V as(19)T=UΣVT,where the left and the right matrices U∈Rqno×qno and V∈Rqni×qni are orthogonal and Σ∈Rqno×qni is diagonal.(20)Σ=diagσ1,σ2,…,σk,where k=minqno,qni, and the singular values are arranged in a descending order:(21)σ1≥σ2≥⋯≥σk≥0.For actual vibration systems, some of the singular values might be very close to zero, and the direct inverse of the T matrix might induce numerical instability, so the truncated SVD or TSVD method [22] is often used to obtain a more stable inverse. In the TSVD method, an approximation of the T matrix is produced by(22)Tr=UΣrVT,Σr=diagσ1,…,σr,0,…,0,where r is the number of reserved singularity values which are larger than a given threshold. The pseudoinverse of Tr matrix can be written as(23)Tr+=VΣr+UT,Σr+=diagσ1−1,…,σr−1,0,…,0.Then,Tr+ can be used as a pseudoinverse to preclude the ill-posed problem, and it can be seen that(24)Tr+T=Ir000.When using this method, the parameterr will influence the accuracy of the pseudoinverse matrix, thus affecting the result of drive signal generation. This parameter can be determined by cross validation methods [22].


## 2.2. A New Test Method

From the multistep prediction model of the test system, the inverse multistep prediction model can be derived and used to generate the drive signals. Combining the inverse multistep prediction model with the control strategy based on matrix power, a new time-domain MIMO sweep vibration test method can be constructed.


### 2.2.1. Inverse Multistep Prediction Model

If the T matrix in equation (13) were always invertible, then the inverse model could be obtained in a quite straight forward way, which is unfortunately not always true for test systems. Therefore, it is necessary to use pseudoinverse derived in equation (23) instead. Left multiply equation (13) by Tr+ on both sides to yield(25)Tr+yqk=Tr+Tuqk+Tr+Bupk−p−Tr+Aypk−p=Ir000uqk+Tr+Bupk−p−Tr+Aypk−p.It can be seen from equation (25) that only the first r elements of uqk are selected. However, as we mentioned in equation (14), uqk comprises drive force data from uk to uk+q−1, so the first r elements might not be meaningful in practice. Therefore, in order to select the drive force data of the length l from uk to uk+l−1, a new T¯ is formed by picking up the first l⋅ni row l⋅ni<r of Tr+. Substituting Tr+ in equation (25) with T¯ yields a new equation:(26)ulk=T¯yqk+B¯ypk−p−A¯upk−p,where yq,yp, and up are the same as in equation (14) and (15) and(27)ulk=ukuk+1⋮uk+l−1,(28)T¯=Il⋅ni0Tr+,B¯=T¯A,A¯=T¯B,in which T¯ is a l⋅ni×q⋅no matrix, A¯ is a l⋅ni×p⋅ni matrix, and B¯ is a l⋅ni×p⋅no matrix. As mentioned before, p is the number of time steps of past drives and responses and q is the number of time steps of future responses. Besides, the new parameter l stands for the time steps of future drives the inverse system generates.Needless to say, equation (26) bears remarkable resemblance to the multistep prediction model in equation (13), except that the T¯ matrix is truncated, which precludes the instability in the inverse system caused by the singularity problem while computing the inverse of T. Hence, the model defined in equation (26) is named as the inverse multistep prediction model hereafter and will be used in the MIMO swept-sine test.


### 2.2.2. Drive Signal Generation

In the real world, test articles start vibrating from rest conditions which cannot achieve stable response immediately; a segment of transient response will always exist. The reference signals for the swept-sine test actually specify the stable response of the system, so if the reference signals are used as the input for the inverse multistep prediction model directly, the drive signals generated will be unrealistic. Hence, a transition segment should be added to the forepart of the original reference signals in order to generate reasonable drive signals as illustrated in Figure1.Figure 1
Modification of the original reference signal.To compose this transition segment, a segment of fixed frequency sine signal which matches the initial amplitude and phase of the original reference signal is generated, and a 1/4 sine window is applied to ensure the amplitude of this segment varies smoothly from zero to the same initial amplitude of the original reference signal. This transition segment is inserted before the original reference signal to modify the reference signal. Obviously, the response of this transition segment is redundant and should be ignored.Since equation (26) produces ulk instead of uqk, using the reference signals with a length of q can only generate drive signals with a length of l, where l < q. Hence, the reference signals will be used in an overlapped manner in order to generate full length drive signals. In the drive signal generating process, if the reference signals are not enough in the last segment, zeros will be supplemented. The process of drive signal generation is illustrated in Figure 2 (for clarity, the figure shows only one channel, but this process is applicable for MIMO systems).Figure 2
Drive signal is generated by the inverse system method and the reference signal is used in an overlapped manner.From the above discussions, the generation of drive signals for the MIMO swept-sine test can be concluded as follows (also illustrated in Figure3):(a)
Estimate the finite difference model of the testing system(b)
ObtainT, A, and B in equation (13) from the finite difference model(c)
Calculate the pseudoinverse of T by the TSVD method to getTr+ and then to form T¯ from Tr+(d)
Establish the inverse multistep prediction model as the inverse system using equation (28)(e)
Modify the reference signals and use them as the inputs of the inverse system to produce drive signalsFigure 3
Generation of drive signals for the MIMO swept-sine test.


### 2.2.3. Control Strategy

Due to the noises in the input and output signals and other errors in the modeling of the system, the responses of the first loop may not be satisfactory without control. An offline control strategy is put forward, in which the response spectra are obtained and compared to the reference counterparts, and thus the errors in frequency domain are obtained. Then, the reference spectra for next control loop are updated according to these errors, and the drive signals are updated afterwards. To avoid damaging the test article and equipment, it advisable to control the swept-sine responses within acceptable range at a lower level and then restore the original level test.The reference spectra is updated as(29)Enω=R0iωSyiω,Rn+1ω=EniωεRniω,i=1,2,…,no,n=0,1,2,…where Syω is the response spectra at present loop, Rnω is the reference spectra after n times of modification (R0 is the actual reference spectra), Enω is the spectral error at present loop, Rn+1ω is the modified reference spectrum for the next loop, and ε is a positive constant and ε<1 [23]. All these spectra take the form of vectors, and they are used in equation (29) in an element-by-element manner, namely, Syiω,Eni, andRniωi=1,2,…,no represent the ith elemnt of the corresponding vectors.The control scheme for the swept-sine test is shown in Figure4.Figure 4
Offline swept-sine control scheme where A andΦ are amplitude and initial phase vectors defined in reference spectra.


## 2.2.1. Inverse Multistep Prediction Model

If the T matrix in equation (13) were always invertible, then the inverse model could be obtained in a quite straight forward way, which is unfortunately not always true for test systems. Therefore, it is necessary to use pseudoinverse derived in equation (23) instead. Left multiply equation (13) by Tr+ on both sides to yield(25)Tr+yqk=Tr+Tuqk+Tr+Bupk−p−Tr+Aypk−p=Ir000uqk+Tr+Bupk−p−Tr+Aypk−p.It can be seen from equation (25) that only the first r elements of uqk are selected. However, as we mentioned in equation (14), uqk comprises drive force data from uk to uk+q−1, so the first r elements might not be meaningful in practice. Therefore, in order to select the drive force data of the length l from uk to uk+l−1, a new T¯ is formed by picking up the first l⋅ni row l⋅ni<r of Tr+. Substituting Tr+ in equation (25) with T¯ yields a new equation:(26)ulk=T¯yqk+B¯ypk−p−A¯upk−p,where yq,yp, and up are the same as in equation (14) and (15) and(27)ulk=ukuk+1⋮uk+l−1,(28)T¯=Il⋅ni0Tr+,B¯=T¯A,A¯=T¯B,in which T¯ is a l⋅ni×q⋅no matrix, A¯ is a l⋅ni×p⋅ni matrix, and B¯ is a l⋅ni×p⋅no matrix. As mentioned before, p is the number of time steps of past drives and responses and q is the number of time steps of future responses. Besides, the new parameter l stands for the time steps of future drives the inverse system generates.Needless to say, equation (26) bears remarkable resemblance to the multistep prediction model in equation (13), except that the T¯ matrix is truncated, which precludes the instability in the inverse system caused by the singularity problem while computing the inverse of T. Hence, the model defined in equation (26) is named as the inverse multistep prediction model hereafter and will be used in the MIMO swept-sine test.


## 2.2.2. Drive Signal Generation

In the real world, test articles start vibrating from rest conditions which cannot achieve stable response immediately; a segment of transient response will always exist. The reference signals for the swept-sine test actually specify the stable response of the system, so if the reference signals are used as the input for the inverse multistep prediction model directly, the drive signals generated will be unrealistic. Hence, a transition segment should be added to the forepart of the original reference signals in order to generate reasonable drive signals as illustrated in Figure1.Figure 1
Modification of the original reference signal.To compose this transition segment, a segment of fixed frequency sine signal which matches the initial amplitude and phase of the original reference signal is generated, and a 1/4 sine window is applied to ensure the amplitude of this segment varies smoothly from zero to the same initial amplitude of the original reference signal. This transition segment is inserted before the original reference signal to modify the reference signal. Obviously, the response of this transition segment is redundant and should be ignored.Since equation (26) produces ulk instead of uqk, using the reference signals with a length of q can only generate drive signals with a length of l, where l < q. Hence, the reference signals will be used in an overlapped manner in order to generate full length drive signals. In the drive signal generating process, if the reference signals are not enough in the last segment, zeros will be supplemented. The process of drive signal generation is illustrated in Figure 2 (for clarity, the figure shows only one channel, but this process is applicable for MIMO systems).Figure 2
Drive signal is generated by the inverse system method and the reference signal is used in an overlapped manner.From the above discussions, the generation of drive signals for the MIMO swept-sine test can be concluded as follows (also illustrated in Figure3):(a)
Estimate the finite difference model of the testing system(b)
ObtainT, A, and B in equation (13) from the finite difference model(c)
Calculate the pseudoinverse of T by the TSVD method to getTr+ and then to form T¯ from Tr+(d)
Establish the inverse multistep prediction model as the inverse system using equation (28)(e)
Modify the reference signals and use them as the inputs of the inverse system to produce drive signalsFigure 3
Generation of drive signals for the MIMO swept-sine test.


## 2.2.3. Control Strategy

Due to the noises in the input and output signals and other errors in the modeling of the system, the responses of the first loop may not be satisfactory without control. An offline control strategy is put forward, in which the response spectra are obtained and compared to the reference counterparts, and thus the errors in frequency domain are obtained. Then, the reference spectra for next control loop are updated according to these errors, and the drive signals are updated afterwards. To avoid damaging the test article and equipment, it advisable to control the swept-sine responses within acceptable range at a lower level and then restore the original level test.The reference spectra is updated as(29)Enω=R0iωSyiω,Rn+1ω=EniωεRniω,i=1,2,…,no,n=0,1,2,…where Syω is the response spectra at present loop, Rnω is the reference spectra after n times of modification (R0 is the actual reference spectra), Enω is the spectral error at present loop, Rn+1ω is the modified reference spectrum for the next loop, and ε is a positive constant and ε<1 [23]. All these spectra take the form of vectors, and they are used in equation (29) in an element-by-element manner, namely, Syiω,Eni, andRniωi=1,2,…,no represent the ith elemnt of the corresponding vectors.The control scheme for the swept-sine test is shown in Figure4.Figure 4
Offline swept-sine control scheme where A andΦ are amplitude and initial phase vectors defined in reference spectra.


## 3. Simulation Example

### 3.1. Simulation Settings

A two-input-two-output swept-sine vibration test is simulated while using an aluminum cantilever beam as the test object. The parameters of the beam are listed in Table1 and shown in Figure 5.Table 1
Physical properties of the beam.
PropertyValueDensity (kg/m3)7850Elastic modulus (MPa)720Length (mm)1000Width (mm)50Height (mm)15Figure 5
Cantilever beam used in simulation.The finite element model (FEM) of the cantilever beam is used in the simulation. The beam is divided evenly into 10 elements and has 20 freedoms, and all modal damping ratios are set as 0.8%. Thus, the mass (M), stiffness (K), and damping (Ξ) matrices of the system can be set up from this finite element model. The continuous state-space model of the system is(30)x.=Acx+Bcu,y=Ccx+Dcu,where u and y are the input and output vectors, x is the state vector, and(31)Ac=0I−M−1K−M−1Ξ,Bc=0M−1Bf,Cc=−CaM−1K−CaM−1Ξ,Dc=CaM−1Bf,where Bf is the indicating matrix used to define the excitation force on each DOF and Ca is used to pick up output DOFs of the state-space model. Responses can be acquired from given forces with this continuous state-space model.In the example, the sampling frequency is 5120 Hz, the multistep prediction model of the system is estimated by a length of 1024 steps, the inverse multistep prediction model is truncated to a length of 512 steps, the sweep mode is logarithm, and the sweep speed is 2 octaves per minute (a sweep from 5 Hz to 250 Hz takes roughly 169 seconds). The reference spectra at two control points are defined in Tables2 and 3. The acceleration between the break points in the table is interpolated linearly on log-log axes, and the phases are formed by linear interpolation with frequency in a logarithmic scale.Table 2
Reference spectrum at control point 1.
Freq. (Hz)Acc. (g)Disp. (mm)Phase (°)15—0.7520225.731—2031501—2042500.5—20Table 3
Reference spectrum at control point 2.
Freq. (Hz)Acc. (g)Disp. (mm)Phase (°)15—0.7540225.731—4032501—40In practice, the vibration at low frequency is usually given in the form of the peak-peak displacement restrictions to avoid damaging test equipment. In addition, the crossover frequency (25.73 Hz in the example) is where both requirements on acceleration and displacement are met. However, since the vibration in the test is often measured by using an accelerometer, the displacement defined in the low frequency must be transformed into acceleration. The transformation formula is(32)A=ω2D2,where A represents the amplitude of acceleration and D is the peak-peak amplitude of displacement.


### 3.2. Simulation Results

The results of the two-input-two-output swept-sine vibration test are shown in Figure6. Here, the amplitudes and phases of the responses at the two control points are obtained by using a Co-Quad analyzer [9].Figure 6
Simulation results of the two-input-two-output swept-sine vibration test example.
(a)(b)From Figure6, it can be seen that the amplitudes and phases of the responses at the controlled points match the references very well.


### 3.3. Comparison with Frequency-Domain Method

The key step in the frequency-domain methods is to get the impedance matrix by inversing the frequency response function matrix H(f) and then use it to obtain the amplitude and initial phase information of the drive signals. The calculation process is as follows:(33)Xf=H−1fYf,where X and Y are complex vectors. Then, the relationship between frequency and time of sweeping process f (t) is substituted for X (f) to obtain the amplitude A and initial phase Φ.(34)At=Xft,Φt=∠Xft.The drive signals can be obtained by substituting the amplitude and the initial phase into equation (6). For the low damping model used in simulation, the condition number of the frequency response function matrix near the formant increases sharply, as shown in Figure 7, which makes the corresponding impedance matrix rather sensitive. Therefore, when calculating the impedance matrix, it is generally necessary to use regularization methods such as TSVD to avoid generating excessive excitation.Figure 7
The condition number of FRF.By comparing the inverse multistep prediction method proposed in this paper with the general frequency-domain method (Figure8), it can be seen that the method proposed in this paper can get better results when the frequency sweeps at a rather fast rate (2 oct/min) in a low damping system. This is because the frequency response function describes the relationship between sinusoidal excitation and response in the steady state, but the sweeping process is transient and deviates from the steady-state hypothesis as the sweeping rate becomes faster. On the other hand, this issue can be avoided by the proposed method since it directly generates the drive signals according to the time-domain reference signals.Figure 8
Comparison between the inverse multistep prediction method (IMSP) and frequency-domain method (FD).
(a)(b)


## 3.1. Simulation Settings

A two-input-two-output swept-sine vibration test is simulated while using an aluminum cantilever beam as the test object. The parameters of the beam are listed in Table1 and shown in Figure 5.Table 1
Physical properties of the beam.
PropertyValueDensity (kg/m3)7850Elastic modulus (MPa)720Length (mm)1000Width (mm)50Height (mm)15Figure 5
Cantilever beam used in simulation.The finite element model (FEM) of the cantilever beam is used in the simulation. The beam is divided evenly into 10 elements and has 20 freedoms, and all modal damping ratios are set as 0.8%. Thus, the mass (M), stiffness (K), and damping (Ξ) matrices of the system can be set up from this finite element model. The continuous state-space model of the system is(30)x.=Acx+Bcu,y=Ccx+Dcu,where u and y are the input and output vectors, x is the state vector, and(31)Ac=0I−M−1K−M−1Ξ,Bc=0M−1Bf,Cc=−CaM−1K−CaM−1Ξ,Dc=CaM−1Bf,where Bf is the indicating matrix used to define the excitation force on each DOF and Ca is used to pick up output DOFs of the state-space model. Responses can be acquired from given forces with this continuous state-space model.In the example, the sampling frequency is 5120 Hz, the multistep prediction model of the system is estimated by a length of 1024 steps, the inverse multistep prediction model is truncated to a length of 512 steps, the sweep mode is logarithm, and the sweep speed is 2 octaves per minute (a sweep from 5 Hz to 250 Hz takes roughly 169 seconds). The reference spectra at two control points are defined in Tables2 and 3. The acceleration between the break points in the table is interpolated linearly on log-log axes, and the phases are formed by linear interpolation with frequency in a logarithmic scale.Table 2
Reference spectrum at control point 1.
Freq. (Hz)Acc. (g)Disp. (mm)Phase (°)15—0.7520225.731—2031501—2042500.5—20Table 3
Reference spectrum at control point 2.
Freq. (Hz)Acc. (g)Disp. (mm)Phase (°)15—0.7540225.731—4032501—40In practice, the vibration at low frequency is usually given in the form of the peak-peak displacement restrictions to avoid damaging test equipment. In addition, the crossover frequency (25.73 Hz in the example) is where both requirements on acceleration and displacement are met. However, since the vibration in the test is often measured by using an accelerometer, the displacement defined in the low frequency must be transformed into acceleration. The transformation formula is(32)A=ω2D2,where A represents the amplitude of acceleration and D is the peak-peak amplitude of displacement.


## 3.2. Simulation Results

The results of the two-input-two-output swept-sine vibration test are shown in Figure6. Here, the amplitudes and phases of the responses at the two control points are obtained by using a Co-Quad analyzer [9].Figure 6
Simulation results of the two-input-two-output swept-sine vibration test example.
(a)(b)From Figure6, it can be seen that the amplitudes and phases of the responses at the controlled points match the references very well.


## 3.3. Comparison with Frequency-Domain Method

The key step in the frequency-domain methods is to get the impedance matrix by inversing the frequency response function matrix H(f) and then use it to obtain the amplitude and initial phase information of the drive signals. The calculation process is as follows:(33)Xf=H−1fYf,where X and Y are complex vectors. Then, the relationship between frequency and time of sweeping process f (t) is substituted for X (f) to obtain the amplitude A and initial phase Φ.(34)At=Xft,Φt=∠Xft.The drive signals can be obtained by substituting the amplitude and the initial phase into equation (6). For the low damping model used in simulation, the condition number of the frequency response function matrix near the formant increases sharply, as shown in Figure 7, which makes the corresponding impedance matrix rather sensitive. Therefore, when calculating the impedance matrix, it is generally necessary to use regularization methods such as TSVD to avoid generating excessive excitation.Figure 7
The condition number of FRF.By comparing the inverse multistep prediction method proposed in this paper with the general frequency-domain method (Figure8), it can be seen that the method proposed in this paper can get better results when the frequency sweeps at a rather fast rate (2 oct/min) in a low damping system. This is because the frequency response function describes the relationship between sinusoidal excitation and response in the steady state, but the sweeping process is transient and deviates from the steady-state hypothesis as the sweeping rate becomes faster. On the other hand, this issue can be avoided by the proposed method since it directly generates the drive signals according to the time-domain reference signals.Figure 8
Comparison between the inverse multistep prediction method (IMSP) and frequency-domain method (FD).
(a)(b)


## 4. Test Analysis

### 4.1. Test Settings

A two-input-two-output swept-sine test is conducted with a three-axis shaker to validate the method proposed in the paper. The test system exhibited in Figure9 consists of a three-axis vibration shaker (Shinken G-6080-3HT-20), a VXI digital signal processor, two acceleration sensors, a computer, and other peripherals.Figure 9
Test system constitution.Due to the limited number of the VXI channels, only the accelerations ofx and y directions of the 3-axis shaker are controlled. The VXI digital signal processor and the computer are used to form the vibration controller. In Figure 9, the first two outputs (OUT 1∼2) are used as drive signal channels and the other two outputs (OUT 3∼4) are used to transfer two swept-sine signals which have same amplitudes but with the phase difference of 90°. These two swept-sine signals are used as “standard signals” in order to obtain the transient frequency and to extract amplitudes and phases from responses. The first four input channels (IN 1∼4) are the direct feedback from the output channels, and the last two input channels (IN 5∼6) are used to acquire responses from the acceleration sensors. The roles of these inputs, outputs, and measurements are also shown in Figure 10.Figure 10
Block diagram for the test system.The reference spectra given in the simulation example in Section3 are used again in the test, but additional abort boundaries are set at ±6 dB. During the test, Hilbert transform [24] is applied to one of the “standard” swept-sine signals to obtain the transient frequency and amplitude. The amplitudes and phases of the responses are obtained by using a Co-Quad analyzer [9].


### 4.2. Test Results

A MIMO vibration test generally starts from a level lower than the normal references and climbs up to the normal reference level after it has been controlled. This is due to the fact that in MIMO cases the test cannot be controlled instantly, so in order to avoid any possible damage to the equipment, a preliminary test at a lower level is necessary. In this test, the preliminary test is conducted at the –6 dB level of the references. After 4 iteration loops, the response spectra are controlled within acceptable tolerances, and then the normal level test can be performed. The response spectra before and after control in the lower level test are shown in Figures11 and 12, respectively. The responses in the normal level test are given in Figure 13. From Figure 11, it can be seen that the response amplitudes are roughly within the limits except for the low frequencies, but the phases have some obtrusive deviations without control. In Figure 12, it is clear that after control both the amplitudes and phases are within the acceptable limits. From Figure 13, it can be seen that after the normal level is applied, the response spectra become slightly worse, but they are still within the alert limits.Figure 11
Responses of the lower level (−6 dB) preliminary test without control.
(a)(b)Figure 12
Responses of the lower level (−6 dB) preliminary test after control (4 iterations).
(a)(b)Figure 13
Responses of the normal level test after control.
(a)(b)


## 4.1. Test Settings

A two-input-two-output swept-sine test is conducted with a three-axis shaker to validate the method proposed in the paper. The test system exhibited in Figure9 consists of a three-axis vibration shaker (Shinken G-6080-3HT-20), a VXI digital signal processor, two acceleration sensors, a computer, and other peripherals.Figure 9
Test system constitution.Due to the limited number of the VXI channels, only the accelerations ofx and y directions of the 3-axis shaker are controlled. The VXI digital signal processor and the computer are used to form the vibration controller. In Figure 9, the first two outputs (OUT 1∼2) are used as drive signal channels and the other two outputs (OUT 3∼4) are used to transfer two swept-sine signals which have same amplitudes but with the phase difference of 90°. These two swept-sine signals are used as “standard signals” in order to obtain the transient frequency and to extract amplitudes and phases from responses. The first four input channels (IN 1∼4) are the direct feedback from the output channels, and the last two input channels (IN 5∼6) are used to acquire responses from the acceleration sensors. The roles of these inputs, outputs, and measurements are also shown in Figure 10.Figure 10
Block diagram for the test system.The reference spectra given in the simulation example in Section3 are used again in the test, but additional abort boundaries are set at ±6 dB. During the test, Hilbert transform [24] is applied to one of the “standard” swept-sine signals to obtain the transient frequency and amplitude. The amplitudes and phases of the responses are obtained by using a Co-Quad analyzer [9].


## 4.2. Test Results

A MIMO vibration test generally starts from a level lower than the normal references and climbs up to the normal reference level after it has been controlled. This is due to the fact that in MIMO cases the test cannot be controlled instantly, so in order to avoid any possible damage to the equipment, a preliminary test at a lower level is necessary. In this test, the preliminary test is conducted at the –6 dB level of the references. After 4 iteration loops, the response spectra are controlled within acceptable tolerances, and then the normal level test can be performed. The response spectra before and after control in the lower level test are shown in Figures11 and 12, respectively. The responses in the normal level test are given in Figure 13. From Figure 11, it can be seen that the response amplitudes are roughly within the limits except for the low frequencies, but the phases have some obtrusive deviations without control. In Figure 12, it is clear that after control both the amplitudes and phases are within the acceptable limits. From Figure 13, it can be seen that after the normal level is applied, the response spectra become slightly worse, but they are still within the alert limits.Figure 11
Responses of the lower level (−6 dB) preliminary test without control.
(a)(b)Figure 12
Responses of the lower level (−6 dB) preliminary test after control (4 iterations).
(a)(b)Figure 13
Responses of the normal level test after control.
(a)(b)


## 5. Conclusion

An inverse multistep prediction model to represent the inverse system for the generation of drive signals for the MIMO continuous swept-sine test is proposed in the paper. TheT¯ matrix in the inverse multistep prediction model is truncated, thus causing the model to be noncausal and able to generate more stable drive signals. Since the length of the generated drive signals is less than the length of the input reference signals, the reference signals are used in an overlapped manner. In addition, in order to prevent the sudden occurrence of large drive damage to the test pieces and other test equipment, a transition section is added to the starting part of the reference signal; an offline control method is proposed in which the reference signals are modified by multiplication, and the drive signals are then updated using the inverse multistep prediction model.A numerical example of a MIMO continuous swept-sine test on a cantilever beam is put forward, and the results show that the proposed method can perform excellently if the inverse multistep prediction model can be obtained precisely. The method is also verified by a MIMO continuous swept-sine test on a 3-axis shaker. The test results show that the response spectra can be controlled within the acceptable limits.



---

*Source: 8869802-2020-09-28.xml*