请问，新版本中复合材料SHELL181结果与老版本中SHELL91结果差别很大的原因是什么？ 20

同一模型，SHELL181结果为-2.33，SHELL91结果为4.06，请问原因是什么？谢谢！

SHELL91命令流：

/PREP7 ET,1,SHELL91,,1,,,1,1,,1,1,1,,
KEYOPT,1,1,16, 
KEYOPT,1,2,1KEYOPT,1,5,1
KEYOPT,1,6,1
KEYOPT,1,8,1 
R,1 
RMODIF,1,1,3,0, , , ,0  
RMODIF,1,13,1,0,1.2,0,0,0,         
RMODIF,1,19,2,0,30,0,0,0,            
RMODIF,1,25,1,0,1.2,0,0,0,
 MP,EX,1,14000 
MP,PRXY,1,0.3 
DENS,1,1.9e-3
 
MP,EX,2,20
MP,PRXY,2,0.3
DENS,2,6e-5r=8500
percent=0.75      
dn=0.9
v=67     jichu=(percent*2-1)*r
jx=sqrt(r*r-jichu*jichu)
K, ,jx,,-jichu,
K, ,,,r,
k,,,,,
LARC,1,2,3,r,  
LDIV,1,dn, ,2,0
AROTAT,1,2, , , , ,3,2,180,4,lsel,s,,,1,9,2
LESIZE,all, , ,50, , , , ,1lsel,s,,,11,18,1
LESIZE,all, , ,10, , , , ,1lsel,s,,,2,10,2
LESIZE,all, , ,5, , , , ,1
 
amesh,allNSEL,S,LOC,Z,-jichu
d,all,,,,,,all,,,,,
NSEL,S,LOC,Y,0D,all, , , , , ,UY,ROTX,ROTZ, , , 
*GET,en,ELEM,,COUNT, , , , 
*do,i,1,en             esel,S,,,i   *GET,XE,ELEM,I,CENT,X*GET,YE,ELEM,I,CENT,Y          *GET,ZE,ELEM,I,CENT,Z  K,100, XE,YE,ZECSYS,2
*GET,zz,KP,100,LOC,y*GET,kyy,KP,100,LOC,Z
yy=(90-kyy)*3.1415926/180zz=zz*3.1415926/180q1=-1.3q2=0.1801*sin(yy)*cos(zz)q3=0.78*sin(yy)*sin(yy)q4=0.78*sin(yy)*sin(yy)*cos(2*zz)q5=0.25*0.5599*sin(yy)*sin(yy)*sin(yy)*cos(3*zz)q6=0.75*0.5599*sin(yy)*sin(yy)*sin(yy)*cos(zz)
f=v*v/16*2.125/100000*(q1+q2+q3+q4+q5+q6)sfe,i,1,pres,,-f,,,
*enddo
ACEL,0,0,9.8e-3,ALLSEL,ALL 
!* 
CSWPLA,11,2,1,1,
EMODIF,all,ESYS,11,FINISH 
/SOLANTYPE,0
 PSTRES,ON  SOLVE  
 
FINISH 
/SOLUTION  
ANTYPE,1
 
BUCOPT,LANB,3,0,0  /STATUS,SOLU
SOLVE  

shell181命令流：

/PREP7 ET,1,SHELL181  
KEYOPT,1,8,1EX,1,14000 
PRXY,1,0.3 
DENS,1,1.9e-3 
 
EX,2,20
PRXY,2,0.3 
DENS,2,6e-5MPLIST,ALL,,,EVLT  
sect,1,shell,, 
secdata, 1.2,1,0.0,3
secdata, 30,2,0.0,3
secdata, 1.2,1,0.0,3
secoffset,MID  
seccontrol,,,, , , ,
r=8500
percent=0.75      
dn=0.9
v=67     jichu=(percent*2-1)*r
jx=sqrt(r*r-jichu*jichu)
K, ,jx,,-jichu,
K, ,,,r,
k,,,,,
LARC,1,2,3,r,  
LDIV,1,dn, ,2,0
AROTAT,1,2, , , , ,3,2,180,4,lsel,s,,,1,9,2
LESIZE,all, , ,50, , , , ,1lsel,s,,,11,18,1
LESIZE,all, , ,10, , , , ,1lsel,s,,,2,10,2
LESIZE,all, , ,5, , , , ,1
 
amesh,allNSEL,S,LOC,Z,-jichu
d,all,,,,,,all,,,,,
NSEL,S,LOC,Y,0D,all, , , , , ,UY,ROTX,ROTZ, , , 
*GET,en,ELEM,,COUNT, , , , 
*do,i,1,en             esel,S,,,i   *GET,XE,ELEM,I,CENT,X*GET,YE,ELEM,I,CENT,Y          *GET,ZE,ELEM,I,CENT,Z  K,100, XE,YE,ZECSYS,2
*GET,zz,KP,100,LOC,y*GET,kyy,KP,100,LOC,Z
yy=(90-kyy)*3.1415926/180zz=zz*3.1415926/180q1=-1.3q2=0.1801*sin(yy)*cos(zz)q3=0.78*sin(yy)*sin(yy)q4=0.78*sin(yy)*sin(yy)*cos(2*zz)q5=0.25*0.5599*sin(yy)*sin(yy)*sin(yy)*cos(3*zz)q6=0.75*0.5599*sin(yy)*sin(yy)*sin(yy)*cos(zz)
f=v*v/16*2.125/100000*(q1+q2+q3+q4+q5+q6)sfe,i,1,pres,,-f,,,
*enddo
ACEL,0,0,9.8e-3,ALLSEL,ALL 
!* 
CSWPLA,11,2,1,1,
EMODIF,all,ESYS,11,FINISH 
/SOLANTYPE,0
 PSTRES,ON  SOLVE  
 
FINISH 
/SOLUTION  
ANTYPE,1
 
BUCOPT,LANB,3,0,0  /STATUS,SOLU
SOLVE