基于近完全時滯補償?shù)碾[式實時混合試驗方法

王 貞， 劉進進， 吳 斌， 許國山

(1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點實驗室， 黑龍江 哈爾濱 150090; 2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院， 黑龍江 哈爾濱 150090; 3.蘇州設(shè)計研究院股份有限公司， 江蘇 蘇州 215021)

基于近完全時滯補償?shù)碾[式實時混合試驗方法

王 貞1,2， 劉進進3， 吳 斌1,2， 許國山1,2

(1.哈爾濱工業(yè)大學(xué)結(jié)構(gòu)工程災(zāi)變與控制教育部重點實驗室， 黑龍江 哈爾濱 150090; 2.哈爾濱工業(yè)大學(xué)土木工程學(xué)院， 黑龍江 哈爾濱 150090; 3.蘇州設(shè)計研究院股份有限公司， 江蘇 蘇州 215021)

實時混合試驗是一種新型結(jié)構(gòu)抗震混合試驗方法。隱式逐步積分算法雖然具有穩(wěn)定性好的特點，但在實時混合試驗中實施困難。在分析兩種多自由度隱式實時混合試驗方法的基礎(chǔ)上，結(jié)合近完全時滯補償方法，提出了一種新型隱式實時混合試驗方法，并分析了該方法的性能。數(shù)值模擬表明，該方法具有較高的收斂速度和計算精度，能夠同時考慮時滯補償，能滿足自由度數(shù)目較多的多自由度結(jié)構(gòu)實時混合試驗的要求。

實時混合試驗； 隱式逐步積分算法； 時滯補償； 迭代

引 言

傳統(tǒng)的抗震試驗方法有三類，即擬靜力、擬動力和振動臺。近些年，擬動力試驗方法得到了快速發(fā)展，出現(xiàn)了子結(jié)構(gòu)化、網(wǎng)絡(luò)化和實時化的趨勢，分別產(chǎn)生了子結(jié)構(gòu)擬動力試驗方法、網(wǎng)絡(luò)分布式擬動力試驗和實時混合試驗(有時也稱實時子結(jié)構(gòu)試驗或者動力子結(jié)構(gòu)試驗等[1-4])。實時混合試驗方法把原型結(jié)構(gòu)分割為試驗子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)，分別通過真實試驗和實時計算進行模擬，二者之間的相互作用通過液壓伺服加載系統(tǒng)的實時加載保證。試驗子結(jié)構(gòu)往往是具有很強非線性或者力學(xué)性能時間相關(guān)的構(gòu)件或子結(jié)構(gòu)，數(shù)值子結(jié)構(gòu)是原型結(jié)構(gòu)的其余部分。實時混合試驗由于實時執(zhí)行加載和數(shù)值仿真，能反應(yīng)試件的速度/加速度相關(guān)力學(xué)性能，彌補了傳統(tǒng)擬動力試驗方法的不足，近些年得到了快速發(fā)展。為了保證實時加載[5-6]和實時計算[7]，尤其前者，使得實時混合試驗具有更大挑戰(zhàn)。

實時混合試驗需要在線實時求解結(jié)構(gòu)的運動方程，對逐步積分算法提出了較高要求。常用的顯式算法，如中心差分法，由于其條件穩(wěn)定性而可能在保證計算實時性時無法滿足對積分步長的限制條件。目前已經(jīng)提出的基于隱式的實時混合試驗方法，主要有兩種，即吳斌等人提出的等效力控制方法[8-9]和Shing等人提出的實時迭代混合試驗方法[10]。此兩種方法取得了一定的成功，但也存在一定的局限性。本文在分析它們局限性的基礎(chǔ)上，結(jié)合作者提出的時滯補償方法[11-12]，提出了一種新的隱式實時混合試驗方法，研究表明該方法對于多自由度結(jié)構(gòu)實時混合試驗具有較明顯的優(yōu)勢。

1 現(xiàn)有隱式實時混合試驗方法

實時混合試驗離散運動方程和基于常加速度假定的位移、速度表達式為：

MNai+1+CNvi+1+RN(di+1)+RE(ai+1，vi+1，di+1)=

Fi+1

(1)

(2)

(3)

式中 角標(biāo)N表示數(shù)值子結(jié)構(gòu)，角標(biāo)E表示試驗子結(jié)構(gòu)；M為結(jié)構(gòu)質(zhì)量矩陣；C為結(jié)構(gòu)阻尼矩陣；R為子結(jié)構(gòu)恢復(fù)力向量；F為外部荷載向量；d，v，a分別為結(jié)構(gòu)位移、速度、加速度向量；Δt為逐步積分時間間隔。

對式(2)和式(3)進行整理，用di+1和上一步結(jié)構(gòu)響應(yīng)來表示vi+1和ai+1，并代入結(jié)構(gòu)的運動方程(1)，可得

RN(di+1)+KPDdi+1+RE(ai+1，vi+1，di+1)=FEQ，i+1

(4)

其中

(5)

(6)

式(4)是一個關(guān)于di+1的非線性方程，在采用牛頓迭代法求解時，有：

(7)

(8)

(1)迭代收斂速度的波動導(dǎo)致作動器運動不平穩(wěn)，容易對作動器造成損傷且得不到試件的真實力反應(yīng)；

(2)不易在線獲得結(jié)構(gòu)剛度，矩陣求逆耗時較大。

針對隱式實時混合試驗存在的問題，Shing等人[10]提出了以實時迭代、插值發(fā)送命令等為特征的方法，如圖1所示。主要內(nèi)容包括：

(9)

(3)考慮到固定迭代次數(shù)及作動器時滯等因素帶來的誤差，在每積分步最后一個迭代步進行修正。取積分步末點處的位移di+1為

(10)

試驗子結(jié)構(gòu)的恢復(fù)力修正為

(11)

Shing方法取得了一定的成功，但仍然存在如下問題：1)通過對迭代位移進行插值發(fā)送命令，雖然解決了作動器加載速度突變的問題，但改變了下一次迭代計算的起點，不僅使得迭代過程變得復(fù)雜，收斂速度也大幅下降，最后一迭代步可能仍然存在較大誤差。2)在積分步長內(nèi)每個采樣時刻點迭代一次，意味著必須在一個控制采樣步長內(nèi)完成迭代計算、位移命令的內(nèi)插并發(fā)送給作動器，當(dāng)結(jié)構(gòu)自由度較多時，很可能無法滿足實時性的要求。3)即使迭代獲得精度很高的位移，位移命令進行加載時不可避免地涉及系統(tǒng)時滯的影響，引入新的誤差。

圖1 位移命令插值示意圖[10]Fig.1 Schematic for displacement command interpolation[10]

等效力控制方法[8-9]把式(4)等號左側(cè)看作是3個并聯(lián)元件，采用力反饋控制環(huán)取代了數(shù)值迭代來求解該方程，其原理如圖2所示。圖中CF稱為力-位移轉(zhuǎn)化矩陣，功能是把力信號轉(zhuǎn)化成位移信號，實現(xiàn)作動器的位移控制加載。采用該方法執(zhí)行實時混合試驗的流程是，首先得到第i+1步的等效力，內(nèi)插得到積分步內(nèi)的等效力命令，并在對應(yīng)的時刻發(fā)送給控制回路，作動器受相應(yīng)的位移命令驅(qū)動；在積分步長結(jié)束時，把試件的反力反饋到式(4)，再次求解位移，以降低控制誤差造成的影響。

可見，等效力控制方法創(chuàng)造性地提出采用控制的思路求解非線性方程，由于引入了控制的思想，使其具有鮮明的個性和突出的特點。分析表明，采用積分控制器并且階躍發(fā)送命令的等效力控制方法類似于傳統(tǒng)的迭代方法。合理設(shè)計的控制器可以較好地考慮等效力命令插值帶來的影響，而數(shù)值迭代則很難實現(xiàn)這一點。此外，由于采用了力反饋控制環(huán)，等效力控制方法可以同時考慮到系統(tǒng)作動性能和迭代過程，在一定程度上減小了系統(tǒng)時滯。不過，該方法存在以下缺點：1)內(nèi)插發(fā)送等效力命令，雖然能保證作動器運動的連續(xù)性，但會導(dǎo)致系統(tǒng)的控制性能變差。2)雖然合理設(shè)計的等效力控制器可以有效減小系統(tǒng)時滯，但是有時仍需要其他時滯補償方法來補償系統(tǒng)時滯。3)積分步內(nèi)固定計算次數(shù)，要求數(shù)值子結(jié)構(gòu)靜反力和等效力偏差的求解等過程必須在一個控制采樣步長內(nèi)完成，尤其是最后一個計算步內(nèi)，還要完成位移的修正計算。當(dāng)數(shù)值子結(jié)構(gòu)自由度較多時，矩陣運算復(fù)雜耗時，很難滿足實時混合試驗實時性的要求。

2 基于近完全時滯補償方法的隱式實時混合試驗方法

縱觀前述兩種隱式實時混合試驗方法，它們的主要缺點是：1)計算耗時：每一計算步復(fù)雜耗時的矩陣計算必須在一個控制采樣步長內(nèi)完成，限制了數(shù)值子結(jié)構(gòu)的復(fù)雜程度；2)計算精度：每一采樣步的計算以當(dāng)前所實現(xiàn)的位移為基礎(chǔ)，計算收斂速度有很大程度的下降，在固定次數(shù)計算后可能仍然存在較大誤差，該誤差在試驗過程中傳播、累積，對試驗結(jié)果可能產(chǎn)生較大影響；3)系統(tǒng)時滯：即使一個積分步長內(nèi)的計算精度很高，將位移命令發(fā)送給作動器進行加載時，不可避免地涉及到時滯，因此會受到時滯補償效果的影響。

鑒于此，提出以近完全時滯補償方法為基礎(chǔ)的隱式實時混合試驗方法。近完全時滯補償方法[11,12]的主要特點是借助過預(yù)測結(jié)構(gòu)位移響應(yīng)，使得試驗子結(jié)構(gòu)的位移超前于數(shù)值子結(jié)構(gòu)，只需在已實現(xiàn)的實測位移數(shù)據(jù)中搜索與期望位移最接近的位移值，同一時刻測得的試件子結(jié)構(gòu)反力即為試驗子結(jié)構(gòu)與期望位移對應(yīng)的反力。

在進行隱式實時混合試驗非線性方程的迭代求解時，若采用近完全時滯補償方法進行時滯補償，由于試驗子結(jié)構(gòu)的位移響應(yīng)超前于數(shù)值子結(jié)構(gòu)，迭代計算得到位移后，無需發(fā)送給作動器進行加載，而是直接從當(dāng)前完成的實測位移數(shù)據(jù)中搜索得到與迭代位移近似相等的試驗子結(jié)構(gòu)位移值，同一時刻的試驗子結(jié)構(gòu)實測反力即為與當(dāng)前迭代位移相對應(yīng)的反力值，將其直接帶入下一迭代步進行計算即可。圖3是基于近完全時滯補償方法[11-12]的隱式實時混合試驗方法原理示意圖。下面以試驗子結(jié)構(gòu)為彈簧試件的單自由度結(jié)構(gòu)體系為基礎(chǔ)，對照圖3詳細介紹基于近完全時滯補償方法的隱式實時混合試驗方法的流程(后文統(tǒng)一稱之為“新方法”)：

圖2 基于等效力控制方法的實時混合試驗原理Fig.2 Principle of real-time hybrid testing based on equivalent force control method

圖3 基于近完全時滯補償方法的隱式實時混合試驗原理Fig.3 Principle of implicit real-time hybrid testing based on nearly-complete delay compensation method

1)位移過預(yù)測：ti時刻預(yù)測ti+1+τp時刻的位移，如采用顯式Newmark表達式，有dti+1+τp=di+vi(τp+Δt) +0.5ai(τp+Δt)2，τp=τa+Δτ，τa為預(yù)估的作動器最大時滯，Δτ為過預(yù)測時間跨度；

2)內(nèi)插生成作動器位移命令：在dti+1+τp與ti時刻的作動器位移命令dti+τp之間，內(nèi)插得到作動器位移命令dc(ti～ti+1)，并發(fā)送給作動器。

3)采集數(shù)據(jù)：采集試驗子結(jié)構(gòu)位移響應(yīng)de(ti～ti+1)及試件反力RE(ti～ti+1)。

4)迭代求解非線性方程：

①每積分步第一次迭代以上一步的位移響應(yīng)為初值，其余次迭代以上次迭代結(jié)果為初值;

④進行多次迭代，直到收斂；更新位移預(yù)測所需的位移di+1、速度vi+1、加速度ai+1；當(dāng)試驗時間達到ti+1，開始下一個積分步的計算。

可見，基于近完全時滯補償方法的隱式實時混合試驗方法具有如下特點：1)采用簡化牛頓迭代，僅需試驗前計算逆矩陣；2)第k+1迭代步的計算完全基于第k迭代步的位移及其響應(yīng)進行，實現(xiàn)了標(biāo)準(zhǔn)的簡化牛頓迭代，其迭代收斂速度要優(yōu)于Shing方法和等效力控制方法，計算效率更高；3)積分步長內(nèi)無固定迭代次數(shù)，位移迭代達到收斂精度即可，而收斂迭代次數(shù)預(yù)期少于前兩種方法；4)在積分步長內(nèi)作動器的加載不依賴于當(dāng)前迭代位移值，只要在一個積分步長內(nèi)完成迭代收斂計算即可，對單步計算時間的限制比前兩種方法寬松，能適應(yīng)自由度數(shù)目更多的結(jié)構(gòu)；5)作動器的加載基于預(yù)測位移內(nèi)插生成位移命令進行，不會出現(xiàn)加載速度失真現(xiàn)象，并能保證試件加載的實時性；6)在迭代求解非線性方程的同時已經(jīng)考慮了作動器加載過程中的時滯，并進行了時滯補償；7)該方法的不足是，采取了位移過預(yù)測方法，預(yù)測的時間跨度更大，帶來了更大的預(yù)測誤差。該誤差可以通過預(yù)測精度更好的位移預(yù)測方法并通過回搜過程來減小。此外，還可通過其他方法先降低系統(tǒng)時滯，以減少預(yù)測跨度過大帶來的預(yù)測誤差。

3 數(shù)值模擬

3.1 單積分步計算精度

數(shù)值分析表明，對于不考慮時滯的小步長分析，三種方法都具有較好的精度?？紤]到多自由度分析中存在高頻成分、實時計算要求采用較大步長，采取如下參數(shù)：采樣步長δt=0.004 s，積分步長Δt=0.04 s;采用積分等效力控制器，控制參數(shù)KI=400；結(jié)構(gòu)參數(shù)為：KE=2×103N/m，MN=1 kg，KN=2×103N/m。三種方法單積分步位移計算結(jié)果如圖4所示?？梢?，在積分步長較大或結(jié)構(gòu)頻率較高時，即使在最后一個計算步對積分點位移進行了修正，Shing方法和EFCM方法的收斂精度也不是很理想。

圖4 單積分步計算精度比較Fig.4 Computational accuracy comparisons in single integration step

時滯補償是實時混合試驗最關(guān)鍵的問題，鑒于時滯可能發(fā)生波動，采用常規(guī)補償方法很難實現(xiàn)實測的完全補償，接下來的分析中假定存在約1 ms的欠補償時滯。作動器模型取為一階傳遞函數(shù)TA(s)=700/(700+s)，對于頻率為1 Hz的信號，該系統(tǒng)時滯約為1 ms。取積分步長Δt=0.01 s，不同方法的單積分步計算精度如圖5所示?？梢?，時滯得不到完全補償時，Shing方法和EFCM的單積分步計算精度低于新方法。

由前述分析可知，數(shù)值子結(jié)構(gòu)自由度較多或結(jié)構(gòu)頻率較高時，Shing方法和EFCM的計算精度會有所下降；實際試驗的欠補償時滯也會使得兩種方法計算精度下降。即使進行了位移修正，計算精度也不如所提出的基于簡化牛頓迭代的新方法。

圖5 時滯欠補償時單積分步位移計算精度比較Fig.5 Computational accuracy comparisons in single integration step with under-compensated delay

3.2 計算耗時

本小節(jié)比較應(yīng)用于多自由度結(jié)構(gòu)實時混合試驗時Shing方法、EFCM和新方法的計算耗時情況。取如圖6所示的5層剪切型結(jié)構(gòu)為分析對象，層間剛度為ki=40×103N/m，集中質(zhì)量為mi=100 kg(i=1～5)，結(jié)構(gòu)自振頻率為0.91，2.64，4.17，5.36，6.11 Hz，結(jié)構(gòu)無阻尼。選底層柱為試驗子結(jié)構(gòu)，剛度取kE=k1/2。采用理想作動器模型進行數(shù)值模擬，積分步長取Δt=0.01 s，每個積分步計算10次。新方法的計算精度取10-7m。計算平臺參數(shù)如表1所示。

圖6 5自由度結(jié)構(gòu)計算示意圖Fig.6 Computational schematic of 5DOF structure

單積分步內(nèi)的各迭代步耗時如圖7所示?？梢?，3種方法單迭代步的耗時并無太大差別，均約為1×10-5s。3種方法的第一迭代步耗時較其余步更大。EFCM在最后一個迭代步要完成位移命令計算和位移修正計算，計算耗時增加。新方法由于收斂速度較快，只進行3次迭代計算即達到精度要求，完成了該積分步內(nèi)方程求解。在實時混合試驗中，Shing方法和EFCM方法要求采樣步長大于單步迭代最大耗時，新方法的迭代次數(shù)遠小于它們，因此新方法大大提高了計算效率，更適用于多自由度結(jié)構(gòu)實時混合試驗。進一步的時程耗時分析見文[13]。需要說明的是，計算耗時統(tǒng)計是在性能較高的計算機上完成，各計算步耗時均較小。實際試驗中，控制板不僅要完成數(shù)值子結(jié)構(gòu)的計算、時滯補償計算，還要完成液壓伺服作動器的控制計算，多任務(wù)的同步實時處理，使得所能完成的實時混合試驗的自由度數(shù)目大大受限。

圖7 單積分步計算耗時Fig.7 Time consumed for computation in single integration step

表1 計算平臺參數(shù)

3.3 時程分析

為了進一步檢驗3種方法的性能，分別完成了無阻尼自由振動、無阻尼地震激勵、非線性自由振動、非線性地震激勵等多種工況的數(shù)值模型。限于篇幅，此處僅介紹最后一種工況的數(shù)值模擬。

分析模型為圖6所示的5層剪切型結(jié)構(gòu)，假設(shè)試驗子結(jié)構(gòu)具有雙線性恢復(fù)力關(guān)系，起始剛度為kE= 20×103N/m，屈服位移為0.008 m，屈服后剛度為起始剛度的0.115倍。首層自由度初始位移為0，初始速度為0；其他自由度的初始狀態(tài)為0。外荷載激勵取El Centro(NS，1940)地震記錄，峰值加速度調(diào)為0.2g。

考慮未能完全補償時滯的影響，即假設(shè)欠補償時滯約為1 ms，取作動器模型為TA(s)=700/(700+s)。等效力控制器參數(shù)不變。首層自由度位移時程及恢復(fù)力-位移關(guān)系見圖8和9所示?？梢娦路椒ǖ挠嬎憔纫糜谄渌麅煞N方法。

圖8 位移時程Fig.8 Displacement time histories

圖9 試件恢復(fù)力-位移曲線Fig.9 Restoring force-displacement relation of specimen

4 結(jié) 論

分析了Shing方法和等效力控制方法(EFCM)在實時混合試驗中可能存在的困難，提出了基于近完全時滯補償方法的隱式實時混合試驗方法，詳細闡述了其原理和實現(xiàn)流程，并從理論上分析了其特點。最后，通過數(shù)值模擬比較了三種方法在收斂速度和計算精度等方面的性能，驗證了理論分析結(jié)果。主要結(jié)論如下：

1)Shing方法和EFCM應(yīng)用于多自由度結(jié)構(gòu)實時混合試驗時，可能會由于自由度數(shù)目的增多，使得運算耗時增加，難以滿足實時性的要求；

2)當(dāng)積分步長增大或試驗子結(jié)構(gòu)的初始剛度較大，以及試驗中可能存在時滯欠補償情況時，Shing方法和EFCM的計算精度會有所下降；

3)基于近完全時滯補償方法的隱式實時混合試驗方法，具有較高的計算精度和收斂速度，能夠同時考慮時滯補償，能滿足較多自由度數(shù)目的多自由度結(jié)構(gòu)實時混合試驗的要求。

[1] Nakashima M, Kato H, Takaoka E. Development of real-time pseudo-dynamic testing[J]. Earthquake Engineering & Structure Dynamics, 1992,21(1):79—92.

[2] 吳斌,王倩穎.實時子結(jié)構(gòu)實驗的研究進展[J].實驗力學(xué),2007,22(5):1—9.

Wu Bin, Wang Qianying. Development of real-time substructure testing[J]. Journal of Experimental Mechanics, 2007,22(5):1—9.

[3] 遲福東,王進廷,金峰,等.土-結(jié)構(gòu)動力相互作用的實時耦聯(lián)動力試驗的時滯穩(wěn)定性[J].工程力學(xué),2012,29(8):1—7.

Chi Fudong, Wang Jinting, Jin Feng, et al. Delay-dependent stability of real-time dynamic hybrid testing for soil-structure interaction analysis[J]. Engineering Mechanics, 2012,29(8):1—7.

[4] 蔡新江,田石柱.振動臺試驗方法的研究進展[J].結(jié)構(gòu)工程師,2011,27(增刊):42—46.

Cai Xinjiang, Tian Shizhu. Research advances of shaking tabel testing method[J]. Structural Engineers, 2011,27(Sup):42—46.

[5] Horiuchi T, Inoue M, Konno T， et al. Namita. 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):1 121—1 141.

[6] 王貞,吳斌.基于最小二乘法的時滯實時在線估計方法[J].振動工程學(xué)報,2009,22(6):625—631.

Wang Zhen, Wu Bin. A real-time approach to delay estimation based on the Least-Square algorithm[J]. Journal of Vibration Engineering, 2009,22(6):625—631.

[7] 賈傳果,李英民,劉立平，等.Rosenbrock耦合積分方法及其收斂性分析[J].重慶大學(xué)學(xué)報,2013,36(2):127—133.

Jia Chuanguo, Li Yingmin, Liu Liping, et al. Rosenbrock-based coupled time integration method and its convergence analysis[J]. Journal of Chongqiong University, 2013,36(2):127—133.

[8] Wu B, Wang Q, Shing P B, et al. Equivalent force control method for generalized real-time substructure testing with implicit integration[J]. Earthquake Engineering & Structural Dynamics, 2007,36:1 127—1 149.

[9] 許國山,吳斌.采用等效力控制方法的非線性結(jié)構(gòu)實時子結(jié)構(gòu)試驗[J].振動工程學(xué)報,2010,23(3):237—242.

Xu Guoshan, Wu Bin. Real-time substructure testing of nonlinear structures with equivalent force control method[J]. Journal of Vibration Engineering, 2010,23(3):237—242.

[10]Jung R, Shing P, Stauffer E, et al. Performance of a real-time pseudo-dynamic test system considering nonlinear structural response[J]. Earthquake Engineering & Structural Dynamics, 2007,36(12):1 785—1 809.

[11]Bin Wu, Zhen Wang, Oreste S. Bursi. Actuator Dynamics Compensation Based on Upper Bound Delay for Real-time Hybrid Simulation[J]. Earthquake Engineering & Structural Dynamics, 2013,42(12):1 749—1 765.

[12]王貞,劉進進,吳斌.實時混合試驗近完全時滯補償方法的參數(shù)確定[J]. 工程力學(xué)，2014,31(10):158—166.

Wang Zhen, Liu Jinjin, Wu Bin. Parameter determination of Nearly-complete compensation scheme for time delay in real-time hybrid simulation[J]. Engineering Mechanics, 2014,31(10):158—166.

[13]劉進進.近完全時滯補償方法及其在實時混合試驗中的應(yīng)用[D].哈爾濱：哈爾濱工業(yè)大學(xué), 2013.

Liu Jinjin. Nearly complete compensation scheme for time delay and its application to real-time hybrid testing[D]. Harbin: Harbin Institute of Technology, 2013.

Implicit method for real-time hybrid testing based on nearly-complete delay compensation scheme

WANGZhen1,2,LIUJin-jin3,WUBin1,2,XUGuo-shan1,2

(1.Key Lab of Structures Dynamic Behavior & Control (Harbin Institute of Technology), Ministry of Education,Harbin 150090, China; 2.School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China;3.Suzhou Institute of Architectural Design Co., LTD., Suzhou 215021, China)

Real-time Hybrid Testing (RHT) is a novel hybrid testing method for evaluating structural seismic performance. Even though implicit time-stepping approaches exhibit better stability, their applications to RHT are limited due to their difficulties. This paper firstly analyses limitations of two existing implicit RHT schemes for Multi-Degree of Freedom (MDoF) structures, and then proposes a new implicit RHT based on the nearly-complete delay compensation developed by the authors formerly. Theoretical performance analysis and numerical simulations are performed following that. Numerical simulations show that the proposed scheme is endowed with better convergence speed and accuracy and incorporates delay compensation, which indicate that it is more suitable to RHT for MDoF systems.

real-time hybrid testing; implicit time-stepping method; delay compensation; iteration

2013-12-01；

2014-07-14

國家自然科學(xué)基金國際合作項目(51161120360)和重點項目(91315301-09);中國博士后科學(xué)基金面上資助項目(2013M531046);中央高?；究蒲袠I(yè)務(wù)費專項資金資助項目(HIT.NSRIF.2014102);哈爾濱工業(yè)大學(xué)海外博士科研啟動項目

TU317； P315.9

A

1004-4523(2015)03-0418-07

10.16385/j.cnki.issn.1004-4523.2015.03.011

王貞(1983—), 男, 博士，講師。電話： 15004548396； E-mail: zhenwang@hit.edu.cn