该怎么得到壳单元的剪力和弯矩？

梁单元求解后可以通过建立单元表smisc，8来得到梁内的剪力。我尝试用shell63建立的模型截面是1*1矩形，长10，加均匀载荷，最后得到的最大应力和手工计算一致，想和beam单元建模一样得到单元的剪力、弯矩，不知道该怎么做，请各位高手赐教！谢谢！
命令流如下：
/PREP7
ET,1,SHELL63
R,1,1, , , , , ,
RMORE, , , ,
RMORE
RMORE, ,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,8e9
MPDATA,PRXY,1,,0.33
RECTNG,0,10,0,1,
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
FLST,5,2,4,ORDE,2
FITEM,5,2
FITEM,5,4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,1, , , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,1
FITEM,5,3
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,50, , , , ,1
!*
MSHAPE,0,2D
MSHKEY,1
!*
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FINISH
/SOL
FLST,2,2,4,ORDE,2
FITEM,2,2
FITEM,2,4
!*
/GO
DL,P51X, ,UX,
FLST,2,2,4,ORDE,2
FITEM,2,2
FITEM,2,4
!*
/GO
DL,P51X, ,UY,
FLST,2,2,4,ORDE,2
FITEM,2,2
FITEM,2,4
!*
/GO
DL,P51X, ,UZ,
FLST,2,1,5,ORDE,1
FITEM,2,1
/GO
!*
SFA,P51X,1,PRES,10
/STATUS,SOLU
SOLVE