A computação numérica como ferramenta para o professor de Física do Ensino Médio

The numerical computation as a tool for high school Physics theacher

NOTAS E DISCUSSÕES

A computação numérica como ferramenta para o professor de Física do Ensino Médio

The numerical computation as a tool for high school Physics theacher

Augusto César de Castro Barbosa^I,1; Cláudio Gonçalves Carvalhaes^I; Marcus Vinicius Tovar Costa^II

^IInstituto de Matemática e Estatística, Universidade do Estado do Rio de Janeiro, Rio de Janeiro, RJ, Brasil

^IIInstituto de Aplicação, Universidade do Estado do Rio de Janeiro, Rio de Janeiro, RJ, Brasil

RESUMO

A proposta deste trabalho é apresentar ao professor de Física do ensino médio uma forma de empregar recursos computacionais para esclarecer e aprofundar conceitos de Física que são explorados de forma limitada por não se poder recorrer ao Cálculo Diferencial e Integral. O pêndulo simples é usado como protótipo para o tratamento sugerido. Mostramos que apenas com o auxílio de uma planilha eletrônica - sem a necessidade, portanto, de conhecimentos avançados de computação - é possível criar um ambiente para o aluno testar o processo de convergência numérica e obter informações do comportamento do sistema a partir de análises gráficas.

Palavras-chave: pêndulo simples, diferenças finitas, planilha eletrônica.

ABSTRACT

The proposal of this work is to present a way to explore computational tools to clarify concepts in the high school Physics classes that are briefly explored since it is not possible to use diferential and integral calculus methods. The simple pendulum is used as a prototype to the suggested treatment. We show that with only the use of a spreadsheet - without requiring advanced computation knowledge - it is possible to create an environment to the student to test the numerical convergence processes and to obtain information about the behavior of the system based on graphical analysis.

Keywords: simple pendulum, finite differences, spreadsheet.

1. Introdução

Uma preocupação do professor de Física do Ensino Médio, em geral, é que no trabalho a ser desenvolvido se busque o entendimento dos fenômenos estudados através da construção de modelos. No entanto, estes levam muitas vezes à uma descrição física complexa. Além disso, uma outra fonte de dificuldade, é que o tratamento matemático necessário para a elaboração de um modelo satisfatório pode estar assentado em um grau de sofisticação muito além daquele usualmente tratado no Ensino Médio [1]. Isto ocorre, por exemplo, no estudo do pêndulo simples, em que o emprego da segunda lei de Newton conduz naturalmente à uma equação diferencial.

A beleza na formulação de um modelo está na busca dos ingredientes fundamentais que, segundo nosso entendimento, devem estar presentes na descrição adequada do fenômeno físico que se deseja estudar [2]. Dessa forma, podemos elaborar duas perguntas: Quais são estes tais ingredientes? Como saber quais são os mais importantes? Em muitos casos, o confronto entre o modelo proposto e a experimentação é a melhor maneira de se proceder. Esta busca na elaboração de um modelo para a descrição de um dado fenômeno não deve deixar de incluir características como ''simplicidade", precisão e poder de previsão [3].

Em geral, o conteúdo de Física apresentado no Ensino Médio passa por uma simplificação que, mesmo limitando as aplicabilidades dos modelos estudados, oferece ao aluno um conjunto de possibilidades que deve ser adequadamente exposto e explorado pelo professor. Nesse sentido, a prática laboratorial ganha uma enorme importância, na medida em que essas práticas possibilitam o aluno perceber com maior clareza os limites dos modelos envolvidos e suas simplificações na descrição dos fenômenos, além de proporcionar a construção dos conceitos básicos.

O computador é um elemento que vem se tornando um poderoso instrumento facilitador no processo ensino/aprendizagem nos últimos anos, principalmente na área científica [1, 4, 5]. Os recursos computacionais têm sido incorporados com sucesso tanto nas práticas de laboratório, como nas demais atividades desenvolvidas em sala de aula [6]. O computador desempenha cada vez mais ativamente um papel de destaque no ensino básico de ciências, principalmente na abordagem de sistemas que apresentam um alto grau de complexidade [7].

Neste trabalho, apresentamos uma forma de se empregar o computador para estudar sistemas físicos que são discutidos de forma restrita no Ensino Médio, devido às estruturas matemáticas envolvidas. Como exemplo, utilizamos o pêndulo simples [8]. Nossa proposta envolverá noções de matemática finita, que podem ser apresentadas com facilidade ao aluno. Quanto ao tratamento computacional, será necessário usar apenas um programa de planilha eletrônica, que faz parte dos pacotes básicos instalados na maioria dos computadores.

2. O pêndulo simples

O pêndulo simples é um dos sistemas protótipos no estudo de oscilações mecânicas, tanto na universidade quanto no Ensino Médio. Além de ter um caráter didático de grande valor, o pêndulo simples é também um ótimo exemplo para se ilustrar os métodos da Física experimental. Isto pode ser observado na estimativa do valor da aceleração da gravidade que pode ser feita facilmente no laboratório didático.

O pêndulo simples é um sistema constituído por uma partícula de massa m suspensa por meio de um fio inextensível de comprimento L e de massa desprezível, como mostra a Fig. 1 [9, 10, 11, 12, 13].

A partícula é afastada de sua posição de equilíbrio (ponto O), de maneira que o fio forme um ângulo q com a vertical, e é então abandonada desta posição. Vamos considerar que a resistência do ar e os efeitos de atrito na articulação entre o fio e o suporte sejam desprezíveis. Neste problema, a inércia é representada pela massa da partícula e a força restauradora é a força gravitacional entre a Terra e a partícula.

As forças que atuam sobre a partícula são o seu peso e a tensão na corda (Fig. 1a). Decompondo-se a força peso em seus componentes radial e tangencial à trajetória (Fig. 1b), a segunda lei de Newton aplicada a este sistema nos fornece:

onde a_r e a_q são as acelerações radial e tangencial, respectivamente, T representa a tensão no fio e g = |g| é a aceleração da gravidade. Neste sistema, não existe variação de velocidade na direção radial, logo a componente da aceleração é nula. A componente tangencial é definida por:

onde x = Lq representa o deslocamento da partícula ao longo do arco. Dessa forma, o sistema de equações toma a forma:

Reescrevemos a segunda equação como:

Neste ponto, em geral, nos livros textos, procura-se calcular a freqüência natural de oscilação, supondo que q seja muito pequeno (q << 1). Apresentamos, no ^apêndice, algumas considerações sobre o movimento do pêndulo para um ângulo qualquer. Dentro dessa aproximação, a Eq. (4) fica:

Esta é uma equação diferencial linear homogênea de segunda ordem, cuja solução geral é:

onde q₀ representa a amplitude máxima de oscilação; a freqüência angular w e o período de oscilação P, independente da massa do pêndulo, são dados por:

Um ponto importante é que o período de oscilação também independe da amplitude q₀, o que é uma conseqüência da aproximação senq » q. Para uma melhor compreensão do erro introduzido por esta aproximação, veja o ^apêndice. Vale mencionar que os livros textos, tanto no Ensino Médio, como no ensino superior, analisam o problema do pêndulo simples a partir dessa aproximação, procurando fazer uma analogia com o problema da mola, deduzindo a Eq. (7) da Lei de Hooke, fazendo:

3. Tratamento computacional

Em diversas áreas da Física os pesquisadores recorrem com freqüência às técnicas de análise numérica quando os problemas matemáticos enfrentados exigem um tratamento analítico intrincado, ou mesmo não permitem um tratamento puramente analítico. De fato, nos anos recentes, nenhuma área da matemática cresceu mais em importância para a Física do que a análise numérica. À medida que os computadores são aperfeiçoados, essa importância aumenta, e aumenta também a rapidez com que a própria análise numérica se desenvolve.

Há anos os métodos numéricos fazem parte do currículo de graduação dos cursos da área tecnológica, onde, normalmente, são oferecidas disciplinas que abordam métodos clássicos para se resolver problemas de Cálculo numericamente.

Nas disciplinas de Cálculo Numérico [14], invariavelmente, o professor apresenta aos alunos métodos para se resolver numericamente integrais e equações diferenciais. Em certos casos, no entanto, é possível se adotar a formulação numérica como via principal para descrever a dinâmica do sistema. Em outras palavras, ao invés de tomar o caminho de descrever a dinâmica analiticamente e, se preciso, resolver as equações resultantes numericamente, às vezes é possível se tratar o problema desde o início empregando-se estruturas numéricas. Esse tratamento evita as equações diferenciais, que são impróprias para o Ensino Médio.

O primeiro passo é considerar que a evolução temporal das variáveis do sistema se dá em intervalos discretos de tempo. Isto significa transformar a variável temporal contínua t Î [0,¥) em uma variável discreta t_i = 0,Dt,2Dt, 3Dt,..., onde cada t_i difere de seus vizinhos mais próximos por um valor fixo Dt,

Esse valor fixo de variação temporal é chamado de passo temporal ou simplesmente passo.

Com o tempo discreto, as demais grandezas dinâmicas também se transformam em variáveis discretas. Por exemplo, a função posição q(t) passa a ser escrita como q_i, sendo q₀ a posição no instante t₀, q₁ a posição no instante t₁, e assim por diante. O mesmo ocorre com as variáveis que medem taxas de variação instantâneas, como velocidade e aceleração, como veremos a seguir.

Na formulação contínua, num instante t¢ qualquer, as derivadas primeira (velocidade instantânea, v_q) e segunda (aceleração instantânea, a_q) são definidas por:

Essas variáveis podem ser transformadas de forma simples em grandezas discretas, tomando-as como valores médios de q e v_q, calculados em intervalos de tempo bastante pequenos (e finitos). Essa forma de definir derivadas é conhecida como Método das Diferenças Finitas (MDF).

Denotando, por conveniência, as derivadas primeira e segunda de q por e , respectivamente, definimos:

e

As definições (12) e (13) podem ser interpretadas como uma aproximação discreta das derivadas contínuas em (10) e (11), em que os seus limites são suprimidos e Dt é muito pequeno. Vale observar aqui que o professor deve utilizar, na apresentação da teoria, a notação padrão adotada nos livros didáticos. Por exemplo, onde encontramos , podemos substituir por a - aceleração angular.

Vamos agora ao problema do pêndulo simples. Seguindo os mesmos passos empregados para se chegar à formulação diferencial da equação do pêndulo linear (Eq. (5)), mas tomando a velocidade instantânea como em (12) e a aceleração como em (13), obtemos que a dinâmica do pêndulo simples linear é descrita por:

Essa é a formulação discreta da equação do pêndulo, obtida via MDF. A Eq. (14) desempenha o mesmo papel que a Eq. (5) desempenhava na formulação contínua. A vantagem é que a Eq. (14) é resolvida para q sem a necessidade de se aplicar conhecimentos matemáticos avançados. Dessa forma, obtemos então:

A relação (15) fornece, recursivamente, os valores de q_i em termos de seus valores nos tempos imediatamente anteriores t_i-1 e t_i-2. O valor q₀ corresponde à amplitude inicial do pêndulo. Se Dt é suficientemente pequeno, q₁ deve ser tomado próximo de q₀, pois não esperamos variações significativas para q num intervalo de tempo pequeno.

Um programa de computador como o OpenOffice ou o Excel é suficiente para se implementar e iterar facilmente a relação (15), obter uma tabela de valores de q_i e o gráfico q×t. Como ilustração, apresentamos a seguir algumas curvas q×t obtidas dessa forma. Para avaliar o resultado da formulação discreta, apresentamos no mesmo par de eixos, o gráfico da solução da formulação contínua fornecida pela relação (6).

A Fig. 2 mostra duas soluções da relação (15) com passos temporais distintos. Na primeira, Dt = 0,03 s e na segunda, Dt = 0,06 s. O eixo vertical representa a amplitude de oscilação em radianos e o horizontal o tempo em segundos. As condições iniciais são q₀ = 0.1 rad e q₁ = 0,1 rad, correspondendo a uma velocidade inicial nula (relação (12)).

Observa-se no gráfico que já no final do primeiro meio ciclo as soluções numéricas analítica são bastante diferentes. Essas diferenças se ampliam à medida que o número de ciclos observados aumenta. Em outras palavras, as soluções numéricas divergem de forma crescente da solução analítica, tanto para Dt = 0,03 s quanto para Dt = 0,06 s, sendo mais acentuado para o Dt maior. Essa discrepância indica que o tratamento numérico é impreciso. É possível melhorar a precisão empregando-se passos temporais que tornem as três curvas indistingüíveis.

Tomando-se valores cada vez menores para Dt como, por exemplo, 0,02, 0,01, pode-se verificar que as curvas referentes a esses valores se aproximam cada vez mais do resultado analítico. Isto evidencia a melhora da precisão numérica. Construindo o gráfico, o leitor constatará que por vários ciclos não haverá diferenças apreciáveis entre as duas formulações. Isto ilustra o papel crítico que a escolha do passo exerce na formulação numérica. Com passos suficientemente pequenos, os valores médios podem ser boas aproximações para as taxas instantâneas.

Vale salientar que há casos em que a não linearidade das equações diferenciais gera instabilidades numéricas que impedem a convergência do MDF, mesmo com passos temporais pequenos. Nesta situação, métodos numéricos mais elaborados ou outras versões do MDF precisam ser considerados.

A planilha eletrônica Excel foi empregada para gerar os gráficos q×t. A planilha Excel ou seu equivalente OpenOffice é um instrumento computacional de fácil aprendizado, onde os cálculos numéricos podem ser executados sem a necessidade de se aprender uma linguagem de programação. Apesar do Excel não ser um software livre, optamos por empregá-lo por ser mais comum encontrá-lo nas escolas do que o OpenOffice.

Os Quadros 1 e 2 abaixo mostram os comandos do Excel que geram o traçado gráfico da relação (15). A célula B₁ armazena o valor de w (fixo) e a célula B₂ o de Dt (também fixo). As colunas C a F irão armazenar os valores de t e da variável angular, nos tempos i-2, i-1 e i. Após o preenchimento das células, como mostradas nos Quadros, deve-se selecionar as células no intervalo C3 a F3, clicando nas células C3 e F3 com a tecla [Shift] pressionada. Após este procedimento deve-se clicar e arrastar a alça de preenchimento da seleção (quadrado escuro no canto direito inferior da célula F3) para estender o cálculo para diversas linhas.

4. Considerações finais

Procuramos mostrar que é possível se empregar o computador para aprofundar o estudo de Física no Ensino Médio, em problemas nos quais o conhecimento de Cálculo Diferencial e Integral é exigido. Esse tratamento propicia ao aluno a possibilidade de estudar um conjunto maior de problemas físicos, como também o introduz, de forma intuitiva, às idéias do Cálculo Numérico e do Cálculo Diferencial e Integral.

Apesar da utilização crescente de computadores como elemento no fornecimento de informações e também como agente auxiliar na construção do conhecimento, existe ainda uma grande resistência dos professores e alunos em utilizá-lo como ferramenta para o desenvolvimento e aprendizagem na área científica. Essa resistência na maioria dos casos se dá em função do desconhecimento de linguagens de computação. A discussão que apresentamos tenta vencer essa resistência, mostrando como se pode usar o computador para resolver um problema importante em Física (o pêndulo simples), sem a necessidade de tais conhecimentos.

Foi mostrado que a determinação do passo temporal adequado à formulação numérica é um exercício de repetição de uma rotina computacional. No entanto, com o devido cuidado, esse exercício pode ser aproveitado para levar o aluno a refletir melhor sobre a diferença entre grandezas médias e instantâneas, preparando-o para no futuro compreender adequadamente o conceito de derivada.

Por fim, queremos ressaltar a importância das atividades laboratoriais para o aluno, já que a Física é uma ciência com profundo caráter experimental. O computador e outras tecnologias são ferramentas auxiliares importantes no processo de ensino/aprendizagem mas, no entanto, não substituem as realizações experimentais.

Recebido em 29/9/2005; Revisado em 2/3/2006; Aceito em 4/4/2006

Quando o fio faz com a vertical um ângulo q₀ qualquer, pode ser provado pelo princípio de conservação de energia [13, 16] que o período de oscilação é dado por:

A Eq. (16) mostra que o período difere cada vez mais da aproximação de pequenas oscilações, à medida que a amplitude de oscilação cresce. Esse resultado pode ser facilmente constatado experimentalmente, com a utilização de dois pêndulos de mesmo comprimento, mas que oscilem com amplitudes diferentes. A Tabela 1 mostra, além dos valores de alguns ângulos selecionados e seus respectivos senos, o erro percentual do cálculo do período considerando a aproximação para pequenos ângulos (Eq. (7)), e a expansão de P(q₀) até o quinto termo (Eq. (16)).

Os dados da Tabela 1 foram obtidos da relação:

A Fig. 5 mostra o gráfico do erro percentual em função do ângulo em graus, destacando o intervalo em que normalmente os experimentos em laboratórios didáticos são realizados.

Este erro percentual, exibido na Tabela 1, pode ser descrito pelo polinômio:

que permite, com boa aproximação, uma avaliação direta do erro que se comete ao se considerar um determinado ângulo q₀Î (0^°,70^°], podendo auxiliar o professor de Ensino Médio na preparação de atividades experimentais envolvendo o pêndulo simples.

• [1] C. Fiolhais e J. Trindade, Rev. Bras. Ens. Fis. 25, 259 (2003).
• [2] I.M. Greca and M.A. Moreira, Science Education 86, 106 (2002).
• [3] I. Halloun, J. Res. Sci. Teach, 33, 1019 (1996).
• [4] R.P. Taylor (ed.), The Computer in the School: Tutor, Too, Tutee (Teachers' College Press, New York, 1980).
• [5] S. Papert, Mind Storms: Children, Computers and Powerful Ideas (Basic Books, New York, 1980).
• [6] E.A. Veit e V.D. Teodoro, Rev. Bras. Ens. Fis. 24, 87 (2002).
• [7] A. Medeiros e C.F. Medeiros, Rev. Bras. Ens. Fis. 24, 77 (2002).
• [8] M. Bazin e P. Lucie, Rev. Bras. Ens. Fis. 3, 3 (1981).
• [9] H. Moysés, Nussenzveig, Curso de Física Básica (Ed. Edgard Blücher, Săo Paulo, 1997), v. 2, 3 ed.
• [10] D. Halliday, R. Resnick e K.S. Krane, Física (LTC, Rio de Janeiro, 1996), v. 2.
• [11] N. Villas Boas, R.H. Doca e G.J. Biscuola, Tópicos de Física (Ed. Saraiva, Săo Paulo, 2003), v. 2.
• [12] L.A. Guimarăes e M. Fonte Boa, Física (Ed. Futura, Săo Pauo, 2004), v. 1.
• [13] C. Kittel, W.D. Knight e M.A. Luderman, Curso de Física de Berkeley (Ed. Blücher, Săo Paulo, 1970), v. 1.
• [14] D.M. Cláudio e J.M. Marins, Cálculo Numérico Computacional (Ed. Atlas S.A., Săo Paulo, 1989).
• [15] J.S. Figueira e E.A. Veit, Rev. Bras. Ens. Fis. 26, 203 (2004).
• [16] C. Medina, S. Velazco e J. Salinas, Rev. Bras. Ens. Fis. 24, 254 (2002).

apêndice

Datas de Publicação

• Publicação nesta coleção
  16 Out 2006
• Data do Fascículo
  Jun 2006

Histórico

• Aceito
  04 Abr 2006
• Revisado
  02 Mar 2006
• Recebido
  29 Set 2005