Multiphysics Modeling of Thermo-Viscoelastic Damage in Functionally Graded Abradable Coatings with Probabilistic Geometric Tolerance Analysis

Resumo Técnico: Modelagem Multifísica de Dano Termo-Viscoelástico em Revestimentos Abradáveis Funcionalmente Graduados com Análise Probabilística de Tolerância Geométrica

Definição do Problema
Revestimentos abradáveis funcionalmente graduados são críticos para o controle da folga da ponta das palhetas em selos de labirinto de motores de aeronaves. No entanto, sua durabilidade é governada por uma interação complexa entre o amolecimento viscoelástico dependente da temperatura, dano progressivo, modulação microestrutural induzida pela deposição e tolerâncias geométricas. Os modelos existentes frequentemente tratam esses efeitos de forma independente ou assumem gradientes de material idealizados e suaves em geometrias nominais. Essa abordagem corre o risco de subestimar o início do dano porque falha em capturar concentrações de tensão locais decorrentes da heterogeneidade microestrutural (por exemplo, variações periódicas de rigidez devido à deposição por spray térmico) e a amplificação desses efeitos pela manufatura e tolerâncias geométricas de montagem. O artigo aborda a lacuna na previsão de como esses fatores acoplados influenciam a confiabilidade, observando que geometrias extremas raras podem dominar o risco, mas são invisíveis para análises determinísticas.

Metodologia
O estudo propõe uma estrutura unificada multifísica–probabilística que integra quatro drivers acoplados:

1. Termo-Viscoelasticidade: Modelada usando uma representação de série de Prony (Maxwell generalizada) com parâmetros dependentes da temperatura. A formulação considera a eigendeformação térmica e o descompasso de expansão revestimento–substrato.
2. Dano Progressivo: Uma lei de dano escalar dirigida por tensão (tipo Kachanov–Lemaitre) governa a degradação da rigidez. O modelo utiliza uma taxa de saturação de dano (expoente positivo) em vez de uma taxa geradora de ruptura, tratando o dano como uma métrica macroscópica de perda de rigidez.
3. Arquitetura do Material: O revestimento é modelado como um material funcionalmente graduado (FGM) com um gradiente contínuo através da espessura, sobreposto a uma modulação senoidal periódica para representar a heterogeneidade microestrutural induzida pela deposição.
4. Propagação de Tolerância Probabilística: As tolerâncias geométricas (desvios de folga/offset) são modeladas como variáveis aleatórias Gaussianas. Elas são propagadas através do solver multifísico completo usando simulações de Monte Carlo (800 realizações) para gerar distribuições de resposta.

Implementação Numérica e Verificação
O domínio resolvido é uma coluna de revestimento unidimensional local através da espessura, impulsionada por histórias prescritas de temperatura e deformação normal equivalente. O solver emprega um esquema de integração temporal incremental explícito.

• Verificação: A implementação é rigorosamente verificada através de estudos de refinamento espacial e temporal (mostrando convergência espacial próxima de segunda ordem e temporal de primeira ordem), testes de unidade de forma fechada (validando o núcleo viscoelástico, expansão térmica e cinética de dano), convergência de Monte Carlo aninhada com intervalos de confiança de Wilson e uma superfície de resposta polinomial validada por hold-out.
• Limitações de Escopo: O artigo distingue explicitamente a verificação numérica da validação experimental. O caso simulado é um "benchmark genérico" usando parâmetros consistentes com a literatura para sistemas abradáveis baseados em AlSi, não sendo um modelo identificado para um revestimento específico. O solver exclui mecânica de contato, fricção, remoção de material e profundidade de desgaste; ele resolve para tensão e dano sob excitação prescrita, não para o processo de abrasão em si.

Principos Resultados

1. Efeito da Modulação Microestrutural: A modulação periódica de propriedades aumenta significamente o dano ao final do ciclo em comparação com um gradiente monotônico. Ela concentra a tensão máxima dentro de cristas de modulação localizadas (especificamente próximo à interface revestimento–substrato, em torno de profundidade normalizada $\zeta \approx 0,96$), acelerando a cinética de dano. O dano médio para o caso modulado foi 58% maior do que o caso monotônico ao final do ciclo.
2. Impacto das Tolerâncias Geométricas: A propagação de tolerâncias geométricas simétricas através do solver acoplado produz uma distribuição de dano fortemente assimétrica à direita. Embora o dano médio permaneça abaixo do nível de classificação ($D_{crit} = 0,10$), a cauda da distribuição leva a uma probabilidade de excedência significativa ($P_e = 0,076$).
3. Não Linearidade e Falha Determinística: A probabilidade de excedência para o gradiente modulado é mais do que o dobro da do gradiente monotônico ($0,076$ vs. $0,035$). Este comportamento não pode ser capturado por uma análise determinística na geometria nominal, que falha em prever a "massa de excedência" na cauda da distribuição.
4. Sensibilidade: A análise de sensibilidade de parâmetros indica que a temperatura de permanência (dwell temperature) é o fator dominante (contribuição de 86%), seguida pelo tempo de relaxação principal (12%). A influência da tolerância geométrica é quantificada diretamente através da dispersão do conjunto em vez de um índice escalar único.

Significância e Alegações
O artigo alega que sua principal contribuição reside no tratamento unificado da modulação microestrutural periódica e das tolerâncias geométricas dentro de um único solver de dano termo-viscoelástico. Ao acoplar esses efeitos, a estrutura revela interações não lineares que análises desacopladas perdem.

• Design Orientado à Confiabilidade: A estrutura vincula os parâmetros de deposição (amplitude/comprimento de onda da modulação) e as faixas de tolerância diretamente às probabilidades quantificadas de excedência de dano. Isso apoia a triagem de arquiteturas de revestimento e especificações de tolerância orientada à confiabilidade.
• Insights Informados pela Manufatura: Os resultados demonstram que a gradação funcional sozinha não garante melhor durabilidade se a heterogeneidade induzida pela deposição e a variabilidade dimensional forem ignoradas. O estudo fornece uma cadeia quantitativa processo–estrutura–propriedade–desempenho, permitindo que projetistas avaliem como a variabilidade da manufatura afeta a confiabilidade em serviço.
• Benchmarking: O trabalho serve como um benchmark numericamente verificado para estudos futuros, distinguindo explicitamente a metodologia computacional verificada da necessidade de futura validação experimental (especificamente a 400°C e para previsão de desgaste).

Os autores concluem que, embora a estrutura seja demonstrada para selos de labirinto de turbinas, sua estrutura matemática é geral o suficiente para ser aplicada a outros sistemas graduados ou em camadas (ex: revestimentos de barreira térmica) onde a heterogeneidade induzida pela deposição e a variabilidade dimensional governam conjuntamente o desempenho.