ElastoPlasticRodMaterial.h

Go to the documentation of this file.

 
#pragma once
#include "materials/ElasticRodMaterial.h"
#include <jem/numeric/algebra/utilities.h>
#include <jive/util/FuncUtils.h>
 
// Forward declarations
class MaterialFactory;
 
using jem::idx_t;
using jem::numeric::Function;
using jive::Ref;
using jive::String;
using jive::util::FuncUtils;
 
class ElastoPlasticRodMaterial : public ElasticRodMaterial
{
 public:
  static const char *TYPE_NAME;
  static const char *YIELD_PROP;       
  static const char *YIELD_DERIV_PROP; 
  static const char *YIELD_TYPE_PROP;  
  static const char *YIELD_STRENGTH_PROP; 
  static const char *YIELD_EXPONENT_PROP; 
  static const char *ISO_HARD_PROP;    
  static const char *KIN_HARD_PROP;    
 
  enum YieldType
  {
    EXPRESSION,   
    SUPER_ELLIPSE 
  };
 
  JEM_DECLARE_CLASS(ElastoPlasticRodMaterial, ElasticRodMaterial);
 
  ElastoPlasticRodMaterial(const String &name,
                           const Properties &conf,
                           const Properties &props,
                           const Properties &globdat);
 
  static Ref<Material> makeNew(const String &name, const Properties &conf,
                               const Properties &props, const Properties &globdat);
 
  static void declare();
 
  virtual void configure(const Properties &props, const Properties &globdat) override;
 
  virtual void getConfig(const Properties &conf, const Properties &globdat) const override;
 
  virtual void getHardVals(const Vector &hardVals, const Vector &hardParams) const;
 
  virtual void getStress(const Vector &stress, const Vector &strain) override;
 
  virtual void getStress(const Vector &stress, const Vector &strain, const idx_t &ielem, const idx_t &ip, const bool inelastic = true) override;
 
  virtual void applyDeform() override;
 
  virtual void rejectDeform() override;
 
  virtual void getTable(XTable &strain_table, const Vector &weights, const String &name, const IdxVector &items) const override;
 
  virtual double getDissipatedEnergy(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual double getHardeningPotential(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual double getPotentialEnergy(const idx_t &ielem, const idx_t &ip) const override;
 
 protected:
  ~ElastoPlasticRodMaterial();
 
  double getYield_(const Vector &args) const;
 
  void getYieldGrad_(const Vector &grad, const Vector &args) const;
 
  double getSuperEllipse_(const Vector &args) const;
 
  void getSuperEllipseGrad_(const Vector &grad, const Vector &args) const;
 
 protected:
  Ref<Function> yieldCond_;         
  FuncUtils::FuncArray yieldDeriv_; 
  idx_t maxIter_;                   
  double precision_;                
  YieldType yieldType_;             
  bool isoHard_;                    
  bool kinHard_;                    
  Vector yieldStrength_;            
  Vector yieldExponent_;            
 
  idx_t dofCount_; 
  idx_t argCount_; 
 
  Matrix materialH_; 
 
  Cubix oldHardParams_;  
  Cubix currHardParams_; 
 
  Cubix oldPlastStrains_;  
  Cubix currPlastStrains_; 
 
  Cubix oldStresses_; 
 
  Vector plastStrain_;   
  Vector hardParams_;    
  Vector hardStress_;    
  Vector elasticStrain_; 
  Vector yieldArgs_;     
  Vector yieldGrad_;     
 
  Matrix currDeltaFlow_; 
  Matrix energyDiss_;    
  Matrix energyHardPot_; 
};