FINISH/CLEAR/FILNAME,DRILLINGRISER/TITLE.ANAKYSISOFDRILLINGRISER/UNIT,SI/PREP7ET,1,PIPE20ET,2,PIPE59R,1,0.5334,0.0254R,2,0.5334,0.0254,0.7,2.0,1030,0R,3,0.2534,0.0254MP,EX,1,2.1E11MP,PRXY,1,0.3MP,DENS,1,7850TBDE,WATE,1!设置波流参数TB,WATE,1,,,0!WATERTABLETBMODIF,1,1,2!选择STOKES五阶波理论TBMODIF,1,2,0!波流耦合作用方式TBMODIF,1,3,30!定义水深TBMODIF,1,4,1030!海水密度TBMODIF,1,5,0!波浪作用方向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,14,1,8!波高TBMODIF,14,2,7.4!波周期TBMODIF,14,3,57!波浪相位角K,1,0,0,18K,2,0,0,K,3,0,0,-30K,4,0,0,-33.2L,1,2L,2,3L,3,4LSEL,S,LINE,,1,2,1LATT,1,2,2LESIZE,ALL,1LMESH,ALLLSEL,S,LINE,,3LATT,1,1,1LESIZE,3,,,5LMESH,ALLALLSELSAVEFINISH/SOLUANTYPE,0NSEL,S,LOC,Z,-33.2D,ALL,ALLALLSELNSEL,S,LOC,Z,18D,ALL,UXD,ALL,UYALLSELNSEL,S,LOC,Z,1,17F,ALL,FX,6675/17ALLSELNSEL,S,LOC,Z,18F,ALL,FZ,-30E4ALLSELACEL,0,0,9,8SOLVESAVEFINISH/POST1SET,LAST!读入计算结果文件/ESHAPE,1!打开模型线形状选项PLNSOL,U,SUM,0,1!绘制结构整体变形图PLDISP,0!绘制结构位移等值图PLNSOL,S,EQV,0,1!绘制结构等效应力分布云图PLNSOL,EPTO,EQV,0,1!绘制结构等效应变分布云图ETABLE,MOMENT_I,SMISC,5!定义I端节点弯矩ETABLE.MOMENT_J,SMISC,11!定义J端节点弯矩PLLS,MOMENT_I,MOMENT_J,1,0!绘制弯矩图*DIM,WLX,ARRAY,361,1,1,,,!定义用于储存相位角的一维数组*DIM,WLFX,ARRAY,361,1,1,,,!定义用于储存波流耦合力值的一维数组*DO,I,0,360,1!定义相位角范围0°~360°,并开始循环*SET,WLX(I+1),I!将相位角值储存于WLX数组/PREP7!进入前处理模块TBDATA,81,I!将相位角取值赋予WATERTABLE/SOLU!进入求解计算模块TIME,1!定义结束结束时间SOLVE!开始求解*GET,RFX,NODE,50,RF,FX!提取该节点X方向的反力*SET,WLFX(I+1),-(RFX)!将所提取的反力取相反数后储存于WLFX数组*ENDDO!结束循环*CREATE,ANSUO!创建一个文件*CFOPEN,'RFORCE','OUT',''!打开文件*VWRITE,WLX(1),WLFX(1),!写入文件(3F15.4)!小数点后保留3位小数,整个输出占15个字符位置。*CFCLOS*END/INPUT,ANSUO!打开目录下此文件,观察相位角