This paper deals with the reentry trajectory optimization problem for lunar return with consideration of entry vehicle’s fore-body shape. Three performance objectives are applied in this work: cross range, peak heat flux and total heat load. Aerothermodynamic models are based on modified Newtonian impact theory and semi-empirical correlations for convective and radiative stagnation-point heat transfer. A population based evolutionary algorithm has been executed to optimize the multidisciplinary problem. At last the numerical example showed the Pareto frontiers for spherical segment and sphere cone respectively, one of optimal trajectory designs selected from the Pareto frontiers are showed in this paper. The mission requirements are satisfied through the aerothermodynamic balance.