Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "EigenDecompositionMaterial.h"
      11             : #include "RankTwoTensor.h"
      12             : 
      13             : registerMooseObject("SolidMechanicsApp", EigenDecompositionMaterial);
      14             : registerMooseObject("SolidMechanicsApp", ADEigenDecompositionMaterial);
      15             : 
      16             : template <bool is_ad>
      17             : InputParameters
      18          96 : EigenDecompositionMaterialTempl<is_ad>::validParams()
      19             : {
      20          96 :   InputParameters params = Material::validParams();
      21          96 :   params.addClassDescription("Emits material properties for the eigenvalues and eigenvectors of a "
      22             :                              "symmetric rank two tensor.");
      23         192 :   params.addRequiredParam<MaterialPropertyName>(
      24             :       "rank_two_tensor",
      25             :       "The name of the symmetric rank two tensor to used in eigen decomposition.");
      26         192 :   params.addParam<std::string>(
      27             :       "base_name",
      28             :       "Optional parameter to allow multiple tensors to be decomposed on the same block.");
      29          96 :   return params;
      30           0 : }
      31             : 
      32             : template <bool is_ad>
      33          72 : EigenDecompositionMaterialTempl<is_ad>::EigenDecompositionMaterialTempl(
      34             :     const InputParameters & parameters)
      35             :   : Material(parameters),
      36         132 :     _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
      37          72 :     _tensor(getGenericMaterialProperty<RankTwoTensor, is_ad>("rank_two_tensor")),
      38          72 :     _max_eigen_vector(
      39          72 :         declareGenericProperty<RealVectorValue, is_ad>(_base_name + "max_eigen_vector")),
      40          72 :     _mid_eigen_vector(
      41          72 :         declareGenericProperty<RealVectorValue, is_ad>(_base_name + "mid_eigen_vector")),
      42          72 :     _min_eigen_vector(
      43          72 :         declareGenericProperty<RealVectorValue, is_ad>(_base_name + "min_eigen_vector")),
      44         144 :     _max_eigen_value(declareGenericProperty<Real, is_ad>(_base_name + "max_eigen_value")),
      45         144 :     _mid_eigen_value(declareGenericProperty<Real, is_ad>(_base_name + "mid_eigen_value")),
      46         216 :     _min_eigen_value(declareGenericProperty<Real, is_ad>(_base_name + "min_eigen_value"))
      47             : {
      48             :   if (LIBMESH_DIM != 3)
      49             :     mooseError("EigenDecompositionMaterial is only defined for LIBMESH_DIM=3");
      50          72 : }
      51             : 
      52             : template <bool is_ad>
      53             : void
      54        1948 : EigenDecompositionMaterialTempl<is_ad>::computeQpProperties()
      55             : {
      56             : 
      57        1948 :   if (!_tensor[_qp].isSymmetric())
      58           2 :     mooseError("EigenDecompositionMaterial will only operate on symmetric rank two tensors.");
      59             : 
      60        1946 :   std::vector<GenericReal<is_ad>> eigval(3, 0.0);
      61        1946 :   GenericRankTwoTensor<is_ad> eigvec;
      62             : 
      63        1946 :   _tensor[_qp].symmetricEigenvaluesEigenvectors(eigval, eigvec);
      64             : 
      65        1946 :   _max_eigen_vector[_qp] = eigvec.column(2);
      66        1946 :   _mid_eigen_vector[_qp] = eigvec.column(1);
      67        1946 :   _min_eigen_vector[_qp] = eigvec.column(0);
      68             : 
      69        1946 :   _max_eigen_value[_qp] = eigval[2];
      70        1946 :   _mid_eigen_value[_qp] = eigval[1];
      71        1946 :   _min_eigen_value[_qp] = eigval[0];
      72        1946 : }