完整word版,FLAC3D 实例命令流1

第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

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) end

hist load

hist gp ydis 0,2,0

step 1000

plot his 1 vs -2

2-5用FISH函数计算体积模量和剪砌模量

new

def derive

s_mod = y_mod / (2.0 * (1.0 + p_ratio))

b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))

end

set y_mod = 5e8 p_ratio = 0.25

derive

print b_mod

print s_mod

2-6 在FLAC输入中使用符号变量

New

def derive

s_mod = y_mod / (2.0 * (1.0 + p_ratio))

b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))

end

set y_mod = 5e8 p_ratio = 0.25

derive

gen zone brick size 2,2,2

model elastic

prop bulk=b_mod shear=s_mod

print zone prop bulk

print zone prop shear

2-7 控制循环

New

def xxx

sum = 0

prod = 1

loop n (1,10)

sum = sum + n

prod = prod * n

end_loop

end

xxx

print sum, prod

new

gen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elas

plot set rotation 0 0 45

plot block group

def install

pnt = zone_head

loop while pnt #null

z_depth = -z_zcen(pnt)

y_mod = y_zero + cc * sqrt(z_depth)

z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio))

z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))

pnt = zone_next(pnt)

end_loop

end

set p_ratio=0.25 y_zero=1e7 cc=1e8

install

2-8 拆分命令行

new ;example of a sum of many things

def long_sum

temp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10 long_sum = temp + v11 + v12 + v13 + v14 + v15

end

2-9 变量类型

new

def haveone

aa = 2

bb = 3.4

cc = ’Have a nice day’

dd = aa * bb

ee = cc + ’, old chap’

end

haveone

print fish

2-10 IF条件语句

new

def abc

if xx > 0 then

abc = 1000

else

abc = -1000

end_if

end

set xx = 10

print abc

set xx = 0

print abc

2-11 索单元自动生成

new

gen zone brick size 10 3 5

plot set rotation 0 0 45

plot block group

def place_cables

loop n (1,5)

z_d = float(n) - 0.5

command

sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_command

end_loop

end

place_cables

plot grid sel geom red

new

gen zone brick size 10 3 5

plot set rotation 15 0 60

plot block group

mod mohr

prop bulk 1e8 shear .3e8 fric 35

prop coh 1e3 tens 1e3

ini dens 1000

set grav 0,0,-10

fix x y z range z -.1 .1

fix y range y -.1 .1

fix y range y 2.9 3.1

fix x range x -.1 .1

fix x range x 9.9 10.1

set large

hist unbal

solve

save cab_str.sav

ini xdis 0 ydis 0 zdis 0

hist gp xdisp 0,1,5

def place_cables

loop n (1,5)

z_d = 5.5 - float(n)

z_t = z_d + 0.5

z_b = z_d - 0.5

command

free x range x -.1,.1 z z_b z_t

solve

sel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7

sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7

sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7

sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 &

gr_k 2e10 gr_coh 1e10 gr_per 1.0

end_command

end_loop

end

place_cables

save cab_end.sav

plot sketch sel cable force red

2-12圆形隧道开挖模拟计算

;建立模型

gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &

size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderock

gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &

size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen zon reflect dip 90 dd 90 orig 0 0 0

gen zon reflect dip 0 dd 0 ori 0 0 0

gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3

gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4

gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1

gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1

gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1

;绘制模型图

plot block group

plot add axes red

;plot set rotation 0 0 45 用于显示三维模型

;设置重力

set gravity 0 0 -10

;给定边界条件

fix z range z -12.01,-11.99

fix x range x -21.01,-20.99

fix x range x 20.99,21.01

fix y range y -0.01 0.01

fix y range y 0.99,1.01

;求解自重应力场

model mohr

ini density 1800 ;围岩的密度

prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度

set mech ratio=1e-4

solve

save Gravsol.sav

plot cont zdisp outl on

plot cont szz

;毛洞开挖计算

initial xdisp=0 ydisp=0 zdisp=0

model null range group insiderock any group concretliner any

plot block group

plot add axes red

set mech ratio=5e-4

solve

save Kaiwsol.sav

plot cont zdisp

plot cont sdisp

plot cont szz

plot cont xzz

;模筑衬砌计算

model elas range group concretliner any

plot block group

plot add axes red

ini density 2500 range group concretliner any ;衬砌混凝土的密度

prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模

set mech ratio=1e-4

solve

save zhihusol.sav

plot cont zdisp

plot cont sdisp

plot cont szz

plot cont xzz

;完成计算分析

完整word版,FLAC3D中应变软化硬化模型参数取值(table)

一FLAC3D中应变软化\硬化模型参数取值（table） 例：model ss prop den 2500 bulk 2e8 shear 1e8 co 2e6 fric 45 ten 1e6 dil 10 prop ftab 1 ctab 2 dtab 3 table 1 0，45 0. 05，42 0. 1，40 table 2 0，2e6 0. 05，1e6 0. 1，5e5 table 3 0，10 0. 05，3 0. 1，0 Table 的功能是，用来定义某一变量的值随另一变量的值而发生变化的变化规律。 例如，对于应变软化模型（strain-softening model ），可使用PROP ftable, PROP ctable 和PROP dtable等命令来设置摩擦角（friction）、黏聚力（cohesion）和剪胀角（dilation），随着累积塑性剪切应变（accumulated plastic shear strain ）而变化的规律； 或者使用PROP ttable命令，来定义抗拉强度（tensile strength ）随着累积塑性拉应变（accumulated plastic tensile strain ）而变化的规律。 Table命令的调用格式为： TABLE n x1 y1 . . . 定义1号table（ftab 1；即定义摩擦角的变化规律）， 累积塑性切应变取值为0 时，摩擦角取值为45； 累积塑性切应变取值为0.05 时，摩擦角取值为42； 累积塑性切应变取值为0.1 时，摩擦角取值为40；

FLAC3D快速入门及简单实例

FLAC3D快速入门 及简单实例 李佳宇编 LJY指南针教程

前言 FLAC及FLAC3D是由国际著名学者、英国皇家工程院院士、离散元的 发明人Peter Cundall博士在70年代中期开始研究的，主要面对岩土工程的通 用软件系统，目前已经在全球70多个国家得到广泛应用，在岩土工程学术界 和工业界赢得了广泛的赞誉。前国际岩石力学会主席 C.Fairhurst(1994)对 FLAC程序的评价是：“现在它是国际上广泛应用的可靠程序。” 我从研二(2010年)开始接触FLAC3D，最初的原因是导师要求每一个人至 少学会一个数值计算软件，而他嘴里每天念叨最多的就是FLAC，自己当时对数 值计算一无所知，便答应老师要学会FLAC3D。第一次打开软件界面，我心里 就凉了大半截，面对着一个操作界面跟记事本无异的所谓“功能强大”的岩土工 程专业软件，半点兴趣也提不起来。年底，从项目工地回到学校准备论文开题， 老师对我的开题报告非常不满意，当着全教研室师生的面，劈头盖脸大批一顿， 第二天又找谈话。在巨大的压力和强烈的自尊心驱使下，我硬着头皮开始啃 FLAC3D，一个半月之后，终于有了初步的计算结果，对老师有个交代，我也能 回家过年了。 前面这一段过程可能是大多数FLAC3D初学者的必经阶段，或者是即将 开始软件学习的人惧怕的事情。毫无疑问，FLAC3D极其不友好的界面是阻碍 初学者前进的很大障碍，当然还包括它是一个全英文的软件。但是当你费尽周折 的走进FLAC3D的世界，你就会发现它独特的魅力，比如简洁的界面，快捷的 命令流操作，高效的计算方法，不易报错等等。另外一个拿不上台面的优点就是 它非常小巧，包括Manual在内一共才几十兆大小，而且已经被破解成绿色版， 只要把它和命令流装进U盘，你就可以随便找一个身边功能最强大的电脑开始 计算了，如果你有过ANSYS、ABAQUS等大型软件痛苦的安装经历，你便能 毕业之后，本以为不用再接触数值计算，但工作需要使得我又一次开始与理解“绿色版”的含义，当然还请大家尊重知识产权，支持正版。 FLAC3D进行亲密接触。我的领导给了我很多新思路和很大的支持，如今我的 水平比研究生时有了不少提高。于是，我想把我的经验总结成文，希望对初学 者起

完整word版,FLAC3D 实例命令流1

第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 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) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear

Flac3D常用命令流Word版

FLAC3D 常用命令流 ;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl contour szz outline on ;在模型中显示位移-应变曲线 hist gp ydisp 0,0,0 hist zone syy 0,1,0 hist zone syy 1,1,0 pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror" ;对模型进行压缩实验的方法 ;即在模型两侧施加相反方向的速度 ini yvel 1e-7 range y -.1 .1 ini yvel -1e-7 range y 1.9 2.1 ;修改模型的坐标值 ini x add -100 y add -100 z add -100 ;显示云图的同时也显示模型网格轮廓 plot add cont disp outline on ;gradient更精确 ;输入角度、弧度方法pi=π，90°为90.0*degrad def set_vals ptA = 25.0 * sin(pi/2);ptA=25.0 ptB = 25.0 * cos( 60.0*degrad );ptB=12.5 ptC=pi;ptC=3.1415926 end set_vals print ptA ptB ptC ;施加结构单元方法 sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) & end2=(0.0,25.0,0.0) radius=24.5 not plot add sel geom black black cid on scale=0.03 sel node init zpos add -25.0 ;如何显示某一平面 plot create name_plane

FLAC模拟隧道开挖支护的实例

;FLAC3D3.0在某隧道工程开挖支护中的应用 ;隧道建模命令流入下： new set log on set logfile yang.log gen zon radcyl p0 0 0 0 p1 9.0 0 0 p2 0 50 0 p3 0 0 8 & size 4 20 6 4 dim 6 5 6 5 rat 1 1 1 1 group 围岩 gen zon cshell p0 0 0 0 p1 6.0 0 0 p2 0 50 0 p3 0 0 5.0 & size 4 20 6 4 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group 初期支护 gen zon cshell p0 0 0 0 p1 5.6 0 0 p2 0 50 0 p3 0 0 4.6 & size 4 20 6 4 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬砌 fill group 原岩gen zon radcyl p0 0 0 0 p1 0 0 -8.0 p2 0 50 0 p3 9.0 0 0 & size 4 20 6 4 dim 3 6 3 6 rat 1 1 1 1 group 围岩2 gen zon cshell p0 0 0 0 p1 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 & size 4 20 6 4 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group 仰拱初期支护 gen zon cshell p0 0 0 0 p1 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 & size 4 20 6 4 dim 2 5 2 5 rat 1 1 1 1 group 仰拱二次衬砌 fill group 仰拱原岩gen zone reflect normal -1 0 0 gen zone radtun p0 0 0 0 p1 45 0 0 p2 0 50 0 p3 0 0 20 & size 3 20 3 12 dim 9 8 9 8 rat 1 1 1 1.1 group 围岩3 gen zon reflect dip 0 ori 0 0 0 range x 0 9 y 0 50 z 8 20 gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z 8 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20 gen zon brick p0 -45 0 -20 p1 -45 0 -40 p2 -45 50 -20 p3 45 0 -20 & size 5 20 6 rat 1.1 1 1 group 围岩4 save tun_model.sav ;假设围岩岩体符合mohr-coulomb本构模型，给围岩赋参数命令流如下， ; mohr-coulomb model model mohr def derive s_mod1=E_mod1/(2.0*(1.0+p_ratio1)) b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) s_mod2=E_mod2/(2.0*(1.0+p_ratio2)) b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2)) end set E_mod1=0.6e9 p_ratio1=0.27 E_mod2=0.8e9 p_ratio2=0.26 derive prop bulk b_mod1 shear s_mod1 cohe 1.8e6 tens 0.8e6 fric 30 range z 4.5 20 prop bulk b_mod2 shear s_mod2 cohe 2.8e6 tens 1.0e6 fric 35 range z -40 4.5 ini dens=2300 set grav 0 0 -10 ; boundary and initial conditions

FLAC3D实例命令流1

FLAC3D实例命令流1 第1部分命令流按照顺序进行 2-1定义一个FISH函数newdefabc abc=25某3+5Endprintabc 2-2使用一个变量newdefabchh=25 abc=hh某3+5EndPrinthhPrintabc 2-3对变量和函数的理解newdefabchh=25 abc=hh某3+5End etabc=0hh=0printhhprintabcprinthhnewdefabc abc=hh某3+5end ethh=25printabc etabc=0hh=0printhhprintabcprinthh 2-4获取变量的历史记录new genzonebrickize121modelmohr prophear=1e8bulk=2e8cohe=1e5ten=1e10 fi某某yzrangey-0.10.1 applyyvel-1e- 5rangey1.92.1plotetrotation0045plotblockgroupdefget_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)endget_addefload

load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfun bal(ad4)endhitload hitgpydi0,2,0tep1000 plothi1v-2 2-5用FISH函数计算体积模量和剪砌模量new defderive _mod=y_mod/(2.0某(1.0+p_ratio)) b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end ety_mod=5e8p_ratio=0.25derive printb_modprint_mod 2-6在FLAC输入中使用符号变量New defderive _mod=y_mod/(2.0某(1.0+p_ratio)) b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end ety_mod=5e8p_ratio=0.25derive genzonebrickize2,2,2modelelatic propbulk=b_modhear=_modprintzonepropbulkprintzoneprophear 2-7控制循环Newdef某某某 um=0prod=1

FLAC3D命令流(整理版)

精心整理1、怎样查看模型？ 答：plotgrid可以查看网格，plotgridnum可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2

flac命令流参考实用

flac命令流参考实用 ;**************************************** group命令 group soil range z 1 2 ;定义group prop bulk 7.8e6 shear 3.0e6 coh 10e3 fric 15 ran group soil ;给group赋值 model null range group soil model elastic range group soil plot block group range group dam ;只显示dam range命令 range name trench x 0 1 y 0 4 z 0 2 model null range trench model null range x=2,4 y=2,6 z=5,10 range name Big_Brick x -3 3 y -2 2 z -1 1 model elastic range Big_Brick prop bulk 1e8 shear 1e8 range Big_Brick range name Layer1 plane dip 0 dd 0 ori 0 0 0 above range name Layer2 plane dip 0 dd 0 ori 0 0 0 below range cylinder end1 x1 y1 z1 end2 x2 y2 z2 radius r cylindrical range with one end of the cylinder axis (end1) at location (x1, y1, z1), the other end (end2) at location (x2, y2, z2), and with a cylinder radius of r '由(x1, y1, z1) 、(x2, y2, z2)两点确定旋转轴 macro命令 macro Sand 'bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35' macro Clay 'bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0' prop sand range Layer1 prop clay range layer2 macro Pt0 'p0 0 0 0'

(完整word版)FLAC3D模拟静水压力的变形模拟

FLAC3D静水压力的变形应变模拟模拟土体存在地下水位的情况下的应力应变。 通过弹塑性求解生成一个简单模型初始地应力场的过程.土体孔隙率为0。5,水位线从模型底部起为1m,计算中土体密度分水上与水下分别设置为饱和密度和天然密度,采用无渗流模式进行计算,其它计算条件见命令。 建模及命令语句如下: ；开始一个新的分析 N ;生成网格模型 gen zon bri size 1 1 2 ；设置强度参数 model m ;设置强度参数 prop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 ;固定z=0平面所有节点z向速度 fix z ran z 0 ;固定x=0平面所有节点x向速度 fix x ran x 0 ;固定x=1平面所有节点x向速度 fix x ran x 1 ；固定y=0平面所有节点y向速度 fix y ran y 0

;固定y=1平面所有节点y向速度 fix y ran y 1 ；设置土体饱和密度(水位线以下） ini dens 2000 ran z 0 1 ;设置土体干密度（水位线以上） ini dens 1500 ran z 1 2 ；设置土体竖向初始应力 ini szz —35e3 grad 0 0 20e3 ran z 0 1；设置土体y向初始应力 ini syy -17。5e3 grad 0 0 10e3 ran z 0 1；设置水下土体x向初始应力 ini sxx —17。5e3 grad 0 0 10e3 ran z 0 1；设置水上土体竖向初始应力 ini szz —15e3 grad 0 0 15e3 ran z 1 2 ;设置水上土体y向初始应力 ini syy —7.5e3 grad 0 0 7。5e3 ran z 1 2；设置水上土体x向初始应力 ini sxx —7.5e3 grad 0 0 7。5e3 ran z 1 2；设置初始孔隙水压力 ini pp 10e3 grad 0 0 -10e3 ran z 0 1 ；设置重力加速度 set grav 0 0 -10 ；按软件默认精度求解

FLAC3D入门基本知识

FLAC3D入门基本知识 FLAC3D一点知识点，仅以参考 4、id,cid的区别 id是指在整个结构中的编号，而cid是指在某一类比如说cable中的编号。拿cable 中的一个单元来说，它既有自己在整个结构中的cd,又有自己在cable中的cid 如果我设置了两个pile sel pile id=1 begin=(10.0, 1.0, 0.0) end=(10.0, 1.0, -10.0) nseg=5 sel pile id=2 begin=(10.0, 3.0, 0.0) end=(10.0, 3.0, -10.0) nseg=5 那么,id=1是不是代表第一根桩? 第一根桩分五段,cid=1~5,那么第二根桩是cid=6~10！ 5、什么情况下使用set large？ 初始应力平衡的时候，不能用large模式。在进行初始应力平衡时一定不要用！在进行大变形计算时，最好要用！！一般硬岩可以使用FLAC默认的小应变,如果是土体和软岩,用大应变 . 在做开挖的时候在进行原始应力平衡计算的时候是用小应变，后面的开挖以 及支护的时候选用大应变. 6、得到初始应力的方法： 方法、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。如在手册中给的第一个示例中就是这样做的。下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis. Gen zone brick size 6 8 8 Mode mohr Prop bulk 1e8 shear 0.3e8 fric 35

(完整word版)FLAC3D滑坡模拟

FLAC 3D滑坡模拟 一、源程序 ;Create Material Zones gen zone brick size 5 5 5 & pO (0,0,0) pl (3,0,0) p2 (0,3,0) p3 (0,0,5) & p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5) gen zone brick size 5 5 5 p0 (0,0,5) edge 5.0 group Material ;Create Bin Zones gen zone brick size 1 5 5 & p0 (4,1,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) & p6 add (3,0,5) p7 add (3,6,5) gen zone brick size 1 5 5 & p0 (6,1,5) p1 add (1,0,0) p2 add (0,5,0) & p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) & p6 add (1,0,5) p7 add (1,6,5) gen zone brick size 5 1 5 & p0 (1,4,0) p1 add (3,0,0) p2 add (0,3,0) & p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) & p6 add (5,2,5) p7 add (6,3,5) gen zone brick size 5 1 5 & p0 (1,6,5) p1 add (5,0,0) p2 add (0,1,0) & p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) & p6 add (5,0,5) p7 add (6,1,5) group Bin range group Material not ;Create n amed range synonyms range n ame=B in group Bin range n ame=Material group Material ;Assig n models to groups model mohr range Material model elas range Bin ;Create in terface eleme nts int 1 face ran pla ne ori (4,0,0) n or (-5,0,2) dist 0.01z (0,5) y (1,6) int 2 face ran pla ne ori (0,4,0) n or (0,-5,2) dist 0.01z (0,5) x (1,6) int 1 face ran x 5.9 6.1 y 1 6 z 5 10 int 2 face ran x 1 6 y 5.9 6.1 z 5 10 int 1 maxedge 0.55 int 2 maxedge 0.55 ;Move bin toward material ini x add -1.0 range Bin ini y add -1.0 range Bin ;Assig n properties

FLAC3D简介(word文档良心出品)

1.FLAC3D知识基本介绍SimWe 岩土工程结构的数值解是建立在满足基本方程（平衡方程、几何方程、本构方程）和边界条件下推导的。由于基本方程和边界条件多以微分方程的形式出现，因此，将基本方程近假发改用差分方程（代数方程）表示，把求解微分方程的问题改换成求解代数方程的问题，这就是所谓的差分法。差分法由来已久，但差分法需要求解高阶代数方程组，只有在计算机的出现，才使该法得以实施和发展。 FLAC3D（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。目前，FLAC 有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存（64K），所以，程序求解的最大结点数仅限于2000个以内。1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。因此，大大发护展了计算规模。FLAC3D是一个三维有限差分程序，目前已发展到V2.1版本。 FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。 FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发变形和移动（大变形模式）。FLAC3D采用的显式拉格朗日算法和混合-离散分区技术能够非常准确发模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。FLAC3D采用ANSI C++语言编写的。 FLAC3D有以下几个优点： 1 对模拟塑性破坏和塑性流动采用的是“混合离散法“。这种方法比有限元法中通常采用的“离散集成法“更为准确、合理。 2 即使模拟的系统是静态的，仍采用了动态运动方程，这使得FLAC3D在模拟物理上的不稳定过程不存在数值上的障碍。 3 采用了一个“显式解“方案。因此，显式解方案对非线性的应力-应变关系的求解所花费的时间，几互与线性本构关系相同，而隐式求解方案将会花费较长的时间求解非线性问题。面且，它没有必要存储刚度矩阵，这就意味着；采用中等容量的内存可以求解多单元结构；

(完整word版)FLAC3D5.0模型及输入参数说明

1.1模型参数代码 可参考manual中各个章节的command命令及说明,注意单位.用prop 赋值。

1.1.12

经典粘弹性模型的材料参数(Classical Viscoelastic （Maxwell Substance) –MODEL mechanical viscous） 1bulk弹性体积模量，K 2shear弹性剪切模量，G 3viscosity动力粘度，η 1.1.13粘弹性模型 粘弹性模型的材料参数（Burgers Model –MODEL mechanical burgers） 1bulk弹性体积模量，K 2kshear Kelvin弹性剪切模量，G K 3kviscosity Kelvin动力粘度，ηK 4mkshear Maxwell切边模量，G M 5mviscosity Maxwell动力粘度,ηM 1.1.14二分幂律模型 二分幂律模型的材料参数（Power Law –MODEL mechanical power） 1a_1常数,A1 2a_2常数,A2 3bulk弹性体积模量，K 4n_1指数，n1 5n_2指数,n2 6rs_1参考应力,σ1ref 7rs_2参考应力，σ2ref 8shear弹性剪切模量，G 1.1.15蠕变模型 蠕变模型材料参数（WIPP Model –MODEL mechanical wipp） 1act_energy活化能，Q 2a_wipp常数，A 3b_wipp常数,B 4bulk弹性体积模量,K 5d_wipp常数，D 6e_dot_star 临界稳定状态蠕变率, 7gas_c气体常数，R 8n_wipp指数,n 9shear弹性剪切模量，G 10temp温度，T 下列参数可以显示、绘图和通过fish访问 1e_prime累积主蠕变应变 2e_rate累积主蠕变应变率 Burger、蠕变组合材料模型的材料参数(Burgers—Creep Viscoplastic Model –MODEL mechanical cvisc) 1bulk弹性体积模量,K 2cohesion内聚力，c 3density密度，ρ 4dilation剪胀角，Ψ 5friction内摩擦角,Φ 6kshear Kelvin弹性剪切模量，G K

FLAC-3D深基坑的开挖与支护的命令流

FLAC D3深基坑的开挖与支护的命令流 一、实例工程 南宁地区地层属于河流阶地二元地层，广泛分布有较厚的圆砾层，国内尚无在类似地层条件下建设地铁基坑的经验，为此，可使用FLAC3D 对基坑开挖的全过程进行三维数值模拟，在对比实测数据的基础上，总结圆砾层中地铁车站深基坑的地下连续墙水平变形及周围地表沉降变形特征。 该基坑位于大学路与明秀路交叉路口处，沿大学东路东西向布置。车站基坑长465m，标准断面宽度为20.7m，为地下两层式结构，底板埋深为15.535m（相对地面），顶板覆土厚度大于3m。本工程主体建筑面积21163.6m2，主要结构形式为双柱三跨框架箱型结构。本工程所处的大学路为南宁市东西向的主要交通枢纽，车流量大，人流密集，地面条件复杂。 基坑施工采用明挖顺作法施工，围护结构为800mm厚地下连续墙+内撑（三道内支撑加一道换撑）的支护体系。第一道支撑采用钢筋混凝土支撑，尺寸为800×900mm，冠梁同时作为第一道钢筋混凝土支撑的围檩。第二、三道支撑及换撑使用钢支撑并施加预加力，直径为609mm，壁厚为t=16mm，斜撑段采用800×1000mm钢筋砼腰梁，其余为2×I45C 钢围檩。 二、模型建立 建模工作由两部分组成，实体模型部分，包括土体和地下连续墙；结构单元部分，包括混凝土支撑和钢支撑。 根据对称性原理，拟选取1/2 的实际工程尺寸进行分析。考虑到实际的基坑长度将近500m，根据以往的经验，选取全部长度的一半虽然能够得到满意的结果，但是由于中间部分的基坑基本处于同样的受力状态，这样会使大部分的计算长度变为重复的计算，降低了计算效率。根据初步计算结果和经验，最终确定的基坑尺寸为，宽度取基坑的最大宽度24m，开挖深度19m，基坑长度36m。根据地勘报告，合并相似土层，模型中共划分了7个土层。 在FLAC3D 中，围护结构可以用衬砌单元(liner)或实体单元模拟。根据Zdravdovi的研究，在二维平面基坑模拟中，分别采用实体单元和梁单元(相当于三维模型中的衬砌单元)计算所产生的墙体变形差别小于4%，而引起地表沉降的主要原因是围护结构变形造成的地层损失，可见上述两种方法计算结果的差别可忽略不计。相比于衬砌单元(liner)，实体单元物理模型清晰，参数相对较少且容易确定，因此本例围护结构采用实体单元模拟。 模型的范围必须要足够大以覆盖基坑开挖可能影响的区域，同时又要兼顾计算效率。根据现行国家和个地方规范的要求，参考已有的研究成果，影响范围为4 倍开挖深度，整个模型的尺寸为96×144×39（长×宽×高m）。为减少不必要的网格，网格划分原则为基坑附近密集，远处稀疏，模型共有节点10890个，单元9408个。模型如图所示。

(完整word版)flac命令流

1、FLAC3D常见命令: 1. FLAC3D是有限元程序吗? 答：不是！是有限差分法。 2。最先需要掌握的命令有哪些？ 答：需要掌握gen, ini, app，plo，solve等建模、初始条件、边界条件、后处理和求解的命令. 3。怎样看模型的样子? 答：plo blo gro可以看到不同的group的颜色分布 4. 怎样看模型的边界情况? 答：plo gpfix red 5。怎样看模型的体力分布？ 答：plo fap red 6。怎样看模型的云图? 答：位移：plo con dis (xdis，ydis，zdis) 应力：plo con sz (sy， sx, sxy, syz，sxz） 7。怎样看模型的矢量图？ 答：plo dis (xdis，ydis, zdis) 8。怎样看模型有多少单元、节点？ 答:pri info 9。怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择File/Print,将保存格式选择为jpe文件 10. 怎样调用一个文件？ 答：File/call或者call命令 10. 如何施加面力？ 答：app nstress

11. 如何调整视图的大小、角度？ 答:综合使用x, y, z, m，Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键 12. 如何进行边界约束？ 答：fix x ran (约束的是速度，在初始情况下约束等效于位移约束) 13. 如何知道每个单元的ID？ 答：用鼠标双击单元的表面，可以知道单元的ID和坐标 14。如何进行切片？ 答:plo set plane ori （点坐标) norm (法向矢量) plo con sz plane （显示z方向应力的切片) 15. 如何保存计算结果？ 答:save +文件名. 16. 如何调用已保存的结果？ 答：rest +文件名；或者File / Restore 17。如何暂停计算？ 答：Esc 18。如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令,用continue进行继续 19. 如何跳过某个计算步？ 答：在计算中按空格键跳过本次计算，自动进入下一步 20。Fish是什么东西? 答:是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学? 答：可以不用，需要的时候查Mannual获得需要的变量就可以了

flac3d5.0常用命令集锦

建 模 1、调用文件： ①文件与工程在同一个文件夹，只写文件名即可：If the called file is located in the same folder as the FLAC3D project file, then only the file name need be entered with the CALL command. ②不在同一个文件夹，全路径：Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat ). Undo ;撤销上一条命令 2、创建 旋转 缩放视图 3、建模命令 model mech mohr ;莫尔库伦模型 model mech elastic ;弹性模型 set grav 0, 0, -9.81 ; 重力加速度 negative z -direction. （垂直向下！常用的） 下下面面这这代代码码，，是是沿沿着着--y y 方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！ gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续

p1 10, -10, -20 ... p2 -10, 10, -20 ... p3 -10, -10, 0 plot zone gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体 p1 10, -10, -20 p2 -10, 10, -20 ... p3 -10, -10, 0 p4 10, 10, -20 ... p5 -10, 10, 10 p6 10, -10, 0 ... p7 10, 10, 10 plot current plot Plot01 plot clear plot zone Undo ;撤销命令 set logfile 127X1001.txt set log on truncate set log off list zone prin range x 0 1 y 0 1 z 0 1 ;显示指定围各单元的主应力，结果如下

(完整word版)FLAC3D原理..

2.2 三维数值模拟方法及其原理 2.2.1 FLAC3D工程分析软件特点 FLAC3D是由美国Itasca Consulting Group, Inc. 为地质工程应用而开发的连续介质显式有限差分计算机软件。FLAC即Fast Lagrangian Analysis of Continua 的缩写。该软件主要适用于模拟计算岩土体材料的力学行为及岩土材料达到屈服极限后产生的塑性流动，对大变形情况应用效果更好。 FLAC3D程序在数学上采用的是快速拉格朗日方法，基于显式差分来获得模型全部运动方程和本构方程的步长解，其本构方程由基本应力应变定义及虎克定律导出，运动平衡方程则直接应用了柯西运动方程，该方程由牛顿运动定律导出。 计算模型一般是由若干不同形状的三维单元体组成，也即剖分的空间单元网络区，计算中又将每个单元体进一步划分成由四个节点构成的四面体，四面体的应力应变只通过四个节点向其它四面体传递，进而传递到其它单元体。当对某一节点施加荷载后，在某一个微小的时间段内，作用于该点的荷载只对周围的若干节点（相邻节点）有影响。利用运动方程，根据单元节点的速度变化和时间，可计算出单元之间的相对位移，进而求出单元应变，再利用单元模型的本构方程，可求出单元应力。在计算应变过程中，利用高斯积分理论，将三维问题转化为二维问题而使其简单化。在运动方程中，还充分考虑了岩土体所具有的粘滞性，将其视作阻尼附加于方程中。 FLAC3D具有一个功能强大的网格生成器，有12种基本形状的单元体可供选择，利用这12种基本单元体，几乎可以构成任何形状的空间立体模型。 FLAC3D主要是为地质工程应用而开发的岩土体力学数值评价计算程序，自身设计有九种材料本构模型： （1）空模型（Null Model） （2）弹性各向同性材料模型（Elastic, Isotropic Model） （3）弹性各向异性材料模型（Elastic, anisotropic Model） （4）德拉克-普拉格弹塑性材料模型（Drucker-Prager Model） （5）莫尔-库伦弹塑性材料模型（Mohr-Coulomb Model） —29 —