large_displacement_beam_element.hpp

Go to the documentation of this file.

//
//   Project Name:        KratosSolidMechanicsApplication $
//   Created by:          $Author:            JMCarbonell $
//   Last modified by:    $Co-Author:                     $
//   Date:                $Date:            November 2015 $
//   Revision:            $Revision:                  0.0 $
//
//
 
#if !defined(KRATOS_LARGE_DISPLACEMENT_BEAM_ELEMENT_H_INCLUDED)
#define  KRATOS_LARGE_DISPLACEMENT_BEAM_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/beam_elements/beam_element.hpp"
 
namespace Kratos
{
 
 
class KRATOS_API(SOLID_MECHANICS_APPLICATION) LargeDisplacementBeamElement
    :public BeamElement
{
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 BeamElement::ElementDataType                ElementDataType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( LargeDisplacementBeamElement );
 
 
    LargeDisplacementBeamElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
    LargeDisplacementBeamElement(IndexType NewId, GeometryType::Pointer pGeometry,  PropertiesType::Pointer pProperties);
 
    LargeDisplacementBeamElement(LargeDisplacementBeamElement const& rOther);
 
    ~LargeDisplacementBeamElement() override;
 
 
 
 
    Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes,  PropertiesType::Pointer pProperties) const override;
 
 
    void Initialize(const ProcessInfo& rCurrentProcessInfo) override;
 
    void InitializeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
    void InitializeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
    void FinalizeNonLinearIteration(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void FinalizeSolutionStep(const ProcessInfo& rCurrentProcessInfo) override;
 
 
    //************* COMPUTING  METHODS
 
 
    void CalculateMassMatrix(MatrixType& rMassMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
 
    //on integration points:
    void CalculateOnIntegrationPoints(const Variable<double>& rVariable,
                      std::vector<double>& rOutput,
                      const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void CalculateOnIntegrationPoints( const Variable< array_1d<double, 3 > >& rVariable,
                                       std::vector< array_1d<double, 3 > >& Output,
                                       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:
 
 
    IntegrationMethod mReducedIntegrationMethod;
 
    IntegrationMethod mFullIntegrationMethod;
 
    int  mIterationCounter;
 
    double mInvJ0;
 
    std::vector<Vector>  mCurrentCurvatureVectors;
 
    std::vector<Vector>  mPreviousCurvatureVectors;
 
    std::vector<QuaternionType>  mFrameQuaternionsReduced;
 
    std::vector<QuaternionType>  mFrameQuaternionsFull;
 
 
    LargeDisplacementBeamElement() {};
 
 
    IntegrationMethod GetReducedIntegrationMethod() const;
 
    IntegrationMethod GetFullIntegrationMethod() const;
 
    void CalculateDynamicSystem( LocalSystemComponents& rLocalSystem,
                 const ProcessInfo& rCurrentProcessInfo ) override;
 
    virtual void MapToSpatialFrame(const ElementDataType& rVariables, Matrix& rVariable);
 
 
    virtual void CalculateCurrentCurvature(ElementDataType& rVariables, const Variable<array_1d<double, 3 > >& rVariable);
 
 
    void CalculateKinematics(ElementDataType& rVariables,
                                     const unsigned int& rPointNumber) override;
 
    virtual void CalculateFrameMapping(ElementDataType& rVariables,
                       const unsigned int& rPointNumber);
 
 
    virtual void UpdateStrainVariables(ElementDataType& rVariables,
                       const unsigned int& rPointNumber);
 
 
    void CalculateConstitutiveMatrix(ElementDataType& rVariables) 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;
 
    virtual void CalculateAndAddKuug2(MatrixType& rLeftHandSideMatrix,
                                     ElementDataType& rVariables,
                                     double& rIntegrationWeight);
 
    virtual void CalculateAndAddKuuf(MatrixType& rLeftHandSideMatrix,
                                     ElementDataType& rVariables,
                                     double& rIntegrationWeight);
 
 
 
    virtual void CalculateAndAddFollowerForces(VectorType& rRightHandSideVector,
                           ElementDataType& rVariables,
                           double& rIntegrationWeight);
 
 
    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 CalculateOperator(MatrixType& rOperator,
               Vector& rN,
               const int& rNode);
 
    void CalculateDiscreteOperator(MatrixType& rDiscreteOperator,
                   ElementDataType& rVariables,
                   const int& rNode);
 
    virtual void CalculateDifferentialOperator(MatrixType& rDifferentialOperator,
                           ElementDataType& rVariables,
                           const int& rNode,
                           double alpha = 0);
 
 
    void CalculateBmatrix(MatrixType& rBmatrix,
              ElementDataType& rVariables);
 
    double& CalculateTotalMass( SectionProperties& Section, double& rTotalMass );
 
    void CalculateInertiaDyadic(SectionProperties& rSection, Matrix& rInertiaDyadic);
 
    virtual void CalculateRotationLinearPartTensor(Vector& rRotationVector, Matrix& rRotationTensor);
 
 
    void GetCurrentNodalMovements(Vector& rValues, const int& rNode, unsigned int PointNumber = 0);
 
    void GetCurrentNodalVelocities(Vector& rValues, const int& rNode, unsigned int PointNumber = 0);
 
    void GetKineticMatrix(Matrix& rKineticMatrix, const double& rMass, const Matrix& rInertia);
 
 
    virtual void CalculateExternalForcesEnergy(double& rEnergy, ElementDataType& rVariables, Vector& rVolumeForce, double& rIntegrationWeight);
 
    virtual void CalculateInternalForcesEnergy(double& rEnergy, ElementDataType& rVariables, double& rIntegrationWeight);
 
 
    virtual void CalculateStrainEnergy(double& rEnergy, ElementDataType& rVariables, const ProcessInfo& rCurrentProcessInfo, double& rIntegrationWeight);
 
    virtual void CalculateKineticEnergy(double& rEnergy, ElementDataType& rVariables, const ProcessInfo& rCurrentProcessInfo, double& rIntegrationWeight);
 
    virtual void CalculateMomentumRelations(array_1d<double,3>& LinearMomentum, array_1d<double,3>& AngularMomentum, ElementDataType& rVariables, const ProcessInfo& rCurrentProcessInfo, double& rIntegrationWeight);
 
 
private:
 
    friend class Serializer;
 
 
    // A private default constructor necessary for serialization
 
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
}; // Class LargeDisplacementBeamElement
 
} // namespace Kratos.
#endif //  KRATOS_LARGE_DISPLACEMENT_BEAM_ELEMENT_H_INCLUDED defined