materialDeck.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 INP_MATERIALDECK_H
#define INP_MATERIALDECK_H
 
#include "util/utilIO.h"
#include <cmath>
#include <string>
#include <vector>
 
namespace inp {
 
struct MatData {
 
  double d_E;
 
  double d_G;
 
  double d_K;
 
  double d_nu;
 
  double d_lambda;
 
  double d_mu;
 
  double d_KIc;
 
  double d_Gc;
 
  MatData()
      : d_E(-1.), d_G(-1.), d_K(-1.), d_nu(-1.), d_lambda(-1.), d_mu(-1.),
        d_KIc(-1.), d_Gc(-1.){};
 
  std::string printStr(int nt = 0, int lvl = 0) const {
 
    auto tabS = util::io::getTabS(nt);
    std::ostringstream oss;
    oss << tabS << "------- MatData --------" << std::endl << std::endl;
    oss << tabS << "Young's modulus = " << d_E << std::endl;
    oss << tabS << "Shear modulus = " << d_G << std::endl;
    oss << tabS << "Bulk modulus = " << d_K << std::endl;
    oss << tabS << "Poisson ratio = " << d_nu << std::endl;
    oss << tabS << "Lame parameter Lambda = " << d_lambda << std::endl;
    oss << tabS << "Lame parameter Mu = " << d_mu << std::endl;
    oss << tabS << "Critical stress intensity factor = " << d_KIc << std::endl;
    oss << tabS << "Critical energy release rate = " << d_Gc << std::endl;
    oss << tabS << std::endl;
 
    return oss.str();
  }
 
  void print(int nt = 0, int lvl = 0) const { std::cout << printStr(nt, lvl); }
 
  double toNu(double lambda, double mu) { return lambda * 0.5 / (lambda + mu); }
 
  double toNuClassical(double K, double G) { return (3*K - 2*G) / (2*(3*K+G)); }
 
  double toE(double K, double nu) { return K * (3. * (1. - 2. * nu)); }
 
  double toK(double E, double nu) { return E / (3. * (1. - 2. * nu)); }
 
  double toLambdaE(double E, double nu) {
    return E * nu / ((1. + nu) * (1. - 2. * nu));
  }
 
  double toLambdaK(double K, double nu) { return 3. * K * nu / (1. + nu); }
 
  double toGE(double E, double nu) { return E / (2. * (1. + nu)); }
 
  double toGK(double K, double nu) {
    return 3. * K * (1. - 2. * nu) / (2. * (1. + nu));
  }
 
  double toELambda(double lambda, double nu) {
    return lambda * (1. + nu) * (1. - 2. * nu) / nu;
  }
 
  double toGc(double KIc, double nu, double E) { return KIc * KIc / E; }
 
  double toKIc(double Gc, double nu, double E) { return std::sqrt(Gc * E); }
 
};
 
struct MaterialDeck {
 
  bool d_isPlaneStrain;
 
  std::string d_materialType;
 
  size_t d_bondPotentialType;
 
  size_t d_statePotentialType;
 
  size_t d_influenceFnType;
 
  std::vector<double> d_bondPotentialParams;
 
  std::vector<double> d_statePotentialParams;
 
  std::vector<double> d_influenceFnParams;
 
  bool d_irreversibleBondBreak;
 
  bool d_stateContributionFromBrokenBond;
 
  double d_checkScFactor;
 
  bool d_computeParamsFromElastic;
 
  inp::MatData d_matData;
 
  double d_density;
 
  bool d_applyContact;
 
  MaterialDeck()
      : d_isPlaneStrain(false), d_bondPotentialType(0), d_statePotentialType(0),
        d_influenceFnType(0), d_irreversibleBondBreak(true),
        d_stateContributionFromBrokenBond(true), d_checkScFactor(1.),
        d_computeParamsFromElastic(true), d_matData(inp::MatData()),
        d_density(1.), d_applyContact(false){};
 
  std::string printStr(int nt = 0, int lvl = 0) const {
 
    auto tabS = util::io::getTabS(nt);
    std::ostringstream oss;
    oss << tabS << "------- MaterialDeck --------" << std::endl << std::endl;
    oss << tabS << "Is plain strain = " << d_isPlaneStrain << std::endl;
    oss << tabS << "Material type = " << d_materialType << std::endl;
    oss << tabS << "Bond potential type = " << d_bondPotentialType << std::endl;
    oss << tabS << "Bond potential params = ["
        << util::io::printStr<double>(d_bondPotentialParams, 0) << "]"
        << std::endl;
    oss << tabS << "State potential type = " << d_statePotentialType
        << std::endl;
    oss << tabS << "State potential params = ["
        << util::io::printStr<double>(d_statePotentialParams, 0) << "]"
        << std::endl;
    oss << tabS << "Influence function type = " << d_influenceFnType
        << std::endl;
    oss << tabS << "Influence function params = ["
        << util::io::printStr<double>(d_influenceFnParams, 0) << "]"
        << std::endl;
    oss << tabS
        << "Irreversible bond breaking enabled = " << d_irreversibleBondBreak
        << std::endl;
    oss << tabS << "State contribution from broken bond enabled = "
        << d_stateContributionFromBrokenBond << std::endl;
    oss << tabS << "Check Sc factor = " << d_checkScFactor << std::endl;
    oss << tabS << "Compute parameters from elastic properties = "
        << d_computeParamsFromElastic << std::endl;
    oss << d_matData.printStr(nt + 1, lvl);
    oss << tabS << "Density = " << d_density << std::endl;
    oss << tabS << "Apply contact on broken bonds = " << d_applyContact << std::endl;
    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_MATERIALDECK_H