CN108583847A - A kind of low reynolds number high power factor aerofoil profile suitable for long endurance unmanned aircraft - Google Patents
A kind of low reynolds number high power factor aerofoil profile suitable for long endurance unmanned aircraft Download PDFInfo
- Publication number
- CN108583847A CN108583847A CN201810382463.7A CN201810382463A CN108583847A CN 108583847 A CN108583847 A CN 108583847A CN 201810382463 A CN201810382463 A CN 201810382463A CN 108583847 A CN108583847 A CN 108583847A
- Authority
- CN
- China
- Prior art keywords
- aerofoil profile
- reynolds number
- power factor
- airfoil
- chord length
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C3/00—Wings
- B64C3/10—Shape of wings
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Aviation & Aerospace Engineering (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
The present invention provides a kind of suitable for long endurance unmanned aircraft low reynolds number high power factor aerofoil profile, disclose maximum camber and its position of the aerofoil profile, maximum gauge and its position, the equation and aerofoil profile coordinate that upper and lower surface meets respectively, f is maximum camber, Xf is the abscissa value at aerofoil profile maximum camber, t is maximum gauge, and Xt is the abscissa value at aerofoil profile maximum gauge, and C is chord length.The origin of coordinate system is airfoil mean line leading edge point where defining aerofoil profile, and X-axis is overlapped with string, and airfoil trailing edge is directed toward in direction by aerofoil profile leading edge, and Y-axis is directed toward the direction of airfoil mean line bending perpendicular to X-axis.By carrying out rational modification to the upper and lower surface gradient distribution of aerofoil profile, delay turning point Forward bubble separation, it may be implemented under the operating Reynolds number of setting, the aerofoil profile has more preferably performance, i.e. larger range of lift resistance ratio, higher power factor coordinates Aircraft Conceptual Design, improves the cruise duration of unmanned plane.
Description
Technical field
The invention belongs to Unmanned Aerial Vehicle Airfoil design fields, and in particular to a kind of low reynolds number high power factor aerofoil profile.
Background technology
Long endurance unmanned aircraft, such as small-sized solar unmanned plane and Altitude Long Endurance Unmanned Air Vehicle, flying speed is small, flight thunder
Promise number is low, and cruise Reynolds number is about in 20-40 ten thousand or so, and mostly referred to as low reynolds number (claims in the industry 100,000-100 ten thousand Reynolds number
For low reynolds number), under the reynolds number range, relate generally to low reynolds number, the Laminar Flow of low speed, laminar flow bubble, laminar flow point
From and turn twist etc. flowing physical phenomenons, with powered by conventional energy middle and high Reynolds number flying condition under of the laminar flow to turbulent flow have very
It is different.Thus determine that such long endurance unmanned aircraft Aerodynamic Design method has obviously with thinking with middle and high Reynolds number aircraft
Difference.And solar energy unmanned plane also needs to consider the compound of photovoltaic module, therefore airfoil profile also suffers restraints.
Although there is part to can be used for low Reynolds number airfoil under low reynolds number state of flight in existing aerofoil profile library, for example originally set
The SD7032 of aerofoil profile as a comparison in meter, under the Reynolds number of 200,000 magnitudes, top airfoil turning point with the increase of the angle of attack and
Rapid forward edge movement, detaches, causing property loss of energy so as to cause air-flow.
Invention content
In view of the deficiencies of the prior art, performance bad situation when being matched with practical work for existing low Reynolds number airfoil,
The present invention provides a kind of low reynolds number high power factor aerofoil profile, can be well matched with the low flight Reynolds number of long endurance unmanned aircraft
Actual condition, meanwhile, upper surface bending degree is small, is suitable for the matching of the photovoltaic module of solar energy unmanned plane.
In order to realize that the above technical purpose, the present invention use following technical scheme:
Specific flying condition is analyzed, such as assumes that flight Reynolds number is 200,000.So under this condition, it first tentatively picks out often
The low Reynolds number airfoil seen, under specified operating mode, with the aerodynamic analysis technology of aircraft industry, analysis different airfoil profiles are being set
Performance under operating mode.
Since the main performance of aerofoil profile is heavily dependent on its shape upper and lower surfaces of, by further grinding
Study carefully performance capabilities preferably several groups of aerofoil profiles in back, by pointedly changing leading-edge radius, upper and lower surface gradient is distributed, most
The parameters such as heavy camber position then find out the best aerofoil profile of performance capabilities, as designing airfoil, and pass through experimental verification and making
Practical unmanned plane Flight.
Above-mentioned aerodynamic analysis is primarily referred to as calculating under setting operating mode, using the lift resistance ratio of the aircraft complete machine of specific aerofoil profile,
The common performance parameter such as power factor.
According to above-mentioned design principle, the present invention specifically provides low reynolds number high power factor aerofoil profile, and is named
For DMXfoil-1, for being described as follows for aerofoil profile:
The ratio between the maximum camber f and chord length C of aerofoil profile f/C=3.88%,
The maximum camber position Xf/C=48.09 of aerofoil profile,
The ratio between the maximum gauge t and chord length C of aerofoil profile t/C=10.50%,
The maximum gauge position of aerofoil profile is Xt/C=28.99%;
When definition chord length is 1, the equation of aerofoil profile upper surface is:
Y=-2.687263401080x6+9.080382811085x5-12.119920168585x4+8.215218339892x3-
3.252320013999x2+0.754048283792x+0.008814393076
When definition chord length is 1, the equation of aerofoil profile lower surface is:
Y=3.108344757931x6-10.324958446486x5+13.345995816048x4-8.586501844602x3+
2.878571621603x2-0.414971953131x-0.005291297688
Wherein, Xf is the abscissa value at aerofoil profile maximum camber, and Xt is the abscissa value at aerofoil profile maximum gauge.Define aerofoil profile institute
It is airfoil mean line leading edge point in the origin of coordinate system, X-axis is overlapped with string, and airfoil trailing edge, Y-axis are directed toward in direction by aerofoil profile leading edge
The direction of airfoil mean line bending is directed toward perpendicular to X-axis.
When the chord length of the aerofoil profile is 1, then the coordinate corresponding to the upper and lower surface of aerofoil profile is as follows:
Upper table areal coordinate:
Following table areal coordinate:
Description of the drawings
Fig. 1 is the composition figure of designing airfoil, and in figure, 1 is aerofoil profile leading edge, and 2 be aerofoil profile upper surface, and 3 be aerofoil profile lower surface, 4
It is airfoil trailing edge for aerofoil profile string, 5, f is aerofoil profile maximum camber, and Xf is the abscissa value at aerofoil profile maximum camber, and C is chord length, t
For aerofoil profile maximum gauge, Xt is the abscissa value at aerofoil profile maximum gauge.
Fig. 2 is the geometry external form comparison diagram of designing airfoil DMXfoil-1 and comparison aerofoil profile SD7032.
Fig. 3 be designing airfoil DMXfoil-1 under design point pressure distribution curve figure (Re=20 ten thousand, angle of attack=
5.6 °, low speed).
Fig. 4 be compare aerofoil profile SD7032 under design point pressure distribution curve figure (Re=20 ten thousand, angle of attack=5.6 °,
Low speed).
Fig. 5 is the lift resistance ratio characteristic curve comparison diagram (Re=20 of designing airfoil DMXfoil-1 and comparison aerofoil profile SD7032
Ten thousand, low speed), in figure, 6 be the aerodynamic characteristic curve graph under design point of designing airfoil DMXfoil-1, and 7 be comparison aerofoil profile
The aerodynamic characteristic curve graph under design point of SD7032.
Fig. 6 is that designing airfoil DMXfoil-1 and comparison aerofoil profile SD7032 is respectively adopted with comparison aerofoil profile in design point
Under power factor Character Comparison figure, in figure, 6 be the aerodynamic characteristic curve under design point of designing airfoil DMXfoil-1
Figure, 7 be the aerodynamic characteristic curve graph under design point for comparing aerofoil profile SD7032.(Re=20 ten thousand, low speed)
Specific implementation mode
Below in conjunction with attached drawing, the present invention is described in detail:
A kind of low reynolds number high power factor aerofoil profile suitable for long endurance unmanned aircraft of the present invention, design principle are:
Specific flying condition is analyzed, such as in the present embodiment, at the beginning of analysis, that is, assumes that flight Reynolds number is 20W.In the condition
Under, common low Reynolds number airfoil is tentatively first picked out, under specified operating mode, with the aerodynamic analysis technology of aircraft industry,
Analyze performance of the different airfoil profiles in the case where setting operating mode.
Since the main performance of aerofoil profile is heavily dependent on its shape upper and lower surfaces of, by further grinding
Study carefully performance capabilities preferably several groups of aerofoil profiles in back, by pointedly changing leading-edge radius, upper and lower surface gradient is distributed, most
The parameters such as heavy camber position then find out the best aerofoil profile of performance capabilities, as designing airfoil.
Above-mentioned aerodynamic analysis is primarily referred to as calculating under setting operating mode, using the lift resistance ratio of the aircraft complete machine of specific aerofoil profile,
The common performance parameter such as power factor.
In the following, being given by the rule characteristic for the low reynolds number high power factor aerofoil profile that the optimization method obtains, give
The designing airfoil regularity of distribution upper and lower surfaces of, the equation and shape of satisfaction, can refer to Fig. 1.
According to above-mentioned design principle, the present invention specifically provides low reynolds number high power factor aerofoil profile, and is named
For DMXfoil-1, for being described as follows for aerofoil profile:
The ratio between maximum camber f and chord length C of aerofoil profile
F/C=3.88%,
The maximum camber position of aerofoil profile
Xf/C=48.09,
The ratio between maximum gauge t and chord length C of aerofoil profile
T/C=10.50%
The maximum gauge position of aerofoil profile is
Xt/C=28.99%;
When definition chord length is 1, the equation of aerofoil profile upper surface is:
Y=-2.687263401080x6+9.080382811085x5-12.119920168585x4+8.215218339892x3-
3.252320013999x2+0.754048283792x+0.008814393076
When definition chord length is 1, the equation of aerofoil profile lower surface is:
Y=3.108344757931x6-10.324958446486x5+13.345995816048x4-8.586501844602x3+
2.878571621603x2-0.414971953131x-0.005291297688
Wherein, Xf is the abscissa value at aerofoil profile maximum camber, and Xt is the abscissa value at aerofoil profile maximum gauge.Define aerofoil profile institute
It is airfoil mean line leading edge point in the origin of coordinate system, X-axis is overlapped with string, and airfoil trailing edge, Y-axis are directed toward in direction by aerofoil profile leading edge
The direction of airfoil mean line bending is directed toward perpendicular to X-axis.
When the chord length of the aerofoil profile is 1, then the coordinate corresponding to the upper and lower surface of aerofoil profile is as follows:
Upper table areal coordinate:
Following table areal coordinate:
It is special to have carried out a series of Pneumatic Calculations and reality with comparison aerofoil profile SD7032 in order to illustrate the progressive of the design aerofoil profile
Verification, it is as follows to obtain result:
Reynolds number is set as 200,000, when the aerofoil profile angle of attack is 5.6 °, the lift resistance ratio of designing airfoil DMXfoil-1 reaches maximum, meter
It is 83.8 to calculate result, and about at 52% chord length, and under same Reynolds number, the lift resistance ratio of SD7032 is only top airfoil turning point
62.3, this is because for SD7032 aerofoil profiles, under this angle of attack, turning point is very close to leading edge locus, since energy subtracts
It is few, it is not enough to be attached to upper surface, airflow diversion has occurred, resistance is made to increase severely, die-off so as to cause lift resistance ratio.
Under above-mentioned condition, the pressure distribution curve figure of DMXfoil-1 and SD7032 are attached to Fig. 3 and Fig. 4 respectively.
And in the angle-of-attack range of bigger (Fig. 5), it is not difficult to find out, the angle of attack is designing airfoil within the scope of -5 ° -3 °
The lift resistance ratio characteristic of DMXfoil-1 and comparison aerofoil profile SD7032 are close, and when the angle of attack is within the scope of 3 ° -20 °, designing airfoil
The lift resistance ratio of DMXfoil-1 is better than comparison aerofoil profile SD7032.
Further, (major parameter is wing chord length 400mm, the span for our designs to same solar energy unmanned plane
5.4m, RECTANGULAR WINGS add wingtip), designing airfoil DMXfoil-1 and comparison aerofoil profile SD7032 is respectively adopted, compares setting flight reappearance
Under the conditions of power factor, the results are shown in Figure 6, also can be seen that, except when the wing angle of attack be 1 ° when power factor having the same
Outside, under remaining angle of attack, the power factor of designing airfoil DMXfoil-1 is above comparison aerofoil profile SD7032.Pay special attention to cruise shape
State, that is, range of angles of attack, between 4 ° -6 °, the difference of power factor is more apparent, the cruising power breath of this value and aircraft
Manner of breathing closes, and directly affects the cruise duration of aircraft.In the case where flight reappearance is certain, high power factor advantageously reduces
Cruising power, and then improve cruise duration.For long endurance unmanned aircraft, there is obvious significance.
Claims (2)
1. a kind of low reynolds number low-power factor aerofoil profile, this aerofoil profile are named as DMXfoil-1, it is characterised in that:Aerofoil profile is most
The ratio between heavy camber f and chord length C f/C=3.88%, maximum camber position Xf/C=48.09%, the maximum gauge t and chord length C of aerofoil profile
The ratio between t/C=10.50%, maximum gauge position be Xt/C=28.99%;
When definition chord length is 1, the equation of aerofoil profile upper surface is:
y = -2.687263401080 x6 + 9.080382811085 x5 - 12.119920168585 x4 +
8.215218339892 x3 - 3.252320013999 x2 + 0.754048283792 x + 0.008814393076
When definition chord length is 1, the equation of aerofoil profile lower surface is:
y = 3.108344757931 x6 - 10.324958446486 x5 + 13.345995816048 x4 -
8.586501844602 x3 + 2.878571621603 x2 - 0.414971953131 x - 0.005291297688
Wherein, Xf is the abscissa value at aerofoil profile maximum camber, and Xt is the abscissa value at aerofoil profile maximum gauge,
The origin of coordinate system is airfoil mean line leading edge point where defining aerofoil profile, and X-axis is overlapped with string, and direction is referred to by aerofoil profile leading edge
To airfoil trailing edge, Y-axis is directed toward the direction of airfoil mean line bending perpendicular to X-axis.
2. a kind of low reynolds number high power factor aerofoil profile suitable for long endurance unmanned aircraft according to claim 1, special
Sign is:When the chord length of the aerofoil profile is 1, then the coordinate corresponding to the upper and lower surface of aerofoil profile is as follows:
Upper table areal coordinate:
Following table areal coordinate:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810382463.7A CN108583847B (en) | 2018-04-26 | 2018-04-26 | Low-Reynolds-number high-power-factor wing section suitable for long-endurance unmanned aerial vehicle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810382463.7A CN108583847B (en) | 2018-04-26 | 2018-04-26 | Low-Reynolds-number high-power-factor wing section suitable for long-endurance unmanned aerial vehicle |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108583847A true CN108583847A (en) | 2018-09-28 |
CN108583847B CN108583847B (en) | 2020-08-11 |
Family
ID=63610024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810382463.7A Active CN108583847B (en) | 2018-04-26 | 2018-04-26 | Low-Reynolds-number high-power-factor wing section suitable for long-endurance unmanned aerial vehicle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108583847B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109969381A (en) * | 2019-04-12 | 2019-07-05 | 中国空气动力研究与发展中心计算空气动力研究所 | A kind of low reynolds number airfoil with high ratio of lift over drag with sharp leading edge heavy camber |
CN110386243A (en) * | 2019-07-22 | 2019-10-29 | 北京耐威科技股份有限公司 | A kind of low-speed unmanned aerial vehicle aerofoil profile of high lift-drag ratio |
CN110435873A (en) * | 2019-08-15 | 2019-11-12 | 西北工业大学 | A kind of half blended wing-body anury formula Unmanned Aerial Vehicle Airfoil race cruised from trim |
CN111498084A (en) * | 2020-04-15 | 2020-08-07 | 成都飞机工业(集团)有限责任公司 | Low-resistance laminar flow airfoil applied to high-altitude high-speed long-endurance unmanned aerial vehicle |
CN112977815A (en) * | 2021-05-10 | 2021-06-18 | 北京三快在线科技有限公司 | Rotor craft, blade of rotor craft and wing section of blade |
CN114194373A (en) * | 2021-12-03 | 2022-03-18 | 中电科芜湖通用航空产业技术研究院有限公司 | Aircraft wing and aircraft |
CN114718903A (en) * | 2022-04-19 | 2022-07-08 | 成都航空职业技术学院 | High-performance wing section for heat dissipation axial flow fan |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0068121B1 (en) * | 1981-07-01 | 1985-08-28 | Dornier Gmbh | Air foil for short take-off and landing aircraft |
CN103482054A (en) * | 2013-08-14 | 2014-01-01 | 西北工业大学 | Low-Reynolds-number wing section matched with full-wing solar unmanned aerial vehicle |
CN104354850A (en) * | 2014-10-30 | 2015-02-18 | 中国人民解放军空军航空大学 | Airfoil profile of airfoil of high-altitude long-endurance fixed-wing aircraft |
CN105351248A (en) * | 2015-12-17 | 2016-02-24 | 新昌县三新空调风机有限公司 | High-performance airfoil for fan |
CN105752314A (en) * | 2016-03-22 | 2016-07-13 | 西北工业大学 | High-altitude, low-speed and high-lift natural laminar flow airfoil |
CN106050739A (en) * | 2016-07-22 | 2016-10-26 | 西安航空制动科技有限公司 | High-performance wing section for cooling fan |
-
2018
- 2018-04-26 CN CN201810382463.7A patent/CN108583847B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0068121B1 (en) * | 1981-07-01 | 1985-08-28 | Dornier Gmbh | Air foil for short take-off and landing aircraft |
CN103482054A (en) * | 2013-08-14 | 2014-01-01 | 西北工业大学 | Low-Reynolds-number wing section matched with full-wing solar unmanned aerial vehicle |
CN104354850A (en) * | 2014-10-30 | 2015-02-18 | 中国人民解放军空军航空大学 | Airfoil profile of airfoil of high-altitude long-endurance fixed-wing aircraft |
CN105351248A (en) * | 2015-12-17 | 2016-02-24 | 新昌县三新空调风机有限公司 | High-performance airfoil for fan |
CN105752314A (en) * | 2016-03-22 | 2016-07-13 | 西北工业大学 | High-altitude, low-speed and high-lift natural laminar flow airfoil |
CN106050739A (en) * | 2016-07-22 | 2016-10-26 | 西安航空制动科技有限公司 | High-performance wing section for cooling fan |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109969381A (en) * | 2019-04-12 | 2019-07-05 | 中国空气动力研究与发展中心计算空气动力研究所 | A kind of low reynolds number airfoil with high ratio of lift over drag with sharp leading edge heavy camber |
CN110386243A (en) * | 2019-07-22 | 2019-10-29 | 北京耐威科技股份有限公司 | A kind of low-speed unmanned aerial vehicle aerofoil profile of high lift-drag ratio |
CN110386243B (en) * | 2019-07-22 | 2024-04-16 | 北京耐威科技股份有限公司 | High lift-drag ratio low-speed unmanned wing type |
CN110435873A (en) * | 2019-08-15 | 2019-11-12 | 西北工业大学 | A kind of half blended wing-body anury formula Unmanned Aerial Vehicle Airfoil race cruised from trim |
CN110435873B (en) * | 2019-08-15 | 2021-04-23 | 西北工业大学 | Cruise self-balancing semi-wing body fusion tailless unmanned aerial vehicle wing type family |
CN111498084A (en) * | 2020-04-15 | 2020-08-07 | 成都飞机工业(集团)有限责任公司 | Low-resistance laminar flow airfoil applied to high-altitude high-speed long-endurance unmanned aerial vehicle |
CN112977815A (en) * | 2021-05-10 | 2021-06-18 | 北京三快在线科技有限公司 | Rotor craft, blade of rotor craft and wing section of blade |
CN112977815B (en) * | 2021-05-10 | 2021-08-27 | 北京三快在线科技有限公司 | Rotor craft, blade of rotor craft and wing section of blade |
CN114194373A (en) * | 2021-12-03 | 2022-03-18 | 中电科芜湖通用航空产业技术研究院有限公司 | Aircraft wing and aircraft |
CN114194373B (en) * | 2021-12-03 | 2024-02-23 | 中电科芜湖通用航空产业技术研究院有限公司 | Aircraft wing and aircraft |
CN114718903A (en) * | 2022-04-19 | 2022-07-08 | 成都航空职业技术学院 | High-performance wing section for heat dissipation axial flow fan |
Also Published As
Publication number | Publication date |
---|---|
CN108583847B (en) | 2020-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108583847A (en) | A kind of low reynolds number high power factor aerofoil profile suitable for long endurance unmanned aircraft | |
CN105752314B (en) | A kind of high-altitude low speed Natural Laminar Flow high-lift airfoil | |
CN103482054B (en) | Low-Reynolds-number wing section matched with full-wing solar unmanned aerial vehicle | |
EP3484764A1 (en) | Airplane wing | |
Rubbert et al. | A general three-dimensional potential-flow method applied to V/STOL aerodynamics | |
CN106050739A (en) | High-performance wing section for cooling fan | |
CN105351248B (en) | A kind of fan high-performance aerofoil profile | |
Prince et al. | Aerodynamic stall suppression on aerofoil sections using passive air-jet vortex generators | |
CN112623254B (en) | Hybrid laminar flow wing air suction energy loss engineering calculation method | |
Morizawa et al. | Effect of Mach number on airfoil characteristics at Reynolds number of 3,000 | |
Eftekhari et al. | Investigation of a cruising fixed wing mini unmanned aerial vehicle performance optimization | |
GölcÁ¼k et al. | Winglet design and analysis for low-altitude solar-powered UAV | |
CN103390109A (en) | Quick prediction method for aerodynamic property | |
Cerón-Muñoz et al. | Experimental analysis of the aerodynamic characteristics adaptive of multi-winglets | |
CN105775159A (en) | Design method for air-blowing ports with function of suppressing separated flow of wings | |
CN106564584B (en) | A kind of unmanned plane | |
Jin et al. | Drag reduction design for a long-endurance electric powered UAV | |
CN113626935A (en) | Design method of transonic crescent wing with high cruising efficiency | |
Woolard | Thin-airfoil theory of an ejector-flapped wing section | |
M Jaffal et al. | Theoretical and Experimental Study of a Forward Swept Wing | |
Yuhara et al. | Drag-Divergence Characteristics of Winglets | |
CN110667820A (en) | Aircraft wing | |
Prasannakumar et al. | Suction System Design for Transonic HLFC Wing | |
CN113665836B (en) | Method for designing steering engine bulge based on sharp Von Karman curve of shock wave intensity of back edge line and the like | |
Sutardi et al. | Flow characteristics around a NACA 4412 airfoil in the ground proximity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |