FLAC3D命令流(整理版)

1、怎样查看模型？

答：plot grid 可以查看网格，plot grid num 可以查看节点号。

2、请问在圆柱体四周如何施加约束条件？

答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 1

3、怎么能把一个PLOT 的图像数据导出来以便用其他软件绘图？答：用set log on 命令，把数据导出来，转到excel 里处理一下，然后用surfer或者什么作图软件绘制就行了。

4、用命令建立模型后，如何显示点的坐标？

答：使用plo bio gro gpnum on 命令

5、关于gauss_dev对性质进行高斯正态分布的问题？答：根据手册上的说明：下面的命令设定一个平均摩擦角为40 度，标准方差是± 5％。则命令如下：prop friction 40 gauss_dev 2 问题：请问gauss_dev 2 中的2 是如何计算的？如果把± 5％改为± 10%，则命令应如何写？40 X 5% = 2

6、reflect 问题

问：gen zone radbrick &p0 （0,0,0） p1 （10,0,0） p2 （0,10,0） p3 （0,0,10） &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill

plot surf

gen zone reflect dip 0 dd 90 （对xy 面做镜像）

gen zone reflect dip 90 dd 90 （对yz 面做镜像）

（1）dd 表示y 轴正向顺时针到那条射线的夹角，dip 表示对称参照

面与xy 平面的夹角，对称参照面与xy 平面的夹角在xy 平面的投影是一条射线。首先应该按照dd 的方向大体确定这个面的朝向，dd 指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy 平面上的投影的夹角，然后再确定dip，dip 指的是从xy 平面转向所要确定的平面的角．（及z 轴负方向转过角度）

（2）命令：reflect normal xyz origin xyz （根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：normal -1 0 0 与normal 1 0 0 的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。

reflect dd a dip b origin xyz 这个就是用里建立斜面镜像的。有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，gen zone reflect••… 7、关于两种

fix 是否有区别？

fix z range z -0.1 0.1

fix x range x -0.1 0.1 any x 1.9 2.1 any

fix y range y -0.1 0.1 any y 1.9 2.1 any

和

fix z range z -0.1 0.1

fix x range x -0.1 0.1

fix x range 1.9 2.1

fix y range y -0.1 0.1

fix y range y 1.9 2.1

答：是一样的，any表示并列连接关系，及and的意思

8、initial apply set 的区别

答：initial 初始化命令，初始化计算应力状态，定义材料参数、重力

等；apply边界条件限制命令，如施加边界的力、位移等约束等；set用于设置计算的模式，控制计算过程。

Ini 施加的初始力在计算过程中是会变的，一般体力需要初始化，而apply 施加的力、边界条件不会随着计算变化。set命令作用于整个网格（赋值全局变量）,ini 则初始化指定范围内网格点、域的偏移量。初始应力不仅仅是重力引起的。

9、T ension-p shear-p tension-n shear-n 表示的含义答：在计算的循环里面,每个循环中,每个zone 都依据failure criterion 处于不同的状

态,shear 和tension 分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态。

10、h istory 保存数据，print 显示

在错误的命令流中“原点位移.dat”文件中只有：“Step 1 Y-Displace 而没有保存下数值，怎么回事？我怎么才能获得table中的数据？

输入'print tabel n'能得到table中的数据

（正确）

hist gp ydisp 0,0,0

hist zone syy 0,1,0

hist zone syy 1,1,0

step 3000

hist write 1 file 原点位移.dat

hist write 2 file 中部单元应力.dat hist write 3 file 边缘单元应力.dat hist write 1 table 1

hist write 2 table 2

hist write 3 table 3

hist write 2 3 vs 1 file 单压.dat （错误）

hist gp ydisp 0,0,0

hist zone syy 0,1,0

hist zone syy 1,1,0

hist write 1 file 原点位移.dat hist write 2 file 中部单元应力.dat hist write 3 file 边缘单元应力.dat hist write 1 table 1

hist write 2 table 2

hist write 3 table 3

hist write 2 3 vs 1 file 单压.dat step 3000 11、plot 显示x、y 坐标图hist n 1

hist gp ydisp 0,0,0

hist ax_str

hist gp xdisp 1,1,0

step 1500

plot hist 2 vs -1 ;axial stress vs axial disp.

plot hist 3 vs -1 ; circumferential disp. vs axial disp

答：那 2 是不是代表hist gp ydisp 0,0,0，而-1 指得是hist gp ydisp

0,0,0 为什么不是1 而是负值？

1）在plot hist m vs n的形式里,m代表y轴,n代表x轴（不管m,n的正负）;

2）"-" 表示对其值作"mirror", 如果大家使用AUTOCAD 的话就明白"mirror"的含义了，比如上面的hist 1 gp ydisp的值是0~100那么vs -1 就变成-100~0.以此类推.

3）要这么做主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内。

12、关于water.fis 的命令注释（设置地下水位）

答：（1）water.fis的功能是设置地下水位面，在二维计算中，通过water table 命令即可设置水位线，在三维计算中由于地下水位是面，因此需要作一些插值运算，以便形成水位面。这个FISH 程序包含三个函数parm,aux1 和aux2。parm 用来设定水位的边值，即设置水位的区域，aux1 插值计算等值线的点，aux2 调用aux1 的计算结果，使用

FLAC3D 命令形成水位面。

( 2)首先设置参数值，然后定义一个table.

def parm

nptab1=7

nprof=10