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

    Experimental validation of a new morphing trailing edge system using Price - Pa?¨doussis wind tunnel tests

    2019-07-01 07:42:30COMMUNIERBOTEZWONG
    CHINESE JOURNAL OF AERONAUTICS 2019年6期

    D. COMMUNIER, R.M. BOTEZ, T. WONG

    Laboratory of Applied Research in Active Controls, Avionics and AeroServoElasticity (LARCASE), ETS, Montre′al H3C 1K3, Canada

    KEYWORDS Aerodynamic performances;Aileron;Morphing wing;Trailing edge;Wing tunnel tests

    Abstract This paper presents the design and manufacturing of a new morphing wing system carried out at the Laboratory of Applied Research in Active Controls, Avionics and AeroServoElasticity (LARCASE) at the ETS in Montre′al. This f irst version of a morphing wing allows the deformation of its trailing edge,denote by Morphing Trailing Edge(MTE).In order to characterize the technical impact of this deformation,we compare its performance with that of a rigid aileron by testing in the LARCASE's price—Pa?¨doussis subsonic wind tunnel. The f irst set of results shows that it is possible to replace an aileron by a MTE on a wing, as an improvement was observed for the MTE aerodynamic performances with respect to the aileron aerodynamic performances.The improvement consisted in the fact that the drag coeff icient was smaller, and the lift-to-drag ratio was higher for the same lift coeff icient.

    1. Introduction

    Due to its monitoring function, an Unmanned Aerial Vehicle(UAV) makes irregular f lights, in which, the f light conditions change frequently and signif icantly. For a conventional aircraft,the surface of the wing is f ixed.That is,airfoils representing wing shapes are chosen and utilized throughout the aircraft's design.As the aircraft passes through different f light phases(climb,descent,and cruise),a compromise is considered to obtain the best eff iciency of the global f light. If this eff iciency is compared to the optimal eff iciency for each f light phase, the eff iciency of each phase will be less optimal than the total eff iciency. This loss of eff iciency would results in higher fuel consumption and/or a shorter operating time.1-3

    The f irst objective is to modify the geometrical shape of the wing according to f light conditions with the aim to reach as much as possible its optimal shape for each f light condition.The morphing wing will be designed to increase its aerodynamic performance by increasing lift-drag ratio, which therefore would be equivalent to the increase of the lift and reduction of the drag.This fact will have the effect of reducing fuel consumption, increasing f light autonomy, etc.

    Nomenclature

    CLlift coeff icient

    CDdrag coeff icient

    e airfoil thickness

    t thickness of the slit

    L distance between slit and trailing edge

    yttotal displacement

    y displacement of the trailing edge

    p depth of slit

    l width of slit

    α Angle of the morphing trailing edge

    n number of slits

    The second objective is to keep the wing weight less than or equal to its current weight. To perform this objective, the structure of the current wing should be analyzed to establish its design criteria (weight, f lexural strength and torsion, maximum permissible load), and to further design a new wing that meets these design criteria.

    The morphing wings are studied at our Laboratory of Applied Research in Active Controls, Avionics and AeroServoElasticity (LARCASE) at the ETS since 2003.Indeed,this team has already explored numerous development strategies for morphing wing technologies. This research was carried out as part of two major projects of the Consortium of Research and Innovation in Aerospace in Quebec(CRIAQ).CRIAQ projects are carried out within the framework of signed agreements between universities, industrial partners and research institutes. Each of these projects has led to the designing and manufacturing of prototypes and their experimental wind tunnel tests, and to publications of results in scientif ic journals and conferences. These projects, CRIAQ 7.1 and CRIAQ MDO 505,were carried out in collaboration with Bombardier, Thales, NRC-IAR and E′ cole Polytechnique in Canada. In addition, the CRIAQ MDO 505 project was realized in collaboration with Italian partners,Alenia,CIRA,and the University of Naples—Frederico II.These two studies carried out under the leadership of LARCASE team aimed at deforming the upper surface of the morphing wing in order to improve its aerodynamic performance. The f irst project was called CRIAQ 7.1, in which the deformation of the upper surface was applied to a wing using ‘‘Smart Material Actuators” (SMA).4,5The resulting deformation made it possible to delay the f low transition on the wing. An ideal rectangular‘‘wing” was used, and therefore no existing structural constraints for a real wing were considered. This concept was experimentally validated in the NRC-IAR wind tunnel. The second CRIAQ project was called the CRIAQ MDO 505 project,6-8where MDO stands for ‘‘Multidisciplinary Design Optimization”. The objective was to design and validate a wing with an aileron. For that wing, it was necessary to move the passage of the air f low along the chord of the wing in order to delay the transition from the laminar to the turbulent f low.In order to perform this motion, the wing tip was equipped with four punctual actuators to change its aileron shape design for a real regional Bombardier aircraft that exhibited structural constraints.The experimental results showed that the system could produce the delay of the airf low transition. A subsequent project with an airfoil of an ATR 42 (Regional Transport Aircraft) made it possible to modify the morphing wing upper surface by means of two oval spars.9When those spars rotated, the surface was pushed outwards. This system therefore makes it possible to change the surface outward but not inward,and therefore the improvement of the aerodynamic performances was limited. These deformation methods allowed modif ication of the airf low,but they proved to be limited in modifying the coeff icients values for aerodynamic performance improvement on the morphing ATR-42 wing (less than 5% for the lift and less than 2% for the drag7). From these three projects,in order to signif icantly improve the aerodynamic coeff icients of a wing10(for example, an increase of more than 50%of the lift with an increase of lift-to-drag ratio L/D), the deformation of the upper surface alone was not suff icient. Therefore, in this work, the design and experimental validation of a morphing wing trailing edge system are analyzed.

    Concerning the state of the art in wing morphing, we can refer to the work of Sof la et al.11

    2. Design of the morphing wing system

    2.1. Deformation of the camber

    Studies have shown that the most promising type of deformation to meet this study objectives would be the deformation of an airfoil's camber.12Indeed, by modifying the camber of an airfoil, it could be possible, for example, to convert a NACA0012 airfoil into a NACA4412 airfoil,thereby attaining a signif icant increase of the lift of the wing while incurring a smaller increase of the drag, thus resulting in a higher L/D.National Advisory Committee for Aeronautics (NACA) airfoils, whose names consist of two digits and a two-digit number (convention for NACA airfoil with 4 digits), follow the following standard: the f irst digit denotes the maximum camber of the airfoil as a percentage of the chord;the second digit gives the position of this maximum camber in tens of percent of the chord; the third number (the last two digits) represents the thickness of the airfoil at a percentage of the chord.13

    Two types of deformations, combined together, made possible to change an airfoil's camber;and these deformations are:the deformation of the leading edge and the deformation of the trailing edge.14However, an often encountered problem when the shape of the airfoil deforms consists in the occurrence of waves on the surface of the wing dues to the interactions between the actuating system and the wing structure.15

    Another diff iculty encountered in the development of a morphing wing mechanism for reduced sized aircraft is its weight and its bulk. Any gain in aerodynamic performance(lift,drag)could thus be cancelled out by the additional weight of the mechanism. This is our motivation for developing a mechanism that can improve the aerodynamic performance of a morphing wing without increasing its weight.

    2.2. Deformation of the trailing edge

    A morphing aileron system was developed as part of the CRIAQ MDO 505 project.16This system made it possible to operate an articulated aileron from an arm that was rotated,which allowed the shape of the aileron to be changed without the need for a mechanism external to the aileron. Research in Germany at Braunschweig presented a concept of the deformation of the ribs at the trailing edge of the wing using an articulated skeleton.17Their report concluded that using the deformation of the trailing edge,it would be possible to replace the ailerons at the end of the wing, but that this deformation alone would not be suff icient to replace the f laps on the wing.

    Many research studies have been carried out on the morphing trailing edge, using different actuators types (mechanical,piezoelectric, smart materials actuators). The mechanical deformations were mainly done using articulated ribs like f ingers18,19using a morphing structure20,21(compliant mechanisms). These mechanisms made possible to control the deformation of the structure, but the surface of the wing did not completely respect the desired deformations (wave formations on the surface).In addition,the necessary actuators were often heavy and bulky.The use of smart material actuators22-24made it possible to replace mechanical actuators. This reduced the weight of the system but required greater power consumption in return.Another method considered the acting on the surface directly using piezoelectric materials.25,26However,these materials only allowed small displacements,and required large electrical consumption.

    In our work, to obtain a deformation of the trailing edge rib, vertical incisions were made. Depending on the length,the number and the width of these incisions, it is possible to modify the f lexibility of the rib and the amplitude of the deformation.To control the deformation of the rib,a servomotor is used that acts directly on the trailing edge as shown in Fig. 1.

    In order to analyze the aerodynamics of the deformed rib,it is f irstly necessary to compute its structural deformation using a Finite Element Analysis (FEA) modelling in CATIA V5 as illustrated in Fig. 2.

    From the results of the Finite Element Analysis, the contour of the deformed rib was extracted using the PROFSCAN tool.27PROFSCAN can draw a curve from an image to create a ‘‘.DAT” f ile that is used further in the XFLR5 software(Fig.3).To compare the aerodynamic eff iciency of the two airfoils,a wing was def ined with the same geometrical dimensions as the reference wing previously studied at the LARCASE28:10 in (254 mm) for the chord and 11.5 in (292.1 mm) for the span. This wing was analyzed for the inclined, rigid aileron,and then, for the Morphing Trailing Edge (MTE), both with a vertical displacement of trailing edge corresponding at 6.8% of the chord.

    Fig. 1 Control of the trailing edge rib deformation.

    Fig. 2 Deformed rib in CATIA V5.

    Fig. 3 NACA0012 airfoil with inclined aileron and MTE.

    The curves presented in Fig.4 compare the performance of a wing with an inclined aileron and that of a wing with a MTE for the airfoils presented in Fig.3(a),and in Fig.3(b).The calculations were carried out for a speed of 20 m/s (65.62 ft/s)using the 3D Panels method of the XFLR5 code.29The 3D Panels method calculates the aerodynamic f low around a wing quickly. This method alone does not allow to obtain the viscous f low characteristics (parameters), but by using the XFoil analysis on the wing airfoils (2D), an extrapolation was made to calculate the viscous f low for the wing (3D). XFLR5 proposes another method of calculation, but we prefer this method because it allow us to recover the distribution of the coeff icients of pressure on the upper and lower surfaces of the wing in order to carry out an aero structural analysis on CATIA V5.28The calculation method does not allow for the aerodynamic coeff icients to be computed after the stall of the air f low, and for this reason, the calculations stop at the angle of 14°.

    Regarding the drag coeff icient reduction, the MTE generates less drag than the inclination of the aileron for the same displacement, which meets our goal. This decrease in drag comes mainly from the fact that the deformation of the camber takes place on a larger portion of the wing chord than the deformation of the camber induced by the aileron. The size of the aileron was chosen accordingly to its required geometry by the aircraft model (~25% of the chord), and the length of deformation according to the chord of the MTE was chosen to be close to a NACA4412 airfoil chord, so that a variation of camber from 40% of the chord to the leading edge was obtained.

    The wing with an inclined aileron generates more lift than the MTE for their same vertical displacement of the trailing edge of 6.8% of the chord. Wing with inclined aileron generates more lift but they incur more drag,while the MTE generates less drag but also less lift. Therefore, it is diff icult to determine which system is the most eff icient from the point of view of these two forces. However, the variation of the L/D with the angle of attack is traced in Fig. 4(c) to indicate which system would have the best L/D.

    Fig. 4 Comparation of the performance of a wing with inclined aileron and MTE.

    The L/D variation with the angle of attack presented in Fig.4(c)shows that the MTE has a better L/D than an inclined aileron for angles of attack higher than-5°;thus, the MTE is more eff icient than the inclined aileron for positive angles of attack.These f irst theoretical results encouraged us to continue in this direction in order to obtain experimental results.

    In the next section,we will discuss on the ways in which the slits in the ribs were sized through design to obtain the MTE that was manufactured using wood.

    3. Sizing design and control of the MTE

    3.1. Sizing design

    To achieve rib f lexibility,the anisotropic property of wood was used.Using the f ibers in the direction of the chord of the airfoil and making a slit section on the thickness of the airfoil, the bending points were created according to the characteristics of the ‘‘compliant mechanisms”.30These bending points behave like ‘‘pivot points”31by allowing the simplif ication of the calculation of the maximum displacement of the MTE.The advantage of this design using bending pivots with respect to a conventional design with mechanical pivots is that there is no need for a mechanism to allow the rotation of the MTE.This fact reduces the weight of the assembly, and facilitates wing maintenance. In order to reduce the force required to bend the rib, the slit (see Fig. 1) must be designed as deep as possible, while keeping enough material for the rib structure to withstand tangential stresses.In order to f ind out the values of these stresses,the aerodynamic forces were f irstly computed on the MTE.

    For the f irst prototype, no sizing calculations were made,and the already acquired experience in wooden wing design and fabrication was suff icient for its manufacturing. In addition, to facilitate the prototype design, the symmetrical slits were dimensioned with symmetrical airfoils that allowed the MTE to change its shape by moving upwards and downwards.The slit was considered as a pivot point, and its displacement for a given width in the rib was calculated.Fig.5 illustrates the parameters required for the rib deformation calculation.

    From Fig.5,the following equations are obtained as the rib is symmetrical:

    The displacement y depends on the dimensions of the slit according to the next equation:

    To obtain this equation, the slits were considered as pivot and so the equation correspond to an angle return.

    As we have n slits in the rib, we get the total displacement given by

    Fig. 5 Deformation of the rib.

    3.2. MTE control

    The deformation of the MTE is controlled by a servomotor via a rod connected to the trailing edge (see Fig. 1). In order for the system to operate properly, the actuator must deliver suff icient torque.A FEA was performed to obtain the needed torque for the control of the wing under aerodynamical pressures(20 m/s=65.62 ft/s) with an angle of attack of 15°, as shown on Fig.6,the surface stresses caused by the aerodynamic pressure were low.The servomotor must mainly counter the elastic resistance of the ribs. The geometry of the ribs makes their resistance to be very low. However, the servomotor provides enough forces to counteract the aerodynamic forces on these ribs. For the prototype, the same servomotor as the one for a conventional aileron was used (Fig. 7). Therefore, the morphing mechanism was transferred on a conventional wing without changing its weight from the point of view of the controls.

    In the next section, the manufacturing of the MTE system and, the results of the wind tunnel test are presented.

    Fig. 6 FEA of MTE control under aerodynamical pressures.

    Fig. 7 Servomotor specif ications.

    4. Manufacturing and wind tunnel testing

    4.1. Manufacturing of the MTE system

    LASER cutting machine was used to manufacture both test wings, and allowed obtaining their components quickly and at a reasonable cost. To validate the functionality of this MTE, a prototype was manufactured to deform three ribs as shown in Fig. 8. Fig. 9 shows the system including the ribs deformation.The trailing edge is shown on the right hand side of Figs. 8 and 9. The structure of the wing is simple as is equipped with a main spar and three ribs. The thin spars have a main purpose of helping during the manufacturing process but are not required in the wing structure because balsa sheeting is added to cover the wing and the main spar is designed to support all the loads on the wing.

    4.2. Wind tunnel tests

    The numerical results were validated by means of wind tunnel tests on the various morphing wings. The morphing wings have the same dimensions as the reference wing (Fig. 10(a)),so that the impact of the system on the morphing wing was compared to the same wing without any deformation system(the reference wing). The reference wing has the NACA0012 airfoil with a chord of 10 in (254 mm) and a wingspan of 11.5 in (292.1 mm). The 10 in (254 mm) chord relates to the chord of the test wing of the ATR42 research project at the LARCASE which also was of 10 in (254 mm). The wing span is 12 in(304.8 mm),which corresponds to the mid-height of the wind tunnel test chamber,but 0.5 in(12.7 mm)is embedded in the base of the wing that is attached to the aerodynamic scale,which is the reason why the analysis is carried out on a wing with a span of 11.5 in (292.1 mm). The aerodynamic scales is used to measure the loads (forces and moments) on the aerodynamic model installed in the wind tunnel.

    Fig. 8 Structure of the MTE system.

    Fig. 9 Prototype with a MTE system.

    Fig. 10 Reference rigid wing dimensions.

    An objective of the design of the reference wing was to validate f irstly the functioning of the LARCASE aerodynamic scales by comparing the measured values with the numerical values of the loads (forces and moments) acting on the wing.Fig. 10(b) shows the dimensions of the reference wing.

    Fig. 11 LARCASE Price-Pa?¨doussis subsonic wind tunnel.

    Wind tunnel tests were performed on this reference wing as part of the Master's thesis of the main author,32in which the experimental aerodynamic performances with respect to their numerical values computed by XFLR5 software28were compared, and further validated. The tests were carried out in the Price-Pa?¨doussis wind tunnel's rectangular test chamber(2 ft (609.6 mm)×3 ft (914.4 mm), Fig. 11) at speeds of 20 m/s (65.62 ft/s), 25 m/s (82.02 ft/s), 30 m/s (98.43 ft/s) and 35 m/s (114.83 ft/s). The Price-Pa?¨doussis wind tunnel can reach a maximum speed of 37 m/s (121.39 ft/s) for this test chamber.

    The angle of attack of the wing was changed from -10°to 20° by 1° steps. The measurements were carried out with an aerodynamic scales equipped with an Omega 160 force sensor from ATI Industrial Automation. Its technical details expressed in term of maximum loads (forces and moments)are presented in Table 1. This scale was designed by the LARCASE team. Following the manufacturer's specif ications, it was possible to improve the resolution of the sensor,as well as its accuracy by f iltering the measurement given by the aerodynamic scales. This f iltering allowed us to obtain resolution of 0.01 N with a maximum error of 0.1 N due to the rotation mechanism of the aerodynamic scales. The applied digital f iltering enable higher accuracy sampling of the drag force. The theoretical values of the lift, the drag and the pitch coeff icients were calculated with the ‘‘3D Panels” method by using XFLR5 code for a wing with the span of 11.5 in (254 mm) and the chord of 10 in(292.1 mm). The wing in the wind tunnel was in contact with the f loor of the test chamber, and this contact was simulated using the symmetry of the span in the XFLR5 code. The calculated wing therefore had a total span of 23 in (584.2 mm).Figs. 12-14 show the variation of the drag, lift and pitch coeff icients of the wing with the angles of attack. The measured values in the wind tunnel are compared to the calculated values using XFLR5 code.

    At a speed of 20 m/s (65.62 ft/s), the wing has a stall angle of 18°, and at a speed of 25 m/s (82.02 ft/s), the stall occurs at 19°. For a speed of 30 m/s (98.43 ft/s) and 35 m/s (114.83 ft/s), there was no stall occurrence before 20°. The stall angles at 20 m/s (65.62 ft/s) and 25 m/s (82.02 ft/s) were not obtained during the calculations because of the fact that the calculation method does not allow the coeff icients to be computed.28

    The objective of the static test is to validate that the mechanism of the MTE is morphed when a command is sent to the system. There is the need to validate that this experimental deformation of the MTE corresponds to the calculated deformation.For the static tests,a command was sent to the servomotor controlling the MTE, and the vertical displacement at the tip of the trailing edge was measured.

    Fig. 15 shows the dimensioned slits necessary to calculate the total deformation of a rib of the wing. An optimization method was not used to determine the number of slots and their dimensions. Instead, an iterative method using FEA module of CATIA V5 was used, by manually changing the sizes of the slits according to the stresses on the ribs obtained using FEA.The main objective of static tests was to obtain the deformation of the inclined aileron presented in Fig. 15 with the constraints, expressed in term of the torque available by the actuator and the mechanical strength of the wood used for the prototype (plywood33).

    The Table 2 lists the dimensions of the slits from left to right in inches.

    According to Eq. (6) and values of ln, Lnand pngiven in Table 2, the maximum displacement ytis calculated as

    Table 1 Range and resolution for ATI Omega 160 F/T sensor.

    Fig. 12 Drag coeff icients variation with angle of attack for the reference wing.

    Fig. 16 illustrates the measurement of the displacement of the trailing edge during a 7° rotation of the actuator.

    The distance between the axis of the servomotor that controls the deformation and the axis that moves the trailing edge is 4.25 in(107.95 mm)and the distance between the axis of the servomotor and the trailing edge is 5.6 in (142.24 mm)(Fig. 17). Thus, for a displacement of the servomotor axis of 7°, the displacement of the trailing edge axis is 0.52 in(13.21 mm)according to the Eq.(3).For a linear displacement of the trailing edge, this displacement would be 0.69 in(17.526 mm). Since the deformation is given by six articulations (joints) that are not aligned on the actuator axis, the MTE does not move along the same axis as the actuator arm.These calculations do not give a precise value for the displacement; a FEA method makes it possible to obtain a very good precision on the deformation value.

    Wind Tunnel Tests measurements were performed to validate the dynamic operation of the system. The main objective of these measurements was to validate that the deformation was well carried out in the presence of aerodynamic loads.These tests allowed us to observe the inf luence of the mechanism on the drag of the wing compared to the drag of the reference wing (without any mechanism). In this section, the MTE is not moved (displaced) by the controller. Wind tunnel measurements with controlled displacement of the MTE will be presented in the next section. Fig. 18 presents the test wing used to obtain the following results, shown in Fig. 19.

    Two causes were identif ied for the increase in drag on the morphing wing with respect to the drag for the reference wing:

    (1) The discontinuity between the f ixed section of the wing and the moving section of the wing

    Fig. 13 Lift coeff icients variation with angle of attack for the reference wing.

    (2) The presence of the slits on the surface of the wing allowing deformations.

    The f irst cause is common to both the aileron and the MTE system. The second cause is due only to the specif ic design of the MTE system. For this reason, it was chosen to work on the second cause in order to improve the design of the MTE system.

    Four experiments were therefore def ined to improve the airf low around the wing. The f irst experiment reproduced a rigid wing by covering all the slits on the wing with tapes (Fig. 20).For the second experiment, the tape covering the intersection between the morphing part and the f ixed part of the wing was removed (Fig. 21). For the third experiment, the tapes were installed on one side of the slits (leading edge side) with the aim to hide them but without hampering the deformation(Fig. 22). For the fourth and last experiment (Fig. 23), tapes were installed on both sides of the slits that allowed some deformation of the trailing edge (a bump or a hollow were formed on the tape during the deformation).

    The f irst experiment was conducted to conf irm that the drag difference was not due to an imperfection on the wing.The second experiment was performed to identify the inf luences of the discontinuities between the f ixed part and the morphing part of the wing. Experiments 3 and 4 were conducted on improved design of the morphing test wing and had the aim of improving its aerodynamic performances. The drag coeff icient measurements obtained with experiments 1-4 were traced versus the drag coeff icients obtained numerically for the morphing and reference test wings (Fig. 24).

    During Experiment 1, difference was recorded in the stall angles (from 9° to 16° for positive angles and from -16° to-20° for negative angles). Further studies are necessary to identify the cause of this difference. A drag coeff icient curve variation similar to that of the reference drag coeff icients was observed (Fig. 25).

    Experiment 2 indicates that the discontinuity would have a higher inf luence for an angle of attack close to 0° than the inf luence for higher angles of attack (±5°). The comparison between results obtained for Experiments 3 and 4 is not obvious. Thus Experiment 3 was found to be more eff icient than Experiment 4 for positive angles but less eff icient for negative angles (Experiment 3 shows an asymmetry in the results).However, an analysis of the L/D variation with the angle of attack makes it possible to determine that one experiment is more effective than another (Fig. 26). In Fig. 26, the L/D variation with the angle of attack between 4° and 8° for the reference wing is due to the maximum error of 0.1 N on the measurement of the drag force.

    5.Comparison of aerodynamics coeff icients for a MTE versus a conventional aileron

    Studies17revealed that elastic trailing edge could be used for adaptive wing. To validate in this paper that the MTE is capable of replacing an aileron, wind tunnel tests were performed on both, the MTE and also on an aileron.

    Fig. 14 Pitch coeff icients variation with angle of attack for the reference wing.

    Fig. 15 Size and location of the slits in the MTE rib (unit: in).

    Table 2 Size of the slits in the deformable rib.

    Fig. 16 0.75 in displacement of the trailing edge tip (19.05 mm)for a 7° rotation.

    Fig. 17 Controlling distance.

    Fig. 18 Test wing with MTE.

    To perform these wind tunnel tests, a wing was designed and manufactured with its basic dimensions given in Fig. 27,on which we placed an aileron. The width of the aileron is 25%of the wing chord and its span is 9.5 in(241.3 mm).There is 1 in(25.4 mm)on each side of the aileron as shown in Fig.27.

    Fig. 19 Inf luence of the MTE on the drag coeff icient at 15 m/s(49.21 ft/s).

    Fig.20 Experiment 1 with all slits covered by tapes(continuous surface).

    Fig.21 Experiment 2 tapes covering the intersection between the morphing f ixed parts were removed.

    Fig.22 Experiment 3 with tapes installed on one side of the slits(def lector).

    Fig. 23 Experiment 4 with tapes installed on both sides of the slits.

    To be able to retrieve data eff iciently, a control interface was added for the aileron to the Labview interface of the aerodynamic balance. Labview software makes it possible to program commands for actuator controls, and to read instrument outputs through a graphical interface that displays results. This interface controls a servomotor by indicating an angle for its control arm. The force was calculated and displayed for each angle of the actuator. To obtain an angle of the servomotor, the forces reading were averaged, and were further saved in a table.As long as this position was kept constant by the servomotor, its average value in the array was updated; when the angle changed, it was the average value of this new angle that was updated. By reviewing all the possible angles, the graph of the forces was drawn as a function of the angle of the servomotor, and the results were directly visualized. A backup function was added to write the values table of to a text f ile as a means to save the data. Then these values were studied in order to analyze the wings' aerodynamic performance.

    Fig. 24 CD variation with the angles of attack at the speed of 15 m/s (49.21 ft/s).

    Fig. 25 Zoom-in of the drag coeff icients variation with the angles of attack between -8° and 8°.

    Fig. 26 Wing L/D variation with the angle of attack.

    5.1. Comparison of aerodynamics coeff icients

    Fig. 27 2D plan of the wing with an aileron.

    In order to determine if the MTE system can replace an aileron on a wing, the corresponding lift and drag coeff icients of the two systems were compared. Since the angles of the two systems are different, their behaviors are also different; their respective motions cannot be directly compared because of the difference in size between moving surfaces (Fig. 15) and the control arms (Figs. 16 and 17). However, a similar evolution of the lift coeff icient was observed for the two systems.As the goal of this paper is to determine if one system creates more drag than another system, the drag of both systems was compared for the same generated lift. The following graphs(Figs. 28 and 29) show the variation of the lift and drag coeff icients CLand CDwith the angle of attack for the aileron and the MTE compared to the CLand CDof the reference wing.From these f igures,it was observed that the aileron can create a CLgreater than the CLof the MTE,but the CDof the aileron was also larger than the CDof the MTE.Analysis of the behavior of the CLallowed to def ine the roll that it could induce,and thus, to determine if the MTE can indeed replace ailerons.34This observation alone made it possible to def ine if one system was more eff icient than another.

    It was observed that the wing with the MTE created less drag(Fig.30),and therefore it had a higher L/D(Fig.31)than the wing with an aileron.It can be therefore concluded that the MTE was more eff icient than the aileron for angles of attack between -10° and 10°. Therefore, by replacing the conventional aileron with the MTE should result in a reduction in fuel consumption and further in an increase in aircraft autonomy and eff iciency. This observation might also be valid for other aircraft control surfaces such as elevator and rudder.

    5.2. Electrical consumption

    Fig. 28 CL and CD variation with the angle of attack. At 15 m/s(49.21 ft/s) for each system.

    Fig. 29 L/D variation with the angle of attack at 15 m/s(49.21 ft/s) for each system.

    Fig. 30 CD variation with angle of attack zoomed-in around 0°.

    Fig.31 L/D variation with angle of attack zoomed-in around 0°.

    As the MTE requires deformation of the ribs actuated by the servomotor, the MTE requires more electrical energy to operate than an aileron.For a conventional aileron,without external aerodynamic forces,the servomotor does not need to act to keep the aileron in a f ixed position (the wing is placed vertically) but for the MTE, the electrical resistance of the ribs induces an increase of the current required by the servomotor when its deformation is increased. In order to ascertain the magnitude of this difference,the current consumed by the servomotors of the test wings was measured. For the wing equipped with an aileron, the current consumption remains constant at 3.5 mA from -40° to 40°. For the wing with the MTE, the current consumption varies from 3.5 mA at 0° to 70 mA at±7° (with a jump from 16 mA at 6° to 70 mA at 7°). The current jump is due to the mechanical limit of the MTE. It is therefore important to accurately measure the travel limits of the MTE to avoid over-consumption of the system that could lead to a breakage of the servomotor. On the other hand, for the two test wings, an increase of the current consumption has not been noticed when an external aerodynamic force was applied. The dimensions of our test wings and the speeds used were too small for the forces, on the control surfaces, to disturb signif icantly the servomotors.

    6. Conclusions

    This paper describes the f irst phase of a new morphing wing project at the LARCASE.The MTE was f irstly studied in this paper because this section of the wing was easier to deform(morph) than its other sections. The next step will be to adapt existing MTE design to modify also the leading edge of the wing. That leading edge section is more complicated than the trailing edge section to deform because of its smaller width(20%width for the leading edge and 50%width for the trailing edge), and thus, it is not evident that the morphing leading edge design will achieve a reduction in the wing drag.

    However, the results obtained in this paper have shown an improvement in the effectiveness of the MTE on the wing drag,in particular in the sections where the ailerons are located.This comparative study between the MTE and the rigid aileron has verif ied that the MTE had the same behavior as an aileron regarding the increase of lift with the angle of attack.The drag analysis showed us that the MTE generated less drag than the aileron. Finally, analysis of the L/D has shown that the MTE gave a lower drag for the same lift as an aileron, which led to less fuel consumption.Moreover,the MTE has the same components (rib and actuator) as an aileron, and thus does not increase the weight of the wing (both test wings weight is a total of 725 g (1.6 lb). The use of the MTE led to the satisfaction of the following three criteria: to replace ailerons, to reduce fuel consumption and to maintain the mass of the wing.

    Acknowledgments

    The authors would like to thank Mr Aure′lien Mons for his help in the preliminary design of the MTE, Mr Antoine Machetto for his help in the design of the rotation system of the aerodynamic balance and Mr The′o Andre′ for his help in carrying out tests in the wind tunnel for the wing with an aileron. In addition, the authors would like to thank very much to the NSERC for the Canada Research Chair in Aircraft Modeling and Simulation New Technologies Funding.

    日韩欧美三级三区| 国产高潮美女av| 一区二区三区高清视频在线| 午夜a级毛片| 久久久久久国产a免费观看| 一夜夜www| 午夜爱爱视频在线播放| 级片在线观看| 国产精品野战在线观看| 精品国产露脸久久av麻豆 | 亚洲人成网站高清观看| 伦理电影大哥的女人| 欧美日本亚洲视频在线播放| 丰满人妻一区二区三区视频av| 久久久精品欧美日韩精品| 男女视频在线观看网站免费| 午夜精品在线福利| 久久久久性生活片| 国产精华一区二区三区| 一级av片app| 秋霞伦理黄片| 两性午夜刺激爽爽歪歪视频在线观看| 少妇人妻一区二区三区视频| 亚洲国产欧洲综合997久久,| 男插女下体视频免费在线播放| 国产亚洲一区二区精品| 国产大屁股一区二区在线视频| 国产探花在线观看一区二区| 免费av不卡在线播放| 日本免费一区二区三区高清不卡| 精品久久久久久电影网 | 超碰97精品在线观看| 成人午夜精彩视频在线观看| 久久这里有精品视频免费| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 搡女人真爽免费视频火全软件| 九九热线精品视视频播放| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 级片在线观看| 夫妻性生交免费视频一级片| 精品久久久久久电影网 | 麻豆成人av视频| 日本免费在线观看一区| 人人妻人人澡欧美一区二区| 在线观看美女被高潮喷水网站| 欧美色视频一区免费| 老司机福利观看| 一个人看视频在线观看www免费| 1000部很黄的大片| 国产欧美另类精品又又久久亚洲欧美| 国产美女午夜福利| 真实男女啪啪啪动态图| 三级经典国产精品| 免费不卡的大黄色大毛片视频在线观看 | 国产精品一区二区三区四区免费观看| 亚洲精品久久久久久婷婷小说 | 美女cb高潮喷水在线观看| 晚上一个人看的免费电影| 国产精品野战在线观看| 在线观看av片永久免费下载| av免费在线看不卡| 国产精品久久视频播放| 亚洲三级黄色毛片| 寂寞人妻少妇视频99o| 日韩中字成人| 最近手机中文字幕大全| 国产av一区在线观看免费| 欧美+日韩+精品| 久久热精品热| 日韩欧美三级三区| 精品无人区乱码1区二区| 亚洲在线自拍视频| 国产一级毛片七仙女欲春2| av又黄又爽大尺度在线免费看 | 26uuu在线亚洲综合色| 国产精品三级大全| 久久久国产成人精品二区| 国产成人freesex在线| 亚洲精品自拍成人| 九色成人免费人妻av| 水蜜桃什么品种好| 亚洲人成网站高清观看| 日本-黄色视频高清免费观看| 国产欧美另类精品又又久久亚洲欧美| 毛片女人毛片| 精品久久久久久久久久久久久| 色网站视频免费| 色噜噜av男人的天堂激情| 久久人人爽人人爽人人片va| 精品人妻熟女av久视频| 午夜免费男女啪啪视频观看| 国产精品久久久久久久电影| 国产伦精品一区二区三区视频9| 国产av不卡久久| 少妇被粗大猛烈的视频| 国产精品乱码一区二三区的特点| 岛国毛片在线播放| av天堂中文字幕网| 久久久亚洲精品成人影院| 午夜日本视频在线| 日本-黄色视频高清免费观看| 亚洲欧美精品自产自拍| 亚洲国产精品合色在线| 春色校园在线视频观看| 亚洲一级一片aⅴ在线观看| 免费电影在线观看免费观看| 亚洲av中文av极速乱| 日韩av在线免费看完整版不卡| h日本视频在线播放| 亚洲精品,欧美精品| 国产精品国产三级专区第一集| 亚洲精品亚洲一区二区| 成人毛片60女人毛片免费| 日韩,欧美,国产一区二区三区 | 色5月婷婷丁香| 日本一本二区三区精品| 久久精品国产鲁丝片午夜精品| 国产精品一区www在线观看| 久久精品夜色国产| 少妇熟女欧美另类| 久久国内精品自在自线图片| 欧美成人精品欧美一级黄| 直男gayav资源| 国产午夜精品论理片| 亚洲成人av在线免费| 日韩人妻高清精品专区| 日韩精品青青久久久久久| 国产极品天堂在线| 国产国拍精品亚洲av在线观看| 欧美一区二区亚洲| 国产亚洲av嫩草精品影院| 在线免费观看不下载黄p国产| 久热久热在线精品观看| 波多野结衣巨乳人妻| 亚洲欧美日韩东京热| 国内少妇人妻偷人精品xxx网站| 亚洲欧洲国产日韩| 嫩草影院新地址| 国产真实乱freesex| 国产精品麻豆人妻色哟哟久久 | 欧美97在线视频| 国产毛片a区久久久久| 国产又黄又爽又无遮挡在线| 国产午夜精品一二区理论片| 亚洲,欧美,日韩| 国产精品一区二区三区四区免费观看| 九九在线视频观看精品| 国产白丝娇喘喷水9色精品| 天堂中文最新版在线下载 | 神马国产精品三级电影在线观看| 国产高潮美女av| 久久精品91蜜桃| 国产成人精品久久久久久| 日日啪夜夜撸| 久久精品人妻少妇| 哪个播放器可以免费观看大片| 三级经典国产精品| 一区二区三区乱码不卡18| 丰满少妇做爰视频| 黑人高潮一二区| 美女黄网站色视频| 舔av片在线| 最近手机中文字幕大全| 国产69精品久久久久777片| 国产精品野战在线观看| 国产成人精品久久久久久| 国产精品一区二区三区四区久久| 又粗又硬又长又爽又黄的视频| 99久国产av精品| 少妇的逼水好多| 黄色配什么色好看| 小说图片视频综合网站| 亚洲色图av天堂| 免费观看在线日韩| 又粗又爽又猛毛片免费看| 亚洲精品亚洲一区二区| 亚洲乱码一区二区免费版| 国产精品熟女久久久久浪| 久久精品久久久久久久性| 老师上课跳d突然被开到最大视频| 日韩一区二区视频免费看| 亚洲成色77777| 搡女人真爽免费视频火全软件| 一边亲一边摸免费视频| 色综合亚洲欧美另类图片| 中文精品一卡2卡3卡4更新| 日本一本二区三区精品| 中文字幕av在线有码专区| 免费看光身美女| 人人妻人人澡欧美一区二区| 成人特级av手机在线观看| 禁无遮挡网站| 日韩 亚洲 欧美在线| 亚洲av.av天堂| 中国国产av一级| 久久6这里有精品| 精品人妻熟女av久视频| 国产又黄又爽又无遮挡在线| 久久久久免费精品人妻一区二区| 岛国毛片在线播放| 国产av码专区亚洲av| 久久国产乱子免费精品| 亚洲av福利一区| 岛国在线免费视频观看| 精品一区二区三区视频在线| 真实男女啪啪啪动态图| 久久久久国产网址| 国产精品综合久久久久久久免费| 国产不卡一卡二| 99九九线精品视频在线观看视频| 熟女人妻精品中文字幕| 国产亚洲av片在线观看秒播厂 | av免费在线看不卡| 超碰97精品在线观看| 成年版毛片免费区| 国产黄a三级三级三级人| 亚洲av不卡在线观看| 国产精品.久久久| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 夫妻性生交免费视频一级片| 亚洲精品日韩在线中文字幕| 成人三级黄色视频| 亚洲在线自拍视频| 一个人看视频在线观看www免费| 国产老妇女一区| 国内精品美女久久久久久| 熟女电影av网| 1000部很黄的大片| 在线观看一区二区三区| 久久精品影院6| 能在线免费看毛片的网站| 干丝袜人妻中文字幕| 亚洲丝袜综合中文字幕| 久久久久久久久久久免费av| 中文天堂在线官网| 国产男人的电影天堂91| 日韩av在线大香蕉| 成人鲁丝片一二三区免费| 国产免费又黄又爽又色| 欧美人与善性xxx| 天堂av国产一区二区熟女人妻| 国产精品乱码一区二三区的特点| 1000部很黄的大片| 亚洲av成人av| 精品一区二区三区视频在线| 国产在线男女| 我的老师免费观看完整版| av国产免费在线观看| 亚洲综合色惰| 亚洲最大成人av| 亚洲婷婷狠狠爱综合网| 五月伊人婷婷丁香| 两个人的视频大全免费| av播播在线观看一区| 婷婷六月久久综合丁香| 免费看美女性在线毛片视频| 国产人妻一区二区三区在| 亚洲国产精品国产精品| 一级黄色大片毛片| 精品久久久久久电影网 | 午夜精品在线福利| 男女那种视频在线观看| 啦啦啦啦在线视频资源| 日本av手机在线免费观看| 亚洲高清免费不卡视频| 桃色一区二区三区在线观看| 亚洲成av人片在线播放无| 亚洲精品aⅴ在线观看| 亚洲精品乱码久久久久久按摩| 国产私拍福利视频在线观看| 久久精品久久精品一区二区三区| 51国产日韩欧美| 国产成人一区二区在线| 国产乱人视频| 久久人人爽人人爽人人片va| 久久久久九九精品影院| 国产 一区 欧美 日韩| 女人久久www免费人成看片 | 99九九线精品视频在线观看视频| 国产成人aa在线观看| 丝袜喷水一区| 美女脱内裤让男人舔精品视频| 99热网站在线观看| 一本久久精品| 免费看日本二区| 久久精品国产亚洲av涩爱| 天天躁夜夜躁狠狠久久av| 校园人妻丝袜中文字幕| 免费一级毛片在线播放高清视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费无遮挡裸体视频| 免费av毛片视频| 午夜福利成人在线免费观看| 国产午夜福利久久久久久| 免费观看人在逋| 亚洲人成网站在线观看播放| 美女国产视频在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲成人av在线免费| 激情 狠狠 欧美| 久久精品熟女亚洲av麻豆精品 | 尾随美女入室| 变态另类丝袜制服| 亚洲久久久久久中文字幕| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 日韩在线高清观看一区二区三区| 亚洲伊人久久精品综合 | 亚洲av二区三区四区| 久久精品影院6| 国产高清三级在线| www日本黄色视频网| 国内精品宾馆在线| 在线观看66精品国产| 国产高清国产精品国产三级 | 最近最新中文字幕免费大全7| 久久久久性生活片| 国产免费一级a男人的天堂| 国模一区二区三区四区视频| 色5月婷婷丁香| 色哟哟·www| 日本三级黄在线观看| 久久久午夜欧美精品| 亚洲欧美中文字幕日韩二区| 亚洲成色77777| 舔av片在线| 国产午夜精品一二区理论片| 亚洲伊人久久精品综合 | 波野结衣二区三区在线| 精品久久国产蜜桃| 国产亚洲av片在线观看秒播厂 | 日本wwww免费看| 国产麻豆成人av免费视频| 亚洲成色77777| 卡戴珊不雅视频在线播放| 人人妻人人看人人澡| 免费电影在线观看免费观看| 国产在视频线在精品| 永久免费av网站大全| 91久久精品国产一区二区成人| 久久精品国产亚洲av涩爱| 变态另类丝袜制服| 成人欧美大片| 亚洲av成人av| 亚洲精品乱码久久久v下载方式| 国产午夜精品久久久久久一区二区三区| 国产精华一区二区三区| 国产精品久久久久久久电影| 69人妻影院| 国产av一区在线观看免费| 日本wwww免费看| 色视频www国产| 免费看a级黄色片| 国产亚洲av嫩草精品影院| www.av在线官网国产| 91在线精品国自产拍蜜月| 波多野结衣巨乳人妻| 国产精品日韩av在线免费观看| 中文在线观看免费www的网站| 色网站视频免费| 成人无遮挡网站| 国产探花极品一区二区| 日韩成人av中文字幕在线观看| 日韩一区二区视频免费看| 日韩欧美三级三区| 久久精品国产鲁丝片午夜精品| 亚洲成人精品中文字幕电影| 成人二区视频| av在线观看视频网站免费| 午夜激情欧美在线| 一区二区三区高清视频在线| 麻豆成人午夜福利视频| 永久网站在线| 国产黄色小视频在线观看| 欧美+日韩+精品| 久久久久久国产a免费观看| 国产激情偷乱视频一区二区| 亚洲av.av天堂| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 一级黄片播放器| 免费看a级黄色片| 最近2019中文字幕mv第一页| 日本-黄色视频高清免费观看| 亚洲av免费高清在线观看| 久久久久久久久久久免费av| 精品久久久久久久末码| 成人高潮视频无遮挡免费网站| 一区二区三区四区激情视频| 亚洲无线观看免费| 九色成人免费人妻av| 日本一本二区三区精品| 精华霜和精华液先用哪个| 久久久久久久国产电影| 五月伊人婷婷丁香| 99久久成人亚洲精品观看| 免费观看a级毛片全部| 亚洲av不卡在线观看| 亚洲乱码一区二区免费版| 国模一区二区三区四区视频| 久久久久免费精品人妻一区二区| 黄片wwwwww| 超碰av人人做人人爽久久| 女人被狂操c到高潮| 国产av码专区亚洲av| 有码 亚洲区| 国产v大片淫在线免费观看| 天堂中文最新版在线下载 | 成年av动漫网址| 国产成人a区在线观看| 中文字幕av成人在线电影| 能在线免费看毛片的网站| 国产在线男女| 久久婷婷人人爽人人干人人爱| 色视频www国产| 欧美日韩综合久久久久久| 亚洲内射少妇av| 老司机福利观看| 国产极品天堂在线| 一级黄片播放器| 亚洲电影在线观看av| 亚洲国产精品成人久久小说| 日本免费在线观看一区| 精品人妻熟女av久视频| 午夜福利高清视频| 亚洲欧美一区二区三区国产| 国产乱人视频| 能在线免费看毛片的网站| 免费av不卡在线播放| 免费大片18禁| 欧美最新免费一区二区三区| 18禁在线无遮挡免费观看视频| 国产精品久久久久久久电影| 日韩高清综合在线| 我的老师免费观看完整版| 大香蕉97超碰在线| 又粗又爽又猛毛片免费看| 狂野欧美白嫩少妇大欣赏| 成年女人看的毛片在线观看| 大香蕉97超碰在线| 免费观看在线日韩| 中文字幕精品亚洲无线码一区| 三级国产精品欧美在线观看| 国产淫片久久久久久久久| 亚洲精品乱码久久久久久按摩| 欧美3d第一页| 亚洲自拍偷在线| 蜜桃久久精品国产亚洲av| 久久99精品国语久久久| 午夜a级毛片| 亚洲国产成人一精品久久久| 日韩中字成人| 国产综合懂色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av成人精品一二三区| 一卡2卡三卡四卡精品乱码亚洲| 我的女老师完整版在线观看| 国产在视频线在精品| 婷婷六月久久综合丁香| 成人一区二区视频在线观看| 麻豆国产97在线/欧美| 中国美白少妇内射xxxbb| 中国国产av一级| 白带黄色成豆腐渣| 欧美成人a在线观看| 成人高潮视频无遮挡免费网站| 天堂√8在线中文| 国产精品爽爽va在线观看网站| 亚洲综合精品二区| 久久久久久久国产电影| 日本免费一区二区三区高清不卡| 丰满少妇做爰视频| 免费av观看视频| 久久婷婷人人爽人人干人人爱| 亚洲一区高清亚洲精品| 在线天堂最新版资源| 久久久精品欧美日韩精品| 丝袜喷水一区| 国产av在哪里看| 最近2019中文字幕mv第一页| 少妇高潮的动态图| 国产高清国产精品国产三级 | 菩萨蛮人人尽说江南好唐韦庄 | 久久精品熟女亚洲av麻豆精品 | 免费人成在线观看视频色| 国产精品一及| 99久久成人亚洲精品观看| 国产成人免费观看mmmm| 成人午夜精彩视频在线观看| 国产亚洲最大av| 黄色配什么色好看| 亚洲婷婷狠狠爱综合网| 国产老妇女一区| 国产精品美女特级片免费视频播放器| 亚洲精品自拍成人| 欧美成人午夜免费资源| 高清日韩中文字幕在线| 国产又色又爽无遮挡免| 久久综合国产亚洲精品| 中国美白少妇内射xxxbb| 最近中文字幕高清免费大全6| 亚洲性久久影院| 性色avwww在线观看| 国产极品精品免费视频能看的| 成年免费大片在线观看| 亚洲精品成人久久久久久| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 麻豆一二三区av精品| 卡戴珊不雅视频在线播放| .国产精品久久| 国产黄色视频一区二区在线观看 | 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产精品sss在线观看| 如何舔出高潮| 五月玫瑰六月丁香| 午夜视频国产福利| 99久久无色码亚洲精品果冻| 中文亚洲av片在线观看爽| 最近的中文字幕免费完整| 国产精品福利在线免费观看| 国产91av在线免费观看| 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 亚洲国产最新在线播放| 精品一区二区三区人妻视频| 18禁动态无遮挡网站| 久久婷婷人人爽人人干人人爱| 国产在视频线在精品| 中国国产av一级| 国内精品一区二区在线观看| 成人性生交大片免费视频hd| 久久99蜜桃精品久久| av在线老鸭窝| 色综合站精品国产| 老司机影院毛片| 国产精品熟女久久久久浪| 91av网一区二区| 别揉我奶头 嗯啊视频| 国产精品日韩av在线免费观看| 搞女人的毛片| 一二三四中文在线观看免费高清| 中文字幕av成人在线电影| 麻豆精品久久久久久蜜桃| 亚洲av.av天堂| 亚洲丝袜综合中文字幕| av卡一久久| 亚洲电影在线观看av| 亚洲精品456在线播放app| 黄片无遮挡物在线观看| 水蜜桃什么品种好| 九草在线视频观看| 在线天堂最新版资源| av视频在线观看入口| 看片在线看免费视频| 亚洲国产欧美在线一区| 亚洲久久久久久中文字幕| 我的老师免费观看完整版| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 91在线精品国自产拍蜜月| 亚洲精华国产精华液的使用体验| 日韩强制内射视频| 可以在线观看毛片的网站| 免费av毛片视频| 高清毛片免费看| 色吧在线观看| 亚洲图色成人| 久久精品夜夜夜夜夜久久蜜豆| 国产乱人偷精品视频| 欧美bdsm另类| 精品久久久久久成人av| 美女黄网站色视频| av黄色大香蕉| 亚洲av成人精品一二三区| 国产成人精品一,二区| av免费在线看不卡| 桃色一区二区三区在线观看| 22中文网久久字幕| 嫩草影院入口| 精品熟女少妇av免费看| 亚洲综合精品二区| 欧美最新免费一区二区三区| 久久99热6这里只有精品| 欧美日韩精品成人综合77777| 国产在线男女| 一夜夜www| 又粗又爽又猛毛片免费看| 91精品国产九色| 色哟哟·www| 久久久久国产网址| av在线播放精品| 亚洲电影在线观看av| 亚洲av电影不卡..在线观看| 97超视频在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品熟女久久久久浪| 搡女人真爽免费视频火全软件| 国产淫语在线视频| 97热精品久久久久久| 日韩人妻高清精品专区| 国产大屁股一区二区在线视频| 国产一区二区在线av高清观看| 黄色日韩在线| 日韩 亚洲 欧美在线| 色吧在线观看| 蜜臀久久99精品久久宅男| 一区二区三区乱码不卡18| 三级经典国产精品|