FLAC3D命令流Word文档格式.docx
《FLAC3D命令流Word文档格式.docx》由会员分享,可在线阅读,更多相关《FLAC3D命令流Word文档格式.docx(57页珍藏版)》请在冰豆网上搜索。
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
3-1
;
-------------------------------------工程信息
ProjectRecordTreeexport
Title:
Simpletest
---------------------------------计算第一步
...STATE:
STATE1....
config
grid10,10
modelelastic
group'
User:
Soil'
notnull
modelelasticnotnullgroup'
propdensity=1500.0bulk=3E6shear=1E6notnullgroup'
fixxyj1
fixxi1
fixxi11
setgravity=9.81
history999unbalanced
savestate1.sav
----------------------------------计算第二步
STATE2....
initialxdisp0ydisp0
initialxvel0yvel0
modelnulli47j810
null'
i47j810
groupdelete'
history1xdispi=4,j=11
savestate2.sav
--------------------------------绘图命令
***plotcommands****
plotname:
syy
plotholdgridsyyfill
Unbalancedforce
plotholdhistory999
grid
plotholdgridmagnify20.0lredgriddisplacement
Xdis-A
plotholdhistory1line
5-1
===============================
定义球体半径和半径方向上单元网格数
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
make_sphere
利用镜像生成完整球体网格
genzoneref
genzonerefdip90
genzonerefdip90dd90
显示球体网格
plotsurf
plsetbackwh
plblgr
5-3
n
genzonbrisize112
groupsoilranz11
grouprockranz01
expgrid1.flac3d
6-1
modelmohr
propbu3e6sh1e6coh10e3fric15
histunbal
setgrav10
solveelastic
save6-1.sav
6-2
rest6-1.sav
inixd0yd0zd0xv0yv0zv0
appnstress-100e3ranz2.93.1x12y12
save6-2.sav
6-3
histid=2gpzdis113
histid=3gpzdis112
histid=4gpxdis113
histid=5gpxdis113
histid=6zoneszz113
histid=7zoneszz1.51.52.5
histid=8zonesxz1.51.52.5
save6-3.sav
6-4
rest6-3.sav
setlogon
setlogfile6-2.log
printzonestress
printgpdis
setlogoff
6-5
plotsetrot20030
plotconszzouonmagf10
plotaddhist1
setmovieavistep1file6-5.avi
moviestart
moviefinish
7-1
propbulk3e7shear1e7
fixzranz0
fixxranx0
fixxranx1
fixyrany0
fixyrany1
ploconsz
7-2
propbulk3e7shear1e7c1e10f15tension1e10
propbulk3e7shear1e7c10e3f15ten0
7-3
genzonebricksize112
propbulk3e7shear1e7coh10e3fri15ten0
solveelas
7-4
fixyran