Universidade de Brasília – UnB

Faculdade UnB Gama – FGA

Engenharia Aeroespacial

Estudo de Sistema de Controle de Atitude para

Simulador de Nanossatélite

Autor: Maria Cristina Mendes Faustino

Orientador: Dr. William Reis Silva

Brasília, DF
2019
 Maria Cristina Mendes Faustino

Estudo de Sistema de Controle de Atitude para

Simulador de Nanossatélite

Monografia submetida ao curso de graduação

em Engenharia Aeroespacial da Universidade
de Brasília, como requisito parcial para ob-
tenção do Título de Bacharel em Engenharia
Aeroespacial.

Universidade de Brasília – UnB

Faculdade UnB Gama – FGA

Orientador: Dr. William Reis Silva

Coorientador: Dr. Renato Alves Borges

Brasília, DF
2019
Maria Cristina Mendes Faustino
Estudo de Sistema de Controle de Atitude para Simulador de Nanossatélite/
Maria Cristina Mendes Faustino. – Brasília, DF, 2019-
61 p. : il. (algumas color.) ; 30 cm.

Orientador: Dr. William Reis Silva

Trabalho de Conclusão de Curso – Universidade de Brasília – UnB

Faculdade UnB Gama – FGA , 2019.
1. Controle de Atitude. 2. Regulador Quadrático Linear. I. Dr. William Reis
Silva. II. Universidade de Brasília. III. Faculdade UnB Gama. IV. Estudo de
Sistema de Controle de Atitude para Simulador de Nanossatélite

CDU 02:141:005.6
 Maria Cristina Mendes Faustino

Estudo de Sistema de Controle de Atitude para

Simulador de Nanossatélite

Monografia submetida ao curso de graduação

em Engenharia Aeroespacial da Universidade
de Brasília, como requisito parcial para ob-
tenção do Título de Bacharel em Engenharia
Aeroespacial.

Trabalho aprovado. Brasília, DF, 03 de Julho de 2019 – Data da aprovação do

trabalho:

Dr. William Reis Silva

Orientador

Dr. Thiago Felippe Kurudez Cordeiro

Convidado 1

Dr. Manuel Nascimento Dias Barcelos

Júnior
Convidado 2

Brasília, DF
2019
Dedico este trabalho aos meu pais Eliana Mendes e Manoel Carlito Alves Faustino, que
sempre me apoiaram, incentivaram e nunca mediram esforços para que eu completasse
minha trajetória na Universidade de Brasília
 Agradecimentos

Agradeço primeiramente a Deus por toda força, ânimo e coragem que me ofereceu
para ter alcançado minha meta.
Agradeço aos meus pais, Eliana e Carlito, cujo o apoio e incentivo foram essencial
para que eu alcançasse meu sonho. Agradeço também à minha irmã Cristiane Mendes
Faustino, por sempre estar ao meu lado, me aconselhando, ouvindo e dando força.
Agradeço ao Guilherme Gonzaga Elizeu, por estar ao meu lado nos 3 últimos
anos de graduação me dando força, me ouvindo, apoiando e ajudando de todas as formas
possíveis.
Agradeço aos professores do curso de graduação em Engenharia Aeroespacial da
Faculdade do Gama por todo conhecimento transmitido ao longo de todos esses anos.
Agradeço ao Profº Dr. William Reis Silva e Profº Dr. Renato Alves Borges pela ori-
entação e ajuda neste trabalho, sem elas, ele não poderia ter sido desenvolvido. Agradeço
também ao Profº Dr. Renato Alves Borges por ter me dado a oportunidade de desenvolver
esse trabalho dentro do Laboratório de Simulação e Controle de Sistemas Aeroespaciais -
LODESTAR, cuja infraestrutura e suporte disponibilizados por ele foram essenciais para
o trabalho.
Agradeço aos alunos Nikson Bernardes Ferreira e Luan Haickel Araújo que ajuda-
ram e compartilharam seus conhecimentos, isso foi extremamente enriquecedor para este
trabalho.
Por fim, agradeço a todos que contribuíram direta ou indiretamente no desenvol-
vimento deste trabalho, muito obrigada!
 “Prefiram a minha instrução à prata,
e o conhecimento ao ouro puro, pois a
sabedoria é mais preciosa do que rubis;
Nada do que vocês possam desejar compara-se a ela.“
(Bíblia Sagrada, Provérbios 8, 10-11)
 Resumo
O apontamento e rotação do satélite é feito pelo, do inglês, Attitude Determination and
Control System (ADCS), ele determina a atitude atual do corpo e altera para o alvo de
desejo. Para que isso funcionem bem após o lançamento do satélite é necessário realizar
um estudo em ambiente simulado. Para tal, plataformas de teste são construídas e colo-
cadas em cima de mancais a ar, de forma a criar um ambiente mais próximo daquele do
espaço. Visando promover esses estudos em terra, o Laboratório de Simulação e Controle
de Sistemas Aeroespaciais (LODESTAR) conta com uma dessas plataformas em suas ins-
talações. Esse trabalho tem como objetivo o estudo e a validação do sistema de controle
em 1-eixo, esse estudo irá contribuir para o aperfeiçoamento do simulador bem como po-
der fornecer dados para futuras missões. No estudo foi utilizado um mockup 1U atuado
por uma roda de reação no eixo Z, e o Linear Quadratic Regulator (LQR) foi utilizado
para efetuação do controle no sistema. Uma análise numérica do problema foi feita e com
os resultados obtidos verificou-se necessidade de se fazer a identificação, a planta do sis-
tema foi modelada e com os valores obtidos por meio de software a análise experimental
foi realizada. Os resultados das duas análises foram comparadas. Foi demonstrado aqui a
ação de um PI sobre um LQR e a sua aplicação em um sistema real.

Palavras-chave: Controle de Atitude. Regulador Quadrático Linear. Simulador de con-

trole de atitude.
 Abstract
The satellite’s pointing and rotation is done by the Attitude Determination and Control
System (ADCS), it determines the current attitude of the body and changes to the target
of desire. For this to work well after the launch of the satellite is necessary to carry out a
study in simulated environment. To do this, test stands are built and placed on top of air
bearings, in order to create an environment closer to that of space. Aiming to promote
these studies on land, the Laboratory of Simulation and Control of Aerospace Systems
(LODESTAR) has one of these platforms in its facilities. This study aims to study and
validate the 1-axis control system, this study will contribute to the improvement of the
simulator as well as being able to provide data for future missions. In the study, a 1U
mockup was used by a reaction wheel in the Z axis, and the Linear Quadratic Regulator
(LQR) was used to carry out the control in the system.A numerical analysis of the problem
was made and with the results obtained it was necessary to make the identification, the
system plan was modeled and with the values obtained through software the experimental
analysis was performed. The results of the two numerical analyzes were compared. It has
been demonstrated here the action of an IP on an LQR and its application in a real
system.

Key-words: Attitude Control. Linear Quadratic Regulator. Attitude control simulator.

 Lista de ilustrações

Figura 1 – Mesa de rolamento a ar do LODESTAR Fonte:(SILVA, 2016) . . . . . . . 17

Figura 2 – Componentes do mancal a ar. Fonte:(SILVA; RODRIGUES, 2015) . . . . . . 18
Figura 3 – Gaiola de Helmholtz. Fonte:(PLOEG, 2015) . . . . . . . . . . . . . . . . . 19
Figura 4 – Sistema Equatorial Terrestre. Fonte: (SILVA, 2016) . . . . . . . . . . . . 22
Figura 5 – Sistema Orbital. Fonte:(SILVA, 2016) . . . . . . . . . . . . . . . . . . . . 23
Figura 6 – Sistema Orbital Local. Fonte:(SILVA, 2016) . . . . . . . . . . . . . . . . . 24
Figura 7 – Sistema Orbital Local visualizado no plano orbital. Fonte:(SILVA, 2016) . 25
Figura 8 – Representação do Sistema de Atitude em verde. Fonte:(SILVA, 2016) . . 26
Figura 9 – Representação da sequência de rotação adotada (𝜑,𝜃,𝜓). Fonte: Modificada
de (SILVA, 2016) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
Figura 10 – CAD Mockup 1U - Vista Externa. Fonte: Autora . . . . . . . . . . . . . 35
Figura 11 – CAD Mockup 1U - Estrutura interna. Fonte: Autora . . . . . . . . . . . 35
Figura 12 – Quadro com os dados de estimação dos momentos de inércia. Fonte:
Autora . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
Figura 13 – Mockup finalizado. Fonte: Autora . . . . . . . . . . . . . . . . . . . . . . 36
Figura 14 – Interior do Mockup. Fonte: Autora . . . . . . . . . . . . . . . . . . . . . 37
Figura 15 – Esquemático de ligação do motor. Fonte: (SILVA, 2013) . . . . . . . . . . 38
Figura 16 – Ângulos Yaw, Pitch e Roll sendo controlados. Fonte:Autor . . . . . . . 40
Figura 17 – Velocidade angular. Fonte:Autor . . . . . . . . . . . . . . . . . . . . . . 40
Figura 18 – Torque das rodas de reação. Fonte:Autor . . . . . . . . . . . . . . . . . 41
Figura 19 – Gráfico da resposta ao degrau (Ângulo). Fonte: Autora . . . . . . . . . . 42
Figura 20 – Gráfico da resposta ao degrau (Velocidade Ângular). Fonte: Autora . . . 42
Figura 21 – Gráfico do Torque necessário. Fonte: Autor . . . . . . . . . . . . . . . . 43
Figura 22 – Planta do sistema real. Fonte: Autor . . . . . . . . . . . . . . . . . . . . 44
Figura 23 – Montagem do experimento. Fonte: Autor . . . . . . . . . . . . . . . . . 45
Figura 24 – Resultado da simulação. Fonte: Autor . . . . . . . . . . . . . . . . . . . 45
Figura 25 – Dados da simulação.(Zoom) Fonte: Autor . . . . . . . . . . . . . . . . . 46
Figura 26 – Percentual de tensão aplicada ao motor ao utilizar sistema de controle
baseado em PI + LQR. Fonte: Autor . . . . . . . . . . . . . . . . . . . . 46
Figura 27 – Mockup com atuadores nos 3 eixos. Fonte: Autora . . . . . . . . . . . . . 49
Figura 28 – Montagem do exterior do Mockup. Fonte: Autora . . . . . . . . . . . . . 57
Figura 29 – Montagem do Interior do Mockup. Fonte: Autora . . . . . . . . . . . . . 58
Figura 30 – Desenho com as especificações do CubeSat 1U. Fonte:(SLO, 2014) . . . . 61
 Lista de tabelas

Tabela 1 – Características técnicas da Roda de Reação . . . . . . . . . . . . . . . 38

Tabela 2 – Parâmetros da simulação . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Tabela 3 – Pesos de Q e R . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
Tabela 4 – Ganhos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
 Lista de abreviaturas e siglas

ADCS Attitude Determination and Control System (Sistema de Determinação

e Controle de Atitude)

INPE Instituto Nacional de Pesquisas Espaciais

LODESTAR Laboratório de Simulação e Controle de Sistemas Aeroespaciais

LQR Linear Quadratic Regulator (Regulador Quadrático Linear)

PI Proporcional Integral (Controlador)

PLA Poliácido Láctico

 Lista de símbolos

O𝑒𝑞 X𝑒𝑞 Y𝑒𝑞 Z𝑎𝑞 Sistema Equatorial Terretre

O𝑜𝑝 x𝑜𝑝 y𝑜𝑝 z𝑜𝑝 Sistema Orbital

O𝑜 x 𝑜 y 𝑜 z 𝑜 Sistema Orbital Local

Oxxz Sistema de Atitude

𝜃 Ângulo de picth

𝜑 Ângulo de roll

𝜓 Ângulo de yaw

A321 Matriz de rotação representada por ângulos de Euler

h Vetor do momento angular

ℎ𝐼 Vetor do momento angular em torno do centro de massa do referencial

inercial

ℎ𝐵 Vetor do momento angular em torno do centro de massa do corpo

𝜔 Vetor da velocidade angular do satélite

𝜔˙ Vetor da aceleração angular do satélite

𝜔𝑥 Matriz antissimétrica

Ω Vetor de velocidade angular da roda de reação

Ω̇ Vetor da aceleração angular da roda de reação

I Vetor do momento de inércia do satélite

𝐼𝑊 Vetor do momento de inércia da Roda de reação

N Torques externos

𝑥 Vetor de estados

𝐴 Matriz que representa a dinâmica do sistema

𝐵 Matriz que representa como o sinal de controle afeta os estados

𝐶 Matriz que representa os estados que estão sendo observados

𝐷 Matriz que representa os sinais de controle que afetam diretamente a
saída

𝑄 Matriz que representa o custo da precisão de posicionamento

𝑅 Matriz que representa o custo da energia gasta pelos atuadores

𝐽 Funcional

𝑢 Sinal de controle

𝑢* Sinal de controle ótimo

𝑡 Tempo

𝐻 Hamiltoniana

𝐾 Ganho
 Sumário

1 INTRODUÇÃO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.1 Contextualização . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.2 Motivação . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.3 Justificativa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.4 Objetivo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.4.1 Objetivo Geral . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.4.2 Objetivos Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.5 Metodologia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
1.6 Apresentação do Trabalho . . . . . . . . . . . . . . . . . . . . . . . . . 21

2 FUNDAMENTAÇÃO TEÓRICA . . . . . . . . . . . . . . . . . . . . 22
2.1 Sistema de Referência . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.1.1 Sistema de Referência Equatorial Terreste . . . . . . . . . . . . . . . . . . 22
2.1.2 Sistema de Referência Orbital . . . . . . . . . . . . . . . . . . . . . . . . 23
2.1.3 Sistema de Referência Orbital Local . . . . . . . . . . . . . . . . . . . . . 24
2.1.4 Sistema de Referência de Atitude . . . . . . . . . . . . . . . . . . . . . . 26
2.2 Atitude de Satélite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.2.1 Representação da Atitude por Ângulos de Euler . . . . . . . . . . . . . . . 27
2.2.2 Dinâmica do Satélite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.2.3 Cinemática do Satélite . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
2.3 Regulador Linear Quadrático (LQR) . . . . . . . . . . . . . . . . . . . 29
2.3.1 Definição das Matrizes em Espaço de Estados . . . . . . . . . . . . . . . . 32

3 PROJETO DO MOCKUP 1U . . . . . . . . . . . . . . . . . . . . . 34
3.1 Padrão CubeSat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.2 Modelagem 3D . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
3.3 Mockup 1U - LODESTAR . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.3.1 Roda de Reação . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.3.2 BeagleBone Blue . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

4 ANÁLISE NUMÉRICA . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.1 Análise em 3-eixos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.2 Análise em 1-eixo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41

5 ANÁLISE EXPERIMENTAL . . . . . . . . . . . . . . . . . . . . . . 44
5.1 Modelando o Sistema . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
6 CONCLUSÃO . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48

7 TRABALHOS FUTUROS . . . . . . . . . . . . . . . . . . . . . . . . 49

REFERÊNCIAS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50

APÊNDICES 52

APÊNDICE A – EQUAÇÃO DA CINEMÁTICA EM ÂNGULOS DE

EULER . . . . . . . . . . . . . . . . . . . . . . . . 53

APÊNDICE B – MATRIZES DO ESPAÇO DE ESTADO . . . . . . 55

APÊNDICE C – MONTAGEM DO MOCKUP . . . . . . . . . . . . 57

APÊNDICE D – CÓDIGO LQR GENÉRICO . . . . . . . . . . . . . 59

ANEXOS 60

ANEXO A – ESPECIFICAÇÕES DO CUBESAT 1U . . . . . . . . . 61

 16

1 Introdução

A descoberta de que foguetes de propulsão química poderiam levar carga útil até
o espaço foi inicialmente utilizada durante a Segunda Guerra Mundial para fins bélicos.
Porém, foi após a guerra que o potencial do foguete para a exploração do espaço foi vista.
Com isso, o início da Era Espacial foi marcado com o lançamento do primeiro satélite
artificial da Terra, Sputnik 1. Desde então os satélites artificiais tornaram-se mecanismos
indispensáveis.
Devido aos grandes avanços tecnológicos a redução dos circuitos embarcados tornou-
se possível, isso proporcionou o surgimento dos nanossatélites. Esta classe de satélites vem
se destacando pelo seu baixo custo e facilidade de produção. Cada vez que um satélite de
grande porte falha os prejuízos econômicos gerados são grandes. Além dos benefícios cita-
dos anteriormente, os nanossatélites podem ser lançados juntamente com outros, grandes
ou pequenos, o que acaba diminuindo o custo de envio. Algumas das suas aplicações são
coletas de dados, teste de tecnologias para posterior uso em satélites maiores, observação
da Terra e acadêmicas.
Durante a construção de um satélite é necessário fazer um estudo, projeto e por
fim a validação dos sistemas nele embarcados com a finalidade de garantir o pleno funcio-
namento após ser lançado. Um desses sistemas é o Sistema de Determinação e Controle da
Atitude (Attitude Determination and Control System - ADCS), este é quem irá garantir
a orientação do satélite em órbita. O controle de atitude é fundamental para a realização
de missões que estão relacionadas a apontamento de equipamentos, como por exemplo
câmeras, painéis solares, antenas e também está relacionada à própria estabilização do
satélite.
Neste trabalho será apresentado um estudo do sistema de controle de atitude,
onde será demonstrado todo o equacionamento da dinâmica e cinemática de atitude para
um CubeSat atuado por roda de reação, em 1 eixo. Para então projetar o controlador
utilizando o Regulador Quadrático Linear (LQR), o estudo será realizados de duas formas:
por meio de uma análise numérica e Análise experimental. A validação do sistema de
controle se dará por meio de simulações em um mockup tamanho 1U (corresponde a
1 unidade de um CubeSat, um cubo de 10 cm de aresta). Os testes foram realizados no
Laboratório de Simulação e Controle de Sistemas Aeroespaciais (LODESTAR), localizado
na Universidade de Brasília.
Capítulo 1. Introdução 17

1.1 Contextualização
O ADCS é o sistema, cuja função é regular o movimento rotacional e o apon-
tamento do satélite. Ele determina a orientação relativa a um corpo ou a um ponto de
interesse, ou seja, determina a atitude atual do corpo e muda para o alvo de desejo.
Para que todo esse sistema funcione com precisão, uma diversidade tecnológica se
faz necessária tanto para a atuação quanto para a detecção. Sensores, como rastreadores
de estrela, Sol e Terra podem definir a orientação colhendo dados de imagem combinando
as posições do satélite e padrões do mapa, é um exemplo. As tecnologias que envolvem
apontamento podem ser citados os magnetorquers e rodas de reação. O primeiro são
eletroímãs que orientam os satélites em relação ao campo magnético da Terra, já o segundo
atuador cria um movimento oposto ao do satélite na direção do giro (RUSSEL; STRAUB,
2017).
O LODESTAR promove contribuição nas atividades de pesquisa e desenvolvimento
no setor aeroespacial na região Centro-Oeste do país, também estabelecendo cooperações
nacionais e internacionais com institutos e empresas nesse domínio (ISHIOKA, 2017). O
Laboratório possui uma plataforma de simulação, na qual é possível estudar e analisar
sistemas de controle de atitude para nanossatélites.
Esse simulador trata-se de um airbearing-platform de macal esférico, que simula
um ambiente espacial com ausência de torques gravitacionais. A plataforma, mostrada na
Fig. 1, possui total liberdade de rotação no eixo Z, porém é limitada em 48° nos eixos X
e Y.

Figura 1 – Mesa de rolamento a ar do LODESTAR Fonte:(SILVA, 2016)

A simulação de um ambiente de microgravidade é praticável devido à existência

de um colchão de ar localizado entre o mancal da base e a esfera da plataforma. O ar-
Capítulo 1. Introdução 18

comprimido chega na base do mancal esférico por meio de seis furos igualmente espaçados,
Fig. 2. A mesa de testes é instrumentada com uma Unidade de Mensuração Inercial (IMU)
de 9 graus de liberdade. Além disso, são embarcados na bancada um microcontrolador,
o rádio XBee para comunicação sem fio, drivers de motores, motores de passo para o
movimento translacional de massas em três graus de liberdade não redundantes e baterias.

(a) Base do Mancal (b) Semiesfera do mancal

Figura 2 – Componentes do mancal a ar. Fonte:(SILVA; RODRIGUES, 2015)

A bancada de teste encontra-se instalada no interior de uma gaiola de Helmholtz

(Fig. 3), isso possibilita fazer simulação do campo magnético em condições semelhantes
às orbitais. A gaiola é formada por um conjunto de seis bobinas quadradas, duas em cada
eixo da gaiola, utilizadas para gerar um campo magnético homogêneo de acordo com a
corrente elétrica aplicada(ISHIOKA, 2017).
Capítulo 1. Introdução 19

Figura 3 – Gaiola de Helmholtz. Fonte:(PLOEG, 2015)

A gaiola tem a capacidade de simular variações de campo magnético em magnitude

e direção com intensidades de ±1,5 Gauss (G) nos três eixos. É importante ressaltar que a
uniformidade do campo magnético no centro da gaiola foi verificada e garante a realização
de testes com atuador magnético para o volume ocupado pela mesa (PLOEG, 2015).
De modo a colaborar com os avanços nas pesquisas na área de simulação de siste-
mas de controle de atitude de nanossatélite dentro do LODESTAR. Há uma necessidade
de projeto e construção de um mockup para efetuação dos testes com as rodas de reação
no simulador, e assim colaborar para o aperfeiçoamento da bancada de testes.

1.2 Motivação
Avanços na área da comunicação, observação da Terra, espaço e acontecimentos
climáticos ocorreram graças ao progresso que ocorreu na área Aeroespacial e tecnológica
no geral. Contudo, este conhecimento não está mais apenas restrito à agencias espaciais
e grandes empresas. Uma vez que, surgiram os nanossatélites o que barateou o processo
possibilitando assim que um maior público pudesse participar de uma missão espacial.
Universidades e até mesmo escolas vem investindo recursos e tempo em desenvolvimento
e pesquisa de Cubesat e plataformas de testes para sistemas aeroespaciais.
Seguindo este viés, o LODESTAR da Universidade de Brasília procura incentivar
o desenvolvimento e a pesquisa de tecnologias espaciais. Ele é equipado para simular o
Capítulo 1. Introdução 20

ambiente espacial para efetivação do estudo do controle. Este laboratório proporciona

inclusão da comunidade acadêmica neste meio promovendo muitos avanços em pesquisas,
trocas de conhecimentos e experiências.
Este trabalho foi motivado a fim de colaborar com os estudos na área de controle
de atitude a serem realizados no simulador do LODESTAR, e também ajudando assim a
aperfeiçoá-lo, para então um dia seu possível uso na validação de missões espaciais reais.

1.3 Justificativa
Quando um satélite é lançado espera-se que seus subsistemas tenham um funcio-
namento perfeito. Caso qualquer um desses venham a falhar, não há possibilidade de se
fazer a manutenção pondo um ponto final na missão.
Dentre todos os subsistemas que estão envolvido no projeto de um nanossatélite,
aqui será destacado o sistema de determinação e controle de atitude. Este é essencial para
a estabilização do satélite após o seu desacoplamento do lançador, e outras aplicações
citadas anteriormente. Um sistema como este deve ser muito bem projetado e pesquisado.
Pensado nisso, plataformas de teste são desenvolvidas. Estas são capazes de possi-
bilitar o estudo aprofundado do controle de atitude do satélite. Estudos como este possi-
bilitam o aperfeiçoamento dos simuladores, bem como validam o controlador fornecendo
dados para possíveis missões futuras.

1.4 Objetivo
1.4.1 Objetivo Geral
Estudo e validação do sistema de controle de atitude em um eixo de um CubeSat
1U atuado por roda de reação para realização de apontamento e frenagem.

1.4.2 Objetivos Específicos

• Fazer um estudo sobre sistema de referência, equações da dinâmica e cinemática de

atitude e definir o modelo matemático do sistema de 1 eixo.

• Projetar o controlador onde serão definidas as matrizes de ganho Q e R do LQR

que satisfarão os requisitos.

• Realizar simulações numéricas para verificação do projeto do simulador.

• Projeto e fabricação da estrutura de um mockup tamanho 1U.

• Implementação do controlador no ambiente simulado para ser feita a validação.

Capítulo 1. Introdução 21

1.5 Metodologia
Considerando todos os objetivos específicos, observa-se uma divisão do trabalho
em três partes: Análise numérica, projeto estrutural e simulação experimental. Para a
realização da análise numérica primeiramente foi necessário conhecer a atitude de um
satélite, que possui roda de reação como atuador, foi feito um estudo das equações da
dinâmica e cinemática.
Após, as matrizes da dinâmica do sistema e do sinal de controle que afetam o estado
foram geradas baseando-se nas equações citadas anteriormente. O software MATLAB foi
utilizado para realizar a análise nesta etapa.
Antes de prosseguir com os testes experimentais, houve a necessidade da constru-
ção de um mockup tamanho 1U. Este foi modelado no CATIA e construído por meio de
manufatura aditiva. Feito isso, as simulações foram realizadas podendo analisar o com-
portamento do controle de atitude em ambiente simulado.

1.6 Apresentação do Trabalho

Este trabalho está organizado em sete capítulos, incluindo este capítulo introdu-
tório.
No Capítulo 2 alguns conceitos básicos são demonstrados. Estes conceitos são
referentes à matemática necessária para a definição formal da atitude, sistema referencial
e ao Regulador Quadrático Linear.
No Capítulo 3 o projeto do mockup 1U utilizado para as simulações é apresentado.
É mostrado desde o processo de modelagem até a sua completa montagem com toda a
eletrônica embarcada.
No Capítulo 4 a análise numérica para o controle de atitude em um eixo é feito.
Nele será apresentada a modelagem matemática para o sistema de um eixo e alguns
gráficos gerados via software.
No Capítulo 5 o teste da simulação do controle de atitude realizado com o mockup
será apresentado, assim como todos os resultados das simulações.
No Capítulo 6 a conclusão obtida a partir dos resultados.
No Capítulo 7 uma sugestão de possível continuação da pesquisa para avanços nos
resultados.
 22

2 Fundamentação Teórica

Neste capítulo será feita uma breve fundamentação teórica referente a sistema de
referência, a dinâmica e controle de atitude e ao Regulador Linear Quadrático (LQR).

2.1 Sistema de Referência

Um sistema de referência é sempre definido por uma origem 𝑂𝑖 e por um triedro
de vetores ortogonais. Segundo (FEHSE, 2003) , três tipos de sistemas de coordenadas
são necessários:

• Sistema de referência orbital: Descreve a orientação da órbita no espaço inercial,

por exemplo, a Terra. Descrevendo o movimento do corpo em sua órbita.

• Sistema de referência orbital local do satélite: Descreve o movimento relativo

para um ponto particular na órbita ou para outro corpo.

• Sistema de referência da atitude do satélite: Descreve a dinâmica e cinemática

de atitude do corpo com relação ao seu centro de massa

2.1.1 Sistema de Referência Equatorial Terreste

Coordenada usada para descrever o movimento orbital ao redor do centro da Terra
com direções fixas e inerciais, representado na Fig. 4. A Terra é admitida como esférica,
ou seja, o centro de massa e o centro geométrico são coincidentes e por conseguinte, o
ponto focal do movimento orbital. (SILVA, 2016)

Figura 4 – Sistema Equatorial Terrestre. Fonte: (SILVA, 2016)

Capítulo 2. Fundamentação Teórica 23

Tal sistema possui:

• Origem O𝑒𝑞 : centro de massa da Terra.

• Eixo X𝑒𝑞 : aponta na direção 𝛾 equinócio vernal.

^ = 𝐼^ × 𝐽^
• Eixo Y𝑒𝑞 : no plano equatorial tal que: 𝐾

• Eixo Z𝑒𝑞 : normal ao plano equatorial, sentido Norte.

2.1.2 Sistema de Referência Orbital

Descreve o movimento no plano orbital, como representado na Fig. 5.

Figura 5 – Sistema Orbital. Fonte:(SILVA, 2016)

Este sistema possui:

• Origem O𝑜𝑝 : centro de massa da Terra.

• Eixo x𝑜𝑝 : no plano orbital, apontado na direção do nódo ascendente.

• Eixo y𝑜𝑝 : no plano orbital tal que: 𝑘^𝑜𝑝 = ^𝑖𝑜𝑝 × ^𝑗𝑜𝑝

• Eixo z𝑜𝑝 : normal ao plano orbital.

A transformação de coordenada do Sistema de Referência Equatorial Terrestre

para o Sistema de Referência Orbital é dada por meio da rotação de um ângulo Ω (ascensão
Capítulo 2. Fundamentação Teórica 24

reta do nódo ascendente) em Z e por uma rotação de ângulo 𝐼 (inclinação) em x𝑜𝑝 , dessa
forma:

⎡ ⎤ ⎡ ⎤⎡ ⎤⎡ ⎤
⎢
𝑥𝑜𝑝 ⎥ ⎢
1 0 0 ⎥⎢
cos Ω sin Ω 0 ⎥⎢
𝑋 ⎥
𝑦𝑜𝑝 = 0 cos 𝐼 sin 𝐼 − sin Ω cos Ω 0 𝑌 (2.1)
⎢ ⎥ ⎢ ⎥⎢ ⎥⎢ ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥⎢ ⎥
⎣ ⎦ ⎣ ⎦⎣ ⎦⎣ ⎦
𝑧𝑜𝑝 0 − sin 𝐼 cos 𝐼 0 0 1 𝑍

2.1.3 Sistema de Referência Orbital Local

Coordenada usada para descrever o movimento com relação a posição do satélite,
como representado na Fig. 6

Figura 6 – Sistema Orbital Local. Fonte:(SILVA, 2016)

A Fig. 7 a seguir, apresenta o Sistema Orbital Local no plano orbital para uma
órbita elíptica.
 Capítulo 2. Fundamentação Teórica 25

Figura 7 – Sistema Orbital Local visualizado no plano orbital. Fonte:(SILVA, 2016)

Este sistema possui:

• Origem O: centro de massa do corpo.

• Eixo x𝑜 : aponta na direção do vetor velocidade orbital.

• Eixo y𝑜 : normal a órbita, oposta em sentido ao vetor momento angular orbital.

• Eixo z𝑜 : radial, no sentido Nadir.

A mudança de coordenada do Sistema de Referência Orbital para o Sistema de

Referência Orbital Local é dada por uma rotação de um ângulo 𝛽 (ângulo medido a partir
do nódo ascendente)1 mais 𝜋/2 em z𝑜𝑝 e por uma rotação de um ângulo 𝜋/2 em x dessa
forma:

𝜋 𝜋
⎡ (︂ )︂ (︂ )︂ ⎤
⎡ ⎤ ⎡ ⎤⎢ cos 𝛽 + sin 𝛽 + 0 ⎥⎡ ⎤
𝑥𝑜 1 0 0 ⎢ 2 2 ⎥ 𝑥
⎥ ⎢ 𝑜𝑝
⎥⎢ 𝜋 𝜋
⎢ ⎥ ⎢ (︂ )︂ (︂ )︂ ⎥
⎥⎢ ⎥⎢
𝑦𝑜 =⎢ 0 0 1 − sin 𝛽 + cos 𝛽 + 0 ⎥ 𝑦𝑜𝑝 (2.2)
⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎥⎢ ⎥⎢ ⎥
⎣ ⎦ ⎣ ⎦⎢
⎢ 2 2 ⎥⎣ ⎦
𝑧𝑜 0 −1 0 ⎢ ⎥ 𝑧
𝑜𝑝
0 0 1
⎣ ⎦

1
Para órbitas elípticas 𝛽 = 𝑤 + 𝑓 , em que 𝑤 representa o argumento do perigeu e 𝑓 a anomalia verdadeira.
Para órbitas circulares 𝛽 = 𝑤* 𝑡, em que 𝑤* é a velocidade angular orbital do corpo
Capítulo 2. Fundamentação Teórica 26

2.1.4 Sistema de Referência de Atitude

Usado para descrever todo o movimento de rotação do corpo ao redor de seu centro
de massa, representado na Fig. 8. Para estudos de satélites estabilizados em 3 eixos geo-
apontado, é prático definir os eixos em roll(𝜑), pitch(𝜃) e yaw(𝜓) (SILVA, 2016).
Tal sistema possui:

• Origem O: centro de massa do corpo

• Eixo de roll em x, nominalmente alinhado com x𝑜 e define o movimento em torno

da direção da velocidade orbital.

• Eixo de pitch em y, nominalmente alinhado com y𝑜 e define o movimento em torno

da direção normal à órbita, com sentido oposto ao vetor momento angular orbital.

• Eixo de yaw em z, nominalmente alinhado com z𝑜 e define o movimento em torno

da direção radial, sentido Nadir.

Figura 8 – Representação do Sistema de Atitude em verde. Fonte:(SILVA, 2016)

2.2 Atitude de Satélite

A direção para qual o satélite aponta, ou seja, a sua orientação é denominada de
atitude. Entender como se dá o movimento de atitude é importante para prever o modo
com que o satélite irá se comportar.
Existe mais de uma parametrização possível para a atitude de um satélite, dentre
elas: Cossenos diretores, Eixos de Euler, Quatérnions e Ângulos de Euler. Neste trabalho
estão apresentadas as parametrizações por quarténions, cuja aplicação mais comum é na
navegação inercial, e por ângulos de Euler. Esta última parametrização é utilizada para
controle de atitude com pequenos ângulos, pois apresenta problema de singularidade.
(WERTZ, 2002).
 Capítulo 2. Fundamentação Teórica 27

2.2.1 Representação da Atitude por Ângulos de Euler

Um modo de representar a atitude é por ângulos de Euler, são eles: 𝜑 (roll), 𝜃
(pitch) e 𝜓 (yaw). Neste trabalho os ângulos de Euler serão utilizados para mostrar a
relação entre dois sistemas de coordenadas, um deles fixo no satélite e outro associado a
um sistema inercial.
Dentre todas as possíveis combinações de rotação existentes para os ângulos de
Euler, a que será utilizada neste trabalho será a sequência de rotação 3-2-1. Esta está
representada na fig. 9

Figura 9 – Representação da sequência de rotação adotada (𝜑,𝜃,𝜓). Fonte: Modificada de

(SILVA, 2016)

A primeira rotação é do ângulo yaw em torno do eixo z0 , a segunda é do ângulo

pitch em torno do eixo y’ e a terceira é do ângulo roll em torno do eixo x. Deste modo,
pode-se calcular a matriz de rotação (SIDI, 2006) por meio de:

⎡ ⎤ ⎡ ⎤⎡ ⎤⎡ ⎤⎡ ⎤
⎢
𝑥 ⎥ ⎢
1 0 0 ⎥⎢
cos 𝜃 0 − sin 𝜃 ⎥⎢
cos 𝜓 sin 𝜓 0 ⎥⎢
𝑥𝑜 ⎥
𝑦 =⎢ 0 cos 𝜑 sin 𝜑 0 1 0 − sin 𝜓 cos 𝜓 0 𝑦𝑜 (2.3)
⎢ ⎥ ⎢ ⎥⎢ ⎥⎢ ⎥⎢ ⎥
⎢ ⎥ ⎥⎢ ⎥⎢ ⎥⎢ ⎥
⎣ ⎦ ⎣ ⎦⎣ ⎦⎣ ⎦⎣ ⎦
𝑧 0 − sin 𝜑 cos 𝜑 sin 𝜃 0 cos 𝜃 0 0 1 𝑧𝑜

Resolvendo a Equação (2.3), tem-se que a matriz de rotação A321 2 é:

⎡ ⎤
⎢
𝑐𝜃𝑐𝜓 𝑐𝜃𝑠𝜓 −𝑠𝜃 ⎥
A321 = ⎢
⎢ 𝑠𝜑𝑠𝜃𝑐𝜓 − 𝑠𝜓𝑐𝜑 𝑠𝜑𝑠𝜃𝑠𝜓 + 𝑐𝜑𝑐𝜓 𝑠𝜑𝑐𝜃 ⎥
⎥
(2.4)
⎣ ⎦
𝑐𝜑𝑠𝜃𝑐𝜓 + 𝑠𝜑𝑠𝜓 𝑐𝜑𝑠𝜃𝑠𝜓 − 𝑠𝜑𝑐𝜓 𝑐𝜑𝑐𝜃
2
c e s significam cos e sin, respectivamente
Capítulo 2. Fundamentação Teórica 28

2.2.2 Dinâmica do Satélite

A equação da dinâmica de atitude será obtida a partir da equação do momento
de Euler Equação 2.5, onde o h é o momento angular em torno do centro de massa e o
subscrito ”𝐼” indica uma derivada no referencial inercial, enquanto o subscrito ”𝐵” uma
derivada no referencial do corpo rotativo (SIDI, 2006). Será considerado a presença de
elementos rotativos dentro do satélite conhecidos como dispositivos de troca de momento.
O dispositivo utilizado neste caso será a roda de reação.

ḣI = ḣB + 𝜔 × h (2.5)

O momento angular de um satélite que é atuado por rodas de reação, segundo

(SOUZA; ARENA, 2014) é dado pela Equação 2.6

h = 𝐼𝜔 + 𝐼 w (Ω + 𝜔) (2.6)

Sendo que I=diag(I11 ,I22 ,I33 ) é a matriz do momento de inércia do satélite, 𝜔 é

a velocidade angular do satélite, Iw =diag(Iw1 ,Iw2 ,Iw3 ) é a matriz do momento de inércia
das rodas de reação e Ω = (Ω1 , Ω2 , Ω3 ) são as velocidades angulares das rodas de reação.
Diferenciando a Equação 2.6, considerando que a velocidade angular do sistema
de referencia do corpo é 𝜔 e que o torque externo é igual a N, tem-se que

ḣ + 𝜔 x h = N (2.7)

Onde a matriz 𝜔 x é chamada de matriz antissimétrica, que representa o produto

vetorial 𝜔×h igual a:

⎡ ⎤
⎢
0 −𝜔 z 𝜔 y ⎥
𝜔x = ⎢ 𝜔z 0 −𝜔 x ⎥ (2.8)
⎢ ⎥
⎣ ⎦
−𝜔 y 𝜔 x 0

Substituindo a Equação 2.6 na 2.7 e rearranjando, tem-se a Equação 2.9 que é a

dinâmica de atitude.

[︁ ]︁
𝜔˙ = (I+Iw )-1 −𝜔 x (I+Iw )𝜔 − 𝜔 x Iw Ω − Iw Ω̇ + N (2.9)

2.2.3 Cinemática do Satélite

A cinemática estuda o movimento, o qual independe das forças que o provocam.
A equação cinemática de um satélite é dada pela relação entre a velocidade angular e a
Capítulo 2. Fundamentação Teórica 29

taxa de variação da atitude do corpo. Essa equação com ângulos de Euler é apresentada
por (SIDI, 2006) como

⎡ ⎤ ⎡ ⎤⎡ ⎤
⎢
𝜑˙ ⎥ ⎢
1 𝑠𝑒𝑛𝜑𝑡𝑎𝑛𝜃 𝑐𝑜𝑠𝜑𝑡𝑎𝑛𝜃 ⎥⎢
𝜔x ⎥
𝜃˙ = 0 𝑐𝑜𝑠𝜑 −𝑠𝑒𝑛𝜑 𝜔y (2.10)
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎣ ⎦ ⎣ ⎦⎣ ⎦
𝜓˙ 0 𝑠𝑒𝑛𝜑𝑠𝑒𝑐𝜃 𝑐𝑜𝑠𝜑𝑠𝑒𝑐𝜃 𝜔z

onde 𝜔𝑥 , 𝜔𝑦 e 𝜔𝑧 são as componentes do vetor velocidade angular do satélite em

roll, pitch e yaw. Sua obtenção está explicitada no Apêndice A.
A equação da cinemática do simulador de atitude que descreve a orientação do sis-
tema de referência do corpo para o de atitude, utilizando os ângulos de Euler na sequência
3-2-1 é obtida rotacionando o eixo 3 em um ângulo 𝜓, o eixo 2 em um ângulo 𝜃 e o eixo
1 em um ângulo 𝜑. A sequência foi escolhida devido a sua singularidade em 𝜃 = 90°. Isso
é fisicamente impossível no simulador, logo a sequência pode ser estudada sem restrições
(GONZALES, 2009).

2.3 Regulador Linear Quadrático (LQR)

O LQR é uma técnica de controle baseada em realimentação de estados que consiste
em calcular um sinal de controle que, quando aplicado em uma planta em que os estados
são diferentes de zero e fazer com que os estados convirjam para zero ou um valor próximo
de zero. Este problema geralmente ocorre quando a planta está sujeita a perturbações
(PINHEIRO, 2013).
No LQR será assumido que todos os estados estão disponíveis para realimentação.
A dinâmica do sistema é dada por

˙
𝑥(𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡) (2.11)

𝑦(𝑡) = 𝐶𝑥(𝑡) + 𝐷𝑢(𝑡) (2.12)

Onde a matriz A é n×n, B uma matriz n×m, C uma matriz k×n e D uma matriz
k×m. O x é o vetor de estados, neste caso ele será composto pelos ângulos de Euler e
velocidades angulares. O u é o sinal de controle que será computado pelo LQR, aqui será
gerado pela roda de reação. E por fim y é a saída do sistema, ou seja, os valores que se
deseja controlar.
No caso estudado neste trabalho será considerado um tempo de otimização infinito,
ou seja, o controlador estará sempre em funcionamento a fim de manter o sistema o mais
Capítulo 2. Fundamentação Teórica 30

próximo possível do equilíbrio. Devido a isso será utilizado o funcional da Equação 2.13.

∫︁ ∞
𝐽= (𝑥′ 𝑄𝑥 + 𝑢′ 𝑅𝑢)𝑑𝑡 (2.13)
𝑡0

O funcional J representa a soma da energia dos estados e do sinal de controle que

deve ser minimizado em relação ao sinal de controle u. A matriz de peso Q é uma matriz
semi-definida positiva que penaliza os estados e a matriz de peso, R é uma matriz definida
positiva e penaliza os sinais de controle. Se a matriz Q aumenta enquanto R permanece
constante, o tempo de acomodação será reduzido assim como os estados os aproximarão
de zero com uma razão maior. Se R é muito maior se comparado com Q, a energia de
controle é penalizada de maneira mais intensa (PINHEIRO, 2013) e, assim, o sinal de
controle obtido possui menor amplitude. Isso resulta em menor consumo de energia pelo
sistema de controle, mas um maior tempo de acomodação.
Não há uma formulação geral que pode ser usada para a seleção das matrizes Q
e R. Uma regra simples é escolher Q e R diagonais, sempre positivos ou zero. Contudo,
é conveniente calcular diferentes controladores, baseando-se em distintos valores para as
matrizes Q e R e assim verificar sua eficiência por um processo iterativo, modificando Q
e R de acordo com resultados obtidos a cada simulação.
Será utilizado o princípio de Pontriagin (PEÑA, 2012; PINHEIRO, 2013) para a
solução da otimização do funcional que caracteriza o controlador LQR. Primeiramente é
preciso definir a Hamiltoniana

1
𝐻(𝑥, 𝜆, 𝑡) = (𝑥′ 𝑄𝑥 + 𝑢′ 𝑅𝑢) + 𝜆′ (𝐴𝑥 + 𝐵𝑢) (2.14)
2
O princípio mínimo estabelece que o controle ótimo e as trajetórias dos estados
precisam satisfazer as condições listadas abaixo

𝑥˙ = 𝜕𝐻
𝜕𝜆
−𝜆˙ = 𝜕𝐻
𝜕𝑥
(2.15)
𝜕𝐻
𝜕𝑢
=0

Aplicando essa condições na Equação 2.14, obtém-se

𝑥˙ = 𝐴𝑥 + 𝐵𝑥
−𝜆˙ = 𝑄𝑥 + 𝐴′ 𝜆 (2.16)
𝑢* = −𝑅-1 𝐵 ′ 𝜆

Em que u* é o sinal de controle ótimo, já a matriz R deverá ser definida positiva

para que sua inversa exista. As equações diferenciais presentes na Equação 2.16 constituem
um problema de valor de contorno em dois pontos que por causa da condição de contorno
Capítulo 2. Fundamentação Teórica 31

mista, torna-se difícil encontrar uma solução numérica. Escrevendo a Equação 2.15 no
formato de espaço de estados tem-se

⎡ ⎤ ⎡ ⎤⎡ ⎤
-1 ′
𝑥˙ ⎦ ⎣ 𝐴 −𝐵𝑅 𝐵 ⎦ ⎣ 𝑥 ⎦
⎣ = (2.17)
𝜆˙ −𝑄 −𝐴′ 𝜆

A matriz a direita é a hamiltoniana. O problema de valor de contorno não precisa

ser resolvido, necessariamente. Para isso a seguinte substituição será feita

𝜆 = 𝑃𝑥 (2.18)

Derivando a Equação 2.18 e substituindo na Equação 2.17 tem-se que

𝑑𝜆 𝑑𝑃
= 𝑥 + 𝑃 𝐴𝑥 − 𝑃 𝐵𝑅-1 𝐵 ′ 𝑃 𝑥 (2.19)
𝑑𝑡 𝑑𝑡
e inserindo 𝜆˙ da Equação 2.15 obtem-se

𝑑𝑃
−𝑄𝑥 − 𝐴′ 𝑃 𝑥 = 𝑥 + 𝑃 𝐴𝑥 − 𝑃 𝐵𝑅-1 𝐵 ′ 𝑃 𝑥 (2.20)
𝑑𝑡
A igualdade apresentada na equação acima deverá ser mantida para qualquer valor
de x, portanto uma condição suficiente para o controle ótimo é que P deve satisfazer

𝑑𝑃
= 𝑄 + 𝐴′ 𝑃 + 𝑃 𝐴 − 𝑃 𝐵𝑅-1 𝐵 ′ 𝑃 (2.21)
𝑑𝑡
A Equação 2.21 é conhecida como a equação diferencial de Riccati. Se o par (A,B)
é controlável, o par (A,C) é observável e as matrizes A, B, R e Q são constantes, P(t)→P
(uma matriz constante) quando tf → ∞. Se estas condições são satisfeitas, a lei de controle
ótimo para um processo de duração infinita é estacionária. Dessa forma a equação dife-
rencial de Riccati (Equação 2.21) passa a ser uma equação algébrica (PINHEIRO, 2013)
representada por

𝑃 𝐴 + 𝐴′ 𝑃 + 𝐶 ′ 𝑄𝐶 − 𝑃 𝐵𝑅-1 𝐵 ′ 𝑃 = 0 (2.22)

uma vez que 𝑑𝑃𝑑𝑡

= 0 na Equação 2.21. Com isso, o sinal de controle que minimiza
o funcional J é representado pela equação a seguir

𝑢(𝑡) = −𝐾𝑥(𝑡) (2.23)

onde K=R-1 B’P.

Capítulo 2. Fundamentação Teórica 32

Fazendo a substituição da Equação 2.23 na Equação 2.11 tem-se o controle da

dinâmica em malha fechada

˙
𝑥(𝑡) = (𝐴 − 𝐵𝐾)𝑥(𝑡) (2.24)

O sistema sendo controlável o, sistema em malha fechada será assintoticamente

estável, isso quer dizer que os autovalores de A-BK possuem parte real negativa.

2.3.1 Definição das Matrizes em Espaço de Estados

Definidos os estados que deseja-se controlar, considerando o sistema da Equação
2.10, o vetor de estados será representado por

𝑥 = [𝜑 𝜃 𝜓 𝜔𝑥 𝜔𝑦 𝜔𝑧 ]′ (2.25)

Em que os ângulos de Euler estão nas três primeiras posições, enquanto as veloci-
dades angulares do satélite estão nas outras três restantes. Agora podemos definir o vetor
de controle (u), o qual será formado pelas acelerações angulares das rodas de reação dado
a seguir

𝑢 = [Ω̇𝑥 Ω̇𝑦 Ω̇𝑧 ]′ (2.26)

Escrevendo o sistema em forma de espaço de estados, tem-se

˙
𝑥(𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡)
(2.27)
𝑦 = 𝐶𝑥 + 𝐷𝑢

Em que:

• A é a dinâmica do sistema.

• B é o matriz de controle que afeta os estados.

• C são os estados que estão sendo realimentados.

• D são os matriz de controle que afetam diretamente a saída.

Capítulo 2. Fundamentação Teórica 33

As matrizes A e B serão resultantes das expressões matriciais descritas nas Equa-

ções 2.9 e 2.10. Considerando os resultados obtidos na Equações 2.25 e 2.26. E igualando
N a zero na Equação 2.9 (PINHEIRO, 2013), tem-se que as matrizes A e B são

⎡ ⎤
0 0 0 1 𝑠𝑒𝑛𝜑𝑡𝑎𝑛𝜃 𝑐𝑜𝑠𝜑𝑡𝑎𝑛𝜃
⎢ ⎥
⎢
⎢ 0 0 0 0 𝑐𝑜𝑠𝜑 −𝑠𝑒𝑛𝜑 ⎥
⎥
⎢ ⎥
⎢
⎢ 0 0 0 0 𝑠𝑒𝑛𝜑𝑠𝑒𝑐𝜃 𝑐𝑜𝑠𝜑𝑠𝑒𝑐𝜃 ⎥
⎥
𝐴= ⎢ 𝐼 22 𝜔𝑧 −𝐼𝑤 Ω𝑧 −𝐼 33 𝜔𝑦 +𝐼𝑤 Ω𝑦
⎥ (2.28)
⎢
⎢ 0 0 0 0 𝐼 11 +𝐼𝑤 𝐼 11 +𝐼𝑤
⎥
⎥
⎢ ⎥
−𝐼 11 𝜔𝑧 +𝐼𝑤 Ω𝑧 𝐼 33 𝜔𝑥 −𝐼𝑤 Ω𝑥
0 0 0 0
⎢ ⎥
𝐼 22 +𝐼𝑤 𝐼 22 +𝐼𝑤
⎢ ⎥
⎣ ⎦
𝐼 11 𝜔𝑦 −𝐼𝑤 Ω𝑦 −𝐼 22 𝜔𝑥 +𝐼𝑤 Ω𝑥
0 0 0 𝐼 33 +𝐼𝑤 𝐼 33 +𝐼𝑤
0
⎡ ⎤
0 0 0
⎢ ⎥
⎢
⎢ 0 0 0 ⎥
⎥
⎢ ⎥
⎢
⎢ 0 0 0 ⎥
⎥
𝐵= ⎢
⎢ −𝐼𝑤
⎥ (2.29)
⎢ 𝐼 +𝐼𝑤 0 0 ⎥
⎥
⎢ 11 ⎥
−𝐼𝑤
0 0
⎢ ⎥
𝐼 22 +𝐼𝑤
⎢ ⎥
⎣ ⎦
−𝐼𝑤
0 0 𝐼 33 +𝐼𝑤

Assumindo-se que todos os estados são observáveis, então a matriz C será igual à
matiz identidade (C = I6×6 ) e a matriz D, visto que não há alimentação direta de entradas,
será igual a zero (D = 06×3 ). Assim o sistema presente na Equação 2.27 está determinado.
Para um melhor entendimento da obtenção das matrizes A e B, ver Apêndice B.
 34

3 Projeto do Mockup 1U

Neste capítulo serão detalhados todos os passos do projeto, construção e montagem

do mockup. Será também falado a respeito do padrão CubeSat.

3.1 Padrão CubeSat

O Projeto CubeSat teve início em 1999, começou como um esforço colaborativo en-
tre o Prof. Jordi Puig-Suari na Universidade Estadual Politécnica da Califórnia, San Luis
Obispo, e o Prof. Bob Twiggs no Laboratório de Desenvolvimento de Sistemas Espaciais
da Universidade de Stanford (SSDL). O objetivo do projeto é fornecer um padrão para o
projeto de nanossatélites para reduzir o custo e o tempo de desenvolvimento, aumentar a
acessibilidade ao espaço e sustentar lançamentos frequentes. Atualmente, o Projeto Cu-
beSat é uma colaboração internacional de mais de 100 universidades, escolas secundárias
e empresas privadas que desenvolvem nanossatélites contendo cargas úteis científicas, pri-
vadas e governamentais. Um CubeSat é um cubo de 10 cm com uma massa de até 1,33
kg (SLO, 2014).

3.2 Modelagem 3D
A modelagem do mockup, Fig. 10, seguiu as especificações que está presente no
Anexo A, que são fornecidas por (SLO, 2014). Porém a única exceção à regra é o topo,
como as simulação foram feitas inicialmente com ele pendurado houve a necessidade de
fazer um ponto de ancoragem.
Capítulo 3. Projeto do Mockup 1U 35

Figura 10 – CAD Mockup 1U - Vista Externa. Fonte: Autora

O software utilizado foi o CATIA V5R19 e por meio dele também foram estimados
os momentos de inércia. Na Fig. 11 é possível observar a organização interna. A roda de
reação fica na parte inferior posicionada no eixo Z, no meio fica a bateria e em cima fica
a BeagleBone Blue.

Figura 11 – CAD Mockup 1U - Estrutura interna. Fonte: Autora

Capítulo 3. Projeto do Mockup 1U 36

Os valores identificados, pelo software, dos momentos de inércia estão presentes

na Fig.12.

Figura 12 – Quadro com os dados de estimação dos momentos de inércia. Fonte: Autora

3.3 Mockup 1U - LODESTAR

A fabricação estrutural se deu por meio de manufatura aditiva e o material uti-
lizado foi o Poliácido Láctico (PLA), este não é um material qualificado para voo. Mas
como o mockup só será utilizado para fins experimentais dentro do laboratório, este não
é um problema. Este material foi escolhido por ser mais barato e por haver uma maior
facilidade na fabricação. A versão finalizada está apresentada na Fig. 13

Figura 13 – Mockup finalizado. Fonte: Autora

Capítulo 3. Projeto do Mockup 1U 37

A parte interna está ilustrada na Fig. 14, a configuração dos componentes manteve-
se a mesma do projeto inicial.

Figura 14 – Interior do Mockup. Fonte: Autora

3.3.1 Roda de Reação

A roda de reação está posicionada de modo que coincide com o eixo Z do mockup.
Ela utiliza motor brushless DC da Faulhaber 2610006BSC e fica acoplada no eixo de
rotação, na Tabela 1 é possível ver as informações técnicas da roda. Esse dispositivo
conta com um driver interno que permite controlar a direção e velocidade de rotação do
motor. O sistema é energizado por uma bateria Lipo Feasso FFB-006 (15C) 1200mAh
7.4v.
O motor possui seis pinos, como mostra o esquemático da Fig. 15. O pino 1 é a
alimentação do driver interno do motor, o pino 2 é a entrada de tensão para as bobinas e o
pino 3 a referência. Os pinos 4, 5 e 6 são pinos digitais sendo o 4 o controle da velocidade
de rotação, o pino 5 a direção de rotação e o pino 6 uma saída que permite verificar a
velocidade angular do motor através do período das ondas quadradas fornecidas (SILVA,
2013).
Capítulo 3. Projeto do Mockup 1U 38

Figura 15 – Esquemático de ligação do motor. Fonte: (SILVA, 2013)

Tabela 1 – Características técnicas da Roda de Reação

Parâmetro Valor
Torque máximo [𝑚𝑁 × 𝑚] 3,7
Velocidade máxima [𝑅𝑃 𝑀 ] 7000
Inércia da Roda de Reação [𝐾𝑔𝑚 ] 32×10-6
2

Momento [𝑚𝑁 × 𝑚 × 𝑠] 23,4

3.3.2 BeagleBone Blue

O 𝐵𝑒𝑎𝑔𝑙𝑒𝐵𝑜𝑛𝑒r Blue é um computador tudo-em-um baseado em Linux para ro-
bótica, integrando em uma única placa pequena o microprocessador Octavo OSD3358
junto com wifi, bluetooth, IMU, barômetro, regulação de potência e estado de LEDs de
carga para um LiPo de 2 células, ponte-H e conectores discretos para 4 motores DC +
codificadores, 8 servos e todos os barramentos comumente necessários para periféricos
adicionais em aplicações embarcadas (BEAGLEBOARD.ORG, ).
Durante o processo de escolha do computador de bordo, haviam 3 opções: o Ar-
duino, o Raspbarry Pi 3 e o BeagleBone Blue. Porém, quando levou-se em consideração
toda a eletrônica que seria necessária embarcar, tal como, IMU 9-eixos, um carregador de
bateria 2-celulas, transmissão de dados, entre outros. O BeagleBone foi o que proporcionou
maiores vantagens, dentre as opções citadas.
 39

4 Análise Numérica

Neste capítulo serão apresentadas as simulações feitas, bem como as matrizes para
o problema de 1-eixo e os resultados obtidos das simulações.

4.1 Análise em 3-eixos

Os espaços de estados utilizados foram calculados no capítulo 2 e os valores dos
estados estão presentes na Tabela 2. Os momentos de inércia foram retirados do quadro
da Fig. 12

Tabela 2 – Parâmetros da simulação

Parâmetro Valor
𝐼 xx = 7.328 × 10−4
Momento de Inércia do simulador [𝐾𝑔𝑚2 ] 𝐼 yy = 7.294 × 10−4
𝐼 zz = 6.332 × 10−4
Velocidade Angular Inicial do Simulador [𝑔𝑟𝑎𝑢𝑠/𝑠] 𝜔𝑥 = 𝜔𝑦 = 𝜔𝑧 = 0
2
Momento de Inércia da Roda de Reação [𝐾𝑔𝑚 ] Iw = 32×10-6
Velocidade Angular Inicial das rodas de reação [𝑔𝑟𝑎𝑢𝑠/𝑠] Ω𝑥 = Ω𝑦 = Ω𝑧 = 0
Ângulos de Euler (1º Caso) [𝑔𝑟𝑎𝑢𝑠] 𝜑 = 1°; 𝜃 = 1°; 𝜓 = 15°
Ângulos de Euler (2º Caso) [𝑔𝑟𝑎𝑢𝑠] 𝜑 = 1°; 𝜃 = 1°; 𝜓 = 90°

Para a seleção dos pesos das matrizes Q e R foram feitos testes com vários valo-
res até que um resultado satisfatório fosse obtido, foram escolhidos dois valores de peso
diferentes para a matriz R. Os resultados escolhidos estão apresentados na Tabela 3.
Como Q e R são matrizes diagonais, tem-se que Q = diag[𝑄1 𝑄2 𝑄3 𝑄4 𝑄5 𝑄6 ] e R =
diag[𝑅1 𝑅2 𝑅3]

Tabela 3 – Pesos de Q e R
1º Caso 2º Caso
Matriz Valor Matriz Valor
Q1...3 0,0001 Q1...3 2000
Q4...6 1 Q4...6 1
R1...3 1100 R1...3 0,46

Nas Fig. 16, 17 e 18 estão os gráficos resultantes das simulações, os que estão locali-
zados a esquerda são referentes aos dados do 1º caso e os que estão a direita são referentes
aos dados do 2º caso. O código utilizado para fazer estas simulações está disponível no
Apêndice D.
Capítulo 4. Análise Numérica 40

Figura 16 – Ângulos Yaw, Pitch e Roll sendo controlados. Fonte:Autor

Nas figuras são observadas duas situações diferentes, a primeira referente ao 1º caso
foi priorizado encontrar uma matrizes de ganho que fizessem com que o torque máximo
fosse bem próximo ao das rodas de reação. Porém, para conseguir isso o tempo de controle
teve um valor considerável (aproximadamente 48 minutos até que todo o sistema fosse
estabilizado). Esse resultado foi possível para ângulos pequena, isso justifica a escolha dos
15° para esta primeira análise.

Figura 17 – Velocidade angular. Fonte:Autor

O 2º caso, tinha como objetivo era obter um sistema que reagisse com uma maior
rapidez. Com isso, o torque aumentou. O ângulo escolhido foi de 90° e a análise dos
dois casos está sendo feita com variação apenas no eixo Z. O tempo conseguido foi 6,2
segundos, obtendo uma resposta bem rápida. Contudo ao avaliar o gráfico de torque da
Capítulo 4. Análise Numérica 41

Fig. 18 o torque vai além da capacidade do motor da roda indicando assim uma possível
saturação.
Como o foco maior deste trabalho é controlar o sistema com rapidez o 2º caso
será estudado com maior profundidade. Para isso, o sistema será simplificado apenas para
análise em 1 eixo (O eixo Z permanecerá) e o torque será então reavaliado.

Figura 18 – Torque das rodas de reação. Fonte:Autor

4.2 Análise em 1-eixo

A análise será feita em torno do eixo Z, a partir das equações encontradas no
Capítulo 2. Com isso vamos considerar que 𝜔𝑥 = 𝜔𝑦 = Ω𝑥 = Ω𝑦 = 0, 𝜃 = 𝜑 = 0. Partindo
disso, temos que as matrizes dos espaço de estado para 1 eixo fica:

⎡ ⎤
𝜓 ⎦
𝑥=⎣ (4.1)
𝜔𝑧

[︁ ]︁
𝑢= Ω̇𝑧 (4.2)

⎡ ⎤
0 1 ⎦
𝐴=⎣ (4.3)
0 0
⎡ ⎤
0
𝐵=⎣ ⎦ (4.4)
−𝐼𝑊 /(𝐼33 + 𝐼𝑊 )
Capítulo 4. Análise Numérica 42

⎡ ⎤
1 0 ⎦
𝐶=⎣ (4.5)
0 1

[︁ ]︁
𝐷= 0 (4.6)

Utilizando os valores de 𝐼𝑧𝑧 e 𝐼𝑊 presentes na Fig 12 e Tabela 1, respectivamente.

Com o auxílio do código presente no Apêndice D, adaptando para este caso, foi analisado
a resposta ao degrau.
As matrizes de ganho Q e R selecionadas para este caso que resultaram em um
tempo de pico igual a 2,45 segundos e tempo de acomodação 4,5 segundos foram

⎡ ⎤
2000 0 ⎦
𝑄=⎣ 𝑅 = [0, 46] (4.7)
0 1

Figura 19 – Gráfico da resposta ao degrau (Ângulo). Fonte: Autora

Figura 20 – Gráfico da resposta ao degrau (Velocidade Ângular). Fonte: Autora

Ao analisar o próximo gráfico (Fig. 21), percebe-se que o torque máximo necessário
para o controle ultrapassa aquele que está na Tabela 1, que é o máximo da roda. Isso quer
dizer que a roda irá saturar antes que o controle seja feito, assim como na análise feita
Capítulo 4. Análise Numérica 43

anteriormente. Contudo o gráfico da Fig. 19 mostra que no momento de pico a roda não
atingiu sua velocidade máxima, que é 7000 RPM.

Figura 21 – Gráfico do Torque necessário. Fonte: Autor

Uma análise experimental será feita para a verificação de como o sistema se com-
portará.
 44

5 Análise Experimental

Neste capítulo será descrito como foi feita e ocorreu a parte de simulação experi-
mental em um eixo, serão também apresentados os resultados.

5.1 Modelando o Sistema

A planta do sistema foi desenhada, levando-se em consideração o sistema real.
Logo terá como entrada uma tensão e não o torque como antes. A planta do sistema está
representada na Fig. 22

Figura 22 – Planta do sistema real. Fonte: Autor

Utilizando a Regra de Mason, a partir da planta é possível obter a Equação 5.1

𝐾 𝑏
𝜃 𝐾 𝑎 𝐽𝑀 𝐾
= 𝐾𝑏 𝐾𝑣 ⇔ (5.1)
𝐸𝑎 𝑠(𝑠 + 𝐾𝑎 𝐽𝑤 ) 𝑠(𝑠 + 𝑎)

Onde 𝐽𝑀 é o momento de inércia do mockup e o 𝐽𝑤 é o momento de inércia da

roda de reação.
Desta equação, pode-se verificar que as matrizes A e B irão ser modificadas em
relação às apresentadas no capítulo anterior.

⎡ ⎤
𝜓 ⎦
𝑥=⎣ (5.2)
𝜔𝑧
⎡ ⎤
0 1 ⎦
𝐴=⎣ (5.3)
0 −𝑎
⎡ ⎤
0 ⎦
𝐵=⎣ (5.4)
𝐾
Capítulo 5. Análise Experimental 45

Para descobrir os valores de a e K, foi realizado um experimento. A malha foi

fechada e aplicou-se uma onda quadrada, no sistema, que ia de 0° até 90°. Para isso, o
mockup foi pendurado como mostrado na Fig. 23. O objetivo era verificar se é possível
controlar o sistema em 10 segundos, sem saturar as rodas. Este limite de tempo e escolha
do ângulo, foram inspirados no 2° CubeDesign promovido pelo Instituto Nacional de
Pesquisas Espaciais (INPE).

Figura 23 – Montagem do experimento. Fonte: Autor

O sinal foi enviado do computador para o Beaglebone Blue, e a coleta dos resultados
foram feitas. Um ajuste de ganho foi feito até se obter um resultado satisfatório que foi o
mostrado na Fig. 24

Figura 24 – Resultado da simulação. Fonte: Autor

No gráfico apresentado na Fig. 24 , os 30 primeiros segundos somente um PI

está atuando sobre o sistema para levar o mockup para a referência 0°. Somente após
Capítulo 5. Análise Experimental 46

ter chegado muito próximo à referência o controlador em espaço de estados começa a

funcionar, isso é para que nenhuma perturbação atrapalhe o experimento. A perturbação
inicial é gerada pelo próprio motor, pelo fato do driver deste ser alimentado antes do
controlador. Com isso o LQR só começou a atuar após 30 segundos, onde iniciou-se o
recebimento da onda quadrada.
Esse resultado foi muito satisfatório, ampliando um pouco a imagem, podemos ver
na Fig 21. O tempo de pico foi de 2,45 segundos e o de acomodação 4,5 segundos, isso
fornece um sobre sinal de 0,5 o que é adequado para o caso

Figura 25 – Dados da simulação.(Zoom) Fonte: Autor

O gráfico da Fig.26 mostra a porcentagem de tensão que está sendo utilizada pelo
motor,o valor máximo observado (durante a ação do LQR) foi de 32,01%. Isso quer dizer
que o motor está trabalhando com 32% da sua capacidade apenas. Consequentemente,
pode-se dizer que o controlador atua mais rápido do que ele já está atuando.

Figura 26 – Percentual de tensão aplicada ao motor ao utilizar sistema de controle baseado

em PI + LQR. Fonte: Autor

Os ganhos K’s utilizados neste experimento estão apresentados na Tabela 4

Capítulo 5. Análise Experimental 47

Tabela 4 – Ganhos
𝐾1 𝐾2
0.4000 0.1962

Feito isso, tendo o valor do sobressinal e do tempo de pico utilizando o software

MATLAB obteve-se os seguintes valores para as matrizes a e K.

⎡ ⎤
0 1
𝐴=⎣ ⎦ (5.5)
0 −0, 8048
⎡ ⎤
0
𝐵=⎣ ⎦ (5.6)
8, 5979
 48

6 Conclusão

Neste trabalho os testes foram concluídos de forma satisfatória. A importância da

validação de um sistema após projetado ficou claro aqui, pois sem a simulação o contro-
lador teria sido rejeitado e tido como inapropriado para o caso. A análise experimental
mostrou exatamente o comportamento do sistema, o que possibilitou o ajuste do contro-
lador que nos deu uma resposta melhor que a esperada, prevista em análise numérica.
Os dados aqui apresentados foram importantes para experimentos futuros dentro
do laboratório LODESTAR, principalmente no que diz respeito ao torque da roda de
reação. Aqui ficou evidente que o modelo não representou o torque requerido para o
controle do sistema. Na análise numérica o torque obtido foi muito além do que está
especificado no datasheet da roda, indicando assim a sua saturação antes da realização do
controle. Isso mostra que há uma necessidade de rever o equacionamento proposto para
investigar o que está causando tamanho erro.
Os testes a priore não foram feitos na plataforma simuladora uma vez que essa
possui momento de inércia muito grande para a capacidade da roda, por isso escolheu-se
fazer o teste dessa maneira.
Este trabalho foi realizado com o objetivo de realizar um estudo dentro do labo-
ratório, e os primeiros passos foi realizar com um mockup que embora esteja construído
dentro dos padrões de dimensão, os materiais utilizados não são os mesmos usados na
construção de um satélite real. Isso indica que quando os testes forem realizados para
um nanossatélite real, os momentos de inércia não serão os mesmos apresentados neste
trabalho. Necessitando assim, uma atenção na escolha dos valores dos ganhos.
 49

7 Trabalhos futuros

Encontrar o real problemas do modelo que represente o torque necessário requerido

pela roda é importante para que os estudos no possam avançar. Com o modelo definido a
simulação precisa ser feita na plataforma simuladora para análise dos resultados e então
avançar os estudos para atuação nos 3-eixos.
Uma configuração de posicionamento das rodas sugerida é posicionar cada uma
em um eixo (X,Y,Z) do mockup, como mostrado na Fig. 27.

Figura 27 – Mockup com atuadores nos 3 eixos. Fonte: Autora

Mas vale ressaltar que um estudo relacionado à configuração do posicionamento e

quantidade de rodas de reação dentro do mockup deve ser feito para verificar a eficiência.
O estudo para outros controladores também deve ser feito, para o comportamento
deles no simulador do LODESTAR contribuindo assim para o seu aperfeiçoamento.
 50

Referências

BEAGLEBOARD.ORG. BeagleBone Blue: What is beaglebone blue? [S.l.]. Disponível

em: <http://beagleboard.org/blue>. Citado na página 38.

FEHSE, W. Automated Rendezvous and Docking of Spacecraft. New York, USA:

Cambridge University Press, 2003. 495 p. Citado na página 22.

GONZALES, R. G. Utilização dos Métodos SDRE e Filtro de Kalman para o controle

de atitude de simuladores de satélite. 120 p. Dissertação (Mestrado do Curso de
Pós-Graduação em Engenharia e Tecnologia Espacial) — Instituto Nacional de Pesquisas
Espaciais, São José dos Campos, 2009. Citado na página 29.

ISHIOKA, I. S. K. Desenvolvimento de um atuador magnético para o sistema de controle

de atitude de nanosatélites. Brasília, p. 104, 2017. Citado 2 vezes nas páginas 17 e 18.

PEÑA, L. G. Dinâmica e Controle de Microsatélites em Órbitas Baixas. 81 p. Tese

(Doutorado) — Instituto Nacional de Pesquisas Espaciais, Curso de Pós-Graduação em
Engenharia e tecnologia Espaciais com concentração em Mecânica Espacial e Controle,
São José dos Campos, 2012. Citado na página 30.

PINHEIRO, E. R. Aplicação do Controlador Misto H2 /H∞ no Controle na atitude de

um microssatélite na Presença de Incertezas Paramétricas. 117 p. Dissertação (Mestrado
do Curso de Pós-Graduação em Engenharia e Tecnologia Espacial/ Mecânica Espacial
e Controle) — Instituto Nacional de Pesquisas Espaciais, São José dos Campos, 2013.
Citado 4 vezes nas páginas 29, 30, 31 e 33.

PLOEG, L. C. V. Desenvolvimento de sistema para simulaÇÃo do campo magnÉtico

terrestre em Órbitas baixas. p. 65, 2015. Trabalho de Graduação em Engenharia
de Controle e Automação, Publicação FT.TG-n010/2017, Faculdade de Tecnologia,
Universidade de Brasília, DF. Citado 2 vezes nas páginas 9 e 19.

RUSSEL, M.; STRAUB, J. Characterization of command software for an autonomous

attitude determination and control system for spacecraft. International Journal of
Computers and Applications, Taylor & Francis Group, Lynchburg, p. 13, 2017. Citado
na página 17.

SIDI, M. J. Spacecraft Dynamics and Control: A pratical engineering approach.

Cambrige: Cambridge University Press, 2006. 409 p. Citado 3 vezes nas páginas 27, 28
e 29.

SILVA, L. M. B. Projeto e implementação de um sistema de rodas de reação para uma

bancada de simulação de atitude de nanosatélites. p. 58, 2013. Trabalho de Conclusão de
Curso em Engenharia Aeroespacial, Faculdade do Gama, Universidade de Brasília, DF.
Citado 3 vezes nas páginas 9, 37 e 38.

SILVA, R. C. Batch identification of the reaction-wheel-based actuation system of the

microsatellite test bed. Brasília, p. 1–6, 2016. Citado 2 vezes nas páginas 9 e 17.
Referências 51

SILVA, R. C. da; RODRIGUES, U. A. Simulador de sistema de determinação e controle

de atitude de pequenos satélites. p. 112, 2015. Trabalho de Graduação em Engenharia
de Controle e Automação, Publicação FT.TG-n 14/2015, Faculdade de Tecnologia,
Universidade de Brasília, DF. Citado 2 vezes nas páginas 9 e 18.

SILVA, W. R. Filtro H∞ estendido de segunda ordem e filtro de partículas regularizado

com Roughening aplicados na estimação de atitude de satélites artificiais. 145 p. Tese
(Doutorado) — Instituto Nacional de Pesquisas Espaciais, Curso de Pós-Graduação
em Engenharia e tecnologias Espaciais / Mecânica Espacial e Controle, São José dos
Campos, 2016. Citado 7 vezes nas páginas 9, 22, 23, 24, 25, 26 e 27.

SLO, C. P. Cubesat design specification rev. 13. p. 42, 2014. Citado 3 vezes nas páginas
9, 34 e 61.

SOUZA, L. C. G. de; ARENA, V. M. R. Design of satellite attitude control algorithm

based on the sdre method using gas jets and reaction wheels. Journal of Engineering,
Rio de Janeiro, p. 9, 2014. Citado na página 28.

WERTZ, J. R. Spacecraft Attitude Determination and Control. Dordrecht, Boston,

London: Klumer Academic Publishers, 2002. v. 73. 858 p. Citado na página 26.

ZIPFEL, P. H. Modeling and Simulation of Aerospace Vehicle Dynamics. 2. ed. Virginia:

AIAA education series, 2007. 567 p. Nenhuma citação no texto.
Apêndices
 53

APÊNDICE A – Equação da Cinemática em

Ângulos de Euler

A velocidade angular do satélite,no sistema de referência de atitude, ao redor do

próprio corpo é dada pela Equação A.1

𝜔 = 𝜔𝑥^i + 𝜔𝑦^j + 𝜔𝑧 k
^ (A.1)

Observando a Fig. 9, a velocidade angular pode ser representada (utilizando a

sequência 3-2-1) por:

′
𝜔 = 𝜓˙ 𝑘^0 + 𝜃˙^j + 𝜑˙^i (A.2)

Por meio da matriz de rotação da Equação 2.4 é possível descobrir os valores de

k^0 e ^j’, uma vez que o valor de ^𝑖 já é conhecido. Substituindo esses valores na Equação
A.2 tem-se que

⎡ ⎤ ⎡ ⎤ ⎡ ⎤
⎢
−𝑠𝑒𝑛𝜃 ⎥ ⎢
0 ⎥ ⎢
1 ⎥
⎥ ˙ ⎥ ˙ ⎢ ⎥ ˙
𝜔=⎢
⎢ 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥𝜓 +⎢
⎢ 𝑐𝑜𝑠𝜑 ⎥𝜃 + ⎢ 0 ⎥𝜑 (A.3)
⎣ ⎦ ⎣ ⎦ ⎣ ⎦
𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 −𝑠𝑒𝑛𝜑 0

Logo,

⎡ ⎤⎡ ⎤
⎢
1 0 −𝑠𝑒𝑛𝜃 ⎥ ⎢ 𝜑˙ ⎥
𝜔=⎢ 0 𝑐𝑜𝑠𝜑 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥ 𝜃˙ (A.4)
⎢ ⎥⎢ ⎥
⎢ ⎥
⎣ ⎦⎣ ⎦
0 −𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 𝜓˙

Rearranjando a Equação A.4 obtêm-se

⎡ ⎤ ⎡ ⎤−1 ⎡ ⎤
⎢
𝜑˙ ⎥ ⎢
1 0 −𝑠𝑒𝑛𝜃 ⎥ ⎢
𝜔𝑥 ⎥
𝜃˙ =⎢ 0 𝑐𝑜𝑠𝜑 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥ 𝜔𝑦 (A.5)
⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ ⎥ ⎢ ⎥
⎣ ⎦ ⎣ ⎦ ⎣ ⎦
𝜓˙ 0 −𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 𝜔𝑧

considerando que,

⎡ ⎤
⎢
1 0 −𝑠𝑒𝑛𝜃 ⎥
Λ=⎢ 0 𝑐𝑜𝑠𝜑 𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥ (A.6)
⎢ ⎥
⎣ ⎦
0 −𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃
APÊNDICE A. Equação da Cinemática em Ângulos de Euler 54

A inversa de Λ pode ser calculada pela seguinte expressão:

(𝑐𝑜𝑓 (Λ))𝑇
Λ−1 = (A.7)
𝑑𝑒𝑡(Λ)

Desenvolvendo os cálculos, têm-se:

𝑑𝑒𝑡(Λ) = 𝑐𝑜𝑠𝜃 (A.8)

⎡ ⎤
⎢
(−1)1+1 (𝑐𝑜𝑠𝜃) (−1)1+2 (0) (−1)1+3 (0) ⎥
𝑐𝑜𝑓 (Λ) = ⎢ (−1)2+1 (−𝑠𝑒𝑛𝜑𝑠𝑒𝑛𝜃) (−1)2+2 (𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃) (−1)2+3 (−𝑠𝑒𝑛𝜑) (A.9)
⎢ ⎥
⎥
⎣ ⎦
(−1)3+1 (𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑) (−1)3+2 (𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃) (−1)3+3 (𝑐𝑜𝑠𝜑)

⎡ ⎤
⎢
𝑐𝑜𝑠𝜃 0 0 ⎥
𝑐𝑜𝑓 (Λ) = ⎢ 𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 𝑠𝑒𝑛𝜑 ⎥ (A.10)
⎢ ⎥
⎣ ⎦
𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑 −𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 𝑐𝑜𝑠𝜑

Então a transposta de A.10 será

⎡ ⎤
⎢
𝑐𝑜𝑠𝜃 𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑 𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑 ⎥
𝑐𝑜𝑓 (Λ)𝑇 = ⎢ 0 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 −𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥ (A.11)
⎢ ⎥
⎣ ⎦
0 𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑

Substituindo as Equações A.11 e A.8 na expressão A.7 é obtido o valor da matriz

inversa de Λ

⎡ ⎤
𝑐𝑜𝑠𝜃 𝑠𝑒𝑛𝜃𝑠𝑒𝑛𝜑 𝑠𝑒𝑛𝜃𝑐𝑜𝑠𝜑 ⎥
1 ⎢
Λ−1 = ⎢
⎢ 0 𝑐𝑜𝑠𝜑𝑐𝑜𝑠𝜃 −𝑠𝑒𝑛𝜑𝑐𝑜𝑠𝜃 ⎥
⎥
(A.12)
𝑐𝑜𝑠𝜃 ⎣ ⎦
0 𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑

Substituindo a Equação A.12 na A.5, têm a Equação da Cinemática do problema

⎡ ⎤ ⎡ ⎤⎡ ⎤
⎢
𝜑˙ ⎥ ⎢
1 𝑠𝑒𝑛𝜑𝑡𝑎𝑛𝜃 𝑐𝑜𝑠𝜑𝑡𝑎𝑛𝜃 ⎥⎢
𝜔x ⎥
𝜃˙ =⎢ 0 𝑐𝑜𝑠𝜑 −𝑠𝑒𝑛𝜑 𝜔y (A.13)
⎢ ⎥ ⎥⎢ ⎥
⎢ ⎥ ⎢ ⎥⎢ ⎥
⎣ ⎦ ⎣ ⎦⎣ ⎦
𝜓˙ 0 𝑠𝑒𝑛𝜑𝑠𝑒𝑐𝜃 𝑐𝑜𝑠𝜑𝑠𝑒𝑐𝜃 𝜔z
 55

APÊNDICE B – Matrizes do Espaço de

Estado

Reescrevendo a Equação matricial B.1, porém considerando os torques externos

igual a zero.

[︁ ]︁
𝜔˙ = (I+Iw )-1 −𝜔 x (I+Iw )𝜔 − 𝜔 x Iw Ω − Iw Ω̇ (B.1)

Obtêm-se,

𝜔˙ = 𝐴𝜔 𝜔 + 𝐵 Ω̇ + 𝐶𝜔 (B.2)

Onde
𝐴𝜔 = −(𝐼 + 𝐼𝑊 )−1 𝜔 𝑥 (𝐼 + 𝐼𝑊 ) (B.3)

𝐵𝜔 = −(𝐼 + 𝐼𝑊 )−1 (𝐼𝑊 ) (B.4)

𝐶𝜔 = −(𝐼 + 𝐼𝑊 )−1 𝜔 𝑥 (𝐼𝑊 )Ω (B.5)

Fazendo as devidas substituições em B.3, tem-se que

⎡ ⎤
𝐼22
0 𝜔
(𝐼11 +𝐼𝑊 ) 𝑧
− (𝐼11𝐼+𝐼
33
𝑊)
𝜔𝑦
⎢ ⎥
⎢ − (𝐼 𝐼+𝐼
𝐴1𝜔 = ⎢ 11
𝜔𝑧 0 𝐼33
𝜔 (B.6)
⎥
𝑊) (𝐼22 +𝐼𝑊 ) 𝑥
⎥
⎣ 22 ⎦
𝐼11
𝜔
(𝐼33 +𝐼𝑊 ) 𝑦
− (𝐼33𝐼+𝐼
22
𝑊)
𝜔𝑥 0

Repetindo as mesma operações para B.4, obtem-se a seguinte matriz

−𝐼𝑊
⎡ ⎤
𝐼 11 +𝐼𝑊
0 0
⎢ ⎥
−𝐼𝑊
𝐵𝜔 = ⎢
⎢ 0 𝐼 22 +𝐼𝑊
0
⎥
⎥ (B.7)
⎣ ⎦
−𝐼𝑊
0 0 𝐼 33 +𝐼𝑊

Por último, fazendo o mesmo processo com B.5

⎡ ⎤
𝐼𝑊
𝐼 11 +𝐼𝑊
(Ω𝑦 𝜔𝑧 − Ω𝑧 𝜔𝑦 )
⎢ ⎥
𝐼𝑊
𝐶𝜔 = ⎢ (Ω𝑧 𝜔𝑥 − Ω𝑥 𝜔𝑧 ) ⎥ (B.8)
⎢ ⎥
⎣ 𝐼 22 +𝐼𝑊 ⎦
𝐼𝑊
𝐼 33 +𝐼𝑊
(Ω𝑥 𝜔𝑦 − Ω𝑦 𝜔𝑥 )
APÊNDICE B. Matrizes do Espaço de Estado 56

Colocando as velocidades angulares do satélite em evidência na Equação B.8, 𝐶𝜔

fica

⎡ ⎤⎡ ⎤
𝐼𝑊 𝐼𝑊
0 (𝐼11 +𝐼𝑊 )
(−Ω𝑧 ) (𝐼11 +𝐼𝑊 )
(Ω𝑦 ) 𝜔𝑥
⎢ ⎥⎢ ⎥
𝐼𝑊 𝐼𝑊
𝐶𝜔 = ⎢ (Ω𝑧 ) 0 (−Ω𝑥 ) 𝜔𝑦 (B.9)
⎢ ⎥⎢ ⎥
(𝐼22 +𝐼𝑊 ) (𝐼22 +𝐼𝑊 )
⎥⎢ ⎥
⎣ ⎦⎣ ⎦
𝐼𝑊 𝐼𝑊
(𝐼33 +𝐼𝑊 )
(−Ω𝑦 ) (𝐼33 +𝐼𝑊 )
(Ω𝑥 ) 0 𝜔𝑧

A matriz B.9 pode então ser somada à matriz B.6, resultando assim na matriz 𝐴𝜔
final

−𝐼33 𝜔𝑦 +𝐼𝑊 Ω2
⎡ ⎤
𝐼22 𝜔𝑧 −𝐼𝑊 Ω𝑧
0 (𝐼11 +𝐼𝑊 ) (𝐼11 +𝐼𝑊 )
⎢ ⎥
⎢ −𝐼11 𝜔𝑧 +𝐼𝑊 Ω𝑧 𝐼33 𝜔𝑥 −𝐼𝑊 Ω𝑥
𝐴𝜔 = ⎢ 0 (B.10)
⎥
(𝐼22 +𝐼𝑊 ) (𝐼22 +𝐼𝑊 )
⎥
⎣ ⎦
𝐼11 𝜔𝑦 −𝐼𝑊 Ω𝑦 −𝐼22 𝜔𝑥 −𝐼𝑊 Ω𝑥
(𝐼33 +𝐼𝑊 ) (𝐼33 +𝐼𝑊 )
0

A Equação B.10 permite reescrever a Equação B.2, ficando

𝜔˙ = 𝐴𝜔 𝜔 + 𝐵 Ω̇ (B.11)

Já que agora o 𝐶𝜔 está inserido em 𝐴𝜔 . Com isso, a partir das Equações B.11 e
2.10 o sistema pode ser descrito na forma de espaço de estados.

˙
𝑥(𝑡) = 𝐴𝑥(𝑡) + 𝐵𝑢(𝑡) (B.12)
 57

APÊNDICE C – Montagem do Mockup

Figura 28 – Montagem do exterior do Mockup. Fonte: Autora

APÊNDICE C. Montagem do Mockup 58

Figura 29 – Montagem do Interior do Mockup. Fonte: Autora

 59

APÊNDICE D – Código LQR Genérico

%I n p u t s
I 1 ; I 2 ; I 3 ; Iw ; O1 ; O2 ; O3 ; t ; w1 ; w2 ; w3 ; p hi ; t h e t a ; p s i ;

%S t a t e v e c t o r
x0 = [ p h i t h e t a p s i w1 w2 w3 ] ;

%M a t r i c e s
A= [ ] ; B = [ ] ; C= [ ] ; N= [ ] ; Q [ ] ; R [ ] ;

%G e t t i n g th e matrix K
[K] = l q r (A, B,Q, R,N ) ;

%Closed−l o o p system
s y s _ c l = s s (A−B∗K, B, C , 0 ) ;

%I n i t i a l f u n c t i o n
x= i n i t i a l ( sys_cl , x0 , t ) ;

%C o n t r o l s i g n a l
u=−K∗x ’ ;
Anexos
 61

ANEXO A – Especificações do CubeSat 1U

Figura 30 – Desenho com as especificações do CubeSat 1U. Fonte:(SLO, 2014)