《FLAC3D命令流(挺实用)》由会员分享，可在线阅读，更多相关《FLAC3D命令流(挺实用)（20页珍藏版）》请在金锄头文库上搜索。
1、1、怎样查看模型？答： plot grid 可以查看网格， plot grid num 可以查看节点号。2、请问在圆柱体四周如何施加约束条件？答：可以用 fix . ran cylinder end1 end2 radius r1 cylinder end1 end2radius r2 not，其中 r2r1，其实就是选择 range 的时候选两个圆柱的差，即得到边界。命令流如下：fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT 的图像数据导出来以便用其他软件绘图？答：用 set lo
2、g on 命令，把数据导出来，转到后用 surfer 或者什么作图软件绘制就行了。excel 里处理一下，然4、用命令建立模型后，如何显示点的坐标？答：使用 plo blo gro gpnum on 命令5、关于 gauss_dev对性质进行高斯正态分布的问题？40 度，答：根据手册上的说明：下面的命令设定一个平均摩擦角为标准方差是 5。则命令如下： prop friction 40 gauss_dev 2问题：请问 gauss_dev 2中的 2 是如何计算的？如果把5改为10，则命令应如何写？405 26、 reflect 问题问： gen zone radbrick &p0 (0,0,0
3、) 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 fillplot surfgen 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 轴正方向按顺时针 （
4、clockwise）方向转向所要确定面的法线方向在 xy 平面上的投影的夹角， 然后再确定 dip，dip 指的是从 xy 平面转向所要确定的平面的角（及 z 轴负方向转过角度）（ 2）命令：reflect normal xyz origin xyz （根据法线和过一点建的面）最常用，至于实际使用过程中有人问道： normal -1 0 0 与 normal 1 0 0 的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表达的是同一个面。reflectdd a dip b origin xyz这个就是用里建立斜面镜像的。有一点注意了， 镜像的时候如果是模型镜像的话， 相当于在建立一个模
5、型，故，在 reflect 的前面要加建模的关键字， gen zone reflect7、关于两种 fix 是否有区别？fix z range z -0.1 0.1fix x range x -0.1 0.1 any x 1.9 2.1 anyfix y range y -0.1 0.1 any y 1.9 2.1 any和fix z range z -0.1 0.1fix x range x -0.1 0.1fix x range 1.9 2.1fix y range y -0.1 0.1fix y range y 1.9 2.1答：是一样的， any 表示并列连接关系，及and 的意思。8
6、、initialapplyset 的区别答：initial初始化命令，初始化计算应力状态，定义材料参数、重力等；apply 边界条件限制命令， 如施加边界的力、位移等约束等； set 用于设置计算的模式，控制计算过程。Ini 施加的初始力在计算过程中是会变的，一般体力需要初始化，而 apply 施加的力、边界条件不会随着计算变化。 set 命令作用于整个网格(赋值全局变量 ),ini 则初始化指定范围内网格点、域的偏移量。初始应力不仅仅是重力引起的。9、Tension-p shear-p tension-n shear-表n示的含义答：在计算的循环里面 ,每个循环中 ,每个 zone 都依据
7、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,0hist zone syy 0,1,0hist zone syy 1,1
8、,0step 3000hist write 1 file 原点位移 .dathist write 2 file 中部单元应力 .dathist write 3 file 边缘单元应力 .dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压 .dat（错误）hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0hist write 1 file原点位移 .dathist write 2 file 中部单元应力
9、.dathist write 3 file 边缘单元应力 .dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压 .datstep 300011、plot 显示 x、y 坐标图hist n 1hist gp ydisp 0,0,0hist ax_strhist gp xdisp 1,1,0step 1500plot hist 2 vs -1 ;axial stress vs axial disp.plot hist 3 vs -1 ; circumferential
10、disp. vs axial disp答：那 2 是不是代表 hist gp ydisp 0,0,0，而-1 指得是 hist gp ydisp 0,0,0 为什么不是 1 而是负值？1) 在 plot hist m vs n 的形式里 ,m 代表 y 轴,n 代表 x 轴(不管 m,n 的正负);2) - 表示对其值作 mirror, 如果大家使用 AUTOCAD 的话就明白 mirror 的含义了 ,比如上面的 hist 1 gp ydisp的值是 0100,那么 vs -1 就变成 -1000.以此类推 .3) 要这么做主要是为了符合岩土工程的习惯需要 ,例如桩载荷试验曲线,沉降曲线等都
11、不是画在第一象限内。12、关于 water.fis的命令注释（设置地下水位）答：（1）water.fis 的功能是设置地下水位面， 在二维计算中，通过 watertable 命令即可设置水位线，在三维计算中由于地下水位是面，因此需要作一些插值运算，以便形成水位面。这个FISH 程序包含三个函数 parm,aux1 和 aux2。parm 用来设定水位的边值，即设置水位的区域， aux1 插值计算等值线的点， aux2 调用 aux1 的计算结果，使用 FLAC3D 命令形成水位面。（2）首先设置参数值，然后定义一个 table.def parmnptab1=7nprof=10endparmta
12、ble 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 语言，这三个变量的含义就不难理解了，只不过
13、这个操作与真正的数据库操作相比，只有两个数据变量 x,y。我用数据库的术语作解释： table 就是产生一个“表” ， xtable,ytable 的作用就是产生一个新表 (Create)以及对表中的数据进行取出(Select)、插入 (Insert)、替换 (Update)操作 .4）FISH 函数 aux1 的作用是在 table 1 的基础上产生出另外个 table， table 的 ID 从 11 到 17，每个 table 有 10 个数据集。 alfa 的计算方法是为适合 (fit) 曲面而采用的坐标点算法。def aux1loop n(1,nptab1) rr=xtable(1,
14、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 2xtable(n+10,k)=rrytable(n+10,k)=0.0end_caseend_loopend_loopendaux1当执行这个 FISH 函数后，生成 (nptab1+10,nprof)组数据，这些数据为 Aux2 调用。注意：本段程序产生的数据是为那个特定的曲面边坡制作的，不是一个 general purpose的 LIB5）FISH 函数 aux2 的作用是使用 WATER TABLE face 命令产生水 面，把 table 1 的 ytable 值作为 z 坐标， aux1 产生的数据作为 x,y 坐标。 def aux2commandset grav 0 0 -10water dens 1000end_com
