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

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

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

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

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

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

1 多尺度建模

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

圖1 實體單元和梁單元界面節(jié)點關(guān)系Fig. 1 Interface node relationship between solid element and beam element

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

2 多尺度模型驗證

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

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

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

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

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

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

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

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

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

采用盈建科軟件設(shè)計一棟簡單規(guī)則的6層4×2跨裝配整體式混凝土空間框架結(jié)構(gòu)，取其中的一榀框架作為研究對象。框架按設(shè)防烈度7度（0.15g）、第一設(shè)計分組、Ⅱ類場地土、三級抗震等級進行設(shè)計。樓面恒荷載為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。框架柱為現(xiàn)澆柱；框架梁為疊合構(gòu)件，預(yù)制部分高350 mm，后澆部分高150 mm。第1層柱和第2~6層柱混凝土強度等級分別為C40、C30；梁和板混凝土強度等級為C30；梁、柱和板的縱筋為HRB400，箍筋為HPB300。該結(jié)構(gòu)的平面布置如圖6所示，取陰影區(qū)域的一榀框架進行計算。

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

使用第2節(jié)的建模方法分別建立了尺寸和配筋相同的裝配式結(jié)構(gòu)和現(xiàn)澆結(jié)構(gòu)的多尺度模型，取框架梁、柱的實體單元計算范圍為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ù)制部分，沒有預(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 主要計算結(jié)果及分析

4.1 模態(tài)分析

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

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

4.2 靜力推覆分析

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

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

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

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

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

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

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

圖10 框架頂層位移時程對比Fig. 10 Comparison of displacement time history of framework top floor

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

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

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

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

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

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

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

5 結(jié) 論

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

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

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

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

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