郭玉榮 龍沐恩

摘? ?要：提出了一種利用鋼筋混凝土柱擬靜力試驗(yàn)數(shù)據(jù)識(shí)別改進(jìn)IMK模型骨架曲線參數(shù)，進(jìn)而提高鋼筋混凝土框架結(jié)構(gòu)非線性模擬精度的方法. 通過(guò)引入可抗差的基于奇異值分解的無(wú)跡卡爾曼濾波算法（抗差SVD-UKF算法），抑制觀測(cè)值粗差對(duì)參數(shù)識(shí)別的影響，采用粒子群算法對(duì)初始協(xié)方差矩陣、過(guò)程噪聲矩陣和測(cè)量噪聲矩陣進(jìn)行自動(dòng)尋優(yōu)，在MATLAB中實(shí)現(xiàn)了柱滯回特征正負(fù)向?qū)ΨQ與非對(duì)稱兩種情況下改進(jìn)IMK恢復(fù)力模型骨架曲線參數(shù)的識(shí)別. 鋼筋混凝土柱實(shí)測(cè)滯回曲線的模型骨架曲線參數(shù)識(shí)別結(jié)果及其在框架結(jié)構(gòu)非線性模擬中的應(yīng)用結(jié)果驗(yàn)證了本文方法的有效性.

關(guān)鍵詞：恢復(fù)力模型;滯回特征;參數(shù)識(shí)別;抗差SVD-UKF算法;粒子群算法

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

文章編號(hào)：1674—2974（2021）01—0126—09

Abstract：A method for identifying the backbone curve parameters of the modified Ibarra-Medina-Krawinkler （IMK） model by using quasi-static test data of reinforced concrete columns and thus improving the simulation accuracy of reinforced concrete frame structures is proposed in this paper. In this method， a robust unscented Kalman filtering algorithm based on singular value decomposition（robust SVD-UKF algorithm） is introduced to suppress the influence of gross error of the observation on the parameter identification， and the particle swarm optimization algorithm is adopted to automatically optimize the initial covariance matrix， the process and measurement noise matrices. The identification of backbone curve parameters of the modified IMK model is realized using MATLAB，in which the symmetric and asymmetric hysteresis behavior of the columns in the positive and negative direction is considered. The effectiveness of the proposed method is verified by model backbone curve parameter identification based on the measured hysteretic curves of reinforced concrete columns and its application in the nonlinear simulation of frame structures.

Key words：hysteretic model;hysteretic behavior;parameter identification;robust SVD-UKF algorithm;particle swarm optimization algorithm

鋼筋混凝土框架結(jié)構(gòu)的地震響應(yīng)混合模擬[1] 及其抗倒塌性能的分析，需要可有效模擬鋼筋混凝土構(gòu)件滯回特征[2]的恢復(fù)力模型及準(zhǔn)確的模型參數(shù). 塑性鉸模型是框架結(jié)構(gòu)非線性模擬常采用的一種模型，它不僅反映構(gòu)件的力學(xué)特征，還與構(gòu)件的材料、約束狀況及空間布局密切相關(guān). 幾十年以來(lái)，塑性鉸模型已有了飛速發(fā)展，Clough等[3]開(kāi)發(fā)了雙線性模型;Wen[4]提出了光滑的塑性鉸模型;Takeda等[5]開(kāi)發(fā)了三線性塑性鉸模型. 但是上述常見(jiàn)模型不能充分考慮構(gòu)件在循環(huán)往復(fù)荷載作用下的剛度和強(qiáng)度退化，影響了整體結(jié)構(gòu)模擬分析的精確性. Ibarra等[6-7]開(kāi)發(fā)了復(fù)雜的塑性鉸模型即改進(jìn)的IMK （Ibarra-Medina-Krawinkler）模型. 改進(jìn)IMK模型作為一種以三折線為骨架曲線的塑性鉸模型，引入了基于能量耗散的退化參數(shù)β，考慮了構(gòu)件在往復(fù)荷載作用下的多種剛度和強(qiáng)度退化，相對(duì)于其他塑性鉸模型，能夠較好地對(duì)鋼筋混凝土梁柱的滯回特征進(jìn)行有效模擬. 但根據(jù)模型經(jīng)驗(yàn)公式直接計(jì)算的改進(jìn)IMK模型骨架曲線參數(shù)可能會(huì)存在明顯的誤差，從而影響框架結(jié)構(gòu)地震響應(yīng)模擬精度或?qū)е陆Y(jié)構(gòu)抗倒塌性能分析產(chǎn)生較大誤差，因此，在選擇恢復(fù)力模型之后，模型參數(shù)的準(zhǔn)確度成為影響結(jié)構(gòu)體系模擬精度的主要因素之一. 為此，本文提出一種結(jié)合擬靜力試驗(yàn)的混合模擬方法，先從框架結(jié)構(gòu)中選取部分關(guān)鍵構(gòu)件進(jìn)行擬靜力試驗(yàn)，然后基于構(gòu)件實(shí)測(cè)滯回曲線對(duì)構(gòu)件恢復(fù)力模型骨架曲線參數(shù)進(jìn)行識(shí)別，最后將識(shí)別的參數(shù)用于更新結(jié)構(gòu)中相同構(gòu)件的恢復(fù)力模型參數(shù)，并進(jìn)行整體結(jié)構(gòu)非線性數(shù)值模擬. 在這一混合模擬方法中，模型參數(shù)識(shí)別是最關(guān)鍵的環(huán)節(jié)之一.

目前在土木工程領(lǐng)域廣泛使用的參數(shù)識(shí)別方法主要有最小二乘估計(jì)法（Least Square Estimation，LSE）、擴(kuò)展卡爾曼濾波法（Extended Kalman Filter，EKF）和無(wú)跡卡爾曼濾波法（Unscented Kalman Filter，UKF）[8]，其中UKF算法以其高效性成為非線性參數(shù)識(shí)別的常用方法. 但UKF方法在實(shí)際運(yùn)用時(shí)仍存在一些局限性，系統(tǒng)狀態(tài)先驗(yàn)信息矩陣會(huì)因?yàn)橛^測(cè)粗差、系統(tǒng)噪聲等不確定的因素失去正定性，中斷濾波進(jìn)程，觀測(cè)粗差等因素還會(huì)影響濾波過(guò)程，影響收斂速度和識(shí)別結(jié)果;此外，在UKF方法中初始協(xié)方差矩陣及過(guò)程噪聲和測(cè)量噪聲矩陣的確定較為繁瑣，但它們的取值卻又決定著參數(shù)識(shí)別的準(zhǔn)確性. 基于此，本文運(yùn)用可抗差的基于奇異值分解的無(wú)跡卡爾曼濾波算法（robust unscented Kalman filtering algorithm based on singular value decomposition，抗差SVD-UKF算法）[9]來(lái)解決狀態(tài)先驗(yàn)信息矩陣失去正定性的問(wèn)題并抑制粗差等對(duì)濾波過(guò)程的影響，并通過(guò)粒子群算法對(duì)初始協(xié)方差矩陣及過(guò)程噪聲矩陣和測(cè)量噪聲矩陣進(jìn)行自動(dòng)尋優(yōu)，采用MATLAB實(shí)現(xiàn)了改進(jìn)IMK恢復(fù)力模型骨架曲線參數(shù)的識(shí)別. 通過(guò)基于陸新征等的框架柱擬靜力試驗(yàn)數(shù)據(jù)[10]的模型參數(shù)識(shí)別以及對(duì)其整體框架結(jié)構(gòu)擬靜力試驗(yàn)結(jié)果[11]的模擬對(duì)比，來(lái)驗(yàn)證本文方法的有效性.

1? ?改進(jìn)的IMK恢復(fù)力模型

改進(jìn)的IMK恢復(fù)力模型是一種用塑性鉸來(lái)模擬梁柱構(gòu)件非線性行為的恢復(fù)力模型，它采用三折線骨架曲線（如圖1所示），其骨架曲線形狀由5個(gè)參數(shù)EIy、θcap，p l、θpc、My、Mc /My確定，參數(shù)的計(jì)算如式（1）～式（5）[12]所示，其循環(huán)能量耗散能力λ的計(jì)算如公式（6）所示.

式中：EIy為割線模量;θcap，p l為塑性轉(zhuǎn)角;θpc為峰值后轉(zhuǎn)角;My為屈服彎矩值;Mc /My為屈服后硬化剛度;P/Ag? f ′ c、ν均為軸壓比;？準(zhǔn)y和ky為由構(gòu)件基本信息計(jì)算出的參數(shù)，在此不作展開(kāi);δ′= d′/d，d′為受壓區(qū)邊緣到受壓鋼筋中心的距離;s為塑性鉸區(qū)的箍筋間距;asl為縱向鋼筋滑移系數(shù)（考慮取1，不考慮取0）;sn為鋼筋屈曲系數(shù)，sn=（s/db）（fy /100）0.5，db為縱筋直徑，fy為縱筋屈服強(qiáng)度;f ′ c為混凝土軸心抗壓強(qiáng)度，單位為MPa;cunits是單位換算量，當(dāng)f ′ c的單位為MPa時(shí)cunits取1;ρsh為柱塑性鉸區(qū)的橫向鋼筋面積比;ρ為縱筋配筋率;b和d分別為柱橫截面的寬度和高度.

2? ?無(wú)跡卡爾曼濾波算法

2.1? ?標(biāo)準(zhǔn)UKF算法

2.1.1? ?時(shí)間預(yù)測(cè)

利用第K-1步狀態(tài)向量構(gòu)造一組總數(shù)為2L+1的sigma點(diǎn)集，L為狀態(tài)量的維度，X在本文中為待識(shí)別的參數(shù) K1、K2、K3、 fy、 fp和恢復(fù)力r組成的狀態(tài)向量：

2.2? ?抗差SVD-UKF算法

在觀測(cè)粗差、計(jì)算機(jī)截?cái)嗾`差等因素的干擾下，狀態(tài)協(xié)方差矩陣可能會(huì)失去正定性，導(dǎo)致無(wú)損變換（Unscented Transformation，UT）中的Cholesky分解失效. 本文采用SVD分解生成sigma點(diǎn)集，之后在UKF算法中通過(guò)抗差因子調(diào)節(jié)噪聲矩陣來(lái)抑制觀測(cè)粗差等因素對(duì)濾波過(guò)程的影響，該算法的流程圖如圖2所示.

2.2.1? ?抗差模型

抗差因子將觀測(cè)值作為考量因素引入濾波過(guò)程中，起到動(dòng)態(tài)調(diào)整測(cè)量噪聲的作用. 采用觀測(cè)殘差構(gòu)建抗差因子，抗差因子θk的表達(dá)式與IGGIII函數(shù)的形式相似，其表達(dá)式如下：

式中：k0、k1為閾值參數(shù)，k0通常取1.0 ～ 2.0，k1取3.0 ～8.5，本文中k0取1，k1取3;vk為當(dāng)前步觀測(cè)值與估計(jì)值的殘差; vk為當(dāng)前步觀測(cè)值與估計(jì)值殘差的均值;svk、σk分別為當(dāng)前步觀測(cè)值與估計(jì)值殘差的標(biāo)準(zhǔn)化量和觀測(cè)值與估計(jì)值殘差的樣本標(biāo)準(zhǔn)差.

2.2.2? ?SVD分解與觀測(cè)協(xié)方差矩陣的調(diào)整

為了解決在異常情況下?tīng)顟B(tài)協(xié)方差矩陣失去正定性而無(wú)法使用Cholesky分解的問(wèn)題，引入SVD分解計(jì)算sigma點(diǎn)集，計(jì)算公式如下：

2.2.3? ?粒子群算法的自動(dòng)尋優(yōu)

粒子群算法是模擬鳥群隨機(jī)搜尋食物的捕食行為，鳥群通過(guò)自身經(jīng)驗(yàn)和種群之間的交流調(diào)整自己的搜尋路徑，從而找到食物最多的地點(diǎn)，其中每只鳥的位置為自變量組合，每次到達(dá)地點(diǎn)的食物密度即為函數(shù)值. 每次搜尋都會(huì)根據(jù)自身經(jīng)驗(yàn)（自身歷史搜尋的最優(yōu)地點(diǎn)）和種群交流（種群歷史搜尋的最優(yōu)地點(diǎn)）調(diào)整自身的搜尋方向和速度，這稱為跟蹤極值，從而找到最優(yōu)解. 粒子群算法的流程圖如圖3所示.

本文中定義粒子的位置為初始協(xié)方差矩陣的每個(gè)元素以及過(guò)程噪聲和測(cè)量噪聲矩陣的每個(gè)元素，即相當(dāng)于構(gòu)造了一個(gè)N維的空間，并讓鳥群在這個(gè)N維空間中搜尋使適應(yīng)度函數(shù)取值最小的坐標(biāo)解，適應(yīng)度函數(shù)為在每一個(gè)粒子自變量組合經(jīng)過(guò)SVD-UKF算法識(shí)別參數(shù)后，由識(shí)別參數(shù)計(jì)算所得的每個(gè)恢復(fù)力值與相應(yīng)試驗(yàn)值的殘差平方和，殘差平方和越大表示識(shí)別的結(jié)果越不可靠，反之則識(shí)別結(jié)果越可靠，也即粒子自身包含的自變量組合的取值越合理. 狀態(tài)轉(zhuǎn)移方程是表示粒子由當(dāng)前時(shí)刻至下一時(shí)刻的位置變化方式，狀態(tài)轉(zhuǎn)移方程如下：

2.2.4? ?抗差SVD-UKF算法抗差效果分析

基于改進(jìn)IMK三折線模型的MATLAB仿真，以表1參數(shù)為仿真參數(shù)值，通過(guò)輸入位移峰值依次為4、10 、15 、20 、25 、30 、40 、45 、55 mm以及各級(jí)位移峰值循環(huán)兩圈的加載序列，計(jì)算得到恢復(fù)力序列，之后為檢驗(yàn)抗差效果將兩類異常觀測(cè)引入仿真結(jié)果.

1）單點(diǎn)異常：N（0，2Z（i）2）（均值為0，方差為2Z（i）2）的單點(diǎn)隨機(jī)誤差，其中Z（i）為當(dāng)前恢復(fù)力值.

2）連續(xù)段異常：N（0，2Z（i）2）的連續(xù)段隨機(jī)誤差.

仿真結(jié)果加入異常后的觀測(cè)值如圖4所示，之后分別采用未抗差的UKF算法及抗差的SVD-UKF算法對(duì)加入異常的觀測(cè)恢復(fù)力值進(jìn)行濾波，濾波過(guò)程如圖5所示.

圖6為在引入單點(diǎn)異常后截取的2 000～4 000加載步的濾波過(guò)程圖. 由圖6可見(jiàn)，抗差算法的濾波過(guò)程跟真實(shí)值基本吻合，而未抗差算法的濾波過(guò)程受到單點(diǎn)異常影響而出現(xiàn)了偏離.

誤差后截取的5 000～7 000加載步的濾波過(guò)程圖. 由圖7可知，抗差算法的濾波過(guò)程跟真實(shí)值基本吻合，而未抗差算法的濾波過(guò)程受到了連續(xù)段異常的影響而出現(xiàn)了較大波動(dòng).

可見(jiàn)，在兩種異常的影響下抗差SVD-UKF算法擁有較強(qiáng)的魯棒性與吻合程度.

3? ?低周反復(fù)荷載下鋼筋混凝土柱的骨架曲線參數(shù)識(shí)別

3.1? ?試驗(yàn)概況

本文試驗(yàn)實(shí)測(cè)數(shù)據(jù)采自陸新征等[10-11]鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究，試驗(yàn)包括一榀3層3跨整體框架的擬靜力試驗(yàn)和一些關(guān)鍵構(gòu)件的擬靜力試驗(yàn). 關(guān)鍵構(gòu)件包括底層邊柱和中柱各兩個(gè)試件，層一的邊節(jié)點(diǎn)和中節(jié)點(diǎn)各一個(gè)試件，本節(jié)參數(shù)識(shí)別選用了其中邊柱A和中柱C的試驗(yàn)數(shù)據(jù). 邊柱A和中柱C的幾何尺寸和配筋如圖8所示，其中中柱角筋直徑為10 mm，鋼筋材料參數(shù)如表2所示，混凝土強(qiáng)度和對(duì)柱子施加的軸力如表3所示. 柱子擬靜力試驗(yàn)的加載路徑如下.

邊柱A試驗(yàn)中首先施加柱端豎向荷載至試驗(yàn)軸壓比，然后以水平荷載控制加載10 kN/級(jí)至30 kN，各級(jí)荷載循環(huán)一圈. 此后以水平加載點(diǎn)位移控制加載，所加水平位移峰值分別為10、15、20、25、30、37.5、45、55 mm，各級(jí)位移循環(huán)兩圈.

中柱C試驗(yàn)中首先施加柱端豎向荷載至試驗(yàn)軸壓比，之后以水平荷載控制加載10 kN/級(jí)至40 kN，各級(jí)荷載循環(huán)一圈. 此后以水平加載點(diǎn)位移控制加載，所加水平位移峰值依次為10、15、20、25、30、37.5、45、55 mm，各級(jí)位移循環(huán)兩圈.

3.2? ?基于柱試驗(yàn)結(jié)果的改進(jìn)IMK對(duì)稱模型參數(shù)識(shí)別

本文采用抗差SVD-UKF算法，在MATLAB中編寫了識(shí)別改進(jìn)IMK對(duì)稱模型骨架曲線參數(shù)的程序，識(shí)別參數(shù)包括彈性剛度K1、硬化剛度K2、退化剛度K3、屈服力fy和極限力fp這5個(gè)參數(shù). 在后續(xù)的滯回曲線對(duì)比中，試驗(yàn)結(jié)果為柱子擬靜力試驗(yàn)的實(shí)測(cè)滯回曲線;初始結(jié)果是指通過(guò)公式（1）～（6）計(jì)算得出改進(jìn)IMK模型參數(shù)值后，依據(jù)柱子的位移加載路徑計(jì)算得到的滯回曲線;以初始結(jié)果中計(jì)算得到的模型參數(shù)為參數(shù)識(shí)別的初始值，識(shí)別出5個(gè)待識(shí)別參數(shù)后，依據(jù)柱子的位移加載路徑計(jì)算得到的滯回曲線稱為識(shí)別結(jié)果.

邊柱A的改進(jìn)IMK模型骨架曲線參數(shù)的初始值和識(shí)別值對(duì)比如表4所示. 用初始值和識(shí)別值分別計(jì)算出柱子的滯回曲線，邊柱A計(jì)算滯回曲線和試驗(yàn)滯回曲線的對(duì)比分別如圖9和圖10所示，與初始結(jié)果相比，識(shí)別結(jié)果與試驗(yàn)結(jié)果吻合較好，也即通過(guò)參數(shù)識(shí)別提高了模型參數(shù)的準(zhǔn)確度.

圖11為邊柱A的試驗(yàn)結(jié)果、識(shí)別結(jié)果和初始結(jié)果的滯回環(huán)面積對(duì)比. 由圖11的滯回環(huán)面積數(shù)據(jù)計(jì)算得到邊柱A的識(shí)別結(jié)果滯回環(huán)總面積與試驗(yàn)結(jié)果滯回環(huán)總面積相差8.15%，邊柱A初始結(jié)果滯回環(huán)總面積與試驗(yàn)結(jié)果滯回環(huán)總面積相差15.77%. 從滯回環(huán)總面積差來(lái)看，與初始參數(shù)相比，采用識(shí)別參數(shù)可以更好地模擬柱構(gòu)件的滯回曲線.

3.3? ?基于柱試驗(yàn)結(jié)果的改進(jìn)IMK非對(duì)稱模型參數(shù)識(shí)別

在上述識(shí)別方法中，通過(guò)修改MATLAB中的程序，實(shí)現(xiàn)了改進(jìn)IMK非對(duì)稱模型骨架曲線參數(shù)的識(shí)別.

中柱C的改進(jìn)IMK模型骨架曲線參數(shù)初始值和識(shí)別值對(duì)比如表5所示. 用初始值和識(shí)別值分別計(jì)算出柱子的滯回曲線，計(jì)算滯回曲線和試驗(yàn)滯回曲線的對(duì)比分別如圖12和圖13所示. 與初始結(jié)果相比，識(shí)別結(jié)果與試驗(yàn)結(jié)果吻合較好，也即通過(guò)參數(shù)識(shí)別提高了模型參數(shù)的準(zhǔn)確度.

圖14為中柱C試驗(yàn)結(jié)果、識(shí)別結(jié)果和初始結(jié)果的滯回環(huán)面積對(duì)比. 由圖14的滯回環(huán)面積數(shù)據(jù)計(jì)算得到識(shí)別結(jié)果滯回環(huán)總面積與試驗(yàn)結(jié)果滯回環(huán)總面積相差26.03%，初始結(jié)果滯回環(huán)總面積與試驗(yàn)結(jié)果滯回環(huán)總面積相差31.79%. 可見(jiàn)，與初始參數(shù)相比，采用識(shí)別參數(shù)可更好地模擬柱構(gòu)件的滯回曲線.

4? ?基于柱構(gòu)件參數(shù)識(shí)別結(jié)果的應(yīng)用

為檢驗(yàn)基于柱實(shí)測(cè)滯回曲線的模型骨架曲線參數(shù)識(shí)別值的有效性，以下將識(shí)別的模型參數(shù)應(yīng)用于整體框架結(jié)構(gòu)的數(shù)值模擬，所選的框架結(jié)構(gòu)來(lái)自陸新征等[10-11]完成的一榀3層3跨鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn). 整體框架的數(shù)值模擬采用OpenSees建模，建模流程如下.

首先計(jì)算得到改進(jìn)IMK模型的輸入?yún)?shù)，然后在OpenSees中調(diào)用改進(jìn)IMK模型并定義塑性鉸材料，最后將材料用于梁柱端部零長(zhǎng)度單元的彎曲特性，其余單元均為彈性梁柱單元. 柱端部改進(jìn)IMK模型參數(shù)按兩種工況取值換算，工況1采用由柱端截面的尺寸、材料、配筋、軸力等參數(shù)按照公式（1）～（6）計(jì)算得到的模型參數(shù)值;工況2采用基于柱實(shí)測(cè)試驗(yàn)數(shù)據(jù)得到的參數(shù)識(shí)別值，其中中柱的參數(shù)值選用中柱C負(fù)向的參數(shù)識(shí)別值，將參數(shù)識(shí)別值換算成模型參數(shù)值后，對(duì)割線剛度EIy、塑性轉(zhuǎn)角θcap，p l和峰值后轉(zhuǎn)角θpc根據(jù)框架的實(shí)際軸力按照公式（1）（3）（4）中軸壓比的分項(xiàng)作了調(diào)整，考慮到屈服彎矩與軸壓比關(guān)系的復(fù)雜性，對(duì)屈服彎矩My和屈服后硬化剛度Mc /My分別按公式（2）和公式（5）計(jì)算取值，對(duì)λ取一較大值，即不考慮構(gòu)件的循環(huán)能量耗散. 之后，根據(jù)擬靜力試驗(yàn)設(shè)定的加載路徑，進(jìn)行了兩種工況下框架結(jié)構(gòu)的數(shù)值模擬. 圖15為兩種工況下基底剪力-頂點(diǎn)位移數(shù)值模擬滯回曲線和試驗(yàn)滯回曲線的對(duì)比. 圖16為工況1、工況2與試驗(yàn)結(jié)果的滯回環(huán)面積對(duì)比.

從圖15曲線對(duì)比可以看出，工況1與工況2均對(duì)試驗(yàn)情況作出了較好的模擬. 由圖16的滯回環(huán)面積數(shù)據(jù)計(jì)算得到基底剪力-頂點(diǎn)位移識(shí)別結(jié)果的滯回環(huán)總面積與試驗(yàn)結(jié)果的滯回環(huán)總面積相差23.71%，初始結(jié)果滯回環(huán)總面積與試驗(yàn)結(jié)果滯回環(huán)總面積相差36.50%. 從滯回環(huán)總面積差來(lái)看，與初始參數(shù)相比，采用識(shí)別參數(shù)可以更好地模擬框架結(jié)構(gòu)的整體響應(yīng).

5? ?結(jié)? ?論

1）將抗差SVD-UKF算法運(yùn)用于鋼筋混凝土柱改進(jìn)IMK模型的參數(shù)識(shí)別，有效抑制了觀測(cè)粗差對(duì)識(shí)別進(jìn)程的影響，并通過(guò)粒子群算法對(duì)初始協(xié)方差矩陣、過(guò)程噪聲矩陣、測(cè)量噪聲矩陣進(jìn)行自動(dòng)尋優(yōu)，省去了大量的人力調(diào)參工作.

2）基于陸新征等[10-11]鋼筋混凝土框架柱構(gòu)件擬靜力試驗(yàn)實(shí)測(cè)滯回曲線數(shù)據(jù)，對(duì)柱子的改進(jìn)IMK模型骨架曲線參數(shù)進(jìn)行了識(shí)別，結(jié)果表明識(shí)別參數(shù)精度較初始參數(shù)精度有明顯提高.

3）將基于柱子試驗(yàn)數(shù)據(jù)識(shí)別的骨架曲線參數(shù)用于框架結(jié)構(gòu)的擬靜力加載試驗(yàn)數(shù)值模擬，結(jié)果表明利用識(shí)別參數(shù)可提高整體結(jié)構(gòu)的數(shù)值模擬精度.

參考文獻(xiàn)

[1]? ? 郭玉榮，劉鐘真. 基于OpenSees的子結(jié)構(gòu)擬動(dòng)力試驗(yàn)方法研究[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，42（9）：1—7.

GUO Y R，LIU Z Z. Research on substructure pseudo-dynamic test method based on OpenSees[J]. Journal of Hunan University （Natural Sciences），2015，42（9）：1—7. （In Chinese）

[2]? ? 過(guò)鎮(zhèn)海. 鋼筋混凝土原理[M]. 北京：清華大學(xué)出版社，1999：147—178，192—212.

GUO Z H. Reinforcement concrete theory[M]. Beijing： Tsinghua University Press，1999：147—178，192—212. （In Chinese）

[3]? ? CLOUGH R W，JOHNSTON S B. Effect of stiffness degradation on earthquake ductility requirements[C]//Proceedings of the 2nd Japan earthquake engineering symposium. Tokyo，1966.

[4]? ? WEN Y K. Method for random vibration of hysteretic systems[J]. Journal of the Engineering Mechanics Division，1976，102（2）：249—263.

[5]? ? TAKEDA T，SOZEN M A，NIELSEN N N. Reinforced concrete response to simulated earthquakes [J]. Journal of the Structural Division，1970，96（12）：2557—2573.

[6]? ? IBARRA L F，KRAWINKLER H. Global collapse of frame structures under seismic excitations [R]. California：Pacific Earthquake Engineering Research Center，2005.

[7]? ? IBARRA L F，MEDINA R A，KRAWINKLER H. Hysteretic models that incorporate strength and stiffness deterioration [J]. Earthquake Engineering & Structural Dynamics，2005，34（12）：1489—1511.

[8]? ? 王濤，吳斌. 基于約束UKF模型更新的混合試驗(yàn)方法[J]. 地震工程與工程振動(dòng)，2013，33（5）：100—109.

WANG T，WU B.Hybrid testing method based on model updating with constrained unscented Kalman filter[J]. Journal of Earthquake Engineering and Engineering Vibration，2013，33（5）：100—109. （In Chinese）

[9]? ? 譚興龍，王堅(jiān)，李增科. 基于SVD的改進(jìn)抗差UKF算法及在組合導(dǎo)航中的應(yīng)用[J]. 控制與決策，2014，29（10）：1744—1750.

TAN X L，WANG J，LI Z K. SVD aided improved robust UKF algorithm and its application for integration navigation [J]. Control and Decision，2014，29（10）：1744—1750.（In Chinese）

[10]? 陸新征，葉列平，潘鵬，等. 鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究及數(shù)值模擬競(jìng)賽Ⅱ：關(guān)鍵構(gòu)件試驗(yàn)[J]. 建筑結(jié)構(gòu)，2012， 42（11）：23—26.

LU X Z，YE L P，PAN P，et al. Pseudo-static collapse experiments and numerical prediction competition of RC frame structure Ⅱ：key elements experiment[J]. Building Structure，2012，42（11）：23—26. （In Chinese）

[11]? 陸新征，葉列平，潘鵬，等. 鋼筋混凝土框架結(jié)構(gòu)擬靜力倒塌試驗(yàn)研究及數(shù)值模擬競(jìng)賽Ⅰ：框架試驗(yàn)[J]. 建筑結(jié)構(gòu)，2012，42（11）：19—22.

LU X Z，YE L P，PAN P，et al. Pseudo-static collapse experiments and numerical prediction competition of RC frame structureⅠ：RC frame experiment[J]. Building Structure，2012，42（11）：19—22. （In Chinese）

[12]? HASELTON C B，LIEL S，TAYLOR L，et al. Beam-column element model calibrated for predicting flexural response leading to global collapse of RC frame buildings[R]. California：Pacific Earthquake Engineering Research Center，2007.