FluidSolidPorousMaterial.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/>.
//----------------------------------------------------------------------------
// $Revision: 1.11 $
// $Date: 2003/10/07 21:20:48 $
// $Source: /usr/local/cvs/OpenSees/SRC/material/nD/soil/FluidSolidPorousMaterial.h,v $
                                                                        
// Written: ZHY
// Created: August 2000
// Revision: A
//
// Description: This file contains the class prototype for FluidSolidPorousMaterial.
//
// What: "@(#) FluidSolidPorousMaterial.h, revA"

#ifndef FluidSolidPorousMaterial_h
#define FluidSolidPorousMaterial_h

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

namespace XC {
class Response;

//
//
class FluidSolidPorousMaterial: public SoilMaterialBase
  {
  private:
    static std::vector<double> combinedBulkModulusx;
    static double pAtm;
    NDMaterial *theSoilMaterial;
    mutable double trialExcessPressure;
    double currentExcessPressure;
    double trialVolumeStrain;
    double currentVolumeStrain;
    mutable double initMaxPress;

    static Vector workV3;
    static Vector workV6;
    static Matrix workM3;
    static Matrix workM6;

    void free_mem(void);
    void alloc(const NDMaterial *);
  protected:
    int sendData(Communicator &);
    int recvData(const Communicator &);
  public:
     // Initialization constructor
     FluidSolidPorousMaterial(int tag= 0);
     FluidSolidPorousMaterial(int tag, int nd,const NDMaterial &soilMat, double combinedBulkModul, double atm=101.);
     FluidSolidPorousMaterial(const FluidSolidPorousMaterial &);
     FluidSolidPorousMaterial &operator=(const FluidSolidPorousMaterial &tro);

     // Destructor: clean up memory storage space.
     virtual ~FluidSolidPorousMaterial();

     // Sets the values of the trial strain tensor.
     int setTrialStrain(const Vector &strain);

     // Sets the values of the trial strain and strain rate tensors.
     int setTrialStrain(const Vector &v, const Vector &r);

     int setTrialStrainIncr(const Vector &v);
     int setTrialStrainIncr(const Vector &v, const Vector &r);

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

     double getRho(void) const;

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

     // Accepts the current trial strain values as being on the solution path, and updates 
     // all model parameters related to stress/strain states. Return 0 on success.
     int commitState(void);

     // Revert the stress/strain states to the last committed states. Return 0 on success.
     int revertToLastCommit(void);

     int revertToStart(void);

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

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

     int sendSelf(Communicator &);  
     int recvSelf(const Communicator &);     
     Response *setResponse(const std::vector<std::string> &argv, Information &matInfo);
     int getResponse(int responseID, Information &matInformation);
     void Print(std::ostream &s, int flag =0) const;

     int updateParameter(int responseID, Information &eleInformation);

     // RWB; PyLiq1 & TzLiq1 need to see the excess pore pressure and initial stresses.
     friend class PyLiq1;
     friend class TzLiq1;
  };
} // end of XC namespace

#endif