Adding the hessian term in variational form (broken for now) and spliting the displacement into multiple terms