完整版)FLAC3D命令流(整理版

1、1、怎样查看模型？答： plot grid 可以查看网格， plot grid num 可以查看节点号。2、请问在圆柱体四周如何施加约束条件？答：可以用 fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius 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 的图像数据导出来以便用其他软件绘图？ 答：用 se

2、t 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，则命令应如何写？ 405 26、reflect 问题问： gen zone radbr

3、ick &p0 （0,0,0） 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 指

4、 的是从 y 轴正方向按顺时针 （ clockwise）方向转向所要确定面的法线 方向在 xy 平面上的投影的夹角， 然后再确定 dip ，dip 指的是从 xy 平面转 向所要确定的平面的角 （及 z 轴负方向转过角度）（ 2）命令：reflect normal xyz origin xyz （根据法线和过一点建的面） 最常用，至于实际使用过程中有人问道： normal -1 0 0 与 normal 1 0 0 的区别，我试过两者的效果是一样的，没区别，虽然方向不同，但表 达的是同一个面。reflect dd a dip b origin xyz 这个就是用里建立斜面镜像的。 有一点注意了

5、， 镜像的时候如果是模型镜像的话， 相当于在建立一个 模型，故，在 reflect 的前面要加建模的关键字， gen zone reflect 7、关于两种 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答：

6、是一样的， any 表示并列连接关系，及 and的意思8、initial apply set 的区别答： initial 初始化命令，初始化计算应力状态，定义材料参数、重力 等；apply 边界条件限制命令， 如施加边界的力、位移等约束等； set 用 于设置计算的模式，控制计算过程。Ini 施加的初始力在计算过程中是会变的，一般体力需要初始化，而 apply 施加的力、边界条件不会随着计算变化。 set命令作用于整个网 格 （赋值全局变量 ）,ini 则初始化指定范围内网格点、域的偏移量。初 始应力不仅仅是重力引起的。9、Tension-p shear-p tension-n shear-n

7、 表示的含义 答：在计算的循环里面 ,每个循环中 ,每个 zone 都依据 failure criterion 处于不同的状态 ,shear 和 tension 分别表示因受剪和受拉而处于塑性 状态 .n 表示 now,p 表示 previous,即分别表示在现在和以前的循环当 中处于塑性状态。10、history 保存数据， print 显示 在错误的命令流中“原点位移 .dat”文件中只有：“ Step 1 Y-Displace” 而没有保存下数值，怎么回事？我怎么才能获得 table 中的数据？ 输入print tabel n能得到 table 中的数据（正确）hist gp ydisp

8、 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0step 3000hist write 1 file 原点位移 .dathist write 2 file 中部单元应力 .dat hist write 3 file 边缘单元应力 .dat hist 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 w

9、rite 1 file 原点位移 .dat hist write 2 file 中部单元应力 .dat hist write 3 file 边缘单元应力 .dat hist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压 .dat step 3000 11、plot 显示 x、 y 坐标图hist n 1hist gp ydisp 0,0,0hist ax_strhist gp xdisp 1,1,0step 1500plot hist 2 vs -1 ;axial stres

10、s 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 而是负值？1）在 plot hist m vs n 的形式里 ,m 代表 y 轴,n 代表 x 轴 （不管 m,n 的正 负）;2）-表示对其值作 mirror, 如果大家使用 AUTOCAD 的话就明白 mirror的含义了,比如上面的 hist 1 gp ydisp的值是 0100,那么vs -1 就变成-1000.

11、以此类推 .3）要这么做主要是为了符合岩土工程的习惯需要 ,例如桩载荷试验曲 线,沉降曲线等都不是画在第一象限内。12、关于 water.fis 的命令注释（设置地下水位）答：（1）water.fis 的功能是设置地下水位面， 在二维计算中， 通过 water table 命令即可设置水位线，在三维计算中由于地下水位是面，因此 需要作一些插值运算，以便形成水位面。这个 FISH 程序包含三个函 数 parm,aux1 和 aux2。 parm 用来设定水位的边值，即设置水位的区 域， aux1 插值计算等值线的点， aux2 调用 aux1 的计算结果，使用 FLAC3D 命令形成水位面。(

12、2)首先设置参数值，然后定义一个 table.def parmnptab1=7nprof=10endparm 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 的所有软件都使用了相同的数据结构储存数据， 称之为链表结 构，类似于关系型

13、数据库中的数据链接方法。 如果有数据结构的知识， 或掌握关系型数据库的 SQL 语言，这三个变量的含义就不难理解了， 只不过这个操作与真正的数据库操作相比，只有两个数据变量x,y。我用数据库的术语作解释： table 就是产生一个“表” ， xtable,ytable 的作用就是产生一个新表 (Create)以及对表中的数据进行取出 (Select)、插入 (Insert)、替换 (Update)操作 .4) FISH 函数 aux1的作用是在 table 1的基础上产生出另外个 table， table的 ID 从11到17，每个 table有 10个数据集。 alfa的计算方法 是为适合

14、(fit) 曲面而采用的坐标点算法。def aux1loop 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 2xtable(n+10,k)=rr ytable(n+10,k)=0.0 end_caseend_loo

15、p end_loop end aux1 当执行这个 FISH 函数后，生成 (nptab1+10,nprof)组数据，这些数据为 Aux2 调用。注意：本段程序产生的数据是为那个特定的曲面边坡制 作的，不是一个 general purpose的 LIB5) FISH函数 aux2的作用是使用 WATER TABLE face命令产生水 面， 把 table 1 的 ytable 值作为 z 坐标， aux1 产生的数据作为 x,y 坐标。 def aux2 commandset grav 0 0 -10water dens 1000 end_command loop n(1,nptab1-1)

16、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=zz1xx3=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=zz3commandWATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3WATER TABLE face

17、xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4 end_commandend_loop end_loop end aux213、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 windowplot set window size 0.8 0.8 position 0.1 0

18、.0 ; graphics windowplot set center 100 0 40根据情况更改14、一段 fish 命令的解释 def aux1 定义函数 loop n(1,nptab1) 开始一个循环，从 1 循环到变量 nptab1 rr=xtable(1,n) 从 table 1 的 x 列中读第 n 个数付给变量 rr zz=ytable(1,n) 从 table 1 的 y 列中读第 n 个数付给变量 zz loop k(1,nprof) 开始第二个循环，从 1 循环到变量 nprof case_of k 下面是一个选择结构，利用变量 k 判断 alfa=0.5*(pi+.1)

19、*float(k-2)/float(nprof-2) xtable(n+10,k)=rr*cos(alfa) ytable(n+10,k)=rr*sin(alfa) case 1 当 k 1 的时候执行下面语句 xtable(n+10,k)=rr table n+10 的 x 列第 k 个值等于变量 rr ytable(n+10,k)=-40.0 table n+10的 y 列第 k 个值等于变量 -40 case 2 当 k 2 的时候执行下面语句 xtable(n+10,k)=rr table n+10 的 x 列第 k 个值等于变量 rr ytable(n+10,k)=0.0 table

20、 n+10的 y 列第 k 个值等于变量 0 end_case 结束选择结构 end_loop 结束第二个循环 end_loop 结束第一个循环 end 函数定义结束 aux1 执行什么定义的函数15、有关命令的解释( 1)apply nvel 0.1 plane dip 60 dd 270 range left_bounapply 施加边界条件， plane 法向速度为 0.1， plane 面的位置：倾角 60(在 xoy 平面从-Z 轴方向逆时针旋转) ， dd 走向 270度，从 +y 方 向绕 z 轴顺时针旋转；( 2) ini szz -5e4 grad 0 0 -1e4ini 初

21、始边界条件， szz=-5e4+x*0+y*0+-z*-1e4, 注： -5e4 是 range范 围起点坐标的 SZZ 值， x,y,z所施加应力网格点的坐标。16、渗流示例; 参数部分bulk modulus, K 390 MPa 体积模量（土质比较硬） shear modulus, G 280 MPa 剪切模量 soil dry density, d 1200 kg/m3 土的密度 water density, w 1000 kg/m3 水的密度 wall density, wal 1500 kg/m3 档土墙的密度 permeability, k 10-12 m2/Pa-s 渗透系数

22、porosity, n 0.3 孔隙率fluid bulk modulus, Kf 2.0 GPa 流体的体积模量; 以下是命令流; Excavation in a saturated soil config fluid ; 设置流体算法; - geometrical model - 建立模型gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1 group soil ;设置土组group excavate range x 0 4 y 0 4 z 0 5设 ; 置开挖范围组group wal1 range x 4

23、 5 y 0 5 z 0 7 ; 设置档土墙 1 组group wal2 range x 0 4 y 4 5 z 0 7 设; 置档土墙 2 组group wall range group wal1 any group wal2 any 档; 土墙 1 和 2 合并为 一个档土墙组; - fluid flow model - 流体模型model fl_iso ; 各项同性流体prop perm 1e-12 poro 0.3 设; 置渗透系数和？比 or ？率ini fdensity 1e3 ; 初始化水密度ini fmod 2.0e9 ftens -1e-3 ; 初始化水的体积模量 及 ？mo

24、del fl_null range group excavate ; 开挖部分没有水流model fl_null range group wall ; 档土墙 不渗水ini pp 0 grad 0 0 1e4 ;初始化孔隙水压力fix pp range z -0.1 0.1 ;约束 应该是整个模型的底部吧？fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ; 开挖部分也给约束住了，而 其他的地方不管; - mechanical model - 力学模型model elas ; 弹塑性prop bul 3.9e6 shea 2.8e6 设; 置体积模量和

25、剪切模量，这两个必须， 有一次我只设置了体积模量，就错的一塌糊涂model null range group excavate 力; 学模型的开挖ini density 1.2e3 ;初始化密度 ，应该是土的密度ini density 1.5e3 range group wall ;初始化墙（档土墙）的密度，fix x range x -.1 .1 ; 施加约束，据说差分法约束 0 的时候，必须取一 个-0.100.1,范围fix x range x 11.9 12.1 ;施加约束，或者叫做边界条件fix y range y -.1 .1 ;施加约束fix y range y 11.9 12.

26、1fix z range z 11.9 12.1; initial total stresses 出世总应力ini szz 0 grad 0 0 -1.5e4 ;z方向的应力大一些ini sxx 0 grad 0 0 -1.2e4 ;水平方向的，小一些ini syy 0 grad 0 0 -1.2e4 ;apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstre

27、ss -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1; - setting - 设置set gravity 0 0 10 ;设置重力加速度; - initial state - ; 初始化状态solve force 1 ; check initial equilibrium 求解平衡; - histories - 记录set hist_rep 40 记; 录步为 40hist fltime ; 记录时间hist gp pp 0 0 6 ;记录某点的孔隙水压力hist gp xdis 4 0 0 ;水平位移hist gp xdis 4 0 2 ;记录水平位移

28、hist gp xdis 4 2 0hist gp zdis 0 0 5 ;记录竖向位移hist gp zdis 2 0 5hist gp zdis 4 0 5hist gp zdis 2 2 5hist gp zdis 4 2 5hist gp zdis 4 4 5hist gp zdis 10 0 1hist gp zdis 10 0 2 ; - excavation - 开挖set fluid off ; 关闭水; apply pore pressure at walls 在墙上施加孔隙水压力apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y

29、 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1 solve ;force 1 ;求解 save exc1.sav ; 保存; - drainage - 排水 apply remove nstress 撤; 掉刚才的压力 def relaxsetup ;定义一个释放 函数的参数 step0 = step end relaxsetup ;调用 def relax ;

30、定义一个释放函数 rstep = step - step0 if rstep gp0_dispgp0_disp=gp_dispendifp_gp = gp_next(p_gp)endloop;找最大值的坐标p_gp=gp_headss=0loop while p_gp # nullgp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp = sqrt(gp_disp

31、)if gp_disp=gp0_dispss=ss+1xtable(1,ss)=gp_xpos(p_gp)ytable(1,ss)=gp_ypos(p_gp)endifp_gp = gp_next(p_gp)endloopendget_gp_maxdisp18、 答：如何显示计算后模型的变形状态？ pl sk magfac *（ 一个放大的数字 ）19、 答：如何查看剖面上任一点的位移值print gp disp range x() y( ) z( )或 print gp disp range idplot gp disp range x() y( ) z( )或 plot gp disp r

32、ange id20、group 命令划分斜面答： group groupname colour range plane above dd * dip * origin * plane below dd * dip * origin *dd 和 dip 就是两个斜面的倾向和倾角 ,origin 是这个斜面通过的某个 点,这个命令就可以 group 两个斜面间的区域21、ansys 导入 flac 的时候，如何划分一下 group 答：在 ansys 中对不同 group 用不同的材料号（ mat）定义，再用 ansys2flac之类的程序转换成 flac 的.dat 文件就可， group 自动

33、生成。22、fix x y z range x -0.1 0.1 afix x range x -0.1 0.1b请问 a式与 b 式有何区别？ 答： a 表示 在 x=0 这个平面 是固定的 b 表示在 x=0 这个平面 只 x 方向不能发生运动23、如何知道每个单元的 ID ？ 答：用鼠标双击单元的表面，可以知道单元的 ID 和坐标24、如何进行切片？答：plo set plane ori （点坐标） norm （法向矢量 ）plo con sz plane （显示 z 方向应力的切片 ） 这两行命令，用来显示指定面上的位移变化云图（把第二行改改，还 可以显示应力云图，等其他） ；plot

34、 plot plot plotcon syycon szz outline on con szz effective ony 向应力云图z 向应力云图加显示网格z 向有效应力con szz inter 1e4 云图显示增量（后面的值越小，图上颜色条带就越细，）plot con szz min -10e4 显示最小值， 图上小于这个数的值都用同一 色表示plot con szz max 10e4 显示最大值，图上大于这个数的值都用同一 色表示plot fip outside on显示力的箭头现在看边界，如果是不动面，用 fix x range 如果边界带速率，或位 移用 ini set grav

35、 0 0， -9.81ini dens 1000设置重力就一定要设置密度， 不考虑重力就不要设计密度。 两者共存 的，设计程序的时候一定要记的。set mech force 50solve用 SET 命令设定了最大不平衡力的上限值为 50N ，即：计算的最大 不平衡力一旦小于此值， 我们便人为模型达到了平衡状态。 求解过程 终止。设置这个目的是在一定精度要求下，加快计算（值越大，算的 量越少，花的时间也越少） 。25、如何在程序中进行暂停，并可恢复计算？答：在命令中加入 pause命令，用 continue 进行继续。 在我们分步求解中想得到某一个过程中的结果， 不用等到全求完， 还 可以在分

36、布求解错误的时候就进行改正，而不是等到结果出来。26、如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步27、Fish 是什么东西？ Fish 是否一定要学？答：是 FLAC3D 的内置语言，可以用来进行参数化模型、完成命令 本身不能进行的功能。 Fish可以不用学，需要的时候查 Mannual 获得 需要的变量就可以了。28、如何调用一些可选模块？ 答： config dyn （fluid, creep, cppudm） 29、range 的命令： 取圆球面的范围，就是同一个中心，不同半径之间夹的那一部分， 用在给球面施加力，约束时用该命令，命令：range annulus

37、 center xy（z 中心坐标） radius r1 r2（两半径），想要表达 的是球面时， r1r2取值相差小点即可， 如 radius 2.9 3.1 这时范围就是 半径为 3 的面。圆柱的环面施加力或约束时也用到该命令。range cylinder end1 xyz end2 xyz radius r(end 是表示柱体断面的中心 点坐标 )range plane above 表示的是 z 大于 0 的上部结构range plane blow表示的是 z 小于 0 的下部结构range dd a dip b origin xyz表示的是走向角度 a，倾向角度 b的面dd 表示，当面与

38、 xoy 面相交时的垂线与 y 轴正向的夹角，顺为正， 说白话，就是面的指向方向面是朝那个方向的意思。dip 面与 xoy 面的夹角，与 z 轴负向为正， origin 表示面经过 xyz 点； range normal xyz 法线方向为 xyz 的面range not 没试用过(范围外的东西，想想用到圆上取范围非常好) 。30、如何调用 FISH 文件？答：新手调用 FISH文件出错的原因：是：文件名为 fishcall.fis.txt( 因 为在 记事本 里这 个 txt 的 后缀是默 认 的)。修 改方 法是：将 fishcall.fis.txt 里面的内容全部复制到新的记事本里，再按

39、保存，保存 时候将保存类型改为所有文件 (*,*), 这时你便可在文件名栏内输入 fishcall.fis ，类型也就正确了。31、问 a 式与 b 式有何区别？答： fix x y z range x -0.1 0.1 afix x range x -0.1 0.1b a 表示 在 x=0 这个平面 是固定的； b 表示在 x=0 这个平面 只 x 方向不能发生运动。32、如何查看各个时段不平衡力的具体数值？ 答：采用 his 来记录计算，包括位移应力等 命令 his unbal； his gp(zone) zdis range (0 0 0) 或者 id=？； 导出数据命令 his wri

40、te n vs m begin 时步 end 时步 + file filename.his n 表示纪录的 id m 表示时步； 要导出不平衡力的具体数值 his unbal step 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his 使用上述命令就可以查看各个时步下的不平衡力的具体数值。33、如何显示变形轮廓线的命令？ plo ske magf 10 其中 10 为放大系数34、怎样查看模型的矢量图？ pl grid vel (disp) red plo blo gro 可以看到不同的 group 的颜色分布(看

41、模型的样子) ； plo gpfix red sk看模型的边界情况plo fap red sk看模型的体力分布plo con dis (xdis, ydis, zdis) 看模型的云图：位移plo con sz (sy, sx, sxy, syz, sxz) 看模型的云图：应力 plo dis (xdis, ydis, zdis)看模型的矢量图plo info 看模型有多少单元、节点 plot grid可以查看网格， plot grid num 可以查看节点号。plot fix查看边界条件plot grid mark 查看预先留设的开挖边界pl syy (sxx xdis ydis sig1.

42、) fill ；查看模型的云图。 注意， 带 fill 关键字 就是查看云图，不带 fill 关键字就是查看等直线图。plo add ske 剖面的沉降云图上增加网格 plo add dis plane 剖面上的矢量 plo add axe 剖面上坐标值35、怎样查看 INTERFACE ？ 答：pl grid iface id red print iface id ; id 是 interface 的编号36、FLAC 怎么调整视图？答： ctrl+z 或 shift+z ,然后点击要放大部位的左上角和右下角就可以 放大某个部位，如果要缩小，双击屏幕。37、判断模型中网格点是一个还是两个的方

43、法。答：显示模型网格点的 ID ，如果有两个 ID 号，说明是两个。节理面 上的网格点应该有两个 ID 号。38、flac 中影响运行时间的因素 答：(1)单元材料，结构面之间的刚度对比很大； (2)划分单元尺寸对比很大， 尽量避免长细比大于 5：1 的细长单元时时显示云图39、简单命令解释：答： 1.plot add cont disp out on2. plot block group gpnum on3. plot inter red ske4. del range group 1 not5. del range group 2显示节点编号 看接触面的效果 . 不删除 group 信息

44、删除 group 信息interface 1 建立接触面 清除水位线 设置背景 设置图片大小 设置图例数字大小 设置图片中心表示在 x=0 这个平面是固定的 表示这个平面只 x 方向不能发生运动6.ini z add * range group * （ name） 7.interface 1 face range x（* ，*）y（* ，*）z（*，*）或者 face range plane orig（x y z 坐标） norm（法线方向）8. water table clear9. pl set back white10. set plot bitmap size （800,510）11.

45、 plot set caption size 3712. plot set center 100 0 4013. fix x y z range x -0.1 0.114. fix x range x -0.1 0.1 40、在固定面上施加位移时， fix 与 apply 不能同时使用，应该用 ini 和 fix, 才能施加上。41、在进行三轴压缩时候，如何施加围压？答：使用 apply sxx =-1e5 szz=-1e5 range cyl end1 0,0,0 end2 0,2,0 radius 1，其中 X 和 Z 轴是在我们眼睛正对的面。42、在进行三轴压缩时候，如何对轴向力进行位移

46、控制？答：使用 ini yvel 1e-7 range y -0.1 0.1 ini yvel 1e-7 range y 1.9 2.1 即 可，数值可以变化。43、在 FLAC3D 计算中，如何记录程序中的值，如果不想时时记录 该怎么弄？答：用 hist wtite n1 vs n2 file.txt ， 想间隔记录可用 set hist_rep n44、想记录接触面上的相关力学参数该怎么办？答： hist inter i nstress /ndis/sstress/sdisp 45、接触面设置参数包括哪几个？ 答： coh dial fric kn ks ten.46、如何输出单元格的数据

47、？ 答； expgrid .flacsd47、建模的过程中， 或同一个文件中有重复应用的命令， 就可以用宏 命令了。答：宏命令： macro 名称 命令的内容例： macro shiping plot con zdip outline on用 shiping 代替后面的引号里的东西。在建模中只要有相同的都可以 用 shiping 代替，缩减了书写命令的时间，整体上也能让人看的更明 白。注意，你命名的符号最好跟 falc 自带的命令不要重复，处理方法 是，名称后带数字。如 shiping1；如果你想修改宏对象，改变命令的内容，在名称上加单 引号就是替换的意思。例 macro shiping1 新

48、的命令。48、group 的用法： group直接加到模型的后面相当于给整个模型命名； group name range x xyy zz给一个范围的体命名 。 例：在 gen zone brick size group name49、plot 中 add 与 set 的区别add 原来有的加新的用 addset 原来没有的，设置用 setplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20plot set center 2.5 4.2 4.0plot add bound behindplot

49、add bcont szz planeplot add axesplot show 。我们想显示别的图像， 又不想把这个弄好的图像丢掉， 则用 plot create name 建立一个图像， 在建立新的图像时， plot create name1（新的名） 显示以前的图像，那么输入如下命令： plot current name。另一种方法，输入： plot show name50、不同的区域定义不同的本构模型model mohr range z=0,5prop ClayeyGravel range z=0,5在定义参数后面加个 range 范围，为了避免在输入命令时总要通过 range全局坐

50、标域的方法来限定该命令执行的区域，我们可以先通 过 range 命令来给不同的区域命名，然后就可以用 range区域名的 方法来限定范围了。 这样就避免了每次输入命令要限定在该区域执行 时，都要输入一遍坐标范围的麻烦。具体如何给区域命名，以及日后 如何用区域名来限定该区域，可以看下面这个例子： range name trench x=-1,1 y=0,4 z=-2,2 model null range trench51、右手法则：建立模型时，各点的排放时有规律的。 方法如下：拿出右手伸直了 4指合并指向就是 p0到 p1的方向，掌心 指向就是 p0到p2的方向，拇指与 4指垂直的方向就是 p0

51、到p3的方 向。很多人说网格后的云图有带齿什么的， 改进的方法就是重复用下面的 命令Gen ADJust用此命令来微调网格 , 使之离散平滑。该命令可以连续使用以增进平 滑效果的。52、如果想看模型在计算过程中的变化情况Set movie avi step 1 file 1.avi；定义动画， 动画记录的频率为 1 时步一 频率（计算量大的时候就没必要 1 时步的记录，可以 10,100 时步的 记录，减少电脑的使用） ：Movie start；动画开始附件里这个命令“ hist write 5 vs 2 file dibiao.txt ” 能将结果按时步讲 就算结果输出 hist n= 5 ：N（Nstep）=5 指定了每迭代计算 5 次记录一 次相关的值（如果不指定，则系统默认值为 10）；hist unbal 命令的第二行指定记录最大不平衡力（ the maximum unbalanced force）；hist gp zdisp 4,4,8 记录坐标值为 （