Method addInertiaLoadSensitivityToUnbalance

SRC/domain/node/Node.cpp:1000–1040

int
Node::addInertiaLoadSensitivityToUnbalance(const Vector &accelG, double fact, bool somethingRandomInMotions)
{
  // simply return if node has no mass or R matrix
  if (mass == 0 || R == 0)
    return 0;

  // otherwise we must determine MR accelG
  if (accelG.Size() != R->noCols()) {
    opserr << "Node::addInertiaLoadToUnbalance - accelG not of correct dimension";
    return -1;
  }

  // if no load yet create it and assign
  if (unbalLoad == 0) {
      unbalLoad = new Vector(numberDOF); 
      if (unbalLoad == 0 || unbalLoad->Size() != numberDOF) {
	  opserr << "FATAL Node::addunbalLoad - ran out of memory\n";
	  exit(-1);
      }  
  }

  // form - fact * M*R*accelG and add it to the unbalanced load
  //(*unbalLoad) -= ((*mass) * (*R) * accelG)*fact;


	Matrix massSens(mass->noRows(),mass->noCols());
	massSens = this->getMassSensitivity();

	Matrix MR(mass->noRows(), R->noCols());

	if (somethingRandomInMotions) {
	  MR.addMatrixProduct(0.0, *mass, *R, 1.0);
	}
	else {
	  MR.addMatrixProduct(0.0, massSens, *R, 1.0);
	}
	unbalLoad->addMatrixVector(1.0, MR, accelG, -fact);

  return 0;
}