Perform analysis on the expected battery consumption based on hypothetical flight path trajectory.
Use the following input parameter:
# Trajectory Study
# AMAT 2024-2025 AIAA DBF
# Outline
# Find the energy spent
# E = integral(Power over time)
# Cruise --> Thrust_x_dir = Drag
# From the Drag, compute required Thrust
# From Thrust, compute Power
# Basic actuator disk/table Look up
# Assumptions for Cruise Flight:
# Cruise speed = 80 ft/s
# Drag from Drag Polar (provided by Esther)
# CL = 2pi(alpha - alpha_0)
# Re = ~5e5
# L/D = ~7
# Thrust = 2lbf (take off = max thrust)
# Thrust --> Power (Table look up)
# Energy = Int(power) = Power*(Distance/Velocity) for cruise
# Intregate power over time
# To find time -->
# find acceleration
# find velocity change
# find position change
# Assumptions for Turns:
# Cruise Speed = 80 ft/s
# Lift = load factor (n) * weight (W)
# For 2g turn: n = 2 (60deg turn)
# For 1.4g turn: n = sqrt(2) (45deg turn)
# L = W/(cos(phi))
# Get Cl from L
# usually CL = ~0.5
# Get Cd from Cl
# Get Drag from Cd
# AR = 6, e = 0.7
# CD = CD0 + CL^2/(pi*e*AR)
# Cl/(7) = CD0 + CL^2/(pi*e*AR)
# CD0 --> Through OpenVSP
# CD0 = ~0.052
# get thrust from drag
# get power from thrust
# get distance from centripetal force
# find Acceleration from bank angle (a = sin(phi))
# from there you can find the radius (a = v^2/R = g*tan(phi))
# r = v^2/(g*tan(phi)) --> pi*r = distance travel in half turn (x4)
# make this turn as energy efficient as possible
# Tasks:
# write the following functions:
# Drag polar function ( given alpha, return CL and CD )
# A solver for the above function, given CL, return alpha and CD
# Output: get a converged alpha based on required CL
# From there, get a CD estimate
# A function for each segment (cruise section, half turn) given velocity
# Find estimated total energy required complete one lap of the course (ideal conditions)