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

    直升機(jī)熱管理與紅外輻射特性耦合分析方法

    2021-03-27 04:48:50李言青宣益民
    航空學(xué)報(bào) 2021年3期
    關(guān)鍵詞:換熱器直升機(jī)燃油

    李言青,宣益民

    南京航空航天大學(xué) 能源與動(dòng)力學(xué)院,南京 210016

    直升機(jī)具有機(jī)動(dòng)靈活、快速機(jī)動(dòng)的優(yōu)點(diǎn),在戰(zhàn)爭(zhēng)中得到越來越廣泛的應(yīng)用。在未來軍事作戰(zhàn)中要求直升機(jī)具有巡航速度更快、續(xù)航能力更強(qiáng)、可靠性更高等性能,這就要求直升機(jī)在攜帶有限燃料的前提下,提高整機(jī)能量利用率,實(shí)現(xiàn)能量的綜合高效利用。同時(shí),未來直升機(jī)將裝備大功率發(fā)動(dòng)機(jī)和電子設(shè)備,這加劇了直升機(jī)系統(tǒng)和整機(jī)的熱管理問題。而紅外抑制和復(fù)合材料的機(jī)身設(shè)計(jì),使得機(jī)身重量和散熱能力受到了極大的限制。因此,面對(duì)未來直升機(jī)的熱管理和紅外抑制需求,從整機(jī)層面統(tǒng)籌管理整機(jī)熱量,合理匹配熱沉梯級(jí)利用方法,在提高能量利用率的同時(shí),保障各子系統(tǒng)安全工作,并降低整機(jī)紅外輻射特性是未來直升機(jī)技術(shù)的發(fā)展方向。

    子系統(tǒng)的熱管理問題已被廣泛地研究并應(yīng)用于不同類型的飛行器中。Balland等[1]研究了高速運(yùn)輸機(jī)的熱管理循環(huán),基于燃油預(yù)熱和電能與機(jī)械能的相互轉(zhuǎn)化關(guān)系設(shè)計(jì)了熱/能量管理系統(tǒng),通過評(píng)估得到了在給定飛行任務(wù)下的可行設(shè)計(jì)方案。Doman[2-3]提出了一種典型噴氣式飛機(jī)的燃油熱管理系統(tǒng),該系統(tǒng)主要用于確定最佳的巡航條件,推導(dǎo)了以燃料作為熱沉?xí)r,飛機(jī)性能參數(shù)和大氣環(huán)境參數(shù)等因素對(duì)燃油溫度的影響。Haselbach等[4]描述了未來大型民用飛機(jī)的發(fā)動(dòng)機(jī)系統(tǒng)、集成方法和潛在的好處。Seki等[5]提出了一種基于空氣/燃油和蒸發(fā)循環(huán)相結(jié)合的綜合熱管理系統(tǒng),該研究的目的是有效地減少未來多電飛機(jī)的沖壓阻力和燃料消耗。Esser等[6]介紹了未來高超聲飛行器中4種熱管理方法(2種被動(dòng)方法和2種主動(dòng)方法),通過模擬仿真和測(cè)試,驗(yàn)證了這4種熱管理方法的可行性。Zilio等[7]通過對(duì)被動(dòng)冷卻技術(shù)和主動(dòng)冷卻技術(shù)的適當(dāng)結(jié)合,提出了一種應(yīng)用于直升機(jī)的航空電子設(shè)備熱管理方案,開發(fā)了原理樣機(jī),通過試驗(yàn)測(cè)試了不同熱載荷下的系統(tǒng)性能。NASA格倫研究中心[8]為商用飛機(jī)設(shè)計(jì)了一套熱管理系統(tǒng),在飛機(jī)內(nèi)部布置熱傳輸管路網(wǎng)絡(luò),通過聲波管提取發(fā)動(dòng)機(jī)廢熱,通過熱泵、熱管和燃油系統(tǒng)提取電子設(shè)備廢熱,將這些低品位熱轉(zhuǎn)化為高品位熱之后,通過管路網(wǎng)絡(luò)將熱量傳遞到全機(jī)身,經(jīng)回收利用之后,回到發(fā)動(dòng)機(jī)燃燒室,用以增加發(fā)動(dòng)機(jī)功率。真正的廢熱最后由排氣系統(tǒng)排出,既解決了整機(jī)的熱管理問題,又提高了發(fā)動(dòng)機(jī)功率,降低了燃油消耗率。歐洲TOICA項(xiàng)目[9]打算從根本上改進(jìn)飛機(jī)器的設(shè)計(jì)方法,研究?jī)?nèi)容包括發(fā)展協(xié)同仿真能力;提高數(shù)據(jù)的生成、管理和成熟度;提出新的熱管理方法;評(píng)估和驗(yàn)證已有的開發(fā)能力和技術(shù)通用性等方面。

    另一方面,隨著紅外探測(cè)技術(shù)的發(fā)展,直升機(jī)受到紅外武器打擊的威脅。目前,直升機(jī)主要采用排氣紅外抑制器和紅外抑制涂層2種手段來降低直升機(jī)的紅外輻射。YAH-64“阿帕奇”直升機(jī)上的“黑洞”紅外抑制器是一種低成本、沒有任何移動(dòng)部件的系統(tǒng)[10],它通過混合管將新鮮的冷空氣和發(fā)動(dòng)機(jī)熱排氣混合,大大降低了排氣系統(tǒng)的紅外輻射信號(hào)。A129“貓鼬”直升機(jī)也采用了類似的紅外抑制器[11]?!昂k唷敝鄙龣C(jī)上的波瓣紅外抑制器由于其優(yōu)越的引射能力和混合效率被人們廣泛地研究[12-17],探討了波瓣的渦旋結(jié)構(gòu)和混合機(jī)理。RAH-66“科曼奇”直升機(jī)是把紅外隱身、雷達(dá)隱身與氣動(dòng)布局三者統(tǒng)籌考慮的直升機(jī)[18-19],其紅外抑制器位于尾梁中,采用獨(dú)特的長(zhǎng)條形排氣口設(shè)計(jì),有足夠的長(zhǎng)度使發(fā)動(dòng)機(jī)的排氣和外界冷卻空氣充分摻混。在直升機(jī)表面涂覆紅外隱身涂層能夠有效降低其紅外輻射,從而增強(qiáng)直升機(jī)的隱蔽性,其具有工藝簡(jiǎn)單、不改變目標(biāo)外形等優(yōu)點(diǎn)。隨著技術(shù)的成熟和發(fā)展,紅外隱身涂層的研究主要集中在如何降低涂層發(fā)射率、如何實(shí)現(xiàn)涂層高溫穩(wěn)定性、涂層自適應(yīng)紅外隱身技術(shù)以及涂層多波段兼容隱身設(shè)計(jì)等方面。

    綜上所述,雖然人們?cè)谥鄙龣C(jī)子系統(tǒng)的熱管理(如電子設(shè)備冷卻系統(tǒng)設(shè)計(jì)[7]、環(huán)控系統(tǒng)設(shè)計(jì)分析[20]以及紅外抑制設(shè)計(jì)[17-18])和整機(jī)紅外輻射特性的研究都已經(jīng)取得了進(jìn)展。但不管是固定翼飛機(jī)還是直升機(jī),在這兩個(gè)方面的研究仍是獨(dú)立分析和獨(dú)立優(yōu)化,沒有考慮耦合的影響。當(dāng)處理系統(tǒng)熱問題時(shí),人們主要關(guān)注子系統(tǒng)的熱傳遞問題,即熱量的產(chǎn)生、提取、傳輸及耗散等過程。當(dāng)設(shè)計(jì)紅外抑制措施時(shí),主要為對(duì)紅外涂層和排氣紅外抑制器的設(shè)計(jì)研究。從直升機(jī)整機(jī)的角度出發(fā),其紅外輻射特性不僅受到紅外涂層和紅外抑制器設(shè)計(jì)的影響,還受到飛機(jī)內(nèi)部各系統(tǒng)部件的散熱量和機(jī)身表面的散熱口布置的影響,這些影響因素相互影響和制約,必須統(tǒng)籌考慮和分析。

    因此,為了應(yīng)對(duì)未來直升機(jī)面臨的熱管理和紅外抑制問題,需要從整機(jī)層面探究熱管理和紅外輻射特性耦合機(jī)制,對(duì)熱/質(zhì)實(shí)現(xiàn)優(yōu)化分配和動(dòng)態(tài)管理,同時(shí)發(fā)展熱量高效轉(zhuǎn)換與利用技術(shù),以實(shí)現(xiàn)能量利用效率的提高,并優(yōu)化整機(jī)紅外輻射性能的目的。本文從直升機(jī)整機(jī)系統(tǒng)的結(jié)構(gòu)特點(diǎn)、系統(tǒng)工作原理和能量平衡關(guān)系出發(fā),基于熱/質(zhì)流產(chǎn)生、收集傳輸、儲(chǔ)存利用和傳熱關(guān)系,建立了以系統(tǒng)部件為內(nèi)部邊界條件、外環(huán)境為外部邊界條件的整機(jī)耦合傳熱模型,通過數(shù)值仿真得到系統(tǒng)工作狀態(tài)、整機(jī)溫度和紅外輻射特性,并以此提出了3個(gè)熱優(yōu)化方向,為直升機(jī)整機(jī)熱管理和紅外抑制提供方法支撐。

    1 直升機(jī)能量平衡關(guān)系

    直升機(jī)整機(jī)溫度和紅外輻射特性分布受到內(nèi)部因素和外部因素的共同影響,能量平衡關(guān)系如圖1所示。在直升機(jī)外部,旋翼產(chǎn)生的下洗流場(chǎng)、前飛流場(chǎng)和飛行環(huán)境共同影響了整機(jī)周圍的對(duì)流換熱,環(huán)境冷空氣與發(fā)動(dòng)機(jī)熱排氣摻混,降低了熱排氣溫度。同時(shí),熱排氣受到外流場(chǎng)的影響,沖擊后機(jī)身、尾槳涵道和垂直尾翼使其表面溫度升高。另外,機(jī)身表面還受到太陽、天空和地面輻射以及環(huán)境溫度、濕度等外部環(huán)境因素的影響。而在直升機(jī)內(nèi)部,各子系統(tǒng)及其部件之間存在能量耦合關(guān)系,且部件作為內(nèi)熱源,其產(chǎn)生的熱量通過固體結(jié)構(gòu)的熱傳導(dǎo)、艙內(nèi)空氣的熱對(duì)流和壁面間熱輻射的形式影響直升機(jī)整機(jī)內(nèi)部的溫度分布。因此,在進(jìn)行直升機(jī)整機(jī)熱管理和紅外抑制設(shè)計(jì)研究時(shí),必須綜合考慮內(nèi)外影響因素的耦合作用,通過對(duì)整機(jī)熱量合理調(diào)控與管理,既提高整機(jī)能量利用率,又降低機(jī)身表面溫度,優(yōu)化紅外輻射性能。

    圖1 直升機(jī)能量平衡關(guān)系Fig.1 Energy balance of helicopter

    根據(jù)直升機(jī)系統(tǒng)特點(diǎn)、能量需求和散熱需求,要在整機(jī)層面建立直升機(jī)熱管理和紅外抑制的耦合分析方法仍需考慮以下困難和挑戰(zhàn):① 直升機(jī)各子系統(tǒng)相對(duì)獨(dú)立且尚未提出整機(jī)熱管理的概念;② 對(duì)整機(jī)子系統(tǒng)的集成熱管理需要重新考慮能量利用策略、部件互聯(lián)和性能的可靠性;③ 發(fā)動(dòng)機(jī)和電子設(shè)備功率的提高需要設(shè)計(jì)更高冷卻效率的冷卻系統(tǒng);④ 發(fā)動(dòng)機(jī)性能的提升使得其散熱量和排氣溫度升高,從而導(dǎo)致整機(jī)紅外輻射強(qiáng)度增大;⑤ 熱管理方法與紅外抑制設(shè)計(jì)相互影響且制約,因此需要統(tǒng)籌建模、分析和優(yōu)化;⑥ 在直升機(jī)整機(jī)層面的系統(tǒng)綜合一體化改造需建立完善的能效評(píng)估體系。

    2 整機(jī)熱模型

    本文所采用的物理模型如圖2所示,結(jié)構(gòu)包括駕駛艙、電子設(shè)備艙、油箱、發(fā)動(dòng)機(jī)艙和主減速器艙等結(jié)構(gòu)。2個(gè)發(fā)動(dòng)機(jī)位于機(jī)身兩側(cè)發(fā)動(dòng)機(jī)艙內(nèi);駕駛艙設(shè)有座椅、駕駛控制電子設(shè)備;電子設(shè)備艙設(shè)有電子設(shè)備箱;油箱位于駕駛艙下側(cè)。模型主要參數(shù)為:起飛總重量4 000 kg,機(jī)身長(zhǎng)11.6 m,機(jī)身寬2.0 m,機(jī)身高3.6 m。主旋翼采用4片槳葉,直徑為6.0 m,轉(zhuǎn)速為350 r/min。

    圖2 典型直升機(jī)結(jié)構(gòu)Fig.2 Typical helicopter structure

    本文所建立的直升機(jī)熱管理和紅外輻射特性耦合仿真和優(yōu)化方法如圖3所示,首先基于能量流動(dòng)和功能的耦合關(guān)系,建立部件級(jí)、子系統(tǒng)級(jí)和整機(jī)級(jí)的能量耦合計(jì)算模型,描述部件及子系統(tǒng)的共同工作關(guān)系和部件散熱情況;然后建立以系統(tǒng)部件為內(nèi)部邊界條件、外環(huán)境為外部邊界條件的直升機(jī)整機(jī)耦合傳熱計(jì)算模型,探究直升機(jī)的流場(chǎng)、溫度場(chǎng)和紅外輻射特性分布;最后根據(jù)仿真結(jié)果提出優(yōu)化方向,如系統(tǒng)布局和能量利用策略、紅外抑制器、隔熱和低紅外發(fā)射率涂層等多方面的改進(jìn)優(yōu)化方法。通過這種方法耦合建模、分析和優(yōu)化,實(shí)現(xiàn)整機(jī)的熱量統(tǒng)籌管理。

    2.1 子系統(tǒng)能量耦合計(jì)算模型

    直升機(jī)艙內(nèi)各個(gè)部件的散熱情況影響了直升機(jī)內(nèi)部空間的溫度分布。要計(jì)算各部件的散熱情況,首先要建立各部件的能量耦合關(guān)系。根據(jù)直升機(jī)內(nèi)部各部件結(jié)構(gòu)特點(diǎn)、工作原理和部件間的能量平衡關(guān)系,可建立如圖4所示的直升機(jī)子系統(tǒng)能量耦合關(guān)系示意圖[21],包括旋翼子系統(tǒng)、發(fā)動(dòng)機(jī)子系統(tǒng)、燃料子系統(tǒng)、滑油子系統(tǒng)和環(huán)境控制子系統(tǒng)。其中,電子冷卻子系統(tǒng)和液壓子系統(tǒng)簡(jiǎn)化為燃油子系統(tǒng)中的熱載荷。

    圖3 熱管理和紅外輻射特性耦合仿真和優(yōu)化方法Fig.3 Coupling simulation and optimization method of thermal management and infrared radiation characteristics

    圖4 直升機(jī)子系統(tǒng)能量耦合關(guān)系[21]Fig.4 Energy coupling relationship of helicopter subsystems[21]

    發(fā)動(dòng)機(jī)子系統(tǒng)是直升機(jī)能量最大來源,軸功由自由渦輪輸出經(jīng)傳動(dòng)系統(tǒng)傳向主旋翼子系統(tǒng)、尾槳等各類附件。在穩(wěn)態(tài)工況下,各部件滿足質(zhì)量平衡、能量平衡和功率平衡關(guān)系,發(fā)動(dòng)機(jī)的共同工作特性可根據(jù)以下方程迭代計(jì)算[21]:

    (1)

    在動(dòng)態(tài)工況下,上述穩(wěn)態(tài)共同工作特性不再滿足,此時(shí)發(fā)動(dòng)機(jī)的動(dòng)態(tài)共同工作特性滿足容積動(dòng)力學(xué)和轉(zhuǎn)子動(dòng)力學(xué)方程[22]:

    (2)

    旋翼子系統(tǒng)的需求功率Pr根據(jù)飛行高度、馬赫數(shù)和飛行環(huán)境確定,則發(fā)動(dòng)機(jī)軸輸出功Pe與旋翼需求功率Pr的關(guān)系為[23]

    Pr=ζ·Pe

    (3)

    式中:ζ為機(jī)械傳遞效率。

    (4)

    式中:θ為空氣利用系數(shù);Qcabin為駕駛艙或者電子設(shè)備艙的熱載荷;cba為引氣比熱容;Tcabin為駕駛艙或者電子設(shè)備艙的控制溫度;Tin為進(jìn)口溫度。

    (5)

    子系統(tǒng)及其各部件的詳細(xì)熱力學(xué)數(shù)學(xué)模型見文獻(xiàn)[21]。通過對(duì)子系統(tǒng)能量耦合性能的仿真可獲得直升機(jī)各部件的工作參數(shù)和散熱情況。

    2.2 流場(chǎng)和溫度場(chǎng)耦合傳熱模型

    直升機(jī)整機(jī)流場(chǎng)、溫度場(chǎng)和紅外輻射特性受到飛行環(huán)境、飛行馬赫數(shù)、旋翼流場(chǎng)、排氣流場(chǎng)、太陽輻射和機(jī)身內(nèi)部部件的散熱情況等因素的綜合影響。因此,流場(chǎng)和溫度場(chǎng)計(jì)算涉及到的計(jì)算區(qū)域包括直升機(jī)內(nèi)部和外部空間,控制方程包括連續(xù)性方程、動(dòng)量守恒方程、能量守恒方程、組分輸運(yùn)方程以及輻射傳輸方程。連續(xù)性方程應(yīng)用于整個(gè)計(jì)算流體中,適用于可壓和不可壓流體:

    (6)

    式中:ρ為密度;v為速度矢量;Sm為質(zhì)量源項(xiàng)。慣性坐標(biāo)系下的動(dòng)量方程為

    (7)

    式中:p為靜壓;τ為應(yīng)力張量;ρg是重力產(chǎn)生的體積力;F為附加體積力或其他自定義源項(xiàng)。計(jì)算域中的能量守恒方程為

    (8)

    式中:E為總能;keff為流體的有效導(dǎo)熱系數(shù);T為溫度;hj和Jj分別為組分j的焓和擴(kuò)散通量;τeff為有效應(yīng)力張量;Sh為能量源項(xiàng)。輻射熱流是作為能量方程中的一個(gè)能量源項(xiàng)考慮,因此,在計(jì)算區(qū)域內(nèi)需要考慮輻射傳輸方程,即

    (9)

    式中:I(r,s)為輻射強(qiáng)度;r為位置矢量;s為方向矢量;s為沿程長(zhǎng)度;α為吸收系數(shù);σs為散射系數(shù);n為折射系數(shù);σ為Stefan-Boltzmann常數(shù);s′為散射方向向量;Φ為散射相函數(shù);Ω′為立體角。

    發(fā)動(dòng)機(jī)高溫排氣中含有的二氧化碳、水蒸氣等氣體也會(huì)參與輻射傳輸,且與環(huán)境冷空氣進(jìn)行摻混,所以在整個(gè)計(jì)算域內(nèi)需要考慮氣體組分輸運(yùn)方程:

    (10)

    式中:Yj為組分的當(dāng)?shù)刭|(zhì)量濃度。

    固體壁面的對(duì)流換熱系數(shù)可計(jì)算如下:

    (11)

    通過運(yùn)用數(shù)值方法求解式(6)~式(10),即可獲得直升機(jī)整機(jī)溫度分布。求解思路又可分為整體耦合求解法和區(qū)域分解求解法。整體耦合求解法是對(duì)直升機(jī)所有求解區(qū)域進(jìn)行流場(chǎng)和溫度場(chǎng)的耦合求解,流固界面轉(zhuǎn)化為求解區(qū)域內(nèi)部,不再需要設(shè)置界面的邊界條件,簡(jiǎn)化了數(shù)值計(jì)算問題的處理,這里不再贅述;區(qū)域分解求解法是將直升機(jī)所有區(qū)域劃分為幾個(gè)子區(qū)域,分別建立各子區(qū)域的數(shù)值計(jì)算模型,再考慮各區(qū)域之間的相互影響,必須給出不同區(qū)域界面處的邊界條件,再建立整體的數(shù)值計(jì)算模型,如下所述:

    2.2.1 直升機(jī)內(nèi)部傳熱數(shù)值計(jì)算模型

    對(duì)于區(qū)域分解的數(shù)值計(jì)算方法,在直升機(jī)艙內(nèi)部流體域中,可運(yùn)用集總參數(shù)法描述艙內(nèi)空氣溫度的變化,則能量平衡方程式(8)可寫為[25]

    (12)

    當(dāng)環(huán)控系統(tǒng)關(guān)閉時(shí),艙內(nèi)流動(dòng)狀況為大空間內(nèi)的自然對(duì)流條件下的換熱,則艙內(nèi)表面的換熱系數(shù)hc可以計(jì)算如下[26]:

    (13)

    式中:Nuc為艙內(nèi)努塞爾數(shù);Gr為格拉曉夫數(shù);Pr為氣體普朗特?cái)?shù);kc為艙內(nèi)氣體導(dǎo)熱系數(shù);l為特征尺寸;系數(shù)C和指數(shù)b與物體幾何形狀及GrPr數(shù)值的大小等因素有關(guān)。艙內(nèi)定性溫度Tm采用邊界層的算數(shù)平均數(shù)計(jì)算,即

    (14)

    式中:Tw為艙內(nèi)壁面表面的溫度;T∞為未受壁面影響流體的溫度。

    當(dāng)環(huán)控系統(tǒng)開啟時(shí),艙內(nèi)流動(dòng)狀況為強(qiáng)迫流動(dòng)條件下的對(duì)流換熱,則艙內(nèi)表面的換熱系數(shù)hc根據(jù)以下公式計(jì)算[24]:

    hc=11.34+5.84vc

    (15)

    式中:vc為艙內(nèi)氣流速度。

    艙內(nèi)熱載荷包括各種部件的散熱和乘員生理散熱等,其中,發(fā)動(dòng)機(jī)是艙內(nèi)的最大熱源,對(duì)整機(jī)溫度分布有很大影響。根據(jù)子系統(tǒng)能量耦合模型可獲得發(fā)動(dòng)機(jī)主流氣體沿發(fā)動(dòng)機(jī)軸向的溫度分布Tg(x),則沿發(fā)動(dòng)機(jī)軸向的機(jī)匣表面溫度分布根據(jù)以下公式計(jì)算:

    (16)

    式中:a為熱擴(kuò)散系數(shù);r為機(jī)匣半徑;r1和r2分別為機(jī)匣內(nèi)、外半徑;h1和h2分為機(jī)匣內(nèi)外側(cè)的對(duì)流換熱系數(shù);kw為機(jī)匣材料導(dǎo)熱系數(shù);Tw1和Tw2分別為機(jī)匣內(nèi)壁面和外壁面的溫度;qg,rad為高溫氣體對(duì)機(jī)匣內(nèi)側(cè)的輻射換熱量;qc,rad為機(jī)匣外側(cè)與艙內(nèi)其他壁面的輻射換熱量。直升機(jī)艙內(nèi)第i個(gè)表面的凈輻射熱流qi可以表示為

    qi=Ji-Gi=εiEbi-αiGi

    (17)

    式中:Ji和Gi分別為第i個(gè)表面的有效輻射和入射輻射;εi為表面發(fā)射率;Ebi為第i個(gè)表面的黑體輻射;αi為表面吸收率。有效輻射Ji可以表示為[27]

    (18)

    式中:Ti為表面溫度;ρi為表面反射率。

    假設(shè)艙內(nèi)無二氧化碳和水蒸氣等參與性介質(zhì),則入射輻射可以用以下公式計(jì)算:

    (19)

    式中:Fi-j為面元i對(duì)面元j的輻射角系數(shù),可通過蒙特卡洛法來計(jì)算。

    2.2.2 蒙皮簡(jiǎn)化處理

    直升機(jī)蒙皮將內(nèi)外區(qū)域的傳熱關(guān)系耦合起來,蒙皮主要采用復(fù)合蜂窩材料,結(jié)構(gòu)和傳熱原理如圖5所示,當(dāng)蜂窩板受到熱載荷作用時(shí),其內(nèi)部存在3種傳熱方式,即蜂窩芯固體和氣體的導(dǎo)熱換熱、蜂窩芯內(nèi)表面間的輻射換熱以及蜂窩芯內(nèi)氣體的對(duì)流換熱。

    為了簡(jiǎn)化計(jì)算,本文將蜂窩蒙皮簡(jiǎn)化為具有一定厚度的固體壁面。因此,在蒙皮中,連續(xù)性方程和動(dòng)量守恒方程式(8)不再求解,且認(rèn)為熱量只在厚度方向上傳遞,則能量方程退化為熱傳導(dǎo)方程,即

    (20)

    式中:ρequ、cequ和kequ分別為蒙皮等效密度、等效比熱容和等效導(dǎo)熱系數(shù);hc和ho分別為蒙皮內(nèi)側(cè)和外側(cè)的對(duì)流換熱系數(shù);Tc,w和To,w分別為艙蒙皮內(nèi)側(cè)和外側(cè)的壁面溫度;Te為蒙皮外側(cè)附面層溫度;qc,rad為蒙皮內(nèi)側(cè)壁面與艙內(nèi)其他壁面之間的輻射換熱量,可根據(jù)式(17)~式(19)計(jì)算;qo,rad為蒙皮外側(cè)壁面與環(huán)境、機(jī)身其他壁面和熱排氣之間的輻射換熱量。

    1-下面板;2-蜂窩芯;3-上面板圖5 蒙皮結(jié)構(gòu)示意圖Fig.5 Schematic diagram of skin structure

    蒙皮等效密度、等效比熱容和等效導(dǎo)熱系數(shù)可分別根據(jù)以下公式計(jì)算[28-29]:

    (21)

    (22)

    (23)

    式中:ρs和ρg分別為蜂窩芯內(nèi)部的固體和氣體的密度;As和Ag分別為蜂窩芯橫截面中固體和氣體所占的面積;cs和cg分別為固體和氣體的比熱;ks和kg分別為固體和氣體的導(dǎo)熱系數(shù);krad為蜂窩芯內(nèi)部輻射導(dǎo)熱系數(shù),可根據(jù)以下公式計(jì)算[28-29]:

    (24)

    式中:ξ為蜂窩芯的高徑比;ε為蜂窩單元內(nèi)表面發(fā)射率;L為蜂窩芯高度;T1、T2分別為蜂窩板上下表面的溫度。

    蒙皮外側(cè)壁面輻射換熱中的環(huán)境入射輻射包括太陽輻射、地面輻射、天空背景輻射等[25,30],則qo,rad可根據(jù)以下方程計(jì)算:

    qo,rad=εiEbi-αi(Gg+Gs)i-

    (25)

    式中:Gs和Gg分別為機(jī)身其他壁面和熱排氣的入射輻射;αi,s為i表面的太陽吸收率;A′i為考慮遮擋后i表面被太陽照射到的部分在垂直于太陽入射光線的平面上的投影面積;Ai為i表面的實(shí)際換熱面積;qD為太陽直射輻射;qd為太陽散射輻射;qR為地面對(duì)太陽的反射輻射;qsky為天空輻射;qgr為地面輻射。太陽直射輻射根據(jù)以下公式計(jì)算:

    qD=rsIscPma

    (26)

    式中:rs為日地間距引起的修正值;Isc為太陽常數(shù);P為大氣透明度;ma為大氣質(zhì)量。水平面和傾斜面的太陽散射輻射qd分別用以下公式計(jì)算:

    (27)

    (28)

    式中:C1和C2為經(jīng)驗(yàn)系數(shù),其值取決于大氣透明情況;hs為太陽高度角;βs為斜面傾角。直升機(jī)接收地面對(duì)太陽的反射輻射qR根據(jù)以下公式計(jì)算:

    (29)

    式中:ρgr為地面反射率。

    對(duì)于晴朗天空,傾斜表面上的天空輻射qsky和地面輻射qgr可分別用以下公式計(jì)算:

    (30)

    (31)

    式中:εsky為大氣等效發(fā)射率;Ta為空氣溫度;εgr為地表表面發(fā)射率;Tgr為地面溫度。

    艙外蒙皮還會(huì)受到其他壁面和高溫排氣的熱輻射,假設(shè)熱排氣無散射,則輻射傳輸方程式(9)可以簡(jiǎn)化為

    (32)

    將該方程在傳輸方向上進(jìn)行離散,即將光線傳遞的行程l分為n層,設(shè)第m層的厚度為Δlm,則其他壁面和高溫排氣的入射輻射Gg+Gs可以用以下方程計(jì)算[31]:

    Ebg(T1)[1-τ(1)]+

    (33)

    式中:Ebg為氣體黑體輻射;Fi-j為面元i對(duì)面元j的輻射角系數(shù);Tm為第m層氣體溫度;τ(m)為第m層氣體的透射率,可用以下公式計(jì)算[31]:

    τ(m)=exp[-Δlmα(m)]

    (34)

    式中:α(m)表示第m層氣體的吸收系數(shù),與氣體的溫度、壓力和組分濃度等參數(shù)有關(guān)。若光線從機(jī)身固體壁面射出時(shí),E(0)為對(duì)應(yīng)固體壁面的有效輻射;當(dāng)光線從天空或地面射出時(shí),E(0)=0。

    2.2.3 直升機(jī)外部傳熱數(shù)值計(jì)算模型

    在直升機(jī)外部,飛行環(huán)境、旋翼氣流、排氣氣流和機(jī)身形狀等因素都會(huì)對(duì)直升機(jī)外部的氣流速度、溫度和排氣組分分布產(chǎn)生影響,從而影響直升機(jī)的紅外輻射特性。因此,需要耦合求解方程式(6)~式(10)來計(jì)算艙外空氣流場(chǎng)、溫度場(chǎng)和氣體組分分布。其中,飛行環(huán)境溫度Th和壓力Ph根據(jù)以下公式確定:

    Ph=

    (35)

    式中:H為飛行高度。

    前飛速度V0根據(jù)以下公式計(jì)算:

    (36)

    式中:Ma為飛行馬赫數(shù)。

    旋翼槳盤處的誘導(dǎo)速度分布v1(r)可根據(jù)旋翼子系統(tǒng)模型計(jì)算,則旋翼槳盤處的速度V1分布為誘導(dǎo)速度v1(r)和前飛速度V0的矢量和[23]:

    V1=V0+v1(r)

    (37)

    發(fā)動(dòng)機(jī)排氣口的溫度、質(zhì)量流量和壓力根據(jù)發(fā)動(dòng)機(jī)模型計(jì)算。排氣口的氣體組分根據(jù)發(fā)動(dòng)機(jī)空氣和燃油的流量,并按照化學(xué)平衡方程式確定(燃油化學(xué)式以C12H23為例)[32]:

    C12H23(l)+17.75·αaO2+17.75×3.76·αaN2→

    12CO2+11.5H2O+17.75×3.76·αaN2+

    17.75·(αa-1)O2

    (38)

    式中:αa為余氣系數(shù)(過量空氣系數(shù))。

    艙外空氣的流動(dòng)一般按照平板強(qiáng)迫對(duì)流的情況處理,則蒙皮外表面的換熱系數(shù)ho可根據(jù)以下方程確定[24]:

    (39)

    式中:Re為氣體雷諾數(shù);ko為艙外氣體導(dǎo)熱系數(shù)。艙外定性溫度計(jì)算中,艙外溫度一般取為附面層的溫度Te,即

    (40)

    式中:re為恢復(fù)系數(shù)。

    2.3 紅外輻射通量計(jì)算模型

    對(duì)于直升機(jī)任一單元表面上的總輻射通量Etotal包括自身輻射通量Eλ1-λ2和反射輻射通量ER兩部分:

    Etotal=Eλ1-λ2+ER

    (41)

    其中,自身紅外輻射通量可以由普朗特公式[26]出發(fā),通過對(duì)其在紅外波段范圍的積分得到,即

    (42)

    式中:λ1和λ2分別為給定的紅外波段范圍的上下限;ελ為單元表面的光譜發(fā)射率;B1為第一輻射常數(shù);B2為第二輻射常數(shù)。反射輻射通量可以用以下公式計(jì)算:

    (43)

    式中:ρλ為單元表面對(duì)第j個(gè)背景的光譜反射率;Hλ為第j個(gè)背景(如太陽、天空、地面等)的光譜輻射能量,分別可根據(jù)式(26)~式(31)計(jì)算。

    3 計(jì)算結(jié)果與討論

    直升機(jī)典型飛行剖面如圖6所示,該飛行剖面包含了水平加速、勻速上升、水平勻速、水平減速、勻速下降等典型飛行過程,總歷時(shí)1 200 s(20 min)。圖7所示為部分部件能量耦合模型仿真計(jì)算結(jié)果。從圖中可見,燃燒室出口溫度和排氣溫度的變化主要受馬赫數(shù)影響,且由于發(fā)動(dòng)機(jī)容積慣性和轉(zhuǎn)子慣性的影響,溫度變化存在一定的響應(yīng)時(shí)間?;偷幕赜蜏囟群腿加偷倪M(jìn)油溫度隨飛行的時(shí)間變化逐漸上升。飛行任務(wù)結(jié)束時(shí),滑油的回油溫度上升至382 K,燃油進(jìn)油溫度上升至358 K。燃油和滑油的溫度變化主要受到飛行狀態(tài)、燃油和滑油循環(huán)中各個(gè)熱源部件的發(fā)熱量以及散熱器布局的影響。

    圖6 典型直升機(jī)飛行剖面Fig.6 Typical flight profile of helicopter

    圖7 能量耦合模型計(jì)算結(jié)果Fig.7 Computational results of energy coupling model

    直升機(jī)流場(chǎng)和溫度場(chǎng)計(jì)算的網(wǎng)格如圖8所示。計(jì)算域包括直升機(jī)內(nèi)部和外部2個(gè)計(jì)算域,面網(wǎng)格采用三角形網(wǎng)格,體網(wǎng)格采用非結(jié)構(gòu)化四面體網(wǎng)格。在計(jì)算得到直升機(jī)各個(gè)子系統(tǒng)的工作參數(shù)和散熱量之后,將旋翼下洗速度分布、發(fā)動(dòng)機(jī)機(jī)匣溫度分布、外環(huán)境溫度和壓力、機(jī)身內(nèi)部熱源散熱量和太陽輻射強(qiáng)度以邊界條件的形式賦到流場(chǎng)和溫度場(chǎng)耦合計(jì)算模型中。圖9所示為網(wǎng)格無關(guān)性驗(yàn)證,共生成了網(wǎng)格數(shù)量為222萬和445萬兩套網(wǎng)格,網(wǎng)格質(zhì)量均大于0.2,仿真時(shí)間間隔采用0.5 s和1 s,監(jiān)控點(diǎn)為發(fā)動(dòng)機(jī)排氣口中心處的氣體溫度。從計(jì)算結(jié)果可知,采用222萬網(wǎng)格、0.5 s時(shí)間間隔的算例和445萬網(wǎng)格、0.5 s時(shí)間間隔的算例的計(jì)算結(jié)果已基本吻合,因此為了節(jié)約計(jì)算資源,本文采用222萬網(wǎng)格、0.5 s時(shí)間間隔的算例。

    仿真計(jì)算工況如圖6所示,圖10所示為直升機(jī)在120 s時(shí)刻的整機(jī)溫度分布,此時(shí)排氣溫度約為895 K。從圖中可以看出,在前飛狀態(tài)中,整機(jī)溫度最高的區(qū)域?yàn)榕艢夤苤車鷧^(qū)域。同時(shí)發(fā)動(dòng)機(jī)整流罩受到發(fā)動(dòng)機(jī)散熱和熱排氣的影響,溫度較高。同時(shí),發(fā)動(dòng)機(jī)熱排氣在前飛流場(chǎng)和旋翼流場(chǎng)的共同作用下沖擊垂直尾翼,使得垂直尾翼溫度升高。

    圖8 計(jì)算域和網(wǎng)格細(xì)節(jié)Fig.8 Computational domains and mesh details

    圖9 網(wǎng)格無關(guān)性驗(yàn)證Fig.9 Mesh independence verification

    直升機(jī)3~5 μm波段和8~14 μm波段的固體壁面紅外輻射通量分布如圖11和圖12所示。從圖中可知,在3~5 μm波段中,紅外輻射通量最大的區(qū)域?yàn)榘l(fā)動(dòng)機(jī)系統(tǒng)影響區(qū)域,可以清楚地看到排氣管和發(fā)動(dòng)機(jī)機(jī)匣壁面的紅外輻射通量。在8~14 μm波段,發(fā)動(dòng)機(jī)系統(tǒng)影響區(qū)域的機(jī)身蒙皮和垂直尾翼的紅外輻射通量明顯增大。

    圖10 直升機(jī)整機(jī)溫度分布Fig.10 Temperature distribution of helicopter

    圖12 8~14 μm波段紅外輻射通量Fig.12 Infrared radiation flux in 8-14 μm band

    4 熱管理優(yōu)化分析

    基于以上計(jì)算結(jié)果,本文提出了在熱管理和紅外抑制的優(yōu)化過程中以空間層面、時(shí)間層面和熱回收利用為出發(fā)點(diǎn)的3個(gè)優(yōu)化方向。在空間層面,分析了部件層面的換熱器設(shè)計(jì)和子系統(tǒng)層面系統(tǒng)布局的優(yōu)化設(shè)計(jì)方法。在時(shí)間層面,分析了影響發(fā)動(dòng)機(jī)和換熱器時(shí)間常數(shù)的因素。在熱回收利用方面,提出了熱回收利用的熱管理方案,分析該改進(jìn)方案對(duì)發(fā)動(dòng)機(jī)性能和紅外輻射特性的優(yōu)化效果。

    4.1 空間層面

    直升機(jī)各個(gè)系統(tǒng)的空間結(jié)構(gòu)呈現(xiàn)3個(gè)層次結(jié)構(gòu),即整機(jī)層、子系統(tǒng)層和部件層。熱源也呈現(xiàn)空間分布的不均勻性和隨時(shí)間變化的特點(diǎn)。因此,要實(shí)現(xiàn)整機(jī)層面的熱利用最優(yōu),就需要在這3個(gè)層次上逐步優(yōu)化。下文以部件層的換熱器和整機(jī)層的系統(tǒng)布局為例分析空間層面的熱優(yōu)化。

    4.1.1 換熱器優(yōu)化

    隨著直升機(jī)性能需求的不斷提高,整機(jī)不同部位的熱載荷增加,這就需要設(shè)計(jì)不同形式的換熱器以適應(yīng)不同的散熱需求,如高效電子設(shè)備換熱器和機(jī)身蒙皮換熱器等。

    傳統(tǒng)電子設(shè)備的冷卻方法是空氣主動(dòng)冷卻,該方法簡(jiǎn)單可靠,但已快到冷卻極限,不能滿足高熱流密度的換熱要求。液體冷卻的換熱系數(shù)及比熱容比空氣大得多,具有更高的冷卻效率,同時(shí)采用泵驅(qū)動(dòng)式的閉式循環(huán)系統(tǒng)的工作能力也更穩(wěn)定可靠。圖13所示為電子設(shè)備液體冷卻系統(tǒng)及其冷板結(jié)構(gòu)的安裝示意圖。

    由于圓角S型流道具有壓力損失小、冷卻效果好、流量分配均勻等優(yōu)點(diǎn),本文設(shè)計(jì)了如圖14所示的4種改型S型流道結(jié)構(gòu)。

    計(jì)算工況和邊界條件為:電子模塊元器件的尺寸和功率見表1所示,冷板的厚度為5 mm,進(jìn)口冷卻液的流速為0.5 m/s,溫度為300 K,系統(tǒng)的環(huán)境設(shè)為空氣,空氣自然對(duì)流的換熱系數(shù)設(shè)為8 W/(m2·K),環(huán)境溫度設(shè)為300 K。

    圖13 電子設(shè)備液體冷卻系統(tǒng)Fig.13 Liquid cooling system for electronic equipment

    圖14 4種S型流道冷板結(jié)構(gòu)Fig.14 Four S-type cold plate structures

    表1 電子元器件尺寸和功率Table 1 Size and power of electronic components

    仿真結(jié)果如圖15和表2所示,在入口流量相同的情況下,流道數(shù)量增加導(dǎo)致單通道內(nèi)流速降低,但流道數(shù)目的增加使得換熱更加充分,所以結(jié)構(gòu)a、b、c的冷卻效果相差不大。隨著流道數(shù)目的增加,沿程的壓力損失降低,在加工工藝和強(qiáng)度允許的前提下,可適當(dāng)增加流道數(shù)目,如結(jié)構(gòu)c。沿著流動(dòng)方向,液體溫度逐漸上升,上游芯片溫度明顯低于下游芯片溫度,均溫性不夠好,故可考慮在上下冷卻液接口均設(shè)置流體進(jìn)出口。在流道中增加強(qiáng)化換熱翅片的設(shè)計(jì)結(jié)構(gòu)d的壓力損失更小,冷卻效果更優(yōu)。

    圖15 溫度分布Fig.15 Temperature distribution

    表2 冷板溫度計(jì)算結(jié)果Table 2 Computational results of cold plate temperature

    直升機(jī)飛行速度低,艙外空氣可作為熱沉使用,如圖16所示為一種蒙皮冷卻液-空氣換熱器,該換熱器將冷卻通道安裝在直升機(jī)蒙皮內(nèi)部,管道中流體的熱量通過管壁釋放到蒙皮上再擴(kuò)散到外部冷卻氣流中。這種散熱器可布置的面積大,不僅將直升機(jī)內(nèi)部的熱量擴(kuò)散到外部環(huán)境中,同時(shí)換熱器翅片產(chǎn)生旋渦使直升機(jī)后體邊界層的內(nèi)外氣流摻混,增加了邊界層內(nèi)氣流的能量,增強(qiáng)了抵抗逆壓梯度、延緩分離的能力,因此該換熱器設(shè)計(jì)還達(dá)到了減阻的目的。

    圖16 直升機(jī)蒙皮換熱器Fig.16 Helicopter skin heat exchanger

    4.1.2 系統(tǒng)布局優(yōu)化

    典型直升機(jī)子系統(tǒng)布局如圖17所示,各個(gè)系統(tǒng)分散,只有在燃滑油換熱器處有熱量交換,燃油系統(tǒng)無回油。各個(gè)子系統(tǒng)的熱量不易統(tǒng)一管理,且燃油流量小,不能作為冷卻流體。因此,為了將各子系統(tǒng)集成和對(duì)熱量的統(tǒng)一管理,采用如圖18所示的子系統(tǒng)設(shè)計(jì),以燃油和沖壓空氣為熱沉。將電子設(shè)備冷卻子系統(tǒng)、滑油子系統(tǒng)、液壓子系統(tǒng)、環(huán)控子系統(tǒng)和燃油子系統(tǒng)集成起來。

    圖17 典型直升機(jī)子系統(tǒng)布局Fig.17 Typical layout of helicopter subsystem

    圖18 直升機(jī)子系統(tǒng)集成熱管理Fig.18 Integrated thermal management of helicopter subsystems

    針對(duì)子系統(tǒng)的集成回路設(shè)計(jì),本文設(shè)計(jì)了以下4種回路布局方案:

    1) 簡(jiǎn)單串聯(lián)帶回路的布局方案

    這種布局方案增加了進(jìn)油路上的燃油流量,依次帶走電子設(shè)備冷卻回路、環(huán)控系統(tǒng)、滑油回路和液壓回路的熱量,如圖19所示,一部分燃油進(jìn)入發(fā)動(dòng)機(jī),其余的燃油經(jīng)回油路返回油箱。

    2) 換熱器并聯(lián)的布局方案

    對(duì)于換熱器數(shù)目較多的系統(tǒng),使用串聯(lián)式布局會(huì)大大增加沿程阻力,為了使系統(tǒng)維持在額定流量下需要選擇較大功率的泵,這將增加系統(tǒng)的熱載荷,增加燃油溫度,降低換熱器的換熱性能。因此,將滑油回路換熱器和液壓回路換熱器并聯(lián)設(shè)置,如圖20所示,2個(gè)支路的燃油流量根據(jù)熱載荷的比例進(jìn)行分配。

    3) 分主輔換熱器的布局方案

    直升機(jī)飛行速度低,可利用沖壓空氣來輔助冷卻。同時(shí),冷卻回油燃油后的沖壓空氣溫度不高,可在冷卻回油燃油之后對(duì)發(fā)動(dòng)機(jī)引氣作初步冷卻。在滑油回路和液壓油回路設(shè)置主輔換熱器,在主換熱器換熱不足的情況下,可開啟輔助換熱器做進(jìn)一步換熱,如圖21所示。

    4) 輔助換熱器冷側(cè)為沖壓空氣的布局方案

    當(dāng)燃油循環(huán)流量較小或輔助換熱器換熱效率較低的情況下,為了提高換熱效率,可將輔助換熱器設(shè)計(jì)為以沖壓空氣為冷側(cè)流體的布局方案,如圖22所示。

    為了比較上述不同回路布局方案的換熱性能,取一般巡航狀態(tài),飛行高度1 km,馬赫數(shù)0.25作為系統(tǒng)設(shè)計(jì)點(diǎn),各回路的熱載荷和循環(huán)流量如表3和表4所示。

    圖19 換熱器串聯(lián)的布局方案1Fig.19 Series layout 1 of heat exchangers

    圖21 分主輔換熱器的布局方案3Fig.21 Layout 3 with main and auxiliary heat exchangers

    圖22 輔助換熱器冷側(cè)為沖壓空氣的布局方案4Fig.22 Layout 4 of auxiliary heat exchangers with ram air on cold side

    不同布局穩(wěn)態(tài)仿真結(jié)果如表5所示,沿進(jìn)油方向,燃油溫度逐漸上升,使得布置在流路下游的換熱器換熱效率降低。布局2方案的整體溫升和布局1一致。因此,當(dāng)在進(jìn)油路上換熱部件較少時(shí)可采用布局1的設(shè)計(jì)。但當(dāng)換熱部件較多時(shí),為避免壓頭損失,采用并聯(lián)式的布局2較合適。采用布局3的設(shè)計(jì)降低了發(fā)動(dòng)機(jī)引氣進(jìn)入環(huán)控次級(jí)換熱器的溫度,使該換熱器燃油出口溫度降低約4.7 K,提高了滑油換熱器和液壓油換熱器的換熱效率?;秃鸵簤河统隹跍囟确謩e降低了11.8 K和10.9 K,換熱量增加,提高了燃油進(jìn)油溫度。當(dāng)輔助換熱器換熱效率較低時(shí),可通過沖壓空氣輔助冷卻,如布局4。

    表3 各回路熱載荷Table 3 Thermal load of each circuit

    表4 各回路循環(huán)油量Table 4 Mass flow rate of each circuit

    表5 監(jiān)控點(diǎn)計(jì)算溫度Table 5 Computational temperature of monitoring points

    4.2 時(shí)間層面

    空間層面的穩(wěn)態(tài)仿真并不能反應(yīng)實(shí)際飛行狀態(tài)變化過程中的時(shí)間響應(yīng)。直升機(jī)飛行工況的改變是通過控制燃油進(jìn)油油量以改變發(fā)動(dòng)機(jī)輸出功率來實(shí)現(xiàn)的,而發(fā)動(dòng)機(jī)的工況變化導(dǎo)致引氣參數(shù)和發(fā)動(dòng)機(jī)散熱量等各種參數(shù)的變化,這些參數(shù)都會(huì)影響環(huán)控子系統(tǒng)、燃油子系統(tǒng)等其他系統(tǒng)中部件的工作狀態(tài)和散熱量。飛行工況變化也影響了換熱器中沖壓空氣的參數(shù)。因此,各系統(tǒng)間的耦合導(dǎo)致整個(gè)系統(tǒng)的時(shí)間響應(yīng)問題非常復(fù)雜。整個(gè)系統(tǒng)中存在以下3種類型的時(shí)間響應(yīng)[33]:

    1) 控制體內(nèi)與流體的可壓縮性相關(guān)的慣性

    (44)

    2)熱部件的熱慣性

    (45)

    式中:cv和cp分別為控制體內(nèi)氣體定容比熱容和定壓比熱容。

    3)機(jī)械轉(zhuǎn)動(dòng)部件的轉(zhuǎn)動(dòng)慣性

    (46)

    式中:J為轉(zhuǎn)動(dòng)慣量;ω為轉(zhuǎn)動(dòng)角速度;Psurplus為剩余功率。

    對(duì)于發(fā)動(dòng)機(jī)來說,各部件容積小,容積慣性和熱慣性相對(duì)于轉(zhuǎn)動(dòng)慣性而言,其趨于平衡的速度更快,因此,轉(zhuǎn)軸是發(fā)動(dòng)機(jī)動(dòng)態(tài)模型中的最大慣性環(huán)節(jié)。圖23所示為直升機(jī)加速過程仿真結(jié)果。仿真時(shí)間為200 s,飛行馬赫數(shù)在100 s時(shí)從0.2加速到0.25。從仿真結(jié)果可見,轉(zhuǎn)動(dòng)慣量越大,發(fā)動(dòng)機(jī)工作參數(shù)達(dá)到穩(wěn)定的時(shí)間越久。當(dāng)轉(zhuǎn)動(dòng)慣量為J時(shí),工作參數(shù)再次達(dá)到穩(wěn)定所需的時(shí)間約為25 s;當(dāng)轉(zhuǎn)動(dòng)慣量為2J時(shí),工作參數(shù)再次達(dá)到穩(wěn)定所需的時(shí)間約為50 s。

    另外,系統(tǒng)中有大量的換熱器設(shè)備,對(duì)系統(tǒng)的動(dòng)態(tài)換熱特性有較大的影響。換熱器時(shí)間常數(shù)的影響因素很多,包括進(jìn)口流量、進(jìn)口溫度、傳熱系數(shù)、擾動(dòng)類型、換熱器整體材料、尺寸結(jié)構(gòu)以及流動(dòng)布置等因素。氣/氣換熱器的主要蓄熱環(huán)節(jié)為金屬壁面的蓄熱,時(shí)間常數(shù)可以表示為金屬壁面的熱慣性時(shí)間常數(shù)。而在氣/液換熱器和液/液換熱器中,冷側(cè)工質(zhì)的比熱容大,其蓄熱能力不能忽略。本文以燃油/滑油換熱器為例,仿真計(jì)算了不同因素對(duì)換熱器動(dòng)態(tài)換熱特性的影響,如圖24~圖26所示。從圖中可知,換熱器質(zhì)量越輕,通道高度尺寸越小,壁厚越小,則時(shí)間常數(shù)越小。

    圖23 轉(zhuǎn)動(dòng)慣量對(duì)發(fā)動(dòng)機(jī)時(shí)間響應(yīng)的影響Fig.23 Effect of moment of inertia on engine time response

    圖24 材料屬性對(duì)換熱器出口溫度的影響Fig.24 Effect of material properties on outlet temperature of heat exchanger

    圖25 流道高度對(duì)換熱器出口溫度的影響Fig.25 Effect of channel height on outlet temperature of heat exchanger

    圖26 壁厚對(duì)換熱器出口溫度的影響Fig.26 Effect of wall thickness on outlet temperature of heat exchanger

    4.3 熱回收利用方法

    目前直升機(jī)上應(yīng)用較多的增壓和防冰供氣源來自發(fā)動(dòng)機(jī)壓氣機(jī)引氣。在飛行高度為1 km,馬赫數(shù)為0.2的工況下,不同引氣量對(duì)發(fā)動(dòng)機(jī)的影響如表6所示。隨著對(duì)發(fā)動(dòng)機(jī)引氣量的減少,發(fā)動(dòng)機(jī)主流氣體溫度均明顯降低,因此,發(fā)動(dòng)機(jī)的散熱量顯著降低。在發(fā)動(dòng)機(jī)輸出軸功率不變的情況下,主流氣體流量隨著引氣量的減少而增大。圖27 說明了引氣量的減少可降低燃油消耗流量,當(dāng)發(fā)動(dòng)機(jī)引氣量為主流流量的4%時(shí),燃油消耗流量比無引氣時(shí)增加了約4%(從0.043 kg/s增加到0.044 8 kg/s);當(dāng)發(fā)動(dòng)機(jī)引氣量為主流流量的10%時(shí),燃油消耗流量比無引氣時(shí)增加了約10%(從0.043 kg/s增加到0.047 5 kg/s)。由此可見,降低引氣量不僅可以降低發(fā)動(dòng)機(jī)的散熱量,還能降低燃油消耗流量,從而降低排氣溫度和紅外輻射強(qiáng)度。

    為了降低發(fā)動(dòng)機(jī)熱影響區(qū)域的固體紅外輻射和熱排氣的氣體輻射,減少引氣量,并利用排氣中的可用熱量,本文設(shè)計(jì)了如圖28所示的熱回收利用熱管理方案,采用非發(fā)動(dòng)機(jī)引氣的環(huán)控系統(tǒng)設(shè)計(jì),并對(duì)尾氣熱、滑油熱和電子設(shè)備熱進(jìn)行回收利用,模式包括制冷模式和制熱模式。在制冷模式中,環(huán)控系統(tǒng)冷媒在座艙的蒸發(fā)器中吸收空氣中的熱量,在電子設(shè)備換熱器中吸收電子設(shè)備液冷循環(huán)中的熱量,通過制冷循環(huán)在冷凝器中將熱量傳遞給外界環(huán)境。此時(shí),熱回收利用循環(huán)不工作。在制熱模式中,熱回收利用循環(huán)開始工作,環(huán)控系統(tǒng)冷媒在蒸發(fā)器和熱回收換熱器中吸收環(huán)境、滑油和尾氣中的熱量,在電子設(shè)備換熱器中吸收電子設(shè)備液冷循環(huán)中的熱量,通過制冷循環(huán)在座艙冷凝器中將熱量傳遞給座艙空氣。熱回收循環(huán)包括滑油熱回收循環(huán)和尾氣熱回收循環(huán),由載冷劑將滑油和尾氣中的熱量輸運(yùn)至熱回收換熱器和防冰部件。這種非發(fā)動(dòng)機(jī)引氣的環(huán)控系統(tǒng)設(shè)計(jì)方案通過耦合多個(gè)換熱器,實(shí)現(xiàn)了環(huán)境控制系統(tǒng)、廢熱回收利用和防冰系統(tǒng)的一體化。

    表6 不同引氣百分比對(duì)發(fā)動(dòng)機(jī)的影響

    圖27 不同引氣百分?jǐn)?shù)對(duì)燃油消耗流量的影響Fig.27 Effect of different percentages of bleed air on fuel consumption mass flow rate

    圖28 熱回收利用熱管理方案Fig.28 Thermal management scheme for heat recovery and utilization

    為了進(jìn)一步說明熱回收利用熱管理方案對(duì)直升機(jī)紅外輻射特性的影響,取原始方案與熱管理方案做比較。其中,原始方案引氣百分比為10%,熱管理方案無引氣,且假設(shè)通過熱管理方案從尾氣中提取50 kW的熱量用于進(jìn)氣道防冰或座艙加溫。此時(shí),原始方案排氣溫度為947 K,熱管理方案排氣溫度為812 K。圖29所示為發(fā)動(dòng)機(jī)熱影響區(qū)域溫度分布,圖30和圖31所示分別為3~5 μm波段和8~14 μm波段機(jī)身固體壁面的紅外輻射通量分布。通過減少引氣和排氣熱回收利用的措施,使得燃油消耗流量降低,從圖中可知發(fā)動(dòng)機(jī)機(jī)匣及其熱影響區(qū)域的溫度明顯降低,紅外輻射通量明顯降低。

    5 結(jié) 論

    1) 根據(jù)直升機(jī)系統(tǒng)特點(diǎn)、熱管理和紅外抑制的需求,分析了影響直升機(jī)溫度和紅外輻射特性的內(nèi)部因素和外部因素,并提出了整機(jī)熱管理和紅外輻射特性耦合仿真和優(yōu)化思路以及需要考慮困難和挑戰(zhàn)。

    2) 基于整機(jī)熱/質(zhì)流和能量平衡關(guān)系,建立了以系統(tǒng)部件為內(nèi)部邊界條件、外環(huán)境為外部邊界條件的直升機(jī)整機(jī)溫度和紅外輻射特性計(jì)算模型。

    3) 直升機(jī)整機(jī)系統(tǒng)呈現(xiàn)3個(gè)層次結(jié)構(gòu),即整機(jī)層、子系統(tǒng)層和部件層,熱源也呈現(xiàn)空間分布的不均勻性和隨時(shí)間變化的特點(diǎn)。本文從空間層面、時(shí)間層面和熱回收利用方法3個(gè)方向出發(fā),以換熱器結(jié)構(gòu)設(shè)計(jì)、系統(tǒng)布局設(shè)計(jì)、部件的時(shí)間響應(yīng)和廢熱回收利用方案設(shè)計(jì)為例分析了整機(jī)熱管理和紅外抑制的優(yōu)化方向。

    圖29 發(fā)動(dòng)機(jī)熱影響區(qū)域的溫度分布Fig.29 Temperature distribution in engine affected region

    圖30 3~5 μm波段紅外輻射通量Fig.30 Infrared radiation flux in 3-5 μm band

    圖31 8~14 μm波段紅外輻射通量Fig.31 Infrared radiation flux in 8-14 μm band

    4) 本研究是指導(dǎo)直升機(jī)整機(jī)系統(tǒng)設(shè)計(jì)、部件選型、熱防護(hù)和紅外抑制設(shè)計(jì)的第一步,可為直升機(jī)綜合熱管理和紅外抑制的設(shè)計(jì)提供依據(jù)。但整機(jī)能量利用策略、部件互聯(lián)的性能可靠性需要進(jìn)一步試驗(yàn)驗(yàn)證。同時(shí),在直升機(jī)整機(jī)層面的一體化系統(tǒng)改造需要建立更全面、更完善的能效評(píng)估方法。

    猜你喜歡
    換熱器直升機(jī)燃油
    直升機(jī)?
    直升機(jī)
    ASM-600油站換熱器的國(guó)產(chǎn)化改進(jìn)
    能源工程(2021年3期)2021-08-05 07:26:14
    燃油泄漏闖了禍
    翅片管式換熱器的傳熱研究進(jìn)展
    制冷(2019年2期)2019-12-09 08:10:36
    310S全焊式板殼型換熱器的焊接
    焊接(2016年10期)2016-02-27 13:05:35
    奔馳S500車燃油表不準(zhǔn)
    邁騰1.8TSI車燃油消耗量大
    奔馳GLA200車燃油表顯示異常
    直升機(jī)取票
    海峽姐妹(2015年3期)2015-02-27 15:09:59
    极品人妻少妇av视频| 久久精品国产a三级三级三级| 波多野结衣av一区二区av| 深夜精品福利| 色精品久久人妻99蜜桃| 狂野欧美激情性xxxx| 亚洲,一卡二卡三卡| 麻豆乱淫一区二区| 熟女av电影| 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 高清视频免费观看一区二区| 日本色播在线视频| 91老司机精品| 妹子高潮喷水视频| 免费看不卡的av| 欧美精品人与动牲交sv欧美| 欧美少妇被猛烈插入视频| 2021少妇久久久久久久久久久| 亚洲专区中文字幕在线| 国产精品.久久久| 国产一区二区三区av在线| 男女高潮啪啪啪动态图| 亚洲欧美一区二区三区国产| 亚洲成人国产一区在线观看 | 秋霞在线观看毛片| 亚洲成人国产一区在线观看 | 中文字幕亚洲精品专区| 国产三级黄色录像| 亚洲欧洲国产日韩| 国产免费一区二区三区四区乱码| 无遮挡黄片免费观看| 人人妻人人澡人人爽人人夜夜| 观看av在线不卡| 黄片播放在线免费| 老汉色∧v一级毛片| 亚洲国产精品999| 久久免费观看电影| 成人18禁高潮啪啪吃奶动态图| 欧美日韩精品网址| 国产91精品成人一区二区三区 | 精品国产一区二区三区久久久樱花| 亚洲人成电影免费在线| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产av蜜桃| 黄色视频不卡| 免费不卡黄色视频| 女性被躁到高潮视频| 国产精品久久久av美女十八| 亚洲一区二区三区欧美精品| 人人澡人人妻人| 国产亚洲一区二区精品| 看免费成人av毛片| av片东京热男人的天堂| 在现免费观看毛片| 成人午夜精彩视频在线观看| 国产精品久久久久久人妻精品电影 | 日韩制服丝袜自拍偷拍| 一级毛片黄色毛片免费观看视频| 18禁裸乳无遮挡动漫免费视频| 久久精品国产综合久久久| 国产深夜福利视频在线观看| 少妇裸体淫交视频免费看高清 | 国产高清videossex| 亚洲欧美日韩另类电影网站| 亚洲久久久国产精品| 亚洲人成电影观看| 日日摸夜夜添夜夜爱| 欧美日韩亚洲高清精品| 午夜福利免费观看在线| 国产成人欧美在线观看 | 中国美女看黄片| 国产精品三级大全| 另类精品久久| 亚洲成色77777| 日韩制服骚丝袜av| 免费在线观看视频国产中文字幕亚洲 | 日韩 欧美 亚洲 中文字幕| 在线看a的网站| 操美女的视频在线观看| 看免费av毛片| 真人做人爱边吃奶动态| 韩国高清视频一区二区三区| 涩涩av久久男人的天堂| videos熟女内射| 在线观看一区二区三区激情| 赤兔流量卡办理| 日本黄色日本黄色录像| 欧美日韩亚洲综合一区二区三区_| 男女边摸边吃奶| 久久精品久久久久久噜噜老黄| 免费久久久久久久精品成人欧美视频| 50天的宝宝边吃奶边哭怎么回事| 日韩免费高清中文字幕av| 狂野欧美激情性bbbbbb| 纯流量卡能插随身wifi吗| 国产免费又黄又爽又色| 国产成人一区二区在线| 欧美激情高清一区二区三区| 亚洲天堂av无毛| 中文字幕色久视频| 夫妻午夜视频| 国产高清videossex| 热re99久久精品国产66热6| 91精品三级在线观看| 日本色播在线视频| 国产一区二区在线观看av| 久久综合国产亚洲精品| 黄色a级毛片大全视频| 天天躁日日躁夜夜躁夜夜| 一级黄色大片毛片| 青青草视频在线视频观看| 亚洲熟女精品中文字幕| 国产精品国产三级专区第一集| 一区二区av电影网| 中文字幕最新亚洲高清| 99久久99久久久精品蜜桃| 97精品久久久久久久久久精品| 伊人久久大香线蕉亚洲五| 久热这里只有精品99| 欧美激情极品国产一区二区三区| av网站免费在线观看视频| 日本wwww免费看| 国产视频首页在线观看| 亚洲欧美清纯卡通| 国产日韩欧美在线精品| 后天国语完整版免费观看| 亚洲精品av麻豆狂野| av在线app专区| 亚洲欧美精品自产自拍| 9色porny在线观看| 亚洲专区中文字幕在线| 国产免费福利视频在线观看| 精品一区二区三区av网在线观看 | 日韩av免费高清视频| 777久久人妻少妇嫩草av网站| 热99久久久久精品小说推荐| 后天国语完整版免费观看| 王馨瑶露胸无遮挡在线观看| 一区福利在线观看| 日本午夜av视频| 亚洲av片天天在线观看| a级片在线免费高清观看视频| 亚洲熟女精品中文字幕| 美女国产高潮福利片在线看| 欧美日本中文国产一区发布| 777米奇影视久久| 97精品久久久久久久久久精品| 国产一区二区三区综合在线观看| 爱豆传媒免费全集在线观看| 欧美日韩黄片免| 国产免费现黄频在线看| 国产激情久久老熟女| 亚洲精品中文字幕在线视频| 满18在线观看网站| 久久 成人 亚洲| videos熟女内射| 日韩一本色道免费dvd| 看免费av毛片| 国产在线观看jvid| 女性生殖器流出的白浆| 成人免费观看视频高清| 精品一区二区三卡| 亚洲中文字幕日韩| 精品一区二区三卡| 欧美人与善性xxx| 国产成人一区二区在线| 99国产精品一区二区蜜桃av | 国产成人精品久久久久久| 母亲3免费完整高清在线观看| 国产成人91sexporn| 精品久久久精品久久久| 免费少妇av软件| 国产成人精品久久二区二区91| 色婷婷av一区二区三区视频| 午夜免费鲁丝| 午夜两性在线视频| 国产日韩欧美在线精品| 叶爱在线成人免费视频播放| 操出白浆在线播放| 人人澡人人妻人| 亚洲欧美一区二区三区国产| 久久久精品区二区三区| 自线自在国产av| 欧美精品高潮呻吟av久久| 精品一区二区三卡| 亚洲精品一二三| 老汉色av国产亚洲站长工具| 五月开心婷婷网| av有码第一页| 叶爱在线成人免费视频播放| 精品免费久久久久久久清纯 | 精品国产乱码久久久久久男人| 免费观看av网站的网址| 国产成人精品久久二区二区91| 美女大奶头黄色视频| 欧美人与性动交α欧美精品济南到| 日本午夜av视频| 国产成人91sexporn| 大片电影免费在线观看免费| 免费久久久久久久精品成人欧美视频| 在线观看国产h片| 亚洲情色 制服丝袜| 91老司机精品| 欧美日韩国产mv在线观看视频| 水蜜桃什么品种好| 亚洲欧美色中文字幕在线| 午夜福利在线免费观看网站| 只有这里有精品99| 欧美另类一区| 最近手机中文字幕大全| 亚洲 国产 在线| 国产成人系列免费观看| 久久精品国产亚洲av涩爱| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美日韩另类电影网站| 无遮挡黄片免费观看| 99re6热这里在线精品视频| 国产伦理片在线播放av一区| 51午夜福利影视在线观看| 日本色播在线视频| 91老司机精品| 国产精品熟女久久久久浪| 欧美 日韩 精品 国产| 国产精品久久久av美女十八| www.999成人在线观看| 伊人亚洲综合成人网| 少妇的丰满在线观看| 国精品久久久久久国模美| 99国产精品一区二区三区| 在线观看www视频免费| 精品人妻一区二区三区麻豆| 青春草视频在线免费观看| 国产欧美日韩一区二区三 | 又黄又粗又硬又大视频| 国产精品一区二区精品视频观看| 韩国高清视频一区二区三区| 自线自在国产av| 天天影视国产精品| 18禁国产床啪视频网站| www.av在线官网国产| 天天操日日干夜夜撸| 欧美日韩福利视频一区二区| 在线av久久热| 精品少妇内射三级| 欧美成人精品欧美一级黄| 熟女av电影| 亚洲七黄色美女视频| 欧美xxⅹ黑人| 国产精品一国产av| 欧美日韩亚洲综合一区二区三区_| 亚洲精品日韩在线中文字幕| 国产成人欧美在线观看 | 免费观看a级毛片全部| 免费不卡黄色视频| 欧美另类一区| 高潮久久久久久久久久久不卡| 国产免费视频播放在线视频| 女人爽到高潮嗷嗷叫在线视频| 午夜福利影视在线免费观看| 在线 av 中文字幕| 又紧又爽又黄一区二区| 美女视频免费永久观看网站| 少妇精品久久久久久久| 国产不卡av网站在线观看| 中文字幕人妻丝袜制服| 亚洲美女黄色视频免费看| 日韩一卡2卡3卡4卡2021年| 国产男女内射视频| 国精品久久久久久国模美| 亚洲午夜精品一区,二区,三区| 性色av一级| 久久久国产欧美日韩av| 男女床上黄色一级片免费看| www.av在线官网国产| 男的添女的下面高潮视频| 两人在一起打扑克的视频| 高清不卡的av网站| 精品少妇一区二区三区视频日本电影| 久久毛片免费看一区二区三区| 久久天躁狠狠躁夜夜2o2o | 亚洲国产精品国产精品| 日韩熟女老妇一区二区性免费视频| 久久精品成人免费网站| 国产精品秋霞免费鲁丝片| 久久久精品94久久精品| 欧美日韩国产mv在线观看视频| 久久精品国产a三级三级三级| 夫妻性生交免费视频一级片| 波多野结衣一区麻豆| 国产男女内射视频| 亚洲中文日韩欧美视频| 精品福利观看| 嫩草影视91久久| 亚洲第一青青草原| 在线看a的网站| 国产97色在线日韩免费| 精品少妇久久久久久888优播| 国产精品免费视频内射| 一二三四社区在线视频社区8| 欧美精品人与动牲交sv欧美| 国产视频首页在线观看| 在线亚洲精品国产二区图片欧美| 90打野战视频偷拍视频| 成年女人毛片免费观看观看9 | 老司机影院毛片| 精品欧美一区二区三区在线| 亚洲自偷自拍图片 自拍| 99国产精品一区二区三区| 悠悠久久av| 国产不卡av网站在线观看| 黄频高清免费视频| 亚洲人成电影免费在线| 久久久久精品国产欧美久久久 | 最黄视频免费看| 国产男女内射视频| 最近手机中文字幕大全| 999久久久国产精品视频| 亚洲,欧美精品.| 日本a在线网址| 久9热在线精品视频| av国产精品久久久久影院| 亚洲伊人色综图| 国产伦人伦偷精品视频| 国产男人的电影天堂91| 少妇粗大呻吟视频| 午夜福利,免费看| 亚洲一码二码三码区别大吗| 首页视频小说图片口味搜索 | 久久亚洲精品不卡| 丝袜美腿诱惑在线| 性少妇av在线| 亚洲欧美清纯卡通| 精品亚洲成a人片在线观看| 国产三级黄色录像| 亚洲国产av新网站| 午夜日韩欧美国产| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品大桥未久av| 亚洲免费av在线视频| netflix在线观看网站| 亚洲精品一区蜜桃| 天堂中文最新版在线下载| 亚洲国产中文字幕在线视频| 午夜两性在线视频| 亚洲久久久国产精品| 亚洲av男天堂| 午夜福利影视在线免费观看| 午夜精品国产一区二区电影| 国产欧美日韩一区二区三 | 婷婷丁香在线五月| 大片电影免费在线观看免费| svipshipincom国产片| 国产视频首页在线观看| 天堂8中文在线网| 啦啦啦在线免费观看视频4| netflix在线观看网站| 狂野欧美激情性bbbbbb| 伊人久久大香线蕉亚洲五| 女人高潮潮喷娇喘18禁视频| 在线观看人妻少妇| 精品国产一区二区三区四区第35| 精品国产一区二区三区四区第35| 一区二区日韩欧美中文字幕| 久久天堂一区二区三区四区| 人妻一区二区av| 黄色视频不卡| 老司机影院毛片| 成年美女黄网站色视频大全免费| cao死你这个sao货| 亚洲精品久久久久久婷婷小说| 69精品国产乱码久久久| 两人在一起打扑克的视频| av一本久久久久| 亚洲欧美一区二区三区久久| 热re99久久精品国产66热6| 欧美老熟妇乱子伦牲交| 国产熟女午夜一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 国产免费又黄又爽又色| 高清黄色对白视频在线免费看| 亚洲国产欧美在线一区| 国产黄色视频一区二区在线观看| 日韩,欧美,国产一区二区三区| 91麻豆av在线| 在线观看www视频免费| 久久亚洲精品不卡| 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 97精品久久久久久久久久精品| 精品人妻在线不人妻| 亚洲黑人精品在线| 国产片内射在线| 欧美成人午夜精品| 亚洲av美国av| 观看av在线不卡| 91精品伊人久久大香线蕉| 亚洲av在线观看美女高潮| 免费高清在线观看日韩| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 久久久久久久国产电影| 久久久国产精品麻豆| 国产麻豆69| 又黄又粗又硬又大视频| 在线观看一区二区三区激情| 纵有疾风起免费观看全集完整版| 一级片'在线观看视频| 下体分泌物呈黄色| 亚洲视频免费观看视频| 欧美成狂野欧美在线观看| 久久久久久久精品精品| 男人舔女人的私密视频| 欧美日韩综合久久久久久| www.精华液| 久久精品久久久久久久性| 国产精品.久久久| www.熟女人妻精品国产| 国产高清视频在线播放一区 | 亚洲,欧美精品.| 极品人妻少妇av视频| 国产成人精品久久二区二区免费| 十八禁高潮呻吟视频| 精品免费久久久久久久清纯 | 一本—道久久a久久精品蜜桃钙片| 免费一级毛片在线播放高清视频 | 日韩一区二区三区影片| 男女床上黄色一级片免费看| 黄色片一级片一级黄色片| 丝袜美足系列| 免费少妇av软件| 这个男人来自地球电影免费观看| 高潮久久久久久久久久久不卡| 亚洲专区中文字幕在线| 国产伦人伦偷精品视频| 男女高潮啪啪啪动态图| 欧美亚洲日本最大视频资源| 国产成人精品在线电影| √禁漫天堂资源中文www| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| 丰满迷人的少妇在线观看| 国产在线一区二区三区精| 久久精品亚洲av国产电影网| 国产欧美日韩一区二区三区在线| 91字幕亚洲| 亚洲欧美精品自产自拍| 国产片内射在线| 国产精品久久久久久精品古装| 看十八女毛片水多多多| 91麻豆av在线| 亚洲av综合色区一区| 国产精品av久久久久免费| 高清av免费在线| 国产精品久久久久久精品古装| 精品国产一区二区三区久久久樱花| 日本猛色少妇xxxxx猛交久久| 亚洲精品国产区一区二| 免费观看av网站的网址| 亚洲情色 制服丝袜| 搡老岳熟女国产| 免费在线观看黄色视频的| 少妇 在线观看| 国产男女超爽视频在线观看| 色视频在线一区二区三区| 亚洲第一av免费看| 日韩伦理黄色片| 久热爱精品视频在线9| 女性生殖器流出的白浆| 国产精品三级大全| 亚洲专区中文字幕在线| 麻豆av在线久日| av国产久精品久网站免费入址| 黄片小视频在线播放| 可以免费在线观看a视频的电影网站| 亚洲精品国产区一区二| 美女大奶头黄色视频| 国产成人精品久久二区二区免费| 人人妻人人爽人人添夜夜欢视频| 国产精品人妻久久久影院| 久久久久久久国产电影| 国产色视频综合| 一边摸一边做爽爽视频免费| 国产一区二区激情短视频 | 在线天堂中文资源库| 大片电影免费在线观看免费| 亚洲精品第二区| 在线观看免费视频网站a站| 欧美成狂野欧美在线观看| 91精品国产国语对白视频| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| 久久ye,这里只有精品| 黑人巨大精品欧美一区二区蜜桃| 久久久亚洲精品成人影院| 国产成人精品久久二区二区免费| 精品人妻在线不人妻| 日韩制服骚丝袜av| 一级黄片播放器| 国产熟女欧美一区二区| 青青草视频在线视频观看| 亚洲欧美色中文字幕在线| 69精品国产乱码久久久| 欧美日韩精品网址| 人人妻人人澡人人爽人人夜夜| 免费av中文字幕在线| 人人妻,人人澡人人爽秒播 | 亚洲欧洲日产国产| 久久人人97超碰香蕉20202| 18在线观看网站| 欧美日韩av久久| 国产成人精品久久久久久| 热re99久久国产66热| 老汉色av国产亚洲站长工具| av一本久久久久| 国产淫语在线视频| 狂野欧美激情性xxxx| 亚洲伊人色综图| 99香蕉大伊视频| 九色亚洲精品在线播放| 欧美日韩亚洲综合一区二区三区_| 久久久久久久精品精品| 国产午夜精品一二区理论片| 91麻豆av在线| 人人妻,人人澡人人爽秒播 | 在线av久久热| 在线亚洲精品国产二区图片欧美| 1024视频免费在线观看| 一本大道久久a久久精品| 男女床上黄色一级片免费看| 一级片'在线观看视频| 久久精品国产a三级三级三级| 免费人妻精品一区二区三区视频| www日本在线高清视频| 色94色欧美一区二区| 久久久久久久久免费视频了| 日韩精品免费视频一区二区三区| 性少妇av在线| 欧美 日韩 精品 国产| 国产黄色免费在线视频| 老司机影院成人| 美女脱内裤让男人舔精品视频| 欧美日韩视频高清一区二区三区二| 欧美黑人精品巨大| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 亚洲人成电影观看| 久久久久网色| 中文字幕人妻丝袜制服| 中文字幕高清在线视频| 人妻一区二区av| 男女国产视频网站| 久久狼人影院| 亚洲成国产人片在线观看| av线在线观看网站| 最黄视频免费看| 天堂俺去俺来也www色官网| 亚洲黑人精品在线| 国产精品熟女久久久久浪| 99国产综合亚洲精品| 高清不卡的av网站| 宅男免费午夜| 亚洲欧洲日产国产| 精品一区二区三卡| 狂野欧美激情性xxxx| 日韩精品免费视频一区二区三区| 免费av中文字幕在线| 日韩中文字幕欧美一区二区 | 日本av免费视频播放| 热re99久久精品国产66热6| 欧美日韩福利视频一区二区| 乱人伦中国视频| 99久久精品国产亚洲精品| 亚洲图色成人| 极品人妻少妇av视频| 老司机深夜福利视频在线观看 | 好男人视频免费观看在线| 人人妻人人爽人人添夜夜欢视频| 一级毛片电影观看| 赤兔流量卡办理| 国产精品三级大全| av视频免费观看在线观看| svipshipincom国产片| 久久精品aⅴ一区二区三区四区| 大码成人一级视频| 男女午夜视频在线观看| 精品久久久久久电影网| www.自偷自拍.com| 国产国语露脸激情在线看| 18禁观看日本| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 十八禁人妻一区二区| 亚洲精品在线美女| 中文欧美无线码| 精品国产超薄肉色丝袜足j| 十八禁人妻一区二区| 在线观看www视频免费| 久久九九热精品免费| 亚洲欧美日韩另类电影网站| 婷婷色麻豆天堂久久| 国产老妇伦熟女老妇高清| 男女免费视频国产| 一二三四在线观看免费中文在| 日韩av不卡免费在线播放| 韩国高清视频一区二区三区| 麻豆国产av国片精品| 韩国精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 精品人妻一区二区三区麻豆| 欧美黄色淫秽网站| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新|