Dynamics and Large Deformation System: not converging after static initialization

[chunhuizhao478]

Hi all, I would like to simulate dynamic system with large deformation, what the boundary achieved is $\sigma_{xx} = 120 MPa$ and $\sigma_{yy} = 135 MPa$, see the attached minimal example:

[Mesh]
    [./msh]
        type = GeneratedMeshGenerator
        dim = 2
        nx = 100
        ny = 100
        xmin = -60000
        xmax = 60000
        ymin = -60000
        ymax = 60000
    []
    [./extranodeset1]
        type = ExtraNodesetGenerator
        coord = '0 -60000 0'
        new_boundary = corner_ptr
        input = msh
    []
    displacements = 'disp_x disp_y'
    use_displaced_mesh = true
[]

[GlobalParams]
    displacements = 'disp_x disp_y'
[]


[Variables]
    [disp_x]
        order = FIRST
        family = LAGRANGE     
    []
    [disp_y]
        order = FIRST
        family = LAGRANGE    
    []
[]

[AuxVariables]
    [vel_x]
        order = FIRST
        family = LAGRANGE
    []
    [accel_x]
        order = FIRST
        family = LAGRANGE
    []
    [vel_y]
        order = FIRST
        family = LAGRANGE
    []
    [accel_y]
        order = FIRST
        family = LAGRANGE
    []
[]

[AuxKernels]
    [accel_x]
        type = NewmarkAccelAux
        variable = accel_x
        displacement = disp_x
        velocity = vel_x
        beta = 0.25
        execute_on = timestep_end
    []
    [vel_x]
        type = NewmarkVelAux
        variable = vel_x
        acceleration = accel_x
        gamma = 0.5
        execute_on = timestep_end
    []
    [accel_y]
        type = NewmarkAccelAux
        variable = accel_y
        displacement = disp_y
        velocity = vel_y
        beta = 0.25
        execute_on = timestep_end
    []
    [vel_y]
        type = NewmarkVelAux
        variable = vel_y
        acceleration = accel_y
        gamma = 0.5
        execute_on = timestep_end
    []
[]

[Kernels]
    [dispkernel_x]
        type = TotalLagrangianStressDivergence
        variable = disp_x
        component = 0
        large_kinematics = true
    []
    [dispkernel_y]
        type = TotalLagrangianStressDivergence
        variable = disp_y
        component = 1
        large_kinematics = true
    []
    [./inertia_x]
        type = InertialForce
        use_displaced_mesh = true
        variable = disp_x
        acceleration = accel_x
        velocity = vel_x
        beta = 0.25
        gamma = 0.5
    []
    [./inertia_y]
        type = InertialForce
        use_displaced_mesh = true
        variable = disp_y
        acceleration = accel_y
        velocity = vel_y
        beta = 0.25
        gamma = 0.5
    []
[]

[Materials]
    [density]
        type = GenericConstantMaterial
        prop_names = 'density'
        prop_values = '2700'
    []
    [strain]
        type = ComputeLagrangianStrain
        large_kinematics = true
        output_properties = 'deformation_gradient'
        outputs = exodus
    []
    [elastic_tensor]
        type = ComputeIsotropicElasticityTensor
        lambda = 1e10
        shear_modulus = 1e10
    []
    [compute_stress]
        type = ComputeStVenantKirchhoffStress
        large_kinematics = true
        output_properties = 'green_lagrange_strain pk2_stress'
        outputs = exodus
    []   
[]  

[Preconditioning]
    [smp]
      type = SMP
      full = true
    []
[]
  
[Executioner]
    type = Transient
    solve_type = 'NEWTON'
    # solve_type = 'PJFNK'
    start_time = 0
    end_time = 1e100
    # num_steps = 1
    l_max_its = 100
    l_tol = 1e-7
    nl_rel_tol = 1e-8
    nl_max_its = 20
    nl_abs_tol = 1e-10
    petsc_options_iname = '-ksp_type -pc_type'
    petsc_options_value = 'gmres     hypre'
    # petsc_options_iname = '-pc_type -pc_factor_shift_type'
    # petsc_options_value = 'lu       NONZERO'
    # petsc_options_iname = '-ksp_gmres_restart -pc_type -sub_pc_type'
    # petsc_options_value = '101                asm      lu'
    # petsc_options_iname = '-ksp_type -pc_type -pc_hypre_type  -ksp_initial_guess_nonzero -ksp_pc_side -ksp_max_it -ksp_rtol -ksp_atol'
    # petsc_options_value = 'gmres        hypre      boomeramg                   True        right       1500        1e-7      1e-9    '
    # automatic_scaling = true
    # nl_forced_its = 3
    line_search = 'bt'
    dt = 1e-8
    [./TimeIntegrator]
        type = NewmarkBeta
        beta = 0.25
        gamma = 0.5
        inactive_tsteps = 1
    [../]
[]

[Controls] # turns off inertial terms for the first time step
  [./period0]
    type = TimePeriod
    disable_objects = '*/vel_x */vel_y */accel_x */accel_y */inertia_x */inertia_y'
    start_time = 0.0
    end_time = 1e-3 # dt used in the simulation
  [../]
[../]

[Outputs] 
    exodus = true
    time_step_interval = 1
[]

[BCs]
    [bc_load_top_y]
        type = DirichletBC
        variable = disp_y
        value = -337.5
        boundary = top
    []
    [bc_fix_bottom_x]
        type = DirichletBC
        variable = disp_y
        value = 0
        boundary = bottom
    []
    [bc_load_left_x]
        type = DirichletBC
        variable = disp_x
        value = 213.75
        boundary = left
    []
    [bc_load_right_x]
        type = DirichletBC
        variable = disp_x
        value = -213.75
        boundary = right
    []   
    # fix ptr
    [./fix_cptr1_x]
        type = DirichletBC
        variable = disp_x
        boundary = corner_ptr
        value = 0
    []
    [./fix_cptr2_y]
        type = DirichletBC
        variable = disp_y
        boundary = corner_ptr
        value = 0
    []
[]

The solution is unable to converge after static initialization, the nonlinear residual is unable to drop below the threshold:

Time Step 0, time = 0

Time Step 1, time = 0.001, dt = 0.001
Pre-SMO residual: 4553.43
 0 Nonlinear |R| = 9.198460e+13
      0 Linear |R| = 9.198460e+13
      1 Linear |R| = 1.832397e+12
      2 Linear |R| = 1.099771e+11
      3 Linear |R| = 8.131518e+09
      4 Linear |R| = 6.987512e+08
      5 Linear |R| = 7.147934e+07
      6 Linear |R| = 7.881670e+06
 1 Nonlinear |R| = 6.985300e+13
      0 Linear |R| = 6.985300e+13
      1 Linear |R| = 1.176999e+12
      2 Linear |R| = 5.558359e+10
      3 Linear |R| = 3.994652e+09
      4 Linear |R| = 3.251245e+08
      5 Linear |R| = 3.553913e+07
      6 Linear |R| = 4.296090e+06
 2 Nonlinear |R| = 3.096218e+13
      0 Linear |R| = 3.096218e+13
      1 Linear |R| = 3.160956e+11
      2 Linear |R| = 1.175575e+10
      3 Linear |R| = 6.269621e+08
      4 Linear |R| = 4.995027e+07
      5 Linear |R| = 4.587469e+06
      6 Linear |R| = 4.920255e+05
 3 Nonlinear |R| = 2.070896e+12
      0 Linear |R| = 2.070896e+12
      1 Linear |R| = 2.259257e+10
      2 Linear |R| = 9.136488e+08
      3 Linear |R| = 5.926474e+07
      4 Linear |R| = 6.822258e+06
      5 Linear |R| = 7.688284e+05
      6 Linear |R| = 7.898013e+04
 4 Nonlinear |R| = 1.485067e+10
      0 Linear |R| = 1.485067e+10
      1 Linear |R| = 1.583929e+08
      2 Linear |R| = 6.614170e+06
      3 Linear |R| = 4.246260e+05
      4 Linear |R| = 5.145280e+04
      5 Linear |R| = 5.621947e+03
      6 Linear |R| = 5.576845e+02
 5 Nonlinear |R| = 9.039025e+05
 Solve Converged!

Time Step 2, time = 0.002, dt = 0.001
Pre-SMO residual: 903903
 0 Nonlinear |R| = 9.039025e+05
      0 Linear |R| = 9.039025e+05
      1 Linear |R| = 5.210092e+04
      2 Linear |R| = 5.679163e+03
      3 Linear |R| = 8.261564e+02
      4 Linear |R| = 1.748659e+02
      5 Linear |R| = 2.812000e+01
      6 Linear |R| = 4.871999e+00
      7 Linear |R| = 1.013468e+00
      8 Linear |R| = 2.009413e-01
      9 Linear |R| = 5.026781e-02
 1 Nonlinear |R| = 6.720198e+03
      0 Linear |R| = 6.720198e+03
      1 Linear |R| = 5.664326e+02
      2 Linear |R| = 1.076791e+02
      3 Linear |R| = 2.050985e+01
      4 Linear |R| = 4.079992e+00
      5 Linear |R| = 1.253128e+00
      6 Linear |R| = 2.933435e-01
      7 Linear |R| = 8.639430e-02
      8 Linear |R| = 2.735438e-02
      9 Linear |R| = 8.229276e-03
     10 Linear |R| = 2.375981e-03
     11 Linear |R| = 6.550166e-04
 2 Nonlinear |R| = 5.010084e+03
      0 Linear |R| = 5.010084e+03
      1 Linear |R| = 4.970719e+02
      2 Linear |R| = 9.544612e+01
      3 Linear |R| = 1.816948e+01
      4 Linear |R| = 3.691453e+00
      5 Linear |R| = 1.078413e+00
      6 Linear |R| = 2.836297e-01
      7 Linear |R| = 8.393670e-02
      8 Linear |R| = 2.520353e-02
      9 Linear |R| = 7.548594e-03
     10 Linear |R| = 2.199123e-03
     11 Linear |R| = 6.257451e-04
     12 Linear |R| = 2.019227e-04
 3 Nonlinear |R| = 4.692583e+03
      0 Linear |R| = 4.692583e+03
      1 Linear |R| = 4.107929e+02
      2 Linear |R| = 7.826963e+01
      3 Linear |R| = 1.593809e+01
      4 Linear |R| = 3.269233e+00
      5 Linear |R| = 9.755990e-01
      6 Linear |R| = 2.459756e-01
      7 Linear |R| = 7.532317e-02
      8 Linear |R| = 2.390892e-02
      9 Linear |R| = 6.912469e-03
     10 Linear |R| = 2.020222e-03
     11 Linear |R| = 5.925084e-04
     12 Linear |R| = 1.894466e-04
 4 Nonlinear |R| = 4.363519e+03
      0 Linear |R| = 4.363519e+03
      1 Linear |R| = 3.687721e+02
      2 Linear |R| = 6.885182e+01
      3 Linear |R| = 1.396455e+01
      4 Linear |R| = 2.993200e+00
      5 Linear |R| = 8.631987e-01
      6 Linear |R| = 2.299164e-01
      7 Linear |R| = 7.264178e-02
      8 Linear |R| = 2.195557e-02
      9 Linear |R| = 6.318704e-03
     10 Linear |R| = 1.880952e-03
     11 Linear |R| = 5.714754e-04
     12 Linear |R| = 1.770981e-04
 5 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
 6 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
 7 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
 8 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
 9 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
10 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
11 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
12 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
13 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
14 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
15 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
16 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
17 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
18 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
19 Nonlinear |R| = 4.324588e+03
      0 Linear |R| = 4.324588e+03
      1 Linear |R| = 3.657097e+02
      2 Linear |R| = 6.782493e+01
      3 Linear |R| = 1.377488e+01
      4 Linear |R| = 2.950402e+00
      5 Linear |R| = 8.503640e-01
      6 Linear |R| = 2.281406e-01
      7 Linear |R| = 7.245050e-02
      8 Linear |R| = 2.160866e-02
      9 Linear |R| = 6.241758e-03
     10 Linear |R| = 1.856704e-03
     11 Linear |R| = 5.682125e-04
     12 Linear |R| = 1.756610e-04
20 Nonlinear |R| = 4.324588e+03
  Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 20
 Solve Did NOT Converge!
Currently Executing
  Finished Solving                                                                       [ 14.10 s] [  177 MB]
Aborting as solve did not converge

Also if I try to add loading to the system during dynamic solving, it turns to give weird velocity and displacements results, I wonder if there are some settings I neglect to cause this issue, thanks a lot!

[GiudGiud]

Hello
For your first time step you are accepting an unconverged solution. You ll want to tighten the tolerance on the nonlinear solver just in case that solution is causing trouble for the second time step

[chunhuizhao478]

Hi Giud, I tried to tighten the tolerance but it is the same

[chunhuizhao478]

Hi @GiudGiud, I'm just curious, in the current 2D setting, how can I make sure it is under plane strain condition? Thank!

[maxnezdyur]

In your input file replace start_time = 0 to
start_time = -1e-12 in the control block and add automatic_scaling = true to the Executioner block. Static Initialization is not well supported for solid dynamics but I think the way you did it was correct.

[chunhuizhao478]

Thanks @maxnezdyur, I think this kind of solves my problem!