命令流如下:
new
genzonecylp1100p2040p3001size123012
genzonereflenorm100
genzonereflenorm001
moss
propden2500bulk2e8she1e8co2e6fric45ten1e6dil10
propftab1ctab2dtab3
table10450.05420.140140
table202e60.051e60.15e515e5
table30100.0530.10
fixxrangeend1100end2140rad1end10.500end20.540rad0.5
回答:
加上的约束不对!
fixxrangeend1100end2140rad1end10.500end20.540rad0.5
回答:
上面的rad是不是相差太大了?
是不是在1左右来限制了
2.关于计算不收敛问题
;初始网格
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_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&
size882rat1.021.021.02
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&
size881rat1.021.021.02
genzonereflectdip90dd90origin(0,0,1)rangez1415
grouplayer4rangez1415
interface4facerangeplanenorm(0,0,1)origin(0,0,13)dist0.01
inizadd-1.0rangegrouplayer4
;No5
genzonebrickp00015p110015p201015p30017&
size882rat1.021.021.02
genzonereflectdip90dd90origin(0,0,15)rangez1517
grouplayer5rangez1517
interface5facerangeplanenorm(0,0,1)origin(0,0,14)dist0.01
inizadd-1.0rangegrouplayer5
;No6
genzonebrickp00017p110017p201017p30018&
size881rat1.021.021.02
genzonereflectdip90dd90origin(0,0,1)rangez1718
grouplayer6rangez1718
interface6facerangeplanenorm(0,0,1)origin(0,0,16)dist0.01
inizadd-1.0rangegrouplayer6
;No7
genzonebrickp00018p110018p201018p30020&
size882rat1.021.021.02
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
;---------------------------------------------
modelmohr
;-----------材料属性--------------------------
;-----------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
histunbal
step5000
applyszz-26e6rangez27.928.1
initialszz-1.08e6grad002.6e4rangez028
setgravity00-10
initialdensity2.6e3
initialsxx-45.8e6grad002.6e4rangez028
initialsyy-45.8e6grad002.6e4rangez028
initialszz-1.08e6grad002.6e4rangez028
initialdensity2.6e3
initialsxx-45.8e6grad002.6e4rangez028
initialsyy-45.8e6grad002.6e4rangez028
这是给的模型加的初始地应力,难道这样加地应力,模型就自动平衡了吗?
回答:
你用ini平衡,为什么还要计算5000步?
问:
我的INI是用来施加初始地应力的,难道用INI那样加应力边界条件就不用step或solve算平衡了吗?
模型就能自动平衡吗?
;-----------初始条件-------------------------
applyszz-26e6rangez27.928.1
initialszz-1.08e6grad002.6e4rangez028
setgravity00-10
initialdensity2.6e3
initialsxx-45.8e6grad002.6e4rangez028
initialsyy-45