def test_TotalLoads(inputs_str, inputs_aeroloads):
npts = np.size(inputs_str['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', TotalLoads(npts), promotes=['*'])
prob = init_IndepVar_add(prob, inputs_str)
prob = init_IndepVar_add(prob, inputs_aeroloads)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, inputs_str)
prob = init_IndepVar_set(prob, inputs_aeroloads)
check_gradient_unit_test(prob, tol=.01, display=False)
def test_TestRootMoment():
data_aero = inputs_aeroloads_operating()
data_str = inputs_str()
data = {}
data['aeroloads_r'] = np.array([1.5375, 1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697, 2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267, 8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998, 18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509, 28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998, 42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998, 56.20419994, 58.9374998, 61.67080026, 63.03135, 63.0375])
data['aeroloads_Px'] = np.array([0., 610.90043568, 469.20341791, 529.97827499, 580.02866009, 625.05962708, 661.78185992, 694.27930979, 830.93330879, 856.02758285, 872.85385017, 1155.30591175, 1278.65677079, 1394.84126103, 1598.03834945, 1726.35834871, 1933.87963881, 2271.78104859, 2748.8320438, 3083.27106663, 3677.62201087, 4343.75465079, 5051.54596049, 5731.9904625, 7004.11264534, 7364.31723123, 8093.26478732, 8599.46800605, 9000.96202336, 9390.16142378, 9555.38019573, 9708.33654591, 9964.10107999, 10160.41657631, 10284.8330047, 10072.19130201, 9644.34268454, 8877.83504026, 7799.39529499, 0.])
data['aeroloads_Py'] = np.array([-0., -144.78495637, -124.72748876, -141.81359748, -155.88095732, -168.54159208, -178.87593197, -188.03239815, -226.67105832, -233.77977156, -238.54367701, -315.58006958, -346.98243408, -376.42732276, -428.40371938, -460.55814935, -514.08253305, -608.45781344, -797.59235326, -945.37444233, -1201.2028286, -1460.69809022, -1699.93141114, -1895.13844916, -2173.05749742, -2231.72002959, -2324.32717937, -2368.68132531, -2391.87757363, -2400.12541632, -2397.52018148, -2390.65909733, -2364.75629832, -2319.62320261, -2178.87813088, -2006.41143031, -1828.46265047, -1590.92929559, -887.78584275, -0.])
data['aeroloads_Pz'] = np.array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])
data['r_pts'] = np.array([1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697, 2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267, 8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998, 18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509, 28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998, 42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998, 56.20419994, 58.9374998, 61.67080026, 63.03135])
data['totalCone'] = np.array([2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5])
data['x_az'] = np.array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])
data['y_az'] = np.array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])
data['z_az'] = np.array([1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697, 2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267, 8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998, 18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509, 28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998, 42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998, 56.20419994, 58.9374998, 61.67080026, 63.03135])
data['s'] = np.array([1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697, 2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267, 8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998, 18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509, 28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998, 42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998, 56.20419994, 58.9374998, 61.67080026, 63.03135])
data['r_pts'] = data['r_pts'] + 0.1
npts = np.size(data['z_az'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', RootMoment(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['form'] = 'central'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['step_calc'] = 'relative'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
# check_gradient_unit_test(prob, tol=5e-3)
check_gradient_unit_test(prob, tol=0.05)
def test_TotalLoads(inputs_str, inputs_aeroloads):
npts = np.size(inputs_str['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', TotalLoads(npts), promotes=['*'])
prob = init_IndepVar_add(prob, inputs_str)
prob = init_IndepVar_add(prob, inputs_aeroloads)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, inputs_str)
prob = init_IndepVar_set(prob, inputs_aeroloads)
check_gradient_unit_test(prob, tol=.01, display=False)
def test_OutputsAero():
data = {}
data['AEP_in'] = 23490703.70603071
data['V_in'] = np.array([0.178275041966, 0.190296788962, 0.202568028297, 0.215071993446, 0.227791807017, 0.240710518196, 0.253811139789, 0.267076684744, 0.280490202059, 0.294034811963, 0.307693740274, 0.32145035184, 0.335288182978, 0.349190972822, 0.363142693506, 0.377127579106, 0.391130153286, 0.405135255571, 0.419128066214, 0.43309412959, 0.447019376098, 0.460890142513, 0.474693190784, 0.488415725243, 0.502045408211, 0.515570374001, 0.528979241305, 0.54226112398, 0.555405640233, 0.568402920224, 0.581243612108, 0.593918886545, 0.606420439704, 0.618740494789, 0.630871802149, 0.642807637993, 0.654541801767, 0.666068612246, 0.67738290239, 0.68848001302, 0.699355785386, 0.710006552662, 0.720429130464, 0.730620806428, 0.740579328925, 0.750302894983, 0.759790137479, 0.76904011166, 0.778052281079, 0.786826502994, 0.795363013307, 0.803662411103, 0.811725642861, 0.819553986387, 0.827149034543, 0.834512678823, 0.841647092838, 0.848554715765, 0.855238235811, 0.861700573747, 0.867944866556, 0.873974451244, 0.879792848857, 0.885403748745, 0.890810993109, 0.896018561866, 0.90103055787, 0.90585119251, 0.91048477172, 0.914935682413, 0.919208379383, 0.923307372661, 0.927237215373, 0.931002492094, 0.934607807711, 0.93805777681, 0.941357013592, 0.944510122313, 0.947521688261, 0.950396269264, 0.95313838773, 0.955752523204, 0.95824310546, 0.960614508091, 0.962871042612, 0.965016953057, 0.967056411054, 0.968993511377, 0.970832267946, 0.972576610282, 0.974230380379, 0.975797329995, 0.977281118336, 0.978685310119, 0.980013374, 0.981268681344, 0.982454505324, 0.98357402033, 0.984630301666, 0.985626325531, 0.986564969246, 0.987449011729, 0.988281134189, 0.98906392103, 0.989799860938, 0.990491348148, 0.991140683869, 0.991750077849, 0.992321650072, 0.992857432566, 0.993359371317, 0.99382932827, 0.994269083401, 0.994680336859, 0.995064711163, 0.995423753433, 0.995758937661, 0.996071667003, 0.996363276084, 0.996635033314, 0.996888143197, 0.997123748634, 0.997342933215, 0.997546723483, 0.997736091178, 0.997911955448, 0.998075185027, 0.998226600365, 0.998366975728, 0.998497041241, 0.998617484891, 0.99872895447, 0.998832059475, 0.998927372946, 0.999015433252, 0.999096745819, 0.999171784802, 0.999240994698, 0.999304791902, 0.999363566204, 0.999417682229, 0.999467480822, 0.999513280372, 0.999555378083, 0.999594051191, 0.999629558123, 0.999662139611, 0.999692019741, 0.999719406969, 0.999744495073, 0.999767464065, 0.999788481055, 0.99980770107, 0.999825267834, 0.999841314496, 0.999855964334, 0.999869331401, 0.999881521151, 0.999892631024, 0.999902750988, 0.999911964064, 0.999920346811, 0.999927969778, 0.99993489794, 0.999941191092, 0.999946904234, 0.999952087915, 0.999956788569, 0.99996104882, 0.999964907768, 0.99996840126, 0.999971562138, 0.999974420473, 0.999977003779, 0.999979337216, 0.999981443776, 0.999983344458, 0.999985058426, 0.999986603162, 0.999987994599, 0.999989247253, 0.999990374339, 0.999991387879, 0.999992298805, 0.999993117053, 0.999993851641, 0.999994510758, 0.999995101829, 0.999995631585, 0.999996106123, 0.999996530963, 0.999996911098, 0.999997251045, 0.999997554883, 0.999997826298, 0.999998068616, 0.999998284835, 0.999998477661, 0.999998649529, 0.999998802632])
data['P_in'] = np.array([22025.3984542, 31942.1185436, 42589.4002304, 53993.2085473, 66179.5085267, 79174.2652013, 93003.4436034, 107693.008766, 123268.925721, 139757.159501, 157183.67514, 175574.437668, 194955.41212, 215352.563527, 236791.856922, 259299.257337, 282900.729806, 307622.23936, 333489.751032, 360529.229855, 388766.640862, 418227.949084, 448939.119554, 480926.117306, 514214.907371, 548831.454781, 584801.724571, 622151.681771, 660907.291415, 701094.518535, 742739.328164, 785867.685334, 830505.555078, 876678.902429, 924413.692418, 973735.890078, 1024671.46044, 1077246.36854, 1131486.57941, 1187418.05809, 1245066.76959, 1304458.67896, 1365619.75124, 1428575.95144, 1493353.24461, 1559977.59577, 1628474.96996, 1698871.33222, 1771192.64757, 1845464.88105, 1921713.99769, 1999965.96251, 2080246.74057, 2162582.29688, 2246998.59648, 2333521.6044, 2422177.28568, 2512991.60535, 2605990.52843, 2701200.01997, 2798646.04499, 2898354.56853, 3000351.55562, 3104662.97129, 3211314.78058, 3320332.94852, 3431743.44013, 3545572.22046, 3661284.29229, 3777883.06101, 3895967.17288, 4015492.83297, 4136416.74759, 4258697.48598, 4382294.05017, 4507165.36302, 4631531.67845, 4757024.71876, 4883638.63072, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0])
data['rated_V_in'] = 11.738606532616528
data['rated_Omega_in'] = 12.
data['rated_pitch_in'] = 0.0
data['rated_T_in'] = 714206.5321080858
data['rated_Q_in'] = 3978873.5772973835
data['hub_diameter_in'] = 3.075
data['diameter_in'] = 125.95500453593273
data['max_chord_in'] = 3.0459289459935825
data['V_extreme_in'] = 70.0
data['T_extreme_in'] = 0.0
data['Q_extreme_in'] = 0.0
data['Rtip_in'] = 63.0375
data['precurveTip_in'] = 0.0
data['precsweepTip_in'] = 0.0
prob = Problem()
prob.root = Group()
prob.root.add('comp', OutputsAero(len(data['P_in'])), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_DamageLoads():
data = {}
data['r'] = np.array([2.8667, 5.6, 8.3333, 11.75, 15.85, 19.95, 24.05, 28.15, 32.25, 36.35, 40.45, 44.55, 48.65, 52.75, 56.1667, 58.9, 61.6333]) # (Array): new aerodynamic grid on unit radius
data['rstar'] = np.array([0.000, 0.022, 0.067, 0.111, 0.167, 0.233, 0.300, 0.367, 0.433, 0.500,
0.567, 0.633, 0.700, 0.767, 0.833, 0.889, 0.933, 0.978]) # (Array): nondimensional radial locations of damage equivalent moments
data['Mxb'] = 1e3*np.array([2.3743E+003, 2.0834E+003, 1.8108E+003, 1.5705E+003, 1.3104E+003,
1.0488E+003, 8.2367E+002, 6.3407E+002, 4.7727E+002, 3.4804E+002, 2.4458E+002, 1.6339E+002,
1.0252E+002, 5.7842E+001, 2.7349E+001, 1.1262E+001, 3.8549E+000, 4.4738E-001]) # (Array, N*m): damage equivalent moments about blade c.s. x-direction
data['Myb'] = 1e3*np.array([2.7732E+003, 2.8155E+003, 2.6004E+003, 2.3933E+003, 2.1371E+003,
1.8459E+003, 1.5582E+003, 1.2896E+003, 1.0427E+003, 8.2015E+002, 6.2449E+002, 4.5229E+002,
3.0658E+002, 1.8746E+002, 9.6475E+001, 4.2677E+001, 1.5409E+001, 1.8426E+000]) # (Array, N*m): damage equivalent moments about blade c.s. y-direction
data['theta'] = np.array([13.308, 13.308, 13.308, 13.308, 11.48, 10.162, 9.011, 7.795, 6.544, 5.361, 4.188, 3.125, 2.319, 1.526, 0.863, 0.37, 0.106])
npts = np.size(data['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', DamageLoads(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=5e-5)
def test_SetupRunVarSpeed(inputs_control):
del inputs_control['control_ratedPower'], inputs_control['control_minOmega']
prob = Problem()
prob.root = Group()
prob.root.add('comp', SetupRunVarSpeed(20), promotes=['*'])
prob = init_IndepVar_add(prob, inputs_control)
prob.setup(check=False)
prob = init_IndepVar_set(prob, inputs_control)
check_gradient_unit_test(prob)
def test_SetupRunVarSpeed(inputs_control):
del inputs_control['control_ratedPower'], inputs_control[
'control_minOmega']
prob = Problem()
prob.root = Group()
prob.root.add('comp', SetupRunVarSpeed(20), promotes=['*'])
prob = init_IndepVar_add(prob, inputs_control)
prob.setup(check=False)
prob = init_IndepVar_set(prob, inputs_control)
check_gradient_unit_test(prob)
def test_CCBladeGeometry():
data = {}
data['Rtip'] = 63.0
data['precone'] = 5.0
data['precurveTip'] = 0.0
prob = Problem()
prob.root = Group()
prob.root.add('comp', CCBladeGeometry(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_CCBladeGeometry():
data = {}
data['Rtip'] = 63.0
data['precone'] = 5.0
data['precurveTip'] = 0.0
prob = Problem()
prob.root = Group()
prob.root.add('comp', CCBladeGeometry(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_ExtremeLoads():
data = {}
data['T'] = np.array([2414072.40260361, 188461.59444074])
data['Q'] = np.array([10926313.24295958, -8041330.51312603])
data['nBlades'] = 3
prob = Problem()
prob.root = Group()
prob.root.add('comp', ExtremeLoads(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=1.5e-4)
def test_ExtremeLoads():
data = {}
data['T'] = np.array([2414072.40260361, 188461.59444074])
data['Q']= np.array([10926313.24295958, -8041330.51312603])
data['nBlades'] = 3
prob = Problem()
prob.root = Group()
prob.root.add('comp', ExtremeLoads(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=1.5e-4)
def test_GustETM():
data = {}
data['V_mean'] = 10.0
data['V_hub'] = 11.7733866478
data['std'] = 3
data['turbulence_class'] = TURBULENCE_CLASS['B']
prob = Problem()
prob.root = Group()
prob.root.add('comp', GustETM(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_CSMDrivetrain(aeroPower, ratedPower):
data = {}
data['aeroPower'] = aeroPower
data['ratedPower'] = ratedPower
data['drivetrainType'] = DRIVETRAIN_TYPE['GEARED']
prob = Problem()
prob.root = Group()
prob.root.add('comp', CSMDrivetrain(len(aeroPower)), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=6e-4)
def test_CSMDrivetrain(aeroPower, ratedPower):
data = {}
data['aeroPower'] = aeroPower
data['ratedPower'] = ratedPower
data['drivetrainType'] = DRIVETRAIN_TYPE['GEARED']
prob = Problem()
prob.root = Group()
prob.root.add('comp', CSMDrivetrain(len(aeroPower)), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=6e-4)
def test_GustETM():
data = {}
data['V_mean'] = 10.0
data['V_hub'] = 11.7733866478
data['std'] = 3
data['turbulence_class'] = TURBULENCE_CLASS['B']
prob = Problem()
prob.root = Group()
prob.root.add('comp', GustETM(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_BladeCurvature():
data = {}
data['r'] = np.array([1.5375, 1.84056613137, 1.93905987557, 2.03755361977, 2.14220322299, 2.24069696719, 2.33919071139, 2.90419974, 3.04095863882, 3.13945238302, 5.637500205, 7.04226699698, 8.370800055, 11.8494282928, 13.8375, 14.7182548841, 15.887499795, 18.5537233505, 19.9875, 22.0564071286, 24.087500205, 26.16236509, 28.187499795, 32.2875, 33.5678634821, 36.387500205, 38.5725768593, 40.487499795, 40.6967540173, 40.7964789782, 40.8975, 42.6919644014, 44.5875, 52.787499795, 56.204199945, 58.937499795, 61.67080026, 63.0375])
data['precurve'] = np.array([0.0, 0.043324361025, 0.0573893371698, 0.0714469497372, 0.0863751069858, 0.100417566593, 0.114452695996, 0.194824077331, 0.214241777459, 0.228217752953, 0.580295739194, 0.776308800624, 0.960411633829, 1.4368012564, 1.7055214864, 1.823777005, 1.98003324362, 2.32762426752, 2.50856911855, 2.76432512112, 3.0113656418, 3.26199245912, 3.50723775206, 4.0150233695, 4.17901272929, 4.55356019347, 4.85962948702, 5.14086873143, 5.17214747287, 5.18708601127, 5.20223968442, 5.47491847385, 5.77007321175, 7.12818875977, 7.7314427824, 8.22913789456, 8.73985955154, 9.0])
data['presweep'] = np.array([0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
data['precone'] = 2.5
npts = np.size(data['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', BladeCurvature(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=5e-5)
def test_RayleighCDF():
data = {}
data['xbar'] = 6.0
data['x'] = np.array([3.0, 3.08805770406, 3.17611540812, 3.26417311218, 3.35223081624, 3.4402885203, 3.52834622436, 3.61640392842, 3.70446163248, 3.79251933654, 3.8805770406, 3.96863474466, 4.05669244872, 4.14475015278, 4.23280785684, 4.3208655609, 4.40892326496, 4.49698096902, 4.58503867308, 4.67309637714, 4.7611540812, 4.84921178526, 4.93726948932, 5.02532719338, 5.11338489744, 5.2014426015, 5.28950030556, 5.37755800962, 5.46561571369, 5.55367341775, 5.64173112181, 5.72978882587, 5.81784652993, 5.90590423399, 5.99396193805, 6.08201964211, 6.17007734617, 6.25813505023, 6.34619275429, 6.43425045835, 6.52230816241, 6.61036586647, 6.69842357053, 6.78648127459, 6.87453897865, 6.96259668271, 7.05065438677, 7.13871209083, 7.22676979489, 7.31482749895, 7.40288520301, 7.49094290707, 7.57900061113, 7.66705831519, 7.75511601925, 7.84317372331, 7.93123142737, 8.01928913143, 8.10734683549, 8.19540453955, 8.28346224361, 8.37151994767, 8.45957765173, 8.54763535579, 8.63569305985, 8.72375076391, 8.81180846797, 8.89986617203, 8.98792387609, 9.07598158015, 9.16403928421, 9.25209698827, 9.34015469233, 9.42821239639, 9.51627010045, 9.60432780451, 9.69238550857, 9.78044321263, 9.86850091669, 9.95655862075, 10.0446163248, 10.1326740289, 10.2207317329, 10.308789437, 10.3968471411, 10.4849048451, 10.5729625492, 10.6610202532, 10.7490779573, 10.8371356614, 10.9251933654, 11.0132510695, 11.1013087735, 11.1893664776, 11.2774241817, 11.3654818857, 11.4535395898, 11.5415972938, 11.6296549979, 11.717712702, 11.8505355749, 11.9833584479, 12.1161813209, 12.2490041939, 12.3818270669, 12.5146499398, 12.6474728128, 12.7802956858, 12.9131185588, 13.0459414318, 13.1787643047, 13.3115871777, 13.4444100507, 13.5772329237, 13.7100557967, 13.8428786696, 13.9757015426, 14.1085244156, 14.2413472886, 14.3741701616, 14.5069930345, 14.6398159075, 14.7726387805, 14.9054616535, 15.0382845265, 15.1711073994, 15.3039302724, 15.4367531454, 15.5695760184, 15.7023988914, 15.8352217644, 15.9680446373, 16.1008675103, 16.2336903833, 16.3665132563, 16.4993361293, 16.6321590022, 16.7649818752, 16.8978047482, 17.0306276212, 17.1634504942, 17.2962733671, 17.4290962401, 17.5619191131, 17.6947419861, 17.8275648591, 17.960387732, 18.093210605, 18.226033478, 18.358856351, 18.491679224, 18.6245020969, 18.7573249699, 18.8901478429, 19.0229707159, 19.1557935889, 19.2886164618, 19.4214393348, 19.5542622078, 19.6870850808, 19.8199079538, 19.9527308267, 20.0855536997, 20.2183765727, 20.3511994457, 20.4840223187, 20.6168451916, 20.7496680646, 20.8824909376, 21.0153138106, 21.1481366836, 21.2809595565, 21.4137824295, 21.5466053025, 21.6794281755, 21.8122510485, 21.9450739215, 22.0778967944, 22.2107196674, 22.3435425404, 22.4763654134, 22.6091882864, 22.7420111593, 22.8748340323, 23.0076569053, 23.1404797783, 23.2733026513, 23.4061255242, 23.5389483972, 23.6717712702, 23.8045941432, 23.9374170162, 24.0702398891, 24.2030627621, 24.3358856351, 24.4687085081, 24.6015313811, 24.734354254, 24.867177127, 25.0])
n = len(data['x'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', RayleighCDF(n), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test( prob, tol=2e-6)
def test_SetupPCModVarSpeed():
data = {}
data['control_tsr'] = 8.0
data['control_pitch'] = 1.0
data['Vrated'] = 12.0
data['R'] = 63.0
data['Vfactor'] = 0.7
prob = Problem()
prob.root = Group()
prob.root.add('comp', SetupPCModVarSpeed(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_SetupPCModVarSpeed():
data = {}
data['control_tsr'] = 8.0
data['control_pitch'] = 1.0
data['Vrated'] = 12.0
data['R'] = 63.0
data['Vfactor'] = 0.7
prob = Problem()
prob.root = Group()
prob.root.add('comp', SetupPCModVarSpeed(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_WeibullWithMeanCDF():
data = {}
data['xbar'] = 8.0
data['k'] = 2.2
data['x'] = np.linspace(1.0, 15.0, 50)
n = len(data['x'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', WeibullWithMeanCDF(n), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_WeibullWithMeanCDF():
data = {}
data['xbar'] = 8.0
data['k'] = 2.2
data['x'] = np.linspace(1.0, 15.0, 50)
n = len(data['x'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', WeibullWithMeanCDF(n), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_MassProperties():
data = {}
data['blade_mass'] = 17288.717087
data['blade_moment_of_inertia'] = 11634376.0531
data['tilt'] = 5.0
data['nBlades'] = 3
prob = Problem()
prob.root = Group()
prob.root.add('comp', MassProperties(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_MassProperties():
data = {}
data['blade_mass'] = 17288.717087
data['blade_moment_of_inertia'] = 11634376.0531
data['tilt'] = 5.0
data['nBlades'] = 3
prob = Problem()
prob.root = Group()
prob.root.add('comp', MassProperties(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_TipDeflection():
data = {}
data['dx'] = 4.27242809591
data['dy'] = -0.371550675139
data['dz'] = 0.0400553989266
data['theta'] = -0.0878099
data['pitch'] = 0.0
data['azimuth'] = 180.0
data['tilt'] = 5.0
data['totalConeTip'] = 2.5
data['dynamicFactor'] = 1.2
prob = Problem()
prob.root = Group()
prob.root.add('comp', TipDeflection(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_TipDeflection():
data = {}
data['dx'] = 4.27242809591
data['dy'] = -0.371550675139
data['dz'] = 0.0400553989266
data['theta'] = -0.0878099
data['pitch'] = 0.0
data['azimuth'] = 180.0
data['tilt'] = 5.0
data['totalConeTip'] = 2.5
data['dynamicFactor'] = 1.2
prob = Problem()
prob.root = Group()
prob.root.add('comp', TipDeflection(), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_BladeCurvature():
data = {}
data['r'] = np.array([
1.5375, 1.84056613137, 1.93905987557, 2.03755361977, 2.14220322299,
2.24069696719, 2.33919071139, 2.90419974, 3.04095863882, 3.13945238302,
5.637500205, 7.04226699698, 8.370800055, 11.8494282928, 13.8375,
14.7182548841, 15.887499795, 18.5537233505, 19.9875, 22.0564071286,
24.087500205, 26.16236509, 28.187499795, 32.2875, 33.5678634821,
36.387500205, 38.5725768593, 40.487499795, 40.6967540173,
40.7964789782, 40.8975, 42.6919644014, 44.5875, 52.787499795,
56.204199945, 58.937499795, 61.67080026, 63.0375
])
data['precurve'] = np.array([
0.0, 0.043324361025, 0.0573893371698, 0.0714469497372, 0.0863751069858,
0.100417566593, 0.114452695996, 0.194824077331, 0.214241777459,
0.228217752953, 0.580295739194, 0.776308800624, 0.960411633829,
1.4368012564, 1.7055214864, 1.823777005, 1.98003324362, 2.32762426752,
2.50856911855, 2.76432512112, 3.0113656418, 3.26199245912,
3.50723775206, 4.0150233695, 4.17901272929, 4.55356019347,
4.85962948702, 5.14086873143, 5.17214747287, 5.18708601127,
5.20223968442, 5.47491847385, 5.77007321175, 7.12818875977,
7.7314427824, 8.22913789456, 8.73985955154, 9.0
])
data['presweep'] = np.array([
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
])
data['precone'] = 2.5
npts = np.size(data['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', BladeCurvature(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=5e-5)
def test_DamageLoads():
data = {}
data['r'] = np.array([
2.8667, 5.6, 8.3333, 11.75, 15.85, 19.95, 24.05, 28.15, 32.25, 36.35,
40.45, 44.55, 48.65, 52.75, 56.1667, 58.9, 61.6333
]) # (Array): new aerodynamic grid on unit radius
data['rstar'] = np.array([
0.000, 0.022, 0.067, 0.111, 0.167, 0.233, 0.300, 0.367, 0.433, 0.500,
0.567, 0.633, 0.700, 0.767, 0.833, 0.889, 0.933, 0.978
]) # (Array): nondimensional radial locations of damage equivalent moments
data['Mxb'] = 1e3 * np.array([
2.3743E+003, 2.0834E+003, 1.8108E+003, 1.5705E+003, 1.3104E+003,
1.0488E+003, 8.2367E+002, 6.3407E+002, 4.7727E+002, 3.4804E+002,
2.4458E+002, 1.6339E+002, 1.0252E+002, 5.7842E+001, 2.7349E+001,
1.1262E+001, 3.8549E+000, 4.4738E-001
]) # (Array, N*m): damage equivalent moments about blade c.s. x-direction
data['Myb'] = 1e3 * np.array([
2.7732E+003, 2.8155E+003, 2.6004E+003, 2.3933E+003, 2.1371E+003,
1.8459E+003, 1.5582E+003, 1.2896E+003, 1.0427E+003, 8.2015E+002,
6.2449E+002, 4.5229E+002, 3.0658E+002, 1.8746E+002, 9.6475E+001,
4.2677E+001, 1.5409E+001, 1.8426E+000
]) # (Array, N*m): damage equivalent moments about blade c.s. y-direction
data['theta'] = np.array([
13.308, 13.308, 13.308, 13.308, 11.48, 10.162, 9.011, 7.795, 6.544,
5.361, 4.188, 3.125, 2.319, 1.526, 0.863, 0.37, 0.106
])
npts = np.size(data['r'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', DamageLoads(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=5e-5)
def test_BladeDeflection():
NINPUT = 5
data = {}
data['dx'] = np.array([0.0, 1.39292987639e-05, 2.42964362361e-05, 3.73611415086e-05, 5.41276445811e-05, 7.25773833888e-05, 9.35928069108e-05, 0.000268596883961, 0.000326596195925, 0.000372259193056, 0.00253285667923, 0.00455502505195, 0.00706548815047, 0.017530829505, 0.0272151750396, 0.0326222002275, 0.0410034329532, 0.0655343595834, 0.0818962168404, 0.109404026845, 0.140873579829, 0.177728457961, 0.218834564559, 0.318921553007, 0.354759106888, 0.442199496033, 0.51925808593, 0.593774567144, 0.602307030873, 0.606403310701, 0.610574803156, 0.688594372351, 0.779039709944, 1.2575016857, 1.49348304065, 1.69400834707, 1.90049599542, 2.00437890947])
data['dy'] = np.array([0.0, -9.11673273998e-07, -1.59100573731e-06, -2.44958708331e-06, -3.55417484862e-06, -4.77185559999e-06, -6.16078184177e-06, -1.77741315033e-05, -2.16387627714e-05, -2.46907649535e-05, -0.000177142122562, -0.00032886393339, -0.000523351402122, -0.00136667950935, -0.00215626107629, -0.00258384925494, -0.00322230578732, -0.00497009730837, -0.00607094164543, -0.00784913963514, -0.00980572470543, -0.0120149227037, -0.0143773305481, -0.0197932419212, -0.0216584171821, -0.0260612290834, -0.0297556792352, -0.0331954636074, -0.0335827012855, -0.0337680683508, -0.0339563913295, -0.0373950218981, -0.0412189884014, -0.0597108654077, -0.0681209435104, -0.0750137913709, -0.0819682184214, -0.085450339495])
data['dz'] = np.array([0.0, 0.000190681386865, 0.000249594803444, 0.000305808888278, 0.000363037181592, 0.000414559803159, 0.000465006717172, 0.00079021115147, 0.000878140045153, 0.000937905765842, 0.00236994444878, 0.003141201122, 0.00385188728597, 0.0059090995974, 0.00721749042255, 0.00775193614485, 0.00838246793108, 0.00965065210961, 0.010256099836, 0.0110346081014, 0.0117147556349, 0.0123373351129, 0.0128834601832, 0.0138847141125, 0.014168689932, 0.0147331399477, 0.0151335663707, 0.015450772371, 0.015480903206, 0.0154926033582, 0.0155017716157, 0.0156218668894, 0.0157089720273, 0.0159512118376, 0.0160321965202, 0.0160695719649, 0.0160814363339, 0.0160814363339])
data['pitch'] = 0.0
data['theta'] = np.array([13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 12.9342399734, 12.4835173655, 11.4807910924, 10.9555299481, 10.2141623046, 9.50473135323, 8.7980712345, 8.12522342051, 6.81383731475, 6.42068577363, 5.58414615907, 4.96394315694, 4.44088253161, 4.38489418276, 4.35829308634, 4.33139992746, 3.86272702658, 3.38639207628, 1.57773054352, 0.953410121155, 0.504987738102, 0.0995174088527, -0.0878099])
data['r_in0'] = np.array([1.5375, 13.2620872, 17.35321975, 40.19535987, 63.0375])
data['Rhub0'] = 1.5375
data['r_pts0'] = np.array([1.5375, 1.84056613137, 1.93905987557, 2.03755361977, 2.14220322299, 2.24069696719, 2.33919071139, 2.90419974, 3.04095863882, 3.13945238302, 5.637500205, 7.04226699698, 8.370800055, 11.8494282928, 13.8375, 14.7182548841, 15.887499795, 18.5537233505, 19.9875, 22.0564071286, 24.087500205, 26.16236509, 28.187499795, 32.2875, 33.5678634821, 36.387500205, 38.5725768593, 40.487499795, 40.6967540173, 40.7964789782, 40.8975, 42.6919644014, 44.5875, 52.787499795, 56.204199945, 58.937499795, 61.67080026, 63.0375])
data['precurve0'] = np.zeros_like(data['dx']) #np.linspace(0.0, 5.0, np.size(data['dx']))
data['bladeLength0'] = 61.5
npts = np.size(data['dx'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', BladeDeflection(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=1e-4)
def test_RayleighCDF():
data = {}
data['xbar'] = 6.0
data['x'] = np.array([
3.0, 3.08805770406, 3.17611540812, 3.26417311218, 3.35223081624,
3.4402885203, 3.52834622436, 3.61640392842, 3.70446163248,
3.79251933654, 3.8805770406, 3.96863474466, 4.05669244872,
4.14475015278, 4.23280785684, 4.3208655609, 4.40892326496,
4.49698096902, 4.58503867308, 4.67309637714, 4.7611540812,
4.84921178526, 4.93726948932, 5.02532719338, 5.11338489744,
5.2014426015, 5.28950030556, 5.37755800962, 5.46561571369,
5.55367341775, 5.64173112181, 5.72978882587, 5.81784652993,
5.90590423399, 5.99396193805, 6.08201964211, 6.17007734617,
6.25813505023, 6.34619275429, 6.43425045835, 6.52230816241,
6.61036586647, 6.69842357053, 6.78648127459, 6.87453897865,
6.96259668271, 7.05065438677, 7.13871209083, 7.22676979489,
7.31482749895, 7.40288520301, 7.49094290707, 7.57900061113,
7.66705831519, 7.75511601925, 7.84317372331, 7.93123142737,
8.01928913143, 8.10734683549, 8.19540453955, 8.28346224361,
8.37151994767, 8.45957765173, 8.54763535579, 8.63569305985,
8.72375076391, 8.81180846797, 8.89986617203, 8.98792387609,
9.07598158015, 9.16403928421, 9.25209698827, 9.34015469233,
9.42821239639, 9.51627010045, 9.60432780451, 9.69238550857,
9.78044321263, 9.86850091669, 9.95655862075, 10.0446163248,
10.1326740289, 10.2207317329, 10.308789437, 10.3968471411,
10.4849048451, 10.5729625492, 10.6610202532, 10.7490779573,
10.8371356614, 10.9251933654, 11.0132510695, 11.1013087735,
11.1893664776, 11.2774241817, 11.3654818857, 11.4535395898,
11.5415972938, 11.6296549979, 11.717712702, 11.8505355749,
11.9833584479, 12.1161813209, 12.2490041939, 12.3818270669,
12.5146499398, 12.6474728128, 12.7802956858, 12.9131185588,
13.0459414318, 13.1787643047, 13.3115871777, 13.4444100507,
13.5772329237, 13.7100557967, 13.8428786696, 13.9757015426,
14.1085244156, 14.2413472886, 14.3741701616, 14.5069930345,
14.6398159075, 14.7726387805, 14.9054616535, 15.0382845265,
15.1711073994, 15.3039302724, 15.4367531454, 15.5695760184,
15.7023988914, 15.8352217644, 15.9680446373, 16.1008675103,
16.2336903833, 16.3665132563, 16.4993361293, 16.6321590022,
16.7649818752, 16.8978047482, 17.0306276212, 17.1634504942,
17.2962733671, 17.4290962401, 17.5619191131, 17.6947419861,
17.8275648591, 17.960387732, 18.093210605, 18.226033478, 18.358856351,
18.491679224, 18.6245020969, 18.7573249699, 18.8901478429,
19.0229707159, 19.1557935889, 19.2886164618, 19.4214393348,
19.5542622078, 19.6870850808, 19.8199079538, 19.9527308267,
20.0855536997, 20.2183765727, 20.3511994457, 20.4840223187,
20.6168451916, 20.7496680646, 20.8824909376, 21.0153138106,
21.1481366836, 21.2809595565, 21.4137824295, 21.5466053025,
21.6794281755, 21.8122510485, 21.9450739215, 22.0778967944,
22.2107196674, 22.3435425404, 22.4763654134, 22.6091882864,
22.7420111593, 22.8748340323, 23.0076569053, 23.1404797783,
23.2733026513, 23.4061255242, 23.5389483972, 23.6717712702,
23.8045941432, 23.9374170162, 24.0702398891, 24.2030627621,
24.3358856351, 24.4687085081, 24.6015313811, 24.734354254,
24.867177127, 25.0
])
n = len(data['x'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', RayleighCDF(n), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=2e-6)
def test_TestRootMoment():
data_aero = inputs_aeroloads_operating()
data_str = inputs_str()
data = {}
data['aeroloads_r'] = np.array([
1.5375, 1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322,
2.24069697, 2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021,
7.042267, 8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998,
18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509,
28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998,
42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998,
56.20419994, 58.9374998, 61.67080026, 63.03135, 63.0375
])
data['aeroloads_Px'] = np.array([
0., 610.90043568, 469.20341791, 529.97827499, 580.02866009,
625.05962708, 661.78185992, 694.27930979, 830.93330879, 856.02758285,
872.85385017, 1155.30591175, 1278.65677079, 1394.84126103,
1598.03834945, 1726.35834871, 1933.87963881, 2271.78104859,
2748.8320438, 3083.27106663, 3677.62201087, 4343.75465079,
5051.54596049, 5731.9904625, 7004.11264534, 7364.31723123,
8093.26478732, 8599.46800605, 9000.96202336, 9390.16142378,
9555.38019573, 9708.33654591, 9964.10107999, 10160.41657631,
10284.8330047, 10072.19130201, 9644.34268454, 8877.83504026,
7799.39529499, 0.
])
data['aeroloads_Py'] = np.array([
-0., -144.78495637, -124.72748876, -141.81359748, -155.88095732,
-168.54159208, -178.87593197, -188.03239815, -226.67105832,
-233.77977156, -238.54367701, -315.58006958, -346.98243408,
-376.42732276, -428.40371938, -460.55814935, -514.08253305,
-608.45781344, -797.59235326, -945.37444233, -1201.2028286,
-1460.69809022, -1699.93141114, -1895.13844916, -2173.05749742,
-2231.72002959, -2324.32717937, -2368.68132531, -2391.87757363,
-2400.12541632, -2397.52018148, -2390.65909733, -2364.75629832,
-2319.62320261, -2178.87813088, -2006.41143031, -1828.46265047,
-1590.92929559, -887.78584275, -0.
])
data['aeroloads_Pz'] = np.array([
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0.
])
data['r_pts'] = np.array([
1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697,
2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267,
8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998,
18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509,
28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998,
42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998,
56.20419994, 58.9374998, 61.67080026, 63.03135
])
data['totalCone'] = np.array([
2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5,
2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5,
2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5
])
data['x_az'] = np.array([
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0.
])
data['y_az'] = np.array([
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
0., 0.
])
data['z_az'] = np.array([
1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697,
2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267,
8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998,
18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509,
28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998,
42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998,
56.20419994, 58.9374998, 61.67080026, 63.03135
])
data['s'] = np.array([
1.54365, 1.84056613, 1.93905988, 2.03755362, 2.14220322, 2.24069697,
2.33919071, 2.90419974, 3.04095864, 3.13945238, 5.63750021, 7.042267,
8.37080006, 10.54495078, 11.7875002, 13.54900997, 15.8874998,
18.55372335, 19.9875, 22.05640713, 24.08750021, 26.16236509,
28.1874998, 32.2875, 33.56786348, 36.38750021, 38.57257686, 40.4874998,
42.58004108, 43.57729024, 44.5875, 46.58134943, 48.6875002, 52.7874998,
56.20419994, 58.9374998, 61.67080026, 63.03135
])
data['r_pts'] = data['r_pts'] + 0.1
npts = np.size(data['z_az'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', RootMoment(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['form'] = 'central'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.root.comp.deriv_options['step_calc'] = 'relative'
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
# check_gradient_unit_test(prob, tol=5e-3)
check_gradient_unit_test(prob, tol=0.05)
def test_OutputsAero():
data = {}
data['AEP_in'] = 23490703.70603071
data['V_in'] = np.array([
0.178275041966, 0.190296788962, 0.202568028297, 0.215071993446,
0.227791807017, 0.240710518196, 0.253811139789, 0.267076684744,
0.280490202059, 0.294034811963, 0.307693740274, 0.32145035184,
0.335288182978, 0.349190972822, 0.363142693506, 0.377127579106,
0.391130153286, 0.405135255571, 0.419128066214, 0.43309412959,
0.447019376098, 0.460890142513, 0.474693190784, 0.488415725243,
0.502045408211, 0.515570374001, 0.528979241305, 0.54226112398,
0.555405640233, 0.568402920224, 0.581243612108, 0.593918886545,
0.606420439704, 0.618740494789, 0.630871802149, 0.642807637993,
0.654541801767, 0.666068612246, 0.67738290239, 0.68848001302,
0.699355785386, 0.710006552662, 0.720429130464, 0.730620806428,
0.740579328925, 0.750302894983, 0.759790137479, 0.76904011166,
0.778052281079, 0.786826502994, 0.795363013307, 0.803662411103,
0.811725642861, 0.819553986387, 0.827149034543, 0.834512678823,
0.841647092838, 0.848554715765, 0.855238235811, 0.861700573747,
0.867944866556, 0.873974451244, 0.879792848857, 0.885403748745,
0.890810993109, 0.896018561866, 0.90103055787, 0.90585119251,
0.91048477172, 0.914935682413, 0.919208379383, 0.923307372661,
0.927237215373, 0.931002492094, 0.934607807711, 0.93805777681,
0.941357013592, 0.944510122313, 0.947521688261, 0.950396269264,
0.95313838773, 0.955752523204, 0.95824310546, 0.960614508091,
0.962871042612, 0.965016953057, 0.967056411054, 0.968993511377,
0.970832267946, 0.972576610282, 0.974230380379, 0.975797329995,
0.977281118336, 0.978685310119, 0.980013374, 0.981268681344,
0.982454505324, 0.98357402033, 0.984630301666, 0.985626325531,
0.986564969246, 0.987449011729, 0.988281134189, 0.98906392103,
0.989799860938, 0.990491348148, 0.991140683869, 0.991750077849,
0.992321650072, 0.992857432566, 0.993359371317, 0.99382932827,
0.994269083401, 0.994680336859, 0.995064711163, 0.995423753433,
0.995758937661, 0.996071667003, 0.996363276084, 0.996635033314,
0.996888143197, 0.997123748634, 0.997342933215, 0.997546723483,
0.997736091178, 0.997911955448, 0.998075185027, 0.998226600365,
0.998366975728, 0.998497041241, 0.998617484891, 0.99872895447,
0.998832059475, 0.998927372946, 0.999015433252, 0.999096745819,
0.999171784802, 0.999240994698, 0.999304791902, 0.999363566204,
0.999417682229, 0.999467480822, 0.999513280372, 0.999555378083,
0.999594051191, 0.999629558123, 0.999662139611, 0.999692019741,
0.999719406969, 0.999744495073, 0.999767464065, 0.999788481055,
0.99980770107, 0.999825267834, 0.999841314496, 0.999855964334,
0.999869331401, 0.999881521151, 0.999892631024, 0.999902750988,
0.999911964064, 0.999920346811, 0.999927969778, 0.99993489794,
0.999941191092, 0.999946904234, 0.999952087915, 0.999956788569,
0.99996104882, 0.999964907768, 0.99996840126, 0.999971562138,
0.999974420473, 0.999977003779, 0.999979337216, 0.999981443776,
0.999983344458, 0.999985058426, 0.999986603162, 0.999987994599,
0.999989247253, 0.999990374339, 0.999991387879, 0.999992298805,
0.999993117053, 0.999993851641, 0.999994510758, 0.999995101829,
0.999995631585, 0.999996106123, 0.999996530963, 0.999996911098,
0.999997251045, 0.999997554883, 0.999997826298, 0.999998068616,
0.999998284835, 0.999998477661, 0.999998649529, 0.999998802632
])
data['P_in'] = np.array([
22025.3984542, 31942.1185436, 42589.4002304, 53993.2085473,
66179.5085267, 79174.2652013, 93003.4436034, 107693.008766,
123268.925721, 139757.159501, 157183.67514, 175574.437668,
194955.41212, 215352.563527, 236791.856922, 259299.257337,
282900.729806, 307622.23936, 333489.751032, 360529.229855,
388766.640862, 418227.949084, 448939.119554, 480926.117306,
514214.907371, 548831.454781, 584801.724571, 622151.681771,
660907.291415, 701094.518535, 742739.328164, 785867.685334,
830505.555078, 876678.902429, 924413.692418, 973735.890078,
1024671.46044, 1077246.36854, 1131486.57941, 1187418.05809,
1245066.76959, 1304458.67896, 1365619.75124, 1428575.95144,
1493353.24461, 1559977.59577, 1628474.96996, 1698871.33222,
1771192.64757, 1845464.88105, 1921713.99769, 1999965.96251,
2080246.74057, 2162582.29688, 2246998.59648, 2333521.6044,
2422177.28568, 2512991.60535, 2605990.52843, 2701200.01997,
2798646.04499, 2898354.56853, 3000351.55562, 3104662.97129,
3211314.78058, 3320332.94852, 3431743.44013, 3545572.22046,
3661284.29229, 3777883.06101, 3895967.17288, 4015492.83297,
4136416.74759, 4258697.48598, 4382294.05017, 4507165.36302,
4631531.67845, 4757024.71876, 4883638.63072, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0,
5000000.0, 5000000.0, 5000000.0, 5000000.0, 5000000.0
])
data['rated_V_in'] = 11.738606532616528
data['rated_Omega_in'] = 12.
data['rated_pitch_in'] = 0.0
data['rated_T_in'] = 714206.5321080858
data['rated_Q_in'] = 3978873.5772973835
data['hub_diameter_in'] = 3.075
data['diameter_in'] = 125.95500453593273
data['max_chord_in'] = 3.0459289459935825
data['V_extreme_in'] = 70.0
data['T_extreme_in'] = 0.0
data['Q_extreme_in'] = 0.0
data['Rtip_in'] = 63.0375
data['precurveTip_in'] = 0.0
data['precsweepTip_in'] = 0.0
prob = Problem()
prob.root = Group()
prob.root.add('comp', OutputsAero(len(data['P_in'])), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob)
def test_CCBlade(bem_case):
data = {}
data['r'] = np.array([
2.8667, 5.6000, 8.3333, 11.7500, 15.8500, 19.9500, 24.0500, 28.1500,
32.2500, 36.3500, 40.4500, 44.5500, 48.6500, 52.7500, 56.1667, 58.9000,
61.6333
])
data['chord'] = np.array([
3.542, 3.854, 4.167, 4.557, 4.652, 4.458, 4.249, 4.007, 3.748, 3.502,
3.256, 3.010, 2.764, 2.518, 2.313, 2.086, 1.419
])
data['theta'] = np.array([
13.308, 13.308, 13.308, 13.308, 11.480, 10.162, 9.011, 7.795, 6.544,
5.361, 4.188, 3.125, 2.319, 1.526, 0.863, 0.370, 0.106
])
data['Rhub'] = 1.5
data['Rtip'] = 63.0
data['hubHt'] = 80.0
data['precone'] = 2.5
data['tilt'] = -5.0
data['yaw'] = 0.0
data['B'] = 3
data['rho'] = 1.225
data['mu'] = 1.81206e-5
data['shearExp'] = 0.2
data['nSector'] = 4
data['precurve'] = np.zeros(len(data['r']))
data['precurveTip'] = 0.0
# airfoils
basepath = os.path.join(
os.path.join(os.path.dirname(os.path.realpath(__file__)), '..'),
'5MW_AFFiles/')
# load all airfoils
airfoil_types = [0] * 8
airfoil_types[0] = basepath + 'Cylinder1.dat'
airfoil_types[1] = basepath + 'Cylinder2.dat'
airfoil_types[2] = basepath + 'DU40_A17.dat'
airfoil_types[3] = basepath + 'DU35_A17.dat'
airfoil_types[4] = basepath + 'DU30_A17.dat'
airfoil_types[5] = basepath + 'DU25_A17.dat'
airfoil_types[6] = basepath + 'DU21_A17.dat'
airfoil_types[7] = basepath + 'NACA64_A17.dat'
# place at appropriate radial stations
af_idx = [0, 0, 1, 2, 3, 3, 4, 5, 5, 6, 6, 7, 7, 7, 7, 7, 7]
n = len(data['r'])
af = [0] * n
for i in range(n):
af[i] = airfoil_types[af_idx[i]]
data['airfoil_files'] = af
prob = Problem()
prob.root = Group()
if bem_case == 'power':
# run_case = 'power'
data['Uhub'] = np.array([
3.0, 4.15789473684, 5.31578947368, 6.47368421053, 7.63157894737,
8.78947368421, 9.94736842105, 11.1052631579, 12.2631578947,
13.4210526316, 14.5789473684, 15.7368421053, 16.8947368421,
18.0526315789, 19.2105263158, 20.3684210526, 21.5263157895,
22.6842105263, 23.8421052632, 25.0
])
data['Omega'] = np.array([
3.43647024491, 4.76282718154, 6.08918411817, 7.41554105481,
8.74189799144, 10.0682549281, 11.3946118647, 12.0, 12.0, 12.0,
12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0
])
data['pitch'] = np.array([
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
])
n2 = len(data['Uhub'])
prob.root.add('comp', CCBladePower(n, n2), promotes=['*'])
elif bem_case == 'loads':
# run_case = 'loads'
data['V_load'] = 12.0
data['Omega_load'] = 10.0
data['pitch_load'] = 0.0
data['azimuth_load'] = 180.0
n2 = 1
prob.root.add('comp', CCBladeLoads(n, n2), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=1e-3, display=True)
def test_BladeDeflection():
NINPUT = 5
data = {}
data['dx'] = np.array([
0.0, 1.39292987639e-05, 2.42964362361e-05, 3.73611415086e-05,
5.41276445811e-05, 7.25773833888e-05, 9.35928069108e-05,
0.000268596883961, 0.000326596195925, 0.000372259193056,
0.00253285667923, 0.00455502505195, 0.00706548815047, 0.017530829505,
0.0272151750396, 0.0326222002275, 0.0410034329532, 0.0655343595834,
0.0818962168404, 0.109404026845, 0.140873579829, 0.177728457961,
0.218834564559, 0.318921553007, 0.354759106888, 0.442199496033,
0.51925808593, 0.593774567144, 0.602307030873, 0.606403310701,
0.610574803156, 0.688594372351, 0.779039709944, 1.2575016857,
1.49348304065, 1.69400834707, 1.90049599542, 2.00437890947
])
data['dy'] = np.array([
0.0, -9.11673273998e-07, -1.59100573731e-06, -2.44958708331e-06,
-3.55417484862e-06, -4.77185559999e-06, -6.16078184177e-06,
-1.77741315033e-05, -2.16387627714e-05, -2.46907649535e-05,
-0.000177142122562, -0.00032886393339, -0.000523351402122,
-0.00136667950935, -0.00215626107629, -0.00258384925494,
-0.00322230578732, -0.00497009730837, -0.00607094164543,
-0.00784913963514, -0.00980572470543, -0.0120149227037,
-0.0143773305481, -0.0197932419212, -0.0216584171821, -0.0260612290834,
-0.0297556792352, -0.0331954636074, -0.0335827012855, -0.0337680683508,
-0.0339563913295, -0.0373950218981, -0.0412189884014, -0.0597108654077,
-0.0681209435104, -0.0750137913709, -0.0819682184214, -0.085450339495
])
data['dz'] = np.array([
0.0, 0.000190681386865, 0.000249594803444, 0.000305808888278,
0.000363037181592, 0.000414559803159, 0.000465006717172,
0.00079021115147, 0.000878140045153, 0.000937905765842,
0.00236994444878, 0.003141201122, 0.00385188728597, 0.0059090995974,
0.00721749042255, 0.00775193614485, 0.00838246793108, 0.00965065210961,
0.010256099836, 0.0110346081014, 0.0117147556349, 0.0123373351129,
0.0128834601832, 0.0138847141125, 0.014168689932, 0.0147331399477,
0.0151335663707, 0.015450772371, 0.015480903206, 0.0154926033582,
0.0155017716157, 0.0156218668894, 0.0157089720273, 0.0159512118376,
0.0160321965202, 0.0160695719649, 0.0160814363339, 0.0160814363339
])
data['pitch'] = 0.0
data['theta'] = np.array([
13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783,
13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783, 13.2783,
12.9342399734, 12.4835173655, 11.4807910924, 10.9555299481,
10.2141623046, 9.50473135323, 8.7980712345, 8.12522342051,
6.81383731475, 6.42068577363, 5.58414615907, 4.96394315694,
4.44088253161, 4.38489418276, 4.35829308634, 4.33139992746,
3.86272702658, 3.38639207628, 1.57773054352, 0.953410121155,
0.504987738102, 0.0995174088527, -0.0878099
])
data['r_in0'] = np.array(
[1.5375, 13.2620872, 17.35321975, 40.19535987, 63.0375])
data['Rhub0'] = 1.5375
data['r_pts0'] = np.array([
1.5375, 1.84056613137, 1.93905987557, 2.03755361977, 2.14220322299,
2.24069696719, 2.33919071139, 2.90419974, 3.04095863882, 3.13945238302,
5.637500205, 7.04226699698, 8.370800055, 11.8494282928, 13.8375,
14.7182548841, 15.887499795, 18.5537233505, 19.9875, 22.0564071286,
24.087500205, 26.16236509, 28.187499795, 32.2875, 33.5678634821,
36.387500205, 38.5725768593, 40.487499795, 40.6967540173,
40.7964789782, 40.8975, 42.6919644014, 44.5875, 52.787499795,
56.204199945, 58.937499795, 61.67080026, 63.0375
])
data['precurve0'] = np.zeros_like(
data['dx']) #np.linspace(0.0, 5.0, np.size(data['dx']))
data['bladeLength0'] = 61.5
npts = np.size(data['dx'])
prob = Problem()
prob.root = Group()
prob.root.add('comp', BladeDeflection(npts), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test(prob, tol=1e-4)
def test_CCBlade(bem_case):
data = {}
data['r'] = np.array([2.8667, 5.6000, 8.3333, 11.7500, 15.8500, 19.9500, 24.0500, 28.1500, 32.2500, 36.3500, 40.4500, 44.5500, 48.6500, 52.7500, 56.1667, 58.9000, 61.6333])
data['chord'] = np.array([3.542, 3.854, 4.167, 4.557, 4.652, 4.458, 4.249, 4.007, 3.748, 3.502, 3.256, 3.010, 2.764, 2.518, 2.313, 2.086, 1.419])
data['theta'] = np.array([13.308, 13.308, 13.308, 13.308, 11.480, 10.162, 9.011, 7.795, 6.544, 5.361, 4.188, 3.125, 2.319, 1.526, 0.863, 0.370, 0.106])
data['Rhub'] = 1.5
data['Rtip'] = 63.0
data['hubHt'] = 80.0
data['precone'] = 2.5
data['tilt'] = -5.0
data['yaw'] = 0.0
data['B'] = 3
data['rho'] = 1.225
data['mu'] = 1.81206e-5
data['shearExp'] = 0.2
data['nSector'] = 4
data['precurve'] = np.zeros(len(data['r']))
data['precurveTip'] = 0.0
# airfoils
basepath = os.path.join(os.path.join(os.path.dirname(os.path.realpath(__file__)),'..'), '5MW_AFFiles/')
# load all airfoils
airfoil_types = [0]*8
airfoil_types[0] = basepath + 'Cylinder1.dat'
airfoil_types[1] = basepath + 'Cylinder2.dat'
airfoil_types[2] = basepath + 'DU40_A17.dat'
airfoil_types[3] = basepath + 'DU35_A17.dat'
airfoil_types[4] = basepath + 'DU30_A17.dat'
airfoil_types[5] = basepath + 'DU25_A17.dat'
airfoil_types[6] = basepath + 'DU21_A17.dat'
airfoil_types[7] = basepath + 'NACA64_A17.dat'
# place at appropriate radial stations
af_idx = [0, 0, 1, 2, 3, 3, 4, 5, 5, 6, 6, 7, 7, 7, 7, 7, 7]
n = len(data['r'])
af = [0]*n
for i in range(n):
af[i] = airfoil_types[af_idx[i]]
data['airfoil_files'] = af
prob = Problem()
prob.root = Group()
if bem_case == 'power':
# run_case = 'power'
data['Uhub'] = np.array([3.0, 4.15789473684, 5.31578947368, 6.47368421053, 7.63157894737, 8.78947368421, 9.94736842105, 11.1052631579, 12.2631578947, 13.4210526316, 14.5789473684, 15.7368421053, 16.8947368421, 18.0526315789, 19.2105263158, 20.3684210526, 21.5263157895, 22.6842105263, 23.8421052632, 25.0])
data['Omega'] = np.array([3.43647024491, 4.76282718154, 6.08918411817, 7.41554105481, 8.74189799144, 10.0682549281, 11.3946118647, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0, 12.0])
data['pitch'] = np.array([0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
n2 = len(data['Uhub'])
prob.root.add('comp', CCBladePower(n, n2), promotes=['*'])
elif bem_case == 'loads':
# run_case = 'loads'
data['V_load'] = 12.0
data['Omega_load'] = 10.0
data['pitch_load'] = 0.0
data['azimuth_load'] = 180.0
n2 = 1
prob.root.add('comp', CCBladeLoads(n, n2), promotes=['*'])
prob = init_IndepVar_add(prob, data)
prob.root.comp.deriv_options['check_step_calc'] = 'relative'
prob.root.comp.deriv_options['check_form'] = 'central'
prob.setup(check=False)
prob = init_IndepVar_set(prob, data)
check_gradient_unit_test( prob, tol=1e-3, display=True)
