流体与动态计算实例分析

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