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

    軟基水閘底板脫空反演的響應(yīng)面模型優(yōu)化

    2022-01-07 05:59:50劉雙平柯賢勇李火坤李輝唐義員涂源方靜
    關(guān)鍵詞:水閘反演底板

    劉雙平,柯賢勇,李火坤,李輝,唐義員,涂源,方靜

    (南昌大學(xué)建筑工程學(xué)院,江西 南昌 330031)

    長期處于激振交替、動靜懸殊的工作環(huán)境中的軟基水閘容易出現(xiàn)底板脫空[1],目前,基于動力學(xué)特性對水閘結(jié)構(gòu)進(jìn)行反分析被廣泛研究[2-4]。尋找準(zhǔn)確可靠的響應(yīng)面模型代替有限元模型,是進(jìn)行軟基水閘底板脫空反演的關(guān)鍵。

    響應(yīng)面方程是決定脫空參數(shù)與模態(tài)參數(shù)之間非線性關(guān)系的重要因素,當(dāng)響應(yīng)面方程過于簡單時,響應(yīng)面模型難以保證其精度;當(dāng)響應(yīng)面方程過于復(fù)雜時,響應(yīng)面模型的效率會降低[5]。于景飛等[6]采用線性多項式函數(shù)構(gòu)建響應(yīng)面方程,對大跨徑混凝土斜拉橋進(jìn)行了優(yōu)化預(yù)測研究;Fang等[7]指出參數(shù)間相互效應(yīng)對響應(yīng)面模型總方差的貢獻(xiàn)非常小;張運濤等[8]基于逐步回歸的方法對變量的選取進(jìn)行顯著性分析,提高了所選取不含交叉項二次多項式響應(yīng)面模型的精度和計算效率,并將其運用于大跨連續(xù)剛構(gòu)橋長期變形預(yù)測;李延強(qiáng)等[9]采用F檢驗法進(jìn)行參數(shù)的顯著性分析,以二次多項式擬合響應(yīng)面模型,對比了考慮交叉項和不考慮交叉項對斜拉橋模型的修正結(jié)果;宗周紅、劉杰等[10-11]采用高階響應(yīng)面模型應(yīng)用于橋梁結(jié)構(gòu)的損傷識別。張建偉等[12]采用基于方差分析的F檢驗法,對響應(yīng)面模型參數(shù)進(jìn)行顯著性檢驗,并將優(yōu)化后的響應(yīng)面模型應(yīng)用于某大型渡槽結(jié)構(gòu)的有限元模型修正;蒙偉等[13]分析了F檢驗和t檢驗的特點,采用了一種能直觀方便地判斷回歸模型和回歸系數(shù)顯著性的p值法檢驗,并將其應(yīng)用于初始地應(yīng)力場的反演分析。

    為獲得軟基水閘底板脫空反演中最佳的響應(yīng)面模型,本文提出了一種基于顯著性檢驗的響應(yīng)面模型優(yōu)選與優(yōu)化方法。首先,采用拉丁超立方抽樣(LHS)獲得脫空參數(shù)樣本集,并將其輸入水閘有限元模型計算相應(yīng)的模態(tài)參數(shù)集,以此構(gòu)建表征脫空參數(shù)與模態(tài)參數(shù)之間非線性關(guān)系的多種不同形式的響應(yīng)面模型;其次,對不同形式的響應(yīng)面模型進(jìn)行F檢驗,綜合考慮精度和效率,優(yōu)選形式最佳的響應(yīng)面模型;然后,對形式最佳的響應(yīng)面模型進(jìn)行t檢驗,優(yōu)化響應(yīng)面模型的參數(shù),剔除對模態(tài)參數(shù)影響不足的脫空參數(shù)及其組合(贅余項)以獲得最佳的響應(yīng)面模型;最后,將最佳的響應(yīng)面模型應(yīng)用于軟基水閘物理模型的底板脫空反演,驗證方法的有效性。

    1 響應(yīng)面模型優(yōu)選與優(yōu)化

    1.1 拉丁超立方抽樣試驗設(shè)計

    試驗樣本的選取方法與數(shù)量關(guān)系到響應(yīng)面模型的精度和效率,以及分析結(jié)果的準(zhǔn)確性[14]。采用拉丁超立方抽樣試驗設(shè)計,能夠以少量的脫空參數(shù)樣本,全面而準(zhǔn)確地反映軟基水閘底板脫空的分布情況。采用LHS生成脫空參數(shù)樣本集思路如下[15-16]:

    然后,根據(jù)脫空參數(shù)的概率分布情況,對每一個子區(qū)間進(jìn)行獨立的等概率抽樣,第j個子區(qū)間抽樣獲取的脫空參數(shù)值dij滿足以下關(guān)系:

    (1)

    式中:dij表示第i個脫空參數(shù)在第j個子區(qū)間的抽樣值;r表示服從[0,1]均勻分布的隨機(jī)數(shù);m表示脫空參數(shù)的個數(shù),即抽樣維度;n表示子區(qū)間數(shù),即抽樣組數(shù)。

    最后,將各子區(qū)間生成的樣本點匯總成樣本集,對各脫空參數(shù)的抽樣值隨機(jī)組合形成多維脫空參數(shù)樣本集。

    1.2 響應(yīng)面模型

    有限元模型結(jié)構(gòu)復(fù)雜,單元龐大,難以進(jìn)行反分析計算。合理構(gòu)建脫空參數(shù)與模態(tài)參數(shù)(固有頻率和節(jié)點振型)之間的響應(yīng)面模型代替有限元模型,是進(jìn)行軟基水閘底板脫空反演的基礎(chǔ)。本文采用的響應(yīng)面方程為基于Taylor展開的多項式函數(shù),響應(yīng)面方程如下:

    (2)

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    式中:α表示頻率響應(yīng)面方程的待定系數(shù)。

    響應(yīng)面模型的構(gòu)建思路:首先,通過LHS獲取脫空參數(shù)樣本集;其次,將脫空參數(shù)樣本集輸入有限元模型中,計算相應(yīng)的模態(tài)參數(shù)集;然后,根據(jù)脫空參數(shù)集和模態(tài)參數(shù)集,以上述7種不同形式的多項式函數(shù)構(gòu)建響應(yīng)面方程組;最后,根據(jù)多元回歸分析確定響應(yīng)面方程組的待定系數(shù),從而建立7種不同形式的響應(yīng)面模型。

    1.3 顯著性檢驗

    為獲得軟基水閘底板脫空反演中最佳的響應(yīng)面模型,采用顯著性檢驗對響應(yīng)面模型進(jìn)行優(yōu)選與優(yōu)化。首先,基于回歸模型的顯著性檢驗(F檢驗),對上述7種不同形式的響應(yīng)面模型進(jìn)行顯著性分析,優(yōu)選形式最佳的響應(yīng)面模型;隨后,基于回歸系數(shù)的顯著性檢驗(t檢驗),對形式最佳的響應(yīng)面模型的回歸系數(shù)進(jìn)行顯著性分析,優(yōu)化響應(yīng)面模型的參數(shù),得到最佳的響應(yīng)面模型。

    1.3.1F檢驗

    H0:αi=αii=αiii=αiiii=αij=αijk=αiij=0

    (10)

    式中:H0表示頻率響應(yīng)面模型的檢驗狀態(tài)。

    利用方差分析法,計算上述7種響應(yīng)面方程的F統(tǒng)計量,公式如下:

    (11)

    (12)

    (13)

    (14)

    1.3.2t檢驗

    H0:αI=0,I∈{i,ii,iii,iiii,ij,ijk,iij}

    (15)

    式中:H0表示頻率響應(yīng)面模型回歸系數(shù)的檢驗狀態(tài)。

    可由以下公式計算響應(yīng)面方程中每一項回歸系數(shù)的TI值:

    (16)

    (17)

    (18)

    (19)

    1.4 響應(yīng)面模型評價指標(biāo)

    在響應(yīng)面方程形式持續(xù)復(fù)雜化的過程中,響應(yīng)面模型的精度必然會一步步提高,與此同時,擬合響應(yīng)面以及進(jìn)行軟基水閘底板脫空反演的效率必然會有所降低。因此,在優(yōu)選形式最佳的響應(yīng)面方程時,需綜合考慮精度和效率兩方面的影響。本文采用3個精度評價指標(biāo)和1個效率評價指標(biāo)來衡量響應(yīng)面模型的優(yōu)劣。

    1.4.1 精度評價指標(biāo)

    (20)

    (21)

    相對均方根誤差RMSE′反映響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的誤差,其表達(dá)式如下:

    (22)

    相對均方根誤差RMSE′絕對值(經(jīng)歸一化之后的振型存在負(fù)值)越小,表明響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的誤差越小,響應(yīng)面模型精度越高。

    以上兩種評價指標(biāo)均表示響應(yīng)面模型的整體精度,而每一組脫空參數(shù)樣本點對應(yīng)的響應(yīng)面模型模態(tài)參數(shù)的輸出值與有限元模型模態(tài)參數(shù)的計算值之間的相對誤差Rsp,能夠有效地評價響應(yīng)面模型點對點的擬合精度,其表達(dá)式如下:

    (23)

    對所有的脫空參數(shù)樣本點對應(yīng)的模態(tài)參數(shù)由式(23)遍歷,設(shè)定閾值為5‰[17]。對LHS的所有樣本點,超出閾值的個數(shù)越少,Rsp max越小,響應(yīng)面模型精度越高。

    1.4.2 效率評價指標(biāo)

    效率是檢驗響應(yīng)面模型的擬合過程以及在代替有限元模型進(jìn)行軟基水閘底板脫空反演過程中耗時是否經(jīng)濟(jì)合理的重要指標(biāo)。考慮到不同計算機(jī)之間的差異性,計算耗時并不適合直接作為效率評價指標(biāo)。在利用多元回歸分析擬合響應(yīng)面方程中,待定系數(shù)越多,擬合效率越低,此外,在進(jìn)行軟基水閘底板脫空反演中,響應(yīng)面模型的參數(shù)越多,反演效率越低,因此,以響應(yīng)面方程展開后的項數(shù)K作為響應(yīng)面模型的效率評價指標(biāo)是合理可行的。響應(yīng)面方程展開后的項數(shù)K越小,響應(yīng)面模型效率越高。

    1.5 響應(yīng)面模型優(yōu)化流程

    基于顯著性檢驗的軟基水閘底板脫空反演中響應(yīng)面模型優(yōu)選與優(yōu)化流程圖如圖1所示。

    圖1 軟基水閘底板脫空反演中響應(yīng)面模型優(yōu)選與優(yōu)化流程圖Fig.1 Flowchart of response surface model selection and optimization in the inversion for the floor void of sluice on the soft foundation

    2 軟基水閘物理模型驗證

    2.1 軟基水閘物理模型基本概況

    以某水閘工程為背景實例,制作一單孔軟基水閘室內(nèi)物理模型[4]。基于水閘物理模型,采用人工掏空的方式對水閘模型底板順?biāo)鞣较蜻M(jìn)行了單側(cè)型和雙側(cè)型隨機(jī)脫空模擬,并測量脫空縱深,如圖2所示。在描述底板脫空范圍的數(shù)學(xué)模型中[3],采用5個參數(shù)描述單側(cè)脫空(即m=5);10個參數(shù)描述雙側(cè)脫空(即m=10)。脫空模擬共分為3種工況:其中,工況一和工況二為單側(cè)脫空,工況三為雙側(cè)脫空。不同工況下的脫空參數(shù)di如表1所示。

    (a) 物理模型脫空

    表1 水閘物理模型底板脫空工況及脫空參數(shù)Tab.1 Void conditions and void parameters of the sluice physical model bottom plate

    2.2 構(gòu)建響應(yīng)面模型

    根據(jù)水閘物理模型尺寸,結(jié)合工程經(jīng)驗,將底板脫空的抽樣范圍確定為上游側(cè)[0,0.65]m,下游側(cè)[0,0.4]m。單側(cè)脫空和雙側(cè)脫空的抽樣維度分別為5和10,采用LHS生成2 000組脫空參數(shù)樣本。將脫空參數(shù)樣本集輸入水閘有限元基準(zhǔn)模型中,計算與之相應(yīng)的模態(tài)參數(shù),提取前四階頻率和振型,分別構(gòu)建脫空參數(shù)與頻率和振型之間的響應(yīng)面模型。

    2.3 優(yōu)選與優(yōu)化響應(yīng)面模型

    2.3.1 優(yōu)選響應(yīng)面模型的形式

    以雙側(cè)脫空,頻率響應(yīng)面模型為例,對響應(yīng)面模型進(jìn)行F檢驗,F(xiàn)統(tǒng)計量及其臨界值如表2所示,可知,7種響應(yīng)面方程的F統(tǒng)計量遠(yuǎn)大于臨界值,表明7種響應(yīng)面模型均能準(zhǔn)確有效地代替有限元模型。此外,第三階響應(yīng)面模型的F統(tǒng)計量比其余三階都小,表明第三階模態(tài)參數(shù)樣本變異性更大,對脫空參數(shù)的敏感性更高。

    表2 7種頻率響應(yīng)面模型的F統(tǒng)計量及其臨界值(雙側(cè)脫空)Tab.2 F statistics of seven frequency response surface models and their critical values (two-sided void)

    表3 7種頻率響應(yīng)面模型的整體精度(雙側(cè)脫空)Tab.3 Overall accuracy of seven frequency response surface models (two-sided void)

    表4 7種頻率響應(yīng)面模型的樣本點精度(雙側(cè)脫空)Tab.4 Sample point accuracy of seven frequency response surface models (two-sided void)

    表5 7種響應(yīng)面模型的效率Tab.5 Efficiency of seven response surface models

    2.3.2 優(yōu)化響應(yīng)面模型的參數(shù)

    (a) 第一階

    根據(jù)圖3可以直觀地看出對頻率響應(yīng)影響不足的脫空參數(shù)及其組合,將各階頻率響應(yīng)面方程中的贅余項剔除,常數(shù)項及保留項所構(gòu)成的響應(yīng)面方程為最佳的響應(yīng)面方程,經(jīng)參數(shù)優(yōu)化的響應(yīng)面模型即為最佳的響應(yīng)面模型。

    2.4 結(jié)果驗證

    為驗證經(jīng)顯著性檢驗的響應(yīng)面模型的優(yōu)越性,將最佳的響應(yīng)面模型應(yīng)用于軟基水閘底板脫空反演,并將反演結(jié)果與參考文獻(xiàn)[4]對比。不同工況下水閘底板脫空控制參數(shù)反演值與實際值如表6所示、識別的脫空面積及趨勢與實際脫空如圖4所示、識別的脫空面積與實際的脫空面積之間的相對誤差如表7所示,同時,為更準(zhǔn)確衡量反演結(jié)果的準(zhǔn)確性,本文提出了一種更為科學(xué)合理的評價指標(biāo)——面積不重合度,即識別脫空與實際脫空不重合部分的面積在實際脫空的面積占比。從反演識別結(jié)果來看,識別脫空與實際脫空能較好地吻合,3種工況下水閘底板脫空反演結(jié)果與實際脫空的面積不重合度分別為14.79%,5.46%,7.36%,相比參考文獻(xiàn)[4]的反演結(jié)果有了較大的提升,驗證了本文方法的有效性。

    表6 不同工況下脫空參數(shù)反演結(jié)果Tab.6 Inversion results of void parameters under different conditions

    (a) 工況一

    表7 識別脫空面積與實際脫空面積的相對誤差Tab.7 Relative error between the identify void area and the actual void area

    3 結(jié)論

    (1) 基于顯著性分析,采用F檢驗,從7種不同形式的響應(yīng)面模型中優(yōu)選出形式最佳的響應(yīng)面模型;隨后,對形式最佳的響應(yīng)面模型的脫空參數(shù)及其組合進(jìn)行t檢驗以刪除贅余項,建立了能夠準(zhǔn)確高效地代替有限元模型的響應(yīng)面模型。

    (2) 提出了采用面積不重合度來衡量反演結(jié)果的準(zhǔn)確性,使得軟基水閘底板脫空識別的精度評價更加科學(xué)合理。

    (3) 將經(jīng)顯著性檢驗獲得的最佳的響應(yīng)面模型應(yīng)用于軟基水閘底板脫空反演,結(jié)果表明:識別脫空與實際脫空能夠較好地吻合,3種工況下的水閘底板脫空反演結(jié)果與實際脫空的面積不重合度分別為14.79%,5.46%,7.36%,提高了反演精度,驗證了本文方法的有效性。

    猜你喜歡
    水閘反演底板
    反演對稱變換在解決平面幾何問題中的應(yīng)用
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    2016年河南省已建成水閘數(shù)量
    全省已建成水閘數(shù)量
    軟土地基上超長底板節(jié)制閘沉降及應(yīng)力分析
    底板巷一巷兩用的回采工作面防火實踐
    底板隔水層破壞的力學(xué)模型及破壞判據(jù)
    河南省2014年已建成水閘數(shù)量
    疊前同步反演在港中油田的應(yīng)用
    成人二区视频| av国产精品久久久久影院| 日韩不卡一区二区三区视频在线| 一区二区三区精品91| 男人舔奶头视频| 午夜福利网站1000一区二区三区| 交换朋友夫妻互换小说| 三级国产精品欧美在线观看| 欧美精品人与动牲交sv欧美| 亚洲av欧美aⅴ国产| 国产乱人视频| 亚洲欧美一区二区三区黑人 | 国产亚洲一区二区精品| 晚上一个人看的免费电影| 国产av一区二区精品久久 | 天美传媒精品一区二区| 男人狂女人下面高潮的视频| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 免费久久久久久久精品成人欧美视频 | 亚洲欧美日韩卡通动漫| 成年免费大片在线观看| 色5月婷婷丁香| 亚洲精品亚洲一区二区| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| a级一级毛片免费在线观看| av在线蜜桃| 亚洲av.av天堂| 亚洲高清免费不卡视频| 纯流量卡能插随身wifi吗| 乱码一卡2卡4卡精品| 人妻一区二区av| 少妇人妻精品综合一区二区| 男女下面进入的视频免费午夜| 1000部很黄的大片| 国产高潮美女av| 久久精品人妻少妇| 五月伊人婷婷丁香| 少妇人妻久久综合中文| 精品一区二区免费观看| 一本色道久久久久久精品综合| 国产国拍精品亚洲av在线观看| 亚洲精品国产色婷婷电影| 日日摸夜夜添夜夜爱| 大片免费播放器 马上看| 亚洲aⅴ乱码一区二区在线播放| 18禁裸乳无遮挡免费网站照片| 国产人妻一区二区三区在| 涩涩av久久男人的天堂| 美女脱内裤让男人舔精品视频| 日本wwww免费看| 男男h啪啪无遮挡| 丝瓜视频免费看黄片| 黄色配什么色好看| 日韩大片免费观看网站| 国产人妻一区二区三区在| 日韩一区二区视频免费看| 寂寞人妻少妇视频99o| 各种免费的搞黄视频| 高清欧美精品videossex| 精品亚洲成国产av| 男的添女的下面高潮视频| 午夜激情久久久久久久| 黄色视频在线播放观看不卡| 欧美老熟妇乱子伦牲交| 少妇丰满av| 久久人人爽av亚洲精品天堂 | 91狼人影院| kizo精华| 国产亚洲最大av| 久久精品久久精品一区二区三区| 精品久久久精品久久久| 国产精品嫩草影院av在线观看| 极品教师在线视频| a级一级毛片免费在线观看| 欧美精品一区二区大全| 黑人高潮一二区| 丰满少妇做爰视频| 女的被弄到高潮叫床怎么办| 久久久久久久久大av| 尾随美女入室| 亚洲国产高清在线一区二区三| 国产精品偷伦视频观看了| 亚洲精品日韩av片在线观看| 成人国产av品久久久| 国产高潮美女av| 少妇人妻久久综合中文| 欧美成人一区二区免费高清观看| 纯流量卡能插随身wifi吗| 免费看光身美女| 天天躁夜夜躁狠狠久久av| 草草在线视频免费看| 国产精品女同一区二区软件| 亚洲aⅴ乱码一区二区在线播放| 六月丁香七月| av不卡在线播放| 最近的中文字幕免费完整| 韩国高清视频一区二区三区| 亚洲图色成人| 80岁老熟妇乱子伦牲交| 亚洲,欧美,日韩| 国产在线免费精品| 亚洲电影在线观看av| 免费人妻精品一区二区三区视频| 最新中文字幕久久久久| 五月玫瑰六月丁香| 深爱激情五月婷婷| 下体分泌物呈黄色| 久久精品国产亚洲av天美| 99热6这里只有精品| 大码成人一级视频| 777米奇影视久久| 久久99热这里只频精品6学生| 亚洲欧美清纯卡通| 国产精品久久久久久精品电影小说 | 国产91av在线免费观看| av专区在线播放| 免费看光身美女| 大话2 男鬼变身卡| 亚洲精品日韩av片在线观看| 一级a做视频免费观看| 久久国产亚洲av麻豆专区| 亚洲国产精品成人久久小说| 香蕉精品网在线| 久久久久久久亚洲中文字幕| 麻豆成人av视频| 你懂的网址亚洲精品在线观看| 亚洲熟女精品中文字幕| 日本vs欧美在线观看视频 | 一级毛片我不卡| 亚洲精品一区蜜桃| 男女边摸边吃奶| 视频中文字幕在线观看| tube8黄色片| 国产成人精品久久久久久| 一本—道久久a久久精品蜜桃钙片| 精品久久国产蜜桃| 欧美日韩一区二区视频在线观看视频在线| 欧美亚洲 丝袜 人妻 在线| 美女国产视频在线观看| 国产高清国产精品国产三级 | 少妇 在线观看| 直男gayav资源| 97超碰精品成人国产| 日本-黄色视频高清免费观看| 久久99精品国语久久久| 少妇人妻久久综合中文| 日本爱情动作片www.在线观看| 亚洲精品日本国产第一区| 老师上课跳d突然被开到最大视频| 国产精品国产av在线观看| 国产熟女欧美一区二区| 日韩伦理黄色片| 成人18禁高潮啪啪吃奶动态图 | 日韩欧美 国产精品| 国产精品av视频在线免费观看| 中文字幕精品免费在线观看视频 | 91久久精品国产一区二区成人| 国产高清三级在线| 精品人妻视频免费看| 欧美xxxx黑人xx丫x性爽| 男女无遮挡免费网站观看| 日日啪夜夜撸| 2021少妇久久久久久久久久久| 久久久欧美国产精品| 久久精品国产自在天天线| 久久国内精品自在自线图片| 国产成人一区二区在线| 国产男人的电影天堂91| 99热这里只有是精品在线观看| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 免费在线观看成人毛片| 成年美女黄网站色视频大全免费 | 一区二区av电影网| 久久人人爽人人片av| 九草在线视频观看| 黑人猛操日本美女一级片| 久热久热在线精品观看| 寂寞人妻少妇视频99o| 偷拍熟女少妇极品色| 久久精品国产a三级三级三级| 两个人的视频大全免费| av国产精品久久久久影院| 亚洲国产毛片av蜜桃av| 亚洲欧美精品专区久久| 国产精品不卡视频一区二区| 这个男人来自地球电影免费观看 | 免费久久久久久久精品成人欧美视频 | 成人综合一区亚洲| 免费播放大片免费观看视频在线观看| 天天躁日日操中文字幕| 国产伦理片在线播放av一区| 国产成人91sexporn| 亚洲国产色片| 看十八女毛片水多多多| 激情 狠狠 欧美| 中文字幕免费在线视频6| 亚洲欧美一区二区三区黑人 | 久久久久精品性色| 亚洲av二区三区四区| 国产男女内射视频| 亚洲精品乱码久久久v下载方式| 高清午夜精品一区二区三区| 啦啦啦中文免费视频观看日本| 亚洲av成人精品一区久久| 午夜福利在线在线| 老司机影院成人| 亚洲性久久影院| 亚洲国产精品专区欧美| 久久久久性生活片| 人妻夜夜爽99麻豆av| 综合色丁香网| 国产又色又爽无遮挡免| 国产亚洲91精品色在线| 亚洲国产精品国产精品| 国产 一区精品| 大香蕉97超碰在线| 亚洲经典国产精华液单| 色视频在线一区二区三区| 亚洲第一av免费看| 亚洲怡红院男人天堂| 午夜激情福利司机影院| 身体一侧抽搐| 欧美日韩国产mv在线观看视频 | 日日摸夜夜添夜夜爱| 国产色婷婷99| 亚洲欧美成人综合另类久久久| 国产片特级美女逼逼视频| 国产在线免费精品| 国产综合精华液| 一级爰片在线观看| 日韩制服骚丝袜av| 久热久热在线精品观看| av在线观看视频网站免费| 一级毛片我不卡| 国产毛片在线视频| 日韩三级伦理在线观看| av黄色大香蕉| 亚洲av成人精品一二三区| 欧美xxxx性猛交bbbb| 亚洲成色77777| 搡老乐熟女国产| 中文字幕人妻熟人妻熟丝袜美| 汤姆久久久久久久影院中文字幕| 丝袜脚勾引网站| 成人免费观看视频高清| 观看美女的网站| 联通29元200g的流量卡| av免费在线看不卡| 性色av一级| 欧美区成人在线视频| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 爱豆传媒免费全集在线观看| 又爽又黄a免费视频| 国产极品天堂在线| 久久久久精品性色| 国产免费福利视频在线观看| 成人无遮挡网站| 亚洲性久久影院| 国产精品免费大片| 成人美女网站在线观看视频| 最近最新中文字幕免费大全7| 成人二区视频| 亚洲精品视频女| 国产在视频线精品| 亚洲av二区三区四区| 欧美一区二区亚洲| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 久久人人爽人人片av| 国产免费福利视频在线观看| 亚洲欧美清纯卡通| 九九在线视频观看精品| av卡一久久| 尾随美女入室| 舔av片在线| 男女无遮挡免费网站观看| 少妇人妻久久综合中文| 国产深夜福利视频在线观看| 五月伊人婷婷丁香| 热re99久久精品国产66热6| 亚洲av欧美aⅴ国产| 国精品久久久久久国模美| 91精品一卡2卡3卡4卡| 日韩中文字幕视频在线看片 | 亚洲国产av新网站| 国产亚洲一区二区精品| 亚洲真实伦在线观看| 国产人妻一区二区三区在| 国产有黄有色有爽视频| 91久久精品国产一区二区成人| 国产伦理片在线播放av一区| 一级a做视频免费观看| av在线老鸭窝| 日日撸夜夜添| 国产乱人偷精品视频| 日本色播在线视频| 少妇熟女欧美另类| 男女啪啪激烈高潮av片| 精品人妻熟女av久视频| 国产精品久久久久久久电影| 18禁裸乳无遮挡免费网站照片| 亚洲中文av在线| 青春草亚洲视频在线观看| 秋霞伦理黄片| 日韩中文字幕视频在线看片 | 成人特级av手机在线观看| 日本欧美视频一区| 亚洲自偷自拍三级| 少妇精品久久久久久久| 亚洲一级一片aⅴ在线观看| 亚洲av不卡在线观看| 亚洲精品,欧美精品| 简卡轻食公司| 国产成人精品久久久久久| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看 | 亚洲av二区三区四区| av国产精品久久久久影院| 亚洲无线观看免费| 久久影院123| 久热这里只有精品99| 国产男人的电影天堂91| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 精品视频人人做人人爽| 一边亲一边摸免费视频| 成人一区二区视频在线观看| 国产精品免费大片| a级一级毛片免费在线观看| 久久人人爽人人爽人人片va| 久久精品熟女亚洲av麻豆精品| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 少妇人妻 视频| 午夜精品国产一区二区电影| 日韩不卡一区二区三区视频在线| av线在线观看网站| 日韩精品有码人妻一区| 内地一区二区视频在线| 亚洲国产精品专区欧美| 极品少妇高潮喷水抽搐| 亚洲av.av天堂| 亚洲国产精品国产精品| 中国三级夫妇交换| 又黄又爽又刺激的免费视频.| 国产精品一区二区在线观看99| 成人漫画全彩无遮挡| 少妇的逼水好多| 2018国产大陆天天弄谢| 国产欧美日韩精品一区二区| 亚洲欧洲日产国产| 另类亚洲欧美激情| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看| 欧美性感艳星| 亚洲国产日韩一区二区| 久久6这里有精品| 一个人看视频在线观看www免费| 国产精品爽爽va在线观看网站| 亚洲人成网站在线观看播放| 欧美国产精品一级二级三级 | 国产永久视频网站| 少妇熟女欧美另类| 91aial.com中文字幕在线观看| 国产精品国产三级国产专区5o| 婷婷色麻豆天堂久久| 最近的中文字幕免费完整| 日韩av在线免费看完整版不卡| 久久久久国产精品人妻一区二区| 国产爱豆传媒在线观看| 最近2019中文字幕mv第一页| 亚洲精品,欧美精品| 久久精品国产亚洲网站| 亚洲天堂av无毛| 如何舔出高潮| 国内揄拍国产精品人妻在线| 99热这里只有精品一区| 日本vs欧美在线观看视频 | 人体艺术视频欧美日本| 热99国产精品久久久久久7| av.在线天堂| 在现免费观看毛片| 交换朋友夫妻互换小说| 一级黄片播放器| 日韩av不卡免费在线播放| 狂野欧美白嫩少妇大欣赏| 国产乱人视频| 亚洲成人一二三区av| 黄片无遮挡物在线观看| 如何舔出高潮| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| 亚洲国产欧美在线一区| av专区在线播放| 一级毛片电影观看| 亚洲久久久国产精品| 天美传媒精品一区二区| 少妇高潮的动态图| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 91精品国产国语对白视频| 亚洲精华国产精华液的使用体验| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有精品一区| 熟女人妻精品中文字幕| 国产欧美亚洲国产| 亚洲av成人精品一二三区| 亚洲国产色片| 成人高潮视频无遮挡免费网站| 久久人人爽人人爽人人片va| 香蕉精品网在线| 久久精品久久久久久久性| 少妇 在线观看| 亚洲国产精品999| 少妇人妻 视频| 国产黄色视频一区二区在线观看| 亚洲最大成人中文| 久久久国产一区二区| 国产爽快片一区二区三区| 自拍偷自拍亚洲精品老妇| 亚洲国产色片| 国产熟女欧美一区二区| 国产成人免费无遮挡视频| 色5月婷婷丁香| 免费av中文字幕在线| 一级爰片在线观看| 亚洲综合色惰| 亚洲精品国产av蜜桃| 久久99蜜桃精品久久| 国产人妻一区二区三区在| 亚洲欧美一区二区三区黑人 | 午夜福利在线在线| 国产 精品1| 亚洲av日韩在线播放| 97热精品久久久久久| 久久精品国产亚洲av涩爱| 在线观看国产h片| 秋霞伦理黄片| 中文精品一卡2卡3卡4更新| tube8黄色片| 国产精品一及| 直男gayav资源| 国内揄拍国产精品人妻在线| 一级av片app| 蜜臀久久99精品久久宅男| 成年人午夜在线观看视频| 欧美日韩综合久久久久久| 少妇人妻久久综合中文| 插阴视频在线观看视频| 少妇高潮的动态图| av.在线天堂| 男女免费视频国产| 国产成人免费无遮挡视频| 高清在线视频一区二区三区| 国产精品无大码| 在线免费十八禁| h日本视频在线播放| 欧美性感艳星| 午夜免费男女啪啪视频观看| 欧美另类一区| 国产淫语在线视频| 在线观看三级黄色| 国产 一区 欧美 日韩| 精品久久久久久电影网| 成人综合一区亚洲| 久久人人爽av亚洲精品天堂 | 丰满迷人的少妇在线观看| 成人高潮视频无遮挡免费网站| 精品久久久精品久久久| 亚洲精品亚洲一区二区| 国产爽快片一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 在线播放无遮挡| 美女福利国产在线 | 在线观看免费高清a一片| 全区人妻精品视频| 最近的中文字幕免费完整| 免费观看a级毛片全部| 尤物成人国产欧美一区二区三区| 亚洲欧洲日产国产| 国产成人精品福利久久| 成年人午夜在线观看视频| 中文欧美无线码| 精品人妻一区二区三区麻豆| a级毛片免费高清观看在线播放| 久久精品国产亚洲av涩爱| 国语对白做爰xxxⅹ性视频网站| 精品人妻偷拍中文字幕| 狠狠精品人妻久久久久久综合| 国产一区亚洲一区在线观看| 麻豆成人av视频| 最后的刺客免费高清国语| 久久久亚洲精品成人影院| 极品教师在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美日韩视频精品一区| 亚洲国产毛片av蜜桃av| 久久综合国产亚洲精品| 国产精品久久久久久久电影| 建设人人有责人人尽责人人享有的 | 国产av精品麻豆| 下体分泌物呈黄色| 精品国产乱码久久久久久小说| 国产色爽女视频免费观看| 国产精品一区二区在线不卡| 成人亚洲精品一区在线观看 | 国产精品久久久久久久电影| 狂野欧美激情性bbbbbb| 26uuu在线亚洲综合色| 大又大粗又爽又黄少妇毛片口| 成人漫画全彩无遮挡| 99精国产麻豆久久婷婷| 在线免费观看不下载黄p国产| 视频区图区小说| 青春草视频在线免费观看| 最近中文字幕高清免费大全6| 欧美xxⅹ黑人| 男人狂女人下面高潮的视频| 亚洲精品色激情综合| 高清视频免费观看一区二区| 99久久精品国产国产毛片| 国产精品成人在线| 少妇精品久久久久久久| 天美传媒精品一区二区| 三级国产精品欧美在线观看| 极品教师在线视频| 2018国产大陆天天弄谢| 女人久久www免费人成看片| 国产淫片久久久久久久久| 亚洲成人中文字幕在线播放| 成年av动漫网址| 一本色道久久久久久精品综合| 亚洲三级黄色毛片| 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 国产大屁股一区二区在线视频| 国产亚洲91精品色在线| 久久久欧美国产精品| 日韩成人av中文字幕在线观看| 亚洲最大成人中文| videossex国产| 我的老师免费观看完整版| 夜夜骑夜夜射夜夜干| 纵有疾风起免费观看全集完整版| 在线免费观看不下载黄p国产| 街头女战士在线观看网站| 晚上一个人看的免费电影| 3wmmmm亚洲av在线观看| 国产 一区 欧美 日韩| 一级a做视频免费观看| 91狼人影院| 午夜激情久久久久久久| 小蜜桃在线观看免费完整版高清| 久久久久久久精品精品| 老司机影院毛片| 国产精品久久久久成人av| 国产黄色免费在线视频| 热re99久久精品国产66热6| 菩萨蛮人人尽说江南好唐韦庄| 国产淫语在线视频| 视频中文字幕在线观看| 欧美人与善性xxx| 亚洲国产最新在线播放| 在线看a的网站| 在线免费十八禁| 亚洲精品,欧美精品| 欧美成人精品欧美一级黄| 各种免费的搞黄视频| 午夜视频国产福利| 欧美性感艳星| 少妇高潮的动态图| 一级毛片我不卡| 亚洲成人av在线免费| 久久久久网色| 欧美精品一区二区免费开放| 中文字幕制服av| 少妇丰满av| 新久久久久国产一级毛片| 亚州av有码| 18+在线观看网站| 最近中文字幕高清免费大全6| 精品国产一区二区三区久久久樱花 | 18禁在线播放成人免费| 国产精品.久久久| 五月天丁香电影| av在线观看视频网站免费| 久久精品国产亚洲av天美| 各种免费的搞黄视频| 成人午夜精彩视频在线观看| av在线播放精品| 国产成人aa在线观看| 国产黄色免费在线视频| 国产v大片淫在线免费观看| 亚洲美女视频黄频| 国产淫语在线视频| 免费大片黄手机在线观看| 观看美女的网站| 永久网站在线| 国产精品久久久久久精品电影小说 | 纵有疾风起免费观看全集完整版| 美女高潮的动态| 亚洲国产精品专区欧美| 欧美变态另类bdsm刘玥| 亚洲,一卡二卡三卡| 亚洲av日韩在线播放| 蜜桃亚洲精品一区二区三区| 久久久久视频综合| 尾随美女入室| 18禁裸乳无遮挡免费网站照片|