ABAQUS模擬滲流要點

1、ABAQUS模擬滲流要點(2011-05-24 22:56:10)轉(zhuǎn)載標(biāo)簽：雜談從其他網(wǎng)站看到一些這樣的介紹，希望對飽和土方面模擬的朋友有幫助：當(dāng)只進行滲流計算時：1、由于Abaqus中缺乏非耦合的孔壓單元，這時可采用耦合單元，但要約束住所有位移的自由度。2、滲流材料參數(shù)選擇。在CAE中都是在(Material-creat-other-pore fluid)選項中。(1)Gel：定義凝膠微粒吸濕膨脹的發(fā)育過程，這在一般的巖土分析中應(yīng)用不多。(2)Moisture swelling：定義由于吸濕飽和所引起的固體骨架體積膨脹(或負(fù)吸力引起的骨架收縮)。(3)Permeability：定義飽和介質(zhì)的

2、滲透系數(shù)，該滲透系數(shù)可以在type選項中定義為各向同性、正交各向異性和各向異性，并且可以根據(jù)Void Ratio定義為孔隙比的函數(shù)。在Suboptions中選擇Saturation Dependent 參數(shù)來指定與飽和度相關(guān)性系數(shù)ks(s)，缺省設(shè)置為kss3，而非飽和介質(zhì)滲透系數(shù)k=ksk。選擇Velocity dependence參數(shù)可以激活Forchheimer定律，缺省的是Darcy定律。(4)Pore Fluid Expansion：定義固體顆粒與流體體積熱變化效應(yīng)。(5)Porous Bulk Moduli：定義固體顆粒與流體體積模量。(6)Sorption：定義負(fù)孔隙壓力與飽和度

3、之間的相關(guān)性。當(dāng)type=Absorption時，定義吸濕曲線，type=Exsorption時定義排水曲線。3、載荷及邊界條件(1)通過(Load-creat-step-fluid-surface pore fluid)選項定義沿著單元表面的外法線方向的滲流速度vn，當(dāng)考慮降雨影響時可采用此載荷(2)邊界條件(Boundary condition-creat-other-pore pressure)選項定義孔壓邊界條件，此時要先假定浸潤面的位置，然后定義浸潤面上的孔壓為零，Abaqus會在后續(xù)的分析計算中自動計算出浸潤面的位置。Abaqus默認(rèn)的是不透水邊界。(3)當(dāng)滲流自由面遇到臨空的自由

4、排水面時，需要定義一個特殊的邊界條件。此時可以通過在inp文件中加入*Flow或*Sflow來定義(4)初始條件的定義。初始條件中一般要定義以下幾種：*initial condition,type=saturation初始飽和度*initial condition,type=pore pressure初始孔壓*initial condition,type=ratio初始孔隙比當(dāng)進行耦合分析時，基本步驟同上，但要去掉除邊界條件之外的約束，同時還要在邊界上加上流體壓力。*這些問題都是很有意義的問題和很關(guān)鍵的技術(shù)，一下子討論和理解的人也不是很多，如能給出一些例子來進一步闡述和說明（越詳細(xì)越好），我想

5、可以作為精華帖。*likuiming編輯上午只看了一個damseep.inp文件，大壩滲流。只是了解其中每個命令的含義，沒有運行看其結(jié)果。*MATERIAL, NAME=soil*DENSITY2000.0,0.0*ELASTIC, TYPE = ISOTROPIC1000.0E9,0.2,0.0*E,r,溫度*PERMEABILITY，TYPE=ISOTROPIC（默認(rèn)），ANISOTROPIC，SATURATION，（飽和）VELOCITY（速度）*PERMEABILITY,SPECIFIC=10000.0*流體比重specific weight（FL3）*2.11E-4,*滲透系數(shù)K，空

6、隙率e,溫度*SORPTION*定義材料吸附作用TYPE=ABSORPTION (default) *-.0,0.04*孔隙壓力（負(fù)壓），飽和度s（0.01=s=1.0）*-.0,0.050.0,1.0*INITIAL CONDITIONS,TYPE=SATURATION*定義初始飽和度，默認(rèn)為1*ALLNODE,1.0*節(jié)點集或節(jié)點號，飽和度（默認(rèn)為1）*定義初始孔隙水壓力（initial pore fluid pressures）for a coupled pore fluid diffusion/stress analysis*INITIAL CONDITIONS,TYPE=PORE P

7、RESSUREPOR,.,0.0,0.0,20.*節(jié)點集，流體第一孔隙壓力，垂直向第一孔隙壓力，流體第二孔隙壓力，垂直向第二孔隙壓力*POR11,0.0POR12,0.POR13,0POR14,0.POR15,0.*INITIAL CONDITIONS,TYPE=RATIO*定義初始孔隙比*ALLNODE,1.0*節(jié)點，第一孔隙比，垂向，第二，垂向*BOUNDARYALLNODE,1ALLNODE,2*RESTART,WRITE*STEP,INC=100,UNSYMM=YES*SOILS,UTOL=0.2*調(diào)用自動時間步*0.2,1.,1.0E-5,0.3*初始時間增量，計算總時間，最小，最大時間增量，孔隙壓力隨時間變化率*DLOADDAM,GRAV,9.8,0.0,-1.,0.*BOUNDARY66,8,0.067,8,19600.68,8,39200.69,8,58800.70,8,78400.71,8,98000.72,8,.73,8,.74,8,.75,8,.6,8,.26,8,58800.027,8,39200.28,8,19600.29,8,0.0*FLOW*定義滲透系數(shù)和孔隙壓力？*FSIDE,Q2D,0.01*單元集，滲流類型，起始孔隙壓力，起