La Clotoide (Grupo 21)

De MateWiki
Saltar a: navegación, buscar
Trabajo realizado por estudiantes
Título La clotoide. Grupo 21
Asignatura Teoría de Campos
Curso 2025-26
Autores Paula Rodríguez Rey, Ignacio Moya Casasola, Adrián Eguilleor Prieto, Mencía Benitez Del Castillo Mateos, Pablo Cortina Gómez.
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción.

Matemáticamente, una clotoide es una curva que parte siendo tangente al eje de abscisas y cuya curvatura aumenta progresivamente, de modo que su radio de curvatura disminuye en proporción inversa a la longitud recorrida sobre la propia curva.
Para estudiar sus características, examinaremos primero los vectores de velocidad y aceleración, junto con los elementos del triedro de Frenet.
Más adelante, relacionaremos estos conceptos con su utilización en ingeniería civil.

2 Dibujo de la curva.

La expresión matemática de la clotoide es:

[math] \gamma (t)=(x(t),y(t))=\left ( \int_{0}^{t}cos(\frac{s^2}{2})ds, \int_{0}^{t}sin(\frac{s^2}{2})ds \right ), t\in (0,5) [/math]


La representación gráfica de la curva se ha obtenido mediante el siguiente código:

Figura 1: Clotoide
clear; clc; clf;
% Definimos los parámetros
 L = 4;       
 n = 500;  
 t = linspace(0, L, n);  

% Definimos los vectores para las coordenadas x y y
 x = zeros(1, n);
 y = zeros(1, n);

% Definimos las funciones
 f1= @(s) cos(s.^2/2);
 f2= @(s) sin(s.^2/2);

% Aproximamos la integral usando el método del rectángulo
for i = 2:n
    % Para x(t), sumamos la función cos(s^2 / 2) de t = 0 hasta t = t(i)
    x(i) = x(i-1) + f1(t(i-1)) * (t(i) - t(i-1));
    
    % Para y(t), repetimos el método usando sin(s^2 / 2) 
    y(i) = y(i-1) + f2(t(i-1))* (t(i) - t(i-1));
end

% Representamos gráficamente la curva
figure;
plot(x, y);
axis equal;
xlabel('eje x');
ylabel('eje y');
title('Curva de la clotoide');
grid on;

3 Velocidad y aceleración.

Para calcular ambos vectores, se han aplicado las siguientes fórmulas de velocidad [math] \dot{\gamma } [/math] y aceleración [math] \ddot{\gamma } [/math]

[math] \vec{{\gamma }'}=cos(\frac{t^2}{2}) \vec{i} +sin(\frac{t^2}{2})\vec{j} [/math]
[math] \vec{{\gamma }''}= -t\cdot sin(\frac{t^2}{2}) \vec{i} +t\cdot cos(\frac{t^2}{2})\vec{j} [/math]


Para representarlo, partiremos del código y gráfica del apartado anterior, añadiendo el siguiente código, y obteniendo:

Figura 2: Vectores velocidad y aceleración junto a la clotoide
% Calculamos las derivadas numéricas de x(t) y y(t) (velocidad)
dx = cos(t.^2/2);  % Derivada primera de x(t)
dy = sin(t.^2/2);  % Derivada primera de y(t)

% Calculamos las derivadas de las velocidades (aceleración)
ddx = -t.*sin(t.^2/2);  % Derivada segunda de x(t)
ddy = t.*cos(t.^2/2);  % Derivada segunda de y(t)

hold on;

% Dibujamos los vectores de velocidad (negro) y aceleración (azul)
for i = 1:4:n  
    % Vectores de velocidad
    quiver(x(i), y(i), dx(i), dy(i), 0.2, 'k', 'LineWidth', 0.5, 'MaxHeadSize',1);
    
    % Vectores de aceleración
    quiver(x(i), y(i), ddx(i), ddy(i), 0.025, 'b', 'LineWidth', 0.5, 'MaxHeadSize',0.5);
end

% Etiquetas y configuración de la gráfica
title('Curva, Vectores de Velocidad y Aceleración');
legend('Curva', 'Velocidad', 'Aceleración','Location','Best');
hold off;


4 Longitud de la curva


La longitud de la curva viene dada por la siguiente expresión:

[math] L(γ'(t))=\int_{0}^{t}|γ'(t)|dt [/math]


Como se ha plasmado en el apartado anterior:

[math] \vec{{\gamma }'}= cos(\frac{t^2}{2}) \vec{i} +sin(\frac{t^2}{2})\vec{j} [/math]


Cuyo módulo es:

[math] |γ′(t)| = \sqrt {cos^2(\frac{t^2}{2})+sin^2(\frac{t^2}{2})} = \sqrt {1} = 1 [/math]


Por tanto la longitud es:

[math] L(γ) = \int_{0}^{4}\sqrt {cos^2(\frac{t^2}{2})+sin^2(\frac{t^2}{2})}dt = \int_{0}^{4}1dt = 4-0 = 4 [/math]


5 Vectores tangente y normal

Los vectores tangente y normal de la clotoide vienen dadas por:

[math]\vec{t}(t)=\frac{\gamma {}'(t)}{\left | \gamma {}'(t) \right |}=\frac{cos(\frac{t^2}{2})\vec{i}+sin(\frac{t^2}{2})\vec{j}}{1} [/math]


[math]\vec{n}(t)={\frac{\gamma'(t) \times \gamma''(t)}{|\gamma'(t) \times \gamma''(t)|}}\times{\frac{cos(\frac{t^2}{2})\vec{i}+sin(\frac{t^2}{2})\vec{j}}{1}}= {-sin(\frac{t^2}{2})\vec{i} + cos(\frac{t^2}{2})\vec{j}} [/math]


Para representarlo, partiremos del código y gráfica del apartado anterior, añadiendo el siguiente código, y obteniendo:

Figura 3: Vectores tangente y normal junto a la clotoide
% Calculamos los vectores tangente x(t) e y(t) 
tx = cos(t.^2/2); 
ty = sin(t.^2/2);  

% Calculamos los vectores normal x(t) e y(t)
nx = -sin(t.^2/2); 
ny = cos(t.^2/2);  

hold on;

% Dibujamos los vector tangente (negro) y normal (azul)
for i = 1:4:n  
    % Vector tangente
    quiver(x(i), y(i), tx(i), ty(i), 0.2, 'k', 'LineWidth', 0.5, 'MaxHeadSize',1);
    
    % Vector normal
    quiver(x(i), y(i), nx(i), ny(i), 0.1, 'b', 'LineWidth', 0.5, 'MaxHeadSize',0.5);
end

% Etiquetas y configuración de la gráfica
title('Curva, Vectores tangente y normal');
legend('Curva', 'Tangente', 'Normal','Location','Best');
hold off;

6 Curvatura k(t).

La curvatura se calcula con la siguiente fórmula:

[math] k(t)=\frac{|\gamma'(t) \times \gamma''(t)|}{|\gamma'(t)|^3}=\frac{\left| \left( \cos\left(\frac{t^2}{2}\right) \vec{i} + \sin\left(\frac{t^2}{2}\right) \vec{j} \right) \times \left( -t \sin\left(\frac{t^2}{2}\right) \vec{i} + t \cos\left(\frac{t^2}{2}\right) \vec{j} \right) \right|}{\left| \cos\left(\frac{t^2}{2}\right) \vec{i} + \sin\left(\frac{t^2}{2}\right) \vec{j} \right|^3} = t [/math]


La gráfica de la curvatura se calcula mediante el siguiente código de Matlab

Figura 4: Curvatura
% Definimos el parámetro t
 t=linspace(0,4,50);
% Definimos la curvatura k(t)
 k=t;
% Representamos la gráfica de la curvatura
 figure;
 plot(k,t);
 title('Curvatura');
 xlabel('Eje x');
 ylabel('Eje y');










7 Circunferencia osculatriz.

La circunferencia osculatriz es una aproximación local de la curva en cada punto de esta, es decir, la circunferencia tiene la misma tangente, curvatura y centro de curvatura que la curva en cada punto.
Dada esta definición y dado P= [math] \gamma (1.5) [/math], es decir, t=1.5, el radio de la circunferencia osculatriz y su centro son:

[math]R(t)=\frac{1}{\kappa(t)}[/math]

[math]Q(t)=\gamma (t)+\frac{1}{\kappa (t)}\bar{n}(t)[/math]



Realizando las operaciones correspondientes, tenemos:

[math]R(1.5)=\frac{1}{1.5}[/math]

[math]Q(1.5) = \left\{\begin{matrix} Q_x(1.5)=\int_{0}^{1.5}cos(\frac{s^2}{2})ds - (\frac{sin(1.5)}{2})\\\ Q_y(1.5)=\int_{0}^{1.5}sin(\frac{s^2}{2})ds + (\frac{cos(1.5)}{2}) \end{matrix}\right.[/math]


Con el centro recientemente calculado, se realiza el gráfico, añadiendo el siguiente código, al anterior de la clotoide, y por tanto, obteniendo la circunferencia osculatriz:

Figura 5: Circunferencia osculatriz y la curva

{{matlab|codigo= % Calculamos las integrales de la curva para t = 1.5 X1 = integral(f1,0,1.5); Y1 = integral(f2,0,1.5);

% Radio de la circunferencia osculatriz R = 1/1.5;

% Vector normal unitario en t = 1.5 nx = -sin(1.5^2/2); ny = cos(1.5^2/2);

% Centro de la circunferencia osculatriz Qx = X1 + R*nx; Qy = Y1 + R*ny;

% Parametrización de la circunferencia theta = linspace(0,2*pi,500); Cx = Qx + R*cos(theta); Cy = Qy + R*sin(theta);

% Representación hold on plot(x,y,'r')  % Clotoide plot(Cx,Cy,'b')  % Circunferencia osculatriz axis equal xlabel('Eje x') ylabel('Eje y') title('Circunferencia osculatriz en t = 1.5') hold off

8 Propiedades para la ingeniería.

La clotoide describe un fenómeno de transición suave entre una trayectoria recta y una curva circular, ya que, como se ha expuesto anteriormente, su curvatura crece de forma lineal. Conociendo esto, en el punto de inicio, el radio de curvatura es infinito, y a medida que avanza, el radio disminuye hasta tomar un valor finito, estableciendo una curvatura más definida.

Su principal aplicación en la ingeniería es el diseño de carreteras y ferrocarriles en el que la clotoide se usa para suavizar la transición entre un tramo recto y una curva circular. Esta transición es crucial, ya que evita cambios abruptos en la aceleración centrípeta y la ajusta gradualmente. Sin una transición suave, se podría generar incomodidad o incluso peligro para los vehículos y pasajeros ya que se enfrentarían a un aumento agresivo de las fuerzas centrípetas, lo que puede afectar a la estabilidad.

Las propiedades de la clotoide ofrecen otras aplicaciones como ayudar a mantener un flujo de agua estable, diseñar rutas de entrada y salida para embarcaciones en los puertos e incluso para construir montañas rusas.