1.FLAC3D是有限元程序吗?
答:不是!是有限差分法。
2.最先需要掌握的命令有哪些?
答:需要掌握gen,ini,app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令。
3.怎样看模型的样子?
答:ploblogro可以看到不同的group的颜色分布
4.怎样看模型的边界情况?
答:plogpfixredsk,
5.怎样看模型的体力分布?
答:plofapredsk
6.怎样看模型的云图?
答:位移:plocondis(xdis,ydis,zdis)
应力:ploconsz(sy,sx,sxy,syz,sxz)
7.怎样看模型的矢量图?
答:plodis(xdis,ydis,zdis)
8.怎样看模型有多少单元、节点?
答:ploinfo
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/Restor
17.如何暂停计算?
答:Esc
18.如何在程序中进行暂停,并可恢复计算?
答:在命令中加入pause命令,用continue进行继续
19.如何跳过某个计算步?
答:在计算中按空格键跳过本次计算,自动进入下一步
20.Fish是什么东西?
答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能
21.Fish是否一定要学?
答:可以不用,需要的时候查手册获得需要的变量就可以了
22.FLAC3D允许的命令文件格式有哪些?
答:无所谓,只要是文本文件,什么后缀都可以
23.如何调用一些可选模块?
答:configdyn(fluid,creep,cppudm)
24.如何在圆柱体四周如何施加约束条件?
可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其中r2<r1,其实就是选择range的时候选两个圆柱的差,即得到边界。
25.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?
回答:用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平面转向所要确定的平面的角.
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.绘制相关曲线命令的一些解释
1)在plothistmvsn的形式里,m代表y轴,n代表x轴(不管m,n的正负);
2)"-"表示对其值作"mirror",如果大家使用AUTOCAD的话就明白"mirror"的含义了,比如上面的hist1gpydisp的值是0~100,那么vs-1就变成-100~0.以此类推.
3)为什么要这么做呢?主要是为了符合岩土工程的习惯需要,例如桩载荷试验曲线,沉降曲线等都不是画在第一象限内.
30.如何查看剖面上任一点的位移值?
printgpdisprangex()y()z()或
printgpdisprangeid
plotgpdisprangex()y()z()或
plotgpdisprangeidA1
31.fixxyzrangex-0.10.1..........................a
fixxrangex-0.10.1................................b
请问a式与b式有何区别?
a表示在x=0这个平面是固定的
b表示在x=0这个平面只x方向不能发生运动
32.initial与apply有何区别?
initial初始化命令,如初始化计算体的应力状态等;apply边界条件限制命令,如施加边界的力、位移等约束等。
initial的应力状态会随计算过程的发生而发生改变,一般体力需要初始化,而apply施加的边界条件不会发生变化。
33.如何查看各个时段不平衡力的具体数值?
采用his来记录计算,包括位移应力等命令hisunbal
hisgp(zone)zdisrange(000)或者id=?
导出数据命令hiswritenvsmbegin时步end时步filefilename.his
n表示纪录的idm表示时步
要导出不平衡力的具体数值
hisunbal
step100000orsolve
hiswrite1vsstepbegin1end1000file123.his
使用上述命令就可以查看各个时步下的不平衡力的具体数值
34.如何显示变形轮廓线的命令?
ploskemagf10其中10为放大系数
plotaddcontdispouton时时显示云图
plotblockgroupgpnumon显示节点编号
plotinterredske看接触面的效果;
delrangegroup1not不删除group信息
delrangegroup2删除group信息
inizadd数值大小rangegroupname顺z方向移动
35dip和dd的问题
dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影是一条射线
dd表示y轴正向顺时针到那条射线的夹角。首先应该按照dd的方向大体确定这个面的朝向,
dd指的是从y轴正方向按顺时针(clockwise)方向转向所要确定面的法线方向在xy平面上的投影的夹角,
然后再确定dip,dip指的是从xy平面转向所要确定的平面的角.
如genzonereflectdip90dd90是以YZ为对称面,genzonereflectdip90dd0是以xz为对称面
建立接触面interface1facerangex(范围)y(范围)z(范围)
或者interface1facerangeplaneorig(xyz坐标)norm(法线方向)
36设置图片
watertableclear清除水位线
plsetbackwhite;设置背景
setplotbitmapsize(800,510);设置图片大小
plotsetcaptionsize37;设置图例数字大小(default=35范围10-50)
plotsetcenter100040
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判断
alfa=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个值等于变量
case2当k=2的时候执行下面语句
xtable(n+10,k)=rrtablen+10的x列第k个值等于变量
ytable(n+10,k)=0.0tablen+10的y列第k个值等于变量
end_case结束选择结构
end_loop结束第二个循环
end_loop结束第一个循环
end函数定义结束
aux1执行定义的函数aux1
fixxyzrangex-0.10.1表示在x=0这个平面是固定的fixedboundary
fixxrangex-0.10.1表示在x=0这个平面只x方向不能发生运动rollerboundary
38判断模型中网格点是一个还是两个的方法
显示模型网格点的ID,如果有两个ID号,说明是两个。节理面上的网格点应该有两个ID号。
39改变模型应力边界时FLAC与3DEC的不同之处
在3DEC中,当改变应力边界时,模型表面的应力直接被更换,FLAC中改变应力边界时,模型表面的应力为两者应力相加。
40对DD,DIP的理解
DD表示节理面正方向在XY面上的投影与Y轴的夹角,也等于节理面与XY面的交线与X轴的夹角;
DIP表示节理面与XY面的夹角。+V0K*r!F)r5_5E'v2}&j
41塑性状态的解释
-p表示某一区域的应力在模型运行过程中进入过屈服状态,但现在已经推出了屈服状态;
-n表示某一区域正处于屈服状态
42对网格点速度的理解
flac中网格点的速度表示的是单位时步网格点发生的位移,查看网格点的速度有两种方法:
(1)绘制整个区域的速度图
(2)记录某一网格点的速度
43flac中影响运行时间的因素
(1)单元材料,结构面之间的刚度对比很大;
(2)划分单元尺寸对比很大,尽量避免长细比大于5:1的细长单元。
43阻尼的选择
(1)对于静力分析,局部阻尼比较有效,是默认状态的(SETMECHANICDANMLOCAL)
(2)组合阻尼,比局部阻尼更容易在大量网格点发生刚体运动时更好的收敛到稳定状态,使用SETMECH
DAMPCOMBINE来调用组合阻尼,组合阻尼对于减小动能方面不如局部阻尼有效。
44FLAC中有多少中材料模型
内置11中本构模型
(1)零模型用于开挖回填的区域
(2)各向同性弹性模型
(3)正交各向异性弹性模型
(4)横观各向同性弹性模型
(5)德鲁克-普拉格塑性模型
(6)莫尔-库仑塑性模型一般岩土体的力学行为
(7)节理化塑性模型用于研究材料破坏后行为(如边破稳定,屈服煤矿开采)
(8)应变硬化/软化模型-库仑塑性模型薄板层状材料破坏后力学行为,
(9)双线性应变硬化/软化节理化塑性模型临近层状土层开挖
(10)双屈服塑性模型
(11)修正的剑桥粘土模型
45初始化应力时为什么用了ini命令后还用apply命令
ini命令给模型所有区域应力赋予给定的值,但是它不能保证应力是平衡的。这导致了至少存在两个潜在的问题,
第一应力可能违反该区域设定的非线性模型屈服准则,step命令运行后,模型会出现塑性流动
第二,网格边界的指定应力可能和假定的初始应力不相等,当step命令执行后网格点可能移动。
为保证应力确实达到初始状态,必须使用ini和apply命令。
46对不平衡力的理解
flac中用最大不平衡力和典型内力大的比值R表示模型的不平衡力的相对大小,该值为百分数,且从不会减小到零,因此,模型不会达到绝对稳定的平横状态,可以根据要求的精度设定R值,模型默认的R值是1e-5。
注意R值小仅表示所有网格点的力都平衡,但模型还可能发生塑性流动,通过plotvel命令查看
网格点速度,一个低振幅随即速度场是没有塑性流动的平衡状态的可靠的指示器。
在固定面上施加位移时,fix与apply不能同时使用,应该用ini和fix,才能施加上。
47.菜单plotitems/add/structureelements/geometry/结构单元的形式,把对话框里的magfac的值改成不同的数,就可以看到位移放大后的结构单元。
48.你可以用键盘上的上、下、左、右键在窗口中移动模型。
49.我理解的dip和dd(所有的面都能很好地解释
首先你要确定你定义面的法线方向,dip是倾角的意思(xoy面和定义面的夹角),法线向上dip是锐角,法线向下dip时钝角。dd是y轴顺时针转到法线在xoy面上的投影线时,所转过的角度。大家可以试一试。
50.如何能把一个PLOT的图像数据导出来以便用其他软件绘图?回答:setlogon命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了请教这个方法的具体应用步骤:我想导出一个点位移的历史曲线,就是不知道怎么样才能把数据导到那个LOG中。