flac3d命令流实例大全

第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 blo

new

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill

gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1 PLOT SURF PAUSE K

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,3

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×

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

实例分析命令：

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: creat

实例分析命令：

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: creat

实例分析命令：

1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转

Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,（重力加速度） 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 图像输出（格式、像素、名称） plot set magf 视图的放大倍数为；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 i

1、模型建好如何检验所建模型是否有问题

注意：第一步建模一定不要建太大的模型，可以建非常小的模型，确保模型的准确性，这样可以尽快开始后边的步骤，等后边计算分析步骤快速做完验证没有大问题，再将准确的模型建好换上，重新计算得到更加准确的结果。 2、如何检验应力边界、位移边界、加载速率等是否准确

Plot add gpfix：看位移边界

Plo add fap：看应力边界，看加载方向是否准确 Plot add vel：看速率，加载速率是否准确 3、如何查看运算结果

Plo con szz outline on：看z向应力云图，并显示边界网格 Plo add disp：看位移

Plo con zdisp outline on magf 5：看z向位移云图，变形扩大5倍 Plo add bo ：添加边界 Plo blo state：看是否屈服

Plo add stensor：查看大主应力方向 Plo con smin：小主应力 Plo con smax：大主应力

Plo con ssi（shear strain increment）：查看剪切应变增量云图 Plo bcont smin：查看单元小主应力 4、如何做切面并查看切面结果

在编辑“Edit“菜