HomogenizationConstraint.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 "ElementUserObject.h"

// Helpers common to the whole homogenization system
namespace Homogenization
{
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 HomogenizationConstraint : public ElementUserObject
{
public:
  static InputParameters validParams();

  HomogenizationConstraint(const InputParameters & parameters);

  virtual void initialize() override;
  virtual void execute() override;
  virtual void threadJoin(const UserObject & y) override;
  virtual void finalize() override;

  virtual const RankTwoTensor & getResidual() const { return _residual; }
  virtual const RankFourTensor & getJacobian() const { return _jacobian; }
  virtual const Homogenization::ConstraintMap & getConstraintMap() const { return _cmap; }

protected:
  virtual RankTwoTensor computeResidual();
  virtual RankFourTensor computeJacobian();

protected:
  const bool _large_kinematics;

  const std::string _base_name;

  const MaterialProperty<RankTwoTensor> & _F;
  const MaterialProperty<RankTwoTensor> & _pk1_stress;
  const MaterialProperty<RankFourTensor> & _pk1_jacobian;

  Homogenization::ConstraintMap _cmap;

private:
  unsigned int _qp;

  RankTwoTensor _residual;
  RankFourTensor _jacobian;
};