Se presenta un modelo refinado del Método de Rigidez para estructuras tridimensionales, el cual considera el comportamiento geométricamente no lineal del sistema hasta el umbral de inestabilidad. El modelo se basa en expansiones de las matrices locales de rigidez truncadas al término cúbico de la solicitación axial local y ha sido implementado computacionalmente mediante un proceso iterativo. Se presentan resultados numéricos.