Method computeReactions

SRC/element/forceBeamColumn/ForceBeamColumn2d.cpp:1067–1127

}
    
void
ForceBeamColumn2d::computeReactions(double *p0)
{
  int type;
  double L = crdTransf->getInitialLength();
  
  for (int i = 0; i < numEleLoads; i++) {
    
    double loadFactor = eleLoadFactors[i];
    const Vector &data = eleLoads[i]->getData(type, loadFactor);

    if (type == LOAD_TAG_Beam2dUniformLoad) {
      double wa = data(1)*loadFactor;  // Axial
      double wy = data(0)*loadFactor;  // Transverse
      
      p0[0] -= wa*L;
      double V = 0.5*wy*L;
      p0[1] -= V;
      p0[2] -= V;
    }
    if (type == LOAD_TAG_BeamUniformMoment) {
      double mz = data(2)*loadFactor;  // About z
      
      p0[1] += mz;
      p0[2] -= mz;
    }    
    else if (type == LOAD_TAG_Beam2dPartialUniformLoad) {
      double waa = data(2)*loadFactor;  // Axial
      double wab = data(3)*loadFactor;  // Axial
      double wya = data(0)*loadFactor;  // Transverse
      double wyb = data(1)*loadFactor;  // Transverse      
      double a = data(4)*L;
      double b = data(5)*L;

      p0[0] -= waa*(b-a) + 0.5*(wab-waa)*(b-a);
      double Fy = wya*(b-a);
      double c = a + 0.5*(b-a);
      p0[1] -= Fy*(1-c/L);
      p0[2] -= Fy*c/L;
      Fy = 0.5*(wyb-wya)*(b-a);
      c = a + 2.0/3.0*(b-a);
      p0[1] -= Fy*(1-c/L);
      p0[2] -= Fy*c/L;      
    }
    else if (type == LOAD_TAG_Beam2dPointLoad) {
      double P = data(0)*loadFactor;
      double N = data(1)*loadFactor;
      double aOverL = data(2);
      
      if (aOverL < 0.0 || aOverL > 1.0)
	continue;
      
      double V1 = P*(1.0-aOverL);
      double V2 = P*aOverL;
      
      p0[0] -= N;
      p0[1] -= V1;
      p0[2] -= V2;