flac怎样输出制定节点的应力值
flac 结构单元用法
（1）cable单元模拟锚杆（全长连接、非全长连接，预应力锚杆），土钉等
（2）beam单元模拟支撑、围檩、冠梁等
（3）liner单元模拟面层，地下连续墙
Fix x range x -0.1 0.1
Fix x range x 5.9 6.1
Fix y range y -0.1 0.1
Fix y range y 7.9 8.1
Fix z range z -0.1 0.1
set mech force=50
Gen zone brick size 6 8 8
Mode mohr
Prop bulk 1e8 shear 0.3e8 fric 35
Prop cohesion 1e10 tens 1e10 ;注意在此这个值给的很大。
Init dens 1000
Set gravity 0 0 -10
1、 得到初始应力的方法：
方法1、可以先给一些材料参数很大的值，进行初始求解，在计算之前再将材料参数设为正常值，即可。如在手册中给的第一个示例中就是这样做的。下面是例子，These are only initial values that are used during the development of gravitational stresses within the body. In effect, we are forcing the body to behave elastically during the development of the initial in-situ stress state.* This prevents any plastic yield during the initial loading phase of the analysis.
flac3应力-位移曲线
flac3应力-位移曲线
要查看FLAC3D中的应力-位移曲线，可以按照以下步骤进行操作：
打开FLAC3D软件并进入模型界面。
在模型界面中，选择需要查看的模型或网格。
在主菜单栏中，点击“Plot”选项卡。
在弹出的窗口中，选择“New Plot”选项。
在“Plot Type”中选择“Stress-Displacement”。
在“DataSource”中选择需要查看的节点或网格。
确保“X-Axis”选择的是“Displacement”或者“U”，而“Y-Axis”选择的是“Stress”。
点击“OK”按钮，FLAC3D将生成应力-位移曲线并显示在界面上。
通过以上步骤，您可以在FLAC3D中查看和分析应力-位移曲线，从而更好地了解模型的应力和位移分布情况。
FLAC若干问题的解答
FLAC若干问题的解答把结果中的信息输出到指定文件分为两种：1、把单元信息输出到指定文件set log onset log on finame.dat;可以输入文件路径，否则按当前目录处理print zoneset log off通过反复使用该命令，可以把不同信息输出到不同的文件。
2、把节点信息输出到指定文件set log onset log on finame.dat;可以输入文件路径，否则按当前目录处理print gpset log off请教一下print gp里的keyword该如何填。
比如想输出某点的zdis谢谢print gp dispprint gp positionlog 文件中有，然后处理，这些论坛里都有的关于flac3d—内置fish语言精讲FISH语言是FLAC3D程序的内置编程语言因为FLAC3D的最佳操作方式是命令流文件方式这一点与ANSYS很相似，而FISH就相当于ANSYS的APDL语言。
它包括循环、判断等结构。
如果你还用过其它高级语言，那么从形式上讲你也可以把它理解为子函数。
FISH语言的引入极大的方便了用户进行复杂的程序建模它不但可以嵌入命令流文件里工作而且还可以引用FLAC3D本身的任何命令所以说它实现了对FLAC3D的完全控制引入FISH语言的作法值得其它通用软件效仿。
以上仅是个人的使用体会请版主与管理员指教也供初学者入门之用。
谢谢！Fish函数增加了以下几种新的特性：1.增加了fish变量来获取结点、单元和界面变量2.fish提供了获取结构单元变量的途径3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量5.提供了fish绘图子程序函数够生用户定义的图形内容6.fish函数已经增加了从文件读、写数据的能力FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说也来凑个热闹，呵呵注意点：1.fish函数可以嵌套使用;2.以save命令保存模型时，fish函数和变量也同时保存；3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish 函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；4.变量或函数名不能以一个数字开头也不能是下列字符：. , * / ^ = > < # ( ) [ ] @ ; “ '5.如果用命令set safe on 指定了编译安全模式，则用户调用fish 函数时，函数名前必须加@；6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定；7.fish函数的调用方法：.可以出项在其他fish函数的单独行中；.可以出现在其他fish函数的表达式中；.出现在flac3d的命令行中；.作为命令set,print,hist的参数。
FLAC命令流
【常规显示】PLOT con syy ;表示Y方向应力云图PLOT con sxy ;表示XY方向的应力云图PLOT con szz ou on ;ou是outline的缩写，用来显示网格PLOT con szz ef on ;ef是effective的缩写，用来显示z方向有效应力的云图PLOT con szz inter 1e4 ; inter是interval的缩写，以改变云图显示中的增量大小，本例中默认为5.0E3PLOT con szz max -10e3 ;max表示显示的应力数值最大不超过-10e3，本例中默认为0 PLOT con szz min -40e3 ;min表示显示的应力数值最小值为-40e3，本例中默认为-60e3 PLOT bcon prop bu ;显示体积模量bulk的赋值块云图【面荷载作用下的变形】rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0 初始应变归零app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 施加面荷载solvesave 6-2.savPLOT con zd 分析面荷载作用下的沉降云图PLOT con zd ou on magf 20 变形后的网格图（20为放大倍数）PLOT con szz ou on magf 20 将变形后的网格与应力云图结合起来plot block state 显示塑性区域plot block state shear ;获得剪切屈服的单元，包括shear-n和shear-pplot block state tension-p ;获得过去拉伸屈服的单元plot block state now ;获得当前处于塑性状态的单元，包括shear-n和tension-nplot block state past ;获得过去处于塑性状态的单元，包括shear-p和tension-p【历史跟踪】rest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2hist id=2 gp zdis 1 1 3 记录此节点的Z方向位移，记录号为2hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 2hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solvesave 6-3.savPLOT hist 2 查看记录的2号记录PLOT hist 7 v 8 ;v是vs的缩写hist write 7 v 8 file 6-3hist.txt 使用该命令，程序会在默认文件夹中生成一个名为6-3hist.txt 的文本文件，读者可以打开这个文本文件观察所得的数据结果。
flac3应力-位移曲线
flac3应力-位移曲线Flac3 应力-位移曲线Flac3 是一种用于土壤力学分析的数值方法，用于模拟土体在加载过程中的应力-位移变化。
本文将从理论分析、图解解释和实例展示三个方面介绍 Flac3 中应力-位移曲线的特点和应用。
一、理论分析在理论分析部分，我们将首先介绍应力-位移曲线的一般形态。
在加载过程中，土体受到外部力的作用而发生变形，应力-位移曲线是描述土体变形特征的重要工具。
曲线通常呈现出刚性阶段、弹性阶段、塑性阶段等不同段落。
其次，我们将详细解释 Flac3 中应力-位移曲线的构建过程。
Flac3 是利用有限元方法对土体进行数值模拟的工具。
它通过离散化土体模型，计算节点处的应力和位移，并绘制应力-位移曲线。
具体而言，Flac3 在加载过程中会考虑土体的弹性参数、塑性指数和受力条件等因素。
二、图解解释在图解解释部分，我们将利用图表对应力-位移曲线进行详细说明。
通过绘制应力-位移曲线，我们可以直观地观察土体在加载过程中的变化。
我们将以某个工程项目中的岩土体模型为例，展示应力-位移曲线的实际应用。
通过对该模型的分析，我们可以清楚地看到初始时的刚性阶段、加载过程中的弹性阶段以及后续的塑性变形阶段。
这一图解解释将使读者更好地理解 Flac3 中应力-位移曲线的意义。
三、实例展示在实例展示部分，我们将提供一个具体的 Flac3 仿真案例。
通过分析该案例中的应力-位移曲线，我们可以深入了解 Flac3 在模拟土壤力学问题上的应用。
在该案例中，我们将研究一个基坑开挖施工过程。
通过加载不同的边界条件，比如土体初始应力和施工阶段的应力变化，我们可以观察到不同环境下的应力-位移曲线的变化。
这个实例展示将帮助读者更好地掌握 Flac3 在实际工程中的应用技巧。
综上所述，本文从理论分析、图解解释和实例展示三个方面介绍了Flac3 中应力-位移曲线的特点和应用。
通过深入探讨 Flac3 的工作原理和仿真案例，我们可以更好地理解和应用这一技术，为土壤力学分析提供有力的支持。
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个字符。
FLAC-3D小技巧
1、如果是外观上看模型，要切换视角的话除了上述的快捷键，命令是：plot set rotation 20 0 60 rotation后面的3个数字是绕x、y、z轴转动的角度。
使用快捷键时，每次操作转动的角度，或者move的大小，可以通过plot set rotincrement和plot set moveincrement来设置。
2、从surfer把数据倒入FLAC3D!不知道各位在做象边坡地形表面这种不规则面的模型时有没有什么好招。
我最近发现了一种方法。
在FLAC3D的fish函数库中有一个叫TOPO的函数，可以生成自由表面模型的网格，关键是它要用到的那个数据表做起来非常麻烦，具体就是在地形图上平行y轴画一系列的等距的直线，然后读出每条与每条直线相交的等高线的x值和高程，然后生成一张表，格式是table y x z。
用surfer可以很方便的生成这张表，在DATA菜单下的convert...可以将GRD文件生成以上要求的那种格式。
3、按找手册中的解释，print gp position 是用来输出gp的位置坐标的；print zone gp是用来输出与zone相关的gp的id值的。
这样我就有一个上面的疑问，两个命令中的gp都有id值，是不是不同命令中的gp的id值相同的话，他们的坐标就可以看作是一样的啊？答：如果要画等值线的话，先输出gp的坐标来，将其另存为一个文件，再输出zone的应力值，将与zone对应的gp的id与前面输出的gp的id值对应起来，找到gp的坐标。
用excel打开上面的文件进行编辑，最后用surfer画等值线。
这样的办法太苯，不知到有没有更好的办法啊？4、plot stensor plane local on out off能够显示的主应力矢量图，但仍有两个问题1）显示的矢量图中应力好像没有变化，如何通过设置符号的大小表示应力的大小变化？2）应力符号好像总在图形没有发生变形前的单元中心，如何使应力符号在发生变形了的单元中心？5、个人经验，无论有限元还是flac，如果支护措施仅仅只考虑其刚度的贡献，岩体的变形和塑性区变化往往都不大，可能需要适当提高岩体的参数，才能反映支护的作用6、＜锚杆支护的数值模拟＞．中国矿业大学岩石力学与土木工程学报前段时间有篇文章专门关于flac锚杆单元应用模型修改的文章，7、如何模拟预应力自由式锚索？首先我想说说我在模拟预应力自由式锚索的做法,然后在说出现的问题,希望大家帮忙解决一下!我采用先整根生成的方法来生成锚索,然后给自由段和锚固段分别赋不同的参数,自由段的灌浆参数统统赋为零,锚固段的按试验值确定.在模拟垫墩时,把自由段外端的交接删除,然后重新生成刚性铰接(之所以这样做是因为模型介质较软,不这样做就收敛不了,介质内发生了流变.而全长粘结式锚索不需这样处理也和试验值拟合的很好),使得该结点不会和周围介质发生相对位移;然后在垫墩作用的面积上再施加一个均布力,大小等于预应力除以垫墩面积.问题是位移值拟合锝还可以,但是锚索的预应力增加值比实测的大的多,有没有其他的方法来模拟自由式锚索???8、Fish函数增加了以下几种新的特性：1.增加了fish变量来获取结点、单元和界面变量2.fish提供了获取结构单元变量的途径3.休单元和面单元性质目前可以通过单元变量名z_prop(i_z string)和界面单元名i_prop(i_z string)分别加以识别4.fish函数可以获取单元应变和应力速率，还提供了全应变增量张量和应变速率张量5.提供了fish绘图子程序函数够生用户定义的图形内容6.fish函数已经增加了从文件读、写数据的能力FLAC3D是一个强大的软件，但是不得不承认的，它的界面没有ansys好用，商业运行不是很强的说注意点：1.fish函数可以嵌套使用;2.以save命令保存模型时，fish函数和变量也同时保存；3.fish函数不支持缩写，这与flac3d命令不同，另外所有的fish函数或变量不区分大小写，程序同意转化为大写进行编译，当然也可以通过执行set case-sensitivity on来区分大小写；4.变量或函数名不能以一个数字开头也不能是下列字符：. , * / ^ = > < # ( ) [ ] @ ; “ '5.如果用命令set safe on 指定了编译安全模式，则用户调用fish函数时，函数名前必须加@；6.如果变量不曾赋值，则系统默认为零（整形），如果赋值，其类型由值的类型决定；7.fish函数的调用方法：.可以出项在其他fish函数的单独行中；.可以出现在其他fish函数的表达式中；.出现在flac3d的命令行中；.作为命令set,print,hist的参数。
FLAC 使用步骤
步骤3：应用公式简略计算
应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。
步骤4：建立x, y坐标与node i, j之间的关系
FLAC 使用步骤
FLAC 使用步骤
1.FLAC程式使用前准备步骤
步骤1：依比例画出所欲分析的资料
在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。
步骤2：换算输入资料成同一单位
将现有地层资料，如 Density, Bulk modulus, Young`s modulus, tension, cohesion, friction Angle 等资料，换算成同一单位。
在图上距离及深度的关系，建立x, y坐标系统，再由x, y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应的x, y 坐标。
建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于1000，以节省时间。
(1)调出起初的应力平衡：re_____ .sav
(2)设定新的材料性质：model,prop
(3)设定新的支撑性质：struct
(4)设定新的外力 apply
(5)设定边界条件 apply,fix
(6)求工程时的应力平衡 solve, step, cycle
(7)储存 save file.sav
(2)设定材料性质：prop
(3)设定外力：Set Grav, Apply Pressure, ini sxx, Syy
flac2d常用图形命令
set log 11.txt
print pl syy ;pr syy line x1 y1 x2 y2 n
set log off
当从giic界面调用命令，网格不显示时，可以输入“！”。
把事先设好的监测点调出命令
plot pen history 5 line
line x1y1 代表第一点 x2y2代表第二点 n代表输出多少个点
x1 y1 为具体X、Y坐标而不是网格
set log on
set log 123.txt
print syy line 1，1 1,2 1
set log off
4、怎样查看模型？
答：plot grid 可以查看网格，plot grid num 可以查看节点号。
输出一条直线上y方向的应力
print syy line x1y1 x2y2 n
syy代表y方向的应力
line x1y1 代表第一点 x2y2代表
print sxx line x1y1 x2y2 n
sxx代表x方向的应力
8、怎样查看模型的矢量图？
答：pl grid vel (disp) red
9、怎样查看INTERFACE？
答：pl grid iface id red
print iface id ; id 是interface的编号
his write 8 vs 4 begin 600 end 800 123.txt
把监测点输出
his dump 1 3 7 vs 2 begin 150 end 375
;shu chu xian shang yingli
set pline x1 y1 x2 y2 ;ju ti dian bu shi i j dian
flac常用命令
flac3d常用命令1、最先需要掌握的命令有哪些？答：需要掌握 ge n, in i, ap p, pl o, so lve 等建模、初始条件、边界条件、后处理和求解的命令。
2、怎样输出模型的后处理图？答：Fil e/Pri nt ty pe/Jp g fil e，然后选择 Fil e/Pri nt，将保存格式选择为 jpe文件。
3、怎样调用一个文件？答：F ile/c all 或者 cal l 命令4、如何施加面力？答：app nstr ess 5、如何调整视图的大小、角度？答：综合使用 x,y, z, m, S hift键，配合使用Ctr l+R，C trl+Z等快捷键。
6、如何进行边界约束？答：fixx ran（约束的是速度，在初始情况下约束等效于位移约束）。
7、如何知道每个单元的ID？答：用鼠标双击单元的表面，可以知道单元的ID 和坐标。
8、如何进行切片？答：plo s et pl ane o ri (点坐标) n orm (法向矢量) plocon s z pla ne (显示 z 方向应力的切片) 9、如何保存计算结果？答：sav e +文件名 10、如何调用已保存的结果？答：r est +文件名；或者 Fil e / R estor 11、如何暂停计算？答：E sc 12、如何在程序中进行暂停，并可恢复计算？答：在命令中加入p ause命令，用conti nue 进行继续。
在我们分步求解中想得到某一个过程中的结果，不用等到全求完，还可以在分布求解错误的时候就进行改正，而不是等到结果出来。
13、如何跳过某个计算步？答：在计算中按空格键跳过本次计算，自动进入下一步14、Fis h 是什么东西？Fi sh 是否一定要学？答：是FLAC3D 的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能。
Flac学习心得
Flac&3dec学习心得快捷键—显示位移等值线：pl cont disp快捷键—显示塑性区：pl bl st快捷键—显示网格：pl grid显示某根锚杆的应力：pl sel cable stress range id xx在显示等直线等时可以用modify中的outline显示单元的边界平移模型的命令：ini x add 30 (x方向平移30m)在三维计算中，想看某一剖面的图，可以用edit菜单中的cutting plane操作任何操作如果不想对某个区域起作用，可以用not关键词，比如gr 5 not gr 6 not 改善锚杆作用的一种操作—把锚杆单元的link删掉，再重新见建一个link与周围的单元连接起来，如下所示：sel del link range id 1;sel del link range id 11 删掉link11;sel del link range id 12;sel del link range id 22;sel link 2 11 target zone 在原来的link 处（此处的2是该处的node号）新建link2 与单元连接;sel link net rang id 11;sel link net rang id 12;sel link 13 22 target zone;sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 11 any id 1 any id 12 any id 22三种方式模拟预应力锚索：方式1、通过删除－建立link模拟。
删除掉锚索端头原来CABLE和ZONE之间自动建立的连接，然后建立钢性连接。
示例命令如下：sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &gr_coh 10e5 gr_k 2e7 range cid 11,17sel delete link range id 1sel link id=100 1 target zonesel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100sel cable id=1 pretension 60e3 range cid 1,10方式2、通过赋三段属性模拟。
abaqus种输出某个节点的应力曲线
abaqus种输出某个节点的应力曲线
1.打开Abaqus/CAE，加载您的模型。
2. 找到您要输出应力的节点。
您可以使用“Find”命令来查找节点。
3. 单击该节点，然后选择“Output Requests”菜单。
4. 在“Output Requests”对话框中，选择“Stress”选项卡。
5. 您可以选择输出哪些应力分量。
如果您想输出全部应力分量，可以选择“All Components”选项。
6. 点击“OK”按钮。
现在您已经设置了应力输出请求。
7. 在Abaqus/Viewer中，选择“Field Output”菜单。
8. 在“Field Output”对话框中，选择“Stress”选项卡。
9. 选择您刚刚设置的应力输出请求。
10. 点击“OK”按钮。
现在您已经成功输出了该节点的应力曲线。
11. 您可以使用“Curve”菜单来修改和导出该曲线。
- 1 -。
Flac3D命令--完整经典版
实例分析命令：1. X ，Y ，Z 旋转 Shift+ X ，Y ，Z 反向旋转Gen zone ……；model ……；prop ……（材料参数）；set grav 0,0,（重力加速度） 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 图像输出（格式、像素、名称） plot set magf 视图的放大倍数为；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 = 定义depth 为10；end 结束对函数的定义；setup 运行函数setup8. （隧道生成）上部圆形放射性圆柱及下部块体单元体的建立，然后镜像。
flac使用说明
FLAC程序简介FLAC是FAST LAGRANGIAN ANALYSIS OF CONTINUA的缩写,是由美国明尼苏达ITASCA软件公司开发的通用程序。
该程序在中国大陆以外已有较多的用户，应用很普遍。
该程序刚引进国内，目前国内尚在推广应用。
FLAC程序的基本原理和算法与离散元法相近，是由P.A．Cundall提出的。
它与离散元法的区别在于它应用了节点位移连续的条件，可用于连续介质的大变形分析。
由于它不必形成像有限元法中那样的整体刚度矩阵，因此可以在内存较小的微机上计算较大规模的题目。
例如对于4M内存的微机可运行大约15000个单元的题目。
FIAC程序可以模拟弹性模型材料，摩尔—库仑模型材料，横观各向同性模型的层状材料，具有软弱夹层的节理材料等六种。
它还可以模拟地应力场的生成、洞室或边坡开挖、回填混凝土、锚杆锚索安设、地下渗流等。
尤其是对锚杆的设置非常方便，可以在任何指定位置设置锚杆而不考虑网格的划分和结点的分布。
FLAC程序的另一特点是，它具有强大的前后处理功能。
网格自动生成，界面美观。
用户可以使用各种命令修正网格以适应各种复杂边界，计算结果均可以有图形输出，并可着色。
这包括各期的主应力分布向量、σx、σy、τxy分布等值线，位移向量，Ux、Uy等值线，塑性区范围，锚杆受力等等。
使用方便快速。
FLAC可以按两种方式运行，既可以通过数据／命令文件运行，又可通过人机对话方式运行。
用户可以在FIAC运行中的任何时候中断它，修改数据后继续运行。
FLAC2.25使用说明书FLAC 的输入和一般的数值模拟的程序不一样, 它可以用交互的方式从键盘输入各个命令, 也可以写成命令文件, 类似于批处理, 由文件来驱动。
FLAC 命令大小写一样。
所有的命令可以附带若干个关键词和有关的数值。
在下面的命令解释中, 只有大写的字母起作用, 小写的字母写不写、写多少个都没有关系。
i,j,m 和n 开始的变量要求整型数, 否则要求实型数。
FLAC教程
FLAC特点
综上可知FLAC有以下特点：
(1)采用混合离散化方法模拟塑性破裂与塑性流动，比 采用归纳积分法更合理；
(2)采用全动态运动方程使FLAC在处理不稳定问题时不 会遇到数值困难；
(3)采用显式解法，在求解非线性应力—应变关系时， 不需要存储任何矩阵进行修改，与普通隐式解法相比， 大大节约了 机时； (4)FLAC按行与列(而不是按顺序)的形式进行单元编号， 这对于某些指定单元的研究很方便；
2.2 FLAC软件的基本原理
FLAC软件的基本原理即是拉格朗日差分法，是一种利 用拖带坐标系分析大变形问题的数值方法。并利用差 分格式按时步积分求解。 拉格朗日元法源于流体力学，在流体力学中有两种主 要的研究方法。一种是定点观察法，亦称欧拉法；另 一种是随机观察法，称为拉格朗日法。后者是研究每 个流体质点随时间而变化的状态，即研究一物体质点 在一段时间内的运动轨迹、速度、压力等特征。
FLAC中可以模拟的模型
(1)零空模型(Mull)：代表网格中的孔洞(开挖单元)； (2)应变硬化/软化模型(SS)：代表非线性，不可逆剪切 破碎与压缩； (3)粘弹性蠕变模型； (4)界面模型(界面为平面，沿界面允许滑动和分开)： 模拟断层、节理和摩擦边界； (5)水利模型：模拟可变形空隙体与粘性流动的全藕荷； (6)结构单元模型：模拟岩土体加固、衬砌、锚杆、混 凝土喷层、可缩支柱及钢拱等。
i=21
Call
文件名
写成的命令文件可以用 Call 命令来调用 , 命令文件的最后一行必须是 RETURN, 以返回到交互方 式。命令文件中不能有 CALL 命令本身。 Config 关键词 FLAC 用以解平面应变问题,但经过配臵命令也 可以用于解平面应力问题或轴对称问题。需要时应在 形成网格之前发。关键词有: P_STR AX CYC n 该命令同 STEP, 为执行 n 个时步的循环运算。 平面应力问题 轴对称问题
Flac3D中文手册
Flac3D中文手册Flac3D 中文手册FLAC3D的计算模式中是否需要做孔压分析取决于是否采用config fluid命令。
1 无渗流模式（不使用config fluid）即使不使用命令config fluid，仍然可以在节点上施加孔压。
这种模式下，孔压将保持为常量。
如果采用塑性本构模型的话，材料的破坏将由有效应力状态来控制。
节点上的孔压分布可由initial pp命令或water table命令来设定。
如果采用water table命令，由程序自动计算水位线以下的静水孔压分布。
此时，必须施加流体密度（water density）和重力（set gravity）。
流体密度值和水位位置可以用命令print water显示。
如果水位线是由face关键字来定义的，则可用命令plot water命令显示水位。
这两种情况，单元的孔压都由节点孔压值平均求出，并在本构模型计算中用作有效应力。
这种计算模式下，体积力中不反映流体的出现：用户必须根据水位线以上或以下相应地指定干密度和湿密度。
使用命令print gp pp和priint zone pp可分别得到节点或单元孔压。
plot contour pp命令可绘出节点孔压云图。
2 渗流模式（使用config fluid）如果使用命令config fluid，则可进行瞬时渗流分析，孔压改变和潜水面的改变都可能出现。
在config fluid模式下，有效应力计算（静态孔压分布）和非排水计算均被执行。
除此之外，还可进行全耦合分析，这种情况下，孔压改变将使固体产生变形，同时体积应变反过来影响孔压的变化。
如果采用渗流模式，单元孔压仍由节点孔压平均求出。
但这种模式，用户只能指定干密度（不论是水位以上还是以下），因为FLAC3D 将流体的影响考虑到了体积力的计算中。
采用渗流模式时，渗流模型必须施加到单元上，使用命令model fl_isotropic模拟各向同性渗流，model fl_anisotropic模拟各向异性渗流，model fl_null模拟非渗透物质。
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的单元为各向同性弹性本构。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
怎样输出制定节点的应力值我计算一个模型有几万个节点，我怎么才能把指定的节点应力输出。
指定节点的应力？到底是节点的力嘛，还是单元的应力哦？两个办法：1. 先plot block model id on，然后找到你想看的那个单元的id，假设是10，然后就plot con szz out on range id 10；2. 用fish，通过gp_head，gp_near，gp_id等等获得你要的节点，然后可以得到gridpoint的各种信息。
通过flac的fish语言，写节点信息和相关的力到文件既然，ansys的格式文件能够导入到flac文件，那么用fish肯定能够导出相关的数据，然后再看处理结果。
如题，我建议，通过flac的fish语言，写节点信息和相关的力，请帮主支持。
如果能够把相关的信息导入到文件，那么就很容易的对数据进行处理。
我先贴一个：用fish写文件的，当然也是从这个网络上找到的，我解释一下，也不算直接抄袭。
;; Initialization;初始化def ini_mesh2tecIO_READ = 0;定义读文件关闭吧？！！IO_WRITE = 1;定义写文件IO_FISH = 0;定义一个？？？IO_ASCII = 1;定义用asci码iarray buf(1);定义一个数组，tec_file = 'tec10.dat';打开文件的名字zone_ngp = z_numgp(zone_head);zone 的标志数endini_mesh2tec;; Write Tecplot File Headdef write_mesh_headbuf(1) = 'TITLE = "Flac3d Mesh to Tecplot Version 10" \n'buf(1) = buf(1) + 'VARIABLES = "X" "Y" "Z" "XDISP" "YDISP" "ZDISP"\n' buf(1) = buf(1) + 'ZONE T="Tecplot v10" \n'buf(1) = buf(1) + ' N=' + string(ngp) + ','buf(1) = buf(1) + ' E=' + string(nzone) + ','if zone_ngp = 4 thenbuf(1) = buf(1) + ' ZONETYPE=FETETRAHEDRON \n'elsebuf(1) = buf(1) + ' ZONETYPE=FEBrick \n'endifbuf(1) = buf(1) + ' DATAPACKING=POINT \n'buf(1) = buf(1) + ' DT=(SINGLE SINGLE SINGLE SINGLE SINGLE SINGLE)' status = write(buf,1)end;; Write Grid Point Coordinates and Displacements;写节点信息def write_gpp_gp = gp_headloop while p_gp # null;这个循环，叫做遍历所有节点buf(1) = string(gp_xpos(p_gp))buf(1) = buf(1) + ' ' + string(gp_ypos(p_gp))buf(1) = buf(1) + ' ' + string(gp_zpos(p_gp))buf(1) = buf(1) + ' ' + string(gp_xdisp(p_gp))buf(1) = buf(1) + ' ' + string(gp_ydisp(p_gp))buf(1) = buf(1) + ' ' + string(gp_zdisp(p_gp));这就是传说中的写status = write(buf,1)p_gp = gp_next(p_gp)endloopend;; Write Zone Connectivity;写单元信息def write_zonep_z = zone_head;获得单元的最高号码loop while p_z # null;循环，或者叫做遍历buf(1) = ''if zone_ngp = 4 thenbuf(1) = buf(1) + string(gp_id(z_gp(p_z, 1))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 2))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 3))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 4))) + ' 'elsebuf(1) = buf(1) + string(gp_id(z_gp(p_z, 1))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 2))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 5))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 3))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 4))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 7))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 8))) + ' 'buf(1) = buf(1) + string(gp_id(z_gp(p_z, 6))) + ' 'endifstatus = write(buf,1)p_z = z_next(p_z)endloopend;; Main Functiondef mesh2tecstatus = closestatus = open(tec_file,IO_WRITE,IO_ASCII)if status = 0 thenwrite_mesh_headwrite_gpwrite_zonestatus = close;关闭文件ii = out('Successfully Write Data Into File ' + tec_file);显示信息elseii=out('Open File Error! Status = ' + string(status))endifendmesh2tec;调用把所有的记录文件写到文件中，名字不一样，那出来给大家共享，如果谁知道谁自动获得最大的记录的编号，那么这个程序就真正完美了，不过，可以凑合一下哦！！！文件名是循环生成的哦，最大数是300，请使用的时候，修改这个值 loop iixx_ (1,300)上面的300，是根据自己的项目变化的哦！！！def write_logloop iixx_ (1,300)logfile11= 'log' + string ( iixx_ ) + '.his'commandhis write iixx_ file logfile11end_commandend_loopendwrite_log自动生成的，乖乖，好多哦！！！！一个一个的看吧！！！一下是异想天开几点：1.一个文件写n行，现在是一个文件写2行2.直接导入到excel，然后就作图3.想办法和上面的那个loop结合起来，导入到能够画等值线的程序中，如ansys，显示结果！！！set log onpr histset log off然后算一下有多少行就知道有多少纪录了如何寻求cable中的应力值试试plot add sel cable grout可以通过直径显示应力大小，加颜色可以表示出应力方向补充一下plot add sel cable grout stress首先谢谢ying111您所说的是，从图像上确定cable的受力性状；如何确定cable中沿轴线方向任何一点的具体应力值，能否说的更清楚一些！请赐教。
可以试试：hist sel cable id=n grout stress x,y,z请问大虾们，print gp disp出来的数据如何保存吗？？使用命令中的文件菜单显示记录功能。
ｌｏｇfile设置log打开刚才试了一下好象不行：仅给出这样一个文件：;****************************************;Log File Started 13:17:48 Thu Apr 15 2004;Using FLAC3D version 2.10-239;By: Itasca Consulting Group, Inc.; Minneapolis, MN USA;Job:Flac3D>;Log File Ended 13:18:47 Thu Apr 15 2004;****************************************不知道问题出在什么地方。
最好的办法是：hist write 1（或你前面记录的任何因素,作为纵坐标；最多可以写12个因素）vs （或你前面记录的任何因素,作为横坐标）（skip 步数）file 你定义的文件名。
这样可以任意输出你需要点的或区域的位移历史或应力历史的数据，保存在记事本文件中很方便应用的。
比print 好多了！log file配合list命令FLAC的缺省计算模式就是平面应变问题.如要进行平面应力分析,设定命令: config P_stressunderstand it?。