摘要
随着社会经济发展以及人民生活水平的提高,能源需求急剧增加,导致煤炭、石油等一次能源的过度开采应用,人类即将面临能源枯竭的威胁以及能源消耗产生的环境问题日益加剧带来的巨大压力,因此,新能源问题成为国家中长期发展规划的战略需求和世界的焦点。风能是最易于开发、发展最成熟的新能源。叶片作为风力发电设备的关键部件,其性能好坏决定了风电设备的发电效率。随着风力机单机容量增高,机组尺寸不断增大,叶片也越来越长,为了使之能够更好的捕捉风能,叶片扭转角也越大,不仅使得设计与制造困难,还增加了叶片对攻角敏感性,导致风力机输出功率更加不稳定,此外,叶片大型化也给运输与维护带来极大地不便,传统叶片越来越难以支持风力机的快速发展。本文利用圆柱形Magnus叶片结合翼型尾翼,提出了一种新型风力机叶片结构。
论文利用FLUENT软件仿真研究了高雷诺数高周速比情况下Magnus叶片绕流流场变化及其升/阻力特性,分析了风速、转速、叶片直径等参数对Magnus叶片升/阻力特性的影响规律,得到了影响Magnus叶片特性的最主要参数可以归纳为雷诺数和周速比,并分析了其影响规律。从理论上分析了叶片升阻比对风力机性能的影响规律,研究阐明了Magnus叶片应用于水平轴风力机所存在的问题。
研究了Magnus叶片在不同周速比和雷诺数条件下的流场变化规律,并分析了其流场形态变化产生的原因,得出了制约Magnus叶片升阻比提升的原因;研究了流场形态对Magnus叶片升/阻力特性的影响规律,得出了升阻比随前驻点和尾涡位置的变化规律,为Magnus叶片改进研究提供了理论支持。
提出了两种针对Magnus叶片的改进方案,增加头部结构和增加尾翼结构。分析了不同头部结构形状及其位置对流场及叶片升/阻力特性的影响;同时研究了尾翼结构参数对叶片气动性能影响规律,确定尾翼形状。根据尾翼结构几何参数对叶片流场影响,对叶片几何参数进行了优化。最后通过仿真计算,获得了风速、转速、攻角、周速比、雷诺数对新型叶片的升阻力系数影响规律。
设计了一套用于测量Magnus型风力机叶片气动性能的测试系统,利用该测试系统对圆柱形Magnus叶片和改进型Magnus叶片气动性能进行了实验研究,分析两种叶片升阻比受周速比、攻角等参数的影响,并与仿真结果进行了比较,分析了实验产生误差的原因并提出了改进建议。
应用改进型Magnus叶片设计了改进型Magnus水平轴风力机几何模型,根据叶素动量定理及其修正理论,建立了改进型Magnus水平轴风力机计算模型,通过该计算模型研究了不同参数下风力机气动特性变化规律,最后分析了改进型Magnus风力机功率特性、转矩特性和推力特性,并与传统翼型NACA0018构成的水平轴风力机和Magnus水平轴风力机进行了比较分析,得出了改进型Magnus风力机相比于另外两种风力机优势所在。
综上所述,本文提出了一种改进型Magnus叶片,升阻比可达到40以上,相比于Magnus叶片提升47%,同时相比于传统翼型叶片,具有攻角不敏感性,在较大攻角范围内不会出现失速现象。利用该叶片设计的改进型Magnus风力机具有较高的风能利用率和较好的自启动性能,功率输出较为稳定,同时叶片可以采用直叶片形式,使得叶片可以模块化设计与生产,大大降低了叶片设计制造难度以及运输与维护成本。改机型适用于大中型水平轴风力发电机建设,同时也可以应用于风速较低、风能资源相对贫乏地区的特殊风力机的开发。
With the social and economic development and the improvement of people'sliving standards, demand for energy increased dramatically, leading toover-exploitation of coal, oil and other primary energy applications, mankind isabout to face the enormous pressure of environmental problems caused by theenergy consumption. Therefore, new energy issues become the focus of the strategicneeds of the world's long-term development planning. Wind energy is the mostmature renewable energy. As a key component of the wind power generationequipment, blades' aerodynamic has important influence to power generationefficiency. As the power of wind turbine increasing, longer wind blades were need,at the same time how to ensure that optimize loads, reduce weight, enhance theadaptability to environment friendly and transport facilitation will becomeresearching hot topic of blade technology. With the increase of unit capacity of windturbine, size increasing of the unit, the blades are getting longer and longer, and inorder to enable it to capture the wind energy well the blade torsion angle becomegreater, which not only makes the design and manufacture of the blade moredifficult but also heighten the blade's sensitivity to the attack angle and lead to theinstability of turbine output power. What's more, the upsizing of blades brings thegreatly inconvenience of transportation and maintenance, which makes thetraditional blades become more and more difficult to support the rapidlydevelopment of the wind turbine. In this thesis, a new type of wind turbine bladewas studied, which integrated Magnus effect theory and advantage of airfoil.
The FLUENT software was used to study the flow field changes of rotatingcylinder in Magnus wind blade and it's lift/resistance characteristic at high Reynoldsnumber, and the rule of influence by flow velocity, rotating speed and cylinderdiameter was also analyzed, which has concluded that the most importantparameters in all those affact the Magnus blades characteristic can be summarized asthe Reynolds number and thetip-speed-ratio, and the rule of influence is alsoanalyzed. Theoretical analysis of the blade lift-drag ratio effect on the performanceof the wind turbine has been investigated. The studies found the problems of the theMagnus blade when applied to a horizontal axis wind turbine.
The variation of the flow field of Magnus wind blade within different Reynoldsnumber and thetip-speed-ratio was studied and the causes of the morphologicalchanges of the flow field was also analyzed; The formation mechanism of the resistance of Magnus bladess was studied and the constraints of its lift-to-drag ratioimprovement was derived; The influence law of flow field morphology to theMagnus blades lift/resistance characteristic was studied, and the variation law oflift-to-drag ratio with changes of the position of the front stagnation point and thetrailing vortex was obtained, which provides theoretical support for the the Magnusblades improvement research.
Analyzed principle of resistance formation when cylinder rotated in flow field,then proposed two programs to reduce the rotating cylindrical resistance; designedtwo programs improving structure characteristic for Magnus wind blade, which wereincreasing the head structure and the tail structure.; analyzed the impact of the flowfield by different head shape structures; then, influence of aerodynamic performanceby different parameters of tail structure were investigated, according to the results,optimized parameters of the new blade; at last, variation of lift and drag coefficientof the new blade by changing flow velocity, rotating speed, attack angle peripheralspeed to flow velocity ratio, Re numbers were investigated.
Designed a testing system used to measure aerodynamic performance of theMagnus type wind turbine blade, then experimental studied lift and drag coefficientof the new wind blades by the system; analyzed the experimental error according tocomparing the simulation results.
An improved Magnus horizontal axis wind turbine geometric model wasdesigned based on improved Magnus blades, and mathematical calculations modelof improved Magnus horizontal axis wind turbine was studied by using bladeelement momentum while the variation law of the turbine aerodynamic performancewithin different parameters was studied based on this model. In addition, the powercharacteristics, torque characteristics and thrust characteristics of improved Magnuswind turbine were analyzed, and a comparative analysis with horizontal axis windturbine consisting of traditional airfoil NACA0018and Magnus horizontal axis windmachine was conducted, which has obtained the advantage of improved Magnuswind turbine compared to the other two wind turbine.
In summary, a improvement wind turbine blade based on the Magnus effect hadbeen investigated in this thesis, the lift-to-drag ratio of the blade was exceed40,47%increased compared to the Magnus wind blade; at the same time, compared to theconventional airfoil blade, the improvement wind turbine blade had a better angle ofinsensitivity of attack. The improved Magnus wind turbine based on this windturbine blade has high wind energy utilization and better self-starting performance,and its power output is more stable, the blades can adopt the form of a straight bladeat the same time, what makes the design and production of the blade moer modular. What is more, the difficulty of blade design and manufacture as well astransportation and maintenance costs are also reduced. Horizontal axis wind turbineconstituted of the improvement bladese has a better aerodynamic performanceaccording to numerical study, it is suitable for large and medium-sized horizontalaxis wind turbine construction, and also can be applied to special wind turbinedevelopment in areas of low wind speed and poor wind energy resource.
引文
[1]李俊峰.2012中国风电发展报告[R].北京:中国环境科学出版社,2012:2-4.
[2]陈绍杰,申振华,徐鹤山.复合材料与风力机叶片[J].可再生能源,2008,26(2):90-91.
[3]谢校臻,王继辉.风机叶片设计制造发展趋势[C].玻璃钢/复合材料增刊(第十七届玻璃钢/复合材料学术年会论文集),2008:353-355.
[4]张希良.风能开发利用[M].北京:化学工业出版社,2004.
[5]倪受元.风力机的类型与种类[J].太阳能,2000,2:6-10.
[6]蒋超奇,严强.水平轴与垂直轴风力发电机的比较研究[J].上海电力,2007,2:163-165.
[7] Erich H. Wind turbines[M].Berlin: Springer,2005.
[8]田海姣,王铁龙,王颖.垂直轴风力发电机发展概述[M].应用能源技术,2006,(11):22-27.
[9] Price T J. James Blyth-Britain’s first modern wind power engineer. Windenergy [J],2005,29(3):191-200.
[10] Hagedorn L. Origins of wind power: Mr. Brush’s windmill dynamo [EB/OL].(2011-07-28)[2012-09-25]. http://www.renewableenergylawinsider.com/2011/07/28/origins-of-wind-power-mr-brushs-windmill-dynamo/.
[11] Shepherd D.G. Historical development of the windmill[M]. ScientificAmerican,1990,63(25):54.
[12] Hau E. Wind turbine [M]. Berlin: Springer,2005:23-64.
[13] Betz A. Introduction to the theory of flow machines [M]. Oxford: PergamonPress,1966:5-85.
[14] Trevor J P. Edward Golding's Iinfluence on Wind Power [J]. Wind Engineering,2005,29(6):513-530.
[15] Vestergaard J, Brandstrup L, Robert D. A brief history of the wind turbineindustries in Denmark and the United States[C]. Academy of InternationalBusiness (Southeast USA Chapter) Conference Proceeding. Florida: AIB,2004:322-327.
[16] Wikipedia. NASA wind turbine (EB/OL).(2012-05-11)[2012-09-25].http://en.wikipedia.org/wiki/NASA_wind_turbines.
[17]赵炜,李涛.国外风力发电机的现状及前景展望[M].电力需求侧管理,2009,11(2):77-80.
[18] Riegler H. HAWT versus VAWT: small VAWTs find a clear niche [J]. Refocus,2003,4(4):44–46.
[19] Shikha B T S, Kothari D P. Early development of modern vertical andhorizontal axis wind turbines: a review [J]. Wind Engineering,2005,29(3):287-299.
[20] Peace S. Another approach to wind. Mechanical Engineering,2004,126(6):28-31.
[21] Darrieus G J M. Turbine having its rotating shaft transverse to the flow of thecurrent [P]. US Patent, No.1.835.018,1931.
[22] Peace S. Why not vertical axis?[J]. Refocus,2003,4(3):30-33.
[23] Windside. WS series of VAWT (EB/OL).(2010-08-11)[2012-09-23].http://www. windside. com/.
[24] Mag-Wing. MW1100(EB/OL).(2009-05-11)[2012-09-23]. http://www.facebook. com/public/Mag-Wing.
[25] PACWIND. Pacwind Vertical Axis Wind Turbines (EB/OL).(2008-03-15)[2012-09-23]. http://en.wikipedia.org/wiki/PacWind.
[26] Stato Eolian. Gual industrie (EB/OL).(2009-07-15)[2012-09-23]. http://www.gual-industrie. com.
[27]陈云程,陈孝耀.朱成名.风力机设计与应用[M].上海科学技术出版社,1990.
[28] Timmer W A. An overview of NACA6-digit airfoil series characteristics withreference to airfoils for large wind turbine blades [C]. Aerospace SciencesMeeting Including the New Horizons Forum and Aerospace Exposition.Florida: AIAA,2009:1-13.
[29] Tangler J L, Somer D M. NREL airfoils families for HAWTs [R]. Washington:AWRA,1995:1-9.
[30] Timmer W A, Rooij V. Summary of the Delft University wind turbinededicated airfoils [R]. Virginia: AIAA,2003:1-10.
[31] Fuglsang P, Bak C. Development of the Ris wind turbine airfoils [J]. WindEnergy,2004,(7):145-162.
[32] Wikipedia. XFOIL [EB/OL].(2012-01-26)[2012-09-25].http://en.wikipedia.org/wiki/XFOIL.
[33] MIT. XFOIL [EB/OL].(2008-04-07)[2012-09-26]. http://web.mit.edu/drela/Public/web/xfoil/.
[34]张石强.风力机专用翼型及叶片关键设计理论研究.重庆大学博士学位论文,2010:124-126
[35]叶枝全,黄继雄,陈严等.风力机新系列翼型气动性能研究[J].太阳能学报,2002,23(2):211-216.
[36]张果宇,冯卫民,刘长陆等.6种风力机叶片翼型的气动性能数值模拟研究[J].可再生能源,2009,27(2):11-15.
[37]陈进,王旭东,沈文忠等.风力机叶片的形状优化设计.机械工程学报.2010,46(3),131-134
[38]刘雄,陈严,叶枝全.水平轴风力机风轮叶片优化设计模型研究[J].汕头大学学报,2006,21(1):44-49
[39]王旭东,陈进,张石强等.风力机叶片翼型型线集成设计理论研究[J].中国机械工程.2009,20(2),211-213
[40]张仲柱,王会社,赵晓路等.水平轴风力机叶片气动性能研究.工程热物理学报.2007,28(5),781-783
[41]刘旺玉,张海全,李明刚.基于仿生设计的风力发电机叶片力学性能的实验研究.实验力学.2009,24(2),121-126
[42]廖猜猜,徐建中,席光.基于一种改进的PSO算法在风力机叶片优化中的应用.工程热物理学报,2008,29(5),773-776
[43]白井艳,杨科,黄哀武等.风力机专用翼型的反设计.工程热物理学报,2012,33(11),1885-1887
[44]郭婷婷,吴殿文,王成荫.风力发电机叶片预弯设计及其数值研究.动力工程学报.2010,30(6),450-455
[45] Modi V J. Moving surface boundary-layer control: A review [J]. Journal ofFluids and Structures,1997,11:627-663.
[46] Bery H, Yao Y X. Numerical Simulation of Unsteady Flow to ShowSelf-starting of Vertical Axis Wind Turbine Using Fluent [J]. Journal ofApplied Sciences,2011,(11):962-970.
[47] Bery H, Yao Y X. Computatinal Analysis of vertical axis wind turbine withopen-able airfoil [C]. Power and Energy Engineering Conference (APPEEC),Wuhan: IEEE,2011:1-5.
[48]申振华,于国亮.翼型弯度对风力机性能的影响[J].动力工程,2007,27(1):136-139.
[49]唐进.提高风力机叶型气动性能的研究[D].清华大学,2004:65-78.
[50]庄月晴,黄典贵.前缘旋转的风力机翼型气动性能研究[J].工程热物理学报,2011,32(1):43-46.
[51] Seifert J. A review of the Magnus effect in aeronautics [J]. Progress inAerospace Sciences,2012,55:17-45.
[52] Swanson W M. The Magnus effect: a summary of investigations to date [J].Journal of Basic Engineering Transactions of ASME,1961,83(3):470.
[53] Alaways L W, Hubbard M. Experimental determination of baseball spin andlift [J]. Journal of Sports Sciences.2001,19(5):349-358.
[54] Fleming P D, Probert S D. The evolution of wind-turbines: An historicalreview. Applied Energy,1984,18(3):163-177.
[55] Mehta R D, Pallis J M. The aerodynamics of a tennis ball [J]. SportsEngineering,2001,4(4):177-189.
[56] Cayzac R, Carette E, Denis P, et al. Magnus Effect: Physical Origins andNumerical Prediction[J]. Transactions of the ASME-E-Journal of AppliedMechanics,2011,78(5):051005.
[57] Thorn A. Experiments on the flow past a ratating cylinder[D]. Glasgow:University of Glasgow.1931:1-26
[58] Mishuris G, Plakhov A. Magnus Effect and dynamics of a spinning disc in ararefied medium[J]. Archives of Mechanics,2009,61(5):391-406.
[59] Glauert M B. The flow past a rapidly rotating circular cylinder [J].Proceedings of the Royal Society of London. Series A, Mathematical andPhysical Sciences,1957,242(1228):108-115
[60] Ludwig G R. An experimental investigation of laminar separation from amoving wall [J]. AIAA Paper,1964,(64):6.
[61] Diaz F, GavaldàJ, Kawall J G, Keffer J G and Giralt F. Vortex shedding from aspinning cylinder [J]. Physics of Fluids,1983,26(12):3454-3460.
[62] Diaz F, GavaldàJ, Kawall J G, Keffer J G and Giralt F. Asymmetrical wakegenerated by a spinning cylinder [J]. AIAA Journal,1985,23(1):49-54.
[63] Coutanceau M and Menard C. Influence of rotation on the near-wakedevelopment behind an impulsively started circular cylinder [J]. Journal ofFluid Mechanics,1985,158:399-446.
[64] Badr H M and Dennis S C R. Time-dependent viscous flow past animpulsively started rotating and translating circular cylinder [J]. Journal ofFluid Mechanics,1985,158:447-488.
[65] Badr H M, Dennis S C R and Young P J S. Steady and unsteady flow past arotating circular cylinder at low Reynolds numbers [J]. Journal of FluidMechanics,1989,17:579-609.
[66] Badr H M, Coutanceau M, Dennis S C R and Ménard C. Unsteady flow past arotating circular cylinder at Reynolds number1e3and1e4[J]. Journal of FluidMechanics,1990,220:459-484.
[67] Chew Y T, Cheng M and Luo S C. A numerical study of flow past a rotatingcircular cylinder using a hybrid vortex scheme [J]. Journal of Fluid Mechanics,1995,299:35-71.
[68] Nair M T, Sengupta T K and Chauhan U S. Flow past rotating cylinders at highReynolds numbers using higher order upwind scheme [J]. Computers&Fluids,1998,27(1):47-70.
[69] Kang S, Choi H, and Lee S. Laminar flow past a rotating circular cylinder [J].Physics of Fluids,1999,11(11):3312-3321.
[70] Mittal S. Flow past rotating cylinders: effect of eccentricity [J]. Transactionsof the ASME. Journal of Applied Mechanics,2001,68(4):543-552.
[71] Mittal S and Kumar B. Flow past a rotating cylinder [J]. Journal of FluidMechanics,2003,476:303-334.
[72] Mittal S. Three-dimensional instabilities in flow past a rotating cylinder [J].Journal of Applied Mechanics,2004,71:89-95.
[73] Karabelas S J. Large Eddy Simulation of High-Reynolds Number flow past arotating cylinder [J]. International Journal of Heat and Fluid Flow.2010,(31):518-527.
[74] Dihj R. Flettner-Rotorschiffe: Alte Technik für neue Schiffe [J]. HANSAInternational Maritime Journal,2007,12:16-22.
[75] Reu H J. Flettner-Rotorschiffe—Alte Technik fur neue Schiffe[J].International Maritime Journal,2007,144(12):16–22.
[76] Calderon A A and Arnold F R. A study of the aerodynamic characteristics of ahigh-lift device based on a rotating cylinder and flap [Technical ReportRCF-1]. Stanford University,1961.
[77] Cichy D R, Harris J W and Mackay J K. Flight tests of a rotating cylinder flapon a North American Rockwell YOV-10aircraft, NASA: CR-2135,1972.
[78] Weiberg J A and Gamse B. Large-scale wind-tunnel tests of an airplane modelwith two propellers and rotating cylinder flaps. NASA Technical Note,TN-D-4489,1968.
[79] Weiberg J A and Dickinson S O. Large-scale wind-tunnel tests of an airplanemodel with four propellers and rotating cylinder flaps. NASA Technical Note,TN-D-5742,1970.
[80]李岩.垂直轴风力机及其发展概况.可再生能源,2009,27(1):121-123.
[81] Hanson T F. Magnus air turbine system. United State:262136(P),1981-5-11.
[82] Bychkov N M. Magnus wind turbine.3. Calculated characteristics of thewindwheel [J]. Thermophysics and Aeromechanics,2008,15(2):321.
[83] Akiyoshi I, Chisachi K and Nobuhiro M. Evaluation of aerodynamic propertiesof Magnus wind turbines with spiral fins [J]. Nihon Kikai Gakkai Ronbunshu,2010,(76):380-382.
[84] Mecaro. Newest Model Installation Report [EB/OL].(2007-09-18)[2012-09-28]. http://www.mecaro.jp/eng/new_repo.html.
[85] Luo D H, Huang D G and Wu G Q. Analytical solution on Magnus windturbine power performance based on the blade element momentum theory [J].Renewable and Sustainable Energy,2011,(3):1-14.
[86]武鑫.风能技术[M].北京:科学出版社,2009.
[87]钱翼稷.空气动力学[M].北京航空航天大学出版社,2004.
[88] Happel J, Brenner H. Low Reynolds number hydro dynamics with specialApplication to particulate media, New Jersey: Prentice-Hall,1965.
[89]陈卓如.工程流体力学[M].北京:高等教育出版社,2004:346-371.
[90]林建忠.流体力学.北京:清华大学出版社,2005.
[91]王福军.计算流体力学分析—CFD软件原理与应用[M].北京:清华大学出版社,2006:13-22.
[92]韩占忠,王敬,兰小平. Fluent流体工程仿真计算实例与应用[M].北京:北京理工大学出版社,2004:20-21.
[93]李进良,李承曦,胡仁喜.精通FLUENT6.3流场分析[M].北京:化学工业出版社,2009:50-80.
[94]李友荣.直接数值模拟(DNS)的应用现状[C].第十四届全算流体力学会议论文集.贵阳:中国空气动力学会,2009:990-993.
[95]万曦.湍流大涡数值模拟的现状和展望[J].航空科学技术,2011,(2):55-56.
[96]江帆. Fluent高级应用与实例分析[M].北京:清华大学出版社,2008:22-26.
[97] Chew Y T, Luo S C, Cheng M. Numerical study of a linear shear flow past arotating cylinder. Journal of Wind Engineering and Industrial Aerodynamics1997,66:107-125.
[98]牛山泉.风能技术[M].北京:科学出版社,2009:43-45.
[99] Manwell J F, Mcgowan J G, Rogers A L. Wind Energy Explained [M]. JohnWiley&Sons,2002.
[100] David A S. Wind Turbine Technology[M]. New York: ASME PRESS,1994.
[101]施鹏飞.风力机的理论和设计[M].北京:机械工业出版社,1987.
[102] Schlichting H. Boundary-layer theory [M]. New York: McGraw-Hill,1997.
[103]苏铭德,黄素逸.计算流体力学基础.北京:清华大学出版社,1997.
[104] Kundu P K, Cohen I M. Fluid Mechanics[M]. California: Elsevier AcademicPress,2004:110-115.
[105] Karabelas S J, Koumroglou B C, Argyropoulos C D. High Reynolds numberturbulent flow past a rotating cylinder[J]. Applied Mathematical Modelling,2012,(36):379-398.
[106] Moore F K. Boundary layer research[M]. Berlin: Springer Verlag,1958:296-310.
[107] Rott N. Unsteady viscous flow in the vicinity of a stagnation point[J].Quarterly Journal of Applied Mathematics,1956,(13):444-451.
[108] Sears W R. Some recent development in airfoil theory[J]. Journal of theAeronautical Sciences,1956,(23):490-499.
[109] Aldoss T K, Abou-Arab T W. Experimental study of the flow around a rotatingcylinder in crossflow[J]. Experimental Thermal and Fluid Science,1990,(3):316-322.
[110] Aldoss T K, Mansour Awad. Theoretical calculations of the flow around arotating circular cylinder placed in a uniform flow[J]. Journal of FluidsEngineering,1988,(110):96-98.
[111] Inoue O. MRS criterion for flow separation over moving walls[J]. AIAAJournal,1981,9(19):1108-1111.
[112] Peller H. Thermofluiddynamic experiments with a heated and rotating circularcylinder in crossflow, Part2.1Boundary layer profiles and location ofseparation points[J]. Experiments in Fluids,1986,(4):223-231.
[113]王勇,田波平.概率论与数理统计[M].北京:科学出版社,2005:93-105.
[114]叶枝全,黄继雄,陈严等.适用于风力机的新翼型气动性能的实验研究[J].太阳能学报,2003,24(4):548-554.
[115]贺德馨.风工程与工业空气动力学[M].北京:国防工业出版社,2006:83-85.
[116] Burton T, Sharpe D, Jenkins N. Wind Energy Handbook[M]. Chichester: JohnWiley&Sons,2001:84-87.
[117]贺德馨,张维智.风力机翼型大攻角气动性能手册.中国空气动力研究与发展中心,1998.