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

    熱金屬顆粒對(duì)地表森林可燃物的傳熱仿真分析

    2021-12-28 13:50:04李東暉胡海清郭妍
    森林工程 2021年6期
    關(guān)鍵詞:傳熱森林火災(zāi)有限元分析

    李東暉 胡海清 郭妍

    摘?要:為研究導(dǎo)線短路熔化形成的熱金屬顆粒對(duì)地表森林可燃物的傳熱機(jī)理,本文通過(guò)有限元仿真技術(shù),基于傳熱學(xué)理論,建立熱金屬顆粒與地表可燃物間的傳熱仿真模型,將墜落在地面的熱金屬顆粒作為熱源,考慮熱傳導(dǎo)、對(duì)流和輻射的作用,利用正交試驗(yàn)分析顆粒材質(zhì)、環(huán)境溫度、落地溫度與熔點(diǎn)的比值、顆粒球心嵌入深度倍數(shù)和顆粒半徑對(duì)可燃物加熱過(guò)程的影響。結(jié)果表明,顆粒材質(zhì)、落地溫度與熔點(diǎn)的比值、顆粒半徑是影響引燃的主要原因,且銅質(zhì)顆粒在墜落后對(duì)可燃物的引燃能力更強(qiáng)。此研究為避免熱金屬顆粒引起森林火災(zāi)提供參考。

    關(guān)鍵詞:熱金屬顆粒;傳熱;森林火災(zāi);有限元分析;引燃

    中圖分類號(hào):S762.1?文獻(xiàn)標(biāo)識(shí)碼:A?文章編號(hào):1006-8023(2021)06-0047-06

    Abstract:In order to study the heat transfer mechanism of hot metal particles formed by wire short-circuit melting to surface forest fuel, in this paper, the finite element analysis technology was used to establish the heat transfer simulation model between hot metal particle and fuel on ground based on the heat transfer theory. Taking hot metal particle as the heat source, considering the effects of heat conduction, convection and radiation, the effects of particle material, ambient temperature, ratio of landing temperature to melting point, multiple of particle center embedding depth and particle radius on the heating process of fuel were analyzed by orthogonal test. The results showed that material, ratio of landing temperature to melting point and radius were the main factors affecting ignition, and copper particles had stronger ignition ability to fuel. This study provides a reference for avoiding forest fire caused by hot metal particle.

    Keywords:Hot metal particles; heat transfer; wildfire; finite element analysis; ignite

    0?引言

    隨著自然環(huán)境的日益惡劣,木材已成為國(guó)際經(jīng)貿(mào)的戰(zhàn)略物資[1]。森林火災(zāi)是造成木材損失的主要原因,在災(zāi)害中釋放出的大量溫室氣體,將進(jìn)一步加劇極端氣候的發(fā)生頻率[2-3]。伴隨著人類社會(huì)基礎(chǔ)建設(shè)的不斷進(jìn)步,由輸電導(dǎo)線短路產(chǎn)生的熱金屬顆粒已成為一類全新的火源形式[4],并已引發(fā)多場(chǎng)森林火災(zāi),造成巨大的損失[5-6]。

    在該類事故中,導(dǎo)線金屬受短路電弧熱量加熱,熔化生成金屬液滴,受重力作用向地面墜落的過(guò)程中冷卻變成熱金屬顆粒。若墜落在地面上的熱金屬顆粒仍存有足夠熱量,將引燃地表森林可燃物,引發(fā)森林大火。因此,通過(guò)對(duì)熱金屬顆粒墜落后對(duì)地表可燃物的加熱過(guò)程進(jìn)行研究,可以更合理地預(yù)防該類森林的火災(zāi)發(fā)生。目前,學(xué)界多通過(guò)室內(nèi)模擬實(shí)驗(yàn)對(duì)各類熱金屬顆粒的引燃能力進(jìn)行研究。Hadden等[7]通過(guò)室內(nèi)模擬實(shí)驗(yàn),熱金屬顆粒的直徑和初始溫度共同決定森林可燃物的引燃概率;Urban 等[8-9]根據(jù)自制的試驗(yàn)裝置,確定在引燃形式相同時(shí),體積較大的熱金屬顆粒所需要的初始溫度更低;王蘇盼[10]通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)熱金屬顆粒在直接引燃可燃物時(shí),既是加熱源也充當(dāng)先導(dǎo)點(diǎn)燃源,而由無(wú)焰燃燒向有焰燃燒的轉(zhuǎn)變過(guò)程中僅充當(dāng)加熱源;此外,李夢(mèng)媛等[11]通過(guò)實(shí)驗(yàn)對(duì)熱金屬顆粒沉入可燃物內(nèi)部時(shí)的傳熱過(guò)程進(jìn)行研究?,F(xiàn)有的研究雖然能夠通過(guò)實(shí)驗(yàn)現(xiàn)象測(cè)定熱金屬顆粒的引燃能力,但無(wú)法從森林可燃物受熱分解的角度定量地描述各種熱金屬顆粒的區(qū)別。

    隨著計(jì)算機(jī)技術(shù)的不斷發(fā)展,有限元仿真技術(shù)[12-13]開(kāi)始被廣泛應(yīng)用于林火研究中,借用仿真技術(shù)可以直觀地獲得熱金屬顆粒對(duì)地表森林可燃物的傳熱情況。為此本文通過(guò)有限元數(shù)值分析方法,墜地后熱金屬顆粒對(duì)地表森林可燃物的熱量傳遞過(guò)程進(jìn)行分析,從而獲得熱金屬顆粒材質(zhì)、直徑、落地溫度、嵌入深度和周邊環(huán)境溫度對(duì)森林可燃物溫度變化的影響。

    1?研究對(duì)象

    由于森林內(nèi)的地表可燃物較為松軟,可將墜落后的熱金屬顆粒視作嵌入到平面內(nèi)部的球體[14]。熱金屬顆粒與地表森林可燃物的傳熱模型如圖1所示。圖1中:邊界①-⑦為模型計(jì)算圓柱形區(qū)域的軸心;邊界①-②-⑤-④-⑦-⑥為大氣環(huán)境;邊界④-⑤-⑨-⑧-⑥-⑦為地表森林可燃物;邊界③-⑤為地面;r為熱金屬顆粒的半徑;h為熱金屬顆粒球心嵌入地面的深度。

    2?研究方法

    在熱金屬顆粒與地表的森林可燃物接觸后,可燃物受熱釋放出可燃性氣體,當(dāng)氣體到達(dá)燃點(diǎn)后發(fā)生燃燒[15]。在分析熱金屬顆粒與地表森林可燃物間的換熱過(guò)程時(shí),首先需要根據(jù)傳熱模型建立仿真區(qū)域的實(shí)體模型,根據(jù)實(shí)際情況將材料屬性和接觸關(guān)系賦予實(shí)體模型。再賦予實(shí)體模型初始溫度和熱金屬顆粒的初始溫度,進(jìn)行初始熱場(chǎng)分析。然后根據(jù)初始熱場(chǎng)分析的結(jié)果,進(jìn)行瞬態(tài)熱場(chǎng)分析,獲得一定時(shí)間內(nèi)熱金屬顆粒與可燃物的溫度變化。如圖2所示。圖2中:T0為熱金屬顆粒的初始溫度;T1為可燃物的初始溫度(環(huán)境溫度);h為金屬顆粒球心嵌入地面的深度;r為熱金屬顆粒的半徑;t為熱場(chǎng)分析時(shí)間;t0為初始分析時(shí)刻;step(step)為設(shè)定仿真步數(shù);ti為時(shí)間步長(zhǎng);te為設(shè)定的傳熱計(jì)算時(shí)間。熱金屬顆粒對(duì)于地表森林可燃物的直接引燃常于12 s內(nèi)結(jié)束[10],本研究中te=12 s。在整個(gè)計(jì)算流程中,可以通過(guò)修改T0、T1、h和r對(duì)計(jì)算進(jìn)行調(diào)整,研究各條件下可燃物的溫度分布情況,實(shí)現(xiàn)穩(wěn)定、可靠和實(shí)用的計(jì)算平臺(tái)。

    3?仿真模型的建立

    3.1?傳熱計(jì)算模型

    根據(jù)傳熱學(xué)相關(guān)理論,熱金屬顆粒與地表森林可燃物間的瞬態(tài)傳熱分布計(jì)算方程為[12]:

    由于森林地表可燃物的主要成分以纖維素為主,本文將纖維素的傳熱系數(shù)作為地表森林可燃物的等效傳熱系數(shù)。地表森林可燃物的傳熱系數(shù)λe和溫度差(T0-T1)間的關(guān)系為[16]:

    熱金屬顆粒的輻射發(fā)射率受溫度變化影響,需要根據(jù)金屬材質(zhì)通過(guò)“發(fā)射率-溫度”函數(shù)表示[17-18]。由于森林可燃物的形貌及顏色較為多樣,為考慮最惡劣情況,將可燃物的表面發(fā)射率設(shè)置為1[19]。導(dǎo)線多以銅或鋁制成,本文以銅顆粒和鋁顆粒作為主要研究對(duì)象,金屬和可燃物的基本物理性能見(jiàn)表1[20-21]。

    3.2?有限元仿真模型

    根據(jù)計(jì)算模型,通過(guò)仿真軟件搭建有限元仿真模型[22],如圖3所示。該模型由空氣區(qū)域、金屬顆粒區(qū)域與森林可燃物區(qū)域組成,模型的外部邊界狀態(tài)與外部環(huán)境相同。為提高計(jì)算精度,該仿真模型為四面體二次網(wǎng)格,共含有200萬(wàn)個(gè)節(jié)點(diǎn),網(wǎng)格平均質(zhì)量大于0.9。金屬顆粒、空氣和森林可燃物依照表1進(jìn)行物理參數(shù)賦值,圖3中的T0、T1、h和r的定義與前面的定義相同。

    4?仿真分析

    4.1?正交實(shí)驗(yàn)仿真方案

    熱金屬顆粒在墜地時(shí)的最高溫度無(wú)法超過(guò)其自身熔點(diǎn)。在描述嵌入深度時(shí),也需要考慮球體的半徑尺寸,才能準(zhǔn)確描述球心的嵌入深度。在設(shè)計(jì)正交實(shí)驗(yàn)方案時(shí),需要以顆粒初始溫度與環(huán)境溫度的比值作為熱金屬顆粒初始溫度的特征參數(shù),以球心嵌入地表的深度與顆粒半徑的比值作為球心嵌入深度的特征參數(shù)。由此,通過(guò)將顆粒的半徑、球心嵌入地表的深度與顆粒半徑的比值、環(huán)境溫度、顆粒落地溫度與其熔點(diǎn)的比值作為變化參數(shù),通過(guò)仿真計(jì)算獲得可燃物的最終溫度。為提高計(jì)算研究效率,本文通過(guò)正交實(shí)驗(yàn)法,分別對(duì)銅、鋁2種材質(zhì)的熱金屬顆粒使用L25(56)正交表進(jìn)行仿真分析,共進(jìn)行50次計(jì)算,見(jiàn)表2。

    4.2?顆粒材質(zhì)對(duì)可燃物加熱過(guò)程的影響

    根據(jù)表2內(nèi)的各梯度布置,計(jì)算得到的部分地表森林可燃物溫度變化云圖如圖4和圖5所示。圖4和圖5中的T0、T1、r和h為對(duì)應(yīng)的各因素取值,2#、5#、21#和25#為對(duì)應(yīng)正交試驗(yàn)中的計(jì)算編號(hào)。

    通過(guò)對(duì)圖4和圖5內(nèi)的部分仿真結(jié)果進(jìn)行分析,可以獲得不同材質(zhì)熱金屬顆粒對(duì)地表森林可燃物的傳熱特性。由于銅的熔點(diǎn)大于鋁,則銅質(zhì)熱金屬顆粒的墜地溫度往往大于鋁質(zhì)顆粒。由于材料自身熱性能的不同,銅顆粒的溫度下降得較快,對(duì)可燃物的加熱過(guò)程較為迅速。而鋁顆粒的溫度降低較慢,對(duì)可燃物的加熱過(guò)程較為持續(xù)。此外,在材質(zhì)相同時(shí),半徑(r)較小的顆粒的降溫過(guò)程較為迅速,能夠在短時(shí)間內(nèi)將自身熱量傳向可燃物。該結(jié)果與Casey D Zak的實(shí)驗(yàn)研究相似[10-11],說(shuō)明仿真計(jì)算的結(jié)果與實(shí)際過(guò)程較為吻合。

    4.3?各因素對(duì)可燃物加熱過(guò)程的影響

    在完成正交試驗(yàn)計(jì)算后,各因素影響下可燃物的溫度均值隨梯度的變化如圖6所示。影響因素方差分析見(jiàn)表3。

    結(jié)合圖6和表3中的方差分析結(jié)果可以發(fā)現(xiàn),由于正交實(shí)驗(yàn)仿真結(jié)果的R2均等于0.99,說(shuō)明各因素變化可以解釋最終溫度99%的變化。當(dāng)P<0.05時(shí),即可認(rèn)為該參數(shù)的變化與最終溫度間存在顯著性差異關(guān)系。因此,在對(duì)相同材質(zhì)的熱金屬顆粒進(jìn)行分析時(shí),熱顆粒半徑(r)與熔點(diǎn)與落地溫度的比值是影響可燃物加熱過(guò)程的主要因素。

    綜合4.2中的計(jì)算結(jié)果,由于銅顆粒的落地溫度較高,且對(duì)可燃物的加熱過(guò)程較為迅速,更易直接造成可燃物發(fā)生有焰燃燒;而鋁顆粒的落地溫度相對(duì)較低,對(duì)可燃物的加熱過(guò)程較為持續(xù),可能會(huì)引起無(wú)焰燃燒轉(zhuǎn)化為有焰燃燒的森林火災(zāi)。當(dāng)材質(zhì)相同時(shí),半徑(r)較小的顆粒質(zhì)量較小,其相對(duì)表面積更大,只有當(dāng)其落地溫度較高時(shí)才會(huì)直接點(diǎn)燃地表可燃物。若落地溫度不能直接點(diǎn)燃地表可燃物,其散熱過(guò)程較快,無(wú)法以陰燃的形式緩慢點(diǎn)燃可燃物。而體積較大的顆粒即使落地溫度不能直接點(diǎn)燃地表可燃物,也能夠通過(guò)持續(xù)加熱的方式使可燃物陰燃,更易引起森林火災(zāi)。

    5?結(jié)論

    本文基于有限元仿真技術(shù)對(duì)熱金屬顆粒對(duì)地表森林可燃物的熱交換過(guò)程進(jìn)行建模,研究熱金屬顆粒半徑、顆粒嵌入深度倍數(shù)、環(huán)境溫度和熔點(diǎn)與落地溫度的比值對(duì)地表森林可燃物溫度變化的影響,得出以下結(jié)論。

    (1)熱金屬顆粒墜地時(shí)的溫度是影響可燃物溫升的主要原因。在其他條件相同時(shí),由于銅的熔點(diǎn)高于鋁,在墜地時(shí)銅質(zhì)金屬顆粒的溫度往往更高,對(duì)可燃物的引燃能力更強(qiáng)。

    (2)半徑較大的金屬顆粒能夠以較高溫度長(zhǎng)時(shí)間對(duì)地表可燃物進(jìn)行加熱,即使無(wú)法直接點(diǎn)燃,也可能通過(guò)緩慢加熱使其發(fā)生陰燃。

    (3)由于熱金屬顆粒的溫度高于其他常見(jiàn)的火源物,環(huán)境溫度并不會(huì)對(duì)引燃過(guò)程產(chǎn)生很大的影響。在對(duì)該類災(zāi)害防控時(shí),不應(yīng)以季節(jié)溫度劃分警戒等級(jí)。

    【參?考?文?獻(xiàn)】

    [1]盧宏亮,李佳恒,李國(guó)源,等.天保工程對(duì)中國(guó)木質(zhì)林產(chǎn)品貿(mào)易的影響[J].林業(yè)經(jīng)濟(jì)問(wèn)題,2020,40(6):634-642.

    LU H L, LI J H, LI G Y, et al. Research on the impact of natural forest protection project on China s forest product trade[J]. Issues of Forestry Economics, 2020, 40(6): 634-642.

    [2]楊波,谷會(huì)巖.大興安嶺3種樹(shù)種木材組分燃燒煙對(duì)珍珠梅種子萌發(fā)的影響[J].東北林業(yè)大學(xué)學(xué)報(bào),2020,48(12):27-37.

    YANG B, GU H Y. Effects of wood components combustion smoke of three tree species in Daxing an mountains on germination of pearl plum seed[J]. Journal of Northeast Forestry University, 2020, 48(12): 27-37.

    [3]王偉,張興龍,竇旭,等.中度火干擾對(duì)興安落葉松林土壤磷和鉀元素的影響[J].森林工程,2020,36(1):10-18.

    WANG W, ZHANG X L, DOU X, et al. Effect of moderate fire disturbance on soil phosphorus and potassium of Dahurian larch (Larxi gmelinii) forest[J]. Forest Engineering, 2020, 36(1): 10-18.

    [4]楊光,袁思博,舒立福,等.森林火災(zāi)中高能量火:飛火研究進(jìn)展[J].世界林業(yè)研究,2020,33(1):20-25.

    YANG G, YUAN S B, SHU L F, et al. Research progress of high-energy forest fire: spotting fire[J]. World Forestry Research, 2020, 33(1): 20-25.

    [5]BADGER S G. Large-loss fires in the United states in 2011[J]. NFPA Journal, 2012, 106(6): 60-67.

    [6]FERNANDEZ-PELLO A C. Wildland fire spot ignition by Sparks and firebrands[J]. Fire Safety Journal, 2017, 91: 2-10.

    [7]HADDEN R M, SCOTT S, LAUTENBERGER C, et al. Ignition of combustible fuel beds by hot particles: an experimental and theoretical study[J]. Fire Technology, 2011, 47(2): 341-355.

    [8]URBAN J L, ZAK C D, FERNANDEZ-PELLO C. Cellulose spot fire ignition by hot metal particles[J]. Proceedings of the Combustion Institute, 2015, 35(3): 2707-2714.

    [9]URBAN J L, ZAK C D, SONG J Y, et al. Smoldering spot ignition of natural fuels by a hot metal particle[J]. Proceedings of the Combustion Institute, 2017, 36(2): 3211-3218.

    [10]王蘇盼.飛火顆粒點(diǎn)燃的實(shí)驗(yàn)及機(jī)理研究[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2016.

    WANG S P. Experimental and theoretical study on spotting ignition[D]. Hefei: University of Science and Technology of China, 2016.

    [11]李夢(mèng)媛,陳海翔,張林鶴.金屬熱顆粒引燃可燃堆積材料的實(shí)驗(yàn)研究[J].火災(zāi)科學(xué),2017,26(3):140-146.

    LI M Y, CHEN H X, ZHANG L H. Experimental study on the ignition of combustible bulk materials by a hot metal particle[J]. Fire Safety Science, 2017, 26(3): 140-146.

    [12]胡海清,郭妍,李東暉.基于計(jì)算流體動(dòng)力學(xué)的燃燒箱熱學(xué)特性建模及優(yōu)化[J].森林工程,2020,36(1):60-67.

    HU H Q, GUO Y, LI D H. Based on computational fluid dynamics modeling and optimization of thermal characteristics of combustion tank[J]. Forest Engineering, 2020, 36(1): 60-67.

    [13]李靜靜,陳萌.鋰動(dòng)力電池電化學(xué)-熱特性建模及仿真研究[J].森林工程,2020,36(6):87-94.

    LI J J, CHEN M. Modeling and simulation study of electrochemical and thermal characteristics of lithium power battery[J]. Forest Engineering, 2020, 36(6): 87-94.

    [14]CASEY D Z. The Effect of particle properties on hot particle spot fire ignition[D]. Berkeley: University Of California Berkeley, 2015.

    [15]郭妍.黑龍江8種喬木可燃物熱解特征及釋放氣體研究[D].哈爾濱:東北林業(yè)大學(xué),2020.

    GUO Y. Study on pyrolysis characteristics and gas emissions of eight tree species in Heilongjiang Province[D]. Harbin: Northeast Forestry University, 2020.

    [16]楊曉曉.基于典型產(chǎn)物的纖維素快速熱解機(jī)理研究及模型構(gòu)建[D].北京:北京林業(yè)大學(xué),2020.

    YANG X X. Mechanism research and modelling of cellulose fast pyrolysis based on the typical products[D]. Beijing: Beijing Forestry University, 2020.

    [17]WEN C D, MUDAWAR I. Mathematical determination of emissivity and surface temperature of aluminum alloys using multispectral radiation thermometry[J]. International Communications in Heat and Mass Transfer, 2006, 33(9): 1063-1070.

    [18]YU K, ZHANG H Y, LIU Y, et al. Study of normal spectral emissivity of copper during thermal oxidation at different temperatures and heating times[J]. International Journal of Heat and Mass Transfer, 2019, 129: 1066-1074.

    [19]方祥,王海暉,陶駿駿,等.林火蔓延過(guò)程中輻射換熱和點(diǎn)燃特性分析[J].林業(yè)科學(xué),2019,55(3):97-105.

    FANG X, WANG H H, TAO J J, et al. Analysis on radiation heat transfer and the associated ignition during forest fire spread[J]. Scientia Silvae Sinicae, 2019, 55(3): 97-105.

    [20]蔣婷.木材水射流輔助激光加工有限元仿真及實(shí)驗(yàn)研究[D].哈爾濱:東北林業(yè)大學(xué),2020.

    JIANG T. Research on finite element simulation and experiment in waterjet assisted laser machining wood[D]. Harbin: Northeast Forestry University, 2020.

    [21]TSE S D, FERNANDEZ-PELLO A C. On the flight paths of metal particles and embers generated by power lines in high winds—a potential source of wildland fires[J]. Fire Safety Journal, 1998, 30(4): 333-356.

    [22]熊世磊,丁賽飛,王啟慧,等.籽瓜破碎取籽分離機(jī)機(jī)架的有限元分析及優(yōu)化[J].森林工程,2021,37(2):86-94.

    XIONG S L, DING S F, WANG Q H, et al. Finite element analysis and optimization of the frame of seed melon crushing separator[J]. Forest Engineering, 2021, 37(2): 86-94.

    猜你喜歡
    傳熱森林火災(zāi)有限元分析
    不同外墻保溫體系框架柱熱橋二維非穩(wěn)態(tài)傳熱比較
    景福鎮(zhèn)預(yù)防森林火災(zāi)方法措施分析
    綠色科技(2016年20期)2016-12-27 16:46:25
    淺談森林火災(zāi)的現(xiàn)代監(jiān)測(cè)技術(shù)
    自錨式懸索橋鋼箱梁頂推施工階段結(jié)構(gòu)分析
    隨機(jī)振動(dòng)載荷下發(fā)射裝置尾罩疲勞壽命分析
    航空兵器(2016年4期)2016-11-28 21:54:01
    有限元分析帶溝槽平封頭的應(yīng)力集中
    飛機(jī)起落架支撐桿強(qiáng)度有限元分析
    科技視界(2016年18期)2016-11-03 22:31:14
    高爐熱風(fēng)爐理論與設(shè)計(jì)
    氣候變化對(duì)森林火災(zāi)的影響及預(yù)防對(duì)策的探討
    科技視界(2016年6期)2016-07-12 19:39:42
    森林火災(zāi)的危害及重要滅火方法的分析
    男女午夜视频在线观看| 国产男女超爽视频在线观看| 欧美色视频一区免费| 亚洲中文日韩欧美视频| 99热网站在线观看| 最新的欧美精品一区二区| 欧美黑人精品巨大| 久久久久久免费高清国产稀缺| 久久青草综合色| 欧美黄色淫秽网站| 嫁个100分男人电影在线观看| 一本一本久久a久久精品综合妖精| 精品一区二区三区视频在线观看免费 | 欧美日韩视频精品一区| 欧美黄色片欧美黄色片| 纯流量卡能插随身wifi吗| 热re99久久国产66热| 中文字幕色久视频| 麻豆av在线久日| 大香蕉久久网| 国产又爽黄色视频| 国产精品1区2区在线观看. | 精品国产国语对白av| 久久草成人影院| 18在线观看网站| 国产不卡av网站在线观看| 久久婷婷成人综合色麻豆| 成人永久免费在线观看视频| 欧美精品av麻豆av| 在线观看日韩欧美| 日韩欧美免费精品| 天堂俺去俺来也www色官网| 我的亚洲天堂| 亚洲成国产人片在线观看| 18在线观看网站| 狠狠狠狠99中文字幕| 91精品三级在线观看| 色在线成人网| 亚洲欧美激情在线| 国产乱人伦免费视频| 怎么达到女性高潮| 精品久久蜜臀av无| 成年人黄色毛片网站| 欧美日韩亚洲综合一区二区三区_| 欧美激情 高清一区二区三区| 自线自在国产av| 欧美激情久久久久久爽电影 | 老熟妇仑乱视频hdxx| 激情在线观看视频在线高清 | 天天影视国产精品| 身体一侧抽搐| 色老头精品视频在线观看| 国产精品电影一区二区三区 | 日本黄色视频三级网站网址 | 天天影视国产精品| 国产精品影院久久| 777久久人妻少妇嫩草av网站| 欧美大码av| 久久这里只有精品19| 亚洲黑人精品在线| 亚洲男人天堂网一区| 久久久国产欧美日韩av| 国产无遮挡羞羞视频在线观看| 成人永久免费在线观看视频| 亚洲情色 制服丝袜| 一级毛片女人18水好多| 亚洲精品国产色婷婷电影| 久99久视频精品免费| www.自偷自拍.com| 久久热在线av| 夜夜夜夜夜久久久久| 十八禁人妻一区二区| 免费女性裸体啪啪无遮挡网站| 国产精品香港三级国产av潘金莲| 欧美最黄视频在线播放免费 | 无遮挡黄片免费观看| 一级毛片高清免费大全| 一进一出抽搐动态| 视频在线观看一区二区三区| 亚洲九九香蕉| 美女高潮喷水抽搐中文字幕| 精品一品国产午夜福利视频| 十八禁高潮呻吟视频| 丰满饥渴人妻一区二区三| 久久精品亚洲精品国产色婷小说| 久久人人爽av亚洲精品天堂| 国产成人av教育| 久久久久久久久久久久大奶| 叶爱在线成人免费视频播放| 午夜免费观看网址| 国产一区在线观看成人免费| 女人被躁到高潮嗷嗷叫费观| 亚洲一区二区三区不卡视频| 黄片播放在线免费| 欧洲精品卡2卡3卡4卡5卡区| 夜夜夜夜夜久久久久| av福利片在线| 国产黄色免费在线视频| 久久中文字幕人妻熟女| 久久久水蜜桃国产精品网| av欧美777| 免费看a级黄色片| 免费看a级黄色片| 午夜福利视频在线观看免费| 日韩免费av在线播放| 欧美+亚洲+日韩+国产| 亚洲精品中文字幕在线视频| 国产高清视频在线播放一区| 侵犯人妻中文字幕一二三四区| 九色亚洲精品在线播放| 亚洲国产精品合色在线| 亚洲va日本ⅴa欧美va伊人久久| 好看av亚洲va欧美ⅴa在| 少妇猛男粗大的猛烈进出视频| 淫妇啪啪啪对白视频| 国产在线精品亚洲第一网站| 黄色a级毛片大全视频| 一区二区三区国产精品乱码| 国产在线精品亚洲第一网站| 制服人妻中文乱码| 丰满迷人的少妇在线观看| 两个人看的免费小视频| 中文字幕av电影在线播放| 人妻丰满熟妇av一区二区三区 | 一级毛片高清免费大全| 看黄色毛片网站| 天天躁日日躁夜夜躁夜夜| 精品国产国语对白av| 老司机亚洲免费影院| 国产免费男女视频| 不卡一级毛片| 国产精品久久久av美女十八| 欧美黄色片欧美黄色片| 天天操日日干夜夜撸| 18禁黄网站禁片午夜丰满| 两人在一起打扑克的视频| 最新在线观看一区二区三区| 最新在线观看一区二区三区| 十分钟在线观看高清视频www| 久久久久视频综合| 亚洲精品乱久久久久久| 久久性视频一级片| 亚洲精品自拍成人| 国产精品久久视频播放| 久久精品成人免费网站| a级毛片在线看网站| 欧美精品人与动牲交sv欧美| 国内毛片毛片毛片毛片毛片| videosex国产| 日韩视频一区二区在线观看| 国产精品自产拍在线观看55亚洲 | 亚洲精品国产色婷婷电影| 午夜成年电影在线免费观看| 黄色片一级片一级黄色片| 免费观看a级毛片全部| 中文字幕人妻熟女乱码| 一进一出抽搐动态| 国产成人一区二区三区免费视频网站| 国产单亲对白刺激| 午夜成年电影在线免费观看| 国产精品一区二区免费欧美| 黄片播放在线免费| 高潮久久久久久久久久久不卡| 久久久精品免费免费高清| 欧美在线黄色| 欧美精品高潮呻吟av久久| 啦啦啦免费观看视频1| 亚洲精品一二三| 日韩欧美三级三区| 欧美日本中文国产一区发布| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人欧美| 老熟女久久久| 天天影视国产精品| 国产蜜桃级精品一区二区三区 | 欧美日韩亚洲国产一区二区在线观看 | 青草久久国产| 91国产中文字幕| 亚洲人成电影免费在线| 操出白浆在线播放| 国产亚洲欧美98| 自拍欧美九色日韩亚洲蝌蚪91| 欧美久久黑人一区二区| 怎么达到女性高潮| 亚洲一区二区三区不卡视频| 黄色成人免费大全| 激情视频va一区二区三区| 欧美 日韩 精品 国产| 色在线成人网| 亚洲精品一二三| 丝袜美腿诱惑在线| 丝袜美腿诱惑在线| 亚洲av日韩在线播放| 国产精品 国内视频| 免费看十八禁软件| 成人手机av| 制服诱惑二区| 久久这里只有精品19| 久久精品国产99精品国产亚洲性色 | 丝袜人妻中文字幕| 色精品久久人妻99蜜桃| 欧美人与性动交α欧美精品济南到| 成年人黄色毛片网站| 国产精品 国内视频| 久9热在线精品视频| 日本撒尿小便嘘嘘汇集6| 成熟少妇高潮喷水视频| 制服人妻中文乱码| 国精品久久久久久国模美| 真人做人爱边吃奶动态| a级毛片黄视频| 国产精品久久久久久精品古装| 精品久久久久久,| 一区二区日韩欧美中文字幕| 久久中文字幕人妻熟女| av网站免费在线观看视频| 一级作爱视频免费观看| 色婷婷av一区二区三区视频| 狂野欧美激情性xxxx| 国产成人欧美在线观看 | 999久久久国产精品视频| 婷婷精品国产亚洲av在线 | 亚洲免费av在线视频| xxxhd国产人妻xxx| 高清欧美精品videossex| 999久久久精品免费观看国产| 精品亚洲成a人片在线观看| 亚洲精品一二三| 俄罗斯特黄特色一大片| 黄片大片在线免费观看| 中文字幕色久视频| 99精品在免费线老司机午夜| 午夜福利影视在线免费观看| 身体一侧抽搐| 久久精品熟女亚洲av麻豆精品| 亚洲精品粉嫩美女一区| 男女之事视频高清在线观看| 成人黄色视频免费在线看| 中出人妻视频一区二区| 99热只有精品国产| 国产欧美亚洲国产| 又大又爽又粗| 首页视频小说图片口味搜索| 久久久久精品国产欧美久久久| 免费日韩欧美在线观看| 国产xxxxx性猛交| 国产亚洲精品久久久久5区| 亚洲成国产人片在线观看| 久久 成人 亚洲| 99国产极品粉嫩在线观看| 大香蕉久久网| 嫁个100分男人电影在线观看| 熟女少妇亚洲综合色aaa.| 在线十欧美十亚洲十日本专区| 欧美最黄视频在线播放免费 | 国产精品 国内视频| 亚洲免费av在线视频| 国产日韩一区二区三区精品不卡| 色婷婷久久久亚洲欧美| videos熟女内射| 午夜影院日韩av| 国产一区二区激情短视频| 久久精品亚洲av国产电影网| 19禁男女啪啪无遮挡网站| 啪啪无遮挡十八禁网站| bbb黄色大片| 18禁裸乳无遮挡动漫免费视频| 免费av中文字幕在线| 天天躁日日躁夜夜躁夜夜| 亚洲国产精品合色在线| 日韩欧美国产一区二区入口| 久久亚洲精品不卡| 大陆偷拍与自拍| 欧美日韩成人在线一区二区| 老司机福利观看| 99国产精品99久久久久| 精品少妇一区二区三区视频日本电影| 精品亚洲成国产av| 久久久久久久午夜电影 | 国产高清激情床上av| 国产精品自产拍在线观看55亚洲 | 岛国在线观看网站| 老熟妇仑乱视频hdxx| 他把我摸到了高潮在线观看| 免费观看人在逋| 日本wwww免费看| 夫妻午夜视频| 真人做人爱边吃奶动态| 激情在线观看视频在线高清 | 桃红色精品国产亚洲av| 高清av免费在线| 亚洲一区二区三区不卡视频| 中文字幕制服av| 宅男免费午夜| 精品福利永久在线观看| 麻豆av在线久日| 色尼玛亚洲综合影院| 又黄又爽又免费观看的视频| 国产深夜福利视频在线观看| 宅男免费午夜| 18在线观看网站| 91麻豆av在线| 在线国产一区二区在线| 少妇 在线观看| 国产欧美日韩一区二区三| 99国产精品一区二区蜜桃av | 在线十欧美十亚洲十日本专区| 老熟妇仑乱视频hdxx| 久久久精品免费免费高清| 中文字幕人妻丝袜制服| 欧美日韩福利视频一区二区| 欧美一级毛片孕妇| 国产成人免费无遮挡视频| 中文字幕精品免费在线观看视频| 国产精品国产高清国产av | 两个人免费观看高清视频| 日韩熟女老妇一区二区性免费视频| 成年人黄色毛片网站| 最近最新中文字幕大全电影3 | 亚洲综合色网址| 99国产精品99久久久久| 好男人电影高清在线观看| 啦啦啦 在线观看视频| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 97人妻天天添夜夜摸| 亚洲美女黄片视频| 午夜福利免费观看在线| 在线观看免费日韩欧美大片| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品国产综合久久久| 国产精品影院久久| 午夜福利免费观看在线| 69精品国产乱码久久久| 成人影院久久| 国产精品一区二区精品视频观看| 搡老岳熟女国产| 俄罗斯特黄特色一大片| 免费不卡黄色视频| 久久国产亚洲av麻豆专区| 老汉色∧v一级毛片| 国产蜜桃级精品一区二区三区 | 一级毛片精品| 国产无遮挡羞羞视频在线观看| 久久天堂一区二区三区四区| 欧美精品人与动牲交sv欧美| 成人三级做爰电影| 高清欧美精品videossex| 日韩成人在线观看一区二区三区| 色94色欧美一区二区| 18禁黄网站禁片午夜丰满| av国产精品久久久久影院| 丝袜在线中文字幕| 老汉色av国产亚洲站长工具| 久久性视频一级片| 丰满迷人的少妇在线观看| 黄网站色视频无遮挡免费观看| 国产一卡二卡三卡精品| 久久精品91无色码中文字幕| 国产精品欧美亚洲77777| 99国产精品一区二区三区| 国产男靠女视频免费网站| 精品一区二区三卡| 大型av网站在线播放| 老司机福利观看| 久久国产精品影院| 国产精品免费一区二区三区在线 | 国产亚洲欧美精品永久| 欧美乱色亚洲激情| 国产精品久久久av美女十八| 国产成人精品无人区| 亚洲精品国产一区二区精华液| 亚洲欧美精品综合一区二区三区| 精品视频人人做人人爽| 久久99一区二区三区| xxx96com| 国产乱人伦免费视频| 国产不卡av网站在线观看| 国产免费现黄频在线看| 1024视频免费在线观看| 亚洲avbb在线观看| 一区二区三区国产精品乱码| 丝袜在线中文字幕| 久久香蕉国产精品| 久久久国产成人免费| 黄色a级毛片大全视频| 在线观看66精品国产| 波多野结衣av一区二区av| 精品国产乱码久久久久久男人| 一级片'在线观看视频| 变态另类成人亚洲欧美熟女 | 国产欧美日韩综合在线一区二区| 怎么达到女性高潮| 欧美 日韩 精品 国产| 国产一卡二卡三卡精品| av超薄肉色丝袜交足视频| 久久午夜综合久久蜜桃| 国产精品综合久久久久久久免费 | 中亚洲国语对白在线视频| 国产精品免费大片| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 亚洲五月婷婷丁香| 黄色a级毛片大全视频| 十八禁高潮呻吟视频| 欧美日韩亚洲综合一区二区三区_| 亚洲第一欧美日韩一区二区三区| 欧美性长视频在线观看| 国产精品香港三级国产av潘金莲| 超色免费av| 香蕉久久夜色| 一区二区日韩欧美中文字幕| cao死你这个sao货| 欧美国产精品一级二级三级| 水蜜桃什么品种好| 午夜福利一区二区在线看| av超薄肉色丝袜交足视频| www.熟女人妻精品国产| 超碰97精品在线观看| 在线永久观看黄色视频| 狠狠狠狠99中文字幕| av一本久久久久| 老汉色∧v一级毛片| 91成年电影在线观看| 在线av久久热| 曰老女人黄片| 亚洲精品一二三| 精品无人区乱码1区二区| 国产精品99久久99久久久不卡| 欧美黑人欧美精品刺激| 亚洲一区高清亚洲精品| 9热在线视频观看99| 99re在线观看精品视频| 国产黄色免费在线视频| 夫妻午夜视频| 老熟妇仑乱视频hdxx| 亚洲九九香蕉| 法律面前人人平等表现在哪些方面| 热re99久久精品国产66热6| 久久久水蜜桃国产精品网| 久久久久久亚洲精品国产蜜桃av| 亚洲一区中文字幕在线| 国产蜜桃级精品一区二区三区 | 91精品国产国语对白视频| 91成人精品电影| 免费看十八禁软件| 亚洲精品久久午夜乱码| 亚洲少妇的诱惑av| av超薄肉色丝袜交足视频| 一级毛片精品| 欧美在线黄色| 18在线观看网站| 免费一级毛片在线播放高清视频 | 国产成人一区二区三区免费视频网站| 在线十欧美十亚洲十日本专区| 国产一卡二卡三卡精品| 美女 人体艺术 gogo| 欧美日韩av久久| 国产欧美日韩精品亚洲av| 中文字幕色久视频| 国产深夜福利视频在线观看| 99精品欧美一区二区三区四区| 老司机午夜十八禁免费视频| 一级a爱视频在线免费观看| 老汉色av国产亚洲站长工具| 国产精品亚洲av一区麻豆| 女人被躁到高潮嗷嗷叫费观| 欧美在线一区亚洲| 在线观看舔阴道视频| 91麻豆精品激情在线观看国产 | 久久精品成人免费网站| 天天躁日日躁夜夜躁夜夜| 亚洲,欧美精品.| 成熟少妇高潮喷水视频| 新久久久久国产一级毛片| 欧美日韩视频精品一区| 激情在线观看视频在线高清 | 欧美大码av| 亚洲国产精品sss在线观看 | 人人妻人人爽人人添夜夜欢视频| 两人在一起打扑克的视频| 欧美成人午夜精品| 久久人妻福利社区极品人妻图片| 久久午夜综合久久蜜桃| 精品视频人人做人人爽| 婷婷丁香在线五月| 久久国产乱子伦精品免费另类| 老司机午夜十八禁免费视频| 国产99白浆流出| 少妇裸体淫交视频免费看高清 | 999精品在线视频| 最近最新中文字幕大全免费视频| 亚洲欧美日韩高清在线视频| 免费黄频网站在线观看国产| 欧美色视频一区免费| 激情在线观看视频在线高清 | 亚洲 欧美一区二区三区| 日韩欧美一区视频在线观看| 欧美黄色淫秽网站| 欧美 日韩 精品 国产| 男人操女人黄网站| 欧美日韩精品网址| 日本黄色日本黄色录像| 精品欧美一区二区三区在线| 大型黄色视频在线免费观看| 青草久久国产| 一边摸一边做爽爽视频免费| 男女免费视频国产| 91九色精品人成在线观看| 国产一区二区三区综合在线观看| 国产精品 欧美亚洲| videos熟女内射| 欧美色视频一区免费| 老熟妇乱子伦视频在线观看| 露出奶头的视频| 9色porny在线观看| 成年动漫av网址| 最新在线观看一区二区三区| 欧美 日韩 精品 国产| 视频在线观看一区二区三区| 国产精品乱码一区二三区的特点 | 91精品国产国语对白视频| www.999成人在线观看| 国内久久婷婷六月综合欲色啪| 涩涩av久久男人的天堂| 欧美乱妇无乱码| 少妇 在线观看| 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 99热只有精品国产| 亚洲精品国产一区二区精华液| 亚洲精品在线观看二区| 大香蕉久久成人网| 国产精品一区二区精品视频观看| 欧美日韩乱码在线| 一区二区三区国产精品乱码| 亚洲熟妇中文字幕五十中出 | 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 亚洲精品自拍成人| 国产精品乱码一区二三区的特点 | 国产99久久九九免费精品| 免费在线观看完整版高清| www.精华液| 国产xxxxx性猛交| 女人精品久久久久毛片| 嫩草影视91久久| 免费观看a级毛片全部| 一级黄色大片毛片| 久久久久久久精品吃奶| 久久精品国产综合久久久| 久久香蕉国产精品| 人人妻人人添人人爽欧美一区卜| 国产高清激情床上av| 久久草成人影院| 免费黄频网站在线观看国产| 日本wwww免费看| 国产欧美日韩一区二区三区在线| 两性夫妻黄色片| 麻豆乱淫一区二区| 午夜福利在线观看吧| 欧美日韩av久久| 国产成人精品久久二区二区91| 亚洲一区二区三区不卡视频| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 老汉色av国产亚洲站长工具| 丰满迷人的少妇在线观看| av视频免费观看在线观看| 欧美成人免费av一区二区三区 | 国产精品二区激情视频| 免费久久久久久久精品成人欧美视频| 搡老岳熟女国产| 久久久久国产一级毛片高清牌| 丝瓜视频免费看黄片| 精品福利永久在线观看| 亚洲av成人一区二区三| 国产精品九九99| 日日爽夜夜爽网站| 岛国毛片在线播放| 亚洲精品av麻豆狂野| 国产精品九九99| 中文字幕av电影在线播放| 一本一本久久a久久精品综合妖精| av天堂久久9| 宅男免费午夜| xxx96com| 亚洲av成人av| 国产1区2区3区精品| 国产免费现黄频在线看| 亚洲成人手机| 成在线人永久免费视频| 亚洲欧美一区二区三区久久| 日韩视频一区二区在线观看| 在线观看舔阴道视频| 精品国产一区二区三区久久久樱花| 99久久99久久久精品蜜桃| 极品教师在线免费播放| 美国免费a级毛片| 狠狠婷婷综合久久久久久88av| 国产免费男女视频| 国内毛片毛片毛片毛片毛片| 国产免费男女视频| 高清黄色对白视频在线免费看| 757午夜福利合集在线观看| 亚洲av欧美aⅴ国产| 亚洲精品一卡2卡三卡4卡5卡| 中国美女看黄片| 中文亚洲av片在线观看爽 | 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | 男女下面插进去视频免费观看| 又黄又爽又免费观看的视频| 老鸭窝网址在线观看| 在线观看免费视频日本深夜|