郭浩

摘 要：擬動力試驗方法是一種研究建筑結(jié)構(gòu)抗震性能的有效方法?？紤]試驗場地和費用等因素，子結(jié)構(gòu)擬動力試驗方法被提出，它是將建筑結(jié)構(gòu)中的關(guān)鍵部位進行試驗，其余部分通過有限元軟件OpenSees進行數(shù)值分析，試驗子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)通過OpenFresco平臺建立連接，實現(xiàn)數(shù)據(jù)的傳遞，不僅可以提高效率，還能降低費用，因此，擬動力試驗技術(shù)具有較好的研究價值和發(fā)展前景。

關(guān)鍵詞：擬動力試驗；有限元分析；試驗子結(jié)構(gòu)；數(shù)值子結(jié)構(gòu)

1 研究目的和意義

隨著我國經(jīng)濟及城市化的快速發(fā)展，越來越多的超高層建筑出現(xiàn)在各個地區(qū)。伴隨著大跨度、復(fù)雜的建筑結(jié)構(gòu)的出現(xiàn)，對工程結(jié)構(gòu)的抗震性能要求也越來越高。我國是世界上多地震國家之一，歷史上曾發(fā)生多次強烈地震，例如2010年的青海省玉樹縣的7.1級地震、2013年的四川省雅安市的7.0級地震、2014年新疆于田的7.3級地震，都造成巨大的經(jīng)濟損失和人員傷亡，因此提高建筑結(jié)構(gòu)的抗震減振性能就顯得尤為關(guān)鍵。結(jié)構(gòu)抗震試驗為評判結(jié)構(gòu)抗震性能、研究結(jié)構(gòu)抗震機理發(fā)揮不可替代的作用，受到科學(xué)研究者的廣大重視[1]。

當前抗震試驗主要有擬靜力試驗、地震模擬振動臺、擬動力試驗三種方式。擬動力試驗綜合擬靜力和地震模擬振動臺的優(yōu)點，可以較好的模擬大型結(jié)構(gòu)的地震反應(yīng)。但是出于實驗設(shè)備和費用的考慮，子結(jié)構(gòu)擬動力試驗方法被提出，建筑物發(fā)生嚴重破壞通常是局部的構(gòu)件，將建筑結(jié)構(gòu)中最易被破壞的部分進行試驗，其余部分進行計算機模擬，通過試驗和模擬的結(jié)合讓整體結(jié)構(gòu)在動力方程中得到統(tǒng)一[2]。本文將對擬動力試驗方法的起源、發(fā)展及研究近況進行簡要的介紹，分析擬動力試驗方法中存在的問題；然后對擬動力試驗方法的數(shù)值計算方法的穩(wěn)定性和精確性進行討論，最后對擬動力實驗方法未來的發(fā)展情況做出簡要的概述。

2 擬動力試驗的研究概況

2.1 擬動力試驗的定義

擬動力試驗方法自1969年由日本學(xué)者M.Hakuno等[3]首次提出后，經(jīng)過多年的發(fā)展，擬動力試驗由早期的單自由度試驗發(fā)展到后來的多自由度擬動力試驗，以及近些年來比較關(guān)注的子結(jié)構(gòu)擬動力試驗。子結(jié)構(gòu)擬動力試驗包括兩個部分：受力復(fù)雜且關(guān)鍵部分的試驗子結(jié)構(gòu)和其余部分的計算子結(jié)構(gòu)，該方法稱為結(jié)構(gòu)混合試驗，它采用開源有限元軟件OpenSees進行數(shù)值部分的計算，為了使計算結(jié)果容易收斂，可以選擇自由度少的構(gòu)件作為數(shù)值部分[4]，例如混凝土結(jié)構(gòu)；然后通過接口軟件OpenFresco實現(xiàn)數(shù)值模擬部分和實驗加載部分的結(jié)合，將有限元軟件與試驗控制系統(tǒng)MTS進行數(shù)據(jù)交換。混合試驗過程就是通過作動器給試驗結(jié)構(gòu)施加荷載獲取恢復(fù)力，利用客戶端軟件實現(xiàn)試驗子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)的連接，通過數(shù)值積分方法求解整體結(jié)構(gòu)的動力方程，獲得各個時間段上的每個自由度的位移、速度、加速度和內(nèi)力。該方法的提出解決試驗結(jié)構(gòu)因尺寸過大限制了加載條件的問題，范云蕾等[5]運用該方法對一座橋梁進行了試驗，驗證該方法的可行性。同時根據(jù)相關(guān)研究結(jié)果 [6]表明采用數(shù)值模擬方法對超高層建筑抗震性能的研究是一個有效的方法，當前開源有限元軟件已在高層和超高層的抗震性能研究中得到廣泛地運用。

2.2 擬動力試驗的應(yīng)用

由以往的歷次地震現(xiàn)象表明，建筑結(jié)構(gòu)在地震中倒塌的主要原因之一是鋼筋混凝土柱抵抗變形的能力不足而造成的，因此，探究鋼筋混凝土結(jié)構(gòu)的抗震性能是十分必要的。近年來研究人員利用你動力實驗對其抗震性能做了大量的試驗分析，以下將進行具體地闡述。

楊濤[7]等利用擬動力試驗方法對單跨兩層的框架縮尺模型施加不同峰值大小的地震作用，該框架結(jié)構(gòu)由鋼筋混凝土梁和T形截面鋼骨混凝土異形柱組成，試驗結(jié)果表明異形柱框架結(jié)構(gòu)具有良好的抗震性能；為研究框架結(jié)構(gòu)底層開間大小對結(jié)構(gòu)抗震性能的影響，唐岱新等[8]對底層大開間框剪組合墻的縮尺結(jié)構(gòu)進行子結(jié)構(gòu)擬動力試驗，通過檢測框架柱的側(cè)向剛度證明底層大開間對框架結(jié)構(gòu)的抗震性能有一定的影響，可以通過調(diào)整層間的彈性剛度比保證結(jié)構(gòu)的抗震性能；袁小龍等[9]利用有限元軟件對單層平面鋼框架結(jié)構(gòu)進行混合實驗，通過接口軟件實現(xiàn)數(shù)值計算部分和物理加載部分的連接，實現(xiàn)軟件和試驗控制系統(tǒng)之間的數(shù)據(jù)交換，結(jié)果表明該方法能夠有效地檢測品面鋼框架結(jié)構(gòu)在地震作用下的性能，結(jié)果的可靠性較好。

伴隨著結(jié)構(gòu)體系地復(fù)雜化，影響實驗結(jié)果的因素也在增加。結(jié)構(gòu)模型的縮尺是一項重要的內(nèi)容，在靜力相似理論的基礎(chǔ)上，陳再現(xiàn)等[10]提出一套基于靜力相似的縮尺模型擬動力試驗方法，降低傳統(tǒng)縮尺模型試驗的相似過程的復(fù)雜性，結(jié)果顯示用該方法得到的結(jié)果與足尺結(jié)構(gòu)的分析結(jié)果具有較高的吻合度，證實該方法的可行性；復(fù)雜的結(jié)構(gòu)進行擬動力試驗時，子結(jié)構(gòu)數(shù)量會對實驗結(jié)果產(chǎn)生影響，為減少影響，郭玉榮等[11]提出最小二乘法模型修正方法，并使用該方法對剪切模型和桿系防屈曲支撐模型進行時程分析，結(jié)果顯示使用該方法可以保證在子結(jié)構(gòu)數(shù)量有所限制的情況下，計算結(jié)果仍然具有較高的精度；為研究超高層剪力墻構(gòu)件的受力狀態(tài)，王麗莎等[12]提出一種高性能四邊形平板殼單元及其幾何非線性列式，將其用于分析剪力墻構(gòu)件在地震作用下的力學(xué)特性，通過對比有限元分析結(jié)果和試驗結(jié)果，驗證了該單元模型的合理性和可靠性。

2.2 擬動力試驗的發(fā)展

纖維模型是結(jié)構(gòu)非線性分析模型的研究熱點，它被廣泛地運用到研究鋼筋混凝土結(jié)構(gòu)的彈塑性分析中，它可以把截面劃分成不同部分，各個纖維根據(jù)用戶的需求來決定本構(gòu)關(guān)系，更接近與實際受力狀態(tài)。杜軻等[13]將開源有限元軟件中的三種非線性單元（基于位移插值的分布塑性單元、基于力插值的分布塑性單元和基于力插值的集中塑性鉸單元），用于計算懸臂柱的基底剪力和頂端位移，對比不同纖維單元的計算結(jié)果，得到各個纖維單元的劃分尺寸，確保得到最佳的計算結(jié)果；為準確描述鋼筋混凝土結(jié)構(gòu)的破壞形式，張行等[14]使用開源有限元軟件中的非線性纖維單元對四根鋼筋混凝土柱的抗震性能進行分析，選擇合適的材料本構(gòu)，對構(gòu)件施加低周反復(fù)荷載，試驗結(jié)果表明使用該類纖維模型可以提高模擬結(jié)果與試驗結(jié)果的吻合度，有效地反映構(gòu)件的加載后期剛度和強度的退化；王向英等[15]將子結(jié)構(gòu)擬動力技術(shù)應(yīng)用于振動臺試驗，提出試驗子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)之間的邊界條件，根據(jù)振動臺試驗的結(jié)果驗證該技術(shù)的可靠度。

伴隨著互聯(lián)網(wǎng)技術(shù)的發(fā)展，擬動力試驗進入到新的階段，通過網(wǎng)絡(luò)將各地的實驗室連接起來形成一個網(wǎng)絡(luò)實驗體系，增強信息的傳遞和資源的共享，可以對結(jié)構(gòu)復(fù)雜的體系進行擬動力試驗[16]，湖南大學(xué)和美國南加州福尼亞大學(xué)之間已經(jīng)、通過模擬實驗平臺開展橋柱和樁模在地震作用下的模擬實驗；加州大學(xué)伯克利分校開發(fā)的接口軟件實現(xiàn)限元軟件和電液伺服作動系統(tǒng)的連接，相關(guān)研究結(jié)果[17]證實該系統(tǒng)的可靠度，為實施多個子結(jié)構(gòu)試驗的遠程協(xié)同擬動力試驗提供保證。徐國山等[18]利用開源有限元軟件、接口軟件和試驗控制裝置進行混合試驗，試驗結(jié)果表明三者之間可以保持良好的通訊，計算結(jié)果具有較高的精度；賈紅星等[19]利用接口軟件實現(xiàn)有開源有限元軟件和電液伺服作動系統(tǒng)的連接，進行框架結(jié)構(gòu)的擬動力實驗，將結(jié)構(gòu)中柱子進行試驗加載，其余部分進行數(shù)值模擬，結(jié)果證實該方法是研究框架結(jié)構(gòu)抗震性能的可靠手段。

2.3 數(shù)值積分法

擬動力試驗中的數(shù)值積分部分通過運動微分方程進行求解，如式1所示，利用逐步積分算法獲得結(jié)構(gòu)在地震作用下各質(zhì)點處的位移，將計算結(jié)果施加給結(jié)構(gòu)，得到各質(zhì)點處的恢復(fù)力并及時反饋給結(jié)構(gòu)，最后根據(jù)恢復(fù)力和地震荷載進行下一時刻位移量的計算。常見的算法分為顯示算法和隱示算法，顯示算法有中央差分法、顯示紐馬克法和修正紐馬克法，隱示算法有Newmark-α法和Newmark-β法。顯示算法對積分步長有所限制，它的穩(wěn)定性和步長與自振周期有關(guān)，為保證算法的穩(wěn)定性，步長盡量取小，但不可過小，以免導(dǎo)致計算誤差增大；隱示算法是無條件穩(wěn)定，需進行迭代控制，但它的精度比顯示算法高。

Mai+Cvi+ri=fi （1）

式中M為結(jié)構(gòu)各質(zhì)點的質(zhì)量矩陣，ai為第i時刻的加速度，C為阻尼矩陣，vi為第i時刻的速度，ri為各質(zhì)點處的恢復(fù)力，fi為結(jié)構(gòu)各質(zhì)點處施加的外力。

（1）中央差分法

中央差分法是基于用有限差分代替位移對時間的求導(dǎo)，通過對位移方程求一階導(dǎo)和二階求導(dǎo)得速度和加速度。結(jié)構(gòu)的運動微分方程如式1所示，第i時刻的速度第i+1時刻和第i-1時刻位移的插值除以兩個時刻間的時間。將第i+1時刻的位移和第i+1時刻的速度代入結(jié)構(gòu)的運動微分方程，即可得到結(jié)構(gòu)下一時刻的位移。

（2）顯示紐馬克法

顯示紐馬克法是將時間離散化，將ti時刻的位移值作為初始值，通過逐步的積分算法，得到下一時刻的位移di+1和速度vi+1，位移di+1和速度vi+1的數(shù)值由兩個參數(shù)β和γ確定。

當β和γ系數(shù)不同時，對應(yīng)三種計算算法。當β=1/4，γ=1/2時，為平均加速度法；當β=1/6，γ=1/2時，為線性加速度法；當β=0，γ=1/2時，為顯示紐馬克法。將下一時刻的位移di+1和速度vi+1帶入動力方程，即可得到下一時刻的加速度。

（3）Newmark-β法

Newmark-β法是一種無條件穩(wěn)定的隱示算法，同樣是將時間離散化，運動方程只要在離散的時間點上滿足即可，通過對加速度的假設(shè)，利用第i時刻的位移獲得第i+1時刻的位移。通過聯(lián)立第i時刻的速度和位移表達式，可以得到第i+1時刻的速度和加速度計算公式，進一步獲得第i+1時刻的位移。

（4）其他

為提高數(shù)值算法的精度和計算效率，相關(guān)研究也在不斷進行，冉田苒等[20]提出一種新新型顯示算法，在隱式算法HHT-α法的基礎(chǔ)上，假定結(jié)構(gòu)作為線性系統(tǒng)進行數(shù)值分析，通過減少時間步長降低計算過程的誤差，使用該方法進行兩層無阻尼框架結(jié)構(gòu)的數(shù)值計算，結(jié)果證實該數(shù)值算法具有良好的可靠性和精度；鞠高云[21]提出一種基于改進分段二次插值的精細積分單步法和多步法，通過內(nèi)插點構(gòu)造基于改進分段二次插值的精細積分單步法和通過外插點構(gòu)造基于改進分段二次插值的精細積分多步法，并通過算例驗證了該方法的可行性。

3 擬動力試驗的缺陷

雖然擬動力實驗技術(shù)為研究建筑結(jié)構(gòu)的抗震性能提供便利，但還存在一些缺陷，例如，地震作用作為一種難以預(yù)測其變化規(guī)律的環(huán)境荷載，擬動力試驗還不能完全地反映建筑結(jié)構(gòu)在地震作用下的受力過程，這是擬動力試驗的一個缺陷，不能真實地反映結(jié)構(gòu)動力反應(yīng)的時間效應(yīng)；其次擬動力試驗中的速度和加速度是由數(shù)值計算得出的，與實際荷載有一定的差距。

4 總結(jié)

綜上所述，擬動力試驗技術(shù)結(jié)果是可靠的，它能真實地反映結(jié)構(gòu)在地震作用下動力特性，因此擬動力試驗方法可以在結(jié)構(gòu)抗震等許多領(lǐng)域發(fā)揮重要的作用。擬動力試驗作為一種新型的試驗方式，仍存在一些局限性，近年來擬動力試驗技術(shù)得到快速的發(fā)展，并出現(xiàn)了實時、混合、遠程等新型試驗技術(shù)，但其還處于初步發(fā)展階段，大量的研究和理論還有待于進一步探討的完善。

參考文獻：

[1]倪志偉，紀金豹，杜陽陽.擬動力試驗方法的發(fā)展與展望[J].災(zāi)害學(xué)，2010，25（b10）：79-84.

[2]邱法維.結(jié)構(gòu)抗震試驗方法進展 [J].土木工程學(xué)報，2004，37（10）：19-27.

[3] Hakuno M，Shidawara M，Hara T.Dynamic Destructive Test of a Cantilever Beam，Controlled by an Analog-Computer[J].2010，1969（171）：1-9.

[4]王瑾，蔡新江，田石柱.基于OpenSees的CFRP加固RC短柱抗震性能數(shù)值模擬[J].世界地震工程，2010（4）：108-111.

[5]范云蕾，肖巖，郭玉榮.異地分布式橋梁結(jié)構(gòu)遠程協(xié)同擬動力試驗[J].鐵道科學(xué)與工程學(xué)報，2013， 10（5）：28-33.

[6]解琳琳，黃羽立，陸新征.基于OpenSees的RC框架-核心筒超高層建筑抗震彈塑性分析[J].工程力學(xué)，2014，31（1）：64-71.

[7]楊濤，張喜德，鐘海牛，鄧志恒.T形截面鋼骨混凝土異形柱框架抗震性能擬動力試驗[J].中南大學(xué)學(xué)報（自然科學(xué)版），2015，46（6）：2711-2717.

[8]唐岱新，朱本全，王鳳來，等.底層大開間框剪組合墻結(jié)構(gòu)1/3比例模型房屋子結(jié)構(gòu)擬動力試驗研究[J].建筑結(jié)構(gòu)，2001（4）：20-22.

[9]袁小龍，韓曉健.基于OpenSees的混合試驗方法[J].混凝土與水泥制品，2014（10）：72-75.

[10]陳再現(xiàn)，陳芍橋，吳斌，王煥定.基于靜力相似的縮尺擬動力試驗方法數(shù)值模擬研究[J].土木工程學(xué)報，2014（s2）：307-311.

[11]郭玉榮，王超.子結(jié)構(gòu)擬動力試驗的模型修正法[J].科學(xué)技術(shù)與工程，2016，16 （22）：88-94.

[12]王麗莎，岑松，解琳琳，等.基于新型大變形平板殼單元的剪力墻模型及其在OpenSees中的應(yīng)用[J].工程力學(xué)，2016，33（3）：47-54.

[13]陳哲明.三層新型裝配式平面鋼框架擬動力試驗研究[D].哈爾濱：哈爾濱工業(yè)大學(xué)，2014：10-30.

[14]張行，張謝東.鋼筋混凝土柱擬靜力試驗數(shù)值模擬分[J].工程抗震與加固改造，2013，35（5）：30-35.

[15]王向英，田石柱.子結(jié)構(gòu)地震模擬振動臺混合試驗原理與實現(xiàn)[J].地震工程與工程振動，2009， 29（4）：46-52.

[16]肖巖，胡慶，郭玉榮.結(jié)構(gòu)擬動力遠程協(xié)同試驗網(wǎng)絡(luò)平臺的開發(fā)研究[J].建筑結(jié)構(gòu)學(xué)報， 2005， 26（3）：122-129.

[17]王亞明，劉煦， 袁國鋒.遠程協(xié)同子結(jié)構(gòu)擬動力試驗研究[J].低溫建筑技術(shù)，2012，34（7）：61-63.

[18]許國山，郝偉，陳永盛，等.基于有限元軟件OpenSees和OpenFresco的混合試驗[J].土木工程學(xué)報，2012（s1）：36-41.

[19]賈紅星，蔡新江，田石柱.應(yīng)用OpenFresco平臺的抗震混合試驗技術(shù)[J].四川建筑科學(xué)研究，2013， 39（4）：209-212.

[20]冉田苒，王濤，周惠蒙，等.一種無條件穩(wěn)定的顯式數(shù)值積分算法[J].防災(zāi)減災(zāi)工程學(xué)報，2016， v.36（1）：50-54.

[21]鞠高云.精細積分方法及其擬動力試驗算法研究[D].中南大學(xué)，2014.

基金項目：華僑大學(xué)研究生科研創(chuàng)新能力培育計劃資助項目（1611304008）