In Phase I, the ZONA unified hypersonic/supersonic aerodynamic method ZONA7U in the aerothermoelastic software development for TPS/TAV design/analysis, was proven a successful tool through feasibility cases studied in the CKEM body, blunt cones and X-34 wing-body. Further upgrade of ZONA7U to a high-fidelity panel method ZONAIR is proposed in order to enhance the universality of the aerodynamic capability with: 1) unified AIC for hypersonic/supersonic/subsonic Mach numbers; 2) exact paneling to model OML surface; 3) high AoA flow; 4) corrected location/strength of impinged shocks; 5)two-body interference aerodynamics. Other proposed improvements of the aerothermoelastic software include: i) incorporation of AEROHEAT to replace SHABP with ZONAIR and ZSTREAM; ii) optimized TPS sizing procedure using complex variable differentiation sensitivity with MINIVER/EXIT; iii) applying temperature mapping onto structural FEM. The structural optimizer ASTROS* is used throughout for sizing TAV/TPS while maintaining the TPS design as a part of the load-carry structure. The TPS design adopts POST/OTIS for trajectory optimization in the outmost design cycle, thereby leading to a final minimized fuel and weight objective. Meanwhile, ZONAIR and ASTROS* FEM will be driven by the automated parametric grid-panel mesh generation of AML, developed by Technosoft. By means of AML, an integrated software framework with a feature-based environment can be established to support an underlying object-orient architecture that will harness and link all software modules in configuration geometry, aerothermoelastic analysis, TPS design, etc. RMLS and HyperX are proposed as the test-bed vehicles for software demonstration in Phase II.