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

    A Compensation Controller Based on a Nonlinear Wavelet Neural Network for Continuous Material Processing Operations

    2019-11-07 03:13:04ChenShenYoupingChenBingChenandJingmingXie
    Computers Materials&Continua 2019年10期

    Chen Shen,Youping ChenBing Chen and Jingming Xie

    Abstract:Continuous material processing operations like printing and textiles manufacturing are conducted under highly variable conditions due to changes in the environment and/or in the materials being processed.As such,the processing parameters require robust real-time adjustment appropriate to the conditions of a nonlinear system.This paper addresses this issue by presenting a hybrid feedforward-feedback nonlinear model predictive controller for continuous material processing operations.The adaptive feedback control strategy of the controller augments the standard feedforward control to ensure improved robustness and compensation for environmental disturbances and/or parameter uncertainties.Thus,the controller can reduce the need for manual adjustments.The controller applies nonlinear generalized predictive control to generate an adaptive control signal for attaining robust performance.A wavelet-based neural network model is adopted as the prediction model with high prediction precision and time-frequency localization characteristics.Online training is utilized to predict uncertain system dynamics by tuning the wavelet neural network parameters and the controller parameters adaptively.The performance of the controller algorithm is verified by both simulation,and in a real-time practical application involving a single-input single-output double-zone sliver drafting system used in textiles manufacturing.Both the simulation and practical results demonstrate an excellent control performance in terms of the mean thickness and coefficient of variation of output slivers,which verifies the effectiveness of this approach in improving the long-term uniformity of slivers.

    Keywords:Continuous material processing,wavelet neural network(WNN),nonlinear generalized predictive control(NGPC),auto-leveling system.

    1 Introduction

    Textiles manufacturing,printing,and other continuous material processing operations involve essentially nonlinear dynamic systems with coupling between the multiple inputs and multiple outputs of the manifold physical process[Valenzuela Bentley and Lorenz(2004);Tsakalis,Dash,Green et al.(2002);Djiev(2016);Moghassem and Fallahpour(2011)].These processing operations differ significantly from metal cutting processing operations according to two primary features.First,the materials being processed(e.g.,cotton,paper,and polymer films)are handled by multiple processing units that must maintain a high level of mechanical synchronization.For example,slivers are spun into yarn after several times of drawing,and the multicolor printing of paper requires multiple printing operations with precise registration and post-processing.Second,the quality of the end product is not only related to the precise control of processing units(e.g.,sliver drafting motors and roller printing motors),but also depends on the continuous physical processes applied to materials,such as fiber drawing processes and the application of tension to paper materials in printing processes.Moreover,these physical processes are closely related to the material characteristics and numerous external environmental factors.Thus,the quality of the final product cannot be assured by simply improving the control precision of processing units.Meanwhile,manual adjustments are unavoidable for addressing these changes in external and internal conditions.

    The issue of control system robustness has been addressed by the development of many control strategies.For example,generalized predictive control(GPC),which was proposed in the 1980's[Clarke,Mohtadi and Tuffs(1987)],has been one of the most frequently used model-based control methods adopted in industry[Camacho and Bordons(1999);Low and Cao(2008)].In most GPC applications,linear models are employed to predict plant behavior over the prediction horizon,and to evaluate future sequences of control signals.Nevertheless,controllers based on linear models perform poorly when applied to nonlinear systems that operate over a wide range of operating conditions like continuous material processing systems.In this respect,numerous studies have analyzed the feasibility of GPC for the modeling and control of nonlinear system dynamic[Gu and Hu(2002)].For example,nonlinear GPC(NGPC)models have been developed successfully using conventional network structures such as artificial neural networks(ANNs),which have been shown to be capable of approximating a wide range of nonlinear functions to any desired degree of accuracy under specified conditions[Geng and Geary(1997);Hamdia,Lahmer,Nguyen-Thoi et al.(2015)].The neural network and hybrid regression models are applied in material design and producing field[Badawy,Msekh,Hamdia et al.(2017);Faizollahzadeh Ardabili,Najafi,Alizamir et al.(2018);Fardad,Najafi,Ardabili et al.(2018)].In addition,wavelets have been incorporated with ANNs to develop wavelet neural networks(WNNs)that combine the capability of ANNs to learn from processes together with the high resolution of wavelet decomposition[Delyon,Juditsky and Benveniste(1995);Chen and Hsiao(1999)].Here,the second layer of a WNN employs a wavelet transform rather than a standard activation function like the sigmoid function employed in conventional ANNs.The use of the wavelet transform allows for exceptional localization in the time domain via translation of the mother wavelet(a shifting process)and also in the frequency domain via dilation of the mother wavelet(a scaling process).Moreover,this time localization capability in particular readily lends the use of wavelet transforms to non-stationary signal analysis.Thus,the structure of WNNs have been demonstrated to provide greater potential than conventional ANNs for enriching the mapping relationships between inputs and outputs[Delyon,Juditsky and Benveniste(1995)].As such,WNNs are ideally suited for the modeling and control of dynamic systems[Billings and Wei(2005);Lu(2009);Abiyev and Kaynak(2008);Sureshbabu and Farrell(1999)].In addition,the training algorithms adopted for WNNs typically converge in a fewer number of iterations than those adopted for conventional ANNs.However,the standard feedforward network structure of WNNs is not the most suitable for solving temporal problems like predicting the behaviors of complex chaotic systems.To address this issue,Yoo et al.[Yoo,Park and Choi(2005)]developed a self-recurrent wavelet neural network(SRWNN)that combined the attractor dynamics of a recurrent neural network(RNN)with the rapid convergence of a WNN.Based on SRWNN,they developed a GPC for stable path tracking of mobile robots[Yoo,Choi and Park(2006)].Lu[Lu(2009)]proposed a stable predictive controller(SPC)for a class of nonlinear discrete time system using GPC with recurrent WNN(RWNN)model.This type of controller has its simplicity in parallelism to conventional GPC design and efficiency to deal with complex nonlinear dynamics.

    In terms of processes requiring robust control,the sliver drafting process employed in textiles manufacturing represents a very good example of a continuous material processing operation with nonlinear dynamics.Here,the dynamics of the sliver drafting process play an essential role in the automatic control systems employed in textiles manufacturing for the reduction of sliver irregularities.Because a textile sliver is composed of thousands of discrete fibers,a sliver drafting system is complex,and includes inherent nonlinearity,which makes the control of such systems a challenging task.The long-term auto-leveling control of such systems is particularly challenging because instability and external disturbances are frequently observed.While a number of modern effectual control methods have been developed for auto-leveling sliver drafting systems[Kim,Lim and Huh(2012);Moghassem and Fallahpour(2011);Han,Youn-Sung,Soon-Yong et al.(2008)],and the problems associated with the modeling and control of the sliver drafting processes have been discussed in control engineering literature[Djiev(2016)],most of these methods can be applied in practical applications only when at least an approximate mathematical model of the process is available.As such,these methods lack robustness and are difficult to implement in nonlinear systems.Some studies have addressed this problem by implementing linearization methods[Guo,Chen and Hu(2003);Huang and Bai(2001)].However,attempts at linearization have not significantly contributed to the robustness of system control.In this sense,the feature of an SRWNN in the absence of a pre-existing mathematical model is uniquely advantageous.

    In light of the above research focused on improving control system robustness in continuous material processing operations like the sliver drafting process,this paper proposes a hybrid feedforward-feedback nonlinear model predictive controller.Block diagram of the proposed control system is shown in Fig.1.The adaptive feedback control strategy of the controller augments the standard feedforward control to ensure improved robustness and compensation for environmental disturbances and/or parameter uncertainties.Thus,the controller can reduce the need for manual adjustments.The controller adopts a wavelet-based neural network model with high prediction precision and time-frequency localization characteristics.Online training is utilized to predict uncertain system dynamics by tuning the wavelet neural network parameters and the controller parameters adaptively.The controller applies NGPC based on a SRWNN to generate an adaptive control signal for attaining robust performance.

    Figure 1:Block diagram of the proposed control system

    This paper is organized as follows.Section 2 presents the adopted SRWNN model structure and adaptive learning algorithm.A brief description of the adaptive predictive control is introduced in Section 3,and the design procedures of the proposed control system are also described in detail.Section 4 presents the self-leveling process and dynamics of a standard sliver drafting system.The effectiveness of the proposed control system is validated by numerical simulations and experimental results involving a single-input single-output double-zone sliver drafting system.Finally,Section 5 concludes the paper.

    2 SRWNN model structure and training algorithm

    To obtain a prediction model for the proposed prediction control scheme that is appropriate to the characteristics of continuous material processing systems,we adopt a class of nonlinear autoregressive moving average(NARX)time-series model at time stepkas the system model:

    wheref(·):R n→Ris a smooth-valued nonlinear function,kindicates the time step index,y(k)is the system output for a system inputu(k),nyandnuare the orders of the output and input,respectively,N=ny+nu,andξ(k)∈Ris a zero-mean Gaussian white noise sequence.

    2.1 SRWNN model structure

    The control strategy is implemented using an SRWNN model developed as an approximation of the nonlinear system given by(1).The proposed SRWNN adopted a mother wavelet layer composed of self-feedback neurons that can store the past information of the network,allowing the SRWNN to capture the dynamic response of a system.As a result,the SRWNN can be applied effectively to complex chaotic systems,even though the SRWNN employs fewer wavelet nodes than a WNN.Thus,the structure of the SRWNN can be simpler than that of a corresponding WNN,which makes the SRWNN more appropriate than WNNs for real-time control applications.

    Figure 2:Schematic of the proposed SRWNN structure

    A schematic illustrating a WNN and the proposed SRWNN structure is given in Fig.2.The structure comprises an input layer(layer I),a self-recurrent wavelet layer(layer II),a wavelet layer(layer III),and an output layer(layer IV)providing an estimated system output?(k).Here,xi(k),i=1,…,n,is theith input of the SRWNN,θil,l=1,…,L,denotes thelth self-feedback weight of theith input directed toward a node of layer II,z-1represents delay,αiis the connection weight between the input nodes of layer I and the output nodes of layer IV,andωilis the connection weight between the product nodes of layer III and its output nodes.The signal propagation and the function expression in each layer are given as follows.

    Layer I:The nodes in this layer receive the input variables and transmit them to the next layer directly.

    Layer II:Each node of this layer has a wavelon and a self-feedback loop.The first derivative of a Gaussian functionis adopted as a mother wavelet.The output of this layer can be represented as

    Layer III:The output of this layer is the product of all the wavelet neurons,i.e.,

    Layer IV:The single-node output of this layer is a summation which combining the outputs of layer III and the adaptive input values from the input of layer I.Therefore,the output of the SRWNN is composed of all incoming signals and the assigned tuning parameters as follows:

    The initial values ofωilandmilare given randomly in the range[-1,1],whiledilis given randomly in the range(0,1](i.e.,dil>0).In addition,the initial values ofθilare given as 0,indicating an absence of feedback initially.

    The aforementioned SRWNN can be used to be a universal uniform approximator for continuous functions over compact sets when satisfies specific conditions.These conditions and a detailed proof of this claim can be obtained elsewhere[Lin and Chen(2006);Lee and Teng(2000)].The leveling process discussed in our paper is a complex parametric nonlinear dynamical system,which has the random-like behavior usually shown in statistical systems although it is associated with deterministic dynamics.It can be seen as continuous function over compact sets and in the calculation process,suitable parameter constraints can make sure it satisfies the conditions.

    2.2 Online training algorithm for the SRWNN

    A back-propagation algorithm was adopted for SRWNN training in the purposed control system,and all weights includingmilanddilare trained via the gradient descent algorithm.The target is to minimize the following error functionJ0(k):

    Applying the gradient descent method provides the following updating laws formil,dil,θilandωil.

    Here,ηis the learning rate,and the following definitions are applied.

    3 Adaptive non-linear control strategy

    The proposed control strategy was implemented here using a nonlinear SRWNN model as a predictor and the schematic of the control strategy is given in Fig.3.The task of the predictor is to predict plant output based on the regressed inputs at each sampling time.This is conducted for all control operations within an established prediction range.The value of the control horizon should always be less than the value of the prediction horizon.A diophantine equation is used to solve and minimize the complex real-time optimization cost function at each sampling time to determine the optimum control inputs that yield the least error between the predicted output and the trajectory reference signals and which minimize the controller efforts[Astrom and Wittenmark(1994)].

    Figure 3:Structure of the predictive control strategy

    To derive the NGPC law and to find the j step-ahead prediction of y(k),the SRWNN model(4)is rewritten as

    whereΔ=1-z-1is the difference operator,

    The system parametersailandbil,that isωil,inA(z -1)andB(z -1)are estimated online with variable forgetting factor recursive least square method(VFFRLS)adaptively.From Eq.(12),

    The parameters are obtained by the following formula:

    λ(0 <λ<1)is the forgetting factor.

    In order to get the optimal predictions,Eq.(12)could be simplified as

    The proposed control law is derived to minimize the expected valueE[·]of the following predictive performance criterion:

    Here,npis the prediction output horizon,andyr(k+j)is a known bounded reference output for the discrete timek+j.In general,npis chosen to encompass all the responses that are significantly affected by the present control.Here,np Tsis typically the same magnitude as the rise time of the controlled system based on the sampling timeTs[Astrom and Wittenmark(1994)].

    In order to optimize the cost functionJ(k),the predictiony(k+j)forj ≥ 1andj ≤ Npwill be obtained.Consider the following Diophantine equation:

    where the following definitions are applied.

    The resulting overall system has been proven to be stable[Lu(2009)].Assuming that the parametersmil,dil,θilandωilin(4)are updated according to(6)-(8)and(14),the proposed SRWNN algorithm is convergent,provided thatηsatisfies the following condition:

    This convergent condition of the purposed process are verified for a chosen Lyapunov function,and the detailed proof can be obtained elsewhere[Yoo,Park and Choi(2005)].To guarantee thatηresides within this stable region,we apply an adaptive learning rate for the SRWNN as follows:

    From the above,the implementation of the adaptive control procedure can be summarized as follows,and the state flow chart of the calculation process is shown in Fig.4

    (2)Measure plant outputy(k)and the reference tracking outputyr(k+j);

    (4)Solve the Diophantine Eq.(17)and obtain;

    (5)Construct the vector R,Y(k),Δ U(k-j)and matrix F1,F2,G;

    (6)Calculate and implementu(k)via(19);

    (7)Repeat Steps 2-6.

    Figure 4:State flow chart of the calculation process

    4 Numerical simulations and experimental results

    4.1 Sliver drafting process dynamics and auto-leveling control

    A frequently-used auto-leveling system for a double-zone sliver drafter is illustrated in Fig.5.The positions of the fibers with respect to each other and the numbers of fibers within cross-sections are affected by the different speeds of the rollers[Hlava(2003)].At each sampling timek,displacement sensors 1 and 2 obtain the thickness of the input sliverx(k)and the thickness of the output slivery(k),respectively,from which the linear densities of the respective slivers are calculated.In this system,variations in the input sliver linear density can be seen as the dominating disturbance,which greatly affects the uniformity of the output sliver.Usually,this variation is stochastic with occasional step-like disturbances owing to nubs in the input sliver,or when changes in the characteristics of the input sliver(e.g.,batch of material)affect.In addition to this measured disturbance,a number of other unmeasured disturbancesv(k)(e.g.,environment moisture change,stacking of cotton fiber)can affect the quality of the output sliver,and these can also be treated as stochastic disturbances.

    A conventional auto-leveling feedforward control algorithm is a simple linear adjustment based on reference valueyr(k)and the measurementsx(k)from sensor1.In existing feedforward strategies,x(k)serves as an independent input that directly affectsy(k)according to a pre-established parametric model used to describe the drawing process.Accordingly,the irregularity of the output sliver is managed by sending the signal from sensor 1 to the controller to determine the control signal,which then adjusts the velocities of the drawing rollers using the servomotor to obtain an appropriate drawing rate.The signal from sensor 2 strictly serves a monitoring purpose,without participating in the automatic control process.This process allows the controller to compensate immediately for the impact of variations inx(k)ony(k)rather than waiting until the effect appears in the output.Thus,feedforward control ensures the short-term uniformity of the output sliver.In addition,the parameters of the parametric model employed in feedforward control are subject to change if the input sliver material or its characteristics are altered.

    Figure 5:Schematic of a standard feedforward controlled auto-leveling sliver drafting system

    The proposed compensation controller consists of a conventional auto-leveling control strategy in the feedforward component in addition to an online trained NGPC controller in the feedback loop.In this control method,the signal of Sensor 2 is used in the controller to automatically compensate for performance degradation.While the feedforward controller compensates immediately for the measured disturbance,the proposed SRWNN-based model predictive controller provides feedback compensation for unmeasured disturbancesv(k).Here,unmeasured disturbances represent an independent input that is not affected by the controller or the plant,and is always potentially present,but is only observable fromy(k).Unmeasured disturbances represent unknown,unpredictable events that are best addressed as un-modeled system dynamics.Therefore,all information based on the outputy(k)obtained from Sensor 2 is fed back to the NGPC controller.Then,the proposed NGPC controller observes the future behavior of the sliver drafting system and compares the actual performance to a desired reference model performance,and accordingly calculates the control input that will optimize plant performance over a specified future time horizon.The learning algorithm modifies the parameters of the NGPC controller online based on the model-following error(MFE)to obtain a match with the desired reference model response.

    4.2 Numerical simulation of the auto-leveling control system

    Simulation results are presented to verify the feasibility of the proposed SRWNN-based NGPC control scheme under various operating conditions.All algorithms were developed using MATLAB/SIMULINK in a control computer.In the simulations,the disturbance rejection capabilities of the auto-leveling system were evaluated under different long-term low-frequency disturbancesv(k)effecting on the plant.

    Simulation 1:

    At first,the two sensors data of feedforward control,x(k)andy0(k)(voltage),collected in advance are used for offline simulation.A nonlinear model for the sliver drafting process was identified by an identification experiment with a sampling period of 0.001 second,and this model serves as the plant model in our simulation[Chun,Bae,Kim et al.(2006)].Referring to the system model(1),the input variables of the SRWNN model are specified by{y(k-1),y(k-2),y(k-3),u(k-1),u(k-2)}.After removing mean,training the network parameters of the SRWNN using the input-output data,selecting the key parameter of the SRWNN asL=3 was found to be effective for this nonlinear system.The prediction horizon of the proposed control law was selected asNp=9.Then the proposed controller was employed to match the system outputy(k)to a reference output yr(k),where yr(k)=0 and the external disturbancev(k)was specified as follows:

    Case 1:v(k)=0.07sign(sin(4πkts)),

    Case 2:v(k)=0.07sin(4πkts),

    where sign(·)represents the signum function andtsis the sample time.These two cases represent long term perturbation under a step disturbance and a sinusoidal disturbance,respectively.The results for case 1 and case 2 are shown in Figs.6(a)and 6(b),respectively.The individual plots show the input thickness signalsx(k),and the output thickness signalsy(k),control errorse(k),and the control signalsu(k)of the conventional feedforward auto-leveling controller(denoted by the subscript 0)and the proposed controller.

    These figures clearly indicate that good dynamic performances,in terms of command tracking and drift restraint,are realized for the auto-leveling system.From the change of input and output sensor datax(k)andy(k),the evenness of the sliver has improved after auto-leveling,in addition,the error with the average sliver thicknessey(k)is smaller thaney0(k),which means that the control effect is enhanced by the proposed method comparing with the traditional approach.Seen fromu(k),the algorithm is convergence with step disturbance(a;case 1)and sinusoidal disturbance(b;case 2).

    Figure 6:The input thickness signals x(k),and the responses of the conventional feedforward auto-leveling controller(denoted by the subscript 0)and the proposed controller for a reference output yr(k)=0 with a step disturbance(a;case 1)and a sinusoidal disturbance(b;case 2)

    Table 1:Simulation results of feedforward control(FF)alone,FF in conjunction with a standard GPC controller(FF+GPC),and the proposed control system(FF+NGPC)in terms of the controlled output y(k)relative to a reference output yr(k)=0

    Simulation 2:

    In order to verify the validity of the algorithm for long segment evenness,20000 sampling of output sliver thickness withTs=0.01 s are analyzed.The simulation results are listed in Tab.1 in terms of the mean and standard deviations ofy0(k)andy(k).We note from the table that,compared with the results obtained using the standard feedforward control,the mean values ofy(k)obtained using the proposed controller were decreased by 87.18% and 94.95% for Cases 1 and 2,respectively,and the corresponding standard deviations ofy(k)were decreased by 30.48% and 19.78%,respectively.Accordingly,we can conclude that the proposed controller can greatly reduce the nonuniformity of the output sliver relative to that obtained using only feedforward control.Tab.I also includes the mean and standard deviations ofy(k)obtained using feedforward control in conjunction with a conventional GPC controller[Clarke,Mohtadi and Tuffs(1987)].As seen from the results,the performance of the proposed controller is significantly better than that get from the conventional GPC controller.

    Simulation 3:

    To test the robustness of the proposed auto-leveling controller,we conducted a simulation in which the system parametersailandbilwere varied by ±5% after plant operation for 300 s.The results are shown in Fig.7,which illustrate that the proposed SRWNN-based NGPC controller can adapt to arbitrary changes of the plant parameters.The control error reveals no significant increase after the parameter changed.

    The simulation results indicate that the SRWNN-based NGPC controller demonstrates satisfactory tracking performance and system robustness.

    Figure 7:Output tracking and error responses in case of system parameter variations operative after 300 s

    4.3 Experiments on a double-zone sliver drafting system

    The experimental results were obtained using an STM32 microcontroller integrated circuit(STM Electronics)and a field programmable gate array(FPGA)control board.The input,self-recurrent wavelet,wavelet,and output layers of the SRWNN included 2,14,7,and 1 neurons,respectively.The output response in terms ofy(k)with the proposed SRWNN-based NGPC controller is shown in Fig.8.Here,the values were estimated from 20,000 consecutively sampledy(k)data with hardware sampling frequency 1000Hz.Tab.II lists the mean thickness and coefficient of variation(CV)of the output slivers respectively obtained using feedforward control alone and using the proposed controller to regulate the speed of the back roller.The values here are recorded by a USTER evenness tester,which is a professional sliver testing equipment in textile factories and laboratories.Obviously,the control method proposed in this paper can regulate the mean thickness at the desired value,and the CV of the sliver can be substantially reduced.

    Figure 8:The measures value of displacement Sensor 2(y(k))obtained over time using the proposed SRWNN-based NGPC controller

    Table 2:Output sliver thickness values obtained during experimental testing

    5 Conclusion

    This paper proposed an online-trained adaptive SRWNN-based model prediction controller for continuous material processing systems.The SRWNN was employed for establishing a discrete-time model for the nonlinear system dynamic,and the NGPC controller functioned as adaptive feedback compensation for augmenting the existing open-loop feedforward control,and for providing improved setting value tracking and external disturbance resisting capabilities.The proposed NGPC algorithm,including the adaptive learning rate for the training of SRWNN model weights,was applied to a sliver drafting process,and the simulation results indicates the better stable tracking ability and adaptability,comparing with the traditional control strategy.The physical experiments verified that the proposed control system is effective for ensuring the long-term uniformity of slivers both with and without input sliver irregularities and external noise disturbances.

    免费人成视频x8x8入口观看| 中国国产av一级| 青春草视频在线免费观看| 免费看av在线观看网站| 日本三级黄在线观看| 国产成人aa在线观看| 天堂√8在线中文| 亚洲av免费在线观看| 亚洲真实伦在线观看| 插逼视频在线观看| 成人高潮视频无遮挡免费网站| 少妇的逼水好多| 国产午夜精品论理片| 国产激情偷乱视频一区二区| 国产亚洲av片在线观看秒播厂 | 看黄色毛片网站| www.av在线官网国产| 黄色配什么色好看| a级一级毛片免费在线观看| 最后的刺客免费高清国语| 亚洲内射少妇av| 九色成人免费人妻av| 国产精品人妻久久久影院| 久久亚洲精品不卡| 此物有八面人人有两片| 久久久久久九九精品二区国产| 亚洲人成网站在线播放欧美日韩| 天天一区二区日本电影三级| 亚洲,欧美,日韩| 亚洲三级黄色毛片| 日韩视频在线欧美| 十八禁国产超污无遮挡网站| 国产黄色视频一区二区在线观看 | 深夜精品福利| 国产精品三级大全| 国产国拍精品亚洲av在线观看| 亚洲欧洲国产日韩| 亚洲欧美精品专区久久| 国产极品天堂在线| 国产熟女欧美一区二区| 尾随美女入室| 欧美日本视频| 联通29元200g的流量卡| 欧美高清性xxxxhd video| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品456在线播放app| 99热6这里只有精品| 一本一本综合久久| 99在线视频只有这里精品首页| 男人和女人高潮做爰伦理| 久久午夜福利片| 丝袜美腿在线中文| 亚洲欧美日韩卡通动漫| 国产成人影院久久av| 国产精品一区www在线观看| 最近2019中文字幕mv第一页| 成人鲁丝片一二三区免费| a级毛色黄片| 国产精品蜜桃在线观看 | 久久精品夜色国产| 狂野欧美白嫩少妇大欣赏| 97热精品久久久久久| 夜夜爽天天搞| 国产黄a三级三级三级人| 久久这里只有精品中国| 欧美又色又爽又黄视频| 日本免费一区二区三区高清不卡| 你懂的网址亚洲精品在线观看 | 亚洲精华国产精华液的使用体验 | 91aial.com中文字幕在线观看| 欧美区成人在线视频| 国产老妇女一区| 在线观看66精品国产| 亚洲最大成人中文| 网址你懂的国产日韩在线| av在线亚洲专区| 日日摸夜夜添夜夜添av毛片| 看黄色毛片网站| 国产午夜精品论理片| 国产毛片a区久久久久| kizo精华| 日韩在线高清观看一区二区三区| 国产欧美日韩精品一区二区| 日韩大尺度精品在线看网址| 亚洲在线观看片| 乱码一卡2卡4卡精品| 日本三级黄在线观看| 九九爱精品视频在线观看| 校园春色视频在线观看| 校园春色视频在线观看| 国产单亲对白刺激| 国产一级毛片在线| 一级av片app| 日本三级黄在线观看| 久久人妻av系列| 亚洲三级黄色毛片| 啦啦啦观看免费观看视频高清| 中文字幕久久专区| 欧美潮喷喷水| 亚洲成av人片在线播放无| 久久亚洲国产成人精品v| 亚洲国产日韩欧美精品在线观看| 亚洲无线观看免费| 99久久久亚洲精品蜜臀av| 欧美最新免费一区二区三区| 国产精品日韩av在线免费观看| 色吧在线观看| 久久久成人免费电影| 日韩精品有码人妻一区| 黄片wwwwww| 男女做爰动态图高潮gif福利片| а√天堂www在线а√下载| 晚上一个人看的免费电影| 伦精品一区二区三区| 欧美日韩国产亚洲二区| 日韩在线高清观看一区二区三区| 欧美变态另类bdsm刘玥| 精品久久久久久久人妻蜜臀av| 99精品在免费线老司机午夜| 久久久欧美国产精品| 看片在线看免费视频| 久久久欧美国产精品| 日韩在线高清观看一区二区三区| 在线播放无遮挡| 在线播放无遮挡| 伦精品一区二区三区| 精品少妇黑人巨大在线播放 | 人妻久久中文字幕网| 亚洲精品456在线播放app| 一卡2卡三卡四卡精品乱码亚洲| 国产91av在线免费观看| 蜜桃久久精品国产亚洲av| 蜜桃久久精品国产亚洲av| 欧美高清性xxxxhd video| 99久久人妻综合| 在线播放国产精品三级| 又爽又黄a免费视频| 成人永久免费在线观看视频| 成人特级黄色片久久久久久久| 免费观看在线日韩| 日本-黄色视频高清免费观看| 国产免费一级a男人的天堂| 91久久精品电影网| 国产乱人视频| 2021天堂中文幕一二区在线观| 看非洲黑人一级黄片| 一区福利在线观看| 久久精品国产亚洲av香蕉五月| 久久精品国产清高在天天线| 一级毛片电影观看 | 国模一区二区三区四区视频| 国产精品永久免费网站| 在线观看免费视频日本深夜| 国产精品国产三级国产av玫瑰| 亚洲四区av| 欧美最黄视频在线播放免费| 九九在线视频观看精品| 2022亚洲国产成人精品| 国产成人91sexporn| 久久久久久伊人网av| 天堂√8在线中文| 禁无遮挡网站| 免费人成在线观看视频色| 神马国产精品三级电影在线观看| 一边亲一边摸免费视频| av在线天堂中文字幕| 韩国av在线不卡| 人人妻人人看人人澡| 久久热精品热| 黄色视频,在线免费观看| 一边摸一边抽搐一进一小说| 婷婷色av中文字幕| 久久久久久久久中文| 国产午夜精品论理片| 黄色一级大片看看| 亚洲人成网站在线播放欧美日韩| 欧美一区二区精品小视频在线| 日本与韩国留学比较| 看非洲黑人一级黄片| av.在线天堂| 久久精品91蜜桃| 亚洲欧美精品自产自拍| 国产精品1区2区在线观看.| www.av在线官网国产| 久久久久国产网址| 国产日韩欧美在线精品| 好男人在线观看高清免费视频| 最近2019中文字幕mv第一页| 日本一二三区视频观看| av在线天堂中文字幕| 久久午夜福利片| 国产亚洲av嫩草精品影院| 亚洲av一区综合| 青青草视频在线视频观看| 精品久久国产蜜桃| 观看美女的网站| 久久九九热精品免费| 简卡轻食公司| 女人被狂操c到高潮| 亚洲精品456在线播放app| 亚洲欧美日韩东京热| 国内揄拍国产精品人妻在线| 青青草视频在线视频观看| 午夜老司机福利剧场| 啦啦啦啦在线视频资源| 国产高清三级在线| 国产成人精品久久久久久| 日本-黄色视频高清免费观看| 国产成人精品久久久久久| 婷婷六月久久综合丁香| 欧美性猛交黑人性爽| 性插视频无遮挡在线免费观看| 51国产日韩欧美| 国产乱人偷精品视频| 欧美变态另类bdsm刘玥| 身体一侧抽搐| 亚洲国产欧洲综合997久久,| 国产精品一区www在线观看| 高清在线视频一区二区三区 | 亚洲精品国产av成人精品| 欧美成人精品欧美一级黄| 国产欧美日韩精品一区二区| 黄色一级大片看看| av福利片在线观看| 天堂√8在线中文| 1024手机看黄色片| 国产精品电影一区二区三区| 九草在线视频观看| 国产三级中文精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本五十路高清| 国产黄a三级三级三级人| 亚洲人与动物交配视频| 男女那种视频在线观看| 婷婷精品国产亚洲av| 男人舔奶头视频| 在线免费观看的www视频| 高清毛片免费观看视频网站| 国产麻豆成人av免费视频| 国产探花在线观看一区二区| 伦精品一区二区三区| 自拍偷自拍亚洲精品老妇| 精品免费久久久久久久清纯| 国产亚洲欧美98| 国产精品永久免费网站| 国产激情偷乱视频一区二区| 不卡视频在线观看欧美| 人妻夜夜爽99麻豆av| 国产精品人妻久久久影院| 午夜亚洲福利在线播放| 毛片一级片免费看久久久久| 亚洲人与动物交配视频| 老师上课跳d突然被开到最大视频| av天堂中文字幕网| 亚洲精品日韩av片在线观看| 18禁在线播放成人免费| 三级毛片av免费| 国产大屁股一区二区在线视频| 日本-黄色视频高清免费观看| 九九久久精品国产亚洲av麻豆| 日韩欧美精品免费久久| 啦啦啦观看免费观看视频高清| 日韩av不卡免费在线播放| 天堂√8在线中文| 91久久精品国产一区二区成人| 精品人妻熟女av久视频| 久久精品国产亚洲av天美| 可以在线观看毛片的网站| 国内精品一区二区在线观看| 22中文网久久字幕| 久久婷婷人人爽人人干人人爱| 中文字幕熟女人妻在线| 国产成人影院久久av| 男人狂女人下面高潮的视频| 99热精品在线国产| 国产精品一区二区三区四区免费观看| 2022亚洲国产成人精品| 激情 狠狠 欧美| 免费看日本二区| 成人高潮视频无遮挡免费网站| avwww免费| 我要看日韩黄色一级片| 亚洲欧美清纯卡通| 国产精品伦人一区二区| 国产av麻豆久久久久久久| 一级黄色大片毛片| 自拍偷自拍亚洲精品老妇| 国产成年人精品一区二区| 免费人成视频x8x8入口观看| 69av精品久久久久久| 午夜爱爱视频在线播放| 观看免费一级毛片| 看十八女毛片水多多多| 亚洲成人精品中文字幕电影| 国产国拍精品亚洲av在线观看| 最后的刺客免费高清国语| 成熟少妇高潮喷水视频| 老师上课跳d突然被开到最大视频| 黄色欧美视频在线观看| 九九在线视频观看精品| 亚洲欧洲国产日韩| 国产黄片视频在线免费观看| 最近的中文字幕免费完整| 国产高清三级在线| 波野结衣二区三区在线| 五月伊人婷婷丁香| 美女xxoo啪啪120秒动态图| 热99re8久久精品国产| 国产精品电影一区二区三区| 免费观看在线日韩| 久久久久久伊人网av| 嫩草影院入口| 亚洲欧美精品综合久久99| 尾随美女入室| 午夜福利在线在线| 日韩三级伦理在线观看| 一本久久精品| 级片在线观看| 久久久久久久久久久免费av| 欧美日韩在线观看h| 噜噜噜噜噜久久久久久91| 亚洲欧美中文字幕日韩二区| 国产成人a区在线观看| 亚洲av二区三区四区| 少妇人妻一区二区三区视频| 欧美日韩乱码在线| 精品无人区乱码1区二区| 免费人成在线观看视频色| 永久网站在线| 美女黄网站色视频| 97热精品久久久久久| 蜜桃久久精品国产亚洲av| 天天一区二区日本电影三级| 久久99蜜桃精品久久| 噜噜噜噜噜久久久久久91| 国产淫片久久久久久久久| 国产亚洲av片在线观看秒播厂 | www.色视频.com| 欧美在线一区亚洲| 久久草成人影院| 熟妇人妻久久中文字幕3abv| 久久99精品国语久久久| 国产视频首页在线观看| 久99久视频精品免费| 国产日韩欧美在线精品| 亚洲国产高清在线一区二区三| av免费在线看不卡| 毛片一级片免费看久久久久| 亚洲人成网站在线播| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产高清在线一区二区三| 亚州av有码| 好男人在线观看高清免费视频| 精品久久久噜噜| 边亲边吃奶的免费视频| 麻豆国产97在线/欧美| 国产美女午夜福利| 国产激情偷乱视频一区二区| 91午夜精品亚洲一区二区三区| 少妇的逼好多水| 国产精品美女特级片免费视频播放器| 亚洲丝袜综合中文字幕| 久久精品人妻少妇| av在线天堂中文字幕| 一级av片app| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 欧美成人一区二区免费高清观看| 久久久久久久久久久丰满| 国产高清不卡午夜福利| 爱豆传媒免费全集在线观看| 国产精华一区二区三区| 99热精品在线国产| 日韩精品青青久久久久久| 美女xxoo啪啪120秒动态图| 嘟嘟电影网在线观看| 老熟妇乱子伦视频在线观看| 九九久久精品国产亚洲av麻豆| 最近手机中文字幕大全| 青春草视频在线免费观看| 亚洲激情五月婷婷啪啪| 99九九线精品视频在线观看视频| 91精品国产九色| 免费无遮挡裸体视频| 高清日韩中文字幕在线| 男人和女人高潮做爰伦理| 日韩一区二区视频免费看| 男人舔奶头视频| 一个人免费在线观看电影| 午夜精品国产一区二区电影 | 国产伦精品一区二区三区四那| 亚洲精品久久久久久婷婷小说 | 日韩欧美精品v在线| 亚洲av第一区精品v没综合| 内地一区二区视频在线| 三级男女做爰猛烈吃奶摸视频| 99热只有精品国产| 成人性生交大片免费视频hd| 亚洲成av人片在线播放无| 97人妻精品一区二区三区麻豆| 99国产精品一区二区蜜桃av| 尤物成人国产欧美一区二区三区| 中文亚洲av片在线观看爽| 国产精品一区二区三区四区久久| 成人性生交大片免费视频hd| 欧美日本视频| 日韩精品有码人妻一区| 男人的好看免费观看在线视频| 成人欧美大片| 欧美日韩在线观看h| 亚洲欧洲日产国产| 久久国内精品自在自线图片| 99久久精品一区二区三区| 神马国产精品三级电影在线观看| 久久久久久大精品| 最后的刺客免费高清国语| 国产精品麻豆人妻色哟哟久久 | 丰满乱子伦码专区| 中出人妻视频一区二区| 久久精品国产亚洲av涩爱 | 99热这里只有是精品在线观看| 国产探花极品一区二区| av免费观看日本| 国产成人影院久久av| 亚洲av中文字字幕乱码综合| 自拍偷自拍亚洲精品老妇| 六月丁香七月| 免费人成在线观看视频色| 亚洲国产精品sss在线观看| 插阴视频在线观看视频| 免费一级毛片在线播放高清视频| av天堂在线播放| 免费搜索国产男女视频| 大又大粗又爽又黄少妇毛片口| 国产午夜精品久久久久久一区二区三区| 久久6这里有精品| 黄色欧美视频在线观看| 中文欧美无线码| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| avwww免费| 亚洲国产色片| 熟妇人妻久久中文字幕3abv| 久久中文看片网| 青春草国产在线视频 | av国产免费在线观看| 春色校园在线视频观看| 欧美日本视频| 99热网站在线观看| 国产麻豆成人av免费视频| 伦理电影大哥的女人| 国产高清视频在线观看网站| 我的女老师完整版在线观看| 免费看光身美女| 欧美一区二区精品小视频在线| av视频在线观看入口| 99精品在免费线老司机午夜| av免费观看日本| 好男人在线观看高清免费视频| 91狼人影院| 国产精品,欧美在线| 成年免费大片在线观看| 高清午夜精品一区二区三区 | 亚洲精品456在线播放app| 亚洲中文字幕一区二区三区有码在线看| 国产精品女同一区二区软件| 丝袜喷水一区| 久久99蜜桃精品久久| 男插女下体视频免费在线播放| 亚洲性久久影院| 日韩欧美三级三区| 色视频www国产| 国产精品一区二区性色av| 亚洲欧美日韩卡通动漫| 美女被艹到高潮喷水动态| 久久亚洲精品不卡| 久久人人爽人人片av| 夜夜看夜夜爽夜夜摸| 亚洲av一区综合| 国产色婷婷99| 欧美bdsm另类| 非洲黑人性xxxx精品又粗又长| 精品不卡国产一区二区三区| 偷拍熟女少妇极品色| 成人综合一区亚洲| 久久精品国产自在天天线| 国产成人aa在线观看| 欧美日韩综合久久久久久| 免费大片18禁| 久久久精品大字幕| 男女下面进入的视频免费午夜| 精品无人区乱码1区二区| 成人特级黄色片久久久久久久| 69人妻影院| 精品久久久久久久人妻蜜臀av| 我的女老师完整版在线观看| 晚上一个人看的免费电影| 久久热精品热| 美女cb高潮喷水在线观看| 岛国毛片在线播放| 亚洲精品国产成人久久av| 丝袜喷水一区| avwww免费| 桃色一区二区三区在线观看| 成人毛片60女人毛片免费| 性插视频无遮挡在线免费观看| 91精品一卡2卡3卡4卡| 亚洲av免费在线观看| 精品99又大又爽又粗少妇毛片| 久久精品人妻少妇| 男人的好看免费观看在线视频| 国产成人freesex在线| 精品久久国产蜜桃| 国产亚洲av片在线观看秒播厂 | 少妇裸体淫交视频免费看高清| 国产乱人偷精品视频| 久久国产乱子免费精品| 亚洲成人久久爱视频| av.在线天堂| 嫩草影院入口| 伊人久久精品亚洲午夜| 天堂√8在线中文| 国产大屁股一区二区在线视频| 国产午夜精品久久久久久一区二区三区| 在线免费十八禁| av天堂中文字幕网| 亚洲av第一区精品v没综合| 国产精品野战在线观看| 亚洲一级一片aⅴ在线观看| 91精品国产九色| 国产三级在线视频| 国产成人aa在线观看| 精品熟女少妇av免费看| 搡老妇女老女人老熟妇| 久久久欧美国产精品| 天堂√8在线中文| 国产 一区 欧美 日韩| 热99在线观看视频| 91精品国产九色| 久久人人精品亚洲av| 三级毛片av免费| 12—13女人毛片做爰片一| 亚洲欧美清纯卡通| h日本视频在线播放| 亚洲av一区综合| 中文资源天堂在线| 国产成人福利小说| 国产又黄又爽又无遮挡在线| 国产成人freesex在线| 欧美性感艳星| 国产av麻豆久久久久久久| 亚洲自偷自拍三级| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品国产av成人精品| 亚洲性久久影院| 偷拍熟女少妇极品色| 尾随美女入室| 乱人视频在线观看| 日韩在线高清观看一区二区三区| 国产成年人精品一区二区| av在线老鸭窝| 久久久久久久久大av| 九九爱精品视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美另类亚洲清纯唯美| 欧美bdsm另类| 精品99又大又爽又粗少妇毛片| 国产不卡一卡二| 亚洲久久久久久中文字幕| 久久99精品国语久久久| 日本黄色片子视频| 亚洲国产精品sss在线观看| 国产单亲对白刺激| 精华霜和精华液先用哪个| 老司机福利观看| 午夜a级毛片| 蜜桃久久精品国产亚洲av| 久久这里只有精品中国| 久久精品国产鲁丝片午夜精品| 亚洲综合色惰| 又黄又爽又刺激的免费视频.| 五月玫瑰六月丁香| 99热精品在线国产| 国产一区二区亚洲精品在线观看| 少妇的逼水好多| 欧美性猛交╳xxx乱大交人| 欧美+日韩+精品| 成人漫画全彩无遮挡| av在线亚洲专区| 狠狠狠狠99中文字幕| 久久久久久久久久久免费av| 欧美xxxx黑人xx丫x性爽| 亚洲中文字幕日韩| 日韩亚洲欧美综合| 国产爱豆传媒在线观看| 久久精品国产亚洲网站| 欧美日韩在线观看h| 成人亚洲精品av一区二区| 国产伦精品一区二区三区视频9| 亚洲欧美精品综合久久99| 亚洲欧洲国产日韩| 天堂影院成人在线观看| 国产精品人妻久久久影院| 国产欧美日韩精品一区二区| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 校园人妻丝袜中文字幕| 久久久久国产网址| 好男人视频免费观看在线| 你懂的网址亚洲精品在线观看 | 伦理电影大哥的女人| videossex国产| 久久精品91蜜桃| 欧美xxxx性猛交bbbb|