flac3d 常用命令师姐整理。
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常用命令
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、怎样查看模型?答: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 5.0常用命令集锦
建 模 【1 】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 ). Undo ;撤消上一条敕令2.创建 扭转 缩放视图3.建模敕令model mech mohr ;莫尔库伦模子model mech elastic ;弹性模子set grav 0, 0, -9.81 ; 重力加快度 negative z -direction. (垂直向下!经常应用的)下下面面这这代代码码,,是是沿沿着着--y y 偏偏向向的的重重力力加加快快度度,,留留意意差差别别!!!!!!!!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 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 type list hist for a listing of the histories and their corresponding numbers.hist nstep 5 ;每5步记载1次.默认是10步记载1次hist add ratio ;记载maximum unbalanced force ratiohist add gp zdisp 4 4 8 ;记载节点的位移model null range x=2,4 y=2,6 z=5,10 ;To excavate the trench,enterThe dimensions of the excavation region (i.e., the null material) are defined by the x-, y- and z-coordinate range.set large ;设置大变形large-strain logic is specifiedini xdis 0 ydis 0 zdis 0 ;位移归零solvesave trench ; 保持盘算成果的文件“TRENCH.F3SAV”不须要后缀4.group敕令(留意新版本的差别 ,group zone)gen zone brick size 6 6 6group zone Tunnel range x 1 6 y 0 6 z 1 5 ;设置groupmodel mech mohr ; 莫尔库伦模子model mech null rangegroup Tunnel ;group zone Tunnel range cylinder end1 100 105 100 end2 100 95 100 radius 3group zone Matrix range group Tunnel not; naming RANGES that are fixed in space and unrelated; to the model geometry - represents stratigraphy...range name Layer1 plane dip 0 dd 0 ori 0 0 0 aboverange name Layer2 plane dip 0 dd 0 ori 0 0 0 belowprop bulk 1e8 shear 0.5e8 coh 0 tens 0 fric 35 range nrange Layer1prop bulk 1e7 shear 0.3e7 coh 1e7 tens 0 fric 0 range nrange Layer2pl zone colorby group4.Interface敕令以下是切割模子法生成interface:寺河模仿用的pl zonegen zone brick p0 0 0 0 p1 200 0 0 p2 0 40 0 p3 0 0 4 size 100 20 2 group 1 ;gen zone brick p0 0 0 4 p1 200 0 4 p2 0 40 4 p3 0 0 8 size 100 20 2 group 2 ;gen separate group 1interface 1 wrap first group 1 second group 2显示interface:plot create plot Plot01plot add surface;显示概况plot add interface red;显示交互面为红色plot show用导来导去法,已经在5.0上运行过,很好,值得借鉴pl zone colorby groupgen zone brick size 3 3 3 &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)group zone Basegen zone brick size 3 3 3 &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)group zone Top range group Base notpl add interfacepauseini zpos add -1.0 range group Top4.hist敕令(5.0版本须要加add)hist nstep=5 //设置采取距离为5 等效于set hist_rep=5 hist add id =2 zone szz 0,0,-1hist unbal //最大非均衡力hist add gp zdisp 4 4 8print histhist add gp zdis 0,0,0hist add gp zvel 0,0,0hist add gp ydisp 0,0,0;监测网格点(0,0,0)点的y 偏向上的位移hist add zone syy 0,1,0;监测单元(0,1,0)点的y 偏向上的应力hist add zone syy 1,1,0;监测单元(1,1,0)点的y 偏向上的应力hist add zone szz 0,0,-1List historyset largehistory nstep 2history add fish tthistory add fish thetahistory add zone sxx (1,1,1)history add zone syy (1,1,1)history add zone szz (1,1,1)plot add his 3 4 5 vs 1HI STORY add fish var;fish var变量的监测5.solve敕令什么意思5.plot敕令经由过程图形操纵界面的方法:pl zone colorby gr ;检讨各个分组经由过程敕令流的方法:plot create plot Plot01plot clearplot add history 1 linestyle color blackplot post plot Plot01 file 127F1002.ps manualplot create plot Plot02plot clearplot add history 2 linestyle color blackplot post plot Plot02 file 127F1003.ps manualplot create plot GravVplot cut add plane dip 90 dd 0 origin 3,4.1,0 ;截面plot set center 3 3 3plot set eye 12 18 12plot add zone plane behindplane on onplane off ;addlabel "Zone" white plot add bcont szz ramp greyplot add axesplot post plot GravV file 127F1004.ps manualres kw50plot create plot szz_contourplot cut add plane dip 90 dd 180 origin 0,18.9,0plot set orient 120 220 0plot set center 45 45 0plot add zone plane behindplane on onplane off ;plot add bcont szz ramp rainbowplot add axesplot showplot create GravVplot set plane dip=90 dd=0 origin=3,4,0plot set rot 15 0 20 ;this would be achieved interactivelyplot set center 2.5 4.2 4.0 ;and is show here simply to illustrateplot add bound behindplot add bcont szz planeplot add axesplot show如许我们就树立了一个叫“GravV”的图形,将它设为当前图层,对此图形设置一个倾角为90°的面,(在x-y面上,商定z的负偏向为下),偏向为0°.(在x-y面上,从y的正偏向开端沿顺时针偏向计量),过x-y面上一点(x=3,y=4,z=0),在平面后加一个网格,垂直偏向的应力Qzz,最后给模子加轴线,整体图形轮廓与添加图形相反,显示每个区域中间的应力盘算值,每个区域的色彩直接对应根本区域的应力,当键入show,我们树立的图形被显示Plot create DispcontPl add cont disp plane behindPl add axesPl showplot sketch;绘出轮廓线plot add vel;添加矢量plot add his 1 2 vs 3;绘出汗青记载n,actLoad , p_solcyc 15000plotCreate TrenchAdd surface yellowAdd axes blackShowplot create szz_contour;创建文件名为szz_contour(z 偏向的应力云图)plot set cent 4 4 5;视图的中间坐标为(4,4,5)plot set rot 20 0 30;视图的扭转角度为(20,0,30)plot add cont szz out on shade on ;打开szz 的等值线图plot add axes;打开坐标轴线plot show;显示窗口;不合的区域,材料赋值不合pl zone colorby modelpl zone colorby group求解前对模子进行检讨网格检讨(如前所述)鸿沟前提检讨速度束缚前提 plot gpfix荷载前提 plot fap模子检讨模子赋值 plo block model参数赋值 plo block prop ***初始应力检讨(如后所述)经常检讨模子的响应plo zcon szz (syy, sxx)应力场plo con zdis (ydis, xdis)位移场plo (add) zone colorby state屈从状况plo gpfix 速度束缚前提plo fap 体力plo hist (unbal) 1不服衡力(必须的)pl bitmap view 1 file myplot.pngpl bitmap view 1 size 1024 768pl create <viewid>pl current <viewid>pl cut add plane dd xxx dip xxxpl cut add plane name xxx normal xxx origin xxx pl dxf <viewid> <file filename.dxf>pl excel <file filename.xlsx>pl svg <viewid> <file filename.svg>pl set background white ;设置布风景plot add his 3 4 5 vs 15.网格的衔接.对称……不合尺寸的zone,网格用attach衔接,粘起来6.plane敕令6.range敕令list rangerange name trench x=-1,1 y=0,4 z=-2,2;界说rangemodel mech null range nrange trench; ; 开挖excavate the trench7.gp敕令list gp position;检讨节点编号.地位坐标print gp disp range id 105hist add gp zdis 0,0,0hist add gp zvel 0,0,0message(1)=string(gp_xpos(pnt))message(1)=message(1)+','+string(gp_zpos(pnt))message(1)=message(1)+','+string(abs(gp_zdisp(pnt)))GP PNT find gp(INT)GP PNT gp find(INT)GP PNT gp headii = ngp total number of gridpoints更多节点变量函数,详见fish reference8.zone敕令list zoneprint zone stress range id 2381 any id 2361 any 2341 any ;输出网格单元的应力hist zone szz 0,0,-1Z PNT find zone(INT)INT z find(INT)Z PNT zone headii = nzone total number of zones更多zone变量函数,详见fish reference9.config9.鸿沟前提Boundary and Initial ConditionsBoundary conditions are normally specified with the APPLY and FIX/FREE commands, and initial conditions are assigned with the INITIAL commandExample: a constant, compressive xx-stress component of 10 MPa can be applied to a boundary located at x = 10 with the command apply sxx -10e6 range x 10; x=10的面上 10MPaapply syy = -20e6 grad 0,0,20e5 range y -20 z 0,10在斜面上施加法向力10 杂项(材料力学性质)在进行模仿时,采取摩尔-库仑(Mohr-Coulomb)强度准则来断定矿体与顶底板的损坏机理.该准则根本的内在为:岩石的剪切损坏是由岩石决裂面产生的损坏剪应力受到岩石材料的内聚力和内摩擦力的抵抗,当决裂面上的损坏剪应力大于岩石的内聚力和内摩擦力时,即 Fs>0 时,材料产生剪切损坏[4]:有限差分法(FDM)Fish函数1. if前提从句2. loop前提从句newdef derives_mod = y_mod / (2.0 * (1.0 + p_ratio))b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio))end@derivelist @b_mod @s_modgen zone brick size (2,2,2)model mech elasticprop bulk=@b_mod shear=@s_modlist zone prop bulklist zone prop shearnewdef xxxsum = 0prod = 1loop n (1,10)sum = sum + nprod = prod * nend_loopend@xxxlist @sum @prodnewdef abcif xx > 0 thenabc = 33elseabc = 11end_ifendset @xx=1list @abcset @xx=-1list @abcnewgen zone brick size 10 3 5def place_cablesloop n (1,5)commandsel cable beg (0.0,1.5,@z_d) end (7.0,1.5,@z_d) nseg 7 end_commandend_loopend@place_cablesplot add zone trans 75plot add sel geomnewgen zone brick size 10 3 5 model mech mohrprop bulk 1e8 shear 0.3e8 fric 35 prop coh 1e3 tens 1e3ini dens 1000set grav (0,0,-10); ; 重力加快度set largehist add unbalsolvesave cab_strini xdis 0 ydis 0 zdis 0hist add gp xdisp (0,1,5)def place_cablesloop n (1,5)z_d = 5.5 - float(n)commandfree x range x 0.0 z @z_b,@z_tsolvesel cable beg (0.0,0.5,@z_d) end (7.0,0.5,@z_d) nseg 7 sel cable beg (0.0,1.5,@z_d) end (7.0,1.5,@z_d) nseg 7 sel cable beg (0.0,2.5,@z_d) end (7.0,2.5,@z_d) nseg 7 sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 ... end_commandend_loopend@place_cablessave cab_endplot zone trans 75 sel cabblock force line width 4newgen zone brick size 1 2 1model mech mohrprop shear=1e8 bulk=2e8 ...cohes=1e5 tens=1e10def 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)end@get_addef loadload=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) endhist add fish loadhist add gp ydis (0,6,0)step 1000plot his 1 vs -2newgen zone brick size 10 3 5model mech mohrprop bulk 1e8 shear 0.3e8 fric 35prop coh 1e3 tens 1e3ini dens 1000set grav (0,0,-10); ; 重力加快度set large; ;大变形hist add unbalsolvesave cab_str; ;保持盘算成果文件,不必加后缀名.Sav ini xdis 0 ydis 0 zdis 0; ; 位移归零hist add gp xdisp (0,1,5)def place_cablesloop n (1,5)z_d = 5.5 - float(n)commandfree x range x 0.0 z @z_b,@z_tsolvesel cable beg (0.0,0.5,@z_d) end (7.0,0.5,@z_d) nseg 7 sel cable beg (0.0,1.5,@z_d) end (7.0,1.5,@z_d) nseg 7 sel cable beg (0.0,2.5,@z_d) end (7.0,2.5,@z_d) nseg 7 sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 ... end_commandend_loopend@place_cablessave cab_endplot zone trans 75 sel cabblock force line width 4def parm global rad=4.0 ; radius of spherical cavityglobal len=10.0 ; length of outer box edgeglobal in_size=6 ; number of zones along outer cube edgeglobal rad_size=10 ; number of zones in radial direction end@parm gen zone radbrick edge @len size @in_size @in_size @in_size @rad_size &rat 1.0 1.0 1.0 1.2 dim @rad @rad @raddef mountainlocal gp = gp_headloop while gp # nulllocal zz = sqrt(gp_xpos(gp)^2 + gp_ypos(gp)^2)local dz = 0.06 * sin(0.2 * zz + 100.0) ; Sum Fourier terms fordz = dz + 0.06 * sin(0.22 * zz - 20.3) ; quasi-random surfacedz = dz - 0.04 * sin(0.33 * zz + 33.3) ; topology.gp_zpos(gp) = 0.5 * gp_zpos(gp) * (1.0 + dz)gp = gp_next(gp)end_loopend@mountainfix x range x -.1 .1fix y range y -.1 .1fix z range z -.1 .1set grav 0,0,-10ini density=2000ini szz=-2.0e6 (grad 0,0,2.0e4)ini sxx=-4.0e6 (grad 0,0,4.0e4) syy=-4.0e6 (grad 0,0,4.0e4)solvedefine afill ; fill matrix with random numbersarray var(4,3)local mloop m (1,4)local nloop n (1,3)var(m,n) = urandend_loopend_loopendnewdefine h_var_1local zi = z_near(1,1,1)h_var_1 = z_sxx(zi)h_var_2 = z_syy(zi)h_var_3 = z_szz(zi)endThe FLAC3D commands to request histories might be hist add fish h var 1hist add fish h var 2hist add fish h var 3degrade第31页,共31页。
FLAC3D命令流(整理版)
1、怎样查看模型?答:plot grid 可以查看网格,plot grid num 可以查看节点号。
2、请问在圆柱体四周如何施加约束条件?答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
命令流如下:fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后,如何显示点的坐标?答:使用plo blo gro gpnum on命令5、关于gauss_dev对性质进行高斯正态分布的问题?答:根据手册上的说明:下面的命令设定一个平均摩擦角为40度,标准方差是±5%。
则命令如下:prop friction 40 gauss_dev 2 问题:请问gauss_dev 2中的2是如何计算的?如果把±5%改为±10%,则命令应如何写?40×5%=26、 reflect问题问:gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill plot surfgen zone reflect dip 0 dd 90 (对xy面做镜像)gen zone reflect dip 90 dd 90 (对yz面做镜像)(1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线。
FLAC3D命令流(整理版)
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常见命令与使用技巧
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常见命令与使用技巧
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常用命令总结
1、模型建好如何检验所建模型是否有问题注意:第一步建模一定不要建太大的模型,可以建非常小的模型,确保模型的准确性,这样可以尽快开始后边的步骤,等后边计算分析步骤快速做完验证没有大问题,再将准确的模型建好换上,重新计算得到更加准确的结果。
2、如何检验应力边界、位移边界、加载速率等是否准确Plot add gpfix:看位移边界Plo add fap:看应力边界,看加载方向是否准确Plot add vel:看速率,加载速率是否准确3、如何查看运算结果Plo con szz outline on:看z向应力云图,并显示边界网格Plo add disp:看位移Plo con zdisp outline on magf 5:看z向位移云图,变形扩大5倍Plo add bo :添加边界Plo blo state:看是否屈服Plo add stensor:查看大主应力方向Plo con smin:小主应力Plo con smax:大主应力Plo con ssi(shear strain increment):查看剪切应变增量云图Plo bcont smin:查看单元小主应力4、如何做切面并查看切面结果在编辑“Edit“菜单中设置切面位置,在上述3命令后加上plane(pl)即可。
如:Plo con smin pl:查看切面小主应力5、后处理将模型导入到tecplot中做切面,加等值线等,对结果进一步处理,得到更漂亮的一些结果。
其他重要命令总结:1、如何在圆柱体四周如何施加约束条件?可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
2、如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer 或者什么作图软件绘制就行了。
(完整版)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> ;绝对收敛标准,体系最大不平衡力小于某一临界值,此局限性较大。
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命令流(整理版)
1、怎样查看模型?答:plot grid 可以查看网格,plot grid num 可以查看节点号。
2、请问在圆柱体四周如何施加约束条件?答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
命令流如下:fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 13、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。
4、用命令建立模型后,如何显示点的坐标?答:使用plo blo gro gpnum on命令5、关于gauss_dev对性质进行高斯正态分布的问题?答:根据手册上的说明:下面的命令设定一个平均摩擦角为40度,标准方差是±5%。
则命令如下:prop friction 40 gauss_dev 2 问题:请问gauss_dev 2中的2是如何计算的?如果把±5%改为±10%,则命令应如何写?40×5%=26、 reflect问题问:gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill plot surfgen zone reflect dip 0 dd 90 (对xy面做镜像)gen zone reflect dip 90 dd 90 (对yz面做镜像)(1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线。
flac3d命令总结
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与FLAC常见命令的整理
FLAC/FLAC3D常规问题的整理1.FLAC3D命令的FAQlakewater整理看到其它板块上都有这个FAQ,也就是常见问题问答,今天抽了时间进行了整理,想到了就写下来了,因为看到很多初学者费了很多的时间,但是还是没有将常用的命令掌握,所以这个也可以作为入门的初级教材,使大家能够快速的上手,而不用为了某个小命令到处求助。
1. FLAC3D是有限元程序吗?答:不是!是有限差分法。
2. 最先需要掌握的命令有哪些?答:需要掌握gen, ini, app, plo, solve等建模、初始条件、边界条件、后处理和求解的命令。
3. 怎样看模型的样子?答:plo blo gro可以看到不同的group的颜色分布4. 怎样看模型的边界情况?答:plo gpfix red sk5. 怎样看模型的体力分布?答:plo fap red sk6. 怎样看模型的云图?答:位移:plo con dis (xdis, ydis, zdis)应力:plo con sz (sy, sx, sxy, syz, sxz)7. 怎样看模型的矢量图?答:plo dis (xdis, ydis, zdis)8. 怎样看模型有多少单元、节点?答:plo 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 / Restor17. 如何暂停计算?答:Esc18. 如何在程序中进行暂停,并可恢复计算?答:在命令中加入pause命令,用continue进行继续19. 如何跳过某个计算步?答:在计算中按空格键跳过本次计算,自动进入下一步20. Fish是什么东西?答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学?答:可以不用,需要的时候查Mannual获得需要的变量就可以了22. FLAC3D允许的命令文件格式有哪些?答:无所谓,只要是文本文件,什么后缀都可以23. 如何调用一些可选模块?答:config dyn (fluid, creep, cppudm)后注:这个工作很繁琐,需要的时间很多,希望广大网友能够将自己曾经遇到的常见问题在后续跟贴,也为了将这个FAQ进行很好的充实。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
pl sel geom cid on ra y 22 24 ;监测衬砌内力,找衬砌cid
pl sel geom cid on node off ra y 21 23 ;监测桩号,找桩cid
pl sel rec sres my surfx 0 1 0 ra y 0 24 scale=0.05
pl add fap ra cyl end1 0 0 -17 end2 0 42 -17 r 3.1
pl add sel geom id on
;查局部坐标
pr sel node local ra id=1
pl bl state shear-n ;看屈服
pl add gpfix ;查固定边界条件
pl add su magf 10 ;变形后位移
pr gp disp ra id=12067;查看某点位移
pl bcon prop bu;无刚度错误错误
pr sel node link ra id ;找node连接的link号
pr sel link tar ra id ;找link的目标节点
pr sel link sou ra id ;找link的源节点
pr sel link attach
pr sel link c lind k;查link的刚度
pr sel rec sr my
;pl sel disp magf 100
;pl sel pile fo fx
;pl sel pile mo my
;pl con zdis
pr sel rec surf
pl sel geom node off nodesys on scale=0.06
pl sel geom node off surfsys on scale=0.06
pl sel pile fo fx; 看桩轴力
pl bl gr id on ra y 20 21;查看zoneID