.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/modeling_features/030-imported-plies.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_modeling_features_030-imported-plies.py: .. _imported_plies_example: Imported ply ============ The definition and use of imported plies is demonstrated in this example. In a nutshell, the difference between :class:`.ImportedModelingPly` and :class:`.ModelingPly` is that the surface mesh of an :class:`.ImportedModelingPly` is defined by an external source, such as a CAD surface, where a :class:`.ModelingPly` is always defined on the initial loaded shell mesh. Therefore, an imported ply can only be used in combination with an :class:`.ImportedSolidModel`. This examples shows hot to: - Load an initial mesh - Add a :class:`.Material` and :class:`.Fabric` - Import geometries which will be used to define the surface of the :class:`.ImportedModelingPly` - Add two imported modeling plies - Create an :class:`.ImportedSolidModel` - Map the imported plies to the solid model - Visualized the mapped plies. .. GENERATED FROM PYTHON SOURCE LINES 49-53 Import modules -------------- Import the standard library and third-party dependencies. .. GENERATED FROM PYTHON SOURCE LINES 53-59 .. code-block:: Python import os import pathlib import tempfile import pyvista .. GENERATED FROM PYTHON SOURCE LINES 60-61 Import the PyACP dependencies. .. GENERATED FROM PYTHON SOURCE LINES 61-69 .. code-block:: Python from ansys.acp.core import CADGeometry, ImportedPlyOffsetType, PlyType, VirtualGeometry, launch_acp from ansys.acp.core.extras import ExampleKeys, get_example_file from ansys.acp.core.material_property_sets import ConstantDensity, ConstantEngineeringConstants CAMERA_POSITION = [(0.0436, 0.0102, 0.0193), (0.0111, 0.0035, 0.0046), (-0.1685, 0.9827, -0.0773)] .. GENERATED FROM PYTHON SOURCE LINES 71-73 Start ACP and load the model ---------------------------- .. GENERATED FROM PYTHON SOURCE LINES 75-76 Get the example file from the server. .. GENERATED FROM PYTHON SOURCE LINES 76-80 .. code-block:: Python tempdir = tempfile.TemporaryDirectory() WORKING_DIR = pathlib.Path(tempdir.name) input_file = get_example_file(ExampleKeys.BASIC_FLAT_PLATE_DAT, WORKING_DIR) .. GENERATED FROM PYTHON SOURCE LINES 81-82 Launch the PyACP server and connect to it. .. GENERATED FROM PYTHON SOURCE LINES 82-84 .. code-block:: Python acp = launch_acp() .. GENERATED FROM PYTHON SOURCE LINES 85-86 Create a model by loading a shell mesh .. GENERATED FROM PYTHON SOURCE LINES 86-88 .. code-block:: Python model = acp.import_model(path=input_file, format="ansys:dat") .. GENERATED FROM PYTHON SOURCE LINES 89-91 Add a material and a fabric with 1mm thickness. The fabric is used for the imported modeling ply. .. GENERATED FROM PYTHON SOURCE LINES 91-122 .. code-block:: Python engineering_constants_ud = ConstantEngineeringConstants.from_orthotropic_constants( E1=5e10, E2=1e10, E3=1e10, nu12=0.28, nu13=0.28, nu23=0.3, G12=5e9, G23=4e9, G31=4e9 ) density_ud = ConstantDensity(rho=2700) ud_material = model.create_material( name="E-Glass UD", ply_type=PlyType.REGULAR, engineering_constants=engineering_constants_ud, density=density_ud, ) engineering_resin = ConstantEngineeringConstants.from_isotropic_constants(E=5e9, nu=0.3) density_resin = ConstantDensity(rho=1200) void_material = model.create_material( name="Void material", ply_type=PlyType.ISOTROPIC, engineering_constants=engineering_resin, density=density_resin, ) filler_material = model.create_material( name="Filler material", ply_type=PlyType.ISOTROPIC, engineering_constants=engineering_resin, density=density_resin, ) fabric = model.create_fabric(name="E-Glass Fabric", material=ud_material, thickness=0.001) .. GENERATED FROM PYTHON SOURCE LINES 123-127 Import CAD geometries --------------------- Import two cad surfaces to define the surface of the imported modeling plies. .. GENERATED FROM PYTHON SOURCE LINES 127-146 .. code-block:: Python def create_virtual_geometry_from_file( example_key: ExampleKeys, ) -> tuple[CADGeometry, VirtualGeometry]: """Create a CAD geometry and virtual geometry.""" geometry_file = get_example_file(example_key, WORKING_DIR) geometry_obj = model.create_cad_geometry() geometry_obj.refresh(geometry_file) # upload and load the geometry file model.update() virtual_geometry = model.create_virtual_geometry( name=os.path.basename(geometry_file), cad_components=geometry_obj.root_shapes.values() ) return geometry_obj, virtual_geometry triangle_surf_cad, triangle_surf_vcad = create_virtual_geometry_from_file( ExampleKeys.RULE_GEOMETRY_TRIANGLE ) top_surf_cad, top_surf_vcad = create_virtual_geometry_from_file(ExampleKeys.SNAP_TO_GEOMETRY) .. GENERATED FROM PYTHON SOURCE LINES 147-149 Definition of Imported Plies ---------------------------- .. GENERATED FROM PYTHON SOURCE LINES 149-170 .. code-block:: Python imported_ply_group = model.create_imported_modeling_group(name="Imported Ply Group") imported_ply_triangle = imported_ply_group.create_imported_modeling_ply( name="Triangle Ply", offset_type=ImportedPlyOffsetType.BOTTOM_OFFSET, ply_material=fabric, mesh_geometry=triangle_surf_vcad, ply_angle=0, rosettes=[model.rosettes["12"]], ) imported_ply_top = imported_ply_group.create_imported_modeling_ply( name="Triangle Ply", offset_type=ImportedPlyOffsetType.MIDDLE_OFFSET, ply_material=fabric, mesh_geometry=top_surf_vcad, ply_angle=45, rosettes=[model.rosettes["12"]], ) model.update() .. GENERATED FROM PYTHON SOURCE LINES 171-175 Imported plies cannot be visualized directly yet but the cad geometries are shown here instead. To visualize the imported plies, you can save the model and load it in ACP standalone. .. GENERATED FROM PYTHON SOURCE LINES 175-191 .. code-block:: Python def plotter_with_all_geometries(cad_geometries): colors = ["green", "yellow", "blue", "red"] plotter = pyvista.Plotter() for index, cad in enumerate(cad_geometries): geom_mesh = cad.visualization_mesh.to_pyvista() plotter.add_mesh(geom_mesh, color=colors[index], opacity=0.1) edges = geom_mesh.extract_feature_edges() plotter.add_mesh(edges, color="white", line_width=4) plotter.add_mesh(edges, color="black", line_width=2) plotter.camera_position = CAMERA_POSITION return plotter plotter = plotter_with_all_geometries([triangle_surf_cad, top_surf_cad]) plotter.show() .. tab-set:: .. tab-item:: Static Scene .. image-sg:: /examples/modeling_features/images/sphx_glr_030-imported-plies_001.png :alt: 030 imported plies :srcset: /examples/modeling_features/images/sphx_glr_030-imported-plies_001.png :class: sphx-glr-single-img .. tab-item:: Interactive Scene .. offlineviewer:: /home/runner/work/pyacp/pyacp/doc/source/examples/modeling_features/images/sphx_glr_030-imported-plies_001.vtksz .. GENERATED FROM PYTHON SOURCE LINES 192-198 Map Imported Plies onto a solid mesh ------------------------------------ An external solid mesh is loaded now to map the imported plies onto the solid model. The next figure shows the imported solid mesh and the imported plies. .. GENERATED FROM PYTHON SOURCE LINES 198-210 .. code-block:: Python local_solid_mesh_file = get_example_file(ExampleKeys.BASIC_FLAT_PLATE_SOLID_MESH_CDB, WORKING_DIR) remote_solid_mesh_file = acp.upload_file(local_solid_mesh_file) imported_solid_model = model.create_imported_solid_model( name="Imported Solid Model", external_path=remote_solid_mesh_file, format="ansys:cdb", ) imported_solid_model.import_initial_mesh() plotter = plotter_with_all_geometries([triangle_surf_cad, top_surf_cad]) plotter.add_mesh(imported_solid_model.solid_mesh.to_pyvista(), show_edges=True, opacity=0.5) plotter.show() .. tab-set:: .. tab-item:: Static Scene .. image-sg:: /examples/modeling_features/images/sphx_glr_030-imported-plies_002.png :alt: 030 imported plies :srcset: /examples/modeling_features/images/sphx_glr_030-imported-plies_002.png :class: sphx-glr-single-img .. tab-item:: Interactive Scene .. offlineviewer:: /home/runner/work/pyacp/pyacp/doc/source/examples/modeling_features/images/sphx_glr_030-imported-plies_002.vtksz .. GENERATED FROM PYTHON SOURCE LINES 211-216 Add a mapping object to link the imported plies with the solid model. In this example, all imported plies are mapped in one go. The remaining elemental volume and elements which do not intersect with the imported plies are filled with a void and filler material, respectively. .. GENERATED FROM PYTHON SOURCE LINES 216-229 .. code-block:: Python imported_solid_model.create_layup_mapping_object( name="Map imported plies", use_imported_plies=True, # enable imported plies select_all_plies=True, # select all plies entire_solid_mesh=True, scale_ply_thicknesses=False, void_material=void_material, delete_lost_elements=False, filler_material=filler_material, rosettes=[model.rosettes["12"]], ) model.update() .. GENERATED FROM PYTHON SOURCE LINES 230-234 Show the imported ply geometries and mapped plies on the solid model. Note that the analysis plies are not yet directly accessible via the API of the imported solid model. Also, elemental data such as thicknesses are not yet implemented for imported plies. .. GENERATED FROM PYTHON SOURCE LINES 234-242 .. code-block:: Python plotter = plotter_with_all_geometries([triangle_surf_cad, top_surf_cad]) for imported_ply in [imported_ply_triangle, imported_ply_top]: for pp in imported_ply.imported_production_plies.values(): for ap in pp.imported_analysis_plies.values(): plotter.add_mesh(ap.solid_mesh.to_pyvista(), show_edges=True, opacity=1) plotter.add_mesh(mesh=imported_solid_model.solid_mesh.to_pyvista(), show_edges=False, opacity=0.2) plotter.show() .. tab-set:: .. tab-item:: Static Scene .. image-sg:: /examples/modeling_features/images/sphx_glr_030-imported-plies_003.png :alt: 030 imported plies :srcset: /examples/modeling_features/images/sphx_glr_030-imported-plies_003.png :class: sphx-glr-single-img .. tab-item:: Interactive Scene .. offlineviewer:: /home/runner/work/pyacp/pyacp/doc/source/examples/modeling_features/images/sphx_glr_030-imported-plies_003.vtksz .. GENERATED FROM PYTHON SOURCE LINES 243-246 The imported solid model can be passed to Mechanical or MAPDL to run an analysis as shown in the examples :ref:`pymechanical_solid_example` and :ref:`pymapdl_workflow_example`. .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 6.284 seconds) .. _sphx_glr_download_examples_modeling_features_030-imported-plies.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: 030-imported-plies.ipynb <030-imported-plies.ipynb>` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: 030-imported-plies.py <030-imported-plies.py>` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: 030-imported-plies.zip <030-imported-plies.zip>` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_