contactDeck.h

Go to the documentation of this file.

/*
 * -------------------------------------------
 * Copyright (c) 2021 - 2024 Prashant K. Jha
 * -------------------------------------------
 * PeriDEM https://github.com/prashjha/PeriDEM
 *
 * Distributed under the Boost Software License, Version 1.0. (See accompanying
 * file LICENSE)
 */
 
#ifndef INP_CONTACTDECK_H
#define INP_CONTACTDECK_H
 
namespace inp {
 
struct ContactPairDeck {
 
  double d_contactR;
 
  bool d_computeContactR;
 
 
  double d_vMax;
  double d_deltaMax;
  double d_Kn;
 
 
  double d_eps;
  double d_betan;
 
 
  double d_mu;
  bool d_dampingOn;
  bool d_frictionOn;
  double d_KnFactor;
  double d_betanFactor;
  double d_kappa;
 
  ContactPairDeck()
      : d_contactR(0.), d_computeContactR(true), d_vMax(0.), d_deltaMax(0.),
        d_Kn(0.), d_eps(0.), d_betan(0.), d_mu(0.), d_dampingOn(true),
        d_frictionOn(true), d_KnFactor(1.), d_betanFactor(1.), d_kappa(1.) {};
 
  ContactPairDeck(const ContactPairDeck &cd)
      : d_contactR(cd.d_contactR), d_computeContactR(cd.d_computeContactR),
        d_vMax(cd.d_vMax), d_deltaMax(cd.d_deltaMax),
        d_Kn(cd.d_Kn), d_eps(cd.d_eps), d_betan(cd.d_betan),
        d_mu(cd.d_mu), d_dampingOn(cd.d_dampingOn),
        d_frictionOn(cd.d_frictionOn), d_KnFactor(cd.d_KnFactor),
        d_betanFactor(cd.d_betanFactor), d_kappa(cd.d_kappa){};
 
  std::string printStr(int nt = 0, int lvl = 0) const {
 
    auto tabS = util::io::getTabS(nt);
    std::ostringstream oss;
    oss << tabS << "------- ContactPairDeck --------" << std::endl << std::endl;
    oss << tabS << "Contact radius = " << d_contactR << std::endl;
    oss << tabS << "v_max = " << d_vMax << ", Delta_max = " << d_deltaMax
              << ", Kn = " << d_Kn << std::endl;
    oss << tabS << "epsilon = " << d_eps << ", Beta_n = " << d_betan << std::endl;
    oss << tabS << "Friction coefficient = " << d_mu << std::endl;
    oss << tabS << "Damping status = " << d_dampingOn << std::endl;
    oss << tabS << "Kn factor = " << d_KnFactor
              << ", Beta n factor = " << d_betanFactor << std::endl;
    oss << tabS << "Bulk modulus = " << d_kappa << std::endl;
    oss << tabS << std::endl;
    return oss.str();
  }
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); }
 
  double getKn(const double &v1, const double &v2) const {
 
    return d_Kn * (v1 * v2) / (v1 + v2);
  }
 
  double getKn(const double &v1, const double &v2) {
 
    return d_Kn * (v1 * v2) / (v1 + v2);
  }
 
  double getWKn(const double &v) const {
 
    return d_Kn * v;
  }
 
  double getWKn(const double &v) {
 
    return d_Kn * v;
  }
 
  double getBetan(const double &v1, const double &v2) const {
 
    return d_betan * std::sqrt((v1 * v2) / (v1 + v2));
  }
 
  double getBetan(const double &v1, const double &v2) {
 
    return d_betan * std::sqrt((v1 * v2) / (v1 + v2));
  }
 
  double getWBetan(const double &v) const {
 
    return d_betan * std::sqrt(v);
  }
 
  double getWBetan(const double &v) {
 
    return d_betan * std::sqrt(v);
  }
};
 
struct ContactDeck {
 
  std::vector<std::vector<ContactPairDeck>> d_data;
 
  const ContactPairDeck &getContact(const size_t &i, const size_t &j) const {
    //return d_data[i < j ? i : j][i < j ? j : i];
    return d_data[i][j];
  }
 
  ContactPairDeck &getContact(const size_t &i, const size_t &j) {
    //return d_data[i < j ? i : j][i < j ? j : i];
    return d_data[i][j];
  }
 
  std::string printStr(int nt = 0, int lvl = 0) const {
 
    auto tabS = util::io::getTabS(nt);
    std::ostringstream oss;
    oss << tabS << "------- ContactDeck --------" << std::endl << std::endl;
    for (size_t i =0; i<d_data.size(); i++) {
      for (size_t j = 0; j < d_data.size(); j++) {
        oss << tabS << "ContactPairDeck id = (" << i << "," << j << ") info:"
                  << std::endl;
        oss << getContact(i,j).printStr(nt+2, lvl);
      }
    }
 
    oss << tabS << std::endl;
 
    return oss.str();
  }
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); }
};
 
} // namespace inp
 
#endif // INP_CONTACTDECK_H