• <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ù)殘差法在有限元模型修正中的應用
    頻響阻抗法診斷變壓器繞組變形
    男人舔女人的私密视频| 亚洲七黄色美女视频| 男人舔女人的私密视频| 又黄又粗又硬又大视频| 丰满人妻熟妇乱又伦精品不卡| 1024视频免费在线观看| 亚洲久久久国产精品| 精品一区二区三卡| 欧美亚洲日本最大视频资源| 久久热在线av| 国产欧美日韩综合在线一区二区| 亚洲性夜色夜夜综合| 999久久久精品免费观看国产| 又黄又粗又硬又大视频| xxx96com| 免费不卡黄色视频| 99精品久久久久人妻精品| 国产国语露脸激情在线看| 亚洲人成77777在线视频| 午夜福利免费观看在线| 日韩欧美国产一区二区入口| 国产欧美日韩一区二区三区在线| 日韩免费av在线播放| 亚洲av片天天在线观看| 精品欧美一区二区三区在线| 他把我摸到了高潮在线观看| 大码成人一级视频| 成人免费观看视频高清| 又黄又爽又免费观看的视频| 夫妻午夜视频| 香蕉国产在线看| 日韩av在线大香蕉| 亚洲人成77777在线视频| 岛国视频午夜一区免费看| 一二三四在线观看免费中文在| 丁香欧美五月| 丝袜在线中文字幕| 国产精品99久久99久久久不卡| 欧美乱色亚洲激情| 91精品国产国语对白视频| 在线播放国产精品三级| 久久久久国产一级毛片高清牌| 亚洲欧美激情综合另类| 成年人黄色毛片网站| 午夜福利在线观看吧| 国产亚洲精品久久久久5区| 国产一区二区在线av高清观看| 久久中文看片网| 亚洲国产看品久久| 91字幕亚洲| 亚洲精品在线美女| 色婷婷久久久亚洲欧美| 精品卡一卡二卡四卡免费| 亚洲欧美日韩高清在线视频| 在线播放国产精品三级| 亚洲精品国产一区二区精华液| 久久久久久亚洲精品国产蜜桃av| 美女 人体艺术 gogo| 中文欧美无线码| 亚洲五月婷婷丁香| 亚洲精华国产精华精| 亚洲av熟女| 天天添夜夜摸| 日韩欧美三级三区| 午夜福利一区二区在线看| 成人三级黄色视频| 亚洲成人免费av在线播放| 亚洲精品中文字幕在线视频| 亚洲熟女毛片儿| 一级a爱视频在线免费观看| 欧美性长视频在线观看| 欧美 亚洲 国产 日韩一| 在线观看免费高清a一片| 黄片播放在线免费| 欧洲精品卡2卡3卡4卡5卡区| 首页视频小说图片口味搜索| 亚洲欧美激情综合另类| 真人一进一出gif抽搐免费| 国产精品偷伦视频观看了| 中亚洲国语对白在线视频| 中文欧美无线码| 真人一进一出gif抽搐免费| 男女下面进入的视频免费午夜 | 成人国语在线视频| 亚洲aⅴ乱码一区二区在线播放 | 老汉色av国产亚洲站长工具| 宅男免费午夜| 日韩欧美一区视频在线观看| 91字幕亚洲| 在线观看www视频免费| 50天的宝宝边吃奶边哭怎么回事| 色尼玛亚洲综合影院| 欧美黄色淫秽网站| 亚洲成a人片在线一区二区| 国产亚洲欧美98| 99久久综合精品五月天人人| 精品国产一区二区久久| 操出白浆在线播放| 香蕉国产在线看| 十八禁人妻一区二区| 国产伦人伦偷精品视频| 69av精品久久久久久| 交换朋友夫妻互换小说| 极品教师在线免费播放| 亚洲av第一区精品v没综合| 国产精品一区二区三区四区久久 | 婷婷六月久久综合丁香| 在线观看免费视频网站a站| av视频免费观看在线观看| 亚洲第一青青草原| 91麻豆精品激情在线观看国产 | 中出人妻视频一区二区| 精品午夜福利视频在线观看一区| 亚洲欧美激情综合另类| 热99re8久久精品国产| 99热国产这里只有精品6| tocl精华| 69精品国产乱码久久久| 男人操女人黄网站| 如日韩欧美国产精品一区二区三区| av超薄肉色丝袜交足视频| 中文字幕另类日韩欧美亚洲嫩草| 国产一区二区三区在线臀色熟女 | 国产亚洲欧美精品永久| 人妻久久中文字幕网| 国产色视频综合| 亚洲人成电影观看| 18美女黄网站色大片免费观看| 国产成人一区二区三区免费视频网站| 一级作爱视频免费观看| 男女下面插进去视频免费观看| 母亲3免费完整高清在线观看| 国产av精品麻豆| 美女高潮到喷水免费观看| 天天添夜夜摸| 亚洲精品中文字幕一二三四区| а√天堂www在线а√下载| 亚洲七黄色美女视频| 老鸭窝网址在线观看| 日韩有码中文字幕| 如日韩欧美国产精品一区二区三区| 欧美日韩视频精品一区| 久久九九热精品免费| 大型av网站在线播放| 久久人妻熟女aⅴ| 视频在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 视频区图区小说| 一进一出抽搐动态| 女人精品久久久久毛片| 久久精品国产99精品国产亚洲性色 | 亚洲成人免费电影在线观看| 999久久久国产精品视频| 三上悠亚av全集在线观看| 自线自在国产av| 婷婷精品国产亚洲av在线| 波多野结衣高清无吗| 久久久久亚洲av毛片大全| 男女午夜视频在线观看| 国产在线观看jvid| 亚洲精品久久成人aⅴ小说| 丰满的人妻完整版| 黄色女人牲交| 亚洲av美国av| 国产一区二区在线av高清观看| 国产精品av久久久久免费| 好看av亚洲va欧美ⅴa在| 夜夜躁狠狠躁天天躁| 亚洲情色 制服丝袜| 99精品在免费线老司机午夜| 丝袜美足系列| 日日夜夜操网爽| 国产精品一区二区免费欧美| 精品久久久久久成人av| 在线播放国产精品三级| 欧美中文日本在线观看视频| 日韩欧美在线二视频| 久久草成人影院| 法律面前人人平等表现在哪些方面| 久久久久久久久免费视频了| 一级毛片高清免费大全| 亚洲av五月六月丁香网| 国产一区二区三区视频了| 国产三级在线视频| 久久精品aⅴ一区二区三区四区| 欧美黄色片欧美黄色片| 香蕉丝袜av| 在线观看日韩欧美| 少妇的丰满在线观看| 精品欧美一区二区三区在线| 亚洲欧美精品综合久久99| 亚洲一卡2卡3卡4卡5卡精品中文| 日日干狠狠操夜夜爽| 国产色视频综合| 日韩一卡2卡3卡4卡2021年| 午夜久久久在线观看| 日韩免费高清中文字幕av| 亚洲精品国产区一区二| 国产一区在线观看成人免费| 亚洲成av片中文字幕在线观看| 亚洲九九香蕉| 国产精品自产拍在线观看55亚洲| 亚洲av成人av| 国产一区二区在线av高清观看| 黄片小视频在线播放| 一级黄色大片毛片| 亚洲 欧美一区二区三区| 国产在线精品亚洲第一网站| 日日摸夜夜添夜夜添小说| 手机成人av网站| 国产成人欧美| 国产真人三级小视频在线观看| 亚洲成av片中文字幕在线观看| 久久久久久大精品| 日韩视频一区二区在线观看| 亚洲五月天丁香| 国产精品美女特级片免费视频播放器 | 欧美黄色片欧美黄色片| 在线天堂中文资源库| 久久久国产精品麻豆| 国产伦人伦偷精品视频| av超薄肉色丝袜交足视频| 在线观看日韩欧美| 亚洲男人的天堂狠狠| 久久青草综合色| 亚洲av成人一区二区三| 午夜日韩欧美国产| 欧美日韩av久久| 国产亚洲精品久久久久久毛片| 97碰自拍视频| 美女大奶头视频| 亚洲久久久国产精品| 丰满的人妻完整版| 午夜福利在线免费观看网站| 视频区欧美日本亚洲| 日韩视频一区二区在线观看| av网站免费在线观看视频| 国产精品av久久久久免费| 在线看a的网站| 亚洲黑人精品在线| 两性夫妻黄色片| 成年女人毛片免费观看观看9| 国产成人av激情在线播放| 脱女人内裤的视频| 一级a爱视频在线免费观看| 成熟少妇高潮喷水视频| 黄色怎么调成土黄色| 亚洲成人免费电影在线观看| 成人av一区二区三区在线看| 国产又色又爽无遮挡免费看| av有码第一页| 亚洲aⅴ乱码一区二区在线播放 | av有码第一页| 18美女黄网站色大片免费观看| 日韩高清综合在线| 久久久久国产一级毛片高清牌| 91老司机精品| 色老头精品视频在线观看| 国产成人欧美| 90打野战视频偷拍视频| 老鸭窝网址在线观看| 中文字幕人妻熟女乱码| 黑人猛操日本美女一级片| 一级毛片高清免费大全| 在线观看66精品国产| 国产亚洲欧美98| 日韩一卡2卡3卡4卡2021年| 丰满人妻熟妇乱又伦精品不卡| 99久久99久久久精品蜜桃| 亚洲欧美激情在线| 亚洲欧美日韩无卡精品| 日本精品一区二区三区蜜桃| 午夜亚洲福利在线播放| 久久久久国内视频| 女人精品久久久久毛片| 我的亚洲天堂| 欧美精品啪啪一区二区三区| 久久婷婷成人综合色麻豆| 国产成人系列免费观看| 天堂√8在线中文| 电影成人av| 国产精品秋霞免费鲁丝片| 久久影院123| 久久久久国产精品人妻aⅴ院| 日韩精品青青久久久久久| 亚洲av成人av| 黄片大片在线免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 日日爽夜夜爽网站| 999久久久精品免费观看国产| 欧美黑人精品巨大| 久久精品国产亚洲av香蕉五月| 天天影视国产精品| 国产在线观看jvid| 村上凉子中文字幕在线| www.自偷自拍.com| 老熟妇仑乱视频hdxx| 悠悠久久av| 亚洲美女黄片视频| 国产三级在线视频| 亚洲少妇的诱惑av| 亚洲伊人色综图| 国产亚洲精品第一综合不卡| 亚洲少妇的诱惑av| 一级a爱片免费观看的视频| 日韩 欧美 亚洲 中文字幕| 精品人妻在线不人妻| 免费av毛片视频| 一个人观看的视频www高清免费观看 | 高清黄色对白视频在线免费看| x7x7x7水蜜桃| 真人做人爱边吃奶动态| 国产一区二区三区视频了| 日韩欧美一区二区三区在线观看| 久久中文字幕人妻熟女| 免费在线观看黄色视频的| 淫妇啪啪啪对白视频| 黑人巨大精品欧美一区二区mp4| 午夜a级毛片| 精品免费久久久久久久清纯| 很黄的视频免费| 国产99久久九九免费精品| 丰满饥渴人妻一区二区三| 国产精品久久久人人做人人爽| 亚洲欧美一区二区三区黑人| 日韩一卡2卡3卡4卡2021年| 中文字幕av电影在线播放| 亚洲精品成人av观看孕妇| 欧美黑人精品巨大| 男女之事视频高清在线观看| 免费在线观看亚洲国产| av国产精品久久久久影院| 无遮挡黄片免费观看| 叶爱在线成人免费视频播放| 成人国语在线视频| 777久久人妻少妇嫩草av网站| 免费少妇av软件| 最近最新中文字幕大全免费视频| 婷婷丁香在线五月| 777久久人妻少妇嫩草av网站| 久久亚洲精品不卡| 亚洲一区二区三区欧美精品| 国产人伦9x9x在线观看| 午夜两性在线视频| 日本精品一区二区三区蜜桃| 久久精品亚洲av国产电影网| 制服人妻中文乱码| 97人妻天天添夜夜摸| 另类亚洲欧美激情| 免费看十八禁软件| 久久人妻av系列| 不卡一级毛片| 亚洲精品美女久久av网站| 久久久国产成人免费| 欧美最黄视频在线播放免费 | 欧美一区二区精品小视频在线| 亚洲三区欧美一区| 午夜视频精品福利| 久久天堂一区二区三区四区| 中亚洲国语对白在线视频| 中文字幕最新亚洲高清| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区三区在线| 日韩欧美三级三区| 99热国产这里只有精品6| 免费一级毛片在线播放高清视频 | 精品福利观看| 国产av在哪里看| 亚洲国产精品一区二区三区在线| 国产成人啪精品午夜网站| 精品国产亚洲在线| ponron亚洲| 国产黄a三级三级三级人| 视频在线观看一区二区三区| √禁漫天堂资源中文www| 国产精品98久久久久久宅男小说| 少妇被粗大的猛进出69影院| 一级片免费观看大全| 国产成人av激情在线播放| 久久久久亚洲av毛片大全| 久久人妻福利社区极品人妻图片| 精品一区二区三区四区五区乱码| 色播在线永久视频| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 人妻丰满熟妇av一区二区三区| 成人特级黄色片久久久久久久| 他把我摸到了高潮在线观看| 久久性视频一级片| 夜夜躁狠狠躁天天躁| 欧美一区二区精品小视频在线| videosex国产| 在线av久久热| 免费日韩欧美在线观看| 国产精品国产av在线观看| 超色免费av| 国产成人免费无遮挡视频| 亚洲成人精品中文字幕电影 | 大香蕉久久成人网| 国产精品野战在线观看 | 18禁裸乳无遮挡免费网站照片 | 欧美不卡视频在线免费观看 | 一进一出好大好爽视频| 久久久久九九精品影院| 99精品在免费线老司机午夜| 国产av在哪里看| 中文字幕人妻丝袜制服| 免费看十八禁软件| 国产av一区在线观看免费| 激情视频va一区二区三区| 国产欧美日韩精品亚洲av| 麻豆av在线久日| 国产激情久久老熟女| 久久人人97超碰香蕉20202| 亚洲美女黄片视频| 香蕉丝袜av| 久久久国产成人精品二区 | 91精品三级在线观看| 精品久久蜜臀av无| 精品国产亚洲在线| 美女高潮喷水抽搐中文字幕| 欧美色视频一区免费| 黑人操中国人逼视频| 色在线成人网| 亚洲aⅴ乱码一区二区在线播放 | 久久草成人影院| 搡老熟女国产l中国老女人| 欧美av亚洲av综合av国产av| 亚洲国产欧美一区二区综合| 丰满人妻熟妇乱又伦精品不卡| 麻豆成人av在线观看| 91精品三级在线观看| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 欧美午夜高清在线| 视频区图区小说| 手机成人av网站| 19禁男女啪啪无遮挡网站| 国产亚洲欧美98| 国产乱人伦免费视频| 久久久水蜜桃国产精品网| www.精华液| 亚洲少妇的诱惑av| 亚洲五月婷婷丁香| 丁香六月欧美| 巨乳人妻的诱惑在线观看| 国产精华一区二区三区| 99精品在免费线老司机午夜| 亚洲欧美日韩无卡精品| 麻豆av在线久日| 嫁个100分男人电影在线观看| 亚洲精品在线美女| 淫秽高清视频在线观看| 亚洲精品国产色婷婷电影| 在线视频色国产色| 欧美日韩视频精品一区| 午夜福利影视在线免费观看| 一个人免费在线观看的高清视频| 久久热在线av| 欧美+亚洲+日韩+国产| 中国美女看黄片| 国产精品久久久久成人av| 国产精品爽爽va在线观看网站 | 夜夜看夜夜爽夜夜摸 | 夜夜夜夜夜久久久久| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 看免费av毛片| 精品人妻在线不人妻| 两人在一起打扑克的视频| 国产高清视频在线播放一区| 在线看a的网站| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区久久 | 黑人操中国人逼视频| 99re在线观看精品视频| 久久天躁狠狠躁夜夜2o2o| 国产av一区在线观看免费| 精品国产亚洲在线| av在线天堂中文字幕 | 91老司机精品| 亚洲少妇的诱惑av| 两性夫妻黄色片| 亚洲av片天天在线观看| 黄色女人牲交| 久久久久久亚洲精品国产蜜桃av| 美女国产高潮福利片在线看| 88av欧美| 精品久久久久久电影网| 久久国产精品影院| 丝袜美腿诱惑在线| 亚洲欧美日韩高清在线视频| 成人三级黄色视频| av福利片在线| 久热爱精品视频在线9| av在线天堂中文字幕 | 精品高清国产在线一区| 亚洲久久久国产精品| 国产激情久久老熟女| 在线观看www视频免费| 国产成人精品久久二区二区91| 免费在线观看日本一区| 午夜a级毛片| 亚洲人成伊人成综合网2020| 久热爱精品视频在线9| 天天添夜夜摸| 欧美久久黑人一区二区| 国产成人系列免费观看| 99香蕉大伊视频| 女同久久另类99精品国产91| 国产亚洲精品久久久久5区| 一级毛片女人18水好多| 久久亚洲精品不卡| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 男女之事视频高清在线观看| 丁香六月欧美| 99精品欧美一区二区三区四区| 两个人看的免费小视频| 国产主播在线观看一区二区| 精品久久久久久电影网| 国产99白浆流出| 日日干狠狠操夜夜爽| 精品卡一卡二卡四卡免费| 欧美激情极品国产一区二区三区| 久久人妻av系列| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| 亚洲美女黄片视频| 日韩一卡2卡3卡4卡2021年| 日韩免费av在线播放| 波多野结衣一区麻豆| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 成年版毛片免费区| 电影成人av| av电影中文网址| 一个人免费在线观看的高清视频| 中文字幕av电影在线播放| 在线观看一区二区三区| 成年女人毛片免费观看观看9| 国产三级在线视频| 老司机午夜十八禁免费视频| 嫩草影视91久久| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 成人手机av| 国产亚洲精品久久久久久毛片| 久久天躁狠狠躁夜夜2o2o| 欧美日韩瑟瑟在线播放| 欧美日韩乱码在线| 精品久久久精品久久久| 91av网站免费观看| 9191精品国产免费久久| 琪琪午夜伦伦电影理论片6080| 国产精品电影一区二区三区| 亚洲久久久国产精品| 精品久久久久久久久久免费视频 | 国产伦人伦偷精品视频| 热re99久久国产66热| 欧美激情极品国产一区二区三区| 久久久久久久久免费视频了| 亚洲激情在线av| 制服人妻中文乱码| 丰满的人妻完整版| 午夜成年电影在线免费观看| 男人舔女人下体高潮全视频| 99久久综合精品五月天人人| 亚洲精品一区av在线观看| 女人精品久久久久毛片| 男女午夜视频在线观看| av片东京热男人的天堂| 韩国精品一区二区三区| 国产极品粉嫩免费观看在线| 黄色 视频免费看| 真人一进一出gif抽搐免费| 9热在线视频观看99| 国产精品久久久久久人妻精品电影| 无限看片的www在线观看| 午夜亚洲福利在线播放| 老司机亚洲免费影院| 人成视频在线观看免费观看| 欧美日韩瑟瑟在线播放| 欧美成人性av电影在线观看| a级毛片黄视频| av网站免费在线观看视频| 欧美日韩精品网址| 视频在线观看一区二区三区| 久久香蕉精品热| 十八禁人妻一区二区| 妹子高潮喷水视频| 午夜精品在线福利| 国产高清视频在线播放一区| 一进一出抽搐动态| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 精品国产一区二区三区四区第35| 99久久综合精品五月天人人| 久久狼人影院| 多毛熟女@视频| 亚洲精品一二三| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 久久影院123| 亚洲男人天堂网一区| 免费观看精品视频网站| 91精品三级在线观看| 欧美最黄视频在线播放免费 | 99国产精品一区二区蜜桃av| 黄色成人免费大全| 国产成人影院久久av| 999久久久精品免费观看国产| 久9热在线精品视频| 99久久人妻综合| 一级片'在线观看视频|