許成祥，胡鴻運

(武漢科技大學(xué)城市建設(shè)學(xué)院，湖北 武漢,430065)

擬靜力、擬動力及地震振動臺試驗是三種主要的結(jié)構(gòu)抗震試驗方法[1]，但通常因試驗設(shè)備能力、試驗室規(guī)模等因素的限制，難以進行大型結(jié)構(gòu)的抗震試驗,不過，由地震作用導(dǎo)致的建筑結(jié)構(gòu)破壞大多為局部結(jié)構(gòu)損壞造成的整體結(jié)構(gòu)失效[2-3]，故學(xué)術(shù)界和工程界常對相關(guān)結(jié)構(gòu)中最薄弱部位進行實際抗震性能試驗，而對結(jié)構(gòu)其它部位則借助數(shù)值模擬，以此實現(xiàn)對大型結(jié)構(gòu)的抗震性能研究并真實反映相應(yīng)結(jié)構(gòu)的地震響應(yīng)，如郭玉榮等[4]基于OpenSees有限元軟件自帶的通訊模塊，通過OpenSees-OpenSees之間的信息交互，完成了復(fù)雜子結(jié)構(gòu)的擬動力試驗，此外，該課題組[5]還利用Netslab平臺開展了遠程協(xié)同擬動力試驗，整合了各試驗室資源，驗證了該試驗平臺的可靠性；蔡新江等[6]利用CSIC和MATLAB將開源試驗平臺OpenFresco與MTS電液伺服加載系統(tǒng)無縫連接，并進行了一榀兩層單跨平面鋼框架的混合試驗驗證；許國山等[7]應(yīng)用OpenSees計算數(shù)值子結(jié)構(gòu)，并通過OpenFresco接口程序?qū)崿F(xiàn)有限元軟件與試驗控制LabVIEW之間的數(shù)據(jù)通信，應(yīng)用dSpace控制作動器加載試件，驗證了試驗系統(tǒng)的穩(wěn)定性。

鋼管混凝土框架組合結(jié)構(gòu)具有良好的受力性能，在實際工程中應(yīng)用廣泛，為進一步提升該類結(jié)構(gòu)的安全性和經(jīng)濟性，本課題組[8]選取T形截面鋼管混凝土異形柱-工字鋼梁框架頂層邊節(jié)點為研究對象，按1∶2的縮尺比例設(shè)計并制作了3個“弱節(jié)點”模型和1個“強節(jié)點”模型，并對其進行加載破壞試驗，分析了節(jié)點荷載特征值、延性、耗能及剛度退化，同時還基于方鋼管混凝土柱-鋼梁組合框架中柱節(jié)點和邊柱節(jié)點的抗震性能試驗，利用ABAQUS進行建模分析，計算結(jié)果與試驗結(jié)果基本吻合[9]。此外，本課題組還建立了鋼管混凝土柱-不等高鋼梁節(jié)點在低周往復(fù)荷載作用下的破壞性試驗?zāi)Ｐ蚚10]，且該模型擬靜力試驗結(jié)果與數(shù)值模擬結(jié)果吻合較好，當(dāng)模擬水平位移增大時，結(jié)構(gòu)出現(xiàn)明顯剛度退化現(xiàn)象，整體變化趨勢和試驗接近[11]，在該數(shù)值模型正確性得到驗證的基礎(chǔ)上，本文采用OpenSees平臺，通過其自帶的通訊模塊實現(xiàn)信息交互，以鋼管混凝土柱-不等高鋼梁框架底層中節(jié)點為試驗子結(jié)構(gòu)，其余部分為數(shù)值子結(jié)構(gòu)，進行子結(jié)構(gòu)擬靜力試驗?zāi)M以驗證子結(jié)構(gòu)方法的可靠性，并進一步分析節(jié)點左右側(cè)梁高比、柱截面尺寸以及軸壓比等參數(shù)對鋼管混凝土柱-不等高鋼梁框架結(jié)構(gòu)抗震性能的影響。

1 基于OpenSees的子結(jié)構(gòu)擬靜力試驗

1.1 模型的建立及子結(jié)構(gòu)的選取

基于本課題組所建且已通過正確性驗證的鋼管混凝土柱-不等高鋼梁節(jié)點在低周往復(fù)荷載作用下的破壞性試驗?zāi)Ｐ蚚10]，利用OpenSees程序構(gòu)建圖1所示模型。該模型結(jié)構(gòu)總高度為9.480 m，高梁跨度為4 m，低梁跨度為2 m，混凝土強度等級均為C40，鋼管柱和H型鋼梁均采用Q235B鋼材并考慮P-Delta效應(yīng)，中柱軸壓比取0.4，左右邊柱軸壓比分別取0.3、0.2，構(gòu)件鋼材力學(xué)性能如表1所示。從模型中一榀三層兩跨框架的梁柱反彎點位置截取底層中節(jié)點作為試驗子結(jié)構(gòu)，其余部分為數(shù)值子結(jié)構(gòu)，選用OpenSees有限元軟件分別建立數(shù)值模型?？紤]計算精度和計算效率，采用基于柔度法的纖維梁柱單元，每個單元選取4個積分點進行插值積分以獲得單元的剛度和抗力。在鋼管混凝土柱-不等高鋼梁組合結(jié)構(gòu)中，因混凝土在鋼管的包圍作用下形成約束區(qū)，核心混凝土強度也隨之提升，故本文選用OpenSees平臺自帶的concrete04混凝土模型及Steel02鋼材模型[11]。

(a)完整結(jié)構(gòu)示意圖 (b)完整結(jié)構(gòu)節(jié)點分布示意圖

(c)試驗子結(jié)構(gòu)示意圖

表1 鋼材力學(xué)性能

1.2 子結(jié)構(gòu)擬靜力試驗

1.2.1 子結(jié)構(gòu)擬靜力試驗方法原理

由數(shù)值子結(jié)構(gòu)和試驗子結(jié)構(gòu)共同組成的整體結(jié)構(gòu)在力學(xué)方程中得到統(tǒng)一，通過計算軟件施加結(jié)構(gòu)整體位移荷載，再利用作動器施加對應(yīng)試驗子結(jié)構(gòu)的位移。模擬試驗參照文獻[12]進行，具體步驟如下：①將結(jié)構(gòu)初始參數(shù)帶入靜力方程，利用有限元軟件求解，得到整體結(jié)構(gòu)的第一步位移；②從整體位移中提取試驗子結(jié)構(gòu)的位移，并利用作動器施加到試驗子結(jié)構(gòu)上，同時由傳感器測得恢復(fù)力；③將得到的恢復(fù)力帶入靜力方程，求解整體結(jié)構(gòu)的下一步位移；④利用作動器施加試驗子結(jié)構(gòu)的下一步位移，并測得恢復(fù)力，如此循環(huán)，完成擬靜力加載的全過程。在本研究中，采用兩個OpenSees軟件分別模擬數(shù)值子結(jié)構(gòu)及試驗子結(jié)構(gòu)，對數(shù)值子結(jié)構(gòu)施加豎向恒定荷載，隨后施加位移控制的水平荷載，加載制度如圖2所示。從每一步加載后的數(shù)值子結(jié)構(gòu)中提取對應(yīng)的節(jié)點位移，通過OpenSees數(shù)據(jù)交互通道將節(jié)點位移施加給試驗子結(jié)構(gòu)，并測得恢復(fù)力回傳給數(shù)值子結(jié)構(gòu)，進行下一步分析。為了實現(xiàn)上述循環(huán)過程，需要借助OpensSees里的if語言，即

if {set i 1} {$i<=$Nsteps} {

…

analyze 1

…

}

從而能夠在循環(huán)中，執(zhí)行每一個分析步內(nèi)的模型修改，將每一步的位移和恢復(fù)力分別作用在試驗子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)上。省略號部分代表具體的修改命令。

圖2 水平加載系統(tǒng)

1.2.2 基于OpenSees的網(wǎng)絡(luò)通信協(xié)議實現(xiàn)

為了實現(xiàn)數(shù)值子結(jié)構(gòu)和試驗子結(jié)構(gòu)的數(shù)據(jù)交互，需要在有限元軟件和試驗加載裝置間建立數(shù)據(jù)傳輸通道。因本研究利用兩個OpenSees軟件分別模擬數(shù)值子結(jié)構(gòu)和試驗子結(jié)構(gòu)，故需要利用網(wǎng)絡(luò)基于TCP/IP通信協(xié)議在兩個OpenSees窗口間實現(xiàn)數(shù)據(jù)交換。命令流如下文所示:

# 服務(wù)端

set count 0

proc Echo_Server {port} {

global count

set s [socket -server EchoAccept $port]

vwait count

}

proc EchoAccept {sock addr port} {

global echo

global count

global s

set s $sock

# Record the client’s information

puts "Accept $sock from $addr port $port"

set echo(addr,$sock) [list $addr $port]

# Ensure that each "puts" by the server

# results in a network transmission

fconfigure $sock -buffering line

# Set up a callback for when the client sends data

# add the event

fileevent $sock readable [list Echo $sock]

incr count 1

}

Echo_Server 9992

# 客戶端

proc Echo_Client {host port} {

set s [socket $host $port]

fconfigure $s -buffering line

return $s

}

set s [Echo_Client localhost 9992]

將兩個OpenSees窗口作為服務(wù)端和客戶端連接，并通過gets和puts命令收發(fā)數(shù)據(jù)，從而實現(xiàn)位移和恢復(fù)力的傳遞交互。

2 有限元結(jié)果分析

2.1 滯回曲線與骨架曲線

完整結(jié)構(gòu)框架和數(shù)值子結(jié)構(gòu)頂層滯回曲線、試驗子結(jié)構(gòu)節(jié)點滯回曲線以及各自的骨架曲線分別如圖3和圖4所示。結(jié)合圖3(a)模擬結(jié)果可見，完整結(jié)構(gòu)框架和數(shù)值子結(jié)構(gòu)的滯回曲線吻合較好，證明了本研究所用模擬試驗方法的準(zhǔn)確性；通過比較圖4(a)中兩者的骨架曲線可見，完整結(jié)構(gòu)的承載力(P)稍高于數(shù)值子結(jié)構(gòu)承載力，差值在8%以內(nèi)，這主要是因為在模擬試驗過程中，只控制了15節(jié)點的水平位移，并未考慮其豎向和轉(zhuǎn)角位移，邊界條件并不充分，加之初始條件存在誤差，子結(jié)構(gòu)試驗第一步是對數(shù)值子結(jié)構(gòu)進行位移加載，此時數(shù)值子結(jié)構(gòu)并未施加邊界條件，因而會對后續(xù)的試驗產(chǎn)生微小誤差；曲線初始剛度偏大，主要是在建模過程中并未考慮鋼管與混凝土之間的黏結(jié)滑移作用，同時曲線的強度退化并不明顯，尤其是在后期同一位移幅度下，承載力降低程度越來越小。此外，結(jié)合圖3(b)和圖4(b)可知，在框架位移加載前期，試驗子結(jié)構(gòu)節(jié)點基本處于彈性工作狀態(tài)。隨著位移增大，節(jié)點滯回環(huán)面積增大，滯回曲線基本呈梭形。當(dāng)位移進一步增大，正向滯回環(huán)面積持續(xù)增大，滯回曲線比較飽滿，未見明顯捏縮現(xiàn)象，表明節(jié)點具有良好的耗能能力和延性。與節(jié)點滯回曲線相似，整體框架的滯回環(huán)面積持續(xù)增大，基本呈梭形，具有較好的耗能能力。

(a)完整框架、數(shù)值子結(jié)構(gòu)頂層滯回曲線

(b)試驗子結(jié)構(gòu)節(jié)點滯回曲線

(a) 完整框架、數(shù)值子結(jié)構(gòu)骨架曲線

(b) 試驗子結(jié)構(gòu)節(jié)點骨架曲線

2.2 破壞形態(tài)

OpenSees雖然沒有應(yīng)力云圖，但是能夠提取纖維截面的彎矩與曲率，并給出具體加載步數(shù)。選擇節(jié)點12～13的梁單元以及節(jié)點1～4柱單元，提取框架梁端及柱腳截面彎矩-曲率如圖5所示。從圖5中可見，在荷載加載初期，框架處于彈性階段，彎矩與曲率成線性關(guān)系；當(dāng)加載至1630步時，框架高梁上翼緣開始屈服，此時框架柱處于彈性階段；繼續(xù)加載至7210步時，框架梁端彎矩到達屈服極限，框架柱仍處于彈性階段；繼續(xù)加載時，梁端截面抗彎剛度下降，隨著曲率增大，彎矩增加緩慢，當(dāng)加載至11 620步時，梁端出現(xiàn)塑性鉸，彎矩不再增大,此時柱腳開始屈服；當(dāng)加載至15 220步時，柱腳彎矩到達屈服極限，而后隨著繼續(xù)加載，柱腳截面抗彎剛度下降，彎矩增加緩慢，當(dāng)加載至21 430步時，柱腳出現(xiàn)塑性鉸，彎矩不再增大。

(a)框架梁端彎矩-曲率

(b)框架柱腳彎矩-曲率

3 參數(shù)分析

3.1 左右側(cè)梁高比

控制節(jié)點和框架對應(yīng)層的高梁高度h2為定值(560 mm)，改變低梁高度h1，在其它相關(guān)參數(shù)保持不變的條件下，分別進行了梁高比(ΔH=h1/h2)為0.3、0.4、0.6、0.8、1的子結(jié)構(gòu)抗震模擬試驗，所得框架及節(jié)點骨架曲線如圖6所示。由圖6可知，隨著節(jié)點梁高比的增大，框架承載力呈上升趨勢，其中當(dāng)梁高比從0.6增至0.8時，承載力最大增幅為4.11%，承載力特征值對應(yīng)的位移值并無明顯變化。而對于節(jié)點自身，當(dāng)梁高比增大時，節(jié)點承載力特征值提高，其對應(yīng)的位移值隨梁高比的增大逐漸減小，其中當(dāng)梁高比從0.6增至0.8時，節(jié)點極限承載力提高了10.68%，極限承載力對應(yīng)的位移減小了8.64%，破壞點對應(yīng)的位移減小了4.31%。模擬結(jié)果表明，框架內(nèi)單一節(jié)點梁高比的改變對框架和節(jié)點自身承載力的影響較為明顯，提高節(jié)點左右側(cè)梁高比可同時提高框架和節(jié)點承載力，但節(jié)點自身延性略有下降。

(a)框架骨架曲線

(b)節(jié)點骨架曲線

3.2 柱截面尺寸

改變節(jié)點和框架對應(yīng)層的柱截面尺寸，考慮到結(jié)構(gòu)數(shù)值模型計算的穩(wěn)定性，分別對柱邊長a為350、400、450、500 mm的子結(jié)構(gòu)模型進行抗震模擬試驗，所得框架及節(jié)點骨架曲線如圖7所示。由圖7可知，當(dāng)節(jié)點柱截面變大時，梁柱線剛度比(高梁截面線剛度與柱截面線剛度比值)逐漸減小，框架和節(jié)點的初始剛度與承載力隨之提升，當(dāng)柱截面邊長由350 mm增至400 mm時，梁柱線剛度比由1.53減至0.90，框架極限承載力最大增幅為5.99%，但承載力特征值對應(yīng)的位移并無明顯變化，節(jié)點極限承載力增幅為61.25%，極限承載力對應(yīng)的位移減小了13.40%，破壞點對應(yīng)的位移減小了5.08%；當(dāng)柱截面邊長不低于450 mm時，梁柱線剛度比小于0.56，框架與節(jié)點承載力提升并不顯著。模擬結(jié)果表明，因框架內(nèi)單一節(jié)點柱截面尺寸的改變而導(dǎo)致的梁柱線剛度比變化對框架承載力的影響較為明顯，柱截面的尺寸增大可同時提高框架與節(jié)點的承載力，雖然節(jié)點延性略有降低，但自身承載力的提升有助于增強結(jié)構(gòu)抗震能力。

(a)框架骨架曲線

(b)節(jié)點骨架曲線

3.3 軸壓比

改變節(jié)點和對應(yīng)框架軸壓比n，在其它相關(guān)參數(shù)保持不變的條件下，分別對節(jié)點軸壓比為0.2、0.4、0.6、0.8 的子結(jié)構(gòu)模型進行抗震模擬試驗，所得框架及節(jié)點骨架曲線如圖8所示。由圖8可知，隨著節(jié)點軸壓比增大，框架的承載力呈下降趨勢，當(dāng)軸壓比由0.2增至0.4時，框架承載力下降并不明顯，極限承載力下降2.89%；當(dāng)軸壓比由0.6增至0.8時，框架承載力骨架曲線形狀發(fā)生變化，極限承載力下降8.72%，破壞點對應(yīng)的位移急劇減小，延性大幅降低。而對于節(jié)點自身，其承載力變化特點與框架近似，其中當(dāng)軸壓比由0.6增至0.8時，節(jié)點極限承載力下降11.27%，承載力達到峰值點后迅速降低，剛度退化加快，延性大幅降低。模擬結(jié)果表明，軸壓比越大，框架和節(jié)點承載力降幅越明顯，剛度退化隨之加快，這將對承受低周往復(fù)荷載的鋼管混凝土柱-不等高鋼梁框架抗震性能產(chǎn)生不利影響。

(a)框架骨架曲線

(b)節(jié)點骨架曲線

4 結(jié)語

為了驗證基于OpenSees的子結(jié)構(gòu)擬靜力試驗的可靠性，建立了鋼管混凝土柱-不等高鋼梁框架與節(jié)點數(shù)值模型，以O(shè)penSees自帶的通訊模塊建立數(shù)據(jù)交互通道，進行有限元模擬，并重點研究了節(jié)點參數(shù)改變對鋼管混凝土柱-不等高鋼梁框架結(jié)構(gòu)抗震性能的影響。模擬計算結(jié)果表明，框架內(nèi)單一節(jié)點左右側(cè)梁高比的改變，對框架和節(jié)點承載力的影響較為明顯，當(dāng)節(jié)點左右側(cè)梁高比從0.6增至0.8時，框架極限承載力提高了4.11%，承載力特征值對應(yīng)的位移值并無明顯變化，節(jié)點極限承載力提高了10.68%，極限承載力對應(yīng)的位移值減小了8.64%，破壞點對應(yīng)的位移值減小了4.31%，節(jié)點在承載力提高的同時延性略有降低；因單一節(jié)點柱截面尺寸的改變而導(dǎo)致的梁柱線剛度比變化對框架和節(jié)點承載力均有明顯影響，當(dāng)柱截面邊長由350 mm增至400 mm時，梁柱線剛度比由1.53降至0.90，框架極限承載力提高了5.99%，承載力特征值對應(yīng)的位移值并無明顯變化，節(jié)點極限承載力提高了61.25%，極限承載力對應(yīng)的位移值減小了13.40%，破壞點對應(yīng)的位移值減小了5.08%，節(jié)點延性略有降低，當(dāng)柱截面邊長不低于450 mm時，梁柱線剛度比小于0.56，框架與節(jié)點承載力的提升均不明顯；軸壓比的存在使得框架和節(jié)點承載力下降，當(dāng)軸壓比由0.6增至0.8時，框架與節(jié)點承載力骨架曲線形狀發(fā)生變化，框架極限承載力下降了8.72%，節(jié)點極限承載力下降了11.27%，承載力到達峰值點后迅速降低，剛度退化加快，延性大幅下降。