En la ciencia física de la dinámica , la dinámica de cuerpos rígidos estudia el movimiento de sistemas de cuerpos interconectados bajo la acción de fuerzas externas . La suposición de que los cuerpos son rígidos (es decir, no se deforman bajo la acción de fuerzas aplicadas) simplifica el análisis, al reducir los parámetros que describen la configuración del sistema a la traslación y rotación de sistemas de referencia adjuntos a cada cuerpo. [1] [2] Esto excluye los cuerpos que muestran un comportamiento fluido , altamente elástico y plástico .
La dinámica de un sistema de cuerpo rígido se describe mediante las leyes de la cinemática y por la aplicación de la segunda ley de Newton ( cinética ) o su forma derivada, la mecánica lagrangiana . La solución de estas ecuaciones de movimiento proporciona una descripción de la posición, el movimiento y la aceleración de los componentes individuales del sistema, y en general del sistema mismo, en función del tiempo . La formulación y solución de la dinámica de cuerpos rígidos es una herramienta importante en la simulación por computadora de sistemas mecánicos .
Si un sistema de partículas se mueve paralelamente a un plano fijo, se dice que el sistema está limitado a un movimiento plano. En este caso, las leyes de Newton (cinética) para un sistema rígido de N partículas, P i , i =1,..., N , se simplifican porque no hay movimiento en la dirección k . Determine la fuerza y el par resultantes en un punto de referencia R , para obtener
donde r i denota la trayectoria plana de cada partícula.
La cinemática de un cuerpo rígido produce la fórmula para la aceleración de la partícula P i en términos de la posición R y la aceleración A de la partícula de referencia, así como el vector de velocidad angular ω y el vector de aceleración angular α del sistema rígido de partículas como ,
Para sistemas que están restringidos a un movimiento plano, los vectores de velocidad angular y aceleración angular se dirigen a lo largo de k perpendicular al plano de movimiento, lo que simplifica esta ecuación de aceleración. En este caso, los vectores de aceleración se pueden simplificar introduciendo los vectores unitarios e i desde el punto de referencia R a un punto ri y los vectores unitarios , entonces
Esto produce la fuerza resultante sobre el sistema como
donde y es el vector unitario perpendicular al plano para todas las partículas P i .
Utilice el centro de masa C como punto de referencia, de modo que estas ecuaciones para las leyes de Newton se simplifiquen y se conviertan en
donde M es la masa total y I C es el momento de inercia alrededor de un eje perpendicular al movimiento del sistema rígido y que pasa por el centro de masa.
Se han desarrollado varios métodos para describir orientaciones de un cuerpo rígido en tres dimensiones. Se resumen en las siguientes secciones.
El primer intento de representar una orientación se atribuye a Leonhard Euler . Imaginó tres sistemas de referencia que podían rotar uno alrededor del otro, y se dio cuenta de que comenzando con un sistema de referencia fijo y realizando tres rotaciones, podía obtener cualquier otro sistema de referencia en el espacio (usando dos rotaciones para fijar el eje vertical y otra para fijar el eje vertical). fijar los otros dos ejes). Los valores de estas tres rotaciones se denominan ángulos de Euler . Comúnmente se utiliza para denotar precesión, nutación y rotación intrínseca.
Estos son tres ángulos, también conocidos como guiñada, cabeceo y balanceo, ángulos de navegación y ángulos cardán. Matemáticamente constituyen un conjunto de seis posibilidades dentro de los doce conjuntos posibles de ángulos de Euler, siendo el orden el que mejor se utiliza para describir la orientación de un vehículo como un avión. En ingeniería aeroespacial se les suele denominar ángulos de Euler.
Euler también se dio cuenta de que la composición de dos rotaciones equivale a una sola rotación alrededor de un eje fijo diferente ( teorema de rotación de Euler ). Por tanto, la composición de los tres primeros ángulos tiene que ser igual a una sola rotación, cuyo eje era complicado de calcular hasta que se desarrollaron las matrices.
Basándose en este hecho, introdujo una forma vectorial de describir cualquier rotación, con un vector en el eje de rotación y un módulo igual al valor del ángulo. Por tanto, cualquier orientación puede representarse mediante un vector de rotación (también llamado vector de Euler) que conduce a ella desde el sistema de referencia. Cuando se utiliza para representar una orientación, el vector de rotación se denomina comúnmente vector de orientación o vector de actitud.
Un método similar, llamado representación eje-ángulo , describe una rotación u orientación utilizando un vector unitario alineado con el eje de rotación y un valor separado para indicar el ángulo (ver figura).
Con la introducción de las matrices se reescribieron los teoremas de Euler. Las rotaciones se describieron mediante matrices ortogonales denominadas matrices de rotación o matrices de cosenos direccionales. Cuando se utiliza para representar una orientación, una matriz de rotación se denomina comúnmente matriz de orientación o matriz de actitud.
El vector de Euler mencionado anteriormente es el vector propio de una matriz de rotación (una matriz de rotación tiene un valor propio real único ). El producto de dos matrices de rotación es la composición de rotaciones. Por lo tanto, como antes, la orientación se puede dar como la rotación desde el cuadro inicial para lograr el cuadro que queremos describir.
El espacio de configuración de un objeto no simétrico en un espacio n -dimensional es SO( n ) × R n . La orientación se puede visualizar adjuntando una base de vectores tangentes a un objeto. La dirección en la que apunta cada vector determina su orientación.
Otra forma de describir las rotaciones es utilizando cuaterniones de rotación , también llamados versores. Son equivalentes a matrices de rotación y vectores de rotación. Con respecto a los vectores de rotación, se pueden convertir más fácilmente desde y hacia matrices. Cuando se utilizan para representar orientaciones, los cuaterniones de rotación suelen denominarse cuaterniones de orientación o cuaterniones de actitud.
Para considerar la dinámica de un cuerpo rígido en un espacio tridimensional, se debe ampliar la segunda ley de Newton para definir la relación entre el movimiento de un cuerpo rígido y el sistema de fuerzas y momentos de torsión que actúan sobre él.
Newton formuló su segunda ley para una partícula como: "El cambio de movimiento de un objeto es proporcional a la fuerza aplicada y se realiza en la dirección de la línea recta en la que se aplica la fuerza". [3] Debido a que Newton generalmente se refería a la masa multiplicada por la velocidad como el "movimiento" de una partícula, la frase "cambio de movimiento" se refiere a la masa multiplicada por la aceleración de la partícula, por lo que esta ley generalmente se escribe como
Si un sistema de N partículas, P i , i=1,..., N , se ensamblan en un cuerpo rígido, entonces la segunda ley de Newton se puede aplicar a cada una de las partículas del cuerpo. Si F i es la fuerza externa aplicada a la partícula P i con masa m i , entonces
Una simplificación importante de estas ecuaciones de fuerza se obtiene introduciendo la fuerza y el par resultantes que actúan sobre el sistema rígido. Esta fuerza y par resultantes se obtienen eligiendo una de las partículas del sistema como punto de referencia, R , donde cada una de las fuerzas externas se aplica con la adición de un par asociado. La fuerza resultante F y el par T vienen dados por las fórmulas,
La segunda ley de Newton para una partícula se combina con estas fórmulas para que la fuerza y el par resultantes cedan,
Las propiedades de masa del cuerpo rígido están representadas por su centro de masa y matriz de inercia . Elija el punto de referencia R para que satisfaga la condición
entonces se le conoce como centro de masa del sistema.
La matriz de inercia [I R ] del sistema con respecto al punto de referencia R está definida por
donde está el vector columna R i − R ; es su transpuesta y es la matriz identidad de 3 por 3.
es el producto escalar de consigo mismo, mientras que es el producto tensorial de consigo mismo.
Usando el centro de masa y la matriz de inercia, las ecuaciones de fuerza y torsión para un solo cuerpo rígido toman la forma
La dinámica de un sistema interconectado de cuerpos rígidos, Bi , j = 1, ..., M , se formula aislando cada cuerpo rígido e introduciendo las fuerzas de interacción. La resultante de las fuerzas externas y de interacción sobre cada cuerpo produce las ecuaciones de fuerza-torque.
La formulación de Newton produce 6 M ecuaciones que definen la dinámica de un sistema de M cuerpos rígidos. [4]
Un objeto en rotación, ya sea bajo la influencia de pares de torsión o no, puede exhibir comportamientos de precesión y nutación . La ecuación fundamental que describe el comportamiento de un cuerpo sólido en rotación es la ecuación de movimiento de Euler :
La solución a esta ecuación cuando no se aplica ningún par se analiza en los artículos Ecuación de movimiento de Euler y Elipsoide de Poinsot .
De la ecuación de Euler se deduce que un par τ aplicado perpendicular al eje de rotación, y por lo tanto perpendicular a L , da como resultado una rotación alrededor de un eje perpendicular tanto a τ como a L. Este movimiento se llama precesión . La velocidad angular de precesión Ω P viene dada por el producto cruzado : [ cita necesaria ]
La precesión se puede demostrar colocando una peonza con su eje horizontal y apoyada sin apretar (sin fricción hacia la precesión) en un extremo. En lugar de caer, como podría esperarse, la peonza parece desafiar la gravedad al permanecer con su eje horizontal, cuando el otro extremo del eje queda sin apoyo y el extremo libre del eje describe lentamente un círculo en un plano horizontal, el resultado giro de precesión. Este efecto se explica por las ecuaciones anteriores. El par en la parte superior es proporcionado por un par de fuerzas: la gravedad que actúa hacia abajo sobre el centro de masa del dispositivo y una fuerza igual que actúa hacia arriba para sostener un extremo del dispositivo. La rotación resultante de este par no es hacia abajo, como podría esperarse intuitivamente, provocando la caída del dispositivo, sino perpendicular tanto al par gravitacional (horizontal y perpendicular al eje de rotación) como al eje de rotación (horizontal y hacia afuera del eje de rotación). punto de apoyo), es decir, alrededor de un eje vertical, haciendo que el dispositivo gire lentamente alrededor del punto de apoyo.
Bajo un par constante de magnitud τ , la velocidad de precesión Ω P es inversamente proporcional a L , la magnitud de su momento angular:
Por convención, estos tres vectores (par, giro y precesión) están orientados entre sí según la regla de la mano derecha .
Se obtiene una formulación alternativa de la dinámica de cuerpos rígidos que tiene varias características convenientes considerando el trabajo virtual de las fuerzas que actúan sobre un cuerpo rígido.
El trabajo virtual de las fuerzas que actúan en varios puntos sobre un solo cuerpo rígido se puede calcular utilizando las velocidades de su punto de aplicación y la fuerza y el par resultantes . Para ver esto, dejemos que las fuerzas F 1 , F 2 ... F n actúen sobre los puntos R 1 , R 2 ... R n en un cuerpo rígido.
Las trayectorias de R i , i = 1, ..., n están definidas por el movimiento del cuerpo rígido. Las velocidades de los puntos R i a lo largo de sus trayectorias son
El trabajo se calcula a partir del producto escalar de cada fuerza con el desplazamiento de su punto de contacto.
o recogiendo los coeficientes de δq j
Para simplificar, considere una trayectoria de un cuerpo rígido que está especificada por una única coordenada generalizada q, como un ángulo de rotación, entonces la fórmula se convierte en
Introduzca la fuerza resultante F y el par T para que esta ecuación tome la forma
La cantidad Q definida por
se conoce como fuerza generalizada asociada con el desplazamiento virtual δq. Esta fórmula se generaliza al movimiento de un cuerpo rígido definido por más de una coordenada generalizada, es decir
Es útil observar que las fuerzas conservativas como la gravedad y las fuerzas del resorte se pueden derivar de una función potencial V ( q 1 ,..., q n ) , conocida como energía potencial . En este caso las fuerzas generalizadas están dadas por
Las ecuaciones de movimiento para un sistema mecánico de cuerpos rígidos se pueden determinar utilizando la forma del principio de trabajo virtual de D'Alembert. El principio de trabajo virtual se utiliza para estudiar el equilibrio estático de un sistema de cuerpos rígidos; sin embargo, al introducir términos de aceleración en las leyes de Newton, este enfoque se generaliza para definir el equilibrio dinámico.
El equilibrio estático de un sistema mecánico de cuerpos rígidos se define por la condición de que el trabajo virtual de las fuerzas aplicadas sea cero para cualquier desplazamiento virtual del sistema. Esto se conoce como el principio del trabajo virtual. [5] Esto equivale al requisito de que las fuerzas generalizadas para cualquier desplazamiento virtual sean cero, es decir Q i =0.
Sea un sistema mecánico construido a partir de n cuerpos rígidos, Bi , i = 1, ..., n , y sea la resultante de las fuerzas aplicadas sobre cada cuerpo los pares fuerza-par , Fi y Ti , i = 1, ..., norte . Observe que estas fuerzas aplicadas no incluyen las fuerzas de reacción donde los cuerpos están conectados. Finalmente, supongamos que la velocidad V i y las velocidades angulares ω i , i = 1, ..., n , para cada cuerpo rígido, están definidas por una única coordenada generalizada q. Se dice que tal sistema de cuerpos rígidos tiene un grado de libertad .
El trabajo virtual de las fuerzas y pares, Fi y Ti , aplicados a este sistema de un grado de libertad está dado por
Si el sistema mecánico está definido por m coordenadas generalizadas, q j , j = 1, ..., m , entonces el sistema tiene m grados de libertad y el trabajo virtual viene dado por,
Estas m ecuaciones definen el equilibrio estático del sistema de cuerpos rígidos.
Considere un solo cuerpo rígido que se mueve bajo la acción de una fuerza resultante F y un par T , con un grado de libertad definido por la coordenada generalizada q . Supongamos que el punto de referencia para la fuerza y el par resultantes es el centro de masa del cuerpo, entonces la fuerza de inercia generalizada Q* asociada con la coordenada generalizada q viene dada por
Esta fuerza de inercia se puede calcular a partir de la energía cinética del cuerpo rígido,
Un sistema de n cuerpos rígidos con m coordenadas generalizadas tiene la energía cinética
La forma de D'Alembert del principio de trabajo virtual establece que un sistema de cuerpos rígidos está en equilibrio dinámico cuando el trabajo virtual de la suma de las fuerzas aplicadas y las fuerzas de inercia es cero para cualquier desplazamiento virtual del sistema. Por tanto, el equilibrio dinámico de un sistema de n cuerpos rígidos con m coordenadas generalizadas requiere que
Si las fuerzas generalizadas Q j son derivables de una energía potencial V ( q 1 , ..., q m ) , entonces estas ecuaciones de movimiento toman la forma
En este caso, introduzca el lagrangiano , L = T − V , de modo que estas ecuaciones de movimiento se conviertan en
El momento lineal y angular de un sistema rígido de partículas se formula midiendo la posición y la velocidad de las partículas en relación con el centro de masa. Sea el sistema de partículas P i , i = 1, ..., n ubicado en las coordenadas r i y velocidades v i . Seleccione un punto de referencia R y calcule la posición relativa y los vectores de velocidad,
Los vectores de momento lineal y angular totales relativos al punto de referencia R son
Si se elige R como centro de masa, estas ecuaciones se simplifican a
Para especializar estas fórmulas en un cuerpo rígido, supongamos que las partículas están conectadas rígidamente entre sí, de modo que P i , i=1,...,n estén ubicadas en las coordenadas r i y las velocidades vi . Seleccione un punto de referencia R y calcule la posición relativa y los vectores de velocidad,
El momento lineal y el momento angular de este sistema rígido medidos con respecto al centro de masa R es
Estas ecuaciones se simplifican para convertirse en,