Sign in to follow this  

Calculating forces, speed,acceleration of a plane

Recommended Posts

I am trying to make the plane fly by calculating forces of thrust gravity drag and lift (but now i am only using now two forces: gravity and thrust)

nvm because i dont have a lift force because of Cl = 2 * pi * attack angle = 0.0ff

the stupidiest thing now i see that: my plane always goes down towards gravity vector even when thrust force is 2 times greater than gravity force and acts in the opposite direction

second issue is that: i cant calculate proper normalized vector (normalized gravity force vector)
in example: i have plane pos: x=50 000 m y = 100 m z = 30 000 m and the normalized vector is something like this: x = 0.7; y = 1.5, z= 0.4; and i donkt know why the y is fck-ing 1.5 here

i always get this effect (lol clouds are 320 km away form the surface - and we are in south america )

i use these equations:
s= v0*t+a*t*t*0.5

i start the simulation with this:

MODEL->t = GetTickCount();
MODEL->pos = position;
MODEL->old = position;

void __fastcall TVehicle::ProcessFrame()
int p = GetTickCount();
float x = abs(float(p-t))/1000.0f;
if (x <=0 ) return;
t3dpoint kk;
kk = YPRangle.rf; //thrust force orientation

float Front_force = Throttle_to_force( Throttle ); //thrust force val
float Q = eGForce * mass; // gravity force val

t3dpoint TForce_vec;
t3dpoint GForce_vec;

TForce_vec = YPRangle.rf; //thrust force direction (movement)

t3dpoint c; c = vectorAB(pos, zero);//center of earth
GForce_vec = Normalize( c ); //g force direction

float drag_force = 0.350*pressure*(V*V/2.0f);

GForce_vec = vector_multiple(GForce_vec, Q);
TForce_vec = vector_multiple(TForce_vec, Front_force);

t3dpoint ALL_FORCES = vectors_add(GForce_vec, TForce_vec);//add forces
//Cl = 2 * pi * pitch angle (in radians)
//L = .5 * Cl * air_density * V^2 * A
t3dpoint VA; //acceleration vector

VA = ALL_FORCES / mass; // a = F / m

float dk = ( x*x ) / 2.0f;

VA = VA * dk; //(a * t^2 )/ 2

float h = V*x;
VA = VA + h; //V0*t + (a * t^2 )/ 2

//s = V0*t + (a * t^2 )/ 2

pos = vectors_add(pos,VA );
t = GetTickCount();
V = n3ddistance(pos,old) / x; // s = V / t so calculate velocity
old = pos;

what am i doing wrong ?

now wypadkowa is ALL_FORCES :] no its not irricht its mine engine

Share this post

Link to post
Share on other sites
I would certainly start with getting rid of GetTickCount() in this kind of use. Resolution of this function is somewhere between 10 and 20 ms, which absolutely isn't enough even for applications running at mere 60 FPS. Try QueryPerformanceCounter() instead.

I didn't have much time to dive into your code, I just noticed that it's quite chaotic and confusing. You should give better names to your variables. For example your "x" is usually called "dt", or "delta_time".
X should be strictly reserved for coordinates.

And short names like VA, h, V, dk tell nothing about their meaning, especially when for example "h" doesn't mean "height".

Share this post

Link to post
Share on other sites
[quote name='___' timestamp='1329954678' post='4915706']
and the GetTickCount() has resolution such as 1 ms not 10 ms
Nope, it doesn't.
It seems to have a 1 ms resolution, because it returns a DWORD value in miliseconds. But if you test it better, you'll notice that when for example you are calling the function repeatedly, it may return something like

In fact, you won't get two succeeding values with difference of 1 ms, it will be returning the same value and after some time the value will suddenly increase in a step.

This limited resolution is mentioned also in the official documentation

Share this post

Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!

Register a new account

Sign in

Already have an account? Sign in here.

Sign In Now

Sign in to follow this