Diferencia entre revisiones de «La Cicloide (Grupo 18)»
(→Cálculo de los vectores velocidad 𝛾 ′ (𝑡) y aceleración 𝛾 ″ (𝑡), y dibujo junto a la curva.) |
(→Bibliografía y Referencias) |
||
| (No se muestran 110 ediciones intermedias del mismo usuario) | |||
| Línea 1: | Línea 1: | ||
| + | |||
{{TrabajoED | La Cicloide (Grupo 18) | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC25/26|2025-26]] | Alejandro Porrúa Perea <br/> Salvador Sánchez Burgos <br/> Natalia Andrés Jiménez <br/> Jaime Colomina López <br/> John Cuenca Uyaguari}} | {{TrabajoED | La Cicloide (Grupo 18) | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC25/26|2025-26]] | Alejandro Porrúa Perea <br/> Salvador Sánchez Burgos <br/> Natalia Andrés Jiménez <br/> Jaime Colomina López <br/> John Cuenca Uyaguari}} | ||
==Introducción== | ==Introducción== | ||
| + | [[Archivo:Geometría de la Cicloide.png|right|]] | ||
La '''Cicloide''' se define como la '''curva''' trazada por un punto contenido en la circunferencia de un círculo que rueda '''sin deslizar''' sobre una línea recta. Aunque también se puede entender como el lugar geométrico descrito por las posiciones de un punto P que pertenece a la circunferencia de un círculo que rueda sin deslizamiento sobre una recta. Dicha recta recibe el nombre de '''directriz''', mientras que la circunferencia recibe el nombre de '''generatriz o ruleta'''. | La '''Cicloide''' se define como la '''curva''' trazada por un punto contenido en la circunferencia de un círculo que rueda '''sin deslizar''' sobre una línea recta. Aunque también se puede entender como el lugar geométrico descrito por las posiciones de un punto P que pertenece a la circunferencia de un círculo que rueda sin deslizamiento sobre una recta. Dicha recta recibe el nombre de '''directriz''', mientras que la circunferencia recibe el nombre de '''generatriz o ruleta'''. | ||
| − | A pesar de su sencilla definición geométrica, esta curva tiene propiedades notables que le consiguen un lugar crucial en la física teórica, como en el estudio de los '''principios variacionales''' que cimentan la '''Teoría de Campos'''.Principalmente se conoce por la aplicación como solución | + | A pesar de su sencilla definición geométrica, esta curva tiene propiedades notables que le consiguen un lugar crucial en la física teórica, así como en el estudio de los '''principios variacionales''' que cimentan la '''Teoría de Campos'''.Principalmente se conoce por la aplicación como solución a los problemas de la '''Braquistócrona''' (La curva de descenso más rápido entre dos puntos) y de la '''Tautócrona''' (la curva en que el tiempo de descenso a un punto inferior es independiente de la posición inicial), estableciendo una relación directa con el formalismo de '''Lagrange y Hamilton'''. |
La '''Cicloide''' se puede representar algebraicamente mediante la siguiente parametrización en coordenadas cartesianas: | La '''Cicloide''' se puede representar algebraicamente mediante la siguiente parametrización en coordenadas cartesianas: | ||
| − | <math>𝛾(𝑡) = (𝑥(𝑡), 𝑦(𝑡)) = (𝑅(𝑡 − sin 𝑡), 𝑅(1 − | + | <math>𝛾(𝑡) = (𝑥(𝑡), 𝑦(𝑡)) = (𝑅(𝑡 − \sin 𝑡), 𝑅(1 − \cos 𝑡)), 𝑡 ∈ (0, 2𝜋)</math>. Donde '''R''' es el radio del círculo y '''t''' es el ángulo de '''desplazamiento angular'''. |
| + | |||
| + | A continuación, se realiza un estudio sobre las distintas partes de la Cicloide, que incluye representaciones y cálculos geométricos. El estudio se desarrolla sobre un ejemplo concreto de una Cicloide descrita a partir de la circunferencia de un círculo de radio R=3. Además, también se facilita el link al póster del trabajo para una mejor comprensión: https://drive.google.com/file/d/1JglUCld_10qRH4lqG34cGiU2vzmZ6BZr/view?usp=sharing | ||
==Dibujo de la curva== | ==Dibujo de la curva== | ||
| − | Para este estudio se toma como origen una circunferencia de radio <math>R = 3</math> | + | Para este estudio se toma como origen una circunferencia de radio <math>R = 3</math>. Este primer apartado sirve como base para definir la curva en el software de '''MATLAB'''. Con la parametrización dada se representa la gráfica. |
| + | |||
| + | [[Archivo:Curva cicloide.png|Curva obtenida en MATLAB|right|]] | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| Línea 40: | Línea 46: | ||
==Cálculo de los vectores velocidad 𝛾 ′ (𝑡) y aceleración 𝛾 ″ (𝑡).== | ==Cálculo de los vectores velocidad 𝛾 ′ (𝑡) y aceleración 𝛾 ″ (𝑡).== | ||
| − | + | Una vez está la curva representada con la '''parametrización''' dada, se deriva en función del parámetro 𝑡, que representa el ángulo de giro, para obtener los vectores pedidos. Como se trata de una función muy sencilla, el vector posición se deriva analíticamente mediante la '''regla de la cadena''', obteniendo así la velocidad y aceleración. | |
| + | [[Archivo:Vel acel18.png|marco|derecha|Gráfica de la Cicloide con los vectores de velocidad y aceleración.]] | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | % | + | %Parámetros de entrada |
R = 3; | R = 3; | ||
| − | % | + | %Cálculo de la trayectoria (Para dibujar la curva completa) |
t = linspace(0, 2*pi,200); | t = linspace(0, 2*pi,200); | ||
% Fórmulas directas de la Cicloide | % Fórmulas directas de la Cicloide | ||
| Línea 51: | Línea 58: | ||
y = R * (1 - cos(t)); | y = R * (1 - cos(t)); | ||
| − | % | + | %Cálculo de vectores en puntos específicos |
t_vec = linspace(0,2*pi,15) | t_vec = linspace(0,2*pi,15) | ||
% Posición de los vectores (origen de las flechas) | % Posición de los vectores (origen de las flechas) | ||
| Línea 63: | Línea 70: | ||
ay = R * cos(t_vec); | ay = R * cos(t_vec); | ||
| − | % | + | %Representación gráfica |
figure('Name','La Cicloide', 'Color', 'w'); | figure('Name','La Cicloide', 'Color', 'w'); | ||
hold on; axis equal; grid on; | hold on; axis equal; grid on; | ||
| Línea 82: | Línea 89: | ||
==Longitud de la curva== | ==Longitud de la curva== | ||
| − | Para este cálculo se hará uso del método del rectángulo. ( | + | Para este cálculo se hará uso del método del rectángulo. (El número escogido no es uno en concreto, pero cuanto mayor sea, mejor será la aproximación) |
| + | |||
| + | El cálculo de la longitud de la curva descrita se resuelve a partir del módulo del vector velocidad de la | ||
| + | mism planteando la siguiente integral en un intervalo (a, b). | ||
| + | |||
| + | <math>l(\gamma) = \int_0^{2\pi} \|\gamma'(t)\| dt = \int_0^{2\pi} \sqrt{x'(t)^2 + y'(t)^2} \, dt = \int_0^{2\pi} \sqrt{(R-R\cos t)^2 +(R\sin^2 t)^2 } \, dt </math> | ||
| + | |||
| + | <math>l(\gamma) = \int_0^{2\pi} 3\sqrt{(1-1\cos t)^2 +(\sin^2 t)^2 } \, dt </math> | ||
| + | |||
| + | <math>l(\gamma) = 12\int_0^{2\pi} \frac{1}{2}\sin(\frac{t}{2})\, dt = 24 </math> unidades | ||
| + | |||
| + | Para comprobar la resolución teórica de la integral, se acude al “Método del rectángulo”. Para ello, se utiliza la plataforma MATLAB que realiza el cálculo a partir del siguiente comando: | ||
| + | |||
{{matlab|codigo= | {{matlab|codigo= | ||
| Línea 113: | Línea 132: | ||
==Cálculo de los vectores tangente 𝑡⃗(𝑡) y normal 𝑛⃗(𝑡)== | ==Cálculo de los vectores tangente 𝑡⃗(𝑡) y normal 𝑛⃗(𝑡)== | ||
| + | Para calcular el vector tangente a la curva <math>𝑡⃗(𝑡)</math> haremos uso del vector velocidad calculado anteriormente y lo normalizamos. | ||
| + | |||
| + | <math>\|\gamma'(t)\| = (R(1-\cos t), R\sin t)</math> | ||
| + | |||
| + | Partiendo de este vector, y obteniendo su valor absoluto <math>\|\gamma'(t)\| = 2Rsin\frac{t}{2}</math> | ||
| + | |||
| + | El vector tangente unitario se obtiene como el cociente entre la velocidad y su norma. Operando y simplificando términos resulta tal que: | ||
| + | |||
| + | <math>𝑡⃗(𝑡)= (\sin(\frac{t}{2}) ,\cos(\frac{t}{2}))</math> | ||
| + | |||
| + | Para calcular el vector normal 𝑛⃗(𝑡)como estamos en R2 será suficiente con cambiar el orden del vector tangente y cambiarle el signo a una coordenada. Entre los dos posibles vectores normales se elegirá el interior a la curva por comodidad visual: | ||
| + | <math>𝑛⃗(𝑡)= (\cos(\frac{t}{2}) ,-\sin(\frac{t}{2}))</math> | ||
| + | |||
| + | Una vez se han obtenido los dos vectores iniciaremos el código de MatLab donde representaremos 25 de estos vectores en la cicloide según el parámetro t: | ||
| + | [[Archivo:Tangente normal.png|450px|marco|derecha|Gráfica con la curva y los vectores tangentes y normales en varios puntos.]] | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | + | ||
R = 3; | R = 3; | ||
| − | |||
t_domain = linspace(0.01, 2*pi - 0.01, 200); | t_domain = linspace(0.01, 2*pi - 0.01, 200); | ||
% Dominio para la curva, evitando cúspides | % Dominio para la curva, evitando cúspides | ||
t_vec = linspace(0.1, 2*pi - 0.1, 15); | t_vec = linspace(0.1, 2*pi - 0.1, 15); | ||
| − | % | + | |
| + | %Cálculo de la Cicloide | ||
X = R * (t_domain - sin(t_domain)); | X = R * (t_domain - sin(t_domain)); | ||
Y = R * (1 - cos(t_domain)); | Y = R * (1 - cos(t_domain)); | ||
| − | % | + | |
| − | % Calculamos t/2 solo para los puntos de la muestra | + | %Cálculo de Vectores Tangente y Normal Unitarios |
| − | + | t_half = t_vec / 2; %Calculamos t/2 solo para los puntos de la muestra | |
| − | + | t(t) = (sin(t/2), cos(t/2)) % Vector Tangente Unitario | |
Tx = sin(t_half); | Tx = sin(t_half); | ||
Ty = cos(t_half); | Ty = cos(t_half); | ||
| − | + | n(t) = (-cos(t/2), sin(t/2)) % Vector Normal Unitario | |
Nx = cos(t_half); | Nx = cos(t_half); | ||
Ny = -sin(t_half); | Ny = -sin(t_half); | ||
| − | % | + | |
| + | %Dibujo de la Curva y Vectores | ||
figure('Name','La Cicloide', 'Color', 'w'); | figure('Name','La Cicloide', 'Color', 'w'); | ||
hold on; | hold on; | ||
| − | |||
plot(X, Y, 'r', 'LineWidth', 2); | plot(X, Y, 'r', 'LineWidth', 2); | ||
| + | |||
% Coordenadas de origen para los vectores (en los puntos t_vec) | % Coordenadas de origen para los vectores (en los puntos t_vec) | ||
X_vec = R * (t_vec - sin(t_vec)); | X_vec = R * (t_vec - sin(t_vec)); | ||
Y_vec = R * (1 - cos(t_vec)); | Y_vec = R * (1 - cos(t_vec)); | ||
| + | |||
% Factor de escala visual para que los vectores sean visibles | % Factor de escala visual para que los vectores sean visibles | ||
escala = 0.8; | escala = 0.8; | ||
| − | |||
quiver(X_vec, Y_vec, Tx, Ty, escala, 'b', 'LineWidth', 0.5, 'MaxHeadSize', 0.1); | quiver(X_vec, Y_vec, Tx, Ty, escala, 'b', 'LineWidth', 0.5, 'MaxHeadSize', 0.1); | ||
| − | |||
quiver(X_vec, Y_vec, Nx, Ny, escala, 'g', 'LineWidth', 0.5, 'MaxHeadSize', 0.1); | quiver(X_vec, Y_vec, Nx, Ny, escala, 'g', 'LineWidth', 0.5, 'MaxHeadSize', 0.1); | ||
title(['Vectores Tangente y Normal de la Cicloide (R=', num2str(R), ')']); | title(['Vectores Tangente y Normal de la Cicloide (R=', num2str(R), ')']); | ||
| Línea 156: | Línea 190: | ||
==Cálculo de la curvatura 𝜅(𝑡)== | ==Cálculo de la curvatura 𝜅(𝑡)== | ||
| + | Para el cálculo de la curvatura se aplica directamente la siguiente fórmula: | ||
| + | <math>\kappa(t) = \frac{\gamma'(t)\times\gamma''(t)}{|\gamma'(t)|^3}</math> | ||
| + | [[Archivo:Curvatura cicloide18.png|250px|marco|left|Representación de la curvatura de la cicloide.]] | ||
| + | |||
| + | Obteniendo como resultado: <math>\frac{1}{4R\sin(\frac{t}{2})} = \frac{1}{12\sin(\frac{t}{2})}</math> | ||
| + | |||
| + | Para el código, es imprescindible definir bien las derivadas parciales para posteriormente evitar problemas en la fórmula de la curvatura. | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | + | ||
| − | + | ||
%Parámetrización de la cicloide con R=3 | %Parámetrización de la cicloide con R=3 | ||
| − | x=3*(t-sin(t)); | + | t=linspace(0,2*pi,500);x=3*(t-sin(t)); |
y=3*(1-cos(t)); | y=3*(1-cos(t)); | ||
| + | |||
%Definición de las derivadas parciales primeras respecto a t | %Definición de las derivadas parciales primeras respecto a t | ||
dxdt=3-3*cos(t); | dxdt=3-3*cos(t); | ||
| Línea 169: | Línea 210: | ||
dx2dt2=3*sin(t); | dx2dt2=3*sin(t); | ||
dy2dt2=3*cos(t); | dy2dt2=3*cos(t); | ||
| + | |||
%Cálculo de la curvatura a lo largo de la cicloide | %Cálculo de la curvatura a lo largo de la cicloide | ||
k=abs((dxdt.*dy2dt2)-(dx2dt2.*dydt))./((dxdt.^2+dydt.^2).^(3/2)); | k=abs((dxdt.*dy2dt2)-(dx2dt2.*dydt))./((dxdt.^2+dydt.^2).^(3/2)); | ||
| − | % | + | |
| − | + | %Grafica de la cicloide en subventana | |
plot(t,k,'b','LineWidth',3); | plot(t,k,'b','LineWidth',3); | ||
| − | |||
title('Curvatura de la Cicloide para R=3') | title('Curvatura de la Cicloide para R=3') | ||
xlabel('Parámetro t'); | xlabel('Parámetro t'); | ||
ylabel('Curvatura \kappa(t)'); | ylabel('Curvatura \kappa(t)'); | ||
| − | %Redimensionado del eje x | + | ax=gca; %Redimensionado del eje x |
| − | + | ||
ax.XTick=[0,pi/2,pi,3*pi/2,2*pi]; | ax.XTick=[0,pi/2,pi,3*pi/2,2*pi]; | ||
ax.XTickLabel={'0','\pi/2','\pi','3\pi/2','2\pi'}; | ax.XTickLabel={'0','\pi/2','\pi','3\pi/2','2\pi'}; | ||
| − | |||
axis manual | axis manual | ||
xlim([0,2*pi]) | xlim([0,2*pi]) | ||
ylim([0,3]) | ylim([0,3]) | ||
grid on | grid on | ||
| + | |||
%Graficación de la curvatura en la segunda subventana | %Graficación de la curvatura en la segunda subventana | ||
subplot(2,1,1); | subplot(2,1,1); | ||
plot(x, y, 'r', 'LineWidth', 2); | plot(x, y, 'r', 'LineWidth', 2); | ||
title('Gráfica de la Cicloide con R=3'); | title('Gráfica de la Cicloide con R=3'); | ||
| − | |||
xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex'); | xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex'); | ||
ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex'); | ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex'); | ||
| Línea 200: | Línea 239: | ||
}} | }} | ||
| − | ==Obtención del centro y el radio de la circunferencia osculatriz== | + | ==Obtención del centro y el radio de la circunferencia osculatriz.== |
| + | |||
| + | Del apartado anterior, se conoce la curvatura de la cicloide, y para las operaciones a realizar en este apartado, vamos a eliminar el valor absoluto del numerador de la definición de curvatura para que al momento de obtener el vector normal (t) salga orientado correctamente, quedando así: <math>\kappa(t) = -\frac{1}{12\sin(\frac{t}{2})}</math> | ||
| + | [[Archivo:Osculatriz18.png |250px|marco|izquierda|Gráfico representación del radio de curvatura y la propia circunferencia osculatriz.]] | ||
| + | Dado que el radio de curvatura es la inversa de la curvatura, el cálculo del radio se obtiene particularizando t=4 en <math>R(t) = \frac{1}{|\kappa(t)|} = 10.9116</math> | ||
| + | |||
| + | Para obtener el centro de curvatura se aplica <math>O(t)=\gamma(t) + R(t)*n(t)</math> y se particulariza en el punto P, obteniendo las coordenadas cartesianas del centro de curvatura. | ||
| + | |||
| + | <math>O(t)= 9,7296i-4,9609j</math> | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | t=linspace(0,2*pi,500); % | + | |
| − | x=@(t) 3*(t-sin(t)); | + | t=linspace(0,2*pi,500); |
| − | y=@(t) 3*(1-cos(t)); | + | %Funciones anónimas que definen las coordenadas x(t) e y (t) de la cicloide |
| + | x=@(t) 3*(t-sin(t)); | ||
| + | y=@(t) 3*(1-cos(t)); | ||
t0=4; %Punto a evaluar (t=4) | t0=4; %Punto a evaluar (t=4) | ||
| + | |||
%Funciónes de las derivadas parciales primeras y segundas | %Funciónes de las derivadas parciales primeras y segundas | ||
dxdt=@(t) 3-3*cos(t); dx2dt2=@(t) 3*sin(t); | dxdt=@(t) 3-3*cos(t); dx2dt2=@(t) 3*sin(t); | ||
dydt=@(t) 3*sin(t); dy2dt2=@(t) 3*cos(t); | dydt=@(t) 3*sin(t); dy2dt2=@(t) 3*cos(t); | ||
| + | |||
%Evaluación de la parametrización de la cicloide y sus derivadas primeras y segundas en t=4 | %Evaluación de la parametrización de la cicloide y sus derivadas primeras y segundas en t=4 | ||
x0=x(t0); dx0=dxdt(t0); dx2_0=dx2dt2(t0); | x0=x(t0); dx0=dxdt(t0); dx2_0=dx2dt2(t0); | ||
y0=y(t0); dy0=dydt(t0); dy2_0=dy2dt2(t0); | y0=y(t0); dy0=dydt(t0); dy2_0=dy2dt2(t0); | ||
k=((dx0.*dy2_0)-(dx2_0.*dy0))./((dx0.^2+dy0.^2).^(3/2)); %Curvatura en el punto | k=((dx0.*dy2_0)-(dx2_0.*dy0))./((dx0.^2+dy0.^2).^(3/2)); %Curvatura en el punto | ||
| − | r=1./abs(k); | + | r=1./abs(k); %Radio de la circunferencia osculatriz |
| + | |||
%Coordenadas cartesianas del centro de la circunferencia en t=4 | %Coordenadas cartesianas del centro de la circunferencia en t=4 | ||
xc=3*4-3*sin(t0)+(1/k).*(-sin(t0))/sqrt(2*(1-cos(t0))); | xc=3*4-3*sin(t0)+(1/k).*(-sin(t0))/sqrt(2*(1-cos(t0))); | ||
yc=3-3*cos(t0)+(1/k)*((1-cos(t0)))/sqrt(2*(1-cos(t0))); | yc=3-3*cos(t0)+(1/k)*((1-cos(t0)))/sqrt(2*(1-cos(t0))); | ||
| − | |||
theta=linspace(0,2*pi,500); | theta=linspace(0,2*pi,500); | ||
xcirc=xc+r*cos(theta); | xcirc=xc+r*cos(theta); | ||
ycirc=yc+r*sin(theta); | ycirc=yc+r*sin(theta); | ||
| − | %Generación | + | |
| + | %Generación de la gráfica | ||
hold on; grid on; axis equal; | hold on; grid on; axis equal; | ||
| − | plot(x(t),y(t),'b','LineWidth',3) % | + | plot(x(t),y(t),'b','LineWidth',3) %Cicloide |
| − | plot(xcirc,ycirc,'g--'); % | + | plot(xcirc,ycirc,'g--'); %Circunferencia |
| − | plot(x0,y0,'ro','MarkerSize',8,'MarkerFaceColor','r'); % | + | plot(x0,y0,'ro','MarkerSize',8,'MarkerFaceColor','r'); %Gamma(t=4) |
plot([xc x0],[yc y0],'k:','LineWidth',1) %Radio de la circunferencia | plot([xc x0],[yc y0],'k:','LineWidth',1) %Radio de la circunferencia | ||
plot(xc,yc,'o','Color', [0.5 0 0.5],'MarkerSize',8,'MarkerFaceColor',' [0.5 0 0.5]') %Centro de la circunferencia | plot(xc,yc,'o','Color', [0.5 0 0.5],'MarkerSize',8,'MarkerFaceColor',' [0.5 0 0.5]') %Centro de la circunferencia | ||
| Línea 233: | Línea 285: | ||
hold off | hold off | ||
%Leyenda y etiquetas | %Leyenda y etiquetas | ||
| − | legend('Gráfica de la cicloide \gamma(t) en t\in(0,2\pi)','Circunferencia osculatriz en t=4', ... | + | legend('Gráfica de la cicloide \gamma(t) en t\in(0,2\pi)','Circunferencia osculatriz en t=4', ...'P=\gamma(4)','Radio de curvatura en P=\gamma(4)','Centro de la circunferencia','Location','southwest') |
| − | + | ||
xlabel('Parámetro t'); | xlabel('Parámetro t'); | ||
ylabel('\gamma(t)'); | ylabel('\gamma(t)'); | ||
| Línea 241: | Línea 292: | ||
==Información adicional de la curva y sus aplicaciones== | ==Información adicional de la curva y sus aplicaciones== | ||
| − | El estudio de esta curva reviste una importancia histórica y práctica considerable dentro de la física y la ingeniería, especialmente en el contexto de la mecánica clásica y el diseño de estructuras. | + | El estudio de esta curva reviste una importancia '''histórica y práctica''' considerable dentro de la física y la ingeniería, especialmente en el contexto de la mecánica clásica y el diseño de estructuras. |
| + | <center>[[Archivo:Geom cic.png|650px|sinmarco]]</center> | ||
| + | La relevancia de la cicloide para la ingeniería se fundamenta en dos propiedades geométricas únicas. En primer lugar, la '''cicloide invertida''' es una curva tautócrona (del griego tautos, mismo, y chronos, tiempo), lo que implica que un cuerpo que se desliza sin fricción a lo largo de esta trayectoria, bajo la acción constante de la gravedad, tardará exactamente el mismo tiempo en alcanzar el punto más bajo, independientemente de su punto de partida en la curva. Esta característica de '''isocronismo''' fue crucial en el desarrollo de la ingeniería de precisión, siendo utilizada por '''Christiaan Huygens''' en el siglo XVII para diseñar el reloj de péndulo isócrono al darse cuenta de que el péndulo que sigue una trayectoria circular no es isócrono. | ||
| − | + | En segundo lugar, la cicloide invertida es la solución al denominado problema de la '''braquistócrona''' (brachistos, más corto, y chronos, tiempo), ya que, de todas las curvas posibles que conectan dos puntos, aquella que permite a un cuerpo descender en el menor tiempo posible, bajo la acción de la gravedad y sin rozamiento, es la cicloide. La resolución de este problema por '''Johann Bernoulli''' en 1696 dio origen a una nueva rama de las matemáticas: el Cálculo de Variaciones. | |
| − | En segundo lugar, la cicloide invertida es la solución al denominado problema de la braquistócrona (brachistos, más corto, y chronos, tiempo), ya que, de todas las curvas posibles que conectan dos puntos, aquella que permite a un cuerpo descender en el menor tiempo posible, bajo la acción de la gravedad y sin rozamiento, es la cicloide. La resolución de este problema por Johann Bernoulli en 1696 dio origen a una nueva rama de las matemáticas: el Cálculo de Variaciones. | + | |
| − | Las aplicaciones de la cicloide en ingeniería son diversas: | + | Las aplicaciones de la cicloide en ingeniería son diversas:<br/> |
| − | Diseño de Mecanismos: Se utiliza en la configuración de engranajes con perfil cicloidal para asegurar una transmisión de movimiento suave y eficiente. | + | *'''Diseño de Mecanismos''': Se utiliza en la configuración de engranajes con perfil cicloidal para asegurar una transmisión de movimiento suave y eficiente.<br/> |
| − | Ingeniería Civil: La forma cicloidal se ha empleado en el diseño de arcos y tejados, como en el Museo de Arte Kimbell, para optimizar la distribución de cargas y la estabilidad estructural. Además, sus propiedades son consideradas en el diseño funcional de puentes y túneles para la distribución eficiente de cargas. | + | *'''Ingeniería Civil''': La forma cicloidal se ha empleado en el diseño de arcos y tejados, como en el Museo de Arte Kimbell, para optimizar la distribución de cargas y la estabilidad estructural. Además, sus propiedades son consideradas en el diseño funcional de puentes y túneles para la distribución eficiente de cargas.<br/> |
| − | Transporte y Ocio: La cicloide es la curva utilizada en toboganes de parques y piscinas, así como en pistas de skate o de salto en la nieve, precisamente para garantizar que se alcanza el punto más bajo en el menor tiempo. | + | *'''Transporte y Ocio''': La cicloide es la curva utilizada en toboganes de parques y piscinas, así como en pistas de skate o de salto en la nieve, precisamente para garantizar que se alcanza el punto más bajo en el menor tiempo.<br/> |
| − | Metrología: Una curva cicloide acortada o alargada se emplea en el diseño de tubos manométricos para medir la velocidad del viento con sensibilidad constante en túneles aerodinámicos. | + | *'''Metrología''': Una curva cicloide acortada o alargada se emplea en el diseño de tubos manométricos para medir la velocidad del viento con sensibilidad constante en túneles aerodinámicos.<br/> |
==Ejemplos de uso== | ==Ejemplos de uso== | ||
| + | |||
| + | Su uso '''no es muy frecuente''', pero cuando aparece, se usa por razones técnicas reales, no solo estéticas, ya que permite una distribución eficiente de esfuerzos por '''compresión'''. Esto la hace adecuada para su empleo en el diseño de arcos y bóvedas. Su curvatura progresiva permite reducir los empujes horizontales y optimizar el espesor de '''elementos portantes'''. Además, gracias a la suavidad de su curvatura, la Cicloide también se utiliza en cubiertas que buscan controlar la iluminación natural de forma uniforme. Por otra parte, es característica la presencia de esta curva en el diseño de puentes al ofrecer soluciones estables, elegantes y funcionales.<br> | ||
| + | |||
| + | *Bóvedas y cubiertas cicloidales: '''Kimbell Art Museum''' (EE.UU.) | ||
| + | [[Archivo:Kimbell3.png|550px|sinmarco|izquierda]] | ||
| + | [[Archivo:Dentro kimbell.png|400px|sinmarco|centro]] | ||
| + | [[Archivo:Kimbell1.png|500px|sinmarco|derecha]] <br> | ||
| + | |||
| + | *Puentes con arco ''aproximado'' a la cicloide (pseudo-cicloidales):<br> | ||
| + | [[Archivo:Puente181.png|500px|sinmarco|izquierda]] | ||
| + | [[Archivo:Puente18.png|700px|sinmarco|derecha]] <br> | ||
| + | |||
| + | *También tiene aplicaciones en '''estructuras dinámicas'''. | ||
==Superficie reglada asociada a la curva== | ==Superficie reglada asociada a la curva== | ||
| + | |||
| + | Se puede ver la cicloide en <math>ℝ3</math> mediante la '''parametrización en cartesianas''' <math>𝛾(𝑡) = (𝑥1 (𝑡), 𝑥2 (𝑡), 𝑥3 (𝑡)) = (0, 𝑅(𝑡 − \sin 𝑡), 𝑅(1 + \cos 𝑡)), 𝑡 ∈ (0, 2𝜋)</math>. | ||
| + | |||
| + | Al cambiar la coordenada 𝑅(1−cos(𝑡)) de signo y sumar 2𝑅, es como considerar la curva ''simétrica'' respecto del eje 𝑥3 y subirla 2𝑅 unidades. | ||
| + | |||
| + | [[Archivo:Superficie18.png|marco|derecha|Superficie reglada de la cicloide]] | ||
{{matlab|codigo= | {{matlab|codigo= | ||
| + | |||
R = 3; | R = 3; | ||
L = 1; | L = 1; | ||
| − | % | + | |
| + | %Mallado | ||
t = linspace(0,2*pi); | t = linspace(0,2*pi); | ||
u = linspace(0, L, 10); | u = linspace(0, L, 10); | ||
[T, U] = meshgrid(t, u); | [T, U] = meshgrid(t, u); | ||
| − | % | + | |
| + | %Ecuaciones paramétricas | ||
X = U; | X = U; | ||
Y = R * (T - sin(T)); | Y = R * (T - sin(T)); | ||
Z = R * (1 + cos(T)); | Z = R * (1 + cos(T)); | ||
| − | % | + | |
| + | %Representación gráfica | ||
figure('Name', 'Superficie Reglada con L=1', 'Color', 'w'); | figure('Name', 'Superficie Reglada con L=1', 'Color', 'w'); | ||
| − | |||
| − | |||
| − | |||
| − | |||
surf(X, Y, Z, 'FaceColor', 'cyan', 'EdgeColor', 'k', 'FaceAlpha', 0.7); | surf(X, Y, Z, 'FaceColor', 'cyan', 'EdgeColor', 'k', 'FaceAlpha', 0.7); | ||
| − | |||
axis equal; | axis equal; | ||
grid on; | grid on; | ||
| Línea 283: | Línea 354: | ||
zlim([0, 2*R + 1]); | zlim([0, 2*R + 1]); | ||
}} | }} | ||
| + | |||
| + | En términos de ingeniería civil y arquitectura, a esta superficie se denomina '''Bóveda de Cañón de Directriz Cicloidal'''. Es muy frecuente usarla, ya que puede cubrir grandes '''luces''' con un arranque vertical, optimizando el volumen interior frente a arcos parabólicos o circulares. y entre algunos de los ejemplos de esta superficie en algunas estructuras tenemos:<br> | ||
| + | |||
| + | *El '''Kimbell Art Museum''', considerado como referencia de uso de la cicloide para resolver cubiertas. Ingeniada por August Komendant y diseñada por el arquitecto Louis Kahn. La solución de la cicloide como cubierta tiene una curvatura leve en la parte superior (Clave). Además, el ingeniero Komendant la diseñó de tal forma que se comporta como una viga de hormigón pretensado gigante. | ||
| + | [[Archivo:Kimbell 2.jpeg|500px|sinmarco]] | ||
| + | [[Archivo:Kimbell 1.jpg|450px|sinmarco]]<br> | ||
| + | *'''Hangares y bodegas''':Para hangares de aviones, estaciones de tren antiguas, bodegas subterráneas… Como se usan geometrías que, si bien a menudo se simplifican como arcos elípticos o parabólicos, se asemejan mucho a la cicloide por una razón geométrica: maximizar el volumen útil. Los camiones, trenes o estanterías se peguen a la pared sin perder espacio, y transmiten la carga verticalmente a la cimentación, reduciendo el empuje horizontal.<br> | ||
| + | *Uso industrial: El uso propio de la cicloide no es muy común en la ingeniería, ya que se suelen emplear superficies y curvas parecidas, aunque no concretamente una cicloide. Entre ellas están la catenaria (muy frecuente en los ferrocarriles) y arcos de elipses y parábolas. | ||
==Distribución de la densidad en la superficie== | ==Distribución de la densidad en la superficie== | ||
| + | |||
| + | Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: <math>𝑓(𝑥1, 𝑥2, 𝑥3) = (1 + 𝑥1 )(1 + 𝑥2 )𝑥3</math> | ||
| + | <br> | ||
| + | |||
| + | *Parametrización de la superficie: <center><math>S(t,u)=(u, R(t-sen t), R(1+cos t))</math> <br></center> | ||
| + | *Densidad de la superficie: <center> <math>𝑓(𝑥1,𝑥2,𝑥3) = (1 + 𝑥1)(1 + 𝑥2)𝑥3</math> <br> </center> | ||
| + | *Cálculo de la masa: <center><math>M=\iint ρ(Φ(u,v))||\vec r_u × \vec r_u||\,du\,dv</math></center><br> | ||
| + | |||
| + | Operando y aplicando el teorema de Fubini se obtiene un resultado de <math>M=\iint ρ(Φ(u,v))||\vec r_u × \vec r_u||\,du\,dv = \int_{0}^{1}\int_{0}^{2π}(1+u)(1+R(t-\sin t))(R(1+\cos t))R\sqrt{(2-2*\cos t)} \ dtdu</math> | ||
| + | |||
| + | <math>M=749.83</math> unidades de forma aproximada. | ||
| + | |||
| + | Como se trata de una integral compleja de resolver a mano, se recurre al método del rectángulo empleado en '''MATLAB'''. | ||
| + | |||
{{matlab|codigo= | {{matlab|codigo= | ||
| − | + | ||
R=3; | R=3; | ||
| + | |||
%Definimos la variable de la masa aproximada | %Definimos la variable de la masa aproximada | ||
Masa_aprox=0; | Masa_aprox=0; | ||
| − | |||
n=1000; | n=1000; | ||
t=linspace(0,2*pi,n); | t=linspace(0,2*pi,n); | ||
| + | |||
%Función (integral) que tenemos que resolver | %Función (integral) que tenemos que resolver | ||
f=@(t) (3/2)*((R^2)+(R^3)*(t-sin(t)))*(1+cos(t)).*sqrt(2-2.*cos(t)); | f=@(t) (3/2)*((R^2)+(R^3)*(t-sin(t)))*(1+cos(t)).*sqrt(2-2.*cos(t)); | ||
| + | |||
%metodo del rectángulo | %metodo del rectángulo | ||
for i=1:n; | for i=1:n; | ||
| Línea 302: | Línea 397: | ||
fprintf('M = %.2f',Masa_aprox); | fprintf('M = %.2f',Masa_aprox); | ||
}} | }} | ||
| + | |||
| + | ==Bibliografía y Referencias== | ||
| + | *https://dibujotecni.com/geometria-plana/curvas-ciclicas/ | ||
| + | *https://es.mathworks.com/help/matlab/math/integration-to-find-arc-length.html | ||
| + | *Corcho Gutiérrez, F. M. (2017). La cicloide [Trabajo Fin de Grado, Universidad de Sevilla]. idUS.https://idus.us.es/items/e2574a35-55ee-4102-9602-62eee5d5ea33 | ||
| + | *https://interiordesign.net/designwire/kimbell-art-museum-rare-works-by-louis-kahn/ | ||
| + | *https://ocw.unican.es/pluginfile.php/1768/course/section/1283/MCP6-Fubini-w.pdf | ||
| + | *Barceló, Juan Antonio, Castro, Carlos, González, David y Schiavone, Nico Michele. (2025). Tema_2_Curso_25_26.pdf. En Teoría de Campos, curso académico 2025–2026. ETSI Caminos, Universidad Politécnica de Madrid. [Documento en línea]. Recuperado el 5 de diciembre de 2025, de moodle.upm.es: https://moodle.upm.es/titulaciones/oficiales/pluginfile.php/13654288/mod_resource/content/9/Tema_2_Curso_25_26.pdf. © 2025 Juan Antonio Barceló, Carlos Castro, David González, Nico Michele Schiavone. | ||
| + | *Google, "Gemini" (Modelo de lenguaje grande), versión 1.5, 2025. [Herramienta de IA]. Disponible en: https://gemini.google.com. [Accedido: Diciembre 2025]. | ||
[[Categoría:TC25/26]] | [[Categoría:TC25/26]] | ||
| + | [[Categoría:Teoría de Campos]] | ||
Revisión actual del 22:43 6 dic 2025
| Trabajo realizado por estudiantes | |
|---|---|
| Título | La Cicloide (Grupo 18) |
| Asignatura | Teoría de Campos |
| Curso | 2025-26 |
| Autores | Alejandro Porrúa Perea Salvador Sánchez Burgos Natalia Andrés Jiménez Jaime Colomina López John Cuenca Uyaguari |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
- 1 Introducción
- 2 Dibujo de la curva
- 3 Cálculo de los vectores velocidad 𝛾 ′ (𝑡) y aceleración 𝛾 ″ (𝑡).
- 4 Longitud de la curva
- 5 Cálculo de los vectores tangente 𝑡⃗(𝑡) y normal 𝑛⃗(𝑡)
- 6 Cálculo de la curvatura 𝜅(𝑡)
- 7 Obtención del centro y el radio de la circunferencia osculatriz.
- 8 Información adicional de la curva y sus aplicaciones
- 9 Ejemplos de uso
- 10 Superficie reglada asociada a la curva
- 11 Distribución de la densidad en la superficie
- 12 Bibliografía y Referencias
1 Introducción
La Cicloide se define como la curva trazada por un punto contenido en la circunferencia de un círculo que rueda sin deslizar sobre una línea recta. Aunque también se puede entender como el lugar geométrico descrito por las posiciones de un punto P que pertenece a la circunferencia de un círculo que rueda sin deslizamiento sobre una recta. Dicha recta recibe el nombre de directriz, mientras que la circunferencia recibe el nombre de generatriz o ruleta.
A pesar de su sencilla definición geométrica, esta curva tiene propiedades notables que le consiguen un lugar crucial en la física teórica, así como en el estudio de los principios variacionales que cimentan la Teoría de Campos.Principalmente se conoce por la aplicación como solución a los problemas de la Braquistócrona (La curva de descenso más rápido entre dos puntos) y de la Tautócrona (la curva en que el tiempo de descenso a un punto inferior es independiente de la posición inicial), estableciendo una relación directa con el formalismo de Lagrange y Hamilton.
La Cicloide se puede representar algebraicamente mediante la siguiente parametrización en coordenadas cartesianas: [math]𝛾(𝑡) = (𝑥(𝑡), 𝑦(𝑡)) = (𝑅(𝑡 − \sin 𝑡), 𝑅(1 − \cos 𝑡)), 𝑡 ∈ (0, 2𝜋)[/math]. Donde R es el radio del círculo y t es el ángulo de desplazamiento angular.
A continuación, se realiza un estudio sobre las distintas partes de la Cicloide, que incluye representaciones y cálculos geométricos. El estudio se desarrolla sobre un ejemplo concreto de una Cicloide descrita a partir de la circunferencia de un círculo de radio R=3. Además, también se facilita el link al póster del trabajo para una mejor comprensión: https://drive.google.com/file/d/1JglUCld_10qRH4lqG34cGiU2vzmZ6BZr/view?usp=sharing
2 Dibujo de la curva
Para este estudio se toma como origen una circunferencia de radio [math]R = 3[/math]. Este primer apartado sirve como base para definir la curva en el software de MATLAB. Con la parametrización dada se representa la gráfica.
%Radio de la circunferencia
R=3;
%Valores entre los que se encuentra (t)
t=linspace(0,2*pi);
%Parametrización en coordenadas cartesianas
x=R*(t-sin(t));
y=R*(1-cos(t));
%Representación de la Cicloide
figure
hold on
plot(x,y, 'r', 'LineWidth',2);
title('Cicloide de R=3');
grid on
xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex');
ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex');
%Tamaño del recuadro en el que estará representada la curva
axis([-1,max(x)+1,0,max(y)+1])
axis equal
hold off
3 Cálculo de los vectores velocidad 𝛾 ′ (𝑡) y aceleración 𝛾 ″ (𝑡).
Una vez está la curva representada con la parametrización dada, se deriva en función del parámetro 𝑡, que representa el ángulo de giro, para obtener los vectores pedidos. Como se trata de una función muy sencilla, el vector posición se deriva analíticamente mediante la regla de la cadena, obteniendo así la velocidad y aceleración.
%Parámetros de entrada
R = 3;
%Cálculo de la trayectoria (Para dibujar la curva completa)
t = linspace(0, 2*pi,200);
% Fórmulas directas de la Cicloide
x = R * (t - sin(t));
y = R * (1 - cos(t));
%Cálculo de vectores en puntos específicos
t_vec = linspace(0,2*pi,15)
% Posición de los vectores (origen de las flechas)
x_vec = R * (t_vec - sin(t_vec));
y_vec = R * (1 - cos(t_vec));
% Velocidad v(t) = r'(t) (Fórmula derivada a mano)
vx = R * (1 - cos(t_vec));
vy = R * sin(t_vec);
% Aceleración a(t) = v'(t) (Fórmula derivada a mano)
ax = R * sin(t_vec);
ay = R * cos(t_vec);
%Representación gráfica
figure('Name','La Cicloide', 'Color', 'w');
hold on; axis equal; grid on;
xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex');
ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex');
title('Cinemática de la Cicloide para R = 3',['r' ...'']);
% Dibujo de la trayectoria de la curva
plot(x, y, 'r', 'LineWidth', 2, 'DisplayName', 'Trayectoria');
% Dibujo de los vectores velocidad (Azul) y aceleración (Rojo)
quiver(x_vec, y_vec, vx, vy, 0.5, 'm', 'LineWidth', 0.5, 'maxheadsize',0.1, 'DisplayName', 'Velocidad');
quiver(x_vec, y_vec, ax, ay, 0.5, 'c', 'LineWidth', 0.5, 'maxheadsize',0.1,'DisplayName', 'Aceleración');
legend('Location', 'best');
xlim([0, R*2*pi + 1]);
ylim([0, 2*R + 2]);
hold off;
4 Longitud de la curva
Para este cálculo se hará uso del método del rectángulo. (El número escogido no es uno en concreto, pero cuanto mayor sea, mejor será la aproximación)
El cálculo de la longitud de la curva descrita se resuelve a partir del módulo del vector velocidad de la mism planteando la siguiente integral en un intervalo (a, b).
[math]l(\gamma) = \int_0^{2\pi} \|\gamma'(t)\| dt = \int_0^{2\pi} \sqrt{x'(t)^2 + y'(t)^2} \, dt = \int_0^{2\pi} \sqrt{(R-R\cos t)^2 +(R\sin^2 t)^2 } \, dt [/math]
[math]l(\gamma) = \int_0^{2\pi} 3\sqrt{(1-1\cos t)^2 +(\sin^2 t)^2 } \, dt [/math]
[math]l(\gamma) = 12\int_0^{2\pi} \frac{1}{2}\sin(\frac{t}{2})\, dt = 24 [/math] unidades
Para comprobar la resolución teórica de la integral, se acude al “Método del rectángulo”. Para ello, se utiliza la plataforma MATLAB que realiza el cálculo a partir del siguiente comando:
R=3;
%Definición de la parametrización derivada
dx=@(t) R*(1-cos(t));
dy=@(t) R*sin(t);
%Definición del intervalo del parámetro
a=0;
b=2*pi;
%Definición del número de rectángulos
n=4000;
%División del intervalo entre el número de rectángulos escogido
dt=(b-a)/n;
%Obtención de vectores que unen los puntos medios de cada subintervalo
t_mid=a+dt/2:dt:b-dt/2;
%Aproximación de la longitud de la curva generada a partir de los vectores anteriores
L=sum(sqrt(dx(t_mid).^2+dy(t_mid).^2))*dt;
%Muestra del resultado obtenido
fprintf("La longitud aproximada de la Cicloide es = %.6f\n",L);
5 Cálculo de los vectores tangente 𝑡⃗(𝑡) y normal 𝑛⃗(𝑡)
Para calcular el vector tangente a la curva [math]𝑡⃗(𝑡)[/math] haremos uso del vector velocidad calculado anteriormente y lo normalizamos.
[math]\|\gamma'(t)\| = (R(1-\cos t), R\sin t)[/math]
Partiendo de este vector, y obteniendo su valor absoluto [math]\|\gamma'(t)\| = 2Rsin\frac{t}{2}[/math]
El vector tangente unitario se obtiene como el cociente entre la velocidad y su norma. Operando y simplificando términos resulta tal que:
[math]𝑡⃗(𝑡)= (\sin(\frac{t}{2}) ,\cos(\frac{t}{2}))[/math]
Para calcular el vector normal 𝑛⃗(𝑡)como estamos en R2 será suficiente con cambiar el orden del vector tangente y cambiarle el signo a una coordenada. Entre los dos posibles vectores normales se elegirá el interior a la curva por comodidad visual: [math]𝑛⃗(𝑡)= (\cos(\frac{t}{2}) ,-\sin(\frac{t}{2}))[/math]
Una vez se han obtenido los dos vectores iniciaremos el código de MatLab donde representaremos 25 de estos vectores en la cicloide según el parámetro t:
R = 3;
t_domain = linspace(0.01, 2*pi - 0.01, 200);
% Dominio para la curva, evitando cúspides
t_vec = linspace(0.1, 2*pi - 0.1, 15);
%Cálculo de la Cicloide
X = R * (t_domain - sin(t_domain));
Y = R * (1 - cos(t_domain));
%Cálculo de Vectores Tangente y Normal Unitarios
t_half = t_vec / 2; %Calculamos t/2 solo para los puntos de la muestra
t(t) = (sin(t/2), cos(t/2)) % Vector Tangente Unitario
Tx = sin(t_half);
Ty = cos(t_half);
n(t) = (-cos(t/2), sin(t/2)) % Vector Normal Unitario
Nx = cos(t_half);
Ny = -sin(t_half);
%Dibujo de la Curva y Vectores
figure('Name','La Cicloide', 'Color', 'w');
hold on;
plot(X, Y, 'r', 'LineWidth', 2);
% Coordenadas de origen para los vectores (en los puntos t_vec)
X_vec = R * (t_vec - sin(t_vec));
Y_vec = R * (1 - cos(t_vec));
% Factor de escala visual para que los vectores sean visibles
escala = 0.8;
quiver(X_vec, Y_vec, Tx, Ty, escala, 'b', 'LineWidth', 0.5, 'MaxHeadSize', 0.1);
quiver(X_vec, Y_vec, Nx, Ny, escala, 'g', 'LineWidth', 0.5, 'MaxHeadSize', 0.1);
title(['Vectores Tangente y Normal de la Cicloide (R=', num2str(R), ')']);
xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex');
ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex');
legend('Cicloide', 'Tangente $\vec{t}(t)$', 'Normal $\vec{n}(t)$', 'Location', 'NorthEast', 'Interpreter', 'latex');
axis equal;
grid on;
hold off;
6 Cálculo de la curvatura 𝜅(𝑡)
Para el cálculo de la curvatura se aplica directamente la siguiente fórmula: [math]\kappa(t) = \frac{\gamma'(t)\times\gamma''(t)}{|\gamma'(t)|^3}[/math]
Obteniendo como resultado: [math]\frac{1}{4R\sin(\frac{t}{2})} = \frac{1}{12\sin(\frac{t}{2})}[/math]
Para el código, es imprescindible definir bien las derivadas parciales para posteriormente evitar problemas en la fórmula de la curvatura.
%Parámetrización de la cicloide con R=3
t=linspace(0,2*pi,500);x=3*(t-sin(t));
y=3*(1-cos(t));
%Definición de las derivadas parciales primeras respecto a t
dxdt=3-3*cos(t);
dydt=3*sin(t);
%Definición de las derivadas parciales segundas respecto a t
dx2dt2=3*sin(t);
dy2dt2=3*cos(t);
%Cálculo de la curvatura a lo largo de la cicloide
k=abs((dxdt.*dy2dt2)-(dx2dt2.*dydt))./((dxdt.^2+dydt.^2).^(3/2));
%Grafica de la cicloide en subventana
plot(t,k,'b','LineWidth',3);
title('Curvatura de la Cicloide para R=3')
xlabel('Parámetro t');
ylabel('Curvatura \kappa(t)');
ax=gca; %Redimensionado del eje x
ax.XTick=[0,pi/2,pi,3*pi/2,2*pi];
ax.XTickLabel={'0','\pi/2','\pi','3\pi/2','2\pi'};
axis manual
xlim([0,2*pi])
ylim([0,3])
grid on
%Graficación de la curvatura en la segunda subventana
subplot(2,1,1);
plot(x, y, 'r', 'LineWidth', 2);
title('Gráfica de la Cicloide con R=3');
xlabel('$x(t) = 3(t - \sin t)$', 'Interpreter', 'latex');
ylabel('$y(t) = 3(1 - \cos t)$', 'Interpreter', 'latex');
xlim([0,2*pi]);
ylim('auto');
axis equal;
grid on;
7 Obtención del centro y el radio de la circunferencia osculatriz.
Del apartado anterior, se conoce la curvatura de la cicloide, y para las operaciones a realizar en este apartado, vamos a eliminar el valor absoluto del numerador de la definición de curvatura para que al momento de obtener el vector normal (t) salga orientado correctamente, quedando así: [math]\kappa(t) = -\frac{1}{12\sin(\frac{t}{2})}[/math]
Dado que el radio de curvatura es la inversa de la curvatura, el cálculo del radio se obtiene particularizando t=4 en [math]R(t) = \frac{1}{|\kappa(t)|} = 10.9116[/math]
Para obtener el centro de curvatura se aplica [math]O(t)=\gamma(t) + R(t)*n(t)[/math] y se particulariza en el punto P, obteniendo las coordenadas cartesianas del centro de curvatura.
[math]O(t)= 9,7296i-4,9609j[/math]
t=linspace(0,2*pi,500);
%Funciones anónimas que definen las coordenadas x(t) e y (t) de la cicloide
x=@(t) 3*(t-sin(t));
y=@(t) 3*(1-cos(t));
t0=4; %Punto a evaluar (t=4)
%Funciónes de las derivadas parciales primeras y segundas
dxdt=@(t) 3-3*cos(t); dx2dt2=@(t) 3*sin(t);
dydt=@(t) 3*sin(t); dy2dt2=@(t) 3*cos(t);
%Evaluación de la parametrización de la cicloide y sus derivadas primeras y segundas en t=4
x0=x(t0); dx0=dxdt(t0); dx2_0=dx2dt2(t0);
y0=y(t0); dy0=dydt(t0); dy2_0=dy2dt2(t0);
k=((dx0.*dy2_0)-(dx2_0.*dy0))./((dx0.^2+dy0.^2).^(3/2)); %Curvatura en el punto
r=1./abs(k); %Radio de la circunferencia osculatriz
%Coordenadas cartesianas del centro de la circunferencia en t=4
xc=3*4-3*sin(t0)+(1/k).*(-sin(t0))/sqrt(2*(1-cos(t0)));
yc=3-3*cos(t0)+(1/k)*((1-cos(t0)))/sqrt(2*(1-cos(t0)));
theta=linspace(0,2*pi,500);
xcirc=xc+r*cos(theta);
ycirc=yc+r*sin(theta);
%Generación de la gráfica
hold on; grid on; axis equal;
plot(x(t),y(t),'b','LineWidth',3) %Cicloide
plot(xcirc,ycirc,'g--'); %Circunferencia
plot(x0,y0,'ro','MarkerSize',8,'MarkerFaceColor','r'); %Gamma(t=4)
plot([xc x0],[yc y0],'k:','LineWidth',1) %Radio de la circunferencia
plot(xc,yc,'o','Color', [0.5 0 0.5],'MarkerSize',8,'MarkerFaceColor',' [0.5 0 0.5]') %Centro de la circunferencia
ylim([-15,8])
yticks(-15:5:8)
hold off
%Leyenda y etiquetas
legend('Gráfica de la cicloide \gamma(t) en t\in(0,2\pi)','Circunferencia osculatriz en t=4', ...'P=\gamma(4)','Radio de curvatura en P=\gamma(4)','Centro de la circunferencia','Location','southwest')
xlabel('Parámetro t');
ylabel('\gamma(t)');
8 Información adicional de la curva y sus aplicaciones
El estudio de esta curva reviste una importancia histórica y práctica considerable dentro de la física y la ingeniería, especialmente en el contexto de la mecánica clásica y el diseño de estructuras.
La relevancia de la cicloide para la ingeniería se fundamenta en dos propiedades geométricas únicas. En primer lugar, la cicloide invertida es una curva tautócrona (del griego tautos, mismo, y chronos, tiempo), lo que implica que un cuerpo que se desliza sin fricción a lo largo de esta trayectoria, bajo la acción constante de la gravedad, tardará exactamente el mismo tiempo en alcanzar el punto más bajo, independientemente de su punto de partida en la curva. Esta característica de isocronismo fue crucial en el desarrollo de la ingeniería de precisión, siendo utilizada por Christiaan Huygens en el siglo XVII para diseñar el reloj de péndulo isócrono al darse cuenta de que el péndulo que sigue una trayectoria circular no es isócrono.
En segundo lugar, la cicloide invertida es la solución al denominado problema de la braquistócrona (brachistos, más corto, y chronos, tiempo), ya que, de todas las curvas posibles que conectan dos puntos, aquella que permite a un cuerpo descender en el menor tiempo posible, bajo la acción de la gravedad y sin rozamiento, es la cicloide. La resolución de este problema por Johann Bernoulli en 1696 dio origen a una nueva rama de las matemáticas: el Cálculo de Variaciones.
Las aplicaciones de la cicloide en ingeniería son diversas:
- Diseño de Mecanismos: Se utiliza en la configuración de engranajes con perfil cicloidal para asegurar una transmisión de movimiento suave y eficiente.
- Ingeniería Civil: La forma cicloidal se ha empleado en el diseño de arcos y tejados, como en el Museo de Arte Kimbell, para optimizar la distribución de cargas y la estabilidad estructural. Además, sus propiedades son consideradas en el diseño funcional de puentes y túneles para la distribución eficiente de cargas.
- Transporte y Ocio: La cicloide es la curva utilizada en toboganes de parques y piscinas, así como en pistas de skate o de salto en la nieve, precisamente para garantizar que se alcanza el punto más bajo en el menor tiempo.
- Metrología: Una curva cicloide acortada o alargada se emplea en el diseño de tubos manométricos para medir la velocidad del viento con sensibilidad constante en túneles aerodinámicos.
9 Ejemplos de uso
Su uso no es muy frecuente, pero cuando aparece, se usa por razones técnicas reales, no solo estéticas, ya que permite una distribución eficiente de esfuerzos por compresión. Esto la hace adecuada para su empleo en el diseño de arcos y bóvedas. Su curvatura progresiva permite reducir los empujes horizontales y optimizar el espesor de elementos portantes. Además, gracias a la suavidad de su curvatura, la Cicloide también se utiliza en cubiertas que buscan controlar la iluminación natural de forma uniforme. Por otra parte, es característica la presencia de esta curva en el diseño de puentes al ofrecer soluciones estables, elegantes y funcionales.
- Bóvedas y cubiertas cicloidales: Kimbell Art Museum (EE.UU.)
- Puentes con arco aproximado a la cicloide (pseudo-cicloidales):
- También tiene aplicaciones en estructuras dinámicas.
10 Superficie reglada asociada a la curva
Se puede ver la cicloide en [math]ℝ3[/math] mediante la parametrización en cartesianas [math]𝛾(𝑡) = (𝑥1 (𝑡), 𝑥2 (𝑡), 𝑥3 (𝑡)) = (0, 𝑅(𝑡 − \sin 𝑡), 𝑅(1 + \cos 𝑡)), 𝑡 ∈ (0, 2𝜋)[/math].
Al cambiar la coordenada 𝑅(1−cos(𝑡)) de signo y sumar 2𝑅, es como considerar la curva simétrica respecto del eje 𝑥3 y subirla 2𝑅 unidades.
R = 3;
L = 1;
%Mallado
t = linspace(0,2*pi);
u = linspace(0, L, 10);
[T, U] = meshgrid(t, u);
%Ecuaciones paramétricas
X = U;
Y = R * (T - sin(T));
Z = R * (1 + cos(T));
%Representación gráfica
figure('Name', 'Superficie Reglada con L=1', 'Color', 'w');
surf(X, Y, Z, 'FaceColor', 'cyan', 'EdgeColor', 'k', 'FaceAlpha', 0.7);
axis equal;
grid on;
xlabel('X [m] (Generatriz)'); ylabel('Y [m] (Directriz)'); zlabel('Z [m] (Altura)');
title('Bóveda Cicloidal (R = 3, L = 1)');
view(3);
xlim([0, 2]);
ylim([-1, 2*pi*R + 1]);
zlim([0, 2*R + 1]);
En términos de ingeniería civil y arquitectura, a esta superficie se denomina Bóveda de Cañón de Directriz Cicloidal. Es muy frecuente usarla, ya que puede cubrir grandes luces con un arranque vertical, optimizando el volumen interior frente a arcos parabólicos o circulares. y entre algunos de los ejemplos de esta superficie en algunas estructuras tenemos:
- El Kimbell Art Museum, considerado como referencia de uso de la cicloide para resolver cubiertas. Ingeniada por August Komendant y diseñada por el arquitecto Louis Kahn. La solución de la cicloide como cubierta tiene una curvatura leve en la parte superior (Clave). Además, el ingeniero Komendant la diseñó de tal forma que se comporta como una viga de hormigón pretensado gigante.
- Hangares y bodegas:Para hangares de aviones, estaciones de tren antiguas, bodegas subterráneas… Como se usan geometrías que, si bien a menudo se simplifican como arcos elípticos o parabólicos, se asemejan mucho a la cicloide por una razón geométrica: maximizar el volumen útil. Los camiones, trenes o estanterías se peguen a la pared sin perder espacio, y transmiten la carga verticalmente a la cimentación, reduciendo el empuje horizontal.
- Uso industrial: El uso propio de la cicloide no es muy común en la ingeniería, ya que se suelen emplear superficies y curvas parecidas, aunque no concretamente una cicloide. Entre ellas están la catenaria (muy frecuente en los ferrocarriles) y arcos de elipses y parábolas.
11 Distribución de la densidad en la superficie
Supongamos que la densidad de la superficie del apartado anterior viene dada por la función: [math]𝑓(𝑥1, 𝑥2, 𝑥3) = (1 + 𝑥1 )(1 + 𝑥2 )𝑥3[/math]
- Parametrización de la superficie:
[math]S(t,u)=(u, R(t-sen t), R(1+cos t))[/math] - Densidad de la superficie:
[math]𝑓(𝑥1,𝑥2,𝑥3) = (1 + 𝑥1)(1 + 𝑥2)𝑥3[/math]
- Cálculo de la masa:
[math]M=\iint ρ(Φ(u,v))||\vec r_u × \vec r_u||\,du\,dv[/math]
Operando y aplicando el teorema de Fubini se obtiene un resultado de [math]M=\iint ρ(Φ(u,v))||\vec r_u × \vec r_u||\,du\,dv = \int_{0}^{1}\int_{0}^{2π}(1+u)(1+R(t-\sin t))(R(1+\cos t))R\sqrt{(2-2*\cos t)} \ dtdu[/math]
[math]M=749.83[/math] unidades de forma aproximada.
Como se trata de una integral compleja de resolver a mano, se recurre al método del rectángulo empleado en MATLAB.
R=3;
%Definimos la variable de la masa aproximada
Masa_aprox=0;
n=1000;
t=linspace(0,2*pi,n);
%Función (integral) que tenemos que resolver
f=@(t) (3/2)*((R^2)+(R^3)*(t-sin(t)))*(1+cos(t)).*sqrt(2-2.*cos(t));
%metodo del rectángulo
for i=1:n;
Masa_aprox=Masa_aprox+((2*pi-0)/n)*f(t(i));
end
fprintf('M = %.2f',Masa_aprox);
12 Bibliografía y Referencias
- https://dibujotecni.com/geometria-plana/curvas-ciclicas/
- https://es.mathworks.com/help/matlab/math/integration-to-find-arc-length.html
- Corcho Gutiérrez, F. M. (2017). La cicloide [Trabajo Fin de Grado, Universidad de Sevilla]. idUS.https://idus.us.es/items/e2574a35-55ee-4102-9602-62eee5d5ea33
- https://interiordesign.net/designwire/kimbell-art-museum-rare-works-by-louis-kahn/
- https://ocw.unican.es/pluginfile.php/1768/course/section/1283/MCP6-Fubini-w.pdf
- Barceló, Juan Antonio, Castro, Carlos, González, David y Schiavone, Nico Michele. (2025). Tema_2_Curso_25_26.pdf. En Teoría de Campos, curso académico 2025–2026. ETSI Caminos, Universidad Politécnica de Madrid. [Documento en línea]. Recuperado el 5 de diciembre de 2025, de moodle.upm.es: https://moodle.upm.es/titulaciones/oficiales/pluginfile.php/13654288/mod_resource/content/9/Tema_2_Curso_25_26.pdf. © 2025 Juan Antonio Barceló, Carlos Castro, David González, Nico Michele Schiavone.
- Google, "Gemini" (Modelo de lenguaje grande), versión 1.5, 2025. [Herramienta de IA]. Disponible en: https://gemini.google.com. [Accedido: Diciembre 2025].











