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

    基于信息熵的動力模型縮聚中主自由度優(yōu)化選擇研究①

    2018-08-03 05:21:10,
    關(guān)鍵詞:振型桁架模態(tài)

    , ,

    (武漢大學(xué)土木建筑工程學(xué)院,湖北 武漢 430072)

    0 引 言

    許多的橋梁工程以及高層建筑等結(jié)構(gòu)在投入使用的過程中都會出現(xiàn)一些局部的破損,隨著時(shí)間的推移,各種因素的影響會使局部破損不斷擴(kuò)大,造成整個(gè)結(jié)構(gòu)的破壞[1]。在工程實(shí)踐中,由于條件限制,為了能更好的對工程結(jié)構(gòu)進(jìn)行檢測和維護(hù),要將有限的傳感器布設(shè)在結(jié)構(gòu)的主要部位,獲取更為全面的信息。

    國內(nèi)外學(xué)者對于傳感器的優(yōu)化布置問題提出了多種優(yōu)化算法,Guyan[2]模型縮聚法是一種常用的傳感器布置方法,Kammer[3]提出了傳感器優(yōu)化布置的有效獨(dú)立法。在國內(nèi),崔飛[4]提出了模態(tài)置信度法,根據(jù)經(jīng)驗(yàn)和初算來擬定一組測點(diǎn),然后從剩余的自由度中每次增加一個(gè)測點(diǎn),以此減小MAC矩陣中非對角元的最大值,來達(dá)到測點(diǎn)優(yōu)化的目的。

    針對桁架橋梁結(jié)構(gòu)模型,運(yùn)用動態(tài)模型縮聚法進(jìn)行主自由度的優(yōu)化選取,以橋梁建模參數(shù)識別結(jié)果的不確定性最小為準(zhǔn)則來優(yōu)化選取主自由度,對該橋梁結(jié)構(gòu)模型進(jìn)行數(shù)值仿真研究以及實(shí)驗(yàn)驗(yàn)證。

    1 動力模型縮聚

    有一類別的模型,假設(shè)有N個(gè)自由度,假如質(zhì)量矩陣由于結(jié)構(gòu)損傷的改變可以忽略,可以得到系統(tǒng)方程如下[5]:

    K(θ)=λjMφj

    (1)

    其中,λj,φj∈RN,j=1,2,…Nt是方程的特征值和特征向量。Nt是測得的模態(tài)階數(shù)。K(θ)是參數(shù)化之后的整體剛度矩陣,其中θ=[θ1,θ2,…,θNθ]T∈RNθ:

    (2)

    當(dāng)θ=θ0時(shí),代表結(jié)構(gòu)處于健康狀況下;K(i),i=1,2,…Nθ是組裝成整體剛度矩陣的單元?jiǎng)偠染仃嚨臄U(kuò)階矩陣;Nθ是待識別的結(jié)構(gòu)參數(shù)的個(gè)數(shù)。

    將上述公式聯(lián)合,并區(qū)分主副自由度,可以得到以下振動方程:

    (3)

    其中{φm}j和{φs}j分別是第j階振型測得的和未測得的部分,維度分別是Nm和Ns,且有Nm+NS=N。

    根據(jù)公式(3)可得:

    {φs}j=Dj{φm}j

    (4)

    (5)

    其中Dj是第j階模態(tài)的動態(tài)縮減矩陣,是關(guān)于λj和θ的函數(shù),如下:

    (6)

    (7)

    因此,對于第j階振型,可以只用主振型表示如下:

    (8)

    其中Tj∈RNNm是第j階模態(tài)的轉(zhuǎn)換矩陣,INm∈RNmNm是單位矩陣。

    將公式(7)帶入公式(3),并且方程兩邊同時(shí)乘以T矩陣的轉(zhuǎn)置矩陣得到:

    (9)

    其中:

    (10)

    (11)

    (12)

    (13)

    (14)

    根據(jù)公式(11)可得[6]:

    A(θ)θ=b(θ)

    (15)

    其中:

    (16)

    (17)

    且有:

    (18)

    zj=zj(λj,{φm}j,θ)=

    (19)

    根據(jù)公式(15)可以得到:

    θ=A(θ)+b(θ)

    (20)

    公式中的A(θ)+是A(θ)矩陣的偽逆矩陣,給定θ初始值θ0,k=1,每步迭代結(jié)果為θk=A(θk-1)+b(θk-1),當(dāng)兩代θ之間的改變值小于一個(gè)指定的極小值,則迭代結(jié)束,θk為識別結(jié)果,否則迭代繼續(xù)。

    2 主自由度優(yōu)化選取理論

    2.1 貝葉斯理論

    假設(shè)選取特定的某一類模型,模型中存在著結(jié)構(gòu)自由參數(shù),假設(shè)模型參數(shù)向量為θ,利用貝葉斯統(tǒng)計(jì)方法根據(jù)時(shí)域內(nèi)的動力實(shí)驗(yàn)數(shù)據(jù),推得結(jié)構(gòu)的合理參數(shù)和不確定性。

    將測得的實(shí)驗(yàn)數(shù)據(jù)記為D={y(m),z(m),m=1,…,N},y(m)∈RN0,z(m)∈RNT分別是輸入和輸出數(shù)據(jù)。N0表示結(jié)構(gòu)觀測自由度的個(gè)數(shù),NT表示輸入自由度的個(gè)數(shù),m表示時(shí)間指數(shù)mΔt,mΔt表示采樣時(shí)間間隔,N表示采樣數(shù)據(jù)的個(gè)數(shù)。當(dāng)選定一類模型后,就可以獲得模型響應(yīng)輸入與輸出之間的關(guān)系[7]:

    y(n)=L0q(m;θ)+n(m;θ)

    (21)

    以下假設(shè)預(yù)估誤差對于不同的測量數(shù)據(jù)滿足∑=σ2I,σ是預(yù)估誤差參數(shù),I∈RN0N0是單位矩陣。根據(jù)貝葉斯系統(tǒng)識別理論,隨機(jī)的參數(shù)向量φ和測量數(shù)據(jù)D,得到模型參數(shù)θ的先驗(yàn)概率密度如下:

    P(ψ,θ|M)=P(ψ|θ,M)P(θ|M)

    (22)

    其中:

    P(θ|M)=N(θ0;Eθ)

    (23)

    先驗(yàn)概率密度函數(shù)如下[8〗:

    (24)

    其中c是常數(shù),使得證據(jù)因子函數(shù)積分為1,而且有:

    (25)

    Σ=[(Σ1)T,(Σ2)T,…,(ΣNt)T]T

    (26)

    公式(26)中的Σj表示振動方程中第j模態(tài)向量的誤差:

    (27)

    根據(jù)上文可得:

    (28)

    同理,可以得到模型修正后的后驗(yàn)概率密度函數(shù):

    P(ψ,θ|D,M)=c2P(D|ψ,θ,M)P(ψ,θ|M)

    (29)

    由此可得:

    (30)

    (31)

    2.2 信息熵理論

    修正的概率密度函數(shù)表示模態(tài)參數(shù)每個(gè)可能值的概率,根據(jù)測得的動力數(shù)據(jù)可以得到參數(shù)值的不確定性。信息熵作為唯一量測結(jié)構(gòu)模型參數(shù)θ的不確定性的工具,表達(dá)方式如下[9〗:

    H(δ,D)=Eθ[-lnp(θ|D)]=

    (32)

    其中Eθ表示關(guān)于θ的數(shù)學(xué)期望。信息熵理論取決于測得的動力數(shù)據(jù)D和選取的主自由度δ。

    根據(jù)C. Papadimitriou提出的理論可以得到信息熵理論[10]:

    -Eθ[lnπθ(θ)]

    (33)

    (34)

    (35)

    (36)

    (37)

    (38)

    s/s0的值越小表明選取的主自由度組合越優(yōu),對于同一結(jié)構(gòu)不同數(shù)量的主自由度組合s/s0也可以進(jìn)行比較。

    3 主自由度優(yōu)化選取

    3.1 MATLAB有限元建模及其驗(yàn)證

    采用自編MATLAB程序建立簡單的桁架橋模型。其中縱梁采用T型截面,橫梁采用工字型截面,彈性模量為2.06×1011kN/m2,剪切模量為7.923×1010MPa,截面面積為0.904×10-3m2,y軸慣性矩為0.145×10-5m4,z軸慣性矩為0.21×10-6m4,泊松比為0.3,材料密度為7850kg/m3。如下圖1所示:

    圖1中所示為桁架橋ANSYS模型圖,4個(gè)端節(jié)點(diǎn)固支。采用ANSYS與MATLAB軟件對該桁架橋建模并進(jìn)行模態(tài)分析,求得前五階模態(tài)參數(shù)。

    表1 MATLAB與ANSYS計(jì)算頻率(Hz)對比

    圖1 桁架橋ANSYS模型圖

    其中表1為頻率結(jié)果,圖2為振型結(jié)果,藍(lán)色線代表的是桁架橋原位曲線,紅色線和黑色線分別為MATLAB和ANSYS提取的振型結(jié)果。從表1以及圖2的模態(tài)結(jié)果對比,可以看出兩種軟件結(jié)果吻合較好,充分說明所建模型的準(zhǔn)確性。

    3.2 主自由度優(yōu)化選取仿真研究

    上圖3為桁架橋模型下部的示意圖,根據(jù)結(jié)構(gòu)對稱性以及振型特點(diǎn),將節(jié)點(diǎn)的豎向自由度作為主要自由度,故可選取的主自由度為1-18號節(jié)點(diǎn)的豎向自由度。待識別建模參數(shù)取為桁架橋下方前八個(gè)縱梁單元的剛度矩陣,共8個(gè)。采用遺傳算法分別獲得6和8個(gè)主自由度組合時(shí)的優(yōu)化結(jié)果,并采用窮舉法驗(yàn)證,因組合數(shù)較多,不列表顯示。

    圖2 MATLAB與ANSYS前四階振型結(jié)果

    圖3 桁架橋主自由度示意圖

    為進(jìn)一步驗(yàn)證,主自由度優(yōu)化結(jié)果的正確性,選取部分不同的主自由度組合與優(yōu)化得到的主自由度組合,進(jìn)行損傷識別結(jié)果對比;一共各選取4種非最優(yōu)組合,以下是6個(gè)主自由度組合識別結(jié)果對比:

    表2 6個(gè)主自由度的各參數(shù)識別結(jié)果的標(biāo)準(zhǔn)差

    將表2中的結(jié)果繪制成圖,如下圖:

    圖4 6個(gè)主自由度的參數(shù)識別結(jié)果標(biāo)準(zhǔn)差

    上圖4中紅色線為優(yōu)化得到的主自由度組合,其余為非最優(yōu)組合。從圖4中可以看出,除了極少數(shù)的參數(shù)標(biāo)準(zhǔn)差偏大之外,優(yōu)化得到的主自由度組合識別結(jié)果標(biāo)準(zhǔn)差絕大部分都是最小的,也就說明,這個(gè)主自由度組合更優(yōu),參數(shù)識別結(jié)果的不確定性最小。其中主自由度的編號為2、22、32、42、47、52〗,也即是下圖中紅色箭頭所標(biāo)注節(jié)點(diǎn)的豎向自由度組合,如圖5所示。

    圖5 6個(gè)優(yōu)化的主自由度組合示意圖

    圖6 8個(gè)主自由度組合的參數(shù)識別結(jié)果標(biāo)準(zhǔn)差

    以上給出主自由度個(gè)數(shù)為6時(shí)的最優(yōu)主自由度組合。為進(jìn)一步說明結(jié)果的正確性,將主自由度的個(gè)數(shù)定為8個(gè),也對應(yīng)選取4種不同的主自由度組合進(jìn)行對比,識別結(jié)果如下:

    表3 8個(gè)主自由度的各參數(shù)識別結(jié)果的標(biāo)準(zhǔn)差

    將表3中的結(jié)果繪制成圖,如下圖:

    上圖6中每條曲線的含義同圖4相同,紅色線為遺傳算法優(yōu)化得到的主自由度組合,其余為任意選取的與之對比的主自由度組合。每個(gè)參數(shù)的標(biāo)準(zhǔn)差越小,表示識別結(jié)果更穩(wěn)定,不確定性更小。根據(jù)上文給出的結(jié)論,可以得到,優(yōu)化得到的主自由度組合識別結(jié)果標(biāo)準(zhǔn)差絕大部分都是最小的,也就說明,這個(gè)主自由度組合更優(yōu);此種主自由度組合下各個(gè)自由度編號為2、12、22、27、47、52、82、87,如圖7所示,紅色箭頭所標(biāo)注節(jié)點(diǎn)的豎向自由度組合。

    圖7 8個(gè)優(yōu)化的主自由度組合示意圖

    圖8 (a)桁架橋模型 (b)實(shí)驗(yàn)裝置系統(tǒng)

    為進(jìn)一步分析主自由度對比選取的結(jié)果,將上文中設(shè)置的兩種工況,分別選定6個(gè)和8個(gè)主自由度的優(yōu)化組合結(jié)果進(jìn)行對比如下表4所示,從表4中可以看出主自由度的選取個(gè)數(shù)增加會使s/s0的值變小,說明隨著主自由度個(gè)數(shù)的增多,計(jì)算中結(jié)構(gòu)的所包含的信息量逐漸增大,這也說明結(jié)構(gòu)中待識別建模參數(shù)識別結(jié)果的不確定性逐漸減小。

    表4 桁架橋兩種工況下計(jì)算結(jié)果表

    3.3 實(shí)驗(yàn)驗(yàn)證

    該桁架橋模型的幾何參數(shù)與第一節(jié)給出的數(shù)值仿真參數(shù)一致,符合工程設(shè)計(jì),以下為實(shí)驗(yàn)室桁架橋模型以及實(shí)驗(yàn)裝置系統(tǒng),如圖8(a)和(b)所示:

    實(shí)驗(yàn)室桁架橋模態(tài)提取試驗(yàn)采用的實(shí)驗(yàn)裝置系統(tǒng)主要有桁架橋模型、數(shù)據(jù)采集儀(KT6016),N-Modal模態(tài)分析軟件、IEPE型加速度傳感器(KT11000L)、沖擊錘、筆記本電腦以及若干導(dǎo)線等。

    圖9 桁架橋?qū)嶒?yàn)流程圖

    圖10 (a) 5號傳感器與12傳感器的加速度時(shí)程響應(yīng)比較

    圖11 (a) 5號節(jié)點(diǎn)與12號節(jié)點(diǎn)的傳感器的FFT圖

    此桁架橋模態(tài)提取實(shí)驗(yàn)首先利用沖擊錘對桁架橋模型施加一個(gè)垂直短時(shí)的沖擊荷載,然后利用吸附在橫梁上的加速度傳感器獲得桁架橋在豎直方向的自由響應(yīng)信號,此時(shí)利用數(shù)據(jù)采集儀采集加速度響應(yīng)歷程數(shù)據(jù)并保存至筆記本電腦中,最后運(yùn)用N-Modal模態(tài)分析軟件提取桁架橋的頻率與振型的參數(shù)信息,具體操作流程如下圖9所示。

    采樣頻率為100Hz,采樣時(shí)間為10s,選取其中的5號節(jié)點(diǎn)、7號節(jié)點(diǎn)以及12號節(jié)點(diǎn)進(jìn)行加速度響應(yīng)分析,各節(jié)點(diǎn)的加速度響應(yīng)曲線如下圖10所示,紅色線是5號節(jié)點(diǎn)上的傳感器采集到的加速度響應(yīng)信號曲線,藍(lán)色線為7號與12號節(jié)點(diǎn)的傳感器采集到的加速度響應(yīng)信號曲線,從圖中可以看出5號節(jié)點(diǎn)的傳感器的響應(yīng)信號與7、12節(jié)點(diǎn)的加速度信號存在一定區(qū)別,相對而言5號傳感器位于振型的中心,故而采集得到的信息量更多。

    為了驗(yàn)證仿真結(jié)果,當(dāng)選取的主自由度個(gè)數(shù)為1個(gè)時(shí),將最優(yōu)主自由度以及其中部分非最優(yōu)主自由度對應(yīng)的目標(biāo)函數(shù)值列于下表5所示。

    表5 不同主自由度對應(yīng)的目標(biāo)函數(shù)值

    從上表5中可知,最優(yōu)主自由度為5號節(jié)點(diǎn)對應(yīng)的豎向自由度,將傳感器分別布置在最優(yōu)的主自由度位置即5號節(jié)點(diǎn)的豎向自由度,與另外非最優(yōu)位置的7號、12號節(jié)點(diǎn)的豎向自由度,得到不同的響應(yīng)FFT圖,如圖11所示。

    圖11中,紅色曲線為5號節(jié)點(diǎn)處的傳感器的時(shí)程響應(yīng)曲線,藍(lán)色曲線分別為7、12號節(jié)點(diǎn)處傳感器對應(yīng)的響應(yīng)信息。從每個(gè)圖中都可以很清晰的看出,所用方法選擇出的最優(yōu)主自由度即5號節(jié)點(diǎn)的豎向自由度得到的頻率成分更加豐富,說明了所采用方法的正確性。

    5 結(jié) 論

    采用文中的方法可得到優(yōu)化的主自由度組合,表明貝葉斯統(tǒng)計(jì)系統(tǒng)識別方法能夠很方便地識別出待識別結(jié)構(gòu)建模參數(shù)的最優(yōu)值以及相應(yīng)的不確定性;信息熵指標(biāo)能夠定量的評估待識別參數(shù)的不確定性程度。增加主自由度的數(shù)量可減小結(jié)構(gòu)待識別參數(shù)的不確定性。以桁架橋?yàn)槔砻髂軌蜻m合于橋梁結(jié)構(gòu)的主自由度優(yōu)化選取問題求解。

    猜你喜歡
    振型桁架模態(tài)
    桁架式吸泥機(jī)改造
    關(guān)于模態(tài)綜合法的注記
    縱向激勵(lì)下大跨鋼桁拱橋高階振型效應(yīng)分析
    擺臂式復(fù)合桁架機(jī)器人的開發(fā)
    塔腿加過渡段輸電塔動力特性分析
    Loader軸在雙機(jī)桁架機(jī)械手上的應(yīng)用
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    矮寨特大懸索橋鋼桁架安裝
    結(jié)構(gòu)振型幾何辨識及應(yīng)用研究
    山西建筑(2015年14期)2015-06-05 09:37:07
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    中文字幕亚洲精品专区| 天天一区二区日本电影三级| 国产精品女同一区二区软件| 亚洲欧美日韩无卡精品| 国产麻豆成人av免费视频| 日本黄大片高清| 麻豆国产97在线/欧美| 三级男女做爰猛烈吃奶摸视频| 又爽又黄a免费视频| 日本爱情动作片www.在线观看| 伦精品一区二区三区| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 91久久精品国产一区二区三区| 国产免费福利视频在线观看| 男人狂女人下面高潮的视频| 亚洲精品日韩av片在线观看| 美女被艹到高潮喷水动态| 欧美激情在线99| 欧美一区二区国产精品久久精品| av免费观看日本| 色综合色国产| 中文乱码字字幕精品一区二区三区 | 可以在线观看毛片的网站| 亚洲欧美精品综合久久99| 99在线人妻在线中文字幕| 又粗又硬又长又爽又黄的视频| 国产亚洲精品av在线| 欧美一区二区国产精品久久精品| 色噜噜av男人的天堂激情| 中文字幕久久专区| 少妇被粗大猛烈的视频| 亚洲国产日韩欧美精品在线观看| eeuss影院久久| 看黄色毛片网站| 国产亚洲一区二区精品| 日本色播在线视频| 你懂的网址亚洲精品在线观看 | 亚洲av成人av| av在线蜜桃| 国产精品99久久久久久久久| 夜夜爽夜夜爽视频| a级毛色黄片| 国产人妻一区二区三区在| 国产黄片视频在线免费观看| 亚洲精品自拍成人| 欧美变态另类bdsm刘玥| 日韩视频在线欧美| 一本一本综合久久| 国产淫语在线视频| 久久人人爽人人爽人人片va| 一区二区三区免费毛片| 少妇猛男粗大的猛烈进出视频 | 欧美三级亚洲精品| 日本五十路高清| 国内揄拍国产精品人妻在线| 精品人妻偷拍中文字幕| 国产精品一区二区性色av| 精品国产三级普通话版| 欧美激情国产日韩精品一区| 国产综合懂色| 一区二区三区高清视频在线| av视频在线观看入口| 91在线精品国自产拍蜜月| 国产午夜福利久久久久久| 人妻系列 视频| kizo精华| 亚洲熟妇中文字幕五十中出| 国产乱人视频| 成年av动漫网址| 天堂影院成人在线观看| 国产 一区精品| 午夜福利视频1000在线观看| 五月玫瑰六月丁香| 十八禁国产超污无遮挡网站| 看免费成人av毛片| 日韩av在线免费看完整版不卡| 日韩欧美三级三区| 精品欧美国产一区二区三| 国产精品久久久久久久久免| 欧美高清成人免费视频www| 国产免费男女视频| 日本五十路高清| 亚洲18禁久久av| 亚洲欧美日韩卡通动漫| 超碰97精品在线观看| 久久精品熟女亚洲av麻豆精品 | h日本视频在线播放| 精品国产一区二区三区久久久樱花 | a级毛色黄片| 亚洲在线观看片| 国产毛片a区久久久久| 亚洲精品自拍成人| 又粗又爽又猛毛片免费看| 97超视频在线观看视频| 桃色一区二区三区在线观看| 2022亚洲国产成人精品| 亚洲天堂国产精品一区在线| 亚洲第一区二区三区不卡| 国产精品一区二区在线观看99 | 免费看日本二区| 大又大粗又爽又黄少妇毛片口| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 久久久精品94久久精品| 国产女主播在线喷水免费视频网站 | 精品熟女少妇av免费看| 免费av毛片视频| 18禁在线播放成人免费| 青春草国产在线视频| 国产又黄又爽又无遮挡在线| 夜夜看夜夜爽夜夜摸| 三级国产精品欧美在线观看| 五月玫瑰六月丁香| 一本久久精品| 天堂影院成人在线观看| 亚洲成人av在线免费| 国产成人91sexporn| 狂野欧美激情性xxxx在线观看| 亚洲四区av| 高清毛片免费看| 中文字幕久久专区| 亚洲人成网站在线播| 国产不卡一卡二| 丰满人妻一区二区三区视频av| 国产视频内射| 免费看a级黄色片| 亚洲国产精品国产精品| 在现免费观看毛片| 五月玫瑰六月丁香| 一个人免费在线观看电影| 男人和女人高潮做爰伦理| 国产91av在线免费观看| 国产精品乱码一区二三区的特点| 久久99蜜桃精品久久| 丰满乱子伦码专区| 少妇熟女欧美另类| 十八禁国产超污无遮挡网站| 日本-黄色视频高清免费观看| 婷婷色综合大香蕉| 视频中文字幕在线观看| 国产中年淑女户外野战色| 午夜精品一区二区三区免费看| 色综合色国产| 国产精品福利在线免费观看| 听说在线观看完整版免费高清| 成人性生交大片免费视频hd| 夜夜爽夜夜爽视频| 精品不卡国产一区二区三区| 男女下面进入的视频免费午夜| 欧美激情在线99| 午夜福利成人在线免费观看| 国产免费男女视频| 国产精品.久久久| 久久欧美精品欧美久久欧美| www日本黄色视频网| a级毛片免费高清观看在线播放| 日本欧美国产在线视频| 国产成人a区在线观看| 欧美最新免费一区二区三区| 精品一区二区免费观看| 久久久久九九精品影院| 日韩视频在线欧美| 国产女主播在线喷水免费视频网站 | 国产在视频线精品| 国产精品麻豆人妻色哟哟久久 | 久久久精品94久久精品| 精品久久久久久久久亚洲| 青春草亚洲视频在线观看| 成人av在线播放网站| 一级爰片在线观看| 一级黄色大片毛片| 国产av在哪里看| 久久综合国产亚洲精品| 久久人人爽人人片av| 伦精品一区二区三区| 亚洲精品亚洲一区二区| 男女下面进入的视频免费午夜| 一二三四中文在线观看免费高清| 久久精品国产鲁丝片午夜精品| 九九在线视频观看精品| 中文字幕熟女人妻在线| 一级爰片在线观看| 男女下面进入的视频免费午夜| 午夜精品在线福利| 久久久久久久久久成人| 床上黄色一级片| 又爽又黄无遮挡网站| 最近中文字幕2019免费版| 亚洲欧洲日产国产| 国产在线一区二区三区精 | 少妇高潮的动态图| 人妻系列 视频| 最近的中文字幕免费完整| 一本一本综合久久| 简卡轻食公司| 欧美成人精品欧美一级黄| 色尼玛亚洲综合影院| 国产淫语在线视频| 天堂网av新在线| 色播亚洲综合网| 国产成人a∨麻豆精品| 欧美成人a在线观看| 欧美色视频一区免费| 国产麻豆成人av免费视频| 欧美3d第一页| 天堂av国产一区二区熟女人妻| 日产精品乱码卡一卡2卡三| 深爱激情五月婷婷| 亚洲人成网站在线播| 最近中文字幕2019免费版| 久久久久精品久久久久真实原创| 毛片女人毛片| 亚洲成av人片在线播放无| 欧美日本亚洲视频在线播放| 亚洲高清免费不卡视频| 大又大粗又爽又黄少妇毛片口| 插逼视频在线观看| 国产欧美日韩精品一区二区| 久久久久精品久久久久真实原创| 又粗又爽又猛毛片免费看| 网址你懂的国产日韩在线| 日韩精品青青久久久久久| 亚洲精品,欧美精品| 波多野结衣高清无吗| 日本与韩国留学比较| 久久久成人免费电影| 青春草亚洲视频在线观看| 一个人免费在线观看电影| 久久久精品欧美日韩精品| 国产精品一区二区性色av| 日本与韩国留学比较| 欧美日本亚洲视频在线播放| 亚洲怡红院男人天堂| 精品久久国产蜜桃| 深爱激情五月婷婷| 日韩亚洲欧美综合| 国产精品伦人一区二区| 日日啪夜夜撸| 国产三级在线视频| 日韩欧美精品v在线| 99久久中文字幕三级久久日本| 国国产精品蜜臀av免费| 亚洲欧美清纯卡通| 国产伦一二天堂av在线观看| 亚州av有码| 国产精品久久久久久精品电影小说 | 午夜激情欧美在线| 久久久久国产网址| 爱豆传媒免费全集在线观看| 免费观看的影片在线观看| 国产成人aa在线观看| 国产精品乱码一区二三区的特点| 亚洲一级一片aⅴ在线观看| 亚洲五月天丁香| 久久久精品大字幕| 国产三级在线视频| 一区二区三区乱码不卡18| 成人欧美大片| 久久99热这里只有精品18| 好男人在线观看高清免费视频| 免费av观看视频| 国产一区二区在线观看日韩| 国产片特级美女逼逼视频| 国产成人免费观看mmmm| 久久久久网色| www.av在线官网国产| 亚洲av免费在线观看| 亚洲精品国产av成人精品| 亚洲av电影不卡..在线观看| 精品久久久久久久久久久久久| av在线亚洲专区| 黄色一级大片看看| av.在线天堂| 精品人妻偷拍中文字幕| 亚洲精品成人久久久久久| 哪个播放器可以免费观看大片| 国产亚洲午夜精品一区二区久久 | 丰满少妇做爰视频| 赤兔流量卡办理| 国产精品无大码| 一本久久精品| 免费人成在线观看视频色| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 极品教师在线视频| av在线蜜桃| 狠狠狠狠99中文字幕| 国产免费又黄又爽又色| 亚洲欧美成人综合另类久久久 | 黄色日韩在线| 成年女人永久免费观看视频| 国产乱来视频区| 熟女电影av网| 欧美性猛交黑人性爽| 色尼玛亚洲综合影院| 永久免费av网站大全| 最近视频中文字幕2019在线8| 久久精品国产亚洲av天美| 日韩欧美 国产精品| 国产黄片视频在线免费观看| 亚洲人成网站高清观看| 国产精品国产三级专区第一集| 亚洲aⅴ乱码一区二区在线播放| 国产色爽女视频免费观看| 成人毛片60女人毛片免费| 丰满少妇做爰视频| 精品久久久久久久久久久久久| 久久精品久久精品一区二区三区| 免费观看人在逋| 嫩草影院新地址| 午夜福利视频1000在线观看| 国产老妇伦熟女老妇高清| 免费黄色在线免费观看| 蜜桃久久精品国产亚洲av| 中文字幕久久专区| 精品少妇黑人巨大在线播放 | 国产综合懂色| 午夜免费男女啪啪视频观看| 色5月婷婷丁香| 国产精品久久久久久久电影| 久久久午夜欧美精品| 亚洲av成人精品一二三区| 少妇的逼好多水| 午夜免费激情av| 麻豆精品久久久久久蜜桃| 亚洲人成网站在线播| 亚洲人成网站在线观看播放| 国产黄片美女视频| 亚洲av中文av极速乱| 日韩三级伦理在线观看| av在线播放精品| 国产精品一二三区在线看| 中文在线观看免费www的网站| 国产午夜精品一二区理论片| 黄色一级大片看看| 亚洲欧美精品自产自拍| 国内揄拍国产精品人妻在线| 97在线视频观看| 色播亚洲综合网| 日本黄色片子视频| 成人午夜精彩视频在线观看| 插逼视频在线观看| 97人妻精品一区二区三区麻豆| 美女国产视频在线观看| 国产亚洲91精品色在线| 亚洲av成人av| 亚洲av电影在线观看一区二区三区 | 国产免费视频播放在线视频 | 精品无人区乱码1区二区| 一级毛片久久久久久久久女| 欧美一区二区亚洲| 成人漫画全彩无遮挡| 久久久精品欧美日韩精品| av在线播放精品| 人妻少妇偷人精品九色| 毛片一级片免费看久久久久| 亚洲高清免费不卡视频| 嫩草影院入口| 国产片特级美女逼逼视频| 青春草国产在线视频| 男人舔女人下体高潮全视频| 久久婷婷人人爽人人干人人爱| 美女国产视频在线观看| 国产成人福利小说| 亚洲国产精品成人综合色| 两个人视频免费观看高清| 免费搜索国产男女视频| videossex国产| 亚洲精品国产av成人精品| 51国产日韩欧美| 麻豆成人午夜福利视频| 精品一区二区三区视频在线| 老司机影院成人| 精品人妻熟女av久视频| 国产淫片久久久久久久久| 精品一区二区三区视频在线| 99热这里只有是精品在线观看| 国产在视频线在精品| 日韩视频在线欧美| 99久久九九国产精品国产免费| 熟女电影av网| 欧美又色又爽又黄视频| 国产精品一及| 国产精品一及| 日本爱情动作片www.在线观看| 禁无遮挡网站| 啦啦啦啦在线视频资源| 99久久精品一区二区三区| 国产成人一区二区在线| 纵有疾风起免费观看全集完整版 | 禁无遮挡网站| 三级国产精品片| 老司机福利观看| 精品无人区乱码1区二区| 女人被狂操c到高潮| 国产精品久久久久久久久免| 亚洲在线自拍视频| 久久久欧美国产精品| 亚洲av二区三区四区| 亚洲自拍偷在线| 亚洲av成人精品一二三区| 97在线视频观看| 日日干狠狠操夜夜爽| 日韩亚洲欧美综合| 国产亚洲一区二区精品| 亚洲不卡免费看| a级一级毛片免费在线观看| 国产精品一区二区三区四区免费观看| 在线a可以看的网站| 午夜福利高清视频| 99久久精品热视频| 精品一区二区免费观看| 日韩强制内射视频| 一个人观看的视频www高清免费观看| 1024手机看黄色片| 麻豆av噜噜一区二区三区| 男人狂女人下面高潮的视频| 日本一本二区三区精品| 国产伦在线观看视频一区| 在线播放无遮挡| 欧美激情国产日韩精品一区| 亚洲经典国产精华液单| 边亲边吃奶的免费视频| 91午夜精品亚洲一区二区三区| 久久精品国产自在天天线| 久久99热这里只有精品18| 菩萨蛮人人尽说江南好唐韦庄 | 在线观看av片永久免费下载| 搡女人真爽免费视频火全软件| 国产亚洲91精品色在线| 亚洲av电影在线观看一区二区三区 | 日产精品乱码卡一卡2卡三| 国产精品av视频在线免费观看| 亚洲真实伦在线观看| 免费观看a级毛片全部| 99热网站在线观看| 亚洲av一区综合| 狂野欧美白嫩少妇大欣赏| 久久久a久久爽久久v久久| 国产探花极品一区二区| 久久精品人妻少妇| 亚洲国产欧美在线一区| 免费搜索国产男女视频| 精品无人区乱码1区二区| 久久精品国产鲁丝片午夜精品| 国产精品不卡视频一区二区| 特级一级黄色大片| 91久久精品电影网| 最新中文字幕久久久久| 九九热线精品视视频播放| 99久国产av精品国产电影| 一级黄片播放器| 成人亚洲精品av一区二区| 91精品一卡2卡3卡4卡| 国内精品美女久久久久久| 只有这里有精品99| 麻豆国产97在线/欧美| 一个人观看的视频www高清免费观看| 国产一区二区在线观看日韩| 亚洲真实伦在线观看| 黑人高潮一二区| 亚洲av熟女| 国产一级毛片在线| 日日摸夜夜添夜夜添av毛片| 午夜福利高清视频| 国产高清视频在线观看网站| 久久午夜福利片| 成人二区视频| 91精品一卡2卡3卡4卡| 国产爱豆传媒在线观看| 黑人高潮一二区| 亚洲精品国产av成人精品| 久久久久久久久久久丰满| 天天一区二区日本电影三级| 亚洲真实伦在线观看| 男女国产视频网站| 一级黄片播放器| 中文字幕久久专区| 美女脱内裤让男人舔精品视频| 日日啪夜夜撸| 淫秽高清视频在线观看| 国产单亲对白刺激| 少妇熟女aⅴ在线视频| 国产三级在线视频| 亚洲av熟女| 午夜精品国产一区二区电影 | av在线蜜桃| 一区二区三区乱码不卡18| 亚洲国产精品合色在线| 一个人免费在线观看电影| 国产午夜精品论理片| 中文字幕熟女人妻在线| 少妇的逼水好多| 久久久精品94久久精品| 国产精品电影一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 老师上课跳d突然被开到最大视频| 久久99热这里只有精品18| 久久精品国产鲁丝片午夜精品| 3wmmmm亚洲av在线观看| 极品教师在线视频| av播播在线观看一区| 国产中年淑女户外野战色| 中文字幕av成人在线电影| 久久久久网色| 两个人的视频大全免费| 久久国内精品自在自线图片| 国产av码专区亚洲av| 亚洲五月天丁香| 国产白丝娇喘喷水9色精品| 一夜夜www| 国产探花极品一区二区| 我要搜黄色片| 麻豆成人av视频| 一个人免费在线观看电影| 日本熟妇午夜| 青青草视频在线视频观看| 97超碰精品成人国产| 久久精品综合一区二区三区| 国产成人午夜福利电影在线观看| 国产亚洲最大av| 尾随美女入室| 狂野欧美激情性xxxx在线观看| 日韩欧美在线乱码| 我要看日韩黄色一级片| 看黄色毛片网站| 乱人视频在线观看| 久久久欧美国产精品| 丰满乱子伦码专区| 亚洲国产日韩欧美精品在线观看| av在线天堂中文字幕| 麻豆成人av视频| 91在线精品国自产拍蜜月| 日韩国内少妇激情av| 校园人妻丝袜中文字幕| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| www.av在线官网国产| 嫩草影院入口| 人体艺术视频欧美日本| 亚洲四区av| 1000部很黄的大片| 简卡轻食公司| 欧美又色又爽又黄视频| 超碰97精品在线观看| 少妇被粗大猛烈的视频| 色尼玛亚洲综合影院| 精品一区二区免费观看| 国产色爽女视频免费观看| АⅤ资源中文在线天堂| 亚洲人成网站在线观看播放| 99九九线精品视频在线观看视频| 久久精品国产亚洲av涩爱| 亚洲国产精品成人综合色| 国产精品国产三级国产av玫瑰| 久久人妻av系列| 国产一区亚洲一区在线观看| 国产欧美日韩精品一区二区| 男女边吃奶边做爰视频| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 成人午夜精彩视频在线观看| 最近视频中文字幕2019在线8| 亚洲性久久影院| 狠狠狠狠99中文字幕| 欧美最新免费一区二区三区| 国产亚洲精品久久久com| 午夜福利高清视频| 桃色一区二区三区在线观看| 亚洲成人精品中文字幕电影| 干丝袜人妻中文字幕| 国内揄拍国产精品人妻在线| 日韩欧美国产在线观看| 国产成人精品婷婷| 可以在线观看毛片的网站| 综合色av麻豆| 亚洲在久久综合| 嫩草影院入口| 国产精品久久久久久久电影| 亚洲国产精品成人久久小说| 日本五十路高清| 欧美另类亚洲清纯唯美| 日韩制服骚丝袜av| 真实男女啪啪啪动态图| 男女边吃奶边做爰视频| 亚洲性久久影院| 99久久无色码亚洲精品果冻| 国产av码专区亚洲av| 在线免费十八禁| 国产亚洲5aaaaa淫片| 成人亚洲精品av一区二区| 亚洲av.av天堂| 欧美潮喷喷水| 深爱激情五月婷婷| 精品少妇黑人巨大在线播放 | 国产白丝娇喘喷水9色精品| 春色校园在线视频观看| av在线蜜桃| 超碰97精品在线观看| 久久久久久伊人网av| 日韩一区二区视频免费看| 观看美女的网站| 99视频精品全部免费 在线| 国产精品一区www在线观看| 日韩欧美精品v在线| 国产久久久一区二区三区| 国产一区亚洲一区在线观看| 午夜免费男女啪啪视频观看| 免费观看性生交大片5| 熟女电影av网| 最近的中文字幕免费完整| 国产精品无大码| 99久久成人亚洲精品观看| 国产亚洲午夜精品一区二区久久 |