• <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è)巷道圍巖加固與巷道底臌的防治
    春色校园在线视频观看| 又黄又爽又刺激的免费视频.| 久久久欧美国产精品| www.av在线官网国产| 天天躁夜夜躁狠狠久久av| 黑人高潮一二区| 99久久中文字幕三级久久日本| 欧美亚洲 丝袜 人妻 在线| 青春草视频在线免费观看| 日韩欧美一区视频在线观看 | 国产一区二区在线观看日韩| 最新中文字幕久久久久| 成人美女网站在线观看视频| 成人综合一区亚洲| 日韩,欧美,国产一区二区三区| 一级毛片黄色毛片免费观看视频| 亚洲av欧美aⅴ国产| av免费观看日本| 夫妻午夜视频| 国产黄片视频在线免费观看| 成人无遮挡网站| 亚洲人与动物交配视频| 中文字幕制服av| 亚洲人成网站高清观看| 亚洲精华国产精华液的使用体验| 岛国毛片在线播放| 国产免费福利视频在线观看| 亚洲精品久久久久久婷婷小说| 久久精品国产自在天天线| 18禁裸乳无遮挡动漫免费视频 | 日日啪夜夜爽| a级一级毛片免费在线观看| 十八禁网站网址无遮挡 | 日本黄色片子视频| 亚洲图色成人| 麻豆成人av视频| 全区人妻精品视频| 国产真实伦视频高清在线观看| 亚洲怡红院男人天堂| 久久精品国产自在天天线| 成人国产麻豆网| 精品人妻视频免费看| eeuss影院久久| 久久99精品国语久久久| 亚洲精品视频女| 免费观看在线日韩| 亚洲一区二区三区欧美精品 | 成年免费大片在线观看| 韩国av在线不卡| 99热全是精品| 国产淫片久久久久久久久| 国产一区二区三区av在线| 在线看a的网站| 欧美精品人与动牲交sv欧美| 波野结衣二区三区在线| 五月开心婷婷网| 国产成人精品福利久久| 天天一区二区日本电影三级| 尾随美女入室| 18+在线观看网站| 亚洲无线观看免费| 亚洲国产高清在线一区二区三| 亚洲精品,欧美精品| 国产高清国产精品国产三级 | 婷婷色综合www| 我要看日韩黄色一级片| 亚洲欧美日韩无卡精品| 欧美日韩在线观看h| 中国三级夫妇交换| 日本av手机在线免费观看| 国产美女午夜福利| 中文字幕人妻熟人妻熟丝袜美| 久久久a久久爽久久v久久| 亚洲久久久久久中文字幕| 亚洲精品国产av蜜桃| 国产欧美日韩一区二区三区在线 | 男女那种视频在线观看| 免费人成在线观看视频色| 精品99又大又爽又粗少妇毛片| 直男gayav资源| 亚洲av中文字字幕乱码综合| 国产毛片在线视频| 欧美日韩精品成人综合77777| 亚洲,一卡二卡三卡| 亚洲久久久久久中文字幕| 亚洲内射少妇av| 亚洲欧美中文字幕日韩二区| 在线亚洲精品国产二区图片欧美 | videos熟女内射| 精品少妇黑人巨大在线播放| 亚洲欧美一区二区三区黑人 | 热re99久久精品国产66热6| 国产色爽女视频免费观看| 欧美+日韩+精品| 91久久精品电影网| 精品国产三级普通话版| 免费电影在线观看免费观看| 亚州av有码| 搡老乐熟女国产| 久久韩国三级中文字幕| 亚洲国产精品成人久久小说| 免费少妇av软件| 成人亚洲欧美一区二区av| 午夜福利网站1000一区二区三区| 欧美亚洲 丝袜 人妻 在线| 日本午夜av视频| 国产精品人妻久久久影院| 亚洲欧洲日产国产| 精品久久久精品久久久| 2021少妇久久久久久久久久久| 国产v大片淫在线免费观看| 色综合色国产| 我的女老师完整版在线观看| 99热全是精品| 国产成人精品福利久久| 日本猛色少妇xxxxx猛交久久| 狂野欧美激情性bbbbbb| 国产探花极品一区二区| 亚洲美女视频黄频| 舔av片在线| 一本一本综合久久| 99热这里只有是精品50| 亚洲国产日韩一区二区| 男人舔奶头视频| 一级片'在线观看视频| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 肉色欧美久久久久久久蜜桃 | 91精品国产九色| 免费黄频网站在线观看国产| 亚洲在久久综合| 日韩一本色道免费dvd| xxx大片免费视频| 国产精品久久久久久av不卡| 人妻少妇偷人精品九色| 蜜桃亚洲精品一区二区三区| 国产色爽女视频免费观看| 日日摸夜夜添夜夜添av毛片| 伦理电影大哥的女人| 你懂的网址亚洲精品在线观看| 午夜日本视频在线| 草草在线视频免费看| 精品久久久久久久久av| 日韩成人av中文字幕在线观看| 最近手机中文字幕大全| 在线亚洲精品国产二区图片欧美 | 色哟哟·www| 久久人人爽av亚洲精品天堂 | 伦精品一区二区三区| 欧美变态另类bdsm刘玥| 国产成人freesex在线| 麻豆久久精品国产亚洲av| 久久久久精品久久久久真实原创| 成人国产av品久久久| 看十八女毛片水多多多| 亚洲av日韩在线播放| 日本午夜av视频| 久久影院123| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 久久人人爽av亚洲精品天堂 | 黄色日韩在线| 久久久精品94久久精品| 99久久人妻综合| 亚洲精品aⅴ在线观看| 美女xxoo啪啪120秒动态图| 亚洲久久久久久中文字幕| 天堂中文最新版在线下载 | 亚洲三级黄色毛片| eeuss影院久久| 久久国内精品自在自线图片| 看免费成人av毛片| 国产黄a三级三级三级人| 成年免费大片在线观看| 三级国产精品欧美在线观看| 综合色丁香网| 夜夜爽夜夜爽视频| 久久久久精品久久久久真实原创| 国产淫片久久久久久久久| 在线精品无人区一区二区三 | 99久久精品国产国产毛片| 国产午夜精品一二区理论片| 中国三级夫妇交换| 国产精品蜜桃在线观看| 日韩人妻高清精品专区| 一区二区三区免费毛片| 在线观看一区二区三区激情| 亚洲精品视频女| 一个人观看的视频www高清免费观看| tube8黄色片| 免费看a级黄色片| av国产久精品久网站免费入址| 日韩欧美一区视频在线观看 | 久久99蜜桃精品久久| 欧美日韩国产mv在线观看视频 | 男的添女的下面高潮视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久国产网址| 五月开心婷婷网| av播播在线观看一区| 我的女老师完整版在线观看| 在线观看美女被高潮喷水网站| 热99国产精品久久久久久7| 如何舔出高潮| 亚洲精品一二三| 1000部很黄的大片| 国内揄拍国产精品人妻在线| 男女边吃奶边做爰视频| 一区二区三区免费毛片| 如何舔出高潮| 色5月婷婷丁香| 亚洲国产精品专区欧美| 美女视频免费永久观看网站| 亚洲精品乱码久久久久久按摩| 亚洲精品影视一区二区三区av| 波多野结衣巨乳人妻| 一区二区av电影网| 蜜桃亚洲精品一区二区三区| 在线看a的网站| 啦啦啦中文免费视频观看日本| 国产高清三级在线| 亚洲美女视频黄频| 嘟嘟电影网在线观看| 日本-黄色视频高清免费观看| 在线观看美女被高潮喷水网站| 国产成人a区在线观看| 日韩一本色道免费dvd| 久久久精品94久久精品| 毛片女人毛片| 日本一本二区三区精品| 欧美日韩综合久久久久久| 欧美精品人与动牲交sv欧美| 日韩制服骚丝袜av| 欧美性猛交╳xxx乱大交人| 亚洲欧美中文字幕日韩二区| 欧美日韩亚洲高清精品| 国产午夜精品一二区理论片| 肉色欧美久久久久久久蜜桃 | 亚洲色图av天堂| 国产 一区精品| 好男人在线观看高清免费视频| www.av在线官网国产| 精品久久国产蜜桃| 在线观看一区二区三区激情| 亚洲国产精品专区欧美| 中文字幕免费在线视频6| 男女啪啪激烈高潮av片| 中文乱码字字幕精品一区二区三区| 日韩国内少妇激情av| 久久人人爽人人爽人人片va| 观看免费一级毛片| 精品久久久久久久久亚洲| 国产精品久久久久久精品电影小说 | 秋霞在线观看毛片| 日本-黄色视频高清免费观看| 国精品久久久久久国模美| 波多野结衣巨乳人妻| 国产精品一及| 熟妇人妻不卡中文字幕| 亚洲真实伦在线观看| 99久久精品国产国产毛片| 欧美激情久久久久久爽电影| 一级毛片电影观看| av专区在线播放| 国产成人91sexporn| 啦啦啦啦在线视频资源| 午夜福利高清视频| 好男人在线观看高清免费视频| 三级经典国产精品| 69人妻影院| 国产欧美日韩精品一区二区| 亚洲成人中文字幕在线播放| 嫩草影院入口| 久久久精品欧美日韩精品| 欧美性猛交╳xxx乱大交人| 男男h啪啪无遮挡| 午夜精品国产一区二区电影 | 亚洲第一区二区三区不卡| 久久99热这里只有精品18| 天天躁日日操中文字幕| 国产免费一区二区三区四区乱码| 亚洲欧洲国产日韩| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 搞女人的毛片| 男女下面进入的视频免费午夜| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 日韩伦理黄色片| 亚洲综合精品二区| 一级毛片aaaaaa免费看小| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 国产精品蜜桃在线观看| 一级毛片黄色毛片免费观看视频| 国产国拍精品亚洲av在线观看| 久久精品久久久久久久性| av福利片在线观看| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站 | 韩国av在线不卡| 日韩国内少妇激情av| 中国三级夫妇交换| 国产乱来视频区| 色网站视频免费| 国产视频首页在线观看| 精品99又大又爽又粗少妇毛片| 人妻夜夜爽99麻豆av| 观看美女的网站| 老师上课跳d突然被开到最大视频| 欧美精品人与动牲交sv欧美| 亚洲内射少妇av| 女人被狂操c到高潮| 99九九线精品视频在线观看视频| 久久久久久久大尺度免费视频| 制服丝袜香蕉在线| 1000部很黄的大片| 欧美日韩视频高清一区二区三区二| 欧美成人精品欧美一级黄| 九草在线视频观看| 精华霜和精华液先用哪个| 中文乱码字字幕精品一区二区三区| 亚洲国产精品999| 精品人妻视频免费看| 婷婷色综合www| 97超碰精品成人国产| 午夜激情久久久久久久| 久久久久久久大尺度免费视频| 亚洲av日韩在线播放| 日韩三级伦理在线观看| 在线看a的网站| 爱豆传媒免费全集在线观看| 色婷婷久久久亚洲欧美| 久久精品熟女亚洲av麻豆精品| 日本黄大片高清| 各种免费的搞黄视频| 亚洲av免费高清在线观看| 男人添女人高潮全过程视频| 久久99蜜桃精品久久| 午夜爱爱视频在线播放| 中文字幕av成人在线电影| 在现免费观看毛片| 欧美xxⅹ黑人| 成年人午夜在线观看视频| 国产熟女欧美一区二区| 日本欧美国产在线视频| 国产高潮美女av| 国产精品嫩草影院av在线观看| 亚洲电影在线观看av| 成人无遮挡网站| av天堂中文字幕网| 18禁裸乳无遮挡免费网站照片| h日本视频在线播放| 少妇猛男粗大的猛烈进出视频 | 嫩草影院精品99| 国产在线男女| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 久久久国产一区二区| 欧美97在线视频| 国产成人精品婷婷| 亚洲av中文字字幕乱码综合| 激情 狠狠 欧美| 国产成人午夜福利电影在线观看| 国产高清国产精品国产三级 | 免费看a级黄色片| 少妇熟女欧美另类| 日韩欧美精品v在线| 超碰97精品在线观看| 国产一区二区三区av在线| 成人鲁丝片一二三区免费| 国产精品偷伦视频观看了| 青春草国产在线视频| 国产极品天堂在线| 国产一级毛片在线| 在线观看一区二区三区| 国产久久久一区二区三区| 国产精品秋霞免费鲁丝片| 久久久久久久国产电影| 亚洲av男天堂| 99热这里只有是精品50| 国产精品福利在线免费观看| av免费在线看不卡| 国产精品久久久久久精品电影| 久久久久精品性色| 91精品一卡2卡3卡4卡| 国产精品久久久久久精品电影| 蜜桃亚洲精品一区二区三区| 亚洲激情五月婷婷啪啪| 2021少妇久久久久久久久久久| 少妇高潮的动态图| 九色成人免费人妻av| 久久精品国产亚洲av天美| 国产真实伦视频高清在线观看| 丰满少妇做爰视频| 国产成人精品一,二区| 2022亚洲国产成人精品| 美女高潮的动态| 国产真实伦视频高清在线观看| 香蕉精品网在线| 久久精品国产自在天天线| 高清毛片免费看| 久久久成人免费电影| 观看免费一级毛片| 国产亚洲午夜精品一区二区久久 | 国产一区亚洲一区在线观看| 80岁老熟妇乱子伦牲交| 午夜免费鲁丝| 老司机影院毛片| 日本一二三区视频观看| 日本黄大片高清| 亚洲成人久久爱视频| 22中文网久久字幕| 超碰97精品在线观看| 精品久久久久久久久亚洲| 国产久久久一区二区三区| 国产精品一及| 日韩欧美一区视频在线观看 | 大片电影免费在线观看免费| 亚洲精品一二三| 亚洲精品久久午夜乱码| 夫妻性生交免费视频一级片| 日本欧美国产在线视频| 亚洲一区二区三区欧美精品 | 少妇丰满av| 国产乱来视频区| 伦理电影大哥的女人| www.av在线官网国产| 777米奇影视久久| 性色av一级| 精品熟女少妇av免费看| 久久精品久久精品一区二区三区| 新久久久久国产一级毛片| 国产午夜精品久久久久久一区二区三区| 99久久精品国产国产毛片| 精品国产一区二区三区久久久樱花 | 99久久精品热视频| 中文字幕久久专区| 久久久精品欧美日韩精品| 天天一区二区日本电影三级| 日韩一本色道免费dvd| 欧美xxxx性猛交bbbb| 亚洲欧美精品专区久久| 狠狠精品人妻久久久久久综合| 精品一区二区免费观看| 男女无遮挡免费网站观看| a级毛片免费高清观看在线播放| 国产精品不卡视频一区二区| 久久99热这里只频精品6学生| 99久久九九国产精品国产免费| 成人无遮挡网站| 成人毛片60女人毛片免费| 国产美女午夜福利| 精品久久久久久久久亚洲| 精品人妻一区二区三区麻豆| 亚洲欧洲日产国产| 欧美人与善性xxx| 午夜福利视频精品| 人妻夜夜爽99麻豆av| 麻豆精品久久久久久蜜桃| 日韩伦理黄色片| av在线老鸭窝| 亚洲精品国产成人久久av| 国产一区二区三区综合在线观看 | 尤物成人国产欧美一区二区三区| 啦啦啦在线观看免费高清www| 欧美高清成人免费视频www| 黄色怎么调成土黄色| 99精国产麻豆久久婷婷| 51国产日韩欧美| 少妇人妻一区二区三区视频| 在线观看人妻少妇| 欧美国产精品一级二级三级 | 91精品一卡2卡3卡4卡| 日韩大片免费观看网站| 99久久精品一区二区三区| 青青草视频在线视频观看| 亚洲最大成人中文| 热99国产精品久久久久久7| 汤姆久久久久久久影院中文字幕| 亚洲美女搞黄在线观看| 啦啦啦在线观看免费高清www| 熟女av电影| 亚洲精品影视一区二区三区av| 大话2 男鬼变身卡| 又爽又黄无遮挡网站| 亚洲国产av新网站| 精品久久久久久久人妻蜜臀av| 听说在线观看完整版免费高清| 三级经典国产精品| 国产淫片久久久久久久久| 欧美日韩在线观看h| 亚洲美女搞黄在线观看| 国产精品秋霞免费鲁丝片| 亚洲久久久久久中文字幕| 午夜福利高清视频| 午夜福利高清视频| 久久国内精品自在自线图片| 不卡视频在线观看欧美| 国产亚洲午夜精品一区二区久久 | 国产午夜精品久久久久久一区二区三区| 国产成人aa在线观看| 国产女主播在线喷水免费视频网站| 男插女下体视频免费在线播放| av天堂中文字幕网| 久久午夜福利片| 夫妻性生交免费视频一级片| 欧美成人a在线观看| 国产黄片美女视频| 国产免费视频播放在线视频| 欧美人与善性xxx| 青青草视频在线视频观看| 国产毛片在线视频| 观看美女的网站| 日本熟妇午夜| 午夜福利视频1000在线观看| 亚洲精品,欧美精品| 涩涩av久久男人的天堂| 欧美性感艳星| 肉色欧美久久久久久久蜜桃 | 免费av毛片视频| 人妻 亚洲 视频| 午夜福利视频精品| 中文字幕av成人在线电影| 免费黄网站久久成人精品| 精品久久国产蜜桃| 成人亚洲精品一区在线观看 | 91精品国产九色| av又黄又爽大尺度在线免费看| 麻豆成人av视频| 免费观看的影片在线观看| 亚洲欧美精品专区久久| 久久久久久久亚洲中文字幕| 国产精品女同一区二区软件| 夜夜看夜夜爽夜夜摸| 成年免费大片在线观看| 日韩三级伦理在线观看| 搞女人的毛片| 男人爽女人下面视频在线观看| 五月伊人婷婷丁香| 亚洲丝袜综合中文字幕| 91精品伊人久久大香线蕉| 日韩人妻高清精品专区| 国产在线一区二区三区精| 久久精品综合一区二区三区| 亚洲成色77777| 亚洲精品成人av观看孕妇| 免费av毛片视频| 中文资源天堂在线| 春色校园在线视频观看| 亚洲天堂av无毛| 国产免费又黄又爽又色| 一区二区三区四区激情视频| 22中文网久久字幕| 久久精品国产鲁丝片午夜精品| 亚洲精品国产色婷婷电影| 波野结衣二区三区在线| 99热这里只有是精品在线观看| 成人国产av品久久久| 午夜免费鲁丝| 两个人的视频大全免费| 亚洲久久久久久中文字幕| a级毛片免费高清观看在线播放| 夫妻性生交免费视频一级片| 国产视频内射| 国产成人91sexporn| 91狼人影院| 色婷婷久久久亚洲欧美| 国内揄拍国产精品人妻在线| 国产精品人妻久久久久久| 亚洲综合精品二区| 国产 一区精品| 成年女人在线观看亚洲视频 | 欧美精品人与动牲交sv欧美| 久久久久久久久久久丰满| 国产亚洲最大av| 国产精品久久久久久久电影| 欧美成人精品欧美一级黄| 亚洲欧洲日产国产| 最近中文字幕2019免费版| 国产乱来视频区| av又黄又爽大尺度在线免费看| 久久鲁丝午夜福利片| 亚洲人与动物交配视频| 日韩大片免费观看网站| 久久鲁丝午夜福利片| 最近2019中文字幕mv第一页| av免费在线看不卡| 日韩,欧美,国产一区二区三区| 日韩 亚洲 欧美在线| 久久久成人免费电影| 欧美人与善性xxx| 国产人妻一区二区三区在| 久久99热这里只有精品18| 欧美人与善性xxx| 黄色视频在线播放观看不卡| 久久99热这里只有精品18| 久久久久久久久久成人| 在线免费十八禁| 久久久成人免费电影| 内射极品少妇av片p| 国产 精品1| 国产精品久久久久久精品电影| 特级一级黄色大片| 亚洲欧美一区二区三区国产| 最近最新中文字幕免费大全7| 免费大片黄手机在线观看| 一级黄片播放器| 97超视频在线观看视频| 白带黄色成豆腐渣| 欧美丝袜亚洲另类| 亚洲内射少妇av| 国产熟女欧美一区二区| 成人亚洲精品av一区二区| av在线天堂中文字幕|