cluster3D.h

Go to the documentation of this file.

// Last Modified by: Salva, latorre@cimne.upc.edu
 
#if !defined KRATOS_CLUSTER3D_H_INCLUDED
#define KRATOS_CLUSTER3D_H_INCLUDED
 
// System includes
#include <string>
#include <iostream>
 
// Project includes
#include "includes/define.h"
#include "includes/node.h"
#include "includes/element.h"
#include "geometries/geometry.h"
#include "includes/properties.h"
#include "includes/process_info.h"
#include "includes/indexed_object.h"
#include "containers/global_pointers_vector.h"
#include "includes/constitutive_law.h"
#include "includes/condition.h"
#include "custom_elements/spheric_particle.h"
#include "custom_utilities/create_and_destroy.h"
#include "utilities/quaternion.h"
#include "custom_elements/rigid_body_element.h"
 
namespace Kratos
{
    class Element;
    class KRATOS_API(DEM_APPLICATION) Cluster3D : public RigidBodyElement3D {
 
    public:
        KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(Cluster3D);
 
        Cluster3D();
        Cluster3D(IndexType NewId, GeometryType::Pointer pGeometry);
        Cluster3D(IndexType NewId, NodesArrayType const& ThisNodes);
        Cluster3D(IndexType NewId, GeometryType::Pointer pGeometry,  PropertiesType::Pointer pProperties);
        Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
 
        virtual ~Cluster3D();
 
        using Element::Initialize;
        void Initialize(const ProcessInfo& r_process_info) override;
        virtual void CreateParticles(ParticleCreatorDestructor* p_creator_destructor, ModelPart& dem_model_part, PropertiesProxy* p_fast_properties, const bool continuum_strategy);
        virtual void GetClustersForce(const array_1d<double,3>& gravity);
        virtual void CollectForcesAndTorquesFromSpheres();
        unsigned int GetNumberOfSpheres() { return mListOfSphericParticles.size(); };
        std::vector<SphericParticle*>  GetSpheres() { return mListOfSphericParticles; };
        virtual void SetContinuumGroupToBreakableClusterSpheres(const int Id);
        virtual void SetInitialConditionsToSpheres(const array_1d<double,3>& velocity);
        virtual void SetInitialNeighbours(const double search_increment);
        virtual void CreateContinuumConstitutiveLaws();
        virtual void Calculate(const Variable<double>& rVariable, double& Output, const ProcessInfo& r_process_info) override;
 
        virtual double SlowGetDensity();
        virtual int SlowGetParticleMaterial();
 
        double GetProgrammedDestructionTime() const {return mProgrammedDestructionTime;}
        void SetProgrammedDestructionTime(const double destruction_time){mProgrammedDestructionTime = destruction_time;}
 
        virtual std::string Info() const override
        {
        std::stringstream buffer;
        buffer << "Discrete Element #" << Id();
        return buffer.str();
        }
 
        virtual void PrintInfo(std::ostream& rOStream) const override
        {
        rOStream << "Discrete Element #" << Id();
        }
 
        virtual void PrintData(std::ostream& rOStream) const override
        {
        //mpGeometry->PrintData(rOStream);
        }
 
    protected:
 
        std::vector<double> mListOfRadii;
        std::vector<SphericParticle*> mListOfSphericParticles;
        double mProgrammedDestructionTime = -1.0; // set to a negative value, so that when marked TO_ERASE, elimination is by default.
 
    private:
 
        friend class Serializer;
        virtual void save(Serializer& rSerializer) const override { KRATOS_SERIALIZE_SAVE_BASE_CLASS(rSerializer, RigidBodyElement3D); }
        virtual void load(Serializer& rSerializer) override { KRATOS_SERIALIZE_LOAD_BASE_CLASS(rSerializer, RigidBodyElement3D); }
 
    }; // Class Cluster3D
 
    inline std::istream& operator >> (std::istream& rIStream, Cluster3D& rThis);
 
    inline std::ostream& operator << (std::ostream& rOStream, const Cluster3D& rThis)
    {
        rThis.PrintInfo(rOStream);
        rOStream << std::endl;
        rThis.PrintData(rOStream);
        return rOStream;
    }
 
}  // namespace Kratos
 
#endif // KRATOS_CLUSTER3D_INCLUDED  defined