Element.hpp

Go to the documentation of this file.

//==============================================================================
//
//                       Seismo Virtual Laboratory
//             Module for Serial and Parallel Analysis of seismic 
//         wave propagation and soil-structure interaction simulation
//         Copyright (C) 2018-2021, The California Institute of Technology
//                         All Rights Reserved.
//
// Commercial use of this program without express permission of the California
// Institute of Technology, is strictly  prohibited. See  file "COPYRIGHT"  in
// main  directory  for  information on  usage  and  redistribution, and for a
// DISCLAIMER OF ALL WARRANTIES.
//
//==============================================================================
//
// Written by:
//   Danilo S. Kusanovic (dkusanov@caltech.edu)
//   Elnaz E. Seylabi    (elnaze@unr.edu)
//
// Supervised by:
//   Domniki M. Asimaki  (domniki@caltech.edu)
//
// References : 
//   [1]
//
// Description:
//------------------------------------------------------------------------------

#ifndef _ELEMENT_HPP_
#define _ELEMENT_HPP_

#include <map>
#include <vector>
#include <memory>
#include <string>
#include <Eigen/Dense>

#include "Node.hpp"
#include "Load.hpp"
#include "Damping.hpp"

class Element{

    public:
        Element(std::string name, const std::vector<unsigned int> nodes, unsigned int ndofs, unsigned int VTKcell, unsigned int SVLcell);

        virtual ~Element() = 0;

        virtual void CommitState() = 0;

        virtual void UpdateState() = 0;

        virtual void ReverseState() = 0;

        virtual void InitialState() = 0;

        virtual void SetDomain(std::map<unsigned int, std::shared_ptr<Node> > &nodes) = 0;

        virtual void SetDamping(const std::shared_ptr<Damping> &damping) = 0;

        virtual std::vector<unsigned int> GetTotalDegreeOfFreedom() const = 0; 

        virtual Eigen::MatrixXd GetStrain() const = 0;

        virtual Eigen::MatrixXd GetStress() const = 0;

        virtual Eigen::MatrixXd GetStrainRate() const = 0;

        virtual Eigen::MatrixXd GetStrainAt(double x3=0.0, double x2=0.0) const = 0;

        virtual Eigen::MatrixXd GetStressAt(double x3=0.0, double x2=0.0) const = 0;

        virtual Eigen::VectorXd GetVTKResponse(std::string response) const = 0;

        virtual double ComputeEnergy() = 0;

        virtual Eigen::MatrixXd ComputeMassMatrix() = 0;

        virtual Eigen::MatrixXd ComputeStiffnessMatrix() = 0;

        virtual Eigen::MatrixXd ComputeDampingMatrix() = 0;

        virtual Eigen::MatrixXd ComputePMLMatrix() = 0;

        virtual Eigen::VectorXd ComputeInternalForces() = 0;

        virtual Eigen::VectorXd ComputeInternalDynamicForces() = 0;

        virtual Eigen::VectorXd ComputeSurfaceForces(const std::shared_ptr<Load> &surface, unsigned int k) = 0;

        virtual Eigen::VectorXd ComputeBodyForces(const std::shared_ptr<Load> &body, unsigned int k=0) = 0;

        virtual Eigen::VectorXd ComputeDomainReductionForces(const std::shared_ptr<Load> &drm, unsigned int k) = 0;

        std::string GetName() const;

        unsigned int GetVTKCellType() const;

        unsigned int GetVTKGroupType() const;

        unsigned int GetNumberOfNodes() const;

        unsigned int GetNumberOfDegreeOfFreedom() const;

        const std::vector<unsigned int>& GetNodes() const;

        bool HasFixedNode(const std::vector<std::shared_ptr<Node> > &nodes) const;

    private:
        std::string Name;

        unsigned int VTKCell;

        unsigned int SVLCell;

        unsigned int NumberOfNodes;

        unsigned int NumberOfDegreeOfFreedom;

        std::vector<unsigned int> Nodes;
};

#endif