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

    梯度校正法燃?xì)鉄嶂禍y量裝置抽氣系統(tǒng)建模

    2021-11-12 02:11:26王玉剛張洪軍趙曉東
    中國測試 2021年10期
    關(guān)鍵詞:抽氣階躍壓差

    劉 偉,王玉剛,張洪軍,趙曉東,潘 江

    (中國計量大學(xué)計量測試工程學(xué)院,浙江 杭州 310018)

    0 引 言

    在天然氣能量計量方法中,以燃燒法測量可燃?xì)怏w熱值是能源領(lǐng)域的基礎(chǔ)問題,也是目前研究熱點。國外很多學(xué)者和計量機構(gòu)已經(jīng)在這方面有較為深入的研究并取得長足進步。20世紀(jì)30年代美國國家標(biāo)準(zhǔn)局的Rossini建立了電能標(biāo)定的等溫式熱值測量計,以電校正的形式將熱值溯源到質(zhì)量、溫度、電流和時間等國際單位上[1-3]。20世紀(jì)80年代,中國計量科學(xué)研究院在Rossini方法的基礎(chǔ)上用氧彈法測量純甲烷的熱值,于2012年開始研制基準(zhǔn)氣體熱量計[4-5]。此后,中國石油西南油氣田分公司天然氣研究院、中國計量大學(xué)等研究機構(gòu)或高校,在天然氣熱值測量技術(shù)、配氣比例、熱值測量主體、熱量計非穩(wěn)態(tài)溫度場優(yōu)化分析、抽氣系統(tǒng)設(shè)計等方面取得較大進展[6-10]。中國計量大學(xué)研制的氣體熱值測量系統(tǒng)是在Rossini法的基礎(chǔ)上改進而來,目前該裝置已初步成型。但該系統(tǒng)在抽氣控制方面的研究還尚待完善,如只能實現(xiàn)恒速抽氣,沒有建立完善的數(shù)學(xué)模型,無法應(yīng)對諸如外界氣壓變化、管路內(nèi)部氣源不穩(wěn)、抽氣過快或過慢等不利因素引起的燃燒不穩(wěn)定情況。本文的研究內(nèi)容就是建立抽氣系統(tǒng)的準(zhǔn)確數(shù)學(xué)模型,并對模型進行驗證,為后期設(shè)計數(shù)字控制器從而實現(xiàn)抽氣系統(tǒng)變速抽氣控制做準(zhǔn)備。

    1 抽氣系統(tǒng)工作原理

    抽氣系統(tǒng)結(jié)構(gòu)如圖1所示,系統(tǒng)由干燥管、緩沖瓶、尾氣收集裝置、壓差傳感器、數(shù)據(jù)采集單元、單片機電路以及控制程序組成。系統(tǒng)工作時,操作人員通過控制程序設(shè)置好抽氣裝置運行參數(shù),發(fā)布抽氣命令,單片機接收到參數(shù)和命令后,控制尾氣收集裝置運行的相關(guān)電路導(dǎo)通,尾氣收集裝置按照設(shè)定的運行速度將燃燒室排氣口出來的尾氣經(jīng)過干燥后抽入尾氣收集裝置的氣缸內(nèi);壓差傳感器檢測到燃燒室進、排氣口兩端壓差的變化并輸出電壓,使用Agilent 34972A數(shù)據(jù)采集單元采集該電壓,通過串口將電壓發(fā)送到控制程序中再換算為氣體壓差值,經(jīng)過控制器的計算后給出下一時刻的運行速度信號,再通過串口自動向單片機發(fā)布新的運行速度信號;如此循環(huán)往復(fù),在控制程序的控制下實現(xiàn)運行速度的自調(diào)整,在一次燃燒實驗完成后,排出氣缸內(nèi)的尾氣。

    圖1 抽氣系統(tǒng)結(jié)構(gòu)示意圖

    2 實驗信號設(shè)計

    2.1 M序列生成

    為了預(yù)先對抽氣系統(tǒng)的動態(tài)特性有一定了解,如確定抽氣系統(tǒng)是否存在滯后現(xiàn)象,確定系統(tǒng)調(diào)節(jié)過程時間等,對抽氣系統(tǒng)進行階躍響應(yīng)測試[11]。階躍響應(yīng)測試時,抽氣系統(tǒng)工作在開環(huán)控制狀態(tài)。設(shè)定尾氣收集裝置活塞運行速度為20 mm/min,待抽氣系統(tǒng)運行處于穩(wěn)定抽氣狀態(tài)時,保持當(dāng)前運行速度,在此運行速度基礎(chǔ)上再突加一個幅值為β的正向速度信號,作為階躍信號輸入抽氣系統(tǒng)。為確定多大幅值的β值能獲得較為理想的階躍響應(yīng)信號,在20 mm/min運行速度基礎(chǔ)上分別施加原有運行速度的5%、10%、15%、20%、25%、50%作為階躍信號輸入系統(tǒng),抽氣系統(tǒng)在不同幅值階躍信號測試下的響應(yīng)曲線如圖2所示。實驗結(jié)果表明,當(dāng)階躍信號幅值β選取為原來速度的25%,即β=5 mm/min時,系統(tǒng)階躍響應(yīng)效果較為理想。

    圖2 不同幅值階躍信號測試下系統(tǒng)的響應(yīng)曲線

    在得到不同幅值階躍信號測試結(jié)果后,將β的幅值設(shè)定為5 mm/min再次進行測試,記錄數(shù)據(jù)如圖3所示。多次重復(fù)實驗結(jié)果表明,抽氣系統(tǒng)的調(diào)節(jié)時間Ts≈21 s,純滯后時間 τ0≈2 s。

    圖3 5 mm/min階躍信號測試下系統(tǒng)的響應(yīng)曲線

    沖激信號作為最理想的系統(tǒng)測試信號,擁有廣泛的頻譜,但讓計算機產(chǎn)生一個理想的沖激信號較為困難。研究發(fā)現(xiàn),白噪聲的自相關(guān)函數(shù)與沖激信號很接近。白噪聲取代沖激信號對一個系統(tǒng)進行測試是最為理想的,但在實際工程應(yīng)用中,執(zhí)行器并不可能完全按照白噪聲的變化規(guī)律來動作,而M序列(即最長線性移位寄存器序列)能較好地模擬執(zhí)行器動態(tài)運行中的輸入信號。通過對M序列自相關(guān)分析發(fā)現(xiàn),M序列自相關(guān)函數(shù)也接近于沖激信號,其自相關(guān)函數(shù)圖形分析如圖4所示。

    圖4 M序列自相關(guān)函數(shù)圖形

    使用M序列產(chǎn)生一系列‘0’和‘1’的碼值,并給碼值‘0’和碼值‘1’賦以相應(yīng)的活塞運行速度值作為抽氣系統(tǒng)激勵信號輸入系統(tǒng),得到一系列M序列的響應(yīng)值。利用移位寄存器產(chǎn)生‘0’和‘1’的碼值,需要確定M序列的時鐘周期T0、序列長度Np、序列幅值A(chǔ)0和試驗長度N。已知抽氣系統(tǒng)的調(diào)節(jié)時間Ts,M序列的時鐘周期可根據(jù)式(1)來確定:

    序列長度根據(jù)式(2)確定:

    代入數(shù)據(jù)計算后取M序列時鐘周期T0=0.5 s,序列長度Np=63。依據(jù)Np=26–1=63,確定M序列由6級移位寄存器發(fā)生,6級移位寄存器連接方式用式(3)表示,寄存器連接示意圖如圖5所示。

    圖5 6級移位寄存器連接示意圖

    根據(jù)抽氣系統(tǒng)階躍響應(yīng)情況可知,系統(tǒng)初始運行時會有較長的調(diào)節(jié)時間,在獲取動態(tài)數(shù)據(jù)時,需連續(xù)發(fā)送兩個M序列實驗信號。第一個M序列用來消除調(diào)節(jié)過程,得到的系統(tǒng)響應(yīng)數(shù)據(jù)不用于參數(shù)辨識,系統(tǒng)的第二個M序列響應(yīng)數(shù)據(jù)用于參數(shù)辨識。依據(jù)6級移位寄存器結(jié)構(gòu)示意圖,利用LabVIEW編寫M序列碼生成程序[12],初始序列設(shè)為101010,產(chǎn)生長度為 126,輸出‘0’和‘1’的 M 序列值。為保證每一個碼值都能有效的輸入抽氣系統(tǒng),程序?qū)序列碼值以一定時間間隔逐點連續(xù)輸出。

    選取M序列幅值A(chǔ)0=±10 mm/min,在原設(shè)定運行速度20 mm/min基礎(chǔ)上,當(dāng)M序列產(chǎn)生碼值‘1’時 賦 30 mm/min,碼 值 ‘0’時 賦 10 mm/min。 在LabVIEW程序中將速度命令按“S3××W”格式以字符串形式發(fā)送至下位機,其中“××”即為賦值的運行速度值。長度為126的M序列01碼值和生成程序如圖6和圖7所示。

    圖6 M序列01碼幅值圖

    圖7 M序列產(chǎn)生程序

    2.2 直流分量處理

    對采樣得到的抽氣系統(tǒng)響應(yīng)信號進行頻譜分析,結(jié)果表明,輸出信號中包含有直流分量。為了使建立的模型能充分反映系統(tǒng)的動態(tài)特性,因此必須把信號中的直流分量去除。圖8顯示了在輸入M序列后抽氣系統(tǒng)響應(yīng)信號在去除直流分量前的波形。將采集得到的輸出時域信號用快速傅里葉變換(fast Fourier transformation,簡稱 FFT)法變換到頻域信號,將信號中頻率為0的分量置為0后,再將信號進行傅里葉逆變換,即可得到不含直流分量的輸出信號,再經(jīng)過濾波和平滑處理得到圖9所示的波形。此時的輸出信號值用作系統(tǒng)數(shù)學(xué)模型的參數(shù)辨識。

    圖8 響應(yīng)信號去除直流分量前波形和FFT波形

    圖9 響應(yīng)信號去除直流分量后波形和FFT波形

    3 梯度校正法

    梯度校正法進行參數(shù)辨識的基本思想是沿著準(zhǔn)則函數(shù)的負(fù)梯度方向,逐漸修正模型參數(shù)的估計值,直至準(zhǔn)則函數(shù)達(dá)到最小值[13-15]。算法結(jié)構(gòu)是在舊的估計值基礎(chǔ)上,將增益矩陣和新息的乘積求和而得到新的估計值。梯度校正法具有算法簡單、計算量小等優(yōu)點。

    已知抽氣系統(tǒng)的模型階次na=nb=2,真時延d=2,在不考慮噪聲的情況下,設(shè)抽氣系統(tǒng)模型方程的差分形式為:

    將式(4)改寫成矩陣形式為:

    其中:

    式中:y(k)——抽氣系統(tǒng)的輸出值;

    u(k)——輸入值;

    k——采樣次數(shù);

    a1,a2,b1,b2——待辨識參數(shù);

    h(k)——測量矩陣;

    θ——參數(shù)矩陣。

    梯度校正法的目標(biāo)就是根據(jù)輸入輸出數(shù)據(jù)來確定參數(shù) θ和h(k),并以此確定估計值使得準(zhǔn)則函數(shù)J(θ)沿著負(fù)梯度方向不斷修正的值,直到J(θ)達(dá)到最小值。

    其中:

    R(k)——N維加權(quán)矩陣,且N=4;

    ΛN(k)——對角陣元素。

    將去除直流分量后的輸出數(shù)據(jù)按梯度校正法進行參數(shù)辨識,辨識算法流程如圖10所示。使用Matlab編寫模型辨識程序,程序開始執(zhí)行后,程序會先給待辨識參數(shù) θ和h(k)賦初值;按照式(8)計算加權(quán)矩陣R(k);下一步根據(jù)J(θ)表達(dá)式和表達(dá)式計算出估計值;之后計算出模型輸出值y(k),再計算系統(tǒng)輸出值y(k)與模型輸出值ym(k)之間的偏差ey;之后再計算準(zhǔn)則函數(shù)的梯度,判斷是否達(dá)到設(shè)定循環(huán)次數(shù),如果沒有則回到權(quán)矩陣計算那一步,直到程序運行到設(shè)定的循環(huán)次數(shù);達(dá)到設(shè)定循環(huán)次數(shù)后,將各項數(shù)值輸出,繪制圖形,結(jié)束程序。

    圖10 梯度校正法參數(shù)辨識算法流程圖

    待辨識參數(shù)梯度校正變化情況如圖11所示,系統(tǒng)輸出值與模型輸出值偏差如圖12所示。圖11表明,當(dāng)辨識過程進行到11次時,待辨識參數(shù)基本達(dá)到穩(wěn)定狀態(tài),圖12也表明此時系統(tǒng)輸出值與模型輸出值的誤差基本達(dá)到穩(wěn)定狀態(tài)。參數(shù)辨識結(jié)果顯示,待辨識參數(shù)a1=0.1185,a2=0.0569,b1= ?0.013 2,b2=0.0031。由此得到抽氣系統(tǒng)的數(shù)學(xué)模型為:

    圖11 待辨識參數(shù)梯度校正變化情況

    圖12 系統(tǒng)輸出值與模型輸出值偏差

    4 數(shù)學(xué)模型的實驗驗證

    抽氣系統(tǒng)是燃?xì)鉄嶂禍y量裝置的組成部分,位于石英玻璃燃燒器的下游,由尾氣收集裝置、管道、測控器件等部分組成,圖13是尾氣收集裝置和石英玻璃燃燒器的實物圖。圖中左側(cè)的尾氣收集裝置,通過伺服電機驅(qū)動運動塊上下移動的方式將燃燒室排出的尾氣抽入缸內(nèi)。系統(tǒng)工作時只需要控制尾氣收集裝置內(nèi)活塞移動速度即可控制玻璃燃燒室進、排氣口兩端壓差。

    圖13 尾氣收集裝置(左)和石英玻璃燃燒器(右)實物圖

    為驗證數(shù)學(xué)模型的準(zhǔn)確性,分別進行了兩組抽氣實驗,第一組是抽氣系統(tǒng)同時通入甲烷、氧氣、氬氣進行抽氣,第二組是供氣接口斷開抽取空氣,之后分別將兩組實驗中抽氣時系統(tǒng)的輸出與數(shù)學(xué)模型仿真的輸出進行對比。圖14是用LabVIEW中控制與仿真模塊下的PID子模塊結(jié)合抽氣系統(tǒng)數(shù)學(xué)模型編寫程序進行仿真。

    圖14 數(shù)學(xué)模型LabVIEW仿真程序

    在抽氣實驗中,設(shè)定條件如下:

    1)以系統(tǒng)開始供氣作為0時刻,供氣時使用甲烷、氧氣、氬氣三種氣體,其體積流量比例為65∶156∶240。

    2)尾氣收集裝置運動塊運行的初始速度設(shè)為20 mm/min。

    3)將期望的壓差值設(shè)為35 Pa,PID控制器中設(shè)置比例系數(shù)Kp=4.2219,積分系數(shù)Ki=0.01407,微分系數(shù)Kd=0.000282。

    三種氣體供氣抽氣實驗和抽空氣實驗輸出壓差變化情況與LabVIEW仿真程序輸出壓差變化情況分別如圖15、圖16所示。從兩圖可以看出,抽氣實驗輸出和數(shù)學(xué)模型仿真輸出均在20 s左右達(dá)到設(shè)定值附近。從抽氣實驗輸出波形和數(shù)學(xué)模型輸出波形達(dá)到較為穩(wěn)定時來看,二者輸出波形較為接近。三種氣體供氣抽氣實驗輸出波形出現(xiàn)波動,是因為供氣不穩(wěn)導(dǎo)致,但輸出波形仍在設(shè)定壓差上下波動。模型輸出和抽氣實驗的輸出結(jié)果對比表明,數(shù)學(xué)模型能比較準(zhǔn)確的反映系統(tǒng)輸入輸出情況。

    圖15 數(shù)學(xué)模型輸出壓差和三種氣體抽氣實驗輸出壓差

    圖16 數(shù)學(xué)模型輸出壓差和抽空氣實驗輸出壓差

    圖17反映抽空氣實驗和三種氣體抽氣實驗工況下的模型誤差變化情況。抽空氣時,數(shù)學(xué)模型輸出結(jié)果與抽氣實驗輸出結(jié)果的平均相對誤差為0.42%,壓差值平均絕對誤差約為0.15 Pa;三種氣體抽氣實驗時,數(shù)學(xué)模型輸出結(jié)果與抽氣實驗結(jié)果的平均相對誤差為0.45%,壓差值平均絕對誤差約為0.16 Pa。

    圖17 模型相對誤差變化圖

    5 結(jié)束語

    以梯度校正法為基本原理,理論和實驗相結(jié)合,在LabVIEW和Matlab等工具軟件的輔助下建立了燃?xì)鉄嶂禍y量裝置的抽氣系統(tǒng)數(shù)學(xué)模型。以梯度校正法的數(shù)學(xué)原理設(shè)計辨識程序,具有代碼量小、計算速度快、通用程度高等優(yōu)點。文中首次通過建立抽氣系統(tǒng)數(shù)學(xué)模型的方式研究燃?xì)鉄嶂禍y量裝置抽氣控制問題;以燃燒室進、排氣口壓差作為反饋信號,實現(xiàn)了抽氣系統(tǒng)的恒壓控制。抽氣系統(tǒng)數(shù)學(xué)模型仿真結(jié)果表明,在PID控制下模型輸出結(jié)果與實際抽氣控制傳感器測得結(jié)果基本吻合。在抽空氣實驗中,數(shù)學(xué)模型與實驗輸出平均相對誤差為0.42%;在抽三種氣體實驗中,數(shù)學(xué)模型與實驗輸出平均相對誤差為0.45%。表明文中所采用梯度校正法建立抽氣系統(tǒng)數(shù)學(xué)模型的正確性,能滿足實際控制需求,為以后設(shè)計抽氣系統(tǒng)的數(shù)字控制器提供了理論依據(jù)。

    猜你喜歡
    抽氣階躍壓差
    基于階躍雙包層光纖的螺旋型光纖傳感器
    抽氣通道性能理論研究
    科技視界(2021年4期)2021-04-13 06:03:58
    在線處理脫丁烷塔壓差高
    探討單位階躍信號的教學(xué)
    某型燃?xì)廨啓C第5級外抽氣流量估算
    機械工程師(2015年9期)2015-02-26 08:38:24
    二維氣體動力學(xué)中壓差方程的特征分解和簡單波
    基于改進經(jīng)濟壓差算法的主站AVC的研究
    淺析濾紙透氣度與初始壓差的關(guān)系
    汽車零部件(2014年2期)2014-03-11 17:46:34
    閉合性氣胸采用深靜脈留置針行閉式引流與胸穿抽氣的療效比較
    一種階躍函數(shù)在矩形時間窗口頻域特性的分析方法
    97在线视频观看| 国产在线免费精品| 午夜久久久在线观看| 免费大片黄手机在线观看| 美女国产高潮福利片在线看| 欧美日韩在线观看h| 国内精品宾馆在线| 久久精品久久久久久噜噜老黄| 午夜视频国产福利| 精品久久久精品久久久| 欧美日韩在线观看h| 久久久久久久国产电影| 久久ye,这里只有精品| 成年人免费黄色播放视频| 午夜91福利影院| 精品国产一区二区久久| 在线观看免费视频网站a站| 大片免费播放器 马上看| av免费观看日本| 亚洲av不卡在线观看| 18禁在线播放成人免费| 成年人免费黄色播放视频| 亚洲成人一二三区av| a级毛片免费高清观看在线播放| 国产高清不卡午夜福利| 狂野欧美白嫩少妇大欣赏| 国产国拍精品亚洲av在线观看| 尾随美女入室| 天堂俺去俺来也www色官网| 免费大片18禁| 亚洲精品,欧美精品| 国产精品人妻久久久久久| 亚洲精品第二区| 国国产精品蜜臀av免费| 欧美xxⅹ黑人| 亚洲第一av免费看| 国产亚洲精品久久久com| 亚洲一级一片aⅴ在线观看| 国产成人精品在线电影| 老司机亚洲免费影院| 免费av不卡在线播放| 最近手机中文字幕大全| 免费播放大片免费观看视频在线观看| 人妻制服诱惑在线中文字幕| 麻豆乱淫一区二区| 亚洲精品国产色婷婷电影| 久久免费观看电影| 国产免费又黄又爽又色| 热re99久久国产66热| 男人操女人黄网站| 熟女av电影| 老司机亚洲免费影院| 五月玫瑰六月丁香| 永久网站在线| 日韩,欧美,国产一区二区三区| 欧美少妇被猛烈插入视频| 母亲3免费完整高清在线观看 | av视频免费观看在线观看| 99国产综合亚洲精品| 2021少妇久久久久久久久久久| 99久久中文字幕三级久久日本| 伦理电影大哥的女人| 国产亚洲精品第一综合不卡 | 成人手机av| 亚洲第一区二区三区不卡| 涩涩av久久男人的天堂| 涩涩av久久男人的天堂| 蜜桃国产av成人99| 亚洲欧美一区二区三区国产| 国精品久久久久久国模美| 制服诱惑二区| 我的老师免费观看完整版| 中国美白少妇内射xxxbb| 亚洲精品乱码久久久v下载方式| 最近中文字幕高清免费大全6| www.色视频.com| 国产精品国产三级专区第一集| 日韩,欧美,国产一区二区三区| 久久这里有精品视频免费| 日韩 亚洲 欧美在线| 中文精品一卡2卡3卡4更新| 成人毛片60女人毛片免费| 免费少妇av软件| 亚洲av综合色区一区| 91在线精品国自产拍蜜月| 日日撸夜夜添| 交换朋友夫妻互换小说| videosex国产| 国产av一区二区精品久久| 国产一级毛片在线| 下体分泌物呈黄色| 成人无遮挡网站| 日韩中文字幕视频在线看片| 精品人妻偷拍中文字幕| 欧美变态另类bdsm刘玥| 中文字幕免费在线视频6| 欧美精品亚洲一区二区| 青春草亚洲视频在线观看| 国产极品天堂在线| 亚洲在久久综合| 十八禁高潮呻吟视频| 成年人午夜在线观看视频| 精品卡一卡二卡四卡免费| 日韩不卡一区二区三区视频在线| 国产高清不卡午夜福利| 日韩电影二区| 亚洲精品乱码久久久v下载方式| 下体分泌物呈黄色| 天美传媒精品一区二区| 午夜av观看不卡| 亚洲精品aⅴ在线观看| 十八禁网站网址无遮挡| 免费黄频网站在线观看国产| 晚上一个人看的免费电影| 日本欧美国产在线视频| 成人影院久久| 国产国语露脸激情在线看| 国产精品.久久久| 99久久人妻综合| 国产精品嫩草影院av在线观看| 日韩av免费高清视频| 视频区图区小说| 精品人妻偷拍中文字幕| 成人亚洲精品一区在线观看| 日韩,欧美,国产一区二区三区| av专区在线播放| 伊人久久国产一区二区| 国产一级毛片在线| 熟女人妻精品中文字幕| 国产在视频线精品| 亚洲精品国产色婷婷电影| 亚洲怡红院男人天堂| 国产精品 国内视频| 国产精品一区二区三区四区免费观看| 亚洲丝袜综合中文字幕| 亚洲av综合色区一区| 久久精品国产亚洲网站| 丰满迷人的少妇在线观看| 黄片无遮挡物在线观看| 日韩,欧美,国产一区二区三区| 亚洲五月色婷婷综合| 黄色怎么调成土黄色| 有码 亚洲区| 国产日韩欧美亚洲二区| 51国产日韩欧美| 91精品伊人久久大香线蕉| 哪个播放器可以免费观看大片| 狠狠婷婷综合久久久久久88av| 亚洲第一av免费看| 下体分泌物呈黄色| 久久影院123| 日韩伦理黄色片| 一级,二级,三级黄色视频| 成年人免费黄色播放视频| 考比视频在线观看| 久热久热在线精品观看| 曰老女人黄片| 最后的刺客免费高清国语| 制服诱惑二区| 我的老师免费观看完整版| 亚洲中文av在线| 亚洲图色成人| 亚洲国产最新在线播放| 亚州av有码| 国产高清不卡午夜福利| 看十八女毛片水多多多| 免费看光身美女| 女性生殖器流出的白浆| 欧美老熟妇乱子伦牲交| 欧美激情国产日韩精品一区| 亚洲精品一二三| 欧美变态另类bdsm刘玥| 熟女av电影| 精品少妇黑人巨大在线播放| 九草在线视频观看| 国产在线免费精品| 久热这里只有精品99| 乱码一卡2卡4卡精品| 99热这里只有精品一区| 夜夜骑夜夜射夜夜干| 美女xxoo啪啪120秒动态图| 亚洲一级一片aⅴ在线观看| 精品久久久久久电影网| 桃花免费在线播放| 日韩中字成人| 日韩不卡一区二区三区视频在线| 国产欧美日韩综合在线一区二区| 一区在线观看完整版| 精品人妻偷拍中文字幕| 插逼视频在线观看| 国产成人精品在线电影| 国产成人精品福利久久| 国产成人一区二区在线| 97超碰精品成人国产| 久久久精品免费免费高清| 午夜影院在线不卡| 夜夜爽夜夜爽视频| 夜夜爽夜夜爽视频| tube8黄色片| 久久久久久久久久久丰满| 久久久久久久大尺度免费视频| 亚洲av不卡在线观看| a级毛色黄片| 久久午夜福利片| 成人漫画全彩无遮挡| 黑人高潮一二区| 女人精品久久久久毛片| 久久久久久久久大av| 欧美精品一区二区免费开放| 免费观看av网站的网址| 一个人看视频在线观看www免费| 精品人妻在线不人妻| 日韩欧美一区视频在线观看| 99久久综合免费| 亚洲国产精品一区三区| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人爽人人爽人人片va| 精品一区二区免费观看| 男女啪啪激烈高潮av片| 国产免费现黄频在线看| 国产精品99久久久久久久久| 亚洲,欧美,日韩| 国产不卡av网站在线观看| 久热久热在线精品观看| 亚洲欧美一区二区三区国产| 精品久久久精品久久久| 人妻人人澡人人爽人人| 在现免费观看毛片| 国产免费一区二区三区四区乱码| 老司机影院成人| 日韩亚洲欧美综合| 国产精品一区www在线观看| 能在线免费看毛片的网站| 久久99热这里只频精品6学生| 丰满迷人的少妇在线观看| 午夜免费鲁丝| 免费高清在线观看视频在线观看| 人人妻人人澡人人爽人人夜夜| 国产在线视频一区二区| 成年人免费黄色播放视频| 亚洲成色77777| 色94色欧美一区二区| 天天躁夜夜躁狠狠久久av| 一本色道久久久久久精品综合| 热re99久久国产66热| av国产久精品久网站免费入址| 热99国产精品久久久久久7| 亚洲国产精品成人久久小说| 国产深夜福利视频在线观看| 国产精品国产av在线观看| 秋霞伦理黄片| 日本-黄色视频高清免费观看| 制服丝袜香蕉在线| 欧美精品人与动牲交sv欧美| 亚洲精品自拍成人| 91国产中文字幕| 久久免费观看电影| 肉色欧美久久久久久久蜜桃| 黑丝袜美女国产一区| 伦理电影大哥的女人| 精品酒店卫生间| av.在线天堂| 亚洲av欧美aⅴ国产| 欧美精品一区二区免费开放| 亚洲精品日韩av片在线观看| 国产成人精品婷婷| 成年人午夜在线观看视频| 国产爽快片一区二区三区| 五月开心婷婷网| 欧美日韩在线观看h| 国产亚洲午夜精品一区二区久久| 简卡轻食公司| 国产色婷婷99| 亚洲av欧美aⅴ国产| √禁漫天堂资源中文www| 国产亚洲一区二区精品| 黄色一级大片看看| 亚洲精品,欧美精品| 国产精品国产av在线观看| 91午夜精品亚洲一区二区三区| 狂野欧美激情性xxxx在线观看| 亚洲av综合色区一区| 精品视频人人做人人爽| 欧美xxⅹ黑人| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久久精品电影小说| 老熟女久久久| 一区二区日韩欧美中文字幕 | 亚洲中文av在线| 欧美精品一区二区大全| www.色视频.com| 久久午夜综合久久蜜桃| 丝袜美足系列| 九色亚洲精品在线播放| 久久精品久久久久久久性| 丝袜喷水一区| kizo精华| 亚洲欧美精品自产自拍| 男人操女人黄网站| 国产精品人妻久久久影院| 日韩强制内射视频| 欧美 亚洲 国产 日韩一| 狂野欧美激情性bbbbbb| 老熟女久久久| 伦精品一区二区三区| 永久免费av网站大全| 美女主播在线视频| 国产精品三级大全| 99久久综合免费| 一边摸一边做爽爽视频免费| 国产色爽女视频免费观看| 国产在线免费精品| 亚洲精品亚洲一区二区| 成人国产av品久久久| 建设人人有责人人尽责人人享有的| 丝瓜视频免费看黄片| 亚洲性久久影院| 91午夜精品亚洲一区二区三区| 久久久精品94久久精品| 成年av动漫网址| 久久久久久久久久成人| 久久免费观看电影| 七月丁香在线播放| 亚洲高清免费不卡视频| 亚洲精品456在线播放app| 精品亚洲成a人片在线观看| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕 | 亚洲国产精品999| 成人毛片60女人毛片免费| 十分钟在线观看高清视频www| 免费大片18禁| 三级国产精品片| 在线观看一区二区三区激情| 黑丝袜美女国产一区| 日本vs欧美在线观看视频| 色视频在线一区二区三区| 美女内射精品一级片tv| 超碰97精品在线观看| 亚洲av成人精品一区久久| 黄色视频在线播放观看不卡| 国产老妇伦熟女老妇高清| 国产在线一区二区三区精| 天天操日日干夜夜撸| 久久ye,这里只有精品| 亚洲av综合色区一区| 国产成人精品久久久久久| 18禁在线播放成人免费| xxxhd国产人妻xxx| 最近手机中文字幕大全| 免费黄网站久久成人精品| 亚洲中文av在线| 99热国产这里只有精品6| 9色porny在线观看| 国产成人aa在线观看| 精品久久久噜噜| 亚洲欧美一区二区三区黑人 | 国产日韩欧美视频二区| 人妻 亚洲 视频| 国产成人精品无人区| 国产黄片视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| freevideosex欧美| av一本久久久久| 欧美成人精品欧美一级黄| 国产成人一区二区在线| 国产熟女午夜一区二区三区 | 久热这里只有精品99| 午夜日本视频在线| 久久午夜综合久久蜜桃| 婷婷色av中文字幕| 插逼视频在线观看| 午夜久久久在线观看| 午夜精品国产一区二区电影| 免费av不卡在线播放| 亚洲欧洲国产日韩| 一区二区三区乱码不卡18| 精品久久久噜噜| 成人毛片a级毛片在线播放| av又黄又爽大尺度在线免费看| 婷婷色综合www| 少妇精品久久久久久久| 永久网站在线| 欧美精品国产亚洲| 中文乱码字字幕精品一区二区三区| 黄色欧美视频在线观看| 高清黄色对白视频在线免费看| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美清纯卡通| av线在线观看网站| 国产白丝娇喘喷水9色精品| 精品国产一区二区久久| 精品国产国语对白av| 91成人精品电影| 男女免费视频国产| 一边亲一边摸免费视频| 久久久久精品性色| 亚洲人成网站在线观看播放| 日产精品乱码卡一卡2卡三| 一级二级三级毛片免费看| 亚洲精品456在线播放app| 少妇高潮的动态图| www.色视频.com| 美女cb高潮喷水在线观看| 亚洲欧美色中文字幕在线| 满18在线观看网站| 多毛熟女@视频| 婷婷成人精品国产| 赤兔流量卡办理| 欧美亚洲 丝袜 人妻 在线| 色视频在线一区二区三区| 亚洲人成网站在线播| 美女脱内裤让男人舔精品视频| 各种免费的搞黄视频| 丝袜脚勾引网站| 80岁老熟妇乱子伦牲交| 久久99精品国语久久久| 色视频在线一区二区三区| 女的被弄到高潮叫床怎么办| 欧美最新免费一区二区三区| 97超视频在线观看视频| 五月天丁香电影| 欧美变态另类bdsm刘玥| 只有这里有精品99| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 亚洲精品成人av观看孕妇| 十八禁网站网址无遮挡| 国产极品粉嫩免费观看在线 | 国产 精品1| 热99久久久久精品小说推荐| 欧美人与性动交α欧美精品济南到 | 狠狠婷婷综合久久久久久88av| 久久久国产一区二区| 国产精品国产三级国产专区5o| 精品酒店卫生间| 国产一区亚洲一区在线观看| av不卡在线播放| 日本91视频免费播放| 国产成人一区二区在线| 国产伦精品一区二区三区视频9| 一本久久精品| 亚洲精品久久成人aⅴ小说 | 精品熟女少妇av免费看| 自拍欧美九色日韩亚洲蝌蚪91| 插阴视频在线观看视频| kizo精华| 最近最新中文字幕免费大全7| 久久精品国产亚洲网站| 久久久久久久久久久免费av| 国产男女内射视频| 国产精品国产三级国产专区5o| 性高湖久久久久久久久免费观看| 亚洲av福利一区| 婷婷成人精品国产| 中文字幕人妻丝袜制服| 热re99久久国产66热| 午夜视频国产福利| 高清av免费在线| 在线观看人妻少妇| 91国产中文字幕| 2021少妇久久久久久久久久久| 国产伦理片在线播放av一区| 日日撸夜夜添| 51国产日韩欧美| .国产精品久久| 狂野欧美激情性xxxx在线观看| 国产在线一区二区三区精| 91午夜精品亚洲一区二区三区| 国产成人免费无遮挡视频| videos熟女内射| 日韩av在线免费看完整版不卡| 91久久精品国产一区二区成人| 久久久久久久亚洲中文字幕| 亚洲色图 男人天堂 中文字幕 | 蜜臀久久99精品久久宅男| 亚洲国产最新在线播放| 日韩欧美一区视频在线观看| 国产国拍精品亚洲av在线观看| 18在线观看网站| 中文字幕av电影在线播放| 日韩一区二区视频免费看| 国产极品粉嫩免费观看在线 | 国产精品麻豆人妻色哟哟久久| 日产精品乱码卡一卡2卡三| 亚洲欧美清纯卡通| 日韩在线高清观看一区二区三区| 少妇人妻 视频| 99热全是精品| 最近中文字幕2019免费版| 成人影院久久| 久久久国产欧美日韩av| 亚洲精品日韩在线中文字幕| 大片免费播放器 马上看| 成人综合一区亚洲| 3wmmmm亚洲av在线观看| 热99国产精品久久久久久7| 久久精品国产亚洲av涩爱| 成人国语在线视频| 97精品久久久久久久久久精品| 黄色配什么色好看| 啦啦啦在线观看免费高清www| 午夜精品国产一区二区电影| 天天躁夜夜躁狠狠久久av| 日韩欧美精品免费久久| 久热这里只有精品99| 少妇人妻久久综合中文| 国产黄色视频一区二区在线观看| 人妻制服诱惑在线中文字幕| 国产亚洲午夜精品一区二区久久| 亚洲国产成人一精品久久久| 亚洲av中文av极速乱| 精品国产国语对白av| 人妻制服诱惑在线中文字幕| 一边亲一边摸免费视频| 午夜福利视频在线观看免费| 五月伊人婷婷丁香| 狂野欧美白嫩少妇大欣赏| 欧美xxⅹ黑人| av国产久精品久网站免费入址| 十八禁网站网址无遮挡| 激情五月婷婷亚洲| 欧美一级a爱片免费观看看| 在线观看三级黄色| 亚洲不卡免费看| av在线观看视频网站免费| 午夜激情福利司机影院| 国产极品天堂在线| 久久99热6这里只有精品| 久久久久久久久久久久大奶| 在线观看美女被高潮喷水网站| 秋霞伦理黄片| xxx大片免费视频| 九色亚洲精品在线播放| 肉色欧美久久久久久久蜜桃| 免费观看无遮挡的男女| 久久av网站| 狂野欧美激情性xxxx在线观看| 秋霞伦理黄片| 少妇的逼好多水| 中文字幕久久专区| 97超视频在线观看视频| 99久久人妻综合| 亚洲欧美成人综合另类久久久| 国产亚洲一区二区精品| 日韩av在线免费看完整版不卡| 亚洲四区av| 国产av精品麻豆| 91久久精品国产一区二区成人| 久久亚洲国产成人精品v| 日韩av免费高清视频| 少妇精品久久久久久久| 各种免费的搞黄视频| 中文字幕人妻熟人妻熟丝袜美| 99久久中文字幕三级久久日本| 国产在线一区二区三区精| 少妇丰满av| 亚洲人成网站在线播| 免费人妻精品一区二区三区视频| 亚洲欧美中文字幕日韩二区| 久久女婷五月综合色啪小说| 国产高清有码在线观看视频| a级片在线免费高清观看视频| 有码 亚洲区| h视频一区二区三区| 桃花免费在线播放| 亚洲不卡免费看| 26uuu在线亚洲综合色| 国产深夜福利视频在线观看| 精品国产一区二区久久| 超碰97精品在线观看| 哪个播放器可以免费观看大片| 特大巨黑吊av在线直播| 国产一区二区三区av在线| 国产69精品久久久久777片| 成人黄色视频免费在线看| 啦啦啦啦在线视频资源| 久久99蜜桃精品久久| 一边摸一边做爽爽视频免费| 亚洲一区二区三区欧美精品| 亚洲五月色婷婷综合| 日产精品乱码卡一卡2卡三| 免费看av在线观看网站| 亚洲熟女精品中文字幕| 999精品在线视频| av免费观看日本| 亚洲精品国产av蜜桃| 国产极品粉嫩免费观看在线 | 又黄又爽又刺激的免费视频.| 少妇人妻久久综合中文| 精品人妻一区二区三区麻豆| 日本vs欧美在线观看视频| 久久精品国产a三级三级三级| av免费观看日本| 18禁动态无遮挡网站| 看非洲黑人一级黄片| 18禁观看日本| 欧美bdsm另类| 免费黄网站久久成人精品| 极品人妻少妇av视频| 老熟女久久久| 桃花免费在线播放| 老女人水多毛片| 免费不卡的大黄色大毛片视频在线观看| videos熟女内射| 国产精品人妻久久久影院| 夜夜爽夜夜爽视频| 精品久久国产蜜桃| 国产午夜精品一二区理论片| 制服人妻中文乱码| 午夜福利,免费看| 午夜福利视频精品| 在线看a的网站|