Increase number of extra output quantities available for shell elemen…

…ts so we can write out all the BladeStiffenedShell DVs
A-CGray · May 19, 2023 · a49fbf4 · a49fbf4
1 parent 3e7a63f
commit a49fbf4
Show file tree

Hide file tree

Showing 3 changed files with 24 additions and 3 deletions.
diff --git a/src/constitutive/TACSBladeStiffenedShellConstitutive.cpp b/src/constitutive/TACSBladeStiffenedShellConstitutive.cpp
@@ -1199,12 +1199,21 @@ void TACSBladeStiffenedShellConstitutive::addFailureDVSens(
 // Retrieve the design variable for plotting purposes
 TacsScalar TACSBladeStiffenedShellConstitutive::evalDesignFieldValue(
     int elemIndex, const double pt[], const TacsScalar X[], int index) {
-  // TODO: Implement this
   switch (index) {
     case 0:
       return this->computeEffectiveThickness();
     case 1:
       return this->computeEffectiveBendingThickness();
+    case 2:
+      return this->panelLength;
+    case 3:
+      return this->stiffenerPitch;
+    case 4:
+      return this->panelThick;
+    case 5:
+      return this->stiffenerHeight;
+    case 6:
+      return this->stiffenerThick;
   }
   return 0.0;
 }

diff --git a/src/elements/TACSElementTypes.cpp b/src/elements/TACSElementTypes.cpp
@@ -85,7 +85,7 @@ int TacsGetOutputComponentCount(ElementType etype, int comp) {
     } else if (comp == TACS_OUTPUT_STRESSES) {
       return 9;
     } else if (comp == TACS_OUTPUT_EXTRAS) {
-      return 4;
+      return 8;
     } else if (comp == TACS_OUTPUT_LOADS) {
       return 6;
     }
@@ -305,6 +305,14 @@ const char *TacsGetOutputComponentName(ElementType etype, int comp, int index) {
           return "dv2";
         case 3:
           return "dv3";
+        case 4:
+          return "dv4";
+        case 5:
+          return "dv5";
+        case 6:
+          return "dv6";
+        case 7:
+          return "dv7";
         default:
           return NULL;
       }

diff --git a/src/elements/shell/TACSShellElement.h b/src/elements/shell/TACSShellElement.h
@@ -1229,7 +1229,11 @@ void TACSShellElement<quadrature, basis, director, model>::getOutputData(
         data[1] = con->evalDesignFieldValue(elemIndex, pt, X, 0);
         data[2] = con->evalDesignFieldValue(elemIndex, pt, X, 1);
         data[3] = con->evalDesignFieldValue(elemIndex, pt, X, 2);
-        data += 4;
+        data[4] = con->evalDesignFieldValue(elemIndex, pt, X, 3);
+        data[5] = con->evalDesignFieldValue(elemIndex, pt, X, 4);
+        data[6] = con->evalDesignFieldValue(elemIndex, pt, X, 5);
+        data[7] = con->evalDesignFieldValue(elemIndex, pt, X, 6);
+        data += 8;
       }
     }
   }