solid_particles.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 SOLID_PARTICLES_H
#define SOLID_PARTICLES_H

#include "elastic_solid.h"
#include "base_particles.h"
#include "base_particles.hpp"

#include "particle_generator_lattice.h"
namespace SPH
{
    //----------------------------------------------------------------------
    //      preclaimed classes
    //----------------------------------------------------------------------
    class Solid;
    class ElasticSolid;

    class SolidParticles : public BaseParticles
    {
    public:
        SolidParticles(SPHBody &sph_body, Solid *solid);
        virtual ~SolidParticles(){};

        StdLargeVec<Vecd> pos0_; 
        StdLargeVec<Vecd> n_;    
        StdLargeVec<Vecd> n0_;   
        StdLargeVec<Matd> B_;    
        virtual Vecd getKernelGradient(size_t index_i, size_t index_j, Real dW_ij, Vecd &e_ij) override;
        virtual StdLargeVec<Vecd> *AverageVelocity() { return &vel_; };
        virtual StdLargeVec<Vecd> *AverageAcceleration() { return &acc_; };

        virtual void initializeOtherVariables() override;

        virtual SolidParticles *ThisObjectPtr() override { return this; };
    };

    class ElasticSolidParticles : public SolidParticles
    {
    public:
        ElasticSolidParticles(SPHBody &sph_body, ElasticSolid *elastic_solid);
        virtual ~ElasticSolidParticles(){};

        StdLargeVec<Matd> F_;     
        StdLargeVec<Matd> dF_dt_; 
        //----------------------------------------------------------------------
        //      for fluid-structure interaction (FSI)
        //----------------------------------------------------------------------
        StdLargeVec<Vecd> vel_ave_; 
        StdLargeVec<Vecd> acc_ave_; 
        // STRAIN
        Matd getGreenLagrangeStrain(size_t particle_i);
        Vecd getPrincipalStrains(size_t particle_i);
        Real getVonMisesStrain(size_t particle_i);
        Real getVonMisesStrainDynamic(size_t particle_i, Real poisson);

        StdLargeVec<Real> getVonMisesStrainVector(std::string strain_measure = "static");
        Real getVonMisesStrainMax(std::string strain_measure = "static");
        Real getPrincipalStrainMax();

        // STRESS
        Matd getStressCauchy(size_t particle_i);
        Matd getStressPK2(size_t particle_i);
        Vecd getPrincipalStresses(size_t particle_i);
        Real getVonMisesStress(size_t particle_i);

        StdLargeVec<Real> getVonMisesStressVector();
        Real getVonMisesStressMax();
        Real getPrincipalStressMax();

        Vecd displacement(size_t particle_i);
        StdLargeVec<Vecd> getDisplacement();
        Real getMaxDisplacement();

        Vecd normal(size_t particle_i);
        StdLargeVec<Vecd> getNormal();

        std::string stress_measure_;

        virtual StdLargeVec<Vecd> *AverageVelocity() { return &vel_ave_; };
        virtual StdLargeVec<Vecd> *AverageAcceleration() { return &acc_ave_; };

        virtual void initializeOtherVariables() override;
        virtual ElasticSolidParticles *ThisObjectPtr() override { return this; };

    protected:
        ElasticSolid *elastic_solid_;
    };

    class ShellParticles : public ElasticSolidParticles
    {
    public:
        ShellParticles(SPHBody &sph_body, ElasticSolid *elastic_solid);
        virtual ~ShellParticles(){};

        Real thickness_ref_;
        StdLargeVec<Matd> transformation_matrix_; 
        StdLargeVec<Real> thickness_;             
        //----------------------------------------------------------------------
        //  extra generalized coordinates in global coordinate
        //----------------------------------------------------------------------
        StdLargeVec<Vecd> pseudo_n_;      
        StdLargeVec<Vecd> dpseudo_n_dt_;  
        StdLargeVec<Vecd> dpseudo_n_d2t_; 
        //----------------------------------------------------------------------
        //  extra generalized coordinate and velocity in local coordinate
        //----------------------------------------------------------------------
        StdLargeVec<Vecd> rotation_;        
        StdLargeVec<Vecd> angular_vel_;     
        StdLargeVec<Vecd> dangular_vel_dt_; 
        //----------------------------------------------------------------------
        //  extra deformation and deformation rate in local coordinate
        //----------------------------------------------------------------------
        StdLargeVec<Matd> F_bending_;     
        StdLargeVec<Matd> dF_bending_dt_; 
        //----------------------------------------------------------------------
        //  extra stress for pair interaction in global coordinate
        //----------------------------------------------------------------------
        StdLargeVec<Vecd> global_shear_stress_; 
        StdLargeVec<Matd> global_stress_;       
        StdLargeVec<Matd> global_moment_;       
        virtual Real ParticleVolume(size_t index_i) { return Vol_[index_i] * thickness_[index_i]; }
        virtual Real ParticleMass(size_t index_i) { return mass_[index_i] * thickness_[index_i]; }

        virtual void initializeOtherVariables() override;
        virtual ShellParticles *ThisObjectPtr() override { return this; };
    };
}
#endif // SOLID_PARTICLES_H