FLAC3D命令流文档格式.docx

FLAC3D命令流文档格式.docx

《FLAC3D命令流文档格式.docx》由会员分享，可在线阅读，更多相关《FLAC3D命令流文档格式.docx（66页珍藏版）》请在冰点文库上搜索。

;

-------------------------------------工程信息

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

iniszz-40e3grad0020e3ranz02

inisyy-20e3grad0010e3ranz02

inisxx-20e3grad0010e3ranz02

7-5

modelm

propbulk3e7shear1e7c10e10f15ten1e10

inidens2000ranz01

inidens1500ranz12

iniszz-35e3grad0020e3ranz01

inisyy-17.5e3grad0010e3ranz01

inisxx-17.5e3grad0010e3ranz01

iniszz-15e3grad0015e3ranz12

inisyy-7.5e3grad007.5e3ranz12

inisxx-7.5e3grad007.5e3ranz12

inipp10e3grad00-10e3ranz01

7-6

inidens2000ranz02

iniszz-50e3grad0020e3ranz01

inisyy-30e3grad0010e3ranz01

inisxx-30e3grad0010e3ranz01

inipp30e3grad00-10e3ranz02

appnstress-10e3ranz2

7-7

genzonebrickp0000p16000p20600p30090&

p460600p506090p6600150p76060150&

size6610

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

step30000

inixdisp0ydisp0zdisp0

inixvel0yvel0zvel0

plocontszz

7-8

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

7-9

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

fixxranx99.9100.1

fixy

modelmohrranz-500

modelnullranz05

propbulk7.8e6shear3.0e6coh10e10tension1e10rangroupsoil

inidens1500rangroupsoil

propbulk3.91e6shear1.5e6coh10e10tension1e10rangroupclay

inidens1800rangroupclay

setgrav00-9.8

histid=1unbal

propbulk7.8e6shear3.0e6coh10e3fric15rangroupsoil

propbulk3.91e6shear1.5e6coh20e3fric20rangroupclay

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

savefill-2.sav

modelelasticranz23

propbulk7.8e6shear3.0e6ranz23

inidens1500ranz23

savefill-3.sav

modelelasticranz34

propbulk7.8e6shear3.0e6ranz34

inidens1500ranz34

savefill-4.sav

modelelasticranz45

propbulk7.8e6shear3.0e6ranz45

inidens1500ranz45

savefill-5.sav

pau

ploblgr

genzonebrickp0000p110000p2050p3005size

genzonebrickp000-50p127.50-50p205-50p300-10size8110groupclay

p5055p62005p72055size815groupdam

7-10

setlogon;

打开log记录

setlogfile1.log;

设置记录文件名为：

1.log

restorefill-1.sav;

调用保存的文件

printgpdisrangeid517anyid533any;

输出两个节点的变形值

restorefill-2.sav

printgpdisrangeid517anyid533any

restorefill-3.sav

restorefill-4.sav

restorefill-5.sav

setlogoff;