来自:台少 > 馆藏分类
配色: 字号:
01
2012-07-14 | 阅:  转:  |  分享 
  
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)=rrcos(alfa)

ytable(n+10,k)=rrsin(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指示:

Genx1,y1x2,y2,x3,y3x4,y4i=i0,i1j=j0,j1详细指令参见使用手册,FLAC程式可自动产生x,y坐标与网格的关系,但由于产生的网格坐标不易控制,将对其它的工作产生负面影响,故依使用前步骤4所建立的关系,将网格依其疏密程度需要的不同,实际控制网格的坐标。

(2)设定材料性质:prop

(3)设定外力:SetGrav,ApplyPressure,inisxx,Syy

(4)设定边界条件:fix,free

(5)求起始的应力平衡:solve

(6)储存:Save

5求工程的影响________

求出区域内的应力分布情况后,再依工程的流程及步骤阶段执行各工程进行过程的影响,建议使用以下的步骤:

(1)调出起初的应力平衡:re_____.sav

(2)设定新的材料性质:model,prop

(3)设定新的支撑性质:struct

(4)设定新的外力

(5)设定边界条件

(6)求工程时的应力平衡

(7)储存
献花(0)
+1
(本文系台少首藏)