updated_lagrangian_U_J_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_ELEMENT_H_INCLUDED )
#define  KRATOS_UPDATED_LAGRANGIAN_U_J_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/solid_elements/large_displacement_element.hpp"
 
#include "pfem_solid_mechanics_application_variables.h"
 
namespace Kratos
{
 
 
 
   class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) UpdatedLagrangianUJElement
      : public LargeDisplacementElement
   {
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
         typedef GeometryData::SizeType SizeType;
         typedef LargeDisplacementElement::ElementDataType ElementDataType;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( UpdatedLagrangianUJElement );
 
 
         UpdatedLagrangianUJElement();
 
         UpdatedLagrangianUJElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         UpdatedLagrangianUJElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         UpdatedLagrangianUJElement(UpdatedLagrangianUJElement const& rOther);
 
 
         virtual ~UpdatedLagrangianUJElement();
 
 
         UpdatedLagrangianUJElement& operator=(UpdatedLagrangianUJElement 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
 
         void SetValuesOnIntegrationPoints(const Variable<double>& rVariable, std::vector<double>& rValues, const ProcessInfo& rCurrentProcessInfo) override;
 
 
         //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 Initialize() override;
 
         /*
          * Called at the beginning of each solution step
          */
 
         void InitializeSolutionStep( ProcessInfo& rCurrentProcessInfo) override;
 
         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;
 
         //************* COMPUTING  METHODS
 
         void CalculateMassMatrix(MatrixType& rMassMatrix, ProcessInfo& rCurrentProcessInfo) override;
 
         void CalculateDampingMatrix(MatrixType& rDampingMatrix, ProcessInfo& rCurrentProcessInfo) override;
 
 
         //************************************************************************************
         //************************************************************************************
         int Check(const ProcessInfo& rCurrentProcessInfo) override;
 
         void CalculateKinematics(ElementDataType& rVariables,
               const double& rPointNumber) override;
 
 
         void CalculateDeformationGradient(const Matrix& rDN_DX,
               Matrix& rF,
               Matrix& rDeltaPosition);
 
         void CalculateDeformationMatrix(Matrix& rB,
               Matrix& rF,
               Matrix& rDN_DX);
 
 
      protected:
 
 
         std::vector< Matrix > mDeformationGradientF0;
 
 
         Vector mDeterminantF0;
 
         /****
          * Some double member variable
          ****/
         double mElementStabilizationNumber;
 
 
 
         virtual SizeType GetDofsSize() const override;
 
         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;
 
         void FinalizeStepVariables(ElementDataType & rVariables, const double& rPointNumber) override;
 
 
 
         void CalculateAndAddKuum(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         void CalculateAndAddKuug(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKuJ (MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         virtual void CalculateAndAddKJu(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddKJJ(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddKJJStab(MatrixType& rK, ElementDataType & rVariables,
               double& rIntegrationWeight
               );
         void CalculateAndAddExternalForces(VectorType& rRightHandSideVector,
               ElementDataType& rVariables,
               Vector& rVolumeForce,
               double& rIntegrationWeight
               ) override;
 
 
         virtual void CalculateAndAddJacobianForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
 
         virtual void CalculateAndAddStabilizedJacobian(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               );
 
         void CalculateAndAddInternalForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) 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;
 
 
         void GetHistoricalVariables( ElementDataType& rVariables,
               const double& rPointNumber ) override;
 
         double& CalculateVolumeChange(double& rVolumeChange, ElementDataType& rVariables) override;
 
         /*
          * Function to modify the deformation gradient to the constitutitve equation
          */
         virtual void ComputeConstitutiveVariables( ElementDataType& rVariables, Matrix& rFT, double& rdetFT);
 
 
      private:
 
 
 
 
 
 
 
 
         friend class Serializer;
 
         // A private default constructor necessary for serialization
 
         void save(Serializer& rSerializer) const override;
 
         void load(Serializer& rSerializer) override;
 
 
 
 
}; // Class UpdatedLagrangianUJElement
 
 
 
} // namespace Kratos
#endif // KRATOS_UPDATED_LAGRANGIAN_U_J_ELEMENT_H_INCLUDED