ComputeElasticityTensorBase.C

Go to the documentation of this file.

//* 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 "ComputeElasticityTensorBase.h"
#include "Function.h"
 
defineLegacyParams(ComputeElasticityTensorBase);
 
InputParameters
ComputeElasticityTensorBase::validParams()
{
  InputParameters params = Material::validParams();
  params.addParam<FunctionName>(
      "elasticity_tensor_prefactor",
      "Optional function to use as a scalar prefactor on the elasticity tensor.");
  params.addParam<std::string>("base_name",
                               "Optional parameter that allows the user to define "
                               "multiple mechanics material systems on the same "
                               "block, i.e. for multiple phases");
  return params;
}
 
ComputeElasticityTensorBase::ComputeElasticityTensorBase(const InputParameters & parameters)
  : DerivativeMaterialInterface<Material>(parameters),
    GuaranteeProvider(this),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _elasticity_tensor_name(_base_name + "elasticity_tensor"),
    _elasticity_tensor(declareProperty<RankFourTensor>(_elasticity_tensor_name)),
    _effective_stiffness(declareProperty<Real>(_base_name + "effective_stiffness")),
    _prefactor_function(isParamValid("elasticity_tensor_prefactor")
                            ? &getFunction("elasticity_tensor_prefactor")
                            : nullptr)
{
}
 
void
ComputeElasticityTensorBase::computeQpProperties()
{
  _effective_stiffness[_qp] = 0; // Currently overriden by ComputeIsotropicElasticityTensor
  computeQpElasticityTensor();
 
  // Multiply by prefactor
  if (_prefactor_function)
  {
    _elasticity_tensor[_qp] *= _prefactor_function->value(_t, _q_point[_qp]);
    _effective_stiffness[_qp] *= std::sqrt(_prefactor_function->value(_t, _q_point[_qp]));
  }
}