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

    基于EMD的某火炮身管模態(tài)參數(shù)識別

    2017-02-09 02:00:28王永剛
    兵器裝備工程學報 2017年1期
    關鍵詞:身管火炮阻尼

    王永剛,尹 強

    (南京理工大學,南京 210094)

    【裝備理論與裝備技術】

    基于EMD的某火炮身管模態(tài)參數(shù)識別

    王永剛,尹 強

    (南京理工大學,南京 210094)

    現(xiàn)階段火炮模態(tài)參數(shù)識別都是以快速傅里葉變換(FFT)為基礎,信號必須平穩(wěn)且是嚴格的周期信號;通常情況火炮工作狀態(tài)下的信號是非平穩(wěn)信號;經(jīng)驗模式分解(EMD)方法是一種新的非平穩(wěn)的信號處理技術,EMD方法在處理非平穩(wěn)信號時相對其他方法有相當大的優(yōu)勢;首先對某火炮身管在有限元分析軟件中建模,分別在有約束和自由狀態(tài)下對身管進行模態(tài)識別;然后用EMD方法對相關數(shù)據(jù)進行處理,得到身管的模態(tài)參數(shù)(頻率、阻尼比、剛度矩陣、阻尼矩陣、質(zhì)量矩陣),并與有限元分析方法得到的結(jié)果進行對比、分析;數(shù)值仿真證明了EMD方法在火炮領域的可行性和有效性。

    火炮;身管;EMD;模態(tài)參數(shù)識別

    火炮作為戰(zhàn)場上的常規(guī)武器,具有火力密集、反應迅速、抗干擾能力強等特點。多年來,我國的火炮機械設計,把功能要求和靜強度作為主要準則,因此常常會出現(xiàn)保守設計的情況,導致火炮的設計質(zhì)量大而且成本比較高[1]。隨著火炮設計技術的不斷發(fā)展,在設計時不但要考慮靜態(tài)特性還要考慮動態(tài)特性。這種火炮設計的需求刺激了模態(tài)分析技術的迅速發(fā)展,加上近些年來計算機技術的高速發(fā)展和快速傅里葉動態(tài)分析技術的發(fā)展,使模態(tài)分析技術成功應用于火炮的結(jié)構設計和性能分析?,F(xiàn)階段的模態(tài)參數(shù)識別方法很多,應用于地震工程和建筑工程上的模態(tài)參數(shù)識別方法主要有分量分析法、加權最小二乘迭代法、有理分式多項式法、正交多項式法、時域復指數(shù)法、Ibrahim時域法(ITD)、節(jié)約時域法(STD)?;趏utput-only的模態(tài)參數(shù)方法有峰值法、頻域分解法、PolyMAX法、隨機子空間法、隨機減量技術(RDT)、自然激勵技術(NExT)、特征系統(tǒng)實驗算法等[2]。

    在火炮方面的應用,現(xiàn)階段識別模態(tài)參數(shù)處理數(shù)據(jù)的方法有蘭索斯法,子空間法,降階法等。應用比較廣泛的是傳統(tǒng)的實驗模態(tài)分析,這種方法通過測量結(jié)構中各個部位在特定激勵下的振動輸出響應,用模態(tài)分析理論及相應的計算軟件分析處理所記錄的激勵信號和響應信號,獲得頻響函數(shù),完成模態(tài)參數(shù)的識別[3]。現(xiàn)在大部分模態(tài)識別方法都是基于FFT理論來完成的。但是FFT分析方法并非適用于所有類型的信號分析,F(xiàn)ourier分析有嚴格的限制條件[4]:信號必須是平穩(wěn)信號而且是嚴格的周期信號;被分析的系統(tǒng)必須是線性系統(tǒng)。否則Fourier分析將失去物理意義。工作狀態(tài)下的火炮,其響應信號多數(shù)情況是非線性非穩(wěn)定信號,基于FFT的模態(tài)分析方法將失去物理意義[5]。

    模態(tài)參數(shù)的識別對火炮結(jié)構設計和性能分析有著重要的意義。根據(jù)火炮的機構設計,要求到工廠加工規(guī)定性能指標的零部件,由于受到加工工序和儀器設備誤差等影響,在加工完成后,得到的零部件往往有一定偏差,這就需要對火炮進行模態(tài)參數(shù)識別,測定加工后的零部件的實際模態(tài)參數(shù),如果沒有達到設計要求,則需要尋找原因,找出哪一環(huán)節(jié)影響了加工結(jié)果。火炮的模態(tài)參數(shù)識別不僅可以加快火炮的設計進程而且可以減小結(jié)構設計成本?;鹋诘慕Y(jié)構參數(shù)還會影響到火炮的戰(zhàn)術技術指標,如射擊準確度,射擊密集度等。所以火炮的模態(tài)參數(shù)識別對火炮的結(jié)構設計和技術分析至關重要。

    1 經(jīng)驗模式分解方法(EMD)

    1998年美籍華人N.E.Huang等[6]提出一種信號處理方法經(jīng)驗模式分解法(Empirical Mode Decompostion),簡稱EMD方法。該方法是將信號進行平穩(wěn)化處理,將信號分解成具有不同尺度波動,不同特征的數(shù)據(jù)序列。每個序列成為一個固有模態(tài)函數(shù)(Intrinsic Mode Function),簡稱IMF。對每個IMF進行Hibert變換,可以得到信號的瞬時振幅和瞬時頻率。與FFT相比,希爾伯特變換(HHT)可以處理非平穩(wěn)和瞬態(tài)問題。HHT得到的每階IMF的振幅和頻率是隨時間變化的,消除了為反映信號的非平穩(wěn)性引入的多余且無物理意義的簡諧波,使信號分析更加靈活方便[7]。

    EMD方法具體過程:

    (1)對任意給定信號x(t),首先確定出x(t)上的所有極值點,用三次樣條曲線連接所有極大值點形成上包絡線和下包絡線。數(shù)據(jù)x(t)與上下包絡線的均值m1的差記為h1,則h1=x(t)-m1。將h1視為新的x(t),重復上述步驟,直到h1滿足IMF的兩個條件(整個事件歷程內(nèi),穿越零點的次數(shù)與極值點數(shù)相等或最多相差1;信號上任意一點,有局部極大值定義的上包絡線和由局部極小值定義的下包絡線的均值為0,即信號關于時間軸局部對稱)時,則其成為原始信號篩選出的第一階IMF,記為C1。

    (2)將C1從x(t)中分離出來得到一個去掉高頻分量的差值信號r1,把r1作為新信號,重復(1)的篩選步驟,直到第n階的殘余信號成為單調(diào)函數(shù)不能再分IMF為止[8]。

    2 希爾伯特變換

    EMD方法基于信號的局部特征時間尺度,將信號自適應地分解為若干個IMF分量之和,對EMD得到的每個特征模態(tài)函數(shù)作希爾伯特變換(HHT)[9]。x(t)為測量信號,則

    (1)

    v(t)的解析信號Y(t)表示為

    (2)

    ω(t)=dθ(t)/dt

    (3)

    對于式(2)中幅值A(t)及相位角θ(t)可進一步表示為

    A(t)=A0e-ζω0t

    θ(t)=ωdt+φ0

    (4)

    對幅值A(t)及相位角θ(t)引入對數(shù)及微分算子,則

    (5)

    3 模態(tài)參數(shù)的識別

    對θ(t)-t圖,用線性最小均方擬合過程進行直線擬合,直線斜率即為ωd。對lnA(t)-t圖,同樣用線性最小均方擬合過程進行直線擬合,直線斜率即為-ζjω0。對于ζ比較小時,即使lnA(t)在擬合直線附近波動,也不會影響識別結(jié)果[10]。

    對于身管的自然頻率和阻尼比只需任意位置的響應信號即可求出,然而要識別模態(tài)質(zhì)量、剛度、阻尼必須得到所有自由度的測量值。

    (6)

    (7)

    由式(6),式(7)可以確定φpj相對于φqj的符號。 這樣,在j階復模態(tài)向量Φj中的所有元素相對于某一特定元素的絕對值和相位角都可以確定。復數(shù)值aj可以導出:

    (8)

    bj=-λjaj

    (9)

    然后利用復模態(tài)的正交特性可以得到質(zhì)量矩陣,剛度矩陣:

    (10)

    (11)

    盡管由式(11)可以確定阻尼矩陣C,但由于估算誤差可能不能滿足運動方程,要從運動方程來估算C??梢缘玫剑?/p>

    K[ΦΦ*] = 0

    (12)

    MΦΛ2+CΦΛ+KΦ=0

    (13)

    C=-{MΦΛ2+KΦ}Λ-1Φ-1

    (14)

    由于對復模態(tài)振型和λj的估算誤差,由等式(11)和(14)估算的剛度和阻尼矩陣可能不是對稱陣。因此,用平均過程處理使得K和C對稱:

    (15)

    由此便可計算出火炮身管的質(zhì)量矩陣,剛度矩陣和阻尼矩陣。

    4 固定約束下身管模態(tài)參數(shù)識別

    算例1:對某火炮身管用有限元軟件abaqus進行建模,然后對身管進行模態(tài)參數(shù)識別。身管模型如圖1所示。

    身管的末端加固定約束。身管的長度為8.060 6 m彈性模量為210 kPa,泊松比為0.3,密度為7.85×10-9t/mm3。在Ls處沿Z方向施加54 N的力,作用時間為0.002 s,分別在L1,L2,L3,L4位置處輸出加速度信號,信號輸出位置之間距離為1.652 0 m。

    然后用有限元法對火炮的模態(tài)參數(shù)進行識別,得到模態(tài)參數(shù)數(shù)據(jù),然后對模型產(chǎn)生的加速度信號用EMD方法進行參數(shù)識別。如圖2所示為L1位置處加速度響應信號。

    圖2 有約束身管L1位置處加速度響應信號

    對實驗得到的位置L1加速度信號進行EMD分解,首先對信號進行FFT變換。如圖3所示。

    圖3 L1位置處信號FFT變換

    如圖4所示身管一階頻率大概為20 Hz,對加速度信號進行15~30 Hz的帶通濾波,得到濾波后的響應信號x1(t)。

    圖4 L1位置處信號濾波后信號x1(t)

    然后對x1(t)進行EMD分解。分解得到IMF如圖5所示。IMF1近似于一階模態(tài)響應信號,然后對IMF1進行HHT變換得到相位角圖和幅值對數(shù)圖。

    圖5 x1(t)進行EMD分解得到的IMF曲線

    如圖6、圖7所示,可以算得ωd=153.18,-ζω0=-2.76。求得f=24.38 Hz,ζ=0.11。由此身管的第一階固有頻率ω0=157.65,f0=25.09。

    圖6 相位角曲線

    圖7 幅值對數(shù)曲線

    身管的前四階固有頻率和阻尼比如表1。

    表1 理論值與EMD識別值

    身管的質(zhì)量,剛度和阻尼矩陣如表2。從數(shù)據(jù)上來看,在固定約束條件下,EMD方法對身管的模態(tài)參數(shù)識別與理論值偏差不大,可用百分比給出各參數(shù)的誤差數(shù)值的最大值。

    5 自由狀態(tài)下身管的模態(tài)參數(shù)識別

    算例2:對身管在自由狀態(tài)下進行模態(tài)參數(shù)識別。方法步驟同固定約束狀態(tài)情況的識別方法。

    自由狀態(tài)下身管的L1位置處響應如圖8所示。

    圖8 自由狀態(tài)下L1位置處響應信號

    然后對自由狀態(tài)下身管的L1位置處響應信號進行FFT變換,估算身管固有頻率值。變換結(jié)果如下:

    如圖9可知,身管的一階固有頻率大概在20 Hz左右,然后對身管進行15~25 Hz的帶通濾波。濾波之后的響應信號x2(t)如圖10所示。

    圖9 自由狀態(tài)下L1位置信號響應FFT結(jié)果

    理論值(有限元方法)識別值(EMD方法)質(zhì)量矩陣/kg915.670000470.250000411.550000392.14901.231754498.45101385399.986386375.98剛度矩陣/kN2907290.66-1116645.5200-1116645.5220831186.55-966541.0300-966541.031856083.03-88954200-8895428895422907123.54-1116625.47250-34-1116625.4720832122.65-1011254.037425-969987.731856978.25-889114.21-105-14-889114.21889114.21阻尼矩陣/(N·s·m-1)540017810-1575009492017810477053430-389870-157500534302310-12134094920-389870-1213406940432014250-1260007594014250382042750-311900126000427501850970707594031900970705550

    圖10 自由狀態(tài)身管L1位置濾波后響應信號x2(t)

    然后對濾波后的響應信號x2(t)進行EMD分解,得到固有模式函數(shù)IMFS,如圖11所示。

    圖11 固有模式函數(shù)IMFS

    對第一個IMF進行HHT變換,求身管的阻尼頻率ωd1和-ζω0的值。

    如圖12、圖13所示,可以算得到ωd=131.19,-ζω0=-2.52,f=20.88 Hz,ζ=0.12。同理計算自由狀態(tài)下身管的其他階次固有頻率和阻尼比如表3。自由狀態(tài)下身管的質(zhì)量矩陣、剛度矩陣阻尼矩陣如表4所示。

    圖12 相位角曲線

    圖13 幅值對數(shù)曲線

    模態(tài)頻率/Hz理論值(有限元方法)識別值(EMD方法)誤差/%阻尼比理論值(有限元方法)識別值(EMD方法)誤差/%一階20.9020.014.40.120.138.3二階58.3458.380.060.180.175.6三階112.65113.510.760.230.224.3四階186.01185.770.130.270.283.7

    表4 自由狀態(tài)下身管的有限元法和EMD法的結(jié)果數(shù)據(jù)

    6 結(jié)論

    本文通過在有限元軟件中建立火炮身管模型,算出了在一端固定約束和自由狀態(tài)兩種情況下火炮身管的模態(tài)參數(shù),用EMD的方法對這兩種情況下的火炮身管進行模態(tài)參數(shù)識別,并分別與有限元法的結(jié)果進行對比。說明EMD識別結(jié)果可靠,EMD相對于傳統(tǒng)的模態(tài)實驗方法數(shù)據(jù)處理方便,結(jié)果可靠,而且適用于工作狀態(tài)下火炮產(chǎn)生的非線性信號。對火炮模態(tài)參數(shù)識別,對火炮的結(jié)構設計,性能指標的分析有著重要的意義。

    [1] 余成寶.火炮系統(tǒng)模態(tài)測試與分析[D].南京:南京理工大學,2007.

    [2] 郭川睿.地震激勵下結(jié)構模態(tài)參數(shù)識別及振動臺試驗驗證[D].長沙:中南大學,2014.

    [3] 任春,張繼承,羅奇峰.HHT方法在結(jié)構模態(tài)參數(shù)時域識別中的應用[J].長江大學學報,2008,5(4):115-118.

    [4] 莫平杰.基于HHT方法的機械系統(tǒng)模態(tài)參數(shù)識別[D].杭州:浙江大學,2011.

    [5] 王慧.HHT方法及其若干應用研究[D].合肥:合肥工業(yè)大學,2009.

    [6] 孫兆偉.基于現(xiàn)代信號處理的結(jié)構模態(tài)參數(shù)識別與損傷識別研究[D].北京:北京郵電大學,2012.

    [7] 任宜春.基于小波分析的結(jié)構參數(shù)識別方法研究[D].長沙:湖南大學,2007.

    [8] 楊永鋒,吳亞鋒.經(jīng)驗模態(tài)分解在振動分析中的應用[M].北京:國防工業(yè)出版社,2013.

    [9] YANG J N,LEI Y,LIN S,et al.Hilbert-Huang Transform Based Approach For Structural Damage Detection,ASCE Journal of Engineering Mechanics[J].ASCE Journal of Engineering Mechanics,2004,130(1):85-95.

    [10]YANG J N,LEI Y,PAN S W,et al.Identification of Linear Structures Based on Hilbert-Huang Spectral Analysis[J].Journal of Earthquake Engineering and Structural Dynamics,2003,32(10):1443-1467.

    [11]王鄭力, 伍鵬宇,羅斌.某型火箭炮身管改進設計[J].兵工自動化,2016(1):93-96.

    (責任編輯周江川)

    Modal Parameters Identification of Gun Barrel Based on EMD Method

    WANG Yong-gang,YIN Qiang

    (Nanjing University of Science and Technology, Nanjing 210094, China)

    The modal parameters identification of artillery is mostly based on the Fast Fourier Transform (FFT). However, the FFT signals must be smooth and strictly periodic. The signal was usually non-stationary in the gun working condition. The EMD method had a great advantage than other methods in dealing with non-stationary signals. In this paper, firstly, we built the model of a gun barrel in the finite element analysis software. Secondly, we identified the gun modal parameters in the constraints and free stated. Thirdly, we dealt the data with the EMD method, and got the gun modal parameters (frequency, damping ratio, stiffness matrix, damping matrix, mass matrix). Finally, we compared the result of the EMD method with the result of the finite element analysis method. The numerical simulation results demonstrate the feasibility and effectiveness of the EMD method for modal parameter identification of artillery.

    gun; barrel; EMD; identify modal parameter

    2016-08-15;

    2016-09-20

    王永剛(1991—),男,碩士研究生,主要從事火炮技術研究。

    10.11809/scbgxb2017.01.017

    王永剛,尹強.基于EMD的某火炮身管模態(tài)參數(shù)識別[J].兵器裝備工程學報,2017(1):69-74.

    format:WANG Yong-gang,YIN Qiang.Modal Parameters Identification of Gun Barrel Based on EMD Method[J].Journal of Ordnance Equipment Engineering,2017(1):69-74.

    TH113

    A

    猜你喜歡
    身管火炮阻尼
    火炮身管壽命評估方法及其標準研究
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關于具有阻尼項的擴散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    火炮也減肥
    基于垂向固有振動仿真的身管口徑倍數(shù)設計研究
    基于數(shù)值仿真的四種類型身管強度比較
    機械工程師(2018年6期)2018-06-14 06:12:28
    具阻尼項的Boussinesq型方程的長時間行為
    身管兵器
    陸戰(zhàn)之神——火炮
    久久青草综合色| 午夜福利一区二区在线看| 国产精品98久久久久久宅男小说| 国产日本99.免费观看| 18禁美女被吸乳视频| www日本黄色视频网| 日本黄色视频三级网站网址| 黄片播放在线免费| 男人操女人黄网站| 丝袜美腿诱惑在线| 日韩欧美一区二区三区在线观看| 亚洲三区欧美一区| av电影中文网址| 国产人伦9x9x在线观看| 最近最新免费中文字幕在线| 亚洲av成人av| 久久久久久人人人人人| 99国产精品99久久久久| www.精华液| 在线观看午夜福利视频| 国产亚洲精品第一综合不卡| 久久久久国产精品人妻aⅴ院| 精品第一国产精品| 中文字幕精品亚洲无线码一区 | 母亲3免费完整高清在线观看| 十八禁网站免费在线| 成人三级黄色视频| 18禁观看日本| 国内少妇人妻偷人精品xxx网站 | 国产成人av教育| 国产成人av教育| 午夜精品久久久久久毛片777| 视频在线观看一区二区三区| 欧美在线一区亚洲| av视频在线观看入口| 淫秽高清视频在线观看| 久久精品国产清高在天天线| 色综合亚洲欧美另类图片| 美女大奶头视频| 大型黄色视频在线免费观看| 国产精品综合久久久久久久免费| 在线十欧美十亚洲十日本专区| 亚洲国产欧美网| 精品国内亚洲2022精品成人| 久久久久久久久免费视频了| 国产片内射在线| 欧美丝袜亚洲另类 | www国产在线视频色| 两人在一起打扑克的视频| 午夜激情福利司机影院| 99riav亚洲国产免费| 满18在线观看网站| 俄罗斯特黄特色一大片| 精品国产乱码久久久久久男人| 日本免费a在线| 黄频高清免费视频| 免费女性裸体啪啪无遮挡网站| 亚洲熟妇中文字幕五十中出| av在线播放免费不卡| 久久久国产欧美日韩av| 日本免费一区二区三区高清不卡| 精品一区二区三区视频在线观看免费| 757午夜福利合集在线观看| 真人做人爱边吃奶动态| 91老司机精品| 国内精品久久久久精免费| 免费在线观看日本一区| 亚洲欧洲精品一区二区精品久久久| 欧美成人免费av一区二区三区| 女警被强在线播放| 1024手机看黄色片| 久久精品国产亚洲av高清一级| a级毛片a级免费在线| 久久亚洲真实| 成人午夜高清在线视频 | 欧美 亚洲 国产 日韩一| 亚洲 欧美一区二区三区| 精品久久久久久久久久久久久 | 亚洲一区二区三区不卡视频| 可以在线观看毛片的网站| 免费在线观看日本一区| 91大片在线观看| e午夜精品久久久久久久| 少妇的丰满在线观看| 99热6这里只有精品| 国内精品久久久久久久电影| 欧美色欧美亚洲另类二区| 久久精品国产99精品国产亚洲性色| e午夜精品久久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 国产极品粉嫩免费观看在线| 亚洲欧美激情综合另类| 国产精华一区二区三区| 国产亚洲精品久久久久5区| 国产三级黄色录像| 精品乱码久久久久久99久播| 亚洲,欧美精品.| 他把我摸到了高潮在线观看| 国产真人三级小视频在线观看| 久久久国产欧美日韩av| 巨乳人妻的诱惑在线观看| 嫩草影视91久久| 动漫黄色视频在线观看| 少妇粗大呻吟视频| 久久精品人妻少妇| 国产又爽黄色视频| 亚洲第一av免费看| 黄色女人牲交| 午夜老司机福利片| 欧美黑人精品巨大| 日韩 欧美 亚洲 中文字幕| 精品电影一区二区在线| 看黄色毛片网站| 在线国产一区二区在线| 午夜日韩欧美国产| 久久精品国产综合久久久| 夜夜躁狠狠躁天天躁| 国产在线精品亚洲第一网站| 两个人看的免费小视频| 精品久久久久久成人av| 成人国语在线视频| 三级毛片av免费| 午夜老司机福利片| 精品无人区乱码1区二区| 高清在线国产一区| 精品久久久久久久久久久久久 | 黑人操中国人逼视频| 国产伦在线观看视频一区| 亚洲av五月六月丁香网| 亚洲五月色婷婷综合| 一级a爱视频在线免费观看| 黑人巨大精品欧美一区二区mp4| 欧美一区二区精品小视频在线| 午夜福利18| 久久狼人影院| 欧美黑人精品巨大| 亚洲国产欧洲综合997久久, | 日本在线视频免费播放| 在线av久久热| √禁漫天堂资源中文www| 久久国产精品人妻蜜桃| 国产成人精品久久二区二区免费| 亚洲男人的天堂狠狠| xxx96com| 国产1区2区3区精品| 看片在线看免费视频| 欧美不卡视频在线免费观看 | 身体一侧抽搐| 日本黄色视频三级网站网址| 黄色 视频免费看| 妹子高潮喷水视频| 欧美一区二区精品小视频在线| 午夜日韩欧美国产| 亚洲一区高清亚洲精品| 51午夜福利影视在线观看| 亚洲国产精品999在线| 少妇熟女aⅴ在线视频| 老司机福利观看| 久久精品91无色码中文字幕| 国产私拍福利视频在线观看| 一边摸一边做爽爽视频免费| 人妻丰满熟妇av一区二区三区| 亚洲成人久久爱视频| videosex国产| 欧美又色又爽又黄视频| 精品卡一卡二卡四卡免费| 久久久精品国产亚洲av高清涩受| 成在线人永久免费视频| 国产av在哪里看| 中文字幕精品免费在线观看视频| av片东京热男人的天堂| 欧美性猛交╳xxx乱大交人| 精品久久久久久,| av中文乱码字幕在线| 国产又爽黄色视频| 亚洲国产看品久久| 啦啦啦免费观看视频1| 精品国产一区二区三区四区第35| 欧美黑人精品巨大| 91麻豆精品激情在线观看国产| 一进一出好大好爽视频| 黄片小视频在线播放| 长腿黑丝高跟| videosex国产| 男女之事视频高清在线观看| 两个人视频免费观看高清| 啪啪无遮挡十八禁网站| 久9热在线精品视频| 老司机午夜十八禁免费视频| 欧美一级毛片孕妇| 久久久久久久久免费视频了| 91在线观看av| 不卡av一区二区三区| 精品少妇一区二区三区视频日本电影| 不卡一级毛片| 午夜福利18| 黄片小视频在线播放| 美女高潮喷水抽搐中文字幕| 19禁男女啪啪无遮挡网站| 黑丝袜美女国产一区| 一个人观看的视频www高清免费观看 | √禁漫天堂资源中文www| 日本三级黄在线观看| 级片在线观看| av电影中文网址| 亚洲 欧美 日韩 在线 免费| 伊人久久大香线蕉亚洲五| 午夜福利欧美成人| 中文字幕人成人乱码亚洲影| av在线播放免费不卡| 国产蜜桃级精品一区二区三区| 亚洲第一青青草原| 亚洲中文av在线| 听说在线观看完整版免费高清| 免费观看人在逋| av超薄肉色丝袜交足视频| 精品久久久久久久末码| 成年人黄色毛片网站| 成年免费大片在线观看| 日本在线视频免费播放| 麻豆一二三区av精品| 欧美日韩精品网址| 国内少妇人妻偷人精品xxx网站 | 久久久精品国产亚洲av高清涩受| 精品欧美一区二区三区在线| 午夜久久久久精精品| 亚洲av成人一区二区三| 黄色成人免费大全| 神马国产精品三级电影在线观看 | 亚洲自偷自拍图片 自拍| 欧美最黄视频在线播放免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一欧美日韩一区二区三区| 热99re8久久精品国产| 久久亚洲真实| 欧美zozozo另类| 丝袜人妻中文字幕| 又大又爽又粗| 黄频高清免费视频| 国产精品日韩av在线免费观看| 欧美三级亚洲精品| 国产伦一二天堂av在线观看| 国内精品久久久久精免费| 国内精品久久久久精免费| www日本在线高清视频| 久久精品91蜜桃| 午夜福利18| 亚洲av成人av| 国产成人精品久久二区二区91| 欧美乱妇无乱码| 久久久久国产一级毛片高清牌| svipshipincom国产片| 看片在线看免费视频| 亚洲一区中文字幕在线| 悠悠久久av| av在线播放免费不卡| 天堂√8在线中文| 啪啪无遮挡十八禁网站| 精品熟女少妇八av免费久了| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 色老头精品视频在线观看| ponron亚洲| 亚洲片人在线观看| 免费在线观看完整版高清| 免费观看人在逋| 成人手机av| 99riav亚洲国产免费| www.www免费av| 精品久久蜜臀av无| 一区二区三区国产精品乱码| 亚洲精品中文字幕在线视频| 一a级毛片在线观看| 丝袜美腿诱惑在线| 久99久视频精品免费| 久久国产精品人妻蜜桃| av免费在线观看网站| 自线自在国产av| 免费电影在线观看免费观看| 久久人妻福利社区极品人妻图片| 精品熟女少妇八av免费久了| 人人澡人人妻人| 欧美精品啪啪一区二区三区| 巨乳人妻的诱惑在线观看| 色综合站精品国产| tocl精华| 日韩欧美国产在线观看| 99久久99久久久精品蜜桃| 亚洲精品国产区一区二| 亚洲美女黄片视频| 国产精品爽爽va在线观看网站 | 国产成人av教育| 精品国产乱子伦一区二区三区| 成人国产综合亚洲| 亚洲国产欧美日韩在线播放| 久久精品夜夜夜夜夜久久蜜豆 | 久久久久亚洲av毛片大全| 男女床上黄色一级片免费看| www.www免费av| 成人三级做爰电影| 欧美 亚洲 国产 日韩一| 搡老岳熟女国产| 久久九九热精品免费| 亚洲人成电影免费在线| 可以在线观看毛片的网站| 欧美日韩黄片免| 国产精品一区二区三区四区久久 | 欧美中文日本在线观看视频| 国产成+人综合+亚洲专区| 欧美一区二区精品小视频在线| 欧美黄色片欧美黄色片| 精品久久久久久久末码| 亚洲av电影不卡..在线观看| 国产亚洲精品一区二区www| 精品电影一区二区在线| 岛国视频午夜一区免费看| 给我免费播放毛片高清在线观看| 精品久久久久久久久久免费视频| 精品少妇一区二区三区视频日本电影| 一区二区日韩欧美中文字幕| 国产成人精品无人区| 精品久久蜜臀av无| 欧美日韩瑟瑟在线播放| 午夜福利18| a级毛片a级免费在线| 很黄的视频免费| 黑人操中国人逼视频| 波多野结衣高清无吗| av有码第一页| 久9热在线精品视频| 午夜免费鲁丝| 国产亚洲精品久久久久久毛片| 国产又色又爽无遮挡免费看| 黄频高清免费视频| xxxwww97欧美| 不卡av一区二区三区| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 男女那种视频在线观看| 丝袜在线中文字幕| 午夜免费观看网址| 在线观看免费日韩欧美大片| 99精品在免费线老司机午夜| 此物有八面人人有两片| 欧美日韩乱码在线| 亚洲人成网站高清观看| 侵犯人妻中文字幕一二三四区| 岛国视频午夜一区免费看| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕日韩| 日韩视频一区二区在线观看| 午夜久久久久精精品| 手机成人av网站| 亚洲精品色激情综合| 可以在线观看的亚洲视频| 欧美色欧美亚洲另类二区| 色av中文字幕| 亚洲五月色婷婷综合| 淫秽高清视频在线观看| 国产极品粉嫩免费观看在线| 欧美中文日本在线观看视频| 亚洲欧美日韩高清在线视频| 国产亚洲精品一区二区www| 少妇 在线观看| 欧美激情 高清一区二区三区| 国产av又大| 久久精品亚洲精品国产色婷小说| 欧美黑人精品巨大| 嫁个100分男人电影在线观看| 99久久99久久久精品蜜桃| 美女午夜性视频免费| 久久精品成人免费网站| 成人欧美大片| 国产一区在线观看成人免费| 亚洲精品久久国产高清桃花| 中文字幕av电影在线播放| 欧美在线黄色| 国产精华一区二区三区| 日日爽夜夜爽网站| av在线播放免费不卡| 国产亚洲av嫩草精品影院| 99国产精品99久久久久| 久久精品夜夜夜夜夜久久蜜豆 | 每晚都被弄得嗷嗷叫到高潮| 色在线成人网| 国产精品自产拍在线观看55亚洲| 十分钟在线观看高清视频www| 亚洲国产精品合色在线| 久久亚洲精品不卡| 一本一本综合久久| 国产亚洲精品综合一区在线观看 | 日韩欧美国产一区二区入口| 欧美日韩乱码在线| 欧美国产日韩亚洲一区| 亚洲电影在线观看av| 国产三级黄色录像| 精品一区二区三区四区五区乱码| 亚洲国产欧美网| 久久人人精品亚洲av| 18禁美女被吸乳视频| 欧美另类亚洲清纯唯美| 亚洲电影在线观看av| 国产精品日韩av在线免费观看| 99久久国产精品久久久| 91在线观看av| 久久中文看片网| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 久久国产精品男人的天堂亚洲| 免费看十八禁软件| 国产精品电影一区二区三区| 亚洲精品在线观看二区| 母亲3免费完整高清在线观看| 亚洲欧美日韩高清在线视频| 欧美又色又爽又黄视频| 免费搜索国产男女视频| 亚洲精品av麻豆狂野| 色在线成人网| 免费看十八禁软件| 国产精品电影一区二区三区| 国产真人三级小视频在线观看| 热99re8久久精品国产| 国产黄片美女视频| 亚洲第一青青草原| 午夜精品久久久久久毛片777| 9191精品国产免费久久| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 午夜免费鲁丝| 男女午夜视频在线观看| 免费搜索国产男女视频| 久久热在线av| 嫩草影视91久久| 欧美成人性av电影在线观看| 一本精品99久久精品77| 国产高清视频在线播放一区| 18禁裸乳无遮挡免费网站照片 | 久久精品国产综合久久久| 久久久久久大精品| 美女 人体艺术 gogo| 国产精品亚洲美女久久久| 国产在线观看jvid| 色精品久久人妻99蜜桃| 亚洲 欧美一区二区三区| 欧美av亚洲av综合av国产av| 男人的好看免费观看在线视频 | 成人精品一区二区免费| 久热这里只有精品99| 两性夫妻黄色片| 日韩国内少妇激情av| 久久久精品欧美日韩精品| 中文字幕最新亚洲高清| 国产男靠女视频免费网站| 18禁裸乳无遮挡免费网站照片 | 成人av一区二区三区在线看| 制服丝袜大香蕉在线| 国产99久久九九免费精品| av在线播放免费不卡| 18禁黄网站禁片午夜丰满| 手机成人av网站| 淫秽高清视频在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲成人精品中文字幕电影| 男人操女人黄网站| 婷婷精品国产亚洲av| 国产精品爽爽va在线观看网站 | 亚洲 欧美 日韩 在线 免费| 亚洲九九香蕉| 最新美女视频免费是黄的| 美女免费视频网站| 国产极品粉嫩免费观看在线| 欧美精品亚洲一区二区| 国产精品自产拍在线观看55亚洲| 中文字幕最新亚洲高清| 亚洲成人精品中文字幕电影| 给我免费播放毛片高清在线观看| 中文字幕精品免费在线观看视频| 亚洲av成人av| 国产成人欧美在线观看| 成人永久免费在线观看视频| 国产欧美日韩一区二区三| 久久精品国产99精品国产亚洲性色| 精品不卡国产一区二区三区| 久久久久久大精品| 国内精品久久久久久久电影| 99国产精品一区二区三区| 国产黄a三级三级三级人| 国产亚洲av嫩草精品影院| 亚洲av片天天在线观看| 在线免费观看的www视频| 久久久久久亚洲精品国产蜜桃av| 不卡一级毛片| 99久久无色码亚洲精品果冻| 淫妇啪啪啪对白视频| 一区二区日韩欧美中文字幕| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 18禁国产床啪视频网站| 在线av久久热| 亚洲精品中文字幕在线视频| 一个人观看的视频www高清免费观看 | 丰满人妻熟妇乱又伦精品不卡| 精品国产亚洲在线| 欧美zozozo另类| 国产伦在线观看视频一区| 哪里可以看免费的av片| 成人18禁在线播放| 美女大奶头视频| 很黄的视频免费| 欧美亚洲日本最大视频资源| 51午夜福利影视在线观看| 男男h啪啪无遮挡| 黄片小视频在线播放| 校园春色视频在线观看| 国产成人系列免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放 | 丁香六月欧美| 亚洲精品美女久久久久99蜜臀| 人人澡人人妻人| 美女高潮喷水抽搐中文字幕| 在线播放国产精品三级| 丁香欧美五月| 亚洲人成伊人成综合网2020| 亚洲中文字幕一区二区三区有码在线看 | 后天国语完整版免费观看| 国产黄片美女视频| 日日干狠狠操夜夜爽| 久久 成人 亚洲| 欧美午夜高清在线| 国产v大片淫在线免费观看| 男人舔奶头视频| 国产欧美日韩一区二区精品| 一二三四在线观看免费中文在| 黑丝袜美女国产一区| 国产激情欧美一区二区| 国产一卡二卡三卡精品| 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 最近在线观看免费完整版| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全免费视频| 亚洲国产日韩欧美精品在线观看 | 国产欧美日韩一区二区三| 伦理电影免费视频| 美女大奶头视频| 久久久精品国产亚洲av高清涩受| 不卡一级毛片| 18禁裸乳无遮挡免费网站照片 | 亚洲av成人不卡在线观看播放网| 久久九九热精品免费| av免费在线观看网站| 亚洲av熟女| 欧美乱码精品一区二区三区| 久久精品影院6| 人人澡人人妻人| 在线免费观看的www视频| 色尼玛亚洲综合影院| 日韩精品免费视频一区二区三区| 精品人妻1区二区| 最近在线观看免费完整版| 亚洲成人免费电影在线观看| 丝袜人妻中文字幕| 亚洲色图av天堂| 亚洲国产精品成人综合色| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久国产精品人妻aⅴ院| 成在线人永久免费视频| 热re99久久国产66热| 国产一区二区激情短视频| 高潮久久久久久久久久久不卡| 日韩欧美三级三区| 观看免费一级毛片| 91麻豆精品激情在线观看国产| 色尼玛亚洲综合影院| 久热这里只有精品99| 日本三级黄在线观看| 亚洲精品久久成人aⅴ小说| 99久久综合精品五月天人人| 波多野结衣av一区二区av| 夜夜夜夜夜久久久久| 日本五十路高清| 日韩一卡2卡3卡4卡2021年| 神马国产精品三级电影在线观看 | 久久午夜亚洲精品久久| 久久久久久免费高清国产稀缺| 嫁个100分男人电影在线观看| 精品不卡国产一区二区三区| 这个男人来自地球电影免费观看| 黄色成人免费大全| 亚洲精品久久国产高清桃花| 麻豆成人av在线观看| 午夜福利在线观看吧| 麻豆成人av在线观看| 午夜福利在线观看吧| 变态另类成人亚洲欧美熟女| 18美女黄网站色大片免费观看| av免费在线观看网站| 最近最新免费中文字幕在线| 亚洲欧美精品综合一区二区三区| 亚洲黑人精品在线| 国产乱人伦免费视频| 亚洲国产欧美一区二区综合| 国产黄片美女视频| 久久精品成人免费网站| 在线观看免费视频日本深夜| 一本久久中文字幕| 亚洲 欧美 日韩 在线 免费| 午夜精品在线福利| ponron亚洲| 久久久精品国产亚洲av高清涩受| 亚洲真实伦在线观看|