ElasticRodMaterial.h

Go to the documentation of this file.

 
#pragma once
#include "materials/Material.h"
#include <jem/numeric/algebra/MatmulChain.h>
 
// Forward declarations
class MaterialFactory;
 
using jem::idx_t;
using jive::Cubix;
using jive::Ref;
using jive::String;
 
class ElasticRodMaterial : public Material
{
 public:
  static const char *TYPE_NAME;
  static const char *STIFF_PREFIX;   
  static const char *YOUNGS_MODULUS; 
  static const char *SHEAR_MODULUS;  
  static const char *POISSON_RATIO;  
  static const char *AREA;           
  static const char *DENSITY;        
  static const char *AREA_MOMENT;    
  static const char *SHEAR_FACTOR;   
  static const char *POLAR_MOMENT;   
  static const char *CROSS_SECTION;  
  static const char *RADIUS;         
  static const char *SIDE_LENGTH;    
  static const char *N_ELEM;      
  static const char *EDGE_FACTOR; 
  static const char *EDGE_ELEMS;  
 
  JEM_DECLARE_CLASS(ElasticRodMaterial, Material);
 
  ElasticRodMaterial(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 Matrix getMaterialStiff() const override;
 
  virtual Matrix getMaterialStiff(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual void getMaterialStiff(const Matrix &stiff) const override;
 
  virtual void getMaterialStiff(const Matrix &stiff, const idx_t &ielem, const idx_t &ip) const override;
 
  virtual Matrix getMaterialMass() const override;
 
  virtual Matrix getMaterialMass(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual void getMaterialMass(const Matrix &mass) const override;
 
  virtual void getMaterialMass(const Matrix &mass, const idx_t &ielem, const idx_t &ip) const override;
 
  virtual Matrix getLumpedMass(const double l) const override;
 
  virtual Matrix getLumpedMass(const double l, const idx_t &ielem) const override;
 
  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 = false) override;
 
  virtual void getTable(XTable &strain_table, const Vector &weights, const String &name, const IdxVector &items) const override;
 
  virtual void applyDeform() override;
 
  virtual void rejectDeform() override;
 
  virtual double getDissipatedEnergy(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual double getPotentialEnergy(const idx_t &ielem, const idx_t &ip) const override;
 
  virtual double getHardeningPotential(const idx_t &ielem, const idx_t &ip) const override;
 
 protected:
  ~ElasticRodMaterial();
 
 protected:
  void calcMaterialStiff_();
  void calcMaterialMass_();
 
  double young_;      
  double shearMod_;   
  double shearParam_; 
 
  double edgeFact_; 
  idx_t edgeElems_; 
  idx_t nElem_;     
 
  double area_;        
  Vector areaMoment_;  
  double polarMoment_; 
 
  String crossSection_; 
  double radius_;       
  Vector sideLength_;   
 
  double density_; 
 
  Matrix materialK_; 
  Matrix materialM_; 
 
  Cubix oldStrains_;  
  Cubix currStrains_; 
  Matrix energyPot_;  
 
  jem::numeric::MatmulChain<double, 1> matvec_;
  Matrix stiffWork_; 
 
  String rodName_; 
};