Fib3DLineSection.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 _FIB3DLINESECTION_HPP_
#define _FIB3DLINESECTION_HPP_

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

#include "Material.hpp"
#include "Section.hpp"

class Fib3DLineSection : public Section{

    public:
        Fib3DLineSection(double h, double b, const std::vector<std::unique_ptr<Material> > &fibers, const std::vector<double> &zi, const std::vector<double> &yi, const std::vector<double> &Ai, double k2=1.0, double k3=1.0, double theta=0.0, unsigned int ip=10);

        ~Fib3DLineSection();

        std::unique_ptr<Section> CopySection();

        Eigen::VectorXd GetStrain();

        Eigen::VectorXd GetStress();

        Eigen::MatrixXd GetDensity();

        Eigen::MatrixXd GetTangentStiffness();

        Eigen::MatrixXd GetInitialTangentStiffness();

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

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

        void CommitState();

        void ReverseState();

        void InitialState();

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

    protected:
        void ComputeSectionCenter();

        unsigned int FiberIndexAt(double x3=0, double x2=0);

    private:
        double h;

        double b;

        double zcm;

        double ycm;

        double kappa2;

        double kappa3;

        double Theta;

        unsigned int InsertPoint;

        unsigned int NumberOfFibers;

        Eigen::VectorXd Strain;

        std::vector<double> zi;

        std::vector<double> yi;

        std::vector<double> Ai;

        std::vector<std::unique_ptr<Material> > theFiber;
};

#endif