10 CLS:KEY OFF
11 PRINT"CALCULO TEMPERATURA PERFIL ACERO SIN PROTECCION EN ESTIMACION INCENDIO
12 'PRINT"AL EI DEL EDIFICIO (Referencia: HEB240, 4 caras, EI=60, uf=0.5)

21 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 del perfil (m2/m) ........... Am = ? 1.38   (obligatorio)"
25 PRINT"Concavidad del perfil de acero (s/n) ......  C = ? S"
26 PRINT"Area expuesta efectiva 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
34 PRINT"Estabilidad estructural del edificio (min)  EI = ? 60"
36 PRINT"Coeficiente grado de utilizacion (Ef/Rf) .. uf = ? 0.5"
37 PRINT"Temperatura critica perfil acero. (低) ... Ofc =   585"
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
42 PRINT"Coefi. emisividad sup.acero\recinto norm. Eres =   0.7":ERES=0.7
43 PRINT"Coefi. transf.termica conv. norm. (W/m2.免) ac =   25":AC=25
44 PRINT"Coefi. regulador activacion incendio .....   r = ? 1     (1 < r < 3)"
45 PRINT"Tiempo de temperatura maxima (min) ...... tmax = ? 20    (10 < tmax < 30)"
40 PRINT"Incremento de tiempo de flujo calor (s) ... At =   5":AT=5
41 PRINT"Iteraciones incremento temperatura perfil..  n =   720
46 PRINT"Temperatura maxima gas incendio (低) ..... Ogt =   735"
47 PRINT"Temperatura maxima perfil acero (低) ....  Omt =   659";
48 'PRINT"Tiempo de temperatura maxima (min) .....  tmax =   90";

49 LOCATE 2,1
51 INPUT"Titulo del calculo (Perfil acero y caras)..  T = ";T$:IF T$="" THEN T$="HEB240 (4 caras)"
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 del perfil (m2) ............. 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=.5:OFC=585:PA=7850:CA=600:AT=5:N=720:ERES=0.7:AC=25:TMAX=20:OGT=735:OMT=659:TMAX=90:LOCATE 25,1:GOTO 240
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 efectiva 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 INPUT"Estabilidad estructural del edificio (min)  EI = ";EI:IF EI=0 THEN EI=60
64 LOCATE CSRLIN-1,51:PRINT EI "    ":INPUT"Coeficiente grado de utilizacion (Ef/Rf) .. uf = ";UF:IF UF=0 THEN UF=0.5
65 OFC=39.19*LOG((1/(0.9674*UF^3.833))-1)+483:LOCATE CSRLIN-1,51:PRINT UF "    ":PRINT"Temperatura critica perfil acero. (低)     Ofc =  " FIX(OFC)
66 PRINT:PRINT:PRINT:PRINT:INPUT"Coefi. regulador activacion incendio ......  r = ";R:IF R=0 THEN R=1
67 INPUT"Tiempo de temperatura maxima (min) ...... tmax = ";TMAX:IF TMAX=0 THEN TMAX=20
68 PRINT:PRINT"Iteraciones incremento temperatura perfil .  n =  " EI*60/AT:N=EI*60/AT

69 LOCATE 23,1
70 OPEN "TEMP.TXT" FOR OUTPUT AS 1
71 PRINT #1,"Titulo del calculo (Perfil acero y caras)..  T =  " T$
72 PRINT #1,"Canto del perfil de acero (m) .............  h = " H
73 PRINT #1,"Ancho del perfil de acero (m) .............  b = " B
74 PRINT #1,"Area expuesta del perfil (m2/m)  .......... Am = " AM
75 PRINT #1,"Concavidad del perfil de acero (s/n) ......  C =  " C$
76 PRINT #1,"Area expuesta efectiva del perfil (m2) .... Ab = " AB
77 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"
84 PRINT #1,"Estabilidad estructural del edificio (min)  EI = " EI
86 PRINT #1,"Coeficiente grado de utilizacion (Ef/Rf) .. uf = " UF
87 PRINT #1,"Temperatura critica perfil acero. (慢) ... Ofc = " FIX(OFC)
88 PRINT #1,"Densidad del perfil acero (kg/m3) ........  Pa =  7850"
89 PRINT #1,"Calor especifico del perfil acero (J/kg.摘) Ca =  600"
91 PRINT #1,"Coefi. emisividad sup.acero\recinto norm. Eres =  0.7
92 PRINT #1,"Coefi. transf.termica conv. norm. (W/m2.摘) ac =  25"
93 PRINT #1,"Coefi. regulador activacion incendio ......  r = " R
94 PRINT #1,"Tiempo de temperatura maxima (min) .....  tmax = " TMAX
95 PRINT #1,"Incremento de tiempo de flujo calor (s) ..  At =  5"
96 PRINT #1,"Iteraciones incremento temperatura perfil .  n = " N

101 PRINT #1,"        I        TH       OMT       OGT     HNETR    HNETC     HNETD        AOT   OMT+AOT"
102 OMT=O1:FOR I=1 TO N
103 TH=(I*AT/3600)*R
104 'OGT=20+345*LOG(8*(TS/60)+1)/LOG(10)
105 IF TH<=TMAX/60 THEN E1=0.324*2.71^(-0.2*TH):E2=0.204*2.71^(-1.7*TH):E3=0.472*2.71^(-19*TH):OGT=20+1235*(1-E1-E2-E3):OGTM=OGT
106 IF TH>TMAX/60 THEN OGT=OGTM-625*(TH-(TMAX/60))

107 HNETR=5.67*(10^(-8))*ERES*(((OGT+273)^4)-(OMT+273)^4)
108 HNETC=AC*(OGT-OMT)
109 HNETD=HNETR+HNETC

115 AOT=SM*KSH*HNETD*AT/(PA*CA):PRINT #1,USING"######.## ";I;TH;OMT;OGT;HNETR;HNETC;HNETD;AOT;OMT+AOT:OMT=OMT+AOT
116 IF AOT>=0 THEN OMTM=OMT
117 IF OMT<=20 THEN THM=TH:GOTO 131
120 NEXT I

131 PRINT"Temperatura maxima gas incendio (低) ..... Ogt =  " FIX(OGTM)
133 PRINT"Temperatura maxima perfil acero (低) ..... Omt =  " FIX(OMTM)
134 IF THM>0 THEN 135 ELSE THM=(20/60)+(-(20/60)+OGTM)/625
135 PRINT"Tiempo total duracion incendio (min) .....  tt =  " FIX(THM*60)

141 PRINT #1,"Temperatura maxima gas incendio (慢) ..... Ogt = " FIX(OGTM)
143 PRINT #1,"Temperatura maxima perfil acero (慢) ..... Omt = " FIX(OMTM)
145 PRINT #1,"Tiempo total duracion incendio (min) .....  tt = " FIX(THM*60)

240 INPUT;"Pulsar tecla para terminar ";Z
300 SYSTEM