model circuitoHidraulico_3 constant Real PI = 2*Modelica.Math.asin(1.0); constant Real g = 10.0; parameter Real rho = 1000.0; parameter Real nu = 1e-6; 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 L = 24.0; parameter Real D = 130e-3; parameter Real A = PI*D^2/4; parameter Real e = 46e-6; parameter Real epsilon = e/D; parameter Real Kentrada = 0.5; parameter Real Ksaida = 1.0; parameter Real Kcurva = 1.3; Real h1, h2; Real q; Real p1, p2; Real v(start=1.0); Real Re; Real f; Real K; Real R; 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; f = fatorDeAtrito(Re,epsilon); K = f*L/D + Kentrada + 2*Kcurva + Ksaida; R = rho*K/(2*A^2); Re = v*D/nu; A * v = q; annotation( Documentation(info = "Sistema com dois reservatórios hidráulicos e fator de carga calculado pela fórmula de Colebrrok.")); end circuitoHidraulico_3;