Method getKg

SRC/element/mixedBeamColumn/MixedBeamColumn3d.cpp:1747–1772

}

Matrix MixedBeamColumn3d::getKg(int sec, double P, double L) {
  double xi[MAX_NUM_SECTIONS];
  beamIntegr->getSectionLocations(numSections, L, xi);

  double temp_x, temp_A, temp_B;

  temp_x = L * xi[sec];

  Matrix kg(NDM_NATURAL,NDM_NATURAL);
  kg.Zero();

  temp_A = 1 - 4 * temp_x / L + 3 * ( temp_x * temp_x ) / ( L * L );
  temp_B = - 2 * temp_x / L + 3 * ( temp_x * temp_x ) / ( L * L );

  kg(0,0) = P / ( L * L );
  kg(1,1) = P * temp_A * temp_A;
  kg(1,3) = P * temp_A * temp_B;
  kg(2,2) = P * temp_A * temp_A;
  kg(2,4) = P * temp_A * temp_B;
  kg(3,1) = P * temp_A * temp_B;
  kg(3,3) = P * temp_B * temp_B;
  kg(4,2) = P * temp_A * temp_B;
  kg(4,4) = P * temp_B * temp_B;

  return kg;
}

Matrix MixedBeamColumn3d::getMd(int sec, Vector dShapeFcn, Vector dFibers, double L) {
  double xi[MAX_NUM_SECTIONS];