ansys搜索波浪相位角的问题

流朦日月 · 发表于 2015-3-6 16:53

这是用ANSYS搜索波浪相位的命令流，请大神帮忙看看哪里有问题？老是报错。
FINISH
/CLEAR
/FILNAME,Drilling riser
/TITLE,Analysis of Drilling riser
/PREP7

ET,1,PIPE59
ET,2,PIPE20

R,1,533.4,25.4,
R,2,533.4,25.4,0.7,2,1.03e-6,0,

MP,EX,1,2.1e5
MP,PRXY,1,0.3
MP,DENS,1,7.85E-6

TB,WATE,1,,,0
TBMODIF,1,1,2
TBMODIF,1,2,0
TBMODIF,1,3,30
TBMODIF,1,4,1030
TBMODIF,1,5,0
TBMODIF,1,6,
TBMODIF,2,1,-30
TBMODIF,2,2,0.96
TBMODIF,2,3,0
TBMODIF,2,4,-15
TBMODIF,2,5,1.47
TBMODIF,2,6,0
TBMODIF,3,1,0
TBMODIF,3,2,1.93
TBMODIF,3,3,0
TBMODIF,3,4,
TBMODIF,3,5,
TBMODIF,3,6,
TBMODIF,4,1,
TBMODIF,4,2,
TBMODIF,4,3,
TBMODIF,4,4,
TBMODIF,4,5,
TBMODIF,4,6,
TBMODIF,5,1,
TBMODIF,5,2,
TBMODIF,5,3,
TBMODIF,5,4,
TBMODIF,5,5,
TBMODIF,5,6,
TBMODIF,6,1,
TBMODIF,6,2,
TBMODIF,6,3,
TBMODIF,6,4,
TBMODIF,6,5,
TBMODIF,6,6,
TBMODIF,7,1,
TBMODIF,7,2,
TBMODIF,7,3,
TBMODIF,7,4,
TBMODIF,7,5,
TBMODIF,7,6,
TBMODIF,8,1,
TBMODIF,8,2,
TBMODIF,8,3,
TBMODIF,8,4,
TBMODIF,8,5,
TBMODIF,8,6,
TBMODIF,9,1,
TBMODIF,9,2,
TBMODIF,9,3,
TBMODIF,9,4,
TBMODIF,9,5,
TBMODIF,9,6,
TBMODIF,10,1,
TBMODIF,10,2,
TBMODIF,10,3,
TBMODIF,10,4,
TBMODIF,10,5,
TBMODIF,10,6,
TBMODIF,11,1,
TBMODIF,11,2,
TBMODIF,11,3,
TBMODIF,11,4,
TBMODIF,11,5,
TBMODIF,11,6,
TBMODIF,12,1,
TBMODIF,12,2,
TBMODIF,12,3,
TBMODIF,12,4,
TBMODIF,12,5,
TBMODIF,12,6,
TBMODIF,13,1,
TBMODIF,13,2,
TBMODIF,13,3,
TBMODIF,13,4,
TBMODIF,13,5,
TBMODIF,13,6,
TBMODIF,14,1,8
TBMODIF,14,2,7.4

csys,0
K,1,,,18,
K,2,,,0,
K,3,,,-30,
K,4,,,-33.2,
lstr,1,2
lstr,2,3
lstr,3,4

LSEL,S,LINE,,1,2,1,,
LATT,1,2,1
LESIZE,ALL,1
LMESH,ALL

ALLS
LSEL,S,LINE,,3
LATT,1,1,2
LESIZE,3,,,5
LMESH,ALL

allsel
save
finish

/PREP7
D,50,ALL
ACEL,0,0,9.8
NPLO

*DIM,WLX,ARRAY,361,1,1,,,       !定义用于储存相位角的一维数组
*DIM,WLFY,ARRAY,361,1,1,,,    !定义用于储存波流耦合力值的一维数组

*DO,I,0,360,1                   !定义相位角范围0°~360°，并开始循环
*SET,WLX(I+1),I                !将相位角值储存于WLX 数组
/PREP7                         !进入前处理模块AFUN,DEG
TBDATA,81,I                   !将相位角取值赋予Water Table
/SOLU                         !进入求解计算模块
TIME,1                         !定义结束结束时间
SOLVE                         !开始求解
*GET,RFX,NODE,50,RF,FX          !提取该节点X 方向的反力
*SET,WLFY(I+1),-(RFX)          !将所提取的反力取相反数后储存于WLFX 数组
*ENDDO                         !结束循环
*CREATE,ANSUO                   !创建一个文件
*CFOPEN,RFORCE,OUT,             !打开文件
*VWRITE,WLX(1),WLFY(1),,,,,,,, !写入文件
(3F15.4)
*CFCLOS
*END
/INPUT, ANSUO

[Ansys] ansys搜索波浪相位角的问题