.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples\10-mesh_operations\05-skin_extraction.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_10-mesh_operations_05-skin_extraction.py: .. _ref_skin_mesh: Extract the skin from a mesh ~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Extracting the skin of a mesh to reduce the amount of data to operate on can be useful for specific results and for performance. .. GENERATED FROM PYTHON SOURCE LINES 11-17 .. code-block:: Python # Import necessary modules from ansys.dpf import core as dpf from ansys.dpf.core import examples from ansys.dpf.core import operators as ops .. GENERATED FROM PYTHON SOURCE LINES 18-20 Create a model object to establish a connection with an example result file and then extract: .. GENERATED FROM PYTHON SOURCE LINES 20-23 .. code-block:: Python model = dpf.Model(examples.download_multi_stage_cyclic_result()) print(model) .. rst-class:: sphx-glr-script-out .. code-block:: none DPF Model ------------------------------ Modal analysis Unit system: MKS: m, kg, N, s, V, A, degC Physics Type: Mechanical Available results: - displacement: Nodal Displacement - stress: ElementalNodal Stress - elastic_strain: ElementalNodal Strain - element_euler_angles: ElementalNodal Element Euler Angles - structural_temperature: ElementalNodal Structural temperature ------------------------------ DPF Meshed Region: 3595 nodes 1557 elements Unit: m With solid (3D) elements ------------------------------ DPF Time/Freq Support: Number of sets: 6 Cumulative Frequency (Hz) LoadStep Substep Harmonic index 1 188.385357 1 1 0.000000 2 325.126418 1 2 0.000000 3 595.320548 1 3 0.000000 4 638.189511 1 4 0.000000 5 775.669703 1 5 0.000000 6 928.278013 1 6 0.000000 .. GENERATED FROM PYTHON SOURCE LINES 24-28 Create the workflow ~~~~~~~~~~~~~~~~~~~~ Maximum principal stress usually occurs on the skin of the model. Computing results only on this skin reduces the data size. .. GENERATED FROM PYTHON SOURCE LINES 28-35 .. code-block:: Python # Create a simple workflow computing the principal stress on the skin # of the model. skin_op = ops.mesh.external_layer(model.metadata.meshed_region) skin_mesh = skin_op.outputs.mesh() .. GENERATED FROM PYTHON SOURCE LINES 36-37 Plot the mesh skin: .. GENERATED FROM PYTHON SOURCE LINES 37-39 .. code-block:: Python skin_mesh.plot() .. image-sg:: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_001.png :alt: 05 skin extraction :srcset: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_001.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 40-41 Compute the stress principal invariants on the skin nodes only: .. GENERATED FROM PYTHON SOURCE LINES 41-50 .. code-block:: Python stress_op = ops.result.stress(data_sources=model.metadata.data_sources) stress_op.inputs.requested_location.connect(dpf.locations.nodal) stress_op.inputs.mesh_scoping.connect(skin_op.outputs.nodes_mesh_scoping) principal_op = ops.invariant.principal_invariants_fc(stress_op) principal_stress_1 = principal_op.outputs.fields_eig_1()[0] principal_stress_2 = principal_op.outputs.fields_eig_2()[0] principal_stress_3 = principal_op.outputs.fields_eig_3()[0] .. GENERATED FROM PYTHON SOURCE LINES 51-53 Manipulate data locally ~~~~~~~~~~~~~~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 56-59 This example goes over the fields, keeping the largest invariant value by node if the averaged value of invariants is large enough. Exploring data allows you to customize it to meet your needs. .. GENERATED FROM PYTHON SOURCE LINES 59-79 .. code-block:: Python node_scoping_ids = principal_stress_1.scoping.ids threshold = 300000.0 field_to_keep = dpf.fields_factory.create_scalar_field( len(node_scoping_ids), location=dpf.locations.nodal ) with field_to_keep.as_local_field() as f: with principal_stress_1.as_local_field() as s1: with principal_stress_2.as_local_field() as s2: with principal_stress_3.as_local_field() as s3: for i, id in enumerate(node_scoping_ids): d1 = abs(s1.get_entity_data_by_id(id)) d2 = abs(s2.get_entity_data_by_id(id)) d3 = abs(s3.get_entity_data_by_id(id)) if (d1 + d2 + d3) / 3.0 > threshold: d = max(d1, d2, d3) f.append(d, id) .. GENERATED FROM PYTHON SOURCE LINES 80-82 Plot result field ~~~~~~~~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 85-86 Plot the result field on the skin mesh: .. GENERATED FROM PYTHON SOURCE LINES 86-88 .. code-block:: Python skin_mesh.plot(field_to_keep) .. image-sg:: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_002.png :alt: 05 skin extraction :srcset: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_002.png :class: sphx-glr-single-img .. GENERATED FROM PYTHON SOURCE LINES 89-91 Plot initial invariants ~~~~~~~~~~~~~~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 94-95 Plot the initial invariants on the skin mesh: .. GENERATED FROM PYTHON SOURCE LINES 95-99 .. code-block:: Python skin_mesh.plot(principal_stress_1) skin_mesh.plot(principal_stress_2) skin_mesh.plot(principal_stress_3) .. rst-class:: sphx-glr-horizontal * .. image-sg:: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_003.png :alt: 05 skin extraction :srcset: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_003.png :class: sphx-glr-multi-img * .. image-sg:: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_004.png :alt: 05 skin extraction :srcset: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_004.png :class: sphx-glr-multi-img * .. image-sg:: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_005.png :alt: 05 skin extraction :srcset: /examples/10-mesh_operations/images/sphx_glr_05-skin_extraction_005.png :class: sphx-glr-multi-img .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 3.330 seconds) .. _sphx_glr_download_examples_10-mesh_operations_05-skin_extraction.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: 05-skin_extraction.ipynb <05-skin_extraction.ipynb>` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: 05-skin_extraction.py <05-skin_extraction.py>` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_