A question about SkyHogg model in Aerospace Blockset
13 次查看(过去 30 天)
显示 更早的评论
in the example model:
asbSkyHogg
in this block (asbSkyHogg/Vehicle System Model/Vehicle/3DOF to 6DOF/3DOF to 6DOF/calc alpha_dot) which is used to calcualate the rate of change of of angle of attack ()
The rate of change of of angle of attack is calculated based on the equation:
So, the acceleration (input #3) in body axes should be transformed to wind axes using the DCM Body-to-Wind as shown in the picture... but why the transpose block is added here???!! In this case, it's actually transformed to body axes again!!
0 个评论
回答(1 个)
Paul
2023-1-28
编辑:Paul
2023-1-28
Hi Islam
I agree with your equation for alphadot, but not the notation. By standard convention (at least what I'm familiary with), u is the component of the velocity vector along the x-axis of the body frame, not the wind frame. By definition, the component of the velocity vector alone the x-axis of the wind frame is the velocity vector. So I would write the equation as
Here, V(t) is the magnitude of the velocity vector and A_zw(t) is the component of the inertial accelaration along the z-axis of the wind axes.
If I followed the model correctly, the signal called Vw should just be [V(t) 0 0].
I agree with you that the the transpose block is an error.
There is also another error. The signal called Ab is NOT the inertial acceleration vector resolved in the body axes. According to the documentation Ab = [udot 0 wdot], which is not what we want. So even if you delete the transpose block, the alpha_dot signal will still be incorrect. Probably not by much, but incorrect. To correct this error you'd have to get the optional output from the 3DOF block, which is the inertial acceleration in inertial (earth) coordinates, and then resolve that into body coordinates and use the result as the input to the "calc alpha_dot" block, instead of Ab.
4 个评论
Paul
2023-2-3
I should have been clearer. When I said Vi is the inertial velocity vector I was assuming that the absSkyHogg model uses a flat earth model, so the velocity relative to an inertial frame would be the same as the velocity relative to coordinate frame "fixed" to the atmosphere.
In this equation: dalpha/dt = q + Z/(mU), the uppercase U in the denominator is the vehicle speed (relative to the atmosphere). That's not same as lowercase u in the other equation for dalpha/dt.
My understanding is that the x-axis of the wind frame is defined by the velocity vector. Therefore, by definition the y- and z-components have to be zero in the wind frame.
Here's one reference: flight dynamics
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Get Started with Aerospace Blockset 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!