Swept Surface

examples/show_swept.py

@@ -0,0 +1,157 @@
+import gustaf as gus
+import numpy as np
+
+import splinepy
+
+if __name__ == "__main__":
+
+    ### TRAJECTORY ###
+
+    # define a hook-trajectory
+    dict_trajectory = {
+        "degrees": [3],
+        "knot_vectors": [
+            [0.0, 0.0, 0.0, 0.0, 0.2, 0.4, 0.6, 0.8, 0.9, 1.0, 1.0, 1.0, 1.0]
+        ],
+        "control_points": np.array(
+            [
+                [0.5, 0],
+                [0.5, 2],
+                [1.0, 3],
+                [2.0, 4],
+                [2.15, 5],
+                [1.8, 5.9],
+                [1.0, 6.2],
+                [-0.25, 6],
+                [-0.5, 5],
+            ]
+        ),
+    }
+
+    # create spline of trajectory dict
+    trajectory = splinepy.BSpline(**dict_trajectory)
+
+    # refine trajectory by inserting knots and control points
+    trajectory.uniform_refine(0, 1)
+
+    ### CROSS SECTIONS ###
+
+    # 1. create a circular 1D-line-cross-section
+    cross_section_circle = splinepy.helpme.create.circle(0.5).nurbs
+
+    # 2. create a circular 2D-surface-cross-section
+    cross_section_disk = splinepy.helpme.create.surface_circle(0.5).nurbs
+
+    # 3. create a rectangular 2D-surface-cross-section
+    cross_section_plate = splinepy.helpme.create.box(1, 1).nurbs
+
+    # Define a custom normal vector for the cross-section when:
+    # a) The cross-section is not planar in the x-y plane, or
+    # b) Crooked sweeping is desired
+    cs_nv = np.array([1, 0, 1])
+
+    ### SWEEP ###
+
+    # create swept surface
+    swept_surface_circle = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_circle,
+        cross_section_normal=None,
+        set_on_trajectory=False,
+        rotation_adaption=None,
+    )
+
+    # create crooked swept solid (circular nr. 1)
+    # the cross-section's normal vector is not default; crooked sweeping
+    swept_solid_disk_1 = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_disk,
+        cross_section_normal=cs_nv,
+        set_on_trajectory=False,
+        rotation_adaption=None,
+    )
+
+    # create swept solid (circular nr. 2)
+    # the cross-sections are set on the trajectory's control points (default)
+    swept_solid_disk_2 = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_disk,
+        cross_section_normal=None,
+        set_on_trajectory=False,
+        rotation_adaption=None,
+    )
+
+    # create swept solid (circular nr. 3)
+    # the cross-sections are set on the trajectory's evaluation points
+    swept_solid_disk_3 = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_disk,
+        cross_section_normal=None,
+        set_on_trajectory=True,
+        rotation_adaption=None,
+    )
+
+    # create swept solid (rectangular nr. 1)
+    swept_solid_plate_1 = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_plate,
+        cross_section_normal=None,
+        set_on_trajectory=False,
+        rotation_adaption=None,
+    )
+
+    # create swept solid (rectangular nr. 2)
+    # rotation adaption with 45 degrees
+    swept_solid_plate_2 = splinepy.helpme.create.swept(
+        trajectory=trajectory,
+        cross_section=cross_section_plate,
+        cross_section_normal=None,
+        set_on_trajectory=False,
+        rotation_adaption=45 * np.pi / 180,
+    )
+
+    ### VISUALIZATION ###
+
+    # first window: swept surface
+    gus.show(
+        ["Trajectory", trajectory],
+        ["1D Cross Section", cross_section_circle],
+        ["Swept Surface", swept_surface_circle],
+        resolution=51,
+        control_mesh=False,
+        control_point_ids=False,
+    )
+
+    # adjust show options
+    swept_solid_disk_2.show_options["alpha"] = 0.3
+    swept_solid_disk_3.show_options["alpha"] = 0.3
+    trajectory.show_options["control_points"] = False
+
+    # second window: swept solids (circular)
+    gus.show(
+        ["2D Cross Section", cross_section_disk],
+        ["Swept Solid - Crooked Sweeping", swept_solid_disk_1],
+        [
+            "Swept Solid - Set on Control Points",
+            swept_solid_disk_2,
+            trajectory,
+        ],
+        [
+            "Swept Solid - Set on Evaluation Points",
+            swept_solid_disk_3,
+            trajectory,
+        ],
+        resolution=51,
+        control_mesh=False,
+        control_point_ids=False,
+    )
+
+    # third window: swept solids (rectangular)
+    gus.show(
+        ["New Cross Section", cross_section_plate],
+        ["Swept Solid without Rotation", swept_solid_plate_1],
+        ["Swept Solid with 45° Rotation", swept_solid_plate_2],
+        resolution=51,
+        control_mesh=False,
+        control_point_ids=False,
+    )

splinepy/bspline.py

@@ -319,7 +319,6 @@ def remove_knots(self, parametric_dimension, knots, tolerance=None):
 
         self._logd(f"Tried to remove {len(knots)} knot(s).")
         self._logd(f"Actually removed {sum(removed)} knot(s).")
-
         return removed
 
     def normalize_knot_vectors(self):