ComputeCrystalPlasticityEigenstrainBase.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 "ComputeCrystalPlasticityEigenstrainBase.h"
#include "RankTwoTensor.h"

InputParameters
ComputeCrystalPlasticityEigenstrainBase::validParams()
{
  InputParameters params = ComputeEigenstrainBase::validParams();

  // The return stress increment classes are intended to be iterative materials, so must set compute
  // = false for all inheriting classes
  params.set<bool>("compute") = false;
  params.suppressParameter<bool>("compute");

  params.addRequiredParam<std::string>(
      "deformation_gradient_name",
      "Material property name for the deformation gradient tensor computed "
      "by this model.");
  return params;
}

ComputeCrystalPlasticityEigenstrainBase::ComputeCrystalPlasticityEigenstrainBase(
    const InputParameters & parameters)
  : ComputeEigenstrainBase(parameters),
    _deformation_gradient_name(_base_name + getParam<std::string>("deformation_gradient_name")),
    _deformation_gradient(declareProperty<RankTwoTensor>(_deformation_gradient_name)),
    _deformation_gradient_old(getMaterialPropertyOld<RankTwoTensor>(_deformation_gradient_name)),
    _crysrot(getMaterialProperty<RankTwoTensor>(
        "crysrot")) // defined in the elasticity tensor classes for crystal plasticity
{
}

void
ComputeCrystalPlasticityEigenstrainBase::initQpStatefulProperties()
{
  ComputeEigenstrainBase::initQpStatefulProperties();

  // Initialize deformation gradient to be identity
  _deformation_gradient[_qp].zero();
  _deformation_gradient[_qp].addIa(1.0);
}

void
ComputeCrystalPlasticityEigenstrainBase::computeQpEigenstrain()
{
  // updates the deformation gradient from the child class
  computeQpDeformationGradient();
  // use the updated deformation gradient for the eigenstrain calculation
  _eigenstrain[_qp] = 0.5 * (_deformation_gradient[_qp].transpose() * _deformation_gradient[_qp] -
                             RankTwoTensor::Identity());
}

void
ComputeCrystalPlasticityEigenstrainBase::setQp(const unsigned int & qp)
{
  _qp = qp;
}

void
ComputeCrystalPlasticityEigenstrainBase::setSubstepDt(const Real & substep_dt)
{
  _substep_dt = substep_dt;
}

const RankTwoTensor
ComputeCrystalPlasticityEigenstrainBase::getDeformationGradient() const
{
  return _deformation_gradient[_qp];
}

const RankTwoTensor
ComputeCrystalPlasticityEigenstrainBase::getDeformationGradientInverse() const
{
  return _deformation_gradient[_qp].inverse();
}