Ver vectores

Añadir carga

Elipse

Catenaria

t0: {0, 0, 0}
t1: {0, 0, 0}
g: (0, -9.81)
A = (3.9, 23.1)
B = (107.9, 16.2)
c1 = Ellipse(A, B, B + (2 * (B - A) / abs(B - A)))
M1 = Point(c1)
tt = 352
dt = 0
anima = 0 (slider from 0 to 1, step 0.001)
i1 = Tangent(M1, c1)
n1 = UnitVector(i1)
gt1: (((g * n1)) * n1)
j1 = Segment(M1 + g, M1 + gt1)
v1 = (0, 0)
vt1: ((sgn((v1 * n1)) * abs(v1)) * n1)
C1 = Point(c1)
tirolina = PolyLine(A, M1, B)
reg1: {(5.26167601809081, 17.68919251188505)}
pol1 = PolyLine(reg1)
c2 = Circle(M1, 4)
M2 = Point(c2)
f = Segment(M1, M2)
C2 = (3.110692653106716, 14.31676451953863)
i2 = Tangent(M2, c2)
n2 = UnitVector(i2)
v2 = (0, 0)
vt2: ((sgn((v2 * n2)) * abs(v2)) * n2)
reg2: {(3.110692653106716, 14.31676451953863)}
pol2 = PolyLine(reg2)
gt2: (((g * n2)) * n2)
j2 = Segment(M2 + g, M2 + gt2)
verVectores = true
verCarga = true
vMax1 = 0
vMax2 = 0
M1': M1
texto1: "Velocidad máxima: " + vMax1 + " m/s "
texto2: "Velocidad máxima
de la carga: " + vMax2 + " m/s"
verElipse = true
verCatenaria = true
L: abs(B - A) + 4
a: abs(x(B - A)) / ((2 * x(Root[sinh(x) - (sqrt(L^(2) - (y(B - A))^(2)) / abs(x(B - A)) * x), 20])))
x0: Min[x(A), x(B)] + (a * ln(((L - (((2 * (x(B - A) > 0)) - 1) * y(B - A))) / a * ℯ^(abs(x(B - A)) / a) / (ℯ^(abs(x(B - A)) / a) - 1))))
y0: y(A) - (a * cosh((x0 - x(A)) / a))
a': x(Root[(((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) + (abs(x(B - A)) * x)) * sqrt(1 + ((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) + (abs(x(B - A)) * x))^(2))) - (((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) - (abs(x(B - A)) * x)) * sqrt(1 + ((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) - (abs(x(B - A)) * x))^(2))) + asinh((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) + (abs(x(B - A)) * x)) - asinh((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) - (abs(x(B - A)) * x)) - ((4 * L) * x), 2000])
b: (((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) - ((2 * Min[x(A), x(B)]) * a') - (a' * abs(x(B - A)))
c: (a' * Min[x(A), x(B)]^(2)) - ((((2 * (x(A - B) < 0)) - 1) * y(B - A) / abs(x(B - A))) * Min[x(A), x(B)]) + ((a' * abs(x(B - A))) * Min[x(A), x(B)]) + If[x(A) < x(B), y(A), y(B)]
cat(x) = If[Min[x(A), x(B)] ≤ x ≤ Max[x(A), x(B)], y0 + (a * cosh((x - x0) / a))]
tPausa = 0
tPausaAux = 0
t: ((t1(3) < t0(3)) * 3600) + ((t1(3) - t0(3)) * 60) + t1(2) - t0(2) + (t1(1) - t0(1)) / 1000 + tPausa
animacion = false