The accumulation of ice on an airplane in flight is one of the leading contributing factors to general aviation accidents. To date, only relatively sophisticated methods based on detailed empirical data and flight data exist for its analysis. This paper develops a methodology and simulation tool for preliminary safety and performance evaluations of airplane dynamic response and climb performance in icing conditions. The important aspect of dynamic response sensitivity to pilot control input with the autopilot disengaged is also highlighted. Using only basic mass properties, configuration, propulsion data, and known icing data from a similar configuration, icing effects are applied to the dynamics of a non-real-time, six degree-of-freedom simulation model of a different, but similar, light airplane. Besides evaluating various levels of icing severity, the paper addresses distributed icing which consists of wing alone, horizontal tail alone, and unequal distributions of combined wing and horizontal tail icing. Results presented in the paper for a series of simulated climb maneuvers with various levels and distributions of ice accretion show that the methodology captures the basic effects of ice accretion on pitch response and climb performance, and the sensitivity of the dynamic response to pilot control inputs.
Nomenclature= arbitrary stability and control derivative C A iced = arbitrary stability and control derivative with icing effectshe has been with Texas A&M University, where he is Associate Professor of aerospace engineering, and Director of the Flight Simulation Laboratory. He teaches courses on digital control, nonlinear systems, flight mechanics, aircraft design, and coteaches a short course on digital flight control for the University of Kansas Division of Continuing Education. Current research interests include control of morphing air and space vehicles, multi-agent systems, intelligent autonomous control, vision-based navigation systems, and fault-tolerant adaptive control. He is an Associate Fellow of AIAA, past Chairman and current member of the Atmospheric Flight c = mean geometric chord D = carry through matrix f ice = icing factor g = gravitational acceleration h = integration step size I yy = airplane moment of inertia about y axis i = index, imaginary component j = imaginary component k 0 C A = coefficient icing factor constant L = roll angular acceleration M = pitch angular acceleration; modal matrix N = yaw angular acceleration nframes = number of time steps P = airplane body-axis roll rate p = perturbed airplane body-axis roll rate Q = airplane body-axis pitch rate q = perturbed airplane body-axis pitch rate _ q = airplane body-axis pitch acceleration q = dynamic pressure R = airplane body-axis yaw rate r = perturbed airplane body-axis yaw rate S = wing area t = time U = airplane velocity in the body-axis x direction U = control input vector u = stability axis airplane velocity in the x direction _ u = incremental change in the stability axis airplane velocity in the x direction V = airplane...