martes, 12 de febrero de 2019

Hormigón

Programa básico para resolver problemas de dimensionamiento de una viga de hormigón armado. Se especifican las unidades de entrada y de salida de los datos

1 PRINT"HORMIGON"
2 INPUT"INTRODUCE DIMENSIONES: BASE(M)=";B,"ALTURA(M)=";H
3 INPUT"FCK";FCK:EC=8500*((FCK+8)^(1/3)):PRINT"EC(MPa)=";OC=0.6*FCK:PRINT"OC(max) (MPa)=";OC:INPUT"MODULO ACERO (GPa)=";ES
4 INPUT"FISURACION Y FLECHAS (1) DIMENSIONAR (2)";H11
5 IF H11=1 GOTO7
6 IF H11=2 GOTO90
7 PRINT"ARMADURA MAYOR":INPUT"DIAMETRO(mm)"; FI1,"NUMERO DE ARM";NR1:A1=NR1*PI*FI1^2/4:PRINT"AS(TRACCION)'(mm2)=";A1
8 PRINT"ARMADURA MENOR":INPUT"DIAMETRO(mm)"; FI,"NUMERO DE ARM";NR:A=NR*PI*FI^2/4:PRINT"AS(COMP)'(mm2)=";A
9 INPUT"RECUBRIMIENTO MECANICO (m)";RCM
10 N=ES*1E3/EC:PRINT "COCIENTE A CERO/HORM=";N
*/Nota, aqui despejamos para obtener n y que cumpla la fisuración/*
11 INPUT"INTRODUCE BETA (1.7)=";BE,"K1(.125)=";K1,"K2(.5)=";K2
14 FCM=0.3*CUR(FCK^2):PRINT"FCT,M(MPa)=";FCM
16 FC1=1.6-H/1000*FCM:IF FC1>FCM THEN FCF=FC1 ELSE FCF=FCM:PRINT"FCTM,FL (MPa)=";FCF
17 IB=1/12*B*H^3:PRINT"IB (m4)=";IB
18 MCR=FCF*IB/(H/2)*1E3:PRINT"MCR(KNm)=";MCR
19 PRINT"MCR>MA FISURACION"
20 PRINT"EMPOTRAMIENT"
21 C1=RCM*1E3-FI1/2:PRINT"c(mm)=";C1
22 AE1=B*1E3*(RCM*1E3+7.5*FI1):PRINT"AEFICAZ EMPOTR(mm2)=";AE1
23 INPUT"HUECOS ENTRE ARMADURAS";P1:S1=(B-2*RCM)/P1*1E3:PRINT"SEPARACION(mm)=";S1
24 SM1=2*C1+.2*S1+.4*K1*FI1*AE1/A1:PRINT"SM EMPOTRAMIENTO(mm)=";SM1
25 PRINT"SECCION CENTRO LUZ"
26 C=RCM*1E3-FI/2:PRINT"c(mm)=";C
27 AE=B*1E3*(RCM*1E3+7.5*FI):PRINT"AEFICAZ C LUZ(mm2)=";AE
28 INPUT"HUECOS ENTRE ARMADURAS";P:S=(B-2*RCM)/P*1E3:PRINT"SEPARACION(mm)=";S
29 SM=2*C+.2*S+.4*K1*FI*AE/A:PRINT"SM EMPOTRAMIENTO(mm)=";SM
*/Nota: ecuacion de segundo grado/*
30 PRINT"ECUACION":D=H-RCM:D1=RCM:PRINT"D(C)=";D."D'(T)=";D1
31 c=-(A1*D1+A*D)*N
32 c1=-(A*D1+A1*D)*N
33 b=(N*(A1+A)/100):a=1/2*B/100
34 x=((-b+sqr(b^2-4*a*c))/2/a)*10
36 x1=((-b+sqr(b^2-4*a*c1))/2/a)*10
37 PRINT"EMPOTRAMIENTO", "AS'(mm)=";A,"AS(mm)=";A1,"d'(m)=";D1,"d(m)=";D
38 PRINT"X(cm)=";x1/10
39 IG1=(1/3*B*1E3*x1^3+N*A*(x1-D1*1E3)^2+N*A1*(x1-D*1E3)^2
40 PRINT"IGH(m4)=";IG1;IG1*1E-12
41 OS1=N*(x1*1E-3-(H-RCM))/(IG1*1E-12)/1000
*/Nota: Tensor en funcion del momento/*
42 PRINT"OS[MA] (MPa)=",OS1
43 INPUT"MOMENTO DE EMPOTRAMIENTO";M1
44 IF M1<-MCR THEN PRINT"FISURA" ELSE PRINT"NO FISURA"
45 OS1=OS1*M1:PRINT"OS=";OS1
45 ESM1I=(OS1/(ES*1E3))*(1-K2*(MCR1/M1)^2):PRINT"ESM(LADOIZQ ECUACION)=";ESM1I