updated_lagrangian_U_J_W_wP_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosPfemSolidMechanicsApplication $
//   Created by:          $Author:                     PNavas $
//   Last modified by:    $Co-Author:               LMonforte $
//   Date:                $Date:                 October 2017 $
//   Revision:            $Revision:                     -0.1 $
//
//
 
#if !defined(KRATOS_UPDATED_LAGRANGIAN_U_J_W_WP_ELEMENT_H_INCLUDED )
#define  KRATOS_UPDATED_LAGRANGIAN_U_J_W_WP_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/updated_lagrangian_U_J_element.hpp"
 
namespace Kratos
{
 
 
 
class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) UpdatedLagrangianUJWwPElement
    : public UpdatedLagrangianUJElement
{
public:
 
    typedef ConstitutiveLaw ConstitutiveLawType;
    typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
    typedef ConstitutiveLawType::StressMeasure StressMeasureType;
    typedef GeometryData::IntegrationMethod IntegrationMethod;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( UpdatedLagrangianUJWwPElement );
 
 
    UpdatedLagrangianUJWwPElement();
 
    UpdatedLagrangianUJWwPElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
    UpdatedLagrangianUJWwPElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
    UpdatedLagrangianUJWwPElement(UpdatedLagrangianUJWwPElement const& rOther);
 
 
    virtual ~UpdatedLagrangianUJWwPElement();
 
 
    UpdatedLagrangianUJWwPElement& operator=(UpdatedLagrangianUJWwPElement const& rOther);
 
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, PropertiesType::Pointer pProperties) const override;
 
    Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override;
 
    //************* STARTING - ENDING  METHODS
 
 
    void GetDofList(DofsVectorType& rElementalDofList, ProcessInfo& rCurrentProcessInfo) override;
 
    void EquationIdVector(EquationIdVectorType& rResult, ProcessInfo& rCurrentProcessInfo) override;
 
    void GetValuesVector(Vector& rValues, int Step = 0) const override;
 
    void GetFirstDerivativesVector(Vector& rValues, int Step = 0) const override;
 
    void GetSecondDerivativesVector(Vector& rValues, int Step = 0) const override;
 
 
    //************************************************************************************
    //************************************************************************************
    int Check(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateMassMatrix(MatrixType& rMassMatrix, 
            ProcessInfo& rCurrentProcessInfo) override;
 
    void CalculateDampingMatrix(MatrixType& rDampingMatrix, 
            ProcessInfo& rCurrentProcessInfo) override;
 
 
protected:
 
    
   /**** 
       the time step (requiered). It shall be somewhere else.
    ****/    
    double mTimeStep;
 
 
 
 
    virtual void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
                                    ElementDataType& rVariables,
                                    double& rIntegrationWeight) override;
 
    virtual void CalculateAndAddKWwP(MatrixType& rK,
                                     ElementDataType & rVariables,
                                     double& rIntegrationWeight
                                    );
 
    virtual void CalculateAndAddKUwP(MatrixType& rK,
                                     ElementDataType & rVariables,
                                     double& rIntegrationWeight
                                    );
    virtual void CalculateAndAddKPPStab(MatrixType& rK,
                                     ElementDataType & rVariables,
                                     double& rIntegrationWeight
                                     );
    virtual void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
                                    ElementDataType& rVariables,
                                    Vector& rVolumeForce,
                                    double& rIntegrationWeight) override;
 
    virtual void CalculateAndAddInternalWaterForces(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    void CalculateAndAddExternalForcesUJWwP(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          Vector & rVolumeForces,
          double& rIntegrationWeight
          );
    void CalculateAndAddFluidLinearMomentum(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
    void CalculateAndAddMassBalanceEquation(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    virtual void CalculateAndAddStabilizationRHS(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
    virtual void CalculateAndAddMassStabilizationMatrix(MatrixType& rMassMatrix,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    virtual void CalculateAndAddDampingStabilizationMatrix(MatrixType& rDampingMatrix,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    virtual void CalculateAndAddHighOrderDampingMatrix(MatrixType& rDampingMatrix,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    virtual void CalculateAndAddHighOrderKPP(MatrixType& rK,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
    virtual void CalculateAndAddHighOrderRHS(VectorType& rRightHandSideVector,
          ElementDataType & rVariables,
          double& rIntegrationWeight
          );
 
 
    virtual void InitializeElementData(ElementDataType & rVariables, const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
    void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
                                  VectorType& rRightHandSideVector,
                                  Flags& rCalculationFlags) override;
 
 
    double & CalculateStabilizationFactor( ElementDataType & rVariables, double & rStabilizationFactor);
 
 
private:
 
 
 
 
 
 
 
 
    friend class Serializer;
 
    // A private default constructor necessary for serialization
 
    virtual void save(Serializer& rSerializer) const override;
 
    virtual void load(Serializer& rSerializer) override;
 
 
 
 
}; // Class UpdatedLagrangianUJWwPElement
 
 
 
} // namespace Kratos
#endif // KRATOS_UPDATED_LAGRANGIAN_U_J_W_wP_ELEMENT_H_INCLUDED