Cálculo ángulo betha
Publicado por Lelar (4 intervenciones) el 24/12/2021 10:18:32
Hola,
Tengo un código en dónde quiero calcular un díámetro de una sección de canal circular, en donde tengo 2 ecuaciones disponibles, pero cuando calculo en ángulo beta necesario me da un error.
Espero que puedan ayudarme.
Saludos
.
Creo que el problema se encuentra en la función np.arccos.
Tengo un código en dónde quiero calcular un díámetro de una sección de canal circular, en donde tengo 2 ecuaciones disponibles, pero cuando calculo en ángulo beta necesario me da un error.
Espero que puedan ayudarme.
Saludos
.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
import numpy as np
from scipy.optimize import fsolveimport matplotlib.pyplot as plt%matplotlib inlineycf= 2.8220646159801364
Q=108b=4.6
m=1.5
g=9.81
def D1(x):
yc2=x[0]
D =x[1]
A1 = (b+m*ycf)*ycf
Ec =ycf+Q**2/(2*g*A1**2)
aux=yc2/Dbetha=np.arccos(1-2*((aux)))
A2 = (betha-np.sin(np.radians(betha))*np.cos((betha)*(D**2)/4)
B2 = D*np.sin(betha)
F=np.zeros((2))
F[0]=yc2+Q**2/(2*g*A2**2)-Ec #Ec Conservación de Energía
F[1]=(1+0.77*(Q**2/(g*D**5)**-3)**-0.085)-yc2 #Condición de crisis igualada a 0
return F
x_inicial=np.array([1,1])
x=fsolve(D1,x_inicial)
print("La altura critica =",np.round(x[0],2),"m")
print("El diámetro es =",np.round(x[1],2),"m")
Creo que el problema se encuentra en la función np.arccos.
Valora esta pregunta
0
