large_displacement_beam_emc_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_LARGE_DISPLACEMENT_BEAM_EMC_ELEMENT_H_INCLUDED)
#define  KRATOS_LARGE_DISPLACEMENT_BEAM_EMC_ELEMENT_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "custom_elements/beam_elements/large_displacement_beam_element.hpp"
 
 
namespace Kratos
{
 
 
class KRATOS_API(SOLID_MECHANICS_APPLICATION) LargeDisplacementBeamEMCElement
    :public LargeDisplacementBeamElement
{
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 LargeDisplacementBeamElement::ElementDataType ElementDataType;
 
    KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION( LargeDisplacementBeamEMCElement );
 
 
    LargeDisplacementBeamEMCElement(IndexType NewId, GeometryType::Pointer pGeometry);
 
    LargeDisplacementBeamEMCElement(IndexType NewId, GeometryType::Pointer pGeometry,  PropertiesType::Pointer pProperties);
 
    LargeDisplacementBeamEMCElement(LargeDisplacementBeamEMCElement const& rOther);
 
    ~LargeDisplacementBeamEMCElement() 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 EMC Element #" << Id();
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "Large Displacement Beam EMC Element #" << Id();
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
      GetGeometry().PrintData(rOStream);
    }
 
protected:
 
 
 
    std::vector<Vector>  mCurrentStrainResultantsVector;
 
    std::vector<Vector>  mPreviousStrainResultantsVector;
 
 
    LargeDisplacementBeamEMCElement() {};
 
 
 
    void InitializeElementData(ElementDataType & rVariables,
                        const ProcessInfo& rCurrentProcessInfo) override;
 
 
    void MapToSpatialFrame(const ElementDataType& rVariables, Matrix& rVariable) override;
 
 
    void CalculateKinematics(ElementDataType& rVariables,
                                     const unsigned int& rPointNumber) override;
 
    Matrix& CalculatePreviousDeltaPosition(Matrix & rDeltaPosition);
 
 
    void CalculateFrameMapping(ElementDataType& rVariables,
                       const unsigned int& rPointNumber) override;
 
 
    void UpdateStrainVariables(ElementDataType& rVariables,
                       const unsigned int& rPointNumber) override;
 
 
    void CalculateAlphaRotationMatrix( const Matrix& rPreviousRotationMatrix,
                       const Matrix& rCurrentRotationMatrix,
                       Matrix& rAlphaRotationMatrix,
                       Matrix& rAlphaRotationMatrixAsterisk,
                       double Alpha);
 
 
    virtual void CalculateCurrentStrainResultantsVector(ElementDataType& rVariables,
                            Vector& rCurrentStrainResultantsVector,
                            double Alpha);
 
    virtual void CalculateCurrentCurvatureVector(ElementDataType& rVariables,
                         Vector& rCurrentCurvatureVector,
                         double Alpha);
 
    void CalculateConstitutiveMatrix(ElementDataType& rVariables) override;
 
    virtual void CalculateStrainResultants(Vector& rStrainResultants, ElementDataType& rVariables, double alpha);
 
    virtual void CalculateStrainCouples(Vector& rStrainCouples, ElementDataType& rVariables, double alpha);
 
 
    void CalculateStressResultants(ElementDataType& rVariables, const unsigned int& rPointNumber) 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 CalculateAndAddInertiaLHS(MatrixType& rLeftHandSideMatrix,
                       ElementDataType& rVariables,
                       const ProcessInfo& rCurrentProcessInfo,
                       double& rIntegrationWeight) override;
 
 
    void CalculateAndAddInertiaRHS(VectorType& rRightHandSideVector,
                       ElementDataType& rVariables,
                       const ProcessInfo& rCurrentProcessInfo,
                       double& rIntegrationWeight) override;
 
    void CalculateDifferentialOperator(MatrixType& rDifferentialOperator,
                           ElementDataType& rVariables,
                           const int& rNode,
                           double alpha) override;
 
    void CalculateRotationLinearPartTensor(Vector& rRotationVector, Matrix& rRotationTensor) override;
 
 
    void CalculateStrainEnergy(double& rEnergy, ElementDataType& rVariables, const ProcessInfo& rCurrentProcessInfo, double& rIntegrationWeight) override;
 
 
private:
 
    friend class Serializer;
 
 
    // A private default constructor necessary for serialization
 
 
    void save(Serializer& rSerializer) const override;
 
    void load(Serializer& rSerializer) override;
 
 
 
}; // Class LargeDisplacementBeamEMCElement
 
} // namespace Kratos.
#endif //  KRATOS_LARGE_DISPLACEMENT_BEAM_EMC_ELEMENT_H_INCLUDED defined