哪位大侠帮我看看这个命令流错错错在哪里

杨叮当

哪位大侠帮我看看这个命令流错错错在哪里 例子是ANSYS高级工程有限元范例精选上 第十六章例子的命令流
/clear
/filenam,ACTUATOR
/PREP7
/GRAPHICS,FULL
!PARAMETER
*AFUN,DEG  !使用角度
ch=1/1000  !用于将MM转化为M
BTOL,1E-4*CH  !定义逻辑操作的公差
VOL=0  !低电位势0V
VOH=20  !高电位势20V

!#########Material  property#################
!3---PZT
et,3,solid5 !UX,UY, UZ,TEMP,VOLT,MAG
mp,nuxy,3,0.27   !?????
mp,dens,3,7.64e3 !kg/m3
mp,perx,3,8.93e-9 !F/m
mp,pery,3.8.93e-9  !F/m
mp,perz,3,6.92e-9  !F/m
mp,kxx,3,1
mp,murx,3,1

tb,anel,3   !Anisotropic  elastic matrix 72a??? 99.3.4
tbdata,1,11.565e10,6.489e10,6.229e10,0,0,0 !定义压电陶瓷的弹性矩阵
tbdata,7,11.565e10,6.229e10,0,0,0
tbdata,12,9.298e10,0,0,0
tbdata,16,1.786e10,0,0
tbdata,19,1.786e10,0
tbdata,21,1.786e10

tb,piez,3  !Piezoelectric matrix  72A??? 99.3.4
tbdata,1,0,0,-12.3082  !定义压电陶瓷压电矩阵
tbdata,4,0,0,-12.3082
tbdata,7,0,0,20.7574
tbdata,10,0,0,0
tbdata,13,0,17.0357,0
tbdata,16,17.0357,0,0
!下面定义其它材料，为看模型和选择方便，为每种材料定义单元，所以很多单元号都是
!*************4----potting adhesive***************
et,4,solid45
mp,ex,4,0.03448E10
mp,nuxy,4,0.3
mp,dens,4,1.1e3  !kg/m3

!******5---stainless steel********
et,5,solid45
mp,ex,5,19.3e10
mp,nuxy,5,0.29
mp,dens,5,8.03e3  !kg/m3

!**8---SLIDER****
et,8,solid45
mp,ex,8,38e10
mp,nuxy,8,0.29
mp,dens,8,4.086e3

!****13---Zirconia******
et,13,solid45
mp,ex,13,206E9
mp,nuxy,13,0.30
mp,dens,13,6.0e3  !kg/m3

!***15---Pt****
et,15,solid45
mp,ex,15,172E9
mp,nuxy,15,0.30
mp,dens,15,21.4e3  !kg/m3

!*****16---GOLD****
!et,16,solid45   
!mp,ex,16,75E9
!mp,nuxy,16,0.42
!mp,dens,16,19.3e3  !kg/m3

!-----AREA ELEMENT------
et,11,43
mp,ex,11,2.4e10
mp,nuxy,11,0.29
mp,dens,11,2.2e3  !kg/m3
R,11,0.001

!-------end--------
!下面定义几何模型参数和建立几何模型
!------parameter-------
W=1.29*CH
L0=1.9*CH
T_B=0.05*CH
T_H=0.05*CH
G=1.05*CH
L_H=0.3*CH
L_B=1.5*CH
T_PT=0.003*CH
T_PZT=0.015*CH
L_START=0.531*CH
L_END=1.7*CH
!----------Zirconia--------------
NUMSTR,KP,1   !为各部件独立，定义起点号
NUMSTR,LINE,1
NUMSTR,AREA,1
LOCAL,11,0, , , , ,90
CLOCAL,12,0,0.6*ch  !切换坐标系
CSYS,12
WPCSYS, ,12
K, ,-L_B,G/2+T_H-0.05*CH  !定义驱动器基体关键点
K, ,-L_B,G/2+T_H
K, ,-L_B+0.1*CH,G/2+T_H
!K, ,-L_B+0.15*CH,G/2+T_H   !?????????
K, ,-L_B+0.2*CH,G/2+T_H
K, ,-L_B+0.5*CH,G/2+T_H
K, ,-L_H,G/2+T_H
K, ,-L_H,G/2
K, , ,G/2
K, , ,G/2+T_H+T_B
K, ,-L_START,G/2+T_H+T_B
K, ,-L_END,G/2+T_H+T_B
K, ,-L0,G/2+T_H+T_B
KSYMM,2,1,12   !对称
FLST,2,24,3    !通过点形成面
FITEM,2,1
FITEM,2,2
FITEM,2,3
FITEM,2,4
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,8
FITEM,2,9
FITEM,2,10
FITEM,2,11
FITEM,2,12
FITEM,2,24
FITEM,2,23
FITEM,2,22
FITEM,2,21
FITEM,2,20
FITEM,2,19
FITEM,2,18
FITEM,2,17
FITEM,2,16
FITEM,2,15
FITEM,2,14
FITEM,2,13
A,P51X
!---------PZT_1----------
NUMSTR,KP,101    !为各部件独立，定义起点号
NUMSTR,LINE,101
NUMSTR,AREA,101
k, ,0.0,G/2+T_H+T_B+T_PT   !定义驱动器驱动部件一陶瓷关键点
k, ,-L_START,ky(101)
k, ,-L_END,ky(101)
k, ,-L0,ky(101)
k, ,-L_START-0.17*CH,G/2+T_H+T_B+4*(T_PZT+T_PT)
k, ,-L_END+0.1*CH,ky(105)
ksymm,y,101,106
a,9,101,102,103,104,12,11,10
a,102,105,106,103
a,21,107,108,109,110,24,23,22
a,108,109,112,111
!-----------PZT_2---------
NUMSTR,KP,201
NUMSTR,LINE,201
NUMSTR,AREA,201
k, , ,G/2+T_H+T_B+T_PT+T_PZT
k, ,-L0,KY(201)
k, , ,KY(201)+T_PT
k, ,-L0,KY(201)+T_PT
k, , ,KY(201)+T_PT+T_PZT
k, ,-L0,KY(201)+T_PT+T_PZT
k, , ,KY(201)+T_PT*2+T_PZT
K, ,-L0,KY(201)+T_PT*2+T_PZT
K, , ,KY(201)+T_PT*2+T_PZT*2
k, ,-L0,KY(201)+T_PT*2+T_PZT*2
k, , ,KY(201)+T_PT*3+T_PZT*2
k, ,-L0,KY(201)+T_PT*3+T_PZT*2
l,201,202
l,203,204
l,205,206
l,207,208
l,209,210
l,211,212
lsymm,y,201,206
FLST,3,6,4,ORDE,2
FITEM,3,201
FITEM,3,-206
ASBL,   102,P51X
FLST,3,6,4,ORDE,2
FITEM,3,207
FITEM,3,-212
ASBL,104,P51X
!--------PZT_3-----------
NUMSTR,KP,301
NUMSTR,LINE,301
NUMSTR,AREA,301
k, ,kx(104),ky(226)
ksymm,y,301
LARC,301,226,1,0.3*ch
LARC,302,202,4,0.3*ch
a,103,104,301,226
a,109,110,302,202
!---------slider---------
NUMSTR,KP,401
NUMSTR,LINE,401
NUMSTR,AREA,401
CSYS,11
WPCSYS,11
RECTNG,-1.235/2*CH,1.235/2*CH,-0.5*CH,0.5*CH  !工作体几何模型
!------------potting adhesine and dimple---------
NUMSTR,KP,501
NUMSTR,LINE,501
NUMSTR,AREA,501
CSYS,12
WPCSYS,12
W_GLUE=0.005  !建立工作体和驱动器间的胶
K, ,-W_GLUE*CH,2*CH
K, ,-W_GLUE*CH,-2*CH
K, ,-0.3*CH+W_GLUE*CH,2*CH
K, ,-0.3*CH+W_GLUE*CH,-2*CH
L,501,502
L,503,504
FLST,2,2,5,ORDE,2
FITEM,2,1
FITEM,2,401
FLST,3,2,4,ORDE,2
FITEM,3,501
FITEM,3,-502
ASBL,P51X,P51X
k, ,2*ch
k, ,
k, ,-2*ch
l,501,502
l,502,503
FLST,3,2,4,ORDE,2
FITEM,3,501
FITEM,3,-502
ASBL,504,P51X
FLST,2,4,3
FITEM,2,512
FITEM,2,510
FITEM,2,505
FITEM,2,513
A,P51X
FLST,2,4,3
FITEM,2,515
FITEM,2,507
FITEM,2,509
FITEM,2,511
A,P51X
!--------potting ashesive_2----------
NUMSTR,KP,601
NUMSTR,LINE,601
NUMSTR,AREA,601
rectng,-L0+0.075*ch,-L_B-0.075*ch,G/2+T_H,-(G/2+T_H)
FLST,2,2,5,ORDE,2
FITEM,2,508
FITEM,2,601
AOVLAP,P51X
!---------modify  elsctrode---------
NUMSTR,KP,701
NUMSTR,LINE,701
NUMSTR,AREA,701
k, ,-L_H,-1.0*ch   !建立电极几何模型
k, ,-L_H,1.0*ch
k, ,-L_B,-1.0*ch
k, ,-L_B,1.0*ch
l,701,702
l,205,229
l,105,111
l,703,704
FLST,2,2,5,ORDE,2
FITEM,2,101
FITEM,2,103
ASBL,P51X,   701
FLST,2,6,5,ORDE,6
FITEM,2,201
FITEM,2,203
FITEM,2,208
FITEM,2,210
FITEM,2,703
FITEM,2,-704
ASBL,P51X,   702
FLST,2,16,5,ORDE,9
FITEM,2,202
FITEM,2,204
FITEM,2,-207
FITEM,2,209
FITEM,2,211
FITEM,2,-214
FITEM,2,708
FITEM,2,710
FITEM,2,-714
FLST,3,2,4,ORDE,2
FITEN,3,703
FITEM,3,-704
ASBL,P51X,P51X
L,    229,    705
L,    205,    707
FLST,2,4,4
FITEM,2,219
FITEM,2,213
FITEM,2,712
FITEM,2,703
AL,P51X
FLST,2,4,4
FITEM,2,233
FITEM,2,201
FITEM,2,714
FITEM,2,704
AL,P51X
FLST,2,4,5,ORDE,4
FITEM,2,705
FITEM,2,708
FITEM,2,-709
FITEM,2,715
AADD,P51X
FLST,2,4,5,ORDE,4
FITEM,2,706
FITEM,2,-707
FITEM,2,710
FITEM,2,716
AADD,P51X
alls
!下面做平面网络划分，为体网络划分做准备，一旦体网络划分完毕，删去面单元。
!############# AREA  MESH  ##########
!CSYS,0
!WPCSYS,,0
MRATIO=0.50 !该参数用于调整网络密度
!------------2D MESH----------
TYPE,11
REAL,11
MAT,11
FLST,5,46,5,ORDE,4
FITEM,5,703
FITEM,5,-704
FITEM,5,717
FITEM,5,-760
ASEL,S, , ,P51X
mshkey,1
ESIZE,0.05*CH*MRATIO   !不同部分采用不同网络密度
AMESH,ALL
FLST,5,15,5,ORDE,14
FITEM,5,301
FITEM,5,-302
FITEM,5,501
FITEM,5,-502
FITEM,5,504
FITEM,5,-506
FITEM,5,511
FITEM,5,602
FITEM,5,701
FITEM,5,-702
FITEM,5,705
FITEM,5,711
FITEM,5,761
FITEM,5,-762
ASEL,S, , ,P51X
ESIZE,0.05*CH*MRATIO  !不同部分采用不同网络密度
mshkey,0
amesh,all
FLST,5,5,5,ORDE,5
FITEM,5,503
FITEM,5,507
FITEM,5,509
FITEM,5,-510
FITEM,5,601
ASEL,S, , ,P51X
ESIZE,0.10*CH*MRATIO   !不同部分采用不同网络密度
amesh,all
ALLS

!下面做体单元
!-----------SLIDER-----------
NUMSTR,KP,1001
NUMSTR,LINE,1001
NUMSTR,AREA,1001
NUMSTR,VOLU,1001     !为新产生的几何元素从一个新号开始
type,8
real,8
mat,8
esize, ,1
VEXT,ALL, , , , ,-0.25*CH  !拉伸出体单元
NUMSTR,KP,1200
NUMSTR,LINE,1200
NUMSTR,AREA,1200
NUMSTR,VOLU,1200
type,8
real,8
mat,8
esize,,1
FLST,5,4,5,ORDE,4
FITEM,5,503
FITEM,5,507
FITEM,5,509
FITEM,5,-510
ASEL,S, , ,P51X
VEXT,ALL, , , , ,(0.3-0.25)*CH  !工作体体单元
!删去面单元
!-------------CLEAR  2D-------------
ASEL,S,area, ,1,999
ACLEAR,ALL

!-----------MODIFY element------------
!-----------Zirconia-----------
FLST,5,6,6,ORDE,6
FITEM,5,1003
FITEM,5,-1004
FITEM,5,1007
FITEM,5,-1008
FITEM,5,1013
FITEM,5,-1014
VSEL,S, , ,P51X
ESLV
EMODIF,ALL,TYPE,13   !将第一层拉出的体单元属于驱动器基体的换单元，材料和常数
EMODIF,ALL,REAL,13
EMODIF,ALL,MAT,13
!-----------PZT-----------
FLST,5,26,6,ORDE,23
FITEM,5,1001
FITEM,5,-1002
FITEM,5,1017
FITEM,5,-1018
FITEM,5,1021
FITEM,5,-1024
FITEM,5,1026
FITEM,5,-1027
FITEM,5,1030
FITEM,5,-1031
FITEM,5,1035
FITEM,5,-1036
FITEM,5,1039
FITEM,5,-1040
FITEM,5,1045
FITEM,5,-1046
FITEM,5,1049
FITEM,5,-1051
FITEM,5,1053
FITEM,5,1055
FITEM,5,1057
FITEM,5,1059
FITEM,5,-1060
VSEL,S, , ,P51X
ESLV
EMODIF,ALL,TYPE,3     !换成驱动器驱动部件一陶瓷
EMODIF,ALL,REAL,3
EMODIF,ALL,MAT,3
!--------------Electrode--------
FLST,5,28,6,ORDE,21
FITEM,5,1015
FITEM,5,-1016
FITEM,5,1019
FITEM,5,-1020
FITEM,5.1025
FITEM,5,1028
FITEM,5,-1029
FITEM,5,1032
FITEM,5,-1034
FITEM,5,1037
FITEM,5,-1038
FITEM,5,1041
FITEM,5,-1044
FITEM,5,1047
FITEM,5,-1048
FITEM,5,1052
FITEM,5,1054
FITEM,5,1056
FITEM,5,1058
FITEM,5,1061
FITEM,5,-1066
VSEL,S, , ,P51X
ESLV
EMODIF,ALL,TYPE,15   !换成电极
EMODIF,ALL,REAL,15
EMODIF,ALL,MAT,15
!------------potting adhesive---------
FLST,5,2,6,ORDE,2
FITEM,5,1006
FITEM,5,1012
VSEL,S, , ,P51X
ESLV
EMODIF,ALL,TYPE,4   !换成相连的胶
EMODIF,ALL,REAL,4
EMODIF,ALL,MAT,4
ALLS
!----------------------------------
csys,0
wpcsys, ,0
!读取后面要使用的节点号
!----------GET THE POINT ON SLIDER---------
KSEL,S,KP, ,1209
NSLK,S
*GET,N01,NODE,,NUM,MIN, , , ,
KSEL,S,KP, ,1207
NSLK,S
*GET,N02,NODE,,NUM,MIN, , , ,
KSEL,S,KP, ,1200
NSLK,S
*GET,N03,NODE,,NUM,MIN, , , ,
KSEL,S,KP, ,1202
NSLK,S
*GET,N04,NODE, ,NUM,MIN, , , ,
ALLS
!定义边界条件
!-------------boundary condition--------------
ASEL,S,AREA, ,1052
NSLA, ,1
D,ALL,UX,0
D,ALL,UY,0
D,ALL,UZ,0
NSEL,S,NODE,,N01
NSEL,A,NODE,,N02
NSEL,A,NODE,,N03
NSEL,A,NODE,,N04
!D,ALL,UY,0
ALLS
!耦合要加电位势的面上的节点，并读取最小节点
!--------------COUPLE THE TWO PACE OF ELECTRODE
FLST,5,4,5,ORDE,4
FITEM,5,1219
FITEM,5,1225
FITEM,5,1230
FITEM,5,1244
ASEL,S, , ,P51X
NSLA, ,1
CM,UPTOP,NODE
CP,1,VOLT,ALL
*get,UP_TOP,node,0 ,num,min   !将加相同电位势的面耦合，获最小节点号

FLST,5,4,5,ORDE,4
FITEM,5,1220
FITEM,5,1226
FITEM,5,1228
FITEM,5,1245
ASEL,S, , ,P51X
NSLA, ,1
CM,UPBOT,NODE
CP,2,VOLT,ALL
*get,UP_BOT,node,0,num,min

FLST,5,4,5,ORDE,4
FITEM,5,1223
FITEM,5,1234
FITEM,5,1236
FITEM,5,1241
ASEL,S, , ,P51X
NSLA, ,1
CM,DNTOP,NODE
CP,3,VOLE,ALL
*get,DN_TOP,node,0,num,min!将加相同电位势的面耦合，获最小节点号

FLST,5,4,5,ORDE,4
FITEM,5,1222
FITEM,5,1233
FITEM,5,1238
FITEM,5,1242
ASEL,S, , ,P51X
NSLA, ,1
CM,DNBOT,NODE
CP,4,VOLT,ALL
*get,DN_BOT,node,0,num,min    !将加相同电位势面耦合，获取最小节点号
ALLS
finish

!------------求解----------------
/solu
antype,static   !静力求解
CSYS,0
ALLS
d,UP_TOP,volt,voh
d,UP_BOT,volt,vol
d,DN_TOP,volt,-voh
d,DN_BOT,volt,vol
d,edge_0,volt,vol

ESEL,S,TYPE,,3
NSLE
D,ALL,TEMP
D,ALL,MAG
ALLS

time,1
!nsubst,10  !10个字部
outres,all,all   !输出所有子步的结果
kbs,0
nsel,all
esel,all
save
fini
/POST1
KSEL,S,,403     !读工作体最前端摆动
NSLK
*GET,NSLD,NODE,,NUM,MIN
STROKE=UZ(NSLD)   !将摆动用参数记录