Flac学习入门.docx
- 文档编号:10540454
- 上传时间:2023-05-26
- 格式:DOCX
- 页数:16
- 大小:23.27KB
Flac学习入门.docx
《Flac学习入门.docx》由会员分享,可在线阅读,更多相关《Flac学习入门.docx(16页珍藏版)》请在冰点文库上搜索。
Flac学习入门
FLAC学习
1.FLAC3D是有限元程序吗?
答:
不是!
是有限差分法。
2.最先需要掌握的命令有哪些?
答:
需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。
3.怎样看模型的样子?
答:
ploblogro可以看到不同的group的颜色分布
4.怎样看模型的边界情况?
答:
plogpfixred
5.怎样看模型的体力分布?
答:
plofapred
6.怎样看模型的云图?
答:
位移:
plocondis(xdis,ydis,zdis)
应力:
ploconsz(sy,sx,sxy,syz,sxz)
7.怎样看模型的矢量图?
答:
plodis(xdis,ydis,zdis)
8.怎样看模型有多少单元、节点?
答:
priinfo
9.怎样输出模型的后处理图?
答:
File/Printtype/Jpgfile,然后选择File/Print,将保存格式选择为jpe文件
10.怎样调用一个文件?
答:
File/call或者call命令
10.如何施加面力?
答:
appnstress
11.如何调整视图的大小、角度?
答:
综合使用x,y,z,m,Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键
12.如何进行边界约束?
答:
fixxran(约束的是速度,在初始情况下约束等效于位移约束)
13.如何知道每个单元的ID?
答:
用鼠标双击单元的表面,可以知道单元的ID和坐标
14.如何进行切片?
答:
plosetplaneori(点坐标)norm(法向矢量)
ploconszplane(显示z方向应力的切片)
15.如何保存计算结果?
答:
save+文件名.
16.如何调用已保存的结果?
答:
rest+文件名;或者File/Restore
17.如何暂停计算?
答:
Esc
18.如何在程序中进行暂停,并可恢复计算?
答:
在命令中加入pause命令,用continue进行继续
19.如何跳过某个计算步?
答:
在计算中按空格键跳过本次计算,自动进入下一步
20.Fish是什么东西?
答:
是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能
21.Fish是否一定要学?
答:
可以不用,需要的时候查Mannual获得需要的变量就可以了
22.FLAC3D允许的命令文件格式有哪些?
答:
无所谓,只要是文本文件,什么后缀都可以
23.如何调用一些可选模块?
答:
configdyn(fluid,creep,cppudm)
24.如何在圆柱体四周如何施加约束条件?
答:
可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其
中r2 25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图? 答: 用setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘 制就行了。 26.关于dd、dip的理解: 答: 首先你要确定你定义面的法线方向,dip是倾角的意思(xoy面和定义面的夹角),法线向上dip是锐角,法线向下dip时钝角。 dd是y轴顺时针转到法线在xoy面上的投影线时,所转过的角度。 27.如何调用FISH文件? 答: 新手调用FISH文件出错的原因: 是: 文件名为fishcall.fis.txt(因为在记事本里这个txt的后 缀是默认的)。 修改方法是: 将fishcall.fis.txt里面的内容全部复制到新的记事本里,再按保存 ,保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入fishcall.fis,类型也就 正确了。 28.如何理解Tension-pshear-ptension-nshear-n表示的含义? 答: 在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和 tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前 的循环当中处于塑性状态. 29.绘制相关曲线命令的一些解释 在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了 30.如何查看剖面上任一点的位移值? 答: printgpdisprangex()y()z()或 printgpdisprangeid plotgpdisprangex()y()z()或 plotgpdisprangeid 31.a: fixxyzrangex-0.10.1 b: fixxrangex-0.10.1请问a式与b式有何区别? 答: a表示在x=0这个平面是固定的;b表示在x=0这个平面只x方向不能发生运动 32.initial与apply有何区别? 答: 初始化命令,如初始化计算体的应力状态等; apply边界条件限制命令,如施加边界的力、位移等约束等。 initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件 不会发生变化。 33.如何查看各个时段不平衡力的具体数值? 答: 采用his来记录计算,包括位移应力等命令hisunbal hisgp(zone)zdisrange(000)或者id=? 导出数据命令hiswritenvsmbegin时步end时步+filefilename.his 要导出不平衡力的具体数值 hisunbal step100000orsolve hiswrite1vsstepbegin1end1000file123.his 使用上述命令就可以查看各个时步下的不平衡力的具体数值 34.如何显示变形轮廓线的命令? ploskemagf10其中10为放大系数 35、怎样开始入门? 答: 初步学习首先看说明书User'sgurid的Section2Gettingstarted和Section3Problem slovingwithFLAC 36、哪些命令需要掌握? 答: gridgenpropinterfaceinifixsethiststepsolvepausecontsav等,在学习的过程 中遇到这些命令和不会的命令应该查看说明书的CommandReference,里面讲的很清楚,特别需要注 意的是一定要仔细看一下CommandReference开头的OptionalRangePhrases,这是讲范围的,因为 很多命令的施加都需要加到模型的一部分,这里就告诉你怎样来表示这一部分。 这一原则同样适合与 其他5个ITASCA软件。 还有一点就是六个软件(FLACFALC3DUDEC3DECPFCPFC3D)的range大不相 37、FLAC和FLAC3D建模有什么区别? 答: FLAC和FLAC3D基本的建模思想大致是一样的,不同的是FLAC是首先用GRID命令把整体的网格画出 来,然后用MODELLNULL命令把网格分成许多部分,再用GEN命令把这些部分一块一块搭接起来。 然后 再用GENLINE,GENARC等命令画出要开挖的边界,在FLAC中产生网格的命令只有一个GRID,而且只 能用一次。 而FLAC3D是直接用GEN命令一块一块搭接,GEN可以用很多次。 38、怎样查看模型? 答: plotgrid可以查看网格,plotgridnum可以查看节点号。 39、怎样查看边界条件? 答: plotfix 40、怎样查看预先留设的开挖边界? 答: plotgridmark。 41、怎么查看模型的云图? 答: plsyy(sxxxdisydissig1..)fill,注意,带fill关键字就是查看云图,不带fill关键字 就是查看等直线图。 42、怎样查看模型的矢量图? 答: plgridvel(disp)red 43、怎样查看INTERFACE? 答: plgridifaceidred printifaceid;id是interface的编号 44、FLAC怎么调整视图? 答: ctrl+z或shift+z,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小, 双击屏幕。 45、判断模型中网格点是一个还是两个的方法 显示模型网格点的ID,如果有两个ID号,说明是两个。 节理面上的网格点应该有两个ID号。 46、改变模型应力边界时FLAC与3DEC的不同之处 在3DEC中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界时,模型 表面的应力为两者应力相加。 47、flac中影响运行时间的因素 (1)单元材料,结构面之间的刚度对比很大; (2)划分单元尺寸对比很大,尽量避免长细比大于5: 1的细长单元。 48FLAC中有多少中材料模型 内置11中本构模型%_4[W/N,nR5q (1)零模型用于开挖回填的区域 (2)各向同性弹性模型 (3)正交各向异性弹性模型 (4)横观各向同性弹性模型 (5)德鲁克-普拉格塑性模型 (6)莫尔-库仑塑性模型一般岩土体的力学行为 (7)节理化塑性模型用于研究材料破坏后行为(如边破稳定,屈服煤矿开采) (8)应变硬化/软化模型-库仑塑性模型薄板层状材料破坏后力学行为 (9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖 (10)双屈服塑性模型 (11)修正的剑桥粘土模型 49、对不平衡力的理解 flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小,该值为百分数,且从 不会减小到零,因此,模型不会达到绝对稳定的平横状态,可以根据要求的精度设定R值,模型默认 的R值是1e-5。 注意R值小仅表示所有网格点的力都平衡,但模型还可能发生塑性流动,通过plotvel命令查看 网格点速度,一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。 50、简单命令解释: plotaddcontdispouton时时显示云图 plotblockgroupgpnumon显示节点编号 plotinterredske看接触面的效果. delrangegroup1not不删除group信息 delrangegroup2删除group信息 inizadd数值大小rangegroupname顺z方向移动 建立接触面interface1facerangex(范围)y(范围)z(范围) 或者interface1facerangeplaneorig(xyz坐标)norm(法线方向) watertableclear清除水位线 plsetbackwhite;设置背景 setplotbitmapsize(800,510);设置图片大小 plotsetcaptionsize37;设置图例数字大小(default=35范围10-50) plotsetcenter100040 fixxyzrangex-0.10.1表示在x=0这个平面是固定的fixedboundary fixxrangex-0.10.1表示在x=0这个平面只x方向不能发生运动rollerboundary 51、fish命令的解释: defaux1定义函数 loopn(1,nptab1)开始一个循环,从1循环到变量nptab1 rr=xtable(1,n)从table1的x列中读第n个数付给变量rr zz=ytable(1,n)从table1的y列中读第n个数付给变量zz loopk(1,nprof)开始第二个循环,从1循环到变量nprof case_ofk下面是一个选择结构,利用变量k判断 iFalfa=0.5*(pi+.1)*float(k-2)/float(nprof-2) xtable(n+10,k)=rr*cos(alfa) ytable(n+10,k)=rr*sin(alfa) case1当k=1的时候执行下面语句 xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rr ytable(n+10,k)=-40.0tablen+10的y列第k个值等于变量-40 case2当k=2的时候执行下面语句 xtable(n+10,k)=rrtablen+10的x列第k个值等于变量rr ytable(n+10,k)=0.0tablen+10的y列第k个值等于变量0 end_case结束选择结构 end_loop结束第二个循环 end_loop结束第一个循环 end函数定义结束 aux1执行什么定义的函数 52.在固定面上施加位移时,fix与apply不能同时使用,应该用ini和fix,才能施加上。 53.菜单plotitems/add/structureelements/geometry/结构单元的形式,把对话框里的magfac的值改成不同的数,就可以看到位移放大后的结构单元。 54.你可以用键盘上的上、下、左、右肩在窗口中移动模型。 55.如何能把一个PLOT的图像数据导出来以便用其他软件绘图? 答: 用setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。 56在进行三轴压缩时候,如何施加围压? 答: 使用applysxx=-1e5szz=-1e5rangecylend10,0,0end20,2,0radius1,其中X和Z轴是在我们眼睛正对的面。 57在进行三轴压缩时候,如何对轴向力进行位移控制? 答: 使用iniyvel1e-7rangey-.1.1iniyvel1e-7rangey1.92.1即可,数值可以变化 58在进行边坡计算的时候,为什么用sol,fos命令运行不行? 答: 应该输入solfosslope3dfos.savassociated(不要也行)再用restore命令调用即可。 59在FLAC3D计算中,如何记录程序中的值,如果不想时时记录该怎么弄? 答: 用histwtiten1vsn2file.txt,想间隔记录可用sethist_repn。 60想记录接触面上的相关力学参数该怎么办? 答: histinterinstress/ndis/sstress/sdisp 61接触面设置参数包括哪几个? 答: cohdialfricknksten. 62如何输出单元格的数据? 答;expgrid.flacsd 63如何在剖面的沉降云图上增加网格、剖面上的矢量及坐标值? DQZ9s1}e4YU%q 答: ploaddske ploadddisplane ploaddaxe 64如何做切片? 答: plotsetplaneorixyznormxyz 65举一个隧道支护开挖的例子? genzoneradbrickp0(0,0,0)p1(10,0,0)p2(0,10,0)p3(0,0,10)& size3557rat1111.5dim142fill modelmohr propbulk1e8shear.3e8fric35 propcoh1e3tens1e3 fixxrangex-0.10.1 fixzrangez-0.10.1 fixyrangey9.910.1 setgrav0100 inidens1000 inisyy=0.0grad0,-10000,0 inisxx=0.0grad0,-5000,0 iniszz=0.0grad0,-5000,0 applysxx=0.0grad0,-5000,0rangex9.9,10.1 applyszz=0.0grad0,-5000,0rangez9.9,10.1 setlarge histunbal solve savetr_eq.sav ;monitortrenchexcavation inixdis=0ydis=0zdis=0 histgpxdisp1,0,0 histgpzdisp0,0,2; ;excavationsteps rangenameexc1x01y0.00.8z02 rangenameexc2x01y0.81.6z02 rangenameexc3x01y1.62.4z02 rangenameexc4x01y2.43.2z02 rangenameexc5x01y3.24.0z02 ;excavationstep1 modelnullrangeexc1 step10000 saveexc1.sav ;excavationstep2 modelnullrangeexc2 step1000 saveexc2.sav ;excavationstep3 modelnullrangeexc3 step1000 saveexc3.sav ;excavationstep4 modelnullrangeexc4 step1000 saveexc4.sav ;excavationstep5 modelnullrangeexc5 step1000 saveexc5.sav ;installcablesupport restexc2.sav selcablebeg1.0,0.4,1.5end5.0,0.4,1.5nseg4 selcablebeg1.0,0.4,0.5end5.0,0.4,0.5nseg4 selcablebeg1.0,1.2,1.5end5.0,1.2,1.5nseg4 selcablebeg1.0,1.2,0.5end5.0,1.2,0.5nseg4 selcablepropemod2e9ytension1e8xcarea1.0gr_coh1e10gr_k2e9gr_per1.0 modelnullrangeexc3 step1000 savecab3.sav selcablebeg1.0,2.0,1.5end5.0,2.0,1.5nseg4 selcablebeg1.0,2.0,0.5end5.0,2.0,0.5nseg4 selcablepropemod2e9ytension1e8xcarea1.0gr_coh1e10gr_k2e9gr_per1.0 modelnullrangeexc4; step1000 savecab4.sav selcablebeg1.0,2.8,1.5end5.0,2.8,1.5nseg4 selcablebeg1.0,2.8,0.5end5.0,2.8,0.5nseg4 selcablepropemod2e9ytension1e8xcarea1.0gr_coh1e10gr_k2e9gr_per1.0 modelnullrangeexc5 step1000 savecab5.sav plotcreateTrenchSupport plotsetrot63.89,334.35,3.4 plotsetcent2.893,3.446,3.899 plotsetsize4.136 plotaddcontdispouton plotaddselcableforce plotaddaxesblack 66举一个三轴压缩的例子? genzonecylp0000p1100p2020p3001size454 genzonereflectnorm1,0,0 genzonereflectnorm0,0,19 modelmohr propbulk1.19e10shear1.1e10 propcoh2.72e5fric44ten2e5 fixxyzrangey-.1.1 fixxyzrangey1.92.1 applysxx=-1e5szz=-1e5rangecylend10,0,0end20,2,0radius1 iniyvel0grad0-1e-40rangey02 sol 1.FLAC程式使用前准备步骤 步骤1: 依比例画出所欲分析的资料 在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。 步骤2: 换算输入资料成同一单位 将现有地层资料,如Density,Bulkmodulus,Young`smodulus,tension,cohesion,frictionAngle等资料,换算成同一单位。 附注: 需谨慎检查输入资料的单位,如因单位不同而造成过大或过小的值,将会造成FLAC无法计算,而产生ERROR信息。 步骤3: 应用公式简略计算 应用公式或依据经验,简略算出FLAC输出资料的范围,以作为Debug及输出资料分析时验证。 步骤4: 建立x,y坐标与nodei,j之间的关系 在图上距离及深度的关系,建立x,y坐标系统,再由x,y坐标系统,转换与网格间系,为了便于以后输出资料的分析,故应确实掌握网格的位置及其相对应的x,y坐标。 建议在敏感区域使用较密的网格,其它地方则使用较疏的网格,刚开始执行程序时,不宜使用网格太大的网格数目,因尽量使网格总数少于1000,以节省时间。 2.FLAC3D程序的编写步骤: 1Config________ 2Grid________ 3Model________ 4求起始的应力平衡________ (1)建立x,y坐标与网格的关系,建议使用Gen指示
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Flac 学习 入门