Postprocessing automation in Icepak

[mjmashayekhi]

Hi,

I'm trying to develop a pyAEDT script to automate Icepak simulation for certain models. My intention is to run the model for different power models and postprocess the results by obtaining the max temperature of an object called "Box1" in the Icepak model.

I used the following script and several variations of this script, but always getting error messages for the postprocessing portion.

[mjmashayekhi]

Script:

import os
from pyaedt import Icepak
from openpyxl import Workbook
import matplotlib.pyplot as plt

#Define the AEDT project and design
project_path = r"E:\Automation\AEDT parametric\20240821-AEDT_Automation.aedt"
design_name = "Design1"
analysis_setup = "Setup2"
output_path = r"E:\Automation\AEDT parametric\output.xlsx"

#List of power values to sweep for $Pbox
pbox_values = [1.0, 2.0, 3.0] # Modify as needed

#Create an Excel workbook to save results
wb = Workbook()
ws = wb.active
ws.title = "Temperature Data"
ws.append(["$Pbox (W)", "Max Temperature on Box1 (°C)"])

#List to store max temperatures for plotting
max_temperatures = []

try:
with Icepak(projectname=project_path,
designname=design_name,
specified_version="2024.2",
non_graphical=False,
new_desktop_session=True) as icepak:
print("Available objects:", icepak.modeler.user_defined_component_names)
print("Native object names:", icepak.modeler.object_names)
icepak.analysis_setup = analysis_setup

    for pbox in pbox_values:
        print(f"\nRunning simulation with $Pbox = {pbox} W...")

        #Update the $Pbox parameter
        icepak["$Pbox"] = pbox

        #Solve
        icepak.analyze()

        #Create a report to extract temperature data
        field_data = icepak.post.create_report(
            expressions=["Temperature"],
            primary_sweep_variable="Distance",
            report_category="Fields",
            polyline_points=500,
        )

        #Extract the maximum temperature from the report data
        #Use the correct method to access the data
        report_data = field_data.get_solution_data()  # Hypothetical method
        distance    = [k[0] for k, _ in report_data.full_matrix_mag_phase[0]["Temperature"].items()]
        temperature_values = [v for _, v in report_data.full_matrix_mag_phase[0]["Temperature"].items()]
        max_temp = max(temperature_values)
        max_temperatures.append(max_temp)

        #Append results to Excel
        ws.append([pbox, max_temp])

#Save the Excel workbook
wb.save(output_path)

#Plot the results
plt.figure()
plt.plot(pbox_values, max_temperatures, marker='o')
plt.title('Max Temperature on Box1 vs. $Pbox')
plt.xlabel('$Pbox (W)')
plt.ylabel('Max Temperature on Box1 (°C)')
plt.grid(True)
plt.show()

except Exception as e:
print(f"An error occurred: {e}")

[Samuelopez-ansys]

Hi @mjmashayekhi ,

Please debug the code in order to find the issue. Does field_data return an object? If not, you probably should also pass the setup to create_report.

Here you can find more examples: https://examples.aedt.docs.pyansys.com/version/dev/examples/low_frequency/multiphysics/maxwell_icepak.html

Please share all the files needed to reproduce the issue, because from my side, it works.

Thanks

CC: @lorenzovecchietti @siva-krishnaswamy

[mjmashayekhi]

Hi Samuel,

Thank you for your feedback.
I tried multiple different approaches to obtain the max temperature on the faces of an object in my Icepak model and each of them give me different errors.

• For the code I shared here, I get the following error:
  AttributeError: 'bool' object has no attribute 'full_matrix_mag_phase'

• I also tried "temperature = icepak.post.get_face_temperature(faces)"
  - I get this error: "AttributeError: 'PostProcessor' object has no attribute 'get_face_temperature'"

• I tried:
  for fid in face_ids:
  # Create a temporary field report for the face
  report_name = f"Temp_Face_{fid}"
  icepak.post.create_report(
  expressions="Temperature",
  domain="Faces",
  context={"Faces": [fid]},
  plotname=report_name,
  plot_type="Field"
  )
  And I got this error message:
  PyAEDT ERROR: AEDT GRPC API CALL ERROR on create
  PyAEDT ERROR: Failed to execute grpc AEDT command: CreateReport

20240821-AEDT_parametric_automation.zip

Thanks,
Mohammad

[lorenzovecchietti]

@mjmashayekhi

I hope I haven’t missed any critical details that would make this approach unsuitable for your use case. Based on what you're aiming to achieve, I believe it can be done quite simply using the get_temperature_extremum method.

Here’s a quick example:

ipk = ansys.aedt.core.Icepak()
out = ipk.post.get_temperature_extremum("Box1", "Max", "Surface")

# out[0] gives the (x, y, z) position of the maximum temperature
# out[1] gives the maximum temperature value

Let me know if this works for you!

[mjmashayekhi]

Hi @lorenzovecchietti, Thank you for your response. I'm not able to find this attribute in the icepak.post. Is it specific to a pyAEDT vertsion?

Thanks

[lorenzovecchietti]

https://github.com/ansys/pyaedt/blame/fea45cf6a73f8346c1d6e7fa01e81faae6439bbd/src/ansys/aedt/core/visualization/post/post_icepak.py#L404

It was added about six months ago, so it's been part of quite a few releases now.