MultiPlasticityRawComponentAssembler.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "TensorMechanicsPlasticModel.h"
#include "UserObjectInterface.h"
 
class MultiPlasticityRawComponentAssembler;
 
template <>
InputParameters validParams<MultiPlasticityRawComponentAssembler>();
 
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 TensorMechanicsPlasticModel *> _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);
};