FLAC3D 高手总结
- 格式:doc
- 大小:34.00 KB
- 文档页数:4
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 end1end2 radius r2 not,其中r22、如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
flac3d常用命令师姐整理。
pl sel rec sres my surfx 0 1 0 ra y 0 24 scale=0.05pl add fap ra cyl end1 0 0 -17 end2 0 42 -17 r 3.1pl add sel geom id on;查局部坐标pr sel node local ra id=1pr sel rec surfpl sel geom node off nodesys on scale=0.06pl sel geom node off surfsys on scale=0.06pl sel pile fo fx; 看桩轴力pl bl gr id on ra y 20 21;查看zoneIDpl sur gp on ra y 20 24 ;监测点位移,找点gpidpl sel geom cid on ra y 22 24 ;监测衬砌内力,找衬砌cidpl sel geom cid on node off ra y 21 23 ;监测桩号,找桩cidpl 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 attachpr 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。
FLAC3D网络高手总结教程FLAC3D是一种广泛应用于地下工程和岩土工程中的数值模拟软件,它能够模拟地下结构的应力、变形、流体流动等物理过程,以便工程师能够更好地理解和预测地下工程的行为。
在FLAC3D中,网络高手是指熟练掌握软件的使用技巧和模拟方法的用户。
本文将总结一些网络高手常用的技巧和教程,并提供一些在FLAC3D中使用3DEC的相关知识。
首先,一个网络高手需要了解FLAC3D的基本操作和模拟流程。
FLAC3D的模拟过程可以分为几个步骤:几何建模、材料参数定义、边界条件设定、数值求解和结果分析。
网络高手需要熟悉软件中各个功能的使用方法,并能够根据实际需求合理地选择参数和设定条件。
其次,网络高手需要学会使用脚本语言来批量处理模型和自动完成重复性工作。
FLAC3D支持FISH和C++两种脚本语言,网络高手需要熟悉这些脚本语言的语法和常用函数,并能够灵活运用它们来自动化地进行模拟和数据处理。
同时,网络高手还需要了解岩土材料的本构模型和参数的选择方法。
FLAC3D提供了多种本构模型,如弹性模型、Mohr-Coulomb模型、Hoek-Brown模型等。
网络高手需要在根据实际工程条件选择合适的本构模型,并合理确定其参数。
此外,网络高手还需要熟悉FLAC3D中的边界条件和荷载设定方法。
边界条件和荷载是模拟中重要的控制因素,网络高手需要根据实际工程情况设置合适的边界条件和荷载,以保证模拟结果的准确性。
最后,作为一个网络高手,了解FLAC3D中的一些常见技巧和优化方法也是必要的。
例如,如何合理地划分网格、如何选择合适的求解器和收敛准则、如何有效地进行后处理等。
这些技巧能够帮助网络高手更好地利用FLAC3D进行模拟,提高工作效率和模拟结果的准确性。
对于使用FLAC3D的网络高手来说,掌握上述技巧和教程是非常重要的。
通过合理地运用这些技巧和方法,网络高手可以更好地解决实际工程中的问题,提高工作质量和效率。
另外,对于一些特殊问题,FLAC3D提供了3DEC作为补充工具,网络高手也可以学习相关的使用方法和技巧,以扩展自己的应用能力。
Flac3D常用命令及常规问题总结Flac 3D 常用命令总结new 不退出flac3d,新建一个新的分析计算任务generate 简写为ge 产生网格zone 简写为zo 指示为3维网格体brick 简写为b 指示为长方形的网格体size x y z 指示长方形网格体在x、y和z轴所划分的网格数plot 显示,输入之后回车会出现plot base/0 可以输入plot的任何子命令(add添加条目(surface 网格体颜色、axes 坐标系颜色)等),若回到flac3d命令状态,则按一次回车或输入qiut 按x ,y,z 小写字母可以旋转当前视图中的X,Y ,Z 轴,相应的大写字母则方向相反。
m 放大,M缩小。
Ctrl+G可以切换彩色图和灰色图。
Model 简写为mo 定义材料模型Mohr 简写为moh 指摩尔-库伦模型Property 简写为pro 定义模型的材料参数,材料参数包括:密度density 简写为de,体积模量bulk简写为bu,切变模量shear 简写为sh,内摩擦角friction 简写为fric,内聚力cohesion 简写为c 抗拉强度tension 简写为ten等Set 用于设置模拟条件或控制条件,如重力加速度gravity,mechanical 设置静态力学分析的有关参数简写为mec,force 最大不平衡力简写为fo等。
比如set mechanical force 50Initial 简写为in 初始化网格的相关值。
Fix 保持网格节点指定参数(速度、压力、温度等)的值不变。
Range 指定范围History 简写为h 采样或记录迭代时模型中变量的值,检验采样记录时可以在plot子命令中输入history 1 、history 2等,忘记编号时可以输入print history就可以看到什么记录号对应什么变量或参数Nstep 简写为n 迭代次数,默认值等于10,ubalance 简写为unb 最大不平衡力,gp 采样网格节点的有关值,zdisplacement简写为zdis z轴方向的位移Solve 简写为so 控制相关过程的自动时间步一般等值线图(包括位移、速度、应力、压力、温度等),可以在一个新的视图中显示,首先要在plot子命令中新创建一个视图,可以输入“create 视图名”,回车,继续输入add contour/bcontour disp/szz…(contour为在当前视图中显示等值线图简写为con ,bcontour 绘制指定区的等值线图简写为bcon,disp 指定为位移量等值线图,szz 指定为应力等值线图)任意剖面上的等值线图,可以在刚创建的视图中再创建一个视图然后输入剖面参数,比如“set plane dip 90 dd 0 origin 3 4 0”(plane 简写为p 设置一个剖平面,dip 为剖面倾角,dd 为剖面的倾向,origin 简写为o 为剖面中的一点)Boundary 简写为bo 在视图中增加面的边界线框Behind 简写为be 为当前视图后面要转换视图则在plot子命令中输入“current 视图名”、show开挖把网格体的材料模型设置成空(null)模型即可,命令如下:Model null range x y zSet large 把材料设置成大变形求解时计算不收敛则用step命令,step 步数节点的ID号和重心的全局坐标可以由下面的命令得到:Print gp position单元体的ID号和重心的全局坐标可由如下的命令得到:Print zone所有命令大小写一样,都可以附带若干个关键字和有关的值,格式如下:Command 〈keyword value…〉[keyword,value…]…命令、关键字、数字可以用数个空格隔开,也可以用(),=隔开,;为注释符,;开始到行末都是注释一个简单的命令行只能至多80个字符(包括注释),如果超过80个字符则在行末输入&,每个命令最多含1024个字符。
学习flac3d的心得:经过几天的学习,我学会了很多关于flac3d的基础知识,我先从基本的符号入手,知道了很多词语的意思,如xdis-x的位移,xvel-x的速度,fri-摩擦角,bu sh coh -代表不同的模量,还有很多基本的词语.然后从最简单的如何生成网格学习,学会了简单网格的生成,如怎么生成立方体(长,宽,高为95 85 75),步鄹如下:gen zone bri p0 0 0 0 p1 95 0 0 p2 0 85 0 p3 0 0 75 p4 95 85 0 p5 0 85 75 p6 95 0 75 p7 95 85 75 size 2 3 4 ratio 1.2 1.2 3 如何在flac3d的软件显示立方体只需要在命令行输入plot block group 这样的命令.还有可以在图形上显示射线的,步鄹如下:gen zone rad 下面的如上所述,但是我对球体的形成还是不懂,再者我学习了解决实际工程问题的步鄹,第一,根据题目中的条件在纸上按比例画出图形,第二,计算出初始的应力平衡,第三,根据画出的图形在flac3d生成,第四,对生成的图形进行赋值还有定义dens ,grav等已知量,第五,进行数据的处理,得出想要的结果.最后我学习对数据的后处理,可以用图形得到各层的受力的情况,还可以得到很多对实际工程有用的数据.如输出xx轴的受力输入命令:plot con sxx(不显示网格) plot con sxx ou on(显示网格).下面例题:NewGen zon bri size 2 2 4Model mohrProp bu 3e7 sh 1e7 coh 1e10 fri 15Fix z ran z 0Fix z ran z 2Fix x ran x 0Fix x ran x 1Fix y ran y 0Fix y ran y 4Ini dens 3000Set grav 9.8Ini sxx -3e4 grad 3.7777e3 0 0 ran z 0 1 Ini syy -4e5 grad 4.6666e4 0 0 ran z 0 1 Ini szz -6e7 grad 5.8888e6 0 0 ran z 0 1 App nstress -4e7 ran z 1Hist id=1 gp zdis 0 0 0Hist id=2 gp zdis 1 0 0Hist id=3 gp xdis 0.5 0 0Hist id=4 gp ydis 2.5 0 0Hist id=5 zone szz 0.5 2.5 1Hist id=6 zone syy 0.5 1 1Hist id=7 zone sxx 0.5 3 1Save 1-1.savSolve。
目录1.FLAC3D的固流耦合计算模式 ------------------------------------------------------------------------------------ 1 2.FLAC3D固流耦合学习小结 --------------------------------------------------------------------------------------- 5 3.关于流固耦合的问题 ----------------------------------------------------------------------------------------------- 6 4.也谈采用FLAC3D对地下采矿的模拟 ------------------------------------------------------------------------- 8 5.FLAC3D本构模型开发 ---------------------------------------------------------------------------------------------- 8 6.FLAC3D自定义本构模型----------------------------------------------------------------------------------------- 11 7.数值计算中初始应力场的模拟-------------------------------------------------------------------------------- 13 8.FLAC3D应变分析--------------------------------------------------------------------------------------------------- 13 9.FLAC3D的调参 ------------------------------------------------------------------------------------------------------ 14 10.开采沉陷垂直剖面等值线的生成--------------------------------------------------------------------------- 15 11.FLAC3D的应变硬化软化模型 -------------------------------------------------------------------------------- 16 12.FLAC3D的塑性流动格式 --------------------------------------------------------------------------------------- 17 13.FLAC3D的动画制作---------------------------------------------------------------------------------------------- 17 14.地下连续墙基坑开挖支护 ------------------------------------------------------------------------------------- 18 15.一个汇的小例子 -------------------------------------------------------------------------------------------------- 21 16.用3DEC生成岩体随机节理网络 ---------------------------------------------------------------------------- 23 17.固结小算例 --------------------------------------------------------------------------------------------------------- 24 FLAC3D的固流耦合计算模式英文原文-------------------------------------------------------------------------- 261.FLAC3D的固流耦合计算模式/blog/static/323428402007102243049387/FLAC3D的计算模式中是否需要做孔压分析取决于是否采用config fluid命令。
FLAC3D学习笔记(自己总结版)一、建立模型笔记1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转、如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转、+号的功能是加大模型旋转的角度,-号的功能是减小模型旋转的角度。
2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动、3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:gen zon tunint p0 0 0 0 p1 010 0 p255 0 p3 0 05 p6 0105 dim sq2 sq21 sq2111 group(编组)z1(组名为:z1)&fill group tunnel4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点、其中一个是原点(origin(xyz))、5镜像部分模型的方法: gen zone reflect normal (x,y,z)origin (x,y,z)range group z1就可以了、6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像、7 group radcyl range group cylint not (编组的时候不包括cylint)7 快捷键总结Ctrl+Z 局部缩放的功能。
Ctrl+R的功能是恢复到初始状态。
8模型建立的方法:方法1:利用默认值生成网格,各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值、方法2:利用4个点的坐标来生成矩形网格、p0~p3 size方法3:利用edge来确定边长生成矩形网格、方法4:利用参数ratio来确定单元体几何边画边生成矩形网格、方法5:利用参数add(相对坐标)来生成矩形网格、方法6:利用八个点的左边来生成矩形网格、9、用户自定义模型的方法我来贡献一点自己的成果:; U" j/ ?7 L+ a0 uFLAC3D的二次开发环境提供了开放的用户接口,在软件安装文件中包含了软件自带所有本构模型的源代码,且给出了Mohr-Coulomb模型和应变软化模型的编译示例,因此可以方便地进行本构模型的修改与开发。
一、建立模型笔记1键盘上的键的功能:X,Y,Z分别控制所建立的模型围绕X,Y,Z三个轴做逆时针旋转.如果打开大小写锁定键,分边按XYZ时建立的模型会做顺时针旋转.+号的功能是加大模型旋转的角度,-号的功能是减小模型旋转的角度。
2键盘上的上下左右四个键可以实现图形的向上向下,向左和向右的移动.3当在一个区域建立多个模型的时候,有时候为了镜像单个模块儿,需要预先编组,编组的方法如下:gen zon tunint p0 0 0 0 p1 0 10 0 p2 5 5 0 p3 0 0 -5 p4 5 10 0 &p5 5 5 -5 p6 0 10 -5 p7 5 10 -5 dim sq2 sq2 1 sq2 1 1 1 group(编组) z1(组名为:z1) &fill group tunnel4镜像整个模型的方法:首先判断即将镜像的模型在原模型的哪个方向(即三个轴的方向),然后在轴所在的方向上去两个点.其中一个是原点(origin(xyz)).5镜像部分模型的方法: gen zone reflect normal (x,y,z) origin (x,y,z) range group z1就可以了.6特别要注意的一点:交叉巷道的巷道充填和巷道的网格是两个组,所以在镜像的时候不要忘了给充填部分镜像.7 group radcyl range group cylint not(编组的时候不包括cylint)7 快捷键总结Ctrl+Z 局部缩放的功能。
Ctrl+R的功能是恢复到初始状态。
8模型建立的方法:方法1:利用默认值生成网格,各个默认值:(1)p0的默认值是(0,0,0),网格的每边的单元格数size默认值是10,网格每边的长度默认值为size的值.方法2:利用4个点的坐标来生成矩形网格.p0~p3 size方法3:利用edge来确定边长生成矩形网格.方法4:利用参数ratio来确定单元体几何边画边生成矩形网格.方法5:利用参数add(相对坐标)来生成矩形网格.方法6:利用八个点的左边来生成矩形网格.9、用户自定义模型的方法我来贡献一点自己的成果:FLAC3D的二次开发环境提供了开放的用户接口,在软件安装文件中包含了软件自带所有本构模型的源代码,且给出了Mohr-Coulomb模型和应变软化模型的编译示例,因此可以方便地进行本构模型的修改与开发。
;模型镜像
gen zone radcylinder size 25 1 25 25
gen 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,0
hist zone syy 0,1,0
hist zone syy 1,1,0
pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror"
;对模型进行压缩实验的方法
;即在模型两侧施加相反方向的速度
ini yvel 1e-7 range y -.1 .1
ini 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*degrad
def set_vals
ptA = 25.0 * sin(pi/2);ptA=25.0
ptB = 25.0 * cos( 60.0*degrad );ptB=12.5
ptC=pi;ptC=3.1415926
end
set_vals
print 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 not
plot add sel geom black black cid on scale=0.03
sel node init zpos add -25.0
;如何显示某一平面
plot create name_plane
plot set plane origin 3 4 0 normal 1 0 0
plot add cont disp plane behind shade on
plot add sel geom black
plot add axes red
plot current name_plane
plot 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 aa
get_id=gp_id(gp_near(3,0,0))
end
aa
print get_id
print gp position range id 3
;保存命令流文件
set log on
set logfile hua.log;文件存放的路径与调用的文本路径相同
......
set log off
;fish函数的说明与解释
def find_add ;定义fish函数find_add
head = null ;给head赋值
p_gp = gp_head ;第一个网格结点的指针赋给p_gp
loop_while p_gp # null ;当p_gp值不为null时作循环
x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_pos
if 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值赋给head
endif
p_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地址链条,方便以后使用。
;FLAC后处理
;输出JPG图片
set plot jpg
set plot quality 100
plot hard file XX.jpg
;输出TXT文件
hist write X file xx.txt ;(x是历史记录ID)
print con zd
;输出单元应力和节点位移
set log on
set logfile xx.log
print zo stress
print gp dis
set log off
;生成影片文件
plot set rot 20 0 30
plot con szz ou on magf 10
plot add hist 1
set movie avi step 1 file xx.avi
movie start
solve
;计算间隔命令在solve命令之前
set hist_rep 1 ;(步数)
;检查模型赋值情况
plot bcon prop bu
;变形云图加网格放大20 倍
plot con zd ou on magf 20
;变形网格
plot sk magf 20
;变形矢量图(加网格)
plot sk dis scale 0.07 ;(scale是箭头长度,默认0.05)
;如何观察塑性区
plot block state
plot block state sh_n
plot block state sh_p
;切片操作
plot set plane ori 1 2 3 norm 0 0 -1
plot con zd plane
plot add dis plane
plot add sk
plot add ax
;切片后方部分显示
plot con zd plane behind shade on
交集就是range group a group b
并集就是range group a any group b any not是差集。