BeamColumnJoint2d.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/>.
//----------------------------------------------------------------------------
/* ****************************************************************** **
**    OpenSees - Open System for Earthquake Engineering Simulation    **
**          Pacific Earthquake Engineering Research Center            **
**                                                                    **
**                                                                    **
** (C) Copyright 1999, The Regents of the University of California    **
** All Rights Reserved.                                               **
**                                                                    **
** Commercial use of this program without express permission of the   **
** University of California, Berkeley, is strictly prohibited.  See   **
** file 'COPYRIGHT'  in main directory for information on usage and   **
** redistribution,  and for a DISCLAIMER OF ALL WARRANTIES.           **
**                                                                    **
** Developed by:                                                      **
**   Frank McKenna (fmckenna@ce.berkeley.edu)                         **
**   Gregory L. Fenves (fenves@ce.berkeley.edu)                       **
**   Filip C. Filippou (filippou@ce.berkeley.edu)                     **
**                                                                    **
** ****************************************************************** */
                                                                        
// $Revision: 1.2 $
// $Date: 2004/10/06 19:21:45 $
// $Source: /usr/local/cvs/OpenSees/SRC/element/special/joint/BeamColumnJoint2d.h,v $
                                                                        
// Written: NM (nmitra@u.washington.edu)
// Created: April 2002
// Revised: August 2004
//
// Description: This file contains the class definition for beam-column joint.
// This element is a 4 noded 12 dof (3 dof at each node) finite area super-element introduced by
// Prof. Lowes and Arash. The element takes in 13 different material types in order to simulate
// the inelastic action observed in a reinforced beam column joint.
// Updates: Several concerning Joint formulation (presently a revised formulation for joints)
                                                                        
#ifndef BeamColumnJoint2d_h
#define BeamColumnJoint2d_h

#include <domain/mesh/element/ElemWithMaterial.h>
#include "domain/mesh/element/utils/physical_properties/BeamColumnJointPhysicalProperties.h"
#include <utility/matrix/Matrix.h>
#include <utility/matrix/Vector.h>

namespace XC {
class Node;
class Channel;
class FEM_ObjectBroker;
class Response;
class Renderer;
class UniaxialMaterial;

//
class BeamColumnJoint2d: public ElemWithMaterial<4,BeamColumnJointPhysicalProperties>
  {
  private:
  
    // private methods
    void getGlobalDispls(Vector&) ;
    void getBCJoint(void);
    void getdg_df(void);
    void getdDef_du(void);
    void matDiag(Vector, Matrix&);
    void getMatResponse(Vector, Vector&, Vector&);
    void formR(Vector);
    void formK(Vector);
    double getStepSize(double,double,Vector,Vector,Vector,Vector,double);
  
    int nodeDbTag, dofDbTag;
  
    // various other element parameters
    double elemActHeight;
    double elemActWidth;
    double elemWidth;
    double elemHeight;
    double HgtFac; 
    double WdtFac; 
  
    Vector Uecommit; 
    Vector UeIntcommit; 
    Vector UeprCommit; 
    Vector UeprIntCommit; 
    Matrix BCJoint; 
    Matrix dg_df; 
    Matrix dDef_du; 
  
    mutable Matrix K; 
    mutable Vector R; 
  public:
    // default constructor
    BeamColumnJoint2d(void); 
    // defined constructor
    BeamColumnJoint2d(int tag,int Nd1, int Nd2, int Nd3, int Nd4, const UniaxialMaterial &theMat1, const UniaxialMaterial &theMat2, const UniaxialMaterial &theMat3, const UniaxialMaterial &theMat4, const UniaxialMaterial &theMat5, const UniaxialMaterial &theMat6, const UniaxialMaterial &theMat7, const UniaxialMaterial &theMat8, const UniaxialMaterial &theMat9, const UniaxialMaterial &theMat10, const UniaxialMaterial &theMat11, const UniaxialMaterial &theMat12, const UniaxialMaterial &theMat13);
  
    BeamColumnJoint2d(int tag,int Nd1, int Nd2, int Nd3, int Nd4, const UniaxialMaterial &theMat1, const UniaxialMaterial &theMat2, const UniaxialMaterial &theMat3, const UniaxialMaterial &theMat4, const UniaxialMaterial &theMat5, const UniaxialMaterial &theMat6, const UniaxialMaterial &theMat7, const UniaxialMaterial &theMat8, const UniaxialMaterial &theMat9, const UniaxialMaterial &theMat10, const UniaxialMaterial &theMat11, const UniaxialMaterial &theMat12, const UniaxialMaterial &theMat13, double Hgtfac, double Wdtfac);
    Element *getCopy(void) const;

  bool  isSubdomain(void) { return false; } ;
  
  // return number of DOFs
  int getNumDOF(void) const;    
  
  // set domain performs check on dof and associativity with node
  void setDomain(Domain *theDomain);
  
  
  // commit state
  int commitState(void);
  
  // revert to last commit
  int revertToLastCommit(void);        
  
  // revert to start
  int revertToStart(void);        
  
  // determine current strain and set strain in material
  int update(void);
  
  
  // returns converged tangent stiffness matrix
  const Matrix &getTangentStiff(void) const;
  const Matrix &getInitialStiff(void) const;   
  
  // not required for this element formulation
  const Matrix &getDamp(void) const;
  const Matrix &getMass(void) const;
  
  // not required for this element formulation
  int addLoad(ElementalLoad *theLoad, double loadFactor);
  int addInertiaLoadToUnbalance(const Vector &accel);
  
  // get converged residual
  const Vector &getResistingForce(void) const;
  
  // get converged residual with inertia terms
  const Vector &getResistingForceIncInertia(void) const;            
  
  // public methods for element output for parallel and database processing
  int sendSelf(Communicator &);
  int recvSelf(const Communicator &);
  
  // print out element data
  void Print(std::ostream &s, int flag =0) const;    
  
  // implemented to print into file
  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