HingeModel.h

Go to the documentation of this file.

 
#pragma once
 
#include <jem/base/Class.h>
#include <jem/util/ArrayBuffer.h>
#include <jem/util/Properties.h>
#include <jem/util/StringUtils.h>
 
#include <jive/Array.h>
#include <jive/algebra/MatrixBuilder.h>
#include <jive/fem/ElementGroup.h>
#include <jive/fem/ElementSet.h>
#include <jive/fem/NodeSet.h>
#include <jive/fem/XElementSet.h>
#include <jive/fem/XNodeSet.h>
#include <jive/implict/Names.h>
#include <jive/model/Actions.h>
#include <jive/model/Model.h>
#include <jive/model/ModelFactory.h>
#include <jive/model/StateVector.h>
#include <jive/solver/Solver.h>
#include <jive/util/Assignable.h>
#include <jive/util/Constraints.h>
#include <jive/util/DofSpace.h>
#include <jive/util/FuncUtils.h>
 
#include "models/SpecialCosseratRodModel.h"
#include "utils/helpers.h"
#include "utils/testing.h"
 
using jem::newInstance;
using jem::Ref;
using jem::String;
using jem::numeric::Function;
using jem::util::Properties;
 
using jive::BoolMatrix;
using jive::StringVector;
using jive::algebra::MatrixBuilder;
using jive::fem::ElementGroup;
using jive::fem::ElementSet;
using jive::fem::NodeSet;
using jive::fem::XElementSet;
using jive::fem::XNodeSet;
using jive::model::Model;
using jive::util::Assignable;
using jive::util::Constraints;
using jive::util::DofSpace;
using jive::util::FuncUtils;
 
typedef jem::util::ArrayBuffer<idx_t> IdxBuffer;
 
class HingeModel : public Model
{
public:
  static const char *TYPE_NAME;  
  static const char *YIELD_PROP; 
 
  JEM_DECLARE_CLASS(HingeModel, Model);
 
  explicit HingeModel
 
      (const String &name,
       const Properties &conf,
       const Properties &props,
       const Properties &globdat);
 
  virtual bool takeAction
 
      (const String &action,
       const Properties &params,
       const Properties &globdat) override;
 
  static Ref<Model> makeNew
 
      (const String &name,
       const Properties &conf,
       const Properties &props,
       const Properties &globdat);
 
  static void declare();
 
private:
  ElementGroup createHinges_(const String &elementName, const Properties &globdat);
 
  void updForces_(const Vector &fint);
 
  bool evalPlastic_(const Vector &disp);
 
  void init_(const Properties &globdat);
 
  void getCons_();
 
private:
  Assignable<ElementGroup> egroup_; 
  Assignable<NodeSet> nodes_;       
  Assignable<ElementSet> elems_;    
  Ref<DofSpace> dofs_;              
  Ref<Constraints> cons_;           
 
  IdxVector jtypes_;    
  StringVector jnames_; 
  String elName_;       
 
  double prec_;   
  idx_t maxIter_; 
  idx_t iiter_;   
 
  Ref<Material> material_;  
  Ref<Function> yieldCond_; 
  Vector ell_;              
  Matrix intForces_;        
  Matrix intForcesOld_;     
  Matrix plasticDisp_;      
};