mpm_particle_base_dirichlet_condition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:        BSD License
//                  Kratos default license: kratos/license.txt
//
//  Main authors:    Bodhinanda Chandra
//
 
 
#if !defined(KRATOS_MPM_PARTICLE_BASE_DIRICHLET_CONDITION_3D_H_INCLUDED )
#define      KRATOS_MPM_PARTICLE_BASE_DIRICHLET_CONDITION_3D_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "custom_conditions/particle_based_conditions/mpm_particle_base_condition.h"
#include "particle_mechanics_application_variables.h"
 
namespace Kratos
{
 
 
 
 
 
 
class MPMParticleBaseDirichletCondition
    : public MPMParticleBaseCondition
{
 
public:
 
    typedef std::size_t SizeType;
 
    // Counted pointer of MPMParticleBaseDirichletCondition
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( MPMParticleBaseDirichletCondition );
 
    using MPMParticleBaseCondition::CalculateOnIntegrationPoints;
    using MPMParticleBaseCondition::SetValuesOnIntegrationPoints;
 
 
    // Constructor void
    MPMParticleBaseDirichletCondition()
    {};
 
    // Constructor using an array of nodes
    MPMParticleBaseDirichletCondition( IndexType NewId, GeometryType::Pointer pGeometry )
        : MPMParticleBaseCondition(NewId,pGeometry)
    {};
 
    // Constructor using an array of nodes with properties
    MPMParticleBaseDirichletCondition( IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties )
        : MPMParticleBaseCondition(NewId,pGeometry,pProperties)
    {};
 
    // Destructor
    ~MPMParticleBaseDirichletCondition() override
    {};
 
 
 
    int Check( const ProcessInfo& rCurrentProcessInfo ) const override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateOnIntegrationPoints(
        const Variable<array_1d<double, 3 > >& rVariable,
        std::vector<array_1d<double, 3 > >& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void SetValuesOnIntegrationPoints(
        const Variable<array_1d<double, 3 > >& rVariable,
        const std::vector<array_1d<double, 3 > >& rValues,
        const ProcessInfo& rCurrentProcessInfo) override;
 
 
protected:
 
 
    array_1d<double, 3> m_imposed_displacement;
    array_1d<double, 3> m_imposed_velocity;
    array_1d<double, 3> m_imposed_acceleration;
    array_1d<double, 3> m_contact_force;
 
 
    virtual void MPMShapeFunctionPointValues(Vector& rResult) const override;
 
 
private:
 
    friend class Serializer;
 
    void save( Serializer& rSerializer ) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, MPMParticleBaseCondition );
        rSerializer.save("imposed_displacement",m_imposed_displacement);
        rSerializer.save("imposed_velocity",m_imposed_velocity);
        rSerializer.save("imposed_acceleration",m_imposed_acceleration);
        rSerializer.save("contact_force",m_contact_force);
 
    }
 
    void load( Serializer& rSerializer ) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, MPMParticleBaseCondition );
        rSerializer.load("imposed_displacement",m_imposed_displacement);
        rSerializer.load("imposed_velocity",m_imposed_velocity);
        rSerializer.load("imposed_acceleration",m_imposed_acceleration);
        rSerializer.load("contact_force",m_contact_force);
    }
 
 
}; // class MPMParticleBaseDirichletCondition.
 
 
} // namespace Kratos.
 
#endif // KRATOS_MPM_PARTICLE_BASE_DIRICHLET_CONDITION_3D_H_INCLUDED  defined