PEER NEES NEESit

• HOME
• USER
• DEVELOPER
• PROJECTS
• SUPPORT
• SITE MAP

• Dev Doc
• API
• Source
• Download
• Builds
• Bug Reports
• Message Board

LagrangeMP_FE.cpp Go to the documentation of this file. /* ****************************************************************** ** ** 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.5 $ // $Date: 2006/02/08 20:20:00 $ // $Source: /usr/local/cvs/OpenSees/SRC/analysis/fe_ele/lagrange/LagrangeMP_FE.cpp,v $ // Written: fmk // Created: 02/99 // Revision: A // // Purpose: This file contains the code for implementing the methods // of the LagrangeMP_FE class interface. #include <LagrangeMP_FE.h> #include <stdlib.h> #include <Element.h> #include <Domain.h> #include <Node.h> #include <DOF_Group.h> #include <Integrator.h> #include <Subdomain.h> #include <AnalysisModel.h> #include <Matrix.h> #include <Vector.h> #include <Node.h> #include <MP_Constraint.h> #include <DOF_Group.h> LagrangeMP_FE::LagrangeMP_FE(int tag, Domain &theDomain, MP_Constraint &TheMP, DOF_Group &theGroup, double Alpha) :FE_Element(tag, 3,(tag, TheMP.getConstrainedDOFs()).Size()+ (TheMP.getRetainedDOFs()).Size() + (TheMP.getRetainedDOFs()).Size()), alpha(Alpha), theMP(&TheMP), theConstrainedNode(0), theRetainedNode(0), theDofGroup(&theGroup), tang(0), resid(0) { const Matrix &constraint = theMP->getConstraint(); int noRows = constraint.noRows(); int noCols = constraint.noCols(); int size = 2*noRows+noCols; tang = new Matrix(size,size); resid = new Vector(size); if (tang == 0 || resid == 0 || tang->noCols() == 0 || resid->Size() == 0) { opserr << "FATAL LagrangeMP_FE::LagrangeMP_FE() - out of memory\n"; exit(-1); } tang->Zero(); resid->Zero(); theRetainedNode = theDomain.getNode(theMP->getNodeRetained()); theConstrainedNode = theDomain.getNode(theMP->getNodeConstrained()); if (theRetainedNode == 0) { opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()"; opserr << "- no asscoiated Retained Node\n"; exit(-1); } if (theConstrainedNode == 0) { opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()"; opserr << "- no asscoiated Constrained Node\n"; exit(-1); } if (theMP->isTimeVarying() == false) { this->determineTangent(); } // set the myDOF_Groups tags indicating the attached id's of the // DOF_Group objects DOF_Group *theConstrainedNodesDOFs = theConstrainedNode->getDOF_GroupPtr(); if (theConstrainedNodesDOFs == 0) { opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()"; opserr << " - no DOF_Group with Constrained Node\n"; exit(-1); } DOF_Group *theRetainedNodesDOFs = theRetainedNode->getDOF_GroupPtr(); if (theRetainedNodesDOFs == 0) { opserr << "WARNING LagrangeMP_FE::LagrangeMP_FE()"; opserr << " - no DOF_Group with Retained Node\n"; exit(-1); } myDOF_Groups(0) = theConstrainedNodesDOFs->getTag(); myDOF_Groups(1) = theRetainedNodesDOFs->getTag(); myDOF_Groups(2) = theDofGroup->getTag(); } LagrangeMP_FE::~LagrangeMP_FE() { if (tang != 0) delete tang; if (resid != 0) delete resid; } // void setID(int index, int value); // Method to set the correMPonding index of the ID to value. int LagrangeMP_FE::setID(void) { int result = 0; // first determine the IDs in myID for those DOFs marked // as constrained DOFs, this is obtained from the DOF_Group // associated with the constrained node if (theConstrainedNode == 0) { opserr << "WARNING LagrangeMP_FE::setID(void)"; opserr << "- no asscoiated Constrained Node\n"; return -1; } DOF_Group *theConstrainedNodesDOFs = theConstrainedNode->getDOF_GroupPtr(); if (theConstrainedNodesDOFs == 0) { opserr << "WARNING LagrangeMP_FE::setID(void)"; opserr << " - no DOF_Group with Constrained Node\n"; return -2; } const ID &constrainedDOFs = theMP->getConstrainedDOFs(); const ID &theConstrainedNodesID = theConstrainedNodesDOFs->getID(); int size1 = constrainedDOFs.Size(); for (int i=0; i<size1; i++) { int constrained = constrainedDOFs(i); if (constrained < 0 || constrained >= theConstrainedNode->getNumberDOF()) { opserr << "WARNING LagrangeMP_FE::setID(void) - unknown DOF "; opserr << constrained << " at Node\n"; myID(i) = -1; // modify so nothing will be added to equations result = -3; } else { if (constrained >= theConstrainedNodesID.Size()) { opserr << "WARNING LagrangeMP_FE::setID(void) - "; opserr << " Nodes DOF_Group too small\n"; myID(i) = -1; // modify so nothing will be added to equations result = -4; } else myID(i) = theConstrainedNodesID(constrained); } } // now determine the IDs for the retained dof's if (theRetainedNode == 0) { opserr << "WARNING LagrangeMP_FE::setID(void)"; opserr << "- no asscoiated Retained Node\n"; return -1; } DOF_Group *theRetainedNodesDOFs = theRetainedNode->getDOF_GroupPtr(); if (theRetainedNodesDOFs == 0) { opserr << "WARNING LagrangeMP_FE::setID(void)"; opserr << " - no DOF_Group with Retained Node\n"; return -2; } const ID &RetainedDOFs = theMP->getRetainedDOFs(); const ID &theRetainedNodesID = theRetainedNodesDOFs->getID(); int size2 = RetainedDOFs.Size(); for (int j=0; j<size2; j++) { int retained = RetainedDOFs(j); if (retained < 0 || retained >= theRetainedNode->getNumberDOF()) { opserr << "WARNING LagrangeMP_FE::setID(void) - unknown DOF "; opserr << retained << " at Node\n"; myID(j+size1) = -1; // modify so nothing will be added result = -3; } else { if (retained >= theRetainedNodesID.Size()) { opserr << "WARNING LagrangeMP_FE::setID(void) - "; opserr << " Nodes DOF_Group too small\n"; myID(j+size1) = -1; // modify so nothing will be added result = -4; } else myID(j+size1) = theRetainedNodesID(retained); } } // finally set the ID corresponding to the ID's at the LagrangeDOF_Group const ID &theGroupsID = theDofGroup->getID(); int size3 = theGroupsID.Size(); for (int k=0; k<size3; k++) myID(k+size1+size2) = theGroupsID(k); return result; } const Matrix & LagrangeMP_FE::getTangent(Integrator *theNewIntegrator) { if (theMP->isTimeVarying() == true) this->determineTangent(); return *tang; } const Vector & LagrangeMP_FE::getResidual(Integrator *theNewIntegrator) { return *resid; } const Vector & LagrangeMP_FE::getTangForce(const Vector &disp, double fact) { opserr << "WARNING lagrangeMP_FE::getTangForce() - not yet implemented\n"; return *resid; } const Vector & LagrangeMP_FE::getK_Force(const Vector &disp, double fact) { opserr << "WARNING lagrangeMP_FE::getK_Force() - not yet implemented\n"; return *resid; } const Vector & LagrangeMP_FE::getKi_Force(const Vector &disp, double fact) { opserr << "WARNING LagrangeMP_FE::getKi_Force() - not yet implemented\n"; return *resid; } const Vector & LagrangeMP_FE::getC_Force(const Vector &disp, double fact) { opserr << "WARNING lagrangeMP_FE::getC_Force() - not yet implemented\n"; return *resid; } const Vector & LagrangeMP_FE::getM_Force(const Vector &disp, double fact) { opserr << "WARNING lagrangeMP_FE::getM_Force() - not yet implemented\n"; return *resid; } void LagrangeMP_FE::determineTangent(void) { const Matrix &constraint = theMP->getConstraint(); int noRows = constraint.noRows(); int noCols = constraint.noCols(); int n = noRows+noCols; tang->Zero(); for (int j=0; j<noRows; j++) { (*tang)(n+j, j) = -alpha; (*tang)(j, n+j) = -alpha; } for (int i=0; i<noRows; i++) for (int j=0; j<noCols; j++) { double val = constraint(i,j) * alpha; (*tang)(n+i, j+noRows) = val; (*tang)(noRows+j, n+i) = val; } } Generated on Mon Oct 23 15:04:57 2006 for OpenSees by  1.5.0

---

opensees-support @ berkeley.edu

Supported by the National Science Foundation

©2006, UC Regents