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