用Tecplot对Flac3D计算结果进行后处理
FLAC3D知识介绍
FLAC 3D基础知识介绍一、概述FLAC（Fast Lagrangian Analysis of Continua）由美国Itasca公司开发的。
目前，FLAC有二维和三维计算程序两个版本，二维计算程序V3.0以前的为DOS版本，V2.5版本仅仅能够使用计算机的基本内存64K），所以，程序求解的最大结点数仅限于2000个以内。
1995年，FLAC2D已升级为V3.3的版本，其程序能够使用护展内存。
因此，大大发护展了计算规模。
FLAC3D是一个三维有限差分程序，目前已发展到V3.0版本。
FLAC3D的输入和一般的数值分析程序不同，它可以用交互的方式，从键盘输入各种命令，也可以写成命令（集）文件，类似于批处理，由文件来驱动。
因此，采用FLAC程序进行计算，必须了解各种命令关键词的功能，然后，按照计算顺序，将命令按先后，依次排列，形成可以完成一定计算任务的命令文件。
FLAC3D是二维的有限差分程序FLAC2D的护展，能够进行土质、岩石和其它材料的三维结构受力特性模拟和塑性流动分析。
调整三维网格中的多面体单元来拟合实际的结构。
单元材料可采用线性或非线性本构模型，在外力作用下，当材料发生屈服流动后，网格能够相应发生变形和移动（大变形模式）。
FLAC3D采用的显式拉格朗日算法和混合-离散分区技术，能够非常准确的模拟材料的塑性破坏和流动。
由于无须形成刚度矩阵，因此，基于较小内存空间就能够求解大范围的三维问题。
三维快速拉格朗日法是一种基于三维显式有限差分法的数值分析方法，它可以模拟岩土或其他材料的三维力学行为。
三维快速拉格朗日分析将计算区域划分为若干四面体单元，每个单元在给定的边界条件下遵循指定的线性或非线性本构关系，如果单元应力使得材料屈服或产生塑性流动，则单元网格可以随着材料的变形而变形，这就是所谓的拉格朗日算法，这种算法非常适合于模拟大变形问题。
三维快速拉格朗日分析采用了显式有限差分格式来求解场的控制微分方程，并应用了混合单元离散模型，可以准确地模拟材料的屈服、塑性流动、软化直至大变形，尤其在材料的弹塑性分析、大变形分析以及模拟施工过程等领域有其独到的优点。
FLAC3d后处理步骤
1、首先rest123.sav，然后再callFLAC3d2Tecplop生成Tec10.dat2、打开Tecplot程序，执行【File】/【Load Data File】命令，打开数据文件对话框，选择刚刚生成的tec10.dat文件，并执行【打开】命令，在出现的Select Initial Plot 对话框中，单击OK按钮。
并在Tecplot左上角显示的3维笛卡尔坐标系下选择【3D Cartasion 】3、在Tecplot中选中【Data】→【Extract】→【Slice from Plane】，打开切面对话框。
在对话框中我们可以选中x=6/y=5/z=7等切片，并单击Extract.4、在Tecplot中选中【File】→【White Data File】，打开【White Data FileOptions】对话框，在对话框中选中ASCII、Point。
并在下面的Zone（s）选项中选中你所要的2 SLc=3，在Variable(s)中同时选中1 x（m）,2 y(m)及你所要的13 szz（pa）等，并单击OK按钮。
保存的文件自己最好独自命名（如tec-90.dat格式）5、以txt格式打开刚才命名的tec-90.dat，把里面的三列数复制到txt里面后再导入Excel里面，同时在Excel里面把szz的值变成相反的值保存在D列，（txt里面注意不是三列的都不要，只复制三列数据，把前面的和后面的都去掉）。
6、打开surfer，【网格】→【数据】，选中刚才生成的123.xls(即刚才生成的Excel数据，选中相应的A、B和D列),并单击打开。
对于【网格化数据】选项单则单击【确定】，则会生成与123.xls相对应的123.grd（即与刚才生成的Excel 文件名相同）。
7、在【地图】→【表面图】单击，在【打开网格对话框】中选中生成的123.grd文件，并单击打开。
8、右击图片选中【属性】并打击打开，可以对网格、颜色进行设置（颜色谱选安装文件夹中的rainbow.CLR即可,具体是在打开的【属性对话框】→【常规】→【材质颜色】→【上部】单击上部中的颜色，在颜色谱中选择【载入】，在打开里面选中安装文件夹里面的rainbow.CLR即可，在多次确认）9、此外也可以在surfer里面对其他参数进行设置。
深部大断面硐室围岩变形及控制技术
深部大断面硐室围岩变形及控制技术苏晓建【摘要】采用数值模拟方法分析了赵固二矿－800 m泵房围岩变形原因及控制技术。
认为深部大断面硐室围岩应力集中区距离硐室中心较远，但底鼓较为明显，支护方案设计时应重点关注。
基于此提出了主动支护与被动支护联合支护方案，首先锚网喷支护，其次双底拱刚性支架联合支护，第三围岩注浆加固。
工程应用表明，锚注＋刚性支架支护能有效控制围岩变形，满足泵房保持长期稳定的要求。
%Deformation mechanism and control technology of surrounding rocks of -800 m water pump house in Zhaogu No. 2 Coal Mine were analyzed by using numerical simulation. It was thought that the distance from stress concentration area to the center of chamber was rela-tively far,so the obvious floor heave should be the key point of supporting scheme. Based on the condition,combined supporting scheme of initiative and passive support was proposed which in-cluded bolt-mesh-spurting supporting,rigid support with double bottom arch and grouting rein-forcement. The result showed bolting-grouting and rigid support could effectively control the surrounding rock deformation and meet the requirement of chamber's long-term stability.【期刊名称】《中国煤炭》【年(卷),期】2017(043)001【总页数】4页(P75-78)【关键词】深部硐室;大断面;围岩控制;联合支护;数值模拟【作者】苏晓建【作者单位】中国矿业大学北京资源与安全工程学院，北京市海淀区，100083; 河南能源化工集团焦煤公司赵固二矿，河南省辉县市，453633【正文语种】中文【中图分类】TD353赵固二矿Ⅰ盘区－800 m泵房布置在二1煤层下方20 m左右的砂质泥岩中,埋深约900 m。
06-边坡变形-flac3d前后处理
• ini sxx=-5.0e6
• fix x range x -0.1 0.1
• fix x range x 5.9 6.1 • pause
• pl con sxx • solve • pl con sxx
• 注意是约束x方向
第三十八页，共84页。
第三十二页，共84页。
赋予初始应力
• gen zone brick size 6 6 6 • model elas • prop bulk 5e9 shear 3e9
• ini sxx=-5e6
• pause
• solve
只用ini，所有单元将有sxx应力， 但solve变为0，因为没有约束， 力引起位移最后平衡.
。 体）和Tetrahedral（四面体）单元
第二十四页，共84页。
第二十五页，共84页。
在ansys选择菜 单:
Input from
选择APDL命令流 文件(在仿真论坛 可以下载此程序)
执行后在电脑的D:\ 下生成ansys单元节 点信息的两个文件
01_node.dat
02_ele.dat
第二十六页，共84页。
• apply sxx=-5e6 range x 5.9 6.1 • pause
• pl con sxx • solve
• pl con sxx
• 注意是约束x方向---------变形参数并不影响应力场，它只影响位移场。如果考虑重力的时候， 变形参数影响泊松比，从而影响侧压力系数，才会影响水平的应力场分布。
平衡力大于容许值时，才会进行cyc。设置不同的变形参数并不影响位移结果，因为系统平衡了。 如果人为设置执行计算流程（如cyc 10000），由于位移边界条件约束，边界上位移为0，内 部则会慢慢随时间产生位移，但应力场是不会发生变化的。
Tecplot画等值线的操作步骤详细介绍00
FLAC3D计算结果—Tecplot画位移等值线，比较详细的操作步骤介绍！！！主要内容：详细介绍在FLAC3D中计算完毕后，把计算结果导入Tecplot中，绘制位移等值线的方法。
第一部分：Tecplot画位移等值线的具体操作步骤。
1. 运行FLAC3D程序，先调用FLAC3D计算结果的保存文件xxxx.sav，命令为：restore xxxx.sav然后，调用dynamax 编制的flac3d到tecplot的转换程序Flac3d2T ecplot.dat，命令为：call Flac3d2Tecplot.dat则在flac3d计算结果文件的当前目录下，生成一个名为tec10.dat文件，这个文件就是将flac3d计算结果转换成符合tecplot 10的数据格式后的文件，可供tecplot 10调用。
2. 在tecplot 10中，调用第一步生成的tec10.dat文件，将FLAC3D计算结果导入到tecplot 中，具体操作为：运行tecplot 10程序，点击File主菜单，在下拉菜单中选择Load Data File(s)选项，在弹出的Load Data File对话框中，选中上一步生成的tec10.dat文件，点击“打开”，之后一路点击默认选项，得到下图，3.），显示模型内部的一个指定剖面。
（相后，图形显示只选择on plane）4. 点击Iinteractive slicing tool*弹出如下对话框，显示切片工具的详细设置，先勾选Show Slices选项（显示切片）；再选择切片的方向（X-Planes和Y-Planes等分别表示法线方向平行于X轴和Y轴等方向的切片）；然后指定切片的位置，在Position后面输入数字，该数字为沿着所选切片的法线方向上的切片位置坐标。
切片设置好后，点击Close，关闭设置对话框。
5. 将切片截面上的数据读出，形成除整体zone外的另一格zone，操作如下，之后弹出如下对话框，点击Extract，弹出如下信息对话框，显示切片上数据提取成功。
tecplot生成等高线的步骤
1、先用FLAC把结果计算出来，可以保存几张图片（zdisp,xdisp,ydisp等），再接着Call“Flac3d2Tecplot位移应力.dat”这个文件，则生成了“tec10.dat”这个文件。
2、打开Tecplot10.3这个程序，点File>>Load data files选择“tec10.dat”文件。
在左边位置的zone下的“Contour”勾上、“Boundary”勾上,zone effects下的“Lighting”勾上,其余的不选。
3、“Plot”>>contour出现对话框，点“Show Lables”，可将等高线的数字显示出来，同时可设置字体的颜色和格式“Color,Font”
3、主菜单的“Data”>>Extract>>Slice from plane，出现对话框，在对话框中“Slice Plane”下的“constant”选择例如“X=0.5”就是X=0.5这个平面。
点“Extract”
4、点在左边的“Zone style”出现一个框，选择第一个“Tecplot v10”，点“zone show”选择“Deactive”；只让需要的面击活，其余的都再点后面的Deactive ，选择需要的那个面，在“Contour Type”选择“Line”，再点后面的“Line color”选择线的颜色。
5、点“Plot”>>Axis出现对话框，把“Show X-Axis”前的勾取消（X、Y、Z 三个都要取消）
6、点“View” >>“Rotate”出现对话框，点“Preset views”可选择视图的方向，点“Fit to full size”则可将图形填充在整个区域。
flac3d后处理
;Sk是草图sketch的缩写，作用使显示模型的外围网格线~
; ··变形矢量图··
;pl sk dis
;变形矢量图可以形象的表征变形的发生方向和相对大小
;pl sk dis scale 0.07
;scale关键字用来调整矢量线头的的大小
;pl block state
;显示塑性区域
;可在pl bl state 添加关键字
;pl bl state shear ; 获得剪切屈服的单元
;pl bl sta tension-p; 获得过去拉伸屈服的单元
;pl bl state now tension-n; 获得当前处于塑性状态的单元
;pl bl state past tension-p; 获得过去处于塑性状态的单元
new
gen zon b size 3 3 3
mod mohr
prop sh 1e6 coh 10e3 fric 15 bu 3e6
fix z ra z -.1 .1
fix x ra x -.1 .1
fix x ra x 2.9 3.1
fix y ra y -.1 .1
;
;save 6-2.sav
;施加额外的边界条件，
;为了更形象的描述模型变形后的结果，常常使用变形后的网格图命令如下
;··变形云图··
;pl con zd ou on magf 20
;magf 为放大倍数，不能缩写为mag，倍数自定，pl bcon不可以跟随 magf
;PL SK MAGF 20
;pl con szz max -10e3
;max 表示显示的应力数值最大不超过-10e3，默认为0
“用tecplot对flac3d计算结果进行后处理”相关讨论帖整理(第二部分)
第二部分：对dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”及相关讨论帖的整理和总结。
1. dynamax大牛的原创主题帖“【原创】用tecplot对flac3d计算结果进行后处理”的链接地址为：SimWe仿真论坛» D01：FLAC3D/FLAC2D » 【原创】用tecplot对flac3d计算结果进行后处理说明：该flac3d-tecplot转换程序，由dynamax大人利用fish编写，总共有三个版本，不断改进。
最初只能导入位移数据，画出位移等值线，并且不能剔除挖空单元（null），只能导入全部网格；第二版不再导出null zone，并且可以指定显示范围；第三版则可进一步同时导出位移和应力结果。
dynamax发表于2004-7-11 17:45[前后处理] 【原创】用tecplot对flac3d计算结果进行后处理fish写的因为flac3d中不能标注等值线图将位移计算结果导入tecplot中进行处理;可以绘出漂亮的等值线图熟悉tecplot的可以试一下不然得到转换后的数据文件还是不知道如何处理给出了简单的测试实例283356-tecplot.rar (45.85 KB) （第一版转换程序）dynamax切片功能，从菜单中选择，操作为：data = >extract = >slice from planedynamax上图显示成2D Cartesian可能更好要显示数值的话，plot-> Contour-> More-> Labels-> Show Labels坐标轴可以关掉plot-> axis-> show %-axisdynamax我给的例子只能导出全部网格的信息(针对第一版转换程序而言)如果要只导出开挖后剩下的部分因为节点编号的关系要另外做处理才行否则没有办法做出像下面这样的图我看代码有点不对，稍加修改应该没有问题，主要是把null部分的网格信息排除我们努力看看行不行。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用Tecplot对Flac3D计算结果进行后处理
第一步：用FLAC3D打开已经完成的*.sav文件；
第二步：然后【File】/【Call】运行Flac3d2Tecplot.dat文件，运行结束时会在命令行里显示写入成功，在命令文件所在的文件夹里会发现tec10.dat数据文件（可以根据自己的需要更改文件名）。
第三步：打开Tecplot程序，执行【File】/【Load Data Files(s)】命令打开tec10.dat数据文件。
以下是画位移等值线图的步骤：
1.三维图形选择3D Cartesian，二维选择2D Cartesian；
2.在边栏左上方Zone Layers栏目中,去掉Mesh前的√，选中Contour前的√。
然后按下
more按钮，在Var选择自己需要显示的量：
3.选择切片工具，并打开下面的Tool Details，可以设置切片的方向和位置等；
4.将切片截面上的数据读出，形成整体zone外的另一个zone；
5.点击zone style，编辑zone；
6.点击surfaces，发现有两个zone，在zoneshow中激活Slice：z=0，关掉tecplot v10；
7.点击contour，contour type选取lines；
8.然后再等值线上显示标注；
9.然后可以进行去掉颜色、网格等编辑。