HomogenizedTotalLagrangianStressDivergenceA.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 "TotalLagrangianStressDivergenceS.h"

// Helpers common to the whole homogenization system
namespace HomogenizationA
{
const MultiMooseEnum constraintType("strain stress none");
enum class ConstraintType
{
  Strain,
  Stress,
  None
};
typedef std::map<std::pair<unsigned int, unsigned int>, std::pair<ConstraintType, const Function *>>
    ConstraintMap;
}

class HomogenizedTotalLagrangianStressDivergenceA : public TotalLagrangianStressDivergenceS
{
public:
  static InputParameters validParams();
  HomogenizedTotalLagrangianStressDivergenceA(const InputParameters & parameters);

protected:
  // Add overrides to base class contributions to only happen for _beta==0, to happen only once
  virtual Real computeQpResidual() override;
  virtual Real computeQpJacobianDisplacement(unsigned int alpha, unsigned int beta) override;

  virtual void computeScalarResidual() override;

  virtual void computeScalarJacobian() override;

  virtual void computeScalarOffDiagJacobian(const unsigned int jvar_num) override;

  virtual Real computeScalarQpOffDiagJacobian(const unsigned int jvar_num) override;

  virtual void computeOffDiagJacobianScalarLocal(const unsigned int svar_num) override;

  virtual Real computeQpOffDiagJacobianScalar(const unsigned int svar_num) override;

  virtual void computeScalarOffDiagJacobianScalar(const unsigned int svar_num) override;

protected:
  const unsigned int _beta;

  const MooseVariableScalar * const _kappao_var_ptr;

  const unsigned int _kappao_var;

  const unsigned int _ko_order;

  const VariableValue & _kappa_other;

  HomogenizationA::ConstraintMap _cmap;

  HomogenizationA::ConstraintType _ctype = HomogenizationA::ConstraintType::None;

  unsigned int _m;
  unsigned int _n;
  unsigned int _a;
  unsigned int _b;
};