• <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
    毛片一级片免费看久久久久 | 男女做爰动态图高潮gif福利片| 国产黄片美女视频| 一本精品99久久精品77| 狂野欧美白嫩少妇大欣赏| 99国产精品一区二区蜜桃av| 淫秽高清视频在线观看| 色综合亚洲欧美另类图片| 国产男靠女视频免费网站| 999久久久精品免费观看国产| 国产精品乱码一区二三区的特点| 国产aⅴ精品一区二区三区波| 欧美中文日本在线观看视频| 成年人黄色毛片网站| 欧美成人a在线观看| 国产黄色小视频在线观看| 日日撸夜夜添| 国产一区二区在线观看日韩| 国产高潮美女av| 婷婷精品国产亚洲av| 中文资源天堂在线| 中文字幕久久专区| 级片在线观看| 精品日产1卡2卡| 免费在线观看成人毛片| 老女人水多毛片| 日日干狠狠操夜夜爽| 免费av毛片视频| 亚洲成人中文字幕在线播放| 国产一区二区三区视频了| 在线a可以看的网站| 婷婷亚洲欧美| 亚洲成av人片在线播放无| 免费观看精品视频网站| 国产精品久久电影中文字幕| 免费看日本二区| 国产日本99.免费观看| 精品人妻熟女av久视频| 精品欧美国产一区二区三| 午夜免费激情av| 久久久久九九精品影院| 日韩欧美精品免费久久| 嫩草影院新地址| 美女xxoo啪啪120秒动态图| 精品久久国产蜜桃| 精品一区二区免费观看| eeuss影院久久| 色在线成人网| 人人妻人人澡欧美一区二区| 国产一级毛片七仙女欲春2| 女同久久另类99精品国产91| 波野结衣二区三区在线| 少妇熟女aⅴ在线视频| 在线观看舔阴道视频| 色综合婷婷激情| 亚洲三级黄色毛片| 乱码一卡2卡4卡精品| 高清在线国产一区| 精品99又大又爽又粗少妇毛片 | 色哟哟哟哟哟哟| 日本成人三级电影网站| 黄色日韩在线| 亚洲经典国产精华液单| 亚洲色图av天堂| 国产探花极品一区二区| 日韩国内少妇激情av| 欧美日本视频| 亚洲午夜理论影院| 成人精品一区二区免费| 亚洲在线观看片| 亚洲人成网站高清观看| 十八禁国产超污无遮挡网站| 国产在线精品亚洲第一网站| 男女那种视频在线观看| 真实男女啪啪啪动态图| 在线播放国产精品三级| 日韩欧美在线二视频| 国产久久久一区二区三区| 97超级碰碰碰精品色视频在线观看| 美女 人体艺术 gogo| 能在线免费观看的黄片| 久9热在线精品视频| 最近在线观看免费完整版| 日韩亚洲欧美综合| 一进一出抽搐动态| 精品无人区乱码1区二区| 日韩 亚洲 欧美在线| 成人性生交大片免费视频hd| av专区在线播放| 老司机福利观看| 如何舔出高潮| 特大巨黑吊av在线直播| 啦啦啦啦在线视频资源| aaaaa片日本免费| 女人被狂操c到高潮| 如何舔出高潮| 国产精品久久电影中文字幕| 别揉我奶头 嗯啊视频| 99热这里只有精品一区| 别揉我奶头 嗯啊视频| 尤物成人国产欧美一区二区三区| 99久久成人亚洲精品观看| 别揉我奶头 嗯啊视频| 亚洲avbb在线观看| 最近最新中文字幕大全电影3| 亚洲不卡免费看| 色尼玛亚洲综合影院| 国内久久婷婷六月综合欲色啪| 亚洲第一电影网av| 老女人水多毛片| 99精品在免费线老司机午夜| a级一级毛片免费在线观看| 真人做人爱边吃奶动态| 人人妻人人看人人澡| 看黄色毛片网站| 国产成人一区二区在线| 又黄又爽又免费观看的视频| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 日韩欧美免费精品| 亚洲国产精品成人综合色| 男女视频在线观看网站免费| 日韩高清综合在线| 亚洲国产色片| 在线免费观看不下载黄p国产 | av天堂中文字幕网| 女同久久另类99精品国产91| 99久久九九国产精品国产免费| 国产探花在线观看一区二区| 免费人成在线观看视频色| 国产又黄又爽又无遮挡在线| 熟女人妻精品中文字幕| 欧美成人一区二区免费高清观看| 制服丝袜大香蕉在线| 色播亚洲综合网| 欧美xxxx黑人xx丫x性爽| 禁无遮挡网站| 人妻少妇偷人精品九色| 美女 人体艺术 gogo| 国产69精品久久久久777片| 亚洲精品456在线播放app | 日本免费一区二区三区高清不卡| 日韩人妻高清精品专区| 波多野结衣高清无吗| 男女下面进入的视频免费午夜| 亚洲av美国av| 成年女人永久免费观看视频| 99九九线精品视频在线观看视频| 免费大片18禁| 国产高清有码在线观看视频| 午夜福利视频1000在线观看| 久99久视频精品免费| 欧美中文日本在线观看视频| 直男gayav资源| 最近最新中文字幕大全电影3| 免费人成在线观看视频色| 在线国产一区二区在线| 成人高潮视频无遮挡免费网站| 精品不卡国产一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 精品一区二区三区人妻视频| 久久久久久九九精品二区国产| 亚洲真实伦在线观看| avwww免费| 亚洲av成人av| 校园人妻丝袜中文字幕| 精品久久久久久久末码| 久久九九热精品免费| x7x7x7水蜜桃| 国产精品永久免费网站| 大又大粗又爽又黄少妇毛片口| 成人精品一区二区免费| 欧美激情久久久久久爽电影| 色哟哟哟哟哟哟| 男女啪啪激烈高潮av片| 亚洲精品日韩av片在线观看| 婷婷六月久久综合丁香| 亚洲精品国产成人久久av| av视频在线观看入口| 午夜免费激情av| xxxwww97欧美| 天堂影院成人在线观看| 免费人成在线观看视频色| 国产精品一区二区三区四区免费观看 | 乱人视频在线观看| 免费搜索国产男女视频| 久久久久九九精品影院| 中文字幕精品亚洲无线码一区| 亚洲不卡免费看| 日日摸夜夜添夜夜添小说| 黄色一级大片看看| ponron亚洲| 精品一区二区免费观看| 一进一出抽搐gif免费好疼| 伊人久久精品亚洲午夜| av.在线天堂| 人人妻人人澡欧美一区二区| 国产精品一及| 综合色av麻豆| 天天躁日日操中文字幕| 在线免费十八禁| 欧美另类亚洲清纯唯美| 联通29元200g的流量卡| 老熟妇乱子伦视频在线观看| 欧美一区二区精品小视频在线| 美女免费视频网站| 亚洲熟妇中文字幕五十中出| 嫩草影院新地址| 一本一本综合久久| 午夜福利成人在线免费观看| 大又大粗又爽又黄少妇毛片口| 欧美+日韩+精品| 国产高清不卡午夜福利| 91午夜精品亚洲一区二区三区 | 日韩亚洲欧美综合| 日本五十路高清| 自拍偷自拍亚洲精品老妇| 九九久久精品国产亚洲av麻豆| 成人特级av手机在线观看| 一区二区三区免费毛片| aaaaa片日本免费| 真实男女啪啪啪动态图| 色5月婷婷丁香| 精品不卡国产一区二区三区| 最好的美女福利视频网| 久久久久久久久大av| 国产黄色小视频在线观看| 日韩欧美精品免费久久| 久久草成人影院| 免费在线观看日本一区| 国产黄片美女视频| 精品人妻1区二区| 99热只有精品国产| 狠狠狠狠99中文字幕| 禁无遮挡网站| 成熟少妇高潮喷水视频| 国产av麻豆久久久久久久| 婷婷亚洲欧美| 欧美日本亚洲视频在线播放| 在线免费观看的www视频| 偷拍熟女少妇极品色| 精品欧美国产一区二区三| 蜜桃久久精品国产亚洲av| 搡老熟女国产l中国老女人| 永久网站在线| 亚洲综合色惰| 亚洲狠狠婷婷综合久久图片| 欧美色视频一区免费| 欧美中文日本在线观看视频| 国产激情偷乱视频一区二区| 免费无遮挡裸体视频| 免费看a级黄色片| 国产精品日韩av在线免费观看| 欧美成人一区二区免费高清观看| 欧美丝袜亚洲另类 | 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 成年人黄色毛片网站| 国产在视频线在精品| 伦精品一区二区三区| 日韩欧美 国产精品| 丰满人妻一区二区三区视频av| 亚洲七黄色美女视频| 一级黄片播放器| 免费在线观看日本一区| av天堂中文字幕网| www日本黄色视频网| 丰满乱子伦码专区| 国产欧美日韩一区二区精品| 尤物成人国产欧美一区二区三区| 欧美激情久久久久久爽电影| av天堂在线播放| 免费高清视频大片| 亚洲国产精品成人综合色| 男人狂女人下面高潮的视频| 麻豆一二三区av精品| 成人特级黄色片久久久久久久| 午夜a级毛片| 69人妻影院| 国产亚洲欧美98| 日日摸夜夜添夜夜添小说| 日本黄大片高清| 国产精品av视频在线免费观看| 观看美女的网站| 国内精品久久久久久久电影| 免费一级毛片在线播放高清视频| 久久香蕉精品热| 久久精品久久久久久噜噜老黄 | 丰满乱子伦码专区| 国产成人影院久久av| 精品99又大又爽又粗少妇毛片 | 99在线人妻在线中文字幕| 亚洲最大成人手机在线| 精品人妻一区二区三区麻豆 | 桃色一区二区三区在线观看| 中文字幕av在线有码专区| 国产免费一级a男人的天堂| 国产伦精品一区二区三区视频9| 欧美在线一区亚洲| 99久国产av精品| 久久久久久久久久成人| 在线观看av片永久免费下载| 欧美精品国产亚洲| 日韩高清综合在线| 久久久久久久久久久丰满 | 色在线成人网| 欧美一区二区亚洲| 久久久久性生活片| 日日摸夜夜添夜夜添av毛片 | 天天躁日日操中文字幕| 99热这里只有精品一区| 少妇高潮的动态图| 自拍偷自拍亚洲精品老妇| 黄色视频,在线免费观看| 亚洲人与动物交配视频| 日日夜夜操网爽| 老熟妇仑乱视频hdxx| 国产av不卡久久| 午夜精品在线福利| 亚洲人成网站高清观看| 成年女人看的毛片在线观看| 在线免费观看的www视频| 午夜免费激情av| 国产成人a区在线观看| 国产午夜福利久久久久久| 婷婷六月久久综合丁香| 嫩草影院入口| 高清在线国产一区| 国产精品亚洲一级av第二区| 国产私拍福利视频在线观看| 嫁个100分男人电影在线观看| 国产精品久久视频播放| 欧美丝袜亚洲另类 | 免费av毛片视频| 欧美黑人巨大hd| 欧美国产日韩亚洲一区| 白带黄色成豆腐渣| 一区福利在线观看| 啪啪无遮挡十八禁网站| 精品久久久久久成人av| 看十八女毛片水多多多| 99久久精品热视频| 午夜老司机福利剧场| 真实男女啪啪啪动态图| 成年女人看的毛片在线观看| 日本精品一区二区三区蜜桃| 久久久国产成人免费| 国产av麻豆久久久久久久| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲| 亚洲成a人片在线一区二区| 亚洲,欧美,日韩| 美女大奶头视频| 欧美zozozo另类| 欧美高清性xxxxhd video| av天堂中文字幕网| 丰满人妻一区二区三区视频av| 毛片女人毛片| 日韩大尺度精品在线看网址| 我的老师免费观看完整版| 免费高清视频大片| 国产精品美女特级片免费视频播放器| 成人精品一区二区免费| 国产男靠女视频免费网站| 超碰av人人做人人爽久久| 国产69精品久久久久777片| 伊人久久精品亚洲午夜| 国产不卡一卡二| 精品久久久久久久末码| 一本精品99久久精品77| 国产午夜福利久久久久久| 国产av不卡久久| 51国产日韩欧美| 欧美绝顶高潮抽搐喷水| 日韩中字成人| 亚洲国产精品成人综合色| 乱码一卡2卡4卡精品| 成年女人看的毛片在线观看| 国产一区二区在线av高清观看| 国产人妻一区二区三区在| 尾随美女入室| 久久精品久久久久久噜噜老黄 | 亚洲七黄色美女视频| 又爽又黄a免费视频| 1000部很黄的大片| 露出奶头的视频| 欧美丝袜亚洲另类 | 乱人视频在线观看| 女的被弄到高潮叫床怎么办 | 免费观看人在逋| 成人亚洲精品av一区二区| 亚洲中文字幕日韩| 99在线人妻在线中文字幕| 久久久国产成人精品二区| 亚洲av二区三区四区| 亚洲综合色惰| 成人性生交大片免费视频hd| 中国美白少妇内射xxxbb| 国产一区二区在线观看日韩| 日韩高清综合在线| av专区在线播放| 国产精品福利在线免费观看| 久久久久久久精品吃奶| 一本精品99久久精品77| 12—13女人毛片做爰片一| 一夜夜www| 男女那种视频在线观看| 精品人妻偷拍中文字幕| 欧美一区二区亚洲| 别揉我奶头 嗯啊视频| 成人av一区二区三区在线看| 极品教师在线免费播放| 校园春色视频在线观看| 亚洲国产高清在线一区二区三| 亚洲av日韩精品久久久久久密| 成人鲁丝片一二三区免费| 国产黄片美女视频| 亚洲av美国av| 色视频www国产| 国产视频一区二区在线看| 国产国拍精品亚洲av在线观看| 成人二区视频| 国产精品久久久久久久电影| 国产一区二区激情短视频| 特大巨黑吊av在线直播| 欧美三级亚洲精品| 夜夜夜夜夜久久久久| 欧美丝袜亚洲另类 | 日日撸夜夜添| 中文字幕精品亚洲无线码一区| 色综合色国产| 精品午夜福利在线看| aaaaa片日本免费| 亚洲经典国产精华液单| 精品久久国产蜜桃| 男人狂女人下面高潮的视频| 日本成人三级电影网站| a级毛片a级免费在线| 中文字幕免费在线视频6| 日本黄大片高清| 制服丝袜大香蕉在线| 九色国产91popny在线| 嫩草影视91久久| 亚洲aⅴ乱码一区二区在线播放| 久久久精品大字幕| 亚洲国产精品sss在线观看| 午夜精品在线福利| 国产伦在线观看视频一区| 国产精品,欧美在线| 国产成人av教育| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 女的被弄到高潮叫床怎么办 | 男女视频在线观看网站免费| 天堂av国产一区二区熟女人妻| 18禁黄网站禁片免费观看直播| 给我免费播放毛片高清在线观看| 欧美+日韩+精品| 色吧在线观看| 我要搜黄色片| 亚洲图色成人| 亚洲精品亚洲一区二区| 久久天躁狠狠躁夜夜2o2o| 国产在线男女| 草草在线视频免费看| 最近中文字幕高清免费大全6 | 国产av一区在线观看免费| 国产白丝娇喘喷水9色精品| 国产精品三级大全| 国产精品久久电影中文字幕| 久久精品人妻少妇| 中亚洲国语对白在线视频| 级片在线观看| 午夜福利在线观看免费完整高清在 | 亚洲黑人精品在线| 成人一区二区视频在线观看| 18禁黄网站禁片免费观看直播| 亚洲av免费高清在线观看| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 亚洲欧美精品综合久久99| 三级国产精品欧美在线观看| 日本爱情动作片www.在线观看 | 国产精品永久免费网站| 免费在线观看影片大全网站| 最近最新免费中文字幕在线| 日韩欧美在线乱码| 亚洲av美国av| 91精品国产九色| 一区二区三区高清视频在线| 欧美黑人巨大hd| 日韩在线高清观看一区二区三区 | 免费搜索国产男女视频| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 韩国av在线不卡| 国产极品精品免费视频能看的| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 免费av毛片视频| 女人被狂操c到高潮| 精品人妻偷拍中文字幕| 亚洲va在线va天堂va国产| 国产精品乱码一区二三区的特点| 国产黄a三级三级三级人| 亚洲美女搞黄在线观看 | 亚洲电影在线观看av| 超碰av人人做人人爽久久| 成人欧美大片| 久久久久久久久大av| 欧美中文日本在线观看视频| 成人毛片a级毛片在线播放| 性插视频无遮挡在线免费观看| 日日摸夜夜添夜夜添av毛片 | 成人亚洲精品av一区二区| 亚洲 国产 在线| 真人一进一出gif抽搐免费| 欧美不卡视频在线免费观看| 日韩强制内射视频| 免费人成在线观看视频色| 九九久久精品国产亚洲av麻豆| 国产视频内射| 男女啪啪激烈高潮av片| 毛片女人毛片| 国产白丝娇喘喷水9色精品| 校园人妻丝袜中文字幕| 国产亚洲精品av在线| 天堂√8在线中文| 又爽又黄a免费视频| 午夜a级毛片| 免费看av在线观看网站| 国产精品综合久久久久久久免费| 久久久久九九精品影院| 免费看光身美女| or卡值多少钱| 亚洲av中文字字幕乱码综合| 99久久中文字幕三级久久日本| 婷婷丁香在线五月| 欧美又色又爽又黄视频| av女优亚洲男人天堂| 在线观看免费视频日本深夜| 亚洲欧美日韩卡通动漫| 成人永久免费在线观看视频| 有码 亚洲区| 男女那种视频在线观看| 久久久久久伊人网av| 精品久久久久久成人av| 亚洲精品日韩av片在线观看| 人妻少妇偷人精品九色| 一区二区三区激情视频| 亚洲国产精品成人综合色| 国产伦在线观看视频一区| a级毛片a级免费在线| 看片在线看免费视频| 九色成人免费人妻av| 一本精品99久久精品77| 草草在线视频免费看| 欧美成人免费av一区二区三区| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| 成人国产一区最新在线观看| 麻豆成人午夜福利视频| 久久久久久九九精品二区国产| 久久国内精品自在自线图片| 小说图片视频综合网站| 啪啪无遮挡十八禁网站| 色综合色国产| 亚洲七黄色美女视频| 亚洲,欧美,日韩| 超碰av人人做人人爽久久| 亚洲成av人片在线播放无| 国产私拍福利视频在线观看| 国产在视频线在精品| 99热网站在线观看| 亚洲不卡免费看| 草草在线视频免费看| 亚洲国产精品久久男人天堂| 国产主播在线观看一区二区| 欧美高清性xxxxhd video| 亚洲精品成人久久久久久| 91精品国产九色| 精品久久久久久久久亚洲 | 日本黄大片高清| 给我免费播放毛片高清在线观看| 中文亚洲av片在线观看爽| 日韩高清综合在线| 熟女电影av网| 99久久精品国产国产毛片| 国产一级毛片七仙女欲春2| 精品一区二区三区视频在线| 国产免费一级a男人的天堂| 国产亚洲精品久久久com| 欧美日韩黄片免| 好男人在线观看高清免费视频| 悠悠久久av| 国产精品一区www在线观看 | 桃红色精品国产亚洲av| 国产毛片a区久久久久| 日本a在线网址| 欧美日韩瑟瑟在线播放| 成人午夜高清在线视频| 国产毛片a区久久久久| 女人十人毛片免费观看3o分钟| 欧美激情久久久久久爽电影| 日韩av在线大香蕉| 国产精品98久久久久久宅男小说| 色精品久久人妻99蜜桃| 日本成人三级电影网站| 男插女下体视频免费在线播放| 国产午夜福利久久久久久| 好男人在线观看高清免费视频| 悠悠久久av| 午夜激情欧美在线| 黄色配什么色好看| 一个人免费在线观看电影|