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

    寒區(qū)寬大結(jié)構(gòu)局部冰荷載特征研究

    2021-09-02 02:59:24王國軍資林欽岳前進張大勇
    船舶力學 2021年8期
    關鍵詞:海冰試件尺寸

    王國軍,彭 鑫,資林欽,許 寧,岳前進,張大勇,5

    (1.大連理工大學工程力學系,遼寧大連 116023;2.大連理工大學海洋科學與技術(shù)學院,遼寧盤錦 124221;3.昆明船舶設備研究試驗中心,昆明 650200;4.國家海洋環(huán)境監(jiān)測中心,遼寧大連 116023;5.大連理工大學寧波研究院,浙江寧波 315000)

    0 引 言

    冰荷載對結(jié)構(gòu)的尺寸和形狀非常敏感,因為不同的形狀與尺寸將導致不同的冰板失效破碎模式。直立結(jié)構(gòu)上,海冰的破碎模式比較多樣,在不同的寬厚比和加載速率組合下,冰板可能發(fā)生包括屈曲、徑向劈裂、環(huán)向開裂、擠壓、彎曲等破碎模式,其中海冰的擠壓破碎產(chǎn)生的荷載最大,對結(jié)構(gòu)的危害較大。一般根據(jù)海洋結(jié)構(gòu)水面處的尺寸大小,可以將其分為窄體結(jié)構(gòu)和寬體結(jié)構(gòu)。對于較寬的結(jié)構(gòu),海冰擠壓破壞形成荷載的過程相比窄體結(jié)構(gòu)更為復雜。寬大的抗冰設施(例如沉箱式、半潛式平臺等結(jié)構(gòu)),水線處剛度較大,研究人員需要關注結(jié)構(gòu)的整體和局部冰載荷,即總體及局部冰壓力,開展平臺的定位能力和結(jié)構(gòu)局部、整體結(jié)構(gòu)強度校核,確保結(jié)構(gòu)滿足強度要求。

    Ashby等[1]提出了海冰與結(jié)構(gòu)相互作用時存在非同時破壞,隨著與結(jié)構(gòu)接觸面積的增加,結(jié)構(gòu)的整體冰壓力隨之減小。Sanderson[2]、Wright[3]和Blanchet[4]等基于現(xiàn)場原型測量數(shù)據(jù)分析發(fā)現(xiàn),海冰存在非同時破壞的現(xiàn)象。目前抗冰結(jié)構(gòu)設計中,由于海冰強度的離散性較大,通常將海冰的破壞強度取值為平均強度(有效值),整個結(jié)構(gòu)受到的總冰荷載為接觸面積與該有效值的乘積,如ISO 19906 中建議的直立結(jié)構(gòu)上的整體冰壓力計算公式(式1)[5],此式得到的冰壓力為整個接觸面積上的平均冰壓力。

    式中:pG是冰與結(jié)構(gòu)作用的名義接觸面積上的單位面積冰壓力;CR是冰壓縮強度系數(shù);w為結(jié)構(gòu)寬度;h為海冰厚度;h1為參考厚度;m為經(jīng)驗系數(shù);n也為經(jīng)驗系數(shù),當h<1.0 m時,n=-0.5+h/5;當h≥1.0 m時,n=-0.3。此公式適用于寬厚比w/h大于2的剛性結(jié)構(gòu)。

    海冰與直立結(jié)構(gòu)相互作用時,有些局部接觸區(qū)域面上的壓力較高,主要原因有:(1)海冰與整個結(jié)構(gòu)界面的接觸存在非同時性以及海冰破碎的非同時性;(2)局部接觸點上的海冰破碎強度不完全相同。Timco 等[6]對整體冰壓力與局部冰壓力進行了定義,將海冰與結(jié)構(gòu)的整個接觸面分為多個較小的面,每個小面元上的平均冰壓力為局部冰壓力,整個接觸面上的平均冰壓力為結(jié)構(gòu)的整體冰壓力,如圖1所示。在海洋結(jié)構(gòu)總體設計和穩(wěn)定性設計方面,通常要考慮整體壓力,當結(jié)構(gòu)設計中需要考慮局部面元的安全性時,需要確立局部面元所承受的極限荷載。因此,結(jié)構(gòu)的局部冰壓力為面元高壓區(qū)的最大冰壓力,本文采用以上的整體冰壓力與局部冰壓力的定義方法。

    圖1 整體冰壓力和局部冰壓力定義示意圖Fig.1 Illustration of the definitions of global and local pressures

    1988 年,Sanderson[7]通過分析室內(nèi)模型試驗、現(xiàn)場監(jiān)測不同類型結(jié)構(gòu)上的冰力數(shù)據(jù),發(fā)現(xiàn)隨著海冰與結(jié)構(gòu)接觸面積的增加,平均冰壓力呈現(xiàn)減小的趨勢,冰壓力在越小的接觸面積上越大;Masterson[8]對海洋平臺現(xiàn)場測量和船舶實驗中的局部冰壓力數(shù)據(jù)進行了分析,得到了隨著接觸面積的增加而冰壓力降低的趨勢,通過分析接觸面積大于100 m2的結(jié)構(gòu)的測量數(shù)據(jù),形成了冰區(qū)結(jié)構(gòu)的局部冰壓力分析方法,同時發(fā)現(xiàn)隨著局部接觸面積的增加局部冰壓力隨之減?。籉rederking等[9]利用Molikpaq平臺的測量數(shù)據(jù),建立了以局部接觸面積為變量的局部冰壓計算方法,獲得了局部冰壓力的垂直和水平分布特征;Frederking[10]根據(jù)日本海洋工業(yè)協(xié)會1999 年的海冰模型試驗數(shù)據(jù),分析了冰壓力分布的詳細信息,發(fā)現(xiàn)局部壓力與整體壓力的變化趨勢相同;K?rn? 等[11]提取了波羅的海燈塔局部壓力極值的現(xiàn)場數(shù)據(jù),通過分析發(fā)現(xiàn)冰破碎引起的冰壓力隨冰厚的增加而減小,通過轉(zhuǎn)換獲得一組標準化的最大壓力,使極值分析中使用的數(shù)據(jù)值不依賴于冰厚;Masterson 等[12]通過對大量現(xiàn)場測量數(shù)據(jù)的分析,提出了一種修正的用于孤立小面積區(qū)域的冰壓力-面積的關系,即局部壓力與面積的負指數(shù)關系;Palmer等[13]對冰擠壓過程的高壓區(qū)應力做了新的假設,重新討論局部區(qū)域的定義,完善了Masterson 的結(jié)論;Timco 等通過對Sanderson 的壓力-面積曲線和數(shù)據(jù)進一步分析,明確了局部冰壓力區(qū)的定義,考慮了包括面積在內(nèi)的冰速、寬厚比、冰破壞模式和冰力學性質(zhì)對壓力的影響,結(jié)果表明并非所有數(shù)據(jù)都服從壓力-面積曲線[6];田喜民等[14]總結(jié)分析了冰區(qū)航行船舶的整體冰力和局部冰力與面積的關系,對比了現(xiàn)有冰荷載理論與數(shù)值分析方法在船舶分析上的應用;Su 等[15]采用數(shù)值方法分析了船舶行駛過程中局部冰載荷特征;龔榆峰[16]基于斷裂力學和損傷力學的理論建立了冰荷載的直接計算方法,采用有限元軟件分析了極地船舶的冰載荷;Timco等[17]通過研究Molikpaq 海域沉箱式寬大結(jié)構(gòu)冰荷載大量測量數(shù)據(jù),總結(jié)了影響冰力的相關參數(shù),提出研究寬大結(jié)構(gòu)的冰壓力隨結(jié)構(gòu)寬度變化是否合理、結(jié)構(gòu)上是否存在冰壓力的上限的疑問;龍雪[18]基于離散元模擬分析了直立結(jié)構(gòu)上的冰壓力分布特性,并深入研究了高壓區(qū)的產(chǎn)生機理;Taylor 等[19]提出了一種模擬局部高壓區(qū)域載荷的概率模型;孔令學[20]通過研究《港口工程荷載規(guī)范》中的冰荷載計算公式,認為計算方法對大尺度結(jié)構(gòu)的局部擠壓系數(shù)考慮比較簡單。目前對直立結(jié)構(gòu)局部冰力分析方法的研究成果主要基于現(xiàn)場結(jié)構(gòu)測量數(shù)據(jù)分析得到,同時也有部分是基于室內(nèi)試驗數(shù)據(jù),由于結(jié)構(gòu)形式比較單一,研究結(jié)論主要還是考慮局部區(qū)域的面積為影響因素,忽略了冰厚與結(jié)構(gòu)寬度相對值的影響,同時不同研究成果也存在一定的差異性。

    造成海冰與直立結(jié)構(gòu)相互作用存在局部高壓區(qū)的原因眾多,海冰力學特性、冰與結(jié)構(gòu)接觸的非同時性是兩大主要原因。本文通過海冰單軸壓縮強度實驗,分析了海冰的壓縮強度與試件尺寸的關系;通過對比現(xiàn)有的局部冰壓力計算方法,分析不同計算方法存在差異性的原因;基于有限元和離散元方法將不同尺寸的結(jié)構(gòu)劃分為多個有限面元,研究冰厚、接觸寬度對局部冰壓力大小的影響,分析海冰與結(jié)構(gòu)相互作用的冰壓力分布特征,最終提出寬大結(jié)構(gòu)局部冰壓力預測模型。

    1 海冰尺寸對壓縮強度的影響分析

    海冰內(nèi)部存在各種各樣的缺陷,其在形成過程中會產(chǎn)生鹽水氣泡、固態(tài)鹽等,成長中晶體粒徑也大小不一,海冰的平均粒徑范圍在1~11 mm 之間[21],同時在海冰內(nèi)部也存在大小不一的裂紋。因此,海冰的尺寸越大,其包含的缺陷越多,在外力作用下越易于各個缺陷的累積從而造成海冰的破壞。Jones 等[22]通過研究試件尺寸和晶粒尺寸對冰壓縮強度的影響實驗,發(fā)現(xiàn)當試件的尺寸為晶粒尺寸的12 倍以上時才不會對冰的壓縮強度產(chǎn)生影響。張麗敏[21]通過分析平均粒徑D與試件尺寸A比值的變換(D/A)對冰壓縮強度的影響,得到當D/A<1.5時,冰的單軸壓縮強度隨D/A的增大而減??;當D/A>1.5時,冰的單軸壓縮強度無明顯變換趨勢。沈梧等[23]通過對葫蘆島和鲅魚圈海域海冰的壓縮強度、拉伸強度及斷裂韌性隨試件尺寸變化的分析,得到海冰強度均反映出不同程度的尺寸效應,海冰的尺寸對其強度有一定的影響。

    1.1 實驗方法及結(jié)果

    為了進一步明確海冰的尺寸效應特征,本文對不同尺寸海冰進行單軸壓縮強度實驗,海冰取自黃海北部莊河附近海域,黃渤海海域的海冰以柱狀冰為主,此次采集的海冰也為柱狀冰。將采集得到的大塊冰坯運送至實驗室,利用油鋸對大塊冰坯進行粗加工,而后放入冷柜中保溫。

    由于柱狀冰表現(xiàn)出各向異性的力學性質(zhì),基于加載方向與冰層生長方向的夾角分為水平加載(加載方向與生長方向垂直)和垂直加載(加載方向與生長方向平行)兩種實驗加載方法。根據(jù)1988年國際水利工程師協(xié)會(IAHR)冰工程會議上提出的冰實驗推薦方法,海冰單軸壓縮實驗標準試件的截面尺寸為70×70 mm,試件高度為175 mm。

    為了分析海冰的尺寸效應,本文定義了兩種海冰尺寸變化限制條件:一是保持試件高寬比為2.5(即標準試件的高寬比),改變試件的體積大??;二是保持試件的高度為175 mm(即標準試件的高度),改變試件的截面尺寸。設計的海冰試件尺寸如表1 所示,每個尺寸制作4 個試件,試件通過實驗室的臺鋸進行精細加工后存儲于低溫柜(-10°C)中24 h以上。

    表1 海冰單軸壓縮強度尺寸設計Tab.1 Size design of sea ice uniaxial compression strength

    本次實驗全程在低溫環(huán)境實驗室內(nèi)進行,室內(nèi)可控環(huán)境溫度精度為±0.1℃,室內(nèi)溫度控制在-10±1℃。采用WDW-50E 型微機控制電子萬能試驗機進行實驗加載,每次加載實驗前采集海冰的溫度、實際長寬高和質(zhì)量數(shù)據(jù)。溫度采用非接觸式激光溫度槍測量試件3處以上不同位置的溫度計算平均值;長寬高采用分辨率為1 mm 的角尺進行測量;質(zhì)量采用分辨率為0.01 g 的電子天平測量(如圖2 所示)。

    圖2 海冰加載實驗Fig.2 Sea ice loading experiment

    本次海冰單軸壓縮實驗中在試件的上下端均放有橡膠墊塊,用于保障加載過程中海冰試件的受力均勻以及防止試件端部的破碎和融化。記錄實驗的力與變形的變化曲線,提取極限應力作為試件的壓縮強度,加載實驗結(jié)束后采集此試件不同位置的碎冰用于測量試件的鹽度,海冰的密度通過質(zhì)量與體積的比值計算,海冰的主要物理參數(shù)如表2所示。

    表2 海冰試件密度、鹽度、溫度的均值Tab.2 Mean values of density,salinity and temperature of sea ice specimens

    海冰試件單軸壓縮實驗分為兩種尺寸,其中高寬比為2.5的海冰壓縮強度實驗結(jié)果見表3,高度為定值的試件主要實驗結(jié)果見表4。

    表3 試件體積因素尺寸效應試驗結(jié)果Tab.3 Test results of size effect of volume factor of specimen

    表4 試件橫截面積因素尺寸效應試驗結(jié)果Tab.4 Test results of size effect of cross-sectional area factors of specimens

    1.2 結(jié)果分析

    (1)試件體積的尺寸效應

    基于實驗數(shù)據(jù)分析冰試件高度與寬度(橫截面邊長)比值保持為2.5 不變的條件下海冰壓縮強度的變化趨勢,此時將試件的體積作為變量,體積的變化范圍是2.85×105~2.56×106mm3。考慮柱狀冰物理力學性質(zhì)的各向異性,分別進行了垂直加載與水平加載的海冰單軸壓縮實驗,試驗結(jié)果如圖3和圖4所示。

    圖3 海冰試件垂直加載壓縮強度與體積的關系Fig.3 Relation between compressive strength and volume of sea ice under vertical loading

    圖4 海冰試件水平加載壓縮強度與體積的關系Fig.4 Relation between compressive strength and volume of sea ice under horizontal loading

    由圖3~4分析可知,海冰壓縮強度與體積呈現(xiàn)線性遞減的關系。無論垂直加載還是水平加載,海冰的壓縮強度都隨試件體積的增大而減小,其中垂直加載的減小速率要大于水平加載的速率,海冰體積的增加影響著整體強度的變化。

    (2)試件橫截面的尺寸效應

    試驗中在冰試件高度保持為175 mm 不變的條件下,改變試件橫截面的邊長,橫截面為正方形,變化因素可看作是試件的橫截面積大小,橫截面積變化范圍是2 500~14 400 mm2。同樣考慮了冰的各向異性,分別進行了垂直加載與水平加載的海冰單軸壓縮實驗,試驗結(jié)果如圖5和圖6所示。

    圖5 海冰試件垂直加載壓縮強度與面積的關系Fig.5 Relation between compressive strength and area of sea ice under vertical loading

    圖6 海冰試件水平加載壓縮強度與面積的關系Fig.6 Relation between compressive strength and area of sea ice under horizontal loading

    由圖5~6 分析可知,海冰壓縮強度與試件截面積呈現(xiàn)出線性遞減的關系。無論垂直加載還是水平加載,海冰的壓縮強度均隨截面積的增大而減小,其中垂直加載的減小速率要大于水平加載的速率,水平加載的離散度要大于垂直加載的試件。海冰體積的增加影響著海冰的整體強度變化。

    基于以上的分析發(fā)現(xiàn),海冰體積和橫截面的變化均對海冰的壓縮強度有較大影響,海冰單軸壓縮強度具有明顯的尺寸效應,海冰在越小的局部面積上具有越大的局部冰壓力。如果結(jié)構(gòu)設計中忽略海冰對局部強度的影響,對于水面處尺寸較小的結(jié)構(gòu)影響較小,但對于水面處的結(jié)構(gòu)尺度較大及水面處對荷載反應明顯的單元(如殼體結(jié)構(gòu))則影響較大,基于整體冰荷載計算方法的平均冰壓力較小,而海冰壓縮強度存在的尺寸效應,接觸面的局部壓力遠大于平均冰壓力,因此需要研究寬大結(jié)構(gòu)上局部冰壓力的影響因素,明確局部冰壓力的預測模型。

    2 寬大結(jié)構(gòu)局部冰壓力研究

    由以上對局部冰壓力的分析可以發(fā)現(xiàn),結(jié)構(gòu)越寬基于整體冰力加載的平均冰壓力就越小。而一些學者通過現(xiàn)場監(jiān)測數(shù)據(jù)分析發(fā)現(xiàn),冰與較寬的結(jié)構(gòu)相互作用時,整個接觸面不同區(qū)域冰的破碎存在非同時性,在小的局部區(qū)域則存在同時破碎,在此局部區(qū)域的冰壓力就會較高,從而冰壓力可能遠大于整個接觸面積上的平均冰壓力,因此在結(jié)構(gòu)的局部強度設計中需要考慮局部冰載荷的作用。

    冰與直立結(jié)構(gòu)相互作用時,擠壓破壞是較為常見的破壞形式。在冰擠壓破碎過程中,相互作用區(qū)具有三個不同的壓力區(qū):臨界區(qū)、基本壓力區(qū)和初始破碎區(qū)。臨界區(qū)可定義為短時間內(nèi)發(fā)生強烈壓力的局部冰區(qū),臨界區(qū)對破碎過程影響較大。臨界區(qū)域的出現(xiàn)與相互作用的規(guī)?;蝾愋蜔o關,盡管臨界區(qū)具有高度隨機性,但基本參數(shù)(如區(qū)域大小、力、壓力和空間密度)是可量化的。臨界區(qū)面積約為0.10 m2,可以施加0.1~4.0 MN 的力。臨界區(qū)空間密度為單位面積的區(qū)域數(shù),范圍約為0.6~0.8 個/m2,受到海冰與結(jié)構(gòu)尺度效應的影響。臨界區(qū)解釋了平均壓力隨接觸面積增加而降低的現(xiàn)象,存在兩種壓力-面積關系,一種是隨時間增加而增加的接觸面積,另一種是較大的、無約束的接觸面積(包含較小的、高度受限的區(qū)域)[24]。

    國際上通過現(xiàn)場測量數(shù)據(jù)和部分室內(nèi)試驗研究,建立了一系列局部冰壓力計算方法,多數(shù)研究成果認為海冰的局部壓力隨著接觸面積的增大而減小。Timco通過研究發(fā)現(xiàn)該結(jié)論也存在一定的不足,認為在一定的局部面積上,由于相互作用情況的不同,壓力不一定隨面積的增大而減小[2];資林欽基于離散元方法對局部冰力進行了一定的研究,初步獲得了局部冰力與冰厚和局部區(qū)域面積的關系[24],但缺乏對現(xiàn)有局部冰力計算方法差異性的分析和對分析方法的驗證分析。因此需要進一步分析局部冰壓力的影響因素,考慮多種因素對局部冰壓力的影響。

    2.1 局部冰力計算方法對比

    目前局部冰力計算方法大多是以局部區(qū)域面積為變量的計算公式,下文列舉出幾個典型的局部冰力公式,同時分析各個計算方法之間的差異性。

    ISO 19906規(guī)范中給出了考慮不同結(jié)構(gòu)條件時的局部冰力的計算方法[1]:

    (1)局部冰力采用如下公式計算:

    式中:ad為局部區(qū)域高度,要求ad大于0.14 m,且ad≤0.4h,h為冰厚;wL為局部區(qū)域?qū)挾?,要求wL/ad≤10。如果ad>0.4h,則采用式(2)計算。

    (2)全冰厚冰壓力計算公式如下:

    式中,當h>0.35 m 時使用該式,當h≤0.35 m 時,pF=4.0 MPa。冰壓力主要分布在中心區(qū)域內(nèi),局部冰壓強可由式(4)計算得到。

    式中,pL為局部冰壓力,γL取2.5。

    (3)對于冰厚大于1.5 m的情況,給出局部冰壓力的計算公式如下:

    式中,A為局部區(qū)域面積,當A≤10 m2時使用該式,當A>10 m2時,pL=1.48 MPa。

    DNV規(guī)范中對于局部冰壓力的計算公式如下[25]:

    式中,σc為海冰壓縮強度,規(guī)范中規(guī)定pL的最大值應小于20 MPa。

    API和CSA采用Masterson和Frederking(1993)的研究結(jié)果[12],建議對于面積不超過19 m2的區(qū)域的冰壓力計算公式為

    式中,A為局部面積,當A>19 m2時,局部冰壓力pL=1.5 MPa。

    結(jié)合Sanderson 的實驗數(shù)據(jù)[2],本文分析了不同局部區(qū)域面積及不同海冰厚度條件下各種冰壓力計算方法的差異性,結(jié)果如圖7~8 所示。其中,圖7 是基于海冰厚度為0.5 m,改變局部區(qū)域的寬度(0.5~20 m)的結(jié)果;圖8 考慮了0.5 m 局部寬度下改變海冰厚度(冰厚0.1~0.3 m)和6 m 局部寬度下改變海冰厚度(冰厚0.4 ~1.6 m)。

    圖7 冰厚不變、改變局部寬度下的局部冰力Fig.7 Local ice force with constant ice thickness and varying local width

    圖8 局部寬度不變、改變冰厚下的局部冰壓力Fig.8 Local ice pressure with constant local width and varying ice thickness

    由圖7分析可以發(fā)現(xiàn):

    (1)ISO規(guī)范中的計算方法在計算冰厚0.35 m 到1.5 m之間的局部強度時,只考慮了冰厚的影響,不能夠反應局部區(qū)域?qū)挾茸兓挠绊?。當局部面積小于2 m2時,計算結(jié)果遠小于Sanderson 實驗數(shù)據(jù)的結(jié)果;當局部面積大于2 m2時,計算結(jié)果略大于實驗數(shù)據(jù)。

    (2)API/CSA 規(guī)范的計算方法是基于局部面積,其計算結(jié)果均大于實驗數(shù)據(jù)的結(jié)果,其中局部面積小于1 m2時計算結(jié)果與實驗數(shù)據(jù)相差更大。

    (3)DNV 規(guī)范的計算方法考慮了冰厚和局部面積,在計算局部面積小于1 m2時計算結(jié)果小于實驗數(shù)據(jù)的結(jié)果,但在局部面積大于1 m2時與實驗數(shù)據(jù)比較接近。

    由圖8分析可知:

    (1)ISO規(guī)范的計算方法在計算冰厚小于0.35 m的局部強度時,局部強度均為4 MPa,不能夠反應局部區(qū)域?qū)挾仍黾拥挠绊懀裨?.35 m 到1.5 m之間時局部冰壓力隨冰厚的增加而減小。

    (2)API/CSA 規(guī)范的計算方法是基于局部面積計算,在計算局部面積較小時計算結(jié)果偏大,例如,對于0.05 m2的局部面積,局部冰壓力計算值可以達到36 MPa。

    (3)DNV 規(guī)范的計算方法考慮了冰厚和局部面積,在局部寬度不變的情況下,隨著冰厚的增大局部冰壓力也在增大,能夠較好地反映海冰與結(jié)構(gòu)相互作用的特征,但其計算結(jié)果在面積小時偏小,面積大時偏大。

    計算局部冰壓力時同時考慮冰厚和局部區(qū)域面積的影響能夠更好地反映局部冰壓力的變化特征。

    2.2 基于離散元的冰壓力分析

    目前現(xiàn)場獲得的局部冰壓力數(shù)據(jù)較少,涉及的結(jié)構(gòu)寬度范圍有限,同時不同的計算方法得到的局部冰壓力計算結(jié)果差異較大。因此,本文通過ANSYS 有限元軟件將不同寬度結(jié)構(gòu)劃分為小的面元,然后基于離散元軟件(采用大連理工大學季順迎團隊開發(fā)的SDEM 離散元分析軟件[26-28])計算海冰與結(jié)構(gòu)的相互作用,獲得局部面元上的冰壓力,研究局部冰壓力的計算方法,軟件通過將海冰離散為具有一定質(zhì)量和大小、具有黏結(jié)-破碎功能的顆粒單元,單元間具有相應的粘接作用力。圖9 為以Norstr?msgrund燈塔的結(jié)構(gòu)特征為例所進行的局部冰壓力分析。

    圖9 Norstr?msgrund燈塔面元劃分和離散元模型Fig.9 FEM model and DEM model of Norstr?msgrund lighthouse

    Norstr?msgrund 燈塔水線處的結(jié)構(gòu)為圓柱體,直徑為7.52 m,燈塔上安裝的壓力盒由17 個大小為1.2 m×1.6 m 的小壓力盒組合而成(編號Panel-1 到Panel-9 和Panel-91 到Panel-98),其中獲得的兩次不同測量時間的冰力數(shù)據(jù)見表5,兩次數(shù)據(jù)對應的冰況均為:冰速0.15 m/s,冰厚0.6 m。

    表5 Norstr?msgrund燈塔壓力盒測量數(shù)據(jù)Tab.5 Measurement data of Norstr?msgrund lighthouse pressure box

    離散元模擬中將單個面元大小劃分為0.5 m×0.5 m,其中海冰運動速度為0.15 m/s,冰厚為0.6 m。提取結(jié)構(gòu)的冰力數(shù)據(jù)(見圖10和圖11),通過后期數(shù)據(jù)處理,利用Tecplot軟件繪制燈塔的局部壓力分布圖(如圖12所示),本次模擬得到的整體結(jié)構(gòu)峰值冰力為1 808.43 kN,與實測的冰力峰值相近(見表5)。

    圖10 SDEM模擬的海冰與燈塔相互作用破碎過程Fig.10 Interaction between sea ice and lighthouse simulated by SDEM

    圖11 SDEM模擬得到的燈塔整體冰力時程Fig.11 Ice force time history of lighthouse simulated by SDEM

    圖12 Norstr?msgrund燈塔模擬的局部冰壓力分布Fig.12 Simulated local ice pressure distribution

    將上述局部冰壓力轉(zhuǎn)換為與壓力盒大?。?.2 m×1.6 m)相當?shù)哪P途植棵娣e上(1.5 m×1.5 m),得到了圖13 的對比結(jié)果,對比Panel 上的實測壓力盒平均壓力與模擬的平均壓力,模擬值要大于實測值。由于模擬值為壓力最大值區(qū)域的平均壓力,而且單個面元的劃分較小(0.5 m 直角邊的等腰直角三角形),從而造成了局部壓力值較大,在實際測量中壓力盒測量的區(qū)域較大可能會造成平均壓力較小。通過對比可知離散元分析方法能夠滿足對局部冰壓力的預測。

    圖13 實測的局部冰壓力與模擬值對比Fig.13 Comparison of measured local ice pressure with simulated value

    2.3 局部冰壓力預測模型研究

    通過模擬不同寬度結(jié)構(gòu)與海冰相互作用,獲得結(jié)構(gòu)在不同大小局部面積區(qū)域的冰壓力,同時對比ISO 19906(2010)中提供的現(xiàn)場實測數(shù)據(jù),分析結(jié)果如圖14所示,離散元計算的結(jié)果與實測數(shù)據(jù)具有相同的分布情況,基于離散元分析計算局部冰壓力具有可行性與可信度。

    圖14 數(shù)值模擬和實測數(shù)據(jù)的局部冰壓力與面積關系Fig.14 Local ice pressure of simulated values and measured data versus area

    由圖14 可知,局部冰壓力隨著接觸面積的增大而減小,兩者成指數(shù)關系,ISO、API 和CSA 等規(guī)范中以及一些學者的研究表明局部冰壓力與局部面積為冪指數(shù)關系。本文基于冪指數(shù)函數(shù)對離散元數(shù)值模擬的結(jié)果進行擬合,分析得到如下公式:

    式中,pL為局部冰壓力,A為局部區(qū)域面積,離散元分析中局部面積A通過壓力較大單元附件多個單元面積相加獲得不同局部面積,局部壓力為這些單元的平均壓力。通過對比可知,以上的擬合結(jié)果與ISO、API和CSA規(guī)范中的局部冰壓力存在一定的差異,但與Pond Inlet的試驗數(shù)據(jù)比較吻合,這也說明僅考慮局部面積作為局部冰壓力的計算變量不夠全面。

    由2.1 節(jié)的分析可知,局部冰壓力由冰厚和局部面積兩個變量共同影響,基于這兩個變量構(gòu)建的預測模型能夠較好地反映局部冰壓力的變化趨勢。因此,需要進一步研究局部壓強、局部面積和冰厚三者之間的關系,同時海冰的壓縮強度也直接影響局部冰壓力的大小。

    分析冰與結(jié)構(gòu)相互作用時,常引入冰厚與結(jié)構(gòu)寬度的比值進行分析。本文引入h/wL無量綱參數(shù)分析局部冰壓力,其中h為冰厚,wL為局部區(qū)域?qū)挾龋梃bDNV的局部冰力分析公式,引入海冰的壓縮強度,從而建立了如下關系式:

    式中,σc為海冰壓縮強度?;谑剑?)對離散元計算得到的局部冰壓力數(shù)據(jù)進行分析,擬合結(jié)果如圖15所示。

    圖15 局部冰壓力與h/wL的關系Fig.15 Relationship between local ice pressure and h/wL

    經(jīng)過擬合得到pL與h/wL的關系式為

    通過對比建立的擬合公式(10)與Pond Inlet 的試驗數(shù)據(jù)[2(]如圖16 所示),可以發(fā)現(xiàn)構(gòu)建的計算公式與試驗的局部冰壓力數(shù)據(jù)吻合較好,擬合結(jié)果與實測值的上限值變化趨勢相近,能夠保守地評價局部冰壓力。

    圖16 離散元結(jié)果與實測數(shù)據(jù)對比Fig.16 Comparison between DEM results and measured data

    當局部區(qū)域的寬度大于冰厚(即wL>h)時,可將式(10)改為

    式中,AL為局部區(qū)域的面積,綜合考慮冰厚及面積因素的局部冰壓力可由式(11)計算。

    2.4 局部冰壓力分布特征分析

    以寬度為10 m 的結(jié)構(gòu)為例,其冰情特征為:冰厚0.5 m,冰速0.2 m/s,研究局部冰壓力在結(jié)構(gòu)寬度和冰厚方向上的分布特征,圖17為不同面元邊長的局部冰壓力分布特征。

    由圖17 可知,在冰與結(jié)構(gòu)的整體接觸區(qū)域內(nèi)結(jié)構(gòu)冰壓力分布是不均勻的,在一些區(qū)域會出現(xiàn)高壓區(qū)。由兩圖對比可看出,不同單元大小所得到的局部應力分布情況是不同的,單元越小,局部應力越集中。可以推知在局部冰力集中位置,冰會發(fā)生局部同時破壞,冰與結(jié)構(gòu)接觸的邊緣位置由于截面突變會導致應力集中,寬結(jié)構(gòu)局部冰力分布隨機性很強。

    圖17 局部冰壓力分布特征Fig.17 Local ice pressure distribution

    選取計算時間內(nèi)某個峰值冰力對應時刻,提取此時刻的局部冰力分布數(shù)據(jù),冰厚中心位置區(qū)域的局部冰力(取面積0.25 m2)沿結(jié)構(gòu)水平方向作圖,如圖18所示,在10 m寬的結(jié)構(gòu)上,高壓區(qū)局部冰力較大,結(jié)構(gòu)整體受到的冰載荷主要由高壓區(qū)傳遞。

    圖18 局部冰壓力在結(jié)構(gòu)寬度方向的分布Fig.18 Distribution of local ice pressure in the width direction of structure

    同時將上述時刻的局部冰力(取面積0.25 m2)沿冰厚方向作圖,如圖19 所示,將縱軸設計為離冰厚中心線的距離與冰厚的比值z/h。通過繪制冰壓力分布發(fā)現(xiàn):在冰厚方向上,冰板中心位置壓力值最大,沿冰厚方向降低,在冰板邊緣處降低為0,在中心位置的海冰擠壓時受到冰板自身的限制,力學強度增大,導致中心位置局部冰壓力較大,同時實際海冰的中心位置比較平整,更容易發(fā)生同時破碎;靠近邊緣的海冰受到擠壓時能夠向上下剝落清除,壓力值會較低。局部冰壓力的冰厚分布與ISO19906中給出的規(guī)律一致。

    圖19 局部冰壓力在冰厚方向的分布Fig.19 Distribution of local ice pressure in the direction of ice thickness

    3 結(jié) 論

    本文通過分析海冰壓縮強度與試件尺寸的關系發(fā)現(xiàn):海冰在垂直加載與水平加載兩種方式下,其壓縮強度都隨海冰尺寸的增大而減??;而海冰垂直加載的減小率要大于水平加載的減小率,海冰試件體積的增加造成的壓縮強度減小率要大于海冰截面積增加造成的壓縮強度減小速率。這表明海冰的壓縮強度具有明顯的尺寸效應,海冰體積越大時其內(nèi)部的缺陷越多,在外力的作用下更加容易發(fā)生破碎,也表明了海冰在小的局部面積上具有更大的局部冰壓力。

    由于海冰與直立結(jié)構(gòu)相互作用時存在局部高壓力區(qū),通過分析現(xiàn)有的局部冰壓力計算方法,發(fā)現(xiàn)不同計算方法對相同局部面積的冰壓力分析結(jié)果差異較大。ISO 規(guī)范的計算方法僅考慮了冰厚的影響,同時對局部面積小于2 m2的計算結(jié)果與實驗數(shù)據(jù)偏差較大;API/CSA 規(guī)范的計算方法僅考慮了局部面積的影響,同時局部面積小于1 m2時計算結(jié)果遠大于實驗數(shù)據(jù);DNV規(guī)范同時考慮了冰厚和局部面積的影響,能夠較好反映局部冰壓力變化趨勢,但是與實驗數(shù)據(jù)相比其計算結(jié)果在面積小時偏小,面積大時偏大。

    本文基于Norstr?msgrund燈塔壓力盒測量的局部冰力數(shù)據(jù),研究了基于離散元軟件的局部冰壓力分析方法,通過計算不同寬度結(jié)構(gòu)上的局部冰壓力,構(gòu)建了考慮冰厚和局部面積雙因素的寬大結(jié)構(gòu)局部冰壓力預測模型,通過與實測數(shù)據(jù)對比,兩者吻合較好。同時局部冰壓力在冰厚方向上的分布呈現(xiàn)出中心位置壓力值最大,沿鉛垂方向降低,到冰板邊緣處降低為0,局部冰壓力在冰厚方向上的分布與ISO19906中給出的規(guī)律一致。

    猜你喜歡
    海冰試件尺寸
    尺寸
    智族GQ(2022年12期)2022-12-20 07:01:18
    復材管纖維纏繞角度對約束混凝土軸壓性能的影響研究
    廣東建材(2022年5期)2022-06-10 08:48:20
    CIIE Shows Positive Energy of Chinese Economy
    末次盛冰期以來巴倫支海-喀拉海古海洋環(huán)境及海冰研究進展
    海洋通報(2021年3期)2021-08-14 02:20:38
    不同參數(shù)對開圓洞鋼板剪力墻抗震性能的影響
    D90:全尺寸硬派SUV
    基于SIFT-SVM的北冰洋海冰識別研究
    國際標準攻堅克難“S試件”美玉漸成
    佳石選賞
    中華奇石(2015年5期)2015-07-09 18:31:07
    應用MODIS數(shù)據(jù)監(jiān)測河北省近海海域海冰
    河北遙感(2014年4期)2014-07-10 13:54:59
    国产黄色视频一区二区在线观看| 国产精品秋霞免费鲁丝片| 女性生殖器流出的白浆| a级片在线免费高清观看视频| 中国三级夫妇交换| 日韩制服骚丝袜av| 成人综合一区亚洲| 蜜桃久久精品国产亚洲av| 女的被弄到高潮叫床怎么办| 久久久久久久大尺度免费视频| 久久久久久久国产电影| 国产日韩欧美视频二区| 免费大片18禁| 欧美精品一区二区免费开放| 午夜福利网站1000一区二区三区| 色5月婷婷丁香| 五月开心婷婷网| 夜夜看夜夜爽夜夜摸| 亚洲欧洲国产日韩| 啦啦啦视频在线资源免费观看| 国产黄色视频一区二区在线观看| 日韩成人伦理影院| 最新的欧美精品一区二区| 一级片'在线观看视频| 午夜福利影视在线免费观看| 丝袜喷水一区| 欧美区成人在线视频| 国产精品福利在线免费观看| www.色视频.com| 中国三级夫妇交换| 亚洲在久久综合| 国产精品人妻久久久影院| 91久久精品国产一区二区三区| 日韩中文字幕视频在线看片| 激情五月婷婷亚洲| 一区二区三区免费毛片| a级毛色黄片| 中文字幕制服av| 极品教师在线视频| 精品久久久噜噜| 日本欧美视频一区| 高清视频免费观看一区二区| 一个人看视频在线观看www免费| 伦精品一区二区三区| 蜜桃久久精品国产亚洲av| 亚洲欧美日韩东京热| 蜜臀久久99精品久久宅男| 久久久国产精品麻豆| 天堂中文最新版在线下载| 午夜日本视频在线| 九九在线视频观看精品| 十八禁网站网址无遮挡 | 日本黄大片高清| 91久久精品电影网| 久久99蜜桃精品久久| 国产一区二区三区av在线| 大码成人一级视频| 久久久久久久精品精品| 性色avwww在线观看| 一级毛片久久久久久久久女| 乱系列少妇在线播放| 国产69精品久久久久777片| 国产精品蜜桃在线观看| 丁香六月天网| 伦精品一区二区三区| 欧美人与善性xxx| 中文天堂在线官网| 日韩在线高清观看一区二区三区| 国产精品一区二区三区四区免费观看| 寂寞人妻少妇视频99o| 男女边摸边吃奶| 久久精品久久久久久噜噜老黄| www.色视频.com| 成人影院久久| 一区二区三区免费毛片| 99热这里只有是精品在线观看| 九九爱精品视频在线观看| 亚洲av在线观看美女高潮| 99re6热这里在线精品视频| 国产精品麻豆人妻色哟哟久久| 18+在线观看网站| 日本免费在线观看一区| 国产av码专区亚洲av| 婷婷色综合大香蕉| 久久99精品国语久久久| 国产精品秋霞免费鲁丝片| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 国产精品久久久久成人av| 亚洲精品第二区| 97在线人人人人妻| 国产一区二区三区综合在线观看 | 九色成人免费人妻av| 一区二区av电影网| 日本午夜av视频| 久久久久精品久久久久真实原创| 久久人人爽人人爽人人片va| av.在线天堂| 涩涩av久久男人的天堂| 免费黄网站久久成人精品| 又爽又黄a免费视频| 中文在线观看免费www的网站| 亚洲欧美日韩另类电影网站| 精品一区二区三区视频在线| 亚洲怡红院男人天堂| 亚洲av中文av极速乱| 国产午夜精品久久久久久一区二区三区| 黄色一级大片看看| 九草在线视频观看| 老司机影院毛片| 伦精品一区二区三区| 日韩欧美一区视频在线观看 | 久久免费观看电影| 啦啦啦啦在线视频资源| 晚上一个人看的免费电影| 在线观看人妻少妇| 亚洲av在线观看美女高潮| 日日摸夜夜添夜夜添av毛片| 人人妻人人澡人人爽人人夜夜| 久久久久久久久久久久大奶| 少妇人妻一区二区三区视频| 久久99精品国语久久久| 亚洲av.av天堂| 在线观看三级黄色| 久久综合国产亚洲精品| 精品久久久久久久久av| 大码成人一级视频| 国产黄片视频在线免费观看| 国产综合精华液| 国产成人aa在线观看| 国产日韩欧美亚洲二区| 人人妻人人看人人澡| av线在线观看网站| 久久韩国三级中文字幕| 国产淫片久久久久久久久| 天美传媒精品一区二区| 青青草视频在线视频观看| 91在线精品国自产拍蜜月| 狂野欧美激情性bbbbbb| 在线观看www视频免费| 午夜福利影视在线免费观看| 久久国内精品自在自线图片| 色婷婷久久久亚洲欧美| 亚洲精品成人av观看孕妇| 99久久精品国产国产毛片| 日本黄色片子视频| 欧美丝袜亚洲另类| 99热这里只有精品一区| 日韩精品免费视频一区二区三区 | 欧美性感艳星| 91久久精品国产一区二区三区| 观看美女的网站| av国产久精品久网站免费入址| 亚洲精品456在线播放app| av线在线观看网站| 午夜免费男女啪啪视频观看| 蜜桃在线观看..| 成人亚洲精品一区在线观看| 汤姆久久久久久久影院中文字幕| 三级国产精品欧美在线观看| 亚洲,欧美,日韩| 日韩精品免费视频一区二区三区 | 一级,二级,三级黄色视频| 少妇人妻一区二区三区视频| 女人精品久久久久毛片| 一级,二级,三级黄色视频| 国产伦精品一区二区三区视频9| 久久人人爽人人爽人人片va| 夫妻午夜视频| 国产精品女同一区二区软件| 我要看日韩黄色一级片| 午夜91福利影院| 免费观看的影片在线观看| 国产老妇伦熟女老妇高清| 国产精品伦人一区二区| 免费在线观看成人毛片| 国产一级毛片在线| 国产色爽女视频免费观看| 中文字幕久久专区| 国产伦精品一区二区三区四那| 中国美白少妇内射xxxbb| 久久女婷五月综合色啪小说| 免费观看a级毛片全部| 欧美最新免费一区二区三区| 涩涩av久久男人的天堂| 少妇人妻精品综合一区二区| 中文乱码字字幕精品一区二区三区| 热99国产精品久久久久久7| 男人添女人高潮全过程视频| 亚洲av中文av极速乱| 中国国产av一级| 亚洲婷婷狠狠爱综合网| 国产精品伦人一区二区| 日韩视频在线欧美| 欧美+日韩+精品| 亚洲精品日韩在线中文字幕| 成人国产av品久久久| 免费黄频网站在线观看国产| 看十八女毛片水多多多| 国产欧美日韩综合在线一区二区 | 女的被弄到高潮叫床怎么办| a级片在线免费高清观看视频| 国产亚洲5aaaaa淫片| 亚洲欧洲日产国产| 中文乱码字字幕精品一区二区三区| 免费人成在线观看视频色| 五月伊人婷婷丁香| 日韩亚洲欧美综合| 久久久久久久久久成人| 2018国产大陆天天弄谢| 国产日韩欧美亚洲二区| 久久精品国产自在天天线| 丰满饥渴人妻一区二区三| 国产成人精品福利久久| 97在线视频观看| 超碰97精品在线观看| 夫妻性生交免费视频一级片| √禁漫天堂资源中文www| 18禁动态无遮挡网站| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 亚洲激情五月婷婷啪啪| 亚洲高清免费不卡视频| 熟女av电影| 99久久中文字幕三级久久日本| 熟女人妻精品中文字幕| 大香蕉久久网| 久久午夜综合久久蜜桃| 亚洲丝袜综合中文字幕| 草草在线视频免费看| 少妇 在线观看| 欧美最新免费一区二区三区| 亚洲国产精品一区三区| 亚洲三级黄色毛片| 亚洲av.av天堂| 欧美丝袜亚洲另类| 在线看a的网站| 国产成人免费无遮挡视频| 免费人成在线观看视频色| 日韩精品有码人妻一区| 日韩 亚洲 欧美在线| 波野结衣二区三区在线| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 亚洲国产精品专区欧美| 成人二区视频| 亚洲国产精品一区二区三区在线| 一区二区三区乱码不卡18| 国产成人午夜福利电影在线观看| 97超视频在线观看视频| 国产欧美日韩一区二区三区在线 | 99久久中文字幕三级久久日本| 欧美日韩av久久| 国产精品一区二区在线观看99| 狂野欧美激情性xxxx在线观看| 久久久久久久国产电影| 亚洲真实伦在线观看| 一本一本综合久久| 欧美日韩综合久久久久久| 久久久久人妻精品一区果冻| a级一级毛片免费在线观看| 日韩不卡一区二区三区视频在线| 久久99热6这里只有精品| 国产精品久久久久久久电影| 王馨瑶露胸无遮挡在线观看| 免费看光身美女| 久久久精品94久久精品| 又爽又黄a免费视频| 国产一区二区三区av在线| 国内精品宾馆在线| 成人免费观看视频高清| 91在线精品国自产拍蜜月| 另类亚洲欧美激情| 精品一区二区免费观看| 在现免费观看毛片| 成人国产麻豆网| 中文精品一卡2卡3卡4更新| 欧美精品亚洲一区二区| 午夜视频国产福利| 成人毛片60女人毛片免费| 一级毛片 在线播放| 久久6这里有精品| 婷婷色av中文字幕| h视频一区二区三区| 国产精品99久久久久久久久| 王馨瑶露胸无遮挡在线观看| 如何舔出高潮| 大陆偷拍与自拍| 大片电影免费在线观看免费| 最近2019中文字幕mv第一页| 久久久欧美国产精品| 日产精品乱码卡一卡2卡三| 国产精品女同一区二区软件| 欧美成人精品欧美一级黄| 国产一区二区三区综合在线观看 | 久久精品国产亚洲av涩爱| 九色成人免费人妻av| 91成人精品电影| 少妇丰满av| 另类精品久久| 久久99蜜桃精品久久| 成年人午夜在线观看视频| 国产男女超爽视频在线观看| 少妇猛男粗大的猛烈进出视频| 欧美老熟妇乱子伦牲交| 我要看日韩黄色一级片| 亚洲欧美成人精品一区二区| 97在线视频观看| 国产亚洲av片在线观看秒播厂| 欧美bdsm另类| 在线观看www视频免费| 男女啪啪激烈高潮av片| 精品久久久精品久久久| h日本视频在线播放| 青青草视频在线视频观看| 男女无遮挡免费网站观看| 美女国产视频在线观看| 久久久久久人妻| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 久久99蜜桃精品久久| 精品一区二区三区视频在线| 国产免费福利视频在线观看| 亚洲成人手机| 亚洲av二区三区四区| av在线播放精品| 国产午夜精品一二区理论片| 色哟哟·www| 在线观看一区二区三区激情| 亚州av有码| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 男人狂女人下面高潮的视频| 69精品国产乱码久久久| 国产91av在线免费观看| 国产熟女午夜一区二区三区 | 十八禁网站网址无遮挡 | 人体艺术视频欧美日本| 午夜福利视频精品| 老司机影院毛片| 亚洲色图综合在线观看| 一个人看视频在线观看www免费| 18禁在线播放成人免费| 人人妻人人爽人人添夜夜欢视频 | 在线观看免费视频网站a站| 97超视频在线观看视频| 国产色婷婷99| 国产亚洲午夜精品一区二区久久| 亚洲人成网站在线观看播放| 免费久久久久久久精品成人欧美视频 | 精品熟女少妇av免费看| 啦啦啦中文免费视频观看日本| 亚洲国产毛片av蜜桃av| 3wmmmm亚洲av在线观看| 欧美精品亚洲一区二区| 少妇高潮的动态图| 99热网站在线观看| 亚洲欧美中文字幕日韩二区| 午夜福利,免费看| www.色视频.com| 国产精品国产三级国产专区5o| 成人亚洲欧美一区二区av| 多毛熟女@视频| 国产91av在线免费观看| 亚洲av国产av综合av卡| 久久韩国三级中文字幕| 精品一区二区三卡| 久久精品国产自在天天线| 亚洲欧美一区二区三区国产| 九草在线视频观看| 亚洲精品乱久久久久久| 秋霞在线观看毛片| 亚洲无线观看免费| 草草在线视频免费看| 国产69精品久久久久777片| 欧美区成人在线视频| 亚洲情色 制服丝袜| av在线老鸭窝| 91aial.com中文字幕在线观看| 午夜91福利影院| 欧美人与善性xxx| 一级毛片 在线播放| 国产精品伦人一区二区| 欧美精品一区二区大全| 人人妻人人添人人爽欧美一区卜| tube8黄色片| 久久国产乱子免费精品| 欧美精品亚洲一区二区| 熟女电影av网| 黄色怎么调成土黄色| 久久久精品免费免费高清| 精品久久久久久电影网| 韩国av在线不卡| 色网站视频免费| 久久久久人妻精品一区果冻| 亚洲国产毛片av蜜桃av| 国产真实伦视频高清在线观看| 在现免费观看毛片| 在线天堂最新版资源| 久久国产亚洲av麻豆专区| 国产在线男女| 少妇裸体淫交视频免费看高清| 伊人久久精品亚洲午夜| 王馨瑶露胸无遮挡在线观看| 寂寞人妻少妇视频99o| 国内精品宾馆在线| 国产女主播在线喷水免费视频网站| 欧美日韩av久久| 国产黄色免费在线视频| 极品教师在线视频| 国产精品秋霞免费鲁丝片| 涩涩av久久男人的天堂| 国产精品免费大片| 中文资源天堂在线| 久久6这里有精品| 人人妻人人澡人人爽人人夜夜| 国语对白做爰xxxⅹ性视频网站| 插逼视频在线观看| 欧美 亚洲 国产 日韩一| 天天操日日干夜夜撸| 久久久久国产精品人妻一区二区| 日本与韩国留学比较| 国产精品蜜桃在线观看| 精品国产一区二区久久| 99热这里只有是精品50| 91久久精品国产一区二区三区| 国产精品人妻久久久久久| 亚洲国产精品一区三区| 看十八女毛片水多多多| 男人和女人高潮做爰伦理| av网站免费在线观看视频| 国产乱来视频区| 在线观看免费日韩欧美大片 | 国产成人一区二区在线| 国产黄片视频在线免费观看| 色94色欧美一区二区| 高清欧美精品videossex| 三级国产精品欧美在线观看| 免费黄色在线免费观看| 国产一区二区三区av在线| 国产精品久久久久久久久免| 搡女人真爽免费视频火全软件| 亚洲人成网站在线播| 女人精品久久久久毛片| 偷拍熟女少妇极品色| 亚洲久久久国产精品| 日日爽夜夜爽网站| 欧美三级亚洲精品| 免费在线观看成人毛片| 又大又黄又爽视频免费| tube8黄色片| 一级a做视频免费观看| 97超碰精品成人国产| 国产高清国产精品国产三级| 女人久久www免费人成看片| 日韩免费高清中文字幕av| 纵有疾风起免费观看全集完整版| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 精品亚洲成a人片在线观看| 少妇的逼好多水| 熟女av电影| 少妇猛男粗大的猛烈进出视频| 日本免费在线观看一区| 国产 精品1| 国产精品成人在线| 亚洲国产精品国产精品| 欧美精品国产亚洲| 久久久久久久久久成人| 一级毛片我不卡| av.在线天堂| 下体分泌物呈黄色| 国产 一区精品| 欧美成人午夜免费资源| 亚洲精品色激情综合| 国产精品人妻久久久久久| 日日摸夜夜添夜夜爱| 黄色一级大片看看| 91久久精品国产一区二区成人| 九草在线视频观看| 夜夜骑夜夜射夜夜干| 成年女人在线观看亚洲视频| 日韩不卡一区二区三区视频在线| 国产精品伦人一区二区| 中文欧美无线码| 中文乱码字字幕精品一区二区三区| 免费观看av网站的网址| 婷婷色综合www| 精品国产一区二区三区久久久樱花| 六月丁香七月| 五月开心婷婷网| 少妇猛男粗大的猛烈进出视频| 成年人免费黄色播放视频 | 超碰97精品在线观看| 亚洲av中文av极速乱| 亚洲国产精品成人久久小说| 高清午夜精品一区二区三区| 免费观看无遮挡的男女| 大香蕉久久网| 国产黄片视频在线免费观看| 蜜桃在线观看..| av一本久久久久| 日韩电影二区| 97在线视频观看| 精品一品国产午夜福利视频| 有码 亚洲区| 国产伦精品一区二区三区四那| 色哟哟·www| 久久毛片免费看一区二区三区| 国产亚洲最大av| 国产午夜精品一二区理论片| av国产精品久久久久影院| 午夜福利在线观看免费完整高清在| 97精品久久久久久久久久精品| 日韩中字成人| 日韩视频在线欧美| 久久综合国产亚洲精品| 老司机影院毛片| 亚洲欧美中文字幕日韩二区| 深夜a级毛片| 国产精品麻豆人妻色哟哟久久| 偷拍熟女少妇极品色| 熟女电影av网| 国产69精品久久久久777片| 最近最新中文字幕免费大全7| av视频免费观看在线观看| 久久6这里有精品| 永久免费av网站大全| 热re99久久精品国产66热6| 18禁裸乳无遮挡动漫免费视频| 午夜福利,免费看| 好男人视频免费观看在线| 色婷婷久久久亚洲欧美| 2021少妇久久久久久久久久久| 成人免费观看视频高清| 赤兔流量卡办理| 男女国产视频网站| av视频免费观看在线观看| 99九九线精品视频在线观看视频| 91精品一卡2卡3卡4卡| a级毛片免费高清观看在线播放| 又黄又爽又刺激的免费视频.| 99热国产这里只有精品6| 乱人伦中国视频| 国产伦精品一区二区三区视频9| 精品亚洲乱码少妇综合久久| 亚洲国产精品专区欧美| 精品亚洲乱码少妇综合久久| 日韩三级伦理在线观看| 亚洲国产毛片av蜜桃av| 一边亲一边摸免费视频| 精品一区二区三卡| 97在线视频观看| 亚洲人成网站在线观看播放| 亚洲欧洲精品一区二区精品久久久 | 内射极品少妇av片p| 观看美女的网站| 女性生殖器流出的白浆| 人妻人人澡人人爽人人| 国产精品一区二区性色av| 国产精品久久久久久精品电影小说| 男女边摸边吃奶| 妹子高潮喷水视频| 国产真实伦视频高清在线观看| 亚洲av综合色区一区| 精华霜和精华液先用哪个| 欧美97在线视频| 观看美女的网站| 边亲边吃奶的免费视频| 纯流量卡能插随身wifi吗| 欧美另类一区| 美女脱内裤让男人舔精品视频| 人妻系列 视频| 日日撸夜夜添| 蜜桃在线观看..| 国产亚洲精品久久久com| 日韩中字成人| 午夜免费男女啪啪视频观看| 国产免费一区二区三区四区乱码| 在线观看免费日韩欧美大片 | 伊人亚洲综合成人网| 伊人久久国产一区二区| 国产综合精华液| 亚洲精品aⅴ在线观看| 如何舔出高潮| 男女边吃奶边做爰视频| 亚洲经典国产精华液单| 久久久久国产精品人妻一区二区| 丰满人妻一区二区三区视频av| 亚洲国产日韩一区二区| 免费黄色在线免费观看| 99热网站在线观看| 91精品一卡2卡3卡4卡| 国产在视频线精品| 欧美日韩在线观看h| 又黄又爽又刺激的免费视频.| 国精品久久久久久国模美| 午夜日本视频在线| 涩涩av久久男人的天堂| 亚洲色图综合在线观看| freevideosex欧美| 国产日韩欧美在线精品| 亚洲色图综合在线观看| 婷婷色麻豆天堂久久| 欧美一级a爱片免费观看看| 大片免费播放器 马上看| 亚洲美女黄色视频免费看| 日韩一区二区视频免费看| 久久韩国三级中文字幕| 亚洲精品视频女| 精品久久久噜噜| 免费观看无遮挡的男女| 免费人妻精品一区二区三区视频| 日韩电影二区|