Interaction.h

Go to the documentation of this file.

/*******************************************************************************
 * Copyright (C) 2017-2026 Theodore Chang
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 ******************************************************************************/
#ifndef INTERACTION_H
#define INTERACTION_H
 
#include <Domain/Tag.h>
#include <set>
 
class DomainBase;
class Element;
 
class InteractionPair {
    const shared_ptr<Element> object_i, object_j;
 
    bool inertial = false;
    unsigned dimension{};
 
public:
    double effective_mass{};
    double effective_radius{};
    double effective_modulus{};
    double effective_damping{};
 
    InteractionPair(const shared_ptr<Element>&, const shared_ptr<Element>&);
 
    void set_dimension(const unsigned dim) { dimension = dim; }
    void set_inertial(const bool flag) { inertial = flag; }
 
    [[nodiscard]] vec position_i() const;
    [[nodiscard]] vec position_j() const;
 
    [[nodiscard]] vec velocity_i() const;
    [[nodiscard]] vec velocity_j() const;
 
    [[nodiscard]] const uvec& dof_i() const;
    [[nodiscard]] const uvec& dof_j() const;
 
    [[nodiscard]] double initial_gap() const;
};
 
class Interaction : public CopyableTag {
protected:
    std::weak_ptr<DomainBase> domain;
 
public:
    using CopyableTag::CopyableTag;
 
    void initialize(const shared_ptr<DomainBase>&);
 
    virtual void apply(bool, const shared_ptr<Element>&) const = 0;
    virtual void apply(bool, const shared_ptr<InteractionPair>&) const = 0;
};
 
class Hertzian final : public Interaction {
    static constexpr double four_third = 4. / 3.;
 
public:
    using Interaction::Interaction;
 
    void apply(bool, const shared_ptr<Element>&) const override {}
    void apply(bool, const shared_ptr<InteractionPair>&) const override;
};
 
class FixedParticle final : public Interaction {
    std::set<unsigned> particles;
    const double multiplier;
 
public:
    FixedParticle(unsigned, double, std::set<unsigned>&&);
 
    void apply(bool, const shared_ptr<Element>&) const override;
    void apply(bool, const shared_ptr<InteractionPair>&) const override {}
};
 
#endif