hola. ahora no te entiendo lo q me dices. te voy a mandar el programa q yo hice. lo q yo quiero hacer es q despues de q me rechace el valor me vuelva a pedir la temperatura hasta q este correcta. si tienes tiempo dime como hacerlo. muchisisisisisisisisimas gracias.
el programa es este:
clc
disp('==============================================')
disp(' ')
disp(' PARCIAL DE DISEÑO DE PLANTAS I ');
disp(' PRODUCCIÓN DE CLORURO DE ALILO A PARTIR DE PROPILENO Y CLORO')
disp(' ')
disp('==============================================')
disp(' ')
p = input('producción de cloruro de alilo al 99.5%wt (Kg/año) = ');
rr=5;
T = input('Temperatura entrada al reactor (ºC) = ');
if T>=440 & T<=520
s1=-3.3333e-6;
s2=0.0054214;
s3=-2.9902;
s4=563.28;
s=(s1*(T^3))+(s2*(T^2))+(s3*T)+s4;
r1=-9.375e-8;
r2=0.00015786;
r3=-0.090105;
r4=17.56;
r=(r1*(T^3))+(r2*(T^2))+(r3*T)+r4;
c1=-8.7261e-20;
c2=-3.2143e-5;
c3=0.036357;
c4=-9.224;
c=(c1*(T^3))+(c2*(T^2))+(c3*T)+c4;
allyl6=((0.995*p)*(100/99))/76.45;
dicloro6=allyl6/((r+1)*s);
doscloro6=dicloro6*r;
hcl6=allyl6+doscloro6+2*dicloro6;
cloro5=hcl6/c;
propileno5=rr*cloro5;
cloro6=cloro5-hcl6;
propileno6=propileno5-allyl6-doscloro6-dicloro6;
propileno8=propileno6;
propileno4=propileno5;
propileno3=propileno4;
propileno1=propileno3-propileno8;
cloro12=cloro6;
cloro4=cloro5;
cloro2=cloro4-cloro12;
hcl9=hcl6;
Hcl9=36.45*hcl9;
corriente9=Hcl9*(100/32);
h2o=0.68*corriente9;
Ganancia_con_recirculado=(1.72*p)+(0.075*corriente9)-(0.20*(dicloro6*110.9))-(0.2*((76.45*doscloro6)-(0.005*p)))-(0.0003*h2o)-(0.47*(propileno1*42))-(0.18*(cloro2*70.9))
plot(T,Ganancia_con_recirculado,'o'), hold on
s1=-3.3333e-6;
s2=0.0054214;
s3=-2.9902;
s4=563.28;
s=(s1*(T^3))+(s2*(T^2))+(s3*T)+s4;
r1=-9.375e-8;
r2=0.00015786;
r3=-0.090105;
r4=17.56;
r=(r1*(T^3))+(r2*(T^2))+(r3*T)+r4;
c1=-8.7261e-20;
c2=-3.2143e-5;
c3=0.036357;
c4=-9.224;
c=(c1*(T^3))+(c2*(T^2))+(c3*T)+c4;
allyl6=((0.995*p)*(100/99))/76.45;
dicloro6=allyl6/((r+1)*s);
doscloro6=dicloro6*r;
hcl6=allyl6+doscloro6+2*dicloro6;
cloro5=hcl6/c;
propileno5=rr*cloro5;
cloro6=cloro5-hcl6;
propileno6=propileno5-allyl6-doscloro6-dicloro6;
propileno8=propileno6;
propileno4=propileno5;
propileno3=propileno4;
propileno1=propileno3-propileno8;
cloro4=cloro5;
cloro2=cloro4;
hcl9=hcl6+(2*cloro6);
Hcl9=36.45*hcl9;
corriente9=Hcl9*(100/32);
h2o=(0.68*corriente9)+(18*cloro6);
Ganancia_reaccion=(1.72*p)+(0.075*corriente9)-(0.20*(dicloro6*110.9))-(0.2*((76.45*doscloro6)-(0.005*p)))-(0.0003*h2o)-(0.47*(propileno1*42))-(0.18*(cloro2*70.9))
plot(T,Ganancia,'o'), hold on
for T=440:0.5:520
s1=-3.3333e-6;
s2=0.0054214;
s3=-2.9902;
s4=563.28;
s=(s1*(T^3))+(s2*(T^2))+(s3*T)+s4;
r1=-9.375e-8;
r2=0.00015786;
r3=-0.090105;
r4=17.56;
r=(r1*(T^3))+(r2*(T^2))+(r3*T)+r4;
c1=-8.7261e-20;
c2=-3.2143e-5;
c3=0.036357;
c4=-9.224;
c=(c1*(T^3))+(c2*(T^2))+(c3*T)+c4;
allyl6=((0.995*p)*(100/99))/76.45;
dicloro6=allyl6/((r+1)*s);
doscloro6=dicloro6*r;
hcl6=allyl6+doscloro6+2*dicloro6;
cloro5=hcl6/c;
propileno5=rr*cloro5;
cloro6=cloro5-hcl6;
propileno6=propileno5-allyl6-doscloro6-dicloro6;
propileno8=propileno6;
propileno4=propileno5;
propileno3=propileno4;
propileno1=propileno3-propileno8;
cloro12=cloro6;
cloro4=cloro5;
cloro2=cloro4-cloro12;
hcl9=hcl6;
Hcl9=36.45*hcl9;
corriente9=Hcl9*(100/32);
h2o=0.68*corriente9;
Ganancia=(1.72*p)+(0.075*corriente9)-(0.20*(dicloro6*110.9))-(0.2*((76.45*doscloro6)-(0.005*p)))-(0.0003*h2o)-(0.47*(propileno1*42))-(0.18*(cloro2*70.9));
plot(T,Ganancia,'k'), hold on
s1=-3.3333e-6;
s2=0.0054214;
s3=-2.9902;
s4=563.28;
s=(s1*(T^3))+(s2*(T^2))+(s3*T)+s4;
r1=-9.375e-8;
r2=0.00015786;
r3=-0.090105;
r4=17.56;
r=(r1*(T^3))+(r2*(T^2))+(r3*T)+r4;
c1=-8.7261e-20;
c2=-3.2143e-5;
c3=0.036357;
c4=-9.224;
c=(c1*(T^3))+(c2*(T^2))+(c3*T)+c4;
allyl6=((0.995*p)*(100/99))/76.45;
dicloro6=allyl6/((r+1)*s);
doscloro6=dicloro6*r;
hcl6=allyl6+doscloro6+2*dicloro6;
cloro5=hcl6/c;
propileno5=rr*cloro5;
cloro6=cloro5-hcl6;
propileno6=propileno5-allyl6-doscloro6-dicloro6;
propileno8=propileno6;
propileno4=propileno5;
propileno3=propileno4;
propileno1=propileno3-propileno8;
cloro4=cloro5;
cloro2=cloro4;
hcl9=hcl6+(2*cloro6);
Hcl9=36.45*hcl9;
corriente9=Hcl9*(100/32);
h2o=(0.68*corriente9)+(18*cloro6);
Ganancia=(1.72*p)+(0.075*corriente9)-(0.20*(dicloro6*110.9))-(0.2*((76.45*doscloro6)-(0.005*p)))-(0.0003*h2o)-(0.47*(propileno1*42))-(0.18*(cloro2*70.9));
plot(T,Ganancia,'r'), hold on
title('Rendimiento del proceso de producción de cloruro de alilo', 'FontSize', 12);
xlabel('Temperatura', 'FontSize',12);
ylabel('Ganancia', 'FontSize', 12)
end
else
disp('el reactor no puede trabajar en estas condiciones')
end