编的APDL ，命令流，施加移动载荷，进行瞬态热分析

shiyan155165

我用APDL，*DO循环，和*IF判断语句，实现在长方体上施加往复移动载荷，进行瞬态热分析。


在X坐标，28,40.和970,982施加热流密度，大小12.
对流载荷，X 50---960施加对流载荷，体温度30，系数h130.
在X150---800施加往复移动热流密度，大小12.每三秒移动50mm。
循环，先设定为12次。


post1 的温度场图都显示不出来。post26的时间历程更别说了，读取不了.rth文件。

请求帮助，我这个是什么问题，万分感谢！！！
附命令流如下：
FINISH
/CLEAR
/FILENAME,Ball-screw
/TITLE,thermal-annalysis
!参数设定
*SET,a,0
*SET,j,0
*SET,tim,0
!设定分析类型
KEYW,PR_THERM,1
!进入前处理
/PREP7
ET,1,SOLID90
R,1,1
!定义材料导热系数
MP,KXX,1,53.6
MP,HF,1,130
BLC4,0,0,980,20,20
/REPLOT
!划分网格
SMRTSIZE,4
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL
/REPLOT
SAVE
FINISH


/SOLU
ANTYPE,4
TRNOPT,full
/PSF,PRES,NORM,2,0,1
!设置面载荷用箭头表示
!求解控制如下
time=0
TIMINT,ON
!施加热流密度载荷
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,28,40
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Y,20
NSEL,R,LOC,X,28,40
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,28,40
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,28,40
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,12

NSEL,S,LOC,Y,0
NSEL,R,LOC,X,965,977
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Y,20
NSEL,R,LOC,X,965,977
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,965,977
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,12
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,965,977
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,12


!施加对流载荷
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,50,960
NSEL,R,LOC,Z,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Y,20
NSEL,R,LOC,X,50,960
NSEL,R,LOC,Z,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,50,960
NSEL,R,LOC,Y,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,50,960
NSEL,R,LOC,Y,0,20
SF,ALL,CONV,130,30
!施加移动载荷
a=0
j=0
*DO,i,0,12,1
tim=tim+3
TIME,tim
!判断移动方向
*IF,a,EQ,0,THEN
j=j+1
*ELSE
j=j-1
*ENDIF
!施加螺母产生的热流密度
NSEL,S,LOC,Y,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,15
NSEL,S,LOC,Y,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SF,ALL,HFLUX,15
NSEL,S,LOC,Z,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,15
NSEL,S,LOC,Z,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SF,ALL,HFLUX,15
ALLS
NSUBST,5
OUTRES,ALL,ALL
EQSLV,JCG
!选择求解器
SOLVE
!删除上面施加的载荷
NSEL,S,LOC,Y,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SFDELE,ALL,HFLUX
NSEL,S,LOC,Y,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SFDELE,ALL,HFLUX
NSEL,S,LOC,Z,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SFDELE,ALL,HFLUX
NSEL,S,LOC,Z,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SFDELE,ALL,HFLUX
!施加这部分区域的对流载荷
NSEL,S,LOC,Y,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Y,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Z,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Z,0
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SF,ALL,CONV,130,30
NSEL,S,LOC,Z,20
NSEL,R,LOC,150+j*50,250+j*50
NSEL,R,LOC,Y,0,20
SF,ALL,CONV,130,30
ALLS
!判断是否转向
*IF,250+j*50,GE,850,THEN
j=12
a=1
*ELSEIF,150+j*50,LE,150,THEN
j=0
a=0
*ENDIF
*ENDDO
SAVE
FINISH



!进入POST1后处理器
/PSOT1
*SET,,,1,,50
PLNSOL,TEMP
/IMAGE,SAVE,wenduchang1,BMP
*SET,LAST
PLNSOL,SET
/IMAGE,SAVE,wenduchang2,BMP
ALLSEL




!进入POST26后处理器
/POST26
/AXLAB,X,TIME,(sec)
/AXLAB,Y,TEMP
!定义坐标轴标题
/GTHK,AXIS,3
/GTHK,CURVE,3
/COLOR,CURVE,MRED,1
!定义坐标轴粗度
WPCSYS,-1
/REPLOT
NSEL,S,LOC,X,150
NSEL,R,LOC,Y,10
!选择节点
*GET,NODE1,NODE,,NUM,MAX
!根据节点坐标读取最大节点编号
NSEL,S,LOC,X,325
NSEL,R,LOC,Y,10
*GET,NODE2,NODE,,NUM,MAX
NSEL,S,LOC,X,500
NSEL,R,LOC,Y,10
*GET,NODE3,NODE,,NUM,MAX
NSEL,S,LOC,X,675
NSEL,R,LOC,Y,10
*GET,NODE4,NODE,,NUM,MAX
NSEL,S,LOC,X,850
NSEL,R,LOC,Y,10
*GET,NODE5,NODE,,NUM,MAX
NSOL,2,NODE1,TEMP
!定义变量2,3,4,5,6
NSOL,3,NODE2,TEMP
NSOL,4,NODE3,TEMP
NSOL,5,NODE4,TEMP
NSOL,6,NODE5,TEMP
PLVAR,2
!绘制A（坐标？？）点温度随时间的变化规律曲线
/IMAGE,SAVE,bianhuaquxian1,BMP
!生成图文件，并保存
PLVAR,3
/IMAGE,SAVE,bianhuaquxian2,BMP
PLVAR,4
/IMAGE,SAVE,bianhuaquxian3,BMP
PLVAR,5
/IMAGE,SAVE,bianhuaquxian4,BMP
PLVAR,6
/IMAGE,SAVE,bianhuaquxian5,BMP
SAVE
FINISH
/EXIST,ALL

shiyan155165

真心求助！！！求大神知道

DONG.X.L

你好，我把你的程序稍微改了一下，但觉得还是有点不妥，有空交流下啊，发你站内信了，加QQ聊

夜魔狂舞

DONG.X.L 发表于 2012-9-26 14:47

                               
登录/注册后可看大图

你好，我把你的程序稍微改了一下，但觉得还是有点不妥，有空交流下啊，发你站内信了，加QQ聊

能不能给我也发一份啊大神，我也在学习，94197029@qq.com,谢谢。

lqchinavib

谢谢分享。。。。