UP_condition.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Lorenzo Gracia
//
 
 
#if !defined(KRATOS_UP_CONDITION_H_INCLUDED )
#define  KRATOS_UP_CONDITION_H_INCLUDED
 
// System includes
#include <cmath>
 
// Project includes
#include "includes/define.h"
#include "includes/condition.h"
#include "includes/serializer.h"
#include "includes/process_info.h"
 
// Application includes
#include "custom_utilities/poro_condition_utilities.hpp"
#include "dam_application_variables.h"
 
namespace Kratos
{
 
template< unsigned int TDim, unsigned int TNumNodes >
class UPCondition : public Condition
{
 
public:
 
    KRATOS_CLASS_POINTER_DEFINITION( UPCondition );
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    // Default constructor
    UPCondition() : Condition() {}
 
    // Constructor 1
    UPCondition( IndexType NewId, GeometryType::Pointer pGeometry ) : Condition(NewId, pGeometry) {}
 
    // Constructor 2
    UPCondition( IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties ) : Condition(NewId, pGeometry, pProperties)
    {
        mThisIntegrationMethod = this->GetGeometry().GetDefaultIntegrationMethod();
    }
 
    // Destructor
    virtual ~UPCondition() {}
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    Condition::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
 
    void GetDofList(DofsVectorType& rConditionDofList, const ProcessInfo& rCurrentProcessInfo) const override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void CalculateLocalSystem(MatrixType& rLeftHandSideMatrix,VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
    void EquationIdVector(EquationIdVectorType& rResult, const ProcessInfo& rCurrentProcessInfo) const override;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
protected:
 
    struct ConditionVariables
    {
        double Density;
 
        double AcelerationCoefficient;
 
        Vector Np;
        Vector NormalVector;
        BoundedMatrix<double,TDim, TNumNodes*TDim> Nu;
        double IntegrationCoefficient;
 
        array_1d<double,TNumNodes> PressureVector;
        Vector AccelerationVector;
 
        BoundedMatrix<double,TNumNodes*TDim,TNumNodes> UPMatrix;
        BoundedMatrix<double,TNumNodes,TNumNodes*TDim> PUMatrix;
        array_1d<double,TNumNodes*TDim> UVector;
        array_1d<double,TNumNodes> PVector;
    };
 
    // Member Variables
 
    GeometryData::IntegrationMethod mThisIntegrationMethod;
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
    void CalculateAll( MatrixType& rLeftHandSideMatrix, VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateLHS( MatrixType& rLeftHandSideMatrix, const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateRHS( VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo);
 
    void CalculateNormalVector(VectorType& rNormalVector,  const Matrix& Jacobian);
 
    void CalculateIntegrationCoefficient(double& rIntegrationCoefficient, const Matrix& Jacobian, const double& weight);
 
    void InitializeConditionVariables(ConditionVariables& rVariables, const GeometryType& Geom, const PropertiesType& Prop, const ProcessInfo& rCurrentProcessInfo);
 
    void GetAccelerationVector( Vector& rValues, int Step );
 
    void CalculateLHSContribution(MatrixType& rLeftHandSideMatrix, ConditionVariables& rVariables);
 
    void CalculateRHSContribution(VectorType& rRightHandSideVector, ConditionVariables& rVariables);
 
//----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
private:
 
    // Serialization
 
    friend class Serializer;
 
    void save(Serializer& rSerializer) const override
    {
        KRATOS_SERIALIZE_SAVE_BASE_CLASS( rSerializer, Condition )
    }
 
    void load(Serializer& rSerializer) override
    {
        KRATOS_SERIALIZE_LOAD_BASE_CLASS( rSerializer, Condition )
    }
 
}; // class UPCondition.
 
} // namespace Kratos.
 
#endif // KRATOS_UP_CONDITION_H_INCLUDED defined