子模型技術(shù)在帶伸臂框架-核心筒結(jié)構(gòu)分析中的應(yīng)用研究

何玉明

(四川建筑職業(yè)技術(shù)學(xué)院， 四川德陽(yáng) 618000)

子模型技術(shù)在帶伸臂框架-核心筒結(jié)構(gòu)分析中的應(yīng)用研究

何玉明

(四川建筑職業(yè)技術(shù)學(xué)院， 四川德陽(yáng) 618000)

研究了ABAQUS子模型技術(shù)在帶伸臂框-筒結(jié)構(gòu)有限元分析中的精確程度，提出用子模型分析復(fù)雜高層建筑關(guān)鍵部位，實(shí)現(xiàn)建筑結(jié)構(gòu)精細(xì)化有限元分析。闡述了子模型技術(shù)的原理、子模型驅(qū)動(dòng)條件獲取的方法、子模型分析計(jì)算流程。計(jì)算了帶伸臂框-筒結(jié)構(gòu)地震作用下的變形、受力特點(diǎn)，對(duì)比了靜力及動(dòng)力作用下子模型的計(jì)算精度。計(jì)算結(jié)果顯示：(1)地震作用下帶伸臂框-筒結(jié)構(gòu)加強(qiáng)層上下樓層層間位移角突變，墻、柱應(yīng)力集中；(2)靜力作用下子模型計(jì)算結(jié)果可靠，采用精細(xì)化子模型能得到更詳細(xì)的計(jì)算結(jié)果；(3)子模型與全局模型動(dòng)力彈性時(shí)程分析結(jié)果高度吻合，子模型技術(shù)能用于高層建筑地震反應(yīng)分析。

加強(qiáng)層； 子模型技術(shù)； 精細(xì)模型； 抗震分析； 仿真模擬

目前，基于有限元法的仿真分析是結(jié)構(gòu)體系復(fù)雜的異形建筑及高層、超高層建筑結(jié)構(gòu)抗震及抗風(fēng)性能預(yù)測(cè)和評(píng)估的重要手段[1]。通用有限元軟件ABAQUS提供了豐富的單元庫(kù)和材料庫(kù)，對(duì)線(xiàn)性分析和非線(xiàn)性分析問(wèn)題都具有很好的適用性，在高層、大跨建筑結(jié)構(gòu)和大型橋梁結(jié)構(gòu)的抗震分析領(lǐng)域的應(yīng)用日趨廣泛[2]。然而，ABAQUS前處理過(guò)程繁瑣，建立精細(xì)化的全局模型困難，制約了在復(fù)雜建筑結(jié)構(gòu)分析中的普遍應(yīng)用。結(jié)構(gòu)分析最終更關(guān)心的是關(guān)鍵部位的分析結(jié)果，采用等效替代的方式建立相對(duì)粗略的整體模型，然后選取關(guān)鍵的部位建立精細(xì)化子模型分析，使結(jié)構(gòu)數(shù)值模擬“主次分明”，同時(shí)在滿(mǎn)足分析精度的前提下，也能降低對(duì)計(jì)算機(jī)硬件的要求，實(shí)現(xiàn)“小機(jī)器解決大問(wèn)題”的目的。

徐偉等[3]采用子模型技術(shù)對(duì)大跨徑斜拉橋橋面結(jié)構(gòu)進(jìn)行分析，結(jié)果表明，獨(dú)立梁段分析橋面板豎向變形結(jié)果與子模型法較接近，認(rèn)為子模型法是進(jìn)行此類(lèi)分析的較有效方法。謝素明等[4]運(yùn)用子模型技術(shù)對(duì)汽車(chē)鑄鋼側(cè)架進(jìn)行分析，計(jì)算所得應(yīng)力分量與實(shí)驗(yàn)結(jié)果對(duì)比，精度很高。張崗等[5]采用3種子模型實(shí)施方法對(duì)混凝土箱梁懸臂板進(jìn)行分析，得出子模型技術(shù)適合于混凝土箱梁懸臂板局部精細(xì)分析，且實(shí)體單元分析結(jié)果比板殼單元更接近實(shí)驗(yàn)結(jié)果。宋良豐等[6]以大崗山拱壩為例，采用子模型技術(shù)對(duì)比分析了不同網(wǎng)格尺寸下孔口區(qū)域應(yīng)力分布及損傷開(kāi)裂情況，計(jì)算結(jié)果顯示子模型與全局模型計(jì)算結(jié)果吻合良好。

帶伸臂框架-核心筒結(jié)構(gòu)體系[7]是從核心筒中伸出剛度較大的伸臂構(gòu)件，并與外部框架鉸接(半剛接)，形成加強(qiáng)層，從而協(xié)調(diào)核心筒和外部框架共同作用，此類(lèi)結(jié)構(gòu)體系能高效地減小結(jié)構(gòu)在水平荷載作用下的位移，實(shí)際工程中被廣泛采用[8]。但是加強(qiáng)層附近樓層剛度不連續(xù)，并有明顯突變，在地震作用下，結(jié)構(gòu)的破壞容易集中在加強(qiáng)層附近，形成薄弱層[9]。本文以帶伸臂桁架的框架-核心筒結(jié)構(gòu)為研究對(duì)象，研究子模型技術(shù)在復(fù)雜建筑結(jié)構(gòu)分析中的應(yīng)用。

1 子模型技術(shù)

1.1 子模型技術(shù)的基本原理

子模型是得到模型部分區(qū)域中更加精確解的有限元技術(shù)，子模型方法又稱(chēng)為切割邊界位移法或特定邊界位移法[10]。子模型獨(dú)立于整體模型，需保證子模型與全局模型中對(duì)應(yīng)部分空間位置保持一致，子模型建模過(guò)程不受全局模型的影響，能選取任意關(guān)鍵樓層、關(guān)鍵的構(gòu)件和關(guān)鍵節(jié)點(diǎn)建立精細(xì)化子模型，然后調(diào)用存儲(chǔ)在全局模型計(jì)算結(jié)果文件中子模型切割邊界上的節(jié)點(diǎn)(面)對(duì)應(yīng)時(shí)刻的位移(應(yīng)力)作為子模型的驅(qū)動(dòng)變量，實(shí)際上就是將全局模型對(duì)子模型部分的作用以驅(qū)動(dòng)變量代替。根據(jù)圣維南原理可知，將全局模型對(duì)子模型的作用以驅(qū)動(dòng)變量等效替代后，子模型的應(yīng)力和應(yīng)變只會(huì)在驅(qū)動(dòng)邊界附近有改變，遠(yuǎn)離驅(qū)動(dòng)邊界的位置的應(yīng)力和應(yīng)變不受影響。因此，通過(guò)適當(dāng)選取子模型的切割邊界，就可以保證子模型中所關(guān)心的部位的分析結(jié)果是精確的(圖1)。

圖1 子模型生成示意

1.2 ABAQUS子模型驅(qū)動(dòng)條件的獲取

ABAQUS子模型技術(shù)提供了兩種驅(qū)動(dòng)方法：(1)基于節(jié)點(diǎn)的子模型，以節(jié)點(diǎn)位移為驅(qū)動(dòng)變量；(2)基于面的子模型，以驅(qū)動(dòng)面的應(yīng)力為驅(qū)動(dòng)變量。后者只能用于采用實(shí)體單元模型，且限于彈性分析，因此，工程上較少使用?；诠?jié)點(diǎn)的子模型，驅(qū)動(dòng)節(jié)點(diǎn)不能位于總體模型中只有一維單元(桁架單元、梁?jiǎn)卧?、軸對(duì)稱(chēng)殼單元等)的區(qū)域，因?yàn)樵谶@些區(qū)域中單元形函數(shù)的插值多項(xiàng)式階次低，無(wú)法獲得必須的插值結(jié)果(圖2)。

圖2 驅(qū)動(dòng)變量插值示意

根據(jù)單元形函數(shù)可以通過(guò)單元節(jié)點(diǎn)位移插值得到驅(qū)動(dòng)節(jié)點(diǎn)的位移值(驅(qū)動(dòng)變量)，以四節(jié)點(diǎn)二維矩形線(xiàn)性單元為例，說(shuō)明插值過(guò)程。

四節(jié)點(diǎn)二維矩形線(xiàn)性單元的形函數(shù)為：

[N]=[(1-ξ)(1-η)/4, (1+ξ)(1-η)/4,

(1+ξ)(1+η)/4, (1-ξ)(1+η)/4]T

(1)

如圖2所示驅(qū)動(dòng)節(jié)點(diǎn)i在t時(shí)刻的位移為：

[Di(t)]=[ui(t),vi(t)]T

(2)

其值通過(guò)母單元(驅(qū)動(dòng)節(jié)點(diǎn)所處的單元)節(jié)點(diǎn)的位移值插值得到，

(3)

[Dj(t)]=[uj(t),vj(t)]T

(4)

其中，Nj為母單元的形函數(shù)，Dj(t)為母單元j節(jié)點(diǎn)t時(shí)刻的位移，由全局模型計(jì)算得到。

有限元是基于位移法求解，平衡方程為：

[K][D]=[P]

(5)

其中，[K]為結(jié)構(gòu)總剛度矩陣，[P]為外荷載向量，[D]為位移向量。對(duì)于子模型，位移向量[D]中，驅(qū)動(dòng)節(jié)點(diǎn)位移[D1]已知，其余節(jié)點(diǎn)位移令為[D2]，將式(5)矩陣分塊得到，

(6)

展開(kāi)式(6)得，

[K22][D2]=[P2]-[K21][D1]

(7)

由式(7)知，子模型驅(qū)動(dòng)節(jié)點(diǎn)的位移，成為求解[D2]的荷載項(xiàng)的一部分，即驅(qū)動(dòng)變量進(jìn)入荷載項(xiàng)。實(shí)際上，驅(qū)動(dòng)節(jié)點(diǎn)的位移成為子模型的附加位移邊界條件。

1.3 子模型分析計(jì)算的流程

子模型的計(jì)算分析以全局模型計(jì)算分析結(jié)果為基礎(chǔ)，計(jì)算分析流程如圖3所示。主要通過(guò)子模型與全局模型切割邊界處節(jié)點(diǎn)位移值、子模型相應(yīng)區(qū)域云紋圖變化趨勢(shì)的吻合程度來(lái)判斷子模型是否正確，同時(shí)通過(guò)控制兩者的吻合程度來(lái)保證子模型分析結(jié)果的精度。

圖3 子模型計(jì)算流程

圖4 加強(qiáng)層示意

2 模型介紹

2.1 工程概況

結(jié)構(gòu)計(jì)算模型共40層，層高為3.0 m，占地面積12.96 m×12.96 m，高度120 m，高寬比5.56。在20層和40層設(shè)置加強(qiáng)層(圖4)?？拐鹪O(shè)防烈度為7度，設(shè)計(jì)基本地震加速度為0.15g，設(shè)計(jì)地震分組為第二組。場(chǎng)地類(lèi)別為Ⅱ類(lèi)，特征周期0.40 s。

結(jié)構(gòu)高度滿(mǎn)足GB 50011-2010《建筑抗震設(shè)計(jì)規(guī)范 》和JGJ 3-2010《高層建筑混凝土結(jié)構(gòu)技術(shù)規(guī)程》規(guī)定框架-核心筒結(jié)構(gòu)在7度設(shè)防時(shí)的130級(jí)高度建筑的最大適用高度140 m限值，屬A級(jí)高度高層建筑。1～10層混凝土強(qiáng)度等級(jí)為C50，11～40層混凝土強(qiáng)度等級(jí)為C40，構(gòu)件尺寸如表1所示。采用中國(guó)建筑科學(xué)院開(kāi)發(fā)的SATWE進(jìn)行結(jié)構(gòu)驗(yàn)算，驗(yàn)算結(jié)果滿(mǎn)足我國(guó)現(xiàn)行結(jié)構(gòu)設(shè)計(jì)規(guī)范。

表1 構(gòu)件尺寸 mm

2.2 計(jì)算模型簡(jiǎn)化

計(jì)算模型共40層，高120 m，全樓結(jié)構(gòu)模型如圖5所示。選取設(shè)有伸臂桁架的一榀(軸③)建立有限元計(jì)算模型，研究伸臂桁架對(duì)結(jié)構(gòu)的受力性能的影響(圖6)。各構(gòu)件選用的單元類(lèi)型見(jiàn)表2。

圖5 整樓模型

圖6 二維模型

表2 構(gòu)件單元類(lèi)型

2.3 簡(jiǎn)化模型驗(yàn)證、分析

采用SATWE和ABAQUS分別對(duì)三維模型和簡(jiǎn)化模型進(jìn)行模態(tài)分析，簡(jiǎn)化模型與整體模型的第一振型均為X向平動(dòng)，周期誤差為3 %(表3)。在設(shè)防烈度下，該簡(jiǎn)化模型彈性時(shí)程分析得出的最大層間位移角滿(mǎn)足規(guī)范要求。因此，該簡(jiǎn)化計(jì)算模型是合理的，能較好地反映整體結(jié)構(gòu)在XZ平面內(nèi)的受力性能。

表3 振型對(duì)比

彈性時(shí)程分析表明：地震作用下，中部加強(qiáng)層處層間位移角發(fā)生明顯突變這是由于伸臂桁架增大了加強(qiáng)層的剛度，導(dǎo)致上下樓層剛度突變，而加強(qiáng)層的柱端約束增強(qiáng)更為明顯(圖7)。

3 子模型計(jì)算精度驗(yàn)證

3.1 子模型選取

設(shè)置加強(qiáng)層后，豎向剛度發(fā)生突變，加強(qiáng)層上下樓層為薄弱層，本文選取17～23層建立細(xì)化網(wǎng)格的子模型，如圖6所示。通過(guò)對(duì)比子模型與全局模型的云紋圖、節(jié)點(diǎn)位移、單元應(yīng)力等指標(biāo)來(lái)研究子模型技術(shù)在建筑結(jié)構(gòu)分析中精確程度。

3.2 靜力分析驗(yàn)證

對(duì)模型施加倒三角形X向水平荷載(模擬風(fēng)荷載)，最大值120 kN/m。子模型計(jì)算得出的右柱角點(diǎn)X向位移誤差在0.05 %以?xún)?nèi)，由此計(jì)算得出的層間位移角誤差在0.5 %以?xún)?nèi)，圖8可以看出，子模型與全局模型靜力分析結(jié)果吻合度高。由圖9、圖10可知，全局模型和子模型的云紋圖變化趨勢(shì)相同，由于子模型網(wǎng)格更細(xì)，云紋圖變化更連續(xù)，局部計(jì)算結(jié)果更精確。

圖7 層間位移角曲線(xiàn)

圖8 樓層-層間位移角曲線(xiàn)

圖9 X向位移云圖對(duì)比

圖10 最小主應(yīng)力云圖對(duì)比(混凝土)

3.3 地震反應(yīng)分析結(jié)果驗(yàn)證

選取El-Centro波(東西向數(shù)據(jù))，加速度峰值修正為110 cm/s2，加速度譜見(jiàn)圖11，對(duì)結(jié)構(gòu)進(jìn)行彈性時(shí)程分析。對(duì)比全局模型和子模型節(jié)點(diǎn)位移和單元應(yīng)力，校驗(yàn)子模型的計(jì)算精確程度。為使選取的節(jié)點(diǎn)具有代表性，本文分別選取墻、柱中心節(jié)點(diǎn)對(duì)比校核位移，應(yīng)力對(duì)比則分別選取混凝土單元及鋼桁架單元。

圖11 El-Centro波加速度時(shí)程曲線(xiàn)

圖12 柱節(jié)點(diǎn)位移-時(shí)間曲線(xiàn)

由圖12節(jié)點(diǎn)位移-時(shí)間曲線(xiàn)可以看出，子模型和全局模型計(jì)算的節(jié)點(diǎn)位移在各時(shí)間點(diǎn)上吻合度很高，在位移均勻變化的區(qū)段，子模型和全局模型結(jié)算結(jié)果基本相同，在位移極值點(diǎn)處(速度方向變化)節(jié)點(diǎn)位移誤差最大。計(jì)算表明：子模型計(jì)算誤差在3 %以?xún)?nèi)。

由圖13、圖14單元應(yīng)力歷程曲線(xiàn)可以看出，子模型與全局模型計(jì)算的單元應(yīng)力吻合度較高，誤差變化規(guī)律與節(jié)點(diǎn)位移誤差變化規(guī)律相同，極值點(diǎn)處誤差相對(duì)較大。應(yīng)力最大值、最小值及誤差見(jiàn)表4。

圖13 混凝土單元應(yīng)力-時(shí)間曲線(xiàn)

圖14 鋼桁架單元應(yīng)力-時(shí)間曲線(xiàn)

表4 單元應(yīng)力對(duì)比 Pa

圖15 最小主應(yīng)力云圖對(duì)比(14.2 s時(shí)刻)

對(duì)比可知，子模型與全局模型計(jì)算的位移誤差比應(yīng)力誤差小，這是由于該子模型是以位移為驅(qū)動(dòng)變量，應(yīng)力是由位移計(jì)算而得，導(dǎo)致應(yīng)力誤差累積。鋼材的彈性模量比混凝土大，因此相同的位移差會(huì)導(dǎo)致鋼桁架單元應(yīng)力曲線(xiàn)比混凝土單元吻合度低。 選取如圖15(a)所示框架柱，對(duì)比全局模型和子模型在14.2 s時(shí)刻的最小主應(yīng)力云圖(圖15)，可以看出，動(dòng)力分析中，子模型和全局模型的應(yīng)力云圖分布趨勢(shì)一致，子模型云圖更細(xì)致，能更為詳細(xì)地顯示框架的應(yīng)力分布。

4 結(jié)論

子模型技術(shù)為復(fù)雜高層建筑結(jié)構(gòu)的關(guān)鍵部位有限元分析提供了新途徑，采用子模型分析關(guān)鍵部位可以提高有限元分析效率。子模型極大降低模型體量，因此可以實(shí)現(xiàn)精細(xì)化模型的建立，提高分析精度。本文通過(guò)對(duì)帶伸臂的框架-核心筒簡(jiǎn)化模型計(jì)算分析表明：(1)加強(qiáng)層伸臂使結(jié)構(gòu)豎向剛度不連續(xù)，地震作用下加強(qiáng)層附近層間位移角突變，加強(qiáng)層上下2～3層為薄弱層，應(yīng)重點(diǎn)計(jì)算、設(shè)計(jì)。(2)子模型靜力分析結(jié)果可靠，運(yùn)用子模型技術(shù)分析風(fēng)荷載作用下高層建筑的側(cè)移，能得到精確的結(jié)果。(3)基于節(jié)點(diǎn)的子模型動(dòng)力彈性時(shí)程分析結(jié)果與全局模型吻合度較高，位移誤差在5 %以?xún)?nèi)，應(yīng)力誤差10 %左右，可以運(yùn)用子模型技術(shù)對(duì)建筑結(jié)構(gòu)進(jìn)行地震反應(yīng)分析。

[1] 盧嘯.超高巨柱-核心筒-伸臂結(jié)構(gòu)地震災(zāi)變及抗震性能研究[D].北京:清華大學(xué),2013.

[2] 聶建國(guó),王宇航.ABAQUS中混凝土本構(gòu)模型用于模擬結(jié)構(gòu)靜力行為的比較研究[J] .工程力學(xué),2013,30(4)：59-67.

[3] 徐偉,李智,張肖寧.子模型法在大跨徑斜拉橋橋面結(jié)構(gòu)分析中的應(yīng)用[J].土木工程學(xué)報(bào),2004, 37(6):30-34.

[4] 謝素明,江渡,兆文忠.提高計(jì)算精度的子模型技術(shù)及側(cè)架實(shí)例[J].大連鐵道學(xué)院學(xué)報(bào),2000,21(3):13-16.

[5] 張崗,王新敏,賀拴海.子模型技術(shù)在混凝土箱梁懸臂板分析中的應(yīng)用[J].鄭州大學(xué)學(xué)報(bào): 工學(xué)版, 2007,28(3):28-32.

[6] 宋良豐,武明鑫,王進(jìn)廷,等.拱壩孔口地震反應(yīng)分析的有限元子模型方法[J].水力發(fā)電學(xué)報(bào)，2014,33(3).

[7] 周穎,呂西林,張翠強(qiáng).消能減震伸臂桁架超高層結(jié)構(gòu)抗震性能研究[J].振動(dòng)與沖擊,2011,30(11):186-189.

[8] 汪大綏,周建龍.我國(guó)高層建筑鋼-混凝土組合結(jié)構(gòu)發(fā)展與展望[J].建筑結(jié)構(gòu)學(xué)報(bào),2010,31(6):62-70.

[9] Kidong Park(2010). A Comparison Study of Conventional Construction Methods and Outrigger Damper System for the Compensation of Differential Column Shortening in High-rise Buildings. International Journal of Steel Structures10(4)： 317-324.

[10] Bazant Z P,Cedolin L.Blunt crack band propagation infinite element analysis[J].J.Eng.Mech.Div.ASCE,1979,105(EM2): 297-315.

何玉明(1990～)，男，碩士，助教，從事鋼筋混凝土結(jié)構(gòu)損傷破壞機(jī)理研究。

TU973.1+7

A

[定稿日期]2016-09-28