Plastic3DBA.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] Multiaxial Cyclic Plasticity Model for Clays, Ronaldo I. Borja and Alexander 
//      P. Amies, Journal of Geotechnical Engineering, Vol. 120, No. 6, June, 1994
//
// Description:
//------------------------------------------------------------------------------

#ifndef _PLASTIC3DBA_HPP_
#define _PLASTIC3DBA_HPP_

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

#include "Material.hpp"

class Plastic3DBA : public Material{

    public:
        Plastic3DBA(const double K, const double G, const double rho, const double H0, const double h, const double m, const double Su, const double beta);

        ~Plastic3DBA();

        //Clone the 'Plastic3DBA' 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);

    protected:
        Eigen::VectorXd ComputeIdentityVector() const;

        Eigen::MatrixXd ComputeDeviatoricTensor() const;

        Eigen::MatrixXd ComputeIdentityTensor() const;

        Eigen::MatrixXd ComputeIdentityOperator() const;

        double ComputeTensorTrace(const Eigen::VectorXd &T);

        double ComputeTensorNorm(const Eigen::VectorXd &T);

        double ComputeInnerProduct(const Eigen::VectorXd &V1, const Eigen::VectorXd &V2);

        Eigen::VectorXd ComputeHardening(Eigen::VectorXd const &DeviatoricIncrStrain, const Eigen::VectorXd &DeviatoricStress);

        Eigen::VectorXd ComputeHardeningBisection(Eigen::VectorXd const &DeviatoricIncrStrain, const Eigen::VectorXd &DeviatoricStress);

    private:
        double K;

        double G;

        double Rho;
    
        double H0;

        double h;
    
        double m;

        double Su;

        double R;

        double beta;

        double psi;

        double kappa;

        double psi_n;

        double kappa_n;
    
        Eigen::VectorXd Strain;

        Eigen::VectorXd Strain_n;

        Eigen::VectorXd Stress;

        Eigen::VectorXd Stress_n;

        Eigen::VectorXd DeviatoricStress0;

        Eigen::VectorXd DeviatoricStress0_n;

        Eigen::MatrixXd TangentStiffness;

        Eigen::MatrixXd TangentStiffness_n;

        int FirstLoadFlag;

        int FirstLoadFlag_n;
        
        int rootFlag;

        int rootFlag_n;
        
        double Hn;
};

#endif