new
;===================生成壳单元==============
gen zo b p0=(-3.5 0 0) p1=(0 0 0) p2=(-3.5 20 0) p3=(-3.5 3.14 7.85) p4=(0 20 0) p5=(-3.5 20 7.85) p6=(0 3.14 7.85) p7=(0 20 7.85) size 7 40 16 group penshehunningtu
pl add surface red
pl show
sel shell id=1 elemtype=dkt crossdiag range plane above normal 0,-1,0.4 origin -1.75,0.108,0 & plane below normal 0,-1,0.4 origin -1.75,0.62,1.55
delete
pl add sel geom node=off shrinkfac=0 magfac=0 red red
sel group shotceng range z=0 7.85
