3D的动力(爆破)分析问题FLAC

admin管理员组

文章数量:1530914

2024年7月25日发(作者：)

3D的动力（爆破）分析问题

我刚在做一个隧道爆破掘进的模拟，想这样做：

q0HG tZ9u!tIl

1.先模拟原岩应力状态

-V0pvT

mo mohr

v)p4FvX

;岩体参数

V%ws~'c r}T!]1r7F+w

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

{+u(gHf T^k

ini density 1650 range z 17.2 13.6

R4kg8ek'_3b

pro bulk 9.9e6 she 2.82e6 fric 29 coh 0 ten 0 range z 13.6 7.1

E l Ny U7T X"K K+Pt

ini density 1900 range z 13.6 7.1

f/IL;I)s1-g.L

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

J-q^!e&c !c

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

HB$O-g}gA#tM&i

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

PS]r0Vu0v4N/{

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自由)

Y$o?W:e!Lm

fix z range z -26.8 -27

0?cN9c Hm

fix x range x -39.9 -40.1

fix x range x 39.9 40.1

,d%l2mJP QX"YR

fix y range y -0.1 0.1

fix y range y 39.9 40.1

9s@u @-e&[N:D3z

set grav 0 0 -9.8

pj;w#^J_?3L

;施加地面面状荷载

#O+L t6v)~g

apply szz -50e3 range z 18 18.2

sol rat 1e-5

save

2.模拟爆破掘进

new

&L c-C1ayu)W,EuM

rest

]]^+d!B2W4@

config dyn

set dyn on

ini xdisp 0 ydisp 0 zdisp 0

ini xvel 0 yvel 0 zvel 0

kTG$[0hNah-~

ini state 0

def setup

#Jt S~R~ {

omega = 2.0 * pi * freq

!F!d0Q(P[j8P

end

setup

|e4GW){3YC

;y向速度时程曲线

z -a8m'k

def wave

wave = 1/2.718282^(1.5*pi*dytime)*ampl*sin(omega * dytime) ;大概震动1.3s

end

set freq = 40.625 ampl = 0.04765

setup

;p1为隧道掘进面

apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01

N%y&l(ue;k1@)}5U/]O-a

;设定粘滞边界

8{VGE$G2S H;b D

apply xquiet yquiet ran z -26.8 -27

]q:A9d*O5g

apply xquiet yquiet ran x -39.9 -40.1

apply xquiet yquiet ran x 39.9 40.1

apply xquiet yquiet ran y -0.1 0.1

apply xquiet yquiet ran y 39.9 40.1

B4Y;L(t*c8I3a1p2QO

hist dytime

hist wave

set dyn damp local 0.15

a;DfR%]V1SH"R

sol age 1.0

m null range group p1 y 0 1

sol age 2.0

:B;z$n*j'i-J%d3U{-k

;之后就是循环的向前掘进，施加爆破波，掘进……

请大家帮忙看看这些命令流有没有什么错误

有这样几个问题：

（用flac没多久，也是刚刚才接触动力模块，所以应该都是很浅显的问题，希望大家不要不

屑于赐教哈）

1.施加动力荷载时，apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01中的yvel=？我看好多

例子似乎都是=1.0，不知道这个赋的值代表啥？

EN2g7@~O o5CA&]l

2.同样是这一句，apply yvel=0.0 hist=wave range gro p1 y -0.01 0.01，最后那施加的位置好像

多是一个面吧，能否施加在一定体积内呢？比如range gro p1 y 0 1？或者一条线上（线装炸

药）？

3.粘滞边界的设定跟之前静力计算中的约束条件有冲突么？比如fix z range z -26.8 -27与

apply xquiet yquiet ran z -26.8 -27？

r1C u"N9w

4.如果同时施加x、y、z三个方向的速度时程，是否就是这样写：

9|)~+oa$F

apply xvel=0.0 hist=wave1 range gro p1 y -0.01 0.01

O"p*|-^2P Bb

apply yvel=0.0 hist=wave2 range gro p1 y -0.01 0.01

apply zvel=0.0 hist=wave3 range gro p1 y -0.01 0.01

3@:^bl(p B

kU6[M(f;S?

主要就是设置粘滞边界条件和apply动力荷载这两个地方就是搞不清楚，希望版主或者明白

的朋友能给些指点，谢谢啦

本文标签： 爆破动力施加