ADComputeInstantaneousThermalExpansionFunctionEigenstrain.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 "ADComputeInstantaneousThermalExpansionFunctionEigenstrain.h"
 
#include "Function.h"
#include "RankTwoTensor.h"
 
registerADMooseObject("TensorMechanicsApp",
                      ADComputeInstantaneousThermalExpansionFunctionEigenstrain);
 
defineADLegacyParams(ADComputeInstantaneousThermalExpansionFunctionEigenstrain);
 
template <ComputeStage compute_stage>
InputParameters
ADComputeInstantaneousThermalExpansionFunctionEigenstrain<compute_stage>::validParams()
{
  InputParameters params = ADComputeThermalExpansionEigenstrainBase<compute_stage>::validParams();
  params.addClassDescription("Computes eigenstrain due to thermal expansion using a function that "
                             "describes the instantaneous thermal expansion as a function of "
                             "temperature");
  params.addRequiredParam<FunctionName>("thermal_expansion_function",
                                        "Function describing the instantaneous thermal expansion "
                                        "coefficient as a function of temperature");
  return params;
}
 
template <ComputeStage compute_stage>
ADComputeInstantaneousThermalExpansionFunctionEigenstrain<compute_stage>::
    ADComputeInstantaneousThermalExpansionFunctionEigenstrain(const InputParameters & parameters)
  : ADComputeThermalExpansionEigenstrainBase<compute_stage>(parameters),
    _temperature_old(coupledValueOld("temperature")),
    _thermal_expansion_function(getFunction("thermal_expansion_function")),
    _thermal_strain(
        declareADProperty<Real>("InstantaneousThermalExpansionFunction_thermal_strain")),
    _thermal_strain_old(
        getMaterialPropertyOld<Real>("InstantaneousThermalExpansionFunction_thermal_strain")),
    _step_one(declareRestartableData<bool>("step_one", true))
{
}
 
template <ComputeStage compute_stage>
void
ADComputeInstantaneousThermalExpansionFunctionEigenstrain<compute_stage>::initQpStatefulProperties()
{
  _thermal_strain[_qp] = 0;
}
 
template <ComputeStage compute_stage>
void
ADComputeInstantaneousThermalExpansionFunctionEigenstrain<compute_stage>::computeThermalStrain(
    ADReal & thermal_strain)
{
  if (_t_step > 1)
    _step_one = false;
 
  const Real & current_temp = MetaPhysicL::raw_value(_temperature[_qp]);
 
  const Real & old_thermal_strain = _thermal_strain_old[_qp];
 
  const Real & old_temp =
      (_step_one ? MetaPhysicL::raw_value(_stress_free_temperature[_qp]) : _temperature_old[_qp]);
  const Real delta_T = current_temp - old_temp;
 
  const Point p;
  const Real alpha_current_temp = _thermal_expansion_function.value(current_temp, p);
  const Real alpha_old_temp = _thermal_expansion_function.value(old_temp, p);
 
  thermal_strain = old_thermal_strain + delta_T * 0.5 * (alpha_current_temp + alpha_old_temp);
  _thermal_strain[_qp] = thermal_strain;
}