NeumannModel.h

Go to the documentation of this file.

 
#pragma once
 
#include <jem/base/Object.h>
#include <jem/io/Writer.h>
#include <jive/Array.h>
#include <jive/fem/NodeGroup.h>
#include <jive/model/Model.h>
#include <jive/model/ModelFactory.h>
#include <jive/util/Assignable.h>
#include <jive/util/XDofSpace.h>
 
using jem::Float;
using jem::idx_t;
using jem::IllegalInputException;
using jem::newInstance;
using jem::Ref;
using jem::SliceFrom;
using jem::String;
using jem::System;
using jem::io::Writer;
using jem::util::Properties;
using jive::IntVector;
using jive::StringVector;
using jive::Vector;
using jive::fem::NodeGroup;
using jive::fem::NodeSet;
using jive::model::Model;
using jive::util::Assignable;
using jive::util::Constraints;
using jive::util::DofSpace;
using jive::util::XDofSpace;
 
//-----------------------------------------------------------------------
//   class NeumannModel
//-----------------------------------------------------------------------
 
class NeumannModel : public Model
{
public:
  static const char *TYPE_NAME;      
  static const char *LOAD_INCR_PROP; 
  static const char *INIT_LOAD_PROP; 
  static const char *MIN_LOAD_PROP;  
  static const char *MAX_LOAD_PROP;  
  static const char *REDUCTION_PROP; 
  static const char *NODES_PROP;     
  static const char *DOF_PROP;       
  static const char *FACTORS_PROP;   
 
  JEM_DECLARE_CLASS(NeumannModel, Model);
 
  explicit NeumannModel(const String &name);
 
  virtual void configure(const Properties &props,
                         const Properties &globdat) override;
 
  virtual void getConfig(const Properties &conf,
                         const Properties &globdat) const override;
 
  virtual bool takeAction(const String &action,
                          const Properties &params,
                          const Properties &globdat) override;
 
  void setLoadIncr(double incr);
 
  inline double getLoadIncr() const;
 
  static Ref<Model> makeNew(const String &name,
                            const Properties &conf,
                            const Properties &props,
                            const Properties &globdat);
 
  static void declare();
 
protected:
  virtual ~NeumannModel();
 
private:
  void init_(const Properties &globdat);
 
  void getExtVector_(const Vector &fext,
                     const Properties &globdat) const;
 
  void getExtVector_(const Vector &fext,
                     const Properties &globdat,
                     const double scale) const;
 
  void advance_(const Properties &globdat);
 
  void commit_(const Properties &params,
               const Properties &globdat);
 
  void reduceStep_(const Properties &params,
                   const Properties &globdat);
 
  void increaseStep_(const Properties &params,
                     const Properties &globdat);
 
private:
  Ref<DofSpace> dofs_;        
  Assignable<NodeSet> nodes_; 
 
  idx_t ngroups_;   
  IntVector idofs_; 
 
  double loadScale0_; 
  double loadScale_;  
  double loadIncr_;   
 
  // Input specification for boundary conditions
  StringVector nodeGroups_; 
  StringVector dofTypes_;   
  Vector factors_;          
 
  // Configuration parameters
  double reduction_;   
  double loadIncr0_;   
  double minLoadIncr_; 
  double maxLoadVal_;  
  double initLoad_;    
 
  String varName_; 
  bool extScale_;  
};