Algoritmo incremental iterativo

El algoritmo de solución es bastante flexible, ya que permite la utilización de los métodos de Newton-Raphson (NR), Newton-Raphson modificado (mNR) o un procedimiento de solución híbrido NR-mNR. Resulta claro que el ahorro en términos de esfuerzo computacional para la formación, ensamblaje y reducción de la matriz de rigidez durante el proceso iterativo puede resultar significativo al utilizar el método mNR en lugar del de NR. Sin embargo, el método mNR suele requerir un mayor número de iteraciones, lo cual puede resultar, en ocasiones, en una demanda computacional excesiva. Por esta razón, el enfoque híbrido, en el cual la matriz de rigidez es actualizada sólo durante las primeras iteraciones de un incremento de carga, suele resultar óptimo.

El procedimiento iterativo sigue los esquemas convencionales utilizados en análisis no lineales, en los cuales las fuerzas internas correspondientes a un incremento de desplazamiento son calculadas y la convergencia es verificada. Si ésta no es alcanzada, se aplican sobre la estructura las fuerzas residuales no equilibradas (resultantes de la diferencia entre el vector de carga aplicado y las fuerzas internas equilibradas) y se computa el nuevo incremento de desplazamiento. Este ciclo continúa hasta alcanzar la convergencia (mensaje del marcador del registro igual a Converg) o el número máximo de iteraciones definido por el usuario (mensaje del marcador del registro igual a Max_Ite).

Por mayores detalles y precisiones respecto de los algoritmos mencionados anteriormente, se recomienda a los usuarios referirse a la bibliografía disponible, como por ejemplo los trabajos de Cook et al. [1988], Crisfield [1991], Zienkiewicz and Taylor [1991], Bathe [1996] and Felippa [2002], por citar sólo algunos.

Nota: Algunos tipos de elementos (infrmFB, infrmFBPH, infrmDBPH & mampostería) requieren un cierto número de iteraciones para alcanzar el equilibrio interno. En algunos casos, el equilibrio interno no puede ser alcanzado, lo cual es señalado por los mensajes del marcador del registro elm_inv y elm_ite. Por mayor información sobre este tema, por favor refierase aquí.