1、流體與動態計算實例分析-堤壩 ;流體與動態計算實例分析-大壩newconf dyn fluid;設置動態與流體算法set dyn off fluid off;關閉動態與流體算法 generate foundation and embankment grids and attach interfaces;生成基礎和堤壩網格并粘貼接觸面gen zone bri p0 0,0,0 p1 40,0,0 p2 0,10,0 p3 0,0,10 size 20 5 5;基礎gen zone bri p0 22,0,10 p1 40,0,10 p2 22,10,10 p3 22,0,20 size 9,5,

2、5;堤壩;在z=10,x=1522,y=010上生成接觸面1interface 1 face range x 15.0 22.0 y 0.0 10.0 z 9.9 10.1;在x=22,z=1021,y=010上生成接觸面2interface 2 face range x 21.9 22.1 y 0.0 10.0 z 10.0 21.0設置最大的邊界長度為1.0interface 1 maxedge 1.0interface 2 maxedge 1.0generate block wall 生成擋水墻gen zone bri p0 15,0.5,11 p1 21,0.5,11 p2 15,9.

3、5,11 p3 15,0.5,20.9 siz 3 5 5. ; name groups and move block on to soil命名群組并移動擋水墻到土壤上, group block range x=16,22 y=0,10 z=10,20group soil embankment range x=22,40 y=0,10 z=10,20group dense soil foundation range x=0,40 y=0,10 z=0,10擋墻沿x方向增加一個單位ini x add 1.0 range group block擋墻沿z方向增加一個單位ini z add -1.0

4、range group block assign models to group;給群組設定計算模型;model mohr range group block not;除擋墻外的群組為摩爾庫侖模型model elas range group block;擋墻為彈性assign mechanical properties設置力學參數prop shear=1e8 bulk=2e8 cohes=1e10 range group soil embankmentprop shear=5e8 bulk=1e9 cohes=1e10 range group dense soil foundationprop

5、 shear=9.15e9 bulk=10e9 range group blockini dens=2100 range group block;初始化質量密度ini dens=1800 range group blockinterface 1 prop coh=0 fric 60. dil 0. kn=1e8 ks=1e8 ten 0.;接觸面參數.interface 2 prop coh=0 fric 60. dil 0. kn=1e8 ks=1e8 ten 0.model fl_iso;各項同性流體模型,計算流體必須的. mechanical boundary and initial c

6、onditions物理邊界和初始條件fix z range z=-.1 .1;固定z=0的面% Q- S; Q8 w. O& t. Q5 X$ Z! Rfix x range x=-.1 .1fix x range x=39.9 40.1;固定x=40的面fix y range y=-.1 .18 fix y range y=9.9 10.1初始應力,垂直應力為水平的2倍,在z方向有梯度變化/ R4 I) l: F* I3 l P: Oini szz -3.6e5 grad 0 0 1.8e4ini sxx -1.8e5 grad 0 0 0.9e4ini syy -1.8e5 grad 0

7、0 0.9e4 set grav 0 0 -10;設置重力加速度;記錄監測數據hist unbal;不平衡力hist gp zdisp 16,5,20;點(16,5,20)的z方向位移hist gp zdisp 30,5,20求解solve. save block1.sav;保存assign realistic strength properties設置現實的強度參數prop cohes=0 tens 0 fric 35 range group soil embankmentprop cohes=0 tens 0 fric 40 range group dense soil foundatio

8、nsave block1.sav流體分析rest block1.sav;調用保存的文件; specify fluid properties設定流體參數prop perm=1e-8 poros=0.3ini fdens=1000;初始化流體密度ini fmod=2e3;流體的體積模量set fluid pcut on;設置流體進程,assign water table設置水位water density 1000.;水的密度.water table face 0 0 20 0 10 20 40 10 20 40 0 20 ;水平面由四個節點創建, ; block wall is impermeab

9、le擋墻不透水fix pp 0 range group block;擋墻區域內孔隙水壓力為0;施加水壓力apply nstress -1e5 range x=0,16 y=0,10 z=9.9,10.1;加在土壩上的apply nstress -2e5 grad 0,0,1e4 range x=15.9,16.1 y=0,10 z=10,20;擋墻上的z向梯度變化,;啟動流體算法set fluid on8 solve1 save block2.sav動態分析rest block2.savset dyn on;啟動動態算法set large;大變形set dyn multi on;設置動態多步計

10、算,有什么區別呢, ;turns multi-stepping on or off. Multi-stepping speeds up calculations in dynamic models which have a largezone size or modulus contrast. Sub-stepping only works when dynamic mode is in operation (SET dyn on), and is effective only when the grid is nonuniform or there is a contrast in mate

11、rial properties. ;初始化速度,位移和狀態ini xvel 0 yvel 0 zvel 0ini xdisp 0 ydisp 0 zdisp 0:ini state 0 ;設置土中水的模量 set fluid modulus for water in soilsini fmod 250e6 range group block not 施加動態邊界條件; apply dynamic boundary conditionsdef wave7 wave = ampl *c2 c sin (2.0*pi*freq*dytime):;釋放邊界z=0+ |& Y+ W$ s7 m* E6

12、r! g6 Efree x y range z -.1 .1施加動態條件apply xvel 1.0 hist wave ran z -.1 .1apply yvel 1.0 hist wave ran z -.1 .11 e9apply ffset freq 10.0 ampl 0.5 ;設置頻率和放大系數dynamic histories 動態歷史記錄2 U)set dyn time 0hist resethist dytimehist gp xvel 20,10,03 hist zone pp 19,5,5 ;記錄點孔隙壓力hist zone pp 30,5,5hist zone pp 30,5,15hist zone pp 19,5,9hist gp xdis 16,5,20set dyn damp local