望哪位大侠不吝赐教，对1/4输油管端施加弯矩后，ANSYS不计算

qiuchengzhi0

小弟最近在做管道凹痕分析，当我在管道端部施加弯矩（用MPC184单元），在管道内侧施加内压（均布压力），在管道另一端施加凹痕（接触）。程序仅不计算我施加的弯矩荷载对模型的影响（应力）。望各位大小不吝赐教，此问题困扰小弟太久，若能解决，不胜感激！！！（这是小弟全部积蓄了，谢）
   以下是我做的部分命令流，我在接触部分暂不会命令，我用GUI路径做的。
      前处理：（除接触）
/prep7
et,1,solid95
mp,ex,1,2.1e5
mp,nuxy,1,0.3
TB,BKIN,1                  
TBDATA,1,485,769
cylind,480,500,0,1000,-90,90
sph4,0,540,0,40
vsbw,2
kwpave,9
wprota,0,90
vsbw,3
vsbw,4
wprota,0,0,90
vsbw,2
vsbw,3
vsbw,4
vsbw,5
vsbw,6
wpcsys,-1,0
wpoffs,0,0,300
vsbw,1
wprota,0,90
vsbw,6
vsbw,11
wprota,0,0,60
vsbw,12
vsbw,13
lesize,6,,,2
lesize,49,,,2
lesize,35,,,2
lesize,8,,,2
lesize,46,,,2
lesize,33,,,2
lesize,1,,,2
lesize,11,,,2
lesize,37,,,2
lesize,59,,,2
lesize,54,,,2
lesize,3,,,2
lesize,36,,,30
lesize,34,,,30
lesize,56,,,30
lesize,58,,,30
lesize,45,,,30
lesize,43,,,30
lesize,31,,,30
lesize,32,,,30
lesize,7,,,30
lesize,5,,,30
lesize,57,,,30
lesize,55,,,30
lesize,65,,,30
lesize,66,,,30
lesize,40,,,15
lesize,41,,,15
lesize,50,,,15
lesize,64,,,15
lesize,38,,,15
lesize,51,,,15
lesize,42,,,15
lesize,39,,,15
lesize,60,,,10
lesize,61,,,10
lesize,62,,,10
lesize,63,,,10
lesize,47,,,10
lesize,48,,,10
lesize,4,,,10
lesize,2,,,10
lesize,12,,,10
lesize,44,,,10
lesize,9,,,10
lesize,10,,,10
lesize,26,,,6
lesize,23,,,6
lesize,22,,,6
lesize,19,,,6
lesize,16,,,6
lesize,24,,,6
lesize,15,,,6
lesize,17,,,6
lesize,20,,,6
lesize,25,,,6
lesize,21,,,6
lesize,13,,,6
lesize,14,,,6
lesize,29,,,6
lesize,18,,,6
lesize,30,,,6
lesize,27,,,6
lesize,28,,,6
vmesh,all
et,2,184
keyopt,2,1,1
n,10000000,0,0,2000
nsel,s,loc,z,1000
type,2
*get,nnum,node,0,count
*get,ND,node,0,num,min
*do,i,2,nnum
E,10000000,ND
ND=NDNEXT(ND)
*enddo
allse
/solu
antype,0
nlgeom,on
autots,on
deltim,1e-4,1e-8,0.25,on
time,1
eqslv,sparse,,,
da,2,symm,,
da,46,symm,,
da,3,symm,,
da,29,symm,,
da,33,symm,,
da,13,symm,,
da,35,symm,,
da,50,ux,0,,
da,50,uy,0,,
da,50,uz,0,,
da,37,ux,0,,
da,37,uy,0,,
da,37,uz,0,,
da,1,ux,0,,
da,1,uy,0,,
da,1,uz,0,,
dk,9,ux,0
dk,9,uz,0
dk,10,ux,0
dk,10,uz,0
dk,14,uy,-10
dk,14,ux,0
dk,14,uz,0
sfa,32,,pres,4
sfa,52,,pres,4
sfa,48,,pres,4
sfa,39,,pres,4
sfa,5,,pres,4
sfa,30,,pres,4

f,10000000,mx,1e9

qiuchengzhi0

这是我做的结果，ANSYS没有计算由弯矩产生的应力