How to get the grid of a scoped meshed_region by named selection?

[rmanno91]

Hi I am trying to get the grid of a meshed region scoped by name selection. Unfortunalty I was not able to find the solution yet would you be able to help?
I have uploaded a zip containg the script and a result from PyMAPDL I am trying to plot
@germa89
files.zip

[germa89]

@PProfizi this is the topic we talked about like a week ago. Could you give us a hand please?

[rmanno91]

Well we found the solution, thanks anyways :)

mesh_op = dpf.operators.mesh.from_scoping()
mesh_scoping = dpf.mesh_scoping_factory.named_selection_scoping("CM_1", model=model)
my_scoping = mesh_op.inputs.scoping.connect(mesh_scoping)
my_inclusive = int(0)
mesh_op.inputs.inclusive.connect(my_inclusive)
mesh_op.inputs.mesh.connect(model.metadata.meshed_region)

[PProfizi]

@rmanno91 Glad you found out how to split the mesh!
I also tried to make a small example on my end, plotting on all the mesh with opacity, and I think the way I plot a deformed mesh is different from what you did originally.
This is what I ended-up with:

from ansys.dpf import core as dpf
import pyvista as pv
import os
###############################################################################
# Post-processing of the results using PyDPF
# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# PyDPF can be used in this context to visualise the crack opening throughout
# the simulation as an animation. The steps for producing that are provided in
# below.

# Generating DPF model
model = dpf.Model(os.path.join(os.getcwd(), "file.rst"))
print(model)

# Getting mesh of the model
meshed_region = model.metadata.meshed_region

# Print available named selections
available_named_selections = meshed_region.available_named_selections
print(available_named_selections)
# ['BOT_NOD', 'CM_1', 'CM_2', 'TOP_NOD', '_CONTACT', '_ELEMCM', '_TARGET', '__NODE__']

# NMISC index of damage results
nmisc_index = 70

# All displacements for deformation and time_scoping on all frequencies
complete_disp_fc = model.results.displacement.on_all_time_freqs.eval()
complete_time_scoping = dpf.time_freq_scoping_factory.scoping_on_all_time_freqs(model)

# Plot displacement and damage for each named selection

for named_selection in available_named_selections:
    scoping = model.metadata.named_selection(named_selection)

    # Displacement
    disp_fc = model.results.displacement.on_all_time_freqs(mesh_scoping=scoping).eval()
    print(disp_fc[-1])
    # Filter out named_selection giving empty fields (fix required)
    if disp_fc[-1].elementary_data_count > 0:
        meshed_region.plot(field_or_fields_container=disp_fc[-1],
                           deform_by=complete_disp_fc[-1],
                           opacity=0.3,
                           text=named_selection)
    else:
        print(f"Named selection {named_selection} gave an empty displacement field.")

    # Damage
    damage_fc = dpf.operators.result.nmisc(time_scoping=complete_time_scoping,
                                           data_sources=model.metadata.data_sources,
                                           mesh_scoping=scoping,
                                           item_index=nmisc_index).eval()
    print(damage_fc[-1])
    # Filter out named_selection giving empty fields (fix required)
    if damage_fc[-1].elementary_data_count > 0:
        meshed_region.plot(damage_fc[-1],
                           deform_by=complete_disp_fc[-1],
                           show_edges=True,
                           scalar_bar_args={"title": "Damage"},
                           clim=[0, 1],
                           opacity=0.3,
                           text=named_selection)
    else:
        print(f"Named selection {named_selection} gave an empty damage field.")

# If one wants a combination of named selections, one can add the fields_containers and then plot it
# cm_1_scoping = model.metadata.named_selection("CM_1")
# cm_1_fc = model.results.displacement.on_all_time_freqs(mesh_scoping=cm_1_scoping).eval()
# cm_2_scoping = model.metadata.named_selection("CM_2")
# cm_2_fc = model.results.displacement.on_all_time_freqs(mesh_scoping=cm_2_scoping).eval()
# cm_fc = dpf.operators.add_fc(cm_1_fc, cm_2_fc).eval()

# Concerning animation of results, the next release (in the coming days) of PyDPF-Core will 
# introduce it as:
# FieldsContainer.animate(save_as="name.gif")

[PProfizi]

As explained in the comments, FieldsContainer.animate is to be released soon and will enable you to save GIFs easily.

[itanfeng]

Hi! I have a problem. I run the following code, it prints nothing but a blank list []
meshed_region.available_named_selections

I actually created several named selections in Workbench. It seems that the created named selections are not saved in the file.rst file or the file.rth file.
Really need help!
Thanks a lot!

[PProfizi]

Hello @itanfeng, if this problem is still there, please raise an issue in the Issues section of the GitHub project.