实验大作业传热主课老师

传热大

6题目0.1m与壁内温度一致，tf1=tf2=t0=5h1=11W/m2K、h2=23W/m2K,壁的导热系数=0.43W/mKa=0.3437×10-6m2/s。如果一侧的环境tf150℃并维持不变，计算在其它参数不变的条件下，平壁内温度分布及分析数值解法个子区域，将网络线的交点作为节点每个节点就代表以它为中心的子区域（元体或称为控解题过程1、划分网首先将问题中的无限大平板沿有限长方向实行时间---x为空间坐标，1001011101为边界点，2~99为内节t为时间坐标，以为时间步长，将时间坐标上的计算区域划分为（I—1）等份，I个时间节点。但时间点个数不能提前给出，由达到稳态所需要的时间决定。选x=1mmBi＝h1x

2.55814102 h

Bi＝ ＝

Fo

2(1Bi Bi1<Bi2,Bi1

x

＝20.343710312.55814102＝1.4187s故选择=0.5s2、建立节点离散方程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、稳态数值计q

tf1tf 1

1

12h12

q

38.85-t(I,1)<0.014流程框5程温度分布程 %计算网格数 %5,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 画出不同时间点温度在平面上的分for %120dt1holdon;forj=1200:600:9000holdon;forj=7200:1200:19200holdon;forj=18000:6000:79200holdon;

holdoff; xlabel(‘x/m’);%x轴 %ygridon;%画出热流密度随时间的变化iq1(i)=11*(50-t(i,1))%q1q2(i)=23*(t(i,101)-5)%q26、计算结果与的进行，物体左边界温度升高，q1逐渐减小，且变化率也在减小。说明在外界温度未变的情况下，物体边界温度变化趋于缓慢，此结论可在温度分布图中得到验证。q1最终在接近度未变的情况下，物体边界温度变化趋于缓慢，此结论可在温度分布图中得到验证。q2最选x=10mm，则Bi＝h1x

h

Bi＝ ＝

Fo

2(1Bi Bi1<Bi2,Bi1

x

＝20.343710310.255814＝1.4187s故需选择=0.1s

实验总体操作过程并掌握了的基本使用。可谓收获良多。我的实验过程可分为三个阶段。第二阶段：的学习。在此之前只听说功能很强大，但没有用过，因此书了循环和画图的基本知识。第三阶段：写程序，画图并调试。最初的程序没有要求显示每次循环计算得到t(i,1)和循环次数，因而在计算过程中看不到任何变化，心里很着急，也不知道程序是否1min