Better handle on negative lubrication height and zero element volume

This pull request is to better handling lubrication height and zero detJ. Before, this was done with EH and it halts the calculation. Now, they are handled in similar manner as negative detJ in mesh equation: Under-relax Newton's method in steady state problem or decrease time step in transient problem