MateWiki - Contribuciones del usuario [es]

La Catenaria

2023-12-11T16:55:30Z

Alejandrojg: /* Introducción */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cate6.jpg|900px|miniaturadeimagen|centro|Cálculo centro y radio]]
{{matlab|codigo=
% Definimos el radio y el centro de la circunferencia
t=0
K=cosh(t)/(1+sinh(t)^2)^(3/2)
radio=1/K
n=1/sqrt(1+sinh(t)^2)*[-sinh(t),1];
Q=[t,cosh(t)]+(1/K)*n

% Crear un conjunto de ángulos para generar puntos en la circunferencia
theta = linspace(0, 2*pi, 100);

% Calcular las coordenadas de los puntos en la circunferencia
x_circunferencia = Q(1) + radio * cos(theta);
y_circunferencia = Q(2) + radio * sin(theta);

% Parametrización de la catenaria
tt =-1:0.09:1;
x_catenaria = tt;
y_catenaria = cosh(tt);

% Dibujo de la circunferencia y la catenaria
figure;
plot(x_circunferencia, y_circunferencia, 'b-', 'LineWidth', 2);
hold on; % Mantener la gráfica actual y agregar nuevas
plot(x_catenaria, y_catenaria, 'r-', 'LineWidth', 2);

axis equal; % Hace que las unidades en x e y sean iguales
grid on;
title('Circunferencia y Catenaria');
xlabel('Eje X');
ylabel('Eje Y');
legend('Circunferencia', 'Catenaria');
hold off; % Liberar la gráfica actual para futuros trazados
}}

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)
[[Archivo:Caten5.png|400px|miniaturadeimagen|centro|Catenaria Madrid-Aranjuez]]
[[Archivo:Cates5.png|400px|miniaturadeimagen|centro|Puente de Gudián (Galicia)]]

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
[[Archivo:Cat7.png|350px|miniaturadeimagen|derecha|Superficie]]
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie conociendo la densidad en la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:41:12Z

Alejandrojg: /* Cálculo de la masa de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
[[Archivo:Cat7.png|350px|miniaturadeimagen|derecha|Superficie]]
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie conociendo la densidad en la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:40:30Z

Alejandrojg: /* Código */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
[[Archivo:Cat7.png|350px|miniaturadeimagen|derecha|Superficie]]
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:40:11Z

Alejandrojg: /* Código */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
[[Archivo:Cat7.png|400px|miniaturadeimagen|derecha|Superficie]]
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:38:27Z

Alejandrojg: /* Cálculo de la masa de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:37:40Z

Alejandrojg: /* Masa de la superficie con dicha densidad */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X1 = cosh(U) .* cos(V);
X2 = cosh(U) .* sin(V);
X3 = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

La Catenaria

2023-12-11T16:36:57Z

Alejandrojg: /* Cálculo de la masa de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|centro|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==
[[Archivo:Cat6.png|900px|miniaturadeimagen|centro|Cálculo centro y radio]]

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

El código para cálcular la masa de la superficie es el siguiente:

{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 1000);
v = linspace(0, 2*pi, 1000);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X1 = cosh(U) .* cos(V);
X2 = cosh(U) .* sin(V);
X3 = U;
% Densidad
densidad = U.^2;
% Calcular la masa de la superficie
M = trapz(v, trapz(u, densidad .* sqrt((sinh(U).^2 + cosh(U).^2) .* U.^2), 1));
disp(['La masa de la superficie es aproximadamente: ' num2str(M)]);
}}

La masa de la superficie seria 5,12

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:33:01Z

Alejandrojg: /* Cálculo de la masa de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==
[[Archivo:Cat5.png|900px|miniaturadeimagen|derecha|Cálculo]]
===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]

==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>

En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:31:43Z

Alejandrojg: /* Distribución de la densidad a lo largo de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
[[Archivo:Cate1.jpg|900px|miniaturadeimagen|centro|Cálculo]]
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

Una vez conocida la distribución de la densidad en la superficie, se puede calcular la masa de la superficie.

===Cálculo de la masa de la superficie===
La masa de la superficie se obtiene integrando la función densidad dada sobre la superficie. El cálculo es el siguiente:
<math>M=∬SfdS</math>
donde <math>S</math> es la superficie parametrizada.
La integral se puede expresar en coordenadas paramétricas <math>(u,v)</math> como:

<math>M=\int_a^b\int_c^d f(u,v)*\sqrt[]{((∂x/∂u)^2+(∂x∂v)^2+(∂x/∂u)^2)}\; dv \; du</math>
En este caso, se puede utilizar la parametrización <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math> y la densidad <math>f(u,v)=u^2</math> en la integral.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:25:37Z

Alejandrojg: /* Distribución de la densidad a lo largo de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a <math>u</math>, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:24:37Z

Alejandrojg: /* Distribución de la densidad a lo largo de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=z^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:24:11Z

Alejandrojg: /* Distribución de la densidad a lo largo de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x,y,z)=x^2(x,y,z)=z^2</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=z^2(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>z</math>, es decir, a mayor distancia del eje <math>z</math>, mayor es la densidad.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:21:58Z

Alejandrojg: /* Distribución de la densidad a lo largo de la superficie */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

La densidad de la superficie es <math>f(x1,x2,x3)=x^23(x1,x2,x3)=x^23</math>. Para obtener cómo se distribuye la densidad en la superficie, se puede considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización queda <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Siendo la densidad <math>f</math> a lo largo de la superficie:
<math>f(u,v)=x^23(u,v)=u^2</math>
De lo que se deduce que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>x3</math>, es decir, a mayor distancia del eje <math>x3</math>, mayor es la densidad.

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:11:09Z

Alejandrojg: /* Código */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
{{matlab|codigo=
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;
}}

==Distribución de la densidad a lo largo de la superficie==

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:10:19Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

{{ TrabajoED | La Catenaria. Grupo 12 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Alejandro Jiménez García Marta García-Moris Fontcuberta Alejandro Seises López Alberto Nuñez Cobo Álvaro Matías Acedo }}
==Introducción==
La catenaria se define como la curva generada por una cuerda o cable sin rigidez, suspendida por sus dos extremos y sometida al campo gravitatorio. En ingeniería de caminos, esta curva se encuentra en los cables principales de los puentes colgantes o en los cables de las líneas eléctricas ferroviarias, las cuales se denominan catenaria, por la curva que describe.

==Representación de la curva ==
Se tiene la siguiente curva plana parametrizada en coordenadas cartesianas que representa una catenaria:
 '''<math> γ(t) = (x(t),y(t)) = (t,cosh(t)), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^2</math>
 Para dibujarla hacemos uso del software de programación y cálculo numérico Matlab/Octave.

[[Archivo:Catenaria.png|380px|miniaturadeimagen|derecha|Representación gráfica de la catenaria]]
{{matlab|codigo=
% Definición de los parámetros
a=-1; b=1;h=0.01;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
plot(x,y,"Color","b");
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Curva Catenaria')
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Vectores velocidad y aceleración==
Siendo <math> γ(t) </math> el vector posición:
 
El vector velocidad es la derivada con respecto al tiempo del vector posición.
 
 
<math> γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(sen(h)t)\vec j </math>
 
 
El vector aceleración es la derivada segunda de la posición con respecto al tiempo (derivada de la velocidad con respecto al tiempo).
 
 
<math> γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + (cosh(t))\vec j</math>
 
 
Los vectores velocidad y aceleración se representan gráficamente usando el siguiente código:
[[Archivo:Catenaria 2.png|750px|miniaturadeimagen|derecha|Vector velocidad y vector acelración]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.1;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Definición de la velocidad
v1=t./t;
v2=sinh(t);
% Definición de la aceleración
a1=0.*t;
a2=cosh(t);
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,v1,v2);
quiver(x,y,a1,a2, "Color","g");
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation ='origin';
% Etiquetas
title('Vector velocidad y aceleración')
legend("Catenaria","Velocidad","Aceleración")
axis("equal")
xlim([-1 1])
ylim([0.8 2])
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Longitud de la curva==
La longitud de la curva se ha calculado a través del Método del Rectángulo usando el siguiente código:
{{matlab|codigo=
function A = rectangulo(f, a, b, n)
f=@(t) sqrt(1+(sinh(t)).^2)
a=-1
b=1
n=200
A = 0;
h = (b-a)/n;
for k=1:n
xk = a + h*k; % Valor de x en el extremo del intervalo
A = A + f(xk)*h; % Sumamos el área del rectángulo
end

end
}}

[[Archivo:Catenaria3.jpg|680px|miniaturadeimagen|centro|Cálculo de la longitud de la curva]]

Dando como resultado queda de la siguiente forma:

L = 2,35 m

==Vector tangente y vector normal==
===Código vector tangente===
[[Archivo:Catenaria4.png|450px|miniaturadeimagen|derecha|Vector tangente]]
{{matlab|codigo=
x=t;
y=cosh(t);
% Vectores tangentes unitarios interiores
t1i=(t./t)./(sqrt(1+(sinh(t)).^2));
t2i=sinh(t)./(sqrt(1+(sinh(t)).^2));
% Vectores tangentes unitarios
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,t1i,t2i);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vector tangente unitario')
legend("Catenaria","Vector tangente unitario")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}
===Código vector normal===
[[Archivo:Catenaria5.png|550px|miniaturadeimagen|derecha|Vector normal]]
{{matlab|codigo=
% Definición de los parámetros
a=-1;
b=1;
h=0.09;
t=a:h:b;
% Definición de la curva
x=t;
y=cosh(t);
% Vectores normales unitarios orientación interior
n1i=sinh(t)./(sqrt(1+(sinh(t)).^2));
n2i=-(t./t)./(sqrt(1+(sinh(t)).^2));
% Vectores normales unitarios orientación exterior
n1e=-sinh(t)./(sqrt(1+(sinh(t)).^2));
n2e=(t./t)./(sqrt(1+(sinh(t)).^2));
hold on
plot(x,y,'LineWidth',2);
quiver(x,y,n1i,n2i);
quiver(x,y,n1e,n2e);
hold off
% Centrado de la gráfica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
% Etiquetas
title('Vectores normales')
legend("Catenaria","Vector normal interior","Vector normal exterior")
axis("equal")
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';
box on
grid minor
}}

==Curvatura y gráfica==

===Código===
{{matlab|codigo=
a=-1;
b=1;
h=0.09
t=a:h:b %Valores del intervalo
Ka=(sqrt(1+(sinh(t)).^2).^3)./cosh(t);
plot(t,Ka);
grid on
xticks([-1 1])
}}
[[Archivo:Catenaria6.png|400px|miniaturadeimagen|centro|Curva y gráfica]]
==Circunferencia osculatriz en t=0==

==Fenómeno descrito por la curva==

==Empleo de la curva en Ingeniería Civil==
La catenaria en ingeniería civil es una curva muy común, al tratarse de una curva natural sometida a su propio peso. Este fenómeno de una cuerda o cable sometido a gravedad se da en los cables principales de los puentes colgantes, que es el cable de donde cuelgan los tensores que se anclan al tablero. Un ejemplo de un puente de esta tipología puede ser el Golden Gate en San Francisco (Estados Unidos). Otro fenómeno en el que interviene esta curva es en las líneas eléctricas ferroviarias, que son los cables de donde se obtiene la energía eléctrica del tren. Estos cables cuelgan sometidos a su propio peso de unos elementos llamados catenarias dando lugar a la curva en cuestión. También es importante mencionar, que esta curva no solo aparece en cables o cuerdas, también están en los puentes de tipo arco con tablero superior, como puede ser el Puente de Gudián (Galicia)

==Superficie de revolución alrededor del eje vertical==
Podemos ver la catenaria en mediante la siguiente parametrización en coordenadas cartesianas en <math> \mathbb{R}^3 </math>:
 '''<math> γ(t) = (x(t),y(t),z(t)) = (0,cosh(t),t), t∈(-1,1)</math>'''
 Donde:
 <math>γ:t\to\mathbb{R}^3</math>
 Para obtener la superficie de revolución hacemos uso del software de programación y cálculo numérico Matlab/Octave.

===Código===
% Parámetros
u = linspace(-1, 1, 100);
v = linspace(0, 2*pi, 100);
[U, V] = meshgrid(u, v);
% Parametrización de la superficie en coordenadas cilíndricas
X = cosh(U) .* cos(V);
Y = cosh(U) .* sin(V);
Z = U;
% Dibujar la superficie de revolución
figure;
surf(X, Y, Z, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje Z');
xlabel('x');
ylabel('y');
zlabel('z');
axis equal;
grid on;

==Distribución de la densidad a lo largo de la superficie==

==Masa de la superficie con dicha densidad==

[[Categoría:Teoría de Campos]]
[[Categoría:TC23/24]]

La Catenaria

2023-12-11T16:06:26Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

La Catenaria

2023-12-11T16:01:49Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T15:20:16Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T15:15:31Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T15:12:02Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T15:11:15Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T15:10:55Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T15:09:19Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T15:04:10Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T14:58:41Z

Alejandrojg: /* Vectores velocidad y aceleración */

La Catenaria

2023-12-11T14:58:25Z

Alejandrojg: /* Vectores velocidad y aceleración */

La Catenaria

2023-12-11T14:58:12Z

Alejandrojg: /* Vectores velocidad y aceleración */

La Catenaria

2023-12-11T14:57:12Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T14:50:11Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T14:49:38Z

Alejandrojg: /* Longitud de la curva */

Archivo:Velocidad Velocidad Aceleración.png

2023-12-11T14:43:37Z

Alejandrojg:

La Catenaria

2023-12-11T14:42:34Z

Alejandrojg: /* Vectores velocidad y aceleración */

La Catenaria

2023-12-11T14:40:08Z

Alejandrojg: /* Representación de la curva */

La Catenaria

2023-12-11T14:39:56Z

Alejandrojg: /* Representación de la curva */

La Catenaria

2023-12-11T14:36:12Z

Alejandrojg: /* Representación de la curva */

La Catenaria

2023-12-11T14:17:47Z

Alejandrojg: /* Empleo de la curva en Ingeniería Civil */

La Catenaria

2023-12-11T14:16:44Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

La Catenaria

2023-12-11T14:14:29Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

La Catenaria

2023-12-11T14:09:10Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

La Catenaria

2023-12-11T14:08:33Z

Alejandrojg: /* Superficie de revolución alrededor del eje vertical */

La Catenaria

2023-12-11T14:03:24Z

Alejandrojg: /* Curvatura y gráfica */

La Catenaria

2023-12-11T14:01:50Z

Alejandrojg: /* Representación de la curva */

La Catenaria

2023-12-11T14:00:23Z

Alejandrojg: /* Código vector normal */

La Catenaria

2023-12-11T13:59:31Z

Alejandrojg: /* Código vector tangente */

La Catenaria

2023-12-11T13:58:30Z

Alejandrojg: /* Vector tangente y vector normal */

La Catenaria

2023-12-11T13:57:32Z

Alejandrojg: /* Código */

La Catenaria

2023-12-11T13:52:17Z

Alejandrojg: /* Código */

La Catenaria

2023-12-11T13:51:47Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T13:51:20Z

Alejandrojg: /* Longitud de la curva */

La Catenaria

2023-12-11T13:48:39Z

Alejandrojg: /* Vectores velocidad y aceleración */