Métodos de solución iterativos para estudios no lineales

Estudios estáticos no lineales

En análisis estáticos no lineales, el conjunto básico de ecuaciones que deben ser resueltas en cualquier paso de "tiempo", t+Δt, es:

t+Δt{R} - t+Δt{F} = 0,

donde:

t+Δt{R} = Vector de cargas nodales de aplicación externa

t+Δt{F} = Vector de fuerzas nodales generadas internamente.

Debido a que las fuerzas nodales internas t+Δt{F} dependen de los desplazamientos nodales en el tiempo t+Δt, t+Δt{U} se debe utilizar un método iterativo. Las siguientes ecuaciones representan el diseño básico de un esquema iterativo para solucionar las ecuaciones de equilibrio en un paso de tiempo determinado, t+Δt,

{ΔR}(i-1) = t+Δt{R} - t+Δt{F}(i-1)

t+Δt[K](i) {ΔU}(i) = {ΔR}(i-1)

t+Δt{U}(i) = t+Δt{U}(i-1) + {ΔU}(i)

t+Δt{U}(0) = t{U}; t+Δt{F}(0) = t{F}

donde:

t+Δt{R} = Vector de cargas nodales de aplicación externa

t+Δt{F}(i-1) = Vector de fuerzas nodales generadas internamente en la iteración (i)

{ΔR}(i-1) = Vector de cargas no equilibradas en la iteración (i)

{ΔU}(i) = Vector de desplazamientos nodales incrementales en la iteración (i)

t+Δt{U}(i) = Vector de los desplazamientos totales en la iteración (i)

t+Δt[K](i) = Matriz jacobiana (rigidez tangente) en la iteración (i).

Existen diferentes esquemas para realizar las iteraciones mencionadas. A continuación, se presenta una breve descripción de dos métodos de tipo Newton:

Métodos de solución iterativos: Esquema Newton-Raphson (NR)

En este esquema, la matriz de rigidez tangencial se forma y se descompone en cada iteración dentro de un paso en particular según se muestra en la figura siguiente. El método NR tiene un alto coeficiente de convergencia y su coeficiente de convergencia es cuadrático. Sin embargo, debido a que la rigidez tangencial se forma y descompone en cada iteración, lo cual puede ser demasiado costoso para modelos grandes, es recomendable utilizar otro método iterativo.

Métodos de solución iterativos: Esquema Newton-Raphson modificado (MNR)

En este esquema, la matriz de rigidez tangencial se forma y se descompone al iniciarse cada paso (o según se especifique al definir las propiedades del estudio) y se utiliza durante las iteraciones, como muestra la siguiente figura.