18、如何能把一个PLOT的图像数据导出来以便用其他软件绘图?
答:
用setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
19、关于dd、dip的理解:
dip-directionangle,measuredintheglobalxy-planeclockwisefromthepositivey-axis这个是从y轴到x顺时针方向所转过的角度;
dipdipangle,measuredinthenegativez-directionfrom theglobalxy-plane这个是xy平面与z轴负方向的夹角;
dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线;
dd指的是从y轴正方向按顺时针(clockwise)方向转向所要确定面的法线方向在xy平面上的投影的夹角,然后再确定dip,dip指的是从xy平面转向所要确定的平面的角。
20、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没试用过(范围外的东西,想想用到圆上取范围非常好)。
21、如何调用FISH文件?
新手调用FISH文件出错的原因:
是:
文件名为fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的)。
修改方法是:
将fishcall.fis.txt里面的内容全部复制到新的记事本里,再按保存,保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis,类型也就正确了。
22、如何理解Tension-pshear-ptension-nshear-n表示的含义?
在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和tension分别表示因受剪和受拉而处于塑性状态;n表示now,
p表示previous,即分别表示在现在和以前的循环当中处于塑性状态。
23、绘制相关曲线命令的一些解释。
1)在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负);
2)"-"表示对其值作"mirror",镜像。
比如上面的hist1gpydisp的值是0~100,那么vs-1就变成-100~0.以此类推;
3)为什么要这么做呢?
主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内。
24、如何查看剖面上任一点的位移值?
printgpdisprangex()y()z()或printgpdisprangeid
plotgpdisprangex()y()z()或plotgpdisprangeid
25、问a式与b式有何区别?
fixxyzrangex-0.10.1..........................a
fixxrangex-0.10.1................................b
a表示在x=0这个平面是固定的;
b表示在x=0这个平面只x方向不能发生运动。
26、initial与apply有何区别?
initial初始化命令,如初始化计算体的应力状态等;
apply边界条件限制命令,如施加边界的力、位移等约束等;
initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件不会发生变化。
27、如何查看各个时段不平衡力的具体数值?
采用his来记录计算,包括位移应力等命令hisunbal;
hisgp(zone)zdisrange(000)或者id=?
;
导出数据命令hiswritenvsmbegin时步end时步+filefilename.his
n表示纪录的idm表示时步;
要导出不平衡力的具体数值
hisunbal
step100000orsolve
hiswrite1vsstepbegin1end1000file123.his
使用上述命令就可以查看各个时步下的不平衡力的具体数值。
28、如何显示变形轮廓线的命令?
ploskemagf10其中10为放大系数
29、怎样查看模型的矢量图?
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剖面上坐标值
30、怎样查看INTERFACE?
答:
plgridifaceidred
printifaceid;id是interface的编号
31、FLAC怎么调整视图?
答:
ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小,双击屏幕。
32、判断模型中网格点是一个还是两个的方法。
显示模型网格点的ID,如果有两个ID号,说明是两个。
节理面上的网格点应该有两个ID号。
33、改变模型应力边界时FLAC与3DEC的不同之处。
在3DEC中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界时,模型表面的应力为两者应力相加。
34、flac中影响运行时间的因素
(1)单元材料,结构面之间的刚度对比很大;
(2)划分单元尺寸对比很大,尽量避免长细比大于5:
1的细长单元。
35、FLAC中有多少中材料模型。
(1)零模型 :
用于开挖回填的区域
(2)各向同性弹性模型(3)正交各向异性弹性模型(4)横观各向同性弹性模型(5)德鲁克-普拉格塑性模型(6)莫尔-库仑塑性模型 ,一般岩土体的力学行为(7)节理化塑性模型 :
用于研究材料破坏后行为(如边破稳定,屈服煤矿开采)(8)应变硬化/软化模型-库仑塑性模型:
薄板层状材料破坏后力学行为(9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖(10)双屈服塑性模型(11)修正的剑桥粘土模型。
36、对不平衡力的理解。
flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小,该值为百分数,且从不会减小到零,因此,模型不会达到绝对稳定的平横状态,可以根据要求的精度设定R值,模型默认的R值是1e-5。
注意R值小仅表示所有网格点的力都平衡,但模型还可能发生塑性流动,通过plotvel命令查看网格点速度,一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。
37、简单命令解释:
plotaddcontdispouton时时显示云图
plotblockgroupgpnumon显示节点编号
plotinterredske 看接触面的效果.
delrangegroup1not不删除group信息
delrangegroup2删除group信息
inizadd数值大小rangegroupname顺z方向移动
建立接触面interface1facerangex(范围)y(范围)z(范围)或者interface1facerangeplaneorig(xyz坐标)norm(法线方向)watertableclear清除水位线
plsetbackwhite设置背景
setplotbitmapsize(800,510)设置图片大小
plotsetcaptionsize37设置图例数字大小(default=35范围10-50)
plotsetcenter100040
fixxyzrangex-0.10.1表示在x=0这个平面是固定的
fixedboundaryfixxrangex-0.10.1表示在x=0这个平面只x方向不能发生运动rollerboundary。
38、在固定面上施加位移时,fix与apply不能同时使用,应该用ini和fix,才能施加上。
39、你可以用键盘上的上、下、左、右肩在窗口中移动模型。
40、在进行三轴压缩时候,如何施加围压?
答:
使用applysxx=-1e5szz=-1e5rangecylend10,0,0end20,2,0radius1,其中X和Z轴是在我们眼睛正对的面。
42、在进行三轴压缩时候,如何对轴向力进行位移控制?
答:
使用iniyvel1e-7rangey-.1.1iniyvel1e-7rangey1.92.1即可,数值可以变化。
43、在进行边坡计算的时候,为什么用sol,fos命令运行?
答:
应该输入solfosslope3dfos.savassociated(不要也行);再用restore命令调用即可。
44、在FLAC3D计算中,如何记录程序中的值,如果不想时时记录该怎么弄?
答:
用histwtiten1vsn2file.txt,想间隔记录可用sethist_repn
45、想记录接触面上的相关力学参数该怎么办?
答:
histinterinstress/ndis/sstress/sdisp
46、接触面设置参数包括哪几个?
答:
cohdialfricknksten.
47、如何输出单元格的数据?
答;expgrid.flacsd
48、如何做切片?
答:
plotsetplaneorixyznormxyz
plotconzdplane
这两行命令,用来显示指定面上的位移变化云图(把第二行改改,还可以显示应力云图,等其他);
plotconsyyy向应力云图
plotconszzoutlineonz向应力云图加显示网格
plotconszzeffectiveonz向有效应力
plotconszzinter1e4云图显示增量(后面的值越小,图上颜色条带就越细,)
plotconszzmin-10e4显示最小值,图上小于这个数的值都用同一色表示
plotconszzmax10e4显示最大值,图上大于这个数的值都用同一色表示
plotfipoutsideon显示力的箭头
现在看边界,如果是不动面,用fixxrange…如果边界带速率,或位移用ini…
setgrav00,-9.81
inidens1000
设置重力就一定要设置密度,不考虑重力就不要设计密度。
两者共存的,设计程序的时候一定要记的。
setmechforce50
solve
用SET命令设定了最大不平衡力的上限值为50N,即:
计算的最大不平衡力一旦小于此值,我们便人为模型达到了平衡状态。
求解过程终止。
设置这个目的是在一定精度要求下,加快计算(值越大,算的量越少,花的时间也越少)。
49、FLAC结果数据的导出
plsetbackwhite设置背景
setplotbitmapsize(800,510)设置图片大小
plotsetcaptionsize37设置图例数字大小(default=35范围10-50)
mainwinsize11position00.8;commandwindow
plotsetwindowsize0.80.8position0.10.0;graphicswindow
plotsetcenter100040根据情况更改。
50、一段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执行什么定义的函数
51、有关命令的解释
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所施加应力网格点的坐标。
52、FLAC程式使用前准备步骤
步骤1:
依比例画出所欲分析的资料
在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。
步骤2:
换算输入资料成同一单位
需谨慎检查输入资料的单位,如因单位不同而造成过大或过小的值,将会造成FLAC无法计算,而产生ERROR信息。
步骤3:
应用公式简略计算
应用公式或依据经验,简略算出FLAC输出资料的范围,以作为Debug及输出资料分析时验证。
步骤4:
建立x,y坐标与nodei,j之间的关系
在图上距离及深度的关系,建立x,y坐标系统,再由x,y坐标系统,转换与网格关系,为了便于以后输出资料的分析,故应确实掌握网格的位置及其相对应的x,y坐标。
建议在敏感区域使用较密的网格,其它地方则使用较疏的网格,刚开始执行程序时,不宜使用网格太大的网格数目,因尽量使网格总数少于1000,以节省时间。
53、FLAC3D程序的编写步骤:
(1)建立x,y坐标与网格的关系,建议使用Gen指示:
Genx1,y1x2,y2,x3,y3x4,y4i=i0,i1j=j0,j1
(2)设定材料性质:
prop
(3)设定外力:
SetGrav,ApplyPressure,inisxx,Syy
(4)设定边界条件:
fix,free
(5)求起始的应力平衡:
solve
(6)储存:
Save
54、建模的过程中,或同一个文件中有重复应用的命令,就可以用宏命令了,
宏命令:
macro名称‘命令的内容’
例:
macroshiping’plotconzdipoutlineon’
用shiping代替后面的引号里的东西。
在建模中只要有相同的都可以用shiping代替,缩减了书写命令的时间,整体上也能让人看的更明白。
注意,你命名的符号最好跟falc自带的命令不要重复,处理方法是,名称后带数字。
如shiping1;如果你想修改宏对象,改变命令的内容,在名称上加单引号就是替换的意思。
例macro’shiping1‘‘新的命令’。
55、很多时候我们要对单一的对象进行操作,例如进行模拟运算后,只想单独的显示一部分的内容,单一的土层,锚杆,等就要用到group。
group给对象命名,
例:
在genzonebrick------size----groupname(name随便取,不要跟flac的已有的命令重合就行)这里的意思就是对这个体进行了命名,接下来想把这一块模型挖掉加:
modelnullrangegroupname
group的用法:
①group直接加到模型的后面相当于给整个模型命名;
②groupnamerangexxyyzz给一个范围的体命名。
56、上面说过建模一般只建一半然后用reflect命令(镜像)。
命令:
reflectnormalxyzoriginxyz(根据法线和过一点建的面)最常用至于
实际使用过程中有人问道:
normal-100与normal100的区别,我试过两者的效果是一样的,没区别,虽然方向不同,但表达的是同一个面。
reflectddadipboriginxyz这个就是用里建立斜面镜像的。
有一点注意了,镜像的时候如果是模型镜像的话,相当于在建立一个模型,故,在reflect的前面要加建模的关键字,genzonereflect……
plotsetrotation20030;设置观测角度的
plotsetcenterauto;设置观测中心点
plotsetmagnification1.5;设置放大倍数的
plotaddsurfaceyellow
plotaddaxes
plotshow;显示
这一整部分就是我开头说的要自己整理出来的句段了,等以后用到的时候直接插入,作用就是显示模型的。
任何模型的显示都可以用他来表示。
很好很强大。
如果还要显示其他功能就把相应的命令在插入上面那部分就行了。
这个就是我常用来查看模型的命令,大家也可以自己整理一份。
57、add与set的区别
add原来有的加新的用add
set原来没有的,设置用set
plotcreateGravV
plotsetplanedip=90dd=0origin=3,4,0
plotsetrot15020
plotsetcenter2.54.24.0
plotaddboundbehind
plotaddbcontszzplane
plotaddaxes
plotshow。
我们想显示别的图像,又不想把这个弄好的图像丢掉,则用plotcreatename建立一个图像,在建立新的图像时,plotcreatename1(新的名)显示以前的图像,
那么输入如下命令:
plotcur