No stateful props used (idaholab#25840)

loganharbour · Feb 20, 2024 · dccd45e · dccd45e
1 parent 9d500a0
commit dccd45e
Show file tree

Hide file tree

Showing 4 changed files with 2 additions and 26 deletions.
diff --git a/...chanics/include/materials/crystal_plasticity/ComputeCrystalPlasticityThermalEigenstrain.h b/...chanics/include/materials/crystal_plasticity/ComputeCrystalPlasticityThermalEigenstrain.h
@@ -27,9 +27,6 @@ class ComputeCrystalPlasticityThermalEigenstrain
 
   ComputeCrystalPlasticityThermalEigenstrain(const InputParameters & parameters);
 
-  /// We need to set initial values for lattice thermal expansion coefficients
-  virtual void initQpStatefulProperties() override;
-
 protected:
   ///Compute the deformation gradient due to thermal expansion
   virtual void computeQpDeformationGradient() override;

diff --git a/modules/tensor_mechanics/include/materials/crystal_plasticity/ComputeElasticityTensorCP.h b/modules/tensor_mechanics/include/materials/crystal_plasticity/ComputeElasticityTensorCP.h
@@ -34,7 +34,6 @@ class ComputeElasticityTensorCP : public ComputeElasticityTensor
    * Defines the constant rotation matrix from the user specified
    * Bunge Euler Angles or user-supplied rotation matrix.
    */
-  virtual void initQpStatefulProperties() override;
   virtual void computeQpElasticityTensor() override;
 
   virtual void assignEulerAngles();

diff --git a/...r_mechanics/src/materials/crystal_plasticity/ComputeCrystalPlasticityThermalEigenstrain.C b/...r_mechanics/src/materials/crystal_plasticity/ComputeCrystalPlasticityThermalEigenstrain.C
@@ -49,17 +49,12 @@ ComputeCrystalPlasticityThermalEigenstrain::ComputeCrystalPlasticityThermalEigen
 }
 
 void
-ComputeCrystalPlasticityThermalEigenstrain::initQpStatefulProperties()
+ComputeCrystalPlasticityThermalEigenstrain::computeQpDeformationGradient()
 {
-  ComputeCrystalPlasticityEigenstrainBase::initQpStatefulProperties();
-  // rotate the thermal deforamtion gradient for crystals based on Euler angles
+  // rotate the thermal deformation gradient for crystals based on Euler angles
   _lattice_thermal_expansion_coefficients[_qp] =
       _thermal_expansion_coefficients.rotated(_crysrot[_qp]);
-}
 
-void
-ComputeCrystalPlasticityThermalEigenstrain::computeQpDeformationGradient()
-{
   // compute the deformation gradient due to thermal expansion
   Real dtheta = (_temperature[_qp] - _temperature_old[_qp]) * _substep_dt / _dt;
   RankTwoTensor residual_equivalent_thermal_expansion_increment =

diff --git a/modules/tensor_mechanics/src/materials/crystal_plasticity/ComputeElasticityTensorCP.C b/modules/tensor_mechanics/src/materials/crystal_plasticity/ComputeElasticityTensorCP.C
@@ -22,9 +22,6 @@ ComputeElasticityTensorCP::validParams()
                                   "GeneralUserObject to read element "
                                   "specific property values from file");
 
-  // WIP figure out why
-  params.set<bool>("_force_stateful_init") = true;
-
   return params;
 }
 
@@ -77,18 +74,6 @@ ComputeElasticityTensorCP::assignEulerAngles()
   _R.update(_Euler_angles_mat_prop[_qp]);
 }
 
-void
-ComputeElasticityTensorCP::initQpStatefulProperties()
-{
-  if (!_user_provided_rotation_matrix)
-  {
-    assignEulerAngles();
-    _crysrot[_qp] = _R.transpose();
-  }
-  else
-    _crysrot[_qp] = _rotation_matrix.transpose();
-}
-
 void
 ComputeElasticityTensorCP::computeQpElasticityTensor()
 {