FLAC3D命令流(整理版)文档格式.doc

FLAC3D命令流(整理版)文档格式.doc

《FLAC3D命令流(整理版)文档格式.doc》由会员分享，可在线阅读，更多相关《FLAC3D命令流(整理版)文档格式.doc（21页珍藏版）》请在冰豆网上搜索。

下面的命令设定一个平均摩擦角为40度，标准方差是±

5％。

则命令如下：

propfriction40gauss_dev2

问题：

请问gauss_dev2中的2是如何计算的？

如果把±

5％改为±

10％，则命令应如何写？

40×

5％＝2

6、reflect问题

问：

genzoneradbrick&

p0（0,0,0）p1（10,0,0）p2（0,10,0）p3（0,0,10）&

size3,5,5,7&

ratio1,1,1,1.5&

dim142fill

plotsurf

genzonereflectdip0dd90（对xy面做镜像）

genzonereflectdip90dd90（对yz面做镜像）

（1）dd表示y轴正向顺时针到那条射线的夹角，dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线。

首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向

在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．（及z轴负方向转过角度）

（2）命令：

reflectnormalxyzoriginxyz（根据法线和过一点建的面）最常用，至于实际使用过程中有人问道：

normal-100与normal100的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。

reflectddadipboriginxyz这个就是用里建立斜面镜像的。

有一点注意了，镜像的时候如果是模型镜像的话，相当于在建立一个模型，故，在reflect的前面要加建模的关键字，genzonereflect……

7、关于两种fix是否有区别？

fixzrangez-0.10.1

fixxrangex-0.10.1anyx1.92.1any

fixyrangey-0.10.1anyy1.92.1any

和

fixxrangex-0.10.1

fixxrange1.92.1

fixyrangey-0.10.1

fixyrangey1.92.1

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

8、initialapplyset的区别

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

apply边界条件限制命令，如施加边界的力、位移等约束等；

set用于设置计算的模式，控制计算过程。

Ini施加的初始力在计算过程中是会变的，一般体力需要初始化，而apply施加的力、边界条件不会随着计算变化。

set命令作用于整个网格（赋值全局变量）,ini则初始化指定范围内网格点、域的偏移量。

初始应力不仅仅是重力引起的。

9、Tension-pshear-ptension-nshear-n表示的含义

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

10、history保存数据，print显示

在错误的命令流中“原点位移.dat”文件中只有：

“Step1Y-Displace”

而没有保存下数值，怎么回事？

我怎么才能获得table中的数据？

输入'

printtabeln'

能得到table中的数据

（正确）

histgpydisp0,0,0

histzonesyy0,1,0

histzonesyy1,1,0

step3000

histwrite1file原点位移.dat

histwrite2file中部单元应力.dat

histwrite3file边缘单元应力.dat

histwrite1table1

histwrite2table2

histwrite3table3

histwrite23vs1file单压.dat

（错误）

histgpydisp0,0,0

histwrite1file原点位移.dat

histwrite3table3

step3000

11、plot显示x、y坐标图

histn1

histax_str

histgpxdisp1,1,0

step1500

plothist2vs-1;

axialstressvsaxialdisp.

plothist3vs-1;

circumferentialdisp.vsaxialdisp

那2是不是代表histgpydisp0,0,0，而-1指得是histgpydisp0,0,0为什么不是1而是负值？

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

2）"

-"

表示对其值作"

mirror"

如果大家使用AUTOCAD的话就明白"

的含义了,比如上面的hist1gpydisp的值是0~100,那么vs-1就变成-100~0.以此类推.

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

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

（1）water.fis的功能是设置地下水位面，在二维计算中，通过watertable命令即可设置水位线，在三维计算中由于地下水位是面，因此需要作一些插值运算，以便形成水位面。

这个FISH程序包含三个函数parm,aux1和aux2。

parm用来设定水位的边值，即设置水位的区域，aux1插值计算等值线的点，aux2调用aux1的计算结果，使用FLAC3D命令形成水位面。

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

defparm

nptab1=7

nprof=10

end

parm

table1（120）（18.2512.5）（2517）（3521）（5024）（7525）（15025）

产生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的作用是在table1的基础上产生出另外个table，table的ID从11到17，每个table有10个数据集。

alfa的计算方法是为适合（fit）曲面而采用的坐标点算法。

defaux1

loopn（1,nptab1）

rr=xtable（1,n）

zz=ytable（1,n）

loopk（1,nprof）

case_ofk

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）

case1

xtable（n+10,k）=rr

ytable（n+10,k）=-40.0

case2

ytable（n+10,k）=0.0

end_case

end_loop

end_loop

aux1

当执行这个FISH函数后，生成（nptab1+10,nprof）组数据，这些数据为Aux2调用。

注意：

本段程序产生的数据是为那个特定的曲面边坡制作的，不是一个generalpurpose的LIB

5）FISH函数aux2的作用是使用WATERTABLEface命令产生水面，把table1的ytable值作为z坐标，aux1产生的数据作为x,y坐标。

defaux2

command

setgrav00-10

waterdens1000

end_command

loopn（1,nptab1-1）

loopk（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）