Computer-assisted planning methods have increasingly contributed to preoperative ablation planning; however, these methods cannot automatically obtain the final optimal solution within a short time and are rarely validated in practice, greatly limiting their clinical applicability. We aimed to propose a full-automatic multi-stage ablation and needle trajectory planning method for CT-guided percutaneous liver ablation to attain the final optimal plans under multiple clinical constraints rapidly. Our proposed method integrates the ablation zone planning fulfilling complete tumor coverage and critical structure avoidance while reaching a trade-off between ablation number and healthy tissue damage, and needle trajectory planning under multiple clinical constraints. Our needle trajectory planning determines feasible skin entry regions based on hard constraints, where the multi-objective optimization (MOO) considering soft constraints is performed using the Pareto Optimality and Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) methods for the final optimal solution. The performance of our proposed method was evaluated on 30 tumors of various characteristics from 23 patients and clinically validated in five clinical cases. Our proposed method achieved 99.8% treatment zone coverage and 40.5% ablation efficiency without involving critical structures, and completely satisfied multiple clinical constraints in all needle trajectory planning results. The average planning time was 23.6 s for tumors of different sizes. All the plans were considered clinically acceptable by the doctors' evaluation. Our method achieved complete tumor coverage without complications in clinical case validation. Our proposed planning method can generate a final optimal plan satisfying multiple clinical constraints within a short time, potentially facilitating preoperative planning for hepatic tumor ablation.