position change in ECEF EoM

Hi
I might have some rust on my mechanics, but I couldn't understand why both the ECEF and the LLA position changes with time under a zero body forces and moments.
It is easy to test, an ECEF quaternion EoM block with zero inputs and arbitrary initial position (not at a pole) and zero initial body translational and rotational velocities and a WGS84 planet.
I know that the position change is due to the transport term (specifically the centripetal acceleration) emerging from solving eom in a rotating frame. However, these transport terms will affect the inertial position (ECI position) not the ECEF and LLA positions. Thus a stationary point on a rotating planet will have a constant ECEF and LLA position, but its ECI position will be changing with time.
Would anyone kindly explain what am I missing? or is the block doing something it shouldn't?
Kind regards,
Mohamad

답변 (1개)

Saif Awwad
Saif Awwad 2015년 5월 8일

0 개 추천

I faced a somehow similar problem some time ago, which I could not answer until today. Please let me know if you find an answer.
Thanks.

댓글 수: 4

MTA
MTA 2015년 5월 11일
OK Saif, I'll let you know if I had a progress, you too please.
Saif Awwad
Saif Awwad 2015년 5월 13일
Hi Mohamad. Thanks for the reply.
Frankly, I doubt I will have any progress. This is an old problem on which I got stuck some time ago, so I decided to leave it aside until a brighter idea comes up. And now you helped me bring it back to the surface again.
I am not sure whether your time allows, but I am kind of hoping to get some help on a part of my work. Is it okay if I send you a file or two to check and confirm?
MTA
MTA 2015년 5월 13일
Hi Saif
Sure, you can personal message it if it is not contributing to the thread directly.
Saif Awwad
Saif Awwad 2015년 5월 16일
With a body of zero forces and moments, is the change in both of ECEF and LLA positions analogical? As you said, only ECI should be changing with time. So I am doubting some sort of a bug in the code converting ECEF to LLA or vise versa.
I will be messaging about further topics outside the scope of this thread.

댓글을 달려면 로그인하십시오.

카테고리

질문:

MTA
2015년 4월 30일

댓글:

2015년 5월 16일

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by