• <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)版圖
    久久久精品免费免费高清| 日本wwww免费看| 日本av手机在线免费观看| 曰老女人黄片| 男女无遮挡免费网站观看| 国产成人精品在线电影| 欧美成狂野欧美在线观看| 午夜激情久久久久久久| 成年版毛片免费区| 一区福利在线观看| 亚洲黑人精品在线| 一级片'在线观看视频| 在线观看舔阴道视频| 三上悠亚av全集在线观看| 久久久欧美国产精品| 亚洲欧美一区二区三区黑人| 亚洲综合色网址| 无限看片的www在线观看| 亚洲伊人色综图| 久久久水蜜桃国产精品网| 欧美日韩成人在线一区二区| 99热这里只有精品一区 | 88av欧美| 日韩欧美三级三区| 免费大片18禁| 欧美日韩一级在线毛片| 亚洲第一电影网av| 亚洲国产日韩欧美精品在线观看 | 中文亚洲av片在线观看爽| 日本与韩国留学比较| 天堂√8在线中文| 亚洲av第一区精品v没综合| 亚洲美女视频黄频| 国产真人三级小视频在线观看| 成人国产一区最新在线观看| 国产一区二区在线观看日韩 | 国产熟女xx| 免费观看精品视频网站| 久久久久久国产a免费观看| 亚洲中文av在线| 亚洲国产色片| 中文字幕av在线有码专区| 午夜福利高清视频| 少妇裸体淫交视频免费看高清| 久久久久性生活片| 国产亚洲精品一区二区www| 欧美又色又爽又黄视频| 亚洲九九香蕉| 亚洲成a人片在线一区二区| 色综合婷婷激情| 国产精品久久久久久亚洲av鲁大| 在线观看美女被高潮喷水网站 | 国产免费av片在线观看野外av| 欧美zozozo另类| 精品日产1卡2卡| 欧美不卡视频在线免费观看| 欧美黄色片欧美黄色片| 亚洲成人久久爱视频| 美女被艹到高潮喷水动态| 国产高潮美女av| 床上黄色一级片| 可以在线观看的亚洲视频| 国产av麻豆久久久久久久| 日韩人妻高清精品专区| 男女之事视频高清在线观看| 国产高清videossex| 中文亚洲av片在线观看爽| 亚洲,欧美精品.| 男人的好看免费观看在线视频| 久久热在线av| 亚洲成人久久性| 麻豆国产97在线/欧美| 精品熟女少妇八av免费久了| 亚洲一区二区三区不卡视频| avwww免费| 成人午夜高清在线视频| 久久人妻av系列| 一区二区三区激情视频| 亚洲18禁久久av| 麻豆成人午夜福利视频| 人人妻人人看人人澡| 日本五十路高清| 叶爱在线成人免费视频播放| svipshipincom国产片| 国产亚洲av高清不卡| 免费电影在线观看免费观看| 熟女人妻精品中文字幕| 无限看片的www在线观看| av国产免费在线观看| 99久久无色码亚洲精品果冻| 天堂网av新在线| 黄色 视频免费看| 女生性感内裤真人,穿戴方法视频| 可以在线观看的亚洲视频| 亚洲av成人一区二区三| 国产一区在线观看成人免费| 欧美黑人巨大hd| 一本精品99久久精品77| 日本免费一区二区三区高清不卡| 一个人免费在线观看电影 | 天堂动漫精品| 亚洲成a人片在线一区二区| 一个人看视频在线观看www免费 | 给我免费播放毛片高清在线观看| 91在线观看av| 99re在线观看精品视频| 黄片大片在线免费观看| 欧美高清成人免费视频www| 最好的美女福利视频网| 18美女黄网站色大片免费观看| 久久热在线av| 亚洲乱码一区二区免费版| 91在线精品国自产拍蜜月 | 久久国产精品影院| 亚洲精品国产精品久久久不卡| 无人区码免费观看不卡| 精品福利观看| 九色成人免费人妻av| 免费在线观看成人毛片| 亚洲欧美精品综合一区二区三区| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| www日本在线高清视频| 久久精品国产亚洲av香蕉五月| 俺也久久电影网| 国产伦精品一区二区三区四那| 亚洲乱码一区二区免费版| 免费搜索国产男女视频| 脱女人内裤的视频| www日本黄色视频网| 国模一区二区三区四区视频 | 日韩欧美免费精品| 制服人妻中文乱码| 午夜免费成人在线视频| 欧美xxxx黑人xx丫x性爽| 精品一区二区三区av网在线观看| 2021天堂中文幕一二区在线观| 亚洲av熟女| 国产又黄又爽又无遮挡在线| 免费看a级黄色片| 色播亚洲综合网| 国产高清三级在线| 日韩成人在线观看一区二区三区| 日韩高清综合在线| 99国产精品99久久久久| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| 国产精品自产拍在线观看55亚洲| av黄色大香蕉| 国产三级在线视频| 欧美乱妇无乱码| 在线国产一区二区在线| avwww免费| 亚洲色图 男人天堂 中文字幕| 亚洲av成人一区二区三| 在线免费观看的www视频| 熟女人妻精品中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一区二区三区四区久久| 97超级碰碰碰精品色视频在线观看| 99国产精品一区二区三区| 日韩欧美一区二区三区在线观看| 美女高潮的动态| 黄色成人免费大全| 极品教师在线免费播放| 亚洲天堂国产精品一区在线| 天堂影院成人在线观看| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久com| 亚洲av片天天在线观看| 亚洲美女视频黄频| 两人在一起打扑克的视频| 亚洲色图av天堂| 亚洲熟妇中文字幕五十中出| 午夜精品在线福利| 国产高清videossex| 国产免费男女视频| 神马国产精品三级电影在线观看| 最新中文字幕久久久久 | 免费看日本二区| 九九在线视频观看精品| 精品免费久久久久久久清纯| 国产在线精品亚洲第一网站| 最新中文字幕久久久久 | 亚洲人成网站在线播放欧美日韩| 综合色av麻豆| 国产av一区在线观看免费| 一进一出抽搐gif免费好疼| 亚洲av第一区精品v没综合| 真实男女啪啪啪动态图| 美女免费视频网站| 最近在线观看免费完整版| 亚洲专区字幕在线| 亚洲精品美女久久久久99蜜臀| 亚洲一区高清亚洲精品| 日韩欧美 国产精品| а√天堂www在线а√下载| 亚洲av成人精品一区久久| 国产三级在线视频| 在线观看免费午夜福利视频| 国产亚洲av高清不卡| 国产一区二区三区在线臀色熟女| 淫妇啪啪啪对白视频| 日本与韩国留学比较| 嫁个100分男人电影在线观看| 黑人操中国人逼视频| 色综合亚洲欧美另类图片| 亚洲人成伊人成综合网2020| 国产精品久久久久久久电影 | 欧美3d第一页| 久久久久久人人人人人| 一边摸一边抽搐一进一小说| 不卡av一区二区三区| 国产欧美日韩精品亚洲av| 人妻久久中文字幕网| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放| 国产探花在线观看一区二区| 亚洲国产看品久久| 国产日本99.免费观看| 亚洲成人免费电影在线观看| av在线天堂中文字幕| 久久午夜综合久久蜜桃| 午夜福利欧美成人| 97超级碰碰碰精品色视频在线观看| 国产视频一区二区在线看| 天天添夜夜摸| 看片在线看免费视频| 国产伦精品一区二区三区视频9 | 午夜久久久久精精品| 男女午夜视频在线观看| 国产成人aa在线观看| 国产乱人伦免费视频| 琪琪午夜伦伦电影理论片6080| 国内久久婷婷六月综合欲色啪| www.自偷自拍.com| 香蕉国产在线看| 天天一区二区日本电影三级| 99国产综合亚洲精品| 久久久精品大字幕| 免费av毛片视频| 国产视频内射| 国产三级中文精品| 欧美日韩一级在线毛片| 不卡一级毛片| 桃红色精品国产亚洲av| 天天一区二区日本电影三级| 99国产综合亚洲精品| 亚洲在线观看片| 99re在线观看精品视频| 欧美中文日本在线观看视频| 国产三级中文精品| 久久性视频一级片| 日本成人三级电影网站| a级毛片在线看网站| 国产午夜精品论理片| 亚洲熟妇中文字幕五十中出| 熟女少妇亚洲综合色aaa.| 日韩精品青青久久久久久| www.www免费av| 免费看光身美女| 色综合站精品国产| 成人欧美大片| 久久精品影院6| 黄色丝袜av网址大全| 亚洲国产精品sss在线观看| 欧美日韩一级在线毛片| 三级国产精品欧美在线观看 | 大型黄色视频在线免费观看| 欧美色视频一区免费| 国产精品一区二区精品视频观看| а√天堂www在线а√下载| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜一区二区| 日韩欧美在线二视频| 亚洲成人精品中文字幕电影| 桃色一区二区三区在线观看| 亚洲 国产 在线| 午夜福利视频1000在线观看| 国产精品98久久久久久宅男小说| 亚洲欧美日韩卡通动漫| 天天一区二区日本电影三级| 男女之事视频高清在线观看| 熟妇人妻久久中文字幕3abv| 免费观看人在逋| 中文字幕av在线有码专区| 中文字幕人妻丝袜一区二区| 精品一区二区三区av网在线观看| 亚洲,欧美精品.| 欧洲精品卡2卡3卡4卡5卡区| 中文字幕人成人乱码亚洲影| 亚洲精品456在线播放app | 久久久久精品国产欧美久久久| 搡老岳熟女国产| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 亚洲五月天丁香| 亚洲欧美日韩东京热| or卡值多少钱| 成人精品一区二区免费| 欧美日韩中文字幕国产精品一区二区三区| 99热精品在线国产| 九色成人免费人妻av| 久久国产精品人妻蜜桃| 国产精品一及| 级片在线观看| 午夜免费观看网址| 国产aⅴ精品一区二区三区波| 久久精品aⅴ一区二区三区四区| 少妇的丰满在线观看| 此物有八面人人有两片| 午夜激情福利司机影院| 97碰自拍视频| 亚洲欧美激情综合另类| 国产三级在线视频| 久久久国产成人精品二区| 午夜福利在线在线| e午夜精品久久久久久久| 伊人久久大香线蕉亚洲五| 久9热在线精品视频| 一区二区三区激情视频| 国产成人影院久久av| 免费av不卡在线播放| 69av精品久久久久久| 国产亚洲精品一区二区www| 一个人免费在线观看的高清视频| 亚洲成av人片在线播放无| 琪琪午夜伦伦电影理论片6080| 免费观看精品视频网站| 国产精品九九99| 亚洲av成人精品一区久久| 宅男免费午夜| 国产高清有码在线观看视频| 搞女人的毛片| 亚洲av第一区精品v没综合| 亚洲一区二区三区不卡视频| 国产黄色小视频在线观看| 久久精品国产综合久久久| 九九热线精品视视频播放| 国产伦精品一区二区三区四那| 老汉色av国产亚洲站长工具| 美女被艹到高潮喷水动态| 久久精品国产清高在天天线| 国产人伦9x9x在线观看| 黄片小视频在线播放| 国内久久婷婷六月综合欲色啪| 国内揄拍国产精品人妻在线| 国产精品 欧美亚洲| 色噜噜av男人的天堂激情| 老司机福利观看| 精华霜和精华液先用哪个| 欧美丝袜亚洲另类 | 午夜福利成人在线免费观看| 搡老岳熟女国产| 成年女人毛片免费观看观看9| 精品久久久久久,| 欧美丝袜亚洲另类 | 制服人妻中文乱码| 天堂动漫精品| 一级毛片精品| or卡值多少钱| 欧美绝顶高潮抽搐喷水| 岛国在线免费视频观看| 国产精品一区二区精品视频观看| 丝袜人妻中文字幕| 日本黄色视频三级网站网址| 欧美成狂野欧美在线观看| 白带黄色成豆腐渣| 午夜福利免费观看在线| 成人性生交大片免费视频hd| 精品熟女少妇八av免费久了| 精品不卡国产一区二区三区| cao死你这个sao货| 午夜激情欧美在线| 国产精品亚洲美女久久久| 亚洲国产欧洲综合997久久,| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 中亚洲国语对白在线视频| 久久久色成人| 夜夜看夜夜爽夜夜摸| 亚洲精品久久国产高清桃花| 激情在线观看视频在线高清| 俄罗斯特黄特色一大片| 黄色丝袜av网址大全| 日韩欧美精品v在线| a级毛片在线看网站| 欧美激情久久久久久爽电影| 国产三级黄色录像| 久久久久国产一级毛片高清牌| 999久久久精品免费观看国产| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 波多野结衣高清作品| 特大巨黑吊av在线直播| 男女午夜视频在线观看| 国产亚洲av嫩草精品影院| 婷婷亚洲欧美| 国产伦人伦偷精品视频| 天堂av国产一区二区熟女人妻| 久久亚洲精品不卡| 热99re8久久精品国产| 国产精品一区二区精品视频观看| 久久精品综合一区二区三区| 两人在一起打扑克的视频| 一区二区三区激情视频| 制服丝袜大香蕉在线| 亚洲国产精品久久男人天堂| 日韩免费av在线播放| 亚洲欧美日韩卡通动漫| av黄色大香蕉| 村上凉子中文字幕在线| 这个男人来自地球电影免费观看| 国产亚洲精品av在线| 色吧在线观看| 免费在线观看成人毛片| 不卡一级毛片| 国模一区二区三区四区视频 | 国产精品日韩av在线免费观看| 少妇的逼水好多| 免费大片18禁| 一二三四在线观看免费中文在| 午夜免费激情av| 12—13女人毛片做爰片一| 狂野欧美白嫩少妇大欣赏| 人人妻人人澡欧美一区二区| 国产精品久久电影中文字幕| 亚洲精品色激情综合| 日日摸夜夜添夜夜添小说| 精品久久久久久久毛片微露脸| 国产免费男女视频| 在线观看午夜福利视频| 日韩免费av在线播放| 国产亚洲欧美98| 巨乳人妻的诱惑在线观看| 国内久久婷婷六月综合欲色啪| 天堂动漫精品| 欧美乱码精品一区二区三区| 久久久久久久精品吃奶| 999久久久精品免费观看国产| 成年人黄色毛片网站| 在线观看免费视频日本深夜| 国产成人精品久久二区二区91| 亚洲人成网站高清观看| 三级国产精品欧美在线观看 | 老司机在亚洲福利影院| 在线观看66精品国产| 狠狠狠狠99中文字幕| 9191精品国产免费久久| 99精品欧美一区二区三区四区| 国产精品永久免费网站| 在线观看一区二区三区| 国产97色在线日韩免费| 婷婷精品国产亚洲av在线| 在线观看免费午夜福利视频| 人人妻人人看人人澡| 精品福利观看| 91在线精品国自产拍蜜月 | 国产视频内射| 神马国产精品三级电影在线观看| 操出白浆在线播放| 亚洲精品在线美女| 狠狠狠狠99中文字幕| 啦啦啦观看免费观看视频高清| 国产又色又爽无遮挡免费看| 欧美中文日本在线观看视频| 嫩草影视91久久| 午夜福利视频1000在线观看| 白带黄色成豆腐渣| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 久久久久久久久久黄片| 丰满人妻熟妇乱又伦精品不卡| 精华霜和精华液先用哪个| 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 狠狠狠狠99中文字幕| 国产亚洲精品一区二区www| 久久久水蜜桃国产精品网| 久久伊人香网站| 国产不卡一卡二| 午夜精品在线福利| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 又黄又粗又硬又大视频| 成在线人永久免费视频| 国产av麻豆久久久久久久| 国产成年人精品一区二区| 校园春色视频在线观看| 美女免费视频网站| 国产极品精品免费视频能看的| 免费看十八禁软件| 欧美一区二区精品小视频在线| 午夜精品久久久久久毛片777| av中文乱码字幕在线| 国产97色在线日韩免费| 免费电影在线观看免费观看| 黑人操中国人逼视频| 亚洲黑人精品在线| 国产精品综合久久久久久久免费| 欧美成人免费av一区二区三区| 丰满人妻一区二区三区视频av | 亚洲自偷自拍图片 自拍| 丝袜人妻中文字幕| 999久久久国产精品视频| 成人鲁丝片一二三区免费| 午夜免费观看网址| 亚洲国产日韩欧美精品在线观看 | 高清在线国产一区| 欧美一区二区国产精品久久精品| 色老头精品视频在线观看| 欧美色视频一区免费| 成人三级黄色视频| 午夜精品一区二区三区免费看| 少妇裸体淫交视频免费看高清| 又紧又爽又黄一区二区| 无人区码免费观看不卡| 男女做爰动态图高潮gif福利片| 丁香六月欧美| 亚洲黑人精品在线| 日韩欧美 国产精品| 亚洲精品456在线播放app | 国产熟女xx| 精品久久久久久久末码| 久久精品人妻少妇| 高清在线国产一区| 九九热线精品视视频播放| 国产精品乱码一区二三区的特点| 亚洲欧美日韩高清专用| www日本黄色视频网| 免费在线观看日本一区| 欧美激情在线99| 国产高清视频在线播放一区| 99久久综合精品五月天人人| 亚洲国产欧洲综合997久久,| 757午夜福利合集在线观看| 久9热在线精品视频| 97人妻精品一区二区三区麻豆| 午夜久久久久精精品| 极品教师在线免费播放| 欧美三级亚洲精品| 狂野欧美激情性xxxx| 在线免费观看的www视频| 欧美乱色亚洲激情| tocl精华| 少妇丰满av| 99国产精品99久久久久| 一边摸一边抽搐一进一小说| 久久欧美精品欧美久久欧美| 在线播放国产精品三级| 琪琪午夜伦伦电影理论片6080| 韩国av一区二区三区四区| 亚洲av熟女| 亚洲av五月六月丁香网| 啪啪无遮挡十八禁网站| 九九在线视频观看精品| 最近在线观看免费完整版| 日本五十路高清| 99久久国产精品久久久| 日本黄大片高清| 在线观看日韩欧美| 丁香欧美五月| 国产精品 欧美亚洲| 热99在线观看视频| 婷婷丁香在线五月| 国产精品久久久久久精品电影| 久久午夜综合久久蜜桃| 波多野结衣高清作品| 一区二区三区国产精品乱码| 噜噜噜噜噜久久久久久91| 国产成+人综合+亚洲专区| 久久久久久久精品吃奶| 午夜两性在线视频| 在线观看日韩欧美| 日本与韩国留学比较| 国产精品久久电影中文字幕| 美女黄网站色视频| 国产激情欧美一区二区| 听说在线观看完整版免费高清| 又黄又粗又硬又大视频| 嫩草影院入口| 国内精品美女久久久久久| 美女 人体艺术 gogo| 美女免费视频网站| 中文在线观看免费www的网站| 亚洲熟妇中文字幕五十中出| 亚洲va日本ⅴa欧美va伊人久久| 久99久视频精品免费| 欧美午夜高清在线| 日本熟妇午夜| 午夜激情欧美在线| 久久久久久久久久黄片| 国产毛片a区久久久久| 狂野欧美白嫩少妇大欣赏| 97超视频在线观看视频| 欧美日韩瑟瑟在线播放| 精品无人区乱码1区二区| 88av欧美| 欧美日韩精品网址| 成年免费大片在线观看| 亚洲精品美女久久久久99蜜臀| 一个人看的www免费观看视频| 国产成人av教育| 岛国视频午夜一区免费看| 成人国产综合亚洲| 99久久精品国产亚洲精品| www.999成人在线观看| 黄色日韩在线| www国产在线视频色| 国产一区二区在线av高清观看| 99视频精品全部免费 在线 | 欧美黄色片欧美黄色片| 白带黄色成豆腐渣| 国产69精品久久久久777片 |