solver_sdf_pw.hpp

Go to the documentation of this file.

#pragma once
 
#include "collision_geometries.hpp"
#include "collision_solver_pw.hpp"
#include "particle.hpp"
#include "shape.hpp"
 
namespace netdem {
 
class SolverSDFPW : public CollisionSolverPW {
public:
  enum PotentialType { linear, hertz };
 
  int potential_type{PotentialType::linear};
 
  bool solve_two_sides{false};
 
  bool use_equivalent_stiffness{true};
 
  double beta{0.65};
 
  SolverSDFPW();
 
  SolverSDFPW(Particle *const p, Wall *const w);
 
  CollisionSolverPW *Clone() const override;
 
  void Init(Particle *const p, Wall *const w) override;
 
  bool Detect() override;
 
  void ResolveInit(ContactPW *const cnt, double timestep) override;
 
  void ResolveUpdate(ContactPW *const cnt, double timestep) override;
 
  void ResolveInit(CollisionGeometries *const cnt_geoms, double timestep);
 
  void ResolveUpdate(CollisionGeometries *const cnt_geoms, double timestep);
 
  void ResolveInitWToP(CollisionGeometries *const cnt_geoms, double timestep);
 
  void ResolveUpdateWToP(CollisionGeometries *const cnt_geoms, double timestep);
 
  void ResolveInitPToW(CollisionGeometries *const cnt_geoms, double timestep);
 
  void ResolveUpdatePToW(CollisionGeometries *const cnt_geoms, double timestep);
 
  double GetPotential(double dist, double node_area, double size_eq,
                      double size_ratio);
 
  static double KiToKn(double ki, double size_eq, double size_ratio);
 
  static double KiFromKn(double kn, double size_eq, double size_ratio);
 
  static double KiToEstar(double ki, double size_eq, double size_ratio);
 
  static double KiFromEstar(double Estar, double size_eq, double size_ratio);
 
private:
  bool solve_1to2{false};
 
  bool solve_2to1{false};
 
  double bound_sphere_radius_1;
 
  double bound_sphere_radius_2;
 
  // particle as reference
 
  Vec3d pos_1;
 
  Vec3d dpos_12;
 
  Vec3d dpos_12_ref;
 
  Vec4d quat_1;
 
  Vec4d dquat_12;
 
  // wall as reference
 
  Vec3d pos_2;
 
  Vec3d dpos_21;
 
  Vec3d dpos_21_ref;
 
  Vec4d quat_2;
 
  Vec4d dquat_21;
 
  Shape *shape_1{nullptr};
 
  Shape *shape_2{nullptr};
 
  VecXT<int> node_id_12_list;
  VecXT<int> node_id_21_list;
 
  VecXT<double> node_dist_12_list;
  VecXT<double> node_dist_21_list;
};
 
} // namespace netdem