MultiPlasticityRawComponentAssembler.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

/* MOOSE - Multiphysics Object Oriented Simulation Environment  */
/*                                                              */
/*          All contents are licensed under LGPL V2.1           */
/*             See LICENSE for full restrictions                */
/****************************************************************/

#include "SolidMechanicsPlasticModel.h"
#include "UserObjectInterface.h"

class MultiPlasticityRawComponentAssembler : public UserObjectInterface
{
public:
  static InputParameters validParams();

  MultiPlasticityRawComponentAssembler(const MooseObject * moose_object);

  virtual ~MultiPlasticityRawComponentAssembler() {}

protected:
  const InputParameters & _params;

  unsigned int _num_models;

  unsigned int _num_surfaces;

  std::vector<std::vector<unsigned int>> _surfaces_given_model;

  MooseEnum _specialIC;

  std::vector<const SolidMechanicsPlasticModel *> _f;

  virtual void yieldFunction(const RankTwoTensor & stress,
                             const std::vector<Real> & intnl,
                             const std::vector<bool> & active,
                             std::vector<Real> & f);

  virtual void dyieldFunction_dstress(const RankTwoTensor & stress,
                                      const std::vector<Real> & intnl,
                                      const std::vector<bool> & active,
                                      std::vector<RankTwoTensor> & df_dstress);

  virtual void dyieldFunction_dintnl(const RankTwoTensor & stress,
                                     const std::vector<Real> & intnl,
                                     const std::vector<bool> & active,
                                     std::vector<Real> & df_dintnl);

  virtual void flowPotential(const RankTwoTensor & stress,
                             const std::vector<Real> & intnl,
                             const std::vector<bool> & active,
                             std::vector<RankTwoTensor> & r);

  virtual void dflowPotential_dstress(const RankTwoTensor & stress,
                                      const std::vector<Real> & intnl,
                                      const std::vector<bool> & active,
                                      std::vector<RankFourTensor> & dr_dstress);

  virtual void dflowPotential_dintnl(const RankTwoTensor & stress,
                                     const std::vector<Real> & intnl,
                                     const std::vector<bool> & active,
                                     std::vector<RankTwoTensor> & dr_dintnl);

  virtual void hardPotential(const RankTwoTensor & stress,
                             const std::vector<Real> & intnl,
                             const std::vector<bool> & active,
                             std::vector<Real> & h);

  virtual void dhardPotential_dstress(const RankTwoTensor & stress,
                                      const std::vector<Real> & intnl,
                                      const std::vector<bool> & active,
                                      std::vector<RankTwoTensor> & dh_dstress);

  virtual void dhardPotential_dintnl(const RankTwoTensor & stress,
                                     const std::vector<Real> & intnl,
                                     const std::vector<bool> & active,
                                     std::vector<Real> & dh_dintnl);

  virtual void buildActiveConstraints(const std::vector<Real> & f,
                                      const RankTwoTensor & stress,
                                      const std::vector<Real> & intnl,
                                      const RankFourTensor & Eijkl,
                                      std::vector<bool> & act);

  unsigned int modelNumber(unsigned int surface);

  bool anyActiveSurfaces(int model, const std::vector<bool> & active);

  void activeModelSurfaces(int model,
                           const std::vector<bool> & active,
                           std::vector<unsigned int> & active_surfaces_of_model);

  void activeSurfaces(int model,
                      const std::vector<bool> & active,
                      std::vector<unsigned int> & active_surfaces);

  bool returnMapAll(const RankTwoTensor & trial_stress,
                    const std::vector<Real> & intnl_old,
                    const RankFourTensor & E_ijkl,
                    Real ep_plastic_tolerance,
                    RankTwoTensor & stress,
                    std::vector<Real> & intnl,
                    std::vector<Real> & pm,
                    std::vector<Real> & cumulative_pm,
                    RankTwoTensor & delta_dp,
                    std::vector<Real> & yf,
                    unsigned & num_successful_plastic_returns,
                    unsigned & custom_model);

private:
  std::vector<unsigned int> _model_given_surface;

  std::vector<unsigned int> _model_surface_given_surface;

  void buildActiveConstraintsRock(const std::vector<Real> & f,
                                  const RankTwoTensor & stress,
                                  const std::vector<Real> & intnl,
                                  const RankFourTensor & Eijkl,
                                  std::vector<bool> & act);

  void buildActiveConstraintsJoint(const std::vector<Real> & f,
                                   const RankTwoTensor & stress,
                                   const std::vector<Real> & intnl,
                                   const RankFourTensor & Eijkl,
                                   std::vector<bool> & act);
};