⚡ electrical engineering

Unified Multilaminate Constitutive Model for Clay and Sand

Este estudo apresenta um modelo constitutivo multilaminado unificado que integra uma formulação de superfície de contorno com uma relação de dilatância dependente do estado para reproduzir com precisão os comportamentos mecânicos complexos, anisotrópicos e dependentes do caminho de tensões de tanto argilas quanto areias dentro de uma estrutura única e computacionalmente eficiente.

Autores originais: Behnam Ghobadi, Ehsan Taheri, Mosleh Eftekhari

Publicado 2026-09-14

📖 6 min de leitura🧠 Leitura aprofundada

CC BY 4.0

Autores originais: Behnam Ghobadi, Ehsan Taheri, Mosleh Eftekhari

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

O solo sob nossos pés raramente está imóvel. Desde o lento deslizamento de uma encosta até o tremor súbito de um terremoto, o solo é um material que se rearranja constantemente sob pressão. Para engenheiros que projetam arranha-céus, pontes ou túneis, prever exatamente como essa terra se comportará é uma questão de segurança e estabilidade. O desafio reside no fato de que o solo não é um bloco uniforme de concreto; é uma mistura complexa de minúsculos grãos, água e ar que reage de forma diferente dependendo de como foi formado, de quão úmido está e da direção em que é empurrado. Os métodos tradicionais para modelar esse comportamento frequentemente tratam o solo como uma substância simples e uniforme que reage da mesma forma em todas as direções. Embora úteis para tarefas simples, essas abordagens mais antigas têm dificuldade quando confrontadas com a realidade desordenada do mundo natural, onde as camadas de solo são inclinadas, os grãos estão alinhados e o caminho do estresse muda conforme uma estrutura se assenta.

Para resolver isso, pesquisadores buscam há muito tempo um conjunto único de regras que pudesse descrever tanto os solos argilosos e pegajosos encontrados em vales de rios quanto os solos arenosos e frouxos encontrados em desertos. O objetivo é um modelo unificado que entenda que uma areia densa e uma argila macia podem compartilhar a mesma física subjacente, mesmo que pareçam e pareçam diferentes. Esse novo entendimento exige reconhecer que o solo possui uma memória de como foi carregado no passado e que sua estrutura interna muda conforme é espremido ou cisalhado. Sem um modelo que capture essas nuances, as simulações computacionais de projetos geotécnicos podem ignorar falhas críticas, como o colapso repentino de uma encosta ou o afundamento inesperado de uma fundação.

Em um estudo recente, uma equipe de pesquisadores da Universidade Tarbiat Modares, no Irã, desenvolveu uma nova estrutura matemática projetada para preencher essas lacunas. Eles criaram um modelo constitutivo unificado capaz de simular o comportamento tanto de argila quanto de areia sob uma ampla variedade de condições. Em vez de depender de um conjunto único e rígido de equações que assume que o solo é o mesmo em todas as direções, a abordagem deles decompõe o problema em muitas partes pequenas e gerenciáveis. Imagine o solo em um único ponto como uma esfera cercada por muitas superfícies planas diferentes, cada uma voltada para uma direção ligeiramente distinta. Os pesquisadores propuseram que a reação do solo ao estresse deve ser calculada independentemente em cada uma dessas superfícies, ou "planos", antes de serem combinadas novamente para formar o quadro geral. Este método, conhecido como estrutura multilaminada, permite que o modelo considere naturalmente o fato de que o solo frequentemente se comporta de maneira diferente dependendo de qual direção está sendo empurrado, um fenômeno conhecido como anisotropia.

A equipe integrou essa abordagem direcional com uma forma sofisticada de lidar com a maneira como o solo cede e se deforma. Eles usaram um conceito chamado superfície de subcarga, que permite que o solo comece a mudar de forma e fluir mesmo antes de atingir seu ponto máximo de ruptura. Isso é crucial porque o solo real raramente espera até estar completamente estressado para começar a se deformar; ele começa a se deslocar gradualmente. Ao permitir esse movimento precoce, o modelo pode reproduzir com maior precisão a transição suave de um estado rígido e sólido para um estado mais macio e fluido. Além disso, os pesquisadores introduziram uma regra que vincula a densidade e a pressão atuais do solo à medida que ele irá expandir ou contrair quando cisalhado. Essa regra "dependente de estado" ajuda o modelo a distinguir entre areia frouxa, que tende a colapsar para dentro quando sacudida, e areia densa, que tende a estufar para fora.

Para testar sua criação, os pesquisadores realizaram simulações computacionais de testes laboratoriais padrão em cinco tipos diferentes de solo: argila de Weald, argila de Londres, argila de Kaolin de Cardiff, areia de Hostun e areia de Ottawa. Esses testes incluíram cenários onde a água podia escapar (drenado) e cenários onde a água estava presa no interior (não drenado), mimetizando condições do mundo real, como um edifício assentando em solo úmido ou um banco de areia sendo carregado rapidamente. Os resultados mostraram que o novo modelo era capaz de reproduzir as curvas complexas vistas em experimentos reais. Ele capturou a maneira como a areia densa atingiria um pico de resistência e depois amoleceria, a maneira como a areia frouxa poderia liquefazer-se subitamente e perder toda a resistência, e a maneira como diferentes argilas mudariam seu volume e pressão conforme fossem comprimidas.

Uma das descobertas mais significativas foi a capacidade do modelo de lidar com os comportamentos contrastantes de argilas e areias dentro da mesma estrutura matemática. No passado, engenheiros frequentemente tinham que usar um conjunto de regras para argila e um completamente diferente para areia. Esta nova abordagem unificada utiliza um único conjunto de parâmetros para descrever ambos, baseando-se em um "parâmetro de estado" que mede o quão longe o solo está de sua condição crítica e estável. Se o solo é frouxo, o parâmetro é positivo; se é denso, é negativo. Essa mudança simples permite que o modelo preveja se o solo irá contrair ou expandir sem a necessidade de alternar entre diferentes teorias. As simulações também previram corretamente a ocorrência de liquefação estática na areia de Ottawa frouxa, uma condição perigosa onde o solo subitamente se torna um estado fluido sob carga, um fenômeno que modelos antigos frequentemente têm dificuldade em capturar com precisão.

Os pesquisadores também demonstraram que seu modelo poderia lidar com os efeitos da sobreconsolidação, que é quando o solo foi espremido no passado e depois a pressão foi removida, deixando-o em um estado duro e denso. Quando simularam testes na argila de Kaolin de Cardiff com diferentes níveis dessa pressão passada, o modelo mostrou corretamente como o comportamento do solo mudou de espremer para dentro para empurrar para fora conforme a sobreconsolidação aumentava. Este nível de detalhe é vital para prever como o terreno existente reagirá a novas construções. O modelo alcançou tudo isso mantendo-se computacionalmente eficiente, o que significa que poderia ser usado em softwares de engenharia de grande escala sem retardar os cálculos excessivamente.

Embora o estudo tenha focado no carregamento monotônico, o que significa empurrar o solo em uma única direção sem reverter a força, os resultados sugerem uma ferramenta poderosa para entender o solo. O modelo reproduziu com sucesso o comportamento não linear de tensão-deformação, a resistência ao cisalhamento de pico e o amolecimento pós-pico observados em materiais reais. Mostrou que, ao observar o solo de muitos ângulos simultaneamente, os engenheiros podem obter uma imagem mais clara e honesta de como a terra se moverá. O trabalho não pretende ter resolvido todos os problemas da mecânica dos solos, particularmente no que diz respeito ao carregamento complexo de ida e volta, como o causado por terremotos, mas fornece uma base robusta e fisicamente significativa para futuros avanços. Ao unificar a descrição de argila e areia e considerar a natureza direcional do solo, esta pesquisa oferece uma maneira mais confiável de garantir que as estruturas construídas sobre a terra permaneçam seguras e estáveis.

Resumo Técnico: Modelo Constitutivo Multilaminado Unificado para Argila e Areia

Definição do Problema
A previsão precisa da resposta mecânica complexa de geomateriais sob diferentes trajetórias de tensão permanece um desafio significativo na modelagem constitutiva de solos. Modelos de estado crítico convencionais, como o modelo Modified Cam-Clay (MCC), frequentemente têm dificuldade em reproduzir a dependência da trajetória de tensão, a anisotropia do solo, o amolecimento por deformação (strain softening) e os comportamentos contrastantes de argilas e areias dentro de uma única estrutura. As formulações existentes frequentemente falham em capturar tensões desviatórias de pico em areias fofas e argilas normalmente consolidadas ou tendem a superestimar as tensões de ruptura sob trajetórias de tensão específicas. Além disso, embora estruturas de superfície delimitadora (bounding surface) e multilaminadas avançadas tenham sido desenvolvidas de forma independente, a integração simultânea de capacidades direcionais intrínsecas com uma formulação constitutiva unificada aplicável tanto a solos coesivos quanto granulares permanece limitada.

Metodologia
O estudo apresenta um modelo constitutivo multilaminado unificado que integra a estrutura multilaminada com uma formulação de superfície delimitadora baseada no conceito de superfície de subcarga (subloading surface). A metodologia é estruturada da seguinte forma:

Estrutura Multilaminada: O modelo aproxima a distribuição direcional de tensão e deformação discretizando uma esfera unitária em 17 planos de amostragem representativos. As respostas constitutivas são avaliadas independentemente em cada plano usando componentes de tensão normal e cisalhante locais. A resposta macroscópica é obtida através da integração numérica ponderada dessas respostas locais, proporcionando uma representação intrínseca da anisotropia induzida pela tensão e da evolução constitutiva direcional sem exigir um grande número de variáveis anisotrópicas macroscópicas adicionais.
Formulação Constitutiva Unificada: Um único conjunto de equações constitutivas é aplicado tanto para argilas quanto para areias. O modelo utiliza:
- Parâmetro de Estado: O parâmetro de estado ( $\psi$ ), definido como a diferença entre a razão de vazios atual e a razão de vazios do estado crítico na mesma pressão de confinamento, serve como uma variável de estado comum para descrever o comportamento mecânico de ambos os tipos de solo.
- Superfície de Subcarga e Superfície Delimitadora: A deformação plástica é permitida desde o início do carregamento via uma superfície de subcarga, que compartilha a mesma forma de uma superfície delimitadora. Uma regra de mapeamento radial conecta o estado de tensão atual à superfície delimitadora.
- Regra de Dilatância Generalizada: Uma regra de fluxo não associada é utilizada, regida por uma relação de dilatância generalizada dependente do estado. Esta formulação considera os efeitos combinados da razão de tensão e do parâmetro de estado, permitindo que o modelo distinga entre comportamentos contrativos e dilatantes.
- Endurecimento e Fluxo: Uma formulação avançada de endurecimento isotrópico governa a evolução do tamanho da superfície delimitadora, enquanto o conceito de superfície de subcarga garante uma transição suave entre o comportamento elástico e elastoplástico.
Implementação Numérica: Um procedimento de integração de tensão explícito incremental é adotado. As equações constitutivas são resolvidas diretamente em cada plano de amostragem usando o critério de carregamento de Hashiguchi. Enquanto o sistema macroscópico de equações lineares é resolvido via eliminação de Gauss sem iterações de Newton–Raphson, o cálculo da razão de rendimento normal ( $\gamma$ ) necessária para a formulação da superfície de subcarga é realizado usando um esquema iterativo Newton–Raphson. As matrizes constitutivas locais resultantes são transformadas e integradas para obter a matriz constitutiva macroscópica.

Principais Contribuições
A principal contribuição deste trabalho é o desenvolvimento de uma estrutura unificada que combina com sucesso as capacidades direcionais da abordagem multilaminada com uma formulação de superfície delimitadora dependente do estado. Os principais avanços técnicos incluem:

Tratamento Unificado: O modelo descreve tanto solos coesivos (argilas) quanto granulares (areias) usando um único conjunto de parâmetros e equações, baseando-se no parâmetro de estado para diferenciar o comportamento.
Anisotropia Intrínseca: Ao avaliar as leis constitutivas em 17 planos de amostragem independentes, o modelo captura naturalmente a rotação do eixo principal de tensão, a não-coaxialidade e a anisotropia induzida pela tensão.
Dilatância Dependente do Estado: A incorporação do parâmetro de estado na regra de dilatância permite a reprodução precisa de transições de fase e os comportamentos volumétricos contrastantes de solos densos/superconsolidados versus fofos/normalmente consolidados.
Integração da Superfície de Subcarga: O uso do conceito de superfície de subcarga permite a simulação de comportamento elastoplástico precoce e amolecimento por deformação sem transições abruptas.

Resultados
O desempenho preditivo do modelo foi validado contra resultados de testes triaxiais drenados e não drenados para cinco tipos distintos de solo: argila de Weald, argila de Londres, argila de Cardiff Kaolin, areia de Hostun e areia de Ottawa. As simulações numéricas demonstraram o seguinte:

Argilas: O modelo reproduziu com sucesso o comportamento não linear tensão-deformação, a deformação volumétrica e a transição do comportamento elástico para o plástico para as argilas de Weald e Londres. Para a argila de Cardiff Kaolin, o modelo capturou com precisão a influência de diferentes razões de sobreconsolidação (OCRs), incluindo a transição do comportamento contrativo para o dilatante e a evolução correspondente da poropressão excessiva.
Areias: Para a areia densa de Hostun, o modelo reproduziu o pico de resistência ao cisalhamento, o amolecimento por deformação pós-pico e a dilatação progressiva. Para a areia fofa de Ottawa, as simulações capturaram o pico de tensão desviatória, a redução subsequente em direção ao estado crítico e a ocorrência de liquefação estática sob condições não drenadas.
Desempenho Geral: Em todos os casos, o modelo reproduziuu efetivamente a evolução da trajetória de tensão, os comportamentos contrativos e dilatantes dependentes do estado e a liquefação estática, mostrando boa concordância com os dados experimentais.

Significância e Alegações
Os autores afirmam que a formulação proposta fornece uma estrutura consistente e fisicamente significativa para descrever uma ampla gama de respostas monotônicas do solo. O estudo afirma que o modelo oferece uma vantagem significativa sobre as formulações de estado crítico isotrópicas convencionais ao capturar simultaneamente o pico de resistência, o amolecimento por deformação, a dilatância dependente do estado e a evolução da trajetória de tensão, mantendo uma estrutura numérica relativamente eficiente adequada para aplicações de elementos finitos.

O artigo observa explicitamente que a validação atual é focada principalmente em carregamento triaxial monotônico. Os autores afirmam que pesquisas adicionais são necessárias para investigar o desempenho do modelo sob carregamentos cíclicos e trajetórias de carga não proporcionais complexas, rotação do eixo principal de tensão e evolução de estrutura (fabric) mais pronunciada. Eles também sugerem que trabalhos futuros devem focar em melhorias na identificação de parâmetros e na extensão do modelo para análises de elementos finitos hidro-mecânicas e cíclicas completas. O estudo conclui que o modelo estabelece uma estrutura promissora para a simulação numérica de argilas e areias, equilibrando realismo físico, versatilidade e eficiência computacional.

Mais como este