FLAC3D学习（SIMWE论坛君之枫整理） 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
结果在“原点位移.dat”文件中只有：“Step 1 Y-Displace” 而没有保存下数值，怎么回事？我怎么才能获得table中的数据？ 修改：
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4 gen zone reflect norm 1,0,0 gen zone reflect norm 0,0,1 model mohr
prop bulk 1.19e10 shear 1.1e10 prop coh 2.72e5 fric 44 ten 2e5 fix x y z range y -.1 .1 fix x y z range y 1.9 2.1 ini yvel 1e-7 range y -.1 .1 ini yvel -1e-7 range y 1.9 2.1 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
save a 'print tabel n' will solve your second qusetion 22.这个plot是什么意思？ hist n 1
hist gp ydisp 0,0,0 hist ax_str hist gp xdisp 1,1,0
36
FLAC3D学习（SIMWE论坛君之枫整理） step 1500
plot hist 2 vs -1 ;axial stress vs axial disp.
plot hist 3 vs -1 ; circumferential disp. vs axial disp plot 显示 hist 历史曲线 曲线2为横坐标，曲线一为纵坐标 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而是负值？
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
1) 在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负);
2) \表示对其值作\如果大家使用AUTOCAD的话就明白\的含义了,比如上面的 hist 1 gp ydisp的值是0~100,那么vs -1就变成-100~0.以此类推.
3) 为什么要这么做呢?主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内.
23.关于water.fis的命令注释（设置地下水位）
1）water.fis的功能是设置地下水位面，在二维计算中，通过water table命令即可设置水位线，在三维计算中由于地下水位是面，因此需要作一些插值运算，以便形成水位面。这个FISH程序包含三个函数parm,aux1和aux2。parm用来设定水位的边值，即设置水位的区域，aux1插值计算等值线的点，aux2调用aux1的计算结果，使用FLAC3D命令形成水位面。 2）首先设置参数值，然后定义一个table. def parm nptab1=7
37
FLAC3D学习（SIMWE论坛君之枫整理） nprof=10 end parm
table 1 (12 0) (18.25 12.5) (25 17) (35 21) (50 24) (75 25) (150 25) 产生table1，包含7组数据点。
3) 在解释的两个FISH函数之前，对table,xtable,ytable作一些总的说明。water.fis程序反复使用了这三个变量，理解了这三个变量的含义，water.fis这段程序就解决了大半问题了。
Itasca的所有软件都使用了相同的数据结构储存数据，称之为链表结构，类似于关系型数据库中的数据链接方法。如果有数据结构的知识，或掌握关系型数据库的SQL语言，这三个变量的含义就不难理解了，只不过这个操作与真正的数据库操作相比，只有两个数据变量x,y。我用数据库的术语作解释：table就是产生一个“表”，xtable,ytable的作用就是产生一个新表(Create)以及对表中的数据进行取出(Select)、插入(Insert)、替换(Update)操作.
4）FISH函数aux1的作用是在table 1的基础上产生出另外个table，table的ID从11到17，每个table有10个数据集。alfa的计算方法是为适合(fit)曲面而采用的坐标点算法。 def aux1 loop n(1,nptab1) rr=xtable(1,n) zz=ytable(1,n) loop k(1,nprof) case_of k
alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2) xtable(n+10,k)=rr*cos(alfa) ytable(n+10,k)=rr*sin(alfa) case 1
xtable(n+10,k)=rr ytable(n+10,k)=-40.0 case 2
xtable(n+10,k)=rr ytable(n+10,k)=0.0 end_case end_loop end_loop end aux1
当执行这个FISH函数后，生成(nptab1+10,nprof)组数据，这些数据为Aux2调用。注意：本段程序产生的数据是为那个特定的曲面边坡制作的，不是一个general purpose的LIB
38
FLAC3D学习（SIMWE论坛君之枫整理）
5）FISH函数aux2的作用是使用WATER TABLE face命令产生水
面，把table 1的ytable值作为z坐标，aux1产生的数据作为x,y坐标。 def aux2 command
set grav 0 0 -10 water dens 1000 end_command loop n(1,nptab1-1) loop k (1,nprof-1) xx1=xtable(n+10,k) yy1=ytable(n+10,k) zz1=ytable(1,n) xx2=xtable(n+10,k+1) yy2=ytable(n+10,k+1) zz2=zz1
xx3=xtable(n+11,k+1) yy3=ytable(n+11,k+1) zz3=ytable(1,n+1) xx4=xtable(n+11,k) yy4=ytable(n+11,k) zz4=zz3 command
WATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3 WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4 end_command end_loop end_loop end aux2
24.产生结构元的几何形状(Geometry Creation)
FLAC3D有6种类型的结构单元，它们分别是梁单元beamSELs,锚杆单元cableSELs,桩单元pileSELs,壳单元shellSELs,地质网格单元geogridSELs和衬砌单元linerSELs.每一种类型的结构由相应分量对象的集合来表示。例如锚杆结构(cable)由一系列的cableSELs组成，而衬砌结构(liner)则由一系列的linerSELs组成。每种结构元都由两个不同的辨识号码来表示，以便与其它结构元相区别。
用Seltype命令和相应的关键字来定义结构单元:Seltype
beam/cable/geogrid/liner/pile/shell,每一种支护类型的性质用SEL 命令的property关
39
FLAC3D学习（SIMWE论坛君之枫整理） 键字定义.
History SEL keywords
ID：用来标识一个结构对象，比如第一根锚杆的ID＝1，第二根锚杆的ID=2，它是对整个结构进行编号；
CID:用来表示一个结构对象内部的分量，比如一根锚杆划分为3段，则这根锚杆的CID分别是1,2,3,用来表示这根锚杆的不同部分。
ID是为结构进行编号，CID是为结构的分量进行编号。比如下面的例子： sel pile id=1 begin=(10.0, 1.0, 0.0) end=(10.0, 1.0, -10.0) nseg=5 sel pile id=2 begin=(10.0, 3.0, 0.0) end=(10.0, 3.0, -10.0) nseg=5
在这个例子中，设置了2根桩，第一根桩以id=1来表示，第二根桩以id=2来表示；第一根桩分为5段，则cid=1~5，第二根桩也分为5段，则cid=6~10。 可以用plot sel geom 命令来查看设置是否正确。 问题: 我设了三根桩，
sel pile id=1 begin=(10.0, 1.0, 0.0) end=(10.0, 1.0, -3.0) nseg=3 sel pile id=2 begin=(10.0, 3.0, 0.0) end=(10.0, 3.0, -3.0) nseg=3 sel pile id=3 begin=(10.0, 5.0, 0.0) end=(10.0, 5.0, -3.0) nseg=3 图中红色的数字1-3-4-2代表什么？是node？还是link？
生成node的时候，如果周围存在zone，会自动生成link，而且link和node位置总是重合的
25. 基于ansys复杂地质体flac模型自动生成
首先将复杂的三维地质模型建于ansys，然后导出结点坐标信息和单元信息，最后运用此程序就可以生成复杂地质体的flac命令，用flac程序运行此命令就可以生成复杂地质体的flac模型。
复杂的三维地质模型一般是在cad中建立的,然后导入ansys中，通过此它建立复杂地形的flac模型。用autocad建好模型后，输出为.sat文件，就可以导入到ansys，不用软件转换。
26. FLAC结果数据的导出 pl set back white ;设置背景
set plot bitmap size (800,510) ;设置图片大小
plot set caption size 37 ; 设置图例数字大小 (default=35 范围10-50) mainwin size 1 1 position 0 0.8 ; command window
plot set window size 0.8 0.8 position 0.1 0.0 ; graphics window
40