FLAC3D命令流.docx
《FLAC3D命令流.docx》由会员分享，可在线阅读，更多相关《FLAC3D命令流.docx（66页珍藏版）》请在冰豆网上搜索。
FLAC3D命令流
命令流
1
genzonebrip0000p11000p20100p30010&
p415150p501515p615010p7202020&
size101010rat1.00.91.1groupbrick_1
genzonebrip02000p1add1000p2add0200p3add0015&
size101010rat1.00.91.1groupbrick_2
genzonebrip04000edge10size101010rat1.00.91.1groupbrick_3
plotsur
2-1
new
genzonbrisize333
modelelas
propbulk3e8shear1e8
inidens2000
fixzranz-.1.1
fixxranx-.1.1
fixxranx2.93.1
fixyrany-.1.1
fixyrany2.93.1
setgrav00-10
solve
appnstress-10e4ranz3x12y12
histgpvel003
histgpvel033
plohist1red
ploaddhist2blue
solve
3-1
;-------------------------------------工程信息
;ProjectRecordTreeexport
;Title:
Simpletest
;---------------------------------计算第一步
;...STATE:
STATE1....
config
grid10,10
modelelastic
group'User:
Soil'notnull
modelelasticnotnullgroup'User:
Soil'
propdensity=1500.0bulk=3E6shear=1E6notnullgroup'User:
Soil'
fixxyj1
fixxi1
fixxi11
setgravity=9.81
history999unbalanced
solve
savestate1.sav
;----------------------------------计算第二步
;...STATE:
STATE2....
initialxdisp0ydisp0
initialxvel0yvel0
modelnulli47j810
group'null'i47j810
groupdelete'null'
history1xdispi=4,j=11
solve
savestate2.sav
;--------------------------------绘图命令
;***plotcommands****
;plotname:
syy
plotholdgridsyyfill
;plotname:
Unbalancedforce
plotholdhistory999
;plotname:
grid
plotholdgridmagnify20.0lredgriddisplacement
;plotname:
Xdis-A
plotholdhistory1line
5-1
new
;===============================
;定义球体半径和半径方向上单元网格数
;===============================
defparm
rad=10.0
rad_size=5
end
parm
;===============================
;建立八分之一球体外接立方体网格
;===============================
genzonepyramidp0rad00p1rad0radp2radrad0p3000&
p4radradradsizerad_sizerad_sizerad_sizegroup1
genzonepyramidp00rad0p1radrad0p20radradp3000&
p4radradradsizerad_sizerad_sizerad_sizegroup2
genzonepyramidp000radp10radradp2rad0radp3000&
p4radradradsizerad_sizerad_sizerad_sizegroup3
;==================================
;利用FISH语言将内部立方体节点调整到球面
;==================================
defmake_sphere
p_gp=gp_head
loopwhilep_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）
ifdist>0then
;节点位置调整
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*pz
end_if
p_gp=gp_next（p_gp）
end_loop
end
make_sphere
;===============================
;利用镜像生成完整球体网格
;===============================
genzoneref
genzonerefdip90
genzonerefdip90dd90
;===============================
;显示球体网格
;===============================
plotsurf
plsetbackwh
plblgr
5-3
n
genzonbrisize112
groupsoilranz11
grouprockranz01
expgrid1.flac3d
6-1
genzonbrisize333
modelmohr
propbu3e6sh1e6coh10e3fric15
fixzranz-.1.1
fixxranx-.1.1
fixxranx2.93.1
fixyrany-.1.1
fixyrany2.93.1
inidens2000
histunbal
setgrav10
solveelastic
save6-1.sav
6-2
rest6-1.sav
inixd0yd0zd0xv0yv0zv0
appnstress-100e3ranz2.93.1x12y12
solve
save6-2.sav
6-3
rest6-1.sav
inixd0yd0zd0xv0yv0zv0
appnstress-100e3ranz2.93.1x12y12
histid=2gpzdis113
histid=3gpzdis112
histid=4gpxdis113
histid=5gpxdis113
histid=6zoneszz113
histid=7zoneszz1.51.52.5
histid=8zonesxz1.51.52.5
solve
save6-3.sav
6-4
rest6-3.sav
setlogon
setlogfile6-2.log
printzonestress
printgpdis
setlogoff
6-5
rest6-1.sav
inixd0yd0zd0xv0yv0zv0
appnstress-100e3ranz2.93.1x12y12
plotsetrot20030
plotconszzouonmagf10
plotaddhist1
setmovieavistep1file6-5.avi
moviestart
solve
moviefinish
7-1
n
genzonbrisize112
modelelas
propbulk3e7shear1e7
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000
setgrav00-10
solve
ploconsz
7-2
n
genzonbrisize112
modelmohr
propbulk3e7shear1e7c1e10f15tension1e10
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000
setgrav00-10
solve
propbulk3e7shear1e7c10e3f15ten0
solve
ploconsz
7-3
n
genzonebricksize112
modelmohr
propbulk3e7shear1e7coh10e3fri15ten0
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000
setgrav00-10
solveelas
ploconsz
7-4
new
genzonebricksize112
modelmohr
propbulk3e7shear1e7coh10e3fri15ten0
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000
iniszz-40e3grad0020e3ranz02
inisyy-20e3grad0010e3ranz02
inisxx-20e3grad0010e3ranz02
setgrav00-10
solve
ploconsz
7-5
n
genzonbrisize112
modelm
propbulk3e7shear1e7c10e10f15ten1e10
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000ranz01
inidens1500ranz12
iniszz-35e3grad0020e3ranz01
inisyy-17.5e3grad0010e3ranz01
inisxx-17.5e3grad0010e3ranz01
iniszz-15e3grad0015e3ranz12
inisyy-7.5e3grad007.5e3ranz12
inisxx-7.5e3grad007.5e3ranz12
inipp10e3grad00-10e3ranz01
setgrav00-10
solve
ploconsz
7-6
n
genzonbrisize112
modelm
propbulk3e7shear1e7c10e10f15ten1e10
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
inidens2000ranz02
iniszz-50e3grad0020e3ranz01
inisyy-30e3grad0010e3ranz01
inisxx-30e3grad0010e3ranz01
inipp30e3grad00-10e3ranz02
appnstress-10e3ranz2
setgrav00-10
solve
ploconsz
7-7
new
genzonebrickp0000p16000p20600p30090&
p460600p506090p6600150p76060150&
size6610
modelelas
probulk10e10she10e10
iniden2500
applysxx-1e9grad001.1111111e7rangex-.1.1
applysxx-1e9grad006.6666666e6rangex59.960.1
applysyy-1e9grad008.3333333e6rangey-.1.1
applysyy-1e9grad008.3333333e6rangey59.960.1
applyszz-1e8grad008.3333333e5ranz0120
setgrav00-10
step30000
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
plocontszz
7-8
new
genzonebrickp0000p16000p20600p30090&
p460600p506090p6600150p76060150&
size6610
modelelas
probulk10e10she10e10
iniden2500
inisxx-1e9grad001.1111111e7rangex-.1.1
inisxx-1e9grad006.6666666e6rangex59.960.1
inisyy-1e9grad008.3333333e6rangey-.1.1
inisyy-1e9grad008.3333333e6rangey59.960.1
iniszz-1e8ranz-.1.1
fixxyzranz-.1.1
setgrav00-10
solve
inixdisp0ydisp0zdisp0
inixvel0yvel0zvel0
plocontszz
7-9
new
genzonebrickp000-50p127.50-50p205-50p300-10size8110groupclay
genzonebrickp027.50-50p11000-50p227.55-50p327.50-10ratio1.111size12110groupclay
genzonebrickp000-10p127.50-10p205-10p3000ratio110.8size814groupsoil
genzonebrickp027.50-10p11000-10p227.55-10p327.500ratio1.110.8size1214groupsoil
genzonebrickp0000p127.500p2050p3005p427.550&
p5055p62005p72055size815groupdam
fixxyzranz-49.9-50.1
fixxranx-.1.1
fixxranx99.9100.1
fixy
modelmohrranz-500
modelnullranz05
propbulk7.8e6shear3.0e6coh10e10tension1e10rangroupsoil
inidens1500rangroupsoil
propbulk3.91e6shear1.5e6coh10e10tension1e10rangroupclay
inidens1800rangroupclay
setgrav00-9.8
histid=1unbal
solve
propbulk7.8e6shear3.0e6coh10e3fric15rangroupsoil
propbulk3.91e6shear1.5e6coh20e3fric20rangroupclay
solve
saveelastic.sav
inixdis0ydis0zdis0;将节点位移清零
inixvel0yvel0zvel0;将节点速度清零
histid=2gpzdis000;记录地基顶部中心点的沉降
histid=3gpzdis27.500;记录路基坡脚处的沉降
histid=4gpxdis27.500;记录路基坡脚处的水平位移
modelelasticranz01;;激活0m~1m的单元
propbulk7.8e6shear3.0e6ranz01
inidens1500ranz01
solve;按软件默认精度求解
savefill-1.sav
modelelasticranz12
propbulk7.8e6shear3.0e6ranz12
inidens1500ranz12
solve
savefill-2.sav
modelelasticranz23
propbulk7.8e6shear3.0e6ranz23
inidens1500ranz23
solve
savefill-3.sav
modelelasticranz34
propbulk7.8e6shear3.0e6ranz34
inidens1500ranz34
solve
savefill-4.sav
modelelasticranz45
propbulk7.8e6shear3.0e6ranz45
inidens1500ranz45
solve
savefill-5.sav
pau
;ploblgr
;genzonebrickp0000p110000p2050p3005size
genzonebrickp000-50p127.50-50p205-50p300-10size8110groupclay
genzonebrickp027.50-50p11000-50p227.55-50p327.50-10ratio1.111size12110groupclay
genzonebrickp000-10p127.50-10p205-10p3000ratio110.8size814groupsoil
genzonebrickp027.50-10p11000-10p227.55-10p327.500ratio1.110.8size1214groupsoil
genzonebrickp0000p127.500p2050p3005p427.550&
p5055p62005p72055size815groupdam
7-10
setlogon;打开log记录
setlogfile1.log;设置记录文件名为：
1.log
restorefill-1.sav;调用保存的文件
printgpdisrangeid517anyid533any;输出两个节点的变形值
restorefill-2.sav
printgpdisrangeid517anyid533any
restorefill-3.sav
printgpdisrangeid517anyid533any
restorefill-4.sav
printgpdisrangeid517anyid533any
restorefill-5.sav
printgpdisrangeid517anyid533any
setlogoff;
升级会员 