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

    Wind tunnel investigation of different engine layouts of a blended-wing-body transport

    2023-10-25 12:12:16ZhengCUIGuojunLAIQifengWANGYuLIANGYihuCAO
    CHINESE JOURNAL OF AERONAUTICS 2023年9期

    Zheng CUI, Guojun LAI, Qifeng WANG, Yu LIANG, Yihu CAO

    a School of Aeronautic Science and Engineering, Beihang University, Beijing 100083, China

    b COMAC Beijing Aircraft Technology Research Institute, Beijing 102211, China

    KEYWORDS PAI;Aerodynamics;BWB;CFD;Civil aviation;Engine layout;Transport aircraft;Wind tunnels

    Abstract A 2% scale, cruising version of a 450-seat class Blended-Wing-Body (BWB) transport was tested in the China Aerodynamic Research and Development Center’s FL-26 2.4-by-2.4-meter subsonic wind tunnel.The focus of the wind tunnel test was to investigate the aerodynamic performance of the latest BWB transport design, which would also aid in choosing a final engine arrangement in the three most potential engine integration layouts.The wind tunnel model can be tested with and without the nacelle and has three sets of different nacelle/tail integration positions.Computational Fluid Dynamics (CFD) simulations were performed in engine-aircraft integration design to find appropriate nacelle installing parameters of each layout.The comparison of CFD with experimental results shows good agreement.Wind tunnel measurements indicate that the tail-mounted engine layout produces the minimum drag penalty, while the fuselage-mounted engine layout increases drag the most.Experimental pressure measurement illustrates the effect of nacelle integration on the wing-body surface pressure distribution.This experimental and numerical research provides a reference for future BWB Propulsion-Airframe Integration (PAI) design.

    1.Introduction

    The United States and Europe have carried out research projects such as Environmentally Responsible Aviation (ERA)1–4 and Clean Sky5intending to maintain the leading edge of the aviation industry and seize the commanding heights of the future aviation market competition, and dozens of unconventional aircraft layout transport configurations such as Blended-Wing-Body (BWB) transport have been proposed and investigated.6–10The BWB layout refers to the shape of a full-lifting surface aircraft in which the wings and the fuselage are highly fused.11–15The BWB layout can reduce the wetted area of the whole aircraft under a given load requiring,and thus reduce the frictional drag.16–20Compared with the traditional conventional layout, the BWB cruise efficiency is improved by 15%-20%,21–24and it has the advantages of reducing noise, emissions, and structural weight, and other potential advantages.These advantages are represented by Boeing’s BWB-450 concept25and X-48 series demonstrator.26,27The academia and industry’s research on BWB is gradually deepening.The US X-48B and X-48C sub-scale BWB flying demonstrator verifies,by flight test,the feasibility,flight control, noise, fuel consumption, and other aspects of BWB layout,and the rich design experience and flight test data have been obtained.In NASA’s ERA project, the ERA-0009H1 is the result of refinement conducted under the ERA-0009A design that was developed under the ERA Advanced Vehicle Concepts (AVC) contract with Boeing.28Later, the Boeing-MIT team proposed BWB layout design H3 for NASA’s N+3 plan.29European ACFA2020 research project also selected BWB layout for next-generation civil aircraft.30Airbus built a demonstrator named MAVERIC(Model Aircraft for Validation and Experimentation of Robust Innovative Controls)for a possible full-scale BWB airliner and claimed that this design can reduce up to 20% of fuel.31The Russian Central Aerohydrodynamic Institute(TsAGI) launched the 900-seat BWB layout airliner program FW-900 back in 1996, and continued to design hybrid wingbody aircraft including Hybrid Wing-Body(TsAGI-IWB),lifting body layout(TsAGI-LF)and Flying Wing layout(TsAGIFW).32A series of research works and wind tunnel experiments have also been performed in Chinese academia, such as BWB-150, BWB-300 and BWB-330 proposed by Northwestern Polytechnical University (NPU).33,34

    The aft engine location of the BWB allows for plenty of installation options.Candidate installation concepts include podded with pylon, upper or lower surface inlet with S-duct,Boundary Layer Ingestion (BLI), or distributed propulsion.The aft engine location on the BWB offers the opportunity to ingest the boundary layer generated on the center body forward of the inlets.Boeing Company under the Ultra Efficient Engine Technology/Propulsion Airframe Integration Project has performed the study to assess the improvements in the BWB’s mission range that occur when podded engines are replaced with the embedded BLI system.35–37Numerical and experimental activities conducted by French Aerospace Lab(ONERA) within the RAPRO2 project allow an in-depth investigation and quantification of the mechanisms of BLI.Also, the experimental results are used to establish and validate a numerical simulation approach for future BLI propulsion system design.38Researchers39–41have been examining distributed propulsion systems that break the traditional twin turbofan engine configuration into arrays of many smaller engines.These arrangements create the opportunity for drag reduction and noise improvement.But the smaller turbofans are generally less efficient than large ones due to the practical limits of the core size and turbine inlet temperature.Distributed propulsion configurations will therefore require innovative electric turbo compound propulsion solutions to minimize these penalties.Simply mounting the engine on pylon is the most practical option, but increased wetted area and weight plus nose-down thrust moment are detractors of this installation.Lockheed Martin (LM) has teamed with both the Air Force Research Laboratory (AFRL) and the National Aeronautics and Space Administration (NASA) to develop a new strategic freighter called the Hybrid Wing Body (HWB) and validate its performance.42A Boeing 757 sized HWB configuration with over-wing engine nacelles was proposed and investigated.The results documented in this report indicate that an HWB-based commercial freighter with over-wing engine layout would use approximately about 7% less fuel than an advanced tube-and-wing configuration using the same technology levels.

    To sum up, many BWB civil transport configurations have been proposed and investigated,covering a variety of research on aerodynamics, structure, stability, and flight control.The previous research has considerably improved the engineering feasibility of the BWB transport concept.However, unlike conventional civil aircraft that generally integrates engines by wing-mounted arrangement, the academia and industry still do not have a general agreement on the question: which is the best way to integrate engines onto the BWB aircraft.And there is a devastating shortage of related wind tunnel data to support the engine integration design.So, in this research,we have investigated the three most popular engine arrangements of BWB civil transport configurations, which are wing-mounted layout, fuselage-mounted layout, and tailmounted layout.We hope that this work will give a good reference for future BWB aircraft design and push the BWB forward as a relatively clear pathway for future aircraft development.

    2.Experimental setup

    2.1.Test description

    Based on the needs of the future civil aircraft market, Commercial Aircraft Corporation of China(COMAC)Beijing Aircraft Technology Research Institute (BATRI) has launched pre-researches on a new generation of passenger aircraft and conducted comprehensive trade-offs of various layouts for future subsonic transport.The prototype of this research is a 450-seat BWB configuration with a cruise Mach number of 0.85 and a cruise altitude of 11000 m(Fig.1).The wind tunnel research is carried out in the FL-26 wind tunnel of the China Aerodynamics Development and Research Center.The FL-26 wind tunnel is the present largest ejector-type semirecirculation transient transonic wind tunnel in China.The size of its test section is 2.4 m×2.4 m,and the maximum Reynolds number per unit chord length test can reach 70 million/m with the capability to test Mach numbers from 0.4 to 1.4.Two primary goals exist for this research.The first is to study the subsonic and transonic aerodynamic characteristics of the BATRI BWB configuration,and the second is to investigate the effects of different engine integration layouts on the aircraft’s overall aerodynamic performance.

    The design of the wind tunnel test model followed the idea of building blocks and modular design.The wind tunnel test starts with the basic wing-body assembly configuration, and then gradually adds other components for comparison,including the very basic wing-body-tail assembly (baseline) configuration, wing-mounted engine configuration, fuselagemounted engine configuration, and tail-mounted engine configuration.

    Fig.1 BATRI 450-seat class BWB transport.

    Conventional static tests are conducted at a dynamic pressure of 11.4×103Pa,corresponding to a Reynolds number of 8.99 × 106based on mean aerodynamic chord.Data are acquired at angle of attack from —2° to 8°.For all tests, sixcomponent force and moment data were acquired using an internal strain-gauge balance.To minimize Reynolds number effects, transition tripping strip is used to fix transition at 7% chord on the upper and lower surface of the wing, Vtail, and nacelle.

    2.2.Model features

    According to the width of the test section of FL-26,the scale of the model is set to 1:50.The model has a wingspan of 1.6 m and the average aerodynamic chord length is 509 mm.The support form of the model is abdominal support.The installation of the test model is shown in Fig.2.

    The maximum thickness position of the BWB test model is about 135.7 mm, and its distance from the nose is about 312 mm.An internal strain-gauge balance is installed inside the central lift body and connected to the abdominal support.The test model has a total of 181 pressure measuring holes with an inner diameter of 0.4–0.6 mm, which is respectively arranged on eight spanwise positions distributed along the wingspan direction.Their relative positions to the symmetry plane of the aircraft are 0%, 17%, 23.5%, 30%, 36.5%,48%, 67.5% and 93.5%, respectively.

    The prototype BWB transport is designed as a twin-engine with V-tail(see Fig.2).The V-tails are installed on the back of the center lift body.Compared with an integrated V-tail at the wingtip, this tail layout can achieve a higher tail capacity and better aerodynamic and structural efficiency of the V-tail itself.That can further benefit the directional trim under the condition of single-engine failure and also increase the directional stability of the aircraft when flying at a high angle of attack.The three engine integration layouts being investigated herein are the fuselage-mounted, wing-mounted, and tail-mounted layouts.

    The fuselage-mounted layout puts nacelles on the back of the fuselage between the V-tails, which is the most common engine arrangement in the present BWB transport design (see Fig.3 and Fig.4).

    The second is the wing-mounted layout.This layout is to put the nacelle over the wing through a pylon that is connected to the wing’s trailing edge and extends backward to avoid the mutual aerodynamic interference between the nacelle and the wing as much as possible (see Fig.5 and Fig.6).

    Fig.2 Installation of wind tunnel model.

    Fig.3 Fuselage-mounted layout geometry.

    Fig.4 Fuselage-mounted layout model.

    Fig.5 Wing-mounted layout geometry.

    Fig.6 Wing-mounted layout model.

    Fig.7 Tail-mounted layout geometry.

    The last one is the tail-mounted layout.This layout integrates the engine nacelle at the root of the V-tail.Then the nacelle-tail assembly is installed as a whole part (see Fig.7 and Fig.8).

    For the above three engine layouts, the fuselage-mounted layout has the acoustic benefits of the shielding effect of the engines,42and neither of the other two layouts will be equally effective since only part of the engine was shielded; the wingmounted layout has an unloading effect on the wing, which can produce a structural benefit; the tail-mounted configuration has structural reinforcement on the root of V-tail; the fuselage-mounted one also has the advantage of increasing less yaw moment at the time of one engine failure.All the above highly complex factors have to be considered comprehensively in the design progress.However, this wind tunnel test only focuses on the aerodynamic characteristics and the trade-offs between aerodynamics,and other aspects are beyond the scope of this paper.

    2.3.CFD results overview

    CFD simulations were performed to aid in Propulsion-Airframe Integration (PAI) design, and an in-house software S-Flow package based on the Reynolds average Navier-Stokes (RANS) equation was used to simulate the Three-Dimensional (3D), compressible, steady flow over the BWB prototype.Roe’s scheme is used for space discretization.The Lower-Upper Symmetric Gauss-Seidel (LU-SGS) method is used for time advancement.And k-ω Shear Stress Transport(SST) equation is used for the turbulence model.

    For wing-body assembly design, we use a Free-Form Deformation (FFD) approach to parametrize the geometry.The FFD volume parametrizes the geometry changes rather than the geometry itself,resulting in a more efficient and compact set of geometry design variables, thus making it easier to handle complex geometric manipulations.In the aspect of the optimization algorithm, we use a gradient-based optimizer combined with adjoint gradient evaluations to solve the problem efficiently.Afterward,the CFD-assisted design progress of the nacelle installation position is carried out for each engine layout.The investigated variables include the nacelle height and chordwise position.For each engine layout,at least 30 sets of installation positions are calculated by CFD at the corresponding spanwise position (see Fig.9 and Fig.10).Finally,the full configuration of the integrated nacelle, pylon, and Vtail are evaluated.The overall design principle of nacelle integration is to avoid the boundary layer on the upper wing surface being sucked into the nacelle, and at the same time to ensure that the drag increase caused by the nacelle/pylon integration is minimized.The selected nacelle installation positions of different engine layouts are shown in Fig.11, and the CFD drag prediction results of each integrated configuration under the cruising condition are shown in Table 1.

    Fig.8 Tail-mounted layout model.

    Fig.9 Nacelle installation position.

    Fig.10 CFD drag prediction of nacelle installation.

    3.Results

    3.1.CFD and experimental comparison

    Fig.11 Top view and back view of nacelle installation position.

    Table 1 CFD drag prediction of different engine layouts.

    For conventional tube-and-wing configuration,the model support system can be attached to the tubular fuselage.Since the support is far away from the wing, the aerodynamic interference of the support on the wing can be well avoided.The BWB configuration, however, uses an integrated lift surface design where the fuselage and wing are highly fused.The aerodynamic interference of the model support on the wing-body assembly is particularly obvious.Therefore, the Reynolds number correction and support interference correction of the wind tunnel test data of the BWB configuration are quite difficult, and the relevant test data correction methods are still immature.Considering the above reasons, it is very difficult to directly compare the CFD data of the prototype with the wind tunnel test data.To verify the accuracy of the CFD numerical method, we have conducted a CFD analysis of the wing-body assembly model that has the model support integrated.The Reynolds number of the numerical simulation is consistent with the actual experimental conditions.Fig.12 shows the structured grid used in the above CFD simulation,the total cell number is 11.7 × 106, and the simulated Mach number and Reynolds number are 0.85 and 8.99 × 106,respectively.

    Fig.13 shows the comparison of CFD and experimental data for lift and drag coefficients, where CLrepresents the lift coefficient,CDrepresents the drag coefficient,AOA represents the angle of attack, respectively.It can be seen that the CFD prediction of the lift coefficient is in good agreement with the test results,and the CFD prediction of the drag coefficient is slightly higher than the test results, but the deviation is within an acceptable range.Fig.14 and Fig.15 show the comparison between the result of CFD simulated and experimental surface pressure measurement, where Cprepresents the pressure coefficient, X/C represents the local chord-wise position divided by the local chord length.The corresponding spanwise positions are 23.5% and 67.5%, respectively.The CFD and experimental comparisons above show that the S-Flow results are in excellent agreement with experimental results both for the force coefficients and surface pressure measurement.This agreement further proved that the grid generation strategy and numerical method can correctly predict overall aerodynamic characteristics and well support the design of BWB configuration.

    Fig.12 Structured grid of wing-body assembly model with support.

    Fig.13 Comparison of lift coefficient and drag coefficient between CFD and experiment.

    Fig.14 Pressure coefficient comparison at spanwise position 23.5%.

    Fig.15 Pressure coefficient comparison at spanwise position 67.5%.

    3.2.Force measurement results

    Fig.16 Lift coefficients comparison.

    Fig.17 Polar curves comparison.

    Fig.18 Pitching moment coefficients comparison.

    The wind tunnel test results of lift coefficient, polar curve,pitch moment coefficient,and lift-to-drag ratio K versus Angle of Attack(AOA)of the baseline configuration and three different nacelle-integrated configurations are shown in Figs.16-19.It should be noted that the geometric irregularity of the BWB fuselage makes it quite challenging to arrange a regular-size internal strain-gauge balance in the test model, so a ponysize internal strain-gauge balance is used in this experiment.However, during the experiment, we found that the stiffness of the pony-size balance was not enough to suppress test model oscillation at angles of attack over 6°.For safety reasons, in the later tests of the tail-mounted configuration and the fuselage-mounted configuration,the test angle of attack is limited to 5.5°.This adjustment resulted in inconsistencies in the angle of attack range of the data.However, we do not think it will affect the validity of the experimental data, because we did not focus on the high angle of attack characteristics in this experiment.

    Fig.19 Lift-to-drag ratio comparison.

    From the comparison of lift coefficients (see Fig.16), we can see that the integration of the engine nacelle will lead to a reduction in the lift of the basic configuration.Among them,the lift loss caused by the fuselage-mounted and tail-mounted engine layouts is about the same,and the lift coefficient reduction caused by the wing-mounted nacelle is the most obvious.It indicates that the wing-mounted engine layout deteriorates wing performance.So, the wing/engine integrative aerodynamic optimization progress should be performed if this engine arrangement is adapted.From the polar curve(Fig.17), we can see an overall drag increase after integrating the engine nacelle.Compared with the baseline configuration at the lift coefficient of 0.1, the drag increment for the fuselage-mounted, wing-mounted, and tail-mounted nacelle configuration is 55.3, 8.7, and 7.3 counts, respectively.So,the fuselage-mounted engine layout has brought a significantly higher drag cost than the other two engine-integrated configurations.The CFD simulation results also indicate that the fuselage-mounted engine arrangement leads to higher wave drag and interference drag because of the very close distances between nacelles and nacelle to the V-tail.From the comparison of the surface pressure coefficient (see Fig.20 and Fig.21), it can be observed that shock waves appear on the nacelle of the wing-mounted configuration.We believe that this phenomenon is caused by the mutual aerodynamic interference of the nacelle and tail and increases overall drag.

    From the comparison of the moment coefficients (see Fig.18), the pitching moment coefficient maintains a linear change in the range of 0°–5° angles of attack, and all the moment coefficient curves have an upturn at a 5° angle of attack.The three nacelle integration methods have different effects on the pitching moment, but the wind tunnel test data is still inadequate in evaluating which engine layout has a better moment character because the engine thrust influence has to be taken into account.The lift-to-drag ratio varies with the angle of attack as shown in Fig.19.Here we must point out that the test results of the lift-to-drag ratio of each configuration are significantly lower than the CFD results at the designed Reynolds number.This is because the Reynolds number of this test is significantly lower than the flight Reynolds number of the prototype.The flight angle of attack corresponding to the max lift-to-drag ratio is also about 2° higher because of the Reynolds number effects.(The test model has its best lift-to-drag performance at an AOA of 5°,but the prototype aircraft was designed for cruising at an AOA of 3°).At an angle of attack of 5°,the tail-mounted nacelle configuration has the best performance on lift-to-drag ratio,compared to the baseline configuration, it has a lift-to-drag deduction of 0.55,and the wing-mounted nacelle configuration and fuselagemounted configuration have a lift-to-drag deduction of 0.7 and 3.35, respectively.

    Fig.20 Pressure coefficient of tail-mounted configuration.

    Fig.21 Pressure coefficient of wing-mounted configuration.

    3.3.Pressure measurement results

    Fig.22 indicates the spanwise position(η)of the pressure measurement wing section.Figs.23-30 show the pressure measurement results of each wing section under Mach number 0.85 at an angle of attack of 5°.Generally speaking, the influence of the nacelle integration on the pressure distribution mainly concentrates on the inner wing upper surface and center lift body upper surface.From spanwise positions 0%-48%(see Figs.23-28), we can find that the nacelle integration will increase the pressure local near the trailing edge.That is because the nacelle has a blocking effect on the incoming flow, and the speed decrease of incoming flow leads to a pressure increase.The pressure distribution of the tail-mounted engine layout is the closest to the baseline pressure distribution,indicating that this nacelle integration approach has the mildest effect on baseline aerodynamics.On the contrary, the fuselage-mounted nacelle configuration has the most obvious effect on the upper surface of the center lift body (see Figs.23-25), but at spanwise positions of 30%-48% (see Figs.26-28), the effects of the wingmounted nacelle configuration exceed the influence of other configurations.It is also worth noting that at the spanwise position 67.5% wing section (see Fig.29) there is the opposite trend.The integration of the engine leads to a reduction in the pressure on the outer wing’s upper surface, which means that the wing-mounted and fuselage-mounted nacelle has increased the lift on the outer wing, and this is quite different compared to its effects on the inner parts.The reason for this inverted change is the airflow on the upper surface of the inner wing and the center lift body is partly blocked by the nacelle,resulting in more air flow being directed to the outer wing part.From this, we can know that the two nacelle integration layouts, the fuselage-mounted and wing-mounted layouts will lead to a lift decrease at the inner wing and lift increase at the outer wing, so the spanwise load distribution moves outward.The pressure distribution of the tail-mounted engine layout is the closest to the baseline configuration, resulting in the least impact on the spanwise load distribution.

    Fig.22 Spanwise positions of pressure measurement.

    Fig.23 Spanwise position: 0%.

    Fig.24 Spanwise position: 17%.

    Fig.25 Spanwise position: 23.5%.

    Fig.26 Spanwise position: 30%.

    Fig.27 Spanwise position: 36.5%.

    Fig.28 Spanwise position: 48%.

    Fig.29 Spanwise position: 67.5%.

    Fig.30 Spanwise position: 93.5%.

    Fig.31 shows the spanwise Cpdistributions near the 60%chord of the wing-body assembly at an AOA of 5°.The background image shows the location of the spanwise pressure taps as well as the engine nacelle to help interpret the pressure distributions.As expected,the suction pressures on the main wing vary with the change of engine layouts.The effect of the engine nacelle is noticeable, especially in the upstream spanwise Cpdistributions.The dip in the Cpdistributions varying with the nacelle/pylon integration position corresponds to the aerodynamic interference of nacelles and pylons on the wing-body assembly.The spanwise Cpdistributions of wing-mounted and fuselage-mounted layout appear to be dramatically suction pressure drop, implying a potential lift loss induced by the engine integration.

    Fig.31 Variation of spanwise Cp distribution along 60%chord.

    Fig.32 Spanwise lift distribution of CFD results.

    The spanwise lift distribution of the CFD results is illustrated in Fig.32,and the CFD calculation conditions are consistent with the wind tunnel test conditions.Comparing Fig.32 with Fig.31, we can see that the changing trend of the spanwise lift is consistent with the changing trend of the pressure coefficient.The nacelle integration does have a significant effect on the spanwise lift distribution.The wing root bending moment will also change according to the different integration methods of the nacelle.

    4.Conclusions

    A 2% scale version of a 450-seat class BWB transport wind tunnel test model was successfully tested in the FL-26 wind tunnel of the China Aerodynamics Development and Research Center.Wind tunnel measurements included force and pressure data as well as pressure measurement by 181 pressure holes.

    This study investigates the three mainstream engine integration layouts around the cruising Mach number 0.85.Force measurements indicate performance degradation after nacelles are integrated with baseline configuration.The lift coefficient reduction caused by the wing-mounted nacelle is the most obvious.An overall drag increase after integrating the engine nacelle is observed.Compared with the baseline configuration at the lift coefficient of 0.1,the drag increment for the fuselagemounted, wing-mounted, and tail-mounted nacelle configuration is 55.3, 8.7, and 7.3 counts, respectively.The tailmounted nacelle configuration has the best performance on lift-to-drag ratio, compared to the baseline configuration, it has a lift-to-drag deduction of 0.55, and the wing-mounted nacelle configuration and fuselage-mounted configuration have a lift-to-drag deduction of 0.7 and 3.35, respectively.

    Pressure measurement results indicate that the influence of the nacelle integration on the pressure distribution mainly concentrates on the inner wing upper surface and center lift body upper surface.The pressure distribution of the tail-mounted engine layout is the closest to the baseline pressure distribution, indicating that this nacelle integration approach has the mildest effect on baseline aerodynamics.However,other undesired effects brought about by tail-mounted layout require further investigation, including the control performance loss,structural strength, and vibration of the V-tail caused by the engine installation.

    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 study was co-supported by the Beijing Key Laboratory of Simulation Technology for Civil Aircraft Design, China (No.11WX08), and Innovation Foundation of COMAC Beijing Aircraft Technology Research Institute, China (No.Y16QT01).

    人人澡人人妻人| 亚洲国产高清在线一区二区三 | 午夜福利,免费看| 欧美黑人欧美精品刺激| 国产亚洲av嫩草精品影院| 国产精品免费视频内射| e午夜精品久久久久久久| 欧美大码av| bbb黄色大片| 少妇粗大呻吟视频| 国产麻豆69| 少妇裸体淫交视频免费看高清 | 99久久99久久久精品蜜桃| 好男人电影高清在线观看| 国产精品电影一区二区三区| 午夜亚洲福利在线播放| 亚洲成人免费电影在线观看| 一级a爱片免费观看的视频| 久久人妻熟女aⅴ| 一边摸一边抽搐一进一出视频| 涩涩av久久男人的天堂| 欧美丝袜亚洲另类 | 视频区欧美日本亚洲| 91精品国产国语对白视频| 免费在线观看亚洲国产| 一区二区三区激情视频| 美女高潮喷水抽搐中文字幕| 侵犯人妻中文字幕一二三四区| 97超级碰碰碰精品色视频在线观看| av超薄肉色丝袜交足视频| 最好的美女福利视频网| 看片在线看免费视频| 久久久水蜜桃国产精品网| 亚洲免费av在线视频| 伊人久久大香线蕉亚洲五| 欧美成人午夜精品| 99精品在免费线老司机午夜| 亚洲国产欧美一区二区综合| 在线av久久热| 制服丝袜大香蕉在线| 男女下面进入的视频免费午夜 | 黑人巨大精品欧美一区二区mp4| 久久午夜综合久久蜜桃| 丝袜美腿诱惑在线| 午夜两性在线视频| av天堂久久9| 黑丝袜美女国产一区| 一区二区三区高清视频在线| 搡老妇女老女人老熟妇| 美女免费视频网站| 操出白浆在线播放| 可以在线观看的亚洲视频| 在线观看舔阴道视频| 欧美老熟妇乱子伦牲交| 久久久国产成人免费| 亚洲欧美激情在线| 美女国产高潮福利片在线看| 亚洲片人在线观看| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频| 亚洲专区字幕在线| 女性生殖器流出的白浆| 一级黄色大片毛片| 久久精品国产综合久久久| 午夜福利成人在线免费观看| 久久九九热精品免费| 日韩大码丰满熟妇| 国产主播在线观看一区二区| 少妇熟女aⅴ在线视频| 日本五十路高清| 制服诱惑二区| 午夜福利18| 一级片免费观看大全| 色综合婷婷激情| 老司机福利观看| 99久久精品国产亚洲精品| 给我免费播放毛片高清在线观看| 美女大奶头视频| 两个人免费观看高清视频| 国产精品98久久久久久宅男小说| av在线天堂中文字幕| 少妇被粗大的猛进出69影院| 久久精品aⅴ一区二区三区四区| 亚洲狠狠婷婷综合久久图片| 亚洲九九香蕉| 精品无人区乱码1区二区| 丁香六月欧美| 午夜视频精品福利| 国产av一区二区精品久久| 叶爱在线成人免费视频播放| 欧美日韩一级在线毛片| 成人18禁高潮啪啪吃奶动态图| 91精品国产国语对白视频| 欧美乱色亚洲激情| 妹子高潮喷水视频| 在线视频色国产色| 长腿黑丝高跟| 国产又色又爽无遮挡免费看| 一二三四在线观看免费中文在| 成人精品一区二区免费| 精品久久久久久成人av| 亚洲国产精品sss在线观看| 香蕉久久夜色| 国产日韩一区二区三区精品不卡| 午夜福利欧美成人| 99国产精品一区二区三区| 欧美国产精品va在线观看不卡| 99香蕉大伊视频| 午夜亚洲福利在线播放| 精品国产美女av久久久久小说| 在线观看www视频免费| 欧美午夜高清在线| 亚洲人成网站在线播放欧美日韩| 免费在线观看日本一区| 久久这里只有精品19| 欧美日韩乱码在线| 国产成人精品久久二区二区91| 欧美黑人欧美精品刺激| 97人妻天天添夜夜摸| 国产一区在线观看成人免费| 亚洲成国产人片在线观看| 国产精品九九99| 欧美精品亚洲一区二区| 亚洲专区中文字幕在线| 亚洲精品在线观看二区| 精品久久久久久久人妻蜜臀av | 亚洲午夜精品一区,二区,三区| av天堂久久9| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| 欧美日韩乱码在线| 91麻豆av在线| 成人亚洲精品av一区二区| 老司机福利观看| 欧美性长视频在线观看| 夜夜看夜夜爽夜夜摸| 国产精品香港三级国产av潘金莲| 午夜亚洲福利在线播放| 国产成人av教育| 免费观看精品视频网站| 亚洲人成网站在线播放欧美日韩| 美女免费视频网站| 国产精品 国内视频| 首页视频小说图片口味搜索| 国产日韩一区二区三区精品不卡| 亚洲熟女毛片儿| 一区福利在线观看| 日本免费一区二区三区高清不卡 | 丰满的人妻完整版| 老司机午夜十八禁免费视频| 久久久国产成人免费| 午夜精品久久久久久毛片777| 国产精品秋霞免费鲁丝片| 精品日产1卡2卡| 欧美黑人精品巨大| 亚洲av美国av| 大香蕉久久成人网| 久久久精品国产亚洲av高清涩受| 国产精品久久电影中文字幕| 9191精品国产免费久久| 电影成人av| 在线永久观看黄色视频| 亚洲国产看品久久| 国产一区在线观看成人免费| 两个人免费观看高清视频| 十八禁网站免费在线| 一级,二级,三级黄色视频| bbb黄色大片| 免费久久久久久久精品成人欧美视频| 亚洲第一av免费看| 美女大奶头视频| 19禁男女啪啪无遮挡网站| 18禁美女被吸乳视频| 51午夜福利影视在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲最大成人中文| av欧美777| 淫秽高清视频在线观看| 热99re8久久精品国产| 可以免费在线观看a视频的电影网站| 久久久久久久久免费视频了| 黄片播放在线免费| 精品不卡国产一区二区三区| 亚洲视频免费观看视频| 欧美日韩乱码在线| 午夜久久久在线观看| 国产免费男女视频| 两个人看的免费小视频| 91老司机精品| 欧美丝袜亚洲另类 | 欧美老熟妇乱子伦牲交| 69av精品久久久久久| 真人一进一出gif抽搐免费| 久久人妻福利社区极品人妻图片| 女人被躁到高潮嗷嗷叫费观| 黄片小视频在线播放| 国产亚洲精品av在线| 超碰成人久久| 999精品在线视频| 19禁男女啪啪无遮挡网站| 69精品国产乱码久久久| 真人做人爱边吃奶动态| 涩涩av久久男人的天堂| 极品教师在线免费播放| 国产亚洲精品久久久久5区| 精品一区二区三区四区五区乱码| 欧美中文综合在线视频| 久久精品亚洲熟妇少妇任你| 亚洲专区中文字幕在线| 午夜福利免费观看在线| 男女之事视频高清在线观看| 一区在线观看完整版| 久久久久国产一级毛片高清牌| 亚洲情色 制服丝袜| 在线国产一区二区在线| 女同久久另类99精品国产91| 午夜精品国产一区二区电影| 巨乳人妻的诱惑在线观看| 19禁男女啪啪无遮挡网站| 亚洲成a人片在线一区二区| 超碰成人久久| 午夜精品国产一区二区电影| 88av欧美| 女性生殖器流出的白浆| 桃色一区二区三区在线观看| 777久久人妻少妇嫩草av网站| 久久精品国产亚洲av香蕉五月| 亚洲精品久久国产高清桃花| 国产片内射在线| 国产精品免费视频内射| 在线视频色国产色| 日韩欧美一区二区三区在线观看| 女生性感内裤真人,穿戴方法视频| 99riav亚洲国产免费| 9色porny在线观看| 悠悠久久av| 国产成人精品在线电影| 欧美黄色淫秽网站| 少妇 在线观看| АⅤ资源中文在线天堂| 国产精品久久久久久人妻精品电影| www.精华液| 免费无遮挡裸体视频| 亚洲激情在线av| 丝袜人妻中文字幕| 久久久久久久久免费视频了| 亚洲av成人不卡在线观看播放网| 国产精品av久久久久免费| 成人亚洲精品av一区二区| 一二三四在线观看免费中文在| 免费av毛片视频| 一级,二级,三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 在线观看www视频免费| 国产一区二区在线av高清观看| 国产精品秋霞免费鲁丝片| 不卡av一区二区三区| 精品一区二区三区视频在线观看免费| 久久久国产精品麻豆| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看 | 亚洲在线自拍视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本视频| 日韩精品青青久久久久久| 色综合站精品国产| 国内精品久久久久久久电影| 天天一区二区日本电影三级 | 老司机午夜福利在线观看视频| 久久草成人影院| 国产乱人伦免费视频| 精品人妻在线不人妻| 两性夫妻黄色片| 啦啦啦韩国在线观看视频| 欧美一区二区精品小视频在线| 亚洲九九香蕉| 最近最新中文字幕大全电影3 | 美国免费a级毛片| 在线观看一区二区三区| 国产精华一区二区三区| 国产成人av激情在线播放| 一级a爱视频在线免费观看| 久久精品国产清高在天天线| 欧美色视频一区免费| 不卡av一区二区三区| 亚洲激情在线av| 99riav亚洲国产免费| 黄色毛片三级朝国网站| 亚洲精品在线观看二区| 天天添夜夜摸| a级毛片在线看网站| 国产精品 欧美亚洲| 黄网站色视频无遮挡免费观看| 久久久久久大精品| 大码成人一级视频| 丝袜在线中文字幕| 男女下面插进去视频免费观看| 亚洲欧洲精品一区二区精品久久久| 侵犯人妻中文字幕一二三四区| www.999成人在线观看| 国产97色在线日韩免费| 制服人妻中文乱码| av天堂在线播放| 国产精品九九99| 久久久久国产一级毛片高清牌| 欧美乱码精品一区二区三区| av欧美777| 欧美日韩福利视频一区二区| 欧美乱码精品一区二区三区| 一级毛片精品| 日韩国内少妇激情av| 国产欧美日韩精品亚洲av| 精品午夜福利视频在线观看一区| 亚洲午夜精品一区,二区,三区| 色综合亚洲欧美另类图片| 美女 人体艺术 gogo| 久久久久国产精品人妻aⅴ院| 国产91精品成人一区二区三区| 99久久国产精品久久久| 侵犯人妻中文字幕一二三四区| av欧美777| 青草久久国产| 男男h啪啪无遮挡| 无遮挡黄片免费观看| 日韩高清综合在线| 亚洲精品av麻豆狂野| 午夜久久久在线观看| 久久草成人影院| 国产视频一区二区在线看| 久久精品影院6| 国产单亲对白刺激| 亚洲精品中文字幕在线视频| 人成视频在线观看免费观看| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三区在线| 满18在线观看网站| 免费在线观看日本一区| 久久亚洲精品不卡| 十分钟在线观看高清视频www| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 中文字幕av电影在线播放| 国产私拍福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| 久久国产精品影院| 黑丝袜美女国产一区| 大码成人一级视频| 天天躁夜夜躁狠狠躁躁| 母亲3免费完整高清在线观看| 国产高清videossex| 亚洲欧美日韩另类电影网站| 一夜夜www| 精品第一国产精品| 午夜免费成人在线视频| 久久精品成人免费网站| 一级毛片女人18水好多| av欧美777| 女人被狂操c到高潮| 99久久久亚洲精品蜜臀av| 夜夜爽天天搞| 精品卡一卡二卡四卡免费| 一二三四社区在线视频社区8| 午夜成年电影在线免费观看| 午夜免费成人在线视频| 国产熟女午夜一区二区三区| 搞女人的毛片| 久久精品影院6| 国产精品 欧美亚洲| 亚洲精品一区av在线观看| 高清毛片免费观看视频网站| 国产亚洲精品av在线| 日韩成人在线观看一区二区三区| 精品国产一区二区三区四区第35| 亚洲一卡2卡3卡4卡5卡精品中文| 搞女人的毛片| 亚洲精品美女久久久久99蜜臀| 黄色 视频免费看| 少妇裸体淫交视频免费看高清 | 欧美不卡视频在线免费观看 | 曰老女人黄片| 国产91精品成人一区二区三区| 久久精品成人免费网站| 久久久久久久久中文| 国产一区在线观看成人免费| 免费在线观看黄色视频的| 激情在线观看视频在线高清| 在线天堂中文资源库| 99久久精品国产亚洲精品| 精品国产超薄肉色丝袜足j| 亚洲精品中文字幕一二三四区| 九色国产91popny在线| 天天躁狠狠躁夜夜躁狠狠躁| 91国产中文字幕| 一本大道久久a久久精品| www.999成人在线观看| 午夜成年电影在线免费观看| 巨乳人妻的诱惑在线观看| 最新在线观看一区二区三区| 国产一区二区三区视频了| 最近最新免费中文字幕在线| 欧美国产精品va在线观看不卡| 国产三级黄色录像| 黄色a级毛片大全视频| 99国产精品免费福利视频| 欧美一级a爱片免费观看看 | 国产一卡二卡三卡精品| 亚洲情色 制服丝袜| 国产一区二区三区综合在线观看| 免费一级毛片在线播放高清视频 | 欧美 亚洲 国产 日韩一| 性色av乱码一区二区三区2| 亚洲一区中文字幕在线| a在线观看视频网站| 欧美日韩黄片免| 如日韩欧美国产精品一区二区三区| 少妇粗大呻吟视频| 97超级碰碰碰精品色视频在线观看| 乱人伦中国视频| x7x7x7水蜜桃| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 国产精品一区二区三区四区久久 | 国产亚洲精品av在线| 国产伦人伦偷精品视频| 欧美日韩精品网址| 一本久久中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 免费久久久久久久精品成人欧美视频| 两个人免费观看高清视频| 久久人妻av系列| 天堂动漫精品| 丝袜美足系列| 91成年电影在线观看| 精品久久久精品久久久| 亚洲国产欧美一区二区综合| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 韩国精品一区二区三区| 99久久精品国产亚洲精品| 日本欧美视频一区| 亚洲国产精品久久男人天堂| 亚洲熟妇中文字幕五十中出| 亚洲中文日韩欧美视频| 满18在线观看网站| 日本黄色视频三级网站网址| 精品国产一区二区久久| 97人妻精品一区二区三区麻豆 | 国产一区二区三区在线臀色熟女| 国产高清videossex| 长腿黑丝高跟| 亚洲专区字幕在线| 一个人观看的视频www高清免费观看 | 久久亚洲精品不卡| 国内精品久久久久精免费| 淫秽高清视频在线观看| 亚洲午夜精品一区,二区,三区| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av在线| 国产单亲对白刺激| 91麻豆精品激情在线观看国产| 久久影院123| 精品国产一区二区三区四区第35| 自拍欧美九色日韩亚洲蝌蚪91| 他把我摸到了高潮在线观看| 日本在线视频免费播放| 男女做爰动态图高潮gif福利片 | 成人三级黄色视频| 高清在线国产一区| 久久久水蜜桃国产精品网| 香蕉久久夜色| 午夜亚洲福利在线播放| 国产亚洲精品久久久久久毛片| 久99久视频精品免费| 免费看a级黄色片| www.精华液| www.999成人在线观看| 精品国产国语对白av| 亚洲av熟女| 亚洲人成网站在线播放欧美日韩| 国产成人av教育| 日韩 欧美 亚洲 中文字幕| xxx96com| 99国产综合亚洲精品| 美女国产高潮福利片在线看| 欧美色视频一区免费| 婷婷精品国产亚洲av在线| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 韩国av一区二区三区四区| 国产精品 国内视频| 午夜视频精品福利| 婷婷丁香在线五月| 午夜免费观看网址| 黄色毛片三级朝国网站| 69av精品久久久久久| 午夜两性在线视频| 69精品国产乱码久久久| 日本欧美视频一区| 高潮久久久久久久久久久不卡| 99国产综合亚洲精品| 久久国产精品影院| 国产麻豆69| 男女做爰动态图高潮gif福利片 | 美女高潮喷水抽搐中文字幕| 亚洲av成人一区二区三| 一个人观看的视频www高清免费观看 | 变态另类丝袜制服| 搞女人的毛片| 欧美成狂野欧美在线观看| 国产精品影院久久| 国产片内射在线| 亚洲av日韩精品久久久久久密| 激情在线观看视频在线高清| 亚洲一卡2卡3卡4卡5卡精品中文| 久久国产精品男人的天堂亚洲| 高清在线国产一区| 成人亚洲精品一区在线观看| 久热爱精品视频在线9| tocl精华| 亚洲第一av免费看| 国产99久久九九免费精品| 欧美日韩亚洲综合一区二区三区_| 精品不卡国产一区二区三区| 欧美色欧美亚洲另类二区 | 精品久久久久久久人妻蜜臀av | 久久精品亚洲精品国产色婷小说| 一级,二级,三级黄色视频| 欧美+亚洲+日韩+国产| 色av中文字幕| 免费高清视频大片| 欧美 亚洲 国产 日韩一| 欧美乱妇无乱码| 亚洲国产毛片av蜜桃av| 国产精品一区二区在线不卡| 久久久精品欧美日韩精品| 精品不卡国产一区二区三区| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| 国产亚洲精品av在线| 亚洲人成电影观看| 午夜福利18| 人人澡人人妻人| 香蕉丝袜av| 精品国产超薄肉色丝袜足j| 亚洲成国产人片在线观看| 黄片小视频在线播放| 亚洲专区字幕在线| 真人一进一出gif抽搐免费| 无遮挡黄片免费观看| 最新美女视频免费是黄的| 免费看美女性在线毛片视频| 亚洲中文字幕日韩| 亚洲免费av在线视频| 国产私拍福利视频在线观看| 久久精品亚洲精品国产色婷小说| 亚洲成av人片免费观看| 美女午夜性视频免费| 日韩大尺度精品在线看网址 | 久久久久久免费高清国产稀缺| 97超级碰碰碰精品色视频在线观看| 亚洲第一青青草原| 亚洲va日本ⅴa欧美va伊人久久| 好男人在线观看高清免费视频 | 最近最新免费中文字幕在线| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 天天添夜夜摸| 日韩中文字幕欧美一区二区| 欧美丝袜亚洲另类 | 成人永久免费在线观看视频| 十八禁人妻一区二区| 999精品在线视频| 久久久国产成人免费| 男女午夜视频在线观看| 午夜免费鲁丝| 一个人观看的视频www高清免费观看 | 91成年电影在线观看| 国产亚洲精品久久久久久毛片| 亚洲成国产人片在线观看| 99在线人妻在线中文字幕| 午夜激情av网站| 99久久国产精品久久久| 国产精品,欧美在线| 亚洲一区二区三区不卡视频| 一区在线观看完整版| 男女做爰动态图高潮gif福利片 | 国产一区二区激情短视频| 亚洲欧美精品综合一区二区三区| 热99re8久久精品国产| 欧美激情 高清一区二区三区| 国产成+人综合+亚洲专区| 亚洲成av人片免费观看| 男女下面进入的视频免费午夜 | 亚洲专区国产一区二区| 欧美丝袜亚洲另类 | 又黄又爽又免费观看的视频| 人人妻人人澡人人看| 国产av一区二区精品久久| 搡老熟女国产l中国老女人| 午夜免费激情av| 精品国产一区二区三区四区第35| 看黄色毛片网站| 国产色视频综合| 国产高清激情床上av| av电影中文网址| 亚洲色图 男人天堂 中文字幕| 国产高清激情床上av| 午夜福利影视在线免费观看| 香蕉国产在线看| 婷婷六月久久综合丁香| 久久精品国产综合久久久| 久久性视频一级片| 久久国产精品男人的天堂亚洲| 久久久久久久久中文| 制服诱惑二区|