identification page- universityof texasat austin（识别页面——texasat奥斯汀大学）（13页）.docVIP

Table of Contents Table of Contents 1 Program Identification 2 Purpose of the Program 3 Mathematical Technique 3 Program Usage 5 Program Listing 5 Data Runs/Test Cases 6 Appendices 6 FLIGHT_ENVELOPE.m 7 NEWTON_RAPHSON.m 9 CDZERO.m 10 AK.m 10 F.m 11 ATMOS.m 12  Program Identification Program Assignment: flight_envelope.m Programmer: Gregory R. Whitney Date: October 15, 2001 Requirements: Computer running MATLAB Storage Required: 6.22KB Special I/O Requirements: None To plot the data, run flight_envelope.m. There is no input required. Purpose of the Program The purpose of the program is to plot the flight envelope for the Lear 23 at a weight of 11,000 pounds. The flight envelope is plotted for values of the altitude ranging from sea level to 45,000 ft. There are 5 curves generated by the program. The 5 curves are the stall velocity (Vstall), the maximum mach number (VM_Max), the maximum dynamic pressure velocity (VQ_max), the low speed solution (V_low), and the maximum speed solution (V_high). Mathematical Technique The mathematical technique of flight_envelope.m is based on the Newton-Raphson technique for solving for the roots of an equation. The equation that we need to find the roots of comes from the conditions of steady flight at full power setting. The equation is as follows: (1) where T = Thrust D = Drag H = altitude (constant) L = Lift (Constant) Substituting the equation for thrrust and drag produces the next equation: (2) where ρ* = Density at sea level ρ = Density M = Mach number Cdo = Zero-lift drag coefficient K = Induced drag factor W = Weight (W=L) V = velocity Sw = Planform Area The function f(V) can be solved for V numerically using the Newton-Raphson method. The Newton Raphson Method yields (3) (4) where ΔV = 0.0001 By solving Eq.[2] for V, we get two solutions: a high speed solution and a low speed solution. As altitude changes, you have to