How to make mapped mesh?

[SadBoen]

Hello, Dear all,
In Workbench, I could use surface meshing function to generate very pretty cells as following pic:

and I found the key setting is 'mapped mesh' setting to yes.
But in Pyprimemesh, Does anyone know how to handle it? I have tried many times, my code and 3D model are as follows：

import os
import ansys.meshing.prime as prime
from ansys.meshing.prime.graphics.plotter import PrimePlotter

prime_client = prime.launch_prime(ip="127.0.0.1", port=50055)
pmodel = prime_client.model

mesh_util = prime.lucid.Mesh(model=pmodel)
step_file_path = os.path.join(os.getcwd(), "7145521xx-lw.scdocx")
mesh_util.read(file_name=step_file_path, cad_reader_route=prime.CadReaderRoute.SPACECLAIM)

mesh_util = prime.lucid.Mesh(model=pmodel)

scope_for_key_faces = prime.lucid.SurfaceScope(
part_expression="*",
entity_expression="m-lw-s",
scope_evaluation_type=prime.ScopeEvaluationType.LABELS
)

size_control_for_lw = pmodel.control_data.create_size_control(sizing_type=prime.SizingType.SOFT)
size_control_for_lw.set_soft_sizing_params(params=prime.SoftSizingParams(model=pmodel,
max=1.0, growth_rate=1.0))

surf_mesh_results = mesh_util.surface_mesh_with_size_controls(
size_control_names=size_control_for_lw,
generate_quads=True,
scope=scope_for_key_faces
)

vol_mesh_results = mesh_util.volume_mesh(
volume_fill_type=prime.VolumeFillType.POLY
)

display = PrimePlotter()
display.plot(pmodel)
display.show()

prime_client.exit()
7145521xx-lw.zip

I hope to receive guidance from everyone, thank you very much!

[waltersma]

Hi @SadBoen ,

I assume that the result you are aiming for is a 100% hexahedral mesh rather than the surface mesh alone being quad? The reason I ask is that in your example script you specify polyhedral elements. There is not currently a quad map control available for surface meshing directly but if you make use of the multizone method you can achieve a 100% hex mesh as you show in the image.

from ansys.meshing import prime

prime_client = prime.launch_prime()
model = prime_client.model
mesh_util = prime.lucid.Mesh(model=model)

geom = r"7145521xx-lw.scdocx"

mesh_util.read(geom)
print(model)
part = model.parts[0]

sizing_params = prime.GlobalSizingParams(model=model, min=1.0, max=1.0)
model.set_global_sizing_params(params=sizing_params)

multi_zone_scope_volume = prime.ScopeDefinition(
    model=model,
    entity_type=prime.ScopeEntity.VOLUME,
    zone_expression="*",
    evaluation_type=prime.ScopeEvaluationType.ZONES,
)
multi_zone_scope_surface = prime.ScopeDefinition(
    model=model, entity_type=prime.ScopeEntity.FACEZONELETS, zone_expression="*"
)

multi_zone_sizing_params = prime.MultiZoneSizingParams(
    model=model,
    min_size=1.0,
    max_size=1.0,
)
multi_zone_sweep_params = prime.MultiZoneSweepMeshParams(
    model=model,
    sweep_type=prime.SweepType.STANDARD,
)

multi_zone_control = model.control_data.create_multi_zone_control()
multi_zone_control.set_volume_scope(multi_zone_scope_volume)
multi_zone_control.set_surface_scope(multi_zone_scope_surface)
multi_zone_control.set_multi_zone_sizing_params(multi_zone_sizing_params)
multi_zone_control.set_sweep_mesh_params(multi_zone_sweep_params)
multi_zone_control.set_map_mesh_params(
    prime.MultiZoneMapMeshParams(model, scope=multi_zone_scope_surface)
)

auto_mesh = prime.AutoMesh(model)
auto_mesh_params = prime.AutoMeshParams(model)
auto_mesh_params.multi_zone_control_ids = [multi_zone_control.id]

result = auto_mesh.mesh(part.id, auto_mesh_params)

I hope this helps!