ComputeConcentrationDependentElasticityTensor.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 "ComputeConcentrationDependentElasticityTensor.h"
#include "RotationTensor.h"

registerMooseObject("SolidMechanicsApp", ComputeConcentrationDependentElasticityTensor);

InputParameters
ComputeConcentrationDependentElasticityTensor::validParams()
{
  InputParameters params = ComputeRotatedElasticityTensorBase::validParams();
  params.addClassDescription("Compute concentration dependent elasticity tensor.");
  params.addRequiredParam<std::vector<Real>>("C0_ijkl",
                                             "Stiffness tensor for zero concentration phase");
  params.addRequiredParam<std::vector<Real>>("C1_ijkl",
                                             "Stiffness tensor for phase having concentration 1.0");
  params.addParam<MooseEnum>(
      "fill_method0", RankFourTensor::fillMethodEnum() = "symmetric9", "The fill method");
  params.addParam<MooseEnum>(
      "fill_method1", RankFourTensor::fillMethodEnum() = "symmetric9", "The fill method");
  params.addRequiredCoupledVar("c", "Concentration");
  return params;
}

ComputeConcentrationDependentElasticityTensor::ComputeConcentrationDependentElasticityTensor(
    const InputParameters & parameters)
  : ComputeRotatedElasticityTensorBase(parameters),
    _Cijkl0(getParam<std::vector<Real>>("C0_ijkl"),
            (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method0")),
    _Cijkl1(getParam<std::vector<Real>>("C1_ijkl"),
            (RankFourTensor::FillMethod)(int)getParam<MooseEnum>("fill_method1")),
    _c(coupledValue("c")),
    _c_name(coupledName("c", 0)),
    _delasticity_tensor_dc(isCoupledConstant("c") ? nullptr
                                                  : &declarePropertyDerivative<RankFourTensor>(
                                                        _elasticity_tensor_name, _c_name))
{
  // Define a rotation according to Euler angle parameters
  RotationTensor R(_Euler_angles); // R type: RealTensorValue

  // Rotate tensors
  _Cijkl0.rotate(R);
  _Cijkl1.rotate(R);
}

void
ComputeConcentrationDependentElasticityTensor::computeQpElasticityTensor()
{
  // Assign elasticity tensor at a given quad point
  _elasticity_tensor[_qp] = _Cijkl0 + (_Cijkl1 - _Cijkl0) * _c[_qp];
  // Define derivative of elasticity tensor with respect to concentration.
  if (_delasticity_tensor_dc)
    (*_delasticity_tensor_dc)[_qp] = (_Cijkl1 - _Cijkl0);
}