Note
Go to the end to download the full example code.
Maxwell 2D Electrostatic analysis#
This example shows how you can use PyAEDT to create a Maxwell 2D electrostatic analysis. It shows how to create the geometry, load material properties from an Excel file and set up the mesh settings. Moreover, it focuses on post-processing operations, in particular how to plot field line traces, relevant for an electrostatic analysis.
Perform required imports#
Perform required imports.
import pyaedt
Set AEDT version#
Set AEDT version.
aedt_version = "2024.1"
Initialize Maxwell 2D#
Initialize Maxwell 2D, providing the version, path to the project, and the design name and type.
setup_name = 'MySetupAuto'
solver = 'Electrostatic'
design_name = 'Design1'
project_name = pyaedt.generate_unique_project_name()
non_graphical = False
Download .xlsx file#
Set local temporary folder to export the .xlsx file to.
file_name_xlsx = pyaedt.downloads.download_file("field_line_traces", "my_copper.xlsx")
Initialize dictionaries#
Initialize dictionaries that contain all the definitions for the design variables.
geom_params_circle = {
'circle_x0': '-10mm',
'circle_y0': '0mm',
'circle_z0': '0mm',
'circle_axis': 'Z',
'circle_radius': '1mm'
}
geom_params_rectangle = {
'r_x0': '1mm',
'r_y0': '5mm',
'r_z0': '0mm',
'r_axis': 'Z',
'r_dx': '-1mm',
'r_dy': '-10mm'
}
Launch Maxwell 2D#
Launch Maxwell 2D and save the project.
M2D = pyaedt.Maxwell2d(projectname=project_name,
specified_version=aedt_version,
designname=design_name,
solution_type=solver,
new_desktop_session=True,
non_graphical=non_graphical
)
C:\actions-runner\_work\_tool\Python\3.10.9\x64\lib\subprocess.py:1072: ResourceWarning: subprocess 6756 is still running
_warn("subprocess %s is still running" % self.pid,
C:\actions-runner\_work\pyaedt\pyaedt\.venv\lib\site-packages\pyaedt\generic\settings.py:383: ResourceWarning: unclosed file <_io.TextIOWrapper name='D:\\Temp\\pyaedt_ansys.log' mode='a' encoding='cp1252'>
self._logger = val
Create object to access 2D modeler#
Create the object mod2D to access the 2D modeler easily.
mod2D = M2D.modeler
mod2D.delete()
mod2D.model_units = "mm"
Define variables from dictionaries#
Define design variables from the created dictionaries.
for k, v in geom_params_circle.items():
M2D[k] = v
for k, v in geom_params_rectangle.items():
M2D[k] = v
Read materials from .xslx file#
Read materials from .xslx file into and set into design.
mats = M2D.materials.import_materials_from_excel(file_name_xlsx)
C:\actions-runner\_work\pyaedt\pyaedt\.venv\lib\site-packages\pyaedt\modules\MaterialLib.py:907: FutureWarning: Series.__getitem__ treating keys as positions is deprecated. In a future version, integer keys will always be treated as labels (consistent with DataFrame behavior). To access a value by position, use `ser.iloc[pos]`
and val[keys.index(prop)]
C:\actions-runner\_work\pyaedt\pyaedt\.venv\lib\site-packages\pyaedt\modules\MaterialLib.py:908: FutureWarning: Series.__getitem__ treating keys as positions is deprecated. In a future version, integer keys will always be treated as labels (consistent with DataFrame behavior). To access a value by position, use `ser.iloc[pos]`
and not (isinstance(val[keys.index(prop)], float) and math.isnan(val[keys.index(prop)]))
C:\actions-runner\_work\pyaedt\pyaedt\.venv\lib\site-packages\pyaedt\modules\MaterialLib.py:910: FutureWarning: Series.__getitem__ treating keys as positions is deprecated. In a future version, integer keys will always be treated as labels (consistent with DataFrame behavior). To access a value by position, use `ser.iloc[pos]`
props[prop] = float(val[keys.index(prop)])
Create design geometries#
Create rectangle and a circle and assign the material read from the .xlsx file. Create two new polylines and a region.
rect = mod2D.create_rectangle(origin=['r_x0', 'r_y0', 'r_z0'],
sizes=['r_dx', 'r_dy', 0],
name='Ground', matname=mats[0])
rect.color = (0, 0, 255) # rgb
rect.solve_inside = False
circle = mod2D.create_circle(position=['circle_x0', 'circle_y0', 'circle_z0'], radius='circle_radius',
num_sides='0', is_covered=True, name='Electrode', matname=mats[0])
circle.color = (0, 0, 255) # rgb
circle.solve_inside = False
poly1_points = [[-9, 2, 0], [-4, 2, 0], [2, -2, 0],[8, 2, 0]]
poly2_points = [[-9, 0, 0], [9, 0, 0]]
poly1_id = mod2D.create_polyline(points=poly1_points, segment_type='Spline', name='Poly1')
poly2_id = mod2D.create_polyline(points=poly2_points, name='Poly2')
mod2D.split([poly1_id, poly2_id], 'YZ', sides='NegativeOnly')
mod2D.create_region([20, 100, 20, 100])
<pyaedt.modeler.cad.object3d.Object3d object at 0x0000022600FCC970>
Define excitations#
Assign voltage excitations to rectangle and circle.
M2D.assign_voltage(rect.id, amplitude=0, name='Ground')
M2D.assign_voltage(circle.id, amplitude=50e6, name='50kV')
<pyaedt.modules.Boundary.BoundaryObject object at 0x0000022600F7B340>
Create initial mesh settings#
Assign a surface mesh to the rectangle.
M2D.mesh.assign_surface_mesh_manual(assignment=['Ground'], surface_deviation=0.001)
<pyaedt.modules.Mesh.MeshOperation object at 0x0000022600F7B1F0>
Create, validate and analyze the setup#
Create, update, validate and analyze the setup.
setup = M2D.create_setup(name=setup_name)
setup.props['PercentError'] = 0.5
setup.update()
M2D.validate_simple()
M2D.analyze_setup(setup_name)
True
Evaluate the E Field tangential component#
Evaluate the E Field tangential component along the given polylines. Add these operations to the Named Expression list in Field Calculator.
fields = M2D.ofieldsreporter
fields.CalcStack("clear")
fields.EnterQty("E")
fields.EnterEdge("Poly1")
fields.CalcOp("Tangent")
fields.CalcOp("Dot")
fields.AddNamedExpression("e_tan_poly1", "Fields")
fields.EnterQty("E")
fields.EnterEdge("Poly2")
fields.CalcOp("Tangent")
fields.CalcOp("Dot")
fields.AddNamedExpression("e_tan_poly2", "Fields")
Create Field Line Traces Plot#
Create Field Line Traces Plot specifying as seeding faces
the ground, the electrode and the region
and as In surface objects only the region.
plot = M2D.post.create_fieldplot_line_traces(seeding_faces=["Ground", "Electrode", "Region"],
in_volume_tracing_objs="Region", plot_name="LineTracesTest")
Update Field Line Traces Plot#
Update field line traces plot. Update seeding points number, line style and line width.
plot.SeedingPointsNumber = 20
plot.LineStyle = "Cylinder"
plot.LineWidth = 3
plot.update()
True
Export field line traces plot#
Export field line traces plot.
For field lint traces plot, the export file format is .fldplt.
M2D.post.export_field_plot(plot_name="LineTracesTest", output_dir=M2D.toolkit_directory, file_format="fldplt")
'D:/Temp/pyaedt_prj_0XX/Project_YJM.pyaedt\\LineTracesTest.fldplt'
Export the mesh field plot#
Export the mesh in aedtplt format.
M2D.post.export_mesh_obj(setup=M2D.nominal_adaptive)
'D:/Temp/pyaedt_prj_0XX/Project_YJM.pyaedt\\Design1\\Mesh_BWOK7T.aedtplt'
Save project and close AEDT#
Save the project and close AEDT.
M2D.save_project()
M2D.release_desktop()
True
Total running time of the script: (1 minutes 12.346 seconds)