def getAppliedForce(self, rocketState, time, environment, CG):
Aref = self.rocket.Aref
Mach = AeroParameters.getMachNumber(rocketState, environment)
# Normal Force --------------------------------------------------------------------------------
# TODO: Account for rate of pitch/yaw rotation in AOA calculation? Or do separate Pitch/Yaw damping moments?
AOA = AeroParameters.getTotalAOA(rocketState, environment)
normalForceCoefficient = AeroFunctions.Barrowman_GetCN(
AOA, Aref, 0, self.baseArea)
# Drag Force ---------------------------------------------------------------------------------------------
#### Skin Friction ####
skinFrictionDragCoefficient, rollDampingMoment = AeroFunctions.getCylindricalSkinFrictionDragCoefficientAndRollDampingMoment(rocketState, environment, self.length, Mach, \
self.surfaceRoughness, self.wettedArea, Aref, self.rocket.finenessRatio, self.rocket.fullyTurbulentBL)
#### Pressure ####
pressureDragCoefficient = self._getCd_Pressure(Mach)
totalDragCoefficient = pressureDragCoefficient + skinFrictionDragCoefficient
# Damping moments --------------------------------------------------------------------------------------
# Roll damping due to skin friction
dampingMoments = rollDampingMoment
# Combine forces and return total -------------------------------------------------------------------------------------
return AeroFunctions.forceFromCdCN(rocketState,
environment,
totalDragCoefficient,
normalForceCoefficient,
self.CPLocation,
Aref,
moment=dampingMoments)
def test_BarrowmanGetCN(self):
AOA = 0.1 #rad
Aref = 1
xArea_top = 1
xArea_bottom = 2
expected = 0.199666833
self.assertAlmostEqual(
AeroFunctions.Barrowman_GetCN(AOA, Aref, xArea_top, xArea_bottom),
expected)