peridigm_boundary_condition.hpp

Go to the documentation of this file.

#pragma once
 
#include "utils_io.hpp"
#include "utils_math.hpp"
#include <fstream>
#include <sstream>
#include <string>
 
namespace netdem {
 
class PeriDigmBoundaryCondition {
public:
  enum Type { Prescribed_Displacement, Body_Force };
 
  Type type{Type::Prescribed_Displacement};
 
  VecXT<int> node_indices;
 
  VecNT<bool, 3> dim_activated{true, true, true};
 
  bool time_depedent{true};
 
  Vec3d disp_rate{0, 0, 0};
 
  Vec3d loading_rate{0, 0, 0};
 
  Vec3d disp{0, 0, 0};
 
  Vec3d loading{0, 0, 0};
 
  double mech_time{0};
 
  void InsertNode(int node_set_id);
 
  void SetActivatedDimensions(bool x, bool y, bool z);
 
  void SetByDisplacementRate(double x, double y, double z);
 
  void SetByUltimateDisplacement(double x, double y, double z, double t);
 
  void SetByLoadingRate(double x, double y, double z);
 
  void SetByUltimateLoading(double x, double y, double z, double t);
 
  void WriteInputFile(std::ostream &os, int node_set_id);
 
  void WriteNodeSetFile(std::string const &result_dir, int node_set_id);
 
private:
  std::string GetNodeSetFileName(int node_set_id);
 
  std::string GetDisplacementString(int dim);
 
  std::string GetLoadingString(int dim);
};
} // namespace netdem