999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

永磁同步直線電機的小波神經網絡控制

2013-09-20 05:50:10黨選舉,徐小平,于曉明
電機與控制學報 2013年1期
關鍵詞:方法模型系統

0 引言

與旋轉電機相比,永磁同步直線電機(permanent magnet linear synchronous motor,PMLSM)不需要任何中間轉換機構,直接將電能轉化為機械能,因而具有結構簡單、高精確度、高速度和低維護等優點,在機器人、半導體制造、高精度數控機床和X-Y驅動設備等領域有著廣闊的應用前景[1-2]。然而,正是由于其直線結構特點和缺少中間緩沖環節,各種干擾,如推力波動、摩擦力等直接作用于PMLSM,使得PMLSM對干擾更為敏感。在使用PMLSM驅動的高精度數控機床中,干擾會使系統的伺服性能降低,產生振蕩增加能量損耗,導致加工的產品表面出現劃痕,甚至使產品報廢[3]。因此,與旋轉電機相比,PMLSM對干擾抑制的質量要求更高,難度也更大。抑制干擾對PMLSM系統的影響,國內外學者進行了大量的研究工作。復合前饋 PID控制方法[4],采用固定參數模型對PMLSM進行速度、加速度補償,可顯著提高系統的跟蹤精度,但該方法的跟蹤效果受模型參數估計誤差影響較大,魯棒性欠佳,且不能消除干擾影響,其跟蹤精度也有限。神經網絡自適應逆模型控制方法[5],采用神經網絡構建控制對象的逆模型,具有較好的跟蹤性和魯棒性,但該方法中神經網絡把PMLSM與干擾作為一個整體進行學習時,往往會忽略對高頻干擾部分的補償,實際上,主要還是對PMLSM系統進行了速度、加速度補償,系統依然會受干擾的影響。自適應前饋控制與PID控制相結合的控制方法[6],通過對誤差信號作頻譜分析得出推力波動所包含的主要諧波,在前饋補償中植入“木馬”的方法來減少推力波動對系統影響,具有較好的控制效果,但該方法采用離線方式對誤差信號進行傅里葉變換分析時,只能得到干擾的頻域信息,而無任何時域信息,因此并不能完全地實現對干擾的在線估計與補償。

文獻[7-8]在建立直線電機干擾觀測器模型的基礎上,提出了一種有效消除干擾的方法,較好地實現了對直線電機位置的精密控制。

該文采用“分而治之”的思想,在復合前饋PID控制的基礎上,提出了一種基于小波神經網絡的控制方法:1)用固定參數模型對PMLSM系統進行速度、加速度補償;2)推力波動、摩擦力等干擾及固定參數模型不精確所造成的誤差主要由小波神經網絡在線估計進行補償。所提出的基于小波神經網絡控制方法有效地提高了PMLSM系統的跟蹤性及魯棒性,較好地消除了干擾的影響。

1 PMLSM系統模型

PMLSM的干擾主要包括由端部效應、齒槽效應、初級電流和反電動勢存在高次諧波及氣隙磁密波形非正弦性所引起的推力波動和摩擦力等[9-10]。

考慮干擾影響,PMLSM的電壓及動力學方程[11]可表示為

其中:uq(t)、iq(t)分別為PMLSM的q軸驅動電壓及q軸電樞電流;f(t)為電磁推力;fr(x) 、ff()為推力波動及摩擦力;fΔ為不確定性干擾,如電機參數擾動和系統噪音所造成的影響等;Ke、R、L、Kf、M、x及x·分別為電機的反電動勢常數、電阻、電感、推力常數、動子質量、電機位移及速度。

推力波動fr(x)是與電機位移x有關的周期性波動函數,其模型[12]可描述為

式中:ai為推力波動的幅值;ωi為以電機位移為變量的角速度;φi為相位角。

摩擦力ff()的模型[12]可描述為

其中:fm為靜態摩擦力;fc為庫侖摩擦力為動子速度;fv為粘滯摩擦系數為經驗系數。

此外,式(1)中,相對于電阻R值,電感L值很小,參考文獻[6]、[11]與[12],文獻將其忽略不計,則PMLSM的動力學方程可表示為

令a=則PMLSM系統模型如圖1所示,表現為一個強非線性、多變量的伺服系統。

圖1 PMLSM系統模型Fig.1 System model of PMLSM

2 基于小波神經網絡的PMLSM控制系統設計

2.1 復合前饋控制

在高精度伺服系統中,復合控制系統可極大地提高系統的跟蹤精度。其基本思想為:在系統中引入前饋環節,使之與閉環系統的傳遞函數之積為1,從而使系統的輸出能夠完全復現輸入[13],其結構如圖2所示。

圖2 復合前饋控制結構Fig.2 Frame diagram of a feedforward composite control

在圖2中,當Gf(s)=1/G(s)時,電機輸出位移x(s)可復現位置輸入信號x*(s),達到理想的跟蹤性能。但實際系統中的模型參數的精確值是不可能得到的,且系統仍然不可避免地受到干擾的影響。

2.2 基于小波神經網絡的PMLSM控制系統

在復合前饋PID控制的基礎上,利用小波函數對非平穩信號和非線性函數逐步精細描述的強逼近特點,采用小波神經網絡在線估計、補償推力波動、摩擦力等干擾以及補償前饋環節中,由于采用固定參數模型的不精確所造成的誤差。該文所設計的PMLSM控制系統結構如圖3所示。

由圖3設計控制系統可得,控制電壓u=up+uf+ud。up為PID控制器的輸出電壓值,為系統的反饋誤差控制器,e為系統的位置跟蹤誤差,可得

其中,P、TI和TD分別為比例系數、積分時間常數和微分時間常數。

圖3 PMLSM控制系統結構Fig.3 Frame diagram of control for PMLSM system

uf為根據電機逆模型,采用固定參數進行速度、加速度補償電壓值,使系統的跟蹤精度得到顯著提高,其中,a*和b*值可根據電機給定的額定參數計算獲得,也可通過實驗辨識獲得,其輸出值為

由式(8)可得,控制系統適應于連續可微的位置輸入跟蹤信號,如正弦波等,如要跟蹤不可微信號,如三角波,可用不完全微分法近似代替。

根據PMLSM的動力學方程,如式(6)所示,當uf為系統進行速度、加速度補償后,系統的跟蹤精確度可得到顯著提高,則影響系統跟蹤性能的主要因素集中于推力波動、摩擦力等干擾以及固定參數模型不精確造成的擾動等,可由小波神經網絡WNN的輸出值ud進行補償,其學習信號為反饋誤差控制器PID的輸出值up,輸入信號為電機位移x,目的是使得系統的反饋誤差值趨于一個較小值,從而達到進一步提高系統跟蹤精度目的。TDL(Tapped Delay Line)為多分頭時延單元。因此,在該系統中,WNN的作用主要有兩個:1)估計推力波動、摩擦力等干擾,實現在線補償;2)補償前饋環節中固定參數模型不精確所造成的誤差。

3 小波神經網絡的綜合問題

3.1 小波函數的選擇及對主要干擾的逼近性能分析

小波神經網絡中小波函數的選擇,目前仍無統一的理論確定哪類小波函數適合做不同網絡中神經元的激勵函數,主要根據經驗和實際情況。該文選取了在函數逼近方面應用較多的一維墨西哥草帽函數[14-15],即

證明:Ψ(x)=(1 -x2)的傅里葉變換為

顯然此函數滿足

則由式(11)可得Ψ(x)滿足母小波函數條件。令

其中,a和b分別為伸縮參數和平移參數。借鑒文獻[16]中定理2,通過適當的選擇a和b,{Ψmn(x)}亦可以構成L2(R)上的框架,即

PMLSM中的推力波動為連續函數,故采用一維墨西哥草帽函數構成的小波神經網絡可以任意精度實現對推力波動信號的逼近。

證畢。

考察式(5)摩擦模型可得,摩擦力在速度x·=0是不連續的。根據定理1可得,采用一維墨西哥草帽函數構成的小波神經網絡可任意精確度逼近摩擦力的連續部分,但對不連續部分則無法有效逼近,其誤差主要通過反饋誤差控制環節來抑制。

3.2 小波神經網絡的結構及學習算法

與徑向基函數神經網絡結構相同,小波神經網絡結構采用3層前向網絡。其結構如圖4所示。

圖4 小波神經網絡結構Fig.4 Frame diagram of wavelet neural network

圖4 中,輸入變量 X=[x1,x2,…,xn]T為電機的位移量,隱含層中的基向量為 ψ =[ψ1,ψ2,…ψm]T,激勵函數使用一維墨西哥草帽小波函數得

式中:bj=[bj1,bj2,…,bjn]T為網絡節點j的移位參數;aj=[a1,a2,…,am]T為節點j的伸縮參數,網絡輸出層的權向量為 Wm=[ω1,ω2,…,ωj,…,ωm]T。

根據圖3中設計的結構圖,WNN的輸入變量為X=[x(k-1),x(k-2),…,x(k-n)]T,網絡輸出為

學習信號為反饋誤差控制器PID的輸出值,即

根據反饋誤差學習方法[17-18],可避免使用 Jacobian信息。同時根據梯度下降法,輸出權值、伸縮參數及平移參數的迭代算法分別如下(η為學習速率,α為動量因子),即

3.3 隱層節點數的確定及參數初始化

與小波函數的選擇相同,小波神經網絡中隱層節點數的確定和參數初始化尚無統一的理論參考。從理論上講,隱層節點數越多,網絡的逼近精度越高,但網絡的結構也越復雜。因此,需要在逼近精度和結構復雜度之間進行折中[16]。參數初始化主要是對伸縮參數、平移參數和權值進行初始化,其大小直接影響網絡的收斂性和逼近精度。文獻[16]提出伸縮參數初始化主要考慮處理函數的極值,平移參數主要取函數時域邊界之間的均值,權值初始化則同一般神經網絡完全一樣,其收斂速度較快。實際上,只是憑借經驗判斷,能夠保證網絡收斂,但并不一定能夠保證網絡達到較好的逼近性能。

論文采用K-均值聚類算法對網絡參數進行初始化訓練,可較好地解決網絡結構對初值敏感問題,并保證網絡達到較好的逼近性能。具體算法步驟可參考文獻[19]。

4 實驗結果及分析

4.1 實驗平臺簡介

PMLSM實驗平臺主要包括某微納科技有限公司自主研制的無鐵芯永磁同步直線電機平臺、控制箱及PC機組成。

圖5 控制系統的硬件原理Fig.5 Hardware block diagram of control system

系統的硬件原理框圖如圖5所示。硬件實物圖如圖6所示。在圖6中,控制箱中的驅動器采用的是以色列 Elmo公司生產的 Harmonica驅動器。PMLSM的主要額定參數為:動子質量M=0.3 kg,線圈電阻R=19.3 Ω,推力常數Kf=11.71 N/A,反電動勢常數Ke=9.6 V·m/s,電感L=2.49 mH,最大運動行程為370 mm。

圖6 控制系統實驗平臺Fig.6 Platform of system experiment

軟件資源方面主要是基于TMS320F2812開發的 cSPACE 控制系統,具有 A/D、D/A、I/O、Encoder和快速控制原型開發和硬件在環仿真功能。設計過程中,只需使用Matlab/SIMULINK環境自動生成代碼便可完成軟件設計,不需要另外編寫代碼。因此,該平臺具有開放性、快速性和易與Matlab聯接等諸多優點。

4.2 基于小波神經網絡的控制方法實現

基于小波神經網絡的控制方法在Matlab/SIMULINK環境下搭建的cSPACE控制算法及運行結果如圖7、8所示。圖7中小波神經網絡補償器NNC采用了效率高、速度快的C-MEX S-function編程,該算法直接編譯下載到cSPACE系統后,便可實現實時控制。

圖7 基于小波神經網絡的cSPACE控制算法Fig.7 cSPACE control algorithm based on wavelet neural network

圖8 基于小波神經網絡的控制方法運行結果Fig.8 Rresults of control based on wavelet neural network

在圖7中,輸入信號為正弦信號,頻率為1Hz,“WM-Write1”用于設定正弦信號的幅值,與圖8中“Write1”對應,設置為20 mm;圖7中“In2”與“WM-Write3”分別為小波神經網絡的學習速率與增益系數,與圖8中“Write2”與“Write3”對應,分別設置為0.45與 0.005;圖 7中“WM -Write4”、“WM -Write5”與“WM-Write6”分別 PID控制器中的比例、微分及積分系數,與圖8中“Write4”、“Write5”與“Write6”對應,分別設置為0.50、0.005與0;圖7中“WM -Write7”與“WM -Write8”分別為系統的速度與加速度補償系數,與圖 8中“Write7”與“Write8”對應,根據PMLSM的額定參數的計算值,設置為0.009 6與0.000 49。

在圖8運行結果圖中,左上角第一個窗口顯示的是位置參考輸入信號,與圖7中“WM-Read1”對應;右上角第二個窗口顯示的是電機實際位移,與圖7中“WM-Read2”對應;左下角第三個窗口顯示的是電機位置跟蹤誤差,與圖7中“WM-Read3”對應;右下角第四個窗口顯示的是小波神經網絡的補償電壓值,與圖7中“WM-Read4”對應。由圖8可得,該方法下的位置跟蹤誤差約為0.15 mm,誤差百分比為0.75%。

圖9與圖10為采用復合前饋PID控制方法下的cSPACE控制算法及運行結果圖。

圖9 復合前饋PID控制的cSPACE控制算法Fig.9 cSPACE control algorithm of combined feedforward plus PID control

圖10 復合前饋PID控制方法運行結果Fig.10 Results of combined feedforward plus PID control

調整PID控制器中比例、微分與積分系數為0.65、0.003 5與 0.10,速度、加速度補償系數為0.009 6與0.000 49。由圖10運行結果可知,采用復合前饋PID控制方法下的位置跟蹤誤差約為0.44mm,誤差百分比為2.2%。

圖11與圖12給出了采用小波神經網絡逆模型控制方法下的cSPACE控制算法及運行結果圖。

圖11 小波神經網絡逆模型控制的cSPACE控制算法Fig.11 cSPACE control algorithm of wavelet neural network adaptive inverse control

圖12 小波神經網絡逆模型控制運行結果Fig.12 Results of wavelet neural network adaptive inverse control

調整PID控制器中比例、微分與積分系數為0.35、0.004與0,小波神經網絡學習速率與增益系數為0.06與0.015。由圖12運行結果可得,采用小波神經網絡逆模型控制方法下的位置跟蹤誤差約為0.56 mm,誤差百分比為2.8%。

為了便于觀察干擾消除情況,對圖8中基于小波神經網絡控制方法與圖12中小波神經網絡逆模型控制方法中的位置跟蹤誤差信號,分別觀察其頻譜圖,其結果如圖13所示。表1給出了3種控制方法下跟蹤性能的匯總結果。

圖13 位置跟蹤誤差的頻譜分析比較Fig.13 Comparisons of frequency analysis for position tracking errors

表1 3種控制方法的跟蹤性能比較Table 1 Comparisons of tracking performance

比較圖7~13及表1數據可知,采用復合前饋PID控制方法對PMLSM系統的補償量只是根據PMLSM模型本身進行了近似補償,并沒有對干擾進行補償,且其跟蹤效果受參數估計值精度影響較大;小波神經網絡自適應逆模型控制方法不需要估計模型參數值,類似于無模型控制思想,可有效應對模型參數變化,魯棒性較好,但經頻譜分析可知,此方法不能有效消除干擾影響。基于小波神經網絡的控制方法采用的是“分而治之”思想,雖然在對PMLSM系統進行速度、加速度補償時需要對模型參數值進行估計,但參數不精確造成的誤差及推力波動、摩擦力等干擾影響可通過小波神經網絡在線估計進行補償,因此受參數估計值影響較小。較前兩種控制方法,基于小波神經網絡控制的位置跟蹤精度能夠得到顯著提高,且經頻譜分析可得該方法能較好地消除干擾影響,從而驗證了基于小波神經網絡控制的有效性。

5 結論

1)采用“分而治之”思想,在復合前饋PID控制的基礎上,引入小波神經網絡在線補償推力波動、摩擦力以及固定參數模型不精確等因素對系統造成的干擾,可使PMLSM系統獲得理想的跟蹤效果,且能夠較好地消除干擾影響。

2)小波神經網絡對摩擦力中非連續部分估計效果不佳,導致PMLSM系統仍然受到靜摩擦力影響,在后續工作中可加入簡單的靜摩擦力補償環節達到完善提高。

3)提出的控制系統適應于連續可微的位置輸入跟蹤信號。如正弦波等,如要跟蹤不可微信號,如三角波,可用不完全微分法近似代替。

4)對小波神經網絡中權值參數的初始化問題,仍需作進一步研究。

[1]王麗梅,武志濤,左濤.永磁直線電機自構式模糊神經網絡控制器設計[J].電機與控制學報,2009,13(5):643 -647.

WANG Limei,WU Zhitao,ZUO Tao.PMLSM self-constructing fuzzy neural network controller design[J].Electric Machines and Control,2009,13(5):643 - 647.

[2]李騰,劉延杰,孫立寧,等.考慮端部效應的永磁直線同步電機鍵合圖建模[J].電機與控制學報,2011,15(2):71 -75.

LI Teng,LIU Yanjie,SUN Lining,et al.Bond graph approach to the modeling of permanent magnet linear synchronous motor with consideration of end-effect[J].Electric Machines and Control,2011,15(2):71 -75.

[3]郭慶鼎,王成元.直線交流伺服系統的精密控制[M].北京:機械工業出版社,2000:50-72.

[4]侯伯杰,李小清,周云飛,等.直線電機伺服系統的復合前饋 P I D 控制[J].機床與液壓,2009,37(2):56 -61.

H0U Bojie,LI Xiaoqing,ZH0U Yunfei,et al.The development of feedforward plus PID controller for linear motor[J].Machine Tool and Hydraulics,2009,37(2):56 -61.

[5]夏加寬,王成元,李皞東,等.高精度數控機床用直線電機端部效應分析及神經網絡補償技術研究[J].中國電機工程學報,2003,23(8):100 -104.

XIA Jiakuan,WANG Chengyuan,LI Haodong,et al.Study of endeffects analysis and NN compensation technique of linear motor for high precision NC machine[J].Proceedings of the Chinese Society for Electrical Engineering,2003,23(8):100 -104.

[6]ZHAO S,TAN K K.Adaptive feedforward compensation of force ripples in linear motors[J].Control Engineering Practice,2005(13):1081-1092.

[7]張代林,陳幼平,艾武,等.基于觀測器模型的直線電機干擾抑制技術的研究[J].中國電機工程學報,2007,27(12):14 -18.

ZHANG Dailin,CHEN Youping,AI Wu,et al.Research on disturbance suppression technology for linear motors based on a disturbance cbserver model[J].Proceedings of the Chinese Society for Electrical Engineering,2007,27(12):14 -18.

[8]ZHANG Dailin,CHEN Youping,AI Wu,et al.Force ripple suppression technology for linear motors based on back propagation neural network[J].Chinese Journal of Mechanical Engineering,2008,21(2):13 -16.

[9]宋亦旭,王春洪,尹文生,等.永磁直線同步電動機的自適應學習控制[J].中國電機工程學報,2005,25(20):151 -156.

SONG Yixu,WANG Chunhong,YIN Wensheng,et al.Adaptivelearning control for permanent-magnet linear synchronous motors[J].Proceedings of the CSEE,2005,25(20):151 -156.

[10]ZHU Y W,JIN S M,CHUNG K,et al.Control-based reduction of detent force for permanent magnet linear synchronous motor[J].IEEE Transactions on Magnetics,2009,45(6):2827-2830.

[11]CHEN S L,TAN K K,HUANG S N,et al.Modeling and compensation of ripples and friction in permanent-magnet linear motor using a hysteretic relay[J].IEEE Transactions on Mechatronics,2010,15(4):586 -594.

[12]AHN H,CHEN Y Q,DOU H F.State-periodic adaptive compensation of cogging and coulomb friction in permanent-magnet linear motors[J].IEEE Transactions on Magnetics,2005,41(1):90-98.

[13]劉金琨.先進PID控制及其MATLAB仿真[M].北京:電子工業出版社,2004:55-58.

[14]YLMAZ S,OYSAL Y.A fuzzy wavelet neural network model for system identification[C]//2009 Ninth International Conference on Intelligent Systems Design and Applications,2009,Nov 30 -Dec 2,2009,Pisa,Italy.2009:1284-1289.

[15]程啟明,王勇浩.基于小波神經網絡的控制方法及其應用研究[J].工業儀表與自動化裝置,2004,4:6 -9.

CHEN Qiming,WANGYonghao.A study on the control method and its application based on wavelet neural network[J].Industrial Instrumentation & Automation,2004,4:6 -9.

[16]張茁生,劉貴忠,張宗平.基于小波框架的自適應徑向基函數網絡[J].自動化學報,2002,28(2):229 -236.

ZHANG Zhuosheng,LIU Guizhong,ZHANG Zongping.Radial basis function networks based on wavelet frames[J].Acta Automatica Sinica,2002,28(2):229 -236.

[17]阮曉鋼,丁名曉,于乃功,等.加入預測信息的反饋誤差學習模型及其仿真研究[J].系統仿真學報,2006,18(11):3227-3229,3246.

RUAN Xiaogang,DING Mingxiao,YU Raigong,et al.Design and simulation of predictive feedback error learning model[J].Journal of System Simulation,2006,18(11):3227 -3229,3246.

[18]黨選舉.壓電陶瓷執行器的神經網絡實時自適應逆控制[J].光學精密工程,2008,16(7)1266 -1272.

DANG Xuanju.Real-time adaptive inverse control based on neural networks for piezoceramic actuator[J].Optics and Precision Engineering,2008,2008,16(7)1266 -1272.

[19]MACQUEEN J.Some methods for classification and analysis of multivariate observations[C]//Proc of the 5th Berkeley Symposium on Mathematical Statistics and Probability.Berkeley,USA:University of California Press,1967:281 -297.

猜你喜歡
方法模型系統
一半模型
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
3D打印中的模型分割與打包
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
主站蜘蛛池模板: 日本www色视频| 激情综合婷婷丁香五月尤物| 亚洲高清资源| 99视频在线免费| 最新亚洲人成无码网站欣赏网 | 国产精品自在拍首页视频8| 手机精品福利在线观看| 91国内视频在线观看| 欧美一级在线看| 天天摸天天操免费播放小视频| 在线va视频| 国产精品视频白浆免费视频| 国产在线视频自拍| 99在线观看国产| 小说 亚洲 无码 精品| 国产精品美女免费视频大全| 91精品福利自产拍在线观看| 国产成人精品一区二区| 国产91在线|日本| 波多野结衣国产精品| 最新亚洲av女人的天堂| 午夜无码一区二区三区| 无码av免费不卡在线观看| 亚洲无码电影| 久久人人爽人人爽人人片aV东京热| 日韩在线影院| 久草青青在线视频| 国产乱视频网站| 亚洲第一区在线| 国产丝袜精品| 精品自拍视频在线观看| 精品1区2区3区| 久久国产V一级毛多内射| 欧美一区精品| 欧美久久网| 亚洲区一区| 欧美午夜精品| 精品国产免费观看| 免费xxxxx在线观看网站| www.亚洲一区| 最新国产你懂的在线网址| 亚洲AV人人澡人人双人| 人人爱天天做夜夜爽| 成人蜜桃网| 中文字幕免费播放| 亚洲av日韩综合一区尤物| 亚洲男人的天堂在线| 午夜人性色福利无码视频在线观看| 亚洲丝袜中文字幕| 国产精品免费p区| 2021精品国产自在现线看| 亚洲人成人无码www| 嫩草国产在线| 亚洲精品在线91| 色丁丁毛片在线观看| 亚洲成人www| 国产人人乐人人爱| 久久人人97超碰人人澡爱香蕉 | 免费一级全黄少妇性色生活片| 亚洲色图欧美视频| 一区二区在线视频免费观看| 色综合五月| 国模视频一区二区| 亚洲自偷自拍另类小说| 美女视频黄频a免费高清不卡| 欧美激情视频一区二区三区免费| 在线看片中文字幕| 国产特级毛片| 国产亚洲精| 波多野结衣二区| 国产精品久久久免费视频| 日韩AV无码一区| 孕妇高潮太爽了在线观看免费| 亚洲人成网站观看在线观看| 亚洲国产中文精品va在线播放| 97久久精品人人做人人爽| 在线观看国产小视频| 一级片免费网站| 青青青国产精品国产精品美女| 福利在线一区| 日韩在线成年视频人网站观看| 免费人成视网站在线不卡|