A robust monolithic solver for phase-field fracture integrated with fracture energy based arc-length method and under-relaxation.