Maxwell 2D: transient winding analysis#

This example shows how you can use PyAEDT to create a project in Maxwell 2D and run a transient simulation. It runs only on Windows using CPython.

The following libraries are required for the advanced postprocessing features used in this example:

Install these libraries with:

pip install numpy pyvista matplotlib

Perform required imports#

Perform required imports.

import os
import pyaedt

Set non-graphical mode#

Set non-graphical mode. "PYAEDT_NON_GRAPHICAL" is needed to generate documentation only. You can set non_graphical either to True or False.

non_graphical = False

Insert Maxwell 2D design and save project#

Insert a Maxwell 2D design and save the project.

maxwell_2d = pyaedt.Maxwell2d(solution_type="TransientXY", specified_version="2023.1", non_graphical=non_graphical,
                              new_desktop_session=True, projectname=pyaedt.generate_unique_project_name())

Create rectangle and duplicate it#

Create a rectangle and duplicate it.

rect1 = maxwell_2d.modeler.create_rectangle([0, 0, 0], [10, 20], name="winding", matname="copper")
added = rect1.duplicate_along_line([14, 0, 0])
rect2 = maxwell_2d.modeler[added[0]]

Create air region#

Create an air region.

region = maxwell_2d.modeler.create_region([100, 100, 100, 100, 100, 100])

Assign windings and balloon#

Assigns windings to the sheets and a balloon to the air region.

maxwell_2d.assign_winding([,], name="PHA")
<pyaedt.modules.Boundary.BoundaryObject object at 0x000001BC1F63F640>

Plot model#

Plot the model.

maxwell_2d.plot(show=False, export_path=os.path.join(maxwell_2d.working_directory, "Image.jpg"), plot_air_objects=True)
Maxwell2D Transient
<pyaedt.generic.plot.ModelPlotter object at 0x000001BC21FDAC70>

Create setup#

Create the transient setup.

setup = maxwell_2d.create_setup()
setup.props["StopTime"] = "0.02s"
setup.props["TimeStep"] = "0.0002s"
setup.props["SaveFieldsType"] = "Every N Steps"
setup.props["N Steps"] = "1"
setup.props["Steps From"] = "0s"
setup.props["Steps To"] = "0.002s"

Create rectangular plot#

Create a rectangular plot.
    "InputCurrent(PHA)", domain="Time", primary_sweep_variable="Time", plotname="Winding Plot 1"
<pyaedt.modules.report_templates.Standard object at 0x000001BC21398730>

Solve model#

Solve the model.


Create output and plot using PyVista#

Create the output and plot it using PyVista.

cutlist = ["Global:XY"]
face_lists = rect1.faces
face_lists += rect2.faces
timesteps = [str(i * 2e-4) + "s" for i in range(11)]
id_list = [ for f in face_lists]
animatedGif =
    intrinsic_dict={"Time": "0s"},
animatedGif.isometric_view = False
animatedGif.camera_position = [15, 15, 80]
animatedGif.focal_point = [15, 15, 0]
animatedGif.roll_angle = 0
animatedGif.elevation_angle = 0
animatedGif.azimuth_angle = 0
# Set off_screen to False to visualize the animation.
# animatedGif.off_screen = False
Maxwell2D Transient

Generate plot outside of AEDT#

Generate the same plot outside AEDT.

solutions ="InputCurrent(PHA)", primary_sweep_variable="Time")
Simulation Results Plot
<Figure size 2000x1000 with 1 Axes>

Close AEDT#

Close AEDT.


Total running time of the script: ( 1 minutes 4.617 seconds)

Gallery generated by Sphinx-Gallery