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

    基于頻響函數(shù)提升小波總能量的桁架結構隨機模型修正

    2021-07-01 06:07:24趙永鵬彭珍瑞
    計算力學學報 2021年3期
    關鍵詞:頻響桁架小波

    趙永鵬, 殷 紅, 彭珍瑞

    (蘭州交通大學 機電工程學院,蘭州 730070)

    1 引 言

    隨著工程技術領域的不斷發(fā)展,有限元技術的應用更加廣泛,但是各種不確定性因素,如建模時的簡化及安裝誤差等對結果的影響很大,因此建立一個可靠有效的有限元模型就變得尤為重要。當前,通過模型修正來獲得準確的有限元模型(FEM)是主流的方法。模型修正方法分為確定性和不確定性兩種[1],確定性方法在近年來有了廣泛發(fā)展。但是,由于各種不確定性的情況,出現(xiàn)的誤差較多。所以,在模型修正中充分考慮不確定性的影響,從而構建一種隨機性的模型修正方法勢在必行。

    參數(shù)不確定性的模型修正方法指通過多次模態(tài)實驗得到相應數(shù)據(jù),計算出具有統(tǒng)計特征的有限元模型,來表明實際結構的靜動力學行為、機理和演化規(guī)律等。方圣恩等[2]充分考慮了修正參數(shù)的不確定性,通過計算參數(shù)的統(tǒng)計特征量對一組鋼板進行了修正。Fang等[3]通過Hermite多項式的顯式多項式混沌展開,把不確定的參數(shù)和響應通過構造一個隨機響應面模型(RSM)對一組金屬板進行了修正。Zhai等[4]以模態(tài)參數(shù)為響應,基于改進的響應面模型和先進的蒙特卡洛(MC)方法,對一航空發(fā)動機定子系統(tǒng)(殼體)進行了修正。Bi等[5]構建了一個完全嵌入Bhattacharyya距離的近似貝葉斯模型的修正框架,對一質量彈簧系統(tǒng)進行了模型修正,證明了方法的有效性。Deng等[6]通過構建代理模型、設置校準參數(shù)以及將距離函數(shù)作為參數(shù)相關性的衡量標準等步驟進行了隨機模型修正。Khodaparast等[7]利用模態(tài)響應(如固有頻率和模態(tài)振型)中得到的可變性來估計隨機修正參數(shù)的第一和第二統(tǒng)計矩,通過一個三自由度的質量彈簧系統(tǒng)對所提方法進行了驗證。

    通過結構模態(tài)參數(shù)進行修正將不可避免地在識別的過程中引入一些誤差,實際工程中由于實驗模態(tài)參數(shù)的識別而出現(xiàn)的誤差有可能會在某些特定情況下比建模精度低引起的誤差更大。而基于頻響函數(shù)(FRFs)的模型修正方法因為可以略去結構模態(tài)參數(shù)識別而得到廣泛應用,該方法減少了各種誤差的出現(xiàn),適于對模態(tài)分布相對比較密的結構進行修正。李偉明等[8]利用頻響函數(shù)為輸出對一個二維桁架進行了修正。王巨濤等[9]利用加速度頻響函數(shù)曲線對應頻率點處的響應值之差構造目標函數(shù),對GARTEUR飛機模型進行了修正。但是頻響函數(shù)包含的數(shù)據(jù)信息量很大,怎樣選出好的頻率點數(shù)據(jù)以得到期望的修正效果,現(xiàn)在還沒有形成統(tǒng)一的標準。因此,可以考慮引入一個特征向量對頻響函數(shù)進行表征,以避免上述問題。張勇等[10]提取頻響函數(shù)的特征向量對一扭力梁模型進行了修正。上述研究大都應用于確定性的模型修正,但鮮見應用于不確定性模型修正。提升小波變換LWT(Lifting Wavelet Transform)屬于第二代小波變換,與傳統(tǒng)小波變換相比,提升小波變換不僅擁有多分辨率的特點,而且具有計算時間少,整個流程相對簡單以及變換后數(shù)據(jù)長度不變等優(yōu)勢[11]。宋萌萌等[12]選擇提升小波變換對得到的初始信號進行降低噪聲的處理,得到的結果滿足要求。袁海英等[13]利用提升小波變換對振動信號進行時頻特性分析和信息預處理,通過預測器和更新器的設計取代了小波基函數(shù)選取過程。提升小波變換不再使用傅立葉變換,也沒有利用函數(shù)的平移或伸縮,而是在時域中就可以完成原始信號不同頻帶上的分離。經過變換后各層分量的小波能量對信號的變化有很高的敏感性,且對噪聲污染、小波種類選擇及小波層數(shù)具有魯棒性。因此,選擇提升小波變換對頻響函數(shù)進行分解,提取總能量值來表征頻響函數(shù)可以最大程度上保證頻響函數(shù)的特征不發(fā)生變化,而且也減少了頻響函數(shù)信息過多引起的各種處理困難。

    本文構建了一種以提升小波總能量為響應,通過篩選響應樣本反求出待修正參數(shù)的隨機模型修正方法。首先對加速度頻響函數(shù)進行提升小波變換,提取對應的提升小波總能量作為輸出,以待修正參數(shù)作為輸入,構造一個多項式響應面代理模型代替有限元模型;然后利用蒙特卡洛抽樣得到大量的響應樣本,并通過設置閾值對響應樣本進行優(yōu)化篩選;最后構造一個優(yōu)化反問題,利用布谷鳥優(yōu)化算法反求出輸入即待修正參數(shù),計算修正后參數(shù)的統(tǒng)計特征。分別對平面桁架和三維鋼桁架結構模型進行修正,證明了所構建方法的有效性。

    2 頻響函數(shù)提升小波總能量

    2.1 加速度頻響函數(shù)

    在具有阻尼的系統(tǒng)中,通過一個簡諧激勵作用所得的動力學方程為

    (1)

    在頻域內的輸入輸出關系可以表示為

    A(ω)=H(ω)F(ω)

    (2)

    式中A(ω)為穩(wěn)態(tài)響應,H(ω)為加速度頻響函數(shù),F(xiàn)(ω)為簡諧激勵,ω為激勵頻率。

    相應的加速度頻響函數(shù)可以表示為

    (3)

    2.2 提升小波變換

    對于得到的頻響函數(shù),選擇提升小波變換對其能量值進行提取。因為常見的小波濾波器的輸出是浮點數(shù),所以對變換后的值選擇壓縮時需要量化才能得出對應的整數(shù)。Swelden[14]選擇了根據(jù)提升方式改進小波特性得到小波基函數(shù)的新模式即提升小波變換。提升小波變換共分為三個過程: 分裂、預測和更新。

    (1) 分裂。根據(jù)原始信號sj奇偶性的不同可以得到兩個沒有相關性的子集。所得子集的長度是初始信號的1/2。如式(4)所示:

    sj=(ej - 1,oj - 1)

    (4)

    式中ej - 1={ej - 1,k=sj,2k}為偶序列,oj - 1={oj - 1,k=sj,2k + 1}為奇序列。

    (2) 預測。根據(jù)第一步提取的偶序列ej - 1來預測奇序列oj - 1。如式(5)所示,用實際值oj - 1與預測值ej - 1的差值dj - 1來表示其接近程度,稱作細節(jié)系數(shù),同時也屬于原始信號的高頻分量。

    dj - 1=oj - 1-P(ej - 1)

    (5)

    式中P為預測函數(shù),用相鄰數(shù)據(jù)的平均值來表示:

    pk(ej - 1)=(ej - 1,k+ej - 1,k + 1)/2

    (6)

    (3) 更新。這一步的目的是要使第二步產生的某些特征(如均值等)與原始信號的特征保持不變。可以通過算子U來實現(xiàn):

    sj - 1=ej - 1+U(dj - 1)

    (7)

    式中sj - 1為原始信號的低頻部分,U為更新算子函數(shù),可以表示為Uk(dj - 1)=dj - 1,k/2。

    通過提升小波變換之后,可以將原始信號sj變?yōu)楦哳l分量dj - 1和低頻分量sj - 1;將低頻部分sj - 1繼續(xù)進行相同的步驟,又可以分解為二層高頻dj - 2和二層低頻sj - 2;以此類推,可以根據(jù)要求進行n層的提升小波分解,原始信號就可以表示為{sj - n,dj - n,…,dj - 1}。其中sj - n為信號的低頻部分,{dj - n,…,dj - 1}為信號的高頻部分。

    對頻響函數(shù)進行提升小波變換后,即可得出各層的高頻部分和低頻部分。設Ej d為原始信號sj在尺度n上的高頻分量總能量:

    (8)

    相應的低頻分量能量Ej s也可以得到,信號的提升小波總能量為各個尺度下的小波能量的總和[15]。

    Ez=E1d+E2d+…+Ej d+Ej s

    (9)

    3 響應面模型

    對于一些大型復雜結構,由于其內部存在很多的零部件,以及很多邊界尺寸數(shù)據(jù)的不確定性,所以在有限元建模的過程中存在困難,從而建立的有限元模型變得十分復雜,在模型修正過程中可能會出現(xiàn)計算無法收斂或者所用時間很長等問題。選擇一個代理模型來代替這個復雜的有限元模型,可以方便計算,提高效率。多項式響應面代理模型具有建立相對簡單且預測精度高的優(yōu)點,因此本文選擇建立響應面模型來代替有限元模型。這個代理模型將一個物理系統(tǒng)的輸入?yún)?shù)x和輸出響應y通過數(shù)學多項式結合起來:

    y=f(x1,x2,…,xn)+ε

    (10)

    式中f為輸入與輸出間的映射,ε為建模誤差,n為輸入?yún)?shù)的數(shù)量。

    對于一般性的工程問題,大多采用二階多項式模型作為一種基本的形式,即

    (11)

    式中bo,bi,bi i和bi j為回歸系數(shù),xi和xj為輸入?yún)?shù)。

    響應面構建完成后,需對其進行精度檢驗,判斷是否滿足后續(xù)計算要求,常用的檢驗方法為計算真實值與預測值間的均方根誤差,即

    (12)

    因此,根據(jù)提取的輸出與輸入?yún)?shù)就可以構建二階響應面代理模型,然后再對其進行精度檢驗。

    4 響應樣本獲取

    4.1 蒙特卡洛抽樣

    蒙特卡洛方法提供了一種可替代分析數(shù)學評估統(tǒng)計量在隨機樣本中行為的方法。根據(jù)樣本的概率分布,用計算機實現(xiàn)統(tǒng)計模擬或抽樣,屬于一種隨機抽樣方法。常見的數(shù)值計算方法由于需要多次迭代,導致最后誤差累積很大,影響結果,而蒙特卡洛抽樣雖然也有迭代,但是沒有像其他方法那樣頻繁迭代,產生的誤差一般可以忽略不計,對最后的抽樣結果也是基本沒有影響。因此,在已知響應概率分布的情況下,運用蒙特卡洛抽樣可以得到大量的響應樣本,認為這些樣本就是試驗值,然后用于隨機模型修正過程。

    4.2 響應樣本篩選

    通過蒙特卡洛抽樣得到大量響應樣本,由于是一種隨機的抽樣方法,所以可能會存在抽取的樣本波動范圍過大,以至于超出了實際試驗所能測得的范圍,樣本稱之為失真樣本。如果將這些數(shù)據(jù)也代入修正過程,就失去了實際的意義,可能導致修正的結果很不理想。因此,可以通過設定一個閾值來對抽樣數(shù)據(jù)進行一次優(yōu)化篩選,使得篩選后的數(shù)據(jù)與實際情況更加接近。本文通過計算抽樣響應值與真實響應值間的均方根誤差進行樣本的篩選,

    (13)

    為了保證所抽取樣本的真實性并且選擇有效的樣本,按照至少保留80%左右樣本的準則進行篩選。計算所有的RMSE值,根據(jù)該準則設定一個閾值,當RMSE值比閾值小時保留,反之則淘汰。

    5 模型修正

    根據(jù)上述理論,通過響應面模型代替有限元模型后,以蒙特卡洛抽樣后優(yōu)選的響應值與響應面模型預測響應值差值最小來構造目標函數(shù)。

    (14)

    建立目標函數(shù)之后,模型修正變成對目標函數(shù)求解的問題。由于布谷鳥優(yōu)化算法具有全局搜索能力強和參數(shù)少等優(yōu)勢,因此本文選擇該算法進行尋優(yōu)以得到修正后的參數(shù)值。

    綜合頻響函數(shù)提升小波總能量的提取,響應面模型的建立,響應樣本的篩選,目標函數(shù)的求解等步驟,可得隨機模型修正流程如圖1所示。

    圖1 模型修正流程

    6 數(shù)值算例

    6.1 平面桁架

    如圖2所示平面桁架,該桁架結構共有14個節(jié)點和25個自由度。桿單元的彈性模量為210 GPa,材料密度為7850 kg/m3。將該平面桁架結構作為試驗模型,選擇彈性模量和密度為待修正參數(shù)。對參數(shù)進行偏離,其中密度減小10%,彈性模量增加10%作為相應的有限元模型列入表1。選擇第11節(jié)點y方向為激勵點位置,第3節(jié)點y方向為測點位置。

    圖2 平面桁架模型

    表1 試驗模型與有限元模型參數(shù)

    選擇拉丁超立方法獲取樣本點。在待修正參數(shù)的±20%區(qū)間內抽取40個樣本點。對于抽樣所得40組樣本,選擇前20組為樣本集擬合響應面模型,后20組為測試集,用來驗證代理模型的精度。通過回歸擬合出響應面模型如圖3所示。

    圖3 平面桁架的響應面模型

    假設試驗值均服從正態(tài)分布。依據(jù)抽樣所得樣本點,可以得到其加速度頻響函數(shù)值,對其進行三層的提升小波變換,提取相應的提升小波總能量。根據(jù)蒙特卡洛方法完成抽樣,共抽取200個樣本。對于抽取的樣本通過設定閾值進行篩選,按照至少保留80%數(shù)據(jù)篩選過后,代入修正過程進行求解。

    對構建的響應面模型進行精度檢驗后,所得均方根誤差RMSE值如圖4所示??梢钥闯?,RMSE值最大不超過0.06。由此可得,所構建的響應面模型精度滿足計算要求,該模型有效。

    圖4 響應面模型RMSE值

    通過上述方法進行模型修正,修正結果列入 表2。由表2可知,修正誤差小于2.1%,達到了較高的修正精度。

    為了進一步說明方法的有效性,修正前后的頻響函數(shù)曲線如圖5所示。可以看出,修正后的頻響函數(shù)曲線與真實值的曲線重合度高,證明用提升小波總能量來表征頻響函數(shù)作為結構響應進行模型修正是可行的。

    圖5 修正前后頻響函數(shù)曲線

    6.2 三維桁架

    如圖6所示三維桁架。該桁架共有28個節(jié)點,66個單元,48個自由度。總長2.8 m,寬0.39 m,高0.27 m。桁架共四個支座固定(節(jié)點編號為1,8,9,16),節(jié)點鉸接。桁架的材料參數(shù)彈性模量為190 GPa,密度為7800 kg/m3。尺寸參數(shù)上弦桿和下弦桿橫截面積為85.5 mm2,直腹桿橫截面積為141.0 mm2,斜腹桿橫截面積為45.0 mm2。結構的激勵位置和響應位置分別如圖6所示。

    圖6 三維桁架模型

    6.2.1 試驗設計與樣本選取

    選擇如圖6所示三維桁架結構作為本算例的試驗模型,相應的有限元模型通過對待修正參數(shù)進行一定范圍的偏離得到。選擇材料參數(shù)的彈性模量和密度以及尺寸參數(shù)的上弦桿橫截面積共三個參數(shù)為待修正參數(shù)。對參數(shù)進行偏離,其中密度減小10%,彈性模量增加10%,上弦桿橫截面積減少10%,有限元模型參數(shù)列入表3。

    表3 試驗模型與有限元模型參數(shù)

    建立響應面模型,樣本點的選擇很重要。常用的試驗設計方法有中心復合法、全因子法和拉丁超立方法等。本文選擇拉丁超立方法獲取樣本點,該方法采樣比較均勻且樣本具有代表性,符合模型建立的要求。在三個待修正參數(shù)的±20%區(qū)間內利用拉丁超立方抽樣抽取40個樣本點。

    假設試驗值均服從正態(tài)分布。根據(jù)蒙特卡洛方法完成抽樣,共抽取2000個樣本。對于抽取的樣本通過設定閾值進行篩選,按照至少保留80%數(shù)據(jù)(即1600個)的要求篩選后,代入修正過程進行求解。

    6.2.2 提升小波總能量的提取

    依據(jù)拉丁超立方抽樣所得樣本點,可以得到其加速度頻響函數(shù)值,如圖7所示。以試驗值的頻響函數(shù)曲線為例,對其進行三層的提升小波變換,可以求得各層高頻分量和低頻分量如圖8所示,然后提取相應的提升小波能量,將提取的各層能量求和,就可以得到提升小波總能量。

    圖7 試驗模型頻響函數(shù)曲線

    圖8 三層提升小波變換分量

    6.2.3 響應面代理模型的構造及檢驗

    對于抽樣所得40組樣本,選擇前20組為樣本集作為響應面模型的輸入x,后20組為測試集,用來驗證代理模型的精度。將20組輸入與其所對應的輸出(即頻響函數(shù)提升小波總能量)通過回歸擬合求出響應面模型的各個參數(shù),得到代理模型表達式為

    y=18.2007-0.0011x1-0.0091x2-215.08x3+

    0.0013x1x3+0.043x2x3

    (15)

    將后20組樣本代入計算得到相應的輸出,然后對構建的響應面模型進行精度檢驗。所得均方根誤差RMSE值如圖9所示。均方根誤差越接近0表明響應面精度越高。可以看出,RMSE值不超過0.07,所構建的響應面模型精度滿足計算要求,該模型有效。

    圖9 響應面模型RMSE值

    6.2.4 模型修正

    通過蒙特卡洛抽樣得到2000個響應樣本如 圖10 所示,通過多次預實驗發(fā)現(xiàn)設定閾值為0.092時滿足篩選準則要求。對所有樣本排序后進行篩選,經過篩選后共有1643個樣本符合要求。如圖11所示為篩選后的樣本。將這些樣本依次代入目標函數(shù)中,通過布谷鳥算法對目標函數(shù)求得最小值,可得到對應的1643組輸入即為修正后的參數(shù)。算出修正后各參數(shù)的均值,修正前后的參數(shù)均值列入表4。

    圖10 蒙特卡洛抽取的響應值

    圖11 優(yōu)化篩選的響應值

    表4 修正前后參數(shù)均值

    由表4可知,修正后參數(shù)均值的誤差相比于修正前有了明顯的減小,說明修正效果良好。修正前后的提升小波總能量值列入表5。由表5可知,提升小波總能量的誤差比修正前也有了大幅的減小,表明該方法完成隨機模型修正是可行的。采用修正后的參數(shù)均值可得其加速度頻響函數(shù)如 圖12 所示。圖中修正后的頻響函數(shù)曲線與真實值的曲線基本重合,相較于修正前的曲線,修正效果達到了預期的目的。

    表5 修正前后提升小波總能量

    圖12 修正前后頻響函數(shù)曲線

    圖13和圖14分別給出修正前后頻響函數(shù)實部和虛部曲線,可以看出,修正值與真實值的實部虛部曲線基本重合,進一步說明修正結果良好。

    圖13 頻響函數(shù)實部曲線

    圖14 頻響函數(shù)虛部曲線

    為了進一步驗證所提方法的有效性,以平面桁架為例,分別將頻響函數(shù)和提升小波總能量作為響應進行模型修正,結果列入表6,所有結果均是在同一臺計算機上運行20組??梢钥闯?,將提升小波總能量作為響應進行模型修正所用時間更短;在修正精度方面,彈性模量修正精度有了很大提高,雖然材料密度修正精度略有降低,但是總體來說效率有了一定提高,計算量也有所減小,能夠達到隨機模型修正的目的。

    表6 誤差與時間對比

    7 結 論

    (1) 本文構建了一種基于加速度頻響函數(shù)提升小波總能量的隨機模型修正方法,提取頻響函數(shù)的提升小波總能量作為響應,應用于隨機模型修正過程,既利用了頻響函數(shù),又減小了計算量,提高了計算效率,修正的結果達到預期目標。

    (2) 對于蒙特卡洛抽樣得到的響應樣本,通過設定閾值進行了一次優(yōu)化篩選,充分考慮了修正過程的隨機性,保證了樣本可靠性的同時也使后續(xù)的修正過程變得更加容易。

    (3) 通過平面桁架和三維桁架模型對所提方法進行驗證,修正后參數(shù)均值的最大誤差小于 3.3%,對應點的頻響函數(shù)曲線吻合良好,達到了隨機模型修正的目的。

    猜你喜歡
    頻響桁架小波
    桁架式吸泥機改造
    構造Daubechies小波的一些注記
    科技風(2021年19期)2021-09-07 14:04:29
    基于分塊化頻響函數(shù)曲率比的砌體房屋模型損傷識別研究
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    擺臂式復合桁架機器人的開發(fā)
    美團外賣哥
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    Loader軸在雙機桁架機械手上的應用
    頻響函數(shù)殘差法在有限元模型修正中的應用
    頻響阻抗法診斷變壓器繞組變形
    国产欧美另类精品又又久久亚洲欧美| 国产精品国产三级专区第一集| 尤物成人国产欧美一区二区三区| 一本一本综合久久| 欧美少妇被猛烈插入视频| 国产成人精品久久久久久| 一本一本综合久久| 九色成人免费人妻av| 不卡视频在线观看欧美| 国产一级毛片在线| 国模一区二区三区四区视频| 在线观看三级黄色| 亚洲国产av新网站| 少妇丰满av| 欧美日韩在线观看h| 黄片无遮挡物在线观看| 国产精品.久久久| 久久精品国产亚洲av涩爱| 成人综合一区亚洲| 一级毛片我不卡| 精品久久久久久电影网| 成年美女黄网站色视频大全免费 | 交换朋友夫妻互换小说| av黄色大香蕉| 一区二区三区乱码不卡18| 国产精品欧美亚洲77777| 一级片'在线观看视频| 精品少妇黑人巨大在线播放| 亚洲av在线观看美女高潮| 如何舔出高潮| 久久99热6这里只有精品| 久久6这里有精品| 夫妻午夜视频| 免费久久久久久久精品成人欧美视频 | 色网站视频免费| 亚洲va在线va天堂va国产| 亚洲美女视频黄频| 视频区图区小说| 99国产精品免费福利视频| av视频免费观看在线观看| 亚洲精品国产av蜜桃| 久久久精品94久久精品| 国产精品国产三级国产av玫瑰| 国产伦在线观看视频一区| 又爽又黄a免费视频| 国产亚洲精品久久久com| 一级二级三级毛片免费看| 777米奇影视久久| www.色视频.com| 亚洲欧洲国产日韩| 亚洲美女黄色视频免费看| 久久99蜜桃精品久久| 欧美亚洲 丝袜 人妻 在线| 亚洲精品亚洲一区二区| 久久精品人妻少妇| 欧美精品一区二区大全| 日韩一区二区视频免费看| 免费看av在线观看网站| 男人和女人高潮做爰伦理| 国产淫片久久久久久久久| 大码成人一级视频| 性高湖久久久久久久久免费观看| 2021少妇久久久久久久久久久| 91aial.com中文字幕在线观看| 日韩一区二区视频免费看| 一级av片app| 水蜜桃什么品种好| 精品久久久精品久久久| 啦啦啦中文免费视频观看日本| 日韩三级伦理在线观看| 国产极品天堂在线| 久久久久久久久久久免费av| 人妻 亚洲 视频| 高清不卡的av网站| 中文精品一卡2卡3卡4更新| 高清黄色对白视频在线免费看 | 少妇熟女欧美另类| 亚洲人与动物交配视频| 在线免费十八禁| 亚洲综合精品二区| 成年人午夜在线观看视频| a 毛片基地| 大片电影免费在线观看免费| av播播在线观看一区| av天堂中文字幕网| 欧美日本视频| 青春草视频在线免费观看| 在线观看一区二区三区激情| 黄色欧美视频在线观看| 极品少妇高潮喷水抽搐| 特大巨黑吊av在线直播| 少妇熟女欧美另类| 久久精品久久久久久久性| av不卡在线播放| 好男人视频免费观看在线| 99久国产av精品国产电影| 99热6这里只有精品| 丰满少妇做爰视频| 91精品国产国语对白视频| 国产成人a∨麻豆精品| 日本猛色少妇xxxxx猛交久久| 中国三级夫妇交换| 久久精品久久久久久噜噜老黄| 国产 一区精品| 亚洲性久久影院| av福利片在线观看| 精品人妻一区二区三区麻豆| 老司机影院成人| 春色校园在线视频观看| 亚洲av福利一区| 国产精品成人在线| 国产熟女欧美一区二区| 亚洲色图av天堂| 99久久精品一区二区三区| 蜜桃久久精品国产亚洲av| 水蜜桃什么品种好| 久久精品国产亚洲网站| 亚洲av综合色区一区| 嘟嘟电影网在线观看| 直男gayav资源| 亚洲图色成人| 如何舔出高潮| 免费黄频网站在线观看国产| 成年女人在线观看亚洲视频| 尾随美女入室| 亚洲精华国产精华液的使用体验| 成人一区二区视频在线观看| 国产精品蜜桃在线观看| 我要看日韩黄色一级片| 人人妻人人爽人人添夜夜欢视频 | 亚洲国产av新网站| 国产爱豆传媒在线观看| 又爽又黄a免费视频| 欧美成人a在线观看| 欧美日韩视频精品一区| 日韩国内少妇激情av| 中文资源天堂在线| 日韩成人伦理影院| 亚洲性久久影院| 一个人看视频在线观看www免费| 日韩视频在线欧美| 精品久久久噜噜| 香蕉精品网在线| 少妇 在线观看| av卡一久久| 午夜福利在线在线| 在线观看人妻少妇| 在现免费观看毛片| 在现免费观看毛片| 性色av一级| 丝袜脚勾引网站| 51国产日韩欧美| 国产精品久久久久久久久免| 亚洲国产成人一精品久久久| av黄色大香蕉| 成年免费大片在线观看| 在线观看国产h片| 亚洲欧洲国产日韩| 寂寞人妻少妇视频99o| 午夜福利高清视频| 午夜福利在线在线| 国产成人freesex在线| 成年av动漫网址| 成年免费大片在线观看| 18禁动态无遮挡网站| 又爽又黄a免费视频| 日韩欧美一区视频在线观看 | 亚洲精品国产av蜜桃| 免费看日本二区| 国产乱人偷精品视频| 久久久久久久久久久丰满| 中文字幕亚洲精品专区| 免费在线观看成人毛片| 欧美 日韩 精品 国产| 久久精品国产亚洲网站| 欧美 日韩 精品 国产| 精品久久久久久久久亚洲| 日韩伦理黄色片| 亚洲第一区二区三区不卡| 国产极品天堂在线| 大陆偷拍与自拍| 国产黄片美女视频| 日韩,欧美,国产一区二区三区| 日本欧美视频一区| 亚洲欧美日韩东京热| 久热这里只有精品99| 香蕉精品网在线| 久久午夜福利片| 春色校园在线视频观看| 国产 精品1| 香蕉精品网在线| 又爽又黄a免费视频| 人妻 亚洲 视频| 伦理电影大哥的女人| 你懂的网址亚洲精品在线观看| 欧美精品人与动牲交sv欧美| av免费观看日本| 亚洲怡红院男人天堂| 国产成人91sexporn| 国产精品免费大片| 天天躁日日操中文字幕| 国产精品嫩草影院av在线观看| 国产视频首页在线观看| 国产精品欧美亚洲77777| 精品国产乱码久久久久久小说| 久久午夜福利片| 免费大片18禁| 天天躁日日操中文字幕| 国产亚洲精品久久久com| a级毛色黄片| 亚洲不卡免费看| 大又大粗又爽又黄少妇毛片口| 麻豆成人av视频| 精品久久久噜噜| 亚洲av成人精品一区久久| 亚洲性久久影院| 国产成人一区二区在线| 欧美激情极品国产一区二区三区 | 亚洲欧美日韩东京热| 午夜福利网站1000一区二区三区| 精品人妻视频免费看| 草草在线视频免费看| 九九久久精品国产亚洲av麻豆| 99久久精品热视频| 大片电影免费在线观看免费| 久久韩国三级中文字幕| 在线免费观看不下载黄p国产| 99久久人妻综合| 国产免费又黄又爽又色| 欧美高清成人免费视频www| 国产亚洲av片在线观看秒播厂| 精品99又大又爽又粗少妇毛片| 国产久久久一区二区三区| 精品视频人人做人人爽| 久久精品国产亚洲av涩爱| 在线观看av片永久免费下载| 中国三级夫妇交换| 99久久中文字幕三级久久日本| 一本色道久久久久久精品综合| 欧美日韩视频高清一区二区三区二| 五月天丁香电影| 免费黄网站久久成人精品| 国产精品不卡视频一区二区| 永久网站在线| 成人黄色视频免费在线看| 国产精品三级大全| 一区二区三区免费毛片| 免费在线观看成人毛片| 成年女人在线观看亚洲视频| 黄色日韩在线| 日韩伦理黄色片| 国产伦精品一区二区三区视频9| 国产精品女同一区二区软件| 成年人午夜在线观看视频| av在线老鸭窝| 五月开心婷婷网| 亚洲精品日韩在线中文字幕| 日韩欧美 国产精品| 寂寞人妻少妇视频99o| 22中文网久久字幕| 只有这里有精品99| 亚洲精品自拍成人| 日韩精品有码人妻一区| 精品久久久久久电影网| 我要看日韩黄色一级片| 免费大片18禁| 国内精品宾馆在线| 国产人妻一区二区三区在| 久热久热在线精品观看| 大又大粗又爽又黄少妇毛片口| 欧美激情国产日韩精品一区| 成人毛片a级毛片在线播放| 日韩不卡一区二区三区视频在线| 亚洲国产精品一区三区| 99热6这里只有精品| 精品一区二区三卡| 男女免费视频国产| 日韩中文字幕视频在线看片 | 午夜视频国产福利| 我的老师免费观看完整版| 一本—道久久a久久精品蜜桃钙片| 久久韩国三级中文字幕| 久久久久久人妻| 人妻制服诱惑在线中文字幕| 国内少妇人妻偷人精品xxx网站| 亚洲精品乱久久久久久| 天天躁日日操中文字幕| 王馨瑶露胸无遮挡在线观看| 麻豆精品久久久久久蜜桃| 国产av一区二区精品久久 | 人妻制服诱惑在线中文字幕| 18+在线观看网站| 久久鲁丝午夜福利片| 亚洲欧美日韩东京热| 亚洲欧美日韩东京热| 久久久久久九九精品二区国产| 国产伦精品一区二区三区视频9| 美女xxoo啪啪120秒动态图| av女优亚洲男人天堂| 乱码一卡2卡4卡精品| 伦精品一区二区三区| 在线观看av片永久免费下载| 寂寞人妻少妇视频99o| 久久久欧美国产精品| 插逼视频在线观看| 久久av网站| 赤兔流量卡办理| 久久国产精品大桥未久av | 亚洲电影在线观看av| av在线观看视频网站免费| 久久久精品94久久精品| 久久久久久久久久人人人人人人| 99久久精品国产国产毛片| 男男h啪啪无遮挡| 亚洲四区av| 日韩一区二区三区影片| 成人美女网站在线观看视频| 91精品一卡2卡3卡4卡| 在线观看av片永久免费下载| 午夜福利在线观看免费完整高清在| 成人毛片60女人毛片免费| 九九久久精品国产亚洲av麻豆| 好男人视频免费观看在线| 老师上课跳d突然被开到最大视频| 国产成人午夜福利电影在线观看| 欧美高清性xxxxhd video| 国产69精品久久久久777片| 亚洲性久久影院| 久久精品久久久久久噜噜老黄| 校园人妻丝袜中文字幕| 久久精品国产鲁丝片午夜精品| 国产精品一区二区在线观看99| 26uuu在线亚洲综合色| 亚洲欧美日韩卡通动漫| 国产深夜福利视频在线观看| 国产熟女欧美一区二区| 精品久久国产蜜桃| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 亚洲真实伦在线观看| 在线观看国产h片| 人妻少妇偷人精品九色| 亚洲精华国产精华液的使用体验| 亚洲美女黄色视频免费看| 久久久精品94久久精品| 自拍欧美九色日韩亚洲蝌蚪91 | 99久久精品热视频| 亚洲四区av| 91精品国产九色| 亚洲电影在线观看av| 国产真实伦视频高清在线观看| 午夜激情久久久久久久| 国产精品久久久久久av不卡| 久久久久国产精品人妻一区二区| 国产 一区 欧美 日韩| 高清视频免费观看一区二区| 日韩人妻高清精品专区| 日韩中文字幕视频在线看片 | 久久久久久久亚洲中文字幕| 最近中文字幕2019免费版| 日韩大片免费观看网站| 国产亚洲精品久久久com| 国产在线男女| av国产免费在线观看| av卡一久久| 免费少妇av软件| 日韩欧美 国产精品| 国产国拍精品亚洲av在线观看| 韩国高清视频一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 精品国产一区二区三区久久久樱花 | 寂寞人妻少妇视频99o| 男女下面进入的视频免费午夜| av国产精品久久久久影院| 中文字幕久久专区| 中文欧美无线码| 超碰av人人做人人爽久久| 国产 精品1| 在线观看人妻少妇| 内射极品少妇av片p| 在线观看免费日韩欧美大片 | 国产乱人偷精品视频| av福利片在线观看| 久久久久久久久大av| 男人和女人高潮做爰伦理| h视频一区二区三区| 亚洲一区二区三区欧美精品| 美女xxoo啪啪120秒动态图| 嘟嘟电影网在线观看| 国产男女超爽视频在线观看| 欧美bdsm另类| 免费观看无遮挡的男女| 精品亚洲乱码少妇综合久久| 一个人看的www免费观看视频| 久久这里有精品视频免费| 成年美女黄网站色视频大全免费 | 国产爽快片一区二区三区| 高清午夜精品一区二区三区| 色哟哟·www| 在线亚洲精品国产二区图片欧美 | av卡一久久| 美女中出高潮动态图| 91aial.com中文字幕在线观看| 精品久久久精品久久久| 熟女av电影| 亚洲精品乱久久久久久| 我要看黄色一级片免费的| 国产91av在线免费观看| 国产av一区二区精品久久 | 国产精品国产三级国产av玫瑰| 国产成人a∨麻豆精品| 日韩视频在线欧美| 欧美亚洲 丝袜 人妻 在线| 毛片一级片免费看久久久久| 国产淫语在线视频| 色网站视频免费| 免费av中文字幕在线| 亚洲成人一二三区av| 国产男女超爽视频在线观看| 亚洲精品国产成人久久av| 伊人久久国产一区二区| 高清黄色对白视频在线免费看 | 成人毛片a级毛片在线播放| 97在线视频观看| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 精品久久国产蜜桃| 18禁动态无遮挡网站| 亚洲av中文av极速乱| 日韩亚洲欧美综合| 99re6热这里在线精品视频| xxx大片免费视频| 亚洲成人一二三区av| 肉色欧美久久久久久久蜜桃| 狂野欧美激情性xxxx在线观看| 久久人人爽人人片av| 高清视频免费观看一区二区| 国产精品秋霞免费鲁丝片| 亚洲欧美中文字幕日韩二区| 免费大片黄手机在线观看| 亚洲真实伦在线观看| 亚洲欧洲国产日韩| 亚洲aⅴ乱码一区二区在线播放| 少妇人妻精品综合一区二区| 精品久久久噜噜| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| 亚洲精品国产av蜜桃| 成人亚洲精品一区在线观看 | 日韩欧美精品免费久久| 久久久久久久久久久丰满| 美女脱内裤让男人舔精品视频| 精品久久久久久久久av| 大香蕉97超碰在线| 91狼人影院| 久久精品国产自在天天线| 在线观看免费高清a一片| 99久久人妻综合| 男女啪啪激烈高潮av片| 在线免费十八禁| 成人无遮挡网站| 内射极品少妇av片p| 亚洲真实伦在线观看| 一级毛片我不卡| 久久av网站| 国产精品人妻久久久影院| 亚洲成人一二三区av| 亚洲av日韩在线播放| 国精品久久久久久国模美| 精品久久久噜噜| 精华霜和精华液先用哪个| 夜夜骑夜夜射夜夜干| 国产成人午夜福利电影在线观看| 国内揄拍国产精品人妻在线| 久久久久网色| 蜜桃在线观看..| 久久久久久伊人网av| 国产乱来视频区| 天堂中文最新版在线下载| 在线观看美女被高潮喷水网站| 国产日韩欧美在线精品| 日本色播在线视频| 国产片特级美女逼逼视频| 国产精品一区二区三区四区免费观看| 男女边吃奶边做爰视频| 亚洲人成网站在线播| 欧美xxⅹ黑人| 91在线精品国自产拍蜜月| 久久婷婷青草| 亚洲国产精品专区欧美| 青春草国产在线视频| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 久久久久性生活片| .国产精品久久| 国产精品女同一区二区软件| 尤物成人国产欧美一区二区三区| 全区人妻精品视频| 一级毛片电影观看| 韩国av在线不卡| 九草在线视频观看| 色吧在线观看| 男人爽女人下面视频在线观看| 久久婷婷青草| av播播在线观看一区| 高清不卡的av网站| 极品少妇高潮喷水抽搐| 国产av码专区亚洲av| 欧美精品国产亚洲| 男人爽女人下面视频在线观看| 在线观看免费视频网站a站| 日韩欧美 国产精品| 国产精品三级大全| 一级黄片播放器| 日韩一区二区视频免费看| 这个男人来自地球电影免费观看 | 91精品国产国语对白视频| 久久精品国产a三级三级三级| 欧美日韩一区二区视频在线观看视频在线| 五月玫瑰六月丁香| 亚洲欧美日韩另类电影网站 | 欧美人与善性xxx| 欧美bdsm另类| 在现免费观看毛片| 网址你懂的国产日韩在线| 人人妻人人添人人爽欧美一区卜 | 日本-黄色视频高清免费观看| 精品久久久久久久末码| 性高湖久久久久久久久免费观看| 一区二区三区四区激情视频| 日日摸夜夜添夜夜添av毛片| 丰满乱子伦码专区| 亚洲丝袜综合中文字幕| 啦啦啦在线观看免费高清www| tube8黄色片| 中国三级夫妇交换| 中文在线观看免费www的网站| 最近最新中文字幕免费大全7| 性色avwww在线观看| 看十八女毛片水多多多| 嫩草影院新地址| 国产欧美日韩精品一区二区| 精品国产乱码久久久久久小说| 欧美丝袜亚洲另类| 伦精品一区二区三区| 精品午夜福利在线看| 国产成人a区在线观看| 婷婷色综合大香蕉| 简卡轻食公司| 国产免费一级a男人的天堂| 99久久综合免费| 国产精品不卡视频一区二区| 久久久国产一区二区| 国产日韩欧美亚洲二区| 亚洲熟女精品中文字幕| 国产成人aa在线观看| 日本与韩国留学比较| 女性被躁到高潮视频| 少妇高潮的动态图| 日本爱情动作片www.在线观看| kizo精华| 国产亚洲午夜精品一区二区久久| 直男gayav资源| 亚洲av成人精品一二三区| 国产v大片淫在线免费观看| 最黄视频免费看| av线在线观看网站| 亚洲综合色惰| 亚洲av成人精品一区久久| 黄片wwwwww| av黄色大香蕉| 国产欧美另类精品又又久久亚洲欧美| 国产高清三级在线| 国产精品国产三级国产av玫瑰| 乱码一卡2卡4卡精品| 七月丁香在线播放| 国产亚洲av片在线观看秒播厂| 18禁裸乳无遮挡动漫免费视频| 精品亚洲成国产av| 99久久精品热视频| 国产大屁股一区二区在线视频| 日日撸夜夜添| 亚洲人成网站高清观看| 777米奇影视久久| 一级毛片aaaaaa免费看小| 久久久国产一区二区| 久久婷婷青草| 搡老乐熟女国产| 国产伦在线观看视频一区| 综合色丁香网| 国产成人freesex在线| 免费观看无遮挡的男女| 久久女婷五月综合色啪小说| 国产精品久久久久久av不卡| 久久鲁丝午夜福利片| 高清不卡的av网站| 精品国产一区二区三区久久久樱花 | 一级毛片 在线播放| 欧美高清性xxxxhd video| 少妇人妻 视频| 尤物成人国产欧美一区二区三区| 亚洲综合色惰| tube8黄色片| 啦啦啦在线观看免费高清www| 国产淫语在线视频| 精品一区在线观看国产| 日日摸夜夜添夜夜添av毛片| 成人18禁高潮啪啪吃奶动态图 | 国产女主播在线喷水免费视频网站| 亚洲国产色片| 亚洲av国产av综合av卡| 干丝袜人妻中文字幕| av国产免费在线观看|