FLAC3D 高手总结
flac3D部分技巧
cohesion 内聚力
tension 抗拉强度
fric 内摩擦角
一个例子里的: prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
pl con zdis 位移、沉降
pl con szz zz应力 stress
pl hist 1 1、2等顺序是之前定义的先后顺序
pl hist 2 3 同时显示
上述命令后可加上 outline on(简写outl on) 表示显示网格线
ห้องสมุดไป่ตู้
几个模量 : bulk 体积模量K=E/3(1-2v) v是泊松比,当v接近0.5是不能盲目使用
2、dip 为剖面的倾角,dd为剖面的倾向!
3、经过几次固定坐标测试,终于搞清楚了,dip是指倾角,即剖切面和xoy平面的锐夹角,dd是倾向,即剖平面倾斜方向在XOY平面投影和Y轴正方向的夹角。以Y的顺时针为正。
总结:1、dip:xoy平面绕ox轴,顺时针(从x正方向看)旋转xx度
2、dd: xoy平面绕oz轴,顺时针(从z正方向看)旋转xx度
3、ori:经以上两次旋转后,将ox轴移动到 坐标(x,y,z)处,y一般为0.
ps:其实先进行第二条选择,再进行第一项旋转更好理解。
call 的文件名字里不能有空格,可以是中文
默认的global ratio limit 是 1e-5。计算过程应该是波浪形的 起伏。
查看结果:
pl bl gr 显示模型
dip和dd:我通过陈育民的那本书了解到一点不知道对不,不成熟的还请指教
FLAC3D常用命令总结
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常用命令师姐整理。
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非常重要的总结的程序语句
plot con szz ou on magf 20 ;变形后网格与应力云图
plot set plane ori x y z (切面过的点的坐标) norm x y z(切面的法向量)
plot set plane ori 2 2 0 norm 0 1 0 做切面
hist write 1 file lishi.txt
plot sketch;绘出轮廓线
plot add vel;添加矢量
plot add his 1 2 vs 3;绘出历史记录 n,actLoad , p_solcyc
15000;运行 15000 时步
plot add interface red 显示断层面红色
plot con zd plane (剖面切面Z向变形云图)
plot block state (显示塑性分布区)
plot block group 显示模型
plot block group id on (显示ID)
使用plo blo gro gpnum on 命令显示坐标
plot add surface
;显示表面
interface i effctive =on/off 接触面上有效应力打开或关闭。默认为on
interface i perm=on/off 打开或关闭渗流,默认on。只用于CONFIG fluid模式下,若perm打开,则不要设置maxedge,否则接触面就自动变为不透水。
print gp position 节点的ID号和全局坐标
prinlane behind shade on 显示切片后部模型
FLAC3D网络高手总结教程
FLAC3D网络高手总结教程FLAC3D是一种广泛应用于地下工程和岩土工程中的数值模拟软件,它能够模拟地下结构的应力、变形、流体流动等物理过程,以便工程师能够更好地理解和预测地下工程的行为。
在FLAC3D中,网络高手是指熟练掌握软件的使用技巧和模拟方法的用户。
本文将总结一些网络高手常用的技巧和教程,并提供一些在FLAC3D中使用3DEC的相关知识。
首先,一个网络高手需要了解FLAC3D的基本操作和模拟流程。
FLAC3D的模拟过程可以分为几个步骤:几何建模、材料参数定义、边界条件设定、数值求解和结果分析。
网络高手需要熟悉软件中各个功能的使用方法,并能够根据实际需求合理地选择参数和设定条件。
其次,网络高手需要学会使用脚本语言来批量处理模型和自动完成重复性工作。
FLAC3D支持FISH和C++两种脚本语言,网络高手需要熟悉这些脚本语言的语法和常用函数,并能够灵活运用它们来自动化地进行模拟和数据处理。
同时,网络高手还需要了解岩土材料的本构模型和参数的选择方法。
FLAC3D提供了多种本构模型,如弹性模型、Mohr-Coulomb模型、Hoek-Brown模型等。
网络高手需要在根据实际工程条件选择合适的本构模型,并合理确定其参数。
此外,网络高手还需要熟悉FLAC3D中的边界条件和荷载设定方法。
边界条件和荷载是模拟中重要的控制因素,网络高手需要根据实际工程情况设置合适的边界条件和荷载,以保证模拟结果的准确性。
最后,作为一个网络高手,了解FLAC3D中的一些常见技巧和优化方法也是必要的。
例如,如何合理地划分网格、如何选择合适的求解器和收敛准则、如何有效地进行后处理等。
这些技巧能够帮助网络高手更好地利用FLAC3D进行模拟,提高工作效率和模拟结果的准确性。
对于使用FLAC3D的网络高手来说,掌握上述技巧和教程是非常重要的。
通过合理地运用这些技巧和方法,网络高手可以更好地解决实际工程中的问题,提高工作质量和效率。
另外,对于一些特殊问题,FLAC3D提供了3DEC作为补充工具,网络高手也可以学习相关的使用方法和技巧,以扩展自己的应用能力。
Flac3D常用命令及常规问题总结
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个字符。
学习Flac心得
学习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。
FLAC3D网络高手总结教程(含部分3DEC)
目录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学习笔记(自己总结版)
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模型和应变软化模型的编译示例,因此可以方便地进行本构模型的修改与开发。
FLAC3D学习笔记(自己总结版)
一、建立模型笔记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模型和应变软化模型的编译示例,因此可以方便地进行本构模型的修改与开发。
(总结452类)FLAC3D命令流挺实用
、怎样查看模型?答:可以查看网格,可以查看节点号。
、请问在圆柱体四周如何施加约束条件?答:可以用... ,其中<,其实就是选择的时候选两个圆柱的差,即得到边界。
命令流如下:、怎么能把一个的图像数据导出来以便用其他软件绘图?答:用命令,把数据导出来,转到里处理一下,然后用或者什么作图软件绘制就行了。
、用命令建立模型后,如何显示点的坐标?答:使用命令、关于对性质进行高斯正态分布的问题?答:根据手册上的说明:下面的命令设定一个平均摩擦角为度,标准方差是±%。
则命令如下:问题:请问中的是如何计算的?如果把±%改为±%,则命令应如何写?×%=、问题问:() () () ()(对面做镜像)(对面做镜像)()表示轴正向顺时针到那条射线的夹角,表示对称参照面与平面的夹角,对称参照面与平面的夹角在平面的投影是一条射线。
首先应该按照的方向大体确定这个面的朝向,指的是从轴正方向按顺时针()方向转向所要确定面的法线方向在平面上的投影的夹角,然后再确定,指的是从平面转向所要确定的平面的角.(及轴负方向转过角度)()命令:(根据法线和过一点建的面)最常用,至于实际使用过程中有人问道:与的区别,我试过两者的效果是一样的,没区别,虽然方向不同,但表达的是同一个面。
这个就是用里建立斜面镜像的。
有一点注意了,镜像的时候如果是模型镜像的话,相当于在建立一个模型,故,在的前面要加建模的关键字,……7、关于两种是否有区别?和答:是一样的,表示并列连接关系,及的意思。
8、的区别答:初始化命令,初始化计算应力状态,定义材料参数、重力等。
边界条件限制命令,如施加边界的力、位移等约束等。
用于设置计算的模式,控制计算过程。
施加的初始力在计算过程中是会变的,一般体力需要初始化,而施加的力、边界条件不会随着计算变化。
命令作用于整个网格(赋值全局变量)则初始化指定范围内网格点、域的偏移量。
初始应力不仅仅是重力引起的。
FLAC3D 软件使用心得体会 张强
FLAC3D 软件使用心得体会张强 20104323 10级水利3班一,初识FLAVC3D软件:1,作用:它是一款数值计算软件,有二维和三维计算程序两个版本,在岩土工程领域用处很大,其特点是简单易懂,能够进行土质.岩石和其他材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
2,主要应用:(1) 承受荷载能力与变形分析: 用于边坡稳定和基础设计(2) 渐进破坏与坍塌反演: 用于硬岩采矿和隧道设计(3) 断层构造的影响研究: 用于采矿设计(4) 施加于地质体锚索支护所提供的支护力研究: 岩锚和土钉的设计(5) 排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究: 挡土墙结构的地下水流动, 和土体固结研究(6) 粘性材料的蠕变特性: 用于碳酸钾盐矿设计(7) 陡滑面地质结构的动态加载: 用于地震工程和矿山岩爆研究(8) 爆炸荷载和振动的动态响应: 用于隧道开挖和采矿活动(9) 结构的地震感应: 用于土坝设计(10) 由于温度诱发荷载所导致的变形和结构的不稳定(11) 大变形材料分析: 用于研究粮仓谷物流动和放矿的矿石流动二.操作对象:Flac3D中为岩土工程问题的求解开发了特有的本构模型, 总共包含了10种材料模型:(1) 开挖模型null(2) 3个弹性模型(各向同性, 横观各向同性和正交各向同性弹性模型)(3) 6个塑性模型(Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型).Flac3D网格中的每个区域可以给以不同的材料模型, 并且还允许指定材料参数的统计分布和变化梯度. 还包含了节理单元, 也称为界面单元, 能够模拟两种或多种材料界面不同材料性质的间断特性. 节理允许发生滑动或分离, 因此可以用来模拟岩体中的断层、节理或摩擦边界.FLAC3D中的网格生成器gen, 通过匹配、连接由网格生成器生成局部网格, 能够方便地生成所需要的三维结构网格. 还可以自动产生交岔结构网格(比如说相交的巷道), 三维网格由整体坐标系x, y, z系统所确定, 这就提供了比较灵活的产生和定义三维空间参数.三.操作命令:结合一些边界条件加上各个实际情况就可以使用该软件对岩石的力学等性质进行计算了。
FLAC3D学习经验积累
FLAC3D学习经验积累Flac 3D 学习经验积累(by seapie)用FLAC3D解决问题时,为了得到最有效的分析使模型最优化是很重要的。
通过资料收集和作者本人经验,对改进模型的运行提供了一些方法建议。
一、学习经验积累1、模型建立1.1 考虑网格划分的密度FLAC3D使用常应变单元。
如果应力/应变曲线倾斜度比较高,那么你将需要许多区域来代表多变的分区。
通过运行划分密度不同的同一个问题来检查影响。
FLAC3D应用常应变区域,因为当用多的少节点单元与用比较少的多节点单元模拟塑性流动时相比更准确。
(见理论卷第二章和实例卷中的塑性例子)。
应尽可能保持网格,尤其是重要区域网格的统一。
避免长细比大于5:1的细长单元,并避免单元尺寸跳跃式变化(即应使用平滑的网格)。
应用GENERA TE命令中的比率关键词,使细划分区域平滑过渡到粗划分区域。
2、本构模型与材料属性2.1 改变模型材料FLAC3D 对一个模拟中所用的材料数没有限制。
这个准则已经尺寸化,允许用户在自己所用版本的FLAC3D中最大尺寸网格的每个区域(假如设定的)使用不同的材料。
3、边界条件3.1 初始应力场的加载岩土工程除常应力场外一般均应加重力场set grav 0 0 -10.0,并应加单元的密度ini density dens1(kg/m3)。
若还考虑构造应力场时,两种情况:1 常构造应力,命令用apply,如apply sxx -1.0e52 递度构造应力,命令用ini,如ini sxx -80 grad 0 0 8.0,注意gz=8.0μ-的倍数关系加水平应力场,所以在通注意:在加载SZZ或重力场时,程序自动按)1(μ常的加载水平应力场时应考虑这部分的值。
3.2 边界约束初始化变量在模拟基坑开挖过程时,在达到目的前通常要初始化网格顶点位移。
因为计算次序法则不要求位移,所以可以初始化位移,这只是由网格顶点的速度决定,并有益于用户初始化速度却是一件难事。
FLAC3D网络高手总结教程
初始条件优化
通过试算和对比分析,对初始条件进行调整和优化,以减小仿真 误差。
求解过程监控与结果
求解过程监控
01
在求解过程中,实时监控模型的计算状态,如收敛情况、计算
时间等。
结果输出与Байду номын сангаас视化
02
将仿真结果以云图、等值线、矢量图等形式进行可视化输出,
了解最新研究进展和应用案例。
利用在线资源和学习平台
03
推荐学员们利用各类在线资源和学习平台,如教程、
案例库、技术论坛等,持续学习和提升自己的技能。
THANKS
感谢观看
大变形问题处理
采用FLAC3D的大变形分析技术,对大变形问题进 行高效求解,如滑坡、泥石流等灾害过程的模拟。
非线性迭代算法优化
通过优化FLAC3D的非线性迭代算法,提高 非线性问题求解的收敛性和计算效率。
并行计算加速技术应用
01
并行计算环境搭建
介绍如何在FLAC3D中搭建并行 计算环境,包括硬件配置和软件 设置等。
结果输出异常或不符合预期原因排查
异常一
输出结果与预期不符,如何排查?
排查
首先检查模型设置和输入参数是否正确。然后,尝试使用不同的算法或调整参数进行计算,观察结果是否有变化 。此外,也可以查看详细的计算过程和日志信息,以便更深入地了解问题所在。
结果输出异常或不符合预期原因排查
异常二
结果中出现明显的数值不稳定现象,如何排查?
经验分享
在案例分析和实战演练过程中,积累了一些宝贵的经验和教训。例如,在建模过程中要 注重细节和质量控制;在选择材料参数时要结合实际情况进行综合考虑;在设置边界条 件和荷载时要确保准确性和合理性等。这些经验和教训对于提高FLAC3D建模水平和解
flac自己总结小命令集合
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 end2radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
网上FLAC3D学习笔记2
网上FLAC3D学习笔记21 请问在圆柱体四周如何施加约束条件?可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
命令流如下:newgen zone cyl p1 1 0 0 p2 0 4 0 p3 0 0 1 size 12 30 12gen zone refle norm 1 0 0gen zone refle norm 0 0 1mo ssprop den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10 prop ftab 1 ctab 2 dtab 3table 1 0 45 0.05 42 0.1 40 1 40table 2 0 2e6 0.05 1e6 0.1 5e5 1 5e5table 3 0 10 0.05 3 0.1 0fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5回答:加上的约束不对!fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 0.5 0 0 end2 0.5 4 0 rad 0.5回答:上面的rad 是不是相差太大了?是不是在1左右来限制了2.关于计算不收敛问题;初始网格gen zone radcyl p0 15 0 -15.75 p1 add 15 0 0 p2 add 0 6 0 p3 add 0 0 15.75 &size 3 4 8 5 dim 2.25 2.25 2.25 2.25 rat 1 1 1 1.1 fillgen zone radtun size 3 4 3 5 dim 0.75 2.25 0.75 2.25 rat 1 1 1 1.1 & p0 15 0 -15.75 p1 add 0 0 -14.25 p2 add 0 6 0 p3 add 15 0 0 fillgen zone reflect dip 90 dd 270 ori 15 0 -15.75; 摩尔库仑模型model mohrmacro niyan 'bulk 3.2825e3 shear 1.5150e3 fric 20 coh 0.07' macro meiceng 'bulk 2.1735e3 shear 8.8917e2 fric 19 coh 0.05'prop niyan range z =-13.5 0prop meiceng range z =-15.5 -13.5prop niyan range z = -30 -15.5;为开挖定义整体部分range name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25range name hangdao x 12.75 17.25 y 0 6 z -16.50 -15.75;进行加载fix x y z range z -30.001 -29.999app sxx=-20 range x -0.001 0.001app sxx=-20 range x 29.999 30.001app syy=-20 range y -0.001 0.001app syy=-20 range y 5.999 6.001app szz=-20 range z -0.001 0.001solve(1)model null range inmodel null range hangdaorange name in z -15.75 0 cyl end1 15 0 -15.75 end2 15 6 -15.75 rad 2.25range name hangdaoliangce1 x 12.749 12.751 y 0 6 z -16.50 -15.75range name hangdaoliangce2 x 17.249 17.251 y 0 6 z -16.50 -15.75range name hangdao x 12.75 17.25 y 0 6 z -16.501 -16.499 apply nstress -0.05 plane dip 90 dd 0 range in; 给巷道底部加压力app nstress -0.05 plane dip 90 dd 0 range hangdao; 给巷道直墙部分(左侧)加压力apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce2 ; 给巷道直墙部分(右侧)加压力apply nstress -0.05 plane dip 90 dd 0 range hangdaoliangce1;对点进行监测hist unbalhist gp xdisp 12.75 0 -15.75:监测左帮hist gp zdisp 15 0 -13.5 ;监测顶板hist gp zdisp 15 0 -16.5 ;监测底板solve问:第一个solve平衡了,是第二个solve不能平衡,这是为什么?如果我把支护力加大,也不能平衡回答:我遇到的不收敛的原因,有两个:1)检查你的模型是不是没有建好,相互之间是不是没有连接好;2)仔细检查你的参数,是不是已经发生了破坏。
flac3d命令总结
6.加载前的位移清零
ini xdisp=0 ydisp=0 zdisp=0
7.开挖命令
m null range x 20 100 y 20 30 z 20 24.50
8.看切面及记录数据
pl set plane ori 0 0 76.85 norm 0 0 1
12.显示塑性区
plot block state
13.显示部分网格情况(位移上的范围)
plot block state/zadisp/szz range z 3.15 7.65
14.显示部分网格情况(id上的范围)
pl block state/zadisp/szz range id 501 1100
21.出图命令
set plot jpg
set plot quality 40
plot hard file name
22.看已知单元id号岩性参数方法
菜单栏-display-set range-"add"-ID-"在low里输入ID号"-确定-再点display-zones-property-bulk
15.修改间隔
pl con szz interval 2.5e06
16.显示网格
pl grid
17.查看节点号
pl grid num
18.查看模型有多少个单元、节点。
pl info
19.查看边界条件
plot fix
20.如何保存/名
1.建模型命令
new
plot block group
gen zone brick size 60 100 2 p0 0 0 0 p1 120 0 0 p2 0 200 0 p3 0 0 8.40 group 1
Flacd学习体会
Flac3d学习体会Flac3d是有限差分数值分析软件,它能够模拟岩土和其它材料的三维力学特性。
因此,在岩土工程力学中取得了普遍的应用,咱们能够用该程序较好的模拟地质材料在达到强度极限或屈服极限时发生的破坏或塑性流动的力学特性,专门适用于分析渐进破坏失稳和模拟大变形。
该软件中提供的壮大丰硕的本构模型能够进行静力、动力、蠕变、渗流、温度等多种计算模式,各类模式间能够相互耦合,也能够模拟多种结构形式,如梁、锚杆、土工格栅、衬砌、桩等。
作为一名初学者,要紧了解了下flac3d中的弹性模型和莫尔库仑模型,通过一段时刻的练习要紧有以下心得;在分析工程案列时,咱们能够依照flac3d的一样求解流程进行剖析案列,将命令流归纳为三个部份,即分析模型、模拟及求解、计算结果输出。
分析模型部份包括生成网格单元、设置初始条件和边界条件和初始应力平稳;在模拟求解部份包括加载及持续建模、后续工程施工、场方程的有限差分求解;输出计算结果部份要紧为图表的绘制、相关数据的输出等;通过这三个部份,咱们能够写出一个完整的命令流,方便后续在运行进程中修改命令流和查出具体在哪个环节犯错了。
那个地址咱们应该注意在flac3d的成立分析模型部份,材料性质的概念、初始条件和边界条件的设置并无明显的前后顺序。
关于工程实例中,咱们在分析模型中,第一看是不是对称,若是对称咱们能够简化模型,然后用gen zone reflect dip 关键字dd 关键字ori 关键字那个命令来镜像模型。
那个地址咱们需要注意的确实是dip和dd的具体含义才能把模型建好;dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线,dd表示y轴正向顺时针到那条射线的夹角。
第一应该依照dd的方向大体确信那个面的朝向,dd 指的是从y轴正方向按顺时针(clockwise)方向转向所要确信面的法线方向在xy平面上的投影的夹角,然后再确信dip,dip指的是从xy平面转向所要确信的平面的角.gen zone reflect dip 90 dd 90是以YZ面.gen zone reflect dip 90 dd 0 才是以xz为对称面的当初始应力平稳后,咱们能够求工程的阻碍,求出区域内的应力散布情形后,再依工程的流程及步骤时期执行各工程进行进程的阻碍,建议利用以下的步骤:(1)调出起初的应力平衡:re_____ .sav(2)设定新的材料性质:model,prop(3)设定新的支撑性质:struct(4)设定新的外力 apply(5)设定边界条件 apply,fix(6)求工程时的应力平衡 solve, step, cycle(7)储存 save在初始应力平稳进程中,咱们要注意是不是收敛,关于不收敛一样要检查模型是不是没有建好,彼此之间是不是没有连接好;认真检查参数,是不是已经发生了破坏;关于计算时刻太长,咱们能够适当修改模型的网格数以减少计算时刻,或使用step命令,那个地址的关键问题是要运行多少步。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
;模型镜像
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是差集。