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

    An Adaptable Walking-skid for Seabed ROV under Strong Current Disturbance

    2014-04-24 02:00:25

    School of Marine Science and Technology,Newcastle University, Newcastle upon Tyne, NE1 7RU, United Kingdom

    An Adaptable Walking-skid for Seabed ROV under Strong Current Disturbance

    Jianting Si and Chengsiong Chin*

    School of Marine Science and Technology,Newcastle University, Newcastle upon Tyne, NE1 7RU, United Kingdom

    This paper proposed a new concept of an adaptable multi-legged skid design for retro-fitting to a remotely-operated vehicle (ROV) during high tidal current underwater pipeline inspection. The sole reliance on propeller-driven propulsion for ROV is replaced with a proposed low cost biomimetic solution in the form of an attachable hexapod walking skid. The advantage of this adaptable walking skid is the high stability in positioning and endurances to strong current on the seabed environment. The computer simulation flow studies using Solidworks Flow Simulation shown that the skid attachment in different compensation postures caused at least four times increase in overall drag, and negative lift forces on the seabed ROV to achieve a better maneuvering and station keeping under the high current condition (from 0.5 m/s to 5.0 m/s). A graphical user interface is designed to interact with the user during robot-in-the-loop testing and kinematics simulation in the pool.

    ROV;hexapod; multi-legged; skid; seabed; tidal current;

    1 Introduction

    Most underwater robots (Corradiniet al., 2011) use propellers for their propulsion and its theoretical mechanism is well-established. However, there exist a few difficulties to conduct pipeline inspection work because of disturbances caused by the high tidal current. Flow from the propellers causes disturbance on the seabed, stirring up soft sediment that impairs visibility. Efficiency of seabed operation is compromised as the operator will have to depend on the ocean current to carry away or clear the suspended particles. This effect magnifies when thrusters are actively used for discrete adjustment of the ROV (Chinet al.,2011; Enget al.,2008) or to provide downward thrust for a positively buoyant vehicle on a soft seabed to achieve desired position. Effects of such propeller action are also detrimental to the marine environment and undesirable to activities such as underwater archaeology as it disturbs the natural state of the subject to be observed (Rife and Rock, 2006; Tsusakaet al., 1986). In addition, shallow waters usually present a plethora of marine growth such as corals and seaweeds, which may be drawn into the suction stream of the thrusters and causefouling on mechanical propeller.

    Furthermore, ROVs have difficulties in getting precise positioning, accurate manipulation and clear acoustic image particularly in high current environment (Junet al., 2011). It is also not easy to achieve fine attitude control due to the highly non-linear dynamic characteristics caused by dead zone, delay and saturation. Although, the propellers have high efficiency when performing underwater flight, it is obviously ineffective in partially submerged or out of water condition. For the tasks that require amphibious capabilities in shallow water, reliance on propellers as the sole mode of propulsion will severely limit the ROV’s range of effectiveness. Besides, thrusters are expensive and have limited service life due to the seals and oil compensation system for watertight integrity.

    In addition to propeller-based propulsion, other methods of ROV locomotion could be found. A tracked-based category is a subset of ROVs commonly known as crawlers. They move on seabed via belted threads, or in combination with the thrusters. The belted threads have a large contact surface area with contact on the seabed to provide sufficient traction for benthic maneuvers for ROV (Inoueet al., 2013; Herzoget al., 2007). However, it caused stirring up of large clouds of sediment, sinkage due to thixotropic behavior of seabed. A wheel-based locomotion is another type of crawlers (Waldmann and Bergenthal, 2010; Nowaket al., 2008) developed to perform underwater intervention by direct contact with the seabed. A new type of compound-mobile mechanism (Yuet al., 2010) achieves amphibious ability by a wheel-propeller-leg driving method was proposed. It is apparent that characteristics of wheel-based locomotion are very much similar to the track-based. Wheels have less direct contact area with the terrain that minimizes the impact of disturbance. But the wheels’ ability to tackle irregularity may be slightly inferior to the tracked design. Compound-driven mechanisms may be novel and effective, but combination of different mechanisms can be quite difficult.

    In addition, a body undulation ROVs based on different modes of locomotion and hyper-redundant mechanisms were also used. Example of this type of vehicles are the AmphiBot I and II (Crespiet al., 2005; Yuet al., 2011; Zuoet al., 2009), which are both snake-like and achieves motion using linked non-linear oscillators to produce travelling waves both on land and in open water (McIsaacet al.,2003).Besides, various swimming modes of fish are enumerated and the analytical method for their propulsive mechanism was reviewed (Sfakiotakiset al., 1999). The advantage of these robots is their dexterity in a highly unstructured environment that makes them ideal for search and rescue operations. The modularity of such robots makes reconfiguration easy and mitigates the risk of the vehicle loss due to damage such as water ingress. However, due to the inherent shape and construction of these robots, they have limited payload capacity and are incapable of tool utilization other than basic camera and lights for search and observation operations. Also, the poor power efficiency and the requirement of complex kinematics and dynamics result in undulating locomotion.

    Recently, a multiple-legged underwater robot (Junet al., 2012) was then designed for high tidal current. The kinematic structure and gait of lobster were analyzed but the studies on the hydrodynamic on a ROV platform during testing were not performed. This spurred a further development of an amphibious six-legged underwater robot (Tanakaet al., 2004; Christinaet al., 2009; Theberge and Dudek, 2006)that they focused on the watertight technologies for the robot in underwater. Unfortunately, the six-legged design does not allow adaptability to other robots platform. The graphical-user-interface (GUI) for simulation and experimental testing of the seabed ROV are not available in the robot-in-the-loop testing. To best of our knowledge, the design and robot-in-the-loop simulation of an adaptable hexapod walking-skid for a seabed ROV has not been shown in a single publication.

    Hence, this paper will develop a low cost biomimetic solution using an attachable hexapod walking skid on a ROV platform. The biomimetic lobster-like seabed ROV is designed to improve the maneuvering and station keeping on seabed under turbulent environment. A simplified kinematic model to develop a ROV’s simulation and prototype control program using the GUI during the robot-in-the-loop testing will be presented. The drag and lift forces (and its coefficients) acting on the seabed ROV’s body and legs compensation postures will be determined using the computer software under the sea current disturbances (from 0.5 m/s to 5.0 m/s).

    This paper is organized as follows: The concept of the adaptable hexapod walking skid is addressed in Section 2. It is followed by kinematic modelling and simulation in Section 3. In Section 4 and 5, the computer simulation and experimental tests on the seabed ROV are presented respectively. Lastly, the conclusions are drawn.

    2 Design of adaptable hexapod walking skid

    Generally, legged locomotion refers to the conventional way of walking on the seabed or in the form of paddling of legs to result in vehicular displacement. Taking inspiration from nature, these bio-mimetic robots strive to emulate nature’s using locomotion to achieve proficiency on the seabed. In order for the body of the vehicle to be entirely decoupled from the irregularities of the seabed environment, the legs should be independently controlled in three degrees of freedom (DOF). A six-legged walking machine (hexapod) (Shimet al., 2013) was used for modeling, based on the criteria of a statically stable means of locomotion, and simplicity. Bottom dwelling marine creatures such as lobsters and crabs have evolved around the conditions of heavy turbulence caused by tides, current and wave action. Lobsters for example use their abdomen and appendages to create hydrodynamic adaptable shapes to cope with directional high velocity water flow (Maude and Williams, 1983). The resultant hydrodynamic translational forces allow them to achieve stability on the seabed even with a buoyant body and avoid negative lift that will sweep those away (Martinez, 2001). In addition, the hydrodynamic thrust achieved by body posture adjustment pushes their legs further into the substrate, increasing traction (Junet al., 2012) without the need to adjust body weight.

    To achieve the proposed legged form of locomotion on the seabed ROV, the attachable hexapod walking skid will be fitted onto existing ROV that has thrusters as its main propulsion. As such, the scope of the conventional propeller-driven ROV can be extended with the added range of motion to operate efficiently on the seabed and achieve amphibious ability to travel from deep to shallow water. The attachable hexapod walking skid as seen in Fig. 1 will be designed and prototyped in a small scale and constructed using readily available components.

    The mechanical system design for the seabed ROV can be seen in Fig. 1 (with no thrusters installed). To support development of the seabed ROV, a test-bed vehicle (weighted 3.5kg in air with a dimension of approximately 0.25 m×0.25 m×0.25 m at its full extended position) dedicated for the test-bedding of the attachable skid and control subsystems was developed by the Newcastle University. The seabed ROV is positively buoyant by design. The robot enclosure design is pressure tested at 50m. The maximum walking speed is not more than 1m/s. For initial prototype design, the payload (i.e. weight carrying capacity) should not be more than 2kg due to the servomotors used and the dimension of the prototype. The attachable skid will be attached to a small-size ROV platform with minimal modification.

    The ROV is equipped with 7.4 V 5.2 Ah lithium polymer battery power supply for up to 50 m for testing and initial deployment purpose. The key components are servomotors and electronics casing containing electronics circuit boards. The electronics casing was used to represent a solely propeller-driven ROV in which the hexapod walking skid will be attached. Four hexapod walking skids (or legs) are shown in Fig.1. Each leg consists of three servomotors, brackets and metal parts that come into contact with the ground. The two remaining legs were used as manipulator that consists of five servomotors, brackets and a gripper mechanism. The completed mechanical subsystems of theadaptable hexapod walking skid can be seen in Fig. 2. Besides, the mechanical subsystems, the buoyancy of the seabed ROV need to be adjusted. Due to the considerable volume of the casing which represents the ROV that is mounted on the skid, the entire prototype is positively buoyant as it displaced more water than its own weight.

    A microcontroller was used to control the actuators for moving the attachable hexapod walking skid. Four Seeeduino Mega microcontrollers were used in the project. The seabed ROV is connected to the computer via a RS232 communication. Each microcontroller takes 5 V 500 mA power from the lithium-polymer battery. The microcontrollers control the 22 servomotors (each sealed with silicone grease, Nitrile rubber o-ring and adhesive) attached to its legs and manipulators of the seabed ROV. The robot-in-the-loop testing and robot’s kinematic simulation is accomplished through the MATLABTMGUI software.

    Fig. 1 3D models of attachable walking skid fixed on a ROV (right)

    Fig. 2 Actual fabricated walking skid (left) and seabed ROV prototype (center & right)

    3 Kinematics of hexapod walking skid

    There are many methods used in robotics to describe the position and orientation of robotic links in Euclidean space(Zhu, 2001; Junet al., 2011). The approach taken in the kinematic derivation of the hexapod walking skid in this project differs slightly and undertakes the use of homogenous coordinates and application of corresponding homogenous transformation matrices in relation to a single absolute reference point, namely the point of origin. The four legs and two manipulators of the hexapod walking skid are each simplified into two-link three degree of freedom (DoF) linkages (see Fig. 3). The manipulators consist of five DoF, the remaining two DoF (the translational opening/closing of the gripper and the rotation of the gripper’s wrist) are omitted to reduce the complexity of the simulation.

    For brevity, only the derivation of one appendage, the front left manipulator (see Fig. 4) is shown as the rest are similar.

    Fig. 3 Overview definition of the coordinate system of the hexapod walking skid

    Fig. 4 Top view of the front left manipulator’s coordinate system and notation

    The point vectors are named with capital letters, which in the case of the front left manipulator – the capital letter ‘A’, and the lengths of legs or manipulator links with the capital letter ‘L’. Joint angles are preceded with the particular vector notation of interest, followed by ‘J’ and then the corresponding angle notation (See Fig.5). CP stands for controlling point that indicates the location of the hexapod skid in three-dimensional (3D) space. The rest of the coordinates of the model are then calculated based on CP. The length and width of the skid are represented by BL and BW, respectively. The lengths of the two link manipulator are assigned asL1andL2, respectively in Fig.5, with one end ofL2in contact with the ground and the other end ofL1to the skid.L3AandL4Aare virtual lengths, used to calculate the joint angles (AJa,AJb,AJcandAJd) betweenL1,L2, ground and skid, respectively.

    Fig. 5 Detailing of the front left manipulator’s coordinate system and notation

    With the knowledge ofCP,BL,BW, the vectors for all of the left manipulator can be found as follows: where OFF represents the linear distance ofA3from the skid’s edge in length and breadth-wise.

    With specification of the skid’s locating pointCPby (1), the ends of the manipulator links can be located by (2) and (3). The required angles for determination of the joint vectorA2between the two manipulators links are determined as follows:

    As shown in (1) to (10), specification of CP takes precedence before the determination ofA1andA3, which ultimately leads to calculation ofA2. The equations are applicable to all the appendages in the hexapod skid. For convenience of handling the equations, vectors representing the connection of the appendages’ link to the side of skid and vectors represent the other end of the appendages’ link to the ground are combined as follows:

    For manipulating the vectors in 3D space, the canonical transformation matrices are used. The below shown the steps applied in MATLABTMto obtain the desired control of the hexapod walking skid. For example, consider having the hexapod walking skid start at the coordinates of (0,0,50). Using (1), the CP corresponds to:

    Applying (2) to all the six legs of the skids from A1 to F1, with BL = 175, BW = 105, L1 = 63 and L2 = 77, the Body matrix is achieved with the following value:

    To achieve a desired elevation of 20 units in the direction of the Z-axis, the following translation operation is used.

    To simplify the mathematical complexities when performing arbitrary rotations that involve various reference frames and orientations, the following is used. If the location of the skid is not at the point of origin and a rotation about an arbitrary axis is required, CP (which represents the location of the skid) is first translated to the point of origin, where the required matrix operations are performed with the canonical rotation matrices. CP is then translated back to its initial position to achieve a similar result as a rotation about an arbitrary axis.

    whereToriis translation back to origin,Tiniis translation back to original position before rotation andRxis the rotation matrix about x-axis. The above equations illustrates the basic kinematics for positioning the components such as the‘spine’, attached legs and manipulators of the hexapod walking skid in 3D space.

    4 Computational fluid dynamic (CFD) simulation tests

    The SolidWorks Flow Simulation was used to simulate the seabed ROV under the strong current speed. It solves the basic equations that govern the fluid motion using the Finite Volume (FV) method. SolidWorks Flow Simulation solves the full Navier-Stokes equations while also applying fluid state equations that define the nature of the fluid, and also by other empirical laws for the dependency of viscosity and thermal conductivity. SolidWorks Flow Simulation begins by analyzing the model geometry and the computational domain. Boundary conditions such as fluid properties, pressure, and temperature must be given to the computational domain before solving for the solution. However, the smaller the mesh size implies greater computational domain or processing power and time to obtain the flow solution. As such, compromises need to be made during the setting up of the boundary conditions and the mesh size. The initial conditions are based on the flow study in 10 m depth of water measuring at 20?C with flow velocity of 1 m/s. The input temperature and pressure are set to 293 K and 199 249 Pa, respectively. The turbulence parameters are automatically calculated by the software and its recommended values are used. The ideal mesh level was determined for every individual flow study. For example, instead of computing the drag force at different mesh levels over the velocity range of 0 to 5 m/s, intermediate values of 1.0 m/s, 3.0 m/s and 5.0 m/s were used at incremental mesh levels. The deviation of values between each level was calculated and compared when the deviation drops below ±5%. This implies that the flow studies at higher mesh level have no significant changes to the simulation results besides having a longer computation time. As such, the full set of velocities from lowest to highest (0.5 m/s to 5.0 m/s) were computed at the mesh level.

    The computational domain is defined as the geometric area around the ROV. To obtain a complete flow study for analyzing the external flow characteristics, the size of the computational domain is shown in Fig. 6. The external flow boundary was set to five times the size of the ROV. In total, around 750 000 tetrahedral elements were built. The region of disturbed flow (usually turbulent) occurred at downstream of the ROV body. This region creates high velocity (or low pressure) region at the rear end, and thus low velocity (or high pressure) in the front to resist the ROV motion.

    Fig. 6 Computational domain for seabed ROV

    SolidWorks Flow Simulation runs the flow studies to derive the forces over the external surface of the 3D model. The aim of the sequence simulations is to determine the magnitude of the additional drag and lift force (and it corresponding drag and lift coefficients) when the hexapod walking skid is attached to the seabed ROV during maneuvering and on the seabed (with/without skid and at different seabed ROV’s body postures).

    4.1 Flow Simulation during underwater maneuvering

    There are many posture that the legs and manipulator of the hexapod walking skid can be positioned during maneuvering. The two main configurations are considered here. Configuration 1 (see Fig. 7 in the center) has the legs of the skid kept as close to the body of the ROV, with the aim of reducing frontal projection area against the current flow. By keeping the legs and manipulators retracted and close to the sides of the ROV, the servomotors at the joints do not have to strain as much against the flow of the current to maintain the position of the appendages. Configuration 2 (see Fig. 7 on the right) has the legs/manipulators of the hexapod walking skid spread out laterally to create a more streamlined vehicle, much like the shape of the fins of a stingray. Both configurations have similar frontal projected area, with Configuration 2’s area greater than Configuration 1 by approximately 8.6 %. Besides, the drag and lift forces in the flow direction were computed. The corresponding drag and lift coefficients (and Reynolds number) were determined using the highest flow velocities. The corresponding graphs were generated.

    (a) Model of the propeller-driven ROV only;

    (b) Model of the propeller-driven ROV with skid attached (Configuration 1);

    (c) Model of the propeller-driven ROV with skid attached in more compact posture (Configuration 2).

    Fig. 7 Underwater flight flow studies for without skid (left), configuration 1(center) and 2 (right)

    The derived values of drag force of underwater flight with/without the skid attached with both configurations are shown in Fig. 8. As observed, the attachment of the hexapod walking skid caused the overall drag of the vehicle to increase by approximately four times at the highest simulation stream velocity ?5 m/s. This is expected as the skid contributes to the higher frontal projected area. However, it is also observed that Configuration 1, which has a slightly smaller frontal projection area, actually incurred 14.3 % more drag than Configuration 2. This may be due to the free stream flow trajectories of the models, which the Configuration 1 presents a much more streamlined exterior with a lesser turbulent vortices around the vehicle than Configuration 2.

    Fig. 8 Drag force of ROV’s underwater flow studies

    4.2Flow simulation on seabed with body compensation postures

    The series of flow studies in this section attempts to determine the drag force under different body compensation postures (see Fig. 9) under the high tidal current on the seabed.

    (a) Model of the propeller-driven ROV with skid-attached in 10°compensation posture;

    (b)Model of the propeller-driven ROV with skid-attached in 20°compensation posture;

    (c)Model of the propeller-driven ROV with skid-attached in 30°compensation posture;

    (d) Model of the propeller-driven ROV with skid-attached in 40°compensation posture.

    The derived values of drag/lift force and coefficients of the model on the seabed at various compensation postures are shown in Fig. 10. As compared with the four compensation postures, it shows that the degree of tilt increases with reference to the flat seabed. The drag force on the model increases due to the negative lift on the model. The result implies that the tilting of the model in the current flow creates a downward force that is significant at compensation postures above 20°. Fig. 10 also included data for the model with the skid attached with legs and manipulators in default position. It can be seen that the incurred drag force on the model is less than the rest of the postures when compared with the default position (except at 40° posture). All the compensation postures results in much greater negative lift than the default position.

    In summary, the numerical flow studies indicated that the hexapod walking skid results in increased drag force, and a higher negative lift force (by the compensation postures) to adhere the ROV onto the seabed in the midst of a considerable current flow (from 0 to 5 m/s).

    Fig. 9 Compensation postures flow studies and actual implementation for configuration (a) to (d) - starting from top

    Fig. 10 Drag and lift coefficients compensation postures from 10 to 40 degrees

    Fig. 11 MATLABTMGraphics User Interface

    5 Experimental tests in air and water

    The movement of the hexapod walking skid was simulated using the GUI in MATLABTM. The object-oriented programming (OOP) design approach was taken in developing the MATLABTMscripts to support the GUI. The OOP allows modularization of the codes for easier management. It also allows concepts such as encapsulation, inheritance and polymorphism to be exploited in order to optimize both the development and execution of the program. As shown in Fig. 11, GUI was created to interact with the user during the dry and water testing. The user has a choice of performing simulation of the hexapod walking skid and/or real- time control of the actual prototype.

    The GUI provides access to methods such as manipulator direction controls together with planar and angular controls of the skid. There are three types of walking gaits (see Fig. 12), each with a different characteristic as follows. Metachronal wave gait is the slowest but most stable out of the three gaits, suited for traversing uneven terrains. Only one leg (usually from the rear) is raised, moved forward and lowered at single time before the next leg on the same side is moved in a similar manner. A five-point contact with the ground is made at all times that is more stable but slow. However, it has a good load carrying capability as the weight is distributed equally on the ground. Ripple gait is similar to the metachronal wave except the two legs are diagonally across each other are lifted to move forward. It is slightly less stable but faster than the metachronal wave gait. On the other hand, tripod gait is the fastest and least stable of the three gaits. It is suitable for traversing on flat and even terrains. The middle leg on one side of the body is lifted together with the front and rear legs of the other side to move forward, before the remaining legs in the first sequence are repeated with the same motions. Load carrying capability is weaker than the previous gaits as the load is distributed among the three legs only.

    Fig. 12 Screenshot of metachronal wave gait (left), ripple gait (center) and tripod gait (right)

    Fig. 13 depicts a block diagram of the MATLABTMGUI interface for the seabed ROV attached with a hexapod walking skid. After the user has selected the choice, the commands go into the serial communications on the four Seeeduino microcontrollers. Each microcontroller uses the arduino interactive development environment (IDE) to provide basic interface between the Arduino microcontroller, as well as to control the servomotors interactively from the MATLABTMcommand window, without having to write, debug, compile, upload and run the C programs. However, the Arduino IO can only control up to two servomotors connected to the microcontroller, whereas a total of 22 servomotors are required for the control of the hexapod walking skid. With that, the Arduino IO package code was modified and extended.

    The constructed seabed ROV prototype will be tested to verify the derived kinematic simulations and software/hardware interfacing. All tests were conducted with the help of the GUI. For clarity, only the GUI for the robot-in-the-loop testing is shown in Fig. 14. Initially, the stationary planar and angular motion were performed under the dry condition. The purpose is to verify basic coordination of the servomotors and identify any faulty servomotors or GUI interfacing error.

    During the pool test, it was observed that seabed ROV prototype and the GUI displayed some correlation. Besides, it was found that the coordination of servomotors is generally acceptable except the jerkiness when performing slow motion. The jerkiness was due to the analogue servomotors programmed to move in a fixed series incremental steps. This can be minimized by decreasing the step size, but it would affect the speed or rotation of the servomotor’s output shaft. With that, the digital servomotors were used and the parameters were programmed into the servomotor’s control board.

    Fig. 13 MATLABTMGUI for seabed ROV with attachable hexapod walking skid

    Fig. 14 Dry test showing correlation of both prototype and simulation on GUI

    Fig. 15 Test of prototype in the midst of yaw rotation underwater (left) and metachronal wave gait (right)

    In addition, the stationary planar, angular motion and metachronal wave gait (see Fig. 15) were conducted in a swimming pool. The prototype performed quite well in water and exhibit stability during the walking gait. The legs were able to achieve traction despite the slippery tiles on the swimming pool. Due to the deliberately slow and minimal accelerated motions, the underwater performance of the prototype barely differed from the one tested in dry condition.

    The testing of feasibility of compensation postures was performed. The ppol test verified the feasibility of implementing the compensation postures on the ROV. The prototype was able to achieve the various degrees of compensation postures (except for the 40?compensation posture). It was due to the fabrication fault that caused the top of the leg collided with the servomotor brackets of the hind legs. As a result, only a maximum of 35° compensation posture could be achieved. In summary, the waling skid on the seabed ROV is able to perform the required postures and maneuvering that give an increase on the drag and negative lift force required to sustain on the seabed under the high tidal current.

    6 Conclusions

    A new concept of an adaptable multi-legged skid design that could be retro-fitted to an existing remotely-operated vehicle (ROV) for high tidal current underwater pipeline inspection was proposed. The advantage of this concept is the ROV is able to have a higher stability during positioning on seabed than the one using propellers only. A simplified kinematic model and hydrodynamic simulation of the seabed ROV was presented. The computational-fluid dynamic (CFD) flow studies shown that the skid attachment provides an increase in overall drag and negative lift force of the vehicle at a maximum speed of 5m/s at certain body compensation postures. Hence, the hexapod walking skid at different posture provides a better maneuvering and station keeping on seabed under the sea current condition. With the graphics user interface using MATLABTM, it allows user to interact with the seabed robot during the robot-in-the-loop testing both in wet and dry condition. For the future works, the sensor such as gyroscope, altimeter, force sensors and accelerometer will be included in the seabed ROV for feedback control and comparisons with the simulated results. Lastly, actual implementation on the existing ROV (mounted with thrusters) for underwater pipeline inspection will be performed.

    Acknowledgment

    The authors would like to thank Newcastle University in United Kingdom (project account number:C0570D2330) and Ngee Ann Polytechnic (in Singapore) for providing laboratory space, and the reviewers for their valuable comments, which improve the readability of this paper significantly.

    Christina G, Meyer N, Martin B (2009). Simulation of an underwater hexapod robot.Ocean Engineering,36, 39-47.

    Chin C, Lau M, Low E (2011). Supervisory cascaded controllers design: experiment test on a remotely-operated vehicle.Proceedings of the IMechE Part C: Journal of Mechanical Engineering Science, 225(3), 584-603.

    Corradini ML, Monteriù A, Orlando G (2011). An actuator failure tolerant control scheme for an underwater remotely operated vehicle.IEEE Transactions on Control Systems Technology, 19(5), 1036-1046.

    Crespi A, Badertscher A, Guignard A, Ijspeert AJ (2005). AmphiBot I: an amphibious snake-like robot.Robotics and Autonomous Systems, 50, 163-175.

    Eng YH, Lau WS, Low E, Seet GL, Chin CS (2008). Estimation of the hydrodynamic coefficients of an ROV using free decay pendulum motion.Engineering Letters, 16(3), 326-331.

    Herzog K, Schulte E, Atmanand MA, Schwarz W (2007). Slip control system for a deep-sea mining machine.IEEE Transactions on Automation Science and Engineering, 4, 282-286.

    Inoue T, Shiosawa T, Takagi K (2013). Dynamic analysis of motion of crawler-type remotely operated vehicles.IEEE Journal of Oceanic Engineering, 38(2), 375-382.

    Jun BH, Lee PM, Baek H, Lim YK (2011). Approximated modeling of hydrodynamic forces acting on legs of underwater walking robot. IEEE OCEANS 2011, Spain, 1- 6.

    JunBH, ShimH, ParkJY, KimB, LeePM (2011). A new concept and technologies of multi-legged underwater robot for high tidal current environment.IEEE Symposium on Underwater Technology (UT) and Workshop on Scientific Use of Submarine Cables and Related Technologies (SSC), Tokyo, Japan, 475-479.

    Jun BH, Shim H, Kim B, Park JY, Baek H, Lee PM, Kim WJ, Park YS (2012). Preliminary design of the multi-legged underwater walking robot CR200. OCEANS 2012, Yeosu, South Korea, 1-4.

    Martinez MM (2001). Running in the surf: hydrodynamics of the shore crab Grapsus tenuicrustatus.Journal of Experimental Biology, 204, 3097-3112.

    Maude SH, Williams DD (1983). Behavior of crayfish in water currents: hydrodynamics of eight species with reference to their distribution patterns in southern Ontario.Canadian Journal of Fisheries and Aquatic Sciences, 40, 68-77.

    McIsaac KA, Ostrowski JP (2003). Motion planning for anguilliform locomotion. IEEE Transactions on Robotics and Automation, 19(4), 637-652.

    Nowak BM, Whitney T, Ackley SF (2008). Analysis of ROV video imagery for krill identification and counting under Antarctic sea ice.IEEE/OES Autonomous Underwater Vehicles, 1-9.

    Rife JH, Rock SM (2006). design and validation of a robotic control law for observation of deep-ocean Jellyfish.IEEE Transactions on Robotics, 22(2), 282-291.

    Sfakiotakis M, Lane, DM, Davies JBC (1999). Review of fish swimming modes for aquatic locomotion. IEEE Journal of Oceanic Engineering, 24, 2, 237-252.

    Shim H, Jun BH, Kang H, Yoo S, Lee GM, Lee PM (2013). Development of underwater robotic arm and leg for seabed robot, CRABSTER200.MTS/IEEE OCEANS, Bergen, Norway, 1-6.

    Tanaka T, Sakai H, Akizono J (2004). Design concept of a prototype amphibious walking robot for automated shore-line survey work. MTS/IEEE Techno-Ocean ′04, 834-839.

    Theberge M, Dudek G (2006).Gone swimmin’[seagoing robots].IEEE Spectrum, 43(6), 38-43.

    Tsusaka Y, Ishidera H, Itoh Y (1986). MURS-300 MK II: a remote inspection system for underwater facilities of hydraulic power plants.IEEE Journal of Oceanic Engineering,11(3), 358-363.

    Waldmann C, Bergenthal M(2010). CMOVE-a versatile underwater vehicle for seafloor studies.OCEANS 2010, Seattle, USA, 1-3.

    Yu JC, Tang YG, Zhang XQ, Liu CJ (2010). Design of a wheel-propeller-leg integrated amphibious robot.11th International Conference on Control Automation Robotics & Vision (ICARCV), 1815-1819.

    Yu S, Ma S, Li B, Wang Y (2011). An amphibious snake-like robot with terrestrial and aquatic gaits.IEEE International Conference on Robotics and Automation (ICRA), 2960-2961.

    Zhu WH (2005). On adaptive synchronization control of coordinated multirobots with flexible/rigid constraints.IEEE Transactions on Robotics, 21(3), 520-525.

    Zuo Z, Wang Z, Li B, Ma S (2009). Serpentine locomotion of a snake-like robot in water environment.IEEE International Conference on Robotics and Biomimetics, 25-30.

    Author biographies

    Jianting Siisthe assistant superintendent for technical services at POSH Fleet Services Pte Ltd. He has been working in maritime industry for few years. He worked in Lloyd’s Register UK for his internship and has a few years of experience as a marine engineer with Neptune Orient Lines (NOL) and American President Lines (APL) Technical Service Department. He has a degree of bachelor of Engineering with First Class Honours in Marine Engineering. He received the school prize for Best Performance in Bachelor of Engineering with Honours in Marine Engineering, RINA-Keppel Student Naval Architect award for the Best Final Year Project, SIT-Keppel Gold Medal, SNAMES award and IMAREST award. His research interests are in the robotics, systems integration and hardware-in-the-loop testing

    Chengsiong Chinisa senior lecturer at Newcastle University (Singapore). He received the Ph.D from Nanyang Technological University (NTU) in 2009 and a M.Sc. from Control System Centre in The University of Manchester in 2001 after completed his B.Eng. at NTU in 2000. He worked in industry for 6 years before moving into academia. He currently holds 3 U.S. Patents, 2 provisional US patent applications and 2 Trade Secrets. In 2013, he obtained a research grant from Singapore Maritime Institute (SMI) on High Performance Lithium-Ion Battery Power System for Long Endurance Deep Water Operation (worth S$0.59M) and EDB-IPP on noise and vibration prediction. He is the FHEA, FIMarEST, SIEEE, CEng, Eur.Ing and MIET. He is also a Member of the Editorial Advisory Board of The Mediterranean Journal of Measurement and Control. Recently, he becomes a member of the IMarEST′s Biofouling Management Expert Group. He had authored a book entitled Computer-aided Control Systems Design: Practical Applications using MATLABTMand SimulinkTMpublished by CRC Press, USA in 2012. His research interests are in the marine mechatronic systems design.

    1671-9433(2014)03-0306-10

    Received date:2014-02-03.

    Accepted date:2014-06-05.

    Foundation:Suuported by Newcastle University in United Kingdom (Project account number:C0570D2330).

    *Corresponding author Email: cheng.chin@ncl.ac.uk

    ? Harbin Engineering University and Springer-Verlag Berlin Heidelberg 2014

    精品人妻偷拍中文字幕| 中文字幕亚洲精品专区| 亚洲人与动物交配视频| 三级国产精品片| 国产亚洲一区二区精品| 国产精品蜜桃在线观看| 亚洲精品色激情综合| 亚洲精品乱久久久久久| 国产免费一区二区三区四区乱码| 精品人妻一区二区三区麻豆| 大陆偷拍与自拍| 看非洲黑人一级黄片| 久久综合国产亚洲精品| 亚洲精品456在线播放app| kizo精华| 欧美精品一区二区大全| 国产亚洲午夜精品一区二区久久| 高清不卡的av网站| 在线亚洲精品国产二区图片欧美| 久久精品人人爽人人爽视色| 99热网站在线观看| 美国免费a级毛片| 高清毛片免费看| 国产探花极品一区二区| 亚洲天堂av无毛| 黑人欧美特级aaaaaa片| 韩国高清视频一区二区三区| 日韩,欧美,国产一区二区三区| 久久99精品国语久久久| 亚洲,一卡二卡三卡| 在线免费观看不下载黄p国产| 男人添女人高潮全过程视频| 91午夜精品亚洲一区二区三区| 精品午夜福利在线看| 99久久综合免费| 久久久久久久国产电影| www.av在线官网国产| 亚洲熟女精品中文字幕| 精品少妇内射三级| 高清在线视频一区二区三区| 午夜福利在线观看免费完整高清在| 久久精品国产自在天天线| 亚洲精品久久午夜乱码| 亚洲美女搞黄在线观看| 插逼视频在线观看| 又粗又硬又长又爽又黄的视频| 99久久精品国产国产毛片| 涩涩av久久男人的天堂| 欧美xxⅹ黑人| 夫妻午夜视频| 国产精品一二三区在线看| 熟女av电影| 亚洲精品视频女| 国产精品欧美亚洲77777| 国产片特级美女逼逼视频| 国产日韩欧美在线精品| 亚洲精品乱码久久久久久按摩| 夫妻午夜视频| 青春草国产在线视频| 国产伦理片在线播放av一区| 免费人妻精品一区二区三区视频| 亚洲国产精品专区欧美| 亚洲欧洲国产日韩| 一级毛片电影观看| 9191精品国产免费久久| 国产在视频线精品| 免费看av在线观看网站| 五月天丁香电影| 麻豆乱淫一区二区| 亚洲国产成人一精品久久久| 国产xxxxx性猛交| 在线观看免费日韩欧美大片| 亚洲欧美成人综合另类久久久| 亚洲欧美日韩卡通动漫| 黑人巨大精品欧美一区二区蜜桃 | 午夜av观看不卡| 狂野欧美激情性xxxx在线观看| 秋霞在线观看毛片| 国产在线免费精品| 街头女战士在线观看网站| 久久久久精品久久久久真实原创| 少妇的丰满在线观看| 男女免费视频国产| 欧美xxⅹ黑人| 大香蕉久久成人网| 日韩av免费高清视频| 国产精品一区二区在线不卡| 国产免费一区二区三区四区乱码| 五月玫瑰六月丁香| 欧美日韩视频高清一区二区三区二| 男男h啪啪无遮挡| 人妻少妇偷人精品九色| 日韩制服丝袜自拍偷拍| 亚洲人成网站在线观看播放| 人体艺术视频欧美日本| 一区二区三区精品91| 五月玫瑰六月丁香| 亚洲色图综合在线观看| 97精品久久久久久久久久精品| 国产又爽黄色视频| 国产精品偷伦视频观看了| 久久人人爽av亚洲精品天堂| 女人被躁到高潮嗷嗷叫费观| 伦理电影大哥的女人| 精品一区二区三区四区五区乱码 | 97超碰精品成人国产| 91国产中文字幕| 国产成人av激情在线播放| 欧美日韩亚洲高清精品| 亚洲第一区二区三区不卡| 久久人人爽人人爽人人片va| 亚洲成人一二三区av| 在线天堂中文资源库| 日韩伦理黄色片| 国产精品久久久久久av不卡| 久久这里有精品视频免费| 亚洲av国产av综合av卡| 美女福利国产在线| h视频一区二区三区| 视频在线观看一区二区三区| 国产成人免费无遮挡视频| av线在线观看网站| 日韩人妻精品一区2区三区| 黑人欧美特级aaaaaa片| 丁香六月天网| 国产福利在线免费观看视频| 久久久久人妻精品一区果冻| 国产一级毛片在线| 色网站视频免费| 少妇人妻久久综合中文| 免费大片18禁| 日韩一区二区三区影片| 91国产中文字幕| 久久久久久伊人网av| av卡一久久| 亚洲精品一区蜜桃| 99国产综合亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| 国产成人av激情在线播放| 超碰97精品在线观看| a级片在线免费高清观看视频| 满18在线观看网站| 啦啦啦在线观看免费高清www| 老司机亚洲免费影院| 国产乱人偷精品视频| 看十八女毛片水多多多| 日韩制服骚丝袜av| 久久国产亚洲av麻豆专区| 久热久热在线精品观看| 久久精品国产亚洲av涩爱| 欧美bdsm另类| 国产免费视频播放在线视频| 999精品在线视频| 人妻 亚洲 视频| 黑丝袜美女国产一区| 波多野结衣一区麻豆| 免费av中文字幕在线| 最近手机中文字幕大全| 久久精品国产亚洲av天美| 国产探花极品一区二区| 高清黄色对白视频在线免费看| 成人影院久久| 新久久久久国产一级毛片| 国产免费一级a男人的天堂| 一本色道久久久久久精品综合| 日韩av免费高清视频| 寂寞人妻少妇视频99o| 五月伊人婷婷丁香| 国产日韩一区二区三区精品不卡| 国产精品久久久久久精品电影小说| 18禁动态无遮挡网站| 十八禁高潮呻吟视频| 少妇的逼水好多| 伊人亚洲综合成人网| 国产国语露脸激情在线看| 国产一区亚洲一区在线观看| 成人无遮挡网站| 黄色 视频免费看| 免费在线观看完整版高清| 最近最新中文字幕大全免费视频 | 大话2 男鬼变身卡| 国产成人91sexporn| 99九九在线精品视频| 高清毛片免费看| 国产免费现黄频在线看| 国产极品粉嫩免费观看在线| a级片在线免费高清观看视频| 精品国产一区二区三区久久久樱花| 男人操女人黄网站| 一级片'在线观看视频| 免费观看a级毛片全部| 青春草视频在线免费观看| 伦精品一区二区三区| 26uuu在线亚洲综合色| 久久 成人 亚洲| 亚洲综合精品二区| 亚洲伊人久久精品综合| 国产成人av激情在线播放| 纯流量卡能插随身wifi吗| 免费大片黄手机在线观看| 国产白丝娇喘喷水9色精品| 日韩在线高清观看一区二区三区| 纵有疾风起免费观看全集完整版| 少妇的丰满在线观看| 麻豆乱淫一区二区| 国产精品不卡视频一区二区| 美女大奶头黄色视频| 日本猛色少妇xxxxx猛交久久| 午夜91福利影院| 精品一区二区三区视频在线| 免费黄频网站在线观看国产| 90打野战视频偷拍视频| 国产成人免费观看mmmm| 午夜老司机福利剧场| 国产高清三级在线| 蜜臀久久99精品久久宅男| 欧美精品国产亚洲| 曰老女人黄片| 高清视频免费观看一区二区| 亚洲精品久久久久久婷婷小说| 久久97久久精品| 乱人伦中国视频| 九九在线视频观看精品| 精品久久久久久电影网| xxxhd国产人妻xxx| 新久久久久国产一级毛片| 天天操日日干夜夜撸| 亚洲情色 制服丝袜| 中文字幕免费在线视频6| 日本免费在线观看一区| 国产亚洲午夜精品一区二区久久| 欧美成人精品欧美一级黄| videosex国产| 日韩大片免费观看网站| av女优亚洲男人天堂| 赤兔流量卡办理| 免费观看av网站的网址| 国国产精品蜜臀av免费| 日韩一区二区三区影片| 在线天堂中文资源库| 菩萨蛮人人尽说江南好唐韦庄| 男男h啪啪无遮挡| 久久精品久久久久久久性| 国产一区有黄有色的免费视频| 99re6热这里在线精品视频| 边亲边吃奶的免费视频| 亚洲av成人精品一二三区| 成人漫画全彩无遮挡| 成年人午夜在线观看视频| 亚洲av在线观看美女高潮| 欧美+日韩+精品| 中文天堂在线官网| 亚洲伊人久久精品综合| 最近中文字幕高清免费大全6| 在线免费观看不下载黄p国产| 国产亚洲av片在线观看秒播厂| 看免费成人av毛片| 免费高清在线观看视频在线观看| 啦啦啦在线观看免费高清www| 午夜91福利影院| 中文字幕制服av| 亚洲性久久影院| 中文字幕人妻熟女乱码| 视频中文字幕在线观看| 亚洲精品456在线播放app| 亚洲中文av在线| 成年动漫av网址| 亚洲国产日韩一区二区| 国产熟女午夜一区二区三区| 观看美女的网站| 另类精品久久| 18禁国产床啪视频网站| 制服诱惑二区| 国产成人a∨麻豆精品| 免费黄频网站在线观看国产| 国产成人精品无人区| 黄片播放在线免费| 老女人水多毛片| 亚洲av欧美aⅴ国产| 人人澡人人妻人| 天天影视国产精品| 久久久久视频综合| 丝袜喷水一区| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 日本av免费视频播放| 日韩 亚洲 欧美在线| av.在线天堂| 国产精品久久久久久av不卡| 女性生殖器流出的白浆| 亚洲欧美日韩卡通动漫| 久久鲁丝午夜福利片| 国产片内射在线| 久久国产精品男人的天堂亚洲 | 侵犯人妻中文字幕一二三四区| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品| av有码第一页| 91精品国产国语对白视频| 久久久久久人人人人人| 国产黄频视频在线观看| 日韩一区二区三区影片| 国产精品国产三级专区第一集| 久久久久精品性色| 成人毛片a级毛片在线播放| 美国免费a级毛片| 女性生殖器流出的白浆| 黄色毛片三级朝国网站| 伦理电影大哥的女人| 观看美女的网站| 国产又色又爽无遮挡免| 一级毛片电影观看| 国产精品国产av在线观看| 亚洲国产欧美在线一区| 精品国产一区二区三区四区第35| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲精品第一综合不卡 | 久久婷婷青草| 午夜激情久久久久久久| 日韩电影二区| 九色成人免费人妻av| videos熟女内射| 成人国产麻豆网| 欧美 日韩 精品 国产| 欧美亚洲 丝袜 人妻 在线| 国产成人aa在线观看| 成人毛片60女人毛片免费| 男男h啪啪无遮挡| 一区二区日韩欧美中文字幕 | 久久婷婷青草| 亚洲人成77777在线视频| 十分钟在线观看高清视频www| 国产精品久久久久成人av| 久久这里有精品视频免费| 如日韩欧美国产精品一区二区三区| 国产成人免费无遮挡视频| 日本欧美视频一区| 免费人妻精品一区二区三区视频| 好男人视频免费观看在线| 精品人妻在线不人妻| 最新中文字幕久久久久| 亚洲国产色片| 岛国毛片在线播放| 毛片一级片免费看久久久久| 亚洲成人av在线免费| 啦啦啦中文免费视频观看日本| 99re6热这里在线精品视频| 国产精品女同一区二区软件| 人人妻人人爽人人添夜夜欢视频| 精品少妇久久久久久888优播| 97超碰精品成人国产| 国产亚洲最大av| 日韩欧美一区视频在线观看| 国产男人的电影天堂91| 秋霞伦理黄片| av天堂久久9| 97超碰精品成人国产| 狂野欧美激情性xxxx在线观看| 久久毛片免费看一区二区三区| 丰满乱子伦码专区| 宅男免费午夜| 在线观看美女被高潮喷水网站| 欧美老熟妇乱子伦牲交| 黄色怎么调成土黄色| 精品少妇黑人巨大在线播放| 一二三四中文在线观看免费高清| 成人免费观看视频高清| 最后的刺客免费高清国语| 成人亚洲欧美一区二区av| 精品国产国语对白av| 成人亚洲欧美一区二区av| 成人无遮挡网站| 国产av精品麻豆| 美国免费a级毛片| 在线天堂最新版资源| 最近手机中文字幕大全| 国产女主播在线喷水免费视频网站| 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区黑人 | 国产xxxxx性猛交| 大片电影免费在线观看免费| 人妻 亚洲 视频| 丰满乱子伦码专区| 人妻一区二区av| 国产综合精华液| 午夜福利,免费看| 免费大片18禁| 性高湖久久久久久久久免费观看| 午夜老司机福利剧场| 欧美日韩一区二区视频在线观看视频在线| 老司机影院成人| 亚洲av欧美aⅴ国产| 国产黄色视频一区二区在线观看| av网站免费在线观看视频| 97人妻天天添夜夜摸| 亚洲精品视频女| 建设人人有责人人尽责人人享有的| 91国产中文字幕| 婷婷成人精品国产| 一区二区av电影网| 欧美精品高潮呻吟av久久| 欧美精品国产亚洲| 国产成人午夜福利电影在线观看| 午夜免费鲁丝| 性高湖久久久久久久久免费观看| 午夜福利影视在线免费观看| 精品福利永久在线观看| 天堂俺去俺来也www色官网| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 欧美精品国产亚洲| kizo精华| 免费观看无遮挡的男女| 国产有黄有色有爽视频| 黄色视频在线播放观看不卡| 一区二区三区精品91| 黄色 视频免费看| 久久久久精品久久久久真实原创| 久久久a久久爽久久v久久| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| 久久99热这里只频精品6学生| 女人久久www免费人成看片| 亚洲高清免费不卡视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美中文字幕日韩二区| 日产精品乱码卡一卡2卡三| 超碰97精品在线观看| 精品国产乱码久久久久久小说| 久久精品国产鲁丝片午夜精品| 亚洲成国产人片在线观看| 青春草视频在线免费观看| 91久久精品国产一区二区三区| 亚洲成av片中文字幕在线观看 | 老司机影院成人| 久久久亚洲精品成人影院| 久久精品久久久久久久性| 26uuu在线亚洲综合色| 精品人妻一区二区三区麻豆| 日本猛色少妇xxxxx猛交久久| 亚洲欧美成人精品一区二区| 久久精品久久精品一区二区三区| 精品国产露脸久久av麻豆| 精品一区二区三区视频在线| 日本与韩国留学比较| 晚上一个人看的免费电影| 色婷婷av一区二区三区视频| 亚洲精品久久成人aⅴ小说| 日韩制服丝袜自拍偷拍| 免费观看在线日韩| 成人漫画全彩无遮挡| 久热这里只有精品99| 一本色道久久久久久精品综合| 18禁观看日本| 亚洲熟女精品中文字幕| 少妇人妻精品综合一区二区| 国产福利在线免费观看视频| 人妻人人澡人人爽人人| 国产精品.久久久| 两个人看的免费小视频| 精品一品国产午夜福利视频| 黑人高潮一二区| 精品人妻一区二区三区麻豆| 又大又黄又爽视频免费| 国产一级毛片在线| 国产精品偷伦视频观看了| 99热全是精品| 满18在线观看网站| 国产精品熟女久久久久浪| 一级,二级,三级黄色视频| 国产一级毛片在线| 成人国产av品久久久| 成人影院久久| 亚洲国产日韩一区二区| 国产免费视频播放在线视频| 高清黄色对白视频在线免费看| 男人爽女人下面视频在线观看| 汤姆久久久久久久影院中文字幕| 成年av动漫网址| 自拍欧美九色日韩亚洲蝌蚪91| 日韩成人伦理影院| av国产久精品久网站免费入址| 成人毛片a级毛片在线播放| 妹子高潮喷水视频| 我的女老师完整版在线观看| 国产精品人妻久久久久久| av国产精品久久久久影院| 人人妻人人澡人人爽人人夜夜| 亚洲精品日本国产第一区| 日本与韩国留学比较| 亚洲伊人色综图| 午夜久久久在线观看| 99国产精品免费福利视频| 一区二区三区乱码不卡18| 一级,二级,三级黄色视频| 女的被弄到高潮叫床怎么办| 国产男人的电影天堂91| 亚洲在久久综合| 国产在视频线精品| av在线观看视频网站免费| 亚洲av在线观看美女高潮| 国产男女内射视频| 精品人妻熟女毛片av久久网站| 夜夜骑夜夜射夜夜干| 国产精品国产三级专区第一集| 黑人巨大精品欧美一区二区蜜桃 | 日本黄色日本黄色录像| 一本大道久久a久久精品| 国产精品国产三级专区第一集| 黑人巨大精品欧美一区二区蜜桃 | 看免费成人av毛片| 在线 av 中文字幕| 69精品国产乱码久久久| 欧美97在线视频| 黄片播放在线免费| 你懂的网址亚洲精品在线观看| 免费黄色在线免费观看| 精品久久久久久电影网| 午夜免费男女啪啪视频观看| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 久久久久久久大尺度免费视频| 国产成人精品婷婷| 国产欧美亚洲国产| 亚洲精品日本国产第一区| 成人黄色视频免费在线看| 精品亚洲成国产av| 美女xxoo啪啪120秒动态图| 插逼视频在线观看| 性色av一级| 免费不卡的大黄色大毛片视频在线观看| kizo精华| 日韩免费高清中文字幕av| 大香蕉久久成人网| 丁香六月天网| 国产综合精华液| 精品熟女少妇av免费看| 亚洲精品456在线播放app| 大码成人一级视频| 在线天堂中文资源库| a级毛片黄视频| 婷婷色综合大香蕉| 精品国产一区二区三区四区第35| 亚洲精品av麻豆狂野| 高清黄色对白视频在线免费看| 欧美成人精品欧美一级黄| 2018国产大陆天天弄谢| 国产精品.久久久| 日韩 亚洲 欧美在线| 草草在线视频免费看| 午夜免费男女啪啪视频观看| 午夜91福利影院| 一区二区三区精品91| 国产高清三级在线| 看免费av毛片| 中文天堂在线官网| av线在线观看网站| 成人漫画全彩无遮挡| 免费观看av网站的网址| 国产精品久久久久久久久免| 各种免费的搞黄视频| 99国产精品免费福利视频| 国产老妇伦熟女老妇高清| 最新中文字幕久久久久| 国产又色又爽无遮挡免| 日韩免费高清中文字幕av| 国产精品无大码| 成人综合一区亚洲| 日日爽夜夜爽网站| 男女边摸边吃奶| 乱码一卡2卡4卡精品| 午夜免费鲁丝| av网站免费在线观看视频| 99久久人妻综合| 美女xxoo啪啪120秒动态图| 国产免费福利视频在线观看| 中文字幕免费在线视频6| 精品一区在线观看国产| 最近中文字幕2019免费版| 超碰97精品在线观看| 国产亚洲欧美精品永久| 国产一区二区激情短视频 | 免费看不卡的av| 91在线精品国自产拍蜜月| 成人亚洲欧美一区二区av| 亚洲欧美一区二区三区国产| 国产成人免费观看mmmm| 精品少妇黑人巨大在线播放| 免费日韩欧美在线观看| 日韩精品有码人妻一区| 亚洲欧美成人精品一区二区| 丰满少妇做爰视频| 男的添女的下面高潮视频| 亚洲成人一二三区av| 女人精品久久久久毛片| 久久久久精品人妻al黑| 国产 精品1| av免费观看日本| 免费黄频网站在线观看国产| 男女国产视频网站| 狠狠婷婷综合久久久久久88av| 亚洲国产成人一精品久久久| 国产亚洲最大av| 永久网站在线| 纯流量卡能插随身wifi吗| 中文字幕最新亚洲高清| 制服丝袜香蕉在线| 黑人巨大精品欧美一区二区蜜桃 | 九九在线视频观看精品| 超碰97精品在线观看| 国产亚洲午夜精品一区二区久久| 99国产综合亚洲精品| 女人被躁到高潮嗷嗷叫费观| 欧美精品av麻豆av| av.在线天堂|