泥石流沖擊荷載下攔擋壩的動力響應分析

第一作者周勇男，博士，教授，1978年10月生

泥石流沖擊荷載下攔擋壩的動力響應分析

周勇1，2，劉貞良1，2，王秀麗1，2，周鳳璽1，2

(1.蘭州理工大學甘肅省土木工程防災減災重點實驗室，蘭州730050；2.西部土木工程防災減災教育部工程研究中心，蘭州730050)

摘要：選取典型的泥石流沖擊荷載，把鋼筋混凝土攔擋壩簡化成懸臂梁，建立其動力偏微分方程，運用數學物理方程中的變量分離法和結構動力學中的振型疊加法求得動力偏微分方程的解析解，最后根據材料力學得到相應的應力以及應變。算例分析表明解析解是可靠的，且在攔擋壩的壩頂處有最大的位移，為實際工程力學計算提供一種思路。

關鍵詞：泥石流攔擋壩；沖擊荷載；動力響應分析

基金項目：國家科技支撐計劃項目(2011BAK12B07)；甘肅省自然科學基金資助項目(1208RJZA158)；蘭州理工大學紅柳青年教師培養計劃資助項目(Q201108)；甘肅省高等學校基本科研業務費項目

收稿日期：2013-11-01修改稿收到日期：2014-05-14

中圖分類號:TH212；TH213.3文獻標志碼：A

Dynamic response analysis for a dam against impact load of debris flow

ZHOUYong1,2,LIUZhen-liang1,2,WANGXiu-li1,2,ZHOUFeng-xi1,2(1. Key Laboratory of Disaster Prevention and Mitigation in Civil Engineering of Gansu Province, Lanzhou University of Technology, Lanzhou 730050, China; 2. Western Engineering Research Center of Disaster Mitigation in Civil Engineering ofMinistry of Education, Lanzhou University of Technology, Lanzhou, 730050, China)

Abstract:Here, the typical impact load of debris flow was chosen, and a reinforced concrete dam against debris flow was simplified into a cantilevered beam, its dynamic partial differential equation was established. Then, the variable separation method in mathematical physics equations and the modal superposition method in the structural dynamics were used to get the analytical solution to the dynamic partial differential equation. According to the mechanics of materials, the corresponding stress and strain of the beam were obtained. The example analysis showed that the analytical solution to the dynamic equation is reliable, and the top of the dam against debris flow has the largest displacement. The study results provided an idea for practical engineering mechanical calculation.

Key words:dam against debris flow; impact load; dynamic response analysis

目前，泥石流重力式攔擋壩是根據《泥石流災害防治工程設計規范》[1]設計的，即采用類似擋土墻方法設計。吳瑋江等[2]在規范的基礎上規范和統一了重力式攔擋壩的穩定性計算方法和公式。喬振華[3]介紹了泥石流地區攔擋壩的作用及適用條件，從多個方面，闡述了攔擋壩的平面布置原則，分析探討了攔擋壩的高度、間距及其構造。王曉明[4]采用綜合對比、逐項分析的方法，對實體重力壩、普通格柵壩和鋼構格柵壩的工程防治效果進行了對比分析。熊道錕等[5]分析了攔擋壩壩體破壞、壩基破壞和壩肩坡壞的主要原因，提出了防止攔擋壩潰決的措施。李冠奇等[6]提出了按有效匯流面積計算泥石流峰值流量的優化設計方法，用該方法設計的攔擋壩的溢流斷面面積較小，從而達到提高設防高度，增大庫容，提高攔擋壩的容災能力。賈世濤[7]在野外調查的基礎上，結合防護體系配置優化的新思路展開試驗研究。在室內二維水槽模型試驗中，通過改變泥石流容重、一次過程體積總量和模型壩的開孔率等控制參數，得出泥石流性能在過壩前、后的變化規律。王根龍等[8]提出用鋼混結構重力式攔擋墻替代原有漿砌塊石重力式攔擋墻，滿足了壩體對沖擊力與強度的關系。何曉英等[9]進行了泥石流漿體沖擊特性模型實驗，采用小波方法揭示泥石流漿體沖擊特性。綜上所述，人們都是將泥石流的沖擊荷載作為一個靜力荷載作用在鋼筋混凝土攔擋壩上，從而進行靜力學分析。而泥石流沖擊荷載是動力荷載，本文采用結構動力學的方法，建立了泥石流沖擊荷載與鋼筋混凝土攔擋壩的動力方程，分析了攔擋壩的動力響應。

1泥石流沖擊荷載的選擇

泥石流沖擊荷載模型很多。何思明等[10]以Hertz接觸理論為基礎，考慮結構的彈塑性特性，給出泥石流大塊石沖擊力的計算方法。何思明[11]又提出了以彈塑性理論為基礎，結合試驗引入合理假設，提出了常見幾類泥石流防治結構的沖擊力計算方法，并給出了相應的計算公式。王強等[12]提出了泥石流防撞墩沖擊力理論計算方法，陳洪凱等[13]提出泥石流兩相沖擊力以及沖擊時間計算方法。

由于陳洪凱提出的泥石流沖擊模型與實際觀測的值比較接近，故本文采用他的模型。

防治結構或岸坡表面單位面積承受的泥石流液相沖擊力與泥石流液相漿體沖擊力、固相顆粒沖擊力等有關，其表達式可以寫為：

P=f(qf,PS,K0)

其顯式為

P=K0(qf+PS)

(1)

式中：P為單位面積承受的泥石流液相沖擊力；K0為沖擊力顯式系數，一般取500～550；qf為泥石流液相漿體沖擊力；PS為固相顆粒沖擊力。

以沖擊形跡為分析對象，運用極限平衡理論建立沖擊形跡的力平衡方程并據此獲得泥石流沖擊時間。

(2)

式中：T為沖擊形跡的形成時間，即泥石流沖擊時間；Sk為坑內表面積；τk為沖擊形跡抗剪強度；Kr為材料承沖系數；Ak為沖擊形跡坑口面積；θ為泥石流沖擊方向與承沖平面的銳交角。

2動力方程的建立

由于攔擋壩的壩體設計比較復雜，有溢流口，排水孔等構造，而且壩體幾何形狀很不規則，加之壩體本身還有自重存在。為了把鋼筋混凝土攔擋簡化成懸臂梁，故引入下面假設。

(1)假設泥石流攔擋壩的壩體受力如懸臂梁，取1m壩寬分析。

(2)假設整個結構在受力過程中始終在線彈性范圍內變化。

(3)不考慮壩體自重對方程建立的影響，即懸臂梁在橫向運動過程中不考慮軸向力的影響。

(4)由于截面是沿著高度連續變化的，假設高度l不變，根據慣性矩等效原則得到相應的矩形截面等效慣性矩以及等效截面。

文獻[14]具體闡述了這一轉換方法，且經計算，劃分成四等分就可以得到足夠的精度。即根據最大位移相等原理，運用辛普生數值積分，將變截面分為四段，各段處的慣性矩分別為I0，I1，I2，I3，I4，由于是梯形截面，所以I0等于零，因此在均布荷載作用下的等效慣性矩可以表示為下式

(3)

梁的理論有兩種，一種是Timshenko梁[15]，一種是Euler梁；Timshenko梁是考慮彎曲和剪切影響，而Euler梁忽略了彎曲和剪切。Euler梁以其形式簡單，便于應用因而在工程上被廣泛采用。根據歐拉伯努利梁的理論可以得到懸臂梁位移的運動方程。

(4)

式中：ω(x,t)為懸臂梁的位移；EIeq為彎曲剛度；ρ為材料的密度；A為橫截面積；c為線性阻尼系數。

初始條件：在初始時刻，位移和速度都為零，即

(5)

邊界條件：在固定端，位移和轉角為零，即

(6)

在自由端，剪力和彎矩為零，即

(7)

3動力方程的解答

方程的解答可以分為兩部分，在泥石流沖擊時間內，鋼筋混凝土攔擋壩作受迫阻尼運動，在泥石流沖擊時間之外作自由阻尼運動。

3.1　自由阻尼運動

即當泥石流沖擊荷載為零，f(x,t)=0時。該方程是線性齊次的偏微分方程，邊界條件是齊次邊界條件，故采用數學物理方程中的變量分量法[16]，將線性齊次偏微分方程化成兩個常微分方程，然后根據結構動力學中的振型疊加法[17]求解。求得解為

ζk[sinh(ukx)-sin(ukx)]}·

(8)

式中

對于工程問題一般k取1，2，3，即用前三個振型，即可滿足精度要求。

3.2　受迫阻尼運動

即當泥石流沖擊荷載為f(x,t)=P，采用特征函數系法[16]，利用攔擋壩作自由阻尼運動中已求的特征函數，經過拉普拉斯變換和拉普拉斯逆變換[18]，即可求得方程的解。即

(9)

式中：

式中：

Xk(x)=cosh(ukx)-cos(ukx)+

ζk[sinh(ukx)-sin(ukx)]

對于工程問題一般k的取法同前。

4動力響應分析

由式(8)已經得到攔擋壩的位移ω(x,t)，又由材料力學中的位移與彎矩、彎矩與應力以及應力與應變關系之間可以得到應力以及應變。

(10)

式中：σ為應力。

(11)

式中：ε為應變。

5算例分析

5.1　工程概況

在甘肅省舟曲縣三眼峪溝泥石流災害治理工程中，攔擋工程共采用了10個重力式泥石流攔擋壩，其中主1號壩的尺寸為：壩頂寬3 m，壩底寬11.25 m，壩高15 m，壩長63.73 m，采用C30混凝土，彈性模量3.0×104N/mm2，配筋按構造配筋，鋼筋混凝土密度2 500 kg/m3，泊松比0.2。沖擊荷載假設采用陳洪凱的沖擊力為7 022 kPa，沖擊時間為8 s。經過計算得EI=6.30×1011N·m2，A=6.33 m2。經過計算可以得到表1中各參數值。

5.2　解析求解

方程的解析解由兩部分組成，即Xk(x)乘以Tk(t)，由于兩者是相乘的關系所以這兩者是獨立的，故分析的時候也把兩部分別求解。采用MATLAB編程計算，當k=1，k=2，k=3分別得到圖1、圖2、圖3。

表1　計算參數

圖1 k=1,Xk(x)和Tk(t)函數圖像Fig.1k=1,ThefunctionimageofXk(x)andTk(t)圖2 當k=2時Xk(x)和Tk(t)函數圖像Fig.2k=2,ThefunctionimageofXk(x)andTk(t)圖3 當k=3時Xk(x)和Tk(t)函數圖像Fig.3k=3,ThefunctionimageofXk(x)andTk(t)

圖中的左邊為Xk(x)函數的圖像，圖中的右邊為Tk(t)函數的圖像。通過對圖1、圖2、圖3中的左圖分析，可以比較得出，當x=15 m，Xk(x)取到最大位移。

將Xk(x)中的最大位移取出，分別為13 mm、7.8 mm、6.4 mm，與對應的Tk(t)函數相乘，然后三項相加，最后通過MATLAB計算得到，當x=15 m，x=10 m，x=5 m處攔擋壩隨時間的位移響應曲線如圖4中虛線所示，采用同樣的方法可以得到其應力響應曲線如圖5中虛線所示。

5.3　數值驗證

用MIDAS NFX軟件進行分析，把鋼筋混凝土材料等效為一種材料，分析參數如表2所示。

采用自動生成3D網格劃分模型，模型的邊界接觸采用固定約束，動力荷載采用時間依存的節點力。

表2　參數特性

圖4　位移動力響應曲線 Fig.4 The displacement of the dynamic response curve

圖5　應力動力響應曲線 Fig.5 The stress of the dynamic response curve

如圖4中實線所示，得到當x=15 m，x=10 m和x=5m時的鋼筋混凝土攔擋壩的位移響應曲線。如圖5中實線所示，得到當x=15 m，x=10 m和x=5 m時的鋼筋混凝土攔擋壩的應力曲線。

通過圖4可以發現，位移動力響應解析解和數值解曲線的波動特征吻合程度一般，但曲線的幅值特征吻合程度比較突出，而且數值解整體比解析解大， 這是由于理論分析的時候，只是分析壩體寬度為1 m的攔擋壩，沒有考慮攔擋壩的空間橫向作用。

通過圖5同樣可以發現，可以得出解析解應力響應曲線與數值解應力響應曲線波動特征吻合程度一般，但曲線的幅值特征吻合程度比較突出，而且數值解整體比解析解大。還可以發現應力響應曲線在某些點出現了交叉，這說明解析解得出的結果與數值解得出的結果有一定的吻合程度。

6結論

(1)將攔擋壩簡化成了懸臂梁，選取了典型的泥石流沖擊荷載，建立了荷載與位移之間動力偏微分方程，得到了解析解。

(2)通過解析求解和數值分析分別得到鋼筋混凝土攔擋壩的位移響應曲線和應力響應曲線，通過對兩者的對比，發現曲線波動特征吻合程度一般，但曲線的幅值特征吻合程度比較突出，表明了這種分析方法是合理可靠的。

(3)由于采用了結構動力學方法分析，攔擋壩會產生位移，需要具有一定的變形能力，故本文分析方法僅適用于鋼筋混凝土攔擋壩，不適用于毛石壘砌的重力式攔擋壩。

參考文獻

[1]DZ/T239-2004.泥石流災害防治工程設計規范[S].北京：中國地質調查局，2004.

[2]吳瑋江，馮樂濤. 重力式泥石流攔擋壩穩定性計算[J].甘肅科學學報, 2012, 24(3): 56-59.

WU Wei-jiang, Feng Le-tao. Calculation of the stability of the gravity dam in front of debris flows [J]. Journal of Gansu Sciences, 2012, 24(3): 56-59.

[3]喬振華. 淺論泥石流地區攔擋壩[J]. 山西建筑, 2005, 31(13): 59-60.

QIAO Zhen-hua. Dam theory of shallow debris flow area [J]. Shanxi Architecture, 2005, 31(13), 59-60.

[4]王曉明. 泥石流防治中的攔擋壩組合設計[J]. 中國水運, 2012, 12(21): 263-264.

WANG Xiao-ming. Debris flow prevention and control of the dam combination design[J]. China Water Transport, 2012, 12(21): 263-264.

[5]熊道錕，徐世民. 泥石流攔擋壩之虞[J]. 中國地質災害與防治學報, 2010, 21(4): 136-138.

XIONG Dao-kun,XU Shi-min. The danger of dam against debris flow[J]. The Chinese Journal of Geological Hazard and Control,2012, 12(21): 263-264.

[6]李冠奇，魏鴻，谷明成，等. 泥石流攔擋壩優化設計[J]. 地質災害與環境保護, 2009, 20(2): 29-32.

LI Guan-qi, WEI Hong, GU Ming-cheng, et al. Dam against debris flow optimization design[J]. Journal of Geological Hazards and Environment Preservation, 2010, 21(4): 136-138.

[7]賈世濤，崔鵬，陳曉清，等. 攔沙壩調節泥石流攔擋與輸移性能的試驗研究[J]. 巖石力學與工程學報,2011,30(11): 2338-2345.

JIA Shi-tao, CUI Peng, CHEN Xiao-qing, et al. Experimental study of regulating barrage and transportation properties of debexperimental study of regulation barrage and transportation properties of debris flow by silt-trap dam[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(11): 2338-2345.

[8]王根龍，張茂省，張新社，等. 舟曲“8．8”特大泥石流治理中的關鍵技術研究[J]. 西北地質, 2011, 44(3)：100-105.

WANG Gen-long, ZHANG Mao-sheng, ZHANG Xin-she, et al. Zhouqu to heavy “8.8” debris flow governance research of key techniques[J]. Northwestern Geology, 2011, 44(3): 100-105.

[9]何曉英,唐紅梅,朱繡竹,等.泥石流漿體沖擊特性實驗研究[J].振動與沖擊,2013,32(24):127-134.

HE Xiao-ying, TANG Hong-mei,ZHU Xiu-zhu,et al.Tests for impacting characteristics of debris flows slurry[J].Journal of Vibration and Shock,2013,32(24):127-134.

[10]何思明，李新坡，吳永. 考慮彈塑性變形的泥石流大塊石沖擊力計算[J]. 巖石力學與工程學報, 2007, 26(6): 1664-1669.

HE Si-ming, LI Xin-po, WU Yong. Calculation of impact force of outrunner blocks in debris flow considering elastoplastic deformation[J].Chinese Journal of Rock Mechanics and Engineering, 2007, 26(6): 1664-1669.

[11]何思明，吳永，沈均. 泥石流大塊石沖沖擊力的簡化計算[J].自然災害學報, 2009, 18(5): 51-56.

HE Si-ming, WU Yong, SHEN Jun. Simplified calculation of impact force of massive stone in debris flow[J]. Journal of Natural Disasters, 2009, 18(5): 51-56.

[12]王強，何思明，張俊云. 泥石流防撞墩沖擊力理論計算方法[J]. 防災減災工程學報, 2009, 29(4): 423-427.

WANG Qang, HE Si-ming, ZHANG Jun-yun.Mudslides anti-collision pier impact theory calculation method[J]. Journal of Disaster Prevention and Mitigation Engineering, 2009, 29(4): 423-427.

[13]陳洪凱，唐紅梅. 泥石流兩相沖擊力及沖擊時間計算方法[J]. 中國公路學報, 2006, 19(3):19-23.

CHEN Hong-kail, TANG Hong-mei. Method to calculate impact force and impact time of two-phase debris flow[J]. China Journal of Highway And Transport, 2006, 19(3)：19-23.

[14]高軒能，朱皓明，沈劍波.截面連續變化懸臂構件撓度計算的等效柱法[J].數學的實踐與認識,2011,41(8):67-72.

GAO Xuan-neng, ZHU Hao-ming, SHEN Jian-bo. Cross section of continuous variation cantilever components deflection calculation of equivalent column method[J].Mathematiocs in Practics and Theoy, 2011,41(8):67-72.

[15]Timoshenko S, Young D H, Weaver W. Vibration problems in engineering [M]. 4th edition, John Wiley & Sons, Inc. 1974.

[16]陳才生，李剛，周繼東，等. 數學物理方程[M]. 北京：科學出版社，2008.

[17]Roy R. Craig，Andrew J. Kurdila. Fundamentals of structural dynamics [M]. 3th edition John Wiley & Sons Inc. 2006.

[18]嚴鎮軍. 數學物理方法[M]. 合肥：中國科學技術大學出版社，2005：169-197.