💻 computer science

A new approach for the determination of through-thickness and free-edge stresses in composite laminates based on structural elements

Este artigo propõe um novo método utilizando elementos coesivos estruturais e uma rigidez de penalidade precisa derivada de camadas ricas em resina para recuperar componentes de tensão fora do plano em elementos de casca fina baseados em Kirchhoff-Love, permitindo, assim, a previsão precisa da iniciação de danos em laminados compostos.

Autores originais: Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

Publicado 2026-07-31

📖 7 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Xiaopeng Ai, Christos Kassapoglou, Boyang Chen

Artigo original sob licença CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Esta é uma explicação gerada por IA do artigo abaixo. Não foi escrita nem endossada pelos autores. Para precisão técnica, consulte o artigo original. Ler aviso legal completo

Imagine que você está construindo um arranha-céu usando milhares de pequenos blocos de Lego coloridos. No mundo real, os engenheiros não usam plástico; eles usam "compósitos laminados", que são como sanduíches superfortes feitos de camadas de tecido embebidas em cola (resina). Esses materiais são o ingrediente secreto por trás de aviões leves e carros rápidos porque são incrivelmente fortes, mas pesam muito pouco. No entanto, existe um problema complicado: quando você dobra ou estica esses sanduíches, as camadas querem deslizar umas sobre as outras ou se separar, especialmente perto das bordas. Isso é chamado de "delaminação", e é um assassino silencioso de estruturas.

Para prever se a asa de um avião resistirá, os engenheiros usam simulações de computador. Para estruturas finas, eles geralmente usam um atalho chamado "elementos de casca" (shell elements). Pense em um elemento de casca como uma folha de papel em um programa de computador. Ele é ótimo para calcular como o papel estica ou dobra em sua superfície, mas é péssimo para entender o que está acontecendo dentro da espessura do papel. É como tentar entender o recheio de um sanduíche olhando apenas para a fatia de pão de cima; você perde o queijo e a carne inteiros. Como esses modelos de computador padrão não conseguem ver as tensões "internas" (as forças que empurram as camadas para longe umas das outras), eles frequentemente perdem os sinais de alerta antes de uma estrutura falhar. Este artigo intervém para corrigir esse ponto cego, permitindo que os engenheiros vejam as forças ocultas dentro do sanduíche sem ter que construir um modelo 3D massivo, lento e complicado de cada grão de cola.

A Grande Ideia do Artigo: Vendo as Camadas Invisíveis

Os autores, Xiaopeng Ai, Christos Kassapoglou e Boyang Chen, criaram uma nova maneira de espiar dentro desses sanduíches compostos usando um truque inteligente. Eles estão construindo sobre um método que trata cada camada do laminado como um "elemento estrutural" (uma versão sofisticada e de alta tecnologia da nossa folha de papel) e as conecta com "elementos coesivos". Você pode pensar nesses elementos coesivos como tiras de cola invisíveis e elásticas entre as camadas.

No passado, essas tiras de cola eram um certo mistério. Os engenheiros tinham que adivinhar o quão rígidas elas deveriam ser. Se eles errassem o palpite, o computador poderia pensar que as camadas estavam coladas com uma epóxi superforte (muito rígida) ou uma fita fraca (muito maleável/elástica), levando a respostas erradas sobre quanta tensão a estrutura poderia suportar. O artigo argumenta que as antigas formas de adivinhar essa rigidez eram falhas porque tratavam a cola como se ela tivesse espessura zero, o que não é verdade no mundo real, onde sempre há um pequeno espaço rico em resina entre as camadas.

A "Receita Secreta" da Camada Rica em Resina

A principal descoberta da equipe é uma nova fórmula para calcular exatamente o quão rígidas essas "tiras de cola" devem ser. Em vez de adivinhar, eles perceberam que a rigidez depende das propriedades reais da camada rica em resina — aquela zona minúscula e carregada de cola entre as camadas de tecido.

Imagine que você está esticando uma pilha de notas adesivas (tipo Post-it). A força que você sente depende de quão grossa é a parte colante e de quão forte é a cola. Os autores derivaram uma regra matemática que diz: "A rigidez da tira de cola é simplesmente a resistência da cola dividida pela espessura dessa minúscula camada de resina".

Eles testaram essa ideia executando simulações em diferentes tipos de sanduíches compostos:

Laminados de cruzamento (cross-ply): Camadas empilhadas como um tabuleiro de xadrez (0° e 90°).
Laminados de ângulo (angle-ply): Camadas empilhadas em ângulos como um zigue-zague (+45° e -45°).
Laminados quase isotrópicos: Pilhas complexas que agem da mesma forma em todas as direções.

Nestas simulações, eles aplicaram um estiramento uniforme nas extremidades do laminado, exatamente como puxar um elástico. Eles então compararam seu novo método contra duas coisas:

Soluções analíticas: Estas são respostas matemáticas perfeitas derivadas por outros cientistas (como Pagano e Mittelstedt) que servem como o "padrão ouro".
Modelos de elementos sólidos: Estes são os modelos de computador pesados e lentos que decompõem o sanduíche em milhões de pequenos blocos 3D para ver cada detalhe.

O Que Eles Descobriram

Os resultados foram impressionantes. Quando os autores usaram sua nova fórmula de rigidez "rica em resina", seus modelos de casca simples previram as tensões internas quase perfeitamente.

Precisão: Seus resultados corresponderam às soluções matemáticas do "padrão ouro" e aos modelos de blocos 3D pesados. Eles puderam prever com precisão a tensão normal interlaminar (a força que tenta separar as camadas, denotada como $\sigma_z$ ) e a tensão de cisalhamento (a força que tenta fazer as camadas deslizarem umas sobre as outras, denotada como $\tau_{xz}$ e $\tau_{yz}$ ).
Efeitos de Borda Livre (Free-Edge Effects): Uma das coisas mais difíceis de prever é o que acontece na borda extrema do material, onde as camadas estão expostas ao ar. É aqui que ocorrem as "singularidades de tensão" (picos extremos de força). O novo método capturou esses picos com precisão, enquanto métodos antigos frequentemente os perdiam ou exigiam malhas incrivelmente finas e lentas para visualizá-los.
Robustez: A equipe também verificou se o método era sensível a pequenas mudanças. Eles perguntaram: "E se a camada de resina for um pouco mais grossa ou mais fina? E se a cola for um pouco mais forte ou mais fraca?". Eles descobriram que seu método era muito estável. Mesmo que a espessura da camada de resina variasse entre 0,6 × 10⁻³ polegadas e 0,9 × 10⁻³ polegadas, ou se a rigidez da cola variasse entre 3 GPa e 5 GPa, os resultados não mudaram muito. Isso significa que os engenheiros não precisam conhecer os detalhes microscópicos exatos da cola para obter uma boa resposta.

O Ganho de Velocidade

Talvez a parte mais emocionante seja a velocidade. Os modelos de blocos 3D pesados (os modelos de "elementos sólidos") são precisos, mas lentos. Para obter uma boa resposta, eles precisavam fatiar o laminado em oito blocos minúsculos para cada camada. O novo método, no entanto, precisava de apenas um elemento de casca por camada.

Em seus testes, o novo método reduziu o tempo de processamento do computador em cerca de 45%. Para um caso de teste específico, o método antigo levou 9,52 segundos por etapa de cálculo, enquanto o novo método levou apenas 5,21 segundos. Isso é uma economia enorme quando você está projetando a asa de um avião inteiro com milhares de partes.

O Que Isso Significa para o Futuro

O artigo não afirma ter resolvido todos os problemas do universo. Eles admitem que, para padrões de tensão muito complexos dentro de uma única camada, o método utiliza uma "interpolação linear" (uma linha reta) para estimar os valores entre as camadas, o que pode perder detalhes minúsculos e ondulados. No entanto, para o panorama geral de como as camadas se separam e deslizam, funciona maravilhosamente bem.

Ao provar que é possível obter dados de tensão 3D precisos usando elementos de casca simples e rápidos, este trabalho abre as portas para que engenheiros simulem danos e delaminação de forma muito mais eficiente. É como atualizar de um telescópio lento e embaçado para um nítido e rápido, permitindo-nos ver os perigos ocultos em nossas estruturas leves antes que elas sequer decolagem. Os autores sugerem que este framework pode ser a base para ferramentas futuras que simulem tanto danos intraplanares (rachaduras no tecido) quanto delaminação (separação de camadas) simultaneamente, tornando nossos aviões e carros mais seguros e leves.

Resumo Técnico: Uma Nova Abordagem para Tensões Transversais e de Borda Livre em Laminados Compósitos

Definição do Problema
Compósitos reforçados com fibras são críticos para estruturas aeroespaciais leves, porém prever sua resistência e tolerância ao dano permanece um desafio. Embora elementos de casca (shell) finos baseados na hipótese de Kirchhoff–Love sejam amplamente utilizados por sua eficiência computacional, eles inerentemente contabilizam apenas três componentes de tensão no plano ( $\sigma_x$ , $\sigma_y$ e $\tau_{xy}$ ). Consequentemente, as componentes de tensão fora do plano necessárias para critérios de dano tridimensionais (3D) convencionais não estão disponíveis. Essa limitação impede a previsão precisa da iniciação de danos, particularmente a delaminação interlaminar e efeitos de borda livre, que são impulsionados por tensões normais e de cisalhamento fora do plano.

Métodos existentes para recuperar essas tensões frequentemente dependem de teorias de casca de Mindlin–Reissner ou técnicas de pós-processamento que ou negligenciam elementos coesivos (limitando a simulação de delaminação) ou enfrentam gargalos de densidade de malha impostos por modelos de zona coesiva. Além disso, formulações anteriores para a rigidez de penalidade em elementos coesivos frequentemente assumem espessura zero ou dependem de suposições cinemáticas inconsistentes com a hipótese de Kirchhoff–Love, levando a níveis de tensão superestimados.

Metodologia
Os autores propõem uma estrutura unificada baseada em elementos estruturais que captura efeitos de borda livre em laminados compósitos multicamadas com empilhamentos arbitrários. A metodologia consiste em três componentes principais:

Formulação do Elemento Estrutural: O laminado é modelado usando elementos de casca cúbicos para cada camada, conectados por elementos coesivos estruturais. Esta abordagem mantém as suposições cinemáticas de Kirchhoff–Love (onde as normais permanecem retas e perpendiculares ao plano médio) enquanto elimina a limitação de densidade de malha da zona coesiva encontrada na modelagem de zona coesiva tradicional.
Mecanismo de Recuperação de Tensão: As componentes de tensão fora do plano ( $\sigma_z$ , $\tau_{xz}$ , $\tau_{yz}$ ) são recuperadas diretamente das trações transmitidas pelos elementos coesivos estruturais. Ao impor o equilíbrio na interface, as trações coesivas são igualadas às tensões fora do plano nas camadas adjacentes. Essas tensões interfaciais são então interpoladas através da espessura para reconstruir o campo de tensão 3D completo.
Rigidez de Penalidade Fisicamente Representativa: Uma inovação crítica é a derivação de uma nova formulação de rigidez de penalidade para os elementos coesivos. Diferente de abordagens anteriores que assumem espessura zero ou deformação de cisalhamento simples, este trabalho define a rigidez de penalidade como uma função das propriedades mecânicas e da espessura da camada rica em resina entre as camadas adjacentes.
- A rigidez de penalidade normal ( $K_n$ ) é derivada como $E_3 / h_{rr}$ , onde $E_3$ é o módulo elástico fora do plano e $h_{rr}$ é a espessura da camada rica em resina.
- A rigidez de penalidade de cisalhamento ( $K_s, K_t$ ) é derivada como $G / h_{rr}$ , onde $G$ é o módulo de cisalhamento da camada rica em resina.
- Esta formulação reconhece que a espessura efetiva do elemento coesivo de casca é metade da soma das espessuras das cascas adjacentes, alinhando-se com as suposições cinemáticas da teoria de Kirchhoff–Love.

Resultos Principais
O método proposto foi validado contra soluções analíticas (ex: Pagano, Mittelstedt) e modelos de elementos sólidos de alta fidelidade (C3D8 no Abaqus) através de várias configurações de laminados:

Laminados de cruzamento (cross-ply) ([0°/90°]s): O método previu com precisão as tensões normais interlaminares ( $\sigma_z$ ) e de cisalhamento ( $\tau_{yz}$ ) na interface e no plano médio, mostrando excelente concordância com as soluções analíticas.
Laminados de ângulo (angle-ply) ([+θ/−θ]s): O método capturou com sucesso a significativa tensão de cisalhamento interlaminar ( $\tau_{xz}$ ) próxima às bordas livres, correspondendo tanto aos resultados analíticos quanto aos de elementos sólidos.
Laminados quase-isotrópicos e empilhamentos gerais: A abordagem foi estendida para sequências de empilhamento complexas (ex: [45/0/−45/90]s e empilhamentos arbitrários) onde não existem soluções analíticas. Nestes casos, os resultados dos elementos estruturais demonstraram concordância próxima com os modelos de elementos sólidos, capturando precisamente singularidades de tensão e valores de pico.
Análise de Sensibilidade: Os resultados mostraram que as distribuições de tensão previstas são relativamente insensíveis a variações no módulo elástico e na razão de Poisson da camada rica em resina dentro de faixas realistas. O método permanece robusto mesmo com propriedades aproximadas da camada rica em resina.
Eficiência Computacional: Ao utilizar um único elemento de casca por camada em vez de múltiplos elementos sólidos através da espessura, o método proposto reduziu o tempo de CPU por incremento em aproximadamente 45%, mantendo uma precisão comparável na previsão de tensões.

Significância e Alegações
O artigo afirma fornecer uma abordagem dedicada para recuperar o campo de tensão 3D completo dentro de elementos de casca de ordem superior, uma capacidade anteriormente indisponível para formulações de casca de ordem elevada. Ao introduzir uma formulação de rigidez de penalidade baseada na camada rica em resina, os autores abordam a discrepância entre as suposições teóricas de espessura zero e a espessura efetiva finita dos elementos coesivos estruturais.

A significância deste trabalho reside na sua capacidade de:

Permitir a previsão precisa de tensões transversais e de borda livre sem recorrer a modelos de sólidos 3D computacionalmente caros.
Facilitar a simulação de danos combinados de plano e delaminação em malhas grossas, removendo as limitações de densidade de malha dos modelos de zona coesiva tradicionais.
Fornecer uma base física para desenvolvimentos futuros em modelos de elementos coesivos estruturais para compósitos.

Os autores concluem que este método oferece uma alternativa computacionalmente eficiente e precisa para a análise de laminados compósitos complexos, particularmente para aplicações que exigem a avaliação de iniciação de delaminação e efeitos de borda livre.

A Grande Ideia do Artigo: Vendo as Camadas Invisíveis

A "Receita Secreta" da Camada Rica em Resina

O Que Eles Descobriram

O Ganho de Velocidade

O Que Isso Significa para o Futuro

Mais como este