傳熱學(xué)實驗--墻角matlab導(dǎo)熱問題

1、二維導(dǎo)熱物體溫度場的數(shù)值模擬姓名小明學(xué)號111111班級能動學(xué)院能動一、問題描述有一墻角模型，尺寸如圖1所示，導(dǎo)熱系數(shù)(mK),墻角內(nèi)外壁為第一類邊界條件。求解該模型的溫度分布及導(dǎo)熱量。圖1F11L1二、計算原理根據(jù)熱平衡法列出節(jié)點方程，各方向?qū)雴卧w的熱量之和為零。內(nèi)節(jié)點和絕熱邊界點(圖1點劃線上的點)的方程形式不同。圖2圖2所示的內(nèi)節(jié)點和絕熱邊界節(jié)點方程如下內(nèi)節(jié)點：+=九NSEW=0TOC o 1-5 h z(tt)Ax(tt)Ax(tt)/y(tt)/y*v*-*v*v*v*ij1I,J+ij1I,J+i+1jI,J+i1jI,JAyAy&絕熱邊界點：+=NSEW(tt)Ax(tt)A

2、x(tji,j+ji,j+0+i1Jt)Ayi,jAx=0三、計算過程用語言編寫計算程序，初取網(wǎng)格步長Ax二Ay二0.1m5HTafunct:L=M;endform=6;16forn=6:12T皿nj=0:end.end1(16,2:5)=24,18,12.6;1(2:5,12)=24,18,12,6;vhile1for1=2:15forj=2:5；temp=I-ijj)；心jJ=(T(i-Lj)+T(1+LjHT(i,j-l)十心j+lD/4；eps=Iijj)-temp:end.endrwr?ifeps1e8bxsak:endendwhile1tenip=rr16j);1(16,j)-Q*

3、T(1.5,ji+TClG,j-D+Ttie,j+l)/4;eps=T16jj-temi);-endife-psle-8break;endendwhile1for1=2:5forj=2:11二:temp=T(ijj);龍jJ+rCi+UjJ+KLj-D+Tti,j+L)/4;epIijj)te皿p;end-endifepsle-8break;end-end-vhile1fori=2:5teHp=r(i.j12J；Ki,L2)=Q*T(l,llHKi-L12)+T(i+1,12)/4;eps=Tij12一七emp;?n.difepsleSbreak;end-endpcolorI1:shading

4、interpfigurejcontour(TJlu);運行結(jié)果：jO.O：jO.O：30.0：30.0：30.0：3C.0：3：.(03：.CD3：.?)3：.?)3：.(O30.0：jO.O：28.9:27.8=26.900526.9：26.02盯門2.79膽.邪24.13乩0524.0:jO.O：r茫25.2甌：21.S72C.3；三門IE.79IE.3IE.22IE.0318.OEjO.O：26.223.4?20.217.3:1E.0：1-.6212.E6二労12.23二1)12.OEjO.O：藥5221.16.612.3：E.6j7.216.59c.ZD5.15c.056.05jO.O

5、：24.4918.812.85F，用C.O：:.CD：.(O0.0：jO.O：24.2718.4512.466.2=：.(O:.CD:0.0：jO.O：14二18.2512.2c6.1:C.O：.(O:.CD：.(O0.0：jO.O：24.OE18.匸二6.0E：.(O:.CD2/0:）0.0：jO.O：24.0518.0?12.0?6.0:：:.CD：.：D2/0:）0.0：jO.O：24.0218.0412.06.02C.O：.(O:.CD：.(O0.0：jO.O：24.0118.012.026.01C.O：:.CD：.(O0.0：jO.O：24.0118.0112.016.01：.(O:.CD:0.0：jO.O：24.0：18.0112.016.0：C.O：.(O:.CD：.(O0.0：jO.O：24.0：18.0212.0：6.0：.(O:.CD2/0:）0.0：jO.O：24.0：18.0：12.0：6.0：:.CD：.：D2/0:）0.0：圖1：各個點的溫度數(shù)值圖2：分層設(shè)色等溫線分布圖3：等溫線分布（每兩條線間隔為三度）四、小結(jié)本次數(shù)值模擬是運用matlab程序用于數(shù)值計算。小組成員共同討論并復(fù)