StrainAdjustedDensity.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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 "StrainAdjustedDensity.h"

registerMooseObject("SolidMechanicsApp", StrainAdjustedDensity);
registerMooseObject("SolidMechanicsApp", ADStrainAdjustedDensity);

template <bool is_ad>
InputParameters
StrainAdjustedDensityTempl<is_ad>::validParams()
{
  InputParameters params = Material::validParams();

  params.addRequiredCoupledVar(
      "displacements",
      "The displacements appropriate for the simulation geometry and coordinate system. If "
      "displacements are not coupled, a different material class such as GenericMaterialProperty "
      "or ParsedMaterial should be used.");

  params.addParam<std::string>("base_name",
                               "Optional parameter that allows the user to define "
                               "multiple material systems on the same block, "
                               "e.g. for multiple phases");
  params.addRequiredParam<MaterialPropertyName>("strain_free_density",
                                                "Material property for strain-free density");
  params.addClassDescription("Creates density material property");

  return params;
}

template <bool is_ad>
StrainAdjustedDensityTempl<is_ad>::StrainAdjustedDensityTempl(const InputParameters & parameters)
  : Material(parameters),
    _coord_system(getBlockCoordSystem()),
    _disp_r(this->template coupledGenericValue<is_ad>("displacements", 0)),
    _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
    _strain_free_density(
        this->template getGenericMaterialProperty<Real, is_ad>("strain_free_density")),
    _grad_disp(this->template coupledGenericGradients<is_ad>("displacements")),
    _density(declareGenericProperty<Real, is_ad>(_base_name + "density"))
{
  if (getParam<bool>("use_displaced_mesh"))
    paramError("use_displaced_mesh",
               "StrainAdjustedDensity needs to act on an undisplaced mesh. Use of a displaced mesh "
               "leads to "
               "incorrect gradient values");

  // fill remaining components with zero
  _grad_disp.resize(3, &genericZeroGradient<is_ad>());
}

template <bool is_ad>
void
StrainAdjustedDensityTempl<is_ad>::initQpStatefulProperties()
{
  computeQpProperties();
}

template <bool is_ad>
void
StrainAdjustedDensityTempl<is_ad>::computeQpProperties()
{
  auto A = GenericRankTwoTensor<is_ad>::initializeFromRows(
      (*_grad_disp[0])[_qp], (*_grad_disp[1])[_qp], (*_grad_disp[2])[_qp]);
  A.addIa(1.0);

  switch (_coord_system)
  {
    case Moose::COORD_XYZ:
      break;

    case Moose::COORD_RZ:
      if (_q_point[_qp](0) != 0.0)
        A(2, 2) = _disp_r[_qp] / _q_point[_qp](0) + 1.0;
      break;

    case Moose::COORD_RSPHERICAL:
      if (_q_point[_qp](0) != 0.0)
        A(1, 1) = A(2, 2) = _disp_r[_qp] / _q_point[_qp](0) + 1.0;
      break;
  }

  _density[_qp] = _strain_free_density[_qp] / A.det();
}

template class StrainAdjustedDensityTempl<false>;
template class StrainAdjustedDensityTempl<true>;