On modeling fracture of soft polymers

Este artigo propõe e valida um modelo de dano baseado em um valor crítico de trabalho de tensão para prever a fratura em polímeros moles rate-dependentes, demonstrando sua eficácia em diversos materiais, geometrias e condições de carregamento ao capturar tanto respostas frágeis quanto dúcteis.

O essencial

Imagine que os polímeros macios (como borrachas, géis e até alguns tecidos biológicos) são como uma multidão de pessoas segurando as mãos em uma grande rede. Às vezes, essas mãos são unidas por um aperto forte e permanente (como em uma borracha de pneu); outras vezes, são unidas por um aperto que pode soltar e se reconectar rapidamente (como em um gel que se auto-cura).

O grande problema que os cientistas enfrentam é: como prever exatamente quando e como essa rede vai se romper?

Até agora, as ferramentas usadas para medir a "resistência à quebra" desses materiais funcionavam bem apenas para borrachas duras, mas falhavam miseravelmente com materiais mais moles e "viscosos" (que escorrem como mel). A resistência parecia mudar dependendo de quão rápido você puxava o material ou de qual formato a peça tinha. Era como tentar medir o peso de uma nuvem usando uma balança de banheiro: o resultado dependia de como você segurava a nuvem.

A Grande Descoberta: O "Trabalho de Estresse" Crítico

Os autores deste artigo, Aditya Konale e Vikas Srivastava, propuseram uma nova maneira de pensar. Em vez de olhar apenas para a força ou para a energia de uma forma complexa, eles sugerem medir o "Trabalho de Estresse" (Stress Work).

Pense nisso como a quantidade total de "suor" ou esforço que o material precisa gastar para se romper completamente.

• A Analogia do Elástico: Imagine que você tem um elástico. Se você puxar devagar, ele estica muito antes de quebrar. Se você puxar rápido, ele estica pouco e estala. O "Trabalho de Estresse" é a soma total da energia que você gastou puxando até o momento da quebra.
• A Descoberta: Eles descobriram que, para uma certa classe desses materiais, esse "suor total" é um número quase fixo. Não importa se você puxou rápido ou devagar, ou se a peça era redonda ou quadrada; o material sempre precisa gastar a mesma quantidade de energia para se romper. É como se o material tivesse um "orçamento de energia" fixo para a sua vida.

O Modelo: Uma Rede com Dois Mecanismos

Para prever como a quebra acontece, eles criaram um modelo matemático (um "simulador" no computador) que vê o material como tendo dois sistemas de defesa:

1. O Sistema de Atrito (Viscoso): Como o mel ou o xarope. Se você puxa devagar, as moléculas têm tempo para deslizar e se reorganizar, absorvendo o estresse.
2. O Sistema de Rede (Elástico): Como as molas. Se você puxa rápido, elas não têm tempo de se reorganizar e esticam até o limite.

O modelo deles combina esses dois mundos. Ele consegue prever dois comportamentos fascinantes que os modelos antigos não conseguiam:

• O Perfil de Trompete (Materiais Viscosos): Quando você estica um gel muito rápido, a ponta da fissura não fica pontiaguda. Ela fica arredondada (como um trompete) e, às vezes, forma uma segunda pequena fissura à frente. É como se o material tentasse "amortecer" o golpe antes de ceder.
• O Perfil Parabólico (Borrachas Elásticas): Quando você estica uma borracha comum, a fissura avança de forma limpa e curva, como uma parábola, sem tentar se reorganizar.

Por que isso é importante?

1. Previsão Universal: Agora, podemos usar uma única "regra" para prever a quebra de borrachas de pneu, géis de cura de feridas, adesivos inteligentes e até tecidos biológicos.
2. Segurança e Design: Isso ajuda engenheiros a criar materiais mais seguros. Se você sabe exatamente quanto "suor" o material aguenta, pode projetar dispositivos médicos ou proteções contra impactos que não falham inesperadamente.
3. Do Macro para o Micro: O modelo é tão bom que, olhando para o número macroscópico (o "orçamento de energia" total), eles conseguiram estimar a energia necessária para quebrar as ligações químicas microscópicas dentro do material. Foi como deduzir o peso de um grão de areia medindo o peso de uma montanha inteira.

Resumo em uma frase

Os autores criaram um "termômetro de energia" universal que mede quanto esforço um material macio precisa gastar para se romper, permitindo prever com precisão como ele vai falhar, seja ele uma borracha elástica ou um gel viscoso, independentemente de quão rápido você o puxa.

Resumo técnico

1. O Problema

Os polímeros moles (incluindo elastômeros, géis, tecidos biológicos e polímeros viscosos) são materiais ubíquos com propriedades que variam de independentes a dependentes da taxa de carregamento (rate-dependent).

• Limitações Atuais: As medidas convencionais de tenacidade à fratura (baseadas na taxa de liberação de energia ou intensidade de tensão na ponta da trinca) não são únicas para materiais dependentes da taxa. Elas variam significativamente com o perfil de carregamento e a geometria do espécime.
• Falta de Modelos Gerais: Trabalhos existentes sobre modelagem de fratura em polímeros moles dependentes da taxa são limitados a geometrias pré-trincadas específicas ou focam apenas em elastômeros. Não há um modelo consensual e geral capaz de prever a iniciação, crescimento e coalescência de trincas em geometrias arbitrárias e complexas para uma ampla classe de polímeros (desde elastômeros até polímeros altamente viscosos).
• Desafio Específico: Polímeros com ligações cruzadas dinâmicas (como o PBS - poliborosiloxano) exibem comportamentos extremos: fluidez em escalas de tempo longas e comportamento elástico em altas taxas de carregamento, com respostas de fratura distintas (ex.: perfis de fratura em "trombeta" vs. parabólicos).

2. Metodologia

Os autores propõem e desenvolvem um modelo unificado baseado em três pilares principais:

A. Nova Medida de Resistência à Fratura ($W_{cr}$)

• Propõem que um valor crítico de trabalho de tensão ($W_{cr}$) pode ser usado como uma medida quase única de resistência à fratura para uma classe específica de polímeros moles.
• $W_{cr}$ é definido como a densidade de energia total (por unidade de volume de referência) necessária para a falha completa de um ponto material arbitrário sob taxas de carregamento suficientes para iniciar o dano.
• Diferente da tenacidade à fratura tradicional (J/m²), $W_{cr}$ tem unidades de J/m³ e é hipoteticamente independente da geometria do espécime e da história de carregamento, dependendo apenas do modo de carregamento (tração vs. compressão).

B. Modelo de Dano Generalizado (Gradient-Damage)

• Desenvolveram um modelo de dano baseado no método do poder virtual e em uma estrutura de dano de gradiente generalizado.
• Decomposição Multiplicativa: O modelo utiliza a decomposição multiplicativa do gradiente de deformação em partes elásticas e plásticas, considerando múltiplos mecanismos de deformação paralelos (M mecanismos).
• Mecanismos Específicos para Polímeros:
  • Mecanismo A (Resistência Intermolecular): Resistência elástica e dissipativa devido a interações intermoleculares (viscosidade).
  • Mecanismo B (Resistência da Rede): Estiramento de subcadeias dentro de ligações cruzadas dinâmicas (reversíveis).
• Critério de Iniciação de Dano: A iniciação do dano é governada por um critério energético baseado na contribuição energética do trabalho de tensão ($\psi^+_{cr}$), que representa a fração de $W_{cr}$ associada à energia armazenada. O dano cresce quando a energia livre deformacional atinge este valor crítico.
• Evolução do Dano: O crescimento do dano é modelado através de uma equação de evolução que inclui um termo de gradiente (parâmetro de comprimento $l$) para regularizar a zona de processo de fratura e evitar dependência de malha.

C. Implementação Numérica

• O modelo foi implementado no software de elementos finitos ABAQUS usando elementos definidos pelo usuário (VUEL para o campo de dano e VUMAT para o comportamento mecânico).
• Utilizou-se uma formulação explícita devido às escalas de tempo curtas envolvidas na fratura rápida.
• Foi desenvolvida uma procedura numérica eficiente para atualizar as variáveis de estado dependentes do histórico (como o estiramento das subcadeias formadas em tempos anteriores) sem armazenar todo o histórico de deformação.

3. Contribuições Principais

1. Proposta de $W_{cr}$: Demonstração experimental e teórica de que o trabalho de tensão crítico é uma medida robusta e independente de geometria para resistência à fratura em polímeros moles.
2. Modelo Unificado: Desenvolvimento de um modelo constitutivo capaz de descrever desde elastômeros (comportamento elástico) até polímeros viscosos (comportamento fluido/viscoelástico) usando a mesma estrutura matemática, apenas ajustando os parâmetros de relaxação de ligações cruzadas.
3. Predição de Fenômenos Físicos Complexos: O modelo consegue prever:
   • Transições entre respostas frágeis e dúcteis.
   • Perfis de fratura distintos: o perfil "trombeta" (com formação de um segundo entalhe) observado em polímeros viscosos e o perfil parabólico simples em elastômeros.
   • Fenômenos de "flow-fracture" (transição entre fluxo e fratura) em materiais com ligações dinâmicas.
4. Conexão Micro-Macro: O modelo permite inferir quantidades microestruturais a partir de parâmetros macroscópicos. Os autores estimaram a energia de dissociação de subcadeias a partir de $\psi^+_{cr}$, validando o resultado com experimentos de reologia oscilatória.

4. Resultados e Validação

O modelo foi validado contra dados experimentais de diversos materiais e configurações:

• Materiais Testados:
  • PBS (Poliborosiloxano): Polímero viscoso com ligações cruzadas dinâmicas.
  • EPDM, SBR, NR: Elastômeros (borracha de etileno-propileno, borracha de estireno-butadieno, borracha natural).
  • EPS25 Vitrimer: Polímero com ligações cruzadas dinâmicas.
• Geometrias e Condições:
  • Ensaios de tração uniaxial homogênea (para calibração de $W_{cr}$).
  • Ensaios de tração com entalhe (single edge notch, center notch).
  • Ensaios de cisalhamento puro (pure shear).
  • Configurações complexas com insertos rígidos.
  • Variações de taxas de carregamento (de 2 mm/s a 60 mm/s).
• Desempenho:
  • O modelo previu com precisão as curvas força-deslocamento e os pontos de falha para todas as geometrias e materiais, sem re-calibração para cada novo caso (apenas usando parâmetros calibrados no ensaio de tração homogênea).
  • Validação de Perfis de Fratura: A simulação reproduziu corretamente a transição do perfil de fratura em "trombeta" (para PBS em alta velocidade) para o perfil parabólico (para elastômeros), explicando o fenômeno pela dependência da taxa de tensão local na ponta da trinca.
  • Predição de Parada de Trinca: No caso do vitrimer EPS25, o modelo previu corretamente a nucleação, propagação e subsequente parada (arrest) da trinca em taxas de carregamento intermediárias, devido à reconfiguração da rede que reduz a energia disponível abaixo do crítico.

5. Significado e Conclusão

Este trabalho oferece uma solução robusta para um problema de longa data na mecânica dos sólidos: a falta de um modelo unificado para fratura em materiais moles dependentes da taxa.

• Unificação: O modelo elimina a necessidade de abordagens separadas para elastômeros e polímeros viscosos.
• Capacidade Preditiva: Ao contrário de métodos variacionais de campo de fase que minimizam energia sem critérios explícitos de iniciação (e que falham em capturar superfícies de falha experimentais), o modelo proposto usa critérios explícitos baseados em tensão/energia, permitindo prever superfícies de falha complexas e geometrias arbitrárias.
• Impacto Microestrutural: A capacidade de extrair parâmetros físicos microestruturais (como energia de ligação) de testes macroscópicos abre novas portas para o design de materiais e a compreensão de mecanismos de falha em escala molecular.
• Aplicabilidade: O modelo é crucial para o desenvolvimento de aplicações de engenharia que utilizam polímeros moles, como proteção contra impactos, adesivos, sensores autocicatrizantes e dispositivos biomédicos, onde a previsão precisa da falha sob diferentes condições de carregamento é essencial.