updated_lagrangian_U_Pressure_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosPfemSolidMechanicsApplication $
//   Created by:          $Author:                  LMonforte $
//   Last modified by:    $Co-Author:                         $
//   Date:                $Date:                    July 2015 $
//   Revision:            $Revision:                      0.0 $
//
//
 
#if !defined(KRATOS_UPDATED_LAGRANGIAN_UP_SECOND_ELEMENT_H_INCLUDED )
#define  KRATOS_UPDATED_LAGRANGIAN_UP_SECOND_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/solid_elements/updated_lagrangian_U_P_element.hpp"
 
namespace Kratos
{
 
 
   // THIS IS THE SECOND VERSION OF THE UP-ELEMENT FOR ANY CONSTITUTIVE EQUATION. 
   // In this version, the pressure is the pressure and it is not some "2D" invented pressure,..
   // the constitutive equation should not be the UP version
 
 
   class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) UpdatedLagrangianUPressureElement
      : public UpdatedLagrangianUPElement
   {
 
      protected:
         typedef struct
         {
            unsigned int voigtsize; 
 
            double NodalMeanStress;
            double ElementalMeanStress;
 
            Vector StressVector;
 
            Matrix DeviatoricTensor; 
            
         } ThisElementData;
      
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
         typedef GeometryData::SizeType SizeType;
         typedef UpdatedLagrangianUPElement::ElementDataType ElementDataType;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( UpdatedLagrangianUPressureElement );
 
 
         UpdatedLagrangianUPressureElement();
 
         UpdatedLagrangianUPressureElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         UpdatedLagrangianUPressureElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         UpdatedLagrangianUPressureElement(UpdatedLagrangianUPressureElement const& rOther);
 
 
         virtual ~UpdatedLagrangianUPressureElement();
 
 
         UpdatedLagrangianUPressureElement& operator=(UpdatedLagrangianUPressureElement 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;
 
 
 
 
         //************************************************************************************
         //************************************************************************************
         int Check(const ProcessInfo& rCurrentProcessInfo) override;
 
         /*
          * Get on rVariable a Matrix Value from the Element Constitutive Law
          */
         void GetValueOnIntegrationPoints(const Variable<Matrix>& rVariable, std::vector<Matrix>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
         
         
         // Get Values defined in order to avoid a clang warning (?)
         void GetValueOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
         virtual void CalculateOnIntegrationPoints(const Variable<Matrix >& rVariable, std::vector< Matrix >& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
 
 
 
      protected:
 
 
         double mElementStabilizationNumber;
 
         double mElementScalingNumber; 
 
 
 
         virtual void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               double& rIntegrationWeight) override;
 
         virtual void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               Vector& rVolumeForce,
               double& rIntegrationWeight) override;
 
         virtual void InitializeElementData(ElementDataType & rVariables, 
               const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
         virtual void CalculateAndAddKuumElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKuugElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKupElemUP (MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKpuElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddKppElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddKppStabElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         /*
          * Calculation of the Internal Forces due to sigma. Fi = B * sigma
          */
         void CalculateAndAddInternalForcesElemUP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddPressureForcesElemUP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddStabilizedPressureElemUP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateThisElementData( ThisElementData& rElementVariables, const ElementDataType & rVariables);
 
 
 
      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 UpdatedLagrangianUPressureElement
 
 
 
} // namespace Kratos
#endif // KRATOS_____