def test_Lam_Wind_28_wind_rad_tan(self):
"""Test machine plot with Slot 28 and winding rad=1, tan=2 and rad=2 and tan=1
"""
print("\nTest plot Slot 28")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=35e-3,
Rext=84e-3,
is_internal=True,
is_stator=False,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.rotor.axial_vent = [
VentilationCirc(Zh=6, Alpha0=pi / 6, D0=15e-3, H0=0.045)
]
test_obj.rotor.slot = SlotW28(
Zs=42, W0=3.5e-3, H0=0.45e-3, R1=3.5e-3, H3=14e-3, W3=5e-3
)
test_obj.rotor.winding = WindingCW2LT(qs=3, p=3)
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=85e-3,
Rext=0.2,
is_internal=False,
is_stator=True,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW28(
Zs=18, W0=7e-3, R1=10e-3, H0=5e-3, H3=30e-3, W3=5e-3
)
test_obj.stator.winding = WindingDW2L(qs=3, p=3, Lewout=60e-3)
test_obj.stator.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.frame = Frame(Rint=0.2, Rext=0.25, Lfra=1)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_1-Machine.png"))
# Rotor + stator + 2 for frame + 1 for Shaft
self.assertEqual(len(fig.axes[0].patches), 133)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_2-Rotor.png"))
# 2 for lam + Zs*2 for wind + 6 vent
self.assertEqual(len(fig.axes[0].patches), 92)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_3-Stator.png"))
# 2 for lam, 2*Zs for wind
self.assertEqual(len(fig.axes[0].patches), 38)
def test_Lam_Wind_22_wind_22(self):
"""Test machine plot with Slot 22 and winding rad=2, tan=2
"""
print("\nTest plot Slot 22")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.8,
Nrvd=4,
Wrvd=0.05,
)
test_obj.rotor.slot = SlotW22(Zs=6,
W0=pi / 20,
W2=pi / 10,
H0=20e-3,
H2=150e-3)
test_obj.rotor.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=0.8,
is_internal=False,
is_stator=True,
L1=0.8,
Nrvd=4,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW22(Zs=18,
W0=pi / 20,
W2=pi / 10,
H0=20e-3,
H2=150e-3)
test_obj.stator.winding = WindingDW2L(qs=3, p=3)
test_obj.stator.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.stator.winding.Lewout = 60e-3
test_obj.frame = Frame(Rint=0.8, Rext=0.8, Lfra=1)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s22_1-Machine.png"))
# Rotor + Stator + 0 for frame + 0 for shaft
self.assertEqual(len(fig.axes[0].patches), 63)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s22_2-Rotor.png"))
# 1 for lam + 4*Zs for wind
self.assertEqual(len(fig.axes[0].patches), 25)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s22_3-Stator.png"))
# 2 for lam + Zs*2 for wind
self.assertEqual(len(fig.axes[0].patches), 38)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s22_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s22_Tooth_out.png"))
def test_Lam_Wind_27_wind_22(self):
"""Test machine plot with Slot 27 and winding rad=2, tan=2
"""
print("\nTest plot Slot 27")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.2,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.rotor.axial_vent = [
VentilationCirc(Zh=6, Alpha0=pi / 6, D0=60e-3, H0=0.35)
]
test_obj.rotor.slot = SlotW27(Zs=6,
H0=0.05,
W0=30e-3,
H1=0.125,
W1=0.06,
H2=0.05,
W2=0.09,
W3=0.04)
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.rotor.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=0.8,
is_internal=False,
is_stator=True,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW27(Zs=18,
H0=0.05,
W0=30e-3,
H1=0.07,
W1=0.05,
H2=0.15,
W2=0.1,
W3=0.05)
test_obj.stator.winding = WindingDW2L(qs=3, p=3)
test_obj.stator.winding.Lewout = 60e-3
test_obj.stator.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.frame = Frame(Rint=0.8, Rext=0.9, Lfra=1)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s27_1-Machine.png"))
# Rotor + Stator + 2 for frame + 1 for shaft
self.assertEqual(len(fig.axes[0].patches), 73)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s27_2-Rotor.png"))
# 2 for lam + 4*Zs for wind + 6 vent
self.assertEqual(len(fig.axes[0].patches), 32)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s27_3-Stator.png"))
# 2 for lam + 2*Zs for wind
self.assertEqual(len(fig.axes[0].patches), 38)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s27_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s27_Tooth_out.png"))
def test_Lam_Wind_14_wind_22(self):
"""Test machine plot with Slot 14 and winding rad=2, tan=2
"""
print("\nTest plot Slot 14")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.2,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.95,
Nrvd=1,
Wrvd=0.05,
)
test_obj.rotor.slot = SlotW14(Zs=6,
W0=100e-3,
W3=200e-3,
H0=15e-3,
H1=25e-3,
H3=140e-3)
test_obj.rotor.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.shaft.mat_type.name = "M270_35A"
test_obj.rotor.mat_type.name = "Load_M400"
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=1,
is_internal=False,
is_stator=True,
L1=0.95,
Nrvd=1,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW14(Zs=6,
W0=100e-3,
W3=200e-3,
H0=15e-3,
H1=25e-3,
H3=140e-3)
test_obj.stator.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.frame = Frame(Rint=1, Rext=1.1, Lfra=1)
test_obj.frame.mat_type.name = "M330_35A"
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s14_1-Machine.png"))
# Rotor + Stator + 2 for frame + 1 for shaft
self.assertEqual(len(fig.axes[0].patches), 55)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s14_2-Rotor.png"))
# 2 for lam + Zs*4 for wind
self.assertEqual(len(fig.axes[0].patches), 26)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s14_3-Stator.png"))
# 2 for lam + Zs*4 for wind
self.assertEqual(len(fig.axes[0].patches), 26)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s14_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s14_Tooth_out.png"))
def test_Lam_Wind_28_wind_22(self):
"""Test machine plot with Slot 28 and winding rad=2, tan=2
"""
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.01,
Rext=0.129,
is_internal=True,
is_stator=False,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.rotor.slot = SlotW28(
Zs=6, W0=20e-3, R1=25e-3, H0=10e-3, H3=50e-3, W3=15e-3
)
test_obj.rotor.winding = WindingUD(
user_wind_mat=wind_mat, qs=4, p=4, Lewout=60e-3
)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.13,
Rext=0.4,
is_internal=False,
is_stator=True,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.stator.axial_vent = [
VentilationCirc(Zh=6, Alpha0=pi / 6, D0=60e-3, H0=0.25)
]
test_obj.stator.slot = SlotW28(
Zs=6, W0=40e-3, R1=50e-3, H0=10e-3, H3=70e-3, W3=85e-3
)
test_obj.stator.winding = WindingUD(
user_wind_mat=wind_mat, qs=4, p=4, Lewout=60e-3
)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_4-Machine.png"))
# Rotor + stator + 0 for frame + 1 for shaft
self.assertEqual(len(fig.axes[0].patches), 59)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_5-Rotor.png"))
# 2 for lam + Zs*4 for wind
self.assertEqual(len(fig.axes[0].patches), 26)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_6-Stator.png"))
# 2 for lam, 4*Zs for wind + 6 vents
self.assertEqual(len(fig.axes[0].patches), 32)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s28_Tooth_out.png"))
def test_Lam_Wind_12_wind_22(self):
"""Test machine plot with Slot 12 and winding rad=2, tan=2
"""
print("\nTest plot Slot 12")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.2,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.rotor.axial_vent = [
VentilationPolar(Zh=6, Alpha0=pi / 6, W1=pi / 6, D0=100e-3, H0=0.3)
]
test_obj.rotor.slot = SlotW12(Zs=6,
R2=35e-3,
H0=20e-3,
R1=17e-3,
H1=130e-3)
test_obj.rotor.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=0.8,
is_internal=False,
is_stator=True,
L1=0.9,
Nrvd=2,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW12(Zs=18,
R2=25e-3,
H0=30e-3,
R1=0,
H1=150e-3)
test_obj.stator.winding.Lewout = 60e-3
test_obj.stator.winding = WindingDW2L(qs=3, p=3)
test_obj.stator.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.frame = Frame(Rint=0.8, Rext=0.9, Lfra=1)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s12_1-Machine.png"))
# Rotor + Stator + 2 for frame + 1 for Shaft
self.assertEqual(len(fig.axes[0].patches), 73)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s12_2-Rotor.png"))
# 2 for lam + Zs*4 for wind + 6 vents
self.assertEqual(len(fig.axes[0].patches), 32)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s12_3-Stator.png"))
# 2 for lam + Zs*2 for wind
self.assertEqual(len(fig.axes[0].patches), 38)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s12_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s12_Tooth_out.png"))
def test_Lam_Wind_10_wind_22(self):
"""Test machine plot with Slot 10 and winding rad=2, tan=2
"""
print("\nTest plot Slot 10")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.2,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.95,
Nrvd=1,
Wrvd=0.05,
)
test_obj.rotor.slot = SlotW10(
Zs=6,
W0=50e-3,
W1=90e-3,
W2=100e-3,
H0=20e-3,
H1=35e-3,
H2=130e-3,
H1_is_rad=False,
)
test_obj.rotor.winding = WindingUD(
user_wind_mat=wind_mat, qs=4, p=4, Lewout=60e-3
)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=0.8,
is_internal=False,
is_stator=True,
L1=0.95,
Nrvd=1,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW10(
Zs=6,
W0=50e-3,
W1=80e-3,
W2=50e-3,
H0=15e-3,
H1=25e-3,
H2=140e-3,
H1_is_rad=False,
)
test_obj.stator.winding = WindingUD(
user_wind_mat=wind_mat, qs=4, p=4, Lewout=60e-3
)
test_obj.frame = Frame(Rint=0.8, Rext=0.9, Lfra=1)
test_obj.frame.mat_type.name = "M330_35A"
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_1-Machine.png"))
# Rotor + Stator + 2 for frame + 1 for Shaft
self.assertEqual(len(fig.axes[0].patches), 55)
test_obj.rotor.plot()
fig = plt.gcf()
self.assertEqual(len(fig.axes[0].patches), 26)
fig.savefig(join(save_path, "test_Lam_Wind_s10_2-Rotor.png"))
# 2 for lam + Zs*4 for wind
self.assertEqual(len(fig.axes[0].patches), 26)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_3-Stator.png"))
# 2 for lam + Zs*4 for wind
self.assertEqual(len(fig.axes[0].patches), 26)
lines = test_obj.stator.slot.build_geometry_half_tooth(is_top=False)
surf = SurfLine(line_list=lines)
surf.plot_lines()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_bottom_out.png"))
lines = test_obj.stator.slot.build_geometry_half_tooth(is_top=True)
surf = SurfLine(line_list=lines)
surf.plot_lines()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_top_out.png"))
lines = test_obj.rotor.slot.build_geometry_half_tooth(is_top=False)
surf = SurfLine(line_list=lines)
surf.plot_lines()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_bottom_in.png"))
lines = test_obj.rotor.slot.build_geometry_half_tooth(is_top=True)
surf = SurfLine(line_list=lines)
surf.plot_lines()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_top_in.png"))
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
mesh_dict = tooth.comp_mesh_dict(5e-3)
for line in tooth.get_lines():
mid = line.get_middle()
plt.text(mid.real, mid.imag, str(mesh_dict[line.label]))
fig.savefig(join(save_path, "test_Lam_Wind_s10_Tooth_out.png"))
def test_Lam_Wind_21_wind_22(self):
"""Test machine plot with Slot 21 and winding rad=2, tan=2
"""
print("\nTest plot Slot 21")
plt.close("all")
test_obj = MachineDFIM()
test_obj.rotor = LamSlotWind(
Rint=0.2,
Rext=0.5,
is_internal=True,
is_stator=False,
L1=0.85,
Nrvd=3,
Wrvd=0.05,
)
test_obj.rotor.slot = SlotW21(
Zs=6,
W0=40e-3,
W1=60e-3,
W2=40e-3,
H0=20e-3,
H1=0,
H2=130e-3,
H1_is_rad=False,
)
test_obj.rotor.axial_vent.append(
VentilationTrap(Zh=6,
Alpha0=pi / 6,
W1=30e-3,
W2=60e-3,
D0=0.05,
H0=0.3))
test_obj.rotor.axial_vent.append(
VentilationTrap(Zh=6,
Alpha0=pi / 6,
W1=60e-3,
W2=90e-3,
D0=0.05,
H0=0.4))
test_obj.rotor.winding = WindingUD(user_wind_mat=wind_mat,
qs=4,
p=4,
Lewout=60e-3)
test_obj.rotor.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.shaft = Shaft(Drsh=test_obj.rotor.Rint * 2, Lshaft=1)
test_obj.stator = LamSlotWind(
Rint=0.51,
Rext=0.8,
is_internal=False,
is_stator=True,
L1=0.85,
Nrvd=3,
Wrvd=0.05,
)
test_obj.stator.slot = SlotW21(
Zs=18,
W0=40e-3,
W1=60e-3,
W2=90e-3,
H0=15e-3,
H1=35e-3,
H2=140e-3,
H1_is_rad=False,
)
test_obj.stator.winding = WindingDW2L(qs=3, p=3)
test_obj.stator.axial_vent.append(
VentilationCirc(Zh=12, Alpha0=pi / 6, D0=50e-3, H0=0.75))
test_obj.stator.mat_type.mag = MatMagnetics(Wlam=0.5e-3)
test_obj.stator.winding.Lewout = 60e-3
test_obj.frame = Frame(Rint=0.8, Rext=1, Lfra=1)
test_obj.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s21_1-Machine.png"))
# Rotor + Stator + 2 for frame + 1 for shaft
self.assertEqual(len(fig.axes[0].patches), 91)
test_obj.rotor.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s21_2-Rotor.png"))
# 2 for lam + 4*Zs for wind + 6 vents + 6 vents
self.assertEqual(len(fig.axes[0].patches), 38)
test_obj.stator.plot()
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s21_3-Stator.png"))
# 2 for lam + Zs*2 for wind + 12 vents
self.assertEqual(len(fig.axes[0].patches), 50)
tooth = test_obj.rotor.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s21_Tooth_in.png"))
tooth = test_obj.stator.slot.get_surface_tooth()
tooth.plot(color="r")
fig = plt.gcf()
fig.savefig(join(save_path, "test_Lam_Wind_s21_Tooth_out.png"))