关于力矩施加问题？？求大侠帮忙解释下

keyking1013

/prep7
block,0,1,0,1,0,2
k,9,0.5,0.5,2.5
mp,ex,1,2e10
mp,prxy,1,0.2   
mp,prxy,1,0.3
r,2,1e-6
et,1,45
et,2,21
keyopt,2,3,0
lesize,all,0.2
vmesh,all
ksel,s,,,9
type,2
real,2
kmesh,all
allsel
nsel,s,loc,z,2,3
NPLOT   
CERIG,node(0.5,0.5,2.5),ALL,ALL, , , ,
allsel  
/SOLU
f,node(0.5,0.5,2.5),my,100e3
FINISH  
/SOL
nsel,s,loc,z,0  
d,all,all   
allsel  
solve
按上述命令出来的图是

keyking1013

按常理应该是第二幅图才对。。。求高手解释下。。。

ajacky007

估计是因为实体单元有三个自由度的原因吧

no_feat

第一主应力看起来还行，UZ的位移也还可以；就是等效应力理解不了，请高手指教一下。