好文档 - 专业文书写作范文服务资料分享网站

FLAC3D命令流(整理版)

天下 分享 时间: 加入收藏 我要投稿 点赞

实用标准文档

3.plot inter red ske 看接触面的效果.

4.del range group 1 not 不删除group信息 5.del range group 2 删除group信息 6.ini z add ** range group ** (name) 7.interface 1 face range x(*,*)y(*,*)z(*,*)或者 interface 1 face range plane orig(x y z 坐标)norm(法线方向) 建立接触面

8.water table clear 清除水位线 9.pl set back white 设置背景

10.set plot bitmap size (800,510) 设置图片大小 11.plot set caption size 37 设置图例数字大小 12.plot set center 100 0 40 设置图片中心

13.fix x y z range x -0.1 0.1 表示在x=0 这个平面是固定的

14.fix x range x -0.1 0.1 表示这个平面只x方向不能发生运动

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、在进行三轴压缩时候,如何对轴向力进行位移控制?

答:使用 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 n

44、想记录接触面上的相关力学参数该怎么办? 答:hist inter i nstress /ndis/sstress/sdisp

45、接触面设置参数包括哪几个?

文案大全

实用标准文档

答:coh dial fric kn ks ten.

46、如何输出单元格的数据? 答;expgrid .flacsd

47、建模的过程中,或同一个文件中有重复应用的命令,就可以用宏命令了。

答:宏命令: macro 名称 ‘命令的内容’

例:macro shiping ’ plot con zdip outline on’

用shiping 代替后面的引号里的东西。在建模中只要有相同的都可以用shiping代替,缩减了书写命令的时间,整体上也能让人看的更明白。注意,你命名的符号最好跟falc自带的命令不要重复,处理方法是,名称后带数字。

如shiping1;如果你想修改宏对象,改变命令的内容,在名称上加单引号就是替换的意思。 例macro ’shiping1‘ ‘新的命令’。

48、group的用法:

①group直接加到模型的后面相当于给整个模型命名; ②group name range x xyy zz给一个范围的体命名 。 例:在gen zone brick ------ size ---- group name

49、plot中add与set的区别 add原来有的加新的用add set原来没有的,设置用set plot create GravV

plot set plane dip=90 dd=0 origin=3,4,0 plot set rot 15 0 20

plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axes plot show 。

我们想显示别的图像,又不想把这个弄好的图像丢掉,则用plot create name建立一个图像,在建立新的图像时,plot create name1(新的名)显示以前的图像, 那么输入如下命令:plot current name。另一种方法,输入:plot show name

文案大全

实用标准文档

50、不同的区域定义不同的本构模型 model mohr range z=0,5

prop ClayeyGravel range z=0,5

在定义参数后面加个range范围,为了避免在输入命令时总要通过range+全局坐标域的方法来限定该命令执行的区域,我们可以先通过range命令来给不同的区域命名,然后就可以用range+区域名的方法来限定范围了。这样就避免了每次输入命令要限定在该区域执行时,都要输入一遍坐标范围的麻烦。具体如何给区域命名,以及日后如何用区域名来限定该区域,可以看下面这个例子: range name trench x=-1,1 y=0,4 z=-2,2 model null range trench

51、右手法则:建立模型时,各点的排放时有规律的。

方法如下:拿出右手伸直了4指合并指向就是p0到p1的方向,掌心指向就是p0到p2的方向,拇指与4指垂直的方向就是p0到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 记录坐标值为(4,4,8)的节点Z方向的位移(这个用的多,查看收敛量时就是这个命令,上面两个个人认为用处不大,也很少用过)。

53、应变硬化/软化模型SS,如何选择参数表

文案大全

实用标准文档

Pro ftab 1 ctab 2 dtab 3

Table 1 0 45 .05 42 .1 40 1 40 Table 2 0 2e6 .05 1e6 .1 5e5 1 5e5 Table 3 0 10 .05 3 .1 0

答:对于SS模型来说,材料需要硬化或者软化,也就是体现在内聚力和摩察角的不断变大或变小。table 1就是在0, 0.05, 0.1, 1的应变/位移值上体现出内聚力的变化。table 2、3同理。点的选取决定了应力—应变曲线的形状,她的数量可以是任意的。

应力与应变之间的关系

塑性应变plastic strain 与参数之间的关系,自己随便定义

gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12

gen zone reflect norm 1,0,0 gen zone reflect norm 0,0,1

model sspro den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10

pro ftab 1 ctab 2 dtab 3table 1 0 45 .05 42 .1 40 1 40 table 2 0 2e6 .05 1e6 .1 5e5 1 5e5

table 3 0 10 .05 3 .1 0fix x y z range y -.1 .1 fix x y z range y 3.9 4.1

ini yvel 2.5e-5 range y -.1 .1 ;施加恒定速度测试反力(应变控制)ini yvel -2.5e-5 range y 3.9 4.1 def ax_str

文案大全

实用标准文档

str = 0

pnt = gp_head

loop while pnt # null

if gp_ypos(pnt) < 0.1 then

str = str + gp_yfunbal(pnt) endif pnt = gp_next(pnt) endloop

ax_str = str / pi ; cylinder radius = 1 end

hist n 1

hist gp ydisp 0,0,0 hist ax_str

hist gp xdisp 1,1,0

plot hist -2 vs 1 ;axial stress vs axial disp. step 5000

save beforeplzones.sav

; Plot of plastic region as zones with strain > 0.2 def ShowPlasticZones zp = zone_head

loop while zp # null

if z_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 plo crea qqq

plo add surf red range group yield plo add axes green plo set rot 123 313 3 plo set mag 1.5 ret

54、查看网格塑性区

文案大全

FLAC3D命令流(整理版)

实用标准文档3.plotinterredske看接触面的效果.4.delrangegroup1not不删除group信息5.delrangegroup2删除group信息6.inizadd**rangegro
推荐度:
点击下载文档文档为doc格式
3dtj43udbf3uh255c6he20sz532aec00c99
领取福利

微信扫码领取福利

微信扫码分享