Apply bending moment on a cylinder meshed with SHELL93

[Lirock]

Hi,
A cylinder with internal diameter d=0.09 mm, external diameter D=0.11 mm, and the length L=100mm.One end is fixed,and apply a bending moment 1mm.N on the other one.We can get the maximum deflection equals about 7.56 mm with formula.
I did it with ANSYS in two ways:
1. built a mass element at the center of the end, rigid the end, and applied the moment on the mass element.
2. average the moment and applied on all the nodes at the end.

The results shows that way NO.1 get a very close value with the formula one.So I think this way is good.
However, when I check the reaction solution, both the ways are not godd. I think the reaction moment should be equal with what I applied. But it is not. I wonder why. Is there anybody know about this.Thank you!

Rock

Here below you can find the commands I used:
/PREP7
ET,1,SHELL93
R,1,0.01, , , , , ,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E5
MPDATA,PRXY,1,,0.3
K, ,,,,
K, ,,,10,
PCIRC,0.1, ,0,360,
LSTR, 1, 2
FLST,2,4,4,ORDE,2
FITEM,2,1
FITEM,2,-4
ADRAG,P51X, , , , , , 5
MSHAPE,0,2D
MSHKEY,0
FLST,5,4,5,ORDE,2
FITEM,5,2
FITEM,5,-5
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
SMRT,6
FLST,5,4,5,ORDE,2
FITEM,5,2
FITEM,5,-5
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
ACLEAR,_Y1
AMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FLST,3,5,5,ORDE,2
FITEM,3,1
FITEM,3,-5
AGEN,2,P51X, , , ,-5, , ,0
ET,2,MASS21
R,2,1,1,1, , , ,
NKPT,0, 2
TYPE, 2
MAT, 1
REAL, 2
ESYS, 0
SECNUM,
TSHAP,LINE
E, 17825
FLST,5,33,1,ORDE,11
FITEM,5,5
FITEM,5,-8
FITEM,5,37
FITEM,5,-43
FITEM,5,782
FITEM,5,-788
FITEM,5,1158
FITEM,5,-1164
FITEM,5,1534
FITEM,5,-1540
FITEM,5,17825
NSEL,S, , ,P51X
FLST,5,32,1,ORDE,10
FITEM,5,5
FITEM,5,-8
FITEM,5,37
FITEM,5,-43
FITEM,5,782
FITEM,5,-788
FITEM,5,1158
FITEM,5,-1164
FITEM,5,1534
FITEM,5,-1540
/NOPR
CM,_NODECM,NODE
*SET,_z1, 17825
NSEL,S, , ,P51X
NSEL,A, , ,_Z1
CM,_CERGCM,NODE
CMSEL,S,_NODECM
/GO
/NOPR
CMSEL,S,_CERGCM
/GO
CERIG,17825,ALL,ALL, , , ,
/NOPR
CMSEL,S,_NODECM
/GO
ALLSEL,ALL
FLST,2,8,4,ORDE,4
FITEM,2,1
FITEM,2,-4
FITEM,2,14
FITEM,2,-17
/GO
DL,P51X, ,ALL,
SBCTRAN
FLST,2,1,1,ORDE,1
FITEM,2,17825
/GO
F,P51X,MX,1
FLST,2,32,1,ORDE,11
FITEM,2,8922
FITEM,2,9292
FITEM,2,-9299
FITEM,2,11520
FITEM,2,11890
FITEM,2,-11896
FITEM,2,13748
FITEM,2,14118
FITEM,2,-14124
FITEM,2,15975
FITEM,2,-15981
/GO
F,P51X,MX,1/32
ALLSEL
FINISH
/SOL
SOLVE