Mad Scientist,
i agree that in a theoric way, standard motion equation shoudn't be used in airplane axis system. But in a real time numeric-model, with a rate of min. 30 updates/sec ,i feel that the reference system's motion is really slight and doesn't really matter. for me, a good solution is :
- compute aerodynamic forces in aerdynamics reference (x on airspeed vector, y in airplane xy plane, z direct)
- compute inertial forces in the local (airplane) reference x longitudinal axis, y right wing...
- compute weights in terrestrial reference. (xy horizontal, z upward)
when it is done, transform all in airplane local reference and apply Newton Law.
I'have already wrote some code with that and plane's behaviour was correct.
However, i hadn't efficient means to quantify the error caused by the reference problem that you had pointed out.