ComputeLagrangianStrainBase.h

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

#pragma once

#include "Material.h"
#include "RankFourTensorForward.h"
#include "RankTwoTensorForward.h"
#include "StabilizationUtils.h"
#include "GradientOperator.h"

template <class G>
class ComputeLagrangianStrainBase : public Material, public G
{
public:
  static InputParameters baseParams();
  static InputParameters validParams();
  ComputeLagrangianStrainBase(const InputParameters & parameters);
  virtual void initialSetup() override;

protected:
  virtual void initQpStatefulProperties() override;
  virtual void computeProperties() override;
  virtual void computeQpProperties() override;
  virtual void computeQpIncrementalStrains(const RankTwoTensor & dL);
  virtual void subtractQpEigenstrainIncrement(RankTwoTensor & strain);
  virtual void computeQpUnstabilizedDeformationGradient();
  virtual void computeDeformationGradient();

  // Displacements and displacement gradients
  const unsigned int _ndisp;
  std::vector<const VariableValue *> _disp;
  std::vector<const VariableGradient *> _grad_disp;

  const std::string _base_name;

  const bool _large_kinematics;

  const bool _stabilize_strain;

  // The eigenstrains
  std::vector<MaterialPropertyName> _eigenstrain_names;
  std::vector<const MaterialProperty<RankTwoTensor> *> _eigenstrains;
  std::vector<const MaterialProperty<RankTwoTensor> *> _eigenstrains_old;

  // The total strains
  MaterialProperty<RankTwoTensor> & _total_strain;
  const MaterialProperty<RankTwoTensor> & _total_strain_old;
  MaterialProperty<RankTwoTensor> & _mechanical_strain;
  const MaterialProperty<RankTwoTensor> & _mechanical_strain_old;

  MaterialProperty<RankTwoTensor> & _strain_increment;

  MaterialProperty<RankTwoTensor> & _spatial_velocity_increment;

  MaterialProperty<RankTwoTensor> & _vorticity_increment;

  MaterialProperty<RankTwoTensor> & _F_ust;

  // The average deformation gradient over the element for F-bar stabilization. Note that the
  // average deformation gradient is undefined if stabilization is not active.
  MaterialProperty<RankTwoTensor> & _F_avg;

  // The deformation gradient. If stabilization is active, this will be the stabilized deformation
  // gradient. Otherwise this will be equal to the unstabilized deformation gradient.
  MaterialProperty<RankTwoTensor> & _F;

  const MaterialProperty<RankTwoTensor> & _F_old;

  MaterialProperty<RankTwoTensor> & _F_inv;
  MaterialProperty<RankTwoTensor> & _f_inv;

  std::vector<MaterialPropertyName> _homogenization_gradient_names;

  std::vector<const MaterialProperty<RankTwoTensor> *> _homogenization_contributions;

  MaterialProperty<RankTwoTensor> & _rotation_increment;
};