Aircraft failures alter the aircraft dynamics and cause maneuvering flight envelope to change. Such envelope variations are nonlinear and generally unpredictable by the pilot as they are governed by the aircraft complex dynamics. Hence, in order to prevent in-flight Loss of Control it is crucial to practically predict the impaired aircraft's flight envelope variation due to any a-priori unknown failure degree. This paper investigates the predictability of the number of trim points within the maneuvering flight envelope and its centroid using both linear and nonlinear least-squares estimation methods. To do so, various polynomial models and nonlinear models based on hyperbolic tangent function are developed and compared which incorporate the influencing factors on the envelope variations as the inputs and estimate the centroid and the number of trim points of the maneuvering flight envelope at any intended failure degree. Results indicate that both the polynomial and hyperbolic tangent function-based models are capable of predicting the impaired fight envelope variation with good precision. Furthermore, it is shown that the regression equation of the best polynomial fit enables direct assessment of the impaired aircraft's flight envelope contraction and displacement sensitivity to the specific parameters characterizing aircraft failure and flight condition.Nomenclature V = total airspeed, knot , = angle of attack, sideslip angle, respectively, deg p, q, r = angular velocity components (roll rate, pitch rate, yaw rate, respectively), deg/s 1 Ph.D. Candidate, Faculty of New Sciences and Technologies; ramin.norouzi@ut.ac.ir. Student Member AIAA. 2 Associate Professor, Faculty of New Sciences and Technologies; kosari_a@ut.ac.ir. 3 Assistant Professor, Faculty of New Sciences and Technologies; sabourmh@ut.ac.ir. , , , = Euler angles (roll, pitch, yaw, respectively), flight path angle, deg , = state vector, control vector, respectively , , , ℎ = deflection angles (elevator, aileron, rudder, respectively), deg, throttle setting (%) ∈ [0, 1] , �, = fitted model's residual, outcome, parameters (coefficients) vectors, respectively , = predictor vector, model parameters (coefficients) vector, respectively = total degree of polynomial , ℍ, , = gradient, Hessian, Jacobian of the objective function , identity matrix, respectively = Newton's method search direction = Marquardt parameter = hyperbolic tangent function 1 , 2 = total number of neurons in the first (hidden) layer, second (last) layer, respectively = backpropagation training iteration ℒ , ℒ , ℒ , ℒ, = (net output, net input, bias) of ℒ ℎ neuron in layer , respectively ℒ,= weight's element of ℒ ℎ neuron in layer due to ℎ neuron of the previous layer .̂ = standard backpropagation sensitivity, Marquardt sensitivity, respectively ℰ, = contribution of the ℰ ℎ neuron of the last (second) layer to the error of the ℎ training sample = statistical expected value of the ℎ input factor Subscripts , = total number of training samples, total number of model parameters, res...