Method getd_hat

SRC/element/mixedBeamColumn/MixedBeamColumn2d.cpp:1518–1552

}

Vector
MixedBeamColumn2d::getd_hat(int sec, const Vector &v, double L, bool geomLinear){
  double xi[MAX_NUM_SECTIONS];
  beamIntegr->getSectionLocations(numSections, L, xi);

  double x, C, E, F;
  Vector D_hat(NDM_SECTION);
  D_hat.Zero();

  x = L*xi[sec];
  C = 1/L;
  E = -4/L + 6*x/(L*L);
  F = -2/L + 6*x/(L*L);

  if (geomLinear) {

    D_hat(0) =  C*v(0);
    D_hat(1) =  E*v(1) + F*v(2);

  } else {

    double A,B;
    A = 1 - 4*(x/L) + 3*pow(x/L,2);
    B =   - 2*(x/L) + 3*pow(x/L,2);

    D_hat(0) =  C*v(0) +
                0.5*A*A*v(1)*v(1) +
                    A*B*v(1)*v(2) +
                0.5*B*B*v(2)*v(2);
    D_hat(1) =  E*v(1) + F*v(2);

  }

  return D_hat;
}

Matrix
MixedBeamColumn2d::getKg(int sec, double P, double L) {