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

    High-speed wind tunnel test of the CAE aerodynamic validation model

    2018-04-19 08:28:58RoyGEBBINKGanglinWANGMinZHONG
    CHINESE JOURNAL OF AERONAUTICS 2018年3期

    Roy GEBBINK,Ganglin WANG,Min ZHONG

    aGerman-Dutch Wind Tunnels,Marknesse 8316 PR,The Netherlands

    bChinese Aeronautical Establishment,Beijing 100012,China

    1.Introduction

    The global demand and visions1for more sustainable means of air transportation bring forth many challenges for all aeronautical disciplines,from the level of aircraft design up to aircraft operations.In the field of aerodynamics,design innovation for commercial civil aircraft commonly2–4focuses on the increase of the aerodynamic efficiency Ma×L/D,where Ma is the Mach number,L the lift,and D the drag.Optimizing Ma×L/D in the transonic flight regime is often challenging,as the drag associated with shockwave formations generally increases progressively when increasing the flight Mach number.In the framework of a design study for an Ma=0.87 long-haul business jet aircraft concept,5the Chinese Aeronautical Establishment(CAE)initiated a project to conduct a dual-purpose wind tunnel experiment.The first purpose of the experiment was to establish the aerodynamic cruise performance for multiple configurations of the aircraft.The second purpose was to establish a reliable dataset for validation of an in-house developed CFD code.To serve both purposes within a single wind tunnel campaign,the CAE Aerodynamic Validation Model(CAE-AVM)was developed,and an experiment was designed that consisted of parallel deployments of various measurement techniques.This paper addresses the design features of the CAE-AVM,along with a description of the wind tunnel experiment and the typical test results at a cruise Mach number.

    2.AVM design and manufacture

    2.1.Configuration

    The Ma=0.87 long-haul business jet aircraft concept is a configuration that consists of a narrow fuselage with a low mounted swept wing of a high aspect ratio,two aft-fuselage mounted engines,and a T-tail(see Fig.1).The overall length of the aircraft is 33 m and the wing span is 33.5 m,including the winglets.The aircraft targets a range of 11000–13000 km with a cruise Mach number of 0.87 at an altitude of 13 km.As the aircraft comprises a slender design,it is anticipated that a conventional metallic scale model would be severely affected by aeroelastic effects when tested in a pressurized wind tunnel.Refs.6–8have shown that for the transonic speeds of interest,aeroelastic model deformation(in particular wing twist)is an important uncertainty in a direct comparison between experimental and computational results.

    In order to establish a comprehensive experimental dataset for CFD validation purposes,the goal was defined to keep the wing deformation of the model within a predefined limit of no more than-1°twist at the wingtip for the envisioned cruise lift condition.To realize this target,the cruise Mach number for the wind tunnel test model was reduced from 0.87 to 0.85,which allowed the new wing to be thickened and the sweep angle to be reduced.In parallel,the winglets were omitted,and modifications were made to the baseline wing sections,in an attempt to remove trailing edge separations at the lower model scale Reynolds number attainable in the envisioned wind tunnel.Similarly,the shape of the engine mounting pylons was adapted to remove shock-induced separations at the model scale Reynolds number.The engines themselves were represented by through- flow-nacelles.The resulting geometry was denominated as the CAE-AVM.

    In collaboration with the German-Dutch Wind Tunnels(DNW)and the Netherlands Aerospace Centre(NLR),a project was defined for building an accurate wind tunnel model of the AVM,suitable for testing in the DNW High-Speed Tunnel(HST).Given the dimensions of the HST test section,the size of the AVM was selected so as to keep the ratio of the model frontal area to the test section’s cross-sectional area less than 1%and to keep the ratio of the wing span to the test section’s width less than 70%.As a result from these two wall interference considerations,9the CAE-AVM has a reference wing area of 0.208 m2,a reference span of 1.37 m,and a mean aerodynamic chord of 0.199 m.The wing has an aspect ratio of 9.0 and a quarter chord sweep angle of 35°.The aerodynamic design point,also referred to as the cruise lift condition,is Ma=0.85,CL=0.50,and Re=4.7×106(where CLis the lift coefficient based on the reference wing area,and Re the Reynolds number based on the mean aerodynamic chord).

    Fig.1 Artistic impression of the long-haul business jet,with its cruise Mach number of 0.87.

    2.2.Model features

    The wind tunnel model is designed for simultaneous measurements of the overall forces and moments,wing pressure distributions,and model attitude.The model design allows for testing of two main aircraft configurations:the fuselage/wing(BW)combination and the complete configuration(BWNVH)that consists of the fuselage/wing/pylon-nacelles/vertical-tail/horizontal-tail(see Fig.2).

    The core of the model is designed as a single central body,which houses an internal strain gauge balance developed by the NLR.The balance features an end-face flange that provides a model-balance interface with minimal hysteresis through the principle of a pre-stressed connection.The design of the model is such that the balance can be mounted on either a ventral z-sting support or a dorsal blade sting support.

    For the assessment of the wing load distribution,a total of 180 pressure taps were incorporated.The taps were distributed in six span-wise stations amongst the left and right wings,at η={0.20,0.35,0.45,0.55,0.65,0.75}.Because the inclusion of the pressure taps weakened the construction,special design considerations were implemented to keep stiffness losses to a minimum.For example,the pockets in the wings(required for installation of the pressure tubes)were provided with load-bearing covers;in the wing root area,special attention was given to the tube and cable routing whilst maintaining a firm fixation interface between the wings and the central balance house.For reasons of strength and stiffness of the wings,the central balance housing and the empennage were all built from Ramax stainless steel.All other model parts were built from an aluminum alloy.

    To verify that the design meets the predefined twist requirement,a finite element analysis(FEA)was conducted.As input for the FEA,a wing load distribution was used that was obtained from a CFD simulation for the nominal,nondeformed,aircraft geometry.The use of a load distribution for a non-deformed wing is judged to yield a conservative estimate of the wing twist,as the loading of the outboard wing sections is expected to be higher compared to a load distribution associated with a deformed wing shape.The FEA results confirmed that the wing twist should not exceed-1°at the tip,at the cruise lift condition.

    Fig.2 As-built CAE-AVM full-model assembly,with its cruise Mach number of 0.85.

    Upon completion of design,the actual manufacturing took place.During the manufacturing process,the accuracy of the model contour was closely monitored at the level of individual parts,as well as in-assembly.Hand finishing was applied to ensure smooth surfaces and no gaps and/or any forward facing steps between joints.After manufacturing,a Coordinate Measuring Machine(CMM)was used as independent verification that the as-built model was within tolerances.In terms of the as-built wing twist,the CMM confirmed that it was within±0.06°along the span.

    3.Test approach

    3.1.Wind tunnel facility

    The test facility is the DNW-HST10in Amsterdam,the Netherlands.The HST is a continuous- flow pressurized wind tunnel that allows for aircraft configuration testing at subsonic and transonic speeds:the Mach number range is from 0.2 up to 1.3.The shell structure of the wind tunnel(Fig.3)allows the facility to be pressurized or evacuated,enabling testing at a constant Reynolds number at variable Mach numbers,or variable Reynolds numbers at a constant Mach number.The attainable total pressure range is from 25 to 390 kPa,and the operational total temperature ranges between 275 and 320 K.The air in the circuit is driven by an axial compressor which comprises of several rows of rotor vanes with adjustable blade pitch;a control system allows the Mach number to be kept constant within±0.002 during operation.The compressor itself is powered by a 20 MW electrical drive system.

    The test section is of rectangular shape with a fixed width of 2.0 m and an adjustable height of either 1.6 m or 1.8 m.The side walls are solid,whilst the top and bottom walls are slotted with an openness ratio of 12%per wall.For the AVM experiment,the 1.8 m by 2.0 m slotted wall test section configuration was used.The test section is surrounded by a plenum chamber to accommodate ventilation of the air through the slotted walls.Downstream in the test section,a permanently installed vertical strut provides interfacing for several types of model support systems.

    3.2.Test set-ups

    The primary test set-up for the CAE-AVM consisted of a ventral z-sting with the so-called straight boom which was connected to the boom-base on the vertical strut downstream in the test section(Fig.4(a)).The boom-base is traversable in height and provides a pivot point for pitch angle variation.

    Fig.3 Artistic impression of the DNW-HST circuit.

    The kinematics was such that the Moments Reference Point(MRP)of the model could be kept on the tunnel centerline during pitch angle sweeps up to 10°.In this ventral set-up,the AVM was mounted upside-up on the z-sting via an internal balance.The z-sting was specifically selected for the primary set-up,because of its modest interference on the model aerodynamics,whilst providing the flexibility of testing the tail-off and tail-on model configurations.

    Aerodynamic interference of the primary z-support(i.e.,the ventral z-sting with the straight boom support),a secondary test set-up was used.For the secondary set-up,the AVM was mounted via the internal balance attached to the dorsal blade sting shown in Fig.4(b).In this set-up,the model was mounted upside-down to allow the dorsal sting to connect to a pitching mechanism underneath the test section floor.Between the dorsal sting and the pitching mechanism,an xtraverse was installed.The x-traverse and pitching mechanism could be operated so as to match the exact though inverted model position and orientation,as occurred during pitch angle sweeps with the primary ventral set-up.

    The aerodynamic influence of the primary set-up was established with two measurements in the dorsal set-up:one measurement with the model solely on the dorsal sting(Fig.4(b));the other measurement with the model mounted to the dorsal sting,whilst a dummy of the ventral z-support was present(Fig.4(c)).The latter consisted of the straight boom plus the so-called dummy z-sting,which entered the ventral model cavity without mechanical contact.The model loads were measured with the internal balance connected to the dorsal sting.In order to prevent any model-internal flow from the dorsal sting cavity to the z-sting cavity,the model was provided with a seal downstream from the balance.For the configuration without the dummy support,the ventral z-sting cavity was closed with an aft-fuselage plug so as to represent the nominal fuselage geometry.Hereby,the difference in balance readings between measurements with and without the dummy zsupport include the external interference effects due to the dummy z-support,plus any potential cavity(pressure)effects.This whole procedure relies on the assumption that the interference from the dorsal blade sting on the test object has no interaction with the(dummy)ventral z-support.This is deemed a reasonable assumption,thanks to the streamlined shape of the dorsal sting and its placement on the opposite side of the model so that there is no direct wake-interaction between the two support systems.

    Fig.4 Side views of the test set-ups.

    3.3.Instrumentation

    For measurements of the overall forces and moments,a sixcomponent strain gauge balance was used.The balance’s longitudinal load ranges are:±14600 N for normal force,±1,360 N for axial force,and ±930 N·m for pitching moment.The balance was calibrated in advance of the test,by means of a multi-component loading scheme,and has a calibrated uncertainty(i.e.,a 95%confidence interval11,12)better than±0.3%of the range of each respective component.

    The model pitch and roll orientations were measured with two Q- flex inclinometers.The inclinometers have an uncertainty of ±0.02°.The output of the pitch angle inclinometer was compensated for model vibrations by means of supporting signals from two linear accelerometers plus two angular rate accelerometers.13

    The 180 wing pressures were measured with three Electronically Scanned Pressure(ESP)modules with digital temperature compensation.The ESP modules have a range of±105 kPa and an uncertainty of±0.15%of their range.Besides the wing pressures,three more pressures were measured in the aft-fuselage cavity where the z-sting enters the model.Furthermore,two rows of 56 pressure taps were measured along the test section’s side wall using ESP modules with a range of±17.5 kPa and an uncertainty of±0.15%of their range.

    For the measurement of wing deformation,the so-called Stereo-Pattern Recognition(SPR)technique was used.SPR is an optical technique that makes use of a stereoscopic camera set-up to record images of markers on the wind tunnel model.For the AVM,ultraviolet fluorescent markers were applied on the wings and fuselages(Fig.5(a)).The markers were illuminated with several light emitting diodes to obtain highcontrast marker images,whilst keeping reflections to a minimum.The SPR technique allowed the position of an individual marker to be determined with an uncertainty of±0.1 mm.

    For the application of infrared thermography(IRT),14,15the port wing was provided with a coating.The coating provided sufficient insulation from conductivity effects,so that surface heat radiation due to aerodynamic phenomena(such as boundary layer transition and/or shock waves)could be detected with infrared cameras.Thermal images were acquired with two FLIR SC3000 cameras,set-up so as to view the top and bottom surfaces of the wing.A 0.1 mm coating was applied so that all pressure taps remained open.On top of the coating,three rows of silver painted tick marks were applied(one mark each 10%of the chord,as shown in Fig.5(b))merely for an easy chord-wise position reference in the raw thermal images.

    Fig.5 Views of the lower sides of the wings.

    3.4.Test conditions and data acquisition

    Primary database measurements for the BW and BWNVH model configurations were conducted on the ventral z-sting support.The measurements were conducted as sweeps of angle of attack at different Mach numbers,whilst the Reynolds number was kept constant.The Mach numbers ranged from 0.4 up to 0.9,and the angle of attack was varied between-4°≤α≤20°depending on the Mach number.The Reynolds number was kept constant at either 4.7×106or 2.0×106,and the total temperature ranged between 290 K and 310 K.

    The majority of the test program was conducted with forced laminar-turbulent transition tripping,using selfadhesive dots.The tripping dots were applied on all relevant surfaces.On the wing,dots with a height of 0.089 mm were applied at 7%chord on the upper surface and at 5%chord on the lower surface.On the horizontal tail,0.079 mm dots were used at 7%chord on both sides.Similarly,0.089 mm dots were used at 7%chord of the vertical tail.The fuselage was provided with 0.114 mm dots applied at 10 mm downstream from the nose,whilst the pylons and nacelle inner and outer surfaces had 0.089 mm dots located 15 mm downstream from their respective leading edges.The sizing of the dots was based on the method in Ref.16.

    Wind-on testing took place within runs,with one or more polars per run.In the ventral set-up,measurements were conducted in a continuous sweep test mode with an α-sweep rate of 0.1(°)/s and a data increment of 0.25°.For polars where wing deformation was measured,a step-by-step approach was used.During these polars,the flow conditions,balance loads,model attitude,wing surface and cavity pressures,wing deformation,thermal images,and wall pressures were all acquired in a parallel and synchronized mode.Fig.6 shows how various techniques were arranged,so as to allow all information to be acquired simultaneously during the primary database test.In addition,colored oil flow visualizations were conducted during separate runs,with the three views indicated in Fig.6.

    Fig.6 Schematic of the optical arrangement.

    In the dorsal set-up,only the model attitude,balance loads,and cavity pressures were measured.Wing surface pressures were not measured,as space requirements associated with the dorsal sting hampered the inclusion of the ESP modules in the narrow fuselage of the AVM.With the dummy zsupport installed,measurements were conducted in a stepby-step test mode,in the range from-2°≤α≤6°with a step of 0.5°.Without the dummy z-support,measurements were conducted in a continuous sweep testing mode again,in a range from-3°≤ α ≤ 7°with an increment of 0.25°.

    4.Data reduction and corrections

    4.1.General

    In the data reduction,the conversion of the measured quantities to relevant engineering units and dimensionless aerodynamic coefficients was handled.In the process,the following corrections were applied to account for fundamental wind tunnel interference effects:

    (1)Blockage correction: a Mach number correction based onresults of dedicated test section centerline calibrations.

    (2)Empty test section buoyancy correction:a correction on the drag to account for the effect of the empty test section’s axial pressure gradient.The buoyancy force is proportional to the product of the pressure gradient and the effective volume of the body.17For our purposes,the correction is formulated as

    where S is the reference wing area,?Cp/?x the centerline axial pressure gradient,A(x)the model cross-sectional area distribution,x the coordinate along the test section centerline,and xmthe length of the model.

    (3)Corrections for the in fluences of the z-support:corrections on the lift,drag,and pitching moment coefficients,as derived from dedicated dummy support measurements(Section 4.2).

    The model’s angle of attack is not specifically corrected for flow angle upwash;neither is a correction applied for the wall interference of the slotted test section.Regular reproduction measurements with a reference model have shown that upwash in the HST is within ±0.01°.As this is smaller than the uncertainty band of the applied inclinometer it was concluded that an upwash correction is not necessary for this test.As for the wall interference,Ref.9suggests that for the specific openness ratio of the slotted test section,no corrections are required for civil aircraft models with a reasonable span.The justification if wall interference is indeed also negligible for the AVM is presently hampered,due to the lack of wall pressure information along the top and bottom walls.

    4.2.Support interference

    As the support system generally imposes one of the largest influences on the test article aerodynamics,18–20the method for assessing and correcting this effect is addressed in this section in some length.Support interference is here defined as the influence of the primary z-support plus the influence of the empty test section with its permanently installed vertical strut.

    The effects due to the empty test section are accounted by means of the blockage and buoyancy corrections described in Section 4.1.The aerodynamic influence of the support is derived via two measurements on the dorsal set-up:one measurement with the model solely on the dorsal sting and the other measurement with the model mounted on the dorsal sting,whilst a dummy of the ventral z-support is present.For each test set-up,the reference flow conditions are set at the model MRP.The flow conditions at the MRP are known through calibrations of the static pressure distribution along the centerline of the test section.Such centerline calibrations are conducted for two situations: firstly,the situation with the straight boom present in the test section;secondly,the bare test section without the straight boom support present.To ensure identical flow conditions(i.e.,identical reference Mach number)between the two set-ups,ΔMa at the MRP due to the presence of the straight boom support is taken into account,on top of the Mach number due to the test section itself.The blockage due to the z-sting is neglected in this procedure,for the practical reason that the z-sting is not included in the centerline static pressure calibrations.The latter is deemed to be a reasonable approximation is deemed to be a reasonable assumption due to the small volume of the sting and its streamlined shape.

    With an identical reference Mach number set for both test set-ups,the difference in forces and moment coefficients is calculated by subtracting the balance results of the situation without the dummy z-support(DOR)from the results of the situation with the dummy z-support(DUM)at an identical angle of attack as follows:

    where X={L,D,m}for respectively the lift,drag,and pitching moment;ΔCXzrepresents the influence of the z-support.At this stage,the empty test section buoyancy correction of Eq.(1)is added,so as to end up with a correction ΔCXwhich covers all support interference effects as follows:

    Through this bookkeeping,ΔCXcomprises of the combined effect of the z-sting plus the straight boom(i.e.,the z-support)and the test section.The interaction between the z-support and the walls is hereby implicitly taken into account.Moreover,because it is verified that the pressure level in the sting-entry cavity is equal between the Fig.4(a)and(c)situations,it can be stated that the in fluences of the z-sting cavity and/or any model-internal cavity pressure effects are included in the support interference correction as well.Because the support system interferes differently for different model configurations,angles of attack,and Mach numbers,the corrections are determined separately for multiple combinations.Corrected results are finally calculated as

    where subscript SIC refers to the support interference corrected results,and subscript VEN refers to the uncorrected results on the Fig.4(a)ventral z-support.

    4.3.Wing deformation

    The SPR technique was used for assessment of the deformation of the wing under aerodynamic loading.Recorded stereo-images serve as input to derive the locations of model markers by means of a volumetric calibration,which links stereo-pairs of pixels to a unique position in a threedimensional space.The reconstruction of model markers in a three-dimensional space is done twice:once for a situation at a wind-off condition;the other for a situation at a wind-on condition.The wind-on results are aligned with the wind-off results by means of a subset of four reference markers on the wing root and fuselage.In turn,model deformation is defined as the difference between the wind-on and wind-off positions of the remaining markers in the model axis system.For each individual marker,the displacement in the z-direction is defined as its bending.To derive a value for the local twist,a pair of markers is used:one marker near the Leading Edge(LE)and the other near the Trailing Edge(TE).The vector connecting the LE and TE markers is determined twice:once for the wind-on situation;the other for the wind-off situation.Finally,the local twist is defined as the angular difference between these two vectors,projected in the zx-plane;nose-up twist is defined as positive.

    5.Results and discussion

    5.1.Tripping validation

    Infrared thermography was used to verify the effectiveness of the applied tripping dots to force laminar-turbulent transition on the wing.Fig.7 shows exemplary IRT results for a freetransition configuration(i.e.without tripping dots),as well as results for a configuration where the tripping dots were installed.

    Fig.7 IRT results for the upper wing surface at Ma=0.85.

    For the free-transition case in Fig.7(a),a clear temperature step can be seen over the observable part of the wing span.This temperature step reflects a sudden change in the local heat transfer coefficient as the boundary layer flow transitions from laminar to turbulent.For the free-transition case in Fig.7(b),again a clear temperature step is visible,but in this case,the temperature step is caused by the upper surface shock wave.Fig.7(c)and(d)shows that with the dots,the radiation images are notably different,as transition is forced at the locations of the dots.This can be deduced from Fig.7(c)where wedge-like shapes occur at the two stations where the dots were removed locally.The wedges represent an area with locally a lower heat transfer coefficient as is typical for laminar boundary layer conditions whilst directly behind the tripping dots,the heat transfer coefficient is larger,implying turbulent boundary layer conditions.Based on Fig.7,along with all other IRT data acquired during the test,it was concluded that the tripping dots are effective in forcing laminar-turbulent transition along the wing span at all Ma,α conditions considered.

    5.2.Data repeatability

    The repeatability of the measurement chain was inspected at various points along the test.This is done with two repeat polars measured within a single run,as well as with four repeat polars that were measured in different runs which had multiple model configurations changes between those runs.Fig.8 presents the lift coefficient versus the differences in the angle of attack Δα,drag coefficient ΔCD,and pitching moment coefficient ΔCm.The presented differences are obtained via linear interpolation of the data points from one polar to the identical lift coefficient value of the second repeat polar.These differences,or residuals,represent the level of variation between the repeat polars.The dashed lines shown on each plot indicate the 95%confidence interval,based on the applied instrumentation.Fig.8 shows that essentially all residuals are within the uncertainty band,which is in compliance with the expected standards for this type of test.18

    5.3.Support interference results

    Fig.8 Repeatability of the angle of attack,drag,and pitching moment coefficients at an identical(interpolated)lift coefficient at Ma=0.85(Dashed lines are the 95%confidence interval based on the applied instrumentation).

    Fig.9 Established support interference corrections for the BW and BWNVH configurations at Ma=0.85.(Error bars show the 95%confidence interval based on the applied instrumentation).

    The established support interferences for the BW and BWNVH configurations at the cruise Mach number are presented in Fig.9.The figure shows the interference levels according to Eq.(3);the error bars reflect the 95%confidence interval,based on the underlying instrumentation.The first observation is that the established interferences for the lift and pitching moment are relatively small for BW compared to BWNVH.Near the design lift condition(α ~ 2.5°),the established interferences are ΔCL=-0.005 and ΔCm=0.001 for BW,compared to ΔCL=-0.015 and ΔCm=0.025 for BWNVH.The signs and magnitudes of these ΔCmvalues imply that the support interferes mostly on the aft-fuselage,pylon-nacelles,and empennage.This is in-line with the expectation that the presence of the support causes the flow to stagnate progressively towards the back-end of the model,causing a decrease of the horizontal stabilizer effectiveness.The axial pressure gradient as induced by the test section and the z-support is further expected to cause a buoyancy force that alleviates the measured drag.The positive sign of the ΔCDcorrection for BW(Fig.9(b))confirms that the support is indeed pushing the model in the upstream direction.For the more complex BWNVH configuration,the ΔCDcorrection is initially also positive(i.e.,at a negative α),but changes smoothly towards negative drag interference levels as α increases.In fact,all interferences are found to vary to some extent with α,with a noticeable trend change near α ~ 2.5°.This change of trend is found to correlate with the point where the flow separation starts on the wing.

    Upon establishing the interference levels at several discrete angles of attack,the data is least-squared fitted with a polynomial so as to obtain a continuous correction function for the considered α-range.When correcting data at angles outside of this α-range(e.g., α≥ 5°),the value of the nearest border is maintained constant.

    Fig.10 SPR results for two repeat α-polars at Ma=0.85(Re=4.7×106,dynamic pressure q=54 kPa).

    5.4.Wing deformation

    An impression of the established twist and bending response along two α-polars is shown in Fig.10.For illustration purposes,the results of the two repeats are combined to generate a least-square fitted polynomial surface response for respectively bending and twist,as a function of the semi-span η and the angle of attack α.The resulting response surfaces show that the wing deforms elastically as a function of the angle of attack.Wing deformation is most pronounced at the tip,where the bending response varies from-9 mm to+18 mm along the α-polar,whilst the tip twist varies between+0.5°and-1.4°.At the design lift condition CL=0.5,the established tip twist is-1.0°,which complies with the initial twist requirement.In terms of uncertainty,the 95%confidence interval of the established twist result is determined as±0.05°in the wing root area and±0.2°at the wing tip.This is primarily due to data-reduction procedures described in Section 4.3,which causes the aforementioned SPR measurement accuracy of±0.1 mm for an individual marker to translate into a twist uncertainty that increases progressively towards the tip of the tapered wing.Because of the swept wing,one may expect that the established twist variation translates into a progressive loss of the total lift and an increase of the nose-up pitching moment along the α-polar.Such effects may be quantified using a CFD approach that considers the adapted geometry,according to the deformation results determined in this wind tunnel experiment.

    6.Conclusions

    Through design and wind tunnel testing of the CAE-AVM,an experimental dataset is established for two configurations of the aircraft.The overall loads and the wing load distributions are established for a variety of flow conditions between Mach numbers of 0.4 and 0.9.The repeatability of the overall loads is found to be within the expected uncertainty band of the underlying instrumentation.The interference imposed by the model support system is quantified on the level of the overall loads.The effectiveness of the applied wing transition trips is verified,and the amount of twist and bending of the wing is determined.Thereby,a comprehensive dataset is established which can be used for validation of computational methods.

    As part of potential future test work,additional wall pressure measurements can be considered for confirmation on the accuracy of the current dataset.

    1.European Commission.Flight Path 2050:Europe’s vision for aviation.Munich:Public Office of the European Union;2011.

    2.Poll DIA.Aerospace and the environment what contribution can aerodynamics make.The Aeronautical J 2000;104(1037):321–3.

    3.Takahashi TT.Optimum aspect ratio for subsonic air vehicles.J of Aircraft 2011;48(6):1984–93.

    4.Torenbeek E.Synthesis of subsonic aircraft design.Netherlands:Delft University Press;1982.p.155–66.

    5.Hua J,Zheng S,Zhong M,Wang GL,Chu L,Liu FL,et al.Design and verification study of an aerodynamic validation model.7th Asia-Pacific international symposium on aerospace technology;2015.

    6.Levy DW,La flin KR,Vassberg JC,Tinoco EN,Mani M,Rider B,et al.Summary of data from the fifth AIAA CFD Drag Prediction Workshop.Reston:AIAA;2013.Report No.:AIAA-2013-0046.

    7.Eberhardt S,Benedict K,Hedges L,Robinson A,Tinoco EN.Inclusion ofaeroelastic twist into the CFD analysis of the twinengine NASA common research model.Reston:AIAA;2014.Report No.:AIAA-2014-0251.

    8.Keye S,Brodersen O,Rivers MB.Investigation of aeroelastic effects on the NASA common research model.J of Aircraft 2014;51(4):1323–30.

    9.Maseland JEJ,Laban M,Ven H,van der Kooi JW.Development of CFD-based interference models for the DNW-HST transonic wind-tunnel.Reston:AIAA;2006.Report No.:AIAA-2006-3639.

    10.Philipsen I,Poel M,Elsenaar A.DNW-HST(High Speed Tunnel)50-year anniversary.Reston:AIAA;2006.Report No.:AIAA-2006-3639.

    11.Bergmann R,Philipsen I.Some contemplations on a proposed definition of uncertainty for balances.Reston:AIAA;2010.Report No.:AIAA-2010-4546.

    12.Coleman HW,Steel WG.Experimentation and uncertainty analysis for engineers.2nd ed.New York:Wiley&Sons;1989.p.24–8.

    13.Fuijkschot PH.Vibration compensation of gravity sensing inclinometers in wind tunnel models.Proceedings ISA 42nd international instrumentation symposium;1996.p.493–503.

    14.Kompenhans J,Agocs J,Egami Y,Engler R,Fey U,Frahnert H,et al.Recent developments of image based measurement methods for application to transonic flows in industrial wind tunnels.Chin J Aeronaut 2006;19(2):115–25.

    15.Crawford BK,Duncan Jr GT,West DE,Saric WS.Laminarturbulent boundary layer transition imaging using IR thermography.Optics and Photonics J 2013;3(3):233–9.

    16.Braslow AL,Knox EC.Simplified method for determination of critical height of distributed roughness for boundary layer transition at mach numbers 0 to 5.Langley Field:National Advisory Committee for Aeronautics,Langley Aeronautical Lab.;1958.Report No.:NACA-TN-4363.

    17.Swanson RS,Gillis,CL.Wind-tunnel calibration and correction procedures for three-dimensional models.Langley Field:National Advisory Committee for Aeronautics,Langley Aeronautical Lab.;1944.Report No.:NACA-WR-L-1.

    18.Wubben FJM,Takara E.Wind tunnel model support and wall interference corrections in DNW-HST-ensuring high data quality standards.5th CEAS air&space conference;2015.

    19.Elsenaar A,Han SOTH.A break-down of sting interference effects.Model support corrections in wind tunnels symposium;1991.

    20.Rivers MB,Hunter CA,Campbell RL.Further investigation of the support system effects and wing twist on the NASA common research model.Reston:AIAA;2012.Report No.:AIAA-2012-3209.

    乱系列少妇在线播放| 中文字幕亚洲精品专区| 国产国拍精品亚洲av在线观看| 久久久国产精品麻豆| 午夜影院在线不卡| 少妇人妻一区二区三区视频| 卡戴珊不雅视频在线播放| 久久久久久久大尺度免费视频| 亚洲成人av在线免费| 精品一区二区免费观看| 毛片一级片免费看久久久久| 久久久国产欧美日韩av| 在线观看国产h片| 丁香六月天网| 六月丁香七月| 夫妻午夜视频| 精品午夜福利在线看| 成人国产麻豆网| 久久久午夜欧美精品| 亚洲电影在线观看av| 岛国毛片在线播放| 午夜福利,免费看| 男人狂女人下面高潮的视频| 国产真实伦视频高清在线观看| 大码成人一级视频| 能在线免费看毛片的网站| 自拍欧美九色日韩亚洲蝌蚪91 | 国产免费又黄又爽又色| 日韩欧美精品免费久久| 久久国产乱子免费精品| 人人妻人人看人人澡| 久久人人爽人人片av| 日韩强制内射视频| 国产永久视频网站| 中文资源天堂在线| 青青草视频在线视频观看| 免费人妻精品一区二区三区视频| 亚洲精品国产av蜜桃| videossex国产| 亚洲av福利一区| 成人黄色视频免费在线看| 另类亚洲欧美激情| 国产精品嫩草影院av在线观看| 亚洲一级一片aⅴ在线观看| 国产精品.久久久| 国产乱人偷精品视频| 少妇猛男粗大的猛烈进出视频| 亚洲国产av新网站| 日韩一本色道免费dvd| 伦理电影大哥的女人| 在线观看一区二区三区激情| 香蕉精品网在线| 日本爱情动作片www.在线观看| 精品一区二区免费观看| 国产av精品麻豆| 亚洲国产精品成人久久小说| 欧美精品一区二区大全| 久久久久网色| 亚洲精品,欧美精品| 久久精品国产亚洲av天美| 少妇人妻一区二区三区视频| 精品国产一区二区三区久久久樱花| 亚洲不卡免费看| 妹子高潮喷水视频| 一本久久精品| 欧美性感艳星| 国内揄拍国产精品人妻在线| 美女大奶头黄色视频| 蜜桃在线观看..| 久久精品夜色国产| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说| 亚洲av成人精品一区久久| 97超视频在线观看视频| 国产av一区二区精品久久| 少妇 在线观看| 国产高清国产精品国产三级| 3wmmmm亚洲av在线观看| 免费大片黄手机在线观看| 欧美97在线视频| 亚洲国产精品国产精品| 丝袜在线中文字幕| 男男h啪啪无遮挡| 久久久久精品久久久久真实原创| 欧美日韩一区二区视频在线观看视频在线| av在线老鸭窝| 99久久精品国产国产毛片| 亚洲自偷自拍三级| 亚洲欧美成人综合另类久久久| 99re6热这里在线精品视频| 国产爽快片一区二区三区| 免费黄色在线免费观看| 啦啦啦中文免费视频观看日本| 亚洲熟女精品中文字幕| 永久网站在线| 欧美最新免费一区二区三区| 国产 精品1| 99热网站在线观看| 噜噜噜噜噜久久久久久91| 亚洲av成人精品一二三区| 亚洲欧美成人综合另类久久久| 国产成人91sexporn| 精品国产一区二区久久| 99久久精品热视频| 精品久久久久久电影网| 少妇人妻久久综合中文| 精品国产一区二区久久| 国产伦精品一区二区三区视频9| 日本wwww免费看| 在线观看免费日韩欧美大片 | 欧美精品高潮呻吟av久久| 精品人妻熟女毛片av久久网站| 免费高清在线观看视频在线观看| 一级二级三级毛片免费看| 乱系列少妇在线播放| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| 日本色播在线视频| 日本av手机在线免费观看| 99热6这里只有精品| 亚洲三级黄色毛片| 永久网站在线| 中文精品一卡2卡3卡4更新| 一区二区三区乱码不卡18| 国产91av在线免费观看| 亚洲经典国产精华液单| 成人黄色视频免费在线看| 一二三四中文在线观看免费高清| 国产欧美日韩一区二区三区在线 | 国产免费福利视频在线观看| 亚洲国产欧美在线一区| 老司机影院成人| a级一级毛片免费在线观看| 亚洲第一区二区三区不卡| 99久久中文字幕三级久久日本| 毛片一级片免费看久久久久| 精品亚洲乱码少妇综合久久| 亚洲经典国产精华液单| 美女大奶头黄色视频| 2018国产大陆天天弄谢| 国产中年淑女户外野战色| 男女无遮挡免费网站观看| 不卡视频在线观看欧美| 亚洲精品乱码久久久久久按摩| 中国三级夫妇交换| 国产亚洲精品久久久com| 乱码一卡2卡4卡精品| .国产精品久久| 成人国产麻豆网| 高清不卡的av网站| 精品亚洲成a人片在线观看| 波野结衣二区三区在线| xxx大片免费视频| 人人妻人人添人人爽欧美一区卜| 婷婷色av中文字幕| 久久久久久久久久久免费av| 精品久久久久久电影网| 一本大道久久a久久精品| 熟妇人妻不卡中文字幕| 精品熟女少妇av免费看| 久久人人爽人人爽人人片va| 日本欧美国产在线视频| 一级毛片aaaaaa免费看小| 国产高清不卡午夜福利| 最近最新中文字幕免费大全7| 嫩草影院入口| 乱人伦中国视频| 国产白丝娇喘喷水9色精品| 国产极品粉嫩免费观看在线 | 爱豆传媒免费全集在线观看| 99久国产av精品国产电影| 在线观看人妻少妇| 两个人的视频大全免费| 精品久久久久久久久av| 亚洲四区av| 日本91视频免费播放| 你懂的网址亚洲精品在线观看| 色婷婷av一区二区三区视频| 国内精品宾馆在线| 夜夜爽夜夜爽视频| av网站免费在线观看视频| 男女边吃奶边做爰视频| 国产亚洲一区二区精品| 国产成人精品福利久久| 777米奇影视久久| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 日韩电影二区| 亚洲精品,欧美精品| av天堂久久9| 嫩草影院入口| 久久婷婷青草| 国产精品久久久久久久电影| xxx大片免费视频| 丝袜脚勾引网站| 国产免费又黄又爽又色| 午夜日本视频在线| 日韩大片免费观看网站| 亚洲精品第二区| 在线播放无遮挡| 国产精品福利在线免费观看| 久久人人爽av亚洲精品天堂| 三上悠亚av全集在线观看 | 十分钟在线观看高清视频www | 亚洲人成网站在线播| 五月开心婷婷网| 少妇人妻 视频| 99久久综合免费| 丁香六月天网| 亚洲精品亚洲一区二区| 国产淫语在线视频| av.在线天堂| 日韩不卡一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 亚洲第一av免费看| 中文欧美无线码| 久久青草综合色| av在线老鸭窝| 亚洲av国产av综合av卡| av卡一久久| 国产视频内射| 欧美最新免费一区二区三区| 五月玫瑰六月丁香| 一级av片app| 国产成人一区二区在线| 亚洲综合色惰| 一级毛片我不卡| 国产免费又黄又爽又色| 国产一区二区三区av在线| 久久精品久久精品一区二区三区| 少妇精品久久久久久久| 亚洲熟女精品中文字幕| 精品国产国语对白av| 国产免费又黄又爽又色| 免费人妻精品一区二区三区视频| 毛片一级片免费看久久久久| 乱码一卡2卡4卡精品| 人妻夜夜爽99麻豆av| 国国产精品蜜臀av免费| 午夜精品国产一区二区电影| 高清黄色对白视频在线免费看 | 久久久国产欧美日韩av| 日本-黄色视频高清免费观看| 国产国拍精品亚洲av在线观看| a级毛色黄片| 五月开心婷婷网| 亚洲精品色激情综合| 亚洲欧美一区二区三区国产| 成年女人在线观看亚洲视频| 一区二区三区精品91| 午夜免费鲁丝| 中文字幕免费在线视频6| 夫妻性生交免费视频一级片| 久久久精品94久久精品| 在线观看www视频免费| 少妇高潮的动态图| 十分钟在线观看高清视频www | 亚洲第一av免费看| 男人舔奶头视频| 观看美女的网站| 一级黄片播放器| 高清视频免费观看一区二区| 亚洲美女视频黄频| 人妻少妇偷人精品九色| 九草在线视频观看| 啦啦啦在线观看免费高清www| 久久人人爽人人片av| 国产成人精品一,二区| 男女边摸边吃奶| 青青草视频在线视频观看| 蜜桃久久精品国产亚洲av| 国产有黄有色有爽视频| 久久午夜综合久久蜜桃| 精品亚洲乱码少妇综合久久| 亚洲精品亚洲一区二区| 一区二区三区乱码不卡18| h视频一区二区三区| av播播在线观看一区| 精品少妇内射三级| 成人综合一区亚洲| 久久国产乱子免费精品| 能在线免费看毛片的网站| av专区在线播放| 国产一区二区三区av在线| 精品卡一卡二卡四卡免费| 精品一品国产午夜福利视频| 久久精品国产自在天天线| 2021少妇久久久久久久久久久| 少妇的逼水好多| 三级经典国产精品| 狂野欧美激情性xxxx在线观看| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 成人漫画全彩无遮挡| 中文字幕亚洲精品专区| 男的添女的下面高潮视频| 视频中文字幕在线观看| 伊人久久国产一区二区| 亚洲av中文av极速乱| 亚洲,一卡二卡三卡| 91成人精品电影| 亚洲av日韩在线播放| 免费大片18禁| 涩涩av久久男人的天堂| 日韩人妻高清精品专区| 嫩草影院新地址| 国产日韩欧美亚洲二区| 久久国产乱子免费精品| 下体分泌物呈黄色| 欧美激情极品国产一区二区三区 | 欧美97在线视频| 成人毛片a级毛片在线播放| 美女xxoo啪啪120秒动态图| 国产在线一区二区三区精| 欧美精品一区二区大全| 日本猛色少妇xxxxx猛交久久| 97精品久久久久久久久久精品| 欧美激情极品国产一区二区三区 | 最近中文字幕2019免费版| 午夜久久久在线观看| 欧美日韩亚洲高清精品| a级毛片在线看网站| 18禁在线无遮挡免费观看视频| 欧美一级a爱片免费观看看| 自拍偷自拍亚洲精品老妇| 91在线精品国自产拍蜜月| 美女视频免费永久观看网站| 你懂的网址亚洲精品在线观看| 激情五月婷婷亚洲| 免费黄频网站在线观看国产| 亚洲精品亚洲一区二区| 亚洲成人手机| 国产成人精品一,二区| 亚洲一区二区三区欧美精品| 一级毛片我不卡| kizo精华| 亚洲国产欧美在线一区| 国产在线一区二区三区精| 亚洲真实伦在线观看| 亚洲综合精品二区| 久久av网站| 极品人妻少妇av视频| 日韩一区二区三区影片| 亚洲欧美一区二区三区黑人 | 中文字幕制服av| av网站免费在线观看视频| 亚洲成人手机| 午夜视频国产福利| 精品一区在线观看国产| 成人午夜精彩视频在线观看| 草草在线视频免费看| 久久久久国产精品人妻一区二区| 国产精品久久久久久久电影| 国产黄片美女视频| 嫩草影院入口| 夫妻性生交免费视频一级片| 日日爽夜夜爽网站| 久久亚洲国产成人精品v| 丰满人妻一区二区三区视频av| 99久久人妻综合| freevideosex欧美| 曰老女人黄片| 亚洲av欧美aⅴ国产| 在线观看免费高清a一片| 涩涩av久久男人的天堂| 看非洲黑人一级黄片| 国产精品伦人一区二区| 久久久久久久久久久久大奶| 乱系列少妇在线播放| 精品国产一区二区三区久久久樱花| 综合色丁香网| av视频免费观看在线观看| 最黄视频免费看| 欧美丝袜亚洲另类| 22中文网久久字幕| videossex国产| 亚洲成色77777| 亚洲av电影在线观看一区二区三区| 日本91视频免费播放| 国产一区二区在线观看日韩| 如何舔出高潮| 国产极品天堂在线| 丝袜脚勾引网站| 国产av码专区亚洲av| 最近手机中文字幕大全| 免费大片黄手机在线观看| 国国产精品蜜臀av免费| av线在线观看网站| 国产 一区精品| av不卡在线播放| 免费观看在线日韩| 美女内射精品一级片tv| 性高湖久久久久久久久免费观看| 精品少妇久久久久久888优播| 欧美变态另类bdsm刘玥| 免费久久久久久久精品成人欧美视频 | 亚洲精品乱码久久久v下载方式| 国产精品福利在线免费观看| 麻豆精品久久久久久蜜桃| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久| 丝袜喷水一区| 亚洲欧洲国产日韩| 国产欧美日韩一区二区三区在线 | 日日啪夜夜撸| 91精品国产国语对白视频| 精品少妇黑人巨大在线播放| 国产欧美另类精品又又久久亚洲欧美| 亚洲av不卡在线观看| 国产黄频视频在线观看| 熟女人妻精品中文字幕| 国产永久视频网站| 十八禁高潮呻吟视频 | 一级毛片久久久久久久久女| 精品国产露脸久久av麻豆| 久久久久国产网址| 成人18禁高潮啪啪吃奶动态图 | 热re99久久精品国产66热6| 国产免费视频播放在线视频| 久久午夜综合久久蜜桃| 国产熟女欧美一区二区| 国产黄片美女视频| 最新中文字幕久久久久| 国产视频首页在线观看| 日韩欧美精品免费久久| 乱系列少妇在线播放| 一级爰片在线观看| 久久精品国产亚洲av涩爱| 中文天堂在线官网| 桃花免费在线播放| 有码 亚洲区| 国语对白做爰xxxⅹ性视频网站| 女的被弄到高潮叫床怎么办| av有码第一页| 欧美精品一区二区大全| 国产高清国产精品国产三级| 亚洲精品日本国产第一区| 少妇被粗大的猛进出69影院 | 精品视频人人做人人爽| 女人久久www免费人成看片| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 高清毛片免费看| 少妇的逼水好多| 欧美日韩视频精品一区| 免费久久久久久久精品成人欧美视频 | 亚洲精品日本国产第一区| a级毛色黄片| 国产色爽女视频免费观看| 国产欧美日韩精品一区二区| 成人免费观看视频高清| 国产精品无大码| 精品久久久久久久久av| av卡一久久| 久久久久网色| 精品久久久久久久久av| 免费大片18禁| 最近2019中文字幕mv第一页| 国产成人精品一,二区| 国产精品麻豆人妻色哟哟久久| 日本av手机在线免费观看| 日本wwww免费看| 97在线视频观看| 能在线免费看毛片的网站| 国产欧美日韩综合在线一区二区 | 久久鲁丝午夜福利片| 99视频精品全部免费 在线| 国产真实伦视频高清在线观看| 中国国产av一级| 国产乱来视频区| 乱系列少妇在线播放| 亚洲国产毛片av蜜桃av| 99视频精品全部免费 在线| 日本午夜av视频| 国产亚洲最大av| 多毛熟女@视频| 少妇裸体淫交视频免费看高清| 久久久久精品性色| 天堂俺去俺来也www色官网| 韩国高清视频一区二区三区| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 在线天堂最新版资源| 国产高清不卡午夜福利| 高清黄色对白视频在线免费看 | 我要看日韩黄色一级片| 午夜福利影视在线免费观看| 黄色毛片三级朝国网站 | 少妇 在线观看| 久久99热6这里只有精品| 成人美女网站在线观看视频| 国产有黄有色有爽视频| 99久久中文字幕三级久久日本| 一级二级三级毛片免费看| 97精品久久久久久久久久精品| 久久婷婷青草| 色5月婷婷丁香| 亚洲欧美日韩东京热| 91久久精品电影网| 我的女老师完整版在线观看| 精品国产国语对白av| 九九爱精品视频在线观看| 国产在线免费精品| 国产成人一区二区在线| freevideosex欧美| av天堂久久9| av国产精品久久久久影院| av线在线观看网站| h日本视频在线播放| 亚洲熟女精品中文字幕| 成人漫画全彩无遮挡| 大片免费播放器 马上看| 中文字幕亚洲精品专区| 久久久国产精品麻豆| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产色片| 另类亚洲欧美激情| 边亲边吃奶的免费视频| 日日啪夜夜撸| 国产成人精品福利久久| av福利片在线| 国产成人精品婷婷| 热99国产精品久久久久久7| 国产免费又黄又爽又色| 国产在线男女| 一级,二级,三级黄色视频| av免费在线看不卡| 嫩草影院新地址| 日韩 亚洲 欧美在线| 十八禁高潮呻吟视频 | 久久狼人影院| 蜜桃在线观看..| 午夜av观看不卡| 人人妻人人添人人爽欧美一区卜| 国产在视频线精品| 丝袜喷水一区| 国产精品国产三级专区第一集| 久久久久久久精品精品| 免费看日本二区| 国产成人精品久久久久久| 最新的欧美精品一区二区| 97在线视频观看| 性色av一级| 看十八女毛片水多多多| 亚洲av成人精品一二三区| 永久免费av网站大全| 日韩伦理黄色片| 国产成人精品婷婷| 日韩一区二区三区影片| 一个人看视频在线观看www免费| 精品国产一区二区久久| 一本久久精品| av在线老鸭窝| 蜜臀久久99精品久久宅男| 卡戴珊不雅视频在线播放| 国产无遮挡羞羞视频在线观看| 国产毛片在线视频| av免费观看日本| 精品久久久久久久久av| 国产视频内射| 亚洲久久久国产精品| 搡女人真爽免费视频火全软件| 一本一本综合久久| 免费看av在线观看网站| 亚洲经典国产精华液单| 国产午夜精品久久久久久一区二区三区| 精品卡一卡二卡四卡免费| 中文字幕制服av| 国产欧美亚洲国产| 国产精品三级大全| 一级av片app| 高清黄色对白视频在线免费看 | 肉色欧美久久久久久久蜜桃| 男女啪啪激烈高潮av片| 最近中文字幕高清免费大全6| 国产成人午夜福利电影在线观看| 亚洲天堂av无毛| 五月天丁香电影| 国产成人精品久久久久久| 最近中文字幕高清免费大全6| 欧美 日韩 精品 国产| 制服丝袜香蕉在线| 国产欧美日韩综合在线一区二区 | 成年人午夜在线观看视频| 蜜桃久久精品国产亚洲av| 国产熟女午夜一区二区三区 | 日本色播在线视频| 最近中文字幕高清免费大全6| 看免费成人av毛片| 噜噜噜噜噜久久久久久91| 国产精品国产av在线观看| 亚洲怡红院男人天堂| 国产精品99久久99久久久不卡 | 秋霞在线观看毛片| av有码第一页| 能在线免费看毛片的网站| 亚洲国产日韩一区二区| 久久国产亚洲av麻豆专区| 人妻一区二区av| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看 | 国产视频首页在线观看| 18禁动态无遮挡网站| 最黄视频免费看| 国产69精品久久久久777片| 美女主播在线视频| 国产精品一区www在线观看| a级毛片在线看网站| 久久午夜福利片| 精品一品国产午夜福利视频| 男女边吃奶边做爰视频| 一级av片app| 哪个播放器可以免费观看大片| 亚洲欧美日韩另类电影网站| 亚洲精品日韩在线中文字幕|