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

基于自適應(yīng)模糊PID的連續(xù)管作業(yè)機(jī)注入頭速度控制研究

2024-03-14 02:15:04湯清源丁悅曾行健康凱武一葦杜宇成白龍
機(jī)床與液壓 2024年4期
關(guān)鍵詞:模型系統(tǒng)

湯清源,丁悅,曾行健,康凱,武一葦,杜宇成,白龍

(1.中石油江漢機(jī)械研究所有限公司,湖北荊州 434000;2.華中科技大學(xué)機(jī)械科學(xué)與工程學(xué)院,湖北武漢 430074)

0 前言

近年來(lái),隨著國(guó)內(nèi)先進(jìn)技術(shù)的發(fā)展,連續(xù)管作業(yè)機(jī)已廣泛用于鉆井、壓裂[1-3],是現(xiàn)代油井開發(fā)的重要裝備。連續(xù)管作業(yè)機(jī)的主要部件有注入頭、滾筒、液壓動(dòng)力單元控制室等[4],工作過(guò)程中核心的關(guān)鍵控制技術(shù)是注入頭起/下管速度精準(zhǔn)控制。目前連續(xù)管作業(yè)機(jī)液壓控制系統(tǒng)中的注入頭速度控制技術(shù)仍存在著一些技術(shù)問(wèn)題,需要根據(jù)經(jīng)驗(yàn)實(shí)時(shí)手動(dòng)調(diào)節(jié)泵的排量從而控制馬達(dá)轉(zhuǎn)速輸出,不僅對(duì)操作人員要求較高,而且注入頭液壓系統(tǒng)很難在較短時(shí)間內(nèi)達(dá)到期望速度,增加了作業(yè)的難度[5]。為此,文中將連續(xù)管作業(yè)機(jī)注入頭液壓模型簡(jiǎn)化為閉式泵控馬達(dá)子系統(tǒng)并對(duì)其進(jìn)行控制。

近年來(lái),泵控馬達(dá)液壓傳動(dòng)子系統(tǒng)引起了廣泛的關(guān)注[6-7],目前對(duì)泵-馬達(dá)液壓傳動(dòng)子系統(tǒng)的研究主要集中在效率特性和速度控制策略上[8]。文中所研究泵控馬達(dá)子系統(tǒng)是變量泵-定量馬達(dá)系統(tǒng)且采用閉式回路,是一種典型的容積調(diào)速系統(tǒng),此系統(tǒng)中液壓泵及液壓馬達(dá)通過(guò)改變流量進(jìn)行動(dòng)力傳遞[9],其主要優(yōu)點(diǎn)是沒(méi)有節(jié)流損失和溢流損失、回路效率高,適用于要求響應(yīng)速度快的高速調(diào)速系統(tǒng)[10]。當(dāng)容積調(diào)速系統(tǒng)負(fù)載的轉(zhuǎn)速需要發(fā)生變化時(shí),通過(guò)調(diào)節(jié)變量泵的排量控制定量馬達(dá)的轉(zhuǎn)速,從而達(dá)到期望的負(fù)載轉(zhuǎn)速。然而,泵控馬達(dá)子系統(tǒng)在變速輸入和負(fù)載擾動(dòng)時(shí)的可控性、快速性和穩(wěn)定性較差[11]。因此研究合適且簡(jiǎn)便的方法對(duì)實(shí)現(xiàn)馬達(dá)轉(zhuǎn)速控制的快速性、準(zhǔn)確性和穩(wěn)定性具有現(xiàn)實(shí)且重要的工程意義[12]。

自適應(yīng)模糊PID控制也廣泛用于液壓回路控制領(lǐng)域中[13-17]。故本文作者以閉式泵控馬達(dá)子系統(tǒng)為研究對(duì)象,以穩(wěn)定快速精確控制馬達(dá)轉(zhuǎn)速為目標(biāo),通過(guò)分析馬達(dá)運(yùn)動(dòng)的工作原理,不再建立液壓系統(tǒng)復(fù)雜數(shù)學(xué)模型,而是以液壓系統(tǒng)的工作原理為基礎(chǔ),通過(guò)AMESim軟件搭建泵控馬達(dá)子系統(tǒng)的液壓模型,并利用MATLAB/Simulink軟件構(gòu)建PID及自適應(yīng)模糊PID 2種控制系統(tǒng)的仿真模型;通過(guò)AMESim-Simulink的軟件接口進(jìn)行聯(lián)合仿真分析,實(shí)現(xiàn)對(duì)泵控馬達(dá)子系統(tǒng)馬達(dá)轉(zhuǎn)速的控制,克服變速輸入時(shí)泵-馬達(dá)液壓傳動(dòng)子系統(tǒng)輸出速度可控性和穩(wěn)定性差的問(wèn)題。

1 泵控馬達(dá)系統(tǒng)工作原理

將連續(xù)管作業(yè)機(jī)注入頭液壓原理簡(jiǎn)化為圖1所示的液壓模型,整個(gè)回路分為執(zhí)行油路部分、負(fù)載油路部分及時(shí)序信號(hào)部分。執(zhí)行油路部分采用典型的閉式泵控馬達(dá)調(diào)速系統(tǒng),它是由電動(dòng)機(jī)1、變量馬達(dá)2、電液比例方向閥6以及定量馬達(dá)7組合而成的一種液壓系統(tǒng),模擬連續(xù)管作業(yè)機(jī)注入頭的控制油路部分,它主要負(fù)責(zé)保證連續(xù)管作業(yè)機(jī)工作過(guò)程中起/下管的速度,并時(shí)刻監(jiān)測(cè)作業(yè)機(jī)的工作狀態(tài),在出現(xiàn)異常狀態(tài)時(shí)根據(jù)需要切換到降速起/下模式、緊急制動(dòng)或蠕動(dòng)模式等。

圖1 泵控馬達(dá)系統(tǒng)液壓回路原理

負(fù)載油路部分由馬達(dá)帶動(dòng)負(fù)載回路里的油泵,油泵出油口連接溢流閥,溢流閥處設(shè)定壓力模擬負(fù)載源,負(fù)載回路是連續(xù)管作業(yè)機(jī)的執(zhí)行對(duì)象,它真實(shí)地模擬作業(yè)機(jī)在不同工況下的負(fù)載情況,當(dāng)負(fù)載發(fā)生變化時(shí),會(huì)及時(shí)反饋,執(zhí)行回路進(jìn)行相應(yīng)的控制調(diào)整,以達(dá)到連續(xù)管作業(yè)機(jī)平穩(wěn)工作的目的。執(zhí)行油路部分控制可以通過(guò)調(diào)節(jié)電動(dòng)機(jī)1的轉(zhuǎn)速,從而改變變量泵的排量,達(dá)到控制定量馬達(dá)轉(zhuǎn)速的目的;同時(shí)也可以通過(guò)調(diào)節(jié)電液比例方向閥6的開度以及換向,對(duì)定量馬達(dá)的轉(zhuǎn)速進(jìn)行控制。而時(shí)序信號(hào)部分主要是控制負(fù)載液壓油流向,根據(jù)液壓馬達(dá)及負(fù)載油泵的正反轉(zhuǎn)情況調(diào)整負(fù)載油的流向。

對(duì)電動(dòng)機(jī)1以及電液比例方向閥6的控制是控制定量馬達(dá)轉(zhuǎn)速的關(guān)鍵環(huán)節(jié)。傳統(tǒng)的開環(huán)控制方法通過(guò)給定電液比例方向閥6恒定開度以及合適的電動(dòng)機(jī)1轉(zhuǎn)速,從而得到期望的馬達(dá)轉(zhuǎn)速輸出;而文中所研究的閉環(huán)控制方法則根據(jù)定量馬達(dá)轉(zhuǎn)速傳感器的反饋值與給定信號(hào)的偏差,按照預(yù)設(shè)的算法規(guī)則進(jìn)行計(jì)算控制,得到電動(dòng)機(jī)的轉(zhuǎn)速以及電液比例方向閥的開度,其原理框圖如圖2所示。

圖2 泵控馬達(dá)調(diào)速系統(tǒng)原理框圖

2 聯(lián)合仿真分析

2.1 液壓模型結(jié)果分析

根據(jù)圖1所示的泵控馬達(dá)系統(tǒng)的液壓回路在AMESim軟件中進(jìn)行仿真分析,設(shè)置期望轉(zhuǎn)速如表1所示,得到圖3所示的結(jié)果。其中圖1中的參數(shù)設(shè)定及輸入信號(hào)列于表2,其余參數(shù)為默認(rèn)值。

表1 期望轉(zhuǎn)速

表2 泵控馬達(dá)系統(tǒng)輸入信號(hào)

圖3 液壓模型定量馬達(dá)轉(zhuǎn)速輸出曲線

由圖3可以看出:無(wú)反饋調(diào)節(jié)僅通過(guò)液壓模型內(nèi)部操作,液壓馬達(dá)的轉(zhuǎn)速存在超調(diào)較大、響應(yīng)速度較慢、穩(wěn)態(tài)誤差較大的問(wèn)題,且無(wú)法準(zhǔn)確穩(wěn)定在期望馬達(dá)轉(zhuǎn)速,動(dòng)態(tài)特性較差,故需要在此泵控馬達(dá)系統(tǒng)液壓模型中加入反饋調(diào)節(jié)進(jìn)行精準(zhǔn)調(diào)控。

2.2 常規(guī)PID控制及結(jié)果分析

常規(guī)PID控制器是一個(gè)經(jīng)典的負(fù)反饋控制回路,它由比例調(diào)節(jié)P、積分調(diào)節(jié)I、微分調(diào)節(jié)D三部分構(gòu)成,因其結(jié)構(gòu)及原理較為簡(jiǎn)單,應(yīng)用已十分普遍。常規(guī)PID控制算法中輸入r(t)與輸出o(t)之間的關(guān)系式為

(1)

式中:Kp為比例系數(shù);Ki(=Kp/Ti)為積分系數(shù);Ti為積分時(shí)間常數(shù);Kd(=KpTd)為微分系數(shù);Td為微分時(shí)間常數(shù)。

將常規(guī)PID控制算法運(yùn)用于泵控馬達(dá)系統(tǒng)的工作原理框圖如圖4所示,使用過(guò)程中僅需要調(diào)節(jié)3個(gè)參數(shù)(Kp、Ki、Kd),通過(guò)對(duì)比期望馬達(dá)轉(zhuǎn)速與輸出馬達(dá)的轉(zhuǎn)速得到偏差r(t)。再通過(guò)PID控制算法得到輸出o(t),將它輸入到液壓系統(tǒng)中進(jìn)行液壓模型內(nèi)部調(diào)節(jié),從而得到期望的馬達(dá)轉(zhuǎn)速輸出。

圖4 常規(guī)PID控制工作原理

聯(lián)合仿真平臺(tái)的搭建需要通過(guò)AMESim與Simulink軟件之間的接口實(shí)現(xiàn),通過(guò)接口可以將AMESim中所建立的液壓模型轉(zhuǎn)換成Simulink中控制模型所需的S函數(shù)。聯(lián)合仿真有2種接口設(shè)置:標(biāo)準(zhǔn)接口(Simulink)、聯(lián)合仿真接口(SimuCosim)。標(biāo)準(zhǔn)接口采用連續(xù)的仿真過(guò)程,由Simulink所選定的求解器計(jì)算,且在AMESim與Simulink之間實(shí)時(shí)傳輸狀態(tài)變量以及輸入/輸出變量。而聯(lián)合仿真接口無(wú)需中間變量,它是一個(gè)不連續(xù)的仿真過(guò)程,2個(gè)軟件采用各自的求解器計(jì)算模型結(jié)果,只交換涉及到的輸入輸出的變量,求解速度更快,故文中采用聯(lián)合仿真接口進(jìn)行聯(lián)合仿真平臺(tái)的搭建。圖5所示為聯(lián)合仿真接口的數(shù)據(jù)交換示意。

圖5 聯(lián)合仿真接口的數(shù)據(jù)交換示意

根據(jù)常規(guī)PID控制泵控馬達(dá)系統(tǒng)原理及所搭建的液壓模型,通過(guò)上述接口聯(lián)合AMESim與Simulink軟件,完成聯(lián)合仿真設(shè)計(jì),建立如圖6—7所示的AMESim-Simulink泵控馬達(dá)系統(tǒng)速度PID控制仿真模型,可以獲得圖8所示的定量馬達(dá)轉(zhuǎn)速輸出曲線。

圖6 泵控馬達(dá)系統(tǒng)常規(guī)PID控制聯(lián)合仿真模型

圖7 泵控馬達(dá)系統(tǒng)常規(guī)PID控制Simulink模型

圖8 常規(guī)PID控制定量馬達(dá)轉(zhuǎn)速n輸出曲線

利用階躍信號(hào)跟蹤仿真結(jié)果分析此液壓模型的調(diào)速效果,設(shè)定馬達(dá)期望轉(zhuǎn)速為nm=200 r/min,控制電動(dòng)機(jī)處PID的參數(shù)分別設(shè)置為:Kp=2、Ki=100、Kd=0.16,控制電磁閥處的PID參數(shù)分別設(shè)置為:Kp=1、Ki=1、Kd=0。由圖8可知:在所受負(fù)載恒定情況下,采用常規(guī)PID控制算法的速度調(diào)節(jié)方式,定量馬達(dá)轉(zhuǎn)速在0.4 s左右到達(dá)200 r/min并保持穩(wěn)定,在速度調(diào)節(jié)過(guò)程中未出現(xiàn)嚴(yán)重超調(diào)現(xiàn)象,整體響應(yīng)較為迅速,并且達(dá)到穩(wěn)態(tài)的時(shí)間較短。

3 自適應(yīng)模糊PID控制器

3.1 自適應(yīng)模糊PID控制理論

模糊控制是基于規(guī)則的專家系統(tǒng)、模糊集理論和控制理論相結(jié)合而產(chǎn)生的成果,由模糊控制器和控制對(duì)象組成,其結(jié)構(gòu)如圖9所示。

圖9 模糊控制系統(tǒng)組成

如圖10所示,文中所設(shè)計(jì)的自適應(yīng)模糊PID控制器是一個(gè)二輸入三輸出的結(jié)構(gòu),以系統(tǒng)輸入與輸出間偏差e及偏差變化率ec為輸入,模糊系統(tǒng)通過(guò)模糊化、模糊推理和解模糊得到PID的調(diào)整量ΔKp、ΔKi和ΔKd,輸出PID 3個(gè)控制參數(shù)的變化量,即在系統(tǒng)工作中實(shí)時(shí)采集輸入與輸出偏差及偏差變化率,由模糊規(guī)則對(duì)控制參數(shù)進(jìn)行實(shí)時(shí)修改,使它滿足當(dāng)下偏差及偏差變化率的要求。自適應(yīng)模糊PID控制器的公式為

(2)

圖10 自適應(yīng)模糊PID控制原理

式中:ΔK′up、ΔK′ui、ΔK′ud為調(diào)整系數(shù);ΔKp、ΔKi、ΔKd為每個(gè)時(shí)間周期后PID參數(shù)的變化量。

文中選用二輸入三輸出的自適應(yīng)模糊PID控制器來(lái)調(diào)節(jié)定量馬達(dá)的轉(zhuǎn)速,由于三角隸屬度函數(shù)計(jì)算簡(jiǎn)單、敏銳度高,可以滿足泵控馬達(dá)系統(tǒng)的精度要求,故輸入輸出均選用三角隸屬度函數(shù)。輸入、輸出量模糊化處理在MATLAB軟件中的Fuzzy Logic Designer工具中設(shè)置。

設(shè)輸入e、ec,輸出ΔKp、ΔKi、ΔKd的模糊子集為{NB,NM,NS,Z,PS,PM,PB},制定模糊控制規(guī)則如表3—5所示,其對(duì)ΔKp、ΔKi、ΔKd的參數(shù)調(diào)整規(guī)則如下:

表3 ΔKp模糊控制規(guī)則

表4 ΔKi模糊控制規(guī)則

表5 ΔKd模糊控制規(guī)則

(1)在系統(tǒng)剛啟動(dòng)時(shí),系統(tǒng)偏差值e最大,為提高系統(tǒng)的響應(yīng)速度應(yīng)適當(dāng)增大比例系數(shù)ΔKp,同時(shí)為防止出現(xiàn)超調(diào)量較大以及初始誤差變化較大造成的積分/微分飽和,ΔKi取在零值附近,ΔKd取中等值大小;

(2)當(dāng)系統(tǒng)運(yùn)轉(zhuǎn)一段時(shí)間后,偏差e在中等范圍內(nèi)變化,選用較小ΔKp以防止出現(xiàn)較大超調(diào)量,同時(shí)ΔKi、ΔKd均在中間值附近變化;

(3)當(dāng)系統(tǒng)趨于穩(wěn)定時(shí),偏差e數(shù)值較小,應(yīng)選擇較大ΔKp及ΔKi以增強(qiáng)整個(gè)系統(tǒng)的穩(wěn)定性,并且應(yīng)當(dāng)在ec較小時(shí)選擇較大ΔKd以減小系統(tǒng)振蕩。

根據(jù)上述規(guī)則完成設(shè)置后,可以獲得馬達(dá)轉(zhuǎn)速偏差e與偏差變化率ec以及控制參數(shù)ΔKp、ΔKi、ΔKd的變化,所設(shè)置的規(guī)則輸入量與輸出量之間的關(guān)系曲面如圖11所示。

圖11 輸入量與輸出量之間的關(guān)系曲面

自適應(yīng)模糊PID控制器的最后一步需要將上述模型推理所產(chǎn)生的模糊值轉(zhuǎn)化為清晰的控制量,故進(jìn)行去模糊化處理。常見的去模糊方法有:重心法、加權(quán)平均法和最大隸屬度法,文中選用重心法作為解模糊方法,取隸屬度函數(shù)曲線與橫坐標(biāo)圍成面積的重心作為模糊推理的最終輸出值,所采用的公式為

(3)

式中:z0為模糊控制器輸出量解模糊后的精確值;zi為模糊控制器論域內(nèi)的值;μc(zi)為zi的隸屬度值。

3.2 自適應(yīng)模糊PID控制結(jié)果分析

按照?qǐng)D10所示的自適應(yīng)模糊PID控制原理框圖,在AMESim中構(gòu)建泵控馬達(dá)系統(tǒng)自適應(yīng)模糊PID控制液壓模型,通過(guò)采用自適應(yīng)模糊PID控制電動(dòng)機(jī)的輸入轉(zhuǎn)速,如圖12所示;同時(shí)在MATLAB/Simulink中搭建自適應(yīng)模糊PID控制Simulink控制模型,如圖13所示。

圖12 泵控馬達(dá)系統(tǒng)自適應(yīng)模糊PID控制液壓模型

圖13 自適應(yīng)模糊PID控制Simulink控制模型

與常規(guī)PID控制模型相同,利用階躍信號(hào)跟蹤仿真結(jié)果分析此液壓模型的速度調(diào)節(jié)效果。同樣設(shè)定馬達(dá)期望轉(zhuǎn)速為nm=200 r/min,模糊PID的參數(shù)設(shè)置如圖13所示。分析圖14可知:在所受負(fù)載恒定情況下,采用自適應(yīng)模糊PID控制算法的速度調(diào)節(jié)方式,定量馬達(dá)轉(zhuǎn)速在0.3 s左右達(dá)到200 r/min并保持穩(wěn)定,在速度調(diào)節(jié)過(guò)程中未出現(xiàn)嚴(yán)重超調(diào)現(xiàn)象,相較于常規(guī)PID控制整體響應(yīng)更為迅速,并且達(dá)到穩(wěn)態(tài)的時(shí)間更短。

圖14 常規(guī)PID與自適應(yīng)模糊PID控制結(jié)果對(duì)比

3.3 變轉(zhuǎn)速馬達(dá)轉(zhuǎn)速控制

在連續(xù)管作業(yè)機(jī)工作過(guò)程中,起/下管轉(zhuǎn)速根據(jù)工況會(huì)隨時(shí)增加或減少,于是要求注入頭處馬達(dá)轉(zhuǎn)速能隨時(shí)跟上突然增加或減少的期望轉(zhuǎn)速。在某固定工況下,利用聯(lián)合仿真模型跟蹤階躍信號(hào)分析其調(diào)速效果,設(shè)定期望轉(zhuǎn)速大小如表1所示,運(yùn)行AMESim-Simulink泵控馬達(dá)系統(tǒng),獲得圖15所示的定量馬達(dá)轉(zhuǎn)速輸出曲線。可知:在固定工況下,馬達(dá)期望轉(zhuǎn)速如果突然增加,PID控制及自適應(yīng)模糊PID控制下的馬達(dá)轉(zhuǎn)速相比普通液壓調(diào)節(jié)穩(wěn)定更快、未出現(xiàn)較大超調(diào),但PID控制在馬達(dá)轉(zhuǎn)速減小時(shí)會(huì)出現(xiàn)較大的振蕩和超調(diào);自適應(yīng)模糊PID控制在馬達(dá)轉(zhuǎn)速陡然減小時(shí),雖然也存在超調(diào)現(xiàn)象,但其達(dá)到穩(wěn)定所需時(shí)間較短。

圖15 常規(guī)PID控制定量馬達(dá)變轉(zhuǎn)速輸出曲線

3.4 變負(fù)載馬達(dá)轉(zhuǎn)速控制

在連續(xù)管作業(yè)機(jī)工作過(guò)程中,馬達(dá)受到的負(fù)載可能會(huì)發(fā)生階段性的突變,因此,可以用方波負(fù)載仿真工況的變化情況。給定負(fù)載信號(hào)為方波信號(hào),最小值為5 N·m,最大值是100 N·m。其中馬達(dá)的期望轉(zhuǎn)速為nm=100 r/min,結(jié)果如圖16所示。

圖16 負(fù)載變化下馬達(dá)轉(zhuǎn)速輸出曲線

由圖16可知:在工況發(fā)生變化的情況下,即加入方波負(fù)載干擾,PID控制相比自適應(yīng)模糊PID控制,存在超調(diào)量較大且再次達(dá)到穩(wěn)態(tài)的時(shí)間較長(zhǎng)等缺點(diǎn)。由此可見,自適應(yīng)模糊PID控制算法對(duì)于負(fù)載變化時(shí)連續(xù)管作業(yè)機(jī)注入頭馬達(dá)轉(zhuǎn)速的控制效果最佳。

4 結(jié)語(yǔ)

文中對(duì)連續(xù)管作業(yè)機(jī)中注入頭起/下管過(guò)程的速度液壓系統(tǒng)進(jìn)行了建模仿真分析,根據(jù)所設(shè)計(jì)的簡(jiǎn)化泵控馬達(dá)液壓子系統(tǒng)原理,提出了馬達(dá)轉(zhuǎn)速PID的控制方法,建立了聯(lián)合仿真模型;利用AMESim軟件建立泵控馬達(dá)子系統(tǒng)液壓仿真模型,同時(shí)利用MATLAB/Simulink軟件設(shè)計(jì)出自適應(yīng)模糊PID控制器,建立了泵控馬達(dá)速度控制聯(lián)合仿真模型,結(jié)合PID控制與模糊控制規(guī)則進(jìn)行速度控制,與常規(guī)液壓控制模型進(jìn)行對(duì)比,結(jié)果表明自適應(yīng)模糊PID控制取得了較好的速度控制效果。

文中通過(guò)構(gòu)建泵控馬達(dá)子系統(tǒng)液壓模型進(jìn)行聯(lián)合仿真,避免了復(fù)雜數(shù)學(xué)模型的構(gòu)建,降低了模型搭建難度并保證了仿真分析的準(zhǔn)確性。在輸入階躍信號(hào)的情況下,自適應(yīng)模糊PID控制將模糊規(guī)則與常規(guī)PID控制進(jìn)行結(jié)合,模型的響應(yīng)速度更快、無(wú)超調(diào)和滯后現(xiàn)象、穩(wěn)態(tài)誤差更小,具有良好的動(dòng)態(tài)特性。

猜你喜歡
模型系統(tǒng)
一半模型
Smartflower POP 一體式光伏系統(tǒng)
WJ-700無(wú)人機(jī)系統(tǒng)
ZC系列無(wú)人機(jī)遙感系統(tǒng)
重要模型『一線三等角』
重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
基于PowerPC+FPGA顯示系統(tǒng)
半沸制皂系統(tǒng)(下)
連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
3D打印中的模型分割與打包
主站蜘蛛池模板: 国产成年无码AⅤ片在线| jizz国产视频| 免费人成视网站在线不卡| 欧美日本激情| 国产区免费精品视频| 欧美人人干| 国产后式a一视频| 国产中文在线亚洲精品官网| 国产亚洲美日韩AV中文字幕无码成人 | 国产精品极品美女自在线看免费一区二区| 中文字幕人妻无码系列第三区| 国产玖玖视频| 一区二区三区在线不卡免费| 91麻豆精品国产91久久久久| 国产视频资源在线观看| 99久久99视频| 国产一区自拍视频| 情侣午夜国产在线一区无码| 无码免费试看| 国产成人无码综合亚洲日韩不卡| 亚洲欧美人成人让影院| 日本不卡视频在线| 九九热视频在线免费观看| 欧美人与性动交a欧美精品| 亚洲AV无码乱码在线观看裸奔| 亚洲美女视频一区| 免费观看男人免费桶女人视频| 国产男女免费视频| 一级毛片在线播放| 欧美精品二区| 国产精品自在自线免费观看| а∨天堂一区中文字幕| 欧美狠狠干| 国产精品三级av及在线观看| 好吊色国产欧美日韩免费观看| 干中文字幕| 色播五月婷婷| 国产福利一区在线| 亚洲午夜福利精品无码不卡| 国产成人亚洲综合a∨婷婷| 亚洲欧美不卡| 激情六月丁香婷婷| 欧美在线国产| 人妻中文字幕无码久久一区| 99人体免费视频| 久久国产拍爱| 亚洲第一视频网| 不卡无码网| 波多野结衣一二三| 欧美中文一区| 国产精品久久久久久久久久久久| 综合五月天网| 99久久精品美女高潮喷水| 久久一本日韩精品中文字幕屁孩| 欧美亚洲日韩中文| 国产精品主播| 免费看美女自慰的网站| 国产精品极品美女自在线网站| 亚洲成人黄色网址| 欧美日本在线观看| 免费不卡视频| 在线观看欧美国产| 国产精品蜜臀| 五月天久久婷婷| 欧美成人看片一区二区三区| 精品国产免费观看| 99久久精品视香蕉蕉| 国产在线日本| 欧美在线三级| 尤物精品国产福利网站| 久久久久久久蜜桃| 麻豆国产精品一二三在线观看| 亚洲AV无码久久精品色欲| 伊人久久大香线蕉aⅴ色| 日韩在线视频网站| 中国黄色一级视频| 色综合热无码热国产| 色综合狠狠操| 欧美日韩精品综合在线一区| 国产免费人成视频网| 国产毛片片精品天天看视频| 久久综合五月婷婷|