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

    不同后進(jìn)氣系統(tǒng)對雙面復(fù)合葉輪的影響

    2022-10-25 05:04:34羅炳亮楊蓓王云張呈
    航空工程進(jìn)展 2022年5期
    關(guān)鍵詞:總溫橢圓型進(jìn)氣管

    羅炳亮,楊蓓,王云,張呈

    (1.南昌航空大學(xué) 飛行器工程學(xué)院,南昌 330063)(2.臺(tái)州學(xué)院 航空工程學(xué)院,臺(tái)州 318000)

    0 引言

    微型燃?xì)廨啓C(jī)在航空航天、地面設(shè)備等動(dòng)力裝置上應(yīng)用廣泛,一直是研究的熱點(diǎn)。國內(nèi)外對微型燃?xì)廨啓C(jī)的各個(gè)方向的研究不勝枚舉,但近幾十年來技術(shù)發(fā)展趨于平緩,亟待有原理和結(jié)構(gòu)上的創(chuàng)新來實(shí)現(xiàn)高性能、高循環(huán)參數(shù)的突破。雙面復(fù)合葉輪具有更高的氣動(dòng)性能,有望為高效費(fèi)比的微燃機(jī)提供必要的技術(shù)儲(chǔ)備。

    有關(guān)常規(guī)葉輪進(jìn)氣畸變的研究已成為近年來的研究熱點(diǎn),國外,I.Ariga等采用實(shí)驗(yàn)方法證明了進(jìn)氣的總壓畸變會(huì)使葉輪的性能惡化,其中周向畸變影響力最大;A.Engeda等發(fā)現(xiàn)在進(jìn)氣管道內(nèi)設(shè)置整流板可以改善葉輪進(jìn)口流場。國內(nèi),周頌東等較早地通過試驗(yàn)對比了有無進(jìn)氣畸變的常規(guī)葉輪,發(fā)現(xiàn)徑向組合畸變拓寬了葉輪工作的穩(wěn)定范圍,改善了其性能,且隨轉(zhuǎn)速上升愈加明顯;李杜等、老大中等分別對帶有90°和180°彎管的常規(guī)葉輪進(jìn)行了大量的試驗(yàn)與數(shù)值模擬研究,得出安裝周向位置、軸向位置的最優(yōu)方案,發(fā)現(xiàn)彎管的總壓畸變與下游的蝸殼傳播上來的周向總壓畸變耦合是造成葉輪性能下降的原因等一系列豐碩成果。

    復(fù)雜的二次流動(dòng)和局部溫升情況嚴(yán)重影響著葉輪進(jìn)氣條件和氣動(dòng)性能,但目前國內(nèi)外對微燃機(jī)的進(jìn)氣畸變研究主要針對常規(guī)葉輪進(jìn)氣通道,對新型雙面復(fù)合葉輪的進(jìn)氣鮮少涉及。由于雙面復(fù)合葉輪采用雙向進(jìn)氣,即除了與傳統(tǒng)葉輪具有相同的前進(jìn)氣通道外,還存在后進(jìn)氣通道,而后進(jìn)氣通道與前進(jìn)氣通道完全不同,因此有必要闡明后進(jìn)氣管對雙面復(fù)合葉輪的影響機(jī)理。

    本文通過對三種形狀后進(jìn)氣管進(jìn)行壓縮特性及流場分析,對比不同后進(jìn)氣管的優(yōu)缺點(diǎn),以期為后進(jìn)氣系統(tǒng)設(shè)計(jì)提供借鑒。

    1 研究模型和研究方法

    1.1 雙面復(fù)合葉輪模型

    本文研究的雙面復(fù)合葉輪模型是由雙面常規(guī)葉輪與空心渦輪葉片徑向連接為一體的,是一個(gè)鏡像幾何對稱的雙面常規(guī)葉輪,而空心渦輪葉片在此對稱基面上并不對稱,如圖1所示,其主要設(shè)計(jì)參數(shù)如表1所示。

    圖1 雙面復(fù)合葉輪3D模型Fig.1 3D model of double-sided composite impeller

    表1 雙面復(fù)合葉輪設(shè)計(jì)參數(shù)Table 1 Design parameters of double-sided composite impeller

    雙面復(fù)合葉輪進(jìn)氣系統(tǒng)一端軸向進(jìn)氣,另一端則經(jīng)過較為彎曲的通道實(shí)現(xiàn)進(jìn)氣,如圖2所示,將軸向進(jìn)氣的葉輪定義為主葉輪,反之為副葉輪。

    圖2 軸輸出功的新型微燃機(jī)Fig.2 New micro gas turbine with shaft output work

    1.2 網(wǎng)格劃分及邊界條件

    不同幾何形狀后進(jìn)氣管的雙面復(fù)合葉輪各內(nèi)流場均采用meshing劃分網(wǎng)格并拼接成一體。網(wǎng)格的參考單元畸變度最大值均在0.94以下,表明該計(jì)算模型可用于計(jì)算仿真。其中雙面復(fù)合葉輪的網(wǎng)格數(shù)目約有266萬,含直管型后進(jìn)氣管的雙面復(fù)合葉輪341萬、對稱翼型后進(jìn)氣管的雙面復(fù)合葉輪765萬、橢圓型后進(jìn)氣管的雙面復(fù)合葉輪759萬。

    采用CFX(Computational Fuid X)對不同幾何形狀后進(jìn)氣管的雙面復(fù)合葉輪進(jìn)行數(shù)值模擬,基于有限體積法與自帶壁面函數(shù)的標(biāo)準(zhǔn)k-ε湍流模型對三維雷諾平均Navier-Stokes方程進(jìn)行近似求解,其中熱傳導(dǎo)模型采用Total Energy,轉(zhuǎn)速設(shè)置為設(shè)計(jì)轉(zhuǎn)速,計(jì)算殘差到10以下確認(rèn)為收斂。對軸向進(jìn)氣及后進(jìn)氣的雙面復(fù)合葉輪采用進(jìn)口總溫總壓出口靜壓邊界條件,對于后進(jìn)氣管氣熱耦合情況,外通道設(shè)置進(jìn)口總溫、流量出口靜壓,通過調(diào)整空心渦輪葉片出口靜壓值得到總流量,從而糾正外通道進(jìn)口條件,計(jì)算出不同工況點(diǎn)數(shù)據(jù)。含后進(jìn)氣管的雙面復(fù)合葉輪計(jì)算模型如圖3所示。

    圖3 含后進(jìn)氣管的雙面復(fù)合葉輪計(jì)算模型Fig.3 Calculation model of double-sided composite impeller with rear intake pipe

    1.3 網(wǎng)格無關(guān)性驗(yàn)證

    為檢驗(yàn)網(wǎng)格是否會(huì)影響數(shù)值模擬的計(jì)算結(jié)果,采用5種不同的網(wǎng)格尺度對雙面復(fù)合葉輪計(jì)算模型進(jìn)行模擬,如圖4所示,可以看出:當(dāng)網(wǎng)格保持在266萬以上時(shí),設(shè)計(jì)工況下的總壓比和等熵效率數(shù)值趨于平穩(wěn),總壓比和等熵效率數(shù)值隨著網(wǎng)格數(shù)量的增大,誤差在5%以內(nèi)。因此綜合計(jì)算精度和計(jì)算量的需要,可采用266萬網(wǎng)格數(shù)的雙面復(fù)合葉輪計(jì)算模型進(jìn)行模擬。

    圖4 雙面復(fù)合葉輪網(wǎng)格無關(guān)性驗(yàn)證Fig.4 Grid independence verification of double-sided composite impeller

    2 后進(jìn)氣管對雙面復(fù)合葉輪的影響

    鑒于雙面復(fù)合葉輪(雙面常規(guī)葉輪共享空心渦輪葉片)不同于常規(guī)葉輪獨(dú)自享有擴(kuò)壓器,取常規(guī)葉輪出口為出口計(jì)算截面進(jìn)行分析。

    不同后進(jìn)氣管道的雙面復(fù)合葉輪性能比較如圖5所示,可以看出:不同幾何形狀的后進(jìn)氣管讓雙面復(fù)合葉輪的流量范圍達(dá)到了理想要求,而對稱翼型和橢圓型后進(jìn)氣管對葉輪的性能影響差異很小,兩者相比直管型后進(jìn)氣管的雙面復(fù)合葉輪,性能在靠近失速點(diǎn)的小流量范圍內(nèi)差異性最大。

    圖5 不同后進(jìn)氣管的雙面復(fù)合葉輪性能比較Fig.5 Performance comparison of double-sided composite impellers with different rear intake pipes

    不同進(jìn)氣畸變的副葉輪對比如圖6所示,可以看出:進(jìn)口總壓隨流量增大而減少,影響葉輪性能,同時(shí)圖6驗(yàn)證了進(jìn)口總溫的畸變對葉輪性能的影響程度更大。

    圖6 不同進(jìn)氣畸變的副葉輪對比Fig.6 Comparison of auxiliary impellers with different inlet distortion

    雙面復(fù)合葉輪質(zhì)量流量比與折合質(zhì)量流量的特性曲線如圖7所示,質(zhì)量流量比定義為主葉輪氣體流量與副葉輪氣體質(zhì)量流量的比值。

    圖7 流量比分布特性曲線Fig.7 Flow ratio distribution characteristic curve

    從圖7可以看出:隨折合流量減少,后進(jìn)氣直管的雙面復(fù)合葉輪流量比逐步增大,說明對于幾何形狀一致的主/副葉輪,由于出口空心渦輪形狀的影響,對副葉輪而言,對應(yīng)的空心渦輪葉片出口面積小,有憋背壓的作用,副葉輪趨于小流量;對稱翼型和橢圓型后進(jìn)氣管的雙面復(fù)合葉輪流量比隨折合流量的減少而顯著增大,即副葉輪更趨于小流量,這不僅是下游空心渦輪葉片的影響,也可能是副葉輪由于徑向進(jìn)氣方式所帶來的進(jìn)氣畸變造成的,進(jìn)氣畸變使得副葉輪流場復(fù)雜。另一方面,流量的減少對應(yīng)著進(jìn)口絕對速度減少,由速度三角形可知,副葉輪進(jìn)口正攻角增大,副葉輪要比主葉輪提前發(fā)生喘振,使得主副葉輪工作模式發(fā)生轉(zhuǎn)變,這在文獻(xiàn)[21]中也有所提及。因此,選取近失速與近堵塞工況點(diǎn)對不同幾何形狀后進(jìn)氣管的雙面復(fù)合葉輪進(jìn)行分析。

    3 后進(jìn)氣管與副葉輪流場結(jié)構(gòu)研究

    3.1 后進(jìn)氣管流場

    鑒于后進(jìn)氣管是一個(gè)軸對稱模型,為方便研究,取對稱翼型與橢圓型管道對稱中心面為剖面進(jìn)行分析,從左到右依次是直管型、對稱翼型、橢圓型,近堵塞點(diǎn)后進(jìn)氣管各參數(shù)云圖如圖8所示。

    圖8 近堵塞點(diǎn)后進(jìn)氣管剖面各參數(shù)云圖Fig.8 Cloud diagram of parameters of inlet pipe profile near the plug point

    從圖8(a)可以看出:由于副葉輪的抽吸作用,使得離副葉輪軸向位置越近的抽吸效果越明顯,且因?yàn)楦比~輪葉尖速度的擾動(dòng)影響,也讓后進(jìn)氣管對應(yīng)出口處速度增大。反之,在后進(jìn)氣管外側(cè)區(qū)域出現(xiàn)低速區(qū),橢圓型后進(jìn)氣管相對對稱翼型后進(jìn)氣管流場較佳。這在圖8(b)中也有所體現(xiàn),此外相對對稱翼型后進(jìn)氣管,橢圓型進(jìn)氣管受進(jìn)口輪轂的滯止氣流影響較弱,管內(nèi)流場均勻,這在圖8(c)、圖8(d)中也得以驗(yàn)證。

    對稱翼型后進(jìn)氣管和橢圓型后進(jìn)氣管外通道(尾噴管通道)的馬赫數(shù)云圖如圖8(e)所示,可以看出:對稱翼型后進(jìn)氣管外通道喉道已出現(xiàn)超聲現(xiàn)象,易出現(xiàn)激波干擾附面層損失。而橢圓型后進(jìn)氣管外通道喉道未出現(xiàn)超聲現(xiàn)象,且流場分布均勻。就近失速點(diǎn)而言,與上同理,不再贅述。

    這兩種后進(jìn)氣管方案對雙面復(fù)合葉輪性能影響并無明顯差異(從圖5可知),加之考慮到尾噴管流動(dòng)損失直接關(guān)系到發(fā)動(dòng)機(jī)推力的好壞。因此考慮采用橢圓型后進(jìn)氣管方案作為后進(jìn)氣系統(tǒng)。

    為進(jìn)一步分析進(jìn)氣方式所帶來的副葉輪流場流動(dòng)的變化,下文均以橢圓型后進(jìn)氣管的雙面復(fù)合葉輪為對象展開分析。

    3.2 副葉輪進(jìn)口流場

    3.2.1 直管進(jìn)氣時(shí)副葉輪進(jìn)口流場

    直管進(jìn)氣時(shí)副葉輪進(jìn)口5個(gè)葉高位置總溫和總壓分布如圖9所示。

    圖9 近堵塞點(diǎn)副葉輪進(jìn)口不同葉高流場參數(shù)分布(直管型后進(jìn)氣管)Fig.9 Distribution of flow field parameters at different blade heights at the inlet of auxiliary impeller near the plug point(straight rear inlet)

    從圖9可以看出:隨葉高的上升,總溫沿周向分布差異性增大,且呈周期性正弦曲線變化更明顯。這種波動(dòng)振幅正好對應(yīng)著副葉輪各主流道位置(R1~R14),說明這是受葉輪葉片葉尖的擾動(dòng)所造成的結(jié)果,總溫正負(fù)峰值也在285~300 K之間波動(dòng)。分析總壓分布時(shí)可見,葉輪進(jìn)口處的總壓沿周向也呈現(xiàn)周期性正弦曲線變化,但由于受葉輪輪轂滯止氣流及葉輪抽吸作用,葉根處的周向總壓比葉中段高,總壓得到了恢復(fù)。同時(shí)受空心渦輪葉片的壓力梯度傳播至上游的影響,使得葉尖處周向總壓波動(dòng)更加明顯。

    綜上,受到葉輪葉片分布位置及空心渦輪葉片壓力梯度的影響,總溫和總壓在周向分布上呈現(xiàn)周期性正弦曲線現(xiàn)象,且隨葉高上升,波動(dòng)情況更加明顯。

    3.2.2 橢圓型后進(jìn)氣管進(jìn)氣時(shí)副葉輪進(jìn)口流場

    為方便研究副葉輪進(jìn)口流場情況,進(jìn)口葉輪近堵塞點(diǎn)和近失速點(diǎn)五種葉高位置的各參數(shù)分布分別如圖10~圖11所示,其中比較同流量下不同進(jìn)氣管給副葉輪進(jìn)口流場帶來的影響。

    圖10 近堵塞點(diǎn)副葉輪進(jìn)口不同葉高流場參數(shù)分布(橢圓型后進(jìn)氣管)Fig.10 Distribution of flow field parameters at different blade heights at the inlet of auxiliary impeller near the plug point(elliptical rear inlet)

    圖11 近失速點(diǎn)副葉輪進(jìn)口流場參數(shù)分布Fig.11 Parameter distribution of flow field at inlet of auxiliary impeller near stall point

    從圖10可以看出:相對于直管進(jìn)氣的副葉輪,橢圓型后進(jìn)氣管的副葉輪進(jìn)口周向與徑向總溫畸變變化明顯,受后進(jìn)氣管對流換熱影響,進(jìn)口除葉中處外其他處的總溫均得到上升且徑向總溫畸變增加,說明在后進(jìn)氣管換熱時(shí),越靠近后進(jìn)氣管壁面位置的后進(jìn)氣氣流換熱越明顯,反之進(jìn)氣管的中部氣流無明顯溫升;而由于后進(jìn)氣管進(jìn)氣周向位置與副葉輪葉片周向位置的疊加效果,使得總溫的正峰值高近325 K。從圖10(b)可以看出:受后進(jìn)氣管彎道的影響,進(jìn)口葉根處的總壓恢復(fù)減少并開始惡化,同時(shí)后進(jìn)氣管的徑向進(jìn)氣結(jié)構(gòu)及彎道的沿程損失讓靠近葉尖處的總壓周向分布波動(dòng)更大,整個(gè)進(jìn)口流場參數(shù)分布不均勻化加劇。

    進(jìn)口的總溫、總壓組合畸變與空心渦輪葉片壓力梯度向上游傳播相耦合,共同改變了副葉輪進(jìn)口流場的結(jié)構(gòu)。

    從圖11可以看出:在近失速點(diǎn)的小流量工況下,徑向總溫畸變變化較小,而在0.7葉高處的周向總溫畸變卻整體升高,說明在小流量工況下,更多的氣流參與了與后進(jìn)氣管的壁面換熱,且管中部的氣流占比小,冷卻效果下降。從總壓分布可知,小流量工況下的各葉高處周向總壓分布波動(dòng)下降,徑向總壓畸變得到緩解。

    3.2.3 副葉輪進(jìn)口流場對其性能的影響

    直管與橢圓型后進(jìn)氣管近堵塞點(diǎn)副葉輪進(jìn)口不同葉軸向速度分布如圖12所示。

    圖12 直管與橢圓型后進(jìn)氣管近堵塞點(diǎn)副葉輪進(jìn)口不同葉軸向速度分布Fig.12 Axial velocity distribution of different blades at the inlet of auxiliary impeller near the plug point of straight pipe and elliptical rear inlet pipe

    從圖12可以看出:直管進(jìn)氣的副葉輪軸向速度在周向分布上周期正弦曲線規(guī)律明顯,同時(shí)葉尖區(qū)域波動(dòng)最大;相比于直管進(jìn)氣,橢圓型后進(jìn)氣管的副葉輪在葉尖區(qū)域的周向分布波動(dòng)平緩,而近葉根部位的軸向速度整體下滑,對應(yīng)著上述副葉輪進(jìn)口總溫總壓組合畸變?nèi)~根處?;谒俣热切卧?,在轉(zhuǎn)速與葉輪進(jìn)口流量不變的條件下,有且僅當(dāng)軸向速度變化時(shí),才能夠改變進(jìn)口相對氣流角。因此,受橢圓型后進(jìn)氣管結(jié)構(gòu)及氣熱耦合的影響,副葉輪葉根處進(jìn)口軸向速度下降,因此進(jìn)口處的相對氣流角增大,攻角減少,副葉輪性能惡化。

    與近堵塞點(diǎn)工況不同,直管與橢圓型后進(jìn)氣管近失速點(diǎn)副葉輪進(jìn)口不同葉軸向速度分布如圖13所示,可以看出:近失速點(diǎn)工況下,直管進(jìn)氣的副葉輪軸向速度依次沿葉高上升逐步下降,說明由于直管壁面的沿程損失產(chǎn)生的邊界層分離渦,使得副葉輪葉尖處相對氣流角增大,影響葉輪效率。與直管的副葉輪軸向速度比較發(fā)現(xiàn),橢圓型后進(jìn)氣管進(jìn)氣會(huì)造成副葉輪進(jìn)口葉根處軸向速度下降并影響至葉中處,同時(shí)也說明,在近失速點(diǎn)附近,后進(jìn)氣管結(jié)構(gòu)引起的總溫總壓組合畸變造成副葉輪進(jìn)口葉根處相對氣流角增大才是導(dǎo)致副葉輪性能下降的原因。

    圖13 直管與橢圓型后進(jìn)氣管近失速點(diǎn)副葉輪進(jìn)口不同葉軸向速度分布Fig.13 Axial velocity distribution of different blades at the inlet of auxiliary impeller near stall point of straight pipe and elliptical rear inlet pipe

    綜上所述,副葉輪進(jìn)口流場受橢圓型后進(jìn)氣管結(jié)構(gòu)帶來的總溫總壓組合畸變的影響,導(dǎo)致葉根處的軸向速度減少,這也是副葉輪性能下降的主要原因。

    3.3 副葉輪流場流動(dòng)研究

    為了更直觀地了解葉輪流動(dòng)的機(jī)理,基于邊界渦量動(dòng)力學(xué),選取經(jīng)典軸向邊界渦量流物理量進(jìn) 行 分 析。近 堵 塞 點(diǎn) 副 葉 輪BVFz(Axial Boundary Vorticity Flux)分布如圖14所示,可以看出:直管進(jìn)氣的副葉輪前緣及尾緣葉尖區(qū)域出現(xiàn)BVFz峰值,說明該區(qū)域壓力梯度過大,產(chǎn)生了較多的壁面分離渦,這也是誘導(dǎo)葉尖泄露渦及二次流的主要原因;與直管進(jìn)氣比較,橢圓型后進(jìn)氣管的副葉輪表面的BVFz峰值已出現(xiàn)在前緣葉根部,說明該處也存在分離流動(dòng)的情況,副葉輪流場不均勻性加劇,也驗(yàn)證了上述葉輪根部出現(xiàn)相對氣流角增大的現(xiàn)象。同時(shí)由于副葉輪的壓力面負(fù)峰值區(qū)域減少,葉輪的壓縮能力下降。

    圖14 近堵塞點(diǎn)雙面常規(guī)葉輪BVFz分布Fig.14 BVFz distribution of double-sided conventional impeller near the plug point

    綜上可知,將邊界渦量動(dòng)力學(xué)與CFD結(jié)合可以更加直接準(zhǔn)確地找到副葉輪流動(dòng)損失的位置,直管型后進(jìn)氣管的副葉輪損失主要位于葉輪前緣與后緣葉尖處,而橢圓型后進(jìn)氣管的副葉輪前緣根部還存在分離流動(dòng)的情況;副葉輪表面的壓力梯度主要是由副葉輪進(jìn)口流場條件所決定的,因此控制BVFz正峰值區(qū)域和改善葉輪進(jìn)口流場條件,是優(yōu)化后進(jìn)氣管的主要研究方向。

    4 結(jié)論

    (1)相比直管進(jìn)氣,不同幾何形狀后進(jìn)氣管使雙面復(fù)合葉輪的工作穩(wěn)定范圍達(dá)到理想值,但由于受到后進(jìn)氣管換熱和結(jié)構(gòu)的影響,在小流量范圍內(nèi)性能差異最大,而橢圓型后進(jìn)氣管與對稱翼型后進(jìn)氣管的雙面復(fù)合葉輪之間性能無明顯差異。相對流量變化的同時(shí)也改變了雙面復(fù)合葉輪的工作模式,通常副葉輪較主葉輪提前出現(xiàn)喘振。

    (2)受葉輪輪轂滯止及抽吸作用影響,對應(yīng)的后進(jìn)氣管流場出口處靜壓得到提高,且離副葉輪進(jìn)口越近抽吸效果越明顯;橢圓型后進(jìn)氣管因外通道流場均勻,可為進(jìn)氣系統(tǒng)設(shè)計(jì)提供參考。

    (3)直管進(jìn)氣的葉輪進(jìn)口流場分布較為均勻且葉尖區(qū)域呈正弦曲線規(guī)律性強(qiáng),而橢圓型后進(jìn)氣管結(jié)構(gòu)造成的總溫總壓組合畸變與空心渦輪葉片耦合共同決定了葉輪進(jìn)口流場的結(jié)構(gòu),同時(shí)也使葉輪進(jìn)口葉根處的相對氣流角增大,是葉輪性能損失的主要原因之一。

    (4)將邊界渦量動(dòng)力學(xué)與CFD融合,可更方便直觀地找到葉輪流動(dòng)損失的根源,橢圓型后進(jìn)氣管的副葉輪前緣的葉尖、葉根部均有分離流動(dòng)的情況,這也是雙面復(fù)合葉輪流動(dòng)損失的主要原因,有待后期對后進(jìn)氣管進(jìn)一步優(yōu)化。

    猜你喜歡
    總溫橢圓型進(jìn)氣管
    一類帶臨界指數(shù)增長的橢圓型方程組兩個(gè)正解的存在性
    3m×2m結(jié)冰風(fēng)洞熱流場品質(zhì)提高及評估
    9E燃機(jī)二級燃料切換噴嘴進(jìn)氣管彎頭裂紋分析及防范措施
    能源工程(2020年5期)2021-01-04 01:29:00
    一種蓄熱式全氧燃燒窯爐
    玻璃(2020年2期)2020-12-11 22:20:04
    大型結(jié)冰風(fēng)洞熱流場符合性驗(yàn)證
    一種汽車發(fā)動(dòng)機(jī)艙與乘客艙隔熱裝置的設(shè)計(jì)
    高超音速氣流總溫傳感器設(shè)計(jì)
    一類擬線性橢圓型方程的正解
    一類完全非線性橢圓型方程組解的對稱性
    RN擬線性橢圓型方程兩個(gè)非負(fù)解的存在性
    色播在线永久视频| 熟女少妇亚洲综合色aaa.| 欧美变态另类bdsm刘玥| 搡老岳熟女国产| 一个人免费看片子| kizo精华| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美在线一区| 欧美日韩国产mv在线观看视频| 99九九在线精品视频| 久久国产亚洲av麻豆专区| 国产精品欧美亚洲77777| 久久久久精品国产欧美久久久 | 亚洲人成77777在线视频| 久久中文字幕一级| 亚洲黑人精品在线| 日本猛色少妇xxxxx猛交久久| 又黄又粗又硬又大视频| 国产在线观看jvid| 国产免费福利视频在线观看| 国产熟女午夜一区二区三区| 久久国产亚洲av麻豆专区| 欧美成人精品欧美一级黄| a级片在线免费高清观看视频| 亚洲黑人精品在线| 国产视频首页在线观看| 久久久久国产精品人妻一区二区| 亚洲视频免费观看视频| 中文欧美无线码| 成人影院久久| 激情视频va一区二区三区| 少妇被粗大的猛进出69影院| 只有这里有精品99| 激情五月婷婷亚洲| 久久久精品区二区三区| 日本黄色日本黄色录像| 精品一区二区三区四区五区乱码 | 国产不卡av网站在线观看| 国产女主播在线喷水免费视频网站| 一本色道久久久久久精品综合| 午夜91福利影院| 欧美+亚洲+日韩+国产| 老汉色av国产亚洲站长工具| 视频区图区小说| 日日摸夜夜添夜夜爱| 三上悠亚av全集在线观看| 91精品三级在线观看| 国产精品久久久久久人妻精品电影 | 波野结衣二区三区在线| 51午夜福利影视在线观看| 人成视频在线观看免费观看| 久久精品成人免费网站| 久久免费观看电影| 新久久久久国产一级毛片| 少妇粗大呻吟视频| 欧美日韩亚洲高清精品| 国产野战对白在线观看| 国产精品久久久av美女十八| 国产午夜精品一二区理论片| 国产淫语在线视频| 国产老妇伦熟女老妇高清| 中文字幕人妻熟女乱码| 久久精品亚洲av国产电影网| 操美女的视频在线观看| 国产精品三级大全| 天堂8中文在线网| 亚洲成人手机| www.自偷自拍.com| 免费观看人在逋| 首页视频小说图片口味搜索 | 国产在线免费精品| 久久这里只有精品19| 精品少妇内射三级| 新久久久久国产一级毛片| 香蕉丝袜av| 亚洲免费av在线视频| 美国免费a级毛片| 欧美av亚洲av综合av国产av| 亚洲av日韩在线播放| 午夜激情久久久久久久| 亚洲欧美精品综合一区二区三区| 蜜桃在线观看..| 欧美日韩成人在线一区二区| 我的亚洲天堂| 亚洲专区中文字幕在线| 欧美激情高清一区二区三区| 亚洲欧洲日产国产| 欧美日韩国产mv在线观看视频| 精品少妇一区二区三区视频日本电影| 久久久久久免费高清国产稀缺| 韩国精品一区二区三区| 国产一区有黄有色的免费视频| 欧美激情 高清一区二区三区| 女人久久www免费人成看片| 亚洲激情五月婷婷啪啪| 精品一区二区三区四区五区乱码 | 少妇 在线观看| 每晚都被弄得嗷嗷叫到高潮| 午夜免费男女啪啪视频观看| 国产免费福利视频在线观看| 精品一区二区三区四区五区乱码 | 亚洲伊人色综图| www.av在线官网国产| 老汉色av国产亚洲站长工具| 免费高清在线观看日韩| 午夜av观看不卡| 99热国产这里只有精品6| 国产一区二区在线观看av| 中文字幕色久视频| 极品人妻少妇av视频| 中文字幕最新亚洲高清| 欧美国产精品va在线观看不卡| 悠悠久久av| 亚洲精品久久成人aⅴ小说| 狠狠婷婷综合久久久久久88av| 大香蕉久久成人网| 精品久久蜜臀av无| 丁香六月欧美| 婷婷色av中文字幕| 亚洲久久久国产精品| 久久精品亚洲av国产电影网| av又黄又爽大尺度在线免费看| 国语对白做爰xxxⅹ性视频网站| 亚洲伊人色综图| 日本黄色日本黄色录像| 国产精品偷伦视频观看了| 一边亲一边摸免费视频| av不卡在线播放| 国产在线一区二区三区精| 亚洲国产日韩一区二区| 国产成人精品久久二区二区91| 亚洲免费av在线视频| 国精品久久久久久国模美| 亚洲精品在线美女| 中文字幕人妻丝袜一区二区| 如日韩欧美国产精品一区二区三区| 男女国产视频网站| 亚洲欧美一区二区三区国产| 在线观看人妻少妇| 久久精品国产亚洲av涩爱| 人成视频在线观看免费观看| www.精华液| a级毛片黄视频| 国产亚洲av片在线观看秒播厂| 黄色视频不卡| 在线观看免费日韩欧美大片| 夫妻午夜视频| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟女毛片儿| 婷婷色综合www| 妹子高潮喷水视频| 亚洲视频免费观看视频| 1024视频免费在线观看| 一级片'在线观看视频| 王馨瑶露胸无遮挡在线观看| 亚洲图色成人| 成人亚洲欧美一区二区av| 日本wwww免费看| 在线观看免费高清a一片| 下体分泌物呈黄色| 国产又爽黄色视频| 高清黄色对白视频在线免费看| e午夜精品久久久久久久| 欧美少妇被猛烈插入视频| 老司机影院成人| 精品亚洲成国产av| 天天添夜夜摸| 观看av在线不卡| www.熟女人妻精品国产| 丝袜脚勾引网站| 无遮挡黄片免费观看| 夫妻性生交免费视频一级片| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性xxxx| av国产精品久久久久影院| 国产老妇伦熟女老妇高清| 人人妻,人人澡人人爽秒播 | 欧美在线一区亚洲| 91国产中文字幕| 日本a在线网址| 老司机影院毛片| 亚洲美女黄色视频免费看| 午夜激情久久久久久久| 亚洲av片天天在线观看| 国产爽快片一区二区三区| 欧美日本中文国产一区发布| 久9热在线精品视频| 青青草视频在线视频观看| 欧美日韩亚洲高清精品| 久久av网站| 大香蕉久久成人网| 激情五月婷婷亚洲| 高潮久久久久久久久久久不卡| 下体分泌物呈黄色| 午夜福利乱码中文字幕| 人成视频在线观看免费观看| 99国产综合亚洲精品| 日本黄色日本黄色录像| 日韩av不卡免费在线播放| 亚洲伊人久久精品综合| 久久人人爽av亚洲精品天堂| 亚洲一区中文字幕在线| 91成人精品电影| 成人免费观看视频高清| 操出白浆在线播放| 久久国产亚洲av麻豆专区| 一边亲一边摸免费视频| 伊人亚洲综合成人网| 女人爽到高潮嗷嗷叫在线视频| 一区二区三区激情视频| 超碰成人久久| 91精品伊人久久大香线蕉| 人人妻人人添人人爽欧美一区卜| 男女边吃奶边做爰视频| 欧美国产精品一级二级三级| 久久久久久亚洲精品国产蜜桃av| 十分钟在线观看高清视频www| 日本av手机在线免费观看| 亚洲精品第二区| 老司机亚洲免费影院| 一边亲一边摸免费视频| 一本大道久久a久久精品| 亚洲精品在线美女| 青春草视频在线免费观看| 中国美女看黄片| 一级毛片电影观看| 一区二区三区四区激情视频| 中文字幕亚洲精品专区| 视频区欧美日本亚洲| 国产国语露脸激情在线看| 国产亚洲午夜精品一区二区久久| 久久久久国产精品人妻一区二区| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线观看99| 99久久精品国产亚洲精品| 亚洲国产最新在线播放| 国产野战对白在线观看| 国产熟女欧美一区二区| 美女扒开内裤让男人捅视频| 亚洲欧美中文字幕日韩二区| 少妇被粗大的猛进出69影院| 大香蕉久久成人网| 好男人视频免费观看在线| 亚洲专区中文字幕在线| 色网站视频免费| 久久精品久久久久久噜噜老黄| 亚洲自偷自拍图片 自拍| 日本91视频免费播放| 看十八女毛片水多多多| 亚洲综合色网址| 久久国产精品影院| 免费在线观看完整版高清| 国产麻豆69| av不卡在线播放| 日韩 欧美 亚洲 中文字幕| 亚洲国产精品一区二区三区在线| 一级片'在线观看视频| 大片免费播放器 马上看| 波多野结衣av一区二区av| 香蕉丝袜av| 桃花免费在线播放| netflix在线观看网站| 成年人黄色毛片网站| 性少妇av在线| 两个人看的免费小视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲图色成人| 国产一区二区激情短视频 | 国产精品秋霞免费鲁丝片| 在现免费观看毛片| 国产免费又黄又爽又色| 久久性视频一级片| 99国产精品99久久久久| 国产高清视频在线播放一区 | 亚洲成色77777| 国产一区二区三区av在线| 91精品国产国语对白视频| videos熟女内射| av网站免费在线观看视频| 亚洲国产精品一区三区| 久久久久久人人人人人| 90打野战视频偷拍视频| 在线av久久热| 色婷婷久久久亚洲欧美| 国产xxxxx性猛交| 丝袜人妻中文字幕| 婷婷丁香在线五月| 国产精品国产三级国产专区5o| 成人免费观看视频高清| 最新在线观看一区二区三区 | 亚洲,一卡二卡三卡| netflix在线观看网站| 麻豆乱淫一区二区| 国产亚洲精品第一综合不卡| 中文乱码字字幕精品一区二区三区| 免费av中文字幕在线| 久久国产亚洲av麻豆专区| 久久久久国产精品人妻一区二区| 亚洲av综合色区一区| 国产精品久久久久久精品古装| 国产精品香港三级国产av潘金莲 | 一区在线观看完整版| 久久久国产一区二区| 黄片播放在线免费| 五月天丁香电影| 亚洲国产精品国产精品| 日韩欧美一区视频在线观看| 亚洲国产欧美一区二区综合| www.熟女人妻精品国产| 久久精品久久久久久久性| 精品一区二区三卡| 久久精品人人爽人人爽视色| 国产成人一区二区三区免费视频网站 | 一个人免费看片子| 波野结衣二区三区在线| 国产精品一国产av| 午夜日韩欧美国产| 久久久国产欧美日韩av| 国产高清不卡午夜福利| 欧美黄色淫秽网站| 亚洲国产欧美网| 精品亚洲乱码少妇综合久久| 大陆偷拍与自拍| 99热国产这里只有精品6| 成人18禁高潮啪啪吃奶动态图| 在线观看www视频免费| 国产精品av久久久久免费| 国产激情久久老熟女| 亚洲中文字幕日韩| 一二三四社区在线视频社区8| 精品久久久久久电影网| 我要看黄色一级片免费的| 亚洲国产看品久久| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 国产精品麻豆人妻色哟哟久久| 欧美精品一区二区大全| 天堂8中文在线网| 亚洲 国产 在线| 性色av乱码一区二区三区2| 免费av中文字幕在线| 黄色 视频免费看| 国产三级黄色录像| 亚洲伊人久久精品综合| 亚洲一区中文字幕在线| 亚洲精品美女久久av网站| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 叶爱在线成人免费视频播放| 成人三级做爰电影| 蜜桃国产av成人99| 你懂的网址亚洲精品在线观看| 丰满迷人的少妇在线观看| 97在线人人人人妻| 欧美精品一区二区免费开放| 王馨瑶露胸无遮挡在线观看| h视频一区二区三区| 波多野结衣一区麻豆| 一级毛片我不卡| 亚洲国产看品久久| 波多野结衣一区麻豆| 高清欧美精品videossex| 久久毛片免费看一区二区三区| 亚洲av成人精品一二三区| 日本欧美国产在线视频| 久久精品国产亚洲av涩爱| 亚洲成国产人片在线观看| 亚洲成人免费电影在线观看 | 国产成人啪精品午夜网站| 美女视频免费永久观看网站| 国产成人精品无人区| 一区二区三区四区激情视频| 日本欧美视频一区| 男人操女人黄网站| 国产精品熟女久久久久浪| 欧美黄色片欧美黄色片| 99热国产这里只有精品6| 亚洲三区欧美一区| 成在线人永久免费视频| 中国国产av一级| 精品一品国产午夜福利视频| 黄色片一级片一级黄色片| 秋霞在线观看毛片| 看免费成人av毛片| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 亚洲欧美色中文字幕在线| 亚洲美女黄色视频免费看| 男男h啪啪无遮挡| 久热这里只有精品99| 最黄视频免费看| 国产成人a∨麻豆精品| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看| 在线观看国产h片| 国产高清视频在线播放一区 | 十八禁网站网址无遮挡| 欧美黄色淫秽网站| a 毛片基地| 久久久欧美国产精品| 波野结衣二区三区在线| 波多野结衣一区麻豆| 中文欧美无线码| 久久久国产一区二区| 国产熟女午夜一区二区三区| 欧美大码av| 亚洲中文日韩欧美视频| 女人爽到高潮嗷嗷叫在线视频| 国产精品99久久99久久久不卡| 一边摸一边做爽爽视频免费| 久久国产精品男人的天堂亚洲| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 在线观看人妻少妇| 成人国产一区最新在线观看 | 免费观看av网站的网址| 亚洲七黄色美女视频| 久久九九热精品免费| 一级毛片电影观看| 欧美 亚洲 国产 日韩一| 午夜影院在线不卡| 午夜福利在线免费观看网站| 欧美在线一区亚洲| √禁漫天堂资源中文www| 人人妻人人添人人爽欧美一区卜| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 国产成人精品久久二区二区91| 一级毛片 在线播放| 欧美国产精品va在线观看不卡| 午夜两性在线视频| 亚洲第一av免费看| 日本a在线网址| 亚洲欧洲日产国产| 亚洲欧美中文字幕日韩二区| 国产免费又黄又爽又色| 久久中文字幕一级| 人人妻人人澡人人看| 国产一区二区 视频在线| 丝袜脚勾引网站| 国产xxxxx性猛交| 婷婷色麻豆天堂久久| 亚洲情色 制服丝袜| 成人国产一区最新在线观看 | 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久人妻精品电影 | 国产成人影院久久av| 国产欧美日韩综合在线一区二区| 国产视频首页在线观看| 国产精品熟女久久久久浪| 精品亚洲成国产av| 老司机影院毛片| 欧美国产精品一级二级三级| 悠悠久久av| 人体艺术视频欧美日本| 亚洲国产精品一区三区| 国产精品久久久久久精品古装| 性高湖久久久久久久久免费观看| 老司机午夜十八禁免费视频| 日本欧美国产在线视频| 久久狼人影院| 国产成人精品久久二区二区91| 亚洲国产中文字幕在线视频| 亚洲五月色婷婷综合| 人妻一区二区av| 精品一区二区三区四区五区乱码 | 悠悠久久av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产爽快片一区二区三区| 婷婷色综合大香蕉| 激情视频va一区二区三区| 99国产综合亚洲精品| 亚洲成人手机| 人人妻人人澡人人爽人人夜夜| 啦啦啦视频在线资源免费观看| 国产av一区二区精品久久| 可以免费在线观看a视频的电影网站| 欧美+亚洲+日韩+国产| 亚洲av成人精品一二三区| 在线观看www视频免费| 久久精品亚洲av国产电影网| 国产精品99久久99久久久不卡| 高清黄色对白视频在线免费看| 亚洲精品av麻豆狂野| 日韩熟女老妇一区二区性免费视频| 日日夜夜操网爽| 精品少妇黑人巨大在线播放| 中国国产av一级| 成人午夜精彩视频在线观看| 777米奇影视久久| 在线观看人妻少妇| 成人手机av| 秋霞在线观看毛片| 最新的欧美精品一区二区| 水蜜桃什么品种好| 亚洲黑人精品在线| av电影中文网址| 日韩熟女老妇一区二区性免费视频| 精品少妇久久久久久888优播| 亚洲第一青青草原| 国产一区二区 视频在线| 高清视频免费观看一区二区| 国产成人免费观看mmmm| www.自偷自拍.com| 国产精品久久久久久人妻精品电影 | 黄色怎么调成土黄色| 欧美+亚洲+日韩+国产| 国产精品偷伦视频观看了| 亚洲精品日本国产第一区| 日韩熟女老妇一区二区性免费视频| 看免费av毛片| 亚洲,欧美精品.| 婷婷丁香在线五月| 999久久久国产精品视频| 如日韩欧美国产精品一区二区三区| 婷婷色综合www| 波野结衣二区三区在线| 精品高清国产在线一区| 国产欧美亚洲国产| 亚洲国产精品一区三区| 欧美变态另类bdsm刘玥| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久大尺度免费视频| 国产精品成人在线| 久久国产精品影院| 麻豆av在线久日| 97人妻天天添夜夜摸| 亚洲av欧美aⅴ国产| 欧美 亚洲 国产 日韩一| 性色av乱码一区二区三区2| 精品高清国产在线一区| 19禁男女啪啪无遮挡网站| 亚洲精品av麻豆狂野| 成年人免费黄色播放视频| 热re99久久国产66热| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 欧美在线一区亚洲| 亚洲欧美精品综合一区二区三区| 一二三四社区在线视频社区8| 精品高清国产在线一区| 亚洲激情五月婷婷啪啪| 亚洲少妇的诱惑av| 国产一区二区 视频在线| 中文字幕人妻丝袜制服| 丝袜在线中文字幕| 精品欧美一区二区三区在线| 在线天堂中文资源库| 手机成人av网站| 别揉我奶头~嗯~啊~动态视频 | 亚洲,一卡二卡三卡| 久久天堂一区二区三区四区| 各种免费的搞黄视频| 亚洲国产看品久久| 亚洲精品国产一区二区精华液| 亚洲国产精品国产精品| 精品久久久久久电影网| 亚洲av在线观看美女高潮| 免费在线观看黄色视频的| 另类精品久久| 亚洲人成电影观看| 久久精品亚洲熟妇少妇任你| 99国产精品一区二区蜜桃av | 一区二区三区激情视频| 欧美精品人与动牲交sv欧美| 久久久久网色| 色精品久久人妻99蜜桃| 国产亚洲一区二区精品| 亚洲天堂av无毛| 欧美日韩黄片免| 国产精品成人在线| 国产老妇伦熟女老妇高清| 欧美在线一区亚洲| 手机成人av网站| 久久精品成人免费网站| 亚洲熟女毛片儿| 亚洲欧美成人综合另类久久久| 精品少妇一区二区三区视频日本电影| 视频在线观看一区二区三区| 99国产精品99久久久久| 亚洲人成77777在线视频| 中文字幕高清在线视频| 波野结衣二区三区在线| 亚洲精品av麻豆狂野| 国产欧美亚洲国产| 99久久综合免费| 成在线人永久免费视频| 日本欧美国产在线视频| 九草在线视频观看| 男女国产视频网站| 免费在线观看完整版高清| 黄色怎么调成土黄色| 国产精品.久久久| 日日摸夜夜添夜夜爱| 五月开心婷婷网| 最近手机中文字幕大全| 一级片免费观看大全| 美女视频免费永久观看网站| 老司机午夜十八禁免费视频| 男人添女人高潮全过程视频| 两个人看的免费小视频| 国产精品一区二区在线观看99| 婷婷色综合www| 久久天躁狠狠躁夜夜2o2o | 亚洲一区二区三区欧美精品| 妹子高潮喷水视频| 日本一区二区免费在线视频| 国产精品久久久av美女十八| 久久ye,这里只有精品| 女人久久www免费人成看片| 热99国产精品久久久久久7| 只有这里有精品99| 91字幕亚洲|