empirical_spring.hpp

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors: Klaus B. Sautter
//
//
//
 
#if !defined(KRATOS_EMPIRICAL_SPRING_H_INCLUDED )
#define  KRATOS_EMPIRICAL_SPRING_H_INCLUDED
 
// System includes
 
// External includes
 
// Project includes
#include "includes/element.h"
#include "includes/define.h"
 
namespace Kratos
{
    class KRATOS_API(CABLE_NET_APPLICATION) EmpiricalSpringElement3D2N : public Element
    {
    protected:
        //const values
        static constexpr int msNumberOfNodes = 2;
        static constexpr int msDimension = 3;
        static constexpr unsigned int msLocalSize = msNumberOfNodes * msDimension;
 
    public:
        KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(EmpiricalSpringElement3D2N);
 
 
        typedef Element BaseType;
        typedef BaseType::GeometryType GeometryType;
        typedef BaseType::NodesArrayType NodesArrayType;
        typedef BaseType::PropertiesType PropertiesType;
        typedef BaseType::IndexType IndexType;
        typedef BaseType::SizeType SizeType;
        typedef BaseType::MatrixType MatrixType;
        typedef BaseType::VectorType VectorType;
        typedef BaseType::EquationIdVectorType EquationIdVectorType;
        typedef BaseType::DofsVectorType DofsVectorType;
 
 
        EmpiricalSpringElement3D2N() {};
        EmpiricalSpringElement3D2N(IndexType NewId,
                        GeometryType::Pointer pGeometry);
        EmpiricalSpringElement3D2N(IndexType NewId,
                        GeometryType::Pointer pGeometry,
                        PropertiesType::Pointer pProperties);
 
 
        ~EmpiricalSpringElement3D2N() override;
 
    Element::Pointer Create(
        IndexType NewId,
        GeometryType::Pointer pGeom,
        PropertiesType::Pointer pProperties
        ) const override;
 
    Element::Pointer Create(
        IndexType NewId,
        NodesArrayType const& ThisNodes,
        PropertiesType::Pointer pProperties
        ) const override;
 
        void EquationIdVector(
            EquationIdVectorType& rResult,
            const ProcessInfo& rCurrentProcessInfo) const override;
 
        void GetDofList(
            DofsVectorType& rElementalDofList,
            const ProcessInfo& rCurrentProcessInfo) const override;
 
        virtual BoundedMatrix<double,msLocalSize,msLocalSize>
         CreateElementStiffnessMatrix(const ProcessInfo& rCurrentProcessInfo);
 
 
        void CalculateLocalSystem(
            MatrixType& rLeftHandSideMatrix,
            VectorType& rRightHandSideVector,
            const ProcessInfo& rCurrentProcessInfo) override;
 
 
        void CalculateRightHandSide(
            VectorType& rRightHandSideVector,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        void CalculateLeftHandSide(
            MatrixType& rLeftHandSideMatrix,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        void AddExplicitContribution(
            const VectorType& rRHSVector, const Variable<VectorType>& rRHSVariable,
            const Variable<array_1d<double, 3>>& rDestinationVariable,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        void AddExplicitContribution(
            const VectorType& rRHSVector,
            const Variable<VectorType>& rRHSVariable,
            const Variable<double >& rDestinationVariable,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        void CalculateMassMatrix(
            MatrixType& rMassMatrix,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        void CalculateLumpedMassVector(
            VectorType &rLumpedMassVector,
            const ProcessInfo &rCurrentProcessInfo) const override;
 
        void CalculateDampingMatrix(MatrixType& rDampingMatrix,
            const ProcessInfo& rCurrentProcessInfo) override;
 
        double EvaluatePolynomial(const Vector& rPolynomial) const;
        double EvaluatePolynomialFirstDerivative(const Vector& rPolynomial) const;
 
        void LocalizeVector(BoundedVector<double,msLocalSize>& rInputVector);
        void GlobalizeVector(BoundedVector<double,msLocalSize>& rInputVector);
        void GlobalizeMatrix(BoundedMatrix<double,msLocalSize,msLocalSize>& rInputMatrix);
 
        void WriteTransformationCoordinates(BoundedVector<double,msLocalSize>
            & rReferenceCoordinates);
 
        void CreateTransformationMatrix(BoundedMatrix<double,msLocalSize,
            msLocalSize>& rRotationMatrix);
 
        double GetElementElongation() const;
 
 
        void GetValuesVector(
            Vector& rValues,
            int Step = 0) const override;
 
        void GetSecondDerivativesVector(
            Vector& rValues,
            int Step = 0) const override;
 
        void GetFirstDerivativesVector(
            Vector& rValues,
            int Step = 0) const override;
 
        int  Check(
            const ProcessInfo& rCurrentProcessInfo) const override;
 
private:
 
    friend class Serializer;
    void save(Serializer& rSerializer) const override;
    void load(Serializer& rSerializer) override;
    };
 
 
}
 
 
#endif