Tesis Quadrotor PDF
Tesis Quadrotor PDF
Tesis Quadrotor PDF
Director
CARLOS BORRAS PINILLA
Ingeniero Mecánico
_________________________________
_________________________________
_________________________________
_________________________________
_________________________________
_________________________________
_________________________________
Presidente de Jurado
_________________________________
Jurado
_________________________________
Jurado
A mi madre Leonor, por haberme dado la vida, por todo su amor y esfuerzo para
hacerme un hombre integro.
A mis familiares, Roberto, Isaías, Soraya, Alexander y Juan Sebastián, por sus
buenos deseos de ver culminada esta etapa de mi vida.
A Yuleidy Álvarez, por estar a mi lado y acompañarme cada día con su amor y
comprensión.
Primero a Dios por darme el alivio y fortaleza en los momentos difíciles de mi vida
A mi madre Mabel Vega por creer siempre en mí, por su apoyo incondicional y
grandes consejos que hacen que cada día sea un mejor ser humano,
A mi abuela Delfa que con su cariño desmedido hace que supere todo obstáculo
A mis amigos que han sido pocos pero suficientes Dairo, Alex, Luza, Yuli, Carolina,
Mauricio, Cristian, Alvaro, porque siempre he recibido las mejores palabras de ánimo
cuando lo he necesitado.
A mis amigos Freddy Augusto Bonilla, Álvaro Romero, Pedro David Medina,
Martín Mojica, Erick Humberto Hernández, Álvaro Celis, Juan Fernando
Arango, Manuel Niño, Nelson Acosta, por su invaluable ayuda en los
momentos difíciles de mi vida.
A mis amigos Álvaro Tasco, Cristian Pérez, Alexander Ríos, Dairo Álvarez,
por su ayuda y colaboración en el desarrollo de este proyecto.
2. MARCO TEORICO 23
de Vuelo 30
3. METODOLOGIA DE DISEÑO 40
9. CONCLUSIONES 254
BIBLIOGRAFIA 259
ANEXOS 263
LISTA DE FIGURAS
Pág.
Figura1. Prototipo de Quadrotor Proto 1 (Virginia Tech) 5
Figura 83. Alarma sonora de baja tensión para baterías LiPo de 3-4S 133
Figura 133. Empleo de pistola de aire caliente sobre la fibra de carbono 209
Figura 134. Sección superior e inferior del fuselaje en fibra de carbono 210
Figura 135. Sección superior e inferior del brazo en fibra de carbono 210
Figura 173. Respuesta del control Kp= 4, Ki=0.020, Kd=10 Prueba 5 248
Figura 174. Respuesta del control Kp= 4, Ki=0.020, Kd=20, prueba 6 249
Pág.
TITULO:
DISEÑO Y CONSTRUCCIÓN DEL PROTOTIPO EXPERIMENTAL DE LA ESTRUCTURA
PARA UN VEHICULO AEREO NO TRIPULADO (UAV) TIPO QUADROTOR*
AUTORES:
JOSE ALBERTO HERRERA VASQUEZ
RICARDO ALEXIS RINCON VEGA**
PALABRAS CLAVE:
Vehículo, Aéreo, Tripulado, Quadrotor, Carbono, Micro-controlador, Sensores, Control
DESCRIPCION:
El objetivo de este proyecto de grado es el de proveer a la Escuela de Ingeniería Mecánica de
la Universidad Industrial de Santander un vehículo aéreo no tripulado experimental tipo
Quadrotor que sirva como herramienta de desarrollo tecnológico en el área de nuevos
materiales y sus procesos de manufactura, así como en la incorporación de sistemas de
control autónomo.
En este documento se plasman las fases de diseño, los procesos de selección de materiales y
manufactura requeridos para la fabricación de los componentes estructurales del vehículo,
donde se hizo empleo de herramientas computacionales de diseño mecánico y maquinado
como SolidWorks, Ansys y PowerMILL.
Así mismo, se desarrolla el análisis dinámico del vehículo a través del modelado matemático
de sus principales movimientos, realizando la simulación respectiva a través de herramientas
de análisis dinámico como Simulink de Matlab para determinar su comportamiento
optimizándolo a través de métodos de control como filtro PID para regular el alto grado de
inestabilidad que estos tipos de vehículos poseen.
________________________
* Trabajo de Grado
** Facultad de Ingenierías Físico-Mecánicas. Escuela de Ingeniería Mecánica. Director: Carlos
Borras
SUMMARY
AUTHORS:
JOSE ALBERTO HERRERA VASQUEZ
RICARDO ALEXIS RINCON VEGA**
KEYWORDS:
Vehicle, aerial, Manned, Quadrotor, Carbon, Microcontroller, Sensors, Control.
DESCRIPTION:
The objective of this degree project is to provide the Mechanical Engineering faculty of The
Santander Industrial University with one experimental Quadrotor type unmaneed aerial
vehicle that serves as a tool for technological development in the new materials and
manufacturing processes field, as well as the incorporation of autonomous control systems.
This document embodies the design phase, selection processes and manufacturing materials
required for the manufacture of the vehicle structural components, where computational
tools for mechanical design and machining as SolidWorks, Ansys and PowerMILL were
used.
In the same way, a vehicle dynamic analysis is developed through dynamic analysis tools,
such as Simulink of Matlab, to determine its behavior and optimizing it through control
methods, for example, a PID filter to regulate the high instability grade possessed by this
kind of vehicles.
The final result is the Quadrotor experimental prototype, fully functional, equipped with
stability and control systems based on micro-controllers, sensors and radio-control that meets
the proposed objectives of this project, which can be employed in future researches.
_______________
* Work Degree
** Faculty of Physico-mechanical Engineerings. School of Mechanical Engineering. Project
director: Carlos Borrás Pinilla.
INTRODUCCION
Para lograr una óptima versatilidad en las formas de vuelo que permite una
plataforma aérea de despegue y aterrizaje vertical como la aquí diseñada, se
requiere la incorporación de sistemas de control basados en arquitecturas
dinámicas, de gran capacidad de procesamiento y provistos con lenguajes de
fácil programación, en entornos amigables al usuario.
2
Esta tecnología no permite ningún tipo de exposición o riesgo para el
operador, constituyéndose así, en uno de los campos de la ingeniería con más
amplio desarrollo en algunos países alrededor del mundo en los últimos
años.
Ref.[1] S.D, Prior, et al. THE FUTURE OF BATTLEFIELD MICRO AIR VEHICLE SYSTEMS,
Middlesex University, School of Engineering and Information Sciences, Department of Product Design
and Engineering, London, United Kingdom, 2009.
3
Como se puede observar, los precios de estos sistemas aéreos oscilan entre
£5000 y £30000 ($15’000.000.oo y $90’000.000.oo pesos colombianos),
dependiendo del tamaño y equipo transportado para navegación y
observación. Su alto valor comercial y servicio postventa hace que estos
dispositivos sean de difícil adquisición para algunos operadores que
disponen de bajo presupuesto, pero requieren contar con esta tecnología.
4
Figura1. Prototipo de Quadrotor Proto 1 (Virginia Tech)
Fuente: http://www.dept.aoe.vt.edu/~mason/Mason_f/IARC-FinalReport-v6.0.pdf
5
Figura 2. Prototipo de Quadrotor (Columbia University)
Fuente: http://me.columbia.edu/seniordesigns/2007/QUAVe/FinalReport.pdf
6
1.2. JUSTIFICACIÓN PARA SOLUCIONAR EL PROBLEMA
7
1.3. OBJETIVOS DEL TRABAJO DE GRADO
8
Las especificaciones técnicas del vehículo aéreo experimental
propuesto son:
- Dimensiones Aproximadas:
Longitud 800 mm (Entre ejes de rotores)
Altura 150 mm
- Pesos:
Peso Estructura 1000 -1500 g
Peso de Carga Útil 500 g
- Autonomía de Vuelo: 15 minutos
9
- Control de Velocidad del Motor:
Circuito Electrónico Voltaje de Entrada 20V
Corriente de Salida 35 A
10
1.4. JUSTIFICACIÓN DE LA SOLUCIÓN
11
1.5. REQUERIMIENTOS DEL CONSUMIDOR
12
Figura 3. Balance de fuerzas en vuelo vertical
13
Figura 4. Daños en la estructura de un Quadrotor
Fuente: http://coe.uncc.edu/orgs/robotics/quadrotor/
14
Figura 5. Estructura de Quadrotor tipo abierta
Fuente: http://www.funnyrobotics.com/2010/10/do-it-yourself-quadrotor-in-pre-order.html
Fuente: http://www.rcgroups.com/forums/showthread.php?t=716870&page=464
15
4. Diseño estético: El atributo estético proporcionado en la fase de diseño a
los componentes estructurales del Quadrotor se contempla como un
requerimiento importante, ya que permite la realización de un producto
atractivo, práctico y funcional, provisto de una excelente terminación y
acabado.
16
1. Elevada capacidad de levante: La ventaja que otorga el transporte de una
adecuada carga útil, permite la adecuación de diferentes tipos de sensores
y dispositivos electrónicos (GPS, altímetro, cámaras de fotografía y video,
etc.) en el interior del Quadrotor, destinados a aplicaciones específicas.
17
5. Vuelo estable y seguro: La estabilidad en el vuelo del Quadrotor se logra
mediante una óptima sintonización de las ganancias en el control PID
(Proporcional Integral Derivativo), función incorporada en la
programación del micro-controlador.
18
4. Espacio reducido de almacenamiento: La facilidad para almacenar y
trasladar cómodamente los componentes estructurales del Quadrotor
requiere que el espacio para realizar estas tareas sea lo más reducido
posible, lo cual se logra con un diseño compacto de su maleta de
transporte.
19
3. Reducción de costos en los procesos de manufactura: La selección de
insumos y materias primas económicas, sumado a un proceso de
manufactura realizado íntegramente a nivel local, contribuye a una
reducción significativa en el costo total de la fabricación del Quadrotor, en
comparación con modelos realizados en el exterior.
20
2. Trazabilidad de componentes: La totalidad de las materias primas e
insumos empleados en la fabricación de la estructura, así como los
equipos y dispositivos electrónicos son productos certificados y
garantizados, adquiridos a proveedores nacionales e internacionales.
Los servicios de soporte técnico que dispone el operador deben contar con los
siguientes requerimientos:
21
3. Capacitación en operación y mantenimiento: Un importante servicio de
soporte técnico son las clases de capacitación y entrenamiento que puede
recibir el operador para la futura ejecución de vuelos. En la Figura 7 se
observa un modelo de Quadrotor denominado QRO (Quad Rotor
Observer) versión 5, desarrollado y construido por JL Naudin Ref. [2],
agregado en el software de simulación para aeromodelismo AeroFly
Professional Deluxe versión 1.8.0.9.
22
2. MARCO TEORICO
23
Figura 8. Torpedo aéreo Kettering Bug
Fuente: Unmanned Aerial Vehicles: Robotic Air Warfare 1917-2007, Steven J Zaloga, Osprey Publishing
Ltd, pág.7, United Kingdom, 2008.
Fuente: Air Force UAV`s The Secret History, Thomas P. Ehrhard, Mitchell Institute Study, pág.5,
United States of America, July 2010.
24
Posteriores desarrollos en este campo produjeron aeronaves más lentas y
económicas, así como funcionales y fáciles de operar en el terreno que fueron
denominadas RPV (Remote Piloted Vehicles) o vehículos pilotados
remotamente que demostraron la aplicabilidad y eficacia del concepto.
Fuente: Enciclopedia Aviones de Guerra “El combate aéreo hoy”, Editorial Planeta-De Agostini S.A.,
Volumen II, Fascículo No.17, pág.323, Barcelona, España, Diciembre 1986.
25
2.1.1 Producción de UAV’s en el mundo
Estadísticas recientes tomadas de Ref. [1] muestran que en los años 2009-10
los sistemas aéreos no tripulados se desarrollaron en 51 fuentes alrededor del
mundo, experimentando un incremento del 22% con respecto a años
anteriores.
26
Los UAV`s con pesos máximos de despegue (Maximun Take-Off Weight,
MTOW) menores a 5 kg ocuparon un 11% del total de desarrollos. La
clasificación de 1190 sistemas implementados en esos años establece que 683
UAV`s (57%) fueron empleados en cometidos militares, 150 (13%) para uso
civil y comercial y 260 (22%) para propósitos mixtos. La mayor proporción de
UAV`s lo constituyen estructuras de ala fija (72%), seguido por vehículos de
ala rotatoria (17%), ducted fan (3%) y sistemas aerostáticos (3%).
Argentina con 13 modelos, México con 4, Brasil con 4 y Chile con 2 son las
principales naciones que lideran proyectos (I+D) surgidos de las necesidades
propias de cada operador en materia de vigilancia y observación. Algunos
países, incluido Colombia han optado por adquirir UAV`s extranjeros,
provenientes de proveedores en Estados Unidos e Israel con el objetivo de
cubrir sus crecientes necesidades de observación y vigilancia más allá del
alcance visual. No obstante, en la actualidad se están desarrollando proyectos
interesantes por parte de la industria nacional en alianza con empresa
privada y algunas instituciones académicas que prometen dar resultados en
el corto y mediano plazo, como el UAV Iris diseñado y construido por la
Corporación de la Industria Aeronáutica de Colombia (CIAC).
27
2.2. CLASIFICACION DE LOS VEHÍCULOS AEREOS NO
TRIPULADOS
Fuente: Unmanned Aerial Vehicles: Robotic Air Warfare 1917-2007, Steven J Zaloga, Osprey Publishing
Ltd, pág.29, United Kingdom, 2008.
28
2.2.1 Según su altitud de vuelo
Fuente: Unmanned Systems Integrated Roadmap 2009-2034, Department of Defense USA, Washington
D.C., pág. 95, March 12, 2009.
Los UAV`s pertenecientes las categorías II y III Ref. [3], requieren certificación
de aeronavegabilidad por parte de las organizaciones de regulación
aeronáutica y su vuelo deber ser monitoreado constantemente por los
servicios de control de tránsito aéreo (Air Traffic Control, ATC), pues operan
a altitudes de vuelo similares de las aeronaves tripuladas. Por su parte, los
UAV’s situados en la categoría I son considerados aeromodelos de vuelo libre
que operan en espacios aéreos no controlados Ref. [4].
Ref.[3] Code of Federal Regulations, Title 14, Aeronautics and Space Part 91,101,103, Right of Way
Rules, FAA, Washington D.C., 2010.
Ref.[4] Advisory Circular AC 91-57 “Model Aircraft Operating Standards”, FAA, Washington D.C.,
June 9, 1981.
29
2.2.2 Según su peso máximo en decolaje, altitud y autonomía de vuelo
Fuente: “Unmanned Aerial Vehicles: An Overview”, María de Fátima Bento, insideGNSS, págs. 54-61,
January-February, 2008.
30
2.2.3 Según sus características de vuelo
Fuente: Ref.[5]
Ref.[5] Unmanned Aircraft Systems: UAV’s Design, Development and Deployment, Reg Austin, Wiley
Publication, págs. 34-42, United Kingdom, 2010.
31
- Tail-aft on Fuselage: Estabilizadores ubicados en la cola del fuselaje.
- Tail-aft on Booms: Fuselaje de cola dividido en dos secciones.
- Flying-Wing: Ala y fuselaje integrados.
- Canard: Estabilizador horizontal ubicado delante del ala.
- Delta: Ala en configuración triangular.
Fuente: Ref.[5]
32
- Single Rotor: Un rotor principal y un rotor de cola anti-torque.
- Co-axial Rotor: Dos rotores principales giran en sentido contrario.
- Tándem Rotor: Dos rotores principales ubicados lado a lado.
- Quad Rotor: Cuatro rotores montados sobre la estructura.
Fuente: Ref.[5]
33
- Tilt-Rotor: Rotores montados en las alas giran en vuelo.
- Tilt-Wing: Alas con rotores giran en vuelo.
- Tilt-Wing-Body: Fuselaje con rotor ubicado verticalmente.
- Ducted Fan: El rotor encapsulado en el interior del fuselaje.
34
Figura 16. Helicóptero Oehmichen No.2
Fuente: http://ingaeronautica.files.wordpress.com/2010/11/helicoptero-2.jpg?w=590&h=272
Fuente: http://illumin.usc.edu/assets/media/152/4589757501_315a2bc282_o.jpg
36
Figura 18. Componentes electrónicos en un Quadrotor
Fuente: http://static.rcgroups.net/forums/attachments/9/8/4/9/7/a2928058-15-OMMmkQuad.jpg
37
Figura 19. Quadrotor AirRobot AR-100-B
Fuente: http://www.minimumsecurity.net/blog/surveils.jpg
Ref.[6] Cotte Corredor Jorge Mario, et al. Diseño de control robusto de velocidad de motores brushless
para robótica aérea, Universidad Nacional de Colombia, 2010.
Ref.[7] Jáuregui Uribe Alejandra Lucía, et al. Diseño de un sistema de control de cabeceo y alabeo de un
vehículo aéreo no tripulado (UAV), Universidad Industrial de Santander, 2010.
Ref.[8] Peña Giraldo Mauricio Vladimir, et al. Modelamiento Dinámico y control LQR de un Quadrotor,
Universidad Libre, 2010.
38
Tabla 4. Clasificación de los Multi-rotores
NOMBRE CARACTERÍSTICAS EJEMPLO
Número de Motores: 3
Número de Rotores: 3
Tricopter
Ángulo de separación: 120º
Configuración: Triangular
Número de Motores: 4
Número de Rotores: 4
Quadcopter
Ángulo de separación: 90º
Configuración: Cruz, Equis
Número de Motores: 6
Número de Rotores: 6
Hexacopter
Ángulo de separación: 60º
Configuración: Hexagonal
Número de Motores: 6
Número de Rotores: 6
Y6copter
Ángulo de separación: 120º
Configuración: Triangular
Número de Motores: 8
Número de Rotores: 8
Octocopter
Ángulo de separación: 45º
Configuración: Octogonal
39
3. METODOLOGIA DE DISEÑO
40
3.1.1 Configuraciones básicas
41
2. Configuración en Equis: Otra forma de estructura, está definida por la
forma en equis (X), en la cual la dirección de vuelo no coincide con los ejes
de los brazos del Quadrotor, encontrándose rotada un ángulo de 45º con
respecto al centro del fuselaje (Figura 22).
1 3 2 4 (1)
M y 0 M x 0
F1 (r ) F3 (r ) 0 (2) F2 (r ) F4 (r ) 0 (3)
F1 F3 F2 F4
43
3.1.3 Direcciones de Vuelo
r ( F1 F3 ) (4)
44
Figura 24. Movimiento de Cabeceo hacia adelante
r ( F3 F1 ) (5)
45
Figura 25. Movimiento de cabeceo hacia atrás
46
2. Alabeo (Roll): La actitud de vuelo en alabeo o roll concede al Quadrotor
la capacidad de efectuar desplazamientos laterales hacia la derecha o
izquierda de su eje central x.
r ( F2 F4 ) (6)
47
Figura 27. Movimiento de alabeo hacia la derecha
48
Figura 28. Movimiento de alabeo hacia la izquierda
49
3. Guiñada (Yaw): El control direccional o movimiento de guiñada en el
Quadrotor se obtiene al girar este sobre su eje vertical z.
4
i (8)
i 1
50
Figura 30. Movimiento de guiñada hacia la derecha
51
Figura 31. Movimiento de guiñada hacia la izquierda
54
Figura 35. Modo de vuelo en estacionario
55
3.2 DISEÑO CONCEPTUAL DEL QUADROTOR
Fuente: Aircraft Design Projects for Engineering Students, Lloyd R. Jenkinson, Chapter 1 Design
Methodology, Fig. 1.1, pág. 2, Butterworth Heinemann.
56
A continuación, se explican cada una de estas fases aplicadas en el diseño del
Quadrotor propuesto:
57
Las características básicas que se tuvieron en cuenta para la realización de los
primeros dibujos de diseño se detallan en la Tabla 12.
58
3.2.2 Diseño del proyecto
59
Al diseñarse el Quadrotor con estructura tipo cerrada, todos estos equipos
deben instalarse en el interior del cuerpo central, protegidos de agentes
externos, así como de golpes y caídas. Con el fin de facilitar los procesos
de manufactura en relación a la fabricación de moldes, contra-moldes,
piezas estructurales y su ensamblado, el cuerpo central fue diseñado en
dos (2) secciones: una superior (Figura 39A) y otra inferior (Figura 39B).
Figura 39. Diseño 3D de las secciones superior e inferior del cuerpo central
61
La sección superior del brazo es un volumen rectangular que presenta
mayor grosor en el plano anterior donde se acopla con el cuerpo central
del Quadrotor, hacia la parte posterior el contorno se adelgaza y termina
en una forma cilíndrica que sirve de carenado al motor eléctrico. La
sección inferior del brazo presenta las mismas características, pero la
forma cilíndrica externa es de menor altura.
62
Figura 43. Ensamble 3D de los componentes superiores e inferiores
63
Obtenido este conjunto tridimensional, el equipo diseñador pudo cerciorar la
forma básica que adopta la estructura del vehículo, prestando atención a su
dimensionamiento, geometría y volumetría general desde diferentes ángulos
y planos de referencia con el fin de validar los trabajos de diseño efectuados
hasta este punto.
64
3.2.3 Diseño detallado
Figura 46. Diseño detallado de las secciones superior e inferior del fuselaje
65
La sección inferior del fuselaje recibió el mismo tratamiento que la sección
superior, pero en el plano circular fue hecha una abertura rectangular que
tiene como objetivo la facilidad de instalación y desinstalación de la
batería Litio-Polímero para propósitos de recarga, configuración que se
observa en la Figura 46B.
66
2. Brazo: De igual forma, la estructura de las secciones superior e inferior del
brazo fueron aligeradas por medio de operaciones de vaciado hasta
obtener un espesor máximo de 2mm. En la Figura 48A se muestra la
forma final de la sección superior en donde se observa el cuerpo
cilíndrico externo que se convierte en una forma carenada en cuyo interior
se aloja el motor eléctrico. A la sección inferior del brazo, visible en la
Figura 48B le fue añadido un elemento cónico que sirve como estructura
de soporte en tierra o tren de aterrizaje.
Figura 48. Diseño detallado de las secciones superior e inferior del brazo
67
Figura 49. Diseño detallado del ensamblaje del brazo
68
Figura 51. Ensamble 3D de la estructura del Quadrotor
69
3. Quadrotor: Una vez se realizados y comprobados los sub-ensamblajes y
ensamblajes de los componentes estructurales, se procedió a instalar
virtualmente los equipos eléctricos y electrónicos internos, spinner`s y
hélices para conformar el conjunto final del Quadrotor, ejercicio que se
observa en la Figura 53.
70
Figura 54. Vista superior e inferior del Quadrotor
71
En el diseño de la disposición interna de los equipos eléctricos y electrónicos
se tuvo en cuenta su agrupación de forma compacta y de fácil instalación en
dos conjuntos principales ubicados en el cuerpo central y los extremos de los
brazos como lo muestra el corte transversal de la estructura del Quadrotor en
la Figura 56.
El cuerpo central se diseñó como una estructura cerrada y hueca que posee
un volumen interno suficientemente grande para la adecuación de equipos
que sirven para su funcionamiento. La Figura 57 detalla los compartimientos
del cuerpo central, los cuales fueron divididos de la siguiente manera:
72
Figura 57. Compartimientos del cuerpo central
73
4. MODELADO MATEMATICO DEL QUADROTOR
Con el fin de realizar la simulacion del sitema dinamico del vehiculo, se hace
necesario la formulacion matematica para el analisis de la implementacion de
los movimientos de control.
Las referencias [9], [10], y [11] describen estas condiciones de diseño, de las
cuales se ha tomado nota para realizar el correcto modelado matemático.
74
4.2 MODELAMIENTO DEL QUADROTOR
75
El momento que genera el movimiento de alabeo o de roll (ángulo ) se
76
En la Tabla 14 se describen tales efectos, donde C representa términos
constantes, Ω es la velocidad del rotor, JR es el momento de inercia rotacional
del rotor alrededor de su eje, l es la distancia del centro de masa a los rotores,
J es el momento de inercia del cuerpo rígido y , son los ángulos de Euler.
centro de masa del Quadrotor con respecto al sistema inercial . Así mismo, la
78
4.3.1 Ángulos de Euler y Matrices de rotación
Los ángulos de Euler son tres coordenadas angulares que permiten relacionar
la orientación de un sistema de ejes respecto a otro así. La Figura 60 muestra
al sistema de referencia determinado mediante estos tres ángulos que
permiten los movimientos del vehículo de la siguiente manera:.
79
Las matrices de rotación son necesarias para convertir de un sistema a otro
mediante rotaciones alrededor de los ejes de coordenadas. Un ejemplo de esto
es cuando en un momento dado representamos la posición de un satélite
artificial en el sistema XYZ. Un tiempo t después, la posición del satélite
habrá variado debido al movimiento en su órbita, pero también la Tierra
habrá rotado, de modo que tendremos un nuevo sistema, que no es igual al
sistema inicial XYZ´, (Ref.[16])
Fuente: http://nacc.upc.es/nacc-libro/img/ejemplo-ecef-satelite.png
80
Considerando los dos sistemas de coordenadas mostrados en la Figura 61, en
donde sus ejes Z y C coinciden y ix, jY, iA, jB corresponden a los vectores
unitarios de su respectivo sistema de coordenadas, es posible representar un
vector P de la siguiente forma:
(9)
(10)
Igualando las ecuaciones 9 y 10, así como realizando el producto punto con ix,
jy, kz en ambos lados de la ecuación y organizando el resultado en una
matricial, se obtiene:
(11)
81
Donde es el vector definido en el sistema , es el vector
).
El determinante es 1 ( = 1).
(13)
82
Figura 63. Rotación alrededor del eje x
(14)
(15)
83
(16)
(17)
por:
(18)
84
A partir de la matriz de rotación (17) generada por las tres rotaciones
sucesivas y su propiedad de orto-normalidad, relacionando la derivada de la
matriz orto-normal con una cierta matriz anti-simétrica (ref. [18]), se puede
obtener las ecuaciones cinemáticas de rotación del vehículo que establecen las
relaciones entre las velocidades angulares.
(19)
(20)
Definiendo:
(21)
Ref. [17] S. Bouabdallah, A.˜ Noth, y R. Siegwart. Dynamic Modelling of UAVs. Technical Report
Version 2.0, Autonomous Systems Laboratory - EPFL, May 2006. http://www.sky-
sailor.ethz.ch/docs/Dynamic_Modeling_of_Fixed-Wing_UAVs_(12.05.2006).pdf
Ref. [18] J. J. Craig. Introduction to Robotics - Mechanics y Control. Addison-Wesley Publishing
Company, Inc., USA, 2nd edn edition, 1989. http://es.scribd.com/doc/52377374/Introduction-to-
robotics-mechanics-and-control-2nd-ed-1986-John-J-Craig
85
Por lo tanto, las ecuaciones cinemáticas para determinar la postura del
Quadrotor, suponiendo la matriz de rotación (17), vienen dadas por:
(22)
(23)
(24)
(25)
86
El movimiento rotacional del Quadrotor viene dado por las componentes de
las velocidades angulares en los tres ejes: velocidad angular de alabeo
Ref. [19] S. Esteban. Control Nolineal en Tres Escalas de Tiempo de un Helicóptero de Radio/ Control
sobre una Plataforma. Technical report, Universidad de Sevilla, Sevilla, ES, 2005.
http://bibing.us.es/proyectos/abreproy/70017/fichero/Tesis_Master_GuilhermeRaffo.pdf
87
4.4 FORMULACIÓN DE NEWTON-EULER
(27)
(28)
88
(29)
Donde y
(30)
Con las ecuaciones de fuerzas y pares (30), el modelo dinámico (29) se puede
reescribir como sigue:
(31)
, donde
89
w es la velocidad del Quadrotor con respecto al aire) (Ref. [20]).
es la constante gravitacional ( )
(32)
constante.
Ref. [20] V. Mistler, A. Benallegue, y N. K. M’Sirdi. Exact linearization and noninteracting control of a 4
rotors helicopter via dynamic feedback. In Proc. IEEE Int. Workshop on Robot and Human Inter.
Communic., 2001.
90
Así, a través de la segunda Ley de Newton se obtiene (Ref. [11]):
(33)
(34)
(35)
91
El eje de rotación del motor se mueve con la velocidad angular del eje de
referencia, lo cual produce los siguientes movimientos giroscópicos:
(36)
Así, usando las ecuaciones (23), (27), (32) y (35), y definiendo el nuevo vector
de estados como (Ref. [20]):
(37)
(38)
92
Este modelo matemático puede asumirse como suficientemente preciso en la
representación de todos los movimientos funcionales de un vehículo aéreo
autónomo (Ref. [16]). Sin embargo, no es adecuado para el diseño del
controlador porque éste depende de fuerzas y momentos aerodinámicos, los
cuales son desconocidos. En consecuencia, se tienen en cuenta las siguientes
consideraciones para justificar la eliminación de estos términos y sean
considerados como perturbaciones externas:
(39)
93
La estabilidad del Quadrotor está ligada a mantener en un punto fijo los
ángulos de Euler y estabilizar su altura. Se ha decidido suprimir las
ecuaciones que describen la posición X e Y para facilitar la
implementación de un filtro de control optimo que cumpla con estas
expectativas.
(40)
Deduciendo que U1, U2, U3, y U4 son las fuerzas de empuje en la aceleración,
alabeo, cabeceo y guiñada respectivamente, tenemos:
94
(41)
(42)
95
4.5 GEOMETRÍAS Y CÁLCULOS DE LOS MOMENTOS DE INERCIA
En esta sección se utiliza los datos que se obtienen al recolectar las medidas
básicas de los componentes principales, teniendo en cuenta la teoría de los
fundamentos del momento de inercia descritos en el Anexo A.
Aquí denota el momento de inercia alrededor del eje cuando los objetos
giran alrededor del eje x, denota el momento de inercia alrededor del eje y
Es por esto que el Quadrotor debe ser modelado como varios componentes
con geometrías más sencillas, (ver Figura 64) con la excepción de la estructura
principal que está compuesta por fuselaje y brazos (objetivo principal en este
proyecto); el cual se modeló con la herramienta de diseño mecánico
SolidWorks 2010, aprovechándose la versatilidad de este software para
calcular los momentos de inercia para esta los estructurales que presentan
una forma geométrica compleja. Para esto, se ejecutó una buena
aproximación definiendo:
96
Cuatro 4 hélices, modeladas como cilindros planos con densidad variable
(aumentando a medida que se acerca al centro de la misma).
Cuatro acoples o núcleos de madera utilizados para la unión del brazo con
el fuselaje central, modelados como un paralelepípedo rectangular sólido.
97
La estructura tiene un espesor de pared de 2mm, distancia entre el centro y
cada brazo de 300 mm y una altura de 84 mm; densidad de la fibra de
carbono de 1000 Kg/m3. Luego de ensamblar sus componentes estructurales
se determinan los momentos de inercia a través de la herramienta de
propiedades físicas de SolidWorks; arrojando los siguientes reultados:
ICx= 0.0127
ICx= 0.0254
98
4.5.2 Momento de inercia de los cuatro motores
Un radio
Una altura de
Una masa
99
Los momentos de inercia del motor alrededor del eje X, , alrededor del
(43)
(44)
(45)
(46)
(47)
(48)
Las cuatro hélices, como los anteriores componentes, muestran una completa
simetría, por tanto, el cálculo se realiza sólo en la parte frontal la hélice y
luego incorporado en los otros. Se han modelado como cilindros planos, con
densidad variable.
100
Esta forma ha sido elegida ya que una placa plana se asemeja a la forma de la
hélice, además, la densidad depende del radio y disminuye a medida que se
aleja del centro. El cilindro plano del modelo de la hélice (Figura 67) presenta
las siguientes características:
Radio
Altura
Una masa .
alrededor del eje Z, , son definidos en las ecuaciones (49), (50) y (51):
101
(49)
(50)
(51)
(52)
(53)
(54)
Altura ,
Longitud ,
Ancho ,
La masa de .
102
Figura 68. Modelado de la batería como una caja rectangular maciza
alrededor del eje Z, , son definidos en las ecuaciones (55), (56) y (57):
(55)
(56)
(57)
Esta parte es un sistema de acople que une el brazo con el fuselaje, fabricado
en madera con una forma irregular, se modela en forma de caja rectangular
maciza de:
103
Altura
Longitud
Ancho
masa de
alrededor del eje Z, , son definidos en las ecuaciones (58), (59) y (60).
(58)
104
(59)
(60)
(61)
(62)
(63)
Un radio
Una altura de
Una masa
105
Figura 70. Geometría del modelado del cilindro de los componentes
electrónicos.
(64)
(65)
(66)
106
Por último es posible calcular el momento de inercia del cuerpo tres , y
definido.
(67)
(68)
(69)
107
5. SELECCIÓN Y CARACTERIZACION DE HARDWARE
Los motores eléctricos más empleados en este campo son los tipo Brushless
(Figura 71), que funcionan en AC debido a que se alimentan con una señal
trifásica (Por este motivo requieren controladores de velocidad), sin embargo
se clasifica como un motor de DC porque al igual que los motores comunes
tienen imanes permanentes. Estos motores poseen la principal característica
de que no emplean escobillas en la conmutación para la transferencia de
energía, ya que estas producen rozamiento, disminuyen el rendimiento,
generan calor, son ruidosos y demandan una sustitución periódica y por
tanto, un mayor mantenimiento.
Ref. [21] Cristóbal Ángel. Técnica: Motores eléctricos brushless, Revista AEROTEC, Número 108, pág.
22, RDM Editorial, Madrid, España, 2003.
108
Figura 71. Componentes internos de un motor eléctrico brushless TM 350-28
Fuente: Revista ELECTRIC FLIGHT UK, pág. 9, Edición 72, Reino Unido, 2003.
109
5.5.1 Parámetro Kv en un motor Brushless:
Esta constante significa la cantidad de vueltas (rpm) que da el motor por cada
voltio aplicado al ESC (A máxima potencia). Por ejemplo, si a un motor que
posee un Kv de 1100 se le aplica una tensión de 11,1v operará a un valor
máximo de revoluciones de:
110
Una forma fácil para determinar el motor adecuado es utilizando
herramientas de software o programas en línea que algunos fabricantes
proporcionan en bases de datos con referencias y especificaciones que
permiten el cálculo de la eficiencia y el empuje aproximado dependiendo de
la variación de los elementos que componen el sistema tales como hélices,
ESC’s y baterías. Una muestra de ello se observa en el software gratuito
llamado Drivecalc versión 3.4, herramienta para el análisis de motores DC y
el cálculo completo de la potencia de sistemas para aeromodelismo con motor
eléctrico. Contiene una base de datos con tablas para baterías, ESC, cajas
reductoras, motores, hélices y estatores. En la Figura 72 se observa la interfaz
del programa y los posibles datos que se incluyeron para el análisis de
equipos.
111
Figura 72. Cálculos en el software Drive Calculator 3.4
Fuente: http://www.drivecalc.de/
112
En primer lugar se debe elegir el tipo de motor a utilizar, que para este caso
corresponde al motor Brushlees tipo outrunner, preferido por sus grandes
ventajas respecto a los motores inrunner que proporcionan un par muy bajo
por lo cual se descartaron.
Fuente: http://www.electrifly.com/powersystem/powersystem-index.html
113
Figura 74. Etapa de ingreso del tamaño de hélice
Fuente: http://www.electrifly.com/config/index.php
Fuente: http://www.electrifly.com/config/index.php
114
El motor seleccionado para equipar el Quadrotor es un motor brushless
outrunner, marca RIMFIRE (30-35-1250) modelo .10 que presenta las
siguientes dimensiones (Figura 75) y especificaciones técnicas.
Fuente: http://www.electrifly.com/motors/gpmg4595.jpg
ESPECIFICACIONES TECNICAS
Kv 1250 rpm/volt
Empuje 3d 2 lb (905g)
115
Para determinar el comportamiento del grupo moto-propulsor se realizaron
mediciones para la obtención de la relación fuerza (empuje) contra la señal
del receptor (PPM). Debido a que la fuerza ejercida por cada motor está
ligada a la señal del receptor, que en este caso particular es la señal de
control, se hizo necesaria la construcción de un banco de pruebas que
facilitara esta toma de medidas, el cual se muestra en la Figura 76.
1. Procedimiento
117
En cada una de las posiciones de incremento en la palanca de aceleración, se
produce una señal PWM que ingresa al control de velocidad electrónica
(ESC). Estas señales, tanto los valores correspondientes a la fuerza de empuje
producida por el grupo moto-propulsor fueron listados en la Tabla 19.
Por medio del software Curve Expert 1.4, se desarrolló la ecuación que
caracteriza el comportamiento de las variables de empuje con respecto a la
señal de control PPM que se mide en la interfaz gráfica de acuerdo a la
posición de la palanca del acelerador. Al ingresar los datos de las mediciones
en el programa se observa un comportamiento casi lineal en la curva, por lo
que una regresión lineal simplifica la ecuación del comportamiento de las
variables con el fin de incluirla en las ecuaciones del modelado matemático y
facilitar el diseño del controlador.
118
Tabla 19. Resultados de las pruebas en el banco de caracterización del motor.
119
Por medio de un ajuste polinomial de tercer orden se obtiene la siguiente
relación:
Datos de coeficientes:
120
Arrojando la siguiente ecuación lineal:
Ajuste lineal:
Datos de coeficientes:
(72)
el modelado matemático.
121
5.2 CARACTERIZACIÓN Y SELECCIÓN DE LAS HELICES
Los perfiles que componen una hélice están sujetos a las mismas leyes y
principios que cualquier otro perfil aerodinámico, por ejemplo un ala. Cada
uno de los cuales tiene un ángulo de ataque, respecto al viento relativo de la
pala que en este caso es cercano al plano de revolución de la hélice, y un paso
igual al ángulo de incidencia (Figura 79). El giro de la hélice, acelera el flujo
de aire hacia el borde de salida de cada perfil, a la vez que deflecta este hacia
atrás (lo mismo que sucede en un ala).
Este proceso da lugar a la aceleración hacia atrás de una gran masa de aire,
movimiento que provoca una fuerza de reacción que es la que propulsa la
aeronave hacia adelante (Ref.[22]).
122
Figura 79. Ángulos de la hélice
Fuente: http://www.manualvuelo.com/GIFS/Fig_322.gif
Fuente:
123
5.3 CARACTERIZACIÓN Y SELECCIÓN DE LOS CONTROLADORES
DE VELOCIDAD
124
La mayoría de los variadores incorporan un dispositivo denominado BEC
(Battery Elimination Circuit), se trata de un regulador lineal de tensión que
suministra al receptor y los servos; en este caso particular alimenta la tarjeta
Arduino con 5V, de este modo se elimina la necesidad de otra batería para el
sistema de control o una etapa de regulación para la misma, representando
un ahorro en el peso del vehículo.
Fuente: http://www.electrifly.com/escs/gpmm1830.jpg
125
En la Tabla 20 se detallan las especificaciones técnicas del controlador
seleccionado:
Fuente: http://manuals.hobbico.com/gpm/gpmm1800-1841-manual-v1_2.pdf
126
5.4 CARACTERIZACIÓN Y SELECCIÓN DE LA BATERIA
Uno de los problemas que han tenido los modelos eléctricos ha sido la corta
duración de tiempo de trabajo y por tanto la baja autonomía de los aviones
así como la poca potencia que entregaban a los motores y su elevado peso. En
la actualidad existen baterías que son capaces de mover aviones de varios
kilos y con autonomías similares a los motores de combustión, siendo su peso
mucho menor.
A pesar de estos grandes beneficios existen puntos negativos para este tipo de
batería, ya que son costosas en comparación con las de NiCad y NiMH.
Además de problemas de seguridad debido a que el electrolito es muy volátil,
produciendo incendios o explosiones.
127
Las diferencias entre las baterías de iones de litio (Li-Ion) y las de polímero de
litio (LiPo) están en cómo las células se empaquetan y el tipo de electrolito
que se utiliza. En la Figura 82 se puede observar la composición y disposición
de las celdas en forma apilada, conectores y plug para el balanceo de cada
celda de una batería de LiPo.
Fuente: http://traxxas.com/sites/default/files/images/2872-lipo-exploded-view_m.jpg
128
5.4.2 Baterías LiPo
• Cada celda de una batería NiMH tiene 1,2 V, por lo que una batería de 8
celdas tiene 1.2x8 = 9.6V
• Cada celda de una batería LiPo tiene 3,7 V, por lo que una batería LiPo
de 3 celdas tiene 3.7x3 = 11.1V
129
La selección de la batería se realizó basándose en dos aspectos, el tiempo de
vuelo y la corriente máxima necesaria, variables relacionadas con el motor.
Según las especificaciones del motor Rimfire .10, la corriente de consumo es
de 18.1 [A], lo cual produce un total de 72.4 A, sumatoria por los 4 motores
que la batería necesita ser capaz de suministrar. La corriente máxima que una
batería puede entregar está determinada por:
(73)
Donde:
La corriente nominal o capacidad de la batería;
La capacidad de descarga.
(74)
130
Para estimar la capacidad mínima nominal es necesario el uso de la corriente
para estado estacionario (43.24 A) y el valor del tiempo de vuelo mínimo
requerida:
(75)
Este análisis indicaría que se necesita una batería mayor a 5000 mA, pero
teniendo en cuenta que al aumentar la capacidad de la batería aumenta el
peso de la misma, para satisfacer la durabilidad estimada se haría necesario
disponer de 2 baterías en paralelo; 1 de 5000mA y otra de 2200 mA.
(76)
Que sería tres veces mayor la corriente que demandan los motores en estado
estacionario.
(77)
131
Para la selección de la batería se tuvo en cuenta los anteriores criterios, pero el
criterio más predominante es el precio y la calidad; ya que se escogió una
batería LiPo marca Venom que posee las siguientes especificaciones técnicas:
132
Este dispositivo visible en la Figura 83 dispone de 4 leds de alta eficiencia,
visibles en exterior, así como dos zumbadores sonoros que el operador puede
escuchar a una distancia considerable.
Figura 83. Alarma sonora de baja tensión para baterías LiPo de 3-4S
Fuente:http://img1.mlstatic.com/alarma-sonora-de-bajo-voltaje-para-bateria-lipo-2s-4s_MCO-O-
19722213_1346.jpg
La alarma está ajustada para que se activen los leds y el zumbador a los
niveles mínimos de la batería que recomiendan los fabricantes: (2S-6,3V) (3S-
9,45V) y 3,15V por cada celda de más de 3S.
133
5.5 CARACTERIZACIÓN Y SELECCIÓN DEL RADIOCONTROL
Fuente: http://www.megastorerc.com/1092-1550-thickbox/receptor-futaba-r-617fs--7ch--24ghz.jpg
134
5.5.1 Amplitud modulada (AM)
135
Figura. 85 Estructura de un pulso PPM
Fuente: http://static.rcgroups.net/forums/attachments/2/4/1/5/3/a4901379-164-frame_ppm.gif
136
Tabla 22. Especificaciones técnicas del sistema radio-control Futaba.
Consumo de
Requisito
corriente 170mA 4,8 - 6V
de energía
Consumo
80 mA (sin señal)
de corriente
137
5.6 CARACTERIZACIÓN Y SELECCIÓN DEL MICROCONTROLADOR
Fuente: http://arduino.cc/en/uploads/Main/ArduinoUno_R3_Front_450px.jpg
138
Las especificaciones técnicas y de operación del micro-controlador Arduino
UNO se pueden ver resumidas en la Tabla 23. El datasheet* del micro-
controlador utilizado, donde se puede observar la asignación de los pines,
características térmicas, eléctricas, además de los planos del micro-
controlador ATMEGA 328.
MICROCONTROLADOR ATMEGA328
Voltaje de operación 5V
Voltaje de entrada (recomendado) 7-12V
Voltaje de entrada (limites) 6-20V
Pines Digitales I/O 14 (de los cuales proveen 6 salidas PWM)
Pines de entradas análogas 6
Pin de corriente DC I/O 40 mA
Pin de corriente DC para 3.3V 50 mA
32 KB (ATmega328) de los cuales 0.5 KB usados
Memoria Flash
para el gestor de arranque
SRAM 2 KB (ATmega328)
EEPROM 1 KB (ATmega328)
Velocidad de reloj 16 MHz
* http://www.atmel.com/Images/8271s.pdf
139
Contiene todo lo necesario para utilizar el micro-controlador; simplemente
conectándolo a un ordenador a través del cable USB o aliméntalo a través de
una batería. El Arduino, tiene un regulador de voltaje, que además de poder
ser alimentado con un mínimo de 6 voltios y un máximo de 20 voltios tiene 2
salidas de voltaje para alimentar sensores, una de 5 V y otra de 3.3 V. La que
se usa para la IMU es la segunda opción.
Serial: 0 (RX) and 1 (TX). Son usados para recibir (Rx) y transmitir (Tx)
datos de señal TTL.
Interruptores externos: 2 y 3. Estos pines pueden ser configurados para
provocar un valor alto o bajo, un flanco ascendente o descendente, o un
cambio de valor.
PWM: 3, 5, 6, 9, 10, and 11. Proveen salidas PWM de 8-bit con la
función analog Write.
SPI: 10 (SS), 11 (MOSI), 12 (MISO), 13 (SCK). Estos pines soportan la
comunicación SPI usando la librería SPI.
6 entradas análogas, etiquetadas desde la A0 hasta la A5, cada una
provee una resolución d 10 bits (valores diferentes 1024 )
TWI: pin A4 o SDA y pin A5 o SCL. Soporta comunicación TWI usando
la librería Wire.
140
2. Comunicación: La tarjeta micro-controladora presenta una serie de
facilidades para comunicarse con un ordenador, otro Arduino y otros
micro-controladores. El ATmega328 ofrece UART TTL (5V) de
comunicación en serie, que está disponible en los pines digitales 0 (Rx) y 1
(Tx). Un ATmega16U2 en los canales de comunicación esta en serie a
través de USB y aparece como un puerto COM virtual con el software en
el ordenador. El software de Arduino incluye un monitor serial que
permite el envío de datos de texto desde y hacia la tarjeta. Los Leds Rx y
Tx en la placa parpadean cuando se están transmitiendo datos a través de
la conexión de USB en serie y la computadora.
141
El entorno más utilizado es el software Processing, que corresponde a la
plataforma utilizada en este proyecto para la interacción entre el
comportamiento de los sensores, calibración y ajuste de las principales
variables del control.
142
Figura 87. Interfaz gráfica del sistema Quadrotor en processing
143
f) Calcular la acción de control de los movimientos de navegación y
estabilidad de los mismos mediante la implementación de un filtro
controlador tipo PID.
g) Comandar a los variadores de velocidad electrónicos (ESC’s) mediante la
acción de control calculada anteriormente a través de las salidas PWM.
h) Enviar y recibir todos los datos obtenidos por los sensores, calibrar
sensores y ajustar las ganancias del controlador a través del puerto serial
en comunicación con la interfaz codificada en Processing.
144
En el mercado se encuentran IMU`s de diferentes configuraciones, desde 3
Grados de libertad (DOF) hasta 9 DOF. El más común es el de 6 DOF que
consta de un acelerómetro de 3 ejes y un giróscopo de 3 ejes; de que se
obtienen mediciones exactas de los ángulos de navegación cabeceo y alabeo;
además de poder tener una medición no referenciada del ángulo de guiñada.
APLICACIÓN RANGO
Navegación 4 a 15 mV/dps
Fuente:http://www.scantec.de/fileadmin/pdf/Invensense/Selection-and-integration-of-
MEMS-based-motion-processing-in-consumer-apps-070809-EE-Times.pdf
145
La selección de la IMU se realizó teniendo en cuenta los parámetros
anteriormente mencionados como el rango de sensibilidad para aplicaciones
de navegabilidad y junto con la limitante de la disponibilidad en el mercado
nacional, información que se muestra en la Tabla 25.
Especificaciones
Dispositivo Rango Interface Ejes Alimentación
especiales
IMU Razor 9
DOF Compatible con
Acel:±16g Acel: 3 Arduino, código
serial 3.3-16VDC
Giro:±300°/s Giro: 3 opcional AHRS,
magnetómetro de 3 ejes.
IMU Fusion –
ADXL345 y
IMU
Acel:±2, 4, 8, 16g Soporte de hardware
3000 Acel: 3
Giro:±250- I2C 3.3-16VDC para los datos Motion
Giro: 3
2000°/s Fusion
Procesador Digital
MPU Motion ™ (DMP ™),
6050 bus I2C auxiliar para
Acel:±2, 4, 8, 16g
Acel: 3 magnetómetros
Giro:±250, 500, I2C 2.3 - 3.4VDC
Giro: 3 externos y otros
1000, 2000°/s
sensores, salida digital
del sensor de
temperatura
146
En este proyecto se seleccionó una IMU Fusion de la marca Sparkfun de 6
grados de libertad que está compuesta por un acelerómetro de referencia
ADXL345 de 3 ejes y un giróscopo IMU3000 de 3 ejes, dispositivo que se
puede observar en la Figura 88.
Fuente: https://dlnmh9ip6v2uc.cloudfront.net/images/products/1/0/2/5/2/10252-04_i_ma.jpg
147
Dentro de los más importantes encontramos los de última tecnología que
son los MEMS (Microelectromechanical Systems), distinguiendo tres
categorías principales de acelerómetros MEMS: el capacitivo de silicio, el
piezorresistivo y, finalmente, los acelerómetros térmicos. Hasta el
momento, los acelerómetros capacitivos de silicio dominaban
ampliamente el mercado. Los fabricantes que usan capacitivos
micromecanizados en la superficie o en masa, como DI, Bosch, Denso,
Freescale, Analog Devices, entre otras.
**http://www.analog.com/static/imported-files/data_sheets/ADXL345.pdf
148
2. Giróscopo: Dispositivo capaz de medir la orientación, basado en los
principios de conservación del momento angular. Cuando se somete a un
momento de fuerza que tienda a cambiar de dirección como lo haría un
cuerpo que no girase, cambia la orientación en una dirección
perpendicular a la dirección intuitiva.
Fuente:http://www.scantec.de/fileadmin/pdf/Invensense/Selection-and-integration-of-MEMS
based-motion-processing-in-consumer-apps-070809-EE-Times.pdf
149
Diapasones vibratorios integran un giróscopo ya que la vibración de las
masas provoca una oscilación de igual magnitud y en sentidos opuestos;
cuando el dispositivo se hace girar, la fuerza de Coriolis crea una fuerza de
vibración ortogonal proporcional a la velocidad de rotación, que se mide " a
lo largo del perímetro de la estructura además estos dispositivos proveen
una alta sensibilidad y requieren menos amplificación, lo que conduce a un
menor ruido.
Los sensores magnéticos son sistemas de medida que utilizan las propiedades
magnéticas o magneto-elásticas de un material, generalmente magnético, que
se suele denominar núcleo del sensor, para obtener una señal eléctrica, de
tensión o corriente, proporcional a la magnitud que se pretende medir,
detectar o controlar. Los sensores magnéticos más comunes son los
magnetorresistivos (MR), que se basa en la propiedad de un material que
cambia su resistividad por la presencia de un campo magnético.
***http://www.sparkfun.com/datasheets/Sensors/IMU/ps-imu-3000a-00-01.1.pdf
150
En la Figura 90 se muestra el principio de funcionamiento básico de un sensor
MR, con una tira de material ferromagnético llamado permalloy (20% de Fe,
80% de Ni). Suponiendo que, cuando no está presente un campo magnético
externo, el permalloy tiene un vector de magnetización interno paralelo al
flujo de la corriente ‘I’ de izquierda a derecha.
Fuente: http://es.scribd.com/doc/3838277/50/Sensores-Magnetorresistivos
151
El magnetómetro seleccionado corresponde al modelo HMC5883L de
Honeywell, el cual es un montaje de la serie magneto-resistivos, que provee
un ADC de 12 bits para mejor precisión. El HMC5883L utiliza un magneto
anisotrópico, tecnología que ofrece ventajas en la precisión ya que caracteriza
una sensibilidad y linealidad en el eje direccional.
Fuente: https://dlnmh9ip6v2uc.cloudfront.net/images/products/1/0/5/3/0/10530-02c_i_ma.jpg
****http://www51.honeywell.com/aero/common/documents/myaerospacecatalog-
documents/Defense_Brochures-documents/HMC5883L_3-Axis_Digital_Compass_IC.pdf
152
Sus características técnicas más importantes son:
153
5.9 CARACTERIZACIÓN Y SELECCIÓN DEL BARÓMETRO
Fuente: http://www.sparkfun.com/tutorial/Barometric/Pressure-Altitude.PNG
154
Figura 93. Sensor de presión tipo piezorresistivo
Fuente: http://tamarugo.cec.uchile.cl/~cabierta/libros/l_herrera/iq54a/especif.htm
*****http://www.adafruit.com/datasheets/BMP085_DataSheet_Rev.1.0_01July2008.pdf
155
6. SELECCIÓN Y CARACTERIZACION DE SOFTWARE
156
Figura 94. Sistema de control de lazo cerrado o retroalimentado.
Fuente: http://www.paginasprodigy.com/5557917723/imagenes/lazo%20cerrado.gif
(79)
Donde:
es la señal de salida del sistema
es la función de transferencia.
; ; y
con
1. Planta de la altitud
(80)
158
2. Planta movimiento de alabeo
(81)
(82)
159
Aplicando la transformada de Laplace, da como resultado:
(83)
Momentos de inercia
Longitud brazo
Constante de empuje Kg
160
Reemplazando los valores de la Tabla 29 en las ecuaciones (80), (81), (82) y
(83), se determinan las plantas de altitud, alabeo-cabeceo-guiñada,
respectivamente:
(84)
(85)
Ahora que se tienen las plantas del vehículo, resta saber cómo es la respuesta
a perturbaciones. En el caso de la planta de los movimientos de navegación
( ) se requiere saber el comportamiento del ángulo (set point=0 grados)
162
6.2 IMPLEMENTACION DE TECNICA FILTRO PID
Dónde:
de transferencia queda:
(86)
(87)
163
La señal de control tiene un valor diferente de cero cuando la señal
(88)
164
Figura 97. Comparación del tiempo de respuesta entre el control P, PI y PID.
Fuente:http://1.bp.blogspot.com/_Hkf1Ma_fPE/R2BuZ77vBI/AAAAAAAAAIU/YYcStloXbdM/s400
/Dibujo4.JPG
Figura 98. Retroalimentación con control PID sin sintonizar para la planta de
altitud.
165
Figura 99. Retroalimentación con PID sin sintonizar para la planta de
movimiento de cabeceo
166
Tabla 29. Métodos usados para la sintonización de un controlador PID.
Ziegler-
Nichols
Los valores para Kp, Ti y Td son determinados como:
Desventaja: Es difícil
Tipo de Kp Ti Td
exactamente establecer el
control
aumento crítico.
p T/L ∞ 0
pi 0.9 T/L L/0.3 0
pid 1.2 T/L 2L 0.5 L
167
Para este modelo se usa una herramienta de auto-sintonización que contiene
Matlab llamada Tune en el bloque PID discrete. Realizando esta operación
que determina automáticamente las ganancias, se cuenta con la posibilidad
de variar el tiempo de respuesta a la que el sistema se estabiliza. En la Figura
100 se observa la aplicación de la auto-sintonización para la planta de altitud.
168
Figura 101. Ganancias de la auto-sintonización de planta altitud en Matlab.
169
La Figura 103 muestra las ganancias para esta respuesta.
170
Se puede observar en el tiempo de 1.5 s como el error proporciona picos de
inestabilidad. En este caso el control logra estabilizarse en un tiempo de 3.5 s;
entendiendo que los errores en la medición afectan enormemente en el
tiempo de respuesta y control de los sistemas. Aplicando un cambio en las
ganancias para este nuevo tipo de error se obtiene (Figura 105):
171
Por medio de la implementación de esta herramienta computacional, se logra
tener una visión más clara acerca del sistema de control que se requiere para
este tipo de sistema.
Debido a que Arduino es basado en C++, se debe seguir una sintaxis básica
para la organización del lenguaje de programación como declarar librerías
que incluyan funciones propias; declarar variables como registros, constantes,
pines, etc., lo cual determina el tipo de datos a trabajar. Además de iniciar la
comunicación serial que es requerida para la transmisión de datos entre la
tarjeta y la computadora.
172
La arquitectura del algoritmo tiene la estructura general de un programa
desplegable en la plataforma Arduino. Contando con dos procedimientos de
entrada que son llamados automáticamente por el planificador del
Framework de Arduino.
173
El programa está basado de una arquitectura por interrupciones, las cuales
son rupturas del flujo del programa principal para la ejecución de una
subrutina de interrupción.
174
1. Módulo QUADMECANICA: En este módulo de definen las variables
principales del control ya que es el encargado de recibir todos los datos de
los demás módulos y de procesar las variables de la estabilización y
calcular la potencia (señal PWM) a cada uno de los motores. En la Figura
107 se ve el despliegue de este modulo en el Framework de Arduino.
Las entradas del controlador PID las forman los tres ángulos calculados
por el módulo IMU, mientras que las salidas son las señales PWM que
comandan a cada uno de los motores. Para cada eje se tiene una única
entrada, que es el ángulo calculado por el módulo IMU, y cuatro salidas
correspondientes a las señales PWM de los motores.
175
La señal de salida a suministrar a los cuatro motores para provocar un
cambio en un eje, se simplifica a un único valor que representa un
diferencial de potencia entre dos parejas de motores. De esta manera,
cuando se calcula un ángulo real (calculado por el módulo IMU) en un eje
concreto, y teniendo en cuenta el ángulo de referencia que seria 0 grados,
se puede calcular el error absoluto entre ambos (que sería su diferencia).
Dicho error es la entrada al controlador, mientras que la salida
representaría el empuje necesario de los motores eliminando dicho error
en forma de diferencia de potencia entre dos pares de motores.
176
4. Módulo Receptor: Este módulo es el encargado de recibir las salidas PPM
del receptor comandas por el transmisor. Dichas lecturas se realizan a
través del uso de Timer’s, para poder obtener todos los canales del
radiocontrol. Se utiliza el Timer (PCINT) en modo de captura para
procesar la información cada vez que exista un flanco de subida en la
señal que indica el final del ciclo de trabajo. Por cada flanco comprobamos
si fue el inicio del marco o si fue el ciclo de trabajo de alguna de las
señales. De esta forma podemos conocer los valores de todos los canales.
177
Además permite la visualización de dicha información en tiempo real. En
el anexo D se cuenta con una guía para el uso de dicha interfaz.
178
7. PROCESOS DE MANUFACTURA
En este capítulo se hace una detallada explicación de los criterios técnicos que
se tuvieron en cuenta para la selección de cada una de las materias primas e
insumos empleados en todas las etapas de fabricación de los componentes de
la estructura del Quadrotor; así como la caracterización de algunas de sus
propiedades físicas, químicas y mecánicas más importantes.
179
1. Madera: Con el propósito de realizar el procedimiento de moldeado
inicial de las piezas estructurales del Quadrotor a través de la técnica
escogida de una forma sencilla, económica y eficaz, fue seleccionada y
empleada madera de cedro dado presenta una excelente estabilidad
dimensional y buenas propiedades mecánicas requeridas para este trabajo
de elaboración de formas y contornos. Para ejecutar esta tarea se cortaron
y adecuaron dos (2) bloques rectangulares y dos (2) bloques cuadrados de
madera de cedro que correspondían a los brazos y fuselaje del Quadrotor.
120 Medio
240 Fino
400 Muy Fino
180
Figura 109. Sellador y laca poliuretano
Fuente: http://www.madecentro.com/modulos/Productos
181
7.1.2 Materiales empleados para la fabricación de contra-moldes
182
Los hilos en fibra de vidrio que están unidos entre sí por medio de una
emulsión soluble en estireno, presentan una óptima compatibilidad con
todo tipo de resinas poliéster, además es fácil de trabajar y permite una
buena supresión de burbujas durante la aplicación de la emulsión. La
ficha técnica de la lámina adquirida se presenta en la Tabla 31.
Fuente: http://tienda.resineco.com/es/PF00029
184
La capa de alcohol polivinílico debe ser aplicada uniformemente sobre
toda la superficie del molde expuesta a la resina y cuando se seca, se
presenta la formación de una fina película que no permite la adherencia
de la emulsión.
185
Figura 112. Resina poliéster CRISTALAN ® Referencia 856
186
Tabla 34. Ficha técnica de la resina poliéster
NOMBRE COMERCIAL Resina Poliéster Cristalán Referencia 856
PRODUCTOR ANDERCOL S.A.
Resina Poliéster: 98%
COMPONENTES
Naftenato de Cobalto: 2%
Resina: 98%
RELACION DE MEZCLA Estireno: La proporción depende de la viscosidad deseada
Peróxido: 1% - 2%
Los tiempos de endurecimiento dependen de la proporción de los
ENDURECIMIENTO agentes acelerantes (Cobalto-MEK) y la temperatura ambiente. No debe
utilizarse a temperaturas inferiores a 15ºC.
Tanto el líquido como sus vapores son inflamables, nocivos si son
ingeridos, inhalados o absorbidos a través de la piel o mucosas. Produce
RIESGOS
irritación en la piel, ojos y tracto respiratorio. Afecta el sistema nervioso
central, hígado y sistema reproductivo.
No es reciclable, ni biodegradable. Se debe evitar su desecho en fuentes
IMPACTO AMBIENTAL
de naturales de agua y cañerías. Su uso debe ser racionalizado.
187
Tabla 35. Ficha técnica del Gel Coat
NOMBRE COMERCIAL Gel Coat
PRODUCTOR Autores
COMPONENTES Resina Poliéster, Peróxido, Talco Mineral y Pigmento.
Resina Poliéster: 50%
Talco Mineral: 40%
RELACION DE MEZCLA
Pigmento: 8%
Peróxido: 1% - 2%
El tiempo de endurecimiento depende de la proporción de acelerante
(Cobalto) que contenga la resina, así como la proporción de MEK
ENDURECIMIENTO
aplicado en la mezcla y la temperatura de trabajo. No debe utilizarse a
temperaturas inferiores a 15ºC.
Tanto el líquido como sus vapores son inflamables, nocivos si son
ingeridos, inhalados o absorbidos a través de la piel o mucosas. Produce
RIESGOS
irritación en la piel, ojos y tracto respiratorio. Afecta el sistema nervioso
central, hígado y sistema reproductivo.
No es reciclable, ni biodegradable. Se debe evitar su desecho en fuentes
IMPACTO AMBIENTAL
de naturales de agua y cañerías. Su uso debe ser racionalizado.
188
Figura 113. Catalizador MEK
190
7.1.3 Materiales empleados para la fabricación de piezas estructurales
Fuente: http://www.aireyespacio.com/2009/07/la-fibra-de-carbono.html
192
Terminado este proceso, se obtienen filamentos individuales de fibra de
carbono que son entrelazados en grupos de 5000 y 12000 filamentos,
conocidos con el nombre de Roving, los cuales se entrelazan formando
diferentes tipos de tejidos, siendo los más importantes los que se pueden
observar en la Figura 116.
Fuente: http://www.metalactual.com
a) Tejido Plano o Plain: En este tipo de tejido cada hilo longitudinal y pasa
por encima y debajo del hilo transversal, proporcionando buena
estabilidad, evitando su distorsión dimensional.
193
En la Tabla 37 se muestra la ficha técnica del tejido de fibra de carbono
seleccionado para la fabricación de las piezas estructurales del Quadrotor.
194
2. Resina Epóxica: La sustancia de impregnación que sirve como matriz
para la fabricación de las piezas compuestas es la resina epóxica,
polímero termoestable que se endurece cuando es mezclado con un
compuesto endurecedor. Este tipo de resina mostrada en la Figura 117A
fue seleccionada por que supera un dureza a otras como las de poliéster y
viniléster.
196
Tabla 39. Comparación entre las resinas epóxica y poliéster
CARACTERISTICAS RESINA EPOXICA RESINA POLIESTER
Precio Elevado Bajo
Facilidad de uso Fácil Difícil
Módulo de Tensión Excelente Muy bueno
Módulo de Flexión Excelente Muy bueno
Tolerancia a ultravioleta Bajo Muy bueno
Uso en compuesto de
Excelente No recomendable
Carbono o Kevlar
197
7.2.1 Empleo de software y equipo de maquinado
198
Figura 119. Máquina CNC Ledwel V-30
199
7.3 FABRICACIÓN DE CONTRA-MOLDES
200
En la Figura 121 se muestra el estado final del molde en madera de la sección
inferior del brazo una vez realizados los trabajos de acabado superficial.
Figura 121. Superficie del molde lijada y tratada con laca de poliuretano
201
Figura 122. Aplicación de cera desmoldante en la superficie del molde
202
3. Aplicación de Gelcoat: Se aplica manualmente con brocha una capa de
Gelcoat hasta recubrir totalmente la superficie del molde como se observa
en la Figura 124.
203
5. Aplicación de resina poliéster: Los trozos de fibra de vidrio se
impregnaron con la mezcla requerida de resina de poliéster y catalizador
con el fin de alcanzar la rigidez requerida una vez haya fraguado. En la
Figura 126 se detalla esta fase del proceso.
La Figura 128 muestra los contra-moldes terminados: (A) Sección inferior del
fuselaje, (B) Sección superior del brazo, (C) Sección inferior del brazo y (D)
Sección superior del fuselaje.
204
Figura 127. Moldes en madera y contra-moldes en fibra de vidrio
205
7.4 FABRICACIÓN DE PIEZAS ESTRUCTURALES
206
2. Aplicación de resina epóxica: Cuando se evaporó completamente el
desmoldante líquido, el contra-molde fue impregnado con una mezcla
de resina-endurecedor (1:1) con el fin de formar una película protectora
(Figura 130A), que se dejó secar por un tiempo de 2 horas.
207
Figura 131. Empleo de bolsa de vacío en el proceso de laminación
208
4. Polimerización con calor: Aunque el proceso de polimerización de
manera industrial se realiza en hornos autoclave, artesanalmente puede
exponerse la resina a condiciones ambiente (≥15ºC) obteniendo excelentes
resultados en cuanto a su resistencia mecánica. Para este trabajo se
incrementó la temperatura en las piezas estructurales por medio de una
pistola de aire caliente (Figura 133). Esta acción permitió una mayor
difusión de la resina a través de la fibra de carbono y una reducción en el
tiempo de fraguado.
209
7.4.2 Componentes estructurales terminados
210
1. Acondicionamiento de las piezas: Los excesos de material en las
secciones de la estructurales fueron retirados manualmente empleando
una herramienta de corte tipo motor-tool como se puede apreciar en la
Figura 136A. El estado final de este acondicionamiento de la pieza
correspondiente al brazo se observa en la Figura 136B.
211
7.5 ENSAMBLAJE DE PIEZAS Y COMPONENTES ESTRUCTURALES
212
Figura 139. Instalación de núcleos en el fuselaje
214
7.5.3 Ensamblaje de componentes estructurales
Como ya se mencionó, los brazos se ensamblan en los extremos del cuerpo
central mediante un sistema de broche mecánico tal como lo muestra la
Figura 143. Al introducir la sección rectangular de aluminio que sobresale del
brazo en el orificio cuadrado presente en el bloque de madera, el broche
superior e inferior también lo hacen hasta que llegan a su tope, quedando
asegurados e impidiendo el movimiento del brazo.
215
Figura 144. Conjunto superior e inferior de la estructura del Quadrotor
216
Figura 145. Instalación de los elementos del grupo moto-propulsor
217
Como se aprecia en la Figura 147 el ensamblaje de los brazos queda
terminado y listo para ser acoplado en los costados del fuselaje del
Quadrotor.
219
Figura 150. Montaje de las tarjetas electrónicas
220
Figura 152. Instalación de cúpulas en el cuerpo central del Quadrotor
221
8. PRUEBAS Y RESULTADOS
222
El resultado de este procedimiento se describe en la Tabla 40, donde se listan
estos pesos individuales.
Según las especificaciones técnicas del motor eléctrico Rimfire 0.10 este es
capaz de generar un empuje máximo 905 gr-f en condiciones de vuelo
vertical, por lo cual el grupo moto-propulsor del Quadrotor produce 3620 gr-f
de empuje máximo. De esta manera, el peso máximo en despegue del
Quadrotor representa un 73.2% de la capacidad máxima de empuje de los
motores, con lo que se cuenta un 26.8% de potencia disponible para ser
empleada en cualquier condición de vuelo.
224
Figura 155. Principales dimensiones del Quadrotor
225
8.1.2 Análisis de Esfuerzos
226
Cumplidas estas condiciones, se ejecutó el análisis estático en cada uno de los
componentes, obteniendo de esta manera los valores de deformación,
esfuerzo y factor de seguridad requeridos, los cuales son tabulados en la
Tabla 41.
227
Figura 156. Aplicación de carga en el extremo del brazo
228
Los resultados obtenidos en las anteriores gráficas permiten determinar qué:
229
8.2 PRUEBAS ELÉCTRICAS Y ELECTRÓNICAS
231
La realización de pruebas de integración en el banco vertical (Figura 161)
sirvió para hacer mover al Quadrotor en sus tres ejes de rotación (pitch, roll,
yaw), así como para afinar la sintonización de las ganancias del control PID y
verificar el correcto funcionamiento del grupo de sensores (Unidad inercial,
barómetro y magnetómetro).
232
Figura 162. Conexión Motor eléctrico – Control de velocidad
233
2. Control de Velocidad – Micro-controlador: Del control de velocidad sale
un cable que se conecta normalmente al receptor del radio-control, de
donde sale la señal de control PWM (Cable naranja), así como los cables
de alimentación (Cables rojo y marrón). Para la integración de los
componentes electrónicos en el Quadrotor se hizo necesario realizar una
modificación para llevar estas conexiones al micro-controlador.
234
3. Micro-controlador – Sensores - Receptor: Para enviar correctamente las
señales de control al ESC, el micro-controlador Arduino debe recibir
información del grupo de sensores (IMU, barómetro y magnetómetro) y
del receptor del equipo de radio-control. La Figura 165A muestra una
gráfica esquemática de las conexiones realizadas entre estos componentes.
235
Al accionar las palancas del radio-control (transmisor), se envía las señales
correspondientes (FM/PPM) al receptor, el cual las recibe y re-envía
gracias a las conexiones realizadas al micro-controlador, encargado de
procesarlas. Cada uno de los canales presentes en el receptor cumple una
función en el movimiento o control del Quadrotor, conexiones que se
explican en la Tabla 43.
236
Figura 166. Integración Batería – Alarma sonora de bajo voltaje
(89)
237
Tabla 44. Medición de la demanda de corriente según señal PPM
238
8.2.3 Desempeño de la batería sin carga útil
Se realizaron pruebas de vuelo estacionario sin carga, es decir, sin equipo útil,
observando el desempeño de la batería a un régimen constante de aceleración
con el fin de estimar la caída de tensión con respecto al tiempo de operación.
En la tabla 45 se listan los resultados de esta prueba.
239
Figura 167. Curva de la caída de voltaje de la batería, prueba sin carga útil.
En la Figura 167, se puede observar el comportamiento del voltaje (eje x), con
respecto al tiempo (eje y); donde se aprecia una caída de tensión casi lineal,
iniciando desde los 12.45 V (máxima carga por pack) hasta un voltaje mínimo
por motivos de seguridad de la batería de 9.6 voltios.
Según esta prueba que se realizó a un 50% de aceleración de los motores del
se estima un tiempo de vuelo del vehículo de 9 minutos.
Se realizaron pruebas de vuelo estacionario con una carga útil de 340 gramos
adosada en la parte inferior del fuselaje. La palanca de aceleración en el radio-
control se mantuvo a un 50% de su recorrido. En la Tabla 46 se listan los
resultados de esta prueba.
240
Tabla 46. Desempeño de la batería con carga útil
241
Figura 168. Curva de la caída de voltaje de la batería, prueba con carga útil.
242
1. Se incrementa la ganancia proporcional hasta que se tenga una respuesta
oscilante.
2. Se recomienda disminuir un poco la ganancia proporcional e incrementar
en pasos pequeños la ganancia integral mientras se evalúa el sistema en
cada caso.
3. Se aumenta la ganancia derivativa en pasos pequeños, hasta obtener un
comportamiento cíclico. Después de este ajuste se puede incrementar la
ganancia proporcional para mejorar la velocidad de respuesta.
4. Otro procedimiento es el siguiente: trabajar primero con una ganancia
proporcional que da una ligera oscilación ante la perturbación; luego se
aumenta la acción derivativa hasta eliminar la oscilación.
5. Finalmente se ajusta la acción integral con aumentos muy pequeños hasta
eliminar el error u offset.
243
1. Ángulo de alabeo con controlador proporcional
244
Prueba 2 Característica de la prueba
Proporcional (Kp) = 7
El sistema es más inestable en el
Integral (Ki) = 0
tiempo, presenta oscilaciones
Derivativo (Kd) = 0
muy violentas.
Tiempo de prueba = 10 s
245
Prueba 3 Característica de la prueba
Proporcional (Kp) = 4
El sistema presenta oscilaciones
Integral (Ki) = 0.010
en el tiempo, pero la respuesta
Derivativo (Kd) = 0
es más lenta.
Tiempo de prueba = 10 s
Figura 171. Respuesta del control proporcional Kp= 4, Ki= 0.010 Prueba 3.
246
Prueba 4 Característica de la prueba
Proporcional (Kp) = 4 El sistema se estabiliza, aunque tiene
Integral (Ki) = 0.020 pequeñas oscilaciones muy lentas
Derivativo (Kd) = 0 una tras otras. Comportamiento
Tiempo de prueba = 10 s sobreamortiguado.
247
Prueba 5 Característica de la prueba
Proporcional (Kp) = 4
El sistema describe un comportamiento
Integral (Ki) = 0.020
amortiguado, un tiempo de estabilización
Derivativo (Kd) = 10
muy corto.
Tiempo de prueba = 10 s
248
Prueba 6 Característica de la prueba
Proporcional (Kp) = 4
Integral (Ki) = 0.020 El sistema describe un comportamiento
Derivativo (Kd) = 20 amortiguado, mejorando la prueba 5.
Tiempo de prueba = 10 s
249
2. Combinación de los movimientos alabeo y cabeceo
250
Como resultado de esta última prueba, se puede concluir que estos valores
son la mejor combinación encontrada para optimizar la estabilidad del
Quadrotor. Sin embargo, este de hecho no sugiere que sea la única
combinación posible que cumpla con los requerimientos de diseño, ya que
queda para futuras tareas la búsqueda de una configuración más óptima del
controlador PID propuesto o simplemente lograrla con la aplicación de
métodos de control más complejos y eficaces.
251
8.4 PRUEBAS DE VUELO
252
Figura 177. Movimiento de alabeo hacia la derecha y hacia la izquierda
253
9. CONCLUSIONES
254
Se diseñó y programó un algoritmo de control basado en lenguaje de
programación C++ compatible con Arduino, con el fin de establecer las
rutinas de control y estabilidad, adaptándolas a los requerimientos del
diseño propuesto.
255
10. RECOMENDACIONES
256
Se debe tener un especial cuidado en la manipulación y operación de la
batería LiPo, ya que esta es susceptible a sufrir daños ocasionados en su
carga y descarga, que pueden incluso a producir incendio.
257
11. DESARROLLOS FUTUROS
258
BIBLIOGRAFIA
259
Bento María de Fátima, “UNMANNED AERIAL VEHICLES: AN
OVERVIEW”, insideGNSS, January-February, 2008.
260
T. Bresciani, MODELLING, IDENTIFICATION AND CONTROL OF A
QUADROTOR HELICOPTER, Department of Automatic Control 2008.
http://roboticsclub.org/redmine/attachments/467/Quadrotor_Bible.pdf
261
S. Bouabdallah, A. Noth, R. Siegwart, DYNAMIC MODELLING OF UAVS,
Technical Report Version 2.0, Autonomous Systems Laboratory - EPFL, May
2006.
http://www.sky-sailor.ethz.ch/docs/Dynamic_Modeling_of_Fixed-
Wing_UAVs_(12.05.2006).pdf
262
ANEXOS
263
ANEXO A. GEOMETRÍAS Y EL CÁLCULO DE LOS MOMENTOS DE
INERCIA
densidad constante de .
264
Cilindro sólido: Con radio , altura , masa y una densidad constante
265
266
Cilindro solido de densidad variable: Con radio , altura , masa y
Figura B.
267
Si el centro de masa del objeto no es igual al origen de los ejes (asumido como
el centro de gravedad del Quadrotor), existe un teorema (teorema de ejes
paralelos) el cual ayuda a determinar el momento de inercia sin necesidad de
recalcular las integrales.
El teorema establece que si el objeto es desplazado una distancia D desde el
eje de rotación del centro de gravedad, el desplazamiento y el momento de
inercia son relacionados con la siguiente ecuación:
268
Figura C. Rectángulo prismático macizo con centro de gravedad desplazado.
269
270
ANEXO B. PROTOCOLO DE COMUNICACIÓN I2C
(INTER-INTEGRATED CIRCUIT) ENTRE LOS SENSORES.
271
ANEXO C. ESTIMACIÓN DE LOS ÁNGULOS DE EULER
http://www.starlino.com/wp-content/uploads/data/imu_guide/01.png
272
Según este modelo de la figura A, el vector es el vector de fuerza que el
(1)
Para obtener las mediciones de estas componentes en los tres ejes, los
acelerómetros brindan información a través de un protocolo en serie como
, y , proporcionando la salida a un nivel de tensión dentro de
Considerando que el módulo ADC de 10 bits genera valores para los tres
canales o ejes del acelerómetro:
273
Cada módulo ADC tendrá un voltaje de referencia, en nuestro caso es de
3.3V. Para convertir un valor ADC de 10 bits a un voltaje se utiliza la
siguiente fórmula:
(2)
decimales)
274
Ahora se tienen las lecturas del acelerómetro en voltios, para convertir este
voltaje en gramos se aplica la sensibilidad que por lo general se expresa en
275
Figura B. Ángulos respecto a la horizontal del vector inercial.
Fuente: http://www.starlino.com/wp-content/uploads/data/imu_guide/02.png
Los ángulos que nos interesan son los ángulos entre del vector de
y similarmente:
276
Donde se deduce de la ecuación (1) que
(4)
(5)
277
Figura C. Modelo vectorial de las medidas del giróscopo
Fuente: http://www.starlino.com/wp-content/uploads/data/imu_guide/03.png
Cada canal del giróscopo mide la rotación alrededor de uno de los ejes. Por
ejemplo, un giróscopo de 2 ejes medirá la rotación alrededor de los ejes X e Y.
Para expresar esta rotación en primer lugar se define lo siguiente:
Pitágoras tenemos:
, y similarmente:
278
También se nota que:
(6)
Un giróscopo (a menos que sea un giróscopo digital especial) rara vez le dará
un valor expresado en . Igual que para el acelerómetro arrojará un
bits).
(7)
279
medir cada salida del eje con un voltímetro luego de soldarlas ya que la
mayoría sufre una ligera compensación)
Como un ejemplo, se supone que nuestro módulo ADC arroja los siguientes
valores:
En otras palabras, el dispositivo gira alrededor del eje Y (o podemos decir que
gira en el plano XZ) con una velocidad de 306 grados/s, y alrededor del eje X
con una velocidad de -94 grados/s.
280
Se observa que el signo negativo significa que el dispositivo gira en la
dirección opuesta a la dirección positiva convencional. Si la hoja de
especificaciones del giróscopo muestra en qué dirección es positivo, de lo
contrario se debe encontrar mediante la experimentación con el dispositivo y
tomando nota de la dirección de los resultados de la rotación de la tensión
creciente en el pin de salida.
281
El filtro complementario es un algoritmo que se inspiró en algunas ideas
utilizadas del filtro de Kalman, sin embargo, es mucho más simple y más fácil
de implementar en dispositivos embebidos. A partir de los valores como
o derivados de los datos del acelerómetro, se
282
Este valor será el resultado de la corrección de los valores sobre la base de
datos de los giróscopos y en base a datos estimados; chequeando siempre la
posición del vector con los datos anteriores de ; generando un
sucesivamente.
componentes:
283
Realizando el cálculo para una componente, en este caso . Observando
284
De la ecuación 1; se puede escribir la magnitud del vector como sigue:
Adoptando una notación temporal más corta para los siguientes cálculos:
interior de la ( ) por .
285
Note que , finalmente:
Nota: es posible simplificar aún más esta fórmula. Al dividir las dos partes de
la fracción por obteniendo:
Donde,
286
Ahora, se puede encontrar:
Donde; cuando , y
cuando
giróscopo.
287
Podemos simplificar esta fórmula, dividiendo el numerador y el
denominador de la fracción por w1.
288
Y luego se hallan los ángulos como se determinaron anteriormente en la
ecuación 4; según los triángulos que forman estas componentes.
289
ANEXOD. GUÍA DEL USUARIO (CONFIGURACIÓN DEL
QUADROTOR A TRAVÉS DE LA INTERFAZ GRÁFICA)
290
5. Indicador de los valores de los comandos del radio control.
6. Indicador de los valores de la señal PPM enviada a los motores.
7. Indicador de los valores que generan los sensores. Aquí se pueden desactivar
o activar en las casillas de los colores que identifican las curvas de los
sensores que aparecen en el recuadro 11.
8. Zona de configuración del controlador PID; para variar los valores de las
casillas del PID se mantiene oprimido y se desliza hacia la derecha para
aumentar el valor o hacia la izquierda para disminuirlo.
9. Indicadores de la activación o desactivación de los sensores. Estos cambian de
acuerdo al comando y configuración que se aplique en la casillas de los
auxiliares 1 y 2.
291
10. Escala de las curvas que grafican el comportamiento de los sensores. Se opera
igual que los botones deslizadores como en las casillas PID.
11. Zona de graficas de los sensores. En esta configuración la velocidad de
muestreo es de 100 ms, es decir, un punto de la gráfica demora 10 s en llegar
al otro extremo del recuadro.
12. Zona de instrumentos principales que indican el movimiento de alabeo,
cabeceo y rumbo del vehículo.
13. Imagen tridimensional que responde en tiempo real con los movimientos del
vehículo, simula el movimiento del quadrotor según la configuración que se
indique.
292
También parpadea cuando se transmiten datos y queda encendido cuando
está el Quadrotor listo para despegar, es decir, ARMADO.
5. Una vez visto lo anterior, oprima el botón INICIAR para observar el tipo de
configuración con que se cuenta, activación de sensores, movimiento del
vehículo en 3D, comandos principales, etc.
6. Antes de despegar el vehículo, es necesario calibrar los sensores; seleccione
CALIBRAR_ACELEROMETRO cuando tenga el vehículo en una zona
horizontal. Luego seleccione CALIBRAR_MAGNETOMETRO para ajustar el
sensor según el norte magnético de la tierra (se cuenta con 30 s para hacer
girar el Quadrotor y se corrija el error de medida).
7. Una vez realizada la calibración se puede continuar con la configuración del
vuelo.
8. Para comenzar a volar solo encienda el transmisor, encienda el interruptor
principal del Quadrotor y póngalo en una zona de despegue plana, calibre
sensores, coloque la palanca del acelerador en el mínimo, con el AUX1 se
arma el sistema; cuando esté encendido el led verde indicando el ARMADO,
ponga en máximo la palanca del acelerador para activar los ESC’s (suenan
dos zumbidos) y luego se pone en la posición mínima (sonará 4 veces). Esto
indica que el Quadrotor está listo para el despegue.
9. Por seguridad se tiene configurado la opción de desactivar los motores
cuando está en la posición mínima, una vez este armado el vehículo.
293
Configuración del comportamiento y modos de vuelo.
294
Con estas casillas se puede fácilmente activar o desactivar opciones a través
de dos interruptores en su radiocontrol. Los dos interruptores corresponden a
canal 5 (AUX1) y el canal 6 (AUX2).
Por medio de las palancas se puede activar funciones tales como armar o
desarmar el sistema activarlo/desactivarlo y calibrar sensores. Los
movimientos de las palancas se muestran a continuación:
295
Armar motores
Desarmar motores
Calibrar acelerómetro
Calibrar giróscopo
Calibrar magnetómetro
296
ANEXO E. ESPECTROS DE LOS ANALISIS MECANICOS EN ANSYS.
297
Figura E4. Aplicación de momento en el núcleo del cuerpo central
298
Figura E7. Aplicación de momentos en el cuerpo central
299
Figura E10. Aplicación de carga en la bancada de la batería
300
Figura E13. Determinación del factor de seguridad en el brazo
301
ANEXO F. ESQUEMA DE INTEGRACION DE LOS EQUIPOS
ELECTRICOS Y ELECTRONICOS.
302
303