⚡ electrical engineering

A discrete‑natured 3D DEM framework for spatial fracture analysis of unidirectional CFRP composites

Este estudo apresenta um novo framework de método de elementos discretos 3D com um esquema de empacotamento cúbico que trata células unitárias como elementos contínuos para modelar eficazmente a fratura transversal em compósitos de CFRP unidirecionais, revelando que, embora a distribuição das fibras afete minimamente a resposta mecânica global, a fração volumétrica e o raio das fibras governam criticamente a rigidez, a resistência e a localização de danos.

Autores originais: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

Publicado 2026-09-11

📖 5 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Dechao Sun, Xiaoxuan Ding, Shuqi Zhang, Jianlin Liu, Zewen Gu, Xiaonan Hou

Artigo original sob licença CC BY 4.0 (https://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

A engenharia moderna depende fortemente de materiais que são fortes, porém leves, como os polímeros reforçados com fibra de carbono. Esses compósitos são feitos ao incorporar finos e rígidos fios de carbono em uma cola plástica mais macia, criando uma estrutura que pode suportar cargas pesadas em aplicações aeroespaciais e automotivas. No entanto, como qualquer material, eles têm uma fraqueza: quando puxados lateralmente, a cola plástica entre as fibras pode rachar. Esse tipo de rachadura, conhecido como fratura transversal, frequentemente começa como pequenas falhas invisíveis que crescem e se conectam, acabando por fazer com que toda a estrutura perca sua resistência. Compreender exatamente como essas rachaduras se formam e se movem através do complexo labirinto de fibras e cola é difícil porque o dano ocorre em uma escala pequena demais para ser vista facilmente com ferramentas padrão, e as rachaduras não viajam em linhas retas. Elas torcem, viram e interagem com as fibras no espaço tridimensional, tornando difícil prever quando uma peça pode falhar.

Para resolver esse quebra-cabeça, os pesquisadores Dechao Sun e sua equipe da Universidade do Petróleo da China, trabalhando com Xiaonan Hou da Universidade de Lancaster, desenvolveram uma nova maneira de simular esses materiais em um computador. Em vez de tratar o material como um bloco sólido e liso, eles o dividiram em uma vasta coleção de partículas individuais minúsculas que se tocam e se empurram umas contra as outras. Essa abordagem, chamada método de elementos discretos, permite que o computador observe como o material se comporta como se fosse feito de milhões de blocos de construção microscópicos. A equipe criou um modelo tridimensional específico de um cubo deste material compósito, preenchendo-o com fibras de carbono virtuais e resina epóxi. Eles então conceberam uma nova regra matemática para conectar a rigidez das pequenas conexões entre essas partículas à resistência geral do material. Uma inovação fundamental em seu método é que o modelo gera automaticamente a tendência do material de estufar quando comprimido — uma propriedade conhecida como razão de Poisson — simplesmente observando como as partículas se movem, em vez de forçar o computador a usar um número pré-estabelecido.

Quando os pesquisadores testaram este novo modelo, descobriram que ele podia reproduzir com precisão o comportamento real do material. A simulação previu que o material se esticaria e quebraria de maneiras que correspondiam a experimentos laboratoriais existentes e a modelos computacionais bidimensionais anteriores. No entanto, a visão tridimensional revelou algo que os modelos antigos perderam: as rachaduras não atravessam o material de uma só vez como uma faca através de um pão de forma. Em vez disso, a frente da rachadura avança de forma irregular. Em alguns pontos, a rachadura avança rapidamente, enquanto em outros, ela fica para trás, criando um caminho irregular e sinuoso através da espessura do material. Essa progressão desigual, que os pesquisadores chamam de tortuosidade da frente da rachadura, é um detalhe crucial que apenas uma visão tridimensional completa pode capturar.

A equipe então usou essa ferramenta para investigar como diferentes arranjos das fibras afetam a resistência do material. Eles testaram se o posicionamento aleatório das fibras dentro da cola plástica alterava a resistência do material. Descobriram que, embora o caminho exato que uma rachadura percorresse dependesse fortemente de onde as fibras estavam localizadas, a resistência e a rigidez geral do material permaneciam quase as mesmas, independentemente do arranjo. Isso sugere que, para prever o desempenho geral do compósito, o padrão aleatório específico das fibras importa menos do que se pensava anteriormente. No entanto, a quantidade de fibra compactada no material fez uma diferença significativa. Aumentar o conteúdo de fibra de 50 por cento para 68 por cento tornou o material mais rígido e forte, mas também o tornou mais quebradiço, fazendo com que falhassem com um estiramento menor. Os pesquisadores observaram que, quando as fibras eram compactadas demais, não havia cola suficiente entre elas para absorver o estresse, levando a rachaduras rápidas.

Finalmente, os pesquisadores observaram o tamanho das próprias fibras, mantendo a quantidade de fibra constante. Eles descobriram que o raio das fibras, ou quão grossas elas eram, influenciava como o material lidava com o estresse. Um raio de fibra de 3,5 micrômetros provou ser a configuração mais eficiente, oferecendo o melhor equilíbrio entre a transferência de carga e a capacidade de a cola esticar. As simulações mostraram que o tamanho das fibras ditava o espaçamento entre elas, o que, por sua vez, controlava onde o estresse se concentrava e onde as rachaduras começariam. Ao mapear as forças entre as partículas virtuais, a equipe pôde ver exatamente como a carga passava da cola para as fibras e onde o material finalmente cedia. Este novo arcabouço tridimensional fornece uma ferramenta poderosa para engenheiros entenderem a mecânica oculta dos materiais compósitos, oferecendo uma imagem mais clara de como essas estruturas complexas falham e de como elas podem ser projetadas para durar mais.

Resumo Técnico: Uma Estrutura de DEM 3D de Natureza Discreta para Análise de Fratura Espacial de Compósitos CFRP Unidirecionais

Definição do Problema
Polímeros reforçados com fibra de carbono (CFRPs) são críticos em aplicações aeroespaciais e automotivas, mas sofrem com uma tolerância limitada ao dano, particularmente no que diz respeito ao trincamento transversal. Este modo de falha inicia microtrincas que se propagam através da espessura da camada, levando à degradação da rigidez e ao comprometimento estrutural. Embora as abordagens numéricas existentes, como a Mecânica do Dano Contínuo (CDM), Modelos de Zona Coesiva (CZM) e o Método de Elementos Finitos Estendido (XFEM), tenham avançado o campo, elas frequentemente dependem de caminhos de trinca predefinidos ou enfrentam dificuldades com a convergência e a complexidade de múltiplas trincas 3D que se interceptam. Além disso, os estudos existentes de Método de Elementos Discretos (DEM) são predominantemente limitados a análises de seção transversal 2D, falhando em capturar a continuidade da trinca através da espessura, a conectividade de dano 3D e a redistribuição de tensão espacial inerente às microestruturas heterogêneas. Existe uma lacuna em modelos de 3D DEM explícitos capazes de simular o trincamento transversal em CFRPs unidirecionais sem depender de propriedades anisotrópicas equivalentes ou estruturas homogeneizadas.

Metodologia
Os autores propõem uma nova estrutura de método de elementos discretos tridimensional (3D DEM) para investigar a fratura transversal em compósitos de epóxi reforçados com fibra de carbono. A metodologia é construída sobre os seguintes pilares:

Modelagem Discreta: Um Elemento de Volume Representativo (RVE) cúbico de 49 µm é construído usando um esquema de empacotamento cúbico 3D. A microestrutura consiste em 38 fibras circulares (raio de 3,5 µm) imersas em uma matriz de epóxi, com interações modeladas via molas de força central ao longo de seis direções ortogonais.
Calibração de Natureza Discreta: Diferente dos modelos tradicionais de partículas ligadas que tratam as ligações como elementos contínuos que exigem calibração simultânea das rigidezes normal ( $k_n$ $k_{n}$ ) e tangencial ( $k_s$ $k_{s}$ ), este estudo introduz uma inovadora relação de equivalência de energia. O modelo trata cada célula unitária como um elemento contínuo onde a razão de Poisson macroscópica é gerada automaticamente pela deformação da célula, em vez de ser prescrita.
- A estrutura deriva relações analíticas explícitas (Eqs. 41 e 42) vinculando a rigidez normal microscópica ( $k_n$ ) ao módulo de Young macroscópico ( $E$ ) e à razão de Poisson ( $\nu$ ).
- Esta abordagem elimina a necessidade de rigidez tangencial ( $k_s$ ) no processo de calibração, baseando-se exclusivamente na rigidez normal para satisfazer a equivalência de energia entre as forças de contato microscópicas e a deformação do contínuo macroscópico.
Validação: O modelo calibrado é validado contra dados experimentais e resultados de DEM 2D. As simulações são conduzidas sob carregamento de tração por deslocamento quase-estático para observar a evolução do dano, a propagação de trincas e as respostas tensão-deformação macroscópicas.

Principais Resultados
O estudo realiza investigações paramétricas sobre os efeitos da distribuição aleatória de fibras, fração volumétrica de fibras ( $V_f$ ) e raio da fibra ( $r$ ):

Validação do Modelo: O 3D DEM reproduz com sucesso o comportamento mecânico macroscópico do epóxi puro e de RVEs de CFRP, com o módulo de Young e a razão de Poisson simulados desviando-se em menos de 2% dos valores alvo. O modelo captura a continuidade da trinca através da espessura e a tortuosidade da frente da trinca, mostrando que a propagação da trinca não é simultânea ao longo do comprimento da fibra, ao contrário das trajetórias planas em modelos 2D.
Distribuição de Fibras: Embora as distribuições aleatórias de fibras alterem significamente os locais de iniciação de trincas locais e os caminhos finais das trincas (por exemplo, flutuando entre fibras densamente compactadas), elas têm uma influência negligenciável nas respostas mecânicas globais. RVEs estatisticamente equivalentes exibem comportamento tensão-deformação e evolução do módulo de Young quase idênticos.
Fração Volumétrica de Fibras ( $V_f$ ): O aumento de $V_f$ (de 50% para 68%) aumenta a rigidez transversal e a resistência de pico, mas reduz a deformação de ruptura, levando a uma resposta pós-pico mais frágil. Um $V_f$ mais alto induz deficiência de resina, promovendo concentração de tensão localizada e dano dominado pela matriz de forma precoce.
Raio da Fibra: Para uma fração volumétrica fixa ( $V_f = 60\%$ ), o raio da fibra governa o equilíbrio entre a transferência de tensão interfacial e a acomodação de deformação da matriz. A configuração com raio de fibra de $r = 3,5$ µm demonstrou a capacidade de carga mais eficiente e a maior rigidez em comparação com $r = 3,2$ µm e $r = 3,8$ µm.

Significância e Alegações
O artigo afirma que a estrutura 3D DEM proposta oferece uma ferramenta eficaz para análise de dano progressivo multiescala em compósitos de CFRP unidirecionais. Sua principal significância reside em:

Calibração de Natureza Discreta: Fornecer um mapeamento rigoroso, baseado em energia, de parâmetros micro para propriedades macro que respeita a natureza descontínua do DEM enquanto gera a razão de Poisson naturalmente.
Análise de Fratura Espacial: Capturar com sucesso mecanismos complexos de falha 3D, incluindo penetração de trinca através da espessura, tortuosidade da frente da trinca e coalescência de danos, que são inalcançáveis com simplificações 2D.
Insight Mecanístico: Elucidar que, embora a distribuição de fibras dite os padrões de fratura locais, a fração volumétrica de fibras é o fator dominante que controla o desempenho mecânico global.

Os autores concluem que, embora o modelo atual aborde efetivamente o trincamento transversal em laminados unidirecionais, trabalhos futuros são necessários para incorporar layups de fibras multidirecionais, relações constitutivas refinadas da interface fibra-matriz e estratégias de acoplamento de dano multiescala.

Mais como este