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

    隧道運(yùn)營期結(jié)構(gòu)健康評價及MATLAB應(yīng)用研究

    2016-07-04 07:04:48金煜皓王桂萱趙杰
    地震研究 2016年1期
    關(guān)鍵詞:模糊層次分析法

    金煜皓 王桂萱 趙杰

    摘要:針對模糊層次綜合評判方法計算的復(fù)雜性,使用MATLAB計算權(quán)重與模糊關(guān)系矩陣,并將兩者綜合得到模糊綜合評判矩陣。最后以實際運(yùn)營隧道為例,分析了模糊層次評價法的MATLAB程序在安全性評價中的應(yīng)用。結(jié)果表明,將模糊層次綜合評判法與MATLAB結(jié)合,能較好地進(jìn)行隧道運(yùn)營期結(jié)構(gòu)健康評價。

    關(guān)鍵詞:模糊層次分析法;MATLAB應(yīng)用;隧道健康評價

    中圖分類號:TU312 文獻(xiàn)標(biāo)識碼:A 文章編號:1000-0666(2016)01-0120-06

    0 引言

    隨著我國高速公路產(chǎn)業(yè)的快速發(fā)展,公路隧道的數(shù)量也呈直線增長態(tài)勢。據(jù)統(tǒng)計,我國已有公路隧道11345座,總長9606km,年均增長率高達(dá)20%,且有逐年增速加快的趨勢,已成為公路隧道工程建設(shè)規(guī)模和建設(shè)速度第一大國。然而,大部分已建隧道卻存在著混凝土強(qiáng)度低、襯砌開裂、襯砌背后空洞、襯砌厚度不足等病害,這些病害不僅會給隧道本身帶來嚴(yán)重的影響,而且有可能危及行車(人)的生命財產(chǎn)安全(王亞瓊等,2014)。因此,對運(yùn)營隧道結(jié)構(gòu)健康狀態(tài)的綜合評價就顯得十分重要。傳統(tǒng)的隧道結(jié)構(gòu)健康評價研究,大多仍停留在隧道施工期,且采用傳統(tǒng)的定性方法,這是因為隧道結(jié)構(gòu)性能受到多種因素的影響,具有很大的模糊性,無法用具體的數(shù)學(xué)公式準(zhǔn)確定量表達(dá)。近年來部分學(xué)者運(yùn)用建立數(shù)學(xué)模型的方法來描述運(yùn)營隧道結(jié)構(gòu)的安全性,洪平和劉鵬舉(2011)用層次分析法模型分析了鐵路運(yùn)營隧道健康情況,但卻無法消除層次分析法中特征值和特征向量復(fù)雜難求的缺陷,不利于工程應(yīng)用。李云等(2010)利用多級模糊評價算法對已建隧道襯砌健康進(jìn)行評價,但模糊評價模型沒有程序化,在面臨大量監(jiān)測數(shù)據(jù)時,往往需要進(jìn)行非常復(fù)雜的模糊計算,且結(jié)果精確度較差。

    綜上所述,已有的隧道運(yùn)營期結(jié)構(gòu)健康評價模型分析方法較為單一,僅有層次分析法或模糊評價法。另外,此類模型均沒有實現(xiàn)全面的程序化,即沒有建立一個全過程的結(jié)構(gòu)評判程序化模型,這就導(dǎo)致了計算過程效率十分低下,計算結(jié)果不穩(wěn)定且不準(zhǔn)確。因此,為了能更好的對隧道結(jié)構(gòu)健康狀態(tài)做出科學(xué)的、合理的評價,有必要選用合適的評價方法,并對該方法的全過程進(jìn)行程序化,以方便工程應(yīng)用。

    本文利用層次分析法求權(quán)重矩陣,模糊綜合評價法求模糊關(guān)系矩陣,并通過建立隧道運(yùn)營期結(jié)構(gòu)健康指標(biāo)體系及結(jié)構(gòu)安全等級劃分,將層次分析和模糊綜合評判法相結(jié)合,建立隧道運(yùn)營期結(jié)構(gòu)模糊層次健康評價模型,并將此評價模型MATLAB程序化(宋飛,趙法鎖,2008),最后以某實際公路隧道工程為例,說明模型的適用性與可靠性。

    1 模糊層次分析法及公路隧道運(yùn)營結(jié)構(gòu)安全評價體系

    1.1 層次分析法

    層次分析法(The Analytic Hierarchy Process,簡稱AHP)在20世紀(jì)70年代中期由美國運(yùn)籌學(xué)家Saaty正式提出。其分析步驟可分為建立層次結(jié)構(gòu)體系、構(gòu)造對比判斷矩陣、計算權(quán)向量,并做一致性檢驗(Harker,Vargas,1987)。

    層次分析法具有嚴(yán)格的上下層隸屬關(guān)系,可運(yùn)用1~9標(biāo)度法,進(jìn)行相對于上一層指標(biāo)重要性的兩兩比較,得判斷矩陣為

    A=(aijm×n. (1)式中:ni是第i個指標(biāo)相對第j個指標(biāo)的重要程度;矩陣A為互反矩陣,滿足aij>0,aij=1/aij,aii=1。對于層次結(jié)構(gòu),從最頂層到最底層有幾個隸屬關(guān)系,就要建立幾個指標(biāo)因素判斷矩陣。在建立了判斷矩陣之后,關(guān)鍵是求出判斷矩陣的最大特征值及該特征值對應(yīng)的特征向量,所得特征向量即為權(quán)重向量。如矩陣A的特征值可由下式求得:

    det(A-λI)=0. (2)式中:I為單位向量,λ為特征值,特征向量可由特征值進(jìn)一步求得。從數(shù)學(xué)的角度來說,隨著影響因素的增加,階數(shù)變高,判斷矩陣也會變得復(fù)雜,傳統(tǒng)特征向量求解不僅越來越困難,而且求出的近似解的精度很差。研究表明,MATLAB可用來精確快速求解復(fù)雜矩陣的特征值和特征向量。

    1.2 模糊綜合評判法

    模糊綜合評判法是模糊數(shù)學(xué)中最常見的方法之一(James et al,2001),運(yùn)用模糊及隸屬度理論,考慮了不同指標(biāo)對某一共同對象的綜合影響,據(jù)此,對該共同對象進(jìn)行等級劃分。當(dāng)某一事物被多種因素影響,但又無法精確給出每個影響因素對該事物的實際影響,就需要用到模糊綜合評判法,而大多數(shù)事物往往受到不確定的各種因素的共同影響,因此模糊綜合評判應(yīng)用十分廣泛。模糊綜合評價具體過程如下:將影響事物的多種因素看成一組模糊集合,稱之為因素集u,每個因素u都需要一組評判值進(jìn)行評判,故設(shè)立評判集v,而每個u對每個v都有相對應(yīng)且唯一的隸屬關(guān)系,此隸屬關(guān)系是模糊的,不是精確的,所有u對所有v的隸屬關(guān)系所組成的集合稱之為隸屬度矩陣(模糊關(guān)系矩陣R),將得到的模糊關(guān)系矩陣R與因素u所對應(yīng)的權(quán)重進(jìn)行模糊合成,就得到評價的定量解值,根據(jù)隸屬度最大原則,按照劃分好的等級,求出事物的等級。上述即為模糊綜合評價過程,其優(yōu)點(diǎn)是將不確定的、模糊的關(guān)系定量化,最大程度減少人為和定性所造成的誤差。其缺點(diǎn)是在對權(quán)重指標(biāo)的確定過程中,主觀性較強(qiáng),針對模糊綜合評價求權(quán)重的缺點(diǎn),可采用層次分析法求權(quán)重,并與模糊綜合評價法相結(jié)合,通過模糊矩陣合成,求出事物等級。

    1.3 隧道運(yùn)營安全評價體系

    參考目前國內(nèi)外已有的研究成果(戴憂華等,2010)和《公路隧道養(yǎng)護(hù)技術(shù)規(guī)范》(JTG H12-2003)以及指標(biāo)選取原則,同時考慮本文的研究范圍,選取僅考慮無特殊地質(zhì)條件,且無深埋的公路隧道洞身標(biāo)準(zhǔn)斷面襯砌結(jié)構(gòu)為研究對象。結(jié)合層次分析法和模糊綜合評價理論,建立公路隧道運(yùn)營期結(jié)構(gòu)安全性評價體系,如圖1所示。由圖1可知,隧道主體結(jié)構(gòu)模糊層次評價模型分為3層:第一層為目標(biāo)層,即隧道運(yùn)營期結(jié)構(gòu)安全狀態(tài)A;第二層為準(zhǔn)則層,因隧道運(yùn)營期結(jié)構(gòu)安全主要受襯砌影響,根據(jù)重點(diǎn)突出,少而精且互不影響原則,可分為襯砌損害A1,襯砌混凝土特性A2,襯砌背后空洞A3,這3個因素組成了隧道結(jié)構(gòu)安全評價體系的制約因素集;第三層為指標(biāo)層,反映了影響上一層即準(zhǔn)則層的最可能的影響因素,各指標(biāo)為:裂縫寬度A11變形量A12碳化深度A13,混凝土厚度指標(biāo)A21、混凝土強(qiáng)度指標(biāo)A22,襯砌空洞直徑A31。

    2 公路隧道運(yùn)營安全等級劃分

    常見的幾種安全等級劃分法中,三級劃分法過于簡單,五級劃分法實質(zhì)是在四級劃分法基礎(chǔ)上對某項進(jìn)行了細(xì)化,十級劃分法又劃分過細(xì),很多狀態(tài)是不需要的,因此筆者參考已有運(yùn)營隧道結(jié)構(gòu)安全等級劃分情況,綜合分析各種劃分法的優(yōu)缺點(diǎn),選用四級劃分法作為隧道運(yùn)營期結(jié)構(gòu)安全性等級劃分的依據(jù):Ⅰ級為隧道結(jié)構(gòu)非常安全,襯砌無損傷;Ⅱ級為隧道基本安全,襯砌輕微受損;Ⅲ級為隧道結(jié)構(gòu)不安全,襯砌存在裂縫,背后出現(xiàn)空洞;Ⅳ級為隧道結(jié)構(gòu)嚴(yán)重破壞,襯砌大面積掉塊,危及行車安全。

    3 模糊層次綜合評價模型的MATLAB程序?qū)崿F(xiàn)

    3.1 特征向量及一致性檢驗的MATLAB實現(xiàn)

    MATLAB可以精確快速地求解層次分析法中復(fù)雜矩陣特征值與特征向量,可調(diào)用eig函數(shù)進(jìn)行求解。MATLAB中具體編程如下,定義M函數(shù)MAX-eigvalvec.m:

    function[eigval,w]=MAX-eigvalvec(A);%求最大特征值和歸一化的特征向量,A為判斷矩陣;

    [eigvec,eigval]=eig(A);

    maxeigval=max(eigval);%最大特征值;

    v=eigvec(:,index);%最大特征值對應(yīng)的特征向量;

    w=v./sum(v);%特征向量歸一化;

    end

    層次分析法要求必須對各層次間的判斷矩陣進(jìn)行一致性檢驗。這是為了避免同一層中A比B重要,B比C重要,而C卻比A重要的錯誤出現(xiàn)。定義s-examine.m函數(shù)實現(xiàn)一致性檢驗:

    function[RI,CI]=s-examine(maxeigval,A);%單排序一致性檢驗;

    n=size(A,1);

    RIT=[0.0 0.0 0.58 0.90 1.12 1.24 1.32 1.411.45 1.49 1.51 1.54 1.56 1.58 1.59];

    RI=RIT(n);%隨機(jī)一致性指標(biāo);

    CI=(maxeigval-n)/(n-1);%一致性指標(biāo);

    CR=CI/RI;

    if CR>=0.1

    disp([input(‘矩陣沒有通過一致性檢驗,請重新調(diào)整判斷矩陣)]);

    else

    disp([input(‘矩陣通過一致性檢驗)]);

    end

    3.2 MATLAB模糊綜合評判模型

    定義模糊綜合評價M函數(shù)Fuzzy_F.m,在MATLAB實現(xiàn)如下:

    function R=Fuzzy_F(A1,A2,~Am);%A1~Am為指標(biāo)層實際監(jiān)測值;

    R=zeros(m,n);%m行,n列空矩陣,作為所求模糊關(guān)系矩陣R的存放空間;

    S=[S11S12…S12S1n;S21S22…S12S2n;Sm1Sm2…Smn];%評價標(biāo)準(zhǔn),即各影響因素的隸屬函數(shù)參數(shù),可結(jié)合評判基準(zhǔn)求得;

    T=zeios(1,i);

    T(1,1)=A1;T(1,2)=A2;T(1,3)=A3;~T(1,m)=Am;%將實際影響因素值賦予T中;

    for i=1:m;

    R(i,1)=fi1(T(i),S(i,)..);

    R(i,2)=fi2(T(i),S(i,)..);

    R(i,3)=fi3(T(i),S(i,)..);

    R(i,4)=fi4(T(i),S(i,)..);

    end

    上述程序是將T(1,i)代入隸屬度函數(shù),計算隸屬度矩陣(模糊關(guān)系矩陣),這4個隸屬度函數(shù),對應(yīng)所劃分的4個等級,該程序?qū)崿F(xiàn)了所列影響因素實際值分別隸屬Ⅰ、Ⅱ、Ⅲ和Ⅳ四個安全等級的隸屬程度,即隸屬度矩陣Ro而隸屬函數(shù)構(gòu)造的好壞直接決定評價系統(tǒng)的優(yōu)劣,一般而言,隸屬函數(shù)有模糊統(tǒng)計方法、三分法、模糊分步法、專家評分法等,考慮到統(tǒng)計方法需要較多樣本,專家評分法主觀性較強(qiáng),故隸屬函數(shù)常常采用模糊分布法。遵循隸屬函數(shù)確定的原則和方法,結(jié)合蘇永華(2007)關(guān)于模糊隸屬函數(shù)的構(gòu)造方法研究以及對各分布模型進(jìn)行對比分析,本文選取模糊分布法中的“降半梯形”作為指標(biāo)層各指標(biāo)的隸屬函數(shù)。

    3.3 模糊層次評判結(jié)果MATLAB實現(xiàn)

    經(jīng)綜合比較分析,本文采用加權(quán)平均模糊算子,進(jìn)而將層次分析法所得權(quán)重矩陣與模糊綜合評價所得隸屬矩陣進(jìn)行加權(quán)模糊結(jié)合,選用加權(quán)平均模糊算子,即先將兩個矩陣?yán)锏臄?shù)對應(yīng)相乘后相加M(·,①),MATLAB程序?qū)崿F(xiàn)如下,定義M函數(shù)jq(w,R):

    function WR=jq(w,R);

    m=size(W,1);n=size(R,2);

    for i-1:W

    for j=1:R

    WR(i,j)=sum(w(i,:).*R(:,j));

    end

    end

    B=WR./sum(WR);%對矩陣歸一化,B為模糊合成矩陣,由最大隸屬原則,確定模型等級。

    4 實例分析

    某運(yùn)營10年隧道為雙向二車道隧道,左線起訖樁號為K195+620~K197+390,全長1700m;右線起訖樁號為YK195+619~YK197+369,全長1740m。本隧道通過地段地層結(jié)構(gòu)簡單,地表均為第四系全新統(tǒng)(Q4)殘坡積層所覆蓋,隧道穿越處最高標(biāo)高為360m。該隧道由于運(yùn)營時間較長,其襯砌結(jié)構(gòu)已出現(xiàn)不同程度的損壞現(xiàn)象,因此有必要對其結(jié)構(gòu)健康進(jìn)行安全性評價。

    另外,由于國內(nèi)外對隧道運(yùn)營期結(jié)構(gòu)健康評估還處于起步階段,更缺少完整的程序化的評價模型,因而文中對該隧道運(yùn)營期結(jié)構(gòu)安全進(jìn)行基于程序化模型的分析,說明模糊層次分析在隧道運(yùn)營期結(jié)構(gòu)健康評價中的MATLAB應(yīng)用。

    4.1 權(quán)重設(shè)計及模糊隸屬度矩陣計算

    針對層次分析法求權(quán)重中的標(biāo)度選用問題,駱正清和楊善林(2004)對常用的幾種層次分析法標(biāo)度進(jìn)行比較,得出對單一準(zhǔn)則下的排序,1~9標(biāo)度較好;對多準(zhǔn)則下的排序,指數(shù)標(biāo)度更好。綜上分析,依據(jù)圖1的安全層次結(jié)構(gòu),參考工程經(jīng)驗,利用1~9標(biāo)度法,對監(jiān)測到的該隧道各個具體影響因素實際數(shù)值進(jìn)行兩兩比較評分,得判斷矩陣,如表1~3所示。

    對背后空洞直徑指標(biāo),因為其只有一項,故A31=[1]。

    對準(zhǔn)則層指標(biāo)進(jìn)行兩兩比較評分,得判斷矩陣,如表3所示。

    根據(jù)上述判斷矩陣,在MATLAB的CommandWindow中鍵入如下命令,可迅速得出各層的權(quán)重矩陣:

    A=[131;1/311/3;131];

    [maxA,wA]=maxeigvalvec(A);

    可得

    [ωA]=[0.4286 0.1429 0.4286];

    同理可得

    [ωA1]=[0.1429 0.7143 0.1429];

    [ωA2]=[0.1667 0.8333];

    [ωA3]=[1].

    基于MATLAB程序求模糊隸屬度矩陣,其中,隧道運(yùn)營期結(jié)構(gòu)安全隸屬函數(shù)參數(shù)為S=[0.050.217 0.05 0.217 0.417 0.217 0.417 0.75 0.417O.75:0.05 0.217 0.05 0.217 0.417 0.217 0.4170.75 0.417 0.75;0.5 2 0.5 2 4 2 4 6 4 6;050 0 50 300 50 300 600 300 600:0.5 2 0.5 26.5 2 6.5 12 6.5 12:0.5 2 0.5 2 3.5 2 3.54.5 3.5 4.5],進(jìn)而根據(jù)現(xiàn)場傳感器采集到的監(jiān)測數(shù)據(jù),并調(diào)用本文所用隸屬度函數(shù)及Fuzzy_F.m函數(shù),在Command Window中鍵入以下代碼:

    R=Fuzzy_FA1(0.5,0.5,0.5);

    得A1所包含的指標(biāo)層影響因素對應(yīng)的隸屬度矩陣:

    將現(xiàn)場監(jiān)測數(shù)據(jù)A21=0.05,A22=0.05,A31=0代入相應(yīng)Fuzzy_F.m函數(shù),得:

    4.2 一級模糊綜合評判

    將由MATLAB編程得到的權(quán)重向量W與模糊關(guān)系矩陣R進(jìn)行加權(quán)平均模糊合成,在CommandWindow中鍵入如下命令:

    C(A1)=jq(WA1,RA1);該語句實現(xiàn)了WA1與RA1加權(quán)模糊合成;

    B(A1)=C(A1)./sum(C(A1));表示對模糊合成的矩陣進(jìn)行歸一化;

    點(diǎn)擊回車鍵,得:

    B(A1)=[0.142 9 0 0.643 5 0.213 6].

    同理可得:

    B(A2)=[1 0 0 0];

    B(A3)=[1 0 0 0].

    綜上可得一級模糊層次綜合評判矩陣為

    4.3 二級模糊綜合評判

    A表示其所包含的準(zhǔn)則層的3個影響因素,將WA與RA進(jìn)行加權(quán)模糊合成:

    C(A)=jq(WA,RA);

    對模糊合成的矩陣進(jìn)行歸一化:

    B(A)=C(A)./sum(C(A));得

    B(A)=[0.6327 0 0.2758 0.0915].

    由以上二級模糊綜合評判計算結(jié)果,并根據(jù)最大隸屬原則,0.6327最大,所以運(yùn)營隧道結(jié)構(gòu)處于Ⅰ級狀態(tài),即結(jié)構(gòu)很安全,這與該運(yùn)營隧道實際情況一致,但相對而言,0.2758(Ⅲ級)也較大,說明隧道安全還是存在風(fēng)險,再加上隧道襯砌結(jié)構(gòu)隨著時間的增長,影響其安全的各個指標(biāo)也會有所下降,必然會對運(yùn)營隧道結(jié)構(gòu)安全帶來威脅,所以應(yīng)當(dāng)采取合適的安全措施。

    5 結(jié)論

    針對運(yùn)營隧道結(jié)構(gòu)健康評價問題,建立基于模糊層次綜合評價法的程序模型,通過分析得出如下結(jié)論:

    (1)將層次分析法和模糊評價法相結(jié)合,建立隧道運(yùn)營期結(jié)構(gòu)安全模糊層次評價模型,克服了評價過程中主觀隨意性的缺點(diǎn),選用“降半梯形”作為指標(biāo)層各指標(biāo)的隸屬函數(shù),使得評價結(jié)果更符合客觀實際。

    (2)利用MATLAB語言實現(xiàn)了評價模型建立全過程的程序化,程序化后模型不僅可以極大縮短傳統(tǒng)模型計算時間,也顯著提高了計算結(jié)果的精度和準(zhǔn)確度,為隧道運(yùn)營期結(jié)構(gòu)健康評價系統(tǒng)開發(fā)提供所需的應(yīng)用程序。

    (3)以實際隧道運(yùn)營期結(jié)構(gòu)健康評價為例,結(jié)果表明,將模糊層次綜合評判法結(jié)合MATLAB應(yīng)用,能較好地進(jìn)行隧道運(yùn)營期結(jié)構(gòu)健康評價,是一種開展運(yùn)營隧道結(jié)構(gòu)安全風(fēng)險分析的有效方法。

    猜你喜歡
    模糊層次分析法
    美國經(jīng)濟(jì)金融化測度研究
    基于業(yè)務(wù)優(yōu)先級的QoS控制技術(shù)研究
    公私合營模式下港口工程建設(shè)風(fēng)險評價
    綠色平衡計分卡
    生鮮農(nóng)產(chǎn)品供應(yīng)鏈績效評價指標(biāo)體系構(gòu)建研究
    基于模糊層次分析法的房地產(chǎn)項目風(fēng)險研究
    商(2016年13期)2016-05-20 10:46:33
    基于模糊層次分析法的管制員班組資源管理分析
    科技視界(2016年7期)2016-04-01 20:00:57
    民進(jìn)軍企業(yè)綜合能力評價體系探究
    QFD與FAHP在軟件開發(fā)項目中的應(yīng)用
    房地產(chǎn)電子商務(wù)交易風(fēng)險評價研究
    亚洲熟女精品中文字幕| 这个男人来自地球电影免费观看 | 免费观看性生交大片5| 精品少妇久久久久久888优播| 校园人妻丝袜中文字幕| 久久热在线av| 综合色丁香网| 亚洲成人手机| 亚洲国产最新在线播放| 精品卡一卡二卡四卡免费| 国产精品久久久久成人av| 99久国产av精品国产电影| 丝袜脚勾引网站| 精品少妇一区二区三区视频日本电影 | 久久国产亚洲av麻豆专区| 嫩草影院入口| 2021少妇久久久久久久久久久| 最近2019中文字幕mv第一页| 日韩一区二区三区影片| 妹子高潮喷水视频| 在线亚洲精品国产二区图片欧美| 亚洲精品一区蜜桃| 免费少妇av软件| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜一区二区 | 啦啦啦在线免费观看视频4| 亚洲婷婷狠狠爱综合网| av不卡在线播放| 亚洲欧美一区二区三区国产| 一区二区三区激情视频| 亚洲国产欧美日韩在线播放| 自线自在国产av| 一级爰片在线观看| 午夜福利乱码中文字幕| 亚洲精华国产精华液的使用体验| 国产探花极品一区二区| 成人亚洲精品一区在线观看| 久久人妻熟女aⅴ| 久久免费观看电影| 日本午夜av视频| 桃花免费在线播放| 成年女人毛片免费观看观看9 | 赤兔流量卡办理| 91精品三级在线观看| 亚洲人成77777在线视频| 国产福利在线免费观看视频| 日韩av在线免费看完整版不卡| 欧美人与性动交α欧美精品济南到| 国产精品秋霞免费鲁丝片| 777米奇影视久久| av天堂久久9| 亚洲精品久久午夜乱码| 99久久99久久久精品蜜桃| 丝袜脚勾引网站| 2021少妇久久久久久久久久久| 各种免费的搞黄视频| www.熟女人妻精品国产| 嫩草影院入口| 97精品久久久久久久久久精品| 超碰97精品在线观看| 丁香六月天网| 亚洲精品久久成人aⅴ小说| 亚洲男人天堂网一区| 黑丝袜美女国产一区| 国产精品av久久久久免费| 国产精品久久久久久人妻精品电影 | 亚洲av日韩在线播放| 亚洲美女黄色视频免费看| 桃花免费在线播放| 成人漫画全彩无遮挡| 美国免费a级毛片| 妹子高潮喷水视频| 极品少妇高潮喷水抽搐| 秋霞伦理黄片| 十八禁网站网址无遮挡| 欧美日韩一区二区视频在线观看视频在线| 午夜日本视频在线| 日韩制服骚丝袜av| 国产精品 欧美亚洲| 久久国产亚洲av麻豆专区| 美女福利国产在线| 97人妻天天添夜夜摸| 久久久久久久大尺度免费视频| 新久久久久国产一级毛片| 熟女av电影| 老司机影院毛片| 天天躁夜夜躁狠狠久久av| 在线观看免费日韩欧美大片| 久久精品久久久久久噜噜老黄| av网站在线播放免费| 成人漫画全彩无遮挡| 夜夜骑夜夜射夜夜干| 激情视频va一区二区三区| 丝袜喷水一区| 久久鲁丝午夜福利片| 亚洲精品aⅴ在线观看| 日本vs欧美在线观看视频| 亚洲欧洲日产国产| 91精品国产国语对白视频| 9热在线视频观看99| 韩国精品一区二区三区| 青春草视频在线免费观看| 亚洲国产日韩一区二区| av网站免费在线观看视频| 成人国产麻豆网| 国产亚洲午夜精品一区二区久久| 久久 成人 亚洲| 亚洲五月色婷婷综合| 亚洲精品一区蜜桃| 国产精品国产av在线观看| 色网站视频免费| 这个男人来自地球电影免费观看 | 日韩制服丝袜自拍偷拍| avwww免费| 天天添夜夜摸| 在线观看免费日韩欧美大片| 性少妇av在线| 成人漫画全彩无遮挡| 老司机深夜福利视频在线观看 | 国产97色在线日韩免费| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美一区二区三区久久| 国产黄色视频一区二区在线观看| 乱人伦中国视频| 免费看不卡的av| 人体艺术视频欧美日本| 成人黄色视频免费在线看| av网站在线播放免费| 在线观看免费高清a一片| 秋霞在线观看毛片| 亚洲成人免费av在线播放| 人人妻人人爽人人添夜夜欢视频| 免费女性裸体啪啪无遮挡网站| 国产熟女午夜一区二区三区| 亚洲欧洲日产国产| 咕卡用的链子| 老司机靠b影院| 久久久欧美国产精品| 欧美日韩综合久久久久久| 国产黄色免费在线视频| 各种免费的搞黄视频| xxxhd国产人妻xxx| 曰老女人黄片| www.自偷自拍.com| 无遮挡黄片免费观看| 精品人妻熟女毛片av久久网站| 亚洲av欧美aⅴ国产| 亚洲美女黄色视频免费看| 桃花免费在线播放| 中文字幕最新亚洲高清| 午夜福利影视在线免费观看| 亚洲三区欧美一区| 一个人免费看片子| 日韩,欧美,国产一区二区三区| 久久精品国产亚洲av高清一级| 亚洲天堂av无毛| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 欧美精品亚洲一区二区| 丝袜喷水一区| 性高湖久久久久久久久免费观看| 亚洲精品视频女| 色吧在线观看| 日韩成人av中文字幕在线观看| 少妇 在线观看| 欧美av亚洲av综合av国产av | 国产精品一区二区在线不卡| 欧美 亚洲 国产 日韩一| 丝瓜视频免费看黄片| 高清视频免费观看一区二区| 久久久久久久久久久久大奶| 精品一区二区免费观看| 久久99热这里只频精品6学生| 人妻一区二区av| 国产一区二区三区综合在线观看| 午夜福利网站1000一区二区三区| 精品久久久久久电影网| 校园人妻丝袜中文字幕| 少妇 在线观看| 日韩视频在线欧美| 国产有黄有色有爽视频| 亚洲欧洲国产日韩| 久久热在线av| 亚洲欧美精品综合一区二区三区| av在线老鸭窝| videosex国产| 一区在线观看完整版| 国产精品一区二区在线不卡| 老司机在亚洲福利影院| 无遮挡黄片免费观看| 人人妻人人澡人人爽人人夜夜| 亚洲一卡2卡3卡4卡5卡精品中文| www.熟女人妻精品国产| 老司机影院成人| videos熟女内射| 久久精品人人爽人人爽视色| 欧美人与性动交α欧美精品济南到| 青春草亚洲视频在线观看| 一边摸一边抽搐一进一出视频| 国产亚洲精品第一综合不卡| 人人妻人人澡人人爽人人夜夜| 日韩伦理黄色片| 久久久久久久久免费视频了| 免费黄色在线免费观看| 欧美黑人精品巨大| 免费观看人在逋| 精品一区二区三区四区五区乱码 | 日韩不卡一区二区三区视频在线| 两个人看的免费小视频| 曰老女人黄片| 国产淫语在线视频| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美成人综合另类久久久| 人人妻人人添人人爽欧美一区卜| 免费黄色在线免费观看| 一区二区日韩欧美中文字幕| 满18在线观看网站| 国产一区亚洲一区在线观看| 热re99久久国产66热| 韩国精品一区二区三区| 国产在线免费精品| 观看av在线不卡| av卡一久久| www.自偷自拍.com| 大陆偷拍与自拍| 久久这里只有精品19| 久久国产亚洲av麻豆专区| 国产一级毛片在线| 欧美黑人欧美精品刺激| 久久性视频一级片| 在线观看国产h片| 亚洲精品第二区| 1024香蕉在线观看| 亚洲精品久久久久久婷婷小说| 欧美另类一区| 亚洲欧美一区二区三区国产| 日韩一区二区视频免费看| 久久久亚洲精品成人影院| 久久久精品94久久精品| 最近2019中文字幕mv第一页| 久久精品aⅴ一区二区三区四区| 亚洲av欧美aⅴ国产| 日韩av免费高清视频| 女性被躁到高潮视频| 精品国产一区二区三区四区第35| 男女下面插进去视频免费观看| 亚洲免费av在线视频| 久久精品国产综合久久久| 精品人妻一区二区三区麻豆| 久久国产精品大桥未久av| 色婷婷av一区二区三区视频| 亚洲国产欧美网| 人人妻,人人澡人人爽秒播 | h视频一区二区三区| 亚洲一码二码三码区别大吗| 国产一区亚洲一区在线观看| 超碰成人久久| 亚洲国产看品久久| 精品免费久久久久久久清纯 | 天天操日日干夜夜撸| 大码成人一级视频| 午夜福利影视在线免费观看| 国产精品麻豆人妻色哟哟久久| 久久精品国产亚洲av高清一级| 久久久国产一区二区| 日韩 亚洲 欧美在线| 国产伦理片在线播放av一区| 亚洲精品,欧美精品| 精品一区二区三区av网在线观看 | 狠狠精品人妻久久久久久综合| 久久久久久久久久久免费av| 色婷婷av一区二区三区视频| 一本大道久久a久久精品| 美女午夜性视频免费| 亚洲av欧美aⅴ国产| 亚洲精品国产色婷婷电影| 国产日韩欧美视频二区| 亚洲国产精品999| 日本午夜av视频| av卡一久久| 精品亚洲成国产av| 成人毛片60女人毛片免费| 黑人猛操日本美女一级片| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 又大又黄又爽视频免费| 美国免费a级毛片| 免费少妇av软件| 久久免费观看电影| 黑人欧美特级aaaaaa片| 国产深夜福利视频在线观看| 国产有黄有色有爽视频| 一边亲一边摸免费视频| 久热爱精品视频在线9| 亚洲,欧美精品.| 99热国产这里只有精品6| 黄色怎么调成土黄色| 午夜激情av网站| 亚洲欧美成人精品一区二区| 成人黄色视频免费在线看| 热re99久久国产66热| 亚洲成人免费av在线播放| 在线看a的网站| 欧美国产精品一级二级三级| 在线免费观看不下载黄p国产| 男女床上黄色一级片免费看| 又大又黄又爽视频免费| 亚洲精品国产av成人精品| 在线观看免费视频网站a站| 久久99精品国语久久久| 亚洲 欧美一区二区三区| 日韩制服丝袜自拍偷拍| 满18在线观看网站| 久久影院123| 91成人精品电影| 国产成人精品无人区| 欧美精品av麻豆av| 悠悠久久av| 精品国产露脸久久av麻豆| 国产成人啪精品午夜网站| 啦啦啦啦在线视频资源| 国产成人精品福利久久| 色吧在线观看| 一二三四在线观看免费中文在| 99re6热这里在线精品视频| 国产麻豆69| 美女国产高潮福利片在线看| 深夜精品福利| 菩萨蛮人人尽说江南好唐韦庄| 国产老妇伦熟女老妇高清| 亚洲欧美成人精品一区二区| 国产极品粉嫩免费观看在线| 亚洲欧美清纯卡通| 久热这里只有精品99| 国产一区二区三区av在线| 久久精品国产综合久久久| 日韩精品免费视频一区二区三区| 街头女战士在线观看网站| 9191精品国产免费久久| 日韩 欧美 亚洲 中文字幕| 国产精品无大码| 午夜福利影视在线免费观看| 考比视频在线观看| 黄色一级大片看看| 国产成人欧美在线观看 | 亚洲欧美成人综合另类久久久| 久久精品国产综合久久久| 日韩欧美一区视频在线观看| 日本色播在线视频| 欧美日韩国产mv在线观看视频| 国产淫语在线视频| 亚洲天堂av无毛| 不卡av一区二区三区| 亚洲国产精品999| 少妇人妻久久综合中文| 日本欧美国产在线视频| 无限看片的www在线观看| 欧美日本中文国产一区发布| 亚洲av男天堂| 国产乱来视频区| 国产男女超爽视频在线观看| 欧美另类一区| 嫩草影视91久久| 国产免费福利视频在线观看| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 高清黄色对白视频在线免费看| 精品少妇久久久久久888优播| tube8黄色片| 黄片无遮挡物在线观看| av电影中文网址| 亚洲第一青青草原| 色吧在线观看| 久久久精品区二区三区| 欧美激情高清一区二区三区 | 国产精品久久久av美女十八| 这个男人来自地球电影免费观看 | 国产成人精品在线电影| 一级毛片电影观看| 日韩制服丝袜自拍偷拍| 美女脱内裤让男人舔精品视频| 少妇人妻 视频| 欧美精品亚洲一区二区| 午夜福利乱码中文字幕| 精品福利永久在线观看| 中文字幕av电影在线播放| 国产又爽黄色视频| 日本午夜av视频| 一本久久精品| 亚洲第一青青草原| 桃花免费在线播放| 丰满饥渴人妻一区二区三| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品一区三区| 无遮挡黄片免费观看| 91aial.com中文字幕在线观看| av不卡在线播放| 国产一区二区激情短视频 | 久久久国产精品麻豆| 欧美少妇被猛烈插入视频| www.自偷自拍.com| 欧美亚洲 丝袜 人妻 在线| 无限看片的www在线观看| 精品少妇黑人巨大在线播放| av又黄又爽大尺度在线免费看| 老司机靠b影院| 亚洲av欧美aⅴ国产| 精品卡一卡二卡四卡免费| 午夜免费男女啪啪视频观看| 七月丁香在线播放| 欧美成人精品欧美一级黄| 电影成人av| 日韩欧美精品免费久久| 99国产精品免费福利视频| 777米奇影视久久| 熟妇人妻不卡中文字幕| 1024视频免费在线观看| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久av网站| 又粗又硬又长又爽又黄的视频| 国产爽快片一区二区三区| 日韩av不卡免费在线播放| 人妻人人澡人人爽人人| 99re6热这里在线精品视频| 男女边吃奶边做爰视频| 国产亚洲精品第一综合不卡| 哪个播放器可以免费观看大片| 狂野欧美激情性bbbbbb| 国产成人a∨麻豆精品| 精品一区二区三卡| 成人国产麻豆网| 青春草国产在线视频| 制服丝袜香蕉在线| 美女福利国产在线| 尾随美女入室| 亚洲精品国产一区二区精华液| 自线自在国产av| 街头女战士在线观看网站| av女优亚洲男人天堂| 久久久国产精品麻豆| 人人澡人人妻人| 一区二区日韩欧美中文字幕| 一本—道久久a久久精品蜜桃钙片| 欧美成人精品欧美一级黄| 色精品久久人妻99蜜桃| 亚洲欧美清纯卡通| 国产成人精品福利久久| 久久午夜综合久久蜜桃| 久久国产精品男人的天堂亚洲| 国产极品粉嫩免费观看在线| 亚洲欧美激情在线| 欧美中文综合在线视频| 欧美精品亚洲一区二区| 制服丝袜香蕉在线| 热re99久久国产66热| 久久精品国产亚洲av涩爱| www.熟女人妻精品国产| 国产成人系列免费观看| 日韩精品有码人妻一区| 青春草国产在线视频| 麻豆av在线久日| 亚洲一区中文字幕在线| 国产有黄有色有爽视频| 丁香六月欧美| 18禁国产床啪视频网站| 精品午夜福利在线看| 无遮挡黄片免费观看| 久久这里只有精品19| 女人被躁到高潮嗷嗷叫费观| 亚洲激情五月婷婷啪啪| 国产精品熟女久久久久浪| 欧美最新免费一区二区三区| 最近中文字幕高清免费大全6| 国产精品国产av在线观看| 啦啦啦在线观看免费高清www| 91成人精品电影| 久久久久精品久久久久真实原创| 蜜桃国产av成人99| 亚洲精品视频女| 丝袜脚勾引网站| 高清黄色对白视频在线免费看| 99九九在线精品视频| 黄色一级大片看看| 天美传媒精品一区二区| 国产成人午夜福利电影在线观看| 一边亲一边摸免费视频| 巨乳人妻的诱惑在线观看| 成人国产av品久久久| 中文字幕亚洲精品专区| 亚洲人成电影观看| 97人妻天天添夜夜摸| 国产熟女欧美一区二区| 久久97久久精品| 日韩,欧美,国产一区二区三区| 9色porny在线观看| 18禁国产床啪视频网站| 伦理电影免费视频| 国产淫语在线视频| 国产亚洲欧美精品永久| 天天影视国产精品| 亚洲伊人久久精品综合| 丝瓜视频免费看黄片| 人人妻人人澡人人看| 中文字幕制服av| 日韩伦理黄色片| 婷婷成人精品国产| 街头女战士在线观看网站| 亚洲三区欧美一区| 国产黄频视频在线观看| 精品少妇一区二区三区视频日本电影 | 天堂中文最新版在线下载| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 日韩欧美精品免费久久| 国产成人精品在线电影| 嫩草影院入口| 成人亚洲欧美一区二区av| 99久久99久久久精品蜜桃| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 日韩熟女老妇一区二区性免费视频| 国产av国产精品国产| 国产成人精品福利久久| 毛片一级片免费看久久久久| 女人久久www免费人成看片| 久久精品久久精品一区二区三区| 亚洲男人天堂网一区| 久久人人爽av亚洲精品天堂| 一本一本久久a久久精品综合妖精| 成人18禁高潮啪啪吃奶动态图| 亚洲第一青青草原| 黑人欧美特级aaaaaa片| 免费在线观看完整版高清| 国产精品女同一区二区软件| 99re6热这里在线精品视频| 午夜免费鲁丝| 捣出白浆h1v1| 国产成人a∨麻豆精品| 90打野战视频偷拍视频| 嫩草影院入口| 叶爱在线成人免费视频播放| 精品久久久精品久久久| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 成人毛片60女人毛片免费| 97人妻天天添夜夜摸| 一级片免费观看大全| 亚洲精品第二区| 久久性视频一级片| 美女视频免费永久观看网站| 欧美日韩福利视频一区二区| www.自偷自拍.com| 色婷婷av一区二区三区视频| 国产亚洲av片在线观看秒播厂| 国产精品人妻久久久影院| avwww免费| 国产视频首页在线观看| 亚洲成人国产一区在线观看 | 亚洲精品av麻豆狂野| 亚洲人成77777在线视频| 一区二区三区乱码不卡18| 天堂中文最新版在线下载| 国产精品久久久久久精品电影小说| av网站免费在线观看视频| 久久人人爽人人片av| 婷婷色麻豆天堂久久| 日韩一区二区三区影片| 国产 一区精品| 亚洲国产成人一精品久久久| 精品一区二区免费观看| 一个人免费看片子| 午夜久久久在线观看| 久久亚洲国产成人精品v| 青青草视频在线视频观看| 黄色视频不卡| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 最近2019中文字幕mv第一页| av在线观看视频网站免费| 婷婷色av中文字幕| 人人妻人人澡人人看| 国产精品麻豆人妻色哟哟久久| 亚洲欧洲日产国产| 亚洲,欧美精品.| 啦啦啦视频在线资源免费观看| 亚洲人成77777在线视频| av不卡在线播放| 亚洲一区二区三区欧美精品| 丝袜在线中文字幕| 极品人妻少妇av视频| 无遮挡黄片免费观看| 夫妻午夜视频| 欧美日本中文国产一区发布| 亚洲欧美一区二区三区久久| 久久99一区二区三区| 国产乱来视频区| 看非洲黑人一级黄片| 无遮挡黄片免费观看| 丰满乱子伦码专区| 深夜精品福利| 亚洲欧美一区二区三区久久| 天天操日日干夜夜撸| 日韩中文字幕欧美一区二区 | 国产精品国产三级专区第一集| 黄色 视频免费看| 日本vs欧美在线观看视频| 在线观看免费视频网站a站| 两个人看的免费小视频| 国产精品一区二区在线观看99| 男人爽女人下面视频在线观看| 超色免费av| 伊人久久大香线蕉亚洲五| 亚洲综合色网址| 国产成人一区二区在线| 国产成人a∨麻豆精品|