摘要：針對(duì)混凝土面板堆石壩施工質(zhì)量空間上的差異性導(dǎo)致的壩料物理力學(xué)參數(shù)具有空間差異性的問(wèn)題，提出了考慮實(shí)際施工質(zhì)量與填筑進(jìn)度（面貌）的精細(xì)有限元分析方法。首先基于試坑試驗(yàn)數(shù)據(jù)，采用空間插值算法對(duì)大壩壓實(shí)質(zhì)量進(jìn)行空間估計(jì)；然后建立各分區(qū)壩料本構(gòu)模型參數(shù)與壓實(shí)質(zhì)量之間的定量關(guān)系，實(shí)現(xiàn)對(duì)任意單元模型參數(shù)的空間估計(jì)與賦值；最后提出考慮填筑進(jìn)度（面貌）的精細(xì)建模方法，并結(jié)合工程實(shí)例對(duì)混凝土面板堆石壩的應(yīng)力變形等性態(tài)進(jìn)行了精細(xì)模擬與分析。結(jié)果表明：相比于基于設(shè)計(jì)參數(shù)模擬得到的沉降量平均誤差10.46 cm，采用考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的精細(xì)分析方法模擬得到的沉降量平均誤差僅為2.41 cm，其結(jié)果與實(shí)測(cè)值具有更好的一致性。研究成果為混凝土面板堆石壩的結(jié)構(gòu)性態(tài)精準(zhǔn)分析提供了新的途徑。

關(guān) 鍵 詞：混凝土面板堆石壩；施工質(zhì)量；填筑進(jìn)度；空間估計(jì)；精細(xì)有限元模型

中圖法分類號(hào)：TV641.43

文獻(xiàn)標(biāo)志碼：ADOI：10.16232/j.cnki.1001-4179.2024.10.019

0 引 言

混凝土面板堆石壩由于具備對(duì)復(fù)雜地形、地質(zhì)和氣候等適應(yīng)性強(qiáng)，施工工期短，成本低等特點(diǎn)，在水利水電工程中得到了大規(guī)模應(yīng)用^［1-3］，其安全問(wèn)題也受到了廣泛關(guān)注^［4］。混凝土面板堆石壩的實(shí)際施工質(zhì)量和填筑進(jìn)度是影響大壩結(jié)構(gòu)性態(tài)及服役安全的重要因素^［5-6］，大壩實(shí)際施工質(zhì)量的空間差異性導(dǎo)致了壩料物理力學(xué)特性在空間分布上存在差異性，其填筑進(jìn)度（面貌）對(duì)壩體的應(yīng)力變形等性態(tài)也會(huì)造成一定影響^［7］。因此，根據(jù)大壩實(shí)際施工質(zhì)量確定壩料物理力學(xué)參數(shù)的空間分布，并基于實(shí)際填筑面貌確定模擬施工步，是精確分析面板堆石壩結(jié)構(gòu)性態(tài)的重要前提。

采用不同本構(gòu)模型的有限單元法是現(xiàn)階段研究大壩結(jié)構(gòu)性態(tài)的重要手段。在壩料本構(gòu)模型的研究方面，張向韜等^［8］基于各向異性狀態(tài)參數(shù)，發(fā)展了考慮各向異性和顆粒破碎的堆石料彈塑性本構(gòu)模型，對(duì)已建的糯扎渡高心墻堆石壩進(jìn)行了應(yīng)力變形有限元計(jì)算分析；何亮等^［9］基于堆石料的靜動(dòng)力大三軸試驗(yàn)結(jié)果，建立了適用于堆石料的廣義塑性本構(gòu)模型，對(duì)高面板堆石壩展開(kāi)了數(shù)值模擬；黃焜等^［10］通過(guò)APDL參數(shù)化語(yǔ)言對(duì)鄧肯-張 E-B模型進(jìn)行了二次開(kāi)發(fā)，分析了猴子巖超高面板堆石壩在蓄水期的應(yīng)力和變形；Xu等^［11］建立了一種能夠在廣義塑性框架內(nèi)描述軸向應(yīng)力-應(yīng)變和體積應(yīng)變-軸向應(yīng)變關(guān)系的新型本構(gòu)模型，并對(duì)200 m高混凝土面板堆石壩進(jìn)行了數(shù)值研究。在本構(gòu)參數(shù)反演方面，陳家琦等^［12］基于改進(jìn)的自適應(yīng)混沌變異粒子群并行算法，對(duì)面板堆石壩不同材料分區(qū)的鄧肯E-B材料參數(shù)進(jìn)行了反演；莊文宇等^［13］提出了一種基于改進(jìn)自適應(yīng)遺傳算法和BP神經(jīng)網(wǎng)絡(luò)的反演方法，可針對(duì)多點(diǎn)監(jiān)測(cè)資料和非線性數(shù)值仿真進(jìn)行力學(xué)參數(shù)反演，有效避免了反演值的失真問(wèn)題；Ma等^［14］提出了一種基于多策略改進(jìn)Sooty Tern優(yōu)化算法的混凝土壩參數(shù)反演優(yōu)化方法，并在實(shí)際工程中得到了良好應(yīng)用。上述研究有效提高了大壩結(jié)構(gòu)性態(tài)的模擬精度，但在研究中均假定同一分區(qū)壩料的物理力學(xué)參數(shù)相同，未考慮實(shí)際施工質(zhì)量對(duì)壩料參數(shù)空間分布的影響以及實(shí)際填筑進(jìn)度對(duì)壩體應(yīng)力變形的影響，不能真實(shí)反映實(shí)際施工情況，存在一定的計(jì)算誤差。

為了解決傳統(tǒng)分析中壩料物理力學(xué)參數(shù)的均一性假設(shè)與實(shí)際壩料參數(shù)的空間差異性不符的問(wèn)題，并綜合考慮現(xiàn)場(chǎng)施工填筑面貌多變的實(shí)際情況，本文提出了一種考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的精細(xì)有限元分析方法，用以更精確地反映大壩的應(yīng)力變形情況，以為混凝土面板堆石壩的結(jié)構(gòu)性態(tài)精確分析提供一條有效途徑。

1 考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的精細(xì)模擬方法1.1 壓實(shí)質(zhì)量空間估計(jì)方法

1.1.1 工程樁號(hào)與實(shí)際坐標(biāo)系的匹配

對(duì)于壩軸線為曲線的大壩，工程樁號(hào)（包括樁號(hào)、壩縱）與實(shí)際坐標(biāo)系（x、y）不匹配。為實(shí)現(xiàn)基于工程樁號(hào)的壓實(shí)質(zhì)量數(shù)據(jù)與實(shí)際坐標(biāo)系的匹配，利用BP神經(jīng)網(wǎng)絡(luò)方法將工程樁號(hào)轉(zhuǎn)換為實(shí)際坐標(biāo)。

BP神經(jīng)網(wǎng)絡(luò)由接收樣本數(shù)據(jù)的輸入層、計(jì)算最終結(jié)果的輸出層以及隱藏層構(gòu)成^［15-16］。訓(xùn)練完成的網(wǎng)絡(luò)在接受輸入（樁號(hào)、壩縱）時(shí)，可以給出準(zhǔn)確的輸出（x、y坐標(biāo)）。對(duì)每個(gè)神經(jīng)元來(lái)說(shuō)，其輸入為

net=XW（1）

式中：X為輸入向量；W為連接權(quán)向量。隱藏層的激活函數(shù)為

f（net）=11+e^-net（2）

神經(jīng)網(wǎng)絡(luò)的訓(xùn)練以誤差（即實(shí)際輸出與理想輸出之差）控制在較低的范圍為目標(biāo)，調(diào)整網(wǎng)絡(luò)的權(quán)值。第p個(gè)樣本的誤差可描述為

E=12∑ni=1（y-O）（3）

式中：O為第p個(gè)樣本在第i個(gè)節(jié)點(diǎn)的實(shí)際輸出；y為第p個(gè)樣本在第i個(gè)節(jié)點(diǎn)的理想輸出；n為輸出節(jié)點(diǎn)數(shù)。

樣本集總誤差為

E=∑mj=1E（4）

式中：m為樣本總數(shù)量。當(dāng)求出滿足誤差要求的樣本時(shí)，模型訓(xùn)練結(jié)束。

1.1.2 基于KNN插值法的壓實(shí)質(zhì)量空間估計(jì)

基于試坑試驗(yàn)獲得的孔隙率、干密度等壓實(shí)質(zhì)量數(shù)據(jù)，利用KNN插值法估計(jì)壩體任意位置處的壓實(shí)質(zhì)量。KNN插值算法是一種基于實(shí)例學(xué)習(xí)的方法^［17-18］，其計(jì)算步驟主要包括以下3個(gè)方面。

（1）訓(xùn)練集X由已知的試坑點(diǎn)坐標(biāo)、干密度與孔隙率等數(shù)據(jù)構(gòu)成，X=（x，x，x;ρ，e）；測(cè)試集中的某點(diǎn)取為X=（x，x，x;ρ，e）。其中， q=1，2，…，l，l為試坑點(diǎn)總數(shù)；（x，x，x）為試坑點(diǎn)坐標(biāo)；ρ為干密度，g/cm³；e為孔隙率，%。

（2）訓(xùn)練集中各點(diǎn)X與測(cè)試集中某點(diǎn)X的距離按（5）式計(jì)算：

L（X，X）=∑lp=1（x-x）²（5）

式中：r=1，2，3，x代表試坑點(diǎn)的3個(gè)坐標(biāo)x，x，x。

（3）對(duì)求得的距離L進(jìn)行排序，選擇訓(xùn)練集中與X最近的k個(gè)鄰近點(diǎn)X（1 ≤ X≤ k），這k個(gè)近鄰點(diǎn)輸出的平均值即為X的輸出預(yù)測(cè)值^［17-18］。

ρ^=∑kt=1ρk（6）

e^=∑kt=1ek（7）

1.2 本構(gòu)模型參數(shù)空間估計(jì)方法

采用鄧肯-張E-B模型作為混凝土面板堆石壩壩料的本構(gòu)模型，并對(duì)其參數(shù)進(jìn)行空間估計(jì)。首先通過(guò)大型室內(nèi)三軸試驗(yàn)獲得不同壓實(shí)質(zhì)量下的鄧肯-張E-B模型參數(shù)值，然后基于三軸試驗(yàn)得到的模型參數(shù)規(guī)律，通過(guò)回歸分析建立壓實(shí)質(zhì)量與鄧肯-張E-B模型參數(shù)的回歸模型^［19］。

X^T=［φ，Δφ，K，n，R，K，m］^T［f（ρ），f（ρ），f（ρ），

f（ρ），f（ρ），f（ρ），f（ρ）］^TY^T（8）

式中：φ，Δφ分別為非線性強(qiáng)度參數(shù)；K為彈性模量數(shù)；n為彈性模量指數(shù)；R為破壞比；K為體積模量數(shù)；m為體積模量指數(shù)；f（ρ）為鄧肯-張E-B模型各參數(shù)與壩料壓實(shí)質(zhì)量的回歸函數(shù)；X^T為鄧肯-張E-B模型參數(shù)組合得到的轉(zhuǎn)置矩陣；Y^T為回歸函數(shù)f（ρ）組合得到的矩陣。

1.3 考慮填筑面貌的精細(xì)模擬方法

（1）三維模型建立與網(wǎng)格劃分。利用AutoCAD、Hypermesh等軟件建立混凝土面板堆石壩的幾何模型，并基于壩體分區(qū)與實(shí)際填筑面貌劃分網(wǎng)格。有限元網(wǎng)格以六面體結(jié)構(gòu)化網(wǎng)格為主，輔以少量四面體、五面體網(wǎng)格，以更精確地反映大壩的實(shí)際面貌，提高模擬精度。

（2）基于實(shí)際填筑面貌進(jìn)行荷載步劃分。將大壩三維填筑面貌模型導(dǎo)入ABAQUS軟件，根據(jù)施工進(jìn)度劃分模擬荷載步T（i=1，2，…，s），其中s為荷載步的總數(shù)量。為每個(gè)荷載步定義各自的單元集合TS，TS由含有壓實(shí)質(zhì)量信息的單元集合組成。

（3）以加載過(guò)程模擬實(shí)際填筑過(guò)程。采用ABAQUS軟件的單元生死功能（MODEL CHANGE），在初始階段將模型所有單元?jiǎng)h除，作為大壩尚未施工時(shí)的狀態(tài)；然后根據(jù)填筑順序依次激活相應(yīng)的單元集合LS，并為激活的單元施加重力荷載，直到壩體竣工，以模擬大壩的填筑施工；最后施加水荷載，以模擬大壩施工完成后蓄水的狀態(tài)。

（4）對(duì)模擬結(jié)果展開(kāi)分析。有限元模擬計(jì)算完成后，提取計(jì)算結(jié)果并與實(shí)測(cè)值進(jìn)行對(duì)比，驗(yàn)證計(jì)算結(jié)果的準(zhǔn)確性與有限元模擬的合理性。

2 實(shí)例分析

2.1 工程概況

五岳抽水蓄能電站的上水庫(kù)大壩為混凝土面板堆石壩，壩頂高程351 m，正常蓄水位347.50 m，最大壩高131 m，沿壩軸線長(zhǎng)1 110 m，主壩從左壩頭至右壩頭壩軸線依次由3段直線段和3段圓弧形線段交替連接組成。具體為：左壩頭起始端直線段，長(zhǎng)17.83 m，右接第一段圓弧線段，圓弧長(zhǎng)188.55 m（曲率半徑151 m，中心角71.546°）；右接第二段直線段，長(zhǎng)44.39 m，右接第二段圓弧線段，圓弧長(zhǎng)415.25 m（曲率半徑321 m，中心角74.118°）；再右接第三段直線段，長(zhǎng)217.10 m，右接第三段圓弧線段至右壩頭，圓弧長(zhǎng)226.88 m（曲率半徑181 m，中心角71.820°）。壩面采用鋼筋混凝土面板防滲，面板厚0.4 m，混凝土強(qiáng)度等級(jí)為C25，按12 m間隔設(shè)有垂直縫。壩體從上游至下游依次為墊層區(qū)、過(guò)渡區(qū)、主堆石區(qū)、下游堆石區(qū)。

2.2 大壩精細(xì)建模與填筑面貌方案

基于多個(gè)大壩斷面圖建立三維模型，并按照實(shí)際填筑進(jìn)度（面貌）來(lái)劃分網(wǎng)格。根據(jù)實(shí)際施工中分級(jí)填筑的層高，并參考類似工程有限元分析中的網(wǎng)格劃分情況^［20-22］，按10 m左右的網(wǎng)格高度模擬壩體填筑面貌，得到能夠反映平行上升面貌與實(shí)際填筑面貌的面板堆石壩三維有限元模型如圖1所示。

為了方便對(duì)比，選用兩種不同的填筑面貌方案進(jìn)行建模：第一種為平行上升的填筑形象面貌，壩體在全剖面同步水平上升，典型橫斷面的填筑面貌如圖2（a）所示，模擬加載步自壩基至壩頂分別為第1～14級(jí)。第二種方案為考慮實(shí)際填筑進(jìn)度的形象面貌，填筑順序?yàn)橄忍钪鲏螀^(qū)，再填筑左、右岸壩區(qū)，主壩區(qū)的典型橫斷面填筑形象面貌如圖2（b）所示，模擬加載步自壩基至壩頂分層交替上升，分別為第1～21級(jí)。部分典型施工階段的三維有限元模型如圖3所示。

2.3 大壩壓實(shí)質(zhì)量空間估計(jì)

依據(jù)DL/T 5128-2021《混凝土面板堆石壩施工規(guī)范》，基于大壩壓實(shí)質(zhì)量試坑試驗(yàn)數(shù)據(jù)，對(duì)大壩主堆石區(qū)、下游堆石區(qū)和過(guò)渡區(qū)進(jìn)行壓實(shí)質(zhì)量（干密度）的空間估計(jì)，結(jié)果如圖4所示。由空間估計(jì)結(jié)果可知，壩體內(nèi)部各斷面的干密度分布具有明顯的空間差異性。

2.4 大壩本構(gòu)模型參數(shù)空間估計(jì)

在鄧肯-張E-B參數(shù)中，彈性模量數(shù)K、彈性模量指數(shù)n、體積模量指數(shù)K、體積模量指數(shù)m對(duì)壓實(shí)質(zhì)量變化的敏感性較高，因此僅考慮上述4個(gè)參數(shù)隨干密度的變化，其余參數(shù)仍采用設(shè)計(jì)值。利用1.2節(jié)所述方法得到各分區(qū)的壩料干密度與K，n，K，m的回歸關(guān)系如表1所列。限于篇幅，僅列出主堆石區(qū)壩料的干密度與模型參數(shù)的回歸效果（圖5）。由圖5可見(jiàn)，各回歸模型的擬合決定系數(shù)R²均高于0.97，干密度與參數(shù)K，n，K，m的擬合關(guān)系良好，可認(rèn)為回歸模型能夠精確反映壩料干密度與本構(gòu)模型參數(shù)的關(guān)系。

通過(guò)上述回歸關(guān)系，并結(jié)合干密度的空間分布，進(jìn)一步計(jì)算得到各分區(qū)的本構(gòu)模型參數(shù)。限于篇幅，僅展示代表參數(shù)K和K的空間估計(jì)結(jié)果，如圖6所示。與主要分區(qū)的設(shè)計(jì)參數(shù)（表2）相比，空間估計(jì)參數(shù)在同一分區(qū)內(nèi)具有明顯的空間差異性，能夠更精細(xì)地表征壩體內(nèi)材料參數(shù)的空間分布。

3 計(jì)算結(jié)果與分析

3.1 施工面貌對(duì)大壩性態(tài)影響分析

基于設(shè)計(jì)參數(shù)分別對(duì)設(shè)計(jì)平行上升填筑與實(shí)際填筑進(jìn)度的兩種工況進(jìn)行模擬，并將模擬結(jié)果與實(shí)際監(jiān)測(cè)值比較，分析不同填筑施工面貌對(duì)大壩性態(tài)的影響。選擇2023年7月（對(duì)應(yīng)竣工期，即壩體填筑完成但未蓄水）的現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)作為實(shí)測(cè)值，豎直沉降量的測(cè)點(diǎn)布置于壩右0+003.35，如圖7所示。圖中，VSA21測(cè)線的測(cè)點(diǎn)從左至右分別為VSA21-1～VSA21-3，其他類似。

在兩種施工面貌工況下，各監(jiān)測(cè)點(diǎn)的沉降量模擬值與實(shí)測(cè)值對(duì)比如圖8所示。由圖8可知，在兩種工況下沉降量模擬值差異較小，多數(shù)測(cè)點(diǎn)的差值為1～3 cm之間，說(shuō)明施工面貌對(duì)大壩的沉降變形不會(huì)帶來(lái)根本性的規(guī)律改變；但在上游處靠近面板的范圍（即施工期的臨時(shí)邊坡處），兩種工況的沉降量模擬值差距相對(duì)較大，為3～5 cm，說(shuō)明實(shí)際填筑施工時(shí)形成的臨時(shí)邊坡會(huì)對(duì)堆石體的變形產(chǎn)生一定影響。

此外，兩種工況下模擬值與實(shí)測(cè)值差距均較大，采用平行上升填筑的模擬結(jié)果與采用實(shí)際填筑的模擬結(jié)果的沉降量誤差平均值分別為10.46 cm和9.81 cm，說(shuō)明采用設(shè)計(jì)參數(shù)模擬得到的沉降量與實(shí)測(cè)值有較大差異。因此，采用空間估計(jì)參數(shù)進(jìn)行精細(xì)有限元分析是必要的。

3.2 考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的精細(xì)分析

3.2.1 堆石體沉降量模擬值與實(shí)測(cè)值的對(duì)比

基于實(shí)際填筑進(jìn)度與空間估計(jì)參數(shù)模擬得到的竣工期沉降量如圖9所示。將沉降量的模擬值與實(shí)測(cè)值進(jìn)行對(duì)比，如表3和圖10所示。由圖10和表3可知，基于空間估計(jì)參數(shù)模擬得到的堆石體沉降量與實(shí)測(cè)值之間差異較小，各測(cè)點(diǎn)沉降量絕對(duì)誤差的平均值為2.41 cm；與基于設(shè)計(jì)參數(shù)的模擬誤差平均值（9.81 cm）相比，考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的模擬精度有了較大提升。

（蓄至正常蓄水位347.50 m）典型剖面的豎直沉降和水平位移分布。由圖可見(jiàn)，剖面的最大沉降量為36.85 cm，位于剖面中軸線300 m高程附近，相比于竣工期的最大沉降量模擬值（34.68 cm）增加了2.17 cm。剖面向上游方向的最大水平位移為0.17 cm，位于上游主堆石區(qū)高程280 m附近；向下游方向的最大水平位移為18.42 cm，位于下游次堆石區(qū)280 m高程附近。水平位移的分布規(guī)律基本為：上游堆石區(qū)位移指向上游，下游堆石區(qū)位移指向下游；受傾斜地形的影響，下游的水平位移總體上大于上游。

3.2.3 混凝土面板蓄水期變形預(yù)測(cè)

圖12為蓄水期面板撓度和沿壩軸向位移的分布。由圖12可見(jiàn)，面板由于水壓力產(chǎn)生的撓度和軸向位移分布符合面板壩的一般變形規(guī)律，其中面板中部約1/2到1/3壩高部位撓度較大，最大值為19.14 cm；向左、右岸方向的軸向位移最大值分別為3.5 cm和1.8 cm。

3.2.4 混凝土面板應(yīng)力預(yù)測(cè)

圖13為蓄水期面板順軸向應(yīng)力和坡向應(yīng)力的分布。由圖13可見(jiàn)，在水壓力的作用下，面板既存在壓應(yīng)力區(qū)又存在拉應(yīng)力區(qū)。軸向拉應(yīng)力最大值為1.5 MPa，壓應(yīng)力最大值為4.5 MPa；坡向拉應(yīng)力最大值為2.0 MPa，壓應(yīng)力最大值為3.0 MPa，二者的拉、壓應(yīng)力均在混凝土的強(qiáng)度范圍之內(nèi)。其中，坡向應(yīng)力表現(xiàn)為靠近趾板范圍受拉，面板中上部受壓，拉應(yīng)力較大的區(qū)域主要處于壩體高度最大的壩段中下部，符合面板壩的一般應(yīng)力規(guī)律（正值為拉應(yīng)力，負(fù)值為壓應(yīng)力）。

3.2.5 周邊縫變形預(yù)測(cè)

圖14為蓄水期周邊縫總變形的矢量分布。由圖14可見(jiàn)，蓄水期的周邊縫合變形（張拉、剪切和沉降的合位移）較小，最大值為5.6 cm，位于壩高最大的壩段范圍附近。周邊縫變形的模擬結(jié)果與同類工程的對(duì)比如表4所列?？梢?jiàn)，面板周邊縫的變形與同類工程相比較小。

4 結(jié) 論

為了研究面板堆石壩壓實(shí)質(zhì)量在空間上的差異所導(dǎo)致的壩料物理力學(xué)性態(tài)的空間差異性，并分析不同施工進(jìn)度（面貌）對(duì)于壩體應(yīng)力變形的影響，本文提出了考慮實(shí)際施工質(zhì)量與填筑進(jìn)度的面板堆石壩精細(xì)分析方法；基于工程設(shè)計(jì)填筑資料與壩料三軸試驗(yàn)，開(kāi)展了壓實(shí)質(zhì)量空間估計(jì)與本構(gòu)模型參數(shù)空間估計(jì)，進(jìn)而實(shí)現(xiàn)了考慮實(shí)際填筑進(jìn)度的有限元精細(xì)建模與模擬，并在實(shí)際工程中得到了應(yīng)用。具體結(jié)論如下：

（1）分別對(duì)平行上升的填筑面貌與考慮實(shí)際填筑進(jìn)度面貌的兩種工況進(jìn)行模擬，結(jié)果表明兩種工況的沉降量模擬值僅在實(shí)際填筑面貌的臨時(shí)邊坡處差異較大，為3～5 cm左右，說(shuō)明實(shí)際填筑時(shí)形成的臨時(shí)邊坡會(huì)對(duì)堆石體的沉降和水平變形產(chǎn)生一定影響；但其總體的變形規(guī)律不會(huì)發(fā)生根本性的改變。

（2）相較于采用同一分區(qū)完全一致的設(shè)計(jì)參數(shù)，由考慮實(shí)際施工質(zhì)量的空間估計(jì)參數(shù)模擬得到的變形結(jié)果與實(shí)測(cè)值具有更好的一致性，反映了考慮實(shí)際施工質(zhì)量的大壩性態(tài)模擬結(jié)果更為合理，精度較高。

（3）在蓄水期，預(yù)測(cè)得到的堆石體變形和面板變形的分布符合面板壩的一般變形規(guī)律；面板應(yīng)力值在混凝土的強(qiáng)度范圍內(nèi)，其分布也符合面板壩的一般應(yīng)力規(guī)律；大壩的周邊縫變形較小。模擬結(jié)果也驗(yàn)證了面板堆石壩在實(shí)際工程應(yīng)用中的優(yōu)越性。

參考文獻(xiàn)：

［1］ 王櫻畯，趙琳，雷顯陽(yáng)，等.抽水蓄能電站高瀝青混凝土面板堆石壩設(shè)計(jì)［J］.人民長(zhǎng)江，2022，53（1）：148-153.

［2］ 任翔，位敏，劉伍，等.極端條件下面板壩破壞規(guī)律及滲漏快速修復(fù)技術(shù)［J］.人民長(zhǎng)江，2019，50（增2）：119-122.

［3］ 殷鳴，胡清義，熊澤斌，等.高混凝土面板堆石壩接縫止水設(shè)計(jì)綜述［J］.人民長(zhǎng)江，2024，55（9）：165-172.

［4］ 常周梅，劉東海，陳輝，等.考慮壓實(shí)質(zhì)量和改進(jìn)本構(gòu)的瀝青心墻壩性態(tài)分析［J］.水力發(fā)電學(xué)報(bào)，2024，43（2）：86-98.

［5］ 劉東海，劉志磊，馮友文.堆石壩料壓實(shí)監(jiān)測(cè)指標(biāo)影響因素及適用性分析［J］.水力發(fā)電學(xué)報(bào)，2019，38（6）：1-10.

［6］ 劉東海，陳輝.數(shù)字化施工下高心墻堆石壩結(jié)構(gòu)性態(tài)精細(xì)數(shù)值模擬研究［J］.巖土工程學(xué)報(bào)，2019，41（3）：448-456.

［7］ 劉東海，張釗，陳輝.填筑形象對(duì)土石壩結(jié)構(gòu)性態(tài)影響的精細(xì)有限元分析［J］.河海大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，47（6）：514-522.

［8］ 張向韜，高溢釗，王翔南，等.基于廣義塑性模型的高心墻堆石壩反演計(jì)算分析［J］.水力發(fā)電學(xué)報(bào)，2020，39（9）：88-98.

［9］ 何亮，李雄威.基于廣義塑性模型的加筋面板堆石壩數(shù)值模擬［J］.水利水電技術(shù)，2020，51（1）：109-114.

［10］黃焜，蔡德所，涂小龍.猴子巖超高面板堆石壩應(yīng)力和位移非線性有限元分析［J］.水電能源科學(xué)，2018，36（5）：57-60.

［11］XU J，NING Z，WANG Y，et al.A novel constitutive model for the concrete face slab of a CFRD and the numerical simulation of its seismic behavior［J］.International Journal for Numerical and Analytical Methods in Geomechanics，2024，48（5）：1160-1196.

［12］陳家琦，岑威鈞，李鄧軍，等.基于改進(jìn)ACMPSO并行算法的土石壩本構(gòu)參數(shù)反演［J］.水利水電科技進(jìn)展，2021，41（3）：66-71.

［13］莊文宇，張如九，徐建軍，等.基于IAGA-BP算法的高拱壩-壩基力學(xué)參數(shù)反演分析［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，62（8）：1302-1313.

［14］MA L，MA F，CAO W，et al.A multi-strategy improved sooty tern optimization algorithm for concrete dam parameter inversion［J］.Water，2024，16（1）：119.

［15］SONG S，XIONG X，WU X，et al.Modeling the SOFC by BP neural network algorithm［J］.International Journal of Hydrogen Energy，2021，46（38）：20065-20077.

［16］袁旭山，劉京會(huì)，宋珂.基于BP神經(jīng)網(wǎng)絡(luò)的洪澇災(zāi)害承災(zāi)體脆弱性評(píng)估［J］.人民長(zhǎng)江，2024，55（2）：26-34.

［17］SAHOO A，GHOSE D K.Imputation of missing precipitation data using KNN，SOM，RF，and FNN［J］.Soft Computing，2022，26（12）：5919-5936.

［18］ZHANG S，ZHOU Y，YAN J，et al.Missing data completion method based on KNN and Random Forest［C］∥.Second LYSF Academic Symposium on Artificial Intelligence and Computer Engineering.SPIE，2021，12079：383-389.

［19］陳輝，劉東海，戚藍(lán).數(shù)字化施工下堆石壩模型參數(shù)空間估計(jì)及賦值［J］.巖土工程學(xué)報(bào)，2018，40（2）：278-286.

［20］江德軍，黃會(huì)寶，柯虎，等.土石壩填筑數(shù)值模擬中分層數(shù)影響研究［J］.人民長(zhǎng)江，2015，46（增1）：130-131，136.

［21］徐澤平，陸希，翟迎春，等.狹窄河谷中混凝土面板壩的應(yīng)力變形規(guī)律及工程措施研究［J］.水利學(xué)報(bào)，2022，53（12）：1397-1409.

［22］楊青坡，鄒德高，李云清，等.兩種殘余變形模型計(jì)算結(jié)果對(duì)比分析［J］.人民長(zhǎng)江，2015，46（4）：74-77.

（編輯：胡旭東）

Fine finite element analysis on panel dam considering actual construction

quality and filling progressTANG Desheng¹，XIE Pengfei¹，MA Ziru²，LIU Donghai²

（1.Sinohydro Bureau 5 Co.，Ltd.，Chengdu 610066，China； 2.State Key Laboratory of Hydraulic Engineering Intelligent Construction and Operation，Tianjin University，Tianjin 300350，China）

Abstract： Aiming at the problem that the spatial differences of construction quality of concrete face rockfill dams lead to the spatial difference of physical and mechanical parameters of dam materials，a fine finite element analysis method considering the actual construction quality and filling progress （dam appearance） was proposed.Firstly，based on the test data of the test pit，the spatial interpolation algorithm was used to estimate the compaction quality of dam in space.Then，the quantitative relationship between the parameters of the constitutive model of the dam material and the compaction quality was established，and the spatial estimation and assignment of the parameters of any element model were realized.Finally，a fine modeling method considering the filling progress （dam appearance） was proposed，and the stress and deformation of the concrete face rockfill dam were simulated and analyzed in combination with the engineering example.The results showed that compared with the average error of settlement （10.46 cm） based on the simulation of designed parameters，the average error of settlement obtained by the fine analysis method considering the actual construction quality and filling progress was only 2.41 cm，and the simulation results were in better agreement with the measured results.The research results provide a new way for the accurate analysis of structural behavior of concrete face rockfill dams.

Key words： concrete face rockfill dam；construction quality；filling progress；spatial estimation；fine finite element model