TotalLagrangianStressDivergenceBase.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 "LagrangianStressDivergenceBase.h"
#include "GradientOperator.h"

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

protected:
  virtual RankTwoTensor gradTest(unsigned int component) override;
  virtual RankTwoTensor gradTrial(unsigned int component) override;
  virtual void precalculateJacobianDisplacement(unsigned int component) override;
  virtual void precalculateResidual() override;
  virtual void precalculateJacobian() override;
  virtual void precalculateOffDiagJacobian(unsigned int jvar) override;
  virtual Real computeQpResidual() override;
  virtual Real computeQpJacobianDisplacement(unsigned int alpha, unsigned int beta) override;
  virtual Real computeQpJacobianTemperature(unsigned int cvar) override;
  virtual Real computeQpJacobianOutOfPlaneStrain() override;

  void computeAverageGradientSpatialTest();

  void computeAverageGradientSpatialPhi(unsigned int beta);

  void computeAvgTestPhiCross(unsigned int beta);

  Real gradXTestComponent(unsigned int component) const;

  Real gradXPhiComponent(unsigned int component) const;

  const MaterialProperty<RankTwoTensor> & _pk1;

  const MaterialProperty<RankFourTensor> & _dpk1;

  const MaterialProperty<RankFourTensor> & _dpk1_d_grad_u;

  const MaterialProperty<RankFourTensor> & _dpk1_bypass_fbar;

protected:
  virtual RankTwoTensor gradTrialUnstabilized(unsigned int component);

private:
  void populateLocalPK1Cache(unsigned int beta);

  void cachePK1ContractionFUst();

  std::vector<std::vector<RankTwoTensor>> _dpk1_grad_trial_cache;

  std::vector<std::vector<RankTwoTensor>> _grad_trial_cache;

  std::vector<std::vector<RankTwoTensor>> _delta_PK1_NL_cache;

  std::vector<std::vector<RankTwoTensor>> _dpk1_total_cache;

  std::vector<Real> _pk1_ddot_F_ust;

  std::vector<std::vector<Real>> _dA_dU_cache;
};