Un algoritm standard Newton-Raphson (NR) complet este utilizat pentru a găsi soluția unei probleme neliniare de tip MEF.
În general, algoritmul NR nu converge adesea atunci când întreaga încărcare este aplicată într-un singur pas. O abordare uzuală, care este utilizată și aici, este de a aplica încărcarea secvențial, în mai multe incremente, și de a folosi rezultatul incrementului de încărcare anterior pentru a începe soluția Newton a celui următor. În acest scop, a fost implementat un algoritm de control al încărcării peste metoda Newton-Raphson. În cazul în care iterațiile NR nu converg, incrementul de încărcare curent este redus la jumătate din valoarea sa, iar iterațiile NR sunt reluate.
Un al doilea scop al algoritmului de control al încărcării este de a găsi încărcarea critică, care corespunde anumitor „criterii de oprire” – în mod specific deformația maximă în beton, alunecarea maximă în elementele de aderență, deplasarea maximă în elementele de ancoraj și deformația maximă în barele de armătură. Încărcarea critică este determinată folosind metoda bisecției. În cazul în care criteriul de oprire este depășit oriunde în model, rezultatele ultimului increment de încărcare sunt eliminate, iar un nou increment, de jumătate din mărimea celui anterior, este calculat. Acest proces este repetat până când încărcarea critică este găsită cu o anumită toleranță de eroare.
Pentru beton, criteriul de oprire a fost stabilit la o deformație de 5% la compresiune (adică, aproximativ un ordin de mărime mai mare decât deformația reală de cedare a betonului) și 7% la întindere, la punctele de integrare ale elementelor de tip shell. La întindere, valoarea a fost stabilită astfel încât să permită atingerea mai întâi a deformației limită în armătură, care este de obicei în jur de 5% fără a lua în considerare participarea betonului întins. La compresiune, valoarea a fost aleasă dintre mai multe alternative, ca fiind suficient de mare pentru ca efectele strivirii să fie vizibile în rezultate, dar suficient de mică pentru a nu cauza prea multe probleme legate de stabilitatea numerică.
Pentru armătură, criteriul de oprire este definit în termeni de tensiuni. Deoarece sunt modelate tensiunile la fisură, criteriul la întindere corespunde rezistenței la întindere a armăturii, ținând cont de coeficientul de siguranță. Aceeași valoare este utilizată pentru criteriul la compresiune.
Criteriul de oprire în elementele de aderență și în resoartele de ancoraj este α·δumax, unde δumax este alunecarea maximă utilizată în verificările conform codului, iar α = 10.
