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

    船-冰碰撞下冰材料有限元數(shù)值方法研究進展

    2021-09-02 03:00:06畢璐澤
    船舶力學(xué) 2021年8期
    關(guān)鍵詞:冰體計算結(jié)果數(shù)值

    蔡 偉,朱 凌,3,畢璐澤

    (1.高性能船舶技術(shù)教育部重點實驗室(武漢理工大學(xué)),武漢 430063;2.武漢理工大學(xué)交通學(xué)院,武漢 430063;3.高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,武漢 430063)

    0 引 言

    極地資源的開發(fā)、航道運輸以及科考成為當今世界各國關(guān)注的熱點話題。隨著全球氣候變暖,北極冰川加速融化,極地生產(chǎn)運輸活動逐漸增多,致使浮冰,特別是冰山逐漸增多,這給航行于冰區(qū)的船舶的結(jié)構(gòu)安全帶來了巨大的挑戰(zhàn)。極地地區(qū)海面上漂浮的冰體主要分為淡水冰和咸水冰,其中淡水冰主要是冰山冰,或者是一些從冰川分離而出的小型冰山。極地船舶在冰區(qū)航行時難免會遭受到浮冰的碰撞作用,在極端冰碰載荷作用下船體結(jié)構(gòu)難免會出現(xiàn)塑性變形,甚至?xí)霈F(xiàn)破壞斷裂,如圖1和圖2所示,這將會導(dǎo)致人員傷亡、貨物泄露以及環(huán)境污染等問題,甚至?xí)斐沙链鹿?。例如?015年,俄羅斯大型拖網(wǎng)漁船“遠東”號與浮冰相撞,導(dǎo)致船舶沉沒,56人死亡以及13人失蹤;2019年1月,“雪龍”號破冰船在南極阿蒙森海與冰山相撞,船艏桅桿被撞倒,部分舷墻受損,對考察任務(wù)造成了嚴重影響。因此,亟需開展冰碰載荷下船舶結(jié)構(gòu)動力響應(yīng)評估方法研究,以探討冰體碰撞下結(jié)構(gòu)的塑性變形損傷。

    船舶與海洋工程結(jié)構(gòu)與浮冰碰撞是結(jié)構(gòu)在短時間內(nèi)遭受巨大沖擊載荷作用下的一種復(fù)雜非線性動態(tài)響應(yīng)過程,存在著材料非線性、幾何非線性、接觸非線性和運動非線性等問題。針對船-冰碰撞結(jié)構(gòu)動力響應(yīng)問題,目前尚未有較為成熟全面的理論方法,尤其是在冰材料力學(xué)模型方面,目前國內(nèi)外學(xué)者主要是采用經(jīng)驗公式法、試驗方法和數(shù)值方法來研究冰體碰撞下結(jié)構(gòu)動力響應(yīng)問題。隨著計算機技術(shù)的不斷提高,計算力學(xué)方法在不斷發(fā)展,數(shù)值模擬方法已成為一種最為經(jīng)濟高效的方法[1],主要包括有限元方法(FEM)、離散元方法(DEM)[2]、光滑粒子動力學(xué)方法(SPH)[3]以及近場動力學(xué)方法(PD)[4]等。有限元法在模擬介質(zhì)間的間斷和離散性方面能力有限,相比之下,離散元方法在模擬冰體離散特性方面具有一定的優(yōu)勢[5]。除此之外,有限元方法在模擬冰體斷裂以及裂紋擴展方面有一定的局限性,而光滑粒子動力學(xué)方法和近場動力學(xué)方法的粒子離散形式使其在海冰破碎以及裂紋擴展方面具有更突出的優(yōu)勢[4]。然而,在船舶結(jié)構(gòu)動力響應(yīng)研究方面,尤其是針對結(jié)構(gòu)塑性大變形以及斷裂失效問題,有限元方法比離散元、光滑粒子動力學(xué)以及近場動力學(xué)方法顯得更有優(yōu)勢,更加成熟。

    由于冰體是一種力學(xué)性質(zhì)較為復(fù)雜的天然復(fù)合材料,其材料力學(xué)特性的影響因素眾多,尤其是對于海冰材料,主要受鹽分、雜質(zhì)等成分含量、溫度、加載速率和形成年份等影響[6]。船-冰碰撞過程是一種瞬態(tài)動力響應(yīng)問題,由于冰材料的固有特性,冰隨著應(yīng)變率的增加會由韌性破壞到脆性破壞[7-8]。由于冰體材料的復(fù)雜力學(xué)特性,不同條件下冰的力學(xué)性質(zhì)差異較大,其材料數(shù)值模型仍然處于研究探索階段。目前,國內(nèi)外學(xué)者提出冰體理想化材料模型,主要為各向同性彈性失效模型[9-14]、各向同性彈塑性失效應(yīng)變率相關(guān)材料模型[15-18]、可壓碎泡沫型材料模型[19-25]、損傷材料模型[26-27]以及一些其它自定義彈塑性材料數(shù)值模型[28-35]等。在船-冰碰撞過程中,冰體的破壞包括了擠壓、彎曲破壞、堆積等動態(tài)過程,其破壞力學(xué)特性不僅與應(yīng)變率、溫度等因素有關(guān),還與其自身形成條件有密切關(guān)系,如年份、鹽分、雜質(zhì)等成分含量,利用數(shù)值模型方法來研究冰體破壞失效過程具有較大的挑戰(zhàn)性。除此之外,通過一些冰體壓縮、彎曲等試驗測得冰體材料參數(shù),如楊氏模量、屈服應(yīng)力、剪切模量和泊松比等,這些結(jié)果具有較大的離散性,開發(fā)一種涵蓋所有條件下的冰體材料模型是不現(xiàn)實的。因此,目前尚未有良好的數(shù)值模型來描述所有條件下的冰體材料力學(xué)特性,上述所提出的材料模型都有各自的局限性。

    因此,本文主要總結(jié)了現(xiàn)有幾種常見的冰體有限元數(shù)值材料模型,并與現(xiàn)有的一些試驗結(jié)果、ISO 推薦的壓力-面積曲線以及單位體積能量吸收經(jīng)驗值進行對比分析,討論現(xiàn)有一些冰體材料模型的適用性和優(yōu)缺點。

    1 冰材料有限元數(shù)值模型研究進展

    有限元數(shù)值計算方法是船-冰碰撞下結(jié)構(gòu)動力響應(yīng)的重要評估手段,但是碰撞過程中冰材料的數(shù)值模擬是其最為關(guān)鍵的難點問題,下面將介紹幾種常見的冰材料數(shù)值模型。

    在大多數(shù)船-冰碰撞情況過程中,船體與冰在接觸過程中高應(yīng)變率使冰體較明顯地體現(xiàn)出脆性性質(zhì),故國內(nèi)外許多學(xué)者將冰體材料模型簡化為各向同性彈性失效材料模型。最早在針對冰雹與復(fù)合材料結(jié)構(gòu)碰撞響應(yīng)進行研究時,Kim 等[9-10]以及Anghileri 等[11]利用ANSYS/LS-DYNA 軟件中的各向同性彈性斷裂模型來模擬在碰撞過程中冰體的力學(xué)特性。該模型采用Von-mises 失效準則,將最大塑性應(yīng)變定義為其破壞模式,把恒定最小壓力模式定義為其分離模式,來模擬冰體失效行為。雖然冰體材料模型比較簡單,可以節(jié)省計算時間,但是僅僅依據(jù)最大塑性應(yīng)變和截斷壓力兩個恒定值來判斷冰體失效模式有所不足。結(jié)合相關(guān)冰體力學(xué)參數(shù),一些學(xué)者[12-14]逐漸將此模型運用于船-冰碰撞分析之中??紤]到冰體材料力學(xué)特性受應(yīng)變率影響較大,Carney 等[15]采用一種各向同性彈塑性失效應(yīng)變率相關(guān)材料模型來模擬冰體失效,此模型廣泛運用于冰雹的高速沖擊數(shù)值模擬[16-17]。胡志寬[18]將各向同性彈塑性失效應(yīng)變率相關(guān)材料模型運用于冰與螺旋槳的碰撞數(shù)值模擬研究中。此種材料模型雖然考慮到了應(yīng)變率變化對冰體力學(xué)特性的影響,但是在船-冰碰撞結(jié)構(gòu)動態(tài)響應(yīng)研究中,碰撞速度有限,其應(yīng)變率變化范圍并不明顯,沒有冰雹高速沖擊的應(yīng)變率大,故該材料模型的適用性具有一定的局限性。

    由于冰體材料參數(shù)試驗測量數(shù)據(jù)離散性較大,具體的數(shù)值大小難以確定,一種可壓縮泡沫型材料模型[19-21]被廣泛運用于船-冰碰撞研究之中,通過運用自定義應(yīng)力-體積應(yīng)變關(guān)系來定義可壓碎泡沫模型的力學(xué)行為??紤]到在碰撞過程中冰體與結(jié)構(gòu)物接觸區(qū)域存在著高壓區(qū),在其接觸區(qū)域周圍存在著低壓區(qū),如圖3所示,但是上述的可壓縮泡沫型材料模型還不能體現(xiàn)出這一區(qū)別。為此,Gagnon[22]對此材料模型所定義應(yīng)力-體積應(yīng)變關(guān)系進行了改進,通過對冰體模型進行分層,在每層賦予不同材料特性的鋸齒形的應(yīng)力-體積應(yīng)變關(guān)系曲線,如圖4所示,來對冰體的高壓區(qū)域和低壓區(qū)域進行模擬。類似于Gagnon 的分層方法,Kim 等[23-25]通過對冰體壓縮試驗現(xiàn)象與結(jié)果的觀察和分析,定義了兩條鋸齒形的應(yīng)力-體積應(yīng)變關(guān)系曲線,并賦予冰體不同區(qū)域處,以表達冰體與結(jié)構(gòu)物接觸時高壓區(qū)域和低壓區(qū)域。此外,Kim 還運用“最大主應(yīng)力”失效準則用來定義冰體破壞,即當冰體網(wǎng)格單元大于最大失效應(yīng)力時,單元自動失效刪除,采用此種材料模型的數(shù)值計算結(jié)果與試驗結(jié)果吻合較好。相比于其它冰材料模型,可壓碎泡沫材料模型在描述冰體力學(xué)特性時在一定程度上有較高的準確性,它可以宏觀地體現(xiàn)出冰體與結(jié)構(gòu)物相互接觸時冰體的變化,可以比較真實地反映實驗過程中的冰體與結(jié)構(gòu)物的接觸力大小,但是其缺乏強有力的物理解釋,無法模擬出冰體脆性失效特性以及冰體在碰撞過程中出現(xiàn)的內(nèi)部缺陷,該模型只反映特定條件下的冰體材料力學(xué)行為。

    冰體在碰撞過程中主要體現(xiàn)出脆性失效性質(zhì),考慮到巖土、混凝土等材料本構(gòu)模型與冰體材料在某種程度上有一定的相似性,而混凝土材料模型研究相對于冰體材料數(shù)值模型研究較為成熟,一些學(xué)者逐漸將一些成熟的混凝土損傷材料模型用于冰體材料模型數(shù)值模擬當中。Polach 和Ehlers[26]使用ANSYS/LS-DYNA 軟件中損傷材料模型(*MAT_DAMAGE_3)來模擬冰體材料,使用Lemaitre 損傷模型來模擬冰體失效,當冰體損傷值大于臨界損傷值時,冰材料就會失效。該模型的數(shù)值計算結(jié)果與實驗結(jié)果吻合較好[37]。此外,Cai 等[27]基于混凝土材料本構(gòu)模型,結(jié)合冰體材料力學(xué)模型,建立了表述應(yīng)力-應(yīng)變關(guān)系的本構(gòu)方程的屈服面公式,以及失效應(yīng)變與截斷壓力之間關(guān)系等式,提出了基于混泥土材料本構(gòu)模型的冰體數(shù)值材料模型,數(shù)值計算結(jié)果與模型試驗結(jié)果吻合較好。由于冰體材料特性的影響因素眾多,如溫度、應(yīng)變率和鹽度,開發(fā)出一種涵蓋所有條件下的冰體材料模型是不現(xiàn)實的。相對于海冰來說,冰山冰為淡水冰,材料屬性較為簡單,可將其視作各向同性材料[38]。Liu 等[28-29]通過開發(fā)LS-DYNA 軟件中的自定義各向同性的彈塑性材料模型來模擬冰山的力學(xué)行為,在所開發(fā)的彈塑性模型中,使用屈服面公式和流動法則來表征冰材料的應(yīng)力-應(yīng)變關(guān)系,建立失效應(yīng)變與壓力之間的關(guān)系曲線來定義冰體的失效準則?;贚iu 等人的自定義冰體彈塑性材料模型,Gao 等[31]提出了一種新的失效準則來模擬冰體失效,Song 等[33]利用Liu 的冰材料模型開展了考慮附連水質(zhì)量的船-冰碰撞有限元計算研究。此外,Shi 等[34]提出了一種與溫度梯度相關(guān)的冰體彈塑性材料模型,用于研究溫度對船舶-冰山碰撞沖擊力的影響。上述的自定義材料模型雖然在冰體失效準則和冰體三向應(yīng)力狀態(tài)模擬上更加準確,但是依然很難反映冰體真實裂紋延伸和斷裂形態(tài),相比于可壓縮材料模型,無法準確模擬計算出冰體在碰撞過程中接觸面壓力大小分布特征。除此之外,Ince等[35]基于KOSORI冰模型的本構(gòu)關(guān)系[36],考慮應(yīng)變速率、溫度和鹽度等參數(shù)的影響,對Johnson-Cook 材料模型進行了修正,通過ABAQUS軟件用戶自定義材料模型(UMAT)來實現(xiàn)冰材料模型的數(shù)值模擬。

    在自然界中,無論是海冰還是淡水冰,不同環(huán)境下冰體的形成機理各不相同[6,36,38-39],如在極地地區(qū)有1年冰和多年冰,不同年份的冰層的材料力學(xué)特性相差較大。冰的應(yīng)力-應(yīng)變特性依賴于應(yīng)變率和溫度,與冰的形成成分如鹽分等也有直接關(guān)系,不同條件下測得的冰體材料參數(shù)也各不相同,如冰體密度、屈服應(yīng)力、楊氏模量、剪切強度、拉伸強度、壓縮強度以及斷裂韌度等參數(shù)。開發(fā)一種涵蓋所有條件下的冰體材料模型是難以實現(xiàn)的,每種冰體數(shù)值材料模型都具有其局限性。

    2 模型試驗與數(shù)值模擬結(jié)果對比分析

    船-冰相互作用試驗研究方法主要包括壓縮和拉伸試驗、彎曲試驗、貫穿與剪切實驗、壓痕試驗以及碰撞試驗[5]。其中,壓縮與拉伸試驗、彎曲試驗、貫穿與剪切實驗以及壓痕試驗主要用于冰體力學(xué)特性研究。在船-冰相互作用過程中,有些工況是低速大質(zhì)量碰撞情形,通常采用準靜態(tài)加載實驗方法來研究結(jié)構(gòu)的彈塑性動力響應(yīng);對于碰撞速度較大時,便采用動態(tài)沖擊實驗方法來研究結(jié)構(gòu)的彈塑性動力響應(yīng)。下面針對三種模型試驗工況來進行數(shù)值模擬分析,并對不同冰材料模型的計算結(jié)果進行分析。值得說明的是,由于所分析的模型實驗工況有各自的特點,所以對下面所選取的冰材料模型參數(shù)與所參考的原文中的冰材料模型進行了微小的改動,以提高數(shù)值計算結(jié)果與模型實驗結(jié)果的可比性。

    2.1 冰體壓縮實驗數(shù)值模擬

    Kim 等[24]開展了不同加載速率圓錐冰體壓縮實驗,實驗裝置如圖5 所示,還采用可壓縮泡沫型冰體材料模型來進行相應(yīng)的數(shù)值仿真研究。在數(shù)值模擬中通過定義兩條鋸齒形的應(yīng)力-體積應(yīng)變關(guān)系曲線,并賦予冰體不同區(qū)域處,采用“最大主應(yīng)力”失效準則來定義冰體破壞失效準則,數(shù)值結(jié)果與實驗結(jié)果吻合較好。將Kim 等人所開展的冰體壓縮實驗其中一組實驗工況來作為有限元數(shù)值模擬對象,其主要實驗參數(shù)如表1所示,有限元數(shù)值模型如圖6所示。圖7對比了模型實驗和6組數(shù)值模型計算結(jié)果,可以發(fā)現(xiàn):各向同性彈性失效模型(MAT_13)和混凝土損傷模型(MAT_78)的計算結(jié)果會出現(xiàn)多個峰值,并出現(xiàn)多次加卸載過程。這是由于這兩種材料模型在加載過程中冰模型會出現(xiàn)網(wǎng)格失效和刪除,導(dǎo)致壓塊與冰體存在多次接觸的過程。因此,在模擬冰體壓縮實驗過程中,采用網(wǎng)格單元失效刪除的方式來模擬冰體失效是不準確的,很難對冰體壓縮力-位移變化機制進行預(yù)測,會導(dǎo)致最終所預(yù)測的載荷值偏小。

    表1 冰體壓縮實驗和數(shù)值模擬主要參數(shù)Tab.1 Main parameters of ice compression test and numerical simulation

    圖7 圓錐冰壓縮實驗與數(shù)值的力-位移曲線結(jié)果對比Fig.7 Comparison of experimental and numerical force-displacement curves for conical ice compression

    2.2 冰體水平碰撞實驗數(shù)值模擬

    Cai 等[27]利用水平?jīng)_擊試驗機開展了楔形冰與船體板模型碰撞實驗,實驗裝置如圖8 所示。本文選取的碰撞實驗工況主要參數(shù)如表2 所示,對應(yīng)的有限元數(shù)值模型如圖9 所示,冰體前端長度為200 mm,前端夾角為90°。下面選取了5種材料模型來進行對比,模型實驗與有限元數(shù)值計算結(jié)果對比如圖10所示。圖10對比了各向同性彈性失效模型(MAT_13)、可壓碎泡沫型材料模型(MAT_63)以及損傷模型(|MAT_78,MAT_153),從對比結(jié)果可以發(fā)現(xiàn):采用低壓區(qū)的低應(yīng)力-體積應(yīng)變曲線(M2)的可壓碎泡沫型材料模型的計算結(jié)果與實驗結(jié)果相比誤差比較大;在碰撞力預(yù)測方面,采用各向同性彈性失效模型和損傷材料模型(MAT_78)的計算結(jié)果與實驗結(jié)果比較接近;在船體板結(jié)構(gòu)變形預(yù)測方面,采用損傷材料中的MAT_78 和MAT_153 單元材料模型的計算結(jié)果與實驗結(jié)果比較接近。在船-冰碰撞數(shù)值模擬過程中,由于可壓碎泡沫型材料模型很難模擬出冰體破壞失效特性,在模擬冰體破壞幾何特征方面存在很大的不足之處。相比較而言,船-冰碰撞過程中冰體與船體局部結(jié)構(gòu)相互作用特性研究方面,由于冰體在碰撞擠壓過程中會出現(xiàn)以顆粒形式剝落的現(xiàn)象,采用如各向同性彈性失效模型(MAT_13)和損傷材料模型(MAT_78、MAT_153)等具有單元失效刪除功能的冰材料模型來進行數(shù)值模擬具有較大的優(yōu)勢。

    表2 冰體水平碰撞實驗主要參數(shù)Tab.2 Main parameters of ice horizontal impact test

    圖10 楔形冰碰撞實驗和數(shù)值的力-位移曲線結(jié)果對比[27]Fig.10 Comparison of experimental and numerical force-displacement curves for wedge ice collision

    2.3 冰體落體碰撞實驗數(shù)值模擬

    Ince等[35]開展了圓錐冰體落體碰撞實驗,實驗裝置如圖11所示,此外還進行相對應(yīng)的有限元數(shù)值模擬,基于KOSORI冰模型本構(gòu)關(guān)系,采用了改進的Johnson-cook材料模型來模擬冰體力學(xué)行為,數(shù)值計算結(jié)果與實驗結(jié)果吻合較好。本文選取的碰撞實驗工況主要參數(shù)如表3 所示,對應(yīng)的有限元數(shù)值模型如圖12 所示。從圖13 的對比結(jié)果可以發(fā)現(xiàn):不同的冰材料模型數(shù)值計算結(jié)果差距較大,在板的變形為60 mm之前時,MAT_63(M2)材料單元模型、MAT_78單元模型與KOSORI冰模型的計算結(jié)果比較接近,與實驗結(jié)果吻合較好:當板的變形超過60 mm時,上述的冰材料模型計算結(jié)果差異較大。

    表3 冰體落體碰撞實驗主要參數(shù)Tab.3 Main parameters of ice drop impact test

    圖13 圓錐冰落體碰撞實驗和數(shù)值模擬的力-位移曲線結(jié)果對比Fig.13 Comparison of the experimental and numerical force-displacement curves for conical ice drop impact

    3 冰材料數(shù)值模型可行性評估方法

    在船-冰碰撞下結(jié)構(gòu)動力響應(yīng)數(shù)值模擬中,冰材料數(shù)值模型的準確性與可行性評估是船-冰碰撞數(shù)值研究中的一個重要環(huán)節(jié)。目前國內(nèi)外學(xué)者主要采用的評估方法有模型實驗驗證方法[20,24,26-27,35]、ISO 規(guī)范等推薦的壓力面積曲線(P-A)對比方法[13,27-28,31,40]以及比吸能(SEA,Specific Energy Absorption)判定方法[41-43]。模型實驗方法是評估數(shù)值模擬方法可行性的最好方法之一,可以通過實驗直接觀察和分析冰體失效過程,并獲得直觀的實驗數(shù)據(jù),但是模型實驗成本較高,耗時較長,且其提供的驗證數(shù)據(jù)有實驗條件限制,具有較大的不確定性。為此,許多學(xué)者基于大量實驗測量得到的數(shù)據(jù),對冰體碰撞過程中接觸壓力和接觸面積進行分析,通過擬合得到壓力-面積關(guān)系,從而為數(shù)值模型驗證提供參考標準。比如,Masterson 等[44]總結(jié)了一系列的試驗數(shù)據(jù),提出了一個由許多研究人員根據(jù)現(xiàn)場測量數(shù)據(jù)擬合得到的P-A關(guān)系式,并將其納入ISO 規(guī)范[45]。Cai等[27]開展了冰體沖擊剛性墻的有限元數(shù)值計算,參考一些現(xiàn)有的有限元計算結(jié)果,總結(jié)了不同冰體材料模型下壓力-面積值,并與ISO/CD 19906[43]以及API/CSA[46]推薦的壓力-面積曲線進行對比,如圖14 所示。利用現(xiàn)有的一些基于實驗數(shù)據(jù)得到的P-A曲線,來對所提出的冰材料數(shù)值計算模型進行標定和校對,是目前比較常用的方法。

    圖14 不同冰體材料模型P-A值與ISO規(guī)范推薦曲線對比[27]Fig.14 Comparison between P-A values calculated by different ice material models and P-A curve recommended by ISO rule[27]

    對于一種破壞失效吸能材料,單位的質(zhì)量材料破壞失效所吸收的能量值(SEA)常用來衡量材料的吸能特性。船-冰碰撞過程中冰體體現(xiàn)出一種較強的擠壓破壞失效特性,因此,SEA也用來表征冰體材料的吸能特性。圖15展示了不同縮尺比實驗測定得到的SEA 值,根據(jù)冰體破壞體積Vc,主要劃分為三個區(qū)域:區(qū)域1,當1.0E-8 m3≤Vc≤1.0E-7 m3時,SEA值的范圍大致在1.0E4 J/kg~1.0E5J/kg之間;區(qū)域2,當1.0E-7 m3≤Vc≤1.0E-5m3時,SEA值的范圍大致在1.0E3 J/kg~1.0E4J/kg之間;區(qū)域3,當1.0E-5m3≤Vc≤1.0E-3m3時,SEA值主要集中在1.0E4 J/kg 左右。在各自區(qū)域內(nèi),隨著冰體破碎體積的增大,SEA 值在逐漸減小。由于圖15 中總結(jié)的實驗數(shù)據(jù)所對應(yīng)的模型實驗冰體尺度不同,所以展現(xiàn)出了SEA 值的三個劃分區(qū)域[42]。根據(jù)Cai等[27]文中的楔形冰與板模型碰撞實驗工況,主要碰撞參數(shù)如表2 工況1 所示,開展不同冰材料模型下有限元數(shù)值計算,得到不同破碎體積下的SEA 值,如圖16所示。由于受到模型實驗尺寸的限制,隨著碰撞時間的變化,冰體的破壞體積范圍主要集中在區(qū)域3。圖16 給出了不同冰材料模型下不同冰體破壞體積所對應(yīng)的SEA數(shù)值計算值,可以發(fā)現(xiàn):圖中的冰材料模型所計算得到的SEA值都小于實驗參考平均值10 000 J/kg,各向同性彈性失效模型(MAT_13)[13]和混凝土損傷模型(MAT_78)[27]計算得到的SEA 值隨著破壞體積的增大有下降的趨勢,而模擬高壓區(qū)的可壓碎泡沫型材料模型(MAT_63,M1)[22]計算得到的SEA 值隨著破壞體積的增大有著明顯的上升趨勢,對于模擬低壓區(qū)的可壓碎泡沫型材料 模 型(MAT_63,M2)[22]和 損 傷 模 型(MAT_153)[26]計算得到的SEA 值沒有很明顯的變化趨勢,相對比較穩(wěn)定。除此之外,從表4 可以看出,各向同性彈性失效模型(MAT_13)和損傷模型(MAT_153)所計算得到的SEA 平均值相對于其它模型的SEA 平均值和實驗參考值較小。每種材料模型所定義的材料參數(shù)都有著各自對應(yīng)的適用條件,與圖15 中的實驗所對應(yīng)的條件有所區(qū)別,因此模型實驗的SEA值和數(shù)值模擬的SEA值有所差別。

    圖15 不同縮尺比實驗測定的SEA值[42]Fig.15 SEA values measured by model tests with different scaled ratios[42]

    圖16 不同冰材料模型下不同冰體破壞體積所對應(yīng)的SEA數(shù)值計算值Fig.16 Numerical calculation values of SEA under different ice material models versus different ice crushed volumes

    表4 不同冰材料模型下的SEA數(shù)值計算值Tab.4 Numerical calculation values of SEA under different ice material models(Experimental reference value:10 000 J/kg)

    4 結(jié) 語

    在自然界中,無論是海冰還是淡水冰,不同環(huán)境下冰體的形成機理各不相同,如在極地地區(qū)有1年冰和多年冰,不同年份的冰層的材料力學(xué)特性相差較大,此外,冰體一些主要的力學(xué)參數(shù)和力學(xué)性質(zhì)還受應(yīng)變率、溫度以及鹽度、雜質(zhì)等冰的形成成分因素的影響。因此,對冰船作用有限元數(shù)值模擬方法來說,開發(fā)一種涵蓋所有條件下的冰體材料數(shù)值模型是目前冰工程領(lǐng)域里一項極具挑戰(zhàn)性的研究工作,其中建立準確的冰材料本構(gòu)關(guān)系是其關(guān)鍵技術(shù)和難點問題,仍處于研究探索階段;除此之外,在船-冰碰撞數(shù)值模擬過程中,冰體變形破壞失效模式和裂紋擴展的真實模擬也是今后亟需解決的有限元數(shù)值模擬難點問題。

    本文總結(jié)了一些現(xiàn)有冰材料有限元數(shù)值計算模型,分析了各自的特點,并選取了幾種模型實驗工況,開展了有限元仿真計算,對比并分析了不同冰材料模型計算結(jié)果的差異性;此外,還總結(jié)了一些現(xiàn)有的冰材料數(shù)值模型計算可靠性驗證方法,針對一些算例,討論并分析了用于冰材料模型驗證的P-A曲線驗證方法和比吸能(SEA)判定方法。從實驗結(jié)果與數(shù)值模型計算結(jié)果對比可知,每種材料模型都有各自的適用范圍及前提條件,目前所開發(fā)的冰材料數(shù)值模型很難適用于各種實驗工況,利用模型實驗驗證方法來對冰材料模型準確性和可行性進行判斷具有一定的局限性。為此,考慮溫度、應(yīng)變率以及冰體形成條件等影響因素,建立不同實驗工況下的實驗結(jié)果數(shù)據(jù)庫,如SEA 值,為冰材料數(shù)值模型提供驗證數(shù)據(jù),是目前船-冰碰撞數(shù)值方法研究的關(guān)鍵之處。

    猜你喜歡
    冰體計算結(jié)果數(shù)值
    用固定數(shù)值計算
    彈體高速侵徹冰體研究
    數(shù)值大小比較“招招鮮”
    高速彈體侵徹冰材料過程數(shù)值模擬研究
    不等高軟橫跨橫向承力索計算及計算結(jié)果判斷研究
    甘肅科技(2020年20期)2020-04-13 00:30:40
    冰體質(zhì)量和撞擊角度對船首結(jié)構(gòu)碰撞性能的影響
    基于船-水-冰耦合技術(shù)的撞擊參數(shù)對船冰碰撞性能的影響
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    超壓測試方法對炸藥TNT當量計算結(jié)果的影響
    噪聲對介質(zhì)損耗角正切計算結(jié)果的影響
    欧美成人一区二区免费高清观看| 真人做人爱边吃奶动态| 内地一区二区视频在线| 啦啦啦观看免费观看视频高清| 丁香六月欧美| 亚洲av二区三区四区| 欧美日韩精品网址| 三级国产精品欧美在线观看| 色噜噜av男人的天堂激情| 国产成年人精品一区二区| 亚洲精品粉嫩美女一区| 亚洲无线在线观看| 欧美成狂野欧美在线观看| 在线观看日韩欧美| 丰满人妻熟妇乱又伦精品不卡| 美女被艹到高潮喷水动态| 国产毛片a区久久久久| 亚洲一区二区三区色噜噜| 国产野战对白在线观看| 日本精品一区二区三区蜜桃| 日日干狠狠操夜夜爽| 桃红色精品国产亚洲av| 国产单亲对白刺激| 麻豆一二三区av精品| 国产精品乱码一区二三区的特点| 老汉色∧v一级毛片| 91久久精品电影网| 欧美精品啪啪一区二区三区| 老汉色∧v一级毛片| 国产av在哪里看| 国产蜜桃级精品一区二区三区| 十八禁网站免费在线| 欧美bdsm另类| 久久99热这里只有精品18| 岛国视频午夜一区免费看| 国产精品日韩av在线免费观看| 亚洲欧美一区二区三区黑人| 老司机午夜十八禁免费视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 老汉色∧v一级毛片| 欧美在线黄色| 老汉色av国产亚洲站长工具| 欧美在线一区亚洲| tocl精华| 日韩欧美国产在线观看| 久久精品国产亚洲av涩爱 | 久久久久久大精品| 国产精品久久视频播放| 国产69精品久久久久777片| 午夜激情福利司机影院| 国产综合懂色| 岛国在线观看网站| 亚洲成人精品中文字幕电影| 国产亚洲欧美在线一区二区| 青草久久国产| 尤物成人国产欧美一区二区三区| 激情在线观看视频在线高清| 神马国产精品三级电影在线观看| 国产av在哪里看| 国产精品99久久久久久久久| 国产精品av视频在线免费观看| 亚洲国产精品合色在线| 五月伊人婷婷丁香| 欧美激情在线99| 欧美丝袜亚洲另类 | 国产精品 国内视频| 99riav亚洲国产免费| 两个人视频免费观看高清| 国产精品久久久久久亚洲av鲁大| 日韩欧美国产一区二区入口| 精品国产美女av久久久久小说| 黄片小视频在线播放| or卡值多少钱| 男人的好看免费观看在线视频| 国产成人av教育| 在线播放国产精品三级| 18+在线观看网站| 午夜福利在线观看免费完整高清在 | 黄片小视频在线播放| 深夜精品福利| 别揉我奶头~嗯~啊~动态视频| 日韩大尺度精品在线看网址| 神马国产精品三级电影在线观看| 午夜亚洲福利在线播放| 色尼玛亚洲综合影院| 女同久久另类99精品国产91| 精品日产1卡2卡| 亚洲中文日韩欧美视频| 国产午夜精品久久久久久一区二区三区 | 嫁个100分男人电影在线观看| 亚洲天堂国产精品一区在线| 亚洲成人免费电影在线观看| 国产一区二区亚洲精品在线观看| 久久精品国产99精品国产亚洲性色| 欧美精品啪啪一区二区三区| 国产成年人精品一区二区| 麻豆成人午夜福利视频| 免费观看精品视频网站| 国产av在哪里看| 搡老妇女老女人老熟妇| 精品久久久久久久久久久久久| 国产老妇女一区| 免费电影在线观看免费观看| 女警被强在线播放| 国产精品国产高清国产av| 日本 av在线| 18美女黄网站色大片免费观看| 级片在线观看| 动漫黄色视频在线观看| 床上黄色一级片| 亚洲人成电影免费在线| 国产午夜精品论理片| 日韩欧美精品v在线| 亚洲欧美一区二区三区黑人| 国产在视频线在精品| 亚洲电影在线观看av| 国产一区二区在线av高清观看| 国产麻豆成人av免费视频| 丁香六月欧美| 久久亚洲精品不卡| 精品欧美国产一区二区三| 国产精品电影一区二区三区| 美女大奶头视频| 日韩av在线大香蕉| 日韩大尺度精品在线看网址| 国产一区二区三区视频了| 变态另类丝袜制服| 成人鲁丝片一二三区免费| www日本在线高清视频| eeuss影院久久| 性色av乱码一区二区三区2| 成人鲁丝片一二三区免费| 亚洲av五月六月丁香网| 国产精品久久久人人做人人爽| 禁无遮挡网站| 欧美高清成人免费视频www| 日本与韩国留学比较| 免费人成视频x8x8入口观看| 级片在线观看| 亚洲美女黄片视频| 午夜免费激情av| 国产精品电影一区二区三区| av欧美777| 最后的刺客免费高清国语| 99热这里只有是精品50| 亚洲成人中文字幕在线播放| 亚洲精华国产精华精| 久久国产精品人妻蜜桃| 亚洲乱码一区二区免费版| 久久亚洲真实| 国产伦精品一区二区三区四那| 女人高潮潮喷娇喘18禁视频| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 男女午夜视频在线观看| 日本一二三区视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 深夜精品福利| 国产视频内射| 99久久九九国产精品国产免费| 美女高潮喷水抽搐中文字幕| 日本撒尿小便嘘嘘汇集6| 少妇熟女aⅴ在线视频| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 久久久久免费精品人妻一区二区| 免费电影在线观看免费观看| 成年女人毛片免费观看观看9| 国产亚洲精品久久久久久毛片| 深夜精品福利| 两人在一起打扑克的视频| 欧美日韩瑟瑟在线播放| 国产国拍精品亚洲av在线观看 | 亚洲中文字幕一区二区三区有码在线看| 综合色av麻豆| 欧美大码av| 此物有八面人人有两片| 国产精华一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩一区二区三| 成人av在线播放网站| 日韩国内少妇激情av| 久久久久久久久大av| 国产中年淑女户外野战色| 韩国av一区二区三区四区| 我的老师免费观看完整版| 高潮久久久久久久久久久不卡| 亚洲人成网站在线播| 高清在线国产一区| 婷婷精品国产亚洲av| 97超视频在线观看视频| 黄色片一级片一级黄色片| 亚洲va日本ⅴa欧美va伊人久久| 97超视频在线观看视频| 午夜激情福利司机影院| 特级一级黄色大片| 久久久久久久久久黄片| 免费看日本二区| 国产午夜福利久久久久久| 中文字幕熟女人妻在线| 国产精品久久电影中文字幕| 亚洲精品在线观看二区| 久久性视频一级片| 99国产精品一区二区三区| 久久久久久久久大av| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 久久香蕉国产精品| 中文字幕人妻熟人妻熟丝袜美 | 99在线视频只有这里精品首页| 老司机深夜福利视频在线观看| 99热只有精品国产| 亚洲国产精品久久男人天堂| 2021天堂中文幕一二区在线观| 成人无遮挡网站| 看黄色毛片网站| 免费av毛片视频| 热99re8久久精品国产| 激情在线观看视频在线高清| 免费看美女性在线毛片视频| 真人做人爱边吃奶动态| 亚洲18禁久久av| 特大巨黑吊av在线直播| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 一个人看视频在线观看www免费 | 99久久九九国产精品国产免费| 国产av不卡久久| 一卡2卡三卡四卡精品乱码亚洲| bbb黄色大片| 国内精品一区二区在线观看| a级毛片a级免费在线| 天堂影院成人在线观看| 少妇丰满av| 在线观看免费午夜福利视频| 变态另类成人亚洲欧美熟女| 欧美成人性av电影在线观看| 99热这里只有是精品50| 他把我摸到了高潮在线观看| 精品无人区乱码1区二区| 可以在线观看毛片的网站| 亚洲精品成人久久久久久| 国产高清有码在线观看视频| 亚洲天堂国产精品一区在线| 国产精品av视频在线免费观看| 久久久久性生活片| 我的老师免费观看完整版| 91久久精品国产一区二区成人 | 岛国视频午夜一区免费看| 国产v大片淫在线免费观看| 99久久久亚洲精品蜜臀av| 久久亚洲真实| 亚洲国产精品合色在线| 欧美日韩乱码在线| 人人妻人人看人人澡| 亚洲av免费在线观看| 成年女人毛片免费观看观看9| 欧美一区二区精品小视频在线| 丝袜美腿在线中文| 无限看片的www在线观看| 天堂√8在线中文| 欧美成人免费av一区二区三区| 91字幕亚洲| 在线观看午夜福利视频| 一个人看的www免费观看视频| 综合色av麻豆| 色尼玛亚洲综合影院| 午夜福利18| 99在线视频只有这里精品首页| 亚洲精品一卡2卡三卡4卡5卡| 天天添夜夜摸| 欧美高清成人免费视频www| 国产高清视频在线观看网站| 国产色爽女视频免费观看| 男人舔女人下体高潮全视频| 国产v大片淫在线免费观看| 欧美+亚洲+日韩+国产| 日韩高清综合在线| 免费无遮挡裸体视频| 免费看美女性在线毛片视频| 国产熟女xx| 久久精品影院6| 日韩欧美国产在线观看| 他把我摸到了高潮在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美av亚洲av综合av国产av| 国产视频内射| 婷婷精品国产亚洲av在线| 麻豆成人午夜福利视频| 亚洲成人久久性| 久久精品91无色码中文字幕| 欧美不卡视频在线免费观看| 久久6这里有精品| 搡老熟女国产l中国老女人| 有码 亚洲区| 制服人妻中文乱码| 色吧在线观看| 尤物成人国产欧美一区二区三区| 国产伦精品一区二区三区视频9 | 亚洲欧美日韩无卡精品| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看| 熟女人妻精品中文字幕| 亚洲成人久久性| 亚洲五月婷婷丁香| 欧美三级亚洲精品| 高清毛片免费观看视频网站| 日韩中文字幕欧美一区二区| 免费看十八禁软件| 好看av亚洲va欧美ⅴa在| 色综合亚洲欧美另类图片| 黄色日韩在线| 日韩欧美国产一区二区入口| 国产亚洲欧美在线一区二区| av黄色大香蕉| 99国产极品粉嫩在线观看| 亚洲av熟女| 熟女少妇亚洲综合色aaa.| 欧美激情久久久久久爽电影| 国产高清视频在线观看网站| 精品人妻一区二区三区麻豆 | 国模一区二区三区四区视频| 欧美日韩综合久久久久久 | 日韩欧美一区二区三区在线观看| 高清在线国产一区| 国产私拍福利视频在线观看| 亚洲天堂国产精品一区在线| 欧美黄色片欧美黄色片| 天堂网av新在线| 国产欧美日韩一区二区三| 日韩欧美在线二视频| 国产不卡一卡二| 又爽又黄无遮挡网站| 男女午夜视频在线观看| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 97超级碰碰碰精品色视频在线观看| 欧美一区二区精品小视频在线| 男人和女人高潮做爰伦理| 亚洲电影在线观看av| 久久久久九九精品影院| 999久久久精品免费观看国产| 男插女下体视频免费在线播放| 婷婷精品国产亚洲av| 母亲3免费完整高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 美女高潮喷水抽搐中文字幕| 欧美中文日本在线观看视频| 丝袜美腿在线中文| 免费看日本二区| 国产精品一区二区三区四区免费观看 | 少妇裸体淫交视频免费看高清| 午夜免费男女啪啪视频观看 | 亚洲色图av天堂| 一夜夜www| 老司机深夜福利视频在线观看| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看| 我要搜黄色片| 久久精品国产亚洲av涩爱 | 亚洲熟妇中文字幕五十中出| 中文字幕av成人在线电影| 国产黄a三级三级三级人| 青草久久国产| 国产一区二区三区视频了| 给我免费播放毛片高清在线观看| 999久久久精品免费观看国产| 亚洲成av人片免费观看| 一进一出抽搐gif免费好疼| 午夜精品在线福利| 精品电影一区二区在线| 午夜激情欧美在线| av中文乱码字幕在线| 成年免费大片在线观看| 亚洲国产精品999在线| 国产精品久久久久久久久免 | 亚洲欧美日韩无卡精品| 亚洲avbb在线观看| 99久久精品热视频| 亚洲人成网站高清观看| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 看黄色毛片网站| 久久久国产精品麻豆| 一区二区三区国产精品乱码| 国产一区二区三区在线臀色熟女| 国产探花极品一区二区| 最近最新免费中文字幕在线| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 欧美+日韩+精品| 我要搜黄色片| 国产乱人伦免费视频| 老鸭窝网址在线观看| 老司机深夜福利视频在线观看| 日本a在线网址| 国产午夜精品论理片| 亚洲av一区综合| 特级一级黄色大片| 男女视频在线观看网站免费| 黄色成人免费大全| tocl精华| 99国产极品粉嫩在线观看| 美女高潮的动态| 国产av在哪里看| 亚洲一区高清亚洲精品| 精品国产三级普通话版| svipshipincom国产片| 亚洲激情在线av| 久久国产精品人妻蜜桃| 久久久久久国产a免费观看| 亚洲午夜理论影院| 噜噜噜噜噜久久久久久91| 欧美乱色亚洲激情| 中亚洲国语对白在线视频| 亚洲中文字幕日韩| 国产成人aa在线观看| 国产三级中文精品| 国产单亲对白刺激| www.色视频.com| 亚洲av成人精品一区久久| 亚洲精品美女久久久久99蜜臀| 欧美最黄视频在线播放免费| 欧美性感艳星| 国产精品国产高清国产av| 亚洲精品久久国产高清桃花| 91久久精品电影网| 成人鲁丝片一二三区免费| 少妇丰满av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产av不卡久久| 色综合婷婷激情| 男女之事视频高清在线观看| 亚洲av二区三区四区| 精品一区二区三区av网在线观看| www.色视频.com| 亚洲av成人精品一区久久| 国产精品日韩av在线免费观看| 国产又黄又爽又无遮挡在线| 变态另类丝袜制服| 国产午夜精品论理片| 亚洲国产精品成人综合色| 一级黄色大片毛片| av女优亚洲男人天堂| 亚洲av成人不卡在线观看播放网| 成熟少妇高潮喷水视频| 日韩欧美一区二区三区在线观看| 成人一区二区视频在线观看| av在线天堂中文字幕| 亚洲最大成人中文| 18禁裸乳无遮挡免费网站照片| 免费在线观看影片大全网站| 欧美3d第一页| 91麻豆精品激情在线观看国产| 午夜精品一区二区三区免费看| 亚洲无线在线观看| 又黄又爽又免费观看的视频| 亚洲最大成人手机在线| 午夜福利在线在线| 亚洲久久久久久中文字幕| 欧美黄色片欧美黄色片| 在线国产一区二区在线| 亚洲av第一区精品v没综合| 在线观看66精品国产| 国产淫片久久久久久久久 | 久久婷婷人人爽人人干人人爱| 国产成人av激情在线播放| 国产一区二区三区在线臀色熟女| 国内精品一区二区在线观看| 成年人黄色毛片网站| 亚洲精品亚洲一区二区| 一卡2卡三卡四卡精品乱码亚洲| 丰满人妻一区二区三区视频av | www.999成人在线观看| 91在线观看av| 一区二区三区免费毛片| 1000部很黄的大片| 色视频www国产| www.色视频.com| 国产精品野战在线观看| 午夜亚洲福利在线播放| 全区人妻精品视频| 最近在线观看免费完整版| 亚洲成人久久性| 最近在线观看免费完整版| 欧美精品啪啪一区二区三区| 网址你懂的国产日韩在线| 超碰av人人做人人爽久久 | 国产老妇女一区| 一本久久中文字幕| 日本成人三级电影网站| 最近最新免费中文字幕在线| h日本视频在线播放| 午夜亚洲福利在线播放| 特级一级黄色大片| 免费人成视频x8x8入口观看| 国产精品99久久99久久久不卡| 露出奶头的视频| 精品久久久久久成人av| 亚洲精品乱码久久久v下载方式 | 欧美成人a在线观看| h日本视频在线播放| 中亚洲国语对白在线视频| 一个人看的www免费观看视频| 国产精品久久视频播放| 亚洲av二区三区四区| 18禁美女被吸乳视频| 少妇熟女aⅴ在线视频| 制服丝袜大香蕉在线| 欧美乱色亚洲激情| 天美传媒精品一区二区| 国产激情偷乱视频一区二区| 亚洲熟妇熟女久久| 成人一区二区视频在线观看| 麻豆国产97在线/欧美| 有码 亚洲区| 久久久久精品国产欧美久久久| av国产免费在线观看| 国产精品亚洲av一区麻豆| 精品久久久久久,| 欧美色欧美亚洲另类二区| 亚洲美女黄片视频| 久久草成人影院| 97人妻精品一区二区三区麻豆| 亚洲欧美精品综合久久99| avwww免费| 国产精品 欧美亚洲| 精品国产美女av久久久久小说| 天堂影院成人在线观看| 88av欧美| av专区在线播放| 亚洲av日韩精品久久久久久密| 亚洲av中文字字幕乱码综合| 噜噜噜噜噜久久久久久91| 9191精品国产免费久久| 天天添夜夜摸| 成人永久免费在线观看视频| 在线观看美女被高潮喷水网站 | 日本黄色片子视频| 色在线成人网| 中文字幕精品亚洲无线码一区| 大型黄色视频在线免费观看| 在线视频色国产色| 三级毛片av免费| 在线观看av片永久免费下载| 国产成人av教育| 91字幕亚洲| 中文字幕精品亚洲无线码一区| 神马国产精品三级电影在线观看| 99精品久久久久人妻精品| 精品电影一区二区在线| 国产淫片久久久久久久久 | 欧美大码av| 88av欧美| 国产精品一区二区免费欧美| 亚洲国产欧洲综合997久久,| 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| 国产久久久一区二区三区| 国产成人系列免费观看| e午夜精品久久久久久久| 村上凉子中文字幕在线| 欧美中文日本在线观看视频| 一级a爱片免费观看的视频| 亚洲精品乱码久久久v下载方式 | 亚洲av二区三区四区| 亚洲人成伊人成综合网2020| 精品国产美女av久久久久小说| 好看av亚洲va欧美ⅴa在| 男人和女人高潮做爰伦理| 中文字幕人成人乱码亚洲影| 小蜜桃在线观看免费完整版高清| 18禁裸乳无遮挡免费网站照片| 国产午夜精品久久久久久一区二区三区 | 99视频精品全部免费 在线| 国产日本99.免费观看| 亚洲av免费在线观看| 国内精品一区二区在线观看| 黄色片一级片一级黄色片| 91麻豆精品激情在线观看国产| 特级一级黄色大片| 日韩有码中文字幕| 午夜免费观看网址| 国产精品乱码一区二三区的特点| 亚洲成人久久爱视频| 久久久久久久久大av| 国产探花极品一区二区| 日韩有码中文字幕| 宅男免费午夜| 欧美日韩瑟瑟在线播放| 99热6这里只有精品| 88av欧美| 精品福利观看| 桃色一区二区三区在线观看| 深夜精品福利| 国产乱人伦免费视频| 国产亚洲欧美98| www.999成人在线观看| h日本视频在线播放| 国产69精品久久久久777片| 国产黄a三级三级三级人| 欧美日韩综合久久久久久 | 欧美成人a在线观看| 欧美日本亚洲视频在线播放| 看免费av毛片| 欧美成人a在线观看| 久久国产乱子伦精品免费另类| 国产激情欧美一区二区| 在线观看午夜福利视频| 精华霜和精华液先用哪个| 午夜日韩欧美国产| 1024手机看黄色片| 精品国产美女av久久久久小说| 99国产精品一区二区三区|