FLAC3D各种命令笔记.doc

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 剪胀角
FISH语言：
Degrade;∏/180
Pi;∏
Ngp;节点总数
Nzone;单元体总数
单元节点遍历：
P_z = zone_head
Loop while p_z # null
………
P_z = z_next(p_z)
End loop
Read(ar，n)；将数组ar的前n个记录写到文件

print zone state 显示塑性区域
编写fish函数时，尽量采用浮点型，例如1.0而不是1
Z_sig1 z_sig2 z_sig3 fish中的三个主应力，由大到小
Print gp poison rang id n 输出节点n的坐标
Xtable(n,s)对ID号为n的表的第s行，x列进行赋值，table表只有x,y两列
Pl table 1 line/mark/both 关键词描述曲线的类型，线、标记、二者综合
File_name=’add_step’+string(n)+’.sav’
Gp_near(x,y,z)得到靠近坐标x,y,z的节点地址
Gp_zdis(p_gp)地址为p_gp的节点的z向变形
String(n)将变量n转换成字符串格式
Dim 看基本网格里面的标准
attach face 常用命令，建模型之后连接面
sxx=syy=(λ/1-λ)szz szz= γH
Sel fix x y z (局部坐标系位移速度) / xrot yrot zrot(局部坐标系转动位移)
Ini xdisp ydisp zdisp(全局坐标系位移) / xpos(节点全局坐标) / xrdisp(全局坐标系旋转位移)
Xv(全局坐标系位移速度) / xrvel (全局坐标系旋转速度)
如果边界是FIX控制，则自动平衡，不用apply 控制应力。
对地下深部开挖的工作，一般忽略顶部到底部的地心引力的差别，set grav被省略
全局坐标系 x y z 局部坐标系 d（倾向） s（走向） n（法向）
 plo fap ran group tunnel ;显示施加在组上的体力
apply remove 施加新的边界条件前，用以移除旧条件
Table  n  x1 y1  <x2 y2>  <x3 y3>
             建成一个点的坐标的表, 以供 FLAC 引用，n为是表号。
材料参数附加关键字：
Gauss_dev s 以输入数值为平均数，s为偏差数，高斯正态分布计算
Gradient gx gy gz 渐变：值=输入值+gx*x+gt*y+gz*z
Uniform_dev s 以输入数值为平均数，s为偏差数，随机正态分布计算
his unbal set mech force 50 检测最大不平衡力到50
M n rang group x x1 x2 开挖
Plot grid 显示网格
Ini z add -6 rang group pile *pile组整体在z方向平移-6个单位
建立网格：实体多少米p0 p1 p2…就设多少米，单元数，在百个左右，各方向与实体成比例
Brick 六面块体 redbrick 放射状 cylinder 柱体网格 radcylinder 渐变放射
Radtunnel 六面体隧道渐变放射 cshell 柱形壳体网格
Gen zon redbrick p0…p8 size n1 n2 n3 n4 ratio r1 r2 r3 r4
Model mohr 赋予本构模型
Prop bu n shear n coh n fric n 赋予材料性质 弹性只赋值弹性、剪切模量
Fix 固定边界
Ini dens n 初始化密度
App 施加力
Reflect 镜像命令 gen zon reflect norm x y z ori x y z
Grad 渐变
ini szz -30e3 grad 0 0 20e3 ran z 0 1
应力ZZ方向的分量按（-30e3+0*x+0*y+20e3*z ）渐变
plot blo gro 显示分组 plot block group range group dam 只显示dam组单元
plotitems_block group_modify_勾选grid point id’s_ok 显示节点ID编号
Hist gp 会自动监测距离输入坐标最近的节点
Hist zone会自动监测距离输入坐标最近的单元
Plot hist 7 vs 8
Set Hist_rap n 监测间隔，越小，曲线越圆滑
hist Write 7 vs 8 file 6-3hist.txt 输出对比信息文件
print 输出结果 print gp dis id rang 517 any 输出节点517的变形值
记录文件：
log文件是FLAC输出文件的格式之一
set log on 建立一个名为flac3d.log的文件，记录窗口中所以信息，直到set log off
继续下一个log文件，
set logfile 6-1.log
set log on
print gp zdisp rang id 2179 a id 2189
set log off
plot con zd / zv 显示z方向位移/速度 plot sxx 显示应力等值阴影图（没有contour）
plot con zd ou on magf 20 附带网格线，变形量放大20
plot sk 显示网格线
plot sk dis scale 0.07 网格上附加矢量图，scale默认0.05越大，箭头越长
后处理：
plot set rot 20 0 30  ;输出前 设置视图角度
plot sk dis scale 矢量图显示

剖面显示：
plot set plane ori 0 1.5 0 norm 0 1 0 设定定一个垂直于Y轴的剖面，ori为切面上一点坐标，norm为切面的法向（即从0 1.5 0 引到0 1 0 连线为法线）
pl add con disp plane behind shade on 剖面后面的区域一样显示
plot con zd plane 在设定的剖面上显示z方向位移
plot add ske 增加网格线 add命令意思是 在绘图区域增加
plot add dis plane 增加位移矢量
plot add axe 增加坐标系
动画制作：
Plot set rot 20 0 30 设置视角
Set movie avi step 1 file 6-5.avi 格式设置
Movie star 动画开关
Solve
求解完成后，自动生成文件。
初始地应力场的生成：弹性求解法，分阶段的弹塑性求解法
存在了null模型，不能采用solve elastic的求解方法获得初始应力，所以采用分阶段的弹塑性求解方法。先将Mohr模型的凝聚力。俏和伉拉强度护赋俏为无穷人进行求解，保证在重力作用卜单儿不至十发生屈服，然后冉将Mohr模型参数赋值为真实值，再进行求解。
接触面：分界面应附在小面区域，与密度大的区域。
胶合模型：不允许张开和滑移
库伦剪切模型：
抗拉粘结模型：粘结型和滑移型。Bslipe off 不滑动（默认） on 滑动
接触面几何模型建立：
移来移去法：首一先将需要建立接触而的两个网格分开建立，然后在一个网格的指定位置建立接触而，然后把另外一部分网格移到特定的位置，这种方法称为移来移去法。要点就是把需要建立接触而的网格分开建立，然后在相应的网格上建立接触而，再把没有建立接触而的网格移动过来，这样做的目的就是使同一坐标点有不同的节点号，从而可以模拟结构物的滑移、分离等情祝。
Gen zon radcyl p0---p11 size a（x） b（y） c（弧度上的单元数）d（放射方向的网格） ratio 1 1 1 1.15
interface 1 face（在指定范围内所有网格的面上，均创建分界面单元） range cyl end1 0 0 0 end2 0 0 -5.1 radius .31 & cyl end1 0 0 0 end2 0 0 -5.1 radius .29 not
radius（半径），后一个半径要小于前一个，即两个圆柱的差，得到边界。命令指的整个圆柱再抠掉后面的
End 圆柱体断面中心点的坐标,两个点确定一条线，再根据半径确定柱体。
Not 就是否认前面的意思
Pl sel geom. 显示创建的sel单元
Ini z add -6 rang group pile *pile组整体在z方向平移-6个单位
导来导去法：1、保存完整模型保存为1.Sav
2、保留需要建接触面的分组，删掉其他所有单元，用到命令：del rang group 2 not，在分组上建立接触面，保存为2.sav；
3、rest 1.sav & del rang group 2 & *导出删除后的模型 expgrid 1.flac3d；
4、rest 2.sav & impgrid 1.flac3d
Interface 1 prop kn(法向刚度) ks(切向刚度) c f ks kn 由实验确定，一般取接触面两侧较小的E的10倍 *注意语法

切割模型法：
此方法是在己建好的模型上，把group 1与周围组分开，然后在group 1 上建立接触面。
Gen separate group 1 *分离出group 1建面
Int 1 wrap 1 2 *int=interface 在group 1 group 2之间建立接触单元。
结构单元：
ID 结构单元编号 CID结构构件编号 node节点 link 结构连接的编号（结构单元间、结构与网格间）
Sel beam beg end nseg
Sel beam prop
Sel pile id=1 cid=i beg 0 0 0 end 0 0 10 nseg=4 *建立桩单元，设立结构单元ID和结构构件ID，结构构件4个，三个节点
Sel node id= *设置节点ID
Pl sel geo id on nod on scale 0.04
Sel shellsel id=I rang cyl end1 end2 radius not *在半径处建一个无厚度的面
Sel node ini zpos add -25 *壳单元的位移，节点初始位置减少25
Delete 删除实体单元 sel delete 删除结构单元
plot add sel geom node = off shrinkfac=0 magfac = 10 red red  显示结构单元 shrinkfac是指结构单元的尺寸，范围是0-1，默认值为0.1
桩单元模拟锚杆支护特性：sel pile prop rockbolt on 可以激活锚杆特性，能够计算支护四周的约束应力、锚杆与周围单元之间的应变软化特性以及锚杆拉断、拉裂等现象。
Print sel node ?
Plot sel node ? *后处理
其中fob表示不平衡力的大小，local表示结构节点的局部坐标系
锚杆：三种非全长锚固模拟方法
1、通过删除锚索端头，即Cable端部的node和zone之间自动建立的连接，然后在他们之间建立刚性链接来模拟托盘。锚杆(锚索)自由段和锚固段通过设置不同的锚固剂参数来模拟，预紧力加在锚杆自由段。
sel del link rang id 1 *删除自动链接
sel link id=100 1 target zone *在cable单元与zone之间设置新链接 target 目标
sel link attach xdir=rigid ydir=rigid zdir=rigid & *六个自由度 rigid 刚性的
xrdir=rigid yrdir=rigid zrdir=rigid rang id=100 *设置链接为刚性
sel node fix x y z xr yr zr ran id=1 *固定六个方向的自由度（此句为桩pile）
sel cable id=1 pretension 60e3 rang cid 1 10 *自由段施加预紧力
2、将锚杆(锚索)的端头、自由段、锚固段赋不同的属性来模拟非全长锚固预应力锚杆〔锚索)，端头的锚固参数设为极大来模拟托盘，这样在锚杆〔锚索)受力时，端头将不会滑动，相当于托盘的作用。预紧力加在锚杆(锚索)自由段。
3、借助别的结构单元模拟托盘