PressureDependMultiYield.h

// -*-c++-*-
//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis C. Pérez Tato
//
//  This program derives from OpenSees <http://opensees.berkeley.edu>
//  developed by the  «Pacific earthquake engineering research center».
//
//  Except for the restrictions that may arise from the copyright
//  of the original program (see copyright_opensees.txt)
//  XC is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or 
//  (at your option) any later version.
//
//  This software is distributed in the hope that it will be useful, but 
//  WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details. 
//
//
// You should have received a copy of the GNU General Public License 
// along with this program.
// If not, see <http://www.gnu.org/licenses/>.
//----------------------------------------------------------------------------
//<<<<<<< PressureDependMultiYield.h
// $Revision: 1.19 $
// $Date: 2004/06/15 18:58:01 $
//=======
// $Revision: 1.19 $
// $Date: 2004/06/15 18:58:01 $
//>>>>>>> 1.7
// $Source: /usr/local/cvs/OpenSees/SRC/material/nD/soil/PressureDependMultiYield.h,v $
                                                                        
// Written: ZHY
// Created: August 2000


// Description: This file contains the class prototype for PressureDependMultiYield.
//
// What: "@(#) PressureDependMultiYield.h, revA"

#ifndef PressureDependMultiYield_h
#define PressureDependMultiYield_h

#include "material/nD/soil/PressureDependMultiYieldBase.h"

namespace XC {

class MultiYieldSurface;

//
class PressureDependMultiYield: public PressureDependMultiYieldBase
  {
  private:
  // user supplied 
     static std::vector<double> liquefyParam4x;

     mutable T2Vector reversalStress;
     mutable T2Vector lockStress;

     T2Vector reversalStressCommitted;
     T2Vector lockStressCommitted;
     
     // Called by constructor
     void setUpSurfaces(double *);  
     void initStrainUpdate(void);

     // Return num_strain_subincre
     int isCriticalState(const T2Vector & stress) const;
     void updatePPZ(const T2Vector & stress) const;
     void PPZTranslation(const T2Vector & contactStress) const;
     double getPPZLimits(int which, const T2Vector & contactStress) const;
     double getPlasticPotential(const T2Vector & stress, const T2Vector & surfaceNormal) const;
     //return 1 if stress locked; o/w return 0.
     int stressCorrection(int crossedSurface) const;

  protected:
    int sendData(Communicator &);
    int recvData(const Communicator &);
  public:
     // Initialization constructor
     PressureDependMultiYield(int tag, 
                   int nd,
                         double rho,
                   double refShearModul,
                   double refBulkModul,
                   double frictionAng,
                   double peakShearStra,
                   double refPress,
                   double pressDependCoe,
                   double phaseTransformAngle, 
                   double contractionParam1,
                   double dilationParam1,
                   double dilationParam2,
                   double liquefactionParam1,
                   double liquefactionParam2,
                   double liquefactionParam4,
                   int   numberOfYieldSurf = 20,
                   double * gredu = 0,
                   double e = 0.6,
                   double volLimit1 = 0.9,
                   double volLimit2 = 0.02,
                   double volLimit3 = 0.7,
                   double atm = 101.,
                   double cohesi = 0.1,
                   double hv = 0.,
                   double pv = 1.);

     PressureDependMultiYield(int tag);

     // Default constructor
     PressureDependMultiYield();

     // Copy constructor
     PressureDependMultiYield(const PressureDependMultiYield &);

     // Calculates current tangent stiffness.
     const Matrix &getTangent(void) const;

     // Calculates the corresponding stress increment(rate), for a given strain increment. 
     const Vector &getStress(void) const;

     // Return an exact copy of itself.
     NDMaterial *getCopy(void) const;

     // Return a copy of itself if "code"="PressureDependMultiYield", otherwise return null.
     NDMaterial *getCopy(const std::string &) const;

     // Return the string "PressureDependMultiYield".
     const std::string &getType(void) const ;

     int sendSelf(Communicator &);  
     int recvSelf(const Communicator &);
     void Print(std::ostream &s, int flag =0) const;
    // RWB; PyLiq1 & TzLiq1 need to see the excess pore pressure and initial stresses.
    friend class PyLiq1;
    friend class TzLiq1;
  };
} // end of XC namespace

#endif