• <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蜜桃| 搡女人真爽免费视频火全软件| 国产日韩欧美亚洲二区| 99热这里只有是精品在线观看| 国产免费福利视频在线观看| 寂寞人妻少妇视频99o| 美女被艹到高潮喷水动态| 亚洲三级黄色毛片| 婷婷色综合www| 亚洲精品影视一区二区三区av| 丰满乱子伦码专区| 哪个播放器可以免费观看大片| 久久99热6这里只有精品| 久久久色成人| 久久午夜福利片| 免费大片18禁| 高清日韩中文字幕在线| 欧美日韩视频精品一区| 夜夜看夜夜爽夜夜摸| 一个人看的www免费观看视频| 热re99久久精品国产66热6| 最近中文字幕高清免费大全6| 搡老乐熟女国产| 亚洲美女视频黄频| 免费观看a级毛片全部| 永久网站在线| 亚洲自拍偷在线| 国产精品一二三区在线看| av国产久精品久网站免费入址| 欧美日韩在线观看h| 亚洲欧美精品自产自拍| 日韩,欧美,国产一区二区三区| 国产免费又黄又爽又色| 极品教师在线视频| 日本欧美国产在线视频| 色婷婷久久久亚洲欧美| 在现免费观看毛片| 热99国产精品久久久久久7| 亚洲欧美日韩卡通动漫| 18禁裸乳无遮挡动漫免费视频 | kizo精华| 精品一区二区三区视频在线| 在线观看免费高清a一片| 嫩草影院精品99| 夫妻午夜视频| 亚洲自偷自拍三级| av免费观看日本| 欧美日韩在线观看h| 三级国产精品欧美在线观看| 日韩一本色道免费dvd| 一个人看的www免费观看视频| 大陆偷拍与自拍| 日韩在线高清观看一区二区三区| 久久这里有精品视频免费| 高清毛片免费看| 亚洲久久久久久中文字幕| 色婷婷久久久亚洲欧美| 在线观看免费高清a一片| 亚洲国产精品成人综合色| 嫩草影院精品99| 亚洲av福利一区| 国产精品av视频在线免费观看| 精华霜和精华液先用哪个| 久久国产乱子免费精品| 婷婷色麻豆天堂久久| 在线天堂最新版资源| 国产毛片a区久久久久| 18+在线观看网站| 欧美人与善性xxx| 国产 一区 欧美 日韩| 蜜桃久久精品国产亚洲av| 联通29元200g的流量卡| 日韩制服骚丝袜av| 久久99精品国语久久久| 欧美日韩一区二区视频在线观看视频在线 | av在线天堂中文字幕| 亚洲精品影视一区二区三区av| 国产精品精品国产色婷婷| 日韩av免费高清视频| 亚洲国产av新网站| 我的女老师完整版在线观看| 日韩中字成人| 亚洲欧美成人综合另类久久久| 久久久久久久久大av| 精品久久久久久久久亚洲| 午夜福利在线在线| 国产亚洲5aaaaa淫片| 成年女人看的毛片在线观看| 亚洲av欧美aⅴ国产| 亚洲精品自拍成人| 午夜福利在线在线| 成人国产麻豆网| 亚洲熟女精品中文字幕| 日韩av在线免费看完整版不卡| 1000部很黄的大片| 久久精品人妻少妇| 纵有疾风起免费观看全集完整版| 免费高清在线观看视频在线观看| 在线看a的网站| 在线天堂最新版资源| 成人亚洲欧美一区二区av| 国内少妇人妻偷人精品xxx网站| 免费观看av网站的网址| 久久久久久久久久人人人人人人| 欧美区成人在线视频| 久久久久久九九精品二区国产| 麻豆久久精品国产亚洲av| 欧美精品人与动牲交sv欧美| 一级毛片黄色毛片免费观看视频| 欧美精品一区二区大全| 亚洲精品视频女| 久久久亚洲精品成人影院| 永久网站在线| 亚洲精品久久久久久婷婷小说| 国产精品一及| 搡女人真爽免费视频火全软件| 久久久久久久久大av| 日韩大片免费观看网站| 国产成人免费观看mmmm| 国产一区亚洲一区在线观看| 国产精品国产三级国产av玫瑰| av国产精品久久久久影院| 国产美女午夜福利| 亚洲精品成人av观看孕妇| 国产亚洲一区二区精品| 免费人成在线观看视频色| 免费少妇av软件| 日本一二三区视频观看| 亚洲国产精品成人久久小说| 美女高潮的动态| 免费在线观看成人毛片| 亚洲av在线观看美女高潮| 精品久久久久久久人妻蜜臀av| 波野结衣二区三区在线| 久久久亚洲精品成人影院| 美女国产视频在线观看| 久久人人爽人人爽人人片va| 亚洲aⅴ乱码一区二区在线播放| 最后的刺客免费高清国语| 一本久久精品| 人妻系列 视频| 男人添女人高潮全过程视频| 亚洲精华国产精华液的使用体验| 岛国毛片在线播放| 国产爱豆传媒在线观看| 日韩精品有码人妻一区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 26uuu在线亚洲综合色| 久久久久久久久大av| 自拍欧美九色日韩亚洲蝌蚪91 | 99re6热这里在线精品视频| 免费大片黄手机在线观看| 日本免费在线观看一区| 1000部很黄的大片| 只有这里有精品99| 亚洲第一区二区三区不卡| 大陆偷拍与自拍| 日日啪夜夜撸| 亚洲av成人精品一二三区| 一本色道久久久久久精品综合| 久久久久久国产a免费观看| 嫩草影院入口| 亚洲人与动物交配视频| av黄色大香蕉| 日韩电影二区| 内射极品少妇av片p| 免费观看av网站的网址| 2022亚洲国产成人精品| 免费av毛片视频| 一级片'在线观看视频| av国产精品久久久久影院| a级毛色黄片| 亚洲欧美成人精品一区二区| 人体艺术视频欧美日本| 国产 一区精品| 国产成人精品婷婷| 午夜福利视频精品| 久久久久国产精品人妻一区二区| 亚洲人成网站高清观看| 久久久久久久久久人人人人人人| 交换朋友夫妻互换小说| 免费看a级黄色片| eeuss影院久久| 日韩av免费高清视频| 69av精品久久久久久| 在线免费十八禁| 极品教师在线视频| av专区在线播放| 国产淫语在线视频| 国产黄片美女视频| 亚洲无线观看免费| 欧美xxⅹ黑人| av福利片在线观看| 久久久久久九九精品二区国产| 亚洲av成人精品一二三区| 欧美性感艳星| 国产人妻一区二区三区在| 亚洲精品一区蜜桃| 秋霞伦理黄片| 岛国毛片在线播放| 亚洲精品日韩在线中文字幕| 又粗又硬又长又爽又黄的视频| 国模一区二区三区四区视频| 日本黄色片子视频| 欧美一级a爱片免费观看看| 精品国产一区二区三区久久久樱花 | 一区二区三区四区激情视频| 好男人在线观看高清免费视频| 日韩 亚洲 欧美在线| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 精品人妻视频免费看| 国产精品人妻久久久影院| 欧美日韩视频高清一区二区三区二| 亚洲精品第二区| 老女人水多毛片| 免费看光身美女| 午夜福利网站1000一区二区三区| 99热国产这里只有精品6| 欧美精品国产亚洲| 久久午夜福利片| 国产视频内射| 精品久久久久久久人妻蜜臀av| 国产av码专区亚洲av| 成人国产麻豆网| 国产综合精华液| 一本一本综合久久| 久久久久久久亚洲中文字幕| 欧美精品一区二区大全| 一级a做视频免费观看| 简卡轻食公司| videos熟女内射| 人妻一区二区av| 国内精品美女久久久久久| 亚洲国产av新网站| 91在线精品国自产拍蜜月| 在线免费十八禁| 真实男女啪啪啪动态图| 内射极品少妇av片p| 精品酒店卫生间| 久久久亚洲精品成人影院| 国内揄拍国产精品人妻在线| 久久精品久久精品一区二区三区| 性色av一级| 天堂网av新在线| 白带黄色成豆腐渣| 亚洲欧美精品专区久久| 久久久久久久午夜电影| 日韩欧美精品免费久久| 99热这里只有精品一区| 97超碰精品成人国产| 美女脱内裤让男人舔精品视频| av在线观看视频网站免费| 色哟哟·www| 久久久色成人| 成人欧美大片| 我的老师免费观看完整版| 国产亚洲91精品色在线| 丰满少妇做爰视频| 又爽又黄无遮挡网站| 精品久久久久久久人妻蜜臀av| 波野结衣二区三区在线| 国产精品秋霞免费鲁丝片| 亚洲,一卡二卡三卡| a级毛色黄片| 欧美日本视频| 国产成人精品久久久久久| 久久久久久久久大av| 少妇猛男粗大的猛烈进出视频 | 一个人看视频在线观看www免费| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影| 少妇裸体淫交视频免费看高清| 国产人妻一区二区三区在| 久久久久久久亚洲中文字幕| 国产黄片美女视频| 亚洲精品乱码久久久久久按摩| 国内揄拍国产精品人妻在线| 噜噜噜噜噜久久久久久91| 丰满少妇做爰视频| 欧美丝袜亚洲另类| 国产午夜精品一二区理论片| 在线a可以看的网站| 街头女战士在线观看网站| 亚洲最大成人av| 国内精品宾馆在线| 国产乱人偷精品视频| 网址你懂的国产日韩在线| 国产欧美另类精品又又久久亚洲欧美| 狂野欧美激情性xxxx在线观看| 在线免费观看不下载黄p国产| 夜夜爽夜夜爽视频| 人体艺术视频欧美日本| 少妇人妻久久综合中文| 国产免费视频播放在线视频| 女人被狂操c到高潮| 天天一区二区日本电影三级| 国产精品一及| 精品久久久精品久久久| 成年av动漫网址| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 久久精品综合一区二区三区| 在线 av 中文字幕| 欧美一级a爱片免费观看看| 男的添女的下面高潮视频| 色哟哟·www| 国产成人aa在线观看| 久久久欧美国产精品| 中文字幕av成人在线电影| 亚洲av一区综合| 国产 精品1| 亚洲人成网站在线观看播放| 婷婷色综合大香蕉| 69av精品久久久久久| 国产 一区 欧美 日韩| 肉色欧美久久久久久久蜜桃 | 国产亚洲av嫩草精品影院| 国产乱来视频区| 亚洲国产日韩一区二区| 3wmmmm亚洲av在线观看| 亚洲aⅴ乱码一区二区在线播放| 一级av片app| 久久久久久久久久久免费av| 91aial.com中文字幕在线观看| 听说在线观看完整版免费高清| 韩国高清视频一区二区三区| 人妻一区二区av| 日本与韩国留学比较| 日韩强制内射视频| 超碰97精品在线观看| 国产毛片a区久久久久| 日日撸夜夜添| 免费观看的影片在线观看| 亚洲精品一区蜜桃| h日本视频在线播放| 在线观看一区二区三区| 欧美一区二区亚洲| 亚洲av在线观看美女高潮| 日韩成人伦理影院| 国产真实伦视频高清在线观看| 观看美女的网站| 亚洲精品成人av观看孕妇| 可以在线观看毛片的网站| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 欧美丝袜亚洲另类| 国产精品99久久久久久久久| 免费看a级黄色片| 91久久精品国产一区二区成人| 欧美 日韩 精品 国产| 亚洲在久久综合| 少妇人妻 视频| 欧美成人a在线观看| videos熟女内射| 久久久亚洲精品成人影院| 一级毛片久久久久久久久女| 亚洲欧美一区二区三区国产| 欧美一级a爱片免费观看看| 在线观看av片永久免费下载| 色哟哟·www| 下体分泌物呈黄色| 欧美激情国产日韩精品一区| 国产精品秋霞免费鲁丝片| 亚洲成人久久爱视频| 男人添女人高潮全过程视频| 婷婷色综合大香蕉| av在线播放精品| av在线app专区| 两个人的视频大全免费| 黄色怎么调成土黄色| 国产免费又黄又爽又色| 国产免费一级a男人的天堂| 一级毛片久久久久久久久女| 国产成人a区在线观看| 免费高清在线观看视频在线观看| 久久精品综合一区二区三区| 成人鲁丝片一二三区免费| 国产毛片在线视频| 亚洲成人中文字幕在线播放| 久久午夜福利片| 18+在线观看网站| 美女主播在线视频| 天天躁日日操中文字幕| 在线观看人妻少妇| 在线观看av片永久免费下载| 久久影院123| 青春草国产在线视频| 欧美+日韩+精品| 美女主播在线视频| 久久国内精品自在自线图片| 人妻 亚洲 视频| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| 国产乱人偷精品视频| 午夜福利网站1000一区二区三区| 夫妻午夜视频| 久久99热6这里只有精品| 伦理电影大哥的女人| 搡女人真爽免费视频火全软件| 嫩草影院入口| 亚洲精品一区蜜桃| 国内精品美女久久久久久| a级毛片免费高清观看在线播放| 晚上一个人看的免费电影| 国产男女超爽视频在线观看| 国产男女内射视频| 国产69精品久久久久777片| 国产黄片美女视频| 久久久久久久午夜电影| 婷婷色av中文字幕| 99视频精品全部免费 在线| 亚洲精品国产色婷婷电影| 18禁动态无遮挡网站| 一级毛片黄色毛片免费观看视频| 亚洲国产欧美人成| 国产美女午夜福利| 人妻制服诱惑在线中文字幕| 丰满人妻一区二区三区视频av| 国产 精品1| 伦理电影大哥的女人| 国产国拍精品亚洲av在线观看| 国产精品久久久久久精品古装| 国产精品.久久久| 亚洲av福利一区| 国产一级毛片在线| 午夜福利在线在线| 青春草亚洲视频在线观看| av在线老鸭窝| 日韩欧美精品v在线| 日韩免费高清中文字幕av| 人妻 亚洲 视频| 免费播放大片免费观看视频在线观看| 亚洲精品日本国产第一区| 精品少妇黑人巨大在线播放| 午夜激情久久久久久久| 亚洲国产成人一精品久久久| 两个人的视频大全免费| 男男h啪啪无遮挡| 成人一区二区视频在线观看| 看非洲黑人一级黄片| 好男人在线观看高清免费视频| .国产精品久久| 国产视频首页在线观看| 亚洲av在线观看美女高潮| 亚洲熟女精品中文字幕| 久久久色成人| av.在线天堂| 在线观看免费高清a一片| 国产伦精品一区二区三区四那| 亚洲av.av天堂| a级毛色黄片| 欧美xxⅹ黑人| 国产精品爽爽va在线观看网站| av在线老鸭窝| 又大又黄又爽视频免费| 丝袜脚勾引网站| 99视频精品全部免费 在线| 国产综合精华液| 久久久欧美国产精品| 国产精品国产三级国产av玫瑰| 又爽又黄无遮挡网站| 欧美精品一区二区大全| 校园人妻丝袜中文字幕| 联通29元200g的流量卡| 亚洲国产av新网站| 三级国产精品片| 日本黄大片高清| 69av精品久久久久久| 赤兔流量卡办理| 我的女老师完整版在线观看| 欧美日韩亚洲高清精品| 国产精品久久久久久久电影| 国产爱豆传媒在线观看| 插阴视频在线观看视频| 亚洲色图av天堂| 精品一区在线观看国产| av黄色大香蕉| 亚洲天堂av无毛| 国产美女午夜福利| 欧美日韩一区二区视频在线观看视频在线 | 成人毛片60女人毛片免费| 啦啦啦中文免费视频观看日本| 欧美日本视频| 黑人高潮一二区| 九九爱精品视频在线观看| 亚洲真实伦在线观看| 日韩亚洲欧美综合| 国产精品一区二区在线观看99| 久久久久久久精品精品| 人人妻人人看人人澡| 国产国拍精品亚洲av在线观看| 免费少妇av软件| 一级毛片 在线播放| 久久久久久久久久久丰满| 国产精品国产三级专区第一集| 精品人妻熟女av久视频| 亚洲av二区三区四区| 国产精品精品国产色婷婷| 午夜视频国产福利| 有码 亚洲区| 亚洲在线观看片| 日韩大片免费观看网站| 国产久久久一区二区三区| 国产精品偷伦视频观看了| 51国产日韩欧美| 纵有疾风起免费观看全集完整版| 最后的刺客免费高清国语| 久久久久久久久久久免费av| 天美传媒精品一区二区| 99久久精品一区二区三区| 国产成年人精品一区二区| 97在线人人人人妻| 91aial.com中文字幕在线观看| 少妇人妻精品综合一区二区| 亚洲av日韩在线播放| 一本久久精品| 国产成人freesex在线| av在线天堂中文字幕| 久久亚洲国产成人精品v| 国产男女内射视频| 久久这里有精品视频免费| 有码 亚洲区| 精品视频人人做人人爽| av在线观看视频网站免费| 精品人妻熟女av久视频| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 精品酒店卫生间| 一二三四中文在线观看免费高清| 高清视频免费观看一区二区| 一级二级三级毛片免费看| 夜夜爽夜夜爽视频| 国产精品人妻久久久影院| 成人美女网站在线观看视频| 七月丁香在线播放| 久久精品国产鲁丝片午夜精品| 成人国产av品久久久| 一级毛片黄色毛片免费观看视频| 精品久久久精品久久久| 波野结衣二区三区在线| 嫩草影院入口| 网址你懂的国产日韩在线| 久久97久久精品| 色综合色国产| av在线播放精品| www.av在线官网国产| 国产精品女同一区二区软件| 亚洲无线观看免费| 黄色怎么调成土黄色| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| 少妇高潮的动态图| 91久久精品国产一区二区成人| 国产精品一区二区三区四区免费观看| 久久6这里有精品| 看免费成人av毛片| 亚洲av一区综合| 十八禁网站网址无遮挡 | 国产男女超爽视频在线观看| 91精品伊人久久大香线蕉| 国产黄片美女视频| a级一级毛片免费在线观看| 亚洲av免费在线观看| 国产 一区 欧美 日韩| av又黄又爽大尺度在线免费看| 国产成人精品一,二区| 深夜a级毛片| 亚洲精品456在线播放app| 日韩中字成人| 日日啪夜夜爽| 高清日韩中文字幕在线| videos熟女内射| 亚洲av免费在线观看| 青春草视频在线免费观看| 欧美精品国产亚洲| 亚洲丝袜综合中文字幕| 日韩强制内射视频| 日韩国内少妇激情av| 午夜福利视频精品| 亚洲国产高清在线一区二区三| 国产淫片久久久久久久久| 亚洲精品自拍成人| 成年免费大片在线观看| 97人妻精品一区二区三区麻豆| 国产欧美日韩一区二区三区在线 | 精品午夜福利在线看| 中文字幕久久专区| 91久久精品国产一区二区三区| 国产成人一区二区在线| av国产精品久久久久影院| 一个人看的www免费观看视频| 嫩草影院精品99| 国产欧美日韩精品一区二区| 久久ye,这里只有精品| 免费看不卡的av| 国产亚洲91精品色在线| 国产色爽女视频免费观看| 国产成人福利小说| 亚洲国产精品999| 久久6这里有精品| 99视频精品全部免费 在线| 亚洲精品日本国产第一区| 午夜日本视频在线| 亚洲精品456在线播放app| 亚洲怡红院男人天堂| 美女视频免费永久观看网站| 亚洲国产av新网站| 成人美女网站在线观看视频| av一本久久久久| 日韩,欧美,国产一区二区三区| 97热精品久久久久久|