TwentyNodeBrick_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:             TwentyNodeBrick_u_p_U.h
// CLASS:            TwentyNodeBrick_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
// PROGRAMMER:        Boris Jeremic, Xiaoyan Wu
// DATE:              Aug. 2001
// UPDATE HISTORY:    Modified from EightNodeBrick_u_p_U.h.   Aug. 2001           
//                    01/16/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 TWENTYNODEBRICK_U_P_U_H
#define TWENTYNODEBRICK_U_P_U_H



#include <domain/mesh/element/ElemWithMaterial.h>
#include "domain/mesh/element/utils/physical_properties/NDMaterialPhysicalProperties.h"
#include "domain/mesh/element/utils/body_forces/BodyForces3D.h"

namespace XC {
  class BJtensor;
  class NDMaterial;

//
class TwentyNodeBrick_u_p_U: public ElemWithMaterial<20,NDMaterialPhysicalProperties>
  {
  private:
    static Matrix K;               // Stiffness
    static Matrix C;           // Damping
    static Matrix M;           // Mass
    static Vector P;           // Residual

    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[3];    // Stores quadrature points
    static const double wts[3];    // Stores quadrature weights
    static BJtensor perm;      // Permeability = k/(rho_f*g)

    BodyForces3D bf; 
    double poro;                   // Porosity
    double alpha;          // Coefficient for soil (approximate equal 1)
    double rho_s;              // Solid density
    double rho_f;              // Fluid density
    double ks;                     // Bulk modulus of solid
    double kf;                     // Bulk modulus of fluid
    double pressure;           // Normal surface traction (pressure) over entire element  //?
    
    Vector *eleQ;                  
    mutable Matrix *Ki;
  private:
    static BJtensor shapeFunction(double, double, double);
    static BJtensor shapeFunctionDerivative(double, double, double);
    BJtensor getGaussPts(void);
    BJtensor getNodesCrds(void) const;
    BJtensor getNodesDisp(void);
    BJtensor Jacobian_3D(BJtensor dh) const;
    BJtensor Jacobian_3Dinv(BJtensor dh) const;
    BJtensor dh_Global(BJtensor dh) const;

    BJtensor getStiffnessTensorKep(void) const;
    BJtensor getStiffnessTensorG12(void) const;
    BJtensor getStiffnessTensorP(void) const;
    BJtensor getMassTensorMsf(void) const;
    BJtensor getDampTensorC123(void) const;
    const Matrix &getStiff(int Ki_flag) const;
    double getPorePressure(double, double, double) const;
    Vector getExForceS();
    Vector getExForceF();
        
  public:
    TwentyNodeBrick_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,
                   int node_numb_9, int node_numb_10, int node_numb_11, int node_numb_12,
                   int node_numb_13, int node_numb_14, int node_numb_15, int node_numb_16,
                   int node_numb_17, int node_numb_18, int node_numb_19, int node_numb_20,
           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); 
   TwentyNodeBrick_u_p_U(void);
   Element *getCopy(void) const;
   ~TwentyNodeBrick_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);                    
};
} // end of XC namespace


#endif