施加扭矩求扭转角，进而求扭转刚度的问题

lideshui

我根据venture和rodge两人所说的加扭矩的方法，试了第二种方法。数据都是和vrnture的一样。
模型简单取1个空心梁，采用实体建模，单元用95，内半径10mm，外半径15mm，长度200mm，矩10E6 Nmm，mass21单元实常数取零。
(1)
梁建模划网
(2)
在梁轴线上生成1节点，偏离端面10mm。
(3)
选定mass21单元属性，通过第二步的节点直接生成单元(E或EN)。
(4)
指定刚性化区域（CERIG），先选取mass21对应的节点，再选取梁端面节点。
(5)
加约束和载荷，梁另一端全约束，在mass21对应节点上加扭矩。

/PREP7
ET,1,95
ET,2,21
KEYOPT,2,3,0
R,1,1E-6
MP,EX,1,2.01e5
MP,PRXY,1,0.3

CYLIND,15,10,0,200,0,360,
wpro,,90,
vsbw,all
wpro,,,90
vsbw,all
WPCSYS,-1,0
K,17 , , ,210
lsel,s,,,13,16,1
lesize,all, , ,8, , , , ,1
lsel,s,,,22
lesize,all, , ,4, , , , ,1
lsel,s, , ,17,20,1
lsel,a, , ,26,27,1
lsel,a, , ,30,31,1
lesize,all, , ,20,0.4, , , ,1
alls
vmesh,all
！！！！！下面一段开始各个方法有所不同，由于前面的建模一样，后面的例子就不再给出
ksel,s,,,17
type,2
real,1
kmesh,all
allsel
nsel,s,loc,z,200,210
npolt
CERIG,node(0,0,210),ALL,ALL, , , ,
!!!!!CERIG命令定义局部刚性区域
allsel
/SOLU
f,node(0,0,210),mz,10e5
FINISH
！！！！！以下一段边界条件的施加各种方法一样，后面例子也不再赘述
/SOL
nsel,s,loc,z,0
d,all,all
allsel
solve



我想得到加扭矩后轴的扭转角，是不是在柱坐标系下看轴的圆周表面和加扭矩表面相交的线的节点的位移。这个位移在柱坐标系下到底是弧长还是弧度，或者是角度，如果是弧度或角度的话，所有圆周表面的节点位移图应该是相等的啊，但是我在柱坐标系下看到的位移图是这样的

这是怎么回事啊。怎样得到正确的扭转角，进而求出扭转刚度，希望各位高手指点下，小弟迷茫了很久。

bjjob1

好像没这个功能吧？