How to calculate contact surface and heat flow with the use of PyDpf
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)
