Acessibilidade / Reportar erro

Descrição temporal de forças de colisão: um modelo didático para laboratório de física assistido por sistema embarcado

Temporal description of colision forces: a didactic model for physics laboratory aided by embedded system

Resumo

Nesse trabalho é modelado um experimento para laboratório de física para descrever temporalmente a força de colisão em uma bola. O modelo emprega um sensor piezoelétrico, um singelo circuito programável (Arduino) ligado a um computador e um suporte para a liberação da bola em queda livre. Um programa em linguagem C foi confeccionado para ler o sinal do sensor e calcular a intensidade da força, bem como para calibrar o sistema. O algoritmo do processamento, empregando as leis da física, é apresentado, assim como o código-fonte. Após a colisão o sistema fornece uma tabela força x tempo com resolução de 120 μs, permitindo o aluno elaborar o gráfico da força ou fazer análise numérica, comprovando as leis da conservação do momento e da energia. Testes de desempenho são discutidos, mostrando que o modelo fornece uma exatidão muito boa, validando assim o método de calibragem. De baixo custo e montagem simples esse sistema embarcado, ou seja, controlado por programa, permite uma exploração minuciosa da força durante a colisão. Pode também ser estendido a outros fenômenos da física, bastando trocar o sensor e ajustar o programa.

Palavras-chave:
força; colisão; sistema embarcado; Arduino; laboratório de física

Abstract

In this study is modelled an experiment for physics laboratory to describe the temporal force in a ball collision. The model employs a piezoeletric sensor, a modest programable circuit (Arduino) connected to a computer and a vertical support to release the ball in free fall. A C-language program was created to read the sensor signal and calculate the force intensity as well as to calibrate the system. The processing algorithm employing the physical laws is presented as the source-code. After the collision a 120 μs resolution force x time table is provided, allowing the student to plot the graphic of force or to do numeric analisys, verifying the momentum and energy conservation laws. Performance tests are discussed showing that the model offers very good accuracy, thus validating the calibration method. At low cost and simple assembly, this embedded system, that is, program-controlled, allows a detailed exploration of force during the colision. As well it can be extended to another phisics effects, merely changing the sensor and ajusting the program.

Keywords:
force; colision; embedded system; Arduino; phisics laboratory

1. Introdução

O estudo de colisões é um tema clássico em experimentos de física. Nesse fenômeno a força de interação começa do zero, atinge um valor máximo e cai a zero novamente. Estudos podem ser concebidos mediante teoria de oscilações e conservação do momento e energia a fim de identificar a forma temporal da força [1[1] M.P. Norton e D.G. Karczub, Fundamentals of Noyse and Vibrations Analysis for Engineers (Cambridge University Press, Cambridge, 2003)., 2[2] D. Halliday, R. Resnick eý J. Walker, Fundamentos de Física (LTC, Rio de Janeiro, 2014), 9ª ed.], F(t). Essa curva pode ser mapeada em laboratório mediante um sensor piezoelétrico, que gera um sinal elétrico quando pressionado, e uma placa Arduino1 1 Nome dado ao dispositivo [3][3] https://www.arduino.cc, acessado em dez. 2017.
https://www.arduino.cc...
, um modesto módulo eletrônico de fácil programação. Ambos baratos e facilmente encontrados no mercado.

Ligado a um computador via cabo USB[4][4] M.A. Cavalcante, C.R.C. Tavolaro e E. Molisani, Rev. Bras. Ensino Fís. 33, 4503 (2011). o Arduino, programado em linguagem Wiring (semelhante à linguagem C++), permite interação com o operador do experimento mediante uma janela de console. Ali são solicitados dados como massa da bola, altura da queda, entre outros, bem como mostra resultados. Uma rotina de calibragem calcula a constante de proporcionalidade entre o sinal elétrico e a força aplicada, gravando o parâmetro na memória EEPROM2 1 Nome dado ao dispositivo do módulo e ali permanecendo mesmo com o sistema desligado. A listagem gerada permite o aluno elaborar o gráfico da força assim como comprovar as leis da mecânica clássica. A seguir serão abordados os requisitos para a compreensão desse trabalho

2. Teoria da conservação da energia e momento linear

A energia mecânica Em em um sistema compreende a energia cinética K e a energia potencial U, ou seja,

(1) E m = K + U

onde

(2) K = m v 2 2 e U = m g h

Em queda livre [5][5] V.R. Coluci, G. Paulino, D.C. Souza e E.P.R. Vasconcelos, Rev. Bras. Ensino Fís. 35, 1 (2013). sua energia potencial gravitacional se converte em cinética. Assim, igualando as Eq. (2) e desprezando as forças dissipativas no ar tem-se

(3) v = 2 g h

e o momento linear p=mv passa a ser

(4) p = m 2 g h

Durante a colisão a força de interação F(t) assume uma forma semelhante à exibida na Figura 1. A bola desce de uma altura h1 e sobe até uma altura menor h2, devido a dissipação de energia no choque.

Figura 1
Perfil aproximado da força de colisão. O detalhe mostra, para uma queda livre, as velocidades inicial e final durante impacto.

Nessa condição a bola recebe um impulso total I[6][6] F.J. Arnold, A.G. Godeny, M.M. Costa, J.F. Viana e R.L. Ximenes, Rev. Bras. Ensino Fís. 40, e2502 (2018). do solo, que equivale à variação da quantidade de movimento mΔv [7][7] F.P. Beer, R.E. Johnston, D.F. Mazurek e E.R. Eisenberg, Mecânica Vetorial para Engenheiros (McGraw-Hill, São Paulo, 2001), 5ª ed., ou seja,

I = m ( v 2 v 1 )

Considerando v2 positivo e v1 negativo tem-se

(5) I = m ( 2 g h 2 + 2 g h 1 ) I = m 2 g ( h 2 + h 1 )

O coeficiente de restituição da bola [7][7] F.P. Beer, R.E. Johnston, D.F. Mazurek e E.R. Eisenberg, Mecânica Vetorial para Engenheiros (McGraw-Hill, São Paulo, 2001), 5ª ed. é

(6) e = p 2 p 1

onde p2 e p1 são respectivamente os momentos lineares da bola imediatamente antes e após o choque. Uma discussão mais detalhada sobre esse coeficiente pode ser vista em [8][8] D.P. Meira Filho, J.K.S. Kamassury e R.C.S. Meira, Rev. Bras. Ensino Fís. 39, e4302 (2017).. Substituindo a Eq. (4) na Eq. (6) [6][6] F.J. Arnold, A.G. Godeny, M.M. Costa, J.F. Viana e R.L. Ximenes, Rev. Bras. Ensino Fís. 40, e2502 (2018).

(7) e = h 2 h 1

A Eq. (5) fica então

(8) I = m 2 g ( e h 1 + h 1 ) I = m 2 g h 1 ( 1 + e )

Esse impulso corresponde à área do gráfico F(t) na Figura 1, ou seja [9][9] J. Santos, B.P.W. Oliveira e O.R. Nelson, Rev. Bras. Ensino Fís. 34, 1 (2012).,

(9) I = F ( t ) d t

Chamando essa área de AF e fazendo uso da Eq. (8), já nomeando h1 de h tem-se

(10) A F = m 2 g h ( 1 + e )

Considerando uma bola caindo de alturas diferentes m passa a ser constante e I na Eq. (8) dependerá apenas da altura h. Portanto

(11) I = α h

onde

(12) α = m 2 g ( 1 + e )

3. Sensor piezoelétrico

Dentre os sensores de força o piezoelétrico figura entre os mais simples. Ele responde apenas a forças dinâmicas [10][10] A. Balbinot e V.J. Brusamarello, Instrumentação e Fundamentos de Medidas (LTC, Rio de Janeiro, 2011), 2ª ed.. Mostrado na Figura 2, na forma de pastilha3 1 Nome dado ao dispositivo , uma camada de cristal piezoelétrico é depositada na região central (branca) da lâmina. Essa região gera uma tensão elétrica quando pressionada.

Figura 2
Sensor Piezoelétrico.

A tensão elétrica gerada é proporcional à força aplicada [11][11] W. B. Lima, Transdutores de deformação à base de sensores piezoelétricos. Dissertação de Mestrado, Universidade Federal Fluminense, Rio de Janeiro (2013).. Assim ela pode ser escrita como

(13) V = λ F

onde λ é uma constante de proporcionalidade e F é a força aplicada no sensor.

Esse tipo de sensor não exige manutenção alguma, não requer fonte de alimentação e é sensível a deformações nanométricas. Responde rápido às deformações na sua rede cristalina, permitindo trabalhar com freqüências da ordem de dezenas de kHz.

Como praticamente não há atraso no sensor a Eq. (13) diz que os gráficos temporais de V e F são proporcionais. Quando a força atinge seu valor máximo a tensão gerada pelo cristal piezoelétrico também será máxima conforme ilustrado na Figura 3.

Figura 3
Perfis da força e da tensão elétrica gerada pelo sensor durante a colisão

4. Plataforma Arduino

Arduino, licenciado pela Ceative Commons, é uma plataforma, de tecnologia aberta, para projetos eletrônicos de pequeno porte. É uma placa programável com um microcontrolador [12[12] F. Pereira, Microcontroladores PIC. Programação em C (Érica, São Paulo, 2002)., 13[13] S. Silveira e M. Girardi, Rev. Bras. Ensino Fís. 39, e4502 (2017)] (Atmega) e várias portas (canais) [14, 15 e 4] para comunicar com outros dispositivos. Sua estrutura pode ser vista na Figura 4.

Figura 4
Arduino. As portas de comunicação estão disponíveis na forma de orifícios nas laterais da placa.

Os dados processados podem ser mostrados no seu console, na forma de uma janela na tela do computador. O diagrama da Figura 5 ilustra as fases de programação e operação desse dispositivo. Para programá-lo basta conectá-lo ao computador via cabo USB e executar o ambiente de desenvolvimento IDE4 1 Nome dado ao dispositivo Arduino, fornecido pelo próprio fabricante. Esse IDE emprega a linguagem Wiring, baseada em C++[16][16] S. Holzner, Programando em C++: um guia prático para a programação profissional (Campus, Rio de Janeiro, 1993)., e é disponível na Internet [3][3] https://www.arduino.cc, acessado em dez. 2017.
https://www.arduino.cc...
. Após ser programado o Arduino pode ser levado ao local onde irá funcionar, necessitando apenas de alimentação elétrica.

Figura 5
(a) Durante a programação o programa é transferido para o Arduino via cabo USB. (b) Em operação os dados processados podem ser exibidos numa janela de console

Quando um mecanismo é controlado por microcontrolador ele é dito sistema embarcado[14][14] S. L. Stevan Junior e R.A. Silva, Automação e instrumentação industrial com Arduino: teoria e projetos (Érica, São Paulo, 2015).. O módulo computacional é destinado a monitorar as propriedades do sistema bem como controlar seus mecanismos.

5. Montagem

A Figura 6 ilustra o suporte usado para a queda da bola. Fabricado em metal, a barra horizontal possui altura ajustável, regulada conforme uma régua anexada à barra vertical. Vibrações na base, decorrentes da própria colisão, são transmitidas ao sensor [10][10] A. Balbinot e V.J. Brusamarello, Instrumentação e Fundamentos de Medidas (LTC, Rio de Janeiro, 2011), 2ª ed., deformando o sinal gerado. Nesse sentido a base deve ser rígida, de massa considerável e estar bem assentada em superfície plana. A base aqui empregada, uma chapa de ferro, possui massa de 1,822 kg e dimensões 8,2 cm x 14,9 cm x 1,9 cm.

Figura 6
Suporte para a queda da bola.

O circuito elétrico está esquematizado na Figura 7. A base do sensor é ligada à porta GND (terra) do Arduino e a região central é ligada à porta A2 (porta analógica) [15] L.A. Dworakowski, A.M. Hartmann, E.M. Kakuno e P.F.T. Dorneles, Rev. Bras. Ensino Fís. 38, e3503 (2016). [14][14] S. L. Stevan Junior e R.A. Silva, Automação e instrumentação industrial com Arduino: teoria e projetos (Érica, São Paulo, 2015).. Quando essa região é pressionada a tensão elétrica gerada é lida pela porta A2, que deve estar aterrada mediante um resistor R, conforme ilustra o circuito elétrico. Isso garante na porta uma tensão de referência V=0. Caso contrário irá manifestar valores arbitrários de leitura, mesmo sem nenhum sinal do sensor. Ao mesmo tempo esse resistor protege a porta de altas correntes, uma vez que ele atua como um divisor de corrente [17][17] R.C. Dorf e J.A. Svoboda, Circuitos Elétricos (LTC, Rio de Janeiro, 2008), 7ª ed..

Figura 7
Circuito elétrico

Percebeu-se, mediante análise em osciloscópio, que o valor de R afeta a intensidade do sinal. Essa relação matemática não foi investigada. Para o tipo de sensor aqui empregado encontrou-se boa resposta para R=10 kΩ.

6. Calibragem

Em funcionamento o sistema deve gerar uma tabela de forças a partir das tensões elétricas lidas pelo Arduino. A força é obtida pela Eq. (13), resultando assim em uma força para cada tensão. Esses dois conjuntos de dados, tensão e força, fornecem dois gráficos e duas áreas de integração, conforme mostrado na Figura 3. Outrossim a Eq. (13) permite calcular a constante que relaciona força e tensão elétrica.

(14) λ = A V A F

onde AV é a área sob a curva no gráfico das tensões e AF é a área sob a curva no gráfico das forças.

Calibrar o aparelho equivale a encontrar o valor de λ. O cálculo dessa constante faz uso da Eq. (14), onde a Eq. (10) fornece um valor teórico para AF. Quanto à área AV ela pode ser obtida por integração numérica na tabela de tensões, ou seja,

(15) A V = i = 0 n = 1 V ( t i ) . ( t i + 1 t i )

Esses cálculos são realizados pelo Arduino. Portanto no processo de calibragem, segundo a Eq. (10), é necessário, antes de soltar a bola, informar a massa m da bola, a altura h da queda e o coeficiente de restituição e da bola.

Uma vez conhecido λ, a Eq. (13) converte o conjunto de tensões {V} no conjunto {F}de forças da seguinte forma

(16) { F } = { V } λ

7. Programação

Os cálculos executados no programa estão diagramados no fluxograma da Figura 8. A cada calibragem a constante λ é recalculada e gravada na memória EEPROM do Arduino[3][3] https://www.arduino.cc, acessado em dez. 2017.
https://www.arduino.cc...
, ali permanecendo, como já mencionado, mesmo após ele desligado. A partir daí, durante as execuções, o programa calcula a força de colisão mediante a Eq. (13), buscando essa constante na memória. Vale comentar que essa memória possui uma vida útil de 100.000 gravações.

Figura 8
Fluxograma do programa no Arduino.

Um material suplementar disponível no final desse documento contém a listagem completa do programa desenvolvido. Nesse material a rotina que faz a leitura do sinal denomina-se loop(), começando na linha 64. Cuidado especial foi tomado no primeiro comando for(). O bloco desse comando rege a taxa de amostragem do sistema. Quanto menos comandos alí houver menor será o tempo entre duas leituras consecutivas. Nesse sentido há apenas duas linhas de comando: uma para ler a tensão e outra para registar o tempo.

A porta analógica do Arduino lê tensões de 0 a 5V numa escala de 0 a 1023. O bloco do segundo comando for faz a conversão da escala 0-1023 para a escala 0-5 [12][12] F. Pereira, Microcontroladores PIC. Programação em C (Érica, São Paulo, 2002)..

A rotina do programa que faz a calibragem, ou seja, calcula λ, denomina-se calculaLambda() e começa na linha 177. A chamada a essa rotina passa como parâmetros m, h e e que são respsctivamente a massa da bola, a altura da queda e o coeficiente de restituição. As equações contempladas nesses cálculos são acompanhados pelo ramo esquerdo da Figura 8. A linha que calcula a variável delta_t divide cada tempo por 106 (linha 190), uma vez que eles são registrados em μs.

Uma captura de tela do console, ilustrada na Figura 9, mostra a interação com o usuário. Os dados solicitados ao operador são digitados no campo ao lado do botão Enviar, dessa forma, não aparecendo na imagem.

Figura 9
Console a caractere do Arduino.

8. Medição

O trabalho executou 5 colisões para 8 alturas diferentes que são 4 cm, 5 cm, 6 cm, 7 cm, 8 cm, 9 cm, 10 cm e 11 cm. Em cada colisão a força foi medida a cada 120 μs, num total de quarenta medidas. A massa da bola é 9,21 g e o diâmetro 25,1 mm. Confeccionada em borracha e apresentando pouca perda de energia na colisão ela é vendida em máquinas de ficha em lanchonetes.

O coeficiente de restituição e, Eq. (7), da bola foi estimado medindo, na régua do próprio aparato (Figura 6), as alturas de descida e subida da bola em queda livre. Foi escolhido h1=27 cm e obtido h2=17 cm resultando em e = 0,9. Esse resultado não é o mesmo para alturas bem maiores. Um método mais preciso para determinar esse coeficiente pode ser obtido mediante [18][18] M.A. Cavalcante, E. Silva, R. Prado e R. Hagg, Rev. Bras. Ensino Fís. 24, 150 (2002)..

A Figura 10 contém a listagem gerada pelo sistema após uma das colisões. É uma tabela força x tempo. Os três primeiros campos, massa da bola, altura da queda, e coeficiente de restituição são informados pelo usuário. O campo Impulso é mostrado para fins didáticos, visto que ele equivale à integral da força, conforme indica a Eq. (9). Essa integral é aqui calculada numericamente como

Figura 10
Captura da tela do console listando o resultado de uma colisão.
(17) I = i = 0 n = 1 F ( t i ) . ( t i + 1 t i )

onde ti e F(ti) correspondem respectivamente às colunas Tempo e Força na listagem. O fator Lambda (calibragem) é buscado na memória EEPROM.

Embora esse estudo considerasse apenas uma determinada bola na análise dos resultados, colisões com outras bolas de tamanhos e massas diferentes também foram efetuadas a título de ilustração. A Figura 11 mostra os gráficos de colisão para bolas diferentes. A tabela tempo x força foi copiada do console para o aplicativo WinPython[19][19] L.E. Borges, Python Para Desenvolveres (Novatec Editora, São Paulo, 2014). que gerou o gráfico. Na legenda m é a massa da bola e h a altura da queda. A curva onde m = 9,21 g corresponde à listagem da Figura 10.

Figura 11
Forças de colisão para três bolas distintas em material e tamanho.

9. Discussão

O impulso I em cada colisão foi calculado e apresentado no gráfico da Figura 12. O cálculo segue a Eq. (17) e o resultado é indicado no cabeçalho da listagem de resultados (campo Impulso na Figura 10).

Figura 12
Impulso recebido pela bola para quedas de 4 cm a 11 cm de altura.

Como a massa da bola é constante o valor teórico do impulso é proporcional à altura da queda h, conforme aponta a Eq. (11). A Eq. (12) permite calcular o valor da constante de proporcionalidade, onde para m = 9,21 g, g = 9,81 m/s2 e e = 0,9 tem-se

α = m 2 g ( 1 + e ) = 0 , 0775

Substituindo esse valor na Eq. (11) chega-se à forma teórica de I como

(18) I = 0 , 0775 h 0 , 5

Como dito anteriormente foram efetuadas cinco colisões para cada uma das alturas estipuladas. Uma regressão potencial [20][20] R. Zentgraf, Estatística Objetiva (ZTG, Rio de Janeiro, 2001). foi realizada com esses valores e foi obtido

(19) I = 0 , 0765 h 0 , 486

As duas curvas, teórica e experimental, dadas respectivamente pelas Eq. (18) e (19), estão representadas no gráfico da Figura 12. A comparaçao entre essas equações manifesta grande semelhança, onde o fator multiplicador e o expoente apresentam erros de 1,29% e 2,8% respectivamente. Essa exatidão é visível também no gráfico, mediante a proximidade entre as duas curvas.

10. Conclusão

O experimento permitiu descrever minuciosamente a força de interação durante a colisão elástica de uma bola, gerando uma listagem tempo x força com uma resolução de 120 μs. Essa tabela permite o aluno elaborar gráficos bem como fazer análise numérica dos dados verificando as leís da mecânica clássica. O modelo oferece notável exatidão, respaldando o método de calibragem criado. A montagem é simples e os componentes, exceto o computador, são de baixo custo. O sensor piezoelétrico custou por volta de R$ 1,00 e a plataforma Arduino, amplamente utilizada e com vasta literatura na Internet, custou aproximadamente R$ 50,00. O suporte de liberação da bola ficou em 100,00, entre material e mão de obra. O algoritmo da calibragem implica conhecimentos de mecânica clássica e cálculo integral. Alunos podem programar a plataforma Arduino mediante outros algoritmos, modificando assim o experimento. O modelo pode ser estendido a outros fenômenos na física, bastando trocar o sensor e ajustar o programa de calibragem.

Supplementary material

O seguinte material suplementar está disponível online:

Material suplementar Programa fonte desenvolvido para a plataforma Arduino.
  • 1
    Nome dado ao dispositivo
  • 2
    Electrically-Erasable Programmable Read-Only Memory
  • 3
    Há vários formatos para esse tipo de sensor. Uma diversidade de licações podem empregar esse dispositivo.
  • 4
    IDE - Integrated Development Environment

Referências

  • [1]
    M.P. Norton e D.G. Karczub, Fundamentals of Noyse and Vibrations Analysis for Engineers (Cambridge University Press, Cambridge, 2003).
  • [2]
    D. Halliday, R. Resnick eý J. Walker, Fundamentos de Física (LTC, Rio de Janeiro, 2014), 9ª ed.
  • [3]
    https://www.arduino.cc, acessado em dez. 2017.
    » https://www.arduino.cc
  • [4]
    M.A. Cavalcante, C.R.C. Tavolaro e E. Molisani, Rev. Bras. Ensino Fís. 33, 4503 (2011).
  • [5]
    V.R. Coluci, G. Paulino, D.C. Souza e E.P.R. Vasconcelos, Rev. Bras. Ensino Fís. 35, 1 (2013).
  • [6]
    F.J. Arnold, A.G. Godeny, M.M. Costa, J.F. Viana e R.L. Ximenes, Rev. Bras. Ensino Fís. 40, e2502 (2018).
  • [7]
    F.P. Beer, R.E. Johnston, D.F. Mazurek e E.R. Eisenberg, Mecânica Vetorial para Engenheiros (McGraw-Hill, São Paulo, 2001), 5ª ed.
  • [8]
    D.P. Meira Filho, J.K.S. Kamassury e R.C.S. Meira, Rev. Bras. Ensino Fís. 39, e4302 (2017).
  • [9]
    J. Santos, B.P.W. Oliveira e O.R. Nelson, Rev. Bras. Ensino Fís. 34, 1 (2012).
  • [10]
    A. Balbinot e V.J. Brusamarello, Instrumentação e Fundamentos de Medidas (LTC, Rio de Janeiro, 2011), 2ª ed.
  • [11]
    W. B. Lima, Transdutores de deformação à base de sensores piezoelétricos Dissertação de Mestrado, Universidade Federal Fluminense, Rio de Janeiro (2013).
  • [12]
    F. Pereira, Microcontroladores PIC. Programação em C (Érica, São Paulo, 2002).
  • [13]
    S. Silveira e M. Girardi, Rev. Bras. Ensino Fís. 39, e4502 (2017)
  • [14]
    S. L. Stevan Junior e R.A. Silva, Automação e instrumentação industrial com Arduino: teoria e projetos (Érica, São Paulo, 2015).
  • [15]
    L.A. Dworakowski, A.M. Hartmann, E.M. Kakuno e P.F.T. Dorneles, Rev. Bras. Ensino Fís. 38, e3503 (2016).
  • [16]
    S. Holzner, Programando em C++: um guia prático para a programação profissional (Campus, Rio de Janeiro, 1993).
  • [17]
    R.C. Dorf e J.A. Svoboda, Circuitos Elétricos (LTC, Rio de Janeiro, 2008), 7ª ed.
  • [18]
    M.A. Cavalcante, E. Silva, R. Prado e R. Hagg, Rev. Bras. Ensino Fís. 24, 150 (2002).
  • [19]
    L.E. Borges, Python Para Desenvolveres (Novatec Editora, São Paulo, 2014).
  • [20]
    R. Zentgraf, Estatística Objetiva (ZTG, Rio de Janeiro, 2001).

Datas de Publicação

  • Publicação nesta coleção
    2019

Histórico

  • Recebido
    30 Jul 2018
  • Revisado
    24 Set 2018
  • Aceito
    12 Nov 2018
Sociedade Brasileira de Física Caixa Postal 66328, 05389-970 São Paulo SP - Brazil - São Paulo - SP - Brazil
E-mail: marcio@sbfisica.org.br