Equivalent Mechanical Models for Sloshing

Resumo Técnico: Modelos Mecânicos Equivalentes para Sloshing (Oscilação de Propelente)

Definição do Problema
O sloshing de propelente permanece um desafio crítico para a estabilidade e o controle de espaçonaves, particularmente em dois regimes distintos: ambientes de microgravidade que exigem alta estabilidade de apontamento (ex: missões interplanetárias) e condições de alta aceleração ($g$) encontradas durante lançamento, reentrada e pouso. Embora modelos baseados em pêndulos e massa-mola-amortecedor sejam ferramentas padrão para analistas de controle avaliarem as perturbações induzidas pelo sloshing, uma formulação matemática rigorosa que conecte essas representações, especificamente contabilizando os efeitos de forças longitudinais nominais em regimes de alto-$g$, é necessária para um projeto de orientação, navegação e controle (GNC) robusto.

Metodologia
O artigo apresenta uma derivação rigorosa de modelos mecânicos equivalentes para sloshing, progredindo de dinâmicas não lineares para representações linearizadas adequadas para análise de controle.

1. Formulação Não Linear:

   • O estudo começa com um pêndulo esférico único de massa $m_P$ e comprimento $\ell$ acoplado a um corpo rígido $B$ de massa $m_B$.
   • Utilizando uma abordagem Newtoniana, os autores derivam as equações de movimento acopladas para o sistema de 8 graus de liberdade (translação/rotação do corpo rígido + ângulos do pêndulo).
   • A derivação lida explicitamente com a força de tensão $N$ e a elimina para reduzir as variáveis do sistema, resultando em uma forma matricial compacta $\mathcal{M}\ddot{\mathbf{x}} = \mathbf{u}$.
   • A formulação é estendida para acomodar múltiplos pêndulos, representando múltiplos tanques ou modos de sloshing, utilizando uma matriz de massa estruturada em blocos para explorar a esparsidade visando uma solução numérica eficiente.
   • Um modelo de amortecimento viscoso é incorporado pela adição de termos proporcionais às taxas angulares, que afetam a dinâmica do pêndulo sem alterar diretamente as equações do corpo rígido.

2. Linearização e Análise Modal:

   • Reconhecendo que os sistemas de controle embarcados tipicamente estimam o centro de massa assumindo uma distribuição de propelente rígida, os autores reformulam a dinâmica em relação a um centro de massa global nominal ($\bar{G}$) em vez do centro de massa do corpo rígido.
   • O sistema é linearizado em torno de um equilíbrio caracterizado por uma atitude constante alinhada com o referencial inercial e uma força longitudinal nominal prescrita $\bar{F}_z$.
   • Dois casos são analisados: forças definidas no referencial inercial e forças definidas no referencial do corpo fixo.
   • Através de transformações de coordenadas e manipulação algébrica, a dinâmica linearizada é convertida em uma "forma modal". Isso envolve a definição de novas coordenadas modais ($\eta_\theta, \eta_\phi$) que desacoplam a dinâmica do pêndulo em dois osciladores planares independentes.
   • O sistema resultante é expresso como um sistema linear MIMO de segunda ordem ($M\delta\ddot{\mathbf{x}} + K\delta\mathbf{x} = \delta\mathbf{u}_{ext}$), identificando explicitamente a matriz de participação modal e os termos de rigidez induzidos pela força nominal.

3. Equivalência com Modelos de Massa-Mola-Amortecedor:

   • O artigo deriva as equações linearizadas para um sistema de massa-mola-amortecedor onde a massa de sloshing é restrita a mover-se lateralmente via molas de rigidez $k$.
   • Ao definir a rigidez da mola como $k = m_P \omega_0^2$ (onde $\omega_0$ é a frequência natural derivada do modelo de pêndulo), os autores demonstram que a dinâmica linearizada de massa-mola-amortecedor é matematicamente idêntica à dinâmica linearizada do pêndulo.
   • A análise destaca que a presença de uma força longitudinal nominal $\bar{F}_z$ é essencial para que o pêndulo exiba comportamento oscilatório análogo a uma mola; na ausência desta força ($\bar{F}_z=0$), o modelo de pêndulo requer modificação (ex: adicionando um terceiro grau de liberdade modal) para representar o sloshing em 0-g.

Principais Resultados

• Equivalência Matemática: O estudo estabelece formalmente a equivalência exata entre o modelo de pêndulo linearizado e o modelo clássico de massa-mola-amortecedor sob condições de alto-$g$, esclarecendo como a força nominal atua como o mecanismo de restauração na formulação do pêndulo.
• Representação Modal: A derivação produz uma representação modal limpa onde os modos de sloshing são desacoplados da dinâmica do corpo rígido na matriz de rigidez, com o acoplamento ocorrendo apenas na matriz de massa via matriz de participação modal.
• Validação:
  • Validação Não Linear: O modelo não linear foi validado contra uma implementação MATLAB Simscape Multibody (cubo rígido com quatro pêndulos). Simulações no domínio do tempo mostraram excelente concordância, com erros máximos abaixo de 0,002% da amplitude da resposta, com discrepâncias atribuídas unicamente a diferenças nos esquemas de integração numérica.
  • Validação Linear: Diagramas de Bode comparando o modelo linearizado analítico com o plano Simscape linearizado demonstraram excelente concordância em toda a faixa de frequência para ambos os inputs de força em referencial inercial e do corpo. Pequenas discrepâncias de baixa frequência no caso do referencial inercial foram atribuídas aos limites de precisão numérica do processo de linearização do Simscape.

Significância e Alegações
O artigo afirma fornecer uma base matemática rigorosa para modelar o sloshing de propelente em ambientes de alta aceleração. Sua principal contribuição é a derivação explícita das equações de movimento não lineares para um sistema de múltiplos pêndulos e a subsequente linearização que contabiliza forças longitudinais nominais. Ao provar a equivalência matemática entre os modelos de pêndulo e de massa-mola-amortecedor, o trabalho valida o uso deste último para o projeto de controle, fundamentando-o na dinâmica de pêndulo mais fundamental. Os autores enfatizam que os modelos propostos são bem postos para aplicações realistas de espaçonaves onde os ângulos de sloshing permanecem dentro de limites razoáveis, e a validação confirma a correção das formulações derivadas tanto para simulação no domínio do tempo quanto para análise no domínio da frequência.