Added max_bond_breaks_per_iteration to original peridynamics_physics

It had been added to all the other peridynamics_physics, but not
the original. Now the user can control how many bonds can break
per iteration. This is helpful in situations where cracks propagate
quickly breaking a lot of bonds. In some cases the linear solves
die if too many bonds break. This is a way to throttle back the damage
per iteration.
Peridynamics_Physics: take dof out of lin sys if all bonds are broken

This seems to resolve the issues with the solver failing in the middle
of a simulation