In this paper the flight envelope partition method is proposed. According to the atmospheric conditions and engine regulation, flight envelope is parted to 83 regions. In each region parameters in region and at boundary points are replaced with parameters at centroid, which are constrained with dynamic quotas of engine dynamic equation. Finally the mathematical simulation of a turbofan engine performance in whole flight envelope is presented, and the results show the method is effective in establishing engine dynamic equation in whole flight envelope.