S?Sr?SB 元/(m?a) (8)
对于平壁：
S'?S'r?S'B 元/(m2?a) (8′)
为了求取保温经济厚度，将(1)或(1′)式及(7)或(7′)式代入(8)或(8′)式，并令其导数为零，即
?S?0?D 管道 1 ?S??0 平壁 ??1
则得下式： 对于管道：
??1?nC?10?3=?K?SB1D1?10-3+SB2? ? ??2?1DD?2?211000?? ?In1?In1?In?2?D2?D(D?2?)?2112?1?2?????
??1??2?21000?? ?(t?t0)????2??62?D?D(D?2?)?(D?2?)?211212?11???元?h/(?10kcal) ?(9)
对于平壁：
2nC?10?6??1?21?????KSB1?1000??11000?2???= 元?h/(?106kcal) (9′)
1(t-t0)?1
4.热损失计算 对于管道：
q?t?t0 kcal/(m?h)
DD?2?2111000In1?In1?2??1D2??2D1?(D1?2?2)? (10)
对于平壁：
q?t?t0?11000?1??21000?2?1 kcal/(m2?h) (10′)
? 二、程序使用说明
根据以上计算原则，用求解非线性方程组的方法可直接求保温经济厚度。本计算程序系在TQ-16机上用ALGOL-60语言编制而成。现已移植到MZ-80B微型计算机上，实现了人机对话，屏幕显示，汉字输出的功能。现仅就程序的使用方法举例说明如下： (一)原始数据
本计算程序的原始数据应按附表1.1及附表1.2的格式要求填写。填写附表的注意事项是：
附表 1.1 简 单 变 量 表 序 号 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 名 称 2 保温对象年运行小时数 保温结构投资年费用率 保温材料密度 保护层材料密度 保温层单价 保护层单价 保温材料导热系数基数 温度系数(一次项) 温度系数(二次项) 保护层厚度 周围环境温度 保温对象周围平均风速 保护层导热系数 计算管道根数 代 号 3 n K G1 G2 SB1 SB2 λ0 单 位 4 h/a — kg/m3 kg/m3 元/m2 元/m2 kcal/(m·h·℃) mm ℃ m/s kcal/(m·h·℃) 根
数 值 5 8000 0.17 350 1000 245 4.82 0.064 0.0001 0 20 20 0 0.2 9 备 注 6 金属保护壳填零 室内管道填零 A A1 δ2 t0 W λ2 N 附表 1.2 数 组 表 管子外径D mm 355.6 508 159 热介质温度t ℃ 540 540 320 热 单 价 C 元/×106kcal 9.98 9.01 6.37 管子长度L m 150 100 80 325 159 2020 426 508 2020 450 100 320 110 540 320 9.64 5 6.37 5 0 0 120 50 1 60 100 1 注：(1)表中行数系管子根数，本例是9根管。
(2)凡管子外径大于2000mm，程序自动按平壁计算有关数据。 (3)热单价填零时，则按表面温度不超过50℃计算有关数据。 (1)各项数据必须按表列的单位填写； (2)附表1.1为简单变量，附表1.2为数组； (3)表列各项必须填满。
(二)TQ-16机数据书写、穿孔及上机操作 1.数据书写及穿孔
本程序原始数据采用N格式输入，按十进制形式书写。这种格式的数据，根据华北电力设计院TQ-16机的编译功能，采用下列两种类型： (1)BS型(分段型)数据
这种类型的数据要分成几个数据段。前面的简单变量作为一个数据段。后面的每一个数组各为一个数据段。因此，本程序原始数据要写五个数据段。各数据段中的数据之间用逗号隔开，最末一个数据之后要书写分号及15个逻辑乘∧∧∧??，在各数据段首末加6*6*，作为起止符合。头一个数据段之前和末一个数据段之后各空两格，数据段之间空半格。 (2)ZN型(紧缩型)数据
ZN型数据一律用八单位纸带穿孔。它不需要分段，程序中所有输入数据严格按上面填写的顺序列在一起，构成一个数据段。这种类型的数据还可采用下面三种缩写形式： 1)数据连续相同时的缩写形式 ＜数据＞(＜重复个数＞)
例如：有五个连续相同的数据540，540，540，540，540，可缩写成540(5)。 2)数据连续等差时的缩写形式
＜等差首项＞(＜等差数据个数＞，＜步长＞)。
例如：对五个连续的数据20，18，16，14，12，可缩写成20(5，-2)。 3)数据连续相同段的缩写形式 ＜重复次数＞［＜数据重复段＞］。
例如：对三个连续相同的数据段50，100，150，50，100，150，50，100，150，可缩写成3［50，100，150］。
数组表(附表1.2)的数据输入顺序是先列后行(横行竖列)。 2.上机操作
输入命令主要靠电传，基本操作是：
(1)将源程序纸带置于某台光电机上，并在控制台上拨好相应的光电机选台板键，将输
入方式板键打向八单位方向，依照顺序按“停机”、“置零”、“启动”按钮，由电传印出DDGZ(等待工作)。
(2)发布编译命令： BX FY
命令执行，源程序输入。编译结束，印出BYJS(编译结束)。
(3)将数据纸带置于同一台光电机上，根据数据穿孔类型采取不同命令： 1)采用BS型数据时，发布命令BS。 2)采用ZN型数据时，发布命令ZN。
当发布ZN命令后，一俟数据输入完毕，即由八十行输出打字机输出已输入的数据，并由电传印出SJJS(数检结束)，停机。按“启动”按钮，继续运行。 (三)MZ-80B微型机数据输入及上机操作
(1)打开软盘驱动器，将软盘放入1号软盘驱动器，接通电源，则屏幕显示如附表1.3。
附表1.3 MZ-80B微型机屏幕显示 **MONTTOR SB-1510** BASIC interpreter SB-6510 COPY right 1981 by SHARP 38636 Bytes Ready (2)发布命令DIR(或按动能键F7)，则屏幕显示本软盘的文件索引。 (3)本程序文件名是：“BWJJS”
发布命令LOAD“BWJJS”，则程序调入内存。
(4)将附表1.1及附表1.2的原始数据从4000语句开始DATA格式输入。数组表（附表1.2）的输入顺序是先行（横行竖列）。
(5)发布命令RUN“BWJJS”，则程序执行。 (四)计算结果表
根据附表1.1及附表1.2所列原始数据进行计算，结果如附表1.4。
附表1.4 保温经济厚度计算结果表 温度 ℃ 热价 元/×10６保长 度 m 温层 厚 度 mm 150 100 80 120 50 190 36.19 200 36.56 100 33.78 170 34.14 40 27.92 394.35 493.07 172.75 313.21 69.38 表面 ℃ 热损失 kcal/(m2·h) 保温层 体 积 m3 保护层 体 积 m3 保温结表面积 m2 构 重 量 kg/m，kg/m2 365.49 297.82 100.28 265.78 43.83 161.46 214.01 52.29 135.57 25.03 管道 外径 mm 温度 kcal/(m·h)， kcal 9.98 9.01 6.37 9.64 5 355.6 540 508 159 325 159 540 320 450 100 48.85 44.49 6.51 31.72 1.25 7.121 5.83 1.905 5.164 0.813 2020 426 508 2020 320 110 540 320 6.37 5 0 0 1 60 100 1 160 50 90 33.9 28.66 49.96 138.98 153.92 598.92 229.62 0.16 4.49 31.01 0.09 0.02 2.058 5.202 0.02 1 106.69 266.41 1 76 60.48 160.55 51.5 150 49.48 168.56 28.136 1448.3 注：表中末行的三个数值分别为保温层体积、保护层体积、表面积的累计值，供统计保温工程量时使用。
附录二 保温经济厚度计算程序参考数据
在进行保温经济厚度计算时，应根据工程具体条件合理地确定有关数据。如果某些数据难以确定时，可参考本附录所列的数据。 一、简单变量
1.保温对象年运行小时数n(h/a)
全年运行的管道 采暖管道 东北地区 华北地区
2.保温结构投资年费用率K
8000 4000 3000
在技术经济比较方法未统一前，投资年费用率暂取0.17。 3.保温材料密度(容重)G1(kg/m3)
微孔硅酸钙 岩棉和矿渣棉 水玻璃珍珠岩 水泥珍珠岩 A级焙烧硅藻土 水泥蛭石
4.保护层材料密度G2(kg/m3)
灰泥抹面保护层 镀锌铁皮保护层
铝或铝合金保护层2850(以LC4为准) 5.保温材料单价SB1
(元/m3)
1000 7850
200～250 100～150 250～300 350～400 400～500 450～500
由电力局颁发的装置性材料预算价格加上安装费用得出。
预算价格包括保温材料的出厂价格、运杂费并考虑安装运输中的损耗率后得出。运杂费包括材料由产地或发货地点运到发电工程现场仓库的运输费、采购费、供销部门手续费以及包装费等。
安装费用(元/m3)：
设备
5万kW以下机组 5～30万kW机组
100 67.9 61.7