Merge pull request #680 from mhscott/mixedBeamColumn

Adding section displacement to element recorder

SRC/element/mixedBeamColumn/MixedBeamColumn3d.cpp

@@ -1417,6 +1417,13 @@ Response* MixedBeamColumn3d::setResponse(const char **argv, int argc,
 
   }
 
+  else if (strcmp(argv[0],"sectionDisplacements") == 0) {
+    if (argc > 1 && strcmp(argv[1],"local") == 0)
+      theResponse = new ElementResponse(this, 1111, Matrix(numSections,3));
+    else
+      theResponse = new ElementResponse(this, 111, Matrix(numSections,3));
+  }
+
     else if (strcmp(argv[0],"xaxis") == 0 || strcmp(argv[0],"xlocal") == 0)
       theResponse = new ElementResponse(this, 201, Vector(3));
 
@@ -1567,6 +1574,53 @@ int MixedBeamColumn3d::getResponse(int responseID, Information &eleInfo) {
 
   }
 
+  else if (responseID == 111 || responseID == 1111) {
+    double L = crdTransf->getInitialLength();
+    double pts[MAX_NUM_SECTIONS];
+    beamIntegr->getSectionLocations(numSections, L, pts);
+    // CBDI influence matrix
+    Matrix ls(numSections, numSections);
+    getCBDIinfluenceMatrix(numSections, pts, L, ls);
+    // Curvature vector
+    Vector kappaz(numSections); // about section z
+    Vector kappay(numSections); // about section y
+    for (int i = 0; i < numSections; i++) {
+      const ID &code = sections[i]->getType();
+      const Vector &e = sections[i]->getSectionDeformation();
+      int order = sections[i]->getOrder();
+      for (int j = 0; j < order; j++) {
+	if (code(j) == SECTION_RESPONSE_MZ)
+	  kappaz(i) += e(j);
+	if (code(j) == SECTION_RESPONSE_MY)
+	  kappay(i) += e(j);
+      }
+    }
+    // Displacement vector
+    Vector dispsy(numSections); // along local y
+    Vector dispsz(numSections); // along local z    
+    dispsy.addMatrixVector(0.0, ls, kappaz,  1.0);
+    dispsz.addMatrixVector(0.0, ls, kappay, -1.0);    
+    beamIntegr->getSectionLocations(numSections, L, pts);
+    static Vector uxb(3);
+    static Vector uxg(3);
+    Matrix disps(numSections,3);
+    static Vector vp(6);
+    vp = crdTransf->getBasicTrialDisp();
+    for (int i = 0; i < numSections; i++) {
+      uxb(0) = pts[i]*vp(0); // linear shape function
+      uxb(1) = dispsy(i);
+      uxb(2) = dispsz(i);
+      if (responseID == 111)
+	uxg = crdTransf->getPointGlobalDisplFromBasic(pts[i],uxb);
+      else
+	uxg = crdTransf->getPointLocalDisplFromBasic(pts[i],uxb);
+      disps(i,0) = uxg(0);
+      disps(i,1) = uxg(1);
+      disps(i,2) = uxg(2);            
+    }
+    return eleInfo.setMatrix(disps);
+  }
+
   else if (responseID >= 201 && responseID <= 203) {
     static Vector xlocal(3);
     static Vector ylocal(3);