Start work on cl vs alpha curves and calculations

1 files changed, 14 insertions, 21 deletions

diff --git a/example_airfoil.py b/example_airfoil.py
index c20bdef..3dc319c 100644
--- a/example_airfoil.py
+++ b/example_airfoil.py
@@ -13,7 +13,7 @@ import resources.materials as mt
 import creator
 import evaluator.evaluator as evaluator
 import generator
-# from generator import
+import matplotlib.pyplot as plt
 
 import time
 start_time = time.time()
@@ -44,14 +44,13 @@ eval = evaluator.Evaluator("eval")
 ac = creator.base.Aircraft(eval, "ac")
 af = creator.wing.Airfoil(ac, 'af')
 af.add_naca(2412)
-# af.info_print(2)
-af.info_save()
 spar1 = creator.wing.Spar(af, 'spar1')
 spar2 = creator.wing.Spar(af, 'spar2', 0.57)
-# af.spar1.info_print(2)
-# af.spar2.info_print(2)
-spar1.info_save()
-spar2.info_save()
+# stringer = creator.wing.Stringer(af, 'stringer')
+# af.stringer.add_coord(af, [af.spar1, af.spar2], NOSE_TOP_STRINGERS,
+#                       TOP_STRINGERS, NOSE_BOTTOM_STRINGERS, BOTTOM_STRINGERS)
+# af.stringer.info_print(2)
+# af.stringer.info_save(SAVE_PATH, 'foo_name')
 
 ac2 = creator.base.Aircraft(eval, "ac2")
 af2 = creator.wing.Airfoil(ac2, 'af2')
@@ -62,6 +61,7 @@ spar4 = creator.wing.Spar(af2, 'spar4', 0.67)
 spar3.info_save()
 spar4.info_save()
 
+
 eval.tree_print()
 eval.tree_save()
 
@@ -73,20 +73,6 @@ print(eval.analysis())
 #     for spar in aircraft.wing.spars:
 #         print(spar, f"is made out of: {spar.material['name']}")
 
-# ac.info_print()
-
-# eval.info_print(2)
-
-# # Create stringer instance
-# af.stringer = wing.Stringer()
-# # Compute the stringer coordinates from their quantity in each zone
-# af.stringer.add_coord(af, [af.spar1, af.spar2], NOSE_TOP_STRINGERS, TOP_STRINGERS,
-# NOSE_BOTTOM_STRINGERS, BOTTOM_STRINGERS)
-# af.stringer.add_area(STRINGER_AREA)
-# af.stringer.add_webs(SKIN_THICKNESS)
-# af.stringer.info_print(2)
-# af.stringer.info_save(SAVE_PATH, 'foo_name')
-
 # Plot components with matplotlib
 # creator.wing.plot_geom(af, [af.spar1, af.spar2], None)
 
@@ -98,6 +84,13 @@ print(eval.analysis())
 # evaluator.plot_geom(eval)
 # evaluator.plot_lift(eval)
 
+# import resources.NACA_2412
+# cl = resources.NACA_2412.cl
+# alpha = resources.NACA_2412.alpha
+# plt.plot(alpha, cl)
+# plt.show()
+
+
 # Final execution time
 final_time = time.time() - start_time
 print(f"--- {round(final_time, 4)}s seconds ---")