An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincaré sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L1, L2, and L3 points as well as Lissajous orbits for L1 and L2 points in the Earth-Moon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.