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

    基于析因設(shè)計(jì)法的隧道圍巖穩(wěn)定性評(píng)價(jià)

    2017-07-06 09:59:38
    關(guān)鍵詞:響應(yīng)值圍巖概率

    高 山

    (中鐵十九局集團(tuán)第七工程有限公司,廣東 珠海 519029)

    ?

    基于析因設(shè)計(jì)法的隧道圍巖穩(wěn)定性評(píng)價(jià)

    高 山

    (中鐵十九局集團(tuán)第七工程有限公司,廣東 珠海 519029)

    由于隧道的賦存環(huán)境復(fù)雜,其圍巖力學(xué)參數(shù)充滿了不確定性,采用單一定值計(jì)算不太嚴(yán)謹(jǐn)。為了給隧道圍巖穩(wěn)定設(shè)計(jì)提供更嚴(yán)謹(jǐn)?shù)囊罁?jù),以隧道圍巖的粘聚力c、內(nèi)摩擦角φ、彈性模量E和側(cè)壓力系數(shù)K4個(gè)因素為研究指標(biāo),根據(jù)析因設(shè)計(jì)表,每個(gè)指標(biāo)以其均值和標(biāo)準(zhǔn)差來進(jìn)行計(jì)算,求得對(duì)應(yīng)響應(yīng)值;根據(jù)線性回歸得到回歸方程后,結(jié)合極限狀態(tài)函數(shù)可獲得圍巖穩(wěn)定性功能函數(shù),最后求得圍巖的失效概率。結(jié)合林嶺隧道,對(duì)實(shí)際工程進(jìn)行了研究,構(gòu)建了林嶺隧道圍巖穩(wěn)定的功能函數(shù),并采用Matlab算得工程隧道的失效概率。結(jié)果表明:析因設(shè)計(jì)法在4因素情況下需進(jìn)行的試驗(yàn)次數(shù)為16次,這為隧道可靠度研究提供了一種更高效而嚴(yán)謹(jǐn)?shù)姆椒ā?/p>

    隧道;圍巖穩(wěn)定;析因設(shè)計(jì);失效概率

    隧道結(jié)構(gòu)體系采用單一安全系數(shù)設(shè)計(jì)時(shí)往往不能衡量工程可靠度大小??煽慷确椒ㄗ鳛橐环N考慮參數(shù)變異或不確定性的方法,多以參數(shù)的均值和方差來計(jì)算結(jié)構(gòu)可靠性,因而用可靠度來表征隧道結(jié)構(gòu)系統(tǒng)的整體性能更為合理[1]。

    可靠度理論發(fā)展至今,形成了一系列方法。一階矩法(FOSM)、二階矩法(SORM)和蒙特卡洛法(Monte Carlo Simulation)最為常見[2-5]。隨著計(jì)算機(jī)的應(yīng)用發(fā)展,數(shù)值模擬方法通過生成樣本點(diǎn)模擬結(jié)構(gòu)參數(shù)的隨機(jī)性來估計(jì)失效概率,衍生了多種數(shù)值模擬算法,廣泛應(yīng)用于實(shí)際工程結(jié)構(gòu)的可靠性分析,如子集模擬法(Subset Simulation)、隨機(jī)有限元法(RFEM)。在隧道結(jié)構(gòu)可靠度計(jì)算時(shí)應(yīng)用響應(yīng)面法非常普遍。趙旭峰和嚴(yán)松宏[2]采用響應(yīng)面法結(jié)合有限元得到響應(yīng)值,對(duì)一座既有隧道進(jìn)行了可靠性分析,認(rèn)為響應(yīng)面法在隧道結(jié)構(gòu)穩(wěn)定分析方面具有良好前景;宋玉香等[4]考慮影響隧道的具有不確定性特征的多種參數(shù), 試驗(yàn)樣本點(diǎn)采用中心復(fù)合設(shè)計(jì)法選取,建立了隧道在地震荷載作用下的響應(yīng)面模型,提供了隧道的震害評(píng)估手段,具有實(shí)際工程應(yīng)用價(jià)值;梁斌[6]等采用析因設(shè)計(jì)法對(duì)隧道有限元模型進(jìn)行了修正,包括步驟樣本設(shè)計(jì)、因素篩選、模型選擇以及數(shù)據(jù)檢驗(yàn),并與實(shí)際工程應(yīng)用進(jìn)行對(duì)比,發(fā)現(xiàn)修正后有限元模型的計(jì)算結(jié)果與實(shí)測(cè)結(jié)果接近。肖志鵬等[7]學(xué)者用響應(yīng)面方法對(duì)隧道靜動(dòng)力有限元模型進(jìn)行了修正,并將其用于解決實(shí)際工程問題,克服了其他有限元模型修正法存在局部解及精度不夠等問題。總結(jié)發(fā)現(xiàn),響應(yīng)面法的樣本點(diǎn)布置較多,而且存在局部解甚至失真等問題,找尋更好的方法并獲取全局解很有必要。

    析因設(shè)計(jì)法是近幾年發(fā)展起來的進(jìn)行結(jié)構(gòu)可靠性分析的一種有效的方法,最早是由數(shù)學(xué)家Box和Wilson于1951年提出來的[8],其原始意思是用一個(gè)合適的修勻函數(shù),即響應(yīng)面,近似表達(dá)一個(gè)未知的函數(shù)。當(dāng)系統(tǒng)的參數(shù)和系統(tǒng)的輸出響應(yīng)之間的關(guān)系以某種隱含的方式存在時(shí),析因設(shè)計(jì)法無疑提供了一種近似表達(dá)這種隱含關(guān)系的合適手段。從本質(zhì)上講,析因設(shè)計(jì)法是一項(xiàng)統(tǒng)計(jì)學(xué)綜合實(shí)驗(yàn)技術(shù),用于處理幾個(gè)變量對(duì)一個(gè)體系或結(jié)構(gòu)的作用問題,也就是結(jié)構(gòu)或體系的輸入(變量)與輸出(響應(yīng))的轉(zhuǎn)換關(guān)系問題。而均勻設(shè)計(jì)是只考慮試驗(yàn)點(diǎn)在試驗(yàn)范圍內(nèi)均勻散布的一種試驗(yàn)設(shè)計(jì)方法,較析因設(shè)計(jì)法,可減少試驗(yàn)次數(shù)。本文將析因設(shè)計(jì)法作為一種非線性系統(tǒng)的可靠性分析方法,結(jié)合均勻設(shè)計(jì)原理,以數(shù)值軟件模擬為試驗(yàn)手段,為隧道工程結(jié)構(gòu)可靠度分析提供了一種新的思路。

    1 析因設(shè)計(jì)法簡介

    析因設(shè)計(jì)法出現(xiàn)于十九世紀(jì),該方法應(yīng)用于數(shù)據(jù)統(tǒng)計(jì)和實(shí)驗(yàn)設(shè)計(jì)中[9],其基本思想是根據(jù)數(shù)學(xué)里的組合原理,將樣本點(diǎn)根據(jù)設(shè)計(jì)表布置,然后對(duì)應(yīng)求得每組數(shù)據(jù)的響應(yīng)值,并采用回歸分析得到回歸方程。

    析因設(shè)計(jì)法的使用包括以下2個(gè)主要內(nèi)容:①參數(shù)設(shè)計(jì);②系數(shù)估計(jì)。

    參數(shù)設(shè)計(jì)即試驗(yàn)樣本點(diǎn)的布置如圖1所示,圖1為一個(gè)包含2個(gè)變量的析因設(shè)計(jì)原理圖。x1和x2包括水平+1和-1,根據(jù)組合原理可知,存在22個(gè)組合,即(+1,+1)、(+1,-1)、(-1,+1)和(-1,-1)。以此類推,包含n個(gè)因素的情況下,即包括2n組試驗(yàn)設(shè)計(jì)。

    參數(shù)設(shè)計(jì)的任務(wù)是在參數(shù)空間中選擇合理的采樣點(diǎn),通過采樣點(diǎn)的輸入輸出數(shù)據(jù)構(gòu)造響應(yīng)面,使其在感興趣的有限區(qū)域內(nèi)能夠有效地逼近真實(shí)響應(yīng)。4因子復(fù)合設(shè)計(jì)見表1,可知其試驗(yàn)總數(shù)為16。參數(shù)設(shè)計(jì)完之后,根據(jù)設(shè)計(jì)表即可對(duì)應(yīng)計(jì)算或做實(shí)驗(yàn)獲取最后的響應(yīng)值;然后求得每個(gè)參數(shù)的回歸系數(shù),一般采用ANOVA,即方差分析法;最后構(gòu)建多項(xiàng)式函數(shù),函數(shù)由一系列因變量x1,x2,…,xn與響應(yīng)值y組成, 如公式(1):

    (1)

    式(1)中,a,bj,bjj分別為常數(shù)項(xiàng)、一次項(xiàng)和二次項(xiàng)的回歸系數(shù)。選擇響應(yīng)面表達(dá)形式時(shí),一方面要盡可能簡單,另一方面要考慮到能夠足夠靈活地反映各種不同的真實(shí)曲面形狀(因?yàn)榍娴恼鎸?shí)形狀并不知道)。如果采用二次超曲面作為響應(yīng)面函數(shù),則對(duì)任意有限維問題,響應(yīng)面只存在著某種意義上的逼近解。

    表1 4因子復(fù)合設(shè)計(jì)

    2 可靠度指標(biāo)求解

    假設(shè)某系統(tǒng)存在抗力R,外加荷載為S,系統(tǒng)功能函數(shù)為G(x)=R-S,即極限狀態(tài)函數(shù)。

    Pf=P[G(x)<0]≈Φ(β)

    (2)

    式(2)中,Pf為失效概率;β為可靠度指標(biāo)。極限狀態(tài)函數(shù)及可靠度指標(biāo)示意圖如圖2所示,可靠度指標(biāo)和失效概率存在一一對(duì)應(yīng)的關(guān)系,可靠度指標(biāo)需滿足式(3)[10]:

    (3)

    式(3)中,Y為響應(yīng)值;C為變量相關(guān)系數(shù);xmin,xmax分別為變量的最小值和最大值。公式(3)求解可以借助Matlab的優(yōu)化工具箱進(jìn)行。可靠度指標(biāo)確定后,即可求得系統(tǒng)的失效概率。

    失效概率求解基本思路如下:①根據(jù)析因設(shè)計(jì)表進(jìn)行樣本點(diǎn)布置,編碼后的數(shù)值需進(jìn)行轉(zhuǎn)碼;②根據(jù)ANOVA進(jìn)行回歸系數(shù)分析,不顯著變量可以適當(dāng)剔除,最后確定回歸方程;③確定極限狀態(tài)方程G(x)=0;④借助Matlab的優(yōu)化工具箱,獲得可靠度指標(biāo)及失效概率。

    3 工程案例分析

    3.1 工程概況

    林嶺隧道是海南瓊中至樂東高速公路的雙修隧道,位于海南省樂東縣樂光農(nóng)場南西側(cè)500 m處。林嶺隧道為小凈距隧道,左線起止樁號(hào)(ZK223+220)~(ZK223+500),長280 m;右線起止樁號(hào)(K223+225)~(K223+510),長285 m。該隧道屬構(gòu)造剝蝕低丘地貌單元,山坡較緩,植被較發(fā)育。隧址區(qū)高程196.0~250.5 m,相對(duì)高差54.5 m。隧道最大埋深38.8 m;隧道穿越山坡地段,山坡坡度一般在15~30°,局部達(dá)45°,左右洞室凈間距9~19 m。隧道進(jìn)出口均采用削竹式洞門,其明洞均為20 m。

    3.2 可靠度分析

    林嶺隧道圍巖力學(xué)參數(shù)彈性模量、粘聚力、內(nèi)摩擦角和側(cè)壓力系數(shù)見表2。4個(gè)參數(shù)的變異系數(shù)取0.2,也就是標(biāo)準(zhǔn)差為均值的0.2倍。確定了4個(gè)參數(shù)的均值和標(biāo)準(zhǔn)差后,就可以按照表1來轉(zhuǎn)換編碼。根據(jù)蘇永華[10]等的研究結(jié)果,一般取x=μ+1.645σ和x=μ- 1.645σ對(duì)應(yīng)的編碼+1和-1。

    表2 隧道力學(xué)參數(shù)

    林嶺隧道建模如圖3所示,取拱頂下沉為響應(yīng)值,采用ABAQUS軟件對(duì)隧道結(jié)構(gòu)的應(yīng)變進(jìn)行了分析。林嶺隧道析因設(shè)計(jì)見表3。根據(jù)表3所示的每個(gè)不同組合,分別建模計(jì)算每個(gè)組對(duì)應(yīng)的位移值,共16組。

    根據(jù)表3及公式(1),可以獲得拱頂下沉值的回歸函數(shù):

    g(x)= 0.0904-0.0058E-0.4788c-0.0044φ+0.292K+0.0046Ec+0.0002φE-0.002EK+0.0095cφ-0.0754cK-0.0035φK

    (4)

    拱頂下沉容許位移u可根據(jù)鐵路隧道監(jiān)控量測(cè)技術(shù)規(guī)程[11]查到,其極限狀態(tài)函數(shù)(即功能函數(shù))為:

    G(x)=g(x)-u

    (5)

    采用Matlab工具箱,可算得圍巖失效概率為6.626%,對(duì)應(yīng)的可靠度指標(biāo)為1.5。

    試驗(yàn)點(diǎn)序編碼變量x1x2x3x4自然變量彈性模量E/GPa粘聚力c/GPa內(nèi)摩擦角φ/°側(cè)壓力系數(shù)K響應(yīng)值位移/mm1-1-1-1-121.220.2329.752.010.15722 1-1-1-142.020.2329.752.010.07993-1 1-1-121.220.4629.752.010.10734 1 1-1-142.020.4629.752.010.05685-1-1 1-121.220.2349.252.010.30966 1-1 1-142.020.2349.252.010.05457-1 1 1-121.220.4649.252.010.02898 1 1 1-142.020.4649.252.010.15649-1-1-1 121.220.2329.752.990.042410 1-1-1 142.020.2329.752.990.215611-1 1-1 121.220.4629.752.990.097712 1 1-1 142.020.4629.752.990.021513-1-1 1 121.220.2349.252.990.109914 1-1 1 142.020.2349.252.990.049715-1 1 1 121.220.4649.252.990.081216 1 1 1 142.020.4649.252.990.0412

    4 結(jié)論

    1) 析因設(shè)計(jì)法需要較多的抽樣次數(shù),抽樣次數(shù)為2n次,可以滿足計(jì)算精度要求,得到的結(jié)果不失真。

    2) 構(gòu)建了隧道圍巖穩(wěn)定的功能函數(shù),并采用Matlab算得工程隧道的失效概率,這為隧道可靠度研究提供了一種更高效而嚴(yán)謹(jǐn)?shù)姆椒ā?/p>

    3) 林嶺隧道的圍巖失效概率為6.626%,對(duì)應(yīng)的可靠度指標(biāo)為1.5。

    [1] 嚴(yán)春風(fēng),劉東燕,張建輝,等.巖土工程可靠度關(guān)于強(qiáng)度參數(shù)分布函數(shù)概型的敏感度分析[J].巖石力學(xué)與工程學(xué)報(bào),1999,18(1):36-39.

    [2] 趙旭峰,嚴(yán)松宏.響應(yīng)面法在隧道襯砌結(jié)構(gòu)可靠度分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2003,22(z2):2853-2856.

    [3] 李宗洋.深埋隧道襯砌結(jié)構(gòu)的可靠性方法研究[D].成都:西南石油大學(xué),2014.

    [4] 宋玉香,劉勇,朱永全.響應(yīng)面方法在整體式隧道襯砌可靠性分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2004,23(11):1847-1851.

    [5] 李奎,趙東平,路軍富.深埋隧道素混凝土襯砌可靠度計(jì)算模型研究[J].鐵道建筑,2014(7):38-42.

    [6] 梁斌.隧道圍巖與支護(hù)結(jié)構(gòu)穩(wěn)定可靠性分析方法研究[D].長沙:湖南大學(xué),2015.

    [7] 肖志鵬,呂慶,趙宇,等.基于移動(dòng)最小二乘法的隧道支護(hù)穩(wěn)定可靠度分析[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2016,49(5):683-689.

    [8] Box GEP.Statistics for experimenters design innovation and discovery for Box George E P[J].Technometrics,2006, 60(48):303-304.

    [9] 駱勇鵬.基于響應(yīng)法的橋梁結(jié)構(gòu)有限元模型靜動(dòng)力修正方法研究[D].蘭州:蘭州理工大學(xué),2013.

    [10] 蘇永華,李翔,丁云,等.基于二次正交試驗(yàn)優(yōu)化的隧道圍巖穩(wěn)定可靠度方法[J].巖土工程學(xué)報(bào),2012,34(2):326-332.

    [11] 中國鐵路總公司.鐵路隧道監(jiān)控量測(cè)技術(shù)規(guī)程Q/CR9218-2015[S].北京:中國鐵路出版社,2015.

    (責(zé)任編輯 吳鴻霞)

    Reliability Assessment of Surrounding Rock in Tunnel Based on Factorial Design

    GaoShan

    (No.7 Engineering Co.Ltd.,China Railway 19th Bureau Group,Zhuhai Guangdong 519029)

    The parameters in tunnel are full of uncertainties due to the complex geology,it is not rigorous to compute the stability of tunnel keep constant value.In order to provide a more rigorous approach for the stability of surrounding rock in tunnels,the mechanical parameters,i.e.cohesioncand internal friction angleφ,elastic modulusEand lateral pressure coefficientKare settled as variables,the response value is obtained with those four variables represented by their mean and standard deviation according to the factorial design table;regression equation can also be obtained according to the linear regression,the probability of failure of surrounding rock in tunnel can be calculated by the performance function by combining the limit state function with regression equation.Combined with the illustrated example,the stability of Linling Tunnel is calculated,the results show that the factorial design method should be tested for 16 times with 4 factors,the probability of failure of the tunnel is obtained with the help of Matlab,which provides a new method for tunnel designation.

    tunnel;stability of surrounding rock;factorial design;probability of failure

    2017-03-23

    國家自然科學(xué)基金項(xiàng)目(項(xiàng)目編號(hào):51378195)。

    高山,工程師,本科。

    10.3969/j.issn.2095-4565.2017.03.011

    U25

    A

    2095-4565(2017)03-0047-05

    猜你喜歡
    響應(yīng)值圍巖概率
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    概率與統(tǒng)計(jì)(二)
    概率與統(tǒng)計(jì)(一)
    基于熒光光譜技術(shù)的不同食用淀粉的快速區(qū)分
    隧道開挖圍巖穩(wěn)定性分析
    提高環(huán)境監(jiān)測(cè)數(shù)據(jù)準(zhǔn)確性初探
    紫外熒光法測(cè)硫各氣路流量對(duì)響應(yīng)值的影響
    山東化工(2019年1期)2019-01-24 03:00:16
    軟弱破碎圍巖隧道初期支護(hù)大變形治理技術(shù)
    江西建材(2018年4期)2018-04-10 12:37:22
    采空側(cè)巷道圍巖加固與巷道底臌的防治
    久久久久国内视频| 国产av在哪里看| 欧美+亚洲+日韩+国产| 国产三级在线视频| 国产av一区二区精品久久| 久久精品国产清高在天天线| 精品国产乱码久久久久久男人| 国产成人精品久久二区二区91| 久久久久久久久久久久大奶| 三级毛片av免费| 色尼玛亚洲综合影院| 可以在线观看毛片的网站| 国产精品久久久久成人av| 狂野欧美激情性xxxx| 丰满人妻熟妇乱又伦精品不卡| 丰满人妻熟妇乱又伦精品不卡| 男男h啪啪无遮挡| 色婷婷久久久亚洲欧美| 亚洲男人的天堂狠狠| 日韩欧美在线二视频| 看免费av毛片| 久久国产精品影院| 国产视频一区二区在线看| 日韩精品青青久久久久久| 欧美老熟妇乱子伦牲交| 交换朋友夫妻互换小说| 日韩一卡2卡3卡4卡2021年| 久久伊人香网站| 久久国产亚洲av麻豆专区| 在线观看免费视频日本深夜| 美女午夜性视频免费| 一二三四社区在线视频社区8| 免费观看人在逋| 成人手机av| 热99re8久久精品国产| 最新在线观看一区二区三区| 脱女人内裤的视频| 在线观看66精品国产| 亚洲精品在线观看二区| 亚洲av成人av| 视频在线观看一区二区三区| 亚洲国产欧美一区二区综合| 国产黄a三级三级三级人| 在线观看一区二区三区激情| 身体一侧抽搐| 97人妻天天添夜夜摸| 久久国产亚洲av麻豆专区| 精品第一国产精品| 精品日产1卡2卡| 丝袜美腿诱惑在线| 国产成人精品在线电影| 五月开心婷婷网| 国产三级黄色录像| 老司机亚洲免费影院| 欧美乱色亚洲激情| 国产精品免费视频内射| 亚洲欧美精品综合一区二区三区| 精品国产超薄肉色丝袜足j| av视频免费观看在线观看| 日韩av在线大香蕉| 国产精品 欧美亚洲| 色婷婷久久久亚洲欧美| 午夜福利一区二区在线看| 最近最新中文字幕大全免费视频| 成熟少妇高潮喷水视频| 欧美精品亚洲一区二区| 少妇被粗大的猛进出69影院| 十八禁人妻一区二区| 在线播放国产精品三级| 国产三级在线视频| 精品一区二区三卡| 99国产精品一区二区蜜桃av| 99久久综合精品五月天人人| 国产伦一二天堂av在线观看| 一进一出抽搐gif免费好疼 | 热99re8久久精品国产| 搡老乐熟女国产| 免费在线观看视频国产中文字幕亚洲| 悠悠久久av| 另类亚洲欧美激情| 精品福利永久在线观看| 中出人妻视频一区二区| 制服人妻中文乱码| 精品一区二区三区视频在线观看免费 | av天堂久久9| 中文字幕精品免费在线观看视频| 欧美日韩亚洲高清精品| 中文字幕人妻丝袜一区二区| 国产av精品麻豆| 波多野结衣一区麻豆| 法律面前人人平等表现在哪些方面| netflix在线观看网站| 身体一侧抽搐| 亚洲av成人av| 国产欧美日韩一区二区精品| 99国产精品99久久久久| 亚洲全国av大片| 黄片小视频在线播放| 大型av网站在线播放| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 一区福利在线观看| av网站免费在线观看视频| 久久久国产成人免费| xxxhd国产人妻xxx| 夜夜看夜夜爽夜夜摸 | 免费女性裸体啪啪无遮挡网站| 国产成人av教育| 9色porny在线观看| 黄网站色视频无遮挡免费观看| 男人舔女人下体高潮全视频| 少妇裸体淫交视频免费看高清 | 在线观看免费日韩欧美大片| 最近最新中文字幕大全电影3 | 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 久久精品国产综合久久久| 12—13女人毛片做爰片一| 80岁老熟妇乱子伦牲交| tocl精华| 最近最新中文字幕大全电影3 | 国产色视频综合| 亚洲成a人片在线一区二区| 色尼玛亚洲综合影院| 制服人妻中文乱码| 欧美日韩视频精品一区| 久久精品亚洲精品国产色婷小说| 韩国av一区二区三区四区| 日韩大尺度精品在线看网址 | 亚洲色图综合在线观看| 亚洲免费av在线视频| 又大又爽又粗| 欧美激情极品国产一区二区三区| 国产不卡一卡二| 丰满饥渴人妻一区二区三| 成人精品一区二区免费| 国产成人一区二区三区免费视频网站| 大型av网站在线播放| 精品电影一区二区在线| 免费观看精品视频网站| av福利片在线| 久久精品影院6| 老汉色av国产亚洲站长工具| 国产精品九九99| 国产av精品麻豆| 大型av网站在线播放| 成年人黄色毛片网站| 午夜老司机福利片| www.熟女人妻精品国产| xxx96com| 成人亚洲精品一区在线观看| 超碰97精品在线观看| 精品国内亚洲2022精品成人| 亚洲欧洲精品一区二区精品久久久| 熟女少妇亚洲综合色aaa.| 成人三级做爰电影| 日本wwww免费看| 国产97色在线日韩免费| 一级片'在线观看视频| 国产精品永久免费网站| 国产成人啪精品午夜网站| 久久久久久久久免费视频了| 亚洲欧美日韩高清在线视频| 欧美大码av| 久久 成人 亚洲| 老司机在亚洲福利影院| 天堂√8在线中文| 精品一区二区三区四区五区乱码| 一本大道久久a久久精品| 国产三级在线视频| 啪啪无遮挡十八禁网站| 久久精品国产综合久久久| 999精品在线视频| 免费高清视频大片| 欧美色视频一区免费| 人妻丰满熟妇av一区二区三区| 黄网站色视频无遮挡免费观看| 亚洲成人国产一区在线观看| 成人黄色视频免费在线看| 成人18禁高潮啪啪吃奶动态图| 久久婷婷成人综合色麻豆| 久久精品影院6| 午夜精品久久久久久毛片777| xxx96com| 亚洲伊人色综图| 岛国视频午夜一区免费看| 嫩草影院精品99| 黄色片一级片一级黄色片| 老司机福利观看| av国产精品久久久久影院| 纯流量卡能插随身wifi吗| 午夜福利在线免费观看网站| 大码成人一级视频| 欧美日韩av久久| 欧美午夜高清在线| 欧美成狂野欧美在线观看| 国产精品永久免费网站| 免费一级毛片在线播放高清视频 | 精品久久久久久成人av| 老司机午夜福利在线观看视频| 国产日韩一区二区三区精品不卡| 欧美日韩瑟瑟在线播放| 变态另类成人亚洲欧美熟女 | 黄色怎么调成土黄色| 欧美乱色亚洲激情| 亚洲伊人色综图| 国产精品九九99| 亚洲人成电影免费在线| xxxhd国产人妻xxx| 一边摸一边抽搐一进一出视频| 波多野结衣av一区二区av| 亚洲欧美日韩无卡精品| 神马国产精品三级电影在线观看 | 12—13女人毛片做爰片一| 久热爱精品视频在线9| 一二三四在线观看免费中文在| 在线观看www视频免费| 日韩成人在线观看一区二区三区| 国产蜜桃级精品一区二区三区| 国产免费现黄频在线看| 天堂√8在线中文| 久久久久久久久久久久大奶| 免费在线观看亚洲国产| 亚洲欧美日韩高清在线视频| 身体一侧抽搐| 亚洲七黄色美女视频| 成熟少妇高潮喷水视频| 女同久久另类99精品国产91| 色综合站精品国产| 色综合欧美亚洲国产小说| 成年版毛片免费区| 国产欧美日韩精品亚洲av| 国产欧美日韩综合在线一区二区| 在线观看免费视频日本深夜| 午夜福利在线免费观看网站| 成人亚洲精品av一区二区 | 亚洲一码二码三码区别大吗| xxx96com| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 咕卡用的链子| 色在线成人网| 午夜免费激情av| 在线观看舔阴道视频| 精品国产一区二区三区四区第35| 午夜免费激情av| 国产成人免费无遮挡视频| 国产又色又爽无遮挡免费看| 午夜精品久久久久久毛片777| 曰老女人黄片| 国产有黄有色有爽视频| 看片在线看免费视频| 亚洲成人免费av在线播放| 波多野结衣av一区二区av| 精品高清国产在线一区| 悠悠久久av| 最好的美女福利视频网| 男女床上黄色一级片免费看| 91成年电影在线观看| а√天堂www在线а√下载| 一边摸一边抽搐一进一小说| 大香蕉久久成人网| 欧美日韩福利视频一区二区| 国产精品美女特级片免费视频播放器 | 欧美在线一区亚洲| 一级片'在线观看视频| 国产视频一区二区在线看| 亚洲精品国产精品久久久不卡| 女人精品久久久久毛片| 日韩中文字幕欧美一区二区| 成在线人永久免费视频| 国产精品永久免费网站| 精品一区二区三卡| 亚洲成a人片在线一区二区| 亚洲精品在线美女| www.www免费av| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女 | 色综合站精品国产| 精品国产乱子伦一区二区三区| 国产不卡一卡二| 亚洲人成伊人成综合网2020| 亚洲专区中文字幕在线| 成人黄色视频免费在线看| 精品福利观看| 亚洲欧美一区二区三区黑人| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 午夜免费鲁丝| 岛国视频午夜一区免费看| 美女扒开内裤让男人捅视频| 一级a爱片免费观看的视频| 女警被强在线播放| 老司机靠b影院| 免费看十八禁软件| 手机成人av网站| 欧美黄色片欧美黄色片| 欧美成人免费av一区二区三区| 国产一区二区三区视频了| 国产伦人伦偷精品视频| 黄色毛片三级朝国网站| 日韩 欧美 亚洲 中文字幕| 久久这里只有精品19| 一边摸一边抽搐一进一出视频| 日本黄色日本黄色录像| 欧美另类亚洲清纯唯美| 成人国产一区最新在线观看| 日韩有码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产成年人精品一区二区 | 久久久久亚洲av毛片大全| 日韩三级视频一区二区三区| 日韩欧美在线二视频| 韩国精品一区二区三区| 男人的好看免费观看在线视频 | 一区福利在线观看| 在线av久久热| 亚洲国产欧美日韩在线播放| 久久久久久久精品吃奶| 久久国产亚洲av麻豆专区| 国产精品二区激情视频| 大码成人一级视频| 男人舔女人的私密视频| 黄片大片在线免费观看| 国产免费男女视频| 99国产精品一区二区三区| 亚洲成av片中文字幕在线观看| 精品久久蜜臀av无| 自线自在国产av| 黄片大片在线免费观看| 成年版毛片免费区| 黑丝袜美女国产一区| 国产精品电影一区二区三区| 亚洲精品一二三| 国产区一区二久久| 亚洲三区欧美一区| 18禁裸乳无遮挡免费网站照片 | 久久中文字幕一级| 久久伊人香网站| 三上悠亚av全集在线观看| 制服诱惑二区| 欧美日韩精品网址| 欧美日韩亚洲综合一区二区三区_| 极品人妻少妇av视频| 在线观看www视频免费| 亚洲精品中文字幕一二三四区| 欧美最黄视频在线播放免费 | 欧美乱妇无乱码| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀| 中文字幕色久视频| 亚洲美女黄片视频| 美女 人体艺术 gogo| 亚洲专区国产一区二区| 午夜亚洲福利在线播放| 91大片在线观看| 好男人电影高清在线观看| 大型av网站在线播放| 久久国产亚洲av麻豆专区| 精品国产美女av久久久久小说| 桃色一区二区三区在线观看| 欧美日韩黄片免| 18禁美女被吸乳视频| 欧美成狂野欧美在线观看| 亚洲国产欧美网| 99精品欧美一区二区三区四区| 欧美另类亚洲清纯唯美| 国产av一区二区精品久久| 美女 人体艺术 gogo| 97碰自拍视频| 欧美色视频一区免费| 久久国产亚洲av麻豆专区| 国产一区二区在线av高清观看| 国产激情久久老熟女| 人人妻,人人澡人人爽秒播| 真人一进一出gif抽搐免费| 欧美激情久久久久久爽电影 | 51午夜福利影视在线观看| 日韩大尺度精品在线看网址 | 五月开心婷婷网| 操出白浆在线播放| 黄色 视频免费看| 久久精品91无色码中文字幕| www日本在线高清视频| 久久天堂一区二区三区四区| 国产成人啪精品午夜网站| 日韩成人在线观看一区二区三区| 很黄的视频免费| 国产成人欧美| 又大又爽又粗| 亚洲va日本ⅴa欧美va伊人久久| 男女之事视频高清在线观看| 一本综合久久免费| 亚洲中文av在线| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 日韩av在线大香蕉| 老司机午夜福利在线观看视频| 一本综合久久免费| 乱人伦中国视频| 极品教师在线免费播放| 中亚洲国语对白在线视频| 成人影院久久| 久久久久久人人人人人| 19禁男女啪啪无遮挡网站| 成人国语在线视频| 水蜜桃什么品种好| 国产xxxxx性猛交| 大香蕉久久成人网| 在线观看免费高清a一片| 一级a爱片免费观看的视频| 亚洲精品国产区一区二| 夫妻午夜视频| 亚洲少妇的诱惑av| 中文字幕人妻熟女乱码| 一级a爱视频在线免费观看| 久久人妻熟女aⅴ| 91在线观看av| 国产精品日韩av在线免费观看 | 性少妇av在线| 亚洲一区二区三区色噜噜 | 欧美黄色淫秽网站| 欧美+亚洲+日韩+国产| 99久久精品国产亚洲精品| 日韩免费av在线播放| 日韩高清综合在线| 一二三四在线观看免费中文在| 男女之事视频高清在线观看| 欧美乱色亚洲激情| av在线天堂中文字幕 | 国产视频一区二区在线看| 日韩av在线大香蕉| 国产精品国产高清国产av| 精品免费久久久久久久清纯| 高清欧美精品videossex| 国产精品一区二区三区四区久久 | 91字幕亚洲| 国产av又大| 国产精品亚洲一级av第二区| 久久久久国内视频| 久久香蕉精品热| 亚洲九九香蕉| 日韩国内少妇激情av| 亚洲性夜色夜夜综合| 看片在线看免费视频| 日韩三级视频一区二区三区| 亚洲 欧美一区二区三区| 91麻豆精品激情在线观看国产 | 国产精品免费一区二区三区在线| 法律面前人人平等表现在哪些方面| 久久香蕉激情| 少妇 在线观看| tocl精华| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| av视频免费观看在线观看| 亚洲在线自拍视频| 久久久精品欧美日韩精品| 免费日韩欧美在线观看| 黄片大片在线免费观看| 97碰自拍视频| 高清av免费在线| 欧美午夜高清在线| 久久国产亚洲av麻豆专区| 少妇裸体淫交视频免费看高清 | 嫩草影院精品99| 99国产精品一区二区三区| 一级,二级,三级黄色视频| 国产99白浆流出| 国产麻豆69| 午夜福利欧美成人| 国产日韩一区二区三区精品不卡| 久久久国产精品麻豆| 亚洲va日本ⅴa欧美va伊人久久| cao死你这个sao货| 亚洲一区二区三区色噜噜 | 国产精品久久久久成人av| 久久精品91无色码中文字幕| 纯流量卡能插随身wifi吗| av片东京热男人的天堂| 男女做爰动态图高潮gif福利片 | 日韩国内少妇激情av| 久久精品91无色码中文字幕| 精品电影一区二区在线| 老司机午夜福利在线观看视频| 欧美午夜高清在线| 看片在线看免费视频| 黄频高清免费视频| 国产精品99久久99久久久不卡| 青草久久国产| 天天躁夜夜躁狠狠躁躁| 久久久国产一区二区| 99热只有精品国产| 精品卡一卡二卡四卡免费| 性色av乱码一区二区三区2| 水蜜桃什么品种好| 久久久久久久精品吃奶| 国产精品久久久久成人av| 99久久人妻综合| 日韩三级视频一区二区三区| 制服诱惑二区| 亚洲全国av大片| 香蕉国产在线看| 9191精品国产免费久久| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 精品卡一卡二卡四卡免费| 亚洲中文av在线| 久久国产精品人妻蜜桃| 免费高清视频大片| xxxhd国产人妻xxx| 99精品欧美一区二区三区四区| 一级毛片精品| 黑丝袜美女国产一区| 成人永久免费在线观看视频| 欧美乱妇无乱码| 丰满饥渴人妻一区二区三| 青草久久国产| 国产欧美日韩一区二区三区在线| 亚洲国产欧美一区二区综合| 女同久久另类99精品国产91| 久久精品亚洲av国产电影网| 亚洲国产欧美日韩在线播放| 波多野结衣一区麻豆| 男人舔女人的私密视频| 国产精品亚洲av一区麻豆| 国产成人系列免费观看| 夫妻午夜视频| 国产精品二区激情视频| 成人精品一区二区免费| 中文字幕人妻丝袜制服| 99国产精品一区二区三区| 国产成人精品无人区| 又紧又爽又黄一区二区| 亚洲人成电影观看| 欧美亚洲日本最大视频资源| 老司机午夜十八禁免费视频| 村上凉子中文字幕在线| 狂野欧美激情性xxxx| 国产xxxxx性猛交| 成人18禁在线播放| 级片在线观看| 一a级毛片在线观看| 免费av中文字幕在线| 免费在线观看完整版高清| 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 9色porny在线观看| 天堂动漫精品| 日日夜夜操网爽| 女性生殖器流出的白浆| 亚洲在线自拍视频| av有码第一页| 欧美国产精品va在线观看不卡| av超薄肉色丝袜交足视频| 国产精品国产av在线观看| 女人被狂操c到高潮| 在线观看一区二区三区| 搡老乐熟女国产| 亚洲精品国产精品久久久不卡| netflix在线观看网站| 亚洲aⅴ乱码一区二区在线播放 | 最好的美女福利视频网| 欧美丝袜亚洲另类 | 日韩高清综合在线| 国产精品免费视频内射| 久久久久九九精品影院| 悠悠久久av| 老司机午夜十八禁免费视频| 18禁裸乳无遮挡免费网站照片 | 另类亚洲欧美激情| 黄色怎么调成土黄色| 欧美日韩中文字幕国产精品一区二区三区 | 香蕉丝袜av| 亚洲午夜理论影院| 叶爱在线成人免费视频播放| 手机成人av网站| 国产色视频综合| 波多野结衣高清无吗| 两人在一起打扑克的视频| 久久国产精品男人的天堂亚洲| 欧美大码av| 国产成人精品在线电影| 久久国产亚洲av麻豆专区| 免费女性裸体啪啪无遮挡网站| 国产精品乱码一区二三区的特点 | 狠狠狠狠99中文字幕| 国产精品免费视频内射| 如日韩欧美国产精品一区二区三区| 久久人人97超碰香蕉20202| 精品一区二区三区四区五区乱码| 波多野结衣av一区二区av| 美女扒开内裤让男人捅视频| 亚洲激情在线av| 国产片内射在线| 免费av毛片视频| xxxhd国产人妻xxx| 国产精品国产高清国产av| 视频区欧美日本亚洲| 真人做人爱边吃奶动态| 99精品欧美一区二区三区四区| 9色porny在线观看| 亚洲第一青青草原| 精品一区二区三区av网在线观看| 久久中文看片网| 夜夜躁狠狠躁天天躁| 精品久久久久久久毛片微露脸| 久久久久国产精品人妻aⅴ院| 最近最新中文字幕大全免费视频| 欧美成人午夜精品| 在线观看日韩欧美| 一级a爱片免费观看的视频| 久久影院123| 午夜福利免费观看在线|