菜鸟问个问题，模型变形一点也不明显，请大家帮忙

hashfinger

/PREP7
!!!!!点
k,1,0,4.720,0.525
k,2,0,4.720,0.265
k,3,0,4.970,0.215
k,4,0,4.970,0.05
k,5,0,4.800,0.05
k,6,0,4.800,0.1
k,7,0,4.720,0.1
k,8,0,4.720,0
k,9,0,5.120,0
k,10,0,5.120,0.375
k,11,0,4.800,0.439
k,12,0,4.800,0.525

!!!!截面
a,1,2,3,4,5,6,7,8,9,10,11,12
!!!!!轴
k,100,0,0,10
k,101,0,0,-10

!!!!!体
vrotat,1,,,,,,100,101,180

!!!!!划分网格
csys,1
lsel,s,loc,y,-45 !从全部线中选中一组线，y坐标-45
lsel,r,loc,x,5.12 !从选中线中选一组线，x坐标5.12
lesize,all,,,50 !所有选中线网格分成50份
allsel !选中所有项
nummrg,all !合并所有项
et,1,solid90,,,
mp,dens,1,7800 !定义材料号及特性，项目dens（质量密度），材料号1，值7800
mp,ex,1,210
mp,kxx,1,70
mp,prxy,1,0.3

!!!!!网格
SMRT,8
MSHAPE,1,3D !指定单元形状，三角形，3D
vmesh,all

!!!!!!对称面约束
/solu
csys,1
asel,s,loc,y,90 !选中Y坐标值90的面
asel,a,loc,y,-90 !加选
da,all,symm,
allsel
asel,s,loc,z,0
da,all,symm !载荷传递
allsel
　
!!!!!!支承约束　
csys,1　　
nsel,s,loc,x,r1
nsel,r,loc,y,-60,-59
!!
finish
/prep7
nrotat,all !所有节点旋转至激活坐标系
finish
/solu
d,all, ,0, , , ,UX,,, , ,
!!
allsel

!!!!!重力
csys,0
acel,,9.8

　 !!!!!内圈载荷
csys,0
Q=4000000/0.75 !Q为每一档轴向单位长度上受的力
pi=3.1416
a=30*pi/180　　　 !a为30度的弧度
W=Q/(4.92*(pi-a+0.5*sin(2*a))) !此公式为《回转窑》上加力公式，4.92为滚圈中径
csys,1　　　　　　 !进入柱座标
nsel,all !选择一组新节点为下面坐准备
nsel,s,loc,x,4.72 !内径
nsel,u,loc,y,60,90 !当前组中不选的部分坐标Y，范围60－90

*get,nmax,node,,num,max
*get,nmin,node,,num,min
*afun,deg
*dim,t1,array,nmax,1,1 !定义数组，名t1,性质array
csys,1
*do,k,nmin,nmax !定义一个循环开始，变量k,起始值nmin，终值nmax
*if,nsel(k),eq,1,then
b=cos(ny(k)-90)
t1(k)=W*(cos(30)-b) !函数加力
*else
t1(k)=0
*endif
*enddo
sffun,pres,t1(1)
sf,all,pres,0
allsel

!!!!!求解
solve
finish

是一个应力分析的问题。
我用上面的命令流作出来后，发现模型变形一点也不明显，几乎没有。如果换成SOLID95的话，又会出现错误提示：The value of UZ at node 1015 is 2262860.95. It is greater than the
current limit of 1000000. This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model is properly constrained.

大家能帮我看看么，谢谢啊

[ 本帖最后由 mjhzhjg 于 2006-12-29 13:12 编辑 ]