MateWiki - Contribuciones del usuario [es]

La catenaria (grupo 24)

2024-12-10T20:00:19Z

Raquel.rserramo:

{{ TrabajoED | La catenaria. Grupo 24| [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/24|2024-25]] |
*David Santafé Palacios
*Pedro Suñé Pérez
*Beatriz Bernal Castañeda
*Raquel Roque Serrano }}
[[Categoría:Teoría de Campos]]
[[Categoría:TC24/25]]
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrización que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:Olympiastadion_Montreal.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Parábola y = A + x^2/A y catenaria. Semejanzas=
En la siguiente imagen podemos ver que la catenaria y la curva que se representa en el dibujo son dos cruvas del mismo tipo, parábolas. El único punto en el que coinciden es en el vértice. Las dos curvas tienen su vértice en el punto (0,2).
[[Archivo:9G24I.jpg|miniaturadeimagen|right|500 px|'''GRÁFICADE LA Y LA CATENARIA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
%Dibujar la parábola
z=linspace(-1,1,1000);
F=a+((z.^2)./a);

% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
hold on
plot(z,F)
hold off
grid on;
}}
 
 

=Catenoide=
===Como se dibuja la superficie===
Para dibujar la catenaria en <math> \mathbb{R}^3 </math> se tiene que parametrizar en coordenadas cartesianas de la siguiente manera:
 

<math>γ(t) = (x1(t), x2(t), x3(t)) = (0, cosh(t), t), t ∈ (−1, 1)</math>. 

[[Archivo:10G24I.jpg|miniaturadeimagen|right|500 px|'''LA CATENARIA EN ℝ3''' ]]

{{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
X1 = cosh(U) .* cos(V);
X2 = cosh(U) .* sin(V);
X3 = U;
% Dibujar la superficie de revolución
figure;
surf(X1, X2, X3, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje x3');
xlabel('x_1');
ylabel('x_2');
zlabel('x_3');
axis equal;
grid on;
}}
 
 

===Información de la superficie y Un ejemplo.===
 Dichas superficie recibe el nombre de Catenoide. Una catenoide es un tipo de superficie que se obtiene por rotación de una catenaria alrededor de un eje coplanario, perpendicular al eje de simetría y que no la corte. Se trata de una superficie mínima o minimal, lo que significa que ocupa el área mínima cuando está delimitada por un espacio cerrado. Fue descrita formalmente en 1744 por el matemático Leonhard Euler.
 Esta superficie se utiliza en estructuras de Ingeniería Civil como:
 
[[Archivo: Mcdonnell planetarium slsc.jpg|miniaturadeimagen|center|600 px|'''El planetario de El St. Louis Science Center en Misouri''' ]]

=Masa de la Catenoide =
La densidad de la superficie proporcionada es <math>f(x1,x2,x3)=x3^2(x1^2+x2^2)</math>. Para describir cómo se distribuye la densidad a lo largo de la superficie, podemos considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización es <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Entonces, la densidad <math>f</math> a lo largo de la superficie está dada por:
<math>f(u,v)=U^2((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math>
Esto significa que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>x3</math>. En otras palabras, a medida que nos alejamos del eje <math>x3</math>, la densidad aumenta.

===Cálculo de la Masa de la Superficie===
La masa de la superficie con la densidad <math>f</math> dada se puede calcular utilizando la integral de la densidad sobre la superficie. La masa M se calcula de la siguiente manera:
<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, puedes 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((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math> en la integral.

{{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 = (X3.^2).*(X1.^2+X2.^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)]);
}}
 
 
===Resultado de la Masa===
El programa de Matlab nos devuelve lo siguiente:
 La masa de la superficie es aproximadamente: 9.7204.

=Referencias=
 '''Bibliografía'''

* Apuntes de clase y de moodle de Teoría de Campos, profesor Carlos Castro.
* Apuntes de informática del año pasado, del profesor Juan Viu Sos.
 '''Enlaces externos'''
* [https://es.wikipedia.org/wiki/Catenaria La catenaria]
* https://www.uv.es/~ivorra/Libros/Catenaria.pdf

La catenaria (grupo 24)

2024-12-09T09:46:42Z

Raquel.rserramo:

{{ TrabajoED | La catenaria. Grupo 24| [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/24|2024-25]] |
*David Santafé Palacios
*Pedro Suñé Pérez
*Beatriz Bernal Castañeda
*Raquel Roque Serrano }}
[[Categoría:Teoría de Campos]]
[[Categoría:TC24/25]]
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:Olympiastadion_Montreal.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Parábola y = A + x^2/A y catenaria. Semejanzas=
En la siguiente imagen podemos ver que la catenaria y la curva que se representa en el dibujo son dos cruvas del mismo tipo, parábolas. El único punto en el que coinciden es en el vértice. Las dos curvas tienen su vértice en el punto (0,2).
[[Archivo:9G24I.jpg|miniaturadeimagen|right|500 px|'''GRÁFICADE LA Y LA CATENARIA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
%Dibujar la parábola
z=linspace(-1,1,1000);
F=a+((z.^2)./a);

% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
hold on
plot(z,F)
hold off
grid on;
}}
 
 

=Catenoide=
===Como se dibuja la superficie===
Para dibujar la catenaria en <math> \mathbb{R}^3 </math> se tiene que parametrizar en coordenadas cartesianas de la siguiente manera:
 

<math>γ(t) = (x1(t), x2(t), x3(t)) = (0, cosh(t), t), t ∈ (−1, 1)</math>. 

[[Archivo:10G24I.jpg|miniaturadeimagen|right|500 px|'''LA CATENARIA EN ℝ3''' ]]

{{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
X1 = cosh(U) .* cos(V);
X2 = cosh(U) .* sin(V);
X3 = U;
% Dibujar la superficie de revolución
figure;
surf(X1, X2, X3, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje x3');
xlabel('x_1');
ylabel('x_2');
zlabel('x_3');
axis equal;
grid on;
}}
 
 

===Información de la superficie y Un ejemplo.===
 Dichas superficie recibe el nombre de Catenoide. Una catenoide es un tipo de superficie que se obtiene por rotación de una catenaria alrededor de un eje coplanario, perpendicular al eje de simetría y que no la corte. Se trata de una superficie mínima o minimal, lo que significa que ocupa el área mínima cuando está delimitada por un espacio cerrado. Fue descrita formalmente en 1744 por el matemático Leonhard Euler.
 Esta superficie se utiliza en estructuras de Ingeniería Civil como:
 
[[Archivo: Mcdonnell planetarium slsc.jpg|miniaturadeimagen|center|600 px|'''El planetario de El St. Louis Science Center en Misouri''' ]]

=Masa de la Catenoide =
La densidad de la superficie proporcionada es <math>f(x1,x2,x3)=x3^2(x1^2+x2^2)</math>. Para describir cómo se distribuye la densidad a lo largo de la superficie, podemos considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización es <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Entonces, la densidad <math>f</math> a lo largo de la superficie está dada por:
<math>f(u,v)=U^2((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math>
Esto significa que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>x3</math>. En otras palabras, a medida que nos alejamos del eje <math>x3</math>, la densidad aumenta.

===Cálculo de la Masa de la Superficie===
La masa de la superficie con la densidad <math>f</math> dada se puede calcular utilizando la integral de la densidad sobre la superficie. La masa M se calcula de la siguiente manera:
<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, puedes 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((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math> en la integral.

{{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 = (X3.^2).*(X1.^2+X2.^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)]);
}}
 
 
===Resultado de la Masa===
El programa de Matlab nos devuelve lo siguiente:
 La masa de la superficie es aproximadamente: 9.7204.

=Referencias=
 '''Bibliografía'''

* Apuntes de clase y de moodle de Teoría de Campos, profesor Carlos Castro.
* Apuntes de informática del año pasado, del profesor Juan Viu Sos.
 '''Enlaces externos'''
* [https://es.wikipedia.org/wiki/Catenaria La catenaria]
* https://www.uv.es/~ivorra/Libros/Catenaria.pdf

La catenaria (grupo 24)

2024-12-09T09:44:11Z

Raquel.rserramo:

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:Olympiastadion_Montreal.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Parábola y = A + x^2/A y catenaria. Semejanzas=
En la siguiente imagen podemos ver que la catenaria y la curva que se representa en el dibujo son dos cruvas del mismo tipo, parábolas. El único punto en el que coinciden es en el vértice. Las dos curvas tienen su vértice en el punto (0,2).
[[Archivo:9G24I.jpg|miniaturadeimagen|right|500 px|'''GRÁFICADE LA Y LA CATENARIA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
%Dibujar la parábola
z=linspace(-1,1,1000);
F=a+((z.^2)./a);

% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
hold on
plot(z,F)
hold off
grid on;
}}
 
 

=Catenoide=
===Como se dibuja la superficie===
Para dibujar la catenaria en <math> \mathbb{R}^3 </math> se tiene que parametrizar en coordenadas cartesianas de la siguiente manera:
 

<math>γ(t) = (x1(t), x2(t), x3(t)) = (0, cosh(t), t), t ∈ (−1, 1)</math>. 

[[Archivo:10G24I.jpg|miniaturadeimagen|right|500 px|'''LA CATENARIA EN ℝ3''' ]]

{{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
X1 = cosh(U) .* cos(V);
X2 = cosh(U) .* sin(V);
X3 = U;
% Dibujar la superficie de revolución
figure;
surf(X1, X2, X3, 'EdgeColor', 'none');
title('Superficie de Revolución de la Catenaria alrededor del eje x3');
xlabel('x_1');
ylabel('x_2');
zlabel('x_3');
axis equal;
grid on;
}}
 
 

===Información de la superficie y Un ejemplo.===
 Dichas superficie recibe el nombre de Catenoide. Una catenoide es un tipo de superficie que se obtiene por rotación de una catenaria alrededor de un eje coplanario, perpendicular al eje de simetría y que no la corte. Se trata de una superficie mínima o minimal, lo que significa que ocupa el área mínima cuando está delimitada por un espacio cerrado. Fue descrita formalmente en 1744 por el matemático Leonhard Euler.
 Esta superficie se utiliza en estructuras de Ingeniería Civil como:
 
[[Archivo: Mcdonnell planetarium slsc.jpg|miniaturadeimagen|center|600 px|'''El planetario de El St. Louis Science Center en Misouri''' ]]

=Masa de la Catenoide =
La densidad de la superficie proporcionada es <math>f(x1,x2,x3)=x3^2(x1^2+x2^2)</math>. Para describir cómo se distribuye la densidad a lo largo de la superficie, podemos considerar la variación de <math>f</math> en la parametrización de la superficie. En este caso, la parametrización es <math>x(u,v)=(cosh(u)cos(v),cosh(u)sin(v),u)</math>.
Entonces, la densidad <math>f</math> a lo largo de la superficie está dada por:
<math>f(u,v)=U^2((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math>
Esto significa que la densidad aumenta cuadráticamente con respecto a u, la coordenada a lo largo del eje <math>x3</math>. En otras palabras, a medida que nos alejamos del eje <math>x3</math>, la densidad aumenta.

===Cálculo de la Masa de la Superficie===
La masa de la superficie con la densidad <math>f</math> dada se puede calcular utilizando la integral de la densidad sobre la superficie. La masa M se calcula de la siguiente manera:
<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, puedes 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((cosh(u)cos(v))^2+(cosh(u)sin(v))^2)</math> en la integral.

{{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 = (X3.^2).*(X1.^2+X2.^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)]);
}}
 
 
===Resultado de la Masa===
El programa de Matlab nos devuelve lo siguiente:
 La masa de la superficie es aproximadamente: 9.7204.

=Referencias=
 '''Bibliografía'''

* Apuntes de clase y de moodle de Teoría de Campos, profesor Carlos Castro.
* Apuntes de informática del año pasado, del profesor Juan Viu Sos.
 '''Enlaces externos'''
* [https://es.wikipedia.org/wiki/Catenaria La catenaria]
* https://www.uv.es/~ivorra/Libros/Catenaria.pdf

2024-12-03T12:24:44Z

Raquel.rserramo: /* Ejemplos de la Catenaria */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:Olympiastadion_Montreal.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=

Archivo:Olympiastadion Montreal.jpg

2024-12-03T12:23:52Z

Raquel.rserramo:

Archivo:Le-stade-olympique.jpeg

2024-12-03T12:20:03Z

Raquel.rserramo:

La catenaria (grupo 24)

2024-12-03T12:17:09Z

Raquel.rserramo: /* Ejemplos de la Catenaria */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:EstadioOlimpico1.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=

La catenaria (grupo 24)

2024-12-03T12:15:43Z

Raquel.rserramo: /* Ejemplos de la Catenaria */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Puente de Vedra, Galicia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Puente de Clifton''' ]]
 
[[Archivo:EstadioOlimpico.jpg|miniaturadeimagen|center|700 px|'''Estadio Olímpico de Canadá''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Cupula de Reichstag''' ]]
 

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=

La catenaria (grupo 24)

2024-12-03T12:13:53Z

Raquel.rserramo: /* Ejemplos de la Catenaria */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:PuentedeClifton.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:EstadioOlimpico.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 
[[Archivo:CupuladeReichstag.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
 

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=

La catenaria (grupo 24)

2024-12-03T12:11:47Z

Raquel.rserramo: /* Ejemplos de la Catenaria */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t/2))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente <math>\vec n(t) \neq \vec 0 </math> , por lo que se tiene que cumplir <math> γ′′(t) = κ(t)\vec n(t) </math>
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''GRÁFICA DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{(cosh(t/2)/2))}{(1+senh(t/2)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}

= Circunferencia osculatriz =
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5)=0 \vec i + 4 \vec j =4 \vec j </math></center>
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:6.1G24I.jpg|miniaturadeimagen|right|700 px|'''GRÁFICADE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=La catenaria =
La catenaria es la curva cuya forma es la que adopta una cuerda de densidad uniforme sujeta por sus dos extremos y sometida únicamente a la fuerza de la gravedad.
 La catenaria tiene como ecuación:
<center><math>y=acosh(\frac{x}{a})</math></center>
 En lo general la ecuación de la catenaria se refiere a cadenas o cuerdas infinitamente flexibles e inextensibles. El requisito de flexibilidad infinita se refiere a que la rigidez flexional sea nula y el requisito de inextensibilidad se refiere a que la longitud de cada tramo de la misma no varíe a pesar de estar sometido a fuerzas. Obviamente en las cuerdas reales estos requisitos se cumplen sólo de forma aproximada. Para cuerdas de gran longitud, la elasticidad de la cuerda las aleja del comportamiento perfectamente inextensible. Si bien la catenaria de una cuerda inextensible es siempre una curva plana, para cables gruesos de pequeña longitud la rigidez flexional finita hace que su deformada no necesariamente esté contenida en un plano.
 Una aplicación de la catenaria se ve por ejemplo en la arquitectura. Dado un elemento lineal sometido solo a cargas verticales, la forma catenaria es precisamente la forma del eje baricéntrico que minimiza las tensiones. Esa propiedad puede aprovecharse para el diseño de arcos. De este modo un arco en forma de catenaria invertida es precisamente la forma en la que se evita la aparición de esfuerzos distintos de los de compresión, como son los esfuerzos cortantes o los flectores.
 Por esa razón, una curva catenaria invertida es un trazado útil para un arco en la arquitectura, forma que fue aplicada, entre otros, y fundamentalmente, por Antonio Gaudí, quien utilizaba maquetas para el estudio de estas estructuras, para luego volcar aquella información en los planos de sus obras.

=Ejemplos de la Catenaria =
[[Archivo:Sagradafamilia-passion.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]
[[Archivo:Ponteulla_Vedra_Galicia_03.jpg|miniaturadeimagen|center|700 px|'''Arcos de la Sagrada Familia''' ]]

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=

Archivo:PuentedeClifton.jpg

2024-12-03T12:10:49Z

Raquel.rserramo:

La catenaria (grupo 24)

2024-12-03T11:29:20Z

Raquel.rserramo: /* Sea P = γ(0,5), es decir t = 0,5. Hallar el centro y el radio de la circunferencia osculatriz en P (se pueden aproximar), y dibujar la circunferencia osculatriz junto a la curva. */

{{ TrabajoED | La Catenaria. Grupo 24 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] | David Santafé Palacios Pedro Suñé Pérez Beatriz Bernal Castañeda Raquel Roque Serrano}}
La catenaria es una curva fundamental en ingeniería civil, ya que describe la forma que adoptan los cables o cadenas flexibles suspendidos bajo su propio peso, sin otras fuerzas externas. Su aplicación es clave en el diseño de estructuras como puentes colgantes, líneas de transmisión eléctrica y cubiertas tensadas, donde la eficiencia estructural y la distribución uniforme de tensiones son esenciales. Este trabajo explora las propiedades matemáticas de la catenaria, su relación con otras curvas como la parábola y sus principales aplicaciones prácticas, destacando su importancia en la optimización y seguridad de proyectos ingenieriles. 

[[Archivo:CatenariaGI.png|miniatura|400 px|centro|'''FORMA DE LA CATENARIA''']]

Siendo la curva representada por:
 '''<math> γ(t) = (x(t),y(t)) = (t,Acosh(t/A)), t∈(-1,1)</math>''' 

Para la representación y cálculos de a continuación usaremos el programa Matlab 

=Dibujar la curva=
[[Archivo:1CatenariaG24I.jpg|500 px||miniaturadeimagen|right|'''REPRESENTACIÓN CATENARIA''' ]]
===Código===
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 1000);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;}}
 
 
 
 
 
 
 

=Calcular los vectores velocidad γ'(t) y aceleración γ"(t), y dibujarlos junto a la curva=

===Definición vector posición, velocidad y aceleración===

El vector posición es el que va desde el origen del sistema de referencia hasta la ubicación de la partícula. El vector velocidad se obtiene al derivar el vector posición con respecto al tiempo. Por su parte, el vector aceleración es la derivada del vector velocidad. En general, el vector aceleración γ′′(t) no tiene por qué ser ortogonal al vector velocidad γ′(t). Sin embargo, sí lo será si la curva γ(t) está parametrizada por la longitud de arco, es decir, cuando la magnitud de la velocidad ∣γ′(t)∣=1|. 

 <center><math>γ(t) = (x(t),y(t)) = (t,2cosh(t/2)) </math></center>

 <center><math>γ´(t) = x′(t)\vec i + y′(t)\vec j = 1\vec i +(senh(t/2))\vec j </math></center>

 <center><math>γ´´(t) = x′′(t)\vec i + y′′(t)\vec j = 0\vec i + \frac{cosh(t/2)}{2}\vec j </math></center>

 

=== Representación de los vectores ===

[[Archivo:2.2G24I.jpg|600px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
a=2;
t = linspace(-1, 1, 20);
x = t;
y = a*cosh(t/a);
% Dibujar la curva
figure;
plot(x, y, 'LineWidth', 2);
title('Curva parametrizada: \gamma(t) = (t, a*cosh(t/a))');
xlabel('x');
ylabel('y');
grid on;

% Velocidad y aceleración
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);

% Gráfica
figure
hold on
plot(x, y, 'r');
quiver(x, y, V1, V2, 1, "Color", "c");
quiver(x, y, A1, A2, 1, "color", "m");
axis equal
hold off;

% Centrado de la grafica en el origen de coordenadas
ax = gca;
ax.XAxisLocation = 'origin';
ax.YAxisLocation = 'origin';

% Etiquetas
xlabel("x", "FontSize", 10);
ylabel("y", "FontSize", 10);}}

=Longitud de curva=

Hemos realizado el calculo de la longitud total de la curva mediante código de MATLAB. Se sabe que la longitud de una curva no depende de la parametrización elegida. La longitud se trata por tanto de un concepto asociado a la propia curva y no a la parametrizacioón que la describe.

=== Código MATLAB longitud de la curva ===
[[Archivo:3.1G24I.jpg|600 px||miniaturadeimagen|right|'''GRÁFICA LONGITUD DE CURVA''' ]]

{{matlab|codigo=
% Calcular longitud de la curva
a=2;
t=linspace(-1,1,50); %valores a t como vector en el intervalo dado
x = t;
y = a*cosh(t/a);
longitud = sum(sqrt(diff(x).^2 + diff(y).^2)); % ponemos la fórmula de la longitud de curva como la integral del módulo del vector velocidad
disp(['Longitud de la curva: ' num2str(longitud)]); % imprime el resultado obtenido como longitud de curva

% Calcular longitud de la curva y longitud acumulativa
longitud_seg = sqrt(diff(x).^2 + diff(y).^2);
longitud_acum = [0, cumsum(longitud_seg)]; % Agregar un cero al inicio para que tengan la misma longitud

% Ajustar t para que tenga la misma longitud que longitud_acum
t_acum = linspace(-1, 1, length(longitud_acum));

% Mostrar la longitud de la curva numéricamente
longitud_total = longitud_acum(end);
disp(['Longitud total de la curva: ' num2str(longitud_total)]);

% Graficar la longitud
figure;
plot(t_acum, longitud_acum, 'LineWidth', 2);
title('Longitud de la curva');
xlabel('t');
ylabel('Longitud ');
grid on;
}}

Se obtiene una longitud total de la curva de 2.0843

=Vectores tangente '''t'''(t) y normal '''n'''(t)=
 <center><math> ℓ(γ) = \int_{a}^{b}|γ′(t)|=\int_{a}^{b}\sqrt {x´(t)^2 +y´(t)^2}dt= \int_{0}^{2π}\sqrt{(1)^2 +(sinh(t))^2}</math></center>

Un vector unitario tangente tiene magnitud uno y es tangente a una curva en un punto, indicando el cambio en dirección sin alterar la rapidez. En contraste, el vector unitario normal, también de magnitud uno, es perpendicular a la curva en ese mismo punto. Este último es clave para calcular la fuerza normal en objetos en movimiento, siendo ambos vectores fundamentales para analizar trayectorias y dinámicas. 
Claramente n (t)≠0 n→(t)≠0, por lo que se tiene que cumplir γ''(t)=κ(t)n(t) γ′′(t)=κ(t) n→(t) 
Se pueden obtener a partir de: 
 <center><math> \vec t(t)=\frac{γ′(t)}{|γ′(t)|}=\frac{(x´(t)\vec i +y´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (1-cost)\vec i +(sent)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}} </math></center>
 <center><math> \vec n(t)=\frac{(-y´(t)\vec i +x´(t)\vec j)}{\sqrt{x´(t)^2 + y´(t)^2}}=\frac{ (-sent)\vec i +(1-cost)\vec j}{\sqrt{(1-cost)^2 +(sent)^2}}</math></center>

===REPRESENTACIÓN===
 
 
[[Archivo:4.1G24I.jpg|550px|miniaturadeimagen|right|'''REPRESENTACIÓN DEL VECTOR TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definición de los vectores normales y tangentes
t= linspace(-1,1,20); %Parámetro
a=2;
x= t;
y= a*cosh(t/a);
% Velocidades/tangentes/normales
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
mod= sqrt(V1.^2+V2.^2);
t1= V1./mod;
t2= V2./mod;
n1= -t2;
n2= t1;
%Representación
figure
axis equal
hold on
plot (x ,y ,'b') ;
quiver(x,y,t1,t2,1,"Color","r") ;
quiver(x,y,n1,n2,1,"Color","g") ;
grid on
hold off;
title('Curva, tangente y normal' )
}}
 
 

=Cálculo de curvatura k(t)=
En este análisis de la curva parametrizada γ(t)=(t,2cosh(t/2)), nos enfocamos en estudiar su curvatura κ, un parámetro que indica cuán alejada está la curva de ser una línea recta en cada punto. Graficando κ(t) , podemos observar cómo varía la "tensión" o el "giro" de la curva a lo largo de su recorrido, lo que nos da una comprensión más profunda de su comportamiento geométrico local.
 Para calcularla se usará la siguiente expresión:
 <center><math>Curvatura: κ(t)=\frac{x´(t)y´´(t)-x´´(t)y´(t)}{(x´(t)^2+y´(t)^2)^(3/2}=\frac{cosh(t)-0senh(t)}{(1^2+senh(t)^2)^(3/2}=\frac{cosh(t)}{(1+senh(t)^2)^(3/2)} </math></center>

=== Código MATLAB de la curvatura ===

[[Archivo:5.1G24I.jpg|500px|miniaturadeimagen|right|'''GRÁFICA DE LA CURVATURA''' ]]
{{matlab|codigo=
% Definir la parametrización
t = linspace(-1, 1, 50);
x = t;
y = 2*cosh(t/2);
% Velocidad y aceleración
V1 = ones(size(t));
V2 = sinh(t/2);
A1 = zeros(size(t));
A2 = cosh(t/2)/2;
% Calcular la curvatura
numerador = V1 .* A2 - V2 .* A1;
denominador = (V1.^2 + A1.^2).^(3/2);
curvatura = abs(numerador ./ denominador);
% Graficar la curvatura en función de t
figure;
plot(t, curvatura, 'LineWidth', 2);
title('Gráfica de la Curvatura \kappa(t)');
xlabel('t');
ylabel('\kappa(t)');
grid on;
}}
= Sea P = γ(0,5), es decir t = 0,5. Hallar el centro y el radio de la circunferencia osculatriz en P (se pueden aproximar), y dibujar la circunferencia osculatriz junto a la curva. =
=Circunferencia osculatriz=
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5) </math></center>Q(0.5)=(0\veci+4\vecj)=(4\vecj)
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:Circunferenci osculatriz.jpg|miniaturadeimagen|left|700 px|'''REPRESENTACIÓN DE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=Circunferencia osculatriz=
La circunferencia osculatriz de una curva en P es la circunferencia contenida en el plano osculador Ω, de centro el centro y radio de la curvatura.
 La circunferencia osculatriz proporciona una buena aproximación a la curva en las proximidades del punto \(P\).
 
En nuestro caso, <math>P=(γ(0.5),0.5).</math>, por lo que el radio y el centro se calculan de la siguiente manera:
 
 <center><math>Q(t) = γ(t) + \frac{1}{κ(t)} \vec{n}(t)= (t,Acosh(t/A))+ \frac{1}{\frac{(cosh(t/2))/2)}{(1+sinh^2(t/2))^\frac{3}{2}}} \frac{(-sinh(t/2)\vec i + 1 \vec j )}{\sqrt{1+sinh(t/2)^2}}</math></center>
 
 <center><math>Q(0.5) = γ(0.5) + \frac{1}{κ(0.5)} \vec{n}(0.5) </math></center>Q(0.5)=(0\veci+4\vecj)=(4\vecj)
 
 <center><math>R(t)=\frac{1}{|κ(t)|}= \frac{1}{|\frac{(cosh(t/2))/2}{(1+sinh^2(t/2))^\frac{3}{2}}|} </math></center>
 
 <center><math>R(0.5)=\frac{1}{|κ(0.5)|} = 2m </math></center>
 
 

=== Código para representar la curva con su circunferencia osulatriz ===
[[Archivo:Circunferenci osculatriz.jpg|miniaturadeimagen|left|700 px|'''REPRESENTACIÓN DE LA CIRCUNFERENCIA Y LA CATENARIA''' ]]
{{matlab|codigo=
%Definir centro y radio
%radio
a=2;
t=0;
V1 = ones(size(t));
V2 = (a/a)*sinh(t/a);
A1 = zeros(size(t));
A2 = (a/a^2)*cosh(t/a);
K=((V1.*A2)-(V2.*A1))./(((V1.^2)+(A1.^2)).^(3/2));
radio=(1/K);

%centro
normal1=(-A1)./(sqrt(V1.^2+A1.^2));
normal2=(V1)./(sqrt(V1.^2+A1.^2));
normal=[normal1,normal2];
centro=[t,a*cosh(t/a)]+(1/K).*normal;

% 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 = centro(1) + radio * cos(theta);
y_circunferencia = centro(2) + radio * sin(theta);

% Definir la parametrización de la catenaria
t = linspace(-1.5, 1.5, 100);
a=2;
x_catenaria = t;
y_catenaria = a*cosh(t/a);

% Dibujar la circunferencia y la catenaria en la misma gráfica
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);

% Configurar el aspecto del gráfico
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
}}
 
 

=Buscar información de la curva dada, y explicar brevemente qué fenómeno describe. Añadir cualquier información que se considere relevante o interesante, especialmente si tiene que ver con el campo de la ingeniería. =

=Mostrar una foto de alguna estructura civil en la que se haya usado la curva dada. =

=Dibujar en la misma figura la catenaria y la parábola de ecuación y = A + x^2/A, y explicar cómo se asemejan ambas. =

=Podemos ver la catenaria en ℝ3 mediante la parametrización en cartesianas=

=Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: =

=Referencias=