MateWiki - Contribuciones del usuario [es]

La Clotoide (grupo 13)

2023-12-14T09:33:08Z

Hugo Gutiérrez Iscar:

{{ TrabajoED | La clotoide. Grupo 13 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Pablo Esteban Coca Hugo Gutiérrez Iscar Nicole Di Natale Berdeal Berta Ramos Domínguez }}

El trabajo correspondiente a nuestro grupo es el número 13, que se centra mayoritariamente en el estudio la clotoide. La clotoide es una curva plana formada por un trozo de espiral, que cumple una serie de condiciones geométricas. Además, desarrollaremos una superficie reglada a partir de una hélice dada. En ambos casos nos enfocaremos en sus estudios matemáticos, así como en su relación con la ingeniería. Para realizar los cálculos con precisión, nos ayudaremos del lenguaje de programación M y de OCTAVE. Los dos programas nos ayudarán a representar gráficamente los elementos pedidos, para así entender los cálculos de manera más visual.
 
 
Consideramos la curva plana dada por la parametrización en coordenadas cartesianas:
 
 
<center> <math> γ(t) = (x(t),y(t)) = (\int_{0}^{t}cos(\frac{s^2}{2})ds, \int_{0}^{t}sin(\frac{s^2}{2})ds), t∈(0,4) </math> </center>
 
 

=='''''La Clotoide'''''==
=== '''Dibujo de la curva''' ===
Comenzaremos el trabajo dibujando la curva dada. Para ello utilizaremos Matlab.
[[File:Dibujo_curva_clotoide1.jpg|410px|miniaturadeimagen|right|'''Dibujo de la curva''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 2000);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);
% Grafica de la clotoide
figure;
plot(x, y);
title('La Clotoide');
% Etiquetado de ejes
xlabel('Eje X');
ylabel('Eje Y');
axis equal;
}}

=== '''Cálculo de vectores Velocidad y Aceleración''' ===

Calcularemos los vectores velocidad y aceleración a partir de la siguiente parametrización:
 
 
<center> <math> γ(t) = (x(t),y(t)) = (\int_{0}^{t}cos(\frac{s^2}{2})ds, \int_{0}^{t}sin(\frac{s^2}{2})ds), t∈(0,4) </math> </center>
 
Para ello usaremos las siguientes fórmulas:
 
*Para el 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>
 
*Para el 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>

[[File:Dibujo_velocidad_aceleracion.jpg|410px|miniaturadeimagen|right|'''Vectores velocidad y aceleración''' ]]

{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 50);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);
% Vectores Velocidad y aceleración
V1 = cos(t.^2/2);
V2 = sin(t.^2/2);
A1 = -t.*sin(t.^2/2);
A2 = t.*cos(t.^2/2);
% Gráfica
figure
hold on
plot (x ,y ,'r') ;
quiver(x,y,V1,V2,"color","r") ; %velocidad
quiver(x,y,A1,A2,"color","g") ; %aceleracion
axis equal
hold off
% Etiqueta de ejes
title('Vectores velocidad y aceleracion');
xlabel("Eje x");
ylabel("Eje y");
}}
 

=== '''Cálculo longitud de la curva''' ===

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

<math> ℓ(γ) = \int_{a}^{b} \left |{\gamma }'(t) \right |= \int_{a}^{b} \left | {\gamma }' (t)\right |= \int_{a}^{b} \sqrt{cos(\frac{t^2}{2})+sin(\frac{t^2}{2})} dt= \int_{0}^{4}1dt= 4-0= 4 </math>
 
 

=== '''Cálculo de los vectores tangente y normal''' ===

Haciendo uso del vector velocidad, calcularemos el vector tangente y normal:
 
 
*El vector tangente:
<math> \vec t (t) =\tfrac{{\gamma}'(t)}{\left |{\gamma}'(t) \right |} = \tfrac{{\gamma}'(t)}{1} = cos(\frac{t^2}{2}) \vec i+sin(\frac{t^2}{2}) \vec j </math>
 
 
*El vector normal:
<math> \vec n (t) =\tfrac{-{y}'(t) \vec i+{x}'(t) \vec j}{\sqrt{{x}'(t)^2 \vec i + {y}'(t)^2 \vec j}} = \tfrac{-sen(\frac{t^2}{2}) \vec i+cos(\frac{t^2}{2}) \vec j}{1} = -sen(\frac{t^2}{2}) \vec i+cos(\frac{t^2}{2}) \vec j </math>
 
 
[[File:Dibujo_tangente_normal.jpg|410px|miniaturadeimagen|right|'''Curva vector, tangente y normal''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 50);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x=arrayfun(x, t);
y=arrayfun(y, t);
%vectores tangente y normal
norma=sqrt(1+4*t.^2);
T1=linspace(0,4,n)./norma;
T2=2*t./norma;
% Grafica
figure;
hold on;
plot(x,y,'r'); %curva
quiver(x,y,T1,T2,'m'); %tangente
quiver(x,y,-T2,T1,'g'); %normal
axis equal
hold off;
title ('Curva, tangente y normal.')
% Etiquetado de ejes
xlabel('Eje X');
ylabel('Eje Y');
axis equal;
}}

=== '''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,4] </math>
[[File:Apartado5curvaturagrupo13.jpg|thumb|miniaturadeimagen|right|'''Dibujo de la curvatura''' ]]
 
 
{{matlab|codigo=
%Representación de la curvatura
t=linspace(0,4,70)
k=t;
figure
plot(t,k,'b');
axis equal
title('Curvatura.');}}

 

=== '''Cálculo de la circunferencia osculatriz''' ===

Dado el punto <math> P=\gamma (1) </math>, es decir <math> t=1 </math>, hallaremos el centro y el radio de la siguiente forma:
 
 
*El radio:
<math> R(t)=\frac{1}{\kappa (t)}=\frac{1}{t} </math>, por lo que el <math> R=1 </math>
 
 

*El centro:
<math> Q(t)=\gamma (t)+\frac{1}{\kappa (t)} \vec n(t) </math>
 
 
<math> Q(t)= (Q_{X},Q_{y})=[\int_{0}^{1}cos(\frac{s^2}{2})ds+\frac{1}{1}\cdot (-sen(\frac{t^2}{2})) , \int_{0}^{1}sin(\frac{s^2}{2})ds+\frac{1}{1}\cdot cos(\frac{t^2}{2})] </math>
 
 
<math> Q(1)=[0+(-sen(\frac{1}{2}))] , [0+cos(\frac{1}{2})]=[-sen(\frac{1}{2}) , cos(\frac{1}{2})] </math>
 

[[File:Apartado6ej11.jpg|410px|miniaturadeimagen|right|'''Circunferencia osculatriz''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 2000);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);

%Calculamos el punto p para t=1
P=[x(1),y(1)];

%Calculamos el vector normal
n=[-sin(1/2),cos(1/2)];

%Curvatura y radio de la curvatura
k=1;
R=1/k;
fprintf('El radio de la curvatura es %d \n',R);

%centro de la circunferencia osculatriz
Q=P+R*n;
fprintf('El centro de la circunferencia es %d \n',Q);

%parametrizamos la circunferencia osculatriz
tt=linspace(0.2*pi,40);
xx=R*cos(tt)+Q(1);
yy=R*sin(tt)+Q(2);

%dibujamos
figure
hold on
plot(x,y,'m','linewidth',3)
plot(P(1),P(2),'*k','linewidth',3) %punto p

plot(xx,yy,'b')
axis equal
hold off
title('Cicunferencia osculatriz.');
}}

=== '''Definiciones e información de interés''' ===

La clotoide, también conocida como espiral de Cornú, es una curva plana, formada con dos extremos de forma parecida a una espiral (aunque no son espirales), enlazadas entre sí por un tramo curvo, formando ambos extremos una simetría central. Su curvatura varía linealmente con la longitud del arco. Debido a esta propiedad la clotoide tiene multitud de aplicaciones.
 
 
La más destacada es el diseño de carreteras y vías ferroviarias. Gracias a sus características, la seguridad del vehículo aumenta notablemente. La curva permite adaptarse al conductor de una forma suave al cambio de trayectoria, por lo que la curva se puede tomar a mayor velocidad.
 
 
Otro campo en el que se aplica la clotoide es en la construcción de montañas rusas. Al construir los bucles verticales con la forma de la curva, la fuerza G se ve reducida, evitando causar daños físicos a los usuarios de la atracción.
 
==== '''Imágenes de estructuras''' ====

La clotoide está presente en múltiples situaciones de nuestra vida cotidiana.
 
 

[[Archivo:clotoide_carretera_1.jpeg|600px|miniaturadeimagen|centro|'''Enlace de carreteras''' ]]
 
[[Archivo:carretera_clotoide_2.jpeg|600px|miniaturadeimagen|centro|'''Nudo vial''' ]]
 
[[Archivo:clotoide_tren_bien.jpg|600px|miniaturadeimagen|centro|'''Vía de tren''' ]]
 
[[Archivo:montaña_rusa.jpeg|600px|miniaturadeimagen|centro|'''Montaña rusa''' ]]

== '''''Superficie reglada. La hélice''''' ==
En esta segunda parte del trabajo nos centraremos en la superficie reglada pedida. Una superficie reglada es aquella superficie generada por una recta de dirección variable que se mueve sobre una curva. En nuestro caso la curva se trata de una hélice.
Consideramos la parametrización de la hélice de <math> \mathbb{R}^3 </math> en coordenadas cartesianas como:
 
 
<center> <math> \gamma (t)=(x_{1}(t),x_{2}(t),x_{3}(t))=(cos(t),sen(t),t), t∈(0,4) </math> </center>
 

=== '''Dibujo de la superficie''' ===
Dibujaremos la superficie reglada asociada a dicha curva mediante segmentos ortogonales de longitud 1 y vector director <math> \vec e_{p} </math>. Para ello seguiremos utilizando Matlab.
 

[[File:Apartado9grupo13.jpg|420px|miniaturadeimagen|right|'''Superficie reglada''' ]]
 
 

{{matlab|codigo=

%Definimos los parámetros
u=linspace(0,1,100);
v=linspace(0,4*pi,100);

%Mallado
[Mu,Mv]=meshgrid(u,v);
Mx=cos(Mv)+(Mu.*cos(Mv));
My=sin(Mv)+(Mu.*sin(Mv));
Mz=Mv;

%Dibujamos
surf(Mx,My,Mz)
shading flat

}}
 
 
 
 
 
Un ejemplo de superficie reglada es:
[[Archivo:parkingclotoide.jpeg|600px|miniaturadeimagen|centro|'''Aparcamiento''' ]]
 
 

=== '''Masa de la superficie reglada''' ===
Suponemos que la densidad de la superficie previamente calculada se rige por la función: <math> f(x_{1},x_{2},x_{3})=10-x_{1}^{2}-x_{2}^{2} </math>
 
 

La masa será igual a: <math> \int_{\phi }f=\int\int f(\phi (u,v))\left |\phi _{u} \times \phi _{v} \right |dvdu </math>
 
 
*La parametrización <math> \phi (u,v) </math> de la superficie es igual a:
 
<math> x=cosv+u\cdot cosv </math>
 
<math> y=sinv+u\cdot sinv </math>
 
<math> z=v </math>
 
 
Para <math> u \epsilon [0,1] </math> y <math> v \epsilon [0,4 \pi] </math>
 
 
*Módulo de los vectores velocidad:
 
Obtenemos los vectores velocidad a partir de la siguiente parametrización:
 
 
<math> \phi _{u}=cosv\vec i+ senv\vec j </math> y <math> \phi _{v}=(-senv-u\cdot senv)\vec i + (cosv+u\cdot cosv)\vec j + \vec k </math>
 
 
Calcularemos el producto vectorial, y posteriormente el módulo:
 
 
<center> <math> \phi _{u} \times \phi _{v}= \begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \\ cosv & senv & 0 \\ -senv-u\cdot senv & cosv+u\cdot cosv & 1 \end{vmatrix}=senv \vec i -cosv \vec j + (1+u) \vec k </math>
 
<math> \left |\phi _{u} \times \phi _{v} \right |= 1+u </math> </center>
 
 
*Densidad en función de la parametrización:
 
<math> f(\phi (u,v))= 10- (cosv+u\cdot cosv)^{2}-(sinv+u\cdot sinv)^{2}=-u^{2}-2u+9 </math>
 
 
Como podemos observar la densidad depende del parámetro <math> v </math>. Previamente, al parametrizar, hemos establecido que <math> v ∈ [0,4\pi] </math>, por lo que v determina la altura de la hélice dada. En consecuencia, podemos asegurar que la densidad se distribuye a lo largo de la altura de nuestra superficie reglada.
 
 
*Cálculo de la masa:
 
<math> masa = \int_{u=0}^{u=1}\int_{v=0}^{v=4 \pi} (9-u^{2}-2u)\cdot (1+u) dvdu= \frac{175}{3}\pi=141.371 </math>

== '''''Bibliografía''''' ==
La Documentación utilizada para la realización de este trabajo es la siguiente:
 
https://www.abc.es/ciencia/abci-clotoide-famosa-curva-hace-segura-carretera-201909210146_noticia.html?ref=https%3A%2F%2Fwww.abc.es%2Fciencia%2Fabci-clotoide-famosa-curva-hace-segura-carretera-201909210146_noticia.html
 
https://trazoide.com/clotoide/
 
https://www.ingenieros-civiles.es/actualidad/actualidad/1/1192/ingenieria-civil-para-dummies-que-son-las-clotoides
 
https://www.microsiervos.com/archivo/ingenieria/loopings-montanas-rusas-circulares-clotoides-fuerzas-g.html

La Clotoide (grupo 13)

2023-12-14T09:21:11Z

Hugo Gutiérrez Iscar:

{{ TrabajoED | La clotoide. Grupo 13 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC23/24|2023-24]] | Pablo Esteban Coca Hugo Gutiérrez Iscar Nicole Di Natale Berdeal Berta Ramos Domínguez }}

El trabajo correspondiente a nuestro grupo es el número 13, que se centra mayoritariamente en el estudio la clotoide. La clotoide es una curva plana formada por un trozo de espiral, que cumple una serie de condiciones geométricas. Además, desarrollaremos una superficie reglada a partir de una hélice dada. En ambos casos nos enfocaremos en sus estudios matemáticos, así como en su relación con la ingeniería. Para realizar los cálculos con precisión, nos ayudaremos del lenguaje de programación M y de OCTAVE. Los dos programas nos ayudarán a representar gráficamente los elementos pedidos, para así entender los cálculos de manera más visual.
 
 
Consideramos la curva plana dada por la parametrización en coordenadas cartesianas:
 
 
<center> <math> γ(t) = (x(t),y(t)) = (\int_{0}^{t}cos(\frac{s^2}{2})ds, \int_{0}^{t}sin(\frac{s^2}{2})ds), t∈(0,4) </math> </center>
 
 

=='''''La Clotoide'''''==
=== '''Dibujo de la curva''' ===
Comenzaremos el trabajo dibujando la curva dada. Para ello utilizaremos Matlab.
[[File:Dibujo_curva_clotoide1.jpg|410px|miniaturadeimagen|right|'''DIBUJO DE LA CURVA''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 2000);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);
% Grafica de la clotoide
figure;
plot(x, y);
title('La Clotoide');
% Etiquetado de ejes
xlabel('Eje X');
ylabel('Eje Y');
axis equal;
}}

=== '''Cálculo de vectores Velocidad y Aceleración''' ===

Calcularemos los vectores velocidad y aceleración a partir de la siguiente parametrización:
 
 
<center> <math> γ(t) = (x(t),y(t)) = (\int_{0}^{t}cos(\frac{s^2}{2})ds, \int_{0}^{t}sin(\frac{s^2}{2})ds), t∈(0,4) </math> </center>
 
Para ello usaremos las siguientes fórmulas:
 
*Para el 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>
 
*Para el 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>

[[File:Dibujo_velocidad_aceleracion.jpg|410px|miniaturadeimagen|right|'''VECTORES VELOCIDAD Y ACELERACIÓN''' ]]

{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 50);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);
% Vectores Velocidad y aceleración
V1 = cos(t.^2/2);
V2 = sin(t.^2/2);
A1 = -t.*sin(t.^2/2);
A2 = t.*cos(t.^2/2);
% Gráfica
figure
hold on
plot (x ,y ,'r') ;
quiver(x,y,V1,V2,"color","r") ; %velocidad
quiver(x,y,A1,A2,"color","g") ; %aceleracion
axis equal
hold off
% Etiqueta de ejes
title('Vectores velocidad y aceleracion');
xlabel("Eje x");
ylabel("Eje y");
}}
 

=== '''Cálculo longitud de la curva''' ===

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

<math> ℓ(γ) = \int_{a}^{b} \left |{\gamma }'(t) \right |= \int_{a}^{b} \left | {\gamma }' (t)\right |= \int_{a}^{b} \sqrt{cos(\frac{t^2}{2})+sin(\frac{t^2}{2})} dt= \int_{0}^{4}1dt= 4-0= 4 </math>
 
 

=== '''Cálculo de los vectores tangente y normal''' ===

Haciendo uso del vector velocidad, calcularemos el vector tangente y normal:
 
 
*El vector tangente:
<math> \vec t (t) =\tfrac{{\gamma}'(t)}{\left |{\gamma}'(t) \right |} = \tfrac{{\gamma}'(t)}{1} = cos(\frac{t^2}{2}) \vec i+sin(\frac{t^2}{2}) \vec j </math>
 
 
*El vector normal:
<math> \vec n (t) =\tfrac{-{y}'(t) \vec i+{x}'(t) \vec j}{\sqrt{{x}'(t)^2 \vec i + {y}'(t)^2 \vec j}} = \tfrac{-sen(\frac{t^2}{2}) \vec i+cos(\frac{t^2}{2}) \vec j}{1} = -sen(\frac{t^2}{2}) \vec i+cos(\frac{t^2}{2}) \vec j </math>
 
 
[[File:Dibujo_tangente_normal.jpg|410px|miniaturadeimagen|right|'''CURVA VECTOR, TANGENTE Y NORMAL''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 50);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x=arrayfun(x, t);
y=arrayfun(y, t);
%vectores tangente y normal
norma=sqrt(1+4*t.^2);
T1=linspace(0,4,n)./norma;
T2=2*t./norma;
% Grafica
figure;
hold on;
plot(x,y,'r'); %curva
quiver(x,y,T1,T2,'m'); %tangente
quiver(x,y,-T2,T1,'g'); %normal
axis equal
hold off;
title ('Curva, tangente y normal.')
% Etiquetado de ejes
xlabel('Eje X');
ylabel('Eje Y');
axis equal;
}}

=== '''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,4] </math>
[[File:Apartado5curvaturagrupo13.jpg|thumb|miniaturadeimagen|right|'''DIBUJO DE LA CURVATURA''' ]]
 
 
{{matlab|codigo=
%Representación de la curvatura
t=linspace(0,4,70)
k=t;
figure
plot(t,k,'b');
axis equal
title('Curvatura.');}}

 

=== '''Cálculo de la circunferencia osculatriz''' ===

Dado el punto <math> P=\gamma (1) </math>, es decir <math> t=1 </math>, hallaremos el centro y el radio de la siguiente forma:
 
 
*El radio:
<math> R(t)=\frac{1}{\kappa (t)}=\frac{1}{t} </math>, por lo que el <math> R=1 </math>
 
 

*El centro:
<math> Q(t)=\gamma (t)+\frac{1}{\kappa (t)} \vec n(t) </math>
 
 
<math> Q(t)= (Q_{X},Q_{y})=[\int_{0}^{1}cos(\frac{s^2}{2})ds+\frac{1}{1}\cdot (-sen(\frac{t^2}{2})) , \int_{0}^{1}sin(\frac{s^2}{2})ds+\frac{1}{1}\cdot cos(\frac{t^2}{2})] </math>
 
 
<math> Q(1)=[0+(-sen(\frac{1}{2}))] , [0+cos(\frac{1}{2})]=[-sen(\frac{1}{2}) , cos(\frac{1}{2})] </math>
 

[[File:Apartado6ej11.jpg|410px|miniaturadeimagen|right|'''CIRCUNFERENCIA OSCULATRIZ''' ]]
{{matlab|codigo=
% Definimos los Parámetros
t = linspace(0, 4, 2000);
% Definimos la funcion
x = @(t) integral(@(s) cos(s.^2/2), 0, t);
y = @(t) integral(@(s) sin(s.^2/2), 0, t);
% Calcular las coordenadas de la clotoide
x = arrayfun(x, t);
y = arrayfun(y, t);

%Calculamos el punto p para t=1
P=[x(1),y(1)];

%Calculamos el vector normal
n=[-sin(1/2),cos(1/2)];

%Curvatura y radio de la curvatura
k=1;
R=1/k;
fprintf('El radio de la curvatura es %d \n',R);

%centro de la circunferencia osculatriz
Q=P+R*n;
fprintf('El centro de la circunferencia es %d \n',Q);

%parametrizamos la circunferencia osculatriz
tt=linspace(0.2*pi,40);
xx=R*cos(tt)+Q(1);
yy=R*sin(tt)+Q(2);

%dibujamos
figure
hold on
plot(x,y,'m','linewidth',3)
plot(P(1),P(2),'*k','linewidth',3) %punto p

plot(xx,yy,'b')
axis equal
hold off
title('Cicunferencia osculatriz.');
}}

=== '''Definiciones e información de interés''' ===

La clotoide, también conocida como espiral de Cornú, es una curva plana, formada con dos extremos de forma parecida a una espiral (aunque no son espirales), enlazadas entre sí por un tramo curvo, formando ambos extremos una simetría central. Su curvatura varía linealmente con la longitud del arco. Debido a esta propiedad la clotoide tiene multitud de aplicaciones.
 
 
La más destacada es el diseño de carreteras y vías ferroviarias. Gracias a sus características, la seguridad del vehículo aumenta notablemente. La curva permite adaptarse al conductor de una forma suave al cambio de trayectoria, por lo que la curva se puede tomar a mayor velocidad.
 
 
Otro campo en el que se aplica la clotoide es en la construcción de montañas rusas. Al construir los bucles verticales con la forma de la curva, la fuerza G se ve reducida, evitando causar daños físicos a los usuarios de la atracción.
 
==== '''Imágenes de estructuras''' ====

La clotoide está presente en múltiples situaciones de nuestra vida cotidiana.
 
 

[[Archivo:clotoide_carretera_1.jpeg|600px|centro]]
 
[[Archivo:carretera_clotoide_2.jpeg|600px|centro]]
 
[[Archivo:clotoide_tren_bien.jpg|600px|centro]]
 
[[Archivo:montaña_rusa.jpeg|600px|centro]]

== '''''Superficie reglada. La hélice''''' ==
En esta segunda parte del trabajo nos centraremos en la superficie reglada pedida. Una superficie reglada es aquella superficie generada por una recta de dirección variable que se mueve sobre una curva. En nuestro caso la curva se trata de una hélice.
Consideramos la parametrización de la hélice de <math> \mathbb{R}^3 </math> en coordenadas cartesianas como:
 
 
<center> <math> \gamma (t)=(x_{1}(t),x_{2}(t),x_{3}(t))=(cos(t),sen(t),t), t∈(0,4) </math> </center>
 

=== '''Dibujo de la superficie''' ===
Dibujaremos la superficie reglada asociada a dicha curva mediante segmentos ortogonales de longitud 1 y vector director <math> \vec e_{p} </math>. Para ello seguiremos utilizando Matlab.
 

[[File:Apartado9grupo13.jpg|420px|miniaturadeimagen|right|'''SUPERFICIE REGLADA''' ]]
 
 

{{matlab|codigo=

%Definimos los parámetros
u=linspace(0,1,100);
v=linspace(0,4*pi,100);

%Mallado
[Mu,Mv]=meshgrid(u,v);
Mx=cos(Mv)+(Mu.*cos(Mv));
My=sin(Mv)+(Mu.*sin(Mv));
Mz=Mv;

%Dibujamos
surf(Mx,My,Mz)
shading flat

}}
 
 
 
 
 
Un ejemplo de superficie reglada es:
[[Archivo:parkingclotoide.jpeg|600px|centro]]
 
 

=== '''Masa de la superficie reglada''' ===
Suponemos que la densidad de la superficie previamente calculada se rige por la función: <math> f(x_{1},x_{2},x_{3})=10-x_{1}^{2}-x_{2}^{2} </math>
 
 

La masa será igual a: <math> \int_{\phi }f=\int\int f(\phi (u,v))\left |\phi _{u} \times \phi _{v} \right |dvdu </math>
 
 
*La parametrización <math> \phi (u,v) </math> de la superficie es igual a:
 
<math> x=cosv+u\cdot cosv </math>
 
<math> y=sinv+u\cdot sinv </math>
 
<math> z=v </math>
 
 
Para <math> u \epsilon [0,1] </math> y <math> v \epsilon [0,4 \pi] </math>
 
 
*Módulo de los vectores velocidad:
 
Obtenemos los vectores velocidad a partir de la siguiente parametrización:
 
 
<math> \phi _{u}=cosv\vec i+ senv\vec j </math> y <math> \phi _{v}=(-senv-u\cdot senv)\vec i + (cosv+u\cdot cosv)\vec j + \vec k </math>
 
 
Calcularemos el producto vectorial, y posteriormente el módulo:
 
 
<center> <math> \phi _{u} \times \phi _{v}= \begin{vmatrix} \vec{i} & \vec{j} & \vec{k} \\ cosv & senv & 0 \\ -senv-u\cdot senv & cosv+u\cdot cosv & 1 \end{vmatrix}=senv \vec i -cosv \vec j + (1+u) \vec k </math>
 
<math> \left |\phi _{u} \times \phi _{v} \right |= 1+u </math> </center>
 
 
*Densidad en función de la parametrización:
 
<math> f(\phi (u,v))= 10- (cosv+u\cdot cosv)^{2}-(sinv+u\cdot sinv)^{2}=-u^{2}-2u+9 </math>
 
 
Como podemos observar la densidad depende del parámetro <math> v </math>. Previamente, al parametrizar, hemos establecido que <math> v ∈ [0,4\pi] </math>, por lo que v determina la altura de la hélice dada. En consecuencia, podemos asegurar que la densidad se distribuye a lo largo de la altura de nuestra superficie reglada.
 
 
*Cálculo de la masa:
 
<math> masa = \int_{u=0}^{u=1}\int_{v=0}^{v=4 \pi} (9-u^{2}-2u)\cdot (1+u) dvdu= \frac{175}{3}\pi=141.371 </math>

== '''''Bibliografía''''' ==
La Documentación utilizada para la realización de este trabajo es la siguiente:
 
https://www.abc.es/ciencia/abci-clotoide-famosa-curva-hace-segura-carretera-201909210146_noticia.html?ref=https%3A%2F%2Fwww.abc.es%2Fciencia%2Fabci-clotoide-famosa-curva-hace-segura-carretera-201909210146_noticia.html
 
https://trazoide.com/clotoide/
 
https://www.ingenieros-civiles.es/actualidad/actualidad/1/1192/ingenieria-civil-para-dummies-que-son-las-clotoides
 
https://www.microsiervos.com/archivo/ingenieria/loopings-montanas-rusas-circulares-clotoides-fuerzas-g.html