马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
毕业设计中的热分析程序:
程序介绍:在平板上(面2)加高斯热源(取1/4平板),平板在煤油中加工,所以想加上对流,所以在面2,4,6上施加了对流载荷(表面上生成了surf152单元,经检查确是生成了),结算后表明高斯热源加上了,但是对流载荷总加不上(运算中出了载荷覆盖的警告)。请您帮我看看程序加载上出了什么问题。最好详细点!谢谢
finish
/clear
/title,transient analysis(gauss heat flux)
/filnam,transientgauss
/TRIAD,RBOT
/UNITS,SI
/PREP7
et,1,solid70
et,2,surf152
keyopt,2,8,2 !要在surf152上施加对流
X1=0.2 !输入平板的能量比例
U=25 !电压
I=3 !电流
Ton=100 !脉冲宽度
pi=3.1415926
Q=3*X1*U*I
lisannum=4 ! 控制载荷布的数量
dividenum=20 !控制网格划分的密度
block_x=0.0005
block_y=0.0005
block_z=-0.00005
BLOCK,0,block_x,0,block_y,block_z,0,
/VIEW,1,1,1,1
save
mptemp,1,20,250,500,750,1000,1500,1700,2500 !定义材料属性
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500
SAVE
FLST,5,2,4,ORDE,2 !分网部分
FITEM,5,4
FITEM,5,-5
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,dividenum,5, , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,1
FITEM,5,8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,dividenum,5, , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,1
FITEM,5,8
CM,_Y1,LINE
LSEL, , , ,P51X
*GET,_z1,LINE,,COUNT
*SET,_z2,0
*DO,_z5,1,_z1
*SET,_z2,LSNEXT(_z2)
*GET,_z3,LINE,_z2,ATTR,NDNX
*GET,_z4,LINE,_z2,ATTR,SPNX
*get,_z6,line,_z2,attr,kynd
*IF,_z3,GT,0,THEN
*IF,_z4,NE,0,THEN
LESIZE,_z2,,,_z3,1/_z4,,,,_z6
*ENDIF
*ENDIF
*ENDDO
CMSEL,S,_Y1
CMDELE,_Y1
!*
FLST,5,3,4,ORDE,3
FITEM,5,9
FITEM,5,-10
FITEM,5,12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,dividenum/5,5, , , ,1
!*
FLST,5,1,4,ORDE,1
FITEM,5,10
CM,_Y1,LINE
LSEL, , , ,P51X
*GET,_z1,LINE,,COUNT
*SET,_z2,0
*DO,_z5,1,_z1
*SET,_z2,LSNEXT(_z2)
*GET,_z3,LINE,_z2,ATTR,NDNX
*GET,_z4,LINE,_z2,ATTR,SPNX
*get,_z6,line,_z2,attr,kynd
*IF,_z3,GT,0,THEN
*IF,_z4,NE,0,THEN
LESIZE,_z2,,,_z3,1/_z4,,,,_z6
*ENDIF
*ENDIF
*ENDDO
CMSEL,S,_Y1
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,';VOLU';
CMSEL,S,_Y
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH !分网结束
*dim,R,array,lisannum,1,1 !高斯热流密度的计算
tempr=0.11*I*I+0.73*I+8.52
*IF,tempr,GE,Ton,then
maxR=2.85*I**0.53*Ton**0.38
*ELSE
maxR=2.85*I**0.53*tempr**0.38
*ENDIF
maxenum=dividenum+1
*dim,flux2,table,maxEnum,maxEnum,lisannum,x,y,time
*DO,k,1,lisannum,1
t=k*Ton/lisannum
R(k)=2.85*I**0.53*t**0.38
*IF,R(k),GE,maxR,then
R(k)=maxR
*ENDIF
*do,j,1,maxEnum,1
flux2(0,j,k)=(j-1)*block_x/dividenum
*enddo
*do,i,1,maxEnum,1
flux2(i,0,k)=(i-1)*block_x/dividenum
*enddo
*do,i,1,maxEnum,1
*do,j,1,maxEnum,1
xcenter=0
ycenter=0
distance=((i-1)*block_x/dividenum-xcenter)**2+((j-1)*block_x/dividenum-ycenter)**2
*if,distance,le,(R(k)*1E-6)**2,then
flux2(i,j,k)=(Q/(pi*R(k)*R(k)*1E-12)*exp(-3*distance/R(k)*R(k)*1E-12)) !????????
*else
flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*Ton/lisannum
*enddo
/SOLU
ANTYPE,trans
TRNOPT,FULL
NROPT,AUTO, ,
EQSLV,
SOLCONTROL,ON
AUTOTS,on
KBC,1
OUTRES,ALL,ALL,
TIME,1e-6
DELTIM,1E-6,1E-6,1E-6
TUNIF,25
kbc,1
solve !先施加25度的开始温度计算,只是第一步起作用
type,2
asel,s,,,2,6,2
NSLA,S,1
ESLN,S,0
ESURF !在2,4,6面上生成表面效应单元
esel,s,type,,2
sf,all,conv,-1,25 !施加对流(hf已经在前面材料属性部分定义了)
allsel,all
nsel,all
eplot
totaltime=Ton*1E-6
inctim=Ton*1E-6/lisannum
autots,1
*do,t,2e-6,totaltime,inctim
time,t
deltim,inctim,inctim,inctim
kbc,0
outres,all,all,
tsres,erase
asel,s,loc,z,0,0
sfa,all,1,hflux,%flux2% !施加高斯热源
solve
OUTRES,ALL,ALL,
/PSF,HFLUX,,2
/REPLOT
*enddo
finish
疑问:表面效应单元生成了,但是并没有生成新的节点,是不是实体载荷和表面效应单元载荷都传到了节点上所以覆盖了?
|