ComputeLagrangianObjectiveStress.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 "ComputeLagrangianStressCauchy.h"
#include "DerivativeMaterialPropertyNameInterface.h"

class ComputeLagrangianObjectiveStress : public ComputeLagrangianStressCauchy,
                                         public DerivativeMaterialPropertyNameInterface
{
public:
  static InputParameters validParams();
  ComputeLagrangianObjectiveStress(const InputParameters & parameters);

protected:
  virtual void initQpStatefulProperties() override;

  virtual void computeQpCauchyStress() override;

  //    This method must provide the _small_stress and _small_jacobian
  virtual void computeQpSmallStress() = 0;

  MaterialProperty<RankTwoTensor> & _small_stress;

  const MaterialProperty<RankTwoTensor> & _small_stress_old;

  MaterialProperty<RankFourTensor> & _small_jacobian;

  const MaterialProperty<RankTwoTensor> & _cauchy_stress_old;

  const MaterialProperty<RankTwoTensor> & _mechanical_strain;

  const MaterialProperty<RankTwoTensor> & _strain_increment;

  const MaterialProperty<RankTwoTensor> & _vorticity_increment;

  const MaterialProperty<RankTwoTensor> & _def_grad;

  const MaterialProperty<RankTwoTensor> & _def_grad_old;

  enum class ObjectiveRate
  {
    Truesdell,
    Jaumann,
    GreenNaghdi
  } _rate;

  const bool _polar_decomp;

  MaterialProperty<RankTwoTensor> * _rotation;

  const MaterialProperty<RankTwoTensor> * _rotation_old;

  MaterialProperty<RankFourTensor> * _d_rotation_d_def_grad;

  MaterialProperty<RankTwoTensor> * _stretch;

private:
  RankTwoTensor objectiveUpdateTruesdell(const RankTwoTensor & dS);

  RankTwoTensor objectiveUpdateJaumann(const RankTwoTensor & dS);

  RankTwoTensor objectiveUpdateGreenNaghdi(const RankTwoTensor & dS);

  std::tuple<RankTwoTensor, RankFourTensor> advectStress(const RankTwoTensor & S0,
                                                         const RankTwoTensor & dQ) const;

  RankFourTensor updateTensor(const RankTwoTensor & Q) const;

  RankFourTensor stressAdvectionDerivative(const RankTwoTensor & S) const;

  RankFourTensor cauchyJacobian(const RankFourTensor & Jinv, const RankFourTensor & U) const;

  void polarDecomposition();
};