林煌斌,王全鳳
(1.華僑大學 土木工程學院,福建 泉州361021;2.同濟大學 航空航天與力學學院,上海200092;3.集美大學 工程技術(shù)學院,福建 廈門 361021)
描述結(jié)構(gòu)動力卸載剛度的指標驗證
林煌斌1,2,3,王全鳳1
(1.華僑大學 土木工程學院,福建 泉州361021;2.同濟大學 航空航天與力學學院,上海200092;3.集美大學 工程技術(shù)學院,福建 廈門 361021)
提出可用于描述結(jié)構(gòu)動力卸載剛度修正系數(shù)和能量變換關(guān)系系數(shù),建立卸載過程的彈性變形能及卸載剛度表達式,并用鋼骨混凝土柱恢復力試驗數(shù)據(jù)對結(jié)論進行可行性驗證 .研究結(jié)果表明:對于在評價結(jié)構(gòu)彈塑性時程分析中引入結(jié)構(gòu)卸載剛度修正系數(shù)和能量變換關(guān)系系數(shù),可以定量評價結(jié)構(gòu)動力卸載,較好地解釋卸載時,結(jié)構(gòu)卸載剛度矩陣隨時間變化構(gòu)成的時變參數(shù)系統(tǒng);彌補結(jié)構(gòu)彈塑性時程分析中未考慮的卸載剛度與運動狀態(tài)的相互影響效應,解決結(jié)構(gòu)切線剛度和割線剛度表示的動力方程不等效的問題.
描述指標;卸載剛度;混凝土結(jié)構(gòu);恢復力模型;彈性變形能
剛度退化模型的剛度變化不連續(xù),存在著拐點或突變點,但剛度計算比較簡單,工程應用上常用剛度退化二線型模型、剛度退化三線型模型、剛度退化四線型模型等.剛度退化三線型(DT3)模型是用三折線表示正反向加載恢復力骨架曲線,并考慮鋼筋混凝土結(jié)構(gòu)或構(gòu)件剛度退化性質(zhì),可更細致地描述鋼筋混凝土構(gòu)件的真實恢復力曲線.根據(jù)是否考慮結(jié)構(gòu)或構(gòu)件屈服后的硬化狀況,剛度退化三線型模型分為坡頂退化三線型模型和平頂退化三線型模型[1-3].結(jié)構(gòu)系統(tǒng)工作時,經(jīng)常是在加載過程中伴隨著卸載的發(fā)生,除了加載對系統(tǒng)產(chǎn)生作用效應,也經(jīng)常伴隨著卸載產(chǎn)生的作用效應.結(jié)構(gòu)系統(tǒng)的加載作用與卸載作用往往相互影響,使結(jié)構(gòu)系統(tǒng)的變形和運動呈現(xiàn)出復雜的非線性[4-5].結(jié)構(gòu)振動中卸載時,結(jié)構(gòu)卸載剛度矩陣隨時間變化構(gòu)成時變參數(shù)系統(tǒng).目前,結(jié)構(gòu)彈塑性時程分析一般未考慮卸載剛度與運動狀態(tài)相互影響效應,往往造成結(jié)構(gòu)切線剛度表示的動力方程和結(jié)構(gòu)割線剛度表示的動力方程不等效,不能完整描述彈塑性結(jié)構(gòu)振動的動力卸載過程.而且,相關(guān)文獻對結(jié)構(gòu)恢復力模型卸載剛度的描述缺少量化指標 .基于此,本文推導恢復力模型卸載剛度參數(shù)的確定方法并用試驗給予驗證.
在增量法中,單元剛度矩陣中的EI0為彈性段內(nèi)的截面抗彎剛度,EI1和EI2為塑性鉸區(qū)的截面抗彎剛度 .該剛度為切線剛度,可以采用條分法精確計算.截面骨架曲線確定,實際計算采用少參數(shù)近似描述骨架曲線.由于內(nèi)力重分布,截面在加卸載時剛度變化,一般情況下取各轉(zhuǎn)化階段彎曲剛度,即
式(1),(2)中:Ku=2.5 Mu/φu.
確定混凝土恢復力模型的卸載剛度采用以上方法,考慮其隨著構(gòu)件位移幅值增加而退化的因素,工程中常采用試驗數(shù)據(jù)考慮軸壓比、配筋率等構(gòu)件物理特征進行擬合卸載剛度.
常用卸載剛度計算方法采用“卸載按彈性”的規(guī)則.在滯回過程中假定計算卸載剛度不變,實際上該過程中Ku(x)矩陣是在不斷改變.變形和運動狀態(tài)在匹配位置附近微小幅度波動時,可以得到穩(wěn)定的計算結(jié)果 .然而,由于計算結(jié)構(gòu)狀態(tài)偏移,必然產(chǎn)生計算誤差 .若變形和運動狀態(tài)偏移匹配位置,則相當于計算中包含“計算擾動”量,偏移程度越大,相應的“計算擾動”也越大;若偏移程度較大或總是往某一方向偏移,則最終可能因誤差積累導致出現(xiàn)計算不穩(wěn)定.
構(gòu)件滯回中,彈性階段變形不會引起結(jié)構(gòu)損傷,體系能量保持守恒,不影響結(jié)構(gòu)的滯回曲線.強化段結(jié)構(gòu)損傷后出現(xiàn)剛度和強度退化損傷,結(jié)構(gòu)對系統(tǒng)外釋放能量.由于結(jié)構(gòu)發(fā)生彈塑性變形后對系統(tǒng)外釋放能量,外荷載做功不等于結(jié)構(gòu)變形能,同時結(jié)構(gòu)損傷導致剛度退化,形函數(shù)矩陣R保持不變,本構(gòu)矩陣K相應變化.結(jié)構(gòu)卸載剛度退化與結(jié)構(gòu)系統(tǒng)能量變化有關(guān).
文中采用的是工程應用較多的平頂退化三線型模型[1-3],如圖1所示.
圖1 剛度平頂退化三線型模型Fig.1 Trilinear model of stiffness degradation
(1)當滯回運動處在上升階段(圖2a),假設(shè)在CY段的X 位置卸載,OC段剛度為Ke,CY段剛度為Kp(Kp=αKe,常數(shù)α由DT3模型確定).
(2)當滯回運動處在下降階段(圖2b),假設(shè)在YM 段的X位置卸載,OC段剛度Ke,YM 段剛度Kd(Kd=βKe,常數(shù)β由DT3模型確定).
圖2 卸載過程彈性變形能彈性變形能示意圖Fig.2 Elastic deformation energy during unloading
當結(jié)構(gòu)處于X位置,卸載過程上升段、下降段的結(jié)構(gòu)累積彈性變形能分別為
式(3),(4)中:∑SR,∑SD分別為上升段、下降段的結(jié)構(gòu)累積彈性變形能;Se為彈性階段變形能,Se=ΔOCXc;Sp為彈塑階段的彈性變形能增量,Sp=ΔCYX1;Sd為下降段的彈性變形能增量,Sd=ΔYXX2;Su為卸載釋放的彈性變形能,Su=ΔXXwXx;γ為能量變換相關(guān)系數(shù),是Se,Sp,Sd的關(guān)系參數(shù),與材料和構(gòu)件物理特征包括軸壓比n0、剪跨比l/d、縱筋配筋率ρl、箍筋配筋率ρw等特征相關(guān).
由式(4)可得能量變換相關(guān)系數(shù)的表達式為
卸載剛度退化應用于表征損傷指數(shù)上,與結(jié)構(gòu)損傷的演化現(xiàn)象相一致.構(gòu)件在循環(huán)受力過程中,應變ε可分解為可逆應變εe和不可逆應變ε′,而ε′包括塑性應變εp和不可閉合損傷殘余應變εu.由損傷原理可得
式(6)中:構(gòu)件破壞時定義為D=1,經(jīng)典損傷模型D采用變形-耗能雙參數(shù)組合[6];ν為僅考慮由最大位移破壞因素的組合系數(shù),0.5≤ν<1.0.
由熱力學內(nèi)變量理論,可得含唯象損傷演化構(gòu)件的力-位移比擬式為
由式(7)可知,采用卸載剛度Ku(x)描述結(jié)構(gòu)損傷值D理論是可行的.
對于延性構(gòu)件(2.0<μ<4.0),設(shè)損傷函數(shù) D=fd(x1,x2,…,xn)為單調(diào)遞增函數(shù),則有
卸載剛度的修正系數(shù)ξ的表達式為
式(9)中:ηi,ηi+1分別為第i次和第i+1次卸載剛度Kui與彈性段剛度比值 .
為驗證結(jié)構(gòu)卸載剛度修正系數(shù)ξ和能量變換關(guān)系系數(shù)γ的指標式(5),(9),采用勁性RC柱和高強混凝土柱變形性能及恢復力特性試驗[6-8]的實測數(shù)據(jù)進行驗證.能量變換關(guān)系系數(shù)γ、卸載剛度Ku和卸載剛度修正系數(shù)ξ與循環(huán)次數(shù)n的關(guān)系,分別如圖3~5所示.
從圖3可知,在未破壞前的循環(huán)段內(nèi),γ在不同構(gòu)件中的變化趨勢相同,表明彈性階段能量與卸載彈性能之間存在關(guān)系,可以通過γ(m0,l/d,ρ1,ρw)擬合.
從圖4可知,各構(gòu)件在循環(huán)過程中,Ku(x)的變化量呈線性單調(diào)變化,不同特征的構(gòu)件卸載剛度Ku在不同循環(huán)中數(shù)值并不保持彈性原則.圖3描述在滯回過程中的卸載剛度并不是保持不變的事實,但最為重要的是dKu(x)基本相等,可以用式(8)表達Kui+1和Kui(x)的關(guān)系.
圖3 能量變換關(guān)系系數(shù)與循環(huán)次數(shù)的關(guān)系Fig.3 Relationship between energy correlation coefficient and load cyclic number
圖4 卸載剛度與循環(huán)次數(shù)的關(guān)系Fig.4 Relationship between unloading stiffness and load cyclic number
圖5 卸載剛度修正系數(shù)與循環(huán)次數(shù)的關(guān)系 Fig.5 Relationship between unloading stiffness correction factor and load cyclic number
從圖5可知,計算得到ξ值穩(wěn)定保持在0.85≤ξ≤0.95區(qū)域內(nèi),可以對試驗數(shù)據(jù)進行擬合得到ξ.
以上試驗充分體現(xiàn)卸載剛度隨著構(gòu)件位移幅值增加而退化現(xiàn)象,可以解釋隨著實際構(gòu)件偏移程度較大或總是往某一方向偏移,最終可能因誤差積累導致出現(xiàn)計算不穩(wěn)定的現(xiàn)象.
引入結(jié)構(gòu)卸載剛度修正系數(shù)ξ和能量變換關(guān)系系數(shù)γ,可以更好地描述和評價結(jié)構(gòu)的動力卸載現(xiàn)象,彌補了結(jié)構(gòu)彈塑性時程分析中一般未考慮卸載剛度與運動狀態(tài)相互影響所帶來的不足,該結(jié)論可供結(jié)構(gòu)彈塑性時程計算參考.由卸載彈性釋放能量、彈性變形能增量和彈塑性階段彈性變形能增量關(guān)系確定卸載剛度關(guān)系,可以評價鋼筋混凝土結(jié)構(gòu)的動力卸載,建立以剛度修正和能量變換關(guān)系的新方法,更接近描述結(jié)構(gòu)構(gòu)件在彈塑性時程的規(guī)律.動態(tài)調(diào)整卸載剛度是簡單又實際操作意義的.
[1]江見鯨,陸新征,葉列平.混凝土結(jié)構(gòu)有限元分析[M].北京:清華大學出版社,2005.
[2]張新培.鋼筋混凝土抗震結(jié)構(gòu)非線性分析[M].北京:科學出版社,2003.
[3]呂西林,金國芳,吳曉涵.鋼筋混凝土結(jié)構(gòu)非線性有限元理論與應用[M].上海:同濟大學出版社,1997.
[4]黃慶豐,王全鳳.彈塑性結(jié)構(gòu)動力系統(tǒng)的時變剛度激勵效應分析[J].土木工程學報,2006,39(12):18-23.
[5]王全鳳,黃慶豐.損傷結(jié)構(gòu)彈塑性抗震性能分析及評價[J].華僑大學學報:自然科學版,2009,30(1):6-11.
[6]PARK Y J,ANG-ALFREDO H S.Mechanistic seismic damage model for reinforced concrete[J].Journal of Structural Engineering,ASCE,1985,111(4):722-739.
[6]張志偉.SRC柱變形能力及恢復力特性試驗研究[D].泉州:華僑大學,2007.
[7]林煌斌,王全鳳.HRB400高強度混凝土柱抗震性能分析[J].建筑結(jié)構(gòu)學報,2008,29(S1):36-41.
[8]林煌斌,王全鳳.基于IDRAC損傷結(jié)構(gòu)彈塑性時程分析與損傷評價[J].建筑結(jié)構(gòu)學報,2008,29(S1):255-259.
Analysis and Experiment of New Indices to Describe Dynamic Unloading Rigidity of Concrete Structure
LIN Huang-bin1,2,3,WANG Quan-feng1
(1.College of Civil Engineering,Huaqiao University,Quanzhou 362021,China;2.School of Aerospace Engineering and Applied Mechanics,Tongji University,Shanghai 200092,China;3.College of Engineering Technology,Jimei University,Xiamen 361021,China)
Unloading stiffness coefficient and energy correlation coefficient to describe unloading stiffness of structure were presented.The expression of elastic deformation energy due to damage during unloading process was established,which was validated by the experimental data of steel reinforce concrete columns.In elastic-plastic time history analysis of concrete structure,the introduction of unloading stiffness coefficient and energy correlation coefficient is conducive to evaluate quantitively the dynamic unloading,to make up the deficiency without considerring the interacting effects between unloading stiffness and motion state,to solve the disequilibrium of dynamic equation between the tangent stiffness and secant stiffness.
description index;unloading stiffness;concrete structure;restoring-force model;elastic deformation energy
黃曉楠 英文審校:方德平)
TU 313
A
1000-5013(2011)06-0676-04
2011-10-27
林煌斌(1981-),男,講師,主要計算機在結(jié)構(gòu)工程中應用的研究.E-mail:linhuangbin@jmu.edu.cn.
國家自然科學基金資助項目(50578066);高等學校博士學科點專項科研基金資助項目(200803850001);福建省自然科學基金資助項目(2011J01320)