bonded_voronois.hpp

Go to the documentation of this file.

#pragma once
 
#include "bond_solver_pp.hpp"
#include "contact_model.hpp"
#include "contact_pp.hpp"
#include "particle.hpp"
#include "shape_trimesh.hpp"
#include "voronoi.hpp"
#include <string>
 
namespace netdem {
 
class Scene;
 
class BondedVoronois {
public:
  VecXT<TriMesh> trimesh_list;
 
  VecXT<Particle> particle_list;
 
  VecXT<ContactPP> contact_list;
 
  VecXT<Vec2i> bond_pair_list;
 
  int cvt_max_iters{1000};
 
  double cvt_tol{1.0e-3};
 
  ContactModel *bond_model{nullptr};
 
  BondedVoronois();
 
  BondedVoronois(BondedVoronois const &bp);
 
  BondedVoronois(BondedVoronois const &&bp);
 
  BondedVoronois &operator=(BondedVoronois const &bp);
 
  BondedVoronois &operator=(BondedVoronois const &&bp);
 
  void SetBondModel(ContactModel *cnt_model);
 
  void Translate(double pos_x, double pos_y, double pos_z);
 
  void RotateByRodrigues(double rot_angle, double rot_axis_x, double rot_axis_y,
                         double rot_axis_z);
 
  Vec3d GetCentroid();
 
  void InitFromSTL(std::string const &filename, int num_voros);
 
  void InitFromSTL(STLModel const &stl_model, int num_voros);
 
  void MakePorosity(double porosity);
 
  void InitBonds();
 
  void SaveAsVTK(std::string const &file_name);
 
  void ImportToScene(Scene *const scene) const;
 
private:
  VecXT<Vec3d> FindSharedVertices(STLModel const &stl_1, STLModel const &stl_2);
 
  Vec3d PolyNormal(VecXT<Vec3d> const &verts);
 
  void PolySortVertices(VecXT<Vec3d> *const pos_vec, Vec3d const &dir_n);
 
  std::tuple<Vec3d, double> PolyCentroid(const VecXT<Vec3d> &verts);
 
  void RefreshPointers();
};
 
} // namespace netdem