塑性区为什么超大？

haiyangzhou

加了动力模块，塑性区变得超大，不知道是为什么？高手指点一下！
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