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

Question

同一模型，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