• <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)之神——火炮
    国产v大片淫在线免费观看| 级片在线观看| 中国国产av一级| 成年av动漫网址| 春色校园在线视频观看| 老司机影院成人| 亚洲精品影视一区二区三区av| 成人三级黄色视频| 一级毛片aaaaaa免费看小| 国产91av在线免费观看| 性插视频无遮挡在线免费观看| 国产成人a区在线观看| 国产黄色视频一区二区在线观看 | 日本免费a在线| 久久久成人免费电影| 国产av一区在线观看免费| 天美传媒精品一区二区| 黑人高潮一二区| 天天躁夜夜躁狠狠久久av| 精品一区二区三区人妻视频| 国产欧美日韩精品一区二区| 麻豆av噜噜一区二区三区| 麻豆乱淫一区二区| 99热这里只有是精品在线观看| 欧美又色又爽又黄视频| 中文在线观看免费www的网站| 卡戴珊不雅视频在线播放| www.色视频.com| 免费在线观看影片大全网站| 一个人看视频在线观看www免费| 在线播放无遮挡| 亚洲精品粉嫩美女一区| 亚洲精品在线观看二区| 亚洲精华国产精华液的使用体验 | 99热这里只有是精品在线观看| 国产一区二区激情短视频| 久久久精品94久久精品| a级毛片a级免费在线| 亚洲无线在线观看| 久久精品国产亚洲av香蕉五月| 久久久a久久爽久久v久久| 亚洲五月天丁香| 露出奶头的视频| 久久久精品94久久精品| 国产一区二区在线观看日韩| 亚洲欧美成人综合另类久久久 | 黄色日韩在线| 国产精品1区2区在线观看.| 变态另类成人亚洲欧美熟女| 夜夜看夜夜爽夜夜摸| 欧美最新免费一区二区三区| 午夜激情福利司机影院| 哪里可以看免费的av片| av免费在线看不卡| 国产视频一区二区在线看| av.在线天堂| 夜夜看夜夜爽夜夜摸| 不卡一级毛片| 精品99又大又爽又粗少妇毛片| 亚洲自偷自拍三级| 在线观看av片永久免费下载| 国产精品不卡视频一区二区| 内地一区二区视频在线| 亚洲aⅴ乱码一区二区在线播放| 内地一区二区视频在线| 亚洲国产精品合色在线| 在线观看66精品国产| 久久精品夜色国产| 99九九线精品视频在线观看视频| 少妇高潮的动态图| 亚洲专区国产一区二区| 国产精品一区二区三区四区免费观看 | 欧美不卡视频在线免费观看| 六月丁香七月| 日韩精品青青久久久久久| 精品一区二区三区视频在线观看免费| 国产免费男女视频| 黄色日韩在线| 午夜福利在线观看免费完整高清在 | 亚洲av第一区精品v没综合| 久久婷婷人人爽人人干人人爱| 国产女主播在线喷水免费视频网站 | 少妇人妻一区二区三区视频| 国产 一区精品| 日本色播在线视频| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 天堂√8在线中文| 精品乱码久久久久久99久播| 亚洲熟妇中文字幕五十中出| 日韩三级伦理在线观看| 国产一区二区亚洲精品在线观看| 亚洲第一区二区三区不卡| 少妇人妻精品综合一区二区 | 欧美+日韩+精品| 人妻丰满熟妇av一区二区三区| 国产片特级美女逼逼视频| 日韩精品有码人妻一区| av天堂在线播放| 波多野结衣高清无吗| 国产在线精品亚洲第一网站| 三级毛片av免费| 日韩高清综合在线| 狂野欧美激情性xxxx在线观看| 亚洲激情五月婷婷啪啪| 蜜臀久久99精品久久宅男| 美女高潮的动态| 亚洲欧美日韩卡通动漫| 亚洲av一区综合| 一级a爱片免费观看的视频| 99久久久亚洲精品蜜臀av| av在线蜜桃| 亚洲成人中文字幕在线播放| 热99在线观看视频| 国产私拍福利视频在线观看| 丰满的人妻完整版| 欧美区成人在线视频| 欧美高清性xxxxhd video| 国产av一区在线观看免费| 亚洲欧美成人精品一区二区| 成年版毛片免费区| 1000部很黄的大片| 在线观看一区二区三区| 国产精品,欧美在线| 日韩中字成人| 白带黄色成豆腐渣| 免费看光身美女| 国产极品精品免费视频能看的| 亚洲av五月六月丁香网| 免费观看在线日韩| 亚洲精品在线观看二区| 免费电影在线观看免费观看| 97超视频在线观看视频| 国产精品,欧美在线| 男女做爰动态图高潮gif福利片| 国产一级毛片七仙女欲春2| 亚洲丝袜综合中文字幕| 看免费成人av毛片| 看黄色毛片网站| 国内久久婷婷六月综合欲色啪| 国产精品一区www在线观看| 99久久精品国产国产毛片| 日本爱情动作片www.在线观看 | 久久久国产成人免费| 俺也久久电影网| 欧美成人一区二区免费高清观看| 欧美一区二区亚洲| 中文字幕熟女人妻在线| 色哟哟·www| 有码 亚洲区| 哪里可以看免费的av片| 十八禁网站免费在线| 美女黄网站色视频| 久久久久久伊人网av| 九九久久精品国产亚洲av麻豆| 中文字幕av成人在线电影| 99久久久亚洲精品蜜臀av| 三级国产精品欧美在线观看| 亚洲电影在线观看av| 亚洲人成网站在线播放欧美日韩| 国产单亲对白刺激| 欧美色视频一区免费| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 国产探花极品一区二区| 日韩一本色道免费dvd| 人人妻,人人澡人人爽秒播| av免费在线看不卡| 男女下面进入的视频免费午夜| 老熟妇仑乱视频hdxx| 亚洲五月天丁香| 长腿黑丝高跟| 一个人观看的视频www高清免费观看| 香蕉av资源在线| 国产成人福利小说| 在线观看66精品国产| 美女被艹到高潮喷水动态| 久久精品国产亚洲网站| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 熟女电影av网| 久久鲁丝午夜福利片| 国产v大片淫在线免费观看| 国产精品久久久久久久久免| 午夜精品在线福利| a级毛色黄片| 亚洲四区av| 亚洲最大成人手机在线| 99热6这里只有精品| 国产精品一二三区在线看| 午夜爱爱视频在线播放| 久久午夜亚洲精品久久| 久久久国产成人免费| 麻豆乱淫一区二区| 欧美成人免费av一区二区三区| 99久国产av精品国产电影| 女人被狂操c到高潮| 夜夜夜夜夜久久久久| 亚洲av二区三区四区| 成人鲁丝片一二三区免费| a级毛片免费高清观看在线播放| 日本五十路高清| 日韩人妻高清精品专区| 日韩精品青青久久久久久| 一个人看的www免费观看视频| 精品日产1卡2卡| 国产精品久久久久久av不卡| 女人十人毛片免费观看3o分钟| 精品一区二区三区人妻视频| 国产乱人偷精品视频| 春色校园在线视频观看| 一级毛片我不卡| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产自在天天线| 丰满的人妻完整版| 波多野结衣高清作品| 久久久久久久亚洲中文字幕| 九九爱精品视频在线观看| 非洲黑人性xxxx精品又粗又长| 欧美一区二区亚洲| 国产 一区精品| 成年版毛片免费区| 伦精品一区二区三区| 国产又黄又爽又无遮挡在线| 在线免费十八禁| av在线播放精品| 亚洲精品一卡2卡三卡4卡5卡| 久久人人爽人人片av| 99久久精品国产国产毛片| 成人无遮挡网站| 国产亚洲欧美98| 中文字幕人妻熟人妻熟丝袜美| 亚洲av中文字字幕乱码综合| 日本三级黄在线观看| 日韩国内少妇激情av| 国产精品亚洲一级av第二区| 国产真实伦视频高清在线观看| av免费在线看不卡| 在线观看66精品国产| 欧美潮喷喷水| 亚洲欧美成人精品一区二区| 国产在线精品亚洲第一网站| 淫妇啪啪啪对白视频| 亚洲,欧美,日韩| 国产午夜精品论理片| 日本黄色视频三级网站网址| 日韩精品有码人妻一区| 亚洲成人久久性| 久久99热6这里只有精品| 午夜福利在线观看吧| 一个人观看的视频www高清免费观看| 亚洲欧美日韩无卡精品| 村上凉子中文字幕在线| 日本黄大片高清| 1000部很黄的大片| 日产精品乱码卡一卡2卡三| 黑人高潮一二区| 国产精品三级大全| 久久久精品欧美日韩精品| 一级毛片电影观看 | 性插视频无遮挡在线免费观看| 国产精品久久久久久av不卡| 久久精品夜色国产| 亚洲av成人精品一区久久| 国产精品1区2区在线观看.| 身体一侧抽搐| av天堂中文字幕网| 不卡一级毛片| 亚洲精品日韩在线中文字幕 | 秋霞在线观看毛片| 日本黄色片子视频| 国产精品永久免费网站| 99久久九九国产精品国产免费| 日本黄大片高清| 人妻丰满熟妇av一区二区三区| 国产成人aa在线观看| 少妇熟女欧美另类| 亚洲图色成人| 国产日本99.免费观看| 色视频www国产| 如何舔出高潮| 欧美成人免费av一区二区三区| 三级毛片av免费| 夜夜夜夜夜久久久久| 色尼玛亚洲综合影院| 日本一本二区三区精品| 国产色爽女视频免费观看| 亚洲图色成人| 亚洲在线自拍视频| 亚洲欧美日韩卡通动漫| 热99在线观看视频| 一a级毛片在线观看| 国内揄拍国产精品人妻在线| 一进一出好大好爽视频| 国产男人的电影天堂91| 免费高清视频大片| 亚洲熟妇熟女久久| 国产av不卡久久| 亚洲一级一片aⅴ在线观看| 香蕉av资源在线| 日韩人妻高清精品专区| 久久久午夜欧美精品| 少妇人妻一区二区三区视频| 久久欧美精品欧美久久欧美| 午夜福利高清视频| 一区二区三区高清视频在线| 精品人妻偷拍中文字幕| 99国产极品粉嫩在线观看| 久久欧美精品欧美久久欧美| 国产一级毛片七仙女欲春2| 在线a可以看的网站| 亚洲无线观看免费| 一级黄片播放器| 欧美xxxx性猛交bbbb| 少妇裸体淫交视频免费看高清| 久久精品国产清高在天天线| 婷婷精品国产亚洲av在线| 22中文网久久字幕| 国产精品福利在线免费观看| 国产三级在线视频| 成人国产麻豆网| 精品国产三级普通话版| 我的女老师完整版在线观看| 日韩欧美三级三区| 有码 亚洲区| 少妇裸体淫交视频免费看高清| 日本免费a在线| 成人午夜高清在线视频| 最后的刺客免费高清国语| 亚洲最大成人手机在线| 夜夜夜夜夜久久久久| 又黄又爽又刺激的免费视频.| 精品久久久久久久久久免费视频| 亚洲欧美精品综合久久99| 亚洲天堂国产精品一区在线| 日本在线视频免费播放| 久久久精品大字幕| 亚洲精华国产精华液的使用体验 | 少妇的逼好多水| 亚洲精品粉嫩美女一区| 乱人视频在线观看| 亚洲欧美清纯卡通| 免费在线观看成人毛片| 亚洲国产精品久久男人天堂| 亚洲欧美日韩高清在线视频| 又黄又爽又免费观看的视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产日韩欧美精品在线观看| 国产色爽女视频免费观看| 天堂网av新在线| 亚洲成av人片在线播放无| 伦精品一区二区三区| 你懂的网址亚洲精品在线观看 | 久久亚洲国产成人精品v| 人妻夜夜爽99麻豆av| www日本黄色视频网| 久久中文看片网| a级毛片a级免费在线| 国产亚洲精品久久久com| 女人被狂操c到高潮| 国产av一区在线观看免费| 在线国产一区二区在线| 变态另类成人亚洲欧美熟女| 欧美一区二区精品小视频在线| 欧美bdsm另类| 日韩成人av中文字幕在线观看 | 最近在线观看免费完整版| 夜夜夜夜夜久久久久| 亚洲精品乱码久久久v下载方式| 蜜桃久久精品国产亚洲av| 联通29元200g的流量卡| 国产色爽女视频免费观看| 老师上课跳d突然被开到最大视频| 国产高清三级在线| 欧美最黄视频在线播放免费| 国产精品永久免费网站| 国内精品一区二区在线观看| 日韩一本色道免费dvd| 一区福利在线观看| 一区福利在线观看| 亚洲无线在线观看| 欧美xxxx黑人xx丫x性爽| 精品一区二区三区视频在线| 国产片特级美女逼逼视频| 欧美国产日韩亚洲一区| 少妇的逼水好多| 日本成人三级电影网站| 波多野结衣高清无吗| 日韩制服骚丝袜av| 国产精品三级大全| 俺也久久电影网| 亚洲欧美精品综合久久99| 成人一区二区视频在线观看| 国产精品国产高清国产av| 午夜福利高清视频| 国产高清视频在线观看网站| av视频在线观看入口| 国产欧美日韩一区二区精品| 女的被弄到高潮叫床怎么办| 国内精品久久久久精免费| 在线观看免费视频日本深夜| 国内少妇人妻偷人精品xxx网站| 91av网一区二区| 国产一区二区亚洲精品在线观看| or卡值多少钱| 成人无遮挡网站| 免费看av在线观看网站| 精品久久久噜噜| 欧美最新免费一区二区三区| 婷婷精品国产亚洲av在线| 国产精品一区二区三区四区久久| 国产精品女同一区二区软件| 日本一二三区视频观看| 欧洲精品卡2卡3卡4卡5卡区| 国产激情偷乱视频一区二区| 亚洲第一区二区三区不卡| 国产精品无大码| АⅤ资源中文在线天堂| 久久久久国内视频| 国产中年淑女户外野战色| av黄色大香蕉| 亚洲av二区三区四区| 亚洲精品一区av在线观看| 国产国拍精品亚洲av在线观看| 99国产精品一区二区蜜桃av| 午夜免费激情av| 99久久成人亚洲精品观看| av在线亚洲专区| 一区二区三区高清视频在线| 免费不卡的大黄色大毛片视频在线观看 | 亚洲一区高清亚洲精品| 国产欧美日韩精品一区二区| 欧美最黄视频在线播放免费| av福利片在线观看| 国产 一区精品| 欧美性猛交╳xxx乱大交人| a级毛片免费高清观看在线播放| av卡一久久| 久久国内精品自在自线图片| 国产av一区在线观看免费| 精品福利观看| 又粗又爽又猛毛片免费看| 欧美精品国产亚洲| 男插女下体视频免费在线播放| 99热6这里只有精品| 狂野欧美白嫩少妇大欣赏| 啦啦啦韩国在线观看视频| 国产精品乱码一区二三区的特点| 永久网站在线| 国产成年人精品一区二区| 日日干狠狠操夜夜爽| 国产av一区在线观看免费| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| 成人鲁丝片一二三区免费| 欧美+日韩+精品| 日本色播在线视频| 日韩成人av中文字幕在线观看 | 免费电影在线观看免费观看| 日本爱情动作片www.在线观看 | 中国美白少妇内射xxxbb| 日韩精品青青久久久久久| 舔av片在线| 人妻少妇偷人精品九色| 午夜久久久久精精品| 黄色日韩在线| 国内揄拍国产精品人妻在线| 日韩av在线大香蕉| 性欧美人与动物交配| 赤兔流量卡办理| 嫩草影视91久久| 99在线人妻在线中文字幕| 欧美xxxx黑人xx丫x性爽| 级片在线观看| 插逼视频在线观看| 久久99热这里只有精品18| 六月丁香七月| 久久久国产成人精品二区| 亚洲成av人片在线播放无| 亚洲精品亚洲一区二区| 亚洲无线在线观看| 一级黄片播放器| 久久久久久久久中文| 三级国产精品欧美在线观看| 免费电影在线观看免费观看| 午夜免费激情av| 亚洲精品乱码久久久v下载方式| 直男gayav资源| 久久久久国产网址| 91久久精品电影网| 精品久久国产蜜桃| 三级毛片av免费| 中文字幕久久专区| 1000部很黄的大片| 国产精品人妻久久久影院| 欧美色视频一区免费| 精品99又大又爽又粗少妇毛片| 小蜜桃在线观看免费完整版高清| 免费在线观看影片大全网站| 国产黄色小视频在线观看| 99热精品在线国产| 男女做爰动态图高潮gif福利片| 亚洲国产色片| 在线观看免费视频日本深夜| 日韩强制内射视频| av在线老鸭窝| 国产高清激情床上av| 女人被狂操c到高潮| 悠悠久久av| 在线观看午夜福利视频| 欧美绝顶高潮抽搐喷水| 天堂√8在线中文| 能在线免费观看的黄片| 国产精品99久久久久久久久| 久久人人精品亚洲av| 欧美一区二区亚洲| 校园春色视频在线观看| 欧美日韩综合久久久久久| 久久精品国产清高在天天线| 国产精品久久视频播放| 成人美女网站在线观看视频| 熟女电影av网| 长腿黑丝高跟| 国产真实乱freesex| 午夜视频国产福利| 亚洲精品日韩在线中文字幕 | 国产亚洲精品综合一区在线观看| 三级国产精品欧美在线观看| 香蕉av资源在线| 麻豆成人午夜福利视频| 看黄色毛片网站| 嫩草影视91久久| 久久6这里有精品| 国产aⅴ精品一区二区三区波| 亚洲av五月六月丁香网| 国产精品免费一区二区三区在线| 一级毛片我不卡| 日韩一区二区视频免费看| 日本色播在线视频| 性插视频无遮挡在线免费观看| 欧美高清性xxxxhd video| 免费av不卡在线播放| 亚州av有码| a级毛色黄片| 九九在线视频观看精品| 十八禁网站免费在线| 最近手机中文字幕大全| 别揉我奶头 嗯啊视频| 麻豆久久精品国产亚洲av| 久久综合国产亚洲精品| 又爽又黄a免费视频| 寂寞人妻少妇视频99o| 欧美性感艳星| 真实男女啪啪啪动态图| 国产一级毛片七仙女欲春2| 丰满乱子伦码专区| 亚洲欧美中文字幕日韩二区| 亚洲三级黄色毛片| 久久久久久九九精品二区国产| 久久精品影院6| 国产精品久久电影中文字幕| 久久热精品热| 中国美女看黄片| 久久久久免费精品人妻一区二区| 婷婷色综合大香蕉| 国产真实乱freesex| 精品福利观看| 色综合站精品国产| 成人鲁丝片一二三区免费| 精品少妇黑人巨大在线播放 | 亚洲一区二区三区色噜噜| 五月玫瑰六月丁香| 日韩 亚洲 欧美在线| 日韩精品中文字幕看吧| 有码 亚洲区| 亚洲图色成人| 自拍偷自拍亚洲精品老妇| 午夜福利视频1000在线观看| 亚洲,欧美,日韩| 美女黄网站色视频| 欧美丝袜亚洲另类| 天堂av国产一区二区熟女人妻| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 男女那种视频在线观看| 精品人妻视频免费看| 亚洲最大成人手机在线| 亚洲国产精品合色在线| 一级毛片电影观看 | 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 国产麻豆成人av免费视频| av天堂中文字幕网| 日韩成人av中文字幕在线观看 | 九九在线视频观看精品| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 两个人的视频大全免费| av视频在线观看入口| 欧美日本视频| 久久久久久久久大av| 国产亚洲欧美98| av在线蜜桃| 五月玫瑰六月丁香| 秋霞在线观看毛片| 亚洲欧美日韩高清在线视频| 国产高清激情床上av| 1024手机看黄色片| 看十八女毛片水多多多| 黑人高潮一二区| 别揉我奶头 嗯啊视频| 午夜福利在线在线| 插阴视频在线观看视频| 中文资源天堂在线| 国产男人的电影天堂91|