ReDim b(1 To cs+1) As Variant '导热系数特征系数b
ReDim z(1 To cs+1) As Variant '导热系数λ
ReDim g(1 To cs+1) As Variant '是否为硅酸铝制品
‘条件输入
tn=val(text2.text) ‘内壁温度
Select Case Combo1.ListIndex
Case 0
Ak = 2.56 '垂直的平壁热面
Case 1
Ak = 3.26 '水平向上的平壁热面
Case Else
Ak = 1.98 '水平向下的平壁热面
End Select
‘赋值数组s、a、b、g
双击MSFlexGrid控件，选择材质，赋值给a、b、g及s,代码略。
'按照公式进行计算
r1=0.03: r2=0.5: r3=1 ‘设定各循环步长及最终温度误差范围
tt=tw+1
For w = tt To tn Step r1
qs = Ak * (w - tw)^1.25 + 4.54 * (((w + 273) / 100)^4 - ((tw + 273) / 100)^4)
For n = 1 To cs
'计算中间各层开始>>
t(1) = w
p = t(n) + 1
For i = p To tn Step r2
If g(n) = 1 Then '若为硅酸铝纤维制品
z(n) = a(n) + 0.000001 * b(n) * ((i + t(n)) / 2) ^ 2
Else
z(n) = a(n) + 0.001 * b(n) * ((i + t(n)) / 2)
End If
t(n + 1) = qs * (s(n) / 1000) / z(n) + t(n) 's（n）表示第n层内衬厚度mm
If Abs(i - t(n + 1)) < r1 Then
Exit For
End If
Next i
'中间各层计算结束<<
Next n
If Abs(tn - t(1 + cs)) < r3 Then
Exit For
End If
Next w
‘输出结果
通过Frame控件内的text控件，显示计算结果（代码略），同时可导出到记事本。
3）界面设计
该程序采用VB进行编写，利用Access作为后台数据库，存放耐火材料导热系数的特征参数及表达式。界面采用Frame、Combo、Text、Command、MSFlexGrid及Option等常用控件，操作简便、直观，见图1。
