Added method for computing effective stiffness but something is wrong in the expression. Also internal force vector differs for nlgeom = 1,-1 - different initial stresses