用mohr——coulomb和hoek——brown做的煤矿采空区

caihuaien

pl add block group
pl add axes black
pl show
set log on
set logfile huangtumodel.log
;建立原始几何模型

gen zone brick p0 0 0 0 p1 400 0 0 p2 0 400 0 p3 0 0 10.5            &
               size 20 20 1
               group  石英砂岩 red

gen zone brick  p0 0 0 10.5 p1 400 0 10.5 p2 0 400 10.5 p3 0 0 12.5   &
                size 20 20 1
                group 煤 black range group 石英砂岩 not
;

gen zone brick p0 0 0 12.5 p1 400 0 12.5 p2 0 400 12.5 p3 0 0 21      &
               size 20 20 1
               group 砂岩 blue  range group 煤 not  group 石英砂岩 not


gen zone brick p0 0 0 21 p1 400 0 21 p2 0 400 21 p3 0 0 71.5          &
               size 20 20 5
               group 页岩 yellow range group 砂岩 not group 煤 not  group 石英砂岩 not
;
gen zone brick p0 0 0 71.5 p1 400 0 71.5 p2 0 400 71.5 p3 0 0 90      &
               size 20 20 3
               group  中细砂岩 orange range group 页岩 not group 砂岩 not group 煤 not &
               group 石英砂岩 not
;
gen zone brick p0 0 0 90 p1 400 0 90 p2 0 400 90 p3 0 0 93.5      &
               size 20 20 1
               group 页岩 green range group  中细砂岩 not group 页岩 not &
               group 砂岩 not group 煤 not group 石英砂岩 not
;
gen zone brick p0 0 0 93.5 p1 400 0 93.5 p2 0 400 93.5 p3 0 0 98      &
               size 20 20 1
               group 中细砂岩 lgreen range group 页岩 not group  中细砂岩 not group 页岩 not &
               group 砂岩 not group 煤 not group 石英砂岩 not
;
gen zone brick p0 0 0 98 p1 400 0 98 p2 0 400 98 p3 0 0 138     &
               size 20 20 4
               group 黄土 yellow range group 中细砂岩 not  group 页岩 not group  中细砂岩 not &
               group 页岩 not group 砂岩 not group 煤 not group 石英砂岩 not
;
;初始属性

model hoekbrown
;range group 中细砂岩 group 页岩 group 砂岩 group 煤 group 石英砂岩
  prop dens 2002 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 黄土
  prop dens 2002 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 中细砂岩
  prop dens 2003 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 页岩
  prop dens 2004 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 中细砂岩
  prop dens 2005 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 页岩
  prop dens 2006 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 砂岩
  prop dens 2007 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 煤
  prop dens 2008 bulk=3.66700006E9 shear=2.2E9 hbs=0.0039 hbmb=1.7 hbsigci=3.0E7 hba=0.5 hbs3cv=0.0 range group 石英砂岩

;model mohr range group 黄土
;  prop dens 2000 tens 1e9 fric 30 coh 1e9 bulk 1e9 shear 1e9 range group 黄土

;固定边界
  fix x range x -0.1   0.1
  fix x range x 399.9  400.1
  fix y range y -0.1   0.1
  fix y range y 399.9  400.1
  fix x y z range z -0.1   0.1
; 加应力变化值
set gravity 0  0 -10
  set large
   set mech r=5e-4
   hist n=5
   hist unbal
solve
save initial_balance.sav
initial xdis 0 ydis 0 zdis 0 state 0
save inti_disp0.sav
def shux

; 密度
   des1=1786
   des2=2630
   des3=2640
   des4=2630
   des5=2640
   des6=2625
   des7=1540
   des8=2745

;弹性模量
   young1=1.0e7
   young2=1.09e10
   young3=5.50e9
   young4=1.09e10
   young5=5.5e9
   young6=1.04e10
   young7=1.00e7
   young8=1.653e10

;泊松比
   psb1=0.30
   psb2=0.20
   psb3=0.25
   psb4=0.19
   psb5=0.25
   psb6=0.20
   psb7=0.36
   psb8=0.15

;内聚力
   cn1=1.9e4

;摩擦角
   fda1=14
   

; 体积模量
  bm1=float(young1/(3*(1-2*psb1)))
  bm2=float(young2/(3*(1-2*psb2)))
  bm3=float(young3/(3*(1-2*psb3)))
  bm4=float(young4/(3*(1-2*psb4)))
  bm5=float(young5/(3*(1-2*psb5)))
  bm6=float(young6/(3*(1-2*psb6)))
  bm7=float(young7/(3*(1-2*psb7)))
  bm8=float(young8/(3*(1-2*psb8)))

;剪切模量
   js1= float(young1/(2*(1+psb1)))
   js2= float(young2/(2*(1+psb2)))
   js3= float(young3/(2*(1+psb3)))
   js4= float(young4/(2*(1+psb4)))
   js5= float(young5/(2*(1+psb5)))
   js6= float(young6/(2*(1+psb6)))
   js7= float(young7/(2*(1+psb7)))
   js8= float(young8/(2*(1+psb8)))


;抗压强度
  _hbsigci2= 9.5e07  
  _hbsigci3= 4.5e07
  _hbsigci4= 9.5e07
  _hbsigci5= 4.5e07
  _hbsigci6= 8.0e07
  _hbsigci7= 0.8e07
  _hbsigci8= 1.5e08
  
;m
  m2= 7
  m3= 4.5
  m4= 7
  m5= 4.5
  m6= 7
  m7= 1
  m8= 7

;s
  s2 = 0.08
  s3 = 0.04
  s4 = 0.08
  s5 = 0.04
  s6 = 0.05
  s7 = 0.004
  s8 = 0.085

;抗拉
  ts1=1.80e3
  
;hbs3cv
   _hbs3cv2=0
   _hbs3cv3=0
   _hbs3cv4=0
   _hbs3cv5=0
   _hbs3cv6=0
   _hbs3cv7=0
   _hbs3cv8=0


a=0.5

     
end
shux

  m  m
  prop dens des1 tension ts1 fri fda1 coh cn1 bulk bm1 shear js1 range group 黄土

  model hoekbrown
  prop dens des2 hbsigci _hbsigci2 hbs3cv _hbs3cv2 hbmb m2 hbs s2 hba a bulk bm2 shear js2  range group 中细砂岩
  prop dens des3 hbsigci _hbsigci3 hbs3cv _hbs3cv3 hbmb m3 hbs s3 hba a bulk bm3 shear js3  range group 页岩
  prop dens des4 hbsigci _hbsigci4 hbs3cv _hbs3cv4 hbmb m4 hbs s4 hba a bulk bm4 shear js4  range group 中细砂岩
  prop dens des5 hbsigci _hbsigci5 hbs3cv _hbs3cv5 hbmb m5 hbs s5 hba a bulk bm5 shear js5  range group 页岩
  prop dens des6 hbsigci _hbsigci6 hbs3cv _hbs3cv6 hbmb m6 hbs s6 hba a bulk bm6 shear js6  range group 砂岩
  prop dens des7 hbsigci _hbsigci7 hbs3cv _hbs3cv7 hbmb m7 hbs s7 hba a bulk bm7 shear js7  range group 煤
  prop dens des8 hbsigci _hbsigci8 hbs3cv _hbs3cv8 hbmb m8 hbs s8 hba a bulk bm8 shear js8  range group 石英砂岩


;开采
;model null range x 150 250 y 100 150 group 煤
step 1000
save excavtemp2.sav
set log off


我这样写代码，老是出错，那位大侠帮我看看，是怎么回事，主要是在一个模型中用两个屈服准则，即黄土用mohr准则，岩体用hoek－brown准则，参数先不管，看看代码过程写的那里有问题！

wolf2005

最好是以附件形式贴上来，这样大家才好运行
还有，出错信息是什么？

zhaojiegj

还是上传附件吧！

caihuaien

问题解决了
model mohr range group 黄土
model hoekbrown  range group 中细砂岩
  model hoekbrown  range group 页岩
  model hoekbrown  range group 砂岩
  model hoekbrown  range group 煤
  model hoekbrown  range group 石英砂岩
  这样就可以了。

caihuaien

另外问一下，附件怎么上传。

zhangboseu_2001

发帖子的时候就有让选择附件的选项，跟发邮件一样

yuanyuekafu

希望lz上传命令附件啊 谢谢！

hanhaishenlun

问一下，hoekbrown模型赋参数值时hbsigci hbs3cv  hbmb  hbs hba 代表什么意思？

luobotou88

楼主，请问 hoekbrown模型赋参数值时hbsigci hbs3cv  hbmb  hbs hba 代表什么