-
Notifications
You must be signed in to change notification settings - Fork 10
/
basic_transformer_center_tapped.py
152 lines (132 loc) · 7.5 KB
/
basic_transformer_center_tapped.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
"""
Basic example to show how to simulate a center-tapped transformer.
After starting the program, the geometry dimensions are displayed. Verify this geometry, close the window, to continue the simulation.
After a short time, B-Field and winding losses simulation results are shown. Winding losses are shown as a colormap.
In the core, the magnitude B-Field in Tesla is shown. With the gmsh window, one can move the picture in the 3D way (not recommended).
If you close this window, the thermal simulation will be continued, if programmed. If true, the thermal heat distribution will be displayed.
To continue with the next simulation (or end the program), you need to close this window. All results are written to the result
folder .../femmt/examples/example_results/simulation_file_name/results/log_electro_magnetic.json. and .../results_thermal.json.
"""
import femmt as fmt
import os
def basic_example_transformer_center_tapped(onelab_folder: str = None, show_visual_outputs: bool = True,
is_test: bool = False):
"""
Run the example code for the center-tapped transformer.
:param onelab_folder: onelab folder path
:type onelab_folder: str
:param show_visual_outputs: True to show visual outputs (simulation results)
:type show_visual_outputs: bool
:param is_test: True for pytest usage. Defaults to False.
:type is_test: bool
"""
def example_thermal_simulation(show_thermal_visual_outputs: bool = True, flag_insulation: bool = True):
# Thermal simulation:
# The losses calculated by the magnetics simulation can be used to calculate the heat distribution of the
# given magnetic component. In order to use the thermal simulation, thermal conductivities for each material
# can be entered as well as a boundary temperature which will be applied on the boundary of the
# simulation (dirichlet boundary condition).
# The case parameter sets the thermal conductivity for a case which will be set around the core.
# This could model some case in which the transformer is placed in together with a set potting material.
thermal_conductivity_dict = {
"air": 0.0263,
"case": { # epoxy resign
"top": 1.54,
"top_right": 1.54,
"right": 1.54,
"bot_right": 1.54,
"bot": 1.54
},
"core": 5, # ferrite
"winding": 400, # copper
"air_gaps": 180, # aluminium nitride
"insulation": 0.42 if flag_insulation else None # polyethylene
}
# Here the case size can be determined
case_gap_top = 0.002
case_gap_right = 0.0025
case_gap_bot = 0.002
# Here the boundary temperatures can be set, currently it is set to 20°C (around 293°K).
# This does not change the results of the simulation (at least when every boundary is set equally)
# but will set the temperature offset.
boundary_temperatures = {
"value_boundary_top": 20,
"value_boundary_top_right": 20,
"value_boundary_right_top": 20,
"value_boundary_right": 20,
"value_boundary_right_bottom": 20,
"value_boundary_bottom_right": 20,
"value_boundary_bottom": 20
}
# Here the boundary sides can be turned on (1) or off (0)
# By turning off the flag a neumann boundary will be applied at this point with heat flux = 0
boundary_flags = {
"flag_boundary_top": 0,
"flag_boundary_top_right": 0,
"flag_boundary_right_top": 1,
"flag_boundary_right": 1,
"flag_boundary_right_bottom": 1,
"flag_boundary_bottom_right": 1,
"flag_boundary_bottom": 1
}
# In order for the thermal simulation to work an electro_magnetic simulation has to run before.
# The em-simulation will create a file containing the losses.
# When the losses file is already created and contains the losses for the current model, it is enough to
# run geo.create_model in order for the thermal simulation to work (geo.single_simulation is not needed).
# Obviously when the model is modified and the losses can be out of date and therefore the
# geo.single_simulation needs to run again.
geo.thermal_simulation(thermal_conductivity_dict, boundary_temperatures, boundary_flags, case_gap_top,
case_gap_right, case_gap_bot, show_thermal_visual_outputs,
color_scheme=fmt.colors_ba_jonas, colors_geometry=fmt.colors_geometry_ba_jonas,
flag_insulation=flag_insulation)
example_results_folder = os.path.join(os.path.dirname(__file__), "example_results")
if not os.path.exists(example_results_folder):
os.mkdir(example_results_folder)
working_directory = os.path.join(example_results_folder, os.path.splitext(os.path.basename(__file__))[0])
if not os.path.exists(working_directory):
os.mkdir(working_directory)
geo = fmt.MagneticComponent(component_type=fmt.ComponentType.Transformer, working_directory=working_directory,
verbosity=fmt.Verbosity.ToConsole, is_gui=is_test)
# This line is for automated pytest running on GitHub only. Please ignore this line!
if onelab_folder is not None:
geo.file_data.onelab_folder_path = onelab_folder
core_dimensions = fmt.dtos.SingleCoreDimensions(window_h=0.025, window_w=0.02, core_inner_diameter=0.015,
core_h=0.05)
core = fmt.Core(core_dimensions=core_dimensions, mu_r_abs=3100, phi_mu_deg=12, sigma=1.2,
permeability_datasource=fmt.MaterialDataSource.Custom,
permittivity_datasource=fmt.MaterialDataSource.Custom)
geo.set_core(core)
air_gaps = fmt.AirGaps(fmt.AirGapMethod.Percent, core)
air_gaps.add_air_gap(fmt.AirGapLegPosition.CenterLeg, 0.0005, 50)
geo.set_air_gaps(air_gaps)
# set_center_tapped_windings() automatically places the conductors
insulation, winding_window = fmt.functions_topologies.set_center_tapped_windings(
core=core,
primary_turns=12,
primary_radius=1.1e-3,
primary_number_strands=50,
primary_strand_radius=0.00011,
secondary_parallel_turns=3,
secondary_thickness_foil=1e-3,
iso_top_core=0.001,
iso_bot_core=0.001,
iso_left_core=0.002,
iso_right_core=0.001,
iso_primary_to_primary=1e-4,
iso_secondary_to_secondary=2e-4,
iso_primary_to_secondary=5e-4,
interleaving_type=fmt.CenterTappedInterleavingType.TypeA,
interleaving_scheme=fmt.InterleavingSchemesFoilLitz.ter_3_4_sec_ter_4_3_sec,
primary_additional_bobbin=100,
winding_temperature=100,
center_foil_additional_bobbin=0e-3,
wrap_para_type=fmt.WrapParaType.FixedThickness,
foil_horizontal_placing_strategy=fmt.FoilHorizontalDistribution.VerticalUpward)
geo.set_insulation(insulation)
geo.set_winding_windows([winding_window])
geo.create_model(freq=200000, pre_visualize_geometry=show_visual_outputs)
geo.single_simulation(freq=200000, current=[20, 120, 120], phi_deg=[0, 180, 180],
show_fem_simulation_results=show_visual_outputs)
example_thermal_simulation(show_visual_outputs, flag_insulation=False)
if __name__ == "__main__":
basic_example_transformer_center_tapped(show_visual_outputs=True)