feat: Integrate geometry imports

doc/changelog.d/952.added.md

@@ -0,0 +1 @@
+feat: Integrate geometry imports

examples/gallery/10_wheel_ground_contact_patch.py

@@ -100,13 +100,13 @@
 
 ###################
 # Visualize results
-# =================
+# ~~~~~~~~~~~~~~~~~
 # .. code-block:: python
 #
 #   display = PrimePlotter()
 #   display.plot(model, scope=prime.ScopeDefinition(model, label_expression="ground, wheel"))
 #   display.show()
-
+#
 print(model)
 
 ###############################################################################
@@ -167,8 +167,8 @@
 #          model, scope=prime.ScopeDefinition(model, label_expression="ground, patch*, wheel")
 #   )
 #   display.show()
-
-###############################################################################
+#
+#######################
 # Wrap the fluid region
 # ~~~~~~~~~~~~~~~~~~~~~
 # The largest internal region in this instance is the fluid region around the wheel.
@@ -192,7 +192,7 @@
 
 #########################
 # Open a pyvistaqt window
-# =======================
+# ~~~~~~~~~~~~~~~~~~~~~~~
 # .. code-block:: python
 #
 #     display = PrimePlotter()
@@ -201,7 +201,7 @@
 #         scope=prime.ScopeDefinition(model, label_expression="ground, patch*, wheel"), update=True
 #     )
 #     display.show()
-
+#
 print(model)
 
 ###############################################################################

examples/gallery/12_import_from_pygeometry.py

@@ -0,0 +1,128 @@
+# Copyright (C) 2024 - 2025 ANSYS, Inc. and/or its affiliates.
+# SPDX-License-Identifier: MIT
+#
+#
+# Permission is hereby granted, free of charge, to any person obtaining a copy
+# of this software and associated documentation files (the "Software"), to deal
+# in the Software without restriction, including without limitation the rights
+# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+# copies of the Software, and to permit persons to whom the Software is
+# furnished to do so, subject to the following conditions:
+#
+# The above copyright notice and this permission notice shall be included in all
+# copies or substantial portions of the Software.
+#
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+# SOFTWARE.
+
+"""
+.. _ref_pygeometry_import:
+
+=================================
+Import a geometry from PyGeometry
+=================================
+
+**Summary**: This example demonstrates how to create a simple geometry using
+PyGeometry and import it into PyPrimeMesh.
+
+
+
+Procedure
+~~~~~~~~~
+#. Launch the Geometry service and create a simple geometry.
+#. Launch Ansys Prime Server.
+#. Create a CAD geometry with PyGeometry.
+#. Import the geometry into PyPrimeMesh.
+#. Mesh and display the imported geometry.
+#. Exit the PyPrimeMesh and PyGeometry sessions.
+
+
+
+Create the geometry
+~~~~~~~~~~~~~~~~~~~
+First we create a simple geometry using PyGeometry. The geometry is a plate
+with a hole in the center:
+
+.. code-block:: python
+
+   from pint import Quantity
+
+   from ansys.geometry.core import launch_modeler
+   from ansys.geometry.core.math import Point2D
+   from ansys.geometry.core.misc import UNITS
+   from ansys.geometry.core.sketch import Sketch
+   from ansys.geometry.core.math import Point2D
+
+   modeler = launch_modeler(hidden=True)
+
+   sketch = Sketch()
+   (
+       sketch.segment(Point2D([-4, 5], unit=UNITS.mm), Point2D([4, 5], unit=UNITS.mm))
+       .segment_to_point(Point2D([4, -5], unit=UNITS.mm))
+       .segment_to_point(Point2D([-4, -5], unit=UNITS.mm))
+       .segment_to_point(Point2D([-4, 5], unit=UNITS.mm))
+       .box(Point2D([0, 0], unit=UNITS.mm), Quantity(3, UNITS.mm), Quantity(3, UNITS.mm))
+   )
+   design = modeler.create_design("ExtrudedPlateNoHoles")
+   body = design.extrude_sketch(f"PlateLayer", sketch, Quantity(2, UNITS.mm))
+   sketch_hole = Sketch()
+   sketch_hole.circle(Point2D([0, 0], unit=UNITS.mm), Quantity(0.5, UNITS.mm))
+
+   hole_centers = [
+       Plane(Point3D([3, 4, 0], unit=UNITS.mm)),
+       Plane(Point3D([-3, 4, 0], unit=UNITS.mm)),
+       Plane(Point3D([-3, -4, 0], unit=UNITS.mm)),
+       Plane(Point3D([3, -4, 0], unit=UNITS.mm)),
+   ]
+   for center in hole_centers:
+       sketch_hole.plane = center
+       design.extrude_sketch(
+           name=f"H_{center.origin.x}_{center.origin.y}",
+           sketch=sketch_hole,
+           distance=Quantity(2, UNITS.mm),
+           cut=True,
+       )
+
+
+Import the geometry into PyPrimeMesh
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+Once the geometry is created, we can pass the Design object to PyPrimeMesh to
+create a mesh:
+
+.. code-block:: python
+
+   import ansys.meshing.prime as prime
+   from ansys.meshing.prime.graphics.plotter import PrimePlotter
+
+   prime_client = prime.launch_prime()
+   model = prime_client.model
+   mesh_util = prime.lucid.Mesh(model=model)
+   mesh_util.from_geometry(design)
+
+
+Mesh the geometry and display it
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+With the geometry imported, we can now mesh it and display the mesh:
+
+.. code-block:: python
+
+   model.set_global_sizing_params(
+       prime.GlobalSizingParams(model, min=0.2, max=10.0, growth_rate=1.2)
+   )
+   mesh_util.surface_mesh(min_size=0.2)
+   mesh_util.volume_mesh()
+
+   # Print the results
+   model.parts[0].print_mesh = True
+   print(model)
+   display = PrimePlotter()
+   display.plot(model)
+   display.show()
+   modeler.close()
+   prime_client.exit()
+"""

pyproject.toml

@@ -32,6 +32,9 @@ dependencies = [
 graphics = [
   "ansys-tools-visualization-interface>=0.4.7,<1",
 ]
+pygeometry = [
+  "ansys-geometry-core>=0.7.6,<1",
+]
 tests = [
   "ansys-tools-visualization-interface==0.9.1",
   "pytest==8.3.5",
@@ -56,6 +59,7 @@ doc = [
 ]
 all = [
   "ansys-tools-visualization-interface>=0.2.6,<1",
+  "ansys-geometry-core>=0.7.6,<1",
 ]
 
 [project.urls]

src/ansys/meshing/prime/lucid/mesh_util.py

@@ -23,6 +23,7 @@
 """Module for meshing utility functions."""
 import enum
 import os
+import tempfile
 from typing import Iterable, List
 
 import ansys.meshing.prime as prime
@@ -31,6 +32,18 @@
 from .utils import check_name_pattern
 
 
+class ImportTypes(enum.Enum):
+    """Allowed geometry import types."""
+
+    PMDB = 0
+    FMD = 1
+    IGES = 2
+    STEP = 3
+    PARA_BIN = 4
+    PARA_TEXT = 5
+    SCDOCX = 6
+
+
 class LabelToZoneMethod(enum.IntEnum):
     """Provides for creating zones from labels."""
 
@@ -69,6 +82,33 @@ def __init__(self, model: prime.Model):
         self._model = model
         self._logger = model.python_logger
 
+    def from_geometry(
+        self, design: "ansys.geometry.core.Design", import_type: ImportTypes = ImportTypes.FMD
+    ):
+        """Import geometry from an Ansys Design object.
+
+        Parameters
+        ----------
+        design : ansys.geometry.core.Design
+            Ansys Design object to import geometry from.
+        import_type : ImportTypes, optional
+            Type of import. The default is ImportTypes.FMD.
+        """
+        with tempfile.TemporaryDirectory() as tmpdir:
+            if import_type == ImportTypes.FMD:
+                file_path = design.export_to_fmd(tmpdir)
+            elif import_type == ImportTypes.IGES:
+                file_path = design.export_to_iges(tmpdir)
+            elif import_type == ImportTypes.STEP:
+                file_path = design.export_to_step(tmpdir)
+            elif import_type == ImportTypes.PARA_BIN:
+                file_path = design.export_to_parasolid_bin(tmpdir)
+            elif import_type == ImportTypes.PARA_TEXT:
+                file_path = design.export_to_parasolid_text(tmpdir)
+            elif import_type == ImportTypes.SCDOCX:
+                file_path = design.export_to_scdocx(tmpdir)
+            self.read(str(file_path))
+
     def read(
         self, file_name: str, append: bool = False, cad_reader_route: prime.CadReaderRoute = None
     ):