给出一段高斯热源程序，有一个问题看似无法解决，大家帮忙看看
/clear/units,si
/config,nres,10000
/prep7
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPDATA,KXX,1,,122
MPDATA,KXX,1,,126
MPDATA,KXX,1,,130
MPDATA,KXX,1,,138
MPDATA,KXX,1,,147
MPDATA,KXX,1,,159
MPTEMP,,,,,,,,
MPTEMP,1,25
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPDATA,C,1,,921
MPDATA,C,1,,1005
MPDATA,C,1,,1047
MPDATA,C,1,,1089
MPDATA,C,1,,1140
MPDATA,C,1,,1200
MPTEMP,,,,,,,,
MPTEMP,1,25
MPDATA,HF,1,,15
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,200
MPTEMP,3,300
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,600
MPTEMP,7,700
MPDATA,DENS,1,,2650
MPDATA,DENS,1,,2616
MPDATA,DENS,1,,2596
MPDATA,DENS,1,,2574
MPDATA,DENS,1,,2550
MPDATA,DENS,1,,2479
MPDATA,DENS,1,,2338
allsel,all
save
FINI
/prep7
/units, si
rb=0.002
length=0.22
width=0.11
height=0.004
et,1,solid70
et,2,solid90
block,0,length,0,width,0,height
Block,0,length,0,3*Rb,0,height
Block,0,length,0,Rb*6,0,height
Block,0,length,0,Rb*12,0,height
VOVLAP, ALL
/view,1,1,1,1
Lesize,4,rb
Lesize,5,rb
Lesize,9,rb
Lesize,10,rb
Lesize,14,rb
Lesize,19,rb
Lesize,23,rb
Lesize,24,rb
Lesize,49,rb
Lesize,50,rb
Lesize,51,rb
Lesize,52,rb
Lesize,31,2*rb
Lesize,26,2*rb
Lesize,43,4*rb
Lesize,38,4*rb
Lesize,2,4*rb
Lesize,7,4*rb
Lesize,53,4*rb
Lesize,54,4*rb
Lesize,55,4*rb
Lesize,56,4*rb
Lesize,11,rb
Lesize,12,rb
Lesize,47,rb
Lesize,48,rb
save
type,1
mshkey,1
mshape,0
vmesh,5
vmesh,6
type,2
mshkey,0
mshape,1
vmesh,7
vmesh,8
Save
/solu
antype,4
tunif,25
nsel,s,loc,x,0
sf,all,conv,15,25
allsel,all
nsel,s,loc,x,length
sf,all,conv,15,25
allsel,all
nsel,s,loc,z,0
sf,all,conv,15,25
allsel,all
nsel,s,loc,z,height
sf,all,conv,15,25
allsel,all
nsel,s,loc,y,width
sf,all,conv,15,25
allsel,all
v=0.002
tinc=rb/v
tm=length/v
U=13
I=120
ef=0.7
p=U*I*ef
r=0.006
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1
max_x=length/rb+1
max_y=width/rb+1
*dim,flux2,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
*do,i,1,max_x,1
flux2(i,0,k)=(i-1)*rb
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*rb
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/exp(3*distance**2/(r**2))
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
*do,t,tinc,tm,tinc
time,t
nsel,s,loc,z,height
sfdele,all,conv 问题在这里，将与空气的热传递载荷去掉了，让后加上了高斯热源载荷，但在实际问题中，与空气的热传递是无法避免的，但是在ANSYS中同时不能加载两个热载荷的。这就是问题所在。
一个解决方法就是把这个平板的表面划分成为两部分，一部分加载高斯热源，另一部分继续加载与空气的热交换，不知道各位有没有更好的解决方案？
sf,all,hflux,%flux2%
allsel,all
nlgeom,on
nropt,full,,on
autots,1
nsubst,25,,,1
kbc,1
lnsrch,1
tsres,erase
tintp,,,,1,,,,,
solve
nsel,s,loc,z,height
sfdele,all,hflux
sf,all,conv,15,25
allsel,all
*enddo
!冷卻階段以1秒為步長冷卻到秒
*DO,I,1,20
allsel,all
TIME,1*I+tm
AUTOTS,1
nsubst,25,,,1
LNSRCH,1
KBC,1
TSRES,ERASE
tintp,,,,1,,,,,
SOLVE
*ENDDO
!冷卻階段以10秒為步長冷卻到秒
*DO,I,1,20
allsel,all
TIME,1*I+tm+20
AUTOTS,1
nsubst,25,,,1
LNSRCH,1
KBC,1
TSRES,ERASE
tintp,,,,1,,,,,
SOLVE
*ENDDO
!冷卻階段以30秒為步長冷卻到秒
*DO,I,1,20
allsel,all
TIME,30*I+tm+20+200
AUTOTS,1
nsubst,25,,,1
LNSRCH,1
KBC,1
TSRES,ERASE
tintp,,,,1,,,,,
SOLVE
*ENDDO
!冷卻階段以60秒為步長冷卻到秒
*DO,I,1,20
allsel,all
TIME,60*I+tm+20+200+600
AUTOTS,1
nsubst,25,,,1
LNSRCH,1
KBC,1
TSRES,ERASE
tintp,,,,1,,,,,
SOLVE
*ENDDO
esel,all
esln,s,1,all
eplot
allsel,all
save
finish The question is that if the air convection is necessary? You should estimate the its magnitude and compare it to the thermal resistances of other heat transfer path and then make juddeif the air convection is required...
页:
[1]