24小时热门版块排行榜

>论坛更新日志 (3725)
>虫友互识 (614)
>文献求助 (214)
>休闲灌水 (213)
>导师招生 (198)
>考博 (129)
>论文投稿 (118)
>博后之家 (90)
>硕博家园 (81)
>招聘信息布告栏 (71)
>论文道贺祈福 (68)
>基金申请 (57)
>教师之家 (49)
>考研 (47)
>找工作 (45)
>公派出国 (21)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yunshan216

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 1300
红花: 4
沙发: 1
帖子: 513
在线: 34.7小时
虫号: 8588912
注册: 2018-04-22
专业: 机械摩擦学与表面技术

[求助] 热分析不收敛，出现负温度已有2人参与

在热分析时不收敛，温度还会比常温还低，不知道怎么回事。热原模型是变半径的高斯热源加载在圆柱面上。没有考虑对流。我是哪里设置不对或是应该怎么修改呢？命令流如下：

/NOPR
KEYW,PR_SET,1
KEYW,PR_STRUC,0
KEYW,PR_THERM,1
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
/GO
/COM,  Thermal

/PREP7
ET,1,SOLID70

ET,2,SURF152

KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,1
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,1
KEYOPT,2,9,0
KEYOPT,2,11,0
KEYOPT,2,13,0

units,si
TOFFST,273
MP,KXX,1,,15
MP,DENS,1,,7930
MP,C,1,,510
MPTEMP,1,0,1424,1426,2849
MPDATA,ENTH,1,,0,5.75906e9,7.92e9,1.3675e10

*set,rj,30e-6

cylind,rj,,0,-2*rj,-90,90
wprota,,90
vsbw,all
wprota,,-90
/pnum,line,1
lplot

mshape,0,3d                !指定划分单元的形状
mshkey,1                   !指定映射网格划分
lsel,s,loc,z,0
lesize,all,,,30
lsel,s,length,,2*rj
lesize,all,,,60,0.1
vmesh,all

asel,s,,,12
nsla,s,1
! /VIEW,1,1,2,3
! /ANG,1
! /REP,FAST
TYPE, 2
MAT,    1
REAL,
ESYS,    0
SECNUM,
esurf,all
alls
FINISH
/SOLU

ANTYPE,4

TRNOPT,FULL
LUMPM,0
ALLSEL,ALL
FLST,2,5876,1,ORDE,2
FITEM,2,1
FITEM,2,-5876
IC,P51X,TEMP,25,       !环境温度25度

*set,u,3
*set,id,3.2
*set,td,20
*set,F,0.6
*set,t_off,500e-6
*set,t_fre,1          !脉冲周期次数

*do,m,1,t_fre,1

*do,n,1,td,1
esel,s,type,,2
*get,nemax_2,elem,,num,max  !获取单元最大序号
*get,nemin_2,elem,,num,min !获取单元最小序号

*do,i,nemin_2,nemax_2,1 !每一个单元循环 i为运算时单元号
*if,esel(i),eq,1,then !如果单元被选择则执行获取中心点坐标
*SET,xsy,centrx(i)
*SET,ysy,centry(i)
*SET,zsy,centrz(i)
*SET,rr,abs(sqrt((xsy/1)*(xsy/1)+(zsy)*(zsy)))
*set,rt,14.823*n**0.2648*id**0.3733*1e-6
*if,ysy,gt,0,then
*if,rr,le,rt,then
*SET,qr,4.55*F*u*id/3.1415926/rt/rt*exp(-4.5*rr*rr/rt/rt)
sfe,i,1,hflux,,qr       !在i单元表面施加面载荷
*endif
*endif
*endif
*enddo

alls
TIME,(m-1)*(t_off+td*1e-6)+n*1e-6 !设置时间
AUTOTS,1
DELTIM,0.5e-6, , ,1
KBC,1
TSRES,ERASE
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2,
OUTRES,ALL,ALL,
ALLSEL,ALL
LNSRCH,1
ALLSEL,ALL
SOLVE
*enddo       !脉冲结束t=td

esel,s,type,,2
sfedele,all,all,hflux
!sfa,3,1,conv,10,25
TIME,(m-1)*(t_off+td*1e-6)+t_off+td*1e-6
AUTOTS,1
DELTIM,100e-6, , ,1
KBC,1
TSRES,ERASE
TIMINT,1
TINTP,0.005, , ,1,0.5,0.2,
OUTRES,ALL,ALL,
ALLSEL,ALL
LNSRCH,1
ALLSEL,ALL
SOLVE
*enddo