tested against febio solutions.
manual newton solver implemented but convergence very slow
petsc nonlinear solver is available...
merged
good. saved kinematics just in case.