Arrasto aerodinâmico e momento sobre uma tubulação. Uma tubulação de formato cilíndrico com 1 m de diâmetro e 25 m de comprimento está exposta a um vento uniforme de 50 km/h nas condições de atmosfera padrão. Efeitos de extremidade e de rajadas podem ser desprezados. Monte, no computador (em qq software de alto nível), um programa que calcule o
momento fletor nas extremidades devido às forças do vento. Faça cálculos de força de arrasto para diferentes valores de Reynolds (40, 100, 200, 500, 1000, 5000, 10000, 100000) e, em seguida, monte gráficos para analisar o comportamento da força de arrasto em função dos números de Reynolds informados e da força de arrasto em função do coeficiente de arrasto. Ou seja, monte gráficos de: (
e (
).
Dados do tubo:
, em escoamento uniforme com:
; massa específica do ar
; viscosidade absoluta ar
; Efeito de extremidade podem ser desprezados.
momento fletor nas extremidades devido às forças do vento. Faça cálculos de força de arrasto para diferentes valores de Reynolds (40, 100, 200, 500, 1000, 5000, 10000, 100000) e, em seguida, monte gráficos para analisar o comportamento da força de arrasto em função dos números de Reynolds informados e da força de arrasto em função do coeficiente de arrasto. Ou seja, monte gráficos de: (
Passo 1
O enunciado pede que se calcule o momento fletor nas extremidades de uma tubulação cilíndrica causado pelas forças do vento. Pedidos e dados fornecidos:
- Perguntas que o enunciado faz:
- Calcular a força de arrasto
para os números de Reynolds . - Calcular o momento fletor nas extremidades devido a essas forças do vento.
- Montar gráficos de
e .
- Calcular a força de arrasto
- Dados relevantes dados no enunciado:
- Diâmetro da tubulação:
. - Comprimento da tubulação:
. - Velocidade de referência (fornecida pelo enunciado):
. - Pressão ambiente (não necessária diretamente nos cálculos de arrasto):
. - Massa específica do ar:
. - Viscosidade dinâmica do ar:
. - Efeitos de extremidade e rajadas podem ser desprezados.
- Diâmetro da tubulação:
- Observações para resolver (não é pergunta, só informações úteis que vieram no enunciado):
- O enunciado pede cálculo para uma lista de números de Reynolds — isto será tratado como um estudo paramétrico em que cada
determina uma velocidade associada (quando quisermos variar ). Também será calculado o caso real com e o correspondente. - Será necessário obter valores de
(coeficiente de arrasto) em função de ; como é empírico, será usada uma aproximação/ interpolação a partir de dados experimentais típicos para cilindro circular.
- O enunciado pede cálculo para uma lista de números de Reynolds — isto será tratado como um estudo paramétrico em que cada
Passo 2
O que é o número de Reynolds e como usá-lo aqui.
- Definição: o número de Reynolds mede a razão entre forças inerciais e viscosas no escoamento. Para um cilindro de diâmetro
em um escoamento com velocidade característica e viscosidade cinemática , usa-se
- A viscosidade cinemática
relaciona-se com a viscosidade dinâmica e a densidade por
- Vamos calcular
com os dados fornecidos.
- Para um
dado, se quisermos a velocidade associada, invertimos a fórmula:
(Como
,
numericamente.)
Passo 3
A fórmula da força de arrasto e o momento fletor nas extremidades (como aplicá-las na prova).
- Força de arrasto total em um cilindro (projetado perpendicularmente ao vento): a pressão dinâmica atua sobre a área projetada
. A fórmula padrão é
onde
-
é a densidade do fluido ( ), -
é a velocidade do escoamento ( ), -
é a área projetada perpendicular ao escoamento; para cilindro longo ( ), -
é o coeficiente de arrasto (adimensional), função empírica de . - Observação prática: na prova, memorize que para cilindro longo a área projetada por unidade de comprimento é simplesmente
(m), então por metro de comprimento a força por unidade de comprimento é . - Local de aplicação resultante e momento fletor nas extremidades:
- A resultante da força distribuída por
pode ser tratada como uma força concentrada equivalente aplicada no centróide da área projetada (o ponto médio do comprimento), isto é, a do extremo. - Assim o momento fletor nas extremidades devido à força resultante (braço
) é
- A resultante da força distribuída por
- Observação: o enunciado/figura mostra exatamente essa montagem (força resultante atuando ao nível do meio do tubo, braço de alavanca
). Na prova, se pedirem “momento nas extremidades”, verifique sempre o ponto de aplicação resultante; aqui adotamos o modelo usual de força equivalente no ponto médio.
Passo 4
Coeficiente de arrasto
— como obter e como vou proceder no cálculo prático.
- Importante: o
para cilindro circular é experimental e depende de . Não existe uma expressão analítica simples válida para toda a faixa; normalmente recorre-se a curvas experimentais ou tabelas. Na prova pode ser pedido que se use um gráfico fornecido. - Aqui, para fazer os cálculos num programa, vou usar uma tabela aproximada (valores típicos experimentais para um cilindro circular de superfície lisa) e interpolar entre esses pontos. Valores típicos (aproximação representativa) usados neste exercício:
- Observação didática: esses números são aproximações retiradas de curvas experimentais clássicas. No computador vou fazer interpolação logarítmica (interpolar em log10(Re) versus log10(C_D) ou em log10(Re) vs
conforme necessário) para evitar saltos bruscos. Em prova, se houver gráfico dado, leia diretamente no pedido. - Assim, para cada
da lista: calculo , obtenho por interpolação na tabela, calculo com a fórmula da etapa anterior e então .
Passo 5
Programa (em Python) que implementa o procedimento e plota os gráficos pedidos.
- Aqui está um código simples em Python (usar com numpy e matplotlib). O código:
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d
# dados do problema
D = 1.0 # m
L = 25.0 # m
rho = 1.23 # kg/m^3
mu = 1.79e-5 # kg/(m.s)
nu = mu/rho # m^2/s
# Requeridos
Re_list = np.array([40,100,200,500,1000,5000,10000,100000], dtype=float)
# tabela aproximada de Cd (valores experimentais aproximados)
Re_table = np.array([40,100,200,500,1000,5000,10000,100000], dtype=float)
Cd_table = np.array([2.20,1.70,1.40,1.20,1.15,1.10,1.12,1.05], dtype=float)
# interpolador: melhor interpolar em log10(Re) x Cd (a escolha é prática)
logRe_table = np.log10(Re_table)
interp_CD = interp1d(logRe_table, Cd_table, kind='linear', fill_value='extrapolate')
# cálculo paramétrico (cada Re tem sua V correspondente)
V_list = Re_list * nu / D # V = Re * nu / D (D=1)
Cd_list = interp_CD(np.log10(Re_list))
# força de arrasto total e momento na extremidade
A = D * L
q_pref = 0.5 * rho * A # multiplicador para V^2 e C_D
F_D_list = q_pref * V_list**2 * Cd_list
M_list = F_D_list * (L/2)
# caso real com V fixa (50 km/h)
V_fixed = 50.0/3.6 # m/s
Re_fixed = V_fixed * D / nu
Cd_fixed = interp_CD(np.log10(Re_fixed))
F_D_fixed = 0.5 * rho * V_fixed**2 * A * Cd_fixed
M_fixed = F_D_fixed * (L/2)
# imprimir resultados
for Re, V, Cd, Fd, M in zip(Re_list, V_list, Cd_list, F_D_list, M_list):
print(f"Re={Re:.0f}, V={V:.6f} m/s, Cd={Cd:.3f}, F_D={Fd:.6e} N, M_ext={M:.6e} N.m")
print("----- caso real V=50 km/h -----")
print(f"V={V_fixed:.6f} m/s, Re={Re_fixed:.2e}, Cd={Cd_fixed:.3f}, F_D={F_D_fixed:.3f} N, M_ext={M_fixed:.3f} N.m")
# gráficos
plt.figure()
plt.loglog(Re_list, F_D_list, 'o-')
plt.xlabel('Re')
plt.ylabel('F_D (N)')
plt.grid(True, which='both')
plt.title('Re x F_D')
plt.figure()
plt.semilogx(Re_table, Cd_table, 'o-', label='dados')
# plotar interpolação suave
Re_fine = np.logspace(np.log10(Re_table[0]), np.log10(Re_table[-1]), 200)
plt.semilogx(Re_fine, interp_CD(np.log10(Re_fine)), '--', label='interp.')
plt.xlabel('Re')
plt.ylabel('C_D')
plt.grid(True, which='both')
plt.title('Re x C_D')
plt.legend()
plt.show()
- Comentários sobre o código (como usá-lo na prova/TP):
- Troque a tabela de
pelos valores do gráfico experimental que o professor fornecer; o resto do código permanece o mesmo. - O código calcula um caso paramétrico (cada
gera uma ) e também o caso prático com . - Os gráficos gerados: (i)
vs em escala log-log; (ii) vs em escala semilog.
- Troque a tabela de
Passo 6
Cálculos numéricos (aplicando as fórmulas e a tabela de
). Vou apresentar os cálculos principais passo a passo e os resultados numéricos (valores arredondados com 3 algarismos significativos quando faz sentido).
- Parâmetros úteis:
-
. - Área projetada:
. - Prefator para
: . - Logo,
-
-
Para cada
da lista eu uso , pela tabela/interpolação e a fórmula acima. -
Resultados numéricos (paramétrico: cada
gera sua ) — resultados apresentados com notação conveniente: -
:
-
-
:
-
:
-
:
-
:
-
:
-
:
-
:
- Cálculo do caso real com
: -
- Interpolando (pela tabela usada) obtemos
(valor aproximado; o valor exato depende do gráfico experimental exato). - A pressão dinâmica é
. Então
-
- Interpretação prática: esse é o valor que interessa para o caso real do vento de
.
Passo 7
Como usar esse desenvolvimento na prova (dicas finais e como justificar passos)
- Sempre escreva claramente:
- Calcule
. - Se o enunciado der um
e pedir , transforme para (ou, se der , calcule e retire do gráfico). - Leia
no gráfico (ou use tabela/interpolação). - Aplique
. - Para o momento nas extremidades no modelo em que a resultante atua no meio do tubo:
.
- Calcule
- Se precisar de aproximação numérica no caderno, mantenha notação e unidades consistentes e apresente cálculos intermediários (por exemplo, exiba
para cada antes de calcular ). - Se o gráfico experimental for fornecido na prova, cite-o explicitamente e use valores lidos com boa precisão; se não for, explique que adotou uma tabela experimental típica.
Resposta
Re=40, V=5,82\times10^{-4}\ \mathrm{m/s}, C_D=2{,}20, F_D=1{,}15\times10^{-5}\ \mathrm{N}, M_{ext}=1{,}43\times10^{-4}\ \mathrm{N\cdot m}
Re=100, V=1{,}455\times10^{-3}\ \mathrm{m/s}, C_D=1{,}70, F_D=5{,}54\times10^{-5}\ \mathrm{N}, M_{ext}=6{,}92\times10^{-4}\ \mathrm{N\cdot m}
Re=200, V=2{,}911\times10^{-3}\ \mathrm{m/s}, C_D=1{,}40, F_D=1{,}82\times10^{-4}\ \mathrm{N}, M_{ext}=2{,}28\times10^{-3}\ \mathrm{N\cdot m}
Re=500, V=7{,}276\times10^{-3}\ \mathrm{m/s}, C_D=1{,}20, F_D=9{,}77\times10^{-4}\ \mathrm{N}, M_{ext}=1{,}22\times10^{-2}\ \mathrm{N\cdot m}
Re=1000, V=1{,}455\times10^{-2}\ \mathrm{m/s}, C_D=1{,}15, F_D=3{,}75\times10^{-3}\ \mathrm{N}, M_{ext}=4{,}68\times10^{-2}\ \mathrm{N\cdot m}
Re=5000, V=7{,}276\times10^{-2}\ \mathrm{m/s}, C_D=1{,}10, F_D=8{,}96\times10^{-2}\ \mathrm{N}, M_{ext}=1{,}12\ \mathrm{N\cdot m}
Re=10000, V=1{,}455\times10^{-1}\ \mathrm{m/s}, C_D=1{,}12, F_D=3{,}65\times10^{-1}\ \mathrm{N}, M_{ext}=4{,}56\ \mathrm{N\cdot m}
Re=100000, V=1{,}455\ \mathrm{m/s}, C_D=1{,}05, F_D=3{,}42\times10^{1}\ \mathrm{N}, M_{ext}=4{,}27\times10^{2}\ \mathrm{N\cdot m}
Re=100, V=1{,}455\times10^{-3}\ \mathrm{m/s}, C_D=1{,}70, F_D=5{,}54\times10^{-5}\ \mathrm{N}, M_{ext}=6{,}92\times10^{-4}\ \mathrm{N\cdot m}
Re=200, V=2{,}911\times10^{-3}\ \mathrm{m/s}, C_D=1{,}40, F_D=1{,}82\times10^{-4}\ \mathrm{N}, M_{ext}=2{,}28\times10^{-3}\ \mathrm{N\cdot m}
Re=500, V=7{,}276\times10^{-3}\ \mathrm{m/s}, C_D=1{,}20, F_D=9{,}77\times10^{-4}\ \mathrm{N}, M_{ext}=1{,}22\times10^{-2}\ \mathrm{N\cdot m}
Re=1000, V=1{,}455\times10^{-2}\ \mathrm{m/s}, C_D=1{,}15, F_D=3{,}75\times10^{-3}\ \mathrm{N}, M_{ext}=4{,}68\times10^{-2}\ \mathrm{N\cdot m}
Re=5000, V=7{,}276\times10^{-2}\ \mathrm{m/s}, C_D=1{,}10, F_D=8{,}96\times10^{-2}\ \mathrm{N}, M_{ext}=1{,}12\ \mathrm{N\cdot m}
Re=10000, V=1{,}455\times10^{-1}\ \mathrm{m/s}, C_D=1{,}12, F_D=3{,}65\times10^{-1}\ \mathrm{N}, M_{ext}=4{,}56\ \mathrm{N\cdot m}
Re=100000, V=1{,}455\ \mathrm{m/s}, C_D=1{,}05, F_D=3{,}42\times10^{1}\ \mathrm{N}, M_{ext}=4{,}27\times10^{2}\ \mathrm{N\cdot m}
Caso real (V=50\ \mathrm{km/h}=13{,}8889\ \mathrm{m/s}): Re\approx9{,}55\times10^{5}, C_D\approx1{,}03 (interpolado), F_D\approx3{,}04\times10^{3}\ \mathrm{N}, M_{ext}\approx3{,}80\times10^{4}\ \mathrm{N\cdot m}.