*DIM,WLY,ARRAY,361,1,1, , , !定义用于储存相位角的一维数组
*DIM,WLFY,ARRAY,361,1,1, , , !定义用于储存波流耦合力值的一维数组
*DO,I,0,360,1 !定义相位角范围0°~360°,并开始循环
*SET,WLY(I+1),I !将相位角值储存于WLX 数组
/PREP7 !进入前处理模块
TBDATA,81,I !将相位角取值赋予Water Table
/SOLU !进入求解计算模块
TIME,1 !定义结束结束时间
SOLVE !开始求解
*GET,RFY,NODE,1,RF,FY !提取该节点X 方向的反力
*SET,WLFY(I+1),-(RFY) !将所提取的反力取相反数后储存于WLFX 数组
*ENDDO !结束循环
*CREATE,ANSUO !创建一个文件
*CFOPEN,'RFORCE','OUT','' !打开文件
*VWRITE,WLY(1),WLFY(1), , , , , , , , !写入文件
(3F15.4) !
*CFCLOS
*END
/INPUT, ANSUO
这是我在某书中看到的程序命令,如果有同问的可以参考下,用了,应为我模型更复杂些,一维数组可能不够,就这么算的时候很多相位角算出来的值是零,不知道什么原因。 |