Flac3D命令--完整经典版 下载本文

实例分析命令:

1. X,Y,Z旋转Shift+X,Y,Z反向旋转

Genzone……;model……;prop……(材料参数);setgrav0,0,-9.81(重力加速度)

plotaddblockgroupredyellow把在group中的部分染成红色和黄色 plotaddaxesblack坐标轴线为黑色;printzonestress%K单元应力结果输出 inidens2000ranzab(设置初始密度,有时不同层密度不同);ini……(设置初始条件);fix……(固定界面)

setplotjpg;setplotquality100;plothardfile1.jpg图像输出(格式、像素、名称)

plotsetmagf1.0视图的放大倍数为1.0;ploconszzz方向应力云图

2. inizadd-1rangegroupone群one的所有单元,在z方向上向下移动1m;然后合并命令

genmerge1e-5rangez0此命令是接触面单元合并成一个整体,1e-5是容差

3. (基坑开挖步骤):Step1:createinitialmodelstate(建立初始模型)

Step2:excavatetrench(开挖隧道)

4. groupToprangegroupBasenot定义(群组Base以外的为)群组Top 5. plotblogro使得各个群组不同颜色显示

6. (两个部分间设置界面;切割法):genseparateTop使两部分的接触网格分离为

两部分;interface1wrapBaseTop在(Base和Top)这两部分之间添加接触单元;plotcreateview_int显示,并创建标题view_int;plotaddsurface显示表面;plotaddinterfacered界面颜色红色

7. (简单的定义函数及运行函数)new;defsetup定义函数setup;numy=8定义常

量numy为8;depth=10.0定义depth为10;end结束对函数的定义;setup运行函数setup

8. (隧道生成)上部圆形放射性圆柱及下部块体单元体的建立,然后镜像。 9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度?— 10. 渐变应力施加:

applynstress-1e6gradient0,0,1e5rangez3.464,0planedip60dd270origin.100;施加法向应力:applynstress-1e6rangeplanedip60dd270origin.100 11. dipdd确定平面位置使用:(纠结) 12. printgppositionrangeid=14647输出节点坐标 applysxx-10e6gradient0,0,1e5rangez-100,0在这个求解方程中,z为变量,所以?xx13.

65为:?xx=-10?10+10?z;原点(0,0,0)

14. freexrangex-.1.1z6.910.1放松x=0平面上,z=7,10这一部分在x方向的约束(可以

在此处产生破坏) 15. 体积模量K和剪切模量G与杨氏模量及泊松比v之间的转换关系如下: 16. 一般而言,大多数问题可以采用FLAC3D默认的收敛标准(或称相对收敛标准),即当体系

-5

最大不平衡力与典型内力的比率R小于定值10;(也可由用户自定义该值,命令:

SETmechratio

17. 所谓体系最大不平衡力,是指每一个计算循环(或称计算时步)中,外力通过网格节点传递

分配到体系各节点时,所有节点的外力与内力之差中的最大值; 所谓典型内力,则是指计算模型所有网格点力的平均值。 18. 要遵循由简单到复杂、由少网格到多网格的思路。 19. 混合离散化法的基本原理是通过适当调整四面体应变率张量中的第一不变量,来

给予单元更多体积变形方面的灵活性。 20. genzoneref以z=0平面为对称面镜像生成网格;genzonerefgenzonerefdip90以y=0面为对称

面镜像生成网格;genzonerefdip90dd90以x=0平面为对称面镜像生成网格 21. 球体模型建立:make_sphere函数循环使用。计算球心至节点距离、获取节点坐

标、调整坐标值得到球面至球心距离、完成八分之一球模型。变换节点内存地址! 22. ①变形云图:(位移、应力、速度等等);plotaddcontszzoutonshadeon打开szz的等

值线图(outline、effective网格、有效的…显示);PLOTconszzouton;PLOTconzdouonmagf20(块云图bcon不能跟magf放大倍数);ploconzdispoutlineonshadeon打开zdisp的等值线云图

②变形矢量图:(变形方向、大小PLOTskdisscale0.07其中scale用来改变箭头大小,默认是0.05);

③塑性区分布printzonestate都可以用来显示模型的塑性区 命令显示那些应力符合屈服准则的区域(或称塑性区)

应力正位于屈服面,或说正处于破坏状态时以shear-n或tension-n标识; 曾进入过屈服状态,但现已经退出以shear-p或tension-p标识 根据塑性区标识判断破坏机制是否在起作用,可按下述步骤进行:

④变量监测(历史跟踪):Histgpzdisp1;Histgpszz2输出2和1的hist关系图用:

plothis-2vs1如果应力值都为负数,则加-2是为了应力坐标轴为正!默认id=1234…

依次(histid=Gpzdisxyz);histwrite7v8file6-3hist.txt使用该命令,程序会在默认文件夹中生成一个名为6-3hist.txt的文本文件,读者可以打开这个文本文件观察所得的数据结果

⑤切片设置:剖面(plotsetplaneori01.50norm010)、云图命令后面加plane、加网格(ske)加矢量(dis)、坐标系(axe)等修饰效果。

23. 输出信息:printzonestress在命令窗口中会输出模型全部单元的6个方向应力数值

printgpdis该命令会显示模型中全部节点三个方向的变形大小

SETlogon在命令行中使用上述命令后,程序在默认目录自动建立一个名为flac3d.log的文件,并开始记录命令窗口中的所有信息,直到用户设置log状态关闭(setlogoff)为止 rest6-3.sav setlogon setlogfile6-2.log printzonestress

printgpdisrangeidaanyidbany setlogoff之后重新打开6-2.log

24. 最大不平衡力不代表没有了塑性流动,只是不平衡力趋于零,并不完全归零。 记录不平衡力随步数的关系。

histunbal监测不平衡力,并保留历史记录;savet1.sav保存到文件t1; histgpzdisp4,4,8监测网格坐标点(4,4,8)在z方向的位移,并保留历史纪录 plotsetrot20030视图的旋转角度为(20,0,30)

PLOTblockgroup查看计算模型Ctrl+G变成灰色图像

25. 采用2种方式观察网格节点速度::HISTgpvel记录某个节点的速度;PLOTvel

命令绘制完整的速度矢量场图;力平衡时,由于网格节点力不为零就存在网格节点速度,仍有可能产生较大位移。

26. 初始地应力场:通常用的是以下三种方法,即弹性求解法、改变参数的弹塑性求解法以及

分阶段弹塑性求解法。

① 由于为弹性求解,在体系达到平衡时,岩、土体中并未有产生屈服的区域modelelas ② 更改强度参数的弹塑性求解法生成是指求解过程中始终采用塑性模型modelmohr

此法与前述弹性求解方法的不同之处在于,计算达到最终平衡时,岩、土体中可能有产生屈服的区域,较合理。

③ (先把基点处总的力算出,然后渐变,一般是递减),高度z=3密度2000;设置初始

应力:iniszz60e3grad0020e3ranz03;

inisxx30e3grad0010e3ranz03;inisyy30e3grad0010e3ranz03;水平为竖向的一半。如果都在水下,则孔隙水压力为inipp30e3grad00-10e3ranz03渐变的过程

饱和密度、干密度、孔隙率以及饱和度间关系式为:?s=?d+ns?w(存在问题)??

27. ID号查询通过标题栏查看或者执行【Plotitems】【/1Blockgroup】【/Modify】命令中grid_pointID’s 28. 建立接触面单元的各类方法(桩的接触面桩端和桩侧最好采用不同id号的接触面单元)

接触面参数的选取见文件155页9.51。接触面有关的常用命令9.7162页。

① 移来移去法:两个模型分开,一个建立接触面单元,然后移动使两个模型合并。 inizadd-1rangegroupone移动命令。

接触面单元interface()facerangecylinderend1(0,0,-4.9)end2(0,0,-5.1)radius.31查看施加结构单元那部分命令(施加接触面单元),—柱体两端面圆心及半径。 ②导来导去法:总模型;删掉其他群组,留下一个群组建立接触面单元(文件1);总模型;删除留下的群组delrangroup2,expgrid导出(文件2);restore接触单元文件1,导入文件2,impgrid文件2。.

interface1facerangex1y12z12 interface1facerangex2y12z12 interface1facerangex12y1z12 interface1facerangex12y12z1 interface1facerangex12y12z2

定义群组方法!

Genzonebricksize333; group2rangex12y12z12; group1rangegroup2not

删除其它群组只保留群组2delrangroup2not

Plotblogrorangegroupdam只显示dam组的单元

③ 切割法:genseparategroup1(表示在1组上建立面);int1wrapgroup1group2,