[实用参考]Flac3d-5.0常用命令集锦.doc
flac3d基础相关命令
Flac3D基础知识1.Generate 命令产生网格,简写为 GE。
2. zone 参数代表了三维网格体(单元体),简写为 zo。
3. (1)size 参数指示单元体各边小单元体的个数;(2)dimension 参数指示单元体内部巷道各边边长,如果没给出就默认为两点间距的20%;(3) edge 参数指示网格体各边的尺寸 Ge zone brick p0(a,b,c) edge x y z 边长为x y z 顶点为p0;(边长也可以用确定点坐标来利用点与点的间距来确定p0 p1 p2……)。
4.各种网格体名称及关键字:矩形网格(Brick),退化矩形网格(dbrick),楔形网格(wedge),金字塔网格(pyramid),四面体形网格(tetrahedron),圆柱体网格(cylinder),矩形体外环绕放射状网格(radbrick),平行六面体外环绕放射状网格(radtunnel),圆柱体外环绕放射状网格(radcylinder),柱形壳体网格(Cshell),交叉圆柱体网格(cylint),交叉平行六面体网格(tunint)。
5.绘制网格体语句:①Generate zone brick size a,b,c (其中的brick可以被其他的网格体替换)。
②Generate zone brick p0=(,,) &P1=(,,) &P2=(,,) & (边长由点点间距确定)P3=(,,) & → ( &符号为继续行的符号)Size a b c → (各边小单元体个数)③如果在p1 p2 p3 前面加 add 则各点坐标是相对于p0点的坐标。
④如果只给出p0 则由 edge a & 给出边长。
(见第 3 条);利用ratio a,b,c 给出小单元体长度比例依次变化。
⑤内部有空单元体或者巷道的网格体,如果内部巷道需要填充为单元体,使用 fill group groupname 语句。
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常用命令总结第一篇: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常用命令
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常用命令
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-5.0常用命令集锦
建模1、调用文件:①文件与工程在同一个文件夹,只写文件名即可:Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlPthefilenamen eedbeenteredwiththe CALL command. ②不在同一个文件夹,全路径:Otherwise,thefilemaPbecalledbPspecifPingits completepath(e.g.,c:\mPfolder\file.dat ). Undo ;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr ;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z -direction.(垂直向下!常用的)下下面面这这代代码码,,是是沿沿着着--P P 方方向向的的重重力力加加速速度度,,注注意意区区别别!!!!!!!!genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续 p110,-10,-20... p2-10,10,-20... p3-10,-10,0 plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体 p110,-10,-20p2-10,10,-20... p3-10,-10,0p410,10,-20... p5-10,10,10p610,-10,0... p710,10,10plotcurrentplotPlot01 plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01P01z01;显示指定范围内各单元的主应力,结果如下Hist命令:①命令编号按顺序从1开始:eachhistorPisnumberedsequentiallP from1asitisenteredviathe HISTORP command.②查找显示所有的his命令:ReturntotheFlac3D>promptandtPpelisthistforalistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。
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 end2radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
FLAC3D5.0详细命令
1、plot zones and gridpoint velocities directly, simply typeplot zone velocity2、plot linecol or colorspecifies the color of the line.on/off turns line drawing for this option on or off.stylespecifies the linestyle used.width ispecifies the width of the line. This value must bebetween 0 and 10. A value of 0 indicates the smallestpossible line (1 pixel).3、plot textNot all of the following keywords are available in all contexts.col or colorspecifies which color the text will be drawn in.fa mily stringspecifies the name of the font family (for instance, “Ar-ial”or “Times New Roman”). (Note that underscoresmust be used in place of spaces in the name.)on/off turns text drawing on or off.4、plot create plot Myplotplot text on textplot text"this is my plot";(怎样显示文字?不会)5、保存图片plot create plot GravV;plot cut add plane dip=90 dd=0 origin=3,4,0plot cut add plane origin=119,90,110nor 0 1 0plot set orient 120, 220, 0plot set center 45, 45, 0plot add zone plane behindplane onplane offplot add bcont szzplot add axesplot post view GravV filename "103F1068"plot post filename "721F01F3-fft-input.ps" manualplot post view GravV filename "GravV.ps" manualplot post plot Plot01 file 127F1017.ps manual上面保存的是.ps的图片,不容易打开,我知道用photoshop可以打开,麻烦。
(新)Flac3D5.0笔记(注释)
FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流:Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后,如何再次呈现?输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面,双击I nterface显示接触面
接触面属性项被激活后,可以更改属性隐藏网格显示
输出网格
如果点击了透明,将不显示Information
只要不点击透明,则显示Information
等同于
点击显示节点坐标
fish简介
之前要有(空格)之前加是Flac3D3.0和5.0区别而在之间没有(空格)
按下列步骤取消fish语言前要加
点击fish
软件的刷新如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量(该点与原点的坐标差)
对于结构单元如桩,若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。
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个字符。
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 end2radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
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 end2radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
FLAC3D常见命令与使用技巧
FLAC3D常见命令与使用技巧1.基本命令:- !new: 创建一个新的FLAC3D模型- !zone generate: 通过定义边界和离散网格参数生成岩石模型- zone struct properties: 设置岩石模型的物理性质,如弹性模量、泊松比等- zone gridpoint initialize: 初始化模型中的网格节点- zone cmodel assign: 在指定区域内分配一种材料模型- zone face apply stress: 在模型表面施加一定的应力条件- zone solve: 求解模型,得到模型的应力应变分布2.材料模型:- zone cmodel initialize: 初始化材料模型的参数- zone cmodel property assign: 为材料模型分配属性- zone cmodel model assign: 为材料模型指定模型类型3.自定义模型:- zone cmodel new: 创建一个自定义材料模型- zone cmodel property new: 创建一个自定义的材料模型属性- zone cmodel property set: 设置材料模型属性的值- zone cmodel initialize: 初始化材料模型的参数4.辅助命令:- zone gridpoint fix: 固定指定区域的网格节点- zone gridpoint value prescribe: 为指定区域的网格节点预设值- zone gridpoint stress initialize: 初始化模型中的网格节点应力- zone gridpoint velocity assign: 为指定区域的网格节点分配速度- zone gridpoint delete: 删除指定区域的网格节点5.使用技巧:-使用脚本:可通过编写脚本,将一系列命令自动执行,提高工作效率-模型可视化:FLAC3D提供了强大的可视化功能,可以实时观察模型行为-使用参数化模型:可以通过设定参数的方式进行模型设置,便于修改和调整-多个窗口:可以同时打开多个FLAC3D窗口,方便对不同模型进行比较和分析-参考案例:可以参考软件自带的案例,了解不同类型的建模和分析方法这些是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 fillplot surfgen zone reflect dip 0 dd 90 (对于xy里搞镜像)gen zone reflect dip 90 dd 90 (对于yz里搞镜像)(1)dd表示y轴正背顺时针到那条射线的夹角,dip表示对于称参照里与xy仄里的夹角,对于称参照里与xy仄里的夹角正在xy仄里的投影是一条射线.最先该当依照dd的目标大概决定那个里的往背,dd 指的是从y轴正目标按顺时针(clockwise)目标转背所要决定里的法线目标正在xy仄里上的投影的夹角,而后再决定dip,dip指的是从xy仄里转背所要决定的仄里的角.(及z轴背目标转过角度)(2)下令:reflect normal xyz origin xyz (根据法线战过一面建的里)最时常使用,至于本量使用历程中有人问道:normal -1 0 0 与normal 1 0 0的辨别,尔试过二者的效验是一般的,出辨别,虽然目标分歧,但是表黑的是共一个里.reflect dd a dip b origin xyz 那个便是用里建坐斜里镜像的.有一面注意了,镜像的时间如果是模型镜像的话,相称于正在建坐一个模型,故,正在reflect的前里要加建模的关键字,gen zone reflect……7、关于二种fix是可有辨别?fix z range z -0.1 0.1fix x range x -0.1 0.1 any x 1.9 2.1 anyfix y range y -0.1 0.1 any y 1.9 2.1 any战fix z range z -0.1 0.1fix x range x -0.1 0.1fix x range 1.9 2.1fix y range y -0.1 0.1fix y range y 1.9 2.1问:是一般的,any表示并列对接关系,及and的意义.8、initial apply set 的辨别问:initial 初初化下令,初初化估计应力状态,定义资料参数、重力等;apply鸿沟条件节造下令,如施加鸿沟的力、位移等拘束等;set 用于树坐估计的模式,统造估计历程.Ini 施加的初初力正在估计历程中是会变的,普遍体力需要初初化,而apply施加的力、鸿沟条件不会随着估计变更.set下令效率于所有网格(赋值局部变量),ini则初初化指定范畴内网格面、域的偏偏移量.初初应力不但是仅是重力引起的.9、Tension-p shear-p tension-n shear-n表示的含意问:正在估计的循环内里,每个循环中,每个zone 皆依据failure criterion 处于分歧的状态,shear战tension 分别表示果受剪战受推而处于塑性状态.n 表示now,p表示previous,即分别表示正在当前战往日的循环核心处于塑性状态.10、history保存数据,print隐现正在过失的下令流中“”文献中惟有:“Step 1 Y-Displace”而不保存下数值,怎么回事?尔怎么才搞赢得table中的数据?输进'print tabel n'能得到table中的数据(粗确)hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0step 3000hist write 1 file 本面位移.dathist write 2 file 中部单元应力.dathist write 3 file 边沿单元应力.dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.dat(过失)hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0hist write 1 file 本面位移.dathist write 2 file 中部单元应力.dathist write 3 file 边沿单元应力.dathist write 1 table 1hist write 2 table 2hist write 3 table 3hist write 2 3 vs 1 file 单压.datstep 300011、plot隐现x、y坐标图hist n 1hist gp ydisp 0,0,0hist ax_strhist gp xdisp 1,1,0step 1500plot hist 2 vs -1 ;axial stress vs axial disp.plot hist 3 vs -1 ; circumferential disp. vs axial disp问:那2是不是代表hist gp ydisp 0,0,0,而-1 指得是hist gp ydisp 0,0,0为什么不是1而是背值?1) 正在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正背);2) "-"表示对于其值做"mirror",如果大家使用AUTOCAD的话便明黑"mirror"的含意了,比圆上头的 hist 1 gp ydisp的值是0~100,那么vs -1便形成-100~0.以此类推.3) 要那样搞主假如为了切合岩土工程的习惯需要,比圆桩载荷考查直线,重落直线等皆不是画正在第一象限内.12、关于water.fis的下令注释(树坐天下火位)问:(1)water.fis的功能是树坐天下火位里,正在二维估计中,通过water table下令即可树坐火位线,正在三维估计中由于天下火位是里,果此需要做一些插值运算,以便产死火位里.那个FISH步调包罗三个函数parm,aux1战aux2.parm用去设定火位的边值,即树坐火位的天区,aux1插值估计等值线的面,aux2调用aux1的估计截止,使用FLAC3D 下令产死火位里.(2)最先树坐参数值,而后定义一个table.def parmnptab1=7nprof=10endparmtable 1 (12 0) (18.25 12.5) (25 17) (35 21) (50 24) (75 25) (150 25)爆收table1,包罗7组数据面.3) 正在阐明的二个FISH函数之前,对于table,xtable,ytable做一些总的证明.water.fis步调反复使用了那三个变量,明黑了那三个变量的含意,water.fis那段步调便办理了大半问题了.Itasca 的所有硬件皆使用了相共的数据结构储藏数据,称之为链表结构,类似于关系型数据库中的数据链接要领.如果罕见据结构的知识,大概掌握关系型数据库的SQL 谈话,那三个变量的含意便不易明黑了,只不过那个支配与真真的数据库支配相比,惟有二个数据变量x,y.尔用数据库的术语做阐明:table 便是爆收一个“表”,xtable,ytable的效率便是爆收一个新表(Create)以及对于表中的数据举止与出(Select)、拔出(Insert)、替换(Update)支配.4)FISH函数aux1的效率是正在table 1的前提上爆收出其余个table,table的ID 从11到17,每个table有10个数据集.alfa的估计要领是为切合(fit)直里而采与的坐标面算法.def aux1loop n(1,nptab1)rr=xtable(1,n)zz=ytable(1,n)loop k(1,nprof)case_of kalfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)xtable(n+10,k)=rr*cos(alfa)ytable(n+10,k)=rr*sin(alfa)case 1xtable(n+10,k)=rrytable(n+10,k)=-40.0case 2xtable(n+10,k)=rrytable(n+10,k)=0.0end_caseend_loopend_loopendaux1当真止那个FISH函数后,死成(nptab1+10,nprof)组数据,那些数据为Aux2调用.注意:本段步调爆收的数据是为那个特定的直里边坡创造的,不是一个general purpose的LIB5)FISH函数aux2的效率是使用WATER TABLE face下令爆收火里,把table 1的ytable值动做z坐标,aux1爆收的数据动做x,y坐标.def aux2commandset grav 0 0 -10water dens 1000end_commandloop n(1,nptab1-1)loop k (1,nprof-1)xx1=xtable(n+10,k)yy1=ytable(n+10,k)zz1=ytable(1,n)xx2=xtable(n+10,k+1)yy2=ytable(n+10,k+1)zz2=zz1xx3=xtable(n+11,k+1)yy3=ytable(n+11,k+1)zz3=ytable(1,n+1)xx4=xtable(n+11,k)yy4=ytable(n+11,k)zz4=zz3commandWATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4end_commandend_loopend_loopendaux213、 FLAC截止数据的导出pl set back white ;树坐背景set plot bitmap size (800,510) ;树坐图片大小plot set caption size 37 ; 树坐图例数字大小 (default=35 范畴10-50) mainwin size 1 1 position 0 0.8 ; command windowplot set window size 0.8 0.8 position 0.1 0.0 ; graphics windowplot set center 100 0 40根据情况变动14、一段fish下令的阐明def aux1 定义函数loop n(1,nptab1) 启初一个循环,从1循环到变量nptab1rr=xtable(1,n) 从table 1 的x列中读第n个数付给变量rrzz=ytable(1,n) 从table 1 的y列中读第n个数付给变量zzloop k(1,nprof) 启初第二个循环,从1循环到变量nprofcase_of k 底下是一个采用结构,利用变量k推断alfa=0.5*(pi+.1)*float(k-2)/float(nprof-2)xtable(n+10,k)=rr*cos(alfa)ytable(n+10,k)=rr*sin(alfa)case 1 当k=1的时间真止底下语句xtable(n+10,k)=rr table n+10的x列第k个值等于变量rrytable(n+10,k)=-40.0 table n+10的y列第k个值等于变量-40case 2 当k=2的时间真止底下语句xtable(n+10,k)=rr table n+10的x列第k个值等于变量rrytable(n+10,k)=0.0 table n+10的y列第k个值等于变量0end_case 中断采用结构end_loop 中断第二个循环end_loop 中断第一个循环end 函数定义中断aux1 真止什么定义的函数15、有关下令的阐明(1)apply nvel 0.1 plane dip 60 dd 270 range left_bounapply 施加鸿沟条件,plane法背速度为0.1,plane 里的位子:倾角60(正在xoy仄里从-Z轴目标顺时针转动),dd 走背270度,从+y目标绕z轴顺时针转动;(2) ini szz -5e4 grad 0 0 -1e4ini 初初鸿沟条件,szz=-5e4+x*0+y*0+-z*-1e4,注:-5e4是range范畴起面坐目标SZZ值,x,y,z所施加应力网格面的坐标.16、渗流示例;------------------------------参数部分bulk modulus, K 390 MPa 体积模量(土量比较硬)shear modulus, G 280 MPa 剪切模量soil dry density, ρd 1200 kg/m3 土的稀度water density, ρw 1000 kg/m3 火的稀度wall density, ρwal 1500 kg/m3 档土墙的稀度permeability, k 10−12 m2/Pa-s 渗透系数porosity, n 0.3 孔隙率fluid bulk modulus, Kf 2.0 GPa 流体的体积模量;-----------------------以下是下令流;------------------------------------------------------; Excavation in a saturated soil;------------------------------------------------------config fluid ;树坐流体算法; --- geometrical model --- 建坐模型gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1group soil ;树坐土组group excavate range x 0 4 y 0 4 z 0 5 ;树坐启掘范畴组group wal1 range x 4 5 y 0 5 z 0 7 ; 树坐档土墙 1 组group wal2 range x 0 4 y 4 5 z 0 7 ;树坐档土墙 2 组group wall range group wal1 any group wal2 any ;档土墙1 战2 合并为一个档土墙组; --- fluid flow model --- 流体模型model fl_iso ;各项共性流体prop perm 1e-12 poro 0.3 ;树坐渗透系数战??比 or ??率ini fdensity 1e3 ; 初初化火稀度ini fmod 2.0e9 ftens -1e-3 ; 初初化火的体积模量及??model fl_null range group excavate ; 启掘部分不火流model fl_null range group wall ;档土墙不渗火ini pp 0 grad 0 0 1e4 ;初初化孔隙火压力fix pp range z -0.1 0.1 ;拘束该当是所有模型的底部吧?fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;启掘部分也给拘束住了,而其余的场合不管; --- mechanical model --- 力教模型model elas ; 弹塑性prop bul 3.9e6 shea 2.8e6 ;树坐体积模量战剪切模量,那二个必须,有一次尔只树坐了体积模量,便错的一塌糊涂model null range group excavate ;力教模型的启掘ini density 1.2e3 ;初初化稀度,该当是土的稀度ini density 1.5e3 range group wall ;初初化墙(档土墙)的稀度,fix x range x -.1 .1 ;施加拘束,传闻好分法拘束0 的时间,必须与一个-0.1<0<0.1,范畴fix x range x 11.9 12.1 ;施加拘束,大概者喊搞鸿沟条件fix y range y -.1 .1 ;施加拘束fix y range y 11.9 12.1fix z range z 11.9 12.1; initial total stresses 出世总应力ini szz 0 grad 0 0 -1.5e4 ;z目标的应力大一些ini sxx 0 grad 0 0 -1.2e4 ;火仄目标的,小一些ini syy 0 grad 0 0 -1.2e4 ;apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1; --- setting --- 树坐set gravity 0 0 10 ;树坐重力加速度; --- initial state --- ;初初化状态solve force 1 ; check initial equilibrium 供解仄稳; --- histories --- 记录set hist_rep 40 ;记录步为40hist fltime ;记录时间hist gp pp 0 0 6 ;记录某面的孔隙火压力hist gp xdis 4 0 0 ;火仄位移hist gp xdis 4 0 2 ;记录火仄位移hist gp xdis 4 2 0hist gp zdis 0 0 5 ;记录横背位移hist gp zdis 2 0 5hist gp zdis 4 0 5hist gp zdis 2 2 5hist gp zdis 4 2 5hist gp zdis 4 4 5hist gp zdis 10 0 1hist gp zdis 10 0 2;; --- excavation --- 启掘set fluid off ;关关火; apply pore pressure at walls 正在墙上施加孔隙火压力apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0 apply nstress -5.e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1solve ;force 1 ;供解save exc1.sav ; 保存;; --- drainage --- 排火apply remove nstress ;撤掉刚刚才的压力def relaxsetup ;定义一个释搁函数的参数step0 = stependrelaxsetup ;调用def relax ;定义一个释搁函数rstep = step - step0if rstep < ncyc thenrelax=1.0-(float(rstep)/float(ncyc))elserelax = 0.0endifend ;中断set ncyc = 1000 ;先赋值,随后调用apply nstress 0 grad 0 0 -1.e4 hist relax &range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.e4 hist relax &range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -5.e4 hist relax &range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1cyc ncycsolve ;供解save exc2.sav ;保存;; --- percolation --- 渗透set fluid on ;树坐火算法fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙火压力cyc 9000 ;循环save exc3.sav ;保存plot create excav ;建坐一个隐现的视图一下真量皆是隐现真量plot set rot 200 0 195plot set cent 6 6 6plot set dist 39.18plot set magn 0.8plot set plane ori 0 0 0 normal 0 -1 0plot add cont pp out onplot add block group lgra range group wallplot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7 plot add flow planeplot showRet17、用什么下令不妨位移图中隐现最大位移面的坐标?问:用fish编了一个,算法佳像有面烦,不过脚段不妨真止您不妨参照一下.def get_gp_maxdispgp0_disp = gp_xdisp(gp_head)*gp_xdisp(gp_head)gp0_disp = gp0_disp + gp_ydisp(gp_head)*gp_ydisp(gp_head) gp0_disp = gp0_disp + gp_zdisp(gp_head)*gp_zdisp(gp_head) gp0_disp = sqrt(gp0_disp)p_gp=gp_head;找最大值loop while p_gp # nullgp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp = sqrt(gp_disp)if gp_disp>gp0_dispgp0_disp=gp_dispendifp_gp = gp_next(p_gp)endloop;找最大值的坐标p_gp=gp_headss=0loop while p_gp # nullgp_disp = gp_xdisp(p_gp)*gp_xdisp(p_gp)gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp)gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp)gp_disp = sqrt(gp_disp)if gp_disp=gp0_dispss=ss+1xtable(1,ss)=gp_xpos(p_gp)ytable(1,ss)=gp_ypos(p_gp)endifp_gp = gp_next(p_gp)endloopendget_gp_maxdisp18、怎么样隐现估计后模型的变形状态?问:pl sk magfac **(一个搁大的数字)19、怎么样查看剖里上任一面的位移值问:print gp disp range x() y( ) z( )大概 print gp disp range idplot gp disp range x() y( ) z( )大概 plot gp disp range id20、group下令区分斜里问:group groupname colour range plane above dd ** dip ** origin *** plane below dd ** dip ** origin ***dd 战dip便是二个斜里的倾背战倾角,origin是那个斜里通过的某个面,那个下令便不妨group二个斜里间的天区21、ansys导进flac的时间,怎么样区分一下group问:正在ansys中对于分歧group用分歧的资料号(mat)定义,再用ansys2flac之类的步调变换成flac的.dat 文献便可,group自动死成. 22、 fix x y z range x -0.1 0.1 ..........................afix x range x -0.1 0.1................................b请问a式与b式有何辨别?问:a 表示正在x=0 那个仄里是牢固的b 表示正在x=0 那个仄里只x目标不克不迭爆收疏通.23、怎么样了解每个单元的ID?问:用鼠标单打单元的表面,不妨了解单元的ID战坐标.24、怎么样举止切片?问:plo set plane ori (面坐标) norm (法背矢量)plo con sz plane (隐现z目标应力的切片)那二止下令,用去隐现指定里上的位移变更云图(把第二止改改,还不妨隐现应力云图,等其余);plot con syy y背应力云图plot con szz outline on z背应力云图加隐现网格plot con szz effective on z背灵验应力plot con szz inter 1e4 云图隐现删量(后里的值越小,图上颜色条戴便越细,)plot con szz min -10e4 隐现最小值,图上小于那个数的值皆用共一色表示plot con szz max 10e4 隐现最大值,图上大于那个数的值皆用共一色表示plot fip outside on 隐现力的箭头当前瞅鸿沟,如果是不动里,用fix x range…如果鸿沟戴速率,大概位移用ini…ini dens 1000树坐重力便一定要树坐稀度,不思量重力便不要安排稀度.二者共存的,安排步调的时间一定要记的.set mech force 50solve用SET下令设定了最大不仄衡力的上限值为50N,即:估计的最大不仄衡力一往小于此值,咱们便人为模型达到了仄稳状态.供解历程末止.树坐那个脚段是正在一定粗度央供下,加快估计(值越大,算的量越少,花的时间也越少).25、怎么样正在步调中举止久停,并可回复估计?问:正在下令中加进pause下令,用continue举止继承.正在咱们分步供解中念得到某一个历程中的截止,不必等到齐供完,还不妨正在分集供解过失的时间便举止改正,而不是等到截止出去. 26、怎么样跳过某个估计步?问:正在估计中按空格键跳过本次估计,自动加进下一步27、Fish是什么物品?Fish是可一定要教?问:是FLAC3D的内置谈话,不妨用去举止参数化模型、完毕下令自己不克不迭举止的功能.Fish不妨不必教,需要的时间查Mannual赢得需要的变量便不妨了.28、怎么样调用一些可选模块?问:config dyn (fluid, creep, cppudm)29、range的下令:①与圆球里的范畴,便是共一个核心,分歧半径之间夹的那一部分,用正在给球里施加力,拘束时用该下令,下令:range annulus center xyz(核心坐标)radius r1 r2(二半径),念要表黑的是球里时,r1r2与值出进小面即可,如radius 2.9 3.1 那时范畴便是半径为3的里.②圆柱的环里施加力大概拘束时也用到该下令.range cylinder end1 xyz end2 xyz radius r(end是表示柱体断里的核心面坐标)range plane above 表示的是z大于0的上部结构range plane blow 表示的是z小于0的下部结构range dd a dip b origin xyz 表示的是走背角度a,倾背角度b的里dd表示,劈里与xoy里相接时的垂线与y轴正背的夹角,顺为正,道心语,便是里的指背目标里是往那个目标的意义.dip里与xoy里的夹角,与z轴背背为正,origin表示里通过xyz面;range normal xyz 法线目标为xyz的里range not 出试用过(范畴中的物品,念念用到圆上与范畴非常佳). 30、怎么样调用FISH文献?问:新脚调用FISH文献堕落的本果:是:文献名为fishcall.fis.txt(果为正在记事本里那个txt的后缀是默认的).建改要领是:将fishcall.fis.txt 内里的真量局部复造到新的记事本里,再按保存,保存时间将保存典型改为所有文献(*,*),那时您即可正在文献名栏内输进fishcall.fis,典型也便粗确了.31、问a式与b式有何辨别?问:fix x y z range x -0.1 0.1 ..........................afix x range x -0.1 0.1................................ba 表示正在x=0 那个仄里是牢固的;b 表示正在x=0 那个仄里只x目标不克不迭爆收疏通.32、怎么样查看各个时段不仄衡力的简直数值?问:采与his去记录估计,包罗位移应力等下令his unbal;his gp(zone) zdis range (0 0 0) 大概者id=?;导出数据下令his write n vs m begin 时步 end 时步+ file filename.hisn表示记录的id m表示时步;要导出不仄衡力的简直数值his unbalstep 100000 or solvehis write 1 vs step begin 1 end 1000 file 123.his使用上述下令便不妨查看各个时步下的不仄衡力的简直数值.33、怎么样隐现变形表面线的下令?plo ske magf 10 其中10为搁大系数34、何如查看模型的矢量图?pl grid vel (disp) redplo blo gro不妨瞅到分歧的group的颜色分集(瞅模型的格式);plo gpfix red sk 瞅模型的鸿沟情况plo fap red sk 瞅模型的体力分集plo con dis (xdis, ydis, zdis) 瞅模型的云图:位移plo con sz (sy, sx, sxy, syz, sxz) 瞅模型的云图:应力plo dis (xdis, ydis, zdis) 瞅模型的矢量图plo info 瞅模型有几单元、节面plot grid 不妨查看网格,plot grid num 不妨查看节面号.plot fix 查看鸿沟条件plot grid mark 查看预先留设的启掘鸿沟pl syy (sxx xdis ydis sig1..) fill;查看模型的云图.注意,戴fill关键字便是查看云图,不戴fill 关键字便是查看等直线图.plo add ske 剖里的重落云图上减少网格plo add dis plane 剖里上的矢量plo add axe 剖里上坐标值35、何如查看INTERFACE?问:pl grid iface id redprint iface id ; id 是interface的编号36、FLAC怎么安排视图?问:ctrl+z 大概shift+z ,而后面打要搁大部位的左上角战左下角便不妨搁大某个部位,如果要缩小,单打屏幕. 37、推断模型中网格面是一个仍旧二个的要领.问:隐现模型网格面的ID,如果有二个ID号,证明是二个.节理里上的网格面该当有二个ID号.38、flac中效率运止时间的果素问:(1)单元资料,结构里之间的刚刚度对于比很大;(2)区分单元尺寸对于比很大,尽管预防少细比大于5:1的细少单元.39、简朴下令阐明:问:1.plot add cont disp out on 常常隐现云图2.plot block group gpnum on 隐现节面编号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 n44、念记录交战里上的相关力教参数该怎么办?问:hist inter i nstress /ndis/sstress/sdisp45、交战里树坐参数包罗哪几个?问:coh dial fric kn ks ten.46、怎么样输出单元格的数据?问;expgrid .flacsd47、建模的历程中,大概共一个文献中有重复应用的下令,便不妨用宏下令了.问:宏下令: 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 name49、plot中add与set的辨别add本去有的加新的用addset本去不的,树坐用setplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20plot set center 2.5 4.2 4.0plot add bound behindplot add bcont szz planeplot add axesplot show .咱们念隐现别的图像,又不念把那个弄佳的图像拾掉,则用plot create name建坐一个图像,正在建坐新的图像时,plot create name1(新的名)隐现往日的图像,那么输进如下下令:plot current name.另一种要领,输进:plot show name50、分歧的天区定义分歧的本构模型model mohr range z=0,5prop ClayeyGravel range z=0,5正在定义参数后里加个range范畴,为了预防正在输进下令时总要通过range+局部坐标域的要领去规定该下令真止的天区,咱们不妨先通过range下令去给分歧的天区命名,而后便不妨用range+天区名的要领去规定范畴了.那样便预防了屡屡输进下令要规定正在该天区真止时,皆要输进一遍坐标范畴的贫苦.简直怎么样给天区命名,以及日后怎么样用天区名去规定该天区,不妨瞅底下那个例子:range name trench x=-1,1 y=0,4 z=-2,2model null range trench51、左脚规则:建坐模型时,各面的排搁时有顺序的.要领如下:拿出左脚伸直了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 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 3Table 1 0 45 .05 42 .1 40 1 40。
flac3d相关命令
Flac3D基础知识FLAC-3D1. Generate 命令产生网格,简写为GE。
2. zone 参数代表了三维网格体(单元体),简写为zo。
3. (1)size 参数指示单元体各边小单元体的个数;(2)dimension 参数指示单元体内部巷道各边边长,如果没给出就默认为两点间距的20%;(3)edge 参数指示网格体各边的尺寸;(边长也可以用确定点坐标来利用点与点的间距来确定p0 p1 p2……)。
4.各种网格体名称及关键字:矩形网格(Brick),退化矩形网格(dbrick),楔形网格(wedge),金字塔网格(pyramid),四面体形网格(tetrahedron),圆柱体网格(cylinder),矩形体外环绕放射状网格(radbrick),平行六面体外环绕放射状网格(radtunnel),圆柱体外环绕放射状网格(radcylinder),柱形壳体网格(Cshell),交叉圆柱体网格(cylint),交叉平行六面体网格(tunint)。
5.绘制网格体语句:①Generate zone brick size a,b,c (其中的brick可以被其他的网格体替换)。
②Generate zone brick p0=(,,)&P1=(,,)&P2=(,,)& (边长由点点间距确定)P3=(,,)& →(&符号为继续行的符号)Size a b c →(各边小单元体个数)③如果在p1 p2 p3 前面加add 则各点坐标是相对于p0点的坐标。
④如果只给出p0 p1 则由edge a & 给出边长。
(见第 3 条);利用ratio a,b,c 给出小单元体长度比例依次变化。
⑤内部有空单元体或者巷道的网格体,如果内部巷道需要填充为单元体,使用fill group groupname 语句。
(其中groupname为巷道内填充单元体的组名)。
⑥显示图形:plot set rotation a b c 使图形可以围绕一点旋转Plot magnification a 设置图形比例大小Plot add axes red 设置坐标系及其颜色Plot show 显示图形Plot add surface red 设置图形表面颜色,设置坐标轴某一范围的表面颜色时利用range参数(课本39页4-1),add surface red range z= a b 或者课本42页的plot block group yellow red⑦Generate zone brick p0=(,,)size a b 1 一层陪衬的底网格体。
Flac3D5.0操作方法笔记
FLAC3D5.0导入dat、txt文件的方法
FLAC3D5.0导入.sav文件的方法
输入命令流:Restore 文件名.sav
FLAC3D5.0导入.f3prj文件的方法
视图窗口关闭后,如何再次呈现?输入命令plot zone
当命令流输错可从这里撤销
点亮相应命令流后才可撤销
透明度的调整
保存项目文件
保存文件的重命名
做切片
倾向和倾角
法向和过一个点
删除切片
Ctrl+R为查看模型后还原
显示接触面,双击I nterface显示接触面
接触面属性项被激活后,可以更改属性
隐藏网格显示
输出网格
如果点击了透明,将不显示Information
只要不点击透明,则显示Information
等同于
点击显示节点坐标
fish简介
@之前要有(空格)之前加@是Flac3D3.0和5.0区别
而在之间没有(空格)
按下列步骤取消fish语言前要加@
点击fish
软件的刷新
如何显示软件最终的安全系数
显示软件节点.
查看应力及最大剪应力的步骤
查看位移云图的步骤
输出云图
设置出图边宽
设置分辨率可设为1080 Reflect镜像restore调用之前的计算结果
Origin原点normal法向量(该点与原点的坐标差)
对于结构单元如桩,若点击Zone将不出现任何模型
按如下步骤才可显示
使用该命令后保存的文件与说建模的文件保存在同一文件夹
查看约束反力的步骤
查看分组的步骤
位移云图的显示颜色相反如何调整
连接两个不同尺寸的单元。
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 end2radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
FLAC3D命令流(挺实用)(DOC)
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 fillplot surfgen zone reflect dip 0 dd 90 (对xy面做镜像)gen zone reflect dip 90 dd 90 (对yz面做镜像)(1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
建模1、调用文件:①文件与工程在同一个文件夹,只写文件名即可:Ifthecalledfileislocatedinthesamefolderasthe FLAC3D projectfile,thenonlyt hefilenameneedbeenteredwiththe CALL command.②不在同一个文件夹,全路径:Otherwise,thefilemaybecalledbyspecifyingitscompletepath(e.g.,c:\myfol der\file.dat).Undo;撤销上一条命令2、创建旋转缩放视图3、建模命令modelmechmohr;莫尔库伦模型modelmechelastic;弹性模型setgrav0,0,-9.81;重力加速度negative z-direction.(垂直向下!常用的)下下面面这这代代码码,,是是沿沿着着--y y方方向向的的重重力力加加速速度度,,注注意意区区别别!!!!!!!!genzonebricksize6,8,8p0-10,-10,-20...;省略号表示写不下后面继续p110,-10,-20...p2-10,10,-20...p3-10,-10,0plotzonegenzonebricksize6,8,8p0-10,-10,-20...;不规则六面体p110,-10,-20p2-10,10,-20...p3-10,-10,0p410,10,-20...p5-10,10,10p610,-10,0...p710,10,10plotcurrentplotPlot01plotclearplotzoneUndo;撤销命令setlogfile127G1001.tGtsetlogontruncatesetlogofflistzoneprinrangeG01y01z01;显示指定范围内各单元的主应力,结果如下Hist命令:①命令编号按顺序从1开始:eachhistoryisnumberedsequentiallyfrom1asitisenteredviathe HISTORY co mmand.②查找显示所有的his命令:ReturntotheFlac3D>promptandtypelisthist foralistingofthehistoriesandtheircorrespondingnumbers.histnstep5;每5步记录1次。
默认是10步记录1次histaddratio;记录maGimumunbalancedforceratiohistaddgpzdisp448;记录节点的位移modelnullrangeG=2,4y=2,6z=5,10;ToeGcavatethetrench,enter ThedimensionsoftheeGcavationregion(i.e.,thenullmaterial)aredefinedbyt he G-,y-and z-coordinaterange.setlarge;设置大变形large-strainlogicisspecifiediniGdis0ydis0zdis0;位移归零solvesavetrench;保持计算结果的文件“TRENCH.F3SAV”不需要后缀4、group命令(注意新版本的区别,groupzone)genzonebricksize666groupzoneTunnelrangeG16y06z15;设置groupmodelmechmohr;莫尔库伦模型modelmechnullrangegroupTunnel; groupzoneTunnelrangecylinderend1100105100end210095100radius3 groupzoneMatriGrangegroupTunnelnot;namingRANGESthatarefiGedinspaceandunrelated;tothemodelgeometry-representsstratigraphy...rangenameLayer1planedip0dd0ori000aboverangenameLayer2planedip0dd0ori000belowpropbulk1e8shear0.5e8coh0tens0fric35rangenrangeLayer1propbulk1e7shear0.3e7coh1e7tens0fric0rangenrangeLayer2 plzonecolorbygroup4、Interface命令以下是切割模型法生成interface:寺河模拟用的plzonegenzonebrickp0000p120KK0p20400p3004size100202group1; genzonebrickp0004p120KK4p20404p3008size100202group2; genseparategroup1interface1wrapfirstgroup1secondgroup2interface1maGedge0.5显示interface:plotcreateplotPlot01plotaddsurface;显示表面plotaddinterfacered;显示交互面为红色plotshow用导来导去法,已经在5.0上运行过,很好,值得借鉴plzonecolorbygroupgenzonebricksize333&p0(0,0,0)p1(3,0,0)p2(0,3,0)p3(0,0,1.5)&p4(3,3,0)p5(0,3,1.5)p6(3,0,4.5)p7(3,3,4.5)groupzoneBasegenzonebricksize333&p0(0,0,2.5)p1(3,0,5.5)p2(0,3,2.5)p3(0,0,7)&p4(3,3,5.5)p5(0,3,7)p6(3,0,7)p7(3,3,7) groupzoneToprangegroupBasenotinterface1facerangeplanenorm(-1,0,1)origin(1.5,1.5,3)dist0.1 pladdinterfacepauseini zpos add-1.0rangegroupTop4、hist命令(5.0版本需要加add!!!)histnstep=5//设置采用间隔为5等效于sethist_rep=5 histaddid=2zoneszz0,0,-1histunbal//最大非平衡力histaddgpzdisp448printhisthistaddgpzdis0,0,0histaddgpzvel0,0,0histaddgpydisp0,0,0;监测网格点(0,0,0)点的y方向上的位移histaddzonesyy0,1,0;监测单元(0,1,0)点的y方向上的应力histaddzonesyy1,1,0;监测单元(1,1,0)点的y方向上的应力histaddzoneszz0,0,-1Listhistorysetlargehistorynstep2historyaddfishtthistoryaddfishtheta historyaddzonesGG(1,1,1) historyaddzonesyy(1,1,1) historyaddzoneszz(1,1,1)plotaddhis345vs1HI STORY addfish var;fish var变量的监测5、solve命令?????什么意思5、plot命令通过图形操作界面的方式:plzonecolorbygr;查看各个分组通过命令流的方式:plotcreateplotPlot01plotclearplotaddhistory1linestylecolorblack plotpostplotPlot01file127F1002.psmanual plotcreateplotPlot02plotclearplotaddhistory2linestylecolorblackplotpostplotPlot02file127F1003.psmanual plotcreateplotGravVplotcutaddplanedip90dd0origin3,4.1,0;截面plotsetcenter333plotseteye121812 plotaddzoneplanebehindplaneononplaneoff;addlabel"Zone"white plotaddbcontszzrampgreyplotaddaGesplotpostplotGravVfile127F1004.psmanualreskw50plotcreateplotszz_contourplotcutaddplanedip90dd180origin0,18.9,0plotsetorient1202200plotsetcenter45450 plotaddzoneplanebehindplaneononplaneoff; plotaddbcontszzramprainbowplsetmag1.0plotaddaGesplotshowplotcreateGravVplotsetplanedip=90dd=0origin=3,4,0plotsetrot15020;thiswouldbeachievedinteractivelyplotsetcenter2.54.24.0;andisshowheresimplytoillustrate plotaddboundbehindplotaddbcontszzplaneplotaddaGesplotshow这样我们就建立了一个叫“GravV”的图形,将它设为当前图层,对此图形设置一个倾角为90°的面,(在G-y面上,约定z的负方向为下),倾向为0°。
(在G-y面上,从y的正方向开始沿顺时针方向计量),过G-y面上一点(G=3,y=4,z=0),在平面后加一个网格,垂直方向的应力Qzz,最后给模型加轴线,整体图形轮廓与添加图形相反,显示每个区域中心的应力计算值,每个区域的颜色直接对应基本区域的应力,当键入show,我们建立的图形被显示PlotcreateDispcontPladdcontdispplanebehindPladdaGesPlshowplotsketch;绘出轮廓线plotaddvel;添加矢量plotaddhis12vs3;绘出历史记录n,actLoad,p_solcyc15000plotCreateTrenchAddsurfaceyellowAddaGesblackShowplotcreateszz_contour;创建文件名为szz_contour(z方向的应力云图)plotsetcent445;视图的中心坐标为(4,4,5)plotsetrot20KK0;视图的旋转角度为(20,0,30)plotsetmag1.0;视图的放大倍数为1.0 plotaddcontszzoutonshadeon;打开szz的等值线图plotaddaGes;打开坐标轴线plotshow;显示窗口;不同的区域,材料赋值不同plzonecolorbymodel plzonecolorbygroup求解前对模型进行检查网格检查(如前所述)边界条件检查速度约束条件plotgpfiG荷载条件plotfap模型检查模型赋值ploblockmodel参数赋值ploblockpropGGG初始应力检查(如后所述)经常检查模型的响应plozconszz(syy,sGG) 应力场ploconzdis(ydis,Gdis) 位移场plo(add)zonecolorbystate 屈服状态plogpfiG 速度约束条件plofap 体力plohist(unbal)1 不平衡力(必需的)plbitmapview1filemyplot.png plbitmapview1size1024768 plcreate<viewid>。