忠玉水電站雙曲拱壩高性能纖維混凝土抗爆性能分析

0 引言

目前，對水庫大壩等大型水利工程安全性的研究主要分為兩方面： ① 正常運行條件下水壩的穩定性，主要包括溫度[1]、滲流[2]、沉降[3]以及裂縫[4]等因素對水壩的影響； ② 抗爆防護，相關研究主要集中在水下爆炸對混凝土壩的結構影響和損傷破壞等方面[5-7]。邵蔚等[8]通過對現場振動數據進行分析，對水下爆破的濾頻效應特點進行了研究。鐘冬望對長江太子磯航道水下爆破工程的振動數據進行研究，發現主振頻率大部分為低頻波而少部分為高頻波。唐海等[10]以現場振動監測數據為基礎，提出了以巖體質點振速峰值作為控制標準。王洋等[1結合大連港大窯灣航道水下爆破工程，采用最小二乘法分析得到了Sadovsk公式中的衰減系數，為改良爆破設計提供了指導。劉釗等[12]為驗證水下爆破施工對沉管存放區干塢是否有影響，開展了一系列爆破試驗。

然而，鉆地武器對高壩大庫抗爆防護安全的威脅也不可忽視。一旦混凝土壩遭受侵徹爆炸打擊，其穩定性可能會大幅下降，發生潰壩將給下游帶來巨大的人員傷亡和財產損失[13-15]。Yang等[16]的研究表明，在混凝土受到侵徹爆炸作用時，初始侵徹損傷會影響混凝土的爆炸損傷;Shu等[1的研究結果顯示，在高速彈丸沖擊下，初始侵徹毀傷對混凝土壩的毀傷機制影響較小;李本平[18]對混凝土重力壩的毀傷效應進行了研究；張海榮[9采用LS-DYNA有限元分析程序研究了混凝土重力壩在侵徹后爆炸的流固耦合模型，探討了彈體整個侵徹和內爆過程對混凝土重力壩造成的破壞。此外，張甲文等2還對混凝土重力壩受到侵徹和爆炸聯合作用下的動力響應進行了仿真分析，并研究了打擊后壩頭部分的毀傷程度及安全性能；陳葉青等[21]指出在侵徹爆炸條件下，大壩易發生毀傷破壞而出現局部潰口。

上述研究表明，混凝土壩的侵徹爆炸作用對壩體的損傷具有不可忽視的影響，對混凝土壩的損傷具有主導作用。目前，對于混凝土壩在侵徹爆炸作用下的損傷機理已有較為豐富的成果，然而對于如何降低混凝土壩在侵徹爆炸作用下的損傷卻少有研究。本文從水壩材料出發，針對高性能纖維混凝土在爆破荷載作用下的損傷分布情況進行分析，以驗證其抗爆性能。

1 工程概況

忠玉水電站位于西藏那曲地區嘉黎縣境內（圖1），主要任務為發電，可研階段推薦的工程樞紐布置方案由混凝土雙曲拱壩、壩身孔口泄洪、壩后水墊塘消能、右岸地下廠房引水發電系統組成。工程規模為I等大（1）型工程。攔河壩采用對數螺旋線型雙曲拱壩，壩頂高程 3225m ，最大壩高 270m 。壩頂中心線弧長 402.08m ，底寬 54m ，厚高比0.2，最大倒懸度0.16，柔度系數8.41。

根據前期預可行報告研究，工程區在大地構造上位于岡瓦拉古陸的岡底斯-騰沖微陸塊，新構造上位于東構造結北側，地貌分區屬藏東南山地河谷區，現今構造活動強烈。庫區第四系土質岸坡存在局部蓄水后庫岸再造、塌岸問題。同時，所在流域具有年較差小、日較差大的氣溫變化特點。為確保忠玉水電站混凝土雙曲拱壩的安全施工及運營，研發了一種高性能纖維混凝土以滿足混凝土的質量和性能要求，并對該高性能纖維混凝土的抗爆性能進行驗證。

2 數值模型

2.1 模型建立

采用LS-DYNA建立高性能纖維混凝土及普通混凝土模型，在混凝土中間進行鉆孔及裝藥，藥包尺寸為 2.1cm×60.0cm 。單位制為 cm-g-μm ，數值計算采用流固耦合算法，混凝土視為固體，采用拉格朗日網格。空氣以及炸藥視為流體，采用歐拉網格。模型中的各邊界均為自由邊界，起爆方式為孔底起爆，模型示意見圖2（圖中 r 為半徑）。

2.2 材料模型及參數

（1）混凝土的RHT模型。高性能纖維混凝土由中熱硅酸鹽水泥、中粗河砂、石灰碎石或花崗巖碎石以及抗裂纖維拌和而成，其中纖維摻用體積率為 2% ，細骨料容積比為0.4，粗骨料容積比為0.5，水灰比為0.3，含氣量為 4%～7% 。與普通混凝土相比，其強度和抗裂性能得到改善，可適應西藏高海拔地區日夜溫差較大的特點。

混凝土材料參數選用RHT本構模型，其在充分考慮了彈性極限面、失效面以及損傷軟化段方程等基礎上將該模型分彈性、線性強化以及損傷軟化3個階段[22-25]，如圖3所示，可對混凝土在爆破作用下的損傷特性進行較好的描述。

在RHT模型中，需要輸入到LS-DYNA中的有38個參數。基本的物理和力學參數，如密度和靜態單軸抗壓強度，由試驗獲得。剪切模量、單軸抗拉強度可根據經驗確定。然后，部分參數可根據一系列理論計算得到（表1）。

（2）炸藥材料參數。炸藥材料采用LS-DYNA內置的MAT_HIGH_EXPLOSIVE_MURN，并加以狀態方程JWL_EOS對炸藥的爆破過程進行描述。

式中： P_cj 為爆炸壓力； V 為相對于未引爆狀態的體積；E₀ 為初始比內能： ;A，B，R₁，R₂ 和 ω 為爆炸常數。炸藥及狀態方程的相關參數如表2所示。

（3）空氣材料參數。空氣材料通常選用MAT_NULL材料及LINEAR_POLYNMIAL狀態方程加以描述，空氣壓力 P 為

P=C₀+C₁μ+C₂μ²+C₃μ³+（C₄+C₅+C₆μ²）E

式中： E 為單位體積內能； μ 為空氣壓縮應變， μ= （ρ/ρ₀）-1，ρ 和 ρ₀ 分別為當前空氣密度和初始空氣密度； C₀，C₁，C₂，C₃，C₄，C₅，C₆ 為空氣物質常數， C₄，C₅ 可用 C₄=C₅=γ-1 計算， γ 為比熱比值。

空氣材料參數及狀態方程的相關參數如下：空氣密度 ρ₀=1.29kg/m³ ， C₄=C₅=0.4，E=0.25MPa，V₀ =1.0 。

3 數值分析

3.1 損傷分析

在對混凝土的爆破數值分析中，混凝土的抗爆性能常以損傷變量來表征。為分析高性能混凝土的抗爆性能，提取模型中不同時刻的損傷云圖并與普通混凝土進行比較，對比結果見圖4。

炸藥引爆時，將產生大量高溫、高壓、高速氣流，激發沖擊波作用于炮孔孔壁，導致炮孔發生變形，并粉碎巖石、擴大爆腔，如圖4所示。隨后爆炸沖擊波衰減為爆炸應力波，在混凝土周圍形成拉伸變形和裂隙；然后爆生氣體填充爆腔，劇烈膨脹并侵入徑向裂隙，促使裂隙進一步擴展。最后彈性能被釋放成卸荷拉伸波，產生環向裂隙，將混凝土分割。

由圖4可知，在相同的爆炸荷載作用下，高性能纖維混凝土的損傷分布較普通混凝土少，兩者損傷分布主要差異體現在模型上部，即遠離起爆點一側。

由于采用的是孔底起爆，根據考慮傳播時間的應力波疊加理論可知，在遠離起爆點一側的應力疊加現象更明顯，因此在模型上部將具有更大的應力，在模型最上側的混凝土將更易造成損傷，形成剝離，進一步導致下面混凝土的損傷加劇。對于普通混凝土，由于其抗爆性能較弱，因此形成圖4（a）所示的損傷分布形態，而對于高性能纖維混凝土，抗爆性能較普通混凝土更優，其最上側混凝土能更好抵抗爆炸作用力造成的損傷，使其對內側混凝土仍能起到一定約束和保護作用，從而減少了整體損傷程度。

為進一步對兩者間的損傷差異進行定量分析，通過 LS-PrePost 后處理軟件獲取 t=1.00ms 時的損傷占比，結果如圖5所示。由圖5可知，普通混凝土的損傷占比為 72.3% ，高性能纖維混凝土的損傷占比為61.6% ，降低了 10.7% 。

由上述分析可知，高性能纖維混凝土較普通混凝土具有更好的抗爆性能。

3.2 有效應力分析

為進一步分析兩種混凝土損傷差異的原因，提取兩模型中的同一測點A（圖6的有效應力數據進行分析，結果如圖7所示。

由圖7可知，在該點處，普通混凝土的有效應力峰值為 75MPa ，高性能纖維混凝土的有效應力峰值為42MPa ，減少了 44% 。這是由于相較于普通混凝土，高性能纖維混凝土的強度有所增強，且拉壓應變率參數增大，使得在爆破作用下，高性能纖維混凝土的有效應力降低，裂紋擴展得到抑制。模型的同一點處，當爆炸荷載相同時，高性能纖維混凝土相較于普通混凝土具有更小的有效應力，對于模型中的一點，有效應力越小意味著會形成的損傷更小，混凝土的抗爆性能更佳。

注：外側圓環為高性能纖維混凝土，內側圓環為普通混凝土。

3.3 振速分析

爆破沖擊波衰減為地震波后，會引起建（構）筑物的振動，若不加以控制，爆破振動同樣會對建（構）筑物造成破壞。為此，以 10cm 為步長，各測點位置如圖8所示，提取兩模型中上部自由面各測點的振動速度峰值數據，分析其衰減規律，結果見圖9。

由圖9可知，在各測點處，高性能纖維混凝土的振動速度峰值均小于普通混凝土。這是由于兩者間的波阻抗不同，高性能纖維混凝土的振速衰減更快。相較于普通混凝土，高性能混凝土的振動峰值速度平均降低 45.93% 。

4結論

為評價和驗證高性能纖維混凝土的抗爆性能，采用LS-DYNA軟件建立了普通混凝土與高性能纖維混凝土模型并進行了對比分析，得到如下結果。

（1）在相同的爆炸荷載作用下，高性能纖維混凝土能更好地抑制損傷分布，尤其是對于模型上部的混凝土。

（2）在模型的同一測點處，高性能纖維混凝土的有效應力相較普通混凝土更小，同時高性能混凝土的振動峰值速度更低。

（3）高性能纖維混凝土具有較好的抗爆性能，能更好地應對突發極端情況，可為忠玉水電站混凝土雙曲拱壩的安全運行提供保障。

參考文獻：

[1］高山，陳建康，陳立寶，等.高溫炎熱地區碾壓混凝土重力壩溫控防裂研究[J].人民長江，2020，51（5）：154－159.

[2]李立平，李程程.河南省寶泉抽水蓄能電站下庫大壩滲流特征分析[J].水利水電快報，2023，44（1）：72-76.

[3]杜成斌，王作彤，黎亞生，等.考慮建基面的瀝青混凝土心墻壩竣工期沉降計算[J].人民長江，2023，54（7）：128-133.

[4]王元鵬，劉濤.楊房溝水電站拱壩9號橫縫二次張開原因及處理措施研究[J].水利水電快報，2023，44（增1）：26-29.

[5]WANGGH，LUWB，YANGGD，etal.A state-of- theartreview onblast resistanceand protection ofhigh dams toblast loads[J]. International Journal of Impact Engineering，2020，139：103529.

[6]周泰安.水下爆破振動對過江輸油管道影響研究[D].南寧：廣西大學，2022.

［7]黃迪.水下爆破振動對周圍建筑物影響研究［D].長沙：湖南大學，2019.

[8］邵蔚，王長柏.水下爆破振動特征及衰減規律研究[J].工程爆破，2018，24（5）：15-22.

[9]鐘冬望.太子磯航道水下鉆孔爆破地震波測試與分析[J].武漢科技大學學報，2011，34（5）：350-353.

[10]唐海，袁超，李俊如.武漢陽邏電廠擴建工程水下爆破安全閥值研究[J].采礦技術，2011，11（5）：186-188

[11]王洋，周明，姜殿科，等.某港區航道擴建水下炸礁爆破振動規律探索[J].中國水運（下半月），2019，19（10）：149 -150.

[12]劉釗，武浩文，孫文豪，等.水下爆破沖擊波對混凝土塢門影響試驗研究[J].中國港灣建設，2022，42（12）：33-36，136.

［13］李宏基.不同混凝土防護工程材料抗侵徹爆炸性能研究及數值模擬[D].南京：南京理工大學，2020.

［14］楊鑫.彈體侵徹與爆炸作用下大體積混凝土動態響應分析［D].哈爾濱：哈爾濱工程大學，2020.

[15]楊石剛，羅澤，許繼恒，等.侵徹爆炸作用下鋼纖維混凝土結構的破壞模式[J].爆炸與沖擊，2024，44（1）：151-163.

[16]YANG G D，WANG G H，LU W B，et al. A SPH -Lagrang-ian-Eulerian approach for the simulation of concrete gravi-ty dams under combined effects of penetration and explosion[J]. KSCE Journal of Civil Engineering，2018，22（8） ：3085-3101.

[17]SHUY Z，WANG GH，LU WB，et al.Damage characteris-tics and failure modes of concrete gravity dams subjected topenetration and explosion[J].EngineeringFailure Analy-sis，2022，134：106030.

[18］李本平.制導炸彈連續打擊下混凝土重力壩的破壞效應[J].爆炸與沖擊，2010，30（2）：220-224.

[19] 張海榮.混凝土重力壩在強沖擊加載下的力學響應[D].北京：北京理工大學，2016.

[20]張甲文，孟會林，盧江仁.混凝土重力壩在侵徹及爆炸加載下的仿真分析[J].彈箭與制導學報，2008，28（3）：126 -130.

[21]陳葉青，呂林梅，汪劍輝，等.爆炸沖擊荷載下的大壩抗爆性能及防護研究進展［J].防護工程，2021，43（2）：1-10.

[22]常銀會，楚京軍，侯榮，等.基于RHT本構模型的鋼渣混凝土SHPB模擬研究［J].混凝土世界，2023（12）：33-37.

[23］蔣霄.超高強和超高韌混凝土高速動態力學性能及其本構模型研究[D].杭州：浙江大學，2023.

[24]黃永輝，孫博，張智宇，等.巖石RHT本構的爆破碎裂判定方法優化及驗證[J].北京理工大學學報，2023，43（6） ：565 -574.

[25］李廣博，羅乙杰.基于RHT模型的沖擊荷載作用下鋼筋混凝土梁動力響應分析［J].吉林建筑大學學報，2022，39（1） ：11-16.

（編輯：李慧）

Analysis of blast resistance ofhigh-performance fiber concrete for double -curved arch dam of Zhongyu Hydropower Station

TIAN Jiangwei，FAN Tian，LI Yanhai，BAIMA Duoji （National Energy Group Xizang Electric Power Company Limited，Lhasa 85oooo，China）

Abstract：Reducing the damage of concrete dams caused bypenetrating explosions is of great significance for the safe constructionand operation of hydropower station dams.Inorder to ensure the safe construction and operationof theconcretehyperbolicarch damof Zhongyu Hydropower Station，the numerical calculation model of three-dimensional highperformance fiber concrete blasting was established by LS-DYNA basedon the three-dimensional non-linear finite elementanalysis method，and numerical simulationanalysis was cariedouttoanalyze the blast-resistantperformanceof high-performance fiberconcrete.Theresultsshowedthatcompared with ordinary concrete，high-performance fiber concrete can effectively inhibit the formation of damage with a reduction of damage ration at 11. 16% .At the same measurement point，the effective stress of high-performance fiber concrete was 42MPa ，which was a reduction of 44% compared with ordinary concrete of 75MPa .The vibration of the peak velocity of high -performance concrete reduced by 45.93% averagely，compared with ordinary concrete. On the whole，the high -performance fiber concrete has a stronger blast -resistant performance.

Key words：double-curved arch dam；high-performance fiber concrete；blast-resistant performance；numerical simulation；Zhongyu Hydropower Station