FLAC3D各种命令笔记

histreset;清空里面所有的监测点nstress=-40e3histramp的意思，很简单nstress=-40e3*ramphist相当于乘法，为了方便控制动荷载的幅值。plotitems——Add——struckturalelements——geometry——cables显示锚杆的几何形状Pl-add-vector-applybodyforce显示施加的力，不过是在求解之后显示Creat创建一新视图plotitems------add-----vetor----applybodyforce显示荷载的施加位置Arrayvar(4,3)二维数组arrayvar(5)一维5个元素的数组Var(m,n)=urand;随机var—数组的名字Pnt=z_next(pnt)查找下一个单元体表Plplas显示塑性区Arrayx_b(x,y)z_b(x,y)创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。Histid=1gpxd700监测700点处的x方向位移，ID=1材料参数dilation剪胀角FISH语言：Degrade;∏/180Pi;∏Ngp;节点总数Nzone;单元体总数单元节点遍历：P_z=zone_headLoopwhilep_z#null………P_z=z_next(p_z)EndloopRead(ar，n)；将数组ar的前n个记录写到文件printzonestate显示塑性区域编写fish函数时，尽量采用浮点型，例如1.0而不是1Z_sig1z_sig2z_sig3fish中的三个主应力，由大到小Printgppoisonrangidn输出节点n的坐标Xtable(n,s)对ID号为n的表的第s行，x列进行赋值，table表只有x,y两列Pltable1line/mark/both关键词描述曲线的类型，线、标记、二者综合File_name=’add_step’+string(n)+’.sav’Gp_near(x,y,z)得到靠近坐标x,y,z的节点地址Gp_zdis(p_gp)地址为p_gp的节点的z向变形String(n)将变量n转换成字符串格式Dim看基本网格里面的标准attachface常用命令，建模型之后连接面sxx=syy=(λ/1-λ)szzszz=γHSelfixxyz(局部坐标系位移速度)/xrotyrotzrot(局部坐标系转动位移)Inixdispydispzdisp(全局坐标系位移)/xpos(节点全局坐标)/xrdisp(全局坐标系旋转位移)Xv(全局坐标系位移速度)/xrvel(全局坐标系旋转速度)如果边界是FIX控制，则自动平衡，不用apply控制应力。对地下深部开挖的工作，一般忽略顶部到底部的地心引力的差别，setgrav被省略全局坐标系xyz局部坐标系d（倾向）s（走向）n（法向）plofaprangrouptunnel;显示施加在组上的体力applyremove施加新的边界条件前，用以移除旧条件Tablenx1y1x2y2x3y3建成一个点的坐标的表,以供FLAC引用，n为是表号。材料参数附加关键字：Gauss_devs以输入数值为平均数，s为偏差数，高斯正态分布计算Gradientgxgygz渐变：值=输入值+gx*x+gt*y+gz*zUniform_devs以输入数值为平均数，s为偏差数，随机正态分布计算hisunbalsetmechforce50检测最大不平衡力到50Mnranggroupxx1x2开挖Plotgrid显示网格Inizadd-6ranggrouppile*pile组整体在z方向平移-6个单位建立网格：实体多少米p0p1p2…就设多少米，单元数，在百个左右，各方向与实体成比例Brick六面块体redbrick放射状cylinder柱体网格radcylinder渐变放射Radtunnel六面体隧道渐变放射cshell柱形壳体网格Genzonredbrickp0…p8sizen1n2n3n4ratior1r2r3r4Modelmohr赋予本构模型Propbunshearncohnfricn赋予材料性质弹性只赋值弹性、剪切模量Fix固定边界Inidensn初始化密度App施加力Reflect镜像命令genzonreflectnormxyzorixyzGrad渐变iniszz-30e3grad0020e3ranz01应力ZZ方向的分量按（-30e3+0*x+0*y+20e3*z）渐变plotblogro显示分组plotblockgrouprangegroupdam只显示dam组单元plotitems_blockgroup_modify_勾选gridpointid’s_ok显示节点ID编号Histgp会自动监测距离输入坐标最近的节点Histzone会自动监测距离输入坐标最近的单元Plothist7vs8SetHist_rapn监测间隔，越小，曲线越圆滑histWrite7vs8file6-3hist.txt输出对比信息文件print输出结果printgpdisidrang517any输出节点517的变形值记录文件：log文件是FLAC输出文件的格式之一setlogon建立一个名为flac3d.log的文件，记录窗口中所以信息，直到setlogoff继续下一个log文件，setlogfile6-1.logsetlogonprintgpzdisprangid2179aid2189setlogoffplotconzd/zv显示z方向位移/速度plotsxx显示应力等值阴影图（没有contour）plotconzdouonmagf20附带网格线，变形量放大20plotsk显示网格线plotskdisscale0.07网格上附加矢量图，scale默认0.05越大，箭头越长后处理：plotsetrot20030;输出前设置视图角度plotskdisscale矢量图显示剖面显示：plotsetplaneori01.50norm010设定定一个垂直于Y轴的剖面，ori为切面上一点坐标，norm为切面的法向（即从01.50引到010连线为法线）pladdcondispplanebehindshadeon剖面后面的区域一样显示plotconzdplane在设定的剖面上显示z方向位移plotaddske增加网格线add命令意思是在绘图区域增加plotadddisplane增加位移矢量plotaddaxe增加坐标系动画制作：Plotsetrot20030设置视角Setmovieavistep1file6-5.avi格式设置Moviestar动画开关Solve求解完成后，自动生成文件。初始地应力场的生成：弹性求解法，分阶段的弹塑性求解法存在了null模型，不能采用solveelastic的求解方法获得初始应力，所以采用分阶段的弹塑性求解方法。先将Mohr模型的凝聚力。俏和伉拉强度护赋俏为无穷人进行求解，保证在重力作用卜单儿不至十发生屈服，然后冉将Mohr模型参数赋值为真实值，再进行求解。接触面：分界面应附在小面区域，与密度大的区域。胶合模型：不允许张开和滑移库伦剪切模型：抗拉粘结模型：粘结型和滑移型。Bslipeoff不滑动（默认）on滑动接触面几何模型建立：移来移去法：首一先将需要建立接触而的两个网格分开建立，然后在一个网格的指定位置建立接触而，然后把另外一部分网格移到特定的位置，这种方法称为移来移去法。要点就是把需要建立接触而的网格分开建立，然后在相应的网格上建立接触而，再把没有建立接触而的网格移动过来，这样做的目的就是使同一坐标点有不同的节点号，从而可以模拟结构物的滑移、分离等情祝。Genzonradcylp0---p11sizea（x）b（y）c（弧度上的单元数）d（放射方向的网格）ratio1111.15interface1face（在指定范围内所有网格的面上，均创建分界面单元）rangecylend1000end200-5.1radius.31&cylend1000end200-5.1radius.29notradius（半径），后一个半径要小于前一个，即两个圆柱的差，得到边界。命令指的整个圆柱再抠掉后面的End圆柱体断面中心点的坐标,两个点确定一条线，再根据半径确定柱体。Not就是否认前面的意思Plselgeom.显示创建的sel单元Inizadd-6ranggrouppile*pile组整体在z方向平移-6个单位导来导去法：1、保存完整模型保存为1.Sav2、保留需要建接触面的分组，删掉其他所有单元，用到命令：delranggroup2not，在分组上建立接触面，保存为2.sav；3、rest1.sav&delranggroup2&*导出删除后的模型expgrid1.flac3d；4、rest2.sav&impgrid1.flac3dInterface1propkn(法向刚度)ks(切向刚度)cfkskn由实验确定，一般取接触面两侧较小的E的10倍*注意语法切割模型法：此方法是在己建好的模型上，把group1与周围组分开，然后在group1上建立接触面。Genseparategroup1*分离出group1建面Int1wrap12*int=interface在group1group2之间建立接触单元。结构单元：ID结构单元编号CID结构构件编号node节点link结构连接的编号（结构单元间、结构与网格间）SelbeambegendnsegSelbeampropSelpileid=1cid=ibeg000end0010nseg=4*建立桩单元，设立结构单元ID和结构构件ID，结构构件4个，三个节点Selnodeid=*设置节点IDPlselgeoidonnodonscale0.04Selshellselid=Irangcylend1end2radiusnot*在半径处建一个无厚度的面Selnodeinizposadd-25*壳单元的位移，节点初始位置减少25Delete删除实体单元seldelete删除结构单元plotaddselgeomnode=offshrinkfac=0magfac=10redred显示结构单元shrinkfac是指结构单元的尺寸，范围是0-1，默认值为0.1桩单元模拟锚杆支护特性：selpileproprockbolton可以激活锚杆特性，能够计算支护四周的约束应力、锚杆与周围单元之间的应变软化特性以及锚杆拉断、拉裂等现象。Printselnode?Plotselnode?*后处理其中fob表示不平衡力的大小，local表示结构节点的局部坐标系锚杆：三种非全长锚固模拟方法1、通过删除锚索端头，即Cable端部的node和zone之间自动建立的连接，然后在他们之间建立刚性链接来模拟托盘。锚杆(锚索)自由段和锚固段通过设置不同的锚固剂参数来模拟，预紧力加在锚杆自由段。seldellinkrangid1*删除自动链接sellinkid=1001targetzone*在cable单元与zone之间设置新链接target目标sellinkattachxdir=rigidydir=rigidzdir=rigid&*六个自由度rigid刚性的xrdir=rigidyrdir=rigidzrdir=rigidrangid=100*设置链接为刚性selnodefixxyzxryrzrranid=1*固定六个方向的自由度（此句为桩pile）selcableid=1pretension60e3rangcid110*自由段施加预紧力2、将锚杆(锚索)的端头、自由段、锚固段赋不同的属性来模拟非全长锚固预应力锚杆〔锚索)，端头的锚固参数设为极大来模拟托盘，这样在锚杆〔锚索)受力时，端头将不会滑动，相当于托盘的作用。预紧力加在锚杆(锚索)自由段。3、借助别的结构单元模拟托盘