model circuitoHidraulico_1
constant Real PI = 2*Modelica.Math.asin(1.0);
constant Real g = 10.0;
parameter Real rho = 1000.0;
parameter Real p0 = 0.1e6;
parameter Real A1 = 1.0;
parameter Real A2 = 1.0;
parameter Real H1 = 10.0;
parameter Real H2 = 8.0;
parameter Real R = 25e6;
Real h1, h2;
Real q;
Real p1, p2;
initial equation
h1 = H1;
h2 = H2;
equation
A1 * der(h1) = -q;
A2 * der(h2) = q;
p1 = p0 + rho * g * h1;
p2 = p0 + rho * g * h2;
if (p1-p2)>0 then
q = sqrt(abs(p1-p2)/R);
else
q = -sqrt(abs(p1-p2)/R);
end if;
annotation(
Documentation(info = "Modelo com dois reservatórios hidráulicos e resistência de perda de carga constante.
Solução da incongruência no cálculo do sentido da vazão.
"));
end circuitoHidraulico_1;