PeriodicBCModel.h

Go to the documentation of this file.

 
#pragma once
 
#include <jem/base/Object.h>
 
#include <cmath>
 
#include <jem/base/Float.h>
#include <jem/util/Properties.h>
#include <jem/util/PropertyException.h>
 
#include <jive/Array.h>
#include <jive/fem/NodeGroup.h>
#include <jive/fem/NodeSet.h>
#include <jive/implict/ArclenActions.h>
#include <jive/model/Actions.h>
#include <jive/model/Model.h>
#include <jive/model/ModelFactory.h>
#include <jive/model/Names.h>
#include <jive/model/StateVector.h>
#include <jive/util/Assignable.h>
#include <jive/util/Constraints.h>
#include <jive/util/DofSpace.h>
#include <jive/util/Globdat.h>
 
#include "modules/PBCGroupInputModule.h"
#include "utils/helpers.h"
#include "utils/testing.h"
 
using jem::util::Properties;
using jive::BoolMatrix;
using jive::IdxMatrix;
using jive::fem::NodeGroup;
using jive::fem::NodeSet;
using jive::implict::ArclenActions;
using jive::implict::ArclenParams;
using jive::model::Model;
using jive::model::StateVector;
using jive::util::Assignable;
using jive::util::Constraints;
using jive::util::DofSpace;
using jive::util::Globdat;
 
using jive_helpers::IdxVectorMatrix;
using jive_helpers::vec2mat;
 
//-----------------------------------------------------------------------
//   class PeriodicBCModel
//-----------------------------------------------------------------------
 
class PeriodicBCModel : public Model
{
public:
  static const char *TYPE_NAME;       
  static const char *MODE_PROP;       
  static const char *GRAD_PROP;       
  static const char *DOF_NAMES_PROP;  
  static const char *ROT_NAMES_PROP;  
  static const char *FIXEDGRAD_PARAM; 
 
  JEM_DECLARE_CLASS(PeriodicBCModel, Model);
 
  enum Mode
  {
    LOAD, 
    UPD,  
    DISP  
  };
 
  explicit PeriodicBCModel(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:
  void init_(const Properties &globdat);
 
  void fixCorners_(const Properties &globdat,
                   const Matrix &currentGrad,
                   const double scale = NAN);
 
  void setConstraints_();
 
  void getExtVec_(const Vector &f,
                  const Properties &globdat,
                  const double scale = 1.);
 
private:
  Assignable<NodeSet> nodes_;      
  Ref<DofSpace> dofs_;             
  Ref<Constraints> cons_;          
  String gradName_;                
  Matrix grad_;                    
  StringVector dofNames_;          
  StringVector rotNames_;          
  IdxVector jdofs_;                
  IdxVectorMatrix masterEdgeDofs_; 
  IdxVectorMatrix slaveEdgeDofs_;  
  IdxMatrix cornerDofs_;           
  IdxVector corner0Dofs_;          
  idx_t pbcRank_;                  
  Mode mode_;                      
  bool ghostCorners_;              
};