Method addLoad

SRC/element/componentElement/ComponentElement3d.cpp:953–1009

}

int 
ComponentElement3d::addLoad(ElementalLoad *theLoad, double loadFactor)
{
  int type;
  const Vector &data = theLoad->getData(type, loadFactor);
  double L = theCoordTransf->getInitialLength();

  if (type == LOAD_TAG_Beam3dUniformLoad) {
    double wy = data(0)*loadFactor;  // Transverse
    double wz = data(1)*loadFactor;  // Transverse
    double wx = data(2)*loadFactor;  // Axial (+ve from node I to J)
    
    double Vy = 0.5*wy*L;
    double Mz = Vy*L/6.0; // wy*L*L/12
    double Vz = 0.5*wz*L;
    double My = Vz*L/6.0; // wz*L*L/12
    double P = wx*L;

    // Reactions in basic system
    p0[0] -= P;
    p0[1] -= Vy;
    p0[2] -= Vy;
    p0[3] -= Vz;
    p0[4] -= Vz;

    // Fixed end forces in basic system
    q0[0] -= 0.5*P;
    if (end1zHinge != 0 && end2zHinge != 0) {
      q0[1] -= Mz;
      q0[2] += Mz;
    }
    if (end1zHinge == 0 && end2zHinge != 0) {
      q0[2] += wy*L*L/8;
    }
    if (end1zHinge != 0 && end2zHinge == 0) {
      q0[1] -= wy*L*L/8;
    }
    
    if (end1yHinge != 0 && end2yHinge != 0) {
      q0[3] += My;
      q0[4] -= My;
    }
    if (end1yHinge == 0 && end2yHinge != 0) {
      q0[4] -= wz*L*L/8;
    }
    if (end1yHinge != 0 && end2yHinge == 0) {
      q0[3] += wz*L*L/8;
    }    
  }

  else {
    opserr << "ComponentElement3d::addLoad()  -- load type unknown for element with tag: " << this->getTag() << endln;
    return -1;
  }

  return 0;
}