Resumos
A parametrização da velocidade angular instantânea de um corpo rígido é usualmente descrita nos cursos de mecânica clássica com o auxílio dos ângulos de Euler. No entanto, para evitar problemas de singularidade na integração numérica das equações de movimento de rotação, outras formas de parametrização são mais interessantes. Dentre as diversas opções, os quatérnions são uma escolha bastante popular e com forte apelo prático na parametrização e integração das equações de Euler. Neste sentido, o presente artigo contribui com uma apresentação didática e comparação da integração numérica das equações de Euler para descrição do movimento de um giroscópio com a cinemática inversa solucionada via ângulos de Euler e por quatérnions, que não são comumente descritas em livros didáticos clássicos em um primeiro curso de mecânica. Busca-se mostrar que além deste trunfo numérico, a integração das equações de movimento com parametrização via quatérnions também pode ser usada com um fim didático mais nobre: introduzir e motivar os estudantes de um curso de mecânica, tanto físicos como engenheiros, no uso de outros sistemas algébricos mais avançados, como a Álgebra Geométrica (Álgebra de Clifford).
Palavras-chave:
Rotações; Giroscópio; Orientação no espaço; Singularidades; Ângulos de Euler; Quatérnions; Álgebra geométrica
The angular velocity of a rigid body is usually described in classical mechanics courses by using Euler angles. However, in order to avoid singularity problems on the numerical integration of the motion equations, another ways are more interesting. Among them, the quaternions are a very popular choice and with a strong practical appeal on the parameterization and integration of the Euler's equations. In this sense, this paper proposes a comparison between the numerical integration of the Euler equations of a gyroscope with the inverse kinematics solved through Euler angles and by quaternions, that are not commonly described in classical mechanics textbooks on a first mechanical course. The goal is to show that apart the numerical asset, the integration of the motion equations with the quaternions parameterization also can be used with a noblest didactical way: to introduce and motivate the students with the use of others more advanced algebraic systems, as Geometric Algebra (Clifford's Algebra).
Keywords:
Gyroscope; Quaternions; Euler Angles; Geometric Algebra; Rigid Body Rotation in Space
1. Introdução
A parametrização de rotações finitas de um corpo rígido no Espaço Euclidiano foi descrita magistralmente por Euler em 1775 a partir do uso de três rotações sequenciais. Este resultado foi um corolário obtido a partir de seu teorema, conhecido hoje como Teorema de Rotação de Euler, que demonstrou que qualquer rotação de um corpo com um ponto fixo pode ser descrita por uma rotação própria em torno de um eixo de rotação, conhecido como Eixo de Euler, solidário ao giro do próprio corpo [1[1] L. Euler, Novi Acad. Sci. Petrop 20, 189 (1775)., 2[2] J. Stuelpnagel, SIAM review 6, 422 (1964).].
Até os dias atuais esta é uma das formas mais populares para representar a orientação de um corpo rígido no espaço. Os três ângulos de Euler parametrizam três rotações sequenciais [3[3] J. Diebel, Matrix 58, 15 (2006).] e são bastante visuais e claros, sendo utilizados para girar um corpo no espaço a partir de uma matriz de rotação com nove elementos, conhecidos como cossenos diretores. Porém, uma vez que o produto de matrizes não é, por via de regra, comutativo, a escolha e ordem destas rotações sequenciais deve ser definida a priori em função do tipo de problema. Em especial, um cuidado adicional deve ser levado em conta no momento desta escolha, pois as velocidades angulares instantâneas também são descritas por estes ângulos e suas respectivas taxas de variação. Um dos graves inconvenientes do uso dos ângulos de Euler para descrição destas velocidades é quando realiza-se algum cálculo envolvendo cinemática inversa, ou seja, quando a partir do conhecimento das velocidades angulares instantâneas, deseja-se integrar numericamente as equações para obter os ângulos que fornecem a orientação do corpo no espaço. Em algumas situações, existe uma clara singularidade que inviabiliza este cálculo de maneira direta [3[3] J. Diebel, Matrix 58, 15 (2006).]. Esse problema será discutido em mais detalhes na seção 4 deste artigo.
Felizmente, existem outras maneiras de parametrização que podem ser usadas com algumas vantagens frente ao uso dos ângulos de Euler [3[3] J. Diebel, Matrix 58, 15 (2006).]. Dentre estas formas, podem ser citados os parâmetros de Cayley-Klein [4[4] J. Neto, Mecânica Newtoniana, Lgrangiana e Hamiltoniana (Editora Livraria da Física, São Paulo, 2004).], parâmetros de Rodrigues e Rodrigues Modificado [5[5] P. Singla, D. Mortari and J.L. Junkins, in: Proceedings of the AAS Space Flight Mechanics Conference, Hawaii, 2004., 6[6] P. Betsch, A. Menzel and E. Stein, Computer Methods in Applied Mechanics and Engineering 155, 273 (1998).], vetor conforme [7[7] M. Trindade and R. Sampaio, Journal of the Brazilian Society of Mechanical Sciences 22, 341 (2000).] e, ainda, os quatérnions [8[8] J. Kuipers, Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (Princeton Paperbacks, Princeton University Press, 1999).]. Uma comparação das vantagens e desvantagens de algumas formas de parametrização pode ser encontrada em Trindade e Sampaio [7[7] M. Trindade and R. Sampaio, Journal of the Brazilian Society of Mechanical Sciences 22, 341 (2000).]. Entretanto estas parametrizações de rotação não são discutidas em detalhes na maioria dos cursos de mecânica e são apresentadas em livros didáticos de mecânica clássica apenas de forma superficial [9[9] K. Vogtmann, V. Arnold and A. Weinstein, Mathematical Methods of Classical Mechanics (Graduate Texts in Mathematics) (Springer, New York, 1989).–12[12] E.T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid Bodies: With an Introduction to the Problem of Three Bodies (CPU Archive, London, 1970).].
Dentre todas estas formas de parametrização, os quatérnions são um tanto quanto especiais, principalmente devido a sua origem histórica. William Rowan Hamilton, em meados do século XIX, sabia que multiplicação de números complexos eram isomorfos de rotação no e buscou de forma persistente estender isto para rotações no [13[13] W.R. Hamilton, Elements of Quaternions (Longmans, Green, & Company, Londres, 1866).]. Suas primeiras tentativas envolveram utilizar elementos no para operar rotações no , mas não obteve sucesso em obter uma álgebra fechada. Somente em 1843, Hamilton constatou que precisaria operar elementos no para rotacionar elementos no . Este novo sistema numérico foi chamado por ele de quatérnion, sendo representado hoje por e com uma propriedade bem interessante: quatérnions podem ser multiplicados, ao contrário de vetores1 1 Um objeto matemático corriqueiro hoje para descrever grandezas físicas, mas que em meados do século XIX não tinha a definição que temos hoje. . Os quatérnions foram a primeira álgebra não-comutativa da história e Hamilton ficou bastante empolgado com o leque de possíveis aplicações em física e engenharia. Tanto que passou o resto de sua vida escrevendo e divulgando sobre os quatérnions, a que ele se referia como uma linguagem matemática universal para descrever fenômenos físicos diversos [14[14] W.R. Hamilton, The Mathematical Papers of Sir William Rowan Hamilton, (CUP Archive, London, 2000)., 15[15] W.R. Hamilton, Lectures on Quaternions, (The Royal Irish Academy, Dublin, 1853).]. Infelizmente, no final do século XIX, o americano Gibbs e o inglês Heaviside popularizaram uma notação e “álgebra” de uma vertente, podendo ser entendida hoje como uma “versão simplificada” de quatérnions, separando o produto de dois quatérnions em duas operações separadas conhecidas hoje como produto vetorial, com o símbolo ; e produto escalar, com o símbolo , na notação de Gibbs [16[16] M. Menon, Revista Brasileira de Ensino de Física 31, 2305 (2009)., 17[17] Gustavo Martínez-Sierra y Pierre Francois Benoit Poirier, Latin-American Journal of Physics Education 2, 201 (2008).]. Apesar de debates ferrenhos na época, a álgebra de quatérnions foi deixada no segundo plano e a “álgebra vetorial” de Gibbs, mesmo que com certos problemas e limitações, foi popularizada rapidamente, em especial devido a inúmeros livros didáticos escritos no início do século XX.
Apesar disto, na segunda metade do século XX, os quatérnions começaram a ser “redescobertos” e aplicados de forma maciça, especialmente em problemas de cinemática inversa em robótica e engenharia aeroespacial [18[18] M.J. Amoruso, Euler Angles and Quaternions in Six Degree of Freedom Simulations of Projectiles - Technical Report (DTIC Document, Fort Belvoir, 1996).–20[20] M.D. Shuster, Navigation 8, 439 (1993).], visão computacional [21[21] K. Shoemake, Quaternion Calculus and Fast Animation, Computer Animation: 3-D Motion Specification and Control (Siggraph, Anaheim, 1987)., 22[22] J.C. Hart, G.K. Francis and L.H. Kauffman, ACM Transactions on Graphics 13, 256 (1994).], biologia molecular [23[23] W.R. Rudnicki, G. Bakalarski and B. Lesyng, Journal of Biomolecular Structure and Dynamics 17, 1097 (2000).–25[25] T. Róg, K. Murzyn, K. Hinsen and G.R. Kneller, Journal of Computational Chemistry 24, 657 (2003).], mecânica dos fluídos [26[26] W. Sprössig, Advances in Applied Clifford Algebras 17, 259 (2007).–28[28] H. Eshraghi and J. Gibbon, Journal of Physics A: Mathematical and Theoretical 41, 344004 (2008).], mecânica quântica [29[29] J. Vaz Jr., Revista Brasileira de Ensino de Física 19, 234 (1997).], entre inúmeras outras áreas. Porém, em grande parte destes artigos e trabalhos, em especial os de cunho técnico e de aplicações, os quatérnions são descritos como uma mera ferramenta matemática, ou método, usados para facilitar a integração numérica e evitar singularidades no problema de cinemática inversa. Carrera [30[30] D.H.Z. Carrera, Movimento de Rotação Sem Restrição de um Corpo Rígido. Tese de Doutorado, PUC, Rio de Janeiro, 2010.] afirma que os valores correspondentes das componentes do quatérnion não correspondem a nenhum significado físico, enquanto Oliveira [31[31] A.A.G. Oliveira, Estudo Estatístico dos Processos Envolvidos em uma Plataforma de Atitude Solidária. Dissertação de Mestrado, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, 2000.] apresenta o uso de quatérnions como uma vantagem de cunho apenas prático na determinação da orientação de um veículo, por meio de um sistema de referência solidário ao corpo, cujo único objetivo é evitar singularidades. Alega-se, ainda, que pouco importa a facilidade em extrair algum significado físico dos quatérnions, já que eles podem ser facilmente convertidos para ângulos de Euler e que, estes sim, possuem a vantagem de permitirem o imediato entendimento físico. Altmann [32[32] S.L. Altmann, Mathematics Magazine 62, 291 (1989).], por sua vez, apresenta os quatérnions como geradores de rotações no espaço. Lima e Mello [33[33] R.S. Lima and M.P. Mello, CNMAC - Congresso Nacional de Matemática Aplicada e Computacional (2004).] tratam os quatérnions como um método de interpretação não tão natural, destacando também a sua dificuldade de interpretação. Outros autores, por outro lado, alegam que a representação de rotações por quatérnions, quando comparada à parametrização por ângulos de Euler, tem a vantagem de ser livre da complexidade analítica que envolve os ângulos de Euler [34[34] V. Pisacane, Fundamentals of Space Systems (Oxford University Press, Oxford, 2005).]. Porém, da mesma forma que números complexos, entidades às quais os estudantes já estão, em muitos casos, habituados e que são intuitivas quando tratam-se rotações no plano, a representação de quatérnions puros no espaço tridimensional também o é, já que por si só fornece uma visualização das trajetórias que são muito mais expressivas do que a decomposição em termos de três ângulos de Euler relacionados a um sistema inercial, fornecendo uma representação matemática e espacial muito simples de rotações. Essa facilidade de visualização, quando combinada com a simplicidade de implementação computacional e ao número menor de operações matemáticas necessárias para realizar rotações, não só traz uma redução no tempo de processamento ao representar trajetórias, como é facilmente compreendida uma vez entendido o formalismo matemático por trás dos quatérnions. Outra questão é que, quando optar-se por expressar os resultados dessas rotações em termos de grandezas angulares, pode-se calculá-las facilmente a partir dos quatérnions.
No entanto, o que alguns destes autores que usam quatérnions não ressaltam é que o uso de quatérnions para parametrizar rotações tem uma função e significado muito mais profundo do que apenas servir de um artifício ou método numérico para evitar singularidades. Neste sentido, este artigo tem como meta descrever esta ideia a partir de um problema concreto envolvendo o estudo clássico de orientação de um giroscópio. A integração numérica das equações de movimento que descrevem o movimento de orientação de um giroscópio pode ter a velocidade angular parametrizada com ângulos de Euler e quatérnions para condição livre e com aplicação de torques externos, com ou sem problema de singularidade numérica dependendo das condições iniciais. A partir este problema, o presente artigo discute como é possível estender e apresentar os quatérnions como uma classe algébrica particular de um sistema mais avançado: a álgebra geométrica (álgebra de Clifford). Assim, em um curso de mecânica clássica é possível introduzir a álgebra geométrica como uma classe especial apresentando os quatérnions como uma álgebra fechada com operação de multiplicação não-comutativa entre os seus elementos que, além de um enorme apelo prático na solução de problemas relevantes e de interesse diverso, serve como porta de entrada para a Álgebra Geométrica.
O artigo está organizado em três partes principais. Inicialmente apresenta-se uma introdução à álgebra dos quatérnions, mostrando as suas principais propriedades e como realizar rotações com estes elementos. Na sequência, se revisa as equações de Euler para orientação de um giroscópio no espaço, seguido do problema da cinemática inversa, solucionado via ângulos de Euler e quatérnions. Em seguida, discute-se como entender os quatérnions como uma álgebra de Clifford. Por fim, algumas considerações finais são pontuadas ao final do texto.
2. A Álgebra dos Quatérnions
Os quatérnions consistem em uma generalização dos números complexos para quatro dimensões e são compostos por uma parte escalar e uma parte vetorial . Pode-se tratar os números reais como números hipercomplexos de dimensão e os números complexos comuns como hipercomplexos de dimensão . Nesses casos, relações de comutatividade tanto para soma quanto para multiplicação são sempre válidas. Por outro lado, para qualquer hipercomplexo de dimensão superior a , essas relações não são necessariamente válidas. Um caso especial são os hipercomplexos de dimensão , conhecidos como quatérnions [8[8] J. Kuipers, Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (Princeton Paperbacks, Princeton University Press, 1999).].
A Fig. 1 descreve uma regra mnemônica para representar a não-comutatividade das multiplicações de quatérnions, onde a multiplicação de dois termos em sentido horário fornece como resultado o terceiro termo com sinal positivo. Por exemplo, se for multiplicados tem-se como resultado , já a multiplicação de dois termos subsequentes em sentido anti-horário também fornece o terceiro termo, porém, com sinal negativo, ou seja quando for multiplicado , tem-se como resultado .
Deste modo, observa-se que os componentes unitários ,, do quatérnion respeitam as relações fundamentais:
representado por:
ou por:
ou ainda, por:
sendo
Assim, os quatérnions são compostos por quatro componentes ortogonais. Quando a parte escalar do quatérnion é nula, tem-se um quatérnion puro, que pode, por exemplo, representar um vetor de três dimensões, do mesmo modo escalares podem ser tratados como quatérnions com parte vetorial nula [35[35] R.A. Wehage, in: Computer Aided Analysis and Optimization of Mechanical System Dynamics (Springer, Berlin, 1984), p. 147-180.].
Um quatérnion possui um conjugado e módulo tal que:
Os quatérnions têm fechamento na soma, assim para , , vale que , onde denota o conjunto dos quatérnions. Estes também apresentam fechamento na multiplicação por escalar, assim e , vale . Também são válidos os seguintes axiomas, que mostram que os quatérnions definem um espaço vetorial:
Assim, o conjunto dos quatérnions é um isomorfo do para quatérnions puros e é isomorfo do para os demais casos.
Observa-se também que, por não satisfazer a comutatividade na multiplicação, o conjunto numérico não forma um corpo. Já a inversa de um quatérnion é definida por:
Kuipers [8[8] J. Kuipers, Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (Princeton Paperbacks, Princeton University Press, 1999).] mostra que se for um quatérnion unitário, ou normalizado, a inversa do quatérnion é próprio complexo conjugado, de modo que:
Interessante ressaltar que os quatérnions são excelentes operadores de rotação. Sejam dois vetores , que podem ser descritos como quatérnions puros, de tal modo que e , onde é após a operação de rotação. Pode-se definir um operador de rotações , associado com o quatérnion e aplicado ao vetor , a partir de:
sendo um operador linear representando geometricamente uma rotação de um ângulo em torno de um vetor diretor no , onde fornece a orientação do eixo de rotação do Teorema de Rotação de Euler e o ângulo de rotação em torno deste eixo. A representação em quatérnions é, então, fornecida por um quatérnion unitário na forma:
A operação de rotações com quatérnions pode ser ainda escrita na forma explícita , de modo que:
Essa representação permite a comparação da matriz de rotações composta por quatérnions com a matriz de rotações envolvendo funções trigonométricas dos ângulos de Euler. Essa relação entre as duas matrizes será realizada na seção 4, nas equações (38), (39) e (40).
3. Orientação de um Giroscópio
Um giroscópio pode ser descrito por um disco com massa e raio , posicionado em e girando em torno de seu eixo de simetria, conforme visto na Fig. 2. O sistema inercial é fixo em , que corresponde ao centro de massa, com os eixos representado pela base ortonormal . O primeiro giro, de precessão, ocorre em torno de com velocidade angular e é solidário ao sistema móvel descrito por representado pela base ortonormal . Já o segundo giro, de nutação, ocorre em torno de com velocidade angular e é solidário ao sistema móvel descrito por com a base . Por fim o terceiro giro, de spin, ocorre em torno de com velocidade angular e é, por sua vez, solidário ao sistema móvel descrito por com a base . Todos estes eixos e giros são descritos na Fig. 2. O movimento da base é solidário ao do corpo, e o tensor de inércia é:
O tensor de inércia, calculado a partir dos valores de parâmetros de massa kg e raio m do giroscópio, é representado na base móvel como:
A matriz de rotação entre os sistemas e pode ser escrita como:
Já a velocidade angular do corpo, em função dos ângulos de Euler e representada na base , é descrita por:
sendo , e os componentes da velocidade angular do giroscópio nos eixos principais de inércia.
O momento angular do giroscópio representado no sistema móvel , pode ser descrito pela relação:
sendo o momento angular do giroscópio.
As equações para descrever a orientação do giroscópio são obtidas a partir do cálculo da taxa de variação do momento angular, descrita por:
sendo a velocidade angular do sistema de referência e a soma dos torques externos, representados no sistema de referência móvel .
Calculando a derivada da amplitude do momento angular e a taxa de variação da direção do momento angular, obtém-se as equações de Euler de primeira ordem:
que podem ser integradas numericamente para obtenção das velocidades instantâneas , e . Os ângulos de Euler , e são obtidos pela integração inversa da equação (15).
4. Cinemática Inversa via Ângulos de Euler versus Quatérnions
O problema da cinemática inversa por meio dos ângulos de Euler consiste em relacionar as velocidades de precessão , nutação e spin com as velocidades angulares nas direções principais de inércia, explicitando , e a partir da inversa da eq. (15):
Estas equações podem ser integradas numericamente. No entanto, para , com tem-se que , de modo que a Eq. (19) pode apresentar singularidade. Este inconveniente de singularidade não ocorre utilizando os quatérnions para representar rotações. Para a solução de um problema deste tipo, inicialmente é necessário estimar os quatérnions correspondentes às condições iniciais do problema, o que pode ser feito com a eq. (11). No caso de rotações sequenciais, sejam , e as rotações em torno dos eixos , e , respectivamente:
A rotação total, resultante das três rotações sequenciais, será:
onde, aplicando as propriedades de multiplicação de quatérnions, pode-se obter as componentes do quatérnion :
Calculados os termos , , e , pode-se integrar também as equações de Euler, eq. (18), e relacionar a taxa de variação de um quatérnion com as velocidades angulares na direção principal, escritas na forma de quatérnions. Essa relação pode ser obtida a partir da ideia de que entre dois quatérnions unitários pode haver um terceiro que represente a transição entre eles [8[8] J. Kuipers, Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (Princeton Paperbacks, Princeton University Press, 1999).]. As formas de obtenção deste quatérnion são discutidas em [36[36] R.A. Mayo, Journal of Guidance, Control, and Dynamics 2, 44 (1979).]. Este quatérnion de transição pode ser, então, escrito como:
sendo:
o incremento da transição entre os quatérnions. Seu ângulo de rotação é e ocorre em torno de um eixo definido pelo vetor diretor . Supondo que o incremento é pequeno, assume-se a velocidade angular como constante no intervalo de tempo , e pode-se assumir e , de forma que a eq. (26) pode ser reescrita como:
substituindo este resultado na eq (25), obtém-se:
se tratarmos o vetor como um quatérnion de parte escalar nula,
e dividindo ambos os lados da eq. (28) por e aplicando o limite para :
e, sendo a taxa de variação de :
sendo a velocidade angular do quatérnion incremental :
finalmente, obtém-se:
Onde a velocidade angular está descrita como um quatérnion de parte escalar nula . A eq. (33) pode, então, ser reescrita como:
ou, ainda, a partir da definição de produto entre dois quatérnions:
resolvendo as operações de produto escalar e produto vetorial descritas na eq. (35), obtém-se:
que pode ser reescrita, na forma matricial, como:
Em suma, para resolver o problema de cinemática inversa por quatérnions, deve-se integrar as equações de Euler, eq. (18), e substituir os valores encontrados para as velocidades angulares nas equações (37) e, em seguida, integrar as equações (37) para a obtenção dos quatérnions em cada instante de tempo. Caso opte-se por determinar os ângulos de Euler, pode-se relacionar a matriz de rotação , eq. (14), à forma explícita da operação de rotações com quatérnions , visto na eq. (12). Igualando-se as matrizes termo a termo, ou seja e realizando-se algumas operações algébricas entre esses termos à fim de isolar os ângulos , e , obtém-se:
Neste ponto são realizadas as integrações numéricas das equações (18) e (19) para o problema da cinemática inversa por ângulos de Euler; e das equações (18) e (37) para o problema da cinemática inversa por quatérnions. Todas as integrações foram feitas utilizando o método de Runge-Kutta de 4° ordem, com amostras temporais e passo de integração s. Foi assumido apenas o caso onde o sistema está livre da ação de momentos externos, com N.m.
A Fig. 3 apresenta os ângulos de Euler obtidos para um caso onde não existia singularidade, de modo que a solução por meio da cinemática inversa por ângulos de Euler e por quatérnions apresentaram resultados idênticos entre si. A solução apresentada é aquela obtida por cinemática inversa via ângulos de Euler e optou-se por omitir a solução por quatérnions para evitar a repetição de resultados. Para esta situação, a velocidade de spin foi definida com valor muito superior às velocidades de precessão e nutação. Nota-se, portanto, um ligeiro aumento na velocidade de spin bem como do próprio ângulo de spin, enquanto os ângulos de precessão e nutação apresentaram comportamento de caráter oscilatório, com tendência ao aumento de suas amplitudes. Para o estudo realizado desprezou-se qualquer presença de atrito ou dissipação de energia.
Ângulos de Euler obtidos por cinemática inversa via ângulos de Euler para as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
A Fig. 4 exibe os resultados para as velocidades angulares nas direções principais de inércia, nota-se que na Fig. 4(a) que é constante ao longo do tempo, o que ocorre pois, se observamos a equação de Euler, Eq. (18), da qual da integração resultam , e , se Nm e , obtém-se que , cujo resultado da integração fornece um valor constante para .
Velocidades angulares na direção principal de inércia , e para as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
A Fig. 5, apresenta as componentes , , e do quatérnion calculado no problema de cinemática inversa utilizando-se os quatérnions. O módulo de , obtido a partir da eq. (7) é igual a 1 para qualquer instante de tempo, garantindo que a operação realizada é apenas de rotação.
Componentes , , e do quatérnion , obtidas para as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
Já para a solução do problema com singularidade, as figuras 6 e 7 expressam os resultados obtidos para os ângulos de precessão , nutação e spin pela cinemática inversa via velocidades angulares instantâneas e por quatérnions, respectivamente. Os resultados expressos nas figuras 6 e 7 deveriam ser idênticos, não fosse a influência das singularidades sob o resultado. A condição inicial do ângulo de nutação foi definida propositalmente com um valor para o qual . A proximidade com um valor onde a eq. (19) é singular faz com que as soluções por meio da cinemática inversa por ângulos de Euler e por quatérnions apresentem resultados diferentes entre si. Não foi utilizado o valor exato devido aos resultados fornecidos serem numericamente indeterminados e impossibilitar a análise comparativa.
Ângulos de Euler obtidos por cinemática inversa via velocidades angulares com as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
Ângulos de Euler obtidos por cinemática inversa via quatérnions com as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
Para os ângulos de precessão e spin, vistos nas figuras 6(a) e 6(c), correspondentes a solução por cinemática inversa através dos ângulos de Euler, é possível notar algumas mudanças bruscas nas trajetórias, enquanto na solução por quatérnions, figuras 7(a) e 7(c), tal comportamento não é percebido. O ângulo de nutação , Fig. 6(c), também apresentou uma mudança repentina em seu comportamento que, na ausência de dissipação de energia ou momentos externos, pode ser facilmente atribuída às singularidades.
O efeito da singularidade no cálculo das velocidades de precessão , nutação e spin, realizado a partir da eq. (19) fica explícito na Fig. 8, onde nota-se que estas velocidades assumem valores muito altos em determinados instantes de tempo, o que ocorre a cada vez que o ângulo de nutação se aproxima muito de °. Esse problema poderia ser parcialmente contornado escolhendo outra sequência de rotações para representar o giroscópio, porém as singularidades seriam apenas “transferidas” para outros valores de ângulos, já que obteríamos um sistema de equações semelhantes à eq. (19) com uma singularidade em , ou seja, para , com , de modo que, utilizando ângulos de Euler para representar as rotações, nunca estaríamos de fato livres do problema. A escolha de outra sequência de rotações seria uma alternativa, por exemplo, na descrição de rotações de um corpo com baixa probabilidade de atingir determinados valores de ângulos por ação de alguma restrição imposta ao problema. Nestes casos, poderíamos admitir singularidades nessas faixas de valores improváveis para tais ângulos. Porém, para o caso onde o sistema pode assumir quaisquer valores de ângulos de nutação, a mudança na sequência das rotações não surtiria grande efeito.
Velocidade de precessão , nutação e spin obtidas por cinemática inversa via velocidades angulares com as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
Outra forma de evitar essas singularidades sem abrir mão do uso dos ângulos de Euler é discutido em [5[5] P. Singla, D. Mortari and J.L. Junkins, in: Proceedings of the AAS Space Flight Mechanics Conference, Hawaii, 2004.], que propõe variar a sequência de rotações adotada de forma dinâmica, alterando a escolha da sequência de rotações entre as doze possíveis quando uma singularidade se aproxima para um determinado valor de ângulo, adotando uma sequência para à qual a singularidade seja deslocada para radianos de distância. Os autores apresentam ainda um algoritmo iterativo para a determinação da sequência de rotações ótima a ser adotada, que requer um número razoável de operações matemáticas. Frente à esta proposta, o uso de quatérnions apresenta um custo computacional muito menor e se mostra muito mais atraente para a solução do problema de cinemática inversa.
Na Fig. 8(b) é possível notar um comportamento inconsistente para a velocidade de nutação, também gerado pela presença de singularidades.
As singularidades não influenciam no cálculo das velocidades angulares nas direções principais de inércia, que são provenientes da integração da eq. (18), como pode ser visto na Fig. 9.
Velocidades angulares na direção principal de inércia , e com as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
Já as componentes , , e do quatérnion são apresentadas na Fig. 10.
Componentes , , e do quatérnion , obtidas para as condições iniciais: °, [rad/s], °, [rad/s], ° e [rad/s].
5. Os Quatérnions Como uma Álgebra de Clifford
Uma característica marcante dos quatérnions é o fato dos mesmos serem uma generalização do conjunto de números complexos . Esta ideia foi estendida e generalizada ainda mais com as ideias de Clifford, introduzidas em 1878, que permitiram trabalhar com elementos de dimensão bem maior. Na verdade, mesmo dimensões menores podem ser descritas álgebras de Clifford. Relações algébricas no espaço Euclidiano podem ser descritas pelos elementos da bases com as propriedades:
Adicionando elementos da forma , sendo , e também dimensões , sendo a dimensão do espaço , pode-se generalizar a álgebra de Clifford.
Tomando um espaço vetorial no corpo dos reais de dimensão com a aplicação bilinear , pode-se definir uma álgebra universal de Clifford2 2 Alguns autores apontam que quando se trabalha no corpo dos reais a álgebra é dita geométrica. com uma álgebra associativa gerada por , satisfazendo a relação fundamental:
Sempre é possível escrever para , sendo e o par de inteiros descrevendo a assinatura da forma quadrática:
sendo uma álgebra de Clifford nomeada de e formada pelos elementos ortogonais entre si, com os elementos com norma e elementos com norma :
sendo números reais e :
e
Assim, pode-se constatar que a álgebra dos reais quando é uma álgebra de Clifford ; já quando , e quando , , mostrando que os quatérnions são uma álgebra de Clifford formada por escalares, vetores e bivetores.
Considerando o espaço vetorial gerado por e , a forma de um elemento qualquer nesta álgebra será descrita por:
Interessante observar que, se for tomado , e , tem-se a relação fundamental da álgebra dos quatérnions , já descrita anteriormente, e o quatérnion , exatamente da forma da eq. (2).
Uma excelente introdução neste tema, acessível em língua portuguesa, é o livro de Jayme Vaz e Roldão da Rocha [37[37] J. Vaz and R. Da Rocha Jr., Álgebras de Clifford e Espinores (Livraria da Física, São Paulo, 2012).].
6. Conclusões
Este artigo exemplificou, a partir de uma aplicação com um giroscópio, como os quatérnions são úteis e interessantes na descrição de rotações e na solução de problemas de cinemática inversa, com diversas vantagens frente a outras opções de parametrização de rotações. Porém, além destes trunfos numéricos, os quatérnions são elementos importantes por formarem a primeira álgebra não-comutativa da história e servirem de porta de entrada para outros sistemas algébricos mais interessantes para descrição matemática de diversos sistemas físicos. Infelizmente, o tema apresentado neste artigo não é, por via de regra, seguido nos cursos de mecânica, tanto para físicos como para engenheiros.
Os autores acreditam ser necessário abordar mais aspectos sobre parametrização de rotação e apresentar os quatérnions não apenas em cursos de mecânica, mas também em disciplinas de álgebra linear e/ou matemática aplicada. Assim, os autores fazem frente aos que defendem a necessidade de um curso de álgebra geométrica no ensino de física, matemática e em cursos de engenharia.
-
1
Um objeto matemático corriqueiro hoje para descrever grandezas físicas, mas que em meados do século XIX não tinha a definição que temos hoje.
-
2
Alguns autores apontam que quando se trabalha no corpo dos reais a álgebra é dita geométrica.
Agradecimentos
O primeiro autor agradece a CAPES por sua bolsa de mestrado. Os autores também agradecem o apoio financeiro do CNPq e da FEPISA, além dos comentários e sugestões do editor e dos revisores da RBEF.
Referências
- [1] L. Euler, Novi Acad. Sci. Petrop 20, 189 (1775).
- [2] J. Stuelpnagel, SIAM review 6, 422 (1964).
- [3] J. Diebel, Matrix 58, 15 (2006).
- [4] J. Neto, Mecânica Newtoniana, Lgrangiana e Hamiltoniana (Editora Livraria da Física, São Paulo, 2004).
- [5] P. Singla, D. Mortari and J.L. Junkins, in: Proceedings of the AAS Space Flight Mechanics Conference, Hawaii, 2004.
- [6] P. Betsch, A. Menzel and E. Stein, Computer Methods in Applied Mechanics and Engineering 155, 273 (1998).
- [7] M. Trindade and R. Sampaio, Journal of the Brazilian Society of Mechanical Sciences 22, 341 (2000).
- [8] J. Kuipers, Quaternions and Rotation Sequences: A Primer with Applications to Orbits, Aerospace, and Virtual Reality (Princeton Paperbacks, Princeton University Press, 1999).
- [9] K. Vogtmann, V. Arnold and A. Weinstein, Mathematical Methods of Classical Mechanics (Graduate Texts in Mathematics) (Springer, New York, 1989).
- [10] H. Goldstein, C. Poole and J. Safko, Classical Mechanics (Addison Wesley, New York, 2002).
- [11] L. Meirovitch, Methods of Analytical Dynamics, Dover Civil and Mechanical Engineering (Dover Publications, Mineola, 2012).
- [12] E.T. Whittaker, A Treatise on the Analytical Dynamics of Particles and Rigid Bodies: With an Introduction to the Problem of Three Bodies (CPU Archive, London, 1970).
- [13] W.R. Hamilton, Elements of Quaternions (Longmans, Green, & Company, Londres, 1866).
- [14] W.R. Hamilton, The Mathematical Papers of Sir William Rowan Hamilton, (CUP Archive, London, 2000).
- [15] W.R. Hamilton, Lectures on Quaternions, (The Royal Irish Academy, Dublin, 1853).
- [16] M. Menon, Revista Brasileira de Ensino de Física 31, 2305 (2009).
- [17] Gustavo Martínez-Sierra y Pierre Francois Benoit Poirier, Latin-American Journal of Physics Education 2, 201 (2008).
- [18] M.J. Amoruso, Euler Angles and Quaternions in Six Degree of Freedom Simulations of Projectiles - Technical Report (DTIC Document, Fort Belvoir, 1996).
- [19] F. Lizarralde and J.T. Wen, IEEE Transactions on Automatic Control 41, 468 (1996).
- [20] M.D. Shuster, Navigation 8, 439 (1993).
- [21] K. Shoemake, Quaternion Calculus and Fast Animation, Computer Animation: 3-D Motion Specification and Control (Siggraph, Anaheim, 1987).
- [22] J.C. Hart, G.K. Francis and L.H. Kauffman, ACM Transactions on Graphics 13, 256 (1994).
- [23] W.R. Rudnicki, G. Bakalarski and B. Lesyng, Journal of Biomolecular Structure and Dynamics 17, 1097 (2000).
- [24] G. Harauz, Ultramicroscopy 33, 209 (1990).
- [25] T. Róg, K. Murzyn, K. Hinsen and G.R. Kneller, Journal of Computational Chemistry 24, 657 (2003).
- [26] W. Sprössig, Advances in Applied Clifford Algebras 17, 259 (2007).
- [27] J.D. Gibbon, D.D. Holm, R.M. Kerr and I. Roulstone, Nonlinearity 19, 1969 (2006).
- [28] H. Eshraghi and J. Gibbon, Journal of Physics A: Mathematical and Theoretical 41, 344004 (2008).
- [29] J. Vaz Jr., Revista Brasileira de Ensino de Física 19, 234 (1997).
- [30] D.H.Z. Carrera, Movimento de Rotação Sem Restrição de um Corpo Rígido Tese de Doutorado, PUC, Rio de Janeiro, 2010.
- [31] A.A.G. Oliveira, Estudo Estatístico dos Processos Envolvidos em uma Plataforma de Atitude Solidária Dissertação de Mestrado, Instituto Nacional de Pesquisas Espaciais, São José dos Campos, 2000.
- [32] S.L. Altmann, Mathematics Magazine 62, 291 (1989).
- [33] R.S. Lima and M.P. Mello, CNMAC - Congresso Nacional de Matemática Aplicada e Computacional (2004).
- [34] V. Pisacane, Fundamentals of Space Systems (Oxford University Press, Oxford, 2005).
- [35] R.A. Wehage, in: Computer Aided Analysis and Optimization of Mechanical System Dynamics (Springer, Berlin, 1984), p. 147-180.
- [36] R.A. Mayo, Journal of Guidance, Control, and Dynamics 2, 44 (1979).
- [37] J. Vaz and R. Da Rocha Jr., Álgebras de Clifford e Espinores (Livraria da Física, São Paulo, 2012).
Datas de Publicação
-
Publicação nesta coleção
Jun 2016
Histórico
-
Recebido
18 Jan 2016 -
Revisado
19 Jan 2016 -
Aceito
30 Jan 2016