Multiscale Calculation of Equivalent Elastic Modulus for Rock Mass Based on Perturbation Theory

Resumo Técnico: Cálculo Multiescala do Módulo de Elasticidade Equivalente para Maciços Rochosos Baseado na Teoria de Perturbação

Definição do Problema
A determinação precisa do módulo de elasticidade equivalente é crítica para a avaliação da segurança e estabilidade em engenharia de maciços rochosos. As metodologias atuais enfrentam limitações significativas: testes laboratoriais em amostras de rocha intacta falham em considerar as juntas e fissuras macro e microestruturais inerentes aos maciços rochosos de engenharia, enquanto testes de campo são frequentemente não representativos devido ao seu escopo espacial limitado. Além disso, os métodos de redução tradicionais usados para escalar resultados laboratoriais para condições de campo dependem de fatores de redução pouco confiáveis e não conseguem capturar o processo de resposta mecânica durante a falha do maciço rochoso. Há uma necessidade reconhecida de uma abordagem de acoplamento multiescala sistemática que estabeleça a ponte entre a mecânica da rocha intacta e o comportamento heterogêneo de rochas danificadas.

Metodologia
Os autores propõem uma abordagem de cálculo multiescala fundamentada na teoria de perturbação, especificamente projetada para considerar o dano inicial em maciços rochosos. A metodologia trata o maciço rochoso com dano inicial como um sistema de perturbação onde um pequeno parâmetro ($\delta$) reflete o processo de expansão multiescala das fissuras internas.

As principais etapas metodológicas incluem:

• Estrutura Teórica: O tensor de deformação causado pela meso-heterogeneidade é utilizado como uma "ponte" conectando as escalas macro e micro. O campo de deslocamento é submetido à expansão assintótica, separando variáveis em coordenadas macroscópicas ($x$) e vizinhanças microscópicas ($y$).
• Equações Governantes: O estudo deriva as equações governantes acoplando a equação diferencial de equilíbrio a uma equação constitutiva que incorpora um fator de dano elástico microestrutural ($D^\delta$). O dano inicial é representado por um tensor de dano elástico inicial.
• Homogeneização: Ao resolver as equações de movimento sob diferentes escalas de dimensões físicas, os autores derivam o módulo de elasticidade macro equivalente homogeneizado ($E_{mnkl}^\delta$). Isso envolve a resolução dos autofunções de tensor ($\chi$) utilizando o método de elementos finitos baseado em princípios variacionais.
• Implementação Numérica: A abordagem foi implementada utilizando MATLAB para cálculos iterativos do coeficiente de dano e COMSOL para modelagem de elementos finitos. O modelo atualiza a carga ($Q$) e o deslocamento ($u$) em cada etapa para refletir o estado de dano evolutivo.
• Validação: A abordagem foi validada através de dois casos primários:
  1. Validação Laboratorial: Testes de compressão uniaxial e tração diametral (Brasil) foram conduzidos em amostras de calcário coletadas de uma mina de metal profunda no sul da China. As amostras foram pré-danificadas em graus variados (0,25$\sigma_{cs}$ a 0,75$\sigma_{cs}$) e caracterizadas via velocidade de onda P e varredura CT.
  2. Validação de Campo: O método foi aplicado a corpos solidificados por injeção em um estaleiro de uma mina de chumbo-zinco, onde o material injetado foi tratado como um maciço rochoso com fissuras de juntas irregulares.

Principais Contribuições

• Modelo de Perturbação Multiescala: O artigo estabelece uma estrutura teórica onde o tensor de deformação decorrente da meso-heterogeneidade serve como o mecanismo de conexão entre as escalas micro e macro, permitindo o cálculo do módulo de elasticidade equivalente considerando o dano inicial.
• Evolução Dinâmica do Dano: Diferente dos modelos convencionais que frequentemente tratam o módulo de elasticidade como uma constante, esta abordagem reflete o desenvolvimento dinâmico do módulo de elasticidade equivalente e a evolução de ambos os danos, macro e micro, durante o processo de carregamento.
• Melhoria da Precisão: O estudo demonstra que o cálculo multiescala proposto produz resultados mais próximos dos dados reais de testes laboratoriais e de campo em comparação com os métodos convencionais de simulação de elementos finitos.

Resultos

• Achados Laboratoriais:
  • Verificou-se que a velocidade da onda P ($V_p$) diminui conforme o dano inicial aumenta, confirmando sua utilidade como um indicador de dano.
  • As varreduras CT revelaram que um maior dano inicial correlaciona-se com o aumento da porosidade e da densidade de fissuras.
  • Teste de Tração Diametral (Brasil): As curvas carga-deslocamento exibiram uma função convexa estritamente inferior. Os resultados do cálculo multiescala alinharam-se mais proximamente com as curvas experimentais do que as simulações convencionais, particularmente na região pré-pico.
  • Teste de Compressão Uniaxial: As curvas tensão-deformação são funções côncavas monotonicamente crescentes. A abordagem multiescala previu com precisão o módulo de elasticidade equivalente, com valores mais próximos dos resultados de teste do que as simulações convencionais. A evolução do dano mostrou um aumento gradual inicialmente, seguido por um desenvolvimento rápido e estabilização.
• Achados de Campo:
  • No teste de compressão de campo dos corpos solidificados por injeção, a abordagem multiescala previu com precisão a curva tensão-deformação até o pico da tensão de compressão uniaxial.
  • O módulo de elasticidade equivalente calculado para os corpos solidificados por injeção coincidiu com os dados de teste de campo, validando a aplicabilidade do método a cenários de engenharia complexos envolvendo agregados de minério quebrado e lama de cimento.
• Limitações Observadas: O estudo observou que, embora os resultados pré-pico fossem precisos, as curvas tensão-deformação pós-pico para os corpos solidificados por injeção divergiram dos resultados de campo. Essa discrepância é atribuída à ausência de um parâmetro de amolecimento por deformação (strain-softening) no modelo, que é necessário para capturar o comportamento pós-falha de materiais que diferem de rochas frágeis profundas.

Significância e Alegações
O artigo alega que a abordagem de cálculo multiescala proposta fornece um método teórico robusto para resolver problemas complexos em engenharia de rochas. Ao considerar o dano inicial e utilizar a teoria de perturbação, o método consegue estabelecer a ponte entre a heterogeneidade microestrutural e o comportamento mecânico macroscópico. Os autores afirmam que esta abordagem oferece uma alternativa mais precisa aos métodos de simulação convencionais para determinar o módulo de elasticidade equivalente de maciços rochosos. Eles sugerem que, se uma ponte adequada entre as escalas macro e micro puder ser identificada, esta metodologia poderá ser estendida para determinar parâmetros mecânicos de outros materiais heterogêneos. O trabalho é apresentado como uma referência para melhorar a avaliação da estabilidade de maciços rochosos em projetos de engenharia, particularmente onde o dano inicial e os efeitos de escala são fatores críticos.