FLAC3D命令流Word文档格式.docx

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