加了动力模块,塑性区变得超大,不知道是为什么?高手指点一下!
new
rest 圆形隧洞.sav
config dynamic
; 本构模型和参数设置
model m
def prop_material
valu_bulk= modul / (3.0*(1-2*p_ratio))
valu_shear= modul / (2*(1+p_ratio))
end
set modul=60e9 ;弹模
set p_ratio=0.25 ;泊松比
prop_material
prop bul valu_bulk shea valu_shear
prop cohesion 30e6 fric 20 dilation 30 ten 1e6
ini dens=2700
;应力边界条件
ini sxx -50e6 syy -40e6 szz -10e6
apply sxx -50e6 range x 99.9 100.1
apply sxx -50e6 range x -99.9 -100.1
apply szz -10e6 range z 99.9 100.1
apply szz -10e6 range z -99.9 -100.1
apply syy -40e6 range y -0.1 0.1
apply syy -40e6 range y 4.9 5.1
fix y
solve
;========清零
ini xdisp 0.0 ydisp 0.0 zdisp 0.0
ini xvel 0 yvel 0 zvel 0
def state_ini
pnt=zone_head
loop while pnt # null
z_state(pnt,0)=0
pnt=z_next(pnt)
endloop
end
state_ini
save 地应力模型.sav
;=====打开动力模块
set dyn on
apply nquiet squiet dquiet range x 99.9 100.1
apply nquiet squiet dquiet range x -99.9 -100.1
apply nquiet squiet dquiet range y 99.9 100.1
apply nquiet squiet dquiet range y -99.9 -100.1
apply nquiet squiet dquiet range z -99.9 -100.1
apply nquiet squiet dquiet range z 99.9 100.1
set dyn damp rayl 0.005 10 ;stiff
model null range group 2
;=====加爆破力
apply nstr -1 range cylinder end1 0 0 0 end2 0 5 0 radius 7.501
solve age=0.1
save 爆破.sav