動(dòng)水荷載作用下粉土圍壩力學(xué)響應(yīng)研究

張宏博，于一凡，宋修廣

（1.山東大學(xué) 土建與水利學(xué)院，濟(jì)南 250061；2.山東省路基安全工程技術(shù)研究中心，濟(jì)南 250061）

水文水資源

動(dòng)水荷載作用下粉土圍壩力學(xué)響應(yīng)研究

張宏博1，2，于一凡1，2，宋修廣1，2

（1.山東大學(xué) 土建與水利學(xué)院，濟(jì)南 250061；2.山東省路基安全工程技術(shù)研究中心，濟(jì)南 250061）

為揭示動(dòng)水荷載作用下粉土圍壩力學(xué)參數(shù)的變化規(guī)律，通過(guò)室內(nèi)試驗(yàn)和數(shù)值模擬，分別研究了干濕循環(huán)及風(fēng)浪影響下粉土圍壩的力學(xué)響應(yīng)。結(jié)果表明：原狀粉土在經(jīng)歷2次以上干濕循環(huán)后，抗剪強(qiáng)度顯著降低，且隨干濕循環(huán)次數(shù)的增加，降低幅度逐漸減??；在周期性波浪荷載作用下，粉土圍壩迎水坡產(chǎn)生震蕩超孔隙水壓力，且隨時(shí)間動(dòng)態(tài)累積，有效應(yīng)力降低嚴(yán)重；孔隙水壓力增大到一定值后將趨于穩(wěn)定，孔壓穩(wěn)定時(shí)間隨波浪周期的增大而增大。該結(jié)論的得出為天然粉土壩坡的短期和長(zhǎng)期穩(wěn)定性評(píng)價(jià)、揭示粉土圍壩滑坡機(jī)理以及建立壩坡坍塌預(yù)測(cè)模型奠定了理論基礎(chǔ)。

粉土圍壩；干濕循環(huán)；波浪荷載；抗剪強(qiáng)度；孔隙水壓力；有效應(yīng)力

水庫(kù)在運(yùn)營(yíng)期間大壩處于開(kāi)放的自然環(huán)境中，長(zhǎng)期受到降雨入滲、波浪荷載沖擊、庫(kù)水位漲落等惡劣環(huán)境條件的影響，使筑壩材料物理力學(xué)性質(zhì)隨運(yùn)營(yíng)時(shí)間的增加而衰減。

黃泛平原區(qū)水庫(kù)大多為就地取土的圍壩型水庫(kù)，筑壩材料為土性不良的粉性土，由于該類(lèi)土具有弱可塑性、低黏結(jié)性、高分散性等性質(zhì)，且具有極強(qiáng)的水敏感性特征，在庫(kù)水位漲落、降雨滲透、波浪等動(dòng)水荷載作用下，壩體強(qiáng)度極易降低，嚴(yán)重影響粉土圍壩的安全運(yùn)營(yíng)，因此有必要對(duì)粉土的干濕循環(huán)特性及波浪荷載作用下的動(dòng)力響應(yīng)進(jìn)行研究。目前，國(guó)內(nèi)外主要研究成果多是針對(duì)膨脹土［1-2］、黏性土的干濕循環(huán)強(qiáng)度衰減規(guī)律及粉土海床［3-5］的波浪荷載動(dòng)力響應(yīng)問(wèn)題，而對(duì)粉土圍壩的水—土耦合作用規(guī)律研究較少，本文將對(duì)動(dòng)水荷載（干濕循環(huán)、波浪荷載）作用下粉土圍壩的力學(xué)響應(yīng)進(jìn)行研究，以期為粉土圍壩壩體坍塌預(yù)測(cè)防護(hù)和治理提供理論依據(jù)。

1 粉土干濕循環(huán)試驗(yàn)

在雨水入滲蒸發(fā)或庫(kù)水位漲落等因素影響下，圍壩長(zhǎng)期處于“濕潤(rùn)—風(fēng)干”的循環(huán)狀態(tài)，干濕循環(huán)效應(yīng)對(duì)土體的強(qiáng)度有較大影響。為研究粉土強(qiáng)度參數(shù)隨干濕循環(huán)次數(shù)的變化規(guī)律，設(shè)計(jì)了無(wú)荷條件下的干濕循環(huán)試驗(yàn)（實(shí)際情況下，除表層土體外，均承受豎向荷載）。

1.1 干濕循環(huán)試驗(yàn)方案

試驗(yàn)所用粉土取自坡面以下1.0m深處，取好的土樣用保鮮膜密封好后放入保濕箱內(nèi)保存，整個(gè)取土、保存過(guò)程中盡量避免對(duì)土樣的擾動(dòng)。

首先采用烘干法測(cè)定了原狀粉土的天然含水率為20.56%，然后制取原狀粉土12個(gè)環(huán)刀直剪試件（h=20mm，φ=61.8mm），將試件分為6組，每組2個(gè)試件，其中第1組做原狀土的直剪試驗(yàn)，其余5組分別進(jìn)行1次，2次，3次，4次，5次干濕循環(huán)，然后按照試驗(yàn)規(guī)程做經(jīng)不同干濕循環(huán)次數(shù)的直剪試驗(yàn)，測(cè)定試件的抗剪強(qiáng)度，然后根據(jù)剪切定律確定粉土的抗剪強(qiáng)度指標(biāo)φ和c，每組取兩個(gè)試件的平均值。試件放在陽(yáng)光下進(jìn)行干燥（在夏季進(jìn)行），待試樣質(zhì)量接近土顆粒質(zhì)量 （根據(jù)含水率確定）時(shí)用滴管滴水使試件逐漸濕潤(rùn)，盡量保證滴入的水分能均勻的滲透到試件中及濕潤(rùn)后的試件與原狀土的含水率相同或相近，因每個(gè)試件的原質(zhì)量基本相同，濕潤(rùn)后只需保證試件的總質(zhì)量與原狀土試件的總質(zhì)量相等即可。

1.2 試驗(yàn)結(jié)果

按照上述試驗(yàn)方案對(duì)原狀粉土經(jīng)不同干濕循環(huán)次數(shù)后的抗剪強(qiáng)度指標(biāo)進(jìn)行室內(nèi)試驗(yàn)，其結(jié)果如表1。

表1 抗剪強(qiáng)度指標(biāo)

從表1可以看出原狀粉土的抗剪強(qiáng)度隨著干濕循環(huán)次數(shù)的增加而逐漸衰減：①黏聚力在第1次干濕循環(huán)之后相對(duì)0次循環(huán)下降了13.96%，內(nèi)摩擦角下降了9.89%；②黏聚力在第2次干濕循環(huán)之后相對(duì)1次循環(huán)下降了12.64%，內(nèi)摩擦角下降了8.86%；③黏聚力在第5次干濕循環(huán)之后相對(duì)4次循環(huán)下降了1.36%，內(nèi)摩擦角下降了1.52%。

可以看出黏聚力和內(nèi)摩擦角在1次、2次干濕循環(huán)后下降的幅度較大，之后隨著干濕循環(huán)次數(shù)的增加變化幅度趨于平緩。經(jīng)歷5次干濕循環(huán)后，黏聚力衰減了32.06%，內(nèi)摩擦角衰減了26.24%。總之，原狀粉土在經(jīng)歷不少于2次干濕循環(huán)后，土顆粒排列順序發(fā)生變化，改變了土體內(nèi)部的應(yīng)力狀態(tài)，致使其內(nèi)部結(jié)構(gòu)遭到疲勞損傷破壞，抗剪強(qiáng)度發(fā)生較大衰減，對(duì)于粉土結(jié)構(gòu)物，其安全系數(shù)將顯著降低。

2 波浪荷載數(shù)值模擬

平原水庫(kù)蓄水后，形成較寬廣的水面，風(fēng)區(qū)擴(kuò)大也使吹程急劇增加，因此庫(kù)區(qū)風(fēng)場(chǎng)所造成的風(fēng)浪波高也得到了較大提高，風(fēng)浪對(duì)邊坡穩(wěn)定的影響已不可忽略。下面將用FLAC軟件計(jì)算波浪作用下粉土圍壩的動(dòng)力響應(yīng)，包括孔隙水壓力響應(yīng)、有效應(yīng)力響應(yīng)、總應(yīng)力響應(yīng)、土體液化及循環(huán)波浪作用下粉土圍壩孔壓穩(wěn)定時(shí)間。

2.1 計(jì)算工況

根據(jù)平原水庫(kù)現(xiàn)場(chǎng)實(shí)際情況，按照莆田試驗(yàn)站波浪要素計(jì)算公式［6］，并考慮波浪破碎時(shí)自由水面與壩坡接觸位置處的波壓力，選取3種計(jì)算工況如表2。

表2 計(jì)算工況

波浪在自由水面處產(chǎn)生的波壓力將沿水深發(fā)生衰減，沿深度方向直線變化，波浪產(chǎn)生的波壓力隨時(shí)間和水平位置的變化規(guī)律為：

式中 p0為波壓力幅值；x為相對(duì)于節(jié)點(diǎn)的水平坐標(biāo)；L為波的長(zhǎng)度；T為波周期；t為時(shí)間。

2.2 模型建立

2.2.1 模型尺寸

壩高15m，壩頂寬6m，壩坡1∶2.5，壩前水位10.5m，壩后水位2m，壩底基礎(chǔ)寬106m，高20m，整個(gè)模型厚度方向取1m，3個(gè)方向的網(wǎng)格點(diǎn)距離均為1m，在壩前坡水面以下0.5，5.5，10.5m處布置了6個(gè)監(jiān)測(cè)點(diǎn)，測(cè)點(diǎn)A、B、C位于泥面下1m，測(cè)點(diǎn)D、E、F位于泥面下2m，如圖1。實(shí)時(shí)輸出6個(gè)位置處的動(dòng)力響應(yīng)。

圖1 模型示意圖

2.2.2 邊界條件

模型基礎(chǔ)底面z方向固定，左右兩側(cè)x方向固定，沿縱向y方向位移為零；模型基礎(chǔ)底面不透水。

2.2.3 粉土參數(shù)

重度20kN/m3，體積模量1.62e6，剪切模量5.4e5N/m2，孔隙率0.4，內(nèi)摩擦角30°。

2.3 數(shù)值模擬

2.3.1 粉土圍壩初始平衡狀態(tài)

施加波浪荷載前，在壩前壩后分別施加靜水壓力，壩前坡腳處?kù)o水壓力105kPa，壩后坡腳處?kù)o水壓力20kPa，沿深度方向的變化率10kPa，使粉土圍壩在靜水條件下達(dá)到平衡狀態(tài)。

粉土圍壩的初始平衡如圖2。靜水條件下，同一縱斷面內(nèi)孔隙水壓力隨深度線性分布，飽和土在每個(gè)測(cè)點(diǎn)的孔隙水壓力等于該位置處的靜水壓力。

從圖2中可以看出，孔隙水壓力由壩前向壩后遞減，總應(yīng)力與有效應(yīng)力沿深度遞增，基本沿壩中心線對(duì)稱(chēng)分布，剪應(yīng)變?cè)隽孔畲笾蛋l(fā)生在壩頂兩側(cè)，說(shuō)明此處受到的剪應(yīng)力最大。

圖2 粉土圍壩的初始平衡狀態(tài)云圖

2.3.2 粉土圍壩波浪響應(yīng)

模型在靜水壓力作用下達(dá)到初始平衡狀態(tài)，之后開(kāi)始計(jì)算動(dòng)力問(wèn)題。打開(kāi)動(dòng)態(tài)模塊和滲流模塊，土體改用Finn模型，然后施加波浪荷載進(jìn)行動(dòng)力計(jì)算。波浪荷載屬于周期動(dòng)荷載，隨時(shí)間變化。FLAC3D不能直接施加波浪荷載，本文使用FISH語(yǔ)言對(duì)模型施加波浪荷載。

下面以工況2為例分析粉土圍壩的波浪響應(yīng)。

圖3 施加波浪荷載后粉土圍壩應(yīng)力云圖

施加波浪荷載后，壩體內(nèi)部應(yīng)力云圖發(fā)生明顯改變，壩前孔隙水壓力明顯增大，壩前有效應(yīng)力明顯減小，甚至出現(xiàn)正值，而總應(yīng)力基本不變，剪應(yīng)變?cè)隽吭龃蠓秶蕡A弧狀，由壩頂延伸至壩底。

圖4 監(jiān)測(cè)點(diǎn)孔壓與總應(yīng)力時(shí)程變化（負(fù)值為總應(yīng)力，正值為有效應(yīng)力）

循環(huán)波浪荷載對(duì)圍壩邊坡的影響主要是通過(guò)改變土體內(nèi)部的孔隙水壓力來(lái)實(shí)現(xiàn)的，從圖3、圖4的孔隙水壓力、有效應(yīng)力、總應(yīng)力變化圖中可以看出，在波浪荷載對(duì)粉土圍壩的循環(huán)波壓力作用下，壩體坡面一定深度范圍內(nèi)的有效應(yīng)力場(chǎng)和孔壓場(chǎng)都隨著時(shí)間和空間發(fā)生變化，壩坡不同深度位置處的孔隙水壓力產(chǎn)生了周期性的震蕩，并隨時(shí)間產(chǎn)生動(dòng)態(tài)累積，以致出現(xiàn)震蕩超孔隙水壓力，而波浪作用過(guò)程中總應(yīng)力無(wú)累積現(xiàn)象，根據(jù)有效應(yīng)力原理，土體中有效應(yīng)力隨時(shí)間逐漸降低，甚至導(dǎo)致邊坡土體液化，從圖2中可以看出，壩坡一定深度范圍內(nèi)的有效應(yīng)力已大于0，土體已部分液化，喪失抵抗剪切的能力，甚至出現(xiàn)壩坡滑塌或潰壩的潛在危險(xiǎn)。當(dāng)孔隙水壓力增大到一定值后將趨于穩(wěn)定，在一定范圍內(nèi)波動(dòng)。

從圖4中可以看出，深度較深處總孔隙水壓力值越大，且孔隙水壓力增量越大，深度較淺處孔隙水壓力震蕩幅值較大。同一深度位置處，靠外的測(cè)點(diǎn)的孔隙水壓力大于靠?jī)?nèi)的測(cè)點(diǎn)。

2.3.3 孔壓穩(wěn)定時(shí)間

波浪荷載作用下，壩坡不同位置處的孔隙水壓力增大到一定值后將趨于穩(wěn)定，為此計(jì)算上述3種工況下粉土圍壩的動(dòng)力響應(yīng)，得出6個(gè)監(jiān)測(cè)點(diǎn)的孔壓穩(wěn)定時(shí)間隨波浪周期的變化規(guī)律如圖5。

圖5 各測(cè)點(diǎn)孔壓穩(wěn)定時(shí)間隨波浪周期的變化

從圖5可以看出，對(duì)于壩坡不同深度位置處，其孔壓穩(wěn)定時(shí)間均隨波浪周期的增大而增大，越靠近水面處，由于波壓力幅值較大，土體的孔隙水壓力很快達(dá)到其峰值，有效應(yīng)力迅速減小，故其孔壓穩(wěn)定時(shí)間越短。由于水平方向測(cè)點(diǎn)間距較小，其孔壓穩(wěn)定時(shí)間相差不大。不同周期波浪作用下，對(duì)于A、D測(cè)點(diǎn)，其孔壓穩(wěn)定時(shí)間為3T～4T（T為周期）；對(duì)于B、C、E、F測(cè)點(diǎn)，其孔壓穩(wěn)定時(shí)間為5T～7T。

3 結(jié)語(yǔ)

（1）原狀粉土在經(jīng)歷不少于2次干濕循環(huán)后，其內(nèi)部結(jié)構(gòu)遭到疲勞損傷破壞，抗剪強(qiáng)度發(fā)生較大衰減，之后隨著干濕循環(huán)次數(shù)的增加衰減幅度趨于平緩。

（2）波浪荷載作用下，粉土圍壩壩坡不同深度位置的孔隙水壓力產(chǎn)生周期性的震蕩，并隨時(shí)間產(chǎn)生動(dòng)態(tài)累積，使土體中有效應(yīng)力逐漸降低，導(dǎo)致邊坡土體液化，土體抗剪強(qiáng)度大幅降低。

（3）波浪荷載作用下，壩坡孔隙水壓力增大到一定值后將趨于穩(wěn)定，此孔壓穩(wěn)定時(shí)間隨波浪周期的增大而增大。

（4）該結(jié)論的得出為天然粉土壩坡的短期和長(zhǎng)期穩(wěn)定性評(píng)價(jià)、揭示粉土圍壩滑坡機(jī)理及建立壩坡坍塌預(yù)測(cè)模型奠定了理論基礎(chǔ)。

（5）實(shí)際情況下，粉土圍壩的干濕循環(huán)效應(yīng)是在有荷條件下發(fā)生的，因此，有荷條件下粉土的干濕循環(huán)試驗(yàn)有待進(jìn)一步補(bǔ)充，以確定真實(shí)情況下粉土的抗剪強(qiáng)度指標(biāo)隨干濕循環(huán)次數(shù)的變化規(guī)律；對(duì)不同坡比、水深等工況下粉土圍壩的波浪響應(yīng)需做進(jìn)一步的數(shù)值模擬和模型試驗(yàn)，以期為不同條件下粉土圍壩壩體坍塌預(yù)測(cè)防護(hù)和治理提供依據(jù)。

［1］楊和平，張銳，鄭健龍.有荷條件下膨脹土的干濕循環(huán)脹縮變形及強(qiáng)度變化規(guī)律［J］.巖土工程學(xué)報(bào)，2006，28（11）:1936-1941

［2］呂海波，曾召田，等.膨脹土強(qiáng)度干濕循環(huán)試驗(yàn)研究［J］.巖土力學(xué)，2009，30（12）:3797-3802.

［3］張民生.波浪作用下黃河三角洲海床穩(wěn)定性研究［D］.大連：中國(guó)海洋大學(xué)，2006.

［4］武伯弢.波浪荷載作用下砂質(zhì)粉土海床的液化研究［D］.上海：同濟(jì)大學(xué)，2009.

［5］吳夢(mèng)喜，樓志剛.波浪作用下海床的穩(wěn)定性與液化分析［J］.工程力學(xué)，2002，19（5）:97-102.

［6］顏宏亮，閆濱.水工建筑物［M］.北京：中國(guó)水利水電出版社，2012.

Mechanical response research of silt dam under hydrodynamic load

ZHANG Hong-bo1,2，YU Yi-fan1,2，SONG Xiu-guang1,2
（1.School of Civil Engineering，Shandong University，Jinan 250061，China；2.Subgrade Security Engineering Research Center of Shandong Province，Jinan 250061，China）

In order to reveal the change rule of mechanical parameters of silt dam under hydrodynamic load,mechanical response of silt dam under cyclic wetting-drying and stormy waves was studied respectively through laboratory test and numerical simulation.The results showed that shear strength of undisturbed silt decreased significantly after more than 2 times of cyclic wetting-drying,and the decreasing amplitude decreased gradually with the increase of cycle number.Excess pore water pressure of silt dam upstream slope oscillated under wave load and accumulated dynamically with time,which led to decrease of effective stress seriously.The pore water pressure would tend to be stable when it increased to a certain value.The pore pressure stabilization time increased as the wave period increased.This conclusion established the theoretical foundation for stability evaluation of natural silt dam slope in short term and long-term,revealing landslide mechanism of silt dam and prediction model establishment of dam slope collapse.

silt dam;cyclic wetting-drying;wave load;shear strength;pore water pressure;effective stress

P341

A

1672-9900（2015）02-0001-05

2015-03-06

山東省科技發(fā)展計(jì)劃（2011GGH21607）；濟(jì)南市科技發(fā)展計(jì)劃（201218001）

張宏博（1978-），男（漢族），山東濟(jì)寧人，副教授，主要從事巖土工程與道路工程的研究工作，（Tel）15552582796。