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

    Investigation of the hydrodynamic performance of crablike robot swimming leg *

    2018-09-28 05:34:04LiquanWang王立權(quán)HailongWang王海龍GangWang王剛XiChen陳曦AskerKhanLixingJin靳勵(lì)行
    關(guān)鍵詞:陳曦王剛海龍

    Li-quan Wang (王立權(quán)), Hai-long Wang (王海龍), Gang Wang (王剛), Xi Chen (陳曦), Asker Khan ,Li-xing Jin (靳勵(lì)行)

    1. College of Mechanical and Electrical Engineering, Harbin Engineering University, Harbin 150001, China

    2. State Key Laboratory of Autonomous Underwater Vehicle, Harbin Engineering University, Harbin 150001,China

    3. College of Mechanical and Electrical Engineering, Heilongjiang Institute of Technology, Harbin 150001,China

    Abstract: The existing amphibious robots cannot usually enjoy a superior adaptability in the underwater environment by replacing the actuators. Based on the bionic prototype of the Portunus trituberculatus, a new leg-paddle coupling crablike robot with a composite propulsion of walking legs and swimming legs is developed, with both the abilities of walking and swimming under water.By simulation and experiment, the effects of the phase difference, the flapping amplitude and the angular bias of the coupling movement, as well as the Strouhal number on the hydrodynamic performance of the swimming legs are studied, and the time dependent tail vortex shedding structure in a cycle is obtained. Both experimental and numerical results indicate that the thrust force with a high propulsion efficiency can be generated by a flapping swimming leg. This work can further be used for analysis of the stability and the maneuverability of the swimming leg actuated underwater vehicles.

    Key words: Leg-paddle coupling crablike robot, swimming leg, hydrofoil propulsion, hydrodynamic performance

    Introduction

    Various kinds of unmanned underwater vehicles were developed for the marine industry. The traditional underwater vehicles equipped with propulsive propellers show low propulsive efficiency and poor maneuverability[1]. The bionic propulsion enjoys many advantages, including low noise and lesser resistance, good concealing ability, strong adaptability,low cost and mass production, therefore, has broad application prospects in the field of underwater vehicles[2]. According to their propulsion methods of imitating underwater creatures, the propulsions of the existing bionic underwater vehicles can be classified as: the bionic fish oscillating propulsion, the bionic squid jet propulsion, the bionic multi-legged walking and the bionic turtle hydrofoil propulsion[3-6].

    The Portunus trituberculatus is a kind of crab living in shoals. It has the abilities of both walking on the seafloor by its 3 pairs of walking legs and swimming in water by its last pair of swimming legs.It has a good adaptability to the sea wave and current in the underwater environment[7-8]. Based on the bionic prototype of the Portunus trituberculatus, a new leg-paddle coupling crablike robot is developed,which can both walk and swim in water. According to various working environments and tasks, it can switch the motion mode independently. There are few studies of swimming leg propulsion of amphibious crustaceans, such as crabs and aquatic insects, but quite many studies of fish pectoral-fin and sea turtle flexible forelimb-hydrofoil propulsion. The propulsion principles and numerical calculation methods of the latter cases can be applied in the research of swimming legs of the crablike robot. The unsteady two-dimensional blade element theory was used to calculate the average thrust and the mechanical efficiency of the fish pectoral fin by using the drag-based swimming mode as well as the lift-based swimming mode. Kato calculated and analyzed the hydrodynamic performance of 2DOF rigid pectoral fins which imitated the black sea bass by using the unsteady vortex lattice method. The effects of the phase difference between the shaking wings and the flapping wings on the average thrust and efficiency were compared with experimental results. The hydrodynamic performances of 3DOF rigid pectoral fins were investigated by numerically solving the three-dimensional unsteady overlapping grid Navier-Stokes equations, and the load characteristics of a mechanical pectoral fin were examined. The results show that the lift-based swimming mode is more suitable for the generation of the propulsive load than the drag-based swimming in a steady flow, and the opposite is the case in a still water[9]. The swimming mechanism of the Labriform fish was simulated. Patar Ebenezer Sitorus et al.[10]designed a fish robot with paired pectoral fins as the thruster and experimentally demonstrated the effect of the flapping frequency, the amplitude and the Strouhal number on the swimming speed. The results indicate that a higher flapping frequency does not necessarily produce a higher swimming speed. Lauder and Lauder.[11]applied the digital particle image velocimetry (DPIV) to study the fish pectoral fin oscillation, obtained the trailing vortex structure and the instantaneous velocity around the pectoral fin and analyzed the turning mechanism of the fish. The hydrodynamics of a 2DOF rigid pectoral fin was calculated by Chiu et al.[12]based on the blade element theory, which shows that the maneuverability of the fish robot can be improved through adjusting the phase difference between the shaking wings and the flapping wings. Gao et al.[13]developed a fish robot propelled by a pair of pectoral fins and carried out experiments in a swimming pool to study its hydrodynamic performance. He found that both the maximum average thrust coefficient and efficiency occur at the same Strouhal number of 0.4, in accordance with the St range of the propulsion of swimming and flying animals. Su et al.[14]studied the hydrodynamics of 2DOF and 3DOF rigid pectoral fins numerically and experimentally and analyzed the influence of the motion parameters on the hydrodynamics. The variations of the hydrodynamic coefficients can be strongly influenced by the rowing biased angle,the feathering biased angle and the phase difference of the coupling motion.

    In this paper, the propulsion mechanism of the Portunus trituberculatus is investigated. Then based on the bionic prototype of the Portunus trituberculatus, a leg-paddle coupling crablike robot, capable of operating in shoals, is established. Finally, numerical calculations and experiments are carried out to analyze the hydrodynamic performance of the swimming leg flapping motion and the complex flow field structure around the hydrofoil of the swimming leg.

    1. Propulsion mechanism of Portunus trituberculatus

    Crabs are crustaceans decapod of the infraorder Brachyura, which typically have a very short projecting tail, usually entirely hidden under the thorax[7]. They live in oceans, in fresh water and on land. They are generally classified into the freshwater crab and the sea crab. The freshwater crabs usually do not have the swimming ability, while the Portunidae is a family of sea crabs and some of them are swimming crabs. Its body is covered with a thick exoskeleton,divided into the cephalothorax and the abdomen, as shown in Fig. 1.

    Fig. 1 (Color online) Shape structure of Portunus trituberculatus

    The crab appendage consists typically of a shaft and 5 pereiopods. The Portunidae is characterized by the flattening of the last pair of legs to become swimming broad paddles, called the swimming legs.This ability, together with their strong chelipeds,allows many species to be fast and aggressive predators.

    Fig. 2 (Color online) 3-D model of crab-like robot

    The experimental study of the Portunidae by Vivo shows that the physiological structure of the swimming leg is similar to that of a walking leg. With the polyarticular skeleton texture, the swimming leg can achieve the flexible flapping under water. The typical swimming legs of the Portunidae consist of six components (the Coxa, the Basis, the Merus, the Carpus, the Propus, the Dacty1) articulated with each other by bicondylar joints[15]. The range of each segment motion is restricted to a single plane. Most of these planes are so positioned as perpendicular to the adjacent joints. 10 groups of stretched muscles are housed in the segments to control the movement of the swimming legs. The pleural muscles consist of 4 types,i.e., those of the Promotor, the Remotor, the Levator and the Depressor. The extensor and flexor musculature of the promotor, originated from the pleural venter and dorsum, generates the up-down motion of the Coxa around the TC joint. The mediolateral motion of the Basis around the CB joint is realized by the Levator and the Depressor. The flexer and extensor muscles in the Merus operate the MC joint to move the legs in the left-right direction. The stretcher and bender housed in the Carpus are responsible for rotating the Propus forward and backward. The opener and closer housed in the Propus operate the PD joint to move the Dacty1 in the up-down direction.Compared with the walking leg, the Propus and the Dacty1 of the swimming legs are muscular. Usually the Dacty1 is unfolded to achieve a greater thrust while the swimming leg is flapping. Both the ratio of the numbers of the nerve fibers and the ratio of the areas of the muscle fibers in the Propus and the Dacty1 are high.

    2. Bionic structure design

    2.1 Crablike robot prototype

    Based on the bionic prototype of the Portunus trituberculatus, a new leg-paddle coupling crablike robot with the composite propulsion of walking legs and swimming legs and the capability of operating in shoals is designed, with the abilities of walking and swimming under water, as shown in Fig. 2.

    The robot consists of three pairs of walking legs,a pair of swimming legs, embedded controller modules, a perception unit, a wireless communication part and an energy system. The joints of the crablike robot are driven by waterproof servos. The control system and the energy system are wrapped in a sealed cabin. The communication between the actuators and the control system is through the watertight connector arranged outside the sealed cabin.

    Fig. 3 Coordinate system of the swimming leg

    In view of the fact that the BM and CP joints with a small angle span during the flapping of the swimming legs can be neglected in the process of the design, a simplified swimming leg model is built, as shown in Fig. 3. Four waterproof servos are combined into an actuator to drive the hydrofoil. These servos are functioned to generate the feathering, lead-lag,flapping and rowing motions. The length of the swimming legs is 280 mm, and the proportions of the Coxa, the Basis, the Propus and the Dacty1 are 0.21,0.11, 0.43 and 0.25, respectively. The rotations of the TC, CB, BP and PD joints are controlled by the servos with angle spans of φTC(-90°≤φTC≤90°),φCB(-45°≤φCB≤90°), φBP(-90°≤φBP≤45°) and φPD(-30°≤φPD≤30°).

    2.2 Kinematic model of swimming leg

    Before analyzing the hydrodynamics of the swimming leg, a kinematic model of the swimming leg should be established. Then the mathematical relationships between the joint angles and the attack angles of the end-hydrofoil should be defined. Further,the joint trajectory plan should be made. In order to analyze the kinematic characteristics of the swimming leg, the inertial reference frame G- XYZ and the′ ′′′ for the Propus and the moving reference frame O-X Y Z moving reference frame O-X Y Z′′ ′′′′′′ for the Dactyl are defined accordingly as shown in Fig. 3.

    Fig. 4 (Color online) Hydrodynamic experiment platform of swimmeret

    The feathering angle, φTC, the lead-lag angle,, the flapping angle, φBP, and the rowing angle,, defined in Fig. 4, are sinusoidally varied:

    where φTCA, φCBA, φBPAand φPDAdenote the amplitudes of the feathering motion, the lead-lag motion, the flapping motion and the rowing motion ,respectively. φTCC, φCBC, φBPCand φPDCare the biased angles of the feathering motion, the lead-lag motion, the flapping motion and the rowing motion,respectively. f and t denote the motion frequency of the swimming leg and the time. ΔφCBis the phase difference between the lead-lag motion and the feathering motion. ΔφBPis the phase difference between the flapping motion and the feathering motion. ΔφPDis the phase difference between the rowing motion and the feathering motion.

    The asymmetric oscillations during the lead-lag motion and the flapping motion of the swimming legs are obtained through the real-time state monitoring and the data analysis of the marine crab movement. In general the down stroke motion takes less time than the upstroke motion. The lead motion takes less time than the lag motion. Due to the asymmetric motion,the effect of the added mass is enhanced and the aggregated vortex rope is merged into the hydrofoil trailing vortex, which increases the pressure and the strength of the jet in the wake. Consequently the thrust force is increased. A mathematical model of the pectoral fin oscillating in an asymmetric motion was set up by Su et al.[14]. The hydrodynamic performance of the pectoral-fin was analyzed when the asymmetric coefficient is 0.5. But in this model. a sudden change of the joint velocity cannot be realized in practice.These problems are solved by developing a new motion-time ratio model of the oscillating asymmetric motion of the swimming legs.

    where subscript i stands for the TC joint, the CB joint, the BP joint and the PD joint, λ stands for asymmetric coefficients. In the drag-based swimming modeand in the lift-based swimming mode, where tF, tB, tUand tD,respectively, stand for the lead motion time, the lag motion time, the upstroke time and the down stroke time.

    By the representation of the Denavit-Hartenberg for the robot manipulator, it is possible to define the homogeneous transformation matrices for adjacent links. The location of the end hydrofoil with respect to the inertial coordinate system is expressed as

    And the homogeneous transformation matrices are

    2.3 Hydrodynamic coefficient and efficiency of swimming legs

    According to the above analysis and the coordinate system, the thrust force coefficient CX,the lateral force coefficient CY, the lift force coefficient CZ, the moment coefficient in XGdirection, the moment coefficient in YGdirection CMYand the moment coefficient in ZGdirection CMZare defined as follows:

    The propulsive efficiency of a swimming leg,which is generally difficult to measure because of the large dynamic loads in both lift and thrust directions,is defined as the power output divided by the power input

    where

    and α is the angle between the lateral axis of the inertial coordinate system and the lateral axis of the body coordinate system.

    The Strouhal number related to the flapping hydrodynamic performance of the swimming leg[16]. It motion has an important influence on the is defined as

    where b is considered here as the estimated length of the wake width, and U is the speed of the swimming leg.

    The Reynolds number is expressed as[16]

    where C is the chord length of the swimming leg and ν is the kinematic viscosity of the fluid.

    3. Numerical calculation principle and method

    3.1 Governing equations

    The effects of the swimming leg oscillations on the fluid flow around such a body are highly unsteady,accompanied with different scale vortices, and a detailed analysis of the fluid-structure interactions is desirable. To simulate the time dependent unsteady viscous flow around an oscillating swimming leg, the unsteady Navier-Stokes equations with an incompressible assumption are used.

    The governing equations for the unsteady incompressible flow are as follows:

    where uidenotes the transient velocity component in i direction,and, respectively, represent the mean velocity components in i and j directions,ρ and μ are the density of the fluid and the dynamic viscosity,is the Reynolds shearing,is the source term,is the mean pressure. The modeled transport equations for k and ε in the practicable model are:

    where Gkis the generation of k due to the mean velocity gradients,

    3.2 Computational domain and grids

    Because for the structured grid method, high quality grids are required, while the distortion of grids is magnificent when the bionic swimming leg is in movement. The precision of the computation using the structured grid method is not good enough. So in this paper the unstructured grid algorithm is used to guarantee the mesh quality near the swimming leg during its movement. In this method the grid zone is divided into two different domains. The interior domain created by the swimming leg and its near domain are moving together and the grid in them will be unchanged. The exterior domain, i.e., the domain between the outer edge of the inner domain and the far field is filled with unstructured grids, and the grid there will be changed to accommodate the swimming leg motion.

    The spring based smoothing and Laplace local re-meshing method is applied to update the dynamic mesh during the motion of the flapping hydrofoil. To smooth the mesh, a value of 0.5 is used for the boundary node relaxation and a standard value of 0.001 is used for the convergence criterion and 0.01 for the Spring Constant Factor.

    The main function of the 4th-DOF of the swimming legs is to expand and recover the hydrofoil,which mainly affects the value of the hydrodynamic coefficient more than its overall trend. This paper investigates the influence of the first three rotational degrees of freedom on the hydrodynamic performance of the swimming legs. The rigid hydrofoils of the swimming legs are simplified in the simulation and experiment, without the end degree of freedom.

    4. Results and discussions

    4.1 Hydrodynamic experiments

    The major part of the hydrodynamic experimental platform is composed of the propulsion system, the control system and the hydrodynamic measurement system, as shown in Fig. 4. In the bionic swimming leg propulsion system, 4 waterproof servos HS-5646WP are used to control the joint movement. The hydrodynamic measurement system mainly consists of a 3-component force sensor, a CAN data acquisition card, and a dSPACE semi-physical simula- tion platform. The 3-axis piezoelectric force sensor is attached on the experimental platform to detect the lift,lateral and thrust forces, to be transferred to the dSPACE through the CAN bus. All signals are digitally filtered at 30 Hz before any calculation and plotting. The embedded control system based on the PC104 on the carriage controls the motions of the swimming legs through a servo motion control card.The hydrodynamic measurement platform is attached to a linear guide located above the free surface of water and can be translated at an adjustable velocity U. The hydrofoil of ABS material is made by the 3-D printing technique.

    4.2 Hydrodynamics of swimming leg

    In view of a close relationship between the reliability of the computation results using the CFD model and the quality of the mesh, in this section, the influence of the grid number in the computational domain on the computation accuracy is investigated.The thrust coefficient tends to be stable when the global grid numbers are over 3.1×106. In this case, the calculation error is less than 1%, within the allowable range, as compared to the cases of 4.66×106and 7.5×106grid numbers. So, with 3.1×106global grid numbers, the basic requirements for the 3-D numerical simulation can be met. More grid numbers may not improve the computation accuracy substantially withmore computation time.

    In order to verify the numerical methods in this study, the computational results of both the blade element method[17]and the CFD method are compared with experimental results. The numerical computation and experiment conditions are: U=0.2 m/s ,φTCA=35°, φCBA=35°, φBPA=40°, ΔφCB=0°,ΔφBP=90°, f=0.5 Hz and Re=1.0× 105, the unassigned parameters are set 0°, and the same is done in what follows.

    Fig. 5 Comparison of hydrodynamic coefficients during one cycle: (a) Time record of the thrust coefficient in one cycle,(b) Time record of the lift coefficient in one cycle

    Figure 5 shows the comparison of results of the FLUENT, the hydrodynamics experiment and the blade element method. It can be seen that the variation trends of the hydrodynamic coefficients of the two calculations are consistent with the experimental results. The extreme value and the inflection point are reached approximately at the same time, but with a certain gap, especially, for the peak value. There are three key factors related to the deviation between the experimental and numerical results. First, the experiment platform vibrates due to the movement of the swimming legs, which increases the interference between the swimming legs and the test platform.Second, the swimming legs are simplified as a rigid plate in the simulation. The motion of the actuator also generates some hydrodynamic force in the experiment besides the hydrofoil, i.e., the surface area of the hydrofoil increases. Third, compared with the viscosity model in the CFD method, the fluid separation and turbulence model in the experiment is more complex.

    4.3 Trail vortex structure

    The 3-D vortex contours of the swimming leg in a cycle are shown in Fig. 6.

    Fig. 6 (Color online) Vortices in an motion period

    Fig. 7 (Color online) The distribution of shedding vortices

    Due to the complexity of the 3-D vortex structure,this paper employs the slice method to analyze the trail vortex of the swimming legs. Figure 7 shows the contours of the vorticity ξyin the direction of Y-axis on the section of Y=100 mm , which depicts the consumption and the shedding of the trail vortex.There are two pairs of vortices on the top and the bottom, which shed sequentially in a motion period.Around 1/4T, an anti-clockwise vortex (ξy>0)and a clockwise vortex (ξy<0) shed below the swimming leg hydrofoil. Around 3/4T, an anticlockwise vortex (ξy>0) and a clockwise vortex(ξy<0) shed below the hydrofoil. By the comparison of the contour structure of the 2-D and 3-D airfoil in the case of flapping, there is a single-row reverse Karman Vortex Street in the wake of the 2-D airfoil edge, while a double-row vortex pair appears in the wake of the 3-D airfoil edge, which significantly improves the propulsion of the hydrofoil[18].

    4.4 Hydrodynamic performance of 2DOF swimming leg

    4.4.1 Effect of phase difference on hydrodynamic performance

    In this part, the effect of the phase difference between the feathering and flapping motions on the swimming leg hydrodynamics is studied by experimental and CFD methods. The motion parameters are as follows: U =0.2 m/s , φTCA=35°,φCBA=40°, f=0.5 Hz and ΔφBP=-120°-120°.

    Figure 8 shows the time variation of the hydrodynamic coefficient CXfor various phase differences ΔφBPbetween the feathering and flapping motions. It is seen that during one flapping period, the variation trends of the hydrodynamic coefficient CXare almost the same for different values of ΔφBP. There are two peak values, with the first peak value at 1/2T and the second peak value at the end of the cycle. We can further find that the change of ΔφBP, results in not only a change of the magnitude of the peak value of the thrust coefficient CX, but also a change of the time when the peak value occurs Meanwhile, the value of the thrust coefficient CXdecreases and the time when the peak value occurs are increased if ΔφBPincreases.

    Fig. 8 Hydrodynamic coefficients with various ΔφBP obtained by experimental method

    Figure 9 shows the curves of the average thrust coefficient CXmand the propulsive η against the phase difference ΔφBPfor the flapping frequency f equal to 0.5 Hz, 1.0 Hz and 1.5 Hz. From Fig. 9 we can see that with the increase of f, the range of the phase difference ΔφBPfor CXm≥0 extends. The average thrust coefficient (CXm≥0) reaches the maximum between the phase differences of 90° and 120° and the maximum value increases as the flapping frequency f increases. When the flapping frequency f keeps constant, the average propulsive efficiency η first rises and then decreases with the increase of ΔφBP. The maxima of the efficiency η at a fixed frequency f can be distinguished for f=0.5 Hz,1.0 Hz and 1.5 Hz, they, respectively, occur at the values of ΔφBPequal to 90°, 120°. The peak value of the efficiency η decreases with f.

    Fig. 9 Influence of ΔφBP on hydrodynamic performance obtained with the numerical method

    4.4.2 Effects of Strouhal number and flapping amplitude on hydrodynamic performance

    The motion parameters of the swimming leg are set as follows: φTCC=0°, φCPC=0°, φTCA=35°and φBPA=10°, 20°, 30° and 40°. The experimental values of the average thrust coefficient CXmand the propulsive efficiency η are plotted as the functions of the Strouhal number St in Figs. 10, 11. Figure 10 shows that the variation trends of the average thrust coefficient CXmversus the Strouhal number St for various φBPAare similar. The average thrust coefficient CXmmonotonically increases with St.Increasing St, one sees a continuous change from a drag-producing motion (CXm<0) to a thrustproducing motion (CXm>0). When St>0.2, the positive average thrust coefficient CXmcan be obtained. The maximum average thrust coefficient CXmis reached when St comes to 1.0, which will benefit the propulsion of the swimming leg. With the same St , when the flapping amplitude φBPAincreasesfrom φBPA=10° to φBPA=40°,CXmdecreases.

    Fig. 10 Influence of St on thrust coefficient obtained with numerical method

    Fig. 11 Influence of St on propulsive efficiency obtained with numerical method

    Figure 11 shows that with various φBPA, on the curves of the efficiency η vs. the Strouhal number St , a maximum efficiency η is observed with Stmaxbeing between 0.3 and 0.5. According to the analysis of the fish swimming body and posterior fin mode by Triantafyllou et al.[19], St is between 0.2 and 0.4.According to Taylor et al.[20], for bird?s wing in the flapping mode, its efficiency reaches the maximum when St is equal to a value in the range of 0.3-0.5, as in a good agreement with the St range of the swimming leg with a high efficiency. When St is increased, the efficiency η increases quickly until it reaches Stmax, then η decreases if St continues to increase. The maximum efficiency η is 53.52%when φCPA=10° and St=0.3. The corresponding average thrust coefficient CXmis very small and as small as 0.35. Through the above analysis, we can see that when the value of St is slightly larger than Stmax, a relatively high efficiency together with a high average thrust can be achieved. In the range of small values of φBPA, the efficiency η is high, but it decreases when φBPAis increased from 10° to 40°.

    The experimental results of the thrust coefficient CXversus time for various flapping amplitudes φBPAare as shown as in Fig. 12. It can be seen that the variation trends of CXwith various φBPAare similar. The time at which the peak value occurs is approximately the same but the magnitude of the peak value will increase if φBPAincreases.

    Fig. 12 Hydrodynamic coefficients for different φBPA obtained with experimental method

    4.4.3 Effect of flapping biased angle on hydrodynamic performance

    The parameters for the numerical simulation and the experiment are: φTCC=φCPC=0°, φTCA=35°,φBPA=40°, f=0.5 Hz, the value of the flapping biased angle φBPCis equal to 5°, 10°, 15° and 20°.

    Variations of the average thrust coefficient CXmand the average lift coefficient CZmare plotted in Fig.13, when the flapping biased angle φBPCvaries from 0° to 40°, St is equal to 0.8 and 1.0, respectively,where the horizontal axis represents CXmand the vertical axis is CZm.

    Fig. 13 Influence of φBPC on hydrodynamic performance obtained with numerical method

    We can see that with a small flapping biased angle φBPC, the lift force has an asymmetrical distributions in positive and negative regions.Approximately a zero average lift force can be obtained when the biased angle is zero. Applying a flapping biased angle will considerably increase the average lift coefficient CZm. The average thrust coefficient CXmdecreases with the increase of φBPC,and negative values of CXmare obtained, when φBPC≥33.4° for St=0.8 and φBPC≥37.8° for St=1.0. Hence, it is indicated that a high St will increase the thrust.

    4.4.4 Effect of feathering amplitude on hydrodynamic performance

    The parameters for the numerical simulation and the experiment are: φTCC=φCPC=0°, φBPA=40°,ΔφBP=90°, f=0.5 Hz and the feathering amplitude φTCAis equal to 10°, 20°, 30° and 40°, respectively.

    Numerical results of the average force coefficient and the propulsive efficiency are plotted as the functions of the feathering amplitude as shown in Fig.14. We can see that when φTCAis increased, the absolute values of CXm, CYmand CZmall increase but in different rates. The rate of CYmis the maximum, the rate of CZmis the minimum. The propulsive efficiency monotonously increases with the increase of φTCAfrom 0° to 40°.

    Fig. 14 Influence of φTCA on hydrodynamic performance obtained with numerical method

    4.4.5 Effect of flow velocity on hydrodynamic performance

    The operating conditions are: φTCC=φCPC=0°,φTCA=35°, φBPA=40°, ΔφBP=90°, f=1.0 Hz,the flow velocity increases from 0.1 m/s to 0.6 m/s with an increment of 0.1 m/s.

    Fig. 15 Influence of U on hydrodynamic performance obtained with numerical method

    Figure 15 shows that the average thrust and lateral force coefficients (absolute values) sharply decrease with the increase of the flow velocity when the flow velocity is less than 0.3 m/s. A very slight decrease of the thrust and lateral force coefficients(absolute values) can be seen when the flow velocity increases further in the large velocity range (0.3 m/s-0.6 m/s). The average lift coefficient approaches zero and remains constant with the increase of the flow velocity. Figure 18 shows that the propulsive efficiency increases with the increase of the flow velocity when the flow velocity is less than 0.3 m/s.After the maximum propulsive efficiency of 42.2% is achieved at U =0.3m/s , it starts to decrease until it is zero, where the average thrust coefficient is very small. Thus the flow velocity 0.2 m/s is selected in the next numerical calculation and experiment.

    4.5 Hydrodynamic performance of 3DOF swimming leg

    In this part the hydrodynamic performance of the 3DOF swimming leg is studied. The leg motion relative to the body consists of feathering, lead-lag and flapping motions.

    4.5.1 Effect of phase difference on hydrodynamic performance

    Figure 16 shows the average hydrodynamic coefficient and the propulsive efficiency against ΔφCBfrom -120° to 120° with various flapping frequencies f and the other parameters, φTCC=φCBC=φBPC=0°, φTCA=35°, φCBC=30°, φBPA=40°, ΔφBP=90°, U=0.2 m/s.

    We can see a slight change of the average hydrodynamic coefficient. The maximum value of CXm, is obtained when ΔφCB=0°. With the increase of ΔφCBin the range -120°≤ΔφCB≤-90°, the propulsive efficiency gradually decreases. The mini-mum propulsive efficiency is obtained at ΔφCB=-90°. With the increase of ΔφCBin the range-90°≤ΔφCB≤0° the propulsive efficiency gradually increases. The highest value of the propulsive efficiency is obtained at ΔφCB=0° and after that it decreases. Comparing the curves of the propulsive efficiency with 3 different frequencies, we can see that with a low flapping frequency (0.5 Hz), the propulsive efficiency is low and with the medium frequency(1.0 Hz), the propulsive efficiency is maximized. A higher frequency (1.5 Hz) will decrease the propulsive efficiency again. A maximum efficiency of around 34%is achieved at f=0.5 Hz and ΔφCB=0°. Comparing with the result from the hydrodynamic experiment of the 2DOF swimming leg, it is shown that adding the lead-lag motion (3DOF) cannot increase the thrust significantly and cannot decrease the propulsive efficiency distinctly either.

    Fig. 16 Influence of ΔφCB on hydrodynamic performance obtained with numerical method

    4.5.2 Effect of lead-lag biased angle on hydrodynamic performance

    Parameters for the numerical simulation and experiment are: φTCC=φBPC=0°, φTCA=35°,φCBA=30°, φBPA=40°, ΔφBP=90°, ΔφCB=0°,U=0.2 m/s , f=0.5 Hz and the range of the lead-lag biased angle is φCBC=10°-30°.

    Figure 17 shows the average hydrodynamic coefficient and propulsive efficiency as functions of the lead-lag biased angle for the swimming leg. The average lift and lateral force coefficients are proportional to the lead-lag biased angle between 0°and 40°. Increasing the lead-lag biased angle φCBC,we also observe that the average thrust coefficient first increases and then decreases. The maximum average thrust coefficient is CXm=5.95 at φCBC=30°. The curve of the pr opu lsive eff icienc y has the same variation trend astheaveragethrustcoefficient.The peak propulsive efficiency occurs at φCBC=20°.

    Fig. 17 Influence of φCBC on hydrodynamic performance obtained with numerical method

    4.5.3 Effect of lead-lag amplitude on hydrodynamic performance

    Figure 18 shows the average hydrodynamic coefficient and propulsive efficiency versus φCBAfor φTCC=φCBC=φBPC=0°, φTCA=35°, φBPC=40°,ΔφBP=90°, ΔφCB=0°, U =2.0 m/s , f=1.0 Hz .

    Fig. 18 Influence of φCBA on hydrodynamic performance obtained with numerical method

    As φCBAincreases from 0° to 40°, the variances of CXmand CYmare very small and the absolute value of CZmincreases considerably. The propulsive efficiency, η, quickly decreases at the beginning(0°<φCBA<5°) until it becomes negative and then decreases gently.

    4.5.4 Effect of foil thickness on hydrodynamic performance

    In order to investigate the influence of the foil thickness on the hydrodynamic performance, three thicknesses (1.5 mm, 3.0 mm and 4.5 mm) are considered in the experiment. The experimental results are as shown as in Figs. 19, 20. When the flapping frequency is 0.5 Hz, the average thrust coefficient is very small for all foil thicknesses. How- ever, when the flapping frequency is high (1.5 Hz), a thin foil(1.5 mm) results in a high thrust coefficient (9.0) and a thick foil (4.5 mm) results in a low thrust coefficient(7.5).

    Fig. 19 Influence of foil thickness on hydrodynamic coefficients obtained with numerical method

    Fig. 20 Influence of foil thickness on propulsive efficiency obtained with numerical method

    Figure 20 shows that the propulsive efficiency decreases with the increase of the foil thickness at a low frequency (0.5 Hz). But the foil thickness has a small effect on the propulsive efficiency when the flapping frequency is high (1.5 Hz). The maximum propulsive efficiency of 19.8% is achieved at f=1.0 Hz when the foil thickness is 1.5 mm. The results suggest that the crab-like robot can swim faster and more efficiently with a thin foil (1.5 mm) than a thick one (4.5 mm).

    5. Conclusions

    In this study, a leg-paddle coupling crablike robot is developed. The hydrodynamic performance of the swimming leg is studied with numerical simulations and experiments. The effects of kinematic parameters and foil thickness on the hydrodynamic performance are also discussed. The following conclusions are drawn:

    (1) The numerical simulation and experiment for the oscillating foil of a swimming leg provide systematic data on thrust, lateral and lift forces generated harmonically. The propulsion tests show the optimum thrust and propulsive efficiency for a certain combination of parameters like the phase difference,the flapping amplitude, the biased angle of the coupling movement and the Strouhal number.

    (2) The double-row vortex pair appears during the flapping motion of the hydrofoil, which significantly improves the propulsion of the hydrofoil.

    (3) The foil thickness has effects on the hydrodynamic performance of the swimming leg. Large thrust and high efficiency can be obtained by a thin foil (1.5 mm).

    猜你喜歡
    陳曦王剛海龍
    Lower bound on the spread of valley splitting in Si/SiGe quantum wells induced by atomic rearrangement at the interface
    封面人物
    天工(2021年2期)2021-03-03 07:29:16
    葉海龍,你別裝啦
    善于總結(jié)化難為易
    The Influence of English Reform in China on English Teaching and Learning
    海龍卷是什么
    Analysis of View of Life and Death of Christianity on Western Culture
    海龍卷是什么
    原野上的花
    世經(jīng)版圖
    色综合色国产| 如何舔出高潮| 亚洲国产欧美在线一区| 身体一侧抽搐| 亚洲人成网站在线播放欧美日韩| 老师上课跳d突然被开到最大视频| 性色avwww在线观看| 99久久精品热视频| 亚洲一区二区三区色噜噜| 精品久久久久久久久久免费视频| 一边摸一边抽搐一进一小说| 亚洲一级一片aⅴ在线观看| 人妻夜夜爽99麻豆av| 免费观看精品视频网站| 亚洲人与动物交配视频| 听说在线观看完整版免费高清| 国产av在哪里看| 日韩,欧美,国产一区二区三区 | 日日摸夜夜添夜夜爱| 亚洲内射少妇av| 国产真实伦视频高清在线观看| 伦理电影大哥的女人| 观看免费一级毛片| 在线免费十八禁| 欧美一区二区亚洲| 婷婷六月久久综合丁香| 亚洲欧美精品自产自拍| 久久这里有精品视频免费| 欧美丝袜亚洲另类| 卡戴珊不雅视频在线播放| 日韩三级伦理在线观看| 精品人妻偷拍中文字幕| 搞女人的毛片| 欧美高清性xxxxhd video| 能在线免费看毛片的网站| 又粗又硬又长又爽又黄的视频 | 免费电影在线观看免费观看| 大型黄色视频在线免费观看| 国产精品乱码一区二三区的特点| 亚洲图色成人| 国产精品久久久久久精品电影小说 | 99热网站在线观看| 大型黄色视频在线免费观看| 国产探花极品一区二区| 中文在线观看免费www的网站| 国产精品久久久久久精品电影小说 | 成人av在线播放网站| 亚洲av第一区精品v没综合| 日韩欧美一区二区三区在线观看| av又黄又爽大尺度在线免费看 | 国产精品永久免费网站| 99在线视频只有这里精品首页| 国产av不卡久久| 亚洲成人中文字幕在线播放| 不卡一级毛片| 久久久久免费精品人妻一区二区| 99热这里只有精品一区| 日韩亚洲欧美综合| 免费大片18禁| 我的女老师完整版在线观看| 国产三级中文精品| 国语自产精品视频在线第100页| 黄色欧美视频在线观看| 91精品国产九色| 春色校园在线视频观看| 噜噜噜噜噜久久久久久91| 久久精品国产自在天天线| 国产一级毛片七仙女欲春2| 尾随美女入室| 老女人水多毛片| 欧美3d第一页| 久久国内精品自在自线图片| 欧美一区二区亚洲| 亚洲自偷自拍三级| 最好的美女福利视频网| 一级黄色大片毛片| 欧美高清成人免费视频www| 精品久久久久久成人av| 在线免费十八禁| 1000部很黄的大片| 人妻制服诱惑在线中文字幕| 99精品在免费线老司机午夜| 女的被弄到高潮叫床怎么办| 麻豆成人午夜福利视频| 91久久精品国产一区二区成人| 麻豆国产97在线/欧美| 搡女人真爽免费视频火全软件| 日韩精品青青久久久久久| 三级经典国产精品| 久久久久久久久中文| 亚洲天堂国产精品一区在线| 韩国av在线不卡| 国产亚洲av嫩草精品影院| 日日干狠狠操夜夜爽| 变态另类丝袜制服| 亚洲欧洲国产日韩| 国产三级在线视频| 黄片无遮挡物在线观看| 成人毛片60女人毛片免费| 网址你懂的国产日韩在线| 国产在线男女| 日韩欧美国产在线观看| 嘟嘟电影网在线观看| 麻豆国产av国片精品| 欧美在线一区亚洲| 丰满乱子伦码专区| 国产精品一及| 只有这里有精品99| 成人特级黄色片久久久久久久| 日韩国内少妇激情av| 免费一级毛片在线播放高清视频| 高清在线视频一区二区三区 | 国产成年人精品一区二区| 婷婷色综合大香蕉| 男女下面进入的视频免费午夜| 中文精品一卡2卡3卡4更新| 亚洲av免费在线观看| 亚洲精品成人久久久久久| 91久久精品电影网| 日本成人三级电影网站| 内射极品少妇av片p| 国产成人午夜福利电影在线观看| 国产日韩欧美在线精品| 蜜桃久久精品国产亚洲av| 麻豆av噜噜一区二区三区| 亚洲成人久久爱视频| 午夜福利在线观看吧| 校园人妻丝袜中文字幕| 九九热线精品视视频播放| 国产精品一及| 草草在线视频免费看| 美女大奶头视频| 国产亚洲5aaaaa淫片| 国产高潮美女av| 午夜福利在线观看免费完整高清在 | 久久久久久九九精品二区国产| 日本成人三级电影网站| 看免费成人av毛片| 91狼人影院| 小蜜桃在线观看免费完整版高清| 国产日韩欧美在线精品| 国产 一区 欧美 日韩| 国产精品永久免费网站| 女的被弄到高潮叫床怎么办| 国产69精品久久久久777片| 久久人人爽人人片av| 国产精品久久久久久精品电影小说 | 国产精品久久久久久久久免| 国产视频内射| 最近的中文字幕免费完整| 一本一本综合久久| 久久久久九九精品影院| 毛片一级片免费看久久久久| 日日干狠狠操夜夜爽| 一夜夜www| 我要搜黄色片| 99九九线精品视频在线观看视频| 男女边吃奶边做爰视频| 国产av不卡久久| 又粗又硬又长又爽又黄的视频 | 欧美成人a在线观看| 少妇猛男粗大的猛烈进出视频 | 美女xxoo啪啪120秒动态图| 91久久精品国产一区二区三区| 99在线人妻在线中文字幕| 亚洲精品色激情综合| 2022亚洲国产成人精品| 国产爱豆传媒在线观看| 99久久久亚洲精品蜜臀av| 国内精品美女久久久久久| 全区人妻精品视频| 日本成人三级电影网站| av天堂中文字幕网| 国产高清不卡午夜福利| 在线观看66精品国产| 欧美日韩综合久久久久久| .国产精品久久| 国语自产精品视频在线第100页| 色综合色国产| 国语自产精品视频在线第100页| 国语自产精品视频在线第100页| 看十八女毛片水多多多| 国产精品1区2区在线观看.| av在线蜜桃| 婷婷六月久久综合丁香| 精品免费久久久久久久清纯| 国产黄片视频在线免费观看| 中文字幕熟女人妻在线| 一级毛片久久久久久久久女| 嫩草影院新地址| 91久久精品国产一区二区成人| 国产成人a∨麻豆精品| 久久午夜福利片| 日韩欧美精品免费久久| 最近手机中文字幕大全| 高清毛片免费观看视频网站| 欧美成人一区二区免费高清观看| 亚洲精品久久国产高清桃花| 美女cb高潮喷水在线观看| 久久午夜亚洲精品久久| 亚洲欧美中文字幕日韩二区| 麻豆成人午夜福利视频| 热99re8久久精品国产| 欧美色欧美亚洲另类二区| 亚洲一区高清亚洲精品| 99久国产av精品| 亚洲图色成人| 真实男女啪啪啪动态图| 国产成人精品一,二区 | 国内揄拍国产精品人妻在线| 久久婷婷人人爽人人干人人爱| 日韩欧美一区二区三区在线观看| 免费电影在线观看免费观看| АⅤ资源中文在线天堂| 综合色av麻豆| 免费观看人在逋| 久久精品国产亚洲网站| 久久精品国产99精品国产亚洲性色| 极品教师在线视频| 亚洲人成网站在线播| 国产伦一二天堂av在线观看| 九草在线视频观看| av天堂中文字幕网| 亚洲一区二区三区色噜噜| 久久精品人妻少妇| 深夜a级毛片| 久久婷婷人人爽人人干人人爱| 黑人高潮一二区| 国产一区二区激情短视频| 国产精品精品国产色婷婷| 两个人视频免费观看高清| 性插视频无遮挡在线免费观看| 最近手机中文字幕大全| 国产成人aa在线观看| 久久久国产成人精品二区| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人妻熟人妻熟丝袜美| 久久久a久久爽久久v久久| 青春草视频在线免费观看| 国产精品不卡视频一区二区| 婷婷亚洲欧美| 国产一级毛片七仙女欲春2| 此物有八面人人有两片| www日本黄色视频网| 丰满人妻一区二区三区视频av| 国产一区二区亚洲精品在线观看| 亚洲av一区综合| 亚洲五月天丁香| 国产人妻一区二区三区在| 一级av片app| 国产 一区 欧美 日韩| 久久久久久久久中文| 欧美不卡视频在线免费观看| 寂寞人妻少妇视频99o| 男人舔奶头视频| 一区福利在线观看| 国产成人精品婷婷| 亚洲精品影视一区二区三区av| 男人和女人高潮做爰伦理| 99热这里只有是精品在线观看| 久久精品国产亚洲av天美| 99热精品在线国产| 成人永久免费在线观看视频| 国产精品久久久久久久电影| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区三区| 非洲黑人性xxxx精品又粗又长| 国产免费男女视频| 欧美一级a爱片免费观看看| 久久久久免费精品人妻一区二区| 变态另类成人亚洲欧美熟女| 亚洲成人久久性| 五月伊人婷婷丁香| 亚洲国产欧洲综合997久久,| 日韩成人av中文字幕在线观看| 禁无遮挡网站| 秋霞在线观看毛片| 亚洲图色成人| 久久久久久久久中文| 欧美最新免费一区二区三区| 男女啪啪激烈高潮av片| 久久这里只有精品中国| 色综合站精品国产| 一个人观看的视频www高清免费观看| 亚洲国产精品成人久久小说 | 亚洲精品亚洲一区二区| 草草在线视频免费看| 日韩强制内射视频| 91久久精品国产一区二区成人| 国产精品一区www在线观看| 成人亚洲精品av一区二区| 可以在线观看毛片的网站| 精品少妇黑人巨大在线播放 | 欧美不卡视频在线免费观看| ponron亚洲| 亚洲欧洲国产日韩| a级毛片a级免费在线| 国产精品久久视频播放| 人妻系列 视频| 99热这里只有是精品50| 日韩,欧美,国产一区二区三区 | 熟女电影av网| 欧美日韩在线观看h| 亚洲自偷自拍三级| 女人十人毛片免费观看3o分钟| 最近2019中文字幕mv第一页| 日本三级黄在线观看| 欧美日本视频| 成年女人看的毛片在线观看| 黄片wwwwww| 日韩强制内射视频| 日本欧美国产在线视频| 晚上一个人看的免费电影| 美女被艹到高潮喷水动态| h日本视频在线播放| 此物有八面人人有两片| 深夜精品福利| av卡一久久| 欧美丝袜亚洲另类| 啦啦啦观看免费观看视频高清| 国产精品.久久久| 亚洲真实伦在线观看| 亚洲自偷自拍三级| 一个人观看的视频www高清免费观看| 99热全是精品| 国产探花在线观看一区二区| 老熟妇乱子伦视频在线观看| 亚洲真实伦在线观看| 又粗又硬又长又爽又黄的视频 | 男人狂女人下面高潮的视频| 天堂av国产一区二区熟女人妻| 全区人妻精品视频| 日韩精品有码人妻一区| 久久综合国产亚洲精品| 禁无遮挡网站| 亚洲图色成人| 久久鲁丝午夜福利片| 国产亚洲5aaaaa淫片| 少妇的逼好多水| 高清毛片免费观看视频网站| 欧美另类亚洲清纯唯美| 成人永久免费在线观看视频| 久久精品人妻少妇| 欧美精品国产亚洲| 最后的刺客免费高清国语| 亚洲av电影不卡..在线观看| 国产成人一区二区在线| 此物有八面人人有两片| 长腿黑丝高跟| 99riav亚洲国产免费| 午夜精品一区二区三区免费看| 免费观看a级毛片全部| 在线观看一区二区三区| 69av精品久久久久久| 欧美3d第一页| 亚洲国产日韩欧美精品在线观看| 欧美精品一区二区大全| 国产一区亚洲一区在线观看| 成人无遮挡网站| 99久久人妻综合| 老熟妇乱子伦视频在线观看| 哪个播放器可以免费观看大片| 免费看美女性在线毛片视频| 乱系列少妇在线播放| 成人一区二区视频在线观看| 成人特级黄色片久久久久久久| 国产精品久久久久久精品电影| 日韩av在线大香蕉| 晚上一个人看的免费电影| 伦精品一区二区三区| 最好的美女福利视频网| 激情 狠狠 欧美| 韩国av在线不卡| 18+在线观看网站| 91午夜精品亚洲一区二区三区| АⅤ资源中文在线天堂| 亚洲欧美日韩高清专用| 真实男女啪啪啪动态图| 亚洲中文字幕日韩| 亚洲一级一片aⅴ在线观看| 亚洲国产最新在线播放| 九九爱精品视频在线观看| 久久鲁丝午夜福利片| 日韩中字成人| 国产精品国产三级国产av玫瑰| 丝袜美足系列| 亚洲精品乱码久久久v下载方式| 伦理电影大哥的女人| 最新中文字幕久久久久| 国产av国产精品国产| 三上悠亚av全集在线观看| 大香蕉久久网| 亚洲美女视频黄频| 精品卡一卡二卡四卡免费| 久久精品国产鲁丝片午夜精品| 国产高清不卡午夜福利| 亚洲av.av天堂| 最后的刺客免费高清国语| 精品亚洲成国产av| 亚洲第一区二区三区不卡| 在线观看三级黄色| av女优亚洲男人天堂| 秋霞伦理黄片| 亚洲av中文av极速乱| 国产精品一区二区在线不卡| 久久99热6这里只有精品| 尾随美女入室| 免费观看a级毛片全部| av在线播放精品| 免费黄网站久久成人精品| 亚洲人成77777在线视频| 我的女老师完整版在线观看| 欧美 亚洲 国产 日韩一| 久久狼人影院| 婷婷色麻豆天堂久久| 免费大片18禁| 午夜精品国产一区二区电影| .国产精品久久| 桃花免费在线播放| 国产精品久久久久成人av| 久久精品熟女亚洲av麻豆精品| 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图 | 日韩欧美一区视频在线观看| 亚洲精品视频女| 国产成人av激情在线播放 | 99久久精品一区二区三区| 18在线观看网站| 久久青草综合色| 高清毛片免费看| 狠狠婷婷综合久久久久久88av| 少妇的逼好多水| 国产乱来视频区| 免费看光身美女| 亚洲国产精品专区欧美| 国产精品嫩草影院av在线观看| h视频一区二区三区| 只有这里有精品99| 亚洲精品av麻豆狂野| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品aⅴ在线观看| 中国国产av一级| 中文字幕制服av| 中文字幕人妻丝袜制服| 色哟哟·www| 蜜桃国产av成人99| 欧美最新免费一区二区三区| 亚洲av福利一区| 国产精品国产三级国产专区5o| 日韩一区二区视频免费看| 日韩人妻高清精品专区| 色视频在线一区二区三区| 99久久人妻综合| 欧美精品一区二区大全| 免费大片18禁| 伊人久久国产一区二区| 亚洲精品aⅴ在线观看| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 在线看a的网站| 人人澡人人妻人| 亚洲图色成人| 91久久精品电影网| 日本爱情动作片www.在线观看| 精品人妻一区二区三区麻豆| 伊人亚洲综合成人网| 久久99热这里只频精品6学生| videos熟女内射| 国产黄色免费在线视频| 少妇猛男粗大的猛烈进出视频| 国产色婷婷99| 国产精品无大码| 日韩免费高清中文字幕av| 免费看不卡的av| 亚洲国产av影院在线观看| 国产 精品1| 国产欧美亚洲国产| 欧美精品国产亚洲| 国产无遮挡羞羞视频在线观看| 九九在线视频观看精品| 中文字幕最新亚洲高清| 久久精品国产亚洲av天美| 欧美精品人与动牲交sv欧美| 国精品久久久久久国模美| 日韩视频在线欧美| 男人操女人黄网站| 国产亚洲精品久久久com| 日韩成人伦理影院| 桃花免费在线播放| 国产日韩欧美视频二区| 国产男女超爽视频在线观看| av一本久久久久| 亚洲av在线观看美女高潮| 亚洲精品乱码久久久久久按摩| 99久久综合免费| 国产精品人妻久久久久久| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 夫妻性生交免费视频一级片| 激情五月婷婷亚洲| 亚洲不卡免费看| 亚洲精品亚洲一区二区| 一区二区三区免费毛片| kizo精华| 97超视频在线观看视频| 亚洲精品中文字幕在线视频| 我的老师免费观看完整版| 日韩成人av中文字幕在线观看| 国产免费又黄又爽又色| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 黄色一级大片看看| 午夜久久久在线观看| 久久久久久久国产电影| 亚洲激情五月婷婷啪啪| 最后的刺客免费高清国语| 久久久久久久久久成人| 欧美精品国产亚洲| 久久久久久人妻| 国产精品国产三级专区第一集| 18禁裸乳无遮挡动漫免费视频| 狂野欧美白嫩少妇大欣赏| 精品国产露脸久久av麻豆| 天堂8中文在线网| 美女中出高潮动态图| 99热国产这里只有精品6| 99热6这里只有精品| 亚洲精品乱码久久久久久按摩| 精品久久久精品久久久| 亚洲久久久国产精品| 国产精品三级大全| 乱码一卡2卡4卡精品| 美女主播在线视频| 午夜激情福利司机影院| 国产精品久久久久久久电影| 免费观看性生交大片5| 国产成人午夜福利电影在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 春色校园在线视频观看| 成年人免费黄色播放视频| 国产视频首页在线观看| 免费少妇av软件| 亚洲国产av影院在线观看| 蜜桃国产av成人99| 亚洲一区二区三区欧美精品| 国产精品久久久久成人av| videosex国产| 视频中文字幕在线观看| 少妇的逼水好多| 国产成人精品婷婷| 在线观看人妻少妇| 少妇人妻 视频| 香蕉精品网在线| 91国产中文字幕| 少妇人妻精品综合一区二区| 制服丝袜香蕉在线| 精品国产国语对白av| www.av在线官网国产| 人人澡人人妻人| 国产成人一区二区在线| 国产成人免费无遮挡视频| 人体艺术视频欧美日本| 国内精品宾馆在线| 老司机影院毛片| 久久午夜福利片| 女性生殖器流出的白浆| 久久99热这里只频精品6学生| a级片在线免费高清观看视频| 多毛熟女@视频| 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久久丰满| 国产极品粉嫩免费观看在线 | 在线观看免费视频网站a站| 女人精品久久久久毛片| 老熟女久久久| 99热国产这里只有精品6| 美女视频免费永久观看网站| 国产黄频视频在线观看| 99热6这里只有精品| 日韩av不卡免费在线播放| 亚洲国产色片| 我的老师免费观看完整版| 亚洲av免费高清在线观看| 你懂的网址亚洲精品在线观看| 欧美日韩国产mv在线观看视频| 亚洲经典国产精华液单| 91精品伊人久久大香线蕉| 黄色一级大片看看| 亚洲欧美日韩卡通动漫| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品999| 亚洲情色 制服丝袜| 尾随美女入室| 国产综合精华液| 99re6热这里在线精品视频| 国产综合精华液| 亚洲在久久综合| 久久99热这里只频精品6学生| 在线精品无人区一区二区三| xxxhd国产人妻xxx| 成人亚洲欧美一区二区av| 午夜老司机福利剧场| 在线天堂最新版资源| 高清欧美精品videossex| 成年美女黄网站色视频大全免费 | 伦理电影大哥的女人| 美女福利国产在线| 国产精品一区二区在线观看99| 精品卡一卡二卡四卡免费| 欧美日韩综合久久久久久| 日韩免费高清中文字幕av|