EightNodeBrick_u_p_U.h

//----------------------------------------------------------------------------
//  XC program; finite element analysis code
//  for structural analysis and design.
//
//  Copyright (C)  Luis Claudio 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/>.
//----------------------------------------------------------------------------
//
// COPYRIGHT (C):     :-))
// PROJECT:           Object Oriented Finite Element Program
// FILE:             EightNodeBrick_u_p_U.h
// CLASS:            EightNodeBrick_u_p_U
// MEMBER FUNCTIONS:
//
// MEMBER VARIABLES
//
// PURPOSE:           Finite Element Class for coupled system
//  "Coupled system": Solid and fluid coexist.
//                    u-- Solid displacement
//                    p-- Pore pressure
//                    U-- Absolute fluid displacement
// RETURN:
// VERSION:
// LANGUAGE:          C++.ver >= 3.0
// TARGET OS:         DOS || UNIX || . . .
// DESIGNER:          Boris Jeremic, Xiaoyan Wu, Chao Cheng (main for last revision)
// PROGRAMMER:        Boris Jeremic, Xiaoyan Wu, Zhaohui Yang, Zhao Cheng (main for last revision)
// DATE:              Sept. 2001
// UPDATE HISTORY:    Modified from EightNodeBrick.h.  Reorganized a lot by Xiaoyan
//                    01/24/2002    Xiaoyan
//                    Add the permeability tensor and ks, kf  to the constructor  Xiaoyan
//
//
//
//                    Clean-up and re-write by Zhao Cheng, 10/20/2004
//
//


#ifndef EIGHTNODEBRICK_U_P_U_H
#define EIGHTNODEBRICK_U_P_U_H

#include <domain/mesh/element/volumetric/BrickBase.h>
#include "domain/mesh/element/utils/body_forces/BodyForces3D.h"

namespace XC {
  class BJtensor;
  class NDMaterial;

//
class EightNodeBrick_u_p_U: public BrickBase
  {
  private:
    static Matrix K; 
    static Matrix C; 
    static Matrix M; 
    static Vector P; 

    static const int  Num_IntegrationPts;
    static const int  Num_TotalGaussPts;
    static const int  Num_Nodes;
    static const int  Num_Dim;
    static const int  Num_Dof;
    static const int  Num_ElemDof;
    static const double pts[2]; 
    static const double wts[2]; 
    static BJtensor perm; 

    BodyForces3D bf; 
    double poro; 
    double alpha; 
    double rho_s; 
    double rho_f; 
    double ks; 
    double kf; 
    double pressure; 

    Vector *eleQ;
    mutable Matrix *Ki;
  public:
    EightNodeBrick_u_p_U(int element_number,
                   int node_numb_1, int node_numb_2, int node_numb_3, int node_numb_4,
                   int node_numb_5, int node_numb_6, int node_numb_7, int node_numb_8,
                   NDMaterial *Globalmmodel,  const BodyForces3D &,
                   double nn, double alf, double rs,double rf,
                   double permb_x, double permb_y, double permb_z,
                   double kks, double kkf, double pp);
    EightNodeBrick_u_p_U(void);
    Element *getCopy(void) const;
    ~EightNodeBrick_u_p_U(void);
    // public methods to obtain information about dof & connectivity
    int getNumDOF(void) const;
    void setDomain(Domain *theDomain);

    // public methods to set the state of the element    
    int update(void);

    // public methods to obtain stiffness, mass, damping and residual information
    const Matrix &getTangentStiff(void) const;
    const Matrix &getInitialStiff(void) const;
    const Matrix &getDamp(void) const;
    const Matrix &getMass(void) const;

    void zeroLoad(void);
    int addLoad(ElementalLoad *theLoad, double loadFactor);
    int addInertiaLoadToUnbalance(const Vector &accel);
    const Vector &getResistingForce(void) const;
    const Vector &getResistingForceIncInertia(void) const;

    int sendSelf(CommParameters &);
    int recvSelf(const CommParameters &);

    void Print(std::ostream &s, int flag =0);

    Response *setResponse(const std::vector<std::string> &argv, Information &eleInfo);
    int getResponse(int responseID, Information &eleInformation);

    //int setParameter(const std::vector<std::string> &argv, Parameter &param);
    //int updateParameter (int parameterID, Information &info);

  private:
    BJtensor shapeFunction(double, double, double) const;
    BJtensor shapeFunctionDerivative(double, double, double) const;
    BJtensor getGaussPts(void);
    BJtensor getNodesCrds(void) const;
    BJtensor getNodesDisp(void) const;
    BJtensor Jacobian_3D(BJtensor dh) const;
    BJtensor Jacobian_3Dinv(BJtensor dh) const;
    BJtensor dh_Global(BJtensor dh) const;

    BJtensor getStiffnessTensorKep() const;
    BJtensor getStiffnessTensorG12() const;
    BJtensor getStiffnessTensorP() const;
    BJtensor getMassTensorMsf() const;
    BJtensor getDampTensorC123(void) const;
    const Matrix &getStiff(int Ki_flag) const;
    double getPorePressure(double, double, double);
    Vector getExForceS();
    Vector getExForceF();
  };
} // end of XC namespace


#endif