Hypersonic vehicle control system design and simulation require models that contain a low number of states. Modeling of hypersonic vehicles is complicated due to complex interactions between aerodynamic heating, heat transfer, structural dynamics, and aerodynamics. Although there exist techniques for analyzing the effects of each of the various disciplines, these methods often require solution of large systems of equations, which is infeasible within a control design and evaluation environment. This work presents an aerothermoelastic framework with reducedorder aerothermal, heat transfer, and structural dynamic models for time-domain simulation of hypersonic vehicles. Details of the reduced-order models are given, and a representative hypersonic vehicle control surface used for the study is described. The methodology is applied to a representative structure to provide insight into the importance of aerothermoelastic effects on vehicle performance. The effect of aerothermoelasticity on total lift and drag is found to result in up to an 8% change in lift and a 21% change in drag with respect to a rigid control surface for the four trajectories considered. An iterative routine is used to determine the angle of attack needed to match the lift of the deformed control surface to that of a rigid one at successive time instants. Application of the routine to different cruise trajectories shows a maximum departure from the initial angle of attack of 8%.corresponding to ith column of A C = correlation matrix c = modal coordinate of thermal proper orthogonal decomposition basis vector Cx = correlation model for kriging c p = specific heat at constant pressure d = structural modal coordinates, kriging sample point E = modulus of elasticity F = thermal load vector of full system in physical space f = generalized thermal load vector of reduced system in modal space F s = structural load vector of full system in physical space f s = generalized structural load vector of reduced system in modal space H i = coefficient matrices for integration of equations of motion h i = thickness of ith layer of thermal protection system K = thermal conductivity matrix of full system in physical space k = generalized thermal conductivity matrix of reduced system in modal space K G = geometric stiffness matrix K s = structural stiffness matrix k s = generalized stiffness matrix of reduced system in modal space k T = thermal conductivity of material K s = modified structural stiffness matrix L = aerodynamic lift, length M = thermal capacitance matrix of full system in physical space, Mach number m = generalized thermal capacitance matrix of reduced system in modal space M s = structural mass matrix of full system in physical space m s = generalized mass matrix of reduced system in modal space n = number of proper orthogonal decomposition snapshots n e = number of snapshots for kriging evaluation cases n k = number of kriging snapshots n s = number of structural parameters in reduced-order aerothermodynamic model n t = number of thermal parameters...