• <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.中國科學院沈陽自動化研究所 機器人學國家重點實驗室,遼寧 沈陽 110016; 2.中國科學院 機器人與智能制造創(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é)合許多學者的后續(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 卡門—錢學森公式

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

    (8)

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

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

    (9)

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

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

    (10)

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

    已有相關(guān)學者對于上述經(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ā)展過程中其力學性能參數(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)變特性,通過力學試驗測定鈦合金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ù)學含義只是一個無量綱的歸一化特征向量,即沒有表征缺陷大小的含義,所以需要將陣型向量的最大振幅值放大一定的倍數(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)系

    為充分驗證極限承載力計算的準確性,在中國科學院沈陽自動化研究所水下試驗室開展水壓爆破試驗,將上述經(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)模型
    一半模型
    《形而上學》△卷的結(jié)構(gòu)和位置
    哲學評論(2021年2期)2021-08-22 01:53:34
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(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ù)成長
    国产野战对白在线观看| 亚洲国产欧美日韩在线播放| 丁香六月天网| 精品国产一区二区三区久久久樱花| 欧美亚洲 丝袜 人妻 在线| 亚洲国产av影院在线观看| 啦啦啦在线观看免费高清www| 激情视频va一区二区三区| 一本大道久久a久久精品| 老司机亚洲免费影院| 下体分泌物呈黄色| a 毛片基地| 亚洲五月色婷婷综合| netflix在线观看网站| 我要看黄色一级片免费的| 下体分泌物呈黄色| 国产精品一二三区在线看| 黄色 视频免费看| 女的被弄到高潮叫床怎么办| 精品亚洲乱码少妇综合久久| 国产在线一区二区三区精| 在线观看一区二区三区激情| 日韩 欧美 亚洲 中文字幕| 尾随美女入室| 国产免费现黄频在线看| 51午夜福利影视在线观看| 99热网站在线观看| 精品久久久精品久久久| 一级爰片在线观看| 欧美老熟妇乱子伦牲交| 妹子高潮喷水视频| 丰满迷人的少妇在线观看| 黄片播放在线免费| 国产精品熟女久久久久浪| 一边摸一边抽搐一进一出视频| 观看美女的网站| 久久久久精品国产欧美久久久 | 性高湖久久久久久久久免费观看| 一级片免费观看大全| 国产精品秋霞免费鲁丝片| 两个人看的免费小视频| 99久久综合免费| 日韩一区二区三区影片| 亚洲国产精品国产精品| 秋霞伦理黄片| 人人妻,人人澡人人爽秒播 | 久久久久人妻精品一区果冻| 亚洲国产精品一区三区| 欧美另类一区| 国产高清国产精品国产三级| 在线观看三级黄色| svipshipincom国产片| 多毛熟女@视频| e午夜精品久久久久久久| 人人妻人人澡人人看| 啦啦啦在线观看免费高清www| 天天躁日日躁夜夜躁夜夜| 麻豆精品久久久久久蜜桃| 乱人伦中国视频| 欧美人与善性xxx| a级毛片在线看网站| 天天躁夜夜躁狠狠久久av| 日本av免费视频播放| 亚洲熟女精品中文字幕| 国产精品免费大片| 亚洲,一卡二卡三卡| 一区二区三区乱码不卡18| 丝袜脚勾引网站| 国产精品久久久久久久久免| 国产伦理片在线播放av一区| 成年人免费黄色播放视频| 国产欧美亚洲国产| 人体艺术视频欧美日本| 久久久久精品人妻al黑| 啦啦啦 在线观看视频| 精品国产露脸久久av麻豆| 久久人妻熟女aⅴ| 国产精品偷伦视频观看了| 少妇的丰满在线观看| 国产日韩一区二区三区精品不卡| 久久韩国三级中文字幕| 美国免费a级毛片| 亚洲欧美清纯卡通| 午夜91福利影院| 欧美日韩国产mv在线观看视频| 一区二区三区乱码不卡18| av在线播放精品| 亚洲一区中文字幕在线| 女人被躁到高潮嗷嗷叫费观| 国产又爽黄色视频| 看十八女毛片水多多多| 中文天堂在线官网| 大码成人一级视频| 咕卡用的链子| 亚洲三区欧美一区| 一级a爱视频在线免费观看| 亚洲国产精品成人久久小说| 五月天丁香电影| 精品少妇一区二区三区视频日本电影 | 69精品国产乱码久久久| 亚洲熟女精品中文字幕| 免费av中文字幕在线| 毛片一级片免费看久久久久| 久久午夜综合久久蜜桃| 国产精品欧美亚洲77777| 国产片特级美女逼逼视频| 男人爽女人下面视频在线观看| 丝袜在线中文字幕| 韩国高清视频一区二区三区| 国产熟女午夜一区二区三区| 亚洲av综合色区一区| 丰满迷人的少妇在线观看| 欧美中文综合在线视频| 成人影院久久| 国产视频首页在线观看| www.精华液| 欧美久久黑人一区二区| 久久久久久久精品精品| 亚洲三区欧美一区| 97在线人人人人妻| 日韩一区二区视频免费看| 久久精品国产亚洲av高清一级| 各种免费的搞黄视频| 天天躁日日躁夜夜躁夜夜| 18禁观看日本| 黄色一级大片看看| 亚洲精品,欧美精品| 久久久精品94久久精品| 一区二区三区四区激情视频| a 毛片基地| 美女中出高潮动态图| 尾随美女入室| 看十八女毛片水多多多| 欧美人与性动交α欧美精品济南到| 精品第一国产精品| 热re99久久精品国产66热6| 另类精品久久| 久久精品国产a三级三级三级| 人妻 亚洲 视频| 你懂的网址亚洲精品在线观看| 操美女的视频在线观看| av在线播放精品| 成年av动漫网址| 五月开心婷婷网| 国产精品偷伦视频观看了| 又大又黄又爽视频免费| 精品第一国产精品| 国产 精品1| 亚洲成人免费av在线播放| 国产精品国产av在线观看| 满18在线观看网站| 精品午夜福利在线看| 高清欧美精品videossex| 亚洲精品国产av蜜桃| 国产成人午夜福利电影在线观看| 美女大奶头黄色视频| 欧美日韩福利视频一区二区| 久久久久人妻精品一区果冻| 日本欧美国产在线视频| 亚洲国产欧美日韩在线播放| h视频一区二区三区| 伦理电影大哥的女人| 如何舔出高潮| 国产淫语在线视频| 久久精品人人爽人人爽视色| 欧美日韩成人在线一区二区| 女人精品久久久久毛片| 观看美女的网站| 丝袜脚勾引网站| 十八禁高潮呻吟视频| 亚洲男人天堂网一区| 中文天堂在线官网| 亚洲国产av新网站| 精品第一国产精品| 老司机影院成人| 九九爱精品视频在线观看| 亚洲成人国产一区在线观看 | 精品亚洲乱码少妇综合久久| 一本一本久久a久久精品综合妖精| 51午夜福利影视在线观看| 日韩精品有码人妻一区| av一本久久久久| 亚洲自偷自拍图片 自拍| 蜜桃国产av成人99| 80岁老熟妇乱子伦牲交| 中文字幕人妻丝袜制服| 久久精品久久久久久噜噜老黄| 国产精品欧美亚洲77777| 一本大道久久a久久精品| 亚洲精品国产av蜜桃| kizo精华| 午夜福利网站1000一区二区三区| 亚洲av成人精品一二三区| 在线观看国产h片| 欧美国产精品va在线观看不卡| 一级爰片在线观看| 人人妻人人澡人人爽人人夜夜| 日本av免费视频播放| 国产精品 欧美亚洲| 欧美人与性动交α欧美软件| 韩国精品一区二区三区| 久久av网站| 日韩人妻精品一区2区三区| 国产成人免费观看mmmm| 亚洲 欧美一区二区三区| 亚洲欧美激情在线| 在线观看人妻少妇| 午夜福利在线免费观看网站| 哪个播放器可以免费观看大片| 国精品久久久久久国模美| 国产成人欧美| 亚洲精品成人av观看孕妇| 一区福利在线观看| 欧美日韩综合久久久久久| 亚洲色图 男人天堂 中文字幕| 欧美少妇被猛烈插入视频| av不卡在线播放| 韩国精品一区二区三区| 老汉色av国产亚洲站长工具| 亚洲av欧美aⅴ国产| 久久久久视频综合| 免费日韩欧美在线观看| 欧美精品高潮呻吟av久久| 99re6热这里在线精品视频| 精品国产露脸久久av麻豆| 十八禁人妻一区二区| 一级毛片电影观看| 欧美精品亚洲一区二区| 午夜福利一区二区在线看| a级毛片在线看网站| 国产片内射在线| 亚洲欧美激情在线| 国产伦理片在线播放av一区| 只有这里有精品99| 综合色丁香网| 女人被躁到高潮嗷嗷叫费观| 十八禁网站网址无遮挡| 亚洲天堂av无毛| 日日摸夜夜添夜夜爱| 久久人人97超碰香蕉20202| 国产免费现黄频在线看| 成人三级做爰电影| 午夜精品国产一区二区电影| 一区二区三区精品91| 久久精品亚洲熟妇少妇任你| 亚洲人成77777在线视频| 精品第一国产精品| 赤兔流量卡办理| 操美女的视频在线观看| 嫩草影视91久久| 亚洲精品久久久久久婷婷小说| 婷婷色综合www| 看免费成人av毛片| 我的亚洲天堂| 欧美激情 高清一区二区三区| 国产一区亚洲一区在线观看| 国产亚洲精品第一综合不卡| 国产精品秋霞免费鲁丝片| 建设人人有责人人尽责人人享有的| 中文字幕人妻丝袜制服| av国产久精品久网站免费入址| 国产一卡二卡三卡精品 | 美女视频免费永久观看网站| 国产极品天堂在线| 国产精品免费视频内射| 交换朋友夫妻互换小说| 国产精品熟女久久久久浪| 美女高潮到喷水免费观看| www.熟女人妻精品国产| 在线观看免费视频网站a站| 视频区图区小说| 秋霞在线观看毛片| 国产精品 国内视频| av电影中文网址| 19禁男女啪啪无遮挡网站| av一本久久久久| 国产成人系列免费观看| 七月丁香在线播放| 精品亚洲乱码少妇综合久久| 午夜老司机福利片| 亚洲av综合色区一区| 亚洲国产欧美一区二区综合| 最近手机中文字幕大全| videosex国产| 国产av国产精品国产| 18禁国产床啪视频网站| 欧美国产精品一级二级三级| h视频一区二区三区| 最新在线观看一区二区三区 | 免费黄频网站在线观看国产| 色婷婷av一区二区三区视频| 日韩制服丝袜自拍偷拍| 日韩伦理黄色片| 国产精品一二三区在线看| 欧美精品一区二区大全| 欧美国产精品va在线观看不卡| 老司机影院毛片| 一区二区三区四区激情视频| 精品人妻熟女毛片av久久网站| 男女国产视频网站| 精品人妻在线不人妻| 成年av动漫网址| 丝袜人妻中文字幕| 韩国av在线不卡| 成人毛片60女人毛片免费| 天天躁日日躁夜夜躁夜夜| av有码第一页| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 黑人欧美特级aaaaaa片| 欧美老熟妇乱子伦牲交| av天堂久久9| 一区二区日韩欧美中文字幕| 高清av免费在线| 国产成人系列免费观看| 成年女人毛片免费观看观看9 | 亚洲少妇的诱惑av| 少妇人妻久久综合中文| 黑人巨大精品欧美一区二区蜜桃| 又大又黄又爽视频免费| www.自偷自拍.com| 亚洲av成人精品一二三区| 欧美成人午夜精品| 亚洲一区二区三区欧美精品| 最近2019中文字幕mv第一页| 久久久久久久久久久久大奶| 国产精品99久久99久久久不卡 | 超碰97精品在线观看| 1024视频免费在线观看| 国产精品久久久久久精品电影小说| 免费观看av网站的网址| 老司机深夜福利视频在线观看 | 19禁男女啪啪无遮挡网站| 久久综合国产亚洲精品| 国产成人一区二区在线| 制服人妻中文乱码| 熟女av电影| 搡老乐熟女国产| 一区二区日韩欧美中文字幕| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 亚洲国产精品一区三区| 欧美人与善性xxx| 两个人看的免费小视频| 精品少妇黑人巨大在线播放| 久久人人爽av亚洲精品天堂| 精品少妇黑人巨大在线播放| 丁香六月欧美| 19禁男女啪啪无遮挡网站| 久久99热这里只频精品6学生| 赤兔流量卡办理| √禁漫天堂资源中文www| 精品视频人人做人人爽| 久久久国产一区二区| 日本欧美视频一区| 亚洲欧洲日产国产| 蜜桃在线观看..| av一本久久久久| 赤兔流量卡办理| 免费av中文字幕在线| 午夜福利乱码中文字幕| 国产亚洲最大av| 欧美变态另类bdsm刘玥| 欧美精品av麻豆av| 黄片小视频在线播放| 在线观看国产h片| 亚洲精华国产精华液的使用体验| 黄片播放在线免费| 久久鲁丝午夜福利片| av在线播放精品| 黄色视频不卡| 91精品伊人久久大香线蕉| 日韩熟女老妇一区二区性免费视频| 日日摸夜夜添夜夜爱| 亚洲中文av在线| 91精品国产国语对白视频| 欧美日本中文国产一区发布| 中文字幕高清在线视频| 国产熟女欧美一区二区| 不卡av一区二区三区| 亚洲欧美一区二区三区黑人| 亚洲一区中文字幕在线| 午夜日本视频在线| 亚洲av综合色区一区| 日本av免费视频播放| 美女视频免费永久观看网站| 少妇人妻 视频| 老熟女久久久| 狂野欧美激情性xxxx| 免费少妇av软件| 亚洲国产欧美在线一区| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 99re6热这里在线精品视频| 久久99热这里只频精品6学生| 老汉色av国产亚洲站长工具| 啦啦啦视频在线资源免费观看| 久久久久久久久久久免费av| 最近手机中文字幕大全| 亚洲av欧美aⅴ国产| 亚洲欧美精品综合一区二区三区| 亚洲人成电影观看| 999精品在线视频| 国产 一区精品| 下体分泌物呈黄色| 天天添夜夜摸| 爱豆传媒免费全集在线观看| 国产精品一区二区精品视频观看| 国产又爽黄色视频| 美女主播在线视频| 久久精品亚洲熟妇少妇任你| 另类亚洲欧美激情| 在线观看免费午夜福利视频| av女优亚洲男人天堂| 黑丝袜美女国产一区| 国产精品欧美亚洲77777| 久久人人97超碰香蕉20202| 免费观看av网站的网址| 国产av国产精品国产| 另类亚洲欧美激情| 狂野欧美激情性xxxx| 18禁观看日本| 制服诱惑二区| 97精品久久久久久久久久精品| 天堂中文最新版在线下载| 久久韩国三级中文字幕| 国产精品二区激情视频| 欧美少妇被猛烈插入视频| 国产精品麻豆人妻色哟哟久久| 观看av在线不卡| 亚洲欧美激情在线| 国产成人精品久久二区二区91 | 无遮挡黄片免费观看| 久久精品久久久久久噜噜老黄| 丰满饥渴人妻一区二区三| 亚洲精品视频女| 性色av一级| 哪个播放器可以免费观看大片| 免费黄网站久久成人精品| 操美女的视频在线观看| h视频一区二区三区| 精品第一国产精品| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 中文天堂在线官网| 国产无遮挡羞羞视频在线观看| 久久久久精品久久久久真实原创| 天天躁狠狠躁夜夜躁狠狠躁| 国产视频首页在线观看| 美女高潮到喷水免费观看| 国语对白做爰xxxⅹ性视频网站| 国产女主播在线喷水免费视频网站| 国产亚洲av高清不卡| 一区在线观看完整版| 天天躁夜夜躁狠狠躁躁| 亚洲国产欧美在线一区| kizo精华| 男男h啪啪无遮挡| 成年美女黄网站色视频大全免费| 国产精品免费视频内射| 最黄视频免费看| 国产亚洲精品第一综合不卡| 丰满乱子伦码专区| 久久鲁丝午夜福利片| 亚洲综合精品二区| 麻豆av在线久日| 悠悠久久av| 十分钟在线观看高清视频www| 中文天堂在线官网| 国产黄色视频一区二区在线观看| 哪个播放器可以免费观看大片| 亚洲国产精品国产精品| 国产激情久久老熟女| 亚洲国产精品999| 国产99久久九九免费精品| 国产精品av久久久久免费| 日韩欧美精品免费久久| 女的被弄到高潮叫床怎么办| 一本色道久久久久久精品综合| 亚洲精品久久成人aⅴ小说| 精品一区二区三区av网在线观看 | 青春草亚洲视频在线观看| 在线天堂中文资源库| av电影中文网址| 亚洲精品视频女| 纯流量卡能插随身wifi吗| 18禁裸乳无遮挡动漫免费视频| 赤兔流量卡办理| 亚洲精品av麻豆狂野| 国产伦人伦偷精品视频| 国产乱来视频区| 精品国产一区二区久久| 秋霞伦理黄片| 最近中文字幕高清免费大全6| 在线天堂最新版资源| 免费黄色在线免费观看| 一区二区av电影网| 日韩欧美一区视频在线观看| 久久青草综合色| 成人亚洲精品一区在线观看| 午夜福利视频精品| 亚洲成色77777| 午夜免费鲁丝| 亚洲男人天堂网一区| 黄色怎么调成土黄色| 午夜福利,免费看| 免费久久久久久久精品成人欧美视频| 日韩电影二区| 在线观看国产h片| 99国产精品免费福利视频| 亚洲国产精品一区三区| 亚洲国产精品999| 国产精品一区二区在线不卡| 精品人妻在线不人妻| 丰满乱子伦码专区| 蜜桃在线观看..| 国产人伦9x9x在线观看| 两个人免费观看高清视频| 免费久久久久久久精品成人欧美视频| 搡老岳熟女国产| 色综合欧美亚洲国产小说| 日韩 欧美 亚洲 中文字幕| 一级片'在线观看视频| 日韩视频在线欧美| av网站在线播放免费| 一区福利在线观看| 中文字幕最新亚洲高清| 国产伦理片在线播放av一区| 亚洲情色 制服丝袜| 中文字幕人妻丝袜一区二区 | 午夜福利,免费看| xxxhd国产人妻xxx| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 波野结衣二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 大陆偷拍与自拍| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 久久精品国产综合久久久| 国产男女超爽视频在线观看| 一本久久精品| 亚洲一码二码三码区别大吗| 国产 精品1| 亚洲精品视频女| 亚洲专区中文字幕在线 | 亚洲国产日韩一区二区| 国产精品久久久久久精品电影小说| 99香蕉大伊视频| 美女中出高潮动态图| 亚洲精品av麻豆狂野| 欧美精品人与动牲交sv欧美| 免费黄网站久久成人精品| 久久久精品94久久精品| 国产男女超爽视频在线观看| 麻豆精品久久久久久蜜桃| 最新在线观看一区二区三区 | 国产精品av久久久久免费| 建设人人有责人人尽责人人享有的| 最近最新中文字幕免费大全7| 国产伦理片在线播放av一区| 2021少妇久久久久久久久久久| 波多野结衣一区麻豆| 欧美97在线视频| 国产免费视频播放在线视频| 亚洲国产欧美一区二区综合| 高清在线视频一区二区三区| 成人18禁高潮啪啪吃奶动态图| 免费观看人在逋| 999精品在线视频| 免费久久久久久久精品成人欧美视频| 国产免费现黄频在线看| 麻豆乱淫一区二区| 一区二区av电影网| 亚洲,欧美精品.| 婷婷色综合www| av免费观看日本| 在线观看免费视频网站a站| 99久久综合免费| 只有这里有精品99| 日本欧美视频一区| 久久韩国三级中文字幕| 18在线观看网站| 亚洲欧美色中文字幕在线| 中文字幕av电影在线播放| 午夜日本视频在线| 亚洲国产毛片av蜜桃av| 日日爽夜夜爽网站| 这个男人来自地球电影免费观看 | 在线观看国产h片| 精品免费久久久久久久清纯 | 大香蕉久久网| 亚洲精品成人av观看孕妇| 少妇的丰满在线观看| 国产一区二区激情短视频 | 国产精品熟女久久久久浪| 免费在线观看黄色视频的| 一区二区日韩欧美中文字幕| 99九九在线精品视频| 亚洲色图 男人天堂 中文字幕| 亚洲精品久久成人aⅴ小说| 少妇猛男粗大的猛烈进出视频| 男女免费视频国产| 韩国av在线不卡| 观看av在线不卡| 国产精品嫩草影院av在线观看| 赤兔流量卡办理| bbb黄色大片| 国产日韩欧美在线精品| 欧美日韩亚洲国产一区二区在线观看 | tube8黄色片| 亚洲国产精品一区三区| av有码第一页| 欧美人与性动交α欧美软件| 亚洲欧美精品自产自拍|