10 CLS:KEY OFF
11 PRINT"CALCULO TEMPERATURA PERFIL ACERO SIN PROTECCION EN COMPROBACION ESTABILIDAD
12 PRINT"AL FUEGO PARA EI DE EDIFICIO (Referencia: HEB240, 4 caras, uf=0.5, EI=60)

20 PRINT"Titulo del calculo (Perfil acero y caras)..  T = ? HEB240 (4 caras)"
22 PRINT"Canto del perfil de acero (m) .............  h = ? 0.24   (indicativo)"
23 PRINT"Ancho del perfil de acero (m) .............  b = ? 0.24   (indicativo)"
24 PRINT"Area expuesta total del perfil (m2/m) ..... Am = ? 1.38   (obligatorio)"
25 PRINT"Concavidad del perfil de acero (s/n) ......  C = ? S"
26 PRINT"Area expuesta envolvente del perfil (m2/m)  Ab = ? 0.96"
27 PRINT"Factor reductor por concavidad ........... Ksh =   0.62"
31 PRINT"Volumen de acero del perfil (m3/m) ........  V = ? 0.0106"
32 PRINT"Factor de seccion (m-1) ................... Sm =   130"
33 PRINT"Temperatura inicial perfil acero. (低) .... O1 =   20":O1=20
38 PRINT"Densidad del perfil acero (kg/m3).......... Pa =   7850":PA=7850
39 PRINT"Calor especifico del perfil acero (J/kg.免) Ca =   600":CA=600
40 PRINT"Coef. emisividad sup.perfil\recinto norm. Eres = ? 0.7"
41 PRINT"Coefi. transf.termica conv. norm. (W/m2.免) ac = ? 25"

42 PRINT"Coeficiente grado de utilizacion (Ef/Rf) .. uf = ? 0.5"
43 PRINT"Temperatura critica perfil acero. (低).... Ofc =   585"
44 PRINT"Estabilidad estructural del edificio (min). EI = ? 60"
45 PRINT"Incremento de tiempo de flujo calor (s).... At =   5":AT=5
46 PRINT"Iteraciones incremento temperatura perfil..  n =   720
47 PRINT"Temperatura final gas incendio (低)....... Ogt =   945"
48 PRINT"Temperatura final perfil acero (低)....... Omt =   937"

49 LOCATE 3,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"Canto del perfil de acero (m) .............  h = ";H:IF H=0 THEN H=0.24
53 LOCATE CSRLIN-1,51:PRINT H "                 ":INPUT"Ancho del perfil de acero (m) .............  b = ";B:IF B=0 THEN B=0.24
54 LOCATE CSRLIN-1,51:PRINT B "                 ":INPUT"Area expuesta total del perfil (m2/m) ..... Am = ";AM:IF AM=0 THEN AM=1.38:C$="S":AB=0.96:KSH=0.62:V=0.0106:SM=130:O1=20:EI=60:UF=0.5:OFC=585:PA=7850:CA=600:AT=5:N=720:ERES=0.7:AC=25:OGT=945:OFT=937:LOCATE 24,1:GOTO 150
55 LOCATE CSRLIN-1,51:PRINT AM "                  ":INPUT"Concavidad del perfil de acero (s/n) ......  C = ";C$:IF C$="S" OR C$="s" THEN C$="S":GOTO 56 ELSE C$="N":AB=AM:KSH=1
56 LOCATE CSRLIN-1,52:PRINT C$ "    ":IF AM<>AB THEN INPUT"Area expuesta envolvente del perfil (m2/m)  Ab = ";AB ELSE LOCATE,51:PRINT AB "     ":PRINT"Factor reductor por concavidad    .......  Ksh =   1.00      ":GOTO 61
57 IF AB=0 THEN AB=0.96:KSH=0.9*AB/AM ELSE KSH=0.9*AB/AM
58 LOCATE CSRLIN-1,51:PRINT AB "    ":PRINT"Factor reductor por concavidad ........... Ksh =  " FIX(100*KSH)/100 "   "
61 INPUT"Volumen de acero del perfil (m3/m).........  V = ";V:IF V=0 THEN V=0.0106
62 LOCATE CSRLIN-1,51:PRINT V "    ":PRINT"Factor de seccion (m-1) ................... Sm =  " FIX(10*AM/V)/10 "   ":SM=AM/V:PRINT

63 PRINT:PRINT:INPUT"Coef. emisividad sup.perfil\recinto norm. Eres = ";ERES:IF ERES=0 THEN ERES=0.7
64 LOCATE CSRLIN-1,51:PRINT ERES "    ":INPUT"Coefi. transf.termica conv. norm. (W/m2.免) ac = ";AC:IF AC=0 THEN AC=25
65 LOCATE CSRLIN-1,51:PRINT AC "    ":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 PERFIL ACERO SIN 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,"Canto del perfil de acero (m)  ............  h = " H
75 PRINT #1,"Ancho del perfil de acero (m)  ............  b = " B
76 PRINT #1,"Area expuesta total del perfil (m2/m) ..... Am = " AM
77 PRINT #1,"Concavidad del perfil de acero (s/n) ......  C =  " C$
78 PRINT #1,"Area expuesta envolvente del perfil (m2/m)  Ab = " AB
79 PRINT #1,"Factor reductor por concavidad    ........ Ksh = " FIX(100*KSH)/100
81 PRINT #1,"Volumen de acero del perfil (m3/m) ........  V = " V
82 PRINT #1,"Factor de seccion (m-1) ................... Sm = " FIX(10*AM/V)/10
83 PRINT #1,"Temperatura inicial perfil acero. (慢)..... O1 =  20"
88 PRINT #1,"Densidad del perfil acero (kg/m3) ........  Pa =  7850"
89 PRINT #1,"Calor especifico del perfil acero (J/kg.摘) Ca =  600"
92 PRINT #1,"Coefi. emisividad sup.acero\recinto norm. Eres = " ERES
93 PRINT #1,"Coefi. transf.termica conv. norm. (W/m2.摘) ac = " AC
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
97 PRINT #1,"Incremento de tiempo de flujo calor (s)...  At =  5"
98 PRINT #1,"Iteraciones incremento temperatura perfil..  n = " N

101 PRINT #1,"        I        TS       OMT       OGT     HNETR    HNETC     HNETD        AOT   OMT+AOT"
102 OMT=O1:FOR I=1 TO N
103 TS=I*AT

105 OGT=20+345*LOG(8*(TS/60)+1)/LOG(10)
106 HNETR=5.67*(10^(-8))*ERES*(((OGT+273)^4)-(OMT+273)^4)
107 HNETC=AC*(OGT-OMT)
108 HNETD=HNETR+HNETC

110 AOT=SM*KSH*HNETD*AT/(PA*CA):PRINT #1,USING"######.## ";I;TS;OMT;OGT;HNETR;HNETC;HNETD;AOT;OMT+AOT
111 OMT=OMT+AOT
120 NEXT I

131 PRINT"Temperatura final gas incendio (低)....... Ogt =  " FIX(OGT)
133 PRINT"Temperatura final perfil acero (低)....... Omt =  " FIX(OMT)

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

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$,H,B,AM,C$,AB,V,ERES,AC,UF,EI
215 CLOSE
220 IF C$="S" OR C$="s" THEN C$="S":GOTO 221 ELSE C$="N":AB=AM:KSH=1
221 IF AM<>AB THEN KSH=0.9*AB/AM
230 SM=AM/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=5
266 N=EI*60/AT
270 LOCATE CSRLIN-1,1
271 PRINT"Titulo del calculo (Perfil acero y caras)..  T =   " T$
272 PRINT"Canto del perfil de acero (m)  ............  h =  " H
273 PRINT"Ancho del perfil de acero (m)  ............  b =  " B
274 PRINT"Area expuesta total del perfil (m2/m) ..... Am =  " AM
275 PRINT"Concavidad del perfil de acero (s/n) ......  C =   " C$
276 PRINT"Area expuesta envolvente del perfil (m2/m)  Ab =  " AB
277 PRINT"Factor reductor por concavidad    ........ Ksh =  " FIX(100*KSH)/100
281 PRINT"Volumen de acero del perfil (m3/m) ........  V =  " V
282 PRINT"Factor de seccion (m-1) ................... Sm =  " FIX(10*AM/V)/10
283 PRINT"Temperatura inicial perfil acero. (低)..... O1 =   20"
288 PRINT"Densidad del perfil acero (kg/m3) ........  Pa =   7850"
289 PRINT"Calor especifico del perfil acero (J/kg.免) Ca =   600"
291 PRINT"Coefi. emisividad sup.acero\recinto norm. Eres =  " ERES
293 PRINT"Coefi. transf.termica conv. norm. (W/m2.免) ac =  " AC
294 PRINT"Coeficiente grado de utilizacion (Ef/Rf) .. uf =  " UF
295 PRINT"Temperatura critica perfil acero. (低).... Ofc =  " FIX(OFC)
296 PRINT"Estabilidad estructural del edificio (min). EI =  " EI
297 PRINT"Incremento de tiempo de flujo calor (s)...  At =   5"
298 PRINT"Iteraciones incremento temperatura perfil..  n =  " N

300 RETURN