麻烦帮我看看这个命令流错在哪里？谢谢！

xinercsm

各位大虾：

麻烦帮我看看这个命令流错在哪里？谢谢！

本命令流是高炳军著作《ANSYS在机械与化工装备中的应用》上的例子E322，为什么运行后显示显示网格错误，希望高手指教！谢谢！

附上：ansys命令流文件

finish
/clear
/filn, E322  
/title, FEA of connecting zone of nozzle to cylinder
! *************参数设定*******************
Rci=1000                 !筒体内半径
tc=30                         !筒体厚度
Rco=Rci+tc                      !筒体外半径
Lc=4000                         !筒体长度
Rno=530                         !接管外半径
tn=15                           !接管厚度
Rni=Rno-tn                 !接管内半径
Hmin=sqrt(Rci**2-Rni**2)        !接管最小内伸长度
Lni=50           !接管内伸长度
Ln=500                         !接管外伸长度
rr1=30                          !焊缝外侧过渡圆角半径
rr2=15                         !焊缝内侧过渡圆角半径
pi=1.2                         !内压
pc=pi*Rci**2/(Rco**2-Rci**2) !筒体端部轴向平衡面载荷
pn=pi*Rni**2/(Rno**2-Rni**2)    !接管端部纵向平衡面载荷
!****************前处理***************************
/prep7
et,1,45                         !定义单元类型
mp,ex,1,2e5                 !定义材料的弹性模量
mp,nuxy,1,0.3                 !定义材料的泊松比
!****************建立模型***************************
cylind,Rco,Rci,0,-Lc/2,90,270, !生成筒体
wpoff,0,0,-Lc/2                 !将工作面沿-Z向移动Lc/2
wprot,0,90,                  !将工作面沿yz旋转90o
cylind,Rno,Rni,-Ln-Rci-tc,-Hmin+Lni,90,180,      !生成接管
vovlap,all                 !体overlap布尔运算
vdel,3,,,1               !删除筒体上开孔失去的部分
vsel,s,,,7               !选择筒体
*afun,deg                 !设定角度函数中单位为角度
ang1=2*nint(asin(Rno/Rci)) !计算接管区切割角度
wprot,0,0,-90+ang1         !旋转坐标系
vsbw,all                 !切割筒体
afillt,21,12,rr1           !筒体与接管外表面圆角
al,33,60,5                 !生成圆角区域所需的面
al,37,61,7                 !生成圆角区域所需的面
va,12,14,16,34,35         !生成圆角区域体
afillt,23,13,rr2         !筒体内表面与接管外表面圆角
al,24,40,68                 !生成圆角区域所需的面
al,25,41,69                 !生成圆角区域所需的面
va,13,23,31,39,40         !生成圆角区域体
alls,all                 !全选
askin,62,70                 !根据内外圆角边界蒙皮生成切割面
vsba,2,36                 !切割筒体
wprot,0,,90-ang1         !旋转坐标系
wprot,0,90,0                 !旋转坐标系
wpoff,0,0,-2*Rni         !移动坐标系至接管区轴向切割位置
vsbw,9                         !切割筒体
vsbw,1                       !切割筒体
!*************划分网格******************************
lsel,s,loc,y,Ln+rci+tc         !选择接管端部各线段
lsel,u,length,,0,tc         !去掉长度为0至tc的线段
lesize,all,,,20,,1            !设定所选线段的剖分数
alls,all
lesize,60,,,4,,1         !设定外圆角某一线段的剖分数
lesize,45,,,12,,1               !设定内接管分段数
lesize,57,,,12,,1
asel,s,,,45,46,1,1         !选择筒体上接管连接影响区外侧两截面
accat,all                 !连接相邻面域
alls,all                        !全选
vsweep,all                      !映射划分
alls
nummrg,all               !合并所有相同项
numcmp,all                 !压缩
fini
!****************求解***************************
/solu
asel,s,loc,z,0                 !选择筒体端面
SFA,all,1,PRES,-pc         !施加端面平衡面载荷
asel,s,loc,x,0                  !选择对称面
asel,a,loc,z,-Lc/2         !选择对称面
DA,all,SYMM                 !施加对称约束
asel,s,loc,y,Ln+Rci+Tc         !选择接管端面
DA,all,UY,                  !约束轴向位移
asel,s,,,6                 !选择所有内表面
asel,a,,,7
asel,a,,,15
asel,a,,,17
asel,a,,,19
asel,a,,,23
asel,a,,,29
asel,a,,,42
asel,a,,,47
asel,a,,,50
SFA,all,1,PRES,pi         !施加内压
alls
solve                         !求解
fini
!****************后处理***************************
/post1
PLNSOL,S,INT,0,1  !显示应力云图
fini