EigenStrain contribution #28558

Emmanuel-324 · 2024-09-06T14:24:44Z

Emmanuel-324
Sep 6, 2024

I have two scripts, one with eingenstrain and the other without contribution.

 [Materials]
    [elasticity_tensor_phase0]
      type = ComputeElasticityTensorCP
      C_ijkl = '2.906e5 1.87e5 1.87e5 2.906e5 1.87e5 2.906e5 1.142e5 1.142e5 1.142e5'
      fill_method = symmetric9
      euler_angle_1 = 0.0
      euler_angle_2 = 0.0
      euler_angle_3 = 0.0
      base_name = phase0
    []
    [stress_phase0]
      type = ComputeMultipleCrystalPlasticityStress
      crystal_plasticity_models = 'trial_xtalpl_phase0'
      tan_mod_type = exact
      rtol = 1e-08
      base_name = phase0
    []
    [trial_xtalpl_phase0]
      type = CrystalPlasticityKalidindiUpdate
      number_slip_systems = 48
      slip_sys_file_name = input_slip_sys_bcc48.txt
      crystal_lattice_type = BCC
      resistance_tol = 0.01
      r = 1.4             
      h = 6000            
      t_sat = 598.5        
      gss_a = 1.5         
      ao = 0.001           
      xm = 0.017             
      gss_initial = 465.5
      base_name = phase0
    []

    [./strain_phase0]
      type = ComputeFiniteStrain
      displacements = 'disp_x disp_y'
      base_name = phase0
      eigenstrain_names = 'eigenstrain_phase0'
    [../]
    
    [./eigen_strain_phase0]
      type = ComputeEigenstrain
      base_name = phase0
      eigen_base = '1 1 1 0 0 0'
      prefactor = '0.2224 0.9493 0.2224'
      eigenstrain_name = 'eigenstrain_phase0'
    [../]

  
    [elasticity_tensor_phase1]
      type = ComputeElasticityTensorCP
      C_ijkl = '2.721e5 1.69e5 1.69e5 2.721e5 1.69e5 2.721e5 1.31e5 1.31e5 1.31e5'
      fill_method = symmetric9
      euler_angle_1 = 0.0
      euler_angle_2 = 0.0
      euler_angle_3 = 0.0
      base_name = phase1
    []
    [stress_phase1]
      type = ComputeMultipleCrystalPlasticityStress
      crystal_plasticity_models = 'trial_xtalpl_phase1'
      tan_mod_type = exact
      rtol = 1e-08
      base_name = phase1
    []
    [trial_xtalpl_phase1]
      type = CrystalPlasticityKalidindiUpdate
      number_slip_systems = 12
      slip_sys_file_name = input_slip_sys.txt
      crystal_lattice_type = FCC
      resistance_tol = 0.01
      r = 1.0             
      h = 6000            
      t_sat = 598.5        
      gss_a = 1.5         
      ao = 0.001           
      xm = 0.017             
      gss_initial = 465.5 
      base_name = phase1
    []
    [./strain_phase1]
      type = ComputeFiniteStrain
      displacements = 'disp_x disp_y'
      base_name = phase1
    [../]
  
    
    # Switching functions for each phase
    [./h0]
      type = SwitchingFunctionMultiPhaseMaterial
      phase_etas = eta0
      all_etas = 'eta0 eta1'
      h_name = h0
    [../]
    [./h1]
      type = SwitchingFunctionMultiPhaseMaterial
      phase_etas = eta1
      all_etas = 'eta0 eta1'
      h_name = h1
    [../]
   
    # Generate the global stress from the phase stresses
    [./global_stress]
      type = MultiPhaseStressMaterial
      phase_base = 'phase0 phase1'
      h          = 'h0     h1'
    [../]   
 []

Below is my ComputeEigenstrain.C with modifications

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "ComputeEigenstrain.h"

registerMooseObject("SolidMechanicsApp", ComputeEigenstrain);
registerMooseObject("SolidMechanicsApp", ADComputeEigenstrain);

template <bool is_ad>
InputParameters
ComputeEigenstrainTempl<is_ad>::validParams()
{
  InputParameters params = ComputeEigenstrainBase::validParams();
  params.addClassDescription("Computes a constant Eigenstrain with anisotropic scaling factors."); 

  params.addRequiredParam<std::vector<Real>>(
      "eigen_base", "Vector of values defining the constant base tensor for the Eigenstrain");

  // Set default values for the prefactor to {1.0, 1.0, 1.0} if not provided, ensuring isotropic scaling
  params.addParam<std::vector<Real>>(
      "prefactor", std::vector<Real>{1.0, 1.0, 1.0}, 
      "Vector of values defining the variable dependence in each principal direction (x, y, z). The vector should normalize to 1.");

  return params;
}

template <bool is_ad>
ComputeEigenstrainTempl<is_ad>::ComputeEigenstrainTempl(const InputParameters & parameters)
  : ComputeEigenstrainBaseTempl<is_ad>(parameters),
   _prefactor(this->template getParam<std::vector<Real>>("prefactor")) // Changed to get prefactor as vector
{
  // Ensure that _prefactor has exactly 3 values (for x, y, z)
  mooseAssert(_prefactor.size() == 3, "Prefactor must have exactly three components."); // Added check for vector size

  // Calculate the norm of the _prefactor vector
  Real norm = std::sqrt(_prefactor[0] * _prefactor[0] + 
                        _prefactor[1] * _prefactor[1] + 
                        _prefactor[2] * _prefactor[2]);

  // Check if the norm of the prefactor vector is approximately 1
  if (!MooseUtils::absoluteFuzzyEqual(norm, 1.0,1e-4))
     this->paramError("prefactor", "Prefactor should have a norm of one.");

  _eigen_base_tensor.fillFromInputVector(this->template getParam<std::vector<Real>>("eigen_base"));
}

template <bool is_ad>
void
ComputeEigenstrainTempl<is_ad>::computeQpEigenstrain()
{
  // Using anisotropic scaling factors from _prefactor
  for (unsigned int i = 0; i < 3; ++i)
  {
    _eigenstrain[_qp](i, i) = _eigen_base_tensor(i, i) * _prefactor[i]; // Apply different scaling to each diagonal component
  }
}

Yet the simulation results at the same timestep for both scripts remain the same in terms of stress or stain.

Below are snippets of one timestep.(The results show timestep-stress in x direction-von misses stress- strain ) In that order.
(eigen strain)
(no eigen strain)

What could be the problem?
@tophmatthews @GiudGiud

GiudGiud · 2024-09-08T14:15:31Z

GiudGiud
Sep 8, 2024
Collaborator

Hello

I dont know. I would check the units, see if there isnt an order of magnitude mistake there, it is fairly common.
The eigenstrain base value of 1 seems small to me, is that what is being used in other examples?

4 replies

Emmanuel-324 Sep 8, 2024
Author

I thought the eigenstrain base works as the voghit notations. Hence my reason for "1 1 1 0 0 0' is the one unit in the x y and z directions.

GiudGiud Sep 8, 2024
Collaborator

could be. Not my field at this point

Emmanuel-324 Sep 9, 2024
Author

I would try without the normalization too. And if you can tag anyone who can help with this I would really appreciate.

GiudGiud Sep 9, 2024
Collaborator

@laagesen ?

laagesen · 2024-09-09T19:25:48Z

laagesen
Sep 9, 2024

you are right that eigenstrain base follows the Voigt convention, so '1 1 1 0 0 0' is dilatational eigenstrain that is uniform in x, y, z directions. However I'm not sure why you are modifying ComputeEigenstrain or why you are modifying to make it possible to specify 3 different prefactors. If you want an eigenstrain that is different in different directions you can just specify how much you want for each direction in eigen_base- are you trying to do something different from that?

2 replies

Emmanuel-324 Sep 9, 2024
Author

Yes, I am trying to get an eigenstrain that is different in each direction. @laagesen

I want to add this eigenstrain for one of my precipitates I found in a paper:

What would be the best way to represent this eigenstrain variant(picture above) in my script.

This thread explains my situation better and how we got here : #28474 (comment)

laagesen Sep 9, 2024

you should just be able to specify that as
eigen_base = '0.0067 0.0286 0.0067 0 0 0'
or
eigen_base = '0.0067 0.0286 0.0067'
within the source code, RankTwoTensor.h lines 654-689 shows a bunch of examples of how the tensor is created based on how eigen_base is specified (though a call to the function fillFromInputVector)

sapitts · 2024-09-09T20:28:41Z

sapitts
Sep 9, 2024
Collaborator

If you want to use an eigenstrain with one of the update-based crystal plasticity models, you should use the crystal plasticity-specific base class, with multiplicative decomposition of the deformation tensor that includes an eigenstrain specific decomposition tensor; see the base class.

You may be able to use the existing ComputeCrystalPlasticityThermalEigenstrain by providing a "dummy" constant auxvariable, set equal to 1, instead of a temperature variable, and using the values from your reference as the arguments to the thermal_expansion_coefficients parameter.

16 replies

Emmanuel-324 Sep 16, 2024
Author

If you want to use an eigenstrain with one of the update-based crystal plasticity models, you should use the crystal plasticity-specific base class, with multiplicative decomposition of the deformation tensor that includes an eigenstrain specific decomposition tensor; see the base class.

You may be able to use the existing ComputeCrystalPlasticityThermalEigenstrain by providing a "dummy" constant auxvariable, set equal to 1, instead of a temperature variable, and using the values from your reference as the arguments to the thermal_expansion_coefficients parameter.

Hi @sapitts, I implemented the eigenstrain model as you suggested but it seems the result after the simulation is no different from my script without the eigenstrain contribution. I have attached my full input script

[Mesh]
  [file]
     type = FileMeshGenerator
     file = Test_5_two_phase_va_out.e-s171
     use_for_exodus_restart = true
   []
 []
 
 [Variables]
  [./disp_x]
  [../]
  [./disp_y]
  [../]

  # order parameter 0
  [./eta0]
    initial_from_file_var = eta0
  [../]
  # order parameter 1
  [./eta1]
    initial_from_file_var = eta1
  [../]
[]

[AuxVariables]
  [./vonmises]
    order = CONSTANT
    family = MONOMIAL
  [../]
  [./stress_xx]
    order = CONSTANT
    family = MONOMIAL
  [../]
 [./dummy_temperature]
   family = MONOMIAL
   order = CONSTANT
 []
[]

[AuxKernels]
  [./vonmises]
    type = RankTwoScalarAux
    rank_two_tensor = stress
    variable = vonmises
    scalar_type = VonMisesStress
    execute_on = timestep_end
  [../]
  [./stress_xx]
    type = RankTwoAux
    rank_two_tensor = stress
    variable = stress_xx
    index_j = 0
    index_i = 0
    execute_on = timestep_end
  [../]
  
   [./set_dummy_temperature]
     type = ConstantAux
     variable = dummy_temperature
     value = 1.0  
   [../]
[]

[BCs]
  [symmy]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [symmx]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [tdisp]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = '0.1*t'
  []
[]

[Materials]
   [elasticity_tensor_phase0]
     type = ComputeElasticityTensorCP
     C_ijkl = '2.906e5 1.87e5 1.87e5 2.906e5 1.87e5 2.906e5 1.142e5 1.142e5 1.142e5'
     fill_method = symmetric9
     euler_angle_1 = 0.0
     euler_angle_2 = 0.0
     euler_angle_3 = 0.0
     base_name = phase0
   []
   [stress_phase0]
     type = ComputeMultipleCrystalPlasticityStress
     crystal_plasticity_models = 'trial_xtalpl_phase0'
     tan_mod_type = exact
     rtol = 1e-08
     base_name = phase0
   []
   [trial_xtalpl_phase0]
     type = CrystalPlasticityKalidindiUpdate
     number_slip_systems = 12
     slip_sys_file_name = input_slip_sys.txt
     crystal_lattice_type = BCC
     resistance_tol = 0.01
     r = 1.4             
     h = 6000            
     t_sat = 598.5        
     gss_a = 1.5         
     ao = 0.001           
     xm = 0.017             
     gss_initial = 600
     base_name = phase0
   []
   [phase0_eigenstrain_0]
     type = ComputeCrystalPlasticityThermalEigenstrain
     eigenstrain_name = phase0_eigenstrain_0
     deformation_gradient_name = phase0_deformation_gradient_0
     temperature = dummy_temperature
     thermal_expansion_coefficients = '0.0067 0.0286 0.0067'
     base_name = phase0
   []
   

   [./strain_phase0]
     type = ComputeFiniteStrain
     displacements = 'disp_x disp_y'
     base_name = phase0
   [../]
 
   [elasticity_tensor_phase1]
     type = ComputeElasticityTensorCP
     C_ijkl = '2.721e5 1.69e5 1.69e5 2.721e5 1.69e5 2.721e5 1.31e5 1.31e5 1.31e5'
     fill_method = symmetric9
     euler_angle_1 = 0.0
     euler_angle_2 = 0.0
     euler_angle_3 = 0.0
     base_name = phase1
   []
   [stress_phase1]
     type = ComputeMultipleCrystalPlasticityStress
     crystal_plasticity_models = 'trial_xtalpl_phase1'
     tan_mod_type = exact
     rtol = 1e-08
     base_name = phase1
   []
   [trial_xtalpl_phase1]
     type = CrystalPlasticityKalidindiUpdate
     number_slip_systems = 12
     slip_sys_file_name = input_slip_sys.txt
     crystal_lattice_type = FCC
     resistance_tol = 0.01
     r = 1.0             
     h = 6000            
     t_sat = 598.5        
     gss_a = 1.5         
     ao = 0.001           
     xm = 0.017             
     gss_initial = 465.5 
     base_name = phase1
   []
   [./strain_phase1]
     type = ComputeFiniteStrain
     displacements = 'disp_x disp_y'
     base_name = phase1
   [../]
 
   
   # Switching functions for each phase
   [./h0]
     type = SwitchingFunctionMultiPhaseMaterial
     phase_etas = eta0
     all_etas = 'eta0 eta1'
     h_name = h0
   [../]
   [./h1]
     type = SwitchingFunctionMultiPhaseMaterial
     phase_etas = eta1
     all_etas = 'eta0 eta1'
     h_name = h1
   [../]
  
   # Generate the global stress from the phase stresses
   [./global_stress]
     type = MultiPhaseStressMaterial
     phase_base = 'phase0 phase1'
     h          = 'h0     h1'
   [../]   
[]

[Kernels]
  [./eta0_dt]
    type = TimeDerivative
    variable = eta0
  [../]
  [./eta1_dt]
    type = TimeDerivative
    variable = eta1
  [../]

  [./TensorMechanics]
    displacements = 'disp_x disp_y'
  [../]
[]

[Postprocessors]
  [./vonmises]
    type = ElementAverageValue
    variable = vonmises
   # block = 'ANY_BLOCK_ID 0'
  [../]
  [./stress_xx]
    type = ElementAverageValue
    variable = stress_xx
   # block = 'ANY_BLOCK_ID 0'
  [../]
[]


[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
   type = Transient
   solve_type = 'PJFNK'
 
   petsc_options_iname = '-pc_type -pc_factor_mat_solver_package -snes_type'
   petsc_options_value = 'lu            superlu_dist          vinewtonrsls'  
 
 #  petsc_options_iname = '-pc_type -pc_asm_overlap -sub_pc_type -ksp_type -ksp_gmres_restart'
 #  petsc_options_value = ' asm      2              lu            gmres     200'
   l_max_its = 20
   nl_max_its = 10
   nl_rel_tol = 1.0e-8
   nl_abs_tol = 1.0e-9
 
   end_time = 40
 
   [./TimeStepper]
     type = IterationAdaptiveDT
     dt = 5e-4
     cutback_factor = 0.75
     growth_factor = 1.2
     optimal_iterations = 20
   [../]
[]

[Outputs]
  exodus = true
  csv = true
[]

sapitts Sep 16, 2024
Collaborator

Remember to use the eigenstrain_names parameter, with the appropriate eigenstrain name, in the ComputeMultipleCrystalPlasticityStress block. See the documentation previously linked for example input files

Emmanuel-324 Sep 16, 2024
Author

Yeah, I forgot to add that. This is the snippet of the block of interest

[elasticity_tensor_phase0]
     type = ComputeElasticityTensorCP
     C_ijkl = '2.906e5 1.87e5 1.87e5 2.906e5 1.87e5 2.906e5 1.142e5 1.142e5 1.142e5'
     fill_method = symmetric9
     euler_angle_1 = 0.0
     euler_angle_2 = 0.0
     euler_angle_3 = 0.0
     base_name = phase0
   []
   [stress_phase0]
     type = ComputeMultipleCrystalPlasticityStress
     crystal_plasticity_models = 'trial_xtalpl_phase0'
     eigenstrain_names = 'phase0_eigenstrain_0'
     tan_mod_type = exact
     rtol = 1e-08
     base_name = phase0
   []
   [trial_xtalpl_phase0]
     type = CrystalPlasticityKalidindiUpdate
     number_slip_systems = 12
     slip_sys_file_name = input_slip_sys.txt
     crystal_lattice_type = BCC
     resistance_tol = 0.01
     r = 1.4             
     h = 6000            
     t_sat = 598.5        
     gss_a = 1.5         
     ao = 0.001           
     xm = 0.017             
     gss_initial = 600
     base_name = phase0
   []
   [phase0_eigenstrain_0]
     type = ComputeCrystalPlasticityThermalEigenstrain
     eigenstrain_name = phase0_eigenstrain_0
     deformation_gradient_name = phase0_deformation_gradient_0
     temperature = dummy_temperature
     thermal_expansion_coefficients = '0.0067 0.0286 0.0067'
     base_name = phase0
   []
   

   [./strain_phase0]
     type = ComputeFiniteStrain
     displacements = 'disp_x disp_y'
     base_name = phase0
   [../]
 
   [elasticity_tensor_phase1]
     type = ComputeElasticityTensorCP
     C_ijkl = '2.721e5 1.69e5 1.69e5 2.721e5 1.69e5 2.721e5 1.31e5 1.31e5 1.31e5'
     fill_method = symmetric9
     euler_angle_1 = 0.0
     euler_angle_2 = 0.0
     euler_angle_3 = 0.0
     base_name = phase1
   []
   [stress_phase1]
     type = ComputeMultipleCrystalPlasticityStress
     crystal_plasticity_models = 'trial_xtalpl_phase1'
     tan_mod_type = exact
     rtol = 1e-08
     base_name = phase1
   []
   [trial_xtalpl_phase1]
     type = CrystalPlasticityKalidindiUpdate
     number_slip_systems = 12
     slip_sys_file_name = input_slip_sys.txt
     crystal_lattice_type = FCC
     resistance_tol = 0.01
     r = 1.0             
     h = 6000            
     t_sat = 598.5        
     gss_a = 1.5         
     ao = 0.001           
     xm = 0.017             
     gss_initial = 465.5 
     base_name = phase1
   []
   [./strain_phase1]
     type = ComputeFiniteStrain
     displacements = 'disp_x disp_y'
     base_name = phase1
   [../]

Emmanuel-324 Sep 16, 2024
Author

Unfortunately I am getting a lot of constitutive failure

Resetting modules to system default. Reseting $MODULEPATH back to system default. All extra directories will be removed from $MODULEPATH.
Run 'mamba init' to be able to run mamba activate/deactivate
and start a new shell session. Or use conda to activate/deactivate.


�[33m
*** Warning, This code is deprecated and will be removed in future versions:
/home/emmanuel324/projects/belson/My_Project/Restart_KKS_mod/9-a16-24/Test2/crystal_plasticity_2000.i:204.3: The 'Kernels/TensorMechanics' syntax is deprecated. Please use 'Physics/SolidMechanics/QuasiStatic' instead.�[39m


�[36m
*** Info ***
SolidMechanics Action: selecting 'total small strain' formulation. Use `incremental = true` to select 'incremental small strain' instead.�[39m
Currently Setting Up
  Initializing
    Finished Initializing Equation Systems                                               [�[33m  0.89 s�[39m] [�[33m  493 MB�[39m]
  Finished Initializing                                                                  [�[33m  0.94 s�[39m] [�[33m  493 MB�[39m]
Finished Setting Up                                                                      [�[33m  1.68 s�[39m] [�[33m  493 MB�[39m]
Framework Information:
MOOSE Version:           git commit 4505a3e on 2024-06-10
LibMesh Version:         1c5369aa05eaae8ebeb7930845fce3b1d59cbb37
PETSc Version:           3.20.3
SLEPc Version:           3.20.1
Current Time:            Mon Sep 16 19:02:16 2024
Executable Timestamp:    Wed Sep 11 14:18:50 2024

Parallelism:
  Num Processors:          2
  Num Threads:             1

Mesh: 
  Parallel Type:           replicated
  Mesh Dimension:          2
  Spatial Dimension:       2
  Nodes:                   
    Total:                 74664
    Local:                 39574
    Min/Max/Avg:           35090/39574/37332
  Elems:                   
    Total:                 70134
    Local:                 35067
    Min/Max/Avg:           35067/35067/35067
  Num Subdomains:          1
  Num Partitions:          2
  Partitioner:             metis

Nonlinear System:
  Num DOFs:                298656
  Num Local DOFs:          158296
  Variables:               { "disp_x" "disp_y" "eta0" "eta1" } 
  Finite Element Types:    "LAGRANGE" 
  Approximation Orders:    "FIRST" 

Auxiliary System:
  Num DOFs:                210402
  Num Local DOFs:          105201
  Variables:               { "vonmises" "stress_xx" "dummy_temperature" } 
  Finite Element Types:    "MONOMIAL" 
  Approximation Orders:    "CONSTANT" 

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

Currently Executing
  Performing Initial Setup
    Setting Up Materials
      Finished Computing Initial Material Values                                         [�[33m  0.75 s�[39m] [�[33m  998 MB�[39m]
    Finished Setting Up Materials                                                        [�[33m  0.75 s�[39m] [�[33m  998 MB�[39m]
  Finished Performing Initial Setup                                                      [�[33m  1.79 s�[39m] [�[33m 1011 MB�[39m]

Time Step 0, time = 0

Postprocessor Values:
+----------------+----------------+----------------+
| time           | stress_xx      | vonmises       |
+----------------+----------------+----------------+
|   0.000000e+00 |   0.000000e+00 |   0.000000e+00 |
+----------------+----------------+----------------+

Warning:  This MeshOutput subclass only supports meshes which are contiguously renumbered!

Time Step 1, time = 0.0005, dt = 0.0005

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 0.00025, dt = 0.00025

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 0.000125, dt = 0.000125

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 6.25e-05, dt = 6.25e-05

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 3.125e-05, dt = 3.125e-05

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.5625e-05, dt = 1.5625e-05

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 7.8125e-06, dt = 7.8125e-06

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 3.90625e-06, dt = 3.90625e-06

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.95313e-06, dt = 1.95313e-06

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 9.76563e-07, dt = 9.76563e-07

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 4.88281e-07, dt = 4.88281e-07

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 2.44141e-07, dt = 2.44141e-07

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.2207e-07, dt = 1.2207e-07

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 6.10352e-08, dt = 6.10352e-08

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 3.05176e-08, dt = 3.05176e-08

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.52588e-08, dt = 1.52588e-08

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 7.62939e-09, dt = 7.62939e-09

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 3.8147e-09, dt = 3.8147e-09

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.90735e-09, dt = 1.90735e-09

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 9.53674e-10, dt = 9.53674e-10

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 4.76837e-10, dt = 4.76837e-10

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 2.38419e-10, dt = 2.38419e-10

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.19209e-10, dt = 1.19209e-10

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 5.96046e-11, dt = 5.96046e-11

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

Nonlinear solve did not converge due to DIVERGED_FUNCTION_DOMAIN iterations 0
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 2.98023e-11, dt = 2.98023e-11

    Computing Residual.                                                                  [�[33m 15.87 s�[39m] [�[33m 1107 MB�[39m]
 0 Nonlinear |R| = �[32m8.895229e+05�[39m
    Computing Jacobian.                                                                  [�[33m 13.15 s�[39m] [�[33m 1149 MB�[39m]
      0 Linear |R| = �[32m8.895229e+05�[39m

    Computing Residual...                                                                [�[33m 27.83 s�[39m] [�[33m 2114 MB�[39m]
      1 Linear |R| = �[32m1.664772e+05�[39m
    Computing Residual...                                                                [�[33m 28.69 s�[39m] [�[33m 2114 MB�[39m]
      2 Linear |R| = �[32m3.334291e+04�[39m
    Computing Residual...                                                                [�[33m 28.28 s�[39m] [�[33m 2114 MB�[39m]
      3 Linear |R| = �[32m1.215332e+04�[39m
    Computing Residual...                                                                [�[33m 27.34 s�[39m] [�[33m 2114 MB�[39m]
      4 Linear |R| = �[32m3.593148e+03�[39m
    Computing Residual...                                                                [�[33m 28.99 s�[39m] [�[33m 2114 MB�[39m]
      5 Linear |R| = �[32m9.162615e+02�[39m
    Computing Residual...                                                                [�[33m 28.37 s�[39m] [�[33m 2114 MB�[39m]
      6 Linear |R| = �[32m2.526824e+02�[39m
    Computing Residual...                                                                [�[33m 29.00 s�[39m] [�[33m 2114 MB�[39m]
      7 Linear |R| = �[32m6.709273e+01�[39m
    Computing Residual...                                                                [�[33m 28.66 s�[39m] [�[33m 2114 MB�[39m]
      8 Linear |R| = �[32m1.763300e+01�[39m
    Computing Residual...                                                                [�[33m 28.01 s�[39m] [�[33m 2115 MB�[39m]
      9 Linear |R| = �[32m4.639385e+00�[39m
    Computing Residual...                                                                [�[33m 28.16 s�[39m] [�[33m 2115 MB�[39m]

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


A MooseException was raised during FEProblemBase::computeResidualTags
ComputeMultipleCrystalPlasticityStress: Constitutive failure
To recover, the solution will fail and then be re-attempted with a reduced time step.

 1 Nonlinear |R| = �[33m8.895229e+05�[39m
Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
  Finished Solving                                                                       [�[33m321.86 s�[39m] [�[33m 1637 MB�[39m]
�[31m Solve Did NOT Converge!�[39m
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 1.49012e-11, dt = 1.49012e-11

    Computing Residual.                                                                  [�[33m 16.99 s�[39m] [�[33m 1637 MB�[39m]
 0 Nonlinear |R| = �[32m8.892258e+05�[39m
    Computing Jacobian.                                                                  [�[33m 13.00 s�[39m] [�[33m 1637 MB�[39m]
      0 Linear |R| = �[32m8.892258e+05�[39m

    Computing Residual.

zuckarm Dec 23, 2024

hi @Emmanuel-324, were you able to resolve this issue? i am facing some similar errors while solving a thermomechanical problem.

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

EigenStrain contribution #28558

{{title}}

{{editor}}'s edit

{{editor}}'s edit

Replies: 3 comments 22 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

EigenStrain contribution #28558

Emmanuel-324 Sep 6, 2024

Replies: 3 comments · 22 replies

GiudGiud Sep 8, 2024 Collaborator

Emmanuel-324 Sep 8, 2024 Author

GiudGiud Sep 8, 2024 Collaborator

Emmanuel-324 Sep 9, 2024 Author

GiudGiud Sep 9, 2024 Collaborator

laagesen Sep 9, 2024

Emmanuel-324 Sep 9, 2024 Author

laagesen Sep 9, 2024

sapitts Sep 9, 2024 Collaborator

Emmanuel-324 Sep 16, 2024 Author

sapitts Sep 16, 2024 Collaborator

Emmanuel-324 Sep 16, 2024 Author

Emmanuel-324 Sep 16, 2024 Author

zuckarm Dec 23, 2024

Emmanuel-324
Sep 6, 2024

Replies: 3 comments 22 replies

GiudGiud
Sep 8, 2024
Collaborator

Emmanuel-324 Sep 8, 2024
Author

GiudGiud Sep 8, 2024
Collaborator

Emmanuel-324 Sep 9, 2024
Author

GiudGiud Sep 9, 2024
Collaborator

laagesen
Sep 9, 2024

Emmanuel-324 Sep 9, 2024
Author

sapitts
Sep 9, 2024
Collaborator

Emmanuel-324 Sep 16, 2024
Author

sapitts Sep 16, 2024
Collaborator

Emmanuel-324 Sep 16, 2024
Author

Emmanuel-324 Sep 16, 2024
Author