擾動淤泥與沉積淤泥的流變特性研究(1)程力學(xué)

《擾動淤泥與沉積淤泥的流變特性研究(1)程力學(xué)》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、擾動淤泥與沉積淤泥的流變特性研究(1)程力學(xué)摘要：為了分析不同沉積狀態(tài)的淤泥流變特性，本文將淤泥作為線性粘彈性體，提出了根據(jù)泥水系統(tǒng)控制方程解答和波浪在泥床上傳播時的實測變化規(guī)律求取淤泥流變參數(shù)的反分析方法，并利用遺傳算法分別求解攪勻(擾動)與自然沉積淤泥的流變參數(shù).結(jié)果表明，淤泥平均密度相同時，自然沉積淤泥的彈性模量明顯大于人工攪勻淤泥的彈性模量，兩者的粘性系數(shù)隨淤泥密度的變化規(guī)律也有明顯區(qū)別，這主要是由于沉積淤泥與攪勻淤泥的內(nèi)部結(jié)構(gòu)不同所造成的.關(guān)鍵詞：自然沉積淤泥擾動淤泥流變特性波浪傳播在河口海岸粘性泥沙運動研究中，淤泥的流變特性和許多物理過程密切相關(guān)，如波浪在淤泥質(zhì)海床上的衰減特性和波

2、浪作用下的泥床質(zhì)量輸移、淤泥液化等[1～3]，是影響波浪與底泥床相互作用的主要因素.淤泥流變特性的變化和水動力作用下床面穩(wěn)定性、底部泥沙的懸揚規(guī)律等也有著直接和間接聯(lián)系.實驗和計算結(jié)果表明，淤泥處于不同的沉積狀態(tài)時，由于流變參數(shù)的變化，波浪在泥床上傳播時波高衰減率可以相差1～2個量級[4].因此，了解淤泥的流變特性是深入研究淤泥質(zhì)海岸泥沙運動規(guī)律的關(guān)鍵問題之一.關(guān)于淤泥的流變關(guān)系，目前已建立了相當(dāng)多的模型[4]，然而由于模型參數(shù)特別是現(xiàn)場條件下的模型參數(shù)難以確定，使得各種模型在應(yīng)用于描述現(xiàn)場條件下波浪和淤泥質(zhì)海床的相互作用規(guī)律時受到限制.由于現(xiàn)有的測量儀器，如旋轉(zhuǎn)同心圓筒流變儀等實際上只能測量

3、擾動淤泥的流變特性，對于自然沉積的底部泥床，目前還沒有很好的方法來測定其流變參數(shù)，對其流變特性的變化規(guī)律有待深入研究.基于以上原因，本文假定淤泥作為線性粘彈性體，在文獻(xiàn)[5]的基礎(chǔ)上，采用根據(jù)實驗結(jié)果反求模型參數(shù)的方法，比較了實驗室內(nèi)人工攪勻和自然沉積淤泥的流變特性，討論了淤泥流變特性變化對波浪衰減規(guī)律的影響.1理論分析1.1關(guān)于淤泥特性的假設(shè)近年來一系列的研究表明，淤泥在波浪作用下表現(xiàn)出復(fù)雜的非線性粘彈性體特征，其流變參數(shù)是應(yīng)變或應(yīng)變率歷史的函數(shù).但為了簡化問題，通過反分析方法確定淤泥的流變參數(shù)，并比較人工攪勻淤泥與沉積淤泥流變特性，這里仍然假定淤泥作為線性粘彈性體.在應(yīng)變較小的情況下，上述

4、假設(shè)精確地描述了淤泥的運動規(guī)律；在應(yīng)變較大的情況下，這種假設(shè)相當(dāng)于對淤泥的本構(gòu)方程進(jìn)行等價線性化.根據(jù)線性粘彈性體假設(shè)，在頻率為σ的循環(huán)(振蕩)荷載作用下，淤泥的本構(gòu)關(guān)系可表示為[6]：對于天然或?qū)嶒炇覂?nèi)形成的自然沉積泥床，沿泥床表面向下的泥密度分布及其流變特性都是變化的，為了簡化問題而方便估計沉積淤泥的流變特性，我們假定泥床具有均勻密度和流變參數(shù)，水波與攪勻或沉積泥床的相互作用就都可以用密度均勻的兩層介質(zhì)線性系統(tǒng)模型來描述.1.2泥—水系統(tǒng)運動控制方程根據(jù)maa和mehta[7]，波浪在泥床上傳播時，上層水體和下層淤泥運動的連續(xù)性方程和線性化運動方程可表示為：假定水面和泥水交界面位移表達(dá)式為

5、：ηj(x,t)=ajexp[i(kx-σt)](5)式中：k=kr+iki為待求的復(fù)波數(shù)，實部與波長l的關(guān)系為kr=2π/l，虛部ki表示波浪衰減率，aj為第j層位移振幅.方程(2)～(4)中uj,uller方法求解.1.3淤泥流變參數(shù)的反求法在淤泥流變參數(shù)確定的情況下，根據(jù)方程(11)可以求出波浪在泥床上傳播時的復(fù)波數(shù)，然后可進(jìn)一步得到速度和壓力分布，這是分析波浪和泥床相互作用問題的一般過程.同樣，如果已知波浪在泥床上傳播的復(fù)波數(shù)而淤泥流變參數(shù)未知，由方程(11)也可以反求出淤泥的表觀復(fù)粘度，從而得到淤泥的彈性系數(shù)和粘性系數(shù).按照把淤泥作為線性粘彈性體的假設(shè)，對于同一種淤泥，泥密度保持不變

6、時，由不同波浪條件的實驗結(jié)果反求出的淤泥粘彈性系數(shù)值應(yīng)該相同.但實際上，由于淤泥的非線性性質(zhì)，其粘彈性特征隨受力狀態(tài)而變化[8，9]，再加上實驗測量誤差的限制，根據(jù)不同實驗條件反求出的淤泥粘彈性系數(shù)是不會相同的.為了求得假定淤泥為線性粘彈性體時代表性的粘彈性參數(shù)值，我們根據(jù)每組淤泥的多組實驗結(jié)果構(gòu)造優(yōu)化問題來反求參數(shù).設(shè)對于某一密度的淤泥共進(jìn)行了m組規(guī)則波浪在淤泥床上傳播的實驗.對于第m(m=1,…，m)組實驗，已知波浪周期t、波高h(yuǎn)、水深d1、水密度ρ1、水的粘性系數(shù)μ1、泥厚d2、泥密度ρ2,根據(jù)實驗測量結(jié)果可求得復(fù)波數(shù)kem.假定粘彈性參數(shù)已知，則根據(jù)式(11)可以求得第m組實驗條件下的

7、復(fù)波數(shù)kcm.由于不存在統(tǒng)一的μe滿足所有實驗條件下(m=1,…，m)的式(11)，我們構(gòu)造一個極值問題來確定最具代表性的淤泥流變參數(shù)值.令：式中：kcm為μe的函數(shù)，則φ最小時的μe值即可認(rèn)為淤泥作為線性粘彈性體時的代表值.文獻(xiàn)[5]采用有約束的外部懲罰函數(shù)方法來反求μe值，存在對初值的依賴性高，最后收斂解難以確定等問題.這里我們采用遺傳算法[10]來求解式(12)，最終得到的μe值與初值無關(guān).