Source code for Tests.Methods.Slot.test_SlotW25_meth

# -*- coding: utf-8 -*-
import pytest

from pyleecan.Classes.SlotW25 import SlotW25
from numpy import ndarray, angle
from pyleecan.Classes.LamSlot import LamSlot
from pyleecan.Classes.Slot import Slot
from pyleecan.Methods.Slot.SlotW25 import S25_HWCheckError

# For AlmostEqual
DELTA = 1e-4

slotW25_test = list()

# Internal Slot
lam = LamSlot(is_internal=True, Rext=0.1325)
lam.slot = SlotW25(Zs=10, H1=3e-3, H2=30e-3, W3=20e-3, W4=40e-3)
slotW25_test.append(
    {
        "test_obj": lam,
        "S_exp": 0.001678213,
        "Ao": 0.32527,
        "Aw": 0.45157,
        "SW_exp": 0.0015532047,
        "H_exp": 0.032848,
    }
)

# External Slot
lam = LamSlot(is_internal=False, Rint=0.1325)
lam.slot = SlotW25(Zs=10, H1=3e-3, H2=30e-3, W3=20e-3, W4=40e-3)
slotW25_test.append(
    {
        "test_obj": lam,
        "S_exp": 0.00236077388,
        "Ao": 0.32527,
        "Aw": 0.49427,
        "SW_exp": 0.00223015199,
        "H_exp": 0.032899,
    }
)


[docs]class Test_SlotW25_meth(object): """pytest for SlotW25 methods"""
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_schematics(self, test_dict): """Check that the schematics is correct""" test_obj = test_dict["test_obj"] point_dict = test_obj.slot._comp_point_coordinate() assert angle(point_dict["Z1"]) < 0 assert angle(point_dict["Z2"]) < 0 assert angle(point_dict["Z3"]) < 0 assert angle(point_dict["Z4"]) < 0 assert angle(point_dict["Z5"]) > 0 assert angle(point_dict["Z6"]) > 0 assert angle(point_dict["Z7"]) > 0 assert angle(point_dict["Z8"]) > 0 # Check height assert abs(point_dict["Z1"] - point_dict["Z2"]) == pytest.approx( test_obj.slot.H1 ) assert abs(point_dict["Z3"] - point_dict["Z4"]) == pytest.approx( test_obj.slot.H2 ) assert abs(point_dict["Z8"] - point_dict["Z7"]) == pytest.approx( test_obj.slot.H1 ) assert abs(point_dict["Z6"] - point_dict["Z5"]) == pytest.approx( test_obj.slot.H2 ) # Radius assert abs(point_dict["Z2"]) == pytest.approx(abs(point_dict["Z3"])) assert abs(point_dict["Z3"]) == pytest.approx(abs(point_dict["Z6"])) assert abs(point_dict["Z6"]) == pytest.approx(abs(point_dict["Z7"]))
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_build_geometry_active(self, test_dict): """Check that the active geometry is correctly split""" test_obj = test_dict["test_obj"] surf_list = test_obj.slot.build_geometry_active(Nrad=3, Ntan=2) # Check label assert surf_list[0].label == "Wind_Stator_R0_T0_S0" assert surf_list[1].label == "Wind_Stator_R1_T0_S0" assert surf_list[2].label == "Wind_Stator_R2_T0_S0" assert surf_list[3].label == "Wind_Stator_R0_T1_S0" assert surf_list[4].label == "Wind_Stator_R1_T1_S0" assert surf_list[5].label == "Wind_Stator_R2_T1_S0" # Check tangential position assert surf_list[0].point_ref.imag < 0 assert surf_list[1].point_ref.imag < 0 assert surf_list[2].point_ref.imag < 0 assert surf_list[3].point_ref.imag > 0 assert surf_list[4].point_ref.imag > 0 assert surf_list[5].point_ref.imag > 0 # Check radial position if test_obj.is_internal: # Tan=0 assert surf_list[0].point_ref.real > surf_list[1].point_ref.real assert surf_list[1].point_ref.real > surf_list[2].point_ref.real # Tan=1 assert surf_list[3].point_ref.real > surf_list[4].point_ref.real assert surf_list[4].point_ref.real > surf_list[5].point_ref.real else: # Tan=0 assert surf_list[0].point_ref.real < surf_list[1].point_ref.real assert surf_list[1].point_ref.real < surf_list[2].point_ref.real # Tan=1 assert surf_list[3].point_ref.real < surf_list[4].point_ref.real assert surf_list[4].point_ref.real < surf_list[5].point_ref.real
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_comp_surface(self, test_dict): """Check that the computation of the surface is correct""" test_obj = test_dict["test_obj"] result = test_obj.slot.comp_surface() a = result b = test_dict["S_exp"] msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg # Check that the analytical method returns the same result as the numerical one b = Slot.comp_surface(test_obj.slot) msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_comp_surface_active(self, test_dict): """Check that the computation of the winding surface is correct""" test_obj = test_dict["test_obj"] result = test_obj.slot.comp_surface_active() a = result b = test_dict["SW_exp"] msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg # Check that the analytical method returns the same result as the numerical one b = Slot.comp_surface_active(test_obj.slot) msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_comp_height(self, test_dict): """Check that the computation of the height is correct""" test_obj = test_dict["test_obj"] result = test_obj.slot.comp_height() a = result b = test_dict["H_exp"] msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg # Check that the analytical method returns the same result as the numerical one b = Slot.comp_height(test_obj.slot) msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_comp_angle_opening(self, test_dict): """Check that the computation of the average opening angle iscorrect""" test_obj = test_dict["test_obj"] a = test_obj.slot.comp_angle_opening() b = test_dict["Ao"] msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg # Check that the analytical method returns the same result as the numerical one b = Slot.comp_angle_opening(test_obj.slot) msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg
[docs] @pytest.mark.parametrize("test_dict", slotW25_test) def test_comp_angle_active_eq(self, test_dict): """Check that the computation of the average angle is correct""" test_obj = test_dict["test_obj"] result = test_obj.slot.comp_angle_active_eq() a = result b = test_dict["Aw"] msg = "Return " + str(a) + " expected " + str(b) assert abs((a - b) / a - 0) < DELTA, msg
[docs] def test_check_error(self): """Check that the check method is correctly raising an error""" lam = LamSlot(is_internal=True, Rext=0.1325) lam.slot = SlotW25(Zs=100, H1=3e-3, H2=30e-3, W3=20e-3, W4=0.0000000000015) with pytest.raises(S25_HWCheckError) as context: lam.slot.check()
[docs] def test_get_surface_active(self): """Check that the get_surface_active works when stator = false""" lam = LamSlot(is_internal=True, Rext=0.1325, is_stator=False) lam.slot = SlotW25(Zs=10, H1=3e-3, H2=30e-3, W3=20e-3, W4=40e-3) result = lam.slot.get_surface_active() assert result.label == "Wind_Rotor_R0_T0_S0" assert len(result.get_lines()) == 6