亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        小球式旋轉直驅壓力伺服閥動態(tài)特性分析優(yōu)化

        2018-10-30 11:54:34陸亮夏飛燕訚耀保原佳陽方向
        航空學報 2018年10期
        關鍵詞:滑閥階躍轉角

        陸亮,夏飛燕,訚耀保,*,原佳陽,方向

        1. 同濟大學 機械與能源工程學院,上海 200092 2. 南京機電液壓工程研究中心 航空機電系統(tǒng)綜合航空科技重點實驗室,南京 210061

        20世紀80年代末期,電機出現(xiàn)了混合勵磁同步結構,其輸出效率高、力矩大,應用廣泛[1-3]?;诟吖β拭芏入姍C的發(fā)展,同期,出現(xiàn)了一批旋轉直接驅動伺服閥的專利[4-6],通過運動轉換機構將電機的旋轉運動轉換成閥芯的直線運動,其目的是取消傳統(tǒng)兩級閥中復雜的前置級。直接驅動電液壓力伺服閥擁有結構簡單、抗污染能力強等優(yōu)勢,可應用在緊湊型液壓系統(tǒng)中替代傳統(tǒng)的兩級電液伺服閥,其應用優(yōu)勢逐漸體現(xiàn)出來,如在飛機電子防滑剎車系統(tǒng)中取代傳統(tǒng)的噴嘴擋板閥作為壓力控制元件,避免易發(fā)的油液污染顆粒堵塞噴嘴口而使壓力控制失效的問題。

        近年來,直接驅動的電液壓力伺服閥發(fā)展迅速。1993年日本學者Urai等[7]提出采用磁致伸縮材料提供閥芯直接驅動的動力,通過通電線圈產(chǎn)生磁場,處于磁場中的磁致伸縮材料發(fā)生變形,進而推動閥芯,此種結構的閥結構緊湊,階躍響應快,但對溫度變化敏感。MOOG公司提出以直線力矩馬達直接驅動閥芯的直動閥結構,如D633/644系列[8],該系列現(xiàn)已形成產(chǎn)品,其動靜態(tài)特性良好,但它將直線力矩馬達與閥芯同軸布置,結構不夠緊湊。2006年以PARK公司為首的音圈驅動[9]單級閥開始發(fā)展,將閥芯與通電線圈相連,在外加磁場作用下,線圈帶動閥芯移動,同時在線圈上內(nèi)置位移傳感器,用于檢測閥芯位置,由于線圈無支承,其慣量小、無摩擦,整閥響應線性度高,但是由于線圈位移有限,其可控的輸出流量或壓力較小。2016年西安交通大學[10]提出一種轉閥結構,在閥芯和閥套上對應開槽,通過閥芯與閥套的相互轉動開啟或閉合對應的通槽,從而控制輸出的流量或壓力。這種結構使得污染顆粒難以在閥內(nèi)堆積,抗污染能力強,且閥芯轉動阻尼相對滑動阻尼小,動態(tài)特性較優(yōu)。

        伺服閥作為飛機剎車系統(tǒng)中的壓力控制元件,其動態(tài)特性的好壞,直接影響剎車系統(tǒng)的控制效果[11]。2002年Yuan和Li[12-13]先后提出通過減小滑閥阻尼長度的方法來增加其響應速度,通過在滑閥節(jié)流邊開槽,將滑閥變?yōu)樨撟枘峥梢杂行У馗纳苹y的動態(tài)特性。2009年吳泊寧等[14]根據(jù)模態(tài)截斷的相關理論,推導了閉環(huán)系統(tǒng)在不同模態(tài)疊加過程中的穩(wěn)定性變化,得出了伺服閥系統(tǒng)在高頻階段失穩(wěn)的結論。2014年黃澄等[15]通過時域仿真結合頻域分析證明了飛機液壓剎車系統(tǒng)管路會影響壓力伺服閥的輸出特性,兩者相耦合會引發(fā)伺服閥振蕩失穩(wěn)。2016年逯九利等[16]利用AMESim仿真驗證了飛機剎車壓力控制系統(tǒng)諧振的原因,并優(yōu)化選取了管路、系統(tǒng)背壓及壓力控制閥參數(shù)。從系統(tǒng)控制的角度來看,壓力伺服閥作為控制系統(tǒng)關鍵環(huán)節(jié),在參數(shù)不匹配的情況下,極易形成振蕩、失穩(wěn)問題[17-19]。

        針對一款具有開發(fā)應用價值的緊湊型旋轉直接驅動電液壓力伺服閥,在動態(tài)特性測試過程中出現(xiàn)的控制壓力超調量大、響應調整時間長等現(xiàn)象,建立了閉環(huán)系統(tǒng)從輸入電信號至輸出壓力的數(shù)學模型,分析傳遞函數(shù)中影響該閥動態(tài)特性的重要結構參數(shù)和電控方法,優(yōu)化選取參數(shù)及改進電控方法后,通過實驗驗證整閥的動態(tài)響應效果。

        1 小球式旋轉直驅閥的工作原理及框圖

        1.1 工作原理

        小球式旋轉直接驅動電液壓力伺服閥(Ball-type Rotary Direct Drive Pressure Servo Valve, BRDDPSV)由有限轉角電機、小球-柱形孔運動副、功率滑閥副及反饋傳感電子電路組成,如圖1所示。

        該旋轉直驅閥采用有限轉角電機驅動,將電機輸出軸上偏置的小球插裝在閥芯末端等徑的柱形孔中,使電機的旋轉運動轉換成閥芯的水平運動和繞軸線的轉動。其中,閥芯的水平運動可控制輸出壓力,而繞軸線的轉動僅為了避免運動干涉。正開口滑閥副為兩凸肩、三通的結構形式,滑閥閥芯兩端接通回油,閥芯左右兩腔壓力p1、p2均為0,負載口置于供油口和回油口之間,負載口中的控制壓力由閥芯滑動位移來調節(jié):滑閥位移在供油口和回油口間形成聯(lián)動的節(jié)流作用,從而使中間的負載口輸出可控的壓力[20]。同時,反饋電機轉角和控制腔中的壓力至輸入信號,構成閉環(huán)控制。

        圖1 BRDDPSV結構示意圖Fig.1 Structure diagram of BRDDPSV

        1.2 數(shù)學模型

        根據(jù)該旋轉直驅壓力伺服閥工作原理,可建立該閥的數(shù)學模型。

        經(jīng)典控制理論為閉環(huán)非線性系統(tǒng)提供了線性化處理方法:若系統(tǒng)中的非線性因素較弱,可直接忽略次要的非線性部分;若系統(tǒng)中的變量只是在工作點附近發(fā)生了微量偏移,可以應用切線法取其中的線性主要部分,即在工作點線性化。線性化后,以成熟的經(jīng)典控制線性理論分析問題。

        1.2.1 電 機

        有限轉角電機采用單向繞組,不需要整流電路,轉角限制在±30°內(nèi)。電機轉子動力學方程為

        (1)

        式中:Jr為電機轉動慣量;α為電機轉角;Tem為電機輸出力矩;Br為轉動阻尼;TL為負載力矩。

        實測電機的輸出電磁力矩Tem和轉角α的關系曲線為拋物線,采用二階函數(shù)模型對該曲線進行擬合得

        Tem=kti0-kmα2

        (2)

        式中:kt為電流-力矩系數(shù);km為轉角-力矩系數(shù);i0為電機輸入電流。

        在電機轉角較小時,輸入電流與輸出力矩近似呈線性關系

        Tem=kti0

        (3)

        1.2.2 小球-柱形孔運動轉換機構

        小球-柱形孔運動轉換機構存在兩處偏心:一是小球球心偏離電機軸線e0;二是小球球心偏離閥芯軸線h0,如圖2所示,建立圖示坐標系。電機繞軸線Y的旋轉運動通過該轉換機構傳遞給閥芯,使閥芯產(chǎn)生兩個方向的運動:一是沿X軸線的水平運動xv;二是繞X軸線的轉動βv。

        由于該壓力伺服閥閥芯位移較小(最大值僅為0.1 mm),所需的電機轉動角度α也較小(小于6°),則有sinα=α,cosα=1-α2/2,驅動接口運動學方程可以簡化為

        xv=e0sinα≈e0α

        (4)

        式中:αx為馬達轉子的任意工作點;kβv為閥芯繞軸線的轉角變化系數(shù)。

        電機所受的阻力矩為

        圖2 運動轉換機構示意圖Fig.2 Structure diagram of movement conversion mechanism

        (5)

        式中:Fx為小球對閥芯X方向的作用力;Tβ為小球對閥芯繞軸線轉動的力矩。

        1.2.3 正開口滑閥

        由圖1可知,滑閥采用二閥肩三通的正開口結構。正常工作時,進、回油口處同時節(jié)流,以控制位于進、回油口中間的控制腔中的壓力,由閥腔內(nèi)液體流動的連續(xù)性可得控制壓力隨閥口開度變化的表達式為

        (6)

        式中:Cd為流量系數(shù);W為節(jié)流窗口面積梯度,W=πD;D為閥肩直徑;ps為系統(tǒng)壓力;pL為整閥控制壓力;U為正開口滑閥預開口量;ρ為油液密度;V為負載容腔體積;E為油液體積彈性模量。

        在負載容腔體積較小(10 mL)時,忽略油液可壓縮性,則控制壓力與閥芯位移的關系可表示為

        (7)

        在任意工作點線性化可表示為

        pL=kpxv

        (8)

        式中:kp為控制壓力增益系數(shù)。

        閥芯在運動過程中受驅動力、穩(wěn)態(tài)液動力作用,正開口滑閥穩(wěn)態(tài)液動力可表示為

        Fs=2CvCdWcosθ[xv(ps-pL)-(U-xv)pL]

        (9)

        式中:Cv為流速系數(shù);θ為射流角。

        穩(wěn)態(tài)液動力與滑閥開度有關,可將其視為彈性力,由式(7)~式(9),有

        (10)

        式中:ks為穩(wěn)態(tài)液動力剛度。

        正開口滑閥所受驅動力和驅動力矩的表達式為

        (11)

        式中:m為閥芯質量;Bv為閥芯沿軸線滑動阻尼系數(shù);Jβv為閥芯繞軸線轉動的轉動慣量;Bβv為閥芯繞軸線轉動阻尼系數(shù)。

        將式(2)~式(11)代入式(1)可得等效在電機轉子上的輸出力矩與電機轉角的關系:

        (12)

        式中:JM、BM、KM為伺服閥中機械運動部件等效到電機轉軸上的轉動慣量、阻尼系數(shù)和彈性系數(shù)。

        1.2.4 電子控制電路

        將實際輸出的控制壓力pL與輸入指令電流i的偏差信號進行比例kp、積分ki運算(PI),輸出電機的控制信號i0。同時小閉環(huán)還反饋電機轉角α以增加電機剛度,提高系統(tǒng)的抗擾動性能。壓力反饋系數(shù)及電機轉角反饋系數(shù)分別為kf1、kf2。

        (13)

        式中:kpwm為功率放大系數(shù);Rc為電機電阻;kb為指令電流放大系數(shù);s為Laplace算子。

        1.3 整閥框圖

        根據(jù)線性化模型式(1)~式(13),可得整閥框圖如圖3所示。i為伺服閥輸入指令信號;ui為調制后的指令信號;u0為輸入電機的電壓信號。

        該閥的動態(tài)響應測試結果表明:現(xiàn)有結構和電控方案下,控制壓力超調量大,響應調整時間長,即該閥穩(wěn)定性較差。由整閥框圖可知,系統(tǒng)響應與結構參數(shù)及PI控制方法均有關系。可通過理論分析獲得結構參數(shù)及控制方法的影響規(guī)律。

        圖3 BRDDPSV在任意工作點線性化框圖Fig.3 Block diagram linearized at any operating point of BRDDPSV

        2 結構參數(shù)的影響及優(yōu)化

        為分析結構參數(shù)對整閥動態(tài)特性的影響,去除PI控制器。由框圖3可知,整閥的傳遞函數(shù)如式(14)所示:它表示無控制器狀態(tài)時,輸出控制壓力和輸入電流之間的關系

        (14)

        式中:ωn為該閥的無阻尼自然頻率;ξ為該閥的阻尼比;k為該閥的閉環(huán)增益。

        2.1 結構參數(shù)的影響及優(yōu)化

        根據(jù)式(14)可以看出,電機轉動慣量、滑閥閥芯質量的增加將導致整閥的慣量增大,會降低系統(tǒng)響應帶寬,使整閥響應變慢;而電機力矩系數(shù)的增加,提高了整閥的自然頻率,降低了阻尼比,整閥動態(tài)響應變快。

        由式(10)可計算得出不同閥芯直徑時的穩(wěn)態(tài)液動力剛度曲線,如圖4所示。由計算結果可知,滑閥在運動至全行程的一半附近時,液動力會出現(xiàn)負剛度,不利于閥芯運動的穩(wěn)定。而減小滑閥直徑可以有效減小系統(tǒng)負向的剛度值。因此,選擇小直徑的閥芯利于削弱液動力對系統(tǒng)穩(wěn)定性的影響。

        圖4 閥芯直徑對穩(wěn)態(tài)液動力剛度的影響Fig.4 Effect of spool diameter on steady flow force stiffness

        圖5 滑閥預遮蓋量對壓力增益的影響Fig.5 Effect of sliding valve pre-cover on pressure gain

        壓力伺服閥采用正開口滑閥結構,負遮蓋量的大小直接影響滑閥的壓力增益kp,如圖5所示。從圖5中可以看出,壓力增益kp隨閥口負遮蓋量的減小而增大。根據(jù)式(14)可知,壓力增益的增加能夠加大整閥的開環(huán)增益,提高控制壓力響應的速度,利于改善整閥的響應性能。因此,在條件允許下應盡量減小正開口滑閥的負遮蓋量。

        2.2 驅動接口的影響及優(yōu)化

        由于安裝或加工誤差,會造成驅動接口處存在微小的間隙或過盈配合,將造成電機轉動過程中的空行程或過大的摩擦力。驅動接口配合偏差如圖6所示,間隙量記為Δx。

        1) 間隙量Δx>0

        Δx>0時,會造成電機轉動的一段空行程。則電機空轉的最大角度為

        (15)

        若Δx=0.05 mm,則Δα=2.05°(e0=1.4 mm)電機轉角產(chǎn)生了2.05°的死區(qū)。

        由于電機轉角的死區(qū),造成了旋轉直驅閥動態(tài)特性的改變,如圖7所示,為Δx>0時的整閥動態(tài)特性。

        圖6 驅動接口配合偏差Fig.6 Deviation of driven interface

        圖7 Δx>0時的整閥動態(tài)特性Fig.7 Valve dynamics at Δx> 0

        可見,驅動接口處的間隙會造成旋轉直驅閥的上升過程死區(qū)增大,同時其額定壓力從8 MPa降低為4.57 MPa。這是由于接口處的間隙造成工作過程中實際推動閥芯的電機轉角減小,導致滑閥驅動力下降,從而造成滑閥位移減小,旋轉直驅閥的控制壓力降低。此外,在驅動小球越過間隙后,會以一定的速度與柱形孔相碰,將使整閥穩(wěn)定性下降,如圖7中間隙為0.05 mm的階躍響應曲線所示。

        2) 間隙量Δx<0

        Δx<0表示驅動接口處存在較大的應力和摩擦力。這里假設摩擦力為常數(shù),設摩擦力產(chǎn)生的阻力矩Tf′=0.01 N·m。如圖8所示為摩擦力矩對旋轉直驅閥動態(tài)特性的影響。與正間隙不同,當存在摩擦力時,閥在階躍響應時不存在死區(qū)。這是由于在這種情況下,驅動小球仍然與柱形孔壁接觸,電機不存在空轉。僅僅由于阻力矩增加,驅動力相對不足,穩(wěn)態(tài)值發(fā)生了變化。

        圖8 摩擦力矩對整閥動態(tài)特性的影響Fig.8 Effect of frictional torque on dynamic characteristics of valve

        由于穩(wěn)態(tài)值的靜差可通過積分作用消除,通過分析可知,驅動接口處的間隙容易引發(fā)閥芯運動的不穩(wěn)定,因此,驅動接口處應選擇小過盈配合,即Δx<0。

        2.3 控制策略的優(yōu)化

        旋轉直驅電液壓力伺服閥通過電機旋轉帶動滑閥平動,從而實現(xiàn)壓力控制。控制策略和控制參數(shù)的選擇是否合理將直接決定旋轉直驅伺服閥動態(tài)特性是否符合使用預期要求。

        2.3.1 積分分離處理

        PID控制器中的積分環(huán)節(jié)的作用是消除穩(wěn)態(tài)值的靜差,保證系統(tǒng)輸出的精度。但對于飛機剎車系統(tǒng)中壓力控制伺服而言,在啟動、停止或大幅度增減設定值時,系統(tǒng)在短時間內(nèi)輸出有很大的偏差,積分環(huán)節(jié)的累積作用會使系統(tǒng)控制壓力出現(xiàn)大幅變化,即表現(xiàn)出超調量大的特征。因此,考慮采用積分分離式PI控制,如圖9所示為旋轉直驅伺服閥的控制框圖。

        積分分離可認為是在積分運算前添加開關環(huán)節(jié),在指定閾值條件時啟動或關閉積分運算[21]。其主要作用是當輸入值變化量較大時,即系統(tǒng)輸入的指令值與系統(tǒng)實際輸出值的偏差大于該積分分離閾值時放棄積分,反之,進行積分消除靜差。設e(k)為偏差值,ε為閾值(ε>0),則積分分離可用分段函數(shù)表示,β為積分項前的系數(shù):

        (16)

        圖9 電反饋控制框圖Fig.9 Block diagram of electric feedback control

        圖10 帶與不帶積分分離的階躍響應曲線(ε=2 mV)Fig.10 Curve of step response with and without integral separation (ε=2 mV)

        將上述積分分離式PI控制器加入仿真模型中,可得圖10所示的階躍輸入情況下RDDPV伺服閥動態(tài)響應。從圖10中可以看出,無積分作用時,穩(wěn)態(tài)值偏離設定值,靜差無法消除;持續(xù)積分作用時,雖然穩(wěn)態(tài)值無靜差,響應速度較快,但會有較大的壓力超調。采用積分分離后,合理的積分閾值可以使控制器在壓力超調處停止積分,避免控制壓力的大幅持續(xù)增加,同時達到消除穩(wěn)態(tài)值的靜差的目的。

        2.3.2 設置動壓反饋

        為了縮短動態(tài)響應的調整時間,考慮增設動壓反饋校正來進一步增大RDDPV伺服閥的阻尼,以提高系統(tǒng)穩(wěn)定性。如圖11所示的電控框圖,在控制壓力的反饋環(huán)節(jié)并聯(lián)加入了控制壓力的微分項(du/dt)。

        動壓反饋的輸入為壓力傳感器的采集信號UpL,反饋至輸入端之后,其傳遞函數(shù)為

        G=kf3s

        (17)

        因此,壓力傳感器的輸出電壓修正為

        UpL=kf2(1+kf3s)pL

        (18)

        式中:kf3為動壓反饋系數(shù)。

        將上述動壓反饋模型加入所建仿真模型中,可得圖12所示的階躍輸入情況下RDDPV伺服閥的動態(tài)響應。可見,動壓反饋將控制壓力的變化率進行反饋,避免了控制壓力的大幅波動,有效縮短了系統(tǒng)響應的調整時間,提高了系統(tǒng)的穩(wěn)定性。

        圖11 采用動壓反饋校正的電控框圖Fig.11 Block diagram of electric feedback control with dynamic pressure correction

        圖12 帶與不帶動壓反饋的階躍響應曲線(kf3=0.000 5 mV·MPa/s)Fig.12 Curve of step response with and without dynamic pressure correction (kf3=0.000 5 mV·MPa/s)

        3 優(yōu)化驗證實驗

        3.1 半物理實驗臺及實驗環(huán)境

        為免去調試電控參數(shù)時頻繁搭建模擬電路的繁瑣過程,本實驗采用半物理仿真平臺進行。所謂的半物理即實驗臺架采用硬件結構,而電控部分采用軟件編寫。旋轉直驅壓力伺服閥測試臺架如圖13所示。圖13中右側的性能測試設備包括測試PC機和信號發(fā)生器,內(nèi)置伺服閥動態(tài)特性測試軟件;左側的數(shù)字控制PC為便攜式筆記本;控制PC機通過信號采集、輸出及放大設備與伺服閥和信號發(fā)生器相連。測試時輸入3.2 mA的階躍電流信號。采用15#液壓油,介質溫度為25 ℃±1 ℃。

        圖13 BRDDPSV測試臺架Fig.13 Test bench of BRDDPSV

        3.2 動態(tài)特性測試結果

        根據(jù)結構和控制策略的分析,優(yōu)化后的參數(shù)選取如表1所示。電控參數(shù)在半物理仿真平臺上進行調試獲得。

        根據(jù)表1中優(yōu)化后的結構和電控參數(shù)進行階躍響應的實驗曲線如圖14所示,階躍響應動態(tài)指標如表2所示??梢?,該閥為欠阻尼系統(tǒng),優(yōu)化前,整閥開啟時的響應超調量大(超調量約為25%穩(wěn)態(tài)值)、調整時間較長(約75 ms),穩(wěn)定性較差。優(yōu)化后,控制壓力階躍響應的超調量明顯降低,約為6.25%穩(wěn)態(tài)值,調整時間縮至30 ms左右。雖然控制壓力階躍響應的上升時間有所增加,但可以控制在20 ms以內(nèi),基本可以滿足飛機剎車壓力控制閥的使用需求。

        表1 BRDDPSV的優(yōu)化參數(shù)Table 1 Optimization parameters of BRDDPSV

        圖14 優(yōu)化前后BRDDPSV的階躍響應(實驗)Fig.14 Step response of BRDDPSV before and after optimization (experiment)

        性能指標優(yōu)化前優(yōu)化后上升時間/ms<25820調整時間/ms<507530控制壓力超調量/%<10256.25

        4 結 論

        1) 提高整閥響應穩(wěn)定性和可靠性的措施:消除整閥負剛度,包括減小閥口負遮蓋量以提高壓力增益,減小閥芯直徑以減小最大負剛度值;提高電機抗干擾能力,如提高電機電流-力矩系數(shù),及在驅動接口處小球和柱形孔間采用小過盈配合。

        2) 優(yōu)化控制方法:采用積分分離式的控制器。理論和實驗表明,適當?shù)姆e分分離閾值可以有效降低控制壓力超調值,使整閥響應趨于穩(wěn)定;設置動壓反饋校正。將控制壓力的變化率反饋到輸入端進行控制,這相當于在系統(tǒng)中增加了阻尼,整閥的階躍響應調整時間明顯縮短。整閥的上升時間稍有加大,但可以控制到20 ms以內(nèi)。

        猜你喜歡
        滑閥階躍轉角
        玩轉角的平分線
        基于階躍雙包層光纖的螺旋型光纖傳感器
        汽輪機啟動滑閥改造經(jīng)驗反饋報告
        中文信息(2021年5期)2021-03-27 12:00:19
        FD1(G)型防滑閥A3修存在問題及對策
        探討單位階躍信號的教學
        三次“轉角”遇到愛
        解放軍健康(2017年5期)2017-08-01 06:27:42
        永春堂贏在轉角
        下一個轉角:邁出去 開啟“智”造時代
        均壓結構在液壓滑閥設計中的應用
        大電機技術(2015年5期)2015-02-27 08:22:32
        一種階躍函數(shù)在矩形時間窗口頻域特性的分析方法
        国产在线观看视频一区二区三区| 亚洲精品美女久久久久久久| 日韩精品一区二区亚洲av| 全免费a级毛片免费看| 国产激情一区二区三区在线蜜臀| 国产一区二区三区护士| 亚洲av午夜精品无码专区| 婷婷久久久亚洲欧洲日产国码av| 亚洲成av人最新无码| 8av国产精品爽爽ⅴa在线观看| 日韩精品成人一区二区在线观看| 国产毛片一区二区三区| 伊人青青草综合在线视频免费播放| 亚洲成av人在线观看网址| 成人精品综合免费视频| 久9热免费精品视频在线观看| 日韩有码中文字幕在线视频| 免费网站内射红桃视频| 成人免费视频在线观看| 91色婷婷成人精品亚洲| 日本一二三区在线不卡| 国产一区二区三区久久悠悠色av| 女人被狂躁高潮啊的视频在线看| 亚洲成a人片在线网站| 一本大道加勒比东京热| 狠狠cao日日橹夜夜十橹| 人人爽久久涩噜噜噜丁香| 亚洲午夜精品久久久久久一区| 农村国产毛片一区二区三区女| 国产91精品在线观看| 国产成人午夜福利在线观看| 最好看的最新高清中文视频| 蜜芽尤物原创AV在线播放| 日本久久精品福利视频| 特级毛片a级毛片100免费播放| 亚洲一区二区三区国产精华液| 亚洲精品一区二区三区播放| 亚洲av一区二区三区色多多| 日本无码人妻波多野结衣| 国产高清国内精品福利99久久| 亚洲第一页在线观看视频网站|