999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

高熱流散熱泵驅兩相流體回路設計及飛行驗證

2017-03-27 06:16:54于新剛苗建印王德偉范含林菅魯京張紅星
宇航學報 2017年2期
關鍵詞:系統

于新剛,徐 侃,苗建印,王德偉,趙 欣,范含林,菅魯京,張紅星

(1.北京空間飛行器總體設計部,北京100094;2.空間熱控技術北京市重點實驗室,北京100094)

高熱流散熱泵驅兩相流體回路設計及飛行驗證

于新剛1,2,徐 侃1,2,苗建印1,2,王德偉1,2,趙 欣1,2,范含林1,菅魯京1,2,張紅星1,2

(1.北京空間飛行器總體設計部,北京100094;2.空間熱控技術北京市重點實驗室,北京100094)

針對激光、雷達等載荷空間應用中存在的超過100 W/cm2高熱流密度散熱問題,采用泵驅兩相流體回路結合金剛石微槽道的方法實現了271 W/cm2熱流密度的散熱。對系統設計及關鍵部件設計進行了研究,整個系統經過地面和在軌飛行驗證,運行平穩。利用儲液器控溫的方式,實現了系統溫度的準確控制。測試結果表明,地面和在軌飛行期間的換熱及流阻相當,驗證了重力無關設計的正確性。

微槽道;兩相流體回路;流動沸騰

0 引 言

隨著技術的發展,雖然電子設備在材料和性能方面不斷改進,但是對大功率、高通量的電子設備的散熱、溫度控制還是變得越來越困難[1-3]。在宇航領域,隨著電子設備的功耗增長,散熱技術由被動熱控技術發展至風冷直至單相液冷,目前宇航領域對于大功率設備散熱最有效的手段仍是單相流體回路技術[4-6],以及以環路熱管為代表的熱管技術。如向艷超等[7]采用重力輔助兩相流體回路解決CE-3月夜問題。近年來,以空間用激光載荷、相控陣雷達為代表的大功率設備的熱流密度已經達到或超過了100 W/cm2,采用傳統的單相流體換熱的方式解決其散熱問題難度很大[8]。這些需求促使研究者發展高性能的兩相換熱方式[9]。熱流密度超過100 W/cm2時,基于相變的流動沸騰換熱(微槽道)、兩相噴霧冷卻、兩相射流冷是常用的幾種散熱方式,這三種換熱方式在換熱熱流密度和換熱溫差上同時滿足實際的工程應用[10]。表1[11]給出了幾種散熱方式的比較,可以看出在系統體積和溫度均勻性上微槽道流動換熱具有優勢,而從換熱的性能來看射流和噴霧冷卻具有優勢。在換熱能力均滿足需求的情況下,采用噴霧冷卻和射流冷卻,系統體積較大,這給空間應用帶來了困難,隨著微加工技術的發展,結構緊湊、輕便、溫度均勻性好的微槽道散熱方式,在未來宇航領域高熱流密度散熱領域應用前景廣闊[12]。

微槽道內流動沸騰是在 1981年第一次被Tuckerman和Pease[13]驗證的。在接下來的三十多年,微通道流動換熱得到了廣泛關注。大量的研究集中在[13-22]以下四個方面:1)微小流動通道里的流動沸騰換熱特性和傳熱機理;2)微槽道內流態的差異以及絕熱和非絕熱流動的特性研究;3)數值模擬以及對流動換熱系數的研究;4)氣泡生長的模擬,不穩定性的克服以及流動的強化。這些研究為微通道流動沸騰傳熱特性的理解提供了指導。

上述的微槽道內流動沸騰換熱的研究集中在地面重力條件下,對于微重力條件下由于氣泡生長、分布與地面不同所引起的微通道內流動沸騰換熱的差異研究還比較少。

Konishi等[23]綜述了微重力流動沸騰換熱的研究,著重從理論上對微重力流動沸騰的特殊性進行分析,包括加熱方向不同對臨界熱流密度的影響,池沸騰和流動沸騰在微重力下的差異以及與重力無關的設計方法,并采用落塔、拋物飛機等開展了少量試驗,但由于試驗條件限制以及試驗時間很短(一般只有幾秒),并未真正實現穩定的微重力條件下的流動沸騰換熱,其結果能否用于空間兩相散熱系統還存在疑問。

表1 不同散熱方式的比較Table 1 Comparison of different cooling scheme

微槽道作為流動沸騰換熱熱沉,需要與泵驅兩相流體回路相結合,即作為泵驅兩相流體回路系統的蒸發器??臻g用泵驅兩相流體回路的研究開始于20世紀80年代,主要經歷了以下幾個階段[24-25]:

1)20世紀80年代,Oren、Stalmach、Haslett開始開展泵驅兩相回路基礎理論研究。

2)1994年重新設計前,國際空間站美國部分、俄羅斯部分采用泵驅兩相回路,后由于熱源數量減少,改為泵驅單相流體回路。

3)1997年 STS-85任務搭載的 NASDA和TOSHIBA研制的TPFLEX為泵驅兩相流體回路的首次搭載。

4)2011年發射的AMS02泵驅兩相流體回路系統,為泵驅兩相流體回路系統的首次在軌應用。

在此之后,空間用泵驅兩相流體回路也被用于解決大型天線以及深空探測任務的飛行器的熱控制[26-28]。

上述空間應用中泵驅兩相流體回路主要解決熱量遠距離傳輸和溫度均勻性問題,超過100 W/cm2空間極高熱流散熱問題的研究還未看到報道。

綜上所述,現有的用于空間極高熱流散熱的微槽道流動沸騰換熱的研究,特別是重力無關設計方法和機理還缺乏穩定的空間飛行驗證,而目前用于航天器熱控的泵驅兩相流體回路系統還未達到極高熱流的范疇,因此為解決熱流密度超過100 W/cm2的空間載荷的散熱問題,需要將二者結合起來,開展進一步研究工作。

本文結合金剛石微槽道散熱技術及泵驅兩相流體回路,解決空間極高熱流密度問題,并對重力無關設計方法進行驗證。

1 兩相流體回路試驗系統

1.1 系統流程

整個泵驅兩相流體回路流程如圖1所示,由蒸發器、預熱器、回熱器、機械泵、儲液器、冷凝器、壓力傳感器和傳輸管路構成氨工質閉合回路。蒸發器內處于氣液兩相狀態的工質,從高熱流密度加熱器處吸收熱量后產生沸騰,干度增大,接著流經回熱器把熱量傳給從泵出來的過冷液,干度降低,再通過相變裝置把熱量排散給相變材料,干度進一步降低,達到過冷狀態。從相變裝置流出的過冷液再順序通過泵、回熱器,經過預熱器加熱至兩相狀態后回流到蒸發器中。工質在蒸發器和相變裝置間的循環流動由泵提供的動力來維持。由于整個系統沒有外部的散熱熱沉,因此采用相變材料作為散熱熱沉。參加飛行試驗的泵驅兩相流體回路系統如圖2所示。

1.2 工質

在目前的設計中采用無水氨作為回路的循環工質。在常用的循環工質中氨的潛熱僅次于水(常壓下分別為1368和2257 kJ/kg。此外氨的三相點為-77℃、臨界溫度為132℃,可以滿足空間應用的需求,也是目前宇航用熱管和環路熱管的常用工質。

1.3 微槽道蒸發器

蒸發器和微槽道的結構如圖3所示,蒸發器由帶有微槽結構的金剛石薄片焊接在金屬本體上組成,微槽的設計寬度為0.2 mm,深度為1 mm。采用的金剛石的導熱系數約為1500 W/m℃,為了減小系統阻力,對微槽道表面進行了處理,處理前后對比如圖4所示。采用金剛石主要是減小蒸發器自身的導熱熱阻,圖5給出了蒸發器上表面施加500 W的熱量,維持下表面溫度為20℃的計算結果,材料為金剛石時,溫差只有7℃。

1.4 循環泵

兩相循環泵是泵驅兩相流體回路的核心,機械泵應用在兩相熱控系統中,與常用的單相泵類產品不同,具有以下特殊性:1)在工質飽和蒸汽壓下工作,耐壓能力強;2)入口工質過冷度較小時,需要具有一定驅動兩相工質的能力;3)泵流量很小,揚程較大。基于上述需求,設計一種軸流式螺紋驅動泵作為流體回路的驅動源。

1.5 高熱流加熱器

模擬超過100 W/cm2的熱源是開展試驗研究的一個關鍵環節。設計中采用了變截面的方式來放大熱流密度,高熱流密度加熱器的外形尺寸如圖6所示,本體材料為紫銅,其頂面上焊接了4路規格相同的4路加熱片,最大功率可達500 W,通過包覆的氣凝膠來隔熱和保溫,加熱器與蒸發器的耦合面的大小為1 cm2,最大可模擬500 W/cm2的熱耗。

1.6 系統控制

整個泵驅兩相流體回路的溫度控制通過儲液器來實現,儲液器內的工質為氣液共存的狀態,通過改變儲液器的溫度,改變系統的壓力,進而控制蒸發和冷凝的溫度,實現對系統的溫度控制。

2 地面測試結果及分析

系統不加熱負荷的情況下運轉,測量了不同轉速條件下,系統的流量數據,作為確定泵工作參數的依據。根據循環泵的測試結果,選擇5000 r/min作為飛行試驗時泵的工作轉速。

表2 系統流量測試結果Table 2 Flow rate of fluid loop

在泵轉速5000 r/min的情況下,地面對系統的運行進行了測試,高熱流密度加熱器分別開1路、2路、3路,測試結果如表3所示,系統工作穩定、正常。從測試結果可以看出:

1)利用儲液器控溫精度較高,儲液器與蒸發器出口溫度差異小于0.2℃,系統處于良好的兩相工作狀態。

2)系統的流阻主要存在于蒸發器兩端,約為60 kPa,而且這一流阻對系統熱流密度(加熱量)不敏感,當加熱量從100.7 W~271.1 W變化時,流阻相當。

3)蒸發器的沸騰換熱系數在7.8~9.1 W/ cm2℃之間。

4)相變裝置與工質之間換熱良好,相變裝置出口工質的溫度與相變裝置側壁的溫差小于2℃。

圖7給出了地面測試時儲液器溫度、蒸發器出入口溫度、泵入口壓力隨時間的變化。從圖7可以看出,在高熱流密度加熱器打開后在蒸發器出口出現一個約4℃的溫度過沖,分析是由于加熱的熱流密度很高,在形成穩定流動沸騰換熱前溫度出現過沖,進入兩相狀態后溫度下降,在后續加熱功率增加時未出現類似的溫度過沖現象。

工作過程中系統的壓力保持穩定,未出現大的壓力波動,只是隨著儲液器的溫度波動而小幅波動。由于泵驅兩相流體回路中循環泵的能力較強,因此小的壓力波動未對系統的工作產生不良影響。

儲液器溫度不斷上升是由于儲液器控溫的方式采用隨動控溫的方式,其控溫的目標是冷凝器出口溫度加上一個增量,由于隨著試驗進行相變裝置溫度不斷上升所以儲液器溫度不斷上升,當相變裝置溫度趨于穩定時儲液器溫度也趨于穩定。

3 在軌飛行結果及其與地面試驗結果比對

2016年6月26日,高熱流散熱泵驅兩相流體回路搭載多功能飛船縮比返回艙進行了在軌飛行。飛行試驗在10∶41開始,13∶41結束,共持續了3個小時,熱源的最高熱流密度為271 W/cm2。

圖8給出了飛行期間關鍵點的溫度、壓力變化。從飛行試驗數據來看,在軌飛行期間溫度和壓力的變化情況與地面試驗類似,通過儲液器控制了蒸發器溫度,系統工作狀態良好。試驗的整體溫度較地面測試低,這主要是由于入軌后,系統的初始溫度較低的緣故。飛行試驗中由于相變材料初始溫度較低(約11℃),而且試驗開展的時間有限,因此相變材料未達到相變點即結束了試驗,圖8中溫度未趨于穩定值。

從試驗過程來看,在高熱流密度加熱器開始加熱后出現了一個1.5℃的溫度過沖,這一過沖較地面試驗時的4℃有所減小,這是由于在高熱流施加的瞬間,蒸發器(見圖9)內還未建立穩定的流動沸騰換熱,因此加熱面溫度瞬間升高,地面試驗時由于氣泡浮力作用難以脫離加熱面,建立穩定流動沸騰換熱時間相對較長,而在軌情況下由于氣泡迅速脫離縮短了建立穩定的流動沸騰換熱的時間,使得加熱面溫度在這一過程中升高較少。

在軌飛行過程中,根據高熱流密度加熱器與工質的溫差推算得到了流動沸騰換熱系數與地面測試結果相當,差別在10%以內,這是由于設計時進行了重力無關設計,根據文獻[23],流速大于1.5 m/s時,流動沸騰換熱受重力影響較小,系統設計時通過流量和槽道尺寸的設計使得槽道內的流速達到2.5 m/s,因此進入了重力無關區域。此外由于進行了重力無關設計在軌飛行期間和地面測試期間系統的流阻相當,均為60 kPa左右。這一結果對于微槽道散熱技術在宇航領域的應用具有顯著的意義。兩相系統應用于空間微重力或其他特殊重力條件下一個主要的障礙就是由于氣泡的浮升力的影響,其在地面和空間行為差異很大,難以在地面進行充分的驗證,因此兩相系統在宇航領域的應用一直比較慎重,需要預留大量的設計余量,或者通過非常復雜的地面試驗方法進行驗證,顯著的制約了其應用,通過此次在軌試驗和地面數據的對比充分驗證了采用合理的重力無關設計,地面試驗可以對在軌工況的流動和換熱特性進行充分的驗證,減少飛行試驗的不確定性。

高熱流密度蒸發器上采用氣凝膠覆蓋的方式進行隔熱,飛行過程中氣凝膠兩側的溫差約為100℃,地面測試時氣凝膠兩側的溫差約為97℃,當量的導熱系數約在10-3W/m℃量級,在軌期間由于氣凝膠中的氣體排出其隔熱效果進一步改善,變化量約為3%,差異不大。

表3 工作性能測試結果Table 3 Performance of fluid loop in terrestrial test

4 結 論

本文利用金剛石微槽道結合泵驅兩相流體回路實現了271 W/cm2熱流密度的散熱,流動沸騰換熱系數達到7.8~9.1 W/cm2℃。在軌和地面數據比對驗證了控制流速的重力無關設計方法的有效性,為后續空間應用奠定了基礎。

[1] Price D C.A review of selected thermal management solutions for military electronic systems[J].IEEE Trans.Compon.Packag.Tech.,2003,26(1):26-39.

[2] Krishnan S,Garimella S V,Chrysler G,et al.Towards a thermal Moore's law[J].IEEE Trans.Adv.Packag.,2007,30(3):462-474.

[3] Mudawar I. Two-phase microchannelheat sinks: theory,applications,and limitations[J]. ASME J. Electron.Packaging,2011,133(4):041002.

[4] 張立,姜軍,趙啓偉.單相流體回路在衛星熱控制中的應用研究[J].中國空間科學技術,2008,28(3):65-71.[Zhang Li,Jiang Jun,Zhao Qi-wei.Application of single-phase fluid loop system on thermal control of satellite[J].Chinese Space Science and Technology,2008,28(3):65-71.]

[5] 范宇峰,黃家榮.熱控流體回路補償器的熱計算方法[J].宇航學報,2008,29(3):826-830.[Fan Yu-feng,Huang Jiarong,Thermal calculation method for the accumulator of the liquid fluid loop in thermal control system[J].Journal of Astronautics,2008,29(3):826-830.]

[6] 吳自帥,黃臻成,莫冬傳,等.多艙耦合流體系統在功耗調配模式下的模擬分析[J].宇航學報,2012,33(12):1860-1865.[Wu Zi-shuai,Huang Zhen-cheng,Mo Dong-chuan,et al.Simulation analysis for a multi-compartment thermal coupled system under power Allocation condition[J].Journal of Astronautics,2012,33(12):1860-1865.]

[7] 向艷超,陳建新,張冰強,等.嫦娥三號“玉兔"巡視器熱控制[J].宇航學報,2015,36(10):1203-1209.[Xiang Yanchao,Chen Jian-xin,Zhang Bing-qiang,et al.Thermal control for Jade Rabbit Rover of Chang’E-3[J].Journal of Astronautics,2015,36(10):1203-1209.]

[8] Anderson T M,Mudawar I.Microelectronic cooling by enhanced pool boiling of a dielectric fluorocarbon liquid[J].ASME J.Heat Transfer,1989,111(3):752-759.

[9] Estes K A,Mudawar I.Comparison of two-phase electronic cooling using free jets and sprays[J].Journal of Electronic Packaging,1995,117(6):323-332.

[10] Mudawar I.Recent advances in high-flux,two-phase thermal management[J].Journal of Thermal Science and Engineering Applications,2013,5(2):021012.

[11] Kandlikar S G.History,advances,and challenges in liquid flow and flow boiling heat transfer in microchannels:a critical review[J].Journal of Heat Transfer,2012,134(3):034001.

[12] 邵保平.星載微型電子系統微尺度傳熱與散熱方法研究[D].哈爾濱:哈爾濱工業大學,2007.[Shao Bao-ping.Study on micro-scale heat transfer and dissipation method for space bome micro-electronics system[D].Harbin;Harbin Institute of Technology,2007.]

[13] Tuckerman D B,Pease R F W.High-performance heat sinking for VLSI[J].Electron Device Letters,IEEE,1981,2(5):126-129.

[14] Park J E,Thome J R.Critical heat flux in multi-microchannel copper elements with low pressure refrigerants[J].Int.J.Heat Mass Transfer,2010,53(1):110-122.

[15] Mudawar I. Two-phase microchannelheatsinks: theory,applications,and limitations[J]. JournalofElectronic Packaging,2011,133(4):041002.

[16] Chen T L,Garimella S V.A study of critical heat flux during flow boiling in microchannel heat sinks[J].Journal of Heat Transfer,2012,134(1):011504.

[17] Kandlikar S G,Widger T T,Kalani A,et al.Enhanced flow boiling over open microchannels with uniform and tapered gap manifolds[J].Journal of Heat Transfer,2013,135(6):061401.

[18] Miner M J,Phelan P E,Odom B A,et al.Optimized expanding microchannel geometry for flow boiling[J].Journal of Heat Transfer,2013,135(4):042901.

[19] Kadam S T,Kumar R.Twenty first century cooling solution: microchannel heat sinks[J].International Journal of Thermal Science,2014,85:73-92.

[20] Adhama A M,GhazaliN M, Ahmad R. Thermaland hydrodynamic analysis of microchannel heatsinks:a review[J].Renewable and Sustainable Energy Reviews,2013,21:614-622.

[21] Prajapati Y K,Pathak M,Kaleem K M.A comparative study of flow boiling heat transfer in three different configurations of microchannels[J].International Journal of Heat and Mass Transfer,2015,85:711-722.

[22] Emran M,Islam M A.Numerical investigation of flow dynamics and heat transfer characteristics in a micro-channel heat sink[J].Procedia Engineering,2014,90:563-568.

[23] Konishi C,Mudawar I.Review of flow boiling and critical heat flux in microgravity[J].International Journal of Heat and Mass Transfer,2015,80(4):469-493.

[24] 孫西輝.機械泵驅CO2兩相流體回路穩定性研究[D].廣州:中山大學,2010.[Sun Xi-hui.Stabilities research of the mechanically pumped two-phase CO2loop[D].Guangzhou:Sun Yat-sen University,2010.]

[25] Delil A A M.Research issues on two-phase loops for space applications[R].Holand:National Aerospace Laboratory,2000.

[26] Armananzasi I S,Hugon J,Ravasso N,et al.Development of a two-phase mechanically pumped loop(2ΦMPL)for the thermal dissipation management of an active antenna:experimental onground results[C].The 43rd International Conference on Environmental Systems,Val.,Co.,July 14-18,2013.

[27] Sunada E,Benjaman F,Bhandari P,et al.A two-phase mechanically pumped fluid loop for thermal control of deep space science missions[C].The 46th International Conference on Environmental Systems,Vienna,Austria,July 10-14,2016.

[28] Elli M C,Kurwitz R C.Development of a pumped two-phase system for spacecraft thermal control[C].The 46th International Conference on Environmental Systems,Vienna,Austria,July 10-14,2016.

通信地址:北京5142信箱364分箱(100094)

電話:(010)68744734

E-mail:27791802@qq.com

(編輯:牛苗苗)

Design and On-Board Validation of Pumped Two-Phase Fluid Loop for High Heat Flux Removal

YU Xin-gang1,2,XU Kan1,2,MIAO Jian-yin1,2,WANG De-wei1,2,ZHAO Xin1,2,FAN Han-lin1,JIAN Lu-jing1,2,ZHANG Hong-xing1,2
(1.Beijing Institute of Spacecraft System Engineering,Beijing 100094,China; 2.Beijing Key Laboratory of Space Thermal Control Technology,Beijing 100094,China)

A cooling system consisting of a pumped two-phase fluid loop and a diamond micro-channel is illustrated to meet the requirement of heat removal for high power payload in space,such as laser and radar.Detailed design of this system with a cooling capability as high as 271 W/cm2is presented,which is verified by the on-board and terrestrial experiments.The system works very well with a two-phase accumulator to control the temperature of the micro-channel heat sinks.The similarity between the on-board and terrestrial results of the heat transfer and flow demonstrates the gravityindependent design which is critical for space applications.

Micro-channel;Two-phase fluid loop;Flow boiling

V416

A

1000-1328(2017)02-0192-06

10.3873/j.issn.1000-1328.2017.02.011

于新剛(1981-),男,博士,主要從事航天器熱控方面的研究。

2016-09-26;

2016-11-27

國家自然科學基金(51406009)

猜你喜歡
系統
Smartflower POP 一體式光伏系統
工業設計(2022年8期)2022-09-09 07:43:20
WJ-700無人機系統
ZC系列無人機遙感系統
北京測繪(2020年12期)2020-12-29 01:33:58
基于PowerPC+FPGA顯示系統
基于UG的發射箱自動化虛擬裝配系統開發
半沸制皂系統(下)
FAO系統特有功能分析及互聯互通探討
連通與提升系統的最后一塊拼圖 Audiolab 傲立 M-DAC mini
一德系統 德行天下
PLC在多段調速系統中的應用
主站蜘蛛池模板: 亚洲人成电影在线播放| 婷婷亚洲最大| 精品第一国产综合精品Aⅴ| 97精品国产高清久久久久蜜芽| 国产丝袜啪啪| 国产免费高清无需播放器| 色婷婷在线影院| 欧美伦理一区| 国产成人精品2021欧美日韩| 国产男女XX00免费观看| 国产成+人+综合+亚洲欧美| 色一情一乱一伦一区二区三区小说 | 久久久久亚洲av成人网人人软件| 四虎永久在线精品影院| 91精品啪在线观看国产91九色| 国产成人亚洲精品蜜芽影院| 亚洲色大成网站www国产| 波多野结衣一区二区三区四区视频| 久久夜色精品| 国产黄色视频综合| 伊人成人在线| 久久一日本道色综合久久| 毛片三级在线观看| 视频二区亚洲精品| 亚洲色图欧美激情| 久久久久久午夜精品| 免费国产不卡午夜福在线观看| 亚洲成肉网| 少妇人妻无码首页| 欧美精品H在线播放| 丁香婷婷激情网| 国产精品蜜臀| 久久国产高清视频| a欧美在线| 青青操国产视频| 亚洲乱码在线视频| 国产区成人精品视频| 国产传媒一区二区三区四区五区| 91久久偷偷做嫩草影院| 国产午夜一级毛片| 中文字幕2区| 国产精品偷伦在线观看| 内射人妻无码色AV天堂| 91人人妻人人做人人爽男同 | 日日噜噜夜夜狠狠视频| 欧美不卡在线视频| 一级黄色网站在线免费看| 成人av专区精品无码国产 | 色天天综合| 99精品国产自在现线观看| 亚洲二三区| 久久黄色毛片| 久久成人18免费| 日韩毛片免费| 毛片在线播放a| 欧美成人看片一区二区三区| 国产女人爽到高潮的免费视频| 91福利片| 亚洲欧美h| 成人免费午间影院在线观看| 久久永久视频| 午夜人性色福利无码视频在线观看| 欧美亚洲国产视频| 亚洲品质国产精品无码| 91久久偷偷做嫩草影院| 日本尹人综合香蕉在线观看| 日韩精品一区二区三区视频免费看| 国产va免费精品观看| 首页亚洲国产丝袜长腿综合| 99re视频在线| 亚洲精品不卡午夜精品| 日韩人妻无码制服丝袜视频| 四虎成人精品| 欧美成人在线免费| 精品少妇人妻av无码久久| 青青草原国产免费av观看| 狠狠干综合| 色妞www精品视频一级下载| 国产福利微拍精品一区二区| 国产欧美高清| 国产精品第一区| 色婷婷久久|