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;
};