updated_lagrangian_U_J_P_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_U_J_P_ELEMENT_H_INCLUDED )
#define      KRATOS_UPDATED_LAGRANGIAN_U_J_P_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/updated_lagrangian_U_J_element.hpp"
 
#include "pfem_solid_mechanics_application_variables.h"
 
namespace Kratos
{
 
 
 
   class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) UpdatedLagrangianUJPElement
      : public UpdatedLagrangianUJElement
   {
 
      protected:
         typedef struct
         {
            unsigned int voigtsize;
 
            double Alpha;
            double Beta; 
 
            double NodalMeanStress;
            double ElementalMeanStress;
 
            double NodalJacobian;
 
            Vector StressVector;
            Vector StressVectorEC;
 
         } UJPElementData;
 
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( UpdatedLagrangianUJPElement );
 
 
         UpdatedLagrangianUJPElement();
 
         UpdatedLagrangianUJPElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         UpdatedLagrangianUJPElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         UpdatedLagrangianUJPElement(UpdatedLagrangianUJPElement const& rOther);
 
 
         virtual ~UpdatedLagrangianUJPElement();
 
 
         UpdatedLagrangianUJPElement& operator=(UpdatedLagrangianUJPElement 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
 
         //SET:
 
         //GET:
 
         void GetValueOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
         void GetValueOnIntegrationPoints(const Variable<Vector>& rVariable, std::vector<Vector>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
         void GetValueOnIntegrationPoints( const Variable<Matrix>& rVariable, std::vector<Matrix>& rValue, const ProcessInfo& rCurrentProcessInfo) 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;
 
 
 
      protected:
 
         double ScalingConstant1;
         double ScalingConstant2;
 
         //std::vector< Matrix > mDeformationGradientF0; // LMV: Crec que no cal
 
         //Vector mDeterminantF0; // LMV: Crec que no cal
 
 
 
 
         void CalculateAndAddLHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               double& rIntegrationWeight) override;
 
         void InitializeElementData(ElementDataType & rVariables, const ProcessInfo& rCurrentProcessInfo) override;
 
         void CalculateAndAddRHS(LocalSystemComponents& rLocalSystem,
               ElementDataType& rVariables,
               Vector& rVolumeForce,
               double& rIntegrationWeight) override;
 
 
         void CalculateAndAddKuumElemUJP(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddKuugElemUJP(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddKuJElemUJP (MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddKup (MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddKJuElemUJP(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddKJJElemUJP(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddKJp(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddKpu(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddKpJ(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddKpp(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
 
         void CalculateAndAddKppStab(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddKJJStabElemUJP(MatrixType& rK,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
 
         void CalculateAndAddJacobianForcesElemUJP( VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
 
 
         void CalculateAndAddPressureForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
 
         void CalculateAndAddStabilizedPressure(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
 
         void CalculateAndAddStabilizedJacobianElemUJP( VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddInternalForcesElemUJP(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               UJPElementData& rElementVariables,
               double& rIntegrationWeight
               );
 
         void InitializeSystemMatrices(MatrixType& rLeftHandSideMatrix,
               VectorType& rRightHandSideVector,
               Flags& rCalculationFlags) override;
 
         //on integration points:
         void CalculateOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
         void CalculateOnIntegrationPoints(const Variable<Vector>& rVariable, std::vector<Vector>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
         void CalculateOnIntegrationPoints(const Variable<Matrix>& rVariable, std::vector<Matrix>& rOutput, const ProcessInfo& rCurrentProcessInfo) override;
 
 
         /*
          * Compute Some integration point variables only once
          */ 
         void CalculateThisElementData( UJPElementData& rElementVariables, const ElementDataType& rVariables);
 
      private:
 
 
 
 
 
 
 
 
         friend class Serializer;
 
         // A private default constructor necessary for serialization
 
         void save(Serializer& rSerializer) const override;
 
         void load(Serializer& rSerializer) override;
 
 
 
 
}; // Class UpdatedLagrangianUJPElement
 
 
 
} // namespace Kratos
#endif // KRATOS_UPDATED_LAGRANGIAN_U_J_P_ELEMENT_H_INCLUDED