FLAC3D学习疑点详解Word文件下载.docx
《FLAC3D学习疑点详解Word文件下载.docx》由会员分享，可在线阅读，更多相关《FLAC3D学习疑点详解Word文件下载.docx（39页珍藏版）》请在冰豆网上搜索。
初始网格
genzoneradcylp0150-15.75p1add1500p2add060p3add0015.75&
size3485dim2.252.252.252.25rat1111.1fill
genzoneradtunsize3435dim0.752.250.752.25rat1111.1&
p0150-15.75p1add00-14.25p2add060p3add1500fill
genzonereflectdip90dd270ori150-15.75
摩尔库仑模型
modelmohr
macroniyan'
bulk3.2825e3shear1.5150e3fric20coh0.07'
macromeiceng'
bulk2.1735e3shear8.8917e2fric19coh0.05'
propniyanrangez=-13.50
propmeicengrangez=-15.5-13.5
propniyanrangez=-30-15.5
为开挖定义整体部分
rangenameinz-15.750cylend1150-15.75end2156-15.75rad2.25
rangenamehangdaox12.7517.25y06z-16.50-15.75
进行加载
fixxyzrangez-30.001-29.999
appsxx=-20rangex-0.0010.001
appsxx=-20rangex29.99930.001
appsyy=-20rangey-0.0010.001
appsyy=-20rangey5.9996.001
appszz=-20rangez-0.0010.001
solve
（1）
modelnullrangein
modelnullrangehangdao
rangenameinz-15.750cylend1150-15.75end2156-15.75rad2.25
rangenamehangdaoliangce1x12.74912.751y06z-16.50-15.75
rangenamehangdaoliangce2x17.24917.251y06z-16.50-15.75
rangenamehangdaox12.7517.25y06z-16.501-16.499
applynstress-0.05planedip90dd0rangein
给巷道底部加压力
appnstress-0.05planedip90dd0rangehangdao
给巷道直墙部分（左侧）加压力
applynstress-0.05planedip90dd0rangehangdaoliangce2
给巷道直墙部分（右侧）加压力
applynstress-0.05planedip90dd0rangehangdaoliangce1;
对点进行监测
histunbal
histgpxdisp12.750-15.75 :
监测左帮
histgpzdisp150-13.5;
监测顶板
histgpzdisp150-16.5;
监测底板
solve
问：
第一个solve平衡了,是第二个solve不能平衡，这是为什么？
如果我把支护力加大，也不能平衡
我遇到的不收敛的原因，有两个：
1）检查你的模型是不是没有建好，相互之间是不是没有连接好；
2）仔细检查你的参数，是不是已经发生了破坏。
力施加有问题，隧道底部隆起破坏，所以无法计算
我这个模型就是该模型就是要出现巷道底臌现象，如果按照你的意思，力应该怎样施加才对？
如果无法收敛，应该用step命令，关键问题是要运行多少步。
3.问：
g1010
me
defdistort
loopi（2,jzones）
loopj（2,jzones）
x（i,j）=x（i,j）+mag*（urand-0.5）
y（i,j）=y（i,j）+mag*（urand-0.5）
end_loop
end
propd1000sh1e8bu2e8
其中：
x（i,j）=x（i,j）+mag*（urand-0.5）中的urand是否是随机生成数的命令？
urand：
randomnumberdrawnfromuniformdistributionbetween 0.0and1.0；
grand：
randomnumberdrawnfromnormaldistribution:
mean=0.0;
standarddeviation=1.0
4.问：
怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？
用setlogon命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。
5.问：
用命令建立模型后，如何显示点的坐标？
使用ploblogrogpnumon命令
6.interface问题
我加了一个层面，初始平衡就出现了这么大的位移？
是怎么原因？
No1
genzonebrickp0000p11000p20100p3008&
size886rat1.021.021.02
genzonebrickp0009p11009p20109p30011&
size882rat1.021.021.02
genzonereflectdip90dd90origin（0,0,6）rangez011
grouplayer1rangez911
interface1facerangeplanenorm（0,0,1）origin（0,0,8）dist0.01
inizadd-1.0rangegrouplayer1
No2
genzonebrickp00011p110011p201011p30013&
genzonereflectdip90dd90origin（0,0,1）rangez1113
grouplayer2rangez1113
interface2facerangeplanenorm（0,0,1）origin（0,0,10）dist0.01
inizadd-1.0rangegrouplayer2
No3
genzonebrickp00013p110013p201013p30014&
size881rat1.021.021.02
genzonereflectdip90dd90origin（0,0,1）rangez1314
grouplayer3rangez1314
interface3facerangeplanenorm（0,0,1）origin（0,0,12）dist0.01
inizadd-1.0rangegrouplayer3
No4
genzonebrickp00014p110014p201014p30015&
genzonereflectdip90dd90origin（0,0,1）rangez1415
grouplayer4rangez1415
interface4facerangeplanenorm（0,0,1）origin（0,0,13）dist0.01
inizadd-1.0rangegrouplayer4
No5
genzonebrickp00015p110015p201015p30017&
genzonereflectdip90dd90origin（0,0,15）rangez1517
grouplayer5rangez1517
interface5facerangeplanenorm（0,0,1）origin（0,0,14）dist0.01
inizadd-1.0rangegrouplayer5
No6
genzonebrickp00017p110017p201017p30018&
genzonereflectdip90dd90origin（0,0,1）rangez1718
grouplayer6rangez1718
interface6facerangeplanenorm（0,0,1）origin（0,0,16）dist0.01
inizadd-1.0rangegrouplayer6
No7
genzonebrickp00018p110018p201018p30020&
genzonereflectdip90dd90origin（0,0,1）rangez1820
grouplayer7rangez1820
interface7facerangeplanenorm（0,0,1）origin（0,0,17）dist0.01
inizadd-1.0rangegrouplayer7
No8
genzonebrickp00020p110020p201020p30024&
size883rat1.021.021.02
genzonereflectdip90dd90origin（0,0,17）rangez2024
grouplayer8rangez2024
interface8facerangeplanenorm（0,0,1）origin（0,0,19）dist0.01
inizadd-1.0rangegrouplayer8
No9
genzonebrickp00024p110024p201024p30029&
size884rat1.021.021.02
genzonereflectdip90dd90origin（0,0,0）rangez2429
grouplayer9rangez2429
interface9facerangeplanenorm（0,0,1）origin（0,0,23）dist0.01
inizadd-1.0rangegrouplayer9
----------边界条件------------------------
fixxrangex-20.1-19.9
fixxrangex19.920.1
fixyrangey-0.10.1
fixyrangey9.910.1
fixzrangez-0.10.1
-----------初始条件-------------------------
applyszz-26e6rangez27.928.1
initialszz-1.08e6grad002.6e4rangez028
setgravity00-10
initialdensity2.6e3
initialsxx-45.8e6grad002.6e4rangez028
initialsyy-45.8e6grad002.6e4rangez028
---------------------------------------------
-----------材料属性--------------------------
-----------No1-----------------------
probulk8e9shear7e9fric25coh7.5e6ten2.6e6rangez08
-----------No2-------------------------
probulk7.8e9shear3.5e9fric30coh6.5e6ten2.4e6rangez810
------------No3---------------------
probulk7e9shear3.0e9fric22coh5.5e6ten2.2e6rangez1012
------------No4-----------------------------
probulk6e9shear2.3e9fric26coh5.5e6ten2.0e6rangez1213
-----------No5-----------------------------------
probulk6.2e9shear2.0e9fric28coh5.2e6ten1.2e6rangez1314
-------------------No6------------------------
probulk5.8e9shear2.5e9fric35coh6.0e6ten1.3e6rangez1416
-------------------No7------------------------
probulk9e9shear3.2e9fric30coh10.3e6ten1.2e6rangez1617
-------------------No8------------------------
probulk8.1e9shear3.1e9fric26coh8.0e6ten1.5e6rangez1719
-------------------No9------------------------
probulk7.6e9shear2.6e9fric23coh8.5e6ten1.8e6rangez1923
-------------------No10------------------------
probulk7.5e9shear2.1e9fric26coh8.0e6ten1.9e6rangez2328
int1propks8e10kn8e10fric30dil6tension2e7bslip=on
int2propks8e10kn8e10fric30dil6tension2e7bslip=on
int3propks8e10kn8e10fric35dil8tension2e7bslip=on
int4propks9e11kn9e11fric32dil5tension2e7bslip=on
int5propks9e11kn9e11fric32dil5tension2e7bslip=on
int6propks8e11kn8e11fric33dil5tension2e7bslip=on
int7propks8e10kn8e10fric33dil5tension2e7bslip=on
int8propks8e10kn8e10fric33dil5tension2e7bslip=on
int9propks8e10kn8e10fric33dil5tension2e7bslip=on
step5000
这是给的模型加的初始地应力，难道这样加地应力，模型就自动平衡了吗？
你用ini平衡，为什么还要计算5000步？
我的INI是用来施加初始地应力的，难道用INI那样加应力边界条件就不用step或solve算平衡了吗？
模型就能自动平衡吗？
把上面的sxx和syy的剃度2.6e4改为9e3试一下，xx和yy方向的梯度不可能与zz方向的梯度一样的。
7.问：
propmat和zone定义的材料特性有什么区别？
前者是对整个模型（如果整个模型只分成一个region的话），后者是对具体的zone的。
8.关于ff边界设置的问题
自由场地边界（FF边界）在动力分析中经常采用，它的原理是采用粘滞阻尼器来模拟静止边界，将自由场节点的不平衡力加到主体的网格边界上。
注意：
模型底部的动力边界条件应当在appff之前，进行appff以后底部的边界条件就自动转为自由场地边界。
现举出一简单例子,比如静力模型的尺寸为10*10*10,有以下几种边界施加方式
1）
freexyzrangz-.1.1
applysquietdquietnquietrangz-.1.1
applynstressvaluehistwaverangz-.1.1
applyff
2）
3）
可以验证,以上的施加方式得到的结果是完全一样的.所谓的自由场边界最重要的特征就是复制特征,现对以上各例子做简要分析.
1）第二三两句就是在10*10的底面范围内分别设置静止边界和动力条件,当加了ff后,ff的网格的条件即从临近的maingrid的zone进行复制.所以ff的底面的位移边界和动力和maingrid完全一样.
2）第二句直接加ff边界,注意这时ff底面是free的,而且动力为0,但是接下来的3.4两句仅仅给定的是z=-.1至.1的范围,导致ff的底面也包括在内了,所以作用方式和命令流1完全相同.
3）同2的分析,同样可以得到结论但是对于以上三个命令流中的两句
applynst
升级会员 