read in vtk data for stress

main.py

@@ -10,5 +10,5 @@ def main(fname):
     return
 
 if __name__ == '__main__':
-    fname = 'torsionbar_hor.stl'
+    fname = 'tensile.vtk'
     pl = Planner(fname, 1, 0.25)

utils/planner.py

@@ -24,10 +24,11 @@ class Planner:
 
         # read mesh data
         self.mesh = pv.read(mesh_file_name)
+        self.boundary = self.mesh.decimate_boundary()
 
         #voxelize
         hwratio = voxel_dim/self.lheight
-        m = self.mesh.copy()
+        m = self.boundary.copy()
         m.points[:,2] *= hwratio
         self.grid = pv.voxelize(m, density=self.voxel_dim)
         self.grid.points[:,2] *= 1/hwratio