public void plotClvsAlpha(String subfolderPath, String name) {
System.out.println("\n \n-----------------------------------------------------");
System.out.println("STARTING PLOT AIRFOIL CL vs ALPHA CURVE " + name);
System.out.println("-----------------------------------------------------");
MyArray alphaArray = new MyArray();
int _numberOfAlpha = 60;
double[] clArray = new double[_numberOfAlpha];
double[] alphaArrayDeg = new double[_numberOfAlpha];
Amount<Angle> alphaActualAmount;
Amount<Angle> alphaStart = Amount.valueOf(toRadians(-6.), SI.RADIAN);
Amount<Angle> alphaEnd = Amount.valueOf(toRadians(25.0), SI.RADIAN);
alphaArray.setDouble(
MyArrayUtils.linspace(
alphaStart.getEstimatedValue(), alphaEnd.getEstimatedValue(), _numberOfAlpha));
for (int i = 0; i < _numberOfAlpha; i++) {
alphaActualAmount = Amount.valueOf(alphaArray.get(i), SI.RADIAN);
clArray[i] = calculateClAtAlpha(alphaActualAmount.getEstimatedValue());
alphaArrayDeg[i] = alphaActualAmount.to(NonSI.DEGREE_ANGLE).getEstimatedValue();
}
MyChartToFileUtils.plotNoLegend(
alphaArrayDeg,
clArray,
null,
null,
null,
null,
"alpha",
"CL",
"deg",
"",
subfolderPath,
"CL vs alpha Airfoil" + name);
System.out.println("-----------------------------------------------------");
System.out.println("\n DONE");
System.out.println("-----------------------------------------------------");
}
// ------------------------------------------------------------------------------------------
// MAIN:
public static void main(String[] args)
throws HDF5LibraryException, NullPointerException, CmdLineException {
// --------------------------------------------------------------------------------------
// Arguments check and initial activities
if (args.length == 0) {
System.err.println("NO INPUT FILE GIVEN --> TERMINATING");
return;
}
System.out.println("------------------------------------------------------------");
System.out.println("PayloadRangeCalc_Test :: TURBOPROP");
System.out.println("------------------------------------------------------------");
System.out.println(" ");
PayloadRange_Test_TP_SA main = new PayloadRange_Test_TP_SA();
System.out.println("Input variable usage:");
main.theCmdLineParser.printUsage(System.out);
System.out.println("-----------------------------------------------------------");
main.theCmdLineParser.parseArgument(args);
String path = main.get_inputFile().getAbsolutePath();
JPADXmlReader reader = new JPADXmlReader(path);
System.out.println("-----------------------------------------------------------");
System.out.println("XML File Path : " + path);
System.out.println("-----------------------------------------------------------");
System.out.println("Initialize reading \n");
// --------------------------------------------------------------------------------------
// Input data
Amount<Angle> sweepLEEquivalent =
(Amount<Angle>) reader.getXMLAmountWithUnitByPath("//ATR72/SweepLE");
Amount<Mass> maxTakeOffMass =
(Amount<Mass>) reader.getXMLAmountWithUnitByPath("//ATR72/Maximum_take_off_mass");
Amount<Mass> operatingEmptyMass =
(Amount<Mass>) reader.getXMLAmountWithUnitByPath("//ATR72/Operating_empty_mass");
Amount<Mass> maxFuelMass =
(Amount<Mass>) reader.getXMLAmountWithUnitByPath("//ATR72/Maximum_fuel_mass");
Amount<Area> surface =
(Amount<Area>) reader.getXMLAmountWithUnitByPath("//ATR72/Planform_surface");
List<String> cd0_property = reader.getXMLPropertiesByPath("//ATR72/CD0");
double cd0 = Double.valueOf(cd0_property.get(0));
List<String> oswald_property = reader.getXMLPropertiesByPath("//ATR72/OswaldFactor");
double oswald = Double.valueOf(oswald_property.get(0));
List<String> ar_property = reader.getXMLPropertiesByPath("//ATR72/AspectRatio");
double ar = Double.valueOf(ar_property.get(0));
List<String> cl_property = reader.getXMLPropertiesByPath("//ATR72/Current_lift_coefficient");
double cl = Double.valueOf(cl_property.get(0));
List<String> tcMax_property = reader.getXMLPropertiesByPath("//ATR72/Mean_maximum_thickness");
double tcMax = Double.valueOf(tcMax_property.get(0));
List<String> altitude_property = reader.getXMLPropertiesByPath("//ATR72/Altitude");
double altitude = Double.valueOf(altitude_property.get(0));
List<String> eta_property = reader.getXMLPropertiesByPath("//ATR72/Propeller_efficiency");
double eta = Double.valueOf(eta_property.get(0));
List<String> currentMach_property = reader.getXMLPropertiesByPath("//ATR72/Mach_number");
double currentMach = Double.valueOf(currentMach_property.get(0));
List<String> byPassRatio_property = reader.getXMLPropertiesByPath("//ATR72/ByPassRatio");
double byPassRatio = Double.valueOf(byPassRatio_property.get(0));
List<String> nPassMax_property =
reader.getXMLPropertiesByPath("//ATR72/Maximum_number_of_passengers");
double nPassMax = Double.valueOf(nPassMax_property.get(0));
List<String> taperRatioEquivalent_property =
reader.getXMLPropertiesByPath("//ATR72/TaperRatio");
double taperRatioEquivalent = Double.valueOf(taperRatioEquivalent_property.get(0));
Amount<Angle> sweepHalfChordEquivalent =
LSGeometryCalc.calculateSweep(
ar,
taperRatioEquivalent,
sweepLEEquivalent.to(SI.RADIAN).getEstimatedValue(),
0.5,
0.0);
// --------------------------------------------------------------------------------------
// the next step is necessary to calculate CL and Speed at point E (or A) of the parabolic polar
ACAerodynamicsManager analysis = new ACAerodynamicsManager();
analysis.calculateDragPolarPoints(
ar,
oswald,
cd0,
AtmosphereCalc.getDensity(altitude),
maxTakeOffMass.getEstimatedValue() * 9.81,
surface.getEstimatedValue());
// --------------------------------------------------------------------------------------
// Assign all default folders
MyConfiguration.initWorkingDirectoryTree();
// --------------------------------------------------------------------------------------
// Setup database(s)
String databaseFolderPath = MyConfiguration.getDir(FoldersEnum.DATABASE_DIR);
String fuelFractionDatabaseFileName = "FuelFractions.h5";
FuelFractionDatabaseReader fuelFractionReader =
new FuelFractionDatabaseReader(databaseFolderPath, fuelFractionDatabaseFileName);
// --------------------------------------------------------------------------------------
// Creating calculator object
PayloadRangeCalc test =
new PayloadRangeCalc(
maxTakeOffMass,
operatingEmptyMass,
maxFuelMass,
sweepHalfChordEquivalent,
nPassMax,
cl,
tcMax,
AirplaneType.TURBOPROP_REGIONAL,
EngineTypeEnum.TURBOPROP,
AirfoilTypeEnum.CONVENTIONAL,
analysis,
fuelFractionReader);
// -----------------------CRITICAL MACH NUMBER CHECK----------------------------
boolean check = test.checkCriticalMach(currentMach);
if (check)
System.out.println(
"\n\n-----------------------------------------------------------"
+ "\nCurrent Mach is lower then critical Mach number."
+ "\nCurrent Mach = "
+ currentMach
+ "\nCritical Mach = "
+ test.getCriticalMach()
+ "\n\n\t CHECK PASSED --> PROCEDING TO CALCULATION "
+ "\n\n"
+ "-----------------------------------------------------------");
else {
System.err.println(
"\n\n-----------------------------------------------------------"
+ "\nCurrent Mach is bigger then critical Mach number."
+ "\nCurrent Mach = "
+ currentMach
+ "\nCritical Mach = "
+ test.getCriticalMach()
+ "\n\n\t CHECK NOT PASSED --> TERMINATING "
+ "\n\n"
+ "-----------------------------------------------------------");
return;
}
// ------------------------MTOM PARAMETERIZATION---------------------------------
test.createPayloadRangeMatrices(
sweepHalfChordEquivalent,
surface,
cd0,
oswald,
cl,
ar,
tcMax,
byPassRatio,
eta,
altitude,
currentMach,
false);
// ------------------------------PLOTTING-----------------------------------------------
// MTOM parameterization:
test.createPayloadRangeCharts_MaxTakeOffMass(test.getRangeMatrix(), test.getPayloadMatrix());
}