An Auto-Adaptive Sub-Stepping Algorithm for Phase-Field Modeling of Brittle Fracture