Capítulo 1 — Cálculo I | Engenharia Segurança de infraestrutura crítica e setores

Ensino gratuito e aberto para quem deseja proteger pessoas e tecnologias. Não é preciso criar uma conta para estudar ou consultar as respostas.

Mapa da aula

Aula 01 publicada. Mapa, exemplo, EX001.

flowchart TD
  N0["1.1 Por que um engenheiro de infraestrutura precisa"]
  N1["1.2 Derivadas que você vai usar"]
  N2["1.3 Integrais e quantidades acumuladas"]
  N3["1.5 Atividade prática: Laboratório 1"]
  N4["Síntese do capítulo"]
  N0 --> N1 --> N2 --> N3 --> N4
  N4 --> EX["Exemplo resolvido"]
  EX --> Q["Exercício EX001"]

Desenho lógico

flowchart LR
  N0["1.1 Por que um engenheiro de infraestrutura precisa"]
  N1["1.2 Derivadas que você vai usar"]
  N2["1.3 Integrais e quantidades acumuladas"]
  N3["1.5 Atividade prática: Laboratório 1"]
  N4["Síntese do capítulo"]
  N0 --> N1 --> N2 --> N3 --> N4 --> Q["EX001"]

[1.1 Por que um engenheiro de infraestrutura precisa]
    |
    v
[1.2 Derivadas que você vai usar]
    |
    v
[1.3 Integrais e quantidades acumuladas]
    |
    v
[1.5 Atividade prática: Laboratório 1]
    |
    v
[Síntese do capítulo]
    |
    v
[Exemplo resolvido] --> [Exercício EX001]

Objetivos de aprendizagem e conexão

Ao fim deste capítulo você será capaz de (1) ler uma taxa de variação como a grandeza física que um operador acompanha — metros cúbicos por hora saindo de um reservatório, megawatts subindo num alimentador; (2) calcular uma derivada e uma integral definida para as curvas suaves de demanda usadas no planejamento de infraestrutura; (3) resolver problemas do tipo "quanto tempo até…" integrando uma taxa e encontrando o instante em que a quantidade acumulada atinge um limite; e (4) explicar por que o atalho da taxa média pode superestimar o tempo de autonomia. O capítulo inicia o caso de projeto do livro e prepara o Capítulo 6, em que o mesmo reservatório se torna uma equação diferencial com vazão de entrada e controle.

1.1 Por que um engenheiro de infraestrutura precisa de cálculo

Serviços críticos são estoques e fluxos. Um reservatório é um estoque de água tratada; bombas e consumidores são fluxos. Um banco de baterias é um estoque de energia; a carga é um fluxo. O tanque de combustível de um hospital, o empacotamento de uma linha de gás, uma fila de mensagens de telemetria não processadas: cada um é uma quantidade que varia a uma taxa. As perguntas de segurança e de resiliência são quase sempre perguntas sobre essas quantidades sob estresse. Se as bombas param porque sua rede de controle caiu, quantas horas de abastecimento restam? Se um atacante consegue manter uma válvula aberta, com que rapidez o nível sobe? A resposta nunca é um único número tirado de uma folha de dados; é a integral de uma taxa ao longo do tempo, comparada com um limite que importa para as pessoas.

A derivada $dV/dt$ é a taxa instantânea com que uma quantidade $V$ varia. Se $V(t)$ é o volume de um reservatório em metros cúbicos e $t$ está em horas, $dV/dt$ está em m³/h. A integral definida $\int_a^b q(t)\,dt$ de um fluxo $q$ é a quantidade total movimentada entre os instantes $a$ e $b$. O teorema fundamental do cálculo liga as duas: se $V' = -q$, então $V(b) - V(a) = -\int_a^b q(t)\,dt$. Todo cálculo de autonomia neste livro é essa frase aplicada a um estoque real.

1.2 Derivadas que você vai usar

Você precisa de um pequeno conjunto de ferramentas: a derivada de um polinômio ($d/dt\, t^n = n t^{n-1}$), da exponencial ($d/dt\, e^{kt} = k e^{kt}$), do seno e do cosseno ($d/dt \sin(\omega t) = \omega\cos(\omega t)$), as regras do produto e da cadeia, e a ideia de ponto estacionário, em que $f'(t) = 0$. A demanda diária de uma rede de água costuma ser aproximada por uma média mais uma senoide com período de 24 horas. Se $d(t) = 150 + 60\sin(2\pi t/24)$ m³/h, então $d'(t) = 60\cdot(2\pi/24)\cos(2\pi t/24) \approx 15.7\cos(2\pi t/24)$ m³/h², de modo que a demanda sobe mais depressa — cerca de 15,7 m³/h a cada hora — quando o cosseno vale 1. Essa inclinação importa para um controlador de bombas: se as bombas só conseguem acelerar a 10 m³/h², elas ficarão atrás da demanda no seu trecho mais íngreme.

1.3 Integrais e quantidades acumuladas

A antiderivada da demanda acima é $D(t) = 150t - \frac{60\cdot 24}{2\pi}\cos(2\pi t/24) + C$. A água retirada do instante 0 até $T$ é $D(T) - D(0) = 150T - \frac{1440}{2\pi}(\cos(2\pi T/24) - 1)$. Observe o que a senoide faz: ao longo de um dia inteiro ela se cancela, então o total diário é simplesmente $150 \times 24 = 3600$ m³; mas ao longo de parte de um dia ela não se cancela, e quando a interrupção começa muda quanto tempo o estoque dura. Para responder "quando o reservatório atinge seu nível mínimo seguro?", iguale o volume retirado ao volume utilizável e resolva para $T$. A maioria dessas equações não tem solução em forma fechada, então os engenheiros as resolvem numericamente. A bissecção divide ao meio um intervalo que sabidamente contém a resposta até que ele fique pequeno o bastante; é lenta, mas não pode falhar se a função for contínua e o intervalo inicial contiver a raiz.

1.4 Exemplo resolvido: autonomia de Hollin Hill (sintético)

Este caso é sintético. O livro acompanha a inventada Lowmarsh Regional Water Authority (LRWA), que opera a estação de tratamento de água de Fenwick, seis estações de bombeamento (PS1–PS6) e o reservatório de distribuição de Hollin Hill, de 4.000 m³, que abastece, entre outros clientes, o St Brenna's Hospital. A LRWA está iniciando um Telemetry Modernisation and Resilience Programme (TMRP), um programa de modernização e resiliência da telemetria; você faz parte da equipe de engenharia. Suponha que uma falha na rede de telemetria pare as bombas que alimentam Hollin Hill às 06:00 (t = 0), quando o reservatório contém 3.400 m³ (85%). A política de operação trata 20% (800 m³) como o piso abaixo do qual a pressão e as reservas para combate a incêndio ficam em risco. A demanda é a $d(t)$ acima, com t = 0 às 06:00. O volume utilizável é de 2.600 m³. Um planejador que divide pela demanda média obtém 2.600/150 = 17,3 horas. Integrar a curva real dá cerca de 14,6 horas, porque a interrupção começa quando a demanda sobe para o pico da manhã e do meio-dia. A diferença de 2,7 horas não é arredondamento; é a diferença entre despachar a tempo uma equipe com gerador e não despachar.

1.5 Atividade prática: Laboratório 1

Execute o laboratório com a biblioteca padrão do Python 3. Ele encontra o instante exato por bissecção, compara a integração passo a passo pelo método de Euler explícito com três tamanhos de passo e imprime a estimativa pela taxa média.

# Lab 1 - how long does Hollin Hill reservoir last after the pumps stop? (synthetic, offline)
import math
V0, VMIN = 3400.0, 800.0          # m3: start at 85% of 4000, alarm floor at 20%
def demand(t):                     # m3/h, t in hours after the pumps stop (t=0 is 06:00)
    return 150.0 + 60.0 * math.sin(2 * math.pi * t / 24)
def drawn(t):                      # exact integral of demand from 0 to t
    return 150.0 * t - (60.0 * 24 / (2 * math.pi)) * (math.cos(2 * math.pi * t / 24) - 1)
lo, hi = 0.0, 48.0
for _ in range(60):                # bisection on drawn(t) = V0 - VMIN
    mid = (lo + hi) / 2
    lo, hi = (mid, hi) if drawn(mid) < V0 - VMIN else (lo, mid)
print(f"exact time to 20%: {lo:.2f} h")
for dt in (1.0, 0.1, 0.01):        # forward Euler
    v, t = V0, 0.0
    while v > VMIN:
        v -= demand(t) * dt; t += dt
    print(f"Euler dt={dt:<5} reaches 20% at {t:.2f} h")
print(f"constant-average model (150 m3/h): {(V0 - VMIN) / 150:.2f} h")
print(f"dV/dt at t=6 h: {-demand(6):.1f} m3/h")

Saída esperada:

exact time to 20%: 14.62 h
Euler dt=1.0   reaches 20% at 15.00 h
Euler dt=0.1   reaches 20% at 14.70 h
Euler dt=0.01  reaches 20% at 14.63 h
constant-average model (150 m3/h): 17.33 h
dV/dt at t=6 h: -210.0 m3/h

O método de Euler com passo de uma hora superestima o tempo em cerca de 0,4 h; com 0,01 h ele fica a 0,01 h da resposta exata. A derivada em t = 6 h é −210 m³/h: o esvaziamento mais rápido da janela.

Antes de confiar no laboratório, verifique se seus números concordam com o exemplo resolvido: o instante exato deve coincidir com a autonomia integrada de Hollin Hill, e a linha da taxa média deve reproduzir a estimativa mais longa do planejador. Se você mudar o horário de início, mantenha o piso em 800 metros cúbicos para que a comparação continue justa. Um resultado que só aparece com um tamanho de passo é uma propriedade do método, não do reservatório.

1.6 Erros comuns e limites

Usar uma taxa média para um fluxo variável no tempo ao longo de menos de um ciclo completo.
Misturar unidades (litros por segundo contra metros cúbicos por hora: 1 L/s = 3,6 m³/h).
Ignorar que o volume "utilizável" exclui o volume morto e as reservas.
Confiar num método numérico sem compará-lo com um passo menor ou com um caso exato.
Tratar o modelo de demanda como fato. Uma senoide é uma aproximação de planejamento; a demanda real tem efeitos de dia da semana, de clima e de eventos, e a demanda durante uma emergência pode cair (avisos de economia) ou subir (pessoas armazenando água).

Exercícios — capítulo 1

EX001 — Inclinação de pico. Objetivo: usar uma derivada para comparar a demanda com os limites dos equipamentos. Entradas: $d(t)$ da §1.2; bombas que conseguem acelerar no máximo 12 m³/h por hora. Passos: (1) encontre $d'(t)$; (2) encontre seu máximo e o horário em que ocorre, se t = 0 é 06:00; (3) decida se as bombas conseguem acompanhar a demanda em todos os instantes. Evidência esperada: a dedução e uma conclusão de engenharia em uma frase. Rubrica: 3 = derivada, máximo e horário corretos, e conclusão que liga a inclinação ao limite de aceleração; 2 = um deslize; 1 = apenas a derivada; 0 = nada.

EX002 — Autonomia com outro horário de início. Objetivo: mostrar que o horário de início importa. Entradas: Laboratório 1. Passos: altere a demanda para que a interrupção comece às 18:00 (substitua sin(2πt/24) por sin(2π(t+12)/24)), execute de novo e compare com o resultado das 06:00 e com a estimativa pela taxa média. Evidência esperada: três números e uma recomendação de 100 palavras sobre qual valor um plano de emergência deve citar. Rubrica: 3 = números corretos e recomendação de planejar pelo pior horário de início plausível, com o motivo; 2 = números certos, recomendação fraca; 1 = não executado; 0 = nada.

Síntese do capítulo

Infraestrutura é feita de estoques e fluxos. Uma derivada é uma taxa que você pode comparar com o que os equipamentos conseguem fazer; uma integral é uma quantidade que você pode comparar com um limite que protege pessoas. O tempo de autonomia é uma integral resolvida para seu limite superior, geralmente de forma numérica, e depende de quando a perturbação começa. Confira os resultados numéricos contra uma solução exata ou mais fina, mantenha as unidades explícitas e nunca deixe uma média ocupar o lugar de uma curva quando é a curva que a interrupção vai encontrar.