How to calculate contact surface and heat flow with the use of PyDpf

Kleanthis · August 9

Kleanthis · August 9

The script below captures the contact for each element by getting the CONTAC174 element and calculates its contact surface. Afterwards it calculates the heat flux of each element.

A prerequisite for the generation of the contact surface named selection is to add the following APDL command under the desired contact:

my_cont=cid

And the following command object before solution:

/PREP7

ESEL, S, TYPE, , my_cont
CM, PCB_CONT, ELEM
ALLSEL

/SOLU

The PyDPF script is as follows:

from ansys.dpf import core as dpf
from ansys.dpf import core as dpf
from ansys.dpf.core.plotter import DpfPlotter

import numpy as np

"""
Load result file and create model database
"""
ds = dpf.DataSources("file.rth")
model = dpf.Model(ds)

# %%
get_NSs = model.metadata.available_named_selections
print(get_NSs)

# %%
mesh_scoping = model.metadata.named_selection("PCB_CONT_2")
print(mesh_scoping)

# %%
print(get_NSs)

mesh_scoping = model.metadata.named_selection("PCB_CONT_2")
print(mesh_scoping)

# %%
contact_status_op = dpf.operators.result.nmisc(
    # time_scoping=-1,
    mesh_scoping=mesh_scoping,
    # fields_container=my_fields_container,
    # streams_container=my_streams_container,
    data_sources=ds,
    # mesh=my_mesh,
    item_index=41,
    # num_components=my_num_components,
)

print(contact_status_op)

# %%
contact_status_res = contact_status_op.outputs.fields_container()
print(contact_status_res)

# %%
scoping_op = dpf.operators.mesh.from_scoping()
scoping_op.inputs.scoping.connect(mesh_scoping)
scoping_op.inputs.mesh.connect(contact_status_res[0].meshed_region)

my_mesh = scoping_op.outputs.mesh()
print(my_mesh)

# %%
plot = DpfPlotter()
plot.add_field(contact_status_res[0], my_mesh, notebook=False)
plot.show_figure(show_axes=True)

# %%
touch_cond = (contact_status_res[0].data >= 2.)
Nelems_touch = sum(touch_cond)

print("Number of Touching Elements: ", Nelems_touch)
print("Number of Total Elements: ", len(contact_status_res[0].data))

# %%
contact_area_op = dpf.operators.result.nmisc(
    # time_scoping=-1,
    mesh_scoping=mesh_scoping,
    # fields_container=my_fields_container,
    # streams_container=my_streams_container,
    data_sources=ds,
    # mesh=my_mesh,
    item_index=58,
    # num_components=my_num_components,
)

print(contact_area_op)

# %%
contact_area_res = contact_area_op.outputs.fields_container()
print(contact_area_res)

# %%
print(contact_area_res[0].data)

# %%
scoping_op = dpf.operators.mesh.from_scoping()
scoping_op.inputs.scoping.connect(mesh_scoping)
scoping_op.inputs.mesh.connect(contact_area_res[0].meshed_region)

my_mesh = scoping_op.outputs.mesh()
print(my_mesh)

# %%
plot = DpfPlotter()
plot.add_field(contact_area_res[0], my_mesh, show_max=True, show_min=True, notebook=False)
plot.show_figure(show_axes=True)

# %%
contact_heatFlux_op = dpf.operators.result.smisc(
    # time_scoping=-1,
    mesh_scoping=mesh_scoping,
    # fields_container=my_fields_container,
    # streams_container=my_streams_container,
    data_sources=ds,
    # mesh=my_mesh,
    item_index=14,
    # num_components=my_num_components,
)

print(contact_heatFlux_op)

# %%
contact_heatFlux_res = contact_heatFlux_op.outputs.fields_container()
print(contact_heatFlux_res)

# %%
scoping_op = dpf.operators.mesh.from_scoping()
scoping_op.inputs.scoping.connect(mesh_scoping)
scoping_op.inputs.mesh.connect(contact_heatFlux_res[0].meshed_region)

my_mesh = scoping_op.outputs.mesh()
print(my_mesh)

# %%
plot = DpfPlotter()
plot.add_field(contact_heatFlux_res[0], my_mesh, show_max=True, show_min=True, notebook=False)
plot.show_figure(show_axes=True)

How to calculate contact surface and heat flow with the use of PyDpf

Comments

Categories