亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        塘口水電站廠房有限元應(yīng)力分析

        2014-07-24 19:58:47何練民沈建玲
        中國建筑科學(xué) 2014年5期
        關(guān)鍵詞:等值線圖洪水位計算結(jié)果

        何練民+沈建玲

        摘 要:塘口水電站新增裝機10MW,為了準(zhǔn)確進(jìn)行施工圖設(shè)計,保證工程安全前提下盡量節(jié)約投資,進(jìn)行有限元應(yīng)力分析,為施工圖設(shè)計提供依據(jù)。

        關(guān)鍵詞:電站廠房;有限元

        1.工程概況

        塘口水電站位于湖北省來鳳縣漫水鎮(zhèn),距來鳳縣城35km,是酉水河上游干流目前已開發(fā)的第二級梯級電站。酉水流域地處鄂、湘、川三省邊界處,是洞庭湖水系沅水下游左岸的一大支流。酉水發(fā)源于湖北省宣恩縣的城墻巖,源頭呈溶洞、泉水及伏流狀態(tài),自羅谷川開始形成地表河道,自北向南縱貫來鳳縣境,從樂坪開始蜿蜒在鄂湘邊界兩側(cè),于湖南沅陵注入沅水。河流全長407km,流域面積18530km2,河床平均坡降1.1‰。本流域地處武陵山區(qū),屬新構(gòu)造運動上升區(qū),地面海拔高程一般在360~370m之間,區(qū)域內(nèi)氣候溫和,雨量充沛,多年平均降雨量1580.1mm,多年平均徑流深1076mm。

        塘口水電站工程位于酉水河上游河段,壩址控制面積3124km2,是一座以發(fā)電為主的中型水利水電樞紐工程。主要建筑物由混凝土重力壩和廠房組成,水庫總庫容3960萬m3,最大壩高37.60m,原電站裝機15MW,3臺機組,單機5MW,電站最大引用流量96.6m3/s,設(shè)計水頭19.10m,原設(shè)計多年平均發(fā)電量8035萬kW·h。本廠房為擴機增容工程,裝機10MW,為左岸引水式電站。廠房結(jié)構(gòu)形式為壩后式,其主要結(jié)構(gòu)尺寸見圖1-2 。

        2.計算分析基本資料

        本次只對廠房水下部分進(jìn)行了五種工況下的靜力計算,計算時未計入結(jié)構(gòu)的動力響應(yīng)。計算程序采用大型通用有限元軟件Ansys9.0。

        本次計算模型建立了高程374.43m以下(發(fā)電層)部分,其以上部分荷載采用等效荷載法施加至模型上。

        (1)計算工況

        根據(jù)設(shè)計資料,計算共考慮5種工況:

        計算工況1:正常蓄水位:上游水位389.6m;下游水位367.40 m;

        計算工況2:正常運轉(zhuǎn)時突然三相短路。在工況1的基礎(chǔ)上加扭矩205t*m;

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61 m;下游水位384.16 m;

        計算工況4:500年一遇校核洪水位:上游水位394.45 m;下游水位387.94m;

        計算工況5:建成時上下游均無水。

        (2)建模

        本模型建模x方向為上下游方向;Y方向為模型高程方向;z方向為與壩軸線平行方向。模型全部采用實體單元,包括一次單元solid45和二次單元solid92。本模型共包含單元281490個,節(jié)點193730個,單元尺寸控制在0.6m左右。

        (3)材料參數(shù)選取

        根據(jù)設(shè)計資料,模型中共采用了3種材料,其計算參數(shù)取值如下表。

        3.邊界條件及荷載施加

        (1)邊界條件施加

        根據(jù)設(shè)計資料,在所有工況下,邊界約束條件為:模型底部約束Y方向;上下游約束x方向;與壩軸線平行方向約束z方向。

        (2)荷載施加

        計算工況1:正常蓄水位上游389.6;下游367.40??紤]內(nèi)水壓力、發(fā)電機下機蓋部分則轉(zhuǎn)換成面壓力、定子基礎(chǔ)荷載轉(zhuǎn)換成面壓力、高程374.43m以上墻體重量轉(zhuǎn)換成面荷載、揚壓力轉(zhuǎn)換成面荷載、出水口水壓力

        計算工況2:正常運轉(zhuǎn)時突然三相短路??紤]在定子基礎(chǔ)處加扭矩205t*m。加載時,扭矩轉(zhuǎn)換成兩對大小相等,方向相反的力加到定子基礎(chǔ)上,其余荷載同計算工況1。

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61m; 下游水位384.16m。考慮內(nèi)水壓力計算、發(fā)電機下機蓋部分、定子基礎(chǔ)荷載、高程374.43m以上墻體重量、防洪墻受下游水壓力、揚壓力。

        計算工況4:500年一遇校核洪水位:上游水位394.45m,下游水位387.94m。考慮防洪墻受下游水壓力,其余荷載同計算工況3。

        計算工況5:建成時上下游均無水,考慮結(jié)構(gòu)重力效應(yīng)。

        4.計算結(jié)果

        五種工況下計算結(jié)果見下表

        5.工況4拉應(yīng)力等值線圖

        對照計算結(jié)果,列出工況4等值線圖。

        (1)控制截面如下圖所示

        (2)圖號說明:

        圖號2-SZ-Y1-1表示:工況2作用下Y1剖面的Z方向應(yīng)力等值線圖第1張。

        (3)等值線圖:

        X方向應(yīng)力等值線圖:X方向最大拉應(yīng)力為1.15Mpa,最大壓應(yīng)力為2.87Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖5-圖10所示。

        Y方向應(yīng)力等值線圖:Y方向最大拉應(yīng)力為2.81Mpa,最大壓應(yīng)力為3.01Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖11-圖15所示。

        Z方向應(yīng)力等值線圖:Z方向最大拉應(yīng)力為1.06Mpa,最大壓應(yīng)力為1.85Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖16-圖20所示。

        6.結(jié)果分析

        (1)在五種工況下,結(jié)構(gòu)的應(yīng)力最大值為:X方向1.18Mpa(工況1和工況2);Y方向2.81Mpa(工況4); Z方向1.45Mpa(工況2)。

        (2)結(jié)構(gòu)在五種工況下,X方向的應(yīng)力基本相同,且最大應(yīng)力發(fā)生的位置完全相同。這主要是因為該方向最大應(yīng)力發(fā)生位置遠(yuǎn)離水壓作用區(qū),受水壓影響可以忽略,其值的產(chǎn)生主要是重力的作用效應(yīng)。

        (3)計算結(jié)果表明,在工況1和工況2下,結(jié)構(gòu)方向應(yīng)力除Z方向略有不同外,其余基本一致,且應(yīng)力最大位置發(fā)生在同一位置。所以可以把這兩種工況合并為一種考慮。

        (4)工況3和工況4均為校核水位,但比較得知結(jié)構(gòu)在工況4時各方向應(yīng)力均大于工況3對應(yīng)的方向應(yīng)力,因此本報告僅以工況4作計算結(jié)果分析,而不考慮工況3下的應(yīng)力結(jié)果,對工況3的計算結(jié)果僅提供應(yīng)力云圖供設(shè)計參考。

        (5)工況5為結(jié)構(gòu)在自重作用下的效應(yīng),應(yīng)力值很小且不控制設(shè)計。

        (6)按照計算應(yīng)力圖形,作為配筋計算依據(jù)。

        摘 要:塘口水電站新增裝機10MW,為了準(zhǔn)確進(jìn)行施工圖設(shè)計,保證工程安全前提下盡量節(jié)約投資,進(jìn)行有限元應(yīng)力分析,為施工圖設(shè)計提供依據(jù)。

        關(guān)鍵詞:電站廠房;有限元

        1.工程概況

        塘口水電站位于湖北省來鳳縣漫水鎮(zhèn),距來鳳縣城35km,是酉水河上游干流目前已開發(fā)的第二級梯級電站。酉水流域地處鄂、湘、川三省邊界處,是洞庭湖水系沅水下游左岸的一大支流。酉水發(fā)源于湖北省宣恩縣的城墻巖,源頭呈溶洞、泉水及伏流狀態(tài),自羅谷川開始形成地表河道,自北向南縱貫來鳳縣境,從樂坪開始蜿蜒在鄂湘邊界兩側(cè),于湖南沅陵注入沅水。河流全長407km,流域面積18530km2,河床平均坡降1.1‰。本流域地處武陵山區(qū),屬新構(gòu)造運動上升區(qū),地面海拔高程一般在360~370m之間,區(qū)域內(nèi)氣候溫和,雨量充沛,多年平均降雨量1580.1mm,多年平均徑流深1076mm。

        塘口水電站工程位于酉水河上游河段,壩址控制面積3124km2,是一座以發(fā)電為主的中型水利水電樞紐工程。主要建筑物由混凝土重力壩和廠房組成,水庫總庫容3960萬m3,最大壩高37.60m,原電站裝機15MW,3臺機組,單機5MW,電站最大引用流量96.6m3/s,設(shè)計水頭19.10m,原設(shè)計多年平均發(fā)電量8035萬kW·h。本廠房為擴機增容工程,裝機10MW,為左岸引水式電站。廠房結(jié)構(gòu)形式為壩后式,其主要結(jié)構(gòu)尺寸見圖1-2 。

        2.計算分析基本資料

        本次只對廠房水下部分進(jìn)行了五種工況下的靜力計算,計算時未計入結(jié)構(gòu)的動力響應(yīng)。計算程序采用大型通用有限元軟件Ansys9.0。

        本次計算模型建立了高程374.43m以下(發(fā)電層)部分,其以上部分荷載采用等效荷載法施加至模型上。

        (1)計算工況

        根據(jù)設(shè)計資料,計算共考慮5種工況:

        計算工況1:正常蓄水位:上游水位389.6m;下游水位367.40 m;

        計算工況2:正常運轉(zhuǎn)時突然三相短路。在工況1的基礎(chǔ)上加扭矩205t*m;

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61 m;下游水位384.16 m;

        計算工況4:500年一遇校核洪水位:上游水位394.45 m;下游水位387.94m;

        計算工況5:建成時上下游均無水。

        (2)建模

        本模型建模x方向為上下游方向;Y方向為模型高程方向;z方向為與壩軸線平行方向。模型全部采用實體單元,包括一次單元solid45和二次單元solid92。本模型共包含單元281490個,節(jié)點193730個,單元尺寸控制在0.6m左右。

        (3)材料參數(shù)選取

        根據(jù)設(shè)計資料,模型中共采用了3種材料,其計算參數(shù)取值如下表。

        3.邊界條件及荷載施加

        (1)邊界條件施加

        根據(jù)設(shè)計資料,在所有工況下,邊界約束條件為:模型底部約束Y方向;上下游約束x方向;與壩軸線平行方向約束z方向。

        (2)荷載施加

        計算工況1:正常蓄水位上游389.6;下游367.40??紤]內(nèi)水壓力、發(fā)電機下機蓋部分則轉(zhuǎn)換成面壓力、定子基礎(chǔ)荷載轉(zhuǎn)換成面壓力、高程374.43m以上墻體重量轉(zhuǎn)換成面荷載、揚壓力轉(zhuǎn)換成面荷載、出水口水壓力

        計算工況2:正常運轉(zhuǎn)時突然三相短路??紤]在定子基礎(chǔ)處加扭矩205t*m。加載時,扭矩轉(zhuǎn)換成兩對大小相等,方向相反的力加到定子基礎(chǔ)上,其余荷載同計算工況1。

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61m; 下游水位384.16m。考慮內(nèi)水壓力計算、發(fā)電機下機蓋部分、定子基礎(chǔ)荷載、高程374.43m以上墻體重量、防洪墻受下游水壓力、揚壓力。

        計算工況4:500年一遇校核洪水位:上游水位394.45m,下游水位387.94m。考慮防洪墻受下游水壓力,其余荷載同計算工況3。

        計算工況5:建成時上下游均無水,考慮結(jié)構(gòu)重力效應(yīng)。

        4.計算結(jié)果

        五種工況下計算結(jié)果見下表

        5.工況4拉應(yīng)力等值線圖

        對照計算結(jié)果,列出工況4等值線圖。

        (1)控制截面如下圖所示

        (2)圖號說明:

        圖號2-SZ-Y1-1表示:工況2作用下Y1剖面的Z方向應(yīng)力等值線圖第1張。

        (3)等值線圖:

        X方向應(yīng)力等值線圖:X方向最大拉應(yīng)力為1.15Mpa,最大壓應(yīng)力為2.87Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖5-圖10所示。

        Y方向應(yīng)力等值線圖:Y方向最大拉應(yīng)力為2.81Mpa,最大壓應(yīng)力為3.01Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖11-圖15所示。

        Z方向應(yīng)力等值線圖:Z方向最大拉應(yīng)力為1.06Mpa,最大壓應(yīng)力為1.85Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖16-圖20所示。

        6.結(jié)果分析

        (1)在五種工況下,結(jié)構(gòu)的應(yīng)力最大值為:X方向1.18Mpa(工況1和工況2);Y方向2.81Mpa(工況4); Z方向1.45Mpa(工況2)。

        (2)結(jié)構(gòu)在五種工況下,X方向的應(yīng)力基本相同,且最大應(yīng)力發(fā)生的位置完全相同。這主要是因為該方向最大應(yīng)力發(fā)生位置遠(yuǎn)離水壓作用區(qū),受水壓影響可以忽略,其值的產(chǎn)生主要是重力的作用效應(yīng)。

        (3)計算結(jié)果表明,在工況1和工況2下,結(jié)構(gòu)方向應(yīng)力除Z方向略有不同外,其余基本一致,且應(yīng)力最大位置發(fā)生在同一位置。所以可以把這兩種工況合并為一種考慮。

        (4)工況3和工況4均為校核水位,但比較得知結(jié)構(gòu)在工況4時各方向應(yīng)力均大于工況3對應(yīng)的方向應(yīng)力,因此本報告僅以工況4作計算結(jié)果分析,而不考慮工況3下的應(yīng)力結(jié)果,對工況3的計算結(jié)果僅提供應(yīng)力云圖供設(shè)計參考。

        (5)工況5為結(jié)構(gòu)在自重作用下的效應(yīng),應(yīng)力值很小且不控制設(shè)計。

        (6)按照計算應(yīng)力圖形,作為配筋計算依據(jù)。

        摘 要:塘口水電站新增裝機10MW,為了準(zhǔn)確進(jìn)行施工圖設(shè)計,保證工程安全前提下盡量節(jié)約投資,進(jìn)行有限元應(yīng)力分析,為施工圖設(shè)計提供依據(jù)。

        關(guān)鍵詞:電站廠房;有限元

        1.工程概況

        塘口水電站位于湖北省來鳳縣漫水鎮(zhèn),距來鳳縣城35km,是酉水河上游干流目前已開發(fā)的第二級梯級電站。酉水流域地處鄂、湘、川三省邊界處,是洞庭湖水系沅水下游左岸的一大支流。酉水發(fā)源于湖北省宣恩縣的城墻巖,源頭呈溶洞、泉水及伏流狀態(tài),自羅谷川開始形成地表河道,自北向南縱貫來鳳縣境,從樂坪開始蜿蜒在鄂湘邊界兩側(cè),于湖南沅陵注入沅水。河流全長407km,流域面積18530km2,河床平均坡降1.1‰。本流域地處武陵山區(qū),屬新構(gòu)造運動上升區(qū),地面海拔高程一般在360~370m之間,區(qū)域內(nèi)氣候溫和,雨量充沛,多年平均降雨量1580.1mm,多年平均徑流深1076mm。

        塘口水電站工程位于酉水河上游河段,壩址控制面積3124km2,是一座以發(fā)電為主的中型水利水電樞紐工程。主要建筑物由混凝土重力壩和廠房組成,水庫總庫容3960萬m3,最大壩高37.60m,原電站裝機15MW,3臺機組,單機5MW,電站最大引用流量96.6m3/s,設(shè)計水頭19.10m,原設(shè)計多年平均發(fā)電量8035萬kW·h。本廠房為擴機增容工程,裝機10MW,為左岸引水式電站。廠房結(jié)構(gòu)形式為壩后式,其主要結(jié)構(gòu)尺寸見圖1-2 。

        2.計算分析基本資料

        本次只對廠房水下部分進(jìn)行了五種工況下的靜力計算,計算時未計入結(jié)構(gòu)的動力響應(yīng)。計算程序采用大型通用有限元軟件Ansys9.0。

        本次計算模型建立了高程374.43m以下(發(fā)電層)部分,其以上部分荷載采用等效荷載法施加至模型上。

        (1)計算工況

        根據(jù)設(shè)計資料,計算共考慮5種工況:

        計算工況1:正常蓄水位:上游水位389.6m;下游水位367.40 m;

        計算工況2:正常運轉(zhuǎn)時突然三相短路。在工況1的基礎(chǔ)上加扭矩205t*m;

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61 m;下游水位384.16 m;

        計算工況4:500年一遇校核洪水位:上游水位394.45 m;下游水位387.94m;

        計算工況5:建成時上下游均無水。

        (2)建模

        本模型建模x方向為上下游方向;Y方向為模型高程方向;z方向為與壩軸線平行方向。模型全部采用實體單元,包括一次單元solid45和二次單元solid92。本模型共包含單元281490個,節(jié)點193730個,單元尺寸控制在0.6m左右。

        (3)材料參數(shù)選取

        根據(jù)設(shè)計資料,模型中共采用了3種材料,其計算參數(shù)取值如下表。

        3.邊界條件及荷載施加

        (1)邊界條件施加

        根據(jù)設(shè)計資料,在所有工況下,邊界約束條件為:模型底部約束Y方向;上下游約束x方向;與壩軸線平行方向約束z方向。

        (2)荷載施加

        計算工況1:正常蓄水位上游389.6;下游367.40??紤]內(nèi)水壓力、發(fā)電機下機蓋部分則轉(zhuǎn)換成面壓力、定子基礎(chǔ)荷載轉(zhuǎn)換成面壓力、高程374.43m以上墻體重量轉(zhuǎn)換成面荷載、揚壓力轉(zhuǎn)換成面荷載、出水口水壓力

        計算工況2:正常運轉(zhuǎn)時突然三相短路??紤]在定子基礎(chǔ)處加扭矩205t*m。加載時,扭矩轉(zhuǎn)換成兩對大小相等,方向相反的力加到定子基礎(chǔ)上,其余荷載同計算工況1。

        計算工況3:50年一遇設(shè)計洪水位:上游水位391.61m; 下游水位384.16m??紤]內(nèi)水壓力計算、發(fā)電機下機蓋部分、定子基礎(chǔ)荷載、高程374.43m以上墻體重量、防洪墻受下游水壓力、揚壓力。

        計算工況4:500年一遇校核洪水位:上游水位394.45m,下游水位387.94m。考慮防洪墻受下游水壓力,其余荷載同計算工況3。

        計算工況5:建成時上下游均無水,考慮結(jié)構(gòu)重力效應(yīng)。

        4.計算結(jié)果

        五種工況下計算結(jié)果見下表

        5.工況4拉應(yīng)力等值線圖

        對照計算結(jié)果,列出工況4等值線圖。

        (1)控制截面如下圖所示

        (2)圖號說明:

        圖號2-SZ-Y1-1表示:工況2作用下Y1剖面的Z方向應(yīng)力等值線圖第1張。

        (3)等值線圖:

        X方向應(yīng)力等值線圖:X方向最大拉應(yīng)力為1.15Mpa,最大壓應(yīng)力為2.87Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖5-圖10所示。

        Y方向應(yīng)力等值線圖:Y方向最大拉應(yīng)力為2.81Mpa,最大壓應(yīng)力為3.01Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖11-圖15所示。

        Z方向應(yīng)力等值線圖:Z方向最大拉應(yīng)力為1.06Mpa,最大壓應(yīng)力為1.85Mpa。截取通過最大應(yīng)力點的三個方向的剖面,其應(yīng)力等值線圖如圖16-圖20所示。

        6.結(jié)果分析

        (1)在五種工況下,結(jié)構(gòu)的應(yīng)力最大值為:X方向1.18Mpa(工況1和工況2);Y方向2.81Mpa(工況4); Z方向1.45Mpa(工況2)。

        (2)結(jié)構(gòu)在五種工況下,X方向的應(yīng)力基本相同,且最大應(yīng)力發(fā)生的位置完全相同。這主要是因為該方向最大應(yīng)力發(fā)生位置遠(yuǎn)離水壓作用區(qū),受水壓影響可以忽略,其值的產(chǎn)生主要是重力的作用效應(yīng)。

        (3)計算結(jié)果表明,在工況1和工況2下,結(jié)構(gòu)方向應(yīng)力除Z方向略有不同外,其余基本一致,且應(yīng)力最大位置發(fā)生在同一位置。所以可以把這兩種工況合并為一種考慮。

        (4)工況3和工況4均為校核水位,但比較得知結(jié)構(gòu)在工況4時各方向應(yīng)力均大于工況3對應(yīng)的方向應(yīng)力,因此本報告僅以工況4作計算結(jié)果分析,而不考慮工況3下的應(yīng)力結(jié)果,對工況3的計算結(jié)果僅提供應(yīng)力云圖供設(shè)計參考。

        (5)工況5為結(jié)構(gòu)在自重作用下的效應(yīng),應(yīng)力值很小且不控制設(shè)計。

        (6)按照計算應(yīng)力圖形,作為配筋計算依據(jù)。

        猜你喜歡
        等值線圖洪水位計算結(jié)果
        湖南省似大地水準(zhǔn)面模型對比研究
        不等高軟橫跨橫向承力索計算及計算結(jié)果判斷研究
        甘肅科技(2020年20期)2020-04-13 00:30:40
        湖南陵水某堤防工程風(fēng)險分析方法研究
        如何來解決等值線問題
        沙雅河護(hù)岸工程防洪堤設(shè)計洪水位的確定
        利用ACAD線型編輯功能實現(xiàn)SURFER繪制等值線圖的修改初探
        江蘇水利(2015年9期)2015-12-12 08:17:18
        超壓測試方法對炸藥TNT當(dāng)量計算結(jié)果的影響
        噪聲對介質(zhì)損耗角正切計算結(jié)果的影響
        水庫相關(guān)水位、庫容(之二)
        Effects of electroacupuncture at Zusanli (ST 36) on neurons in the colonic myenteric plexus in rats with irritable bowel syndrome with constipation*★
        亚洲国产日韩欧美高清片a| 亚洲av一二三区成人影片| 比较有韵味的熟妇无码| 久久精品国产亚洲av大全| 精品中文字幕久久久久久| 日韩av不卡一二三区| 亚洲国产一区二区中文字幕| 新婚少妇无套内谢国语播放| 又大又粗又爽18禁免费看| 国产黄a三级三级三级av在线看| 亚洲av永久无码精品水牛影视| 在线视频免费自拍亚洲| 亚洲伦理第一页中文字幕| 宅男天堂亚洲一区二区三区| 日本伦理精品一区二区三区 | 99无码熟妇丰满人妻啪啪| 国产熟女露脸大叫高潮| 视频在线观看一区二区三区| 亚洲视频在线视频在线视频| 黄页免费人成网址大全| 国产亚洲超级97免费视频| 日韩精品一区二区三区中文| 亚洲色在线v中文字幕| 中文字幕无码人妻丝袜| 少妇极品熟妇人妻高清| 久久精品伊人久久精品| 人妻少妇进入猛烈时中文字幕| 成年女人粗暴毛片免费观看| 一二区视频免费在线观看| 亚洲av无码之国产精品网址蜜芽| 欧美粗大无套gay| 国产一区二区欧美丝袜| 精品人妻一区二区三区不卡毛片| 国产日本精品一二三四区| 久久精品国产久精国产果冻传媒 | 久久99热久久99精品| 亚洲区偷拍自拍29p| 精品人妻一区二区视频| 香蕉成人伊视频在线观看| 乱色熟女综合一区二区三区| 久久精品国产色蜜蜜麻豆|