基于遺傳算法的臨水加筋土擋墻擬靜力內(nèi)部穩(wěn)定分析

侯超群，檀 奕，阮曉波，陳 歡

(合肥工業(yè)大學(xué) 汽車與交通工程學(xué)院，安徽 合肥 230009)

0 引言

加筋土擋墻由于具有經(jīng)濟(jì)、安全以及施工方便等優(yōu)點(diǎn)，在水利、交通、城建等領(lǐng)域得到了廣泛的工程應(yīng)用。在地震荷載的作用下，相比傳統(tǒng)的擋土結(jié)構(gòu)，加筋土擋墻能表現(xiàn)出良好的抗震性能[1]。臨水邊坡由于受到水位的影響，土體強(qiáng)度存在較大的變異性，在地震作用比一般邊坡更易發(fā)生失穩(wěn)，因此在地震頻發(fā)地區(qū)采用加筋土擋墻作為臨水護(hù)岸結(jié)構(gòu)，可以利用其抗震性能，達(dá)到工程應(yīng)用的安全性。

地震作用下加筋土擋墻的內(nèi)部穩(wěn)定分析是其工程設(shè)計中的必要環(huán)節(jié)。地震作用下加筋土擋墻的計算方法有：擬靜力法、擬動力法以及數(shù)值分析法等[2-5]。擬靜力法是將地震慣性力視為自重與地震加速度系數(shù)的乘積，并將其施于墻后不穩(wěn)定體上進(jìn)行極限平衡分析[6]。擬動力法考慮地震頻率等因素，基于彈性波理論來求解土體中的剪切波和壓縮波的波速，從而得到水平和垂直

加速度沿墻高的分布方程，并通過建立平衡方程得到地震力[7]。數(shù)值模擬法用計算機(jī)模擬地震作用下土體本身強(qiáng)度變化和破壞過程，具有易于建模以及能夠反復(fù)操作的優(yōu)點(diǎn)[8]。當(dāng)前的設(shè)計規(guī)范一般簡化滑裂面形式，基于極限平衡原理采用擬靜力法進(jìn)行加筋土擋墻的抗震設(shè)計。Choudhury和Ahmad[9]、黃睿等[10]基于擬靜力法和水平層分析法假定滑裂面為直線型和折線型對臨水加筋土擋墻進(jìn)行內(nèi)部穩(wěn)定分析。Terzaghi[11-12]研究表明，假設(shè)滑動面為直線和折線型，所計算得到的筋材拉力系數(shù)會偏大，計算誤差有時可達(dá)到2～3倍，甚至更大。通常，擋墻計算中簡化破裂面為直線型和折線型等形式，但采用對數(shù)螺旋線破壞機(jī)制更加符合實(shí)際。

在已知破裂面形式的基礎(chǔ)上，建立平衡方程，求得結(jié)構(gòu)破壞的最危險滑裂面，也就意味著獲得了安全系數(shù)(對于邊坡穩(wěn)定分析)、土壓力(對于擋土墻設(shè)計)或者所需筋材拉力(對于加筋土擋墻內(nèi)部穩(wěn)定分析)的極值。采用智能搜索方法是確定最危險滑裂面的一種重要途徑。人工智能算法搜索效率高、不易陷入局部極值，并且易與其他算法相結(jié)合。目前在臨界滑動面搜索中的人工智能方法有遺傳算法、模擬退火算法和仿生算法等[13]。其中，遺傳算法通用性強(qiáng)、全局搜索能力強(qiáng)，且不需要梯度信息、不要求函數(shù)連續(xù)，適合于并行處理和大型復(fù)雜優(yōu)化問題的求解。

本研究假定破裂面是對數(shù)螺旋線，基于極限平衡原理和擬靜力方法，建立地震荷載作用下臨水加筋土擋墻內(nèi)部穩(wěn)定分析的力矩平衡方程；進(jìn)而采用遺傳算法確定臨水加筋土擋墻的最危險滑裂面，獲得其內(nèi)部穩(wěn)定設(shè)計的所需筋材拉力；最后分析地震荷載作用下超靜孔壓、水位高度、填土內(nèi)摩擦角等因素對臨水加筋土擋墻筋材拉力的影響。

1 基本力學(xué)模型

1.1 基本假定

(1)墻后填土為各向同性的均質(zhì)無黏性土。

(2)假設(shè)破裂面是對數(shù)螺旋線狀。

(3)假設(shè)破裂面通過墻趾。

(4)已知沿著擋墻高度的筋材拉力的分布函數(shù)。

(5)墻后填土具有較高的滲透性，水可以在土體中穩(wěn)定滲流，在發(fā)生地震時填土側(cè)產(chǎn)生動水壓力。

圖1 臨水加筋土擋墻力學(xué)模型Fig.1 Mechanical model of waterfront reinforced retaining wall

1.2 受力分析

臨水加筋土擋墻的力學(xué)模型如圖1所示。墻后填土處于浸水狀態(tài)，滑動土楔體受到的力有：水平地震慣性力Qh，豎直地震慣性力Qv，土體自身重力W，以及地震引起的動水壓力Pd1，Pd2。

1.2.1 地震慣性力

本研究運(yùn)用Mononobe-Okabe法[14]計算地震慣性力，用土體的自重和地震加速度系數(shù)的乘積表示滑動土楔體受到的地震慣性力。地震慣性力Qh，Qv包括水位線以上的干土受到的地震慣性力Qhd，Qvd和水位線以下的濕土受到的地震慣性力Qhs，Qvs。

(1)

式中，kh，kv分別是水平、豎直地震慣性力系數(shù)；Wd為干土質(zhì)量。

根據(jù)Matsuzawa.et al[15]的研究，在土的滲透性較強(qiáng)時，水平地震慣性力只考慮作用于填土顆粒上的部分。因此，濕土部分的水平地震慣性力為Qhs=khWd；豎向地震慣性力為Qvs=kvW*，其中W*為用于計算浸水部分豎向地震慣性力的填土重度[2]。

W*=(γsat-γw)(1-ru)Ar，

(2)

式中，γsat為填土的飽和重度；γw為水重度；Ar為滑動土楔體浸水部分的面積；ru為地震超靜孔壓與豎向有效應(yīng)力的比值。

1.2.2 水壓力

浸水部分的填土受到的水壓力有：靜水壓力和動水壓力。由于本研究假設(shè)填土兩側(cè)水位高度相同，故靜水壓力相互抵消，不加以考慮。

根據(jù)Westergaard法[16]，作用于滑動土楔體臨水側(cè)與填土側(cè)的動水壓力為：

(3)

式中hw為水位高度。

1.3 力矩平衡方程

Leshchinsky和San在前人的基礎(chǔ)上，通過變分極限平衡法研究了土質(zhì)邊坡擬靜力條件下的對數(shù)螺旋線破壞機(jī)制及穩(wěn)定性設(shè)計[17]。本研究采用該破壞機(jī)制建立臨水加筋土擋墻內(nèi)部穩(wěn)定受力計算模型，如圖1所示，其中對數(shù)螺旋線半徑為：

R=Ae-ψmβ，

(4)

通過幾何關(guān)系，可以確定滑裂面在直角坐標(biāo)系中的方程，如下：

(5)

式中，A=H/(e-ψm β1cosβ1-e-ψm β2cosβ2)，是對數(shù)螺旋線常數(shù)；β1，β2分別為滑裂面上E，F(xiàn)點(diǎn)對應(yīng)的極角；ψm=tanφ/FS，F(xiàn)S=1.0，為安全系數(shù)；φ為墻后填土內(nèi)摩擦角；β為對應(yīng)計算點(diǎn)的極角；(xC,yC)是對數(shù)螺旋線極點(diǎn)C點(diǎn)所對應(yīng)的直角坐標(biāo)。

(6)

滑動土楔體在極限平衡狀態(tài)下，外力繞極點(diǎn)C產(chǎn)生的力矩平衡，即：

M∑Ti=MG+Mkh-Mkv+MPd1+MPd2，

(7)

式中，MG，Mkh，Mkv，MPd1及MPd2分別為滑動土楔體受到的重力、水平地震慣性力、豎直地震慣性力、臨水側(cè)動水壓力及填土側(cè)動水壓力繞C點(diǎn)的力矩。

(1)重力力矩：滑動土楔體重力繞(xC,yC)點(diǎn)的力矩分為干土和濕土兩部分進(jìn)行積分；力臂為Ae-ψmβsinβ。求得滑動土楔體重力力矩MG見附錄。

(2)地震慣性力力矩：水平和豎直地震慣性力繞(xC,yC)點(diǎn)的力臂分別為：

(8)

滑動土楔體水平、豎直地震慣性力力矩Mkh，Mkv見附錄。

(3)動水壓力力矩

(9)

臨水側(cè)動水壓力繞(xC,yC)點(diǎn)力矩MPd1如上式所示，動水壓力合力作用在0.4hw處[18]。

填土側(cè)動水壓力繞(xC,yC)點(diǎn)力矩MPd2由填土側(cè)動水壓力水平分力和豎直分力積分求得，

MPd2=MPdx+MPdy，

(10)

式中，MPdx，MPdy分別為填土側(cè)動水壓力的水平分力和豎直分力繞(xC,yC)點(diǎn)的力矩。MPdx,MPdy的表達(dá)式見附錄。即可求得筋材拉力合力：

(11)

式中，D為筋材拉力合力作用點(diǎn)的高度，當(dāng)筋材線性分布時：D=H/3,筋材均勻分布時：D=H/2。其中H為擋墻高度。

1.4 筋材拉力系數(shù)

采用無量綱化方式定義臨水加筋土擋墻的拉力合力系數(shù)為：

KT=∑Ti/(Wd+W*)。

(12)

2 內(nèi)部穩(wěn)定的遺傳算法分析

傳算法(Genetic Algorithm，GA)是一種通過模擬自然進(jìn)化過程搜索最優(yōu)解的方法，采用概率化的尋優(yōu)方法自動獲取和指導(dǎo)優(yōu)化的搜索空間，自適應(yīng)地調(diào)整搜索方向，通過對當(dāng)前群體施加選擇、交叉及變異等操作產(chǎn)生新一代群體，并逐步使群體進(jìn)化到最優(yōu)解狀態(tài)。

2.1 遺傳算法的計算過程

應(yīng)用遺傳算法解決實(shí)際問題時，首先要建立適應(yīng)度函數(shù)，然后根據(jù)遺傳算法的基本步驟求解問題的最優(yōu)解。

(1)目標(biāo)函數(shù)的確定。以(11)式為目標(biāo)函數(shù)，即Rt=∑Ti。Rt的值越大，表示其適應(yīng)性能越好。

(2)初始化種群。采用二進(jìn)制編碼，隨機(jī)產(chǎn)生一定數(shù)量的染色體，第j條染色體代表向量[β1j,β2j]，并且β1j，β2j滿足約束條件：

(13)

(3)適應(yīng)度函數(shù)。適應(yīng)度函數(shù)反映每一代個體對環(huán)境的適應(yīng)能力，即適應(yīng)度值越高的個體被遺傳到下一代的可能性越大。本研究求解擋墻筋材拉力的最大值，即設(shè)適應(yīng)度函數(shù)為：

F=Rt=f(β1j,β2j)，

(14)

(4)選擇操作。輪盤賭選擇是一種回放式隨機(jī)采樣方法，是遺傳算法中最常用的一種選擇方法，指某一個體被遺傳到下一代的概率等于它的適應(yīng)度值與整個種群中個體適應(yīng)度值之和的比值。[β1j,β2j]被選中作為父代的概率為：fj=Fj/ΣFj。

(5)交叉操作。按交叉概率pc從群體中隨機(jī)選擇兩個個體進(jìn)行某個或某些位的基因交換，從而產(chǎn)生新個體。交叉概率決定了遺傳算法的全局搜索能力，其值越大，收斂速度越快，但當(dāng)取值過大時，會破壞優(yōu)良個體，因此將交叉概率pc取0.7～0.9。常用的有單點(diǎn)交叉、多點(diǎn)交叉及均勻交叉等。本研究采用的單點(diǎn)交叉如圖2所示，隨機(jī)地在個體編碼串中設(shè)置一個交叉點(diǎn)，然后兩個配對的個體在該點(diǎn)相互交換部分染色體信息以產(chǎn)生新的個體。

(6)變異操作。變異操作是產(chǎn)生新個體的輔助方法，但也是必不可少的一個運(yùn)算步驟，它決定了遺傳算法的局部搜索能力。將交叉后的子代進(jìn)行變異操作，以變異概率pm對每個子代重復(fù)如下過程：在[0,1]內(nèi)取隨機(jī)數(shù)r，如果r

圖2 交叉操作Fig.2 Crossover operation

2.2 遺傳算法的計算結(jié)果

運(yùn)用遺傳算法的基本原理和算法程序?qū)Φ卣鹱饔孟屡R水加筋土擋墻筋材拉力的決策變量進(jìn)行編碼，將式(11)作為目標(biāo)函數(shù)，取案例計算中基本參數(shù)，設(shè)置初始種群數(shù)為30，迭代次數(shù)為30代，交叉概率pc=0.9，變異概率pm=0.01，對線性布筋的垂直擋墻筋材拉力進(jìn)行智能尋優(yōu)。從圖3可知，第1代初始種群的最大筋材拉力已達(dá)到88.45，與全局最優(yōu)解89.62已經(jīng)很接近，在第13代后進(jìn)化曲線趨于平滑。說明遺傳算法具有很好的收斂性。

圖3 尋優(yōu)性能的跟蹤Fig.3 Trace of finding best solution

3 案例計算

基本參數(shù)：擋墻墻高H=5 m，填土飽和重度γsat=19 kN/m3，填土干土重度γd=16 kN/m3，填土內(nèi)摩擦角φ=30°，擋墻垂直(α=90°)，kv=kh/2。力學(xué)模型如圖1所示。

文獻(xiàn)[9]假設(shè)填土具有較高的滲透性，假定直線型和折線型兩種滑裂面，只考慮臨水側(cè)靜水壓力，對浸水部分土體受到的地震慣性力進(jìn)行折減，計算墻背垂直的臨水加筋土擋墻筋材拉力系數(shù)。文獻(xiàn)[10]根據(jù)土的滲透性高低分為水土分算、水土合算兩種情況，考慮直線型和折線型兩種滑裂面計算墻背垂直的臨水加筋土擋墻筋材拉力系數(shù)。本研究取相同的參數(shù)分別與文獻(xiàn)[9]、文獻(xiàn)[10]的計算結(jié)果進(jìn)行對比，如表1所示。

表1 本研究方法計算得出的筋材拉力系數(shù)KT與其他方法的比較Tab.1 Comparison of tensile coefficient KT calculated by proposed method with those obtained by other methods

從表1可知，在kh=0.1，0.2及0.3時，在擋墻線性布筋時本研究方法分別比文獻(xiàn)[9]的筋材拉力計算結(jié)果大33.42%，33.62%及38.93%；均勻布筋下分別比文獻(xiàn)[9]計算結(jié)果大41.24%，37.00%及39.82%。本研究同時考慮臨水側(cè)靜水壓力和填土側(cè)靜水壓力，文獻(xiàn)[9]僅考慮臨水側(cè)靜水壓力，而忽略了填土側(cè)靜水壓力；由于假設(shè)墻后填土具有較高的滲透性，滑動土楔體的浸水部分所受到的水平地震慣性力應(yīng)僅作用于填土顆粒上[15]，而文獻(xiàn)[9]以填土干土重度與平均重度的比值對水平地震慣性系數(shù)進(jìn)行折減后與飽和重度相乘表示浸水部分受到的水平地震慣性力。因此文獻(xiàn)[9]的計算結(jié)果較小，本研究的計算結(jié)果更加合理。

在kh=0.05，0.1，0.15及0.2時，在擋墻線性布筋時本研究方法分別比文獻(xiàn)[10]的筋材拉力系數(shù)計算結(jié)果大2.0%，0.2%，-2.4%及-4.4%；均勻布筋下分別比文獻(xiàn)[10]計算結(jié)果大9.9%，5.8%，1.8%及-0.9%?？梢钥闯霰狙芯糠椒ㄓ嬎闼媒Y(jié)果與文獻(xiàn)[10]非常接近，驗證了本研究方法的可行性。

4 參數(shù)分析

取上文中的基本參數(shù)，考慮墻背垂直(α=90°)、墻背傾斜(α=70°)，筋材線性分布(D=H/3)、筋材均勻分布(D=H/2)，水位高度為hw/H=0.25，hw/H=0.75，填土內(nèi)摩擦角為φ=30°，φ=40°，超靜孔壓與填土豎向有效應(yīng)力的比值為ru=0.3，ru=0.7，水平地震加速度系數(shù)kh=0.05～0.2。

4.1 水位高度與墻高的比值

在實(shí)際情況中，水位高度是隨時發(fā)生變化的，但為了簡化模型，本研究假設(shè)擋墻兩側(cè)水位高度相同。取hw/H=0.25，0.75，水位高度對筋材拉力系數(shù)的影響如圖4所示。在相同墻背傾角、相同的布筋模式和相同地震力作用下，水位高度對筋材拉力系數(shù)有增大作用，且對墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)的增大作用較墻背垂直的擋墻更加明顯；對均勻布筋的臨水加筋土擋墻的筋材拉力系數(shù)的增大作用較線性分布更加明顯。如在kh=0.1時，從hw/H=0.25到hw/H=0.75，線性布筋和均勻布筋下相對應(yīng)的墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)KT分別增長66.88%，69.27%；垂直擋墻KT分別增長63.73%，64.58%。此外，hw/H=0.25時，從kh=0.05到kh=0.2，在線性布筋和均勻布筋下相對應(yīng)的墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)KT分別增長44.71%，34.13%，垂直擋墻的筋材拉力系數(shù)KT分別增長26.13%，19.01%；hw/H=0.75時，從kh=0.05到kh=0.2，在線性布筋和均勻布筋下相對應(yīng)的墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)KT分別增長54.16%，44.38%，垂直擋墻的筋材拉力系數(shù)KT分別增長33.48%，25.63%。在其他條件相同的情況下，水位較高時，地震慣性力對筋材拉力的影響更加顯著；地震慣性力對墻背傾斜(α=70°)的臨水加筋土擋墻的筋材拉力系數(shù)的影響較墻背垂直的擋墻更加顯著；地震慣性力對線性布筋的臨水加筋土擋墻的筋材拉力的影響較均勻分布更加顯著。

圖4 水位高度對KT的影響Fig.4 Influence of water level on KT

4.2 超靜孔壓與豎向有效應(yīng)力的比值

超靜孔壓與地震作用時的循環(huán)震動有關(guān)，是臨水加筋土擋墻地震穩(wěn)定性的重要影響因素。取ru=0.3，0.7，超靜孔壓對筋材拉力的影響如圖5所示。在相同墻背傾角、相同布筋模式和相同地震力作用下，超靜孔壓對筋材拉力有增大作用。如在kh=0.1時,從ru=0.3到ru=0.7，墻背傾斜的擋墻(α=70°)在線性布筋和均勻布筋下相對應(yīng)的筋材拉力系數(shù)KT分別增長24.92%，24.83%；垂直擋墻KT增長24.81%，24.63%。此外，ru=0.3時，從kh=0.05到kh=0.2，在線性布筋和均勻布筋下相對應(yīng)的墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)KT分別增長53.19%，44.38%，垂直擋墻的筋材拉力系數(shù)KT分別增長52.92%，43.49%；ru=0.7時，從kh=0.05到kh=0.2，在線性布筋和均勻布筋下相對應(yīng)的墻背傾斜的擋墻(α=70°)的筋材拉力系數(shù)KT分別增長33.49%，25.63%，垂直擋墻的筋材拉力系數(shù)KT分別增長32.77%，24.46%。在其他條件相同的情況下,超靜孔壓較小時，地震慣性力對筋材拉力的影響更加顯著；地震慣性力對線性布筋的臨水加筋土擋墻的筋材拉力的影響較均勻分布更加顯著。

圖5 超靜孔壓對KT的影響Fig.5 Influence of excess pore pressure on KT

4.3 填土內(nèi)摩擦角

圖6 填土內(nèi)摩擦角對的影響Fig.6 Influence of internal friction angle of filler on KT

墻土內(nèi)摩擦角反映了填土抗剪破壞的極限能力，φ值是由于土粒之間的滑動摩擦以及凹凸面間的鑲嵌作用產(chǎn)生的摩阻力的摩擦角，墻土內(nèi)摩擦角是工程設(shè)計中需要考慮到的重要參數(shù)。取φ=30°，40°，填土內(nèi)摩擦角對筋材拉力的影響如圖6所示。在相同墻背傾角、相同布筋模式和相同地震力作用下，填土內(nèi)摩擦角對筋材拉力有減小作用。如在kh=0.1時，從φ=30°到φ=40°，墻背傾斜的擋墻在線性布筋和均勻布筋下相對應(yīng)的筋材拉力系數(shù)KT約減小38.15%，39.82%；垂直擋墻KT約減小29.48%，30.73%。此外，在kh=0.05時，線性布筋模式下墻背傾斜的擋墻(α=70°)的筋材拉力在φ=30°，40°時分別比均勻布筋大15.74%，12.23%；線性布筋模式下垂直擋墻的筋材拉力在φ=30°，40°時分別比均勻布筋大7.18%，5.04%。在kh=0.2時，線性布筋模式下墻背傾斜的擋墻(α=70°)的筋材拉力在φ=30°，40°時分別比均勻布筋大8.34%，8.46%；線性布筋模式下垂直擋墻的筋材拉力在φ=30°，40°時分別比均勻布筋大0.91%，1.83%。在地震慣性力系數(shù)較大時，線性布筋模式下的擋墻筋材拉力系數(shù)與均勻布筋相近，且墻背垂直的擋墻較墻背傾斜的擋墻(α=70°)更加明顯。

5 結(jié)論

本研究假定滑裂面的型式為對數(shù)螺旋線，基于擬靜力方法建立極限平衡方程，通過合理地選擇遺傳算法控制參數(shù)和設(shè)置搜索范圍，有效地計算得出地震作用下臨水加筋土擋墻內(nèi)部穩(wěn)定設(shè)計所需的筋材拉力系數(shù)。根據(jù)結(jié)果對比和參數(shù)分析得出以下結(jié)論：

(1)通過計算對比分析驗證了本研究方法的可行性，表明了填土側(cè)靜水壓力對臨水加筋土擋墻筋材拉力系數(shù)計算的重要性。

(2)水位高度、超靜孔壓的增大對筋材拉力有增大作用，填土內(nèi)摩擦角的增大對筋材拉力有減小作用，且這些參數(shù)對墻背傾斜的擋墻(α=70°)的影響較墻背垂直的擋墻(α=90°)更加明顯；筋材拉力系數(shù)隨著地震慣性力的增大而增大，且在水位高度較大、超靜孔壓較小時其增大幅度更加顯著。

(3)在水位高度、墻背傾斜角度、填土性質(zhì)、超靜孔壓、地震荷載等條件相同時，臨水加筋土擋墻在線性布筋模式下的筋材拉力系數(shù)比均勻布筋模式小，即在相同的工況下，均勻布筋模式的臨水加筋土擋墻較線性布筋而言存在一定的設(shè)計冗余。該結(jié)論對工程應(yīng)用中的布筋方式有一定的參考價值。

(4)在水位高度、布筋模式、填土性質(zhì)、超靜孔壓、地震荷載等條件相同時，墻背傾斜的臨水加筋土擋墻的筋材拉力系數(shù)比墻背垂直的臨水加筋土擋墻小。在工程應(yīng)用中，在場地許可條件下，合理地將臨水加筋土擋墻墻背設(shè)計為傾斜狀可提升設(shè)計的經(jīng)濟(jì)性。