模型 边界 受力对称，位移结果却不对称！求解！！

hewitthy

rest chushi.sav
plot block group
plot add axes red
plot set rotation 0 0 45 ;用于显示三维模型
gen separate 3
gen separate 1
ini x add -100 rang group 3
ini x add 100 rang group 1
;
;;interface施加
int 1 face range plane ori -185 0 0 nor 1.732 0 1 dist 0.1
int 2 face range plane ori 45 0 0 nor -1.732 0 -1 dist 0.1
plot add interface yellow
;
ini x add 100 rang group 3
ini x add -100 rang group 1
attach face range cylinder end1 0 0 0 end2 -160 0 0 radius 7.6
;设置重力
set gravity 0 0 -10
;给定边界条件
fix z range z -19.59,-19.61
fix y range y -39.61,-39.59
fix y range y 39.61,39.59
fix x range x -0.01,0.01
fix x range x -159.99,-160.01
;求解自重应力场
;参数设置
def para_1
   bulk01_1=elas01/(3*(1-2*v01))
   shear01_1=elas01/(2*(1+v01))
   k_s01_1=10*(bulk01_1+4*shear01_1/3)/z
end
   set elas01=1.5e9
   set v01=0.31
   set z=1.2
para_1
def para_2
   bulk01_2=elas02/(3*(1-2*v02))
   shear01_2=elas02/(2*(1+v02))
end
   set elas02=0.25e9
   set v01=0.45
para_2
model mohr  
hist unbal
;int 1 prop ks=k_s01_1 kn=k_s01_1 fric 20         ;ks=kn=10*max[(K+4G/3)/z],z=1.2
;int 2 prop ks=k_s01_1 kn=k_s01_1 fric 20
ini density 2000  ;围岩的密度
prop bulk=shear01_1 shear=shear01_1 fric=30 coh=2.5e5 tension=1.0e5 &    ;体积、剪切、摩擦角、凝聚力、抗拉强度
                   range group 2 not
prop bulk=shear01_2 shear=shear01_2 fric=22 coh=1.5e5 tension=7e4 range group 2 any
set mech ratio=4e-4
step 7000
save inisolve.sav
plot cont zdisp outl on
pause
initial xdisplacement 0 ydisplacement 0 zdisplacement 0
initial xvelocity 0 yvelocity 0 zvelocity 0
define excav
loop cut (1,45)
      i=cut
      y0=-160+(i-1)*4
      y1=-160+4+(i-1)*4
     y0p=y0-20   
     y1p=y1-20
;
;开挖
if i>=1 then
    if i<=40  then
         command
             model  null  range  group 5 x @y0 @y1
             model  null  range  group 6 x @y0 @y1
         end_command
    endif
endif
;加固
if i>=1 then  
if i<=40 then
      command
       model mohr range group 4 x @y0 @y1
       ini density 2500  range group 4 x @y0 @y1   ;注浆圈的密度
         prop bulk=5e9 shear=1.969e9 fric=30 coh=2.5e5 tension=1.0e5 range group 4 x @y0 @y1 ;注浆圈的体积弹模、剪切弹模
      end_command
endif
endif

;支护
if i>=6 then  
     if i<=45 then
      command
          model elastic range group 6 x @y0p @y1p
          ini density 2700  range group 6 x @y0p @y1p   ;衬砌混凝土的密度
           prop bulk=0.77e9,shear=0.61e9 range group 6 x @y0p @y1p ;衬砌混凝土的体积弹模、剪切弹模
      end_command
     endif
endif
;
if i<=32 then
     command
    step 2000
     end_command
end if   
if i>35 then
     command
     step 2500
     end_command
end if
    ;  
;存储指定步数文件
if i=21 then
   command
   save meter_84.sav
   end_command
   endif   
  if i=45 then
   command
   save meter_erchen_160.sav
   end_command
   endif
end_loop
end
@excav