Method addResistingForceToNodalReaction

SRC/element/Element.cpp:637–680

int 
Element::addResistingForceToNodalReaction(int flag) 
{
  int result = 0;
  int numNodes = this->getNumExternalNodes();
  Node **theNodes = this->getNodePtrs();

  static Vector theVector(48);

  //
  // now determine the resisting force
  //

  const Vector *theResistingForce =0;
  if (flag == 0)
    theResistingForce = &(this->getResistingForce());
  else if (flag == 1) {
    theResistingForce = &(this->getResistingForceIncInertia());
  }
  else if (flag == 2)
    theResistingForce = &(this->getRayleighDampingForces());


  //
  // iterate over the elements nodes; determine nodes contribution & add it
  //

  int nodalDOFCount = 0;

  for (int i=0; i<numNodes; i++) {
    Node *theNode = theNodes[i];

    int numNodalDOF = theNode->getNumberDOF();
    theVector.resize(numNodalDOF);

    for (int j=0; j<numNodalDOF; j++) {
      (theVector)(j) = (*theResistingForce)(nodalDOFCount);
      nodalDOFCount++;
    }
    result +=theNode->addReactionForce(theVector, 1.0);
  }

  return result;
}

double Element::getCharacteristicLength(void)
{