BeamColumnJoint3d.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/>.
//----------------------------------------------------------------------------
/* ****************************************************************** **
**    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/BeamColumnJoint3d.h,v $
                                                                        
// Written: NM (nmitra@u.washington.edu)
// Created: Feb 2003
// Updated: September 2004
//
// Description: This file contains the class defination for beam-column joint.
// This element is a 4 noded 24 dof (6 dof at each node) finite area super-element, being a slight
// variation of the 2d one. The element takes in 13 different material types in order to simulate
// the inelastic action observed in a reinforced beam column joint. Though it has 6 dof per node 
// the out of the plane nodal dof are constrained or fixed and the inplane nodal dof are activated.
                                                                        
#ifndef BeamColumnJoint3d_h
#define BeamColumnJoint3d_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 BeamColumnJoint3d : 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);
  void formTransfMat();
  double getStepSize(double,double,Vector,Vector,Vector,Vector,double);
  
  int nodeDbTag, dofDbTag;
  
  // various other element parameters
  Vector Node1; Vector Node2; Vector Node3; Vector Node4;
  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; 
  
  // static transformation matrices
  static Matrix Transf;
  static Matrix Tran;
  
  public:
    // default constructor
    BeamColumnJoint3d(); 
  
  // defined constructor
  BeamColumnJoint3d(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);
  
  BeamColumnJoint3d(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(CommParameters &);
  int recvSelf(const CommParameters &);
  
  // print out element data
  void Print(std::ostream &s, int flag =0);    
  
  // 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