import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 LPC map from NPSS"""
LPCMap = MapData()
# Map design point values
LPCMap.defaults = {}
LPCMap.defaults['alphaMap'] = 0.0
LPCMap.defaults['NcMap'] = 1.00
LPCMap.defaults['PR'] = 1.969
LPCMap.defaults['RlineMap'] = 2.150
LPCMap.RlineStall = 1.0
LPCMap.alphaMap = np.array([0.000, 90.000])
LPCMap.NcMap = np.array([
0.300, 0.400, 0.500, 0.600, 0.700, 0.750, 0.800, 0.850, 0.900, 0.950,
1.000, 1.050, 1.100, 1.150
])
LPCMap.RlineMap = np.array([
1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600, 2.800, 3.000
])
LPCMap.WcMap = np.array([[[
17.907, 19.339, 20.749, 22.136, 23.498, 24.833, 26.141, 27.420, 28.669,
29.887, 31.011
],
[
24.951, 26.742, 28.485, 30.177, 31.815, 33.397,
34.921, 36.385, 37.788, 39.128, 40.405
],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 LPT map from NPSS"""
LPTMap = MapData()
# Map design point values
LPTMap.defaults = {}
LPTMap.defaults['alphaMap'] = 1.0
LPTMap.defaults['NpMap'] = 100.0
LPTMap.defaults['PRmap'] = 6.0
#effMapDes = 0.9276 # = effMap for no scaling
LPTMap.alphaMap = np.array([1.0, 2.0])
LPTMap.NpMap = np.array([60., 70.0, 80.0, 90.0, 100., 110., 120.])
LPTMap.PRmap = np.array([
3.00, 3.25, 3.5, 3.75, 4.0, 4.25, 4.5, 4.75, 5.0, 5.25, 5.5, 5.75, 6.0,
6.25, 6.5, 6.75, 7.0, 7.25, 7.5, 8.0
])
LPTMap.effMap = np.array(
[[[
0.8560, 0.8532, 0.8497, 0.8463, 0.8432, 0.8402, 0.8371, 0.8341, 0.8313,
0.8288, 0.8263, 0.8240, 0.8214, 0.8188, 0.8165, 0.8142, 0.8121, 0.8101,
0.8082, 0.8046
],
[
0.8855, 0.8853, 0.8844, 0.8832, 0.8818, 0.8799, 0.8782, 0.8764,
0.8748, 0.8726, 0.8703, 0.8682, 0.8661, 0.8642, 0.8624, 0.8606,
0.8589, 0.8573, 0.8557, 0.8527
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 LPC map from NPSS"""
LPCmap = MapData()
# Map design point values
LPCmap.defaults = {}
LPCmap.defaults['alphaMap'] = 0.0
LPCmap.defaults['NcMap'] = 1.100
LPCmap.defaults['PR'] = 1.800
LPCmap.defaults['RlineMap'] = 2.200
LPCmap.RlineStall = 1.0
LPCmap.alphaMap = np.array([0.000, 1.000])
LPCmap.NcMap = np.array([0.300, 0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 1.000, 1.100, 1.200, 1.250])
LPCmap.RlineMap = np.array([1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600, 2.800, 3.000, 3.200])
LPCmap.WcMap = np.array([[[ 38.0744, 42.9399, 47.7510, 52.5016, 57.1863, 61.7994, 66.3359, 70.7905, 75.1584, 76.5663, 76.5663, 76.5663],
[ 54.0383, 60.0388, 65.9233, 71.6816, 77.3038, 82.7808, 88.1038, 93.2648, 98.2565, 101.0545, 101.0545, 101.0545],
[ 70.3200, 77.5153, 84.4949, 91.2421, 97.7419, 103.9805, 109.9459, 115.6273, 121.0156, 124.6409, 124.6409, 124.6409],
[ 87.4860, 95.6896, 103.5393, 111.0128, 118.0907, 124.7566, 130.9971, 136.8019, 142.1633, 146.2312, 146.2312, 146.2312],
[105.8588, 114.8071, 123.2285, 131.0978, 138.3948, 145.1045, 151.2169, 156.7268, 161.6340, 165.7319, 165.7319, 165.7319],
[125.1164, 134.6062, 143.3572, 151.3454, 158.5548, 164.9773, 170.6127, 175.4677, 179.5554, 182.8951, 183.0717, 183.0717],
[144.4910, 154.5703, 163.6243, 171.6346, 178.5959, 184.5149, 189.4099, 193.3090, 196.2491, 198.2745, 198.4155, 198.4155],
[165.9141, 176.2228, 185.1849, 192.7986, 199.0806, 204.0644, 207.7979, 210.3410, 211.7638, 212.1506, 212.1506, 212.1506],
[188.5677, 198.3532, 206.6834, 213.5745, 219.0613, 223.1942, 226.0370, 227.6647, 228.1611, 228.1611, 228.1611, 228.1611],
[214.1402, 222.1943, 228.9021, 234.2963, 238.4220, 241.3359, 243.1030, 243.7959, 243.8124, 243.8124, 243.8124, 243.8124],
[227.8569, 234.5820, 240.1193, 244.5040, 247.7802, 250.0000, 251.2213, 251.5216, 251.5216, 251.5216, 251.5216, 251.5216]],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of axi-5.map Compressor map from NPSS"""
AXI5 = MapData()
# Map design point values
AXI5.defaults = {}
AXI5.defaults['alphaMap'] = 0.0
AXI5.defaults['NcMap'] = 1.00
AXI5.defaults['PRmap'] = 5.2000
AXI5.defaults['RlineMap'] = 2.000
AXI5.RlineStall = 1.0
AXI5.alphaMap = np.array([0.000, 90.000])
AXI5.NcMap = np.array(
[0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100])
# , 4.8430, 5.1909, 5.5289, 5.8564, 6.1729, 6.4780, 6.7714, 7.0525, 7.3212])
AXI5.RlineMap = np.array(
[1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600])
AXI5.WcMap = np.array([[[
4.8430, 5.1909, 5.5289, 5.8564, 6.1729, 6.4780, 6.7714, 7.0525, 7.3212
], [
6.8115, 7.1360, 7.4477, 7.7462, 8.03130, 8.30260, 8.56000, 8.80330, 9.03230
], [
8.9765, 9.2855, 9.5780, 9.8535, 10.1121, 10.3536, 10.5781, 10.7855, 10.9760
],
[
11.0367, 11.4245, 11.7795, 12.1017, 12.3913,
12.6485, 12.8738, 13.0679, 13.2314
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 HPT map from NPSS"""
HPTmap = MapData()
# Map design point values
HPTmap.defaults = {}
HPTmap.defaults['alphaMap'] = 1.0
HPTmap.defaults['NpMap'] = 100.0
HPTmap.defaults['PRmap'] = 5.0
#effMapDes = 0.9276 # = effMap for no scaling
HPTmap.alphaMap = np.array([1.0, 2.0])
HPTmap.NpMap = np.array([60., 70.0, 80.0, 90.0, 100., 110.])
HPTmap.PRmap = np.array([
3.000, 3.250, 3.500, 3.750, 4.000, 4.250, 4.500, 4.750, 5.000, 5.250,
5.500, 5.750, 6.000, 6.250, 6.500, 6.750, 7.000, 7.250, 7.500, 8.000
])
HPTmap.effMap = np.array(
[[[
0.8460, 0.8405, 0.8349, 0.8296, 0.8249, 0.8206, 0.8165, 0.8127, 0.8092,
0.8060, 0.8030, 0.8002, 0.7976, 0.7953, 0.7931, 0.7911, 0.7892, 0.7875,
0.7858, 0.7826
],
[
0.8879, 0.8842, 0.8804, 0.8769, 0.8735, 0.8701, 0.8670, 0.8640,
0.8614, 0.8590, 0.8568, 0.8548, 0.8529, 0.8511, 0.8495, 0.8479,
0.8460, 0.8436, 0.8414, 0.8373
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of HBTF Fan map from NPSS"""
FanMap = MapData()
# Map design point values
FanMap.defaults = {}
FanMap.defaults['alphaMap'] = 0.0
FanMap.defaults['NcMap'] = 1.0
FanMap.defaults['PR'] = 1.39
FanMap.defaults['RlineMap'] = 2.00
FanMap.RlineStall = 1.0
FanMap.alphaMap = np.array([0.000, 5.000])
FanMap.NcMap = np.array([0.500, 0.600, 0.700, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100])
FanMap.RlineMap = np.array([1.000, 1.250, 1.500, 1.750, 2.000, 2.250, 2.500, 2.750, 3.000, 3.200, 3.250, 3.500])
FanMap.WcMap = np.array([[[1104.3380, 1247.2229, 1384.7280, 1514.7731, 1635.3656, 1742.0557, 1829.9553, 1866.5565, 1866.5565, 1866.5565, 1866.5565, 1866.5565],
[1369.0782, 1519.3645, 1661.7891, 1793.9161, 1913.4485, 2012.9872, 2084.3818, 2125.0225, 2126.4814, 2126.4814, 2126.4814, 2126.4814],
[1644.4929, 1791.9211, 1929.5121, 2054.8865, 2165.8337, 2252.1497, 2302.7083, 2315.6926, 2315.6926, 2315.6926, 2315.6926, 2315.6926],
[1919.7177, 2054.9785, 2179.5896, 2291.5830, 2389.1497, 2460.2048, 2491.7344, 2493.1619, 2493.1619, 2493.1619, 2493.1619, 2493.1619],
[2191.8401, 2304.7476, 2407.9941, 2500.2908, 2580.4556, 2636.5833, 2655.9902, 2655.9902, 2655.9902, 2655.9902, 2655.9902, 2655.9902],
[2324.3926, 2422.7034, 2512.5652, 2593.0525, 2663.3154, 2712.5415, 2729.2681, 2729.2681, 2729.2681, 2729.2681, 2729.2681, 2729.2681],
[2446.8806, 2530.6360, 2607.3428, 2676.3726, 2737.1450, 2780.9336, 2798.5239, 2798.8225, 2798.8225, 2798.8225, 2798.8225, 2798.8225],
[2564.0947, 2632.0720, 2694.6523, 2751.4587, 2802.1399, 2840.2256, 2858.9136, 2860.8069, 2860.8069, 2860.8069, 2860.8069, 2860.8069],
[2675.1184, 2726.3958, 2774.0383, 2817.8599, 2857.6853, 2889.2947, 2908.2825, 2914.2979, 2914.2979, 2914.2979, 2914.2979, 2914.2979]],
[[ 883.6700, 1019.1972, 1150.5194, 1275.3372, 1391.4280, 1494.5356, 1580.0111, 1639.0737, 1639.0737, 1639.0737, 1639.0737, 1639.0737],
[1103.0485, 1250.1055, 1390.6215, 1521.7380, 1640.7301, 1740.2568, 1812.2133, 1853.6100, 1860.6985, 1860.6985, 1860.6985, 1860.6985],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of HBTF Fan map from NPSS"""
FanMap = MapData()
# Map design point values
FanMap.defaults = {}
FanMap.defaults['alphaMap'] = 0.0
FanMap.defaults['NcMap'] = 0.99
FanMap.defaults['PR'] = 1.754
FanMap.defaults['RlineMap'] = 2.200
FanMap.RlineStall = 1.0
FanMap.alphaMap = np.array([0.000, 90.000])
FanMap.NcMap = np.array([0.300, 0.400, 0.500, 0.600, 0.700, 0.750, 0.800, 0.850, 0.900, 0.950, 1.000, 1.050, 1.100, 1.150])
FanMap.RlineMap = np.array([1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600, 2.800, 3.000])
FanMap.WcMap = np.array([[[121.797, 150.895, 179.422, 207.275, 234.359, 260.582, 285.861, 310.120, 333.291, 355.315, 369.552],
[194.867, 227.417, 258.872, 289.101, 317.981, 345.407, 371.284, 395.535, 418.097, 438.921, 457.973],
[265.640, 302.320, 337.109, 369.834, 400.351, 428.543, 454.320, 477.621, 498.412, 516.684, 532.453],
[330.650, 373.304, 412.750, 448.759, 481.163, 509.851, 534.771, 555.926, 573.366, 587.188, 597.526],
[380.373, 433.991, 481.896, 523.736, 559.314, 588.580, 611.621, 628.642, 639.948, 645.923, 647.172],
[408.374, 467.955, 520.262, 564.886, 601.649, 630.583, 651.910, 666.008, 673.384, 674.865, 674.865],
[432.494, 499.118, 556.491, 604.137, 641.921, 670.010, 688.838, 699.045, 701.536, 701.536, 701.536],
[470.522, 540.248, 599.144, 646.786, 683.165, 708.637, 723.856, 729.705, 729.841, 729.841, 729.841],
[527.146, 594.695, 650.583, 694.581, 726.881, 748.023, 758.831, 760.796, 760.796, 760.796, 760.796],
[593.025, 653.568, 702.700, 740.438, 767.140, 783.445, 790.213, 790.533, 790.533, 790.533, 790.533],
[643.809, 697.196, 739.568, 771.126, 792.350, 803.950, 806.892, 806.892, 806.892, 806.892, 806.892],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of HBTF HPC map from NPSS"""
HPCMap = MapData()
# Map design point values
HPCMap.defaults = {}
HPCMap.defaults['alphaMap'] = 0.0
HPCMap.defaults['NcMap'] = 0.976
HPCMap.defaults['PR'] = 10.894
HPCMap.defaults['RlineMap'] = 2.050
HPCMap.RlineStall = 1.0
HPCMap.alphaMap = np.array([0.000, 90.000])
HPCMap.NcMap = np.array([0.500, 0.600, 0.700, 0.750, 0.800, 0.850, 0.900, 0.925, 0.950, 0.975, 1.000, 1.025, 1.050, 1.150])
HPCMap.RlineMap = np.array([1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600, 2.800, 3.000])
HPCMap.WcMap = np.array([[[ 7.267, 7.543, 7.802, 8.046, 8.273, 8.484, 8.678, 8.855, 9.017, 9.162, 9.292],
[ 9.809, 10.081, 10.331, 10.558, 10.764, 10.948, 11.111, 11.253, 11.374, 11.475, 11.557],
[13.119, 13.369, 13.593, 13.791, 13.964, 14.111, 14.235, 14.335, 14.412, 14.467, 14.501],
[15.487, 15.754, 15.990, 16.196, 16.372, 16.519, 16.638, 16.729, 16.794, 16.834, 16.850],
[18.810, 19.118, 19.386, 19.614, 19.804, 19.956, 20.073, 20.155, 20.203, 20.219, 20.219],
[24.726, 25.163, 25.538, 25.851, 26.105, 26.302, 26.444, 26.532, 26.570, 26.572, 26.572],
[32.579, 33.147, 33.628, 34.025, 34.340, 34.576, 34.736, 34.824, 34.844, 34.844, 34.844],
[37.318, 37.919, 38.426, 38.843, 39.173, 39.418, 39.582, 39.669, 39.685, 39.685, 39.685],
[42.541, 43.049, 43.481, 43.840, 44.126, 44.343, 44.493, 44.579, 44.603, 44.603, 44.603],
[47.740, 48.152, 48.505, 48.801, 49.040, 49.225, 49.358, 49.440, 49.472, 49.473, 49.473],
[53.232, 53.466, 53.672, 53.849, 53.998, 54.120, 54.216, 54.285, 54.330, 54.349, 54.350],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 HPC map from NPSS"""
HPCmap = MapData()
# Map design point values
HPCmap.defaults = {}
HPCmap.defaults['alphaMap'] = 0.0
HPCmap.defaults['NcMap'] = 1.00
HPCmap.defaults['PR'] = 23.00
HPCmap.defaults['RlineMap'] = 2.0
HPCmap.RlineStall = 1.0
HPCmap.alphaMap = np.array([0.000, 90.000])
HPCmap.NcMap = np.array([
0.500, 0.600, 0.700, 0.750, 0.800, 0.850, 0.900, 0.925, 0.950, 0.975,
1.000, 1.025, 1.050
])
HPCmap.RlineMap = np.array([
1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600, 2.800, 3.000
])
HPCmap.WcMap = np.array(
[[[
22.7411, 24.0487, 25.1548, 26.0615, 26.7738, 27.2992, 27.6470, 27.8286,
27.8634, 27.8634, 27.8634
],
[
31.7548, 33.1181, 34.2670, 35.2054, 35.9397, 36.4783, 36.8308,
37.0085, 37.0362, 37.0362, 37.0362
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 LPT map from NPSS"""
LPTmap = MapData()
# Map design point values
LPTmap.defaults = {}
LPTmap.defaults['alphaMap'] = 1.0
LPTmap.defaults['NpMap'] = 100.0
LPTmap.defaults['PRmap'] = 7.5
#effMapDes = 0.9276 # = effMap for no scaling
LPTmap.alphaMap = np.array([1.0, 2.0])
LPTmap.NpMap = np.array([60.0, 70.0, 80.0, 90.0, 100.0, 110.0, 120.0])
LPTmap.PRmap = np.array([
3.000, 3.250, 3.500, 3.750, 4.000, 4.250, 4.500, 4.750, 5.000, 5.250,
5.500, 5.750, 6.000, 6.250, 6.500, 6.750, 7.000, 7.250, 7.500, 8.000
])
LPTmap.effMap = np.array(
[[[
0.8388, 0.8309, 0.8234, 0.8159, 0.8091, 0.8030, 0.7975, 0.7924, 0.7876,
0.7832, 0.7791, 0.7753, 0.7717, 0.7684, 0.7652, 0.7623, 0.7595, 0.7568,
0.7542, 0.7495
],
[
0.8878, 0.8813, 0.8745, 0.8685, 0.8629, 0.8577, 0.8528, 0.8484,
0.8443, 0.8404, 0.8368, 0.8334, 0.8302, 0.8272, 0.8242, 0.8210,
0.8179, 0.8150, 0.8122, 0.8071
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of CFM56 HPT map from NPSS"""
HPTMap = MapData()
# Map design point values
HPTMap.defaults = {}
HPTMap.defaults['alphaMap'] = 1.0
HPTMap.defaults['NpMap'] = 100.0
HPTMap.defaults['PRmap'] = 6.0
#effMapDes = 0.9276 # = effMap for no scaling
HPTMap.alphaMap = np.array([1.0, 2.0])
HPTMap.NpMap = np.array([60., 70.0, 80.0, 90.0, 100., 110.])
HPTMap.PRmap = np.array([
3.00, 3.25, 3.5, 3.75, 4.0, 4.25, 4.5, 4.75, 5.0, 5.25, 5.5, 5.75, 6.0,
6.25, 6.5, 6.75, 7.0, 7.25, 7.5, 8.0
])
HPTMap.effMap = np.array(
[[[
0.8449, 0.8394, 0.8358, 0.8330, 0.8306, 0.8285, 0.8242, 0.8201, 0.8163,
0.8128, 0.8094, 0.8063, 0.8026, 0.7992, 0.7959, 0.7928, 0.7899, 0.7872,
0.7846, 0.7798
],
[
0.8801, 0.8767, 0.8743, 0.8723, 0.8698, 0.8659, 0.8622, 0.8589,
0.8557, 0.8520, 0.8480, 0.8443, 0.8408, 0.8376, 0.8345, 0.8316,
0.8288, 0.8262, 0.8233, 0.8178
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of msfan3_3 fan map from NPSS"""
FanMap = MapData()
# Map design point values
FanMap.defaults = {}
FanMap.defaults['alphaMap'] = 0.0
FanMap.defaults['NcMap'] = 1.0
FanMap.defaults['PR'] = 3.3
FanMap.defaults['RlineMap'] = 2.00
RlineStall = 1.0
FanMap.alphaMap = np.array([0.000, 90.000])
FanMap.NcMap = np.array(
[0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100])
FanMap.RlineMap = np.array(
[1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600])
FanMap.WcMap = np.array([
# alpha = 0.
[[5.929, 6.562, 7.174, 7.763, 8.328, 8.866, 9.378, 9.862, 10.317],
[8.176, 8.853, 9.495, 10.101, 10.670, 11.201, 11.693, 12.144, 12.556],
[10.456, 11.200, 11.890, 12.523, 13.099, 13.617, 14.077, 14.479, 14.824],
[12.911, 13.755, 14.510, 15.176, 15.753, 16.241, 16.642, 16.959, 17.194],
[16.180, 17.144, 17.969, 18.653, 19.200, 19.615, 19.904, 20.074, 20.133],
[21.567, 22.428, 23.124, 23.661, 24.046, 24.286, 24.391, 24.398, 24.398],
[25.381, 26.000, 26.486, 26.844, 27.081, 27.203, 27.225, 27.225, 27.225],
[29.313, 29.579, 29.778, 29.912, 29.984, 30.000, 30.000, 30.000, 30.000],
[30.535, 30.639, 30.705, 30.735, 30.737, 30.737, 30.737, 30.737, 30.737],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of msfan3_3 fan map from NPSS"""
HPCmap = MapData()
# Map design point values
HPCmap.defaults = {}
HPCmap.defaults['alphaMap'] = 0.0
HPCmap.defaults['NcMap'] = 1.0
HPCmap.defaults['PR'] = 9.3
HPCmap.defaults['RlineMap'] = 2.00
HPCmap.defaults['effDes'] = 0.851
RlineStall = 1.0
HPCmap.alphaMap = np.array([0.000, 90.000])
HPCmap.NcMap = np.array(
[0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100])
HPCmap.RlineMap = np.array(
[1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600])
HPCmap.WcMap = np.array([
# alpha = 0.
[[3.278, 3.470, 3.656, 3.834, 4.004, 4.167, 4.321, 4.467, 4.604],
[4.457, 4.687, 4.903, 5.105, 5.293, 5.465, 5.623, 5.766, 5.895],
[6.137, 6.334, 6.515, 6.682, 6.833, 6.969, 7.090, 7.196, 7.287],
[8.080, 8.271, 8.443, 8.595, 8.728, 8.842, 8.938, 9.017, 9.078],
[11.324, 11.568, 11.780, 11.960, 12.111, 12.231, 12.323, 12.388, 12.427],
[18.702, 19.249, 19.705, 20.074, 20.357, 20.559, 20.684, 20.735, 20.737],
[23.752, 24.269, 24.699, 25.046, 25.311, 25.499, 25.614, 25.658, 25.659],
[29.110, 29.361, 29.574, 29.751, 29.893, 30.000, 30.074, 30.116, 30.128],
import numpy as np
from pycycle.maps.map_data import MapData
battery = MapData()
# Map design point values
battery.defaults = {}
battery.defaults['T_bp'] = 20.0
battery.defaults['SOC_bp'] = 0.5
battery.defaults['tU_oc'] = 3.75
battery.defaults['tC_Th'] = 20000.
battery.defaults['tR_Th'] = 0.002
battery.defaults['tR_0'] = 0.009
# Temperature break points (degC)
battery.T_bp = np.array([5., 20., 40])
battery.SOC_bp = np.array(
[0., 0.1, 0.25, 0.5, 0.75, 0.9, 1.]) # SOC break points
battery.tU_oc = np.array([[3.5, 3.55, 3.65, 3.75, 3.9, 4.1, 4.25],
[3.5, 3.55, 3.65, 3.75, 3.9, 4.1, 4.25],
[3.5, 3.55, 3.65, 3.75, 3.9, 4.1, 4.25]]) # volts, Open Circuit Voltage
# Thevenin components for Polarization Losses
battery.tC_Th = np.array([[2000., 2500., 22200., 1200., 20000., 14000., 10000.],
[1200., 20000., 40000., 20000., 30000., 20000., 25000.],
[30000., 35000., 49000., 27500., 50000., 29000., 30000.]]) # farads
battery.tR_Th = np.array([[0.0110, 0.0070, 0.0050, 0.0040, 0.0040, 0.0040, 0.0038],
[0.0030, 0.0029, 0.0029, 0.0020, 0.0029, 0.0025, 0.0024],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of axi-3_2.map Compressor map from NPSS"""
AXI3_2 = MapData()
# Map design point values
AXI3_2.defaults = {}
AXI3_2.defaults['alphaMap'] = 0.0
AXI3_2.defaults['NcMap'] = 1.00
AXI3_2.defaults['PRmap'] = 3.2000
AXI3_2.defaults['RlineMap'] = 2.000
AXI3_2.RlineStall = 1.0
AXI3_2.alphaMap = np.array([0.000, 90.000])
AXI3_2.NcMap = np.array(
[0.400, 0.500, 0.600, 0.700, 0.800, 0.900, 0.950, 1.000, 1.050, 1.100])
AXI3_2.RlineMap = np.array(
[1.000, 1.200, 1.400, 1.600, 1.800, 2.000, 2.200, 2.400, 2.600])
AXI3_2.WcMap = np.array([[[6.656, 7.378, 8.080, 8.763, 9.423, 10.060, 10.673,
11.260, 11.820],
[9.491, 10.188, 10.860, 11.505, 12.122,
12.710, 13.268, 13.796, 14.293],
[11.897, 12.714, 13.487, 14.214, 14.893,
15.524, 16.105, 16.637, 17.120],
[14.089, 15.192, 16.204, 17.124, 17.948,
18.676, 19.309, 19.847, 20.292],
import numpy as np
from pycycle.maps.map_data import MapData
"""Python version of npc01.map from NPSS"""
NCP01 = MapData()
# Map design point values
NCP01.defaults = {}
NCP01.defaults['alphaMap'] = 0.0 # Unscaled map design point vane angle
NCP01.defaults[
'NcMap'] = 1.0 # Unscaled map design point corrected speed. Max eff point at 100% speed
NCP01.defaults[
'RlineMap'] = 2.0 # Unscaled map design point Rline. Max eff point at 100% speed
NCP01.RlineStall = 1.0 # Rline at stall
# Values along map axes
NCP01.alphaMap = np.array([0.0, 90.0])
NCP01.NcMap = np.array(
[0.5, 0.6, 0.7, 0.8, 0.85, 0.9, 0.95, 1.0, 1.05, 1.1, 1.15])
NCP01.RlineMap = np.array(
[1.0, 1.2, 1.4, 1.6, 1.8, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2])
# Create map arrays
#Nc_data, alpha_data, Rline_data = np.meshgrid(Nc_vals, alpha_vals, Rline_vals, sparse=False)
NCP01.WcMap = np.array(
[[[
1129.0825, 1258.8474, 1382.1909, 1498.6384, 1607.7749, 1709.2917,
1802.8665, 1888.3794, 1965.6499, 2034.6738, 2095.4690, 2148.1226
],
[