摘要: Starting from a variational formulation based on Hamilton’s Principle, the paper exploits finite element technique in time domain order to solve orbital dynamic problems characterised by constrained boundary value rather than initial problems. The solution is obtained assembling suitable number of elements inside interval interest, imposing desired constraints, and solving resultant set non-linear algebraic equations means Newton-Raphson method. In particular, this work general strategy applied periodic orbits determination. effectiveness approach finding unhomogeneous gravity field Moon assessed relevant examples, results are compared with those standard marching techniques as well analytical results.