Utiliza-se um algoritmo padrão de Newton-Raphson (NR) completo para obter a solução de um problema não linear do Método dos Elementos Finitos.
Em geral, o algoritmo NR não converge com frequência quando a carga total é aplicada num único passo. Uma abordagem habitual, que é também a utilizada neste caso, consiste em aplicar a carga sequencialmente em vários incrementos e utilizar o resultado do incremento de carga anterior como ponto de partida da solução de Newton do incremento seguinte. Para este fim, foi implementado um algoritmo de controlo de carga sobre o Newton-Raphson. Caso as iterações NR não convirjam, o incremento de carga atual é reduzido para metade do seu valor e as iterações NR são repetidas.
Um segundo objetivo do algoritmo de controlo de carga é determinar a carga crítica, que corresponde a determinados "critérios de paragem" – nomeadamente a deformação máxima no betão, o deslizamento máximo nos elementos de aderência, o deslocamento máximo nos elementos de ancoragem e a deformação máxima nas barras de armadura. A carga crítica é determinada através do método da bissecção. Caso o critério de paragem seja excedido em qualquer ponto do modelo, os resultados do último incremento de carga são descartados e é calculado um novo incremento com metade da dimensão do anterior. Este processo é repetido até que a carga crítica seja determinada com uma determinada tolerância de erro.
Para o betão, o critério de paragem foi definido como uma deformação de 5% à compressão (ou seja, cerca de uma ordem de grandeza superior à deformação de rotura real do betão) e de 7% à tração nos pontos de integração dos elementos de casca. À tração, o valor foi definido de modo a permitir que a deformação limite da armadura, que é normalmente de cerca de 5% sem considerar o enrijecimento à tração, seja atingida primeiro. À compressão, o valor foi escolhido de entre várias alternativas por ser suficientemente elevado para que os efeitos do esmagamento sejam visíveis nos resultados, mas suficientemente reduzido para não causar demasiados problemas de estabilidade numérica.
Para a armadura, o critério de paragem é definido em termos de tensões. Uma vez que as tensões são modeladas na fenda, o critério à tração corresponde à resistência à tração da armadura, tendo em conta o coeficiente de segurança. O mesmo valor é utilizado para o critério à compressão.
O critério de paragem nos elementos de aderência e nas molas de ancoragem é α·δumax, em que δumax é o deslizamento máximo utilizado nas verificações normativas e α = 10.
