實(shí)驗(yàn)大作業(yè)傳熱主課老師

傳熱大

6題目0.1m與壁內(nèi)溫度一致，tf1=tf2=t0=5h1=11W/m2K、h2=23W/m2K,壁的導(dǎo)熱系數(shù)=0.43W/mKa=0.3437×10-6m2/s。如果一側(cè)的環(huán)境tf150℃并維持不變，計(jì)算在其它參數(shù)不變的條件下，平壁內(nèi)溫度分布及分析數(shù)值解法個(gè)子區(qū)域，將網(wǎng)絡(luò)線的交點(diǎn)作為節(jié)點(diǎn)每個(gè)節(jié)點(diǎn)就代表以它為中心的子區(qū)域（元體或稱為控解題過(guò)程1、劃分網(wǎng)首先將問(wèn)題中的無(wú)限大平板沿有限長(zhǎng)方向?qū)嵭袝r(shí)間---x為空間坐標(biāo)，1001011101為邊界點(diǎn)，2~99為內(nèi)節(jié)t為時(shí)間坐標(biāo)，以為時(shí)間步長(zhǎng)，將時(shí)間坐標(biāo)上的計(jì)算區(qū)域劃分為（I—1）等份，I個(gè)時(shí)間節(jié)點(diǎn)。但時(shí)間點(diǎn)個(gè)數(shù)不能提前給出，由達(dá)到穩(wěn)態(tài)所需要的時(shí)間決定。選x=1mmBi＝h1x

2.55814102 h

Bi＝ ＝

Fo

2(1Bi Bi1<Bi2,Bi1

x

＝20.343710312.55814102＝1.4187s故選擇=0.5s2、建立節(jié)點(diǎn)離散方程htf，為第三類邊界條件。采取

ti1ti(12Fo 2Fo)2Fo(tiBit 1f

ti1

)

Bi

2fnti1n

12Fo

nn0t1nn0

＝5℃（n＝1,2,3…99，100,101）qih ti 1f qih(ti 2 f3、穩(wěn)態(tài)數(shù)值計(jì)q

tf1tf 1

1

12h12

q

38.85-t(I,1)<0.014流程框5程溫度分布程 %計(jì)算網(wǎng)格數(shù) %5,n %設(shè)置時(shí)間節(jié)點(diǎn)初值，iwhile38.85-t(i- %0.01t(i,1)=t(i-1,1)*(1-2*Fo*Bi1-2*Fo)+2*Fo*(t(i-1,2)+Bi1*50); forn=2:100t(i,n)=Fo*(t(i-1,n+1)+t(i-1,n-1))+(1-2*Fo)*t(i-1,n); t(i,101)=t(i-1,101)*1-2*Fo*Bi2-2*Fo+2*Fo*t(i-1,100)Bi2*5);% %1 i 畫(huà)出不同時(shí)間點(diǎn)溫度在平面上的分for %120dt1holdon;forj=1200:600:9000holdon;forj=7200:1200:19200holdon;forj=18000:6000:79200holdon;

holdoff; xlabel(‘x/m’);%x軸 %ygridon;%畫(huà)出熱流密度隨時(shí)間的變化iq1(i)=11*(50-t(i,1))%q1q2(i)=23*(t(i,101)-5)%q26、計(jì)算結(jié)果與的進(jìn)行，物體左邊界溫度升高，q1逐漸減小，且變化率也在減小。說(shuō)明在外界溫度未變的情況下，物體邊界溫度變化趨于緩慢，此結(jié)論可在溫度分布圖中得到驗(yàn)證。q1最終在接近度未變的情況下，物體邊界溫度變化趨于緩慢，此結(jié)論可在溫度分布圖中得到驗(yàn)證。q2最選x=10mm，則Bi＝h1x

h

Bi＝ ＝

Fo

2(1Bi Bi1<Bi2,Bi1

x

＝20.343710310.255814＝1.4187s故需選擇=0.1s

實(shí)驗(yàn)總體操作過(guò)程并掌握了的基本使用?？芍^收獲良多。我的實(shí)驗(yàn)過(guò)程可分為三個(gè)階段。第二階段：的學(xué)習(xí)。在此之前只聽(tīng)說(shuō)功能很強(qiáng)大，但沒(méi)有用過(guò)，因此書(shū)了循環(huán)和畫(huà)圖的基本知識(shí)。第三階段：寫(xiě)程序，畫(huà)圖并調(diào)試。最初的程序沒(méi)有要求顯示每次循環(huán)計(jì)算得到t(i,1)和循環(huán)次數(shù)，因而在計(jì)算過(guò)程中看不到任何變化，心里很著急，也不知道程序是否1min