黃 亮 黃慎江 王靜峰 徐偉杰 郭 彤

(1合肥工業(yè)大學(xué)土木與水利工程學(xué)院，合肥 230009)(2合肥工業(yè)大學(xué)土木工程結(jié)構(gòu)與材料安徽省重點實驗室，合肥 230009)(3東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點試驗室，南京 211189)

抗震試驗是檢驗土木工程領(lǐng)域新理論、新結(jié)構(gòu)和新材料有效性的重要手段，也為抗震設(shè)計提供了重要的實踐支撐[1].由Nakashima等[2]提出的實時混合模擬試驗是一種在有限空間內(nèi)進行大尺度動力試驗的有效方法.該方法繼承了擬動力試驗思路，采用子結(jié)構(gòu)交互技術(shù)，實現(xiàn)了有限元計算分析和動力加載試驗的同步交互耦合，可同時反映整體結(jié)構(gòu)抗震性能和局部構(gòu)件的變形損傷，是今后土木工程實驗的發(fā)展方向之一.

實時混合模擬對實驗的軟硬件設(shè)備要求較高.硬件方面更新了高性能計算機、高速信息通訊設(shè)備及動態(tài)加載控制系統(tǒng)；軟件方面開發(fā)了新的積分算法和誤差修正方法[3].上述技術(shù)革新使得該試驗方法具有在真實時間尺度內(nèi)同步計算、同步加載、同步信息交互的能力.

實時混合模擬要求子結(jié)構(gòu)間的信息同步交互.然而，積分計算延遲、數(shù)字/模擬信號轉(zhuǎn)換、數(shù)據(jù)通訊、控制器采樣頻率和液壓作動器動力響應(yīng)延遲等[4]均會導(dǎo)致子結(jié)構(gòu)間信息交互不同步.其中，積分計算延遲和液壓作動器響應(yīng)延遲分別是導(dǎo)致數(shù)值子結(jié)構(gòu)有限元分析和試驗子結(jié)構(gòu)動力加載延遲的最主要原因，其余延遲極小，可忽略不計.

積分計算和作動器響應(yīng)延遲對實時混合模擬試驗的影響截然不同.前者會導(dǎo)致試驗停頓，破壞試驗的連續(xù)性；后者則會影響試驗結(jié)果的準(zhǔn)確性和試驗系統(tǒng)的穩(wěn)定性.目前的研究大多關(guān)注作動器延遲的影響評估、時滯誤差控制和補償方法等[5-7].

本文研究了積分計算延遲的產(chǎn)生機理和影響，提出減少計算延遲的方案，并采用試驗加以驗證.

1 實時混合模擬試驗

1.1 基本原理

實時混合模擬試驗將整體結(jié)構(gòu)視為研究對象，并劃分為數(shù)值與試驗子結(jié)構(gòu)2個部分.參考文獻[8]，試驗步驟如下：

① 劃分?jǐn)?shù)值與試驗子結(jié)構(gòu)，建立運動方程

(1)

② 選擇積分算法，求解運動方程，獲得第i+1步位移ui+1.

為了保證試驗的實時性，要求每一步試驗的執(zhí)行時間等于該步的積分步長.

1.2 積分算法和子步技術(shù)

MTS公司提供的實時混合模擬解決方案中，數(shù)值子結(jié)構(gòu)采用有限元軟件OpenSees進行分析，并提供Newmark[9-10]、HHT-α[11]、Generalized-α[12]和TRBDF2[13]等積分算法用于數(shù)值積分.為了兼顧算法精度和穩(wěn)定性，推薦使用隱式Generalized-α算法通過Newton迭代求解運動方程.一旦殘差收斂到閾值，便可獲得第i+1步動力響應(yīng).為了平衡計算效率和準(zhǔn)確性，積分步長Δt通常取控制步長δt的整數(shù)倍，且不超過0.02 s，殘差閾值通常取10-6～10-4.

上述方案中，計算耗時與殘差閾值相關(guān).閾值越小，隱式算法的收斂速度越慢，計算耗時越長.然而，實時試驗對同步性要求較高，不允許出現(xiàn)因計算延遲所導(dǎo)致的試驗停頓現(xiàn)象.因此，方案中引入了預(yù)測-修正子步技術(shù)，用于增大積分步長，延長計算時間.子步技術(shù)緩解了計算資源緊張的壓力，但無法徹底解決計算延遲問題.

圖1 預(yù)測-修正過程

(2)

(3)

式中

上述預(yù)測-修正過程可滿足絕大多數(shù)計算和加載的需求，但在某些時刻，隱式算法收斂速度緩慢導(dǎo)致計算耗時超限.當(dāng)控制系統(tǒng)在0.6Δt內(nèi)仍無法獲得下一步指令時，預(yù)測-修正模塊被迫進入等待階段，直至計算機完成積分計算，才重新進入修正狀態(tài).這一特殊過程即為預(yù)測-等待-修正過程.

2 滑移支座試驗

本節(jié)以在東南大學(xué)混凝土及預(yù)應(yīng)力混凝土結(jié)構(gòu)教育部重點實驗室完成的滑移支座實時混合模擬試驗為例，研究試驗中偶發(fā)計算延遲產(chǎn)生的原因及影響.

東晉醫(yī)學(xué)家葛洪就曾寫過一本醫(yī)書《肘后備急方》，書名的“肘后”一詞，就是說這本書可以放在肘后的口袋中隨身攜帶。“袖珍”一詞也是來源于此。而我們常用來比喻官員作風(fēng)清廉的“兩袖清風(fēng)”，其實原意就是說這個官員袖子里沒有裝錢。

2.1 試驗

滑移支座實時混合模擬試驗[15]研究了含新型滑移支座與鉛芯橡膠支座的混合隔震體系對4層鋼筋混凝土框架的隔震作用(見圖2).圖中，B柱下的滑移支座分割為試驗子結(jié)構(gòu)，其余支座和上部框架定義為數(shù)值子結(jié)構(gòu).梁柱均采用桿件單元進行有限元建模，數(shù)值子結(jié)構(gòu)共包含20個節(jié)點，56個自由度，僅底部4個節(jié)點約束了轉(zhuǎn)動自由度，其余16個節(jié)點均包含水平、垂直和轉(zhuǎn)動自由度.

(a) 基礎(chǔ)隔震結(jié)構(gòu)示意圖

(b) 試驗裝置照片

試驗使用的計算機核心設(shè)備包括Intel Core i5-680 3.6 GHz 芯片和2 048 MB 內(nèi)存.數(shù)據(jù)采集設(shè)備包括8通道數(shù)據(jù)采集器，采樣頻率為1 024 Hz.試驗系統(tǒng)采用Simulink平臺進行流程控制，通過閉環(huán)控制系統(tǒng)指揮作動器工作.

試驗共完成了El-Centro、Taft和Nanjing地震波作用下各5組實時混合模擬，試驗時長分別為60、60和20 s.

2.2 結(jié)果分析

經(jīng)15組重復(fù)試驗(編號分別為EL-1～EL-5、Taft-1～Taft-5和NJ-1～NJ-5)，發(fā)現(xiàn)其中9組試驗實現(xiàn)了完全實時同步，其余6組試驗出現(xiàn)計算延遲(見表1).統(tǒng)計結(jié)果顯示，計算延遲具有偶發(fā)性和隨機性.最嚴(yán)重的計算延遲出現(xiàn)在EL-2的第1 044步，共超時55δt；最輕微的計算延遲出現(xiàn)在Taft-2的第806步，計算超時僅δt.試驗時間尺度(即真實耗時與理想耗時的比值)均小于1.001.本文重點分析Taft-3.

表1 計算延遲統(tǒng)計

Taft-3的理想試驗耗時為60 s，實際超時71.3 ms，時間尺度為1.001.試驗中存在4處短暫停頓，分別在第131、1 768、1 869和2 483步，實際耗時分別為41δt、48δt、35δt和29δt.圖3顯示了第131步的試驗停頓現(xiàn)象.

(a) 計算延遲現(xiàn)象

(b) 第131步細(xì)節(jié)圖

(4)

式中

(5)

式中

3 計算延時的影響

計算延遲會對指令位移信號和作動器響應(yīng)產(chǎn)生影響.

3.1 指令位移信號的影響

3.2 作動器響應(yīng)的影響

受液壓伺服系統(tǒng)加載誤差影響，測量反饋力不僅與終點狀態(tài)有關(guān)，還與加載過程密切相關(guān).

Zhao等[16]對作動器位移追蹤進行掃頻試驗，發(fā)現(xiàn)位移追蹤效果與指令信號相關(guān)，信號頻率越高，位移追蹤效果越差.出現(xiàn)計算延遲時，指令位移信號呈S形，導(dǎo)致信號中含有高頻成分.由于作動

圖4 不同計算耗時情況下和

器在短時間內(nèi)被命令先減速等待再加速追趕，動力加載系統(tǒng)需克服管道內(nèi)液壓油、作動器加載頭、試驗構(gòu)件等慣性力作用，使得作動器追蹤出現(xiàn)顯著誤差.因此，S形指令位移信號對位移追蹤效果極為不利.

圖5(a)～(b)比較了Taft-3試驗中t1 766～t1 769與t1 533～t1 536時步的位移追蹤效果.由圖5(a)可知，t1 766～t1 769時步的位移追蹤延遲分別為0、-1、7.8、0 ms，其中負(fù)值表示該時步位移過補償.對于計算延遲出現(xiàn)的t1 768步，由于作動器被命令在48δt內(nèi)先減速再加速，導(dǎo)致該時步的位移追蹤延遲從-1 ms急劇增加到7.8 ms，出現(xiàn)了8.8 ms的響應(yīng)延遲，表明該時步作動器位移追蹤效果極差.由圖5(b)可知，t1 533～t1 536時步的位移追蹤延遲分別為1.9、2.3、2.2、1.5 ms，未出現(xiàn)計算延遲，故作動器位移追蹤效果較好且波動很小.對比發(fā)現(xiàn)，由于液壓作動器在短時間內(nèi)被要求急停急轉(zhuǎn)，其響應(yīng)會出現(xiàn)較大延遲，位移追蹤效果也大大降低.因此，需徹底避免計算延遲，防止作動器出現(xiàn)急停急轉(zhuǎn)的現(xiàn)象.

4 克服計算延遲的改進方案

為了防止偶發(fā)計算延遲、保障試驗的實時性和連續(xù)性，可通過使用顯式積分算法和簡化數(shù)值子結(jié)構(gòu)模型提高計算效率，降低計算耗時.

(a) 計算延遲步響應(yīng)

(b) 標(biāo)準(zhǔn)步響應(yīng)

降低計算耗時的最有效方法是引入顯式積分算法，用于代替需迭代計算的隱式積分算法.實時混合試驗中常用的顯式算法包括CR算法[17]、Chang算法[18]等.顯式積分算法通過遞推直接求解運動方程，計算時間可控，大幅降低了積分計算時間.

其次，簡化數(shù)值子結(jié)構(gòu)模型同樣可減少計算耗時.通過模型簡化可減少計算自由度，減少非線性單元數(shù)量，降低運動方程復(fù)雜度，進而減少計算耗時.

上述改進方案已在鉛芯橡膠支座實時混合模擬試驗[19]中得到驗證.該試驗采用相同的硬件設(shè)備，研究了二層鋼框架隔振結(jié)構(gòu)的隔震性能和支座的變形、耗能能力.試驗中，梁柱均采用線彈性桿件單元建模，數(shù)值子結(jié)構(gòu)僅保留樓板處3個水平自由度，積分算法采用CR顯式積分算法.選用EL-Centro地震波進行了12次重復(fù)試驗，每次試驗時長60 s，未發(fā)現(xiàn)計算延遲現(xiàn)象，所有試驗均保證嚴(yán)格實時同步.該試驗驗證了使用顯式積分算法和簡化數(shù)值子結(jié)構(gòu)模型可有效降低計算耗時，避免出現(xiàn)由計算延遲導(dǎo)致的實時試驗停頓現(xiàn)象.

5 結(jié)論

1) 在實時混合模擬中使用隱式積分算法時，試驗可能出現(xiàn)偶發(fā)計算延遲，該現(xiàn)象具有偶然性和隨機性，無法預(yù)先補償.

2) 計算延遲產(chǎn)生的原因是隱式積分算法通過迭代計算求解運動方程，其殘差收斂速度緩慢.

3) 計算延遲的影響是使作動器在短時間內(nèi)急停急轉(zhuǎn)，導(dǎo)致其位移追蹤效果大大降低.

4) 通過使用顯式積分算法、簡化模型等方法可避免出現(xiàn)計算延遲.