Modelos numéricos geodinâmicos 2D de zonas de subducção: reologias "thermal" e "strain-rate dependent" (dependência não linear). Implementação de condições de fronteira laterais periódicas e a topo do tipo "free slip". Análise de resultados (outputs de strain rate, temperatura e viscosidade) e discussão.