Análise gravimétrica por meio de smartphones e aprimoramentos nas equações do experimento do pêndulo simples como subterfúgio para o ensino de Física Clássica

Gravimetric analysis using smartphones and improvements in the equations of the simple pendulum experiment as a subterfuge for teaching Classical Physics

Resumos

A Gravimetria contemporânea vale-se de gravímetros para encontrar o módulo da aceleração gravitacional (g) de forma pontual na superfície terrestre. Esse trabalho considerou a pertinência da determinação de g no laboratório didático por intermédio do pêndulo, fazendo atenciosa análise dos fenômenos que afetam o dispositivo em oscilação, trazendo contribuições para o ensino da Física Clássica. A primeira, embora modesta, é eficaz, e propõe medir o período de oscilação empregando dois smartphones. A segunda, uma análise teórica acerca do efeito da amplitude angular correlacionada ao arrasto, sobre o período de oscilação do pêndulo. A terceira, refere-se à validação de duas equações que fornecem o valor de g em função da latitude e da altitude, com os dados gravimétricos disponíveis no site do IBGE. No campus Londrina da UTFPR, realizaram-se ensaios independentes, permitindo que se obtivesse o valor $g_{UTFPR-Ld}=(\mathrm{9,7873}\pm \mathrm{0,0012})m\cdot s^{-2}$. Para a validação da metodologia e verificação da acuracidade, as duas equações que fornecem g em função da altitude e da latitude foram confrontadas com dados gravimétricos para definir qual seria mais adequada para gerar o valor de referência. Assim feito, verificou-se que o valor obtido experimentalmente divergia em $\mathrm{0,001}m\cdot s^{-2}$ (0,01%).

Palavras-chave:
Gravimetria; pêndulo; smartphone; ensino de Física

Contemporary gravimetry uses gravimeters to find the magnitude of gravitational acceleration (g) in a punctual manner on the Earth’s surface. This work considered the relevance of determining g in the didactic laboratory using the pendulum, making a careful analysis of the phenomena that affect the oscillating device, bringing contributions to the teaching of Classical Physics. The first, although modest, is effective, and proposes measuring the oscillation period using two smartphones. The second, a theoretical analysis of the effect of angular amplitude correlated to drag, on the period of pendulum oscillation. The third refers to the validation of two equations that provide the value of g like a function of latitude and altitude, with gravimetric data available on the IBGE website. At the Londrina campus of UTFPR, independent tests were carried out, allowing the value $g_{UTFPR-Ld}=(9.7873\pm 0.0012)m\cdot s^{-2}$ to be obtained. To validate the methodology and verify its accuracy, the two equations that provide g as a function of altitude and latitude were compared with gravimetric data to define which would be most suitable for generating the reference value. This done, it was found that the experimentally obtained value differed by $0.001m\cdot s^{-2}$ (0.01%).

Keywords:
Gravimetry; pendulum; smartphone; Physics teaching

1. Introdução

A Gravimetria é o ramo da ciência que se ocupa essencialmente da determinação da aceleração gravitacional, g, de seu gradiente, e do conjunto de procedimentos e técnicas correlatos a estes [¹, ²]. O aprimoramento das metodologias, sobretudo a partir da segunda metade do século XX, alavancou o emprego dos valores gravimétricos em diversos setores. Os dados mais recentes, obtidos por múltiplos meios, valeram-se de levantamentos em terra, no mar, no ar e até mesmo no espaço, por meio de satélites, gerando prodigioso conjunto de dados, permitindo elevar a precisão dos dados gravimétricos. Valores mais confiáveis de g foram empregados em arqueologia, cartografia, ciências ambientais, geofísica, mineralogia, no setor petrolífero, e, também, nas calibrações de precisão para a padronização do quilograma no sistema internacional [³].

O conhecimento da variação temporal da aceleração gravitacional apresenta vasta aplicação tecnológica e científica. O estudo das taxas de variação relativas ao tempo permite o entendimento de detalhes das camadas de gelo polar, do nível do oceano e suas correntes, do ciclo da água no planeta, e da estrutura da Terra. Para ilustrar, entre os anos 2003 e 2006, o satélite GRACE04.16.07 forneceu dados que apontavam para a redução da camada de gelo sobre a Groenlândia. Também possibilitou a detecção de água em aquíferos abaixo da superfície da Terra e mudanças na crosta oriundas daquele terremoto que gerou o tsunami ocorrido no Oceano Índico, em 2004. Considere-se também que a aceleração gravitacional afeta a pressão nos oceanos e a pressão atmosférica, cujos gradientes são empregados para estimativas da velocidade dos ventos e das correntes marítimas profundas, elementos muito importantes para a previsão do tempo. Olhando para o futuro com a mesma audácia dos grandes navegadores lusitanos, na perspectiva de colonizar outros planetas, o conhecimento da aceleração gravitacional nestes novos mundos será fundamental para a análise da viabilidade da vida humana [⁴].

A determinação de g, com precisão da ordem de ${10}^{-9}m\cdot s^{-2}$, contemporaneamente é feita com gravímetros, por meio da queda livre de um corpo, medindo-se a distância percorrida por este em intervalos de tempo extremamente pequenos, modalidade conhecida por determinação absoluta [³]. Embora atualmente em desuso, historicamente a determinação de g iniciou empregando-se o pêndulo, do latim pendulus (pendurado, suspenso), cujo uso para fins diversos é milenar. Na China, no século I A.D., há registro do emprego do dispositivo como sismógrafo primitivo [⁵, ⁶]. Considera-se que Galileu Galilei, por volta de 1602, tenha sido o primeiro cientista a examinar detidamente o pêndulo, constatando seu isocronismo e verificando que, respeitados os limites de precisão de que dispunha, o período de oscilação era proporcional à raiz quadrada do comprimento do fio e independente da massa e da amplitude de oscilação [⁶, ⁷, ⁸]. O pêndulo foi imperioso para a Física Clássica suplantar a Aristotélica. Investindo no empirismo, Galileu empregou o modesto instrumento para medir o tempo e analisar quantitativamente o movimento dos corpos. Inovando, o cientista italiano combina Matemática e Física fundamentada nos ensaios experimentais, inaugurando uma nova forma de fazer ciência, a qual permitia a validação das leis naturais, em contraponto com o pensamento aristotélico que restringia o empirismo. Apoiado no ombro do gigante italiano, Newton desenvolve os pilares da Mecânica. Assim, a Revolução Científica do século XVII, embasada nos paradigmas estabelecidos por Galileu e Newton, esteve plenamente sujeita ao pêndulo. O historiador Richard Westfall chega a afirmar que sem o artefato não haveria o Principia [⁹].

Em 1657, Christiaan Huygens confeccionou o primeiro relógio de pêndulo, aprimorando a cronometragem. Dois anos depois, o cientista vale-se do pêndulo simples para determinar a aceleração gravitacional, estimando-a em $\mathrm{9,5}m.s^{-2}$. Em 1667, aprimorou o cálculo, utilizando-se de sua recém-descoberta equação do período [¹⁰, ¹¹]. Daniel Bernoulli, em 1747, aperfeiçoa o modelo matemático, demonstrando que, em verdade, o período do pêndulo dependia da amplitude angular [¹²]. A partir do início do século XIX, estudos sistemáticos para determinar a aceleração gravitacional foram conduzidos em diversos institutos internacionais, merecendo destaque, em 1818, os trabalhos com o pêndulo reversível desenvolvidos pelo capitão Henry Kater [¹³, ¹⁴, ¹⁵]. De robustos pêndulos demolidores ao impactante pêndulo de Jean Bernard Léon Foucault [¹⁶, ¹⁷], as aplicações científicas, tecnológicas e didáticas multiplicaram-se com o desenrolar dos anos [¹⁸, ¹⁹, ²⁰, ²¹, ²², ²³, ²⁴, ²⁵, ²⁶, ²⁷, ²⁸, ²⁹, ³⁰, ³¹]. No campo do ensino, este dispositivo tão simples agrega múltiplas possibilidades para a construção de subsunçores conceituais da Mecânica Clássica, enquanto oferece belíssima introdução às metodologias investigativas [¹⁸, ¹⁹, ²⁰, ²¹, ²², ²⁴, ³², ³³, ³⁴, ³⁵].

Considerando a abrangência da gravimetria, a importância histórica do pêndulo para a ciência e para a tecnologia, aliada à riqueza didática do dispositivo, entendeu-se que seria válido o presente aporte. Enveredou-se pela trilha focada no emprego de equipamentos amplamente disponíveis e de baixíssimo custo, investindo na análise teórica para suprimir a carência de recursos técnicos. Destarte, a prática com a metodologia proposta pode ser replicada em qualquer instituição de ensino superior. Os trabalhos sobre a determinação da aceleração gravitacional na perspectiva didática são bastante abundantes. Os estudos valem-se de métodos diversos, desde a medida direta do período com cronômetro de acionamento manual, a dispositivos que empregam ultrassom, filmagem e vídeo análise pelo Tracker, ou ainda valendo-se do próprio smartphone como acelerômetro [³⁶, ³⁷, ³⁸, ³⁹, ⁴⁰, ⁴¹]. Entretanto, a prática que se vale do pêndulo apresenta generoso potencial didático que, aliado à simplicidade do dispositivo, torna-se forte argumento a favor de seu emprego, motivando a metodologia aqui desenvolvida,

O objeto de estudo do presente trabalho não está centrado em analisar as amplas aplicações didáticas do pêndulo, contudo, ressalta-se que a proposta permite ao docente a possibilidade de desenvolver, junto aos alunos, o processo de demonstração seguido de verificação. Também poderá ser proveitoso explorá-lo um pouco mais em um viés de demonstração/observação aberta [¹⁹, ²⁰], abrindo espaço para se discutirem conceitos e hipóteses, sendo excelente oportunidade para se desenvolverem os aspectos inerentes da construção do conhecimento científico.

2. Algumas considerações Teóricas acerca do Pêndulo empregado no Cálculo de g

2.1. Pêndulo simples (idealizado)

O pêndulo simples, quando apresentado inicialmente aos estudantes, está idealizado. Constituído por fio muito fino, de massa desprezível, com uma das extremidades presa a um suporte e a outra a um corpo de massa não nula, também com dimensões desprezíveis, de modo que comprimento total do pêndulo coincida com a distância do centro de massa do sistema ao ponto de fixação (L_CM, Figuras 1a e 1b). Hipoteticamente, o dispositivo interagindo exclusivamente com o campo gravitacional uniforme é abandonado ligeiramente fora da posição de equilíbrio, de tal forma que a aproximação $sen\theta =\theta =x/L_{CM}$ não comprometa significativamente a modelagem do problema. Assim, o arco $x=\widehat{BC}=L_{CM}\cdot \theta$ poderia ser considerado um segmento de reta sobre o eixo das abscissas, e o ponto de equilíbrio coincidiria com a origem do eixo cartesiano (x = 0). Neste caso, a componente tangencial da força resultante torna-se $m\cdot a_t=-m\cdot g\cdot \theta$. Com $a_t=\ddot{x}$ [⁴²], chega-se a uma equação diferencial linear de segunda ordem, bastante simples, $\ddot{x}=-{\omega }_0^2\cdot x$, ou, $\ddot{\theta }=-{\omega }_0^2\cdot \theta$, onde ${\omega }_0^2=g_0/L_{CM}$, cuja solução geral é $x(t)=A\cdot \cos ({\omega }_0\cdot t+{\phi }_0)$, sendo $A={\theta }_0\cdot L_{CM}$ a amplitude do movimento e ${\phi }_0$ sua fase inicial, ambos definidas pelas condições iniciais do problema [⁴³]. Como a periodicidade da função cosseno é $2\pi$, implica que o período de oscilação do pêndulo será T, tal que ${\omega }_0\cdot T=2\pi$, resultando em

(1) $$g_0=4{\pi }^2\cdot \frac{L_{CM}}{T^2}$$

Figura 1
O pêndulo simples. (a) Visão tridimensional do pêndulo e suas coordenadas espaciais; (b) Representação das forças que atuam no caso ideal.

2.2. Pêndulo real: considerações cerca do valor do cálculo da aceleração gravitacional

A equação (1) é obtida assumindo-se diversas hipóteses simplificadoras. Contudo, implicitamente, também se considera que o pêndulo oscile no vácuo; que o suporte, onde se prende seja absolutamente rígido e possua inércia extremamente elevada; e ainda, que não haja forças inerciais. Excetuando-se a última conjectura, pois seu efeito não é significativo para a presente proposta [⁴⁴, ⁴⁵, ⁴⁶], as demais hipóteses serão analisadas buscando avaliar o impacto que exercem sobre o cálculo do valor da aceleração. Seguir-se-á uma sequência gradativa quanto à complexidade dos fatores influentes, tratando cada efeito independentemente do outro [⁴⁷, ⁴⁸, ⁴⁹, ⁵⁰].

2.2.1. Empuxo

O efeito do empuxo sobre o pêndulo é relativamente simples de ser analisado e quantificado [⁴²], bastando considerá-lo como um vetor antiparalelo ao peso, atuando no centro de empuxo do sistema. A massa de ar deslocado pelo pêndulo é $m_{ar}={\rho }_{ar}\cdot V_{ar}$, sendo ${\rho }_{ar}$ a densidade e $V_{ar}=V_{fio}+V_{mi\text{ç}}+V_{esf}$ o volume de ar deslocado. Assim concebido, aplicando-se a Segunda Lei de Newton ao sistema, encontra-se $\ddot{x}=-{\omega }^2\cdot x$, com

(2) $${\omega }^2=\left(1-\frac{m_{ar}}{m}\right)\cdot {\omega }_E^2$$

sendo $m=m_{fio}+m_{mi\text{ç}}+m_{esf}$, $\omega =2\pi /T$, ${\omega }_E^2=g_E/L_{CM}$, e o subscrito “E” em g_E indica que este é o valor da aceleração gravitacional obtida considerando-se o efeito do empuxo. Observa-se que a frequência angular do pêndulo na presença do ar é menor que aquela que ocorreria no vácuo, consequentemente, o período de oscilação é maior, mostrando que o empuxo leva o pêndulo a oscilar mais lentamente, e que

(3) $$g_E=g_0\cdot {\left(1-\frac{m_{ar}}{m}\right)}^{-1}$$

sendo $g_E>g_0$, onde g₀ é o valor de g no caso ideal (equação (1)). Aplicando-se a aproximação de primeira ordem $1/(1-\epsilon )=1+\epsilon$, para $\epsilon =m_{ar}/m=m_{ar}/(m_{esf}+m_{mi\text{ç}}+m_{fio})\ll 1$ [⁵¹], tem-se $g_E=g_0+\mathrm{\Delta }{g}_E$, sendo

(4) $$\mathrm{\Delta }{g}_E=g_0\cdot \frac{m_{ar}}{m}$$

2.2.2. Momento de inércia (pêndulo físico)

A próxima correção no cálculo de g considerará que o pêndulo seja um pêndulo físico. O pêndulo empregado no presente estudo possui suas particularidades: constituído por uma esfera de massa m_esf, raio $R=D/2$, conectada a uma esfera pequenina (miçanga) com diâmetro d e massa m_mic, conectada a um fio, de massa m_fio e comprimento L_fio, preso pela outra extremidade à calha de aço (Figura 2). A localização do centro de massa do pêndulo (L_CM, distância do centro de massa ao ponto de fixação do pêndulo) considera que o comprimento total do pêndulo $(L)$ seja dado pela soma $L=L_{fio}+d+D$ (Figuras 2a e 2b), e sua massa total (m) seja dada pela soma das massas dos seus componentes de modo que

(5) $$L_{CM}=\frac{\begin{array}{c}{m}_{fio}\cdot \mathrm{0,5}\cdot L_{fio}+m_{mi\text{ç}}\cdot (L_{fio}+\mathrm{0,5}\cdot d)\hfill \\ +m_{esf}\cdot \left(L_{fio}+d+\mathrm{0,5}\cdot D\right)\hfill \end{array}}{m}$$

Figura 2
(a) Pêndulo constituído por fio (linha de algodão), miçanga e esfera de aço; (b) Detalhe da miçanga que conecta o fio à esfera.

O momento de inércia do pêndulo em relação ao ponto de fixação do fio, no plano do movimento, $I=m_{esf}\cdot L_{CMesf}^{}{}_{}{}^2+(2/5)m_{esf}\cdot R^2+m_{mi\text{ç}}\cdot L_{CMmi\text{ç}}^{}{}_{}{}^2+(1/3)m_{fio}\cdot L_{fio}^{}{}_{}{}^2$[¹⁶, ⁵², ⁵³, ⁵⁴], e o comportamento escalar do torque resultante, para pequenos ângulos, será tal que $I\cdot \ddot{\theta }=-m\cdot g\cdot \theta \cdot L_{CM}$, ou $\ddot{\theta }=-{\omega }^2\cdot \theta$, onde

(6) $${\omega }^2={\omega }_I^2\cdot \frac{m\cdot L_{CM}^{}{}_{}{}^2}{I}$$

com $\omega _I^{}{}_{}{}^2=g_I/L_{CM}$. Uma vez que $m\cdot L_{CM}^{}{}_{}{}^2/I<1$, segue-se que $\omega <{\omega }_0$, mostrando que o período do pêndulo físico também é maior que o do pêndulo ideal. Como $\omega =2\pi /T$, o valor de g corrigido, considerando-se o efeito do momento de inércia será:

(7) $$g_I=g_0\cdot \frac{I}{m\cdot L_{CM}^{}{}_{}{}^2}$$

e definindo-se $g_I=g_0+\mathrm{\Delta }{g}_I$, encontra-se

(8) $$\mathrm{\Delta }{g}_I=g_0\cdot \left(\frac{I}{m\cdot L_{CM}^{}{}_{}{}^2}-1\right)$$

2.2.3. Elasticidade do fio

Admitindo-se a elasticidade do fio, seu comprimento $\mathrm{\ell }(\theta )$, Figura 3, varia durante a oscilação do pêndulo porque o fio de algodão está sujeito à força tensora $\overrightarrow{F}(\theta )$, cujo módulo varia enquanto o pêndulo oscila. O Princípio Fundamental da Dinâmica aplicado ao pêndulo idealizado de comprimento $\mathrm{\ell }(\theta )+R$, sendo R o raio da esfera suspensa, mostra que a componente centrípeta é dada por $m\cdot v^2\cdot {[\mathrm{\ell }(\theta )+R]}^{-1}=F(\theta )-m\cdot g\cdot \cos \theta$. Considerando-se que $\mathrm{\ell }({\theta }_0)={\mathrm{\ell }}_0$ seja o menor comprimento do fio, ocorrendo quando a altura do pêndulo for máxima (h₀, na Figura 3), e que o fio se comporte como uma mola ideal de constante elástica k, de tal modo que $\mathrm{\ell }(\theta )-{\mathrm{\ell }}_0=\mathrm{\Delta }\mathrm{\ell }(\theta )$ seja a variação do seu comprimento. Com $F(\theta )=k\cdot \mathrm{\Delta }\mathrm{\ell }(\theta )$, a energia potencial elástica associada ao fio será $k\cdot \mathrm{\Delta }{\mathrm{\ell }}^2/2$, e tomando-se $h\cong {\mathrm{\ell }}_{M\text{Á}X}-{\mathrm{\ell }}_0\cdot \cos \theta$, a conservação da energia mecânica exigirá que $v^2=(k/m)\cdot (\mathrm{\Delta }\mathrm{\ell }_0^{}{}_{}{}^2-\mathrm{\Delta }{\mathrm{\ell }}^2)+2\cdot g\cdot (h_0-h)$. Então, combinando-se estas relações, encontra-se a variação no comprimento do fio em função do ângulo $\theta (t)$ que o pêndulo forma com a vertical: $\mathrm{\Delta }\mathrm{\ell }=(m\cdot v^2+m\cdot g\cdot {\mathrm{\ell }}_0\cdot \cos \theta )\cdot k^{-1}\cdot {\mathrm{\ell }}^{-1}$. No interior dos parêntesis, o primeiro termo é igual ao dobro da energia cinética, implicando que pode ser substituído pelo dobro da diferença entre a energia mecânica e a potencial, sendo esta última essencialmente gravitacional, de modo que $m\cdot v^2\cong 2\cdot m\cdot g\cdot (h_0-h)$ ou $m\cdot v^2\cong 2\cdot m\cdot g\cdot {\mathrm{\ell }}_0\cdot (\cos \theta -\cos {\theta }_0)$, e assim, $\mathrm{\Delta }\mathrm{\ell }\cong (m\cdot g/k)\cdot (3\cos \theta -2\cos {\theta }_0)$, que permite determinar o comprimento do fio em função de $\theta$:

(9) $$\mathrm{\ell }(\theta )\cong {\mathrm{\ell }}_0+\frac{m\cdot g}{k}\cdot (3\cos \theta -2\cos {\theta }_0)$$

Figura 3
Correlação entre o comprimento do pêndulo e sua altura relativa ao ponto mais baixo.

Expressando-se $\mathrm{\ell }(\theta )$ em função do comprimento do fio quando o pêndulo estiver em repouso no ponto mais baixo da trajetória, $L_{fio}={\mathrm{\ell }}_0+m\cdot g/k$, e ainda, para pequenos ângulos, empregando-se $\cos \theta =1-{\theta }^2/2$, tem-se

(10) $$\mathrm{\ell }(\theta )\cong L_{fio}+\frac{m\cdot g}{K}\cdot \left(\theta _0^{}{}_{}{}^2-\frac{3{\theta }^2}{2}\right)$$

¹de modo que se pode estimar um valor médio para o comprimento do fio em ¼ de ciclo, 0$\le \theta \le {\theta }_0$:

(11) $$L_{m\text{é}dio,fio}=\frac{1}{{\theta }_0}\cdot {\int }_0^{{\theta }_0}\mathrm{}\left[L_{fio}+\frac{m\cdot g}{k}\cdot \left(\theta _0^{}{}_{}{}^2-\frac{3{\theta }^2}{2}\right)\right]\cdot \mathit{d}\theta$$

ou seja, $L_{m\text{é}dio,fio}=L_{fio}+0.5m\cdot g\cdot \theta _0^{}{}_{}{}^2/k$, que levada à equação (2) e posteriormente à equação (1), implica que o valor de g, considerando-se o efeito da elasticidade do fio, será aproximadamente:

(12) $$g_K=g_0\cdot \left(1+\frac{m\cdot g_0}{2k\cdot L_{CM}}\cdot \theta _0^{}{}_{}{}^2\right)$$

ou seja, para $g_k=g_0+\mathrm{\Delta }{g}_k$, encontra-se

(13) $$\mathrm{\Delta }{g}_K=g_o\cdot \mathrm{0,5}\cdot \frac{m\cdot g_0}{k\cdot L_{CM}}\cdot \theta _0^{}{}_{}{}^2$$

sendo que outros autores chegaram a resultados muito próximo ao da equação (13) considerando pequenas perturbações, ou expansões de harmônicos em primeira ordem [⁵⁵, ⁵⁶, ⁵⁷].

2.2.4. Arrasto

O efeito do arrasto, diminuindo a amplitude do movimento a cada ciclo, é facilmente perceptível. A redução suscita dúvidas quanto ao seu efeito sobre o isocronismo do pêndulo, então se deve analisá-las. Primeiramente considere-se que a intensidade da força de arrasto para uma esfera de diâmetro D é dada por $F_{arr}=C\cdot {\rho }_f\cdot A\cdot v^2/2$, onde C é o coeficiente de arrasto, ${\rho }_f$ é a densidade do fluido que a envolve, $A=\pi D^2/4$ é a área de sua seção transversal na direção normal à velocidade relativa ao fluido, $\overrightarrow{v}$ [³⁶, ³⁷, ³⁸]. O coeficiente de arrasto depende do número de Reynolds ($Re={\rho }_f\cdot v\cdot D/\eta$, para um fluido de viscosidade $\eta$), e no intervalo $0<Re<2\cdot {10}^5$, a expressão empírica $C\cong 24/Re+C^{\prime }$ fornece o valor do coeficiente com precisão de 90%. O fator $C^{\prime }$ não é constante, de fato $C^{\prime }=6/(1+\sqrt{R_e})+\mathrm{0,4}$, mas tratá-lo como tal nas condições do presente estudo mostra-se adequado [³⁷, ⁵⁸].

Enquanto o pêndulo oscila, $0\le \theta \le {\theta }_0$, sua velocidade varia de 0 a um máximo, de modo que o mesmo ocorrerá para o número de Reynolds, afetando o arrasto. Desta forma, é útil avaliar-se o intervalo de valores que Re pode ocupar durante um ciclo. Admita-se um pêndulo constituído por uma esfera de diâmetro $D=\mathrm{2,5}cm$, suspensa por um fio de comprimento $L=\mathrm{2,5}m$, oscilando no ar, ${\rho }_{ar}=\mathrm{1,19}kg/m^3$, ${\eta }_{ar}=\mathrm{1,82}\cdot {10}^{-5}\mathrm{}Pa\cdot s$, a $293K$ e altitude $566m$ [³⁶, ³⁷, ³⁸, ⁵⁹], $g=\mathrm{9,79}m\cdot s^{-2}$, com ${\theta }_0={\mathrm{3,50}}^{\circ }$, valores coerentes com os parâmetros desta proposta, de forma que, para esta esfera, encontra-se $0\le Re\le 500$. No intervalo $0\le R_e\le 1$, correspondente às menores velocidades do pêndulo, a intensidade da força de arrasto deve ser predominantemente da forma $F_{arr}\cong -b\cdot v$, com $b=3\cdot \eta \cdot \pi \cdot D$, mas para $R_e\approx {10}^3$, o arrasto será proporcional ao quadrado do módulo da velocidade, $F_{arr}\cong -c\cdot v\cdot |v|$, com $c=C^{\prime }\cdot \pi \cdot D^2\cdot \rho /8$. Assim, durante o movimento, $F_{arr}=-b\cdot v-c\cdot v\cdot |v|$, mas a pequena velocidade da esfera contribuirá para que o arrasto seja preponderantemente linear ($F_{arr}\cong -b\cdot v$), correspondendo à lei de Stokes [⁴⁰, ⁶⁰, ⁶¹].

O módulo do trabalho da força de arrasto durante um quarto de ciclo, $0\le t\le T/4$, ou $0\le \theta \le {\theta }_0$, é $|\tau _{arr}^{}{}_{1/4}{}^{}|={\int }_0^{{\theta }_0}\mathrm{}|F_{arr}|\cdot L_{CM}\cdot \mathit{d}\theta$, e tomando-se a expressão da força de arrasto, $F_{arr}=-b\cdot v-c\cdot v\cdot |v|$, com $\theta ={\theta }_0\cdot \cos \phi$, $\phi ={\omega }_0\cdot t$, $v=L_{CM}\cdot \dot{\theta }$ e integrando-se no intervalo $[0;\pi /2]$ encontra-se o módulo da energia dissipada por ciclo, $|{\tau }_{arr}|\cong 4\cdot |\tau _{arr}^{}{}_{1/4}{}^{}|$, ou seja

(14) $$|{\tau }_{arr}|=\frac{b\cdot \pi }{{\omega }_0}\cdot {({\omega }_0\cdot {\theta }_0\cdot L_{CM})}^2+\frac{8\cdot c}{3{\omega }_0}\cdot {({\omega }_0\cdot {\theta }_0\cdot L_{CM})}^3$$

A energia dissipada pelo arrasto correlaciona-se a uma queda na energia potencial gravitacional dada por $|\mathrm{\Delta }\epsilon |=m\cdot g\cdot L_{CM}\cdot [\cos {\theta }_0-\cos {\theta }_1]$, sendo ${\theta }_1$ o valor da amplitude angular no início do segundo ciclo. Considerando-se $\cos {\theta }_1=1-\theta _1^{}{}_{}{}^2/2$, ${\theta }_0\cong {\theta }_1$, e $\cos {\theta }_0-\cos {\theta }_1\cong {\theta }_0\cdot \mathrm{\Delta }{\theta }_0$, $\mathrm{\Delta }{\theta }_0={\theta }_0-{\theta }_1$, e dividindo-se $|{\tau }_{arr}|$, equação (14), e $|\mathrm{\Delta }\epsilon |$ por $\mathrm{\Delta }t=T=2\pi /{\omega }_0$, e igualando-se as razões encontra-se

(15) $$\begin{array}{l}m\cdot {({\omega }_0\cdot L_{CM})}^2\cdot {\theta }_0\cdot {\displaystyle \frac{\mathrm{\Delta }{\theta }_0}{\mathrm{\Delta }t}}\\ =-{\displaystyle \frac{{\omega }_0}{2\pi }}\cdot [{\displaystyle \frac{b\cdot \pi }{{\omega }_0}}\cdot {({\omega }_0\cdot {\theta }_0\cdot L_{CM})}^2\\ +{\displaystyle \frac{8\cdot c}{3{\omega }_0}}\cdot {({\omega }_0\cdot {\theta }_0\cdot L_{CM})}^3]\end{array}$$

no limite para $\mathrm{\Delta }t\to 0$, $-dt=d{\theta }_0\cdot {[b\cdot {\theta }_0/2m+(4c\cdot {\omega }_0\cdot L_{CM})\cdot \theta _0^{}{}_{}{}^2/3\pi m]}^{-1}$. Integrando-se no intervalo $[0;T]$ com $r=b/(2m)$ e $s=4c\cdot {\omega }_0\cdot L_{CM}/(3\pi m)$ encontra-se

(16) $${\theta }_0(T)={\theta }_{0,m\text{á}x}\cdot \frac{r\cdot e^{-r\cdot T}}{r+s\cdot {\theta }_{0,m\text{á}x}\cdot (1-e^{-r\cdot T})}$$

Apesar da complexidade dos fatores envolvidos no arrasto, a hipótese exigida pela equação (16) é que este contenha duas componentes, uma proporcional à velocidade do pêndulo e outra proporcional ao quadrado desta velocidade, de modo que, posteriormente, os fatores r e s serão obtidos através do ajuste da equação (16) aos dados experimentais. Contudo, espera-se que o arrasto seja predominante linear, $F_{arr}\cong -b\cdot v$, de modo que a equação de movimento do pêndulo corresponderá, com boa aproximação, a $\ddot{\theta }+\gamma \cdot \dot{\theta }+{\omega }_0^2\cdot \theta =0$, sendo $\gamma =b/2m$ e ${\omega }_0^2=g_{arr}/L_{CM}$. A solução geral para o caso subcrítico é

(17) $$\theta (t)={\theta }_0\cdot e^{-\frac{b}{2m}\cdot t}\cdot \cos (\omega \cdot t+{\phi }_0)$$

com $\omega =2\pi /T=\sqrt{{\omega }_0^2-\gamma }<{\omega }_0$, apontando que, também, o arrasto aumenta o período de oscilação [⁴², ⁶²], de modo que o valor corrigido para g devido ao este efeito será:

(18) $$g_{arr}=g_0+\frac{b}{2m}\cdot L_{CM}$$

ou seja, definindo $g_{arr}=g_0+\mathrm{\Delta }{g}_{arr}$, encontra-se

(19) $$\mathrm{\Delta }{g}_{arr}=\frac{b}{2m}\cdot L_{CM}$$

sendo que o fator $b/(2m)$, como dito anteriormente, pode ser determinado pela curva de ajuste da equação (17) ao gráfico de $\theta (t)versust$, obtido experimentalmente.

2.2.5. Amplitude angular

Ao se considerar válida a identidade $sen\theta =\theta$ introduziu-se um erro no cálculo de g, então se faz necessária a respectiva correção. Admitindo que a energia mecânica inicial do pêndulo apresentado na Figura 3, é aproximadamente $m\cdot g\cdot L_{CM}\cdot (1-\cos {\theta }_0)$, tomando-se o nível zero de energia potencial na posição mais baixa do centro de massa pêndulo, sendo θ₀ o ângulo máximo que o fio forma com a vertical, de forma que:

(20) $$\frac{m\cdot v^2}{2}+m\cdot g\cdot L_{CM}\cdot (1-\cos \theta )=m\cdot g\cdot L_{CM}\cdot (1-\cos {\theta }_0)$$

e tomando-se $v=L_{CM}\cdot \dot{\theta }$, é possível escrever a equação diferencial para o tempo em função do ângulo $\theta$: $dt=d\theta /[{\omega }_{\theta }\cdot \sqrt{2(\cos \theta -\cos {\theta }_0)}]$, $\omega _{\theta }^{}{}_{}{}^2=g/L_{CM}$. Integrando esta expressão em um quarto de ciclo, $\theta \in [0;{\theta }_0]$, ou $t\in [0,T/4]$, encontra-se:

(21) $$T=T_0\cdot \frac{2}{\pi }\cdot {\int }_0^{\pi /2}\frac{1}{\sqrt{1-{\delta }^2\cdot {(sen\phi )}^2}}\cdot \mathit{d}\phi$$

Para $T_0=2\pi \sqrt{L_{CM}/g}$, $\delta =sen({\theta }_0/2)$, e $sen(\theta /2)=\delta \cdot sen\phi$. A integral dada na equação (21) é uma integral elíptica do primeiro tipo completa [²⁰, ³⁹, ⁶² pp. 45–49, ⁶³, ⁶⁴, ⁶⁵] e, uma vez solucionada, encontra-se:

(22) $$T=T_0\cdot \left[1+{\left(\frac{1}{2}\delta \right)}^2+{\left(\frac{1\cdot 3}{2\cdot 4}\delta \right)}^4+{\left(\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6}\delta \right)}^6\mathrm{\dots }\right]$$

de modo que, para ${\theta }_0<\mathrm{0,05}rad$, pode-se considerar apenas os dois primeiros termos da série para se obter $T=T_0\cdot \{1+{[sen({\theta }_0/2)]}^2/4\}$, ou ainda $T\cong T_0\cdot \{1+\theta _0^{}{}_{}{}^2/16\}$, de onde

(23) $$g_{{\theta }_0}\cong g_0\cdot \left(1+\frac{\theta _0^{}{}_{}{}^2}{8}\right)$$

Ocorre que, na presença do ar, após n ciclos, a equação (23) deve ser aprimorada considerando a variação da amplitude angular oriunda do arrasto, para não superestimar a correção. A equação (17) mostra que, a cada ciclo, a amplitude angular é aproximadamente reduzida pelo fator de decaimento

(24) $$q=\frac{\theta (T)}{{\theta }_0}=e^{-\frac{b}{2m}\cdot T}$$

Considerando-se que o intervalo de tempo necessário para que a primeira oscilação ocorra seja $T_1=T_0\cdot [1+\theta _0^{}{}_{}{}^2/16]$, no segundo ciclo, como a amplitude reduziu para ${\theta }_2={\theta }_0\cdot q$, o intervalo de tempo necessário para que esta oscilação complete será $T_2=T_0\cdot [1+{({\theta }_0\cdot q)}^2/16]$, e assim, sucessivamente, até a $k-\text{é}sima$ oscilação, onde o período seria dado por $T_k=T_0\cdot [1+{({\theta }_0\cdot q^{k-1})}^2/16]$. Então, o intervalo de tempo necessário para n oscilações completas ocorrerem, $\mathrm{\Delta }{t}_n=T_1+T_2+\mathrm{...}+T_n$, $n\ge 1$, será $\mathrm{\Delta }{t}_n=T_0\cdot \{n+[\theta _0^{}{}_{}{}^2/16]\cdot [1+q^2+q^4+\mathrm{...}+q^{2(n-1)}]\}$, com $\mathrm{\Delta }{t}_n=n\cdot T_n$, encontra-se

(25) $$T_n=2\pi \sqrt{\frac{L}{g_0}}\cdot \left(1+\frac{\theta _0^{}{}_{}{}^2}{16\cdot n}\cdot \frac{1-q^{2n}}{1-q^2}\right)$$

de modo que o valor de g considerando-se a queda na amplitude angular após n ciclos seja:

(26) $$g_{\theta ,n}=g_0\cdot \left(1+\frac{\theta _0^{}{}_{}{}^2}{8n}\cdot \frac{1-q^{2n}}{1-q^2}\right)$$

ou seja, $g_{\theta ,n}\cong g_0+\mathrm{\Delta }{g}_{\theta ,n}$, e

(27) $$\mathrm{\Delta }{g}_{\theta ,n}=g_0\cdot \frac{\theta _0^{}{}_{}{}^2}{8n}\cdot \frac{1-q^{2n}}{1-q^2}$$

3. Valor de Referência para Validação do Resultado Experimental

O estudo precedente forneceu as equações para que o valor da aceleração gravitacional possa ser determinado experimentalmente empregando-se um pêndulo. Contudo, para avaliar a exatidão do valor obtido através destas equações, será necessário um valor de referência $(g_{ref})$. Com este intuito, para se obter g_ref serão analisadas as equações [⁶⁶, ⁶⁷]

(28) $$g1(\lambda ,h)=\mathrm{9,80613}\cdot (1-{\mathrm{2,6325.10}}^{-3}\cdot \cos 2)\cdot (1-{\mathrm{3,92.10}}^{-7}\cdot h)$$

(29) $$\begin{array}{l}g2(\lambda ,h)=\mathrm{9,78032677}\\ \cdot {\displaystyle \frac{[1+\mathrm{0,0052790414}\cdot {(\mathrm{sen}\lambda )}^2+\mathrm{0,0000232718}\cdot {(\mathrm{sen}\lambda )}^4]}{1+{[h/(\mathrm{6,371}\cdot {10}^6)]}^2}}\end{array}$$

as quais fornecem, em unidades do Sistema Internacional, os valores da aceleração em função da latitude e da altitude h. Para se determinar qual destas equações seria a mais adequada para o presente estudo, suas previsões foram confrontadas com os dados gravimétricos das estações da Rede Nacional de Estações Gravimétricas Absolutas (RENEGA) [⁶⁸] e do Relatório das Estações Geodésicas 8061171 do IBGE [⁶⁹, ⁷⁰], cujos dados estão expressos na Tabela 1. No gráfico da Figura 4 estão os valores obtidos com as equações (28) e (29) e os valores gravimétricos absolutos (Tabela 1).

Tabela 1
Valores da aceleração gravitacional (gravimétricos absolutos) em algumas cidades brasileiras.

Cidade	Latitude sul (°)	h (m)	g (m/s2)
1-Teresina (PI)	5,060000	70	9,78016343
2-Brasília (DF)	15,660000	1100	9,78048798
3-Viçosa (RJ)	20,760000	653	9,78460230
4-Valinhos (SP)	23,004167	850	9,78563778
5-Vassouras (MG)	22,400000	400	9,78637581
6-Londrina (PR)	23,305833	539,7	9,7864001
7-Curitiba (PR)	25,454167	910	9,78760387
8-Sta. Maria (RS)	29,720000	85	9,79261636

Figura 4
Valores da aceleração gravitacional: valor gravimétrico absoluto (Tabela 1) e valores obtidos com as equações (28) e (29), para algumas cidades brasileiras.

Objetivando enfatizar a discrepância entre os valores teóricos de g, previstos pelas equações (28) e (29), em relação ao valor gravimétrico absoluto, elaborou-se o gráfico apresentado na Figura 5. A partir desta figura, observa-se que a equação (29) apresenta melhor exatidão na maioria dos casos, contudo, para Londrina, a equação (28) apresenta melhor previsão. Considerando também que, na tabela 1, o valor de g referente a Londrina é fornecido pela estação 8061171 (IBGE), a qual dista somente algo em torno de $\mathrm{2,5}km$ do campus

Figura 5
Diferenças (“erros”) entre os valores obtidos com as equações (28) e (29) e o valor gravimétrico absoluto (Tabela 1).

da UTFPR, entendeu-se que a equação (28) seria mais acertada para a determinação de g_ref para este ensaio. No campus, com $h=568m$ e $={\mathrm{23,3071}}^{\circ }$ [⁶⁹], a equação (28) fornece o valor $\mathrm{9,786226}m\cdot s^{-2}$, e considerando-se que o valor fornecido pela equação (28) divergiu em $\mathrm{0,000074}m\cdot s^{-2}$ do valor gravimétrico absoluto, acrescentou-se este valor a $\mathrm{9,786226}m\cdot s^{-2}$ para se adotar como o valor de referência

$${\mathit{g}}_{\mathit{r}ef}=\mathbf{9,786300}\mathbf{m}\cdot {\mathbf{s}}^{\mathbf{-}\mathbf{2}}$$

para a aceleração gravitacional no campus Londrina da UTFPR.

4. Procedimento Experimental

Foram realizados cinco ensaios independentes com pêndulos distintos que divergiam somente no comprimento por alguns milímetros. Os pêndulos foram confeccionados com um fio de algodão, densidade linear ${\mu }_f=(\mathrm{0,095326}\pm \mathrm{0,000008})g/m$, e uma esfera de aço, diâmetro $D=(\mathrm{25,40}\pm \mathrm{0,05})mm$ e massa $m_{esf}=(\mathrm{67,26}\pm \mathrm{0,03})g$. Para conectar-se o fio à esfera, utilizou-se uma miçanga com diâmetro d e massa m_miç, a qual foi colada à esfera com adesivo instantâneo (“supercola”), sendo $D+d=(\mathrm{29,30}\pm \mathrm{0,05})mm$ e massa $m_{esf}+m_{mi\text{ç}}=(\mathrm{67,26}\pm \mathrm{0,03})g$, (Figura 2). Posteriormente, com auxílio de um ímã de neodímio, fixou-se o pêndulo em uma calha de aço rígida (Figura 6). O comprimento total do pêndulo (L) foi medido com uma trena de aço, distendida desde o ímã até a lâmina nivelada, que resvalava delicadamente a porção inferior da esfera (Figura 6). Para avaliar a calibração da trena empregada, quatro trenas de fabricantes distintos foram esticadas horizontalmente e, até o comprimento de 5,0000 m, a olho nu, não se observaram divergências, sendo a calibração das trenas definida em norma estabelecida [⁷¹]. É importante apontar que se deixou os pêndulos pendurados por cerca de $24h$ para que as fibras dos fios se acomodassem. Surpreendentemente, o comprimento inicial, medido imediatamente à fixação dos pêndulos, era cerca de $\mathrm{1,8}cm$ inferior àquele medido no dia seguinte. Após a acomodação, mesmo depois dos ensaios, não se observaram alterações nos comprimentos, dentro dos limites de precisão das medições com a trena.

Figura 6
Esquema da fixação do pêndulo auxiliada por um ímã de neodímio, e detalhe da lâmina metálica (nivelada) empregada na medição do comprimento total do pêndulo $(L)$.

Uma régua milimetrada foi posicionada horizontalmente atrás do pêndulo (Figura 7a), coincidindo sua posição de equilíbrio com a origem da escala (Figura 7a). Dois smartphones foram utilizados, um operava como cronômetro, o outro filmava o conjunto experimental (Figura 7b). A câmera foi posicionada de forma que não houvesse paralaxe (horizontal ou vertical) em relação à origem (Figuras 8a e 8b). O pêndulo foi posto para oscilar, e a filmagem só foi concluída após 150 ciclos serem completados. A seguir, o vídeo foi projetado quadro a quadro para se determinar o período de oscilação, com incerteza de $\mathrm{0,02}$ s. Objetivando-se otimizar a acuracidade na medição do período, considerou-se a posição de equilíbrio para a determinação de cada ciclo. Evitou-se empregar os pontos extremos, onde a pequena velocidade do pêndulo dificultaria a obtenção precisa do instante em que este pararia e inverteria o sentido do movimento. A medição das amplitudes de oscilação também foi feita por intermédio do vídeo, e considerando o esquema apresentado na Figura 8, onde A_P é o valor da amplitude segundo a leitura direta do vídeo, sujeita a erro de paralaxe, e A é a amplitude corrigida, dada pela relação $A=(h/H)\cdot A_P$.

Figura 7
(a) Posição de equilíbrio do pêndulo coincidindo com a origem da escala; (b) Filmagem do conjunto (pêndulo e cronômetro).

Figura 8
(a) Vista de perfil do pêndulo, da câmera e da régua; (b) Vista superior dos mesmos elementos.

5. Análise de Dados

A seguir serão apresentados e analisados os dados referentes a um dos pêndulos de comprimento total $L=(\mathrm{2,4775}\pm \mathrm{0,0003})m$, e $L_{CM}=(\mathrm{2,4605}\pm \mathrm{0,0005})m$ (equação (5)). Na Figura 9 estão os valores do tempo e das amplitudes de oscilação $(A)$ deste pêndulo a cada 10 ciclos, até que fossem completadas 150 oscilações, $\mathrm{\Delta }{T}_{150}=\mathrm{472,960}s$. A curva de decaimento da amplitude será utilizada posteriormente para se obter o fator de decaimento definido na equação (24). Com os dados apresentados na Figura 9, empregando-se a equação (1), foram gerados os valores de g₀, os quais são apresentados na Figura 10, com os respectivos desvios $(\sigma g_0)$ e diferenças em relação ao valor de referência $(g_0-g_{ref})$. O desvio padrão em g₀, $\sigma g_0=\sqrt{\sigma g_L^{}{}_{}{}^2+\sigma g_T^{}{}_{}{}^2}$, depende da incerteza em L_CM, $\sigma g_L=g_0\cdot {\sigma }_{L_{CM}}/L_{CM}$, a qual permanece constante, $\sigma g_L=\mathrm{2,0}\cdot {10}^{-3}m\cdot s^{-2}$, e também àquela associada a medição do intervalo de tempo, $\sigma g_{Tn}$. Nas primeiras 10 oscilações, $\sigma g_{Tn}$ é $\mathrm{1,8}\cdot {10}^{-2}m\cdot s^{-2}$ e, no último ciclo $(n=150)$, reduz-se a $\mathrm{1,2}\cdot {10}^{-3}m\cdot s^{-2}$.

Figura 9
(a) Valores do intervalo de tempo ($\mathrm{\Delta }{T}_n$) em função do número de ciclos (n); (b) Amplitude de oscilação (A) em função do número de ciclos.

Figura 10
(a) Valores de g₀ versus n (número de ciclos); (b) Comportamento do desvio padrão associado a g₀$(\sigma g_0)$ e do erro $(g_0-g_{ref})$ em função de n.

O motivo de a precisão e a exatidão melhorarem com o incremento de ciclos, conforme se observa na Figura 10, deve-se ao fato de a parcela de incerteza associada à medida do tempo $(\sigma g_T=g_0\cdot 2{\sigma }_{Tn}/T_n)$ estar reduzindo. Sendo $\mathrm{\Delta }{T}_n=⟨\mathrm{\Delta }{T}_n⟩\pm {\sigma }_{\mathrm{\Delta }T}$, segue-se que $T_n=⟨\mathrm{\Delta }{T}_n⟩/n\pm {\sigma }_{\mathrm{\Delta }T}/n$, e para uma incerteza de 0,02 s no início e no término dos n ciclos, ${\sigma }_{\mathrm{\Delta }T}=\mathrm{0,02}\sqrt{2}s$, implica que ${\sigma }_{Tn}=\mathrm{0,02}s\cdot \sqrt{2}/n$, e que $\sigma g_T=(g_0\cdot 2{\sigma }_{\mathrm{\Delta }T}/T_n)/n$ $\cong (\mathrm{0,175}m\cdot s^{-2})/n$. Desta forma, o desvio em g₀ devido à medição do período torna-se menor que a incerteza devida à medida localização do centro de massa $(\sigma g_L)$ após 90 iterações.

O experimentador poderá utilizar manualmente um cronômetro para determinar o intervalo de tempo ($\mathrm{\Delta }{T}_n$). Acaso o faça, a incerteza na medida seria oriunda essencialmente de três elementos: incertezas devidas ao próprio instrumento; ao tempo de reação do operador; e ainda, acerca da avaliação de onde se inicia e onde termina o ciclo por parte do operador. Quando se analisam os vídeos, fica evidente que é bem difícil de definir quando o ciclo se inicia ou termina ao se escolher um dos extremos da elongação como referência. Acaso se opte pelo acionamento manual,

$$\begin{array}{l}\sigma _{\mathrm{\Delta }t}^{}{}_{}{}^2=2\cdot \left(\sigma _{cr}^{}{}_{}{}^2+\sigma _R^{}{}_{}{}^2+\sigma _j^{}{}_{}{}^2\right)\\ \approx 2\cdot \{{[\mathrm{0,01}s]}^2+{(\mathrm{0,02}s)}^2+{(\mathrm{0,02}s)}^2\}\\ \to {\sigma }_{\mathrm{\Delta }t}\approx \mathrm{0,04}s\end{array}$$

[⁷², ⁷³], ou seja, maior que o dobro do que se encontra valendo-se da filmagem.

O valor de g₀ obtido para o ensaio número 1, (Figura 10a, n = 150) considerados os desvios, é

$${\mathit{g}}_{\mathbf{0}\mathrm{,1}(\mathit{n}\mathbf{=}\mathbf{150})}=(\mathbf{9,7706}\mathbf{\pm }\mathbf{0,0023})\mathbf{m}\cdot {\mathbf{s}}^{\mathbf{-}\mathbf{2}},$$

o qual que diverge do valor de referência em $-\mathrm{0,016}$ $m\cdot s^{-2}$, ou seja, está 0,16% abaixo de $g_{ref}=\mathrm{9,786300}m\cdot s^{-2}$. A partir do gráfico mostrado na Figura 10b, é possível observar que, a partir de quarenta ciclos, n = 40, os valores de g₀ não diferem significativamente. Também é possível observar, pelo gráfico, como o desvio padrão reduz conforme o número de oscilações aumenta, assim como o erro relativo ao valor de referência $(g_0-g_{ref})$ também diminui. Note-se que os valores de g₀ sempre estão abaixo do valor de referência ($g_0-g_{ref}<0\to g_0<g_{ref}$), como havia sido previsto.

6. Efetivando as Correções

Todas as correções, detalhadas a seguir, serão referentes a 150 oscilações do ensaio número 1, incluindo-se o orçamento de todos os desvios correlacionados aos dados experimentais. Iniciando pelo empuxo, sendo $V_{fio}=(\mathrm{1,871}\pm \mathrm{0,005})\cdot {10}^{-7}{m}^3$, $V_{mi\text{ç}}=(\mathrm{2,48}\pm \mathrm{0,06})\cdot {10}^{-7}{m}^3$ e $V_{esf}=(\mathrm{8,58}\pm \mathrm{0,02})\cdot {10}^{-6}{m}^3$, assumindo para a densidade do ar ${\rho }_{ar}=(\mathrm{1,1946}\pm \mathrm{0,0003})kg/m^3$, para a pressão atmosférica de $\mathrm{1,0174}\cdot {10}^5\mathrm{}Pa$, com umidade relativa de 45,5%, temperatura de ${\mathrm{22,2}}^{\circ }C$ [⁷⁴, ⁷⁵, ⁷⁶], encontra-se $m_{ar}=(\mathrm{1,0770}\pm \mathrm{0,0029})\cdot {10}^{-5}\mathrm{}kg$. Com a equação (4)

$$\mathrm{\Delta }{g}_E=(\mathrm{0,0015656}\pm \mathrm{0,0000043})m\cdot s^{-2}$$

A correção devida ao momento de inércia do pêndulo considera que $I=(\mathrm{0,40740}\pm \mathrm{0,00012})kg\cdot m^2$ e $m\cdot L_{CM}^2=(\mathrm{0,40692}\pm \mathrm{0,00017})kg\cdot m^2$, substituídos na equação (8), encontra-se

$$\mathrm{\Delta }{g}_I=(\mathrm{0,0114}\pm \mathrm{0,0049})m\cdot s^{-2}$$

Com a linha de algodão utilizada, $k=(\mathrm{21,7}\pm \mathrm{2,8})N/m$, com ${\theta }_0=(\mathrm{0,08106}\pm \mathrm{0,00020})rad$ e $m=(\mathrm{0,0672408}\pm \mathrm{0,0000005})g$, com a equação (13) encontra-se

$$\mathrm{\Delta }{g}_k=(\mathrm{0,000393}\pm \mathrm{0,000040})m\cdot s^{-2}$$

Enfatiza-se que, embora o efeito da elasticidade do suporte tenha sido desconsiderado neste ensaio, acaso o pêndulo fosse fixado em suportes flexíveis, como hastes de aço presas em tripés costumeiramente empregados nos laboratórios didáticos, a magnitude do efeito da flexibilidade está em torno de ${6.10}^{-6}m\cdot s^{-2}$ [⁵⁶, ⁷⁷].

A determinação da correção devida ao arrasto inicia-se com o ajuste das equações (16) e (17) aos dados apresentados na Figura 9b, $A(t)/A_0=\theta (t)/{\theta }_0$, sendo que a equação (16) melhor representa o decaimento porque considera o arrasto de forma mais abrangente (Figura 11a, $R^2=\mathrm{0,9990}$), contudo, o ajuste com a equação (17) fornece

Figura 11
(a) Ajuste pela equação (16) ($R^2=\mathrm{0,9990}$) aos dados experimentais. (b) Ajuste pela equação (17) ($R^2=\mathrm{0,9982}$).

(30) $$\frac{\theta (T)}{{\theta }_0}=\mathrm{0,982965}\cdot e^{-\mathrm{0,00154016}.T}$$

que também oferece uma boa concordância com os dados experimentais, (Figura 11b, $R^2=\mathrm{0,9982}$), apontando que o arrasto seja predominante linear, principalmente nos últimos ciclos, quando a velocidade máxima do pêndulo é menor, conforme se cogitara precedentemente. Assim, o fator de decaimento definido pela equação (24) será $q=\mathrm{0,995513}\pm \mathrm{0,0026}$, para um desvio quadrático médio de 0,026% entre os dados experimentais e a equação (30). Portanto, a correção pelo efeito do arrasto, equação (19) resulta em

$$\mathrm{\Delta }{g}_{arr}=(\mathrm{0,003813}\pm \mathrm{0,000076})m\cdot s^{-2}$$

A última correção refere-se ao efeito da amplitude angular, $\theta (t)$, dada pela equação (27). Para n = 150, ${\theta }_0=(\mathrm{0,081061}\pm \mathrm{0,000068})rad$, encontra-se:

$$\mathrm{\Delta }{g}_{\theta \mathrm{,150}}=(\mathrm{0,00423}\pm \mathrm{0,00041})m\cdot s^{-2}.$$

Todas as correções referentes ao centésimo quinquagésimo ciclo do estudo apresentado estão reunidas na Tabela 2. Constata-se que a correção relacionada ao momento de inércia é a de maior impacto, correspondendo a cerca de 53% de toda a correção. O ângulo de abertura, o arrasto, o empuxo e a elasticidade do fio são responsáveis por de 20%, 18%, 7% e 2% da correção, respectivamente. Deste modo, se o docente considerar somente o efeito do momento de inércia e do empuxo, bastante simples de serem calculados, efetivará cerca de 60% da correção.

Tabela 2
Detalhamento dos fatores de correção para a aceleração gravitacional, n = 150, ensaio n° 1.

Fator de correção	$\mathrm{\Delta }$g (m.s−2)	$\mathrm{\Delta }$g/g0(%)	$\mathrm{\Delta }$g/$\mathrm{\Delta }$gTOTAL(%)
Empuxo	0,001566	0,016	7
Pêndulo físico	0,0114	0,116	53
Elasticidade do fio	0,000393	0,004	2
Arrasto	0,003813	0,039	18
Ângulo (teta)	0,00423	0,043	20
Total	0,021	0,21	100

A correção total ($\mathrm{\Delta }{g}_n=\mathrm{\Delta }{g}_{E,n}+\mathrm{\Delta }{g}_{I,n}+\mathrm{\Delta }{g}_{k,n}+\mathrm{\Delta }{g}_{arr,n}+\mathrm{\Delta }{g}_{\theta \mathrm{,150}}$) apresentada no gráfico da Figura 12, reduz ciclo a ciclo. Isto era esperado porque elevando o número de ciclos, reduz-se a amplitude angular, implicando que o erro na aproximação $sen\theta =\theta$, assumido na equação (1), tornar-se cada vez menor, conforme a correção apresentada na equação (27). Acaso não se levasse o decaimento em consideração, a correção estaria superestimada em $+\mathrm{0,0035}m\cdot s^{-2}$, correspondendo a um erro da ordem de $\mathrm{0,036\%}$ no valor de g. A coerência desta análise fica reforçada quando se observa a Figura 13a, para $n>70$, que o valor de g₀ sofre ligeiro acréscimo, enquanto g tende a permanecer constante, salvo ligeiras flutuações. A Figura 13b detalha as diferenças entre o valor de g_n obtido a cada n ciclos e o valor de referência, mostrando também que, para $n>70$, as diferenças tendem a se estabilizarem, sofrendo flutuações estatisticamente esperadas.

Figura 12
Fatores de correção $\mathrm{\Delta }g$ ciclo a ciclo.

Figura 13
(a) Valores de $g_{o1}n$, (situação ideal) e de g_1,n (valor corrigido para o $n-\text{é}simo$ ciclo do primeiro ensaio apresentado); (b) diferenças percentuais de g_1,n em relação ao valor de referência.

A partir do valor de $g_{\mathrm{0,1}(n=150)}=(\mathrm{9,7706}\pm \mathrm{0,0023})$ $m\cdot s^{-2}$, com todas as correções consideradas (Figura 13) encontra-se o valor da aceleração gravitacional:

$${\mathit{g}}_{\mathbf{1}\mathbf{(}\mathit{n}\mathbf{=}\mathbf{150}\mathbf{)}}\mathbf{=}\left(\mathbf{9,7901}\mathbf{\pm }\mathbf{0,0054}\right)\mathit{m}\mathbf{\cdot }{\mathit{s}}^{\mathbf{-}\mathbf{2}}$$

divergindo do valor de referência em $+\mathrm{0,0038}m\cdot s^{-2}$, ou seja, está 0,039% acima. As incertezas em g são maiores que aquelas em g₀, visto que cada correção carrega a incerteza associada a seus respectivos parâmetros.

Finalmente, seguindo todos os procedimentos adotados para o primeiro ensaio, determinaram-se os valores de g, para os outros quatro ensaios, apresentados na Figura 14. Com os resultados de todos os estudos, encontrou-se a média, com seu respectivo desvio padrão, para o valor de g no campus Londrina da UTFPR:

Figura 14
(a) Valores de g_i para os cinco ensaios $(i=\mathrm{1,2}\mathrm{\dots }5)$, sendo a linha tracejada correspondente ao valor de referência. (b) Diferença percentual (relativa ao valor de referência) para os cinco ensaios realizados.

$${\mathit{g}}_{\mathit{U}\mathit{T}\mathit{F}\mathit{P}\mathit{R}\mathbf{-}\mathit{L}\mathit{d}}\mathbf{=}\left(\mathbf{9,7873}\mathbf{\pm }\mathbf{0,0012}\right)\mathit{m}\mathbf{\cdot }{\mathit{s}}^{\mathbf{-}\mathbf{2}}$$

o qual está $\mathrm{0,001}m\cdot s^{-2}$ abaixo do valor de referência, ou seja, em 0,01%. Vale Reforçar que a maior dificuldade na coleta de dados experimentais e, também, aquela que esteve mais sujeita aos erros humanos, tornando-se o calcanhar de Aquiles do trabalho, foi a medição do comprimento total do pêndulo (L, Figura 2) e, portanto, aquela que potencialmente é a maior responsável pelas flutuações apresentadas nos gráficos da Figura 14. Para ilustrar, um erro de $\mathrm{0,5}mm$ em $L\approx 2500mm$ acarreta um erro $\mathrm{\Delta }{g}_L\approx \mathrm{9,8}m\cdot s^{-2}\cdot \mathrm{0,5}mm/2500mm\approx \mathrm{0,002}m\cdot s^{-2}$, ou seja, 0,02%, o dobro da divergência de $\mathrm{0,001}m\cdot s^{-2}$ entre o valor de g_UTFPR-Ld e o valor de referência.

7. Conclusão

O estudo desenvolvido objetivava aprimorar a acuracidade do valor da aceleração gravitacional obtido por meio de ensaios com o pêndulo simples, no laboratório didático. Para a execução da atividade, empregaram-se aparatos muito simples e amplamente disponíveis. Para a validação do resultado, duas equações que forneciam o valor da aceleração em função da latitude e da altitude foram confrontadas com o valor da aceleração gravitacional mensurada através de gravímetros absolutos, e concluiu-se que o valor de referência para a aceleração gravitacional no Campus Londrina da UTFPR, localidade onde ocorreram os ensaios, seria $g_{ref}=\mathrm{9,786300}m\cdot s^{-2}$. A meta foi atingida, considerando-se que, através da metodologia apresentada, o valor obtido para a aceleração gravitacional foi $g=(\mathrm{9,7873}\pm \mathrm{0,0012})m\cdot s^{-2}$, divergindo em $\mathrm{0,0001}m\cdot s^{-2}$ do valor de referência, ou seja, em 0,01%.

Ficou evidente que, acaso se considere somente o pêndulo idealizado, a aceleração gravitacional obtida tem sua exatidão comprometida por desconsiderar diversos fenômenos suprimidos na simplificação. Entretanto, a boa precisão de g₀, o valor encontrado no caso ideal, foi essencial para que, em se considerando as correções correlacionadas ao empuxo, momento de inércia, elasticidade do fio, arrasto e ângulo de abertura, g fosse obtido com boa acuracidade. A medição cuidadosa do comprimento do fio e utilizar-se da filmagem simultânea do pêndulo oscilando, e do display do smartphone que operava como cronômetro, foi essencial. O procedimento de medição do intervalo de tempo empregando a filmagem elimina os erros humanos associados ao julgamento e eventuais erros do operador, acaso operasse manualmente um cronômetro para definir o intervalo de tempo. Dentre todas as correções, correspondendo a pouco mais de 50%, aquela associada ao momento de inércia é a mais significativa, e o efeito da elasticidade do fio poderá ser desconsiderada acaso se empregue um fio de algodão. O empuxo corresponde a apenas 5% das correções, mas pode ser incorporado por sua simplicidade aliada à oportunidade de incremento na fundamentação teórica. Os demais efeitos serão analisados de acordo com grupo de alunos e da ementa do curso. Para os pêndulos deste trabalho, $L\cong \mathrm{2,5}m$ e esfera de aço com diâmetro $D\cong \mathrm{2,5}cm$, os dados também sugeriram que cerca de oitenta ciclos, algo em torno de quatro minutos, serão suficientes para se obter bons resultados.

O trabalho apresentado norteou-se a fim de que os procedimentos desenvolvidos pudessem ser plenamente replicáveis em qualquer laboratório didático de um curso de graduação em Ciências Exatas. A revisão desta prática tradicional e lúdica gerou procedimentos e análises conceituais que, além de permitirem obter a aceleração gravitacional com boa com boa acuracidade, também ofereceu ricas potencialidades de ensino de Física embasada em um dispositivo extremamente simples. Então, valendo-se desta prática tão popular e versátil, o docente sempre terá uma excelente ferramenta para desenvolver seu trabalho, podendo abordar grande espectro de fundamentos de Física Clássica. Em suma, valer-se do pêndulo como atividade didática é excelente oportunidade para construção do saber científico, reforçando a importância da proposta apresentada.

Referências

• ^[1]
  C.A.C. Castro Junior, Contribuição ao estabelecimento de um sistema gravimétrico para a América do Sul Dissertação de Mestrado, Universidade de São Paulo, São Paulo (2005).
• ^[2]
  R.C. Remache, Ajuste de la red de control básico vertical en funcion de números geopotenciales Tese de Doutorado, Universidad de las Fuerzas Armadas – ESPE, Sangolquí (2010).
• ^[3]
  C.A.C. Castro Junior, G.N. Guimarães e N.C. Ferreira, Geociências 37, 361 (2018).
• ^[4]
  M.R.C. Silva, Ensaio de História 1, 37 (1996).
• ^[5]
  J. Needham, Science and Civilisation in China: Mathematics and the Sciences of the Heavens and the Earth (Cambridge University Press, Cambridge, 1959), v. 3.
• ^[6]
  C. Gauld, em: The Pendulum: Scientific, Historical, Philosophical and Educational Perspectives, editado por M.R. Matthews, C.F. Gauld e A. Stinner (Springer, Dordrecht, 2005).
• ^[7]
  DICIONÁRIO PRIBERAM DA LÍNGUA PORTUGUESA, Pêndulo, disponível em: https://dicionario.pri beram.org/pêndulo, acessado em 11/09/2022.
  » https://dicionario.pri beram.org/pêndulo
• ^[8]
  P. Murdin, Full Meridian of Glory: Perilous Adventures in the Competition to Measure the Earth (Copernicus, New York, 2008).
• ^[9]
  M.R. Matthews, em: The Pendulum: Scientific, Historical, Philosophical and Educational Perspectives, editado por M.R. Matthews, C.F. Gauld e A. Stinner, (Springer, Dordrecht, 2005).
• ^[10]
  A. Stinner, Science & Education 10, 323 (2001).
• ^[11]
  F.L. Silveira e L.O.Q. Peduzzi, Caderno Brasileiro de Ensino de Física 23, 26 (2006).
• ^[12]
  H. Poynting e J.J. Thomson, A Textbook of Physics (Legare Street Press, London, 2022).
• ^[13]
  H. Kater, Philos Trans R Soc Lond. 108, 364 (1818).
• ^[14]
  V.F. Lenzen e R.P. Multihauf, Development of Gravity Pendulums in the 19th Century (Smithsonian Institution, Washington, 1966).
• ^[15]
  HARVARD UNIVERSITY, Reversible (Kater’s) Pendulum, disponível em: https://sciencedemonstrations.fas. harvard.edu/presentations/reversible-katers-pendulum, acessado em 25/11/2022.
  » https://sciencedemonstrations.fas. harvard.edu/presentations/reversible-katers-pendulum
• ^[16]
  H.M. Nussenzveig, Curso de física básica, 1: mecânica (Blucher, São Paulo, 2013), 5 ed.
• ^[17]
  C.H. Furukawa e W. Sano, Revista de Ensino de Física 13, 11 (1991).
• ^[18]
  T.C. Luchese, A. Mergen e R. Andrighetto, Brazilian Journal of Development 6, 4891 (2020).
• ^[19]
  M.S.T. Araújo e M.L.V.S. Abib, Revista Brasileira de Ensino de Física 25, 176 (2003).
• ^[20]
  M.D. Quirino, O período de oscilação do movimento do pêndulo simples e sua relação com seu comprimento e sua amplitude inicial Dissertação de Mestrado, Universidade Federal do Rio de Janeiro, Rio de Janeiro (2020).
• ^[21]
  A.C.M. Silva, M.V.F. Tomaselli e C.P. Barbosa, Arquivos Do MUDI 21, 20 (2017).
• ^[22]
  A.C.P. Fernandes, L.T.S. Auler, J.A.O. Huguenin e W.F. Balthazar, Revista Brasileira de Ensino de Física 38, e3504 (2016).
• ^[23]
  C. Gauld, em: The Pendulum, editado por: M.R. Matthews, C.F. Gauld e A. Stinner (Springer, Dordrecht, 2005).
• ^[24]
  T.D.C. Luchese e M.N. Schneider, Revista Brasileira de Ensino de Física 42, e20200056 (2020).
• ^[25]
  G.B. Santos, D.D.A. Santos, T.D.O. Coghi, J.M. Oliveira, N. Aranha e W. Bonventi, Revista Brasileira de Ensino de Física 42, e0121 (2020).
• ^[26]
  M.M. Almeida, Revista Brasileira de Ensino de Física 43, e20200360 (2021).
• ^[27]
  A.R. Romero Castellanos, C.E. Alvarez-Salazar e P.C. Holanda, Revista Brasileira de Ensino de Física 42, e20190154 (2020).
• ^[28]
  M.C. Agarana e S.A. Iyase, International Journal of Physical Sciences 10, 364 (2015).
• ^[29]
  C.M. Falco, Am J Phys. 44, 733 (1976).
• ^[30]
  A.H. MacDonald e M. Plischke, Phys Rev B. 27, 201 (1983).
• ^[31]
  C.B.C. Silva, B.B. Fagundes, J.D. Ribeiro, A.P. Tonel, F.L. Silveira e P.F. Dorneles, Revista Brasileira de Ensino de Física 42, e20190085 (2020).
• ^[32]
  D.C. Zanotta, E. Cappelletto e M.T. Matsuoka, Revista Brasileira de Ensino de Física 33, 2313 (2011).
• ^[33]
  O.P. Lavor e A.N. Oliveira, REMAT: Revista Eletrônica da Matemática 7, e3001 (2021).
• ^[34]
  A. Original, Caminhos Da Educação Matemática em Revista 11, 230 (2021).
• ^[35]
  K.T. Boni, C.E. Laburú e P.S. Camargo Filho, Investigações em Ensino de Ciências 26, 97 (2021).
• ^[36]
  M.R. Rossini, M.J. Alves, L.H. Amorin e P.S. Camargo Filho, Revista Brasileira de Ensino de Física 42, e20200290 (2020).
• ^[37]
  R.A. Nelson e M.G. Olsson, Am J Phys. 54, 112 (1986).
• ^[38]
  F.J. Arnold, R. Arthur, L.L. Bravo-Roger, M.S. Gonçalves e M.J.G. Oliveira, Revista Brasileira de Ensino de Física 33, 4311 (2011).
• ^[39]
  C.G. Carvalhaes e P. Suppes, Revista Brasileira de Ensino de Física 31, 2701 (2009).
• ^[40]
  W. Bonventi e N.N. Aranh, Revista Brasileira de Ensino de Física 37, 2504 (2015).
• ^[41]
  J. Kuhn e P. Vogt, European J of Physics Education 4, 16 (2013).
• ^[42]
  H.M. Nussenzveig, Curso de Física Básica. Fluidos, Oscilações e Calor (Blucher, São Paulo, 2014), 5 ed.
• ^[43]
  W.E. Boyce, R.C. Diprima e D.B. Meade, Equações Diferenciais Elementares e Problemas de Valores de Contorno (LTC, Rio de Janeiro, 2020), 11 ed.
• ^[44]
  E. Herrera e S. Morett, Revista Brasileira de Ensino de Física 38, e3304 (2016).
• ^[45]
  M. Alonso e E.J. Finn, Física: um Curso Universitário – Mecânica (Blucher, São Paulo, 2014), v. 1, 2 ed.
• ^[46]
  J.F. Price, A Coriolis tutorial, Part 1: 1 the Coriolis force, inertial and geostrophic motion, 2018, disponível em: www.whoi.edu/science/PO/people/jpricejprice@w hoi.edu
  » www.whoi.edu/science/PO/people/jpricejprice@w hoi.edu
• ^[47]
  R. Bellman, Perturbation Techniques in Mathematics, Physics, and Engineering (Dover, New York, 1972).
• ^[48]
  E.J. Routh, The effects of finite amplitude, linear damping, and quadratic dampingare considered as simultaneous perturbations (Dover, New York, 1955), 6 ed.
• ^[49]
  E.J. Routh, A treatise on the dynamics of a system of rigid bodies. Being part II of a treatise on the whole subject (Macmillan, London, 1905).
• ^[50]
  L. Meirovitch, Elements of vibration Analysis (McGraw-Hill, New York, 1975).
• ^[51]
  G.B. Thomas, Cálculo (Pearson Prentice Hall, São Paulo, 2012), v. 2, 12 ed.
• ^[52]
  J.A. Rizcallah, Phys Teach. 53, 482 (2015).
• ^[53]
  S.C. Hong e S.I. Hong, Phys Teach. 51, 139 (2013).
• ^[54]
  J.F. Streib, Am J Phys. 57, 181 (1989).
• ^[55]
  M.G. Olsson, Am J Phys. 44, 1211 (1976).
• ^[56]
  K. Colvill, Phys Educ. 5, 119 (1970).
• ^[57]
  D.D. Novaes, Revista Brasileira de Ensino de Física 35, 1314 (2013).
• ^[58]
  H. Beyer e S. Kempfle, ZAMM – Journal of Applied Mathematics and Mechanics 75, 623 (1995).
• ^[59]
  J. Śmietańska, Kinematic Viscosity of Air Calculator, disponível em: https://www.Omnicalculator.Com/Phy sics/Kinematic-Viscosity-of-Air.
  » https://www.Omnicalculator.Com/Phy sics/Kinematic-Viscosity-of-Air
• ^[60]
  G.D. Quiroga e P.A. Ospina-Henao, Eur J Phys. 38, 065005 (2017).
• ^[61]
  A. Arora, R. Rawat, S. Kaur e P. Arun, arXiv:physics/0608071 (2006).
• ^[62]
  J. Chasnov, em: Scientific Computing (The Hong Kong University of Science and Technology, Hong Kong, 2022).
• ^[63]
  M.L. Vieira, Integrais elípticas – funções hiperbólicas e comprimento de arco, disponível em: https://matematicasimplificada.com/integrais-elipticas-funcoes-hiperbolicas-e-comprimento-de-arco/, acessado em 25/10/2022.
  » https://matematicasimplificada.com/integrais-elipticas-funcoes-hiperbolicas-e-comprimento-de-arco/
• ^[64]
  J. Wang, Q. Xue, L. Li, B. Liu, L. Huang e Y. Chen, Alexandria Engineering Journal 61, 10563 (2022).
• ^[65]
  F.S. Lopes, R.N. Suave e J.A. Nogueira, Revista Brasileira de Ensino de Física 40, e3313 (2018).
• ^[66]
  R. Boynton, em: 60th Annual Conference of the Society of Allied Weight Engineers (Arlington, 2001).
• ^[67]
  W. Lopes, Cad. Bras. Ens. Fís. 25, 561 (2008).
• ^[68]
  C. Gemael, S.R.C. Freitas, P.L. Faggion, J.S. Silva Junior e K. Simões, em: Anais do Simpósio Brasileiro de Geomática (Presidente Prudente, 2002).
• ^[69]
  Mapa Topográfico Londrina, disponível em: https://Pt-Br.Topographic-Map.Com/Map-S24z4/Londrina/?Center=-23.30712%2C-51.11225&zoom=16&popup=-23.30748%2C-51.11414.
  » https://Pt-Br.Topographic-Map.Com/Map-S24z4/Londrina/?Center=-23.30712%2C-51.11225&zoom=16&popup=-23.30748%2C-51.11414
• ^[70]
  IBGE, Relatório de Estação Geodésica, disponível em: http://Www.Bdg.Ibge.Gov.Br/Bdg/Pdf/Relatorio.As p?L1=8061159.
  » http://Www.Bdg.Ibge.Gov.Br/Bdg/Pdf/Relatorio.As p?L1=8061159
• ^[71]
  ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCNICAS, NRB 10123 – Especificação Trena de Fita de Aço (Associação Brasileira de Normas Técnicas, Rio de Janeiro, 2012).
• ^[72]
  K. West, A. Santos, F.J. Almeida e G.M. Dantas, A Física no trânsito: Uma discussão para sala de aula, disponível em: https://editorarealize.com.br/editor a/anais/conedu/2018/TRABALHO_EV117_MD1_SA16_ID10341_14092018013613.pdf
  » https://editorarealize.com.br/editor a/anais/conedu/2018/TRABALHO_EV117_MD1_SA16_ID10341_14092018013613.pdf
• ^[73]
  A. Jain, R. Bansal, A. Kumar e K. Singh, Int J Appl Basic Med Res. 5, 124 (2015).
• ^[74]
  ENGINEERING TOOLBOX, Air – Density, Specific Weight and Thermal Expansion Coefficient vs. Temperature and Pressure, disponível em: https://www.engine eringtoolbox.com/air-density-specific-weight-d_600.h tml, acessado em 11/11/2022.
  » https://www.engine eringtoolbox.com/air-density-specific-weight-d_600.h tml
• ^[75]
  L. Hoyos, Air Density Calculator, disponível em: https://Www.Calctool.Org/Atmospheric-Thermodynamics/Air-Density
  » https://Www.Calctool.Org/Atmospheric-Thermodynamics/Air-Density
• ^[76]
  SIMEPAR. Sistema de Tecnologia e Monitoramento Ambiental do Paraná, Estações Meteorológicas do Simepar no Paraná, disponível em: http://previsao.simepar.br/prognozweb/simepar/dados_estacoes/25425241, acessado em 26/10/2023.
  » http://previsao.simepar.br/prognozweb/simepar/dados_estacoes/25425241
• ^[77]
  O.H.M. Silva, S.D.M. Arruda, C.E. Laburú e E.A.S. Bueno, Caderno Brasileiro de Ensino de Física 30, 409 (2013).

• 1
  Há ainda um detalhe bem interessante revelado pela expressão $\ell \left(\theta \right)\cong L_{fio}+m\cdot g\cdot \left({\theta }_0^2-3{\theta }^2/2\right)/k$: considerando $\theta \left(t\right)={\theta }_0\cdot \cos \omega t$, com ${\left(\cos \omega t\right)}^2=0,5\cdot \left(1+\cos 2\omega t\right)$, encontra-se

  $$\ell \left(\theta \right)\cong L_{fio}+\frac{m.\text{}g}{k}\cdot {\theta }_0^2\left[1-0,75\cdot \cos \left(2\omega .t\right)\right]$$

  revelando que o corpo suspenso oscila radialmente com o dobro da frequência angular do pêndulo.

Datas de Publicação

• Publicação nesta coleção
  26 Jan 2024
• Data do Fascículo
  2024

Histórico

• Recebido
  27 Out 2023
• Aceito
  18 Dez 2023