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

    Numerical and experimental studies of hydrodynamics of flapping foils *

    2018-05-14 01:42:52KaiZhou周凱JunkaoLiu劉軍考WeishanChen陳維山
    關(guān)鍵詞:劉軍

    Kai Zhou (周凱), Jun-kao Liu (劉軍考), Wei-shan Chen (陳維山)

    1. College of Mechanical and Electronic Engineering, Shandong Agricultural University, Taian 271018, China

    2. State Key Laboratory of Robotics and System, Harbin Institute of Technology, Harbin 150001, China

    Introduction

    The conventional screw propellers have the disadvantages of low efficiency, poor maneuvering performance, and large noise , which greatly limit their applications in narrow, complex and dynamic environments. However, after millions of years of evolution,creatures are capable of swimming or flying with high speed and efficiency. Therefore, the biomimetic propulsor attracts research attentions[1].

    Most fast-swimming creatures undulate their fins to obtain thrust. Similarly, flying creatures undulate their wings to resist gravity and move forward. They have similar characters in shape and motion mode, in reducing resistance and obtaining valuable force. The flapping foil based on bionics is a sort of simplified models, which can imitate the motion of the wings or fins of fish and birds[2-6]. With different motion parameters, the hydrodynamic performance varies correspondingly. In addition, a simple imitation of motion is not enough, its mechanism is a more important issue.The hydrodynamic mechanism would provide a theoretical reference for the design of underwater vehicles based on the flapping propulsion.

    In the movements of fish and birds, one sees the symmetric mode and the asymmetric mode. Most fish flap their caudal fins in the symmetric mode, with equal forces generated from both the upstroke and down stroke. On the other hand, with an oscillation parallel to the advance direction, the asymmetrical mode consists of a powerful down stroke generating large force and a weak upstroke with a weak force, as evidenced by the motion patterns of turtles and birds[7-11].

    The symmetric mode and the asymmetric mode were usually studied separately, focusing on certain bionic prototypes, lack of generality. The kinematic parameters of living fish, birds and bionic robots were widely studied, but without enough attention to the hydrodynamic mechanism. In this paper, a united motion model containing the symmetric mode and the asymmetric mode is developed. By analyzing the flow field generated in the process of flapping and forcing of the foil, the propulsion mechanism of the bionic foil is studied. And then the influence of the motion parameters on the hydrodynamic performance of the bionic foil is analyzed. The simulation results are compared with the existing related studies to assess the reliability of the numerical method. In addition, an experimental method is adopted to further validate the numerical approach.

    1. Materials and methods

    1.1 Geometric model

    Based on observations of fish and birds, the sections of fins or wings parallel to the flow direction take streamlined shapes, similar to the NACA airfoils.Thus the hydrodynamic force can be viewed as the resultant force of a series of forces acting on the streamlined sections. In this paper, the NACA0013 airfoil is chosen to construct the section. Figure 1 presents the geometry of the three-dimensional model,wherec=0.1m is the chord length,L=0.3m is the length along the span direction, andOis the reference point.Ois located 0.25cfrom the leading edge and the projected area (one-sided)Sis 3.2×10-2m2. The foil makes a periodic motion along the vertical direction and the horizontal direction,meanwhile, the foil makes a pitching movement around the reference axis.

    Fig.1 (Color online) Geometric model of the foil

    1.2 Motion model

    In this study, the foil is allowed to move with three degrees of freedom. Specifically, the foil is towed forward at a constant speedUand there are three types of movements: the heave motion transverse to the direction of towing, the angular motion around a span wise axis, and the surge motion parallel to the direction of towing. The motion of the foil can be described by the following equations:

    wherefis the flapping frequency,tis the time,and β is the stroke angle[12]. In contrast to the traditional symmetric motion model,x(t) is adopted to describe the oscillatory motion parallel to the forward direction. With the stroke angle β, different trajectories can be formed to represent different motion models of different bionic prototypes. As shown in Fig.2, the case of β=45° describes the motion mode of birds. It can be seen that birds direct their wings forward during the down stroke, creating a highly asymmetric flap to generate a net lift force. The case of β=90° describes the motion mode of fish.In the symmetric mode, a symmetric trajectory is formed without the surge motion parallel to the direction of towing. The case of β=135° describes the motion mode of turtles. It can be seen that turtles direct their flippers backward during the down stroke,creating a highly asymmetric flap to generate a net trust force.

    In addition, the pitching angle θ is given by

    where

    The angle α(t) is the angle of attack. In the case of β =90°

    and in the case of β≠90°

    The anglemθ represents the tangential direction of the trajectory, and the angle of attack α(t) is also an independent variable. The definition makes it easier to develop an optimization method for the motion trajectory or the angle of attack. In the nature, fish and birds adjust the angle of attack and the trajectory constantly to eliminate the unwanted forces and keep the speed and the efficiency at a high level.

    1.3 Numerical method

    Using the commercial software Fluent (ANSYS Inc., Canonsburg, USA), the governing equations are discretized by the finite volume method (FVM) with a first-order discretization in the time and a second order discretization in the space. The fluid motion is computed based on an unstructured, tetrahedral grid.The simulation is a dynamic process and the boundaries of the foil are changing constantly. Therefore,in this paper, the dynamic mesh method is adopted to deal with the moving boundary problem. The motion of the foil is realized through user-defined functions(UDFs) linked to the Fluent. In the calculation, the SIMPLE method of the pressure-speed amending method and thek-εturbulence mode are adopted.

    The simulations start from the static state of the foil, and then the foil moves according to the motion model. The cubic computational domain (2.5 m×1.2 m×1.2 m) is discretized into 8.1×105grid cells.Besides, refined meshes are used near the boundaries of the foil. To verify the grid in dependency, the simulation with a mesh of smaller cell size is carried out and only a difference of 4% between the coarse grid and the fine grid is found.

    To evaluate the hydrodynamic performance of the foil, the force and the torque acting on the foil are extracted and converted to dimensionless forms:

    where ρ is the density of the fluid,Fyis the transverse force, -Fxis the thrust force,Mθis the torque, andSis the projected area of the foil. The propulsive efficiency η is defined as the ratio of the output poweroPto the input powereP, which can be defined as:

    whereTis the period of the motion. The chord Reynolds number and the Strouhal numberStare defined as:

    where ν is the fluid kinematic viscosity.

    2. Experiment approach and setup

    With the developments of computer technology and computational fluid dynamics (CFD), numerical methods are having more and more extensive applications. Compared with the traditional experimental method, the numerical method has the advantage of high repeatability, great maneuverability, and less input requirement. Moreover, it is easier to obtain the details of the flow field with the numerical method,thus, many difficulties of measurement in the experimental method can be avoided. Therefore, the numerical method becomes an important method in the mechanism research. However, due to the limitations of the theoretical model, the reliability of the numerical results still needs to be improved. In addition, the computation parameters of the numerical models need to be modified and confirmed according to the experimental results. So usually both the numerical method and the experimental method are adopted in studies.

    To validate the reliability of the numerical method used in this paper, an experiment platform is designed and the verification experiments are carried out. Figure 3 shows the design schematic diagram of the experiment platform.

    To realize the movements of three degrees of freedom, three servo motors equipped with reducers are used. The synchronous belt transmission is adopted in the towing direction and the gear-gear stripe transmission is adopted in the heave direction.All servo motors are controlled by computer.

    Fig.3 (Color online) Design schematic diagram of the experiment platform

    Fig.4 (Color online) Measuring module of the experiment platform

    Fig.5 (Color online) Actual photo of the experiment platform

    Except the power module and the transmission module, the measuring module is an important part of the system. As is shown in Fig.4, two slide rails are installed on the moving platform along the heave direction. The force sensors act as the transmission components to supply power. According to the mechanical analysis of the moving platform in the horizontal plane, there are the hydrodynamic force,the frictional force, the traction force and the inertia force. Because the coefficients of friction of the slide rails are very small, the frictional force can be ignored.The inertia force can be calculated through the motion equations. The traction force can be provided by the force sensor. Through mathematical calculations, the hydrodynamic force acting on the foil can be obtained.Figure 5 shows the actual photo of the experiment platform.

    3. Results

    Firstly, for comparing with the existing related studies[13, 14], the motion parameters of the foil are selected as follows:Re=11000,St=0.3,h=c=0.055 m, αmax=25° and the stroke angle β=45°,90° and 135°. Figure 6 shows the force coefficients at different stroke angles.

    Fig.6 Force coefficients at different stroke angles

    The curves see a periodic behaviour, and the curves in a cycle are shown in Fig.6. In the case of β =90°, the curve of the thrust force coefficient has the same shape in the upstroke and the down stroke because of the symmetrical motion. Note that the mean value ofCx, which is equal to 0.58, is positive,while the mean value ofCyis almost zero. In other words, the foil obtains the thrust force without the net transverse force. In the case of β=45°, one sees small fluctuations at the beginning and the end of the upstroke when the foil rotates quickly. The mean thrust force coefficient is 0.21 and the mean transverse force coefficient is 2.12. In the case of β=135°, the mean thrust force coefficient in a cycle is 0.12. In the cases of β=45° and β=135°, the mean transverse forces are not equal to zero because of the asymmetrical motion. In the case of β=45°, it is the transverse force that increases significantly rather than the thrust force. The thrust force in the case of β =135° is relatively small because of the negative mean thrust force of the upstroke, and other experiments show similar phenomena during the upstroke,called “the memory effects” in the wake[13]. Besides,the propulsive efficiency of symmetrical mode is about 47.0% which is well within the range of fish.The propulsive efficiency in the case of β=45° is equal to 11.3%, which is low because the for wards biased flap increases the transverse force significantly rather than the thrust force. The propulsive efficiency in the case of β=135° case is equal to 31.0%,which is also relatively low because of the unwanted energy expenditure during the upstroke.

    Figure 7 shows the vortex pattern in the wake. It can be seen that there are vortex rings in the wake,which is similar to the wake of fish and birds. The jet stream can be induced among the vortex rings and the foil obtains a reaction force correspondingly. Specifically, Figure 7(a) shows the 3-D vortex pattern in the case of β=90°. The vortex rings consist of a train of inverted hairpin-like vortices braided together such that the legs of each vortex are attached to the head of the preceding vortex. The vortex structure is similar to those found in the nature and the simulations in other studies[15-18]. Figures 7(b), 7(c) show the comparisons of the vortex patterns in different motion modes viewed fromzdirection. It is shown that the vortex patterns of the asymmetrical mode are characterized by irregular and discontinuous structures. As a result,because of the asymmetrical motion , a net transverse force is produced and more useless energy is expended.

    In the nature, birds direct their wings forward during the down stroke, creating a highly asymmetric flap that can generate a net lift force to support the weight. Turtles direct their flippers backward during the down stroke, creating a highly asymmetric flap that can reduce the oscillatory force. The fins of fish can form a symmetric trajectory without the surge motion parallel to the direction of towing and keep a high propulsive efficiency. In the following parts, the influence of the motion parameters on the hydrodynamic performance is analyzed. For comparing with the experimental results, the chord length is chosen to be 0.1 m and the length along the span direction is 0.3 m. The other parameters are selected as follows: αmax=25°,h=0.05m,U=0.2 m/s ,St=0.3. The force coefficients represent the mean value in a cycle.

    Fig.7 (Color online) Vortex patterns for different stroke angles

    3.1 Influence of the angle of attack

    To study the influence of the angle of attack, its value is chosen to vary from 10° to 45°. The definition of the angle of attack is different in the symmetric and asymmetric modes, so the influence of the angle of attack for different stroke angles should be analyzed,separately. Figures 8-10 show the mean force coefficients and the propulsive efficiency of the foil with different stroke angles. In the symmetrical mode,the mean thrust force keeps a relatively large value when the angle of attack αmaxvaries from 20° to 30°,as is consistent with fish in the nature. Furthermore, as the angle of attack increases, the mean thrust force increases correspondingly and reaches the peak point,then starts to decrease. Besides, it is found that the propulsive efficiency sees a similar trend as compared with the mean thrust force. Larger or smaller αmaxwill bring a large transverse force and a reduced thrust force, and then decrease the propulsive efficiency. In the case of β=45°, the trajectory of the down stroke is more gradual than in other cases, so large αmaxadopted in the powerful down stroke can increase the transverse force significantly and have little effect on the thrust force. As a result, the propulsive efficiency decreases significantly. In contrast, in the case of β= 135°, the trajectory of the down stroke is steeper,so large αmaxis allowed to keep the thrust force.However, the strong angular motion with a large αmaxwill increase the useless energy consumption, to reduce the propulsive efficiency.

    Fig.8 Performance curves against αmax with β=45°

    Fig.9 Performance curves against αmax with β=90°

    Fig.10 Performance curves against αmax with β=135°

    3.2 Influence of the stroke angle

    The performance curves against the stroke angle are shown in Fig.11. The dotted lines in Fig.11 represent the values in the case of β=90°. It is shown that the stroke angle greatly affects the hydrodynamic performance of the foil. By choosing a large stroke angle, the foil can obtain a large thrust force with a small transverse force. In contrast, with a small stroke angle, the foil can obtain both a large transverse force and a large thrust force, but the propulsive efficiency is relatively low. In addition, a stroke angle around 90° can achieve a smaller thrust force with a high propulsive efficiency. However, in the symmetrical mode, a higher efficiency can be achieved as com-pared with the asymmetrical mode. A reasonable stroke angle should be adopted in different engineering applications.

    Fig.11 Performance curves against stroke angles

    3.3 Influence of the heave motion

    Here, the case of β=90° is chosen as a typical example to show the influence of the heave motion on the hydrodynamic performance. Figure 12 shows the force coefficients and the propulsive efficiency versus the amplitude of the heave motion. The Strouhal number is chosen to be 0.3 to keep a high efficiency.It is shown that the impact of the amplitude on the mean thrust force is small. Note that the propulsive efficiency reaches the peak value aroundh/c≈ 0.8.As is known, a large heave motion obviously will lead to a large thrust force but a larger heave motion will lead to a lower frequency to keep the Strouhal number constant. So the change of the mean thrust is insignificant when the amplitude of the heave motionhincreases.

    Fig.12 Performance curves against amplitudes of heave motion

    3.4 Influence of Strouhal number

    The Strouhal number is an important parameter in the biological swimming. Furthermore, most fish are shown to swim near a “universal” optimal valueSt=0.3 because of the high propulsive efficiency[19,20].Figure 13 shows the performance curves against the Strouhal number when β=90°. Here, the frequency is chosen to vary to obtain different Strouhal numbers. It is shown that the propulsive efficiency keeps a relatively large value within the range of 0.2-0.4, as is consistent with the range of fish in the nature. Furthermore, a large Strouhal number can lead to a large thrust force, but the transverse force increases correspondingly.

    Figure 14 shows the vortex pattern in the wake for different values ofStwhen β=90°. At a low value ofSt, the 3-D structure of the wake is a single and continuous vortex ring. At a high value ofSt, a more complex and laterally diverging structure is observed in the wake. Smaller-scale vertical structures attaching to the hairpin vortices are observed in the wake. The energy dissipation is greater in this case, so the propulsive efficiency is lower at a high value ofSt.In addition, the more complex and laterally diverging structure impacts the transverse force and a considerable net transverse force is produced. The results of asymmetric motion models are similar to those of symmetric motion models.

    Fig.13 Performance curves against Strouhal numbers

    Fig.14 (Color online) Performance curves against amplitudes of heave motions

    3.5 Experimental results

    To validate the reliability of the numerical method used in this paper, verification experiments are carried out. Table 1 shows the numerical results and experimental results with same parameters. Figure 15 shows the photo of the experiment site. The selected motion parameters of the foil are as follows:Re=20000,St=0.3,h=0.05 m and αmax=25°.Through a comparison, it is found that the present results compare well with those in the previous studies.Therefore, the adopted numerical method is reliable.

    Table 1 Comparison of results

    Fig.15 (Color online) Photo of experiment site

    4. Conclusions

    In this paper, a universal kinematic model with three degrees of freedom is adopted. The motion parallel to the flow direction considered in the model can produce trajectories of different prototypes. Based on this kinematic model, numerical simulations of the flapping foil with different trajectories are conducted.In certain parameter ranges, the mean thrust forces of different motion modes in a cycle could have positive values, that is, the foil can obtain a thrust force. And the propulsive efficiency of the symmetric motion model is significantly greater than those of the asymmetric motion models. In addition, the influences of the motion parameters are analyzed. It is found that the motion parameters play important roles in the hydrodynamic performance of the flapping foil.Besides, to validate the reliability of the numerical method, an experiment platform is designed and verification experiments are carried out. The consistence between the numerical results and the experimental results indicates that the numerical method of the flapping foil is reasonable.

    The results can help the bionic propulsor design and improvement. Furthermore, through a reasonable optimization algorithm of the motion parameters, the flapping foil can either reduce the oscillatory force or increase the thrust force and the propulsive efficiency.The optimization algorithm will be completed in the later work.

    Acknowledgement

    This work was supported by the State Key Laboratory of independent topics, Harbin Institute of Technology (Grant No. SKLRS200801C).

    [1] Wu T. Y. Fish swimming and bird/insect flight [J].Annual Review of Fluid Mechanics, 2011, 43: 25-58.

    [2] Tang C., Lu X. Y. Self-propulsion of a three-dimensional flapping flexible plate [J].Journal of Hydrodynamics,2016, 28(2): 1-9.

    [3] Han R., Zhang J., Cao L. et al. Propulsive performance of a passively flapping plate in a uniform flow [J].Journal of Hydrodynamics, 2015, 27(4): 496- 501.

    [4] Hu J., Xiao Q. Three-dimensional effects on the translational locomotion of a passive heaving wing [J].Journal of Fluids and Structures, 2014, 46: 77-88.

    [5] Cheng H. Y., Long X. P, Ji B. et al. Numerical investigation of unsteady cavitating turbulent flows around twisted hydrofoil from the Lagrangian viewpoint [J].Journal of Hydrodynamics, 2016, 28(4): 709-712.

    [6] Liu W., Li N., Zhao J. et al. Wake structure and hydrodynamic performance of flapping foils mimicking fish fin kinematics [J].Saudi Journal of Biological Sciences, 2017,24(6): 1344-1354.

    [7] Szymik B. G., Satterlie R. A. Changes in wing stroke kinematics associated with a change in swimming speed in a ptero pod mollusk, Clione limacine [J].The Journal of Experimental Biology, 2011, 214(23): 3935-3947.

    [8] Licht S. C., Wibawa M. S., Hover F. S. et al. In-line motion causes high thrust and efficiency in flapping foils that use power down stroke [J].The Journal of Experimental Biology, 2010, 213(1): 63-71.

    [9] Tobalske B. W., Warrick D. R., Clark C. J. et al. Three dimensional kinematics of hummingbird flight [J].The Journal of Experimental Biology, 2007, 210(13):2368-2382.

    [10] Dileo C., Deng X. Design of and experiments on a dragonfly-inspired robot [J].Advanced Robotics, 2009, 23:1003-1021.

    [11] Berg A. M., Bie wener A. A. Wing and body kinematics of takeoff and landing flight in the pigeon (Columba livia) [J].The Journal of Experimental Biology, 2010, 213(10):1651-1658.

    [12] Lindhe Norberg U. M., Winter Y. Wing beat kinematics of a nectar-feeding bat, Glossophaga soricina, flying at different flight speeds and Strouhal numbers [J].The Journal of Experimental Biology, 2006, 209(19): 3887-3897.

    [13] Izraelevitz J. S., Triantafyllou M. S. Adding in-line motion and model-based optimization offers exceptional force control authority in flapping foils [J].Journal of Fluid Mechanics, 2014, 742: 5-34.

    [14] Zhou K., Liu J., Chen W. Numerical study on hydrodynamic performance of bionic caudal fin [J].Applied Sciences-Basel, 2016, 6(1): 15.

    [15] Xia D., Chen W., Liu J. et al. The three-dimensional hydrodynamics of thunniform swimming under self-propulsion [J].Ocean Engineering, 2015, 110: 1-14.

    [16] Bhalla A. P. S., Griffith B. E., Patankar N. A. A forced damped oscillation framework for undulatory swimming provides new insights into how propulsion arises in active and passive swimming [J].PLOS Computational Biology,2013, 9(6): e1003097.

    [17] Neveln I. D., Bale R., Bhalla A. P. S. et al. Undulating fins produce off-axis thrust and flow structures [J].Journal of Experimental Biology, 2014, 217(2): 201-213.

    [18] Bhalla A. P. S., Bale R., Griffith B. E. et al. A unified mathematical framework and an adaptive numerical method for fluid–structure interaction with rigid, deforming, and elastic bodies [J].Journal of Computational Physics, 2013, 250: 446-476.

    [19] Taylor G. K., Nudds R. L., Thomas A. L. R. Flying and swimming animals cruise at a Strouhal number tuned for high power efficiency [J].Nature, 2003, 425: 707-711.

    [20] Thomas A. L. R., Taylor G. K., Srygley R. B. et al.Dragonfly flight: Free-flight and tethered flow visualizations reveal a diverse array of unsteady lift-generating mechanisms, controlled primarily via angle of attack [J].Journal of Experimental Biology, 2004, 207(24):4299-4323.

    猜你喜歡
    劉軍
    劉軍作品
    大眾文藝(2022年10期)2022-06-08 02:33:28
    彩泥大變身
    Improving robustness of complex networks by a new capacity allocation strategy
    The energy-saving advantages of burst-and-glide mode for thunniform swimming *
    Using spanwise flexibility of caudal fin to improve swimming performance for small fishlike robots *
    映像畜牧業(yè)
    搶電視
    愛你(2017年35期)2017-11-24 18:30:41
    搶電視
    圖說
    看圖作文寫作指導(dǎo)
    日本熟妇午夜| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清| 精品一区二区三区视频在线 | 九九久久精品国产亚洲av麻豆| 制服丝袜大香蕉在线| 亚洲人与动物交配视频| 看免费av毛片| av在线蜜桃| 美女黄网站色视频| 国模一区二区三区四区视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av二区三区四区| 亚洲真实伦在线观看| 女生性感内裤真人,穿戴方法视频| 欧美成狂野欧美在线观看| 级片在线观看| 日本黄色视频三级网站网址| 男女床上黄色一级片免费看| 久99久视频精品免费| 国产精品一区二区免费欧美| 亚洲av二区三区四区| 色噜噜av男人的天堂激情| 亚洲av一区综合| x7x7x7水蜜桃| 久久精品国产综合久久久| 好男人电影高清在线观看| 长腿黑丝高跟| 狂野欧美白嫩少妇大欣赏| 免费人成在线观看视频色| 亚洲中文字幕一区二区三区有码在线看| 国产精品久久久人人做人人爽| 久久久久久久精品吃奶| 老司机午夜福利在线观看视频| 欧美乱码精品一区二区三区| 波多野结衣高清无吗| 国产一区二区在线av高清观看| 欧美乱妇无乱码| 国产综合懂色| 人人妻人人看人人澡| 亚洲专区国产一区二区| 精品一区二区三区视频在线观看免费| 亚洲精品在线美女| 三级国产精品欧美在线观看| 好男人电影高清在线观看| 丁香欧美五月| 18禁国产床啪视频网站| 亚洲精品美女久久久久99蜜臀| 成人特级av手机在线观看| 久久久精品欧美日韩精品| 啪啪无遮挡十八禁网站| 久久午夜亚洲精品久久| 国产日本99.免费观看| 欧美中文综合在线视频| 国内精品久久久久精免费| 无限看片的www在线观看| av在线天堂中文字幕| 最后的刺客免费高清国语| 国内久久婷婷六月综合欲色啪| 婷婷亚洲欧美| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 在线播放国产精品三级| 精品一区二区三区人妻视频| 亚洲av电影在线进入| 18禁黄网站禁片午夜丰满| 制服人妻中文乱码| 91在线观看av| 男人舔女人下体高潮全视频| 美女cb高潮喷水在线观看| 国产又黄又爽又无遮挡在线| 国产一区二区三区视频了| 久久久久久久久久黄片| 99热6这里只有精品| 中文字幕高清在线视频| 此物有八面人人有两片| 女人十人毛片免费观看3o分钟| 精品久久久久久成人av| 日韩中文字幕欧美一区二区| 久久久久久久久久黄片| 国产中年淑女户外野战色| 在线国产一区二区在线| 亚洲精品成人久久久久久| 女人十人毛片免费观看3o分钟| 内地一区二区视频在线| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 日韩亚洲欧美综合| 可以在线观看毛片的网站| 亚洲 国产 在线| 嫩草影院精品99| 性欧美人与动物交配| 99国产精品一区二区三区| 久久精品国产自在天天线| 中亚洲国语对白在线视频| 亚洲天堂国产精品一区在线| av在线蜜桃| 欧美性猛交黑人性爽| 免费在线观看成人毛片| 在线十欧美十亚洲十日本专区| 午夜福利欧美成人| 欧美xxxx黑人xx丫x性爽| 国产又黄又爽又无遮挡在线| 国产精品影院久久| 日本 欧美在线| 制服人妻中文乱码| 国产精品影院久久| 欧美成人性av电影在线观看| 亚洲第一电影网av| 久久久久久九九精品二区国产| 综合色av麻豆| 首页视频小说图片口味搜索| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 桃色一区二区三区在线观看| 波野结衣二区三区在线 | 久久草成人影院| 99久久九九国产精品国产免费| 他把我摸到了高潮在线观看| 午夜精品久久久久久毛片777| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| 91麻豆精品激情在线观看国产| 日本a在线网址| 美女高潮的动态| 国产毛片a区久久久久| 国产免费av片在线观看野外av| 禁无遮挡网站| 欧美日韩国产亚洲二区| 午夜免费观看网址| 三级毛片av免费| 国产免费男女视频| 久久99热这里只有精品18| 青草久久国产| 在线观看舔阴道视频| 美女免费视频网站| 欧美成狂野欧美在线观看| 人人妻人人看人人澡| 一级a爱片免费观看的视频| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 男女下面进入的视频免费午夜| 国产成人影院久久av| 国产激情欧美一区二区| 村上凉子中文字幕在线| 桃红色精品国产亚洲av| 亚洲国产日韩欧美精品在线观看 | 国产aⅴ精品一区二区三区波| 三级国产精品欧美在线观看| 成人午夜高清在线视频| 啦啦啦韩国在线观看视频| 国产男靠女视频免费网站| 亚洲欧美日韩卡通动漫| 可以在线观看毛片的网站| 色在线成人网| 亚洲精华国产精华精| 国产色婷婷99| 免费看a级黄色片| 免费搜索国产男女视频| 人人妻人人看人人澡| 国产不卡一卡二| 亚洲最大成人手机在线| 一个人看视频在线观看www免费 | 国产探花极品一区二区| 国产私拍福利视频在线观看| 色视频www国产| 黄色日韩在线| 最近视频中文字幕2019在线8| 蜜桃久久精品国产亚洲av| 亚洲精品在线美女| 久久6这里有精品| 国产色爽女视频免费观看| 亚洲精品亚洲一区二区| 一边摸一边抽搐一进一小说| 熟妇人妻久久中文字幕3abv| 国产av不卡久久| 国产97色在线日韩免费| 99久久九九国产精品国产免费| 色在线成人网| or卡值多少钱| 热99re8久久精品国产| 国产久久久一区二区三区| 日韩欧美在线乱码| 日韩精品中文字幕看吧| 亚洲美女黄片视频| 一二三四社区在线视频社区8| 母亲3免费完整高清在线观看| 亚洲一区高清亚洲精品| 国产精品一区二区免费欧美| 免费av不卡在线播放| 欧美bdsm另类| 国内毛片毛片毛片毛片毛片| 黄色片一级片一级黄色片| 午夜精品一区二区三区免费看| 亚洲专区国产一区二区| 亚洲片人在线观看| 日韩av在线大香蕉| 免费在线观看日本一区| 日韩精品青青久久久久久| 老鸭窝网址在线观看| 免费一级毛片在线播放高清视频| 国产免费av片在线观看野外av| 国产成人av教育| 精品福利观看| 叶爱在线成人免费视频播放| 看免费av毛片| 香蕉av资源在线| 在线免费观看不下载黄p国产 | 一级毛片高清免费大全| 亚洲av美国av| 亚洲久久久久久中文字幕| 91久久精品电影网| 观看美女的网站| 国产高清三级在线| 久久这里只有精品中国| 国产亚洲欧美98| 亚洲人成网站在线播放欧美日韩| 中文字幕熟女人妻在线| 久久久久久久亚洲中文字幕 | 免费大片18禁| a级一级毛片免费在线观看| 午夜福利欧美成人| 亚洲成a人片在线一区二区| 精品不卡国产一区二区三区| 性色av乱码一区二区三区2| 亚洲av成人精品一区久久| 成人三级黄色视频| 成人性生交大片免费视频hd| 99视频精品全部免费 在线| 亚洲av一区综合| 夜夜躁狠狠躁天天躁| 午夜影院日韩av| 亚洲最大成人手机在线| 婷婷精品国产亚洲av| 国产不卡一卡二| 成人无遮挡网站| 人人妻人人看人人澡| 亚洲精品亚洲一区二区| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 国产av不卡久久| 法律面前人人平等表现在哪些方面| 午夜福利高清视频| 美女cb高潮喷水在线观看| 99国产精品一区二区三区| 国产av不卡久久| 久久伊人香网站| 日本与韩国留学比较| 午夜精品在线福利| 精品不卡国产一区二区三区| 母亲3免费完整高清在线观看| netflix在线观看网站| 黄片小视频在线播放| www国产在线视频色| av在线蜜桃| 母亲3免费完整高清在线观看| 91在线观看av| 亚洲无线在线观看| 免费观看精品视频网站| 国产精品av视频在线免费观看| 桃色一区二区三区在线观看| 亚洲精品亚洲一区二区| 成人av在线播放网站| 免费无遮挡裸体视频| 九九热线精品视视频播放| 九色国产91popny在线| 国产一区二区在线观看日韩 | 亚洲精品在线美女| 免费搜索国产男女视频| 国产熟女xx| 国产一区二区在线av高清观看| av黄色大香蕉| 男女做爰动态图高潮gif福利片| 亚洲人成伊人成综合网2020| 91字幕亚洲| 给我免费播放毛片高清在线观看| h日本视频在线播放| 免费大片18禁| 一级毛片高清免费大全| 有码 亚洲区| 国产 一区 欧美 日韩| 午夜福利欧美成人| 国产精品久久久久久久久免 | 热99在线观看视频| 人妻久久中文字幕网| 国产av在哪里看| 久久久国产成人免费| 国产黄片美女视频| 岛国在线免费视频观看| 成人欧美大片| 少妇的逼好多水| 黄色丝袜av网址大全| 丁香欧美五月| 中出人妻视频一区二区| 国产真实伦视频高清在线观看 | 亚洲五月婷婷丁香| 国产亚洲精品久久久久久毛片| 亚洲av第一区精品v没综合| 欧美高清成人免费视频www| 亚洲精品美女久久久久99蜜臀| 尤物成人国产欧美一区二区三区| 黄色女人牲交| 在线播放国产精品三级| 九九久久精品国产亚洲av麻豆| 蜜桃亚洲精品一区二区三区| 日韩欧美精品v在线| 美女 人体艺术 gogo| 国产午夜福利久久久久久| 中亚洲国语对白在线视频| 成人高潮视频无遮挡免费网站| 黄色视频,在线免费观看| 99热精品在线国产| 国产精品一区二区三区四区久久| 制服人妻中文乱码| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 少妇人妻一区二区三区视频| 日本黄色片子视频| 最好的美女福利视频网| av片东京热男人的天堂| 日本与韩国留学比较| 美女 人体艺术 gogo| 日韩欧美 国产精品| a在线观看视频网站| 老熟妇乱子伦视频在线观看| а√天堂www在线а√下载| 99久久99久久久精品蜜桃| 久久久久久国产a免费观看| 精品久久久久久久久久免费视频| 美女高潮喷水抽搐中文字幕| 国产成人a区在线观看| 宅男免费午夜| 国产在线精品亚洲第一网站| 国产高清视频在线观看网站| 高潮久久久久久久久久久不卡| 美女黄网站色视频| 久久久久久久亚洲中文字幕 | 性色av乱码一区二区三区2| 国产av一区在线观看免费| 日韩精品青青久久久久久| 18禁在线播放成人免费| 国产精品久久久久久久久免 | 国产精品自产拍在线观看55亚洲| svipshipincom国产片| 亚洲一区二区三区不卡视频| 男女床上黄色一级片免费看| 成熟少妇高潮喷水视频| 日韩欧美 国产精品| 日韩av在线大香蕉| 国产伦精品一区二区三区四那| www日本黄色视频网| 亚洲国产色片| 在线十欧美十亚洲十日本专区| 99在线人妻在线中文字幕| 又爽又黄无遮挡网站| 欧美在线一区亚洲| 午夜精品久久久久久毛片777| 中文字幕高清在线视频| 久久99热这里只有精品18| 操出白浆在线播放| 99久久九九国产精品国产免费| 国产国拍精品亚洲av在线观看 | 在线观看舔阴道视频| 欧美午夜高清在线| 无人区码免费观看不卡| 黄色成人免费大全| 亚洲av电影在线进入| 高清在线国产一区| 精品电影一区二区在线| 欧美在线一区亚洲| 免费在线观看日本一区| 村上凉子中文字幕在线| 男人的好看免费观看在线视频| 中文资源天堂在线| 午夜免费观看网址| 夜夜夜夜夜久久久久| 日本成人三级电影网站| 丰满乱子伦码专区| 久久伊人香网站| 午夜影院日韩av| 久99久视频精品免费| 国产精品99久久99久久久不卡| 久久精品综合一区二区三区| 99久久综合精品五月天人人| 久久精品91蜜桃| 亚洲欧美一区二区三区黑人| 欧美一级a爱片免费观看看| 一级a爱片免费观看的视频| 久久久色成人| 搞女人的毛片| 免费高清视频大片| 国产97色在线日韩免费| 欧美三级亚洲精品| 久久久久久久精品吃奶| 99热这里只有是精品50| 人人妻人人澡欧美一区二区| 桃红色精品国产亚洲av| 日本熟妇午夜| 亚洲精品粉嫩美女一区| 欧美国产日韩亚洲一区| 97碰自拍视频| 黄片大片在线免费观看| 亚洲一区二区三区不卡视频| 亚洲av免费在线观看| 中出人妻视频一区二区| 伊人久久大香线蕉亚洲五| 中国美女看黄片| 国产亚洲精品综合一区在线观看| 他把我摸到了高潮在线观看| 欧美极品一区二区三区四区| 日韩欧美精品v在线| 51国产日韩欧美| 国产亚洲精品久久久久久毛片| 国内毛片毛片毛片毛片毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产又黄又爽又无遮挡在线| 18禁黄网站禁片免费观看直播| 操出白浆在线播放| 欧美区成人在线视频| 国产精品女同一区二区软件 | 夜夜爽天天搞| 欧美乱码精品一区二区三区| 亚洲熟妇中文字幕五十中出| 午夜福利成人在线免费观看| 精品电影一区二区在线| 亚洲精品久久国产高清桃花| 久久久久免费精品人妻一区二区| 99热只有精品国产| 亚洲第一欧美日韩一区二区三区| 91久久精品电影网| 国模一区二区三区四区视频| 国产成+人综合+亚洲专区| 国产黄a三级三级三级人| 精品人妻一区二区三区麻豆 | 国产高清激情床上av| 99热6这里只有精品| 身体一侧抽搐| 天天躁日日操中文字幕| 亚洲国产欧美网| 天堂√8在线中文| 国产高清videossex| 3wmmmm亚洲av在线观看| 欧美日韩福利视频一区二区| 欧美三级亚洲精品| 亚洲五月天丁香| 国产精品乱码一区二三区的特点| 欧美黑人欧美精品刺激| 国内久久婷婷六月综合欲色啪| 日本五十路高清| 亚洲人成电影免费在线| 国产视频一区二区在线看| 亚洲国产精品成人综合色| 亚洲午夜理论影院| АⅤ资源中文在线天堂| eeuss影院久久| 女人十人毛片免费观看3o分钟| 免费一级毛片在线播放高清视频| 日韩欧美一区二区三区在线观看| 国产97色在线日韩免费| 少妇高潮的动态图| 午夜亚洲福利在线播放| netflix在线观看网站| 久久香蕉国产精品| 亚洲人成网站高清观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美3d第一页| 色哟哟哟哟哟哟| 国产成人影院久久av| 黄色日韩在线| 在线观看美女被高潮喷水网站 | 色视频www国产| 午夜福利在线在线| 欧美色视频一区免费| 夜夜看夜夜爽夜夜摸| 欧美xxxx黑人xx丫x性爽| 少妇高潮的动态图| 日本a在线网址| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 免费观看精品视频网站| 99精品欧美一区二区三区四区| 九色国产91popny在线| 18禁黄网站禁片免费观看直播| 国产日本99.免费观看| 午夜福利在线在线| 可以在线观看的亚洲视频| 亚洲成人久久性| 91久久精品国产一区二区成人 | 99久久精品热视频| 欧美中文综合在线视频| 最新中文字幕久久久久| 精品99又大又爽又粗少妇毛片 | 观看免费一级毛片| 欧美成狂野欧美在线观看| 熟女电影av网| 男女视频在线观看网站免费| 欧美+日韩+精品| 欧美性感艳星| netflix在线观看网站| 天天添夜夜摸| 18+在线观看网站| 听说在线观看完整版免费高清| netflix在线观看网站| 黄色视频,在线免费观看| 小蜜桃在线观看免费完整版高清| 床上黄色一级片| 淫妇啪啪啪对白视频| 国产不卡一卡二| 欧美3d第一页| 无遮挡黄片免费观看| 国产单亲对白刺激| 日韩成人在线观看一区二区三区| 国产单亲对白刺激| 日韩欧美精品v在线| 国产成人欧美在线观看| 欧美av亚洲av综合av国产av| 免费看日本二区| 欧美在线黄色| 亚洲avbb在线观看| 亚洲美女视频黄频| 特大巨黑吊av在线直播| 国产真人三级小视频在线观看| 91麻豆av在线| 3wmmmm亚洲av在线观看| 午夜福利成人在线免费观看| 久久精品国产亚洲av香蕉五月| 国产在视频线在精品| 97超视频在线观看视频| 国产色爽女视频免费观看| 波多野结衣高清作品| 国产精品野战在线观看| 久久久国产精品麻豆| 欧美乱码精品一区二区三区| 极品教师在线免费播放| 久久精品亚洲精品国产色婷小说| 在线观看美女被高潮喷水网站 | 亚洲久久久久久中文字幕| 男女做爰动态图高潮gif福利片| 又爽又黄无遮挡网站| h日本视频在线播放| 岛国在线免费视频观看| 老司机深夜福利视频在线观看| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| 欧美黄色片欧美黄色片| 国产色婷婷99| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 极品教师在线免费播放| 免费看光身美女| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 国模一区二区三区四区视频| 亚洲 欧美 日韩 在线 免费| 亚洲av免费高清在线观看| 国产精品乱码一区二三区的特点| 丝袜美腿在线中文| 麻豆成人av在线观看| 亚洲精品影视一区二区三区av| 激情在线观看视频在线高清| 母亲3免费完整高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| www.www免费av| 国产午夜福利久久久久久| 精品一区二区三区视频在线 | 国产精品美女特级片免费视频播放器| 中文字幕av成人在线电影| 很黄的视频免费| 日韩欧美精品v在线| 亚洲人成网站在线播| 啦啦啦免费观看视频1| 欧美国产日韩亚洲一区| 国产欧美日韩一区二区三| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久| 国产精品电影一区二区三区| 国产一区二区激情短视频| 久久久久久久久中文| 我要搜黄色片| h日本视频在线播放| 给我免费播放毛片高清在线观看| 色哟哟哟哟哟哟| 免费在线观看亚洲国产| 午夜福利18| 99热这里只有精品一区| 91麻豆精品激情在线观看国产| 九色国产91popny在线| 亚洲精品粉嫩美女一区| 国产在线精品亚洲第一网站| 国产一区二区激情短视频| 搡老熟女国产l中国老女人| 观看美女的网站| 网址你懂的国产日韩在线| 给我免费播放毛片高清在线观看| 老熟妇仑乱视频hdxx| xxxwww97欧美| 欧美av亚洲av综合av国产av| 日本三级黄在线观看| 黄片小视频在线播放| 麻豆一二三区av精品| 欧美一区二区亚洲| 一个人看的www免费观看视频| 午夜精品在线福利| 久久精品人妻少妇| 性欧美人与动物交配| 国产精品久久久久久精品电影| 高潮久久久久久久久久久不卡| 精品久久久久久久毛片微露脸| 欧美+日韩+精品| 欧美xxxx黑人xx丫x性爽| 两个人的视频大全免费| 国产日本99.免费观看| 日本 av在线| 免费在线观看日本一区| 在线观看美女被高潮喷水网站 |