用壳单元建立薄壁圆筒，这样有什么问题吗？

善良的喵

new
  
gen zone cshell p0 0 0 575 p1 0.2 0 575 p2 0 0 401 p3 0 0.2 575 p4 0.2 0 401 p5 0 0.2 401 p8 0.1 0 575  p9 0 0.1 575 p10 0.1 0 401 p11 0 0.1 401 size 1 58 2 1 group neiguan
gen zone cshell p0 0 0 401 p1 0.2 0 401 p2 0 0 310 p3 0 0.2 401 p4 0.2 0 310 p5 0 0.2 310 p8 0.1 0 401  p9 0 0.1 401 p10 0.1 0 310 p11 0 0.1 310 size 1 10 2 1 group neiguan

gen zone cylinder p0 0 0 575 p1 0.1 0 575 p2 0 0 401 p3 0 0.1 575 p4 0.1 0 401 p5 0 0.1 401  size 1 58 2 group zhongxinshang
gen zone cylinder p0 0 0 401 p1 0.1 0 401 p2 0 0 310 p3 0 0.1 401 p4 0.1 0 310 p5 0 0.1 310  size 1 10 2 group zhongxinxia

attach face range group neiguan any group zhongxinshang any group zhongxinxia any


gen zone radc p0 0 0 275 p1 50 0 275 p2 0 0 101 p3 0 50 275 p4 50 0 101 p5 0 50 101 p6 50 50 275 p7 50 50 101 p8 3.2 0 275 p9 0 3.2 275 p10 3.2 0 101 p11 0 3.2 101 size  16 58 25 10
gen zone radc p0 0 0 101 p1 50 0 101 p2 0 0 10 p3 0 50 101 p4 50 0 10 p5 0 50 10 p6 50 50 101 p7 50 50 10 p8 3.2 0 101 p9 0 3.2 101 p10 3.2 0 10 p11 0 3.2 10 size 16 10 25 10
gen zone radc p0 0 0 10 p1 50 0 10 p2 0 0 0 p3 0 50 10 p4 50 0 0 p5 0 50 0 p6 50 50 10 p7 50 50 0 p8 3.2 0 10 p9 0 3.2 10 p10 3.2 0 0 p11 0 3.2 0 size 16 1 25 10

                    
gen zone cshell p0 0 0 275 p1 3.2 0 275 p2 0 0 101 p3 0 3.2 275 p4 3.2 0 101 p5 0 3.2 101 p8 0.5 0 275  p9 0 0.5 275 p10 0.5 0 101 p11 0 0.5 101 size  14 58 25 3
gen zone cshell p0 0 0 101 p1 3.2 0 101 p2 0 0 10 p3 0 3.2 101 p4 3.2 0 10 p5 0 3.2 10 p8 0.5 0 101  p9 0 0.5 101 p10 0.5 0 10 p11 0 0.5 10 size 14 10 25 3
gen zone cshell p0 0 0 10 p1 3.2 0 10 p2 0 0 0 p3 0 3.2 10 p4 3.2 0 0 p5 0 3.2 0 p8 0.5 0 10  p9 0 0.5 10 p10 0.5 0 0 p11 0 0.5 0 size 14 1 25 3


gen zone cshell p0 0 0 275 p1 0.5 0 275 p2 0 0 101 p3 0 0.5 275 p4 0.5 0 101 p5 0 0.5 101 p8 0.2 0 275  p9 0 0.2 275 p10 0.2 0 101 p11 0 0.2 101 size  2 58 4 1
gen zone cshell p0 0 0 101 p1 0.5 0 101 p2 0 0 10 p3 0 0.5 101 p4 0.5 0 10 p5 0 0.5 10 p8 0.2 0 101  p9 0 0.2 101 p10 0.2 0 10 p11 0 0.2 10 size 2 10 4 1
gen zone cshell p0 0 0 10 p1 0.5 0 10 p2 0 0 0 p3 0 0.5 10 p4 0.5 0 0 p5 0 0.5 0 p8 0.2 0 10  p9 0 0.2 10 p10 0.2 0 0 p11 0 0.2 0 size 2 1 4 1


gen zone cshell p0 0 0 10 p1 0.2 0 10 p2 0 0 0 p3 0 0.2 10 p4 0.2 0 0 p5 0 0.2 0 p8 0.1 0 10  p9 0 0.1 10 p10 0.1 0 0 p11 0 0.1 0 size  1 1 2 1
gen zone cylinder p0 0 0 10 p1 0.1 0 10 p2 0 0 0 p3 0 0.1 10 p4 0.1 0 0 p5 0 0.1 0  size  1 1 2                                                                                   


group top range z 131 275  
group gas range cylinder end1 0 0 101 end2 0 0 111 radius 3.2  
group ungas range z 101 131 group gas not  
group bottom range z 0 101  


plo blo gro

ini z add -300 range  group neiguan any group zhongxinshang any group zhongxinxia any

gen zone reflect norm 0 -1 0 origin 0 0 0

gen zone reflect norm -1 0 0 origin 0 0 0



attach face range group neiguan any group zhongxinshang any group zhongxinxia any group top any group gas any group bottom any group ungas any  

ini x add 100 range  group top any group gas any group bottom any group ungas any group neiguan any group zhongxinshang any group zhongxinxia any



def make_sphere
    p_gp=gp_head
   loop while p_gp#null
     px=gp_xpos(p_gp)
  pz=gp_zpos(p_gp)  
   if px#0
     if pz<275.7 then
      gp_zpos(p_gp)=pz-px*5.2/100
     end_if
   end_if
    p_gp=gp_next(p_gp)
  end_loop
end
make_sphere

plot set plane ori 0 0 0 norm 0 1 0

model elas
prop density 1800 bulk 4e7 shear 1.43e7 range group top
prop density 1700 bulk 1e8 shear 4e7 range group gas
prop density 1700 bulk 1e8 shear 4e7 range group ungas
prop density 4500 bulk 3e13 shear 1e13 range group bottom

prop density 1800 bulk 4e7 shear 1.43e7 range group zhongxinshang
prop density 4500 bulk 4e7 shear 1.43e7 range group zhongxinxia

model null range group neiguan
sel shell id=1 range cylinder end1 0 0 10 end2 0 0 275 rad 0.2 group zhongxinshang not any group zhongxinxia not any
sel shell id=1 prop iso=(2.1e11,0.3) thick 0.1 density 7850



fix x y z range z -0.1 0.1
fix y range y -50.1 -49.9
fix y range y 49.9 50.1
fix x range x -0.1 0.1
fix x range x 99.9 100.1

set grav 0 0 -10
water dens 1000
water table ori 0 0 1000 norm 0 0 1

hist unbal
plo add hist 1
set mech ratio 1e-4
solve

plo con pp


app nstress -2.784e6 range  radius 0.5
hist unbal
plo add hist 1
set mech ratio 1e-4
solve

save 弹性平衡