粘性泥石流流量的研究

粘性泥石流流量的研究

1泥石流體作用力粘性泥石流流動的形式主要是矩陣流，其脈沖波形為鋸齒狀。在流動中，脈沖波表示為幾個三角形脈沖，每個脈沖的作用時間非常短?，F(xiàn)在應用于計算泥石流沖擊力的公式主要有:(1)C.M.弗萊施曼公式Ρ動=Κγcαv22g(1)P動=Kγcαv22g(1)(2)伊茲巴什與哈爾德拉公式Ρ動=Κγcv22g(2)P動=Kγcv22g(2)(3)И.И.赫爾赫烏利澤公式Ρ總=0.1γc(5Η0+v2)(3)P總=0.1γc(5H0+v2)(3)上述三式中P動為泥石流對攔擋壩的沖擊動水壓強,T/m2;P總為攔擋壩承受的靜水壓強和動水壓強之總值,T/m2;v為擋壩處的泥石流流速,m/s;γc為泥石流體平均密度,T/m3;g為自由落體加速度,m/s2;K為由互撞體特征決定的系數(shù),(1)式取2.0,(2)式取1.3;α為由動量所決定的動量改正系數(shù),1.33;H0為泥石流泥深,m。(4)章書成,袁建模公式σ=ΚρV2(4)σ=KρV2(4)式中σ為建筑物單位面積作用的沖擊動壓力,N/m2;K為表征泥石流體不均值的系數(shù),用于類同蔣家溝泥石流時,取3.0;ρ為泥石流體密度,kg/m3;V為泥石流體流速,m/s。以上公式的共同之處在于,均假定泥石流體為連續(xù)體,同一流動面上各點的流速和密度是一致的,利用泥石流體流速和密度來推算泥石流沖擊力,但是粘性泥石流的運動狀態(tài)表現(xiàn)為陣流,而并非連續(xù)流,泥石流流動狀態(tài)復雜,其內(nèi)部瞬時速度,流速分布規(guī)律未知,從觀測數(shù)據(jù)上可知,沖擊力也是時時變化的,直接利用上述公式,需同時對泥石流密度和流速兩個變量進行假設,計算結(jié)果與實際情況之間的誤差較大。但如果根據(jù)實測的沖擊力波形簡單取最大值或是平均值作為工程計算參數(shù),也是不合理的。本文提出基于動量守恒原理,采用Hilbert變換對實測數(shù)據(jù)進行信號處理,通過計算求得泥石流體瞬時沖擊力代表值,作為防治工程的設計依據(jù)。2受泥石流沖擊的建筑結(jié)構微元沖擊行為首先引入泥石流流體微元的概念,即沿著泥石流流動方向,按照能夠在攔擋物上產(chǎn)生一個三角形脈沖形式?jīng)_擊力的特征,將泥石流體劃分為若干個流體微元,微元體的面積等于結(jié)構受沖擊作用面積。假定每一流體微元在運動過程中,具有相同的速度和相同的密度,當微元撞擊到建筑結(jié)構上時,其內(nèi)部各點速度同時降為零。對于受泥石流沖擊作用的建筑結(jié)構(建筑物、橋墩、攔擋壩等),在受到泥石流流體微元瞬間作用時,存在著動、沖量守恒mv=Ρ(5)式中m為微元體質(zhì)量,m=ρAlv,A為沖擊力作用面積,即微元體面積,而l=vT為微元長度,T為沖擊力作用時間;P為微元產(chǎn)生的沖量,Ρ=Τ∫0f(t)dt,f(t)為沖擊力曲線函數(shù)。式(5)經(jīng)變換后,得到σ=Τ∫0f(t)dtAΤ=ρv2(6)式(6)可以改寫為F=σA=Τ∫0f(t)dtΤ(7)σ為結(jié)構單位面積上作用的沖擊力,即沖擊壓強;F為沖擊壓力。從(6)式可知,如果可獲得微元沖擊力作用時間,以及在該段時間內(nèi)的作用沖量,利用已知的力作用面積即可求得微元的沖擊壓強。3ft變換泥石流沖擊力作用下的動態(tài)反應信號屬于非平穩(wěn)信號,即統(tǒng)計量(相關函數(shù)、功率譜等)隨時間變化的信號,因此時頻分析方法是分析這類信號的有效手段,Hilbert變換即數(shù)字濾波器中一個重要方法,它是1946年伽波(Gobor)提出的。信號f(t)的Hilbert變換,記為?f(t),是指?f(t)=1πt×f(t)=1π∞∫-∞f(τ)t-τdτ(8)Hilbert變換可以看成是一個積分濾波器,又稱為90°相移濾波器或垂直濾波器,信號經(jīng)過Hilbert變換后產(chǎn)生90°相移。采用Hilbert變換的本質(zhì)是把實測得到的所有有效信號移相90度,由f(t)+i×Hilbert(f(t))組成的復信號的幅值正好精確等于窄帶調(diào)幅信號的包絡,給定信號的包絡為A(t)=√[f(t)]2+[?f(t)]2。利用Hilbert變換,對沖擊力數(shù)據(jù)進行移相處理,求出信號的包絡線,利用包絡線求得沖擊力作用時間,如圖3所示。由于實測沖擊力波型圖上沖擊力開始和停止的時間不易確定,做如下規(guī)定:若記錄數(shù)據(jù)包絡線的最大振幅為A,則從振幅αA開始到波型衰減到αA為止的一段時間作為振動持續(xù)時間T,通常取α=(15～13),本文取=α=15,在T范圍內(nèi)對沖擊力信號做積分,就可以得到?jīng)_量值。4角形脈沖形式的計算結(jié)果在兩個國家自然科學基金項目中采用上述計算方法對實測數(shù)據(jù)進行了分析計算。本文選擇了其中一組計算數(shù)據(jù),圖4為其中一個有代表性的實測微元沖擊力的三角形脈沖形式。計算結(jié)果見表1。從表中可見,除了個別值外,波形中最大沖擊力值與計算沖擊壓力值的比值在3-4之間。5建立泥石流沖擊檢測設施研究泥石流體沖擊力的不穩(wěn)定性,源自于流體內(nèi)攜帶的大量規(guī)格不同的固體物質(zhì),水流和固體物質(zhì)對結(jié)構的撞擊過程不同,引發(fā)的結(jié)構動態(tài)響應歷程不同,這些問題還有待于進一步深入研究。在現(xiàn)有條件下,可以在泥石流溝谷中建立沖擊力測試站,或可以直接利用建在溝谷中的橋墩、攔擋壩等構筑物安裝沖擊力傳感器,采集泥石流沖擊作用歷程信息,通過數(shù)據(jù)處理,利用公式(6)和