Convergence problem of creep model

[lucuo]

Hello Moose Development Team,
I am trying a creep model, but when the creep strain is large, the program always one step converge and then one step error Encountered inf or nan in material return mapping iterations. At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0 In 1000 iterations the residual went from -11.8316 to -nan in 'creep'.. Could you please advise on how to improve convergence under large strains? Below are the relevant settings from my input file.

[Physics/SolidMechanics/QuasiStatic]
  [all]
    strain = FINITE
    incremental = true
    add_variables = true
    generate_output = 'stress_yy strain_yy creep_strain_yy '
  []
[]
[Materials]
  [elasticity_tensor]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 5e11
    poissons_ratio = 0.3
  []
 [./creep]
  type = Creepmodel
  temperature = temp
  coefficients = '2.5e-5 6e-8'
  activation_energys = '18000 10000'
[../]
[./creep_plas]
  type = ComputeMultipleInelasticStress
  tangent_operator = elastic
  inelastic_models = 'creep '
  max_iterations = 50
  absolute_tolerance = 1e-02
[../]
[]
[BCs]

  [u_bottom_fix]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0.0
  []
  [u_yz_fix]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  []
  [u_xy_fix]
    type = DirichletBC
    variable = disp_z
    boundary = back
    value = 0.0
  []
  [temp_fix]
    type = DirichletBC
    variable = temp
    boundary = 'bottom top'
    value = 1000.0
  []
  [./Pressure]
    [./pressure]
      boundary = 'top'
      factor = -100e6
    [../]
  [../]
[]

Here is part of the program output.

Framework Information:
MOOSE Version:           git commit 4107adf72d on 2024-08-28
LibMesh Version:         
PETSc Version:           3.21.4
SLEPc Version:           3.21.1
Current Time:            Wed Jan  8 22:48:45 2025
Executable Timestamp:    Wed Jan  8 22:24:35 2025

Checkpoint:
  Wall Time Interval:      Every 3600 s
  User Checkpoint:         Disabled
  # Checkpoints Kept:      2
  Execute On:              TIMESTEP_END 

Parallelism:
  Num Processors:          1
  Num Threads:             1

Mesh: 
  Parallel Type:           replicated
  Mesh Dimension:          3
  Spatial Dimension:       3
  Nodes:                   8
  Elems:                   1
  Num Subdomains:          1

Nonlinear System:
  Num DOFs:                32
  Num Local DOFs:          32
  Variables:               { "temp" "disp_x" "disp_y" "disp_z" } 
  Finite Element Types:    "LAGRANGE" 
  Approximation Orders:    "FIRST" 

Auxiliary System:
  Num DOFs:                11
  Num Local DOFs:          11
  Variables:               "fast_neutron_flux" { "stress_yy" "strain_yy" "creep_strain_yy" } 
  Finite Element Types:    "LAGRANGE" "MONOMIAL" 
  Approximation Orders:    "FIRST" "CONSTANT" 

Execution Information:
  Executioner:             Transient
  TimeStepper:             ConstantDT
  TimeIntegrator:          ImplicitEuler
  Solver Mode:             Preconditioned JFNK
  PETSc Preconditioner:    lu mat_superlu_dist_fact: SamePattern  mat_superlu_dist_replacetinypivot: true 


Time Step 0, time = 0

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
+----------------+----------------+----------------+


Time Step 1, time = 100000, dt = 100000
 0 Nonlinear |R| = 5.000000e+07
      0 Linear |R| = 5.000000e+07
      1 Linear |R| = 3.944267e+07
      2 Linear |R| = 1.637957e-06
 1 Nonlinear |R| = 1.038691e+08
      0 Linear |R| = 1.038691e+08
      1 Linear |R| = 4.734065e+07
      2 Linear |R| = 3.953094e-04

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -1654.13 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 1, time = 50000, dt = 50000
 0 Nonlinear |R| = 5.000000e+07
      0 Linear |R| = 5.000000e+07
      1 Linear |R| = 3.826928e+07
      2 Linear |R| = 9.195750e-07
 1 Nonlinear |R| = 1.086685e+08
      0 Linear |R| = 1.086685e+08
      1 Linear |R| = 9.936241e+06
      2 Linear |R| = 2.807456e-04

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -74.0775 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 1, time = 25000, dt = 25000
 0 Nonlinear |R| = 5.000000e+07
      0 Linear |R| = 5.000000e+07
      1 Linear |R| = 3.640671e+07
      2 Linear |R| = 2.675547e-06
 1 Nonlinear |R| = 1.173146e+08
      0 Linear |R| = 1.173146e+08
      1 Linear |R| = 2.555997e+06
      2 Linear |R| = 1.435902e-04

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.46939 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 1, time = 12500, dt = 12500
 0 Nonlinear |R| = 5.000000e+07
      0 Linear |R| = 5.000000e+07
      1 Linear |R| = 3.393029e+07
      2 Linear |R| = 3.391693e-07
 1 Nonlinear |R| = 1.920391e+05
      0 Linear |R| = 1.920391e+05
      1 Linear |R| = 1.317353e+05
      2 Linear |R| = 6.661266e-08
 2 Nonlinear |R| = 4.183909e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 3.346897e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
+----------------+----------------+----------------+


Time Step 2, time = 37500, dt = 25000
 0 Nonlinear |R| = 2.762531e+07
      0 Linear |R| = 2.762531e+07
      1 Linear |R| = 7.143196e+03
      2 Linear |R| = 6.244047e-06
 1 Nonlinear |R| = 1.173032e+08
      0 Linear |R| = 1.173032e+08
      1 Linear |R| = 1.345613e+06
      2 Linear |R| = 9.048020e-02

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 2, time = 25000, dt = 12500
 0 Nonlinear |R| = 2.087095e+07
      0 Linear |R| = 2.087095e+07
      1 Linear |R| = 2.704644e+03
      2 Linear |R| = 5.444010e-06
 1 Nonlinear |R| = 1.556647e+05
      0 Linear |R| = 1.556647e+05
      1 Linear |R| = 5.730100e+04
      2 Linear |R| = 7.099486e-05
 2 Nonlinear |R| = 3.029730e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.424607e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
+----------------+----------------+----------------+


Time Step 3, time = 50000, dt = 25000
 0 Nonlinear |R| = 2.762283e+07
      0 Linear |R| = 2.762283e+07
      1 Linear |R| = 7.140029e+03
      2 Linear |R| = 5.080629e-05
 1 Nonlinear |R| = 1.172927e+08
      0 Linear |R| = 1.172927e+08
      1 Linear |R| = 1.345614e+06
      2 Linear |R| = 1.282506e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 3, time = 37500, dt = 12500
 0 Nonlinear |R| = 2.086908e+07
      0 Linear |R| = 2.086908e+07
      1 Linear |R| = 2.702685e+03
      2 Linear |R| = 4.871648e-05
 1 Nonlinear |R| = 1.556537e+05
      0 Linear |R| = 1.556537e+05
      1 Linear |R| = 5.729420e+04
      2 Linear |R| = 2.581781e-04
 2 Nonlinear |R| = 3.029417e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.424116e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
+----------------+----------------+----------------+


Time Step 4, time = 62500, dt = 25000
 0 Nonlinear |R| = 2.762035e+07
      0 Linear |R| = 2.762035e+07
      1 Linear |R| = 7.140039e+03
      2 Linear |R| = 7.507184e-05
 1 Nonlinear |R| = 1.172822e+08
      0 Linear |R| = 1.172822e+08
      1 Linear |R| = 1.345618e+06
      2 Linear |R| = 9.754059e-02

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 4, time = 50000, dt = 12500
 0 Nonlinear |R| = 2.086721e+07
      0 Linear |R| = 2.086721e+07
      1 Linear |R| = 2.702690e+03
      2 Linear |R| = 6.907979e-05
 1 Nonlinear |R| = 1.556427e+05
      0 Linear |R| = 1.556427e+05
      1 Linear |R| = 5.728751e+04
      2 Linear |R| = 5.998405e-04
 2 Nonlinear |R| = 3.029105e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.423635e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
+----------------+----------------+----------------+


Time Step 5, time = 75000, dt = 25000
 0 Nonlinear |R| = 2.761787e+07
      0 Linear |R| = 2.761787e+07
      1 Linear |R| = 7.140123e+03
      2 Linear |R| = 1.735079e-04
 1 Nonlinear |R| = 1.172717e+08
      0 Linear |R| = 1.172717e+08
      1 Linear |R| = 1.345624e+06
      2 Linear |R| = 1.619276e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 5, time = 62500, dt = 12500
 0 Nonlinear |R| = 2.086534e+07
      0 Linear |R| = 2.086534e+07
      1 Linear |R| = 2.702732e+03
      2 Linear |R| = 1.555194e-04
 1 Nonlinear |R| = 1.556318e+05
      0 Linear |R| = 1.556318e+05
      1 Linear |R| = 5.728098e+04
      2 Linear |R| = 1.631024e-03
 2 Nonlinear |R| = 3.028796e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.423142e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
+----------------+----------------+----------------+


Time Step 6, time = 87500, dt = 25000
 0 Nonlinear |R| = 2.761539e+07
      0 Linear |R| = 2.761539e+07
      1 Linear |R| = 7.140077e+03
      2 Linear |R| = 1.715144e-02
 1 Nonlinear |R| = 1.172612e+08
      0 Linear |R| = 1.172612e+08
      1 Linear |R| = 1.345625e+06
      2 Linear |R| = 1.265038e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 6, time = 75000, dt = 12500
 0 Nonlinear |R| = 2.086347e+07
      0 Linear |R| = 2.086347e+07
      1 Linear |R| = 2.702736e+03
      2 Linear |R| = 3.683948e-04
 1 Nonlinear |R| = 1.556208e+05
      0 Linear |R| = 1.556208e+05
      1 Linear |R| = 5.727429e+04
      2 Linear |R| = 3.772965e-03
 2 Nonlinear |R| = 3.028476e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.422648e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
|   7.500000e+04 |   1.078495e-03 |   1.278495e-03 |
+----------------+----------------+----------------+


Time Step 7, time = 100000, dt = 25000
 0 Nonlinear |R| = 2.761292e+07
      0 Linear |R| = 2.761292e+07
      1 Linear |R| = 7.140013e+03
      2 Linear |R| = 8.744307e-03
 1 Nonlinear |R| = 1.172507e+08
      0 Linear |R| = 1.172507e+08
      1 Linear |R| = 1.345627e+06
      2 Linear |R| = 2.042160e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 7, time = 87500, dt = 12500
 0 Nonlinear |R| = 2.086160e+07
      0 Linear |R| = 2.086160e+07
      1 Linear |R| = 2.702695e+03
      2 Linear |R| = 4.448420e-03
 1 Nonlinear |R| = 1.556099e+05
      0 Linear |R| = 1.556099e+05
      1 Linear |R| = 5.726742e+04
      2 Linear |R| = 5.655096e-03
 2 Nonlinear |R| = 3.028167e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.422181e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
|   7.500000e+04 |   1.078495e-03 |   1.278495e-03 |
|   8.750000e+04 |   1.258245e-03 |   1.458245e-03 |
+----------------+----------------+----------------+


Time Step 8, time = 112500, dt = 25000
 0 Nonlinear |R| = 2.761044e+07
      0 Linear |R| = 2.761044e+07
      1 Linear |R| = 7.140005e+03
      2 Linear |R| = 1.648070e-03
 1 Nonlinear |R| = 1.172402e+08
      0 Linear |R| = 1.172402e+08
      1 Linear |R| = 1.345631e+06
      2 Linear |R| = 1.761336e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 8, time = 100000, dt = 12500
 0 Nonlinear |R| = 2.085973e+07
      0 Linear |R| = 2.085973e+07
      1 Linear |R| = 2.702700e+03
      2 Linear |R| = 8.285261e-03
 1 Nonlinear |R| = 1.555989e+05
      0 Linear |R| = 1.555989e+05
      1 Linear |R| = 5.726074e+04
      2 Linear |R| = 4.378199e-03
 2 Nonlinear |R| = 3.027848e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.421676e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
|   7.500000e+04 |   1.078495e-03 |   1.278495e-03 |
|   8.750000e+04 |   1.258245e-03 |   1.458245e-03 |
|   1.000000e+05 |   1.437994e-03 |   1.637994e-03 |
+----------------+----------------+----------------+


Time Step 9, time = 125000, dt = 25000
 0 Nonlinear |R| = 2.760797e+07
      0 Linear |R| = 2.760797e+07
      1 Linear |R| = 7.140087e+03
      2 Linear |R| = 2.019849e-03
 1 Nonlinear |R| = 1.172297e+08
      0 Linear |R| = 1.172297e+08
      1 Linear |R| = 1.345636e+06
      2 Linear |R| = 1.212490e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 9, time = 112500, dt = 12500
 0 Nonlinear |R| = 2.085786e+07
      0 Linear |R| = 2.085786e+07
      1 Linear |R| = 2.702714e+03
      2 Linear |R| = 1.968096e-03
 1 Nonlinear |R| = 1.555880e+05
      0 Linear |R| = 1.555880e+05
      1 Linear |R| = 5.725409e+04
      2 Linear |R| = 8.983486e-03
 2 Nonlinear |R| = 3.027538e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.421197e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
|   7.500000e+04 |   1.078495e-03 |   1.278495e-03 |
|   8.750000e+04 |   1.258245e-03 |   1.458245e-03 |
|   1.000000e+05 |   1.437994e-03 |   1.637994e-03 |
|   1.125000e+05 |   1.617743e-03 |   1.817743e-03 |
+----------------+----------------+----------------+


Time Step 10, time = 137500, dt = 25000
 0 Nonlinear |R| = 2.760550e+07
      0 Linear |R| = 2.760550e+07
      1 Linear |R| = 7.140022e+03
      2 Linear |R| = 3.671282e-03
 1 Nonlinear |R| = 1.172193e+08
      0 Linear |R| = 1.172193e+08
      1 Linear |R| = 1.345639e+06
      2 Linear |R| = 2.377454e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 10, time = 125000, dt = 12500
 0 Nonlinear |R| = 2.085599e+07
      0 Linear |R| = 2.085599e+07
      1 Linear |R| = 2.702691e+03
      2 Linear |R| = 4.810857e-03
 1 Nonlinear |R| = 1.555771e+05
      0 Linear |R| = 1.555771e+05
      1 Linear |R| = 5.724729e+04
      2 Linear |R| = 1.579306e-02
 2 Nonlinear |R| = 3.027224e+01
 Solve Converged!

Outlier Variable Residual Norms:
  disp_y: 2.420726e+01

Postprocessor Values:
+----------------+----------------+----------------+
| time           | creep_strain_yy| strain_yy      |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
|   1.250000e+04 |   1.797493e-04 |   3.797494e-04 |
|   2.500000e+04 |   3.594985e-04 |   5.594986e-04 |
|   3.750000e+04 |   5.392477e-04 |   7.392478e-04 |
|   5.000000e+04 |   7.189969e-04 |   9.189970e-04 |
|   6.250000e+04 |   8.987461e-04 |   1.098746e-03 |
|   7.500000e+04 |   1.078495e-03 |   1.278495e-03 |
|   8.750000e+04 |   1.258245e-03 |   1.458245e-03 |
|   1.000000e+05 |   1.437994e-03 |   1.637994e-03 |
|   1.125000e+05 |   1.617743e-03 |   1.817743e-03 |
|   1.250000e+05 |   1.797492e-03 |   1.997492e-03 |
+----------------+----------------+----------------+


Time Step 11, time = 150000, dt = 25000
 0 Nonlinear |R| = 2.760303e+07
      0 Linear |R| = 2.760303e+07
      1 Linear |R| = 7.139995e+03
      2 Linear |R| = 2.400820e-02
 1 Nonlinear |R| = 1.172088e+08
      0 Linear |R| = 1.172088e+08
      1 Linear |R| = 1.345641e+06
      2 Linear |R| = 2.017448e-01

Encountered inf or nan in material return mapping iterations.
At element 0 _qp=0 Coordinates (x,y,z)=(0.211325, 0.211325, 0.211325) block=0
In 1000 iterations the residual went from -5.4689 to -nan in 'creep'.

To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 11, time = 137500, dt = 12500
 0 Nonlinear |R| = 2.085412e+07
      0 Linear |R| = 2.085412e+07
      1 Linear |R| = 2.702677e+03
      2 Linear |R| = 1.138983e-02
 1 Nonlinear |R| = 1.555662e+05
      0 Linear |R| = 1.555662e+05
      1 Linear |R| = 5.724054e+04
      2 Linear |R| = 2.217140e-02
 2 Nonlinear |R| = 3.026916e+01
 Solve Converged!

[GiudGiud]

Hello

We cannot take huge time steps sometimes. Simply because the element would be so deformed after such a step that it no longer makes geometrical sense.

So my recommendation is to maintain reasonable time steps through the simulation.

Setting preset=false on displacements BCs can help if the deformed elements are near the boundary conditions