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 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 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 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 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 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 LPT map from NPSS"""
HPT1269 = MapData()
# Map design point values
HPT1269.defaults = {}
HPT1269.defaults['alphaMap'] = 1.0
HPT1269.defaults['NpMap'] = 100.0
HPT1269.defaults['PRmap'] = 6.0
#effMapDes = 0.9276 # = effMap for no scaling
HPT1269.alphaMap = np.array([1.0, 2.0])
HPT1269.NpMap = np.array([60., 70.0, 80.0, 90.0, 100., 110.])
HPT1269.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
])
HPT1269.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 lpt2269.map Turbine map from NPSS"""
LPT2269 = MapData()
# Map design point values
LPT2269.defaults = {}
LPT2269.defaults['alphaMap'] = 1.0
LPT2269.defaults['NpMap'] = 100.0
LPT2269.defaults['PRmap'] = 6.0 # Design Pressure Ratio
#effMapDes = 0.9276 # = effMap for no scaling
LPT2269.alphaMap = np.array([1.0, 2.0])
LPT2269.NpMap = np.array([60., 70.0, 80.0, 90.0, 100., 110., 120.])
LPT2269.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
])
LPT2269.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
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
],
[