contact_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/contact_forces_data.hpp>
#include <meltpooldg/particles/dem_util.hpp>
#include <meltpooldg/particles/obstacle_field.hpp>
#include <meltpooldg/particles/particle_accessor.hpp>
#include <meltpooldg/time_integration/time_iterator.hpp>
 
#include <memory>
 
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  class SphericalParticleContactForce
  {
  public:
    explicit SphericalParticleContactForce(
      const SphericalParticleContactData<number>              &contact_data,
      const MeltPoolDG::TimeIntegration::TimeIterator<number> &time_iterator);
 
 
    void
    attach_wall(std::unique_ptr<dealii::Function<dim>> &&wall_signed_distance_function);
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const;
 
  private:
    const SphericalParticleContactData<number> &contact_data;
 
    const number damping_prefactor;
 
    mutable std::map<int, std::map<int, dealii::Tensor<1, dim, number>>> tangential_gaps;
 
    mutable std::map<int, std::map<int, dealii::Tensor<1, dim, number>>> tangential_gaps_with_walls;
 
    const MeltPoolDG::TimeIntegration::TimeIterator<number> &time_iterator;
 
    std::map<int, std::unique_ptr<dealii::Function<dim>>> wall_signed_distance_functions;
 
    struct ContactConfiguration
    {
      ContactConfiguration(const DEMParticleAccessor<dim, number> &self,
                           const DEMParticleAccessor<dim, number> &other,
                           const number                            youngs_modulus,
                           const number                            poisson_ratio);
 
      ContactConfiguration(const DEMParticleAccessor<dim, number> &particle,
                           const dealii::Function<dim>            *wall,
                           const number                            youngs_modulus,
                           const number                            poisson_ratio);
 
      number effective_mass;
 
      number effective_radius;
 
      number effective_youngs_modulus;
 
      number effective_shear_modulus;
 
      dealii::Tensor<1, dim, number> normal_vector;
 
      number normal_overlap;
 
      dealii::Tensor<1, axial_dim<dim>, number> relative_angular_velocity;
 
      struct
      {
        dealii::Tensor<1, dim, number> value;
 
        dealii::Tensor<1, dim, number> normal_component;
 
        dealii::Tensor<1, dim, number> tangential_component;
      } relative_velocity;
    };
 
    dealii::Tensor<1, dim, number>
    normal_contact_force(const ContactConfiguration &contact_configuration) const;
 
    dealii::Tensor<1, dim, number>
    tangential_contact_force(const ContactConfiguration           &contact_configuration,
                             const dealii::Tensor<1, dim, number> &normal_force,
                             dealii::Tensor<1, dim, number>       &tangential_gap) const;
 
    dealii::Tensor<1, axial_dim<dim>, number>
    tangential_contact_torque(const ContactConfiguration           &contact_configuration,
                              const dealii::Tensor<1, dim, number> &tangential_force,
                              const number                          particle_radius) const;
 
    dealii::Tensor<1, axial_dim<dim>, number>
    rolling_resistance_torque(const ContactConfiguration           &contact_configuration,
                              const dealii::Tensor<1, dim, number> &contact_normal_force) const;
 
    number
    compute_damping_prefactor(const number restitution_coefficient) const;
  };
} // namespace MeltPoolDG