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

    低溫極端載荷作用下船艏結(jié)構(gòu)損傷演化過(guò)程研究

    2023-10-18 02:32:28趙南汪高飛葛辛辛李飛張占陽(yáng)
    裝備環(huán)境工程 2023年9期
    關(guān)鍵詞:冰體海冰冰山

    趙南,汪高飛,葛辛辛,李飛,張占陽(yáng)

    (1.中國(guó)船舶科學(xué)研究中心,江蘇 無(wú)錫 214082;2.深海技術(shù)科學(xué)太湖實(shí)驗(yàn)室,江蘇 無(wú)錫 214082)

    北極圈油氣和各類礦產(chǎn)資源十分豐富,目前已經(jīng)成為各國(guó)索取資源的供給基地[1]。同時(shí),隨著全球氣候變暖和北極冰蓋融化進(jìn)程的加速[2],北極航道的通行窗口期將越來(lái)越長(zhǎng),人類在北極的各種商業(yè)、科學(xué)考察等活動(dòng)也將更加頻繁。中國(guó)作為鄰近北極圈的重要大國(guó)和國(guó)際貿(mào)易強(qiáng)國(guó),政治、經(jīng)濟(jì)、科技等方面都受到北極環(huán)境變化的重大影響。針對(duì)這一現(xiàn)狀,我國(guó)也在積極應(yīng)對(duì),并于2018 年1 月發(fā)布了《中國(guó)的北極政策》白皮書(shū)[3],目前已經(jīng)承接部分冰區(qū)船訂單,并對(duì)大型冰區(qū)船舶進(jìn)行了技術(shù)上的研制開(kāi)發(fā)。

    自20 世紀(jì)90 年代至21 世紀(jì)初,主要針對(duì)船-船碰撞、船-橋碰撞等開(kāi)展了大量研究。Paik 等[4]開(kāi)展了觸底狀態(tài)下船舶結(jié)構(gòu)的耐撞性研究。李江濤等[5]開(kāi)展了擱淺于剛性斜坡下船舶結(jié)構(gòu)的耐撞性研究。趙南等[6-7]進(jìn)行了補(bǔ)給作業(yè)船碰撞場(chǎng)景確定及損傷環(huán)境仿真分析。近10 余年來(lái),針對(duì)極地環(huán)境,國(guó)內(nèi)外學(xué)者也開(kāi)展了大量實(shí)驗(yàn)和仿真分析。劉俊杰等[8]開(kāi)展了油船結(jié)構(gòu)與浮冰碰撞仿真分析,獲得了油船結(jié)構(gòu)響應(yīng)。Bruce[9]、Kim[10]、Abdullah[11]、Martin[12]等開(kāi)展了冰載荷作用下的板架結(jié)構(gòu)損傷變形試驗(yàn)和仿真分析,研究了冰載荷作用下加筋板的損傷變形情況等。袁光奇[13]開(kāi)展了導(dǎo)管架平臺(tái)附近碎冰堆積過(guò)程模型試驗(yàn)研究,揭示了碎冰堆積過(guò)程機(jī)理。黃焱等[14-15]、蔡柯等[16]、Antonio[17]、Holm[18]、Marnix[19]開(kāi)展了船冰碰撞載荷模型試驗(yàn)研究以及仿真技術(shù)研究等,得到了冰載荷時(shí)空分布規(guī)律及統(tǒng)計(jì)規(guī)律等。閆巖[20]基于數(shù)值仿真方法開(kāi)展了冰載荷計(jì)算,并獲得了該載荷作用下加筋板結(jié)構(gòu)的承載能力。Bahar[21]基于IACS 規(guī)范開(kāi)展了FPSO 結(jié)構(gòu)響應(yīng)及失效機(jī)理。Abraham[22]、Mihkel[23]、Liu[24]等基于數(shù)值仿真方法開(kāi)展了船冰相互作用下船體結(jié)構(gòu)響應(yīng)分析及承載能力研究。

    總體來(lái)說(shuō),目前對(duì)于船舶在浮冰中航行的冰載荷特性、結(jié)構(gòu)響應(yīng)研究以及冰體與板架結(jié)構(gòu)相互作用的模型試驗(yàn)研究較多,大型以及實(shí)船船冰相互作用模型試驗(yàn)研究較少,且目前對(duì)于如冰山、冰脊極端載荷作用下的船舶結(jié)構(gòu)損傷演化過(guò)程的研究較少。因此,本文將針對(duì)低溫以及極端載荷作用下的船舶艏部結(jié)構(gòu)損傷演化過(guò)程進(jìn)行仿真分析,掌握各主要性能參數(shù)對(duì)損傷演化過(guò)程的影響規(guī)律,為極端環(huán)境下船舶結(jié)構(gòu)設(shè)計(jì)及評(píng)估方法提供技術(shù)支撐。

    1 典型船冰相互作用場(chǎng)景

    船舶航行在極地水域時(shí),其艏部及涉冰帶結(jié)構(gòu)可能與浮冰、冰脊、冰山等發(fā)生接觸碰撞。對(duì)于船舶與浮冰發(fā)生的常規(guī)碰撞,主要發(fā)生在船舶回轉(zhuǎn)運(yùn)動(dòng)時(shí)舷側(cè)涉冰帶結(jié)構(gòu)、船舶直線航行時(shí)艏尖結(jié)構(gòu)等。對(duì)于船舶與冰脊或冰山發(fā)生碰撞情況,主要考慮到氣候條件的影響(如陽(yáng)光輻射、霧氣影響)、船員偶然操作失誤、探測(cè)儀器設(shè)備失效或誤判等因素。如雪龍?zhí)栐?019 年時(shí)曾因濃霧遮擋,能見(jiàn)度極地的環(huán)境下與冰山相撞,泰坦尼克號(hào)在天氣晴朗情況下與冰山相撞,導(dǎo)致舷側(cè)進(jìn)水,船體斷裂而沉沒(méi)[25]。對(duì)于極端環(huán)境條件,主要的船冰相互作用可以概括為:船舶艏部撞擊冰山(見(jiàn)圖1a)、船舶應(yīng)急避險(xiǎn)時(shí)舷側(cè)撞擊冰山(見(jiàn)圖1b)以及船舶快速行駛騎冰脊(見(jiàn)圖1c)等3 種場(chǎng)景。本文主要針對(duì)前兩種場(chǎng)景開(kāi)展艏部結(jié)構(gòu)在低溫極端環(huán)境下的損傷演化過(guò)程分析。

    圖1 極端船冰作用場(chǎng)景Fig.1 Extreme ship-ice action scenarios:a) Scenario 1;b) Scenario 2;c) Scenario 3

    2 仿真模型建立及材料本構(gòu)

    2.1 船體結(jié)構(gòu)建模及材料本構(gòu)模型

    文中采用等效船體梁方法[26]開(kāi)展局部立體結(jié)構(gòu)失效機(jī)理研究,對(duì)可能與海冰接觸并發(fā)生變形的艏部進(jìn)行了詳細(xì)的有限元建模,細(xì)化區(qū)域的甲板、外板、縱桁、橫梁及骨材等所有結(jié)構(gòu)均采用殼單元建立,定義為彈性體。對(duì)于遠(yuǎn)離碰撞發(fā)生區(qū)域的船體結(jié)構(gòu)進(jìn)行簡(jiǎn)化建模,定義為剛體材料,以提升計(jì)算效率,并通過(guò)改變不同站位材料密度方式調(diào)整重心位置、重量等與實(shí)船一致。其中,遠(yuǎn)離碰撞區(qū)域長(zhǎng)度為141 m,艏部細(xì)化區(qū)域長(zhǎng)度為46 m,寬度為16.5 m,型深為17 m,有限元模型見(jiàn)圖2 所示。模型中殼單元總數(shù)為49 631個(gè),節(jié)點(diǎn)數(shù)為44 011 個(gè)。數(shù)值仿真分析過(guò)程根據(jù)工況選擇不同的初始速度,無(wú)其他約束條件。針對(duì)船舶結(jié)構(gòu)與海冰碰撞場(chǎng)景,其船體的附件質(zhì)量取為0.02Ms(Ms為船體排水量)[27],海冰的附件質(zhì)量取為0.1Mi(Mi為海冰質(zhì)量)。

    圖2 船體有限元模型Fig.2 Finite element model of hull

    本文針對(duì)船用材料開(kāi)展了–40 ℃和常溫環(huán)境下的材料力學(xué)性能測(cè)試,獲得了船用高強(qiáng)度鋼應(yīng)力-應(yīng)變曲線,如圖3 所示。在保證材料曲線可以較好地模擬主要特征和變化趨勢(shì)的前提下,將材料曲線進(jìn)行合理簡(jiǎn)化,以減少材料動(dòng)態(tài)力學(xué)性能試驗(yàn)工作量。圖3中給出了真實(shí)材料應(yīng)力應(yīng)變曲線以及用于數(shù)值仿真分析中的簡(jiǎn)化應(yīng)力應(yīng)變曲線,彈性模量為206 GPa,泊松比為0.3,常溫下對(duì)應(yīng)的屈服應(yīng)力為472 MPa,–40 ℃環(huán)境下對(duì)應(yīng)的屈服應(yīng)力為500 MPa。在–60 ℃時(shí),該材料進(jìn)入韌脆轉(zhuǎn)變狀態(tài)。由于本文未開(kāi)展相應(yīng)的材料力學(xué)性能試驗(yàn),且文中的碰撞速度變化區(qū)間較?。?~15 kn),接近文獻(xiàn)[28]中的速度,屬于低應(yīng)變率范圍,因此本文采用Cowper-Symonds 模型[29]來(lái)描述應(yīng)變率對(duì)鋼材料影響:

    圖3 高強(qiáng)度鋼應(yīng)力-應(yīng)變曲線Fig.3 Stress-strain curve of high strength steel

    為了更加真實(shí)地模擬船冰碰撞,本文采用應(yīng)變準(zhǔn)則模擬鋼材失效,并根據(jù)文獻(xiàn)[8]中參數(shù)取鋼材最大塑性應(yīng)變值為0.28。

    2.2 海冰建模及材料本構(gòu)模型

    海冰模型采用Hex8 體單元建立,并考慮計(jì)算效率,在碰撞可能發(fā)生區(qū)域進(jìn)行網(wǎng)格細(xì)化,遠(yuǎn)離碰撞區(qū)域采用大尺度網(wǎng)格。2 種碰撞場(chǎng)景的冰體有限元模型如圖4 所示,海冰長(zhǎng)度方向?yàn)?0 m,寬度方向?yàn)?4 m,高度方向?yàn)?40 m。文中計(jì)算環(huán)境溫度包括常溫和–40 ℃,且為保證分析中考慮單一變量的影響,海冰密度統(tǒng)一取900 kg/m3。

    圖4 冰山有限元模型Fig.4 Finite element model of iceberg:a) finite element model of iceberg in ship bow impact;b) finite element model of icebergs in bow side impact

    在數(shù)值仿真中,將海冰材料性質(zhì)假定為各向同性。海冰失效準(zhǔn)則采用Von-Mises 準(zhǔn)則,海冰的破壞模式采用最大塑性應(yīng)變失效[30],海冰的主要力學(xué)性能參數(shù)見(jiàn)表1。

    表1 冰體材料參數(shù)Tab.1 Ice material parameters

    3 結(jié)構(gòu)損傷及參數(shù)敏感性分析

    2 種碰撞場(chǎng)景下,船冰相對(duì)高度方向的位置根據(jù)冰體受到的浮力以及船體結(jié)構(gòu)吃水來(lái)確定,以保證碰撞發(fā)生的位置更加接近真實(shí)情況。對(duì)于球鼻艏及艏柱位置首先與冰山發(fā)生碰撞,主要研究球鼻艏及艏柱位置結(jié)構(gòu)損傷,并開(kāi)展航速、材料本構(gòu)以及溫度影響分析。對(duì)于舷側(cè)與冰山發(fā)生碰撞情形,主要研究舷側(cè)板架結(jié)構(gòu)損傷情況,并開(kāi)展不同撞擊角度下的損傷分析。通過(guò)對(duì)上述2 種碰撞場(chǎng)景的分析,獲得不同碰撞場(chǎng)景狀態(tài)下主要失效部位的失效模式,以此掌握垮塌過(guò)程失效機(jī)理。

    3.1 船舶正撞冰山場(chǎng)景

    3.1.1 航速影響分析

    航速取6、8、10、12 kn 等4 種,撞擊場(chǎng)景見(jiàn)圖1a。通過(guò)圖5 可以看出,隨著航速的增大,碰撞力量值隨之增大,冰體變形能量不斷增加。4 s 后目標(biāo)船航速基本為0,對(duì)應(yīng)的撞深基本保持不變,且發(fā)現(xiàn)對(duì)應(yīng)的冰體變形能也基本保持不變。

    圖5 航速影響Fig.5 Effect of speed:a) collision force;b) energy

    3.1.2 材料本構(gòu)模型影響分析

    本節(jié)主要分析Cowper-Symonds 模型和真實(shí)材料模型對(duì)損傷過(guò)程及損傷范圍影響。本節(jié)仿真對(duì)應(yīng)的環(huán)境溫度為常溫,船體航速為10 kn。通過(guò)圖6 可知,2種材料模型下,碰撞力、能量的變化趨勢(shì)基本一致,僅在局部峰值處存在一定差異。因此,后續(xù)分析中均采用真實(shí)材料本構(gòu)模型。

    圖6 材料本構(gòu)模型影響Fig.6 Effect of material constitutive model:a) collision force;b) energy

    3.1.3 溫度影響分析

    本節(jié)主要針對(duì)常溫和–40 ℃低溫環(huán)境下材料性能變化對(duì)目標(biāo)船垮塌過(guò)程及損傷范圍等的影響進(jìn)行分析,仿真中船體航速為10 kn。通過(guò)圖7 可知,2種材料模型下碰撞力的變化趨勢(shì)基本一致,僅在局部峰值處存在一定差異。對(duì)于常溫環(huán)境下的材料,由于材料的屈服應(yīng)力為472 MPa,小于–40 ℃環(huán)境下的材料屈服應(yīng)力。因此,采用常溫環(huán)境下材料性能參數(shù)的目標(biāo)船艏部結(jié)構(gòu)的變形能要大于采用–40 ℃環(huán)境下材料性能參數(shù)的目標(biāo)船艏部結(jié)構(gòu)的變形能。

    圖7 環(huán)境溫度影響Fig.7 Effect of ambient temperature:a) collision force;b) energy

    3.1.4 船艏部撞擊冰山場(chǎng)景

    本節(jié)主要針對(duì)表2 中的4 個(gè)工況開(kāi)展仿真分析,研究不同冰體塑性應(yīng)變和航速組合情況下結(jié)構(gòu)損傷過(guò)程、失效范圍,計(jì)算工況環(huán)境溫度均為–40 ℃。

    表2 船艏部撞擊冰山場(chǎng)景計(jì)算工況Tab.2 Calculation condition for bow impact iceberg scenario

    本次計(jì)算過(guò)程中的艏部結(jié)構(gòu)損傷情況如圖8 所示。通過(guò)對(duì)比分析可以得出,隨著冰體塑性應(yīng)變的增大,冰體破碎范圍逐漸減小,結(jié)構(gòu)損傷范圍逐漸增大。主要結(jié)構(gòu)損傷發(fā)生在球鼻艏處,主要表現(xiàn)為球鼻艏板架結(jié)構(gòu)的垮塌,加強(qiáng)筋結(jié)構(gòu)的屈服失效。

    圖8 球艏損傷云圖Fig.8 Damage cloud of bow:a) condition 1;b) condition 2;c) condition 3;d) condition 4

    3.2 舷側(cè)部撞擊冰山場(chǎng)景仿真分析

    本節(jié)主要針對(duì)表3 中的6 個(gè)工況開(kāi)展仿真分析,研究不同撞擊角度、不同撞擊位置下的結(jié)構(gòu)損傷過(guò)程及失效模式等。其中強(qiáng)框位置、板格中心1 和板格中心2 位置如圖9 所示。

    表3 舷側(cè)部撞擊冰山場(chǎng)景仿真計(jì)算工況Tab.3 Calculation conditions for side impact iceberg scene

    圖9 3 個(gè)典型撞擊位置Fig.9 Three typical impact locations

    由于本場(chǎng)景為舷側(cè)與冰山發(fā)生碰撞,其航速的確定根據(jù)敞水狀態(tài),未及時(shí)發(fā)現(xiàn)冰山,保持15 kn 航速行駛確定。為分析其舷側(cè)板架結(jié)構(gòu)的失效機(jī)理,采取了不同航向角、不同撞擊位置進(jìn)行分析,同時(shí)分析不同冰體失效塑性應(yīng)變對(duì)結(jié)構(gòu)垮塌過(guò)程的影響,其中初始角度為船舶正向撞擊冰山為0°,向左舷轉(zhuǎn)彎為正值。

    不同工況的船冰相對(duì)位置初始狀態(tài)如圖10 所示,各工況在碰撞結(jié)束后的船冰相對(duì)位置、冰體破損情況、板架結(jié)構(gòu)損傷情況如圖11 所示??梢缘贸?,隨著初始角度的增大,碰撞結(jié)束后偏離的角度越大。主要是由于碰撞過(guò)程中冰山將給船體艏部一個(gè)側(cè)向力,使得船體發(fā)生整體的偏轉(zhuǎn),而隨著初始角度的增大,產(chǎn)生轉(zhuǎn)向的趨勢(shì)越明顯。隨著初始相對(duì)角度的增大,由于船體型線沿船長(zhǎng)的變化,使得船體與冰山接觸的面積增大,碰撞結(jié)束后,冰體破損范圍也越大。同時(shí),隨著冰體塑性失效應(yīng)變的增大,冰山破損范圍逐漸減小。碰撞過(guò)程呈現(xiàn)局部集中的特點(diǎn),產(chǎn)生塑性變形區(qū)域主要集中在船冰接觸位置。隨著冰體塑性失效應(yīng)變的增大,其損傷范圍逐漸增大,撞擊位置發(fā)生在板格2 附近時(shí),損傷范圍最大,板格1附近次之,最小損傷范圍為強(qiáng)框附近。發(fā)生該狀態(tài)的原因在于,強(qiáng)框位置剛度較大,而對(duì)板格而言,隨著跨長(zhǎng)的變大,其相對(duì)剛度降低,使得損傷變形范圍更大。但由于承受載荷的范圍更廣,應(yīng)力重新分布的范圍區(qū)域更大,使得產(chǎn)生的塑性應(yīng)變最大值在3 種狀態(tài)時(shí)基本一致。

    圖11 船冰相對(duì)位置及冰體破損情況Fig.11 Relative position of ship and ice damage:a) condition 1;b) condition 2;c) condition 3;d) condition 4;e) condition 5;f) condition 6

    6 種工況下的船冰碰撞力曲線如圖12 所示??梢钥闯?,碰撞過(guò)程僅持續(xù)了5 s 左右。對(duì)于撞擊強(qiáng)框、板格1 和板格2 不同位置,結(jié)構(gòu)較弱的首先產(chǎn)生較大的碰撞力,隨著結(jié)構(gòu)增加,產(chǎn)生極大碰撞力的時(shí)間向后推移。船舶與冰山發(fā)生碰撞時(shí)的碰撞力與船和浮冰、碎冰等產(chǎn)生的短暫脈沖形式的碰撞力明顯不同,呈現(xiàn)峰值數(shù)量更少、量值更大的特點(diǎn)。

    圖12 碰撞力曲線Fig.12 Collision force curve:a) conditions 1—3;b) conditions 4—6

    4 結(jié)論

    1)當(dāng)船體材料C-S 模型和真實(shí)材料模型均采用同樣的屈服應(yīng)力、彈性模量等基本參數(shù)時(shí),對(duì)碰撞力、船體結(jié)構(gòu)損傷以及結(jié)構(gòu)變形能量等的影響較小。

    2)隨著航速的增加,船體撞擊區(qū)域的損傷范圍增大,船體結(jié)構(gòu)變形能、撞擊深度以及碰撞力等增大。隨著環(huán)境溫度的降低,船用鋼材料屈服應(yīng)力增大,船體結(jié)構(gòu)損傷范圍減小,在材料達(dá)到韌脆轉(zhuǎn)變溫度前,低溫使得結(jié)構(gòu)更加安全。

    3)對(duì)于緊急避險(xiǎn)舷側(cè)碰撞情況,隨著撞擊角度的增加,船冰接觸面積增大,船體結(jié)構(gòu)損傷區(qū)域增加。隨著碰撞發(fā)生位置的剛度變大,碰撞區(qū)域損傷范圍變小。

    4)對(duì)于本文所涉及的2 種極端場(chǎng)景,艏部結(jié)構(gòu)失效模式主要表現(xiàn)為橫框架屈服失效、加強(qiáng)筋屈曲及側(cè)傾失效、局部外板大變形等。

    猜你喜歡
    冰體海冰冰山
    彈體高速侵徹冰體研究
    末次盛冰期以來(lái)巴倫支海-喀拉海古海洋環(huán)境及海冰研究進(jìn)展
    高速?gòu)楏w侵徹冰材料過(guò)程數(shù)值模擬研究
    冰體質(zhì)量和撞擊角度對(duì)船首結(jié)構(gòu)碰撞性能的影響
    基于船-水-冰耦合技術(shù)的撞擊參數(shù)對(duì)船冰碰撞性能的影響
    不可思議!冰山變綠了
    基于SIFT-SVM的北冰洋海冰識(shí)別研究
    崩塌的冰山
    危險(xiǎn)的冰山
    危險(xiǎn)的冰山
    啦啦啦啦在线视频资源| 一级,二级,三级黄色视频| 亚洲成av片中文字幕在线观看 | 99热这里只有是精品在线观看| 熟女人妻精品中文字幕| 欧美成人精品欧美一级黄| 欧美人与性动交α欧美软件 | 99久久人妻综合| 欧美人与善性xxx| 亚洲三级黄色毛片| 久久精品熟女亚洲av麻豆精品| 日韩成人伦理影院| 成年人免费黄色播放视频| 亚洲精品自拍成人| 日韩av在线免费看完整版不卡| 五月伊人婷婷丁香| 一级毛片 在线播放| 亚洲精品乱久久久久久| 极品人妻少妇av视频| 成人国语在线视频| 男女边吃奶边做爰视频| 狂野欧美激情性xxxx在线观看| 99九九在线精品视频| 少妇 在线观看| 免费黄网站久久成人精品| 夜夜爽夜夜爽视频| av有码第一页| 在线看a的网站| 国产黄频视频在线观看| 欧美精品国产亚洲| 久久青草综合色| 国产亚洲精品久久久com| a级毛色黄片| 亚洲av在线观看美女高潮| 成人综合一区亚洲| 大香蕉久久网| 国产福利在线免费观看视频| 18在线观看网站| 国产国语露脸激情在线看| 777米奇影视久久| 校园人妻丝袜中文字幕| 国产成人aa在线观看| 99热全是精品| 国产爽快片一区二区三区| 又黄又爽又刺激的免费视频.| 伊人亚洲综合成人网| 成人二区视频| 熟女人妻精品中文字幕| av免费观看日本| 黄片无遮挡物在线观看| 蜜臀久久99精品久久宅男| 天天操日日干夜夜撸| 国产又色又爽无遮挡免| 亚洲av电影在线观看一区二区三区| 丰满乱子伦码专区| 黑人猛操日本美女一级片| 精品第一国产精品| 国产成人免费观看mmmm| 80岁老熟妇乱子伦牲交| 婷婷色综合www| 亚洲欧美日韩卡通动漫| 国产视频首页在线观看| 激情五月婷婷亚洲| 亚洲av电影在线观看一区二区三区| 天堂8中文在线网| 国产av国产精品国产| 国产亚洲最大av| 校园人妻丝袜中文字幕| 久久青草综合色| 久久久久视频综合| a级毛片在线看网站| 国产熟女欧美一区二区| 免费在线观看完整版高清| 精品亚洲成a人片在线观看| 日韩一区二区三区影片| 五月天丁香电影| 日日爽夜夜爽网站| 亚洲,一卡二卡三卡| 99久久综合免费| 一级黄片播放器| 男女午夜视频在线观看 | av有码第一页| 精品国产露脸久久av麻豆| 一区二区三区精品91| 亚洲成国产人片在线观看| 成人影院久久| 9热在线视频观看99| 妹子高潮喷水视频| 日本91视频免费播放| 国产欧美亚洲国产| 国产精品 国内视频| 亚洲av.av天堂| 观看美女的网站| 国产精品嫩草影院av在线观看| 欧美成人精品欧美一级黄| 少妇的逼水好多| 国产精品久久久久久久电影| 又大又黄又爽视频免费| 国产成人av激情在线播放| 热re99久久精品国产66热6| 欧美97在线视频| 女性生殖器流出的白浆| 亚洲色图 男人天堂 中文字幕 | 久久久国产一区二区| 欧美bdsm另类| 一级a做视频免费观看| 久久 成人 亚洲| 熟女人妻精品中文字幕| 久久国产精品男人的天堂亚洲 | 国产精品秋霞免费鲁丝片| 观看av在线不卡| 国产精品一国产av| 在线天堂最新版资源| 欧美变态另类bdsm刘玥| 午夜福利乱码中文字幕| 国产午夜精品一二区理论片| 国产免费又黄又爽又色| 晚上一个人看的免费电影| av片东京热男人的天堂| 免费在线观看黄色视频的| 超色免费av| 久久免费观看电影| 国产免费福利视频在线观看| av有码第一页| 亚洲国产日韩一区二区| 亚洲精华国产精华液的使用体验| 人妻 亚洲 视频| 久久精品aⅴ一区二区三区四区 | 亚洲精品色激情综合| 午夜福利视频精品| 三级国产精品片| 精品一区二区三卡| 成年美女黄网站色视频大全免费| 精品午夜福利在线看| 黄色毛片三级朝国网站| 热re99久久精品国产66热6| 午夜福利乱码中文字幕| 亚洲一区二区三区欧美精品| av一本久久久久| 大香蕉97超碰在线| 亚洲国产精品国产精品| 免费看光身美女| 久久久亚洲精品成人影院| 精品久久久久久电影网| 国产免费一级a男人的天堂| 在线天堂中文资源库| 不卡视频在线观看欧美| 国产毛片在线视频| 久久精品aⅴ一区二区三区四区 | 色婷婷av一区二区三区视频| 99热国产这里只有精品6| 国产 一区精品| 哪个播放器可以免费观看大片| 免费观看a级毛片全部| 只有这里有精品99| 香蕉国产在线看| 99香蕉大伊视频| 久久午夜综合久久蜜桃| 在线观看www视频免费| 91成人精品电影| 成人免费观看视频高清| 日韩av不卡免费在线播放| 亚洲精品乱久久久久久| a级毛片黄视频| 97在线人人人人妻| 狠狠婷婷综合久久久久久88av| 看非洲黑人一级黄片| 亚洲av电影在线观看一区二区三区| a 毛片基地| 在线观看人妻少妇| 成人亚洲精品一区在线观看| 日本午夜av视频| av一本久久久久| 国产片特级美女逼逼视频| 久久久国产欧美日韩av| 亚洲精品aⅴ在线观看| 日韩三级伦理在线观看| 又黄又粗又硬又大视频| 中文乱码字字幕精品一区二区三区| 久久婷婷青草| 十八禁网站网址无遮挡| av免费观看日本| 日本av手机在线免费观看| 欧美xxxx性猛交bbbb| 极品少妇高潮喷水抽搐| 大话2 男鬼变身卡| 国产精品国产三级国产av玫瑰| 国产1区2区3区精品| 欧美日韩国产mv在线观看视频| 国产一级毛片在线| 亚洲精品美女久久av网站| 午夜福利网站1000一区二区三区| 国产精品秋霞免费鲁丝片| 午夜福利视频精品| 色5月婷婷丁香| 国产麻豆69| 天美传媒精品一区二区| 精品一品国产午夜福利视频| av福利片在线| 麻豆乱淫一区二区| 国产免费福利视频在线观看| 丝袜美足系列| 国产成人午夜福利电影在线观看| 成人国语在线视频| 18禁观看日本| av播播在线观看一区| 99九九在线精品视频| 免费看不卡的av| 亚洲av日韩在线播放| 久久久国产一区二区| 亚洲欧美日韩卡通动漫| 毛片一级片免费看久久久久| 五月天丁香电影| 久久人人爽人人片av| 激情视频va一区二区三区| 熟妇人妻不卡中文字幕| 亚洲欧美色中文字幕在线| 亚洲成av片中文字幕在线观看 | 欧美日韩精品成人综合77777| 中文天堂在线官网| 亚洲欧美精品自产自拍| 波多野结衣一区麻豆| 高清黄色对白视频在线免费看| 欧美xxⅹ黑人| 久久久久久久国产电影| 999精品在线视频| 国产欧美日韩一区二区三区在线| 秋霞伦理黄片| 国产日韩欧美在线精品| 97超碰精品成人国产| 在线天堂中文资源库| 寂寞人妻少妇视频99o| 亚洲一码二码三码区别大吗| 成人黄色视频免费在线看| 又粗又硬又长又爽又黄的视频| 精品少妇久久久久久888优播| 日本欧美国产在线视频| 999精品在线视频| 国产乱人偷精品视频| 国产不卡av网站在线观看| 最近中文字幕2019免费版| 亚洲一区二区三区欧美精品| 亚洲av成人精品一二三区| √禁漫天堂资源中文www| 大片免费播放器 马上看| 在线观看www视频免费| 日本欧美视频一区| 亚洲av电影在线进入| 亚洲一区二区三区欧美精品| 99久久中文字幕三级久久日本| 丰满迷人的少妇在线观看| 国产精品国产三级国产av玫瑰| 国产精品嫩草影院av在线观看| 午夜免费观看性视频| 欧美日韩一区二区视频在线观看视频在线| 波多野结衣一区麻豆| 中文精品一卡2卡3卡4更新| 欧美亚洲日本最大视频资源| 汤姆久久久久久久影院中文字幕| 午夜福利在线观看免费完整高清在| 18在线观看网站| 中国美白少妇内射xxxbb| 亚洲国产精品一区二区三区在线| 国产精品 国内视频| 波多野结衣一区麻豆| 十分钟在线观看高清视频www| 成人国语在线视频| 免费观看性生交大片5| 亚洲精品成人av观看孕妇| 午夜影院在线不卡| 亚洲av日韩在线播放| 国产国语露脸激情在线看| 国产男女内射视频| 夫妻午夜视频| 在线观看三级黄色| 在现免费观看毛片| 欧美亚洲日本最大视频资源| 午夜精品国产一区二区电影| 精品人妻在线不人妻| 国产av一区二区精品久久| 日本黄大片高清| 亚洲av日韩在线播放| 五月玫瑰六月丁香| 性色avwww在线观看| 精品99又大又爽又粗少妇毛片| 亚洲av在线观看美女高潮| 久久久久精品人妻al黑| 汤姆久久久久久久影院中文字幕| 国产成人精品婷婷| 在线观看美女被高潮喷水网站| 少妇 在线观看| 国产成人精品一,二区| 波野结衣二区三区在线| 街头女战士在线观看网站| 日韩三级伦理在线观看| 国产成人a∨麻豆精品| 亚洲欧美日韩卡通动漫| 国产精品一二三区在线看| 国产成人午夜福利电影在线观看| 亚洲经典国产精华液单| 日韩精品有码人妻一区| av片东京热男人的天堂| 亚洲国产精品专区欧美| 日本av手机在线免费观看| 女人被躁到高潮嗷嗷叫费观| 精品酒店卫生间| 少妇人妻 视频| 热re99久久国产66热| 国产深夜福利视频在线观看| 成人亚洲精品一区在线观看| 国产欧美日韩综合在线一区二区| 日日爽夜夜爽网站| 熟妇人妻不卡中文字幕| 日韩熟女老妇一区二区性免费视频| 国产精品一国产av| 久久久精品94久久精品| 国产免费一级a男人的天堂| 欧美国产精品一级二级三级| 久久久久网色| 精品一区二区三卡| 国产有黄有色有爽视频| 中文字幕人妻熟女乱码| 久久人人爽人人爽人人片va| 亚洲精品国产色婷婷电影| 日本欧美国产在线视频| 国产乱来视频区| 青青草视频在线视频观看| 97在线视频观看| 成人18禁高潮啪啪吃奶动态图| 国产黄色视频一区二区在线观看| 亚洲欧美色中文字幕在线| 欧美激情极品国产一区二区三区 | 在线观看www视频免费| 亚洲国产精品999| 国产极品天堂在线| 免费女性裸体啪啪无遮挡网站| 欧美bdsm另类| 少妇被粗大猛烈的视频| 日本欧美国产在线视频| 在线精品无人区一区二区三| 一级毛片 在线播放| av在线播放精品| 美女中出高潮动态图| 日韩精品有码人妻一区| 大香蕉97超碰在线| 国产伦理片在线播放av一区| 欧美 亚洲 国产 日韩一| 少妇 在线观看| 亚洲,欧美,日韩| 亚洲综合色惰| 人妻少妇偷人精品九色| 80岁老熟妇乱子伦牲交| 国产xxxxx性猛交| 国产福利在线免费观看视频| xxxhd国产人妻xxx| 免费少妇av软件| 美女主播在线视频| 五月天丁香电影| 日韩一区二区视频免费看| 亚洲图色成人| 久久精品熟女亚洲av麻豆精品| 亚洲天堂av无毛| 久久国内精品自在自线图片| 蜜桃国产av成人99| 国产一区二区三区综合在线观看 | 欧美日韩视频高清一区二区三区二| 熟妇人妻不卡中文字幕| 人成视频在线观看免费观看| 精品一区二区三区四区五区乱码 | 久久国内精品自在自线图片| 99re6热这里在线精品视频| 色哟哟·www| 亚洲图色成人| 成人黄色视频免费在线看| 最近的中文字幕免费完整| 黑丝袜美女国产一区| av天堂久久9| 国产日韩欧美亚洲二区| 日本与韩国留学比较| 我要看黄色一级片免费的| 亚洲在久久综合| 中国三级夫妇交换| 婷婷色av中文字幕| 精品人妻偷拍中文字幕| 99久久综合免费| 少妇被粗大的猛进出69影院 | 一本大道久久a久久精品| av有码第一页| 免费大片18禁| 国产成人午夜福利电影在线观看| 我的女老师完整版在线观看| 欧美成人精品欧美一级黄| 1024视频免费在线观看| 中文字幕精品免费在线观看视频 | 精品久久久精品久久久| 深夜精品福利| 国产精品三级大全| 乱码一卡2卡4卡精品| 成人手机av| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 国产精品一区二区在线不卡| 乱人伦中国视频| 国产成人精品久久久久久| 天堂8中文在线网| 精品视频人人做人人爽| 亚洲欧美成人综合另类久久久| av女优亚洲男人天堂| 国产成人精品福利久久| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网| 欧美最新免费一区二区三区| 91久久精品国产一区二区三区| 欧美成人精品欧美一级黄| 最近最新中文字幕大全免费视频 | 亚洲婷婷狠狠爱综合网| 亚洲国产精品999| 国产男人的电影天堂91| 国产日韩欧美亚洲二区| 国产极品天堂在线| 国产精品嫩草影院av在线观看| 免费av不卡在线播放| 看免费av毛片| 亚洲国产精品999| 国产成人精品婷婷| 老熟女久久久| 日产精品乱码卡一卡2卡三| 五月天丁香电影| 日韩,欧美,国产一区二区三区| 2018国产大陆天天弄谢| 精品亚洲乱码少妇综合久久| av播播在线观看一区| 国产 精品1| 久久久久久久久久久久大奶| 青青草视频在线视频观看| 亚洲精品久久成人aⅴ小说| av又黄又爽大尺度在线免费看| 日日啪夜夜爽| 国产欧美日韩综合在线一区二区| 少妇 在线观看| 人体艺术视频欧美日本| 一边亲一边摸免费视频| 免费黄网站久久成人精品| 午夜日本视频在线| 欧美国产精品va在线观看不卡| 国产午夜精品一二区理论片| 最近最新中文字幕免费大全7| 精品亚洲乱码少妇综合久久| 日产精品乱码卡一卡2卡三| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 婷婷色综合大香蕉| av女优亚洲男人天堂| 久久女婷五月综合色啪小说| 欧美性感艳星| 国产又爽黄色视频| 日日撸夜夜添| 极品少妇高潮喷水抽搐| 亚洲国产av新网站| 日韩免费高清中文字幕av| 日韩,欧美,国产一区二区三区| 狂野欧美激情性xxxx在线观看| 最近的中文字幕免费完整| 亚洲精品aⅴ在线观看| 美国免费a级毛片| 成年人免费黄色播放视频| 久久影院123| 久久精品aⅴ一区二区三区四区 | 在线天堂中文资源库| 久久97久久精品| 韩国高清视频一区二区三区| 国产乱人偷精品视频| 久久久久精品久久久久真实原创| 婷婷色av中文字幕| 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 国产在线免费精品| 久久青草综合色| kizo精华| 国产精品欧美亚洲77777| 老熟女久久久| 美女大奶头黄色视频| 国产在视频线精品| 精品国产一区二区三区久久久樱花| 免费av中文字幕在线| 国产精品久久久久久久久免| 国产黄色免费在线视频| 18禁在线无遮挡免费观看视频| 热99久久久久精品小说推荐| 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 欧美bdsm另类| 激情视频va一区二区三区| 母亲3免费完整高清在线观看 | 91国产中文字幕| 我要看黄色一级片免费的| 大香蕉久久网| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频| 少妇的逼好多水| 久久狼人影院| 亚洲综合色惰| 日日啪夜夜爽| 婷婷色麻豆天堂久久| 夜夜骑夜夜射夜夜干| 黄色一级大片看看| 久久99热这里只频精品6学生| 男女午夜视频在线观看 | 三级国产精品片| 亚洲精品自拍成人| 三级国产精品片| 22中文网久久字幕| 一二三四中文在线观看免费高清| 久热这里只有精品99| 成人国产av品久久久| 丝袜美足系列| 国产av一区二区精品久久| 大码成人一级视频| 成人综合一区亚洲| 最近最新中文字幕免费大全7| 午夜精品国产一区二区电影| 国产白丝娇喘喷水9色精品| 精品一区在线观看国产| 亚洲图色成人| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的| 精品熟女少妇av免费看| 亚洲精品久久午夜乱码| 热99久久久久精品小说推荐| 亚洲av.av天堂| 老熟女久久久| 亚洲国产精品专区欧美| 麻豆乱淫一区二区| 一级,二级,三级黄色视频| 亚洲国产色片| 亚洲精品自拍成人| 99热6这里只有精品| 精品人妻偷拍中文字幕| 欧美精品亚洲一区二区| 男女免费视频国产| 9热在线视频观看99| 国产精品 国内视频| 久久99热6这里只有精品| av播播在线观看一区| 国产一级毛片在线| 最新中文字幕久久久久| 三级国产精品片| 国产 一区精品| 纯流量卡能插随身wifi吗| 亚洲美女黄色视频免费看| 超碰97精品在线观看| 成人国语在线视频| 亚洲性久久影院| 大香蕉97超碰在线| 国产精品一国产av| 久久国产精品男人的天堂亚洲 | 国产白丝娇喘喷水9色精品| 97人妻天天添夜夜摸| 日本wwww免费看| 久久狼人影院| 热re99久久国产66热| 亚洲一区二区三区欧美精品| 精品亚洲乱码少妇综合久久| 两个人看的免费小视频| 国产免费福利视频在线观看| 亚洲国产欧美日韩在线播放| 有码 亚洲区| 日韩熟女老妇一区二区性免费视频| 水蜜桃什么品种好| 日韩不卡一区二区三区视频在线| 18+在线观看网站| 午夜老司机福利剧场| 免费人妻精品一区二区三区视频| 丝瓜视频免费看黄片| 久久久久人妻精品一区果冻| 精品国产露脸久久av麻豆| 青春草国产在线视频| 日韩大片免费观看网站| 高清视频免费观看一区二区| 午夜老司机福利剧场| 国产一区二区三区av在线| 亚洲三级黄色毛片| 黄色怎么调成土黄色| 成人毛片a级毛片在线播放| 亚洲精品国产av蜜桃| 国产精品欧美亚洲77777| 国产精品一区二区在线观看99| 97精品久久久久久久久久精品| 日本vs欧美在线观看视频| 亚洲,欧美精品.| 色婷婷久久久亚洲欧美| 免费在线观看黄色视频的| 久久久精品免费免费高清| av天堂久久9| 久久久久精品人妻al黑| 99九九在线精品视频| 欧美日韩一区二区视频在线观看视频在线| 欧美人与性动交α欧美软件 | 好男人视频免费观看在线| 老司机影院毛片| 免费日韩欧美在线观看| 母亲3免费完整高清在线观看 | 高清不卡的av网站| 亚洲人成77777在线视频| 国产精品熟女久久久久浪| 日韩一本色道免费dvd| 国产69精品久久久久777片| 婷婷色综合大香蕉| 高清在线视频一区二区三区| 99热网站在线观看| 男人添女人高潮全过程视频| 青春草亚洲视频在线观看| 久久精品国产亚洲av涩爱| 一区二区三区乱码不卡18| 久久99一区二区三区|