deformation_analysis.hpp

Go to the documentation of this file.

#pragma once
 
#include "deformable_particle.hpp"
#include "modifier.hpp"
#include "particle_group.hpp"
#include <unordered_map>
 
namespace netdem {
 
class DeformationAnalysis : public ParticleGroup {
private:
  enum SolveBy { cycles, time };
 
  // 0: cycles, 1: time
  SolveBy solve_by{SolveBy::cycles};
 
  int solve_cycle_interval{0};
  double solve_time_interval{0};
 
  int solve_cycle_previous{0};
  double solve_time_previous{0};
 
  using SaveBy = SolveBy;
 
  SaveBy save_by{SaveBy::cycles};
 
  int save_cycle_interval{0};
  double save_time_interval{0};
 
  int save_cycle_previous{0};
  double save_time_previous{0};
 
public:
  double neo_k{6.94e5}, neo_mu{5.21e5}, density{500.0};
  double damp_coef{0.7};
 
  Vec3d gravity_coef{0.0, 0.0, -9.81};
  double timestep{1.0e-4};
  int mesh_res{20};
 
  std::string root_path{"tmp/out/deformation_analysis/"};
 
  std::unordered_map<Particle *, std::pair<bool, DeformableParticle>>
      particle_map;
 
public:
  DeformationAnalysis();
 
  Modifier *Clone() const override;
 
  void Init(Simulation *sim) override;
 
  void SetSolveByTime(double interval);
 
  void SetSolveByCycles(double interval);
 
  void SetSaveByTime(double interval);
 
  void SetSaveByCycles(double interval);
 
  void Execute() override;
 
  void Update() override;
 
private:
  void SolveDeformation();
 
  void SaveFEMAsVTK();
 
  void EvaluateBCForce(DeformableParticle *const p_deformable_ptr,
                       ContactPP *const cnt);
 
  void EvaluateBCForce(DeformableParticle *const p_deformable_ptr,
                       ContactPW *const cnt);
 
  void SetSettings(FEMSimulator *const fem_sim);
 
  std::string GetFEMResultFileName(Particle *const p_ptr);
 
  bool CheckIfToExecute();
};
 
} // namespace netdem