Study of dynamic simulators with ADAMS and Matlab/Simulink for control systems

== Resumen ==
La co-simulación ADAMS-MATLAB se usa ampliamente en el desarrollo de sistemas complejos para controlar modelos de mecanismos en lazo cerrado. Una ventaja de ADAMS es que facilita la configuración de parámetros como la fricción. Sin embargo, se desconoce su modelo matemático de fricción. Existen alternativas como Simscape Multibody que facilitan la configuración de cualquier modelo de fricción. Este artículo presenta un caso de estudio de simulación dinámica usando el modelo matemático de un mecanismo para compararlo con ADAMS y Simscape Multibody. El modelo de espacio de estado utiliza los parámetros b y fc: los coeficientes de fricción viscosa y de Coulomb. Sin embargo, ADAMS usa solo un parámetro de fricción dinámica µ. Dado que ambos parámetros de fricción no son conocidos en ADAMS, se empleó un procedimiento de identificación paramétrica para determinar la equivalencia entre los coeficientes de fricción b y fc con respecto a µ. Para identificar el modelo de fricción de ADAMS, se propusieron simulaciones en lazo abierto y cerrado. Para comparar numéricamente los resultados de la simulación se utilizó la norma L2. La referencia fue el modelo resuelto con la función ODE45. 

'''Palabras clave''': ADAMS, identificación de la fricción, Simscape Multibody, simulación.

==Abstract==

ADAMS-MATLAB co-simulation is widely used in complex system development to control models of closed-loop mechanisms. An advantage of ADAMS is to facilitate the configuration of parameters such as friction. However, its mathematical friction model is unknown. There are alternatives such as Simscape Multibody that facilitate the configuration of any friction model. This paper presents a case study of dynamic simulation using the mathematical model of a mechanism in order to compare it with ADAMS and Simscape Multibody. The state-space model uses the parameters b and fc – viscous and Coulomb friction coefficients. However, ADAMS uses only one dynamic friction parameter µ. Since both friction parameters do not appear in ADAMS, a parametric identification procedure was used to determine the equivalency between the friction coefficients b and fc compared to µ. To identify the friction model of ADAMS, open-loop and closed-loop simulations were proposed. To compare the simulation results numerically, the L2 norm was used. The reference was the model solved with the ODE45 function.

'''Keywords''': ADAMS, friction identification, Simscape Multibody, simulation<span id="PointTmp"></span>

==1. Introducción==

La pandemia ocasionada por el COVID-19 no solo trajo un reto enorme al sistema educativo, sino que también a todos los investigadores y educadores que realizan investigación a distancia o que imparten materias a distancia relacionadas con el diseño en ingeniería [1]. La educación a distancia es complicada debido a que es difícil mantener la calidad e interés de los alumnos cuando se necesitan abordar prácticas de laboratorio [2]. Los simuladores basados en modelos virtuales como ADAMS y Simscape Multibody (MATLAB) son una interesante alternativa al desarrollo de prácticas en un laboratorio. Sin embargo, la falta de pericia y experiencia de los alumnos son barreras para usar estos programas. Para impartir una clase de diseño en ingeniería a distancia se requiere tener un conocimiento multidisciplinario para configurar correctamente los parámetros de la simulación ya que si no se hace bien se pueden obtener resultados distintos [3]. Para disminuir la complejidad de los simuladores se pueden utilizar laboratorios virtuales para realizar prácticas en los cuales los parámetros de simulación están precargados para facilitar el aprendizaje y reducir el margen de error [4].

Un ejemplo muy particular es cuando pueden ser utilizados para el diseño e implementación de sistemas de control interactivo en bancos de pruebas de motores de CD para demostrar los conceptos fundamentales en ingeniería mecatrónica [5] o también en simuladores de bola y plato para demostrar la dinámica del sistema bajo la acción de un sistema de control como parte de la materia de robótica [6]. Desafortunadamente, los laboratorios virtuales están limitados porque solamente resuelven la tarea para la que fueron creados. Además, para reproducir un experimento en simulación es necesario contar con el modelo matemático de la plataforma experimental y los valores numéricos de sus parámetros.

En comparación con un laboratorio virtual, el modelo virtual permite diseñar, adaptar y modificar los parámetros de simulación de acuerdo con las necesidades de cada práctica o investigación [7]. Eldirdiry et al. [8] estudiaron un modelo de un robot bípedo el cual puede simular el pie caído de una persona al caminar. Milanowski y Pilat [9] realizaron una simulación de un exoesqueleto basado en datos adquiridos de la marcha humana y Urrea et al. [10] simularon un modelo dinámico de un robot hexápodo. Similarmente, Beaber et al. [11], analizaron un hexápodo y Zhou et al. [12] analizaron la dinámica de un exoesqueleto para verificar su viabilidad antes de construirlo. Otros autores [13]–[18] también han reportado el uso de modelos virtuales en sistemas robóticos. Xia et al. [19] llevaron a cabo la simulación de un robot cuadrúpedo. Liu et al. [20] diseñaron, simularon y controlaron un mecanismo complejo de varias patas en un robot.

También en la literatura están reportadas las co-simulaciones de ADAMS con MATLAB/Simulink. Sosa-Méndez et al. [21] propusieron un análisis cinemático y dinámico de la plataforma Stewart-Gough. Parthasarathy et al. [22] presentaron el modelo dinámico y la simulación de un robot Mitsubishi RM-501. Nair y Ezhilarasi [23], analizaron el desempeño de un controlador de modos deslizantes en el seguimiento del movimiento de un exoesqueleto para miembro inferior. Liao et al. [24] modelaron y simularon un exoesqueleto para rehabilitación de miembros inferiores. Pandey et al. [25] simularon el movimiento de nado de una rana. La simulación fue usada para calcular las fuerzas propulsivas en las patas de la rana. Brezina et al. [26] presentaron una eficiente técnica para el desarrollo de sistemas mecatrónicos basada en la metodología del modelo V. Otros autores [27]–[29] también han reportado simulaciones de mecanismos avanzados, automatización de procesos y análisis dinámicos.

Aunque existe una gran cantidad de trabajos de investigación relacionados con simulaciones y co-simulaciones de modelos virtuales, todavía no existe un estudio que compare y analice ADAMS y Simscape Multibody a detalle, teniendo como referencia un modelo matemático para determinar el nivel de error de cada software de manera independiente. En la búsqueda del estado del arte no se encontró información similar o experimental que aborde este problema por esta razón este es el primer trabajo que aborda esta problemática. El principal objetivo de este trabajo es presentar un caso de estudio de simulación dinámica usando el modelo matemático de un mecanismo para compararlo con ADAMS y Simscape Multibody, debido a que estos paquetes de software son una referencia en el estudio de sistemas multicuerpo. La contribución de esta investigación es el análisis de retos y oportunidades para el desarrollo de simuladores de mecanismos con ADAMS y Simscape Multibody, tomando un caso de estudio de un modelo virtual de un grado de libertad (GDL), el cual puede ser una referencia para futuras investigaciones en el área de control de mecanismos. Esto es importante porque la mayoría de los investigadores y educadores utilizan Matlab para impartir clases o desarrollar investigación. La presentación de este artículo se estructura de la siguiente manera. La sección 2 contiene la descripción de los materiales y métodos usados en este estudio. La sección 3 discute sus resultados en lazo abierto y cerrado. Finalmente, la sección 4 concluye el artículo.

==2. Materiales y métodos==

La [[#img-1|Figura 1]] muestra la metodología. La sección 2.1 muestra el desarrollo del modelo CAD usando SolidWorks, especificando sus materiales y dimensiones. La sección 2.2 se divide en tres partes: 1) obtener las ecuaciones de movimiento de Euler-Lagrange para plantearlas como una ecuación diferencial ordinaria (ODE, por sus siglas en inglés), calculando los parámetros de distribución de masa en SolidWorks, 2) exportar el modelo CAD obtenido en SolidWorks a ADAMS, configurando la fricción de la articulación en ADAMS, así como sus entradas y salidas. Después se exporta la planta no lineal de ADAMS a Simulink, y 3) exportar el modelo de SolidWorks a Simulink usando el Simscape Multibody Link. En la sección 2.3 se realiza la identificación paramétrica de la fricción. En la sección 2.4 se describe el procedimiento para validar el modelo identificado, mientras que en la sección 2.5 se describe la forma de cuantificar los resultados. En este estudio se propone un péndulo rígido de un GDL sin elasticidad en la articulación. Por su modelo dinámico se considera que es un caso particular de un manipulador serial.

<div id='img-1'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Metodología.png|centre|400px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 1'''. Metodología 
|}

===2.1 Modelo CAD en SolidWorks===

En ingeniería es esencial el uso de programas para el diseño y simulación de prototipos antes de su implementación. En el caso de los mecanismos no lineales, como los manipuladores, es de gran interés contar con el modelo dinámico para el diseño y evaluación de esquemas de control. Aunque existen varias herramientas que permiten obtener el modelo dinámico de un manipulador serial a partir de la descripción cinemática basándose en el algoritmo de Denavit-Hartenberg [30], es necesario que el usuario realice la medición de los parámetros a partir de un modelo CAD. Para este estudio, el modelo propuesto es un péndulo que comprende la base y el eslabón. El material seleccionado fue aluminio 6061 por que en la industria es muy común utilizar este material para el prototipado de mecanismos en sistemas de control.

===2.2 Desarrollo de modelos para simulación===

La [[#img-2|Figura 2]] muestra el diagrama cinemático del péndulo. El modelo basado en las ecuaciones de movimiento Euler-Lagrange se presenta en la ec.(1); la [[#tab-1|Tabla 1]] describe su nomenclatura

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>{\theta }_{1}\ddot{q}+{\theta }_{2}\mathrm{sin}\,\left( q\right) +{\theta }_{3}\mathrm{cos}\,\left( q\right) +{\theta }_{4}\overset{\cdot}{q}+{\theta }_{5}sgn\left( \overset{\cdot}{q}\right) =\tau </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}

donde los parámetros agrupados <math display="inline">{\theta }_{i},\, i=</math><math>1,\ldots ,5</math> son los siguientes:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\begin{array}{l} \theta_{1} =m\left( {r}_{x}^{2}+{r}_{y}^{2}\right) +{I}_{zzc}\\ \theta_{2}=mg{r}_{x}\\ \theta_{3}=mg{r}_{y}\\ \theta_{4}=b\\ \theta_{5}={f}_{c}
\end{array}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(2)
|}

mientras que <math display="inline">sgn(\alpha )</math> es la función signo, dada por:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>sgn\left( \alpha \right) =\left\{ \begin{matrix}-1&si\, \alpha <0\\0&si\, \alpha =0\\1&si\, \alpha >0\end{matrix}\right. \quad</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

<div id='img-2'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[File:Diagrama Péndulo actualizado.png|centre|250px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 2'''. Diagrama cinemático del péndulo
|}


<div class="center" style="font-size: 75%;">'''Tabla 1.''' Nomenclatura</div>

<div id='tab-1'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
!'''Símbolo''' !! style="text-align: left;"| '''Descripción'''
|-
|  style="text-align: center;"|<math>{x}_{0}{y}_{0}{z}_{0}</math>
|  style="text-align: left;"|Sistema fijo o inercial
|-
|  style="text-align: center;"|<math>{x}_{1}{y}_{1}{z}_{1}</math>
|  style="text-align: left;"|Sistema unido al eslabón sobre el eje de rotación
|-
|  style="text-align: center;"|<math>{x}_{c}{y}_{c}{z}_{c}</math>
|  style="text-align: left;"|Sistema unido al centro de masa del cuerpo rígido y alineado con el sistema  <math display="inline">{x}_{1}{y}_{1}{z}_{1}</math>
|-
|  style="text-align: center;"|<math>g</math>
|  style="text-align: left;"|Constante gravitacional de la tierra
|-
|  style="text-align: center;"|<math>q,\overset{\cdot}{q},\ddot{q}</math>
|  style="text-align: left;"|Variables de posición, velocidad y aceleración articular
|-
|  style="text-align: center;"|<math>\tau</math> 
|  style="text-align: left;"|Es el par de torsión inyectado por el actuador
|-
|  style="text-align: center;"|<math>m</math>
|  style="text-align: left;"|Masa del cuerpo rígido (eslabón)
|-
|  style="text-align: center;"|<math>{\left[ {r}_{x},{r}_{y},{r}_{z}\right] }^{T}</math>
|  style="text-align: left;"|Vector de posición del centro de masa del cuerpo rígido con respecto al sistema  <math display="inline">{x}_{1}{y}_{1}{z}_{1}</math>
|-
|  style="text-align: center;"|<math>{I}_{zzc}</math>
|  style="text-align: left;"|Momento de inercia en el eje z, calculado con respecto al sistema  <math display="inline">{x}_{c}{y}_{c}{z}_{c}</math>
|-
|  style="text-align: center;"|<math>b</math>
|  style="text-align: left;"|Coeficiente de fricción viscosa
|-
|  style="text-align: center;"|<math>{f}_{c}</math>
|  style="text-align: left;"|Coeficiente de fricción de Coulomb
|}


La ecuación diferencial ordinaria (4) describe el modelo en espacio de estados, empleando el vector de estados <math display="inline">{\left[ q,\overset{\cdot}{q}\right] }^{T}</math>. La solución de esta ecuación puede resolverse numéricamente con Matlab usando el método de Runge-Kutta 4/5 con la función ODE45. La intención de un laboratorio virtual debería ser llevar a experimental lo realizado en simulación, por esta razón se seleccionó un método de paso variable debido a que presenta una solución que se acerca a un sistema analógico

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>\frac{d}{dt}\left[ \begin{matrix}q\\\overset{\cdot}{q}\end{matrix}\right] =\left[ \begin{matrix}\overset{\cdot}{q}\\{\theta }_{1}^{-1}\left[ \tau -{\theta }_{2}{S}_{q}-{\theta }_{3}{C}_{q}-{\theta }_{4}\overset{\cdot}{q}-{\theta }_{5}sgn(\overset{\cdot}{q})\right] \end{matrix}\right] </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

donde <math display="inline">{S}_{q}=\mathrm{sin}\,(q)</math>  y <math display="inline">{C}_{q}=\mathrm{cos}\,(q)</math>.

La [[#img-3|Figura 3]] muestra la captura de pantalla de las propiedades físicas del eslabón. Los parámetros de distribución de masa se indican con recuadros de color, los cuales se obtuvieron con respecto al sistema de coordenadas unido al eslabón <math display="inline">S1</math>, el cual corresponde al sistema <math display="inline">{x}_{1}{y}_{1}{z}_{1}</math>. El color rojo indica la masa, el azul indica el centro de masa, el verde indica el momento de inercia en <math display="inline">{z}_{1}</math> y  el magenta indica los momentos principales de inercia. La [[#tab-2|Tabla 2]] muestra los parámetros de simulación que se emplearon en este estudio.

<div id='img-3'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[File:Parámetros Solidworks.png|centre|500px]] 
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 3'''. Parámetros de distribución de masa en SolidWorks 
|}


Con el modelo CAD desarrollado en SolidWorks se realizó un estudio de movimiento con la posición inicial del péndulo en <math display="inline">{q}_{0}\, =\, 90\mbox{°}</math>. Una vez realizado esto, se exportó a ADAMS.

<div class="center" style="font-size: 75%;">'''Tabla 2'''.Parámetros de distribución de masa</div>

<div id='tab-2'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! Parámetro !! Valor !! Unidad
|-
|  style="text-align: center;vertical-align: top;"|<math>m</math>
|  style="text-align: center;vertical-align: top;"|10.909762
|  style="text-align: center;vertical-align: top;"|kg
|-
|  style="text-align: center;vertical-align: top;"|<math>{r}_{x}</math>
|  style="text-align: center;vertical-align: top;"|0.243926 
|  style="text-align: center;vertical-align: top;"|m
|-
|  style="text-align: center;vertical-align: top;"|<math>{r}_{y}</math>
|  style="text-align: center;vertical-align: top;"|0.000000
|  style="text-align: center;vertical-align: top;"|m
|-
|  style="text-align: center;vertical-align: top;"|<math>{I}_{zzc}</math>
|  style="text-align: center;vertical-align: top;"|0.437764
|  style="text-align: center;vertical-align: top;"|kg·m<sup>2</sup>
|}


La [[#img-4|Figura 4]] muestra los parámetros de distribución de masa en ADAMS. La [[#img-4|Figura 4]](a) muestra la posición del centro de masa con respecto a <math display="inline">S1</math>. La [[#img-4|Figura 4]](b) muestra el valor de la masa y los momentos principales de inercia.

<div id='img-4'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|[[File:Chavez-Olivares_et_al_2022a_6087_300px-Parámetros_ADAMS_2(a).png|350px]]
|-
|  style="text-align: center;font-size: 75%;"|(a) 
|-
|style="padding:10px;"|[[File:Chavez-Olivares_et_al_2022a_7118_300px-Parámetros_ADAMS_2(b).png|350px]]
|-
|  style="text-align: center;font-size: 75%;"|(b)
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 4'''. Parámetros de distribución de masa en ADAMS
|}


<!--[[File:Parámetros ADAMS 2.png|centre|300px]]-->

=== 2.3 Configuración de la fricción en la articulación ===

Debido a que la ecuación de movimiento (1) considera únicamente fricción dinámica, se seleccionó esta opción en el cuadro de diálogo de articulación en ADAMS. El parámetro de fricción seleccionado para la articulación del péndulo fue de <math display="inline">\mu =</math><math>0.03</math> Nm s/rad, la selección se realizó de forma arbitraria para tener un punto de partida, ya que el modelo de fricción de la articulación es desconocido para el usuario ([[#img-5|Figura 5]]). Sin embargo, el modelo en espacio de estados (4) utiliza los parámetros <math display="inline">b</math> y <math display="inline">{f}_{c}</math>. Debido a que no aparecen ambos parámetros de fricción en este cuadro de diálogo, se consideró inicialmente que el parámetro <math display="inline">\mu</math>  es equivalente al coeficiente de fricción viscosa <math display="inline">b</math>, mientras que el coeficiente de fricción de Coulomb <math display="inline">{f}_{c}</math> se consideró nulo.

<div id='img-5'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Parámetro fricción ADAMS.png|centre|350px]] 
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 5'''. Configuración de la fricción en ADAMS
|}


Para obtener la co-simulación ADAMS-Simulink fue necesario configurar entradas y salidas. Como única entrada se tiene el torque aplicado en la articulación <math display="inline">\tau</math>  y como salidas la posición y velocidad articulares, <math display="inline">q</math> y <math display="inline">\overset{\cdot}{q}</math>, respectivamente. Finalmente, con el complemento ''ADAMS controls'' se exportó el modelo como una planta no lineal. Los estados del sistema no lineal fueron resueltos en ADAMS, mientras que en Simulink se realizó la solución del algoritmo de control. La [[#img-6|Figura 6]](a) muestra el diagrama de bloques en lazo abierto. El modelo de SolidWorks se exportó a Simscape Multibody (Simulink) usando el complemento Simscape Multibody Link.

El proceso de exportación de SolidWorks a Simscape Multibody se realizó en dos pasos. El primero fue generar archivos intermedios y el segundo convertirlos al modelo de Simulink. Los archivos generados fueron de tres tipos: a) los sólidos 3D (STL o STEP), b) un archivo de extensión *.m que contiene todos los parámetros cinemáticos y dinámicos del mecanismo (masas, momentos de inercia, desplazamientos y rotaciones), y c) un archivo XML que contiene los datos necesarios para formar el ensamble en Simulink. El segundo paso fue crear el modelo de Simulink con un comando de importación en Matlab. La [[#img-6|Figura 6]](b) muestra el modelo en lazo abierto de Simscape Multibody.

<div id='img-6'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:ADAMS lazo abierto.png|centre|500px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 6'''. Modelo en lazo abierto 
|}


<!--[[File:ADAMS lazo abierto.png|centre|thumb]]<span style="text-align: center; font-size: 75%;">'''Figura 6.'''  Modelo en lazo abierto.</span> </div>-->

=== 2.4 Identificación del modelo de fricción de la articulación ===

Los tres modelos se excitaron con el mismo torque de entrada para comparar sus respuestas ([[#img-7|Figura 7]]), con <math display="inline">b\, = \mu</math>, considerando:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
|<math>\left[ \begin{matrix}{q}_{0}\\{\overset{\cdot}{q}}_{0}\end{matrix}\right] = \left[ \begin{matrix}90\mbox{°}\\0\mbox{°}/s\end{matrix}\right] </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

<div id='img-7'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Comparación lazo abierto.png|centre|376px]] 
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 7'''. Esquema de simulación en lazo abierto
|}


La [[#img-8|Figura 8]] muestra los resultados de simulación en lazo abierto considerando sólo fricción viscosa con el modelo de fricción del cuadro de diálogo de la articulación.

<div id='img-8'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto considerando sólo fricción viscosa.png|centre|376px]] 
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 8'''. Simulación en lazo abierto considerando sólo fricción viscosa. <br> Torques de entrada usados : (a) <math display="inline">\tau = 0</math> N·m. (b) <math display="inline">\tau =5 \sin\left( \frac{2\pi }{3}\cdot t\right)</math>N·m
|}


Los resultados obtenidos de ADAMS no coincidieron con los de Simscape Multibody y ODE45. Tomando en cuenta que los parámetros de distribución de masa son iguales, es posible deducir que el modelo de fricción no fue propuesto correctamente ya que como se mencionó anteriormente este modelo es desconocido para el usuario. Por esta razón, se exploró la posibilidad de usar el modelo de fricción tradicional y caracterizar este modelo. Para encontrar la relación de <math display="inline">\mu</math>  con los coeficientes de la fricción tradicional <math display="inline">b</math> y <math display="inline">{f}_{c}</math>, se realizó una identificación paramétrica en dos pasos: 1) se aplicó un par de excitación en lazo abierto al modelo ADAMS-Simulink y se registraron entradas y salidas como el par aplicado, las posiciones y velocidades articulares resultantes, es decir, <math display="inline">\tau</math> , <math display="inline">q</math> y <math display="inline">\overset{\cdot}{q}</math>, respectivamente, y 2) los valores obtenidos en el paso 1 se usaron en el modelo de estimación paramétrico desarrollado en Matlab/Simulink, utilizando la herramienta ''Parameter estimation'' del toolbox ''Simulink Design Optimization''. La [[#img-9|Figura 9]] muestra el esquema de identificación paramétrica.

<div id='img-9'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Identificacion.png|centre|446x446px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 9'''. Esquema de identificación paramétrica. (a) Modelo ADAMS-Simulink. (b) Modelo de estimación Simulink
|}


Para realizar la identificación paramétrica, se propuso el torque de excitación en lazo abierto mostrado en la ec.(6)

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>{\tau }_{id}=5 \sin\left( \frac{2\pi }{3}\cdot t\right) + 3\sin\left( 2\pi \cdot t+45\mbox{°}\right) \, N\cdot m </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

El resultado obtenido fue:

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>b=0.278421\, N\cdot m\, s/rad {f}_{c}=2.978608\, N\cdot m</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

Tomando en cuenta que el valor configurado en ADAMS fue considerado como <math display="inline">\mu  = 0.03</math> N·m s/rad, se estableció la siguiente relación:

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>b=10\cdot \mu {f}_{c}=100\cdot \mu</math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

=== 2.5 Validación del modelo identificado ===

Para verificar si el modelo identificado es válido, entonces fue necesario aplicar además de la entrada usada en la identificación; otras entradas diferentes con una excitación más compleja para provocar más cambios en velocidad y aceleración en el péndulo con el objetivo de verificar que el seguimiento de trayectorias sea correcto. La validación contempló usar diferentes valores de <math display="inline">\mu</math>  considerando las mismas condiciones iniciales. Los resultados de simulación en lazo abierto se muestran en la sección 3.1, usando las ecs.(9)-(11)

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|
<math>{\tau }_{A1}=5\sin\left( \frac{2\pi }{3}\cdot t\right) +3\sin\left( 2\pi \cdot t+45\mbox{°}\right) </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}<div class="center" style="width: auto; margin-left: auto; margin-right: auto;"></div>

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\tau }_{A2}=8\sin\left( \frac{2\pi }{3}\cdot t\right) +5\sin\left( 2\pi \cdot t+45\mbox{°}\right) \, +3\sin\left( \pi \cdot t-55\mbox{°}\right) \, \,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\tau }_{A3}=10\sin\left( \frac{2\pi }{3}\cdot t\right) +6\sin\left( 2\pi \cdot t+45\mbox{°}\right) +4\sin\left( \pi \cdot t-55\mbox{°}\right) \, \,</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

Cuando se comparan dos o más sistemas, en este caso un péndulo excitado en lazo abierto usando la misma entrada, pero con ligeras variaciones en sus coeficientes de fricción, entonces la posición articular entre estos sistemas puede separarse a lo largo del tiempo. Esta separación puede ser más notable cuando el sistema termina sobrepasando su posición vertical inestable. Sin embargo, en lazo cerrado la variación de los parámetros es compensada por el controlador. Para complementar la validación, se propuso una comparación en lazo cerrado usando un controlador de posición <math>\tanh</math> ([[#img-10|Figura 10]]) [31]. El controlador usado fue el siguiente:

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math>{\tau }_{Ci}={K}_{pi}\tanh\left( {\alpha }_{i}\cdot \tilde{q}\right) -{K}_{vi}\tanh\left( {\beta }_{i}\cdot \overset{\cdot}{q}\right) +g(q) </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

donde <math display="inline">\tilde{q} = qd-q (t)</math> es el error de posición, <math display="inline">K_p > 0</math> es la ganancia de la acción proporcional, <math display="inline">K_v > 0</math> es la ganancia de la acción derivativa, <math display="inline">\alpha  > 0</math> es la constante de sensibilidad al error, <math display="inline">\beta  > 0</math> es la constante de sensibilidad a la velocidad, mientras que el torque de carga gravitacional es <math display="inline">g(q) = {\theta }_{2}\, \sin(q) + {\theta }_{3}\, \cos(q)</math>. El subíndice <math display="inline">i</math> se emplea para especificar el número de la simulación.

<div id='img-10'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Compara lazo cerrado.png|centre|550x550px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 10'''. Esquema de simulación en lazo cerrado
|}


La [[#tab-3|Tabla 3]] muestra las constantes de sintonía del controlador considerando que el error de posición articular se midió en radianes y la velocidad articular en rad/s. Las condiciones iniciales fueron las mismas que en lazo abierto. La sintonía se realizó de forma empírica en base a la experiencia de los autores.

<div class="center" style="font-size: 75%;">'''Tabla 3'''. Ganancias del controlador</div>

<div id='tab-1'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! <math>i</math> !! <math>{K}_{p}</math> !! <math>{K}_{v}</math> !! <math>\alpha</math> !!<math>\beta</math> 
|-style="text-align:center"
|  <math>1</math>
|  7
|  4
|  6
|  1
|-style="text-align:center"
|  <math>2</math>
|  15
|  4 
|  6
|  1
|-style="text-align:center"
|  <math>3</math>
|  15
|  4
|  6
|  1
|}


La trayectoria de referencia <math display="inline">{q}_{d}\, (t)</math> fue configurada para un controlador punto a punto como sigue:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math display="inline">{q}_{d}\, (t)=\left\{ \begin{array}{cl} 0\mbox{°} & 0\, s\, \leq t\leq 5\, s\\ -45\mbox{°} & 5\, s\, <t\leq 10s\\ 0\mbox{°}& 10\, s\, <t\leq 15s\\ 45\mbox{°}&\, 15\, s\, <t\leq 20s\end{array}\right. \quad</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

=== 2.6 Cuantificación de la diferencia ===
Para comparar los resultados de simulación de forma numérica se utilizó la norma <math display="inline">{L}_{2}</math>, aplicada a la diferencia entre las señales. La referencia fue el modelo matemático resuelto con la función ODE45.

''Error de posición en lazo abierto y cerrado'':

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math display="block">{\tilde{q}}_{S}={q}_{Simscape}-{q}_{ODE45}\quad </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math display="block">{\tilde{q}}_{A}={q}_{ADAMS}-{q}_{ODE45} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}

''Error de torque en lazo cerrado'':

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math display="block">{\tilde{\tau }}_{S}={\tau }_{Simscape}-{\tau }_{ODE45} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}
{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math display="block">{\tilde{\tau }}_{A}={\tau }_{ADAMS}-{\tau }_{ODE45} </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}


Sea <math display="inline">\tilde{x}(t)</math> la señal de diferencia para cuantificar, entonces

{| class="formulaSCP" style="width: 100%; text-align: center;"  
|-
|<math display="block">{L}_{2}\left\{ \tilde{x}\right\} =\sqrt{\frac{1}{T}}\int_{0}^{T}{\left\| \tilde{x}(t)\right\| }^{2}dt </math>
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}

El valor más pequeño de la norma <math display="inline">{L}_{2}</math> permite determinar el grado de aproximación entre los sistemas que se comparan.

==3. Resultados==

===3.1 Simulación en lazo abierto===

La Figura 11 muestra 20 segundos de simulación en lazo abierto usando la entrada <math display="inline">{\tau }_{A1}</math> con <math display="inline">\mu =0.03</math>. La [[#img-11|Figura 11]](a) muestra el torque de excitación que va desde -6.74 Nm hasta 6.74 Nm con dos componentes sinusoidales. La [[#img-11|Figura 11]](b) muestra la posición articular del modelo la cual parte de 90°.

<div id='img-11'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto usando la entrada TA1.png|centre|450px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 11'''. Simulación en lazo abierto usando la entrada  <math display="inline">{\tau }_{A1}</math> con <math display="inline">\mu =0.03</math> 
|}


La [[#img-12|Figura 12]] muestra 20 segundos de simulación en lazo abierto usando la entrada <math display="inline">{\tau }_{A2}</math> con <math display="inline">\mu =0.03</math>. La [[#img-12|Figura 12]](a) muestra el torque de excitación que va desde -13.28 Nm hasta 13.09 Nm con tres componentes sinusoidales. La [[#img-12|Figura 12]](b) muestra la posición articular del modelo la cual parte de 90°.

<div id='img-12'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto usando la entrada TA2.png|centre|339x339px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 12'''. Simulación en lazo abierto usando la entrada  <math display="inline">{\tau }_{A2}</math> con <math display="inline">\mu =0.03</math> 
|}


La [[#img-13|Figura 13]] muestra 20 segundos de simulación en lazo abierto usando la entrada <math display="inline">{\tau }_{A3}</math> con <math display="inline">\mu =</math><math>0.05</math>. La [[#img-13|Figura 13]](a) muestra el torque de excitación que va desde -16.55 Nm hasta 16.34 Nm con tres componentes sinusoidales. La [[#img-13|Figura 13]](b) muestra la posición articular del modelo la cual parte de 90°.

<div id='img-13'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto usando la entrada TA3.png|centre|339x339px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 13'''. Simulación en lazo abierto usando la entrada  <math display="inline">{\tau }_{A3}</math> con <math display="inline">\mu =0.05</math> 
|}


Los resultados de simulación en lazo abierto aplicando diferentes señales de excitación muestran un comportamiento similar de la trayectoria realizada por la posición articular.

=== 3.2 Simulación en lazo cerrado ===

La [[#img-14|Figura 14]] muestra 20 segundos de simulación en lazo cerrado usando la entrada <math display="inline">{\tau }_{C1}</math> con <math display="inline">\mu =0.01</math>. La [[#img-14|Figura 14]](a) muestra el torque de control. La [[#img-14|Figura 14]](b) muestra la posición angular. los resultados del sistema en lazo cerrado con el coeficiente <math display="inline">\mu  = 0.01</math>, es decir, un valor por debajo del usado para estimar los parámetros de fricción. En este caso se observa también que los resultados del torque de control y la posición son iguales en los tres modelos, por lo que la relación (8) se considera válida.

<div id='img-14'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto usando la entrada TC1.png|centre|330x330px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 14'''. Simulación en lazo cerrado usando la entrada  <math display="inline">{\tau }_{C1}</math> con <math display="inline">\mu =0.01</math> 
|}


La [[#img-15|Figura 15]] muestra los resultados del sistema configurado en lazo cerrado y el coeficiente <math display="inline">\mu  =0.03</math>, es decir, con el que se realizó la estimación de los parámetros de fricción. Como se puede observar, el controlador todavía puede sintonizarse mejor para evitar el sobretiro en el seguimiento de la referencia punto a punto, sin embargo, no es el propósito de este trabajo. Lo importante de este resultado es que las soluciones obtenidas en los tres modelos son prácticamente iguales, tanto en el torque de control <math display="inline">\tau</math>, como en la posición controlada <math>q</math>, por lo que se considera que la relación (8) es válida.

<div id='img-15'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"|  [[File:Simulación en lazo abierto usando la entrada TC2.png|centre|349x349px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 15'''. Simulación en lazo cerrado usando la entrada  <math display="inline">{\tau }_{C2}</math> con <math display="inline">\mu =0.03</math> 
|}


La [[#img-16|Figura 16]]  muestra los resultados de la simulación en lazo cerrado con el coeficiente <math display="inline">\mu  = 0.05</math>, es decir, un valor superior al usado para estimar los parámetros de fricción. En este caso se observa que las soluciones de los tres modelos para la posición controlada son prácticamente iguales, mientras que los torques de control presentan ligeras diferencias, apenas notables, por lo que la relación propuesta para los coeficientes de fricción se considera válida.

<div id='img-16'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[File:Simulación en lazo abierto usando la entrada TC3.png|centre|339x339px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 16'''. Simulación en lazo cerrado usando la entrada  <math display="inline">{\tau }_{C3}</math> con <math display="inline">\mu =0.05</math> 
|}

===3.3 Comparación numérica===

La [[#tab-4|Tabla 4]] muestra la diferencia entre simulaciones en lazo abierto usando la norma <math display="inline">{L}_{2}</math>. La diferencia de posición se obtuvo usando las ecuaciones (14) y (15). Como se observa, el modelo de Simscape Multibody presentó un error de posición menor en comparación con el modelo de ADAMS. La [[#tab-5|Tabla 5]] muestra la diferencia entre simulaciones en lazo cerrado usando la norma <math display="inline">{L}_{2}</math>.

<div class="center" style="font-size: 75%;">'''Tabla 4'''. Diferencia entre simulaciones en lazo abierto</div>

<div id='tab-4'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
!Entrada !! <math display="inline">{L}_{2}\left\{ {\tilde{q}}_{S}\right\}</math> [rad] !! <math display="inline">{L}_{2}\left\{ {\tilde{q}}_{A}\right\}</math> [rad]
|-style="text-align:center"
|  <math>{\tau }_{A1}</math>
|  0.0022
|  0.0163
|-style="text-align:center"
|  <math>{\tau }_{A2}</math>
|  0.0044
|  0.0187
|-style="text-align:center"
|  <math>{\tau }_{A3}</math>
|  >0.0036
|  0.0326
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<span style="text-align: center; font-size: 75%;">'''Tabla 5.'''</span> <span style="text-align: center; font-size: 75%;">Diferencia entre Simulaciones en Lazo Cerrado</span></div>

{| style="width: 69%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Entrada</span>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math display="inline">{L}_{2}\left\{ {\tilde{\tau }}_{S}\right\}</math> [N<span style="text-align: center; font-size: 75%;">·</span>m]
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math display="inline">{L}_{2}\left\{ {\tilde{\tau }}_{A}\right\}</math> [N·m]
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math display="inline">{L}_{2}\left\{ {\tilde{q}}_{S}\right\}</math> [rad]
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<math display="inline">{L}_{2}\left\{ {\tilde{q}}_{A}\right\}</math> [rad]
|-
|  style="text-align: center;vertical-align: top;"|<math>{\tau }_{C1}</math>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.2804</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.5814</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0043</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0143</span>
|-
|  style="text-align: center;vertical-align: top;"|<math>{\tau }_{C2}</math>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0357</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.1626</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0022</span>
|  style="text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0073</span>
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\tau }_{C3}</math>
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.2948</span>
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.7098</span>
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0047</span>
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.0208</span>
|}


Las diferencias de posición y torque se obtuvieron usando las ecuaciones (14)-(15) y (16)-(17), respectivamente. Como se observa, el modelo de Simscape Multibody presentó errores de posición y torque menores en comparación con el modelo de ADAMS. A partir de los resultados obtenidos en las simulaciones, tanto en lazo abierto como en lazo cerrado, se puede establecer que la regla obtenida para los coeficientes de fricción (8) puede usarse para obtener simuladores equivalentes con los tres modelos. Sin embargo, el modelo de ADAMS-Simulink es muy restrictivo, ya que los coeficientes <math display="inline">b</math> y <math display="inline">{f}_{c}</math> dependen de <math display="inline">\mu</math> . Por ejemplo, si se desea desarrollar un simulador considerando los coeficientes de fricción viscosa y de Coulomb de manera independiente, entonces lo más adecuado sería usar la función ODE45 o Simscape Multibody, y no utilizar el cuadro de diálogo de la fricción de la articulación en ADAMS. Si se utiliza la función ODE45, entonces es necesario contar con la ecuación de movimiento para su implementación. Sin embargo, cuando se trata de obtener las ecuaciones de movimiento de mecanismos con más GDL se vuelve un proceso más complejo. La recomendación para solucionar ambos problemas es desarrollar modelos en Simscape Multibody, debido a que no es necesario el desarrollo de ecuaciones de movimiento y sus coeficientes de fricción se pueden configurar de manera independiente.

Para complementar este trabajo de investigación existe otra manera de agregar amortiguamiento viscoso en ADAMS. Para ello es necesario utilizar un resorte de torsión. El resorte de torsión permite configurar la rigidez y el amortiguamiento viscoso. Por lo cual para este trabajo se desactivó la rigidez permitiendo el análisis solo con amortiguamiento viscoso. En la Figura 17 se muestran los resultados de la simulación en lazo abierto utilizando esta metodología. 

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[File:Simulación en lazo abierto con amortiguamiento viscoso.jpg|centre|thumb]]<span style="text-align: center; font-size: 75%;">'''Figura 17.'''  Simulación en lazo abierto considerando un modelo con amortiguamiento viscoso. Los torques de entrada usados fueron: (a) τ = 0 N·m y (b) </span><math display="inline">\tau =5sin\left( \frac{2\pi }{3}\cdot t\right)</math> <span style="text-align: center; font-size: 75%;">N·m</span> </div>El modelo de ADAMS con amortiguamiento viscoso es una buena opción para simular sistemas que no requieren una configuración de fricción compleja, ya que es un software muy utilizado. A diferencia de SimScape  Multibody construir el modelo en ADAMS requiere más pasos. Con ambos sistemas es posible obtener resultados similares siempre y cuando no se utilice la fricción de la articulación. Es recomendable utilizar el amortiguamiento viscoso que viene incluido en el resorte de torsión.

=4. Conclusiones=

En este artículo se presenta un caso de estudio de simulación dinámica usando el modelo matemático de un mecanismo para compararlo con ADAMS y Simscape Multibody. Se identificó el modelo de fricción de ADAMS y después se validó mediante simulaciones en lazo abierto y cerrado. El modelo de Simscape Multibody no requiere del desarrollo de ecuaciones de movimiento y los parámetros de simulación como condiciones iniciales y fricción se ajustan directamente. Por otro lado, el modelo ADAMS-Simulink requiere de varios pasos para reconfigurar parámetros de simulación y el modelo de espacio de estados requiere de un análisis para obtener la ecuación de movimiento. Simscape Multibody es una buena opción para el desarrollo de simuladores de mecanismos en lazo cerrado, sobre todo para estudiantes e investigadores que no cuentan con una plataforma experimental para el desarrollo de sistemas de control aplicado a robots manipuladores o cualquier otro sistema mecánico. Este estudio puede contribuir al desarrollo de laboratorios virtuales y apoyar el estudio a distancia. Las posibilidades de este estudio van desde un mecanismo simple hasta sistemas más complejos. Como trabajo futuro se tiene contemplado cubrir mecanismos con más GDL incluyendo movimientos restringidos usando fuerzas de contacto.

==Agradecimientos==

Los autores agradecen a la Universidad Autónoma de Aguascalientes por los apoyos otorgados en los proyectos PII21-3, PII19-2 y PII22-4.

==Referencias==

[1] Z. H. Khan and M. I. Abid, Distance learning in engineering education: Challenges and opportunities during COVID-19 pandemic crisis in Pakistan, Int. J. Electr. Eng. Educ., 2021.

[2] F. Ghayoor, A MATLAB-based virtual robotics laboratory: Demonstrated by a two-wheeled inverted pendulum, Int. J. Electr. Eng. Educ., vol. 57, no. 4, pp. 301–320, 2020.

[3] K. Odry, R. Fullér, I. J. Rudas, and P. Odry, Fuzzy control of self-balancing robots: A control laboratory project, Comput. Appl. Eng. Educ., vol. 28, no. 3, pp. 512–535, 2020.

[4] V. M. Petrović, B. Nikolić, K. Jovanovic, and V. Potkonjak, Development of virtual laboratory for mechatronic systems, in Advances in Robot Design and Intelligent Control (A. Rodic and T. Borangiu, eds.), (Cham), pp. 622–630, Springer International Publishing, 2017.

[5] K. S. Alli, A LabVIEW-based online DC servomechanism control experiments incorporating PID controller for students’ laboratory, Int. J. Electr. Eng. Educ., 2019.

[6] E. Fabregas, S. Dormido-Canto, and S. Dormido, Virtual and Remote Laboratory with the Ball and Plate System, IFAC-PapersOnLine, vol. 50, no. 1, pp. 9132–9137, 2017. 20th IFAC World Congress.

[7] V. Potkonjak, M. Vukobratovic, K. Jovanović, and M. Medenica, Virtual mechatronic/robotic laboratory - a step further in distance learning, Comput. Educ., vol. 55, no. 2, pp. 465–475, 2010.

[8] O. Eldirdiry, R. Zaier, A. Al-Yahmedi, I. Bahadur, and F. Alnajjar, Modeling of a biped robot for investigating foot drop using MATLAB/ Simulink, Simul. Model Pract. Theory, vol. 98, p. 101972, 2020.

[9] H. Milanowski and A. K. Pilat, Comparison of Identified and SimScape Model of Human Leg Motion, in 2020 International Conference Mechatronic Systems and Materials (MSM), pp. 1–6, 2020.

[10] C. Urrea, L. Valenzuela, and J. Kern, Design, Simulation, and Control of a Hexapod Robot in Simscape Multibody, in Applications from Engineering with MATLAB Concepts (J. Valdman, ed.), Chapters, IntechOpen, September 2016.

[11] Sameh I. Beaber, Abdelrahman S. Zaghloul, Mohamed A. Kamel, Wessam M. Hussein. Dynamic Modeling and Control of the Hexapod Robot Using Matlab SimMechanics. Proceedings of the ASME 2018 International Mechanical Engineering Congress and Exposition, Volume 4A: Dynamics, Vibration, and Control, Pittsburgh, Pennsylvania, USA, November 2018. V04AT06A036.

[12] L. Zhou, W. Chen, J. Wang, J. Liu, W. Chen, and S. Bai, Optimization design of a bionic lower limb rehabilitation robot with dynamic analysis, in 2017 IEEE Int. Conf. Robot. Biomim (ROBIO), pp. 718–723, 2017.

[13] V. Fedak, F. Durovsky, and R. Uveges, Analysis of Robotic System Motion in SimMechanics and MATLAB GUI Environment, in MATLAB (K. Bennett, ed.), ch. 20, Rijeka: IntechOpen, 2014.

[14] A. Othman, K. Belda, and P. Burget, Physical modelling of energy consumption of industrial articulated robots, in 2015 15th Int. Conf. Control Autom. Syst. (ICCAS), pp. 784–789, 2015.

[15] M. Papananias, S. Sztendel, and C. Pislaru, Development of a Novel MultiBody Mechatronic Model for Five-Axis CNC Machine Tool, in Laser Metrology and Machine Performance XI, LAMDAMAP 2015, pp. 379–389, Huddersfield, UK: EUSPEN, March 2015.

[16] V. H. Andaluz, F. A. Chicaiza, C. Gallardo, W. X. Quevedo, J. Varela, J. S. Sánchez, and O. Arteaga, Unity3D-MatLab Simulator in Real Time for Robotics Applications, in Augmented Reality, Virtual Reality, and Computer Graphics (L. T. De Paolis and A. Mongelli, eds.), (Cham), pp. 246–263, Springer International Publishing, 2016.

[17] Mariappan S.M. and A. Veerabathiran, Modelling and simulation of multi spindle drilling redundant SCARA robot using SolidWorks and MATLAB/SimMechanics, Rev Fac Ing-Univ Ant, pp. 63–72, Dec. 2016.

[18] M. F. Velásquez-Lobo, J. M. Ramirez-Cortés, J. de Jesus Rangel-Magdaleno, and J. L. Vázquez-González, Modeling a biped robot on Matlab/SimMechanics, in CONIELECOMP 2013, 23rd Int. Conf. Electron. Commun. Comput., pp. 203–206, 2013.

[19] D. Xia, S. Guo, L. Shi, H. Xing, X. Hou, Y. Liu, H. Liu, Y. Hu, and Z. Li, Application of ADAMS User-Written Subroutine to Simulation of Multi-gait for Spherical Robot, in 2019 IEEE Int. Conf. Mechatron. Autom. (ICMA), pp. 1436–1441, 2019.

[20] C. Liu, D. Zhao, W. Huang, X. Meng, and W. Du, Modeling and simulation of multi-footed wheel robot based on ADAMS, J. Phys. Conf. Ser., vol. 1601, p. 062045, Aug 2020.

[21] D. Sosa-Méndez, E. Lugo-González, M. Arias-Montiel, and R. A. García-García, ADAMS-MATLAB co-simulation for kinematics, dynamics, and control of the Stewart-Gough platform, Int. J. Adv. Robot. Syst., vol. 14, no. 4, p. 1729881417719824, 2017.

[22] Tejaswin Parthasarathy, Vignesh Srinivasaragavan and Soundarapandian Santhanakrishnan, ADAMS-MATLAB Co-Simulation of a Serial Manipulator, MATEC Web Conf., vol. 95, p. 08002, 2017.

[23] A. S. Nair and D. Ezhilarasi, Performance Analysis of Super Twisting Sliding Mode Controller by ADAMS-MATLAB Co-simulation in Lower Extremity Exoskeleton, Int. J. Precis. Eng. Manuf. - Green Technol., vol. 7, no. 3, pp. 743–754, 2020.

[24] C. Liao, J. Zhang, W. Chen, and M. Lv, Union simulation on lower limbs rehabilitation robot based on MATLAB and ADAMS, in 2015 IEEE 10th Conf. Ind. Electron. Appl. (ICIEA), pp. 557–562, 2015.

[25] J. Pandey, N. S. Reddy, R. Ray, and S. N. Shome, Multi-body dynamics of a swimming frog: A co-simulation approach, in 2013 IEEE Int. Conf. Robot. Biomim. (ROBIO), pp. 842–847, 2013.

[26] T. Brezina, Z. Hadas, and J. Vetiska, Using of Co-simulation ADAMSSIMULINK for development of mechatronic systems, in 14th Int. Conf. Mechatron., pp. 59–64, 2011.

[27] Y. Huang, J. Wang, and D. M. Tham, Kinematics modelling and ADAMS MATLAB/Simulink co-simulation for automated aerobridge docking process, in 2017 3rd Int. Conf. Control Autom. Robot. (ICCAR), pp. 304–309, 2017.

[28] P. Chattopadhyay, H. Dikshit, A. Majumder, S. Ghoshal, and A. Maity, Dynamic analysis of a bio-inspired climbing robot using Adams-Simulink co-simulation, AIP Conference Proceedings, vol. 1952, no. 1, p. 020015, 2018.

[29] S. Xiao, H. Deng, X. Zhang, and X. Ding, Optimization design of flapping mechanism of micro air vehicle based on Matlab and Adams, in 2018 IEEE Int. Conf. Robot. Biomim. (ROBIO), pp. 312–317, 2018.

[30] P. I. Corke, Robotics, Vision & Control: Fundamental Algorithms in MATLAB. Springer, second ed., 2017. ISBN 978-3-319-54413-7.

[31] F. Reyes, Robótica: Control de Robots Manipuladores. Alfaomega Grupo Editor, 2011.