楊俊鴿

（洛陽水利勘測設(shè)計有限責任公司，河南 洛陽 471000）

1 引言

堤防作為水工設(shè)施中重要控流、 防洪建筑，其安全穩(wěn)定性與滲流場變化乃是結(jié)構(gòu)設(shè)計中必須考慮的因素[1，2]。 擋墻作為一種防護結(jié)構(gòu)，其可用作堤防防洪、防滲及加固，故開展擋墻結(jié)構(gòu)設(shè)計優(yōu)化對推動堤防運營安全水平具有重要意義。舒天白等[3]、茹秋瑾等[4]采用ANSYS、ABAQUS 等數(shù)值仿真軟件開展了擋墻的靜力場分析，從應力、位移特征評價結(jié)構(gòu)設(shè)計的技術(shù)優(yōu)勢。 華中等[5]、劉波[6]認為水工設(shè)計不應僅研究靜力穩(wěn)定性，利用FLUENT、MIKE21等滲流場模擬平臺開展了流場特征計算， 分析了不同設(shè)計方案下水工建筑的滲流場演變， 為工程設(shè)計評價提供了依據(jù)。 不論是靜力場計算或是滲流場分析，上述方法集中在數(shù)值計算過程，陳林[7]采用模型試驗手段在室內(nèi)原型復制建立水工模型，通過施加相似比例尺寸的荷載工況，分析模型在運營過程中滲流參數(shù)、應力與位移參數(shù)等變化，也可為工程設(shè)計提供更為高精度的參照。 本文在水庫堤防運營現(xiàn)狀的前提下， 采用FLAC 3D 與FLUENT 仿真平臺開展了加筋格柵擋墻設(shè)計優(yōu)化分析，為堤防選擇最適配擋墻結(jié)構(gòu)提供參考。

2 工程設(shè)計模擬

2.1 工程概況

衛(wèi)家磨水庫乃是靈寶市重要水利中樞， 原設(shè)計蓄水庫庫容為320 萬m3，正常蓄水位19.5m。 該水庫堤防工程最大高度為21.6m，無防浪墻及二次加高擋墻等結(jié)構(gòu)設(shè)施，原堤防基礎(chǔ)乃是原位地基，僅進行防滲墊層及注漿加固處理， 堤防斷面形狀如圖1 所示。 由于衛(wèi)家磨水庫工程堤防兩側(cè)岸坡穩(wěn)定性欠佳，且堤防地基土體以軟弱黏土為主，沉降變形較大，特別是岸坡內(nèi)易形成潛在滑移面，對堤防工程帶來較大威脅。 因而，考慮采用加筋擋土邊墻結(jié)構(gòu)形式對該岸坡進行加固， 其設(shè)計形式如圖2 所示，而加筋方式采用土工格柵結(jié)構(gòu)[8]，不同于其他擋墻結(jié)構(gòu)設(shè)計， 格柵設(shè)計勢必需要契合邊墻，與之受力、變形及防滲保持一致性，才可最大程度發(fā)揮邊墻對衛(wèi)家磨水庫堤防的加固防護作用。 因而，本文重點對帶格柵結(jié)構(gòu)的擋墻結(jié)構(gòu)設(shè)計優(yōu)化分析。

圖1 堤防斷面圖

圖2 加筋擋土墻設(shè)計幾何圖

2.2 工程設(shè)計

采用FLAC 3D 平臺建立加筋帶格柵擋土邊墻結(jié)構(gòu)設(shè)計[9]，如圖3 所示，本模型覆蓋防護范圍為堤防上、下游各15m，而本文研究模型的流場及靜力穩(wěn)定性均考慮在該區(qū)域內(nèi)。 經(jīng)有限元網(wǎng)格劃分后共獲得節(jié)點單元121 652 個， 節(jié)點數(shù)142 685個， 該模型中墻后岸坡巖土層物理力學參數(shù)直接以土工試驗實測取值，包括土體滲透系數(shù)等。 該擋墻的扶壁、面板厚度參數(shù)分別為1.2m、1m，土工格柵分層鋪設(shè)， 經(jīng)本文計算表明格柵層間距參數(shù)應處于合理區(qū)間， 否則不利于邊墻發(fā)揮防滲與加固作用，而格柵層間距又不應超過扶肋厚度，故本文對比研究方案中格柵層間距設(shè)定為0m （無加筋）、0.2m、0.4m、0.6m、0.8m、1m、1.2m，擋墻結(jié)構(gòu)其他設(shè)計參數(shù)保持一致， 探討加筋擋墻格柵設(shè)計參數(shù)對結(jié)構(gòu)靜力場及流場影響特性。

圖3 擋墻有限元模型

不論擋墻是否加筋設(shè)計，其結(jié)構(gòu)安全穩(wěn)定性反映在擋墻上則是各重要部位的應力、 位移特征，圖4 中即是擋墻結(jié)構(gòu)的面板、踵板及扶壁等部位。 在該模型中分別改變格柵設(shè)計參數(shù)，典型方案的模型中格柵布設(shè)如圖5 所示。 模型中設(shè)定X、Y、Z 正向分別為擋墻剖面下游向、墻后土壓力方向及結(jié)構(gòu)自重應力向。 滲流場模擬FLUENT 開展計算[10]，通過圖2 幾何模型的導入， 采用三維滲流場分析單元模塊完成流域內(nèi)滲流場特征參數(shù)計算。

圖4 擋墻結(jié)構(gòu)重要組成部分

圖5 格柵布設(shè)圖

3 格柵設(shè)計對結(jié)構(gòu)靜力場影響

3.1 應力特征

依據(jù)格柵層間距各設(shè)計方案， 獲得擋墻上重要部位拉應力受層間距設(shè)計參數(shù)影響變化特征，如圖6 所示。 分析圖中拉應力變化可知，不論格柵層間距設(shè)計參數(shù)水平如何， 擋墻結(jié)構(gòu)上拉應力最大乃是面板迎水側(cè)， 各方案中該部位拉應力分布為2.25—5.25MPa，而頂板、扶肋等部位拉應力較前者分別減少了16.1%—32.2%、29.6%—50.1%。從抗拉安全設(shè)計考慮， 面板迎水側(cè)迎水側(cè)最易受到張拉應力集中效應， 所選用的結(jié)構(gòu)材料應滿足較大抗拉要求。 當格柵層間距增大后， 面板迎水側(cè)、扶肋及頂板上的最大拉應力均隨之遞減，但不可忽視層間距過大，拉應力降幅減小。 分析表明，格柵布設(shè)間距愈大， 對結(jié)構(gòu)張拉應力的分散調(diào)整具有促進作用， 更好提升擋墻抗拉效果； 另一方面，格柵層間距不應過大，不利于加筋格柵在擋墻內(nèi)分布，應控制層間距設(shè)計參數(shù)在合理區(qū)間即可。當格柵層間距在方案0—0.8m 內(nèi)時， 在各方案間，頂板、 扶肋及面板迎水側(cè)拉應力隨之平均降幅為22.7%、25%、18.7%，而層間距超過0.8m 后，前三個部位拉應力的平均降幅分別為0.7%、1.2%、0.9%。依據(jù)此三部位拉應力變化可初步得知， 格柵層間距最優(yōu)方案不需要超過0.8m。 而擋墻上另一踵底板部位拉應力在各方案中均處于較穩(wěn)定狀態(tài)，維持在1MPa，各方案中拉應力最大波幅不超過1%。從拉應力的計算結(jié)果表現(xiàn)分析， 認為當格柵層間距0.8m 時結(jié)構(gòu)抗拉技術(shù)優(yōu)勢最顯著。

圖6 擋墻各部位最大拉應力特征

根據(jù)結(jié)構(gòu)靜力場計算， 同樣獲得了重點部位壓應力受層間距參數(shù)影響特征，如圖7 所示。 分析壓應力變化可知，踵底板部位處受壓最大，各方案中其最大壓應力分布為11.2—22.3MPa，即結(jié)構(gòu)自重應力乃是決定受壓效果的最主要因素。 當格柵層間距參數(shù)遞增后， 各部位最大壓應力均隨之先增后減變化， 而所有方案中壓應力最大均位于間距0.8m 方案中，在該方案下踵底板、扶肋及面板迎水側(cè)壓應力分別為22.3MPa、19.1MPa、16.3MPa。當格柵層間距過低時，如位于0—0.8m 方案中，當間距每增大0.2m，各部位壓應力平均可增長19.1%、19.5%、20%，該區(qū)間內(nèi)方案對結(jié)構(gòu)張拉應力的重分布調(diào)整效果達到較佳；但當層間距超過0.8m 后，與拉應力受之逆轉(zhuǎn)性變化類似， 各部位壓應力隨之遞減，且降幅較大，其中扶肋壓應力在0.8—1.2m區(qū)間內(nèi)最大降幅超過21.9%， 平均降幅為18.5%，而踵底板與面板迎水側(cè)在該區(qū)間內(nèi)的平均降幅分別也達14.9%、19.9%，故結(jié)構(gòu)設(shè)計時不應調(diào)整層間距超過0.8m，此與拉應力受之影響相契合。結(jié)合擋墻結(jié)構(gòu)拉、壓應力影響變化特征，筆者認為格柵層間距0.8m 時靜力場穩(wěn)定性最優(yōu)。

圖7 擋墻各部位最大壓應力特征

3.2 位移特征

依據(jù)各設(shè)計方案獲得三向位移隨格柵層間距設(shè)計參數(shù)變化特征，如圖8 所示。 分析位移特征觀察可知，Z 向位移乃是各方案中最大，分布為7.9—14.8mm， 此也印證了擋墻結(jié)構(gòu)受力來源主要為結(jié)構(gòu)自重。 隨格柵層間距增大，各向位移均隨之為先減后增變化，表明格柵層間距參數(shù)的改變，會導致結(jié)構(gòu)位移變化出現(xiàn)階段性改變。 當層間距設(shè)計參數(shù)位于0—0.8m 方案區(qū)間內(nèi)時， 在方案0.8m 處Y向位移較之方案0.2m、0.6m 下分別減少了46.2%、14.9%，當層間距設(shè)計參數(shù)每增大0.2m，則Y 向位移平均降低20.2%； 當層間距參數(shù)位于0.8—1.2m時，Y 向位移隨之具有平均增幅40.9%， 方案間最大增幅可達44.2%， 不利于擋墻結(jié)構(gòu)運營穩(wěn)定性。X、Z 向位移變化特征與Y 向位移有所類似， 層間距設(shè)計參數(shù)在0—0.8m 時，X、Z 向位移在各方案間具有平均降幅24.4%、17.2%，而在0.8—1.2m 方案內(nèi)分別具有平均增幅60%、31.8%， 最大增幅分別可達85.6%、33.6%。從三向位移受格柵層間距影響來看，以X 向位移敏感度最大[11]。綜合分析認為，格柵層間距0.8m 方案在各向位移中均處于較佳狀態(tài)，有利于擋墻結(jié)構(gòu)穩(wěn)定性。

圖8 擋墻各向位移特征

4 格柵設(shè)計對堤防滲流場影響

本文以FLUENT 仿真平臺模擬計算獲得堤防滲流場特征， 圖9 為各層間距方案下斷面流速變化特征。

圖9 斷面流速變化特征

從圖中可看出，格柵層間距愈大，則堤防斷面流速水平愈大，在層間距0.2m 方案內(nèi)，全斷面上平均流速為0.34m/s，而層間距每增大0.2m 后，斷面平均流速遞增24.6%。流速水平的高低決定了堤防區(qū)間斷面內(nèi)滲流場穩(wěn)定性， 過大的流速有利于流場內(nèi)沖淤排沙， 但會導致水工結(jié)構(gòu)受水力勢能沖刷影響，同樣過低的流速會限制泥沙沉降，不利于控流[12]， 故合理穩(wěn)定的流速才是堤防運營最需要的。 從斷面流速的穩(wěn)定變化來看，以層間距0.6m、0.8m 更為穩(wěn)定， 兩者方案內(nèi)斷面流速分別穩(wěn)定在0.56m/s、0.62m/s，而層間距在0.2m、0.4m 內(nèi)流速在上、下游分別具有較大波幅，斷面間流速最大波幅分別超過14.7%、15.5%， 而層間距1m、1.2m 方案內(nèi)全斷面流速參數(shù)均出現(xiàn)波動狀態(tài)，穩(wěn)定性欠佳，極易出現(xiàn)局部紊流、渦旋等非穩(wěn)定滲流現(xiàn)象。 綜合滲流場變化與靜力場安全穩(wěn)定性， 認為格柵層間距0.8m 時擋墻運營綜合技術(shù)優(yōu)勢最大。

5 結(jié)論

（1）擋墻上拉應力最大位于面板迎水側(cè)，面板迎水側(cè)、 扶肋及頂板上的最大拉應力均隨格柵層間距參數(shù)遞減，但降幅減小，在層間距0.8m 方案后降幅處于停滯狀態(tài)； 踵底板拉應力穩(wěn)定在1MPa；擋墻最大壓應力隨設(shè)計參數(shù)為先增后減變化，在0.8m 方案中最大。

（2）擋墻位移最大在Z 向；各向位移值均為先減后增變化，以兼具0.8m 方案下位移值最低。

（3） 層間距設(shè)計參數(shù)與斷面流速具有正相關(guān)關(guān)系；流速不應過大或過小，同時也應保持較為穩(wěn)定，層間距在0.6m、0.8m 方案內(nèi)流速穩(wěn)定性最佳。

（4） 綜合擋墻結(jié)構(gòu)設(shè)計靜力穩(wěn)定性與滲流場變化，認為格柵層間距0.8m 為技術(shù)優(yōu)勢最顯著。□