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

    An experimental study of elastic properties of dragonfly-like flapping wings for use in biomimetic micro air vehicles(BMAVs)

    2017-11-20 12:07:30PrveenNirSivsnkrnThomsArthurWrdErfnSlmiRubentherenViypuriChristopherFerdyMohdRfieJohn
    CHINESE JOURNAL OF AERONAUTICS 2017年2期

    Prveen Nir Sivsnkrn,Thoms Arthur Wrd,Erfn Slmi,Rubentheren Viypuri,Christopher J.Ferdy,Mohd Rfie John

    aDepartment of Mechanical Engineering,University of Malaya,Kuala Lumpur 50603,Malaysia

    bSchool of Engineering and Physical Sciences,Heriot-Watt University Malaysia,Putrajaya 62200,Malaysia

    cDepartment of Electrical Engineering,University of Malaya,Kuala Lumpur 50603,Malaysia

    An experimental study of elastic properties of dragonfly-like flapping wings for use in biomimetic micro air vehicles(BMAVs)

    Praveena Nair Sivasankarana,Thomas Arthur Wardb,*,Erfan Salamia,Rubentheren Viyapuria,Christopher J.Feardayc,Mohd Rafie Johana

    aDepartment of Mechanical Engineering,University of Malaya,Kuala Lumpur 50603,Malaysia

    bSchool of Engineering and Physical Sciences,Heriot-Watt University Malaysia,Putrajaya 62200,Malaysia

    cDepartment of Electrical Engineering,University of Malaya,Kuala Lumpur 50603,Malaysia

    ABS;Acrylic;Biomimetic micro air vehicle;Flapping mechanism;PLA;Wing structure

    This article studies the elastic properties of several biomimetic micro air vehicle(BMAV)wings that are based on a dragonfly wing.BMAVs are a new class of unmanned micro-sized air vehicles that mimic the flapping wing motion of flying biological organisms(e.g.,insects,birds,and bats).Three structurally identical wings were fabricated using different materials:acrylonitrile butadiene styrene(ABS),polylactic acid(PLA),and acrylic.Simplified wing frame structures were fabricated from these materials and then a nanocomposite film was adhered to them which mimics the membrane of an actual dragonfly.These wings were then attached to an electromagnetic actuator and passively flapped at frequencies of 10–250 Hz.A three-dimensional high frame rate imaging system was used to capture the flapping motions of these wings at a resolution of 320 pixels?240 pixels and 35000 frames per second.The maximum bending angle,maximum wing tip deflection,maximum wing tip twist angle,and wing tip twist speed of each wing were measured and compared to each other and the actual dragonfly wing.The results show that the ABS wing has considerable flexibility in the chordwise direction,whereas the PLA and acrylic wings show better conformity to an actual dragonfly wing in the spanwise direction.Past studies have shown that the aerodynamic performance of a BMAV flapping wing is enhanced if its chordwise flexibility is increased and its spanwise flexibility is reduced.Therefore,the ABS wing(fabricated using a 3D printer)shows the most promising results for future applications.

    1.Introduction

    Micro air vehicles(MAVs)are a relatively new and rapidly growing area of aerospace research.They were first defined by the US Defense Advanced Research Projects Agency(DARPA)in 1997 as unmanned aircraft that are less than 15 cm in any dimension.Later in 2005,the DARPA defined aircraft with all dimensions less than 7.5 cm and lighter than 10 g(carrying 2 g payload)as nano air vehicles(NAVs).MAVs(or NAVs)generally fit into one of the three categories:fixed wing,rotorcraft,or biomimetic.Biomimetic MAVs(BMAVs)mimic the flapping wing motion of flying organisms(e.g.,insects,birds,bats,etc.).This allows lift and thrust to be achieved from a relatively small wing surface area.This allows BMAVs to be potentially smaller and more lightweight than the other two types.These characteristics make BMAVs ideally suited for flight missions in con fined areas(e.g.,around power lines,narrow streets,indoors,etc.).Therefore,BMAV structural components must be ultra-lightweight,compact,and flexible.Most past MAV research has focused on fixed wings,which are essentially scaled-down versions of wings on conventional fixed-wing aircraft.These wings are unsuitable for BMAVs due to their lack of flexibility,so a new type of structural wing design is required for BMAVs.In this work,a dragonfly wing structure is mimicked to construct a new BMAV wing design.A dragonfly(Odonata)was selected for biomimicry,because they are highly maneuverable flyers,capable of hovering,rapid forward flight,and reverse flight.Therefore,structurally analyzing their wings could yield results that bioinspire the design of more effective wings for BMAVs.This article follows on from research discussed in a previous article(written by the authors)that analyzed the static strength of dragonfly-like wing frames fabricated from common materials used in unmanned aircraft(balsa wood,black graphite carbon fiber,and red pre-impregnated fiberglass).1

    Several past research studies have been conducted on flying insect wing structures to understand their elastic properties.Wootton et al.2conducted numerical investigations on a tethered desert locust(Schistocerca gregaria).They concluded that the wings must undergo an appropriate elastic wing deformation(through the course of a wing beat)in order to achieve an efficient aerodynamic flow suitable for lift and thrust generation.Several studies showed that flexible wings,capable of changing their camber,generate higher peak lift forces than those of rigid wings.2,3Wing flexibility also prevents small tears or warping from occurring.Young et al.4suggested that dragonfly wings appear to be adapted for reversible failures in response to excess loads,enabling them to avoid permanent structural damages.Zhu et al.5conducted a study on the effect of flexibility on flapping wing performance during forward flight.A two-dimensional numerical simulation was done by solving the unsteady incompressible Navier-Stokes equations,coupled with the structural dynamic equation for the motion of a wing.The results showed that the flexibility of a flapping wing can largely in fluence its aerodynamic characteristics.The flexibility can simultaneously increase both the propulsive and lifting ef ficiencies of the wing.Kei et al.6conducted a study in which deformation of wings was modeled to examine the effects of bending and torsion on the aerodynamic forces.Their numerical simulations demonstrated that flexible torsion reduces flight instability.They concluded that living butterflies have structurally flexible wings that improve both the aerodynamic efficiency and flight stability.Their experimental measurements showed that a uniformly flexible wing generates lower aerodynamic forces than those of rigid wings under steady-state conditions.However,the presence of wing veins can substantially enhance aerodynamic performance to match or improve the rigid airfoil.These observations agree with those of Zhao et al.7who concluded that flexible insect wings generate greater forces due to an enhanced camber in flight.

    Luo and Tian et al.8,9found that the chordwise deformation of an elastic wing is greater during upstroke than during downstroke.In a study conducted by Ha et al.,10the asymmetric bending of anAllomyrina dichotomabeetle’s hind wing was investigated.Five differently cambered wings were modeled using the ANSYS finite element analysis software.These models were subjected to loads and pressures from the dorsal and ventral sides.The results revealed that both the stressed stiffening of the membrane and the wing camber affect the bending asymmetry of insect wings.In particular,increasing the chordwise camber increased the rigidity of the wing when a load was applied on the ventral side.Alternatively,increasing the spanwise camber increased the rigidity of the wing when a load was applied on the dorsal side.These results explained the bending asymmetry behavior of flapping insect wings.Yang et al.11conducted research on the effects of chordwise and spanwise flexibility on the aerodynamic performance of micro-sized flapping wings.Four flapping motions were described:pure rigid flapping(no deformation),pure spanwise flapping,pure chordwise flapping,and combined chord-spanwise flapping motions.Their results showed that a large spanwise de flection reduces the aerodynamic performance(e.g.,lift and thrust generation)and a large chordwise de flection increases the performance.They further suggested that the design of a flexible flapping wing should incorporate characteristics that will create a suitable chordwise deformation angle(25?and above)and limit the spanwise deformation angle(5?and below).

    Mountcastle and Combes12conducted an experiment using artificially stiffened bumblebee wings(in vivo)by applying a micro-splint to a single flexible vein joint.The bees were then subjected to load-lifting tests.Bees with stiffened wings showed an 8.6 percent maximum lift reduction.This reduction cannot be accounted for by changes in gross wing kinematics,since the stroke amplitude and flapping frequency were unchanged.The results revealed that flexible wing design and the resulting passive wing deformations enhance the load-lifting capacity in bumblebees.Wu et al.13presented a multidisciplinary experiment that correlated a flapping wing’s elasticity and thrust production,by quantifying and comparing overallthrust,structural deformation, and airflow. Six pairs of hummingbird-shaped membrane wings of different properties were examined.The results showed that,for a specific spatial distribution of flexibility,there is an effective frequency range in thrustproduction.Thewing deformation at thrustproducing wing beat frequencies indicated the importance of flexibility.Both bending and twisting motions interact with aerodynamic loads to enhance wing performance.

    Most past research,that is similar to the objectives of this article,examined the effects of wing flexibility on aerodynamic performance by either using numerical models or experimentation.However,very few researchers have attempted to mimic the detailed structure of an actual insect wing.In this article,biomimicry of a dragonfly wing(frame structure and membrane)is done by fabricating them with different materials:acrylonitrile butadiene styrene(ABS),polylactic acid(PLA),and acrylic.The focus of this article is solely on the flexibility of the fabricated wing structures,not the resulting aerodynamic forces that are generated.The wings were fixed to a flapping mechanism and flapped at variable wing beat frequencies.An actual dragonfly has a natural frequency of 120–170 Hz and a wing beat frequency of 30 Hz.The mechanism used in this study was able to flap up to a maximum wing beat frequency of 250 Hz.This allowed us to study the deformation of wing motions at frequencies beyond the ability of an actual dragonfly.The resulting wing tip deflection,twisting angles,twisting speed,and bending angles were measured using imagery generated by two high frame rate cameras.Comparisons were made with a real dragonfly wing in passive flapping motion.

    2.Materials and methodology

    2.1.Wing design and fabrication

    Fig.1 shows the comparison between an actual dragonfly wing(Diplacodes Bipunctata)and the simplified wing frame structure used in this study.The simplified frame structure was designed based on spatial network analysis,which has been described in a previous article written by the authors.14This analysis utilizes geometric objects within a region specified by vertices or edges.Although this method is commonly used in geographical information systems(GIS)to explore geographic spatial patterns,the idea of applying this algorithm to a biological structure was first introduced in this article.It was inspired by observing the compactly arranged geometrical patterns inherent to dragonfly wings.This method allows this complex biological structure to be mimicked by a simplified frame structure that can be fabricated by machining or 3D printing.

    All of the simplified frame structures were fabricated to be approximately 55 mm in length and 0.05 mm thick.As previously mentioned,they were constructed of three different materials:ABS,PLA,and acrylic(see Fig.2 and Table 1).The ABS and PLA wings were fabricated using a Maker Bot Replicator 2X 3D printer.The acrylic wings were fabricated using micro laser machining.Acrylic or PLA is generally known for its resistance to breakage,elasticity,and flexibility.15,16ABS and PLA are the two most dominant plastics used for 3D printing.ABS is chosen due to its strength,flexibility,and machinability,10while PLA is chosen for its biodegradability,lightweight,flexibility,and elasticity.17The densities of ABS,PLA,and acrylic are 1.05,1.19,and 1.18 g/cm3,respectively.A finite element analysis on von Mises stress was conducted to simulate the flexibility of the materials tested.

    A chitosan nanocomposite film was bonded to the wing frames to serve as a thin(3 mm),ultra-lightweight wing membrane.This chitosan nanocomposite film was developed by our research team for this specific purpose and has been the subject of another article.18It has similar properties to those of the chitin membranes of real dragonflies.It is formed by reinforcing a chitosan suspension with nanometer-scaled nanocrystalline cellulose(NCC)particles and tannic acid.This allows both the mechanical properties and water resistivity of the chitosan film to be controlled to achieve suitable design values.The use of NCC as a filler material elevates the film’s mechanical properties(e.g.,rigidity).The addition of tannic acid as a cross-linking agent reduces the swelling behavior,solubility,and rigidity of the nanocomposite film.The film was adhered to a wing frame by firstly submerging the frame into the nanocomposite solution.This procedure also ensured that the film membrane would have a prescribed,uniform thickness and that both sides of the frame structure were evenly coated.The suspension was then transformed into a film by the casting evaporation method.Once cured,the film created a shiny,transparent film layer that adhered firmly to the frame structure.

    Fig.1 Dragonfly wing structure comparison.

    Fig.2 Wing frame materials of PLA,acrylic,and ABS,respectively.

    2.2.Wing flapping mechanism

    The wing flapping mechanism used in this study was an electromagnetic flapping wing actuator.The power supply usedin this flapping wing drive was 9 V DC.An LM555 crystal clock oscillator integrated circuit(shown in Fig.3)was used to generate a stable oscillation.The free running frequency and duty cycle were accurately controlled with two resistors and one capacitor.The generated oscillation was fed to a Power MOSFET fast switch.The output of the Power MOSFET was used to actuate the miniature PC Board Relay.The frequency of the switch(corresponding to the wing beat frequency)can be adjusted by a 22 kX potentiometer.Each of the different wings was attached to a flat iron plate(2 mm long and 2.75 mm thick)using super glue.This plate(wing platform)was oscillated by an electromagnetic actuator(3 mm?3 mm).Fig.3 shows a wing structure attached to the actuator.The plate was attached to the hinge of the wing to mimic the joint of an actual dragonfly.This flapping mechanism is able to create a linear up-down stroke motion at variable wing beat frequencies,up to a maximum frequency of 250 Hz.The flapping degree was set to be 60?which corresponds to an actual dragonfly wing flapping angle during hovering flight.15,19

    Table 1 Mechanical properties of frame structure materials.15,17

    Fig.3 Flapping mechanism used in this study.

    2.3.Experimental set-up

    Two Phantom Miro310(Vision Research)high frame rate cameras were used to view a flapping wing from two different directions.The cameras’high frame rate enables a precise sequence of images to be captured of the flapping wing motion within a single wing beat.Two cameras were necessary in order to determine the three-dimensional shape and orientation of the wing surface(see Fig.4).The cameras were placed perpendicular to each other following the procedures established by Gui et al.20Both cameras were equipped with a Nikon F lens.A multiple LED lighting system was used to provide sufficient illumination.Imagery wasrecorded ata resolution of 320 pixel?240 pixel and a frame rate of 35000 per second,which allowed the wing beat motion to be precisely captured.The motion videos were stored to a computer via two highspeed Ethernet cables.They were played-back and analyzed using the Vision Research Phantom Camera Control Software(Version 2.6.749.0).

    Measurements were taken of each of the three wings while flapping at varying frequencies:10–250 Hz.Fig.5 shows the front and side views of the wing motions that were measured and recorded from captured imagery.Fig.5(a)illustrates the bending angle(h)and the displaced distance or deflection(d).Fig.5(b)defines the wing tip angle(a)and the wing tip rotational twist speed(x).

    3.Results and discussion

    3.1.Stress simulation results(without a membrane)

    A stress simulation analysis was done on the wing frame materials(without and with a membrane)tested in this experiment using Autodesk Simulation Multiphysics 2015.These results directly relate to the flexibilities of the materials tested in this experiment.The results are shown in Figs.6 and 7.

    Fig.6 shows the von Mises stress results of all the three different frame structures.The highest stresses in the forewing recorded for PLA,acrylic,and ABS are 13,17,and 23 N/mm2,respectively.This shows that ABS is the least flexible material among all three materials tested without a membrane.

    3.2.Stress simulation results(with a membrane)

    Fig.7 shows the forewing models of all three materials used in this experiment.Based on Fig.6,the maximum von Mises stress occurs at approximately the same location for all three materials.The highest stresses occur in regions where the surface-to-area ratio is minimum.The maximum stresses recorded are 14.77,17.29,and 24.23 N/mm2for PLA,acrylic,and ABS,respectively.Both Figs.6 and 7 show that ABS exhibits the maximum stress among all three materials.

    Fig.4 Experimental set-up:two high-speed cameras perpendicular to each other along with multi LED lighting.

    Fig.5 Front and side views of wing motion captured(and measurement axes).

    Fig.6 Stress simulation results for ABS,PLA and acrylic(without a membrane).

    Fig.7 Stress simulation results for ABS,PLA,and acrylic(with a membrane).

    3.3.Dragonfly wing flapping motion

    The experiment was conducted on each of the three types of wings(both with and without the chitosan membrane).This was done to study the flexibility of each wing frame material and to determine the best material for use in a BMAV.An actual dragonfly wing(Diplacodes Bipunctata)was also tested to study its motion during passive flapping at different frequencies and compare it with the fabricated wings.The nomenclature for wing rotation about different axes is shown in Fig.8.21Figs.9 and 10 show a sequence of images,illustrating the wing motion of an actual flapping dragonfly wing during one complete flapping cycle.The wing beat frequency for these images was 30 Hz,which is the nominal wing beat frequency of this species of dragonfly.

    Dragonfly wings greatly deform during flight.This was observed in our experiment as well as by others.22Despite having a certain degree of rigidity,dragonfly wings undergo a considerable amount of bending,twisting,and rotational motions.Figs.9 and 10 show the motion of a flapping wing in one complete cycle at 30 Hz(side and front views).It was shown that in both directions(chord and spanwise),an asymmetric twist-bend motion was observed.Figs.9(d),(f),and 10(d)clearly show these asymmetric motions mentioned.At the end of an upstroke(observed in Fig.9(f)),the wing momentarily exhibited a symmetrical twisting motion.A large feathering rotation range of 154?–179?of the entire wing was observed at the beginning of the downstroke and the end of the upstroke(for all frequencies)(Figs.10(a)and(f)).Even during the steady phase(passive moment occurring when the flapping angle was zero),the wing was observed to undergo internal torsion.This corresponds well to previous findings made by Wootton et al.2,23

    Besides the nominal 30 Hz wing beat frequency,the dragonfly wing was also flapped at frequencies ranging from 10 to 250 Hz.The pattern of deformations was similar for all of the frequencies observed.The measured bending angle,wing tip deflection,wing tip twist angle,and speed for the different wing frames(without and with a membrane)were plotted in comparison to the results obtained from an actual dragonfly wing in Figs.11–14.

    Fig.8 Degrees of freedom for wings of a flying insect.21

    3.4.Bending angle versus flapping frequency

    Fig.9 Side view of dragonfly flapping wing(gray scale)captured by high speed camera during one flapping cycle at 30 Hz.

    The bending angle is directly proportional to the flexibility of a wing.Both inertial and aerodynamic loads influence it.Wootton23found that most insect wings have relatively stiff supporting zones near the wing base and leading edge.Adding to this in a later article,Wootton24wrote that the wing veins taper in diameter from base to tip.The resulting reduction in stiffness reduces the inertial load at the wing tip,reducing the energy expenditure and stress at the wing base.Ennos and Wootton25showed that wings having a tapered stiffness distribution from base(high)to tip(low)are well suited to with stand torques.This article also showed that spanwise bending moments due to the inertia of flapping wings are approximately two times larger than those due to aerodynamic forces.A structural finite element analysis by Jongerius and Lentink26of a dragonfly wing model also showed that the inertial forces along the wingspan are 1.5–3 times higher than the aerodynamic forces.Similarly,Combes and Daniel27modeled dragon fly and hawkmoth wings,and found that the flexural stiffnessdeclined exponentially from wing baseto tip.Although inertial loading dominates,Young et al.4showed that aerodynamic forces(e.g.,lift and thrust)generated by a flapping wing also has an in fluence on wing flexibility.

    This study focuses only on the chordwise flexibility of a passive flapping wing.Bending angles were measured along the chordwise direction.Kang and Shyy22also investigated chordwise flexibility,but for simple,non-anisotropic wing structures.They presented a detailed assessment of the effects of structural flexibility on the aerodynamic performance of flap-ping wings.The Reynolds number(Re=100)considered in this study is relevant to small insect flyers,such as fruit flies.However,this study only includes the roles of chordwise flexibility and passive pitch in two-dimensional plunging motions.

    Fig.10 Front view of dragonfly flapping wing captured by high speed camera(gray scale)during one flapping cycle at 30 Hz.

    Fig.11 Bending angles of different wing frames.

    Fig.12 Wing tip deflection of different wing frames.

    Fig.13 Wing twist angle of different frames versus flapping frequency.

    Our study involves a much more complex wing design than those in many past studies.However,tapering the thickness(declination from base to tip)of the veins in our physical models(similar to actual insect wings)was not possible due to fabrication limitations.Ourwingshave tapered flexibility(declining from base to tip and from leading to trailing edge)solely due to a reduction in the frame planform width sizes(mimicking veins)in these directions.Fig.11 shows the bend-ing angles as the wing beat frequency is varied for the three fabricated wing frames(without and with a membrane)in comparison to that of an actual dragonfly wing.Fig.11 shows that the maximum bending angle(hmax)for all the wings occurs during the upstroke.This was observed for both frames without and with a membrane.This agrees with the results of previous research done by Jongerius and Lentink,26in which this asymmetry(difference in the bending angle between the upstroke and downstroke)was attributed to the directional bending stiffness in a wing structure(e.g.,one-way hinge or a pre-existing camber in the wing surface).

    The maximum bending angle of dragonfly wings at 30 Hz was recorded to be about 6?.The wings were observed to have a maximum bending angle of 10.7?at 120 Hz(natural frequency of an actual dragonfly).This is an increase of 78.3%from 30 Hz.ABS shows a high level of flexibility compared to the other two materials used.Fig.11 shows that the bending angle curves of the fabricated ABS wings are more similar to that of the actual dragonfly wing than those of the other two types.Fig.11(a)shows that the bending angle of an ABS wing(without a membrane)at 30 Hz is 8.5?and 5.9?at 120 Hz.At 30 Hz,the percentage difference between an ABS wing(without a membrane)and an actual dragonfly wing is about 42%.The PLA and acrylic wings recorded reduced percentage differences of 30%and 70%,respectively.In Fig.11(b),ABS exhibited much larger bending angles at 30 Hz when the membrane was added.The value of the ABS wing(with a membrane)is 20.1?at 30 Hz and 34.9?at 120 Hz.This angle is much larger than that of the actual dragonfly wing.The percentage increase between the ABS and the actual dragonfly wing is 233%.The other two materials(PLA and acrylic)exhibited much smaller bending angles than that of the actual dragonfly wing.The percentage reductions in PLA and acrylic(in comparison to the actual dragonfly wing)are 83%and 75%,respectively.

    These observations confirm that the overall flexibility of a wing decreases after a membrane is attached,except for ABS wings.At a frequency of 120–170 Hz,the dragonfly wing bends at a very large angle.Previous research showed that dragonflies do not flap at their natural frequency(120–170 Hz).28Therefore,this result is likely due to a resonance effect caused by the wing beat frequency being proximate to the natural frequency of the wing.This result confirms that dragonflies have a maximum wing beat frequency limitation in this range.The ABS wing frame shows a similar trend at 120 Hz.The bending angle is reduced at frequencies greater than120 Hz for both the actual dragonfly wing and the three fabricated wings.

    Fig.14 Wing tip twist speed of different frames versus flapping frequency.

    3.5.Wing tip deflection versus flapping frequency

    Fig.12 shows the wing tip deflection for varying wing beat frequencies of the three fabricated wing frames(without and with a membrane)in comparison to that of an actual dragonfly wing.Similar to the bending angle,deflection is another measurement that can be used to assess a flapping wing’s flexibility.As mentioned earlier,past studies have shown that wing flexibility has a significant effect on the wing’s ability to generate a suitable time-averaged lift or thrust.7Fig.12 shows that the maximum deflection(dmax)occurs during the upstroke.This again was observed for both frames without and with a membrane.This agrees with the results of previous research done by Luo et al.8

    Fig.12(a)shows that all of the fabricated wing frames(without a membrane)deflect at magnitudes that are similar(only slightly reduced)to that of the actual dragonfly wing at 30 Hz which is about 7.1 mm.At 30 Hz,ABS has a percentage increase of 24%.PLA and acrylic both have percentage reductions of 48%and 62%,respectively.However,Fig.12(b)shows that the fabricated wing frames(with a membrane)have very different de flections from that of the actual dragonfly wing.Only the ABS wing showed a comparable level of de flection,however the dragon fly wing is 41%higher than the ABS wing.The PLA and acrylic wings have percentage reductions of 94%and 66%,respectively,compared to the dragonfly wing.The actual dragonfly wing is able to undergo a large deflection at the tip region.This supports the findings of previous studies,which explain that the difference between the deflections at the tip and the surface is created by the difference in the rigidity(due to the vein and corrugations)along the wing surface.29

    The difference in deflection between the wing frames without and with a membrane shows that the attachment of a membrane causes an increase in rigidity.This increase in rigid-ity was observed to be the highest in the PLA wing.Only the ABS wing shows a curvature trend similar to that of the actual dragonfly wing around 120 Hz.At 120 Hz,a reduction in percentage of 45%(without a membrane)or 70%(with a membrane)is seen in the ABS wing frame.Compared to the PLA wing,there is a percentage reduction of 83%(without a membrane)or 95%(with a membrane).The acrylic wing has a percentage reduction of 85%,both without and with the membrane attached.The trend of the graph again shows that there is a decrease in flexibility after the membrane has been attached.Two high peaks were observed for the actual dragonfly wing(30 and 120 Hz).As already stated,the natural frequency of dragon fly wings has been reported to be between 120 to 170 Hz.28The extreme fluctuation observed in this range confirms the reporting.

    3.6.Wing twist angle versus flapping frequency

    Fig.13 shows the maximum wing tip twist angle of the three fabricated wing frames in comparison to that of an actual dragonfly wing.The maximum twist angle was recorded during the stroke reversal(transition from upstroke to downstroke).The twist angle for the actual dragonfly wing at 30 Hz is 154.58?.Untwisted wings have large,drag producing wing surfaces that are exposed to flow,and hence the importance of twisting in wings is justified.Wing tip twist also plays an important role in enhancement of flight performance.The mid-stroke timing of wing deformation in a butterfly,examined by Zheng et al.,29suggests that the deformation is not due to wing inertia,because the acceleration of the wing is small at this point in the stroke.They suggested that this is instead due to elastic effects,since the aerodynamic forces are very large at mid-stroke.

    Figs.13(a)and(b)shows that both the PLA and acrylic wing frames(both without and with a membrane)closely match the performance of an actual dragonfly wing.At 30 Hz,the ABS wing(without and with a membrane)has percentage reductions of 20%and 1%respectively in comparison to that of the actual dragonfly wing.The PLA wing(without and with a membrane)has percentage increases of 5%and 10%,respectively.The acrylic wing(without and with a membrane)has percentage increases of 7%and 12%,respectively.At 120 Hz,the ABS and acrylic wings(without a membrane)have percentage reductions of 10%and 3%,respectively,compared to that of the dragonfly wing,while the PLA wing(without a membrane)has a percentage increase of 3%.The ABS wing(with a membrane)has a percentage reduction of 36%compared to that of the dragonfly wing,while the PLA and acrylic wings have percentage increases of 5%and 3%,respectively.Based on these results,the PLA and acrylic wings are more similar to the actual dragonfly wing than the ABS wing.The large fluctuation of the ABS wing across varying flapping frequencies(10–250 Hz)makes it a more complicated BMAV option.

    Another trend observed from Fig.13 is that the wing tip twist angle of the dragonfly wing does not vary significantly as the flapping frequency is varied.This matches the finding of a previous study by Zhao et al.8(mentioned earlier)which showed that the flexibility of insect wings increases more chordwise than spanwise,due to the rigid leading edge vein.

    This is true for both categories of wing frames(with and without a membrane).

    3.7.Wing tip twist speed versus flapping frequency

    Fig.14 shows the wing tip twist speed for the three wing frames(without and with a membrane)in comparison to that of an actual dragonfly wing.The wing tip twist speed was measured using the Vision Research Phantom Camera Control Software associated with our high frame rate cameras.Vogel30stated that the wing tip twist speed varies according to size and must exceed a ratio to flight speed(wing tip twist speed:flight speed)by 3.7 or more to enable forward flight.Fig.14 shows that the PLA and acrylic wing frames(both without and with a membrane)show a curvature trend similar to that of the actual dragonfly wing.The wing tip twist speed of the actual dragonfly wing at 30 Hz is 9.2 revolutions per second(r/s).At 30 Hz,the PLA wing shows a percentage increase of 32%(without a membrane)and a percentage reduction of 51%(with a membrane)in comparison with that of the dragonfly wing.The acrylic wing shows a percentage reduction of 4%(without a membrane)and an increase of 48%(with a membrane).The ABS wing shows percentage reductions of 16%(without a membrane)and 13%(with a membrane).At 120 Hz,all of the fabricated wing frames without a membrane attached,showed percentage increases in comparison to an actual dragonfly wing.The ABS,PLA,and acrylic wings show percentage increases of 5%,64%,and 13%,respectively,while the ABS,PLA,and acrylic wing frames without a membrane have percentage increases of 32%,88%,and 36%,respectively.The ABS wing frame shows a much different curvature trend from those of the other two,both with and without a membrane.Fig.14 shows that the wing tip twist speed is highly dependent on the flapping frequency and is less influenced by changes in the frame’s flexibility.This can be confirmed by observing the curves of the wing frames with a membrane.The observed trend is the same across varying flapping frequencies(10–250 Hz)for both types of wing frames.

    Combes and Daniel31conducted a finite element analysis study on the wing structures of several different insects(including dragon flies).In all of the species that they tested,spanwise flexure stiffness was one to two orders of magnitude higher than chordwise flexure stiffness.They concluded that stiff leading edge veins played a primary role in generating this anisotropy.Moreover,as previously mentioned,the study conducted by Yang et al.11concluded that spanwise flexible deformation should be limited to a small range(by use of stiff leading edges)in order to achieve a higher aerodynamic performance for a flapping MAV.Alternatively,a larger chordwise deformation could serve to enhance the aerodynamic performance(e.g.,lift and thrust generation).

    The results of our experiments in flapping an actual dragonfly wing support this observation,by showing that the chordwise deformation is very signi ficant compared to the spanwise deformation.These results suggest that BMAV wings should be designed with a stiff leading edge to limit the spanwise deformation and flexible ribs to keep the chordwise deformation within a significant but suitable range.This indicates that the ABS wing design is better suited for use in a BMAV than the PLA and acrylic wing designs.

    4.Conclusions

    One challenge in constructing a working BMAV involves the need to fabricate a highly deformable and flexible wing that has a large load-bearing capacity.An experimental study was conducted to assess elastic properties of flapping wings fabricated from three different materials(ABS,PLA,and acrylic).The structural design of each of these wings is identical and based on biomimicry of an actual dragonfly wing.The experimental results were compared to those of the actual dragonfly wing,on which they were based,in order to assess their potential applications in a BMAV design.A flapping mechanism that utilizes an electromagnetic actuator was used to flap the wings at various frequencies from 10 to 250 Hz.A high frame rate imaging system that includes two cameras was used to capture the three-dimensional motion of a flapping wing.Several different elastic parameters were measured:bending angle,wing tip deflection,wing tip twist angle,and wing tip twist speed.Analysis of the wing bending angle and wing tip deflection indicates the flexibility of the wing in the chordwise direction,while the wing tip twist angle and speed show the flexibility of the wing in the spanwise direction.The ABS wings exhibited the highest chordwise flexibility(indicated by their large bending angles and wing tip de flections).Although the PLA and acrylic fabricated wings exhibited much lower chordwise flexibility than that of the ABS fabricated wings and the dragonfly wing,their spanwise flexibility(indicated by their wing tip twist angles and speeds)closely matched that of the dragonfly wing.

    These experimental results show that an actual dragonfly wing has a highly deformable structure despite its rigidity.The materials examined in this study(ABS,PLA,and acrylic)were selected due to their high flexibility,low density,and low fabrication costs.This study shows that each of these materials is able to perform like an actual dragonfly wing to varying degrees.However,the ABS wing design gave better results in matching the chordwise flexibility of the actual dragonfly wing,while limiting the spanwise flexibility to a much greater degree than those of the other two designs.

    Acknowledgements

    This research was done under the auspices of the Centre for Transportation Research at the Faculty of Engineering,University of Malaya.It was primarily funded by the High Impact Research Grant from the Malaysian Ministry of Higher Education(UM.C/625/1/HIR/MOHE/ENG/12,H-16001-D000012)and a secondarily by a University of Malaya Research Grant(RG155-12AET).

    1.Sivasankaran PN,Ward TA,Viyapuri R,Johan MR.Static strength analysis of dragonfly inspired wings for biomimetic micro air vehicles.Chin J Aeronaut2016;29(2):411–23.

    2.Wootton RJ,Evans KE,Herbert R,Smith CW.The hind wing of the desert locust(Schistocerca gregaria Forskal).I.Functional morphology and mode of operation.J Exp Biol2000;203(19):2921–31.

    3.Mountcastle AM,Daniel TL.Aerodynamic and functional consequences of wing compliance.Exp Fluids2009;46(5):873–82.

    4.Young J,Walker SM,Bomphrey RJ,Taylor GK,Thomas AL.Detailsofinsectwingdesignanddeformationenhanceaerodynamic function and flight ef ficiency.Science2009;325(5947):1549–52.

    5.Zhu JY,Zhou CY,Wang C,Jiang L.Effect of flexibility on flapping wing characteristics under forward flight.Fluid Dyn Res2014;46(5):055515.

    6.Kei S,Takuya O,Masahiko K,Naoto Y,Norio H,Makoto I.Effects of structural flexibility of wings in flapping flight of butter fly.Bioinspiration Biomimetics2012;7(2):025002.

    7.Zhao L,Huang QF,Deng XY,Sane SP.Aerodynamic effects of flexibility in flapping wings.J R Soc Interface2010;7(44):485–97.

    8.Luo HX,Tian FB,Song JL,Lu XY.Aerodynamic cause of the asymmetric wing deformation of insect wings.In:DFD12 meeting of the American Physical Society(APS);2012 Nov 18–20;San Diego(CA).College Park,MD:APS Division of Fluid Dynamics;2012.

    9.Tian FB,Luo HX,Song JL,Lu XY.Force production and asymmetric deformation of a flexible flapping wing in forward lf ight.J Fluids Struct2013;36:149–61.

    10.Ha NS,Truong QT,Goo NS,Park HC.Biomechanical properties of insect wings:The stress stiffening effects on the asymmetric bending of the Allomyrina dichotoma beetle’s hind wing.PLoS One2013;8(12):e80689.

    11.Yang WQ,Song BF,Song WP,Wang LG.The effects of spanwise and chord-wise flexibility on the aerodynamic performance of micro flapping-wing.Chin Sci Bull2012;57(22):2887–97.

    12.Mountcastle AM,Combes SA.Wing flexibility enhances loadlifting capacity in bumblebees.Proc R Soc London B:Biol Sci2013;280(1759):1–8.

    13.Wu P,Stanford BK,Sa¨llstro¨m E,Ukeiley L,Ifju PG.Structural dynamics and aerodynamics measurements of biologically inspired flexible flapping wings.BioinspirationBiomimetics2011;6(1):016009.

    14.Sivasankaran PN,Ward TA.Spatial network analysis to construct simplified wing structural models for biomimetic micro air vehicles.Aerospace Sci Technol2016;49:259–68.

    15.Srigrarom S,Chan WL.Ornithopter type flapping wings for autonomous micro air vehicles.Aerospace2015;2(2):235–78.

    16.Sudo S,Tsuyuki K,Yano T,Takagi K.A magnetic fluid microdevice using insect wings.J Phys:Condensed Matter2008;20(20):204142.

    17.Chilson L.The difference between ABS and PLA for 3D printing.ProtoParadigm[Internet].[updated 2013 Jan 26;cited 2016 Oct].Availablefrom:http://www.protoparadigm.com/news-updates/the-difference-between-abs-and-pla-for-3d-printing/.

    18.Rubentheren V,Ward TA,Chee YC,Praveena N.Physical and chemical reinforcement of chitosan film using nanocrystalline cellulose and tannic acid.Cellulose2015;22(4):2529–41.

    19.Shyy W,Aono H,Chimakurthi SK,Trizila P,Kang CK,Cesnik CES,et al.Recent progress in flapping wing aerodynamics and aeroelasticity.Prog Aerospace Sci2010;46(7):284–327.

    20.Gui L,Fink T,Cao Z,Sun D,Seiner JM,Streett DA.Fire ant alate wing motion data and numerical reconstruction.J Insect Sci2010;10(19):1–17.

    21.Ward TA,Rezadad M,Fearday CJ,Rubentheren V.A review of biomimetic air vehicle research:1984–2014.Int J Micro Air Vehicles2015;7(3):375–94.

    22.Kang CK,Shyy W.Scaling law and enhancement of lift generation of an insect-size hovering flexible wing.J Royal Soc Interface2013;10(85):20130361.

    23.Wootton RJ.Support and deformability in insect wings.J Zool1981;193(4):447–68.

    24.Wootton RJ.Functional morphology of insect wings.Ann Rev Entomol1992;37:113–40.

    25.Ennos AR,Wootton RJ.Functional wing morphology and aerodynamics of Panorpa germanica(Insecta:Mecoptera).J Exp Biol1989;143(1):267–84.

    26.Jongerius SR,Lentink D.Structural analysis of a dragonfly wing.Exp Mech2010;50(9):1323–34.

    27.Combes SA,Daniel TL.Flexural stiffness in insect wings II.Spatial distribution and dynamic wing bending.J Exp Biol2003;206(17):2989–97.

    28.Mingallon M,Ramaswamy S.The architecture of the dragonfly wing:A study of the structural and fluid dynamic capabilities of the anisoptera’s forewing.ASME 2011 international mechanical engineering congress and exposition;2011 Nov 11-17;Denver(CO).New York:ASME;2011.

    29.Zheng LX,Hedrick TL,Mittal R.Time-varying wing twist improves aerodynamic efficiency of forward flight in butterflies.PLoS One2013;8(1):e53060.

    30.Vogel S.Comparative biomechanics:Life’s physical world.2nd ed.Princeton(NJ):Princeton University Press;2013.p.251–63.;Vogel S.Comparative biomechanics:life’s physical world 2nd ed..Princeton(NJ):Princeton University Press;2013.p.251–63.

    31.Combes SA,Daniel TL.Flexural stiffness in insect wings I.Scaling and the influence of wing venation.J Exp Biol2003;206(17):2979–87.

    29 January 2016;revised 25 April 2016;accepted 27 October 2016

    Available online 15 February 2017

    *Corresponding author.

    E-mail address:T.Ward@hw.ac.uk(T.A.Ward).

    Peer review under responsibility of Editorial Committee of CJA.

    如日韩欧美国产精品一区二区三区| 好看av亚洲va欧美ⅴa在| 一二三四在线观看免费中文在| 嫩草影视91久久| 国产精品二区激情视频| 男人操女人黄网站| 成年女人毛片免费观看观看9| av国产精品久久久久影院| 9热在线视频观看99| 国产三级在线视频| 国产一区二区三区在线臀色熟女 | 午夜免费成人在线视频| 亚洲av电影在线进入| 久9热在线精品视频| 久久热在线av| 一区二区三区激情视频| 国产精品一区二区三区四区久久 | 一区福利在线观看| 亚洲国产精品一区二区三区在线| 精品久久蜜臀av无| 男男h啪啪无遮挡| 久久婷婷成人综合色麻豆| 亚洲av美国av| 91国产中文字幕| 正在播放国产对白刺激| 欧美精品亚洲一区二区| 丰满的人妻完整版| 一区二区三区精品91| 在线播放国产精品三级| 一区二区三区激情视频| 国产午夜精品久久久久久| 男人操女人黄网站| 看黄色毛片网站| 老熟妇仑乱视频hdxx| 亚洲熟妇中文字幕五十中出 | 亚洲欧美激情在线| 欧美午夜高清在线| 成年人黄色毛片网站| 久久热在线av| 成人av一区二区三区在线看| 女生性感内裤真人,穿戴方法视频| 男人操女人黄网站| 一级毛片精品| 新久久久久国产一级毛片| 国产一区二区三区在线臀色熟女 | 80岁老熟妇乱子伦牲交| 午夜免费激情av| 在线播放国产精品三级| 麻豆久久精品国产亚洲av | 精品久久久精品久久久| 村上凉子中文字幕在线| 国产成人精品无人区| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕精品免费在线观看视频| 亚洲国产看品久久| www.www免费av| 天堂影院成人在线观看| 亚洲成人久久性| 老汉色av国产亚洲站长工具| 三级毛片av免费| 在线国产一区二区在线| 国产精品久久久av美女十八| 欧美日韩亚洲综合一区二区三区_| aaaaa片日本免费| 99久久99久久久精品蜜桃| 久久久久亚洲av毛片大全| 一个人免费在线观看的高清视频| 欧美中文日本在线观看视频| 欧美精品啪啪一区二区三区| 1024香蕉在线观看| 国产极品粉嫩免费观看在线| 一级毛片精品| av天堂在线播放| 淫秽高清视频在线观看| 成人国语在线视频| 在线十欧美十亚洲十日本专区| 国产精品久久久av美女十八| 亚洲黑人精品在线| 免费久久久久久久精品成人欧美视频| 国产亚洲精品久久久久久毛片| 亚洲视频免费观看视频| 久久久久九九精品影院| 国产1区2区3区精品| 亚洲 国产 在线| 成年人免费黄色播放视频| 欧美日韩视频精品一区| 麻豆一二三区av精品| 18禁观看日本| 欧美日韩乱码在线| 欧美成人免费av一区二区三区| 国产精品久久久久成人av| 美女大奶头视频| 欧美日韩乱码在线| 在线播放国产精品三级| 国产一区二区三区视频了| 国产精品免费一区二区三区在线| 国产精品美女特级片免费视频播放器 | 欧美另类亚洲清纯唯美| 日本黄色视频三级网站网址| 日本黄色日本黄色录像| 精品一区二区三区四区五区乱码| 久久久久国产精品人妻aⅴ院| 午夜福利欧美成人| 99riav亚洲国产免费| 国产日韩一区二区三区精品不卡| 欧美一区二区精品小视频在线| videosex国产| av网站在线播放免费| 午夜91福利影院| 纯流量卡能插随身wifi吗| 国产成年人精品一区二区 | 久久人人精品亚洲av| 国内毛片毛片毛片毛片毛片| 中出人妻视频一区二区| 又紧又爽又黄一区二区| 99国产精品免费福利视频| 桃红色精品国产亚洲av| 亚洲中文日韩欧美视频| 亚洲成人免费av在线播放| 日日干狠狠操夜夜爽| 免费久久久久久久精品成人欧美视频| 亚洲精品在线美女| 国产一区二区三区在线臀色熟女 | 国产蜜桃级精品一区二区三区| 在线观看免费日韩欧美大片| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看| 这个男人来自地球电影免费观看| 亚洲精品国产色婷婷电影| 夜夜看夜夜爽夜夜摸 | 老司机午夜福利在线观看视频| 老熟妇乱子伦视频在线观看| 免费高清视频大片| 国产成人欧美在线观看| 一级a爱视频在线免费观看| 宅男免费午夜| 久久久国产成人免费| 精品高清国产在线一区| 啦啦啦免费观看视频1| 久久天躁狠狠躁夜夜2o2o| 波多野结衣一区麻豆| 中文字幕人妻熟女乱码| 80岁老熟妇乱子伦牲交| 亚洲熟妇熟女久久| a级毛片在线看网站| 欧美成狂野欧美在线观看| 首页视频小说图片口味搜索| 91老司机精品| 午夜91福利影院| 久久婷婷成人综合色麻豆| 男女下面进入的视频免费午夜 | 亚洲第一欧美日韩一区二区三区| 看免费av毛片| 少妇 在线观看| 久久人人精品亚洲av| 日韩三级视频一区二区三区| 国产有黄有色有爽视频| 亚洲精华国产精华精| 久久久久国产一级毛片高清牌| 麻豆成人av在线观看| 美女高潮到喷水免费观看| 久久精品国产99精品国产亚洲性色 | 看黄色毛片网站| 国产真人三级小视频在线观看| a在线观看视频网站| 国产成人免费无遮挡视频| 搡老岳熟女国产| 免费看十八禁软件| 在线观看舔阴道视频| 在线观看免费视频网站a站| 成人黄色视频免费在线看| 美国免费a级毛片| 国产精品国产av在线观看| 久久香蕉激情| 精品福利观看| 黄色视频不卡| 免费av中文字幕在线| 悠悠久久av| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区在线不卡| 国产男靠女视频免费网站| av超薄肉色丝袜交足视频| 窝窝影院91人妻| 巨乳人妻的诱惑在线观看| 在线观看免费日韩欧美大片| 欧美精品啪啪一区二区三区| 国产一区二区激情短视频| 久久久久国产一级毛片高清牌| 亚洲精品在线美女| 久久亚洲真实| 精品久久久精品久久久| 国产精品一区二区三区四区久久 | 日本免费a在线| 国产av一区在线观看免费| ponron亚洲| 久久99一区二区三区| 精品一区二区三区av网在线观看| 精品久久久久久,| 桃红色精品国产亚洲av| 无限看片的www在线观看| 欧美激情 高清一区二区三区| 免费不卡黄色视频| 身体一侧抽搐| 电影成人av| 91精品国产国语对白视频| 精品国产一区二区三区四区第35| 午夜精品在线福利| 欧美精品亚洲一区二区| 久久中文字幕一级| 啦啦啦在线免费观看视频4| 午夜福利欧美成人| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| 嫩草影视91久久| 黄色丝袜av网址大全| 午夜日韩欧美国产| 国产精品98久久久久久宅男小说| 91字幕亚洲| 久久人人精品亚洲av| 欧美激情极品国产一区二区三区| 一个人观看的视频www高清免费观看 | 中文字幕av电影在线播放| 999久久久精品免费观看国产| 欧美午夜高清在线| 麻豆av在线久日| e午夜精品久久久久久久| 国产亚洲精品一区二区www| 中文亚洲av片在线观看爽| 精品国产美女av久久久久小说| 69av精品久久久久久| 午夜成年电影在线免费观看| 午夜精品国产一区二区电影| 丝袜人妻中文字幕| 不卡av一区二区三区| 99国产精品免费福利视频| 国产精品野战在线观看 | 亚洲精品av麻豆狂野| 精品一品国产午夜福利视频| 亚洲一区高清亚洲精品| 国产精品1区2区在线观看.| 视频区欧美日本亚洲| 国产成人啪精品午夜网站| 成人国产一区最新在线观看| 午夜精品在线福利| 国产亚洲av高清不卡| 久久中文字幕一级| 在线国产一区二区在线| 亚洲精品粉嫩美女一区| 黄色a级毛片大全视频| 国产又色又爽无遮挡免费看| 最新在线观看一区二区三区| 麻豆久久精品国产亚洲av | 一进一出抽搐gif免费好疼 | 乱人伦中国视频| 国产成人影院久久av| 咕卡用的链子| 久久精品影院6| 夫妻午夜视频| 日韩 欧美 亚洲 中文字幕| bbb黄色大片| 亚洲五月色婷婷综合| 丝袜美腿诱惑在线| 在线观看免费高清a一片| 亚洲片人在线观看| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 制服诱惑二区| 欧美精品亚洲一区二区| 99国产精品一区二区蜜桃av| 色播在线永久视频| 欧美黑人欧美精品刺激| 夫妻午夜视频| 精品一品国产午夜福利视频| 色综合欧美亚洲国产小说| 精品电影一区二区在线| 欧美成狂野欧美在线观看| 97人妻天天添夜夜摸| 国产精品av久久久久免费| 99精品在免费线老司机午夜| 国产蜜桃级精品一区二区三区| 国产精品亚洲av一区麻豆| 99久久久亚洲精品蜜臀av| 18美女黄网站色大片免费观看| 国产成+人综合+亚洲专区| 美女高潮到喷水免费观看| 国产精品一区二区免费欧美| 黄片小视频在线播放| 亚洲avbb在线观看| 午夜福利欧美成人| 精品午夜福利视频在线观看一区| 亚洲国产精品999在线| 中文字幕人妻熟女乱码| 日韩中文字幕欧美一区二区| 一级,二级,三级黄色视频| 一本综合久久免费| 极品教师在线免费播放| 精品卡一卡二卡四卡免费| 日本黄色日本黄色录像| 精品久久久久久久久久免费视频 | 国产97色在线日韩免费| 人人妻,人人澡人人爽秒播| 如日韩欧美国产精品一区二区三区| 九色亚洲精品在线播放| 午夜亚洲福利在线播放| 国产伦一二天堂av在线观看| 日韩视频一区二区在线观看| 亚洲国产欧美网| 99在线人妻在线中文字幕| 久久人妻av系列| 日韩欧美一区视频在线观看| 熟女少妇亚洲综合色aaa.| 日韩成人在线观看一区二区三区| 国产亚洲欧美98| 天天躁夜夜躁狠狠躁躁| av福利片在线| 成人av一区二区三区在线看| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 欧美成狂野欧美在线观看| 精品少妇一区二区三区视频日本电影| 老鸭窝网址在线观看| 不卡av一区二区三区| 极品教师在线免费播放| 美国免费a级毛片| 露出奶头的视频| 中文字幕av电影在线播放| 亚洲久久久国产精品| 国产高清国产精品国产三级| 国产91精品成人一区二区三区| 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| 桃色一区二区三区在线观看| 久久国产乱子伦精品免费另类| 午夜老司机福利片| 韩国精品一区二区三区| 日韩有码中文字幕| 大型黄色视频在线免费观看| 精品一区二区三区四区五区乱码| 亚洲少妇的诱惑av| 国产97色在线日韩免费| 国产男靠女视频免费网站| 国产成人一区二区三区免费视频网站| 亚洲av成人不卡在线观看播放网| 精品久久久久久电影网| 涩涩av久久男人的天堂| 国产单亲对白刺激| 成人三级黄色视频| 熟女少妇亚洲综合色aaa.| 欧美在线黄色| 国产精品久久电影中文字幕| 黄色女人牲交| 伦理电影免费视频| 国产一卡二卡三卡精品| 淫妇啪啪啪对白视频| 岛国视频午夜一区免费看| 国产人伦9x9x在线观看| 深夜精品福利| 国产激情欧美一区二区| 久久久水蜜桃国产精品网| 亚洲精品久久午夜乱码| 国产一区二区激情短视频| 日日摸夜夜添夜夜添小说| 欧美日本中文国产一区发布| 精品高清国产在线一区| 叶爱在线成人免费视频播放| 在线十欧美十亚洲十日本专区| 男女之事视频高清在线观看| 老司机在亚洲福利影院| 亚洲国产欧美一区二区综合| 国产亚洲av高清不卡| www.www免费av| 亚洲人成网站在线播放欧美日韩| 桃红色精品国产亚洲av| 悠悠久久av| 在线观看日韩欧美| 亚洲av成人不卡在线观看播放网| 欧美丝袜亚洲另类 | 一夜夜www| 欧美乱色亚洲激情| 国产伦人伦偷精品视频| a级毛片黄视频| 一个人观看的视频www高清免费观看 | 国产三级黄色录像| 亚洲一区二区三区欧美精品| 久久精品亚洲精品国产色婷小说| 十八禁网站免费在线| 人妻久久中文字幕网| 神马国产精品三级电影在线观看 | 搡老乐熟女国产| 国产又色又爽无遮挡免费看| 日本撒尿小便嘘嘘汇集6| 亚洲精品中文字幕在线视频| 久久亚洲真实| 久久久国产成人精品二区 | 18美女黄网站色大片免费观看| 午夜福利影视在线免费观看| 黄片小视频在线播放| 精品久久久久久成人av| 亚洲av成人不卡在线观看播放网| 18禁裸乳无遮挡免费网站照片 | 亚洲欧美一区二区三区黑人| 91老司机精品| 日本三级黄在线观看| 超色免费av| 搡老乐熟女国产| 国产不卡一卡二| 午夜老司机福利片| 精品国产超薄肉色丝袜足j| 亚洲七黄色美女视频| 中文字幕av电影在线播放| 在线免费观看的www视频| 中出人妻视频一区二区| 大型黄色视频在线免费观看| 黄片大片在线免费观看| 色婷婷久久久亚洲欧美| 嫩草影院精品99| 国产成人av激情在线播放| av在线播放免费不卡| bbb黄色大片| 精品国产超薄肉色丝袜足j| 亚洲黑人精品在线| 国产伦一二天堂av在线观看| 日韩三级视频一区二区三区| 制服人妻中文乱码| 天天影视国产精品| 午夜老司机福利片| 精品一区二区三区视频在线观看免费 | 精品国产乱码久久久久久男人| 人人妻人人添人人爽欧美一区卜| 国产一区在线观看成人免费| 美女 人体艺术 gogo| 欧美日韩黄片免| 亚洲国产欧美一区二区综合| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线| 在线av久久热| 女人被躁到高潮嗷嗷叫费观| a在线观看视频网站| 成年版毛片免费区| 日本免费a在线| netflix在线观看网站| 午夜免费观看网址| 日韩人妻精品一区2区三区| 久久青草综合色| 在线免费观看的www视频| 日本三级黄在线观看| 国产精品久久久av美女十八| x7x7x7水蜜桃| 欧美一区二区精品小视频在线| 后天国语完整版免费观看| 成人特级黄色片久久久久久久| 久久天躁狠狠躁夜夜2o2o| 精品福利观看| 亚洲av熟女| 99国产精品一区二区三区| 丝袜美腿诱惑在线| 91在线观看av| 久久精品人人爽人人爽视色| 夜夜躁狠狠躁天天躁| 欧美黄色片欧美黄色片| 黄色丝袜av网址大全| 日韩欧美一区视频在线观看| 婷婷丁香在线五月| 久久人人精品亚洲av| 日韩一卡2卡3卡4卡2021年| 纯流量卡能插随身wifi吗| 如日韩欧美国产精品一区二区三区| 欧美成人午夜精品| 满18在线观看网站| 99久久99久久久精品蜜桃| 91精品三级在线观看| 国产日韩一区二区三区精品不卡| 日韩精品免费视频一区二区三区| e午夜精品久久久久久久| 人人妻,人人澡人人爽秒播| 99久久国产精品久久久| ponron亚洲| 一进一出抽搐动态| 999久久久精品免费观看国产| 日韩大码丰满熟妇| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 亚洲自偷自拍图片 自拍| 精品久久久久久久毛片微露脸| 窝窝影院91人妻| 亚洲精品美女久久久久99蜜臀| 国产高清videossex| 久久亚洲精品不卡| 午夜久久久在线观看| 久久狼人影院| 成年女人毛片免费观看观看9| 新久久久久国产一级毛片| 99久久精品国产亚洲精品| 波多野结衣高清无吗| 日本vs欧美在线观看视频| 伊人久久大香线蕉亚洲五| 18禁黄网站禁片午夜丰满| 亚洲成人久久性| 欧美色视频一区免费| 国产99白浆流出| 两性夫妻黄色片| 久久精品91蜜桃| 黄片大片在线免费观看| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 亚洲色图av天堂| av网站免费在线观看视频| 久久久久国产精品人妻aⅴ院| 久久人妻熟女aⅴ| 女人爽到高潮嗷嗷叫在线视频| 丝袜美腿诱惑在线| 国产亚洲精品一区二区www| 少妇裸体淫交视频免费看高清 | 精品久久久精品久久久| av在线天堂中文字幕 | 日韩成人在线观看一区二区三区| 夜夜夜夜夜久久久久| 日韩成人在线观看一区二区三区| 国产97色在线日韩免费| 亚洲美女黄片视频| 日本撒尿小便嘘嘘汇集6| 免费看a级黄色片| 夜夜夜夜夜久久久久| 亚洲国产精品sss在线观看 | 亚洲av五月六月丁香网| 夜夜爽天天搞| 亚洲五月色婷婷综合| 日韩大码丰满熟妇| 成人特级黄色片久久久久久久| 欧美乱妇无乱码| 12—13女人毛片做爰片一| 中文字幕人妻熟女乱码| ponron亚洲| 精品国产一区二区久久| 新久久久久国产一级毛片| 首页视频小说图片口味搜索| 丁香欧美五月| 成人三级做爰电影| 男人舔女人下体高潮全视频| 国产区一区二久久| 亚洲成人免费av在线播放| 国产欧美日韩一区二区三区在线| 国产欧美日韩一区二区精品| 女人被躁到高潮嗷嗷叫费观| 人人妻人人澡人人看| 91字幕亚洲| 夜夜夜夜夜久久久久| 丁香欧美五月| √禁漫天堂资源中文www| 1024香蕉在线观看| 精品国产乱子伦一区二区三区| 亚洲精品美女久久av网站| 日韩大尺度精品在线看网址 | 欧美日韩乱码在线| 手机成人av网站| 最近最新免费中文字幕在线| 午夜激情av网站| 桃红色精品国产亚洲av| 99久久久亚洲精品蜜臀av| 精品国产一区二区久久| 欧美日本亚洲视频在线播放| 91麻豆av在线| 在线永久观看黄色视频| 黄色a级毛片大全视频| 在线天堂中文资源库| netflix在线观看网站| 精品免费久久久久久久清纯| 老司机在亚洲福利影院| 国产成人影院久久av| 老司机靠b影院| 精品国产一区二区久久| 在线视频色国产色| 99久久综合精品五月天人人| 国产主播在线观看一区二区| 成人影院久久| 在线天堂中文资源库| 在线国产一区二区在线| 18禁美女被吸乳视频| 亚洲av五月六月丁香网| 在线观看免费日韩欧美大片| 少妇粗大呻吟视频| 桃色一区二区三区在线观看| 国产精品香港三级国产av潘金莲| 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻丝袜一区二区| 国产精品电影一区二区三区| 成年人免费黄色播放视频| 色播在线永久视频| а√天堂www在线а√下载| 变态另类成人亚洲欧美熟女 | 欧美人与性动交α欧美精品济南到| 黑人巨大精品欧美一区二区蜜桃| 18禁黄网站禁片午夜丰满| 精品人妻1区二区| 亚洲avbb在线观看| 精品乱码久久久久久99久播| 精品国内亚洲2022精品成人| 淫秽高清视频在线观看| 欧美性长视频在线观看| 国产黄色免费在线视频| 国产精品综合久久久久久久免费 | 黄片大片在线免费观看| 成人国产一区最新在线观看| 亚洲欧美精品综合久久99| 18禁国产床啪视频网站| 免费一级毛片在线播放高清视频 | 欧美日韩乱码在线| 欧美乱色亚洲激情| 久久亚洲真实| 999久久久精品免费观看国产| 很黄的视频免费| 国产精品香港三级国产av潘金莲| 最新美女视频免费是黄的| 午夜福利免费观看在线| 日韩精品中文字幕看吧|