00001 /* ****************************************************************** ** 00002 ** OpenSees - Open System for Earthquake Engineering Simulation ** 00003 ** Pacific Earthquake Engineering Research Center ** 00004 ** ** 00005 ** ** 00006 ** (C) Copyright 1999, The Regents of the University of California ** 00007 ** All Rights Reserved. ** 00008 ** ** 00009 ** Commercial use of this program without express permission of the ** 00010 ** University of California, Berkeley, is strictly prohibited. See ** 00011 ** file 'COPYRIGHT' in main directory for information on usage and ** 00012 ** redistribution, and for a DISCLAIMER OF ALL WARRANTIES. ** 00013 ** ** 00014 ** Developed by: ** 00015 ** Frank McKenna (fmckenna@ce.berkeley.edu) ** 00016 ** Gregory L. Fenves (fenves@ce.berkeley.edu) ** 00017 ** Filip C. Filippou (filippou@ce.berkeley.edu) ** 00018 ** ** 00019 ** ****************************************************************** */ 00020 00021 // $Revision: 1.1.1.1 $ 00022 // $Date: 2000/09/15 08:23:16 $ 00023 // $Source: /usr/local/cvs/OpenSees/SRC/analysis/fe_ele/penalty/PenaltyMP_FE.h,v $ 00024 00025 00026 #ifndef PenaltyMP_FE_h 00027 #define PenaltyMP_FE_h 00028 00029 // File: ~/analysis/fe_ele/penalty/PenaltyMP_FE.h 00030 // 00031 // Written: fmk 00032 // Created: 11/96 00033 // Revision: A 00034 // 00035 // Description: This file contains the class definition for PenaltyMP_FE. 00036 // PenaltyMP_FE is a subclass of FE_Element which handles MP_Constraints 00037 // using the penalty method. 00038 // 00039 // What: "@(#) PenaltyMP_FE.h, revA" 00040 00041 #include <FE_Element.h> 00042 #include <ID.h> 00043 #include <Matrix.h> 00044 #include <Vector.h> 00045 00046 #include <Element.h> 00047 #include <Integrator.h> 00048 #include <AnalysisModel.h> 00049 #include <Domain.h> 00050 #include <MP_Constraint.h> 00051 #include <Node.h> 00052 00067 class PenaltyMP_FE: public FE_Element 00068 { 00069 public: 00071 PenaltyMP_FE(Domain &theDomain, MP_Constraint &theMP, double alpha); 00072 00074 virtual ~PenaltyMP_FE(); 00075 00077 virtual int setID(void); 00079 virtual const Matrix &getTangent(Integrator *theIntegrator); 00081 virtual const Vector &getResidual(Integrator *theIntegrator); 00083 virtual const Vector &getTangForce(const Vector &x, double fact = 1.0); 00084 00085 protected: 00086 00087 private: 00089 void determineTangent(void); 00090 00092 MP_Constraint *theMP; 00094 Node *theConstrainedNode; 00096 Node *theRetainedNode; 00097 00099 Matrix *tang; 00101 Vector *resid; 00103 Matrix *C; 00104 double alpha; 00105 00106 00107 }; 00108 00109 #endif 00110 00111