24小时热门版块排行榜

返回列表

隐月sea

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 17.5
帖子: 3
在线: 5.8小时
虫号: 1841494
注册: 2012-05-31
性别: MM
专业: 机械动力学

[求助] ANSYS 循环移动热载荷结果不收敛求助

我做的是在丝杆上的循环移动热载荷  不知道哪里错了  结果不收敛  速度8000mm/min  r移动块250mm  热流密度6500W/m2 两端有轴承发热 37W/m2
移动距离是1.6m 求高人指点啊  其余部分的对流载荷 7W/M2,K  来回循环12次
下面是我做的命令流
FINISH
/CLEAR
/FILENAME,sea
/TITLE,thermal-structure annalysis
*SET,a,0
*SET,j,0
*SET,tim,0

/PREP7                   !进入前处理
ET,1,SOLID90
R,1,1
MP,KXX,1,48                !定义材料导热系数
Mp,dens,7900
Mp,c,1,460
BLC5,0,0,0.027,0.027,-0.02
cyl4,0,0,0.045/2,360,,,0.225
wpoff,0,0,0.225
cyl4,0,0,0.063/2,360,,,2.251
wpoff,0,0, 2.251
cyl4,0,0,0.025,360,,,0.215

/REPLOT
Wprota,,90
vsbw,all
Wprota,,,90
vsbw,all
nummrg,all
numcmp,all


csys,0                !划分网格
MSHAPE,0,3D
MSHKEY,1
Lsel,s,loc,z,0
Lsel,a,loc,z,-0.02
Lsel,a,loc,z,0.225
Lsel,a,loc,z,0.225+2.251
Lsel,a,loc,z,0.225+2.251+0.215
lesize,all,0.005
lsel,inve
lesize,all,0.02
Vsweep,all
/REPLOT
SAVE
FINISH

/SOLU
ANTYPE,4
TRNOPT,full
/PSF,PRES,NORM,2,0,1    !设置面载荷用箭头表示
                     !求解控制
nropt,auto,,
Eqslv,
Solcontrol,on
Autots,on
Kbc,1
Outres,all,all,
Time,0.001
deltim,0.001,0.001,0.001
TIMINT,off
Tunif,20

Csys,0
Nsel,s,loc,z, -0.02
SF,ALL,CONV,7,20
Nsel,s,loc,z, 0.225+2.251+0.215+0.005
SF,ALL,CONV,7,20

NSEL,S,LOC,x,-0.027/2,0.027/2          !施加对流载荷
NSEL,R,LOC,y,0.027/2
NSEL,R,LOC,z, -0.02,0
SF,ALL,CONV,7,20
NSEL,S,LOC,x,- 0.027/2, 0.027/2
NSEL,R,LOC,y,- 0.027/2
NSEL,R,LOC,z, -0.02,0
SF,ALL,CONV,7,20
NSEL,S,LOC,y,- 0.027/2, 0.027/2
NSEL,R,LOC,x, 0.027/2
NSEL,R,LOC,z, -0.02,0
SF,ALL,CONV,7,20
NSEL,S,LOC,y,- 0.027/2, 0.027/2
NSEL,R,LOC,x,- 0.027/2
NSEL,R,LOC,z, -0.02,0
SF,ALL,CONV,7,20

Csys,1
NSEL,S,LOC,x,0.045/2
NSEL,R,LOC,z,0,0.225-0.15
SF,ALL,CONV,7,20
NSEL,S,LOC,x,0.025
NSEL,R,LOC,z, 0.225+2.251+0.15, 0.225+2.251+0.215
SF,ALL,CONV,7,20

NSEL,S,LOC,x,0.063/2
NSEL,R,LOC,z, 0.225,0.225+2.251
SF,ALL,CONV,7,20
allsel,all
Solve

Timint,on             !施加移动载荷

*DO,i,0,12*12,1             ! tim=0          time,tim
Outres,all,all,
tim=tim+1
TIME,tim

Sfdele,all,hflux

Csys,1
NSEL,S,LOC,x,0.045/2          !施加热流密度载荷
NSEL,R,LOC,z,0.225-0.15,0.225
SF,ALL,HFLUX,37
NSEL,S,LOC,x,0.025
NSEL,R,LOC,z,0.225+2.251,0.225+2.251+0.15
SF,ALL,HFLUX,37

NSEL,S,LOC,x,0.063/2       !施加螺母产生的热流密度
NSEL,R,LOC, z,0.435+j*0.134,0.685+ j*0.134
SF,ALL,HFLUX,7540
ALLS
NSUBST,5
EQSLV,JCG                !选择求解器
SOLVE

NSEL,S,LOC,x,0.063/2       !删除上面添加的载荷
NSEL,R,LOC, z,0.435+ j*0.134,0.685+ j*0.134
SFDELE,ALL,HFLUX,7540

NSEL,S,LOC,x,0.063/2       !施加这部分的对流热载荷
NSEL,R,LOC, z,0.435+ j*0.134,0.685+j*0.134
SF,ALL,conv,7,20
alls

*IF, 0.435+j*0.134,GE,0.435+1.55,THEN             !判断是否转向
j=12
a=1
*ELSEIF, 0.435+j*0.134,Lt,0.435+0.134,THEN
j=0
a=0
*ENDIF

*IF,a,EQ,0,THEN          !判断移动方向
j=j+1
*ELSE
j=j-1
*ENDIF

*ENDDO
SAVE
FINISH

回复此楼