LagrangianStressDivergenceBase.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 "KernelScalarBase.h"
#include "DerivativeMaterialInterface.h"
#include "JvarMapInterface.h"
#include "StabilizationUtils.h"
#include "RankFourTensorForward.h"
#include "GuaranteeConsumer.h"
 
class LagrangianStressDivergenceBase
  : public JvarMapKernelInterface<DerivativeMaterialInterface<KernelScalarBase>>,
    public GuaranteeConsumer
{
public:
  static InputParameters validParams();
  LagrangianStressDivergenceBase(const InputParameters & parameters);
 
  enum class FBarMode
  {
    Total,
    Incremental
  };
 
protected:
  // Helper function to return the test function gradient which may depend on kinematics and
  // stabilization
  virtual RankTwoTensor gradTest(unsigned int component) = 0;
 
  // Helper function to return the trial function gradient which may depend on kinematics and
  // stabilization
  virtual RankTwoTensor gradTrial(unsigned int component) = 0;
 
  virtual void initialSetup() override;
 
  virtual void precalculateJacobian() override;
  virtual void precalculateOffDiagJacobian(unsigned int jvar) override;
 
  virtual void precalculateJacobianDisplacement(unsigned int component) = 0;
 
  virtual Real computeQpJacobian() override;
  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
 
  // Derivatives of the residual w.r.t. the displacement dofs
  virtual Real computeQpJacobianDisplacement(unsigned int alpha, unsigned int beta) = 0;
 
  // Derivatives of the residual w.r.t. the temperature dofs through eigenstrain
  virtual Real computeQpJacobianTemperature(unsigned int cvar) = 0;
 
  // Derivatives of the residual w.r.t. the out-of-plane strain
  virtual Real computeQpJacobianOutOfPlaneStrain() = 0;
 
  RankTwoTensor deltaPK1NonLocalFBar(const RankTwoTensor & delta_F_avg) const;
 
protected:
  bool _large_kinematics;
 
  const bool _stabilize_strain;
 
  const FBarMode _F_bar_mode;
 
  const std::string _base_name;
 
  const unsigned int _alpha;
 
  const unsigned int _ndisp;
 
  std::vector<unsigned int> _disp_nums;
 
  // Averaged trial function gradients for each displacement component
  // i.e. _avg_grad_trial[a][j] returns the average gradient of trial function associated with
  // node j with respect to displacement component a.
  std::vector<std::vector<RankTwoTensor>> _avg_grad_trial;
 
  std::vector<Real> _avg_grad_spatial_test;
 
  std::vector<std::vector<Real>> _avg_grad_spatial_phi;
 
  std::vector<std::vector<std::vector<std::vector<Real>>>> _avg_test_phi_cross;
 
  const MaterialProperty<RankTwoTensor> & _F_ust;
 
  const MaterialProperty<RankTwoTensor> & _F_ust_old;
 
  const MaterialProperty<RankTwoTensor> & _F_avg;
 
  const MaterialProperty<RankTwoTensor> & _f_inv;
 
  const MaterialProperty<RankTwoTensor> & _F_inv;
 
  const MaterialProperty<RankTwoTensor> * _F_ust_inv = nullptr;
  const MaterialProperty<Real> * _F_ust_det = nullptr;
 
  const MaterialProperty<RankTwoTensor> & _F;
 
  const MaterialProperty<RankTwoTensor> & _F_actual;
 
  const MaterialProperty<RankFourTensor> & _d_deformation_gradient_increment_d_F;
 
  const MaterialProperty<RankFourTensor> & _d_F_d_grad_u;
 
  const MaterialProperty<RankFourTensor> & _d_F_stab_d_F_ust;
  const MaterialProperty<RankFourTensor> & _d_F_stab_d_F_avg;
 
  const MooseVariable * _temperature;
 
  const MooseVariable * _out_of_plane_strain;
 
  std::vector<std::vector<const MaterialProperty<RankTwoTensor> *>> _deigenstrain_dargs;
 
  const MaterialProperty<RankFourTensor> * _dcauchy_stress_d_eigenstrain = nullptr;
 
  const MaterialProperty<RankFourTensor> * _d_nl_fbar = nullptr;
};