ansys 模拟管材轧制,j结果出错，请指点？

davile123

本人 ansys 12.0 版本， 在网上找到一个ansys模拟管材轧制的例子（转帖他人成果，感谢原创者），按照相应的命令流做的时候，结果出错，请指点？不胜感激！以下是命令流！
*SET,r1,59.2
*SET,r2,280
*SET,r3,50  
*SET,ex,4   
*SET,alpha1,47.74   
*SET,alpha2,42.37
*SET,sp,10  
*SET,b,151.9
*SET,h,110.4
*SET,loc1x,-148.5
*SET,loc1y,-167.4
*SET,loc2x,0
*SET,loc2y,-4
*SET,loc3x,95.3067
*SET,loc3y,54.9922   
*SET,deltax1,5
*SET,circum,390.07
*SET,thin1,11.5
*SET,dia1,130
*SET,r1out,dia1/2
*SET,r1in,dia1/2-thin1
*SET,biteangle1,23.1155
*SET,dist,66.6611

/PREP7  
wpoff,0,loc2y,0
cyl4, , ,r1,alpha1,r1,90
csys,0  
wpave,0,0,0
csys,0  
wpoff,loc3x,loc3y,0
pcirc,r3, ,alpha2+180,270,  
csys,0  
wpave,0,0,0
asel,s,,,1,2,1
adele,all
lplot
lsel,s,,,2,3,1
lsel,a,,,5,6,1
ldele,all,,,1
allsel
k,,loc1x,loc1y
larc,1,4,3,r2
kdele,3
lplot
lsel,s,,,1,2,1
lsel,a,,,4
lcomb,all,,0
lplot

/pnum,kp,1
numcmp,kp
k,,loc3x+deltax1,sp/2  
k,,loc3x+deltax1,75
k,,loc3x+deltax1,225
k,,loc3x,225
k,,0,225
k,,646,75   
k,,651,80   
k,,651,225
l,7,6
l,6,2
l,2,3
l,3,4
l,4,5
l,5,6
l,4,8
l,8,9
l,9,10
l,10,5
lfillt,5,8,10
l,4,11

lsel,all
lsymm,x,all,,,,0,0

lsel,s,,,2
lsel,a,,,15
lglue,all

lsel,s,,,1
lsel,a,,,14
lglue,all


lsel,s,,,1,3,1
lsel,a,,,15,16,1
lsel,a,,,27
al,all

lsel,s,,,21,26,1
lsel,a,,,19
al,all

lsel,s,,,8,13,1
lsel,a,,,6
al,all

lsel,s,,,3,7,1
lsel,a,,,13
al,all

lsel,s,,,16,20,1
lsel,a,,,26
al,all

aplot

asel,s,,,1,5,1
vrotat,all,,,,,,10,24,360

vsel,s,,,1,6,5
vsel,a,,,11,16,5
vadd,all
vplot

vsel,s,,,3,8,5
vsel,a,,,13,18,5
vadd,all
vplot

vsel,s,,,2,7,5
vsel,a,,,12,17,5
vadd,all
vplot

vsel,s,,,5,10,5
vsel,a,,,15,20,5
vadd,all
vplot

vsel,s,,,4,9
vsel,a,,,14,19,5
vadd,all
vplot

vsel,s,,,2,5,3
vsel,a,,,21
vadd,all
vplot

numcmp,all

vsel,s,,,1,3,1
vsymm,y,all,,,,0,0


wpoff,,,-100
CYL4, , ,49.7, , , ,700
   
wpro,,,90.000000
wpro,,,90.000000
cone,49.7,20,0,100,0,360,  
wpro,,,-90.000000
wpro,,,-90.000000
wpave,0,0,0
csys,0

vadd,7,8


wpoff,0,0,dist+100
cyl4, , ,r1in, ,r1out, ,400
csys,0
wpave,0,0,0
csys,0
numcmp,all
allsel,all

et,1,solid164
keyopt,1,1,1
mp,ex,1,104.4e9
mp,dens,1,7870
mp,nuxy,1,0.4878
tb,plaw,1,,,1
tbdat,1,245e6
tbdat,2,35.08e9
tbdat,3,0
tbdat,4,40.0
tbdat,5,5.0
tbdat,6,0.75

mp,ex,2,170e9
mp,nuxy,2,0.25
mp,dens,2,7096
edmp,rigid,2,7,5

mp,ex,3,170e9
mp,nuxy,3,0.25
mp,dens,3,7096
edmp,rigid,3,7,5

mp,ex,4,ex,205e9
mp,nuxy,4,0.28
mp,dens,4,7830
edmp,rigid,4,4,7

type,1
mat,1
exto,vswe,auto,on
vsweep,7

type,1
mat,2
smrt,6
mshape,1,3d
mshkey,0
vmesh,1,2

smrt,5
mshape,1,3d
mshkey,0
vmesh,3

type,1
mat,3
smrt,6
mshape,1,3d
mshkey,0
vmesh,4,5

smrt,5
mshape,1,3d
mshkey,0
vmesh,6

type,1
mat,4
exto,vswe,auto,on
vsweep,5

numcmp,all
allsel,all

edpart,create

esel,s,mat,,1
nsle,s
cm,tube,node
allsel,all

  
edcgen,fnts, 1, 2,0.2,0,0,0,0, , , , ,0,10000000,0,0
edcgen,fnts, 1, 3,0.2,0,0,0,0, , , , ,0,10000000,0,0
edcgen,fnts, 1, 4,0.2,0,0,0,0, , , , ,0,10000000,0,0000,0,0

edcontact,0.1,0,2,1,4,1,1,4,0   
edsp,off,,,1

edve,velo,tube,0,0,-1000,0,0,0,,,,,,,

/solu

*dim,time,,2
time(1)=0,1
*dim,rotrolling,,2
rotrolling(1) = 69.5,79.05

*dim,velmandril,,2
velmandril(1,1,1) = -1600,-1610

wpoff,0,225,0  

edload,add,rbox,0,2,time,rotrolling,0,,,,,

wpoff,0,-450,0  

edload,add,rbox,0,3,time,rotrolling,0,,,,,

csys,0

wpave,0,0,0

csys,0

edload,add,rbvz,0,4,time,velmandril,0,,,,,

edenergy,1,1,0,0
time,1
edrst,100
edhtime,1000
eddump,1
solve
finish

davile123

???
怎么没有人回答啊？

davile123

???
期待回答啊！谢谢大牛了？

davile123

谁能回答，本人不甚感激！