10 PRINT"UNE-EN ISO 7730:2006 11 PRINT"PROGRAMA INFORMÁTICO PARA EL CALCULO DE LOS PMV (Voto Medio Estimado) Y PPD (Porcentaje Estimado de Insatisfechos) 20 ' Voto Medio Estimado (PMV) y Porcentaje Estimado de Insatisfechos (PPD) de acuerdo con la Norma Internacional ISO 7730 21 ' 1 clo = 0.155 m2.§K/W (aislamiento confortable vestimenta a ta=20 §C y hr=50 % durante 8 horas) 22 ' 1 met = 50 kcal/m2.h (indice metabolico o calor emitido por una persona media en descanso por metro cuadrado de piel en una hora) 23 ' 1 met = 1 kcal/kg.h (relacion aproximada) 24 ' Peso medio de una persona 80 kg 25 ' Superficie aproximada media por persona 80/50=1,6 m2 26 ' Tasa Metabolica Basal media TMB= 1*80*24=1900 kcal al dia 27 ' Tasa Metabolica Total media TMT= 1,5*80*24=3000 kcal al dia 28 ' 1 W = 1 J/s = 1/4187 kcal/s = 3600/4187 kcal/h = 0,86 kcal/h 29 ' 1 met = 50/0,86 W/m2 = 58,15 W/m2 30 CLS:KEY OFF:PRINT"INTRODUCCIÓN DE DATOS" 31 PRINT"INTRODUCIR hr (%) O PRESIÓN DE VAPOR DE AGUA pa (Pa), PERO NO AMBAS" 41 PRINT"Vestimenta (clo) ...............................Icl = ? 1.0 (de 0 a 2)" 42 PRINT"Tasa Metabolica Total (met) ................... TMT = ? 1.2 (de 0.8 a 4)" 43 PRINT"Trabajo externo, normalmente alrededor de 0 (met) W = ? 0" 44 PRINT"Temperatura del aire (§C) .................... ta = ? 19 (de 10 a 30)" 45 PRINT"Temperatura radiante media (§C) ............... tr = ? 18 (de 10 a 40)" 46 PRINT"Velocidad relativa del aire (m/s) ............. var = ? 0.1 (de 0 a 1)" 48 PRINT"Humedad relativa (%) .......................... hr = ? 40 (de 0 a 100)" 49 PRINT"Presion del vapor de agua (Pa) ................ pa = ? (de 0 a 2700)" 50 LOCATE 3,1 55 INPUT"Vestimenta (clo) ...............................Icl = "; CLO:IF CLO=0 THEN CLO=1.0 56 LOCATE CSRLIN-1,56:PRINT CLO " " 60 INPUT"Tasa Metabolica Total (met) ................... TMT = "; MET:IF MET=0 THEN MET=1.2 61 LOCATE CSRLIN-1,56:PRINT MET " " 65 INPUT"Trabajo externo, normalmente alrededor de 0 (met) W = "; WME 66 LOCATE CSRLIN-1,56:PRINT WME " " 70 INPUT"Temperatura del aire (§C) .................... ta = "; TA:IF TA=0 THEN TA=19 71 LOCATE CSRLIN-1,56:PRINT TA " " 75 INPUT"Temperatura radiante media (§C) ............... tr = "; TR:IF TR=0 THEN TR=18 76 LOCATE CSRLIN-1,56:PRINT TR " " 80 INPUT"Velocidad relativa del aire (m/s) ............. var = "; VEL:IF VEL=0 THEN VEL=0.1 81 LOCATE CSRLIN-1,56:PRINT VEL " " 85 INPUT"Humedad relativa (%) .......................... hr = "; RH 86 LOCATE CSRLIN-1,56:PRINT RH " ":IF RH<>0 THEN 140 90 INPUT"Presion del vapor de agua (Pa) ................ pa = "; PA 95 LOCATE CSRLIN-1,56:PRINT PA " " 140 DEF FNPS(T) = EXP(16.6536-4030.183/(T+235)) :'+presión de vapor saturado, kPa a determinada temperatura T 150 IF PA=0 THEN PA=RH*10*FNPS(TA) ELSE GOTO 160 :'+presión de vapor de agua, Pa para humedad RH % 151 PRINT"Presion del vapor de agua (Pa) ................ pa = " FIX(PA) 160 ICL = .155 * CLO :'+aislamiento térmico de la ropa en mK/W 170 M = MET * 58.15 :'+tasa metabólica en W/m 180 W = WME * 58.15 :'+trabajo externo, en W/m 190 MW = M - W :'+produccion de calor interno en el cuerpo humano, en W/m 200 IF ICL <= .078 THEN FCL = 1 + 1.29 * ICL ELSE FCL = 1.05 + 0.645*ICL :'+factor de superficie de la ropa 201 PRINT"Factor de superficie de la ropa ............... fcl = " FIX(100*FCL)/100 210 HCF = 12.1*SQR(VEL) :'+coef. transf. calor por convección forzada 220 TAA = TA + 273 :'temperatura del aire en Kelvin 230 TRA = TR + 273 :'temperatura radiante media en Kelvin 240 '-----CALCULAR LA TEMPERATURA DE SUPERFICIE DE LA ROPA tcl POR ITERACIÓN---- 250 TCLA = TAA + (35.5-TA) / (3.5*(6.45*ICL+.1)) :'primera estimación de la temperatura superficial de la ropa 260 P1 = ICL * FCL :'etapa de cálculo 270 P2 = P1 * 3.96 :'etapa de cálculo 280 P3 = P1 * 100 :'etapa de cálculo 290 P4 = P1 * TAA :'etapa de cálculo 300 P5 = 308.7 - .028 * MW + P2 * (TRA/100)^4 310 XN = TLCA / 100 320 XF = XN 330 N = 0 :'N, número de iteraciones 340 EPS = .00015 :'criterio de fin de iteración 350 XF = (XF+XN)/2 360 HCN =2.38 * ABS(100 * XF-TAA)^0.25:'coeficiente de transf. de calor por convecc. natural 370 IF HCF>HCN THEN HC=HCF ELSE HC=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 :'temperatura superficial de la ropa 421 PRINT"Coeficiente transm. calor por convec. (W/m2.§K) hc = " FIX(10*HC)/10 422 PRINT"Temperatura superficie de la ropa (§C) ........ tcl = " FIX(10*TCL)/10 430 '-------------------------------- COMPONENTES DE LA PÉRDIDA DE CALOR ----------------------------------- 431 TS = .303 * EXP(-.036*M) + .028 :'+coef. trans. de la sensación térmica 432 PRINT"Coeficiente transmis. de la sensacion termica . Ts = " FIX(1000*TS)/1000 440 HL1 = 3.05 * .001*(5733-6.99*MW-PA) :'pérdida de calor por difusión a través de la piel 450 IF MW > 58.15 THEN HL2 = .42 * (MW-58.15) ELSE HL2 = 0 :'pérdida de calor por sudor (comodidad) 460 HL3 = 1.7 * .00001 * M * (5867-PA) :'+pérdida de calor latente por respiración 470 HL4 = .0014 * M * (34-TA) :'+pérdida de calor seco por la respiración 480 HL5 = 3.96 * 10^(-8) * FCL * ((TCL+273)^4 - (TRA)^4) :'pérdida de calor por radiación 490 HL6 = FCL * HC * (TCL-TA) 491 PRINT"Gasto ener. calor interno cuerpo humano (W/m2) Gi = " FIX(10*MW)/10 492 PRINT"Perdida calor por difusion traves piel (W/m2) . hl1 = " FIX(10*HL1)/10 493 PRINT" '' de calor por sudor (W/m2) ............. hl2 = " FIX(10*HL2)/10 494 PRINT" '' de calor latente en respiracion (W/m2) hl3 = " FIX(10*HL3)/10 495 PRINT" '' de calor seco por respiracion (W/m2) . hl4 = " FIX(10*HL4)/10 496 PRINT" '' de calor por radiacion (W/m2) ......... hl5 = " FIX(10*HL5)/10 497 PRINT" '' de calor por conveccion (W/m2) ........ hl6 = " FIX(10*HL6)/10:HL=HL1+HL2+HL3+HL4+HL5+HL6 498 PRINT"Total perdidas de calor (W/m2) ................ hl = " FIX(10*HL)/10:IF MW <= HL THEN 500 499 PRINT"Calentamiento corporal (W/m2) ................. Gc = " FIX(10*(MW-HL))/10:GOTO 520 500 PRINT"Enfriamiento corporal (W/m2) .................. Gf = " FIX(10*(MW-HL))/10 510 '--------------------------------CÁLCULO DE LOS ÍNDICES PMV Y PPD ----------------------------------- 520 PMV = TS * (MW - HL1 - HL2 - HL3 - HL4 - HL5 -HL6) :'+voto medio estimado 530 PPD = 100 - 95 * EXP(-.03353*PMV^4 - .2179*PMV^2) :'+porcentaje estimado de insatisfechos 540 GOTO 570 550 PMV = 999999! 560 PPD = 100 570 'PRINT "SALIDA" :'salida 580 PRINT"Voto Medio Estimado .......................... PMV = ";:PRINT USING " ##.##"; PMV;:PRINT" [de +3 a -3] 590 PRINT"Porcentaje Estimado de Insatisfechos .......... PPD = ";:PRINT USING "###.##"; PPD 600 INPUT "NUEVO CALCULO (S/N)"; R$ 610 IF (R$ = "S" OR R$ = "s") THEN RUN 620 END