当前位置：文档之家 › flac命令流

flac命令流

  • 格式：docx
  • 大小：19.39 KB
  • 文档页数：13

flac命令流

new

gen zone brick size 1 1 2

model mohr

prop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0

fix z ran z 0

fix x ran x 0

fix x ran x 1

fix y ran y 0

fix y ran y 1

ini dens 2000 ran z 0 2

ini szz -50e3 grad 0 0 20e3 ran z 0 1

ini syy -30e3 grad 0 0 10e3 ran z 0 1

ini sxx -30e3 grad 0 0 10e3 ran z 0 1

ini pp 30e3 grad 0 0 -10e3 ran z 0 2

app nstress -10e3 ran z 2

set grav 0 0 -10

solve

；深埋工程初始应力场的S-B法生成

new

gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &

p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &

size 6 6 10

model elas

prop bulk 10e10 shear 10e10

ini dens 2500

apply sxx -1e9 grad 0 0 1.1111111e7 ran x -0.1 0.1

apply sxx -1e9 grad 0 0 7.6666666e6 ran x 59.9 60.1

apply syy -1e9 grad 0 0 8.3333333e6 ran y -0.1 0.1 apply syy -1e9 grad 0 0 8.3333333e6 ran y 59.9 60.1

apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120

set grav 0 0 -10

step 30000

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

plot cont szz

常规方法生成

new

gen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 &

p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 &

size 6 6 10

model elas

prop bulk 10e10 shear 10e10

ini dens 2500

apply sxx -1e9 grad 0 0 1.1111111e7 ran x -0.1 0.1

apply sxx -1e9 grad 0 0 7.6666666e6 ran x 59.9 60.1

apply syy -1e9 grad 0 0 8.3333333e6 ran y -0.1 0.1

apply syy -1e9 grad 0 0 8.3333333e6 ran y 59.9 60.1

apply szz -1e8 z -0.1 0.1

fix x y z ran z -0.1 0.1

set grav 0 0 -10

solve

ini xdisp 0 ydisp 0 z disp 0

ini xvel 0 yvel 0 zvel 0

plot cont szz

7.3 路基施工过程模拟

gen zone brick p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group clay

gen zone brick p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 &

ratio 1.1 1 1 size 12 1 10 group clay

gen zone brick p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0

ratio 1 1 0.8 &

size 8 1 4 group soil

gen zone brick p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3

27.5 0 0 &

ratio 1 1 0.8 size 12 1 4 group soil

gen zone brick p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5

5 0 p5 0 5 5 p6 20 0 5 &

p7 20 5 5 size 8 1 5 group dam

；边界条件

fix x y z ran z -49.9 -50.1

fix x ran x -.1 .1

fix x ran x -99.9 -100.1

fix y

model mohr ran z -50 0

model null ran z 0 5

prop bulk 7.8e6 shear 3.0e6 coh 10e10 tension 1e10 ran

group soil

ini dens 1500 ran group soil

prop bulk 3.91e6 shear 1.5e6 coh 10e10 tension 1e10 ran

group clay

ini dens 1800 ran group clay

set grav 0 0 -9.8

hist id=1 unbal

solve

prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil

prop bulk 3.91e6 shear 1.5e6 coh 20e3 fric 20 ran group clay

solve save

plot con szz outline on

plot con sxx outline on

;第一次填筑

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

hist id=2 gp zdis 0 0 0

hist id=3 gp zdis 27.5 0 0

hist id=4 gp xdis 27.5 0 0

model elastic ran z 0 1

prop bulk 7.8e6 shear 3.0e6 ran z 0 1

ini dens 1500 ran z 0 1

solve

save fill-1.sav

plot con zdis outline on

plot con xdis outline

on

model elastic ran z 1 2

prop bulk 7.8e6 shear 3.0e6 ran z 1 2

ini dens 1500 ran z 1 2

solve

save fill-2.sav

model elastic ran z 2 3

prop bulk 7.8e6 shear 3.0e6 ran z 2 3

ini dens 1500 ran z 2 3

solve

save fill-3.sav

model elastic ran z 3 4

prop bulk 7.8e6 shear 3.0e6 ran z 3 4

ini dens 1500 ran z 3 4 solve

save fill-4.sav

model elastic ran z 4 5

prop bulk 7.8e6 shear 3.0e6 ran z 4 5

ini dens 1500 ran z 4 5

solve

save fill-5.sav

plot con zdis outline on

plot con xdis outline on

set log on

set logfile 1.log

restore fill-1.sav

print gp dis range id 517 any id 533 any

restore fill-2.sav

print gp dis range id 517 any id 533 any

restore fill-3.sav

print gp dis range id 517 any id 533 any

restore fill-4.sav

print gp dis range id 517 any id 533 any

restore fill-5.sav

print gp dis range id 517 any id 533 any

set log off

让土体的模量随小主应力变化的fish程序

例子

gen zon bri size 3 3 3

model elas

prop bulk 3e7 shear 1e7

ini dens 2000 fix x y z ran z -0.1 0.1

fix x ran x -0.1 0.1

fix x ran x 2.9 3.1

fix y ran y -0.1 0.1

fix y ran y 2.9 3.1

set grav 10

solve

ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0

save 8-4.sav

rest 8-4.sav

def E_modify

p_z=zone_head

d_k=704

d_n=0.38

d_pa=101325.0

loop while p_z# null

sigma_3 = -1.0*z_sig1(p_z)

E_new = d_k*d_pa*(sigma_3/d_pa)^d_n

z_prop(p_z,'young')=E_new

p_z=z_next(p_z)

endloop

end

E_modify

plot cont prop young

zone_head 与 gp_head:分别表示单元和节点的头指针

z-next（）与gp_next():分别表示下一个单元（）节点

单元遍历的程序框架为：

p_z=zone head

loop while p_z # null语句

P_z=z_next(p_z)

合集下载

FLAC3D 实例命令流1

FLAC3D 实例命令流1

文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.

1文档来源为:从网络收集整理.word版本可编辑. 第1部分 命令流按照顺序进行

2-1定义一个FISH函数

new

def abc

abc = 25 * 3 + 5

End

print abc

2-2使用一个变量

new

def abc

hh = 25

abc = hh * 3 + 5

End

Print hh

Print abc

2-3对变量和函数的理解

new

def abc

hh = 25

abc = hh * 3 + 5

End

set abc=0 hh=0

print hh

print abc

print hh

new

def abc

abc = hh * 3 + 5

end

set hh=25

print abc

set abc=0 hh=0

print hh

print abc

print hh

2-4获取变量的历史记录

new

gen zone brick size 1 2 1

model mohr

prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10

fix x y z range y -0.1 0.1

apply yvel -1e-5 range y 1.9 2.1

plot set rotation 0 0 45 文档来源为:从网络收集整理.word版本可编辑.欢迎下载支持.

2文档来源为:从网络收集整理.word版本可编辑. plot block group

def get_ad

ad1 = gp_near(0,2,0)

ad2 = gp_near(1,2,0)

ad3 = gp_near(0,2,1)

ad4 = gp_near(1,2,1)

end

get_ad

def load

load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4)

FLAC渗流命令流共11页

FLAC渗流命令流共11页

第 1 页 ini pp 0 j=0

fix pp j=0

这样设置后是透水边界了

底边界不透水,你可以不用管它了,默认为不透水边界

我现在在三维中模拟基坑开挖，x(0~42) y(0~30) z(-40~0),初始水位在－3米处，我的初始水位命令这样写对不对啊？

ini pp -3e4 grad 0 0 -1e4 range x 0 42 y 0 30 z -40 -3

range -3.1 -2.9

;;

在开挖中需要降水，降至－7米，需开挖的区域为 x(0~8) y(0~30)

fix pp range x 8 42 y o 30 z -3.1 -2.9

fix pp 0 range x 0 8 y 0 30 z -7.1 -6.9 pl contour pp

可以显示孔隙水压力场，一般就认为就是渗流场吧！

set flow on 即选择是否要进行渗流分析，其前提是要set fluid on方可以！off就是相反了。但是，set fluid on以后不一定进行set flow on的计算，这点我的理解（有待证实）是，仅有set fluid

on的话，就是把里面的pp当作是静水压力来处理的，而当set flow on 打开以后，就可以进行真正意义上的渗流分析了！再一个，set flow on的一个前提必须设定perm值（渗透系数），否则程序无法继续进行！而set fluid on就没有那么严格，只要是config fluid 以及设定模型以及fmod以后即可以开始计算了。

set mech on就表示把力学计算这部分打开，在默认的情况下，是打开的。在流固耦合有一节中讲到，当只要使系统产生水压力的分布，而不需要由于机械变形（mechanical deformation）而产生的pp的时候，就必须要set mech off，也就是说这个时候只有fluid的计算，而没有mechanical的计算了。

FLAC若干问题的解答

FLAC若干问题的解答

FLAC若⼲问题的解答

把结果中的信息输出到指定⽂件

分为两种：1、把单元信息输出到指定⽂件

set log on

set log on finame.dat;可以输⼊⽂件路径，否则按当前⽬录处理

print zone

set log off通过反复使⽤该命令，可以把不同信息输出到不同的⽂件。

2、把节点信息输出到指定⽂件

set log on

set log on finame.dat;可以输⼊⽂件路径，否则按当前⽬录处理

print gp

set log off

请教⼀下print gp

⾥的keyword该如何填。⽐如想输出某点的zdis

谢谢print gp disp

print gp position

log ⽂件中有，然后处理，这些论坛⾥都有的

关于flac3d—内置fish语⾔精讲FISH语⾔是FLAC3D程序的内置编程语⾔

因为FLAC3D的最佳操作⽅式是命令流⽂件⽅式

这⼀点与ANSYS很相似，⽽FISH就相当于ANSYS的APDL语⾔。它包括循环、判断等结构。如果你还⽤过其它⾼级语⾔，那么从形式上讲你也可以把它理解为⼦函数。FISH语⾔的引⼊极⼤的⽅便了⽤户进⾏复杂的程序建模

它不但可以嵌⼊命令流⽂件⾥⼯作

⽽且还可以引⽤FLAC3D本⾝的任何命令

所以说它实现了对FLAC3D的完全控制

引⼊FISH语⾔的作法值得其它通⽤软件效仿。

以上仅是个⼈的使⽤体会请版主与管理员指教

也供初学者⼊门之⽤。

谢谢！Fish函数增加了以下⼏种新的特性：

1.增加了fish变量来获取结点、单元和界⾯变量2.fish提供了获取结构单元变量的途径

3.休单元和⾯单元性质⽬前可以通过单元变量名z_prop(i_z string)和界⾯单元名i_prop(i_z string)分别加以识别

4.fish函数可以获取单元应变和应⼒速率，还提供了全应变增量张量和应变速率张量

5.提供了fish绘图⼦程序函数够⽣⽤户定义的图形内容

6.fish函数已经增加了从⽂件读、写数据的能⼒

matlab读flac命令流

matlab读flac命令流

ssxx1=importdata('1sxx.txt');

ssxx2=importdata('2sxx.txt');

ssxx3=importdata('3sxx.txt'); ssxx4=importdata('4sxx.txt');

ssxx5=importdata('5sxx.txt');

ssxx6=importdata('6sxx.txt');

ssxx7=importdata('7sxx.txt');

ssxx8=importdata('8sxx.txt'); ssxx9=importdata('9sxx.txt');

ssxx10=importdata('10sxx.txt');

ssxx11=importdata('11sxx.txt');

ssxx12=importdata('12sxx.txt');

ssxx13=importdata('13sxx.txt'); ssxx14=importdata('14sxx.txt');

ssxx15=importdata('15sxx.txt');

ssxx16=importdata('16sxx.txt');

ssxx17=importdata('17sxx.txt');

ssxx18=importdata('18sxx.txt'); ssxx19=importdata('19sxx.txt');

ssxx20=importdata('20sxx.txt');

ssxx21=importdata('21sxx.txt');

ssxx22=importdata('22sxx.txt');

ssxx23=importdata('23sxx.txt'); ssxx24=importdata('24sxx.txt');

ssxx25=importdata('25sxx.txt');

ssxx26=importdata('26sxx.txt');