geometrically_exact_rod_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosSolidMechanicsApplication $
//   Created by:          $Author:            JMCarbonell $
//   Last modified by:    $Co-Author:                     $
//   Date:                $Date:              August 2017 $
//   Revision:            $Revision:                  0.0 $
//
//
 
#if !defined(KRATOS_GEOMETRICALLY_EXACT_ROD_ELEMENT_H_INCLUDED)
#define  KRATOS_GEOMETRICALLY_EXACT_ROD_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/beam_elements/large_displacement_beam_emc_element.hpp"
 
namespace Kratos
{
 
 
class KRATOS_API(SOLID_MECHANICS_APPLICATION) GeometricallyExactRodElement
    :public LargeDisplacementBeamEMCElement
{
public:
 
    typedef ConstitutiveLaw                         ConstitutiveLawType;
    typedef ConstitutiveLawType::Pointer     ConstitutiveLawPointerType;
    typedef ConstitutiveLawType::StressMeasure        StressMeasureType;
    typedef GeometryData::IntegrationMethod           IntegrationMethod;
    typedef BeamMathUtils<double>                     BeamMathUtilsType;
    typedef Quaternion<double>                           QuaternionType;
    typedef GeometryData::SizeType                             SizeType;
    typedef LargeDisplacementBeamEMCElement::ElementDataType ElementDataType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( GeometricallyExactRodElement );
 
 
    GeometricallyExactRodElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
    GeometricallyExactRodElement(IndexType NewId, GeometryType::Pointer pGeometry,  PropertiesType::Pointer pProperties);
 
    GeometricallyExactRodElement(GeometricallyExactRodElement const& rOther);
 
    ~GeometricallyExactRodElement() override;
 
 
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes,  PropertiesType::Pointer pProperties) const override;
 
 
    //************* STARTING - ENDING  METHODS
 
    void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    //************************************************************************************
    //************************************************************************************
    int Check(const ProcessInfo& rCurrentProcessInfo) const override;
 
 
    std::string Info() const override
    {
      std::stringstream buffer;
      buffer << "Large Displacement Beam Element #" << Id();
      return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
      rOStream << "Large Displacement Beam Element #" << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
      GetGeometry().PrintData(rOStream);
    }
 
protected:
 
 
    std::vector<Matrix> mInitialLocalDirectors;
 
    std::vector<Matrix> mCurrentLocalDirectors;
 
    std::vector<Matrix> mPreviousLocalDirectors;
 
    std::vector<Matrix> mInitialLocalDirectorsVelocities;
 
    std::vector<Matrix> mCurrentLocalDirectorsVelocities;
 
    std::vector<Matrix> mPreviousLocalDirectorsVelocities;
 
 
    GeometricallyExactRodElement() {};
 
 
    void CalculateElementalSystem( LocalSystemComponents& rLocalSystem,
                   const ProcessInfo& rCurrentProcessInfo ) override;
 
 
    void CalculateDynamicSystem( LocalSystemComponents& rLocalSystem,
                 const ProcessInfo& rCurrentProcessInfo ) override;
 
 
 
    void InitializeElementData(ElementDataType & rVariables, const ProcessInfo& rCurrentProcessInfo) override;
 
 
 
    void CalculateKinematics(ElementDataType& rVariables,
                                     const unsigned int& rPointNumber) override;
 
    void CalculateConstitutiveMatrix(ElementDataType& rVariables) override;
 
    void CalculateFrameMapping(ElementDataType& rVariables,
                       const unsigned int& rPointNumber) override;
 
 
    void UpdateRotationVariables(ElementDataType& rVariables,
                 const unsigned int& rPointNumber);
 
    virtual void CalculateDirectorsMappingTensor(Matrix& rMappingTensor,
                         ElementDataType& rVariables,
                         const int& rNode,
                         double alpha);
 
 
    void CalculateCurrentStrainResultantsVector(ElementDataType& rVariables,
                            Vector& rCurrentStrainResultantsVector,
                            double Alpha) override;
 
    void CalculateCurrentCurvatureVector(ElementDataType& rVariables,
                         Vector& rCurrentCurvatureVector,
                         double Alpha) override;
 
 
    void CalculateStressResultants(ElementDataType& rVariables, const unsigned int& rPointNumber) override;
 
 
    void CalculateAndAddKuum(MatrixType& rLeftHandSideMatrix,
                 ElementDataType& rVariables,
                 double& rIntegrationWeight) override;
 
 
 
    void CalculateAndAddKuug(MatrixType& rLeftHandSideMatrix,
                 ElementDataType& rVariables,
                 double& rIntegrationWeight) override;
 
 
    void CalculateAndAddKuuf(MatrixType& rLeftHandSideMatrix,
                 ElementDataType& rVariables,
                 double& rIntegrationWeight) override;
 
 
 
    void CalculateAndAddFollowerForces(VectorType& rRightHandSideVector,
                       ElementDataType& rVariables,
                       double& rIntegrationWeight) override;
 
 
    void CalculateAndAddExternalForces(VectorType& rRightHandSideVector,
                       ElementDataType& rVariables,
                       Vector& rVolumeForce,
                       double& rIntegrationWeight) override;
 
    void CalculateAndAddInertiaLHS(MatrixType& rLeftHandSideMatrix,
                   ElementDataType& rVariables,
                   const ProcessInfo& rCurrentProcessInfo,
                   double& rIntegrationWeight) override;
 
 
    void CalculateAndAddInertiaRHS(VectorType& rRightHandSideVector,
                   ElementDataType& rVariables,
                   const ProcessInfo& rCurrentProcessInfo,
                   double& rIntegrationWeight) override;
 
 
    void CalculateAndAddInternalForces(VectorType& rRightHandSideVector,
                       ElementDataType & rVariables,
                       double& rIntegrationWeight) override;
 
 
    void CalculateDiscreteOperatorN(MatrixType& rDiscreteOperator,
                    ElementDataType& rVariables,
                    const int& rNodeI,
                    const int& rNodeJ,
                    const int& rComponent);
 
    void CalculateDiscreteOperatorM(MatrixType& rDiscreteOperator,
                    ElementDataType& rVariables,
                    const int& rNodeI,
                    const int& rNodeJ,
                    const int& rComponent);
 
    void CalculateDifferentialOperator(MatrixType& rDifferentialOperator,
                       ElementDataType& rVariables,
                       const int& rNode,
                       double alpha) override;
 
 
     void CalculateAlphaDirectors(Matrix& rDirectors, ElementDataType& rVariables, const int& rNode, double alpha);
 
 
     void CalculateAlphaDirectorVector(Vector& rDirectorVector, ElementDataType& rVariables, const int& rNode, const int& rDirection, double alpha);
 
     void CalculateAlphaDirectorSkewSymTensor(Matrix& rDirectorSkewSymTensor, ElementDataType& rVariables, const int& rNode, const int& rDirection, double alpha);
 
 
    void CalculateAlphaDirectorVector(Vector& rDirectorVector, ElementDataType& rVariables, const int& rDirection, double alpha);
 
    void CalculateAlphaDirectorSkewSymTensor(Matrix& rDirectorSkewSymTensor, ElementDataType& rVariables, const int& rDirection, double alpha);
 
    void CalculateDirectorDerivativesVector(Vector& rDirectorDerivativesVector, ElementDataType& rVariables, const int& rDirection, double alpha);
 
    void CalculateDirectorDerivativesSkewSymTensor(Matrix& rDirectorDerivativesSkewSymTensor, ElementDataType& rVariables, const int& rDirection, double alpha);
 
 
    void CalculateAlgorithmicInertia(Matrix & rAlgorithmicInertia,
                     const Matrix& rInertiaDyadic,
                     ElementDataType & rVariables,
                     const int& NodeJ,
                     const int& NodeI,
                     double alpha);
 
 
private:
 
    friend class Serializer;
 
 
    // A private default constructor necessary for serialization
 
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
 
}; // Class GeometricallyExactRodElement
 
} // namespace Kratos.
#endif //  KRATOS_GEOMETRICALLY_EXACT_ROD_ELEMENT_H_INCLUDED defined