mirror of
https://github.com/dcs-retribution/dcs-retribution.git
synced 2025-11-10 15:41:24 +00:00
88b9ed29ba
Previously we were trying to make every potential flight plan look just like a strike mission's flight plan. This led to a lot of special case behavior in several places that was causing us to misplan TOTs. I've reorganized this such that there's now an explicit `FlightPlan` class, and any specialized behavior is handled by the subclasses. I've also taken the opportunity to alter the behavior of CAS and front-line CAP missions. These no longer involve the usual formation waypoints. Instead the CAP will aim to be on station at the time that the CAS mission reaches its ingress point, and leave at its egress time. Both flights fly directly to the point with a start time configured for a rendezvous. It might be worth adding hold points back to every flight plan just to ensure that non-formation flights don't end up with a very low speed enroute to the target if they perform ground ops quicker than expected.
189 lines
7.1 KiB
Python
189 lines
7.1 KiB
Python
from __future__ import annotations
|
|
|
|
import logging
|
|
import math
|
|
from datetime import timedelta
|
|
from typing import Optional, TYPE_CHECKING
|
|
|
|
from dcs.mapping import Point
|
|
from dcs.unittype import FlyingType
|
|
|
|
from game.utils import meter_to_nm
|
|
from gen.flights.flight import Flight
|
|
|
|
if TYPE_CHECKING:
|
|
from gen.ato import Package
|
|
|
|
|
|
class GroundSpeed:
|
|
|
|
@classmethod
|
|
def for_flight(cls, flight: Flight, altitude: int) -> int:
|
|
if not issubclass(flight.unit_type, FlyingType):
|
|
raise TypeError("Flight has non-flying unit")
|
|
|
|
# TODO: Expose both a cruise speed and target speed.
|
|
# The cruise speed can be used for ascent, hold, join, and RTB to save
|
|
# on fuel, but mission speed will be fast enough to keep the flight
|
|
# safer.
|
|
|
|
c_sound_sea_level = 661.5
|
|
|
|
# DCS's max speed is in kph at 0 MSL. Convert to knots.
|
|
max_speed = flight.unit_type.max_speed * 0.539957
|
|
if max_speed > c_sound_sea_level:
|
|
# Aircraft is supersonic. Limit to mach 0.8 to conserve fuel and
|
|
# account for heavily loaded jets.
|
|
return int(cls.from_mach(0.8, altitude))
|
|
|
|
# For subsonic aircraft, assume the aircraft can reasonably perform at
|
|
# 80% of its maximum, and that it can maintain the same mach at altitude
|
|
# as it can at sea level. This probably isn't great assumption, but
|
|
# might. be sufficient given the wiggle room. We can come up with
|
|
# another heuristic if needed.
|
|
mach = max_speed * 0.8 / c_sound_sea_level
|
|
return int(cls.from_mach(mach, altitude)) # knots
|
|
|
|
@staticmethod
|
|
def from_mach(mach: float, altitude: int) -> float:
|
|
"""Returns the ground speed in knots for the given mach and altitude.
|
|
|
|
Args:
|
|
mach: The mach number to convert to ground speed.
|
|
altitude: The altitude in feet.
|
|
|
|
Returns:
|
|
The ground speed corresponding to the given altitude and mach number
|
|
in knots.
|
|
"""
|
|
# https://www.grc.nasa.gov/WWW/K-12/airplane/atmos.html
|
|
if altitude <= 36152:
|
|
temperature_f = 59 - 0.00356 * altitude
|
|
else:
|
|
# There's another formula for altitudes over 82k feet, but we better
|
|
# not be planning waypoints that high...
|
|
temperature_f = -70
|
|
|
|
temperature_k = (temperature_f + 459.67) * (5 / 9)
|
|
|
|
# https://www.engineeringtoolbox.com/specific-heat-ratio-d_602.html
|
|
# Dependent on temperature, but varies very little (+/-0.001)
|
|
# between -40F and 180F.
|
|
heat_capacity_ratio = 1.4
|
|
|
|
# https://www.grc.nasa.gov/WWW/K-12/airplane/sound.html
|
|
gas_constant = 286 # m^2/s^2/K
|
|
c_sound = math.sqrt(heat_capacity_ratio * gas_constant * temperature_k)
|
|
# c_sound is in m/s, convert to knots.
|
|
return (c_sound * 1.944) * mach
|
|
|
|
|
|
class TravelTime:
|
|
@staticmethod
|
|
def between_points(a: Point, b: Point, speed: float) -> timedelta:
|
|
error_factor = 1.1
|
|
distance = meter_to_nm(a.distance_to_point(b))
|
|
return timedelta(hours=distance / speed * error_factor)
|
|
|
|
|
|
class TotEstimator:
|
|
# An extra five minutes given as wiggle room. Expected to be spent at the
|
|
# hold point performing any last minute configuration.
|
|
HOLD_TIME = timedelta(minutes=5)
|
|
|
|
def __init__(self, package: Package) -> None:
|
|
self.package = package
|
|
|
|
def mission_start_time(self, flight: Flight) -> timedelta:
|
|
takeoff_time = self.takeoff_time_for_flight(flight)
|
|
startup_time = self.estimate_startup(flight)
|
|
ground_ops_time = self.estimate_ground_ops(flight)
|
|
return takeoff_time - startup_time - ground_ops_time
|
|
|
|
def takeoff_time_for_flight(self, flight: Flight) -> timedelta:
|
|
travel_time = self.travel_time_to_rendezvous_or_target(flight)
|
|
if travel_time is None:
|
|
logging.warning("Found no join point or patrol point. Cannot "
|
|
f"estimate takeoff time takeoff time for {flight}")
|
|
# Takeoff immediately.
|
|
return timedelta()
|
|
|
|
from gen.flights.flightplan import FormationFlightPlan
|
|
if isinstance(flight.flight_plan, FormationFlightPlan):
|
|
tot = flight.flight_plan.tot_for_waypoint(
|
|
flight.flight_plan.join)
|
|
if tot is None:
|
|
logging.warning(
|
|
"Could not determine the TOT of the join point. Takeoff "
|
|
f"time for {flight} will be immediate.")
|
|
return timedelta()
|
|
else:
|
|
tot = self.package.time_over_target
|
|
return tot - travel_time - self.HOLD_TIME
|
|
|
|
def earliest_tot(self) -> timedelta:
|
|
return max((
|
|
self.earliest_tot_for_flight(f) for f in self.package.flights
|
|
)) + self.HOLD_TIME
|
|
|
|
def earliest_tot_for_flight(self, flight: Flight) -> timedelta:
|
|
"""Estimate fastest time from mission start to the target position.
|
|
|
|
For BARCAP flights, this is time to race track start. This ensures that
|
|
they are on station at the same time any other package members reach
|
|
their ingress point.
|
|
|
|
For other mission types this is the time to the mission target.
|
|
|
|
Args:
|
|
flight: The flight to get the earliest TOT time for.
|
|
|
|
Returns:
|
|
The earliest possible TOT for the given flight in seconds. Returns 0
|
|
if an ingress point cannot be found.
|
|
"""
|
|
time_to_target = self.travel_time_to_target(flight)
|
|
if time_to_target is None:
|
|
logging.warning(f"Cannot estimate TOT for {flight}")
|
|
# Return 0 so this flight's travel time does not affect the rest
|
|
# of the package.
|
|
return timedelta()
|
|
startup = self.estimate_startup(flight)
|
|
ground_ops = self.estimate_ground_ops(flight)
|
|
return startup + ground_ops + time_to_target
|
|
|
|
@staticmethod
|
|
def estimate_startup(flight: Flight) -> timedelta:
|
|
if flight.start_type == "Cold":
|
|
if flight.client_count:
|
|
return timedelta(minutes=10)
|
|
else:
|
|
# The AI doesn't seem to have a real startup procedure.
|
|
return timedelta(minutes=2)
|
|
return timedelta()
|
|
|
|
@staticmethod
|
|
def estimate_ground_ops(flight: Flight) -> timedelta:
|
|
if flight.start_type in ("Runway", "In Flight"):
|
|
return timedelta()
|
|
if flight.from_cp.is_fleet:
|
|
return timedelta(minutes=2)
|
|
else:
|
|
return timedelta(minutes=5)
|
|
|
|
@staticmethod
|
|
def travel_time_to_target(flight: Flight) -> Optional[timedelta]:
|
|
if flight.flight_plan is None:
|
|
return None
|
|
return flight.flight_plan.travel_time_to_target
|
|
|
|
@staticmethod
|
|
def travel_time_to_rendezvous_or_target(
|
|
flight: Flight) -> Optional[timedelta]:
|
|
if flight.flight_plan is None:
|
|
return None
|
|
from gen.flights.flightplan import FormationFlightPlan
|
|
if isinstance(flight.flight_plan, FormationFlightPlan):
|
|
return flight.flight_plan.travel_time_to_rendezvous
|
|
return flight.flight_plan.travel_time_to_target
|