• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Aerodynamic optimization of a high-lift system with adaptive dropped hinge flap

    2022-12-04 08:03:56QingJIYufeiZHANGHixinCHENJunkeYE
    CHINESE JOURNAL OF AERONAUTICS 2022年11期

    Qing JI, Yufei ZHANG, Hixin CHEN, Junke YE

    a School of Aerospace Engineering, Tsinghua University, Beijing 100084, China

    b 63870 PLA Troops, Huayin 714200, China

    c Shanghai Aircraft Design and Research Institute, Shanghai 200236, China

    KEYWORDS Adaptive dropped hinge flap;Aerodynamic design;High-Lift aerodynamics;Multidisciplinary optimization;Optimization of aircraft design

    Abstract The Adaptive Dropped Hinge Flap(ADHF)is a novel trailing edge high-lift device characterized by the integration of downward deflection spoiler and simple hinge flap, with excellent aerodynamic and mechanism performance. In this paper, aerodynamic optimization design of an ADHF high-lift system is conducted considering the mechanism performance. Shape and settings of both takeoff and landing configurations are optimized and analyzed, with considering the kinematic constraints of ADHF mechanism, and the desired optimization results were obtained after optimization. Sensitivity analysis proves the robustness of the optimal design. Comparison shows that the ADHF design has better comprehensive performance of both mechanism and aerodynamics than the conventional Fowler flap and simple hinge flap design.

    1. Introduction

    The design of an efficient high-lift system is a challenging task in the aerospace industry. Modern civil aircrafts require complex multi-element high-lift systems to increase the lift coefficient to compensate for the low velocity during the takeoff and landing phases.1Takeoff and landing performances for civil airplanes are constrained by strict requirements, such as FAR(Federal Aviation Regulations) Part 25,2to meet the safety demands. It has also been evidenced that even minor improvements of high-lift systems’aerodynamic characteristics can produce significant gains in the aircraft’s overall weight and performance. An increase in the lift coefficient can reduce the approach attitude and shorten the landing gear;an increase in the takeoff lift-to-drag ratio(L/D)can increase the payload or range.3

    High-lift system design in the early years was mainly targeted to improve aerodynamic performance. Based on the analysis of Smith,4a high-lift system that has more elements can generate a larger lift coefficient.Consequently,some complex high-lift configurations, e.g., triple-slotted or doubleslotted Fowler flap systems, were commonly adopted to increase the maximum lift coefficient in the 1960 s to 1990s.5However, the complex flap mechanism also results in heavy weight, high manufacturing cost, and low reliability. It was not surprising that the high-lift system accounts for approximately 10%of the manufacturing cost of typical jet transport,and an alteration from a triple-slotted flap to a double-slotted flap led to an approximately 1/3 flap cost reduction.6

    In recent years, the attention of the high lift design has turned to reduce the complexity, weight, and maintenance costs while keeping an acceptable lift level, which gives a beneficial effect on the overall cost of the aircraft. Rudolph’s report6pointed out that the change from multi-slotted flaps to single-slotted flaps delivered the largest weight and cost reduction of the high-lift system and suggested the singleslotted flap as a recommendation for future development.Both Airbus and Boeing simplified their trailing edge configuration from multi-slotted flaps to a single-slotted flap in their latest airliner design.5

    The high-lift system mechanism is used to support and deploy high-lift movables to the desired position at landing or takeoff phases. There are many different types of trailing edge flap mechanisms used on commercial aircraft,which were reviewed in detail by Rudolph7and van Dam.8Two types of commonly used flap mechanisms are simple hinge flap and fowler flap, as shown in Fig. 16.

    The simple hinge mechanism was commonly adopted in McDonnell-Douglas series aircraft. As shown in Fig. 1(a), a simple hinge flap is characterized by a simple mechanism of a flap connected to a fixed hinge and pure rotational movement. It is mechanically simple and lightweight but has a very restricted rotation trajectory. To obtain more extension of the flap, the pivot point of the hinged flap has to be located far from the wing, which results in a large flap fairing and cruise drag penalty. Moreover, because of the kinematic constraint of the simple hinge mechanism, the optimal configurations of the landing and takeoff may not be realized at the same trajectory.

    Fig. 1 Conventional trailing-edge high-lift devices on commercial aircraft (based on Ref. 6).

    Fowler flaps are the dominant trailing-edge high-lift mechanism in civil aircraft. It is characterized by a nearly translational movement to extend the wing area in takeoff and a rotational movement to increase the wing camber in landing,which provides optimum aerodynamic efficiencies for both takeoff and landing.However,the mechanism is more complex than a simple hinge flap. Three types of mechanisms are commonly used on the Fowler flap: linkage system, track system,and hybrid system; however, all these types of systems tend to be heavy and costly, as well as with lower reliability and maintainability.

    A particularly novel trailing edge device called the Adaptive Dropped Hinge Flap(ADHF)was first introduced on the Airbus A350XWB in 2013.9–10As shown in Fig.2,it is characterized by the integration of a downward deflection spoiler and a simple hinge flap. The ADHF is mechanically simpler than a conventional Fowler flap and requires fewer moving parts,resulting in a considerable weight reduction. The multifunctional spoiler allows regular operation as an airbrake during landing and is able to perform downward deflection to adjust the gap of the single slotted flap. It can significantly increase the flap efficiency without increasing the complexity or the structure weight of the high-lift system. Wang,11Gubsky,12Wang13–14and Liu15further studied the aerodynamic of the dropped spoiler, and pointed out the dropped spoiler will increase the lift coefficient in the linear section, but also decreases the stall Angle of attack.When the spoiler deflection Angle is large, it is easy to separate the upper surface of the spoiler first. Drooped spoilers coupled with a low-setting flap are also used for variable camber technology16, which can reduce the drag coefficient by up to two percent, resulting in considerable fuel economies.17

    Fig. 2 Multifunctional ADHF of A350XWB.10

    The high-lift mechanism not only has an important influence on the aerodynamic performance but also has a significant impact on the weight, cost, maintainability, and reliability of an airplane. With the increasing attention of the comprehensive performance in practical engineering applications, the design methodology of high-lift systems has evolved progressively from an aerodynamic-focused design to a multidisciplinary design in the last few decades.In 1996,Rudolph6–7developed weight and cost modules of different types of highlift devices based on statistical data and recommended a multidisciplinary design methodology that incorporates aerodynamic performance, weight, and cost for future aircraft development. van Dam18proposed an aeromechanical design methodology to evaluate the aerodynamic performance of a high-lift device considering the kinematic constraints of four different trailing edge flap mechanisms. Recently, more and more studies considered the additional kinematic constraints in aerodynamic optimization.19–22Franke23–25improved the aerodynamic performance and reduced the takeoff and landing distance while fulfilling the kinematic constraints, leading to a more realistic design. Pires et al.26–27developed a combined methodology to take the weight, structure, aerodynamic improvements, and cost into account at an early design stage to improve the design quality. Niu et al.28conducted a multipoint aerodynamic optimization of trailing-edge variablecamber by deflecting the flaps and spoiler synchronously,and effectively improved flight efficiency in a wide range of CLby reducing wave drag and form drag. Liu et al.29conducted an aerodynamics/mechanism coupling optimization to a new configuration of link/straight track mechanism coupling with downward deflection of spoiler. The present study constructs an aerodynamic optimization platform considering the mechanism performance to explore the potential of ADHF design.The multi-element airfoil shape,deflection setting, and mechanism parameters of the ADHF are optimized for the best comprehensive performance of both mechanism and aerodynamics. Several aerodynamic and kinematic constraints are included in the optimization. Parametric analysis is conducted for the optimum configuration to test the sensitivity of the ADHF design. A comparison study between the ADHF optimization and the conventional design of Fowler flaps and simple hinge flaps was conducted to test the advantages of the ADHF configuration.

    2. Numerical method and validation

    A Reynolds-Averaged Navier-Stokes (RANS) solver is used for Computational Fluid Dynamics (CFD) simulation. The Spalart-Allmaras (S-A) turbulence model is adopted for the Reynolds stress closure. The 30P30N multi-element airfoil is adopted as the validation case. It has plenty of experimental data from the Low Turbulence Pressure Tunnel (LTPT) at NASA Langley Research Center in a cooperative effort with McDonnell-Douglas Aerospace.30Three sets of grids of the 30P30N airfoil with different grid numbers are adopted to test the grid convergence.The grid information is listed in Table 1,in which Nxis the number of grid points along the circumference of the solid surface. The computational domain size is[-50c, 50c] in the x-direction and [-40c, 40c] in the ydirection to fulfill the requirement of the far-field boundary condition. c is the chord length of the retracted airfoil. The mesh is refined in close proximity to the slat and the flap, as well as the upper surface of the main element, to ensure that the wakes can be well captured. The grid spacing of the first grid layer normal to the wall is 5×10-6c to ensure Δy+<1.The medium grid is shown in Fig. 3. The flow conditions are Ma∞=0.2 and Re=5×106.

    The lift curves of the computational results are compared with the experimental data, as shown in Fig.4.Both the medium and fine grids agree well with the experimental data,except that the stall angle is slightly later. The relative error of the maximum lift coefficient is less than 2.2%. Fig. 5 demonstrates that the Cpand Cfdistribution predicted at α=21°is in good agreement with the experimental data.This shows that the present computational method is reliable and the grid convergence is achieved. Consequently, the present CFD solver with the S-A turbulence model and the medium grid are adopted in the following study.

    3. Optimization of ADHF airfoil

    The design of the present ADHF airfoil is based on the retracted geometry of the 30P30N airfoil. The maximum thickness-to-chord ratio t/c of the airfoil is 0.106. The ADHF high-lift configuration is a simple three-element airfoil that consists of a slat, a main wing with a downward deflectable spoiler, and a single-slotted flap. The optimization of ADHF is to find the optimal shape,deflection setting,and mechanism parameters of an ADHF for the best comprehensive performance of both mechanism and aerodynamic.The aerodynamic and kinematic constraints are also included in the optimization. This approach is general and can also be extended to three-dimensional configurations. The present optimization is performed under both takeoff and landing conditions of freestream Mach number Ma∞= 0.2 and Reynolds number Re=9×106based on the retracted chord length.

    3.1. Parametrization and optimization settings

    3.1.1. Aerodynamic and kinematic constraints

    To obtain a realistic design, the optimization process must handle the aerodynamic constraints and the kinematic constraints of the ADHF mechanism. The following constraints are considered in the optimization:

    (1) The chord lengths of the slat and flap are limited by the front spar and rear spar of the wing to ensure the stiffness and fuel tank capacity, as shown in Fig. 6. For a conventional commercial transport aircraft,the position of the wing front spar is between 0.15c-0.20c,and the position of the wing rear spar is between 0.65c -0.75c.12–13For the ADHF, the trailing edge high-lift device of the Airbus A350XWB is realized by a small flap chord of 0.19c,which is smaller than a conventional high-lift system. It is beneficial for wing stiffness and weight reduction.In this study,the wing front spar position is located at x=0.145c, and the wing rear spar is located at x=0.75c.The chord of the flap is a variable to be optimized between 0.15c -0.25c.

    Table 1 Grid information of 30P30N airfoil.

    Fig. 3 Medium grid of 30P30N multielement airfoil.

    (2) The spoiler of ADHF is downward deflectable following the flap. Due to the kinematic constrains of the mechanism,the spoiler is connected to a fixed hinge and moved in a pure rotational trajectory, and should not interfere with the flap.The geometry interference check is automatized in the optimization process.The spoiler chord and the flap chord together determine the clean overlap of the flap, which will also be further determined by optimization design.

    (3) Mechanism depth of flap is also an important constraint, because the mechanism depth determines the flap fairing size and weight, and also fairing drag. For commercial transport aircraft, the typical fairing depth of the Fowler flap is 0.09c-0.16c.26The ADHF design of A350XWB and B787 further reduces the fairing size and achieves excellent mechanism performance. In this paper, the hinge depth Yhof ADHF is constrained to be less than 0.075c, and also be taken as the optimization objective for the optimization design. In addition,the ADHF mechanism also brings strong kinematic constraints to the design of the high-lift device,both the slat and the flap moves along the circular trajectories around the hinge points.The slat and flap positions of the takeoff and landing configurations have to be compatible in the same circular trajectories.

    (4) To avoid a sudden drop in the lift coefficient in the linear section, which often occurs due to flow separation occurring on the flap designs at large deflection (≥30°)at lower angle of attack,31the lift coefficient CLat α=4° of the landing configuration is constrained to be larger than 2.7 as aerodynamic constraint.

    Fig. 4 Variation of experimental and computational CL with α.

    Fig. 5 Variation of experimental and computational Cp and Cf at α=21°.

    3.1.2. Design parameters

    The shape of both the slat and flap and the settings at both takeoff and landing are used as design parameters. Twentysix design variables are defined in the optimization process,as shown in Fig. 7 and listed in Table 2:

    (1) Slat shape parameters: the shape of the slat is cut from the baseline clean airfoil by a cubic spline curve. The curve forms both the leading edge of the main element and the cove shape of the slat. The shape of the slat is controlled by 5 control points and has 6 variables. As shown in Fig. 7(a), Ps0and Ps4control the chord of the slat. Ps1is the leading-edge point of the mainelement.Ps2and Ps3control the shape of the nose shape of the main element.

    Fig. 6 General constrains on high-lift design.

    (2) Flap shape parameters: the geometry of the flap is controlled by 6 control points, including 10 variables. As shown in Fig.7(b),Pf5controls the chord of the spoiler.Pf1determines the chord of the flap. Pf3dominates the maximum thickness and maximum thickness position of the flap. Pf2controls the leading edge shape of the flap.Pf4is used to adjust the shape of flap upper surface.

    (3) Takeoff and landing settings: including gaps s, overlaps d, and deflection angles δ of the slat, flap, and spoiler.Because the deflection trajectories of the takeoff and the landing are the same,the takeoff and landing parameters are dependent.In this paper,the hinge locations of the slat and flap are determined by the landing configuration.Once the gaps,overlaps,and deflection angles of the slat or flap of the landing configuration are provided,the hinge locations can be computed by rigid body rotation,as shown in Fig.7.Then,the takeoff configuration can be generated only by the deflection angles of the slat,flap,and spoiler.Consequently,the design parameters of the landing configuration have 7 variables,and the takeoff configuration has 3 variables, as listed in Table 2.

    3.1.3. Optimization objectives

    The present optimization aims to increase the aerodynamic performance of both takeoff and landing conditions and the mechanism performance of the high-lift device. Tri-objective optimizations with several engineering constraints are adopted in the current study to cover the following design requirements:

    (1) The first objective is to increase the maximum lift coefficient CL,maxin the landing configuration since the increase in CL,maxcan increase the maximum landing weight and decrease the approach speed and landing range.32

    (2) The second objective is to improve the lift-to-drag ratio L/D under the takeoff condition(CL=2.5).Because the takeoff climb gradient is dependent on the lift-to-drag ratio, increasing L/D at a given lift coefficient can improve the climb performance during the second segment climb, which is a key component in measuring the total range of an aircraft.

    (3) The third objective is the hinge depth Yh, which represents the distance of the flap hinge to the wing. In this study, the hinge depth Yhis chosen as the key indicator for the mechanism performance of ADHF, because a lower hinge depth Yhleads to a smaller support structure and a lighter mechanism weight, a lower manufacturing cost, a smaller and lighter flap fairing, and also lower fairing drag. The hinge depth Yhis computed based on the settings of the flap. It is worth noting that the hinge depth of the slat is not optimized in this study because the slat is usually driven by an internal gear rack, and no external faring of the mechanism is required.

    Fig. 7 Sketch of design parameters.

    3.2. Optimization modules

    3.2.1. Geometry generation and grid automation

    An in-house-developed high-lift geometry generator is used to generate new configurations. The geometry generator takes design parameters as input and generates high-lift geometry quickly and robustly. The geometry constraints, such as the trajectory constraint of the landing and takeoff and the gapand overlap feasibility, are implemented in the geometry generator.

    Table 2 Range of design parameters.

    The computational grid is automatically generated using a grid deformation script. The baseline grid is a multiblock structured grid based on the medium grid of the 30P30N three-element airfoil. The grid deformation script contains a sequence of operations to readapt the baseline block topology around the new configuration. By carefully designing the grid topology and grid deformation operations,the deformed mesh preserves a good grid quality even at the upper or lower limits of the design variables.

    3.2.2. Optimization algorithm

    Because there are multiple optimization variables and multiple optimization objectives, and the optimization can be highly nonlinear and difficult to converge. Genetic-type algorithms are often used in high-lift design to achieve a globally optimized solution in either two-dimensional33or threedimensional problems.34

    An in-house developed hybrid surrogate-aided differential evolution optimization algorithm called HSADE35is chosen to find the global optimal solution in the current study. The parallel combination of differential evolution and the radial basis function response surface improves the elitist performance and guarantees population diversity at the same time.36The optimization efficiency and high robustness of the HSADE algorithm have been verified by our previous studies.28,37–39

    3.2.3. Optimization process

    The optimization flow chart is illustrated in Fig. 8. An initial set of design variables representing the configuration designed manually is provided to the optimizer.A vector of design variables represents each design candidate. Then, the vector is transferred to the geometry generation module. The landing configuration and takeoff configuration are optimized in a parallel process. The geometry generation module translates the design parameters into landing and takeoff configurations with respect to the kinematic law of the ADHF mechanism. The grid automation module separately generates the grid of landing and takeoff automatically.The CFD analysis module computes the flow fields and evaluates the aerodynamic performance of the configurations.Constant lift coefficient calculation of the takeoff configuration is performed by automatically iterating the angle of attack.Several discrete angles near the stall angle of attack are calculated and compared to determine the maximum lift coefficient of the landing configuration.The aerodynamic performances of landing and takeoff configurations and the hinge depth of the ADHF form the objective functions. If an optimal solution is not achieved, the function values of the objectives are fed back to the HASDE optimizer to generate new populations for the next design generation.The process is repeated until the Pareto front of the optimal solutions is converged.

    Fig. 8 Flow chart of optimization procedure.

    3.3. Optimization results

    The optimization was performed on an Intel 32-core 2.4 GHz workstation. Thirty individuals are computed in each generation.It takes approximately 1 core hour to evaluate an individual. Sixty generations and 1800 individuals are computed for each optimization.

    3.3.1. Selection of optimal solutions

    The tri-objective optimization leads to a three-dimensional Pareto front. Three projections of the three-dimensional (3D)Pareto front are shown in Fig. 9 for better display. The gray dots represent the valid individuals generated in the optimization. The color dots present the individuals of the third objective from the 3D Pareto front of the optimization.As shown in Fig.9(a),the color of dots marks out the Yhvalue of the optimization. As shown in the figure, a sharp Pareto front or an‘‘L-shaped” Pareto front (refer to Refs. 26,16) is obtained in the L/D vs CL,maxplane, indicating that the two objectives are almost independent. It also indicates that the takeoff and landing settings have a more significant influence than the common parameters of the shape on the aerodynamic objectives.A unique optimal solution of the aerodynamic objectives could be obtained at the corner of the ‘‘L-shaped” Pareto front.

    Fig. 9 Individual distributions of tri-objective optimization.

    The smooth Pareto front in Fig.9(b)and 9(c)indicates that both the aerodynamic objectives of L/D and CL,maxare influenced by the mechanism objective of Yh. There is a trade-off between the mechanism size and the aerodynamic performance.

    Three individuals on the Pareto front are chosen for comparison to decide the final optimal individual. The first one has the shortest Yhnamed as ‘‘Yhbest”. The second one is a trade-off of the three objectives named the ‘‘Trade-off”. The third one has the best aerodynamic performance named‘‘Aero best”.

    The landing and takeoff geometries of the three designs,as well as the original multi-element airfoil, are shown in Fig. 10 and Fig. 11. The design parameters are also shown in Table 3 and Table 4.The landing configuration and the takeoff configuration share the same shape and hinge point. The leading edge shapes of the main element and flap of the 3 alternatives almost coincide and show a plump leading edge of both the main element and flap and a shorter spoiler and longer flap chord than the original configuration,which indicates the convergence of the shape parameters.

    The candidate landing configurations have larger deflection angles of the flap and spoiler,and smaller Fowler motions than the original configuration. By comparing the three candidate designs,one can see that a larger CL,maxof the landing configuration, corresponding to a larger slat deflection angle and smaller spoiler and flap deflection angles. And a smaller flap deflection angle leads to a longer hinge depth Yh.

    The candidate Takeoff configurations have larger deflection angles of the flap and spoiler,a smaller angle and nearly sealed slat, and smaller Fowler motions than the original configuration. By contrast, the three takeoff configurations are quite similar.This demonstrates that a larger L/D of the takeoff configuration corresponds to a smaller slat deflection angle and bigger Fowler motions. The nearly sealed slat leadings to better L/D performance and reduces the slat noise,9–10the larger Fowler motion leads to longer hinge depth Yh.

    To further assess the performance at different angles of attack, the lift curves and lift-to-drag curves of the three designs are computed and shown in Fig. 12.

    The‘‘Aero best”configuration has the largest CL,maxin the landing configuration, but the lift curve is lower in the linear section. The ‘‘Aero best” has the best L/D when the takeoff CLis above 2.2.

    The‘‘Yhbest”configuration has the highest CLin the linear section, but it falls suddenly below α=3°. The streamline at α=2° is shown in the cyan frame in Fig. 12(a). The flow detached at the upper surface of the flap.

    Compared with the ‘‘Yhbest” and ‘‘Trade-off” configurations,the‘‘Trade-off”configuration has a moderate maximum lift coefficient CL,max,and lift-to-drag ratio L/D.It also has the largest stall angle of attack (α=24°). Therefore, the ‘‘Tradeoff” configuration is selected as the optimal individual of this optimization. In the following study, this configuration is denoted as the ‘‘ADHF opt” for convenience.

    3.3.2. Performance and flow analysis of landing configuration

    Compared with the original configuration, the performance improvement of the‘‘ADHF opt”is the result of the combined action of the optimized shape and settings.To separately study the contributions of the shape and the settings, a new landing configuration is generated by the original airfoil shape,and the settings of the selected optimal configuration ‘‘ADHF opt”.This configuration is denoted as the ‘‘opt setting” in this section. The lift curves of the original airfoil, the ‘‘ADHF opt”configuration, and the ‘‘opt setting” configuration are shown in Fig. 13.

    (1) Contribution of settings.

    By comparing three CLcurves, it can be seen that the settings have the most important influence on the landing configuration. The lift coefficient of the ‘‘opt setting” improves considerably by approximately 0.6 at a low angle of attack,which is beneficial for reducing the landing gear length.3With the increase in the angle of attack, the slope of the lift curve decreases. The stall lift coefficient of the ‘‘opt setting” is slightly lower than that of the ‘‘original” configuration.

    The pressure coefficients Cpof the three landing configurations at α=4°are shown in Fig.14.The larger flap deflection angle of the optimized solution tends to increase the loading on the main element and the slat by enhancing the upwash of the upstream element.

    The drooped spoiler increases the loading on the main element by increasing the aft camber of the main element and increases the loading on the slat by enhanced upwash. It also reduces the leading-edge suction peak on the large deflected flap by the downwash, thus reducing the pressure gradient and delaying separation, increasing the lift coefficient.

    As the angle of attack increases,the load on the main wing increases rapidly, which increases the downwash to the flap and results in the decrease of the local angle of attack of the flap. The drooped spoiler further reduced the local angle of attack of the flap,which reduces the adverse pressure gradient and delays the trends of flap separation. The reduced loading on the flap also weakens the upwash to the main wing in turn,resulting in the decrease of lift coefficient increment and the decrease in the lift curve slope, as shown in Fig. 13.

    Fig. 10 Sketches of alternative landing configurations.

    Fig. 11 Sketches of alternative takeoff configurations.

    Table 3 Parameters of alternative landing configurations.

    Table 4 Parameters of alternative takeoff configurations.

    (2) Contribution of shapes.

    The difference between the‘‘ADHF opt”and the‘‘Opt setting” comes from the different nose shapes of the two configurations. This shows that the optimized shape is also beneficial for the improvement of the lift coefficient. The stall lift coefficient of the optimized configuration is higher than that of the ‘‘original” configuration.

    However, the nose shape has a limited influence, mainly affecting the leading edge suction peak, as shown in Fig. 14.The larger curvature shape just behind the slat and flap slot accelerates the slot jet flow and increases the leading-edge suction peak of the main element and the flap,as shown in Figs.15 and 16.When the attack angle increases,the jet flow speeds up,and the leading-edge suction peak increases simultaneously.

    Fig.12 Aerodynamic performance of alternative configurations.Consequently,the difference in CLcurves of the‘‘ADHF opt”and ‘‘opt setting” is also increasing, as shown in Fig. 13.

    3.3.3. Performance and flow analysis of takeoff configuration

    The takeoff configuration and the L/D curves of the‘‘original,‘‘the‘‘ADHF opt”and the‘‘opt setting”are shown in Fig.17.

    (1) Contribution of settings

    The improvement between the ‘‘original” and the ‘‘opt setting” comes from the different settings of the two configurations. The larger flap of the ‘‘opt setting” increases the aft camber of the high-lift system, and the spoiler deflection increases the downwash to the flap, delaying the trends of trailing-edge separation,40–41which reduces the angle of attack for the takeoff angle of attack from 10.3° to 5.1°. This is beneficial for the reduction of the drag coefficient and the improvement of the L/D.

    (2) Contribution of shapes

    Fig. 13 Landing configuration and lift curves comparison.

    The difference between the ‘‘ADHF opt” and the ‘‘opt setting” comes from the different nose shapes and the different trajectories of the two configurations. Fig. 17(a) shows that the ‘‘ADHF opt” configuration has an almost sealed slat gap in the takeoff configurations. The trajectory of the slat gap and deflection angle is shown in Fig. 18. The gap size of the‘‘original”and‘‘opt setting”configurations increases homogeneously as the slat angle increases, while ‘‘ADHF opt” maintains an almost sealed slat gap when the slat angle is less than 10° but increases rapidly as the slat angle continues to increase. This results in enough gap for landing but a sealed slat gap for takeoff. Fig. 19 illustrates the Cpcontour of the‘‘ADHF opt” and the ‘‘opt setting” configurations at α=4°.The greater suction at the leading edge of the ‘‘ADHF opt”slat, and larger pressure in the sealed slat cove improves the CLand reduces the pressure drag, and the sealed slat also reduces the friction drag caused by jet flow, as shown in Fig. 20. Consequently, it improves the L/D of the ‘‘ADHF opt”configuration by 11.87%than that of the‘‘original”configuration. Moreover, the non-slotted slat at takeoff also benefits the slat noise reduction of modern civil aircraft.42

    3.3.4. Sensitivity study

    Fig. 14 Cp landing configuration comparison.

    Fig. 15 Slot flow of ‘‘ADHF opt” landing configuration(α=24°).

    A practical design should have both good performance and robustness. It is worth studying the influence of each parameter on the performance of the optimized design.When an optimum design is obtained, the shape and the mechanism of the high-lift system would not change. The hinge depth and the trajectory of the slat and flap would not change either. However, the gaps, overlaps, and deflection angles of the slat and flap might be slightly changed during flight because of elastic deformation. The disturbance of these parameters may interfere with aerodynamic performance. Therefore, the aerodynamic performance under the disturbance of slat and flap settings is analyzed in this section. Because of the kinematic constraints of the ADHF mechanism, the changes in the gap and overlap are correlated with the deflection angle of the slat and flap.

    (1) Landing configuration

    The optimal slat deflection angle for the landing configuration is 25.99°. The landing performance of 21 different slat deflection angles from 24.00°to 28.00°are assessed and shown in Fig. 21. The baseline value is marked out in the figure. The CLat α=4° decreases gently as the slat angle increases. The optimal configuration has the maximum landing CLat α=25.99°, and the CLdecreases slightly near the baseline value.The lift coefficient drops rapidly when the slat deflection angle is too large or too small. The streamlines of the critical deflection angles are also shown in the figure. When the slat angle is larger than 27.2°,the increased slat deflection weakens the jet flow, which leads to massive flow separation above the spoiler,as marked by the blue frame in Fig.21.If the slat angle continues to increase,the flow separation expands to the main wing and destroys the circulation of the airfoil. By contrast, if the slat angle is smaller than 24.6°, the decrease in the slat angle increases the suction peak on the leading edge, which leads to the presence of the separation bubble on the slat, as marked by the red frame in Fig. 21. The formation and shedding of the separation bubble leads to the instability of the flow.

    For the flap and spoiler deflections, the optimal values are δf=40.63°and δsp=12.56°.441 configurations of δfranging from 38.6° to 42.6° and δspranging from 10.6° to 14.6° are assessed and shown in Fig.22.For α=24°,CLincreases when the flap deflection increases and the spoiler deflection decreases. Interestingly, CLremains the same value when the flap and spoiler deflections increase or decrease simultaneously.When α=4°,there is a sudden drop in CLwhen the flap deflection increases and the spoiler deflection decreases.

    Fig. 16 Slot flow of ‘‘Opt setting” landing configuration(α=24°).

    Changing the deflection angles of the flap and spoiler also changes the gap and overlap of the flap.To further investigate the effect of the design parameters, the result of Fig. 22 is redrawn against the changing gap and overlap, as shown in Fig.23.This reveals that the gap is the most important parameter on the lift coefficient CL. When the flap deflection increases and the spoiler deflection decreases simultaneously,the flap gap size increases. The increased gap leads to an increasing CLat 24° but a sudden drop in CLat 4°. Because the increasing flap deflection and the decreasing spoiler deflection increase the suction peak of the flap, the wider flap gap increases the mass flow in the gap and enhances the upwash of the main wing, resulting in a higher CLat 24°. However,a high suction peak on the flap at 4°increases the adverse pressure gradient and tends to flow separation. When the gap size is smaller than 0.0125c, CLincreases gently as the gap size increases. When the gap size exceeds the size of 0.0125c, the CLfalls suddenly because of flap separation, as shown in Fig.23(b).The optimal configuration is at the edge of the flap separation, which proves the convergence of the optimization.For an ADHF high-lift system, the gap size of the flap can be easily controlled by the flexible spoiler. Reducing the gap size improves the robustness at a low angle of attack.However,the maximum lift coefficient of the landing configuration decreases as the penalty increases.

    (2) Takeoff configuration

    Fig. 17 Takeoff configuration and lift curve comparison.

    Fig. 18 Trajectory of slat gap and deflection angle.

    The slot of the flap is determined by both the flap and the drooped spoiler. Therefore, the flap deflection and the spoiler deflection are assessed simultaneously. The coupling effects of the two parameters are also studied. The optimal values of the flap deflection angle and the spoiler deflection angle are δf=29.42°and δsp=7.75°.441 configurations of δfranging from 27.4° to 31.4° and δspranging from 5.8° to 9.8° are computed and shown in Fig. 24. This shows that the CLand the L/D of the takeoff configuration are sensitive to the flap angle. The increase in flap deflection leads to an increase in CLbut a decrease in L/D. The spoiler deflection seems optimum in the baseline configuration.Both CLand L/D decrease as the spoiler deflection deviates from the baseline value.

    Fig. 19 Cp contour map of ‘‘ADHF opt” and ‘‘Opt setting”(α=4°).

    Fig. 20 Cf of ‘‘ADHF opt” and ‘‘Opt setting” (α=4°).

    Fig. 21 Landing configurations with different δs.

    The contours in Fig.24 also demonstrate a strong coupling effect between the two parameters acting on the CLand L/D.This coupling effect decreases the lift coefficient CLif the flap deflection decreases and the spoiler deflection increases simultaneously.This is expected since this movement closes the flap gap, reducing the flux of jet flow and reducing the suction on the flap’s upper side. In contrast, the coupling effect decreases the L/D when the flap deflection increases and the spoiler deflection decreases simultaneously. The opposite tendency indicates that a compromise has to be made between the lift coefficient CLand the lift-to-drag ratio L/D.

    3.4. Comparison with conventional designs

    Two more optimizations of the conventional configuration of the Fowler flap and simple hinge flap are carried out to evaluate the efficiency of the ADHF high-lift device.Both optimizations follow a conventional optimization process, as shown in Fig. 25. A landing configuration optimization of maximizing CL,maxruns first to obtain the best shape and settings. Then,a takeoff configuration optimization runs to maximize the L/D, which has the same shape and kinematic constraint provided by the landing optimization.Both conventional configurations of the Fowler flap and simple hinge flap have the same chord length of the slat and flap and have no spoiler deflection.

    Two optimal configurations of the Fowler flap and simple hinge flap are obtained from the conventional optimization process. Both configurations are compared with the optimal design of the ADHF optimization. The geometries of three designs are shown in Fig. 26, and the parameters are shown in Table 5 and Table 6.

    The simple hinge flap has the same landing configuration as the Fowler flap configuration because of the same objectives and the same optimization process. The conventional design has an almost coincidentally plump shape as the ADHF design to obtain a higher lift. The landing configuration of the conventional design has a smaller flap angle and a larger Fowler motion than the ADHF design.

    For the Fowler flap, the takeoff configuration is optimized separately from the landing configuration. Both the best landing and takeoff configurations have large Fowler motions but a rather different flap deflection angles. A mechanism of at least 2 degrees of freedom is needed to reach the different flap locations. The most commonly used Fowler flap mechanisms are linkage systems, track systems, and hybrid systems.6

    Fig. 22 Landing configurations with different δsp and δf.

    Fig. 23 Landing performance with respect to sf and df.

    Fig. 24 Takeoff configurations with different δsp and δf.

    Fig. 25 Flow chart of conventional optimization process.

    For the simple hinge flap, the flap is connected to a fixed hinge and moved in a pure rotational trajectory.When the best landing configuration is obtained, the hinge point and the trajectory of the flap are determined.The takeoff configuration is a result of the intermediate configuration.The desired optimal configurations in landing and takeoff cannot be met at the same time.The simple hinge flap is mechanically simple,lightweight, and reliable; however, the hinge point of the flap is quite far from the wing, resulting in a large flap fairing and cruise drag penalty.By contrast, the ADHF has a relatively short extension motion. The downwash of the drooped spoiler ensures that the flap deflects to a larger angle without flow separation.This helps to reduce the hinge bar length Yhby 50.6% compared with the simple hinge system. The size and weight of the trailing edge mechanisms can be reduced considerably to reduce the fairing size and the fairing drag, which makes the ADHF an efficient design.

    Table 5 Parameters of 3 types of high-lift designs at landing configurations.

    Table 6 Parameters of 3 types of high-lift designs at takeoff configurations.

    3.4.1. Comparison of aerodynamic performance

    The aerodynamic performance of three designs is compared in Fig. 27.

    Fig. 26 Geometries of the 3 types of high-lift designs.

    Fig. 27 Aerodynamic performance of 3-type high-lift design.

    Table 7 Estimated weights of trailing-edge high-lift mechanisms (lb/ft2 of stowed flap area).

    For landing configuration,the drooped spoiler and the larger flap deflection angle increase the overall camber of the ADHF, resulting in an increment of approximately 0.5 for the landing CLat small angles of attack. The CL,maxof the ADHF landing configuration is slightly decreased by 2.4%,but it still maintains a higher CLbefore 20.5°. The lift coefficient increment at the linear section results in a reduction in the approach attitude. This might lead to a shortening of the landing gear due to lower ground clearance requirements and thus saving weight.

    For takeoff configuration, because of the kinematic constraint of the trailing edge mechanism, the L/D of the simple hinge is approximately 2.6%lower than that of the Fowler flap for a given takeoff CLof 2.5. The L/D of ADHF is considerably increased by 10.3%, and the takeoff angle of attack is reduced from 7.7° to 4.1°, compared with that of the Fowler flap design. The increase in takeoff L/D leads to an increase in the payload or an increase in range. The reduction of the takeoff angle of attack is also conducive to reducing the demand of the plane’s taildownangle and improving the comfort of passengers during takeoff.

    3.4.2. Comparison of mechanism performance

    The weight of different trailing edge mechanisms is estimated by empirical equations based on statistical data.Pepper et al.43divided the weight of a high-lift mechanism into 3 parts: supports & linkages, actuation systems, and fairings, and proposed a fairly accurate prediction method of the weights of the trailing edge high-lift mechanisms. In this paper, the trailing edge weights are evaluated based on Pepper’s method.43The weights estimation of the 4 different types of trailing edge high-lift mechanisms for a representative single slotted flap is presented in Table 7. The support and fairing weight of the ADHF is estimated as half of the simple hinge flap. A reference airplane with a maximum gross weight of 113400 kg(250000 lb) is used for evaluation. For a trailing-edge flap stowed area of 30.67 m2(330 ft2), the evaluated weights are estimated in Table 8.

    Table 8 Estimated weights for trailing-edge high-lift mechanisms of 113400 kg (250000 lb) gross weight airplane.

    It shows that the hooked track and link/track systems for the conventional Fowler mechanisms are quite heavy.The simple hinge mechanism is 36.8%lighter than the link/track mechanism due to its simplicity. The ADHF is even lighter because of its shorter hinge. The weight of the ADHF trailing-edge mechanism is reduced by 315.9 kg compared to that of the link/track mechanism.Since the reduction of structural weight cycles into the airplane takeoff gross weight with a sensitivity factor of 1.5 to 2.5,6the 315.9 kg weight reduction could make the airplane 473.9 kg to 789.8 kg lighter at the same payload or would allow a significant increase in the payload or fuel. This weight estimation demonstrates that the ADHF has an outstanding mechanism performance.

    4. Conclusions

    In this paper, an aerodynamic optimization of ADHF considering the mechanism performance is conducted based on the retracted clean airfoil of the 30P30N airfoil.The optimal takeoff and landing configurations with good comprehensive performance of mechanism and aerodynamics are obtained through optimization.

    (1) The optimization demonstrates the downward deflection of the spoiler,and the optimized settings lead to lift coefficient increments of approximately 0.6 over most of the lift curve.The larger curvature shape behind the slat and flap slots accelerates the slot jet flow and increases the leading-edge suction peak of the main element and the flap.The sealed slat of the ADHF takeoff configuration improves the takeoff L/D by 11.87%.

    (2) A parametric analysis was conducted to test the sensitivity of the ADHF design. This reveals that the gap is the most important parameter on the landing lift coefficient CL. The performance of the landing configuration is robust when keeping a flap gap size smaller than 0.0125c. The robustness of the takeoff configuration is also validated.

    (3) Comparisons between the ADHF optimization and the conventional designs of Fowler flaps and simple hinge flaps shows that the ADHF design has excellent aerodynamic performance of 10.3% higher in takeoff L/D and approximately 0.5 higher in landing CLover most of the lift curve than that of the Fowler flap design, With mechanism of 50.6% shorter and 34.2% lighter than that of simple hinge flap design. Demonstrating that the ADHF design has outstanding comprehensive performance.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Nos. 11872230, 91852108, 91952302,92052203),and the Aeronautical Science Foundation of China(No. 2020Z006058002).

    又黄又爽又免费观看的视频| 日韩熟女老妇一区二区性免费视频| 黄色成人免费大全| av一本久久久久| 在线av久久热| 日本vs欧美在线观看视频| e午夜精品久久久久久久| 久久久久视频综合| 国产高清国产精品国产三级| 男女高潮啪啪啪动态图| 国产精品久久久久久人妻精品电影| 少妇被粗大的猛进出69影院| 亚洲全国av大片| www日本在线高清视频| 国产精品一区二区精品视频观看| 在线观看www视频免费| 好男人电影高清在线观看| 精品久久久久久,| avwww免费| 少妇猛男粗大的猛烈进出视频| xxx96com| 欧美精品av麻豆av| 日韩欧美在线二视频 | 自拍欧美九色日韩亚洲蝌蚪91| 中文欧美无线码| 日本一区二区免费在线视频| 狂野欧美激情性xxxx| 久久国产精品人妻蜜桃| 91精品国产国语对白视频| 亚洲欧美激情综合另类| a级毛片在线看网站| 俄罗斯特黄特色一大片| 日本精品一区二区三区蜜桃| 啦啦啦免费观看视频1| 午夜福利欧美成人| 亚洲视频免费观看视频| av有码第一页| 久久中文字幕人妻熟女| 另类亚洲欧美激情| 一a级毛片在线观看| 人人妻,人人澡人人爽秒播| 欧美成人午夜精品| 欧美日韩av久久| 女同久久另类99精品国产91| 色播在线永久视频| 亚洲午夜理论影院| 国产成人av激情在线播放| 欧美最黄视频在线播放免费 | 一级a爱片免费观看的视频| 欧美午夜高清在线| 亚洲第一av免费看| 另类亚洲欧美激情| 大型av网站在线播放| 后天国语完整版免费观看| 女警被强在线播放| 黄频高清免费视频| 天堂√8在线中文| 在线免费观看的www视频| 久热这里只有精品99| 久久狼人影院| 丰满的人妻完整版| 亚洲男人天堂网一区| 巨乳人妻的诱惑在线观看| 日韩成人在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 国产高清视频在线播放一区| 侵犯人妻中文字幕一二三四区| 伊人久久大香线蕉亚洲五| 国产精品影院久久| netflix在线观看网站| 欧美成人免费av一区二区三区 | 欧美日韩福利视频一区二区| 在线观看免费午夜福利视频| 亚洲成a人片在线一区二区| 高清毛片免费观看视频网站 | 亚洲专区国产一区二区| 成年人午夜在线观看视频| 日韩中文字幕欧美一区二区| 老熟妇仑乱视频hdxx| 中文字幕精品免费在线观看视频| av国产精品久久久久影院| 校园春色视频在线观看| 伦理电影免费视频| 男男h啪啪无遮挡| 9热在线视频观看99| 欧美日韩国产mv在线观看视频| 日日夜夜操网爽| 欧美久久黑人一区二区| 午夜91福利影院| 免费日韩欧美在线观看| 黄色成人免费大全| 99久久人妻综合| 无遮挡黄片免费观看| 搡老熟女国产l中国老女人| 亚洲情色 制服丝袜| 久久久久久久久免费视频了| 亚洲 欧美一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲国产欧美日韩在线播放| 国产精品98久久久久久宅男小说| 日日夜夜操网爽| www.精华液| 免费久久久久久久精品成人欧美视频| 国产精品亚洲av一区麻豆| 我的亚洲天堂| 国产人伦9x9x在线观看| 别揉我奶头~嗯~啊~动态视频| 最新的欧美精品一区二区| 欧美精品人与动牲交sv欧美| 国产av精品麻豆| 亚洲五月天丁香| av福利片在线| 少妇裸体淫交视频免费看高清 | av网站免费在线观看视频| 精品国产一区二区三区久久久樱花| a级片在线免费高清观看视频| 国产av一区二区精品久久| 巨乳人妻的诱惑在线观看| 99re在线观看精品视频| 精品一品国产午夜福利视频| 亚洲熟妇熟女久久| 18禁美女被吸乳视频| 99在线人妻在线中文字幕 | 久久久久国产一级毛片高清牌| 久久草成人影院| 精品一区二区三区四区五区乱码| 免费av中文字幕在线| 变态另类成人亚洲欧美熟女 | 免费看a级黄色片| 国产精品二区激情视频| 热99久久久久精品小说推荐| 成人国产一区最新在线观看| 亚洲熟女毛片儿| 一进一出抽搐gif免费好疼 | 每晚都被弄得嗷嗷叫到高潮| 亚洲熟妇熟女久久| 老鸭窝网址在线观看| 国产成人精品无人区| 99精品在免费线老司机午夜| 国产成人免费无遮挡视频| 亚洲av片天天在线观看| 国产人伦9x9x在线观看| 精品国内亚洲2022精品成人 | 91九色精品人成在线观看| 国产男靠女视频免费网站| 午夜亚洲福利在线播放| 99久久综合精品五月天人人| 精品久久久久久久久久免费视频 | 日本黄色日本黄色录像| 丰满的人妻完整版| 91av网站免费观看| 1024视频免费在线观看| 午夜老司机福利片| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看 | 9色porny在线观看| 50天的宝宝边吃奶边哭怎么回事| 在线播放国产精品三级| 精品一区二区三区av网在线观看| 免费一级毛片在线播放高清视频 | 黄网站色视频无遮挡免费观看| 精品卡一卡二卡四卡免费| 女性被躁到高潮视频| 亚洲五月色婷婷综合| 亚洲少妇的诱惑av| 日本黄色日本黄色录像| 亚洲中文av在线| 日韩欧美免费精品| 亚洲av第一区精品v没综合| 麻豆国产av国片精品| x7x7x7水蜜桃| 欧美大码av| 757午夜福利合集在线观看| 黄色 视频免费看| 亚洲成人免费av在线播放| 亚洲av成人av| 色精品久久人妻99蜜桃| 中亚洲国语对白在线视频| 国产熟女午夜一区二区三区| 在线av久久热| 一区二区三区国产精品乱码| 国产一卡二卡三卡精品| 国产精品亚洲av一区麻豆| 欧美成人午夜精品| 亚洲av第一区精品v没综合| av一本久久久久| 少妇 在线观看| 久久亚洲精品不卡| 久久天堂一区二区三区四区| 亚洲男人天堂网一区| 国产欧美亚洲国产| 高清黄色对白视频在线免费看| 欧美成狂野欧美在线观看| 熟女少妇亚洲综合色aaa.| 亚洲成人免费电影在线观看| 亚洲国产精品一区二区三区在线| 搡老乐熟女国产| 精品福利观看| 精品福利观看| 午夜精品国产一区二区电影| 亚洲av成人不卡在线观看播放网| av在线播放免费不卡| 午夜福利,免费看| 久久精品国产亚洲av高清一级| 亚洲avbb在线观看| 国产精品久久久久成人av| 国产精品国产高清国产av | 香蕉丝袜av| 满18在线观看网站| 如日韩欧美国产精品一区二区三区| 性少妇av在线| 精品第一国产精品| 国产成+人综合+亚洲专区| 女人被躁到高潮嗷嗷叫费观| 久久天堂一区二区三区四区| 日韩欧美一区二区三区在线观看 | 亚洲中文av在线| 视频在线观看一区二区三区| 十分钟在线观看高清视频www| 黄色怎么调成土黄色| 久久国产精品影院| 亚洲国产看品久久| 国产又爽黄色视频| 欧美黑人欧美精品刺激| 久久国产精品大桥未久av| 夜夜爽天天搞| 亚洲一区高清亚洲精品| 国产激情欧美一区二区| 在线观看www视频免费| 涩涩av久久男人的天堂| 亚洲国产精品sss在线观看 | 欧美精品啪啪一区二区三区| 午夜精品国产一区二区电影| 亚洲成av片中文字幕在线观看| 久久性视频一级片| av片东京热男人的天堂| 欧美中文综合在线视频| 人妻丰满熟妇av一区二区三区 | 中文亚洲av片在线观看爽 | bbb黄色大片| 丰满迷人的少妇在线观看| 国产无遮挡羞羞视频在线观看| 啦啦啦免费观看视频1| 亚洲久久久国产精品| 欧美日韩乱码在线| 18禁裸乳无遮挡免费网站照片 | 一区二区三区精品91| 精品午夜福利视频在线观看一区| 91老司机精品| 午夜日韩欧美国产| 欧美国产精品va在线观看不卡| 国产欧美日韩一区二区三区在线| videosex国产| 国产精品偷伦视频观看了| av天堂久久9| 伊人久久大香线蕉亚洲五| 狠狠婷婷综合久久久久久88av| 色综合欧美亚洲国产小说| 国产成人av激情在线播放| 日日爽夜夜爽网站| 成年人午夜在线观看视频| 亚洲精品在线观看二区| 欧美激情 高清一区二区三区| 99国产精品一区二区蜜桃av | 女人被躁到高潮嗷嗷叫费观| 18禁国产床啪视频网站| 大陆偷拍与自拍| 久久性视频一级片| 亚洲综合色网址| www日本在线高清视频| ponron亚洲| 久久久久精品人妻al黑| 精品一区二区三区四区五区乱码| 国产成人精品久久二区二区免费| 黑人巨大精品欧美一区二区mp4| 91成人精品电影| 日韩欧美一区视频在线观看| 99精国产麻豆久久婷婷| 后天国语完整版免费观看| 亚洲精品久久午夜乱码| 免费在线观看影片大全网站| 在线观看一区二区三区激情| 国产精品九九99| 亚洲五月天丁香| 国产精品欧美亚洲77777| 国产亚洲精品一区二区www | 国产成人影院久久av| av一本久久久久| 国产精品乱码一区二三区的特点 | 中亚洲国语对白在线视频| 18禁裸乳无遮挡动漫免费视频| 午夜91福利影院| 欧美精品高潮呻吟av久久| 欧美精品人与动牲交sv欧美| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 久久精品亚洲精品国产色婷小说| 久久久久久人人人人人| 国产精品免费一区二区三区在线 | 久久精品国产99精品国产亚洲性色 | 精品人妻熟女毛片av久久网站| 亚洲国产欧美一区二区综合| 精品午夜福利视频在线观看一区| 天天躁夜夜躁狠狠躁躁| 在线观看日韩欧美| 老熟妇仑乱视频hdxx| 久久精品91无色码中文字幕| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 久久国产乱子伦精品免费另类| 亚洲av成人av| 国产精品久久久久久人妻精品电影| 国产成人免费观看mmmm| 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 亚洲精品一二三| 麻豆国产av国片精品| 一级毛片高清免费大全| 动漫黄色视频在线观看| 欧美最黄视频在线播放免费 | 黄色片一级片一级黄色片| 久久久久久久久免费视频了| 久久久精品免费免费高清| 91麻豆av在线| 亚洲国产看品久久| 美女扒开内裤让男人捅视频| 在线观看免费视频网站a站| 久久精品国产99精品国产亚洲性色 | 国产精品亚洲av一区麻豆| av线在线观看网站| 精品无人区乱码1区二区| 悠悠久久av| 青草久久国产| 老汉色∧v一级毛片| 青草久久国产| 婷婷成人精品国产| 在线观看日韩欧美| 自线自在国产av| 在线观看免费午夜福利视频| 自线自在国产av| 免费在线观看黄色视频的| 丝袜在线中文字幕| 午夜精品久久久久久毛片777| 国产一区二区三区视频了| 国产成人影院久久av| 一级黄色大片毛片| 男人的好看免费观看在线视频 | 宅男免费午夜| 精品久久久久久久久久免费视频 | 大陆偷拍与自拍| 亚洲av片天天在线观看| 俄罗斯特黄特色一大片| 99国产精品一区二区三区| 久久中文字幕一级| 免费在线观看视频国产中文字幕亚洲| 高清黄色对白视频在线免费看| 99国产精品一区二区三区| 天堂中文最新版在线下载| 国产亚洲精品一区二区www | 欧美日韩黄片免| 久久性视频一级片| 成人av一区二区三区在线看| 亚洲欧美精品综合一区二区三区| 亚洲av熟女| 欧美性长视频在线观看| 欧美在线一区亚洲| 亚洲五月色婷婷综合| 在线观看午夜福利视频| 妹子高潮喷水视频| 欧美日韩福利视频一区二区| 丝瓜视频免费看黄片| 国产麻豆69| 国产精品综合久久久久久久免费 | 国产精品九九99| 一级毛片精品| 久久久久精品人妻al黑| 日日爽夜夜爽网站| 成人亚洲精品一区在线观看| 美国免费a级毛片| 欧美激情高清一区二区三区| 色老头精品视频在线观看| 高清黄色对白视频在线免费看| 成年人午夜在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 午夜免费观看网址| 高潮久久久久久久久久久不卡| 十八禁网站免费在线| 精品国内亚洲2022精品成人 | 国产aⅴ精品一区二区三区波| 欧美老熟妇乱子伦牲交| 精品国内亚洲2022精品成人 | 9色porny在线观看| 国产午夜精品久久久久久| 亚洲色图综合在线观看| 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 亚洲国产精品一区二区三区在线| 一二三四在线观看免费中文在| 久久久久国内视频| 悠悠久久av| 黄色片一级片一级黄色片| 婷婷精品国产亚洲av在线 | 宅男免费午夜| av国产精品久久久久影院| 五月开心婷婷网| 午夜91福利影院| 免费女性裸体啪啪无遮挡网站| 老司机深夜福利视频在线观看| 最近最新中文字幕大全免费视频| 一级黄色大片毛片| 免费观看a级毛片全部| 91字幕亚洲| 夜夜夜夜夜久久久久| 午夜福利视频在线观看免费| 一级片'在线观看视频| 久久久精品免费免费高清| 制服人妻中文乱码| 岛国毛片在线播放| av电影中文网址| 亚洲精品中文字幕一二三四区| a级毛片黄视频| av超薄肉色丝袜交足视频| 一二三四在线观看免费中文在| 一区二区三区精品91| 免费在线观看影片大全网站| 日本撒尿小便嘘嘘汇集6| 亚洲国产毛片av蜜桃av| av一本久久久久| 国产精品二区激情视频| 午夜日韩欧美国产| 国产成人精品无人区| 麻豆成人av在线观看| 成人手机av| 18禁观看日本| 狂野欧美激情性xxxx| 精品人妻在线不人妻| 久久ye,这里只有精品| www.999成人在线观看| 久久国产亚洲av麻豆专区| 人人妻人人澡人人爽人人夜夜| 国产成+人综合+亚洲专区| 国产精品久久久久久人妻精品电影| 在线观看免费视频日本深夜| 日韩有码中文字幕| 亚洲国产毛片av蜜桃av| 欧美精品一区二区免费开放| 午夜两性在线视频| 一级作爱视频免费观看| 久久国产乱子伦精品免费另类| 亚洲精品国产一区二区精华液| 午夜福利一区二区在线看| 黄片播放在线免费| 国产精品自产拍在线观看55亚洲 | 青草久久国产| 国产在线一区二区三区精| 久久久久国产一级毛片高清牌| 国产免费av片在线观看野外av| 国产成人av激情在线播放| 人人妻人人澡人人爽人人夜夜| 午夜91福利影院| cao死你这个sao货| 热re99久久精品国产66热6| 视频区图区小说| 国产在线观看jvid| 最近最新中文字幕大全电影3 | 国产在线精品亚洲第一网站| 成人亚洲精品一区在线观看| 色老头精品视频在线观看| 一区二区三区精品91| 自线自在国产av| 精品午夜福利视频在线观看一区| 精品乱码久久久久久99久播| 看黄色毛片网站| 99国产精品一区二区蜜桃av | 女人精品久久久久毛片| 午夜福利一区二区在线看| 精品少妇久久久久久888优播| 成在线人永久免费视频| av免费在线观看网站| 欧美色视频一区免费| 国产野战对白在线观看| 国产片内射在线| 亚洲熟妇熟女久久| av不卡在线播放| 亚洲人成电影观看| 国产深夜福利视频在线观看| 午夜久久久在线观看| 午夜福利,免费看| 亚洲午夜理论影院| 中亚洲国语对白在线视频| 飞空精品影院首页| 久久久久视频综合| 一本一本久久a久久精品综合妖精| 色综合婷婷激情| av网站免费在线观看视频| 999精品在线视频| 午夜精品久久久久久毛片777| 999精品在线视频| 一级作爱视频免费观看| 久久精品国产综合久久久| 变态另类成人亚洲欧美熟女 | 亚洲欧美精品综合一区二区三区| 51午夜福利影视在线观看| 午夜日韩欧美国产| 19禁男女啪啪无遮挡网站| 日日夜夜操网爽| 啦啦啦 在线观看视频| 脱女人内裤的视频| 久久久久精品人妻al黑| 高清在线国产一区| 国产精品久久久人人做人人爽| 国产精品久久电影中文字幕 | 欧美日韩成人在线一区二区| 国产精品免费视频内射| 一级,二级,三级黄色视频| 亚洲三区欧美一区| 久久性视频一级片| 欧美日韩亚洲高清精品| 中出人妻视频一区二区| 免费在线观看完整版高清| 黑人欧美特级aaaaaa片| 久久久水蜜桃国产精品网| 一本综合久久免费| 国产精品成人在线| 天天躁狠狠躁夜夜躁狠狠躁| 淫妇啪啪啪对白视频| 中文字幕人妻丝袜制服| 久久久精品区二区三区| 国产单亲对白刺激| 亚洲av成人av| 亚洲一区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三卡| 丁香六月欧美| a级毛片在线看网站| 韩国精品一区二区三区| 男女高潮啪啪啪动态图| 天天添夜夜摸| 丰满饥渴人妻一区二区三| 亚洲一区高清亚洲精品| 亚洲成av片中文字幕在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 一级作爱视频免费观看| av天堂在线播放| 中文字幕最新亚洲高清| 国产有黄有色有爽视频| 999久久久精品免费观看国产| 成年女人毛片免费观看观看9 | 日韩欧美在线二视频 | www日本在线高清视频| 亚洲国产欧美一区二区综合| 亚洲熟女毛片儿| 人成视频在线观看免费观看| 三级毛片av免费| 三上悠亚av全集在线观看| 国产av又大| 久久久久久免费高清国产稀缺| 久久久国产成人免费| 日韩中文字幕欧美一区二区| 在线观看免费日韩欧美大片| 在线观看www视频免费| 高清毛片免费观看视频网站 | 两性夫妻黄色片| 精品少妇一区二区三区视频日本电影| 国产成人影院久久av| 九色亚洲精品在线播放| 亚洲人成电影观看| 国产99白浆流出| 欧美精品一区二区免费开放| 免费在线观看影片大全网站| 一进一出抽搐动态| 亚洲少妇的诱惑av| 中出人妻视频一区二区| 午夜福利乱码中文字幕| 国产高清国产精品国产三级| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 女性被躁到高潮视频| 这个男人来自地球电影免费观看| 身体一侧抽搐| 国产精品秋霞免费鲁丝片| 天天躁狠狠躁夜夜躁狠狠躁| 久久天躁狠狠躁夜夜2o2o| 黑人欧美特级aaaaaa片| 精品人妻熟女毛片av久久网站| 免费一级毛片在线播放高清视频 | 久久人妻av系列| 一a级毛片在线观看| 大陆偷拍与自拍| 亚洲成人手机| 91av网站免费观看| 热99国产精品久久久久久7| 老鸭窝网址在线观看| 久久香蕉精品热| 香蕉国产在线看| 亚洲专区中文字幕在线| 老司机午夜十八禁免费视频| 色在线成人网| 一级片'在线观看视频| 在线永久观看黄色视频| 亚洲欧洲精品一区二区精品久久久| 99精品欧美一区二区三区四区| 精品少妇久久久久久888优播| 亚洲av日韩精品久久久久久密| 91在线观看av| 狠狠婷婷综合久久久久久88av| 超碰成人久久| av网站免费在线观看视频| 亚洲五月色婷婷综合| 久久精品亚洲av国产电影网| 欧美 日韩 精品 国产| 首页视频小说图片口味搜索| av在线播放免费不卡| 成人18禁高潮啪啪吃奶动态图|