rnpBond.h

//  Copyright (c) 2019 Prashant K. Jha
//  Copyright (c) 2019 Patrick Diehl
//
//  Distributed under the Boost Software License, Version 1.0. (See accompanying
//  file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
 
#ifndef MATERIAL_PD_BONDMATERIAL_NEW_H
#define MATERIAL_PD_BONDMATERIAL_NEW_H
 
#include "baseMaterial.h"
#include "util/point.h"
#include "util/matrixBlaze.h"
#include "geometry/neighbor.h"
#include "fe/mesh.h"
 
#include <vector>
 
#include <hpx/include/parallel_algorithm.hpp>
 
// forward declaration
namespace inp {
struct MaterialDeck;
struct OutputDeck;
}
 
// forward declaration
namespace data {
class DataManager;
}
 
// forward declaration of class
namespace fe {
class Mesh;
} // namespace fe
 
namespace material {
namespace pd {
class BaseInfluenceFn;
} // namespace pd
} // namespace material
 
namespace material {
 
namespace pd {
 
class RNPBond : public BaseMaterial {
 
public:
  RNPBond(inp::MaterialDeck *deck,data::DataManager* dataManager);
 
  std::pair<util::Point3, double> getBondEF(size_t i, size_t j);
 
  double getS(const util::Point3 &dx, const util::Point3 &du);
 
  double getS(size_t i, size_t j);
 
  double getSc(const double &r);
 
  double getSc(size_t i , size_t j);
 
  util::Point3 getBondForceDirection(const util::Point3 &dx,
                                     const util::Point3 &du) const;
 
  double getInfFn(const double &r) const ;
 
private:
  void computeParameters(inp::MaterialDeck *deck, const double &M);
 
  void computeMaterialProperties(inp::MaterialDeck *deck, const double &M);
 
  double d_C;
 
  double d_beta;
 
  double d_rbar;
 
  double d_invFactor;
 
  double d_factorSc;
 
  bool d_irrevBondBreak;
 
  double d_contact_Kn = 0.;
 
  double d_contact_Rc = 0.;
 
  const inp::MaterialDeck *d_deck;
 
  data::DataManager* d_dataManager_p;
 
  material::pd::BaseInfluenceFn *d_baseInfluenceFn_p;
};
 
} // namespace pd
 
} // namespace material
 
#endif // MATERIAL_PD_BONDMATERIAL_NEW_H