Download An Integrated, Modular Simulation System for Education and

% |
sAcB -sAsB
cA | Fzb
% Test data
FxWind = 100;
FyWind = 30;
FzWind = 300;
u=100;
v=50;
w=50;
%***** Trig functions first *****
alpha = atan(w/u);
beta = asin(v/AirSpeed);
%alpha = asin(w/AlphaSpeed)*57.3
%alpha = acos(u/AlphaSpeed)*57.3
%beta = asin(v/AirSpeed)*57.3
%***** just lift and drag *****
FxBody = FxWind*cos(alpha)*cos(beta)-FzWind*sin(alpha);
FyBody = FxWind*sin(beta);
Fzbody = FxWind*sin(alpha)*cos(beta)+FzWind*cos(alpha);
%***** with side force *****
FxBody = FxWind*cos(alpha)*cos(beta)-FyWind*cos(alpha)*sin(beta)-FzWind*sin(alpha)
FyBody = FxWind*sin(beta)+FyWind*cos(beta)
Fzbody = FxWind*sin(alpha)*cos(beta)-FyWind*sin(beta)*sin(alpha)+FzWind*cos(alpha)
%***** Without trig functions *****
% computers are VERY SLOW at trig functions and fast at mult ok at divide
% sA = w/AlphaSpeed
% cA = u/Alphaspeed
% sB = v/Airspeed
% cB = AlphaSpeed/Airspeed
vSquared = u*u+w*w;
154