基于地質統(tǒng)計學的土石壩滲流分析

第一章緒論 21.1研究背景 21.2國內外研究進展 21.3研究內容及技術路線 4第二章：土石壩滲流基本理論及ABAQUS軟件 42.1土石壩滲流隨機理論 42.2滲流基本理論 52.3滲流基本方程 62.4滲流問題的邊界條件和初始條件的確定 72.5滲流的有限元分析方法 82.6ABAQUS有限元軟件 92.7滲流應力耦合基本原理及其在ABAQUS中的應用 2.8ABAQUS二維均質土壩的穩(wěn)定滲流分析案例 第三章：研究內容 3.1邊界條件的確定 3.2壩體滲流分析數(shù)值模擬 3.3降雨入滲結果分析 第四章：總結和展望 4.1結論 4.2展望 參考文獻 第一章緒論土石壩是水工建筑物中最常見的一種，不僅僅是世界上最古老的一種壩型。也是世界大壩工程的建設中的最廣為人知的一種壩型。它集應用最為廣泛且發(fā)展最快等優(yōu)點于一身，基于土石壩可以就地取材的優(yōu)點，所以能夠大大節(jié)約重要的建筑材料例如鋼材﹑水泥﹑木材等等，因為其構造較為簡單，所以利于進行維修處理，以及加高并擴建的工作，因此經(jīng)常作為擋水建筑物來應用于各項水利工程中。但因為其壩體土石材料的限制，壩身不能溢流，由于滲流的影響，使得壩體使壩體滑移，另一方面，由于滲流問題持續(xù)嚴重，壩體內的土石材料自身容重也會產(chǎn)生變化，致使土體顆粒內部的有效應力的急劇減小現(xiàn)象，從而使壩體的抗剪強度也降低，嚴重時甚至導致土石壩的潰壩。土石壩作為一種擋水建筑物，由于發(fā)生的一系列土石壩失事事故往往都是滲流破壞的產(chǎn)物。論：保證土石壩安全穩(wěn)定的前提條件是要對水體的滲流進行一定的控制。潰壩的是一個巨大的威脅隱患。造就邊坡失穩(wěn)的因素很多，這其中的幾何條件包括邊坡有斷層，裂隙發(fā)育水平等邊坡裂隙要素，面對這以上種種因素，最重要的還是水有效解決，因此研究土石壩滲流及控制安全問題就顯得格外必要。驗法等等。最先研究土石壩滲流影響工作的方法代表主要為水力學法和圖解法。這兩種方法都是基于一種科學的模擬，這種方法主要以達西定律為基礎，通過假設某一種情況來進行分析計算?；谕潦瘔螡B流分析問題情況的復雜性，邊界條件的不同，壩體材料的差異等使得通過水力學法來計算土石壩的滲流問題基本是個不可能事件，且伴隨著計算機信息水平的迅猛發(fā)展，涌現(xiàn)了大批的科學方法，人們最早是借助有限元的方法來分析土壩中與之相關的問題，最先應用于滲流問向異性的飽和土體滲流-耦合的有限元公式。在滲流計算分析工作中，最大的重點進行處理，所以大多數(shù)情況下是先假定一個自由面，然后不斷地迭代，使其逼近原理,提出了自己對于地基滲透破壞現(xiàn)象包括其中的發(fā)生原理和影響因素的獨特見過研究砂卵石在滲透變形方面上的形成條件,并在計算方法的創(chuàng)新上，創(chuàng)造性的提出了臨界坡降計算方法的這一理念，進一步對滲透變形的形式做了一個總結。葛說明。我國在滲流分析的數(shù)值模擬方面，成功開發(fā)了一些較為成熟的軟件，這其中包括有南京水利科學研究院的土壩滲流計算軟件UNSST2，河海大學工程力學所三.搜集土石壩的滲流基礎理論資料，并進行有效分類，著重研討數(shù)值模擬方法在土石壩滲流工作中的適用條件及應用范圍。了解三種基本數(shù)值理論，特別是有第二章：土石壩滲流基本理論及ABAQUS軟件以推測巖土體各種性質都具有空間變異性，研究表明，這種不確定性并不等同于和各向同性，所以不能夠表現(xiàn)出空間局部范圍和特定方向上的巖土參數(shù)的變化特征及其規(guī)律?？紤]到巖土成因條件及各項要素都會對巖土的性質或多或少產(chǎn)生一定的作用，所以用隨機場的模型理念來反映巖土參數(shù)的變化。為條件的干涉，上下游的水位不可能固定在某一個位置不發(fā)生變化，而對于滲透系數(shù)的變化就更多了，壩體的任一位置的滲透系數(shù)都有著差別，而這些重要參數(shù)主要是根據(jù)它的數(shù)學參數(shù)例如分布類型，均值，方差，協(xié)方差等來進行描述。而Q=KA（h2-h1）LJ=（h2-h1）LA根據(jù)（2-12-22-3）推導出（2-4）達西定律揭示出了土體中滲流多數(shù)的土例如砂土黏土中的滲流問題都能應用層流規(guī)律進行解決，因此對于簡化計算來說，可以應用達西定律處理。達西定律是進行所有的滲流分析工作求解的FIG.2-1bodyofdifferentialelement?z?z?x?y?z?x?y?z-?x?y?z?x?y?t?y?z?x滲流場根據(jù)運動要素是否與時間有關劃分為兩部分：穩(wěn)定滲流場和非穩(wěn)定滲2.3.2.1穩(wěn)定滲流基本方程 ?(Kxx?h+Kxy?h+Kxz?h)+?(Kyx?h+Kyy?h+Kyz?h)+?(Kzx?h+Kzy?h+Kzz?h)=0?x?x?y?z?y?x?y?z?z?x?y?z上述為普通土體的穩(wěn)定滲流基本方程，若假設土2.3.2.2非穩(wěn)定滲流基本方程 ?(Kxx?h+Kxy?h+Kxz?h)+?(Kyx?h+Kyy?h+Kyz?h)+?(Kzx?h+Kzy?h+Kzz?h)=S?h（2-7）?x?x?y?z?y?x?y?z?z?x?y?z?t此為普通土體的非穩(wěn)定滲流場的基本方程。釋放出來的儲水量。單位為??1?x?x?x?t表達式：〔（Kxx?h+Kxy?h+Kxz?h）lx+(Kyx?h+Kyy?h+Kyz?h)ly+(Kzx?h+Kzy?h+Kzz?h)lz〕?x?y?z?x?y?z?x?y?z?nk?h?nα,β：邊界點上的已知數(shù)有限元法在求解滲流時，會將待求解的滲流區(qū)劃分為多個子區(qū)域，這些子區(qū)域之間能夠相互連接，然后再用這些子區(qū)域內連續(xù)的分區(qū)近似水頭來替代待定的情況復雜的求解域進行模型化處理。其中各子區(qū)域內的近似函數(shù)一般通過未知場以由已經(jīng)知道的若干內插值函數(shù)通過迭加的方法求得[12]。而這些待定系數(shù)就可以近似的為單元邊界各節(jié)點的水頭函數(shù)，相互連接的有限元他們所有的連接節(jié)點上求滲透區(qū)域內的節(jié)點水頭函數(shù)就反映了該滲透區(qū)域整體的水頭函數(shù)。其求解實質就是講微分求解變成有限個待求值的代數(shù)方程，未知量求出后就可以利用插值函數(shù)求得各有限單元的滲流場函數(shù)近似值，從而求得整個滲流場的函數(shù)解[13]。在求解過程中，有限元的數(shù)量越多或者是有限單元的自由度越多那么近似解就越接近有限元法的求解基礎為變分泛函數(shù)積分，把微分方程和邊界條件轉變成線性方程組再求得泛函數(shù)極值，從而求得滲流結果。運用有限元法對滲流場中的水頭：［：［］：①把求解偏微分方程轉變?yōu)榍蠼庾兎?。①把求解偏微分方程轉變?yōu)榍蠼庾兎?。②離散化研究區(qū)域，劃分為若干個有限單元，并確定插值函數(shù)。③利用各有限單元插值函數(shù)逼近變分水頭函數(shù)求得代數(shù)方程組，從而獲得單元滲④將各個子區(qū)域的滲透矩陣整合為總體的滲透矩陣，并聯(lián)系定解條件求得總體的⑤求解得出的線性方程組，獲得各節(jié)點未知水頭值和其它的滲流要素。實際工程中的土石壩滲流問題幾乎都隸屬于三維滲流問題，因此研究土石壩滲流的前提條件為研究三向空間滲流問題。三向空間問題的基本有限元單元可以元法求解土石壩的滲流問題時，滲流自由面的確定是重中之重，因此，如何確定②作為一個合格的巖土工程分析軟件，它必須具備能夠處理復雜邊界及載荷條件模擬填土或開挖造成的邊界條件改變，值得一提的是，該軟件還提供了無限元，模擬填土或開挖造成的邊界條件改變，值得一提的是，該軟件還提供了無限元，③巖土工程數(shù)值分析步中須將初始應力的作用考慮在內，而該軟件為此專門提供了相應的分析步，這樣就可以靈活又精準的打建初始應力狀態(tài)。④土體是典型的三向體，就目前的理論而言，學者們普遍認為土體的強度和變形土體的本構模型可以反映出土體的力學性能，因此具有重要作用。土體的本利用兩個材料參數(shù)：應力E和應變V就能描述。局限性在于該模型并不能反映土鄧肯-張模型的建立的基礎是土體的三軸實驗，康納根據(jù)三軸實驗得到了土體a+bεaEt=d(?1—?3)=a?ε1(a+bε1)2摩爾-庫倫模型把土體的極限平衡理論作為基礎，認為土體破壞時土體的抗剪強度?1—?3?1—?32.7.2.1應力場影響下二維穩(wěn)定滲流數(shù)學模型: ?x?x?y?y?n2?n2?n3-?n32.7.2.2滲流場影響下的應力場數(shù)學模型［K］[△δ]=[△F］+[△Fs］;(x,y)∈Ω2.7.2.3滲流場和應力場耦合的數(shù)學模型［K］[△δ]=[△F］+[△Fs］2.8.2.1模型建立與求解/【Create】，創(chuàng)建一個材料，并命名為/【PoreFluid】/【Permeability】選項中定義飽和時的滲透系數(shù)為/【PoreFluid】/【Sorption】選項中僅定義吸濕曲線。/【Section】命令，把材料賦給相應的區(qū)域。圖2-3計算終止時的孔壓等值線圖FIG.2-3calculationofporepressurecontourmapattermination圖2-5網(wǎng)格積分點流速矢量圖FIG.2-5flowvelocityvectordiagramofgridintegrationpoints圖2-6飽和度的等值線圖FIG.2-6contourmapofsaturation圖2-7下游無水時浸潤線圖2-11滲透系數(shù)隨機無下游水位浸潤線圖2-12滲透系數(shù)隨機下游水位1.5米浸圖2-13滲透系數(shù)隨機下游水位3米FIG.2-13permeabilitycoefficientofrandomdownstreamwaterlevel3m第三章：研究內容浸潤面以下的部分為飽和土，浸潤線以上的部分屬于非飽和土。本文采用固定網(wǎng)圖3-1邊界條件示意圖FIG.3-1schematicdiagramofboundaryconditionsS2壩體與水平基準面交接觸，由于該邊界隸料參數(shù)如下所示：粘聚力c=15kPa，內摩擦角φ=30°,彈性模量E=100Mpa，泊松圖3-2算例的有限元表格FIG.3-2finiteelementtableofanexampleKw=awKws/[aw+（bw*(ua-uw）^cw]3.2.3.1初始狀態(tài)結果分析FIG.3-3porepressuredistributionbeforerainfall圖3-4降雨之前的飽和度分布Figure3-4saturationdistributionbeforerainfallFIG.3-5stressdistributionbeforerainfallH