-1

I'm trying to calculate the angle of attack (AOA) using aircraft some data. From what I understand the AOA is just the angle between the aircraft's forward unit vector (out the nose vector) and the x-component of the velocity vector (vx), as long as both vectors are in the same coordinate frame.

My current plan is to find the aircraft's forward unit vector, and then find the angle between that forward vector and the x-component of the aircraft's velocity vector (vx) using the dot product.

All of my data is in the ECEF frame and consists of the position (x,y,z), velocity components (vx,vy,vz), and the orientation components (psi,theta,phi).

There are two stack exchange answers that I think are on the right track, but they are not using the ECEF frame.

In the first answer, https://stackoverflow.com/a/1568687/12131175 the unit vectors of the aircraft can be found using the rotation matrices, but I think they are using the NED coordinate frame and not ECEF.

In the second answer, https://aviation.stackexchange.com/a/67213/52074 the AOA is found using the arctan function, but vectors are in the aircrafts body frame.

So, my questions are:

  1. Am i on the right track on how to calculate the AOA? Is this the best way?
  2. What are the correct rotation matrices needed to solve for the unit vectors in the ECEF frame?

It seems like there's a lot of different sets of rotation matrices out there for the different orders of rotations (x-y-z vs z-y-z and extrinsic vs intrinsic). Some use different greek letters for the same rotation (i.e. rotation about z-axis, some use φ and some use ψ). I'm not certain which rotation matrices are correct for describing an aircrafts orientation in the ECEF frame.

Thanks in advance!

  • 1
    AOA is actually to the wing chord line, not “out the nose”, so make sure you factor in any angle of incidence built into the airframe... – Michael Hall Aug 20 '20 at 16:59
  • 1
    What's the ECEF frame? – JZYL Aug 20 '20 at 17:33
  • @MichaelHall Thanks. – Marshin8726 Aug 20 '20 at 17:41
  • @JZYL Earth Centered Earth Fixed frame. Also known as ECR. https://en.wikipedia.org/wiki/ECEF – Marshin8726 Aug 20 '20 at 17:42
  • 2
    No, you can't find AOA simply by looking at inertial velocity. You need to know winds too (expressed in whichever frame you're doing your calculations in). – JZYL Aug 20 '20 at 17:43
  • @JZYL Understood. I guess that's why the AOA is usually measured using an onboard sensor? I just wanted a theoretical value by assuming the wind vector was in the opposite direction of the velocity vector. I suppose this may be evolving into more of a math/physics question if it's not the actual AOA that I'm calculating since I don't have any wind data. – Marshin8726 Aug 20 '20 at 18:02
  • 1
    I think you misunderstand the definition of AoA. It is defined between mean chord (not zero angle with the fuselage) and relative wind (not necessarily constant in time). See skybrary for more information. – Manu H Aug 20 '20 at 18:07
  • @ManuH I think the data is actually for a missile and not a plane, so in that case could I assume the chord line is on the same axis as the out the nose vector? There may be fins, but I assume those are also aligned with the missile's center axis. – Marshin8726 Aug 20 '20 at 18:16
  • It doesn't make sense to do this in ECEF. Convert the speed vector to the body reference frame and calculate the angle between the body x axis and the projection of the speed vector on the xz plane. Add the angle of incidence and you're done. Assuming no wind. – DeltaLima Aug 20 '20 at 21:24
  • @DeltaLima Do you mean to convert the speed vector to the NED frame? Exactly which frame should I convert it to? – Marshin8726 Aug 21 '20 at 15:17
  • You should convert everything to the body frame – DeltaLima Aug 21 '20 at 15:20
  • @DeltaLima Thanks for the info. I'm still confused, doesn't the body frame just move with the object? So, wouldn't it have to be relative to some fixed frame? Do you have any links on how to do what you're referring to? Thanks! – Marshin8726 Aug 21 '20 at 15:43

1 Answers1

0

Regarding your first question

  • The angle of attack is the angle between the chord line of the airfoil and the relative airflow, that's first and foremost.

    There is a longitudinal datum that could be what you are referring to as the "aircraft's forward unit vector" but that longitudinal datum shouldn't be measured directly to the velocity vector in your case because that's not what the air considerably is passing through, the air passes through the airfoil which has the chord line in it, the angle between the chord line and the longitudinal datum is called the angle of incidence and it's a fixed angle.

  • The velocity vector could be the relative airflow which is a combination of more than one direction of wind.

Regarding your second question

  • I think the concentration should be towards the lateral axis of the airfoil which that airflow is angled toward and rotates on to create a pitch angle whether it is a positive or a negative angle of attack.

This video about the propeller's angle of attack would help you understand the combination of the velocity vector:

This would explain your second question further (Link is http only 'not secure') : http://www.prepar3d.com/SDKv2/LearningCenter/simobjects/weapon_systems/weapon_attachments.html