10 CLS:KEY OFF
11 PRINT"CALCULO TEMPERATURA DE PERFIL ACERO CON PROTECCION EN COMPROBACION ESTABILIDAD
12 'PRINT"AL FUEGO PARA EI DE EDIFICIO (Referencia: HEB240, 4 caras, EI=60, uf=0.5, 
13 'PRINT"Placa de yeso)
18 PRINT"Titulo del calculo (Perfil acero y caras).   T = ? HEB240 (4 caras)"
19 PRINT"Material protector (denominacion) .........  P = ? Placa de yeso"  
20 PRINT"Canto del perfil de acero (m) .............  h = ? 0.24  (indicativo)"
21 PRINT"Ancho del perfil de acero (m) .............  b = ? 0.24  (indicativo)"
22 PRINT"Area de la proteccion del perfil (m2/m) ..  Ap = ? 0.96  (obligatorio)"
23 PRINT"Volumen de acero del perfil (m3/m) .......   V = ? 0.0106"
24 PRINT"Factor de seccion (m-1) ..................  Sp =   90.5"
25 PRINT"Temperatura inicial perfil acero. (低) ...  O1 =   20":O1=20
29 PRINT"Densidad del perfil acero (kg/m3) ........  Pa =   7850":PA=7850
30 PRINT"Calor especifico del perfil acero (J/kg.免) Ca =   600":CA=600
31 PRINT"Densidad caracteris. del protector (kg/m3) Ppk = ? 800"
32 PRINT"Calor espe. carac. del protector (J/kg.免) Cpk = ? 1700"
34 PRINT"Conductividad cara. del protector (W/m.免) lpk = ? 0.20"
35 PRINT"Coeficiente parcial seguridad del protector Yp = ? 1.00  (1.00\1.10\1.25)"
36 PRINT"Espesor del material protector (m) .......  dp = ? 0.015"
27 PRINT"Coeficiente grado de utilizacion (Ef/Rf) .  uf = ? 0.5"
28 PRINT"Temperatura critica perfil acero. (低) ... Ofc =   585"
26 PRINT"Estabilidad estructural del edificio (min)  EI = ? 60"

37 PRINT"Incremento de tiempo de flujo calor (s) ..  At =   30":AT=30
38 PRINT"Iteraciones incremento temperatura perfil.   n =   120":N=120
40 PRINT"Resistividad ter.ef. protector (m2.免/W) rpefd =   0.084"
41 PRINT"Temperatura final gas incendio (低) ...... Ogt =   945"
42 PRINT"Temperatura final perfil acero (低) ...... Omt =   472";

49 LOCATE 2,1
50 INPUT"Titulo del calculo (Perfil acero y caras)..  T = ";T$:IF T$="" THEN T$="HEB240 (4 caras)"
51 IF T$="1" THEN GOSUB 200:GOTO 69

52 INPUT"Material protector (denominacion) .........  P = ";P$:IF P$="" THEN P$="Placa de yeso"
53 INPUT"Canto del perfil de acero (m) .............  h = ";H:IF H=0 THEN H=0.24
54 LOCATE CSRLIN-1,51:PRINT H "                 ":INPUT"Ancho del perfil de acero (m) .............  b = ";B:IF B=0 THEN B=0.24
55 LOCATE CSRLIN-1,51:PRINT B "                 ":INPUT"Area de la proteccion del perfil (m2/m) ..  Ap = ";AP:IF AP=0 THEN AP=0.96:V=0.0106:SM=90.5:O1=20:EI=60:UF=.5:OFC=585:PA=7850:CA=600:PPK=800:CPK=1700:LPK=0.20:YP=1.00:DP=0.015:AT=30:N=120:RPEFD=0.084:AC=OGT=945:OMT=472:LOCATE 23,1:GOTO 150
56 LOCATE CSRLIN-1,51:PRINT AP "                  ":INPUT"Volumen de acero del perfil (m3/m) .......   V = ";V:IF V=0 THEN V=0.0106
57 LOCATE CSRLIN-1,51:PRINT V "     ":PRINT"Factor de seccion (m-1) ..................  Sp =  " FIX(10*AP/V)/10 "   ":SP=AP/V:PRINT

60 PRINT:PRINT:INPUT"Densidad caracteris. del protector (kg/m3) Ppk = ";PPK:IF PPK=0 THEN PPK=800
61 LOCATE CSRLIN-1,51:PRINT PPK "    ":INPUT"Calor espe. carac. del protector (J/kg.免) Cpk = ";CPK:IF CPK=0 THEN CPK=1700
62 LOCATE CSRLIN-1,51:PRINT CPK "    ":INPUT"Conductividad cara. del protector (W/m.免) lpk = ";LPK:IF LPK=0 THEN LPK=0.20
63 LOCATE CSRLIN-1,51:PRINT LPK "    ":INPUT"Coeficiente parcial seguridad del protector Yp = ";YP:IF YP=0 THEN YP=1.00
64 LOCATE CSRLIN-1,51:PRINT YP "    ":INPUT"Espesor del material protector (m) .......  dp = ";DP:IF DP=0 THEN DP=0.015
65 LOCATE CSRLIN-1,51:PRINT DP "     ":INPUT"Coeficiente grado de utilizacion (Ef/Rf) .  uf = ";UF:IF UF=0 THEN UF=0.5
66 LOCATE CSRLIN-1,51:PRINT UF "    ":OFC=39.19*LOG((1/(0.9674*UF^3.833))-1)+483:PRINT"Temperatura critica perfil acero. (低) ... Ofc =  " FIX(OFC)
67 INPUT"Estabilidad estructural del edificio (min)  EI = ";EI:IF EI=0 THEN EI=60
68 LOCATE CSRLIN-1,51:PRINT EI "    ":PRINT:PRINT"Iteraciones incremento temperatura perfil..  n =  " EI*60/AT:N=EI*60/AT

69 LOCATE 22,1

70 OPEN "TEMP.TXT" FOR OUTPUT AS 1
71 PRINT #1,"CALCULO TEMPERATURA DE PERFIL ACERO CON PROTECCION EN COMPROBACION ESTABILIDAD AL FUEGO PARA EI DE EDIFICIO (EAE)
72 PRINT #1, D$
73 PRINT #1,"Titulo del calculo (Perfil acero y caras)..  T =  " T$
74 PRINT #1,"Material protector (denominacion) .........  P =  " P$ 

75 PRINT #1,"Canto del perfil de acero (m) .............  h = " H
76 PRINT #1,"Ancho del perfil de acero (m) .............  b = " B
77 PRINT #1,"Area de la proteccion del perfil (m2/m) ... Ap = " AP
78 PRINT #1,"Volumen de acero del perfil (m3/m) ........  V = " V
79 PRINT #1,"Factor de seccion (m-1) ................... Sp = " FIX(10*SP)/10
82 PRINT #1,"Temperatura inicial perfil acero. (慢) .... O1 =  20"
86 PRINT #1,"Densidad del perfil acero (kg/m3) ......... Pa =  7850"
87 PRINT #1,"Calor especifico del perfil acero (J/kg.摘) Ca =  600"
88 PRINT #1,"Densidad caracteris. del protector (kg/m3) Ppk = " PPK
90 PRINT #1,"Calor espe. carac. del protector (J/kg.摘) Cpk = " CPK
91 PRINT #1,"Conductividad cara. del protector (W/m.摘) lpk = " LPK
92 PRINT #1,"Coeficiente parcial seguridad del protector Yp = " YP
93 PRINT #1,"Espesor del material protector (m) ........ dp = " DP
94 PRINT #1,"Coeficiente grado de utilizacion (Ef/Rf) .. uf = " UF
95 PRINT #1,"Temperatura critica perfil acero. (慢) ... Ofc = " FIX(OFC)
96 PRINT #1,"Estabilidad estructural del edificio (min). EI = " EI
98 PRINT #1,"Incremento de tiempo de flujo calor (s) ... At =  30"
99 PRINT #1,"Iteraciones incremento temperatura perfil..  n =  120"

102 CPD=CPK:PPD=PPK:F=PPD*CPD*DP*SP/(PA*CA):RPEFK=(1+F/3)*DP/LPK:RPEFD=RPEFK/YP:FP=(2.71^(F/10))-1
103 OAT=O1:OGT0=20
104 PRINT #1,"Resistividad ter.ef. protector (m2.摘/W) rpefd = " FIX(1000*RPEFD)/1000
105 PRINT #1,"Coeficiente reduc.temp. por cap.calorifica  fp = " FIX(10000*FP)/10000
106 PRINT #1,"        I        TS       OAT       OGT       AOT   OMT+AOT"

107 FOR I=1 TO N
108 TS=I*AT
109 OGT=20+345*LOG(8*(TS/60)+1)/LOG(10):AOGT=OGT-OGT0:OGT0=OGT

110 AOAT=((SP/RPEFD)*(OGT-OAT)*AT)/((PA*CA)-((2.71^(F/10))-1)*AOGT)
115 PRINT #1,USING"######.## ";I;TS;OAT;OGT;AOAT;OAT+AOAT:OAT=OAT+AOAT
120 NEXT I

130 PRINT"Resistividad ter.ef. protector (m2.免/W) rpefd =  " FIX(1000*RPEFD)/1000 " "
131 PRINT"Temperatura final gas incendio (低)....... Ogt =  " FIX(OGT)
133 PRINT"Temperatura final perfil acero (低)....... Oat =  " FIX(OAT)

141 PRINT #1,"Temperatura final gas incendio (慢)....... Ogt = " FIX(OGT)
143 PRINT #1,"Temperatura final perfil acero (慢)....... Oat = " FIX(OAT)

150 INPUT;"Pulsar tecla para terminar ";Z
160 SYSTEM

200 OPEN "DATOS.TXT" FOR INPUT AS 1:LINE INPUT #1, D$:CLOSE
201 OPEN "DATOS.TXT" FOR INPUT AS 1
210 INPUT #1, T$,P$,H,B,AP,V,PPK,CPK,LPK,YP,DP,UF,EI
215 CLOSE
230 SP=AP/V
240 O1=20
250 PA=7850
255 CA=600
260 OFC=39.19*LOG((1/(0.9674*UF^3.833))-1)+483
265 AT=30
266 N=EI*60/AT
269 LOCATE CSRLIN-1,1
270 PRINT"Titulo del calculo (Perfil acero y caras)..  T =   " T$
271 PRINT"Material protector (denominacion) .........  P =   " P$
272 PRINT"Canto del perfil de acero (m)  ............  h =  " H "  "
273 PRINT"Ancho del perfil de acero (m)  ............  b =  " B "  "
274 PRINT"Area de la proteccion del perfil (m2/m) ... Ap =  " AP "  "
281 PRINT"Volumen de acero del perfil (m3/m) ........  V =  " V "  "
282 PRINT"Factor de seccion (m-1) ................... Sp =  " FIX(10*SP)/10 "  "
283 PRINT:PRINT:PRINT
290 PRINT"Densidad caracteris. del protector (kg/m3) Ppk =  " PPK "  "
291 PRINT"Calor espe. carac. del protector (J/kg.免) Cpk =  " CPK "  "
292 PRINT"Conductividad cara. del protector (W/m.免) lpk =  " LPK "  "
293 PRINT"Coeficiente parcial seguridad del protector Yp =  " YP "  "
294 PRINT"Espesor del material protector (m) .......  dp =  " DP "   "
295 PRINT"Coeficiente grado de utilizacion (Ef/Rf) .  uf =  " uf " "
296 PRINT"Temperatura critica perfil acero. (低).... Ofc =  " FIX(OFC) "  "
297 PRINT"Estabilidad estructural del edificio (min). EI =  " EI "  "
298 PRINT
299 PRINT"Iteraciones incremento temperatura perfil..  n =  " N " "

300 RETURN