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

    非線性燃燒不穩(wěn)定振蕩幅值演化特征信號(hào)辨識(shí)方法研究①

    2023-04-26 02:04:50張翔宇金秉寧魏少娟劉佩進(jìn)
    固體火箭技術(shù) 2023年1期
    關(guān)鍵詞:適應(yīng)度幅值模態(tài)

    甘 林,張翔宇,金秉寧*,魏少娟,劉佩進(jìn)

    (1.西北工業(yè)大學(xué) 燃燒、熱結(jié)構(gòu)與內(nèi)流場(chǎng)重點(diǎn)實(shí)驗(yàn)室,西安 710072;2.中國(guó)航天科技集團(tuán)有限公司四院四十一所,西安 710025;3.空軍工程大學(xué) 裝備管理與無(wú)人機(jī)工程學(xué)院,西安 710000)

    0 引言

    非線性燃燒不穩(wěn)定問(wèn)題一直是固體火箭發(fā)動(dòng)機(jī)研制過(guò)程中常遇到的棘手問(wèn)題之一[1-5],其典型特征是在發(fā)動(dòng)機(jī)燃燒中產(chǎn)生非線性壓強(qiáng)振蕩現(xiàn)象[6-10]。通常,這種非線性壓強(qiáng)振蕩中包含一個(gè)或多個(gè)聲振蕩模態(tài),且各階模態(tài)受增益阻尼因素的影響,能量從低階模態(tài)向高階模態(tài)傳遞,激發(fā)出高階模態(tài),使得不穩(wěn)定問(wèn)題愈加嚴(yán)重[11-15]。然而,目前仍存在對(duì)能量傳遞演化機(jī)理認(rèn)識(shí)不清的問(wèn)題,無(wú)法對(duì)非線性燃燒不穩(wěn)定的預(yù)估和抑制提供理論指導(dǎo)。因此,需要開(kāi)展能量傳遞演化機(jī)理的研究,而進(jìn)行這項(xiàng)研究的前提則是從非線性振蕩數(shù)據(jù)中準(zhǔn)確的得到各階模態(tài)的關(guān)鍵參數(shù),如振蕩頻率和幅值等。

    對(duì)于固體火箭發(fā)動(dòng)機(jī)非線性不穩(wěn)定振蕩數(shù)據(jù),科研人員多采用商業(yè)軟件平臺(tái)進(jìn)行分析,其特點(diǎn)是功能強(qiáng)大,操作簡(jiǎn)單,經(jīng)過(guò)長(zhǎng)期探索,可以形成一套數(shù)據(jù)分析流程。但在使用過(guò)程中發(fā)現(xiàn),商業(yè)軟件中很多模塊方法和參數(shù)不明確,使得數(shù)據(jù)處理結(jié)果出現(xiàn)失真;此外,受主觀經(jīng)驗(yàn)判斷的影響,分析結(jié)果也存在一定的誤差。在處理類(lèi)似信號(hào)如電力系統(tǒng)低頻振蕩信號(hào)、機(jī)械系統(tǒng)振動(dòng)衰減信號(hào)等方面,更多的是采用時(shí)頻分析方法,如經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)[16]和變分模態(tài)分解(VMD)[17-18]等,并且取得了較好的效果。綜上,有必要借鑒處理類(lèi)似信號(hào)的分析方法建立一種方便而準(zhǔn)確的分析固體火箭發(fā)動(dòng)機(jī)非線性燃燒不穩(wěn)定振蕩數(shù)據(jù)的方法,以獲得非線性不穩(wěn)定振蕩數(shù)據(jù)的關(guān)鍵參數(shù)。

    固體火箭發(fā)動(dòng)機(jī)非線性燃燒不穩(wěn)定振蕩數(shù)據(jù)是典型的非線性信號(hào)。在非線性信號(hào)處理方面,典型的時(shí)頻方法有經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)[19]、局部均值分解(LMD)[20]、本征時(shí)間尺度分解(ITD)[21]和變分模態(tài)分解(VMD)[22]等。其中,前三種基于遞歸迭代的思想,能夠?qū)π盘?hào)模態(tài)進(jìn)行分解。然而,EMD存在著一些問(wèn)題,如缺乏嚴(yán)格的數(shù)學(xué)公式,易受模態(tài)混疊、端點(diǎn)效應(yīng)、過(guò)擬合和欠采樣的影響;LMD在抑制端點(diǎn)效應(yīng)和模態(tài)混疊方面優(yōu)于EMD,但其主要問(wèn)題是計(jì)算復(fù)雜度高,仍容易產(chǎn)生模態(tài)混疊;ITD存在末端效應(yīng)和時(shí)間分辨率低的問(wèn)題;而VMD方法基于求解一個(gè)完全非遞歸變分問(wèn)題的思想,具有理論上嚴(yán)謹(jǐn)?shù)臄?shù)學(xué)原理和對(duì)噪聲及采樣的魯棒性等優(yōu)點(diǎn)[23]。但VMD方法依賴于經(jīng)驗(yàn),需要人工預(yù)先設(shè)定模態(tài)分解個(gè)數(shù)K、二次懲罰項(xiàng)α和帶寬τ,參數(shù)設(shè)置不當(dāng)可能導(dǎo)致模態(tài)混疊,影響分解結(jié)果的準(zhǔn)確性[24]。

    綜合考慮,本文引入了VMD[24]算法,并采用標(biāo)準(zhǔn)粒子群優(yōu)化算法(SPSO 2011,以下簡(jiǎn)稱SPSO)[25]對(duì)VMD算法所需的預(yù)置參數(shù)進(jìn)行優(yōu)化,建立了SPSO-VMD方法。將SPSO-VMD方法應(yīng)用于典型特征的燃燒不穩(wěn)定聲振蕩信號(hào)處理,準(zhǔn)確地獲得了各階模態(tài)的頻率和幅值,驗(yàn)證了該方法的可行性。并將其應(yīng)用在了真實(shí)發(fā)動(dòng)機(jī)振蕩數(shù)據(jù)的處理上,準(zhǔn)確地分離出了各階模態(tài)。

    1 SPSO-VMD方法的建立

    VMD方法將原始信號(hào)X(t)分解為K(模態(tài)分解個(gè)數(shù))個(gè)緊繞其中心頻率ω的本征模態(tài)函數(shù)。在該算法中,定義本征模態(tài)函數(shù)(IMF)為一個(gè)調(diào)幅-調(diào)頻信號(hào),VMD分解過(guò)程實(shí)質(zhì)上是變分問(wèn)題的構(gòu)造和求解過(guò)程,詳見(jiàn)文獻(xiàn)[24]。

    SPSO算法是粒子群算法(Particle Swarm Optimization,PSO)[26]的一種改進(jìn)。PSO算法作為一種群體智能優(yōu)化算法,具有原理簡(jiǎn)單、易于實(shí)現(xiàn)和搜索速度快等優(yōu)勢(shì),在PSO算法中,群體中的每個(gè)個(gè)體稱為一個(gè)粒子,代表優(yōu)化問(wèn)題的一個(gè)潛在解。對(duì)于一個(gè)n維優(yōu)化問(wèn)題,每個(gè)粒子被賦予一個(gè)n維速度向量和n維位置向量。在每一次迭代也即搜索的過(guò)程中,每個(gè)粒子根據(jù)自身的搜索經(jīng)驗(yàn)和群體中其他粒子的搜索經(jīng)驗(yàn)動(dòng)態(tài)調(diào)整各自的運(yùn)動(dòng)狀態(tài),從而保證每個(gè)粒子逐步向問(wèn)題的最優(yōu)解方向靠近[27-28]。但是,基本PSO算法存在兩個(gè)缺陷:(1)不能很好地動(dòng)態(tài)調(diào)整粒子的全局和局部搜索能力;(2)粒子容易陷入搜索迭代停滯[29]。RATNAWEERA等[30]為克服PSO算法的第一個(gè)缺陷,提出了線性自適應(yīng)參數(shù)調(diào)整法則,即根據(jù)當(dāng)前迭代次數(shù)和最大迭代次數(shù),線性地調(diào)整參數(shù)的值;CLERC等[25]為克服基本PSO算法過(guò)早陷入迭代停滯的缺陷,提出了SPSO算法。

    在SPSO算法中,首先根據(jù)每個(gè)粒子的當(dāng)前位置Xi(t)、粒子的個(gè)體最優(yōu)位置Pbesti(t)以及群體的全局最優(yōu)位置Gbest,為每個(gè)粒子定義一個(gè)超球面;然后,在定義的超球面中,為每個(gè)粒子隨機(jī)產(chǎn)生一個(gè)新的位置點(diǎn);最后,將這個(gè)隨機(jī)產(chǎn)生的位置視為擾動(dòng)添加到粒子的速度更新公式中,所以在迭代初期,SPSO算法可從最大程度上保證每個(gè)粒子的飛行速度不為0??梢员M量保證在迭代的初期每個(gè)粒子都能以非零速度在解空間中一直搜索,從而避免算法過(guò)早陷入迭代停滯。

    本文按照線性自適應(yīng)參數(shù)調(diào)整法則的SPSO算法優(yōu)化VMD分解,建立了SPSO-VMD方法,其思路如下:影響VMD算法的三個(gè)需要提前輸入的參數(shù)為模態(tài)分解個(gè)數(shù)K、二次懲罰因子α和帶寬τ。在處理的過(guò)程中,模態(tài)分解個(gè)數(shù)K是通過(guò)對(duì)振蕩信號(hào)進(jìn)行傅里葉變換并濾波從而確定的;二次懲罰因子α和帶寬τ是通過(guò)SPSO算法來(lái)優(yōu)化的。

    SPSO算法優(yōu)化需要一個(gè)適應(yīng)度函數(shù)f(x),其定義如下:首先,對(duì)振蕩信號(hào)進(jìn)行傅里葉變換并濾波,得到K個(gè)模態(tài)頻率及對(duì)應(yīng)的K個(gè)幅值A(chǔ)k;再使用VMD算法對(duì)振蕩信號(hào)進(jìn)行分解得到K個(gè)模態(tài)分量,對(duì)這些模態(tài)分量進(jìn)行傅里葉變換,得到K個(gè)模態(tài)分量對(duì)應(yīng)的K個(gè)幅值A(chǔ)k,VMD(k=1,2,…,K)。各模態(tài)幅值差值ΔA計(jì)算公式為

    ΔAk=|Ak-Ak,VMD|

    ΔA={ΔA1,ΔA2,…,ΔAk}

    (1)

    將Max(ΔA)定義為SPSO算法的適應(yīng)度函數(shù)f(x),也即SPSO優(yōu)化VMD算法是朝著值Max(ΔA)變小的方向去優(yōu)化的。

    應(yīng)用SPSO算法優(yōu)化VMD算法流程如下:首先設(shè)定兩個(gè)參數(shù)α和τ優(yōu)化的具體范圍,粒子在這個(gè)范圍內(nèi)進(jìn)行全局尋優(yōu),α和τ的值由粒子的位置來(lái)表征,每個(gè)粒子每次循環(huán)將自身的位置代入VMD算法對(duì)振蕩信號(hào)進(jìn)行一次分解,并計(jì)算一次適應(yīng)度值,通過(guò)對(duì)比新粒子的適應(yīng)度值進(jìn)行位置的更新,使每個(gè)粒子向著最優(yōu)解方向更新位置,通過(guò)不斷的迭代更新,最終找到使整個(gè)種群全局適應(yīng)度最優(yōu)的位置,也即所需優(yōu)化參數(shù)的最佳值。

    SPSO-VMD方法流程圖如圖1所示。

    2 SPSO-VMD方法的驗(yàn)證

    2.1 典型特征信號(hào)構(gòu)建

    固體火箭發(fā)動(dòng)機(jī)不穩(wěn)定振蕩數(shù)據(jù)中包含振蕩增長(zhǎng)段、極限環(huán)段以及振蕩衰減段。并且在整個(gè)振蕩數(shù)據(jù)中包含多階模態(tài),各模態(tài)之間存在能量傳遞,各模態(tài)幅值變化率也隨時(shí)間變化。因此,構(gòu)造符合上述特征的典型的振蕩信號(hào)用于校驗(yàn)SPSO-VMD模態(tài)分解的準(zhǔn)確性。此時(shí),各階模態(tài)壓強(qiáng)振蕩幅值Rm滿足的非線性微分方程[13]為

    (2)

    式中αm為m階模態(tài)壓強(qiáng)振蕩的總增益系數(shù);A為非線性系數(shù);ωm為m階聲振模態(tài)的角頻率;Emnl為表示各階模態(tài)之間的非線性耦合作用。

    按上式構(gòu)造含20階模態(tài)的仿真信號(hào),前5階模態(tài)的幅值以及幅值變化率如表1所示,模態(tài)數(shù)m=20。6~20階初始幅值均為0,其幅值變化率滿足αm=(m-3)×α4。在表1中,正值αm表示能量通過(guò)線性作用從燃燒場(chǎng)流入振蕩系統(tǒng)中的第m階模態(tài),負(fù)值αm則表示能量從振蕩系統(tǒng)中流出,高階模態(tài)的αm的絕對(duì)值越大,表明高階模態(tài)越難被激發(fā)以及越容易被耗散。按上式構(gòu)造含20階模態(tài)的仿真信號(hào)如圖2所示。從圖2可看出,上升段振蕩幅值出現(xiàn)了明顯線性向非線性轉(zhuǎn)變的過(guò)程,即各階模態(tài)的幅值并不是以線性的方式持續(xù)增長(zhǎng),當(dāng)各階模態(tài)能量傳遞達(dá)到動(dòng)態(tài)平衡過(guò)程時(shí),出現(xiàn)了極限環(huán)狀態(tài);下降段各階模態(tài)的幅值變化率也是隨時(shí)變化各不相同;圖2(b)為整個(gè)振蕩段頻譜分析的前五階結(jié)果,前三階振蕩模態(tài)的幅值占主導(dǎo),所以這里僅分解前三階模態(tài)來(lái)驗(yàn)證分解結(jié)果的準(zhǔn)確性。

    表1 前五階模態(tài)的幅值以及幅值變化率

    (a)Amplitude of the simulated signal

    2.2 典型特征信號(hào)的SPSO-VMD方法分解

    對(duì)上述信號(hào)在給定α和τ范圍為(0,3000]和[0,0.2]內(nèi)進(jìn)行SPSO-VMD方法的分解,優(yōu)化所得α和τ分別為2983.78和0.000 4。各階模態(tài)幅值分解結(jié)果如圖3(b)所示。從圖3(b)可以看分解后的各階模態(tài)幅值與圖2(c)幅值輪廓線高度一致;表2展示了SPSO-VMD方法分解該特征信號(hào)的頻率和幅值結(jié)果。

    表2 SPSO-VMD方法分解特征信號(hào)頻率、幅值結(jié)果

    (a)Iterative process

    這里三個(gè)模態(tài)的分解所得極限環(huán)幅值是通過(guò)統(tǒng)計(jì)的方式計(jì)算所得,其與標(biāo)準(zhǔn)極限環(huán)幅值誤差均低于0.5%。從表中結(jié)果可以看出,對(duì)于非線性燃燒不穩(wěn)定典型特征信號(hào),SPSO-VMD方法能夠準(zhǔn)確的分離各個(gè)模態(tài)、辨識(shí)各階頻率和幅值。比較了以下兩種方法:

    (1)采用VMD算法對(duì)特征信號(hào)進(jìn)行分解,分別取值α和τ為16和0.01;

    (2)采用PSO算法以同樣思路在同樣參數(shù)范圍內(nèi)對(duì)特征信號(hào)進(jìn)行優(yōu)化VMD參數(shù)的分解,對(duì)比結(jié)果如表3和圖4所示。

    圖4 Mode 1極限環(huán)段振蕩某點(diǎn)誤差示意圖

    表3 三種算法對(duì)比

    為求出全局最小適應(yīng)度函數(shù)值f(x),采用VMD法對(duì)特征信號(hào)進(jìn)行分解。其中,α和τ的取值范圍分別為(0,3000]和[0,0.2],α的步長(zhǎng)為1,τ的步長(zhǎng)為0.01,全局求解共耗時(shí)約10 h,結(jié)果如圖5所示。

    圖5 f(x)的全局求解

    如表3和圖4所示,PSO-VMD方法和SPSO-VMD方法分解的誤差遠(yuǎn)小于VMD方法。從圖4中可知,SPSO-VMD方法的誤差略小于PSO-VMD方法,但SPSO-VMD方法可以全局優(yōu)化到與最小適應(yīng)度函數(shù)值相對(duì)應(yīng)的參數(shù),且SPSO-VMD方法所需時(shí)間更短,效率更高,如圖3(a)所示,在本次計(jì)算中,SPSO-VMD方法僅用3步就取得了比PSO方法更好的效果。通過(guò)兩者誤差的對(duì)比分析,可以認(rèn)為在適應(yīng)度函數(shù)為10-5數(shù)量級(jí)時(shí),分解便已經(jīng)達(dá)到了較好的效果。

    3 SPSO-VMD方法的應(yīng)用

    某發(fā)動(dòng)機(jī)在工作過(guò)程中發(fā)生了非線性燃燒不穩(wěn)定,其振蕩幅值數(shù)據(jù)如圖6(a)所示,對(duì)其進(jìn)行FFT處理。頻譜結(jié)果表明,該發(fā)動(dòng)機(jī)在工作過(guò)程中出現(xiàn)了典型的非線性燃燒不穩(wěn)定現(xiàn)象,存在多階模態(tài),且前8階模態(tài)幅值明顯,通過(guò)計(jì)算分析,分解前8階模態(tài)重構(gòu)該信號(hào)的精度能達(dá)到95%,于是這里使用SPSO-VMD方法對(duì)其振蕩數(shù)據(jù)進(jìn)行前8階模態(tài)分解,所用計(jì)算資源CPU型號(hào)為intel Xeon(R)Platinum 8259CL、48核,迭代在第9步時(shí),適應(yīng)度函數(shù)便達(dá)到了7.763 5×10-5,通過(guò)上一節(jié)的誤差分析,可以認(rèn)為該適應(yīng)度值表明分解已經(jīng)達(dá)到了較好的效果,對(duì)應(yīng)用時(shí)約為1880 s。

    (a)Oscillation signal

    8階模態(tài)的分解結(jié)果如圖7(a)~(h)所示,各分解模態(tài)的FFT結(jié)果表明,各階模態(tài)被準(zhǔn)確分離。

    (a)Mode 1 and its FFT result (b)Mode 2 and its FFT result

    針對(duì)于類(lèi)似第6階模態(tài),頻譜圖上出現(xiàn)兩個(gè)臨近頻率峰值的情況,這里作一定的分析:對(duì)第6階模態(tài)進(jìn)行STFT處理,所得結(jié)果如圖8所示,結(jié)果表明第6階模態(tài)在發(fā)動(dòng)機(jī)發(fā)生燃燒不穩(wěn)定的時(shí)間內(nèi)出現(xiàn)了頻率的變化,頻譜圖中兩個(gè)峰值恰好對(duì)應(yīng)STFT結(jié)果中振蕩幅值最大的0.29 s和0.46 s兩個(gè)時(shí)刻下的頻率,也即這兩個(gè)頻率是第6階模態(tài)不同時(shí)刻的振蕩頻率,它們屬于一個(gè)模態(tài)。這種頻率變化的現(xiàn)象主要是由燃燒室裝藥燃燒,燃面退移,導(dǎo)致燃燒室內(nèi)聲腔構(gòu)型發(fā)生變化造成的。

    圖8 第6階模態(tài)STFT處理結(jié)果 圖9 前8階模態(tài)重構(gòu)信號(hào)

    這里不考慮9階及其高階模態(tài),將得到的前8階模態(tài)重新疊加,與原始的振蕩幅值進(jìn)行對(duì)比,如圖9所示,計(jì)算前8階模態(tài)重構(gòu)振蕩幅值的精度為95.2%。如圖10所示,將前8階模態(tài)繪制在同一張圖中,在0.16~0.4 s時(shí)間段范圍內(nèi),壓強(qiáng)信號(hào)雖然保持極限環(huán)振蕩現(xiàn)象,但是各階模態(tài)的幅值卻在不斷變化,可看到奇數(shù)階模態(tài)幅值衰減到了一個(gè)極小值,而偶數(shù)階模態(tài)出現(xiàn)了一定的增長(zhǎng),仿佛出現(xiàn)了偶數(shù)階增長(zhǎng)抑制奇數(shù)階的現(xiàn)象,據(jù)文獻(xiàn)[15]所提:“非線性能量可從一階聲模態(tài)傳遞至所有高階聲模態(tài),而從二階聲模態(tài)傳遞只能向偶數(shù)階傳遞,奇數(shù)階模態(tài)均不發(fā)生能量傳遞過(guò)程,且高階聲模態(tài)不能向低階聲模態(tài)傳遞能量。”在該發(fā)動(dòng)機(jī)振蕩信號(hào)中,由于某種原因,一階模態(tài)幅值衰減到了接近0的值,傳遞向高階聲模態(tài)的能量極少,而二階聲模態(tài)幅值較大,不斷向偶數(shù)階模態(tài)傳遞著非線性能量,導(dǎo)致偶數(shù)階維持一定幅值的振蕩,而奇數(shù)階高階模態(tài)因?yàn)檩斎氲姆蔷€性能量較少,不足以維持自身耗散的作用,導(dǎo)致幅值一直處在一個(gè)較低的值。這種現(xiàn)象恰好也提示研究人員可以從能量傳遞的角度對(duì)發(fā)動(dòng)機(jī)非線性燃燒不穩(wěn)定進(jìn)行抑制,比如切斷能量傳輸?shù)穆窂交蛘吒淖兡芰總鬏數(shù)牧康鹊?這些都有待研究。

    圖10 前8階模態(tài)奇偶數(shù)階幅值演化對(duì)比

    4 結(jié)論

    (1)SPSO-VMD方法能夠準(zhǔn)確地分解發(fā)動(dòng)機(jī)非線性燃燒不穩(wěn)定典型特征信號(hào)的各階模態(tài),獲得各模態(tài)的頻率和幅值,可應(yīng)用于真實(shí)發(fā)動(dòng)機(jī)燃燒不穩(wěn)定振蕩信號(hào)的處理,可為研究非線性不穩(wěn)定能量傳遞演化機(jī)制提供數(shù)據(jù)支撐。

    (2)SPSO-VMD方法改進(jìn)了PSO-VMD方法陷入局部最優(yōu)解的缺陷,能夠準(zhǔn)確地找到全局最優(yōu)解,分解模態(tài)誤差較小,分解效率較高。

    (3)SPSO-VMD方法自適應(yīng)強(qiáng),不依賴主觀經(jīng)驗(yàn),集成性好,使用簡(jiǎn)單方便,可適用于工業(yè)部門(mén)對(duì)發(fā)動(dòng)機(jī)燃燒不穩(wěn)定問(wèn)題的定位和分析。

    猜你喜歡
    適應(yīng)度幅值模態(tài)
    改進(jìn)的自適應(yīng)復(fù)制、交叉和突變遺傳算法
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    基于空調(diào)導(dǎo)風(fēng)板成型工藝的Kriging模型適應(yīng)度研究
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    由單個(gè)模態(tài)構(gòu)造對(duì)稱簡(jiǎn)支梁的抗彎剛度
    少數(shù)民族大學(xué)生文化適應(yīng)度調(diào)查
    日本三级黄在线观看| 国产三级中文精品| 亚洲无线观看免费| 婷婷亚洲欧美| 亚洲色图av天堂| 一区二区三区激情视频| 一个人看的www免费观看视频| 一区二区三区国产精品乱码| 国产精品野战在线观看| 99精品在免费线老司机午夜| 99久久精品一区二区三区| 真人做人爱边吃奶动态| 国产精品野战在线观看| 亚洲一区二区三区不卡视频| 亚洲色图av天堂| 在线天堂最新版资源| 999久久久精品免费观看国产| 亚洲成人久久性| 女人高潮潮喷娇喘18禁视频| 亚洲精品一卡2卡三卡4卡5卡| 日本黄大片高清| 日韩 欧美 亚洲 中文字幕| 乱人视频在线观看| e午夜精品久久久久久久| 国产激情偷乱视频一区二区| 日韩人妻高清精品专区| 久久久国产精品麻豆| 麻豆国产av国片精品| 最新美女视频免费是黄的| 国产乱人伦免费视频| 丁香欧美五月| 欧美激情在线99| 日本撒尿小便嘘嘘汇集6| 久久性视频一级片| 两性午夜刺激爽爽歪歪视频在线观看| a级一级毛片免费在线观看| 亚洲天堂国产精品一区在线| 亚洲第一电影网av| 欧美大码av| 一卡2卡三卡四卡精品乱码亚洲| 久久午夜亚洲精品久久| 国产精品久久视频播放| 天堂√8在线中文| 搡女人真爽免费视频火全软件 | 久久精品夜夜夜夜夜久久蜜豆| 久久久久久人人人人人| 日本成人三级电影网站| 国产激情欧美一区二区| 中文字幕人妻熟人妻熟丝袜美 | 在线观看免费视频日本深夜| 搡老熟女国产l中国老女人| 一级毛片高清免费大全| 成人精品一区二区免费| 国产伦在线观看视频一区| 日韩精品中文字幕看吧| 亚洲av二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 无人区码免费观看不卡| 免费搜索国产男女视频| 免费观看人在逋| 欧美黄色片欧美黄色片| 天天一区二区日本电影三级| 免费观看精品视频网站| 在线看三级毛片| 十八禁网站免费在线| 成年女人看的毛片在线观看| 亚洲av不卡在线观看| 最近最新中文字幕大全免费视频| 嫩草影院精品99| 国产一区二区亚洲精品在线观看| 亚洲国产欧美人成| 丰满的人妻完整版| 男女下面进入的视频免费午夜| 久久久久亚洲av毛片大全| av在线蜜桃| h日本视频在线播放| 亚洲在线观看片| 日韩欧美在线二视频| 精品欧美国产一区二区三| 久久久国产精品麻豆| 亚洲精品成人久久久久久| 在线天堂最新版资源| 国产高清videossex| 一级作爱视频免费观看| 草草在线视频免费看| 91在线观看av| 成年女人毛片免费观看观看9| 色播亚洲综合网| 中文字幕熟女人妻在线| 亚洲美女视频黄频| 999久久久精品免费观看国产| 国产综合懂色| 美女cb高潮喷水在线观看| 一个人免费在线观看的高清视频| 国产午夜精品久久久久久一区二区三区 | www.999成人在线观看| 亚洲av五月六月丁香网| 99热精品在线国产| 91麻豆av在线| 久久久色成人| 噜噜噜噜噜久久久久久91| 亚洲国产高清在线一区二区三| 在线播放无遮挡| 欧美中文综合在线视频| 天天躁日日操中文字幕| 国产高清激情床上av| 亚洲最大成人中文| 成人18禁在线播放| 国产成人av教育| 亚洲av第一区精品v没综合| 97人妻精品一区二区三区麻豆| 夜夜躁狠狠躁天天躁| 麻豆国产av国片精品| 午夜福利在线观看吧| 国产精品精品国产色婷婷| 国产精品,欧美在线| 国产午夜福利久久久久久| 久久精品91蜜桃| 老司机深夜福利视频在线观看| a级一级毛片免费在线观看| 欧美成狂野欧美在线观看| 免费电影在线观看免费观看| 亚洲 国产 在线| 欧美区成人在线视频| 国产极品精品免费视频能看的| 国产三级中文精品| 村上凉子中文字幕在线| 日本与韩国留学比较| 国产欧美日韩精品一区二区| 亚洲午夜理论影院| 看免费av毛片| h日本视频在线播放| 亚洲成人久久性| 国产真实伦视频高清在线观看 | 国产麻豆成人av免费视频| 亚洲欧美日韩高清在线视频| 无限看片的www在线观看| 老汉色∧v一级毛片| 国产精品免费一区二区三区在线| 舔av片在线| 中亚洲国语对白在线视频| 午夜精品久久久久久毛片777| www国产在线视频色| 看片在线看免费视频| 一级黄色大片毛片| 午夜福利18| 极品教师在线免费播放| 国产高清有码在线观看视频| 色精品久久人妻99蜜桃| 一级黄片播放器| 亚洲精品久久国产高清桃花| 免费人成视频x8x8入口观看| 天天添夜夜摸| 国产亚洲av嫩草精品影院| 国产精品自产拍在线观看55亚洲| 中文字幕人妻丝袜一区二区| 99在线人妻在线中文字幕| 18禁美女被吸乳视频| 一本久久中文字幕| eeuss影院久久| 香蕉av资源在线| 欧美激情在线99| 国内精品久久久久久久电影| 一本精品99久久精品77| 黄色丝袜av网址大全| 国产伦一二天堂av在线观看| 国产精品亚洲一级av第二区| 日本五十路高清| 亚洲av美国av| 高清日韩中文字幕在线| 岛国在线免费视频观看| 欧美丝袜亚洲另类 | 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 欧美性猛交╳xxx乱大交人| 性色avwww在线观看| av黄色大香蕉| 少妇熟女aⅴ在线视频| 精品不卡国产一区二区三区| 成人18禁在线播放| 欧美+亚洲+日韩+国产| 午夜亚洲福利在线播放| 中文在线观看免费www的网站| 国产高清视频在线播放一区| 欧美一区二区亚洲| 亚洲人成网站高清观看| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 亚洲在线自拍视频| 国产真人三级小视频在线观看| 黄片大片在线免费观看| 国产伦一二天堂av在线观看| 国产一区在线观看成人免费| 熟女人妻精品中文字幕| 国产乱人伦免费视频| 国产精品一区二区三区四区久久| 老司机午夜十八禁免费视频| 成人特级黄色片久久久久久久| 久久精品综合一区二区三区| 色视频www国产| 亚洲va日本ⅴa欧美va伊人久久| 一a级毛片在线观看| 午夜福利在线观看吧| 脱女人内裤的视频| 亚洲 国产 在线| 狠狠狠狠99中文字幕| 国产色爽女视频免费观看| 麻豆久久精品国产亚洲av| 日韩欧美免费精品| 成人18禁在线播放| 欧美一区二区精品小视频在线| 免费大片18禁| 99久久无色码亚洲精品果冻| 校园春色视频在线观看| 亚洲精品在线观看二区| 久久九九热精品免费| 麻豆一二三区av精品| 热99在线观看视频| 午夜激情福利司机影院| 女同久久另类99精品国产91| 听说在线观看完整版免费高清| 黄色丝袜av网址大全| 国产爱豆传媒在线观看| 国模一区二区三区四区视频| 中文资源天堂在线| 国产精品一及| 精品无人区乱码1区二区| 手机成人av网站| 男女床上黄色一级片免费看| 国产成人欧美在线观看| 不卡一级毛片| 欧美成人a在线观看| 免费观看精品视频网站| 男女做爰动态图高潮gif福利片| 久久性视频一级片| 一个人看视频在线观看www免费 | 午夜久久久久精精品| 好男人电影高清在线观看| 亚洲av熟女| 中文字幕久久专区| 搡老岳熟女国产| 青草久久国产| av天堂中文字幕网| 国产黄色小视频在线观看| 高潮久久久久久久久久久不卡| 亚洲欧美日韩无卡精品| 桃红色精品国产亚洲av| 欧美日韩精品网址| 淫秽高清视频在线观看| 国产色婷婷99| 精品国产美女av久久久久小说| 国产成人av教育| 欧美日韩瑟瑟在线播放| 亚洲av日韩精品久久久久久密| 国产高清videossex| 久久欧美精品欧美久久欧美| 久久性视频一级片| 国产精品99久久99久久久不卡| 亚洲内射少妇av| 精品久久久久久成人av| 亚洲成人久久性| 91久久精品国产一区二区成人 | 一进一出好大好爽视频| 性欧美人与动物交配| 国产老妇女一区| 香蕉丝袜av| 日本一本二区三区精品| 91av网一区二区| 一本一本综合久久| 女同久久另类99精品国产91| 99久久无色码亚洲精品果冻| 内射极品少妇av片p| 九色成人免费人妻av| 色av中文字幕| 97人妻精品一区二区三区麻豆| 99在线视频只有这里精品首页| 激情在线观看视频在线高清| 少妇裸体淫交视频免费看高清| 久久久久性生活片| 亚洲成人免费电影在线观看| 欧美zozozo另类| 啦啦啦韩国在线观看视频| av欧美777| 成人无遮挡网站| 天堂影院成人在线观看| 日本与韩国留学比较| 亚洲av熟女| 国产久久久一区二区三区| 国产三级在线视频| 一进一出抽搐动态| 国产一区二区三区视频了| 91在线精品国自产拍蜜月 | а√天堂www在线а√下载| 成人av一区二区三区在线看| 精品免费久久久久久久清纯| 91九色精品人成在线观看| 欧美黄色片欧美黄色片| 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 一级黄片播放器| 午夜福利在线观看免费完整高清在 | 露出奶头的视频| 一本精品99久久精品77| 国产精品av视频在线免费观看| 在线看三级毛片| 国内精品美女久久久久久| 久久精品91无色码中文字幕| 一个人看视频在线观看www免费 | 亚洲人成网站在线播放欧美日韩| 女人十人毛片免费观看3o分钟| 真人做人爱边吃奶动态| 哪里可以看免费的av片| 99久久成人亚洲精品观看| 亚洲成人中文字幕在线播放| 非洲黑人性xxxx精品又粗又长| 日本黄大片高清| 天天躁日日操中文字幕| 日韩精品中文字幕看吧| 岛国在线免费视频观看| 国产探花极品一区二区| 偷拍熟女少妇极品色| 国产精品电影一区二区三区| 哪里可以看免费的av片| 欧美一级a爱片免费观看看| 国产欧美日韩一区二区三| 国产精品亚洲av一区麻豆| 亚洲精品亚洲一区二区| 亚洲成人久久爱视频| 国产极品精品免费视频能看的| 日本 欧美在线| 欧美日韩国产亚洲二区| 免费大片18禁| 搡老熟女国产l中国老女人| 亚洲av五月六月丁香网| 搡老妇女老女人老熟妇| 国产精品乱码一区二三区的特点| 午夜福利免费观看在线| 午夜日韩欧美国产| 九九热线精品视视频播放| 看免费av毛片| 国产欧美日韩精品亚洲av| 中文字幕人妻熟人妻熟丝袜美 | 中文字幕久久专区| 搡老岳熟女国产| 国产精品精品国产色婷婷| 中文字幕人妻熟人妻熟丝袜美 | bbb黄色大片| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 久久九九热精品免费| 中文字幕久久专区| 嫩草影视91久久| 国产伦人伦偷精品视频| 国产高清videossex| 一区二区三区激情视频| 久久久久久久午夜电影| 欧美色欧美亚洲另类二区| 国产熟女xx| 国产成人系列免费观看| 免费看十八禁软件| 中文字幕高清在线视频| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 久久久久久久久中文| 久久精品夜夜夜夜夜久久蜜豆| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 一级作爱视频免费观看| 嫩草影视91久久| 亚洲人成电影免费在线| 久久久久久人人人人人| 可以在线观看的亚洲视频| ponron亚洲| 欧美乱妇无乱码| 久久香蕉精品热| 18美女黄网站色大片免费观看| 亚洲av五月六月丁香网| 欧美成人a在线观看| 婷婷精品国产亚洲av| 久久6这里有精品| 成人高潮视频无遮挡免费网站| 18禁在线播放成人免费| 亚洲午夜理论影院| 制服丝袜大香蕉在线| 国产免费男女视频| 51午夜福利影视在线观看| 欧美绝顶高潮抽搐喷水| 欧美成人性av电影在线观看| 国语自产精品视频在线第100页| 欧美zozozo另类| 老熟妇乱子伦视频在线观看| 久久亚洲真实| 精品不卡国产一区二区三区| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 午夜精品在线福利| 亚洲专区中文字幕在线| 无人区码免费观看不卡| 村上凉子中文字幕在线| 亚洲欧美精品综合久久99| 亚洲成人中文字幕在线播放| 亚洲avbb在线观看| 九九久久精品国产亚洲av麻豆| 亚洲不卡免费看| 男人舔女人下体高潮全视频| 国产在视频线在精品| 国产精品久久久人人做人人爽| 精品人妻一区二区三区麻豆 | 免费看日本二区| 亚洲国产精品sss在线观看| 久久久国产成人精品二区| 97人妻精品一区二区三区麻豆| 久久草成人影院| 一a级毛片在线观看| 熟女少妇亚洲综合色aaa.| 国产高清激情床上av| 日本在线视频免费播放| 亚洲av免费在线观看| 精品一区二区三区av网在线观看| 成年女人永久免费观看视频| 中文字幕高清在线视频| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看| 黄片大片在线免费观看| 国产高潮美女av| 可以在线观看的亚洲视频| 大型黄色视频在线免费观看| e午夜精品久久久久久久| 午夜福利18| 成熟少妇高潮喷水视频| 久久亚洲精品不卡| 国产精品一区二区免费欧美| 长腿黑丝高跟| 国产熟女xx| 搡老妇女老女人老熟妇| 麻豆成人av在线观看| 亚洲五月婷婷丁香| 中文字幕av成人在线电影| 18禁黄网站禁片午夜丰满| 亚洲精品一区av在线观看| 欧美中文日本在线观看视频| 久久99热这里只有精品18| 高清在线国产一区| 99精品欧美一区二区三区四区| 中文字幕久久专区| 欧美性猛交╳xxx乱大交人| 国产免费男女视频| 一边摸一边抽搐一进一小说| 国产99白浆流出| bbb黄色大片| 又粗又爽又猛毛片免费看| 欧美色欧美亚洲另类二区| 观看美女的网站| 精品久久久久久久人妻蜜臀av| 久久久久久久久大av| 欧美国产日韩亚洲一区| 亚洲18禁久久av| 久久久久久久久中文| 免费看美女性在线毛片视频| 欧美不卡视频在线免费观看| 男人的好看免费观看在线视频| 国产一级毛片七仙女欲春2| 1024手机看黄色片| 久久欧美精品欧美久久欧美| 99久国产av精品| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频| 欧美性猛交╳xxx乱大交人| 日本熟妇午夜| 禁无遮挡网站| 97碰自拍视频| 亚洲男人的天堂狠狠| 国产免费av片在线观看野外av| 九色国产91popny在线| 免费高清视频大片| 国产亚洲精品av在线| 69人妻影院| 国产亚洲精品久久久久久毛片| 欧美色视频一区免费| ponron亚洲| 国产97色在线日韩免费| 国产精品久久久久久精品电影| xxx96com| 最近视频中文字幕2019在线8| 欧美+日韩+精品| 久久精品91无色码中文字幕| 亚洲成人久久性| 丁香欧美五月| 一级黄片播放器| 亚洲精品成人久久久久久| 看片在线看免费视频| 精品99又大又爽又粗少妇毛片 | 久久精品91蜜桃| 少妇裸体淫交视频免费看高清| 一级毛片高清免费大全| 十八禁网站免费在线| 观看免费一级毛片| 国产不卡一卡二| 熟女人妻精品中文字幕| 动漫黄色视频在线观看| 亚洲精华国产精华精| 一个人看视频在线观看www免费 | 国产伦精品一区二区三区四那| 岛国视频午夜一区免费看| 最新美女视频免费是黄的| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 少妇丰满av| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 日韩有码中文字幕| 我的老师免费观看完整版| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| 青草久久国产| 夜夜看夜夜爽夜夜摸| 日韩高清综合在线| 网址你懂的国产日韩在线| 老汉色av国产亚洲站长工具| 最近最新中文字幕大全电影3| 久久这里只有精品中国| 国产精品久久久久久久久免 | 波多野结衣高清作品| 国产亚洲精品av在线| 一区二区三区免费毛片| 91字幕亚洲| 午夜福利成人在线免费观看| 波多野结衣高清作品| 久久精品国产亚洲av涩爱 | 久久久久精品国产欧美久久久| 一二三四社区在线视频社区8| 欧美绝顶高潮抽搐喷水| 欧美乱妇无乱码| 亚洲真实伦在线观看| 黄色丝袜av网址大全| 久久久久久国产a免费观看| 老熟妇乱子伦视频在线观看| 精品一区二区三区视频在线 | 亚洲在线自拍视频| 免费看美女性在线毛片视频| 亚洲人成网站在线播放欧美日韩| 两人在一起打扑克的视频| 国产主播在线观看一区二区| 麻豆一二三区av精品| 欧美在线一区亚洲| 悠悠久久av| 国产精品亚洲av一区麻豆| 亚洲精品色激情综合| 国产伦一二天堂av在线观看| 免费观看的影片在线观看| 国产三级在线视频| 国产在线精品亚洲第一网站| 一夜夜www| 亚洲中文字幕一区二区三区有码在线看| 女人被狂操c到高潮| 黄片小视频在线播放| 欧美日韩国产亚洲二区| 男人的好看免费观看在线视频| 村上凉子中文字幕在线| 淫妇啪啪啪对白视频| 丰满人妻熟妇乱又伦精品不卡| 夜夜躁狠狠躁天天躁| 香蕉丝袜av| 国产极品精品免费视频能看的| bbb黄色大片| 草草在线视频免费看| 人人妻人人看人人澡| 一边摸一边抽搐一进一小说| 久久久久亚洲av毛片大全| 一本精品99久久精品77| 香蕉丝袜av| 不卡一级毛片| 校园春色视频在线观看| 久久久久久久午夜电影| 国产老妇女一区| 女人被狂操c到高潮| 五月玫瑰六月丁香| 国内精品一区二区在线观看| 一级作爱视频免费观看| 国产高清激情床上av| 18+在线观看网站| 黄片大片在线免费观看| 美女大奶头视频| 国产色爽女视频免费观看| 国产成年人精品一区二区| 麻豆一二三区av精品| 99久久精品国产亚洲精品| 日本成人三级电影网站| 欧美色欧美亚洲另类二区| 99热只有精品国产| 免费看光身美女| 熟女人妻精品中文字幕| 精品国产美女av久久久久小说| 欧美日韩国产亚洲二区| 真人一进一出gif抽搐免费| 欧美在线黄色| 一级毛片高清免费大全| 日韩高清综合在线| 亚洲国产精品999在线| 美女被艹到高潮喷水动态| 久久精品91无色码中文字幕| 午夜福利18| 啦啦啦韩国在线观看视频| 亚洲精华国产精华精| 亚洲精品一区av在线观看| 欧美一级a爱片免费观看看| 99久久无色码亚洲精品果冻| 精品人妻1区二区| 午夜精品一区二区三区免费看| 国产高清有码在线观看视频| 成人国产一区最新在线观看| 级片在线观看| 91在线精品国自产拍蜜月 | 热99re8久久精品国产| 国产精品,欧美在线|