Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://www.mooseframework.org
       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 "ComputeEigenstrain.h"
      11             : 
      12             : registerMooseObject("TensorMechanicsApp", ComputeEigenstrain);
      13             : registerMooseObject("TensorMechanicsApp", ADComputeEigenstrain);
      14             : 
      15             : template <bool is_ad>
      16             : InputParameters
      17          72 : ComputeEigenstrainTempl<is_ad>::validParams()
      18             : {
      19          72 :   InputParameters params = ComputeEigenstrainBase::validParams();
      20          72 :   params.addClassDescription("Computes a constant Eigenstrain");
      21         144 :   params.addRequiredParam<std::vector<Real>>(
      22             :       "eigen_base", "Vector of values defining the constant base tensor for the Eigenstrain");
      23         144 :   params.addParam<MaterialPropertyName>(
      24         144 :       "prefactor", 1.0, "Name of material property defining the variable dependence");
      25          72 :   return params;
      26           0 : }
      27             : 
      28             : template <bool is_ad>
      29          54 : ComputeEigenstrainTempl<is_ad>::ComputeEigenstrainTempl(const InputParameters & parameters)
      30             :   : ComputeEigenstrainBaseTempl<is_ad>(parameters),
      31          54 :     _prefactor(this->template getGenericMaterialProperty<Real, is_ad>("prefactor"))
      32             : {
      33         108 :   _eigen_base_tensor.fillFromInputVector(this->template getParam<std::vector<Real>>("eigen_base"));
      34          54 : }
      35             : 
      36             : template <bool is_ad>
      37             : void
      38       17888 : ComputeEigenstrainTempl<is_ad>::computeQpEigenstrain()
      39             : {
      40             :   // Define Eigenstrain
      41       17888 :   _eigenstrain[_qp] = _eigen_base_tensor * _prefactor[_qp];
      42       17888 : }