Method getInitialDeformations

SRC/element/forceBeamColumn/ForceBeamColumn3d.cpp:2545–2611

  }

int
ForceBeamColumn3d::getInitialDeformations(Vector &v0)
{
  v0.Zero();
  if (numEleLoads < 1)
      return 0;

  double L = crdTransf->getInitialLength();
  double oneOverL = 1.0 / L;

  double xi[maxNumSections];
  beamIntegr->getSectionLocations(numSections, L, xi);

  double wt[maxNumSections];
  beamIntegr->getSectionWeights(numSections, L, wt);

  for (int i = 0; i < numSections; i++) {

      int order = sections[i]->getOrder();
      const ID &code = sections[i]->getType();

      double xL = xi[i];
      double xL1 = xL - 1.0;
      double wtL = wt[i] * L;

      static Vector sp;
      sp.setData(workArea, order);
      sp.Zero();

      this->computeSectionForces(sp, i);

      const Matrix &fse = sections[i]->getInitialFlexibility();

      static Vector e;
      e.setData(&workArea[order], order);

      e.addMatrixVector(0.0, fse, sp, 1.0);

      double dei, tmp;
      for (int ii = 0; ii < order; ii++) {
          dei = e(ii)*wtL;
          switch (code(ii)) {
          case SECTION_RESPONSE_P:
              v0(0) += dei;
              break;
          case SECTION_RESPONSE_MZ:
              v0(1) += xL1*dei; v0(2) += xL*dei;
              break;
          case SECTION_RESPONSE_VY:
              tmp = oneOverL*dei;
              v0(1) += tmp; v0(2) += tmp;
              break;
          case SECTION_RESPONSE_MY:
              v0(3) += xL1*dei; v0(4) += xL*dei;
              break;
          case SECTION_RESPONSE_VZ:
              tmp = oneOverL*dei;
              v0(3) += tmp; v0(4) += tmp;