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

    碟式太陽(yáng)能鏡面組件面形誤差影響分析

    2017-07-01 16:28:31黃振軍梁偉青周順軍唐志雄
    電力與能源 2017年3期
    關(guān)鍵詞:檢測(cè)

    張 琛,黃振軍,梁偉青,周順軍,唐志雄

    (上海齊耀動(dòng)力技術(shù)有限公司,上海 201203)

    碟式太陽(yáng)能鏡面組件面形誤差影響分析

    張 琛,黃振軍,梁偉青,周順軍,唐志雄

    (上海齊耀動(dòng)力技術(shù)有限公司,上海 201203)

    運(yùn)用Matlab高效的數(shù)值計(jì)算能力進(jìn)行蒙特卡洛光線追跡模擬太陽(yáng)光的反射聚焦,實(shí)現(xiàn)編程計(jì)算結(jié)果的可視化,進(jìn)而分析碟式拋物面聚光器焦斑處的能流密度分布。通過(guò)分析入射太陽(yáng)張角、反射表面斜率誤差及其散射誤差對(duì)單塊鏡面組件焦斑處能量分布的影響,并通過(guò)Soltrace軟件論證了算法的準(zhǔn)確性。最終根據(jù)光斑實(shí)際形狀以及熱頭位置靶點(diǎn)尺寸,確立合理的面形精度誤差指標(biāo),為鏡面組件制造商的生產(chǎn)加工及面形誤差檢測(cè)提供理論依據(jù)及檢測(cè)標(biāo)準(zhǔn)。

    碟式拋物面聚光器;斜率誤差;蒙特卡洛法;Matlab;Soltrace

    2016年9月1日,國(guó)家發(fā)展改革委下發(fā)《關(guān)于太陽(yáng)能熱發(fā)電標(biāo)桿上網(wǎng)電價(jià)政策的通知》,核定太陽(yáng)能熱發(fā)電標(biāo)桿上網(wǎng)電價(jià)為1.15元/kWh,同時(shí)鼓勵(lì)地方政府對(duì)太陽(yáng)能熱發(fā)電企業(yè)采取稅費(fèi)減免、財(cái)政補(bǔ)貼、綠色信貸、土地優(yōu)惠等措施,促進(jìn)太陽(yáng)能熱發(fā)電產(chǎn)業(yè)發(fā)展。目前主流的太陽(yáng)能熱發(fā)電技術(shù)路線主要包括槽式、塔式、碟式以及線性菲涅爾4種。碟式斯特林熱發(fā)電系統(tǒng)光電轉(zhuǎn)換效率最高,適用于分布式發(fā)電。拋物面型鏡面組件的反射聚焦集熱作為碟式太陽(yáng)能光熱發(fā)電的關(guān)鍵技術(shù)之一,光學(xué)性能尤其是面形精度對(duì)發(fā)電系統(tǒng)的整體聚光性能影響很大,因此研究分析鏡面組件面形誤差對(duì)光斑分布的影響具有重要意義。

    1 數(shù)學(xué)建模

    根據(jù)碟式太陽(yáng)能的旋轉(zhuǎn)拋物面反射聚焦特性及蒙特卡羅光線追跡法,生成隨機(jī)的入射點(diǎn)坐標(biāo),組建單塊鏡面組件的數(shù)學(xué)模型。將接收器平面網(wǎng)格化,統(tǒng)計(jì)各網(wǎng)格的能量分布從而得出整個(gè)接收面的能流密度分布,最后通過(guò)邊界尺寸分析給出不同規(guī)格鏡面組件的斜率誤差精度要求。

    1.1 鏡面組件參數(shù)

    旋轉(zhuǎn)拋物面理論方程為x2+y2=4fz,其中f為拋物面焦距。考慮加工制造難度以及鏡面組件獨(dú)立調(diào)光聚焦的可行性,在實(shí)際工程應(yīng)用中通常簡(jiǎn)化為多圈等角度均勻分布的曲面反射鏡組件。建立數(shù)學(xué)模型時(shí)首先定義單塊鏡面組件內(nèi)徑為r,外徑為R,扇形張角為α,再利用Matlab軟件自帶的unifrnd函數(shù)生成兩組指定范圍內(nèi)均勻分布的隨機(jī)數(shù),隨機(jī)數(shù)的取值范圍分別定義為rand1∈(1,(R/r)2),rand2∈(1,απ/180),計(jì)算出各入射點(diǎn)空間坐標(biāo)(x0,y0,z0)。為在滿足仿真計(jì)算精度的同時(shí)盡量減少程序計(jì)算所需時(shí)間和內(nèi)存消耗,追跡入射光線的總數(shù)取106~107為宜。大量隨機(jī)分布的入射點(diǎn)定義出鏡面組件的三維空間拋物面形如圖1所示。

    (1)

    圖1 鏡面組件定義

    1.2 太陽(yáng)形狀分布

    (2)

    圖2 歸一化太陽(yáng)光強(qiáng)分布

    1.3 面形誤差計(jì)算

    鏡面組件生產(chǎn)加工過(guò)程中,聚光器實(shí)際表面輪廓與理論表面輪廓不一致引起的誤差,可分為位置誤差和傾斜誤差[2]。其中入射點(diǎn)位置與理論值不一致定義為位置誤差,入射點(diǎn)表面的斜率與理論值不一致定義為傾斜誤差,又稱斜率誤差。位置誤差對(duì)鏡面組件聚光性能的影響遠(yuǎn)遠(yuǎn)小于斜率誤差造成的影響[3]。因此,以表面斜率誤差來(lái)表征鏡面組件面形誤差,對(duì)于理論分析面形誤差對(duì)反射聚光性能的影響更為直接準(zhǔn)確,表面斜率誤差近似服從高斯分布,標(biāo)準(zhǔn)差用σslope表示。另外入射光線經(jīng)實(shí)際鏡面反射后并非理想鏡面反射,實(shí)際在各出射方向均有光強(qiáng)分布,在光束截面上以高斯形式分布,散射誤差標(biāo)準(zhǔn)差用σspec表示,各類型誤差分布如圖3所示。

    圖3 光學(xué)誤差分布

    1.4 空間坐標(biāo)系轉(zhuǎn)換

    在初始空間直角坐標(biāo)系OXiYiZi中,將生成的正態(tài)分布夾角數(shù)據(jù),以旋轉(zhuǎn)拋物面上各入射點(diǎn)位置的切平面法線方向重新建立空間直角坐標(biāo)系OXjYjZj??臻g坐標(biāo)旋轉(zhuǎn)轉(zhuǎn)換如圖4所示,初始空間直角坐標(biāo)系OXiYiZi先繞Z軸旋轉(zhuǎn)角度θ可得空間直角坐標(biāo)系OXmYmZm,再繞Y旋轉(zhuǎn)角度α,可得到目標(biāo)空間直角坐標(biāo)系OXjYjZj。代入θ及α角度值計(jì)算方向余弦矩陣[iRj(θ,α)],即可實(shí)現(xiàn)空間直角坐標(biāo)系轉(zhuǎn)換計(jì)算出絕對(duì)反射光線方向矢量。

    (3)

    圖4 空間坐標(biāo)旋轉(zhuǎn)轉(zhuǎn)換

    1.5 能流密度分布

    文獻(xiàn)[4,5]在計(jì)算旋轉(zhuǎn)拋物面型聚光器接收面的能流密度分布時(shí),均采用將接收面按照徑向等分的同心圓來(lái)進(jìn)行網(wǎng)格劃分的方法,通過(guò)統(tǒng)計(jì)各圓環(huán)內(nèi)的反射光線分布數(shù)量以及不同角度反射光線能量大小得出接收面的能流密度分布。碟式太陽(yáng)能鏡面組件在實(shí)際工程應(yīng)用中安裝角度位置需要專業(yè)調(diào)整,并非以理論焦平面或焦點(diǎn)作為投射靶面、靶點(diǎn)。由此可知,實(shí)際反射聚焦后的能量光斑中心并非和焦點(diǎn)重合,同心圓環(huán)形網(wǎng)格劃分不能準(zhǔn)確表述光斑形狀分布,并且徑向環(huán)寬相等造成各圓環(huán)的面積均不等,得出的能流密度分布與實(shí)際能流密度誤差較大。采用正方形網(wǎng)格劃分,將統(tǒng)計(jì)得到的各網(wǎng)格內(nèi)光線能量總和加載到矩陣內(nèi),再調(diào)用Imagesc函數(shù)進(jìn)行Flux Map能流密度分布圖的繪制。

    2 算法驗(yàn)證

    為了驗(yàn)證Matlab算法的準(zhǔn)確性,選取單塊外圈鏡面進(jìn)行光學(xué)仿真分析,再與Soltrace軟件進(jìn)行仿真比對(duì)。Soltrace光學(xué)仿真軟件由美國(guó)國(guó)家可再生能源實(shí)驗(yàn)室(NREL)研究開發(fā),可用于對(duì)太陽(yáng)能光學(xué)系統(tǒng)建模、仿真及性能分析[6]。其采用Delphi高級(jí)語(yǔ)言編寫,基于蒙特卡羅光線追跡法,考慮太陽(yáng)能量分布及元素的光學(xué)誤差分布,計(jì)算效率及精度較高。用戶可根據(jù)不同需求建立塔式、碟式及槽式等各類聚光系統(tǒng)并進(jìn)行聚光光斑能流密度分析。

    各圈鏡面組件尺寸參數(shù)如表1所示。考慮以焦平面位置作為目標(biāo)靶面,取單塊外圈鏡面組件作為反射元件,Soltrace軟件中Slope error設(shè)為0.95 m·r,Specularity error設(shè)為0.2 mard。能流密度分布結(jié)果比對(duì)如圖5所示。Soltrace與Matlab光學(xué)仿真結(jié)果對(duì)比表明,單塊外圈鏡面組件能流密度峰值約為290 kW/m2,光斑短軸處于±80 mm坐標(biāo)范圍內(nèi),光斑長(zhǎng)軸處于±120 mm坐標(biāo)范圍內(nèi),兩者基本吻合。通過(guò)長(zhǎng)短軸截面能流密度分布比對(duì),證明該算法分析鏡面組件面形誤差影響是準(zhǔn)確可行的。

    圖5 能流密度分布結(jié)果比對(duì)

    表1 鏡面組件尺寸參數(shù)

    3 面形誤差分析

    在距焦平面后特定距離的目標(biāo)靶平面上,僅考慮太陽(yáng)張角忽略其他光學(xué)誤差,根據(jù)旋轉(zhuǎn)拋物面鏡面組件規(guī)格,接收器錐角、深度以及加熱器內(nèi)外徑等關(guān)鍵尺寸確定好各鏡面組件的靶點(diǎn)位置,具體分布如圖6所示。下面取右側(cè)水平方向A、B、C三塊鏡面組件作為分析對(duì)象。

    圖6 靶點(diǎn)位置分布

    中空環(huán)形光斑分布如圖7所示。完整光斑分布確定單塊鏡面組件光斑的內(nèi)外邊緣界限,在圖7中以實(shí)線圓環(huán)表示。以A、B、C各靶點(diǎn)為圓心的與實(shí)線圓環(huán)最大內(nèi)切圓確定各圈鏡面組件的最大允許光斑直徑分別為191.0 、270.8 、293.4 mm。

    圖7 中空環(huán)形光斑分布

    考慮目標(biāo)靶上不同安裝位置鏡面組件的反射光斑分布角度各異,并且光斑呈橢圓形高斯分布,光斑長(zhǎng)軸基本與目標(biāo)靶徑向方向一致。取各光斑長(zhǎng)軸尺寸作為特征尺寸,定義為單塊鏡面組件反射的光斑等效直徑。通過(guò)Matlab光學(xué)仿真計(jì)算出不同面形誤差設(shè)置下光斑的等效直徑,繪制出等效光斑直徑的變化曲線,再擬合出面形誤差與鏡面組件聚焦光斑等效直徑的線性關(guān)系,如圖8所示。為保證各圈鏡面組件光斑都處于靶點(diǎn)位置分布中兩實(shí)線圓環(huán)內(nèi),根據(jù)各圈最大允許光斑直徑,可反推得出內(nèi)、中、外圈鏡面允許的最大面形誤差值分別為0.64 、1.09 、0.88 mrad。

    圖8 不同面形誤差下光斑有效直徑

    4 面形誤差檢測(cè)

    目前成功應(yīng)用于太陽(yáng)能光熱發(fā)電領(lǐng)域反射鏡面面形誤差檢測(cè)的技術(shù)方法,可大致分為兩類:一類為基于激光束掃描原理的VSHOT測(cè)量法;二類為基于圖像處理的檢測(cè)技術(shù),包括攝影法和偏折法(或稱為條紋反射法)等[7]。其中,應(yīng)用較廣泛的3D掃描技術(shù),通過(guò)在曲面鏡表面布置一定間距的標(biāo)記點(diǎn)獲取若干離散點(diǎn)云坐標(biāo),擬合生成面形數(shù)據(jù),但只能反映玻璃表面的面形誤差,無(wú)法測(cè)得折射率、厚度不均及鍍銀反射層表面的面形誤差等因素對(duì)鏡面組件聚光性能的影響[8]。

    文獻(xiàn)[9]以長(zhǎng)程輪廓儀為設(shè)計(jì)基礎(chǔ)研制出反射鏡面形檢測(cè)裝置,經(jīng)標(biāo)定傾斜誤差RMS值為1.70 μrad,位置誤差PV值為56 nm。許文斌[10]等對(duì)高架橋式三坐標(biāo)測(cè)量平臺(tái)進(jìn)行定日鏡面形精度檢測(cè),檢測(cè)精度達(dá)到0.1 mrad。這些檢測(cè)裝置均滿足目前鏡面組件的加工生產(chǎn)過(guò)程中的檢測(cè)精度要求,但曲面鏡在線定位、系統(tǒng)測(cè)量速度等問(wèn)題仍有待提高,以滿足鏡面生產(chǎn)廠商對(duì)每片反射鏡進(jìn)行在線檢測(cè)的要求。

    5 結(jié)語(yǔ)

    采用蒙特卡洛光線追跡模擬太陽(yáng)光的反射聚焦,建立鏡面組件模型,考慮太陽(yáng)張角以及鏡面面形誤差分析碟式拋物面聚光器焦斑處的能流密度分布,得出以下結(jié)論:

    (1)運(yùn)用Matlab進(jìn)行編程計(jì)算,通過(guò)與Soltrace軟件對(duì)外圈鏡面組件進(jìn)行算法的對(duì)比驗(yàn)證,仿真結(jié)果基本吻合,表明該算法準(zhǔn)確可行。

    (2)通過(guò)線性擬合以及反推計(jì)算,得出各圈鏡面組件的面形誤差允許偏差,并總結(jié)了目前幾種主要的面形誤差檢測(cè)方法,為鏡面生產(chǎn)廠商的加工制造提供檢測(cè)方法及標(biāo)準(zhǔn)。

    [1]BUIE D, MONGER A G, DEY C J. Sunshape distributions for terrestrial solar simulations[J]. Solar Energy, 2003, 74(2):113-122.

    [2]聚光型太陽(yáng)能熱發(fā)電術(shù)語(yǔ): GB/T 26972—2011 [S].

    [3]JONES S A, GRUETZNER J K, HOUSER R M, et al. VSHOT measurement uncertainty and experimental sensitivity study[C]//Proceedings of the 32th Intersociety Energy Conversion Engineering Conference. Honolulu, 1997.

    [4]劉穎, 戴景民, 孫曉剛. 拋物面型聚光器聚焦光斑能流密度分布的計(jì)算[J]. 太陽(yáng)能學(xué)報(bào), 2007, 28(10):1049-1054.

    LIU Ying, DAI Jingmin, SUN Xiaogang. Calculation method of flux distribution on focal plane in parabolic concentrators[J]. Acta Energiae Solaris Sinica , 2007, 28(10): 1049-1054.

    [5]張付行. 聚光器斜率誤差對(duì)能流分布影響的研究[D]. 長(zhǎng)沙:湖南大學(xué), 2014.

    [6]WENDELIN T, DOBOS A, LEWANDOWSKI A. SolTrace: A Ray-Tracing Code for Complex Solar Optical Systems[R]. NREL/TP-5500-59163, 2013.

    [7]肖君. 太陽(yáng)能熱發(fā)電曲面反光鏡面形檢測(cè)技術(shù)研究[D]. 北京:中國(guó)科學(xué)院大學(xué), 2015.

    [8]許文斌, 盧振武, 張紅鑫, 等. 塔式太陽(yáng)能熱發(fā)電站定日鏡面形檢測(cè)專用儀器的研制[J]. 儀器儀表學(xué)報(bào), 2009, 30(7): 1390-1394.

    XU Wenbin, LU Zhenwu, ZHANG Hongxin, et al. Development of the heliostat shape measuring apparatus of solar tower thermal power plant[J]. Chinese Journal of Scientific Instrument , 2009, 30(7): 1390-1394.

    [9]崔天剛. Wolter I型反射鏡面形檢測(cè)及成像性能檢測(cè)研究[D]. 北京:中國(guó)科學(xué)院研究生院,2011.

    [10]許文斌, 盧振武, 張紅鑫. 定日鏡面形精度檢測(cè)中的三坐標(biāo)掃描系統(tǒng)的精度分析[J]. 光學(xué)技術(shù), 2009, 35(5): 751-753.

    XU Wenbin, LU Zhenwu, ZHANG Hongxin. Accuracy analysis of coordinate measuring machine for the surface shape measurement of heliostats[J]. Optical Technique, 2009, 35(5): 751-753.

    (本文編輯:趙艷粉)

    Shape Error Effect Analysis of Dish Solar Mirror Components

    ZHANG Chen,HUANG Zhenjun,LIANG Weiqing, ZHOU Shunjun, TANG Zhixiong

    (Shanghai Micropowers Co., Ltd.,Shanghai 201203,China)

    In this paper, efficient numerical calculation function of Matlab is used to simulate the reflection and concentration of sunlight with the Monte Carlo method, realize the visualization of programming calculation results and analyze the flux density distribution at the focal plane of dish parabolic concentrator. Then simulation results show the effect analysis of incident angle, the slope error and scattering error of reflective surface on single piece of mirror component of the focal spot and prove the accuracy of the algorithm by Soltrace. Finally, this research establishes the reasonable indexes of surface shape error in accordance with the actual shape of focus spot and target size at the focal position, provides theoretical basis and testing standards for the production and processing of mirror components manufacturers as well as the surface shape error detection.

    dish parabolic concentrator; slope error; Monte Carlo method; Matlab; Soltrace

    10.11973/dlyny201703017

    張 琛(1988—),男,主要從事太陽(yáng)能光熱發(fā)電領(lǐng)域的研發(fā)工作。

    F416.61

    A

    2095-1256(2017)03-0293-05

    2017-02-16

    猜你喜歡
    檢測(cè)
    QC 檢測(cè)
    “不等式”檢測(cè)題
    “一元一次不等式”檢測(cè)題
    “一元一次不等式組”檢測(cè)題
    “幾何圖形”檢測(cè)題
    “角”檢測(cè)題
    “有理數(shù)的乘除法”檢測(cè)題
    “有理數(shù)”檢測(cè)題
    “角”檢測(cè)題
    “幾何圖形”檢測(cè)題
    黄色视频,在线免费观看| 50天的宝宝边吃奶边哭怎么回事| 精品福利观看| 91成人精品电影| 老汉色av国产亚洲站长工具| 男女下面插进去视频免费观看| 精品一区二区三区四区五区乱码| 啦啦啦观看免费观看视频高清 | 99国产精品免费福利视频| 中亚洲国语对白在线视频| 中出人妻视频一区二区| 在线国产一区二区在线| 乱人伦中国视频| 亚洲熟妇熟女久久| 99在线视频只有这里精品首页| 国产aⅴ精品一区二区三区波| 9热在线视频观看99| 欧美av亚洲av综合av国产av| 嫩草影院精品99| 久99久视频精品免费| 欧美一级a爱片免费观看看 | 欧美日韩瑟瑟在线播放| 在线观看免费视频网站a站| 亚洲九九香蕉| 女人精品久久久久毛片| 久久精品亚洲熟妇少妇任你| 日韩精品免费视频一区二区三区| 两个人视频免费观看高清| 国产亚洲欧美98| 亚洲国产欧美网| 亚洲 国产 在线| 国产午夜精品久久久久久| 一区在线观看完整版| 夜夜躁狠狠躁天天躁| 亚洲五月色婷婷综合| avwww免费| 日日干狠狠操夜夜爽| 女人被狂操c到高潮| 国产熟女xx| 午夜免费成人在线视频| 一夜夜www| 十分钟在线观看高清视频www| 一级毛片高清免费大全| 久久人人精品亚洲av| 精品久久久久久久毛片微露脸| av福利片在线| 亚洲国产日韩欧美精品在线观看 | 亚洲av电影在线进入| 69精品国产乱码久久久| 国产99久久九九免费精品| 色哟哟哟哟哟哟| 亚洲第一青青草原| 午夜精品久久久久久毛片777| 精品一区二区三区av网在线观看| 一区二区三区高清视频在线| 国产成人精品在线电影| 老司机午夜十八禁免费视频| 女人爽到高潮嗷嗷叫在线视频| 99国产极品粉嫩在线观看| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 国产亚洲精品av在线| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站 | 国产激情久久老熟女| 老熟妇乱子伦视频在线观看| 久久中文看片网| 禁无遮挡网站| 看黄色毛片网站| 亚洲 国产 在线| 91成年电影在线观看| 非洲黑人性xxxx精品又粗又长| 久久国产亚洲av麻豆专区| 午夜福利高清视频| 精品国产乱码久久久久久男人| 最新在线观看一区二区三区| 成年版毛片免费区| 国产成年人精品一区二区| 真人一进一出gif抽搐免费| 50天的宝宝边吃奶边哭怎么回事| 12—13女人毛片做爰片一| 19禁男女啪啪无遮挡网站| 女警被强在线播放| 国产熟女xx| 亚洲国产精品999在线| 欧美激情高清一区二区三区| 亚洲avbb在线观看| 亚洲人成网站在线播放欧美日韩| 中文字幕人妻丝袜一区二区| 久久香蕉激情| 成人欧美大片| 丁香欧美五月| 高清毛片免费观看视频网站| 久久久久久久久中文| 亚洲全国av大片| 午夜福利免费观看在线| 欧美另类亚洲清纯唯美| 妹子高潮喷水视频| 天天添夜夜摸| 可以免费在线观看a视频的电影网站| 亚洲va日本ⅴa欧美va伊人久久| 可以在线观看毛片的网站| 国产成人欧美| 女性生殖器流出的白浆| 午夜福利成人在线免费观看| 久久香蕉精品热| 亚洲av成人av| 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 美女免费视频网站| 国产伦人伦偷精品视频| 亚洲国产精品久久男人天堂| 一区二区三区激情视频| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 久久婷婷成人综合色麻豆| 啦啦啦 在线观看视频| 老司机午夜十八禁免费视频| 一二三四在线观看免费中文在| av视频免费观看在线观看| 十八禁人妻一区二区| 久久青草综合色| 乱人伦中国视频| 搡老妇女老女人老熟妇| 久久热在线av| 最好的美女福利视频网| 午夜福利在线观看吧| 人妻丰满熟妇av一区二区三区| 亚洲 国产 在线| 香蕉丝袜av| 成人亚洲精品av一区二区| 在线观看免费视频日本深夜| 亚洲精品国产精品久久久不卡| 国内精品久久久久久久电影| 又紧又爽又黄一区二区| 国产精品日韩av在线免费观看 | 国产一卡二卡三卡精品| 亚洲人成伊人成综合网2020| 变态另类成人亚洲欧美熟女 | 韩国av一区二区三区四区| 亚洲国产精品999在线| 欧美成人一区二区免费高清观看 | av在线天堂中文字幕| 99久久综合精品五月天人人| 一二三四在线观看免费中文在| 91精品三级在线观看| 三级毛片av免费| 欧美国产精品va在线观看不卡| 一进一出好大好爽视频| 欧美色欧美亚洲另类二区 | 午夜精品在线福利| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 91国产中文字幕| 看片在线看免费视频| 国产成人系列免费观看| 麻豆国产av国片精品| 久久精品国产综合久久久| 国产又爽黄色视频| 久99久视频精品免费| 搞女人的毛片| 久久久水蜜桃国产精品网| 国产精品亚洲av一区麻豆| 黄色成人免费大全| 国产精品av久久久久免费| 亚洲avbb在线观看| www.自偷自拍.com| 中亚洲国语对白在线视频| 精品国内亚洲2022精品成人| 午夜a级毛片| 久久九九热精品免费| 久久精品国产综合久久久| 精品一品国产午夜福利视频| 很黄的视频免费| 欧美日韩精品网址| 日韩三级视频一区二区三区| www.精华液| 老汉色∧v一级毛片| 99国产综合亚洲精品| 日本在线视频免费播放| 999久久久国产精品视频| 两个人视频免费观看高清| 国产精品一区二区免费欧美| 亚洲欧美日韩无卡精品| 男人操女人黄网站| 亚洲色图av天堂| 99在线人妻在线中文字幕| 老司机午夜福利在线观看视频| 久久香蕉精品热| 淫秽高清视频在线观看| 最近最新中文字幕大全电影3 | 国产不卡一卡二| 欧美色视频一区免费| 亚洲自偷自拍图片 自拍| 亚洲,欧美精品.| 性少妇av在线| 亚洲国产欧美日韩在线播放| 夜夜爽天天搞| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 在线观看66精品国产| 久久久久久久午夜电影| 国产精品久久电影中文字幕| 女人高潮潮喷娇喘18禁视频| 日韩大码丰满熟妇| 久久久久国产精品人妻aⅴ院| 免费高清视频大片| 午夜成年电影在线免费观看| 亚洲精品在线观看二区| 99在线视频只有这里精品首页| 日韩国内少妇激情av| 天堂动漫精品| 国产av一区在线观看免费| 神马国产精品三级电影在线观看 | 久久人妻熟女aⅴ| 成人18禁高潮啪啪吃奶动态图| 91成人精品电影| 国产色视频综合| 一本综合久久免费| 成人免费观看视频高清| 50天的宝宝边吃奶边哭怎么回事| 搡老岳熟女国产| 欧美人与性动交α欧美精品济南到| 熟妇人妻久久中文字幕3abv| 久久久国产精品麻豆| 欧美中文综合在线视频| 色在线成人网| 国产精品秋霞免费鲁丝片| 一二三四在线观看免费中文在| 精品久久久久久久毛片微露脸| 久久久久精品国产欧美久久久| 午夜免费成人在线视频| 后天国语完整版免费观看| 少妇熟女aⅴ在线视频| 啦啦啦 在线观看视频| 欧美中文综合在线视频| 久久久久亚洲av毛片大全| 久久久久国内视频| 又黄又爽又免费观看的视频| 俄罗斯特黄特色一大片| 国产成人av激情在线播放| 午夜福利欧美成人| 国产区一区二久久| 老熟妇仑乱视频hdxx| 母亲3免费完整高清在线观看| 午夜日韩欧美国产| 久久香蕉国产精品| 人人澡人人妻人| 亚洲国产看品久久| 亚洲一区高清亚洲精品| 亚洲成人久久性| 欧美一级毛片孕妇| 亚洲av成人不卡在线观看播放网| 给我免费播放毛片高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| 一区二区三区国产精品乱码| 亚洲精品美女久久久久99蜜臀| 自拍欧美九色日韩亚洲蝌蚪91| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 亚洲国产日韩欧美精品在线观看 | 别揉我奶头~嗯~啊~动态视频| 一边摸一边抽搐一进一出视频| 美女高潮到喷水免费观看| 日本在线视频免费播放| 久久狼人影院| www国产在线视频色| 午夜精品国产一区二区电影| 好男人在线观看高清免费视频 | 国产蜜桃级精品一区二区三区| 欧美黄色淫秽网站| 国产男靠女视频免费网站| 后天国语完整版免费观看| 亚洲色图av天堂| 丝袜美腿诱惑在线| 人妻丰满熟妇av一区二区三区| 国产高清激情床上av| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| a级毛片在线看网站| 一本大道久久a久久精品| 19禁男女啪啪无遮挡网站| 国产人伦9x9x在线观看| 精品一品国产午夜福利视频| 欧美成人免费av一区二区三区| 久久久久亚洲av毛片大全| 真人一进一出gif抽搐免费| 亚洲专区国产一区二区| 国产精品 国内视频| 又黄又粗又硬又大视频| 久久中文字幕一级| 波多野结衣巨乳人妻| 亚洲男人天堂网一区| 日韩欧美在线二视频| 波多野结衣av一区二区av| 99久久国产精品久久久| 午夜视频精品福利| 中文字幕人妻熟女乱码| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 亚洲专区国产一区二区| 中亚洲国语对白在线视频| 亚洲va日本ⅴa欧美va伊人久久| 国产成人精品在线电影| 黄色毛片三级朝国网站| 啦啦啦 在线观看视频| 九色国产91popny在线| 久久精品国产亚洲av香蕉五月| 亚洲一区高清亚洲精品| 久久久国产欧美日韩av| 女人被狂操c到高潮| 亚洲av成人不卡在线观看播放网| 99久久国产精品久久久| 嫩草影视91久久| 欧美色视频一区免费| 中文亚洲av片在线观看爽| 免费人成视频x8x8入口观看| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看 | 欧美亚洲日本最大视频资源| 午夜免费激情av| 免费无遮挡裸体视频| 亚洲欧美一区二区三区黑人| 日韩国内少妇激情av| av在线天堂中文字幕| 精品高清国产在线一区| 亚洲av电影在线进入| 久久 成人 亚洲| 天堂动漫精品| 欧美最黄视频在线播放免费| 长腿黑丝高跟| 中文字幕久久专区| 国产一区二区在线av高清观看| 久久久久久久久免费视频了| 亚洲久久久国产精品| 九色国产91popny在线| 精品国产美女av久久久久小说| 国产精品久久电影中文字幕| 男女做爰动态图高潮gif福利片 | 国产又爽黄色视频| 欧美日韩亚洲国产一区二区在线观看| av欧美777| 亚洲成国产人片在线观看| 男女午夜视频在线观看| 女人被狂操c到高潮| 人人妻人人澡欧美一区二区 | 国产不卡一卡二| 亚洲欧美精品综合久久99| 成人三级黄色视频| 人妻久久中文字幕网| 亚洲国产毛片av蜜桃av| 成人三级做爰电影| 久久狼人影院| 国产精品 欧美亚洲| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美精品永久| 午夜a级毛片| 校园春色视频在线观看| av欧美777| 免费久久久久久久精品成人欧美视频| 午夜福利一区二区在线看| 午夜久久久在线观看| av在线播放免费不卡| 中出人妻视频一区二区| 国产精品香港三级国产av潘金莲| 亚洲avbb在线观看| 国产欧美日韩一区二区三| 黄色丝袜av网址大全| 少妇被粗大的猛进出69影院| 国产片内射在线| 久久天堂一区二区三区四区| 久久热在线av| 亚洲熟妇熟女久久| 色播在线永久视频| 国产亚洲精品久久久久5区| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| 又黄又爽又免费观看的视频| 久久精品91无色码中文字幕| 女同久久另类99精品国产91| tocl精华| 欧美一区二区精品小视频在线| 操出白浆在线播放| 涩涩av久久男人的天堂| 欧美成人免费av一区二区三区| 久久狼人影院| 亚洲天堂国产精品一区在线| 亚洲久久久国产精品| 脱女人内裤的视频| 免费高清在线观看日韩| 不卡av一区二区三区| а√天堂www在线а√下载| www日本在线高清视频| 日韩av在线大香蕉| 中文亚洲av片在线观看爽| 青草久久国产| 亚洲 国产 在线| 性少妇av在线| 精品国产一区二区三区四区第35| 国产片内射在线| 免费在线观看完整版高清| 午夜免费鲁丝| 国产一区二区激情短视频| 色播在线永久视频| 亚洲七黄色美女视频| 日韩精品中文字幕看吧| 18禁裸乳无遮挡免费网站照片 | 亚洲成人精品中文字幕电影| 免费久久久久久久精品成人欧美视频| 亚洲成人国产一区在线观看| 亚洲色图综合在线观看| av在线天堂中文字幕| 99国产精品一区二区三区| 国产成人精品久久二区二区91| 亚洲一区中文字幕在线| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 热99re8久久精品国产| 午夜影院日韩av| 大码成人一级视频| 久久天堂一区二区三区四区| 99精品在免费线老司机午夜| 嫩草影院精品99| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美av亚洲av综合av国产av| 久久伊人香网站| 免费在线观看完整版高清| 精品国内亚洲2022精品成人| 亚洲性夜色夜夜综合| 国产免费男女视频| a级毛片在线看网站| 久久精品91蜜桃| av福利片在线| 亚洲国产毛片av蜜桃av| 黄色丝袜av网址大全| 日本精品一区二区三区蜜桃| 色播在线永久视频| 亚洲aⅴ乱码一区二区在线播放 | 久久人人97超碰香蕉20202| 日本vs欧美在线观看视频| 亚洲自拍偷在线| 欧美中文综合在线视频| 99国产精品免费福利视频| 成人精品一区二区免费| 满18在线观看网站| 亚洲电影在线观看av| av福利片在线| 午夜福利18| 多毛熟女@视频| 国产一卡二卡三卡精品| 精品久久久久久久人妻蜜臀av | 国产成人免费无遮挡视频| 国产精品av久久久久免费| 亚洲av五月六月丁香网| www.熟女人妻精品国产| 深夜精品福利| 精品国产国语对白av| 午夜亚洲福利在线播放| 欧美成人午夜精品| 国产成人av教育| 97人妻精品一区二区三区麻豆 | 国内精品久久久久精免费| 91成年电影在线观看| 色婷婷久久久亚洲欧美| 久久性视频一级片| 日本五十路高清| 美女扒开内裤让男人捅视频| 中文字幕人妻丝袜一区二区| 欧美日本视频| 久久久久久久久久久久大奶| 免费在线观看影片大全网站| 亚洲成a人片在线一区二区| 在线视频色国产色| 国产aⅴ精品一区二区三区波| 久久久久久国产a免费观看| 国产免费男女视频| 大香蕉久久成人网| 国产av在哪里看| 亚洲国产精品久久男人天堂| 国产一区二区三区综合在线观看| 老司机深夜福利视频在线观看| 这个男人来自地球电影免费观看| 嫩草影视91久久| 老司机福利观看| 欧美黑人欧美精品刺激| 亚洲一区二区三区不卡视频| 日韩高清综合在线| 韩国精品一区二区三区| 法律面前人人平等表现在哪些方面| www.999成人在线观看| 午夜日韩欧美国产| 一级,二级,三级黄色视频| 少妇粗大呻吟视频| 老司机福利观看| 88av欧美| 又黄又粗又硬又大视频| 午夜福利视频1000在线观看 | 成人18禁在线播放| 成人免费观看视频高清| 国产乱人伦免费视频| 久久香蕉激情| 中文字幕人成人乱码亚洲影| 免费观看精品视频网站| 性色av乱码一区二区三区2| 亚洲五月色婷婷综合| 国产成人一区二区三区免费视频网站| avwww免费| 脱女人内裤的视频| 亚洲一码二码三码区别大吗| 精品一区二区三区视频在线观看免费| 国产97色在线日韩免费| 精品乱码久久久久久99久播| 精品久久久精品久久久| 久久久久国内视频| 午夜免费成人在线视频| 国产精品一区二区在线不卡| 免费高清在线观看日韩| 亚洲欧美日韩另类电影网站| 很黄的视频免费| 他把我摸到了高潮在线观看| 看黄色毛片网站| 一级,二级,三级黄色视频| 久久人人97超碰香蕉20202| 精品久久久久久成人av| 国产亚洲欧美在线一区二区| 欧美激情久久久久久爽电影 | 国产精品久久久av美女十八| 欧美一级a爱片免费观看看 | 黄色视频,在线免费观看| 亚洲五月色婷婷综合| 日韩一卡2卡3卡4卡2021年| 他把我摸到了高潮在线观看| 欧美日本视频| 午夜福利视频1000在线观看 | 最近最新中文字幕大全电影3 | 丝袜人妻中文字幕| 女人高潮潮喷娇喘18禁视频| 久久国产精品人妻蜜桃| 国产伦人伦偷精品视频| 国产成人欧美| 免费久久久久久久精品成人欧美视频| av有码第一页| 一级黄色大片毛片| 夜夜夜夜夜久久久久| 亚洲男人天堂网一区| 欧美性长视频在线观看| 在线观看日韩欧美| 男女做爰动态图高潮gif福利片 | 久久久水蜜桃国产精品网| 色老头精品视频在线观看| 久99久视频精品免费| 欧美黑人精品巨大| 99久久99久久久精品蜜桃| a在线观看视频网站| 亚洲自拍偷在线| 国内精品久久久久精免费| 亚洲一码二码三码区别大吗| 天天躁夜夜躁狠狠躁躁| 国产av一区在线观看免费| av在线天堂中文字幕| 精品一区二区三区视频在线观看免费| 中出人妻视频一区二区| 国产精品九九99| 国产精品1区2区在线观看.| 亚洲欧美激情综合另类| 午夜福利,免费看| 少妇熟女aⅴ在线视频| 在线观看免费午夜福利视频| 一区二区日韩欧美中文字幕| 亚洲午夜理论影院| 国产色视频综合| 午夜福利欧美成人| 亚洲中文日韩欧美视频| 免费看十八禁软件| 久久国产乱子伦精品免费另类| 人妻丰满熟妇av一区二区三区| cao死你这个sao货| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 搡老妇女老女人老熟妇| 日本免费一区二区三区高清不卡 | 视频区欧美日本亚洲| 在线观看免费视频日本深夜| www国产在线视频色| 女生性感内裤真人,穿戴方法视频| 一二三四在线观看免费中文在| 桃色一区二区三区在线观看| 男女床上黄色一级片免费看| 亚洲国产精品合色在线| 久久久精品国产亚洲av高清涩受| 亚洲,欧美精品.| 岛国在线观看网站| x7x7x7水蜜桃| 99国产精品免费福利视频| 精品福利观看| 十分钟在线观看高清视频www| 亚洲熟妇熟女久久| 麻豆成人av在线观看| 成年女人毛片免费观看观看9| 国产亚洲欧美在线一区二区| 亚洲欧美激情在线| 成熟少妇高潮喷水视频| 精品久久久久久久毛片微露脸| 亚洲中文字幕一区二区三区有码在线看 | 日本免费一区二区三区高清不卡 | 成在线人永久免费视频| 欧美亚洲日本最大视频资源| 国内久久婷婷六月综合欲色啪| 制服丝袜大香蕉在线| 纯流量卡能插随身wifi吗| 91大片在线观看| 欧美在线一区亚洲| 人人妻人人澡欧美一区二区 |