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

    球柱殼耐壓艙體極限承載力研究

    2022-04-02 06:07:50張洪彬徐會希
    海洋工程 2022年2期
    關(guān)鍵詞:承載力結(jié)構(gòu)模型

    張洪彬,徐會希

    (1.中國科學(xué)院沈陽自動化研究所 機器人學(xué)國家重點實驗室,遼寧 沈陽 110016; 2.中國科學(xué)院 機器人與智能制造創(chuàng)新研究院,遼寧 沈陽 110169; 3.遼寧省水下機器人重點實驗室,遼寧 沈陽 110169)

    耐壓作為深潛技術(shù)的重要環(huán)節(jié),是決定水下機器人[1-2]下潛深度的關(guān)鍵。而大壓力條件下的耐壓艙體結(jié)構(gòu)設(shè)計,一直是直接關(guān)乎水下機器人安全性與總體指標的重要研究內(nèi)容。但是單純依賴線性靜強度失效計算,即假設(shè)耐壓艙體結(jié)構(gòu)不會發(fā)生失穩(wěn)破壞,遠遠不能滿足設(shè)計安全性要求。經(jīng)驗表明,大深度鈦合金薄殼及中厚殼類耐壓艙體結(jié)構(gòu),其破壞形式往往是先發(fā)生失穩(wěn)失效,即極限承載力是決定水下機器人大深度下潛的核心指標。

    有限元的特征值屈曲分析只能計算出理想結(jié)構(gòu)的臨界失穩(wěn)載荷,結(jié)果非保守。針對球殼結(jié)構(gòu),傳統(tǒng)經(jīng)驗公式的計算誤差一般可以比試驗值大3倍左右,參考價值更小。圍繞極限承載力的計算,國內(nèi)外已有一些科研成果產(chǎn)出[3-9],但目前的研究成果主要集中在針對載人潛水器的球殼體結(jié)構(gòu),針對無人潛水器常用的球殼封蓋加圓柱形耐壓殼體結(jié)構(gòu)的極限承載力研究少之又少。

    文中研究基于有限元理論的非線性屈曲分析方法,為了提升計算的準確性,首先采用裝配體整體建模策略,搭建復(fù)合屈曲失穩(wěn)的初始受力模型。然后在同時考慮材料非線性、幾何非線性、接觸非線性等方面因素的條件下,評估承壓結(jié)構(gòu)的最大極限承載力。最后針對初始缺陷導(dǎo)入尺度這一關(guān)鍵參數(shù)展開深入研究,結(jié)合試驗結(jié)果驗證,確定柱形耐壓殼體結(jié)構(gòu)的初始缺陷導(dǎo)入尺度參數(shù),為未來的耐壓艙體結(jié)構(gòu)設(shè)計、計算分析及壓力試驗提供試驗指導(dǎo)與技術(shù)支撐。

    1 傳統(tǒng)經(jīng)驗公式計算方法

    1.1 耐壓球殼的經(jīng)典失穩(wěn)公式

    早在1915年,R.Zoelly就基于小變形假設(shè),導(dǎo)出受外壓的耐壓球殼經(jīng)典失穩(wěn)公式[10],這是最早的穩(wěn)定性校核公式。

    (1)

    其中,pE為臨界失穩(wěn)壓力,E為彈性模量,μ為泊松比,t為球殼厚度,R為球殼中面半徑。此經(jīng)典公式包含4個基本假設(shè):1)對材料各向均勻同性;2)有完善的幾何球形;3)無初始內(nèi)部應(yīng)力;4)滿足線彈性本構(gòu)方程。

    該經(jīng)驗公式的局限性和誤差來源有:幾何形狀的完善只能是理想化球殼結(jié)構(gòu),無法考慮結(jié)構(gòu)初始缺陷和結(jié)合非線性的影響;基于線彈性的基本假設(shè),忽略了材料和幾何非線性及內(nèi)部殘余應(yīng)力等因素,使得分析過程過于簡化。

    當球形殼體承受均勻外壓時,球殼的計算應(yīng)力有:

    (2)

    其中,p為外壓力。

    將式(2)代入式(1),臨界失穩(wěn)應(yīng)力有:

    (3)

    1.2 BERCH公式

    在經(jīng)典公式的基礎(chǔ)上,結(jié)合許多學(xué)者的后續(xù)實踐,又有新的經(jīng)驗公式產(chǎn)生。代表性的有BERCH公式:

    (4)

    其中,σs為屈服強度,pcr為殼體的臨界壓力。

    1.3 德川公式

    基于20面體近似球殼推導(dǎo)的德川公式為:

    (5)

    其中,D為球殼直徑。

    將理論模型與試驗?zāi)P瓦M一步結(jié)合后,根據(jù)厚度直徑比的不同具體有:

    (6)

    (7)

    上述經(jīng)驗公式都在經(jīng)典理論的基礎(chǔ)上,結(jié)合模型的試驗結(jié)果進行一定程度的優(yōu)化,式(4)引入并適當考慮了材料的特性,式(6)和(7)針對厚徑比進一步將公式詳細區(qū)分,但是都沒有考慮初始結(jié)構(gòu)缺陷對于非線性穩(wěn)定性的影響。

    1.4 卡門—錢學(xué)森公式

    考慮球殼微小變形的能量法求解得到卡門—錢學(xué)森公式,公式考慮了結(jié)構(gòu)初始缺陷對于結(jié)構(gòu)穩(wěn)定性的影響。

    (8)

    1.5 美國海軍泰勒水池耐壓殼設(shè)計公式

    美國海軍泰勒水池耐壓殼設(shè)計公式為:

    (9)

    1.6 俄羅斯深海潛水器耐壓殼體的臨界壓力公式

    俄羅斯深海潛水器耐壓殼體的臨界壓力公式為:

    (10)

    1.7 經(jīng)驗公式總結(jié)

    已有相關(guān)學(xué)者對于上述經(jīng)驗公式進行研究,對上述部分經(jīng)驗公式的計算結(jié)果進行對比評價,這里不做深入分析。只給出基本結(jié)論:經(jīng)典計算公式和部分經(jīng)驗公式,由于沒有考慮材料非線性、幾何非線性、接觸非線性和初始幾何缺陷諸多因素中的一種或者多種,使得計算結(jié)果誤差較大,工程指導(dǎo)意義有限。

    美國海軍泰勒水池耐壓殼設(shè)計公式和俄羅斯深海潛水器耐壓殼體的臨界壓力公式,充分考慮結(jié)構(gòu)局部缺陷和材料特性對于最大極限承載力的影響,計算公式經(jīng)過大量試驗數(shù)據(jù)修正,有較強的參考價值。但是由于試驗和公式都是基于球殼或者類球殼結(jié)構(gòu)完成的,同時試驗所用材料在漫長的發(fā)展過程中其力學(xué)性能參數(shù)也有較大差異,所以針對鈦合金材料廣泛采用的圓柱筒體加球殼封蓋式承壓結(jié)構(gòu),經(jīng)驗公式計算方法的指導(dǎo)價值有很大的局限性。

    2 非線性屈曲分析流程

    2.1 不同類型缺陷結(jié)構(gòu)分析流程

    非線性屈曲分析,可以細分為確定性缺陷模型分析方法和非確定性缺陷模型分析方法兩種。針對兩種模型建立分析流程如圖1所示。其中確定性缺陷模型分析方法針對已知可測量的缺陷模型,進行三維缺陷建模后開展極限承載力分析;非確定性缺陷模型分析方法針對缺陷無法精確測定的幾何結(jié)構(gòu),借助特征值屈曲的模態(tài)陣型進行缺陷導(dǎo)入后分析。

    圖1 非線性屈曲分析流程

    2.2 材料非線性特性

    極限承載力分析技術(shù)需要充分考慮材料的塑性應(yīng)力應(yīng)變特性,通過力學(xué)試驗測定鈦合金TC4棒料的材料應(yīng)力應(yīng)變曲線和塑性應(yīng)力應(yīng)變曲線分別如圖2和3所示。鋁合金7075棒料的塑性應(yīng)力應(yīng)變曲線如圖4所示。

    圖3 鈦合金TC4單軸拉伸試驗塑性應(yīng)力—應(yīng)變曲線

    圖4 鋁合金7075單軸拉伸試驗塑性應(yīng)力—應(yīng)變曲線

    2.3 幾何非線性特性

    當結(jié)構(gòu)的變形使體系的受力發(fā)生顯著變化,以至不能采用線性體系的分析方法時就稱為幾何非線性,即力與位移不再是直線關(guān)系。常見分類有大位移小應(yīng)變問題、大位移大應(yīng)變問題、大轉(zhuǎn)角問題。直觀理解為當結(jié)構(gòu)的位移相對于結(jié)構(gòu)本體的尺度不再是高階小量時,需要充分考慮幾何非線性問題。

    2.4 接觸非線性特性

    接觸是狀態(tài)改變的非線性,系統(tǒng)剛度隨接觸狀態(tài)的變化而實時更新,是一個高度非線性的過程。對于研究的球殼封蓋加圓柱形耐壓殼體的耐壓艙體結(jié)構(gòu),半球殼與圓柱筒體在接觸面之間力的精確傳遞也是影響全局結(jié)構(gòu)極限承載力的關(guān)鍵環(huán)節(jié)。同時由于球殼與圓柱殼體耐壓能力的差別,等壁厚條件下交界面上的圓柱筒體邊緣剛度遠遠小于球殼體邊緣剛度。所以如果采用零件一體化建?;蛘吖灿霉?jié)點的方式處理,會使得球殼邊界將圓柱殼體邊界進一步強化,邊界傳力不精確,計算結(jié)果誤差偏大。

    所以充分考慮耐壓艙體密封形式及其深海環(huán)境下的實際工作受力狀態(tài),采用裝配體整體建模策略,結(jié)合面采用非線性摩擦接觸建模,將力精確傳遞,同時充分考慮結(jié)合面上不同零件結(jié)構(gòu)剛度的差異。但接觸的引入往往會使得極限承載力的求解收斂變得異常困難,所以文中引入穩(wěn)定能這一求解技術(shù)解決這個問題,使得這一超強非線性的極限承載力分析能夠排除數(shù)值不穩(wěn)定的失穩(wěn)點,獲得物理不穩(wěn)定的數(shù)值解。

    3 初始缺陷的確定

    3.1 非確定性缺陷模型

    針對非確定性缺陷模型,此類缺陷主要針對設(shè)計階段的結(jié)構(gòu),以及無法測定缺陷大小和尺寸的結(jié)構(gòu)。此類缺陷的引入需要通過線性屈曲分析求解出屈曲模態(tài),計算中的陣型需要選擇對于結(jié)構(gòu)最為不利的一階陣型,該階陣型的形態(tài)即為在此結(jié)構(gòu)振動方向上剛度最弱,最容易發(fā)生失效的形態(tài)。對于球殼類結(jié)構(gòu),因其結(jié)構(gòu)形式具有先天的高抗壓能力,經(jīng)驗的屈曲模態(tài)常取第一階陣型。對于圓柱筒體類結(jié)構(gòu),需要通過輸出前6階模態(tài)陣型后進行人工選擇,挑選最不利的一階陣型進行導(dǎo)入。

    由于屈曲模態(tài)陣型的數(shù)學(xué)含義只是一個無量綱的歸一化特征向量,即沒有表征缺陷大小的含義,所以需要將陣型向量的最大振幅值放大一定的倍數(shù)。這個倍數(shù)叫做屈曲形狀導(dǎo)入因子,其具有長度上的量綱,大小為屈曲形狀缺陷(長度量綱)與屈曲模態(tài)陣型最大振幅(無量綱)的比值。

    針對屈曲形狀缺陷的取值,不同的文獻資料有所不同。這里只列舉相關(guān)的參考標準及具體的屈曲形狀缺陷取值。以某型無人潛水器的耐壓艙耐壓筒體為例進行詳細研究,針對不同的屈曲形狀缺陷0.115 mm、0.274 2 mm、0.6 mm分別計算,并給出最大極限承載力。

    表1 屈曲形狀缺陷導(dǎo)入值

    3.2 非線性屈曲有限元模型網(wǎng)格要求

    非線性屈曲分析,是一項要考慮材料非線性、幾何非線性和邊界(接觸)非線性要素,同時加入幾何初始缺陷的高度非線性、復(fù)雜有限元分析,這樣的高級分析計算需要高質(zhì)量的網(wǎng)格劃分技術(shù)支撐。具體的單元選擇與網(wǎng)格質(zhì)量要求為:

    1)采用多節(jié)點高精度單元對結(jié)構(gòu)進行計算;

    2)單元類型殼體單元盡量選擇四邊形殼體單元離散,實體單元盡量選擇六面體單元離散;

    3)網(wǎng)格質(zhì)量盡量劃分到0.7以上,最好不要低于0.5,如果存在局部不規(guī)則特征,可以嘗試局部細化手段處理。

    4 非線性屈曲具體算例校核

    首先針對帶有確定性缺陷的厚殼耐壓艙體結(jié)構(gòu)進行極限承載力分析,針對某型潛水器的鋁合金7075材質(zhì)厚殼耐壓艙體,對其進行表面缺陷尺寸的測繪,將此測繪的缺陷建模埋入到幾何模型中,進行極限承載力分析;將缺陷去除,即使用理想無缺陷厚殼模型進行極限承載力分析;再將計算結(jié)果與基于船級社標準[11]缺陷容差模型的最大極限承載力對比分析。針對某型潛水器的薄殼鈦合金材質(zhì)耐壓艙體模型,不同缺陷尺度的導(dǎo)入標準,分別進行非確定性缺陷最大極限承載力分析。計算得出相應(yīng)的穩(wěn)定性安全系數(shù),對相應(yīng)的耐壓艙安全性進行評估。

    4.1 確定性缺陷算例

    4.1.1 某型潛水器點蝕缺陷厚殼耐壓艙極限承載力校核

    建立某型潛水器耐壓筒體模型,此耐壓艙圓柱筒體20 mm等壁厚,球蓋頂部最薄厚度10.2 mm。預(yù)埋缺陷尺寸為R=5.7 mm,深度2.83 mm,兩個位于中心位置,另外一個沿母線方向80 mm對稱布置,具體如圖5所示,圖中為二分之一對稱建模,只顯示兩個缺陷凹坑。

    圖5 確定性缺陷三維建模

    所加邊界約束條件如圖6所示,約束對稱面上的x(法向)位移和球蓋頂部的切平面內(nèi)自由度,即釋放軸向移動自由度。耐壓艙體接觸面施加摩擦系數(shù)0.1的摩擦約束。后續(xù)章節(jié)邊界條件與本章節(jié)相同,不做贅述。極限承載載荷為89 MPa。

    圖6 邊界條件

    4.1.2 某型潛水器無蝕缺陷厚殼耐壓艙極限承載力校核

    使用4.1.1節(jié)中去除缺陷后的模型,即建立無缺陷理想化耐壓艙有限元模型,按照確定性缺陷的極限承載力流程計算,獲得最大極限承載力為88.8 MPa,與帶有小尺度點蝕凹坑缺陷的耐壓艙體結(jié)構(gòu)最大極限承載力相差不大。

    4.1.3 非確定性缺陷厚殼耐壓艙體結(jié)構(gòu)極限承載力分析

    使用非確定性缺陷模型屈曲計算方法,執(zhí)行船級社規(guī)范[11]要求,缺陷尺度為中性層半徑的0.5%,中性層半徑為125 mm,導(dǎo)入缺陷為0.625 mm。進一步對比評估0.5%的整體模態(tài)屈曲陣型缺陷對于某型耐壓艙體最大極限承載力的影響情況。經(jīng)過計算獲得臨界失穩(wěn)壓力為79 MPa。

    4.1.4 小結(jié)

    1)帶有2.83 mm深度的確定性缺陷厚殼耐壓艙體極限承載力為89 MPa,理想無缺陷某型厚殼耐壓艙體臨界失穩(wěn)壓力為88.8 MPa,由此可得該尺度點蝕類缺陷對于厚殼結(jié)構(gòu)的極限承載力影響微乎其微,工程上可以忽略處理。

    2)真球度和圓柱度等形狀缺陷和耐壓殼體厚度對于極限承載的影響權(quán)重較高。

    4.2 非確定性缺陷算例

    針對某型薄殼耐壓艙體,分別進行屈曲形狀導(dǎo)入缺陷為0.115 mm、0.274 2 mm、0.6 mm 時的最大極限承載力計算。線性屈曲最危險陣型云圖如圖7所示。

    圖7 耐壓艙體最危險線性屈曲陣型云圖

    由圖7可知屈曲模態(tài)陣型最大振幅為1.000 3,所以屈曲缺陷導(dǎo)入因子為0.115 mm/1.000 3=0.115 0 mm、0.274 2 mm/1.000 3=0.274 1 mm、0.6 mm/1.000 3=0.599 8 mm。將會以此導(dǎo)入因子借助軟件命令流將模態(tài)缺陷輸入到計算模型中。不同模態(tài)缺陷導(dǎo)入尺度與極限承載力的關(guān)系如表2所示。

    表2 模態(tài)缺陷導(dǎo)入尺度與極限承載力關(guān)系

    為充分驗證極限承載力計算的準確性,在中國科學(xué)院沈陽自動化研究所水下試驗室開展水壓爆破試驗,將上述經(jīng)三坐標檢測的耐壓艙體放入壓力罐內(nèi)開展爆破試驗,在緩慢升壓至54 MPa后保壓約5 min時,發(fā)生耐壓艙艙體的壓潰碎裂、球殼的變形、法蘭的斷裂,如圖8所示。

    圖8 某型耐壓艙壓力試驗后照片

    艙體的三坐標檢測缺陷最大值檢測結(jié)果詳見表1第二組數(shù)據(jù),實測的缺陷參數(shù)宏觀表現(xiàn)為整體的不圓度,最大值為艙體結(jié)構(gòu)的0.23%,這與有限元中的模態(tài)缺陷定義契合。因此,通過水壓爆破試驗,驗證了極限承載力計算結(jié)果準確可靠。

    5 結(jié) 語

    在裝配體整體建模策略的基礎(chǔ)上,同時考慮材料非線性、幾何非線性、接觸非線性等方面因素,評估承壓結(jié)構(gòu)的最大極限承載力。并針對初始缺陷導(dǎo)入尺度這一關(guān)鍵參數(shù)展開深入研究,研究得出不同類型的缺陷對承壓結(jié)構(gòu)的削弱有所不同。最后給出缺陷導(dǎo)入?yún)?shù)建議值,并對球柱殼體耐壓艙的結(jié)構(gòu)分析提出要求。

    1)帶有2.83 mm深度的確定性缺陷厚殼耐壓艙體極限承載力為89 MPa,理想無缺陷某型厚殼耐壓艙體臨界失穩(wěn)壓力為88.8 MPa,二者相差0.225%,因含有確定性缺陷的艙體結(jié)構(gòu)無法全局采用結(jié)構(gòu)化網(wǎng)格,由此引入一定的數(shù)值計算誤差,但計算結(jié)果差異較小。由此可得該尺度點蝕缺陷對于厚殼結(jié)構(gòu)的極限承載力影響較小,工程上可以忽略處理。

    2)真球度和圓柱度等形狀缺陷和耐壓殼體厚度對于極限承載的影響權(quán)重較高。

    3)針對無人潛水器的球殼封蓋加圓柱形耐壓殼體結(jié)構(gòu),0.5%的缺陷尺度對于結(jié)構(gòu)的穩(wěn)定性削弱過大,存在過安全的問題,對于耐壓艙體數(shù)量眾多的無人潛水器而言,過大的質(zhì)量總體指標也無法接受。

    4)基于研究成果,球殼封蓋加圓柱形耐壓殼體結(jié)構(gòu),缺陷導(dǎo)入尺度為0.23%,該缺陷導(dǎo)入尺度即為現(xiàn)有結(jié)構(gòu)臨界失穩(wěn)極限缺陷值。

    5)鈦合金類薄殼耐壓艙體結(jié)構(gòu)分析的非線性屈曲分析重要性遠遠大于靜強度失效校核分析。所以一個嚴謹?shù)膹姸刃:朔治霰仨毷庆o強度失效校核與非線性屈曲分析兩相結(jié)合。

    猜你喜歡
    承載力結(jié)構(gòu)模型
    一半模型
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    論《日出》的結(jié)構(gòu)
    3D打印中的模型分割與打包
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    創(chuàng)新治理結(jié)構(gòu)促進中小企業(yè)持續(xù)成長
    国产高清有码在线观看视频 | 91九色精品人成在线观看| 又紧又爽又黄一区二区| 亚洲国产精品sss在线观看| 久久精品影院6| 久久天躁狠狠躁夜夜2o2o| 成年免费大片在线观看| 91在线观看av| 男人舔女人的私密视频| 天堂√8在线中文| 亚洲精品久久国产高清桃花| 午夜久久久在线观看| 国产亚洲精品一区二区www| 真人做人爱边吃奶动态| 曰老女人黄片| 草草在线视频免费看| 18禁美女被吸乳视频| 很黄的视频免费| 国产99久久九九免费精品| 国产视频一区二区在线看| 在线观看一区二区三区| 满18在线观看网站| 欧美中文综合在线视频| 老司机在亚洲福利影院| 1024香蕉在线观看| 欧美成狂野欧美在线观看| 精品久久久久久久久久久久久 | 黑丝袜美女国产一区| 操出白浆在线播放| 男人舔奶头视频| 亚洲精品粉嫩美女一区| 女警被强在线播放| 女人爽到高潮嗷嗷叫在线视频| 日本免费一区二区三区高清不卡| 美女午夜性视频免费| 一a级毛片在线观看| 男女午夜视频在线观看| 动漫黄色视频在线观看| 在线观看免费视频日本深夜| 亚洲人成电影免费在线| 国产色视频综合| 高清毛片免费观看视频网站| 久99久视频精品免费| 国产亚洲欧美98| 午夜久久久在线观看| 黄片播放在线免费| 国产亚洲精品av在线| 手机成人av网站| 黄色 视频免费看| 国产三级在线视频| 在线免费观看的www视频| 90打野战视频偷拍视频| 欧美中文日本在线观看视频| 国产精品国产高清国产av| 97人妻精品一区二区三区麻豆 | 国产在线精品亚洲第一网站| 欧美一区二区精品小视频在线| 一二三四社区在线视频社区8| 午夜福利视频1000在线观看| 黄色毛片三级朝国网站| 精品久久久久久,| 国内毛片毛片毛片毛片毛片| 国产日本99.免费观看| 女人爽到高潮嗷嗷叫在线视频| 女人高潮潮喷娇喘18禁视频| 99国产精品一区二区蜜桃av| 九色国产91popny在线| 国产亚洲精品久久久久久毛片| 国产精品久久久人人做人人爽| 国产精品 国内视频| 高清在线国产一区| 悠悠久久av| 国产一区二区在线av高清观看| 久久香蕉激情| 亚洲专区字幕在线| 久久99热这里只有精品18| 麻豆av在线久日| 国产片内射在线| 午夜两性在线视频| 亚洲精品国产精品久久久不卡| 日本一区二区免费在线视频| 日韩国内少妇激情av| 精华霜和精华液先用哪个| 亚洲精品中文字幕一二三四区| 色av中文字幕| 亚洲人成网站在线播放欧美日韩| 国产精品,欧美在线| 禁无遮挡网站| 午夜日韩欧美国产| 成人三级黄色视频| 俺也久久电影网| 淫妇啪啪啪对白视频| 999久久久精品免费观看国产| 天堂动漫精品| 国产区一区二久久| 一级a爱视频在线免费观看| 91麻豆精品激情在线观看国产| a在线观看视频网站| 黄片播放在线免费| 国产三级在线视频| 狠狠狠狠99中文字幕| 日韩欧美在线二视频| 欧美激情久久久久久爽电影| 十八禁网站免费在线| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看 | 国产黄色小视频在线观看| 特大巨黑吊av在线直播 | 麻豆成人午夜福利视频| 99久久99久久久精品蜜桃| 天天躁狠狠躁夜夜躁狠狠躁| 高清毛片免费观看视频网站| 日韩大尺度精品在线看网址| www.自偷自拍.com| 亚洲黑人精品在线| 亚洲精品av麻豆狂野| 国产高清视频在线播放一区| 一夜夜www| 精品久久久久久成人av| 丝袜在线中文字幕| 18禁美女被吸乳视频| svipshipincom国产片| 麻豆成人av在线观看| 岛国视频午夜一区免费看| 午夜精品久久久久久毛片777| 国产成人精品久久二区二区免费| 国产主播在线观看一区二区| 两性夫妻黄色片| 成人精品一区二区免费| 在线观看一区二区三区| 国产av又大| 中文字幕高清在线视频| 久久国产精品男人的天堂亚洲| 麻豆av在线久日| 啦啦啦观看免费观看视频高清| 大型黄色视频在线免费观看| 亚洲国产欧美一区二区综合| 看片在线看免费视频| av中文乱码字幕在线| av视频在线观看入口| 精品久久久久久久久久久久久 | 这个男人来自地球电影免费观看| 黑人操中国人逼视频| 久久国产精品男人的天堂亚洲| 久久久久国内视频| 免费观看精品视频网站| 97人妻精品一区二区三区麻豆 | 欧美激情久久久久久爽电影| 久久精品影院6| e午夜精品久久久久久久| 757午夜福利合集在线观看| 香蕉av资源在线| 久久久国产成人免费| 中文资源天堂在线| 在线观看午夜福利视频| 亚洲第一欧美日韩一区二区三区| 久久国产精品影院| 国产午夜精品久久久久久| 一进一出抽搐动态| 亚洲欧美日韩高清在线视频| 免费在线观看视频国产中文字幕亚洲| 日韩av在线大香蕉| 黄色片一级片一级黄色片| 欧美成人免费av一区二区三区| 黄色a级毛片大全视频| 脱女人内裤的视频| 亚洲精品国产区一区二| 亚洲专区中文字幕在线| 欧美国产精品va在线观看不卡| 国产一级毛片七仙女欲春2 | 99精品久久久久人妻精品| 亚洲第一电影网av| 夜夜夜夜夜久久久久| 一本综合久久免费| 一区二区三区高清视频在线| 国产精品爽爽va在线观看网站 | 老司机在亚洲福利影院| 亚洲国产看品久久| 黄频高清免费视频| 午夜久久久在线观看| 久热爱精品视频在线9| 亚洲成人精品中文字幕电影| 亚洲久久久国产精品| 男女视频在线观看网站免费 | 男人的好看免费观看在线视频 | 麻豆av在线久日| 夜夜夜夜夜久久久久| 日韩国内少妇激情av| 国产主播在线观看一区二区| 嫁个100分男人电影在线观看| svipshipincom国产片| 757午夜福利合集在线观看| 日本 av在线| 欧美zozozo另类| 国产成人一区二区三区免费视频网站| 亚洲av成人一区二区三| 此物有八面人人有两片| 精品国产乱码久久久久久男人| 免费在线观看成人毛片| 丰满人妻熟妇乱又伦精品不卡| 国产一级毛片七仙女欲春2 | www国产在线视频色| 级片在线观看| 国产高清视频在线播放一区| 亚洲久久久国产精品| 三级毛片av免费| 欧美精品啪啪一区二区三区| 丁香欧美五月| 听说在线观看完整版免费高清| 亚洲黑人精品在线| 国产亚洲欧美精品永久| 看黄色毛片网站| 久久久久久免费高清国产稀缺| 国产精品免费一区二区三区在线| 欧美在线一区亚洲| 婷婷丁香在线五月| 黑人操中国人逼视频| 天天添夜夜摸| 午夜精品在线福利| 夜夜躁狠狠躁天天躁| 亚洲成人久久性| 国产av不卡久久| 欧美日韩福利视频一区二区| 精品高清国产在线一区| 韩国精品一区二区三区| 免费看美女性在线毛片视频| 91麻豆av在线| 搡老熟女国产l中国老女人| 欧美又色又爽又黄视频| 一个人免费在线观看的高清视频| 国产视频内射| 国产成人精品久久二区二区免费| 最好的美女福利视频网| 国产欧美日韩一区二区精品| 侵犯人妻中文字幕一二三四区| 亚洲性夜色夜夜综合| 色尼玛亚洲综合影院| 久久精品人妻少妇| 久久久久久大精品| 好男人在线观看高清免费视频 | 国产精华一区二区三区| 欧美国产日韩亚洲一区| 黄色毛片三级朝国网站| 午夜视频精品福利| 91麻豆精品激情在线观看国产| 久久性视频一级片| 窝窝影院91人妻| 夜夜看夜夜爽夜夜摸| 久久精品国产清高在天天线| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲欧美精品永久| 激情在线观看视频在线高清| 国产免费av片在线观看野外av| 黄片大片在线免费观看| 看免费av毛片| 嫩草影视91久久| 中文字幕另类日韩欧美亚洲嫩草| 欧美一级毛片孕妇| 国产精品1区2区在线观看.| 精品国产美女av久久久久小说| 欧美成人一区二区免费高清观看 | 国产精品av久久久久免费| 国产亚洲欧美精品永久| 老熟妇乱子伦视频在线观看| 一二三四社区在线视频社区8| 一卡2卡三卡四卡精品乱码亚洲| 伊人久久大香线蕉亚洲五| 亚洲国产欧美一区二区综合| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久亚洲av毛片大全| 18禁裸乳无遮挡免费网站照片 | 视频区欧美日本亚洲| 级片在线观看| 欧美黑人精品巨大| 99久久国产精品久久久| 成人av一区二区三区在线看| 久久性视频一级片| 亚洲熟女毛片儿| 成人午夜高清在线视频 | 精品国产乱码久久久久久男人| 精品国产乱子伦一区二区三区| 久久狼人影院| xxxwww97欧美| 天天躁夜夜躁狠狠躁躁| 国产真实乱freesex| 日本成人三级电影网站| a级毛片在线看网站| 老司机深夜福利视频在线观看| 欧美激情极品国产一区二区三区| 脱女人内裤的视频| 国产欧美日韩精品亚洲av| 日本在线视频免费播放| 在线视频色国产色| 午夜老司机福利片| 91字幕亚洲| 在线观看午夜福利视频| 国产精品一区二区精品视频观看| 欧美日韩亚洲综合一区二区三区_| 欧美一级a爱片免费观看看 | 又大又爽又粗| 免费无遮挡裸体视频| 久久久久久大精品| 窝窝影院91人妻| 岛国视频午夜一区免费看| 欧美国产日韩亚洲一区| 男女做爰动态图高潮gif福利片| 免费看日本二区| 久久欧美精品欧美久久欧美| 精品无人区乱码1区二区| 999久久久精品免费观看国产| 1024视频免费在线观看| 亚洲欧美精品综合一区二区三区| a级毛片在线看网站| 女性被躁到高潮视频| 黄色视频,在线免费观看| 精品国产乱子伦一区二区三区| 少妇被粗大的猛进出69影院| 老司机福利观看| 日本五十路高清| 日本一区二区免费在线视频| 男人舔女人的私密视频| 国产成人一区二区三区免费视频网站| 亚洲精品中文字幕在线视频| tocl精华| 成人特级黄色片久久久久久久| 99久久无色码亚洲精品果冻| 欧美不卡视频在线免费观看 | 中出人妻视频一区二区| 成人三级黄色视频| 特大巨黑吊av在线直播 | 九色国产91popny在线| 久久精品国产亚洲av香蕉五月| 久久草成人影院| 男女午夜视频在线观看| 成年免费大片在线观看| 久久人人精品亚洲av| 欧美日韩福利视频一区二区| 亚洲 国产 在线| 麻豆成人午夜福利视频| 黄色 视频免费看| 成年人黄色毛片网站| 亚洲五月色婷婷综合| 国产精品 国内视频| 色婷婷久久久亚洲欧美| av视频在线观看入口| 99国产精品一区二区三区| 香蕉丝袜av| 91九色精品人成在线观看| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 成人亚洲精品一区在线观看| 黑人操中国人逼视频| 亚洲人成77777在线视频| 亚洲中文日韩欧美视频| 精品久久久久久久毛片微露脸| 麻豆一二三区av精品| 十八禁网站免费在线| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 国产片内射在线| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久久久久 | 老司机福利观看| 免费看a级黄色片| e午夜精品久久久久久久| 午夜福利在线在线| 欧美丝袜亚洲另类 | 999久久久精品免费观看国产| 两人在一起打扑克的视频| 国产高清视频在线播放一区| 午夜老司机福利片| 久久久国产欧美日韩av| 亚洲av电影不卡..在线观看| 在线永久观看黄色视频| 精品第一国产精品| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 久久久久久国产a免费观看| 99国产精品99久久久久| 精品久久久久久久毛片微露脸| 精品国产国语对白av| 久久中文字幕人妻熟女| 亚洲一区中文字幕在线| 丝袜人妻中文字幕| 国产一级毛片七仙女欲春2 | 国产不卡一卡二| 欧美另类亚洲清纯唯美| 在线观看一区二区三区| 亚洲国产精品久久男人天堂| 99热只有精品国产| 国产激情欧美一区二区| 久久久精品欧美日韩精品| 久久精品国产综合久久久| 亚洲性夜色夜夜综合| 无限看片的www在线观看| а√天堂www在线а√下载| 成人国产综合亚洲| 久久久久九九精品影院| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 欧美日本亚洲视频在线播放| 亚洲精品国产一区二区精华液| 男女做爰动态图高潮gif福利片| 午夜久久久久精精品| 岛国视频午夜一区免费看| 一级毛片女人18水好多| 19禁男女啪啪无遮挡网站| 欧美最黄视频在线播放免费| 亚洲国产看品久久| 女人被狂操c到高潮| 欧美三级亚洲精品| 久久中文看片网| 午夜老司机福利片| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 精品一区二区三区av网在线观看| 美女大奶头视频| 一级作爱视频免费观看| 国产精品一区二区精品视频观看| 听说在线观看完整版免费高清| 黄片播放在线免费| 日韩精品免费视频一区二区三区| 非洲黑人性xxxx精品又粗又长| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 美女高潮喷水抽搐中文字幕| 校园春色视频在线观看| 日韩欧美国产一区二区入口| 91国产中文字幕| 亚洲色图av天堂| 黄色视频,在线免费观看| 人妻丰满熟妇av一区二区三区| 久久99热这里只有精品18| 精品久久蜜臀av无| 最好的美女福利视频网| 女警被强在线播放| 亚洲 欧美 日韩 在线 免费| 可以在线观看的亚洲视频| 一级a爱片免费观看的视频| e午夜精品久久久久久久| 国产一区在线观看成人免费| 九色国产91popny在线| 啦啦啦 在线观看视频| 欧美日本视频| 久久亚洲真实| 99久久久亚洲精品蜜臀av| 免费在线观看亚洲国产| 午夜a级毛片| 欧美久久黑人一区二区| 色综合站精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | av超薄肉色丝袜交足视频| 精品国产亚洲在线| 久久这里只有精品19| 村上凉子中文字幕在线| 日韩欧美在线二视频| 久久婷婷人人爽人人干人人爱| 亚洲在线自拍视频| 极品教师在线免费播放| 白带黄色成豆腐渣| 久9热在线精品视频| 人人妻人人澡欧美一区二区| 色av中文字幕| 亚洲熟妇中文字幕五十中出| 国产亚洲精品一区二区www| 欧美在线一区亚洲| 麻豆久久精品国产亚洲av| av片东京热男人的天堂| 国产高清激情床上av| 国产在线精品亚洲第一网站| 亚洲无线在线观看| 婷婷亚洲欧美| 1024手机看黄色片| 亚洲成人精品中文字幕电影| 久久午夜综合久久蜜桃| or卡值多少钱| 国产精品一区二区精品视频观看| 亚洲人成电影免费在线| 巨乳人妻的诱惑在线观看| 最新美女视频免费是黄的| 亚洲av日韩精品久久久久久密| 天天一区二区日本电影三级| 两性夫妻黄色片| 99国产极品粉嫩在线观看| 亚洲av电影不卡..在线观看| 长腿黑丝高跟| 亚洲专区字幕在线| 亚洲五月婷婷丁香| 免费看a级黄色片| 丰满人妻熟妇乱又伦精品不卡| 国产黄a三级三级三级人| 老司机深夜福利视频在线观看| 日韩高清综合在线| 日韩欧美国产一区二区入口| svipshipincom国产片| 在线天堂中文资源库| av电影中文网址| 久久久久国内视频| 亚洲全国av大片| 色av中文字幕| 久久久水蜜桃国产精品网| 国产伦在线观看视频一区| 一边摸一边抽搐一进一小说| 亚洲av成人不卡在线观看播放网| 日韩欧美在线二视频| 久久香蕉国产精品| 白带黄色成豆腐渣| 一进一出抽搐动态| 午夜a级毛片| 99精品在免费线老司机午夜| 国产成人精品久久二区二区91| 色哟哟哟哟哟哟| 狂野欧美激情性xxxx| 精品乱码久久久久久99久播| 欧美乱码精品一区二区三区| 免费看十八禁软件| 国产精品久久久久久亚洲av鲁大| 黄色丝袜av网址大全| 久久久久久大精品| 日本成人三级电影网站| 国产日本99.免费观看| 久久精品国产亚洲av香蕉五月| 久久国产亚洲av麻豆专区| 欧美成人性av电影在线观看| 悠悠久久av| 99久久综合精品五月天人人| 亚洲av成人av| 亚洲一区中文字幕在线| 亚洲成人久久爱视频| 午夜a级毛片| 黄色毛片三级朝国网站| 男女视频在线观看网站免费 | 亚洲狠狠婷婷综合久久图片| 一级黄色大片毛片| 精品一区二区三区四区五区乱码| 两个人免费观看高清视频| 可以在线观看毛片的网站| 欧美zozozo另类| 一进一出抽搐gif免费好疼| 一卡2卡三卡四卡精品乱码亚洲| 在线天堂中文资源库| 欧美不卡视频在线免费观看 | 精品卡一卡二卡四卡免费| 国产精品 国内视频| 欧美+亚洲+日韩+国产| 国产成人啪精品午夜网站| 亚洲无线在线观看| 欧美成人午夜精品| 亚洲成a人片在线一区二区| 免费电影在线观看免费观看| 国产精品亚洲美女久久久| 欧美乱妇无乱码| 亚洲国产精品sss在线观看| 国产在线精品亚洲第一网站| 亚洲七黄色美女视频| 国产成人欧美| 亚洲精品国产一区二区精华液| 成人免费观看视频高清| 亚洲五月色婷婷综合| 国产区一区二久久| 亚洲国产欧洲综合997久久, | 后天国语完整版免费观看| 国产一区在线观看成人免费| 黄色毛片三级朝国网站| 18美女黄网站色大片免费观看| 午夜免费成人在线视频| 精品第一国产精品| 欧美成人性av电影在线观看| 精品电影一区二区在线| 日本三级黄在线观看| 国产精品自产拍在线观看55亚洲| 熟女电影av网| 国内精品久久久久精免费| 51午夜福利影视在线观看| 亚洲国产欧美日韩在线播放| 午夜久久久在线观看| 国产一区二区三区在线臀色熟女| 在线十欧美十亚洲十日本专区| 51午夜福利影视在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品999在线| 午夜免费激情av| 日本五十路高清| 亚洲av电影不卡..在线观看| 在线免费观看的www视频| 免费在线观看日本一区| 黑丝袜美女国产一区| 成人永久免费在线观看视频| 欧美日韩乱码在线| 亚洲国产精品合色在线| 成人午夜高清在线视频 | 欧美日本亚洲视频在线播放| 性欧美人与动物交配| 亚洲全国av大片| 久久久国产成人免费| 神马国产精品三级电影在线观看 | 久久午夜综合久久蜜桃| 欧美激情极品国产一区二区三区| 久久精品亚洲精品国产色婷小说| 久久精品91无色码中文字幕| 99久久久亚洲精品蜜臀av| 性欧美人与动物交配| 国产一卡二卡三卡精品| 日日爽夜夜爽网站| а√天堂www在线а√下载| 熟女电影av网| 中文字幕精品亚洲无线码一区 | 亚洲精华国产精华精| 99久久久亚洲精品蜜臀av| 狂野欧美激情性xxxx| 国产91精品成人一区二区三区| 日日爽夜夜爽网站|