張?zhí)?徐 磊

（河海大學水利水電學院，南京 210098）

合理確定邊坡巖體力學參數(shù)是正確評價其穩(wěn)定安全性的關(guān)鍵前提.傳統(tǒng)室內(nèi)、現(xiàn)場試驗與工程類比等方法確定邊坡巖體力學參數(shù)存在固有缺陷.近年來，根據(jù)位移監(jiān)測信息來反演邊坡工程巖體力學參數(shù)的智能反演分析方法得到快速發(fā)展，已經(jīng)成為解決邊坡工程巖體力學參數(shù)反演的有效方法.

目前，在巖土工程參數(shù)反演中應用較多的智能方法有人工神經(jīng)網(wǎng)絡、粒子群算法等［1－6］.人工神經(jīng)網(wǎng)絡反演方式需要進行樣本設計，其反演參數(shù)結(jié)果的精度直接取決于樣本的準確性.粒子群算法可進行直接反演，該算法易于實現(xiàn)全局收斂，效率較高，但是最優(yōu)粒子信息不能及時共享.鑒于此，為了使得邊坡巖體力學參數(shù)的反演更具高效性，采用一種比標準粒子群算法更高效、收斂更快的異步粒子群智能算法［7－14］，該算法的優(yōu)點在于提高信息共享的效率和加強信息交換的能力.

邊坡巖體力學參數(shù)實現(xiàn)方法的提出及其相應的數(shù)學描述還不足以使反演順利開展，除非采用功能強大的、高效的正分析工具.隨著計算機技術(shù)和數(shù)值計算方法的發(fā)展，有限元商業(yè)軟件在工程設計和科研領(lǐng)域得到了越來越廣泛的重視和應用，如ABAQUS具有強大的非線性分析功能、豐富的可模擬任意形狀的單元庫和與之相應的材料庫、高效求解器、更有力的網(wǎng)格劃分和強大的結(jié)果可視化技術(shù)，并已在包括邊坡工程在內(nèi)的巖土工程領(lǐng)域得到廣泛的應用.為此，將異步粒子群算法和ABAQUS兩者結(jié)合起來，開發(fā)出一套程序用于邊坡巖體力學參數(shù)的反演，并通過開展算例分析，驗證了所提出的方法和程序編制的可行性和高效性.

1 邊坡工程巖體力學參數(shù)反演模型

本文在已有研究成果的基礎(chǔ)上，建立高效反演模型用于反演邊坡工程中的巖體力學參數(shù).該反演模型采用ABAQUS作為正分析工具，建立以位移量為自變量的目標函數(shù)，構(gòu)建優(yōu)化反演數(shù)學模型，基于異步粒子群智能優(yōu)化算法進行反演.

在建立反演分析計算方法的研究中，變形觀測資料一般用作建立反演計算方程的輸入量，因而通常是進行反演計算的主要依據(jù).變形觀測資料主要包括位移量、應力量或應力增量的量測值，以及描述這些物理量隨時間而變化的規(guī)律的曲線等，采用位移反分析.

1.1 反演數(shù)學模型

在理想狀態(tài)下，監(jiān)測點的有限元位移計算值應與位移監(jiān)測值相一致，而目標函數(shù)值則用于判斷這種一致性.由于理論上和監(jiān)測上的偏差，這使得監(jiān)測點的位移計算值與監(jiān)測值存在誤差，因此優(yōu)化反演的目的致力于尋找使兩者之間的誤差為最小的解答，這類目標一般通過建立目標函數(shù)來實現(xiàn)，以測點位移實測值與有限元位移計算值之間的誤差作為目標函數(shù).

對于邊坡工程巖體力學參數(shù)反演，建立如下反演數(shù)學模型：

其中J：D→R為一個線性或非線性映射，數(shù)學表達式如下：

式中，Uci（s）為測點的有限元計算值；Uti（s）為監(jiān)測點實測值.為了求解最優(yōu)邊坡力學參數(shù)，要求在可行集合T

尋找最優(yōu)解sopt使得

在邊坡工程力學參數(shù)反演中，計算模型通常都很復雜，且計算規(guī)模很大，選擇一個高效的優(yōu)化算法十分必要.異步粒子群優(yōu)化算法（PSO）是一種基于標準粒子群算法的高效智能仿生進化算法，具有易于實現(xiàn)全局優(yōu)化、收斂速度快等優(yōu)點，為此，采用該算法來開展邊坡工程巖體力學參數(shù)反演.

1.2 異步粒子群算法

粒子群算法（PSO）是1995年由社會心理學家Kennedy和電氣工程師Eberhart提出的一種群體智能優(yōu)化算法，它通過個體間的協(xié)作與競爭來完成復雜搜索空間中最優(yōu)解的搜索.在協(xié)作過程中，每個粒子通過信息交流，共享最優(yōu)粒子的信息，進而更新自己來完成粒子群最優(yōu)解的搜索；在競爭過程中，粒子群通過每個粒子的適度值尋找到最優(yōu)粒子，粒子間通過競爭產(chǎn)生最優(yōu)粒子，不是最優(yōu)的粒子被群體丟棄，進而每個粒子向最優(yōu)粒子移動.在PSO算法中，優(yōu)化問題的潛在解都可以認為是d維搜索空間上的一個點，稱之為“粒子”（Particle），若干個粒子構(gòu)成一個粒子群.對于一個粒子群而言，它的每個粒子都有一個被目標函數(shù)決定的的適應度，每個粒子都具有位置和速度兩個特征，位置決定粒子的適應度，而速度決定下一代粒子運動的方向和距離.在每次迭代過程中，粒子通過跟蹤個體極值pBest和全局極值gBest來實現(xiàn)代際更新，直到在整個搜索空間中找到最優(yōu)解或達到最大迭代次數(shù)為止.

假設在一個n維的搜索空間中，由m個粒子構(gòu)成一個粒子群X＝［x1，…，xi，…，xm］，其中，第i個粒子的位置為xi＝［xi1，xi2，…，xin］T，其速度為vi＝［vi1，vi2，…，vin］T，其個體極值為pi＝［pi1，pi2，…，pin］T，全局極值為pg＝［pg1，pg2，…，pgn］T，則粒子xi將按下式改變其位置和速度以完成代際更新.

式中，j＝1，2，…，m，d＝1，2，…，n，k為粒子更新代數(shù)，wk為慣性權(quán)重，控制PSO算法的搜索能力，wk取值較大，則全局尋優(yōu)能力強，局部尋優(yōu)能力弱，wk通常依代數(shù)的增加而線性遞減.rand（）代表均勻分布于（0，1）之間的隨機數(shù)，c1、c2為加速常數(shù)，一般取2.0.式（5）中，第1部分為粒子先前行為的慣性，第2部分為“認知（cognition）”部分，表示粒子本身的思考；第3部分為“社會（social）”部分，表示粒子間的信息共享與相互合作.

對于粒子群算法，在上述粒子間協(xié)作過程中，每個粒子共享信息時的行為表現(xiàn)出同步性，即在每次迭代中，每個粒子通過自己代際間的適度值更新個體極值，粒子群通過同一個迭代步中的粒子適度值更新粒子群全局最優(yōu)值，最優(yōu)粒子的信息不能及時共享.為了提高信息共享的效率和加強信息交換的能力，從而減少計算機運行時間，提高優(yōu)化效率，在粒子間競爭過程中，當每個粒子的適度值小于共享信息的最小值時，粒子群的每個粒子把共享信息的最優(yōu)值作為個體最優(yōu)值，來進行粒子群更新，這使得粒子在協(xié)作競爭過程中表現(xiàn)出異步性，所以將它稱為異步粒子群算法（簡稱為APSO算法，下同）.異步粒子群算法（APSO）實現(xiàn)流程圖見圖1.

1.3 反演程序研制

提出了邊坡巖體力學參數(shù)的聯(lián)合反演法，該方法基于異步粒子群算法和通用有限元軟件ABAQUS，并考慮邊坡的系統(tǒng)錨桿等效模擬，編制了相應的計算程序.該程序在獲取有限元計算模型信息和反演參數(shù)信息的基礎(chǔ)上，首先采用隨機方法初始化粒子群，并基于每個粒子的位置，自動生成可以被ABAQUS求解器調(diào)用的INP文件來進行邊坡工程的有限元正算分析.進而計算粒子的適應度（即目標函數(shù)值）并基于異步粒子群算法完成粒子更新，直至找到符合預定最小閾值的最優(yōu)粒子或達到最大迭代次數(shù)為止，程序?qū)崿F(xiàn)如下.

圖1 異步粒子群算法流程圖

步驟1：根據(jù)反演參數(shù)的取值范圍和個數(shù)隨機初始化種群.

步驟2：調(diào)用ABAQUS進行有限元正演計算，讀入計算位移值，計算目標函數(shù)值.

步驟3：依據(jù)每個粒子的目標函數(shù)值，計算第一代粒子群的適應度、個體極值及全局極值.

步驟4：判斷適度值是否小于給定閥值.如果是程序結(jié)束并輸出結(jié)果，否則進入步驟5.

步驟5：按式（1）更新粒子的速度，按式（2）更新粒子的位置，調(diào)用ABAQUS進行有限元正演計算，計算目標函數(shù)值.然后進行異步粒子群更新，根據(jù)粒子適度值找到全局極值和全局最優(yōu)位置.

步驟6：判斷適度值是否小于給定閥值或者小于給定迭代次數(shù)，若不滿足則返回步驟5，否則程序結(jié)束并輸出結(jié)果.

2 算例分析

為了驗證本文所提出的APSO－ABAQUS聯(lián)合反演法在邊坡工程中巖體力學參數(shù)反演的可行性和高效性，進行如下算例分析.假定由兩種巖體構(gòu)成的邊坡工程分3步開挖，高為50m，上兩級坡比為1∶1，開挖高度為15m；第三級開挖坡比為1∶0.5，開挖高度為20m.上層巖體高為30m，下層巖體高為50m.本構(gòu)模型采用D－P彈塑性模型，巖體材料參數(shù)見表1.

表1 巖體力學參數(shù)

對于巖體中的初始地應力，取自重應力場作為初始地應力場進行計算.對于實測變形資料，根據(jù)工程經(jīng)驗，首先給定一組邊坡巖體力學參數(shù)（彈性模量為4 200MPa，摩擦系數(shù)為0.8，粘聚力為0.4MPa）進行有限元數(shù)值計算，并以測點的計算位移值作為反演所需的實測變形資料.邊坡工程分3次開挖，整體有限元計算模型如圖2所示，邊坡開挖體有限元計算模型如圖3所示.

圖2 整體有限元計算模型（圖中數(shù)字代表材料號）

圖3 邊坡開挖體有限元計算模型（圖中數(shù)字代表開挖步序）

隨后應用所研制的程序反演邊坡巖體力學參數(shù)，最后將反演所得的巖體邊坡材料參數(shù)與給定的材料參數(shù)進行對比分析，驗證了所提出方法及相關(guān)程序的可行性.

根據(jù)以上有限元模型，把材料號標為1的材料作為反演材料組，把材料號標為2的材料作為定值材料組.設置種群規(guī)模為20，學習因子c1＝c2＝2.0，慣性權(quán)重w1＝0.9，w2＝0.4，最大迭代次數(shù)取為100，迭代收斂精度為0.01.優(yōu)化算法程序迭代了19次，通過程序反演上層巖體的彈性模量、摩擦系數(shù)和粘聚力，其最終結(jié)果與初始值對比見表2，反演求解迭代過程曲線如圖4～5所示.

表2 反演參數(shù)對比結(jié)果

圖4 異步粒子群反演迭代過程曲線

圖5 粒子群反演迭代過程曲線

根據(jù)表2由邊坡實測開挖變形進行的彈性分析可以看出，反分析得出的巖體彈性模量與實測值相當接近，誤差可以控制在5%以內(nèi)，結(jié)果表明巖體對材料的彈性模量敏感，反分析得出的巖體摩擦系數(shù)和粘聚力與實測值有較大誤差，相對誤差可控制在15%左右，結(jié)果表明巖體對材料的摩擦系數(shù)和粘聚力不敏感，從而說明程序的實用性.

為了說明程序的高效性，給出了粒子群反演迭代過程曲線，對比可知異步粒子群算法的迭代次數(shù)明顯小于粒子群算法的迭代次數(shù).另外，基于前人在評價優(yōu)化算法上做的研究，對于兩種算法，選擇平均截止代數(shù)和截止代數(shù)分布熵來評估這兩種算法的收斂速度和收斂不穩(wěn)定性［15］.然后，為了評估優(yōu)化效率，構(gòu)建笛卡爾直角坐標系，以平均截止代數(shù)為X軸，以截止代數(shù)分布熵作為Y軸，分別確定異步粒子群算法和粒子群算法的坐標，即X值和Y值，以它們據(jù)原點的距離作為評價算法優(yōu)化效率的指標，距離越小，優(yōu)化效率越高.兩種算法的優(yōu)化效率如圖6所示.

對于多次獨立運行的PSO算法和APSO算法，從圖6可以明顯看出，APSO算法比PSO算法優(yōu)化效率要高，而且平均截止代數(shù)要小，收斂穩(wěn)定性較好.

圖6 兩種算法優(yōu)化效率對比

3 結(jié) 語

合理確定邊坡巖體力學參數(shù)是正確評價其安全性的關(guān)鍵前提.傳統(tǒng)試驗方法和工程類比方法和實際情況有一定誤差，在已有研究成果的基礎(chǔ)上，以功能強大的ABAQUS為正分析工具，建立以位移量為自變量的目標函數(shù)，構(gòu)建優(yōu)化反演數(shù)學模型；基于異步粒子群智能優(yōu)化算法，建立能夠反映主要影響因素作用的基于變形觀測資料的高效反演模型.研制的程序是在粒子群算法的基礎(chǔ)上，將粒子群算法的同步更新改進為異步更新的方式，從而每個粒子的更新都可以用到當代其他已經(jīng)更新的粒子的信息，提高信息共享的效率和加強信息交換的能力，可保證搜索到全局最優(yōu)解的同時，使算法更易于收斂.結(jié)合邊坡工程的算例，采用彈塑性計算模型，較好地模似邊坡工程結(jié)構(gòu)的開挖，成功地反演巖石的力學參數(shù)（彈性模量、摩擦系數(shù)、粘聚力）.程序運行結(jié)果表明，其計算精度是滿意的.聯(lián)合反演程序能夠?qū)崿F(xiàn)基于邊坡工程施工期變形觀測資料的巖體力學參數(shù)反演，驗證所提出的方法和程序編制的可行性和高效性.

［1］ 田明俊.智能反演算法及其應用［D］.大連：大連理工大學，2005.

［2］ Kavangh K，Clough R W.Finite Element Applications in the Characterization of Elastic Solids［J］.Int.J.Solids Structure，1971（7）：11－23.

［3］ Sakuali S，Takeuchi K.Back Analysis of Measured Displacement of Tunnel［J］.Rock Mech.and.Rock Eng，1983，16（3）：173－180.

［4］ 王芝銀，楊志法，王思敬.巖石力學位移反演分析問題及進展［J］.力學進展，1998，28（4）：488－498.

［5］ 蔣樹屏.巖體工程反分析研究的新進展［J］.地下空間，1995，15（1）：25－33.

［6］ 韓 峰，徐 磊，張?zhí)?壩基巖體力學參數(shù)的PSOABAQUS聯(lián)合反演［J］.河海大學學報：自然科學版，2013，41（4）：321－325.

［7］ 張 磊，高 尚.基于異步粒子群優(yōu)化算法的圖像分割方法［J］.微電子與計算機科學，2009，26（4）：174－177.

［8］ 劉蒙蒙.騾坪隧道施工監(jiān)測及其圍巖參數(shù)的智能反分析研究［D］.成都：西華大學，2009.

［9］ 沈永飛.邊坡位移反分析及其工程應用研究［D］.重慶：重慶大學，2010.

［10］張麗平.粒子群優(yōu)化算法的理論及實踐［D］.杭州：浙江大學，2005.

［11］賈善坡，伍國軍，陳衛(wèi)忠.基于粒子群算法與混合罰函數(shù)法的有限元優(yōu)化反演模型及應用［J］.巖土力學，2011，32（S2）：598－603.

［12］常曉林，喻勝春，馬 剛，等.基于粒子遷徙的粒子群優(yōu)化算法及其在巖土工程中的應用［J］.巖土力學，2011，32（4）：1077－1082.

［13］雷金波，陸曉敏，任青文，等.索風營水電站右岸地下廠房洞室群圍巖穩(wěn)定數(shù)值模擬［J］.長江科學院院報，2011，28（1）：47－52.

［14］王 偉，改進粒子群優(yōu)化算法在邊坡工程力學參數(shù)反演中的應用［D］.南京：河海大學，2007.

［15］孫瑞祥，屈梁生.遺傳算法優(yōu)化效率的定量評價［J］.自動化學報，2000，26（4）：552－556.