complex_body.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
 *                              SPHinXsys                                   *
 * --------------------------------------------------------------------------*
 * SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle    *
 * Hydrodynamics for industrial compleX systems. It provides C++ APIs for   *
 * physical accurate simulation and aims to model coupled industrial dynamic *
 * systems including fluid, solid, multi-body dynamics and beyond with SPH  *
 * (smoothed particle hydrodynamics), a meshless computational method using *
 * particle discretization.                                                 *
 *                                                                          *
 * SPHinXsys is partially funded by German Research Foundation              *
 * (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1                *
 * and HU1527/12-1.                                                         *
 *                                                                           *
 * Portions copyright (c) 2017-2020 Technical University of Munich and      *
 * the authors' affiliations.                                               *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may   *
 * not use this file except in compliance with the License. You may obtain a *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
 *                                                                           *
 * --------------------------------------------------------------------------*/
#ifndef COMPLEX_BODY_H
#define COMPLEX_BODY_H

#include "base_body.h"
#include "base_body_part.h"
#include "base_geometry.h"

namespace SPH
{
    class SecondaryStructure
    {
    public:
        explicit SecondaryStructure(){};
        virtual ~SecondaryStructure(){};

        virtual void buildParticleConfiguration(ParticleConfiguration &particle_configuration) = 0;
    };

    class TreeBody : public SecondaryStructure, public RealBody
    {
    public:
        class Branch;

    private:
        UniquePtrKeepers<Branch> branches_ptr_keeper_;

    public:
        StdVec<Branch *> branches_;    
        IndexVector branch_locations_; 
        size_t last_branch_id_;
        Branch *root_;

        explicit TreeBody(SPHSystem &sph_system, SharedPtr<Shape> shape_ptr);
        virtual ~TreeBody(){};

        Branch *createANewBranch(size_t parent_id)
        {
            return branches_ptr_keeper_.createPtr<Branch>(parent_id, this);
        };
        size_t BranchLocation(size_t particle_idx);
        Branch *LastBranch() { return branches_[last_branch_id_]; };

        virtual void buildParticleConfiguration(ParticleConfiguration &particle_configuration) override;
        size_t ContainerSize() { return branches_.size(); };
    };

    class TreeBody::Branch : public Edge<size_t, IndexVector>
    {
    public:
        explicit Branch(TreeBody *tree);
        Branch(size_t parent_id, TreeBody *tree);
        virtual ~Branch(){};

        Vecd end_direction_; 
        IndexVector inner_particles_;
        bool is_terminated_; 
    };

    class TreeTerminates : public BodyPartByParticle
    {
    public:
        TreeBody &tree_;

        explicit TreeTerminates(SPHBody &sph_body);
        virtual ~TreeTerminates(){};
    };
}
#endif // COMPLEX_BODY_H