This paper focuses on the optimization of two-dimensional transfer trajectories for solar sails using analytical approximation and modified analytical approximation methods. Instead of assuming circular orbits [29-31,37], the proposed approach allows for the analysis of solar sails launched from general heliocentric elliptical orbits. By employing linear perturbation theory, analytical propelled trajectories with fixed pitch angles are derived, where the perturbation coefficient serves as the reference lightness coefficient. To obtain the optimal transfer orbit between general heliocentric elliptical orbits, the analytical segments corresponding to different pitch angles are spliced together. A nonlinear optimization algorithm is then utilized to optimize the multi-segment trajectory, aiming to achieve the optimal performance index while satisfying boundary and intermediate constraints. Remarkably, this method requires the optimization of only the pitch angle for each segment, resulting in a small number of optimization variables and significantly improved efficiency in transfer orbit optimization. Additionally, the analytical approximation approach offers substantial time savings compared to numerical integration, particularly when repeated calculations are involved. Furthermore, the modified approximate trajectory exhibits high accuracy, with the optimal analytical transfer trajectory closely approximating the optimal actual trajectory acquired through the shooting method. This paper extensively validates and analyzes the precision of the proposed analytical transfer orbit through numerical verification. Additionally, when contrasted with actual transfer orbit, the suggested analytical approach can significantly decrease the computation time for multi-object challenges by a minimum of 80%. This makes the proposed method highly suitable for rapid design of optimal transfer orbits in multi-target mission scenarios for solar sail-based probes.