ParticleRelaxtionOnMesh.h

Go to the documentation of this file.

 
 
#pragma once
 
#include "Topology/PointSet.h"
#include "ParticleSystem/ParticleSystem.h"
#include "ParticleSystem/Emitters/ParticleEmitter.h"
#include "./ParticleSystem/Module/IterativeDensitySolver.h"
#include "./ParticleSystem/Module/ImplicitViscosity.h"
#include "./ParticleSystem/Module/ParticleIntegrator.h"
#include "./ParticleSystem/Module/NormalForce.h"
#include "Collision/NeighborTriangleQuery.h"
#include "TriangularMeshConstraint.h"
#include "Collision/NeighborPointQuery.h"
#include "Topology/TriangleSet.h"
 
#include "Samplers/PointsBehindMesh.h"
 
namespace dyno
{
 
    /*
    *   @brief£ºGenerate points benhind triangle meshes, 
    *   @biref: and slightly shift the particles benhind triangle meshes to obtain more regular distribution.
    *
    */
 
 
    template<typename TDataType>
    class ParticleRelaxtionOnMesh : public PointsBehindMesh<TDataType>
    {
        DECLARE_TCLASS(ParticleRelaxtionOnMesh, TDataType)
 
    public:
        typedef typename TDataType::Real Real;
        typedef typename TDataType::Coord Coord;
 
        ParticleRelaxtionOnMesh();
 
        ~ParticleRelaxtionOnMesh() override;
 
        DEF_VAR(int, IterationNumber, 30, "");
 
        DEF_VAR(Real, NormalForceStrength, 0.05, "");
 
        DEF_VAR(Real, MeshCollisionThickness, 0.003, "");
 
        DEF_VAR(Real, PointNeighborLength, 0.006, "");
 
        DEF_VAR(Real, MeshNeighborLength, 0.012, "");
 
        DEF_VAR(Real, ViscosityStrength, 50.0f, "");
 
        DEF_VAR_STATE(Real, Delta, 0.03, "");
 
        DEF_VAR(int, DensityIteration, 5, "");
 
    protected:
 
        void resetStates() override;
 
        void preUpdateStates();
 
        void particleRelaxion();
 
        void updatePositions();
    private:
 
        void loadInitialStates();
 
        DEF_ARRAY_STATE(Coord, Velocity, DeviceType::GPU, "Particle velocity");
 
        DEF_ARRAY_STATE(Coord, Force, DeviceType::GPU, "Force on each particle");
 
    
        std::shared_ptr<NeighborPointQuery<TDataType>> ptr_nbrQuery;
        std::shared_ptr<IterativeDensitySolver<TDataType>> ptr_density;
        std::shared_ptr<ImplicitViscosity<TDataType>> ptr_viscosity;
        std::shared_ptr<NeighborTriangleQuery<TDataType>> ptr_nbrQueryTri;
        std::shared_ptr<TriangularMeshConstraint<TDataType>> ptr_meshCollision;
        std::shared_ptr<NormalForce<TDataType>> ptr_normalForce;
        
    };
 
    
}