mcssr.dat

一个基于Itasca Flac3D据强度折减法进行岩土边坡稳定分析的强度折减法fish命令

new
;====================================建模
def _geoparm
size1=20
size2=30
filecal0=string(size1)+'_'+string(size2)+'_cal'+'.sav'

bf=75   ;      改变此处坡角
af=bf*pi/180    ; 
h1=20
;
x00=30.0
x01=105.0
x02=x00
x03=x00+h1/tan(af)
x04=x01
x05=x03
x06=105.0
x07=x06
;
y00=0.0
y01=y00
y02=1.0
y03=y00
y04=y02
y05=y02
y06=y00
y07=y02
;
z00=20.0
z01=z00
z02=z00
z03=100.0
z04=z00
z05=z03
z06=z03
z07=z03
end
_geoparm

gen zone brick p0 0 0 0 p1 30 0 0 p2 0 1 0 p3 0 0 20 size 6 1 6 
gen zone brick p0 30 0 0 p1 105 0 0 p2 30 1 0 p3 30 0 20 size size1 1 6
gen zone brick p0 x00 y00 z00 p1 x01 y01 z01 p2 x02 y02 z02 p3 x03 y03 z03 &
               p4 x04 y04 z04 p5 x05 y05 z05 p6 x06 y06 z06 p7 x07 y07 z07 size size1 1 size2

;===============================================
   ;*************强度折减法*******************
;===============================================
def SSR

;------------------
ait1=0.02         ;精度
k11=1.00          ;上限值
k12=2.50          ;下限值
ks=(k11+k12)/2

loop while (k12-k11)>ait1

;------------------

  ; 弹性模量(Pa)
		flag=0
  E1=3.0e10   
  ; 泊松比
  poi1=0.27
  ;
  ; 粘结力(Pa)
  coh1=0.54e6/ks
  ;
  ; 容重(N/m3)
  weight1=27.0e3
  ; 膨胀角(度)
  dila1=0.0
  ;
  ; 内摩擦角(度)
  fri1=(atan((tan(50.1*pi/180))/ks))*180/pi
  ;
  ;抗拉强度
  ten1=0.05e6
  grav0=-9.80            ; 重力加速度(N/kg)
  dens1=-weight1/grav0   
  K1=E1/(3*(1-2*poi1)) 
  G1=E1/(2*(1+poi1))  
		
;==================

command

model null       ;挖空
model mohr  ;采用莫尔库仑模型
pro bulk K1 she G1 dens dens1 coh coh1 & 
friction fri1 dilation 0. tens ten1

fix x range x -0.1 0.1
fix x range x 104.9 105.1
fix y
fix x y z range z -0.1 0.1
set grav 0 0 grav0

set mech ratio 9.8e-6
solve step 30000
endcommand
;==================

;收敛
if mech_ratio<1.0e-5 
  k11=ks
		k12=k12

else 
  k12=ks
		k11=k11
endif
ks=(k11+k12)/2

endloop
fosfile0=string(size1)+'_'+string(size2)+'_fos'+'.sav'

command 
save fosfile0
endcommand
end
SSR

pr ks