peridigm_settings.hpp

Go to the documentation of this file.

#pragma once
 
#include "utils_math.hpp"
#include <cmath>
#include <fstream>
#include <sstream>
#include <string>
 
namespace netdem {
 
class PeriDigmSettings {
public:
  std::string result_dir{"tmp/out/peridigm/"};
 
  std::string peridigm_exe{"Peridigm"};
 
  double horizon_factor{3.01};
 
  bool omit_bonds_between_blocks{false};
 
  bool use_auto_timestep{true};
 
  double timestep{1.0e-6};
 
  double timestep_factor{0.95};
 
  double mech_time{0.0};
 
  double loading_radius_factor{1.5};
 
  double constrain_radius_factor{1.5};
 
  int output_freqency{10};
 
  void WriteInputFile(std::ostream &os) {
    // solver settings
    os << "  Solver: " << std::endl;
    os << "    Verbose: false" << std::endl;
    os << "    Initial Time: 0.0" << std::endl;
    os << "    Final Time: " + IO::ToString(mech_time) << std::endl;
    os << "    Verlet: " << std::endl;
    if (!use_auto_timestep) {
      os << "      Fixed dt: " + IO::ToString(timestep) << std::endl;
    } else {
      os << "      Safety Factor: " + IO::ToString(timestep_factor)
         << std::endl;
    }
    os << std::endl;
 
    // compute settings, for fraguting the stored strain energy in the particle
    os << "  Compute Class Parameters: " << std::endl;
    os << "    Grain Stored Elastic Energy: " << std::endl;
    os << "      Compute Class: \"Block_Data\"" << std::endl;
    os << "      Calculation Type: \"Sum\"" << std::endl;
    os << "      Block: \"block_1\"" << std::endl;
    os << "      Variable: \"Stored_Elastic_Energy\"" << std::endl;
    os << "      Output Label: \"Grain_Stored_Elastic_Energy\"" << std::endl;
    os << std::endl;
 
    // output settings
    os << "  Output Data: " << std::endl;
    os << "    Output File Type: \"ExodusII\"" << std::endl;
    os << "    Output Format: \"BINARY\"" << std::endl;
    os << "    Output Filename: \"output\"" << std::endl;
    os << "    Output Frequency: " + IO::ToString(output_freqency)
       << std::endl;
    os << "    Parallel Write: true" << std::endl;
    os << "    Output Variables: " << std::endl;
    os << "      Displacement: true" << std::endl;
    os << "      Velocity: true" << std::endl;
    os << "      Element_Id: true" << std::endl;
    os << "      Dilatation: true" << std::endl;
    os << "      Weighted_Volume: true" << std::endl;
    os << "      Volume: true" << std::endl;
    os << "      Damage: true" << std::endl;
    os << "      Radius: true" << std::endl;
    os << "      Horizon: true" << std::endl;
    os << "      Contact_Force: true" << std::endl;
    os << "      Number_Of_Neighbors: true" << std::endl;
    os << "      Neighborhood_Volume: true" << std::endl;
    os << "      Force: true" << std::endl;
    os << "      Force_Density: true" << std::endl;
    os << "      Stored_Elastic_Energy: true" << std::endl;
    os << "      Grain_Stored_Elastic_Energy: true" << std::endl;
    os << "      Kinetic_Energy: true" << std::endl;
    os << std::endl;
  }
};
 
} // namespace netdem