public static float AerodynamicCentre(this VerticalStabilizer ld)
{
var lambda = ld.Sweep * Mathf.Deg2Rad;
var cBar = (0.066667f) * ld.RootChord * (1 + lambda + Mathf.Pow(lambda, 2)) / (1 + lambda);
var YBar = (ld.Height / 3) * ((1 + 2 * lambda) / (1 + lambda));
var centre = Mathf.Tan(lambda) * YBar + 0.25f * cBar;
return(centre + ld.GetXPosition());
}
private static float CalculateVSCD_0(VerticalStabilizer vs, float S_ref, float Re, float C_f)
{
// Calculate form factor
var gamma_quarter = vs.CalculateSweep(0.25f) * Mathf.Deg2Rad;
var S_wet = vs.Area / 100f;
var xc = 0.25f;
var tc = vs.ThicknessChordRatio;
var FF = (1 + (0.6f / xc * tc) + 100 * Mathf.Pow(tc, 4)) * (1.34f * Mathf.Pow(Constants.Metrics.DesignMachNumber, 0.18f) * Mathf.Pow(Mathf.Cos(gamma_quarter), 0.28f));
return((C_f * FF * S_wet) / S_ref);
}