Pos-Graduacao em Engenharia Mecanic

See discussions, stats, and author profiles for this publication at: https://www.researchgate.
net/publication/331797090
II Simpósio de Pós-Graduação em Engenharia Mecânica da Universidade

Federal de Campina Grande: Coletânea de artigos
Book · January 2019

DOI: 10.5935/978-85-7042-048-0
CITATIONS READS
0 567
35 authors, including:
Raissa Queiroga Talles Oliveira

Universidade Federal da Paraíba Universidade Federal de Campina Grande (UFCG)
4 PUBLICATIONS 1 CITATION 17 PUBLICATIONS 10 CITATIONS
SEE PROFILE SEE PROFILE
Romulo Pierre Batista Dos Reis Antonio Almeida

Universidade Federal Rural do Semi-Árido - UFERSA Universidade Federal de Campina Grande (UFCG)
25 PUBLICATIONS 115 CITATIONS 55 PUBLICATIONS 258 CITATIONS
SEE PROFILE SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Development of computational tools for the modeling and analysis of corrosion defects in pressure vessels View project
Anaerobic sludge, Anaerobic reactor, Wastewater View project
All content following this page was uploaded by Wanderley Ferreira Amorim Jr on 21 October 2019.
The user has requested enhancement of the downloaded file.

Raíssa Alves Queiroga
Kaline Ventura Batista
(Organizadoras)
II Simpósio de Pós-Graduação em Engenharia

Mecânica da Universidade Federal de Campina
Grande: Coletânea de artigos
1ª Edição
Belo Horizonte
Poisson
2019
Editor Chefe: Dr. Darly Fernando Andrade
Conselho Editorial
Dr. Antônio Artur de Souza – Universidade Federal de Minas Gerais
Ms. Davilson Eduardo Andrade
Dr. José Eduardo Ferreira Lopes – Universidade Federal de Uberlândia
Dr. Otaviano Francisco Neves – Pontifícia Universidade Católica de Minas Gerais
Dr. Luiz Cláudio de Lima – Universidade FUMEC
Dr. Nelson Ferreira Filho – Faculdades Kennedy
Ms. Valdiney Alves de Oliveira – Universidade Federal de Uberlândia
Dados Internacionais de Catalogação na Publicação (CIP)

S612
II Simpósio de Pós-Graduação em Engenharia
Mecânica da Universidade Federal de
Campina Grande: coletânea de artigos /
Organização Raíssa Alves Queiroga, Kaline
Ventura Batista - Editora Poisson – Belo
Horizonte - MG : Poisson, 2019
130p
Formato: PDF
ISBN: 978-85-7042-048-0
DOI: 10.5935/978-85-7042-048-0
Modo de acesso: World Wide Web

Inclui bibliografia
1. Engenharia 2. Engenharia Mecânica.

I. Título
CDD-620
O conteúdo dos artigos e seus dados em sua forma, correção e confiabilidade são
de responsabilidade exclusiva dos seus respectivos autores.
www.poisson.com.br
contato@poisson.com.br
II Simpósio de Pós-Graduação em Engenharia Mecânica da
Universidade Federal de Campina Grande
Realização
Programa de Pós-Graduação em Engenharia Mecânica - UFCG/CCT
Organizadoras do livro
Comitê Organizador SiPGEM 2018

Eng. Petróleo - Raíssa Alves Queiroga – (Presidente)
Eng. Mecânica - Kaline Ventura Batista – (Vice-Presidente)
Eng. Mecânico - Marcélio Ronnie Dantas de Sá – (Secretário)
Eng. Mecânico - João Manoel de Oliveira Neto – (Secretário)
Eng. Mecânica - Aline Michelly do Nascimento – (Tesoureira)
Eng. Mecânico - Armando Wilmans Nunes da Fonseca Júnior – (Colaborador)
Eng. Mecânico -Felipe Silva Lima – (Colaborador)
Eng. Mecânico - Renan Di Pace Arruda – (Colaborador)
Comitê Científico SiPGEM 2018

Prof. Dr. Antônio Almeida Silva – UAEM
Prof. Dr. Celso Rosendo Bezerra Filho – UAEM
Prof. Dr. Renato Alexandre Costa de Santana – UABQ
Prof. Dr. Walman Benício de Castro – UAEM
Prof. Dr. Wanderley Ferreira de Amorim Júnior – UAEM
Eng. Mecânico Me. Armando Wilmans Nunes da Fonseca Júnior
Eng. Mecânico Me. Carlos Eduardo da Silva Albuquerque
Eng. Mecânico Me. José Jorge da Silva Júnior
Eng. Mecânico Me. José Ricardo Ferreira Oliveira
Eng. Mecânico Me. Paulo César Sales da Silva
APRESENTAÇÃO
Este livro foi concebido a partir de uma série de trabalhos que foram apresentados no II
SiPGEM - Simpósio de Pós-Graduação em Engenharia Mecânica - a partir de pesquisas em
desenvolvimento nos laboratórios da Unidade Acadêmica de Engenharia Mecânica da
Universidade Federal de Campina Grande (UAEM/UFCG). A proposta de melhorar a
visibilidade e inserção social do Programa de Pós-Graduação em Engenharia Mecânica teve
sua primeira edição em 2017, por iniciativa dos alunos de mestrado que se dedicaram,
contribuíram e evoluíram para essa nova formatação, que nesse ano adotou a temática
“expandindo o horizonte do jovem engenheiro”. É destinado a ser um evento científico para a
comunicação e divulgação da comunidade acadêmica local e regional e pretende envolver os
discentes e docentes dos cursos de Graduação e Pós-Graduação em Engenharia Mecânica da
UFCG e de outras IES parceiras.
Tem como objetivos fomentar a integração entre alunos, a interdisciplinaridade, a divulgação

dos resultados das pesquisas realizadas no âmbito das linhas de pesquisas em andamento na
UFCG e colaborar na formação do pesquisador em aspectos relacionados à formação de
recursos humanos de alto nível. Na sua programação inclui palestras, minicursos e seções
técnicas com artigos científicos, bem como a participação de alunos de graduação, por meio
das equipes de projetos de extensão (Grupos PET, ParahyBaja, Labem, Parahyasas e
Scuderia), através de Stands com protótipos e painéis para divulgar os vários trabalhos em
desenvolvimento na Engenharia Mecânica. Também conta com a participação de especialistas
e profissionais da indústria em outras áreas trazendo novidades em temas relevantes e atuais
visando motivar os alunos para a sua formação nesse ambiente acadêmico. Como parte dessa
jornada tem-se buscado apoios em empresas locais e órgãos de fomento e divulgação visando
aumentar o seu alcance e contribuir de maneira formal com os principais temas de pesquisa
em andamento envolvendo aspectos teóricos (métodos computacionais e analíticos) e
experimentais (pesquisas em laboratórios e trabalhos de campo), visando a expansão do
conhecimento e potencial de aplicação de novas ferramentas com interesses compartilhados
pela indústria, fortalecendo, sempre que possível, a interação da universidade com o setor
produtivo.
Enfim agradecemos à UAEM/UFCG e a Editora Poisson que se propôs publicar os resultados

desse evento na forma de capítulos de livro, além das agências que fomentam as pesquisas
(CNPq, CAPES, Fapesq-PB), às empresas parceiras e a todos os alunos de graduação e pós-
graduação envolvidos, especialmente a equipe da comissão organizadora, nesse ano
representada pelas alunas Raissa Alves Queiroga e Kaline Ventura Batista, pelo empenho e
dedicação na formatação do livro.
Antonio Almeida Silva

Coordenador do Programa de Pós-Graduação em Engenharia Mecânica - UFCG/CCT
Sumário
Capítulo 1: Mineração de criptomoedas utilizando energias renováveis.......................................
8
Agripino Pereira da Silva Filho, Heleno da Costa Neto, Maria Aparecida Souza de Andrade,
Talles Caio Linhares de Oliveira
Capítulo 2: Análise numérica de protótipo de ponte estaiada com elementos

com memória de forma superelásticos em vibração forçada ........................................................
19
Armando Wilmans Nunes da Fonseca Júnior, Marcelio Ronnie Dantas de Sá, Rômulo Pierre
Batista dos Reis, Antonio Almeida Silva
Capítulo 3: Influência dos elementos Berílio, Titânio e Boro na dureza de ligas

do sistema Cu-Al passíveis do efeito memória de forma ...............................................................
31
Carlos Cássio de Alcântara, Lucas Daniel Martins Fernandes, José Joelson de Melo Santiago
Capítulo 4: Análise numérica de uma prancha de skate fabricado em

compósito estrutural submetida a esforços transversais ...............................................................
38
Álvaro Barbosa Rocha, Armando Wilmans Nunes da Fonseca Júnior, Felipe Silva Lima,
Wanderley Ferreira de Amorim Júnior
Capítulo 5: Caracterização microestrutural e mecânica de uma liga memória

de forma do sistema Cu-Al-Mn ......................................................................................................
49
Carlos Cássio de Alcântara, Jéssica Menezes Alves, José Joelson de Melo Santiago
Capítulo 6: Análise experimental quanto aos efeitos da dilatação térmica na

gasolina automotiva e no etanol hidratado ....................................................................................
58
João Manoel de Oliveira Neto, Marcelo Bezerra Grilo
Capítulo 7: Análise numérica da influência da profundidade de operação no

63
coeficiente de arrasto do casco de um veículo autônomo submarino...........................................
João Victor Nunes de Sousa, Antonio Gilson Barbosa de Lima
Capítulo 8: Detecção dos danos de impacto em materiais compósitos via

Termografia Infravermelha .............................................................................................................
74
Kaline Ventura Batista, Renan Di Pace Arruda, Gustavo Lima do Nascimento, Wanderley
Ferreira de Amorim Júnior, Juscelino de Farias Maribondo
Sumário
Capítulo 9: Simulação de um dispositivo giroscópio absorvedor do tremor da
Doença de Parkinson ....................................................................................................................
87
Marcelio Ronnie Dantas de Sá, Armando Wilmans Nunes da Fonseca Júnior, Aline Michelly
do Nascimento, Antonio Almeida Silva
Capítulo 10: Revestimentos metálicos obtidos por eletrodeposição: Uma

revisão da literatura .......................................................................................................................
99
Raíssa Alves Queiroga, Mikarla Baía de Sousa, Josiane Dantas Costa, Arthur Filgueira de
Almeida, Renato Alexandre Costa de Santana
Capítulo 11: Estudo de caracterização mecânica em materiais processados

através de modelagem por fusão e deposição (FDM) ..................................................................
111
Ramon Sales de Araújo Batista, Carlos José de Araújo
Capítulo 12: Simulação numérica do processo de separação gás/ óleo

empregado em um vaso separador ..............................................................................................
121
Tuana Maria de Medeiros Costa, Endyara de Morais Cabral, Severino Rodrigues de Farias
Neto
Capítulo 1
Agripino Pereira da Silva Filho

Heleno da Costa Neto
Maria Aparecida Souza de Andrade
Talles Caio Linhares de Oliveira.
Resumo: O cenário mundial ultimamente tem se inclinado para o desenvolvimento

de tecnologias limpas e sustentáveis, com isso essas novas tecnologias requerem
cuidados especiais para entrar com força total no mercado. A mineração de
criptomoedas é uma prática muito comum em todo o mundo, porém no Brasil,
devido às altas taxas cobradas das concessionárias de abastecimento de energia
elétrica, não é viável. A partir disso, foi desenvolvido projeto para utilizar energia
elétrica através de placas solares que conseguem suprir a necessidade de
fornecimento de energia no processo de mineração. Esta pesquisa está dividida em
duas partes principais. A primeira parte abrange todo um estudo teórico e
comparativo e a segunda parte traz um levantamento dos equipamentos utilizados
para o uso de energia sustentável e mineração e seus respectivos valores,
indicando se o projeto é viável ou não. Após todos os cálculos realizados foi
identificado que, aqui no Brasil, a prática de mineração de criptomoedas mostrou-
se ser viável apenas com a utilização de um meio de energia sustentável, que no
presente trabalho foi a energia solar.
Palavras-Chave: Energias renováveis. Criptomoedas. Mineração. Sustentabilidade.

9
1 INTRODUÇÃO sua utilização, reduzirá as agressões ao meio

ambiente e visam retornos financeiros a curtos
As preocupações com novas tecnologias e
prazos.
aplicações estão presentes em praticamente
todos os segmentos, principalmente no setor Deste modo, o projeto visa a utilização de
industrial. Crescentemente o homem foi energia solar para a mineração de
desenvolvendo máquinas de acordo com criptomoedas, sendo realizados estudos de
suas necessidades e aumentando suas viabilidade de sua implementação e melhor
capacidades industriais. Assim, o ritmo de localização. Um comparativo entre as cidades
crescimento populacional cresceu e em de Campina Grande, Paraíba, e Cidade do
paralelo a isso o consumo de energia, Leste, localizado no Paraguai. Identificando
ignorando o cuidado com recursos e qual obterá menor prazo de retorno financeiro
consequências futuras. Há pouco tempo, os mediante as condições geográficas e
benefícios que a tecnologia vinha climáticas qual será melhor localizada para a
proporcionando como praticidade, conforto, implementação da fazenda mineradora.
redução de custos, aumento na produtividade
etc., era até então as justificativas para
realizar qualquer ação na engenharia e 2 REFERENCIAL TEÓRICO
produção. Muito pouco se falava nos danos e
2.1 ENERGIAS SUSTENTÁVEIS
impactos que poderiam ocorrer no meio
ambiente e suas consequências irreversíveis O cenário mundial tem se preocupado
(GOLDEMBERG, 2007). bastante com a questão de sustentabilidade,
devido a essa preocupação tem se falado
Portanto, a forma com que a energia vem
bastante em energias sustentáveis, fazendo
sendo produzida e consumida é discordante
com que as grandes indústrias se preocupem
com o desenvolvimento sustentável. Exemplo
em adquirir selos de sustentabilidade,
disso, são as dimensões que as criptmoedas
valorizando sua marca atraindo consumidores
vem crescendo, e há a necessidade de
que aderem a utilização de energias limpas,
máquinas com hardwares potentes e estão
algumas fontes de sustentabilidade podem
diretamente ligadas ao consumo excessivo de
estar entrelaçadas a produção de energias
energia, pois a máquina necessita estar
renováveis. Muitas fontes de energias são
sempre ligada para ser realizada a
fornecidas pelo sol, uma parte é destinada
mineração. Assim foi notado que a mineração
para a irradiação solar que fornece calor,
de criptmoedas consome mais energia do
outra forma o vento, os potenciais hidráulicos
que alguns pequenos e médios países. Com
dos rios (pela evaporação e condensação), as
essa postura, também crescem as
correntes marinhas e outras (GOLDEMBERG,
preocupações com o impacto ambiental, que
2007).
além da eletricidade, encadeia-se os
materiais para sua geração, já que os Todas as energias renováveis, como por
combustíveis fosseis são hoje seu principal exemplo, energias solares, eólicas,
fornecedor de energia elétrica (MORAIS, geotérmicas, biomassas e oceânicas são
2008). sustentáveis, pois são limpas, abundantes e
estáveis. A grande maioria de países
Adjunto ao grande consumo de energia
desenvolvidos já está buscando medidas
elétrica, está o segmento financeiro, pois no
para utilizar fontes alternativas de energia,
Brasil as taxas de energia são uma das mais
cerca de 20% das carências mundiais de
caras no mundo, principalmente no quesito
energia são provenientes de fontes
impostos, sem mencionar que as estiagens de
renováveis. A energia sustentável entrou em
chuvas baixam o nível dos reservatórios nas
cena principalmente para poder atender não
hidrelétricas que é o principal meio de
só a demanda atual de energia, mas também
geração de energia no país, e é necessário
fornecer opções para o futuro.
utilizar as termoelétricas, que é um meio caro
e altamente poluente, alterando as tarifas com Energias solares tem se destacado bastante,
bandeira vermelha e elevando drasticamente pois é uma fonte interminável de energia,
o valor final das contas de energia. devido ao aquecimento global os níveis de
irradiação solar tem aumentado
Visto as carências encontradas é de
gradativamente fazendo com que até os
fundamental importância à aplicação de
países com os mais baixos níveis de
energias sustentáveis para suprir este
irradiação se tornem grandes produtores
crescente consumo que não só estabiliza a
dessa tipo de energia. Um excelente exemplo
II Simpósio de Pós-Graduação em Engenharia Mecânica da Universidade Federal de Campina Grande: Coletânea de Artigos
10
é a Alemanha que segundo: Bruno Pereira et fundamental para o utilização das energias
al. (2017) tem um poder de irradiação 40% renováveis é a diminuição da liberação de
menor que o brasil, entretanto é um dos gases responsáveis pelo efeito estufa. Assim,
países líderes em produção de energia solar, as energias renováveis ainda tem um custo
e vem sendo utilizado há cerca de 30 anos, elevado comparado com as fontes
indicando confiabilidade e durabilidade de convencionais, mas as taxas sobre a poluição
geração. provenientes do combustíveis fosseis criaram
um nível de preços onde ambas seriam
A participação das fontes renováveis no Brasil
competitivas. Porém, a hidráulica se sobre sai
vem aumentando gradualmente, porem sua
como principal fonte de fornecimento
diversidade ainda é baixa com pode ser
energético (EPE, 2017).
notado na Figura 1. Pois, a questão
Figura 1. Estrutura da oferta interna de energia elétrica no Brasil.
65,20%
8,20% 6,80% 10,50%

0,20% 3% 2,50% 3,60%
Fonte: adaptado EPE (2017).
2.2 CRIPTOMOEDAS modelo original foi o meio termo de

criptografia e bancos de dados espalhados e
As criptomoedas têm chamando bastante
distribuídos por pontos de redes P2P
atenção por vir utilizando as criptografias em
registrando as transações e utilizando funções
controle das maiores transações mundiais,
básicas de segurança, para se certificar que
para a criação de novas moedas digitais.
os bitcoins só serão gastos pelo próprio dono
Assim, as criptomoedas usam a criptografia
eliminando gastos duplicados (SILVA e CRUZ,
encravada em alguns algoritmos
2014).
proporcionando a formação de novas
moedas. Onde, o sistema utilizado é Bitcoin tem seu funcionamento com a rede de
basicamente composto por duas chaves, uma pagamentos peer-to-peer e a moeda virtual,
pública e outra privada, na qual a privada é onde opera basicamente como dinheiro
apenas do dono e se essa chave for roubada online, não dependendo se quer de alguma
ou perdida não poderá ser resgatado um dos autoridade central, um dos principais atrativos
prejuízos da moeda. Pois, não existe um do bitcoin fazendo do mesmo um dos
órgão regulador e consequentemente com pioneiros sistemas de pagamento global
quem reclame a respeito de fraudes e roubos. completamente descentralizado. Portanto, a
Já o ponto positivo é total anonimato do criptomoeda bitcoin utiliza o sistema de
usuário (MOIA, 2014). criptografia de chave privada e pública para
registrar as efetuações de transações
Uma das mais conhecidas é a criptomoeda
realizadas pelo usuário, desta forma as
bitcoin que foi desenvolvida em 2008 por um
transações que ocorreram ficam registradas
programador chamado Satoshi Nakamoto, o
em um gigante banco de dados da rede da
11
criptomoeda denominado de blockchain como retribuição, dependendo diretamente

(SILVA e CRUZ, 2014). da máquina utilizada e que quanto melhor e
mais rápida a capacidade de processamento
As criptomoedas ou moedas digitais sem
do equipamento melhores as chances de se
dúvidas são uma das tecnologias mais
receber bitcoins (PIRES, 2017).
inovadoras e complexas abrangendo tanto
usuários leigos até profissionais na área da Com a popularidade da mineração para se
informática, e por envolver várias áreas de obter bitcoins tornou-se algo muito
diferentes conhecimentos humanos, porém já competitivo e complicado, pois cresceu
é notório que esta tecnologia vem consideravelmente o número de usuários
conquistando cada vez mais seu espaço e minerando e a rede de criptomoedas
permanecera no mercado. começou a dificultar a quantidade de geração
de moedas para não ter um aumento da
inflação da moeda, e um dos pontos que faz a
2.3 MINERAÇÃO moeda tenha grande valor é sua escassez e a
dificuldade em obtê-la. E pensando-se nisso
As maiorias das moedas aderirão ao design
foram criadas máquinas especialmente para
do bitcoin, já que é o modelo de criptomoeda
executar a mineração bitcoin, um exemplo
mais antigo. E para adquirir bitcoin, por
disso é o minerador Bitmain Antminer S9i.
exemplo, é utilizado o método de minerar
Uma gigantesca e poderosa máquina criada
bitcoin, ou seja, é utilizado um programa de
exclusivamente com a finalidade de
computador gratuito ao qual coordena a
mineração de criptomoedas.
mineração executada por meio do hardware
de um computador aos quais computadores A grande desvantagem da mineração de
no mundo inteiro realizam complexos estorvos criptmoedas é a necessidade de máquinas
matemáticos. extremamente potentes ligadas
initerruptamente para realizar a mineração,
É utilizado o termo minerar, pois descreve
consumindo excessivamente energia elétrica.
perfeitamente o que ocorre no processo. A
Segundo o Bitcoin Energy Consumption, em
mineração funciona com a seguinte
novembro de 2017 toda rede de bitcoin usou
configuração, ao iniciar o programa de
76 milhões de quilowatts-hora, um total
mineração das criptomoedas à rede inicia a
descontrole no consumo desenfreado de
utilização dos hardwares dos usuários e os
energia elétrica. Como pode ser visto na
utiliza para realizar funções matemáticas
Figura 2, em apenas um mês.
complexas, assim, o usuário contribuirá com
cada problema resolvido com a rede
aumentando a chance de adquirir bitcoins
Figura 2. Gráfico de índice de consumo de energia entre os meses de outubro e novembro de 2018.
66
Estimativa TWh por ano
64
63
62
60
09/OUT 13/OUT 17/OUT 21/OUT 25/OUT 29/OUT 02/NOV 06/NOV
Fonte: adaptado BitcoinEnergyConsumption.com
12
3 METODOLOGIA pelo serviço de fornecimento do último ano da

concessionaria local.
Diante das definições, este estudo, baseia-se
em tipos de pesquisas e trabalhos, que Diante disto, foi realizado o dimensionamento
dividem-se aos critérios de fins e meios. da quantidade necessária de placas solares
fotovoltaicas para suprir o fornecimento de
De modo que, os fins devem ser viáveis a
energia elétrica no minerador, analisando
partir de pesquisas para a melhor viabilidade
aspectos comerciais e técnicos. A escolha do
de implementação de sistemas solares a fim
equipamento solar se deu pela melhor
de obter o melhor retorno financeiro e um
eficiência de geração, reduzindo a
curto espaço de tempo, aceitando ou
quantidade de placas e a facilidade de
rejeitando o projeto. Tanto quanto aos meios,
adquirir o produtos em ambas as situações.
em virtude de pesquisas bibliográficas para a
utilização de dados retirados de, matérias de Por fim, foram utilizados métodos e analises
portais de órgãos regulamentadores, livros, de viabilidade dos projetos para a melhor
revistas cientificas, outras formas de fontes e decisão de investir ou não no projeto e em
um estudo de caso, realizando uma pesquisa qual das cidades há um real retorno e lucro
detalhada sobre o sistema a ser implantado após sua instalação.
com análises técnicas e sua real viabilidade.
A escolha da região se deu através das
necessidades e condições favoráveis que o 3.1 MATERIAIS UTILIZADOS
sistema de geração de energia requeria, no
3.1.1 MINERADOR
caso placas solares, que estão diretamente
ligadas a intensidade de irradiação solar, e Para ser ter uma viabilidade real, foi
consequentemente, realizando um selecionado um dos equipamentos de
comparativo entre cidades do Brasil e mineração mais eficientes fabricado no
Paraguai, devido as suas diferentes mercado até o presente momento, o modelo
localizações, climas, intensidade solar e da Bitmain Antminer S9i com fonte integrada,
preços de tarifas energéticas. Figura 3. Onde seu desempenho é superior
devido a sua versatilidade de minerar vários
No primeiro instante, foi escolhido o
tipos de moedas, sua velocidade de
equipamento de mineração mais eficiente no
processamento de 14 Teras hashrate por
mercado, para obter competitividade e
segundo, que aumenta a estimativa de
valores de consumo reais. Em seguida, foram
adquirir a criptmoeda mais rapidamente, tem
levantados dados referentes ao consumo de
um consumo de energia melhorado em
energia elétrica que o equipamento minerador
relação ao seu antecessor, de 850 W de
atingiria no final do mês em kWh, a partir de
potência real, e fornece a interface de
sua potência real, para a zona rural de
conectividade com outras mineradoras
Campina Grande e a Cidade do Leste, no
multiplicando o poder de processamento.
Paraguai, e seus respectivos valores pagos
Figura 3. Equipamento de mineração Bitmain Antminer S9i com fonte original.
Fonte: safemagazine.com.br/bitmain
13
3.1.2 PAINÉIS SOLARES de potência nas condições padrão, com área

externa do módulo de 1,64 m² e eficiência
A energia solar fotovoltaica é obtida através
energética de 15,9. Da fabricante Seraphin
do efeito fotovoltaico. Esse efeito se dá devido
Solar System, com certificado do INMETRO.
a alguns materiais semicondutores, a partir de
Onde a partir de comparativos há uma grande
incidência da luz solar onde os fótons da luz
variação de fatores como, eficiência,
excitam os elétrons a saltar da camada de
fabricante, local de fabricação, dimensões e
condução, que acontece em condições
entre outros.
controladas que favorecem a criação de uma
tensão e corrente elétrica. A irradiação
proveniente do sol fornece energia necessária
3.2 PROCEDIMENTOS UTILIZADOS
para que os elétrons saltem para a banda de
condução, com esse movimento entre a Após ser selecionado os equipamentos
lacuna e a banda de condução que a energia principais e ideias no mercado, foi realizado o
elétrica é gerada e armazenada pelos cálculo de consumo mensal de energia
condutores da célula (PORTOLAN, 2009). elétrica em kWh consumido pelo minerador.
Para este projeto foi utilizado placas solares Logo para encontrar o consumo mensal será
com um módulo solar fotovoltaico de 265 W utilizada a seguinte equação 1:
C = P ∗ t ∗ nd (1)
Onde:
C = consumo mensal (kWh/mês);
P = potência real (W);
t = tempo de funcionamento do equipamento (h);
nd = número de dias utilizados.
Em seguida será necessário realizar o cálculo concessionaria utilizada. Como pode ser
do valor mensal em reais com a tarifa expressa na equação 2:
energética de acordo com a região e
Gm = C ∗ Tf (2)
Onde:
Gm = Gasto mensal em reais;
C = consumo mensal (kWh/mês);
Tf = tarifa energética por kWh.
Portanto, após identificar quanto um diretamente na geração de energia elétrica
equipamento consumira mensalmente em por meio da energia solar.
kWh/mês, é preciso o dimensionamento das
A partir disso, é indispensável a obtenção de
placas solares que seriam necessárias para
dados com precisão e segurança de quantos
suprir a demanda de consumo do minerador.
painéis solares serão necessários para suprir
Porém, é preciso levar em consideração que
corretamente a exigência do consumo do
o Brasil é um país de vasta extensão territorial
equipamento.
e com diferentes regiões climáticas, assim
como, o Paraguai tem suas divergências e Foram realizadas estimativas de geração
inconformidades geográficas. Logo, alguns média que um módulo fotovoltaico, teria em
locais têm grande incidência de chuvas, ambas as regiões. Exemplo disto, em
outros já são mais áridos e possuem Campina Grande na região nordeste, como
diferentes variações de temperatura, onde pode ser visto na Figura 4. Segundo
todos esses fenômenos influenciam parâmetros estabelecidos pelo Atlas do
14
potencial energético solar e eólico, o módulo perdas totais do sistema fotovoltaico, seja por
solar fotovoltaico mais eficiente e utilizado tem aquecimento ou sujeira, em 20% da potência
potência de 265 W, sendo utilizado no nominal.
nordeste, inclinado 20º e considerado as
Figura 4. Média anual de insolação diária no Brasil.
Fonte: ATLAS Solarimétrico do Brasil. Editora

Universidade da UFPE, 2000.
Dessa maneira, na região nordeste, onde através de estações meteorológicas, Atlante

abrange a cidade Campina Grande no estado do Governo de Canárias a Espanha,
da Paraíba, tem-se em média, 4,93 horas determinando potencial energético solar,
diárias de irradiação útil à geração onde um minuto de insolação é uma medida
fotovoltaica. de energia recebida por metro quadrado por
minuto, assim o valor de insolação em horas é
Com tal característica, na Cidade do Leste no
aproximadamente 5,88 horas de irradiação
Paraguai, foi utilizada com base de
útil. Como pode ser notado na Figura 5
informação e dados globais de irradiação
(ATLAS, 2014).
solar gerados a partir de medições feitas
15
Figura 5. Quantificação da energia solar média anual por unidade de superfície, em função da
irradiação solar.
Fonte: simulação metereologica meso-escala, (WRF-ARW) realizado

em supercomputador Atlante do Governo de Canárias (Espanha).
Com dados provenientes da National Center of Atmospheric
Research dos Estados Unidos.
Após a observação isolada de irradiação solar Paraguai, não se mantém uma constância. Ao
média por ano em cada país, foi produzido inverso do Brasil, que o gráfico de irradiação
um gráfico comparando ao longo dos meses se mantem o mais constante possível, o mais
o comportamento e as intensidades em cada desejável em um sistema solar, como é
época do ano, facilitando observar que por observado na Figura 6.
mais acentuados os picos de irradiação no
Figura 6. Comparativo da irradiação solar entre Brasil e Paraguai.
Irradiação solar global

7
6
kWh/m² - dia
5
4
3
2
1
0
Brasil Paraguai
Fonte: Dados extraídos do ATLAS do potencial energético solar – ITAIPU

adaptado. 2014.
16
Portanto, considerando os parâmetros citados de insolação diária útil, é defino quanto é

em ambas as situações, um módulo gerado por dia a partir da Equação 3 a seguir:
fotovoltaico com consequência da intensidade
E=P∗t∗η (3)
Onde:
E = energia gerada;
P = Potencia nas condições padrão (W);
t = tempo médio de radiação solar útil (h);
η = rendimento do módulo com perdas totais.
4 RESULTADOS E DISCUSSÕES S9i, A Tabela 1 mostra a relação entre o

tempo de utilização do minerador e sua
4.1 CÁLCULO DE GASTOS MENSAIS DO
potência real.
CONSUMO DE ENERGIA
Portanto, através da equação 1 o consumo
4.1.1 CÁLCULO DO CONSUMO MENSAL DE
mensal de um minerador utilizado
ENERGIA
diariamente, 24 horas por dia será de 612
O cálculo do consumo de energia, foi kWh/mês.
baseado no consumo mensal de uma unidade
do minerador, modelo da Bitmain Antminer
Tabela 1: Dados do consumo mensal do minerador

Potência real (Watts) 850
Tempo de funcionamento (horas) 24
Número de dias 30
4.1.2 CÁLCULO DO VALOR MENSAL EM nos últimos 4 meses do ano devido as

REAIS estiagens será alterada para o adicional B
vermelho, encarecendo seu valor final. A
O consumo mensal em reais será realizado a
Tabela 2 mostra as diferentes taxas entre
partir das tarifas energéticas de cada região,
Campina Grande e Cidade do Leste.
onde o minerador será utilizado, notando que
na zona rural da cidade de Campina Grande,
contará com duas tarifas energéticas, pois
Tabela 2: Valores tarifários de energia em Campina Grande e Cidade do Leste

Tarifas Sem adicional Com adicional B
Campina Grande (R$/kWh) 0,34920 0,70356
Cidade do Leste(R$/kWh) 0,2251 -
A partir dos dados coletados, utilizando a uma tarifa energética baixa, obtivemos um
equação 2, teremos o gasto mensal na zona valor mensal de R$ 139,7871, assim se o
Rural de Campina Grande de R$ 213,7104 equipamento for instalado no Paraguai
nos oito primeiros meses e nos últimos 4 possuirá um gasto anual de R$ 1677,4452.
meses R$ 430,58, assim adquira um gasto Aproximadamente a metade do valor
anual de R$ 3.432,00. Consequentemente, na encontrado no Brasil.
Cidade do Leste, no Paraguai, que possui
17
4.2 CÁLCULO DA ENERGIA GERADA altera diretamente na geração de energia. Na

Tabela 3 fornece os dados necessários.
Como em ambas as regiões tem tempos
médios de irradiação solar úteis diferentes,
Tabela 3: Dados médios de irradiação solar útil

Cidades analisadas Campina G. C. Leste
Potência do módulo fotovoltaico (Watts) e rendimento com perdas (%) 250/80 250/80
Tempo médio de radiação solar (horas) 4,93 5,88
Deste modo, utilizando a equação 3, necessários 20 módulos, e na Cidade do

fornecera que a energia gerado em um dia na Leste seriam necessários 17 módulos.
cidade de Campina grande é de 1,045
kWh/dia e gerada por mês será de 31,35
kWh/mês. Seguidamente, a mesma placa 4.2.1 CUSTO DE INSTALAÇÃO E
solar instalada na Cidade do Leste fornecera EQUIPAMENTOS
1,2466 kWh/dia e gerada por mês será de
37,40 kWh/mês.
Após adquirimos todos os valores necessários
Por fim, visto que o equipamento minerador
de instalação em um equipamento de
consome 612 kWh/mês, é utilizado uma regra
mineração utilizando energia solar, será
de Três simples para saber aproximadamente
realizado os cálculos para uma fazenda de
quantos painéis solares são necessários para
mineração com 20 mineradores de
gerar e suprir o consumo elétrico do
criptomoedas. A Tabela 4 fornece as
equipamento. Logo, o painel fotovoltaico
quantidades que serão utilizadas e seus
instalado em Campina Grande seriam
respectivos valores.
Tabela 4: Valores e quantidades necessárias de equipamentos para instalação

Equipamento Quantidade (unid.) Valor p/ unidade (R$) Valor Total(R$)
Placas solares 400 650,00 260.000,00
Minerador 20 5.880,00 117.600,00
Com base em comparativos no setor de como a fazenda de mineração foi

instalação, foi adquirido o valor da mão de dimensionada com 20 mineradoras o lucro
obra e acessórios dos equipamentos cotados mensal será de R$ 11.047,40 e anualmente de
em média no valor de R$ 60.000,00. Portanto, R$ 132.568,88.
o custo total inicial com, os mineradores,
placas solares e instalação foi de
aproximadamente R$ 437.600,00. 5 CONCLUSÃO
Neste estudo, foram analisados duas regiões
distintas com a finalidade de obter um
4.3 RECEITA EM CRIPTOMOEDAS
comparativo entre ambas para a melhor
Os retornos financeiros não são instalação de uma fazenda de mineração de
necessariamente apenas com a obtenção das criptmoedas, como a zona rural de Campina
criptomoedas, e sim através da mineração. Grande, no estado da Paraíba, e a Cidade do
Exemplo disso, é a mineração da moeda Leste, localizada no Paraguai. Após analises
Bitcoin. Atualmente, essa moeda escolhida iniciais nas tarifas de energia elétrica, foi visto
consegue gerar aproximadamente 0.01966 que nesse quesito há dois extremos
Bitcoin por mês, utilizando o equipamento diferentes, onde o Brasil têm uma das tarifas
utilizado na pesquisa, convertendo fornecerá de energia elétrica mais caras no mundo e o
aproximadamente R$ 552,37 por mês, Paraguai entre as mais baixas na américa do
considerando a taxa de conversão. Ou seja, sul, sem mencionar que nos últimos quatro
18
meses do ano o Brasil sofre com as estiagens, concessionaria em Campina Grande, o

alterando a tarifa com adicional B Vermelho, retorno financeiro do projeto seria em
encarecendo ainda mais seus valores. Em aproximadamente 2,6 anos.
seguida, foram analisadas as questões de Consecutivamente, utilizando o sistema de
instalação de placas solares, que estão energia solar, o retorno financeiro do projeto
diretamente ligadas as condições seria em aproximadamente 3,2 anos.
geográficas, climáticas e especialmente a
Por fim, há uma grande viabilidade de
irradiação solar. Onde, foi visto que em cada
instalação do sistema solar para a mineração
região a um tempo médio de irradiação útil,
de criptomoedas na zona rural de Campina
que determina efetivamente o desempenho
grande, por mais que a estimativa de tempo
da instalação e utilização do sistema solar.
utilizando a concessionaria elétrica seja
Portanto, na Cidade do Leste, no Paraguai, o menor, consequentemente irá variar, se
investimento com as mineradoras e a observado o histórico do crescente aumento
utilização da energia elétrica da das tarifas de energia. Além de vários pontos
concessionaria o retorno financeiro do projeto positivos como, nos últimos quatro meses do
será em 1,9 anos. Assim, fica descartada a ano, onde é a época crítica de escassez nos
instalação de um sistema de energia solar, já níveis dos reservatórios das hidrelétricas e há
que o retorno do investimento é breve, e os o aumento da tarifa, é observado os maiores
impostos são baixos, e não há histórico de índices de irradiação solar, favorecendo sua
grandes variações nos valores de taxas e geração, o equipamento tem vida útil de 25
tributos. anos com baixíssimas manutenções e além
de benefícios financeiros é a importante
À vista disso, com a compra das mineradoras
contribuição do uso da energia solar a
e a utilização da energia elétrica da
consciência socioambiental.
REFERÊNCIAS
[1] ATLAS. “Solar mede energia do sol [6] Pires, H. F. Bitcoin: a moeda do
disponível em qualquer ponto do estado”. ciberespaço. Geousp – Espaço e Tempo (Online),
Disponível em: (https://www.itaipu.gov.br/sala-de- v. 21, n. 2, p. 407-424, agosto. 2017. ISSN 2179-
imprensa/noticia/atlas-solar-mede-energia-do-sol- 0892.
disponivel-em-qualquer-ponto-do-estad).
[7] Portolan, I. S; “Integração de painéis
Acessado em 30 de Outubro de 2018.
solares fotovoltaicos em edificações residenciais e
[2] Daniel P Et Al. “Paradigmas da energia sua contribuição em um alimentador de energia de
solar no Brasil e no mundo”. Revista Eletrônica Em zona urbana mista”. Projeto de mestrado em
Gestão, Educação E Tecnologia Ambiental, V. 20, Engenharia Civil – Universidade Federal de Santa
N. 1, P. 241-247, 2016. Catarina, UFSC. Florianópolis - SC. PPGEC. 2009.
[3] Goldemberg, José; Lucon, Oswaldo. [8] Portal Solar LTDA.” Quanto Custa para
“Energias renováveis: um futuro sustentável”. instalar energia solar”. Disponível em:
REVISTA USP, São Paulo, n.72, p. 6-15, (https://www.portalsolar.com.br/quanto-custa-para-
dezembro/fevereiro 2006-2007. instalar-energia-solar.html) Acessado em 24 de
Outubro de 2018.
[4] Moia, G. H. V; Henquires, M. A. A.
“Avaliação da segurança de protocolos [9] Silva, M. W. S; Cruz, A. R. Bitcoin a moeda
criptográficos usados em moedas virtuais”. digital que se tornou realidade. Publicação em
Disponível em: revista cientifica da UNESC. Curso de Análise e
(http://www.fee.unicamp.br/sites/default/files/depart Desenvolvimento de Sistemas – UNESC.
amentos/dca/artigos/a264v0_IESS.pdf) Acesso em: Rondônia/Cacoal. v. 12, n. 15 (2014).
26 de setembro de 2018.
[10] Smaal, B. Techtudo. Bitcoin: “A mineração
[5] Morais, J. F; Silva, W. S; Rabelo, P. F. R. de moedas”. Disponível em:
“Impactos ambientais provenientes das novas (http://www.techtudo.com.br/artigos/noticia/2014/01
tecnologias de telecomunicações”. Projeto de /bitcoin-a-mineracao-de-moedas.html) Acessado
mestrado em Engenharia Civil – Universidade em: 30 de Agosto de 2018.
Federal de Fluminense, UFF, Niterói – RJ, 2008.
19
Capítulo 2
Armando Wilmans Nunes da Fonseca Júnior

Marcelio Ronnie Dantas de Sá
Rômulo Pierre Batista dos Reis
Antonio Almeida Silva.
Resumo: Entre os vários materiais utilizados atualmente para supressão de

vibrações, estão as ligas com memória de forma (LMF). Estas vêm ganhando
popularidade graças a sua capacidade de sofrer grandes deformações reversíveis,
aliadas às suas propriedades de dissipação de energia. Neste contexto, esse
trabalho tem como objetivo avaliar o uso de uma metodologia para estimar a
resposta numérica em vibração forçada de uma estrutura do tipo ponte estaiada,
em escala reduzida, com elementos com memória de forma superelásticos atuando
como elementos dissipadores de energia. Os resultados mostraram boa coerência
na determinação dos parâmetros modais da estrutura. Com maior erro associado
ao segundo modo de vibrar do sistema. Concluiu-se que os modelos numéricos
permitem uma boa previsão da resposta da estrutura.
Palavras-chave: Análise numérica, Elementos Finitos, Ligas com Memória de

Forma, Ponte Estaiada.
20
1 INTRODUÇÃO superelasticidade (SE) (OTSUKA; WAYMAN,

1998). O comportamento superelástico (ou
Ligas com memória de forma (LMF) são ligas
pseudoelástico) das LMFs está associado à
metálicas que sofrem transformações entre
capacidade destas de sofrer altas
fases sólidas induzidas por mudanças de
deformações induzidas por tensões a partir
temperatura e/ou tensões apropriadas e que
de carregamento aplicado e recuperar-se
podem recuperar deformações
completamente destas mediante
aparentemente permanentes. Tais ligas
descarregamento. A Figura 1 ilustra o
possuem uma série de propriedades
comportamento pseudoelástico típico de
desejáveis, destacando-se o efeito de
LMFs através de um diagrama tensão-
memória de forma (EMF) e a
deformação.
Figura 1: Carregamento cíclico típico de uma LMF superelástica (Lagoudas, 2008).
Pode-se notar que, durante o ciclo que desenvolvem modelos numéricos com o
superelástico, as tensões de transformação propósito de prever o comportamento dessas.
direta (σMs → σMf) e reversa (σAs → σAf) são Um dos modelos mais usados atualmente é
diferentes. Isto resulta em uma histerese de modelo constitutivo proposto por Auricchio e
tensão mecânica. O laço da histerese Taylor (1997), usado em análises de
depende do material da LMF e das condições estruturas de LMF pelo método dos elementos
de teste utilizadas. finitos. Entretanto, esse modelo não leva em
consideração as propriedades dinâmicas do
Fisicamente, essa histerese representa a
material, sendo apenas utilizado para
energia recebida pelo material para sofrer a
aplicações estáticas ou quase estáticas.
deformação. Durante o carregamento, a liga
absorve energia para deformar-se e, mediante Neste contexto, esse trabalho tem como
descarregamento, apenas parte dessa objetivo avaliar o uso de uma metodologia
energia retorna ao sistema, sendo o restante que combina simulação e validação da
dissipado para o meio na forma de calor. Isto análise experimental para estimar a resposta
implica em dizer que o material tem numérica em vibração forçada de uma
capacidade de dissipar uma parcela da estrutura do tipo ponte estaiada, em escala
energia utilizada para realizar o ciclo reduzida, com elementos com memória de
superelástico. Logo, quanto maior a histerese, forma superelásticos atuando como
maior a energia dissipada por cada ciclo. É elementos dissipadores de energia utilizando
devido a esta propriedade que as LMF o método dos elementos finitos. Comparando
superelásticas estão sendo cada vez mais os resultados obtidos numericamente com
utilizadas em dispositivos de controle passivo aqueles obtidos através do protótipo
de vibrações. experimental.
Devido as essas propriedades, as LMF são
objetos de estudo de muitos pesquisadores,
21
2 METODOLOGIA
A metodologia para a realização deste
trabalho é ilustrada na Figura 2.
Figura 2: Metodologia para obtenção da resposta numérica do protótipo de ponte estaiada com
elementos com memória de forma em vibração forçada.
NECESSIDADE: Obter resposta da estrutura

com elementos de LMF em vibração forçada
utilizando métodos numéricos
INÍCIO DA
PESQUISA
MODELAR ESTRUTURA ADICIONAR ESTAIS À

FASE 1 FASE 2
SEM ESTAIS MODELAGEM
Caracterizar simulação com
Realizar ensaios com
ETAPA 1.1 parâmetros geométricos e de ETAPA 2.1
elementos dissipadores
materiais
NÃO NÃO
Alimentar modelagem com Alimentar simulação com
ETAPA 1.2 fatores de amortecimento ETAPA 2.2 parâmetros dos elementos
obtidos experimentalmente dissipadores
Resposta da estrutura Resposta da estrutura

Adequadas? sem elementos SIM Adequadas? com elementos
dissipadores dissipadores
SIM
FASE 3 COMPARAR RESULTADOS
As atividades foram desenvolvidas em três forçada. Essa foi executada utilizando o

fases que serão melhores detalhadas nos método dos elementos finitos através do
tópicos a seguir. software comercial Ansys® Mechanical APDL
(Ansys Parametric Design Language) Release
16.0, mais especificamente por meio da
2.1 FASE 1 – MODELAGEM DA ESTRUTURA ferramenta de análise harmônica completa,
SEM ESTAIS que consiste na resolução da Eq. 1.
Na Fase 1 é realizada a modelagem numérica
inicial da estrutura sem estais em vibração
[𝑀]{𝑢̈ } + [𝐶]{𝑢̇ } + [𝐾]{𝑢} = {𝐹(𝑡)} = 𝐹0 𝑐𝑜𝑠⁡(𝜔𝑡) (1)
Onde:
• [C] – matriz de
• [M] – matriz de massa; • [K] – matriz de rigidez;
amortecimento;
• {𝑢̈ } – vetor nodal de • {𝑢̇ }– vetor nodal de • {𝑢} – vetor nodal de
aceleração; velocidade; deslocamento;
• {F(t)} – vetor de força.
Esta fase se subdivide em duas etapas que geométricos e de material similares ao do

serão discutidas abaixo. protótipo experimental, como ilustrado na
Tabela 1.
Na Etapa 2.1 o modelo numérico é
caracterizado com os parâmetros
22
Tabela 1 Parâmetros geométricos e de material utilizados na modelagem numérica.

Parâmetros Valores
Comprimento 1,5 m
Largura 0,25 m
Espessura 0,002 m
Modulo de Elasticidade (E) 71 GPa
Coeficiente de Poisson (ν) 0,33
Densidade (ρ) 2770 kg/m³
A estrutura que representa o tabuleiro da realizado, uma força harmônica foi aplicada
ponte foi modelada através de uma viga bi no ponto de 0,375 m do comprimento da
apoiada, utilizando elementos do tipo casca estrutura, na lateral da mesma, como ilustrado
(SHELL181). Além disso, com o propósito de na Figura 3.
simular condições similares ao experimento
Figura 3: Modelo numérico do tabuleiro da estrutura de ponte.
Aplicação
da força
Na Etapa 2.2, com o intuito de realizar uma via método de largura de banda em meia
melhor caracterização baseado no protótipo potência. O modelo foi também alimentado
experimental, o modelo numérico, sem com os diferentes níveis de força excitadora
elementos dissipadores, foi alimentando com associada a cada modo de vibrar da
os fatores de amortecimento viscosos estrutura. Estes valores podem ser vistos na
encontrados experimentalmente, calculados Tabela 2.
Tabela 2: Fatores de amortecimento e força excitadora associados aos 3 primeiros modos de vibrar
da estrutura.
1º Modo (Flexão) 2º Modo (Flexão) 3º Modo (Flexão)
Fator de amortecimento viscoso (ζ) 0,039 0,0034 0,0036
Força excitadora (N) 0,195 0,120 0,240
Apenas os três primeiros modos de vibrar da vibração. A partir do quarto modo de vibrar, a
estrutura foram escolhidos devido a estrutura já não apresentava níveis
apresentarem maiores amplitudes de significativos de vibração.
23
A conclusão desta fase permite avaliar o numérica realizada na Fase 1, mais

comportamento do modelo numérico sem os especificamente sob a forma do elemento
elementos dissipadores e comparar este com mola/amortecedor COMBIN14, que possui
o protótipo experimental. associado a si uma constante de rigidez e
amortecimento viscoso. O novo modelo
numérico é ilustrado na Figura 4. Esta Fase se
2.2 FASE 2 – ADIÇÃO DOS ESTAIS À subdivide em duas etapas que serão
MODELAGEM melhores discutidas nos próximos parágrafos.
Na Fase 2 os estais (atuando como elementos
dissipadores) são adicionados à modelagem
Figura 4: Modelo numérico com elementos COMBIN14.
Na Etapa 2.1 são realizados os ensaios de impostos aos elementos dissipadores.

força – deslocamento dos elementos Levando em consideração a pré-carga de
dissipadores em um analisador dinâmico 6,25 mm que tinha como intuito fazer com que
(Dynamic Mechanical Analyzer, DMA) modelo os elementos se encontrassem no centro de
Q800 da fabricante TA instruments. O seu laço histerético, foi possível também
propósito destes foi o de caracterizar os calcular a elongação que estes irão sofrer em
elementos dissipadores e identificar suas regime de trabalho. Todos estes valores estão
propriedades de rigidez e dissipação de ilustrados na Tabela 3.
energia para diferentes níveis de
Os ensaios de força – deslocamento dos
deslocamentos impostos.
elementos dissipadores são ilustrados na
Três diferentes deslocamentos, um associado Figura 5.
a cada modo de vibrar da estrutura, foram
Tabela 3: Níveis de elongação sofridos pelos elementos dissipadores associado a cada modo de
vibrar da estrutura.
Elongação do elemento (mm) 2,5 3,6 1,38
Elongação do elemento com pré-carga (mm) 6,25 ± 1,25 6,25 ± 1,80 6,25 ± 0,69
24
Figura 5: Ensaios de força - deslocamento dos elementos dissipadores associados aos diferentes
modos de vibrar da estrutura.
Na Etapa 2.2, com a conclusão dos ensaios (hLMF) dos elementos dissipadores utilizando
de força – deslocamento, é possível então as Eq. 2 e 3, respectivamente. Estes valores
obter os valores de rigidez equivalente (kLMF) são apresentados na Tabela 4.
e constante de amortecimento histerético
𝐹𝑚𝑎𝑥 − ⁡ 𝐹𝑚𝑖𝑛
𝑘𝐿𝑀𝐹 = ⁡ (2)
𝑋𝑚𝑎𝑥 − ⁡ 𝑋𝑚𝑖𝑛
𝐴𝑙𝑎ç𝑜
ℎ𝐿𝑀𝐹 = ⁡ (3)
𝜋𝑋 2
Tabela 4: Valores de rigidez e constante de amortecimento histerético dos elementos dissipadores.

Rigidez equivalente (kLMF) 367,38 N/m 334,89 N/m 463,37 N/m
Constante de amortecimento histerético (hLMF) 30,48 N/m 27,25 N/m 25,90 N/m
Com a conclusão dessa fase é possível então primeiramente foi realizada uma análise modal
analisar e validar a resposta numérica da numérica no software Ansys APDL. Os modos
estrutura com elementos dissipadores em de vibrar e as frequências naturais da
vibração forçada. estrutura associadas a cada um destes pode
ser visualizado na Figura 6.
3 RESULTADOS E DISCUSSÕES
3.1 ANÁLISE MODAL DO MODELO SEM
ELEMENTOS DISSIPADORES
Com o objetivo de identificar a semelhança
entre o modelo numérico e experimental,
25
Figura 6: Modos de vibrar da estrutura sem elementos dissipadores.
Devido à ausência do modo de vibrar em Foi possível então realizar uma comparação
torção no modelo experimental, este modo entre as frequências naturais experimentais
não foi levado em consideração nas análises da estrutura e aquelas obtidas via modelagem
subsequentes. numérica, como ilustrado na Tabela 5.
Tabela 2: Frequências naturais numéricas e experimentais da estrutura sem elementos dissipadores.

Modelo numérico 2,044 Hz 8,216 Hz 18,585 Hz
Protótipo experimental 2,063 Hz 8,563 Hz 19,250 Hz
Nota-se uma boa coerência entre os valores elementos dissipadores, utilizando a rigidezes
numéricos e experimentais, com o maior erro dos elementos encontrados por meio dos
associado ao terceiro modo de vibrar em ensaios de força-deslocamento, citados na
flexão da estrutura, com erro relativo de 4,8%. seção 2.2 deste trabalho. Tendo em vista que
Isto ocorre devido ao modelo numérico não a análise modal numérica só leva em
possuir uma representação adequada da consideração um valor de rigidez do elemento
rigidez da estrutura para os modos de vibrar dissipador, o valor escolhido foi aquele
mais elevados. associado ao primeiro modo de vibrar do
sistema (kLMF = 367,38 N/m), pois este era o
mais suscetível a quaisquer mudanças de
3.2 ANÁLISE MODAL DO MODELO COM parâmetros.
ELEMENTOS DISSIPADORES
Os modos de vibrar e as frequências naturais
Um procedimento similar ao realizado no do novo sistema podem ser visualizados na
tópico anterior foi feito para a estrutura com os Figura 7.
26
Figura 7: Modos de vibrar e frequências naturais da estrutura com elementos dissipadores.
Mais uma vez, é possível realizar uma modelo numérico e do protótipo experimental.
comparação entre as frequências naturais do Estes podem ser visualizados na Tabela 6.
Tabela 3: Frequências naturais numéricas e experimentais da estrutura com elementos dissipadores.

Modelo numérico 3,76 Hz 9,44 Hz 18,89 Hz
Protótipo experimental 3,81 Hz 9,50 Hz 19,6 Hz
Pode-se perceber que o maior erro relativo foi 3.3 RESPOSTA EM VIBRAÇÃO FORÇADA DO
de 4,13%, associado ao terceiro modo de MODELO SEM ELEMENTOS DISSIPADORES
vibrar da estrutura, proveniente da dificuldade
Na análise harmônica completa (Eq. 1), a
dos modelos numéricos de representar os
matriz de amortecimento tem a seguinte
modos de vibrar mais elevados.
definição:
27
𝑁𝑚𝑎 𝑁𝑚𝑏 𝑁𝑒
1 2 1
[𝐶] = ⁡𝛼[𝑀] + (𝛽 + 𝑔) [𝐾] + ⁡ ∑ 𝛼𝑖𝑚 [𝑀𝑖 ] ⁡ + ⁡ ∑ [(𝛽𝑗𝑚 + 𝑚𝑗 + 𝑔𝑗𝐸 )[𝐾𝑗 ]] + ⁡⁡ ∑[𝐶𝑘 ]⁡ (4)
𝛺 𝛺 𝛺
𝑖=1 𝑗=1 𝑘=1
Onde:
α – amortecimento α de
[C] – matriz de amortecimento; [M] – matriz de massa;
Rayleigh;
g – coeficiente de
Ω – frequência de excitação (em
β – amortecimento de β Rayleigh; amortecimento estrutural
rad/s);
constante;
Nma – número de materiais 𝛼𝑖𝑚 – amortecimento α de

[K] – matriz de rigidez; com amortecimento α de Rayleigh para material i;
Rayleigh;
Nmb – número de materiais 𝛽𝑖𝑚 – – amortecimento de β

[Mi] – parcela da matriz de massa Rayleigh para o material j;
com amortecimento β de
estrutural associada ao material i;
Rayleigh;
𝑔𝑗𝐸 – coeficiente de
mj – coeficiente de amortecimento amortecimento estrutural para o [Kj] – parcela da matriz de
estrutural constante para o material material j; rigidez estrutural associada ao
j; material j;
Ne – número de materiais com [Ck] – matriz de amortecimento

amortecimento especificado; de elemento.
Logo, os amortecimentos encontrados na transformar em amortecimentos β de

seção 2.1 não puderam ser alimentados Rayleigh, associados as frequências naturais
diretamente na simulação, sendo necessário do modelo, utilizando a Eq. 5.
𝜁
𝛽 =⁡ (5)
𝜋𝑓
Utilizando as frequências naturais obtidas na cada modo de vibrar da estrutura, como

seção 3.1, foi possível obter os valores de ilustrado na Tabela 7.
amortecimento β de Rayleigh associados a
Tabela 7: Amortecimentos β de Rayleigh associados aos modos de vibrar do modelo numérico.

Amortecimento β de Rayleigh 6 x 10-3 s 1,25 x 10-4 s 5,91 x 10-5 s
Com posse dessas informações e utilizando resposta em frequência da estrutura sem

os valores de força excitadora encontrados na elementos dissipadores em vibração forçada,
seção 2.1 deste trabalho, é possível obter a como ilustrado na Figura 7.
28
Figura 3: Reposta do modelo numérico sem elementos dissipadores em vibração forçada.
A presença do terceiro pico na faixa de 17 Hz primeiros modos de vibrar em flexão do

ocorre devido a primeira frequência natural modelo numérico e protótipo experimental,
em torção do modelo numérico. Este modo sem a presença de elementos dissipadores.
não foi observado no protótipo experimental. Estes valores podem ser visualizados na
Tabela 8.
É possível então realizar uma comparação
das amplitudes de deslocamento nos três
Tabela 8: Amplitudes de deslocamento do modelo numérico e protótipo experimental associados

aos três primeiros modos de vibrar em flexão.
Modelo numérico 6,726 mm 5,730 mm 1,344 mm
Protótipo experimental 7, 082 mm 5,048 mm 1,150 mm
Nota-se uma boa coerência entre os valores 3.4 RESPOSTA EM VIBRAÇÃO FORÇADA DO
da modelagem numérica e do protótipo MODELO COM ELEMENTOS DISSIPADORES
experimental, com o maior erro relativo de
Tendo em vista que o elemento COMBIN14
16% associado ao terceiro modo de vibrar em
tem associado a si uma constante de
flexão. Isto acontece devido a problemas da
amortecimento viscoso, é necessário
modelagem numérica de representar de
transformar o amortecimento histerético dos
forma fiel os modos de vibrar mais elevados.
elementos dissipadores encontrados na
seção 2.2 deste trabalho. Isso é realizado
utilizando a Eq. 6.
ℎ𝐿𝑀𝐹
𝑐𝑒𝑞 = ⁡ (6)
2𝜋𝑓
Onde ceq é a constante de amortecimento elementos dissipadores através da análise

viscoso equivalente, hLMF é a constante de modal numérica (Seção 3.2), é possível
amortecimento histerético e f é a frequência realizar o cálculo do amortecimento viscoso
de interesse (em Hz). equivalente destes últimos utilizando a Eq. 6 e
os dados da Tabela 4. Estes valores podem
Com a obtenção das frequências naturais que
ser visualizados na Tabela 9.
a estrutura terá com a implementação dos
29
Tabela 9: Coeficiente de amortecimento viscoso equivalente dos elementos dissipadores associados

a cada frequência natural.
Constante de amortecimento viscoso (ceq) 1,287 N.s/m 0,459 N.s/m 0,219 N.s/m
Com a aquisição desses valores, o modelo

numérico é alimentado e são realizadas novas
análises harmônicas, ilustradas na Figura 8.
Figura 4: Resposta do modelo numérico com elementos dissipadores em vibração forçada.
Com esta análise é possível mais uma vez modos de vibrar em flexão do modelo
realizar a comparação das amplitudes de numérico e do protótipo experimental, como
deslocamento associadas aos três primeiros ilustrado na Tabela 10.
Tabela 10: Amplitudes de deslocamento do modelo numérico e protótipo experimental com

elementos dissipadores.
Modelo numérico 1,70 mm 1,42 mm 1,06 mm
Protótipo experimental 1,66 mm 2,35 mm 0,98 mm
Utilizando os valores de rigidez e No entanto, pode-se observar uma grande

amortecimento dos elementos dissipadores, é defasagem dos valores de amplitude no
possível notar boa coerência do primeiro e segundo modo de vibrar, com o maior erro
terceiro modo de vibrar em flexão da relativo de 39%. Em todas as modelagens
modelagem numérica com a resposta do realizadas, pôde-se notar que o segundo
protótipo experimental, com um erro relativo modo de vibrar é o mais suscetível a
máximo de 8% no terceiro modo de vibrar. quaisquer variações de parâmetros, isso se
deve, em parte, à como o experimento foi
realizado. Os pontos de captação de sinal de
30
resposta e de aplicação de força têm o maior dissipadores se mostrou muito próxima

impacto no segundo modo de vibrar. Além daquela obtida a partir do protótipo
disso, as aproximações e linearizações que experimental. Os resultados numéricos
foram realizadas, torna a modelagem deste utilizando os parâmetros dos elementos
modo mais complexa. dissipadores se mostraram muito próximos do
comportamento apresentado pelo protótipo
experimental, para o primeiro e terceiro modo
4 CONCLUSÕES de vibrar, com maior erro relativo de 10%. Isto
indica a capacidade que os modelos
A modelagem do sistema e comparação dos
numéricos têm de prever a resposta da
resultados numéricos iniciais descreveram de
estrutura, mediante conhecimento dos
forma coerente os parâmetros de frequência
parâmetros de trabalho dos elementos
natural e modos de vibrar do protótipo
dissipadores. Houve uma dificuldade de
experimental, com um erro máximo da ordem
modelar o segundo modo de vibrar da
de 4,8%, associado ao terceiro modo de
estrutura, decorrente de como o experimento
vibrar da estrutura.
em vibração forçada foi executado e das
A resposta numérica em vibração forçada do aproximações e linearizações que foram
modelo da estrutura sem elementos realizadas.
REFERÊNCIAS
[1] AURICCHIO F., TAYLOR R.L., Shape- [3] FONSECA JÚNIOR, A. W. N.;
memory alloy modeling and numerical simulations DAMASCENO, F. O.; SILVA, A. A., Modelagem e
of the finite-strain superelastic behavior, Computer Análise Dinâmica de Estruturas do Tipo Ponte de
Methods of Application Mechanics and Vigas Incorporando Elementos com Memória de
Engineering, 143, 175-194, 1997. Forma. In: X Congresso Nacional de Engenharia
Mecânica, 2018. Anais do X Congresso Nacional
[2] FONSECA JÚNIOR, A. W. N.. Modelagem
de Engenharia Mecânica.
e Análise de Protótipo de Ponte Estaiada sob
Cargas Dinâmicas Incorporando Molas de Nitinol [4] LAGOUDAS, D. C. Shape Memory Alloys:
Superelásticas para Supressão de Vibrações. Modeling and Engineering Applications. New York:
2018, 182 f. Dissertação. (Mestrado em Engenharia Springer, 2008. v. 1
Mecânica). Universidade Federal de Campina
[5] OTSUKA, K.; WAYMAN, C. M. Shape
Grande, Campina Grande, Paraíba, 2018.
Memory Materials. Cambridge: Cambridge
University Press, 1998.
31
Capítulo 3
Carlos Cássio de Alcântara

Lucas Daniel Martins Fernandes
José Joelson de Melo Santiago
Resumo: As ligas memória de forma têm potencial para diversas aplicações, desde
as mais simples até aplicações complexas. Dentre os vários tipos de aplicações
destacam-se os dispositivos utilizados nas áreas médicas e odontológicas. Neste
contexto o objetivo desde trabalho foi a elaboração e caracterização de ligas
memória de forma do sistema Cu-Al contendo berílio e titânio-boro nos seguintes
percentuais: Cu-11,8%Al (Liga 1), Cu-11,8%Al-0,6% Be (Liga 2) e Cu-11,8%Al-
0,6%Be-0,5%Ti-0,1%B (Liga 3) (% em peso). As ligas foram elaboradas em forno
sob atmosfera ambiente, tratadas termicamente e então caracterizadas por ensaios
de microscopia óptica e ensaio de dureza Brinell. Constatou-se que a adição de
pequenos teores de Berílio, altera as temperaturas de transformação devido à
mudança da composição da matriz, visto que ocorre uma redução do tamanho
médio do grão, e consequentemente da dureza. Enquanto que os precipitados da
liga mãe Titânio-Boro inibem o crescimento do grão.
Palavras-chave: Ligas memória de forma. Liga Cu-Al-Be-Ti-B. Efeito memória de

forma.
32
1 INTRODUÇÃO através da variação das temperaturas dos

materiais, visto que dependendo do
É notório perceber que o desenvolvimento e o
tratamento térmico elas podem apresentar
avanço da humanidade desde os primórdios
memória de forma, fazendo das LMFs.
vêm atrelado às habilidades dos seus
integrantes de manipular ferramentas para As LMFs são metais que têm a capacidade de
satisfazer suas necessidades. Segundo retornar à sua forma original, após sofrerem
Salluste (1971), os primeiros povos tinham uma deformação pseudoplástica, quando
disponíveis uma pequena quantidade de aquecidas a uma determinada temperatura
materiais, tendo como destaque os materiais (CALLISTER, 2008).
encontrados na natureza, tais como pedra,
Wilkes, Liaw e Wilkes (2000), dizem que as
argila e pele de animais.
princiapsi LMFs são as ligas de Cu-Al-Ni e Ni-
Os primeiros relatos históricos da utilização Ti, devido as suas propriedades, embora
do cobre datam da época de transição do existam várias outras LMF, tais como a liga
período Paleolítico para o Neolítico, no qual o Cu-Al-Be e Cu-Al-Mn. Sun e Huang (2009)
homem deixou de ser um indivíduo nômade e dizem que as LMFs podem existir em duas
passou a ser sedentário, ou seja, o homem fases diferentes, com três diferentes
deixou de mudar de localização estruturas cristalinas, martensita maclada,
periodicamente e passou a fixar-se em martensita e austenita.
determinadas regiões. O cobre foi um dos
Nesse sentido, para Padilha (1997) os
fatores primordiais para que essa fixação
materiais metálicos que possuem essa
tenha tornando-se possível, visto que esse
interessante propriedade funcional, podem
elemento acabou substituindo a pedra na
ser empregados nos mais diversos
fabricação de ferramentas, armas, recipientes
segmentos da Indústria metalúrgica, que vão
e utensílios.
desde aplicações aeroespaciais,
Inicialmente a utilização do cobre e alumínio eletroeletrônica, construção civil, robótica a
era restrita apenas ao estado sólido, visto que bioengenharia.
o ponto de fusão dos mesmos é em torno de
Nas últimas décadas, as ligas memória de
1083ºC e 660,3ºC, respectivamente. A fusão
forma têm sido aplicadas na confecção de
de tais materiais tornou-se possível apenas
vários dispositivos comerciais, tais como:
com a invenção das fornalhas. Com o tempo o
acoplamento de tubos, fixadores mecânicos
homem descobriu que as propriedades de um
(ZAK; KNEISSL; ZATULSKIJ, 1996), válvulas,
material podem ser alteradas por meio de
implantes, aparelhos médicos e ainda
tratamentos térmicos e pela adição de outros
atuadores mecânicos (CZECHOWICZ;
elementos. Com isso foi possível a fundição
LANGBEIN, 2015).
de cobre adicionando outros metais, obtendo
as ligas metálicas à base de cobre com o Neste trabalho foram analisadas as influências
intuito de melhorar algumas características da adição de Berilio e da liga mãe de Ti-Bo no
como resistência mecânica, ductilidade, tamanho do grão, dureza e fases presentes
estabilidade térmica, tenacidade e dureza. na matriz de ligas a base de Cu-Al visto que
de acordo com Garcia, Spim e Santos (2012),
a adição de elementos em ligas metálicas
Segundo Heinz (1995), as ligas à base de Cu
pode modificar as propriedades mecânicas
possuem baixo custo de obtenção e
das mesmas.
conformação em relação às outras ligas
memória de forma (LMF). Dentre as várias
combinações de ligas à base de Cu,
2 MATERIAIS E MÉTODOS
destacam-se as ligas Cu-Zn (latões), Cu-Sn
(bronze), Cu-Ni, Cu-Ni-Zn, Cu-Al, Cu-Al-Be. As atividades desenvolvidas foram
inteiramente realizadas nas instalações da
Para Nava (2013), tais ligas possuem uma
Universidade Federal Rural do Semi-Árido
vasta gama de aplicações, como juntas,
(UFERSA). A metodologia adotada está
engrenagens, parafusos, tubulações,
apresentada no fluxograma da Figura 1.
instrumentos cirúrgicos, elementos funcionais
de sensores e atuadores que funcionam
33
Figura 1- Fluxograma da sequência metodológica
Preparação
Corte do corpo Tratamentos Caracterização Ensaio de
das
de prova Térmicos das ligas Dureza
ligas
Foram elaboradas ligas de Cu-Al (liga 1), Cu- Tabela 1 é apresentada a porcentagem de
Al-Be (liga 2) e Cu-Al-Be-Ti-B (liga 3). Na elementos usados nas ligas.
Tabela 1 - Valores das Composições nominais das ligas fundidas.

Cu Al Be Ti B
Percentuais [%]
Liga 1 88,2 11,8 0 0 0
Liga 2 87,6 11,8 0,6 0 0
Liga 3 87 11,8 0,6 0,5 0,1
Após a determinação das porcentagens de 1200ºC e os elementos são colocados em um

elementos a serem utilizados os materiais cadinho de material inerte a base de argila-
foram pesados em uma balança analítica da grafite de forma sequencial, na qual os
marca Shimadzu modelo TX323L. A fundição primeiros materiais a serem fundidos são os
foi realizada em forno da marca Lucadema e com maiores pontos de fusão. O processo de
modelo Luca 200. O forno foi mantido a fundição pode ser visualizado na Figura 2.
Figura 2 - Processo de Fundição. (a) Elementos de liga sendo misturados. (b) Posicionamento do
cadinho no forno.
Após a fundição, as ligas foram vazadas em durante 12 horas em um forno da marca

moldes cerâmicos, com a solidificação do Lucadema e modelo Luca 200, com
material as mesmas foram usinadas com o subsequente tempera em água a 25ºC para
intuito de preparar corpos de provas (CPs) obtenção do efeito memória de forma.
para serem analisados através dos ensaios de
Os CPs foram preparados para metalografia
metalográficos e de dureza Brinell. Os CPs
com a sequência de lixamento de granulação
foram submetidos a dois tratamentos térmicos
80, 120, 220, 300, 600 e 1200. Polimento com
inicialmente foram homogeneizadas a 850ºC
alumina e finalizando com ataque químico
34
com reagente Cloreto de ferro III (5g FeCl3 , 30 3 RESULTADOS E DISCUSSÕES

ml de HCl fumegante 37% e 100 ml de H2 O
3.1 METALOGRAFIA.
destilada).
De acordo com Barbosa (2014) as ligas de
A análise foi realizada com auxilio do software
Cu-Al com 10% de alumínio, o aquecimento a
IPWIN32 na qual foram determinados os
uma temperatura superior a 850ºC resulta em
diâmetros dos furos que o penetrador fez
uma microestrutura constituída unicamente
sobre a superfície de cada amostra. A esfera
por fase β, e o resfriamento rápido leva à
do penetrador utilizada foi de 2,5 mm com
formação de uma microestrutura totalmente
uma carga de 187,5 N como determina a
martensítica. Por isso a presença de
norma ASTM (2014) para determinação da
plaquetas de agulhas no qual sugere a
dureza.
presença da fase martensítica à temperatura
ambiente nos corpos de prova temperados
observadas na imagem de micrografia
apresentadas na Figura 3.
Figura 3- Micrografia da liga Cu-Al da liga Cu-Al temperada com aumento de 100x.
O tamanho de grão foi determinado pela liga Cu-Al foi de aproximadamente 1,17mm. A
medição dos grãos nas imagens com imagem da micrografia com aumento de 5x
aumento de 5x. O tamanho médio do grão da pode ser visualizado na Figura 4.
Figura 4 – Micrografia da liga Cu-Al temperada com aumento de 5x
35
Por meio da análise da micrografia da liga Cu- mm. O Berílio atuou como refinador de grão,
Al-Be mostrada na Figura 5, identifica-se reduzindo o tamanho médio de grão em torno
apenas a fase austenítica CCC do tipo β2 com de 35%, em relação à liga Cu-Al.
tamanho de grão de aproximadamente 0,760
Figura 5 - Micrografia da Liga Cu-Al-Be temperada. Aumento de 5x.
Por meio da análise da Figura 6 percebe-se presença de aglomerados com pigmentação

que a liga Cu-Al-Be-Ti-B apresentou a fase preta, que estão circunscritos, tais
martensítica representada por (a), devido à aglomerados são precipitados de elementos
presença de agulhas, e a fase austenítica químicos que pode estar relacionado ao
representado por (b), devido à ausência de Titânio que devido à baixa solubilidade não
agulhas. Em certos pontos do grão ocorre a dissolve na matriz do grão.
Figura 6 - Micrografia da Liga Cu-Al-Be-Ti-B temperada. Aumento de 100x
O tamanho de grão da liga Cu-Al-Be-Ti foi redução percentual em relação as ligas Cu-Al
determinado através da Figura 7, na qual é de e Cu-Al-Be de aproximadamente 41% e 10%
aproximadamente 0,6813mm ocorrendo uma respectivamente.
36
Figura 7 - Micrografia da Liga liga Cu-Al-Be-Ti-B, com aumento de 5x.
3.2 DUREZA BRINELL. porcentagem de alumínio utilizado na liga foi

de 11%.
Na análise de dureza Brinell das Ligas CuAl,
Cu-Al-Be e Cu-Al-Be-Ti-B os valores obtidos Os valores de dureza encontrados para liga
estão apresentados na Tabela2. Cu-Al-Be são coerentes com a literatura, o
acréscimo de Be na liga Cu-Al gerou aumento
De acordo com Barbosa (2014), as ligas de
de aproximadamente 14,5% na dureza em
Cu-Al com 10% de alumínio, o aquecimento a
relação a liga Cu-Al. De acordo com Peng
uma temperatura superior a 850ºC e o
(1997), os valores de dureza com a adição de
resfriamento rápido torna essa liga com
Ti na liga Cu-Al deveriam melhorar as
dureza Brinell de cerca de 250 HB. Para
propriedades mecânicas do material, como
Sutou (2008), alumínio aumenta o
por exemplo, a dureza, porém a dureza para
ordenamento da fase β, promovendo maior
a liga 3 (Al-Cu-Be-TiB) encontrada foi
rigidez a liga. Com isso, os valores de dureza
aproximadamente 11,5% menor do que da
média encontrado de 268,66 HB nos ensaios
liga 1 (Cu-Al), e isso deu-se possivelmente
são coerentes com a literatura, visto que o
pela retirada de Ti da matriz por parte do
desvio padrão é de aproximadamente 6,2 e a
precipitado.
Tabela 2 - Valores de Dureza Brinell.
Desvio
Ligas Dureza [HB]
Padrão
Liga 1 268,70 6,2
Liga 2 307,33 7,5
Liga 3 237,88 5,2
4 CONCLUSÕES O acréscimo de elementos de liga, promoveu

uma redução do tamanho do grão e um
Neste trabalho ligas Cu-Al, Cu-Al-Be e Cu-Al-
aumento de dureza para a liga Cu-Al-Be,
Be-Ti-B foram fabricadas com sucesso pelo
processo de fundição sob atmosfera ambiente O acréscimo de elementos de liga promoveu
e caracterizadas por microscopia óptica e uma diminuição na dureza na liga CuAlBeTiB,
ensaios de dureza. que pode estar relacionado com precipitação
do TiB na matriz da liga.
A partir dos resultados obtidos é possível
concluir que a metalografia foi suficiente para
determinar as fases presentes em cada uma
AGRADECIMENTOS
das ligas, o tamanho dos grãos e identificar
precipitados. Foram realizados tratamentos Os autores agradecem a Universidade
térmicos Federal Rural do Semi-Árido (UFERSA) pelo
apoio a realização da pesquisa.
Os resultados de dureza encontrados estão
de acordo com a literatura mencionada no
decorrer do trabalho.
37
REFERÊNCIAS
[1] ALCÂNTARA, C. C. Influência do [15] FERNANDES, FRANCISCO MANOEL
crescimento anormal de grão nas propriedades BRAZ. Ligas com Memória de Forma, Universidade
mecânicas de ligas superelásticas cualmntib e Nova Lisboa CENIMAT. Abril de 2003.
cualmn. 2017. 68 f. Dissertação (Mestrado) - Curso
[16] FERNANDES JUNIOR, Paulo. Ensaios
de Engenharia Mecanica, Universidade Federal da
metalograficos. São Paulo: ., 2011.
Paraíba, João Pessoa - Paraíba, 2017.
[17] GARCIA, Amauri; SPIM, Jaime Alvares;
[2] ARAÚJO FILHO, O. O. Obtenção e
SANTOS, Carlos Alexandre dos. Ensaios dos
caracterização da liga quaternária β-cu-zn-al-ni
Materiais. 2. ed. Rio de Janeiro: Ltc Editora, 2012.
com efeito memória de forma. Dissertação de
mestrado. DEMEC/ UFPE, 2000. [18] GONZALEZ, C. H. Etude des
comportements electro-thermomécaniques et de la
[3] AMERICAN ASSOCIATION STATE
stabilisation martensitique d’alliages monocristallins
HIGHWAY. ASTM E10: Standard Test Method for
à mémoire de forme base cuivre. Tese D. Sc., INSA
Brinell Hardness of Metallic Materials. Eua: Astm,
de Lyon, France. 2002.
2014.
[19] GUNTHER, Eggeler, E. Hornbogen:
[4] AWAN, I. Z. fascinating shape memory
Werkstoffe mit Formgedächtnis. In: Das Magazin.
alloys. Institute Charles Gerhard, Ecole Nationale
Jg. 9, Nr. 1, 1998
Superieure de Chimie, Montpellier, France, 2018.
02p. [20] Kuhn G, Jordan L. Fatigue and mechanical
properties of nickel – titanium endodontic
[5] BARON, M. P. Etude du comportement
instruments. J Endod. 2002.
des interfaces austenite/martensite et
martensite/martensite de deux alliages a effet [21] LAGOUDAS, Dimitris C.. Shape Memory
memoire de forme: le CuCnAl et le CuAlBe. Tese Alloys: Modeling and Engineering Applications. Tx,
de doutoro, pp. 189. L’ Institut National Dês Usa: Springer, 2008.
Sciences Appliqees de Lyon, França. 1998.
[22] LANZINI, F., ROMERO, R., CASTRO, M. L.
[6] BARBOSA, C. Metais não Ferrosos e suas Influence of Be addition on order-disorder
Ligas: Microestrutura, Propriedades e Aplicações. transformations in β Cu-Al. Journal Intermetallics.
Rio de Janeiro: E-papers, 2014. 2008.
[7] BELKAHLA S. Elaboration et [23] LECCE, Leonardo; CONCILIO, Antonio.
caracterisation de nouveaux alliages à mémoire de Shape Memory Alloy ENGINEERING for Aerospace,
forme basse temperature type CuAlBe. Tese de Structural and Biomedical Applications. Naples:
doutorado, pp. 132, L’Institut National Dês Elsevier Ltd, 2015.
Sciences Appliqees de Lyon, França. 1990.
[24] LEXCELLENT, Christian. Shape-memory
[8] CAHN, R.W.; HAASEN, P. Physical Alloys Handbook. ISTE Ltd, Great Britain, 2013.
Metallurgy. North Holland (1983).
[25] MONTENCINOS, S., CUNIBERTI, A., R.,
[9] CÂMARA, Alexandre Sandri. Avaliação CASTRO, M. L. Kinetics of isothermal
das alterações geométricas e dimensionais dos decomposition in polycrystalline b CuAlBe alloys.
instrumentos de niti do sistema protaper para o Journal Intermetallics. 2009.
sistema protaper universal e o efeito destas sobre a
flexibilidade e a resistência torcional. Universidade [26] MONTECELLI, C. A. A competitividade da
Estadual Paulista. ARARAQUARA, 2008. 25p indústria brasileira de fundição. Dissertação –
Universidade Estadual de Campinas. Campinas,
[10] CZECHOWICZ, Alexander; LANGBEIN, 1994.
Sven. Shape Memory Alloy Valves. New York:
Springer, 2015. [27] NAVA, Marcelo. Estudo das Cinéticas de
Recristalização e Crescimento de Grãos de Ligas
[11] CALLISTER, Willian D.. Ciência e com efeito memória de forma. UNIVERSIDADE DE
Engenharia de Materiais. Utah: Ltc Editora, 2002. BRASÍLIA, Brasília. 2013.
[12] COLPAERT; Hubertus. Metalografia dos [28] NISHIYAMA, Z. Martensitic Transformation,
produtos siderúrgicos comuns, 3ª Edição,Editora Mats. Academoc Press, New York, 1978.
Edgarg Blücher Ltda, São Paulo – 1974.
[29] Otsuka K, Wayman CM. Shape memory
[13] FERREIRA, R.A.S. Transformação de fase. materials. Cambridge: Cambridge Univ. Press;
Aspectos cinéticos e morfológicos, Editora da 1998. 284p.
UFPE, 2002.
[30] Otsuka K, Ren X. Physical metallurgy of Ti-
[14] GUNTHER, E. Arco de E. Horn: Materiais Ni-based shape memory alloys. Progress in
com memória de forma . In: A revista . Vol. 9, No. 1 Materials Science. 2005.
, 1998.
38
Capítulo 4
Álvaro Barbosa Rocha

Felipe Silva Lima
Resumo: O crescente desenvolvimento tecnológico em aplicações de engenharia

tem demandado cada vez mais o emprego de materiais sofisticados com
propriedades excepcionais, as quais apenas os materiais convencionais não têm
sido capazes de suprir. Uma boa resposta a este cenário, são os compósitos,
devido à variedade de combinação que seus microconstituintes podem assumir,
derivando em produtos com economia de peso da ordem de 20% a 30%
comparado com os materiais tradicionais. Face a isto, este trabalho apresenta uma
análise numérica do comportamento mecânico de configurações arquitetadas em
materiais compósitos como solução ao arranjo original em madeira para importante
produto de desporto, o skate. Para a realização das análises foi utilizado o método
de elementos finitos com o auxílio do software comercial ANSYS®, especificamente
a extensão ACP (Ansys Composite Pre and Post) e ANSYS® Mechanical™
Academic Student Release 19.1. O método cumpriu as etapas de definição da
geometria inicial, geração da malha do modelo e aplicação de uma análise estática
para verificação do comportamento da estrutura frente as solicitações imputadas. O
modelo numérico atestou que as distintas associações de materiais influenciam
significativamente na seleção da constituição do compósito estrutural permeando
uma análise mais objetiva da estrutura erguida.
Palavras-chave: Análise numérica, Skate, Materiais compósitos
39
1 INTRODUÇÃO constituem, com a prerrogativa de suprir o

ávido mercado consumidor da categoria,
O desenvolvimento tecnológico
reservado a competidores e indivíduos, que
contemporâneo é cada vez mais significativo,
usufruem do esporte para lazer e diversão,
subsistindo a possibilidade de alçar produtos
totalizando cerca de 10 milhões de
para diversas áreas de aplicação, que
praticantes só nos EUA (MAIA,2015).
incorporam características dispares como,
por exemplo elevada resistência mecânica e Poucos estudos centralizaram como tema o
facilidade de fabricação adjunta a baixa skate e suas características mecânicas,
massa específica. Estrategicamente, em entretanto algumas pesquisas se sobressaem
qualquer das situações almeja-se lograr como Maia (2015), que desenvolveu e testou
produtos de qualidade e suficientemente inúmeros protótipos de pranchas (Shape)
confiáveis, fruto da combinação de novos fabricados em material compósito avaliando
materiais e tecnologias, sendo um diferencial parâmetros de rigidez, peso e som, a fim de
competitivo a aposta no desempenho estético encontrar uma configuração competitiva a um
e funcional, sem se desvencilhar de aspectos modelo comercial já existente. Já, Munshi e
relativos aos custos produtivos. Walame (2017) estimam o comportamento
mecânico de um shape concebida mediante o
O desporto avulta-se como um campo, em
emprego de compósitos laminados
que a Engenharia é levada ao limite de suas
reforçados por fibra de bambo. Outros
potencialidades, compelindo centros de
trabalhos enfocam, a tarefa de conceber
referência a desenvolver pesquisas, na
novos modelos de skate, como Andrade
expectativa que atletas arregimentem vitórias
(2015), que colabora com a criação de um
acompanhadas de marcas extraordinárias.
produto específico que atenda as demandas
Notadamente, o avanço destes objetivos
das modalidades donwhill e freeride.
passa pela tecnologia dos materiais,
despontando como uma opção de destaque No âmbito deste trabalho será efetuada uma
os compósitos, materiais cuja constituição análise numérica de diferentes disposições de
detém ao menos dois componentes um compósito sanduíche para prancha de
escolhidos para trabalhar em sinergia skate, contemple a Fig. 1, com faces
predispondo-se a obter uma peça com as constituídas por laminados reforçado por
propriedades requeridas para a aplicação fim fibras e um núcleo em material polimérico
(CHAWLA,1998). desenvolvido no software ANSYS®, que é
baseado no Método dos Elementos Finitos
A indústria de esportes radicais deposita
(MEF). Esta vertente é acolhida, de maneira a
grande expectativa no lançamento de
compreender o comportamento físico de
produtos funcionais que abarquem as
materiais compósitos, quando submetidos a
características supracitadas acima. Um item
esforços inerentes a aplicação de destino.
esportivo pertinente a esta discussão é o
Este estudo se faz necessário, uma vez que
skate, que desde o seu surgimento a partir do
na atualidade existe uma grande demanda da
surf, vem continuamente modificando não só
utilização de materiais compósitos em
suas formas, assim como os materiais que as
diversos setores da indústria.
Figura 1: Diferentes formatos de pranchas (Shape)
Fonte: Autoria própria (2018)
40
2 METODOLOGIA 2.1 MODELAGEM DO MATERIAL

EMPREGADO
O desenvolvimento do presente estudo
empregou o software comercial de elementos A prancha é idealizada como uma viga
finitos ANSYS Workbench 19.1 em sua versão sanduíche, desta forma optou-se por faces
Acadêmica, especificamente o módulo compostas por laminados cuja matriz incide
ANSYS Composite Pre/Post (ACP) e em um polímero do tipo epóxi reforçado com
Structural. O plano de consecução principiou fibras de carbono ou fibras de vidro, visto que
com a construção e especificação da garantem altas propriedades mecânicas de
geometria, material, malha de elementos resistência e rigidez, baixo peso, grande
finitos e, por fim as condições de contorno, durabilidade, somado a facilidade de
que segmentam em etapas os procedimentos fabricação. Por sua vez, o núcleo
mais importantes adotados nas simulações (enchimento) destas estruturas é formado por
numéricas. uma espuma polimérica de PVC entremeada
aos laminados citados, seu uso deve-se a sua
baixa densidade e compatibilidade química
com o material das faces. Acompanhe na
Tabela 1, as principais propriedades dos
constituintes.
Tabela 1: Propriedades mecânicas dos materiais em estudo

Propriedades Fibra de Vidro/Epóxi Fibra de Carbono/Epóxi Espuma de PVC
3
Densidade (g/cm ) 1,85 1,63 0,06
Módulo de Young (GPa) 39,5 215 0,01
Resistência a Tração (MPa) 965 1240 0,53
2.2 MODELAGEM DO ELEMENTO aplicações lineares, de grande rotação e/ ou

ESTRUTURAL de grande tensão não linear de estruturas não
maciças (elementos de casca) tendo as
O elemento utilizado para as simulações da
mesmas, estruturas finas ou moderadamente
prancha tem natureza bidimensional, em
grossas. O elemento possui oito nós com seis
razão da eficiência computacional e sua
graus de liberdade em cada nó: translações
melhor compatibilidade com o módulo de
nos eixos x, y e z e rotações em torno dos
compósitos do software ANSYS®. O elemento
eixos x, y e z, observe sua representação
bidimensional, SHELL281, apropriado para
abaixo.
Figura 2: Prancha elaborada em casca.
Fonte: ANSYS (2013).
41
A discretização do domínio baseia-se no interno e externo do sistema. A igualdade

principios dos deslocamentos virtuais, entre o trabalho interno e externo pode ser
tomando a igualdade entre o trabalho virtual expresso pela equação (1).
𝜹. 𝑾 = ⁡ ∫𝑽𝒆 𝜹𝒖𝒕 𝒃𝒅𝒗 + ⁡ ∫𝑽 𝒖𝒕 𝒎𝒖̈ 𝒅𝒗 (1)
Onde  corresponde ao trabalho virtual vetor de tensões reais para um elemento
externo,  ao vetor das deformações virtuais, genérico e com volume Ve e, são resolvidas
b ao vetor da força de volume real, u o vetor equações matriciais para a rigidez (K),
de deslocamento nodal da estrutura, m a amortecimento (C), massa (M) e forças nodais
matriz de densidade de massa e ü o vetor da (F), de modo que a equação geral do
aceleração do deslocamento nodal da movimento, equação (6) é satisfeita.
estrutura.
Em função da ocorrência de deformações ou
deflexões virtuais, provocadas por um dado
𝑲 = ⁡ ∑𝒏𝒆
𝒏=𝟏 𝑲𝒆 (2)
𝑴 = ⁡ ∑𝒏𝒆
𝒏=𝟏 𝑴𝒆 (3)
𝑪 =⁡∝. 𝑲 + ⁡𝜷. 𝑴 (4)
𝑭 = ⁡ ∑𝒏𝒆
𝒏=𝟏 𝑭𝒆 (5)
𝐊. 𝐮 + 𝐂. 𝐮̇ + 𝐌𝐮̈ = 𝐅 (6)
Onde ne é a matriz função de forma do 2.3 CONDIÇÕES DE CONTORNO

elemento ou matriz de interpolação dos
As condições de contorno empregadas para
deslocamentos correspondentes, Ke como a
a modelagem da problemática, aproximando
matriz de rigidez de cada elemento
o sistema real do sistema computacional,
empregado na discretização da geometria o,
foram formulados tomando a prancha com as
Me a matriz de massa de cada elemento
condições de fixação esboçada na Fig. 3. O
empregado na discretização da geometria, Fe
apoio simples (simply supported) introduz
a matriz de força correspondente a cada
vínculo na estrutura, impedindo o
elemento empregado na discretização da
deslocamento somente em uma direção,
geometria, 𝑢̇ como o vetor da velocidade de
produzindo reações equivalentes a uma força
deslocamento nodal da estrutura e  e 
com linha de ação conhecida, portanto
como coeficientes de ponderação.
remete a situação de uma viga em bi apoiada.
Figura 3: Ilustração da localização dos vínculos.
O valor de referência para o carregamento, região central da prancha, situação extrema

corresponde ao emprego de uma força de esforço para a estrutura em apreço, vide
dinâmica de impacto equivalente ao peso de Fig.4.
80 kg, equivalente a 787,8 N, localizado na
42
Figura 4: Ilustração da localização do carregamento sobre a estrutura.
2.4 MODELO GEOMÉTRICO popularmente conhecido por shape, que

corresponde a interface de contato direto com
Atualmente, os skates são formados por
o usuário, observe a Fig. 5. Quanto maior for o
diversos componentes. O enfoque,
shape mais estabilidade proporcionará em
especificamente, dado neste estudo é sobre o
altas velocidades, porém apresenta menor
elemento denominado de prancha ou mais
agilidade para curvas e rápidos desvios.
Figura 5: Prancha de Skate com seus respectivos elementos.
Na análise proposta assumisse, que a modelagem da estrutura é baseada na

prancha possui dimensões equivalentes a criação de uma estrutura em sanduiche,
81,3 cm (32”) de comprimento e 20,3 cm (8”) conforme figura 3, vindo a apresentar 3
de largura, e aproximadamente 1,3 cm (1/2”) camadas intercaladas com fibra de carbono
de espessura considerado o modelo mais ou de vidro e espuma PVC, com camadas de
comercializado pelos fabricantes de skates. A espessuras média de 0,4 cm.
43
Figura 4: Composição da prancha de Skate, ilustrando as camadas intercaladas com espessura

média de 0,4 cm.
A execução da simulação em ambiente deslocamento e não um dado gráfico,

ANSYS®, especialmente caso de materiais permitindo que a maior fração do poder
compósitos, exige modificações na geometria computacional seja direcionado a resolução
prévia, o que implica na transformação da das equações.
estrutura sólida 3D, para uma superfície As geometrias foram erguidas em um software
correspondente a uma casca de espessura CAD (Computer-Aided Design) de
virtualizada, implicando em uma superfície modelagem 3D, denominado SpaceClaim,
com o perfil da estrutura a ser simulada. O encarregado da criação, edição e porventura
processo de virtualização da espessura reparo da geometria, vide a Fig. 5. O emprego
permite que as propriedades das múltiplas desta plataforma sucedeu-se, em função da
camadas de compósitos seja um dado flexibilidade e facilidade nos comandos,
qualitativo inserido na matriz de equações aliado a integração com o software CAE
que avaliam o quadro de tensão e (Computer-Aided Engineering) ANSYS®.
Figura 5: Prancha elaborada em casca
Uma vez definida a geometria da estrutura, ACP (ANSYS Composite Pre and Post), a fim
para o desenvolvimento deste trabalho é de modelar o compósito, material participe do
utilizado um pré-processador denominado componente estrutural em julgamento.
44
2.5 GERAÇÃO DA MALHA (span angle center), foram configurados para

valores demasiadamente estreitos, a fim de
A malha retrata a partição da estrutura em
garantir a acurácia na convergência de
elementos, que descrevem aproximadamente
resultados. Ademais foi imputado a condição
o comportamento real do sistema em apreço.
de agressive mechanical, com perspectiva a
O ANSYS® permite ao usuário alterar
avaliar deformações apreciáveis e materiais
parâmetros de malha dependendo da
com comportamento intricado (não-linear). Ao
simulação em específico. O presente estudo
final, somado a esses procedimentos foi
transcorreu com a utilização de uma malha
estabelecido um controle de malha que
com relevância + 100, considerada o extremo
associou o face sizing e mapped face
valor – grosseira (-100) e refinada (+100) – em
meshing, almejando uma malha estruturada
escala ofertada pelo Software.
nas superfícies da peça.
Os comandos para dimensionamento (sizing)
O arranjo final contou com 6558 Nós e 6379
da malha, singularmente o centro de
elementos, observe a Fig. 6 que expõe o
relevância (relevance center), suavização
refinamento gerado da configuração inicial
(smoothing), centro do ângulo de expansão
para a configuração final.
Figura 6: Aspectos da malha criada
3 RESULTADOS E DISCUSSÕES somente os critérios de deformação devido a

formação de flecha na estrutura.
3.1 MODELAGEM NUMÉRICA POR
ELEMENTOS FINITOS
Por meio da modelagem da geometria 3.1.1 FIBRA DE CARBONO
empregando materiais distintos e, fazendo
A prancha reforçada por fibra de carbono
uso de um ambiente multífisico do ANSYS®,
embebida a uma matriz polimérica de epóxi,
foi possível determinar a partir da análise do
foi idealizada como uma viga sanduíche, no
material compósito, com número de camadas
qual entre as faces laminadas a um material
pré-determinadas, qual seria melhor
de enchimento a base de espuma polimérica.
empregado para a carga máxima teórica de
Ao presente arranjo constata-se uma
80 kgf. Para as análises apresentadas, a
deformação apresentada segundo as Fig. 7 e
seguir foram analisados e apresentados
8.
45
Figura 7 – Deformação sofrida pela prancha confeccioanda em fibra de carbono.
Ao confrontar a Fig. 7, percebe-se que o contrapartida na Fig. 9 a), constata-se a partir

maior deslocamento ocorre ao centro da do parâmetro denominado de Safety Margin
estrutura, local de aplicação do esforço, (margem de segurança) diretamente ligado
resultando em uma flexa máxima de ao fator de segurança da estrutura (Safety
1,2023mm. Por outro lado, os pontos onde se Factor), o fato que o corpo em apreço resiste
localizam os mancais de fixação das rodas a uma carga superior a duas vezes o valor
tem deformação nesta área de inicial, exprimindo na região em contraste
aproximadamente 0,3000 mm. exibida,em detalhes, na Fig. 9 b) uma
potencial área de colapso conforme os
A formação da flexa pela prancha pode ser
criterios de falha adotados (Tsai-Wu).
facilmente visualizada na Fig. 8. Em
Figura 8 – Formação de flecha sofrida pela geometria confeccionada em fibra de carbono.
46
Figura 9 – a) Esboço da margem de segurança por região; b) detalhe de um local de atenção para o
projetista.
a) b)
3.1.2 FIBRA DE VIDRO deformação é apresentada segundo as

Figuras 10 e 11.
Já a prancha em apreço foi modelada da
mesma maneira, entretanto o reforço passa a
ser fibra de vidro. Para o presente material a
Figura 10 – Deformação sofrida pela prancha confeccioanda em fibra de vidro.
Segundo a Fig. 10, o deslocamento máximo A formação da flecha pode ser vislumbrada
ocorre ao centro da estrutura originando uma na Figura 11. Em outra perspectiva na Fig. 12
flexa máxima para este caso de 1,5824 mm. A a), desprende-se que a estrutura reforçada
deformação na região próxima aos mancais em fibra de vidro apresenta um
de fixação das rodas beira 0,3515 mm. comportamento muito semelhante a reforçada
Diferentemente do caso anterior, a fibra de por fibra de carbono com uma margem de
vidro promove uma deformação maior, segurança um pouco menor, por conseguinte
justificavel pelas propriedades intrísecas do fica a cargo do executor do projeto a tarefa de
próprio reforço inferiores mecanicamente as ponderar o fatores como resistência mecânica
propriedades da fibra de carbono. e custo, a principio nesta situação
conflitantes.
47
Figura 11 – Formação de flecha sofrida pela geometria confeccionada em fibra de vidro.
Figura 12 – Formação de flecha sofrida pela geometria confeccionada em fibra de carbono.
a) b)
4 CONCLUSÕES opção o fator orçamentário necessitando a

posteriori uma análise de otimização, do qual
A partir da análise dos resultados obtidos,
conjuntamente são ponderados resistência e
pode-se concluir que as modificações no
custo em uma mesma função matemática.
comportamento estático repercutem uma
grande flexibilidade no que diz respeito as As zonas de fixação dos mancais das rodas
faixas de utilização do material compósito, expuseram deformações semelhantes para
entretanto, por uma ótica análoga, essa ambos os casos, com uma deformação de
característica pode se transformar em uma aproximadamente 0,3000 mm para a estrutura
grande desvantagem, sendo esta a em fibra de carbono e de 0,3515 mm para a
dificuldade na fabricação de materiais que estrutura em fibra de vidro, sendo que a
reproduzam o comportamento esperado. diferença entre as deformações são
justificavéis pelo comportamento intríseco do
Pode-se afirmar que dentro da análise
material usado como reforço e da presença
estática, tem-se um material mecanicamente
de furos, que atuam como concentradores de
superior, a disposição com fibra de carbono
tensão e promovem um quadro de
como material de reforço, exibiu uma deflexão
deformação considerável.
ao ser exposta a um carregamento transversal
inferior a cerca de 30 % a deflexão enfrentada
pela fibra de vidro, ou seja, promovendo o
AGRADECIMENTOS
direcionamento dos eixos adequadamente de
forma a se obter uma maior resistência do Os autores agradecem o apoio financeiro da
sistema segundo o critério de falha adotado CAPES.
(Tsai-Wu), entretanto sopesa sobre essa
48
REFERÊNCIAS
[1] ANSYS. ANSYS Composite PrepPost [5] MAIA, G.C., 2015, “Tábua de Skate em
User's Guide. Canonsburg: Sas Ip, Inc., 2013. resina poliéster reforçada com fibra de vidro”.
(Dissertação de Mestrado) – Instituto Superior de
[2] ANDRADE, P.N.B., 2015, “Chimpá
Engenharia do Porto – Departamento de
longboard”. (Trabalho de Conclusão de Curso) –
Engenharia Mecânica, Porto.
Universidade Estadual Paulista – Faculdade de
Arquitetura, Artes e Comunicação – Departamento [6] MALLICK, P.K. Composites Engineering
de Design, Bauru. Handbook. Marcel Dekker, Inc., NY, 1997.
[3] CHAWLA, K. K. Composite Materials, [7] MUNSHI, I.A; WALAME, M.V., 2017 “Finite
Science and Engineering. Springer Verlag, 1998. Element Analysis of Skate Board Made of Bamboo
Composite”. International Research Journal of
[4] LIU, H.; COOTE, T.; AIOLOS, C. 2017
Engineering and Technology, v.4, 2017.
“Skateboard deck materials selection” IOP
Conference Series: Earth and Environmental
Science, v.128, 2018.
49
Capítulo 5
Carlos Cássio de Alcântara

Jéssica Menezes Alves
José Joelson de Melo Santiago
Resumo: Neste trabalho foi realizado a caracterização mecânica e microestrutural

de uma liga do sistema Cu-Al-Mn passível de efeito memória de forma. As ligas
foram elaboradas em forno sob atmosfera ambiente, tratadas termicamente por
homogeneização e têmpera, e analisadas por ensaios de microscopia óptica,
difração de raios-x, microdureza e tração uniaxial. Os resultados obtidos dos
ensaios de DRX mostraram que a liga bruta de fundição a temperatura ambiente
(Tamb ≈25 ºC) apresentou a fase β, fase γ e Al8Mn5. A liga tratada termicamente
apresentou a temperatura ambiente a fase Austenítica β e após os ensaios de
tração exibiu fases martensiticas γ’ de fórmula estrutural AlCu3, fase β e fase γ. O
ensaio de microdureza determinou a dureza das fases presentes contribuindo para
a classificação das fases encontradas. O ensaio de tração mostrou um
comportamento linear associado à deformação da fase austenítica presente na liga
à temperatura ambiente.
Palavras-chave: Ligas memória de forma. Liga Cu-Al-Mn. Efeito memória de forma.
50
1 INTRODUÇÃO adequadas podem potencializar e aumentar

suas aplicações em diversas áreas. Visando
As ligas com memória de forma (LMF) são
modificar as propriedades térmicas e
metais que, após sofrerem uma deformação,
mecânicas de ligas à base de cobre, é
têm a capacidade de retornar à sua forma
realizada a inserção de um terceiro ou quarto
original quando aquecidas a determinada
elemento na composição deste sistema,
temperatura. Tais ligas possuem
modificando assim: temperaturas de
propriedades características, tais como:
transformação, grau de deformação,
transformação martensítica, efeito memória de
ductilidade, tamanho de grão entre outras
forma (EMF) e superelasticidade (SE)
propriedades metalúrgicas. As vantagens das
(CALLISTER, 2008). A transformação
ligas ternárias Cu-Al-Mn, comparadas com
martensítica é uma transformação de fase no
metais comuns e outras ligas memória de
estado sólido que ocorre sem difusão de
forma, são a boa resistência mecânica e alta
átomos. Essa transformação é responsável
capacidade de amortecimento. As ligas Cu-
pela mudança de forma entre a fase de alta
Al-Mn têm baixo custo de produção em
temperatura denominada austenita (A), e a
relação a ligas NiTi, devido a não
fase de baixa temperatura chamada
necessidade de fundição a vácuo. As
martensita (M) (CANBAY; GENC; SEKERCI,
temperaturas de transformação podem ser
2014). A transformação de fase é reversível e
ajustadas de -180ºC até 230ºC, através da
não dependente do tempo, uma vez que só
composição química, adição de elementos
distorce a rede atômica sem causar migração
quaternários, processos de fabricação ou
de átomos de longo alcance. Outra forma de
pelos tratamentos térmicos subsequentes.
obter a transformação martensítca é através
Sabe-se que para baixos teores de Al, as
da aplicação de tensão. Essa tensão
ligas exibem boa ductilidade e capacidade de
necessária para promover a transformação de
trabalho a frio. Adicionalmente, o Mn promove
fase é chamada de Tensão de Indução
o aumento da fluidez da liga no processo de
Martensítica (CANBAY; GENC; SEKERCI,
fundição (MIELCZAREK et al., 2008). Neste
2014).
contexto o objetivo deste trabalho é elaborar
As LMF têm sido extensivamente estudadas, uma liga Cu-Al-Mn e avaliar as características
nos últimos anos, para aplicação em diversas microestruturais através de ensaios de
áreas (Rao, et al., 2015). Pode-se destacar microscopia óptica, difração de raio X,
aqui as ligas, Ni-Ti, Cu-Al-Ni, Cu-Al-Be e Cu- microdureza e tração uniaxial.
Al-Mn, como as ligas mais utilizadas. Dentre
estas as ligas Ni-Ti, são as que possuem
melhores propriedades termoelásticas e que 2 MATERIAIS E MÉTODOS
em particular, são amplamente utilizadas em
As atividades desenvolvidas foram
muitas aplicações de engenharia e
inteiramente realizadas nas instalações da
biomédicas. Algumas das desvantagens que
Universidade Federal Rural do Semi-Árido
essas ligas sofrem são baixas temperaturas
(UFERSA). A metodologia adotada está
de transformação, dificuldade na produção e
apresentada no fluxograma da Figura 1.
processamento e alto custo. Por outro lado, as
ligas com memória de forma baseadas em Neste trabalho foi elaborada a liga
cobre são mais fáceis para produzir, policristalina Cu-Al-Mn a partir de elementos
processar e são menos onerosas (CALUÊTE, de alta pureza (99,99%) na composição Cu-
2012). Portanto desenvolver ligas memória de 17%Al-11,2%Mn (% at.)
forma do sistema Cu-Al com propriedades
51
Figura 1- Fluxograma da sequência metodológica
ELABORAÇÃO TRATAMENTOS CARACTERIZAÇÃO PROPRIEDADES

DA LIGA TERMOMECÂNICOS MICROESTRUTURAL MECÂNICAS
Confecção dos Homogeneizaç

Metalografia Microdureza
CPs ão
Bruta de Fusão Têmpera DRX Tração
Temperada
Tracionada
Os elementos foram previamente pesados em As amostras do grupo Tracionada além da

uma balança analítica da marca Shimadzu têmpera e homogeneização foram
modelo TX323L. O manganês, antes da tracionadas até a ruptura.
pesagem, foi submetido ao processo de
O tratamento de homogeneização foi
decapagem para retirada da camada
realizado a 850ºC por um período de 12 horas
superficial de óxido. Para isso, uma certa
com resfriamento lento dentro do próprio
quantidade do elemento foi imersa em uma
forno. A têmpera foi realizada com
solução aquosa com ácido nítrico 10%, de
aquecimento a 850ºC, mantido por 30min,
composição H2O-5%HNO3 (% em volume). A
seguido de resfriamento em água a
liga foi fundida em forno tipo mufla da marca
temperatura ambiente (Tamb ≈25 ºC).
Lucadema, modelo LUCA 2000, em cadinho
de argila-grafite, sob atmosfera ambiente e A análise microestrutural foi realizada com o
vazada em um molde de seção retangular auxílio de um microscópio óptico Axiotech 30
com dimensões 120mm de comprimento x da Carl Zeiss. As amostras foram submetidas
22mm de largura x 40mm de altura. O lingote a lixamento manual com lixas de numeração
bruto de fundição foi homogeneizado a 850ºC granulométrica 80, 100, 320, 600 e 1200,
por um período de 12 horas com resfriamento seguidas de polimento mecânico em alumina
lento dentro do próprio forno. Após a 1µm. Para a visualização dos contornos de
homogeneizado, foram confeccionados grão foi utilizado o método de ataque por
corpos de prova (CPs) adaptados. As esfregação, onde uma porção de algodão foi
dimensões dos CPs estão dispostas na Figura imersa em um reativo de cloreto de ferro III
2. No processo de confecção foi utilizada uma (5g FeCl3, 30 ml de HCl 37% e 100ml de H2O
máquina de eletroerosão a fio, modelo destilada) e esfregada contra a superfície das
Electrocut AR 1300. amostras. O tamanho de grão das amostras
foi realizado através da determinação do
Os corpos de prova foram divididos em 3
comprimento médio do intercepto (DCMI) na
grupos: Bruto de Fundição (BF), Temperado e
direção longitudinal pelo processo de Heyn,
Tracionado. As amostras do grupo BF foram
descrito pela Asociación Mercosur de
apenas homogeneizadas. As amostras do
Normalización (2001).
grupo Temperada foram submetidas a
tratamentos de homogeneização e têmpera.
52
Figura 2 - Dimensões dos CPs.
A identificação das fases foi realizada através 10s. Os ensaios foram realizados à
de difração de raios-x. O ensaio foi realizado temperatura de 20ºC.
com o auxílio de um difratômetro da marca
As propriedades mecânicas da liga foram
Shimadzu, modelo Lab X/XRD-6000, passo de
determinadas através de um ensaio de tração
0,02°, varredura de 10 a 90º e tempo de
uniaxial, com auxílio de uma máquina da
passo de 1s, utilizando a radiação K-alfa do
AUTOGRAPH modelo AG-X 10kN, nos CPs
cobre como fonte de radiação
Temperado. A velocidade de aplicação da
monocromática. Os dados gerados no ensaio
carga no ensaio foi de 0,5 mm/min, até a
foram interpretados através do software X’Pert
ruptura. Desse ensaio foram obtidos a tensão
HighScore Plus, sendo possível a coleta do 2θ
de indução martensítica através do método
e do espaçamento interplanar (di) referentes
das tangentes, o limite de resistência à tração
aos picos de difração. Essas informações
e a deformação máxima.
coletadas foram comparadas com o banco de
dados de fases e estruturas do aplicativo
Crystmet, permitindo a identificação das fases
presentes nos CPs. Todas as análises foram
realizadas em temperatura ambiente. 3.1 METALOGRAFIA.
O ensaio de microdureza foi realizado com A micrografia das amostras BF pode ser
auxílio de um microdurômetro da marca visualizada na Figura 3 que revela a presença
Shimadzu modelo HMV-2 Series, com carga de três regiões distintas: a região clara, a
0.100 kgf e tempo de aplicação da carga de região escura ao fundo e os precipitados sob
a forma de esferas no interior dos grãos.
Figura 3 - Micrografia das amostras BF: a) Micrografia com aumento 5x. b) Micrografia com aumento
de 20x.
53
A micrografia das amostras do grupo plaquetas ou agulhas características da fase

Temperada, ilustrada na Figura 4, revela a martensítica. O tamanho de grão mensurado
presença da fase austenítica β. Um fator que nas amostras foi de 662 µm.
corrobora com essa ideia é a ausência de
Figura 4 - Micrografia das amostras tratadas termicamente
Nas micrografias das amostras do grupo transformação martensítica induzida por

Tracionada, Figura 5, pode-se observar a tensão
presença de agulhas provenientes da
Figura 5 - Micrografia dos CPs temperados após ensaios de tração: a) Micrografia com aumento 5x.
b) Micrografia com aumento de 10x
As micrografias obtidas através do amostras do grupo Temperada podem ser

resfriamento em nitrogênio líquido das visualizadas na Figura 6.
54
Figura 6 - Micrografia com Nitrogênio Líquido dos Corpos de prova tratados termicamente: a)
Micrografia no início da transformação de fase. b) Micrografia no fim da transformação de fase.
3.2 DIFRAÇÃO DE RAIOS-X a têmpera. A amostra exibe apenas a fase

austenítica β, de estrutura L21 e grupo
A Figura 7 apresenta os resultados obtidos do
espacial Fm3m, essa com plano de difração
ensaio de DRX realizado nas amostras BF,
(4 2 2). A amostra Tracionada exibe a fase β
Temperada e Tracionada.
de plano (2 2 0), a fase γ de plano (3 3 0) e a
Na amostra BF é possível identificar a fase β, fase martensítica ortorrômbica γ’, de fórmula
de estrutura L21 e grupo espacial Fm3m, cujo estrutural AlCu3 e grupo espacial P2mm, sob
plano de difração é (2 0 0). Os planos (3 0 0) os planos (5 1 1) e (1 2 10), que representa
e (3 3 2) são referentes à fase γ, que uma das duas fases martensíticas que podem
apresenta estrutura Cu9Al4 e grupo espacial ser formadas em ligas Cu-Al-Mn. Sua
P-43m. Uma terceira fase, cujo plano de sequência de empilhamento, segundo a
difração é o (4 4 4), é identificada sob a forma notação de Ramsdell, é 2H (OBRADÓ et al.,
de Al8Mn5, de estrutura CCC e grupo 1998, p. 14245).
espacial I-43m. Na amostra Temperada é
possível notar a supressão das fases
existentes na amostra bruta de fundição após
Figura 7 - Difratograma das amostras.

Tracionada
Temperada
16000
Bruta de Fusão
14000 '
 
12000
10000

Intensidade
8000
6000
4000
 Al8Mn5
2000
 
10 20 30 40 50 60 70 80 90
2
55
3.3 MICRODUREZA do ensaio de DRX para a classificação das

fases. A Tabela 1 mostra o valor de dureza
Ensaios de microdureza foram realizados para
dessas fases.
a obtenção da dureza das regiões
contrastadas na micrografia. Esses resultados
serão utilizados juntamente com os resultados
Tabela 1 - Dureza das fases presentes na amostra bruta de fusão

FASE CARGA (GF) DUREZA (HV)
HV0.100 395
HV0.100 419
HV0.100 380
Escura Média 398
Desvio 19,67
HV0.100 110
HV0.100 115
HV0.100 98,1
Clara Média 107,70
Desvio 8,68
3.4 ENSAIO DE TRAÇÃO pseudoplásticas de ordenamento das

diversas variantes da martensita. Logo após,
Na Figura 8 está mostrada uma curva obtida
um novo comportamento linear, agora, com
do ensaio de tração. Nesta curva, verifica-se
menor inclinação, é desenvolvido. Este se
em um primeiro estágio um comportamento
associa ao menor módulo de elasticidade que
linear associado à deformação da fase
a variante martensítica predominante
austenítica presente na liga à temperatura
apresenta em relação à sua respectiva fase
ambiente. Em seguida, observa-se um
austenítica. Por fim, termina-se com a queda
decaimento na inclinação da curva. Esse
repentina do valor de tensão, caracterizando
decaimento está relacionado às deformações
a ruptura do material.
Figura 8 - Gráfico Tensão x Deformação da Liga Cu71,8Al17,0Mn11,2
Na Tabela 2 estão os valores das formação está próxima da temperatura na

propriedades obtidos no ensaio de tração. Os qual foi realizado o ensaio de tração. Os
baixos valores da Tensão de Indução dados referentes ao Limite de Resistência à
sugerem que a temperatura de trans-
56
Tração da liga apresentaram pouco desvio e alcançaram uma média de 253,841 MPa.
Tabela 2 - Gráfico Tensão x Deformação da Liga Cu71,8Al17,0Mn11,2

CORPO DE TENSÃO DE INDUÇÃO LIMITE DE RESISTÊNCIA À DEFORMAÇÃO
PROVA MARTENSÍTICA (MPA) TRAÇÃO (MPA) MÁXIMA (%)
1 113 260,9 18,6
2 127 210,3 7,6
3 100 208,0 8,3
4 118 252,7 16,0
5 120 269,0 13,9
6 123 292,5 14,5
7 69 266,0 11,8
8 115 248,9 12,4
9 50 276,2 15,2
Média 103,9 253,8 13,2
Desvio Padrão 20,6 28,4 3,6
4 CONCLUSÕES Do ensaio de tração uniaxial foram obtidos a

Tensão de Indução Martensítica, o Limite de
Os CPs temperados apresentaram tamanho
Resistência à Tração e a Deformação Máxima
de grão médio de 662 µm.
dos CPs Temperados.
O processo de têmpera promoveu a
supressão das fases γ e Al8Mn5, e exibe
apenas a fase β. AGRADECIMENTOS
Houve a formação da fase martensítica γ’ nos Os autores agradecem a Universidade
CPs tracionados. Federal Rural do Semi-Árido (UFERSA) pelo
apoio a realização da pesquisa.
Através do ensaio de microdureza foi possível
determinar as durezas das fases na amostra
Bruta de Fusão e na Temperada.
REFERÊNCIAS
[1] ALCÂNTARA, C. C. de. INFLUÊNCIA DO [5] CALLISTER, W. D. Ciência e engenharia
CRESCIMENTO ANORMAL DE GRÃO NAS de materiais: uma introdução. 7. ed. LTC, 2008.
PROPRIEDADES MECÂNICAS DE LIGAS ISBN 978-85-216-1595-8.
SUPERELÁSTICAS CuAlMnTiB e CuAlMn. 2017. 67
[6] CANBAY, C. A.; GENC, Z. K.; SEKERCI, M.
p. Dissertação (Mestrado em Engenharia
Thermal and structural characterization of Cu–Al–
Mecânica) — Universidade Federal da Paraíba.
Mn–X (Ti, Ni) shape memory alloys. Applied
[2] ASOCIACIÓN MERCOSUR DE Physics: Materials Science & Processing, A, n. 115,
NORMALIZACIÓN. NM241:2001 - Determinação do p. 371 – 377, Apr 2014.
tamanho de grão em materiais metálicos. 2001.
[7] CZECHOWICZ, A.; LANGBEIN, S. Shape
[3] AZEVEDO, M. R. D. ESTUDO DO SISTEMA Memory Alloy Valves: Basics, Potentials, Design.
DE LIGAS Cu-Al PASSÍVEIS DO EFEITO MEMÓRIA Springer, 2015. ISBN 978-3-319-19081-5.
DE FORMA. 2011. Monografia (Bacharel em
[8] LECCE, L.; CONCILIO, A. Shape Memory
Ciência e Tecnologia) — Universidade Federal
Alloy ENGINEERING for Aerospace, Structural and
Rural do Semi-Árido - UFERSA.
Biomedical Applications. : Elsevier, 2015. ISBN
[4] BUEHLER, W. J.; GILFRICH, J. W.; 978-0-08-099920-3.
WILLEY, R. C. Effect of low temperature phase
[9] MALLIK, U.; SAMPATH, V. Influence of
changes on the mechanical properties of alloys
quaternary alloying additions on transformation
near composition TiNi. J. Appi. Phys, n. 34, p. None
temperatures and shape memory properties of Cu-
– 1475, 1963.
Al-Mn shape memory alloy. Journal of Alloys and
57
Compounds, v. 469, p. 156 – 163, Feb 2009. ISSN MEMÓRIA DE FORMA. 2006. 69 p. Dissertação
0925-8388. (Mestrado em Engenharia Mecânica) —
Universidade Federal de Pernambuco.
[10] MIELCZAREK, A. et al. Mechanical and
fatigue properties of Cu - Al - Mn shape memory [14] SUTOU, Y. et al. Ductile Cu–Al–Mn based
alloys with influence of mechanical cycling on shape memory alloys: general properties and
amplitude dependence of internal friction at room applications. Materials Science and Technology, v.
temperature. Solid State Phenomena, v. 137, p. 145 24, n. 8, p. 896 – 901, 2008.
– 154, Mar 2008.
[15] Rao, A., Sirivasa, A. R., Reddy, J.N., 2015.
[11] OBRADÓ, E. et al. Order-disorder Design of Shape Memory Alloy (SMA) Actuators.
transitions of Cu-Al-Mn shape-memory alloys. Springer Internationa Publishing, Texas, 1ft edition.
PHYSICAL REVIEW B, v. 58, n. 21, p. 14 245 – 14
[16] ZAK, G.; KNEISSL, A.; ZATULSKIJ, G.
255, 1998.
SHAPE MEMORY EFFECT IN CRYOGENIC Cu-Al-
[12] OTSUKA, K.; WAYMAN, C. M. Shape Mn ALLOYS. Scripta Materialia, v. 34, n. 3, p. 363 –
memory materials. University Press, 1998. 367, 1996. Software CRYSMET
[13] PINA, E. A. C. de. ESTUDO DA
ESTABILIZAÇÃO NA LIGA Cu-Al-Mn COM EFEITO
58
Capítulo 6
João Manoel de Oliveira Neto

Marcelo Bezerra Grilo
Resumo: O objetivo deste estudo é verificar experimentalmente a influência da

dilatação térmica sobre a gasolina automotiva e o etanol hidratado combustível,
comercializados no município de Campina Grande (PB). O universo pesquisado
consistiu em coletar as amostras de combustível em vinte e quatro (24) postos de
comercialização, que foram visitados nos períodos noturno e diurno. Para avaliar os
efeitos da dilatação térmica, no momento da coleta, em janeiro de 2016, foi
realizado o teste do densímetro e aferida a temperatura das amostras coletadas. Os
resultados obtidos foram, com uma variação média relativa entre os períodos diurno
e noturno, de 0,004 g/ml para o etanol e 0,005 g/ml para a gasolina. Estes
resultados mostram que a temperatura tem influência sobre os volumes
comercializados, especialmente para a gasolina, por possuir um maior coeficiente
de dilatação térmica. Este problema pode ser resolvido quando o combustível for
comercializado em massa (kg) ao invés de volume (litro).
Palavras-Chave: Dilatação térmica, etanol hidratado, gasolina automotiva,

biocombustíveis e temperatura.
59
1 INTRODUÇÃO maior temperatura respectivamente. Para o

período noturno, o horário escolhido foi após
A diferença de temperatura entre o período
as 20:00 h, devido a maioria dos postos de
mais quente e mais frio, ao longo do dia, faz
comercialização fecharem as 22:00 h. E pela
com que o consumidor final acabe adquirindo
manhã, boa parte desses abre após as 7:00
quantidades diferentes de massa de
h, o que dificulta os experimentos.
combustível. [1]. Isto acontece devido ao
fenômeno da dilatação térmica dos líquidos. Na sequência, as amostras coletadas foram
Fenômeno esse mais notável no munícipio de acondicionadas em vidros âmbar (cor
Campina Grande-PB, caracterizado por dias marrom) de um litro, tamponados com
quentes e noites frias, ou seja, com grande batoques, e etiquetadas com o código
amplitude térmica diária. É esperado que numérico para identificação da amostra
devido esta influência de temperatura ao referente ao posto, de acordo com a relação
longo do dia, o consumidor final acabe previamente registrada no Laboratório
adquirindo uma menor massa de combustível, Experimental de Máquinas Térmicas (LEMT)
seja gasolina ou etanol, no período diurno da Universidade Federal de Campina Grande
(período em que são apresentados maiores (UFCG). Sendo identificados por códigos alfa
valores de temperatura). Uma análise numéricos de PG1, PE1, PG2, PE2, ....., PG24.
experimental, no que diz respeito a aferição
Logo em seguida, após a coleta de cada
da massa dos combustíveis nos postos de
amostra, o combustível foi despejado em um
comercialização (simulação da atividade de
recipiente graduado em 500 ml. Foi aferida a
abastecimento), pode ser comprovado se
temperatura do combustível com o auxílio de
realmente essa variação é significante.
um termômetro de mercúrio e realizado o
Motivado pelo crescimento da frota teste do densímetro em uma superfície plana
automotiva nas cidades e rodovias brasileiras por mais ou menos dois minutos, ainda no
e pelo respectivo aumento no preço e no posto de distribuição visando simular o
consumo do etanol e da gasolina, faz-se abastecimento e evitar a perda de energia
necessário este estudo, referente a influência dos combustíveis, onde foi obtida a leitura da
da variação de temperatura sob estes massa específica, similar ao que é exigido
combustíveis, de modo que a comercialização pelas NBRs 5992 e 7148. Após realizadas as
ocorra da forma que o consumidor final venha aferições de massa específica e temperatura,
a ser o mais beneficiado. foi verificado nas NBRs 7148 e 5992 se as
amostras atendiam as requisitos exigidos
mediante a qualidade.
2 PROCEDIMENTO EXPERIMENTAL
Foram coletadas vinte e quatro amostras de
cada combustível, no que diz respeito a
gasolina automotiva comum tipo C e etanol Foi constatado de acordo com a Tabela 1 que
hidratado no volume padrão de um litro. Para para a gasolina houve uma variação média de
seguir os procedimentos recomendados pelas temperatura e densidade de 5,04 °C e 0,005
NBRs 5992 e 7148, as coletas para o período g/ml entre os períodos diurno e noturno e as
noturno ocorreram no horário após 20:00 h, e amostras encontravam-se em conformidade
para o período diurno, o horário de coleta foi no que diz respeito a qualidade referente ao
entre as 14:00 e 16:00 h, períodos de menor e teste do densímetro.
60
Tabela 1: Variação da densidade (g/ml) da gasolina tipo C comum em função da variação da

temperatura, com as coletas simulando o abastecimento noturno e diurno, nas regiões central e
periférica – janeiro 2016.
Massa
T (ºC) específica T (ºC) Massa
Posto (noturno) específica Δρ (g/ml)
(Noturno) (Diurno) (diurno) (g/ml)
(g/ml)
PG1 26.9 0,740 32,2 0,735 0,005
PG2 26,0 0,741 31,2 0,736 0,005
PG3 26,2 0,741 31,6 0,736 0,005
PG4 26,1 0,741 31,3 0,736 0,005
PG5 28,3 0,739 33,4 0,734 0,005
PG6 27,4 0,740 32,1 0,736 0,004
PG7 27,6 0,740 32,8 0,735 0,005
PG8 26,0 0,741 30,7 0,736 0,005
PG9 26,8 0,742 31,6 0,738 0,004
PG10 27,1 0,741 31,1 0,737 0,004
PG11 26,2 0,741 31,0 0,736 0,004
PG12 26,3 0,745 30,4 0,741 0,004
PG13 26,0 0,741 31,4 0,736 0,005
PG14 26,3 0,741 31,8 0,736 0,005
PG15 26,4 0,740 31,6 0,735 0,005
PG16 27,0 0,741 32,6 0,736 0,005
PG17 26,6 0,742 31,9 0,737 0,005
PG18 26,3 0,741 31,1 0,737 0,004
PG19 26,2 0,741 31,4 0,735 0,005
PG20 26,1 0,741 31,7 0,736 0,005
PG21 26,4 0,742 32,0 0,736 0,005
PG22 26,9 0,740 31,8 0,735 0,005
PG23 27,1 0,740 30,2 0,737 0,003
PG24 27,0 0,741 33,2 0,736 0,005
Para o etanol, de acordo com a Tabela 2, encontravam-se em conformidade no que diz

houve uma variação média de temperatura e respeito a qualidade referente ao teste do
densidade de 4,55 °C e 0,004 g/ml entre os densímetro.
períodos diurno e noturno e as amostras
61
Tabela 2: Variação da densidade (g/ml) do etanol comum em função da variação da temperatura,

com as coletas simulando o abastecimento noturno e diurno, nas regiões central e periférica –
janeiro 2016.
T (ºC) Densidade T (ºC)
Posto Densidade (diurno) Δρ (g/ml)
(Noturno) (noturno) (Diurno)
PE1 26,1 0,802 30,9 0,798 0,004
PE2 26,2 0,802 30,4 0,798 0,004
PE3 25,7 0,804 32,2 0,798 0,006
PE4 26,0 0,802 30,3 0,798 0,004
PE5 26,4 0,802 31,8 0,798 0,004
PE6 26,3 0,802 30,2 0,799 0,003
PE7 26,2 0,802 32,1 0,798 0,004
PE8 27,3 0,800 30,3 0,798 0,002
PE9 27,1 0,802 31,2 0,798 0,004
PE10 27,1 0,802 31,3 0,798 0,004
PE11 26,6 0,802 30,1 0,798 0,004
PE12 27,2 0,802 31,3 0,798 0,004
PE13 27,0 0,800 31,5 0,796 0,004
PE14 26,3 0,802 30,2 0,798 0,004
PE15 27,8 0,802 32,4 0,798 0,004
PE16 26,4 0,802 31,1 0,798 0,004
PE17 27,2 0,800 30,2 0,798 0,002
PE18 27,8 0,802 32,7 0,798 0,004
PE19 27,1 0,802 31,4 0,798 0,004
PE20 26,4 0,800 30,9 0,796 0,004
PE21 27,1 0,802 32,2 0,798 0,005
PE22 26,5 0,802 31,0 0,798 0,004
PE23 27,7 0,802 32,8 0,798 0,004
PE24 26,2 0,802 31,9 0,796 0,006
4 CONCLUSÕES levando uma maior quantidade em massa (kg)

para um volume fixo de combustível, quando
De acordo com os dados obtidos quanto a
realiza o abastecimento no período noturno.
variação da massa específica da gasolina e
Essa variação é significativa quando é
do etanol, foi obtido uma média de 0,005 e
comprado combustível em grandes
0,004 g/ml para uma variação média de
quantidades. Para uma logística mais
temperatura de 5,04 e 4,55°C, entre os
eficiente, gerando maior satisfação ao
períodos diurno e noturno. Esses dados
consumidor final e máxima eficiência
comprovam que o efeito da dilatação térmica
energética dos combustíveis, os fatores
sobre ambos os combustíveis, havendo uma
explanados no presente trabalho devem ser
certa sensibilidade decorrente de um
levados em consideração, ou seja, o
gradiente de temperatura, que influencia
combustível deveria ser vendido em massa, e
diretamente no aumento do seu volume.
não em volume. Enquanto forem vendidos em
No momento da compra dos combustíveis, os volume, preferencialmente realizar o
mesmos são vendidos em volume (litros). abastecimento no período noturno, onde os
Como a massa varia em função da efeitos da dilatação térmica são menores.
temperatura, o consumidor final, acaba
REFERÊNCIAS
[1] ALMEIDA SILVA, Felipe Augusto. Estudo [2] ASSOCIAÇÃO BRASILEIRA DE NORMAS
sobre a Qualidade da Gasolina Comercializada em TÉCNICAS. NBR 5992: Álcool etílico e suas
Campina Grande e da Variação de sua Massa misturas com água - Determinação do teor de
Específica em função do Horário de álcool etílico anidro combustível (AEAC) – Método
Abastecimento. Trabalho de Conclusão de Curso do densímetro, Rio de Janeiro, 2008.
(Monografia), Unidade Acadêmica de Engenharia
[3] ASSOCIAÇÃO BRASILEIRA DE NORMAS
Mecânica, Universidade Federal de Campina
TÉCNICAS. NBR 7148: Petróleo e produtos de
Grande, Campina Grande, PB. 2014.
petróleo: Determinação da massa específica,
62
densidade relativa e °API - Método do densímetro. [6] GONÇALVES, Dalton. Física: Termologia,
Rio de Janeiro, 2013. ótica, ondas. Vol. 1, 3ª edição. Capítulo 2: dilatação
térmica, São Paulo: editora moderna, 1979.
[4] CASULO DE ARAÚJO, Helder Giuseppe.
Estudo experimental sobre a qualidade da gasolina [7] OLIVEIRA NETO, J. M., EMERENCIANO, L.
comercializada no estado da Paraíba com o uso da F., GRILO, M. B., Análise experimental da gasolina
norma ABNT 13992:2008 e do planejamento comercializada em Campina Grande, incluindo a
fatorial. 2013. 91 p. Dissertação de Mestrado. influência do horário de abastecimento e a
Programa de Pós-Graduação em Engenharia localização do posto de combustível. In: I
Mecânica (PPGEM), Universidade Federal de CONGRESSO NACIONAL DE ENGENHARIA DE
Campina Grande, Campina Grande. 2013. PETRÓLEO, GÁS NATURAL E BIOCOMBUSTÍVEIS,
número 1. 2015, Campina Grande, PB. Anais...
[5] DOMINGOS, I. S. et al., Determinação do
Campina Grande: Editora Realize, 2015. 8p.
aspecto, massa específica e teor de etanol de
amostras de gasolina tipo C comum, [8] OLIVEIRA NETO, J. M., GRILO, M. B.,
comercializadas em Patu-RN. In: I CONGRESSO Análise experimental da variação de massa da
NACIONAL DE ENGENHARIA DE PETRÓLEO, GÁS gasolina comercializada na cidade de Campina
NATURAL E BIOCOMBUSTÍVEIS, número 1. 2015, Grande-PB devido os efeitos da dilatação térmica.
Campina Grande, PB. Anais... Campina Grande: In: II CONGRESSO NACIONAL DE ENGENHARIA
Editora Realize, 2015. 7p. DE PETRÓLEO, GÁS NATURAL E
BIOCOMBUSTÍVEIS, número 2. 2016, Natal, RN.
Anais... Natal: Editora Realize, 2016. 8p.
63
Capítulo 7
João Victor Nunes de Sousa

Antonio Gilson Barbosa de Lima
Resumo: Os AUVs formam uma classe avançada de veículos submarinos, capazes

de realizar missões pré-estabelecidas, sem comunicação física com a terra ou
intervenção humana. Durante os últimos anos a pesquisa e o desenvolvimento
dessa classe de veículos vêm crescendo, devido às suas excelentes
características, ideais para a sua aplicação em setores como o militar, o científico
ou o industrial. Diante disso, esse artigo tem como objetivo avaliar a influência da
profundidade de operação de um AUV no coeficiente de arrasto do seu casco,
visando obter dados que auxiliem pesquisadores e projetistas. Para essa análise foi
utilizado o software ANSYS-CFX 15.0, que é uma poderosa ferramenta CFD para a
solução de problemas envolvendo a Mecânica dos Fluidos.
Palavras-chave: Veículo Autônomo Submarino, Coeficiente de Arrasto, Simulação

Numérica
64
1 INTRODUÇÃO em pouco tempo e com baixo custo quando

comparado às técnicas experimentais.
Define-se como Veículo Autônomo Submarino,
ou AUV (sigla em inglês para Autonomous Atualmente, as ferramentas de Fluido-
Underwater Vehicle), um veículo que viaja Dinâmica Computacional, ou CFD (sigla em
submerso, sem comunicação física com a inglês para Computational Fluid Dynamics),
terra e sem a necessidade de operador são integradas com outras ferramentas
humano. Os AUVs estão inseridos no grupo numéricas, criando um ambiente de projeto
dos Veículos Submarinos Não-Tripulados, completo, onde os experimentos são feitos
mais conhecidos como UUVs (sigla em inglês apenas para ajustes finais e testes.
para Unmanned Underwater Vehicles).
A maioria dos AUVs existente utiliza baterias
Durante as últimas décadas, vários AUVs têm como fonte de energia para o seu sistema de
sido desenvolvidos e pesquisas na área são propulsão. Quanto maior o arrasto gerado
cada vez mais frequentes, tendo em vista as durante o deslocamento do veículo, mais
características extremamente favoráveis que energia esse sistema necessitará e,
esses veículos possuem, tais como a consequentemente, menor será a autonomia
possibilidade de operar de forma autônoma do AUV, o que é indesejável por qualquer
em ambientes hostis, como áreas projetista.
inexploradas, águas de território inimigo,
Nesse contexto, esse artigo avalia,
áreas contaminadas ou águas profundas.
numericamente, a influência da profundidade
Todas essas características tornam o uso dos
de operação, que é um parâmetro
AUVs muito interessante para os setores
operacional básico, no coeficiente de arrasto
militar, científico e industrial.
do casco de um AUV, visando obter dados
A aplicação de AUVs vem sendo que auxiliem pesquisadores e projetistas em
impulsionada em virtude dos constantes trabalhos/concepções futuras.
avanços na eletrônica e na robótica,
permitindo a execução de missões de alta
precisão, a diminuição do custo e o aumento 2 METODOLOGIA
do nível de processamento dos computadores
2.1 GEOMETRIA DO CASCO
embarcados, bem como devido às novas
tecnologias de baterias e gerenciamento de O casco de AUV analisado é do tipo torpedo,
energia, possibilitando o aumento da ou seja, corpo cilíndrico com alta razão entre
autonomia, manutenção e segurança na seu comprimento e diâmetro, com perfis
operação desses veículos (Dantas, 2014). suaves em sua proa e popa, devido às suas
boas características de projeto, tais como
Visando manter os AUVs como ferramentas
baixo arrasto hidrodinâmico, bom volume
de alto desempenho, seus projetos devem
interno, simplificação do acesso a toda a
passar por constantes evoluções, refletindo
aparelhagem existente e redução do
necessidades presentes e futuras. Para tanto,
custo de fabricação (Madham et al., 2006),
se faz necessário o desenvolvimento teórico e
esse tipo de casco é utilizado atualmente
tecnológico nas mais diversas áreas de
pelos principais fabricantes comerciais dessa
projeto desses veículos, tais como na
classe de veículos.
mecânica, hidrodinâmica, navegação e
controle. Para modelar os perfis da proa e popa do
veículo recorreu-se às equações de Myring.
O uso de técnicas numéricas para a solução
Essas equações teóricas descrevem
de problemas de engenharia é hoje uma
contornos para as extremidades de corpos do
realidade graças ao desenvolvimento de
tipo torpedo que geram o menor coeficiente
computadores de alta capacidade de
de arrasto possível (Myring, 1976). As
processamento. Essas técnicas apresentam
equações são da seguinte forma:
inúmeras vantagens, podendo resolver
problemas complexos, sem soluções
analíticas exatas, e que muitas vezes não
podem ser reproduzidos em experimentos,
65
a) Proa:
1⁄
1 𝑥1 −𝑎 2 𝑛
𝑟1 (𝑥1 ) = 𝐷 [1 − ( ) ] (1)
2 𝑎
b) Popa:
1 3𝐷 𝑡𝑔𝜃 2 𝐷 𝑡𝑔𝜃 3
𝑟2 (𝑥2 ) = 𝐷 [ − ] 𝑥2 + [ − ] 𝑥2 
2 2𝑐 2 𝑐 𝑐3 𝑐2
onde todos os parâmetros de alimentação das equações, com exceção do parâmetro potencial n,
são geométricos, e estão mostrados na Figura 1.

Figura 1: Representação dos parâmetros geométricos necessários para alimentação das equações
de Myring
O casco analisado tem suas características de estarem disponíveis na literatura dados

geométricas baseadas nas do casco experimentais que possibilitaram a validação
Afterbody1 (Tab. 1), estudado por Jagadeesh da modelagem numérica-matemática utilizada
et al. (2009). Essa escolha se deu em virtude no presente artigo.
Tabela 1: Parâmetros geométricos do casco de AUV estudado

Lt = a + b + c (mm) 1.400
a (mm) 250
b (mm) 700
c (mm) 450
D (mm) 140
n (-) 2
θ (°) 20
2.2 DOMÍNIO FÍSICO E MALHA NUMÉRICA As dimensões desse domínio (Fig. 2) foram
escolhidas baseadas no trabalho de Sousa et.
Optou-se por um domínio fluido de forma
al. (2014), que estudou o escoamento ao
semicilíndrica, visando a obtenção de uma
redor do casco de um AUV tipo torpedo com
malha numérica com menor quantidade de
Lt ⁄ D=9, e são informadas na Tabela 2. Para
volumes de controle, quando comparado a
geração da geometria que representa o
um domínio paralelepípedo, reduzindo assim
domínio estudado, bem como para
o custo computacional para solução do
elaboração da malha numérica (Fig. 3), foi
problema.
utilizado o software comercial ICEM-CFD 15.0.
66
Já para as simulações, foi utilizado o software elementos finitos para solução do problema
comercial ANSYS-CFX 15.0, que faz uso do em estudo.
método dos volumes finitos baseado em
Figura 2: Domínio fluido de estudo
Tabela 2: Parâmetros geométricos que definem as dimensões do domínio fluido estudado

A (mm) 0,5 Lt = 700
B (mm) 3 Lt = 5.200
R (mm) 3 D = 420
Figura 3: Malha numérica criada sobre o domínio de estudo, que contempla 454.950 volumes de
controle, dentre tetraedros e prismas
Muitos cuidados devem ser tomados ao se equação que estima essa espessura
construir malhas numéricas para resolver (Equação 3), sendo recomendada a
problemas de escoamento externo, sendo construção de uma malha com grande
necessário levar em conta diversos refinamento nessa região (Fig. 4), visando
parâmetros nessa construção, sendo um capturar precisamente os efeitos viscosos
deles a espessura aproximada da camada- atuantes ao redor do casco do AUV.
limite hidrodinâmica, 𝛿. Ansys (2014) traz uma
67
−1⁄
𝛿 = 0,035𝐿𝑅𝑒 7 ⁡⁡⁡ (3)
onde:
𝐿 é um comprimento de referência; e
𝑅𝑒 é o número de Reynolds, Re, que pode ser definido como:
𝜌𝑈𝐿
𝑅𝑒 = (4)
µ
onde:
𝜌 é a massa específica;
𝑈 é a velocidade de referência do escoamento; e
µ é a viscosidade dinâmica.
O número de Reynolds pode ainda ser chamado de número de Reynolds

definido em base volumétrica, passando a ser volumétrico, 𝑅𝑒𝑉 , calculado pela expressão:
1
𝜌𝑈𝑉 ⁄3
𝑅𝑒𝑉 = (5)
µ
onde
V é o volume do corpo sujeito a escoamento externo.
2.3 MODELAGEM MATEMÁTICA incompressível, isotérmico e turbulento,

apresentando camada-limite hidrodinâmica
Para investigar o escoamento ao redor do
totalmente turbulenta.
casco do AUV analisado, foi considerado que
esse escoamento é monofásico, As equações gerais de conservação que
tridimensional, em regime permanente, regem o problema são:
Figura 4: Detalhe da proa do casco de AUV analisado (em vermelho), com destaque para o
refinamento da malha na região da camada-limite hidrodinâmica
68
A) Equação da conservação da massa:
𝜕𝜌
⃗ ) = 0⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡
+ ∇ ∙ (𝜌𝑈 (6)
𝜕𝑡
B) Equação da conservação da quantidade de movimento linear:
⃗)
𝜕(𝜌𝑈 T
⃗ ⊗U
+ ∇ ∙ (𝜌U ⃗ ) − ∇ ∙ (µ𝑒𝑓𝑓 ∇U
⃗ ) = −∇𝑝′ + ∇ ∙ (µ𝑒𝑓𝑓 ∇U
⃗ ) + 𝜌𝑔⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡(7)
𝜕𝑡
onde:
µ𝑒𝑓𝑓 é a viscosidade efetiva, que contabiliza os efeitos turbulentos, e é definida como:
µ𝑒𝑓𝑓 = µ + µ𝑡 (8)
Onde
⁡µ𝑡 ⁡é a viscosidade turbulenta;
𝑝′ é a pressão corrigida; e
𝑔⁡é o vetor aceleração da gravidade.
Além das equações básicas supracitadas, se baseado na modelagem k-ω. Utilizou-se esse
faz necessário adicionar ao modelo equações modelo devido ao seu bom tratamento de
que regem o fenômeno da turbulência problemas envolvendo escoamento externo
presente no escoamento. A turbulência com altos números de Reynolds e à sua boa
consiste de flutuações no campo de fluxo no precisão em prever a separação do fluxo
tempo e no espaço. É um processo complexo quando o mesmo é submetido a gradientes
e pode ter um efeito significativo sobre as de pressão adversos (Ansys, 2014).
características do escoamento. A turbulência
O modelo SST introduz duas variáveis ao
ocorre quando as forças de inércia agindo no
problema, que são a variável k, que é a
fluido se tornam significativamente maiores
energia cinética turbulenta, e a variável ω, que
que às forças viscosas, e é caracterizada por
é a frequência turbulenta. Essas variáveis são
um elevado número de Reynolds do
calculadas através das equações:
escoamento.
Foi utilizado o modelo de turbulência SST
(Shear Stress Transport). Esse modelo é
C) Equação da energia cinética turbulenta:

𝜕(𝜌𝑘) µ𝑡
⃗⃗ 𝑘) = ∇ ∙ [(µ +
+ ∇ ∙ (𝜌U ) ∇𝑘] + 𝑃𝑘 − 𝐶𝑘2 𝜌𝑘𝜔⁡⁡⁡ (9)
𝜕𝑡 𝐶𝑘1
69
D) Equação da frequência turbulenta:
𝜕(𝜌𝜔) µ𝑡 𝜔
⃗⃗ 𝜔) = ∇ ∙ [(µ +
+ ∇ ∙ (𝜌U ) ∇𝜔] + 𝐶𝜔2 𝑃𝑘 − 𝐶𝜔3 𝜌𝜔2 ⁡⁡⁡ (10)
𝜕𝑡 𝐶𝜔1 𝑘
onde:
𝐶𝑘1 , 𝐶𝑘2 , 𝐶𝜔1 , 𝐶𝜔2 e 𝐶𝜔3 são constantes, e valem, respectivamente, 2; 0,09; 2; 5/9; e 0,075; e
𝑃𝑘 é a produção de turbulência, que é modelada por meio da equação:
T 2
𝑃𝑘 = µ𝑡 ∇U
⃗⃗ ∙ (∇U
⃗⃗ + ∇U ⃗⃗ (3µ ∇ ∙ U
⃗⃗ ) − ∇ ∙ U ⃗⃗ + 𝜌𝑘) (11)
3 𝑡
O modelo SST assume que a viscosidade turbulenta e à frequência turbulenta através

turbulenta está associada à energia cinética da relação:
𝑘
µ𝑡 = 𝜌 𝜔 (12)
O modelo SST assume também que a pressão

corrigida, p', é modelada através da equação:
2
𝑝′ = 𝑝 + 𝜌𝑘 (13)
3
Onde
p é a pressão.
Sempre que existir uma velocidade relativa do movimento relativo entre o sólido e o fluido
entre um corpo sólido e o fluido que o é definida como força de arrasto,⁡𝐹𝑑 .
circunda surgirá uma força resultante 𝐹 sobre
Essa força de arrasto é comumente expressa
este corpo (Fox et al., 2011). O módulo,
em termos de um número adimensional, o
direção e sentido dessa força dependerão de
coeficiente de arrasto, 𝐶𝑑 . Em base
diversos fatores, tais como a velocidade
volumétrica, esse coeficiente passa a se
relativa entre o sólido e o fluido, a forma e
chamar coeficiente de arrasto volumétrico, e
dimensões do sólido e as propriedades do
pode ser definido da seguinte forma:
fluido. A componente de 𝐹 paralela à direção
𝐹𝑑
𝐶𝑑𝑣 = 1 2 (14)
𝜌𝑉 3 𝑈 2
⁄
2
2.4 CONDIÇÕES DE CONTORNO E Foram utilizados os seguintes fluidos nas

PROPRIEDADES DOS FLUIDOS simulações:
A Figura 5 mostra informações relativas às  Água pura (Tab. 4): Para as simulações
condições de contorno adotadas nas de validação; e
simulações.
70
 Água do mar: Para as demais simulações, de operação no coeficiente de arrasto do

que avaliam a influência da profundidade casco de AUV estudado.
Tabela 4: Propriedades da água pura utilizadas nas simulações*

Densidade (kg/m³) 997
Viscosidade dinâmica (mPa.s) 0,8899
Fonte: Ansys, 2014
Figura 5: Condições de contorno do domínio
Schmiegelow (2004) informa que, para c) Camada de fundo: Região de

regiões equatoriais e tropicais, o profundidades superiores a 1.000 m, onde
comportamento médio da temperatura da não há grandes variações na temperatura da
água do mar em função da profundidade água (região praticamente isotérmica).
oceânica pode ser dividido em três regiões:
Baseado no perfil vertical médio de
a) Camada de mistura: Região de até temperatura das águas oceânicas, para
500 m de profundidade, onde não há grandes regiões equatoriais e tropicais, e
variações na temperatura da água (região considerando a salinidade média dessas
isotérmica); águas como sendo 35 g/kg, chegou-se às
propriedades da água do mar descritas na
b) Termoclina principal: Região entre 500
Tabela 5, que foram utilizadas nas
e 1.000 m, onde há grandes gradientes de
simulações.
temperatura; e
Tabela 5: Propriedades da água do mar utilizadas nas simulações*

Profundidade (m) Temperatura (°C) Densidade (kg/m³) Viscosidade dinâmica (mPa.s)
500 25 1.023 0,9600
750 15 1.026 1,2300
1.000 3 1.028 1,7200
Fonte: The Engineering Toolbox, 2018
71
RESULTADOS E DISCUSSÃO simulada computacionalmente e a utilizada

nos experimentos.
3.1 VALIDAÇÃO DAS SIMULAÇÕES
Para validação da modelagem numérica-
matemática utilizada no presente trabalho, 3.2 AVALIAÇÃO DA INFLUÊNCIA DA
comparou-se resultados experimentais PROFUNDIDADE DE OPERAÇÃO NO
obtidos por Jagadeesh et al. (2009) com as COEFICIENTE DE ARRASTO
simulações executadas no presente artigo,
A Tabela 6, a Figura 7 e a Equação 15 trazem
obtendo uma aproximação satisfatória entre
os dados referentes ao coeficiente de arrasto
os coeficientes de arrasto volumétricos para
volumétrico em função da profundidade de
em ambas as situações, com um desvio
operação do casco de AUV analisado, para
máximo inferior a 10% em todos os casos
uma velocidade de escoamento de 1,0 m/s.
(Fig. 6). Baseado nesses resultados, confirma-
Verifica-se um comportamento linearmente
se que a modelagem numérica-matemática
crescente do coeficiente de arrasto
adotada representa bem o fenômeno do
volumétrico com o aumento da profundidade.
escoamento ao redor do casco de AUV
Esse fenômeno pode ser explicado em virtude
analisado, tendo sido observado um desvio
do aumento significativo da viscosidade da
médio de 7,2% entre os resultados numéricos
água do mar, que se dá devido à queda
e os experimentais. Esse pequeno desvio
considerável da temperatura à medida que se
pode ser justificado em virtude de erros
avança para as camadas oceânicas mais
experimentais, erros numéricos e/ou devido à
profundas.
pequena variação entre a geometria do casco
Tabela 6: Coeficiente de arrasto volumétrico em função da profundidade de operação do casco de

AUV
Profundidade (m) Re x 10-5 (-) ReV x 10-5 (-) Cdv (-)
500 14,92 2,71 0,0399
750 11,68 2,12 0,0421
1.000 8,37 1,52 0,0447
Figura 6: Comparativo entre os coeficientes de arrasto volumétricos em função do número de

Reynolds volumétrico, obtidos experimentalmente e numericamente
0,055
Experimental - Jagadeesh et al. (2009)

Numérico
0,050
0,045
Cdv (-)
0,040
0,035
0,030
1 2 3 4
Rev x 10-5 (-)
72
Figura 7: Coeficiente de arrasto volumétrico em função da profundidade de operação do casco de

AUV analisado
0,050
Equação 15
Cdv (-) 0,045
0,040
0,035
400 500 600 700 800 900 1000 1100
Profundidade (m)
𝐶𝑑𝑣 = 0,0351 + 0,0956(𝑝𝑟𝑜𝑓𝑢𝑛𝑑𝑖𝑑𝑎𝑑𝑒. 10−4 ) (15)
onde o coeficiente de correlação foi de 99,61%.
4 CONCLUSÃO se passa de 500 m para 1.000 m de

profundidade de operação.
Baseado nos resultados expostos, conclui-se
que a modelagem numérica-matemática
utilizada no presente artigo representou bem
AGRADECIMENTOS
o fenômeno do escoamento ao redor do
casco de AUV estudado. Conclui-se ainda Os autores agradecem ao Programa de Pós-
que a profundidade de operação do AUV, na Graduação em Engenharia Mecânica
faixa entre 500 e 1.000 m, tem forte influência (PPGEM) da Universidade Federal de
sobre o arrasto do seu casco, e Campina Grande (UFCG), ao Laboratório
consequentemente, sobre o arrasto total do Computacional de Térmica e Fluidos (LCTF) e
veículo, tendo sido verificado um incremento a todos os autores das Referências
de cerca de 12% no arrasto do casco quando Bibliográficas por viabilizarem a produção do
presente trabalho.
REFERÊNCIAS
[1] Ansys, “CFX-Theory Manual”. Ansys Inc. [4] Jagadeesh, P.; Murali, K.; Idichandy, V.G.
2014. “Experimental Investigation of Hydrodynamic Force
Coefficients over AUV Hull Form”. Ocean
[2] Dantas, J.L.D. “Metodologia Numérica
Engineering, 36, 113-118, 2009.
para Estimativa da Manobrabilidade de Veículos
Autônomos Submarinos”. Tese de doutorado em [5] Madham, R.; Desa, Elgar.; Prabhudesai,
Engenharia Mecânica. Universidade de São Paulo, S.; Sebastião, L.; PascoaL, A.; Desa, Ehrlich;
2014. Mascarenhas, A.; MAURYA, P.; Navelkar, G.;
Afzulpurkar, S.; Khalap, S. “Mechanical Design and
[3] Fox, R.W.; McDonald, A.T.; Pritchard, P.J.
Development Aspects of a Small AUV – MAYA”.
“Introdução à Mecânica dos Fluidos”. 7ª edição.
National Institute of Oceanography (India) and
LTC Editora, 2011.
Institute for Systems and Robotics (Portugal), 2006.
73
[6] Myring, D. F. “A Theoretical Study of Body Turbulent Fluid Flow and Drag Coefficient for
Drag in Subcritical Axisymmetric Flow”. Optimizing the AUV Hull Design”. Open Journal of
Aeronautical Quarterly, 27(3):186–94, August 1976. Fluid Dynamics, 2014, 4, 263-277.
[7] Schmiegelow, J.M.M. “O Planeta Azul – [9] The Engineering Toolbox. “Sea Water
Uma Introdução às Ciências Marinhas”. 1ª edição. Properties”. Disponível em
Editora Interciência, 2004. <https://www.engineeringtoolbox.com/sea-water-
properties-d_840.html>. Acesso em: 28 de set. de
[8] Sousa, J.V.N.S.; Macêdo, A.R.L.; Amorim
2018.
Junior, W.F.; Lima, A.G.B. “Numerical Analysis of
74
Capítulo 8

Renan Di Pace Arruda
Gustavo Lima do Nascimento
Juscelino de Farias Maribondo
Resumo: O objetivo deste trabalho é detectar danos de impacto em materiais

compósitos via Termografia Infravermelha com o auxílio de uma câmara de
isolamento térmico. Os métodos utilizados no desenvolvimento desta pesquisa
foram a fabricação de compósitos poliméricos laminados e posterior ensaio de
impacto utilizando uma máquina por queda de peso, a fabricação de câmara de
isolamento térmico, aquecimento por exposição à lâmpadas halógenas e execução
do ensaio não destrutivo através da termografia ativa. Os resultados encontrados
mostram que a Termografia Infravermelha é capaz de detectar os danos de
impacto em materiais compósitos poliméricos e que existe influência no uso do
isolamento térmico e na forma geométrica na percepção dos danos de impacto.
Palavras-chave: Compósitos, Dano de impacto, Termografia Infravermelha,

Aplicação Aeronáutica.
75
1 INTRODUÇÃO fibra de vidro ou com fibra de carbono.

Segundo (MEOLA et al., 2016) no mundo de
O avanço da utilização de materiais
hoje há uma demanda crescente de
compósitos na indústria tem demandado um
estruturas leves e de alto desempenho que
amplo estudo em torno de suas propriedades.
explica o crescente interesse em relação aos
Em certas aplicações, esses materiais se
materiais compósitos. Sua combinação pode
apresentam como uma alternativa
oferecer propriedades especificas adequadas
interessante em relação aos materiais
ao uso, em termos de resistências mecânica e
metálicos. Dentre os pontos fortes da sua
rigidez por peso específico. Dessa forma, o
utilização se destaca a alta resistência
uso de materiais compósitos na indústria
específica, ou seja, esses materiais são
aeronáutica aumentou exponencialmente nos
capazes de fornecerem alta resistência
últimos anos. Grande parte dos equipamentos
mecânica associado a uma baixa densidade,
utilizados na fuselagem, ailerons e spoilers,
promovendo assim uma redução do uso de
por exemplo, são feitos a partir de materiais
combustível e um aumento na capacidade de
compósitos. Em casos como o do Boeing 787,
carga das aeronaves.
cerca de 50% da sua estrutura é feita através
Esses compósitos são, em geral, formados a de compósitos avançados (HALE, 2008).
partir de duas fases distintas. A primeira é
Porém, a combinação carbono-epóxi e vidro-
denominada como reforço, que em geral são
epóxi apresenta resistência baixa aos
fibras, e têm a função de melhorar o
impactos e a fadiga. Isso influi nos danos
desempenho mecânico. A segunda fase do
causados por baixa velocidade, que acabam
compósito é dada por um material ligante,
por gerar danos de difícil verificação nos
conhecida como matriz, que é responsável
laminados, o que pode acarretar em reduções
pela união dos demais constituintes e de
significativas em suas propriedades
transmitir os esforços entre as fibras. Além
mecânicas (ALMEIDA, 2010).
disso, a matriz tem a função de proteger as
fibras individuais contra danos superficiais Devido a esses problemas, geralmente são
que são capazes de causar a formação de utilizados critérios de resistência tanto para
trincas (CALLISTER; RETHWISCH, 2007). seleção do material quanto para o
carregamento, o que acaba acarretando em
Além das propriedades específicas já
custos extras (CUARTAS et al., 2014). Esses
declaradas, esses materiais possuem, em
defeitos geralmente geram danos superficiais
geral, elevada resistência à fadiga, alta
ou subsuperficiais e são causados tanto no
resistência a corrosão, elevada estabilidade
processo de construção, como na utilização e
térmica e baixa condutibilidade, além de
manutenção de aeronaves. Estes danos
serem dados como excelentes amortecedores
variam em relação ao módulo de energia e ao
para vibração (JUNIOR, 2014). Entretanto,
tipo de impacto. A detecção visual desses
apresentam uma grande perda de
defeitos é complicada devido à sua pequena
propriedades na presença de algum defeito,
dimensão.
que são muitas vezes microscópicos. Essa
incerteza no seu comportamento, portanto, faz As aeronaves, quando estão em condições
com que tenham grandes dispersões reais de funcionamento, seja em voo ou em
associadas a durabilidade de estruturas que manutenção, estão sujeitas a diversas fontes
utilizam esses materiais. Além disso, esses de impacto em sua estrutura. (ALMEIDA,
materiais são sensíveis em relação a seu 2010) mostra que estes danos podem ser
processo de produção, possuem uma análise originados desde choque com pássaros,
complexa e uma alta sensibilidade em relação chuva de granizo, fragmentos levantados na
ao dano sofrido. O dano de impacto pode decolagem ou aterrisagem, queda de
reduzir significativamente as propriedades ferramenta até o contato com outra aeronave.
mecânicas residuais de compósitos, embora Na Figura 1 pode-se ver situações onde a
dificilmente visível. (ZHANG et al., 2018). aeronave está sujeita a danos de impacto e
também as consequências dos mesmos.
Grande parte da utilização de materiais
compósitos na indústria aeronáutica é dada
pela utilização do laminado reforçado com
76
Figura 1 – (a) Choque com pássaros durante uma decolagem; (b,c) Efeitos do impacto de chuva de
granizo e consequências dos danos de impacto
a) b) c)
Fonte - (ALMEIDA, 2010)
Nesse contexto, o setor aeronáutico utiliza feitos através da utilização de uma fonte de
cada vez mais de ensaios não-destrutivos calor externa. Esse processo é denominado
para a visualização rápida e precisa dos por termografia ativa e é visto que o objeto de
materiais compósitos. A utilização do ensaio estudo deve ser energeticamente estimulado
termográfico é ainda um ensaio em (por exemplo, através de fontes térmicas
desenvolvimento na indústria aeronáutica. simples como lâmpadas, flashes e ar quente,
Utilizando de radiação térmica para detecção ou então, por meios mais elaborados, tais
de falhas nos compósitos, essa técnica vem como ondas ultrassônicas, correntes
sendo cada vez mais estudada pela sua parasitas, micro-ondas e laser e células de
eficiência e facilidade de operação. Peltier), de forma que um fluxo interno de
calor seja gerado na parte inspecionada.
A termografia infravermelha é baseada na
medida do fluxo térmico ou radiação Estes tipos de estratégias de aquecimento
eletromagnética emitida por uma peça por possibilitaram o acompanhamento dos
ondas infravermelhas (FRAMEZELLE et al., fenômenos térmicos desenvolvidos na parte
2011). Este fluxo é diretamente proporcional a inspecionada, permitindo a identificação e a
temperatura da superfície da peça e, com o caracterização de danos impingidos por
auxílio de uma câmera infravermelha é impacto aos laminados. Uma vez que, a
possível fazer a representação térmica do presença de defeitos na região inspecionada
corpo em questão. A energia infravermelha é irá causar uma perturbação desse fluxo,
proveniente da vibração e agitação das fazendo com que a câmera termográfica
moléculas, ou seja, um corpo com maior acuse o aparecimento de um gradiente
temperatura terá uma maior agitação térmico na descontinuidade do componente
molecular emitindo assim uma maior inspecionado (TARPANI et al., 2009).
quantidade de energia infravermelha.
Dessa forma, esse trabalho tem o objetivo de
As principais vantagens da utilização da analisar danos de impacto em materiais
termografia infravermelha são: a rapidez da compósitos, através da técnica da termografia
inspeção, a interpretação simples das ativa. Os corpos de prova foram fabricados e
imagens, a análise pode ser feita em tempo impactados com danos provenientes de
real, a radiação emitida pela peça é não letal, diferentes níveis de energia de impacto. As
a geometria da peça não se apresenta imagens obtidas foram analisadas a partir de
sensível na detecção do dano, pode ser diferentes tipos de câmaras de isolamento
utilizado em grandes superfícies e pode ser térmico, identificando as diferenças nas
feita também a inspeção total da peça. As imagens, a partir de um aquecimento
poucas desvantagens do método incluem: a semelhante para todos os casos.
falta de uniformidade de aquecimento eo
custo elevado do equipamento. (ALMEIDA,
2010). 2 MATERIAIS E MÉTODOS
Devido a esses fatores, a utilização da Segundo (GIL, 2002) a pesquisa se classifica
termografia tem se tornado bastante segundo a natureza como aplicada em função
abrangente, visando uma promissora redução de ter como objetivo gerar conhecimento para
de custos na manutenção. Na indústria uma aplicação prática e buscar a solução do
aeronáutica, por exemplo, os ensaios são problema em questão. Neste caso, trata-se da
77
avaliação da capacidade da técnica da Também (GIL, 2002) classifica a pesquisa

Termografia detectar danos em materiais com relação aos procedimentos técnicos,
compósitos. Por tratar-se de uma pesquisa este estudo trata-se de uma pesquisa
quantitativa, as informações são expressas experimental, pois a partir da determinação
em valores numéricos, onde serão então de um objeto de estudo, como a detecção de
classificadas e analisadas para a obtenção danos em materiais compósitos, são
dos resultados. Para este trabalho, são selecionadas as variáveis que as influenciam,
obtidos valores de temperatura que indicam a suas formas de controle e observados seus
presença do dano a partir da temperatura do efeitos gerados.
corpo de prova estudado. Portanto, quanto
Todas as fases de realização deste trabalho
aos objetivos, esta pesquisa se enquadra
estão apresentadas em um fluxograma único
como descritiva, pois busca esclarecer um
(Fluxograma 1), o qual é detalhado em
determinado assunto já conhecido.
seguida.
Fluxograma 1 – Fluxograma da metodologia utilizada
Fonte - Autoria própria
78
A partir da necessidade esclarecida que Na Fase 2, realiza-se o ensaio de impacto.

levou ao desenvolvimento deste trabalho, Para tanto, utilizou-se as diretrizes
utilizou-se o processo metodológico estabelecidas pela norma ASTM – D7136M
apresentado no Fluxograma1, o qual é (Standard Test Method for Measuring the
composto por seis fases assim denominadas: Damage Resistance of a Fiber-Reinforced
Fabricar os corpos de prova (Fase 1), Realizar Polymer Matrix Composite to a Drop-Weight
ensaios de impacto (Fase 2), Fabricar a Impact Event), a máquina utilizada foi
câmara de isolamento térmico (Fase 3), projetada para ensaios por queda de peso
Preparar corpos de prova para os testes desenvolvida por (JUNIOR, 2014), conforme
termográficos (Fase 4), Registrar os mostrado na Figura 2 e encontra-se disponível
termogramas (Fase 5) e Analisar os para realização de pesquisas cientificas e
termogramas (Fase 6). acadêmicas na Unidade Acadêmica de
Engenharia Mecânica da Universidade
A Fase 1 se destina a fabricar os corpos de
Federal de Campina grande (UFCG) no
prova em material compósito. Para executa-la
Laboratório de Projetos e Instrumentação
utilizou-se a norma ASTM 7136D/D7136M e
laboratório (LPI).
fez-se uso da técnica de laminação manual.
As placas apresentam como reforço fibra de A máquina possui as seguintes propriedades:
vidro-E e como matriz a resina cristalina,
ambas obtidas pela Empresa EquiFiber  5 níveis de energia de impacto - 20, 39,
Equipamentos de Fibra LTDA., localizada no 58, 77 e 96 J;
distrito industrial de Campina Grande, Estado  Altura de queda dos pesos de 1,3 metros;
da Paraíba, Brasil.
 Impactador de ponta semiesférica com
diâmetro de 16 mm e dureza 57 HRC.
Figura 2 – Máquina de impacto por queda de peso.
Fonte – (JUNIOR, 2014)
Na Fase 3 fabrica-se a câmara de isolamento para isso inicialmente utilizou-se as diretrizes

térmico, sendo assim foi necessário projetar indicadas pelo manual da câmera
uma câmara de isolamento térmico, que termográfica na qual indica os parâmetros
comportasse o objeto de estudo com o intuito que devem ser ajustados antes dos registros,
de que o livrasse das interferências de tais como: Emissividade, Temperatura
radiações térmicas do ambiente. Aparente Refletida, Distância, Umidade
relativa e Temperatura atmosférica. Por optar
Na Fase 4 foram preparados os corpos de
pela técnica da termografia ativa (Figura 3) na
prova ensaiados para os testes termográficos,
79
qual o objeto de estudo é energeticamente para o de aquecimento dos espécimes, assim

estimulado de forma que um fluxo interno de como o tempo de aquecimento e a distância
calor seja gerado na parte inspecionada, se da fonte para o corpo de prova impactado.
fez necessário a escolha da fonte de calor
Figura 3 – Termografia Ativa
Fonte - (TARPANI et al., 2009)
Para determinação da temperatura aparente mensuração das condições ambientais, tais

refletida utilizou-se o método direto como umidade relativa e temperatura
apresentado na norma ABNT NBR 16292-14 atmosférica utilizou-se a Sonda AMI 300
para precisar a distância entre o CP e a (Figura 4), esse equipamento foi
câmera termográfica utilizou-se uma fita disponibilizado pelo Laboratório de
métrica convencional. A emissividade é Ergonomia e Engenharia de Segurança do
determinada pelo próprio manual da câmera Trabalho (LEEST) da UFCG.
para vários tipos de materiais, para
Figura 4 – Sonda AMI 300.
Fonte – Autoria própria.
Em seguida monta-se a bancada do laboratório LEEST. Um croqui da montagem

experimento. Os ensaios foram realizados no do sistema pode ser visto na Figura 5.
80
Figura 5 – Ilustração da bancada de ensaios.
Fonte – Autoria própria.
Na Fase 5 registra-se os termogramas, para pertence ao Laboratório de Ergonomia e

isso se fez a utilização da câmera térmica Engenharia de Segurança do Trabalho
T620 da FLIR SYSTEMS (Figura 6), a qual (LEEST) da UFCG.
Figura 6 – câmera térmica T620 da FLIR SYSTEMS.
Fonte – Autoria própria
Na Fase 6 analisa-se os termogramas, para Na Fase 2, realizou-se o ensaio de impacto de

isso procura-se conferir, utilizando o software queda de peso onde os níveis de energia
FLIR tools e gráficos por meio do Microsoft foram selecionados de acordo com as
Excel, os dados obtidos com os parâmetros pesquisas realizadas na literatura. Sabe-se
de referência para o estudo a fim de verificar que uma aeronave pode sofrer um dano
a existência de anomalias durante o devido a queda de um martelo, isso pode
acionamento do sistema. Após a finalização ocorrer quando o avião está em processo de
da Fase 6 realizou-se a redação deste manutenção, até choque com pássaro
trabalho. quando a aeronave está em voo. Para simular
impactos reais, determinou-se quatro níveis
de energia a partir dos dados da Tabela 1 e
3 RESULTADOS E DISCUSSÕES das limitações da máquina de impacto por
queda de peso utilizada para realizar os
Na Fase 1 obteve-se a placa confeccionada
ensaios, realizou-se com os parâmetros da
com o intuito de utilização na máquina de
máquina dados da seguinte maneira:
ensaio por queda de peso. Portanto as suas
dimensões foram baseadas segundo a norma  4 níveis de energia de impacto – 20, 39,
ASTM D7136/D7136M-12. O corpo de prova 58 e 77 Joules;
utilizado apresenta as seguintes
características: 100 mm de largura; 150 mm  Altura de queda dos pesos de 1,3 metros;
de comprimento; e 2 mm de espessura.  Impactador de ponta semiesférica com
diâmetro de 16mm e dureza 57 HRC.
81
Tabela 1 – Tipos comuns de danos de impacto na estrutura de uma aeronave

Descrição Energia (J) Massa (g) Velocidade (m/s) Circunstâncias
Queda de equipamento 6 330 >6 Manutenção
Queda de elementos removíveis 4 220 >6 Movimento de carga
Componente de manutenção 16 910 >6 Manutenção
Granizo (até 51 mm) 43 62 37,3 Decolagem e voo
Choque com pássaros 3,8-81 (kJ) 1800 65-300 Decolagem e voo
Fragmentos de pista 2-40 9 20-94 Decolagem e voo
Carga concentrada 50 - Estática Movimento de carga
Fonte - (ALVES; CHAVES; BIRCH, 2003)
Os corpos de prova analisados foram o nível de energia que o objeto foi impactado.
nomeados com códigos específicos, sendo Assim, obteve-se a Tabela 2 que apresenta a
CP o prefixo que significa corpo de prova nomenclatura utilizada.
seguido por dois números, os quais significam
Tabela 2 – Nomenclatura utilizada para os corpos de prova

Corpos de Prova Nomenclatura
Corpo de prova impactado com 20 J CP20
Fonte – Autoria própria
A Figura 7 apresenta os danos obtidos de do defeito causado pelo impacto não foram
acordo com o aumento do nível de energia. analisadas, uma vez que o objetivo do
trabalho é verificar se é possível a detecção
É possível verificar que a delaminação cresce
dos danos de impacto através dos
conforme a energia que foi realizada no
termogramas.
ensaio aumenta. Vale ressaltar que, devido ao
caráter qualitativo do trabalho, as dimensões
Figura 7 – Região impactada após o ensaio. (a) CP20; (b) CP39; (c) CP58; (d) CP77 J.
Fonte - Autoria Própria
82
Na Fase 3, a concepção da câmara de A câmara possui uma abertura na face frontal

isolamento térmico se deu a partir de para o encaixe da câmera termográfica
processo manual e a geometria selecionada [Figura 8(a)] e, a uma distância de 41 cm da
foi retangular devido a facilidade construtiva. face frontal se localiza o corpo de prova em
Para o dimensionamento da estrutura se levou análise, pois está é a distância focal mínima
em consideração o comprimento focal mínimo que o fabricante da câmera sugere [Figura
estabelecido pelo fabricante da câmera 8(b)].
termográfica que é de 42 cm, a altura do CP é
Em relação ao projeto da câmara térmica,
de 15 cm e a largura do CP é de 15 cm.
vale destacar que foi projetada com a parte
Diante dessas condições determinou-se uma
superior móvel. Dessa forma, foi possível
câmara retangular com as seguintes
realizar, de forma mais fácil, o posicionamento
dimensões: 62 cm de comprimento, 17 cm de
do corpo de prova aquecido no local
largura e 22 cm de altura, pois essa
especificado para que a região do dano fosse
concepção possibilitou o encaixe e retirada
centralizada. Além disso, essa parte móvel foi
do CP facilmente da estrutura.
projetada para ser facilmente realocada após
o posicionamento do corpo de prova.
Figura 8 – Câmara de isolamento térmico: (a) Câmara fechada; (b) Câmara aberta.
Fonte - Autoria Própria
O material utilizado para sua concepção foi o das partes da câmara foi utilizado cola quente
Poliestireno Expandido (conhecido como de silicone com o auxílio de pequenos
Isopor), pois dificulta a dissipação de calor e bastonetes de madeira assim como, a
permite manter a temperatura ambiente. Para inserção do isolamento interno e externo.
o revestimento interno foi utilizado o Etil-Vinil-
Na Fase 4, as condições ambientais
Acetato (conhecido por EVA) de cor preta e
encontradas no laboratório foram mensuradas
emborrachado (fosco), pois esse material atua
pelo aparelho descrito na metodologia,
como uma barreira radiante amortecendo o
obteve-se a temperatura ambiente média de
fluxo de calor e conserva a temperatura
23°C e umidade relativa média de 57%. Ainda
interna de um equipamento em que os
foi medido a velocidade da corrente de ar em
materiais empregados apresentam baixa
torno dos compósitos expostos ao ambiente
condutividade. De forma a reduzir o ganho de
que foi de 0,0 m/s por se encontrar
calor pelos corpos circunvizinhos foi então
enclausurado dentro da câmara de
necessário o uso de Papel Alumínio, pois é
isolamento térmico. Uma propriedade
um material de alta refletância e baixa
importante para a medição termográfica é a
emissividade, (SILVA, 2017). Para a junção
Emissividade, para este estudo utilizou-se um
83
valor de 0,94 de emissividade, pois segundo aquecer a superfície dos corpos de prova. O
o manual da máquina utilizada FLIR T620 CP foi aquecido por duas lâmpadas
esse é o valor que o material compósito halógenas de 500W (totalizando 1.000W), por
laminado com fibra de vidro apresenta. um tempo de 60 segundos. Após isto, o CP foi
transferido para o interior da câmara de
Como fonte de calor utilizou-se lâmpadas
isolamento térmico de forma a não ocorrer a
halógenas, as quais foram usadas para
perda de calor para o ambiente.
Figura 9 – Posicionamento metodológico de aquecimento do corpo de prova.
CP
Fonte - Autoria própria.
Na Fase 5 realizou-se o registro dos estava posicionada na câmara de isolamento

termogramas, uma vez que, a câmera já como apresenta a Figura 10.
Figura 10 – Posicionamento da câmera na câmara de isolamento térmico.
Fonte - Autoria própria
Os termogramas foram obtidos para os Tabela 3 pode-se ver os termogramas

impactos obtidos pelos quatro níveis de qualitativos e quantitativos registrados.
energia estabelecidos anteriormente. Na
84
Tabela 3 – Termogramas dos corpos de prova danificados por diferentes níveis de energia de
impacto.
Fonte - Autoria Própria.
Na Fase 6 se analisou os termogramas, é visto Uma análise qualitativa nos termogramas

nas imagens térmicas na Tabela 3 que o observa-se que a região do impacto quando
aquecimento por radiação térmica através da estimulada dissipa calor mais rápido ou
exposição da placa compósita às lâmpadas absorve menos energia, portanto há um
halógenas obteve-se uma distribuição decréscimo na leitura da temperatura em
semelhante em todos os casos. Os corpos de função do maior dano, ou seja, quanto maior o
prova atingiram uma temperatura máxima de impacto, maior o dano e menor a temperatura
aproximadamente 45°C com o tempo de ensaiada como pode ser observado no
exposição ao calor, e então o dano originado Gráfico 1. É possível visualizar nos
anteriormente por impacto é mostrado no termogramas da Tabela 3 que ao redor da
centro da placa por diferentes níveis de calor localização do impacto apresenta cores de
emitido (observado pela diferente coloração). temperatura maior do que na região
danificada.
85
Gráfico 1 – Temperatura na região de dano versus a energia de impacto.
TEMPERATURA X ENERGIA
46 44,6
44
TEMPERATURA (°C)
42
41
42
40 38,2
38
36
34
20 39 58 77
ENERGIA DE IMPACTO (J)
Fonte - Autoria Própria.
O dano originado no compósito pode ser Para tanto realizou-se uma revisão na
visto, segundo esta metodologia, a partir do literatura envolvendo os seguintes temas:
nível de energia de 39J, este nível de energia termografia, danos em materiais compósitos
corresponde ao mesmo nível de energia aeronáuticos e ensaio não-destrutivo para
quando uma aeronave sofre choque com um materiais compósitos.
pássaro de 0,5kg voando a 12m/s.
Após estabelecer a metodologia do trabalho,
Os termogramas obtidos do CP39, CP58 e apresentar e discutir os resultados, pode-se
CP77 mostram claramente a extensão da área concluir que com a realização do ensaio não-
danificada. O maior dano é mostrado no destrutivo de termografia é possível visualizar
último termograma da Tabela 3, onde a danos de impacto em materiais compósitos
variação de temperatura no local do dano é mesmo com pouco tempo de aquecimento
de alguns graus de diferença. dos corpos de prova, provando ser um ensaio
rápido e seguro para o cientista, também foi
possível observar que quanto maior for a
4 CONCLUSÕES energia do impacto melhor e mais rápida será
a detecção através deste ensaio e que
Ao longo deste trabalho buscou-se estudar o
quando as áreas de inspeção forem
seguinte problema: “analisar danos de
ensaiadas, as regiões onde houver danos
impacto em materiais compósitos, através da
serão detectadas por apresentarem cores de
técnica da termografia ativa”.
menor temperatura.
REFERÊNCIAS
[1] ALMEIDA, E. G. R. DE. Inspeção [5] FRAMEZELLE, G. et al. Contrôle Non
termográfica de danos por impacto em laminados Destructif De Structures Composites Pour
compósitos sólidos de matriz polimérica fortalecida L’Aéronautique. NDT.net, p. 1–15, 2011.
com fibras de carbono. [s.l.] Universidade de São
[6] GIL, A. C. Como Elaborar Projetos de
Paulo, 2010.
Pesquisa. 4a Edição ed. São Paulo - SP: [s.n.].
[2] ALVES, M.; CHAVES, C. E.; BIRCH, R. S.
[7] HALE, J. Boeing 787 from the Ground Up.
Impact on aircraft. 17th International Congress of
AeroMagazine, p. 9, 2008.
Mechanical Engineering, p. 1–8, 2003.
[8] JUNIOR, F. G. Desenvolvimento de uma
[3] CALLISTER, W.; RETHWISCH, D. Materials
máquina de impacto por queda de peso para
science and engineering: an introduction. 7th. ed.
materiais compósitos. [s.l.] Universidade Federal
New York, NY: [s.n.]. v. 94
de Campina Grande, 2014.
[4] CUARTAS, V. M. et al. Utilisation de la
[9] MEOLA, C. et al. New perspectives on
thermographie infrarouge et de l’emission
impact damaging of thermoset- and thermoplastic-
acoustique pour l’identification de
matrix composites from thermographic images.
l’endommagement d’un composite stratifié
Composite Structures, v. 152, p. 746–754, 2016.
carbone-époxy. OATAO, 2014.
86
[10] SILVA, J. F. DA. Estudo do mapeamento de matriz polimérica reforçados por fibras de
térmico das mãos para o auxílio ao diagnóstico da carbono. Polímeros, v. 19, p. 318–328, 2009.
LER. [s.l.] Universidade Federal de Campina
[12] ZHANG, H. et al. Thermographic Non-
Grande, 2017.
Destructive Evaluation for Natural Fiber-Reinforced
[11] TARPANI, J. R. et al. Inspeção Composite Laminates. Applied Sciences, v. 8, n. 2,
termográfica de danos por impacto em laminados p. 240, 2018.
[13]
87
Capítulo 9
Marcelio Ronnie Dantas de Sá

Aline Michelly do Nascimento
Antonio Almeida Silva.
Resumo: Neste trabalho é proposto a simulação de um aparato capaz de absorver

parte dos tremores em pacientes com a doença de Parkinson. Para isso é sugerido
a utilização de um sistema giroscópio, uma vez que esses dispositivos são
bastantes simples de serem construídos, além de seu alto poder de absorção de
vibração quando dimensionados da forma correta. Essa necessidade surge a partir
da carência de novos métodos de controle dos tremores em pessoas que sofrem
por conta desse sintoma, já que os medicamentos apresentam fortes efeitos
colaterais e as cirurgias envolvem grandes riscos aos pacientes. Todo o
desenvolvimento é realizado a partir de um ambiente virtual, iniciado com a
modelagem matemática do sistema e utilizando um software de interação
matemática para se obter as respostas do sistema sob diversos aspectos. Os
resultados apresentados sugerem boas perspectivas de sucesso em relação ao
aparelho, alcançando uma redução de 70,8% para um paciente com frequência de
tremor em 7Hz. Portanto, o dispositivo giroscópio surge como mais uma alternativa
de controle dos tremores, de modo que possibilite aos pacientes parkinsonianos a
realização de atividades manuais simples do dia a dia.
Palavras-chave: Doença de Parkinson, Tremores, Absorvedor giroscópio.
88
1 INTRODUÇÃO meio de medicamentos. Isso foi possível

graças ao desenvolvimento de dispositivos
O tremor é classificado como um movimento
capazes de absorver a vibração transmitida
rítmico e involuntário que pode afetar
pelo membro afetado.
qualquer parte do corpo, sendo causado por
fenômenos físicos ou psicológicos (ABBASI; Nos últimos anos diversos trabalhos foram
SAFAIE; AFSHARFARD, 2018). Quando esse desenvolvidos, a exemplo do bracelete
ocorre em pessoas saudáveis, é considerado absorvedor de Buki et al. (2018), os
tremor psicológico, por outro lado, uma absorvedores dinâmicos de vibração
manifestação descontrolada e permanente é desenvolvidos por (GEBAI et al., 2016) e
associada a uma patologia, denominando (LANA et al., 2007). Sendo todos os
assim de tremor patológico (GEBAI et al., dispositivos anteriores elaborados utilizando o
2016). Dentre as doenças mais comuns princípio passivo de absorver o tremor. Por
relacionadas com o tremor, a doença de outro lado, exemplos de atenuadores ativos
Parkinson (DP) e o Tremor Essencial (TE) são estão nos trabalhos de (KAZI; MAILAH; ZAIN,
observadas mais frequentemente em 2014) que utilizam atuadores piezelétricos; a
pacientes. Elas surgem principalmente em colher desenvolvida por (ABBASI; SAFAIE;
pessoas com idade mais avançada, AFSHARFARD, 2018) ao qual consegue ler,
geralmente acima de 65 anos, tendo seus interpretar e atuar para reduzir a vibração
sinais amplificados com o passar dos anos causada pela mão e, por fim, o trabalho
(CICHACZEWSKI et al., 2014). realizado por Taheri, Case e Richer (2011),
que foi desenvolvido com a intenção de criar
Segundo Pinheiro et al. (2017), a DP é a
um algoritmo que possa ser utilizado em uma
segunda doença neurodegenerativa mais
órtese, possibilitando a leitura e supressão do
prevalente no mundo, representando cerca de
tremor.
2% da população mundial acima de 60 anos.
Só no Brasil são 200 mil pessoas, onde 75% Dentro desse contexto surge a proposta deste
delas possuem o tremor como sintoma trabalho, ao qual será simulado a utilização de
principal. um aparelho giroscópico para atenuar o
tremor das mãos. Esse dispositivo se baseia
A DP não se caracteriza como uma doença
na conservação do momento angular, o que
primordialmente ligada ao tremor, porém 80%
permite que ele reaja a qualquer movimento
dos pacientes sofrem com esse sintoma.
que o “perturbe”. Ele é composto
Possui uma faixa de frequência entre 3 a 7 Hz,
basicamente por um disco, que pode ser
podendo chegar a um máximo de 12Hz em
acoplado a um motor elétrico para que possa
alguns casos, podendo aparecer como
fornecer a energia cinética, e uma gaiola que
flexão-extensão ou pronação-supinação no
é responsável por alocar o dispositivo e
antebraço (BUKI et al., 2018).
permitir o seu movimento (HE; LI; CHEN,
Apesar do tremor não ser um sintoma fatal, 2014).
ele é responsável por dificultar bastante a
A reação giroscópica é a responsável pela
vida de qualquer pessoa, uma vez que as
absorção da energia do sistema, uma vez
atividades diárias são realizadas
que, quando o dispositivo está girando em
principalmente com as mãos, seja vestir-se,
altas velocidades e sofre um torque
comer, beber e escrever (ABBASI; SAFAIE;
perturbador, esse reage precessando. A
AFSHARFARD, 2018). Portanto, para amenizar
precessão ocorre quando o eixo de rotação é
essa situação, se faz necessário o tratamento
deslocado do seu eixo fixo do estado inicial,
desse sintoma, que pode ser por medicações
fazendo com que o mesmo gire ao redor
ou cirurgias. As drogas são as mais utilizadas,
desse eixo fixo. Dessa forma, essa reação é a
porém somente cerca de 50% dos pacientes
responsável pelo torque estabilizador, uma
respondem bem a esse método. Outra
vez que ela atua como uma resposta à
maneira de tratar é através de cirurgias,
alteração de seu estado inicial, respeitando,
porém mais indicadas em casos mais severos
portanto, a terceira lei de newton. O torque
da doença, já que é um procedimento
reativo gerado pelo giroscópio será sempre
invasivo que envolvem grandes riscos
igual e oposto ao torque perturbador e foi
(ZHANG et al., 2011).
utilizando esse princípio que surgiram os
Embora a medicina tenha se esforçado primeiros absorvedores giroscópicos
bastante em prol de novas tecnologias, o (SCARBOROUGH, 1958).
tremor patológico se tornou um problema
multidisciplinar, indo além do controle por
89
Dentro desse contexto, este trabalho busca 2 METODOLOGIA

utilizar os princípios giroscópios para
Com o intuito de cumprir de forma metódica
desenvolver uma tecnologia barata e eficaz
todos os objetivos necessários a obtenção
na redução dos tremores causados pela
dos resultados, foi estabelecido o fluxograma
doença de Parkinson, uma vez que ela é a
da Figura 1.
principal patologia neurodegenerativa que
envolve tremuras e responsável pela queda
drástica na qualidade de vida dos pacientes.
Figura 1 – Fluxograma da metodologia de atividades.
NECESSIDADE: Atenuar
INÍCIO DA
tremores no braço de
PESQUISA
pacientes Parkinsonianos
SITES: Google acadêmico,

FASE 1 PROJETO INFORMACIONAL
Periódicos Capes, sciencedirect
NÃO
Informações sobre frequências,
modo de vibrar, tratamentos da
Informações
doença. O comportamento dos
suficientes?
giroscópios e sua relação com a
redução de vibrações
SIM
Software interativo de cálculo

FASE 2 SIMULAÇÃO DO SISTEMA ATENUADOR
numérico (MATLAB)
NÃO
Resposta do sistema atenuador ao

Adequado?
tremor
SIM
Ferramentas gráficas e
FASE 3 ESTABELECER A CAPACIDADE DE ATENUAÇÃO DO SISTEMA
matemáticas
NÃO
Adequados? Avaliação e discussão dos resultados
Redigir texto.
FIM
Fonte: Autoria Própria.
2.1 PROJETO INFORMACIONAL com a doença, os principais métodos de

tratamento, as informações acerca dos
A pesquisa surgiu a partir da necessidade de
giroscópios e a utilização do mesmo como um
novos métodos de controle dos tremores em
atenuador. Sendo todos os pontos citados e
pacientes com Parkinson, dessa forma foi
fundamentados na Seção 1.
realizada a primeira fase dessa pesquisa que
se baseia na obtenção das principais
informações que serão necessárias ao
2.2 SIMULAÇÃO DO SISTEMA ATENUADOR
estudo, permitindo assim, obter os parâmetros
elementares no que diz respeito a doença de Para a realização da simulação e obtenção
Parkinson e ao seu controle. Portanto, consta das respostas do sistema, foi seguido o
nessa primeira fase as informações sobre as fluxograma da Figura 2.
características do tremor, a relação desse
90
Figura 2 – Fluxograma de metodologia de simulação do sistema atenuador.
FASE 2 SIMULAÇÃO DO SISTEMA ATENUADOR
ETAPA 2.1 Modelar matematicamente o sistema atenuador

Artigos periódicos, livros,
equações elementares
ETAPA 2.2 Definir e quantificar os principais parâmetros
NÃO ETAPA 2.3 Obter resposta do sistema principal

Software interativo de cálculo
numérico (MATLAB)
ETAPA 2.4 Obter resposta do sistema atenuador
Resposta do sistema atenuador ao

Adequado?
tremor
2.2.1 MODELAGEM MATEMÁTICA DO toda essa estruturação matemática e física

SISTEMA DE CONTROLO GIROSCÓPICO desenvolvida a partir do livro de Arnold
(1961). Na Figura 3 é possível visualizar o
A modelagem parte do princípio da
diagrama esquemático da configuração.
acoplagem do sistema atenuador ao
segmento representativo do braço, sendo
Figura 3 – Diagrama esquemático do sistema giroscópico atenuador de vibrações
Fonte: Arnold (1961).
Os principais elementos que configuram o 𝜓̇ – Velocidade angular do disco;

sistema são denominados abaixo, onde:
𝜃(𝑡) – Deslocamento angular em torno do eixo
𝐼 – Momento de Inércia do sistema atenuado; X;
𝑘1 – Rigidez do sistema a ser atenuado; 𝜙(𝑡) – Deslocamento angular em torno do
eixo Y;
𝑐1 – Amortecimento do sistema a ser
atenuado; 𝑂𝑥𝑦𝑧 – Conjunto de eixos móveis;
𝑇0 𝑠𝑒𝑛(𝜔𝑡) – Torque desenvolvido pelo sistema 𝑂𝑋𝑌𝑍 – Conjunto de eixos fixos;
a ser atenuado;
Ao ser perturbado pelo torque 𝑇0 𝑠𝑒𝑛(𝜔𝑡)
k𝑔 – Rigidez do sistema giroscópio; deslocando o sistema atenuado em torno do
eixo⁡𝑥, o sistema giroscópio precessa
𝑐𝑔 – Amortecimento do sistema giroscópio;
deslocando-se em torno do eixo 𝑦, sendo
𝐶 – Momento de inércia do disco do sistema esse deslocamento responsável pela criação
giroscópio; do torque atenuador ao qual é igual e oposto
ao torque fornecido pelo braço.
91
Assumindo que os momentos de inércias da velocidades angulares e as quantidades de

gaiola e rotor são definidos por (A,B,C) e movimento sobre os eixos móveis Oxyz são
(A’,B’,C’) respectivamente sobre os eixos dadas por:
móveis Oxyz. As componentes das
𝜔1 = 𝜃̇ 𝑐𝑜𝑠(𝜙)
𝜔2 = 𝜙
𝜔3 = 𝜃̇ 𝑠𝑒𝑛(𝜙)
ℎ1 = (𝐴 + 𝐴′)𝜃̇𝑐𝑜𝑠(𝜙) (1)
ℎ2 = (𝐵 + 𝐵′)𝜙̇ (2)
ℎ3 = (𝐶 + 𝐶 ′ )𝜃̇𝑠𝑒𝑛(𝜙) + 𝐶𝜓⁡ (3)
Substituindo os termos nas equações seja, A = B e n = (𝜓̇ + 𝜃̇𝑠𝑒𝑛𝜙) seja a

diferenciais do movimento sob eixos móveis e velocidade angular total sobre o eixo Oz, tem-
assumindo que o disco giroscópico tenha se então que:
momentos de inércia iguais nos eixos x e y, ou
𝑇1 = (𝐴 + 𝐴′)𝜃̈ + 𝐶𝑛𝜙̇ (4)

𝑇2 = (𝐴 + 𝐵′)∅̈ − 𝐶𝑛𝜃̇ (5)
𝑇3 = 0 (6)
A equação do movimento do sistema sobre o

eixo fixo OX é dada por:
𝐼𝜃̈ = ⁡ 𝑇0 𝑠𝑒𝑛(𝜔𝑡) − ⁡ 𝑇1 𝑐𝑜𝑠(𝜙) − ⁡ 𝑐1 𝜃̇ − ⁡ 𝑘1 𝜃⁡ (7)

Analogamente sobre o eixo Oy, tem-se que:
𝑇2 = −⁡𝑐2 𝜙̇ − ⁡ 𝑘1 𝜙⁡ (8)
Substituindo os termos das Equações (4), (5)

e (6) nas equações (7) e (8), tem-se então:
𝐼1 𝜃̈ + 𝑐1 𝜃̇ + 𝑘1 𝜃 + 𝐶𝑛𝜙̇ = 𝑇0 𝑠𝑒𝑛(𝜔𝑡)
{ (9)
𝐼𝑔 ∅̈ + 𝑐2 𝜙̇ + 𝑘1 𝜙 − 𝐶𝑛𝜃̇ = 0
Onde: 𝐼1 = (𝐼 + ⁡𝐴 + 𝐴′ ); 𝐼𝑔 = (⁡𝐴 + 𝐵′ )
Assumindo a solução da equação (9) sendo longitudinal do sistema a ser atenuado, obtém
da forma 𝜃 = 𝜃1 𝑠𝑒𝑛(𝜔𝑡+∈)⁡em relação ao eixo que a sua solução geral é dada por:
𝑑12 +𝑒12 𝑒1 𝑒2
𝜃(𝑡) = √ 𝑇 sen (𝜔𝑡 + 𝑡𝑔−1 − 𝑡𝑔−1 ) (10)
𝑑22 +𝑒22 0 𝑑1 𝑑2
92
Onde:
𝑑1 = (𝑘𝑔 − 𝐼𝑔 𝜔2 )
𝑑2 = (𝑘1 − 𝐼1 𝜔2 )(𝑘𝑔 − 𝐼𝑔 𝜔2 ) − 𝜔2 (𝑐1 𝑐𝑔 + 𝐶 2 𝑛2 )
𝑒1 = 𝑐𝑔
𝑒2 = 𝑐1 (𝑘𝑔 − 𝐼𝑔 𝜔2 ) + 𝑐𝑔 (𝑘1 − 𝐼1 𝜔2 )
Dessa forma, pode-se assumir uma série de sobre o eixo longitudinal do conjunto a ser
fatores adimensionais a fim de simplificar a atenuado e o segundo correspondente ao
solução da Equação (10), sendo esses sistema giroscópico que se movimenta em
relacionados em dois grupos, o primeiro torno do eixo 𝑦 da Figura 3. O Quadro 1
relativo ao sistema principal que movimentará refere-se aos sistemas citados.
Quadro 1 – Parâmetros adimensionais relativos ao sistema principal e ao sistema giroscópico.

Parâmetros Sistema Principal Sistema Giroscópico Torque aplicado
Amplitude de vibração 𝜃1 ⁄𝜃0 ∅1 ⁄∅0
Fator de amortecimento 𝛽1 𝛽𝑔
Relação de frequência 𝑓 𝑟
Efetividade dinâmica ∆
Onde 𝑓 é a razão de frequência do sistema Utilizando as relações adimensionais do

giroscópico, definida pela relação da Quadro 1 na Equação (10), é possível obter a
frequência natural do giroscópio 𝜔𝑛𝑔 e 𝑟 é a relação de amplitude do sistema que é
razão de frequência da aplicação do torque definida por:
ao sistema, que por definição é considerado a
relação entre a frequência de excitação 𝜔 e a
frequência natural do sistema principal 𝜔𝑛1 .
2 2
𝜃1 (𝑓2 −𝑟 2 ) +(2𝛽𝑔 𝑟)
=√ 2 2 (11)
𝜃0 [(1−𝑟 2 )(𝑓2 −𝑟 2 )−4𝛽1 𝛽𝑔 𝑟 2 −∆𝑟 2 ] +4𝑟 2 [𝛽1 (𝑓2 −𝑟 2 )+𝛽𝑔 (1−𝑟 2 )]
Os termos definidos por 𝜃1 ⁄𝜃0 e ∅1 ⁄∅0 da escolha do tamanho do giroscópio e sua

sugerem as relações de amplitudes tomadas relação com a velocidade angular que o
em relação as condições 𝜃0 , ∅0 . A efetividade mesmo tem. Seu parâmetro é definido por:
dinâmica fornecerá dados importantes acerca
(𝐶𝑛)2
∆⁡= ⁡ , 𝑜𝑛𝑑𝑒⁡𝑝1 ⁡é⁡𝑎⁡𝑓𝑟𝑒𝑞𝑢ê𝑛𝑐𝑖𝑎⁡𝑛𝑎𝑡𝑢𝑟𝑎𝑙⁡𝑑𝑜⁡𝑠𝑖𝑠𝑡𝑒𝑚𝑎⁡𝑝𝑟𝑖𝑛𝑐𝑖𝑝𝑎𝑙. (12)
𝐼1 𝐼𝑔 𝑝12
A Equação (12) retrata os principais Equação (12) por 𝐶 e a velocidade do mesmo

elementos a serem manipulados de tal forma dada por 𝑛.
que mesmo ao alterar qualquer parâmetro,
O sistema a ser atenuado possui uma região
será possível ainda obter a mesma efetividade
em que sua resposta terá um pico máximo,
dinâmica ∆. A ênfase é dada principalmente
conhecido como região de ressonância,
em relação ao tamanho do disco, definida na
93
representado pela linha tracejada 𝑎 no gráfico sistema passe por qualquer um desses novos
da Figura 4. Porém a adição do dispositivo pontos de ressonância. Além do mais, existem
giroscópico acrescentará mais um grau de dois pontos⁡𝑃 e 𝑄 significantes que
liberdade ao sistema e, portanto, sendo independente do amortecimento adicionado
possível observar mais um pico no gráfico esses irão passar. Portanto, é necessário
designado pela linha tracejada 𝑏. Além do encontrar o ajuste ideal de modo que os
mais, a adição de amortecimento ao sistema pontos 𝑃 e 𝑄 no gráfico estejam no mesmo
fará com que a amplitude nesses dois picos nível de amplitude e evitando que
diminua, em contrapartida, a região de independente do sistema passar por qualquer
antirressonância onde antes possuía região de ressonância, seu pico não atingirá
amplitude zero, agora passa a ser maior, níveis elevados. Esse ajuste é encontrado
conforme visto no gráfico através pela linha definindo os valores a partir do fator de
contínua 𝑐. Apesar de se ter uma redução na amortecimento 𝛽𝑔 e a relação de frequência
absorção da vibração, ganha-se bastante em do giroscópio 𝑓.
relação a redução das amplitudes caso o
Figura 4 – Curvas de respostas de sistema com 1GDL sem amortecimento e com 2GDL sem e com
amortecimento.
Fonte: Arnold (1961).
2.2.2 MELHOR RAZÃO DE FREQUÊNCIA DO a Equação (11) de modo a deixá-la na forma

SISTEMA GIROSCÓPIO 𝒇 da Equação (13), assumindo que o termo 𝛽1
seja muito pequeno e que a razão de
Partindo do princípio que é necessário
amplitude independa do amortecimento 𝛽𝑔 .
encontrar a melhor condição em que os
pontos 𝑃 e 𝑄 tenham picos iguais, foi utilizado Nessas condições, é obtido que:
2
𝐸1 +𝐹1 𝛽𝑔
𝜃1
=√ 2 (13)
𝜃0 𝐸2 +𝐹2 𝛽𝑔
Logo, a condição imposta só será satisfeita soluções apresentadas serão as abcissas dos
quando 𝐸1 ⁄𝐸2 + 𝐹1 ⁄𝐹2 , que conduzirá a uma dois pontos 𝑃 e 𝑄 definidas por:
equação do segundo grau em 𝑟 2 cuja as
𝑟 4 − (𝑓 2 + 1 + ∆⁄2)𝑟 2 + 𝑓 2 = 0 (14)
As ordenadas desses pontos são obtidas ao

considerar 𝛽𝑔 = ⁡∞ na Equação (11), logo:
94
1 1
=− (15)
1−𝑟 2 1−𝑟 2
Sendo assim, (𝑟12 + 𝑟22 ) = 2 é a soma das

raízes da Equação (14):
𝑓 2 + 1 + ∆⁄2 = 2 (16)
𝑓 = 𝑓0 = √1 − ∆⁄2 (17)
Esse resultado permite não só definir a razão Uma vez ajustado o absorvedor para a
de frequência do sistema giroscópio como frequência definida pela Equação (17) na
também limitar o dispositivo a uma efetividade Equação (14), tem-se que:
dinâmica ∆< 2.
2.2.3 MÁXIMA RAZÃO DE AMPLITUDE
(𝜽𝟏 ⁄𝜽𝟎 )𝒎 PARA ∆< 𝟐
𝑟 4 − 2𝑟 2 + (1 − ∆⁄2) = 0 (18)
𝑟 2 = 1 + √∆⁄2 (19)
Utilizando a Equação (19) na Equação (11)

quando 𝛽𝑔 = ⁡∞, obtém-se que:
𝜃 1
( 1) = ± = ⁡ √2⁄∆ (20)
𝜃0 𝑚 1−𝑟 2
2.2.4 MELHOR CONDIÇÃO DO FATOR DE sendo esse o valor médio entre os pontos 𝑃 e
AMORTECIMENTO 𝜷𝒈 𝑄. Portanto, ao definir 𝑟𝑚 = 1 e 𝑓 = √1 − ∆⁄2
A melhor condição do fator de amortecimento na Equação (11), é encontrado as seguintes
pode ser obtida através da escolha da relações:
amplitude de √2⁄∆ no ponto referente a 𝑟𝑚 ,
2
(∆⁄2)2 +4𝛽𝑔
2
= (21)
∆ ∆2
1
𝛽𝑔 = 𝛽0 = √∆(8 − ∆) (22)
4
3 RESULTADOS E DISCUSSÕES 3.1 INFORMAÇÕES DO BRAÇO DO

PACIENTE
Os resultados foram obtidos seguindo o
procedimento metodológico apresentado na Utilizando a fundamentação deste artigo, foi
Figura 5, sob as melhores condições elaborarado o Quadro 2 que define os
descritas a partir da modelagem do sistema. principais parâmetros do braço de um
paciente parkinsoniano. Foi estabelecido o
momento de inércia a partir do estudo
antropométrico realizado pela NASA (2010)
para um antebraço e mão esquerda de um
homem possuindo uma frequência
característica do tremor.
95
Quadro 2 – Parâmetros relativo a um antebraço de um paciente parkinsoniano.

Parâmetro Valor
Momento de Inércia 2,43. 10−3 ⁡(𝑘𝑔. 𝑚2 )
Frequência do tremor 7⁡(𝐻𝑧)
3.1.1 RESPOSTA DO MOVIMENTO DO apresentada na Figura 5 ao qual é

BRAÇO DE UM PACIENTE SOB TREMOR DE interpretada como um sistema com um grau
7HZ de liberdade cuja ressonância ocorre no
ponto r=1.
A curva de resposta em frequência do
antebraço do paciente com Parkinson é
Figura 6 – Curva de resposta em frequência do antebraço do paciente com Parkinson.
A razão de frequência relacionada ao torque relevante. Outro ponto importante é a

aplicado 𝑟 mostrada no gráfico refere-se ao quantificação do diâmetro do disco que
ponto em que a frequência de excitação compõe o giroscópio, que não deve ser maior
alcança a frequência natural do braço, ou que 120 mm sendo fabricado em alumínio e
seja, em 7Hz e consequentemente é obtida a pesando 50g para manter as condições de
amplitude máxima de 7,53 para esse ponto. tamanho e peso mínimo aliado a uma boa
capacidade de absorção. Além disso, o motor
que deverá fornecer a energia cinética não
3.2 INFORMAÇÕES QUANTO AO SISTEMA deverá ultrapassar os 5000 rpm de modo que
GIROSCÓPICO ABSORVEDOR sua aquisição não seja dispendiosa. De uma
forma mais simplificada, é expressa na Tabela
É apresentado um dispositivo cujo momento
1 a quantificação dos parâmetros:
de inércia seja 30% do momento de inércia do
antebraço do paciente, uma vez que o peso
do aparelho é um parâmetro bastante
Tabela 1 – Parâmetros relativos ao sistema giroscópio

Momento de inércia do dispositivo giroscópio 7.29.10−4 (𝑘𝑔. 𝑚2 )
Momento de inércia do disco giroscópio 8.10−5 (𝑘𝑔. 𝑚2 )
Rotação do disco 4200 (rpm)
96
Dessa forma é obtido a efetividade dinâmica

do dispositivo é obtido a partir da Equação
(12), em que:
∆⁡= ⁡0,3613
Ou seja, respeitando o pré-requisito da seção A escolha da melhor frequência do sistema

2.2.2 onde sugere o ∆⁡< 2. giroscópio é encontrada a partir da Equação
(17), portanto:
𝑓0 = 0,9052
Dessa maneira, essa é a melhor condição amortecimento, de modo que os pontos 𝑃 e 𝑄

levando em consideração o dimensionamento tenham picos iguais, predefinidos na seção
do dispositivo. 2.2.1, logo:
Por fim, baseando-se na Equação (22), é
possível estabelecer o melhor fator de
𝛽0 = 0,4153
3.2.1 RESPOSTA DO SISTEMA efeitos do fator de amortecimento e para uma

ABSORVEDOR análise mais adequada, é apontado a melhor
curva aplicando os parâmetros calculados na
Com a utilização da Equação (11) é possível
seção 3.2, sendo esses dados apresentados
plotar a curvas das respostas sob diversos
na Figura 7.
Figura 7 - Curvas das respostas sob diversos efeitos do fator de amortecimento.
A partir da análise da Figura 7 é perceptível a especialmente onde a curva apresenta o

redução da amplitude no braço devido ao mínimo local, sendo esse a uma razão de
dispositivo como também a atenuação nos amplitude da ordem de 2,199.
dois picos de ressonância do sistema,
97
Uma análise mais sucinta acerca da frequências de 4,2Hz a 10.5Hz. Portanto,

efetividade do absorvedor é apresentada na qualquer paciente que possuir a frequência
Figura 8 que apresenta na região azul a faixa do tremor dentro dessa região, será
de frequência relativa a absorção do tremor, beneficiado pelo dispositivo. A condição mais
tomando como referência a faixa de tremor da severa do sistema ocorre quando o paciente
doença de Parkinson, sendo essas apresentar uma frequência de 8,89Hz, onde a
apresentadas por 0.6 ≤ 𝑟 ≤ 1,5, ou seja, nas curva tem o seu maior pico de amplitude.
Figura 5 – Região de atuação giroscópica relativo a frequência do tremor de Parkinson.
Para uma análise comparativa, a Tabela 2

fornece os dados sobre o quanto o dispositivo
é efetivo:
Tabela 2 – Tabela comparativa da efetividade do giroscópio.

Sem o Absorvedor Giroscópico Com o Absorvedor Giroscópico Redução
Pico 7,53 2,199 70,8%
Figura 9 - Gráfico de comparação do sistema sem o dispositivo e com o dispositivo atuando.
Considerando a Figura 9, é fácil observar a sistema com dois graus de liberdade. A

capacidade de absorção da amplitude que o Tabela 3 define a amplificação nesses dois
dispositivo consegue alcançar mesmo nos pontos tomando como referência o valor
dois pontos de ressonância comum em um mínimo de amplitude:
98
Tabela 3 – Tabela de análise da amplificação do sistema.

Parâmetro Valor Amplificação
Pico de referência 2.199 -
Pico 1 2,467 12,19%
Pico 2 2,93 33,24%
Percebe-se então a grande importância do giroscópio é bastante simples de ser
sistema ajustado segundo os parâmetros de construído além de permitir uma variedade de
razão de frequência e amortecimento do desenhos para a sua construção.
sistema giroscópico, já que garante uma boa
Embora o desenvolvimento do dispositivo
absorção de energia pelo dispositivo dentro
tenha sido realizado a partir de um ambiente
de uma larga faixa de frequência.
simulado, os elementos que compõe o
dispositivo são bastante simples de ser
adquirido, além de fácil construção. Porém,
4 CONCLUSÕES
fatores como elevados ruídos devido ao motor
A partir do estudo baseado na simulação do elétrico que fornece energia cinética ao disco.
dispositivo, foram obtidos os resultados
De uma forma geral, a proposta desse
relativos a sua capacidade de absorção de
dispositivo absorvedor giroscópio é uma
amplitude do tremor. A análise permite
tecnologia relativamente nova, mas altamente
observar bons resultados, chegando na casa
eficiente na redução de vibração, por isso
dos 70% de redução para um paciente com
surge como mais uma alternativa capaz de
tremor a uma frequência de 7Hz. A grande
oferecer melhoria na qualidade de vida dos
vantagem em relação a outros dispositivos é o
pacientes que possuem a Doença de
seu custo, uma vez que o dispositivo
Parkinson.
REFERÊNCIAS
[1] ABBASI, M.; SAFAIE, J.; AFSHARFARD, A. [9] KAZI, S.; MAILAH, M.; ZAIN, Z.
Design of a Noninvasive Hand Tremor Suppression Suppression of Hand Postural Tremor Via Active
System with Active Control. p. 1–7, 2018. Force Control Method 3 Dynamic Model of Human
Hand 2 Background of Hand Tremor. p. 76–82,
[2] BUKI, E. et al. Vib-bracelet: a passive
2014.
absorber for attenuating forearm tremor. Medical
and Biological Engineering and Computing, v. 56, [10] LANA, R. C. et al. Percepção da qualidade
n. 5, p. 923–930, 2018. de vida de indivíduos com Doença de Parkinson
através do PDQ-39. Revista Brasileira de
[3] CARPENTER, M. D.; PECK, M. A.
Fisioterapia, v. 11, n. 5, p. 397–402, 2007.
Reducing Base Reactions With Gyroscopic
Actuation of Space-Robotic Systems. v. 25, n. 6, p. [11] NASA Handbook (2010) NASA Human
1262–1270, 2009. Integration Design Handbook (HIDH)-NASA, Vol.
3407. SP-2010.
[4] CICHACZEWSKI, E. et al. Características
eletrofisiológicas do tremor na doença de [12] PINHEIRO, W. C. et al. Parkinson’s
parkinson e tremor essencial. Arquivos de Neuro- Disease Tremor Suppression A Double Approach
Psiquiatria, v. 72, n. 4, p. 301–306, 2014. Study - Part 1. Vol 1: Biodevices, n. Biostec, p.
149–155, 2017.
[5] GAGNE, J. et al. Active heart stabilization
using adaptive noise cancelling techniques with [13] SCARBOROUGH, J. B. The gyroscope;
gyroscopic actuation. 2010, BioRob 2010, p. 802– Theory and applications, 1958.
807, 2010.
[14] TAHERI, B.; CASE, D.; RICHER, E. Active
[6] GEBAI, S. et al. Tremor Reduction at the Tremor Estimation and Suppression in Human
Palm of a Parkinson’s Patient Using Dynamic Elbow Joint. ASME 2011, Volume 2, p. 115–120,
Vibration Absorber. Bioengineering, v. 3, n. 3, p. 2011.
18, 2016.
[15] TOWNSEND, N. C.; SHENOI, R. A. Control
[7] HE, H.; LI, R.; CHEN, K. High-rise Structure strategies for marine gyrostabilizers. IEEE Journal
Vibration Control with Gyroscopes. n. Icmce, p. of Oceanic Engineering, v. 39, n. 2, p. 243–255,
227–232, 2014. 2014.
[8] KARNOPP, D. Tilt Control for Gyro- [16] ZHANG, D. et al. Neural oscillator based
Stabilized Two-Wheeled Vehicles. Vehicle System control for pathological tremor suppression via
Dynamics, v. 37, n. 2, p. 145–156, 2002. functional electrical stimulation. Control Engineering
Practice, v. 19, n. 1, p. 74–88, 2011.
99
Capítulo 10

Mikarla Baía de Sousa
Josiane Dantas Costa
Arthur Filgueira de Almeida
Renato Alexandre Costa de Santana.
Resumo: A eletrodeposição é considerada uma técnica rápida, fácil e econômica

para a obtenção de revestimentos metálicos. Uma ampla variedade de
propriedades para revestimentos pode ser obtida selecionando diferentes
parâmetros de eletrodeposição. Os revestimentos metálicos compostos de cromo,
cobre, tungstênio, zinco, cobalto e molibdênio são empregados em muitos campos,
como corrosão, desgaste, aplicações magnéticas e elétricas e materiais
eletrocatalíticos. Neste estudo, apresenta-se uma revisão da literatura de alguns
dos metais mais utilizados como revestimentos, sendo de grande valia para guiar
futuras pesquisas. O processo de eletrodeposição, parâmetros efetivos,
propriedades e aplicação de tais ligas foram resumidas considerando os estudos
anteriores.
Palavras-Chave: Eletrodeposição. Ligas Metálicas. Revestimentos Metálicos.

Corrosão.
100
1 INTRODUÇÃO utilizado para a proteção dos metais (NOCE et

al., 2014).
Devido à variedade de materiais (metais, ligas
metálicas, polímeros, etc.) o papel de escolha A eletrodeposição é o processo pelo qual um
destes, por um profissional, para uso em revestimento metálico é aplicado sobre uma
estruturas ou equipamentos se torna difícil. O superfície através de uma corrente elétrica,
material ideal terá que apresentar três geralmente contínua. A deposição pode ser
características importantes: propriedades feita em superfície condutora (metais e ligas
desejadas, menor custo possível e metálicas: aço carbono, latão) ou não
durabilidade. Porém, quando se relaciona a condutora (plásticos, couro, madeira). Através
durabilidade do material com a sua da eletrodeposição é possível obter materiais
resistência à corrosão, prever como será seu metálicos com propriedades satisfatórias
comportamento se torna trabalhoso. Isto como dureza, resistência à corrosão,
porque a natureza do meio corrosivo e as propriedades catalíticas, boa ductilidade, etc.
condições de exposição do material a este
A eletrodeposição possibilita o controle de
meio têm que ser consideradas (AZEVEDO et
parâmetros experimentais, como composição
al., 2015).
da solução eletrolítica, temperatura de
Neste aspecto, os problemas causados pela deposição, densidade de corrente aplicada,
corrosão de materiais se tornam frequentes carga a ser depositada e pH da solução. O
em diversos setores da atividade humana controle destes parâmetros favorece o estudo
(construção civil, automobilística, ferroviária, das características físicas e químicas dos
marítima, odontológica, indústria química, revestimentos (WHITE; BOCKRIS; CONWAY,
etc.) proporcionando grandes gastos 1995).
econômicos diretos (substituição e prevenção
Nos últimos anos o interesse pela
dos materiais) e indiretos (consequência da
eletrodeposição de ligas tem aumentado
corrosão dos materiais). Destas perdas, as
devido as suas propriedades, tais como alta
mais prejudiciais são as indiretas, pois são de
condutividade elétrica, resistência mecânica,
difícil avaliação e na maioria das vezes de
alta condutividade térmica e boas
custos bem mais elevados. Algumas dessas
propriedades catalíticas (SILVA, 2006). Essas
perdas indiretas são: paralisações acidentais,
ligas são aplicadas nas mais diversas áreas,
perdas de produtos, contaminações de
como por exemplo, nas indústrias químicas,
produtos, etc. (GENTIL, 2007).
petrolíferas, navais, de construções civis e
Deve-se considerar a importância social com automobilísticas em decorrência de algumas
relação a fatores de segurança à saúde características especiais, como alta
humana, em que acidentes devido à corrosão resistência à corrosão e desgaste.
são imprevisíveis. Mediante estes problemas
que abrangem diversas áreas, causando
gastos econômicos e riscos de vida, o estudo 2 CONCEITOS E APLICAÇÕES
da corrosão, no que se refere à melhor forma
A eletrodeposição de metais e ligas permite a
de evitá-la ou até mesmo minimizá-la, tem
obtenção de materiais que podem ser
crescido e se mostrado importante do ponto
utilizados numa variedade de aplicações
de vista tecnológico, econômico e social na
(HAMID e HASSAN, 2012), desde que os
preservação da natureza e do homem.
depósitos obtidos apresentem propriedades
Uma forma de se minimizar a corrosão em convenientes, tais como: resistência à
materiais metálicos é revesti-los com outros corrosão, resistência ao desgaste,
materiais que apresentam características mais propriedades magnéticas e ópticas
nobres. Dessa forma, o revestimento metálico desejadas, etc. Todas estas características
é uma das formas de pretensão anticorrosiva. dependem da composição química e da
Estes revestimentos podem ser obtidos microestrutura obtida em revestimentos e/ou
através de deposição em fase gasosa, da topografia da superfície do metal tratado.
eletrodeposição, metalização, imersão a Deste modo, os processos de
quente, redução química, etc. Em eletrodeposição vêm sendo intensivamente
consequência do desenvolvimento dos estudados com o objetivo de controlar a
geradores de corrente contínua e por ser microestrutura e composição dos depósitos,
economicamente viável, a prática da visando à otimização de propriedades que
eletrodeposição tornou-se o processo mais melhorem seu desempenho para as diferentes
aplicações.
101
A eletrodeposição é o método mais submersos em um eletrólito contendo o banho

importante para a produção comercial de de deposição (GAMBURG e ZANGARI, 2011;
revestimentos protetores. É um processo ZOLDAN, 2007).
comumente utilizado, pois se consegue
As correntes dentro do eletrólito provêm da
revestimento muito fino e relativamente livre
movimentação das espécies iônicas, que
de poros, evitando-se excesso do metal
pode ocorrer via migração, convecção e
eletrodepositado. A partir do momento em
difusão, e subsequente oxidação/redução nos
que a técnica da eletrodeposição passou a
eletrodos. A migração é um mecanismo de
ser utilizada na produção de revestimentos
movimentação iônica originada a partir da
finos e multicamadas, ela passou a ter um
influência que o campo elétrico gerado pela
papel de destaque na indústria eletrônica com
diferença de potencial existente entre os
as deposições de materiais condutores ou
eletrodos submersos na solução causa nas
camadas magnéticas. Além do mais, a
partículas carregadas. A convecção é o
eletrodeposição se apresenta como sendo
transporte hidrodinâmico que ocorre devido à
uma técnica de baixo custo de implantação e
diferença nas densidades naturais de
de produção (BELTOWSKA-LEHMAN, 2012).
concentração causadas por flutuações na
As coberturas por eletrodeposição são de
temperatura e na concentração da solução. A
espessuras mais uniformes do que as obtidas
difusão ocorre devido à diferença de
por imersão a quente, menos porosas do que
concentração da espécie a ser depositada
as vaporizadas e de maior pureza. O método
em função da posição (gradiente). Para
permite facilmente o controle da espessura da
eletrólitos não muito concentrados, o
camada depositada.
crescimento do depósito na superfície do
Devido à sua facilidade de aplicação em eletrodo é controlado pelo transporte de
substratos com as mais variadas formas e massa difusivo (regime difusivo). Para
tamanhos possíveis, este processo é eletrólitos muito concentrados, o crescimento
adequado quando se tem peças da camada depositada é controlado pela taxa
relativamente pequenas como parafusos e das reações de oxirredução (regime cinético)
rebites, chapas e fios produzidos em (ZOLDAN, 2007).
processos contínuos, para uso em ambiente
De acordo com a lei de Faraday para a
interno e em ambientes externos pouco
eletrólise, a quantidade de todos os
poluídos. Para que seja possível obter
elementos liberados no cátodo ou no ânodo
eletrodepósitos com as características físicas
durante a eletrólise é proporcional à
e químicas desejadas é necessário
quantidade de eletricidade que se passa
inicialmente entender o processo de
através da solução. Desde que no cátodo íons
eletrodeposição; portanto, é importante
hidrogênio sejam descarregados ao mesmo
lembrar que um grande número de variáveis
tempo em que os íons metálicos, a eficiência
deve ser considerado no estudo de
do processo pode ser, e geralmente é menor
eletrodeposição. Isto faz com que diferentes
que 100%. De qualquer modo, sob quaisquer
técnicas eletroquímicas possam ser
condições, a espessura do depósito é
empregadas nestes estudos e que haja uma
proporcional ao tempo de deposição.
grande inter-relação de diferentes áreas para
estudar a cinética de formação de O ânodo desta célula é, geralmente, do
recobrimentos metálicos e em seguida mesmo metal a ser depositado, mas pode ser
caracterizar o material formado. também de algum material altamente
condutor (inerte como a platina), insolúvel no
O processo de eletrodeposição é,
banho nas condições atuantes. Dessa forma,
essencialmente, uma reação eletroquímica
devido à diferença de potencial, os íons
que ocorre na área catódica da célula
metálicos entram na solução por dissolução,
eletroquímica para produzir depósitos sólidos
movem-se na direção do cátodo, onde são
na superfície de substratos (TOMACHUK,
depositados na forma metálica (MOURA,
1998). A formação dos depósitos está
2005).
associada com reações de oxirredução de
espécies químicas (na reação de oxidação a
espécie perde elétrons, tornando-se mais
3 ELETRODEPOSIÇÃO DE LIGAS
positiva, e na reação de redução recebe
elétrons, tornando-se mais negativa). Tais No caso de eletrodeposição de ligas deve-se
reações são produzidas pela passagem de analisar a possibilidade da codeposição dos
uma corrente elétrica entre dois eletrodos metais. A codeposição de dois metais A e B
102
ocorre a velocidades apreciáveis quando em zinco e cobre-estanho em soluções de

uma determinada densidade de corrente os cianeto.
potenciais de redução de dois metais se
igualam. A pré-condição para formação de  Codeposição de equilíbrio:
qualquer liga é a codeposição dos metais em caracteriza-se pela deposição a partir de uma
velocidades apreciáveis. A modificação de solução que está em equilíbrio químico com
parâmetros como densidade de corrente, a os metais eletrodepositados. Este sistema é o
concentração de íons metálicos e a alteração único em que a razão molar dos íons
da composição das soluções fazem com que metálicos no banho é igual à razão dos metais
seja induzida a codeposição, mesmo no caso no eletrodepósitos. Poucas ligas que
em que os potenciais reversíveis sejam pertencem a este grupo são conhecidas.
diferentes. Observa-se então que uma das Dentre elas pode-se citar: cobre-níquel,
primeiras etapas na investigação do processo eletrodepositada a partir de banhos de
de eletrodeposição de ligas é desenvolver o tiossulfato; estanho-chumbo e cobre-bismuto
banho de eletrodeposição, ou seja, as eletrodepositadas a partir de banhos ácidos.
soluções que serão utilizadas no processo.  Codeposição anômala: é o tipo de
Isto requer um conhecimento prático da deposição onde o metal menos nobre é
eletroquímica dos elementos que constituirão depositado preferencialmente. Esta
a liga, da solubilidade de seus sais e da característica ocorre somente sob condições
química de seus complexos (SANTANA et al., bem determinadas de concentração e
2003). variáveis de operação. Pode ocorrer com íons
Brenner (1963) classificou as ligas livres ou complexados e está frequentemente
eletrodepositadas em cinco classes distintas associada à eletrodeposição de ligas
de codeposição, levando em consideração a contendo um dos três metais do grupo do
composição do depósito, sendo estas: ferro (cobalto, níquel e ferro).
 Codeposição regular: ocorre quando  Codeposição induzida: caracteriza-se

o processo está sendo controlado por difusão. pela deposição de ligas contendo metais
Os efeitos das variáveis do processo na como molibdênio, tungstênio ou germânio.
composição dos depósitos são determinados Estes metais não podem ser
pelas mudanças na concentração dos íons eletrodepositados na forma metálica
metálicos na camada de difusão catódica e individualmente, mas se depositam formando
são, consequentemente, previstos pela teoria ligas com metais do grupo do ferro.
de difusão. A porcentagem do metal mais Os três primeiros são denominados de
nobre eletrodepositado é aumentada pela codeposição normal, pois se caracterizam por
diminuição da densidade de corrente, apresentarem a deposição preferencial do
elevação da temperatura, aumento da metal mais nobre. Os dois últimos são
concentração total do metal e agitação. Este denominados de codeposição anormal
tipo de codeposição ocorre geralmente em porque há deposição preferencial do metal
banhos contendo íons metálicos livres, mas menos nobre.
pode ocorrer também com íons metálicos
complexados. Nesta classe podem-se Em 1994, Landolt publicou um artigo
encontrar ligas de estanho-chumbo, bismuto- propondo três principais mecanismos de
cobre etc. eletrodeposição de ligas que explicam os
principais tipos de deposição propostos por
 Codeposição irregular: caracteriza-se Brenner. Landolt na década de 90 se
por ser controlada em maior extensão pelo mostrava envolvido em uma série de estudos
potencial de ativação do que pelo fenômeno relacionados a processos de eletrodeposição
de difusão. Os efeitos de algumas variáveis de ligas, assim como Brenner na década de
do processo influenciam na composição dos 60.
eletrodepósitos de acordo com a teoria de
difusão, entretanto outras variáveis podem ter A partir deste trabalho os sistemas de
comportamento contrário a esta teoria. eletrodeposição passaram a ser classificados
Também se observa um efeito menor das como:
variáveis do processo na composição dos  Sistemas acoplados por transporte de
eletrodepósitos em relação a um sistema massa (onde a corrente parcial de um metal
regular. A codeposição ocorre geralmente em depende do transporte de massa do outro);
soluções contendo íons complexos. Nesta
classe, podem ser incluídas ligas de cobre-
103
 Sistemas acoplados por transferência depósitos. Estes parâmetros agem direta ou

de carga (onde as correntes de deposição de indiretamente no processo de
cada metal são dependentes uma da outra); eletrodeposição, tendo influência nas
propriedades morfológicas e estruturais dos
 Sistemas não interativos (onde as depósitos. Com base no exposto acima,
correntes de deposição de cada metal são pode-se concluir que estudar o processo de
independentes uma da outra). eletrodeposição de uma determinada liga
A teoria envolvida nestes três processos está envolve, além de um conhecimento de
baseada na teoria de potencial misto e nos diversas áreas, um critério para escolha do
processos que envolvem a reação de material que irá compor esta liga.
eletrodeposição, como o fenômeno de
transferência de elétrons e a difusão dos íons
até a superfície do eletrodo. Há necessidade 4 ALGUNS TIPOS DE LIGAS
em compreender estes fenômenos tanto do Nos últimos anos várias ligas obtidas por
lado da solução quanto da parte da região eletrodeposição vêm sendo estudadas para
interfacial no sólido, para que se possa obter serem aplicadas como revestimentos e assim
eletrodepósitos com as características diminuírem os problemas causados pela
desejadas. Por razões ligadas com a corrosão. A seguir algumas dessas ligas
evolução histórica dos modelos teóricos, a serão apresentadas.
região interfacial na solução é conhecida
como a região da dupla camada elétrica.
Nesta região, o valor do potencial eletrostático 4.1 LIGAS DE CROMO
difere do potencial dos íons no interior da
solução. A espessura desta camada é da O cromo (Cr) possui um alto ponto de fusão
ordem de nanômetros, e é nela que ocorrem de 2.136 K, alta resistência mecânica e é
as reações de transferência de carga, quimicamente estável sob várias condições.
responsáveis pela formação dos depósitos Devido ao alto ponto de fusão seu processo
metálicos. de obtenção à temperatura ambiente torna-se
difícil, e desta forma, as ligas de cromo são
produzidas geralmente por um processo
3.1 FATORES QUE AFETAM O PROCESSO pirometalúrgico (processo sob a condição de
DE ELETRODEPOSIÇÃO alta temperatura e vácuo em uma câmara
fechada). (OHGAI; TANAKA; FUJIMARU,
Após a elaboração de um banho de 2012).
eletrodeposição conveniente, a próxima etapa
a ser investigada é o efeito de todas as Uma solução viável para obtenção dos
variáveis que influenciam na composição dos revestimentos de cromo à temperatura
eletrodepósitos. Estas variáveis são de dois ambiente e pressão atmosférica se dá através
tipos: i) variáveis de composição, tais como a do processo convencional de
natureza dos complexos, aditivos, pH, eletrodeposição. Onde filmes metálicos de Cr
concentração e razão molar entre os íons podem ser sintetizados por eletrodeposição a
metálicos e ii) variáveis de operação, como partir de uma solução aquosa contendo íons
densidade de corrente, temperatura e de cromo hexavalente (Cr6+) ou íons de
agitação. Estas variáveis exercem influência cromo trivalente (Cr3+) (OHGAI; TANAKA;
direta sobre o caráter do depósito e estão FUJIMARU, 2012). Considerações ambientais
relacionadas entre si de tal forma que se uma têm sido responsáveis pelo maior interesse na
delas for alterada, haverá alteração nas outras eletrodeposição de cromo trivalente que é um
(HALIM et al., 2012). processo menos tóxico. Então, nas últimas
décadas, muitas tentativas foram feitas para
A preparação da superfície dos substratos produzir ligas de cromo depositadas a partir
sobre os quais se deseja efetuar a de eletrólitos trivalentes de cromo como novos
eletrodeposição é extremamente importante materiais para aplicação, a fim de substituir o
para a obtenção de coberturas aderentes, processo de galvanoplastia de cromo
compactas e homogêneas. Inicialmente é hexavalente. Como exemplo, podemos citar
fundamental para o êxito da eletrodeposição as ligas Fe-Cr, Cr-Ni, Fe-Cr-Ni, Cr-P e Cr-C (LI
que a peça esteja perfeitamente limpa. Muitos et al., 2008).
pesquisadores acrescentam aditivos
orgânicos na solução para o controle das O grande interesse na obtenção dos
propriedades morfológicas e mecânicas dos revestimentos de cromo se deve a sua alta
104
resistência à corrosão devido a sua necessidade de ligas com características

propriedade de passivação em alguns especiais, novos desenvolvimentos foram
ambientes corrosivos, como ácidos sulfúrico e obtidos visando o mercado cativo da indústria
nítrico, soluções alcalinas e outros. Nestes microeletrônica e de telecomunicações.
ambientes o filme fino que se forma na
A eletrodeposição de cobre e suas ligas a
superfície do cromo impede a dissolução do
partir de soluções alcalinas cianetadas é
metal. No entanto, a desvantagem
amplamente empregada na indústria para
significativa dos revestimentos de cromo é a
obtenção de camadas de cobre de alta
sua instabilidade na presença de íons cloreto.
qualidade sobre substratos ferrosos e não
Isso pode ser superado através do
ferrosos. O interesse por novos processos de
desenvolvimento de ligas com pequenas
eletrodeposição de cobre isento de íon
quantidades de fósforo, molibdênio,
cianeto tem mostrado um aumento nas
tungstênio e nióbio (KUZNETSOV et al., 2015).
pesquisas desde 1980, devido à implantação
As ligas de Co-Cr chamam a atenção dos de novas normas ambientais e imposições de
pesquisadores devido à boas propriedades novas condições de segurança, até os dias
mecânicas e alta resistência à corrosão. Em atuais.
particular, as ligas Co-Cr-W e Co-Cr-Mo
Ligas de Cu-Zn têm sido amplamente
contendo uma pequena quantidade de
aplicadas devido à sua elevada resistência
carbono, as chamadas ligas “Stellite”, tendo
mecânica, boa plasticidade e excelente
uso extensivo tanto na indústria quanto em
resistência à corrosão (ALMEIDA, DE et al.,
aplicações biomédicas. Isto é devido às suas
2011). Por exemplo, as ligas de Cu-Zn têm
boas propriedades mecânicas, resistência à
sido amplamente utilizadas como material
corrosão e biocompatibilidade. A inibição do
elétrico em dispositivos microeletrônicos,
processo de corrosão está relacionado à
conectores automotriz, máquinas elétricas e
formação de um filme constituído por óxido e
instrumentos eletrônicos, etc. Os materiais
hidróxido de cromo na superfície da liga sob
utilizados no contato elétrico são geralmente
exposição a meios corrosivos (KUZNETSOV et
constituídos de ligas de cobre-zinco
al., 2018).
interessantes por sua excelente condutividade
Outro grupo de resvestimentos de cromo que elétrica e um baixo custo.
vem ganhando destaque são os
Os revestimentos eletrodepositados de Cu-Zn
revestimentos amorfos obtidos por
são amplamente utilizados para fins
eletrodeposição, pois estes têm maior
decorativos, bem como para promover a sua
capacidade de passivar espontaneamente e
boa aderência à borracha, que é empregue
maior resistência à corrosão em comparação
na indústria de pneus (VREESE, DE et al.,
aos revestimentos cristalinos. Dentre os
2013). A eletrodeposição comercial de ligas
revestimentos de cromo amorfos os de Fe-Ni-
de Cu-Zn em banhos de cianeto produz
Cr são os mais conhecidos, estes
depósitos de alta qualidade, mas os
revestimentos são amplamente utilizados para
problemas ambientais surgem na utilização e
aplicações aeronáuticas devido à suas
eliminação de cianeto. Assim, uma extensa
propriedades desejáveis de resistência
pesquisa tem sido feita para encontrar
mecânica, resistência à corrosão e efeitos
eletrólitos alternativos satisfatórios (SENNA;
magnéticos (TAVOOSI; BARAHIMI, 2017).
DÍAZ; SATHLER, 2005; CARLOS; ALMEIDA,
DE, 2004; KARAHAN; ÖZDEMIR 2014).
4.2 LIGAS DE COBRE As ligas de Cu-Ni são amplamente utilizadas
na mineração e obras metalúrgicas, na
O cobre é um metal mole, maleável, dúctil,
indústria química devido à sua elevada
tenaz com alta condutividade elétrica, alta
resistência à corrosão, propriedades
condutividade térmica, boa resistência
mecânica e catalítica (WANG et al., 2014;
mecânica, boa trabalhabilidade e boa
GORANOVA; AVDEEV; RASHKOV, 2014).
resistência à corrosão (SILVA, 2006). Por
Além disso, as ligas de Cu-Ni têm recebido
esses motivos, o cobre tem sido cada vez
muita atenção por suas propriedades
mais utilizado em escala industrial e,
magnéticas e termofísicas e também são
principalmente na indústria de
usadas para fins decorativos (BASKARAN et
eletroeletrônicos. O cobre, como elemento
al., 2006). São aplicadas em válvulas,
principal, com a adição de outros elementos,
ferragens, condensadores e trocadores de
apresenta-se como liga para fins específicos
calor, pela sua boa usinabilidade, excelentes
(SPIEGELBERG, 1990). Devido à crescente
105
propriedades térmicas e elétricas e são vidro “Pyrex” (SANTANA; PRASAD; DE

particularmente resistentes à SANTANA, 2003).
biocontaminação (EL-CHIEKH et al., 2005).
Diante do exposto o interesse no
Nos últimos anos, o interesse pela desenvolvimento das ligas de tungstênio
eletrodeposição da liga Cu-Ni em diferentes obtidas por eletrodeposição vem crescendo
composições de banho eletrolítico (variando nos últimos anos. Os revestimentos de
de Cu-05Ni até Cu-64Ni) tem aumentado, por tungstênio possuem propriedades excelentes
suas atribuições mecânicas, tais como força como: alta estabilidade térmica, elevada
de tração, maleabilidade, ductilidade e dureza, resistência à corrosão e ao desgaste
propriedades magnéticas únicas (Ni é (TSYNTSARU et al., 2012), tornando-os ótimos
ferromagnético e Cu é diamagnético). A liga substitutos do cromo duro (revestimento que
Cu-Ni, principalmente as que contêm cerca apresenta o cromo hexavalente em seus
de 30% a 40% de Cu, são resistentes à banhos eletrolíticos, e esse por sua vez é
corrosão em meio ácido, meio básico e tóxico e tem restrições ambientais)
especialmente em soluções contendo cloreto (BROOMAN, 2004; CAPEL; SHIPWAY;
(água do mar). Esta liga tem sido utilizada na HARRIS, 2003; OLIVEIRA et al., 2015).
construção de embarcações marítimas e em
A eletrodeposição direta do tungstênio a partir
trocadores de calor que usem a água do mar
de soluções de tungstato é impedida pela
para refrigeração. Outro aspecto que tem
formação de uma camada de óxido sobre o
motivado o estudo da eletrodeposição de Cu-
cátodo e neste caso toda a corrente utilizada
Ni é sua aplicação na forma de pó em tintas,
no processo é consumida para evolução do
utilizando banhos amoniacais, também
hidrogênio (LASSNER; SCHUBERT, 1999).
empregadas contra a agressividade da
Desta forma, o tungstênio não pode ser
atmosfera marinha (WANG et al., 2014).
eletrodepositado sozinho e sua redução só
acontece na presença de um metal indutor,
este fenômeno é conhecido como "Co-
4.3 LIGAS DE TUNGSTÊNIO
deposição induzida” e foi relatado por
O tungstênio possui um considerável Brenner (BRENNER, 1963). A partir deste
interesse para ser eletrodepositado em forma conhecimento alguns pesquisadores
de ligas, pois este metal apresenta desenvolveram banhos eletrolíticos contendo
propriedades bastante interessantes como: o outros metais e verificaram que a
mais alto ponto de fusão (3.410 °C), o que eletrodeposição do tungstênio apresentou
torna impossível depositálo termicamente sucesso na presença dos metais do grupo do
sobre a superfície de qualquer outro metal, ferro como Fe, Ni, Co e etc.
uma vez que este último muda de fase; o mais
O interesse em ligas binárias e ternárias ricas
baixo coeficiente de dilatação linear térmica
em tungstênio eletrodepositadas com metais
(4,3x10-6 °C- ¹); uma condutividade térmica
do grupo de ferro tem aumentado nos últimos
muito alta (0,487 cal/cm°C) e é um dos metais
anos devido à sua combinação única de
mais densos (19,3 g/cm³). Tem propriedades
propriedades tribológicas, magnéticas,
mecânicas fora do comum, pois conserva sua
elétricas e de eletro-erosão (CHEN et al.,
dureza mesmo em altas temperaturas, tem a
2016; ĆIROVIĆ et al., 2016; MUN et al., 2010;
mais alta resistência à tração (410 kg/mm²) e
TSYNTSARU et al., 2009). Suas aplicações
é muito resistente à corrosão (SANTANA;
são bem documentadas e exemplos podem
PRASAD; DE SANTANA, 2003).
ser encontrados em proteção contra corrosão,
Além disso, o tungstênio não é atacado por eletrônica para camadas de barreira e
nenhum ácido mineral simples na temperatura catálise. A introdução de um terceiro elemento
ambiente. Devido a esta insólita combinação em ligas binárias de tungstênio, em alguns
de propriedades, este metal tem inúmeras casos, pode melhorar as propriedades de
aplicações na indústria e na engenharia, não corrosão, por exemplo, a liga Ni-W-Co teve
obstante o seu elevado preço e a sua relativa melhor resistência à corrosão do que as ligas
raridade. Por exemplo, suas boas amorfas Co-W similares a eletrodos. Além
propriedades mecânicas em elevadas disso, eletrocatalisadores eficientes e de
temperaturas o tornam útil para ser usado em baixo custo são importantes para eletrodos de
filamento de lâmpadas incandescentes, células de combustível e ligas de tungstênio
bocais de motores de foguetes e em aços podem ser uma alternativa possível à platina e
para ferramentas. O seu baixo coeficiente de às ligas caras que contêm metais do grupo da
expansão térmica permite a sua inserção no
106
platina usados atualmente (CESIULIS; XIE; apresentam resistência à corrosão superior

PODLAHA-MURPHY, 2009). àquelas com teores de níquel mais baixos. A
composição do depósito é dependente do
Algumas ligas binárias de tungstênio (Ni-W e
teor de Ni2+ no eletrólito (GNANAMUTHU et
Fe-W) podem apresentar fissuras causadas
al., 2012; GAVRILA et al., 2000; HAMMAMI;
por tensões internas quando o conteúdo de
DHOUIBI; TRIKI, 2009). Em relação ao zinco
tungstênio é demasiadamente alto ou quando
convencional, as ligas Zn-Ni apresentam
a densidade de corrente é alta. Muitos
dureza significativamente maior e melhor
pesquisadores acreditam que a liga ternária
resistência à corrosão, principalmente depois
Ni-Fe-W pode exibir as propriedades chave
de tratamento térmico, o que confere
desses depósitos binários, eliminando as
propriedades adequadas para utilização em
propriedades indesejáveis (defeitos e tensões
componentes de motores dos automóveis
internas, por exemplo). Principalmente se a
(CONDE; ARENAS; DE DAMBORENEA, 2011).
liga ternária possuir estrutura amorfa, pois
esse tipo de estrutura apresenta A eletrodeposição de ligas Zn-Co é de grande
homogeneidade e falta de sítios interesse uma vez que estas ligas apresentam
eletroquimicamente ativos que geram resistência à corrosão superior à do zinco
revestimentos com baixa taxa de corrosão puro (BAJAT et al., 2010). Revestimentos com
(COSTA et al., 2018). baixos teores de cobalto são menos nobres
que o aço e, portanto, apresentam proteção
sacrificial a substratos de aço. Os que
4.4 LIGAS DE ZINCO possuem alto teor de Co são mais nobres que
o aço e proporcionam proteção por efeito
O zinco é um metal pouco nobre que vem
barreira. Comparado com o zinco puro, as
sendo utilizado há mais de 170 anos como
camadas de ligas zinco-cobalto possuem
revestimento de sacrifício. A eletrodeposição
outras propriedades protetoras como dureza,
de zinco é um dos processos de aplicação
ductilidade, tensão interna, soldabilidade e
deste tipo de revestimento sobre superfícies
habilidade para receber revestimento de
de aço-carbono mais utilizados
tintas.
industrialmente, isto se deve principalmente
ao baixo custo do processo. As propriedades O efeito de proteção das ligas Zn-Mn é
das camadas de zinco obtidas são dependente do teor de Mn na liga. A elevada
dependentes das características dos resistência à corrosão das ligas Zn-Mn tem
parâmetros do processo de eletrodeposição. sido atribuída ao efeito duplo de proteção
observado, primeiramente devido à
Os revestimentos obtidos por galvanização e
dissolução do manganês protegendo o zinco
os eletrodepósitos de ligas de zinco, entre
e, em segundo lugar, por assegurar a
outros, são de grande interesse para as
formação de um produto de corrosão de
indústrias, principalmente automobilísticas,
baixa solubilidade. Para ligas Zn-Mn contendo
onde se necessitam de camadas superficiais
11% de Mn em meio de sulfato de sódio, tem
com alta resistência à corrosão, podendo
sido proposto que a corrosão da liga leva à
permitir a aderência de tintas para o
formação de sulfato básico de zinco como
acabamento final (PUSHPAVANAM, 1991;
produto de corrosão. Revestimentos de Zn-Mn
SWATHIRAJAM, 1986).
compactos e com diferentes morfologias
Em anos recentes, algumas ligas de zinco podem ser obtidos a baixas densidades de
depositadas sobre aço por eletrodeposição, corrente. Corrente pulsada pode também ser
especificamente Zn-Ni, Zn-Co e Zn-Mn, foram usada na eletrodeposição deste tipo de
propostas como revestimentos protetores revestimento e os resultados mostram que,
contra a corrosão. Estas ligas são para uma mesma faixa de densidade de
termodinamicamente mais estáveis e, corrente média, todos os métodos levaram à
portanto, mais duradouras que revestimentos formação de camadas com teores de Mn
de zinco de mesma espessura. maiores ou similares aos obtidos por corrente
direta, porém as eficiências de corrente foram
Os eletrodepósitos de ligas Zn-Ni têm atraído
menores (MULLER; SARRET; ANDREU, 2002,
grande interesse para aplicação em
2003).
substituição a revestimentos de cádmio por
apresentarem excelente resistência à
corrosão, ductibilidade, baixa fragilização
pelo hidrogênio e alta estabilidade química.
As ligas com 10% a 15% de níquel
107
4.5 LIGAS DE COBALTO 4.6 LIGAS DE MOLIBDÊNIO

As ligas de cobalto são conhecidas na As ligas de molibdênio e níquel vêm sendo
literatura por sua aplicação em diversas estudadas nos últimos anos como uma
tecnologias (CHU; LIANG; HAO, 2014; WEI et alternativa para os banhos de cromo. O uso
al., 2017; YAPONTSEVA; DIKUSAR; do molibdênio se deve às suas propriedades
KYBLANOVSKII, 2014), as quais estão e com isso suas aplicações. O molibdênio
diretamente associadas aos outros possíveis pode ser usado em aplicações metalúrgicas,
elementos que se ligam ao cobalto formando tais como aço inoxidável e ligas de ferro
a liga. fundido. A adição de molibdênio, no caso das
ligas de aço, melhora as características
Ligas de Co-W, por exemplo, recaem em
comerciais da liga, aumentando a força e a
duas conhecidas aplicações largamente
resistência térmica, bem como a redução do
discutidas: elas podem ser desenvolvidas
potencial de corrosão do material
para fins de revestimento ou até mesmo como
(ALLAHYARZADEH; ROOZBEHANI; ASHRAFI,
alternativa na produção de energias limpas,
2011a).
em células combustíveis (BERA;
SEENIVASAN; RAJAM, 2013; SHENG et al., Atualmente, há um grande interesse por
2018). Além de possuírem elevada atividade revestimentos de molibdênio, devido às suas
eletrocatalítica (SHENG et al., 2018), devido boas propriedades mecânicas, elétricas e
às propriedades excelentes do tungstênio, as ótima resistência à corrosão. Como este metal
ligas de Co-W possuem elevadas dureza, não pode ser eletrodepositado em seu estado
resistências elétrica, térmica, alta resistência puro é realizada uma codeposição do tipo
contra choques mecânicos, assim como induzida. Este processo foi apresentado por
elevada resistência à corrosão (BERA; Brenner, onde o mesmo explica que quando
SEENIVASAN; RAJAM, 2013). um metal não pode ser eletrodepositado a
partir de uma solução aquosa na sua forma
Outras ligas de cobalto, como Co-Fe,
pura, ele pode ser codepositado com algum
possuem diferentes aplicações,
dos metais do 8º grupo como o Ni e o Fe,
principalmente em sistemas microeletrônicos
formando assim uma liga (BELTOWSKA-
e microeletromecânicos na produção de
LEHMAN; INDYKA, 2012).
sensores e atuadores (ZHANG; IVEY, 2018). O
que grande parte das ligas de cobalto Devido à alta resistência à corrosão, ao
possuem, no processo de obtenção por desgaste e ao baixo sobrepotencial de
eletrodeposição, é fato de óxidos serem evolução do hidrogênio da liga Ni-Mo, foram
formados na superfície do cátodo, o que faz dedicadas muitas investigações para
com que o cobalto, geralmente menos nobre, deposição de tal liga. Estudos mostraram que
seja preferencialmente reduzido. Como esta liga apresenta brilho, homogeneidade,
consequência, o pH em torno do cátodo aderência e caráter metálico, além de uma
torna-se mais elevado, fazendo com que a ótima eficiência catódica, chegando a
taxa de produção de hidrogênio seja também aproximadamente 92% (ABDEL HAMID;
elevada, o que comprova, em muitos casos, a HASSAN, 2012; ALLAHYARZADEH;
baixa eficiência de deposição (KIELING, ROOZBEHANI; ASHRAFI, 2011b;
1997). BELTOWSKA-LEHMAN et al., 2012).
REFERÊNCIAS
[1] ABDEL HAMID, Z.; HASSAN, H. B. [3] ALLAHYARZADEH, M. H.; ROOZBEHANI,
Influence of electrodeposition parameters on the B.; ASHRAFI, A. Electrodeposition of high Mo
characteristics of NiMoP film. Surface and Coatings content amorphous/nanocrystalline Ni-Mo alloys
Technology v. 212, p. 37–45 , 2012. using 1-ethyl-3-methyl-imidazolium chloride ionic
liquid as an additive. Electrochimica Acta v. 56, n.
[2] ALLAHYARZADEH, M. H.; ROOZBEHANI,
27, p. 10210–10216 , 2011b.978-1-56677-986-9.
B.; ASHRAFI, A. Electrodeposition of high Mo
content amorphous/nanocrystalline Ni-Mo alloys [4] ALMEIDA, M. R. H. DE; BARBANO, E. P.;
using 1-ethyl-3-methyl-imidazolium chloride ionic CARVALHO, M. F. DE; CARLOS, I. A.; SIQUEIRA, J.
liquid as an additive. Electrochimica Acta v. 56, n. L. P.; BARBOSA, L. L. Electrodeposition of copper-
27, p. 10210–10216 , 2011a. Disponível em: zinc from an alkaline bath based on EDTA. Surface
<http://dx.doi.org/10.1016/j.electacta.2011.09.011> and Coatings Technology, v. 206, n. 1, p. 95–102,
.978-1-56677-986-9. 2011.
108
[5] AZEVEDO, C. A. C.; WORMSBECKER, L. [17] CHEN, H. et al. Wear and Corrosion
E.; ADAIME, L. M.; MILEK, J. L.; SCHEER, S. Properties of Crystalline Ni-W Alloy Coatings
Corrosão em estruturas metálicas. Disponível em: Prepared by Electrodeposition. Materials Science
<http://www.cesec.ufpr.br/metalica/patologias/corr Forum, v. 849, p. 671–676, 2016.
osao/corrosao.htm> Acesso em: 19 jan. 2015.
[18] CHU, Qingwei; LIANG, Jun; HAO,
[6] BAJAT, J. B.; STANKOVIC, S.; JOKIC, B. Jingcheng. Electrodeposition of zinc-cobalt alloys
M.; STEVANOVIC, S. I. Corrosion stability of Zn-Co from choline chloride-urea ionic liquid.
alloys deposited from baths with high and low Co Electrochimica Acta v. 115, p. 499–503 , 2014.
content – The influence of deposition current Disponível em:
density. Surface and Coatings Technology. 204 <http://dx.doi.org/10.1016/j.electacta.2013.10.204>
(16-17) (2010) 2745-2753 .0013-4686.
[7] BASKARAN, I.; SANKARA NARAYANAN, [19] ĆIROVIĆ, N. et al. Synthesis, structure and
T. S. N.; STEPHEN, A. Pulsed electrodeposition of properties of nickel-iron-tungsten alloy
nanocrystalline Cu-Ni alloy films and evaluation of electrodeposits: Part II: Effect of microstructure on
their characteristic properties. Materials Letters, v. hardness, electrical and magnetic properties.
60, p. 1990–1995, 2006. Science of Sintering, v. 48, n. 1, p. 1–16, 2016.
[8] BELTOWSKA-LEHMAN, E. et al. [20] CONDE, A.; ARENAS, M. A.; DE
Electrodeposition and characterisation of DAMBORENEA, J. Electrodeposition of Zn–Ni
nanocrystalline Ni-Mo coatings. Surface and coatings as Cd replacement for corrosion
Coatings Technology v. 211, p. 67–71 , 2012. protection of high strength steel, Corrosion
Disponível em: Science. 53 (4) (2011) 1489-1497
<http://dx.doi.org/10.1016/j.surfcoat.2011.10.011>.
[21] COSTA, J. D. et al. Effect of
[9] BELTOWSKA-LEHMAN, E.; INDYKA, P. electrochemical bath composition on the
Kinetics of Ni-Mo electrodeposition from Ni-rich preparation of Ni-W-Fe-P amorphous alloy.
citrate baths. Thin Solid Films v. 520, n. 6, p. 2046– International Journal of Electrochemical Science, v.
2051 , 2012. Disponível em: 13, n. 3, p. 2969–2985, 2018.
<http://dx.doi.org/10.1016/j.tsf.2011.10.024>.
[22] EL-CHIEKH, F.; EL-HATY, M. T.;
[10] BERA, Parthasarathi; SEENIVASAN, H; MINOURA, H.; MONTASER, A. A. Electrodeposition
RAJAM, K S. Studies on surface structure, and characterization of Cu-Ni-Zn and Cu-Ni-Cd
morphology and composition of Co–W coatings alloys. Electrochimica Acta, v. 50, p. 2857– 2864,
electrodeposited with direct and pulse current 2005.
using gluconate bath. Surface Review and Letters
[23] GAMBURG, Y. D.; ZANGARI, G. Theory
v. 20, n. 01, p. 1350006 , 2013.
and Practice of Metal Electrodeposition, 2011.
[11] BOSHKOV, N; PETROV, K; RAICHEVSKI, (Boston, New York, Heidelberg: Springer)
G. Corrosion behavior and protective ability of
[24] GAVRILA, M.; MILLET, J. P.; MAZILLE, H.;
multilayer galvanic coatings of Zn and Zn-Mn alloys
MARCHANDISE, D.; CUNTZ, J. M. Corrosion
in sulfate containing médium. Surface and Coatings
behaviour of zinc-nickel coatings electrodeposited
Technology. 200 (20-21) (2006) 5995-6001.
on steel. Surface and Coatings Technology. 123 (2-
[12] BRENNER, A. Electrodeposition of Alloys: 3) (2000) 164-172.
Principles and Practice. First Edit ed. New York:
[25] GNANAMUTHU, R. M.; MOHAN, S.;
Academic Press Inc., 1963. v. I
SARAVANAN, G.; LEE, C. W. Comparative study on
[13] BROOMAN, E. W. Wear behavior of structure, corrosion and hardness of Zn– Ni alloy
environmentally acceptable alternatives to deposition on AISI 347 steel aircraft material.
chromium coatings: Cobalt-based and other Journal of Alloys and Compounds. 513 (2012) 449-
coatings. Metal Finishing, v. 102, n. 10, p. 42–54, 454.
2004.
[26] GORANOVA, D.; AVDEEV, G.; RASHKOV,
[14] CAPEL, H.; SHIPWAY, P. H.; HARRIS, S. J. R. Electrodeposition and characterization of Ni-Cu
Sliding wear behaviour of electrodeposited cobalt- alloys. Surface and Coatings Technology, v. 240, p.
tungsten and cobalt-tungsten-iron alloys. Wear, v. 204–210, 2014.
255, n. 7–12, p. 917–923, 2003.
[27] GENTIL, Vicente. Corrosão. 5. ed. Rio de
[15] CARLOS, I. A.; ALMEIDA, M. R. H. DE. janeiro: LTC, 2007.
Study of the influence of the polyalcohol sorbitol on
[28] HALIM J; ABDEL-KARIM R; EL-RAGHY S;
the electrodeposition of copper-zinc films from a
NABIL M; WAHEED A. Electrodeposition and
non-cyanide bath. Journal of Electroanalytical
characterization of nanocrystalline Ni-Mo catalysts
Chemistry, v. 562, p. 153–159, 2004.
for hydrogen production. Journal of Nanomaterials,
[16] CESIULIS, H.; XIE, X.; PODLAHA- 2012.
MURPHY, E. Electrodeposition of Co-W alloys with
[29] HAMMAMI, O.; DHOUIBI, L.; TRIKI, E.
P and Ni. Medziagotyra, v. 15, n. 2, p. 115–122,
Influence of Zn–Ni alloy electrodeposition
2009.
109
techniques on the coating corrosion behaviour in morphological and magnetic characterization of

chloride solution. Surface and Coatings electrodeposited Co-Fe-W alloys. Journal of Alloys
Technology. 202 (19) (2009) 2863-2870. and Compounds, v. 611, p. 243–248, 2014.
[30] HAMID ZA; HASSAN HB. Influence of [43] OHGAI, T.; TANAKA, Y.; FUJIMARU, T.
electrodeposition parameters on the characteristics Soft magnetic properties of Ni-Cr and Co-Cr alloy
of NiMoP film. Surface and Coatings Technology, thin films electrodeposited from aqueous solutions
212:37-45, 2012. containing trivalent chromium ions and glycine.
Journal of Applied Electrochemistry, v. 42, n. 10, p.
[31] KARAHAN, I. H.; ÖZDEMIR, R. Effect of Cu
893–899, 2012.
concentration on the formation of Cu1-x Znx shape
memory alloy thin films. Applied Surface Science, v. [44] OLIVEIRA, A. L. M. et al. Studies on
318, p. 100–104, 2014. electrodeposition and characterization of the Ni–W–
Fe alloys coatings. Journal of Alloys and
[32] KIELING, Virginia Costa. Parameters
Compounds, v. 619, p. 697–703, 2015.
influencing the electrodeposition of Ni-Fe alloys.
Surface and Coatings Technology v. 96, n. 2–3, p. [45] PUSHPAVANAM, M.; NATARAJAN, K.;
135–139 , 1997. BALAKRISHNAN, K.; SHARMA, L.R. Corrosion
behavior of electrodeposited zinc-nickel alloys. J.
[33] KUZNETSOV, V. V. et al. Corrosion
Appl. Electrochem., v. 21, p. 642 – 645, 1991.
resistance of Cr–C–W alloys produced by
electrodeposition. Journal of Solid State [46] SANTANA, R. A. C.; PRASAD, S.; DE
Electrochemistry, v. 19, n. 9, p. 2545–2553, 2015. SANTANA, F. S. M. Revestimento eletrolítico com
uma liga amorfa de Ni-W-B, resistente à corrosão e
[34] KUZNETSOV, V. V. et al. Corrosion
ao desgaste. Ecletica Quimica, v. 28, n. 1, p. 69–
resistance of Co−Cr−W coatings obtained by
76, 2003.
electrodeposition. Journal of Solid State
Electrochemistry, v. 22, n. 7, p. 2267–2276, 2018. [47] SENNA, L. F. DE; DÍAZ, S. L.; SATHLER, L.
Hardness analysis and morphological
[35] LANDOLT, D. Electrochemical and
characterization of copper-zinc alloys produced in
materials science aspects of alloy deposition.
pyrophosphate-based electrolytes. Materials
Electrochimica Acta. 39: 1075, 1994.
Research, v. 8, n. 3, p. 275–279, 2005.
[36] LASSNER, E.; SCHUBERT, W.-D. Tungsten
[48] SHENG, Minqi et al. Co-W/CeO2
Properties, Chemistry, Technology of the Element,
composite coatings for highly active
Alloys, and Chemical Compounds. first edit ed.
electrocatalysis of hydrogen evolution reaction.
New York: Kluwer Academic/Plenum Publishers,
Journal of Alloys and Compounds v. 743, p. 682–
1999.
690 , 2018.
[37] LI, B. et al. Electrodeposition and
[49] SPIEGELBERG, W. D., "Elastic resilience
characterization of Fe-Cr-P amorphous alloys from
and related properties in electronic connector
trivalent chromium sulfate electrolyte. Journal of
alloys selection", Conference: ASM
Alloys and Compounds, v. 453, n. 1–2, p. 93–101,
INTERNATIONAL'S 3RD ELECTRONIC MATERIALS
2008.
& PROCESSING CONGRESS, San Francisco,
[38] MOURA, E. S. Desenvolvimento de ligas California, 20-23 Aug. 1990, p 215-223.
de Ni-Fe-Mo, resistentes à corrosão para o
[50] SWATHIRAJAN, S. Potentiodynamic and
revestimento de equipamentos ligados ao
galvanostatic stripping methods for
transporte e obtenção do petróleo. Monografia
characterization of alloy electrodeposition process
(Graduação). Universidade Federal de Campina
and product. Journal of the Electrochemical
Grande, Campina Grande, 2005.
Society. 133, 4, 671-80, 1986.
[39] MULLER, C; SARRET, M; ANDREU, T.
[51] TAVOOSI, M.; BARAHIMI, A. Corrosion
Electrodeposition of Zn-Mn alloys at low current
behavior of amorphous–nanocrystalline Fe–Ni–Cr
densities. Journal Electrochemistry Society. 149
electrodeposited coatings. Surfaces and Interfaces,
(2002) C600-C609.
v. 8, p. 103–111, 2017.
[40] MULLER, C; SARRET, M; ANDREU, T.
[52] TSYNTSARU, N. et al. Iron-tungsten alloys
ZnMn alloys obtained using pulse reverse and
electrodeposited under direct current from citrate-
superimposed current modulations. Electrochimica
ammonia plating baths. Surface and Coatings
Acta. 48 (17) (2003) 2397-2404.
Technology, v. 203, n. 20–21, p. 3136–3141, 2009.
[41] MUN, S.-J. et al. Mechanical and
[53] TSYNTSARU, N. et al. Modern trends in
Structural Characteristics of Electrodeposited Ni-
tungsten alloys electrodeposition with iron group
Fe-W Alloy after Heat-Treatment. Journal of The
metals. Surface Engineering and Applied
Electrochemical Society, v. 157, n. 3, p. D177–
Electrochemistry, v. 48, n. 6, p. 491–520, 2012.
D180, 2010.
[54] TOMACHUK Santos, C. R. Codepósitos de
[42] NOCE, R. DELLA; BENEDETTI, A. V.;
ZnCo - Estudo do processo de sua eletrodeposição
MAGNANI, M.; PASSAMANI, E. C.; KUMAR, H.;
CORNEJO, D. R.; OSPINA, C. A. Structural,
110
por corrente pulsada. Tese (Doutorado). [58] WHITE, R. E.; BOCKRIS, J. O’M.;
Universidade Estadual de Campinas, 1998. CONWAY, B. E. Modern Aspects of
Electrochemistry, New York, Plenum Press, 1995.
[55] VREESE, P. DE; SKOCZYLAS, A.;
MATTHIJS, E.; FRANSAER, J.; BINNEMANS, K. [59] YAPONTSEVA, Yu. S.; DIKUSAR, A. I.;
Electrodeposition of copper-zinc alloys from an KYBLANOVSKII, V. S. Study of the composition,
ionic liquid-likecholine acetate electrolyte. corrosion, and catalytic properties of Co-W alloys
Electrochimica Acta, v. 108, p. 788–794, 2013. electrodeposited from a citrate pyrophosphate
electrolyte. Surface Engineering and Applied
[56] WANG, S.; GUO, X.; YANG, H.; DAI, J.;
Electrochemistry v. 50, n. 4, p. 330–336 , 2014.
ZHU, R.; GONG, J.; PENG, L.; DING, W.
Disponível em:
Electrodeposition mechanism and characterization
<http://link.springer.com/10.3103/S1068375514040
of Ni-Cu alloy coatings from a eutecticbased ionic
139>.1068375514040.
liquid. Applied Surface Science, v. 288, p. 530–
536, 2014. [60] ZHANG, Yahui; IVEY, Douglas G.
Electrodeposition of nanocrystalline CoFe soft
[57] WEI, Yongsheng et al. Synthesis of
magnetic thin films from citrate-stabilized baths.
bifunctional non-noble monolithic catalyst Co-W-
Materials Chemistry and Physics v. 204, p. 171–178
P/carbon cloth for sodium borohydride hydrolysis
, 2018. Disponível em:
and reduction of 4-nitrophenol. International Journal
<https://doi.org/10.1016/j.matchemphys.2017.10.04
of Hydrogen Energy v. 42, n. 41, p. 25860–25868 ,
3>.1362-4962 (Electronic)r0305-1048 (Linking).
2017. Disponível em:
<https://doi.org/10.1016/j.ijhydene.2017.08.148>.03 [61] ZOLDAN, V. C. Eletrodeposição de Fe em
60-3199. Si monocristalino. Dissertação (Mestrado).
Universidade Federal de Santa Catarina,
Florianópolis, 2007.
111
Capítulo 11
Ramon Sales de Araújo Batista

Carlos José de Araújo
Resumo: Dado o panorama atual das indústrias em todo o mundo e o imperativo

pela sobrevivência econômica, a necessidade pela otimização dos processos de
fabricação e obtenção de mais agilidade na criação de novos produtos tem
crescido bastante. A Prototipagem Rápida (PR) surgiu com a intenção de auxiliar
nesse quadro de necessidade pelo avanço na área de projeto e fabricação, através
da produção de protótipos, moldes, ou até peças. Dentre as tecnologias de PR
tem-se a técnica de Modelagem por Fusão e Deposição (FDM) que se baseia na
construção de um modelo tridimensional por deposição de camada por camada. O
método de FDM se utiliza de uma máquina especifica, popularmente conhecida
como “impressora 3D”, que possui alguns parâmetros de processamento para
realizar a manufatura. Esses parâmetros podem vir a ser considerados críticos
quando alterados de forma indevida, assim podendo causar avarias à peça final,
danos ao equipamento de trabalho e custos desnecessários com matéria-prima. A
partir dessa constatação, surge a evidente necessidade de analisar esses
parâmetros e suas variações durante o processo de fabricação. Nesse contexto,
este trabalho teve por objetivo propor um projeto de caracterização mecânica em
materiais processados por FDM, através do estudo dos parâmetros críticos de
Largura da Parede, Espessura da Camada Construtiva e a Porcentagem de
Preenchimento, para um equipamento MakerBot Replicator 2x, se utilizando das
técnicas de planejamento e análise de experimentos. Assim, por meio da pesquisa
experimental foi possível determinar as condições de processo em que se obtêm
aumento da tensão máxima de tração suportada, encontrando o valor de 34,80
Mpa (experimental) e em comparação com o mesmo observou-se o valor de 34,30
Mpa (predito). Constatando assim que os modelos estatísticos apresentados foram
adequados para a análise de resistência mecânica à tração.
Palavras-chave: FDM, Impressão 3D, Planejamento Experimental, Otimização,

Resistência Mecânica.
112
1 INTRODUÇÃO (AZEVEDO et al, 2016). Lopes (2011) destaca

que os ensaios mecânicos permitem a
O grande desenvolvimento tecnológico
determinação de propriedades mecânicas
vivenciado pelo mercado atual exigiu das
que se referem ao comportamento do material
empresas um maior empenho para inserir
quando submetido à ação de esforços, e que
novos produtos nas suas linhas de produção
são expressas em função de tensões e/ou
em menos tempo e mantendo o nível de
deformações.
flexibilidade e inovação elevados. Segundo
Mello et al (2010), as industrias podem valer- Nesse cenário, o objetivo deste trabalho foi,
se de artifícios para a otimização do tempo de através da variação de três parâmetros do
desenvolvimento de produtos, e uma dessas processo, avaliar a resposta mecânica de
táticas é o uso da prototipagem rápida (PR). corpos de prova produzidas por FDM
De acordo com Oliveira et al (2009), é utilizando polímero ABS e uma impressora 3D
possível conceituar a prototipagem rápida MakerBot Replicator 2x.
como um processo que simboliza a união de
tecnologias de construção física de peças a
partir de um modelo original digital em CAD. 2 MATERIAIS E MÉTODOS
Atualmente, para Dizon et al (2018), a PR Com o intuito de se possibilitar a elaboração
também muito nomeada como Manufatura de um procedimento de ensaios mecânicos
Aditiva (MA) está sendo seriamente do tipo tração, em corpos produzidos através
considerada como produtora de materiais de manufatura aditiva, por meio do processo
para várias aplicações, a saber, construção, de FDM, e utilizando como matéria prima o
vestuário, odontologia, medicina, eletrônica, ABS (Acrilonitrilo Butadieno Estireno), a
automotivo, robôs, militares, oceanografia, princípio o estudo fundamentou-se em uma
aeroespacial e outros. consistente revisão da literatura cientifica
pertinente ao assunto, também se realizando
A aplicação de manufatura aditiva ou
o estudo das normas regulamentadoras.
impressão 3D na medicina, por exemplo,
pode fornecer muitos benefícios, incluindo: a Segundo Schaeffer (2017) um experimento é
personalização de produtos médicos, um procedimento no qual alterações
equipamentos, custo-efetividade, aumento de propositais são feitas nas variáveis de entrada
produtividade, a democratização do design e de um processo ou sistema, de modo que se
fabricação, e a colaboração aprimorada. No possa avaliar as possíveis alterações sofridas
entanto, deve ser advertido que, apesar dos pela variável resposta. Portanto, para modelar
recentes avanços médicos significativos e o presente trabalho foi feito uso das técnicas
empolgantes envolvendo impressão 3D, de planejamento e análise de experimentos,
desafios científicos e regulamentares notáveis ou mais precisamente um DOE – Design of
permanecem (VENTOLA, 2014). Experiments.
Dentro do contexto de manufatura aditiva, um
dos processos de prototipagem rápida
2.1 PLANEJAMENTO FATORIAL
bastante conhecido é o de modelagem por
fusão e deposição (FDM). De acordo com Este estudo buscou analisar a relação entre a
Bakar et al (2010), cada modelo FDM variação de determinados parâmetros críticos
disponível comercialmente tem diferentes do processo de FDM e a resposta de algumas
tipos de parâmetros de processo para propriedades físicas e mecânicas do protótipo
aplicações diferentes. Em conformidade com construído. Dentre as diversas variáveis que
Anitha et al (2001), uma vez que este tem relação com a manufatura de corpos e
processo de FDM pode ser usado para uma peças através de FDM, foram selecionadas as
variedade de aplicações e o custo do seguintes: largura da parede (núm. de
protótipo pode ser geralmente elevado, existe cascas), espessura da camada construtiva
uma necessidade de aperfeiçoar os (precisão do equipamento) e porcentagem de
parâmetros do processo, tanto do ponto de preenchimento. E em um segundo passo foi
vista econômico quanto tecnológico. delimitado os níveis de variação desses
fatores, baseando-se na capacidade da
Ao desenvolver um projeto de Engenharia, é
impressora 3D, como pode ser visto na
de grande importância que o engenheiro
Tabela 1
tenha conhecimento de propriedades
mecânicas fundamentais dos materiais que .
serão empregados na execução do projeto
113
Tabela 1. Seleção dos parâmetros de processo a serem estudados e seus respectivos níveis de
variação.
Níveis dos fatores
Fatores
Mín Central Máx
Largura da Parede (núm. de cascas) 2 6 10
Espessura da Camada Construtiva (mm) (Precisão) 0.3 0.2 0.1
Porcentagem de Preenchimento (%) 10 55 100
Ainda trabalhando com base nas técnicas de ao planejamento fatorial, utilizando o software
planejamento e análise de experimentos, foi STATISTICA 7.0. Na Tabela (2) apresenta-se a
gerada uma matriz de delineamento referente matriz do planejamento adotado.
Tabela 2. Matriz de delineamento do planejamento fatorial.

Ensaios Aleatoriedade Largura da Parede Precisão Preenchimento
1 4 2 0.3 10%
2 1 10 0.3 10%
3 6 2 0.1 10%
4 9 10 0.1 10%
5 7 2 0.3 100%
6 3 10 0.3 100%
7 5 2 0.1 100%
8 12 10 0.1 100%
9 2 6 0.2 55%
10 11 6 0.2 55%
11 10 6 0.2 55%
12 8 6 0.2 55%
A partir da geração da matriz de 2.2.1 CORPOS DE PROVA

delineamento do planejamento fatorial, foi
Ahn et al. (2002) afirmam que a norma padrão
possível identificar que seriam necessários
para ensaios de tração de materiais
doze ensaios diferentes, para que assim fosse
poliméricos é a ASTM D638 (Método Padrão
possível abarcar de forma confiável a
de Teste para Propriedades de Tração de
combinação dos parâmetros e seus
Plásticos). Com base nesta norma, foi criado e
respectivos níveis, que tinham sido
desenhado um modelo em ambiente CAD,
estabelecidos no início do estudo. Esses
utilizando o código comercial Autodesk
ensaios obedeceram a uma ordem de
Inventor Professional 2016. Na sequência,
aleatoriedade, necessária para evitar a
esse modelo foi exportado para o formato
propagação ou influência de possíveis efeitos
STL, que é o formato comumente aceito nos
de outros parâmetros não controláveis.
programas que gerenciam as máquinas de
impressão 3D.
2.2 MANUFATURA ADITIVA POR FDM
Tendo sido definida a matriz do planejamento
fatorial, foi iniciada a construção dos corpos
de prova através do processo de FDM.
114
2.2.2 IMPRESSORA 3D E MATERIAL DE copolímero ABS, original da MakerBot. A

IMPRESSÃO Figura 1 permite visualizar a máquina de
impressão e o material utilizado na
A máquina de impressão 3D utilizada para a
construção.
fabricação dos corpos de prova foi uma
MakerBot Replicator 2 Desktop 3D Printer. O
material escolhido para o projeto foi o
Figura 1. Impressora 3D e Material utilizado na Impressão. (a) Impressora MakerBot Replicator 2

Desktop 3D Printer. (b) Rolo de filamento ABS original MakerBot.
a) b)
Fonte: Adaptado de MakerBot Replicator 2 – User Manual.
Foram então produzidos os doze corpos de trabalho, já que possuindo três réplicas de
prova definidos na Tab. (2), variando os cada configuração é possível extrair uma
parâmetros do processo para cada corpo. média mais representativa das respostas às
Esse procedimento foi repetido por mais duas quais se desejava avaliar. Na Figura 2 é
vezes, para assim totalizar 36 corpos de possível visualizar uma das réplicas dos doze
prova. Essa repetição foi realizada para corpos de prova.
atribuir uma garantia estatística mínima ao
Figura 2. Doze corpos de prova identificados segundo o planejamento fatorial.
2.3 ETAPA DE OBTENÇÃO DAS ensaios de tração para obtenção das curvas
PROPRIEDADES tensão x deformação que caracterizam o
comportamento mecânico desses corpos. Os
Após realizar a impressão dos corpos de
ensaios em questão foram baseados na
prova, os mesmos foram submetidos a
norma ASTM D638 – 2ª. Para realizar os
115
ensaios de tração foi utilizada uma máquina hidráulica, modelo MTS 810, a mesma pode
universal de ensaios de materiais, servo- ser visualizada na Figura 3.
Figura 3. Máquina de ensaios de materiais servo-hidráulica modelo MTS 810.
3 RESULTADOS E DISCUSSÕES configuração. Ou seja, para cada

configuração foram realizados três ensaios e
3.1 ENSAIOS DE TRAÇÃO
em seguida calculada a curva média. Assim,
Na Figura 4 são apresentados os resultados foram geradas doze curvas tensão x
dos ensaios de tração, mostrando as curvas deformação para as doze respectivas
médias referentes às três réplicas de cada configurações.
Figura 4. Curvas tensão x deformação para as doze configurações de corpos de prova.
Com base nesses ensaios foi possível obter foram usados como entradas nas análises
os dados de Tensão Máxima e Deformação estatísticas realizadas. A Tabela 3 permite
na Tensão Máxima para cada um dos doze visualizar um resumo dos dados obtidos nos
corpos de prova. Posteriormente esses dados ensaios.
116
Tabela 3. Matriz de dados resultantes dos Ensaios mecânicos de tração.

Dados Obtidos
Propriedades Mecânicas
Corpos de Prova (Mín - Máx)
Tensão Máxima (MPa) 17.1 - 34.8
Deformação na Tensão Máx. (%) 3.1 - 5.35
Castelo Branco (2016) afirma em sua orientação da impressão e o ângulo de raster

pesquisa que ao se analisar as taxas de torna-se possível indicar que o filamento mais
deposição, percebe-se que quanto maior a grosso oferece melhores propriedades
área de preenchimento maior será a tensão mecânicas se impresso nas direções x e z e
máxima e a tensão de ruptura, filamento mais fino é ideal para impressões na
consequentemente mais rígida será a peça, e direção y. Com isso, frente a essa análise e ao
ele destaca aquelas com taxa de deposição presente estudo, percebesse que o parâmetro
de 100% (peça maciça). No caso do presente da espessura da camada construtiva
trabalho, foi possível constatar que o ensaio necessita ser estudado em paralelo com a
que apresentou melhor índice de resposta direção de impressão, para que assim se
(ensaio 7), foi também um dos ensaios que obtenha dados mais significativos de sua
possuía taxa de deposição de 100%. Porém influência nas propriedades mecânicas das
percebeu-se que nem todos os ensaios que peças. Logo, sabendo que nessse trabalho a
possuíam taxa de preenchimento de 100% direção de impressão não foi considerada
garantiram índices aproximados do máximo como parâmetro de entrada, percebe-se que
encontrado que foi 34,8 Mpa, como foi o caso, não foi possivel avaliar em termos
por exemplo, do ensaio 6 que só alcançou comparativos o efeito da variação da
28,3 Mpa. Isso reforça a importância da espessura da camada construtiva. Podemos
interação como os demais parâmetros do observar melhor comparando o ensaio 3
processo. (menor índice de resposta) com o ensaio 7
(maior índice de resposta), pois ambos
Os resultados que apresentaram os menores
tiveram o mesmo nível de espessura da
índices no geral foram os ensaios com a taxa
camada construtiva, constatando assim a
de preenchimento de apenas 10% e com a
pouca influência desse parâmetro frente aos
largura da parede de apenas 2 cascas, foram
demais considerados no estudo.
eles os ensaios 1 e 3.
Hossain (2013) relata que a diminuição da
largura de parede resulta em aumento da 3.2 ANÁLISE DOS RESULTADOS
resistência à tração. Porém foi possível avaliar ESTATÍSTICOS
dentro do corrente estudo que essa afirmativa
Tendo em mãos os dados referentes à
só se torna válida se associada a uma alta
variável de reposta, Tensão Máxima, que foi
taxa de preenchimento, como podemos
selecionada como aquela que seria analisada,
observar comparando os ensaios 1 e 3
foi possível gerar as tabelas do projeto de
(baixos níveis de largura da parede e
experimento (DOE). Na Tabela 4 são
preenchimento) com os ensaios 5 e 7 (baixos
mostrados os valores utilizados no DOE para
níveis de largura da parede porém altos níveis
os três parâmetros críticos inicialmente
de preenchimento).
apresentados, e na Tabela 5 são
Segundo Bagsik (2011) ao estudar os efeitos apresentados valores codificados e a
da espessura da camada construtiva, a resposta de Tensão Máxima (Y).
Tabela 4. Valores utilizados no DOE para os três parâmetros previamente definidos.

Variáveis Código -1 0 1
Largura da Parede x1 2 6 10
Espessura (mm) (Precisão) x2 0.3 0.2 0.1
Preenchimento (%) x3 10 55 100
117
Tabela 5. Valores codificados e resposta de Tensão Máxima (Y).

Ensaios x1 x2 x3 Y (Mpa)
1 -1 -1 -1 17.7
2 1 -1 -1 29.4
3 -1 1 -1 17.1
4 1 1 -1 33.3
5 -1 -1 1 29.1
6 1 -1 1 28.3
7 -1 1 1 34.8
8 1 1 1 29.7
9 0 0 0 27.9
10 0 0 0 30.7
11 0 0 0 29.1
12 0 0 0 28.8
A partir dos resultados obtidos foram porcentagem de preenchimento (x3) e a

calculados os coeficientes de regressão interação entre essas duas variáveis (x1-x3)
apresentados na Tabela 6. Observa-se que foram estatisticamente significativas, ao nível
para a resposta Y, referente a Tensão de significância 5% (p < 0.05).
Máxima, apenas a largura da parede (x1), a
Tabela 6. Coeficientes de regressão para a resposta Y (MPa).
Coeficiente de Erro Estimativas por intervalo (95%)

Fatores t(5) p-valor
Regressão Padrão Limite Inferior Limite Superior
MÉDIA 27.99167 0.598250 46.78927 0.000000 26.45382 29.529517
x1 2.75000 0.732703 3.75322 0.013250 0.86653 4.633474
x2 1.30000 0.732703 1.77425 0.136200 -0.58347 3.183474
x3 3.05000 0.732703 4.16267 0.008801 1.16653 4.933474
x1 - x2 0.02500 0.732703 0.03412 0.974102 -1.85847 1.908474
x1 - x3 -4.22500 0.732703 -5.76632 0.002204 -6.10847 -2.341526
x2 - x3 0.47500 0.732703 0.64828 0.545392 -1.40847 2.358474
A Equação 1 descreve a Tensão Máxima contém apenas os termos estatisticamente

prevista pelo modelo em função das variáveis significativos.
codificadas, no modelo reparametrizado que
𝑌 = 27.992 + 2.750⁡𝑥1 + 3.050⁡𝑥3 − 4.225⁡𝑥1𝑥3 (1)
A análise de variância (ANOVA) indicou uma largura da parede e a porcentagem de

porcentagem de variação explicada de preenchimento. Com base na Tabela 7 é
93,17%, revelando que o modelo possível observar a análise de variância dos
representado pela Equação 1 é satisfatório dados calculados.
para analisar a Tensão Máxima frente à
118
Tabela 7. ANOVA para resposta Y.
Fatores Soma dos Quadrados Grau de Liberdade Quadrado Médio F calc p-valor
x1 60.5000 1 60.5000 14.08669 0.013250

x2 13.5200 1 13.5200 3.14797 0.136200
x3 74.4200 1 74.4200 17.32780 0.008801
x1 - x2 0.0050 1 0.0050 0.00116 0.974102
x1 - x3 142.8050 1 142.8050 33.25042 0.002204
x2 - x3 1.8050 1 1.8050 0.42027 0.545392
Erro Puro 21.4742 5 4.2948
Total SQ 314.5292 11
A Figura 5 permite observar mais algumas regressão foi considerável para a variável x1 e
informações relevantes com relação à bem significativa para a variável x3 e a
resposta Y, baseando-se no fato de que a interação x1-x3.
Figura 5. Superfície de resposta para a Tensão Máxima (MPa) em funçãoda Largura da Parede e da
Porcentagem de Preenchimento.
3.3 VALIDAÇÃO EXPERIMENTAL tensão máxima suportada, como pode ser

visto na Tabela 8. E também foi possível
A partir da análise da superfície de resposta
observar o valor predito em comparação com
foi possível determinar as condições de
o valor experimental, como está apresentado
processo em que se obtêm aumento da
na Tabela 9.
Tabela 8. Condição ótima de processo.

Cond. Ótima x1 x2 x3
Código -1 1 1
Nível 2 0.1 100
119
Tabela 9. Resultados preditos e experimentais.
Predito Experimental
Variável de Resposta
Condição Ótima
Tensão Máxima (MPa) 34.30 34.80
4 CONCLUSÕES combinação de máximo percentual de

preenchimento e mínima largura de parede
Após os experimentos realizados, e todas as
(número de cascas) é a interação responsável
analises estatísticas foi possível constatar que
por proporcionar maiores índices de
os modelos estatísticos apresentados foram
resistência mecânica à tração nas peças
adequados para a resposta de Tensão
impressas por FDM, e que o parâmetro da
Máxima (Y), pois apresentaram termos
Espessura da Camada Construtiva não é
estatisticamente significativos. Também foi
estatisticamente significativo para a resposta
possível comprovar através da metodologia
da Tensão Máxima (Y).
de superfície de resposta, que dentre os três
parâmetros críticos do processo,
selecionados ao inicio do estudo, a
REFERÊNCIAS
[1] ANITHA, R.; ARUNACHALAM, S.; performance do material PLA através da
RADHAKRISHNAN, P., 2001. Critical parameters manufatura aditiva por meio do processo FDM.
influencing the quality of prototypes in fused 2016. 71f. Dissertação (Mestrado Profissional em
deposition modelling. Journal of Materials Ciência e Tecnologia em Saúde), Centro de
Processing Technology, Vol.118, pp.385-388. Ciências e Tecnologia, Universidade Estadual da
Paraíba - UEPB, Campina Grande - PB.
[2] ANH, S. H.; MONTERO, M.; ODELL, D.;
ROUNDY, S.; WRIGHT, P.K. Anisotropic material [7] DIZON, J. R. C.; ESPERA JR, A. H.; CHEN,
properties of fused deposition modeling ABS. A.; ADVINCULA, R. C.; Mechanical characterization
Rapid Prototyping Journal. Vol 7, No. 4, pp. 248- of 3D-printed polymers. Journal Additive
257. 2002. Manufacturing, Vol. 20, 2018, pp. 44 – 67.
[3] AZEVEDO, C. M. P. de; GODINHO, L. H.; [8] HOSSAIN, M. S., J. RAMOS, D. ESPALIN,
SARAIVA, M. A. L.; ALVARENGA, S. D.; RIBEIRO, M. PEREZ and R. WICKER, Improving tensile
M. V. A. A. Ensaio Mecânico de Tensão- mechanical properties of FDM-manufactured
Deformação sob Tração: Um Estudo Integrado com specimens via modifying build parameters, in Solid
a disciplina de cálculo. Revista Perspectivas Fabrication Symposium, Austin, TX, 2013.
Online. v. 6, n. 15, 2016. Disponível em: <
[9] LOPES, Jorge Teófilo de Barros. Estrutura
http://www.
e propriedades dos materiais. Apostila. Disponível
seer.perspectivasonline.com.br/index.php/exatas_e
em:
_engenharia/article/view/1030 >. Acesso em
http://jorgeteofilo.files.wordpress.com/2011/08/epm
25/10/2018.
-apostila-capc3adtulo09-propriedades.pdf. Acesso
[4] BAGSIK, A., and V. SCHOPPNER, em 23/10/2018.
Mechanical properties of fused deposition
[10] MELLO, C. H. P.; SILVA, C. E. S. da;
modeling parts manufactured with Ultem∗9085. In
COSTA, S. C. da. Comparação de três diferentes
Proceedings of the 69th Annual Technical
tecnologias de prototipagem rápida em relação a
Conference of the Society of Plastics Engineers
critérios de custo e tempo. In: ENCONTRO
(ANTEC ’11), pp. 1294–1298, Boston, MA, May
NACIONAL DE ENGENHARIA DE PRODUÇÃO, 26,
2011.
2006, Fortaleza. Anais do 26° Encontro Nacional de
[5] BAKAR, N. S. A.; ALKAHARI, M. R.; Engenharia de Produção, 2006.
BOEJANG, E., 2010. Analysis on fused deposition
[11] OLIVEIRA, M. F.; MAIA, I. A.; NETO, P. I.;
modelling performance. Journal of Zhejiang
PEREIRA, F. D. A. S.; SILVA, J. V. L. Explorando os
University-SCIENCE A (Applied Physics &
recursos das tecnologias de prototipagem rápida -
Engineering), Vol.11 (12), pp.972-977.
SLS e FDM - em aplicações especiais. In:
[6] CASTELO BRANCO, R. R. Elaboração de CONGRESSO BRASILEIRO DE POLÍMEROS, 10,
protocolo de ensaios mecânicos para avaliação da
120
2009, Foz do Iguaçu. Anais do 10° Congresso processo-de-forjamento-a-quente/# > Acesso em:
Brasileiro de Polímeros, 2009. 03/12/2018.
[12] SCHAEFFER, L.; Projeto de Experimentos [13] USA. ASTM D638-02a - Standard Test
(DOE) aplicado a processo de Forjamento a Method for Tensile Properties of Plastics.
Quente. Revista Ferramental, Jun/2017. Disponível
[14] VENTOLA, C. L. Medical Applications for
em:
3D Printing: Current and Projected Uses. Journal
<https://www.revistaferramental.com.br/?cod=artig
P&T Community, Vol. 39, No. 10, Out. 2014, pp.
o/projeto-de-experimentos-doe-aplicado-a-
704-711
121
Capítulo 12
Tuana Maria de Medeiros Costa

Endyara de Morais Cabral
Severino Rodrigues de Farias Neto
Resumo: Em um campo petrolífero é comum se observar uma produção simultânea

de óleo, gás, água e contaminantes, dependendo das características da reserva
petrolífera. O principal interesse econômico na cadeia produtiva de petróleo, é a
produção de óleo e gás, isso torna de suma importância a existência de um
processo para separar os componentes contidos no fluido extraído. A separação é
feita por ação da gravidade em equipamentos denominados vasos separadores.
Tal processo faz parte da etapa chamada processamento primário de petróleo.
Aprofundar o estudo dos vasos separadores é extremamente necessário, pois o
processamento primário de petróleo é fundamental na indústria petroquímica.
Conhecer a dinâmica do equipamento torna possível fazer melhorias no processo,
buscando as melhores condições de funcionamento para que haja um maior
aproveitamento no processo de separação. No presente trabalho, foi utilizado o
software ANSYS/CFX® para analisar os campos de fração volumétrica, pressão e
velocidade superficial do óleo sob diferentes condições de entrada. Os resultados
mostram o que acontece com a dinâmica do líquido ao longo do tempo
apresentando os campos de pressão fração volumétrica, pressão e velocidade
superficial do óleo. O conjunto de resultados obtidos representa a dinâmica do óleo
em um vaso separador bifásico.
Palavras-chave: simulação numérica, separação, vaso separador
122
1 INTRODUÇÃO  Separação das fases oleosa, gasosa e

aquosa, nos equipamentos conhecidos
Ao longo da vida produtiva de um campo
como
petrolífero é comum se observar uma
produção simultânea de óleo, gás, água e  separadores;
contaminantes, dependendo das
características da reserva petrolífera. A  Tratar a fase oleosa para redução da
presença de água, em geral, é decorrente da água emulsionada e dos sais dissolvidos;
própria água conata ou da água injetada na  Tratar a fase gasosa para redução do teor
recuperação secundária da reserva de de água, e se necessário de outros
petróleo (Brasil et al., 2011; Zimmermann,
2014). O gás, por sua vez, costuma ser  contaminantes;
produzido misturado ao óleo (gás natural
associado) ou isolado (gás natural não
 Tratar a fase aquosa para descarte e/ou
reinjeção nos poços produtores.
associado) dependendo de sua razão de
solubilidade.
O principal interesse econômico nesta cadeia A primeira etapa do processamento primário
produtiva é produção de óleo e gás, sendo de petróleo e gás consiste na separação das
necessário dotar os campos marítimos e fases água, óleo e gás. Essa separação é
terrestres de facilidades de produção, com realizada com o auxílio de equipamentos
instalações destinadas a efetuar, sob denominados vasos separadores que,
condições controladas a separação do óleo, dependendo da quantidade relativa das fases
do gás e da água. Dessa forma, o presentes no fluido, podem ser bifásicos ou
processamento primário do petróleo, trifásicos, atuando em série ou paralelo,
conforme Figura 1, abrange (BRASIL et al., orientados vertical ou horizontalmente
2011): (SANT'ANNA, 2005). Nos separadores
bifásicos ocorre a separação gás/líquido,
enquanto que nos separadores trifásicos
ocorre, também, a separação gás/óleo/água.
Figura 1: Esquemas simplificado do Processamento Primário de Fluidos (Kunert, 2007).
Os separadores horizontais são normalmente Por outro lado, possuem desvantagens em

mais eficientes sob o ponto de vista da relação aos verticais quanto à remoção dos
separação gás/líquido, pois oferecem uma sólidos produzidos, os verticais apresentam
área superficial de interface maior que uma geometria que permite a deposição
permite uma maior decantação das gotículas localizada no fundo do vaso facilitando sua
de óleo presentes na fase gasosa, além de remoção. Outra desvantagem dos
favorecer o desprendimento do gás da fase separadores horizontais é que eles
líquida separada (SILVA, 2013). apresentam uma menor capacidade de
123
absorver grandes variações de fluxo (Kunert, 2 METODOLOGIA

2007).
2.1 DOMÍNIO FÍSICO DE ESTUDO
O presente trabalho tem como propósito
O domínio de estudo consiste de um vaso
avaliar a dinâmica do escoamento de uma
separador bifásico equipado com uma seção
mistura bifásica (gás/óleo) em um vaso
de entrada e duas seções de saída como
separador horizontal durante o processo de
mostra o esquema da Figura 2.
separação, utilizando um modelo matemático
para descrever o processo.
Figura 2: Esquema de um separador bifásico horizontal (Requena e Rodríguez, 2006).
A geometria adotada para o estudo numérico vaso separador. As seções de saída estão
é formada por um domínio simplificado e localizadas no topo e na base do separador,
bidimensional, como ilustrado a Figura 3. Os uma para gás e outra para óleo,
fluidos entram por uma tubulação situada na respectivamente.
parte superior na extremidade esquerda do
Figura 3: Geometria criada para representar o domínio físico de estudo
124
2.2 MALHA NUMÉRICA e foi cedida pelo professor Dr. Severino

Rodrigues de Farias Neto.
A malha numérica que está representada na
Figura 4, possui 87633 elementos tetraédricos
Figura 4 – Representação da malha: (a) visão geral, (b) região de entrada e (c) uma das regiões de
saída.
a)
b) c)
2.3 MODELO MATEMÁTICO  O regime de escoamento é turbulento;

A modelagem matemática definida para  O efeito da gravidade é levado em
descrever a dinâmica dos fluidos no interior consideração;
do vaso separador é baseada nas equações
de conservação de massa, momento linear e  As superfícies foram consideradas
as do modelo de turbulência k-ε (maiores sem rugosidade.
informações podem ser obtidas no trabalho
de Paiva, 2012), em conjunto com as
condições iniciais e de contorno, adotando as O modelo de superfície livre foi utilizado para
seguintes considerações: prever o escoamento e assume que as fases
presentes são contínuas e a força devido à
 O escoamento é isotérmico e os tensão superficial é modelada como uma
fluidos incompressíveis; força concentrada na interface ao invés de
uma força de superfície. O modelo de tensão
 Regime transiente;
superficial utilizado pelo CFX baseia-se no
 Não há transferência de quantidade modelo da Força de Superfície Contínua de
de movimento devido à transferência de Brackbill et al. (1992).
massa na interface gás/óleo;
 Não há reações químicas;
125
Assim, com estas considerações as equações reduzem-se a:

de conservação de massa e momento linear
Equação de conservação de massa:


t

   f   U   0 (1)
onde fa, ra e Ua são, respectivamente, a Equação da quantidade de movimento:

fração volumétrica, viscosidade e vetor
velocidade da fase a.
 U 
t
   
 
  
   f  U   U    f p    f  ef U   U     S
T
M  F (2)
 de momento devido às forças de corpo

onde p é a pressão, ef a viscosidade efetiva externas, tais como a força de empuxo e força
(viscosidade dinâmica mais a viscosidade
F
S rotacional.  representa a força de tensão
turbulenta), o termo M representa as forças superficial dada pela Equação 3:
F  r   (3)
onde rαβ e δαβ são dados pelas Equações 4 e 5:
r      n   S  (4)
   f (5)
onde σ é o coeficiente de tensão superficial, volumétrica), Ñs é o operador gradiente na

nαβ é o vetor normal à interface apontando do interface e κ é a curvatura da superfície
fluido primário ao secundário (calculado a definida pela Equação 6:
partir do gradiente de uma fração
     n
(6)
A área de contato superficial por unidade de

A
volume,  , para o modelo de superfície
livre é definida pela Equação 7:
A  f (7)
onde fα representa a fração volumétrica da Adotou-se o modelo de turbulência .

fase α
 Condições de contorno:
 Modelo de turbulência:
126
Para finalizar o modelo matemático foram Para simular um regime de escoamento

definidas as seguintes condições de contorno estratificado da mistura gás/óleo na região de
sobre as fronteiras do domínio de estudo entrada do vaso, a seção de entrada foi
(Figura 2) para o escoamento dentro do vaso subdividida em duas partes, uma para
separador. entrada de gás e outra para a entrada de
óleo, como ilustra a Figura 5.
Figura 5: Detalhe da região de entrada do vaso.
Inicialmente foi considerado que o vaso Foi adotado um valor uniforme e não nulo para
separador era preenchido com metade de as componentes de velocidade da mistura.
óleo (fgás = 0) e metade de gás (fóleo = 0) e
Na região 0 < r < R (entrada gás), as
ambos estavam estagnados, ou seja, u = v =
seguintes condições foram adotadas:
0 m.s-1.
u  U gás  x, y 
e v = 0 m.s-1
f gás  1,0 f óleo  0

e
Já na região 0 < r < -R (entrada óleo), foram

adotadas as seguintes condições:
u  U óleo e v = 0 m.s-1 x, y 
f gás  0 f óleo  1,0

e
Nas paredes da tubulação, da placa defletora As seções de saída de gás e óleo foram
e do vaso, foi adotada a condição de não- definidas como Opening, com pressão
deslizamento (No Slip Wall) para o óleo e para relativa igual a 0 Pa e Opening Pressure,
o gás (u = v = 0 m.s-1). Essa condição respectivamente.
permite que os fluidos assumam a mesma
velocidade da parede do tubo.
127
Quadro 1 - Considerações adotadas para a resolução das equações governantes.

Características Considerações
Escoamento Bifásico (gás-óleo), bidimensional, incompressível e isotérmico
Modelo de transferência de interface Modelo de superfície livre
Força gravitacional Considerada
Critério de Convergência para a massa
10-4 kg.s-1
(RMS)
Força interfacial Tensão superficial
Propriedades físico-químicas dos fluidos: mistura gás/óleo estão apresentadas na

Tabela 1. O óleo usado como referência foi
As propriedades dos fluidos utilizados na
obtido do trabalho desenvolvido por Buriti
simulação do escoamento estratificado da
(2013).
Tabela 1 - Propriedades físico-químicas dos fluidos usados nas simulações deste estudo.
Propriedades físicas Gás Óleo
Densidade (kg.m-³) 1,185 925,5
Viscosidade dinâmica (Pa.s) 1.831E-05 0,1
Tensão superficial (N.m-1) 0,072
As simulações foram realizadas em um As condições de entrada para a simulação

computador Asus, processador Intel Core i5 foram:
1.7 GHz, 8GB RAM e HD de 1TB. Para a
obtenção da simulação numérica foi utilizado
o software ANSYS CFX® Release 15.0.
Tabela 2 - Dados e condições utilizado na simulação.

Velocidade Óleo Velocidade Gás Tempo Total
-1 -1
(m.s ) (m.s ) (s)
2,5 0,5 8
3 RESULTADOS E DISCUSSÃO interior do vaso separador em diferentes

instantes de tempo (0; 0,1; 0,2; 0,4; 0,8; 1,2;
Os campos de fração volumétrica de óleo
1,6; 2,2; 4,0; 5,5; 6,5 e 8,0 s).
para a simulação estão representados na
Figura 6. Onde pode ser observado o
comportamento das fases gás e óleo no
128
Figura 6: Campos de fração volumétrica de óleo no plano XY em diferentes instantes de tempo.
Nas Figuras 7 e 8, estão representados os velocidade superficial variam com o tempo,

campos de pressão e de velocidade em consequência das modificações causadas
superficial do óleo sobre um plano pela entrada da mistura no vaso. Porém, a
longitudinal. Observa-se de uma maneira partir do instante 5,5 s observa-se um
geral que os campos de pressão e de comportamento constante.
Figura 7: Campos de pressão do óleo no plano XY em diferentes instantes de tempo.
129
Figura 8: Campos de velocidade superficial do óleo no plano XY em diferentes instantes de tempo.
4 CONCLUSÃO escoamento. Além disso, verificou-se como

uma mistura gás/óleo, ao escoar em regime
Com base nos resultados obtidos e
estratificado no tubo de entrada, comporta-se
apresentados acerca do presente estudo,
dinamicamente no interior do vaso separador
pode-se concluir que a modelagem
ao longo de um determinado intervalo de
matemática desenvolvida para o estudo do
tempo.
escoamento bifásico gás/óleo no interior do
vaso separador foi capaz de representar os
fenômenos físicos envolvidos neste
REFERÊNCIAS
[1] Brackbill, J. U.; Kothe, D.B.; Zemach, C. [6] Requena, G.; Rodríguez, M. “Diseño y
“A Continuum Method for Modelling Surface Evaluación de Separadores Bifásicos y Trifásicos”.
Tension”. Journal of Computational Physics Tese de Licenciatura. Engenharia de Petróleo.
100:335-354, 1992. Universidad Central de Venezuela, UCV. Caracas,
2006.
[2] Brasil, N.I; Araújo, M.A.S.; Sousa, E.C.M.
“Processamento Primário de Petróleo e Gás”. Rio [7] Sant'anna, A.A. “Simulação de
de Janeiro: LTC, 2011. processamento de gás natural em plataforma off-
shore.” Monografia de TCC. Escola de Química.
[3] Buriti, C. J. O. “Análise
UFRJ, 2005.
termofluidodinâmica de escoamento água-óleo em
conexões na presença de vazamento: modelagem [8] Silva, A. M. Análise do Equilíbrio Líquido-
e simulação”. Dissertação de mestrado. Programa Vapor em Vasos Separadores de Petróleo de
de Pós-Graduação em Engenharia Química. UFCG, Campos Maduros”. Monografia. Engenharia do
2013. Petróleo. Universidade Federal do Rio Grande do
Norte. Natal, 2013.
[4] Kunert, R (E.). “Processamento Primário de
Petróleo”. Petrobrás, Rio de Janeiro, 2007. [9] Zimmermann, M.S., “Modelagem e
simulação do processo de separação gás/óleo
[5] Paiva, M.K.N. “Avaliação da estabilidade
empregando um separador ciclônico”, Monografia
do escoamento anular (core-flow) numa tubulação
de TCC, Curso de Graduação em Engenharia
horizontal na presença de vazamento”. Trabalho de
Química, UFCG, 2014.
Conclusão de Curso. Curso de graduação em
Engenharia Química. UFCG, 2012.
View publication stats

Pos-Graduacao em Engenharia Mecanic

Enviado por

Direitos autorais:

Formatos disponíveis

Pos-Graduacao em Engenharia Mecanic

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Pos-Graduacao em Engenharia Mecanic

Enviado por

Direitos autorais:

Formatos disponíveis

See discussions, stats, and author profiles for this publication at: https://www.researchgate.

II Simpósio de Pós-Graduação em Engenharia Mecânica da Universidade

Book · January 2019

Raissa Queiroga Talles Oliveira

SEE PROFILE SEE PROFILE

Romulo Pierre Batista Dos Reis Antonio Almeida

SEE PROFILE SEE PROFILE

Anaerobic sludge, Anaerobic reactor, Wastewater View project

The user has requested enhancement of the downloaded file.

II Simpósio de Pós-Graduação em Engenharia

Dados Internacionais de Catalogação na Publicação (CIP)

Modo de acesso: World Wide Web

1. Engenharia 2. Engenharia Mecânica.

Comitê Organizador SiPGEM 2018

Comitê Científico SiPGEM 2018

Tem como objetivos fomentar a integração entre alunos, a interdisciplinaridade, a divulgação

Enfim agradecemos à UAEM/UFCG e a Editora Poisson que se propôs publicar os resultados

Antonio Almeida Silva

Capítulo 2: Análise numérica de protótipo de ponte estaiada com elementos

Capítulo 3: Influência dos elementos Berílio, Titânio e Boro na dureza de ligas

Capítulo 4: Análise numérica de uma prancha de skate fabricado em

Capítulo 5: Caracterização microestrutural e mecânica de uma liga memória

Capítulo 6: Análise experimental quanto aos efeitos da dilatação térmica na

Capítulo 7: Análise numérica da influência da profundidade de operação no

Capítulo 8: Detecção dos danos de impacto em materiais compósitos via

Capítulo 10: Revestimentos metálicos obtidos por eletrodeposição: Uma

Capítulo 11: Estudo de caracterização mecânica em materiais processados

Capítulo 12: Simulação numérica do processo de separação gás/ óleo

Agripino Pereira da Silva Filho

Resumo: O cenário mundial ultimamente tem se inclinado para o desenvolvimento

Palavras-Chave: Energias renováveis. Criptomoedas. Mineração. Sustentabilidade.

1 INTRODUÇÃO sua utilização, reduzirá as agressões ao meio

Figura 1. Estrutura da oferta interna de energia elétrica no Brasil.

8,20% 6,80% 10,50%

Fonte: adaptado EPE (2017).

2.2 CRIPTOMOEDAS modelo original foi o meio termo de

criptomoeda denominado de blockchain como retribuição, dependendo diretamente

09/OUT 13/OUT 17/OUT 21/OUT 25/OUT 29/OUT 02/NOV 06/NOV

Fonte: adaptado BitcoinEnergyConsumption.com

3 METODOLOGIA pelo serviço de fornecimento do último ano da

Figura 3. Equipamento de mineração Bitmain Antminer S9i com fonte original.

3.1.2 PAINÉIS SOLARES de potência nas condições padrão, com área

Figura 4. Média anual de insolação diária no Brasil.

Fonte: ATLAS Solarimétrico do Brasil. Editora

Dessa maneira, na região nordeste, onde através de estações meteorológicas, Atlante

Fonte: simulação metereologica meso-escala, (WRF-ARW) realizado

Figura 6. Comparativo da irradiação solar entre Brasil e Paraguai.

Irradiação solar global

Fonte: Dados extraídos do ATLAS do potencial energético solar – ITAIPU

Portanto, considerando os parâmetros citados de insolação diária útil, é defino quanto é

4 RESULTADOS E DISCUSSÕES S9i, A Tabela 1 mostra a relação entre o

Tabela 1: Dados do consumo mensal do minerador

4.1.2 CÁLCULO DO VALOR MENSAL EM nos últimos 4 meses do ano devido as

Tabela 2: Valores tarifários de energia em Campina Grande e Cidade do Leste

4.2 CÁLCULO DA ENERGIA GERADA altera diretamente na geração de energia. Na

Tabela 3: Dados médios de irradiação solar útil

Deste modo, utilizando a equação 3, necessários 20 módulos, e na Cidade do

Tabela 4: Valores e quantidades necessárias de equipamentos para instalação

Com base em comparativos no setor de como a fazenda de mineração foi