Série de Taylor: a aproximação de pequenos ângulos escondida dentro da análise de segunda ordem

A série de Taylor que você conheceu em Cálculo II costuma chegar como uma forma de aproximar uma função difícil por um polinômio: escreva sin, ou cos, ou 1/(1 menos x) como uma soma infinita de derivadas, e pare quando estiver perto o suficiente. Parece uma ferramenta para gráficos e limites de erro, longe de um pilar ou de uma chapa de base. Não é. Toda vez que um engenheiro escreve o equilíbrio na forma indeformada de uma estrutura, ele está silenciosamente truncando uma série de Taylor depois do primeiro termo, e toda vez que roda uma análise de segunda ordem, está mantendo o termo seguinte.

Essa é toda a promessa deste artigo. Pegamos um pilar de aço real, o inclinamos sob carga, e mostramos que a escada inteira, da primeira ordem à segunda ordem à análise completa de grandes deslocamentos, é uma única expansão de Taylor, lida em três profundidades de truncagem diferentes.

• A análise de primeira ordem mantém apenas o primeiro termo: sin theta é theta, cos theta é 1, e a geometria nunca se move.
• A análise de segunda ordem (P-Delta) mantém o termo seguinte: a carga axial agindo através da flecha lateral, o momento P vezes Delta.
• O famoso amplificador 1/(1 menos P/Pcr) revela-se uma série geométrica, um termo de Taylor por passagem da carga pela flecha crescente.

Todo número aqui bate com o motor FEM do CalcSteel até o decimal. É esse o sentido de ter um motor real por trás das palavras: a conta à mão e a solução da máquina têm que se encontrar, e se encontram. Tudo roda grátis no navegador, então você pode reconstruir o pilar por conta própria enquanto lê.

Três leitores tiram coisas diferentes disto. O estudante de engenharia ganha um motivo concreto de a série de Taylor merecer seu lugar na ementa de Cálculo I e II. O calculista atuante ganha um modelo mental limpo de onde o P-Delta realmente vem. O calculista autônomo ganha uma forma gratuita, nativa do navegador, de ler uma carga crítica Pcr e formar o amplificador sem comprar assento de nada.

Uma série de Taylor reconstrói uma função suave perto de um ponto a a partir de suas derivadas naquele único ponto. Escreva-a como f(x) = soma em n de f⁽ⁿ⁾(a)/n! vezes (x menos a)ⁿ: o valor, mais a primeira derivada vezes (x menos a), mais a segunda derivada sobre 2! vezes (x menos a)², e assim por diante. Quando o ponto-base é a = 0 a série ganha nome próprio, a série de Maclaurin, e é essa a versão que usaremos para sin, cos e 1/(1 menos x).

A imagem a guardar é esta: cada termo acrescentado abraça a curva verdadeira sobre um intervalo mais largo. Mantenha um termo e você tem uma reta horizontal tocando a curva em a; some o seguinte e você tem a tangente; some o próximo e o polinômio começa a curvar junto com a função. Parar cedo dá uma aproximação deliberada, não um erro.

O que diz o quão longe você está é o resto, R_n. A forma de Lagrange limita o erro de parar depois de n termos pelo tamanho do primeiro termo que você largou, avaliado em algum ponto do intervalo. Esse resto é a contabilidade honesta de uma truncagem: diz exatamente quando manter um termo é seguro e quando não é. Segure-o, porque o mesmo resto vai depois avisar um engenheiro de estruturas de quando uma aproximação de pequenos ângulos deixa de ser inofensiva.

Esta é a sequência natural do Cálculo I. Lá você conheceu a derivada ordinária, o único número f'(a). A série de Taylor simplesmente empilha todas as derivadas em a numa só expressão que reconstrói a função em torno daquele ponto. Nada de novo é pedido de você, apenas que guarde mais das derivadas que já sabe tomar.

Aponte a máquina de Maclaurin para as funções trigonométricas e as famosas regras de pequenos ângulos caem prontas. As séries são:

• sin θ = θ − θ³/6 + θ&sup5;/120 − ..., então para θ pequeno o primeiro termo domina e sin θ ≈ θ.
• cos θ = 1 − θ²/2 + θ⁴/24 − ..., então cos θ ≈ 1 (ou, um termo melhor, 1 − θ²/2).
• tan θ = θ + θ³/3 + 2θ⁵/15 + ..., então tan θ ≈ θ.

Estas não são regras de bolso tiradas de uma tabela. São séries de Taylor com todo termo depois do primeiro jogado fora, e o ângulo θ está em radianos, que é a única razão de o coeficiente principal ser um 1 limpo.

Quão boa é a truncagem? Aqui está o erro de largar tudo além do primeiro termo, em porcentagem, nas rotações que um pórtico real de fato vê:

Leia a primeira coluna. A 5° a aproximação sin θ ≈ θ erra por cerca de uma parte em mil, e mesmo a 10° é meio por cento. O pórtico em serviço gira no máximo alguns graus sob carga de serviço, então θ genuinamente é sin θ até uma fração de por cento, e o cosseno, mantido a mais um termo, é ainda melhor. Esse erro minúsculo é exatamente o espaço em que vive a análise estrutural de primeira ordem.

Agora o eixo em torno do qual todo o artigo gira. Quando você escreve as equações de equilíbrio de uma estrutura, soma de momentos e soma de forças, sobre sua geometria original, indeformada, você faz duas suposições silenciosas sem nomeá-las. Ambas são truncagens de pequenos ângulos de uma série de Taylor.

Primeira, você toma cos θ ≈ 1. Fisicamente isso diz que a carga axial permanece vertical e seu braço de alavanca até o nó não muda conforme o elemento se inclina, porque a geometria sobre a qual você escreve o equilíbrio não se moveu. Segunda, você toma sin θ ≈ θ onde quer que uma rotação entre num balanço de forças. Juntas elas mantêm apenas os termos de Taylor de ordem zero e um da trigonometria e descartam todo o resto. A estrutura nas equações está congelada em sua pose descarregada.

Esse é o significado preciso da análise de primeira ordem: equilíbrio na forma indeformada, que é equilíbrio com a série de pequenos ângulos truncada depois do primeiro termo. A flecha que você calcula em seguida é consequência das cargas, mas nunca é realimentada de volta na geometria sobre a qual o equilíbrio foi escrito. A carga e a forma ficam desacopladas pela truncagem.

O importante é que a análise de primeira ordem não é errada. É uma truncagem deliberada cujo erro você pode limitar com o mesmo termo de resto que limita qualquer série de Taylor, e para as pequenas rotações de um pórtico rígido e bem contraventado esse erro é desprezível. O problema só começa quando os termos descartados crescem o bastante para importar, que é exatamente o caso para o qual a análise de segunda ordem foi construída. Então a forma honesta de descrever a análise de primeira ordem não é "a aproximada", mas "a que mantém o primeiro termo de Taylor".

A mesma truncagem se esconde um nível mais fundo, dentro da própria equação da viga. A curvatura exata de um elemento fletido não é a segunda derivada y″. É

κ = y″ / (1 + y′²)^3/2

onde y′ é a inclinação da forma defletida. Esse denominador é o que uma série de Taylor pode achatar. Expanda (1 + y′²)^−3/2 como uma série binomial, que é só a série de Taylor de (1 + u)^−3/2 em u = y′²:

(1 + y′²)^−3/2 = 1 − (3/2) y′² + (15/8) y′⁴ − ...

Para um elemento cuja inclinação y′ é pequena, isto é, cuja rotação é um ângulo pequeno, todo termo além do primeiro é minúsculo, e o fator inteiro é ≈ 1. A curvatura exata colapsa em

κ ≈ y″

que é a relação linear de Euler-Bernoulli M = EI y″ que toda fórmula de flecha, toda matriz de rigidez e toda verificação de norma silenciosamente supõe. É a mesma truncagem de pequenos ângulos que a do equilíbrio na seção anterior, da mesma família de séries de Taylor, aplicada a uma quantidade diferente. Uma vez que você vê que M = EI y″ é uma série truncada, o resto da escada é fácil de subir: mantenha mais termos e você obtém uma estrutura mais honesta.

Até aqui toda truncagem jogou fora os termos de Taylor mais altos. A análise de segunda ordem faz o oposto: mantém o próximo. O truque é simplesmente deixar a geometria se mover antes de escrever o equilíbrio.

Imagine nosso pilar em balanço. Uma carga lateral empurra o topo livre para o lado por uma quantidade Δ. Na análise de primeira ordem a carga axial P desliza por um elemento que, no que as equações sabem, ainda está perfeitamente vertical, então P não contribui com momento algum na base. Mas o pilar real se inclinou. A carga axial P agora fica pendurada sobre a base pela flecha Δ, e acrescenta um momento

P × Δ

na base engastada, além do que a carga lateral já produziu. Esse momento extra é precisamente o termo que a análise de primeira ordem descartou quando congelou a geometria. É o termo seguinte na expansão do equilíbrio em torno da forma deformada: em vez de parar na pose indeformada, você leva o balanço um passo de Taylor adiante e deixa a carga agir através do deslocamento que ela causou.

O nome P-Delta vem direto desse produto. A carga P vezes o deslocamento lateral Δ é o momento de segunda ordem, e qualquer análise que o inclua é uma análise de segunda ordem. Intuitivamente, o balanço inclinado é a imagem mais limpa possível: quanto mais se inclina, mais a carga axial ultrapassa a base, mais momento faz, o que o inclina ainda mais. Esse laço de realimentação é a história toda, e a próxima seção o transforma numa série que você já conhece.

Siga o laço de realimentação uma passagem por vez e o amplificador se escreve sozinho. Aplique a carga lateral H e o pilar deflete a algum valor de primeira ordem Δ₀. Agora a carga axial P, pendurada sobre a base por Δ₀, acrescenta um momento P vezes Δ₀ que flete o pilar mais ainda, por um extra Δ₁ = Δ₀ × ρ, onde ρ = P/Pcr é a fração do caminho até a carga crítica. Mas essa nova flecha deixa P ultrapassar ainda mais, acrescentando Δ₂ = Δ₀ × ρ², depois Δ₃ = Δ₀ × ρ³, e assim por diante para sempre.

A flecha total é a soma de cada passagem:

Δ = Δ₀ (1 + ρ + ρ² + ρ³ + ...)

Esse parêntese é uma série geométrica, a série de Maclaurin de 1/(1 menos x) que você conheceu em Cálculo II, avaliada em x = ρ. Para ρ abaixo de 1 ela soma a uma forma fechada limpa, 1/(1 menos ρ), então

Δ = Δ₀ / (1 − ρ)

e o momento amplifica pelo mesmíssimo fator. O amplificador é

B = 1 / (1 − P/Pcr)

e agora você pode lê-lo pelo que ele é: cada passagem da carga axial pela flecha crescente é mais um termo da série de Taylor. A análise de primeira ordem mantém apenas o 1 inicial. A análise de segunda ordem, em sua forma de primeira ordem amplificada, soma a série inteira. O Pcr no denominador é a carga crítica de Euler do pilar: a carga em que ρ alcança 1, a série para de convergir, e o amplificador explode. Essa divergência não é um defeito na álgebra, é flambagem, contada na linguagem de uma série de Taylor que ficou sem espaço.

Aqui está a estrutura que carregaremos pelo resto do artigo, com números reais e um motor real por trás de cada um deles. É um pilar de aço em balanço, perfil HEB 200, altura L = 4.0 m, engastado na base e livre no topo. Flexão e flambagem ambas governam em torno do eixo fraco, cujo momento de inércia é I_y = 2001.54 cm⁴, e o aço tem E = 20000 kN/cm².

Passo 1, a resposta de primeira ordem. Empurre o topo livre para o lado com uma carga lateral H = 10 kN. O momento na base é simplesmente H vezes a altura:

M₀ = H · L = 10 × 4.0 = 40.0 kNm.

O motor FEM do CalcSteel, o mesmo que roda dentro do app do navegador, retorna 40.0000 kNm, uma diferença de 0.000%. A flecha de primeira ordem no topo de um balanço é a cúbica clássica:

Δ₁ = H · L³ / (3 E I) = 53.29 mm,

e o motor retorna 53.29 mm em cada dígito, porque o elemento de viga cúbico é exato para este carregamento.

Passo 2, a carga crítica de Euler. A análise de flambagem linear do motor dá Pcr = 617.32 kN. A forma fechada para um pilar engastado-livre é Pcr = π² E I / (2L)², com fator de comprimento efetivo K = 2 porque o balanço livre flamba numa quarta parte de onda senoidal. Essa fórmula dá 617.32 kN. Motor e mão concordam a 0.000%.

Passo 3, adicione carga axial e amplifique. Agora pressione uma carga axial de serviço P = 247 kN por este mesmo pilar. A razão de cargas é

ρ = P / Pcr = 247 / 617.32 = 0.400,

então o amplificador é B = 1 / (1 − 0.400) = 1.667. Aplicando-o aos resultados de primeira ordem:

M₂ = B · M₀ = 1.667 × 40.0 = 66.68 kNm, e Δ₂ = B · Δ₁ = 88.84 mm.

Leia isso sem rodeios: a apenas 40% do caminho até Pcr, os efeitos de segunda ordem acrescentaram 67% ao momento de primeira ordem. Os 26.68 kNm extras (isto é 66.68 menos 40.0) são nada além da soma de todo termo de Taylor além do primeiro. Você pode ver as somas parciais geométricas subirem rumo ao amplificador, um termo por passagem da carga:

A análise de primeira ordem mantém apenas o 1.0000 inicial. Toda linha extra é mais um termo da mesma série convergindo para 1.667.

Uma palavra de honestidade sobre o que o motor fez e o que não fez. O motor FEM dá a resposta exata de primeira ordem e o Pcr exato de Euler, batidos até o decimal acima. O amplificador é a estimativa clássica de segunda ordem em forma fechada, uma análise elástica de primeira ordem amplificada (o método B1/B2 do AISC), não uma resolução separada. É a série geométrica em ação, e é precisa exatamente enquanto a série converge, que é o assunto de uma seção posterior.

Você acabou de ver uma série geométrica virar um momento real no papel. Agora faça você mesmo, sem pegar um lápis. A calculadora embutida aqui mesmo roda uma análise de flambagem num pilar que você define, então o Pcr abstrato vira um número que você pode manipular.

Coloque sua própria altura, perfil e condições de extremidade, leia a carga crítica Pcr direto, depois escolha uma carga axial de serviço P e forme a razão ρ = P/Pcr e o amplificador B = 1/(1 − ρ) você mesmo. Mude a altura, veja Pcr cair e o amplificador subir. É o mesmo laço do pilar resolvido: a série geométrica 1 + ρ + ρ² + ... tornada visível.

É genuinamente grátis, não há login para a matemática, e não há limite de quantas vezes você a roda. Explore-a como exploraria uma série em Cálculo II: mude uma entrada, segure o resto fixo, e leia quão rápido o amplificador cresce conforme você se aproxima de Pcr.

Toda truncagem tem uma etiqueta de aviso, e para a série de Taylor é o termo de resto. Olhe de volta a tabela de pequenos ângulos: o erro de sin θ ≈ θ era 0.127% a 5° mas havia crescido para 4.720% a 30°. O resto cresce com o ângulo, e cresce também a pergunta honesta de se a truncagem ainda vale. Há dois modos distintos de ela falhar.

Falha um, as próprias rotações ficam grandes. Cabos esbeltos, elementos muito flexíveis e problemas de busca de forma podem girar o bastante para que sin θ não seja mais θ: o primeiro termo de Taylor largado é agora uma fração real da resposta. Aqui nenhum amplificador salva você, porque a geometria linearizada sobre a qual o amplificador é construído já expirou. Você precisa de uma análise verdadeira de grande deslocamento, geometricamente não linear que mantenha sin θ inteiro em vez de truncado, e atualize a geometria conforme a estrutura se deforma. Este é o regime em que a imagem de Taylor deixa de ser um atalho e o motor tem que carregar a trigonometria completa.

Falha dois, a própria série do amplificador para de convergir. Conforme a carga axial P sobe rumo a Pcr, a razão ρ = P/Pcr se aproxima de 1, e a série geométrica 1 + ρ + ρ² + ... diverge: B = 1/(1 − ρ) dispara para o infinito. Muito antes de literalmente divergir, o atalho de primeira ordem amplificada perde precisão, e é por isso que as normas o cercam. O AISC limita o uso do amplificador B1/B2, e o Eurocode EN 1993-1-1 pede uma análise completa de segunda ordem sempre que a razão crítica elástica α_cr = Pcr/P cai abaixo de 10. Esse α_cr é só 1/ρ: α_cr ≥ 10 significa ρ ≤ 0.1, a carga está no máximo a um décimo do caminho até a flambagem, a série converge rápido, e a imagem de pequenos ângulos está confortável. Quando α_cr < 10, a carga está perto o bastante de Pcr para que os termos mais altos importem e uma resolução rigorosa de segunda ordem seja exigida. Nosso pilar resolvido ficou em ρ = 0.400, então α_cr = 2.5, bem dentro da região onde uma análise rigorosa de segunda ordem, não um simples amplificador, é a escolha honesta.

O termo de resto é o que torna ambos os limites precisos em vez de uma questão de gosto. É a etiqueta de aviso honesta impressa na aproximação de pequenos ângulos: segura enquanto os termos largados são minúsculos, perigosa no instante em que não são.

Dê um passo atrás e o arco inteiro de Cálculo I e II pousa em aço real, uma ferramenta por vez. A derivada ordinária dá a taxa de variação ao longo de uma viga, dM/dx é igual à força cortante V. A segunda derivada dá a curvatura da linha elástica, M = EI y″. A integral definida recupera áreas e propriedades de seção. A derivada parcial dá uma flecha direto da energia armazenada, como mostra o teorema de Castigliano. E agora a série de Taylor dá o amplificador de segunda ordem, e revela que a análise de primeira ordem por baixo de todas as outras era uma truncagem dessa mesma série.

É um único cálculo coerente, não uma prateleira de truques desconexos. A série de Taylor é a ferramenta que finalmente nomeia a aproximação que as demais faziam sem alarde, e lhe dá o termo seguinte quando você precisa. Se você quer ver de onde vem a carga crítica Pcr no denominador, nosso post sobre flambagem de Euler a deriva de ponta a ponta. Mesmo aço, mais um membro da família do cálculo.

O amplificador não é uma curiosidade de sala de aula. Está escrito, quase letra por letra, nas normas de projeto que um engenheiro assina. Cada uma delas pega a mesma série geométrica e a transforma numa verificação exigida.

• AISC 360 permite uma análise elástica de primeira ordem amplificada com dois fatores: B1 para o efeito P-δ do elemento (uma carga axial agindo através da flecha dentro de um elemento) e B2 para o efeito P-Δ do pavimento (a carga agindo através do deslocamento lateral do pavimento inteiro). Ambos são da forma 1/(1 − α P_r/P_e), a mesma série 1/(1 − ρ) com P_e a carga de flambagem elástica. O Capítulo C expõe o método de análise direta e o Apêndice 8 a rota B1/B2.
• Eurocode EN 1993-1-1 atrela toda a decisão a α_cr: se α_cr ≥ 10 (elástico) os efeitos de segunda ordem são pequenos o bastante para ignorar ou tratar por amplificação, e abaixo disso uma análise de segunda ordem mais completa é exigida. Esse é o teste de convergência da série, promovido a cláusula.
• ABNT NBR 8800 carrega seus próprios coeficientes B1 e B2 no mesmo espírito de primeira ordem amplificada, então um calculista brasileiro está somando a série de Taylor idêntica sob notação brasileira.

Dois vizinhos na mesma família de normas valem um link. O deslocamento lateral que B2 amplifica precisa ser disparado antes de mais nada, que é o trabalho das cargas nocionais, pequenas forças laterais que expõem a instabilidade de um pórtico imperfeito sem uma rodada completa de P-Delta. E o momento amplificado não age sozinho: ele se combina com a carga axial na verificação de interação de flexão composta com força axial. Em todo caso o ponto é o mesmo: quantos termos de Taylor você mantém é uma decisão de norma, não uma vibração. A norma diz quando o primeiro termo basta e quando você precisa somar o resto.

A ideia é curta, que é exatamente por que os erros são fáceis de cometer. Aqui está a lista de checagem que mantém uma estimativa de segunda ordem honesta.

• "Pequeno ângulo" não é "pequena flecha". A truncagem que a análise de primeira ordem faz é sobre a rotação, não sobre o deslocamento. Um balanço esbelto pode mover o topo por uma distância grande (nosso pilar resolvido alcançou 88.84 mm) enquanto seus elementos ainda giram apenas alguns graus, então sin θ ≈ θ segue válido mesmo que a flecha seja visivelmente grande. Confundir as duas leva as pessoas a desconfiar de uma análise perfeitamente boa, ou a confiar numa ruim.
• Usar ρ = P/Pcr com o Pcr errado. Este é o erro real de sempre. A carga crítica depende do fator de comprimento efetivo K e de qual eixo flamba; ponha o Pcr do eixo forte, ou K = 1 para um pilar que é de fato engastado-livre (K = 2), e seu amplificador não tem sentido. Note que nosso pilar resolvido flambou em torno do eixo fraco, com I_y = 2001.54 cm⁴, que é o eixo que governa. Acerte Pcr antes de formar ρ.
• Segunda ordem não é o mesmo que não linearidade material. A análise inteira neste artigo é linear elástica: o aço nunca escoa, a lei de Hooke vale por todo lado. A não linearidade de segunda ordem (geométrica) é sobre equilíbrio na forma deformada, coisa diferente de uma curva tensão-deformação não linear. Um problema pode ser uma, a outra, ambas ou nenhuma.
• B explode perto de Pcr por um motivo real. Quando alguém vê 1/(1 − ρ) disparar para o infinito às vezes supõe uma falha numérica. É física: em ρ = 1 a carga alcançou a flambagem de Euler, e a amplificação infinita é a face matemática honesta de um pilar instável. Respeite a divergência, não a corte.
• Primeira ordem é uma truncagem, não um erro. Enquanto α_cr = Pcr/P estiver confortável (10 ou mais em termos do Eurocode), manter apenas o primeiro termo de Taylor é uma escolha legítima, sancionada por norma. A análise de primeira ordem merece seu lugar precisamente porque os termos largados são comprovadamente minúsculos ali.

Um pilar inclinado, uma série de Taylor, lida em três profundidades de truncagem, e todo número de primeira ordem e de Euler batido com um motor FEM real até o decimal.

• A análise de primeira ordem é uma truncagem de pequenos ângulos. O equilíbrio na forma indeformada mantém sin θ como θ e cos θ como 1, apenas o primeiro termo de Taylor, e a geometria nunca se move.
• O amplificador é uma série geométrica. B = 1/(1 − P/Pcr) é exatamente 1 + ρ + ρ² + ρ³ + ..., um termo por passagem da carga axial pela flecha crescente.
• O pilar HEB 200 resolvido. A primeira ordem deu M₀ = 40.0 kNm e Δ₁ = 53.29 mm; em ρ = 0.400 o amplificador B = 1.667 os elevou a M₂ = 66.68 kNm e Δ₂ = 88.84 mm, 67% mais momento a apenas 40% do caminho até Pcr.
• O motor bateu com a conta à mão até o decimal. Sua análise de flambagem retornou Pcr = 617.32 kN, igual à forma fechada de Euler π²EI/(2L)² com K = 2 a 0.000%.
• O resto é a etiqueta de aviso. Perto de Pcr a série diverge, e em grandes rotações sin θ deixa de ser igual a θ. Ali a truncagem falha e você precisa de análise completa de segunda ordem ou de grandes deslocamentos.

Agora ponha as mãos nisso. Construa o pilar, leia Pcr, e forme o amplificador você mesmo no CalcSteel, uma ferramenta nativa do navegador para projeto e análise de estruturas de aço com um motor FEM genuíno. A matemática e as calculadoras são grátis, sem login, ilimitadas, e se você é estudante a oferta /education é grátis para você. A prova não é um teste com cronômetro; é uma ferramenta que você pode abrir agora mesmo e conferir contra cada número desta página.