PROYECTO DE Robotica
PROYECTO DE Robotica
PROYECTO DE Robotica
Facultad de Ingeniería
Programa de Ingeniería Mecatrónica
ROBÓTICA
CICLO: 2021 - I
Trujillo, Perú
2021
UNIVERSIDAD NACIONAL DE
TRUJILLO Resumen
El presente informe del proyecto final del curso de Robótica, consiste en el análisis del
modelo cinemático, estático, dinámico y planificación de trayectorias del robot de 6 grados de
libertad EPSON C12XL el cual utilizo un efector final tipo ventosas para formar cajas de
cartón.
Primeramente, se realizó una revisión bibliográfica de trabajos de investigación científica en
revisas, publicaciones de tesis avocadas al análisis de la cinemática y dinámica del robot. En
esta revisión se encontró muchos trabajos dirigidos a robot de pocos grados de libertad, por lo
que posteriormente se planteó en la etapa de metodología desarrollar el análisis del modelo
cinemático, dinámico y planificación de trayectorias del robot EPSON C12XL de 6 grados de
libertad.
En la etapa de metodología se utilizó el método de desacoplo cinemático para la cinemática
del robot, teniendo en cuenta el método Denavit H. se analizo la cinemática diferencial, así
como el análisis estático, el análisis dinámico se utilizó el método de Lagrange-Euler.
Finalmente se plasmó en un algoritmo el cual se logró hacer una conexión entre los softwares
de CoppeliaSim, Matlab y SolidWorks.
INGENIERÍA
MECATRÓNICA
Abstract
In this project, the design of a.
Capítulo 1 Introducción..............................................................................................................7
Realidad Problemática.......................................................................................................7
Objetivos............................................................................................................................8
Capítulo 4 Metodología.......................................................................................................36
Capítulo 6 Conclusiones......................................................................................................38
Capítulo 7 Recomendaciones..............................................................................................39
Referencias Bibliográficas...................................................................................................40
Anexos.................................................................................................................................41
Lista de tablas
Lista de figuras
UNIVERSIDAD NACIONAL DE TRUJILLO
Capítulo 1
Introducción
El presente trabajo de investigación formativa aborda el análisis del modelado cinemático,
estático, dinámico y planificación de trayectorias del robot EPSON C12XL de 6 grados de
libertad. Para el estudio de la cinemática directa e inversa, cinemática diferencial y para el
comportamiento dinámico del robot, se usarón los softwares de CoppeliaSim y Matlab.
La característica principal de este tipo de robot paletizador es coger un objeto de un
determinado lugar y llevarlo a otro para su organización en una estiba (pick and place). Para
ello es necesario que se establezca una sincronización espacio-temporal entre el brazo y el
objeto en movimiento, para que este último sea trasladado.
El área de trabajo es un factor que influye en el movimiento de este robot, por eso también se
debe tener en cuenta el posicionamiento del brazo robótico para que pueda desplazar los
objetos sin que haya obstáculos a su alrededor y pueda cumplir con el tiempo requerido de
sincronización con el objeto en movimiento.
La presente investigación formativa que tiene como finalidad el análisis del modelo
cinemático, estático, dinámico y planificación de trayectorias del robot EPSON C12XL para
formar cajas de cartón con ventosas. se realizó con el interés de poder aplicar los
conocimientos adquiridos en clase para poder dar una solución que permita reducir el tiempo
de ciclo del robot. En el marco de la teoría de control, el proceso que se realizo fue primero el
de obtener la morfología, localización espacial, cinemática directa e inversa y la cinemática
diferencial de robot.
Por último, se analiza el comportamiento dinámico del robot para evaluar una mejor
planificación de trayectoria en el software de CoppeliaSim.
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO
Realidad Problemática
El BantamPro ELSP de la firma Brenton (Pro Mach) es una empresa que brinda
soluciones compactas de empacado en cajas, de bajo costo, con funciones de máquinas
previamente separadas. Este sistema acomoda cajas RSC y automáticamente realiza múltiples
operaciones de armado de las cajas, carga, sellado y paletizado. También cuenta con la
integración de los equipos adicionales, tales como: etiquetadora, apilador de palés
semiautomático y envoltura encogible. Utilizan un brazo robótico FANUC M-710iC el cual
puede cargar cajas y también paletizar. Este sistema es compacto y a la vez flexible.
INGENIERÍA MECATRÓNICA
UNIVERSIDAD NACIONAL DE TRUJILLO
Por todo lo expuesto anteriormente, podemos darnos cuenta que la falta de inserción
de robot en la industrial genera pérdidas en las empresas ya sea porque los trabajos manuales
demandan mucho tiempo, así como el agotamiento de los trabajadores.
Objetivos
Objetivo General
análisis del modelo cinemático, estático, dinámico y planificación de trayectorias del robot
Objetivos específicos
Planificación de trayectorias
INGENIERÍA MECATRÓNICA
Capítulo 2
Antecedentes (Estado del Arte)
ALONZO y otros (2014). En quito, Ecuador; presenta la tesis sobre “Diseño,
construcción y control de un brazo robótico con cuatro grados de libertad”. En este trabajo
usan conceptos y teorías sobre diseño y hacen posible la simulación del funcionamiento del
estructura. Nos muestran la secuencia detallada de la construcción y montaje del brazo, con
trabajo utilizan como controlador una tarjeta electrónica de control programable Arduino
Mega esta misma permite comandos inalámbricos utilizando su módulo Wii, además utiliza
también un control sencillo con sus módulos tipo puente H para controlar los movimientos.
Programa también una interfaz gráfica para visualizar los movimientos en pantalla y realizar
el control de los movimientos vía Web de 2 tipos, emplea un circuito para controlar
Para culminar utiliza modelos cinemáticos y modelos geométricos para calcular la cinemática
inversa que será necesaria para traducir las coordenadas geométricas de cada motor en
simular el control de un Brazo Robótico Industrial automatizado. Por ese motivo y para lograr
esa finalidad cuenta con un brazo Robot VELLEMAN KSR10 de 5 GDL que en su forma
2560. Dicho Robot cuenta con 5 articulaciones que a su vez son los grados de libertad(GDL)
y que por lo consiguiente el brazo cuenta con 5 motores que cada motor controlara a
determinada articulación, este Brazo no cuenta con control de posición según el trabajo a
realizar, por ello a base de los estudios realizados se pretende modificarlo para que pueda
la carrera e investigando otros no conocidos con el fin de ser un buen profesional al servicio
de la Industria en General”.
Dinámica de un brazo Robótico de 4 grados de libertad para tareas sobre un plano Horizontal
“ nos dice, esta investigación solo es un segmento de todo un proyecto para lograr el
una mejora de los procesos 20 aquí en el Perú, para ello debemos utilizar los conocimientos
adquiridos y poder llevar acabo nuestras ideas que nos permitan mediante herramientas lograr
un plano”. El plano a que se refiere es el A3, este diseño realizado en forma parcial por Luis
Pontificia Universidad Católica del Perú, este trabajo fue concluido y rediseñado por el Sr.
Alberto Orihuela técnico Egresado de SENATI el que fue participe en la fabricación del
brazo robótico que se llevó a cabo en su taller. Se dividió en tres etapas la realización del
diseño de todo el sistema que controle el brazo robot, la primera etapa de interfaz de usuario
que haga posible que el mismo pueda manipular en conjunto el sistema con la ayuda de un
procesador individual y hacer posible compartir información con el controlador, este por su
parte se encargara de recibir toda la información que definan todos los valores de la
trayectoria del procesador personal que sean ingresados por el operador para producir las
señales que permitirán el control requerido en el accionar de los actuadores que harán posible
el movimiento, de la misma manera recibe información de todos los sensores que están
acoplados al sistema, como etapa final se tiene la interfaz de potencia que se encargara de
recibir las señales de control producidas por el micro controlador , que permitirán a los
que nos propone realizar los cálculos del diseño mecánico y electrónico requeridos para que
la maquina funciones correctamente. Para que esto se logre obtienen el modelo cinemático
método geométrico. Logran el modelo dinámico del robot solucionando las ecuaciones de
elementos mecánicos, ensamblaje y planos mecánicos del robot, exportan todos los datos
CAD del robot al software Matlab con la finalidad de comprobar lo que se propone en el
diseño. Utilizan el software Electrónico Eagle (atreves de este programa podemos ver si los
circuitos funcionan correctamente) para poder llevar a cabo el esquema de todos los circuitos.
Para obtener los elementos electrónicos para el proyecto se buscó información de diferentes
Sobre los resultados de diseño de este brazo robótico fueron favorables pues se tuvo
en cuenta la seguridad para su uso ya que cuenta con sensores de corriente esto evitaría
posibles sobre cargas en los motores, y además se le instalo un sistema para poder parar el
funcionamiento del robot en caso que sea necesario. Cabe señalar que realizando diferentes
los movimientos para poder utilizar el brazo robótico para diferentes tareas o trabajos que se
requieran realizar. “Este proyecto tubo fines didácticos para su uso educacional y también de
investigación.
Capítulo 3
Marco Teórico
Estructura Mecánica:
Mecánicamente, un robot está formado por una serie de elementos o eslabones unidos
mediante articulaciones que permitan movimiento relatico entre cada dos eslabones
una combinación de ambos. De este modo es posible los 6 tipos de diferentes articulaciones
Barrientos ]
“estructura paralela”, son menos frecuentes, a pesar que en los últimos años, algunos
paralela conocida como plataforma de Steward Gougth y un robot industrial comercial con
estructura paralela.
estructura Cartesiana y SCARA. Los robots de estructura esférica y cilíndrica, más frecuentes
Figura 14-3.- Estructuras paralelas: a) Plataforma de Stewart-Gouth. b) Robot industrial paralelo [Fuente:
Transmisiones y Reductores:
Las transmisiones son los elementos encargados de transmitir el movimiento desde los
actuadores hasta las articulaciones. Los reductores, junto con las transmisiones, son
encargados de adaptar el par y la velocidad de la salida del actuador a los valores adecuados
para el movimiento de los elementos del robot.
Transmisiones
Un robot, al mover su efector final, lo hace con aceleraciones elevadas, por lo que es de
gran importancia reducir al máximo los momentos de inercia que este posee. Del mismo
modo los pares estáticos que deben vencer los actuadores dependen directamente de la
distancia de las masas al actuador. Por este motivo es que se procura que los actuadores que
son generalmente pesados, estén lo más cerca posible de la base del robot. Esta circunstancia
obliga a utilizar sistemas de transmisión que trasladen el movimiento hasta las articulaciones,
especialmente las que se localizan en el extremo del robot. También son utilizadas para
convertir movimiento circular en lineal o viceversa, lo que en ocasiones puede ser necesario.
Actualmente, en el mercado de robots industriales, existen con acoplamiento directo entre
accionamiento y articulación ventajosos en numerosas ocasiones. Se trata, sin embargo, de
casos particulares dentro de la generalidad que, en los robots industriales actuales, supone la
existencia de sistemas de transmisión junto con reductores para el acoplamiento entre
actuadores y articulaciones.
1.2. Sistema Robótico
Los sistemas robóticos se crean con un propósito y es el de aumentar un volumen de
producción de una compañía manufacturera, objetivo de investigación en este documento.
Muchas cuentan con maquinaria y personal humano que se encargan de realizar tareas
manuales, con un alto costo para la empresa en tiempo y dinero.
Los sistemas basados en robótica efectúan tareas complejas como el embalaje, paletizado o
la preparación de envío con mucha rapidez y eficacia de manera totalmente automatizada
durante las 24 horas del día, (ver Figura 15).
paletizado-automatico-y-robot-despaletizador-de-cajas/]
1.3. Cinemática Del Robot Manipulador
La cinemática se encarga de estudiar el movimiento del robot con respecto al sistema de
robot sin considerar las fuerzas o momentos que originan el movimiento. La cinemática de un
x=f x (∅1 , ∅ 2 , ∅ 3 , … , ∅n )
y=f y ( ∅1 , ∅ 2 , ∅ 3 , … , ∅n ) (1)
z=f z (∅ 1 , ∅ 2 , ∅3 , … , ∅ n)
Rotacion Traslacion R P3 x1
T= [ Perspectiva Escalado ][
= 3x3
f 1 x 3 w 1 x1 ] (2)
Las dos primeras matrices son las más importantes en el tema de la robótica, la matriz de
perspectiva es nula y la submatriz de escalado global es una constante de valor uno.
R P3 x1
T = Rotacion Traslacion = 3 x 3
[ 0 Escalado 0 ][
1 ] (3)
Por lo general un robot de n grados de libertad está formado por n eslabones unidos por n
articulaciones, de tal modo que cada par eslabón-articulación constituye un grado de libertad.
A cada eslabón se le asocia un sistema de referencia solidario a él y, utilizando matrices de
transformación homogénea, es posible representar las rotaciones y traslaciones relativas entre
los distintos eslabones que componen el robot.
La matriz de transformación homogénea representa la posición y orientación relativa entre
[
T 1= sen (θ 1) cos (θ1)
0
0
0
0
0 L 1 sen (θ1 )
1
0
0
1
] (4)
La ventaja que ofrecen las matrices de transformación homogénea es su composición, es
T 2=0¿ 0T 1∗1T 2 ¿
Para hallar los parámetros de Denavit-Hartenberg se hace uso de los siguientes pasos:
Definir el origen O ien la intersección z icon la normal común entre los ejes z i−1y z i
Escoger el eje x i a lo largo de la normal común a los ejes z i−1y z i con dirección de la
articulación i+1.
Escoger el eje Yi, para que con Xi se cumpla la regla de la mano derecha.
Escoger el sistema de coordenadas del efector final donde x n sea normal a Zn-1
θi=¿ Ángulo que forman los ejes x i−1 y x i medido en un plano perpendicular al eje z i−1,
giratorias.
d i=¿Distancia a lo largo del eje z i−1 desde el origen del sistema de coordenadas
(i–1)-ésimo hasta la intersección del eje z i−1 con el eje xi. Se trata de un parámetro variable
en articulaciones prismáticas.
a i=¿ Distancia a lo largo del eje xi que va desde la intersección del eje z i−1 con el eje xi
de articulaciones prismáticas, se calcula como la distancia más corta entre los ejes z i−1
y z i.
α i=¿ Ángulo de separación del eje z i−1 y el eje z i, medido en un plano perpendicular al
El problema cinemático inverso tiene como objetivo encontrar los valores de las
Para resolver el problema cinemático inverso es mejor encontrar una solución cerrada. Si
q k =f k ( x , y , z , φ , θ , ψ) (5)
k =1 … n(GDL)
consideren sólo los primeros grados de libertad, que posicionan el efector final.
Los datos de partida son las coordenadas ( p x , p y , p z) referidas al sistema base. Este robot
posee una estructura planar, quedando este plano definido por el ángulo de la primera
q 1=arctg ( p y / p x ) (6)
Considerando sólo los eslabones 2 y 3, utilizando el teorema del coseno:
r 2= p 2x + p 2y
± √ 1−cos 2 q 3
sen q3 = (8)
cos q3
Como se ve, existen 2 posibles soluciones para q 3, según se tome el signo positivo o
pz pz
β=arctg ( )
r
=arctg
(
± √ p2x + p2y )
(9)
l 3 sen q 3
α =arctg
( l +l 23 cos q 3
2
2
)
Finalmente.
pz l 3 sen q3
q 2=arctg
(√
± p +p 2
x
2
y
) −arctg
( l +l 23 cos q3
2
2
) (10)
T = n o a p =[t ij ]
[
0 0 0 1 ] (11)
T relaciona el efector final con la base del robot.
¿¿
(12)
¿¿
Dado que ❑0T 3 es conocida se pueden realizar las operaciones descritas en la ecuación (10)
e ir comparando los términos, con el fin de poder obtener el valor de las coordenadas
necesarios 6 GDL, por lo general los ejes de las 3 últimas articulaciones se cortan en un
2) Obtener las matrices de Transformación Homogénea que relacione el eslabón i-ésimo con
3) Se tiene un vector que relaciona la base del robot y el punto que intercepta los 3 últimos
sistemas coordenadas.
5) La orientación del extremo del robot 0 R6 es conocida, además0 R6 =0R 3 3 R6 , del cálculo
3
R6=( 0 R3 )−1 0 R 6 (13)
Además:
3
R 6 = 3R 4 4 R 5 5 R 6 (14)
Igualando términos de ambos lados de la ecuación 14, se puede obtener los valores de
q 4 , q5 y q6. Con esto queda resuelto el problema cinemático inverso.
(Espacio articular) y de las velocidades de posición y orientación del efector final (Espacio
X =f x (q 1 , q 2 , q3 , … , q n−1 , qn )
Y =f y ( q 1 , q 2 , q 3 , … , q n−1 , q n)
Z=f z ( q 1 , q 2 , q 3 , … , q n−1 , q n)
(15)
α =f α ( q1 , q2 , q3 , … , qn−1 ,q n )
β=f β ( q1 , q2 , q3 , … , qn−1 , qn )
θ=f θ ( q 1 , q 2 , q 3 , … , q n−1 , q n )
n
∂ f z ∂ qi
Ż= ∑( ∂ qi ∂ t )
i=n
(16)
n
∂ f α ∂ qi
α̇ = ∑( ∂ qi ∂ t )
i=n
n
∂ f β ∂ qi
β̇= ∑( ∂ qi ∂ t )
i=n
n
∂ f θ ∂q i
θ̇= ∑( ∂ qi ∂t )
i=n
Expresado en forma matricial queda de la siguiente forma.
∂f x ∂f x ∂f x
…
[ ]
∂ q1 ∂ q2 ∂ qn
∂f y ∂f y ∂f y
…
Ẋ ∂ q1 ∂ q2 ∂ q2 q̇1
[] []
Ẏ ∂fz ∂fz ∂fz q̇ 2
…
Ż ∂ q1 ∂ q2 ∂ q2 q̇3
α̇
= (17)
∂fα ∂fα ∂fα q̇4
β̇ …
∂ q1 ∂ q2 ∂ q2 q̇5
θ̇ ∂f β ∂f β ∂f β q̇6
…
∂ q1 ∂ q2 ∂ q1
∂f θ ∂fθ ∂fθ
…
∂ q1 ∂ q2 ∂ q1
J(q)
Por lo tanto:
q̇1
0 ❑
❑V n
[ ]
0 ❑
❑ω n
6x1
=J ( q )
[]
:
:
.
q̇n
6xn nx1
(18)
q˙1
0
❑ Ẃ n =J ω ( q )
[]
:
:
.
q˙n
(20)
0
❑ Ẃ n=❑0Ẃ 1 + ❑0Ẃ 2 + ❑0Ẃ 3 +…+ ❑0Ẃ n (21)
0
i−1
❑
[]
Ẃ i= q̇i . ḱ , ḱ = 0
1
(22)
0 0 0 1 0 2 0 n−1
❑Ẃ n=❑Ẃ 1 + ❑R1 ❑Ẃ 2+ ❑ R2 ❑Ẃ 3 +…+ ❑ Rn −1 ❑ Ẃ n
(23)
0
❑Ẃ n=q˙1 . ḱ + ❑0 R1 q̇2 . ḱ + ❑0 R2 q̇3 . ḱ +…+ ❑0 Rn−1 q̇n . ḱ
0
❑Ẃ n=ε i q˙1 Ź 0+ ε i q̇2 ❑0 Ź 1+ ε i q̇3 ❑0Ź 2 +…+ ε i q̇n ❑0Ź n−1
(24)
ε i ε i=0 ; Prismática
{ ε i=1 ; Revoluta
n
0
❑Ẃ n=∑ ε i q̇ i ❑0 Ź i−1 (25)
i=1
q̇1
3x1
0
❑
0 0
J ω ( q3xn
)
0
Ẃ n=[ ε 1 Ź 0 ε 2 ❑Ź 1 ε 3 ❑Ź 2 … ε n ❑ Ź n−1 ]
[]
q̇2
:
q̇n
nx1
(26)
R . R T =I
Derivando en el tiempo:
Ṙ . R T + R . Ṙ T =0
Ṙ . R T =−R . ṘT
Ṙ . R T =−( Ṙ . R T )T
Ω=Ṙ . RT
Capítulo 4
Metodología
Cinemática Directa:
# De eslabones θ d a
1 q1 4724 100 −π / 2
2 q 2−π /2 0 650 0
3 q3 0 30 −π / 2
4 q4 647.5 0 π /2
5 q5 0 0 −π / 2
6 q6 84 0 0
Cinemática Inversa:
Pm =[ x , y , z ]
y
tan ( q1 )=
x
q 1=arctan ( xy )
Sabemos que:
z ' 0=z−L1
r =√ x 2 + y 3
2
R= ( r −t ) + ( z −L1) 2
√
Para q 2:
π
q 2= −γ −β
2
Donde:
z 0−L1
tan ( β )=
r −t
β=arctan ( zr−L
−t )
0 1
y
L23=R2 + L22−2 R L2 cos ( γ )
R2 + L22−L23
cos ( γ ) =
2 R L2
sen ( γ )
γ =arctan
( cos ( γ ) )
Para q 3:
π
q 3= −η
2
Como 0 R6 será la orientación deseada del efector final del robot, se tiene:
0
R6 =0R 3∗3 R6
Donde:
3
R6= 3R 4∗4 R5∗5 R 6
Por lo que:
3
R6= 0R 3T∗0R 6
3
R6= 3R 0∗0 R6
cos (q 4 )∗cos (q 5)∗cos(q6 )−sin( q4 )∗sin(q 6) −cos(q6 )∗sin (q 4)−cos(q4 )∗cos( q5 )∗sin( q6 ) −cos (q 4 )∗
3
❑
[
R6 = cos (q 4 )∗sin (q6 )+cos (q5 )∗cos(q 6 )∗sin(q4 ) cos(q4 )∗cos( q6 )−cos (q5 )∗sin(q 4 )∗sin( q6 ) −sin(q 4 )∗s
Como ya se obtuvo
cos( q6 )∗sin( q5 )
0
−sin (q5 )∗sin(q 6)
R3 y ❑0 R6 son datos numéricos, se obtendría la matriz numérica:
cos(q
❑
r 11 r 12 r 13
[
r ij = r 21 r 22 r 23
r 31 r 32 r 33 ]
Y comparando con la matriz ❑3 R6 , se escogen los siguientes términos para poder obtener
los valores de q 4 , q5 y q 6
r 23=−sin(q 4 )∗sin(q5 )
r 31=cos (q 6)∗sin(q 5)
r 33=cos (q 5)
De donde obtenemos:
r 23
q 4=arctan ( )
r 13
q 5=arccos ( r 33)
−r 32
q 6=arctan ( )r 31
El cálculo de la estática de este robot, emplea las siguientes fórmulas por recursividad:
Eslabón 3
2
(Fuerza
F2,3 =2 R3∗3 F 3,4−m3∗2R 0∗0 ǵ resultante
ejercida sobre el
eslabón 3 por el
eslabón 2)
Eslabón 2
1
(Fuerza
F1,2 =1R 2∗2 F2,3 −m2∗1 R0∗0ǵ resultante
ejercida sobre el
eslabón 2 por el
eslabón 1)
Eslabón 1
0 (Fuerza
F 0,1=0 R1∗1 F 1,2−m1∗0 R 0∗0 ǵ resultante
ejercida sobre el
eslabón 1 por el
eslabón base)
Para los Momentos:
6 '
M 6,7= [ M x M y M z ]
entonces las posiciones de los centros de masa respecto del sistema base son:
0.03
Ṕ1= 0
[ ]
0.32
0.08
Ṕ2= 0
[ ]
0.78
0.09
Ṕ3= 0
1.14 [ ]
0.46
Ṕ4 = 0
1.15[ ]
0.76
Ṕ5=
[ 0
1.15 ]
0.81
Ṕ6=
[ 0
1.15 ]
Hallamos la matriz de inercia D ( q ) :
3
D ( q ) =∑ ( mi∗J iv T∗J iv + J iωT 0R i∗~
I i∗0 RTi ∗J iω )
i=1
∂ Ṕ1 ∂ Ṕ 1 ∂ Ṕ1
J 1v =
[ ∂ q1 ∂ q2 ∂ q3 ]
0 0 00 0 0
[
J 1v = 0 0 0 0 0 0
0 0 00 0 0 ]
Jacobiano de velocidad angular para el eslabón 1
0 0 00 0 0
[
J 1ω= 0 0 0 0 0 0
1 0 00 0 0 ]
Jacobiano de velocidad lineal para el eslabón 2
∂ Ṕ2 ∂ Ṕ 2 ∂ Ṕ2
J 2v =
[ ∂ q1 ∂ q2 ∂ q3 ]
Jacobiano de velocidad angular para el eslabón 2
0 sen q 1 0 0 0 0
1 [
J 2ω= 0 −cos q 1 0 0 0 0
0 00 0 0 ]
Jacobiano de velocidad lineal para el eslabón 3
0 sen q 1 sin q 1 0 0 0
3
ω
1[
J = 0 −cos q 1 −cos q 1 0 0 0
0 0 0 0 0 ]
Matrices inerciales del centro de masa de cada eslabón
Ix x Ix y Ixz
0
❑ 1
[
I = I yx
I zx
1
1
Iyy
I yz
1
1
1
I yz
Iz z
1
1
]
I xx Ix y Ixz
0
❑ 2
[
I = I yx
I zx
2
2
Iyy
I yz
2
2
2
I yz
Iz z
2
2
]
I xx Ix y Ixz
0
❑ 3
[
I = I yx
I zx
3
3
Iyy
I yz
3
3
3
Iyz
Iz z
3
3
]
Reemplazamos en la matriz de inercia
~ 0 T 1
D ( q ) =m1∗J 1T 1 1T 0
v ∗J v + J ω ❑ R1∗ I 1∗❑ R1 ∗J ω +¿
~ 0 T 2
m 2∗J 2T 2 2T 0
v ∗J v + J ω ❑ R 2∗ I 2∗❑ R 2 ∗J ω + ¿
~ 0 T 3
m3∗J 3T 3 3T 0
v ∗J v + J ω ❑ R 3∗ I 3∗❑ R3 ∗J ω
3 3
C ( q , q̇ )=∑ (∑ (
j=1 k=1
)
C ijk ∗q̇ k ) ∗q̇ j
c 11 c 12 c 13
[
C (q , q̇)= c 21 c 22 c 23
c 31 c 32 c 33 ]
*Los términos reales no se colocan debido a que son demasiado extensos.
3
G ( q ) =−∑ ( J kv T ∗mk ∗ǵ0 )
i
k=1
[]
ǵ0= 0
−g
[]
g2
g
G(q)= 3
g4
g5
g6
Finalmente, hallados todas las matrices necesarias, la dinámica del robot de 3GDL
sería:
τ =D ( q ) q̈+C ( q , q̇ ) +G(q)
Capítulo 5
Resultados y Discusiones
Barrientos, Antonio., Peñín, Luis Felipe., Balaguer, Carlos., Aracil, Rafael., 2007.
Fundamentos de Robótica, 2 nd ed. Mc Graw Hill.
Craig, John J., 2006. Robótica, 3rd ed. Pearson
Siciliano, Bruno., Sciavicco, Lorenzo., Villani, Luigi., Oriolo, Giuseppe., 2009. Robotics:
Modelling, Planning and Control. Springer.
SOTO (2015). Un brazo Robótico de 5GDL con sistema de control modificable por el