axisym_updated_lagrangian_U_Pressure_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosSolidMechanicsApplication $
//   Last modified by:    $Author:              LMonforte $
//   Date:                $Date:                July 2013 $
//   Revision:            $Revision:                 -0.1 $
//
//
 
#if !defined(KRATOS_AXISYM_UPDATED_LAGRANGIAN_UPRESSURE_ELEMENT_H_INCLUDED )
#define  KRATOS_AXISYM_UPDATED_LAGRANGIAN_UPRESSURE_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/updated_lagrangian_U_Pressure_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) AxisymUpdatedLagrangianUPressureElement
      : public UpdatedLagrangianUPressureElement
   {
 
 
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( AxisymUpdatedLagrangianUPressureElement );
 
 
         AxisymUpdatedLagrangianUPressureElement();
 
         AxisymUpdatedLagrangianUPressureElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         AxisymUpdatedLagrangianUPressureElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         AxisymUpdatedLagrangianUPressureElement(AxisymUpdatedLagrangianUPressureElement const& rOther);
 
 
         virtual ~AxisymUpdatedLagrangianUPressureElement();
 
 
         AxisymUpdatedLagrangianUPressureElement& operator=(AxisymUpdatedLagrangianUPressureElement 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;
 
 
         virtual void Initialize() override;
 
 
 
 
 
      protected:
 
 
 
 
 
         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
               ) override;
 
         virtual void CalculateAndAddKuugElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKupElemUP (MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKpuElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
 
         virtual void CalculateAndAddKppElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
 
         virtual void CalculateAndAddKppStabElemUP(MatrixType& rK,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddPressureForcesElemUP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
 
         virtual void CalculateAndAddStabilizedPressureElemUP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               ThisElementData & rElementVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateThisElementData( ThisElementData& rElementVariables, const ElementDataType & rVariables) override;
 
 
 
         virtual void CalculateKinematics(ElementDataType& rVariables,
               const double& rPointNumber) override;
 
 
         void CalculateRadius(double & rCurrentRadius,
               double & rReferenceRadius,
               const Vector& rN);
 
         void CalculateDeformationGradient(const Matrix& rDN_DX,
               Matrix& rF,
               Matrix& rDeltaPosition,
               double & rCurrentRadius,
               double & rReferenceRadius);
 
         void CalculateDeformationMatrix(Matrix& rB,
               Matrix& rDN_DX,
               Vector& rN,
               double & rCurrentRadius);
 
 
         void CalculateGreenLagrangeStrain(const Matrix& rF,
               Vector& rStrainVector) override;
 
         void CalculateAlmansiStrain(const Matrix& rF,
               Vector& rStrainVector) override;
 
 
 
 
      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 AxisymUpdatedLagrangianUPressureElement
 
 
 
} // namespace Kratos
#endif // KRATOS_____