TFM Helen Andreina de Lima Briceño
TFM Helen Andreina de Lima Briceño
TFM Helen Andreina de Lima Briceño
Autora:
Helen Andreina De Lima Briceño
Tutora:
Dra. Cecilia Elisabet García Cena
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Agradecimientos
Especialmente quiero dar las gracias a mi familia y amigos por estar siempre presentes,
dándome apoyo y fuerzas para alcanzar mis metas, sin ellos no hubiera sido posible culminar
mis estudios de maestría.
También agradezco a mi tutora Cecilia García Cena, por el apoyo, consejos y dedicación, así
como también a Marie André Destarac Eguizabal y a todo el equipo de trabajo de la empresa
Aura Innovative Robotics por brindarme su ayuda y disposición en todo momento.
Y por último, y no menos importante, deseo agradecer al profesor Antonio Barrientos, por
dedicar su tiempo, disposición, conocimientos y experiencia en ayudar con cualquier duda
referente al máster.
ii
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Resumen
Este trabajo de fin de maestría se basa en el análisis cinemático de un nuevo exoesqueleto para
rehabilitación de miembro superior, desarrollado en Aura Innovative Robotics, ésta órtesis es
ajustable a diferentes estaturas, presenta materiales ligeros y cuenta con 6 servomotores que
proporcionan un sistema de 6 grados de libertad. Dentro de este proyecto, se evalúa el espacio
de trabajo del exoesqueleto a partir del estudio de la manipulabilidad del mismo.
Para ello, se ha realizado una investigación sobre los tipos de lesiones del plexo braquial, los
posibles movimientos articulares del miembro superior, el estudio de diferentes órtesis para
rehabilitación existentes en el mercado y las principales características del nuevo prototipo a
utilizar, como base para calcular el nuevo desempeño cinemático.
Se hizo el análisis cinemático directo por el método de Denavit – Hartenberg y por Matrices de
Desplazamiento, la cinemática inversa se efectuó por el procedimiento de desacople cinemático
y se desarrolló una interfaz en MATLAB, para obtener gráficamente estos resultados en
conjunto con el índice de Yoshikawa, y proporcionar el espacio de trabajo de manera visual,
que delimita el rango de movimiento del robot.
Debido a que el nuevo prototipo es ajustable, se realizaron una serie de simulaciones para
apreciar la manipulabilidad del exoesqueleto de miembro superior, variando las estaturas de los
individuos, con el fin de comparar como influye el tamaño del paciente al momento de
determinar el espacio de trabajo.
Por último, se compararon los rangos de movimiento fijados al nuevo exoesqueleto con respecto
a órtesis comerciales y se evaluaron las ventajas e inconvenientes, con el objetivo de elaborar
propuestas de mejoras para trabajos futuros.
iii
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Índice general
Capítulo 1 ............................................................................................................................................... 1
Introducción ....................................................................................................................................... 1
1.1. Objetivos ............................................................................................................................. 2
Capítulo 2 ............................................................................................................................................... 3
Lesión del plexo braquial .................................................................................................................. 3
1.1. Etiología .............................................................................................................................. 4
1.2. Epidemiología ..................................................................................................................... 5
1.3. Clasificación de la lesión ..................................................................................................... 6
1.4. Cuadro clínico ................................................................................................................... 12
1.5. Tratamientos ...................................................................................................................... 14
Capítulo 3 ............................................................................................................................................. 17
Estado del arte ................................................................................................................................. 17
3.1. Prototipos de exoesqueletos de miembro superior: ........................................................... 17
3.2. Exoesqueletos comerciales de rehabilitación para miembro superior: .............................. 28
Capítulo 4 ............................................................................................................................................. 35
Descripción del Exoesqueleto de Rehabilitación a utilizar .......................................................... 35
4.1. Diseño del Exoesqueleto de rehabilitación de miembro superior ..................................... 35
4.2. Hardware del Exoesqueleto de rehabilitación de miembro superior ................................. 36
4.3. Software del Exoesqueleto de rehabilitación de miembro superior .................................. 37
Capítulo 5 ............................................................................................................................................. 38
Cinemática del Exoesqueleto de rehabilitación de miembro superior ....................................... 38
5.1. Cinemática Directa del Exoesqueleto de rehabilitación de miembro superior .................. 38
5.2. Cinemática Inversa del Exoesqueleto de rehabilitación de miembro superior .................. 43
5.3. Manipulabilidad del sistema .............................................................................................. 45
5.4. Espacio de trabajo ............................................................................................................. 46
5.5. Plataforma software para el cálculo cinemático: ............................................................... 47
Capítulo 6 ............................................................................................................................................. 50
Análisis de Resultados ..................................................................................................................... 50
Capítulo 7 ............................................................................................................................................. 64
Conclusiones y Trabajos Futuros .................................................................................................. 64
Bibliografía .......................................................................................................................................... 66
iv
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Anexos .................................................................................................................................................. 70
Anexo A: código Matlab ............................................................................................................... 71
Anexo B: Diagrama de Gantt ........................................................................................................ 83
v
Capítulo 1
Introducción
Recientemente es muy común hablar de exoesqueletos del miembro superior diseñados para
ayudar al hombre en situaciones de manipulación y rehabilitación, pero cada día el alcance de
estos dispositivos aumenta y ofrece sistemas más confiables. Es por ello, que los grandes
avances tecnológicos han permitido un auge en el diseño y desarrollo de dichos exoesqueletos.
Por otro lado, el mejorar la calidad de vida de personas con ciertas discapacidades es también
un punto de partida que ha motivado el presente trabajo. El cual, parte del estudio de las lesiones
de plexo braquial, como se clasifican y cuáles son los tratamientos para luego evaluar el
funcionamiento de diversos exoesqueletos del miembro superior (para adultos) existentes en el
mercado, que faciliten y optimicen el proceso de rehabilitación tanto para los pacientes con
lesiones del plexo braquial como para el terapeuta o especialista.
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
1.1. Objetivos
Objetivo General:
Objetivos específicos:
1. Localizar e identificar el tipo de lesión del plexo braquial con el fin de evaluar los
movimientos de rehabilitación adecuados mediante un exoesqueleto.
2. Comprender y evaluar los aportes de diferentes exoesqueletos para rehabilitación de
miembro superior existentes en el mercado.
3. Realizar el desempeño cinemático directo e inverso de un exoesqueleto para
rehabilitación de miembro superior en adultos desarrollado en la empresa Aura
Innovative Robotics.
4. Indicar la manipulabilidad del exoesqueleto de 6 grados de libertad para definir nuestro
espacio de trabajo.
5. Señalar los movimientos que dicho exoesqueleto puede realizar para optimizar la terapia
de rehabilitación.
6. Comparar el rango de movimiento del exoesqueleto con respecto a otros comerciales.
2
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 2
Lesión del plexo braquial
El plexo braquial es una red de nervios que provienen de la columna vertebral, formado por los
ramos anteriores de los nervios espinales C5, C6, C7, C8 y T1 (Figura 2.1), localizado en la
base del cuello, la región axilar y el brazo. Es el encargado de la inervación muscular y cutánea
del miembro torácico. Como se aprecia en la Figura 2.2, se distinguen tres troncos primarios:
el tronco superior (C5 y C6), el tronco medio (C7) y el tronco inferior (C8 y T1) [1].
La lesión del plexo braquial es una lesión neurológica que afecta a las extremidades superiores,
acarrea en un deterioro funcional, discapacidad física, ansiedad psicológica, entre otros. Las
lesiones del plexo braquial son ocasionadas por muchas razones, a raíz de caídas,
complicaciones obstétricas, tracción traumática, tumores, heridas penetrantes, síndrome de la
abertura torácica superior, complicaciones de las fracturas de clavícula, etc. Por ello, la
detención temprana de las lesiones del plexo braquial es esencial para aumentar la recuperación
funcional del miembro superior [2].
3
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
1.1. Etiología
Las lesiones en el plexo braquial suelen producirse por traumatismos violentos en el hombro o
en el miembro superior, ocasionando un estiramiento o desgarro de las fibras nerviosas, como
consecuencia de: aplicación de fuerza o tracción sobre el plexo braquial del feto durante el parto
por vía vaginal, accidentes laborales, accidentes de tráfico, problemas vasculares, caídas sobre
el hombro, accidentes deportivos, entre otras [3].
Además, las lesiones del plexo braquial que ocasionan daños en los nervios cervicales y
torácicos a razón de estiramientos excesivos, fractura del húmero, lesiones penetrantes, etc.,
pueden inducir al deterioro neurológico que pueden conllevar a parálisis del brazo y la
hipoestesia (disminución patológica de la sensibilidad) [1].
Mientras que esta patología en los niños, caso obstétrico, se presenta frecuentemente por la
distocia del hombro, cuando el feto queda trabado detrás de la pelvis materna luego de la salida
de la cabeza (ver figura 2.3). Lo cual, genera la lesión del plexo a raíz del estiramiento del
cuello del feto y el aumento del ángulo entre la cabeza y el cuello [4].
4
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Entre algunos de los factores de riesgo asociados con las lesiones del plexo braquial tenemos:
diabetes materna, aumento excesivo del peso materno, anomalías uterinas, antecedentes de
parálisis del plexo braquial, trabajo de parto prolongado, distocia del hombro, pelvis pequeña,
entre otros [4].
Los accidentes vehiculares (en su mayoría de dos ruedas) representan el 90% de los
casos.
Trauma industrial, a causa de caídas sobre los hombros de cierta altura.
Grandes caídas con estiramientos de cuello.
Asaltos con objetos punzantes.
Herida de bala
Lesiones iatrogénicas (alteración del estado del paciente producido por el médico), por
ejemplo, en cirugías donde pueden estar involucradas ciertas raíces nerviosas.
Donde la fisiopatología de la causa de las tres primeras, se atribuye a la tracción del plexo
causada por el ángulo anormal entre el hombro y el cuello después del impacto [5].
1.2. Epidemiología
Desde comienzos del siglo XXI, estudios han revelado la aparición de lesiones del plexo
braquial congénito de forma habitual, con una relación aproximada de 1,5 por cada 1.000
nacidos vivos, variando en un amplio rango (0,3-3 cada 1.000 nacidos vivos) según la región
geográfica, pese a los avances existentes en la práctica obstétrica [6].
5
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Otros estudios han revelado que la incidencia de la parálisis no ha disminuido a lo largo de las
últimas 3 décadas, e indican que el riesgo de que exista la lesión aumenta de forma proporcional
al peso del bebe al nacer. La frecuencia de aparición se establece en 0,9 por 1.000 en los
neonatos que pesaban menos de 4.000g, 1,8 por 1.000 en los que pesaban entre 4.000 y 4.500g
y 2,6 por 1.000 entre los que pesaban más de 4.500g [4].
Las lesiones del plexo braquial suelen afectar comúnmente a la población joven y trabajadora
que conduce a daños socioeconómicos alrededor de la familia. [7] Existen múltiples estudios
enfocados en adultos, los cuales indican que un 70% de las lesiones de plexo braquial
traumáticas se producen por accidentes de tráfico, siendo un 70% de éstos provenientes de
motocicletas (para un total de incidencia de 4,2% en accidentes de motocicleta que afectan al
1,2% de los poli traumatizados) [8].
La lesión del plexo braquial se clasifica a partir de la localización anatómica y del tipo de lesión.
Es por ello, que mediante un examen clínico se determina la ubicación de la lesión que puede
ser superior, inferior o total. Algunas posturas anormales y pérdidas de movimiento en
determinados músculos se les atribuyen a patrones de localización específicos, lo que conlleva
adicionalmente, tras la supervisión del tiempo de recuperación del paciente, a determinar el tipo
de lesión que vendría a ser por avulsión o ruptura [12].
En la Figura 2.4 se observan los troncos primarios del plexo braquial, conformados por los
ramos ventrales de los nervios cervicales C5, C6, C7, C8 y torácico T1:
6
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Localización de la lesión:
1. Lesión Superior: Esta lesión afecta al tronco superior y en ocasiones al tronco medio, el
paciente presenta incapacidad para la abducción del hombro, la rotación externa del
brazo y la supinación del antebrazo, además es posible una parálisis diafragmática. Se
suele denominar “mal hombro, buena mano” debido a que se asocia a una deficiencia
del hombro y un compromiso variable de la mano. Los grupos musculares involucrados
en esta lesión son los rotadores externos y abductores del hombro, los flexores del codo,
los supinadores del antebrazo y en ocasiones los extensores de la muñeca, a pesar de
que la movilidad de la muñeca y los dedos no se ve afectada [4].
2. Lesión Inferior: Esta lesión afecta al tronco inferior, el individuo presenta una
incapacidad para flexionar la muñeca o tomar objetos con la mano. Se suele denominar
“buen hombro, mala mano” debido a que la función del hombro suele ser buena y la de
la mano se ve comprometida. Los grupos musculares involucrados son los flexores de
la muñeca, los flexores de los dedos y los músculos intrínsecos de la mano [4].
3. Lesión Total: Esta lesión afecta a todas las raíces nerviosas del plexo braquial por lo
que encontramos una deficiencia funcional de todo el brazo y la mano. Es una de las
lesiones supraclaviculares más frecuentes, se caracteriza por un brazo flácido y por una
mano flácida, en forma de copa, o en “mano en garra” o “mano del predicador” [4].
7
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Posible parálisis
diafragmática
Buen funcionamiento
general del hombro y
deficiencia funcional de
la mano
Tracción del brazo Incapacidad para
Inferior Raíces C8 y T1
hacia arriba felxionar la muñeca o
tomar objetos con la
mano
Posible síndrome de
Honer
Traumatismo muy Deficiencia funcional en
Raíces C5, C6,
Total violento del brazo todo el brazo
C7, C8 y T1
completo Brazo y mano flácidos
Lesión por fractura del
Tracción nerviosa por hombro o de la clavícula,
Fascículos y sus
Infraclaviculares traumatismo del luxación glenohumeral o
ramificaciones
hombro fractura de la extremidad
proximal del húmero
Lesión radicular o de un
tronco primario, junto
Raíces, troncos Traumatismo grave a
con rotura del nervio
Doble nivel primarios y nivel supraclavicular e
musculocutáneo a su
ramificaciones infraclavicular
entrada en el
coracobraquial
Tabla 2.1. Clasificación por localización de la lesión del plexo braquial [4].
8
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Tipos de Lesión:
1. Lesión por avulsión: conocida también como avulsión radicular o lesión preganglionar
(antes del agujero de conjunción), arrancamiento de las raíces de la médula con muerte
de las neuronas medulares correspondientes (ver Figura 2.5). Es considerada como la
lesión de mayor gravedad, pueden afectar a ramas dorsales (sensitivas), ventrales
(motoras) o ambas (el nervio raquídeo o espinal incluye fibras sensitivas y motoras).
Frecuente en raíces C8-T1. (ver Figura 2.4) [9].
2. Lesión por ruptura: también conocida como lesión postganglionar (después del agujero
de conjunción), es una lesión nerviosa periférica que consiste en la rotura, disección o
compresión en algún punto a lo largo del nervio que constituye el plexo braquial. Esta
lesión de divide en supraclaviculares (afecta a troncos primarios) e infraclaviculares
(afecta a troncos secundarios o nervios terminales del plexo) (ver Figura 2.5). Frecuente
en raíces C5-C6. (ver Figura 2.4) [9].
9
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Los movimientos de la articulación del hombro se desarrollan en conjunto con la acción de los
músculos, tendones y ligamentos, que permiten orientar el miembro superior en los tres planos
del espacio (ver Figura 2.6). A continuación, se detallan los diferentes tipos de movimientos
existentes. Como se pueden aprecia en las Figuras 2.7 y 2.8.
10
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Flexión: disminución del ángulo formado entre dos huesos, dobla o pliega una parte sobre otra.
Extensión: aumento del ángulo formado entre dos huesos, devuelve a la posición anatómica
una parte desde su posición en flexión. Mientras que la denominada hiperextesión se presente
cuando de estira una parte más allá de su posición anatómica.
Abducción: movimiento lateral que separa el miembro superior del tronco (aleja una parte del
plano medio del cuerpo).
Aducción: movimiento que acerca el miembro superior hacia el tronco (mueve una parte hacia
el plano medio). Debido a la presencia del tronco, la aducción absoluta es mecánicamente
imposible por lo que, se precisa combinar este movimiento con la flexión / extensión como se
observa en la Figura 2.8.
Pronación: movimiento del antebrazo y mano que rota el radio medialmente, donde la palma
queda hacia atrás.
Supinación: movimiento del antebrazo y la mano que rota al radio lateralmente, donde la palma
queda hacia adelante.
Rotación: consiste en hacer girar un hueso sobre su propio eje Existen dos tipo de rotación, la
externa, lateral o rotación hacia fuera y la interna, rotación medial o rotación hacia adentro.
11
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a b
Figura 2.8. a) Aducción absoluta del brazo combinada con flexión/extensión. b) Supinación y pronación del
brazo
El diagnóstico de la parálisis del plexo braquial puede realizarse al momento del nacimiento o
poco después, partiendo tanto de la historia clínica como del examen físico realizado por el
médico. El neonato presenta brazo flácido, asimetrías en la región superior del brazo y ausencia
de los reflejos bicipital (a nivel medular en los segmentos C5-C6) y de Moro (o reflejo de
sobresalto en el bebé) en el brazo afectado, que vienen a ser respuestas involuntarias a la
estimulación. El examen clínico es de gran importancia en la detección de lesiones en el plexo
braquial, el especialista debe llevar un control de la evolución del neonato para evaluar el de
forma completa los hombros, codos, muñecas y demás dedos. Para tener una evaluación
completa, también se requiere analizar la flexión y extensión del hombro, codo, pulgar y dedos,
así como también evaluar en el hombre la flexión anterior, abducción, rotación interna y
rotación externa. Durante el período posnatal, varios casos de parálisis del plexo braquial
muestran cuadros clínicos semejantes, y de persistir estos casos deben ser derivados a clínicas
especializadas en PBO (Parálisis Obstétricas Braquiales). Es por ello que, según la extensión y
la localización de la lesión del plexo braquial se determina la gravedad del paciente y se plantea
la necesidad de intervención quirúrgica de ser necesario [4].
Marrero et.al. [10] describen el cuadro clínico para cada tipo de lesión del plexo braquial de la
siguiente forma:
Lesión Superior
Están paralizados [11] (ver Figura 2.9):
12
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Deltoides
Supraespinoso e infraespinoso
Redondo menor
Bíceps braquial
Braquiorradial (o braquirradial)
Braquial
Supinador corto y largo
Figura 2.9. Músculos del brazo. (a) Vista anterior, (b) Vista posterior. [12]
13
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Niño mayor:
Acortamiento del miembro afectado.
Atrofia muscular.
Hombro en aducción y rotación interna.
Limitación de la abducción del hombro.
Codo en abducción parcial y extensión.
Pliegue del codo hacia adentro y la mano hacia atrás.
Antebrazo pronado o supinado.
Ausencia de los reflejos bicipital y radial.
Anestesia en la región deltoidea y disminución de la sensibilidad en la cara lateral del
brazo y del antebrazo.
Lesión inferior
Quedan paralizados:
Flexores del carpo.
Flexores largos de los dedos.
Intrínsecos de la mano.
Ausencia temprana del reflejo de prensión.
Mano en garra por parálisis de los intrínsecos.
Atrofia de la eminencia tenar e hipotenar, que constituye la masa muscular de la base
del dedo pulgar.
La sensibilidad suele ser normal.
Lesión total
Parálisis sensitivo motora de toda la extremidad.
Miembro completamente fláccido e insensible.
1.5. Tratamientos
14
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
corporal mediante actividades lúdicas y de integración, así como reeducar la posición del
miembro superior del tronco. Posteriormente, se debe fortalecer la musculatura mediante el
desarrollo de diversas actividades a pesar de la ausencia de una función neuronal y muscular.
En algunos casos, se toma en consideración la inmovilización intermitente con la meta de
atenuar las fuerzas de gravedad que actúan sobre el brazo y contribuyen a aumentar la tracción
sobre el plexo braquial [4].
En 1903 se realizó la primera cirugía del plexo braquial y no fue hasta la década de los años
sesenta del siglo XX que la microcirugía causo un impacto significativo en el tratamiento de
dicha lesión [10]. Actualmente, las técnicas microquirúrgicas, los pegamentos biológicos y las
técnicas de diagnóstico intraoperatorias han proporcionado grandes mejoras en las fases
tempranas de la lesión e incluso se están haciendo reimplantaciones radiculares en el caso de
arrancamiento medular con buenos resultados.
15
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Supraespinoso, infraespinoso,
Nervio Supraescapular Deltoides y redondo menor
bíceps, braquial y braquiorradial
Supraespinoso, infraespinoso,
Combinado Redondo menor deltoides, redondo menos, bíceps
y braquial
Tabla 2.3. Casos más comunes de injertos de nervios para tratar la LPBS [14] [15].
Los dispositivos robóticos brindan tanto al terapeuta como al paciente optimizar el proceso dela
rehabilitación del miembro superior, ya que, comúnmente consisten en la repetición gradual de
movimientos según la articulación a ejercitar. Es por ello que, los exoesqueletos para
rehabilitación asisten a los pacientes según el grado de la lesión del plexo braquial a forma
práctica, continua, con plataformas interactivas, que, a su vez son capaces de mostrar el
desarrollo del proceso de cada individuo.
16
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 3
Estado del arte
Este capítulo tiene como finalidad evaluar el aporte de diversos exoesqueletos de rehabilitación
utilizados en investigaciones de los últimos años, como base para el posterior análisis
cinemático de la órtesis desarrollada en la empresa Aura Innovative Robotics.
17
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Como se muestra en la Figura 3.1, en la década de los setenta y principios de los ochenta, los
primeros dispositivos para terapia de rehabilitación como las máquinas Lido y Biodex
consistían en palancas y barras conectadas al motor con el fin de trabajar diferentes miembros,
como por ejemplo la flexión y extensión del codo, casos de abducción y aducción, permitiendo
que los pacientes ejerciten una articulación específica. El dinamómetro que tienen, detecta el
momento de rotación y la rotación del brazo del paciente para mostrar la información y
documentar los resultados [16].
Figura 3.1. Precursores de los dispositivos de terapia robótica, los tres de la izquierda: la honda de Suecia
(Swedish sling), el brazo monopatín (arm skateboard) y el soporte de brazo móvil JAECO. El dispositivo de la
derecha es el Dinamómetro Biodex Activo (1 GDL que se ajusta dependiendo de la articulación) [16].
18
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 3.2. Algunos de los primeros dispositivos utilizados en terapias robóticas, de izquierda a derecha:
MITMANUS, MIME, la Guía ARM [16].
Algunos médicos expresaron su escepticismo en relación a la idea de que los robots pudieran
ayudar a cumplir las metas de la rehabilitación. Pensaban que los robots no podían igualar ni la
experiencia ni las habilidades de los terapeutas, consideraban que las terapias se tenían que
realizar con movimientos complejos, en estado de alerta y con una determinada percepción
acorde al problema de cada paciente. De igual manera, los catalogaban como inseguros, no
inteligentes y no les trasmitía confianza el hecho de que en un futuro pudieran llegar a sustituir
el trabajo de los especialistas [16].
Gracias al incentivo de dispositivos como MITMANUS, MIME, la Guía ARM (ver Figura 3.2)
y a la colaboración de grandes especialistas de la rehabilitación este campo de trabajo motivo a
los ingenieros a reforzar la robótica, las teorías de control y los enfoques computacionales para
abarcar otros problemas de rehabilitación como son las lesiones neurológicas [16].
Finalmente, con la meta de obtener mejores resultados clínicos a través de la terapia robótica
se estudia de forma meticulosa diferentes lesiones que afectan al miembro superior. Es por ello,
que la rehabilitación va a depender del tipo de problema a atacar y de la articulación implicada
[16].
19
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En 2005, el instituto de investigación MBEC (Medical & Biological Engineering & Computing)
[17], describe de forma general una visión de diferentes dispositivos robóticos para
rehabilitación del brazo con nuevas tecnologías, que mediante movimientos repetitivos
orientados pueden mejorar la fuerza muscular y la coordinación de movimientos en pacientes
con patologías neurológicas. De igual forma, se realiza una comparación en cuanto a sus
funciones técnicas, aplicabilidad clínica y si existen o no resultados clínicos.
20
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Swedish Helparm (Figura 3.3), es un sistema pasivo basado en el contrapeso, que funciona con
un sistema de cables y poleas que sostienen a los brazos de forma independiente para realizar
las tareas, brinda asistencia funcional y de entrenamiento muscular, se puede utilizar para
practicar actividades de la vida diaria en la clínica o en el domicilio del paciente [17].
En la Figura 3.4, tenemos los dispositivos BTE (Baltimore Therapeutic Equipment) que de igual
forma ayudan al paciente a realizar entrenamientos para actividades de la vida cotidiana y a su
vez permiten realizar terapias para las extremidades inferiores y superiores mediante una gran
variedad de movimientos, con la diferencia de ser un sistema activo [17].
El sistema de rehabilitación asistida ARM (Figura 3.5), puede ser empleado en pacientes con
lesiones cerebrales crónicas debido a que al no implicar ninguna impulsión de movimiento se
clasifica como sistema pasivo. La guía ARM tiene un grado de libertad pasivo de traslación
(para los movimientos de la mano hacia el hombro) y otro semi activo rotatorio (corresponde a
la rotación del hombro en el plano vertical) [17].
Para desarrollar un dispositivo que ayude a la rehabilitación del brazo, se requiere considerar a
modo clínico aspectos psicológicos, médicos y ergonómicos. La presencia del terapeuta durante
la rehabilitación con ayuda del robot de asistencia proporciona confianza en el paciente y por
ende optimiza la terapia. El diseño y aplicación del dispositivo robótico se deje adaptar a cada
individuo, ajustar las longitudes de segmentos, el rango de movimiento y el número de grados
de libertad (GDL) según la articulación a tratar. Al aumentar el número de GDL tendremos
mayor variedad de los movimientos a realizar, pero a su vez implicaría un dispositivo más
complejo y costoso. Es imprescindible que el robot de rehabilitación sea compatible con las
terapias tradicionales y procese la información tanto de los sensores del dispositivo como del
progreso de evolución del paciente [17].
21
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 3.4. BTE device Simulator II for training of ADL (photo courtesy of BTE Technologies, Inc.)
Figura 3.5. ARM guide (1999a) Figura 3.6. Arm rehabilitation robot ARMin (2004)
22
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Por ello, se estudiaron diversas técnicas de control para obtener soluciones que garanticen un
alto seguimiento dinámico ante la presencia de no linealidades e incertidumbres. Se trabajó con
un robot de asistencia que simula los movimientos más importantes de las extremidades
superiores con 5 grados de libertad (GDL): abducción/aducción del hombro, flexión/extensión
del hombro, flexión/extensión del codo, flexión/extensión y rotación interna/externa de la
muñeca. Con el fin de que las articulaciones ejecuten movimientos complejos de tareas
cotidianas [18].
Figura 3.7. Diagrama esquemático y fotografía real del exoesqueleto de miembro superior de 5 GDL.
Debido a que los movimientos del brazo humano son no lineales se tiene que recurrir a técnicas
de control complejas, ya que, los modos convencionales de control cumplen escasamente con
los requerimientos de dicha tarea. Varios trabajos suelen emplear estrategias de control como:
control por métodos deslizantes (SLV), esquemas de control por aprendizaje iterativo, siendo
uno de los más conocidos es el control por modelos dinámicos simplificados. Las cuales, en la
mayoría de los casos, necesitan que se tome en cuenta el modelado del error y de las
incertidumbres, así como también se deben considerar factores de calibración inicial, masa y
tamaño del brazo según el paciente, ya que varían de persona a persona [18].
Por tanto, el presente trabajo muestra tres estrategias de control adaptativo para un exoesqueleto
de miembro superior de 5 GDL que brinde mejores garantías de seguridad:
1. Seguimiento de salidas bajo incertidumbres, discrepancia de los parámetros de varianza
y mal funcionamiento de hardware.
23
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Finalmente, basándonos en la meta de obtener movimientos con una posición específica y una
cierta protección de seguridad, el método del control adaptativo para un exoesqueleto de 5 GDL
facilita la planificación de trayectorias a partir de la intuición, de forma simple y con un bajo
coste tecnológico [18].
Por otro lado, en 2013, Riaan Stopforth en conjunto con el grupo de investigación de
mecatrónica y robótica (MR2G) unidad de bio-ingeniería de la Universidad de KwaZulu-Natal
(Durban, South Africa) desarrolló un control personalizado y ajustable de un exoesqueleto de
miembro superior [19].
Este trabajo estudia los aspectos mecánicos, cinemáticos y biológicos necesarios para el diseño
y desarrollo de un exoesqueleto de miembro superior personalizado y ajustable. El exoesqueleto
propuesto toma en cuenta la capacidad de ajuste, seguridad, funcionabilidad y soporte de
objetos, en base a los movimientos de flexión/extensión, abducción/aducción y
supinación/pronación en un determinado tiempo según pequeñas alteraciones [19].
24
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En la figura 3.9, encontramos la interacción básica del sistema con la PC, donde las líneas azules
muestran la alimentación del sistema y las rojas el flujo de la señal electrónica. El PC envía
órdenes al microcontrolador (en este caso un ATMEGA 1.280), para crear la realimentación
del sistema y se colocan interruptores de seguridad, en caso de fallas, a dicho microcontrolador
a modo de seguridad y así evitar que el exoesqueleto sobrepase los límites humanos y
mecánicos. Se decidió implementar un control PID a modo de prueba, donde la ganancia
proporcional redujo el tiempo de subida, la ganancia integral el error en estado estacionario, y
la ganancia derivativa permite mejorar la estabilidad del sistema, la respuesta transitoria y
reducir el sobre impulso [19].
25
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Se implementó un control PID adaptativo que varía según los diferentes pesos y fuerzas de los
brazos, así como también influyen los movimientos sobre el sistema. Por tanto, se debe realizar
un control de impedancia ajustado a cada individuo para alcanzar un mejor desempeño del
exoesqueleto del miembro superior. Para el desarrollo del exoesqueleto se consideró un factor
de bajo costo y un peso específico para poder ser transportado fácilmente a lugares remotos
[19].
De igual forma, en 2013 [20], Yu, W., & Rosen, J. en conjunto con la organización IEEE
(Institute of Electrical and Electronics Engineers) implementaron un control neuronal PID aun
exoesqueleto de miembro superior. Con el fin minimizar el error en estado estacionario ante
incertidumbres, en el control PID se necesita insertar una ganancia elevada en la componente
integral del control, otra solución, es la de agregar un compensador neuronal en el control PD
añadiendo un valor alto de ganancia en la componente derivativa de control, pero ambas
soluciones deterioran la respuesta transitoria en el control del robot. Por ello, el presente trabajo
le añade un compensador neuronal a un control PID lineal clásico, con la intención de combinar
los beneficios de cada control. Por una parte, la estabilidad asintótica semi-global que presenta
el control PID neuronal y la estabilidad local al agregar observadores para estimar la velocidad
(teorema que permite el entrenamiento de los pesos neuronales y un método de selección
explícita de las ganancias PID).
El control PID es ampliamente conocido y utilizado en diversos sistemas, debido a que tiene
características que facilitan su aplicación, no necesita modelar dinámicamente la planta y
permite ajustar las componentes de las ganancias por separado. Uno de los problemas más
significativos de dicho controlador, está en la componente integral, debido a que se pueden
presentar fallos en la estabilidad ante la presencia de perturbaciones de ruido. Por ello, es que
muchos robots manipuladores prefieren emplear solamente un control PD o añadir ganancias
muy pequeñas en el control PID [20].
Una solución propuesta para garantizar la estabilidad asintótica del sistema en casos no ideales
es cambiar el control PID lineal por uno no lineal. Se han venido implementando métodos
alternativos como son: el control PD con ganancias variables, el control PD con ganancias no
lineales, el control PD con compensación en modo deslizante, etc., que a pesar de que puedan
brindar cierta estabilidad son muy complejos [20].
26
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En la figura 3.10, se aprecia el control del sistema, donde encontramos por un lado el
controlador PID y por otro el compensador neuronal que exoesqueleto (actuador). De la misma
forma, tenemos un observador en velocidad que posteriormente se suma al algoritmo de
impedancias (luego de medir previamente las fuerzas) a modo de realimentar del sistema y
determinar los cálculos del control [20].
Como idea principal, el trabajo propone extender un control PD a un nuevo control PID
neuronal, aplicado en el robot exoesqueleto de 7 GDL de la Universidad de California (Santa
Cruz, UCSC) (Figura 3.11). Donde se demuestra que no se necesitan valores elevados en las
ganancias de las componentes integral y derivativa, se garantiza estabilidad semi-global
asintótica y estabilidad local asintótica a través de un observador de velocidad. Tenemos que,
el exoesqueleto de la UCSC está conformado por 3 GDL para el hombro (J1-J3), uno para el
codo (J4) y tres para la muñeca (J5-J7). Las articulaciones esféricas J1, J2 y J3 son responsables
de la flexión-extensión, aducción y rotación interna-externa del hombro; J4 es la superficie
articular responsable de la flexión-extensión del codo y las articulaciones esféricas J5, J6 y J7
son las responsables de la flexión-extensión, de la muñeca, la pronación-supinación y la
desviación radial-cubital [20].
27
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En conclusión, a partir de los resultados obtenidos (ver gráfica 3.2.6) se obtuvieron excelentes
resultados en el control PID neuronal, en relación a otros, donde a pesar de existir una leve
oscilación en la respuesta del sistema se corrige inmediatamente y se alcanza el estado
estacionario [20].
28
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
ARMin
Es un robot de rehabilitación de brazo con estructura de exoesqueleto (se asemeja a la anatomía
del brazo humano) y equipado con sensores de posición y de fuerza, el cual es aplicable en
clínicas para ayudar, mejorar, evaluar y documentar las terapias de rehabilitación. El dispositivo
puede funcionar en modo de terapia de movilización pasiva, terapia de juego y de formación
orientado a la tarea [21].
Como se aprecia en la Figura 3.13, la versión ARMin II tiene seis grados de libertad, donde el
Eje 1 corresponde a la rotación vertical del hombro, el Eje 2 representa la rotación horizontal
del hombro, el Eje 3 es responsable de la rotación interna / externa del hombro, el Eje 4 está
relacionado a la flexión / extensión del codo, el Eje 5 corresponde a la pronación / supinación
del antebrazo y por último el Eje 6 la flexión /extensión de la mano. En la Tabla 21, se muestra
el rango de movimiento del exoesqueleto ARMin [21].
Tabla 3.3. Rango de movimiento del robot ARMin para rehabilitación de miembro superior. [21]
29
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 3.15. Interfaz para la terapia del brazo, de apoyo (juego con pelota y el laberinto), escenario para
entrenamiento de actividades cotidianas [22].
Figura 3.16. Robot para rehabilitación de brazo ARMin III (Octubre 2007), con 4 GDL accionados para
controlar el hombro y flexión/extensión del codo [23].
30
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 3.17. Relación de los tres exoesqueletos robóticos de la Terapia Armeo (ArmeoPower, ArmeoSpring y
ArmeoBoom) [22].
31
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a modo de juego que hace que la terapia sea divertida. Según los datos técnicos, el exoesqueleto
ArmeoPower de la compañía Hocoma AG, permite realizar los siguientes movimientos: [25]
Abducción horizontal del hombro con respecto al plano sagital: -169º a +50º
Flexión/extensión del hombro en relación al plano frontal: +40º a +120º
Rotación interna/externa del hombro: 0º a 90º
Flexión/extensión del codo: 0º a 100º
Pronación/supinación del antebrazo: -60º a 60º
Flexión/extensión de la muñeca: -60º a 60º
Y adicionalmente, agarre de la mano por sensor de presión [25].
El software Armeo brinda a los pacientes una plataforma que los motiva mientras realizan el
entrenamiento para actividades del día a día y a su vez, aumenta la eficiencia del trabajo del
terapeuta (ver Figura 3.19). Proporciona una gran variedad de juegos, aumenta el rango de
movimiento, permite coordinar los movimientos, realizar funciones de agarre, fuerza,
resistencia y funciones cognitivas. Donde la dificultad y los niveles de los juegos se adaptan a
las capacidades y progresos de cada individuo [24].
32
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Artromot S3
El dispositivo ASTROMOT–S3 (ver Figura 3.20) consiste en una tablilla de movilización
monitorizada que contribuye al tratamiento médico-terapéutico, a través de movimientos
pasivos continuos (MPM) de la articulación del húmero. El objetivo principal está enfocado en
la movilidad de las articulaciones, para evitar daños por inmovilización (prevención de
trombosis, embolias, entre otros), agilizar el proceso curativo y obtener mejoras en los
resultados de la rehabilitación. Características técnicas del rango de movimiento y otras
generales: [26]
33
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
muchos otros beneficios, integra uno de los dispositivos tecnológicos comerciales empleados
actualmente [27].
34
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 4
Descripción del Exoesqueleto de Rehabilitación a utilizar
Figura 4.1. Exoesqueleto de miembro superior para rehabilitación, desarrollado en la empresa Aura Innovative
Robotics.
Uno de los objetivos principales para el diseño del exoesqueleto de miembro superior consistía
en un dispositivo ligero, confiable, de fácil uso, eficiente, entre otros, que según el número de
grados de libertad pudiese brindar una gran gama de movimientos para mejorar la rehabilitación
del paciente, facilitando y optimizando el trabajo del terapeuta.
35
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Es por ello, que en la Figura 4.2 se aprecia uno de los primeros prototipos de la órtesis de
hombro diseñados en Autodesk Inventor 2015 (software de diseño mecánico y CAD 3D). Se
puede apreciar que es un dispositivo ajustable según el tamaño del paciente a tratar bien sea por
la altura, tamaño del brazo o del antebrazo. Así como también, a pesar de ser un diseño sencillo
para la articulación del hombro, permite efectuar los siguientes tipos de desplazamiento: la
flexión, abducción, aducción y extensión. Que en conjunto con un buen sistema de control se
logra alcanzar los objetivos deseados con respecto a diversos tipos de rehabilitaciones acorde a
las carencias específicas existentes.
Figura 4.2. Diseño en inventor de las piezas del exoesqueleto de miembro superior [29].
La mayor parte del exoesqueleto para rehabilitación del miembro superior está conformada por
una estructura a base de ABS (plástico muy duro y resistente a altas temperaturas), partiendo
de los planos realizados previamente en inventor se procedió a imprimir cada pieza en una
impresora 3D. Posteriormente, se fueron agregando los elementos de soporte como son los
tornillos, turcas, cierres mágicos para sostener el brazo del paciente junto al dispositivo, el usillo
para la primera articulación prismática (a modo de tornillo sin fin), etc. Adicionalmente, y no
menos importante, se colocaron los 6 motores en cada una de las articulaciones pertinentes
como se contempla en las figuras 4.1 y 4.2 [30]. La órtesis utiliza servomotores Dynamixel RX-
64 (Figura 4.3), los cuales son actuadores inteligentes con la capacidad de realizar un
seguimiento de la velocidad, temperatura, posición del eje, voltaje y carga. Estos servomotores,
según la especificación técnica del producto [31], poseen una relación de reducción 200:1,
también tienen un microcontrolador incorporado que gracias al algoritmo de control se puede
ajustar la posición del eje, la velocidad y la fuerza de la respuesta de cada motor de forma
individual. Por ejemplo, entre las herramientas de configuración del servomotor se puede
modificar la limitación del par, limitaciones de entrada o velocidad de movimiento, entre otros.
36
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
37
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 5
Cinemática del Exoesqueleto de rehabilitación de miembro
superior
Figura 5.1. Primer prototipo del exoesqueleto para rehabilitación de miembro superior realizado en Aura
Innovative Robotics.
38
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 5.2. Diagrama del Exoesqueleto para realizar el cálculo cinemático del robot.
Figura 5.3. Dos barras rígidas consecutivas unidas por una articulación (un grado de libertad)
39
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Encontramos,
i-1
Ai = Rotz(θi) T(0,0, di) T(ai,0,0) Rotx(αi)
Por tanto,
Matriz de D-H
Es por ello, que siguiendo los pasos de D-H al diagrama de la Figura 5.2, obtuvimos los
resultados de la tabla 5.1.
Articulación Ɵ d a α
1 π/2 L0 + q1 0 π/2
2 q2 + π/2 L1 -a -π/2
3 q3 0 a π/2
4 q4 L2 -a -π/2
5 q5 0 a π/2
6 q6 L3 -a 0
Tabla 5.1. Análisis Cinemático por el método D-H.
Se desarrolló un código en Matlab (ver Anexo A), donde se obtuvo la matriz de transformación
homogénea T0 = 0A1 1A2 2A3 3A4 4A5 5A6. Si,
40
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
0 0 1 𝐿1 + 𝐿2 + 𝐿3
0 −1 0 0
𝐓ο = [ ]
1 0 0 𝐿0 − 𝑎
0 0 0 1
A diferencia del método de DH que sitúa los ejes xyz y el origen del sistema de coordenadas,
éste método por Matrices de Desplazamiento precisa solo identificar el eje de cada articulación
(ejes de desplazamiento), simplificando de forma notable el procedimiento. Luego se calculan
las matrices Di ligadas a cada una de las articulaciones que vienen a ser los desplazamientos
generalizados en el espacio, bien sean rotaciones o traslaciones, que van a depender de los
parámetros que definen los ejes de la articulación. Después, al efectuar el producto de las
matrices Di se obtiene el modelo cinemático que se le añade a la matriz de transformación
homogénea T0, siendo este último la localización del extremo del robot cuando todos los de
libertad adoptan el valor cero [35].
Es por ello que, primero situamos el robot en su posición cero (qi = 0, ∀i), y se obtiene la matriz
de transformación homogénea T0, luego se obtienen las matrices de desplazamiento para lo
cual:
41
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Siendo el operador skew el que a cada vector k = (kx,ky,kz) asocia la matriz asimétrica:
Al igual que en el caso de D-H, tomando el esquema de la Figura 5.2 para calcular la cinemática
directa, obtenemos los resultados de la tabla 5.2.
42
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
0 0 1 𝐿1 + 𝐿2 + 𝐿3
0 −1 0 0
𝐓ο = [ ]
1 0 0 𝐿0 − 𝑎
0 0 0 1
Se busca determinar la configuración que deben adoptar las coordenadas articulares qi = {q1,
q2, q3, q4, q5, q6} del exoesqueleto de rehabilitación de miembro superior para obtener una
posición y orientación deseada del extremo del robot. Según la configuración de la órtesis
pueden existir múltiples soluciones [34].
Al momento de buscar la cinemática inversa del robot es más conveniente encontrar una
solución cerrada, buscar la relación: qk = fk (x, y, z, φ, θ, ψ), k = 1…n (GDL). Debido a que una
solución de tipo iterativo retrasa el resultado en tiempo real y en los casos de solución de casos
cerrado se pueden fijar restricciones que aseguren una respuesta adecuada, ya que a diferencia
del caso del cálculo de la cinemática directa podemos encontrar diferentes soluciones [q1, . . . ,
qn] T que posicionen y orienten al efector final al mismo nodo [34].
Para el presente trabajo se utilizó la Robotics Toolbox de Peter Corke para Matlab, para resolver
el problema de cinemática inversa del exoesqueleto de rehabilitación de miembro superior [36].
Puesto que T = 0A1 1A2 2A3 , tendremos (0A1)-1 T = 1A2 2A3 y (1A2)-1 (0A1)-1 T = 2A3
43
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
0 −𝑆2 𝐶2 0 𝑛𝑥 𝑜𝑥 𝑎𝑥 𝑝𝑥 𝐶3 0 𝑆3 𝑎𝐶3
−1 0 0 0 𝑛𝑦 𝑜𝑦 𝑎𝑦 𝑝𝑦 𝑆3 0 −𝐶3 𝑎𝑆3
[ ][ ]=[ ]
0 −𝐶2 −𝑆2 0 𝑛𝑧 𝑜𝑧 𝑎𝑧 𝑝𝑧 0 1 0 0
𝐿1 𝑎𝑆2 𝐿𝑜 + 𝑞1 − 𝑎𝐶2 1 0 0 0 1 0 0 0 1
Se tiene finalmente,
𝐿1 𝑝𝑥 𝑎 𝑆2 𝑝𝑦 𝑝𝑦 (𝐿1−𝑝𝑥)
q1 = 𝑎𝐶2 − 𝐿0 − 𝑝𝑧
− 𝑝𝑧
q2 = 𝑡𝑎𝑛−1 (− 𝑝𝑧
) q3 = 𝑠𝑖𝑛−1 𝑎
Es por ello que, partiendo de la premisa de que los tres últimos grados de libertad se cortan en
un punto, luego se analizaron las variables articulares del robot (q4, q5, q6) con el propósito de
poder orientar libremente el robot en el espacio de trabajo [34].
Denominamos 0R6 (orientación deseada del extremo del robot) a la submatriz de rotación de
0
T6 y obtenemos,
0
R6 = [n o a] = 0R3 3R6
Donde i-1Ri es la submatriz de rotación de i-1
Ai. Y a partir de los valores conocidos de q1, q2
y q3 tenemos la siguiente relación:
3
R6 = [rij] = (0R3)-1 0R6 = (0R3)T [n o a]
0
Siendo R3 = 0R1 1R2 2R3 y 3
R6 = 3R4 4R5 5R6.
𝐶4 0 −𝑆4 𝐶5 0 𝑆5 𝐶6 −𝑆6 0
3 4 5
R4 = [ 𝑆4 0 𝐶4 ] R5 = [ 𝑆5 0 −𝐶5] R6 = [ 𝑆6 𝐶6 0]
0 −1 0 0 1 0 0 0 1
44
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Por tanto,
𝑟23
r13 = C4 S5 q4 = 𝑡𝑎𝑛−1 ( 𝑟13)
r23 = S4 S5
Por tanto, dada una posición y orientación deseada para resolver el problema de cinemática
inversa por desacoplamiento, primero se procede a calcular las 3 primeras variables articulares
(q1, q2, q3) en el punto, bien sea por métodos geométricos o a partir de matrices de
transformación homogénea y después con la orientación se determina el valor de las variables
articulares restantes (q4, q5, q6) [34].
Consiste en determinar las zonas del espacio de trabajo donde el robot puede moverse en todas
las direcciones. El índice de maniobrabilidad, es una medida escalar de la dificultad del
movimiento del robot. En el presente trabajo, se estudiaran algunos índices de desempeño
definidos sobre el espacio de trabajo. La mayoría de estos índices están definidos en función de
la matriz Jacobiana del manipulador, que establecen una relación entre las velocidades y fuerzas
de las articulaciones del extremo del robot [37].
Uno de los índices más comunes que utilizaremos es el propuesto por Yoshikawa (1985b), cuyo
propósito es medir la capacidad del robot para generar velocidades en el extremo del robot
según una determinada configuración. En el caso de robots redundantes la manipulabilidad está
definida por 𝑤 = √det(𝑱 𝑱𝑇 ) . Mientras, que en robots no redundantes w = |det(J)|, por tanto
la manipulabilidad es equivalente al producto de los valores singulares, w = σ1σ2 . . . σm
(corresponden a las magnitudes de los ejes principales del elipsoide de velocidad, ver Figura
5.4). Por otro lado, en el caso de tener mayores valores de la manipulabilidad existirá una mejor
capacidad de movimiento del efector final, lo cual implica una conversión eficiente de las
velocidades articulares. Inversamente, al tener valores pequeños implica una conversión
ineficiente y en el caso en que la manipulabilidad es cero estaremos en presencia de una
45
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
De igual forma, en el código desarrollado en Matlab (ver Anexo A), se aprecian las jacobianas
del sistema y por ende la representación de la dirección del vector velocidad tanto lineal como
angular del efector final, dado por el giro de cada articulación. Al utilizar “Jacobn(robot,q)”
obtenemos la Jacobiana geométrica en el sistema de referencia del extremo del robot (que nos
permite conseguir la relación entre las velocidades articulares y las velocidades lineales y
angulares del extremo del robot) [34].
46
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 5.5. Jacobiana geométrica que relaciona las velocidades articulares y las velocidades lineales y angulares
del extremo del robot.
Existen dos tipos de singularidades, en los límites del espacio de trabajo cuando el extremo del
robot está en un punto límite hacia donde el mismo no se podrá desplazar o en el interior del
espacio de trabajo por el alineamiento de dos o más ejes de las articulaciones del robot [34].
El desarrollo del cálculo cinemático del exoesqueleto para rehabilitación del miembro superior
se realizó en la plataforma de MATLAB versión R2015b, que facilita y optimiza el trabajo
matemático computacional. De igual forma, la herramienta permite realizar gráficos para una
mejor visualización de la información, brinda una gran gama de librerías, entre otras cosas.
Para el presente trabajo se utilizó la Toolbok de Peter Corke [40] referente a robótica, la cual
proporciona una gran gama de funciones que simplifican y corroboran algunos cálculos
tediosos.
47
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 5.6. Interfaz en Matlab para el análisis del espacio de trabajo del exoesqueleto de miembro superior.
Como se puede observar, código desarrollado permite ajustar las coordenadas articulares del
robot y el tamaño de los eslabones, estos últimos van a depender del tamaño del brazo de cada
paciente. Es por ello, que al ajustar esos parámetros obtenemos la matriz de transformada
48
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Posteriormente, se fijarán los grados a los cuales puede llegar o no el exoesqueleto de miembro
superior para rehabilitación al realizar los diversos movimientos de las articulaciones
explicadas anteriormente. De igual manera, se tomará en cuenta en el siguiente capítulo el
estudio del espacio de trabajo para pacientes de diferentes edades, pesos y tamaños. Con el
objetivo de definir el área del entorno donde la órtesis puede trabajar sin encontrarse con puntos
singulares, se evita que se presenten inconvenientes de desplazamiento.
49
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 6
Análisis de Resultados
Es este capítulo, se muestran los rangos de movimiento del nuevo exoesqueleto de miembro
superior. A partir del estudio de la cinemática directa e inversa, se procedió a analizar el
determinante de la matriz jacobiana y calcular el índice de manipulabilidad de Yoshikawa para
diferentes pacientes, con el fin de tener una comparación detallada de los factores que pueden
afectar al espacio de trabajo del nuevo exoesqueleto para rehabilitación de miembro superior.
Para acotar el rango de movimiento de la nueva órtesis se requiere, en primer lugar, ajustar los
valores de las coordenadas articulaciones según el tope físico previamente fijado. Luego, es
necesario evaluar el espacio de trabajo en dicho rango de movimiento dependiendo del tamaño
de cada paciente. Y por último, se debe realizar una comparación de que tan factible es ejecutar
o no ciertos movimientos al momento de una terapia de rehabilitación. En la tabla 6.1, se
muestran los topes mecánicos existentes para cada articulación implicada del nuevo
exoesqueleto.
Hay que destacar y tomar en cuenta que según el paciente pueden presentarse situaciones
particulares que afecten este rango de movimiento previamente fijado. Al presentarse ocasiones
donde los individuos son de estatura muy reducida o por el contrario su altura les impide entrar
50
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
en la órtesis, o cuando por ejemplo existe un paciente que tiene desarrollado el músculo del
bíceps significativamente el rango de movilidad de una flexión / extensión del codo se va a
reducir, etc.
a b
Figura 6.1. a) Posición horizontal (q3= π/2) y b) Posición de reposo del exoesqueleto.
En la Tabla 6.2, se observan las características pertinentes de seis individuos sanos con
diferentes edades y estaturas para comparar resultados en relación a la posición y orientación
alcanzada por el efector final, el índice de manipulabilidad en ciertos casos, así como también
el determinante de la matriz jacobiana, entre otros, al variar las coordenadas articulares. Es por
ello que, se fijó como distancia a la base un L0=70cm haciendo referencia a que la medida se
toma desde el piso hasta la base de la articulación prismática. Luego, se varían los parámetros
correspondientes a la primera articulación traslacional (coordenada articular q1) y las distancias
L1 (acorde a la separación entre la coordenada articular q2 y q3), L2 (responsable de la medida
del antebrazo y L3 (según el tamaño del brazo).
51
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
L0 o base q1 o L2 o L3 o
Edad Altura
Nº del sistema Traslacional L1 (cm) Antebrazo Brazo
(años) (cm)
(cm) (cm) (cm) (cm)
1 10 73 70 10 18 10,66 13,58
2 15 79 70 10 20 11,53 14,69
3 20 85 70 15 20 12,41 15,81
4 25 161 70 20 22 23,5 29,95
5 30 179 70 20 22 26,13 33,29
6 35 183 70 20 22 26,72 34,04
Tabla 6.2. Rangos de movimientos de las coordenadas articulares.
Para los diferentes análisis del cálculo de la manipulabilidad del exoesqueleto se propone variar
hasta tres coordenadas articulares (del espacio articular), debido a que con un mayor número el
programa se hace muy lento e ineficiente. Para lo cual, se hicieron tres iteraciones que
cumplieran con la resolución deseada de una circunferencia de -π a π (-180º a 180º), como se
detalla en el Anexo A, lo que permite trabajar con una, dos o tres articulaciones al momento de
implementar varias pruebas que muestren de forma gráfica cómo varía el índice de
manipulabilidad de la nueva órtesis. Y por ende, concluir como afectan estas variaciones en el
efector final de nuestro robot para así, definir el espacio de trabajo. Además, al describir un
volumen en forma de elipse hay que tener presente que una cosa es la manipulabilidad del robot
de 6 GDL y otra es la del brazo humano, ya que, algunos tipos de movimientos aunque abarquen
un gran área en las gráficas, en la realidad uno se debe ajustar a los límites colocados en la
Tabla 6.1 para realizar un análisis de los resultados acorde a los objetivos planteados.
Tomando como punto de partida los diferentes tipos de movimientos (Tabla 6.1) que el
exoesqueleto puede realizar se implementaron una serie de pruebas donde se describe una esfera
para que se vea más claro el espacio de trabajo, a pesar de que no todas las coordenadas
articulares pueden o deben alcanzar ciertos grados de movimiento. Por tanto, se busca tener una
gran gama de experimentos donde se aprecie el entorno abarcable por el exoesqueleto para
casos particulares.
Experimento 1:
Entonces, para comenzar las pruebas con el exoesqueleto de miembro superior para
rehabilitación se propone variar sólo las coordenadas articulares q2, q3 y q5 en un rango entre -
180º a 180º. Se tomaron a los individuos sanos 1 y 6, para tomar como punto de partida la menor
52
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
y mayor estatura, con el fin de observar la posición y orientación del efector final del
exoesqueleto a partir del cálculo de la cinemática directa y la gráfica del índice de Yoshikawa.
Figura 6.2. Posición final del efector final del exoesqueleto para el individuo sano 1.
Figura 6.3. Posición final del efector final del exoesqueleto para el individuo sano 6.
53
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En ambas imágenes se observa como la posición y orientación del extremo del robot varía según
al tamaño de cada individuo, también vemos un cambio significativo en el estudio del índice
de manipulabilidad debido a que el área de trabajo aumenta su rango de movimiento al existir
un mayor radio de ejecución por parte del individuo (a mayor distancia el radio aumenta).
a b c
Figura 6.4. Movimientos terapéuticos de: a) abducción/aducción, flexión/extensión del hombro y
flexión/extensión del codo, b) abducción/aducción del hombro y de rotación interna y externa del antebrazo, c)
flexión/extensión del hombro y el codo.
Experimento 2:
Tomando las medidas de los individuos 4 y 5, se propone realizar un barrido de las coordenadas
articulares q2 y q4 entre (-180º a 180º) correspondientes de la abducción/aducción del hombro
y de rotación interna y externa del antebrazo. Las demás coordenadas articulares q1, q3, q5 y q6
se colocan en la posición de reposo (Figura 6.1.b).
54
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a b
Figura 6.5. Movimientos de abducción/aducción del hombro y de rotación interna y externa del
antebrazo para: a) Individuo 4, y para el b) individuo 5.
Experimento 3:
Para este caso en particular solo se tomó en consideración los datos del individuo 2 y 5, con la
meta de realizar ejercicios terapéuticos como los de la Figura 6.4.c. Para lo cual, se propone
hacer un barrido en torno a las coordenadas articulares q3 y q5 dentro de los rangos (-90º a 90º)
y (-180º a 180º), mientras que las demás coordenadas articulares q1, q2, q4 y q6 se colocan en la
posición de reposo.
55
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a b
Figura 6.6. Individuo 2, barrido de las coordenadas articulares q 3 y q5 dentro de los rangos: a) de -90º a
90º y b) de -180º a 180º.
a b
Figura 6.7. Individuo 5, barrido de las coordenadas articulares q 3 y q5 dentro de los rangos: a) de -90º a
90º y b) de -180º a 180º.
Partiendo de las figuras 6.6 y 6.7, observamos que al hacer los movimientos de
flexión/extensión del hombro y el codo el exoesqueleto de miembro superior se mueve a lo
largo del plano sagital, por tanto el área dibujado por el efector final recorrerá una media luna
que es como se aprecia en la figura donde nuestro robot podrá desplazarse sin problema. De
igual manera, estaremos en presencia de puntos singulares si sobrepasamos los límites
propuestos al inicio. Una de las razones del ¿por qué es importante ajustarse a los límites
56
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
propuestos para cada articulación? es para evitar que ocurra, por ejemplo, un estiramiento del
brazo involuntario con fuerza que puede conducir a una lesión traumática (hiperextensión).
Experimento 4:
Para esta prueba se propone variar las coordenadas articulares q2, q3 y q4 de (-180º a 180º),
responsables del movimiento de abducción/aducción, flexión/extensión del hombro y rotación
interna/externa del antebrazo respectivamente, para los valores de los individuos 2 y 6.
57
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Experimento 5:
Para este experimento se emplean los datos de los individuos 1 y 6, se propone variar las
coordenadas articulares q3, q5 y q6 de (-180º a 180º), mientras las demás coordenadas
articulares se quedan en posición de reposo.
58
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Figura 6.11. Movimientos Terapéuticos de flexión/extensión del hombro, flexión/extensión del codo y
supinación/pronación.
a b
Figura 6.12. Movimientos Terapéuticos de flexión/extensión del hombro, flexión/extensión del codo y
supinación/pronación, para: a) individuo 1 y para el b) individuo 6.
Donde q3, q5 y q6 son las responsables del movimiento de flexión/extensión del hombro,
flexión/extensión del codo y supinación/pronación. Obtenemos una respuesta muy diferente
al variar bastante el tamaño de los individuos, encontramos un índice de manipulabilidad
mayor en el primer caso cuando el individuo es mucho menor pero es necesario considerar
que los topes físicos para cada articulación deben ser respetados y en realidad tendríamos
la mitad de una circunferencia como área de trabajo, al simular un ejercicio de rehabilitación
como el de la Figura 6.11.
Experimento 6:
Se cogieron los datos de los individuos 2 y 5, se varió las coordenadas articulares q4 y q5
referentes a la rotación interna/externa del antebrazo y a la flexión/extensión del codo
respectivamente y se obtuvieron los siguientes resultados:
59
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a b
Figura 6.13. Movimientos de rotación interna/externa del antebrazo y a la flexión/extensión del codo,
para: a) individuo 2 y para el b) individuo 5.
En este caso tenemos puntos singulares muy visibles, por ejemplo, al momento de rotar q4 entre
(-30º y 30º) cuando la articulación llega a su límite inferior (-30º) no permite que la coordenada
articular q5 sea mayor de 0º a causa del choque del brazo con el tronco en el plano trasversal.
Experimento 7:
Empleamos al individuo3, se varían las articulaciones q1 y q4 mientras que las otras están en la
posición de reposo. Y ahora al mismo individuo se le se varían las articulaciones q 1 y q5
mientras que las otras están en la posición de reposo.
60
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
a b
Figura 6.15. Individuo 3, al variar: a) q1 y q4, b) q1 y q5.
Para modificar el valor de la primera coordenada articular se consideró en el ajuste inicial del
exoesqueleto a partir de la altura del individuo, pero en este caso particular se sugiere cambiar
el parámetro de q1 en relación a un movimiento reducido de elevación del hombro. A diferencia
de las gráficas obtenidas anteriormente, encontramos una cierta variación en el eje ‘z’ a razón
del cambio en la primera articulación traslacional.
Entre los inconvenientes existentes, durante la realización del código, para evaluar diferentes
pruebas, tenemos: un retraso significativo en el cálculo del índice de manipulabilidad del efector
final de la órtesis al modificar un gran número de los parámetros de entrada del sistema. Al
intentar variar más de tres coordenadas articulares el programa es un poco ineficiente. Otro
inconveniente es que al ejecutar iteraciones entre el rango de -90º y los 90º no se aprecia
detalladamente el espacio de trabajo real del exoesqueleto, por lo que en algunos casos se
evaluaron las iteraciones entre -180º y 180º y se fue calculando paso a paso las singularidades
existentes.
61
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
En la Tabla 6.3, se pueden comparar los tipos de movimientos con sus respectivos rangos de
acción para cada exoesqueleto de miembro superior. Siendo el primero el prototipo estudiado
en el presente trabajo, desarrollado en la empresa Aura Innovative Robotics, mientras que los
otros tres son órtesis comerciales existentes actualmente en el mercado.
62
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Por tanto, luego de realizar una serie de experimentos en base a los cálculos matemáticos
previamente realizados, se tiene que el nuevo prototipo tiene un amplio espacio de tarea que va
a brindar un buen nivel en cuanto a calidad de trabajo tanto al paciente como al terapeuta. Lo
que da pie, para continuar el estudio de mejoras que optimicen este innovador prototipo de
exoesqueleto para rehabilitación de miembro superior que pueda abarcar a un mayor número
de personas según su edad y tamaño.
63
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Capítulo 7
Conclusiones y Trabajos Futuros
En principio, el presente trabajo parte del estudio de conceptos clínicos relacionados con los
tipos de lesiones del plexo braquial, con el fin de ubicar e identificar el problema; para así,
adecuar los movimientos de rehabilitación que deben realizar los pacientes con ayuda de un
especialista. También se ha podido reflexionar acerca de la importancia que se le debe dar a
mejorar la calidad de vida de las personas, lo que hace indispensable realizar un análisis
completo de cuál es el espacio de trabajo de un nuevo e innovador exoesqueleto para
rehabilitación de miembro superior que se acomode a cada paciente.
Una órtesis debe aportar al paciente y al terapeuta una ayuda significativa en cuando a la
capacidad de manipulabilidad que puede llegar a alcanzar el mismo, es por ello que partiendo
del estudio detallado de una serie de exoesqueletos de miembro superior existentes en el
mercado se crean los cimientos para documentar un nuevo prototipo creado en la empresa Aura
Innovative Robotics.
Antes que nada, lo primordial era detallar el desempeño cinemático del exoesqueleto para
obtener una descripción analítica del movimiento espacial en función del tiempo, la relación
entre la posición y la orientación del efector final y los valores de las coordenadas articulares.
Luego, en MATLAB se comprobaron los cálculos matemáticos para posteriormente poder
definir el espacio de trabajo de nuestro exoesqueleto de 6 grados de libertad.
Es por eso que, al conocer el concepto de manipulabilidad como la capacidad que tiene el
extremo del robot para moverse y aplicar ciertas velocidades y/o fuerzas en distintas
direcciones, era necesario acotar el espacio de las entradas (espacio articular) al evaluar los
posibles movimientos de la órtesis en base a las limitaciones del cuerpo humano. En la mayoría
de los casos, al realizar un barrido tridimensional de al menos dos de las coordenadas articulares
del sistema nos encontramos con un volumen en forma de elipse o circunferencia que representa
la manipulabilidad de nuestro robot. Y a partir de una serie de experimentos se verifican los
rangos de movimientos previamente fijados, tomando en cuenta los puntos singulares; que en
su mayoría están conformados por las configuraciones donde el propio brazo humano no puede
alcanzar.
64
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Una forma de corroborar los puntos singulares es al anularse la matriz Jacobiana (determinante
de | J |), el índice de manipulabilidad de Yoshikawa (obtenido gracias a la ayuda de la Toolbok
de Peter Corke a modo de optimizar el tiempo de cálculo) colorea de azul el exoesqueleto de
miembro superior cuando se le hace imposible alcanzar libremente cualquier velocidad del
punto de referencia (valor muy pequeño del índice de Yoshikawa) y colorea de amarillo el área
alcanzable por el extremo del robot.
También, al comparar la nueva órtesis con los exoesqueletos comerciales se observa una gran
similitud en relación a los rangos de movimiento de las coordenadas articulares. En conclusión,
se aprecia una variación del espacio de trabajo según el tamaño del paciente porque, como se
observó a lo largo de las pruebas la amplitud del área de operación del entorno va a aumentar o
disminuir en consecuencia de este factor. De la misma manera, es necesario seguir los límites
fijados para cada articulación y considerar los puntos singulares para no perder grados de
libertad.
Al ser este nuevo prototipo algo original, desconocido y moderno para el mercado, se han
considerado muchas mejoras en cuando a materiales, dispositivos tecnológicos, controladores,
entre otros, que en conjunto con el aporte del presente trabajo optimizan la etapa de desarrollo
e implementación. Y de este modo, alcanzar una posición en el mercado actual que optimice el
proceso de rehabilitación para pacientes con lesiones del plexo braquial y surja para mejoras
futuras.
65
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Bibliografía
[1] Chang, C. Y., Wu, Y. T., Chen, L. C., Chan, R. C., Chang, S. T., & Chiang, S. L. (2015).
Massage-induced brachial plexus injury. Physical therapy, 95(1), 109-116.
[2] Bonham, C., & Greaves, I. (2011). Brachial plexus injuries. Trauma, 13(4), 353-363.
[3] García-López A, y D. Peréa Tortosa. Lesiones del plexo braquial. Manual Secot de cirugía
ortopédica y traumatología. Segunda edición. Editorial Médica Panamericana, cap. 69, pp. 701-
710. Madrid, 2009.
[4] Sutcliffe, T. Brachial plexus injury in the newborn. Neo Reviews: An official Journal of the
American Academy of Pediatrics, Vol.8, No.6, pp. 239-246. Junio, 2007.
[5] Thatte, M. R., Babhulkar, S., & Hiremath, A. (2013). Brachial plexus injury in adults:
Diagnosis and surgical treatment strategies. Annals of Indian Academy of Neurology, 16(1), 26.
[8] López, A. García; Tortosa, D. Perea. Lesiones del plexo braquial. En: Manual de cirugía
ortopédica y traumatología. Tomo 2. Sección XI: Miembro superior. Hombro. Sociedad
Española de Cirugía Ortopédica y Traumatología. 2o ed., Editorial Médica Paramericana S.A.,
2010, vol. 2, p. 701-710. ISBN Tomo 2: 978-84-9835- 324-2. ISBN Obra completa: 978-84-
9835-325-9
[9] Formato de referencia electrónica recomendado por Universidad Complutense de Madrid
(2012, 05 de octubre), [en línea]. Madrid, MADRID: Universidad Complutense de Madrid.
Disponible en: https://www.ucm.es/data/cont/docs/420-2014-02-26-13%20Patologia%20
plexo%20braquial.pdf [2016, 20 de mayo].
66
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
[13] Dogliotti, A.A. Conceptos actuales en la parálisis braquial perinatal. Parte 1: etapa
temprana. Archivos Argentinos de Pediatría, 2011, vol. 109, no 4, p. 347-353.
[14] Goubier, T. Technique of the double nerve transfer to recover elbow exion in C5, C6 or
C5 to C7 brachial plexus palsy. Techniques in Hand & Upper Extremity Surgery, Vol.11,
pp.15.17, 2007.
[15] Terzis, K. Supraescapular nerve reconstruction in 118 cases of adult post-traumatic
brachial plexus. Plastic and Reconstructive Surgery, Vol.117, pp.613-629, 2006.
[16] Reinkensmeyer, D. J. (2009). Robotic assistance for upper extremity training after stroke.
Stud Health Technol Inform, 145, 25-39.
[17] Riener, R., Nef, T., & Colombo, G. (2005). Robot-aided neurorehabilitation of the upper
extremities. Medical and Biological Engineering and Computing, 43 (1), 2-10.
[18] Kang, H. B., & Wang, J. H. (2013). Adaptive control of 5 DOF upper-limb exoskeleton
robot with improved safety. ISA transactions, 52(6), 844-852.
[20] Yu, W., & Rosen, J. (2013). Neural PID control of robot manipulators with application to
an upper limb exoskeleton. Cybernetics, IEEE Transactions on, 43(2), 673-684.
[21] Nef, T., Mihelj, M., Kiefer, G., Perndl, C., Muller, R., & Riener, R. (2007, June). ARMin-
Exoskeleton for arm therapy in stroke patients. In 2007 IEEE 10th International Conference on
Rehabilitation Robotics (pp. 68-74). IEEE.
[22] Nef, T., Mihelj, M., & Riener, R. (2007). ARMin: a robot for patient-cooperative arm
therapy. Medical & biological engineering & computing, 45(9), 887-900.
[23] Formato de referencia electrónica recomendado por The Catholic University Of America
(2010, 15 diciembre), [en línea]. Washington, DC: The Catholic University of. Disponible en:
http://cabrr.cua.edu/devicegallery.cfm [2016, 30 de julio].
[24] Hocoma, C. (2014). Armeo Therapy Concept. Catálogo de productos, [en línea].
Disponible en: http://www.hocoma.com/fileadmin/user/Dokumente/Armeo/bro_Armeo_
Therapy_Concept_140226_en.pdf [2016, 11 de Junio].
[25] Hocoma, C. (2012). Armeo Therapy Concept. Catálogo de productos, 1-2.
[26] ORMED GmbH & Co.KG (2007). ASTROMOT-S3 user manual, [en línea]. Disponible
en: http://storage.djoglobal.eu/en_US/Documents/Production_documents/ARTROMOT-S3_
user_manual.pdf [2016, 27 de Julio].
67
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
[29] García Montaño, J. Diseño mecánico de un robot serial para rehabilitación. Trabajo de
Fin de Grado en Ingeniería Mecánica. Universidad Politécnica de Madrid (UPM), Escuela
Técnica Superior de Ingeniería y Diseño Industrial. 2016.
[30] Monge Chamorro, L. J. Modelado y Simulación de Servomotores para Rehabilitación.
Análisis de la influencia en el movimiento biomecánico. Trabajo de Fin de Grado en Ingeniería
Electrónica Industrial y Automática. Universidad Politécnica de Madrid (UPM), Escuela
Técnica Superior de Ingeniería y Diseño Industrial. Curso 2014-2015.
[31] Dynamixel RX-64 Product Manual. Robotis Inc.
[34] Barrientos, A., & Barrientos, A. (2007). Fundamentos de robótica (No. 681.5 629.892). e-
libro, Corp..
[35] Barrientos, A., Álvarez, M., Hernández, J. D., del Cerro, J., & Rossi, C. (2012). Modelado
de Caden as Cinemáticas mediante Matrices de Desplazamiento. Una alternativa al método de
Denavit-Hartenberg. Revista Iberoamericana de Automática e Informática Industrial
RIAI, 9(4), 371-382.
[36] Corke, P. (2011). Robotics, vision and control: fundamental algorithms in MATLAB (Vol.
73). Springer.
[37] Moreno, H. A., Saltaren, R., Carrera, I., Puglisi, L., & Aracil, R. (2012). Índices de
Desempeño de Robots Manipuladores: una revisión del Estado del Arte. Revista
Iberoamericana de Automática e Informática Industrial RIAI, 9(2), 111-122.
[38] García, J. A. P., Cóbix, M. N., & García, J. J. P. Análisis del Desempeño de un robot en
una Tarea de Transferencia de Piezas.
68
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
[40] Corke, P. (2011). Robotics, vision and control: fundamental algorithms in MATLAB (Vol.
73). Springer.
69
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
Anexos
70
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
71
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
% --- Outputs from this function are returned to the command line.
function varargout = gui_exo_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
set(hObject,'Min', 0);
set(hObject,'Max', 20);
72
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
set(hObject,'Min', -pi);
set(hObject,'Max', pi);
set(hObject,'Min', -pi);
set(hObject,'Max', pi);
73
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
set(hObject,'Min', -pi);
set(hObject,'Max', pi);
set(hObject,'Min', -pi);
set(hObject,'Max', pi);
74
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
set(hObject,'Min', -pi);
set(hObject,'Max', pi);
L0= links(1,1);
L1= links(2,1);
L2= links(3,1);
L3= links(4,1);
75
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
% Denavit Hartenberg:
theta = [pi/2; q2+pi/2; q3; q4; q5; q6];
d = [L0+q1; L1; 0; L2; 0; L3];
a = [0; -10; 10; -10; 10; -10];
alpha = [pi/2; -pi/2; pi/2; -pi/2; pi/2; 0];
r_p = [1; 0; 0; 0; 0; 0];
DH = [theta, d, a, alpha];
% function T=MCD_MD(kk,pp,T0,tipo,n)
kk = [0 0 1; 1 0 0; 0 -1 0; 1 0 0; 0 -1 0; 1 0 0];
pp = [0 0 0; 0 0 20; 22 0 10; 0 0 20; 52 0 10; 0 0 20];
% pp = [0 0 0; 0 0 L0; L1 0 (L0-a); 0 0 L0; (L1+L2) 0 (L0-a); 0 0 L0];
T0 = [0 0 1 76; 0 -1 0 0; 1 0 0 10; 0 0 0 1];
% T0 = [0 0 1 (L1+L2+L3); 0 -1 0 0; 1 0 0 (L0-a); 0 0 0 1];
tipo = [1 0 0 0 0 0]; % prismatico= 1; rotacional = 0
n = 6; %GDL
qq = q;
T = T0;
76
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
I = eye(3);
solucion=1;
for i=1:1:n
k=kk(i,:)'; p=pp(i,:)'; qn=qq(i);
if tipo(i)==1;
matrices(i).D1=[I,k*qn;[0 0 0 1]]
else
K=[0 -k(3) k(2);k(3) 0 -k(1);-k(2) k(1) 0];
R=I*cos(qn)+k*k'*(1-cos(qn))+K*sin(qn);
matrices(i).D2=[R,(I-R)*p;[0 0 0 1]]
end
if(i==1)
solucion=matrices(i).D1;
else
solucion=solucion*matrices(i).D2;
end
end
T = solucion*T0; % Transformada Homogénea
jn = R.jacobn(qq');
detJ = abs(det(jn)); % Determinante para evaluar los puntos singulares
set(handles.text18, 'String', detJ);
set(handles.uitable4,'data', Jn1);
set(handles.uitable5,'data', Jn2);
% Gráfica Jv Lineales
axes(handles.axes2)
cla(handles.axes2)
color = ['m' 'b' 'g' 'k' 'y' 'r' 'c'];
for (i = 1 : n)
xx = norm([Jn1(1:3,i)]);
plot3([0,Jn1(1,i)]./xx,[0,Jn1(2,i)]./xx,[0,Jn1(3,i)./xx],color(i)); hold on
referencia(i) = {['Jv', num2str(i)]};
end
77
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
% Gráfica Jw Angulares
axes(handles.axes3)
cla(handles.axes3)
color = ['m' 'b' 'g' 'k' 'y' 'r' 'c'];
for (i = 1 : n)
xx = norm([Jn2(1:3,i)]);
plot3([0,Jn2(1,i)]./xx,[0,Jn2(2,i)]./xx,[0,Jn2(3,i)./xx],color(i)); hold on
referencia(i) = {['Jw', num2str(i)]};
end
X2(k) = R.maniplty(vecmani2,'yoshikawa');
78
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
79
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
80
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
81
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
82
DEPARTAMENTO DE AUTOMÁTICA,
INGENIERÍA ELÉCTRICA Y ELECTRÓNICA E
INFORMÁTICA INDUSTRIAL
CENTRO DE AUTOMÁTICA Y
ROBÓTICA (CAR-UPM-CSIC)
83