updated_lagrangian_U_P_wP_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosSolidMechanicsApplication $
//   Last modified by:    $Author:              LMonforte $
//   Date:                $Date:                July 2015 $
//   Revision:            $Revision:                 -0.1 $
//
//
 
#if !defined(KRATOS_UPDATED_LAGRANGIAN_U_P_wP_ELEMENT_H_INCLUDED )
#define  KRATOS_UPDATED_LAGRANGIAN_U_P_wP_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/updated_lagrangian_U_Pressure_element.hpp"
 
 
namespace Kratos
{
 
 
 
   class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) UpdatedLagrangianUPwPElement
      : public UpdatedLagrangianUPressureElement
   {
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( UpdatedLagrangianUPwPElement );
 
 
         UpdatedLagrangianUPwPElement();
 
         UpdatedLagrangianUPwPElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         UpdatedLagrangianUPwPElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         UpdatedLagrangianUPwPElement(UpdatedLagrangianUPwPElement const& rOther);
 
 
         virtual ~UpdatedLagrangianUPwPElement();
 
 
         UpdatedLagrangianUPwPElement& operator=(UpdatedLagrangianUPwPElement 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;
 
         //************* GETTING METHODS
 
 
         //************* 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;
 
         void FinalizeSolutionStep(ProcessInfo& rCurrentProcessInfo) override;
 
         //************************************************************************************
         //************************************************************************************
         int Check(const ProcessInfo& rCurrentProcessInfo) override;
 
 
      protected:
 
         double mTimeStep; // because I don't kwnow how to do it better 
 
 
         void CalculateElementalSystem(LocalSystemComponents& rLocalSystem,
               ProcessInfo& rCurrentProcessInfo) override;
 
 
         void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               double& rIntegrationWeight) override;
 
         void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               Vector& rVolumeForce,
               double& rIntegrationWeight) override;
 
         void InitializeElementData(ElementDataType & rVariables, const ProcessInfo& rCurrentProcessInfo) override;
 
         virtual void CalculateAndAddUnconsideredKuuTerms(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKuwP(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKwPu(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKwPP(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKwPwP(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKwPwPStab(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddExternalForces(VectorType& rRightHandSideVector,
               ElementDataType& rVariables,
               Vector& rVolumeForce,
               double& rIntegrationWeight
               ) override;
 
 
         void CalculateAndAddPressureForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
 
         void CalculateAndAddStabilizedPressure(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         void CalculateAndAddInternalForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddWaterPressureForces( VectorType& rRightHandSideVector,
               ElementDataType& rVariables,
               double& rIntegrationWeight
               );
         virtual void CalculateAndAddStabilizedWaterPressure( VectorType& rRightHandSideVector, 
               ElementDataType& rVariables,
               double& rIntegartionWeight
               );
 
         void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
               VectorType& rRightHandSideVector,
               Flags& rCalculationFlags) override;
 
 
         double& CalculateVolumeChange(double& rVolumeChange, ElementDataType& rVariables) override;
 
 
         void GetConstants( double& rScalingConstant, double& rWaterBulk, double& rDeltaTime, double& rPermeability);
 
         virtual double GetElementSize( const Matrix& rDN_DX);
 
 
      private:
 
 
 
 
 
 
 
 
         friend class Serializer;
 
         // A private default constructor necessary for serialization
 
         void save(Serializer& rSerializer) const override;
 
         void load(Serializer& rSerializer) override;
 
 
 
 
}; // Class UpdatedLagrangianUPwPElement
 
 
 
} // namespace Kratos
#endif // KRATOS_UPDATED_LAGRANGIAN_U_P_wP_ELEMENT_H_INCLUDED