寧西占 黃偉

摘? ?要：為解決實(shí)時(shí)混合試驗(yàn)系統(tǒng)中時(shí)變時(shí)滯嚴(yán)重影響試驗(yàn)穩(wěn)定性和精度的問題，以基于混合靈敏度的H∞魯棒控制方法設(shè)計(jì)反饋控制器，使試驗(yàn)系統(tǒng)具有穩(wěn)定的動(dòng)態(tài)特性;以多項(xiàng)式外插進(jìn)一步消除時(shí)滯，并將其應(yīng)用于線性Benchmark問題.為充分探討所提方法性能，采用線性和非線性加載系統(tǒng)模型，開展了物理子結(jié)構(gòu)為線性和非線性的虛擬實(shí)時(shí)混合試驗(yàn)研究.結(jié)果顯示：作動(dòng)器的實(shí)測(cè)位移均與期望位移幾乎完全重合;與已有方法的對(duì)比顯示，采用魯棒時(shí)滯補(bǔ)償?shù)膶?shí)時(shí)混合試驗(yàn)具有更小的誤差;而對(duì)非線性物理子結(jié)構(gòu)，試驗(yàn)系統(tǒng)的均方根誤差和峰值誤差均在0.6%～1.5%之間，且均方根誤差大于10%和峰值誤差大于6%的概率遠(yuǎn)小于0.05.研究結(jié)果表明，魯棒時(shí)滯補(bǔ)償方法可大幅提高實(shí)時(shí)混合試驗(yàn)的模擬精度，并具有較強(qiáng)的魯棒性.

關(guān)鍵詞：實(shí)時(shí)混合試驗(yàn);時(shí)滯;H∞控制;魯棒性;不確定性

中圖分類號(hào)：TU317? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

Numerical Study of Robust Time-delay Compensation

Method for Real-time Hybrid Simulation

NING Xizhan1，2，HUANG Wei1

（1. College of Civil Engineering，Huaqiao University，Xiamen 361021，China;

2. Key Laboratory for Intelligent Infrastructure and Monitoring of Fujian Province（Huaqiao University），Xiamen 361021，China）

Abstract：To address the problem that the variable time delay seriously affects the stability and accuracy of the real-time hybrid simulation（RTHS），a robust time-delay compensation method was studied. A feedback controller? based on H∞ robust control method of mixed sensitivity was designed to stabilize the dynamics of the testing system，while the polynomial extrapolation method was employed to compensate for the time-delay. The proposed method was applied to a linear Benchmark problem. To further investigate the performance of the proposed method，virtual RTHSs were carried out with linear and nonlinear physical substructure（PS） employing linear and nonlinear models of the loading system. The results showed the measured displacements were almost identical to the desired displacements. Compared with the existing method，RTHS with the robust time-delay compensation method had a smaller error. For nonlinear PS，the root-mean-square error（RMSE） and peak error（PE） were both between 0.6% and 1.5%，and the probability of RMSE greater than 10% and PE greater than 6% was far less than 0.05. The study demonstrated that the robust delay compensation method can improve the accuracy of RTHS greatly and behavior strong robustness.

Key words：real-time hybrid simulation;time delay;H∞ control;robustness;uncertainty

混合試驗(yàn)源于1969年日本學(xué)者[1]提出的擬動(dòng)力試驗(yàn)方法，它將數(shù)值計(jì)算和物理試驗(yàn)通過在線的方式有機(jī)結(jié)合在一起，是評(píng)估結(jié)構(gòu)在動(dòng)荷載作用下性能的經(jīng)濟(jì)、實(shí)用的方法. 通常，混合試驗(yàn)方法將待評(píng)估結(jié)構(gòu)中具有較強(qiáng)非線性的部分取出在實(shí)驗(yàn)室進(jìn)行物理加載，剩余部分在計(jì)算機(jī)中模擬. 前者稱為物理子結(jié)構(gòu)，后者稱為數(shù)值子結(jié)構(gòu)，子結(jié)構(gòu)間通過傳遞系統(tǒng)（通常是作動(dòng)器）聯(lián)系在一起. 混合試驗(yàn)提出之初采用了慢速加載方式，并在逐步積分方法、誤差傳遞等方面取得了一系列研究進(jìn)展[2-4].

Nakashima等在進(jìn)行混合試驗(yàn)時(shí)對(duì)試件采用動(dòng)力加載的方式，提出了實(shí)時(shí)混合試驗(yàn)，解決了速率相關(guān)型試件的混合試驗(yàn)問題[5]. 在實(shí)時(shí)混合試驗(yàn)中，數(shù)值計(jì)算、邊界條件實(shí)現(xiàn)及數(shù)據(jù)采集等均需在十分短的時(shí)間步長(zhǎng)內(nèi)完成，以實(shí)現(xiàn)子結(jié)構(gòu)間邊界位移協(xié)調(diào)和力的平衡. 受數(shù)據(jù)交互、模擬信號(hào)與數(shù)字信號(hào)間相互轉(zhuǎn)換、作動(dòng)器的動(dòng)態(tài)性能以及作動(dòng)器與試件間相互作用的影響，作動(dòng)器接收的命令和實(shí)際響應(yīng)之間往往不同步，這種現(xiàn)象通常稱為時(shí)滯.時(shí)滯的存在將減弱實(shí)時(shí)混合試驗(yàn)的精度，甚至嚴(yán)重影響其穩(wěn)定性[6]. 因此，一系列時(shí)滯補(bǔ)償措施不斷被用于抑制時(shí)滯的影響，如基于常時(shí)滯假設(shè)的預(yù)測(cè)方法[6]，基于經(jīng)典控制理論的前饋、反饋及逆控制[7-9]，基于現(xiàn)代控制理論的滑動(dòng)模態(tài)控制[10]、H∞控制[11]，以及基于反饋力修正的方法[12]等.近年來，自適應(yīng)時(shí)滯補(bǔ)償策略[13-16]也得到了廣泛的關(guān)注. 然而，基于控制理論的時(shí)滯補(bǔ)償方法大多數(shù)具有模型依賴性，而反饋控制無法完全消除時(shí)滯的影響. 此外，考慮到物理加載的動(dòng)態(tài)特性、試件的非線性、噪聲以及電壓變化導(dǎo)致的供油壓力變化等因素，前述時(shí)滯補(bǔ)償方法效果被削弱.

針對(duì)以上問題，本文提出了一種魯棒時(shí)滯補(bǔ)償方法，并通過Benchmark模型進(jìn)行了初步驗(yàn)證[15]. 為深入分析魯棒時(shí)滯補(bǔ)償方法的性能，本文首先基于作動(dòng)器的Benchmark模型，采用線性物理子結(jié)構(gòu)，對(duì)比研究了魯棒時(shí)滯補(bǔ)償方法和反饋修正的自適應(yīng)時(shí)滯補(bǔ)償方法;之后，討論了魯棒時(shí)滯補(bǔ)償方法處理非線性物理子結(jié)構(gòu)的能力;最后，基于考慮流量非線性的作動(dòng)器模型，以磁流變阻尼器為非線性物理子結(jié)構(gòu)，研究了考慮試驗(yàn)系統(tǒng)（作動(dòng)器-試件系統(tǒng)）非線性時(shí)魯棒時(shí)滯補(bǔ)償方法的性能.

1? ?魯棒時(shí)滯補(bǔ)償方法

在已有的時(shí)滯補(bǔ)償方法中，大多數(shù)具有模型依賴性，而數(shù)學(xué)模型不能考慮物理子結(jié)構(gòu)的制造誤差、觀測(cè)噪聲、外界擾動(dòng)以及加載系統(tǒng)動(dòng)態(tài)性能改變等因素.而在控制領(lǐng)域中，上述因素通常歸為未建模動(dòng)態(tài)和模型參數(shù)不確定性. 因此，筆者提出了魯棒時(shí)滯補(bǔ)償方法，并將其應(yīng)用于線性Benchmark問題[15].該方法的原理如圖1所示，其中dD是期望位移，dm是測(cè)量位移，dc是命令位移，d是外插預(yù)測(cè)位移.該方法中，時(shí)滯補(bǔ)償器由H∞控制器和多項(xiàng)式外插兩部分組成，前者用于改善被控對(duì)象的動(dòng)態(tài)性能，得到一個(gè)穩(wěn)定且具有魯棒性的閉環(huán)系統(tǒng)，后者則用于進(jìn)一步消除系統(tǒng)的時(shí)滯.

1.1? ?基于混合靈敏度的H∞控制器

H∞控制是通過建立從廣義輸入到性能輸出的性能指標(biāo)傳遞矩陣F，并在Hardy空間中以該性能指標(biāo)的無窮范數(shù)為衡量標(biāo)準(zhǔn)進(jìn)行優(yōu)化，從而得到控制器的一種方法.標(biāo)準(zhǔn)H∞控制問題可通過對(duì)一般反饋控制問題進(jìn)行性能加權(quán)獲取，本研究中加權(quán)的反饋控制框圖如圖2所示.這里，廣義輸入為位移d，性能輸出為z = [z1;z2;z3]，性能指標(biāo)為：

|| F ||∞ = W1 SW2 KSW3 T ∞? ? ? ? ? ?（1）

式中：S = 1/（1 + PK）為靈敏度函數(shù);T = PK/（1 + PK）為補(bǔ)靈敏度函數(shù);P是控制對(duì)象（包含作動(dòng)器和物理子結(jié)構(gòu)）;W1、W2和W3是性能加權(quán)函數(shù)，直接影響系統(tǒng)的追蹤性能和魯棒性. 圖2所示控制問題的廣義被控對(duì)象及狀態(tài)空間實(shí)現(xiàn)為：

設(shè)計(jì)H∞控制器，即尋找一個(gè)控制器K，使得閉環(huán)系統(tǒng)是內(nèi)穩(wěn)定的，并且性能指標(biāo)滿足

|| F ||∞ = || G11 + G12 K（I - G22 K）-1G21||∞ < γ? ? （3）

式中：γ是一個(gè)給定的正數(shù)，下標(biāo)“∞”表示無窮范數(shù).設(shè)計(jì)H∞控制器有多種方法，如線性矩陣不等式法、DGKF法等[17]. 本研究采用DGKF法，即當(dāng)系統(tǒng)（A，B2）為可鎮(zhèn)定的，（C2，A）為可檢測(cè)時(shí)，可通過求解兩個(gè)黎卡提方程得到H∞控制器.

采用基于混合靈敏度的H∞控制設(shè)計(jì)時(shí)，控制器的追蹤性能和魯棒性直接受到三個(gè)權(quán)函數(shù)的影響.吳旭東和解學(xué)書[18]給出了權(quán)函數(shù)選取的一般性原則，這里不再贅述.需要強(qiáng)調(diào)的是，在進(jìn)行控制器設(shè)計(jì)時(shí)，權(quán)函數(shù)W1應(yīng)具有近似積分特性的形式，以實(shí)現(xiàn)在所關(guān)心的頻率范圍內(nèi)系統(tǒng)具有較強(qiáng)的追蹤性能和抗干擾能力;權(quán)函數(shù)W3在低頻段遠(yuǎn)小于1并且在高頻段應(yīng)具有近似s2特性，以實(shí)現(xiàn)對(duì)命令信號(hào)的追蹤和對(duì)高頻噪聲的抑制. 因此，本研究中權(quán)函數(shù)W1和W3具有如下形式：

式中：a1、a2、b1、b2和b3是常數(shù)，且a2是一個(gè)很小的正數(shù)，以使控制器具有積分特性.

1.2? ?多項(xiàng)式外插

H∞控制本質(zhì)上仍是反饋控制，控制后系統(tǒng)必然仍有時(shí)滯存在.但采用H∞控制器對(duì)被控對(duì)象的動(dòng)態(tài)性能進(jìn)行修正后，系統(tǒng)將具有穩(wěn)定的動(dòng)態(tài)性能，具備較強(qiáng)的追蹤性能和抗干擾能力，且系統(tǒng)在較寬頻率范圍內(nèi)將具有穩(wěn)定的時(shí)滯[11，15]. 通常，采用H∞控制修正后的新系統(tǒng)可看作常時(shí)滯系統(tǒng)，可采用經(jīng)典的多項(xiàng)式外插法將該時(shí)滯進(jìn)一步消除.

本研究中，多項(xiàng)式外插采用性能較好的三階外插方法[19]，其表達(dá)式為：

d（ti + τ） = q1 dDi - q2 dDi-1 + q3 dDi-2 - q4 dDi-3? ? （5）

式中：q1 = 1 + 11η/6 + η2 + η3/6;q2 = 3η + 3η2/2 + η3/2;q3 = 3η/2 + 2η2 + η3/2;q4 = η/3 + η2/2 + η3/6;η = τ/Δt，τ是新系統(tǒng)的時(shí)滯，Δt是積分步長(zhǎng).

2? ?作動(dòng)器數(shù)值模型

2.1? ?線性模型

為推進(jìn)實(shí)時(shí)混合試驗(yàn)的發(fā)展，由美國自然科學(xué)基金資助的混合試驗(yàn)多災(zāi)害工程研究體發(fā)布了實(shí)時(shí)混合試驗(yàn)的Benchmark問題[20]，用于評(píng)估不同時(shí)滯補(bǔ)償方法的性能，尤其是魯棒性.該Benchmark問題建立在實(shí)際試驗(yàn)系統(tǒng)之上，待模擬結(jié)構(gòu)是一個(gè)三層兩跨的鋼框架結(jié)構(gòu)，剛度矩陣為KS=[2.605 5，-2.313 4，0.593 7;-2.313 4，3.256 1，-1.442 0;0.593 7，-1.442 0， 0.926 7] × 107 N/m，采用集中質(zhì)量矩陣，每層質(zhì)量相同，并以底層左跨為物理子結(jié)構(gòu);同時(shí)，Benchmark提供了電液伺服作動(dòng)器的線性模型，將液壓伺服閥、液壓作動(dòng)筒和力/位移傳感器等元件以集中參數(shù)表示，物理子結(jié)構(gòu)是線性試件，如圖3所示. 圖中α1、 β0、 β1和β2與液壓伺服閥相關(guān)，α2與自然速率反饋相關(guān)，α3與液壓作動(dòng)筒相關(guān)，mp、cp和kp是物理子結(jié)構(gòu)的質(zhì)量、阻尼和剛度系數(shù). 為考慮液壓傳遞系統(tǒng)、傳感器與控制實(shí)現(xiàn)硬件和試驗(yàn)子系統(tǒng)等的不確定性，參數(shù)β1、β2、α3和kp采用服從標(biāo)準(zhǔn)正態(tài)分布的隨機(jī)變量來描述. 上述參數(shù)取值及其分布情況見表1.

2.2? ?考慮流量非線性的模型

作動(dòng)器線性模型無法考慮加載系統(tǒng)使用過程中非線性引起的時(shí)滯變化，而流量引起的非線性在電液伺服系統(tǒng)中占主要成分. Zhao等[21]提出的作動(dòng)器模型可以有效考慮上述非線性，如圖4所示. 圖中，xv是伺服閥位移，Q是流量，p是負(fù)載油壓，Ka是油的壓縮系數(shù)，Cl是油的泄露系數(shù)，A是活塞面積. 伺服閥的傳遞函數(shù)為：

式中：TPID是模擬控制器傳遞函數(shù)，通常只取比例增益;Kvp是閥的增益系數(shù);xv max是最大閥芯位移;Av是滑閥閥芯面積;K3是內(nèi)置位移傳感器靈敏度系數(shù);τ是滑閥時(shí)間常量. 流量非線性關(guān)系及作動(dòng)器相關(guān)參數(shù)分別見表2和表3.

3? ?模擬結(jié)果與分析

3.1? ?作動(dòng)器線性模型工況

為考察子結(jié)構(gòu)劃分對(duì)實(shí)時(shí)混合試驗(yàn)穩(wěn)定性和精度的影響，Benchmark問題給出了4種子結(jié)構(gòu)劃分方式，樓層質(zhì)量分別為1 000 kg、1 100 kg、1 300 kg和1 000 kg，相應(yīng)于每一種樓層質(zhì)量，對(duì)應(yīng)的阻尼比分別為5%、4%、3%和3%.對(duì)每一種子結(jié)構(gòu)劃分，考慮了El Centro波、Kobe波和Morgan Hill波等3種地震動(dòng)記錄，地震波的調(diào)幅均為0.4g，并分別完成了21次虛擬實(shí)時(shí)混合試驗(yàn)，共計(jì)252次虛擬實(shí)時(shí)混合試驗(yàn)，以驗(yàn)證魯棒時(shí)滯補(bǔ)償方法的有效性和魯棒性.同時(shí)，本文還研究了物理子結(jié)構(gòu)進(jìn)入非線性，并考慮恢復(fù)力模型為Bouc-Wen模型時(shí)補(bǔ)償器的性能.需要說明的是，本文采用了Kalman濾波器對(duì)含噪聲的位移信號(hào)進(jìn)行濾波.

3.1.1? ?補(bǔ)償器設(shè)計(jì)

以圖3所示系統(tǒng)為控制對(duì)象，完成了魯棒時(shí)滯補(bǔ)償器的設(shè)計(jì). 經(jīng)試算，設(shè)計(jì)H∞控制器的3個(gè)權(quán)函數(shù)為：

W3 = 1 × 10-5 s2 + 4 × 10-3 s + 5 × 10-2? ? ? ?（7）

最終，經(jīng)過29次迭代得到可行解，其中最優(yōu)γ值為1.522，相應(yīng)的控制器為：

式中：n5 = 3.94 × 106，n4 = 1.703 × 109，n3 = 1.139 ×1012，n2 = 3.151 × 1014，n1 = 7.355 × 1016，n0 = 2.933 ×1018，d5 = 2.615 × 105，d4 = 2.935 × 108，d3 = 2.136 ×1011，d2 = 9.397 × 1013，d1 = 2.427 × 1016，d0 = 2.369 ×1014.

經(jīng)H∞控制器修正后系統(tǒng)的動(dòng)態(tài)性能如圖5所示，為方便比較，無控系統(tǒng)的性能也在圖中給出.從圖5中可以看出，在關(guān)心的頻率范圍內(nèi)，被控后系統(tǒng)具有穩(wěn)定的幅值，且該幅值接近于0 dB，極大地改善了系統(tǒng)的追蹤性能;從時(shí)滯圖中可以看出，無控系統(tǒng)的時(shí)滯較大且變化明顯，而被控后系統(tǒng)在較寬的頻率范圍內(nèi)具有穩(wěn)定且較小的時(shí)滯.以上分析表明，被控后系統(tǒng)具有穩(wěn)定的動(dòng)態(tài)性能，且可被看作常時(shí)滯系統(tǒng).被控后系統(tǒng)時(shí)滯大小為0.008 4 s，因此，在多項(xiàng)式外插中η = 34.406 4.

3.1.2? ?線性物理子結(jié)構(gòu)

圖6給出了在El Centro地震激勵(lì)下樓層質(zhì)量和阻尼比分別為1 000 kg和3%的實(shí)時(shí)混合試驗(yàn)中作動(dòng)器的追蹤性能，其中參考位移是由數(shù)值子結(jié)構(gòu)計(jì)算得到的位移.從圖6（a）中可以看出，實(shí)測(cè)位移與參考位移吻合良好;從圖6（b）可以看出，實(shí)測(cè)位移與參考位移幾乎完全重合，表明本文所提時(shí)滯補(bǔ)償方法具有良好的追蹤性能.

本文采用魯棒時(shí)滯補(bǔ)償方法的虛擬實(shí)時(shí)混合試驗(yàn)結(jié)果與文獻(xiàn)[22]的結(jié)果進(jìn)行了對(duì)比，如圖7所示.圖中結(jié)果來自在El Centro地震激勵(lì)下虛擬實(shí)時(shí)混合試驗(yàn)的統(tǒng)計(jì)值.其中，誤差指標(biāo)J2和J3反映了作動(dòng)器的追蹤性能，分別是均方根誤差和峰值追蹤誤差;誤差指標(biāo)J4、J5和J6分別是第一、第二和第三層位移的均方根誤差;J7、J8和J9分別是第一、第二和第三層的峰值誤差. J2 ～ J9的定義參見文獻(xiàn)[20]. 從圖7（a）中可以看出，在8個(gè)誤差指標(biāo)的統(tǒng)計(jì)均值中，本文方法的最大值均小于5%，而文獻(xiàn)[22]的方法則均大于5%，且J4的最大值高達(dá)30%，表明本文方法具有較好的追蹤性能，可以大幅提高模擬精度.從圖7（b）可以看出，除誤差指標(biāo)J6略大于文獻(xiàn)[22]的方法外，本文方法的方差整體上是最小的，表明本文方法具有較強(qiáng)的魯棒性.

3.1.3? ?非線性物理子結(jié)構(gòu)

在虛擬實(shí)時(shí)混合試驗(yàn)中，記錄了非線性物理子結(jié)構(gòu)的實(shí)測(cè)位移和實(shí)測(cè)反力，得到其滯回曲線，如圖8所示. 從圖中可看出，物理子結(jié)構(gòu)進(jìn)入了較強(qiáng)的非線性.同時(shí)，本文對(duì)252次虛擬實(shí)時(shí)混合試驗(yàn)所得的均方根誤差和峰值追蹤誤差進(jìn)行了統(tǒng)計(jì)分析，得到了誤差的累積分布函數(shù)，如圖9所示. 從圖中可看出，均方根誤差的范圍為0.6%～1.3%，而峰值追蹤誤差則為0.7%～1.5%;且當(dāng)概率小于0.8時(shí)，均方根誤差和峰值追蹤誤差的值分別小于1%和1.25%.值得說明的是，采用魯棒時(shí)滯補(bǔ)償方法后，系統(tǒng)的計(jì)算時(shí)滯為0 ms.由于試件非線性的影響可歸為模型的參數(shù)不確定性和未建模動(dòng)態(tài)特性，而H∞控制可以很好地處理這些不確定性，因此本文方法對(duì)非線性試件依然具有較強(qiáng)的追蹤能力.

本節(jié)采用結(jié)構(gòu)第三層位移響應(yīng)來討論實(shí)時(shí)混合試驗(yàn)的模擬精度.由樣本統(tǒng)計(jì)和核密度估計(jì)[23]得到的均方根誤差和峰值誤差的統(tǒng)計(jì)直方圖與概率密度函數(shù)曲線如圖10所示，其中，在對(duì)概率密度進(jìn)行估計(jì)時(shí)，核密度估計(jì)方法的核函數(shù)為正態(tài)分布函數(shù)，窗寬由樣本的數(shù)目確定.從圖中可以看出，采用核密度估計(jì)得到的兩種誤差指標(biāo)的概率密度函數(shù)與由樣本統(tǒng)計(jì)得到的直方圖具有較好的一致性，且均方根誤差和峰值誤差的分布范圍較小，分別集中在0%～7%和0%～5%.同時(shí)，從圖中可以看出，均方根誤差大于10%和峰值誤差大于6%的概率遠(yuǎn)小于0.05.結(jié)果表明，采用魯棒時(shí)滯補(bǔ)償方法的實(shí)時(shí)混合試驗(yàn)，模擬結(jié)果的精度得到了很大改善，且表現(xiàn)出較強(qiáng)的魯棒性.

3.2? ?作動(dòng)器非線性模型工況

在本工況中，參考結(jié)構(gòu)是一個(gè)三層單跨的鋼框架結(jié)構(gòu)，并在底層安裝有一磁流變阻尼器.參考結(jié)構(gòu)每層的質(zhì)量為20 250 kg，剛度矩陣為KS = [9.933 3，-5.662，0;-5.662，1.134，-5.662;0，-5.662 ，5.662]×106 N/m，阻尼矩陣C = [724 2，-206 9，0;-206 9，

4 139，-206 9;0，-2 069，2 069] Ns/m.物理子結(jié)構(gòu)為磁流變阻尼器，采用Bouc-Wen模型描述，即：

f = αb z + cb x? ? ? ? （9）

式中：αb = 2 888.78 N;cb = 6 400 N/A2;γz = 200/m;βz = 200/m;n = 2，p = 5 000/m.

3.2.1? ?補(bǔ)償器設(shè)計(jì)

在進(jìn)行虛擬混合之前，對(duì)試驗(yàn)系統(tǒng)輸入幅值為0.01 m、頻率為0～40 Hz的正弦掃頻信號(hào)，識(shí)別得到試驗(yàn)系統(tǒng)的線性數(shù)值模型為：

之后，以式（11）為控制對(duì)象，進(jìn)行補(bǔ)償器設(shè)計(jì).經(jīng)試算，設(shè)計(jì)H∞控制器的3個(gè)權(quán)函數(shù)為：

經(jīng)過26次迭代得到可行解，其中最優(yōu)γ值為0.783，相應(yīng)的控制器為：

此時(shí)，新系統(tǒng)的時(shí)滯大小為0.007 8 s. 因此，在多項(xiàng)式外插中，η = 7.8.

3.2.2? ?模擬結(jié)果分析

在El Centro地震激勵(lì)下，作動(dòng)器的追蹤性能如圖11所示. 從圖11（a）中可以看出，作動(dòng)器的測(cè)量位移與期望位移吻合良好;從11（b）的局部放大圖可以看出，二者完全重合，表明作動(dòng)器具有較強(qiáng)的追蹤性能;從圖11（c）給出了期望位移與測(cè)量位移的關(guān)系圖，其中直線斜率為45°，進(jìn)一步表明，魯棒時(shí)滯補(bǔ)償方法在處理作動(dòng)器非線性和物理子結(jié)構(gòu)非線性時(shí)依然具有較強(qiáng)的追蹤性能.

4? ?結(jié)? ?論

為改善實(shí)時(shí)混合試驗(yàn)的穩(wěn)定性，提高模擬精度，提出了魯棒時(shí)滯補(bǔ)償方法，并將其應(yīng)用于解決線性問題的時(shí)滯補(bǔ)償.基于前期工作，本文針對(duì)非線性問題開展了深入細(xì)致的研究，主要結(jié)論如下：

1）經(jīng)H∞控制后試驗(yàn)系統(tǒng)的時(shí)滯明顯減小，動(dòng)態(tài)性能得到了極大改善;在關(guān)心的頻率范圍內(nèi)，新系統(tǒng)可近似看作常時(shí)滯系統(tǒng)，且該時(shí)滯可采用多項(xiàng)式外插進(jìn)一步消除.

2）考慮液壓傳遞系統(tǒng)、傳感器與控制硬件以及物理子結(jié)構(gòu)等的不確定性和非線性時(shí)，經(jīng)本文方法修正后的試驗(yàn)系統(tǒng)依然可幾乎完全實(shí)現(xiàn)期望位移，使實(shí)時(shí)混合試驗(yàn)的模擬精度大幅提高，且表現(xiàn)出較強(qiáng)的魯棒性.

3）雖然魯棒時(shí)滯補(bǔ)償方法具有較強(qiáng)的追蹤性能和魯棒性，但本文只完成了數(shù)值研究，該方法在真實(shí)試驗(yàn)中的性能有待進(jìn)一步驗(yàn)證.

參考文獻(xiàn)

[1]? ? HAKUNO M，SHIDAWARA M，HARA T. Dynamic destructive test of a cantilever beam，controlled by an analog-computer[J]. Proceedings of the Japan Society of Civil Engineers，1969，1969（171）：1—9.

[2]? ? TAKANASHI K，NAKASHIMA M. Japanese activities on on-line testing[J]. Journal of Engineering Mechanics，1987，113（7）：1014—1032.

[3]? ? NAKASHIMA M.Hybrid simulation：an early history[J]. Earthquake Engineering & Structural Dynamics，2020，49（10）：949—962.

[4]? ? 郭玉榮，黃強(qiáng). 子結(jié)構(gòu)擬動(dòng)力試驗(yàn)邊界條件模擬方法[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，44（3）：68—75.

GUO Y R，HUANG Q. Boundary condition simulation method of substructure pseudo-dynamic test[J]. Journal of Hunan University （Natural Sciences），2017，44（3）：68—75. （In Chinese）

[5]? ? NAKASHIMA M，KATO H，TAKAOKA E. Development of real-time pseudo dynamic testing[J]. Earthquake Engineering & Structural Dynamics，1992，21（1）：79—92.

[6]? ? HORIUCHI T，INOUE M，KONNO T，et al. Real-time hybrid experimental system with actuator delay compensation and its application to a piping system with energy absorber[J]. Earthquake Engineering & Structural Dynamics，1999，28（10）：1121—1141.

[7]? ? CHEN C，RICLES J M. Improving the inverse compensation method for real-time hybrid simulation through a dual compensation scheme[J]. Earthquake Engineering & Structural Dynamics，2009，38（10）：1237—1255.

[8]? ? CHEN P C，CHANG C M，SPENCER B F，et al. Adaptive model-based tracking control for real-time hybrid simulation[J]. Bulletin of Earthquake Engineering，2015，13（6）：1633—1653.

[9]? ? PHILLIPS B M，SPENCER B F Jr. Model-based feedforward-feedback actuator control for real-time hybrid simulation[J]. Journal of Structural Engineering，2013，139（7）：1205—1214.

[10]? WU B，ZHOU H. Sliding mode for equivalent force control in real-time substructure testing[J]. Structural Control and Health Monitoring，2014，21（10）：1284—1303.

[11]? 寧西占，周惠蒙，吳斌，等.基于內(nèi)環(huán)H∞控制的實(shí)時(shí)混合試驗(yàn)[J]. 振動(dòng)與沖擊，2017，36（15）：57—63.

NING X Z，ZHOU H M，WU B，et al. Real-time hybrid simulation based on inner-loop H∞ control[J]. Journal of Vibration and Shock，2017，36（15）：57—63. （In Chinese）

[12]? WU B，WANG Z，BURSI O S. Actuator dynamics compensation based on upper bound delay for real-time hybrid simulation[J]. Earthquake Engineering & Structural Dynamics，2013，42（12）：1749—1765.

[13]? WANG Z，WU B，BURSI O S，et al. An effective online delay estimation method based on a simplified physical system model for real-time hybrid simulation[J]. Smart Structures and Systems，2014，14（6）：1247—1267.

[14]? 李寧，周子豪，李忠獻(xiàn). 基于時(shí)滯追蹤的實(shí)時(shí)混合試驗(yàn)自適應(yīng)補(bǔ)償方法[J]. 工程力學(xué)，2019，36（7）：38—47.

LI N，ZHOU Z H，LI Z X. Time-delay tracing based adaptive compensation algorithm for real-time hybrid testing[J]. Engineering Mechanics，2019，36（7）：38—47. （In Chinese）

[15]? NING X Z，WANG Z，ZHOU H M，et al. Robust actuator dynamics compensation method for real-time hybrid simulation[J]. Mechanical Systems and Signal Processing，2019，131：49—70.

[16]? 李騰飛，隋龑，蘇明周，等. 偏心支撐框架子結(jié)構(gòu)實(shí)時(shí)混合仿真試驗(yàn)研究[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2018，45（11）：46—53.

LI T F，SUI Y，SU M Z，et al. Study on real time hybrid simulation test of an eccentrically braced frame as test sub-structure[J]. Journal of Hunan University （Natural Sciences），2018，45（11）：46—53. （In Chinese）

[17]? ZHOU K，DOYLE J C，GLOVER K. Robust and optimal control [M]. Englewood Cliffs：Prentice Hall，1996：442—445.

[18]? 吳旭東，解學(xué)書. H∞魯棒控制中的加權(quán)陣選擇[J]. 清華大學(xué)學(xué)報(bào)（自然科學(xué)版），1997，37（1）：27—30.

WU X D，XIE X S. Weighting function matrix selection in H∞ robust control [J]. Journal of Tsinghua University （Sci & Tech），1997，37（1）：27—30. （In Chinese）

[19]? BONNET P A，LIM C N，WILLIAMS M S，et al. Real-time hybrid experiments with newmark integration，MCSmd outer-loop control and multi-tasking strategies[J]. Earthquake Engineering & Structural Dynamics，2007，36（1）：119—141.

[20]? SILVA C E，GOMEZ D，MAGHAREH A，et al. Benchmark control problem for real-time hybrid simulation[J]. Mechanical Systems and Signal Processing，2020，135：106381.

[21]? ZHAO J，F(xiàn)RENCH C，SHIELD C，et al. Considerations for the development of real-time dynamic testing using servo-hydraulic actuation[J]. Earthquake Engineering & Structural Dynamics，2003，32（11）：1773—1794.

[22]? TAO J J，MERCAN O. A study on a benchmark control problem for real-time hybrid simulation with a tracking error-based adaptive compensator combined with a supplementary proportional-integral-derivative controller[J]. Mechanical Systems and Signal Processing，2019，134：106346.

[23]? RUPPERT D，CLINE D B H. Bias reduction in kernel density estimation by smoothed empirical transformations[J]. The Annals of Statistics，1994，22（1）：185—210.