axisym_updated_lagrangian_U_J_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosSolidMechanicsApplication $
//   Last modified by:    $Author:              LMonforte $
//   Date:                $Date:                July 2015 $
//   Revision:            $Revision:                  0.0 $
//
//
 
#if !defined(KRATOS_AXISYM_UPDATED_LAGRANGIAN_U_J_ELEMENT_H_INCLUDED )
#define  KRATOS_AXISYM_UPDATED_LAGRANGIAN_U_J_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/updated_lagrangian_U_J_element.hpp"
 
namespace Kratos
{
 
 
 
   class KRATOS_API(PFEM_SOLID_MECHANICS_APPLICATION) AxisymUpdatedLagrangianUJElement
      : public UpdatedLagrangianUJElement
   {
      public:
 
         typedef ConstitutiveLaw ConstitutiveLawType;
         typedef ConstitutiveLawType::Pointer ConstitutiveLawPointerType;
         typedef ConstitutiveLawType::StressMeasure StressMeasureType;
         typedef GeometryData::IntegrationMethod IntegrationMethod;
 
         KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( AxisymUpdatedLagrangianUJElement );
 
 
         AxisymUpdatedLagrangianUJElement();
 
         AxisymUpdatedLagrangianUJElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
         AxisymUpdatedLagrangianUJElement(IndexType NewId, GeometryType::Pointer pGeometry, PropertiesType::Pointer pProperties);
 
         AxisymUpdatedLagrangianUJElement(AxisymUpdatedLagrangianUJElement const& rOther);
 
 
         virtual ~AxisymUpdatedLagrangianUJElement();
 
 
         AxisymUpdatedLagrangianUJElement& operator=(AxisymUpdatedLagrangianUJElement 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 Initialize() override;
 
         //************* COMPUTING  METHODS
 
         void CalculateMassMatrix(MatrixType& rMassMatrix, 
               ProcessInfo& rCurrentProcessInfo) override;
 
 
         //************************************************************************************
         //************************************************************************************
         int Check(const ProcessInfo& rCurrentProcessInfo) override;
 
         virtual void CalculateKinematics(ElementDataType& rVariables,
               const double& rPointNumber) override;
 
 
         void CalculateDeformationGradient(const Matrix& rDN_DX,
               Matrix& rF,
               Matrix& rDeltaPosition, 
               double & rCurrentRadius, 
               double & rReferenceRadius);
 
         virtual void CalculateRadius(double & rCurrentRadius, double & rReferenceRadius, const Vector& rN);
 
 
         void CalculateDeformationMatrix(Matrix& rB,
                                            Matrix& rF,
                                            Vector& rN,
                                            double& rCurrentRadius);
 
         
         void CalculateGreenLagrangeStrain(const Matrix& rF,
               Vector& rStrainVector) override;
 
         void CalculateAlmansiStrain(const Matrix& rF,
               Vector& rStrainVector) 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 CalculateAndAddKuum(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKuug(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKuJ (MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKJu(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
 
 
         virtual void CalculateAndAddKJJ(MatrixType& rK,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddKJJStab(MatrixType& rK, ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
 
         virtual void CalculateAndAddJacobianForces(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
         virtual void CalculateAndAddStabilizedJacobian(VectorType& rRightHandSideVector,
               ElementDataType & rVariables,
               double& rIntegrationWeight
               ) override;
 
 
 
 
         virtual void GetHistoricalVariables( ElementDataType& rVariables, 
               const double& rPointNumber ) override;
 
 
         /*
          * Function to modify the deformation gradient to the constitutitve equation
          */
         virtual void ComputeConstitutiveVariables( ElementDataType& rVariables, Matrix& rFT, double& rdetFT) 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 AxisymUpdatedLagrangianUJElement
 
 
 
} // namespace Kratos
#endif // KRATOS_UPDATED_LAGRANGIAN_U_J_ELEMENT_H_INCLUDED