潘毅 ,陳齊 ,王騰 ,周祎

（西南交通大學(xué) a.土木工程學(xué)院；b.抗震工程技術(shù)四川省重點(diǎn)實(shí)驗(yàn)室，成都 610031）

目前，結(jié)構(gòu)減震設(shè)計(jì)大多采用被動(dòng)控制技術(shù)[1]，且多采用耗能減震設(shè)計(jì)，即在結(jié)構(gòu)的合理位置布置適宜數(shù)量的阻尼器來(lái)消耗地震能量，減小結(jié)構(gòu)地震反應(yīng)。適宜的阻尼器數(shù)量、位置是減震設(shè)計(jì)的關(guān)鍵。減震設(shè)計(jì)流程大致可分為阻尼器豎向布置和阻尼器平面布置。在阻尼器豎向布置方面，將阻尼器布置在合理樓層位置，能增強(qiáng)耗能效果、提高抗震能力[2]。為實(shí)現(xiàn)減震性能和經(jīng)濟(jì)性的平衡，應(yīng)將盡量少的阻尼器布置在最優(yōu)的樓層位置，以最大限度發(fā)揮阻尼器的耗能減震性能并降低經(jīng)濟(jì)成本。但目前的減震設(shè)計(jì)布置方法既不能保證阻尼器布置是最優(yōu)方案，又要進(jìn)行繁瑣的迭代，其效率有待提高[3]。在阻尼器豎向布置方面，高效、合理地確定各樓層阻尼器數(shù)量是結(jié)構(gòu)減震設(shè)計(jì)亟須解決的問(wèn)題。在阻尼器平面布置方面，基于阻尼器豎向布置得到各樓層阻尼器的合理數(shù)量之后，需要為各樓層阻尼器選擇合理的安裝位置，確定阻尼器平面位置是一個(gè)迭代、調(diào)整和耗時(shí)耗力的過(guò)程。因此，在滿足業(yè)主使用功能與阻尼器平面布置規(guī)則的前提下，快速確定阻尼器平面安裝位置是結(jié)構(gòu)減震設(shè)計(jì)亟須解決的問(wèn)題。為解決上述問(wèn)題，眾多學(xué)者進(jìn)行了阻尼器豎向布置優(yōu)化方法和阻尼器平面布置優(yōu)化方法的探索與研究。

在阻尼器豎向布置優(yōu)化方面，因阻尼器布置與結(jié)構(gòu)地震響應(yīng)之間無(wú)明顯、有規(guī)律的顯式關(guān)系，故目前的阻尼器豎向優(yōu)化布置大多采用隨機(jī)搜索優(yōu)化布置方法和遺傳算法。例如，Agrawal 等[4]采用3種智能搜索算法，即序列、WOBI 和ESPS，研究在不同類(lèi)型目標(biāo)函數(shù)下的阻尼器最優(yōu)位置，結(jié)果表明，通過(guò)最小化的結(jié)構(gòu)最大層間位移目標(biāo)函數(shù)得到的阻尼器最佳位置取決于具體地震動(dòng)，WOBI 和ESPS 在改善序列搜索最佳位置方面有效。李鋼等[5]基于遺傳算法，提出一種應(yīng)用于耗能減震結(jié)構(gòu)參數(shù)設(shè)計(jì)的新型優(yōu)化數(shù)學(xué)模型，結(jié)果表明，新模型是一種有效的優(yōu)化設(shè)計(jì)方法。李宏男等[6]基于遺傳算法，提出阻尼器新型優(yōu)化數(shù)學(xué)模型，對(duì)比了位移型和速度型阻尼器優(yōu)化布設(shè)的差異，發(fā)現(xiàn)高層建筑宜選用速度型阻尼器。燕樂(lè)緯等[7]基于數(shù)字序列編碼遺傳算法，提出高層結(jié)構(gòu)黏滯阻尼器優(yōu)化布置方法，結(jié)果表明，該算法能加速種群進(jìn)化，在高層結(jié)構(gòu)阻尼器優(yōu)化布置上有效。馬宏偉等[8]基于粗粒度并行遺傳算法，提出通過(guò)Matlab-ABAQUS-Python 的交互使用來(lái)優(yōu)化被動(dòng)控制結(jié)構(gòu)中阻尼器位置的方法，結(jié)果表明，該算法能顯著提高結(jié)構(gòu)減震率。金波等[9]基于改進(jìn)的遺傳算法，提出以替換桿件模態(tài)應(yīng)變能百分比之和最大為適應(yīng)度函數(shù)來(lái)優(yōu)化大跨網(wǎng)架結(jié)構(gòu)上的黏滯阻尼器數(shù)量和位置的方法，結(jié)果表明，該方法能明顯改善結(jié)構(gòu)受力狀況，減震效果良好。王曙光等[10]采用遺傳算法進(jìn)行多種不同類(lèi)型阻尼器的同步優(yōu)化布置，結(jié)果表明，該算法可以在較短時(shí)間內(nèi)同時(shí)優(yōu)化不同類(lèi)型的阻尼器位置。陳豐收等[11]采用多目標(biāo)并行遺傳算法進(jìn)行減震結(jié)構(gòu)的阻尼器布置數(shù)量和布置位置的多目標(biāo)同步優(yōu)化，結(jié)果表明，該算法收斂性強(qiáng)，阻尼器主要布置在結(jié)構(gòu)中下部且中部較多。周玉嫻等[12]采用基因?qū)傩员Ａ暨z傳算法進(jìn)行阻尼器的雙目標(biāo)同步優(yōu)化布置，結(jié)果表明，該方法可高效處理雙目標(biāo)阻尼器同步優(yōu)化問(wèn)題。盡管以上研究在阻尼器豎向優(yōu)化布置方面已經(jīng)很成熟，但是較少同時(shí)考慮阻尼器的豎向布置與平面布置，即在確定各樓層的阻尼器數(shù)量后，需要考慮在不妨礙業(yè)主使用功能需求且不增大剛度與質(zhì)心差異、結(jié)構(gòu)扭轉(zhuǎn)效應(yīng)的前提下，在各樓層平面上合理布置阻尼器。

在阻尼器平面優(yōu)化布置方面，目前處于起步階段，相關(guān)研究較少。例如，Erdem ?er?evik 等[13]采用蝙蝠算法和蜻蜓算法進(jìn)行黏性墻式阻尼器數(shù)量、阻尼器周邊布置的優(yōu)化，并以樓層加速度和層間位移角為約束條件，結(jié)果表明，該方法有利于改進(jìn)結(jié)構(gòu)響應(yīng)、最小化經(jīng)濟(jì)成本。程衛(wèi)紅等[14]通過(guò)模型算例針對(duì)墻式連接位移型阻尼器進(jìn)行了平面布置分析，結(jié)果表明，阻尼器平面布置宜規(guī)則對(duì)稱(chēng)，且盡量布置在框架結(jié)構(gòu)周邊，從而減少不規(guī)則布置帶來(lái)的附加扭轉(zhuǎn)效應(yīng)。目前針對(duì)阻尼器智能選擇平面位置方面的研究較少，為實(shí)現(xiàn)和加速減震設(shè)計(jì)的智能化，有必要開(kāi)展相關(guān)研究。

針對(duì)以上不足，筆者基于雙目標(biāo)優(yōu)化算法（Dual Objective Optimization Algorithm，DOOA）和生成對(duì)抗網(wǎng)絡(luò)（Generative Adversarial Network，GAN），提出阻尼器豎向優(yōu)化布置與智能選擇阻尼器平面位置的方法。采用雙目標(biāo)優(yōu)化算法進(jìn)行阻尼器豎向布置，并與原結(jié)構(gòu)、工程師設(shè)計(jì)、逐層逼近法進(jìn)行對(duì)比分析；采用生成對(duì)抗網(wǎng)絡(luò)進(jìn)行阻尼器平面布置，并與工程師設(shè)計(jì)結(jié)果進(jìn)行相似性差異度對(duì)比分析。通過(guò)工程案例的應(yīng)用與分析，驗(yàn)證方法的有效性，為減震結(jié)構(gòu)的智能化設(shè)計(jì)提供參考。

1 阻尼器豎向布置優(yōu)化目標(biāo)

1.1 減震結(jié)構(gòu)運(yùn)動(dòng)方程

在地震作用下，未設(shè)置阻尼器的框架結(jié)構(gòu)的運(yùn)動(dòng)方程為

式中：M、C0、K0分別為未布置阻尼器的框架結(jié)構(gòu)的質(zhì)量矩陣、阻尼矩陣和剛度矩陣；(t)(t)、u(t)分別為框架結(jié)構(gòu)的樓層加速度向量、速度向量和位移向量；P(t)為阻尼器所提供的阻尼力矩陣；(t)為地面加速度，即地震動(dòng)加速度；I為單位列向量。

為兼容多種速度型阻尼器的力學(xué)模型，考慮3種常見(jiàn)的速度型阻尼器力學(xué)模型，包括線性模型、Kelvin 模型和Maxwell 模型，3 種模型均假設(shè)阻尼指數(shù)為1。線性模型的阻尼力表達(dá)式為

式中：Cd為阻尼器的阻尼矩陣；(t)為第i層樓面相對(duì)基礎(chǔ)的變形速度。將式（2）代入式（1），可得到布置線性阻尼器模型的框架結(jié)構(gòu)運(yùn)動(dòng)方程。

由式（3）可得出，線性阻尼器模型只改變了結(jié)構(gòu)的阻尼矩陣。

Kelvin 模型采用的是剛度與阻尼并聯(lián)的形式，如圖1 所示。在實(shí)際減震設(shè)計(jì)中，一般設(shè)置剛度Kd=0，以此來(lái)讓阻尼裝置在地震作用下運(yùn)動(dòng)起來(lái)，故Kelvin 模型的阻尼力矩陣表達(dá)式為

圖1 Kelvin 模型Fig.1 Kelvin model

式中：Kd為阻尼器的剛度矩陣；Cd為框架結(jié)構(gòu)的層間側(cè)移剛度矩陣；u(t)為樓面相對(duì)基礎(chǔ)的變形向量。將式（4）代入式（1），可得到布置Maxwell 阻尼器模型的框架結(jié)構(gòu)運(yùn)動(dòng)方程同式（3）。

Maxwell 模型采用的是剛度與阻尼串聯(lián)的形式，如圖2 所示。在實(shí)際減震設(shè)計(jì)中，一般設(shè)置剛度為阻尼系數(shù)的100～1 000 倍，從而讓阻尼裝置消耗地震能量，故項(xiàng)可忽略不計(jì)，Maxwell 模型的單個(gè)阻尼器阻尼力表達(dá)式為

圖2 Maxwell 模型Fig.2 Maxwell model

式中：p(t)為單個(gè)阻尼器提供的阻尼力(t)為單個(gè)阻尼器兩端相對(duì)變形速度大小。將式（5）代入式（1），可得到布置Maxwell 阻尼器模型的框架結(jié)構(gòu)運(yùn)動(dòng)方程同式（3）。

綜上所述，當(dāng)阻尼指數(shù)為1 時(shí)，3 種阻尼器模型的結(jié)構(gòu)運(yùn)動(dòng)方程在形式上統(tǒng)一，均為式（6）形式。

1.2 減震結(jié)構(gòu)響應(yīng)求解

在布置阻尼器后，框架結(jié)構(gòu)的阻尼體系變?yōu)榉墙?jīng)典阻尼體系，即非對(duì)角矩陣，運(yùn)動(dòng)方程組無(wú)法解耦為N個(gè)獨(dú)立的單自由度體系運(yùn)動(dòng)方程。為解決此類(lèi)非經(jīng)典阻尼解耦問(wèn)題，采用復(fù)振型分解法進(jìn)行方程組解耦，由線性定常動(dòng)力體系的微分方程經(jīng)變量代換得到減震結(jié)構(gòu)體系的狀態(tài)方程[15]。具體表達(dá)式為

因矩陣D為非對(duì)稱(chēng)矩陣，故應(yīng)采用廣義模態(tài)分析法進(jìn)行解耦，通過(guò)特征向量的正交性計(jì)算框架結(jié)構(gòu)的地震響應(yīng)。假設(shè)非對(duì)稱(chēng)矩陣D的左、右特征向量與特征值分別為θi、φi、λi，其中i=1,2…2N。

由于振型的正交特性，框架結(jié)構(gòu)的任何響應(yīng)均可采用振型的線性組合進(jìn)行表示。例如，采用右特征向量進(jìn)行振型的線性組合，如式（8）所示。

式中：φ=[φ1,φ2…φi…φ2N]，η(t)為單模態(tài)坐標(biāo)向量。

將式（8）代入式（7）可得

式（9）兩側(cè)同時(shí)左乘進(jìn)行解耦

在確定的地震激勵(lì)作用下，減震結(jié)構(gòu)單模態(tài)響應(yīng)為

因此，框架結(jié)構(gòu)的樓層位移響應(yīng)向量可表示為

框架結(jié)構(gòu)的樓層加速度響應(yīng)向量可表示為

1.3 豎向布置優(yōu)化目標(biāo)函數(shù)

在布置阻尼器后，一般需減小框架結(jié)構(gòu)的樓層加速度峰值與層間位移角峰值，故采用多目標(biāo)優(yōu)化方法，優(yōu)化目標(biāo)函數(shù)為

2 阻尼器平面布置評(píng)價(jià)指標(biāo)

2.1 布置規(guī)則

基于框架結(jié)構(gòu)的平面建筑圖，進(jìn)行阻尼器的平面布置。在已確定框架結(jié)構(gòu)阻尼器豎向布置的基礎(chǔ)上，對(duì)需要布置阻尼器的樓層進(jìn)行智能化平面布置。依據(jù)《全國(guó)民用建筑工程設(shè)計(jì)技術(shù)措施——結(jié)構(gòu)（結(jié)構(gòu)體系）》規(guī)范附錄C.4，阻尼器平面布置應(yīng)符合表1 所示的規(guī)則。

表1 平面布置規(guī)則Table 1 Rules of plane placement

2.2 評(píng)價(jià)指標(biāo)

基于阻尼器平面布置規(guī)則，以較優(yōu)的目標(biāo)設(shè)計(jì)為設(shè)計(jì)標(biāo)準(zhǔn)，提出用于評(píng)價(jià)生成設(shè)計(jì)的相似性差異度綜合評(píng)價(jià)指標(biāo)E，以及兩個(gè)子指標(biāo)：阻尼器質(zhì)心差異度、阻尼器分散性差異度。質(zhì)心差異度可以反映阻尼器質(zhì)心與結(jié)構(gòu)質(zhì)量中心之間的差異程度，該值越小則阻尼器布置越對(duì)稱(chēng)，其附加扭轉(zhuǎn)越小。分散性差異度可以反映阻尼器是否分散布置，該值越小則結(jié)構(gòu)抗扭轉(zhuǎn)能力越強(qiáng)。將生成的平面布置與目標(biāo)設(shè)計(jì)的平面布置進(jìn)行逐樓層、逐柱跨對(duì)比，通過(guò)式（15）～式（17）計(jì)算兩者的相似性。E值越小，說(shuō)明生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)之間的差異越小，相似程度越高。當(dāng)相似性差異度綜合評(píng)價(jià)指標(biāo)E小于臨界值α?xí)r，則生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)相似性極高，結(jié)構(gòu)的抗扭能力較好。

式中：E為生成設(shè)計(jì)的相似性綜合評(píng)價(jià)指標(biāo)；為生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)在阻尼器質(zhì)心上的距離大小；為生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)在第j層樓的阻尼器分散程度上的差異大?。籲為框架結(jié)構(gòu)樓層數(shù)量；xg、yg分別為生成設(shè)計(jì)第j層樓X向、Y向阻尼器質(zhì)心與結(jié)構(gòu)質(zhì)心x0、y0之間的距離；xe、ye分別為目標(biāo)設(shè)計(jì)第j層樓X向、Y向阻尼器質(zhì)心與結(jié)構(gòu)質(zhì)心x0、y0之間的距離；L、B分別為框架結(jié)構(gòu)的長(zhǎng)度與寬度；xgi、xei分別為生成設(shè)計(jì)、目標(biāo)設(shè)計(jì)第j層樓第i個(gè)阻尼器與結(jié)構(gòu)質(zhì)量中心x0之間的距離；ygi、yei分別為生成設(shè)計(jì)、目標(biāo)設(shè)計(jì)第j層樓第i個(gè)阻尼器與建筑質(zhì)量中心y0之間的距離；n1、n2分別為生成設(shè)計(jì)、目標(biāo)設(shè)計(jì)的X向與Y向阻尼器數(shù)量之和。

需要說(shuō)明的是，生成設(shè)計(jì)是基于提取的可布置阻尼器位置圖像進(jìn)行的圖像生成，如圖3 所示，即需要對(duì)圖紙做前處理工作。

圖3 生成設(shè)計(jì)圖像Fig.3 The figure of generate design

3 阻尼器布置智能設(shè)計(jì)算法

3.1 阻尼器豎向智能布置

采用雙目標(biāo)優(yōu)化算法來(lái)實(shí)現(xiàn)阻尼器的豎向布置。雙目標(biāo)優(yōu)化算法是借鑒了麻雀搜索算法的優(yōu)化思想,通過(guò)模擬麻雀搜索食物的行為來(lái)尋優(yōu)，從而尋找最優(yōu)解。在一個(gè)D維搜索空間內(nèi)，存在pop只個(gè)體，第i只個(gè)體在搜索空間中的位置可表示為Xi=[xi1…xij…xiD]，i=1,2…N，其中xij為第i只個(gè)體在第j維空間的位置?；陔p目標(biāo)優(yōu)化算法的阻尼器豎向布置流程圖見(jiàn)圖4，具體步驟如下。

圖4 阻尼器豎向智能布置流程圖Fig.4 Flow chart of damper vertical intelligent placement

Step1：初始化數(shù)據(jù)，種群規(guī)模pop，優(yōu)化目標(biāo)函數(shù)維數(shù)D，迭代次數(shù)T，初始值上下界lb、ub。

Step2：初始化種群，生成pop 個(gè)D維向量Xi，維度對(duì)應(yīng)樓層數(shù)量。

Step3：計(jì)算個(gè)體的目標(biāo)函數(shù)值f1和f2，選出優(yōu)化目標(biāo)函數(shù)f1和f2均相對(duì)較小時(shí)所對(duì)應(yīng)的個(gè)體作為最優(yōu)解。

Step4：循環(huán)執(zhí)行step2、step3，迭代T次，每次迭代后更新阻尼器豎向布置最優(yōu)解。

Step5：當(dāng)?shù)螖?shù)t=T，程序結(jié)束，輸出最后一次迭代所得到的阻尼器豎向布置最優(yōu)解。

3.2 阻尼器平面智能布置

采用GAN 來(lái)實(shí)現(xiàn)阻尼器的平面布置。GAN是自動(dòng)化建筑設(shè)計(jì)中應(yīng)用最廣泛的框架之一[16]，是一種深度生成模型，該模型框架包含兩個(gè)相互競(jìng)爭(zhēng)的神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)：判別器D 和生成器G。生成器用來(lái)捕獲數(shù)據(jù)分布并生成樣本，而判別器用于將生成的樣本與真實(shí)樣本區(qū)分開(kāi)來(lái)[17]，二者通過(guò)對(duì)抗學(xué)習(xí)提升各自的能力。生成對(duì)抗網(wǎng)絡(luò)可通過(guò)對(duì)既有設(shè)計(jì)圖像或文本的學(xué)習(xí)，掌握已有設(shè)計(jì)數(shù)據(jù)中的潛在設(shè)計(jì)規(guī)律，實(shí)現(xiàn)新設(shè)計(jì)圖像的生成[18]?；贕AN 的阻尼器平面布置流程圖見(jiàn)圖5，算法具體步驟如下。

圖5 阻尼器平面智能布置流程圖Fig.5 Flow chart of damper plane intelligent placement

Step1：將用于訓(xùn)練、驗(yàn)證和測(cè)試的建筑平面圖進(jìn)行前處理，即用藍(lán)色標(biāo)出所有能布置阻尼器的位置，用紅色標(biāo)出已布置阻尼器的位置。

Step2：將建筑平面圖紙的所有合理阻尼器布置用于訓(xùn)練生成器，在與判別器博弈過(guò)程中，不斷學(xué)習(xí)阻尼器布置規(guī)律。當(dāng)達(dá)到納什平衡時(shí)，生成器訓(xùn)練完成，可用于阻尼器平面布置。

Step3：將訓(xùn)練好的生成器用于生成與工程師設(shè)計(jì)相似的阻尼器平面布置，從而得到平面布置的生成設(shè)計(jì)。

Step4：采用相似性差異度綜合評(píng)價(jià)E進(jìn)行生成設(shè)計(jì)評(píng)價(jià)。若E≥α，返回到step3 中重新進(jìn)行訓(xùn)練和生成，重新評(píng)價(jià)，直至E＜α。

4 阻尼器布置智能設(shè)計(jì)工程應(yīng)用

為研究雙目標(biāo)優(yōu)化算法和生成對(duì)抗網(wǎng)絡(luò)的可行性和通用性，采用2 個(gè)實(shí)際減震設(shè)計(jì)工程案例進(jìn)行阻尼器布置設(shè)計(jì)，并開(kāi)展性能分析、對(duì)比?；谧枘崞髫Q向布置雙目標(biāo)優(yōu)化算法和阻尼器平面布置生成對(duì)抗網(wǎng)絡(luò)算法，對(duì)地震動(dòng)作用下的框架結(jié)構(gòu)進(jìn)行阻尼器豎向布置與平面布置，在豎向布置上與逐層逼近設(shè)計(jì)、工程師設(shè)計(jì)進(jìn)行對(duì)比，在平面布置上與工程師平面設(shè)計(jì)進(jìn)行對(duì)比，從而驗(yàn)證阻尼器布置方法的有效性。

4.1 工程案例一

4.1.1 結(jié)構(gòu)模型

某5 層鋼筋混凝土框架結(jié)構(gòu)教學(xué)樓，首層層高4.9 m，標(biāo)準(zhǔn)層層高4 m，總高20.9 m；柱截面為600 mm×600 mm，邊主梁截面為500 mm×800 mm，內(nèi)主梁截面為400 mm×800 mm，次梁均為300 mm×600 mm，具體結(jié)構(gòu)平面布置如圖6 所示，各樓層的質(zhì)量和剛度見(jiàn)表2。抗震設(shè)防烈度為8 度（0.2g），地震分組為第2 組，場(chǎng)地類(lèi)別為Ⅱ類(lèi)，設(shè)防類(lèi)別為乙類(lèi)。

表2 某5 層RC 框架結(jié)構(gòu)的質(zhì)量和剛度Table 2 The mass and stiffness of a five story reinforced concrete frame structure

圖6 案例一框架結(jié)構(gòu)平面布置Fig.6 Plane placement of frame structure of case one

按照工程師試算，案例一結(jié)構(gòu)需安裝30 個(gè)黏滯阻尼器，具體阻尼器參數(shù)見(jiàn)表3。優(yōu)化并求解阻尼器合理位置，提高框架結(jié)構(gòu)在地震波作用下的抗震性能，使其層間位移角和樓層加速度滿足一定的減震目標(biāo)，具體減震指標(biāo)為：在8 度罕遇地震作用下，層間位移角控制在1/200 以下，樓層加速度控制在5 m/s2以下。

表3 阻尼器參數(shù)Table 3 Parameters of dampers

4.1.2 阻尼器豎向優(yōu)化布置

依據(jù)場(chǎng)地特征周期Tg1=0.45 s、結(jié)構(gòu)基本周期Ta1=1.25 s，選取與其周期相近的3 條地震動(dòng)：SAN FERNANDO（臺(tái)站110）波、COALINGA（臺(tái)站1624）波、NORTHRIDGE（臺(tái)站5366）波，分別記為A 波、B波、C 波，并將地震波的峰值加速度調(diào)整至400 cm/s2。

在不同地震波作用下，雙目標(biāo)優(yōu)化算法能保持目標(biāo)函數(shù)值f1、f2變動(dòng)幅度較小，且在同一方向上的阻尼器布置近似相同，反映了該算法對(duì)不同地震動(dòng)的魯棒性能較好。同一樓層依據(jù)出現(xiàn)次數(shù)較多的阻尼器數(shù)目、接近平均數(shù)和阻尼器總數(shù)量為15 的原則，雙目標(biāo)優(yōu)化算法在X向與Y向的最終阻尼器豎向布置如表4 所示。相比于逐層逼近法，雙目標(biāo)優(yōu)化算法的目標(biāo)函數(shù)值與其幾乎相同，側(cè)面說(shuō)明該算法在阻尼器豎向布置上的優(yōu)化能力較好。相比于工程師設(shè)計(jì)，雙目標(biāo)優(yōu)化算法的目標(biāo)函數(shù)值f1、f2均小于工程師設(shè)計(jì)，說(shuō)明該優(yōu)化算法在阻尼器數(shù)量固定的前提下能輔助工程師進(jìn)行豎向布置修正，充分發(fā)揮阻尼器的減震性能，減小各樓層的最大加速度與最大層間位移角。

表4 案例一的阻尼器布置方案與目標(biāo)函數(shù)值Table 4 Damper layout schemes and objective function values of case one

與原結(jié)構(gòu)相比，采用不同阻尼器豎向布置方法得到的結(jié)構(gòu)最大層間位移角和最大樓層加速度均有所降低。對(duì)于最大層間位移角，雙目標(biāo)優(yōu)化算法與逐層逼近法的豎向布置均能滿足減震目標(biāo)1/200的要求，而工程師設(shè)計(jì)不能滿足減震目標(biāo)，如圖7 所示。主要原因在于工程師設(shè)計(jì)的首層與頂層未布置阻尼器，導(dǎo)致其首層、頂層成為樓層薄弱處，層間位移角未能減小至減震目標(biāo)。對(duì)于最大樓層加速度，雙目標(biāo)優(yōu)化算法與逐層逼近法的最大樓層加速度均小于減震目標(biāo)5 m/s2，但工程師設(shè)計(jì)未能達(dá)標(biāo)，如圖8 所示。這是因?yàn)槭讓?、頂層為薄弱層，?dǎo)致首層與頂層的最大樓層加速度未減小，甚至增大。

圖7 最大層間位移角包絡(luò)圖Fig.7 Envelope diagram of maximum inter-story displacement angle

圖8 最大樓層加速度包絡(luò)圖Fig.8 Envelope diagram of maximum floor acceleration

4.1.3 阻尼器平面優(yōu)化布置

在雙目標(biāo)優(yōu)化算法確定的案例一阻尼器豎向布置的基礎(chǔ)上，分別讓訓(xùn)練好的生成對(duì)抗網(wǎng)絡(luò)模型和工程師進(jìn)行阻尼器平面布置，并以工程師設(shè)計(jì)的阻尼器平面布置為目標(biāo)設(shè)計(jì)。案例一的目標(biāo)設(shè)計(jì)與生成設(shè)計(jì)如圖9、圖10 所示。

圖9 案例一目標(biāo)設(shè)計(jì)Fig.9 Objective design of case one

圖10 案例一生成設(shè)計(jì)Fig.10 Generative design of case one

在該工程案例中，生成設(shè)計(jì)的相似性綜合評(píng)價(jià)值為E1=0.05，各樓層的質(zhì)心差異度和分散性差異度見(jiàn)圖10。由于相似性綜合評(píng)價(jià)E比較接近0，生成設(shè)計(jì)與工程師設(shè)計(jì)具有較高的相似度。

4.2 工程案例二

4.2.1 結(jié)構(gòu)模型

某6 層鋼筋混凝土框架結(jié)構(gòu)辦公樓，首層層高6 m，標(biāo)準(zhǔn)層層高4 m，總高26 m；柱截面為650 mm×650 mm，邊主梁截面為600 mm ×800 mm，內(nèi)主梁截面為400 mm×700 mm，次梁均為300 mm ×600 mm，結(jié)構(gòu)平面布置如圖11所示，各樓層的質(zhì)量和剛度見(jiàn)表5。抗震設(shè)防烈度為8度（0.3g），地震分組為第3組，場(chǎng)地類(lèi)別為Ⅱ類(lèi)，設(shè)防類(lèi)別為乙類(lèi)。

表5 某6 層RC 框架結(jié)構(gòu)的質(zhì)量和剛度Table 5 The mass and stiffness of a six story reinforced concrete frame structure

圖11 案例二框架結(jié)構(gòu)平面布置Fig.11 Plane placement of frame structure for case two

依據(jù)工程師試算，案例二需安裝48 個(gè)黏滯阻尼器，具體阻尼器參數(shù)見(jiàn)表6。在工程師設(shè)計(jì)的基礎(chǔ)上，優(yōu)化阻尼器布置位置，提升框架結(jié)構(gòu)的抗震性能，使其層間位移角和樓層加速度滿足一定的減震目標(biāo)，具體減震指標(biāo)為：在8 度罕遇地震作用下，層間位移角控制在1/200 以下，樓層加速度控制在6 m/s2以下。

表6 阻尼器參數(shù)Table 6 The parameters of dampers

4.2.2 阻尼器豎向優(yōu)化布置

根據(jù)案例二場(chǎng)地特征周期Tg2=0.5 s 與結(jié)構(gòu)基本周期Ta2=1.27 s，選取出與其周期相近的3 條地震 動(dòng)：SAN FERNANDO（臺(tái) 站110）波、SAN FERNANDO（臺(tái)站160）波、NORTHRIDGE（臺(tái)站5366）波，分別記為A 波、D 波、C 波，并將各地震波的峰值加速度調(diào)整至510 cm/s2。

案例二結(jié)構(gòu)在X向與Y向的質(zhì)量、剛度和平面布置完全相同，故在各樓層兩主軸方向上布置同等數(shù)量的阻尼器。同一樓層依據(jù)出現(xiàn)次數(shù)較多的阻尼器數(shù)目、接近平均數(shù)和阻尼器總數(shù)量為24 的原則，雙目標(biāo)優(yōu)化算法在X向與Y向的最終阻尼器豎向布置如表7 所示。在3 條地震動(dòng)作用下，雙目標(biāo)優(yōu)化算法的目標(biāo)函數(shù)值與逐層逼近法、工程師設(shè)計(jì)非常接近，但其f1值均優(yōu)于逐層逼近法，f2值均優(yōu)于工程師設(shè)計(jì)，再次說(shuō)明該算法在阻尼器豎向布置上的有效性、合理性。而且，布置的阻尼器集中于中下部樓層，符合以往的阻尼器布置經(jīng)驗(yàn)。

表7 案例二的阻尼器布置方案和目標(biāo)函數(shù)值Table 7 Damper placement schemes and objective function values of case two

與原結(jié)構(gòu)相比，不同阻尼器布置方法的最大層間位移角和最大樓層加速度均有所降低。對(duì)于最大層間位移角，雙目標(biāo)優(yōu)化算法與逐層逼近法的豎向布置均能滿足減震目標(biāo)1/200 的要求，而在A 波作用下，工程師設(shè)計(jì)的1、2 層的層間位移角未能滿足減震目標(biāo)，如圖12 所示。主要原因在于首層的樓層剪力較大和阻尼器數(shù)量不足，導(dǎo)致首層的層間變形過(guò)大，使其層間位移角未能滿足減震目標(biāo)。對(duì)于最大樓層加速度，雙目標(biāo)優(yōu)化算法、逐層逼近法和工程師設(shè)計(jì)的最大樓層加速度均小于加速度減震目標(biāo)6 m/s2，如圖13 所示。

圖12 最大層間位移角包絡(luò)圖Fig.12 Envelope diagram of maximum inter-story displacement angle

圖13 最大樓層加速度包絡(luò)圖Fig.13 Envelope diagram of maximum floor acceleration

4.2.3 阻尼器平面優(yōu)化布置

在雙目標(biāo)優(yōu)化算法確定的案例二阻尼器豎向布置的基礎(chǔ)上，分別讓訓(xùn)練好的GAN 模型和工程師進(jìn)行阻尼器平面布置，并以工程師設(shè)計(jì)的阻尼器平面布置為目標(biāo)設(shè)計(jì)。案例二的目標(biāo)設(shè)計(jì)與生成設(shè)計(jì)如圖14、圖15 所示。

圖14 案例二目標(biāo)設(shè)計(jì)Fig.14 Objective design of case two

圖15 案例二生成設(shè)計(jì)Fig.15 Generative design of case two

在該工程案例中，生成設(shè)計(jì)的綜合相似性評(píng)價(jià)值為E=0.02，各樓層的質(zhì)心差異度和分散性差異度見(jiàn)圖15。由于相似性綜合評(píng)價(jià)E非常接近0，生成設(shè)計(jì)與工程師設(shè)計(jì)具有較高的相似度。雖然生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)的阻尼器平面位置不完全相同，但阻尼器質(zhì)心與分散性相同，導(dǎo)致其實(shí)際減震效果完全相同，例如該案例的5F、6F。

4.3 雙目標(biāo)優(yōu)化算法效率分析

阻尼器智能布置在豎向布置上耗時(shí)占比最大，對(duì)算法效率影響最大，故只針對(duì)雙目標(biāo)優(yōu)化算法進(jìn)行效率分析。

在EI Centro 地震波作用下，以案例二的標(biāo)準(zhǔn)層平面布置、阻尼器參數(shù)為依據(jù)，開(kāi)展2、4、6、8、10 和12 層框架結(jié)構(gòu)的雙目標(biāo)優(yōu)化算法效率分析，并與逐層逼近法、工程師設(shè)計(jì)進(jìn)行對(duì)比。采用雙目標(biāo)優(yōu)化算法、逐層逼近法和工程師設(shè)計(jì)對(duì)框架結(jié)構(gòu)分別進(jìn)行5 次設(shè)計(jì)，取其所用平均時(shí)間，其統(tǒng)計(jì)結(jié)果如表8所示。

表8 計(jì)算時(shí)間Table8 The time of calculation

由表8 可知，雙目標(biāo)優(yōu)化算法在阻尼器豎向布置上所用時(shí)間隨樓層數(shù)量與阻尼器數(shù)量增大而增長(zhǎng)，且遠(yuǎn)小于逐層逼近法與工程師設(shè)計(jì)。雙目標(biāo)優(yōu)化算法計(jì)算時(shí)間較短，分析原因主要有兩點(diǎn)：1）算法內(nèi)嵌入了樓層地震響應(yīng)有限元程序，可快速計(jì)算最大樓層加速度、最大層間位移角，避免了調(diào)用有限元軟件計(jì)算速度慢的問(wèn)題；2）調(diào)整雙目標(biāo)優(yōu)化算法的種群數(shù)量為5，迭代次數(shù)為100，既能保證優(yōu)化結(jié)果較好，又能極大降低優(yōu)化計(jì)算所用時(shí)間。而逐層逼近法計(jì)算時(shí)間最長(zhǎng)，主要原因在于需要工程師使用有限元軟件迭代計(jì)算并確定每個(gè)阻尼器所在樓層，直至布置完所有阻尼器。工程師設(shè)計(jì)所用時(shí)間位于兩者之間，是由于工程師可以利用積累的減震設(shè)計(jì)經(jīng)驗(yàn)減少迭代次數(shù)，從而較快地完成阻尼器豎向布置。

4.4 平面布置評(píng)價(jià)指標(biāo)臨界值α

目前，尚未明確平面評(píng)價(jià)指標(biāo)E的臨界值，即：當(dāng)E小于某個(gè)臨界值時(shí)，生成設(shè)計(jì)與目標(biāo)設(shè)計(jì)相似度較高，平面布置可靠，結(jié)構(gòu)力學(xué)性能較好。因此，以抗扭力學(xué)性能的優(yōu)劣來(lái)搜尋平面布置評(píng)價(jià)指標(biāo)臨界值α。針對(duì)阻尼器平面布置的相似性綜合評(píng)價(jià)指標(biāo)E的臨界值，以工程案例二為試算模型，基于工程師設(shè)計(jì)的阻尼器豎向布置，開(kāi)展12 種不同平面布置的扭轉(zhuǎn)力學(xué)分析。值得說(shuō)明的是，雖然工程案例采用的黏滯阻尼器在理論上不提供剛度，但在實(shí)際工程中存在一定的剛度，會(huì)對(duì)結(jié)構(gòu)的扭轉(zhuǎn)產(chǎn)生影響。經(jīng)過(guò)試算得到不同平面布置方式下的框架結(jié)構(gòu)扭轉(zhuǎn)周期比和考慮偶然偏心作用下的X、Y方向最大位移比，以及相似性綜合評(píng)價(jià)E值，如表9 所示。其中，考慮偶然偏心的位移比采用規(guī)范反應(yīng)譜方法進(jìn)行計(jì)算，偶然偏心取0.05[19]。已知工程師設(shè)計(jì)的扭轉(zhuǎn)周期比為0.785 9，最大扭轉(zhuǎn)位移比均值為1.173。

表9 結(jié)構(gòu)扭轉(zhuǎn)周期比、位移比和E 值Table 9 The torsion period ratio,displacement ratio and E value of the structure

將12 種平面布置方式的分析結(jié)果繪制在圖16中，可以看出：當(dāng)相似性差異度綜合評(píng)價(jià)E＜α=0.1時(shí)，可保證生成設(shè)計(jì)與工程師設(shè)計(jì)的扭轉(zhuǎn)周期比差異在10%以內(nèi)，最大扭轉(zhuǎn)位移比差異在5%以內(nèi)，且扭轉(zhuǎn)周期比、扭轉(zhuǎn)位移比分別滿足規(guī)范要求的0.9、1.6。案例一、案例二的相似性差異度綜合指標(biāo)E均小于0.1，間接說(shuō)明生成對(duì)抗網(wǎng)絡(luò)生成的平面布置滿足結(jié)構(gòu)抗扭力學(xué)性能要求。

圖16 力學(xué)性能與評(píng)價(jià)指標(biāo)相關(guān)性Fig.16 The correlation between mechanical properties and evaluation indexes

5 結(jié)論

基于雙目標(biāo)優(yōu)化算法和生成對(duì)抗網(wǎng)絡(luò)算法進(jìn)行阻尼器豎向和平面布置的智能設(shè)計(jì)，并進(jìn)行了工程應(yīng)用，得到以下結(jié)論。

1）基于結(jié)構(gòu)地震響應(yīng)公式與雙目標(biāo)優(yōu)化算法進(jìn)行阻尼器豎向智能布置。在推導(dǎo)出線性黏滯阻尼減震結(jié)構(gòu)彈性階段地震響應(yīng)公式后，采用雙目標(biāo)優(yōu)化算法求解優(yōu)化目標(biāo)函數(shù)值f1、f2最小時(shí)的各樓層阻尼器數(shù)量，可實(shí)現(xiàn)線性黏滯阻尼器的豎向智能布置。

2）基于生成對(duì)抗網(wǎng)絡(luò)進(jìn)行阻尼器水平智能布置。在完成阻尼器豎向布置后，即在確定各樓層阻尼器數(shù)量后，采用生成對(duì)抗網(wǎng)絡(luò)算法，對(duì)其生成器進(jìn)行訓(xùn)練，利用生成器完成阻尼器平面智能布置，可實(shí)現(xiàn)線性黏滯阻尼器的平面智能布置。

3）對(duì)兩個(gè)框架結(jié)構(gòu)減震設(shè)計(jì)案例進(jìn)行雙目標(biāo)優(yōu)化算法與生成對(duì)抗網(wǎng)絡(luò)的工程應(yīng)用。通過(guò)在豎向布置上同逐層逼近法、工程師設(shè)計(jì)進(jìn)行對(duì)比，在平面布置上同工程師平面布置進(jìn)行對(duì)比，可知：雙目標(biāo)優(yōu)化算法生成的阻尼器豎向布置方案滿足減震目標(biāo)、生成對(duì)抗網(wǎng)絡(luò)算法生成的阻尼器平面布置方案均與工程師的布置方案基本一致，直接說(shuō)明了所提智能設(shè)計(jì)方法在阻尼器布置上的有效性。