pro5

Código de programación Matlab

% PROBLEMA 5
clc
clear
format compact
% Para la lectura manométrica en el manómetro de Bourdon PA=14.65 Kpa,
% determinar:
% a)La elevación [m] en la ramas abiertas de los piezómetros E,F,G.
% b) la lectura H [m] del manómetro de U de mercurio (S=13.6)
disp(‘Problema 6 para calcular la elevación y la lectura de h en un manómetro Bourdon y de U de mercurio’)
p1=input(‘¿Cuál es la presión en KPa? ‘);
%p1=17.65
E1=input(‘¿Cuál es la elevación 1 [m]? ‘);
%E1=15
E2=input(‘¿Cuál es la elevación 2 [m]? ‘);
%E2=12
E3=input(‘¿Cuál es la elevación 3 [m]? ‘);
%E3=8
E4=input(‘¿Cuál es la elevación 4 [m]? ‘);
%E4=4
dr1=input(‘Dame el valor de la densidad relativa de la sustancia 1 = ‘);
%dr1=.7
dr2=input(‘Dame el valor de la densidad relativa de la sustancia 2 = ‘);
%dr2=1
dr3=input(‘Dame el valor de la densidad relativa de la sustancia 3 = ‘);
%dr3=1.6
dr4=input(‘Dame el valor de la densidad relativa de la sustancia 4 = ‘);
%dr4=13.6
pea=input(‘Dame el valor del peso específico del agua en [kg/m^3]= ‘);
%pea=999.64
a=101.9716;
h1=(p1*a)/(dr1*pea)
p2=0;
p3=0;
disp(‘Altura en la columna E [m]’)
He=E1-h1
p4=-(p1*a)+(E1-E2)*(dr1*pea)
p5=p4;
h2=p5/pea
disp(‘Altura en la columna F [m]’)
Hf=E2+h2
p6=p4+(E2-E3)*pea
h3=p6/(dr3*pea)
disp(‘Altura en la columna G [m]’)
Hg=h3+E3
p8=p6;
p9=p6+(E3-E4)*pea
h4=p9/(dr4*pea)
disp(‘Altura en la columna de mercurio en [m]’)
Hhg=E4+h4