Статус документа
Статус документа

ГОСТ Р ИСО 7730-2009 Эргономика термальной среды. Аналитическое определение и интерпретация комфортности теплового режима с использованием расчета показателей PMV и PPD и критериев локального теплового комфорта

Приложение D
(справочное)

     
Компьютерная программа для вычисления и


Приведенная компьютерная программа на BASIC позволяет вычислить и для заданного набора вводимых переменных. Другие языки программирования также могут быть использованы, но результат должен быть проверен с помощью данной программы или по примерам, приведенным в таблице D.1.


Таблица D.1 - Примеры результатов

N п/п

Температура воздуха, °С

Средняя температура излучения, °С

Скорость движения воздуха, м/с

, %

Скорость обмена веществ, мет

Изоляция одежды, кло

1

22,0

22,0

0,10

60

1,2

0,5

-0,75

17

2

27,0

27,0

0,10

60

1,2

0,5

0,77

17

3

27,0

27,0

0,30

60

1,2

0,5

0,44

9

4

23,5

25,5

0,10

60

1,2

0,5

-0,01

5

5

23,5

25,5

0,30

60

1,2

0,5

-0,55

11

6

19,0

19,0

0,10

40

1,2

1,0

-0,60

13

7

23,5

23,5

0,10

40

1,2

1,0

0,50

10

8

23,5

23,5

0,30

40

1,2

1,0

0,12

5

9

23,0

21,0

0,10

40

1,2

1,0

0,05

5

10

23,0

21,0

0,30

40

1,2

1,0

-0,16

6

11

22,0

22,0

0,10

60

1,6

0,5

0,05

5

12

27,0

27,0

0,10

60

1,6

0,5

1,17

34

13

27,0

27,0

0,30

60

1,6

0,5

0,95

24

          

Динамические воздействия на изоляцию одежды должны быть скорректированы до ввода результирующей термоизоляции одежды ().

Переменные

Обозначения в программе

Одежда, кло

CLO

Скорость обмена веществ, мет

МЕТ

Наружная работа, мет

WME

Температура воздуха, °С

ТА

Средняя температура излучения, °С

TR

Относительная скорость движения воздуха, м/с

VEL

Относительная влажность, %

RH

Парциальное давление водяного пара, Па

РА

10

Компьютерная программа (BASIC) для вычисления

20

Прогнозируемой средней оценки (PMV) и прогнозируемого процента недовольных (PPD)

30

в соответствии с Международным стандартом ИСО 7730

40

CLS:

PRINT

"ВВОД ДАННЫХ"

ввод данных

50

INPUT

"

Одежда

(кло)

"

CLO

60

INPUT

"

Скорость обмена веществ

(мет)

"

МЕТ

70

INPUT

"

Наружная работа, обычно рядом с 0

(мет)

"

WM

80

INPUT

"

Температура воздуха

(°С)

"

ТА

90

INPUT

"

Средняя температура излучения

(°С)

"

TR

100

INPUT

"

Относительная скорость движения воздуха

(м/с)

"

VEL

110

INPUT

"

"ВЫБЕРИТЕ ИЛИ RH, ИЛИ ДАВЛЕНИЕ ВОДЯНОГО ПАРА, НО НЕ ОБА"

120

INPUT

"

Относительная влажность

(%)

"

RH

130

INPUT

"

Давление водяного пара

(Па)

"

РА

140

DEF FNPS (Т) = ЕХР (16.6536-4030.183/Т+235))

: давление насыщенного пара, кПа

150

IF PA = 0 THEN PA = RH * 10 * FNPS (TA)

: давление водяного пара, Па

160

ICL = .155 * CLO

: термоизоляция одежды, м·К/Вт

170

М = МЕТ * 58.15

: скорость обмена веществ, Вт/м

180

W = WME * 58.15

: наружная работа, Вт/м

190

MW= M-W

: наружное производство тепла в теле человека

200

IF ICL u .078 THEN FCL = 1+1.29 * ICL

ELSE FCL = 1.05 + 0.645 * ICL

: множитель площади одежды

210

HCF = 12.1 * SQR (VEL)

: коэффициент теплообмена при принудительной конвекции

220

ТАА = ТА + 273

: температура воздуха, К

230

TRA = TR + 273

: средняя температура излучения, К

240

- - - - - ВЫЧИСЛИТЕ ТЕМПЕРАТУРУ ПОВЕРХНОСТИ ОДЕЖДЫ С ПОМОЩЬЮ ПОВТОРЕНИЯ - - -

250

TCLA = TAA +(35.5-ТА)/(3.5 * ICL + .1)

: первое предположение о температуре поверхности одежды

260

Р1 = ICL * FCL

: период вычисления

270

Р2 = Р1 * 3.96

: период вычисления

280

Р3 = Р1 * 100

: период вычисления

290

Р4 = Р1 * ТАА

: период вычисления

300

Р5 = 308.7 - .028 * MW + Р2 * (TRA/100) * 4

310

XN = TLCA / 100

320

XF = XN

330

N = 0

: N: количество повторений

340

EPS = .00015

: критерии остановки при повторе

350

XF = (XF + XN)/2

360

HCN = 2.38 * ABS (100 * XF - ТАА) ^ .25

: коэффициент теплообмена при естественной конвекции

370

IF HCF > HCN THEN НС = HCF ELSE НС = HCN

380

XN = (P5 + P4 * HC - P2 * XF ^ 4) / (100 + P3 * HC)

390

N = N + 1

400

IF N > 150 THEN GOTO 550

410

IF ABS (XN-XF) > EPS GOTO 350

420

TCL = 100 * XN - 273

: температура поверхности одежды

430

- - - - - - - - - - - - - - - - - - - КОМПОНЕНТЫ ПОТЕРИ ТЕПЛА - - - - - - - - - - - - - - - - - - - - -

440

HL1 = 3.05 * .001 (5733-6.99 * MW-PA)

: потеря тепла через кожу

450

IF MW > 58.15 THEN HL2 = .42 * (MW - 58.15)

ELSE HL2 = 0!

: потеря тепла за счет потения (комфортного)

460

HL3 = 1.7 * .00001 * m * (5867-РА)

: скрытая потеря тепла из-за дыхания

470

HL4 = .0014 * m * (34 - ТА)

: потеря тепла из-за сухого дыхания

480

HL5 = 3.96 * FCL * (XN^4 - (TRA/100^4)

: потеря тепла из-за излучения

500

- - - - - - - - - - - - - - - - - - ВЫЧИСЛИТЕ PMV И PPD - - - - - - - - - - - -

510

TS = .303 * ЕХР (- .036 * m) + .028

: коэффициент переноса температурной чувствительности

520

PMV = TS * (MW - HL1 - HL2 - HL3 -
- HL4 - HL5 - HL6)

: прогнозируемая средняя оценка

530

PPD = 100 - 95 * ЕХР (- .03353 * PMV ^ 4 -
- .2179 * PMV^ 2)

: прогнозируемый процент недовольных

540

GOTO 570

550

PMV = 999999!

560

PPD = 100

570

PRINT: PRINT "OUTPUT"

: результат

580

PRINT  "Predicted Mean Vote (PMV):"

:PRINT USING "# # . #": PMV

590

PRINT "Predicted Percent of Dissatisfied (PPD):"

:PRINT USING "# # # . #": PPD

600

PRINT: INPUT "NEXT RUN (Y/N)"; RS

610

IF (RS = "Y" OR RS = "y") THEN RUN

620

END

Пример - Ввод данных

Одежда

(кло) ? 1.0

Скорость обмена веществ

(мет) ? 1.2

Наружная работа, обычно около 0

(мет) ? 0

Температура воздуха

(С) ? 19

Средняя температура излучения

(С) ? 18

Относительная влажность воздуха

(м/с) ? 0.1

ВЫБЕРИТЕ ИЛИ RH ИЛИ ДАВЛЕНИЕ ВОДЯНОГО ПАРА, НО НЕ ОБА

Относительная влажность

(%) ? 40

Давление водяного пара

(Па) ?

РЕЗУЛЬТАТ

Прогнозируемая средняя оценка

(PMV): -0.7

Прогнозируемый процент недовольных

(PPD): 15.3