Section.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 _SECTION_HPP_
#define _SECTION_HPP_

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

class Section {

    public:
        Section(std::string name);

        virtual ~Section() = 0;

        virtual std::unique_ptr<Section> CopySection() = 0;

        virtual Eigen::VectorXd GetStrain() = 0;

        virtual Eigen::VectorXd GetStress() = 0;

        virtual Eigen::MatrixXd GetDensity() = 0;

        virtual Eigen::MatrixXd GetTangentStiffness() = 0;

        virtual Eigen::MatrixXd GetInitialTangentStiffness() = 0;

        virtual Eigen::VectorXd GetStrainAt(double x3=0, double x2=0) = 0;

        virtual Eigen::VectorXd GetStressAt(double x3=0, double x2=0) = 0;

        virtual void CommitState() = 0;

        virtual void ReverseState() = 0;

        virtual void InitialState() = 0;

        virtual void UpdateState(Eigen::VectorXd strain, unsigned int cond=0) = 0;

        std::string GetName();

    protected:
        Eigen::MatrixXd GetLineRotationMatrix(double theta);

        Eigen::MatrixXd GetAreaRotationMatrix(double Theta);

        Eigen::MatrixXd GetLineTranslationMatrix(double h, double b, double zc, double yc, unsigned int ip);

        void InsertionPointCoordinates(double &x3, double &x2, double h, double b, double zc, double yc, unsigned int ip);

        Eigen::MatrixXd ComputeLineLocalAxes(double h, double b, double zcm, double ycm, double angle, unsigned int ip);

    private:
        std::string Name;
};

#endif