FLAC3D 3.0 C 常用命令

合集下载

flac3d 常用命令师姐整理。

flac3d 常用命令师姐整理。
pl sur gp on ra y 20 24 ;监测点位移,找点gpid
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

Flac3D 常用命令流

Flac3D 常用命令流

FLAC3D 常用命令流;模型镜像gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一点);沿X轴镜像,通过对称平面法线向量确定对称面gen zone reflect normal 0 0 -1 ;沿z轴镜像;绘图控制pl <add> contour szz outline on;在模型中显示位移-应变曲线hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror";对模型进行压缩实验的方法;即在模型两侧施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1;修改模型的坐标值ini x add -100 y add -100 z add -100;显示云图的同时也显示模型网格轮廓plot add cont disp outline on <gradient/average>;gradient更精确;输入角度、弧度方法pi=π,90°为90.0*degraddef set_valsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint ptA ptB ptC;施加结构单元方法sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) &end2=(0.0,25.0,0.0) radius=24.5 notplot add sel geom black black cid on scale=0.03sel node init zpos add -25.0;如何显示某一平面plot create name_planeplot set plane origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom blackplot add axes redplot current name_planeplot show;如何施加法向应力并显示所施加的力apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力plot grid green fap black;显示施加的力;获取某网格点的ID,坐标值def aaget_id=gp_id(gp_near(3,0,0))endaaprint get_idprint gp position range id 3;保存命令流文件set log onset logfile hua.log;文件存放的路径与调用的文本路径相同......set log off;fish函数的说明与解释def find_add ;定义fish函数find_addhead = null ;给head赋值p_gp = gp_head ;第一个网格结点的指针赋给p_gploop_while p_gp # null ;当p_gp值不为null时作循环x_pos = gp_xpos(p_gp) ;将指针为p_gp的结点的x坐标值赋给x_posif x_pos = 10.0 then ;如果x_pos = 10.0 则(执行)new = get_mem(2) ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址mem(new) = head ;将head类型和数值置于地址为new的fish变量mem(new+1)=p_gp ;第二个变量存储符合条件的zone地址head = new ;将new值赋给headendifp_gp = gp_next(p_gp) ;将结点指针为p_gp的下一个结点的指针赋给p_gp endloop;结束循环end ;结束fish函数;实际上这个fish函数为满足条件(x坐标为10的)的zone的地址开辟一定的地址空间,各地址之间存在一定的联系;找到第一个符合条件的zone地址后,用new = get_mem(2) 从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址,并用;下面的mem(new) = head ,将第一次开辟的两个变量的第一个变量存储地址head(注意第一次head=null),第二个变量存储第一;个符合条件的zone地址,并将第一个变量的地址赋予head(head=new),第一次循环结束;下次循环,同样开辟两个变量对象,第;一个变量对象记录上次循环开辟的第一个变量的地址,第二个记录第二个符合条件的zone地址,其余循环依次类推,这样子就建立了;一个符合条件的zone地址链条,方便以后使用。

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 end2 radius r2 not,其中r225.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。

(完整word版)FLAC3D命令流

(完整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常用命令总结第一篇: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常用命令

(完整版)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常用命令师姐整理。

flac3d常用命令师姐整理。

flac3d常用命令师姐整理。

pl sel rec sres my surfx 0 1 0 ra y 0 24 scale=0.05pl add fap ra cyl end1 0 0 -17 end2 0 42 -17 r 3.1pl add sel geom id on;查局部坐标pr sel node local ra id=1pr sel rec surfpl sel geom node off nodesys on scale=0.06pl sel geom node off surfsys on scale=0.06pl sel pile fo fx; 看桩轴力pl bl gr id on ra y 20 21;查看zoneIDpl sur gp on ra y 20 24 ;监测点位移,找点gpidpl sel geom cid on ra y 22 24 ;监测衬砌内力,找衬砌cidpl sel geom cid on node off ra y 21 23 ;监测桩号,找桩cidpl 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 attachpr 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。

FLAC3D与FLAC常见命令的整理

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. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

gen zone brick gro 1 ratio 1.1 1 1 size 15 1 10 p0 0 0 19.5 p1 20 0 19.5 p2 0 1 19.5 p3 0 0 29.5 建模.
mo null range group 7 开挖
attach face检查节点
gen merge 1e-3融合节点
plo blo gro;按组画图
plo blo den;按密度画图
plo con zdis;Z向位移画图
plo con sz;Z向应力画图
set mech ratio 1e-3;设置计算收敛
plo con ssi disp vel ;剪应变增加滑移线
pl add disp;位移矢量图合成位移
plo sel recoer;可以显示应力
ini y add 2
在地壳每增加100m深度,压力将增加约2.75Mpa
his gp disp id 1851
pl con disp;
his write 7 vs step begin 1 end 10000000000 file 顶板.txt;输出记录数据
pl bl st n;
pl bl st
plo fap red sk;显示体力
plotitems------add-----vetor----apply body force
plo add ve apply body force
ini szz 1000 range z 13 14 x -1.5 1.5 y 0 15 施加局部调整力;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;出图;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
set plot bmp
plot set mag 2.44 rotation 0 0 0
plot set back white
plo con sz
plo add sel geom pile
plot contour ssi ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;出图;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
保存及调用
sol
save kaiwa1.sav
new
res kaiwa1.sav
初始沉降归零
ini state 0
ini xdis 0 ydis 0 zdis 0
ini xve 0 yve 0 zve 0
自重应力场fish函数:
gen zone reflect normal -1 0 0 origin 0 0 0 左右对称镜像;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;第一层;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
_grad = -1910*(-9.81)
_szzl = -_grad*(24-18)
_szz0 = -_grad*24
Command
ini szz add _szz0 grad 0 0 _grad range z 18 24
ini szz add _szzl range z 0 18
end_command ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;最后层;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
_grad = -1900*(-9.81)
_szz0 = -_grad*18
Command
ini szz add _szz0 grad 0 0 _grad range z 0 18
end_command
End
锚杆设置:;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;(1)整体设置锚杆;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel cable id=1 begin 38.7771 0.5 24 end 25.7371 0.5 20.5059 nseg 54 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;(1)预紧力自由段;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel cable id=1 pre 15e4 range x 33.9474 38.7771 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;(1)自由段参数;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel cable id=1 prop emod 2.05e11 gr_coh 0 gr_fric 0.0 gr_k 0 ycomp 1e5 ytens 15.34e7 xcarea 5.54e-3 range x 33.9474 38.7771
;sel cable prop&
;emod 2.05e11&弹性模量
;gr_coh 1e20& 单位长度上水泥浆粘结力
;gr_fric 0.0&水泥浆摩擦角
;gr_k 1e10& 单位长度上水泥浆刚度
;ycomp 1.0e5 & 抗压强度
;ytens 15.34e5& 抗拉强度
;xcarea 5.54e-3& 横截面积
;range x 33.25 40 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;(1)锚固段参数;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel cable id=1 prop emod 2.05e11 gr_coh 1.5e10 gr_fric 25.0 gr_per 0.264 gr_k 10e9 ycomp 1.0e5 ytens 30e5 xcarea 5.54e-3 range x 25.7371 33.9474 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;记录锚杆2应力;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
his sel cables force cid=177
输出记录文件
def weiyi
loop n (8,62)
s=8478
e=8488
command
his write n vs step begin s end e file 底板位移记录.txt endcommand
endloop
end
weiyi。

相关文档
最新文档