A modelagem de escavações subterrâneas no FLAC3D, especialmente quando se trata de sistemas de suporte complexos, apresenta desafios significativos. Este artigo detalha um processo de modelagem passo a passo para a escavação de um túnel ou galeria, incorporando suporte combinado de concreto projetado (shotcrete) e ancoragens, utilizando comandos do FLAC3D.

Configuração Inicial do Modelo Geométrico e Material

Iniciamos com a criação de uma malha tridimensional que representa o maciço rochoso. É crucial garantir que as regiões próximas à futura escavação possuam uma discretização mais refinada para capturar com precisão os gradientes de tensão.

zone generate brick 25 20 15        ; Cria um bloco de solo/rocha com dimensões 25x20x15 unidades
zone group 'MaciçoRocha'
zone cmodel assign mohr to 'MaciçoRocha' ; Atribui o modelo constitutivo de Mohr-Coulomb
zone property bulk 6e9 shear 2.5e9 friction 32 cohesion 1e6 tension 0 to 'MaciçoRocha'

Recomenda-se iniciar a simulação com um modelo elástico linear para verificar a estabilidade e a corretude das condições de contorno, antes de prosseguir com modelos constitutivos mais complexos.

Definição das Condições de Contorno

As condições de contorno são fundamentais para simular o comportamento do maciço. O fundo do modelo é restrito verticalmente, e as faces laterais são restritas nas direções perpendiculares a elas, permitindo deformação no plano horizontal.

zone face apply velocity-z 0 range z = 0
zone face apply velocity-x 0 range x = 0 or x = 25
zone face apply velocity-y 0 range y = 0 or y = 20

Simulação da Escavação em Fases

A escavação do túnel é realizada em etapas para simular o processo construtivo e permitir que o maciço responda progressivamente à remoção de material. Este método é essencial para evitar instabilidades numéricas e obter uma representação mais realista.

def escavacao_progressiva
    local profundidade_escavada = 0.0
    local largura_tunel = 8.0
    local altura_tunel = 6.0
    local x_centro = 12.5
    local y_centro = 10.0
    
    loop passo from 1 to 5
        profundidade_escavada = 2.0 * passo
        zone delete range x = @x_centro - @largura_tunel/2 to @x_centro + @largura_tunel/2 ...
                          y = @y_centro - @altura_tunel/2 to @y_centro + @altura_tunel/2 ...
                          z = 15.0 - @profundidade_escavada to 15.0 - (@profundidade_escavada - 2.0)
        solve
        output message 'Etapa de escavação ' + string(passo) + ' concluída.'
    end_loop
end
call escavacao_progressiva

Implementação do Sistema de Suporte

Concreto Projetado (Shotcrete)

O revestimento de concreto projetado é modelado utilizando elementos de casca (struct shell). A conexão entre o revestimento e o maciço rochoso é crítica e é feita através de interfaces.

struct shell create by-zone-face range group 'ExcavatedSurface' ; Associa à superfície do túnel
struct shell property type=isotropic E=3.0e10 nu=0.25 thick=0.10 group 'RevestimentoConcreto'
interface create type=node-to-node range group 'RevestimentoConcreto'
interface property normal-stiffness=1e12 shear-stiffness=1e12 cohesion=1.2e6 friction=35

Ancoragens (Chumbadores e Cabos)

As ancoragens são representadas por elementos de cabo (struct cable). A correta definição de suas propriedades e sua interface com o maciço é vital para simular seu efeito de confinamento.

; Criação de uma ancoragem simples
struct cable create begin (10.0, 10.0, 14.5) end (8.0, 10.0, 14.5) nseg 5 group 'AncoraModelo'
struct cable property young=2.0e11 area=5e-4 bond-shear-stiff=5e10 bond-cohesion=1e6 yield-force=250e3 group 'AncoraModelo'

; Exemplo de padrão de ancoragens
struct cable create begin (10.0, 8.0, 14.0) end (8.0, 7.0, 13.0) group 'Ancoragens laterais'
struct cable create begin (10.0, 12.0, 14.0) end (8.0, 13.0, 13.0) group 'Ancoragens laterais'
struct cable create begin (12.5, 10.0, 14.0) end (12.5, 10.0, 12.0) group 'Ancoragens teto'
struct cable property young=2.1e11 area=4e-4 bond-shear-stiff=6e10 bond-cohesion=1.2e6 yield-force=300e3 group 'Ancoragens laterais' 'Ancoragens teto'

Pós-Processamento e Verificação de Resultados

Após a simulação, a visualização e análise dos resultados são cruciais para avaliar o desempenho do suporte e a estabilidade do túnel. É importante verificar a aderência entre os elementos de suporte e o maciço para garantir que as interfaces estejam funcionando corretamente.

plot create 'VisualizacaoDeslocamento'
plot add contour displacement magnitude
plot add struct cable force 'Axial' color 'red' ; Visualiza as forças axiais nas ancoragens

plot create 'TensoesAncoragem'
plot add contour struct cable stress 'Axial' range 'AncoraModelo'
plot add surface contour stress 'SigmaZ'

plot add struct shell force 'MomentX' ; Visualiza momentos no concreto projetado

A análise dos campos de deslocamento e das forças nas ancoragens e no concreto projetado permite identificar pontos de maior solicitação e subsidiar ajustes no projeto de suporte. Por exemplo, grandes deformações ou forças de ancoragem próximas ao limite de esocamento podem indicar a necessidade de aumentar o comprimento, a densidade ou o diâmetro das ancoragens.

Análise Paramétrica com FISH

O uso de scripts FISH permite automatizar análises paramétricas, como a otimização do comprimento ou da inclinação das ancoragens, testando diferantes configurações de forma eficiente.

def otimizar_comprimento_ancora
    local comprimentos_testar = (3.0, 4.0, 5.0, 6.0) ; Comprimentos de ancoragem em metros
    local ponto_fixo_x = 10.0
    local ponto_fixo_y = 10.0
    local ponto_fixo_z = 14.0

    foreach comp_atual in comprimentos_testar
        struct cable delete group 'AncorasOtimizacao' ; Remove ancoragens anteriores

        ; Recria as ancoragens com o novo comprimento
        struct cable create begin (@ponto_fixo_x, @ponto_fixo_y, @ponto_fixo_z) ...
                             end (@ponto_fixo_x - @comp_atual, @ponto_fixo_y, @ponto_fixo_z) ...
                             nseg 5 group 'AncorasOtimizacao'
        
        struct cable property young=2.1e11 area=4e-4 bond-shear-stiff=6e10 bond-cohesion=1.2e6 yield-force=300e3 group 'AncorasOtimizacao'
        
        solve
        save 'resultados_len_' + string(comp_atual) + '.sav'
        output message 'Simulação para comprimento de ' + string(comp_atual) + 'm concluída.'
    end_foreach
    output message 'Análise paramétrica de comprimento de ancoragem finalizada.'
end
call otimizar_comprimento_ancora

Para modelos complexos, é prudente salvar o modelo frequentemente e executar as simulações em pequanos incrementos de tempo, verificando a convergência. O FLAC3D 6.0 e versões posteriores oferecem ferramentas de depuração robustas, como o recurso de retrocesso automático, que facilitam a gestão de instabilidades. A combinação de modelagem cuidadosa, análise de resultados e experimentação paramétrica é a chave para o sucesso na simulação de escavações subterrâneas.