On Thu, Jan 4, 2018 at 1:12 AM, Avtar Singh <[email protected]> wrote:
> Hello Petsc Users, > > I am solving a fortran framework to simulate phase-field fracture problem. > Upto the failure point the code is working fine. But as the crack-phase > starts to propagate, there is abrupt changes in the Stiffness Matrix, so > determinant becomes nearly equal to zero. Hence, the code crashes. > I am not sure what formulation you are using, but that does not happen with the Francfort-Marrigo-Bourdin formulation I am familiar with: https://bitbucket.org/litianyi/dynamic-gradient-damage Matt > > I tried with mumps package with pc lu. Also tried, superlu, gmres and > fgmres with jacobi, bjacobi, lu. But the problem still persist. > > Can anyone suggest, Which solver and preconditioner should i use? > > Thank you > > -- > Avtar Singh > Research Scholar > Multiscale Mechanics and Multiphysics Lab > IIT Roorkee, Roorkee > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
