- 积分
- 0
- 注册时间
- 2010-5-20
- 仿真币
-
- 最后登录
- 1970-1-1
|
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(调大),可以初始平衡,但是开挖后结果不收敛。
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|