Estabilização ativa de sensores em UAVs utilizando rodas de reação (reaction wheels)

Modelagem 3-DOF e simulação de rodas de reação para estabilizar o payload de sensores em VANT de asa fixa: PID adaptativo com anti-windup e testes de viabilidade para fotografia de longa exposição sob turbulência e vibração. A pergunta é se o apontamento nadir permanece estável o bastante para evitar motion blur inaceitável.

Neste post documentamos a dinâmica das rodas, o controlador e os ensaios numéricos de longa exposição — fora das séries A–D, mas relevante para imageamento aéreo.

---

1. O Problema da Estabilização de Payloads em VANTs

Sensores embarcados em aeronaves de asa fixa, como câmeras hiperespectrais ou sistemas LiDAR, sofrem com a dinâmica caótica do voo atmosférico. Turbulências, vibrações de alta frequência provenientes do trem de propulsão (motores e hélices) e manobras ativas de pilotagem induzem variações rápidas de atitude em roll (rolagem), pitch (arfagem) e yaw (guinada).

Para capturas de longa exposição ou imageamento panorâmico de alta precisão, esse jitter angular introduz um desfoque de movimento (motion blur) inaceitável.

O Princípio de Conservação do Momento Angular

Diferente de gimbals motorizados tradicionais, que usam motores nos eixos de articulação para girar a câmera fisicamente em relação à fuselagem (atuadores de torque externo), as rodas de reação atuam como atuadores de torque interno.

Uma roda de reação é composta por um disco inercial simétrico acoplado ao eixo de um motor elétrico brushless de alta precisão. Quando o motor aplica um torque para acelerar a roda em uma direção, o corpo da plataforma recebe um torque de reação de igual magnitude e direção oposta, obedecendo à Terceira Lei de Newton e à conservação do momento angular total do sistema:

$$\mathbf{H}_{\text{total}} = \mathbf{H}_{\text{corpo}} + \mathbf{H}_{\text{rodas}} = \text{constante}$$

Este método de atuação oferece vantagens únicas: é mecanicamente selado, silencioso, livre de folgas de engrenagens e extremamente responsivo a perturbações de alta frequência.

---

2. Modelagem Dinâmica Tridimensional (3-DOF)

Modelamos a plataforma do sensor como um corpo rígido com três graus de liberdade rotacionais atuando no espaço tridimensional.

Dinâmica da Plataforma (Corpo Rígido)

A atitude do corpo é representada pelos ângulos de Euler $\boldsymbol{\theta} = [\phi, \theta, \psi]^T$ (Roll, Pitch, Yaw). A evolução temporal da velocidade angular do corpo $\boldsymbol{\omega} = [\omega_x, \omega_y, \omega_z]^T$ no referencial do corpo é regida pelas Equações de Euler para Corpo Rígido:

$$\mathbf{J} \dot{\boldsymbol{\omega}} + \boldsymbol{\omega} \times (\mathbf{J} \boldsymbol{\omega}) = \boldsymbol{\tau}_{\text{ctrl}} + \boldsymbol{\tau}_{\text{dist}}$$

onde:

• $\mathbf{J} \in \mathbb{R}^{3 \times 3}$ é o tensor de inércia da plataforma do sensor.
• $\boldsymbol{\tau}_{\text{ctrl}}$ é o vetor de torque líquido aplicado pelas rodas de reação sobre o corpo.
• $\boldsymbol{\tau}_{\text{dist}}$ representa os torques externos indesejados (vibrações e rajadas aerodinâmicas).

No nosso projeto, consideramos uma plataforma de sensor com massa de $1{,}0$ kg e tensor de inércia diagonal típico:

$$\mathbf{J} = \begin{bmatrix} 0{,}01 & 0 & 0 \\ 0 & 0{,}01 & 0 \\ 0 & 0 & 0{,}015 \end{bmatrix} \ \mathrm{kg \cdot m^2}$$

Física da Roda de Reação

Para cada roda $i$, o torque aplicado pelo motor $\tau_{\text{motor}, i}$ acelera o disco com momento de inércia $I_{w}$ vencendo o atrito viscoso com coeficiente $b$. A dinâmica da velocidade angular da roda $\omega_{w, i}$ é descrita por:

$$I_{w} \dot{\omega}_{w, i} = \tau_{\text{motor}, i} - b \omega_{w, i}$$

O torque de reação exercido sobre a plataforma é oposto ao torque eletromagnético gerado:

$$\tau_{\text{reacao}, i} = -I_{w} \dot{\omega}_{w, i}$$

Saturação de Momento Angular

A limitação fundamental das rodas de reação é a saturação. Como os motores têm uma velocidade de rotação física máxima $\omega_{w, \text{max}}$, a roda possui uma capacidade limite de armazenar momento angular:

$$H_{\text{max}} = I_{w} \omega_{w, \text{max}}$$

Se um torque de perturbação externo atua de forma contínua em um sentido (por exemplo, vento constante inclinando a asa), a roda precisa acelerar constantemente na direção oposta para compensar. Eventualmente, ela atinge $\omega_{w, \text{max}}$ e satura ($\text{Nível de Saturação} \ge 95\%$). Uma vez saturada, a roda perde a capacidade de aplicar torques naquela direção, exigindo mecanismos de desaturação (como transferir momento para superfícies de controle aerodinâmicas da aeronave ou usar gimbals auxiliares).

---

3. Controle PID Adaptativo com Anti-Windup

Para manter a atitude tridimensional da plataforma travada no setpoint $\boldsymbol{\theta}_{\text{desejado}} = [0, 0, 0]^T$ (sensor apontando na vertical em direção ao centro da Terra), projetamos duas variações de controle de malha fechada de alta taxa (operando de $1000$ Hz a $2000$ Hz).

Controlador PID Clássico com Anti-Windup

O erro de atitude $\mathbf{e}(t) = \boldsymbol{\theta}_{\text{desejado}} - \boldsymbol{\theta}_{\text{atual}}$ é normalizado no intervalo $[-\pi, \pi]$. O sinal de controle de torque calculado para cada eixo é:

$$\boldsymbol{\tau}_{\text{calc}}(t) = K_p \mathbf{e}(t) + K_i \int_{0}^{t} \mathbf{e}(\tau) d\tau + K_d \frac{d\mathbf{e}(t)}{dt}$$

Como o motor elétrico possui limites físicos de torque ($\pm 5$ mN·m no projeto), o sinal enviado é saturado:

$$\boldsymbol{\tau}_{\text{ctrl}}(t) = \text{sat}(\boldsymbol{\tau}_{\text{calc}}(t), \tau_{\text{lim}})$$

Para prevenir o fenômeno de integral windup (onde o termo integral acumula valores massivos enquanto o atuador está fisicamente saturado, causando oscilações destrutivas e instabilidade quando o sistema sai da saturação), foi implementado um loop de realimentação de retrocalculação (back-calculation anti-windup):

$$\mathbf{e}_{\text{sat}}(t) = \boldsymbol{\tau}_{\text{ctrl}}(t) - \boldsymbol{\tau}_{\text{calc}}(t)$$ $$\int \mathbf{e} \, dt \leftarrow \int \mathbf{e} \, dt + K_{\text{aw}} \mathbf{e}_{\text{sat}}(t) \cdot dt$$

Adaptação Dinâmica por Saturação de Momento

Uma inovação do projeto é o Controlador PID Adaptativo. À medida que a velocidade da roda de reação se aproxima do limite de rotação física (saturação de momento angular), a autoridade de controle diminui. Para evitar respostas abruptas que levariam a oscilações violentas (overshoots), o controlador reduz dinamicamente os ganhos com base no nível de saturação instantâneo da roda de cada eixo:

$$\text{sat}_{i} = \frac{|H_i(t)|}{H_{\text{max}}}$$ $$K_i(t) = K_{i, \text{base}} \cdot \max\left(0{,}2, \ 1{,}0 - \gamma \cdot \text{sat}_i\right)$$

onde $\gamma \in [0, 1]$ é a taxa de atenuação. O termo integral é atenuado de forma ainda mais agressiva para mitigar o windup inercial.

---

4. Resultados de Simulação Computacional

As simulações foram executadas no Python utilizando um integrador numérico com passo fixo de simulação robusto para dinâmicas rápidas.

Cenário A: Estabilização Básica sob Turbulência Suave

• Condições Iniciais: Erro angular acentuado: $\text{Roll} = 5{,}0^\circ$, $\text{Pitch} = 3{,}0^\circ$, $\text{Yaw} = -2{,}0^\circ$.
• Perturbações: Vibração estrutural do motor a $50$ Hz ($\pm 0{,}5$ mN·m), rajadas de vento senoidais com período de $2$ s ($\pm 1{,}0$ mN·m), e ruído branco de medição.
• Frequência do Loop: $1000$ Hz ($dt = 1$ ms).

Abaixo estão os resultados estatísticos obtidos após $10$ segundos de simulação:

Sob condições de vento moderado e perturbações ordinárias de voo, as rodas de reação mantêm o sensor alinhado à vertical nadir com erro RMS inferior a $0{,}06^\circ$. O nível de saturação máximo das rodas permaneceu abaixo de $5\%$, demonstrando folga de momento angular e excelente estabilização.

---

Cenário B: Teste Crítico de Viabilidade para Longa Exposição

Para fotografia aérea ou imageamento hiperespectral com tempos de integração longos (por exemplo, tempo de exposição $t_{\text{exp}} = 0{,}5$ s), o jitter angular instantâneo deve ser extremamente reduzido para evitar borrões na imagem.

O Limite Físico de Estabilidade

Definimos que o borrão geométrico causado pela rotação durante o tempo de exposição deve ser menor ou igual a 1 pixel. Para uma lente de distância focal $f = 50$ mm e um tamanho de pixel no detector de $p = 5\ \mu\text{m}$, o limite de jitter angular aceitável é dado por:

$$\theta_{\text{lim}} = \frac{p}{f} = \frac{5 \times 10^{-6}\ \text{m}}{50 \times 10^{-3}\ \text{m}} = 10^{-4}\ \text{rad}$$

Convertendo para graus e segundos de arco:

$$\theta_{\text{lim}} = 10^{-4} \times \frac{180}{\pi} \approx 0{,}00573^\circ \approx \mathbf{20{,}63\ \text{arcsec}}$$

Resultados sob Perturbação Agressiva

Simulamos um voo com turbulências atmosféricas pesadas e manobras ativas da aeronave (oscilações de pitch lentas de $\approx 2$ mN·m induzidas por curvas e correção de altitude) em um loop de $2000$ Hz.

O resultado indica que o desvio RMS instantâneo ($\approx 150 - 340$ arcsec) excede o limite estrito de blur de 1 pixel por um fator de até $17\times$. O jitter máximo ultrapassa o limite por mais de $27\times$. A taxa de sucesso para capturas viáveis com tempo de exposição de $0{,}5$ segundos sob este nível de perturbação foi de $0\%$.

---

5. Comparação e Discussão de Atuadores

A análise dinâmica e os dados simulados trazem importantes trade-offs de engenharia para o projeto do payload:

Conclusão e Diretrizes de Projeto

1. Reaction Wheels Isoladas: São excelentes estabilizadores inerciais contra vibrações de motores e jitter de média frequência em voos calmos. Porém, são insuficientes para compensar rajadas constantes ou manobras da asa durante longas exposições fotográficas, pois saturam rapidamente ou exigem motores e discos massivos e ineficientes em consumo.
2. Solução Híbrida Recomendada: Para atender ao rigoroso requisito de $20{,}63\ \text{arcsec}$ sob turbulência, a arquitetura ideal consiste em um Gimbal motorizado ativo de 2 ou 3 eixos (para compensar os movimentos de baixa frequência e grande amplitude da aeronave) acoplado a um pod estabilizado internamente por pequenas reaction wheels de alta frequência (para isolar os harmônicos de vibração do motor e vento rápido).

---

Referências

• Hughes, P. C. (2012). Spacecraft Attitude Dynamics. Courier Corporation.
• Franklin, G. F., Powell, J. D., Emami-Naeini, A. (2014). Feedback Control of Dynamic Systems. 7th ed. Pearson.
• Spong, M. W., Hutchinson, S., Vidyasagar, M. (2005). Robot Modeling and Control. Wiley.