Script for creating HAZ layers without using HAZ named selection creation option in GUI

kiran.sai

I have set a smaller mesh size and I don't want to create HAZ layers that disturb the mesh near the weld region. The number of HAZ layers I can create is 6 as per the current version, so it's not possible for me to take more HAZ layers in order to capture my interest region with the smaller mesh element size

kiran.sai

You can use the following script for this-

import math

SEARCH_DISTANCE = 10e-3  # Distance in m

def get_node_coords(mesh_data, node_id):

    node = mesh_data.NodeById(node_id)

    return (node.X, node.Y, node.Z)

mesh_data=Model.Analyses[0].MeshData

all_bodies = DataModel.GetObjectsByType(DataModelObjectCategory.Body)

weld_bodies = [b for b in all_bodies if b.Name.lower().startswith("weld")]

other_bodies = [b for b in all_bodies if not b.Name.lower().startswith("weld")]

for weld_body in weld_bodies:

    weld_geo = weld_body.GetGeoBody()

    weld_node_ids = set(mesh_data.MeshRegionById(weld_geo.Id).NodeIds)

    neighbor_count = 0

    for other_body in other_bodies:

        other_geo = other_body.GetGeoBody()

        other_node_ids = set(mesh_data.MeshRegionById(other_geo.Id).NodeIds)

        shared_node_ids = weld_node_ids.intersection(other_node_ids)

        if shared_node_ids:

            neighbor_count += 1

            shared_coords = [get_node_coords(mesh_data, nid) for nid in shared_node_ids]

            target_element_ids = []

            other_elem_ids = mesh_data.MeshRegionById(other_geo.Id).ElementIds

            for e_id in other_elem_ids:

                element = mesh_data.ElementById(e_id)

                e_nodes = set(element.NodeIds)

                if not e_nodes.isdisjoint(shared_node_ids):

                    target_element_ids.append(e_id)

                    continue

                ex, ey, ez = element.Centroid[0], element.Centroid[1], element.Centroid[2]

                for sx, sy, sz in shared_coords:

                    dist = math.sqrt((ex-sx)**2 + (ey-sy)**2 + (ez-sz)**2)

                    if dist <= SEARCH_DISTANCE:

                        target_element_ids.append(e_id)

                        break

            if target_element_ids:

                sel_name = "HAZ_Zone_{0}_on_{1}".format(weld_body.Name, other_body.Name)

                selection = ExtAPI.SelectionManager.CreateSelectionInfo(SelectionTypeEnum.MeshElements)

                selection.Ids = target_element_ids

                ns = Model.AddNamedSelection()

                ns.Name = sel_name

                ns.Location = selection