在热分析时不收敛,温度还会比常温还低,不知道怎么回事。热原模型是变半径的高斯热源加载在圆柱面上。没有考虑对流。我是哪里设置不对或是应该怎么修改呢?命令流如下:
/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@wuming524 |