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 :用R无量纲值确定收敛标准（默认10e-5）

格式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文件夹中