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

    A Novel Method Based on UNET for Bearing Fault Diagnosis

    2021-12-10 11:54:20DileepKumarSootherImtiazHussainKalwarTanweerHussainBhawaniShankarChowdhrySanaullahMehranUjjanandTayabDinMemon
    Computers Materials&Continua 2021年10期

    Dileep Kumar Soother,Imtiaz Hussain Kalwar,Tanweer Hussain,Bhawani Shankar Chowdhry,Sanaullah Mehran Ujjan and Tayab Din Memon

    1National Centre of Robotics and Automation,HHCMS Lab,Mehran University of Engineering&Technology,Jamshoro,76020,Sindh,Pakistan

    2Department of Electrical Engineering,DHA SUFFA University,Karachi,Sindh,Pakistan

    3Department of Electronic Engineering,Mehran University of Engineering and Technology,Jamshoro,76020,Sindh,Pakistan

    Abstract:Reliability of rotating machines is highly dependent on the smooth rolling of bearings.Thus,it is very essential for reliable operation of rotating machines to monitor the working condition of bearings using suitable fault diagnosis and condition monitoring approach.In the recent past,Deep Learning(DL)has become applicable in condition monitoring of rotating machines owing to its performance.This paper proposes a novel bearing fault diagnosis method based on the processing and analysis of the vibration images.The proposed method is the UNET model that is a recent development in DL models.The model is applied to the 2D vibration images obtained by transforming normalized amplitudes of the time-series vibration data samples into the corresponding vibration images.The UNET model performs pixel-level feature learning using the vibration images owing to its unique architecture.The results demonstrate that the model can perform dense predictions without any loss of label information,generally caused by the sliding window labelling method.The comparative analysis with other DL models confirmed the superiority of the UNET model which has achieved maximum accuracy of 98.91%and F1-Score of 99%.

    Keywords:Condition monitoring;deep learning;fault diagnosis;rotating machines;vibration

    1 Introduction

    Bearings are critical components in rotating machines used in various applications such as auto-manufacturing,chemical industries,railways,water pumping stations,etc.Their health has significant impacts on the operation of industrial systems and economy of industries,with 45%–55% of faults in industrial motors are caused by bearing faults[1,2].Therefore,effective fault diagnosis of bearings can significantly increase the reliability and life span of machines in various industrial applications,which in turn reduces maintenance and operation costs[3,4].Considering the importance of bearings in rotating machines,various methods have been applied by researchers to effectively monitor the bearing condition and to avoid downtime caused by the bearing failure.

    Vibration is considered as an early and highly accurate indicator of various faults in mechanical equipment,and its measurement and analysis are widely employed with Artificial Intelligence(AI)based fault diagnosis and prognosis methods[5–7].Typically,a fault diagnosis procedure comprises of three steps:data acquisition,feature extraction,and fault diagnosis.Firstly,the data acquisition step for machinery fault detection involves signals such as vibration,current,voltage,temperature,and acoustic emissions.Second step includes the extraction of time-domain features(root mean square,skewness,kurtosis,and gap factor),frequency domain(Fourier transform),and time-frequency domain features(wavelet transform,short-time Fourier transform,empirical mode decomposition,and Hilbert-Huang transform)[8–13].Lastly,fault diagnosis step which has been widely explored by researchers using different model-based techniques and AI-based techniques.Machine Learning(ML)based methods are the foundation for effective fault diagnosis in the arena of AI.Various researches in literature have employed these techniques for fault diagnosis,for instance,Support Vector Machines(SVM),K-Nearest Neighbors(KNN),and Random Forest(RF)[14–16].

    ML methods have yet proven to be an effective diagnosis method;however,these methods have some inherent limitations and their performance extremely relies on designed features.For combined bearing fault detection,it often requires advanced signal processing techniques[17].Moreover,these methods are not effective in a practical industrial environment where background noise and interference of signal components are inevitable[18].Thus,the progress in AI led to the emergence of DL which has been widely exploited by researchers in various fields of this era.DL as a subfield of AI and has provided an efficient way of automatic representative feature learning,even from raw and noisy input data without any human intervention[19,20].Compared to ML methods,DL methods automatically learn rich features even from raw data using high-performance computing hardware[7,21–24].Various DL architectures have been exploited for efficient bearing fault detection including Convolutional Neural Networks(CNN),Long Short-Term Memory(LSTM),Auto-Encoders(AE),Deep Belief Networks(DBN),and Deep Boltzmann Machines(DBM)[3,25–28].Among these DL models,CNNs has been the most representative model and has demonstrated robust and remarkable performance as a supervised learning approach[29,30].Furthermore,Full Convolutional Network(FCN)can perform pixellevel semantic segmentation of image data owing to their rich hierarchical architecture.It has also made substantial progress in achieving higher level of feature hierarchy by processing of lower level features.These models learn features through complex mapping which extracts spectral information from individual pixels[3,6].Considering the merits of CNN model,it has been widely used to diagnose rolling bearing fault using 1D(one dimensional)and 2D(two dimensional)data with 1D-CNN and 2D-CNN configurations,respectively.Guo et al.[31]have used an adaptive CNN(ADCNN)to classify the 4 conditions of the rolling bearing with the Case Western Reserve University(CWRU)bearing data.The initial layers of the model comprise of LeNet-5 and then stacked hierarchical layers which determine two output components:fault type and fault size.The adaptive configuration of the learning rate allowed to maintain a trade-off between accuracy and training speed.This model yielded competitive performance to the earlier models.

    In[32],a deep CNN model with multisensory vibration data for classifying the nine bearing conditions have been experimented.Their model was able to learn features from the raw data and avoided overfitting through dropout regularization.Li et al.[33]have employed deep CNN with data augmentation.The model was trained with 400 data samples of the CWRU bearing dataset with the ten bearing conditions.Raw vibration data was the input to the CNN layers then residual layers were stacked for better feature learning.The network was trained with batch-normalization and dropout layers for fast training and avoiding overfitting through data augmentation.Although these DL architectures have demonstrated effective performance in bearing fault diagnosis,but these models face multiclass window problem and loss of information owing to sliding window labelling technique.Thus,these models pose limitations in performing dense predictions.

    Unlike the above-cited works,a novel bearing fault diagnosis method based on UNET and vibration images is presented in this paper.The proposed approach overcomes the multiclass window problem that occur in DL models and can predict label for each input data sample in the bearing data.The UNET model achieves pixel-level dense predictions through the downsampling and up-sampling layers in the architecture.This model is successfully applied in various domains and has revealed excellent performance in terms of dense predictions[34–36].The main contribution of this research article is twofold.The primary contribution is UNET based novel rolling bearing fault diagnosis method,which can perform dense predictions by overcoming the multiclass window problem,caused by sliding window-labelling technique.To the best of authors’knowledge,UNET network is applied to diagnose bearing faults for the first time.The second is the dense labelled 2D vibration images without any data-augmentation given as input to the UNET model for sample-based bearing condition classification.

    The remaining paper is organized as follows:Section 2 discusses dense labelling and fault classification using UNET,Section 3 explains the UNET architecture and training process,Section 4 discusses the dataset,experiment configuration,and experiment evaluation,Section 5 reports and discusses the results,and Section 6 concludes the investigation.

    2 Dense Predictions Using UNET

    Conventionally,time-series or sequential signals are divided into fixed-length windows using sliding window techniques.Then,each window is assigned with the same label and this label assigning technique leads to a multiclass window problem.There are two sliding window techniques:the first technique assigns a label by selecting the most frequent sample class and the second technique assigns labels by selecting last sample class of time step.Using both the techniques,classifier generates incorrect output information and causes a decrease in accuracy owing to the multiclass window problem.This problem is solved by dense labelling technique which increases classifiers’accuracy by providing correct label information.This method assigns label to each sample of dataset rather than labelling based on sliding windows.Hence,it upholds all the label information in dataset.Fig.1 describes both the labelling techniques.Current sample and next sample are described as S1 and S2.In the case of sliding window labelling method,the most frequent appearing class is assigned with label S1.However,owing to sliding windows it contains information of both classes.Thus,it causes learning incorrect information which in turn brings down the recognition accuracy.

    The traditional CNN restricts implementation of dense prediction,whereas UNET model reduces resolution of top-level output by using the max-pooling operation.For image recognition task,UNET reduces sensitivity of image shift.In the case of time-series data recognition,CNN causes a mismatch between input data length and output label length.To overcome this problem,UNET architecture has been successfully implemented by adding up-convolution layers or upsampling layers in CNN architecture.It allows the use of dense labelling and preserves same length for input series and output label.Thus,UNET can be applied for dense predictions owing to addressing the multiclass window problem which is often posed by traditional CNN.

    Figure 1:Sliding-window labelling and dense labelling

    3 Bearing Fault Classification Using UNET

    3.1 UNET Architecture

    UNET is the recent development on CNN that utilizes deep layer architecture for automatic feature learning from input data.The multilayer architecture allows to learn more abstract features and can perform classification at pixel-level.To achieve pixel-level classification,an end-to-end UNET architecture was proposed by Ronneberger et al.[37]for semantic segmentation of the biomedical images.The main objectives were to improve the precise segmentation and localization of neuronal composition in the microscopic images.The goal behind the development of UNET model was to realize pixel-level classification through an in-depth feature processing.

    UNET architecture consists of two paths:a contractive path and an expansive path and both are symmetric to each other and yield an architecture like U-shape as shown in Fig.2(Thus called as UNET).The network on the left side is contracting path that is like a traditional CNN involving convolution and pooling layers with activations.Thus,learning the image contents.On the right side of Fig.2 is an expansive path which includes stacked up-sampling layers and the corresponding convolutional layers on the left side.Both the network paths are merged to compensate for the loss of information caused by pooling operation.As a result,the architecture preserves the same resolution of images as in input network layer.

    Fig.2 shows that in the contracting network,in which each block consists of the two repeated and unpadded 3 × 3 convolution layer followed by a ReLU(Rectified Linear Unit)activation and 2 × 2 max-pooling layer with a stride of 2 for down-sampling.Whereas,the expansive network consists of the same number of up-sampling layers as the down-sampling layers.In the end,a full-connection layer and softmax classifier are employed to map input feature vectors with the corresponding output class.In a contractive network,convolutional and pooling layers transform high-resolution images to low-resolution images also called down-sampling.While in an expansive network,the transposed convolutional layers transform low-resolution images to high-resolution images also known as up-sampling.Transposed convolution functions exactly the opposite way of normal convolution thus known as de-convolution.It up-samples the images by learning parameters through backpropagation.

    In this paper,a 9-level UNET architecture is employed for efficient pixel-level feature learning in the bearing fault classification.The architecture comprises of three parts:encoding network,decoding network,and a bridge that connects both the networks as given in Tab.1.The encoding network transforms vibration images into compact representations and the decoding network recovers the transformed representations as pixel-wise categorizations.The complete network is constructed using 3×3 convolution layers,pooling layers,and a ReLU activation function.

    Figure 2:UNET architecture

    3.2 Network Training

    The employed network architecture is shown in Fig.2.The network receives input vibration images with size(N,N,C),where N denotes the number of sampling points and C denotes the number of channels.The architecture consists a total of 18 layers including contracting and expanding network layers.The input nodesxi,jreceive three-dimensional inputhl×wl×fl,wherehlandwlare the height and width of the input vector at a location(i,j)andflis number of the input feature maps in layerl.The output of layerlis represented byyi,jand can be found by the formula given in Eq.(1):

    wherekiandkjrepresents convolution kernels andf(·)denotes the type of layer(non-linear operation with activation,matrix multiplication of convolution layer,max-pooling operation of maximum cell layer,etc.).The output of the layerlbecomes input of the layerl+1.Its size is given byhl+1×wl+1×fl+1;wherehl+1=(hl?ki)/si+1 andwl+1=(wl?kj)/sj+1.The stride of movement is expressed bysiandsj.Application of the filling operation on input feature map allows to maintain resolution of output feature map same as the input feature map.

    Table 1:Structure of the UNET network

    To achieve effective and accurate dense predictions using UNET,the network training is performed for estimating appropriate parameters(W,b)from a given dataset and corresponding labels.Through minimizing loss of all the samples in training dataset,accurate dense predictions can be made possible.The loss function is expressed in Eq.(2):

    4 Experimental Study

    This section reports experiments performed on a rolling bearing dataset using UNET model.Furthermore,the investigation and evaluation process of performance of UNET is reported in this section.In addition,this section presents a comparison of UNET model performance with other DL architectures,which are employed in the recent research for bearing fault diagnosis.Firstly,the dataset employed in this research is explained,then the experimental configuration and evaluation metrics are discussed.

    4.1 Dataset and Vibration Image Construction

    In this investigation,the CWRU bearing dataset is used to evaluate performance of the proposed model[38].The reason behind choosing this dataset is that it is the most widely employed by the researchers as a benchmark dataset for bearing fault diagnosis[4].Furthermore,this dataset is provided with open access for research community which allows them to evaluate the performance of their proposed algorithms.The experimental setup includes a dynamometer(right),a 2 HP Reliance electric motor(left),and a torque encoder(center)as depicted in Fig.3.The setup also includes a control system for proper operation of the system;however,it is not shown in the figure.The test bearings are attached to support the motor shaft.Various single point artificial faults of different sizes are introduced to the rolling bearings using electric discharge machining.The vibration is collected using accelerometers attached to the fan-end and derive-end,and base of the motor.These accelerometers are attached using magnetic bases with position at the 12’o clock.The data is acquired using a 16 channel data recorder.The setup is operated under the four different load conditions including 0–3 HP load applied by the dynamometer within range of the speed between 1720 and 1797.

    Figure 3:Experimental setup

    In this paper,a balanced dataset is used,which includes different conditions of the CWRU bearing dataset with 1 HP load.These conditions of bearing include healthy condition and nine bearing faults of different sizes in mils.The dataset includes a total of 10 classes of the vibration data that is collected from the drive-end of the setup with a sampling rate of 48 kHz.The size of each condition is given in Tab.2 with the appropriate class labels.

    The vibration signatures of bearings obtained as 1D data can be transformed into 2D form as images.The process includes normalization of each sample in vibration signal into the range[?1,1].Then,the normalized amplitude of each sample of the signal transformed into the corresponding pixel in the vibration image[39].This transformation between normalized amplitude of vibration signal and pixels can be expressed as given in Eq.(3):

    where P[i,j]is the intensity of the corresponding pixels in M * N vibration image,i = 1:N,and j= 1:M.A[·]represents normalized amplitude of each vibration sample.The number of samples in a vibration signal are transformed into an equal number of pixels in the vibration image.

    Table 2:Bearing fault class labels and data samples

    The dataset including these bearing conditions is transformed into the vibration images with the minimum resolution.The transformation resulted in a total of 4600 vibration images of size 32 × 32 × 1.Fig.4 shows the transformed vibration images of the ten bearing conditions with the appropriate class labels.

    4.2 Experiment Configuration

    To train the UNET network,the vibration images of different bearing conditions are given as input.The size of the each vibration image is 32 × 32 × 1 including height,width,and the number of channels,respectively.The model is trained with the Adam optimizer,batch-size of 128,and the dropout rate of 0.2.The learning rate is set adaptively with the minimum value of 1 × 108.The adaptive learning rate callback allowed to reduce learning by 0.1,if the model validation loss does not improve for five iterations.The early stopping callback allowed to avoid overfitting of the model through stopping the training,if the validation loss does not improve for ten continuous iterations.Subsequently,it saves the learned model parameters.In the end,the softmax loss function is used to classify the ten different conditions of the bearing.The dataset is randomly divided into training,test,and validation set with a percentage of 70%,20%,and 10%,respectively.

    Figure 4:Vibration images of the bearing conditions

    To compare the performance of the proposed model,three DL models are employed which were investigated in the recent researches on this dataset.These models include FCN[40],LeNet-5[41],and ResNet-50[42].These models are described briefly in the following paragraphs.

    ?FCN:This model consists of 6 layers including:two convolutional layers,two max-pooling layers,and two Full Connection(FC)layers.The model is employed for bearing condition diagnosis with a dropout of 0.2,batch size of 128,and zero padding.At the final stage,a softmax classifier is stacked.

    ?LeNet-5:It is a 2D network comprises of nine layers including three convolutional layers,three max-pooling layers,and three FC layers.A softmax layer is added at the end of the model to classify the bearing conditions.The dropout is set at 0.2 and batch-size of 128.

    ?ResNet-50:ResNet-50 model can optimize model parameters through a “shortcut connection or residual unit” in each convolutional block.It also assists in avoiding overfitting without loss of important information.This model consists of two blocks:stacked structure block and two FC blocks.The stacked structure blocks achieve low training parameters through reducing the size of feature map which in turn minimizes hardware requirement.Each stacked block consists of convolution layer,max-pooling layer,and batch-normalization.Meanwhile,FC blocks allow classifying the bearing conditions.Here,this model is utilized with a dropout rate of 0.2 and a batch size of 128.

    The hardware used in this investigation is provided with a GPU model GeForce RTX-2060.The system is programmed with Python 3.7.1 and DL frameworks including Keras and Tensorflow.

    4.3 Experiment Evaluation

    To compare the performance of the UNET model and comparative models,following evaluation indexes are used as given in Eqs.(4)–(7):

    where TP denotes true positives,TN denotes true negatives,FP denotes the false positives,FN denotes false negatives,and m denotes number of examples in data.

    5 Results and Discussions

    In this section,the performance of the four models is compared and summarized based on accuracy and F1-Score indexes.The highest results obtained with the UNET and comparative methods are shown in Tab.3.The results clearly indicate effectiveness of the UNET method,which outperforms the comparative DL methods owing to its capability of dense predictions.The UNET has achieved the best accuracy of 98.91% and the highest F1-Score of 99% among all the employed DL models on the dataset in this research.FCN has achieved 97.61% of accuracy and 97.6% of F1-Score.While,LetNet-5 has demonstrated the second-best performance with accuracy and F1-Score of 96.74% and 96.8%,respectively.However,ResNet-50 has demonstrated the lowest performance in diagnosing the rolling bearing conditions.The comparative models have yielded lower accuracies than the UNET owing to loss of useful information caused by sliding window labelling technique.Contrarily,UNET model preserves useful information from the vibration images owing its hierarchical network architecture which in turn allowed it to perform better than the comparative models.The achieved results conclude the UNET model as an effective model among all the used DL models in this investigation.

    Table 3:The evaluation indexes of the UNET and the comparative methods on CWRU bearing dataset

    Fig.5 shows the F1-Score of the UNET model and comparative DL models in terms of the each bearing condition.Comparatively,UNET model has the best F1-Score for each bearing condition.

    Figure 5:F1-Score of the UNET and comparative methods for each bearing condition

    Fig.6 shows the confusion matrices of the UNET and comparative DL methods in terms of the each bearing condition.As shown in Fig.6a,UNET model accuracy for the individual classes varies within the range of 95.65% to 100%.Similarly,individual class accuracy of the LeNet-5 is shown in Fig.6b,where the accuracy has the minimum value of 84.75% and the maximum of 100%.In Fig.6c,individual class accuracy of the ResNet-50 model is shown which varies within the range of 82.61% and 100%.Lastly,confusion matrix of the FCN model shows the individual class accuracy changes between 82.61% and 100%.These confusion matrices confirm the UNET model as the best model for the individual bearing condition classification among the employed DL models.

    Figure 6:Confusion matrix of the(a)UNET,(b)LeNet-5,(c)ResNet-50,and(d)FCN

    Figure 7:Accuracy and loss of the(a)UNET,(b)LeNet-5,(c)ResNet-50,and(d)FCN

    Fig.7 shows training accuracy and loss function of the employed DL models.These models were trained with early stopping callback which stops the training when a DL model achieves maximum accuracy without overfitting.Thus,number of the epochs varies for each DL model.Fig.7a shows that the accuracy of the UNET model tends to be stable after the 30th epoch and achieves the accuracy of 98.91% and its loss function value improved after the 25th epoch.The LeNet-5 model achieves the accuracy of 96.74%,but the model response keeps fluctuating as shown in Fig.7b.Similarly,the ResNet-50 model achieves the accuracy of 95.43% with the fluctuating response for both the accuracy and the loss function as depicted in Fig.7c.Lastly,Fig.7d shows accuracy of the FCN that is around 97.61% and its accuracy and the loss function response improve after 20 epochs.It can be noted from the figures that the UNET model depicts the most stable classification performance with the highest classification accuracy.

    Although the UNET model was fed with the vibration images of resolution of 32 × 32 ×1,yet it yielded the excellent performance in terms of classifying the bearing conditions owing to the pixel-level learning capabilities.Moreover,with the adaptive learning rate the model learned features hierarchically thorough adjusting the learning rate,if the performance of the model did not improve for certain iterations.Another technique which is known as early stopping,helped to avoid the overfitting,if the model did not improve the performance.The UNET model has depicted robustness in the training process as shown in Fig.7a.

    The UNET model was remarkably able to predict sample-based rolling bearing conditions owing to its inherent property of dense predictions.Its comparison with the other popular DL models revealed that the UNET model has superior and robust performance in terms of rolling bearing fault diagnosis.It was observed from the results that the proposed algorithm possesses the excellent potential for fault diagnosis of rotating machines in various industrial applications.

    6 Conclusion

    In this paper,a novel DL method namely UNET is proposed for rolling bearing fault classification based on the vibration images.Compared to the existing DL methods,this model overcomes the multiclass window problem that is inherent in the sliding window labelling method.The UNET model used in this investigation has effectively performed dense predictions of the bearing conditions with an accuracy of 98.91% and F1-Score of 99%.The obtained results confirm that the UNET model outperforms the comparative DL models such as CNN,LeNet-5,and ResNet-50.The model yielded robust and better results than the comparative methods on short-term feature recognition.The robust classification results of the UNET model on this bearing dataset indicates its excellent potential for applications in other domains.

    Funding Statement:Authors would like to acknowledge the support of the ‘Haptics,Human Robotics,and Condition Monitoring Lab’established in Mehran University of Engineering and Technology,Jamshoro under the umbrella of the National Center of Robotics and Automation.This work was supported by the Higher Education Commission Pakistan(Grant No.2(1076)/HEC/M&E/2018/704).

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding this research work.

    久久精品亚洲熟妇少妇任你| 高潮久久久久久久久久久不卡| 啦啦啦在线免费观看视频4| 一本综合久久免费| 妹子高潮喷水视频| 成年动漫av网址| 少妇裸体淫交视频免费看高清 | 日韩制服丝袜自拍偷拍| 成人免费观看视频高清| 一区二区av电影网| 日本av手机在线免费观看| 侵犯人妻中文字幕一二三四区| 亚洲第一av免费看| 天天操日日干夜夜撸| 十分钟在线观看高清视频www| 丁香六月欧美| 黑人巨大精品欧美一区二区mp4| 国产淫语在线视频| 男女边摸边吃奶| 自拍欧美九色日韩亚洲蝌蚪91| 黑人操中国人逼视频| 老熟妇仑乱视频hdxx| 亚洲av美国av| 男人操女人黄网站| 脱女人内裤的视频| 悠悠久久av| 亚洲av美国av| 日韩有码中文字幕| 亚洲精品久久午夜乱码| 中文字幕制服av| 99国产精品免费福利视频| 中国美女看黄片| 男女高潮啪啪啪动态图| 午夜两性在线视频| 国产高清激情床上av| 国产免费福利视频在线观看| 多毛熟女@视频| 日韩免费高清中文字幕av| 欧美日韩黄片免| 欧美激情 高清一区二区三区| 午夜福利在线观看吧| 2018国产大陆天天弄谢| 变态另类成人亚洲欧美熟女 | 91大片在线观看| 国产精品免费视频内射| 久久中文看片网| 麻豆av在线久日| 国产精品欧美亚洲77777| 欧美激情久久久久久爽电影 | 人人澡人人妻人| 久热爱精品视频在线9| 精品国产亚洲在线| 欧美激情久久久久久爽电影 | 欧美日韩亚洲综合一区二区三区_| 国产成人av教育| 久久 成人 亚洲| 午夜视频精品福利| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 99九九在线精品视频| 嫩草影视91久久| 欧美国产精品va在线观看不卡| 乱人伦中国视频| 90打野战视频偷拍视频| 日韩中文字幕视频在线看片| 9热在线视频观看99| 美女午夜性视频免费| 日本a在线网址| 亚洲成a人片在线一区二区| 在线观看舔阴道视频| 99国产极品粉嫩在线观看| 午夜福利免费观看在线| av又黄又爽大尺度在线免费看| 中文字幕精品免费在线观看视频| 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 久久 成人 亚洲| 亚洲av成人不卡在线观看播放网| 国产伦人伦偷精品视频| 成人18禁在线播放| 国产日韩欧美视频二区| 国产欧美亚洲国产| 亚洲免费av在线视频| 一本综合久久免费| 久久久欧美国产精品| 人成视频在线观看免费观看| 亚洲avbb在线观看| 在线 av 中文字幕| a级片在线免费高清观看视频| 欧美精品人与动牲交sv欧美| 国产高清videossex| 在线观看免费视频日本深夜| 国产一区二区三区综合在线观看| 国产97色在线日韩免费| 侵犯人妻中文字幕一二三四区| 欧美激情 高清一区二区三区| 后天国语完整版免费观看| 精品国产乱子伦一区二区三区| 丝袜在线中文字幕| 国产精品av久久久久免费| 成人国产av品久久久| 欧美激情极品国产一区二区三区| 国产成人av教育| 精品久久久精品久久久| 考比视频在线观看| 欧美日韩一级在线毛片| 欧美性长视频在线观看| 久久毛片免费看一区二区三区| 亚洲中文日韩欧美视频| 日日摸夜夜添夜夜添小说| 18禁观看日本| 免费女性裸体啪啪无遮挡网站| 欧美黄色淫秽网站| 国产人伦9x9x在线观看| 超碰97精品在线观看| 99在线人妻在线中文字幕 | 两个人看的免费小视频| 国产成人欧美| 99国产综合亚洲精品| av又黄又爽大尺度在线免费看| 777久久人妻少妇嫩草av网站| 在线观看免费视频网站a站| 欧美精品一区二区大全| 精品少妇久久久久久888优播| 亚洲av成人不卡在线观看播放网| 日韩人妻精品一区2区三区| 黄色丝袜av网址大全| 变态另类成人亚洲欧美熟女 | 天天操日日干夜夜撸| 老司机午夜十八禁免费视频| 香蕉久久夜色| 亚洲国产看品久久| 久9热在线精品视频| 成人精品一区二区免费| 波多野结衣av一区二区av| 欧美精品一区二区大全| 天堂8中文在线网| 午夜视频精品福利| 母亲3免费完整高清在线观看| 搡老乐熟女国产| 每晚都被弄得嗷嗷叫到高潮| 黄色视频,在线免费观看| 精品欧美一区二区三区在线| 欧美日韩亚洲高清精品| 亚洲一码二码三码区别大吗| 亚洲精品av麻豆狂野| 女人高潮潮喷娇喘18禁视频| 交换朋友夫妻互换小说| 日日摸夜夜添夜夜添小说| 波多野结衣av一区二区av| 久久免费观看电影| 一级毛片精品| 久久久国产欧美日韩av| 麻豆av在线久日| 国产99久久九九免费精品| 国产精品影院久久| 满18在线观看网站| 亚洲色图综合在线观看| 国产又色又爽无遮挡免费看| 黄片大片在线免费观看| 欧美性长视频在线观看| 久久久水蜜桃国产精品网| 99久久精品国产亚洲精品| 99国产精品一区二区三区| 国产欧美日韩一区二区三| 后天国语完整版免费观看| 精品久久久精品久久久| 亚洲男人天堂网一区| 午夜老司机福利片| 五月天丁香电影| 久久天堂一区二区三区四区| 一区二区av电影网| 日本黄色视频三级网站网址 | 女人精品久久久久毛片| 亚洲精品一卡2卡三卡4卡5卡| 99re6热这里在线精品视频| 国产麻豆69| 成人av一区二区三区在线看| 一个人免费在线观看的高清视频| 免费久久久久久久精品成人欧美视频| 久久精品国产a三级三级三级| 国产成人精品久久二区二区91| 91老司机精品| 国产主播在线观看一区二区| 国产成人免费观看mmmm| 妹子高潮喷水视频| 国产精品久久久久久人妻精品电影 | 免费在线观看视频国产中文字幕亚洲| 十八禁网站免费在线| 国产成人欧美在线观看 | 久久中文看片网| 亚洲人成伊人成综合网2020| 精品一区二区三区视频在线观看免费 | 一级毛片电影观看| 老汉色∧v一级毛片| 女性生殖器流出的白浆| 日本vs欧美在线观看视频| 夜夜骑夜夜射夜夜干| 欧美变态另类bdsm刘玥| 亚洲成国产人片在线观看| 侵犯人妻中文字幕一二三四区| 国产精品.久久久| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产区一区二| 亚洲精华国产精华精| 一级毛片精品| 亚洲第一青青草原| 一本色道久久久久久精品综合| 精品国产亚洲在线| 夜夜爽天天搞| 亚洲专区字幕在线| 欧美激情极品国产一区二区三区| 亚洲avbb在线观看| 久久久国产欧美日韩av| 黄色片一级片一级黄色片| 老司机在亚洲福利影院| 亚洲九九香蕉| 窝窝影院91人妻| 日韩 欧美 亚洲 中文字幕| 制服诱惑二区| 国产精品1区2区在线观看. | 一边摸一边抽搐一进一小说 | 涩涩av久久男人的天堂| 18在线观看网站| netflix在线观看网站| av视频免费观看在线观看| 国产精品秋霞免费鲁丝片| 亚洲色图av天堂| 午夜福利免费观看在线| 一级毛片女人18水好多| 18禁国产床啪视频网站| www.精华液| 婷婷成人精品国产| 在线看a的网站| 中文字幕av电影在线播放| 国产在线观看jvid| 日韩有码中文字幕| 一边摸一边做爽爽视频免费| 大型av网站在线播放| 国产免费视频播放在线视频| 久久av网站| 国产单亲对白刺激| 99久久人妻综合| 这个男人来自地球电影免费观看| a级毛片黄视频| 色94色欧美一区二区| 成人影院久久| 午夜福利欧美成人| 精品久久久久久电影网| 在线亚洲精品国产二区图片欧美| a在线观看视频网站| 99国产精品一区二区三区| 成年人免费黄色播放视频| 欧美国产精品一级二级三级| av天堂在线播放| 国产色视频综合| 中文字幕精品免费在线观看视频| 99re6热这里在线精品视频| 2018国产大陆天天弄谢| 丁香欧美五月| 最新在线观看一区二区三区| 大型黄色视频在线免费观看| 国产熟女午夜一区二区三区| 黑人操中国人逼视频| 亚洲成av片中文字幕在线观看| 欧美另类亚洲清纯唯美| 波多野结衣av一区二区av| 别揉我奶头~嗯~啊~动态视频| 亚洲色图av天堂| 99国产综合亚洲精品| 精品国内亚洲2022精品成人 | 日韩欧美三级三区| 大片电影免费在线观看免费| 丁香欧美五月| 亚洲精品乱久久久久久| 久久婷婷成人综合色麻豆| 国产免费福利视频在线观看| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全电影3 | 99精品欧美一区二区三区四区| 熟女少妇亚洲综合色aaa.| 日本五十路高清| 后天国语完整版免费观看| 一二三四在线观看免费中文在| 肉色欧美久久久久久久蜜桃| 精品国产国语对白av| 99精品欧美一区二区三区四区| 亚洲欧洲日产国产| 精品一品国产午夜福利视频| 精品久久蜜臀av无| 国产在线观看jvid| 中文字幕色久视频| 精品亚洲乱码少妇综合久久| 免费观看a级毛片全部| 水蜜桃什么品种好| 亚洲一码二码三码区别大吗| 757午夜福利合集在线观看| 亚洲熟妇熟女久久| 精品一区二区三区av网在线观看 | 18禁裸乳无遮挡动漫免费视频| 男人舔女人的私密视频| 操出白浆在线播放| 成人特级黄色片久久久久久久 | 女人爽到高潮嗷嗷叫在线视频| 国产午夜精品久久久久久| 国产精品 欧美亚洲| 午夜福利视频在线观看免费| e午夜精品久久久久久久| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 纯流量卡能插随身wifi吗| 成人手机av| av免费在线观看网站| 久久精品国产亚洲av高清一级| 午夜91福利影院| 国产欧美亚洲国产| 丁香六月欧美| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲熟女毛片儿| 男女下面插进去视频免费观看| 日韩免费高清中文字幕av| 日韩熟女老妇一区二区性免费视频| 日本av手机在线免费观看| 亚洲成人国产一区在线观看| 成人影院久久| a级片在线免费高清观看视频| 国产精品 欧美亚洲| 国产不卡一卡二| 首页视频小说图片口味搜索| 亚洲欧美色中文字幕在线| 日本wwww免费看| 日本欧美视频一区| 成人国产av品久久久| 国产成人一区二区三区免费视频网站| 三上悠亚av全集在线观看| 黄片大片在线免费观看| 在线观看www视频免费| 日韩大码丰满熟妇| 激情视频va一区二区三区| 这个男人来自地球电影免费观看| 午夜两性在线视频| 两性夫妻黄色片| 建设人人有责人人尽责人人享有的| 国产精品 国内视频| 亚洲午夜精品一区,二区,三区| 少妇裸体淫交视频免费看高清 | 国产成人av教育| 日韩免费av在线播放| 亚洲色图 男人天堂 中文字幕| 日本黄色视频三级网站网址 | 9色porny在线观看| 日本精品一区二区三区蜜桃| 久久青草综合色| 一二三四社区在线视频社区8| 91大片在线观看| 热99久久久久精品小说推荐| 两个人看的免费小视频| www.999成人在线观看| 青青草视频在线视频观看| 午夜福利,免费看| 久久热在线av| 久久久久精品人妻al黑| 99国产精品99久久久久| 免费高清在线观看日韩| 高清欧美精品videossex| 国产av精品麻豆| netflix在线观看网站| 亚洲人成77777在线视频| 看免费av毛片| 国产免费现黄频在线看| 一本综合久久免费| 国产成人精品久久二区二区91| 女性被躁到高潮视频| 亚洲av日韩在线播放| 国产精品熟女久久久久浪| 久久精品国产亚洲av高清一级| www日本在线高清视频| 天堂中文最新版在线下载| 99国产精品一区二区三区| 日本a在线网址| 欧美激情久久久久久爽电影 | 午夜视频精品福利| 久久这里只有精品19| 成人国产一区最新在线观看| av片东京热男人的天堂| 国产xxxxx性猛交| 大码成人一级视频| 日韩 欧美 亚洲 中文字幕| 免费久久久久久久精品成人欧美视频| 亚洲视频免费观看视频| 精品国产乱子伦一区二区三区| 欧美日本中文国产一区发布| 少妇猛男粗大的猛烈进出视频| 久久精品人人爽人人爽视色| 搡老熟女国产l中国老女人| 精品国产乱码久久久久久男人| 高清欧美精品videossex| 亚洲七黄色美女视频| 一边摸一边抽搐一进一出视频| 纵有疾风起免费观看全集完整版| 一二三四在线观看免费中文在| 日韩大码丰满熟妇| 女人被躁到高潮嗷嗷叫费观| 性少妇av在线| 亚洲国产看品久久| 国产片内射在线| a级片在线免费高清观看视频| 天堂8中文在线网| 在线天堂中文资源库| 国产精品电影一区二区三区 | 亚洲伊人久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 精品欧美一区二区三区在线| 久久久欧美国产精品| 久久久久久久久免费视频了| 视频区欧美日本亚洲| 亚洲av成人一区二区三| 91av网站免费观看| 国产精品99久久99久久久不卡| a在线观看视频网站| 亚洲精品国产区一区二| 亚洲国产av新网站| 巨乳人妻的诱惑在线观看| 国产亚洲欧美精品永久| 黄片小视频在线播放| 亚洲 欧美一区二区三区| 91成年电影在线观看| 夜夜骑夜夜射夜夜干| 纵有疾风起免费观看全集完整版| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 亚洲 国产 在线| 色婷婷av一区二区三区视频| 久久狼人影院| 黄片播放在线免费| 成人免费观看视频高清| 欧美性长视频在线观看| 日韩大码丰满熟妇| 天天躁夜夜躁狠狠躁躁| 亚洲第一欧美日韩一区二区三区 | 亚洲精品乱久久久久久| 在线永久观看黄色视频| 这个男人来自地球电影免费观看| 精品国产亚洲在线| 亚洲国产av新网站| 久久精品熟女亚洲av麻豆精品| 国产高清videossex| 一二三四社区在线视频社区8| 国产成人精品久久二区二区免费| 一夜夜www| 亚洲色图av天堂| 国产精品成人在线| 日韩制服丝袜自拍偷拍| 久久久国产一区二区| 亚洲视频免费观看视频| 久久国产亚洲av麻豆专区| 香蕉丝袜av| 国产精品1区2区在线观看. | 亚洲国产欧美在线一区| 国产精品免费一区二区三区在线 | 高清欧美精品videossex| av网站免费在线观看视频| 精品一区二区三区视频在线观看免费 | 丝袜喷水一区| av一本久久久久| 亚洲熟妇熟女久久| 我要看黄色一级片免费的| 久久天躁狠狠躁夜夜2o2o| 黄色 视频免费看| e午夜精品久久久久久久| 亚洲精品久久成人aⅴ小说| 天堂8中文在线网| 国产精品一区二区在线观看99| 欧美性长视频在线观看| 久久精品亚洲av国产电影网| 97在线人人人人妻| 国产精品 欧美亚洲| 十八禁网站网址无遮挡| 女人爽到高潮嗷嗷叫在线视频| 老鸭窝网址在线观看| 精品一区二区三区视频在线观看免费 | 国产在线观看jvid| 最近最新中文字幕大全电影3 | 精品福利永久在线观看| 老司机靠b影院| 日本vs欧美在线观看视频| 亚洲精品久久成人aⅴ小说| 亚洲欧美日韩高清在线视频 | 国产精品免费视频内射| 国产精品美女特级片免费视频播放器 | 99精品久久久久人妻精品| 一个人免费看片子| 欧美久久黑人一区二区| 亚洲人成77777在线视频| 欧美日韩精品网址| 亚洲色图 男人天堂 中文字幕| videos熟女内射| 亚洲精品国产色婷婷电影| 国产精品.久久久| 成年人午夜在线观看视频| 欧美成狂野欧美在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲欧洲精品一区二区精品久久久| 亚洲熟妇熟女久久| 岛国在线观看网站| 如日韩欧美国产精品一区二区三区| 新久久久久国产一级毛片| www.熟女人妻精品国产| 国产亚洲欧美在线一区二区| 女人被躁到高潮嗷嗷叫费观| 精品欧美一区二区三区在线| 精品少妇黑人巨大在线播放| 亚洲va日本ⅴa欧美va伊人久久| 久久精品亚洲av国产电影网| 久久久国产欧美日韩av| 国产一区二区 视频在线| 免费高清在线观看日韩| 精品少妇一区二区三区视频日本电影| 精品国产一区二区久久| 午夜福利乱码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 中文字幕另类日韩欧美亚洲嫩草| 丰满饥渴人妻一区二区三| 日韩一卡2卡3卡4卡2021年| 日韩欧美一区二区三区在线观看 | av不卡在线播放| 两性夫妻黄色片| 少妇精品久久久久久久| 午夜免费鲁丝| 在线看a的网站| 悠悠久久av| 午夜激情久久久久久久| 啦啦啦 在线观看视频| 久久久久网色| 国产欧美日韩综合在线一区二区| av欧美777| 日日夜夜操网爽| 欧美在线一区亚洲| 久久久久久人人人人人| 亚洲欧洲日产国产| 露出奶头的视频| 亚洲五月色婷婷综合| 久久久精品国产亚洲av高清涩受| 精品人妻在线不人妻| 亚洲中文日韩欧美视频| 这个男人来自地球电影免费观看| 国产日韩欧美亚洲二区| 91大片在线观看| 香蕉丝袜av| 另类精品久久| 午夜激情久久久久久久| 五月开心婷婷网| 一级a爱视频在线免费观看| 精品一区二区三区四区五区乱码| 国产一卡二卡三卡精品| 无限看片的www在线观看| 免费人妻精品一区二区三区视频| 精品亚洲成国产av| 亚洲情色 制服丝袜| 国产有黄有色有爽视频| 精品国产超薄肉色丝袜足j| av又黄又爽大尺度在线免费看| 亚洲精品美女久久av网站| 一边摸一边抽搐一进一出视频| 一进一出抽搐动态| 国产不卡av网站在线观看| 久久精品国产综合久久久| 国产精品偷伦视频观看了| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 最近最新免费中文字幕在线| 老熟女久久久| 精品第一国产精品| 久久99一区二区三区| 久久久久久人人人人人| www.熟女人妻精品国产| 香蕉丝袜av| 成人av一区二区三区在线看| 欧美国产精品一级二级三级| 精品国产一区二区三区久久久樱花| 激情视频va一区二区三区| 精品午夜福利视频在线观看一区 | 青青草视频在线视频观看| 亚洲av日韩在线播放| 久久人人97超碰香蕉20202| 国产高清videossex| 国产亚洲精品久久久久5区| 欧美 日韩 精品 国产| 国产精品久久电影中文字幕 | 一级,二级,三级黄色视频| 涩涩av久久男人的天堂| 精品少妇久久久久久888优播| 久久亚洲真实| 97在线人人人人妻| 怎么达到女性高潮| 中亚洲国语对白在线视频| 免费黄频网站在线观看国产| 久久青草综合色| 午夜福利在线观看吧| 国产野战对白在线观看| 淫妇啪啪啪对白视频| 精品福利永久在线观看| 免费看a级黄色片| 国产99久久九九免费精品| 中国美女看黄片| 国产在线观看jvid| 成人精品一区二区免费| 青草久久国产| 欧美黑人精品巨大| 一本综合久久免费| 精品国产乱码久久久久久小说| 国产91精品成人一区二区三区 | 欧美日韩视频精品一区| 久久国产精品男人的天堂亚洲|