BodyDelaunayMassSpringDamper.hpp

Go to the documentation of this file.

#ifndef BODYMASSSPRINGDAMPER_HPP
#define BODYMASSSPRINGDAMPER_HPP
 
#define MOVABLE
 
#include <mecacell/mecacell.h>
#include <mecacell/movable.h>
#include "PluginDelaunayMassSpringDamper.hpp"
#include "../../../../src/core/BaseBody.hpp"
#include <math.h>
 
namespace DelaunayMassSpringDamper {
 
    template<typename cell_t, class plugin_t>
    class BodyDelaunayMassSpringDamper : public MecaCell::Movable, public virtual BaseBody<plugin_t> {
 
    private:
        double baseRadius; 
        double radius; 
        double density; 
        double adhesion; 
    public:
        explicit BodyDelaunayMassSpringDamper(const MecaCell::Vector3D &pos = MecaCell::Vector3D::zero()) {
            baseRadius = 10.0; // 10 µm radius
            radius = baseRadius;
            density = 0.001; // 0.001 ng/µm³ density of water ~ 1000 kg/m³
            adhesion = 0.75;
            this->setMass(density * 4 * M_PI / 3 * pow(radius, 3)); // mass = density * 4π/3 * radius³ ng
            this->setPosition(pos);
        }
 
        inline double getBoundingBoxRadius() const { return radius; }
 
        inline void setRadius(double rad) {
            radius = rad;
            this->setMass(density * 4 * M_PI / 3 * pow(radius, 3)); // mass = density * 4π/3 * radius³ ng
        }
 
        inline double getRadius() { return radius; }
 
        inline double getBaseRadius() const { return baseRadius; }
 
        inline void setBaseRadius(double _baseRadius) { baseRadius = _baseRadius; }
 
        inline void growth(double delta) {
            this->setRadius(baseRadius * std::cbrt(1 + delta));
        }
 
        inline void setDensity(double d) {
            density = d;
            this->setMass(density * 4 * M_PI / 3 * pow(radius, 3)); // mass = density * 4π/3 * radius³ ng
        }
 
        inline double getAdhesion() const { return adhesion; }
 
        inline void setAdhesion(double a) { adhesion = a; }
 
        inline void moveTo(const MecaCell::Vec &v) { this->setPosition(v); }
    };
}
 
#endif // BODYMASSSPRINGDAMPER_HPP