四年级上册数学课件 6.《参观花圃》北师大版(2014秋) (共12张PPT)

精心整理1、怎样查看模型？
答：plotgrid可以查看网格，plotgridnum可以查看节点号。

2、请问在圆柱体四周如何施加约束条件？
答：可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2<r1，其实就是选择range的时候选两个圆柱的差，即得到边界。

命令流如下：fixxrangeend1100end2140rad1end1100end2140rad1
3、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？
答：用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

4、用命令建立模型后，如何显示点的坐标？
答：使用ploblogrogpnumon命令
5、关于gauss_dev对性质进行高斯正态分布的问题？
答：根据手册上的说明：下面的命令设定一个平均摩擦角为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&di m142fill
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
和
fixzrangez-0.10.1
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
histzonesyy0,1,0
histzonesyy1,1,0
histwrite1file原点位移.dat
histwrite2file中部单元应力.dat
histwrite3file边缘单元应力.dat
histwrite1table1
histwrite2table2
histwrite3table3
histwrite23vs1file单压.dat
step3000
11、plot显示x、y坐标图
histn1
histgpydisp0,0,0
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的话就明白"mirror"的含义了,比如上面的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
xtable(n+10,k)=rr
ytable(n+10,k)=0.0
end_case
end_loop
end_loop
end
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)
xx4=xtable(n+11,k)
yy4=ytable(n+11,k)
zz4=zz3
command
WATERTABLEfacexx1,yy1,zz1xx2,yy2,zz2xx3,yy3,zz3 WATERTABLEfacexx1,yy1,zz1xx3,yy3,zz3xx4,yy4,zz4
end_command
end_loop
end_loop
end
aux2
13、FLAC结果数据的导出
plsetbackwhite;设置背景
setplotbitmapsize(800,510);设置图片大小plotsetcaptionsize37;设置图例数字大小(default=35范围10-50) mainwinsize11position00.8;commandwindow plotsetwindowsize0.80.8position0.10.0;graphicswindow plotsetcenter100040
根据情况更改
14、一段fish命令的解释
defaux1定义函数
loopn(1,nptab1)开始一个循环，从1循环到变量nptab1
rr=xtable(1,n)从table1的x列中读第n个数付给变量rr
zz=ytable(1,n)从table1的y列中读第n个数付给变量zz loopk(1,nprof)开始第二个循环，从1循环到变量nprof
case_ofk下面是一个选择结构，利用变量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)
case1当k＝1的时候执行下面语句
xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rr ytable(n+10,k)=-40.0tablen+10的y列第k个值等于变量-40 case2当k＝2的时候执行下面语句
xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rr ytable(n+10,k)=0.0tablen+10的y列第k个值等于变量0
end_case结束选择结构
end_loop结束第二个循环
end_loop结束第一个循环
end函数定义结束
aux1执行什么定义的函数
15、有关命令的解释
（1）applynvel0.1planedip60dd270rangeleft_boun
apply施加边界条件，plane法向速度为0.1，plane面的位置：倾角60（在xoy平面从-Z轴方向逆时针旋转），dd走向270度，从+y方向绕z轴顺时针旋转；
（2）iniszz-5e4grad00-1e4
ini初始边界条件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4是range范围起点坐标的SZZ 值，x,y,z所施加应力网格点的坐标。

16、渗流示例
;------------------------------参数部分
bulkmodulus,K390MPa体积模量（土质比较硬）
shearmodulus,G280MPa剪切模量
soildrydensity,ρd1200kg/m3土的密度
waterdensity,ρw1000kg/m3水的密度
walldensity,ρwal1500kg/m3档土墙的密度
permeability,k10?12m2/Pa-s渗透系数
porosity,n0.3孔隙率
fluidbulkmodulus,Kf2.0GPa流体的体积模量
;-----------------------以下是命令流
;------------------------------------------------------
;Excavationinasaturatedsoil
;------------------------------------------------------
configfluid;设置流体算法
;---geometricalmodel---建立模型
genzonebrickp11200p20120p30012size121212rat111
groupsoil;设置土组
groupexcavaterangex04y04z05;设置开挖范围组
groupwal1rangex45y05z07;设置档土墙1组
groupwal2rangex04y45z07;设置档土墙2组
groupwallrangegroupwal1anygroupwal2any;档土墙1和2合并为一个档土墙组
;---fluidflowmodel---流体模型
modelfl_iso;各项同性流体
propperm1e-12poro0.3;设置渗透系数和？？比or？？率
inifdensity1e3;初始化水密度
inifmod2.0e9ftens-1e-3;初始化水的体积模量及？？
modelfl_nullrangegroupexcavate;开挖部分没有水流
modelfl_nullrangegroupwall;档土墙不渗水
inipp0grad001e4;初始化孔隙水压力
fixpprangez-0.10.1;约束应该是整个模型的底部吧？
fixpprangex-0.14.1y-0.14.1z4.95.1;开挖部分也给约束住了，而其他的地方不管
;---mechanicalmodel---力学模型
modelelas;弹塑性
propbul3.9e6shea2.8e6;设置体积模量和剪切模量，这两个必须，有一次我只设置了体积模量，就错的一塌糊涂
modelnullrangegroupexcavate;力学模型的开挖
inidensity1.2e3;初始化密度，应该是土的密度
inidensity1.5e3rangegroupwall;初始化墙（档土墙）的密度，
fixxrangex-.1.1;施加约束，据说差分法约束0的时候，必须取一个-0.1<0<0.1,范围fixxrangex11.912.1;施加约束，或者叫做边界条件
fixyrangey-.1.1;施加约束
fixyrangey11.912.1
fixzrangez11.912.1
;initialtotalstresses出世总应力
iniszz0grad00-1.5e4;z方向的应力大一些
inisxx0grad00-1.2e4;水平方向的，小一些
inisyy0grad00-1.2e4;
applynstress0grad00-1.2e4rangex0.04.0y3.94.1z0.05.0
applynstress0grad00-1.2e4rangex3.94.1y0.04.0z0.05.0
applynstress-7.5e4rangex0.04.0y0.04.0z4.95.1
;---setting---设置
setgravity0010;设置重力加速度
;---initialstate---;初始化状态
solveforce1;checkinitialequilibrium求解平衡
;---histories---记录
sethist_rep40;记录步为40
histfltime;记录时间
histgppp006;记录某点的孔隙水压力
histgpxdis400;水平位移
histgpxdis402;记录水平位移
histgpxdis420
histgpzdis005;记录竖向位移
histgpzdis205
histgpzdis405
histgpzdis225
histgpzdis425
histgpzdis445
histgpzdis1001
histgpzdis1002
;
;---excavation---开挖
setfluidoff;关闭水
;applyporepressureatwalls在墙上施加孔隙水压力applynstress0grad00-1.e4rangex0.04.0y3.94.1z0.05.0 applynstress0grad00-1.e4rangex3.94.1y0.04.0z0.05.0 applynstress-5.e4rangex0.04.0y0.04.0z4.95.1 solve;force1;求解
saveexc1.sav;保存
;
;---drainage---排水
applyremovenstress;撤掉刚才的压力defrelaxsetup;定义一个释放函数的参数
step0=step
end
relaxsetup;调用
defrelax;定义一个释放函数
rstep=step-step0
ifrstep<ncycthen
relax=1.0-(float(rstep)/float(ncyc))
else
relax=0.0
endif
end;结束
setncyc=1000;先赋值，随后调用
applynstress0grad00-1.e4histrelax&
rangex0.04.0y3.94.1z0.05.0
applynstress0grad00-1.e4histrelax&
rangex3.94.1y0.04.0z0.05.0
applynstress-5.e4histrelax&
rangex0.04.0y0.04.0z4.95.1
cycncyc
solve;求解
saveexc2.sav;保存
；
;---percolation---渗透
setfluidon;设置水算法
fixpp0rangex-0.14.1y-0.14.1z4.95.1;施加孔隙水压力
cyc9000;循环
saveexc3.sav;保存
plotcreateexcav;建立一个显示的视图一下内容都是显示内容plotsetrot2000195
plotsetcent666
plotsetdist39.18
plotsetmagn0.8
plotsetplaneori000normal0-10
plotaddcontppouton plotaddblockgrouplgrarangegroupwall
plotaddcontppint10000max110000rangex04y04z57 plotaddflowplane
plotshow
Ret
17、用什么命令可以位移图中显示最大位移点的坐标？答：用fish编了一个，算法好像有点烦，不过目的可以实现你可以参考一下。

defget_gp_maxdisp
gp0_disp=gp_xdisp(gp_head)*gp_xdisp(gp_head)
gp0_disp=gp0_disp+gp_ydisp(gp_head)*gp_ydisp(gp_head) gp0_disp=gp0_disp+gp_zdisp(gp_head)*gp_zdisp(gp_head) gp0_disp=sqrt(gp0_disp)
p_gp=gp_head
;找最大值
loopwhilep_gp#null
gp_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)
ifgp_disp>gp0_disp
gp0_disp=gp_disp
endif
p_gp=gp_next(p_gp)
endloop
;找最大值的坐标
p_gp=gp_head
ss=0
loopwhilep_gp#null
gp_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)
ifgp_disp=gp0_disp
ss=ss+1
xtable(1,ss)=gp_xpos(p_gp)
ytable(1,ss)=gp_ypos(p_gp)
endif
p_gp=gp_next(p_gp)
endloop
end
get_gp_maxdisp
18、如何显示计算后模型的变形状态？
答：plskmagfac**(一个放大的数字)
19、如何查看剖面上任一点的位移值
答：printgpdisprangex()y()z()或printgpdisprangeid
plotgpdisprangex()y()z()或plotgpdisprangeid
20、group命令划分斜面
答：groupgroupnamecolourrangeplaneabovedd**dip**origin***planebelowdd**dip**origin* **
dd和dip就是两个斜面的倾向和倾角,origin是这个斜面通过的某个点,这个命令就可以group两个斜面间的区域
21、ansys导入flac的时候，如何划分一下group
答：在ansys中对不同group用不同的材料号（mat）定义，再用ansys2flac之类的程序转换成flac的.dat文件就可，group自动生成。

22、fixxyzrangex-0.10.1..........................a
fixxrangex-0.10.1................................b
请问a式与b式有何区别？
答：a表示在x=0这个平面是固定的
b表示在x=0这个平面只x方向不能发生运动。

23、如何知道每个单元的ID？
答：用鼠标双击单元的表面，可以知道单元的ID和坐标。

24、如何进行切片？
答：plosetplaneori(点坐标)norm(法向矢量)
ploconszplane(显示z方向应力的切片)
这两行命令，用来显示指定面上的位移变化云图（把第二行改改，还可以显示应力云图，等其他）；
plotconsyyy向应力云图
plotconszzoutlineonz向应力云图加显示网格
plotconszzeffectiveonz向有效应力
plotconszzinter1e4云图显示增量（后面的值越小，图上颜色条带就越细，）plotconszzmin-10e4显示最小值，图上小于这个数的值都用同一色表示plotconszzmax10e4显示最大值，图上大于这个数的值都用同一色表示
plotfipoutsideon显示力的箭头
现在看边界，如果是不动面，用fixxrange…如果边界带速率，或位移用ini…setgrav00，-9.81
inidens1000
设置重力就一定要设置密度，不考虑重力就不要设计密度。

两者共存的，设计程序的时候一定要记的。

setmechforce50
solve
用SET命令设定了最大不平衡力的上限值为50N，即：计算的最大不平衡力一旦小于此值，我们便人为模型达到了平衡状态。

求解过程终止。

设置这个目的是在一定精度要求下，加快计算（值越大，算的量越少，花的时间也越少）。

25、如何在程序中进行暂停，并可恢复计算？
答：在命令中加入pause命令，用continue进行继续。

在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正，而不是等到结果出来。

26、如何跳过某个计算步？
答：在计算中按空格键跳过本次计算，自动进入下一步
27、Fish是什么东西？Fish是否一定要学？
答：是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能。

Fish可以不用学，需要的时候查Mannual获得需要的变量就可以了。

28、如何调用一些可选模块？
答：configdyn(fluid,creep,cppudm)
29、range的命令：
①取圆球面的范围，就是同一个中心，不同半径之间夹的那一部分，用在给球面施加力，约束时用该命令，命令：
rangeannuluscenterxyz（中心坐标）radiusr1r2（两半径），想要表达的是球面时，r1r2取值相差小点即可，如radius2.93.1这时范围就是半径为3的面。

②圆柱的环面施加力或约束时也用到该命令。

rangecylinderend1xyzend2xyzradiusr(end是表示柱体断面的中心点坐标) rangeplaneabove表示的是z大于0的上部结构
rangeplaneblow表示的是z小于0的下部结构
rangeddadipboriginxyz表示的是走向角度a，倾向角度b的面
dd表示，当面与xoy面相交时的垂线与y轴正向的夹角，顺为正，说白话，就是面的指向方向面是朝那个方向的意思。

dip面与xoy面的夹角，与z轴负向为正，origin表示面经过xyz点；rangenormalxyz法线方向为xyz的面
rangenot没试用过（范围外的东西，想想用到圆上取范围非常好）。

30、如何调用FISH文件？
31、问a式与b式有何区别？
答：fixxyzrangex-0.10.1..........................a
fixxrangex-0.10.1................................b
a表示在x=0这个平面是固定的；
b表示在x=0这个平面只x方向不能发生运动。

32、如何查看各个时段不平衡力的具体数值？
答：采用his来记录计算，包括位移应力等命令hisunbal；
hisgp(zone)zdisrange(000)或者id=？；
导出数据命令hiswritenvsmbegin时步end时步+filefilename.his
n表示纪录的idm表示时步；
要导出不平衡力的具体数值
hisunbal
step100000orsolve
hiswrite1vsstepbegin1end1000file123.his
使用上述命令就可以查看各个时步下的不平衡力的具体数值。

33、如何显示变形轮廓线的命令？
ploskemagf10其中10为放大系数
34、怎样查看模型的矢量图？
plgridvel(disp)red
ploblogro可以看到不同的group的颜色分布（看模型的样子）；
plogpfixredsk看模型的边界情况
plofapredsk看模型的体力分布
plocondis(xdis,ydis,zdis)看模型的云图：位移
ploconsz(sy,sx,sxy,syz,sxz)看模型的云图：应力
plodis(xdis,ydis,zdis)看模型的矢量图
ploinfo看模型有多少单元、节点
plotgrid可以查看网格，plotgridnum可以查看节点号。

plotfix查看边界条件
plotgridmark查看预先留设的开挖边界
plsyy(sxxxdisydissig1..)fill；查看模型的云图。

注意，带fill关键字就是查看云图，不带fill关键字就是查看等直线图。

ploaddske剖面的沉降云图上增加网格
ploadddisplane剖面上的矢量
ploaddaxe剖面上坐标值
35、怎样查看INTERFACE？
答：plgridifaceidred
printifaceid;id是interface的编号
36、FLAC怎么调整视图？
答：ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位，如果要缩小，双击屏幕。

37、判断模型中网格点是一个还是两个的方法。

?答：显示模型网格点的ID，如果有两个ID号，说明是两个。

节理面上的网格点应该有两个ID号。

38、flac中影响运行时间的因素
答：（1）单元材料，结构面之间的刚度对比很大；
（2）划分单元尺寸对比很大，尽量避免长细比大于5：1的细长单元。

39、简单命令解释：
答：1.plotaddcontdispouton时时显示云图
2.plotblockgroupgpnumon显示节点编号
3.plotinterredske????????看接触面的效果.
4.delrangegroup1not不删除group信息??
5.delrangegroup2删除group信息
6.inizadd**rangegroup**（name）
7.interface1facerangex（*，*）y（*，*）z（*，*）或者interface1facerangeplaneorig(xyz 坐标）norm（法线方向）建立接触面
8.watertableclear清除水位线
9.plsetbackwhite设置背景
10.setplotbitmapsize(800,510)设置图片大小
11.plotsetcaptionsize37设置图例数字大小
12.plotsetcenter100040设置图片中心
13.fixxyzrangex-0.10.1表示在x=0这个平面是固定的
14.fixxrangex-0.10.1表示这个平面只x方向不能发生运动
40、在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。

41、在进行三轴压缩时候，如何施加围压？
答：使用applysxx=-1e5szz=-1e5rangecylend10,0,0end20,2,0radius1，其中X和Z轴是在我们眼睛正对的面。

42、在进行三轴压缩时候，如何对轴向力进行位移控制？
答：使用iniyvel1e-7rangey-0.10.1iniyvel1e-7rangey1.92.1即可，数值可以变化。

43、在FLAC3D计算中，如何记录程序中的值，如果不想时时记录该怎么弄？答：用histwtiten1vsn2file.txt，想间隔记录可用sethist_repn
44、想记录接触面上的相关力学参数该怎么办？
答：histinterinstress/ndis/sstress/sdisp
45、接触面设置参数包括哪几个？
答：cohdialfricknksten.
46、如何输出单元格的数据？
答；expgrid.flacsd
47、建模的过程中，或同一个文件中有重复应用的命令，就可以用宏命令了。

答：宏命令：macro名称‘命令的内容’
例：macroshiping’plotconzdipoutlineon’
用shiping代替后面的引号里的东西。

在建模中只要有相同的都可以用shiping代替，缩减了书写命令的时间，整体上也能让人看的更明白。

注意，你命名的符号最好跟falc自带的命令不要重复，处理方法是，名称后带数字。

如shiping1；如果你想修改宏对象，改变命令的内容，在名称上加单引号就是替换的意思。

例macro’shiping1‘‘新的命令’。

48、group的用法：
①group直接加到模型的后面相当于给整个模型命名；
②groupnamerangexxyyzz给一个范围的体命名。

例：在genzonebrick------size----groupname
49、plot中add与set的区别
add原来有的加新的用add
set原来没有的，设置用set
plotcreateGravV
plotsetplanedip=90dd=0origin=3,4,0
plotsetrot15020
plotsetcenter2.54.24.0
plotaddboundbehind
plotaddbcontszzplane
plotaddaxes
plotshow。

我们想显示别的图像，又不想把这个弄好的图像丢掉，则用plotcreatename建立一个图像，在建立新的图像时，plotcreatename1（新的名）显示以前的图像，
那么输入如下命令：plotcurrentname。

另一种方法，输入：plotshowname
50、不同的区域定义不同的本构模型
modelmohrrangez=0,5
propClayeyGravelrangez=0,5
在定义参数后面加个range范围，为了避免在输入命令时总要通过range＋全局坐标域的方法来限定该命令执行的区域，我们可以先通过range命令来给不同的区域命名，然后就可以用range＋区域名的方法来限定范围了。

这样就避免了每次输入命令要限定在该区域执行时，都要输入一遍坐标范围的麻烦。

具体如何给区域命名，以及日后如何用区域名来限定该区域，可以看下面这个例子：
rangenametrenchx=-1,1y=0,4z=-2,2
modelnullrangetrench
51、右手法则：建立模型时，各点的排放时有规律的。

方法如下：拿出右手伸直了4指合并指向就是p0到p1的方向，掌心指向就是p0到p2的方向，拇指与4指垂直的方向就是p0到p3的方向。

很多人说网格后的云图有带齿什么的，改进的方法就是重复用下面的命令GenADJust
用此命令来微调网格,使之离散平滑。

该命令可以连续使用以增进平滑效果的。

52、如果想看模型在计算过程中的变化情况
Setmovieavistep1file1.avi；定义动画，动画记录的频率为1时步一频率（计算量大的时候就没必要1时步的记录，可以10,100时步的记录，减少电脑的使用）：Moviestart；动画开始
附件里这个命令“histwrite5vs2filedibiao.txt”能将结果按时步讲就算结果输出histn=5：N(Nstep)=5指定了每迭代计算5次记录一次相关的值（如果不指定，则系统默认值为10）；
histunbal命令的第二行指定记录最大不平衡力（themaximumunbalancedforce）；histgpzdisp4,4,8记录坐标值为（4，4，8）的节点Z方向的位移（这个用的多，查看收敛量时就是这个命令，上面两个个人认为用处不大，也很少用过）。

53、应变硬化/软化模型SS，如何选择参数表
Proftab1ctab2dtab3
Table1045.0542.140140
Table202e6.051e6.15e515e5
Table3010.053.10
答：对于SS模型来说，材料需要硬化或者软化，也就是体现在内聚力和摩察角的不断变大或变小。

table1就是在0,0.05,0.1,1的应变/位移值上体现出内聚力的变化。

table2、3同理。

点的选取决定了应力—应变曲线的形状，她的数量可以是任意的。

应力与应变之间的关系
塑性应变plasticstrain与参数之间的关系，自己随便定义
genzonecylp0000p1100p2040p3001size123012
genzonereflectnorm1,0,0
genzonereflectnorm0,0,1
modelssproden2500bulk2e8she1e8co2e6fric45ten1e6dil10
proftab1ctab2dtab3table1045.0542.140140
table202e6.051e6.15e515e5
table3010.053.10fixxyzrangey-.1.1
fixxyzrangey3.94.1
iniyvel2.5e-5rangey-.1.1；施加恒定速度测试反力（应变控制）iniyvel-2.5e-5rangey3.94.1 defax_str
str=0
pnt=gp_head
loopwhilepnt#null
ifgp_ypos(pnt)<0.1then
str=str+gp_yfunbal(pnt)endif
pnt=gp_next(pnt)
endloop
ax_str=str/pi;cylinderradius=1
end
histn1
histgpydisp0,0,0
histax_str
histgpxdisp1,1,0
plothist-2vs1;axialstressvsaxialdisp.
step5000
savebeforeplzones.sav
;Plotofplasticregionaszoneswithstrain>0.2
defShowPlasticZones
zp=zone_head
loopwhilezp#null
ifz_prop(zp,'es_plastic')>0.2；塑性剪应变es_plastic;塑性拉应变et_plastic z_group(zp)='yield'
else
z_group(zp)='other'
endif
zp=z_next(zp)
endLoop
end
ShowPlasticZones
plocreaqqq
ploaddsurfredrangegroupyield
ploaddaxesgreen
plosetrot1233133
plosetmag1.5
ret
54、查看网格塑性区
Plblockstate
55、如何建立马蹄形隧道
new
genzoneradcylsize58810&
p0000p11002p20100p30010&
p410102p501010p610010p7101010
;第一部分
genzoneradcylsize58810&
p0000p1100-10p20100p31002&
p41010-10p510102p6100-4p71010-4
;第二部分
genzone brick size5810&
p000-10p1100-10p2010-10p300-2p41010-10&
p5010-2p620-2p7210-2
;第三部分
genzonereflectorig000norm100
;镜像
plocreaqqq
ploaddsurfred
plosho
56、如何显示锚杆Plotselgeomnodeoncidon。