《flac3d命令流资料》由会员分享，可在线阅读，更多相关《flac3d命令流资料（22页珍藏版）》请在金锄头文库上搜索。
1、命令流1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1.0 0.9 1.1 group brick_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1.0 0.9 1.1 group brick_2gen zone bri p0 40 0 0 edge 10 size 10 10
2、10 rat 1.0 0.9 1.1 group brick_3 plot sur 2-1newgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2hist gp vel 0 0 3hist gp vel 0 3 3
3、plo hist 1 redplo add hist 2 bluesolve3-1;-工程信息;Project Record Tree export;Title:Simple test;-计算第一步;. STATE: STATE1 .configgrid 10,10model elasticgroup User:Soil notnullmodel elastic notnull group User:Soilprop density=1500.0 bulk=3E6 shear=1E6 notnull group User:Soilfix x y j 1fix x i 1fix x i 11se
4、t gravity=9.81history 999 unbalancedsolvesave state1.sav;-计算第二步;. STATE: STATE2 .initial xdisp 0 ydisp 0initial xvel 0 yvel 0model null i 4 7 j 8 10group null i 4 7 j 8 10group delete nullhistory 1 xdisp i=4, j=11solvesave state2.sav;-绘图命令;* plot commands *;plot name: syyplot hold grid syy fill;plot
5、 name: Unbalanced forceplot hold history 999;plot name: gridplot hold grid magnify 20.0 lred grid displacement;plot name: Xdis-Aplot hold history 1 line5-1new; =; 定义球体半径和半径方向上单元网格数; =def parm rad=10.0 rad_size=5endparm; =; 建立八分之一球体外接立方体网格; =gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0
6、0 0 &p4 rad rad rad size rad_size rad_size rad_size group 1gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 &p4 rad rad rad size rad_size rad_size rad_size group 2gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 &p4 rad rad rad size rad_size rad_size rad_size group 3; =;
7、利用FISH语言将内部立方体节点调整到球面; =def make_sphere p_gp=gp_head loop while p_gp#null; 获取节点点坐标值：P=(px,py,pz) px=gp_xpos(p_gp) py=gp_ypos(p_gp) pz=gp_zpos(p_gp) dist=sqrt(px*px+py*py+pz*pz)if dist0 then; 节点位置调整 maxp=max(px,max(py,pz) k=(maxp/rad)*(rad/dist) gp_xpos(p_gp)=k*px gp_ypos(p_gp)=k*py gp_zpos(p_gp)=k*p
8、z end_if p_gp=gp_next(p_gp) end_loopendmake_sphere; =; 利用镜像生成完整球体网格; =gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90 ; =; 显示球体网格; =plot surfpl set back whpl bl gr5-3ngen zon bri size 1 1 2group soil ran z 1 1group rock ran z 0 1expgrid 1.flac3d6-1gen zon bri size 3 3 3model mohrprop bu 3e6
9、 sh 1e6 coh 10e3 fric 15fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbalset grav 10solve elasticsave 6-1.sav6-2rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2solvesave 6-2.sav6-3rest 6-1.savi
10、ni xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solvesave 6-3.sav6-4rest 6-3.savset log onset lo
11、gfile 6-2.logprint zone stressprint gp disset log off6-5rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie avi step 1 file 6-5.avimovie startsolvemovie finish7-1ngen zon bri size 1 1 2model e
12、lasprop bulk 3e7 shear 1e7fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveplo con sz7-2ngen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 1e10 f 15 tension 1e10fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set
13、 grav 0 0 -10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solveplo con sz7-3ngen zone brick size 1 1 2model mohrprop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve elasplo con sz7-4newgen zone brick size 1 1 2mod
