FLAC3D命令流.docx

FLAC3D命令流.docx

《FLAC3D命令流.docx》由会员分享，可在线阅读，更多相关《FLAC3D命令流.docx（57页珍藏版）》请在冰豆网上搜索。

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

fixyran