Steel1DFiber.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] Filippou, F. C., Popov, E. P., Bertero, V. V. (1983). "Effects of Bond 
//       Deterioration on Hysteretic Behavior of Reinforced Concrete Joints". 
//       Report EERC 83-19, Earthquake Engineering Research Center, University 
//       of California, Berkeley. 
//
// Description:
//------------------------------------------------------------------------------

#ifndef _STEEL1DFIBER_HPP_
#define _STEEL1DFIBER_HPP_

#include "Material.hpp"

class Steel1DFiber : public Material {

    public:
        Steel1DFiber(double fy, double E, double b, double R0=15.00, double cR1=0.925, double cR2=0.150, double a1=0.0, double a2=1.0, double a3=0.0, double a4=1.0, double nu=0.33, double rho=0.0);

        ~Steel1DFiber();

        //Clone the 'Steel1DFiber' material.
        std::unique_ptr<Material> CopyMaterial();

        double GetDensity() const;

        double GetPoissonRatio() const;

        double GetBulkModulus() const;

        double GetShearModulus() const;

        double GetElasticityModulus() const;

        double GetEnergy() const;

        Eigen::MatrixXd GetDamping() const;

        Eigen::VectorXd GetStrain() const;

        Eigen::VectorXd GetStress() const;

        Eigen::VectorXd GetStrainRate() const;

        Eigen::VectorXd GetTotalStress() const;

        Eigen::MatrixXd GetTangentStiffness() const;

        Eigen::MatrixXd GetInitialTangentStiffness() const;

        void CommitState();

        void ReverseState();

        void InitialState();

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

    private:
        //Steel properties:
        double nu;
        double fy; 
        double E;
        double b;
        double R0;
        double cR1;
        double cR2;
        double a1;
        double a2;
        double a3;
        double a4;
        double Rho;

        //Steel Commited history variables:
        double newStrain;
        double newStress;
        double newTangentStiffness;
        double newMinStrain;
        double newMaxStrain;
        double newPlasticStrain;
        double newStrainInterception;
        double newStressInterception;
        double newStrainLastInversion;
        double newStressLastInversion;
        int newLoadUnloadIndex;

        //Steel Updated history variables:
        double oldStrain;
        double oldStress;
        double oldTangentStiffness;
        double oldMinStrain;
        double oldMaxStrain;
        double oldPlasticStrain;
        double oldStrainInterception;
        double oldStressInterception;
        double oldStrainLastInversion;
        double oldStressLastInversion;
        int oldLoadUnloadIndex;
};

#endif