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 blo gro gpnum on命令

5、关于gauss_dev对性质进行高斯正态分布的问题？

答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是±5％。则命令如下：prop friction 40 gauss_dev 2

问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×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、Tension-p shear-p tension-n shear-n表示的含义

答：在计算的循环里面,每个循环中,每个zone 都依据failure criterion 处于不同的状态,shear和tension 分别表示因受剪和受拉而处于塑性状态.n 表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态。

10、history保存数据，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