Revisión actual del 14:13 12 dic 2024

En este trabajo vamos a exponer la curva conocida como clotoide y sus numerosas propiedades en el ámbito civil. Una clotoide es una curva cuya característica principal es que la tasa de cambio de la curvatura es constante a lo largo de su longitud, es decir, aumenta o disminuye de manera progresiva y suave, sin cambios bruscos.

Consideramos la curva plana dada por la parametrización en coordenadas cartesianas:

1 La Clotoide

1.1 Dibujo de la curva

Comenzaremos el trabajo dibujando la curva dada. Para ello utilizaremos matlab. (L=5)

Dibujo de la curva

t = linspace(0, 5, 200);
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
x = arrayfun(x, t);
y = arrayfun(y, t);
figure;
plot(x, y);
title('La Clotoide');
xlabel('X');
ylabel('Y');
axis equal;

1.2 Cálculo de vectores velocidad y aceleración

Calcularemos los vectores velocidad y aceleración a partir de la siguiente parametrización:

Calculo vector velocidad: [math] {\gamma }'(t)={x}'(t) \vec i + {y}'(t) \vec j \rightarrow {\gamma }'(t)=cos(\frac{t^2}{2}) \vec i+sin(\frac{t^2}{2}) \vec j [/math]
Calculo vector aceleración: [math] {\gamma }''(t)={x}''(t) \vec i + {y}''(t) \vec j \rightarrow {\gamma }'(t)=-t\cdot sin(\frac{t^2}{2}) \vec i+t\cdot cos(\frac{t^2}{2}) \vec j [/math]

Vectores velocidad y aceleración

t = linspace(0, 5, 150);
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
x = arrayfun(x, t);
y = arrayfun(y, t); 
V1 = cos(t.^2/2);
V2 = sin(t.^2/2);
A1 = -t.*sin(t.^2/2);
A2 = t.*cos(t.^2/2);
figure
hold on
plot (x ,y ,'b') ; 
quiver(x,y,V1,V2,"color","g") ; 
quiver(x,y,A1,A2,"color","r") ; 
axis equal
hold off
title('Vectores velocidad y aceleracion');
xlabel("X");
ylabel("Y");

1.3 Cálculo longitud de la curva

Utilizando la siguiente fórmula calcularemos la longitud de la curva:

1.4 Cálculo de los vectores tangente y normal

Haciendo uso del vector velocidad, calcularemos el vector tangente y normal:

El vector tangente:

El vector normal:

Curva vector tangente y normal

t = linspace(0, 5, 100);
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
x=arrayfun(x, t);
y=arrayfun(y, t);
norma=1;
T1 = cos(t.^2/2)./norma;
T2 = sin(t.^2/2)./norma;
N1= -sin((t.^2)./2);
N2= cos ((t.^2)./2);
figure;
hold on;
plot(x,y,'b'); %curva
quiver(x,y,T1,T2,"color",'r');
quiver(x,y,N1,N2,"color",'g');
axis equal
hold off;
title ('Curva, tangente y normal.')
xlabel('Eje X');
ylabel('Eje Y');

1.5 Cálculo de la curvatura

Estudiaremos la curvatura en el punto [math] γ(t) [/math] que viene dada por la siguiente fórmula:

[math] \kappa (t)=\frac{{x}'(t)\cdot {y}''(t)-{x}''(t)\cdot {y}'(t)}{({x}'(t)^2 + {y}'(t)^2)^{\frac{3}{2}}} [/math] [math] =\frac{cos(\frac{t^2}{2})\cdot t\cdot cos(\frac{t^2}{2})-(-t\cdot sin(\frac{t^2}{2})\cdot sin(\frac{t^2}{2}))}{\sqrt{[(cos(\frac{t^2}{2}))^2+(sin(\frac{t^2}{2}))^2]^3}}=\frac{t}{1}=t , t ∈ [0,5] [/math]

Dibujo de la curvatura

t=linspace(0,5,70)
k=t;
figure
plot(t,k,'b');
axis equal
title('Curvatura.');
xlabel('Eje X');
ylabel('Eje Y');

1.6 Cálculo de la circunferencia osculatriz

Dado el punto [math] P=\gamma (2) [/math], es decir [math] t=2 [/math], hallaremos el centro y el radio de la siguiente forma:

El radio:

El centro:

Parametrizamos el centro de la circunferencia osculatriz con la siguiente fórmula:

derecha

% Parámetros iniciales
t = linspace(0, 5, 2000);
x = @(t) arrayfun(@(t_val) integral(@(s) cos(s.^2 / 2), 0, t_val), t);
y = @(t) arrayfun(@(t_val) integral(@(s) sin(s.^2 / 2), 0, t_val), t);
% Cálculo de la clotoide
xc = x(t);
yc = y(t);
% Selección de un punto t1 para la circunferencia osculatriz
t1 = 1;
x1 = integral(@(s) cos(s.^2 / 2), 0, t1);
y1 = integral(@(s) sin(s.^2 / 2), 0, t1);
P = [x1, y1];
fprintf('El punto de curvatura es (%f, %f)\n', P(1), P(2));
% Cálculo de la normal y radio de curvatura
n = [-sin(t1^2 / 2), cos(t1^2 / 2)]; % Vector normal unitario
k = 1;
R = 1 / k; % Radio de curvatura
fprintf('El radio de curvatura es %f\n', R);
% Centro de la circunferencia
Q = P + R * n;
fprintf('El centro de la circunferencia es (%f, %f)\n', Q(1), Q(2));

% Coordenadas de la circunferencia
tt = linspace(0, 2 * pi, 40);
xx = R * cos(tt) + Q(1);
yy = R * sin(tt) + Q(2);
% Gráfica
figure;
hold on;
plot(xc, yc, 'm', 'LineWidth', 1); % Clotoide
plot(x1, y1, '*k', 'LineWidth', 1); % Punto de curvatura
plot(xx, yy, 'b'); % Circunferencia osculatriz
hold off;
title('Circunferencia osculatriz');
xlabel('Eje X');
ylabel('Eje Y');
axis equal;

1.7 Definición e información relevante sobre la clotoide

La clotoide es una curva matemática, cuya principal característica es que su curvatura cambia de manera gradual. En otras palabras, una clotoide es una curva cuya tasa de cambio de curvatura es constante, lo que significa que la transición entre una curva recta y una curva circular que ocurre de manera progresiva y no abrupta.

Esta curva es usada en el ambiente de la ingeniería destacando carreteras, ferrocarriles, y diseño de pistas de aeropuertos. El concepto de clotoide fue desarrollado en el siglo XVIII por el matemático alemán Leonhard Euler quien la uso en caminos y ferrocarriles

1.8 Fotos de algunas estructuras en el ámbito civil

carretera

carretera

2 Helicoide cónico

2.1 Dibujo de la superficie

En este apartado se va a dibujar la superficie reglada asociada a la parametrización de una curva mediante segmentos ortogonales de longitud 1 y vector director eρ. Esta superficie se conoce habitualmente como helicoide cónico

• La parametrización de la superficie en cartesianas es la siguiente:

En función de u y v:

Para que la gráfica muestre la superficie reglada hay que extender los segmentos desde la hélice para cada punto de la curva base, la superficie se expresa como:

Porque el vector director eρ en cartesianas es [math] cos(\theta) \vec i + sin(\theta) \vec j [/math] y ρ controla la magnitud de los segmentos en la dirección de eρ.

Dibujo de la superficie

u=linspace(2*pi, 6*pi, 100); % Valores de u
v=linspace(0, 1, 100);     % Valores de v
[U, V] = meshgrid(u, v);   % Malla para parametrización
% Coordenadas de la superficie reglada
X=U.*cos(U)+V.*cos(U);
Y=U.*sin(U)+V.*sin(U);
Z=U;
% Gráfica de la superficie 
surf(X, Y, Z, 'EdgeColor', 'none');
% Configuración de la gráfica
colormap('winter');
c=colorbar; 
c.Label.String='Valores en Z'; 
axis equal;
xlabel('X');
ylabel('Y'); 
zlabel('Z');
title('Superficie Reglada de la Hélice Cónica'); 
grid on;

2.2 Aplicaciones en ingeniería

Se aplica esta geometría en la fabricación de engranajes cónicos helicoidales, que se emplean en diferentes vehículos de motor, tienen una mayor superficie de contacto que los engranajes rectos. También hay aplicaciones en la construcción:

Torre de observación en el parque Killesberg, Stuttgart.

3 Masa de la superficie reglada.

Suponiendo que la densidad de la superficie obtenida en el apartado anterior está descrita por la función:

Una vez conocida la superficie y la densidad podemos calcular su masa utilizando la siguiente expresión:

Gracias a los cálculos realizados anteriormente podemos calcular la masa :

El vector de posición esta dado por:

Así que deberemos aplicar la siguiente fórmula para obtener rᵤ y rᵥ:

Podemos calcular el producto vectorial de rᵤ y rᵥ como:

De modo que el módulo de dicho producto vendrá dado por:

La expresión en función del vector [math]\vec{{r}}(u,v)[/math] queda:

Por último, aproximamos el resultado de la masa mediante el método del rectángulo, haciendo uso de Matlab.

% Pedimos el número de rectángulos al usuario
n = input('Número de rectángulos: ');

% Definimos los extremos de los intervalos
a = 0; % Límite inferior en u
b = 1; % Límite superior en u
c = 2 * pi; % Límite inferior en v
d = 6 * pi; % Límite superior en v

% Calculamos los pasos en cada eje
h_u = (b - a) / n;
h_v = (d - c) / n;

% Creamos las particiones para u y v
u = a:h_u:b;
v = c:h_v:d;

% Inicializamos el volumen total
vol = 0;

% Iteramos sobre cada rectángulo
for i = 1:n
    for j = 1:n
        % Calculamos el valor de u y v en el centro del rectángulo
        u_c = u(i) + h_u / 2;
        v_c = v(j) + h_v / 2;

        % Calculamos la función en el centro del rectángulo
        f = (-u_c^2 - 2 * u_c + 99) * sqrt(1 + (1 + u_c^2));

        % Calculamos el área del rectángulo proyectado
        area_rect = h_u * h_v;

        % Calculamos el volumen del rectángulo
        vol_rect = f * area_rect;

        % Acumulamos el volumen total
        vol = vol + vol_rect;
    end
end

% Mostramos el resultado
fprintf('Para %d rectángulos, el volumen bajo la curva es: %.3f\n', n, vol);