add a scaled identity to the rod stiffness matrix if there is no Dirichlet boundary: this regularizes the problem

[[Imported from SVN: r6886]]