流体与动态计算实例分析.doc

流体与动态计算实例分析---堤坝 ;流体与动态计算实例分析---大坝 new conf 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,5;堤坝;在z=10,x=15~22,y=0~10上生成接触面1 interface 1 face range x 15.0 22.0 y 0.0 10.0 z 9.9 10.1 ;在x=22,z=10~21,y=0~10上生成接触面2 interface 2 face range x 21.9 22.1 y 0.0 10.0 z 10.0 21.0 \设置最大的边界长度为1.0\ interface 1 maxedge 1.0 interface 2 maxedge 1.0 generate block wall 生成挡水墙 gen zone bri p0 15,0.5,11 p1 21,0.5,11 p2 15,9.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,20 group 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 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 embankment prop shear=5e8 bulk=1e9 cohes=1e10 range group dense soil foundation prop shear=9.15e9 bulk=10e9 range group block ini dens=2100 range group block;初始化质量密度 ini dens=1800 range group block interface 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 conditions物理边界和初始条件 fix z range z=-.1 .1;固定z=0的面% Q- S; Q8 w. O t. Q5 X$ Z! R fix x range x=-.1 .1 fix 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: O ini szz -3.6e5 grad 0 0 1.8e4 ini sxx -1.8e5 grad 0 0 0.9e4 ini syy -1.8e5 grad 0 0 0.9e4 set grav 0 0 -10;设置重力加速度 ;记录监测数据 hist unbal;不平衡力 hist gp zdisp 16,5,20;点(16,5,20)的z方向位移 hist gp zdi