張望喜，王冠杰，龐博，易偉建

（湖南大學(xué) a. 土木工程學(xué)院； b. 工程結(jié)構(gòu)損傷診斷湖南省重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)沙 410082）

裝配式混凝土結(jié)構(gòu)由于具有產(chǎn)品質(zhì)量好、施工效率高、能源消耗小、耗費(fèi)人工少等方面的優(yōu)點(diǎn)，成為我國(guó)在新的工業(yè)時(shí)代重點(diǎn)發(fā)展的主要結(jié)構(gòu)體系之一[1]。節(jié)點(diǎn)的構(gòu)造形式?jīng)Q定了裝配式結(jié)構(gòu)的抗震性能以及施工效率，根據(jù)裝配工藝的不同可分為濕式連接節(jié)點(diǎn)和干式連接節(jié)點(diǎn)[2]。

目前，對(duì)裝配式混凝土結(jié)構(gòu)的試驗(yàn)研究集中在構(gòu)件層面，如節(jié)點(diǎn)、單獨(dú)的梁柱構(gòu)件等[3-4]。然而，構(gòu)件層面的試驗(yàn)難以準(zhǔn)確確定其邊界條件，特別是地震等復(fù)雜的自然災(zāi)害荷載作用下其邊界約束條件比較復(fù)雜，得到的構(gòu)件抗震試驗(yàn)結(jié)果不能完全等同于其在整個(gè)結(jié)構(gòu)中的抗震性能[5]。有限元理論方法的快速發(fā)展及其在現(xiàn)代計(jì)算機(jī)的應(yīng)用，使利用有限元軟件對(duì)大型混凝土結(jié)構(gòu)進(jìn)行輔助抗震分析成為了可能。國(guó)內(nèi)外學(xué)者通常采用實(shí)體單元或桿系單元的建模方式對(duì)裝配式混凝土結(jié)構(gòu)進(jìn)行數(shù)值模擬分析。Kaya等[6]基于ANSYS軟件進(jìn)行了預(yù)應(yīng)力連接裝配式節(jié)點(diǎn)的抗震性能分析；Hawileh等[7]采用ANSYS軟件的接觸單元來(lái)模擬結(jié)合面的接觸行為，考察了裝配式節(jié)點(diǎn)性能對(duì)框架整體性能的影響；張銳[8]基于Marc軟件，采用彈簧單元模擬混凝土界面連接行為和界面處的連接鋼筋，分析了樓板對(duì)框架結(jié)構(gòu)的影響；譚平等[9]基于OpenSees軟件中的宏觀節(jié)點(diǎn)模型分析了裝配式隔震節(jié)點(diǎn)的滯回性能；祝輝慶等[10]使用SAP2000軟件對(duì)采用雙線性彈簧單元模擬節(jié)點(diǎn)連接的裝配式混凝土框架進(jìn)行彈性時(shí)程分析，驗(yàn)證了這種簡(jiǎn)化節(jié)點(diǎn)連接模型的計(jì)算精度能滿足一般工程要求。

實(shí)體單元能精確、直觀地反映結(jié)構(gòu)失效的微觀破壞過(guò)程，但建模過(guò)程復(fù)雜，計(jì)算成本較大，難以用于復(fù)雜結(jié)構(gòu)的分析；桿系單元使用方便，運(yùn)算速度快，但只能得到結(jié)構(gòu)的宏觀反應(yīng)，不能直觀反映結(jié)構(gòu)關(guān)鍵區(qū)域的受力狀態(tài)。多尺度建模方法是在模型的關(guān)鍵區(qū)域采用精細(xì)模型，在非關(guān)鍵區(qū)域設(shè)置宏觀模型，通過(guò)在界面處耦合兩部分來(lái)模擬受力復(fù)雜的構(gòu)件或結(jié)構(gòu)[11]。裝配整體式混凝土框架結(jié)構(gòu)節(jié)點(diǎn)在地震中處于復(fù)雜受力狀態(tài)，往往是結(jié)構(gòu)破壞的薄弱位置，且存在預(yù)制混凝土與現(xiàn)澆混凝土之間接觸摩擦等局部高度非線性失效行為[12]，而其他采用如鋼板焊接、螺栓連接、預(yù)應(yīng)力等裝配工藝的裝配式節(jié)點(diǎn)也存在節(jié)點(diǎn)區(qū)域受力復(fù)雜的問(wèn)題。因此，多尺度建模的策略適合用來(lái)分析裝配式混凝土結(jié)構(gòu)在地震中的破壞模式。近年來(lái)，多尺度計(jì)算開(kāi)始逐漸應(yīng)用于鋼筋混凝土結(jié)構(gòu)[13]、木結(jié)構(gòu)[14]、鋼框架結(jié)構(gòu)[15]、大跨度懸索橋[16-17]等的結(jié)構(gòu)分析。

基于多尺度模型，本研究從裝配整體式混凝土梁柱子結(jié)構(gòu)算例出發(fā)，借助ABAQUS有限元軟件，驗(yàn)證了模型的有效性和準(zhǔn)確性，進(jìn)而對(duì)裝配整體式框架結(jié)構(gòu)進(jìn)行抗震分析，并與現(xiàn)澆框架結(jié)構(gòu)進(jìn)行對(duì)比。通過(guò)求2種結(jié)構(gòu)在靜力推覆作用和動(dòng)力時(shí)程作用下的地震反應(yīng)和損傷情況，量化了地震作用對(duì)裝配整體式框架結(jié)構(gòu)和現(xiàn)澆框架結(jié)構(gòu)的影響。

1 多尺度建模

不同尺度單元界面的有效連接是結(jié)構(gòu)整體協(xié)同工作的關(guān)鍵。常用的連接類別[18]有：1）梁?jiǎn)卧c殼單元；2）梁?jiǎn)卧c實(shí)體單元；3）殼單元與實(shí)體單元。由于連接界面處宏觀模型的節(jié)點(diǎn)少，而精細(xì)模型的節(jié)點(diǎn)多，連接的關(guān)鍵是通過(guò)尋找合適的方法耦合不同單元類型節(jié)點(diǎn)的自由度來(lái)實(shí)現(xiàn)多尺度模型間的變形協(xié)調(diào)和力的傳遞。以圖1中實(shí)體單元和梁?jiǎn)卧倪B接為例，設(shè)Ai(i=1,2,3,…,n-1,n)為實(shí)體單元在連接界面的任意節(jié)點(diǎn)，B為梁?jiǎn)卧谶B接界面的節(jié)點(diǎn)。節(jié)點(diǎn)Ai的位移為

圖1 實(shí)體單元和梁?jiǎn)卧缑婀?jié)點(diǎn)關(guān)系Fig. 1 Interface node relationship between solid element and beam element

可見(jiàn)，實(shí)體單元和梁?jiǎn)卧谧杂啥壬洗嬖诓町?，梁?jiǎn)卧?jié)點(diǎn)B除了3個(gè)平動(dòng)自由度外還存在3個(gè)方向的轉(zhuǎn)動(dòng)自由度。ABAQUS軟件中提供了一點(diǎn)控制多點(diǎn)的多點(diǎn)約束（MPC）方式，選擇其中的梁?jiǎn)卧愋图s束，可以自動(dòng)生成一個(gè)節(jié)點(diǎn)與多個(gè)節(jié)點(diǎn)之間的剛性梁，將第一個(gè)節(jié)點(diǎn)的位移和旋轉(zhuǎn)約束作為其余多個(gè)節(jié)點(diǎn)的位移和旋轉(zhuǎn)，同時(shí)耦合接觸節(jié)點(diǎn)的自由度，從而實(shí)現(xiàn)不同尺度單元之間的變形協(xié)調(diào)，并確保軸力、彎矩和扭矩的正確傳遞[19]。截面采用MPC梁?jiǎn)卧愋图s束的軸向變形和彎曲變形滿足以下n個(gè)方程：

2 多尺度模型驗(yàn)證

為驗(yàn)證裝配式框架結(jié)構(gòu)建模方法以及不同尺度單元連接方法的正確性，基于裝配式梁柱子結(jié)構(gòu)靜力試驗(yàn)數(shù)據(jù)[20]建立了有限元模型。該試驗(yàn)節(jié)點(diǎn)取自一棟6層框架結(jié)構(gòu)，采用2/3的縮尺比例，對(duì)其中錨固方式為90°彎鉤的PC-0-B試件進(jìn)行數(shù)值模擬分析。圖2為PC-0-B試件的尺寸及配筋，梁和柱的混凝土強(qiáng)度等級(jí)為C30；構(gòu)件截面尺寸：柱350 mm×350 mm，梁200 mm×380 mm；縱筋為HRB400，箍筋為HPB300；具體材性見(jiàn)表1。試驗(yàn)裝置見(jiàn)圖3，在梁端分2個(gè)加載階段進(jìn)行單調(diào)加載：第1階段為荷載控制加載，每級(jí)加載5 kN直至試件屈服；第2階段為位移控制加載，每級(jí)加載3 mm直至試件破壞。

表1 鋼筋及混凝土材性特性Table 1 Material properties of rebars and concrete

圖2 試件PC-0-B尺寸及配筋（單位: mm）Fig. 2 Size and reinforcement of PC-0-B specimen（unit: mm）

圖3 試驗(yàn)加載裝置Fig. 3 Test setup

基于第1節(jié)的界面連接原理，在梁柱子結(jié)構(gòu)破壞的關(guān)鍵節(jié)點(diǎn)區(qū)域用實(shí)體單元模擬，非關(guān)鍵區(qū)域的梁柱部分用梁?jiǎn)卧M，通過(guò)MPC梁?jiǎn)卧愋图s束在連接界面處將梁?jiǎn)卧膯我还?jié)點(diǎn)與實(shí)體單元的所有節(jié)點(diǎn)進(jìn)行多點(diǎn)耦合?；炷敛捎镁€性縮減積分單元C3D8R模擬，材料本構(gòu)采用混凝土CDP模型，本構(gòu)方程與損傷因子的計(jì)算采用《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB 50010-2010）[21]附錄C中的相關(guān)公式，泊松比取0.2。鋼筋采用桁架單元T3D2模擬，本構(gòu)采用雙折線模型，泊松比均取0.3，通過(guò)ABAQUS提供的Embedded功能將鋼筋嵌入到整個(gè)模型當(dāng)中。預(yù)制梁端面和柱之間以及預(yù)制梁和現(xiàn)澆部分的疊合面之間采用接觸對(duì)進(jìn)行處理，在法線方向上采用硬接觸，切線方向采用庫(kù)倫摩擦模型。預(yù)制梁端面的摩擦系數(shù)取1.4，預(yù)制梁-現(xiàn)澆層疊合面的摩擦系數(shù)取1.0[22]。宏觀單元采用梁?jiǎn)卧狟31模擬，通過(guò)在inp文件中添加*REBAR語(yǔ)句在梁?jiǎn)卧胁迦脘摻?，與混凝土形成統(tǒng)一截面。梁?jiǎn)卧谢炷梁弯摻畹谋緲?gòu)模型通過(guò)調(diào)用PQ-Fiber[23]子程序，采用其中的UConcrete02和USteel03兩組材料本構(gòu)。

采用位移加載的方式對(duì)實(shí)體單元模型（entity unit model，EUM）和多尺度模型（multi-scale model，MSM）進(jìn)行計(jì)算。文獻(xiàn)[24]表明：在地震作用下，混凝土結(jié)構(gòu)的梁、柱構(gòu)件塑性鉸長(zhǎng)度約為0.4h~1.4h（h為構(gòu)件截面高度）。因此，本文中建立了4個(gè)模型：EUM和梁端計(jì)算破壞區(qū)ld（連接界面到梁端距離）為h、1.5h和2.0h的MSM-1、MSM-2和MSM-3，進(jìn)行模擬計(jì)算，并與試驗(yàn)結(jié)果進(jìn)行對(duì)比。裝配整體式混凝土梁柱子結(jié)構(gòu)的EUM和MSM有限元模型如圖4所示，根據(jù)試驗(yàn)實(shí)際約束條件，縮短柱的長(zhǎng)度以節(jié)省計(jì)算時(shí)間。

圖4 有限元模型Fig. 4 Finite element models

圖5為模擬獲得的荷載-位移曲線與試驗(yàn)的實(shí)測(cè)曲線對(duì)比?？梢?jiàn)，EUM的計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合良好，屈服荷載稍大于試驗(yàn)值；MSM-1、MSM-2和MSM-3三個(gè)模型的峰值荷載與EUM較接近，在彈性階段，隨著梁端計(jì)算破壞區(qū)減小，MSM的剛度有逐漸增大的趨勢(shì)。

圖5 試件PC-0-B數(shù)值模擬與試驗(yàn)曲線對(duì)比Fig. 5 Comparison of numerical simulation and experimental curves of PC-0-B specimen

表2為不同模型的計(jì)算結(jié)果及計(jì)算效率的對(duì)比分析。由表2可以看出，MSM相對(duì)于EUM能夠在保證足夠的計(jì)算精度的前提下，盡可能地提高有限元分析的計(jì)算效率。MSM-2花費(fèi)時(shí)間較少，計(jì)算精度較高，這說(shuō)明了當(dāng)計(jì)算破壞區(qū)域合理時(shí)，多尺度模型可以在計(jì)算效率和計(jì)算精度之間取得平衡。由此可見(jiàn)，隨著有限元計(jì)算模型增大，多尺度計(jì)算的優(yōu)勢(shì)越來(lái)越明顯。

表2 計(jì)算結(jié)果與計(jì)算效率對(duì)比Table 2 Comparison of calculation results and calculation efficiency

3 裝配式混凝土框架結(jié)構(gòu)多尺度模型

采用盈建科軟件設(shè)計(jì)一棟簡(jiǎn)單規(guī)則的6層4×2跨裝配整體式混凝土空間框架結(jié)構(gòu)，取其中的一榀框架作為研究對(duì)象?？蚣馨丛O(shè)防烈度7度（0.15g）、第一設(shè)計(jì)分組、Ⅱ類場(chǎng)地土、三級(jí)抗震等級(jí)進(jìn)行設(shè)計(jì)。樓面恒荷載為5.0 kN/m2，樓面活荷載為2.0 kN/m2；屋面恒荷載為7.0 kN/m2，活荷載0.5 kN/m2；梁上線荷載取16.0 kN/m。結(jié)構(gòu)縱橫向跨度為5.4 m，首層層高4.0 m，其余5層層高3.6 m，建筑總高度22.0 m；梁截面尺寸為300 mm×500 mm，柱截面尺寸為500 mm×500 mm，板厚為150 mm?？蚣苤鶠楝F(xiàn)澆柱；框架梁為疊合構(gòu)件，預(yù)制部分高350 mm，后澆部分高150 mm。第1層柱和第2~6層柱混凝土強(qiáng)度等級(jí)分別為C40、C30；梁和板混凝土強(qiáng)度等級(jí)為C30；梁、柱和板的縱筋為HRB400，箍筋為HPB300。該結(jié)構(gòu)的平面布置如圖6所示，取陰影區(qū)域的一榀框架進(jìn)行計(jì)算。

圖6 結(jié)構(gòu)平面圖（單位: mm）Fig. 6 Structure layout plan（unit: mm）

使用第2節(jié)的建模方法分別建立了尺寸和配筋相同的裝配式結(jié)構(gòu)和現(xiàn)澆結(jié)構(gòu)的多尺度模型，取框架梁、柱的實(shí)體單元計(jì)算范圍為ld=1.5h，框架結(jié)構(gòu)模型尺寸和配筋如圖7。裝配式框架結(jié)構(gòu)模型中混凝土與鋼筋本構(gòu)、單元類型、相互作用等均與第2節(jié)所述相同；現(xiàn)澆混凝土結(jié)構(gòu)的模擬方式與裝配式框架結(jié)構(gòu)的不同之處在于現(xiàn)澆混凝土不存在預(yù)制部分，沒(méi)有預(yù)制部分和現(xiàn)澆部分之間的結(jié)合面，不需要設(shè)置接觸行為。

圖7 裝配整體式框架結(jié)構(gòu)模型尺寸和配筋（單位: mm）Fig. 7 Dimensions and reinforcement of monolithic precast reinforced concrete frame structure model（unit: mm）

4 主要計(jì)算結(jié)果及分析

4.1 模態(tài)分析

裝配整體式混凝土框架結(jié)構(gòu)和現(xiàn)澆混凝土框架結(jié)構(gòu)的動(dòng)力特性如表3所示，對(duì)2個(gè)模型的前6階自振頻率進(jìn)行對(duì)比。從表3中可以看出，裝配整體式混凝土框架結(jié)構(gòu)和現(xiàn)澆混凝土框架結(jié)構(gòu)的前6階模態(tài)基本一致，頻率相差最大值為3.7%，最小僅0.6%。裝配式框架的自振頻率略小于現(xiàn)澆框架，這是因?yàn)檠b配整體式節(jié)點(diǎn)存在新舊混凝土之間的連接接縫，導(dǎo)致結(jié)構(gòu)的剛度偏小。

表3 框架結(jié)構(gòu)自振頻率對(duì)比Table 3 Comparison of natural frequencies of frame structures

4.2 靜力推覆分析

為了解結(jié)合面接縫對(duì)裝配整體式框架抗側(cè)剛度的影響，采用圖8所示的倒三角形荷載分布模式和均勻荷載分布模式分別對(duì)裝配整體式混凝土框架結(jié)構(gòu)和現(xiàn)澆混凝土結(jié)構(gòu)進(jìn)行單向推覆分析。圖9為分析得到的2種荷載分布形式下各模型的能力曲線。

圖8 荷載分布模式Fig. 8 Load distribution mode

圖9 基底剪力-頂點(diǎn)位移曲線Fig. 9 Curves of base shear force to top displacement

可見(jiàn)，荷載較小時(shí)，裝配整體式混凝土框架結(jié)構(gòu)和現(xiàn)澆混凝土框架結(jié)構(gòu)的基底剪力與頂點(diǎn)位移呈線性關(guān)系；隨著荷載增加，2種結(jié)構(gòu)的靜力抗側(cè)剛度出現(xiàn)明顯退化，裝配整體式框架的靜力抗側(cè)剛度始終小于現(xiàn)澆框架。裝配式框架在屈服點(diǎn)處的基底剪力更小，相應(yīng)的頂點(diǎn)位移也更小。裝配式結(jié)構(gòu)和現(xiàn)澆結(jié)構(gòu)在均勻荷載分布模式下的基底剪力最大值相較于倒三角分布模式降低了14.95%和15.57%。

4.3 結(jié)構(gòu)動(dòng)力時(shí)程分析

為進(jìn)一步研究罕遇地震作用下2種框架結(jié)構(gòu)的抗震性能，選取2條天然波（Chi-Chi地震波和Imperial Valley地震波）和1條人工波進(jìn)行7度罕遇地震（310 gal）作用下的動(dòng)力彈塑性時(shí)程分析，動(dòng)力時(shí)程求解算法為Newmark-β法。各地震波的時(shí)程曲線均按照比例將峰值轉(zhuǎn)換成310 cm/s2，結(jié)構(gòu)阻尼比取0.05。提取各地震波輸入20 s內(nèi)2種框架模型頂層的位移，頂層位移時(shí)程曲線見(jiàn)圖10，框架結(jié)構(gòu)各層最大側(cè)向位移見(jiàn)表4。

表4 框架結(jié)構(gòu)各層最大側(cè)向位移Table 4 Maximum lateral displacement of each floor of frame structure mm

圖10 框架頂層位移時(shí)程對(duì)比Fig. 10 Comparison of displacement time history of framework top floor

可見(jiàn)，在地震輸入的前8 s內(nèi)，裝配整體式框架的頂點(diǎn)位移與現(xiàn)澆結(jié)構(gòu)的頂點(diǎn)位移基本一致，隨著地震時(shí)長(zhǎng)增加，裝配整體式框架的頂點(diǎn)位移響應(yīng)與現(xiàn)澆結(jié)構(gòu)的頂點(diǎn)位移響應(yīng)差距不斷加大，這是因?yàn)樵诘卣鹣卵b配整體式框架的損傷更大，其側(cè)向剛度降低，變形增大。裝配整體式框架結(jié)構(gòu)相對(duì)于現(xiàn)澆結(jié)構(gòu)，框架頂層最大位移平均值增大了3.8%。

各地震波作用下2個(gè)框架結(jié)構(gòu)模型的各樓層最大層間位移角平均值的對(duì)比如圖11所示。框架結(jié)構(gòu)最大層間位移角平均值均出現(xiàn)在第2層，裝配整體式結(jié)構(gòu)的最大層間位移角平均值與現(xiàn)澆結(jié)構(gòu)相比在第1、2、3層差距較大，在其他層數(shù)差距不太明顯，最大層間位移角平均值在第1、2、3層相差達(dá)到了7.18%、3.80%、3.08%。

圖11 最大層間位移角對(duì)比Fig. 11 Comparison of maximum inter-story drift ratios

為進(jìn)一步研究節(jié)點(diǎn)的破壞程度，選取裝配整體式框架結(jié)構(gòu)模型與現(xiàn)澆結(jié)構(gòu)模型損傷明顯的薄弱層節(jié)點(diǎn)（第2層中節(jié)點(diǎn)）進(jìn)行分析。圖12和圖13分別為2種框架在Chi-Chi波作用20 s后的混凝土受拉損傷和受壓損傷分布云圖?？梢?jiàn)，在相同地震荷載作用下不論是受拉還是受壓，裝配整體式節(jié)點(diǎn)的損傷程度均大于現(xiàn)澆節(jié)點(diǎn)，且節(jié)點(diǎn)破壞主要集中在梁端部，說(shuō)明計(jì)算模型能夠滿足“強(qiáng)柱弱梁、強(qiáng)節(jié)點(diǎn)弱構(gòu)件”的設(shè)計(jì)要求。現(xiàn)澆結(jié)構(gòu)的混凝土損傷分布比較均勻，與之相比，裝配整體式框架結(jié)構(gòu)的混凝土損傷更多地集中在預(yù)制梁后澆部位。此外，從2種框架整體上看，1~3層節(jié)點(diǎn)的損傷程度大于4~5層節(jié)點(diǎn)，這說(shuō)明需要特別重視1~3層節(jié)點(diǎn)的施工質(zhì)量。

圖12 裝配式混凝土節(jié)點(diǎn)損傷云圖Fig. 12 Damage cloud of monolithic precast concrete joint

圖13 現(xiàn)澆混凝土節(jié)點(diǎn)損傷云圖Fig. 13 Damage cloud of cast-in-place concrete joint

隨著裝配式混凝土節(jié)點(diǎn)構(gòu)造復(fù)雜性提高，形式越發(fā)多樣，采用多尺度模型可以有效模擬不同節(jié)點(diǎn)形式的裝配式混凝土結(jié)構(gòu)整體抗震性能，同時(shí)可以直觀反映裝配式混凝土結(jié)構(gòu)損傷演變、裂縫發(fā)展等過(guò)程，得到結(jié)構(gòu)的薄弱部分以便于進(jìn)一步改進(jìn)設(shè)計(jì)方法。

5 結(jié) 論

基于多尺度建模對(duì)裝配整體式混凝土節(jié)點(diǎn)和框架結(jié)構(gòu)進(jìn)行了數(shù)值模擬，分析了裝配整體式框架結(jié)構(gòu)和現(xiàn)澆混凝土框架結(jié)構(gòu)在單向推覆作用和7度罕遇地震作用下的抗震性能，得出以下結(jié)論：

1）通過(guò)算例驗(yàn)證可知多尺度建模方法可以正確模擬裝配整體式混凝土結(jié)構(gòu)的受力性能，能同時(shí)滿足計(jì)算效率和計(jì)算精度的要求。

2）靜力彈塑性分析表明，裝配整體式框架結(jié)構(gòu)的靜力抗側(cè)剛度小于現(xiàn)澆結(jié)構(gòu)；在不同荷載分布模式作用下，結(jié)構(gòu)的基底剪力最大值差值為15%左右。

3）動(dòng)力彈塑性時(shí)程分析表明，在相同地震波作用下裝配整體式框架結(jié)構(gòu)的基本周期和動(dòng)力特性與現(xiàn)澆結(jié)構(gòu)相近，薄弱部位都位于第2層。裝配整體式框架結(jié)構(gòu)相對(duì)于現(xiàn)澆結(jié)構(gòu)頂層最大位移增加了3.8%，節(jié)點(diǎn)損傷更加嚴(yán)重。

4）多尺度建模方法在裝配式混凝土結(jié)構(gòu)分析中具有較好的適用性。