cohesive_forces.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/parameter_handler.h>
#include <deal.II/base/tensor.h>
 
#include <meltpooldg/particles/cohesive_forces_data.hpp>
#include <meltpooldg/particles/obstacle_field.hpp>
#include <meltpooldg/particles/particle_accessor.hpp>
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  class SphericalParticleCohesiveForce
  {
  public:
    explicit SphericalParticleCohesiveForce(
      const SphericalParticleCohesiveForceData<number> &cohesive_force_data);
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const;
 
  private:
    SphericalParticleCohesiveForceData<number> cohesive_force_data;
 
    struct CohesiveContactConfiguration
    {
      CohesiveContactConfiguration(
        const DEMParticleAccessor<dim, number>           &self,
        const DEMParticleAccessor<dim, number>           &other,
        const SphericalParticleCohesiveForceData<number> &cohesive_force_data);
 
      number normal_distance;
 
      number effective_radius;
 
      number pull_off_force_magnitude;
 
      number cut_off_distance; // denoted as g* in literature
 
      number pull_off_force_limit_distance; // denoted as g_0 in literature
 
      dealii::Tensor<1, dim, number> normal_vector;
    };
  };
 
} // namespace MeltPoolDG