混凝土收缩模拟求助！

命令如下：/PREP7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!*
ET,1,SOLID65 !定义单元
!*
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MP,EX,1,4.3e4 !混凝土弹性模量
MP,PRXY,1,0.17
MP,ALPX,1,0.7e-5
MP,kxx,1,1.63
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
block,0,325,0,100,0,100 !建立长方块
wpoff,162.5,50
CYLIND,9, ,0,85,0,360, !建立混凝土中的栓钉约束
vsbv,1,2
numcmp,all
VSEL,S, , ,1 !选择1号体（混凝土）
VATT,1, , 1,0 !指定1体材料属性1（混凝土）
allsel
MSHAPE,1,3D
MSHKEY,0 !自由划分
VMESH, 1 !自由划分混凝土体
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
finish
/PREP7
asel,s,,,7,8
NSLA,S,1
FINISH
/PREP7
CSWPLA,11,1,1,1, !创建一个新的局部坐标系（柱坐标）
FLST,2,568,1,ORDE,4 !将局部坐标系面上的node，全转换为柱坐标
FITEM,2,1
FITEM,2,-24
FITEM,2,56
FITEM,2,-599
NROTAT,P51X
FINISH
/SOL
FLST,2,568,1,ORDE,4 !将这些柱坐标下的node的x方向约束，即形成了径向约束
FITEM,2,1
FITEM,2,-24
FITEM,2,56
FITEM,2,-599
!*
/GO
D,P51X, ,0, , , ,UX, , , , ,
csys,4
allsel
asel,s,,,9
nsla,s,1
d,all,uz,0
asel,s,,,6
nsla,s,1
d,all,uz,0
allsel
/SOLU
tref,40 !定义温度荷载
bfunif,temp,20
ALLSEL
SOLVE !求解
这个命令式模拟混凝土收缩的，在混凝土中有一个栓钉作为约束，我将栓钉处的混凝土径向进行了约束，并将混凝土底部UY也进行了约束。
荷载为温度从40降到20
但是结果却好像不对，按道理，混凝土纯收缩的话只会沿着XYZ轴均匀的产生位移。
但是分析下来， 整个混凝土却产生了旋转，导致结果不符合预期，请教各位是哪里出错了，谢谢！