Estruturas Estaticamente Indeterminadas: Por Que o Método Manual Para e o Método Matricial Começa

Pegue uma viga biapoiada, um pino, um rolete e uma carga uniforme. Você consegue achar as duas reações com nada além de uma folha de papel e as três equações da estática: somatório das forças horizontais, somatório das forças verticais, somatório dos momentos, cada um igual a zero. Todo esforço cortante e momento interno decorre dessas reações. A viga é estaticamente determinada (isostática): o número de incógnitas é igual ao número de equações de equilíbrio, e só o equilíbrio já lhe diz tudo.

Agora solde a extremidade esquerda dessa mesma viga a uma parede, ou coloque um terceiro apoio sob uma viga longa, ou dobre duas colunas para cima para carregá-la como um pórtico. Você adicionou uma reação, mas não adicionou uma equação, a estática continua lhe dando apenas três no plano. A estrutura agora tem mais incógnitas do que o equilíbrio consegue resolver. Ela é estaticamente indeterminada, também chamada de hiperestática, e as forças extras não podem ser obtidas da estática de jeito nenhum. Isso não é um defeito; é a condição normal de quase toda edificação real. Este artigo é sobre o que muda quando você cruza essa linha: o que é a indeterminação, como contá-la, como os engenheiros a resolvem à mão há mais de um século, precisamente onde esses métodos manuais deixam de ser práticos, e como o método matricial que começa ali virou o motor dentro de todo solver estrutural.

Toda a distinção repousa sobre uma contagem. Em duas dimensões um corpo rígido tem exatamente três equações de equilíbrio independentes, então ele consegue resolver três grandezas de força incógnitas e nada mais. Alinhe suas incógnitas, as reações de apoio mais quaisquer redundâncias de esforço interno, contra essas equações:

• Incógnitas em número menor que as equações: a estrutura é um mecanismo, ela se move, é instável. Não é estrutura coisa nenhuma.
• Incógnitas iguais às equações: estaticamente determinada (isostática). O equilíbrio a resolve exatamente, e os esforços internos não dependem da rigidez das barras, só da geometria e da carga.
• Incógnitas em número maior que as equações: estaticamente indeterminada. O equilíbrio é necessário mas já não é suficiente, e as equações que faltam têm de vir de outro lugar.

Esse outro lugar é a compatibilidade: a exigência de que as flechas e rotações da estrutura encaixem umas nas outras, de que uma viga contínua permaneça contínua sobre o apoio central, de que um engaste realmente não gire. As condições de compatibilidade são equações extras, e são escritas em termos de como a estrutura se deforma. No instante em que uma estrutura é indeterminada, seus esforços internos passam a depender da rigidez relativa de suas barras, de E, de I, da seção e do comprimento, porque a rigidez é o que decide como a força redundante é repartida. Esse único fato, o de que forças em estruturas indeterminadas dependem da rigidez enquanto as de estruturas determinadas não, é a raiz de tudo que se segue.

O grau de hiperestaticidade é uma subtração simples, incógnitas menos equações disponíveis. Para uma viga plana ou uma estrutura simples cujas únicas incógnitas são reações:

GHE = r − 3 − c

onde r é o número de componentes de reação, 3 são as equações de equilíbrio no plano, e c é o número de equações de condição liberadas por rótulas internas (uma rótula interna adiciona uma equação, a de que o momento ali é zero). Um apoio de pino dá dois componentes de reação, um rolete dá um, um engaste dá três.

Para um pórtico plano rígido, conte também as barras e os nós, porque cada barra ligada rigidamente carrega suas próprias três grandezas de esforço interno:

GHE = (3m + r) − 3n − c

com m barras, r componentes de reação, n nós (incluindo os apoios), e c liberações. Aplique aos casos da tabela. Uma viga engastada e apoiada é um engaste (3) mais um rolete (1), então r = 4 e GHE = 1. Uma viga contínua de dois vãos sobre um pino e dois roletes também tem r = 4 e GHE = 1. Um pórtico sobre duas bases engastadas tem m = 3, r = 6, n = 4, então GHE = 9 + 6 − 12 = 3. Cada unidade de GHE é uma redundante, uma força que você pode escolher livremente antes que a compatibilidade a fixe.

Se a indeterminação torna a análise mais difícil, por que quase todo pórtico real é indeterminado? Porque a mesma continuidade que complica a matemática compra três coisas que você quer.

Rigidez. A continuidade e o engastamento combatem a flecha com força. Uma viga biengastada sob carga uniforme flete apenas um quinto do que a mesma viga biapoiada (wL⁴/384EI contra 5wL⁴/384EI). Para barras governadas por serviço, que é a maioria das vigas de piso, essa é a diferença entre uma seção que funciona e uma que não funciona.

Momentos de pico menores. O engastamento espalha a flexão. Aquela viga biengastada tem pico de wL²/12 nas extremidades em vez de wL²/8 no meio do vão, um terço a menos de momento para o mesmo vão e carga. A continuidade sobre um apoio faz o mesmo truque, trocando um momento de vão menor por um momento negativo no apoio.

Redundância, ou seja, segurança. O sentido literal de uma redundante é um caminho de carga alternativo. Se uma seção de uma estrutura indeterminada é sobrecarregada e escoa, a estrutura não necessariamente cai; o momento se redistribui para as partes que ainda têm capacidade, e ela continua de pé enquanto avisa você. Uma estrutura determinada não tem caminho reserva: perca uma barra ou forme uma rótula e ela vira um mecanismo. É por isso que as normas premiam a redundância e por que o raciocínio de colapso progressivo se apoia nela. A indeterminação não é um custo que você tolera; é uma propriedade que você projeta de propósito.

A maneira mais antiga e mais intuitiva de resolver uma estrutura indeterminada é o método das forças, também chamado de método da flexibilidade ou método das deformações consistentes. A ideia é desarmadoramente simples. Você tem forças incógnitas demais, então escolhe GHE delas, as redundantes, e imagina removê-las. O que sobra, a estrutura primária, é determinada, então a estática a resolve de imediato.

Claro, remover um apoio deixa a estrutura se deformar de um jeito que a real nunca faria: um apoio liberado deixa a viga afundar onde o apoio deveria tê-la segurado. Então você calcula essa flecha ilegal na liberação sob a carga aplicada, chame-a de δ₁₀. Depois você recoloca a redundante como uma força incógnita X₁ e calcula a flecha que ela produz por unidade, δ₁₁. A estrutura real não flete no apoio de forma alguma, então as duas têm de se cancelar. Essa é a equação de compatibilidade:

δ₁₀ + X₁·δ₁₁ = 0

Uma redundante, uma equação, uma força incógnita recuperada. Com mais redundantes vira um sistema pequeno, [δ]{X} = −{δ₀}, onde [δ] é a matriz de flexibilidade, uma linha e coluna por redundante. Resolva-o, recoloque as redundantes, e termine com a estática. O método das forças é o retrato mais limpo possível do que a indeterminação é: exatamente GHE condições, cada uma um enunciado geométrico de que a estrutura encaixa.

O método das forças trata as forças como incógnitas. A família rival, os métodos dos deslocamentos (da rigidez), trata os movimentos dos nós como incógnitas, e para pórticos reais em geral ganha, porque um pórtico costuma ter menos rotações nodais livres do que forças redundantes.

O método dos deslocamentos (1915) escreve os momentos de extremidade de cada barra como uma fórmula nas rotações de suas duas extremidades e em qualquer deslocamento lateral entre elas, e depois impõe o equilíbrio de momentos em cada nó. Isso dá uma equação por rotação incógnita, um sistema que você resolve simultaneamente. A distribuição de momentos, o método de Cross de 1930, foi o avanço que colocou a análise de pórticos nas mãos de qualquer um com uma régua de cálculo: ela trava todos os nós, calcula os momentos de engastamento, depois libera os nós um de cada vez e deixa o momento desequilibrado se distribuir para as barras conectadas conforme sua rigidez relativa, iterando até que as sobras sejam desprezíveis. Nenhum sistema de equações simultâneas, apenas escrituração disciplinada, e converge rápido para vigas contínuas e pórticos contraventados. A equação dos três momentos de Clapeyron (1857) é a ferramenta especialista para vigas contínuas, relacionando os momentos sobre três apoios consecutivos em uma única equação por apoio interno.

Esses métodos são genuinamente poderosos, e todo engenheiro estrutural deveria saber rodar a distribuição de momentos em uma viga de dois vãos ou em um pórtico de um único vão. Mas repare no que cada um custa à medida que a estrutura cresce, porque esse custo é a razão inteira de o método matricial existir.

Pegue a primeira estrutura indeterminada canônica: uma viga engastada e apoiada (propped cantilever), engastada numa extremidade e apoiada sobre um rolete na outra, vencendo L = 6 m sob uma carga uniforme w = 24 kN/m. É GHE 1, então o método das forças precisa de exatamente uma equação de compatibilidade.

Escolha a reação do apoio como a redundante X₁ e libere-a. A estrutura primária é um console simples engastado à esquerda. Sob a carga w, a extremidade livre (apoiada) de um console afunda δ₁₀ = wL⁴/8EI. Sob uma força unitária para cima nessa extremidade ela sobe δ₁₁ = L³/3EI. O apoio real segura a extremidade em zero, então:

X₁ = −δ₁₀ / δ₁₁ = (wL⁴/8EI) / (L³/3EI) = 3wL/8 = 54 kN

Repare que EI se cancela: para uma única barra prismática a redundante afinal não depende da rigidez, então a resposta é exata e limpa. Agora termine com a estática. A reação do engaste é wL − X₁ = 144 − 54 = 90 kN (isto é 5wL/8), e o momento de engastamento é wL²/2 − X₁·L = 432 − 324 = 108 kN·m de momento negativo, exatamente wL²/8. O momento no vão tem pico em x = 5L/8 = 3,75 m com valor de 9wL²/128 = 60,75 kN·m, e o momento fletor passa por zero (o ponto de inflexão) em x = L/4 = 1,5 m.

Agora a verificação. Montamos a mesma viga no motor de EF do CalcSteel, o solver real que está em produção, e ele devolve R no engaste = 90,000 kN, R no apoio = 54,000 kN, o momento de engastamento = −108,000 kN·m, e o pico no vão = +60,750 kN·m em x = 3,750 m. Cada número bate com a solução manual até a terceira casa decimal, que é o ponto: para uma redundante o método manual e o motor dão a mesma resposta, e o motor apenas fez a compatibilidade por você.

Acrescente um apoio em vez de um engaste e você obtém o outro caso cotidiano de GHE 1, uma viga contínua de dois vãos: um pino, um rolete interno e um rolete de extremidade, dois vãos de L = 6 m cada, o mesmo w = 24 kN/m ao longo de tudo. Aqui a ferramenta manual natural é a equação dos três momentos de Clapeyron, que para dois vãos iguais com carga uniforme dá o momento no apoio interno diretamente:

M_B = −wL²/8 = −108 kN·m

O apoio interno faz exatamente o que a continuidade compra: puxa um momento negativo para fora da viga e, em troca, mantém os momentos de meio de vão baixos, em 9wL²/128 = 60,75 kN·m, exatamente o mesmo momento de campo da viga engastada e apoiada, o que não é coincidência, já que cada vão desta viga se comporta como uma viga engastada e apoiada, com o apoio interno fazendo o papel do quase-engaste. A estática então dá as reações: 3wL/8 = 54 kN em cada extremidade e 10wL/8 = 1,25wL = 180 kN no apoio interno, que carrega a maior parte.

O motor concorda em toda a linha: reações de extremidade 54,000 kN, reação interna 180,000 kN, momento no apoio interno −108,000 kN·m, picos de vão +60,750 kN·m. Duas estruturas indeterminadas diferentes, uma resolvida pelo método das forças e outra pela equação dos três momentos, e em ambas o solver de EF em produção reproduz o resultado clássico exatamente. Essa é a tranquilidade que os métodos manuais lhe dão: em GHE 1 você ainda consegue manter a solução inteira na cabeça, e o software é auditável contra ela.

Em GHE 1 o método manual é um prazer. O problema é que o esforço não cresce gentilmente. No método das forças a matriz de flexibilidade é GHE por GHE, então uma estrutura com três redundantes precisa de um sistema 3 por 3, seis precisam de 6 por 6, e cada um daqueles coeficientes δ_ij é ele mesmo um cálculo de flecha separado que você tem de resolver e acertar. Nos métodos dos deslocamentos a contagem de incógnitas é o número de deslocamentos nodais livres, que sobe igualmente rápido num edifício real, e o método dos deslocamentos os transforma num sistema simultâneo enquanto a iteração da distribuição de momentos desacelera e se multiplica assim que o deslocamento lateral (sidesway) entra, porque um pórtico livre para pender de lado precisa de ciclos extras de deslocamento lateral empilhados sobre o balanceamento dos nós.

Então o muro é real, e ele chega cedo. Um pórtico de um único vão (GHE 3) já é um cálculo manual sério, uma tarde de método dos deslocamentos com uma equação de deslocamento lateral, ou distribuição de momentos feita duas vezes. Um pórtico de dois pavimentos e dois vãos está além do que qualquer pessoa sã faz à mão. Um edifício inteiro, uma grelha, qualquer coisa em três dimensões, simplesmente não é um problema para fazer à mão, e não há método esperto esperando para resgatá-lo, porque a dificuldade não é a técnica, é o número absoluto de incógnitas acopladas. Os métodos manuais não falharam; eles esbarraram no tamanho do sistema linear que um humano consegue carregar.

O método matricial não supera os métodos manuais em esperteza; ele os supera em organização. É o método dos deslocamentos, a mesma ideia de movimento-do-nó-como-incógnita por trás do método dos deslocamentos, escrito numa forma que um computador consegue montar mecanicamente e resolver em qualquer escala. Ele é chamado de método matricial da rigidez direta, e é a espinha dorsal da análise por elementos finitos.

A receita nunca muda com o tamanho. Cada barra ganha uma pequena matriz de rigidez do elemento k que relaciona as forças em suas duas extremidades aos deslocamentos em suas duas extremidades, pura geometria e EA, EI, L. Cada uma dessas é rotacionada para os eixos globais e somada a uma matriz de rigidez global K nas linhas e colunas dos nós que ela conecta; barras que se sobrepõem simplesmente somam onde compartilham um nó. Aplique as cargas como um vetor F e as condições de contorno, e a estrutura inteira agora é um único sistema linear:

K d = F

Resolva-o uma vez para os deslocamentos nodais d, depois substitua de volta através da k de cada elemento para recuperar os esforços cortantes e momentos de cada barra. Não há tratamento separado para determinada ou indeterminada, nenhuma escolha de redundantes, nenhuma equação de deslocamento lateral, nenhuma iteração; a indeterminação é tratada implicitamente no momento em que as rigidezes são montadas, porque K já codifica como cada barra reparte carga com suas vizinhas. Dobrar o tamanho do pórtico só torna K maior, e um computador não se importa. É exatamente isto que o motor do CalcSteel roda, o mesmo método, na precisão da máquina; os três casos resolvidos acima foram resolvidos por ele, e o pórtico a seguir também.

A calculadora abaixo é o método da rigidez direta, rodando no seu navegador sobre o mesmo motor de EF dos casos resolvidos. Monte uma viga engastada e apoiada ou uma viga contínua, engaste uma extremidade ou adicione um terceiro apoio, coloque uma carga uniforme nela, e leia as reações e o diagrama de momento fletor direto de K d = F. Teste nos casos GHE 1 acima e veja os números caírem em 90, 54, −108 e 60,75 sem você jamais escolher uma redundante.

É gratuita, sem login para a análise, e nunca para em GHE 1, adicione apoios e barras até ter uma estrutura que nenhum método manual conseguiria tocar, e ela resolve no mesmo instante.

Termine com o caso que não tem fórmula em manual nenhum: um pórtico, uma viga de 8 m sobre duas colunas de 5 m engastadas em suas bases, carregando uma carga de gravidade de 20 kN/m na viga e um empurrão horizontal de vento de 15 kN no topo. É GHE 3, e o vento o faz pender, então as três redundantes se acoplam através do deslocamento lateral. Não há forma fechada; este é um problema só matricial.

O motor devolve a solução inteira de uma vez. As duas bases das colunas, embora o pórtico e sua carga de gravidade sejam perfeitamente simétricos, carregam momentos de engastamento bem diferentes, 4,25 kN·m na base a barlavento e 34,45 kN·m na a sotavento. Os dois joelhos de ligação viga-coluna são igualmente desiguais, 12,07 kN·m a barlavento contra 48,37 kN·m a sotavento, e a viga afunda até um pico de 129,8 kN·m no meio do vão. As reações verticais se dividem em 75,46 kN e 84,54 kN à medida que o pórtico pende a sotavento, e os 15 kN de vento são repartidos em 1,56 kN e 16,56 kN entre as bases. Cada uma dessas assimetrias é o deslocamento lateral falando, e nenhuma equação manual isolada as produz; você precisaria do método dos deslocamentos com uma incógnita de deslocamento lateral, com cuidado, para chegar sequer perto.

As verificações de sanidade ainda fecham à mão, que é como você confia no solver: as reações verticais somam 75,46 + 84,54 = 160 kN, exatamente os 20 kN/m sobre 8 m de gravidade, e as reações horizontais somam 1,56 − 16,56 = −15 kN, equilibrando exatamente o vento. O equilíbrio é sempre necessário; aqui o método matricial forneceu a compatibilidade que só o equilíbrio jamais alcançaria.

Confundir indeterminada com instável. Ambas têm um descompasso entre incógnitas e equações, mas em direções opostas. Poucas incógnitas (ou mal arranjadas) é um mecanismo, ele se move. Muitas é indeterminada, ela é mais rígida do que a estática consegue descrever. Sempre verifique a estabilidade e o arranjo, não só a contagem; três reações todas paralelas ou todas passando por um único ponto é instável mesmo que o número pareça certo.

Achar que as forças em estruturas indeterminadas independem da seção. Não independem. No momento em que uma estrutura é indeterminada, os esforços internos dependem da rigidez relativa, de E e I e do comprimento da barra. Mude uma seção e os momentos se redistribuem. Forças em estruturas determinadas não ligam para a rigidez; as em indeterminadas sempre ligam.

Esquecer as equações de condição das rótulas internas. Uma rótula interna libera um momento e adiciona uma equação, baixando o GHE em um. Uma viga que parece superapoiada pode ser determinada, ou até um mecanismo, uma vez contadas suas rótulas. Não leia o GHE só dos apoios.

Uma estrutura estaticamente indeterminada é mais forte? Em geral mais rígida e mais redundante, não automaticamente mais forte no sentido último, mas a redundância dá caminhos de carga alternativos e deixa os momentos se redistribuírem antes do colapso, o que é uma margem de segurança real.

Ainda preciso dos métodos manuais se o software resolve tudo? Sim. Você não consegue auditar uma resposta que não consegue estimar. A distribuição de momentos numa viga de dois vãos ou um wL²/8 de cabeça mantêm o solver honesto e pegam o erro de modelagem, o apoio errado, a liberação faltante, que transforma um método certo num número errado.

Qual é a diferença entre o método da rigidez e os elementos finitos? Nenhuma, em espírito. O método da rigidez direta para pórticos é análise por elementos finitos com um elemento de viga por barra; o MEF é a mesma montagem de K d = F generalizada para placas, cascas e sólidos.

A indeterminação é o estado normal das estruturas reais, e cruzar para dentro dela muda tanto a matemática quanto a ferramenta a que você recorre.

• Uma estrutura é estaticamente indeterminada quando as incógnitas superam as equações de equilíbrio; o excedente é o grau de hiperestaticidade, GHE = r − 3 para uma viga plana, (3m + r) − 3n para um pórtico plano.
• Além do equilíbrio você precisa da compatibilidade, e a partir desse ponto os esforços internos dependem da rigidez da barra, de E, I e L, não só da geometria.
• Os métodos manuais adicionam exatamente as condições que faltam: o método das forças resolve uma equação por redundante; o método dos deslocamentos e a distribuição de momentos resolvem para os movimentos dos nós. Todos escalam mal e batem num muro prático em poucos graus de hiperestaticidade.
• Verificado contra o motor até a terceira casa decimal: uma viga engastada e apoiada dá −108 kN·m (wL²/8) no engaste e +60,75 kN·m (9wL²/128) no vão; uma viga contínua de dois vãos gera −108 kN·m sobre o apoio interno; a redundante do método das forças é a reação do apoio, 3wL/8 = 54 kN.
• O método matricial (rigidez direta) escreve a mesma física como K d = F, montado barra por barra e resolvido por computador, sem teto. Ele roda o motor do CalcSteel, e é a única forma prática de resolver um pórtico sob vento (GHE 3) ou qualquer coisa maior.