水利工程論文-巖石粘彈性模型辨識及參數(shù)反演摘要：本文回顧了常用的幾種巖石粘彈性本構模型，依據(jù)蠕變柔量的概念，用位移和蠕變柔量兩步反分析法從粘彈性本構模型的一般表達式中辨識出巖石的本構模型.先由相應的位移實測值用解析法反演巖石的蠕變柔量，再由巖石的蠕變柔量運用非線性優(yōu)化技術辨識出巖石本構模型表達式，最后進一步由模型參數(shù)反算出巖體的粘彈性參數(shù).文中最后給出了工程實例.關鍵詞：粘彈性本構模型蠕變柔量兩步反分析法反問題巖石巖體屬于各向異性流變介質，由于其客觀復雜性，在巖體理論分析和數(shù)值模擬方面，參數(shù)和模型的正確給定是巖石力學研究中的兩大難題，它們的合理性將大大增強巖體工程分析設計的可靠性.新奧法施工及目前廣泛應用的現(xiàn)場監(jiān)控法或信息反饋施工法是把巖石的變形觀測作為對原設計及施工方案進行修正的依據(jù)1，2.這些方法緊密結合工程實際指導修改原設計，收到了良好的效果，已廣泛應用于巖土工程實踐領域.參數(shù)取值，傳統(tǒng)的方法大多憑經(jīng)驗、工程類比來進行，常常不能作出科學的判斷.參數(shù)反演作為參數(shù)辨識的一種方法，是基于實測位移反求系統(tǒng)某種參數(shù)的一種逆問題，對確定認識系統(tǒng)和進一步正演分析相當重要.參數(shù)的辨識，包括模型參數(shù)的辨識，是在模型的結構式給定的情況下進行的，模型的給定與實際相符合的程度顯然相當重要.工程巖體是復雜的不確定系統(tǒng)，巖體的流變力學模型較多，模型識別的問題有著重要的理論意義和實用價值.在巖體本構模型的辨識，尤其是與時間有關的粘彈性本構關系辨識方面當前還有許多問題值得研究.本文回顧了工程應用和科學研究中幾種常用的巖石粘彈性本構模型，然后依據(jù)蠕變柔量概念，用位移和蠕變柔量兩步反分析法從粘彈性本構模型的一般表達式中辨識巖石本構模型，第一步由相應的位移實測值用解析法反演巖石的蠕變柔量，第二步由巖石的蠕變柔量運用非線性優(yōu)化技術辨識出巖石本構模型的表達式，最后進一步由模型參數(shù)與粘彈性參數(shù)的關系反算出巖石的粘彈性參數(shù).1粘彈性本構模型回顧工程巖體的長期穩(wěn)定性是當今巖土工程領域的一個十分重大的前沿問題3.工程巖體的變形破壞是有時間過程的，具有時間效應，巖體流變是巖體的重要基本力學特性和行為，對于某些實際巖體工程，研究變形的時間效應有著極其重要的意義8,9.真實巖體是一種非連續(xù)、非均質、各向異性的流變介質，目前很難用數(shù)學手段作出適當精確描述.在一定情況下，可將巖體視為似連續(xù)、均勻化、類各向同性介質，采用經(jīng)驗與理論相結合，定量分析定性使用的原則，對工程設計與施工起重要的定性指導作用10,11.在這一前提下，在巖土工程應用和科學研究中提出的粘彈性本構模型主要包括經(jīng)驗模型和組合模型5(彈性和粘性元件的不同組合而成).表1列出了常用的幾種微分型組合模型的公式及其特征(H代表彈性元件，N代表粘性元件).在一般狀態(tài)下，這些組合模型的本構方程微分形式的一般式可表示為5表1常用巖石粘彈性微分型組合本構模型式中：式(1)即線粘彈性微分型本構關系的一般表達式.2粘彈性模型識別2.1巖體介質蠕變柔量的確定在彈性巖體中開挖任意形狀斷面的隧道，應用平面復變影射、保角變換方法可得圍巖內任一點在曲線坐標中的位移為7由彈性理論確定.其中：()為垂直隧道軸線的物理平面Z上的非圓形洞室外域到數(shù)學平面內單位圓外域的映射函數(shù)；()和()為滿足應力邊界條件的復勢函數(shù).應用彈性粘彈性對應原理，對照式(2)，得出任意斷面的隧洞內任一點因開挖引起的粘彈性位移在曲線坐標中的表達式為8D(x)為荷載逐步釋放系數(shù)，若計算斷面遠離開挖面，則可作為平面應變問題處理；若接近開挖面，則嚴格說來是空間問題.為使問題簡化，考慮了開挖面空間效應后，仍作為平面問題來處理.開挖面空間效應通過洞周釋放荷載逐步釋放代替瞬間完全釋放來體現(xiàn).據(jù)研究資料4表明，其取如下形式：D(x)=1-0.7exp(-3.15x/2a).其中，a為隧洞半徑；p0為垂直向初始地應力；G為剪切模量；r為所考慮點到隧道軸線的距離；為側壓力系數(shù)；k為體積模量，D(x)為應力釋放系數(shù)，x為計算斷面距開挖面距離.在流變巖體中開挖任意形狀的隧道，在開挖之前布置量測儀器，可測到由于開挖引起的全部位移，在ti時刻可測到n個測點由于開挖引起的相對徑向全部位移，記為Uk(ti),k=1,2,，n,i=1,2,，l(假定量測了l個時段)，代入式(4)，則有蠕變柔量Jl(ti)和廣義蠕變柔量J2(ti)的線性方程組2.2粘彈性應力應變關系的確定線粘彈性微分型本構關系的一般表達式如式(1)所示，對其進行拉氏變換，考慮光滑化假定，得7對于(s)=0有重根或虛根的情況，J(t)中含有時間t的階次項和正余弦項，且未知參量個數(shù)增加，給反演優(yōu)化大大增加了難度，收斂性差，不利于程序的統(tǒng)一編制，這里作了一些簡化，在工程實際中，復雜的高階次模型也不必要.式(10)中，在模型階次確定時，參數(shù)Ai(i=12m+1)與模型式(1)中的模型參數(shù)p0,p1,p2,pm,q0,q1,q2,qn有一定的對應關系，由式(9)和式(10)可推導得出，模型參數(shù)均可用Ai簡單運算表達.表2常用巖石粘彈性模型參數(shù)表利用實測位移由式(6)可求出眾多不同時段ti的蠕變柔量Ji(ti)的最佳估計值，又由式(10)可得相對應的J(ti)(i=1L)，因而構造如下的非線性優(yōu)化目標函數(shù)：（11）式(10)是Ai的非線性函數(shù)，該問題是一非線性最小二乘優(yōu)化問題，帶約束的隱式非線性優(yōu)化問題，這里采用了可變容差優(yōu)化方法，它是在單純形法和復合形法的基礎上變化而來，具體見文獻6.上述模型參數(shù)估計是在模型已經(jīng)確定的情