flac2d用修正剑桥试算一个小模型遇到问题

wuyong505

Project Record Tree export

;... STATE: GRID ....
config
grid 20 10;;;;;;;attention to 20!
gen 0,0 0,5 8,5 8,0 1 i 1 9 j 1 6
mo elas i=1,8 j=1,5
gen 8,0 8,5 10,5 10,0 i 10 12 j 1 6
mo elas i=10,11 j=1,5
gen 10,0 10,5 18,5 18,0 i 13 21 j 1 6
mo elas i=13,20 j=1,5
gen 0,5 0,10 8,10 8,5 i 1 9 j 6 11
mo elas i=1,8 j=6,10
gen 8,5 8,10 10,10 10,5 i 10 12 j 6 11
mo elas i=10,11 j=6,10
gen 10,5 10,10 18,10 18,5 i 13 21 j 6 11
mo elas i=13,20 j=6,10
;;;;;;;;;;;;;;;;;;;;;
struct node 1 8,10
struct node 2 8,5
struct beam begin node 1 end node 2 prop 1001 seg 5
struct prop 1001 e 2.4E10 area 1.0 I 0.0833
struct node 7 10,10
struct node 8 10,5
struct beam begin node 7 end node 8 prop 1001 seg 5
struct prop 1001 e 2.4E10 area 1.0 I 0.0833
;;;;;;;;;;;;
;;;;;;;;;;;
attach aside from 9 6 to 9 1 bside from 10 6 to 10 1
attach aside from 12 6 to 12 1 bside from 13 6 to 13 1
;;;;;;;;;;;;
interface 1 aside from 9,11 to 9,6 bside from node 1 to node 2
interface 2 aside from 10,11 to 10,6 bside from node 1 to node 2
interface 1 unglued kn=4.645E8 ks=4.645E8 cohesion=0.0 dilation=0.0 &
friction=12.5 tbond=0 sbratio=20 bslip=off
interface 2 unglued kn=4.645E8 ks=4.645E8 cohesion=0.0 dilation=0.0 &
friction=12.5 tbond=0 sbratio=20 bslip=off
;;;;;;;;;;;;
interface 3 aside from 12,11 to 12,6 bside from node 7 to node 8
interface 4 aside from 13,11 to 13,6 bside from node 7 to node 8
interface 3 unglued kn=4.645E8 ks=4.645E8 cohesion=0.0 dilation=0.0 &
friction=12.5 tbond=0 sbratio=20 bslip=off
interface 4 unglued kn=4.645E8 ks=4.645E8 cohesion=0.0 dilation=0.0 &
friction=12.5 tbond=0 sbratio=20 bslip=off
;;;;;;;;;;;
fix x i 1
fix x i 21
fix x y j 1
save grid.sav

;... STATE: INI ....
model elastic notnull
prop density=1800.0 bulk=1.66667E7 shear=5555560.0 notnull j 1 10;;考虑到接触面参数和收敛速度
;;;;;;;;;;;;;;
ini syy -180e3 var 0,180e3 j 1 10
;;;;;;;;;
def ini_k0
loop i (1,izones)
    loop j (1,jzones)
      k0=0.5
      sxx(i,j)=k0*syy(i,j)
      szz(i,j)=sxx(i,j)
    endloop
endloop
end
ini_k0
;;;;;;;;;;;;;;;
set grav=10
history 999 unbalanced
solve
save ini.sav

;... STATE: MODELCAM ....
;;;;;;;;;;;;;;;;;;
ini xdisp 0 ydisp 0
ini xvel 0 yvel 0
;;;;;;;;;;;;;;;;;;;
stru node range 1 12 ini xdis 0 ydis 0
stru node range 1 12 ini xvel 0 yvel 0
;;;;;;;;;;;;;;;;;;
model cam-clay notnull
prop density=1800.0 bulk=1e8 kappa=0.01 lambda=0.14 mpc=180e3 mm=1.2 mp1=1000.0 &
mv_l=3.00 poiss=0.35 notnull
save modelcam.sav

;... STATE: SOLVE ....
;;;;;;;;;;;;;;
def ini_k0
loop i (1,izones)
    loop j (1,jzones)
      k0=0.5
      sxx(i,j)=k0*syy(i,j)
      szz(i,j)=sxx(i,j)
    endloop
endloop
end
ini_k0
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
def camclay_ini_p
loop i (1,izones)
   loop j (1,jzones)
     if model(i,j)=11 then
       OCR=1.0
       s1=abs(sxx(i,j))
       s2=abs(syy(i,j))
       s3=abs(szz(i,j))
       p0=(s1+s2+s3)/3.0
       cam_p(i,j)=p0
       q0=sqrt(((s1-s2)^2+(s2-s3)^2+(s3-s1)^2)*0.5)
       temp1=q0/(mm(i,j)*p0)
       pc0=p0*(1.0+temp1^2)*OCR
       mpc(i,j)=pc0
     endif
   endloop
endloop
end
camclay_ini_p
;;;;;;;;;;;;;
solve
save solve.sav

;*** plot commands ****
;plot name: syy
plot hold syy fill
;plot name: sxx
plot hold sxx fill
;plot name: szz
plot hold szz fill
;plot name: Unbalanced force
plot hold history 999
;plot name: g
plot hold  grid magnify 20.0 lred grid
;plot name: Plot 2
plot hold num
;plot name: mpc
plot hold mpc block
;plot name: state
plot hold state block

结果在连续墙接触面位置出现这样的情况：

到底是什么原因呢？？
是由于模型改变导致的瞬间效应破坏吗？
调整剑桥模型的mpc（调大），可以初始平衡，但是开挖后结果不收敛。

gnik77

楼主您好，能否将FLAC5.0发予我一份？
我安装的版本老出现“the hardware key was not found ,nor FLACS.CFG file. also, information that FLAC was run with a key on this computer was not found..”这个问题，不知道对以后的计算使用会不会产生影响~~
谢谢！