FLAC3D扫盲50问答Word下载.docx

FLAC3D扫盲50问答Word下载.docx

《FLAC3D扫盲50问答Word下载.docx》由会员分享，可在线阅读，更多相关《FLAC3D扫盲50问答Word下载.docx（8页珍藏版）》请在冰豆网上搜索。

i#X0z4{）y8U）I$m（12.如何进行边界约束？

6v'

@&

{0M2S9d5n%A2~6w9W答：

fixxran（约束的是速度，在初始情况下约束等效于位移约束）

A（x;

x5U3g 

q13.如何知道每个单元的ID？

-T:

T+~1s&

B 

y8g9M（z%[*c.G答：

用鼠标双击单元的表面，可以知道单元的ID和坐标

14.如何进行切片？

3Z1y0F%l/W（@#C,x6a6_答：

plosetplaneori（点坐标）norm（法向矢量）

!

M1J&

z1e+H+p&

_7`）m 

ploconszplane（显示z方向应力的切片）

%H7o2K'

E1Z#_.u+l2_*J%` 

15.如何保存计算结果？

save+文件名

16.如何调用已保存的结果？

0h!

o/R!

f'

Z答：

rest+文件名；

或者File/Restor

9P'

A;

\0j$r8I-K!

H1H（W17.如何暂停计算？

4S（^;

h&

B%D4t'

j答：

Esc

2w）W*h8G7H.?

$_18.如何在程序中进行暂停，并可恢复计算？

在命令中加入pause命令，用continue进行继续

-{#B（p7d7Q3H/g%V-D"

L#?

a$~/Z"

_9\,w:

L:

z!

Q19.如何跳过某个计算步？

在计算中按空格键跳过本次计算，自动进入下一步

20.Fish是什么东西？

N5I+]6B&

o.x答：

是FLAC3D的内置语言，可以用来进行参数化模型、完成命令本身不能进行的功能

21.Fish是否一定要学？

可以不用，需要的时候查手册获得需要的变量就可以了

4j-B'

h4?

J*B8n22.FLAC3D允许的命令文件格式有哪些？

4t 

h8r9u&

N-r!

e9z9S#E7x+p答：

无所谓，只要是文本文件，什么后缀都可以

23.如何调用一些可选模块？

1}&

R*l4f8d4t3b答：

configdyn（fluid,creep,cppudm）

24.如何在圆柱体四周如何施加约束条件？

/e,D/?

6B7M6h%O:

J可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not，其中r2<

r1，其实就是选择range的时候选两个圆柱的差，即得到边界。

;

M,[,m4N&

T'

c%y'

G525. 

如何能把一个PLOT的图像数据导出来以便用其他软件绘图？

'

d&

|"

T"

a8}6C7R8P回答：

用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

26.关于dd、dip的理解

dip-directionangle,measuredintheglobalxy-planeclockwisefromthepositivey-axis这个是从y轴到x顺时针方向所转过的角度；

dipdipangle,measuredinthenegativez-directionfrom　theglobalxy-plane这个是xy平面与z轴负方向的夹角。

dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线，dd表示y轴正向顺时针到那条射线的夹角。

首先应该按照dd的方向大体确定这个面的朝向，dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．

*}3x:

L&

o/g,l 

T,D9t&

]/~5F'

P8`2N3F27. 

如何调用FISH文件？

新手调用FISH文件出错的原因：

是：

文件名为fishcall.fis.txt（因为在记事本里这个txt的后缀是默认的）。

修改方法是：

将fishcall.fis.txt里面的内容全部复制到新的记事本里，再按保存，保存时候将保存类型改为所有文件（*,*）,这时你便可在文件名栏内输入fishcall.fis，类型也就正确了。

3M3M 

b;

}3v4h+m&

u!

t（T-P/A28. 

如何理解Tension-pshear-ptension-nshear-n表示的含义？

A6o-C2M"

S7^#G6^6^*@）s在计算的循环里面,每个循环中,每个zone都依据failurecriterion处于不同的状态,shear和tension分别表示因受剪和受拉而处于塑性状态.n表示now,p表示previous,即分别表示在现在和以前的循环当中处于塑性状态.

l8d4B6i!

j,B6h!

i8]:

@!

d.v29. 

绘制相关曲线命令的一些解释

c3K#u;

Y$X4O1）在plothistmvsn的形式里,m代表y轴,n代表x轴（不管m,n的正负）;

2）"

-"

表示对其值作"

mirror"

如果大家使用AUTOCAD的话就明白"

的含义了,比如上面的hist1gpydisp的值是0~100,那么vs-1就变成-100~0.以此类推.

0y;

A0]"

K 

B*a8L3\3）为什么要这么做呢?

主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内.

30. 

如何查看剖面上任一点的位移值？

printgpdisprangex（）y（）z（）或

f/h;

j9J*?

#a6G6Rprintgpdisprangeid

0R#a:

J!

`"

I 

Oplotgpdisprangex（）y（）z（）或

+u 

_-B）F6@'

P!

L7Y4h3M%Qplotgpdisprangeid 

A1

.i!

?

.X4j&

i7d9m"

{*]31. 

fixxyzrangex-0.10.1..........................a

y3h5J1u.O.i6b/C/C 

fixxrangex-0.10.1................................b

/A 

O$C!

l3Z%f%N-k9f#| 

请问a式与b式有何区别？

:

n7I;

O（l7n6u&

Fa表示在x=0这个平面是固定的

5O#c8N&

[&

X'

K9Qb表示在x=0这个平面只x方向不能发生运动

B%]1b 

c.f!

32. 

initial与apply有何区别？

.r（J"

q.]1y（b,w4T#Kinitial初始化命令，如初始化计算体的应力状态等；

apply边界条件限制命令，如施加边界的力、位移等约束等。

%M0f0G'

{+m1w5e）ainitial的应力状态会随计算过程的发生而发生改变，一般体力需要初始化，而apply施加的边界条件不会发生变化。

e6j5@4S 

p3d2w（@"

o9}33. 

如何查看各个时段不平衡力的具体数值？

6H+K#X9P6^&

B:

B3Z-B采用his来记录计算，包括位移应力等命令hisunbal

hisgp（zone）zdisrange（000）或者id=？

导出数据命令hiswritenvsmbegin时步 

end 

时步filefilename.his

%]:

s$c7`7g（q;

In表示纪录的idm表示时步

7e4G;

E6Z'

d 

q#H7g要导出不平衡力的具体数值

hisunbal

*?

K（A+s$M6O0]step100000orsolve

*I*T&

j*Z7G:

x*e）d（s0ghiswrite1vsstepbegin1end1000file123.his

1p+a!

W9o*F;

s/h6w使用上述命令就可以查看各个时步下的不平衡力的具体数值

%i&

`*V,b）s'

{.A/Q34.如何显示变形轮廓线的命令？

Y-M）S9U%k,Xploskemagf10 

其中10为放大系数 

plotaddcontdispouton 

时时显示云图

plotblockgroupgpnumon 

显示节点编号

$V&

.A8x1q/]8?

.v7y-N5];

j8S"

d6Nplotinterredske 

看接触面的效果;

delrangegroup1not不删除group信息 

delrangegroup2删除group信息

9^:

M!

\4f8?

#O$f"

finizadd数值大小rangegroupname顺z方向移动

35dip和dd的问题

dip表示对称参照面与xy平面的夹角，对称参照面与xy平面的夹角在xy平面的投影是一条射线

dd表示y轴正向顺时针到那条射线的夹角。

首先应该按照dd的方向大体确定这个面的朝向，

#R,\/S9[*l1j1W'

~dd指的是从y轴正方向按顺时针（clockwise）方向转向所要确定面的法线方向在xy平面上的投影的夹角，

然后再确定dip，dip指的是从xy平面转向所要确定的平面的角．

如genzonereflectdip90dd90是以YZ为对称面，genzonereflectdip90dd0是以xz为对称面

.n）t;

\5Y 

k"

E4_7j（K建立接触面interface1facerangex（范围）y（范围）z（范围）

#^0}#Z+b4|/P1W%l 

或者interface1facerangeplaneorig（xyz坐标）norm（法线方向）

36设置图片

watertableclear清除水位线

（u'

Z%b8h）G*Gplsetbackwhite;

设置背景

.]2S&

f.a"

e8k#r$wsetplotbitmapsize（800,510）;

设置图片大小

plotsetcaptionsize37;

设置图例数字大小（default=35范围10-50）

2V8{3M1Z9~1M+Dplotsetcenter100040

37fish命令的解释（网上摘录）

defaux1定义函数aux1

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判断

#S）g'

I9l/E0l5Jalfa=0.5*（pi+.1）*float（k-2）/float（nprof-2）

0\9m3J'

c,^*?

xtable（n+10,k）=rr*cos（alfa）

%D!

Q%D1k4Q4Y3m（G9[ 

w#fytable（n+10,k）=rr*sin（alfa）

case1当k＝1的时候执行下面语句

7[（z6n6I9U7[xtable（n+10,k）=rrtablen+10的x列第k个值等于变量rr

4?

）K7D:

@#K7Nytable（n+10,k）=-40.0tablen+10的y列第k个值等于变量

case2当k＝2的时候执行下面语句

;

I0A;

E5R（N1b1o$A.Q!

^xtable（n+10,k）=rrtablen+10的x列第k个值等于变量

ytable（n+10,k）=0.0tablen+10的y列第k个值等于变量

end_case结束选择结构

（B!

I2m,I/C:

_.o 

l6Tend_loop 

结束第二个循环

\${'

o*Z3`/Q-C%f5q2jend_loop结束第一个循环

）b1Q:

X*U）S4G9y 

`8E0A8hend函数定义结束

aux1执行定义的函数aux1

fixxyzrangex-0.10.1表示在x=0这个平面是固定的fixedboundary

fixxrangex-0.10.1表示在x=0这个平面只x方向不能发生运动rollerboundary

38判断模型中网格点是一个还是两个的方法

3W$f"

["

} 

G-h 

显示模型网格点的ID，如果有两个ID号，说明是两个。

节理面上的网格点应该有两个ID号。

39改变模型应力边界时FLAC与3DEC的不同之处

%Q,T 

~6D'

q6^ 

在3DEC中，当改变应力边界时，模型表面的应力直接被更换，FLAC中改变应力边界时，模型表面的应力为两者应力相加。

40,F8X'

R;

|7W#Q4040对DD,DIP的理解

DD表示节理面正方向在XY面上的投影与Y轴的夹角，也等于节理面与XY面的交线与X轴的夹角；

（o）S5v6~9I"

X9o1[&

m$?

$t 

DIP表示节理面与XY面的夹角。

+V0K*r!

F）r5_5E'

v2}&

j

41塑性状态的解释

7W2u 

a 

n!

J.U$f;

[）K 

-p 

表示某一区域的应力在模型运行过程中进入过屈服状态，但现在已经推出了屈服状态；

－n 

表示某一区域正处于屈服状态 

42对网格点速度的理解

&

I#r;

B6Y*t"

u1F'

M5`3} 

flac中网格点的速度表示的是单位时步网格点发生的位移，查看网格点的速度有两种方法：

#S!

O6]2y$O）{,m&

p

（1）绘制整个区域的速度图

l.`"

F5c$E.J4Q

（2）记录某一网格点的速度

4k6W2s#|8z434b43flac中影响运行时间的因素

（1）单元材料，结构面之间的刚度对比很大；

4{&

N）n-Q*K

（2）划分单元尺寸对比很大，尽量避免长细比大于5：

1的细长单元。

43阻尼的选择

+u*Y（N9n1f7U 

（1）对于静力分析，局部阻尼比较有效，是默认状态的（SETMECHANICDANMLOCAL）

（2）组合阻尼，比局部阻尼更容易在大量网格点发生刚体运动时更好的收敛到稳定状态，使用SETMECH

DAMPCOMBINE来调用组合阻尼，组合阻尼对于减小动能方面不如局部阻尼有效。

5c,Q3M1Y-^ 

t/e 

L4444FLAC中有多少中材料模型

内置11中本构模型

6O2[2C:

],y9N 

（1）零模型 

用于开挖回填的区域

（2）各向同性弹性模型

（3）正交各向异性弹性模型

（4）横观各向同性弹性模型

（5）德鲁克－普拉格塑性模型

）?

5{ 

y'

M*~'

l*q 

c 

（6）莫尔－库仑塑性模型 

一般岩土体的力学行为

.A,c#r*T）V3d!

J1x（H 

（7）节理化塑性模型 

用于研究材料破坏后行为（如边破稳定，屈服煤矿开采）

（8）应变硬化/软化模型－库仑塑性模型 

薄板层状材料破坏后力学行为,

（9）双线性应变硬化/软化节理化塑性模型 

临近层状土层开挖

（10）双屈服塑性模型

（11）修正的剑桥粘土模型

{7],y9r.^2s;

[4545初始化应力时为什么用了ini命令后还用apply命令

ini命令给模型所有区域应力赋予给定的值，但是它不能保证应力是平衡的。

这导致了至少存在两个潜在的问题，

第一应力可能违反该区域设定的非线性模型屈服准则，step命令运行后，模型会出现塑性流动

#I8B,S/d'

e;

T第二，网格边界的指定应力可能和假定的初始应力不相等，当step命令执行后网格点可能移动。

#I*M"

E;

G1k0U4d.j/H为保证应力确实达到初始状态，必须使用ini和apply命令。

3o 

r,C'

i4u 

~#\;

l"

k;

B.}）T+Y3V&

y$Q&

H4646对不平衡力的理解

flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小，该值为百分数，且从不会减小到零，因此，模型不会达到绝对稳定的平横状态，可以根据要求的精度设定R值，模型默认的R值是1e-5。

y&

G2Q5W"

z&

o6[/e 

注意R值小仅表示所有网格点的力都平衡，但模型还可能发生塑性流动，通过plotvel命令查看

网格点速度，一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。

在固定面上施加位移时，fix与apply不能同时使用，应该用ini和fix,才能施加上。

47.菜单plotitems/add/structureelements/geometry/结构单元的形式，把对话框里的magfac的值改成不同的数，就可以看到位移放大后的结构单元。

48.你可以用键盘上的上、下、左、右键在窗口中移动模型。

49.我理解的dip和dd（所有的面都能很好地解释

首先你要确定你定义面的法线方向，dip是倾角的意思（xoy面和定义面的夹角），法线向上dip是锐角，法线向下dip时钝角。

dd是y轴顺时针转到法线在xoy面上的投影线时，所转过的角度。

大家可以试一试。

50.

25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图？

26.

27.

28.

29.

30.回答：

31.

32.

33.setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了

请教这个方法的具体应用步骤：

我想导出一个点位移的历史曲线，就是不知道怎么样才能把数据导到那个LOG中。

1、形成切面的操作;

以0为中心取切面plotsetplaneori000norm010$U;

输出szz方向的云图3o7cplotconszzplane;

添加网格线jplotadddisplane-|8t_#X9I;

在切面上添加位移矢量2、生产变形矢量图plotskdis;

输出变形网格plotskmagf20;

（变形倍数）plotconszzouonmagef20;

变形后的网格与应力云图结合起来3、显示单元的id的网格模型+r5?

/|3plotblockgroupidon:

x2C6x54、指定单元的应力，两个办法：

/先plotblockmodelidon，然后找到你想看的那个单元的id，假设是10，然后就plotconszzoutonrangeid10；

再者，用fish，通过gp_head，gp_near，gp_id等等获得你要的节点，然后可以得到gridpoint的各种信息。

5、要输出应力则应监测单元，要输出位移则应监测节点！

对于如何输出，是比较简单的，在计算前你要根据自己的需要，设置监测单元及监测点，这样计算完成后才能输出相应的应力或位移值。

设置监测单元应力historyid=2zonesxxid=776historyid=3zonesxxid=794historyid=4zonesxxid=1889q0H9u/j+o*W"

qhistoryid=6zonesxxid=1835输出：

setlogon+Z）{1L-q+Esetlogfileoutp_str.txt-p'

`7u&

R+~*m/`1orestoriginal_str.savprzonestressranid776anyid794anyid1889anyid1835anyprzoneprincipalranid776anyid794anyid1889anyid1835any7`6setlogoff6、在运