flac3d常用命令.docx

上传人:b****5 文档编号:12188267 上传时间:2023-04-17 格式:DOCX 页数:15 大小:25.41KB
下载 相关 举报
flac3d常用命令.docx_第1页
第1页 / 共15页
flac3d常用命令.docx_第2页
第2页 / 共15页
flac3d常用命令.docx_第3页
第3页 / 共15页
flac3d常用命令.docx_第4页
第4页 / 共15页
flac3d常用命令.docx_第5页
第5页 / 共15页
点击查看更多>>
下载资源
资源描述

flac3d常用命令.docx

《flac3d常用命令.docx》由会员分享,可在线阅读,更多相关《flac3d常用命令.docx(15页珍藏版)》请在冰豆网上搜索。

flac3d常用命令.docx

flac3d常用命令

1、最先需要掌握的命令有哪些?

答:

需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。

2、怎样输出模型的后处理图?

答:

File/Printtype/Jpgfile,然后选择File/Print,将保存格式选择为jpe文件。

3、怎样调用一个文件?

答:

File/call或者call命令

4、如何施加面力?

答:

appnstress

5、如何调整视图的大小、角度?

答:

综合使用x,y,z,m,Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键。

6、如何进行边界约束?

答:

fixxran(约束的是速度,在初始情况下约束等效于位移约束)。

7、如何知道每个单元的ID?

答:

用鼠标双击单元的表面,可以知道单元的ID和坐标。

8、如何进行切片?

答:

plosetplaneori(点坐标)norm(法向矢量)

ploconszplane(显示z方向应力的切片)

9、如何保存计算结果?

答:

save+文件名

10、如何调用已保存的结果?

答:

rest+文件名;或者File/Restor

11、如何暂停计算?

答:

Esc

12、如何在程序中进行暂停,并可恢复计算?

答:

在命令中加入pause命令,用continue进行继续。

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

13、如何跳过某个计算步?

答:

在计算中按空格键跳过本次计算,自动进入下一步

14、Fish是什么东西?

Fish是否一定要学?

答:

是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能。

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

15、FLAC3D允许的命令文件格式有哪些?

答:

无所谓,只要是文本文件,什么后缀都可以。

16、如何调用一些可选模块?

答:

configdyn(fluid,creep,cppudm)

17、如何在圆柱体四周如何施加约束条件?

可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其中r2

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

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 解决方案 > 学习计划

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1