(完整版)FLAC3D常用命令
1. apply（缩写：app）
可用来定义边界条件及初始条件：
1）添加应力
格式1：apply szz -0.3395e6 range z -0.1 0.1
格式1：apply szz -0.3395e6 range group pile
格式3：apply szz -0.3395e6 range z -0.1 0.1 group pile
格式4：apply nstress 数值range z 2.9 3.1（或3）x 1 2 y 1 2
2）以一定速度施加位移边界
格式1：apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step
3）添加边界条件
格式：apply szz <数值> grad <梯度> range <范围>
示例1：apply szz -1e9 grad 0 0 8.3e5 range z 0 120
示例2：apply szz -0.6e6 range z 0.05 0.15 group pile
注：<数值>是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。
2. range(缩写：ran)
通过range功能，可以使命令作用在一定指定范围的目标上；如果一个命令没有使用range来确定范围，则命令对整个模型有效。
1)利用坐标指定一定的范围
格式1：range z 0 1
格式2：range z 2.9 3.1 x 1 2 y 1 2
2)利用分组来指定范围
格式：range group 1
3)以上两种的复合
格式：range z -0.1 0.1 group pile
4）利用id号来指定一定的范围
格式：range id 0 10
该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、
接触面和节点的编号。
例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。
5）命名一个范围（需要先命名这个范围）
格式：range name <自己起一个名字> <范围>
示例：range name intersected_zones x 5 8 y 3 7
引用这个范围的格式：内容+range nrange <自己起的那个名字>
示例：range nrange intersected_zone
6）命名一个柱形范围
格式：同4）
示例1：range name bhole cylinder end1 0 0 0 end2 -9.14 0 0 rad 0.46
示例2：range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 &
cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not 7）any和not
Any跟在range命令后表示“和”的意思，即使若干个定义的范围同时起效；not跟在range后表示“除了”的意思，即除此之外的范围
示例1：group wall range group 1 any group 2 any ;把组1和组2定义为新组，组名为wall
示例2：group wall range group 1 not ;把除了组1之外的所有部分定义为新组，组名为wall.
示例3：range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 &
cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not
3. set
用来设置重力、收敛标准、给某些变量赋值
1）设置重力
格式：set gravity 0 0 -10
2）设置收敛标准
格式1：set mech ratio
格式2：set mech force
以上两种格式也可以进行耦合：
示例：set mech force 0 ratio 1e-4
格式3：自定义收敛标准，比如某一关键点的位移或速度达到某一临界值（遇到这种情况进行补充）
3）在FLAC3D中给FISH语言中的变量赋值
格式：set @变量=value 如：set @hh=-22
4）设置大变形
格式：set large
4）设置监测历史变量时的监测间隔
格式1：set hist_rep <数目> ；默认的应该是10步一个间隔。
示例：set hist_rep 1
格式2 hist nstep <数目>
示例：hist nstep 5
4. history（缩写：hist）
设置监测点要进行监测的项目，不设置编号的话会从1自动进行编号。监测点必须用id或坐标这两种方法之一指定，新版本后需加add，不加add其实也可以。
1）监测整个模型的不平衡力
格式：hist add id <编号>unbal
2）设置节点监测项目
格式1：hist add gp zdis 0 0 6（监测节点z方向位移）
hist add gp zvel 0 0 6（监测节点z方向速度）
格式2：hist add id <给定该监测变量的编号> gp zdis <坐标>或者监测点id 示例1：hist add id 2 gp zdis 1 1 3
示例2：hist add id 1 gp zdis id 1 ；监测1号节点z方向的位移，这个监测
点编号是1
3）设置单元监测项目
格式1：hist add zone szz 0 0 -1(监测节点zz方向应力)
格式2：hist add id <给定该监测变量的编号> zone szz <坐标><坐标>或者监测单元id
示例1：hist add id 7 zone szz 1.5 1.5 2.5
4）设置监测历史变量时的监测间隔
格式1：set hist_rep <数目> ；默认的应该是10步一个间隔
示例：set hist_rep 1
格式2 hist nstep <数目>
示例：hist nstep 5
5）将计算结果写入到txt文件夹中
格式：hist write <监测变量的编号> skip <记录间隔> begin <监测开始的步数> end <监测结束步数> file 路径/文件名.txt
示例1：hist write 7 vs 8 file C:\Users\zhangzhao\Desktop\7-8.txt
示例2：hist write 7 file C:\Users\zhangzhao\Desktop\7.txt
示例3：hist write 1 2 3 4 5 skip 20 begin 8300 end 8400 file C:\Users\zhangzhao\Desktop\7.txt ;把1-5号监测点在8300到8400步之间的数据，按照每20步一个的记录间隔输出到7.txt文件中
5.model
用来指定模型的本构模型，以下格式中mech可以省略
格式1：model mech mohr ;定义全模型的本构模型
格式2：model mech mohr range group <组名>；
格式3：model mech null range group <组名>；
格式4：model mech null range z 0 1
6.property(缩写：pro)
给模型材料赋予参数
格式：pro bulk <值> shear <值> coh <值> fric <值> ten <值> range <范围>
示例1：prop bulk 5.6420e6 shear 1.3789e6 coh 9.77e3 fric 15.36 range group 1示例2：prop bulk 6e6 she 2.2e6 fric 10 coh 5000 range z 17.3 20.1
7.group
用来进行模型的分组
1）建模时直接将这部分命名为一个分组
格式1：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1
格式2：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 fill group 1 ；(对填充区域进行有别于外围材料的分组命名)
格式3：gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1 fill group 2 ；(外围区域命名为1，填充区域进行有别于外围材料的分组命名2)
2）对一定范围内的数据命名分组
格式：group 1 range < > ;用range定义范围
3）分开命令为一个组
示例：
group dam range x 5 7 z -5 0
group dam range id 201 any id 211 any id 231 any id 241 any
group dam range id 202 any id 212 any id 222 any id 232 any
以上区域命名了一个分组dam。
可以不同的区域命名为同一个组名，最后就会是一个组了。
8.attach
进行网格的连接，使用该命令时对各网格连接面上的单元尺寸有限制，要求它们之间的比率成整数倍
格式：attach face range < > ;range后跟定义范围的关键词
9.reflcet(缩写：ref)
以某对称轴生成镜像网格
格式1：gen zone ref ;以z=0平面为对称轴镜像生成网格
格式2：gen zone ref dip 90 ; 以y=0平面为对称轴镜像生成网格
格式3：gen zone ref dip 90 dd 90 ; 以x=0平面为对称轴镜像生成网格
格式4：gen zone ref dip 90 dd 90 range z 20 35 ; 可以加上范围，在z为20～35范围内以x=0平面为对称轴镜像生成网格。
其中：dip为平面在空间坐标系下沿着z轴负向与xy平面所成的夹角；dd为平面的dip方向角，即坐标系的xy平面从y轴正向顺时针旋转到指定方向的夹角。如果不注明代表值为0。
格式5：gen zone reflect normal 1 0 0 ori 0 0 0
格式：normal x y z origin x y z:其中normal为指定平面的法向方向，origin 为法线起点位置。这种方法比较直观，相比于用dd和dip定义平面，更容易记住，也不易出错。
10.fix固定/free放松
进行模型边界条件的设置
1）约束x/y/z方向的位移，其实约束的是速度
格式；fix <限制变形的方向> range <范围>
示例1：fix z range z -0.1 0.1
示例2：fix y range z -0.1 0.1
示例3：fix x z range z -0.1 0.1 group 1
示例4：fix x ；也可以不加范围，把所有的地方都固定
2)fix pp
格式：fix pp
给定孔隙水压力边界条件
3）给定温度边界条件
格式：fix temperature
11.initial(简写ini)
用来设置初始值
1）用来设置密度
格式：ini dens <值> range <范围>
示例1：ini dens 1830 range group 1ceng
示例2：ini dens 2000 range z 0 1
2）初始量归零
格式：ini <量> <值> range <范围>
示例1：ini xdis 0 ydis 0 zdis 0 xvel 0 yvel 0 zvel 0
示例2：ini szz 0 sxx 0 syy 0
3）初始状态为0
格式：ini state <状态>
示例：ini state 0
一般的程序地应力平衡后，初始化一般要进行塑性区、位移、速度清零
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
4）设置初始孔隙水压力
格式：ini pp <初始值> grad <梯度> range <范围>
示例：ini pp 30e3 grad 0 0 -10e3 range z 0 1
注：<初始值>是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。
5）设置初始应力
格式：ini szz <初始值> grad <梯度> range <范围>
示例：ini szz -35e3 grad 0 0 20e3 range z 0 1
注：<初始值>是梯度方向坐标0点的数值，可通过定义坐标范围的上下值与梯度计算得到。
6）将单元体移动
格式：ini <方向> add <距离> range <范围>
示例：ini zposition add -6.0 range group pile
12.save
将计算状态进行保存，会保存在当前文件下
1.常规
格式：save 文件名.sav
示例1：save 6-1.sav
2.自动循环
rest 8-4.sav
table 1 name load_settlement ;可要可不要，不要会自动从1编号
def add_load
p_gp=gp_near(2,1,3)
loop n (1,5)
app_load=n*(-1000e3)
filename='7-6_add_step'+string(n)+'.sav'
command
app nstress @app_load range z 2.9 3.1 x 1 2 y 1 2
solve
save @filename
endcommand
xtable(1,n)=app_load
ytable(1,n)=gp_zdisp(p_gp)
endloop
end
@add_load
save 8-6.sav
13.solve
对模型进行计算
1、格式1：solve ;直接计算，默认最大平衡应力比小于或等于1e-5时，计算结果获得稳定值
2、格式2：solve ratio <计算精度值>；设置最大不平衡应力比，计算达到这个条件即终止
示例：solve ratio 1e-5
3、格式3：step <步数> ；设置计算终止步数
4、渗流求解
格式1：solve ；得到渗流平衡状态下的孔压场，渗流平衡的概念类似于
力学平衡，表示节点上的不平衡流量比达到了程序设置
的最小值。
格式2：step <步数> ；在计算中完成相应的步数以后程序自动结束。
格式3：solve age <时间> ；时间是流动时间，这个时间并不一定是真实的
时间，只有当所有的流体参数都为真实时才是
真实的时间。
14.plot
详见《FLAC3D快速入门及简单实例》第四章FLAC3D图形显示
进行绘图，进行绘图前最好指定一个绘图文件
格式：plot create view(或者plot) <名字>
1）绘制模型单元、建立的界面等
格式1：plot zone
示例1：plot zone trans 75 ;绘制透明度为75的单元体
示例2：plot zone fill off ;单元体不进行填充
示例3：plot zone trans 75 fill off ；以上两种的组合
示例4：plot zone colorlist red
示例5：plot zone colorlist red range nrange <范围名称>
示例6：plot interface colorlist red
示例程序：
new
gen zone brick size 10 10 10
range name intersected_zones x 5 8 y 3 7
range name union_zones x 5 8 y 3 7
plot create view intersected_zones ;（创建一个intersected_zones的绘图纸）plot add zone colorlist yellow trans 80
plot add zone colorlist red range nrange intersected_zone
plot add axes
plot create view union_zones ;(也可通过plot create plot union_zones) plot add zone colorlist yellow trans 80
plot add zone colorlist red range nrange union_zones
plot add axes
格式2：plot interface
2）绘制云图（应力、变形，节点与单元的量都可绘制）
格式1：plot zcon <表示应力的量>
示例1：plot zcon szz
Plot zcon sxy
格式2：plot con <表示变形的量>
格式3：plot zcon <表示的量> ramp grey ;灰度显示
示例1：plot con zdis
3）绘制块云图（一般只表示应力的量）
格式：plot bcon <量>
示例1：plot bcon szz
示例2：plot bcon pro bulk
4）绘制矢量图（最好结合1）内容带上网格）
格式：plot dis；绘制变形矢量图（）
示例：plot dis scale 0.1 ;scale后跟的数值越大，矢量图的箭头越大。
5）绘制监测变量曲线图
格式1：plot hist <给监测变量设置的id编号>
示例：plot hist 2
格式2：plot hist <编号1> vs <编号2> ；两个监测变量之间的关系
示例：plot hist 7 vs 8 (历史变量7位横坐标、8位纵坐标)
6）plot add 的应用
格式：plot add +内容；把要绘制的图形绘制在前面的图形上示例：plot zone trans 75 fill off
Plot add dis
相当于在单元图形上绘制了变形矢量图
7）按组显示
格:1：plot zgroup range <范围>
示例：plot zgroup range z 0 5
格式2：plot zone range group 1
其他常见的plot命令
Plot gpfix
Plot gp
Plot boundary
Plot zone trans 75
Plot clear
Plot interface colorlist red
8）绘制表格数据
格式1：plot table id style both
格式2：plot table id style line
格式3：plot table id style mark
示例1：plot table 1 style both
9）绘制孔压云图
节点：plot con gpp
单元：plot zcon zpp
15.restore(简写为rest)
调用保存好的sav文件，并在此状态的基础上进行下面的操作格式：rest 文件名.sav
示例：rest 6-1.sav
16.state
单元状态
1）显示单元状态
格式：plot zone colorby state
2）输出单元状态
格式：print zone state
17.print
输出计算结果
1）输出节点、单元数据
格式1：print zone/gp ;输出单元数据
格式2：print gp position
格式3：print zone stress
格式4：print gp dis
格式5：print gp dis range id 517 any 533 any
格式6：print gp pos range id 59 ;id号为59的节点坐标信息
2）输出接触面的数据
格式1：print interface ;接触面数据输出
格式2：print interface <编号> stress ；输出接触面上的nstress/stress等格式3：print interface <编号> dis
格式4：print interface <编号> position
格式5：print interface <编号> prop
3）输出监测变量
格式1：print hist ;输出监测点信息
格式2：
2）将计算结果输入到指定的文件(也可参考4hist)
格式：（set log on
Set logfile C:\Users\zhangzhao\Desktop\6-3.log
Print zone stress
Set log off）
以上的程序将单元应力数据输出到6-3文件中
3）输出表格数据
格式：print table id
示例：print table 1
18.plane
创建切片，显示剖面图
Oringin/normal和dip/dd在定义平面上的作用是等效的
1）建立剖面
最后先建立一个绘图纸
格式：plot create view <绘图纸名字>
2）创建剖面
格式：plot cut add plane ori <坐标> norm <方向矢量>
示例：plot cut add plane ori 1.5 1.5 1.5 norm 0 1 0
3）绘制剖面上的内容
格式：plot add <内容> (plane on) ;括号内可省略
示例1：plot add zone fill off ;绘制剖面上的网格，建议先绘制这个。
示例2：plot add zcon szz
示例3：plot add dis
3）绘制剖面外的内容
格式：plot add <内容> plane b(f) b代表与norm反向一侧，f代表正向侧。
示例1：plot add zone plane b
示例2：plot add zone fill off plane f
示例程序：
plot create view plane4
plot cut add plane ori 1.5 1.5 1.5 norm 0 1 0
plot add zone fill off plane o
plot add zcon szz plane on
plot add zone plane b
plot add zone fill off plane f
19.输出图片
将界面上显示的内容以图片形式显示出来
1）输出ps照片
格式：plot post view <要输出的绘图纸名称> file 文件路径/名称.格式
示例：plot post view plot04 file C:\Users\zhangzhao\Desktop\3.bmp
2）输出bitmap(位图)图片注：常用格式bmp png
格式：plot bitmap view <要输出的绘图纸名称> file 文件路径/名称.格式
示例：plot bitmap view plot04 file C:\Users\zhangzhao\Desktop\7.bmp
以上示例中的view plot04是可以省略的，但是输出的就为当前界面的内容。
20.delete
删除指定范围内的网格
格式：delete zone range <>
Delete和model null是不一样的，delete是将网格信息删除，model null是只是将网格赋予了空属性，即乘以一个很小的数，后面还可以通过model将这些属性恢复。
21.interface(int)
建立接触关系，详见《FLAC3D快速入门及简单实例》第9章9.8节、《FLAC/FLAC3D基础与工程实例》第6章
1）移来移去法
格式：Interface
示例1：interface 5 face range cylinder end1 0 0 0 end2 0 0 -15.1 radius 0.16 &
cylinder end1 0 0 0 end2 0 0 -15.1 radius 0.14 not 多加上0.1,本来是0-15
示例2：interface 2 face range cylin end1 (2.,2.,-19.9) end2 (2.,2.,-20.1) radius .26多加上0.1,本来半径是0.25
2）切割模量法
关键命令：gen separation group 1
Int 1 wrap group1 group2
示例程序：
new
gen zone brick size 3 3 3
group 1 range x 1 2 y 1 2 z 1 2
group 2 range group 1 not
gen seperate group 1 ；将group 1 与主体模型分开int 1 wrap 1 2 ；建立接触面
plot interface
3)导来导去法
示例程序：
new
gen zone brick size 3 3 3
group 2 range x 1 2 y 1 2 z 1 2
group 1 range group 2 not
save 1.sav
delete zone range group 1
interface 1 face range x 0.9 1.1 y 1 2 z 1 2
interface 1 face range x 1.9 2.1 y 1 2 z 1 2
interface 1 face range x 1 2 y 0.9 1.1 z 1 2
interface 1 face range x 1 2 y 1.9 2.1 z 1 2
interface 1 face range x 1 2 y 1 1 z 0.9 1.1
interface 1 face range x 1 2 y 1 1 z 1.9 2.1
save 2.sav
rest 1.sav
delete zone range group 2
expgrid C:\Users\zhangzhao\Desktop\group2.flac3d rest 2.sav
impgrid C:\Users\zhangzhao\Desktop\group2.flac3d
plot interface
plot add zone fill off
4）设置界面参数
格式：interface
22.expgrid与impgrid
网格的导出与导入
格式：expgrid <存储路径>\文件名.flac3d
impgrid <存储路径>\文件名.flac3d
示例：expgrid C:\Users\zhangzhao\Desktop\group2.flac3d
impgrid C:\Users\zhangzhao\Desktop\group2.flac3d
23.list
1)打印节点、单元状态
格式1：list zone或者print zone
格式2：list gp position 或者print gp position
2）打印单元数据等
格式1：list zone stress
格式2：list zone prop bulk/shear……
3）打印有关节点数据
格式1：list gp
格式2：list gp position
格式3：list gp dis
4）打印有关fish函数
格式1：list @var ;如果是函数变量，计算函数，并输出数值
格式2：list fish ;按其名称的字母顺序输出所有的当前值，可用来检查5）输出表格数据
格式：list table id
示例：list table 1
6）孔压的输出
节点：print/list gp pp
单元：print/list zone pp
后面可以加上范围
24.generate
详细参考《FLAC3D快速入门及简单实例》第40页、《FLAC/FLAC3D基础与工程实例》第5章
25.pause/continue
这是两个控制命令流（.txt文件）读入进程的命令。Pause的作用是让程序在此处暂停，可以在命令行中输入continue让程序继续运行。
当命令流有问题时但却不知道在哪里的时候，可以通过在不同的位置插入pause来确定问题的所在。
26.structure elements
结构单元，建立方法详见《FLAC/FLAC3D基础与工程实例》第157页。
1）创建线性结构单元（beam cable pile）
线性结构单元是通过sel来创建的
格式1：sel pile id 1 begin 0 0 0 end 0 0 10 nseg 4
格式2：示例程序如下
new
sel node id 1 0 0 0
sel node id 2 2 0 0
sel node id 3 4 0 -1
sel node id 4 5 0 -2
sel beamsel id 1 cid 1 node 1 2
sel beamsel id 1 cid 2 node 2 3
sel beamsel id 1 cid 3 node 3 4
plot sel geometry
plot add ax
2）创建面结构单元（shell geogrid liner）
面结构单元也是sel来创建的
格式：
3）显示
格式1：plot sel geometry ;将所有单元都显示出来
格式2：plot sel pile ;显示某种特定的单元
格式3：plot sel geometry shrink 0.9 marktype node colorby id ;这种格式只能应用在格式1geometry中，因为marktype只有在这种条件下有。
格式4：plot sel pile shrink 0.9 line width 4 colorby id ;这种格式只能应用在格式2pile中，因为line width 4只有在这种条件下有。
4）结构单元之间的连接
通过sel link 命令进行设置
在实体单元内部或表面创建结构单元时，程序会自动对所有结构节点创建了node-zone连接。
5）结构单元的边界与初始条件
结构单元所有的边界条件和初始条件都是通过结构节点来设置的，采用sel node 命令来实现，详见《FLAC/FLAC3D基础与工程实例》第162页
下面6）7）8）内容详见《FLAC3D快速入门及简单实例》7.8节。
6）结构单元的赋值问题
格式：sel type peop <参数> ；type可以是beam/cable/pile/shell/geogrid/liner 示例：sel pile prop xcarea 6e-3 xciz 200e-6 xciy 200e-6 xcj 0
7）监测变量
详见《FLAC3D快速入门及简单实例》7.8节
8）输出各个节点的变量
格式1：list sel beam type ;type为具体的监测项目
格式2：list sel cable type ;type为具体的监测项目
格式1：list sel pile type ;type为具体的监测项目
格式1：list sel recover type ;shell 单元，type为具体的监测项目
格式1：list sel geogrid type ;type为具体的监测项目
格式1：list sel liner type ;type为具体的监测项目
27.config
格式1：config dyn ;打开动力响应
格式2：config fluid ;打开渗流模式
28.显示某些赋值的量
格式1：plot zcon/bcon prop bulk
格式2：list zone prop bulk/shear
29.call
格式1：call find-add.f3fis
格式2：
注：1. FLAC3D及FISH语言中是不分大小写的，全部会被转换为小写；
2. FLAC3D中，速度以位移除以时步表示；
3. 在FLAC3D总不进行设置的量代表为0，不设置范围的代表全模型。
常用的计算步骤（也可进行局部的调整，如开挖与回填问题）：
1、画网格
2、材料本构模型与材料参数（先本构，再参数）
3、边界条件与初始条件
4、计算（初始平衡与后续计算）
5、后处理（处理数据、作图）
可以用来检查编程正确性的命令
1、attach
2、gen merge
3、plot bcon szz
FLAC3D基础知识介绍
FLAC 3D 基础知识介绍 一、概述 FLAC（Fast Lagrangian Analysis of Continua ）由美国Itasca 公司开发的。目前，FLAC 有二维和三维计算程序两个版本，二维计算程序V3.0 以前的为DOS 版本，V2.5 版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。1995 年，FLAC2D 已升级为V3.3 的版本，其程序能够使用护展内存。因此，大大发护展了计算规模。FLAC3D是一个三维有限差分程序，目前已发展到V3.0 版本。 FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。 FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。FLAC3D 采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围 的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析 方法，它可以模拟岩土或其他材料的三维力学行为。三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所 谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。三维快速 拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。 FLAC-3D(Three Dimensional Fast Lagrangian Analysis of Continua)是美国Itasca Consulting Goup lnc 开发的三维快速拉格朗日分析程序，该程序能较好地模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学行为，特别适用于分析渐进破坏和失稳以及模拟大变形。它包含10种弹塑性材料本构模型，有静力、动力、蠕变、渗流、温度五种计算模式，各种模式间可以互相藕合，可以模拟多种结构形式，如岩体、土体或其他材料实体，梁、锚元、桩、壳以及人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩、界面单元等，可以模拟复杂的岩土工程或力学问题。 FLAC3D采用ANSI C++语言编写的。 二、FLAC3D的优点与不足 FLAC3D有以下几个优点: 1对模拟塑性破坏和塑性流动采用的是混合离散法。这种方
Flac3D命令--完整经典版
实例分析命令： 1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转 Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,-9.81（重力加速度） plot add block group red yellow 把在group 中的部分染成红色和黄色 plot add axes black 坐标轴线为黑色；print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b （设置初始密度，有时不同层密度不同）；ini ……（设置初始条件）；fix ……（固定界面） set plot jpg ；set plot quality 100 ；plot hard file 1.jpg 图像输出（格式、像素、名称） plot set magf 1.0视图的放大倍数为1.0；plo con szz z 方向应力云图 2. ini z add -1 range group one 群one 的所有单元，在z 方向上向下移动1m ；然后合并 命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体，1e-5是容差 3. (基坑开挖步骤)：Step 1: create initial model state （建立初始模型）Step 2: excavate trench （开挖隧道） 4. group Top range group Base not 定义（群组Base 以外的为）群组Top 5. plot blo gro 使得各个群组不同颜色显示 6. （两个部分间设置界面；切割法）：gen separate Top 使两部分的接触网格分离 为两部分；interface 1 wrap Base Top 在（Base 和Top ）这两部分之间添加接触单元；plot create view_int 显示，并创建标题view_int ；plot add surface 显示表面；plot add interface red 界面颜色红色 7. （简单的定义函数及运行函数）new ；def setup 定义函数setup ；numy = 8定义常 量numy 为8；depth = 10.0 定义depth 为10；end 结束对函数的定义；setup 运行函数setup 8. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。 9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度？— 10. 渐变应力施加：apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0；施加法向应力：apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0 11. d ip dd 确定平面位置使用：（纠结） 12. p rint gp position range id=14647 输出节点坐标 13. a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中，z 为变量，所以xx σ为：65=-1010+10xx z σ?? ；原点（0,0,0） 14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上，z=7，10 这一部分在x 方向的约 束（可以在此处产生破坏） 15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下： =3(1-2v)E K G=2(1+v) E 16. 一般而言，大多数问题可以采用FLAC 3D 默认的收敛标准（或称相对收敛标准），即当体 系最大不平衡力与典型内力的比率R 小于定值10-5；（也可由用户自定义该值，命令：
FLAC3D基础知识介绍
FLAC 3D基础知识介绍 一、概述 FLAC(Fast Lagrangian Analysis of Continua)由美国Itasca公司开发的。目前,FLAC有二维与三维计算程序两个版本,二维计算程序V3、0以前的为DOS版本,V2、5版本仅仅能够使用计算机的基本内存64K),所以,程序求解的最大结点数仅限于2000个以内。1995年,FLAC2D已升级为V3、3的版本,其程序能够使用护展内存。因此,大大发护展了计算规模。FLAC3D就是一个三维有限差分程序,目前已发展到V3、0版本。 FLAC3D的输入与一般的数值分析程序不同,它可以用交互的方式,从键盘输入各种命令,也可以写成命令(集)文件,类似于批处理,由文件来驱动。因此,采用FLAC程序进行计算,必须了解各种命令关键词的功能,然后,按照计算顺序,将命令按先后,依次排列,形成可以完成一定计算任务的命令文件。 FLAC3D就是二维的有限差分程序FLAC2D的护展,能够进行土质、岩石与其它材料的三维结构受力特性模拟与塑性流动分析。调整三维网格中的多面体单元来拟合实际的结构。单元材料可采用线性或非线性本构模型,在外力作用下,当材料发生屈服流动后,网格能够相应发生变形与移动(大变形模式)。FLAC3D采用的显式拉格朗日算法与混合-离散分区技术,能够非常准确的模拟材料的塑性破坏与流动。由于无须形成刚度矩阵,因此,基于较小内存空间就能够求解大范围的
三维问题。 三维快速拉格朗日法就是一种基于三维显式有限差分法的数值分析方法,它可以模拟岩土或其她材料的三维力学行为。三维快速拉格朗日分析将计算区域划分为若干四面体单元,每个单元在给定的边界条件下遵循指定的线性或非线性本构关系,如果单元应力使得材料屈服或产生塑性流动,则单元网格可以随着材料的变形而变形,这就就是所谓的拉格朗日算法,这种算法非常适合于模拟大变形问题。三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程,并应用了混合单元离散模型,可以准确地模拟材料的屈服、塑性流动、软化直至大变形,尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。 FLAC-3D(Three Dimensional Fast Lagrangian Analysis of Continua)就是美国Itasca Consulting Goup lnc开发的三维快速拉格朗日分析程序,该程序能较好地模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学行为,特别适用于分析渐进破坏与失稳以及模拟大变形。它包含10种弹塑性材料本构模型,有静力、动力、蠕变、渗流、温度五种计算模式,各种模式间可以互相藕合,可以模拟多种结构形式,如岩体、土体或其她材料实体,梁、锚元、桩、壳以及人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩、界面单元等,可以模拟复杂的岩土工程或力学问题。 FLAC3D采用ANSI C++语言编写的。 二、FLAC3D的优点与不足
FLAC3D常用命令
1. apply（缩写：app） 可用来定义边界条件及初始条件： 1）添加应力 格式1：apply szz -0.3395e6 range z -0.1 0.1 格式1：apply szz -0.3395e6 range group pile 格式3：apply szz -0.3395e6 range z -0.1 0.1 group pile 格式4：apply nstress 数值range z 2.9 3.1（或3）x 1 2 y 1 2 2）以一定速度施加位移边界 格式1：apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step 3）添加边界条件 格式：apply szz <数值> grad <梯度> range <围> 示例1：apply szz -1e9 grad 0 0 8.3e5 range z 0 120 示例2：apply szz -0.6e6 range z 0.05 0.15 group pile 注：<数值>是梯度方向坐标0点的数值，可通过定义坐标围的上下值与梯度计算得到。 2. range(缩写：ran)
通过range功能，可以使命令作用在一定指定围的目标上；如果一个命令没有使用range来确定围，则命令对整个模型有效。 1)利用坐标指定一定的围 格式1：range z 0 1 格式2：range z 2.9 3.1 x 1 2 y 1 2 2)利用分组来指定围 格式：range group 1 3)以上两种的复合 格式：range z -0.1 0.1 group pile 4）利用id号来指定一定的围 格式：range id 0 10 该命令后跟起始id和结束id，这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。 例：model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。 5）命名一个围（需要先命名这个围） 格式：range name <自己起一个名字> <围> 示例：range name intersected_zones x 5 8 y 3 7
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
[实用参考]Flac3d-5.0常用命令集锦.doc
建模 1、调用文件： ①文件与工程在同一个文件夹，只写文件名即可：Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlyt hefilenameneed beenteredwiththe CALL command. ②不在同一个文件夹，全路径：Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfol der\file.dat). Undo;撤销上一条命令 2、创建旋转缩放视图 3、建模命令 modelmechmohr;莫尔库伦模型 modelmechelastic;弹性模型 setgrav0,0,-9.81;重力加速度negative z-direction.（垂直向下！常用的） 下下面面这这代代码码，，是是沿沿着着--y y方方向向的的重重力力加加速速度度，，注注意意区区别别！！！！！！！！ genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续 p110,-10,-20... p2-10,10,-20... p3-10,-10,0 plotzone
genzonebricksize6,8,8p0-10,-10,-20...;不规则六面体 p110,-10,-20p2-10,10,-20... p3-10,-10,0p410,10,-20... p5-10,10,10p610,-10,0... p710,10,10 plotcurrentplotPlot01 plotclear plotzone Undo;撤销命令 setlogfile127G1001.tGt setlogontruncate setlogoff listzoneprinrangeG01y01z01;显示指定范围内各单元的主应力，结果如下 Hist命令： ①命令编号按顺序从1开始：eachhistoryisnumberedsequentiallyfrom1asitisenteredviathe HISTORY co mmand. ②查找显示所有的his命令：ReturntotheFlac3D>promptandtype listhist foralistingofthehistoriesandtheircorrespondingnumbers. histnstep5;每5步记录1次。默认是10步记录1次
FLAC3D常见命令与使用技巧
FLAC3D常见命令与使用技巧 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.如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步
Flac3D常见问题整理
1.1常见问题及其解答Gen separate 不能被识别答：原因是FLAC3D版本不行，我用3.0的版本不能。 1. FLAC3D是有限元软件吗？答：不是，是有限差法软件。 2. FLAC3D最先需要掌握的命令有哪些？答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3. 怎样看模型的样子？答：plo blo gro可以看到不同的group的颜色分布。 4. 怎样看模型的边界情况？答：plo gpfix red sk 5. 怎样看模型的体力分布？答：plo fap red sk 6. 怎样看模型的云图？答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz) 7. 怎样看模型的矢量图？答：plo dis (xdis, ydis, zdis) 8. 怎样看模型有多少单元、节点？答：print info 9. 怎样输出模型的后处理图？答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpg文件。 10. 怎样调用一个文件？答：使用菜单File/call 或者call 命令。 11. 如何施加面力？答：app nstress ran 12. 如何调整视图的大小、角度？答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键。 13. 如何进行边界约束？答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束） 14. 如何知道每个单元的ID？答：使用鼠标双击单元的表面，可以知道单元的ID和坐标。 15. 如何进行切片？答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 16. 如何保存计算结果？答：save filename(文件名可自定义) 17. 如何调用已保存的结果？答：使用菜单File/call或者命令rest filename(文件名可自定义)。 18. 如何暂停计算？答：运行中使用Esc命令。 19. 如何在程序中进行暂停，并可恢复计算？答：在命令中加入pause命令，键入continue命令后可恢复计算。 20. 如何跳过某个计算步？答：在计算中按空格键可跳过本次计算，自动进入下一步。 21. FISH是什么？答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能。 22. FISH是否一定要学？答：可以不用，需要的时候查Manual获得需要的变量就可以了。 23. FLAC3D允许的命令文件格式有哪些？答：只要是符合FLAC3D格式要求的文本文件，无论是什么后缀名，都可以为FLAC3D调用。 24. 如何调用一些可选模块？答：使用命令config dyn (fluid, creep, cppudm)。 25. 如何使用gauss_dev对符合高斯正态分布的材料参数进行赋值？答：假定某材料的摩擦角均值为40度，标准差是2，则命令如下：prop friction 40 gauss_dev 2 26. FISH函数中是否能调用“.sav”文件？答：不能。FLAC3D中规定，new和restore命令不允许出现在FISH函数中，因为new和restore 命令会将原有存储信息清除掉。 27. initial 与apply 有何区别？答：initial初始化命令，如初始化计算体的应力状态等；apply边界条件限制命令，如施加边界的力、位移等约束等。initial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。 28. FLAC3D动力分析中是如何计算永久变形的？答：FLAC3D采用动态运动方程求解动力方程，因此采用弹塑性本构模型可以计算永久变形。而土动力学常用的粘弹性模型由于没有考虑土体的塑性，因此不能计算永久变形。 29. 对于初学者而言，是学习FLAC还是FLAC3D？答：FLAC有较好的图形化操作界面，而FLAC3D目前只能通过命令流来操作，从学习难度上来说，FLAC要简单一些，不过复杂的三维问题还是需要使用FLAC3D才能解决。FLAC和FLAC3D的某些命令和分析方法类似，读者在学习过程中可以相互借鉴。 30. interface建模命令中的dist关键词是否表示接触面的厚度？答：FLAC3D 中的interface 是没有厚度的，dist 关键词表示的是接触面建模时选择范围时的容差，表示该范围内的“面”上将被赋予interface 单元。 31. 初始应力场计算中位移场和速度场是否都要清零？答：是的。一般，FLAC和FLAC3D中位移场和速度场的清零命令都是同时使用的。 32. 加了fix边界，再使用apply施加应力边界有效吗？答：无效。fix和apply都是边界条件，两者不能混用，fix的作用是固定节点的速度，只要用户不更改这个速度，在计算中都会保持不变。 33. solve age后面跟随的时间是真实的时间吗？答：FLAC和FLAC3D在动力、渗流、流变模式下才有真实的时间，时间的单位默认为秒，也可以根据读者使用的量纲进行调整。
flac3d常用命令
1、最先需要掌握的命令有哪些？ 答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 2、怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件。 3、怎样调用一个文件？ 答：File/call或者call命令 4、如何施加面力？ 答：app nstress 5、如何调整视图的大小、角度？ 答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键。 6、如何进行边界约束？ 答：fix x ran （约束的是速度，在初始情况下约束等效于位移约束）。 7、如何知道每个单元的ID？ 答：用鼠标双击单元的表面，可以知道单元的ID和坐标。 8、如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 9、如何保存计算结果？ 答：save +文件名 10、如何调用已保存的结果？ 答：rest +文件名；或者File / Restor 11、如何暂停计算？ 答：Esc 12、如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令，用continue进行继续。 在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正，而不是等到结果出来。 13、如何跳过某个计算步？ 答：在计算中按空格键跳过本次计算，自动进入下一步 14、Fish是什么东西？Fish是否一定要学？
答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能。Fish可以不用学，需要的时候查Mannual获得需要的变量就可以了。 15、FLAC3D允许的命令文件格式有哪些？ 答：无所谓，只要是文本文件，什么后缀都可以。 16、如何调用一些可选模块？ 答：config dyn (fluid, creep, cppudm) 17、如何在圆柱体四周如何施加约束条件？ 可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 第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 实用标准文档 1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 1、怎样查看模型？ 答：plot grid 可以查看网格，plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件？ 答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2 1.1常见问题及其解答 Gen separate 不能被识别 答：原因是FLAC3D版本不行，我用3.0的版本不能。 1. FLAC3D是有限元软件吗？ 答：不是，是有限差法软件。 2. FLAC3D最先需要掌握的命令有哪些？ 答：需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。 3. 怎样看模型的样子？ 答：plo blo gro可以看到不同的group的颜色分布。 4. 怎样看模型的边界情况？ 答：plo gpfix red sk 5. 怎样看模型的体力分布？ 答：plo fap red sk 6. 怎样看模型的云图？ 答：位移：plo con dis (xdis, ydis, zdis) 应力：plo con sz (sy, sx, sxy, syz, sxz) 7. 怎样看模型的矢量图？ 答：plo dis (xdis, ydis, zdis) 8. 怎样看模型有多少单元、节点？ 答：print info 9. 怎样输出模型的后处理图？ 答：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpg文件。 10. 怎样调用一个文件？ 答：使用菜单File/call或者call命令。 11. 如何施加面力？ 答：app nstress ran <……> 12. 如何调整视图的大小、角度？ 答：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键。 13. 如何进行边界约束？ 答：fix x ran <……>（约束的是速度，在初始情况下约束等效于位移约束） 14. 如何知道每个单元的ID？ 答：使用鼠标双击单元的表面，可以知道单元的ID和坐标。 15. 如何进行切片？ 答：plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 16. 如何保存计算结果？ 答：save filename(文件名可自定义) 17. 如何调用已保存的结果？ 答：使用菜单File/call或者命令rest filename(文件名可自定义)。 18. 如何暂停计算？ 答：运行中使用Esc命令。 19. 如何在程序中进行暂停，并可恢复计算？ 答：在命令中加入pause命令，键入continue命令后可恢复计算。 gen zone brick gro 1 ratio 1.1 1 1 size 15 1 10 p0 0 0 19.5 p1 20 0 19.5 p2 0 1 19.5 p3 0 0 29.5 建模. mo null range group 7 开挖 attach face检查节点 gen merge 1e-3融合节点 plo blo gro;按组画图 plo blo den;按密度画图 plo con zdis;Z向位移画图 plo con sz;Z向应力画图 set mech ratio 1e-3;设置计算收敛 plo con ssi disp vel ;剪应变增加滑移线 pl add disp;位移矢量图合成位移 plo sel recoer;可以显示应力 ini y add 2 在地壳每增加100m深度，压力将增加约2.75Mpa his gp disp id 1851 pl con disp; his write 7 vs step begin 1 end 10000000000 file 顶板.txt；输出记录数据 pl bl st n; pl bl st plo fap red sk;显示体力 plotitems------add-----vetor----apply body force plo add ve apply body force ini szz 1000 range z 13 14 x -1.5 1.5 y 0 15 施加局部调整力;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;出图;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; set plot bmp plot set mag 2.44 rotation 0 0 0 plot set back white plo con sz plo add sel geom pile plot contour ssi ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;出图;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 保存及调用 sol save kaiwa1.sav new res kaiwa1.sav 2-1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x -.1 .1 fix x ran x 2.9 3.1 fix y ran y -.1 .1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3-1 ;-------------------------------------工程信息 ;Project Record Tree export ;Title:Simple test ;---------------------------------计算第一步 ;... STATE: STATE1 .... config grid 10,10 model elastic group 'User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9.81 history 999 unbalanced solve save state1.sav ;----------------------------------计算第二步 ;... STATE: STATE2 .... initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete 'null' history 1 xdisp i=4, j=11 solve save state2.sav ;--------------------------------绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis-A plot hold history 1 lineFLAC3D 实例命令流1
FLAC3D命令流(整理版)
FLAC3D命令流(整理版)
FLAC3D错误提示解决办法
FLAC3D 3.0 C 常用命令
flac3D基础与工程实例命令