peridigm_dem_coupler.hpp

Go to the documentation of this file.

// Originates from Fan Zhu (HKUST). Modified by Zhengshou Lai to be coupled with
// NetDEM by Sep 12, 2021
 
#pragma once
 
#include "dem_fragment.hpp"
#include "particle.hpp"
#include "peridigm_simulator.hpp"
#include "peridigm_strength.hpp"
#include "simulation.hpp"
#include <string>
 
namespace netdem {
 
class PeriDigmDEMCoupler {
public:
  std::string base_dir{"tmp/out/"};
 
  int sub_dir_index{0};
 
  Particle *particle{nullptr};
 
  PeriDigmSimulator pd_sim;
 
  STLModel surface_stl;
 
  int mesh_res{20};
 
  double node_size_ave{0.0};
 
  VecXT<int> fixed_nodes;
 
  VecXT<int> boundary_force_nodes;
 
  VecXT<double> boundary_force_node_vols;
 
  VecXT<Vec3d> boundary_force_values;
 
  VecXT<int> unbalanced_force_nodes;
 
  Vec3d unbalanced_force_values{0, 0, 0};
 
  double contact_force_max{0.0};
 
  VecXT<double> contact_force_list;
 
  bool use_customized_loading_rate{false};
 
  double loading_rate{1.0e5};
 
  int loading_steps{1000};
 
  double mech_time{0.0};
 
  bool is_broken{false};
 
  double damage_fraction_limit{0.05};
 
  double fragment_vol_limit{0.001};
 
  bool ignore_fines{true};
 
  bool use_alpha_shape{true};
 
  double fragment_alpha{0.0};
 
  PeriDigmStrength strength_params;
 
  PeriDigmMaterial material_params;
 
  PeriDigmDEMCoupler();
 
  void Init(Particle *p);
 
  void Solve();
 
  void ApplyBoundaryForce(Vec3d const &pos, Vec3d const &force);
 
  bool CheckBreakage();
 
  VecXT<DEMFragment> GetFragments();
 
private:
  double damage_limit{0.5};
 
  VecXT<VecXT<double>> damage_data;
 
  void UpdateMaterials();
 
  void UpdateMechTime();
 
  void UpdateCriticalStretch();
 
  std::string GetResultDirectory();
 
  VecXT<int> SeperateFragments(const VecXT<Vec2i> &bond_list);
 
  VecXT<int> GetFragmentNodeIndices(const VecXT<Vec2i> &bond_list,
                                    const VecXT<int> &frag_id_list,
                                    int frag_id);
 
  VecXT<Vec3d> GetFragmentNodes(const VecXT<Vec3d> &node_list,
                                const VecXT<int> &node_ids);
 
  VecXT<Vec3d> GetFragmentNodeVelocities(const VecXT<Vec3d> &velocity_list,
                                         const VecXT<int> &node_ids);
 
  VecXT<double> GetFragmentNodevolumes(const VecXT<double> &node_vol_list,
                                       const VecXT<int> &node_ids);
 
  STLModel GetAlphaShape(const VecXT<Vec3d> &point_list, double alpha);
 
  void ResolveFragmentOverlap(VecXT<DEMFragment> *const frag_list);
 
  DEMFragment GetFragmentCombined(const VecXT<DEMFragment> &frag_list);
 
  void WriteLogFileDEM();
};
 
} // namespace netdem