(完整word版)FLAC3D常见命令与使用技巧
[实用参考]Flac3d-5.0常用命令集锦.doc
![[实用参考]Flac3d-5.0常用命令集锦.doc](https://img.taocdn.com/s3/m/dee971461711cc7930b7163b.png)
建模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次。
(完整word版)FLAC3D命令流

(完整word版)FLAC3D命令流命令流1gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 &p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 &size 10 10 10 rat 1。
0 0。
9 1.1 group brick_1gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 &size 10 10 10 rat 1。
0 0.9 1.1 group brick_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0。
9 1.1 group brick_3plot sur2-1newgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z —。
1 。
1fix x ran x -。
1 。
1fix x ran x 2。
9 3。
1fix y ran y —.1 。
1fix y ran y 2.9 3.1set grav 0 0 —10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve3-1;——-—-—————--——-—--—-———--————-——-———-工程信息;Project Record Tree export;Title:Simple test ;—-———-—————---—-—--———————--—---—计算第一步;。
.。
STATE: STATE1 。
.configgrid 10,10model elastic(完整word版)FLAC3D命令流group 'User:Soil’ notnullmodel elastic notnull group 'User:Soil'prop density=1500.0 bulk=3E6 shear=1E6 notnull group ’User:Soil’fix x y j 1fix x i 1fix x i 11set gravity=9.81history 999 unbalancedsolvesave state1.sav;—-------—---——----—----------—---—计算第二步;... STATE:STATE2 ..。
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、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命令流(整理版)

1、怎样查看模型答:plotgrid可以查看网格,plotgridnum可以查看节点号。
2、请问在圆柱体四周如何施加约束条件答:可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
命令流如下:fixxrangeend1100end2140rad1end1100end2140rad13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图答:用setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后,如何显示点的坐标答:使用ploblogrogpnumon命令5、关于gauss_dev对性质进行高斯正态分布的问题答:根据手册上的说明:下面的命令设定一个平均摩擦角为40度,标准方差是±5%。
则命令如下:propfriction40gauss_dev2问题:请问gauss_dev2中的2是如何计算的如果把±5%改为±10%,则命令应如何写40×5%=26、reflect问题问:genzoneradbrick&p0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)&size3,5,5,7&rati o1,1,1,&dim142fillplotsurfgenzonereflectdip0dd90(对xy面做镜像)genzonereflectdip90dd90(对yz面做镜像)(1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线。
首先应该按照dd 的方向大体确定这个面的朝向,dd指的是从y轴正方向按顺时针(clockwise)方向转向所要确定面的法线方向在xy平面上的投影的夹角,然后再确定dip,dip指的是从xy平面转向所要确定的平面的角.(及z轴负方向转过角度)(2)命令:reflectnormalxyzoriginxyz(根据法线和过一点建的面)最常用,至于实际使用过程中有人问道:normal-100与normal100的区别,我试过两者的效果是一样的,没区别,虽然方向不同,但表达的是同一个面。
[实用参考]Flac3d-5.0常用命令集锦
![[实用参考]Flac3d-5.0常用命令集锦](https://img.taocdn.com/s3/m/0be3be44f46527d3240ce044.png)
建模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次。
(完整word版)FLAC3D5.0模型及输入参数说明

1.1模型参数代码可参考manual中各个章节的command命令及说明,注意单位.用prop 赋值。
1.1.12经典粘弹性模型的材料参数(Classical Viscoelastic (Maxwell Substance) –MODEL mechanical viscous)1bulk弹性体积模量,K2shear弹性剪切模量,G3viscosity动力粘度,η1.1.13粘弹性模型粘弹性模型的材料参数(Burgers Model –MODEL mechanical burgers)1bulk弹性体积模量,K2kshear Kelvin弹性剪切模量,G K3kviscosity Kelvin动力粘度,ηK4mkshear Maxwell切边模量,G M5mviscosity Maxwell动力粘度,ηM1.1.14二分幂律模型二分幂律模型的材料参数(Power Law –MODEL mechanical power)1a_1常数,A12a_2常数,A23bulk弹性体积模量,K4n_1指数,n15n_2指数,n26rs_1参考应力,σ1ref7rs_2参考应力,σ2ref8shear弹性剪切模量,G1.1.15蠕变模型蠕变模型材料参数(WIPP Model –MODEL mechanical wipp)1act_energy活化能,Q2a_wipp常数,A3b_wipp常数,B4bulk弹性体积模量,K5d_wipp常数,D6e_dot_star临界稳定状态蠕变率,7gas_c气体常数,R8n_wipp指数,n9shear弹性剪切模量,G10temp温度,T下列参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率Burger、蠕变组合材料模型的材料参数(Burgers—Creep Viscoplastic Model –MODEL mechanical cvisc)1bulk弹性体积模量,K2cohesion内聚力,c3density密度,ρ4dilation剪胀角,Ψ5friction内摩擦角,Φ6kshear Kelvin弹性剪切模量,G K7kviscosity Kelvin粘度,ηK8shear弹性剪切模量,G9tension抗拉强度,σt10mviscosity Maxwell动力粘度,ηM下列计算参数可以显示、绘图和通过fish访问1es_plastic累积塑性切应变2et_plastic累积塑性拉应变1.1.17幂律模型幂律模型的材料参数(Power-Law Viscoplastic Model –MODEL mechanical cpower)1a_1常数,A12a_2常数,A23bulk弹性体积模量,K4cohesion内聚力,c5dilation剪胀角,Ψ6friction内摩擦角,Φ7n_1指数,n18n_2指数,n29rs_1参考应力,σ1ref10rs_2参考应力,σ2ref11shear弹性剪切模量,G12tension抗拉强度,σt1.1.18粘塑形模型粘塑形模型的材料参数(WIPP—Creep Viscoplastic Model –MODEL mechanical pwipp)1act_energy活化能,Q2a_wipp常数,A3b_wipp常数,B4bulk弹性体积模量,K5d_wipp常数,D6e_dot_star临界稳定状态蠕变率,7gas_c气体常数,R8kshear材料参数,KΦ9n_wipp指数,n10kdil材料参数,q k11kvol材料参数,qΦ12shear弹性切变模量,G13temp温度,T14tension抗拉强度,σt以下计算参数可以显示、绘图和通过fish访问1e_prime累积主蠕变应变2e_rate累积主蠕变应变率3es_plastic累积塑性切应变4et_plastic累积塑性拉应变碎盐变形模型的材料参数(Crushed—Salt Model –MODEL mechanical cwipp)1act_energy活化能,Q2a_wipp常数,A3b_f最终体积模量,K f4b_wipp常数,B5b0蠕变压实系数,B06b1蠕变压实系数,B17b2蠕变压实系数,B28bulk弹性体积模量,K9d_f最终密度,ρf10d_wipp常数,D11e_dot_star临界稳定状态蠕变率,12gas_c气体常数,R13n_wipp指数,n14rho密度,ρ15s_f最终切变模量,G f16shear弹性切变模量,G17temp温度,T以下计算参数可以显示、绘图和通过fish访问1frac_d当前碎片密度,ρd2s_g1蠕变压实参数,G3s_k1蠕变压实参数,K均质流体模型的材料参数1permeability等方向渗透性,k2porosity孔隙率,n(默认时,n=0。
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平面的投影是一条射线。
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或者什么作图软件绘制就行了。
完整word版,FLAC3D 实例命令流1

第1部分命令流按照顺序进行2-1定义一个FISH函数newdef abcabc = 25 * 3 + 5Endprint abc2-2使用一个变量newdef abchh = 25abc = hh * 3 + 5EndPrint hhPrint abc2-3对变量和函数的理解newdef abchh = 25abc = hh * 3 + 5Endset abc=0 hh=0print hhprint abcprint hhnewdef abcabc = hh * 3 + 5endset hh=25print abcset abc=0 hh=0print hhprint abcprint hh2-4获取变量的历史记录newgen zone brick size 1 2 1model mohrprop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10fix x y z range y -0.1 0.1apply yvel -1e-5 range y 1.9 2.1plot set rotation 0 0 45plot block groupdef get_adad1 = gp_near(0,2,0)ad2 = gp_near(1,2,0)ad3 = gp_near(0,2,1)ad4 = gp_near(1,2,1)endget_addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) endhist loadhist gp ydis 0,2,0step 1000plot his 1 vs -22-5用FISH函数计算体积模量和剪砌模量newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25deriveprint b_modprint s_mod2-6 在FLAC输入中使用符号变量Newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))endset y_mod = 5e8 p_ratio = 0.25derivegen zone brick size 2,2,2model elasticprop bulk=b_mod shear=s_modprint zone prop bulkprint zone prop shear2-7 控制循环Newdef xxxsum = 0prod = 1loop n (1,10)sum = sum + nprod = prod * nend_loopendxxxprint sum, prodnewgen zone brick p0 (0,0,0) p1 (-10,0,0) p2 (0,10,0) p3 (0,0,-10) model elasplot set rotation 0 0 45plot block groupdef installpnt = zone_headloop while pnt #nullz_depth = -z_zcen(pnt)y_mod = y_zero + cc * sqrt(z_depth)z_prop(pnt, ’shear’) = y_mod / (2.0*(1.0+p_ratio))z_prop(pnt, ’bulk’) = y_mod / (3.0*(1.0-2.0*p_ratio))pnt = zone_next(pnt)end_loopendset p_ratio=0.25 y_zero=1e7 cc=1e8install2-8 拆分命令行new ;example of a sum of many thingsdef long_sumtemp = v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 + v9 + v10 long_sum = temp + v11 + v12 + v13 + v14 + v15end2-9 变量类型newdef haveoneaa = 2bb = 3.4cc = ’Have a nice day’dd = aa * bbee = cc + ’, old chap’endhaveoneprint fish2-10 IF条件语句newdef abcif xx > 0 thenabc = 1000elseabc = -1000end_ifendset xx = 10print abcset xx = 0print abc2-11 索单元自动生成newgen zone brick size 10 3 5plot set rotation 0 0 45plot block groupdef place_cablesloop n (1,5)z_d = float(n) - 0.5commandsel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7 end_commandend_loopendplace_cablesplot grid sel geom rednewgen zone brick size 10 3 5plot set rotation 15 0 60plot block groupmod mohrprop bulk 1e8 shear .3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav 0,0,-10fix x y z range z -.1 .1fix y range y -.1 .1fix y range y 2.9 3.1fix x range x -.1 .1fix x range x 9.9 10.1set largehist unbalsolvesave cab_str.savini xdis 0 ydis 0 zdis 0hist gp xdisp 0,1,5def place_cablesloop n (1,5)z_d = 5.5 - float(n)z_t = z_d + 0.5z_b = z_d - 0.5commandfree x range x -.1,.1 z z_b z_tsolvesel cable beg 0.0,0.5,z_d end 7.0,0.5,z_d nseg 7sel cable beg 0.0,1.5,z_d end 7.0,1.5,z_d nseg 7sel cable beg 0.0,2.5,z_d end 7.0,2.5,z_d nseg 7sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 &gr_k 2e10 gr_coh 1e10 gr_per 1.0end_commandend_loopendplace_cablessave cab_end.savplot sketch sel cable force red2-12圆形隧道开挖模拟计算;建立模型gen zon radcyl p0 0 0 0 p1 6 0 0 p2 0 1 0 p3 0 0 6 &size 4 2 8 4 dim 3 3 3 3 rat 1 1 1 1.2 group outsiderockgen zone cshell p0 0 0 0 p1 3 0 0 p2 0 1 0 p3 0 0 3 &size 1 2 8 4 dim 2.7 2.7 2.7 2.7 rat 1 1 1 1 group concretliner fill group insiderock gen zon reflect dip 90 dd 90 orig 0 0 0gen zon reflect dip 0 dd 0 ori 0 0 0gen zon brick p0 0 0 6 p1 6 0 6 p2 0 1 6 p3 0 0 13 size 4 2 6 group outsiderock1 gen zon brick p0 0 0 -12 p1 6 0 -12 p2 0 1 -12 p3 0 0 -6 size 4 2 5 group outsiderock2 gen zon brick p0 6 0 0 p1 21 0 0 p2 6 1 0 p3 6 0 6 size 10 2 4 group outsiderock3 gen zon reflect dip 0 dd 0 orig 0 0 0 range group outsiderock3gen zon brick p0 6 0 6 p1 21 0 6 p2 6 1 6 p3 6 0 13 size 10 2 6 group outsiderock4gen zon brick p0 6 0 -12 p1 21 0 -12 p2 6 1 -12 p3 6 0 -6 size 10 2 5 group outsiderock5 gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z 6.1 13.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x -0.1 6.1 z -6.1 -12.1gen zon reflect dip 90 dd 90 orig 0 0 0 range x 6.1 21.1 z -12.1 13.1;绘制模型图plot block groupplot add axes red;plot set rotation 0 0 45 用于显示三维模型;设置重力set gravity 0 0 -10;给定边界条件fix z range z -12.01,-11.99fix x range x -21.01,-20.99fix x range x 20.99,21.01fix y range y -0.01 0.01fix y range y 0.99,1.01;求解自重应力场model mohrini density 1800 ;围岩的密度prop bulk=1.47e8 shear=5.6e7 fric=20 coh=5.0e4 tension=1.0e4 ;体积、剪切、摩擦角、凝聚力、抗拉强度set mech ratio=1e-4solvesave Gravsol.savplot cont zdisp outl onplot cont szz;毛洞开挖计算initial xdisp=0 ydisp=0 zdisp=0model null range group insiderock any group concretliner anyplot block groupplot add axes redset mech ratio=5e-4solvesave Kaiwsol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;模筑衬砌计算model elas range group concretliner anyplot block groupplot add axes redini density 2500 range group concretliner any ;衬砌混凝土的密度prop bulk=16.67e9,shear=12.5e9 range group concretliner any ;衬砌混凝土的体积弹模、剪切弹模set mech ratio=1e-4solvesave zhihusol.savplot cont zdispplot cont sdispplot cont szzplot cont xzz;完成计算分析。
(完整版)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的单元为各向同性弹性本构。
5)命名一个范围(需要先命名这个范围)格式:range name <自己起一个名字> <范围> 示例:range name intersected_zones x 5 8 y 3 7 引用这个范围的格式:内容+range nrange <自己起的那个名字> 示例:range nrangeintersected_zone6)命名一个柱形范围格式:同4)示例1:range name bhole cylinder end10 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 not 7)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> ;绝对收敛标准,体系最大不平衡力小于某一临界值,此局限性较大。
FLAC3D各种命令笔记(2021年整理精品文档)

FLAC3D各种命令笔记编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望(FLAC3D各种命令笔记)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为FLAC3D各种命令笔记的全部内容。
hist reset ;清空里面所有的监测点nstress=-40e3 hist ramp的意思,很简单 nstress = —40e3 * ramp hist相当于乘法,为了方便控制动荷载的幅值。
plotitems--Add--strucktural elements——geometry ——cables 显示锚杆的几何形状Pl-add-vector-apply body force 显示施加的力,不过是在求解之后显示Creat 创建一新视图plotitems—--—--add——-—-vetor-—-—apply body force 显示荷载的施加位置Array var(4,3)二维数组 array var(5)一维5个元素的数组Var(m,n)=urand ;随机 var—数组的名字Pnt=z_next(pnt) 查找下一个单元体表Pl plas 显示塑性区Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格,可以给其赋值,以便调用.Array就是数组的意思。
Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移,ID=1材料参数dilation 剪胀角FISH语言:Degrade;∏/180Pi;∏Ngp;节点总数Nzone;单元体总数单元节点遍历:P_z = zone_headLoop while p_z # null………P_z = z_next(p_z)End loopRead(ar,n);将数组ar的前n个记录写到文件print zone state 显示塑性区域编写fish函数时,尽量采用浮点型,例如1。
Flac 3D 常用命令总结

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.基本命令:- !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--5.0常用命令集锦

建模1、调用文件:①文件与工程在同一个文件夹,只写文件名即可:If the called file is located in the same folder as the FLAC3D project file, then only the file name needbe entered with the CALL command.②不在同一个文件夹,全路径:Otherwise, the file may be called by specifying its complete path (e.g., c:\myfolder\file.dat).Undo ;撤销上一条命令2、创建旋转缩放视图3、建模命令model mech mohr ;莫尔库伦模型model mech elastic ;弹性模型set grav 0, 0, -9.81 ; 重力加速度 negative z-direction. (垂直向下!常用的)gen zone brick size 6,8,8 p0 -10, -10, -20 ... ;省略号表示写不下后面继续p1 10, -10, -20 ...p2 -10, 10, -20 ...p3 -10, -10, 0plot zonegen zone brick size 6,8,8 p0 -10, -10, -20 ... ;不规则六面体p1 10, -10, -20 p2 -10, 10, -20 ...p3 -10, -10, 0 p4 10, 10, -20 ...p5 -10, 10, 10 p6 10, -10, 0 ...p7 10, 10, 10plot current plot Plot01plot clearplot zoneUndo ;撤销命令set logfile 127X1001.txtset log on truncateset log offlist zone prin range x 0 1 y 0 1 z 0 1 ;显示指定范围内各单元的主应力,结果如下Hist命令:①命令编号按顺序从1开始:each history is numbered sequentially from 1 as it is entered via the HISTORY command.②查找显示所有的his命令:Return to the Flac3D> prompt and typelist histfor a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记录1次。
flac3d常用命令

1、最先需要掌握的命令有哪些?答:需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
2、怎样输出模型的后处理图?答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件。
3、怎样调用一个文件?答:File/call或者call命令4、如何施加面力?答:app nstress5、如何调整视图的大小、角度?答:综合使用x, y, z, m, Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键。
6、如何进行边界约束?答:fix x ran (约束的是速度,在初始情况下约束等效于位移约束)。
7、如何知道每个单元的ID?答:用鼠标双击单元的表面,可以知道单元的ID和坐标。
8、如何进行切片?答:plo set plane ori (点坐标) norm (法向矢量)plo con sz plane (显示z方向应力的切片)9、如何保存计算结果?答:save +文件名10、如何调用已保存的结果?答:rest +文件名;或者File / Restor11、如何暂停计算?答:Esc12、如何在程序中进行暂停,并可恢复计算?答:在命令中加入pause命令,用continue进行继续。
在我们分步求解中想得到某一个过程中的结果,不用等到全求完,还可以在分布求解错误的时候就进行改正,而不是等到结果出来。
13、如何跳过某个计算步?答:在计算中按空格键跳过本次计算,自动进入下一步14、Fish是什么东西?Fish是否一定要学?答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能。
Fish可以不用学,需要的时候查Mannual获得需要的变量就可以了。
15、FLAC3D允许的命令文件格式有哪些?答:无所谓,只要是文本文件,什么后缀都可以。
16、如何调用一些可选模块?答:config dyn (fluid, creep, cppudm)17、如何在圆柱体四周如何施加约束条件?可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
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
Flac3D命令--完整经典版(2)

Flac3D命令--完整经典版(2)预览说明:预览图片所展示的格式为文档的源格式展示,下载源文件没有水印,内容可编辑和复制实例分析命令:1. X ,Y ,Z 旋转 Shift+ X ,Y ,Z 反向旋转Gen zone ……;model ……;prop ……(材料参数);set grav 0,0,-9.81(重力加速度)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 1.jpg 图像输出(格式、像素、名称) plot set magf 1.0视图的放大倍数为1.0;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 = 10.0 定义depth 为10;end 结束对函数的定义;setup 运行函数setup8. (隧道生成)上部圆形放射性圆柱及下部块体单元体的建立,然后镜像。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
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或者什么作图软件绘制就行了。
26.关于dd、dip的理解:答:首先你要确定你定义面的法线方向,dip是倾角的意思(xoy面和定义面的夹角),法线向上dip 是锐角,法线向下dip时钝角。
dd是y轴顺时针转到法线在xoy面上的投影线时,所转过的角度。
27.如何调用FISH文件?答:新手调用FISH文件出错的原因是:文件名为fishcall.fis.txt(因为在记事本里这个txt的后缀是默认的)。
修改方法是:将fishcall.fis.txt里面的内容全部复制到新的记事本里,再按保存,保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis,类型也就正确了。
28.如何理解Tension-p shear-p tension-n shear-n表示的含义?答:在计算的循环里面,每个循环中,每个zone都依据failure criterion处于不同的状态,shear和tension 分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.29.绘制相关曲线命令的一些解释在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负);"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了30.如何查看剖面上任一点的位移值?答:print gp disp range x() y( ) z( )或print gp disp range idplot gp disp range x() y( ) z( )或plot gp disp range id31.a: fix x y z range x -0.10.1 b: fix x range x -0.1 0.1 请问a式与b式有何区别?答:a表示在x=0这个平面是固定的;b表示在x=0这个平面只x方向不能发生运动32.initial与apply有何区别?答:初始化命令,如初始化计算体的应力状态等;apply边界条件限制命令,如施加边界的力、位移等约束等。
initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件不会发生变化。
33.如何查看各个时段不平衡力的具体数值?答:采用his来记录计算,包括位移应力等命令his unbal his gp(zone) zdis range (0 0 0)或者id=?导出数据命令his write n vs m begin时步end时步+file filename.his要导出不平衡力的具体数值his unbal step 100000 or solvehis write 1 vs step begin 1 end 1000 file.123.his使用上述命令就可以查看各个时步下的不平衡力的具体数值34.如何显示变形轮廓线的命令?plo ske magf 10其中10为放大系数35、怎样开始入门?答:初步学习首先看说明书User's gurid的Section2 Gettingstarted和Section 3 Problemsloving with FLAC36、哪些命令需要掌握?答:grid gen prop interface ini fix set hist step solve pause cont sav等,在学习的过程中遇到这些命令和不会的命令应该查看说明书的Command Reference,里面讲的很清楚,特别需要注意的是一定要仔细看一下Command Reference开头的Optional Range Phrases,这是讲范围的,因为很多命令的施加都需要加到模型的一部分,这里就告诉你怎样来表示这一部分。
这一原则同样适合与其他5个ITASCA软件。
还有一点就是六个软件(FLAC FALC3D UDEC 3DECPFC PFC3D)的range大不相同37、FLAC和FLAC3D建模有什么区别?答:FLAC和FLAC3D基本的建模思想大致是一样的,不同的是FLAC是首先用GRID命令把整体的网格画出来,然后用MODELL NULL命令把网格分成许多部分,再用GEN命令把这些部分一块一块搭接起来。
然后再用GEN LINE,GEN ARC等命令画出要开挖的边界,在FLAC中产生网格的命令只有一个GRID,而且只能用一次。
而FLAC3D是直接用GEN命令一块一块搭接,GEN可以用很多次。
38、怎样查看模型?答:plot grid可以查看网格,plot grid num可以查看节点号。
39、怎样查看边界条件?答:plot fix40、怎样查看预先留设的开挖边界?答:plot grid mark。
41、怎么查看模型的云图?答:pl syy (sxx xdis ydis sig1..) fill,注意,带fill关键字就是查看云图,不带fill关键字就是查看等直线图。
42、怎样查看模型的矢量图?答:pl grid vel (disp) red43、怎样查看INTERFACE?答:pl grid iface id redprint iface id ; id是interface的编号44、FLAC怎么调整视图?答:ctrl+z或shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小,双击屏幕。
45、判断模型中网格点是一个还是两个的方法显示模型网格点的ID,如果有两个ID号,说明是两个。
节理面上的网格点应该有两个ID号。
46、改变模型应力边界时FLAC与3DEC的不同之处在3DEC中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界时,模型表面的应力为两者应力相加。
47、flac中影响运行时间的因素(1)单元材料,结构面之间的刚度对比很大;(2)划分单元尺寸对比很大,尽量避免长细比大于5:1的细长单元。
48 FLAC中有多少种材料模型内置11中本构模型(1)零模型用于开挖回填的区域(2)各向同性弹性模型(3)正交各向异性弹性模型(4)横观各向同性弹性模型(5)德鲁克-普拉格塑性模型(6)莫尔-库仑塑性模型一般岩土体的力学行为(7)节理化塑性模型用于研究材料破坏后行为(如边破稳定,屈服煤矿开采)(8)应变硬化/软化模型-库仑塑性模型薄板层状材料破坏后力学行为(9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖(10)双屈服塑性模型(11)修正的剑桥粘土模型49、对不平衡力的理解flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小,该值为百分数,且从不会减小到零,因此,模型不会达到绝对稳定的平横状态,可以根据要求的精度设定R值,模型默认的R值是1e-5。
注意R值小仅表示所有网格点的力都平衡,但模型还可能发生塑性流动,通过plot vel命令查看网格点速度,一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。
50、简单命令解释:plot add cont disp out on时时显示云图plot block group gp num on显示节点编号plot inter red ske 看接触面的效果.del range group 1 not不删除group信息del range group 2删除group信息ini z add数值大小range group name顺z方向移动建立接触面interface 1 face range x(范围)y(范围)z(范围)或者interface 1 face range plane orig(x y z坐标)norm(法线方向)water table clear清除水位线pl set back white ;设置背景set plot bitmap size (800,510) ;设置图片大小plot set caption size 37 ;设置图例数字大小(default=35范围10-50)plot set center 100 0 40fix x y z range x -0.1 0.1表示在x=0这个平面是固定的fixed boundaryfix x range x -0.1 0.1表示在x=0这个平面只x方向不能发生运动roller boundary51、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判断i Falfa=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执行什么定义的函数52.在固定面上施加位移时,fix与apply不能同时使用,应该用ini和fix,才能施加上。