FLAC3D常用命令

合集下载

flac3d 常用命令师姐整理。

flac3d 常用命令师姐整理。
pl sur gp on ra y 20 24 ;监测点位移,找点gpid
pl sel geom cid on ra y 22 24 ;监测衬砌内力,找衬砌cid
pl sel geom cid on node off ra y 21 23 ;监测桩号,找桩cid
pl sel rec sres my surfx 0 1 0 ra y 0 24 scale=0.05
pl add fap ra cyl end1 0 0 -17 end2 0 42 -17 r 3.1
pl add sel geom id on
;查局部坐标
pr sel node local ra id=1

pl bl state shear-n ;看屈服
pl add gpfix ;查固定边界条件
pl add su magf 10 ;变形后位移
pr gp disp ra id=12067;查看某点位移
pl bcon prop bu;无刚度错误错误
pr sel node link ra id ;找node连接的link号
pr sel link tar ra id ;找link的目标节点
pr sel link sou ra id ;找link的源节点
pr sel link attach
pr sel link c lind k;查link的刚度
pr sel rec sr my
;pl sel disp magf 100
;pl sel pile fo fx
;pl sel pile mo my
;pl con zdis
pr sel rec surf
pl sel geom node off nodesys on scale=0.06

Flac3D 常用命令流

Flac3D 常用命令流

FLAC3D 常用命令流;模型镜像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一点);沿X轴镜像,通过对称平面法线向量确定对称面gen zone reflect normal 0 0 -1 ;沿z轴镜像;绘图控制pl <add> contour szz outline on;在模型中显示位移-应变曲线hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror";对模型进行压缩实验的方法;即在模型两侧施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1;修改模型的坐标值ini x add -100 y add -100 z add -100;显示云图的同时也显示模型网格轮廓plot add cont disp outline on <gradient/average>;gradient更精确;输入角度、弧度方法pi=π,90°为90.0*degraddef set_valsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint ptA ptB ptC;施加结构单元方法sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &end2=(0.0,25.0,0.0) radius=24.5 notplot add sel geom black black cid on scale=0.03sel node init zpos add -25.0;如何显示某一平面plot create name_planeplot set plane origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom blackplot add axes redplot current name_planeplot show;如何施加法向应力并显示所施加的力apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力plot grid green fap black;显示施加的力;获取某网格点的ID,坐标值def aaget_id=gp_id(gp_near(3,0,0))endaaprint get_idprint gp position range id 3;保存命令流文件set log onset logfile hua.log;文件存放的路径与调用的文本路径相同......set log off;fish函数的说明与解释def find_add ;定义fish函数find_addhead = null ;给head赋值p_gp = gp_head ;第一个网格结点的指针赋给p_gploop_while p_gp # null ;当p_gp值不为null时作循环x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_posif x_pos = 10.0 then ;如果x_pos = 10.0 则(执行)new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址mem(new) = head ;将head类型和数值置于地址为new的fish变量mem(new+1)=p_gp ;第二个变量存储符合条件的zone地址head = new ;将new值赋给headendifp_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp endloop;结束循环end ;结束fish函数;实际上这个fish函数为满足条件(x坐标为10的)的zone的地址开辟一定的地址空间,各地址之间存在一定的联系;找到第一个符合条件的zone地址后,用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址,并用;下面的mem(new) = head ,将第一次开辟的两个变量的第一个变量存储地址head(注意第一次head=null),第二个变量存储第一;个符合条件的zone地址,并将第一个变量的地址赋予head(head=new),第一次循环结束;下次循环,同样开辟两个变量对象,第;一个变量对象记录上次循环开辟的第一个变量的地址,第二个记录第二个符合条件的zone地址,其余循环依次类推,这样子就建立了;一个符合条件的zone地址链条,方便以后使用。

(完整版)FLAC3D常用命令

(完整版)FLAC3D常用命令

1. apply(缩写:app)可用来定义边界条件及初始条件:1)添加应力格式1:apply szz -0.3395e6 range z -0.1 0.1格式1:apply szz -0.3395e6 range group pile格式3:apply szz -0.3395e6 range z -0.1 0.1 group pile格式4:apply nstress 数值range z 2.9 3.1(或3)x 1 2 y 1 22)以一定速度施加位移边界格式1:apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step3)添加边界条件格式:apply szz <数值> grad <梯度> range <范围>示例1:apply szz -1e9 grad 0 0 8.3e5 range z 0 120示例2:apply szz -0.6e6 range z 0.05 0.15 group pile注:<数值>是梯度方向坐标0点的数值,可通过定义坐标范围的上下值与梯度计算得到。

2. range(缩写:ran)通过range功能,可以使命令作用在一定指定范围的目标上;如果一个命令没有使用range来确定范围,则命令对整个模型有效。

1)利用坐标指定一定的范围格式1:range z 0 1格式2:range z 2.9 3.1 x 1 2 y 1 22)利用分组来指定范围格式:range group 13)以上两种的复合格式:range z -0.1 0.1 group pile4)利用id号来指定一定的范围格式:range id 0 10该命令后跟起始id和结束id,这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例:model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

Flac3D命令--完整经典版

Flac3D命令--完整经典版

实例分析命令:1. X ,Y ,Z 旋转 Shift+ X ,Y ,Z 反向旋转Gen zone ……;model ……;prop ……(材料参数);set grav 0,0,(重力加速度) plot add block group red yellow 把在group 中的部分染成红色和黄色plot add axes black 坐标轴线为黑色;print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b (设置初始密度,有时不同层密度不同);ini ……(设置初始条件);fix ……(固定界面)set plot jpg ;set plot quality 100 ;plot hard file 图像输出(格式、像素、名称) plot set magf 视图的放大倍数为;plo con szz z 方向应力云图2. ini z add -1 range group one 群one 的所有单元,在z 方向上向下移动1m ;然后合并命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体,1e-5是容差3. (基坑开挖步骤):Step 1: create initial model state (建立初始模型)Step 2: excavatetrench (开挖隧道)4. group Top range group Base not 定义(群组Base 以外的为)群组Top5. plot blo gro 使得各个群组不同颜色显示6. (两个部分间设置界面;切割法):gen separate Top 使两部分的接触网格分离为两部分;interface 1 wrap Base Top 在(Base 和Top )这两部分之间添加接触单元;plot create view_int 显示,并创建标题view_int ;plot add surface 显示表面;plot add interface red 界面颜色红色7. (简单的定义函数及运行函数)new ;def setup 定义函数setup ;numy = 8定义常量numy 为8;depth = 定义depth 为10;end 结束对函数的定义;setup 运行函数setup8. (隧道生成)上部圆形放射性圆柱及下部块体单元体的建立,然后镜像。

flac3d快捷命令

flac3d快捷命令

1、FLAC3D常见命令:1. FLAC3D是有限元程序吗?答:不是!是有限差分法。

2. 最先需要掌握有哪些?答:需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解。

3. 怎样看模型的样子?答:plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况?答:plo gpfix red5. 怎样看模型的体力分布?答:plo fap red6. 怎样看模型的云图?答:位移:plo con dis (xdis, ydis, zdis)应力:plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图?答:plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点?答:pri info9. 怎样输出模型的后处理图?答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件10. 怎样调用一个文件?答:File/call或者call命令10. 如何施加面力?答:app nstress11. 如何调整视图的大小、角度?答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键12. 如何进行边界约束?答:fix x ran (约束的是速度,在初始情况下约束等效于位移约束)13. 如何知道每个单元的ID?答:用鼠标双击单元的表面,可以知道单元的ID和坐标14. 如何进行切片?答:plo set plane ori (点坐标) norm (法向矢量)plo con sz plane (显示z方向应力的切片)15. 如何保存计算结果?答:save +文件名.16. 如何调用已保存的结果?答:rest +文件名;或者File / Restore17. 如何暂停计算?答:Esc18. 如何在程序中进行暂停,并可恢复计算?答:在命令中加入pause命令,用continue进行继续19. 如何跳过某个计算步?答:在计算中按空格键跳过本次计算,自动进入下一步20. Fish是什么东西?答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学?答:可以不用,需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许文件格式有哪些?答:无所谓,只要是文本文件,什么后缀都可以23. 如何调用一些可选模块?答:config dyn (fluid, creep, cppudm)24 .如何在圆柱体四周如何施加约束条件?答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。

(完整版)FLAC3D常用命令

(完整版)FLAC3D常用命令

1. apply(缩写:app)可用来定义边界条件及初始条件:1)添加应力格式1:apply szz -0。

3395e6 range z -0.1 0。

1格式1:apply szz -0.3395e6 range group pile格式3:apply szz —0.3395e6 range z —0。

1 0.1 group pile格式4:apply nstress 数值 range z 2.9 3.1(或3) x 1 2 y 1 22)以一定速度施加位移边界格式1:apply yvel —1e5 range y -1.9 2.1 ;施加y方向速度—1e5/step3)添加边界条件格式:apply szz <数值> grad <梯度〉 range 〈范围>示例1:apply szz -1e9 grad 0 0 8。

3e5 range z 0 120示例2:apply szz -0。

6e6 range z 0.05 0。

15 group pile注:〈数值〉是梯度方向坐标0点的数值,可通过定义坐标范围的上下值与梯度计算得到.2. range(缩写:ran)通过range功能,可以使命令作用在一定指定范围的目标上;如果一个命令没有使用range 来确定范围,则命令对整个模型有效。

1)利用坐标指定一定的范围格式1:range z 0 1格式2:range z 2。

9 3.1 x 1 2 y 1 22)利用分组来指定范围格式:range group 13)以上两种的复合格式:range z -0.1 0。

1 group pile4)利用id号来指定一定的范围格式:range id 0 10该命令后跟起始id和结束id,这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例:model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

FLAC3D命令流(整理版)

FLAC3D命令流(整理版)

1、怎样查看模型?答:plot grid 可以查看网格,plot grid num 可以查看节点号。

2、请问在圆柱体四周如何施加约束条件?答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。

命令流如下:fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。

4、用命令建立模型后,如何显示点的坐标?答:使用plo blo gro gpnum on命令5、关于gauss_dev对性质进行高斯正态分布的问题?答:根据手册上的说明:下面的命令设定一个平均摩擦角为40度,标准方差是±5%。

则命令如下:prop friction 40 gauss_dev 2 问题:请问gauss_dev 2中的2是如何计算的?如果把±5%改为±10%,则命令应如何写?40×5%=26、 reflect问题问:gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill plot surfgen zone reflect dip 0 dd 90 (对xy面做镜像)gen zone reflect dip 90 dd 90 (对yz面做镜像)(1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线。

FLAC3D常用命令总结

FLAC3D常用命令总结

1、模型建好如何检验所建模型是否有问题注意:第一步建模一定不要建太大的模型,可以建非常小的模型,确保模型的准确性,这样可以尽快开始后边的步骤,等后边计算分析步骤快速做完验证没有大问题,再将准确的模型建好换上,重新计算得到更加准确的结果。

2、如何检验应力边界、位移边界、加载速率等是否准确Plot add gpfix:看位移边界Plo add fap:看应力边界,看加载方向是否准确Plot add vel:看速率,加载速率是否准确3、如何查看运算结果Plo con szz outline on:看z向应力云图,并显示边界网格Plo add disp:看位移Plo con zdisp outline on magf 5:看z向位移云图,变形扩大5倍Plo add bo :添加边界Plo blo state:看是否屈服Plo add stensor:查看大主应力方向Plo con smin:小主应力Plo con smax:大主应力Plo con ssi(shear strain increment):查看剪切应变增量云图Plo bcont smin:查看单元小主应力4、如何做切面并查看切面结果在编辑“Edit“菜单中设置切面位置,在上述3命令后加上plane(pl)即可。

如:Plo con smin pl:查看切面小主应力5、后处理将模型导入到tecplot中做切面,加等值线等,对结果进一步处理,得到更漂亮的一些结果。

其他重要命令总结:1、如何在圆柱体四周如何施加约束条件?可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。

2、如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer 或者什么作图软件绘制就行了。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1. apply(缩写:app)可用来定义边界条件及初始条件:1)添加应力格式1:apply szz -0.3395e6 range z -0.1 0.1格式1:apply szz -0.3395e6 range group pile格式3:apply szz -0.3395e6 range z -0.1 0.1 group pile格式4:apply nstress 数值range z 2.9 3.1(或3)x 1 2 y 1 22)以一定速度施加位移边界格式1:apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step 3)添加边界条件格式:apply szz <数值> grad <梯度> range <围>示例1:apply szz -1e9 grad 0 0 8.3e5 range z 0 120示例2:apply szz -0.6e6 range z 0.05 0.15 group pile注:<数值>是梯度方向坐标0点的数值,可通过定义坐标围的上下值与梯度计算得到。

2. range(缩写:ran)通过range功能,可以使命令作用在一定指定围的目标上;如果一个命令没有使用range来确定围,则命令对整个模型有效。

1)利用坐标指定一定的围格式1:range z 0 1格式2:range z 2.9 3.1 x 1 2 y 1 22)利用分组来指定围格式:range group 13)以上两种的复合格式:range z -0.1 0.1 group pile4)利用id号来指定一定的围格式:range id 0 10该命令后跟起始id和结束id,这里的id可以是实体单元、网格、结构单元、接触面和节点的编号。

例:model elastic range id 1 10 ;指定id为1到10的单元为各向同性弹性本构。

5)命名一个围(需要先命名这个围)格式:range name <自己起一个名字> <围>示例:range name intersected_zones x 5 8 y 3 7引用这个围的格式:容+range nrange <自己起的那个名字>示例:range nrange intersected_zone6)命名一个柱形围格式:同4)示例1:range name bhole cylinder end1 0 0 0 end2 -9.14 0 0 rad 0.46 示例2:range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 &cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not7)any和notAny跟在range命令后表示“和”的意思,即使若干个定义的围同时起效;not跟在range后表示“除了”的意思,即除此之外的围示例1:group wall range group 1 any group 2 any ;把组1和组2定义为新组,组名为wall示例2:group wall range group 1 not ;把除了组1之外的所有部分定义为新组,组名为wall.示例3:range name interface1 cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.36 &cylinder end1 0 0 0 end2 0 0 -0.9 radius 0.34 not3. set用来设置重力、收敛标准、给某些变量赋值1)设置重力格式:set gravity 0 0 -102)设置收敛标准格式1:set mech ratio <value> :用R无量纲值确定收敛标准(默认10e-5)格式2:set mech force <value> ;绝对收敛标准,体系最大不平衡力小于某一临界值,此局限性较大。

以上两种格式也可以进行耦合:示例:set mech force 0 ratio 1e-4格式3:自定义收敛标准,比如某一关键点的位移或速度达到某一临界值(遇到这种情况进行补充)3)在FLAC3D中给FISH语言中的变量赋值格式:set 变量=value 如:set hh=-224)设置大变形格式:set large4)设置监测历史变量时的监测间隔格式1:set hist_rep <数目> ;默认的应该是10步一个间隔。

示例:set hist_rep 1格式2 hist nstep <数目>示例:hist nstep 54. history(缩写:hist)设置监测点要进行监测的项目,不设置编号的话会从1自动进行编号。

监测点必须用id或坐标这两种方法之一指定,新版本后需加add,不加add其实也可以。

1)监测整个模型的不平衡力格式:hist add id <编号>unbal2)设置节点监测项目格式1:hist add gp zdis 0 0 6(监测节点z方向位移)hist add gp zvel 0 0 6(监测节点z方向速度)格式2:hist add id <给定该监测变量的编号> gp zdis <坐标>或者监测点id示例1:hist add id 2 gp zdis 1 1 3示例2:hist add id 1 gp zdis id 1 ;监测1号节点z方向的位移,这个监测点编号是13)设置单元监测项目格式1:hist add zone szz 0 0 -1(监测节点zz方向应力)格式2:hist add id <给定该监测变量的编号> zone szz <坐标><坐标>或者监测单元id示例1:hist add id 7 zone szz 1.5 1.5 2.54)设置监测历史变量时的监测间隔格式1:set hist_rep <数目> ;默认的应该是10步一个间隔示例:set hist_rep 1格式2 hist nstep <数目>示例:hist nstep 55)将计算结果写入到txt文件夹中格式:hist write <监测变量的编号> skip <记录间隔> begin <监测开始的步数> end <监测结束步数> file 路径/文件名.txt示例1:hist write 7 vs 8 file C:\Users\zhangzhao\Desktop\7-8.txt示例2:hist write 7 file C:\Users\zhangzhao\Desktop\7.txt示例3:hist write 1 2 3 4 5 skip 20 begin 8300 end 8400 file C:\Users\zhangzhao\Desktop\7.txt ;把1-5号监测点在8300到8400步之间的数据,按照每20步一个的记录间隔输出到7.txt文件中5.model用来指定模型的本构模型,以下格式中mech可以省略格式1:model mech mohr ;定义全模型的本构模型格式2:model mech mohr range group <组名>;格式3:model mech null range group <组名>;格式4:model mech null range z 0 16.property(缩写:pro)给模型材料赋予参数格式:pro bulk <值> shear <值> coh <值> fric <值> ten <值> range <围>示例1:prop bulk 5.6420e6 shear 1.3789e6 coh 9.77e3 fric 15.36range group 1示例2:prop bulk 6e6 she 2.2e6 fric 10 coh 5000 range z 17.3 20.17.group用来进行模型的分组1)建模时直接将这部分命名为一个分组格式1:gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1格式2:gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 fill group 1 ;(对填充区域进行有别于外围材料的分组命名)格式3:gen zone radcylinder p0 0 0 -18 p1 8 0 -18 p2 0 0 -19.1 p3 0 8 -18 size 8 4 8 15 ratio 1 1 1 1.15 group 1 fill group 2 ;(外围区域命名为1,填充区域进行有别于外围材料的分组命名2)2)对一定围的数据命名分组格式:group 1 range < > ;用range定义围3)分开命令为一个组示例:group dam range x 5 7 z -5 0group dam range id 201 any id 211 any id 231 any id 241 anygroup dam range id 202 any id 212 any id 222 any id 232 any以上区域命名了一个分组dam。

可以不同的区域命名为同一个组名,最后就会是一个组了。

8.attach进行网格的连接,使用该命令时对各网格连接面上的单元尺寸有限制,要求它们之间的比率成整数倍格式:attach face range < > ;range后跟定义围的关键词9.reflcet(缩写:ref)以某对称轴生成镜像网格格式1:gen zone ref ;以z=0平面为对称轴镜像生成网格格式2:gen zone ref dip 90 ; 以y=0平面为对称轴镜像生成网格格式3:gen zone ref dip 90 dd 90 ; 以x=0平面为对称轴镜像生成网格格式4:gen zone ref dip 90 dd 90 range z 20 35 ; 可以加上围,在z 为20~35围以x=0平面为对称轴镜像生成网格。

相关文档
最新文档