Вестник МГТУ, 2023, Т. 26, № 1.

Вестник МГТУ. 2023. Т. 26, № 1. С. 45-56. DOI : https://doi.org/10 .21443/1560-9278-2023-26-1-45-56 необходимости произвести или проверить аналогичные вычисления следует скопировать программу из рис. 4, вставить в поле редактора https://yohan.es/swbasic (или аналогичного интерпретатора языка программирования) и ввести контрольные значения. 10 REM компьютернаяпрограмма (BASIC) длявычисления 20 REM прогнозируемойсреднейоценки (РМѴ) ипрогнозируемогопроцентанедовольных (PPD) 30 REM всоответствиисмеждународнымстандартомISO 7730 40 CLS : PRINT "ВВОДДАННЫХ" : REM вводданных 50 INPUT 60 INPUT 70 INPUT 80 INPUT 90 INPUT 100 INPUT 110 PRINT 120 INPUT 130 INPUT . 140 REM давление насыщенного водяного пара при температуре Т в кПа = ЕХР (16.6536-4030.183/(Т+235) ) 145 е=2.718281828459045 150 IF РА=0 THEN PA=RH*10* (еА (16.6536-4030.183/(ТА+235))): REM давление водяного пара, Па 160 ICL=.155*CLO : REM термоизоляция одежды, кв.м*К/Вт 170 М=МЕТ*58.15 : REM скорость обмена веществ, Вт/кв.м 180 W=WME‘58.15 : REM наружная работа, Вт/кв. м 190 MW=M-W : REM внутреннее производство тепла в теле человека 200 IF ICL<=0.078 THEN FCL=1+1.29*ICL ELSE FCL=1.05+0.645*ICL 201 REM ICL - множительплощадиодежды 210 HCF=12.1*SQR(VEL) : REM коэффициент теплообмена при принудительной кон­ векции 220 ТАА=ТА+273 : REM температуравоздуха, К 230 TRA=TR+273 : REM средняятемператураизлучения, К 240 R E M .................. ВЫЧИСЛЕНИЕТЕМПЕРАТУРЫПОВЕРХНОСТИОДЕЖДЫПОСЛЕДОВАТЕЛЬНЫМИПРИБЛИЖЕНИЯМИ........... Одежда (кло) Н CLO Скорость обмена веществ (мет) ", МЕТ Наружная работа, обычно рядом с 0 (мет) п 1 WM температура воздуха (град.С) ”г ТА Средняя температура излучения (град.С) г TR Средняя скорость движения воздуха (м/с) ", VEL ВВЕДИТЕ ИЛИ ВЛАЖНОСТЬ ИЛИ ДАВЛЕНИЕ ВОДЯНОГО ПАРА, НО НЕ Относительная влажность (%> ( RH Давление водяного пара (Па) й РА 250 TCLA=TAA+(35.5-ТА)/(3.5*ICL+.1) 310 XN=TCLA 320 XF=XN 330 N=0 340 EPS=.0000015 350 XF=(XF+XN)/2 360 HCN=2.38*(ABS(XF-TAA)A.25) rem первоепредположениеотемпературеповерхности одежды : REM номеритерации : REM критерийзавершенияциклапоискаприближенного значения : REM коэффициенттеплообменаприестественнойконвек­ ции 370 IF HCF>HCN THEN HC=HCF ELSE НС= HCN 380 P1=ICL*FCL*((3.96*(10A-8)*((XFA4)-(TRAA4)))+(HC*(XF-TAA)) 381 XN=35.7-(.028*MW)-Pl+273 390 N=N+1 400 IF N>150 THEN GOTO 550 410 IF ABS(XN-XF)>EPS THEN GOTO 350 420 TCL=XN-273 : REM температураповерхностиодежды 430 R E M ........................ ВЫЧИСЛЕНИЕКОМПОНЕНТОВПОТЕРИТЕПЛА......................... 440 HL1=3.05*.001*(5733-6.99*MW-PA) 450 IF ММ>58.15 THEN HL2=.42*(MW-58.15) ELSE HL2=0! 451 460 HL3=1.7*.00001*M*(5867-PA) 470 HL4=.0014*M*(34-TA) 480 HL5=3.96*(10A-8)*FCL*(XNA4-TRAA4) 490 HL6 = FCL*HC*(TCL-TA) 500 R E M ------------------------ ВЫЧИСЛЕНИЕ PMV И PPD- 510 TS=.303*(eA (-.036*M))+.028 REM потерятеплазасчетпотения (комфортного) REM скрытаяпотерятеплачерездыхание REM потерятеплаиз-засухогодыхания REM потерятеплаиз-заизлучения 520 PMV=TS*(MW-HL1-HL2-HL3-HL4-HL5-HL6) 530 PPD=100-95*(еА(-.03353*(РМѴА4)-.2179*(РМѴА2))) 540 GOTO 570 550 PMV=999999! 560 PPD=100 570 PRINT : PRINT "РЕЗУЛЬТАТ" 580 PRINT "прогнозируемаясредняяоценка (PMV):", рмѵ 590 PRINT "прогнозируемыйпроцентнедовольных (PPD):", : rem коэффициентпереносатемпературнойчувствительно­ сти : REM прогнозируемаясредняяоценка : REM прогнозируемыйпроцентнедовольных : rem выводрезультата PPD Рис. 4. Исправленный текст компьютерной программы на языке BASIC для вычисления PMV и PPD Fig. 4. Computer program in BASIC for calculating PMV and PPD 51