Thank you very much for your help. I remember reading your post before making this thread. One thing I haven't understood well is if bodies in KSP generate lift too (like cylindrcal bodies of rockets).
I think there's a bit of confusion here, but you pointed out I, indeed, forgot the squared term.
I derived that formula from the one you wrote in your thread:
In a 2D plane you could rewrite that formula as (by looking at the image, assuming a dragcube with only +x and +y faces)
D = ½ ρ ( V² cos2(α) Ax Cd,x + V² sin2(α) Ay Cd,y ) = ½ ρ V² (cos2(α) Ax Cd,x + sin2(α) Ay Cd,y)
Regarding real life, in my aerodynamics course I never heard about a sine-squared tendency of Cd with respect to angle of attack. For airfoils the behaviour is pseudo-polynomial, similar to x2 . For instance for the NACA 0015:
Other arifoils may have a low near-constant drag region (called laminar bucket) followed by a steep increase of drag outside of it.
While induced drag (which afaik is not modeled by stock aerodynamics, but Ferram does) is proportional to the square of the generated lift.