abaco seção simetrica

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clear
clc
//dados do problema
fc=12.0000 //fc=0,85fck/gamac
fyd=434.0000 // fyd = fyk/gamas
Es=210000.0000 //Módulo de Elasticidade do aço CA-50
deformacaocu=0.0035 //Deformação máxima do concreto
deformacaoyd=fyd/Es //Deformação máxima do aço
alfa=0.0700 //Eq II.4
//Inicialização da tolerância e da taxa mecânica de armadura desejados
tole= input ("Digite a tolerância desejada: ")
w= input ("Digite qual valor de w você deseja calculcar: ")
w= w/2.0000 //O w utilizado é só de um dos bordos, e como o cálculo é feito para armadura simétrica, devemos dividi-lo por 2.
for i=1:20 //Laço para variar o v, que será o eixo x do ábaco.
 v=0.1000*i; //varia de 0.1 a 2.0 para a construção do ábaco 2. 
 v1=0.0000; 
 dif=10.0000; 
 j=0.0000;
 while dif>tole // laço para testar qual beta resultará o valor mais próximo, em relação a tolerância, do v em questão.
 j=j+1.0000;
 bet=alfa+0.0001*j;
 deformacaosl=(1.0000-(0.8000*alfa/bet))*deformacaocu; //Eq. II.41
 deformacaos=((0.8000*(1.0000-alfa)/bet)-1.0000)*deformacaocu; //Eq. II.42
 if deformacaosl<deformacaoyd then //Laço pra descobrir se a armadura do bordo "superior" está no patamar de escoamento ou na fase elástica
 tsl=deformacaosl*Es; //Eq. II.43
 else
 tsl=fyd;
 end
 if deformacaos<deformacaoyd then //Laço pra descobrir se a armadura do bordo "inferior" está no patamar de escoamento ou na fase elástica
 ts=deformacaos*Es; //Eq.II.44
 else
 ts=fyd;
 end
 vl=bet+(w*(tsl-ts))/fyd; //Eq. II.56
 dif=abs(vl-v);
 end
 mi=(bet*(1.0000-bet)+((w*(1.0000-2.0000*alfa))*(tsl+ts))/fyd)/2.0000 //Eq. II.66
 disp(mi, "mi=")
end