李雙龍 梁曉東 胡守營 王傳甲
(奧意建筑工程設(shè)計有限公司,廣東深圳 518031)
隨著我國經(jīng)濟和社會的發(fā)展,城市人口密度持續(xù)增長,新增用地日益緊張,高層建筑的發(fā)展已經(jīng)成為歷史的必然和時代的潮流[1]。由相關(guān)數(shù)據(jù)統(tǒng)計可知,我國已建的高層建筑中80%是住宅,而在這些高層住宅中鋼筋混凝土剪力墻結(jié)構(gòu)占比在90%以上。所以,剪力墻結(jié)構(gòu)體系的優(yōu)化,比如剪力墻的數(shù)量、布置以及尺寸的優(yōu)化工作,即布置多少剪力墻,以及在什么位置才能最大限度地發(fā)揮剪力墻的作用,這些應(yīng)當作為當前的重點工作予以關(guān)注。有數(shù)據(jù)表明,通過優(yōu)化計算后的結(jié)構(gòu)設(shè)計方案造價一般會比傳統(tǒng)設(shè)計方案節(jié)約5% ~30%[2],因此,研究高層剪力墻結(jié)構(gòu)體系的結(jié)構(gòu)優(yōu)化設(shè)計工作,可使資源更加合理優(yōu)化配置,具有良好的社會和經(jīng)濟效益。
國外學(xué)者在結(jié)構(gòu)優(yōu)化方面的研究起步較早。從二十世紀六七十年代起,國外學(xué)者就進行過剪力墻結(jié)構(gòu)體系彈塑性力學(xué)性能的研究,其中以Kumarapillai K和Pekau C A等為代表,取得了卓著的成績[3,4]。之后隨著計算機信息技術(shù)的發(fā)展,計算能力日益提高,可用于結(jié)構(gòu)優(yōu)化的手段更加多樣化。我國在剪力墻結(jié)構(gòu)優(yōu)化方面的研究有:傅學(xué)怡[5]給出了樓層斷面剪力墻的計算公式;邵志堅[6]利用多組模型進行計算,得出了剪力墻的混凝土強度等級、墻率以及墻洞數(shù)量等因素對計算模型內(nèi)力、剛度、位移等結(jié)果的影響;姚琦[7]通過優(yōu)化計算,更改剪力墻布置獲得了合理的結(jié)構(gòu)剛度,實現(xiàn)了對整個剪力墻結(jié)構(gòu)的優(yōu)化。
以上研究,普遍是從單個剪力墻構(gòu)件模型出發(fā),在一定的條件下,得出構(gòu)件響應(yīng)與變量的影響關(guān)系,大多只實現(xiàn)了單目標下的剪力墻數(shù)量的選擇,對剪力墻哪些位置設(shè)置洞口能達到最優(yōu)目前還未見報道。對于實際的剪力墻結(jié)構(gòu),由于存在設(shè)計變量多、約束條件多、變量離散、受建筑功能限制較大等因素,因此前述文獻中給出的公式對實際結(jié)構(gòu)定量化優(yōu)化缺少可操作性,對實現(xiàn)多目標下的高層剪力墻結(jié)構(gòu)的剪力墻位置和數(shù)量的分析還有待研究。多目標結(jié)構(gòu)優(yōu)化設(shè)計將遺傳算法等先進的方法引入到結(jié)構(gòu)領(lǐng)域,是建筑結(jié)構(gòu)設(shè)計領(lǐng)域未來的發(fā)展方向。
建筑結(jié)構(gòu)設(shè)計中,在滿足建筑功能和國家規(guī)范的前提下,基于“低碳、節(jié)能、環(huán)?!崩砟?,對建筑結(jié)構(gòu)進行優(yōu)化設(shè)計,通過合理布置剪力墻位置以及合理選擇各類構(gòu)件的尺寸,不僅能使結(jié)構(gòu)構(gòu)件受力更加合理,而且可以節(jié)省材料,減小企業(yè)的資金壓力。本文基于優(yōu)化分析軟件Isight和通用有限元軟件Abaqus,通過二次開發(fā)將自編的C#和Python腳本程序嵌入優(yōu)化平臺,運用有限元生死單元技術(shù),搭建了完整的可用于高層剪力墻結(jié)構(gòu)體系優(yōu)化設(shè)計的多目標優(yōu)化流程,將傳統(tǒng)被動設(shè)計變?yōu)榫毣鲃訉ふ蚁鄬ψ顑?yōu)解的過程,可以對設(shè)計師提供強有力的幫助。
數(shù)學(xué)上,多目標優(yōu)化問題可以描述為D個設(shè)計變量、N個目標函數(shù)、m+k個約束條件組成的一個求解問題。決策者根據(jù)具體問題要求在非劣解集中選擇比較滿意的一個非劣解作為最終解。多目標優(yōu)化問題的數(shù)學(xué)形式可以如下描述[8]:
其中,x為一個維度為D的設(shè)計向量;y為目標函數(shù)向量,由多個目標函數(shù)組成;N為優(yōu)化目標的數(shù)量;gi(x)≤0和hj(x)=0確定了解的可行域,xd_max和xd_min為每維向量搜索的上下限。
基于高層剪力墻初始結(jié)構(gòu)設(shè)計方案,在不違反建筑功能和規(guī)范的前提下(比如避免出現(xiàn)短肢剪力墻),在可能的位置上預(yù)先布置剪力墻,然后通過自編接口程序?qū)⒛P娃D(zhuǎn)為Abaqus軟件的輸入文件,通過優(yōu)化分析軟件Isight搭建分析流程,并將自編腳本程序嵌入到流程中,對分組剪力墻單元進行敏度分析。之后根據(jù)敏度分析結(jié)果,剔除不敏感位置墻元(不敏感意味著此位置增加或減少墻對結(jié)構(gòu)指標比如層間位移角影響不大),從而達到剪力墻準拓撲相對最優(yōu)分布。
基于上述結(jié)果,將多目標算法和結(jié)構(gòu)、構(gòu)件控制指標腳本嵌入優(yōu)化平臺,設(shè)定優(yōu)化多目標、約束條件和設(shè)計變量,通過優(yōu)化算法多次迭代后將優(yōu)化結(jié)果回代常規(guī)設(shè)計軟件進行配筋計算比較,取相對優(yōu)化結(jié)果作為最終剪力墻構(gòu)型。這種方法實際上屬于一定約束條件下的結(jié)構(gòu)準拓撲優(yōu)化,是對材料滿足約束條件下的合理分布,實現(xiàn)物盡其用,進而達到拓撲層次的優(yōu)化。在剪力墻拓撲分布優(yōu)化的前提下,對結(jié)構(gòu)構(gòu)件進行尺寸優(yōu)化。結(jié)構(gòu)敏度分析和整體優(yōu)化流程如圖1所示。
合肥置地廣場·栢悅公館8號樓位于安徽省合肥市政務(wù)文化新區(qū),項目總占地面積為51 188 m2,總建筑面積為380 000 m2,結(jié)構(gòu)三維模型如圖2所示。
圖1 結(jié)構(gòu)敏度和優(yōu)化分析流程
圖2 結(jié)構(gòu)三維模型圖
本結(jié)構(gòu)為高層住宅,塔樓地上41層,地下2層,結(jié)構(gòu)屋面高度122.75 m,采用剪力墻結(jié)構(gòu)體系。以地下室頂板面為計算嵌固端,重力荷載代表值為4.43×104t,以平地震荷載為控制荷載,結(jié)構(gòu)第一周期3.13 s,振型為Y向平動。
本工程所在場地的特征周期為0.35 s,選用彈性動力時程分析,取地震荷載時長為40 s,時間步長為0.01 s,地震加速度的峰值為 344.5 cm/s2。
計算選取地下室頂板以上部分作為彈性動力時程分析的對象,從地下室頂板處豎向構(gòu)件的端點輸入計算選取的地震波。有限元分析模型包含彈性殼單元和彈性梁單元,梁、柱用可考慮線性剪切效應(yīng)的一階三維梁單元B31模擬,墻、樓板等用4節(jié)點縮減積分曲殼單元S4R和S3模擬。
本工程實例計算條件設(shè)定如下:
多目標函數(shù):結(jié)構(gòu)的造價和第一階周期(假設(shè)鋼筋價格為3 500元/t,混凝土價格為450元/m3)。
約束條件:樓層層間位移角不大于1/1 000;墻軸壓比限值0.5;扭轉(zhuǎn)平動周期比不大于0.85。
設(shè)計變量:分組剪力墻單元的材料屬性。
優(yōu)化前后剪力墻標準層結(jié)構(gòu)平面圖對比如圖3,圖4所示。
圖3 初始設(shè)計標準層結(jié)構(gòu)平面圖
在層間位移角、軸壓比、扭轉(zhuǎn)平動位移比等約束都滿足的條件下,優(yōu)化前后結(jié)構(gòu)造價對比見圖5。結(jié)果表明,與傳統(tǒng)結(jié)構(gòu)設(shè)計方法相比,多目標優(yōu)化方法在滿足各項約束指標的前提下,能進一步優(yōu)化,實現(xiàn)了自動尋優(yōu)精細化設(shè)計,減少了建筑的材料用量,達到了建筑節(jié)能的目標。
圖4 優(yōu)化設(shè)計標準層結(jié)構(gòu)平面圖
圖5 結(jié)構(gòu)優(yōu)化前后造價對比
本文開發(fā)了一套完整的高層剪力墻結(jié)構(gòu)多目標優(yōu)化流程,可以實現(xiàn)高層剪力墻結(jié)構(gòu)的準拓撲優(yōu)化和墻尺寸優(yōu)化,并通過工程算例驗證了其實用性,初步實現(xiàn)了對高層剪力墻結(jié)構(gòu)體系和尺寸的優(yōu)化,彌補了傳統(tǒng)工程設(shè)計的不足,為實現(xiàn)自動化、智能化結(jié)構(gòu)設(shè)計提供了切實可行的思路。