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

    Classifying wakes produced by self-propelled fish-like swimmers using neural networks

    2020-07-01 05:13:16BinglinLiXiangZhangXingZhang

    Binglin Li, Xiang Zhang, Xing Zhang,*

    a State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China

    b School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China

    Keywords:Flow sensing Machine learning Wake classification Self-propelled swimming Fluid-structure interaction

    ABSTRACT We consider the classification of wake structures produced by self-propelled fish-like swimmers based on local measurements of flow variables. This problem is inspired by the extraordinary capability of animal swimmers in perceiving their hydrodynamic environments under dark condition. We train different neural networks to classify wake structures by using the streamwise velocity component, the crosswise velocity component, the vorticity and the combination of three flow variables, respectively. It is found that the neural networks trained using the two velocity components perform well in identifying the wake types, whereas the neural network trained using the vorticity suffers from a high rate of misclassification. When the neural network is trained using the combination of all three flow variables, a remarkably high accuracy in wake classification can be achieved. The results of this study can be helpful to the design of flow sensory systems in robotic underwater vehicles.

    The long-lasting wake structures left behind active swimmers and inanimate objects contain useful hydrodynamic cues about the sources that produce them. Such information can be extracted and exploited by other swimmers to accomplish a variety of tasks [1]. For example, with the help of the lateral line system, even blind fish can school [2] and avoid obstacles in the water [3, 4]. By using the whiskers, harbor seals can track the preys by detecting the wake signals [5]. This remarkable capability in animal swimmers also inspired the development of artificial flow sensory systems [6–10].

    The wake structures behind a bluff body and an active swimmer are often characterized by the von Karman and reverse von Karman vortex streets, respectively. Based on the number of vortices that are shed per flapping cycle, the wakes produced by an active swimmer can also be further classified into different subtypes. The “2S” and “2P” wakes are the two most commonly seen sub-types. Here “2S” refers to the situation in which two vortices of opposite sign are shed per cycle; whereas “2P” refers to the situation in which two pairs of vortices are shed per cycle. Besides, some higher-order exotic wakes, with up to 16 vortices shed per cycle have also been observed [11]. The types of wake produced by the swimmers are intimately associated with the swimming characteristics.

    For an animal (or artificial) swimmer, the capability of distinguishing between different wake types lays the groundwork for motion planning and navigation. However, this seemingly easy task is actually an inverse problem which poses great challenges to the flow sensory system. The reason is that the data of measurements are only available at discrete locations within a small sub-domain of the entire flow field. This is in contrast with the situation faced by an outside observer to whom the globally spacial-temporal flow patterns are accessible.

    With the rapid progress of machine learning, fundamental methodologies in machine learning have been used to handle the tasks of fluid mechanics [12]. Recently, some pioneering works have been conducted to solve this inverse problem by using machine learning algorithms. Wang and Hemati [13] generated “2S” and “2P” wakes with the point vortices by using the ideal-flow model. The wake-induced velocities were measured on the surface of a fish-like body that was immersed in the vortical wakes. The K-nearest neighbors (KNN) algorithm was then used to classify the wake types. Colvert et al. [14] generated the wakes behind a tethered rigid pitching foil placed in an otherwise steady stream by solving the Navier?Stokes equations numerically. The neural networks were then trained to distinguish between three distinct wake types from the time histories of vorticity measured at discrete locations. Alsalman et al. [15] used the same wakes as those generated by Colvert et al. [14], and went one step further to investigate how the accuracy of wake classification was affected by the types of sensors (i.e., sensors measuring different flow variables). They also trained the neural networks to classify wake types from the data of “one time”measurements at a spatially-distributed array of sensors. Motivated by the aforementioned studies, in this work, we train the neural networks to classify wakes behind a self-propelled fish-like swimmer. Similar to the studies of Refs. [14, 15], in the present study the wakes are also produced by high-fidelity numerical simulations. However, the wakes considered here and those in the two aforementioned works are very dissimilar, in terms of the spatial-temporal patterns. Clearly, the wakes considered here bear a much closer resemblance to those produced by an animal (or an artificial) swimmer under free-swimming condition.

    We consider continuous and intermittent self-propelled swimming of a flexible fish-like swimmer which is driven by fluid-structure interaction (FSI) [16, 17]. The localised actuation consists of the coupled heaving and pitching motions that are imposed at the leading edge (see Fig. 1). For the continuous swimming, the driving motions are:

    where ys=0and θs=0are the instantaneous heaving and pitching displacement of the leading edge; A and Aθare the heaving and pitching amplitudes, respectively; f is the actuation frequency,and ? is the phase lag between the heaving and pitching motions. In this study, the phase lag is set to 90°. For the intermittent swimming, certain time duration of passive coasting is interspersed between two half-periods of active bursting(flapping) [16]. The degree of intermittency is characterised by the duty cycle (DC), which is defined as the ratio of the time duration of active flapping to the time duration of one full cycle of burst-and-coast locomotion.

    The incompressible Naiver?Stokes equations and nonlinear dynamics equations constitute the governing equations of the system. These equations can be written in a dimensionless form as:

    where u is the fluid velocity; p is the pressure; X is the position vector for the Lagrangian (marker) point on the flexible swimmer and s is the Lagrangian coordinate along the arc; f and F are the dimensionless Eulerian and Lagrangian forces which represent the interaction between the flow and the structure. We select the chord length L and Uref=L f as the reference length and reference velocity for scaling the system.The four dimensionless quantities that arise in the governing equations are: flapping Reynolds number Ref, mass ratio β,dimensionless tension ξ and dimensionless rigidity γ. The definitions of these four dimensionless quantities as:

    Fig. 1. Schematic diagram of the wake classification problem: an array of sensors is placed in the wake behind a self-propelled fish-like swimmer.

    where ν is the kinematic viscosity of the fluid; ρfand ρsare the densities of the fluid and the filament, respectively; δ is the thickness of the filament; T and B are the tension and bending ri gidity of the structure, respectively. Here ξ is a Lagrangian multiplier which is determined to enforce the inextensibility condition (Eq. (6)) [18]. β and γ are two exponential functions with respect to s [16, 17]. In this study, R efis fixed to 6 000, which is the same as the highest Reynolds number considered in Ref.[16]. The Reynolds numbers based on the swimming velocity lie in the range of 3000 to 4000. The coefficients in the two aforementioned exponential functions are tunable.

    This FSI problem is simulated by using a loosely-coupled solver in which the flow equations and structure equations are advanced sequentially in time. In solving the flow equations, the immersed boundary method is used to impose the no-slip condition on the surface of the swimmer. In solving the structure equations, a finite difference method is used for the spatial discretization and a three-time-level scheme is used for the temporal advancement. This FSI solver has been extensively validated in the earlier works of our group [16–20]. After solving for the velocity field u (x,t), the scalar vorticity field is computed by ωz=(?×u)·k , where k is the unit normal vector which points towards the “outside” of the two-dimensional plane. Here the second-order central-difference scheme is used to discretize the curl operator ( ?×).

    We conducted nine cases of simulation by tuning some control parameters. These parameters include: coefficients of the exponential functions that govern the distributions of β and γ,dimensionless heaving amplitude A /L , pitching amplitude Aθ,and duty cycle D C. Among all simulation cases, three distinct wake types, namely, “2S”, “2P”, “2P+2S”, are identified. The values of the control parameters and wake type corresponding to each simulation case are listed in Table 1. It is seen that each wake type can be observed in three different simulation cases.The representative vortex structures for the “2S”, “2P” and“2P+2S” wakes are shown in Fig. 2.

    The wake detection zone is of rectangular shape and has the dimensions of 5 L×2L. It is symmetric with respect to the horizontal centerline and is placed downstream of the swimmer with a distance of L from the trailing edge. Unlike the previous studies in Refs. [14, 15], where the time histories of flow variables were used for wake classification, here we only use the data of“one time" measurements. The reasons for this choice in the present study are twofold. First, in the wake produced by a selfpropelled swimmer, “decaying with time” is the salient feature in the time histories of flow variables, while temporal periodicity does not show up in the signals. Second, the wake classification based on the data of “one time” measurements implies that there is no time delay. This is certainly more favorable in the design of flow sensory systems in animal (or artificial) swimmers. For each simulation case listed in Table 1, the data sampling frequency is 20 times per flapping cycle, and the duration of sampling is 5 flapping cycles. As a result, the data of 100 snapshots are collected. For each snapshot, the vorticity ( ωz) and two velocity components (uxand uy) at the grid points of the computational mesh are saved. Before later usage in training and testing the neural networks, all flow variables are preprocessed by scaling with the root-mean-square (RMS) values of the collected data.

    A spatially distributed array of senors with the length of L is then used to measure the flow variables (see Fig. 1). The sensor array is in align with the x-axis and the number of sensors is set to 6. For each snapshot, the measurements are conducted 861 times, by placing the first (leftest) sensor on the vertices of a 41×21 lattice which occupies part of the detection zone. The vertices of the lattice are uniformly distributed and the interspace between two neighboring vertices is 0 .1L. Thus, for each simulation case, 861 × 20 × 5 = 86100 measurements are conducted by the sensor array. Since the positions of sensors are not necessarily in coincidence with the positions of grid points, bilinear interpolations are performed to obtain the values of flow variables at the positions of sensors.

    We train fully connected neural networks (FCNs) to classify the wake structures. The hyperparameters of the neural networks are listed in Table 2. Specifically, the FCNs consist of 3 hidden layers and 1 output layer and each hidden layer has 32 units (neurons). We test the effects of the numbers of hidden layers and units on the performance of the networks. It is found that the networks exhibit little variation in performance if the number of hidden layers or the number of units is further increased.The output of each unit is the input of ReLU activation function.The output layer is a softmax function which can map the output of the last hidden layer to a probability distribution vector.This vector has three components which indicate the probability of the three predicted wake types. Thus, the values of thethree components lie in the range of 0 to 1. The loss function of the FCNs is the cross-entropy function which can be written as:

    Table 1 Control parameters of the fish-like swimmers used to produce the wakes.

    where pi(x) is the true probability distribution vector by using one-hot encoding that identifies the wake type and qi(x) is the predicted distribution vector obtained from the FCNs. The training of the FCNs is performed by minimizing the loss function to approach the real distribution of the target. In this work, the neural networks are implemented using the open source machine learning package TensorFlow.

    The FCNs are trained and tested by using the approach of three-fold cross-validation. More specifically, for each wake type, the measurement data obtained from two simulation cases are used for training and the data from the remaining case is reserved for testing. Altogether, three distinct data sets (Data set 1,Data set 2 and Data set 3) are created by different combinations of the simulation cases for training and testing (see Table 3).

    Fig. 2. Representative vortex structures for the three types of wake:a “2S” type; b “2P” type; c “2P+2S” type. The vortex structures are visualized using the contours of dimensionless vorticity (vorticity scaled by f).

    Fig. 3. Values of loss function and network accuracy versus the number of training epochs in training the FCN using the combination of all three flow variables: a loss function; b network accuracy.

    Four types of FCNs are trained using the measurements of different flow variables as the input data. The first three FCNs are trained by using ux, uyand ωz, respectively, while the fourth FCN is trained by using the combination of all three flow variables. In Fig. 3, we show the progress to convergence in training the FCN using all three flow variables as the input. It is seen that 100 epochs are sufficient for the network to reach its best performance, since both the loss function and accuracy exhibit little variation with the increase of epoch number. Similar behaviors are also observed in the training of other three FCNs (the progress to convergence for them are not shown here for brevity).

    We then evaluate the performances of the four FCNs in classifying the wake structures. The performances of the FCNs are visualized using the normalized confusion matrices. The definition of confusion matrix is as follows. Each row of the confusion matrix represents the instances in an actual class, while each column represents the instances in a predicted class. Thus, the diagonal and off-diagonal entries of the matrix show the instances of correct classification and misclassification, respectively. A normalized confusion matrix is constructed by normalizing the entries of each row of the confusion matrix using the total instances in the actual class (i.e., the summation of all entries in that row). The normalized confusion matrices for the three FCNs trained by using ux, uyand ωzare shown in Fig. 4. The ac-curacies shown in the normalized confusion matrices are the averaged results obtained by using three different data sets. It is seen that the FCN trained using uxperforms very well. For all three wake types, the percentage of correct identification are higher than 9 0%. For the “2S” wake, the rate of correct identification can even reach 9 8.3%. The FCN trained by using uyperforms comparably with the FCN trained using ux. For these two FCNs, the highest rate of misclassification is around 6 % to 7 %.The performance of the FCN trained using ωzis much worse. In identifying the “2S”, ”2P” and “2P+2S” wakes, the percentages of misclassification are 6.1%, 58.6%, and 24.4%, respectively.

    Table 2 Hyperparameters of the FCNs used in the wake classification.

    Table 3 Different combinations of simulation cases in the threefold cross-validation for training and testing the neural networks.

    Fig. 4. Normalized confusion matrices for the three FCNs trained using a ux; b uy; c ωz.

    Fig. 5. Contours of the two velocity components in the “2S”, “2P” and “2P+2S” wakes. a and b (dimensionless) streamwise and crosswise velocity components for the “2S” wake; c and d (dimensionless) streamwise and crosswise velocity components for the “2P” wake; e and f (dimensionless) streamwise and crosswise velocity components for the “2P+2S” wake.

    The poor performance of the FCN trained using ωz(especially in identifying the “2P” wakes) can be explained as follows.By examining the vorticity contours shown in Fig. 2, it is seen that for the “2P” wake, in a large portion of the wake detection zone, the vorticity magnitude is very small (close to zero). This makes it very difficult for the FCN to extract the intrinsic feature of the wake from local measurements of ωz. In contrast with the contours of ωz, the contours of uxand uylook very different (see Fig. 5). It is seen that the contours with high magnitudes of uxand uyare less concentratedly distributed in the wake detection zone. This makes it easy for the FCNs to uncover the salient feature of each wake type.ux, uyand ωz) as the input data. It is seen that remarkably high Figure 6 shows the normalized confusion matrix for the FCN trained by using the combination of all three flow variables (i.e.,accuracies in identifying all three wake types have been achieved. The highest rate of misclassification is reduced to 3%(which is observed in identifying the “2P+2S” wakes). The superior performance of the FCN trained using all three variables can be explained as follows. The spatial distributions of the three flow variables measured by the sensor array can complement each other and result in reduction of confusion among the three wake types.

    Fig. 6. Normalized confusion matrix for the FCN trained using the combination of ux, uy and z.

    In this work, we trained four types of FCNs to distinguish among three types of wakes produced by a self-propelled fishlike swimmer. We used the “one time” measurements of some flow variables (including the two velocity components and the vorticity) as the input data in training the neural networks. It turned out that the FCN trained using streamwise or crosswise component of velocity performs much better than the FCN trained using vorticity. The FCN trained by using the combination of all three flow variables outperforms the other three FCNs and achieves a remarkably high accuracy in wake classification.

    There are a number of interesting avenues for future research. First, the sensitivity of the performance of the neural networks to the number and arrangement of sensors needs to be studied. Second, the neural networks can be used to discern more useful information of the swimmer from the wake, such as the body size, species, swimming speed and swimming style.This is achievable only if a solid link between such information and the localized flow measurements can be established. Third,one limit of the current study is the assumption that the placement of sensors in the detection zone does not disturb the flow.The two-way interactions between the sensors and the flow need to be considered.

    Acknowledgement

    This work was supported by the National Natural Science Foundation of China (Grants 11772338 and 11372331), Chinese Academy of Sciences (Grants XDB22040104 and XDA22040203).

    在线观看免费午夜福利视频| 岛国在线观看网站| 99国产极品粉嫩在线观看| 亚洲久久久国产精品| 国产成人精品无人区| 精品国产一区二区久久| 久久国产亚洲av麻豆专区| 亚洲精品在线美女| 国产一卡二卡三卡精品| 色精品久久人妻99蜜桃| 欧美中文综合在线视频| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| 亚洲av熟女| 亚洲免费av在线视频| 一级毛片女人18水好多| 午夜精品国产一区二区电影| av有码第一页| 人人妻人人澡人人看| 人人妻人人添人人爽欧美一区卜| 飞空精品影院首页| 欧美激情极品国产一区二区三区| 无人区码免费观看不卡| 欧美日韩国产mv在线观看视频| 久久久久久人人人人人| 色94色欧美一区二区| 久久性视频一级片| 精品第一国产精品| 又紧又爽又黄一区二区| 国产蜜桃级精品一区二区三区 | 精品福利观看| 新久久久久国产一级毛片| 国产精品免费一区二区三区在线 | 91在线观看av| 涩涩av久久男人的天堂| 露出奶头的视频| 欧美日韩福利视频一区二区| 9191精品国产免费久久| 亚洲国产精品sss在线观看 | 18禁裸乳无遮挡免费网站照片 | 亚洲国产精品sss在线观看 | 人人妻人人澡人人看| 免费一级毛片在线播放高清视频 | 亚洲五月色婷婷综合| 亚洲精品久久成人aⅴ小说| 欧美在线一区亚洲| www.精华液| 亚洲午夜精品一区,二区,三区| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 女人被狂操c到高潮| 国产在视频线精品| 国产高清激情床上av| 久久久久久久精品吃奶| 国内久久婷婷六月综合欲色啪| 9热在线视频观看99| 精品久久久久久,| 黄频高清免费视频| 免费av中文字幕在线| 丝袜美足系列| 国产乱人伦免费视频| 999久久久精品免费观看国产| 伊人久久大香线蕉亚洲五| 国产av一区二区精品久久| 久久精品人人爽人人爽视色| 亚洲精品自拍成人| e午夜精品久久久久久久| 国产高清videossex| 在线国产一区二区在线| 女同久久另类99精品国产91| 999久久久国产精品视频| 亚洲片人在线观看| 视频区欧美日本亚洲| 麻豆av在线久日| 18禁黄网站禁片午夜丰满| 久久午夜综合久久蜜桃| 免费日韩欧美在线观看| 精品福利观看| 很黄的视频免费| 人人澡人人妻人| 午夜日韩欧美国产| 精品国产亚洲在线| 精品人妻1区二区| a级片在线免费高清观看视频| а√天堂www在线а√下载 | 50天的宝宝边吃奶边哭怎么回事| 真人做人爱边吃奶动态| 久久久久久久国产电影| 狠狠婷婷综合久久久久久88av| 18禁裸乳无遮挡动漫免费视频| 好看av亚洲va欧美ⅴa在| 国产在线观看jvid| av超薄肉色丝袜交足视频| 久久ye,这里只有精品| 动漫黄色视频在线观看| 国产av精品麻豆| 精品久久蜜臀av无| 亚洲中文字幕日韩| 国产男女超爽视频在线观看| 成人三级做爰电影| 国产精品一区二区精品视频观看| 久久久水蜜桃国产精品网| 精品国产一区二区久久| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美 日韩 精品 国产| 99re6热这里在线精品视频| 免费观看a级毛片全部| 丁香六月欧美| 婷婷精品国产亚洲av在线 | 久久国产精品影院| 国产欧美日韩精品亚洲av| 精品熟女少妇八av免费久了| 久久精品亚洲熟妇少妇任你| 亚洲成人国产一区在线观看| 久99久视频精品免费| 9191精品国产免费久久| 国产欧美日韩一区二区三| 丝袜人妻中文字幕| 美女扒开内裤让男人捅视频| 精品久久久精品久久久| 精品国产超薄肉色丝袜足j| 一本一本久久a久久精品综合妖精| 久久久国产成人免费| 99国产精品99久久久久| 人妻 亚洲 视频| 啪啪无遮挡十八禁网站| 久久久久国产一级毛片高清牌| 50天的宝宝边吃奶边哭怎么回事| 亚洲va日本ⅴa欧美va伊人久久| 亚洲专区中文字幕在线| a级片在线免费高清观看视频| 亚洲精华国产精华精| 中出人妻视频一区二区| 看黄色毛片网站| 12—13女人毛片做爰片一| 他把我摸到了高潮在线观看| 男女之事视频高清在线观看| 久久这里只有精品19| 精品国内亚洲2022精品成人 | 亚洲av成人av| 久久午夜亚洲精品久久| 啦啦啦视频在线资源免费观看| 国产野战对白在线观看| 丝袜人妻中文字幕| 中文字幕av电影在线播放| 亚洲伊人色综图| 大型av网站在线播放| 欧美大码av| 老司机靠b影院| 精品国产一区二区三区久久久樱花| 成人三级做爰电影| 午夜视频精品福利| 亚洲国产精品一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 欧美精品高潮呻吟av久久| 一边摸一边抽搐一进一小说 | 精品国产乱码久久久久久男人| 国产精品 欧美亚洲| 国产精品久久久久久精品古装| 精品久久久久久久毛片微露脸| 熟女少妇亚洲综合色aaa.| 超碰97精品在线观看| 亚洲性夜色夜夜综合| 久久久久久亚洲精品国产蜜桃av| 99国产精品99久久久久| 飞空精品影院首页| 午夜福利一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 青草久久国产| √禁漫天堂资源中文www| 免费日韩欧美在线观看| 国产又爽黄色视频| 久久中文字幕人妻熟女| 免费在线观看黄色视频的| 欧美久久黑人一区二区| 一级a爱视频在线免费观看| 日韩制服丝袜自拍偷拍| 国产一区二区三区综合在线观看| 后天国语完整版免费观看| 国产精品一区二区免费欧美| 天堂俺去俺来也www色官网| 另类亚洲欧美激情| 久久久久视频综合| 纯流量卡能插随身wifi吗| 成在线人永久免费视频| 桃红色精品国产亚洲av| 国产在线观看jvid| 一个人免费在线观看的高清视频| 亚洲色图 男人天堂 中文字幕| 在线观看午夜福利视频| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 女人被躁到高潮嗷嗷叫费观| 亚洲av美国av| 亚洲五月婷婷丁香| 久久国产精品人妻蜜桃| 国产视频一区二区在线看| 日韩欧美一区二区三区在线观看 | 人人妻人人澡人人看| 国产1区2区3区精品| 纯流量卡能插随身wifi吗| 国产午夜精品久久久久久| 美女 人体艺术 gogo| 99久久99久久久精品蜜桃| av中文乱码字幕在线| 极品教师在线免费播放| 伦理电影免费视频| 免费少妇av软件| 欧美色视频一区免费| 热99国产精品久久久久久7| www.999成人在线观看| 成年女人毛片免费观看观看9 | 一级,二级,三级黄色视频| 精品少妇久久久久久888优播| 亚洲第一av免费看| 午夜老司机福利片| 国产成人欧美| 麻豆av在线久日| 久久精品国产综合久久久| 欧美老熟妇乱子伦牲交| 国产精品免费大片| 亚洲国产毛片av蜜桃av| 国产麻豆69| 国产精品1区2区在线观看. | 久久精品国产综合久久久| 国产一卡二卡三卡精品| 美女扒开内裤让男人捅视频| 精品亚洲成国产av| 欧美人与性动交α欧美精品济南到| av不卡在线播放| 亚洲精品美女久久久久99蜜臀| 大码成人一级视频| 丝瓜视频免费看黄片| 色精品久久人妻99蜜桃| 在线十欧美十亚洲十日本专区| 99国产精品免费福利视频| 大码成人一级视频| 91成人精品电影| 香蕉国产在线看| 成年人午夜在线观看视频| 女性生殖器流出的白浆| 午夜免费鲁丝| 亚洲精品一二三| 亚洲第一欧美日韩一区二区三区| 午夜免费成人在线视频| 欧美精品一区二区免费开放| 亚洲精品成人av观看孕妇| av福利片在线| 欧美精品亚洲一区二区| 中文字幕人妻丝袜制服| 国产成人精品久久二区二区免费| 欧美老熟妇乱子伦牲交| 成人18禁在线播放| av欧美777| 99热网站在线观看| 悠悠久久av| 久久亚洲真实| a级毛片黄视频| 亚洲av电影在线进入| 最近最新中文字幕大全电影3 | x7x7x7水蜜桃| 天天躁夜夜躁狠狠躁躁| 麻豆国产av国片精品| 9热在线视频观看99| 大陆偷拍与自拍| 脱女人内裤的视频| 91精品三级在线观看| 国产精品 欧美亚洲| 色在线成人网| av片东京热男人的天堂| 自线自在国产av| 超色免费av| 一边摸一边抽搐一进一出视频| 天天添夜夜摸| 在线观看免费高清a一片| 99热只有精品国产| 91大片在线观看| av在线播放免费不卡| 亚洲成a人片在线一区二区| 亚洲精品国产区一区二| 国产精华一区二区三区| 午夜激情av网站| 午夜亚洲福利在线播放| 成人永久免费在线观看视频| 精品久久久久久久毛片微露脸| 我的亚洲天堂| 日韩欧美三级三区| 成人18禁高潮啪啪吃奶动态图| 成年动漫av网址| 成人国语在线视频| 免费看a级黄色片| 亚洲色图 男人天堂 中文字幕| 国产精品98久久久久久宅男小说| av网站免费在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 老汉色av国产亚洲站长工具| 50天的宝宝边吃奶边哭怎么回事| 人人妻人人爽人人添夜夜欢视频| 国产在线精品亚洲第一网站| 欧美老熟妇乱子伦牲交| 国产日韩一区二区三区精品不卡| 欧美午夜高清在线| 高清在线国产一区| 国产欧美日韩一区二区三| 校园春色视频在线观看| 精品视频人人做人人爽| 十分钟在线观看高清视频www| 伊人久久大香线蕉亚洲五| 波多野结衣一区麻豆| 国产视频一区二区在线看| 亚洲av欧美aⅴ国产| 99国产精品一区二区三区| 天天添夜夜摸| 国产亚洲精品一区二区www | 一边摸一边做爽爽视频免费| 亚洲五月天丁香| 成年人午夜在线观看视频| 免费观看人在逋| 91在线观看av| 香蕉丝袜av| 亚洲av片天天在线观看| 亚洲av电影在线进入| 日日爽夜夜爽网站| 99国产精品免费福利视频| 最新美女视频免费是黄的| 日韩欧美一区二区三区在线观看 | 波多野结衣av一区二区av| 国产精品秋霞免费鲁丝片| 日韩免费av在线播放| 成年人黄色毛片网站| 中文亚洲av片在线观看爽 | 中文字幕制服av| 亚洲熟女精品中文字幕| 精品国产超薄肉色丝袜足j| 18禁裸乳无遮挡免费网站照片 | 交换朋友夫妻互换小说| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 欧美人与性动交α欧美软件| 男女免费视频国产| 精品少妇久久久久久888优播| 极品少妇高潮喷水抽搐| 99国产精品一区二区蜜桃av | 午夜福利一区二区在线看| 国产欧美日韩综合在线一区二区| 亚洲色图综合在线观看| 精品人妻1区二区| 最新的欧美精品一区二区| 咕卡用的链子| 久久天堂一区二区三区四区| x7x7x7水蜜桃| 国产亚洲av高清不卡| 亚洲九九香蕉| 国产高清videossex| 精品人妻1区二区| 桃红色精品国产亚洲av| 久久精品亚洲熟妇少妇任你| 日本黄色日本黄色录像| 亚洲国产中文字幕在线视频| 久久人人97超碰香蕉20202| 久久亚洲真实| 精品少妇一区二区三区视频日本电影| 嫁个100分男人电影在线观看| 国产亚洲av高清不卡| 99riav亚洲国产免费| 欧美成人午夜精品| 如日韩欧美国产精品一区二区三区| 精品国产超薄肉色丝袜足j| 国产精品成人在线| 一进一出抽搐gif免费好疼 | 看黄色毛片网站| 国产免费现黄频在线看| 一本综合久久免费| 亚洲av日韩在线播放| 757午夜福利合集在线观看| 男女下面插进去视频免费观看| 亚洲一区中文字幕在线| 中文字幕人妻丝袜制服| 看片在线看免费视频| 国产91精品成人一区二区三区| 欧美大码av| 日韩成人在线观看一区二区三区| 日本一区二区免费在线视频| 丝瓜视频免费看黄片| 日韩一卡2卡3卡4卡2021年| 国产在线一区二区三区精| 精品无人区乱码1区二区| 亚洲精品在线观看二区| 亚洲第一欧美日韩一区二区三区| 日韩欧美一区视频在线观看| 欧美久久黑人一区二区| av有码第一页| 女人爽到高潮嗷嗷叫在线视频| 91av网站免费观看| 777久久人妻少妇嫩草av网站| 1024视频免费在线观看| 超碰97精品在线观看| 欧美人与性动交α欧美精品济南到| 午夜视频精品福利| 亚洲成a人片在线一区二区| tube8黄色片| 久久久国产一区二区| 久久精品国产亚洲av香蕉五月 | 日韩视频一区二区在线观看| 国产精品乱码一区二三区的特点 | 欧美激情高清一区二区三区| 岛国在线观看网站| 18禁裸乳无遮挡免费网站照片 | 在线观看午夜福利视频| 黄色 视频免费看| 女警被强在线播放| 村上凉子中文字幕在线| 人妻久久中文字幕网| 欧美黑人精品巨大| 成人18禁高潮啪啪吃奶动态图| 亚洲久久久国产精品| 窝窝影院91人妻| 国产一区二区三区在线臀色熟女 | 在线视频色国产色| 色尼玛亚洲综合影院| 国产一卡二卡三卡精品| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 国产精品99久久99久久久不卡| a在线观看视频网站| 不卡一级毛片| 日本vs欧美在线观看视频| 午夜福利视频在线观看免费| 午夜成年电影在线免费观看| 99国产精品免费福利视频| 99国产精品99久久久久| 亚洲av第一区精品v没综合| 欧美精品亚洲一区二区| 成人黄色视频免费在线看| 亚洲专区中文字幕在线| 精品久久蜜臀av无| 国产91精品成人一区二区三区| 久久青草综合色| 日韩视频一区二区在线观看| 手机成人av网站| 亚洲五月婷婷丁香| 国产深夜福利视频在线观看| 精品久久蜜臀av无| 欧美激情 高清一区二区三区| 欧美最黄视频在线播放免费 | 久久天躁狠狠躁夜夜2o2o| 黄网站色视频无遮挡免费观看| 淫妇啪啪啪对白视频| 日日爽夜夜爽网站| 亚洲av美国av| 国产精品一区二区在线不卡| 美女午夜性视频免费| 亚洲人成伊人成综合网2020| 女人被躁到高潮嗷嗷叫费观| 超色免费av| 麻豆成人av在线观看| 成人黄色视频免费在线看| 人人妻,人人澡人人爽秒播| 老司机福利观看| 他把我摸到了高潮在线观看| 国产1区2区3区精品| 黄频高清免费视频| 国产乱人伦免费视频| 最新的欧美精品一区二区| 老司机靠b影院| 免费在线观看亚洲国产| 久久精品aⅴ一区二区三区四区| www.熟女人妻精品国产| 欧美激情 高清一区二区三区| 一级毛片女人18水好多| 香蕉丝袜av| av不卡在线播放| 又大又爽又粗| 嫩草影视91久久| 国产亚洲欧美精品永久| 日韩欧美一区二区三区在线观看 | av免费在线观看网站| 两性夫妻黄色片| 午夜免费成人在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产欧美日韩一区二区三| 黄色a级毛片大全视频| 午夜福利免费观看在线| 欧美日韩av久久| 亚洲人成电影观看| 嫩草影视91久久| av网站在线播放免费| 人妻一区二区av| 在线观看舔阴道视频| 最近最新免费中文字幕在线| 亚洲精品在线美女| 成人国语在线视频| 国产不卡av网站在线观看| 91精品国产国语对白视频| 亚洲一码二码三码区别大吗| 国产亚洲一区二区精品| 国产97色在线日韩免费| 亚洲视频免费观看视频| 王馨瑶露胸无遮挡在线观看| 母亲3免费完整高清在线观看| 老司机靠b影院| 成年女人毛片免费观看观看9 | x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看 | 中文字幕人妻丝袜一区二区| 搡老熟女国产l中国老女人| 亚洲精品自拍成人| 亚洲伊人色综图| 国产精品久久电影中文字幕 | 最新美女视频免费是黄的| 90打野战视频偷拍视频| 久热爱精品视频在线9| 99久久国产精品久久久| 久久精品aⅴ一区二区三区四区| 日日爽夜夜爽网站| 欧美精品人与动牲交sv欧美| 纯流量卡能插随身wifi吗| 久久久国产一区二区| 国产一区二区三区综合在线观看| www日本在线高清视频| 久久午夜亚洲精品久久| 在线观看一区二区三区激情| 一级毛片女人18水好多| 99国产精品一区二区蜜桃av | 精品午夜福利视频在线观看一区| 十八禁网站免费在线| 热99国产精品久久久久久7| 制服诱惑二区| 啦啦啦 在线观看视频| 亚洲自偷自拍图片 自拍| 成年人黄色毛片网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久国产精品人妻蜜桃| 免费女性裸体啪啪无遮挡网站| 久久久久久免费高清国产稀缺| 如日韩欧美国产精品一区二区三区| 飞空精品影院首页| 十分钟在线观看高清视频www| 亚洲,欧美精品.| 一区二区三区国产精品乱码| 欧美激情高清一区二区三区| 国产精品国产高清国产av | 精品亚洲成国产av| 黄片播放在线免费| 亚洲av日韩在线播放| 国产精品综合久久久久久久免费 | 久久精品国产亚洲av高清一级| 国产成人欧美| 久久久国产成人精品二区 | 99热国产这里只有精品6| 国产精品香港三级国产av潘金莲| 一级毛片高清免费大全| 这个男人来自地球电影免费观看| 人人妻,人人澡人人爽秒播| 91大片在线观看| 国产亚洲精品久久久久5区| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 日韩熟女老妇一区二区性免费视频| 国产精品免费视频内射| 免费在线观看亚洲国产| 欧美黄色片欧美黄色片| www日本在线高清视频| 操出白浆在线播放| 国产精品国产av在线观看| 国产亚洲精品一区二区www | 国产在线观看jvid| 啦啦啦在线免费观看视频4| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 久久人妻av系列| 狠狠狠狠99中文字幕| 18禁观看日本| 国产精品一区二区在线不卡| 国产伦人伦偷精品视频| 悠悠久久av| 男女高潮啪啪啪动态图| 香蕉丝袜av| 欧美 日韩 精品 国产| 黑人猛操日本美女一级片| 日韩视频一区二区在线观看| 国产99白浆流出| 老司机福利观看| 成人亚洲精品一区在线观看| 麻豆乱淫一区二区| 欧美一级毛片孕妇| 午夜福利免费观看在线| 动漫黄色视频在线观看| 国产亚洲欧美精品永久| 午夜福利在线免费观看网站| 国产99白浆流出| 亚洲国产中文字幕在线视频| 国产av精品麻豆| 12—13女人毛片做爰片一| 高清在线国产一区| 天堂中文最新版在线下载| 老熟妇仑乱视频hdxx| 国产高清国产精品国产三级| 亚洲七黄色美女视频| 韩国精品一区二区三区| 精品一区二区三区av网在线观看| 国产精品久久久久久人妻精品电影| 亚洲少妇的诱惑av| 91老司机精品| 亚洲av欧美aⅴ国产| 国产日韩欧美亚洲二区| 国产精品电影一区二区三区 | av线在线观看网站| 我的亚洲天堂| 久久天躁狠狠躁夜夜2o2o| 亚洲国产精品sss在线观看 | 九色亚洲精品在线播放| 最新在线观看一区二区三区| 国产精品九九99| 一边摸一边做爽爽视频免费|