3D动力分析（隧道爆破）求教

chaylo

求教。
1.原岩应力状态
mo mohr
;岩体参数
pro bulk 4.2e6 she 9.23e5 fric 12 coh 1e4 ten 9.6e3 range z 18.1 17.2
ini density 1800 range z 18.1 17.2
pro bulk 7.02e6 she 1.83e6 fric 7.2 coh 4.7e3 ten 4.7e3 range z 17.2 13.6
ini density 1650 range z 17.2 13.6
pro bulk 9.9e6 she 2.82e6 fric 29 coh 0 ten 0 range z 13.6 7.1
ini density 1900 range z 13.6 7.1
pro bulk 4.37e7 she 1.46e7 fric 30 coh 3e4 ten 3e4 range z 7.1 6.4
ini density 1990 range z 7.1 6.4
prop bulk 4.95e7 she 1.65e7 fric 38 coh 5e4 ten 5e4 range z 6.4 -0.4
ini density 2100 range z 6.4 -0.4
prop bulk 6.52e7 she 2.67e7 fric 40 coh 4e4 ten 4e4 range z -0.4 -9.1
ini density 2480 range z -0.4 -9.1
prop bulk 6.78e7 she 3.12e7 fric 43 coh 5e4 ten 5e4 range z -9.1 -17.6
ini density 2510 range z -9.1 -17.6
prop bulk 6.53e8 she 3.37e8 fric 44 coh 5e6 ten 5e6 range z -17.6 -26.9
ini density 2520 range z -17.6 -26.9
;设定边界条件(上边界z自由)
fix z range z -26.8 -27
fix x range x -39.9 -40.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 39.9 40.1
set grav 0 0 -9.8
;施加地面面状荷载
apply szz -50e3 range z 18 18.2
sol rat 1e-5
save geom.sav
2.模拟爆破掘进
new
rest geom_dyn.sav
;pl bl s
set log on
ini state 0
ini xdisp 0 ydisp 0 zdisp 0
ini xvel  0 yvel  0 zvel  0
;hist dis
;地表
hist id 1 gp zdisp 0  0  18.12
;拱顶
hist id 2 gp zdisp 0  0  -2.36
;侧帮
hist id 3 gp xdisp -3.2  0 -5.573
;拱底
hist id 4 gp zdisp 0  0 -9.08
;hist vel & acc
;地表
hist id 5 gp xvel 0 0 18.12
hist id 6 gp yvel 0 0 18.12
hist id 7 gp zvel 0 0 18.12
hist id 8 gp xacc 0 0 18.12
hist id 9 gp yacc 0 0 18.12
hist id 10 gp zacc 0 0 18.12
;洞顶
hist id 11 gp zvel 0 0 -2.36
hist id 12 gp zacc 0 0 -2.36
;侧帮
hist id 13 gp xvel -3.2 0 -5.573
hist id 14 gp xacc -3.2 0 -5.573
;动力模式:开
set dyn on
;设置粘滞(吸收)边界
apply nquiet squiet dquiet ran z -26.78 -26.98
apply nquiet squiet dquiet ran x -39.9 -40.1
apply nquiet squiet dquiet ran x 39.9 40.1
apply nquiet squiet dquiet ran y -0.1 0.1
apply nquiet squiet dquiet ran y 39.9 40.1
;定义应力时程曲线(爆破波)
;第一个爆破波
def wave
if dytime > 0.1
wave = 0.0
else
  if dytime < 0.000
  wave = 0
  else
  wave = 6e6*(2.718281828^(-74.5*dytime)-2.718281828^(-117.1*dytime))
  endif
endif
end
hist unbal
;整理时程曲线
hist dytime
hist wave
;设置局部阻尼
set dyn damp local 0.1571
;pl add hist 1 6 vs 72
;超前预加固
;prop bu 1.537e8 sh 6.9e7 f 42 coh 5e5 ten 4.86e5 density 2510 range group en
;第一个进尺
m null range group p1 y 0 2
m null range group li1 y 0 2
;假定爆破荷载以压力形式的均布荷载作用在隧道洞壁上,方向垂直于洞壁和掌子面
apply ns -1.0 hist=wave range plane normal ( -0.998 0 0.070 ) origin ( -3.20 0 -5.57 ) x -3.20 -3.17 z -5.57 -5.12 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.978 0 -0.208 ) origin ( -3.17 0 -5.12 ) x -3.17 -3.08 z -5.12 -4.68 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.940 0 -0.342 ) origin ( -3.08 0 -4.68 ) x -3.08 -2.92 z -4.68 -4.26 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.883 0 0.469 ) origin ( -2.92 0 -4.26 ) x -2.92 -2.70 z -4.26 -3.84 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.809 0 -0.588 ) origin ( -2.70 0 -3.84 ) x -2.70 -2.41 z -3.84 -3.50 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.707 0 0.707 ) origin ( -2.41 0 -3.46 ) x -2.41 -2.17 z -3.46 -3.22 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.602 0 -0.799 ) origin ( -2.07 0 -3.12 ) x -2.07 -1.73 z -3.12 -2.86 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.485 0 -0.875 ) origin ( -1.73 0 -2.86 ) x -1.73 -1.34 z -2.86 -2.66 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.358 0 0.934 ) origin ( -1.34 0 -2.66 ) x -1.34 -0.94 z -2.66 -2.50 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.225 0 0.974 ) origin ( -0.94 0 -2.56 ) x -0.94 -0.47 z -2.50 -2.40 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.070 0 0.998 ) origin ( -0.47 0 -2.40 ) x -0.47 -0.00 z -2.40 -2.36 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.035 0 -0.999 ) origin ( -3.20 0 -5.57 ) x -3.20 -2.90 z -5.57 -5.56 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.267 0 0.961 ) origin ( -2.41 0 -5.70 ) x -2.90 -2.41 z -5.56 -5.70 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.225 0 0.974 ) origin ( -1.91 0 -5.82 ) x -2.41 -1.91 z -5.70 -5.82 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.174 0 0.985 ) origin ( -1.40 0 -5.91 ) x -1.91 -1.40 z -5.82 -5.91 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.1045 0 0.9945 ) origin ( -0.70 0 -5.99 ) x -1.40 -0.70 z -5.91 -5.99 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.035 0 -0.999 ) origin ( -0.00 0 -6.01 ) x -0.70 0.00 z -5.99 -6.01 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.998 0 0.07 ) origin ( 3.2 0 -5.57 ) x 3.2 3.17 z -5.57 -5.12 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.978 0 -0.208 ) origin ( 3.17 0 -5.12 ) x 3.17 3.08 z -5.12 -4.68 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.94 0 -0.342 ) origin ( 3.08 0 -4.68 ) x 3.08 2.92 z -4.68 -4.26 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.883 0 0.469 ) origin ( 2.92 0 -4.26 ) x 2.92 2.7 z -4.26 -3.84 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.809 0 -0.588 ) origin ( 2.7 0 -3.84 ) x 2.7 2.41 z -3.84 -3.5 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.707 0 0.707 ) origin ( 2.41 0 -3.46 ) x 2.41 2.17 z -3.46 -3.22 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.602 0 -0.799 ) origin ( 2.07 0 -3.12 ) x 2.07 1.73 z -3.12 -2.86 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.485 0 -0.875 ) origin ( 1.73 0 -2.86 ) x 1.73 1.34 z -2.86 -2.66 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.358 0 0.934 ) origin ( 1.34 0 -2.66 ) x 1.34 0.94 z -2.66 -2.5 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.225 0 0.974 ) origin ( 0.94 0 -2.56 ) x 0.94 0.47 z -2.5 -2.4 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.07 0 0.998 ) origin ( 0.47 0 -2.4 ) x 0.47 0 z -2.4 -2.36 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.035 0 -0.999 ) origin ( 3.2 0 -5.57 ) x 3.2 2.9 z -5.57 -5.56 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.267 0 0.961 ) origin ( 2.41 0 -5.7 ) x 2.9 2.41 z -5.56 -5.7 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.225 0 0.974 ) origin ( 1.91 0 -5.82 ) x 2.41 1.91 z -5.7 -5.82 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.174 0 0.985 ) origin ( 1.4 0 -5.91 ) x 1.91 1.4 z -5.82 -5.91 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( -0.1045 0 0.9945 ) origin ( 0.7 0 -5.99 ) x 1.4 0.7 z -5.91 -5.99 y 0 1.9
apply ns -1.0 hist=wave range plane normal ( 0.035 0 -0.999 ) origin ( 0.00 0 -6.01 ) x 0.7 0 z -5.99 -6.01 y 0 1.9
apply ns -1.0 hist=wave range gro p1 y 1.9 2.1
apply ns -1.0 hist=wave range gro li1 y 1.9 2.1
;pl add hist 15 vs 16
solve age 0.1
save dyn0.1.sav
;sol
求教。

lovemm

有点大！

chaylo

有点大！
lovemm 发表于 2009-4-1 10:09

啥子有点大哦？

gcchfut

什么问题，请指明

caikewei1984

学习学习

chaylo

4# gcchfut 新发了个帖子：
http://forum.simwe.com/thread-873406-1-1.html