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

    基于新型智能探頭的大氣數(shù)據(jù)系統(tǒng)算法研究

    2022-09-23 01:50:10駿,鄭偉,王
    測控技術(shù) 2022年9期
    關(guān)鍵詞:馬赫數(shù)迎角靜壓

    劉 駿,鄭 偉,王 禹

    (航空工業(yè)沈陽飛機(jī)設(shè)計(jì)研究所,遼寧 沈陽 110035)

    大氣數(shù)據(jù)系統(tǒng)是一種具有綜合性和高精度特點(diǎn)的大氣數(shù)據(jù)測量解算系統(tǒng)[1-3],可以提供高度Hp、空速Vc、馬赫數(shù)Ma、迎角αt和側(cè)滑角βt等飛行器飛行所需的重要飛行參數(shù),對于保障飛行器的飛行安全起著至關(guān)重要的作用。隨著我國科學(xué)技術(shù)的不斷發(fā)展進(jìn)步,相關(guān)的技術(shù)工藝水平也得到了進(jìn)一步的提高[4],大氣數(shù)據(jù)系統(tǒng)已由最早期的多個(gè)分立式傳感器組合,發(fā)展到由多個(gè)智能探頭組成的分布式大氣數(shù)據(jù)系統(tǒng)[5]。

    在基于智能探頭的分布式大氣數(shù)據(jù)系統(tǒng)中,每一個(gè)或每一組智能探頭都相當(dāng)于一套獨(dú)立的大氣數(shù)據(jù)測量輸出系統(tǒng),可以根據(jù)飛機(jī)測量余度的需要,在機(jī)身不同位置處分別安裝數(shù)只,組成一套分布式大氣數(shù)據(jù)系統(tǒng)。根據(jù)測量形式的不同,目前常用的智能探頭可分為3種[6],分別是風(fēng)標(biāo)式智能探頭、空速管式智能探頭和錐形智能探頭,如圖1所示。

    本文主要設(shè)計(jì)了一種基于半圓柱形智能探頭(下文簡稱為半圓柱探頭)的分布式大氣數(shù)據(jù)系統(tǒng)的算法模型(該種半圓柱探頭由測壓傳感器、智能電路和數(shù)據(jù)接口等部分組成),并開展了基于此算法模型的仿真驗(yàn)證工作,為后續(xù)大氣數(shù)據(jù)系統(tǒng)的設(shè)計(jì)研究工作提供一種新的思路。

    1 半圓柱探頭外形及安裝信息

    該種半圓柱探頭整體呈半圓柱形,圓柱端為安裝端,自由端為半圓柱型,外側(cè)靠近自由端處開4個(gè)測壓孔,其中一個(gè)測壓孔在半圓柱體臺(tái)階平面上,另外3個(gè)測壓孔均勻分布在背部曲面上,半圓柱探頭外形示意圖如圖2所示。

    圖2 半圓柱探頭布局示意圖

    半圓柱探頭在機(jī)上安裝后,通過CFD(Computational Fluid Dynamics,計(jì)算流體力學(xué))軟件計(jì)算得到機(jī)身附近局部流場信息,選擇半圓柱探頭安裝位置[7-8],將每兩個(gè)半圓柱探頭為一組,垂直蒙皮對稱安裝在飛行器前機(jī)身左下和右下腹部處,探頭圓柱面正對來流方向,臺(tái)階平面處于背風(fēng)面。

    2 半圓柱探頭的大氣數(shù)據(jù)系統(tǒng)算法建模

    2.1 算法模型基本原理

    分布式大氣數(shù)據(jù)系統(tǒng)測量大氣數(shù)據(jù)的基本原理[9-10]是:通過風(fēng)洞試驗(yàn)和CFD計(jì)算,標(biāo)定得到在不同狀態(tài)下智能探頭感受到的局部來流壓力、角度和飛機(jī)遠(yuǎn)端自由來流的壓力、角度之間的定量關(guān)系,之后再由大氣解算部件利用局部壓力、角度計(jì)算得到表征自由流的大氣參數(shù)。按照基于半圓柱探頭的大氣數(shù)據(jù)系統(tǒng)方案構(gòu)型[11],系統(tǒng)應(yīng)包括對稱安裝于前機(jī)身左右表面的4個(gè)半圓柱探頭(兩兩對稱),每個(gè)探頭(橫截面外形如圖3所示)感受并測量輸出4個(gè)壓力給大氣解算部件。

    圖3 半圓柱探頭橫截面示意圖

    具體來說,機(jī)頭前遠(yuǎn)方自由來流Pt(QC)、PS、Ma、αt、βt受機(jī)身擾動(dòng),氣流流經(jīng)左、右探頭所在位置處時(shí),其相對左、右探頭的局部來流分別為PLt(QLC)、PLS、MaL、θL和PRt(QRC)、PRS、MaR、θR。左、右局部來流再流經(jīng)左、右探頭時(shí),左、右探頭上各4個(gè)測壓孔感受的壓力分別為PL1、PL2、PL3、PL4和PR1、PR2、PR3、PR4。

    大氣數(shù)據(jù)系統(tǒng)算法軟件模塊駐留在大氣解算部件的每個(gè)通道中,以其中一個(gè)通道為例,其接收對稱位置兩個(gè)探頭的8個(gè)壓力信號為PL1、PL2、PL3、PL4、PR1、PR2、PR3、PR4,先基于半圓柱體探頭本體的氣動(dòng)模型,利用左、右半圓柱探頭各自4個(gè)點(diǎn)壓力解算得到各自所在位置機(jī)身局部流場的QLC、PLS、MaL、θL和QRC、PRS、MaR、θR。

    再根據(jù)探頭所在位置局部流場的氣動(dòng)模型,利用左、右探頭局部來流角均值求解真迎角αt,利用左、右探頭局部來流角差值求解真?zhèn)然铅聇,取左、右探頭局部流場動(dòng)壓均值為指示動(dòng)壓QCi,取左、右探頭局部流場靜壓均值為指示靜壓PSi,利用QCi和PSi之比求解指示Mai。再以αt、βt和Mai為基礎(chǔ),根據(jù)CFD計(jì)算獲取的半圓柱探頭所在位置局部流場的氣動(dòng)模型,完成由指示值到真實(shí)值的迭代計(jì)算,得到當(dāng)前飛行狀態(tài)下的QC、PS、Ma、αt和βt。

    2.2 算法模型相關(guān)參數(shù)定義

    本節(jié)給出氣動(dòng)建模所需部分樣本數(shù)據(jù)涉及到的各參數(shù)定義。

    ① 左、右探頭上各測壓點(diǎn)的壓力系數(shù)定義:

    CpLi(MaL,θL)=(PLi-PLS)/QLC,i=1,2,3,4

    (1)

    CpRi(MaR,θR)=(PRi-PRS)/QRC,i=1,2,3,4

    (2)

    ② 左、右探頭位置處機(jī)身局部流場靜壓誤差壓力系數(shù):

    CpLS(Ma,αt,βt)=(PLS-PS)/QC

    (3)

    CpRS(Ma,αt,βt)=(PRS-PS)/QC

    (4)

    ③ 左、右探頭位置處機(jī)身局部流場總壓誤差壓力系數(shù):

    CpLt(Ma,αt,βt)=(PLt-Pt)/QC

    (5)

    CpRt(Ma,αt,βt)=(PRt-Pt)/QC

    (6)

    ④ 左、右機(jī)身局部流場來流角與馬赫數(shù)Ma、真迎角αt、真?zhèn)然铅聇關(guān)系:

    θL=θL(Ma,αt,βt)

    (7)

    θR=θR(Ma,αt,βt)

    (8)

    2.3 氣動(dòng)模型設(shè)計(jì)與分析

    從2.1節(jié)的內(nèi)容中可知,大氣數(shù)據(jù)系統(tǒng)算法建模的核心分為兩方面,一是獲取半圓柱探頭本體4個(gè)測壓點(diǎn)壓力系數(shù)樣本數(shù)據(jù),并建立相應(yīng)測壓點(diǎn)間壓力系數(shù)差與局部來流動(dòng)、靜壓和角度之間關(guān)系;二是獲取半圓柱探頭所在位置處局部流場的氣動(dòng)模型樣本數(shù)據(jù),并建立局部來流和自由來流之間的映射關(guān)系[12-13]。在確定探頭外形和探頭在機(jī)上的安裝位置和安裝角度之后,上述壓力系數(shù)樣本數(shù)據(jù)即可通過風(fēng)洞試驗(yàn)和CFD計(jì)算獲取。

    2.3.1 半圓柱探頭本體氣動(dòng)模型

    以飛機(jī)左側(cè)探頭為例,首先根據(jù)CFD計(jì)算數(shù)據(jù),做出探頭上P1、P2、P3和P4四個(gè)測壓點(diǎn)處不同馬赫數(shù)下壓力系數(shù)隨局部來流角變化的曲線,圖像如圖4和圖5所示。

    圖4 0.4馬赫數(shù)下各測壓點(diǎn)處壓力系數(shù)隨局部迎角變化曲線

    圖5 0.8馬赫數(shù)下各測壓點(diǎn)壓力系數(shù)隨局部迎角變化曲線

    通過對圖4和圖5中的曲線進(jìn)行分析,可以發(fā)現(xiàn)在不同馬赫數(shù)下,P2點(diǎn)處壓力系數(shù)對局部流場來流角變化不太敏感,在不同馬赫數(shù)下數(shù)值雖隨之變化但基本能保持穩(wěn)定,結(jié)合圖3所示流場來流和探頭處位置關(guān)系示意圖,可以確定用P2測壓點(diǎn)處測量得到的壓力值可以作為探頭所在處的局部靜壓的指示靜壓使用。同時(shí),P1、P3和P4三個(gè)測壓點(diǎn)處的壓力系數(shù)在不同馬赫數(shù)下的隨局部來流角變化時(shí),總有一個(gè)點(diǎn)的壓力系數(shù)接近1,說明此時(shí)該點(diǎn)近似于正對來流方向,該測壓點(diǎn)感受到的壓力值可以作為探頭所在處的局部總壓的指示總壓使用,所以可得到局部流場指示靜壓:

    PLSi=PL2

    (9)

    局部流場指示動(dòng)壓:

    QLCi=PLmax-PLSi

    (10)

    式中:PLmax為PL1、PL3、PL4中的最大值。

    局部流場指示馬赫數(shù):

    MaLi=f(QLCi/PLSi)

    (11)

    式中:MaLi按HB6127—1986第3.9節(jié)公式依據(jù)動(dòng)靜壓比值求解(當(dāng)Ma<1時(shí),PS/PT=(1+0.2Ma2)-3.5)。

    在單獨(dú)一個(gè)測壓點(diǎn)的壓力系數(shù)不具備明顯的隨來流角變化的規(guī)律的情況下,接著研究不同馬赫數(shù)下3個(gè)測壓點(diǎn)間壓力系數(shù)的差值隨來流角變化的規(guī)律,結(jié)果如圖6、圖7所示。

    圖6 Ma=0.4時(shí)各測壓點(diǎn)間壓力系數(shù)差

    圖7 Ma=0.8時(shí)各測壓點(diǎn)間壓力系數(shù)差

    觀察圖6和圖7內(nèi)各測壓點(diǎn)間壓力系數(shù)差隨來流角變化的曲線,可以發(fā)現(xiàn)CpL3-L1在隨來流角變化的過程中呈一條曲線,不具備明顯的規(guī)律性;而CpL4-L1和CpL4-L3在對應(yīng)的角度范圍之內(nèi),即CpL4-L1在來流角處于0°~60°之間時(shí)以及CpL4-L3在來流角處于-30°~30°之間時(shí),與來流角呈近似線性相關(guān)的關(guān)系。所以根據(jù)上述描述可以按照式(12)、式(13)完成探頭本體氣動(dòng)建模設(shè)計(jì)和分析,處理得到左、右探頭兩測壓孔間壓力系數(shù)之差與局部來流角θ的關(guān)系曲線。

    ΔCpL4-Lj(MaL,θL)=CpL4(MaL,θL)-CpLj(MaL,θL)

    =(PL4-PLj)/QLC=KθL+θ0

    (12)

    ΔCpR4-Rj(MaR,θR)=CpR4(MaR,θR)-CpRj(MaR,θR)

    =(PR4-PRj)/QRC=KθR+θ0

    (13)

    式中:當(dāng)θ≤30°時(shí),j=3;當(dāng)θ>30°時(shí),j=1;K和θ0按氣動(dòng)樣本數(shù)據(jù)處理得到,隨Ma不同取具體值。

    但在后續(xù)算法建模中使用上述關(guān)系式時(shí)可以發(fā)現(xiàn),本節(jié)建立的半圓柱體探頭本體的氣動(dòng)模型,是建立探頭上4個(gè)測壓點(diǎn)間的壓力系數(shù)差與半圓柱探頭所在安裝位置處機(jī)身局部來流角的之間的映射關(guān)系,即在這個(gè)模型中,局部來流角作為關(guān)系式中的因變量不能作為上式中j取值的判斷依據(jù),應(yīng)采取Cp1、Cp2、Cp3、Cp4或者CpL4-L3和CpL4-L1作為j取值的判斷依據(jù),由于不同馬赫數(shù)下,測壓點(diǎn)間壓力系數(shù)差變化規(guī)律極為類似,故以Ma=0.4時(shí)的情況為例進(jìn)行分析,仔細(xì)分析圖6可以看出,在CpL4-L1≥-1時(shí),CpL4-L1與局部來流角θL呈線性關(guān)系;而在CpL4-L1<-1時(shí),CpL4-L3與局部來流角θL呈線性關(guān)系。故綜合考慮之后,將CpL4-L1≥-1作為式(12)和式(13)中j值的取值判斷依據(jù),即在式(12)和式(13)中,當(dāng)CpL4-L1≥-1時(shí),j=1;當(dāng)CpL4-L1<-1時(shí),j=3;K和θ0按氣動(dòng)樣本數(shù)據(jù)處理得到,隨Ma不同取具體值。

    將式(10)帶入到式(12)和式(13)中,則可得到探頭所在處局部來流角指示值的表達(dá)式:

    θLi=[(PL4-PLj)/QLCi-θ0]/K

    (14)

    θRi=[(PR4-PRj)/QRCi-θ0]/K

    (15)

    至此,根據(jù)探頭上4個(gè)測壓點(diǎn)感受到的壓力值即可求解得到探頭所在位置處的局部動(dòng)壓指示QLCi、局部靜壓指示值PLSi、局部馬赫數(shù)指示值MaLi和局部來流角指示值θLi共4個(gè)指示值。接下來再進(jìn)行迭代求解計(jì)算,使各參數(shù)收斂至系統(tǒng)要求的精度,即可求得左側(cè)探頭所在位置處的局部動(dòng)壓QLC、局部靜壓PLS、局部馬赫數(shù)MaL和局部來流角θL。同理可解得右側(cè)半圓柱探頭所在位置處的局部動(dòng)壓QRC、局部靜壓PRS、局部馬赫數(shù)MaR和局部來流角θR。

    2.3.2 探頭所在位置局部流場氣動(dòng)模型

    本節(jié)要完成半圓柱探頭所在位置局部流場氣動(dòng)模型設(shè)計(jì)分析工作。前文已述,按照半圓柱探頭大氣數(shù)據(jù)系統(tǒng)方案構(gòu)型,半圓柱探頭大氣數(shù)據(jù)系統(tǒng)包含對稱安裝于前機(jī)身左右表面的4個(gè)半圓柱體探頭(兩兩對稱),但在本節(jié)進(jìn)行半圓柱探頭所在位置局部流場氣動(dòng)模型時(shí),只選擇一對對稱安裝的半圓柱體探頭作為研究對象即可。

    首先研究探頭所在位置局部流場的總、靜壓誤差壓力系數(shù)模型,左、右靜壓誤差壓力系數(shù)隨迎角變化的曲線如圖8所示。

    圖8中,圖8(a)和圖8(b)是馬赫數(shù)為0.4的情況下,局部流場左右靜壓誤差壓力系數(shù)隨迎角變化的曲線,圖8(c)和圖8(d)則是馬赫數(shù)為0.8的情況下的曲線。可以明顯看出,不同馬赫數(shù)下,其局部流場左、右靜壓誤差壓力系數(shù)的變化趨勢幾乎相同,所以在研究其變化規(guī)律時(shí),只以馬赫數(shù)為0.4的情況為例進(jìn)行分析。而對比圖8(a)和圖8(b)兩圖可以發(fā)現(xiàn),在側(cè)滑角為5°的時(shí)候,局部流場左靜壓誤差壓力系數(shù)的變化曲線和側(cè)滑角為-5°時(shí)的局部流場右靜壓誤差壓力系數(shù)變化曲線幾乎重合;類似地,側(cè)滑角為-5°的時(shí)候,局部流場左靜壓誤差壓力系數(shù)的變化曲線和側(cè)滑角為5°時(shí)的局部流場右靜壓誤差壓力系數(shù)變化曲線也幾乎重合;而側(cè)滑角為0°時(shí)的變化曲線則處于側(cè)滑角為5°和-5°時(shí)變化曲線中間,此時(shí)考慮取局部流場左、右靜壓誤差壓力系數(shù)的均值,做出圖像如圖9所示。

    圖8 不同馬赫數(shù)下左右靜壓誤差壓力系數(shù)隨迎角變化曲線

    圖9 局部流場左右靜壓誤差壓力系數(shù)均值隨迎角變化曲線

    觀察圖9可以發(fā)現(xiàn),此時(shí)3條曲線幾乎重合,而且都貼近側(cè)滑角為0°時(shí)的曲線,說明局部流場左、右靜壓誤差壓力系數(shù)的均值與側(cè)滑角的變化沒有較大的關(guān)系,只與馬赫數(shù)和迎角變化有關(guān)。類似地,局部流場左、右總壓誤差壓力系數(shù)的變化規(guī)律也有一樣的特性,如圖10所示。

    圖10 不同馬赫數(shù)下局部流場左右總壓誤差壓力系數(shù)隨迎角變化曲線

    圖11為左右總壓誤差壓力系數(shù)均值隨迎角變化曲線。從圖11中可以看出,局部流場左右總壓誤差壓力系數(shù)均值的變化規(guī)律與局部流場左右靜壓誤差壓力系數(shù)均值的變化規(guī)律類似,局部流場總壓誤差壓力系數(shù)的變化與馬赫數(shù)和迎角變化有關(guān),同時(shí)與側(cè)滑角的絕對值有關(guān)。以探頭所在位置的局部流場總、靜壓誤差壓力系數(shù)樣本數(shù)據(jù)集為輸入,可以按照以下定義和式(16)、式(17)完成探頭所在位置局部流場總、靜壓氣動(dòng)建模設(shè)計(jì)和分析。

    圖11 不同馬赫數(shù)下左右總壓誤差壓力系數(shù)均值隨迎角變化曲線

    (1)取左、右探頭局部流場靜壓均值為指示靜壓PSi,則靜壓誤差壓力系數(shù)為

    CpS(Ma,αt,βt)=(PSi-PS)/QC

    =[CpLS(Ma,αt,βt)+CpRS(Ma,αt,βt)]/2

    (16)

    (2)取左、右探頭局部流場動(dòng)總均值為指示總壓Pti,則總壓誤差壓力系數(shù)為

    Cpt(Ma,αt,βt)=(Pti-Pt)/QC

    =[CpLt(Ma,αt,βt)+CpRt(Ma,αt,βt)]/2

    (17)

    確定總壓誤差壓力系數(shù)的變化規(guī)律和靜壓誤差壓力系數(shù)的變化規(guī)律之后,要進(jìn)一步分析左、右局部來流角θL、θR和迎角、側(cè)滑角及馬赫數(shù)之間的映射關(guān)系。

    通過對CFD流場數(shù)據(jù)進(jìn)行分析可以發(fā)現(xiàn),在同一迎角、側(cè)滑角狀態(tài)下,馬赫數(shù)改變但θL和θR值保持不變,所以可以認(rèn)為θL和θR均只為αt、βt兩個(gè)變量的函數(shù)。此時(shí)以Ma=0.4的情況為例,做出在不同側(cè)滑角狀態(tài)下,θL和θR隨迎角變化的曲線,如圖12、圖13所示。

    圖12 θL隨迎角變化曲線

    圖13 θR隨迎角變化曲線

    對比圖12和圖13中的曲線可以發(fā)現(xiàn)以下幾點(diǎn):

    ① 不同側(cè)滑角情況下的局部來流角隨迎角變化的曲線,其斜率都是相同的;

    ②θL在側(cè)滑角為5°的曲線和θR在側(cè)滑角為-5°的曲線基本重合,類似地,θL在側(cè)滑角為-5°的曲線和θR在側(cè)滑角為5°的曲線基本重合;

    ③ 在同一組曲線中,以θL的變化曲線為例,側(cè)滑角分別為5°和-5°時(shí)的θL曲線,基本以側(cè)滑角為0°時(shí)的曲線為對稱軸,呈對稱關(guān)系。

    基于上述特點(diǎn),可以發(fā)現(xiàn),若是求得同一側(cè)滑角情況下的θL和θR均值,應(yīng)為3條基本重合的曲線,這說明,迎角的變化可以只與θL和θR的均值有關(guān),與馬赫數(shù)和側(cè)滑角均無關(guān),曲線如圖14所示。

    圖14 不同側(cè)滑角情況下θL和θR均值隨迎角變化曲線

    綜上所述,根據(jù)θL值(θR值)對應(yīng)αt、βt兩個(gè)變量的二維樣本數(shù)據(jù)表格,可以按照式(18)~式(20)處理得到相應(yīng)關(guān)系曲線。

    (3)左、右局部來流角與真迎角αt關(guān)系:

    θL(αt,βt)=K1αt+αL0

    (18)

    θR(αt,βt)=K1αt+αR0

    (19)

    (4)左、右局部來流角均值與真迎角αt關(guān)系:

    [θL(αt,βt)+θR(αt,βt)]/2=K1αt+α0

    式中:K1、αL0和αR0按樣本數(shù)據(jù)處理得到,其中K1不隨Ma數(shù)和βt變化,為定值,αL0和αR0隨βt不同取值不同。

    通過上面的分析已經(jīng)驗(yàn)證了θL+θR的均值的變化只與迎角的變化有關(guān),與側(cè)滑角還有馬赫數(shù)無關(guān),接下來進(jìn)一步分析θL和θR的差θL-R與迎角、側(cè)滑角之間的關(guān)系。左右局部來流角差值在不同馬赫數(shù)下隨側(cè)滑角變化的情況如圖15所示。

    圖15 不同迎角下θL-R隨側(cè)滑角變化曲線

    分析圖15可以發(fā)現(xiàn),在不同的馬赫數(shù)和迎角狀態(tài)下,θL-R與側(cè)滑角呈線性關(guān)系,說明θL-R與側(cè)滑角存在一種斜率與迎角有關(guān)的線性關(guān)系,據(jù)此可以得到式(20)。

    (5)左、右局部來流角之差與βt關(guān)系:

    ΔθL-R(αt,βt)=θL(αt,βt)-θR(αt,βt)=K2βt

    (20)

    式中:K2按樣本數(shù)據(jù)處理得到,隨αt和Ma不同取值不同,但是根據(jù)圖15可以看出,K2的取值并不是隨αt以及Ma線性變化,所以最終K2的取值可以隨著αt和Ma的變化進(jìn)行線性插值計(jì)算。

    綜上所述,可以得到如下公式。

    指示動(dòng)壓:

    QCi=(QLC+QRC)/2

    (21)

    指示靜壓:

    PSi=(PLS+PRS)/2

    (22)

    指示總壓:

    Pti=PSi+Qci

    (23)

    指示馬赫數(shù):

    Mai=f(QCi/PSi)

    (24)

    真迎角:

    αt=[(θL+θR)/2-α0]/K1

    (25)

    真?zhèn)然牵?/p>

    βt=(θL-θR)/K2

    (26)

    與2.3.1節(jié)內(nèi)容類似,至此按照上述描述,根據(jù)左右局部來流角、局部總/靜壓,即可求得遠(yuǎn)端自由流場的指示靜壓PSi、指示動(dòng)壓QCi、指示總壓Pti、指示馬赫數(shù)Mai四個(gè)指示值和真迎角αt、真?zhèn)然铅聇。接下來進(jìn)行迭代求解計(jì)算,求得真實(shí)靜壓PS、真實(shí)動(dòng)壓QC、真實(shí)馬赫數(shù)Ma。

    2.3.3 逆算法模型

    由于仿真計(jì)算驗(yàn)證的需要設(shè)計(jì)了逆算法,逆算法是基于同樣的氣動(dòng)建模設(shè)計(jì)輸入條件,由PS、QC、Ma、αt、βt計(jì)算得到左、右探頭8個(gè)測壓點(diǎn)壓力?;舅悸肥峭ㄟ^對Ma、αt、βt進(jìn)行查表和插值計(jì)算,得到左右探頭處局部流場信息,再由局部流場的局部總壓、局部靜壓、局部馬赫數(shù)和局部來流角信息進(jìn)一步通過模型公式和插值計(jì)算得到半圓柱探頭上4個(gè)測壓點(diǎn)處的壓力。

    3 算法仿真驗(yàn)證與分析

    3.1 仿真結(jié)果

    針對算法仿真驗(yàn)證的基本思路是,截取一段真實(shí)的飛參數(shù)據(jù)作為輸入,以動(dòng)壓QC、靜壓PS、馬赫數(shù)Ma、迎角αt、側(cè)滑角βt為輸入值,通過逆向算法計(jì)算得到左右探頭上各4個(gè)測壓孔處壓力值。再將8個(gè)測壓孔處壓力值帶入到半圓柱探頭大氣數(shù)據(jù)系統(tǒng)解算模型中來,得到解算出的解算動(dòng)壓QC′、解算靜壓PS′、解算馬赫數(shù)Ma′、解算迎角αt′、解算側(cè)滑角βt′,對比解算前后對應(yīng)數(shù)據(jù)的差值來驗(yàn)證算法模型的準(zhǔn)確性。

    仿真的輸入數(shù)據(jù)截取了一段某飛機(jī)試飛的飛參數(shù)據(jù),在該段時(shí)間內(nèi),飛機(jī)呈先上升再平飛,先加速再勻速的飛行狀態(tài),迎角和側(cè)滑角均有較大的變化幅度,具體的仿真結(jié)果如圖16~圖20所示。

    圖16 迎角解算結(jié)果

    圖17 側(cè)滑角解算結(jié)果

    圖18 馬赫數(shù)解算結(jié)果

    圖19 靜壓解算結(jié)果

    從圖16~圖20的仿真結(jié)果曲線來看,每組仿真曲線的圖(a)為各大氣參數(shù)的輸入值和解算值的對比圖,圖(b)為各大氣參數(shù)的輸入值和解算值的差值曲線圖,迎角的誤差值基本控制在0.2°以內(nèi),側(cè)滑角的誤差值基本控制在0.01°以內(nèi),馬赫數(shù)的誤差值基本控制在0.001以內(nèi),靜壓的誤差值控制在30 Pa以內(nèi),動(dòng)壓控制在85 Pa以內(nèi)。仿真結(jié)果表明,本文基于半圓柱體智能探頭設(shè)計(jì)的大氣數(shù)據(jù)系統(tǒng)算法模型解算精度達(dá)到了預(yù)期要求,且迭代收斂穩(wěn)定,能夠基本滿足工程使用上的要求。

    圖20 動(dòng)壓解算結(jié)果

    3.2 誤差分析

    筆者設(shè)計(jì)的基于半圓柱體智能探頭的大氣數(shù)據(jù)系統(tǒng)算法模型,主要誤差來源于在正向建模時(shí),由于樣本數(shù)量的限制帶來的數(shù)據(jù)擬合誤差,在模型設(shè)計(jì)的后續(xù)階段會(huì)采取繼續(xù)通過CFD計(jì)算或風(fēng)洞試驗(yàn)增加樣本數(shù)量[14],提高氣動(dòng)建模精度。

    4 結(jié)束語

    本文設(shè)計(jì)了一種基于新型半圓柱形智能探頭的大氣數(shù)據(jù)系統(tǒng)算法模型,介紹了半圓柱探頭外形信息和安裝要求,詳細(xì)說明了算法模型的建模思路和過程,并通過仿真驗(yàn)證證明了該算法模型精度滿足設(shè)計(jì)要求,迭代收斂穩(wěn)定。該種模型算法的研究設(shè)計(jì)方法可以為后續(xù)相關(guān)設(shè)計(jì)工作提供參考,具有較高的理論參考價(jià)值和實(shí)際應(yīng)用價(jià)值[15]。

    猜你喜歡
    馬赫數(shù)迎角靜壓
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    靜壓法沉樁對周邊環(huán)境影響及質(zhì)量控制
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    載荷分布對可控?cái)U(kuò)散葉型性能的影響
    靜壓托換樁在某濕陷性黃土場地地基加固中的應(yīng)用
    超精密液體靜壓轉(zhuǎn)臺(tái)裝配技術(shù)
    一種基于空氣靜壓支承的自調(diào)心裝置
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    半柔壁噴管初步實(shí)驗(yàn)研究
    99久久九九国产精品国产免费| 国产伦理片在线播放av一区| 国产伦在线观看视频一区| 网址你懂的国产日韩在线| 日韩av在线免费看完整版不卡| 亚洲av男天堂| 亚洲综合精品二区| 午夜久久久久精精品| 高清视频免费观看一区二区 | 国产精品日韩av在线免费观看| 在线观看av片永久免费下载| 日韩精品青青久久久久久| 国产一区亚洲一区在线观看| 午夜福利视频1000在线观看| 国产大屁股一区二区在线视频| 亚洲丝袜综合中文字幕| 免费观看av网站的网址| xxx大片免费视频| 亚洲高清免费不卡视频| 亚洲图色成人| 精品国产露脸久久av麻豆 | 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| 别揉我奶头 嗯啊视频| 在线观看人妻少妇| 天堂网av新在线| 欧美一级a爱片免费观看看| 人人妻人人看人人澡| 黑人高潮一二区| 国产黄色视频一区二区在线观看| 午夜福利视频1000在线观看| 国产精品美女特级片免费视频播放器| 久久久久精品性色| 少妇人妻精品综合一区二区| 男人狂女人下面高潮的视频| 久久精品国产鲁丝片午夜精品| 真实男女啪啪啪动态图| 精品人妻一区二区三区麻豆| 精品国产一区二区三区久久久樱花 | 老司机影院成人| 国产精品麻豆人妻色哟哟久久 | 熟女人妻精品中文字幕| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 日韩欧美三级三区| 男女视频在线观看网站免费| 亚洲最大成人手机在线| 男人狂女人下面高潮的视频| 免费黄频网站在线观看国产| 午夜免费男女啪啪视频观看| 91狼人影院| 国产精品一区二区在线观看99 | 国产成年人精品一区二区| 国产精品精品国产色婷婷| 色播亚洲综合网| 久久午夜福利片| 国产麻豆成人av免费视频| 国产精品伦人一区二区| 精品亚洲乱码少妇综合久久| 老女人水多毛片| 国产精品人妻久久久久久| 国产av码专区亚洲av| 美女内射精品一级片tv| 国产高潮美女av| 可以在线观看毛片的网站| 一级av片app| 午夜久久久久精精品| 精品酒店卫生间| 亚洲欧美一区二区三区黑人 | 国产精品一区www在线观看| 干丝袜人妻中文字幕| 色尼玛亚洲综合影院| 国产日韩欧美在线精品| 欧美日本视频| 国产精品1区2区在线观看.| 边亲边吃奶的免费视频| 日韩人妻高清精品专区| 男女下面进入的视频免费午夜| 一级a做视频免费观看| 久久国内精品自在自线图片| 80岁老熟妇乱子伦牲交| 亚洲乱码一区二区免费版| 国产精品女同一区二区软件| 天堂俺去俺来也www色官网 | 18禁在线播放成人免费| 亚洲美女搞黄在线观看| 九九在线视频观看精品| 久久99热6这里只有精品| 亚洲av成人av| 人妻少妇偷人精品九色| 国产成人一区二区在线| 一区二区三区高清视频在线| 国产伦精品一区二区三区视频9| 亚洲熟妇中文字幕五十中出| 成人特级av手机在线观看| 日日摸夜夜添夜夜爱| 国产在视频线精品| 欧美xxxx性猛交bbbb| 日韩一区二区三区影片| 久久精品国产鲁丝片午夜精品| 在线播放无遮挡| 欧美日韩综合久久久久久| 我的老师免费观看完整版| 在线观看av片永久免费下载| 免费av毛片视频| 国产精品爽爽va在线观看网站| 最后的刺客免费高清国语| 久久精品久久久久久久性| 免费黄频网站在线观看国产| 久久热精品热| 啦啦啦中文免费视频观看日本| 午夜福利成人在线免费观看| 免费观看av网站的网址| 久久综合国产亚洲精品| videos熟女内射| 永久网站在线| 青春草亚洲视频在线观看| 国产精品久久久久久久久免| 你懂的网址亚洲精品在线观看| 国产男人的电影天堂91| 亚洲一区高清亚洲精品| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕| 在线观看免费高清a一片| 精品久久久久久久久久久久久| 全区人妻精品视频| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 国产一区二区三区综合在线观看 | 亚洲欧美精品自产自拍| 亚洲最大成人av| 欧美丝袜亚洲另类| 国产精品一区二区在线观看99 | 国产久久久一区二区三区| 最近中文字幕2019免费版| 国产欧美另类精品又又久久亚洲欧美| 国产黄色视频一区二区在线观看| 久久久久国产网址| 国产欧美另类精品又又久久亚洲欧美| a级毛色黄片| 国产毛片a区久久久久| 色综合亚洲欧美另类图片| 男人和女人高潮做爰伦理| 极品少妇高潮喷水抽搐| 亚洲在线观看片| 欧美日韩一区二区视频在线观看视频在线 | 久久久午夜欧美精品| 久久精品久久精品一区二区三区| 真实男女啪啪啪动态图| 深夜a级毛片| 毛片一级片免费看久久久久| 欧美潮喷喷水| 一级黄片播放器| 免费黄色在线免费观看| 精品少妇黑人巨大在线播放| 亚洲av成人av| 又爽又黄无遮挡网站| 亚洲在久久综合| 蜜臀久久99精品久久宅男| 综合色av麻豆| 又爽又黄a免费视频| 身体一侧抽搐| 亚洲国产精品成人综合色| 老女人水多毛片| 日韩一本色道免费dvd| 少妇裸体淫交视频免费看高清| 内地一区二区视频在线| 日韩一本色道免费dvd| 一区二区三区四区激情视频| 日韩视频在线欧美| 婷婷色综合大香蕉| 好男人在线观看高清免费视频| 91精品一卡2卡3卡4卡| 日韩伦理黄色片| 国产在视频线精品| 天天一区二区日本电影三级| 亚洲av免费在线观看| 免费黄色在线免费观看| 高清午夜精品一区二区三区| 嫩草影院新地址| 中文字幕人妻熟人妻熟丝袜美| 1000部很黄的大片| 一二三四中文在线观看免费高清| 99热这里只有是精品50| 亚洲精品aⅴ在线观看| 少妇丰满av| 免费大片18禁| 久久精品国产亚洲网站| 又爽又黄a免费视频| 国产亚洲精品av在线| 一级毛片黄色毛片免费观看视频| 我要看日韩黄色一级片| 九九爱精品视频在线观看| 国内精品美女久久久久久| 乱人视频在线观看| 国产一区有黄有色的免费视频 | 大片免费播放器 马上看| 青春草视频在线免费观看| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久av不卡| 91精品一卡2卡3卡4卡| 亚洲怡红院男人天堂| 亚洲人与动物交配视频| 色网站视频免费| or卡值多少钱| 久久久久久久久久久丰满| 22中文网久久字幕| 高清欧美精品videossex| 欧美高清成人免费视频www| 22中文网久久字幕| 国产精品一区二区三区四区免费观看| 精品久久久噜噜| 欧美激情久久久久久爽电影| 免费看a级黄色片| av一本久久久久| 99久久精品国产国产毛片| 一级av片app| 日韩国内少妇激情av| 日本与韩国留学比较| 日本猛色少妇xxxxx猛交久久| 国内精品宾馆在线| 尤物成人国产欧美一区二区三区| 成人高潮视频无遮挡免费网站| 国产精品99久久久久久久久| 大香蕉久久网| 亚洲精品aⅴ在线观看| 欧美zozozo另类| 国产探花在线观看一区二区| www.色视频.com| 三级毛片av免费| h日本视频在线播放| 狂野欧美激情性xxxx在线观看| 极品少妇高潮喷水抽搐| 亚洲乱码一区二区免费版| 99久国产av精品国产电影| 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 国产一区有黄有色的免费视频 | 久久精品夜色国产| av.在线天堂| 在线观看美女被高潮喷水网站| 欧美一级a爱片免费观看看| 黄片wwwwww| 男人爽女人下面视频在线观看| 少妇熟女aⅴ在线视频| 久久久国产一区二区| 欧美高清性xxxxhd video| 亚洲精品乱码久久久久久按摩| 亚洲高清免费不卡视频| 精品酒店卫生间| 国产成人福利小说| 真实男女啪啪啪动态图| 美女大奶头视频| 国产成人aa在线观看| 精品少妇黑人巨大在线播放| 日韩三级伦理在线观看| 免费观看性生交大片5| 亚洲精品乱码久久久v下载方式| 黄色日韩在线| 亚洲熟女精品中文字幕| 夫妻午夜视频| 天天躁日日操中文字幕| 亚洲精品国产av蜜桃| 小蜜桃在线观看免费完整版高清| 天天躁夜夜躁狠狠久久av| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| 尤物成人国产欧美一区二区三区| 九九在线视频观看精品| 七月丁香在线播放| 毛片女人毛片| 久久久久性生活片| 亚洲在久久综合| 国产三级在线视频| 免费观看性生交大片5| 婷婷色av中文字幕| 永久网站在线| 亚洲va在线va天堂va国产| 久久鲁丝午夜福利片| 夫妻午夜视频| 久久这里有精品视频免费| 高清日韩中文字幕在线| 久热久热在线精品观看| 国产黄频视频在线观看| 日韩 亚洲 欧美在线| 精品不卡国产一区二区三区| 青青草视频在线视频观看| 青春草视频在线免费观看| 亚洲国产精品专区欧美| 熟妇人妻不卡中文字幕| 欧美日韩精品成人综合77777| 亚洲人与动物交配视频| 天天躁夜夜躁狠狠久久av| 禁无遮挡网站| 欧美区成人在线视频| 欧美激情在线99| 成人一区二区视频在线观看| 欧美 日韩 精品 国产| 在线免费观看不下载黄p国产| 肉色欧美久久久久久久蜜桃 | 午夜精品在线福利| 亚洲电影在线观看av| 日日摸夜夜添夜夜爱| 亚洲性久久影院| 夫妻午夜视频| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| av黄色大香蕉| 亚洲av.av天堂| 成人亚洲精品av一区二区| 久久久午夜欧美精品| 在线免费十八禁| 国产精品久久久久久久电影| 国模一区二区三区四区视频| 亚洲av二区三区四区| 综合色丁香网| 十八禁网站网址无遮挡 | 久久久久久久久久人人人人人人| 美女内射精品一级片tv| av福利片在线观看| 丰满少妇做爰视频| av天堂中文字幕网| 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 久久6这里有精品| 99re6热这里在线精品视频| 国产高清国产精品国产三级 | 国产一区有黄有色的免费视频 | 中文字幕亚洲精品专区| 哪个播放器可以免费观看大片| 天堂俺去俺来也www色官网 | 成人性生交大片免费视频hd| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 亚洲成人av在线免费| 午夜福利视频精品| 白带黄色成豆腐渣| 午夜福利视频精品| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩卡通动漫| 不卡视频在线观看欧美| 一个人免费在线观看电影| 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 日韩电影二区| a级毛片免费高清观看在线播放| 亚洲丝袜综合中文字幕| 又大又黄又爽视频免费| 偷拍熟女少妇极品色| 亚洲av中文字字幕乱码综合| 男女国产视频网站| av在线亚洲专区| 日韩强制内射视频| 成人综合一区亚洲| 久久人人爽人人爽人人片va| 青春草视频在线免费观看| 国产亚洲91精品色在线| 国产探花在线观看一区二区| 男女边吃奶边做爰视频| 成人午夜精彩视频在线观看| 舔av片在线| 波多野结衣巨乳人妻| 国产av不卡久久| 免费观看精品视频网站| 少妇的逼好多水| 成人一区二区视频在线观看| 亚洲精品国产成人久久av| 国产午夜精品一二区理论片| 深爱激情五月婷婷| 内射极品少妇av片p| 国产精品女同一区二区软件| 看免费成人av毛片| 日韩av在线大香蕉| 亚洲av二区三区四区| 精品久久久久久久末码| 亚洲成人av在线免费| 18禁裸乳无遮挡免费网站照片| 美女脱内裤让男人舔精品视频| 性色avwww在线观看| 国产熟女欧美一区二区| 亚洲精品日本国产第一区| 亚洲欧美一区二区三区黑人 | 亚洲国产精品sss在线观看| 成人二区视频| 欧美日韩精品成人综合77777| 18禁动态无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品乱码久久久v下载方式| 99久久精品热视频| 久久久a久久爽久久v久久| 亚洲av男天堂| 国产探花在线观看一区二区| 国产精品爽爽va在线观看网站| 国产精品不卡视频一区二区| 日本免费a在线| 久久久久国产网址| 97热精品久久久久久| 亚洲国产高清在线一区二区三| 亚洲四区av| 一区二区三区高清视频在线| 麻豆乱淫一区二区| 久久国产乱子免费精品| 黄色欧美视频在线观看| 精品久久久久久电影网| 国产亚洲av片在线观看秒播厂 | 欧美日韩亚洲高清精品| 免费av毛片视频| 乱码一卡2卡4卡精品| 最近最新中文字幕免费大全7| 国产大屁股一区二区在线视频| 久久这里只有精品中国| 国产高清三级在线| 欧美另类一区| 色5月婷婷丁香| 亚洲精品久久午夜乱码| 欧美不卡视频在线免费观看| 婷婷色综合www| 麻豆久久精品国产亚洲av| 97热精品久久久久久| 亚洲精华国产精华液的使用体验| 97人妻精品一区二区三区麻豆| 国产精品国产三级专区第一集| 哪个播放器可以免费观看大片| 欧美xxxx性猛交bbbb| 国产有黄有色有爽视频| 成人特级av手机在线观看| 亚洲三级黄色毛片| 婷婷六月久久综合丁香| 久久久久网色| 国产白丝娇喘喷水9色精品| 搡老乐熟女国产| 赤兔流量卡办理| 午夜福利在线在线| 秋霞在线观看毛片| 免费看a级黄色片| 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 精品久久久久久久末码| 成人亚洲精品一区在线观看 | 成人欧美大片| 少妇丰满av| 亚洲精品国产av蜜桃| 日韩精品有码人妻一区| 免费看av在线观看网站| av网站免费在线观看视频 | 99久久精品一区二区三区| 免费播放大片免费观看视频在线观看| 一级av片app| 久久久午夜欧美精品| 水蜜桃什么品种好| 欧美xxⅹ黑人| 2021少妇久久久久久久久久久| 中文精品一卡2卡3卡4更新| 亚洲自偷自拍三级| 亚洲图色成人| 国产不卡一卡二| 国产亚洲一区二区精品| 亚洲精品自拍成人| 五月天丁香电影| 成年女人看的毛片在线观看| 一区二区三区免费毛片| 人妻一区二区av| 最近视频中文字幕2019在线8| 免费观看精品视频网站| 亚洲精品色激情综合| 欧美 日韩 精品 国产| 日日啪夜夜撸| 国产片特级美女逼逼视频| 亚洲成人av在线免费| av卡一久久| 国产精品99久久久久久久久| or卡值多少钱| 可以在线观看毛片的网站| 国产成人精品婷婷| 在线观看人妻少妇| www.色视频.com| 九九爱精品视频在线观看| 欧美日韩精品成人综合77777| 欧美丝袜亚洲另类| 亚洲精品自拍成人| 国产精品美女特级片免费视频播放器| 成人毛片60女人毛片免费| 久久国内精品自在自线图片| 欧美区成人在线视频| 亚洲精品影视一区二区三区av| 九草在线视频观看| 日日摸夜夜添夜夜爱| 国产一区二区三区综合在线观看 | 一级爰片在线观看| 日韩不卡一区二区三区视频在线| 免费在线观看成人毛片| 国产高清不卡午夜福利| 好男人视频免费观看在线| 国内精品一区二区在线观看| 国产av国产精品国产| 欧美激情久久久久久爽电影| 观看美女的网站| 日韩国内少妇激情av| 亚洲av一区综合| 国产精品人妻久久久久久| 综合色av麻豆| 赤兔流量卡办理| 一级毛片电影观看| 久久久欧美国产精品| 亚洲国产欧美人成| 国产午夜精品一二区理论片| 舔av片在线| av在线亚洲专区| 一个人免费在线观看电影| 日韩国内少妇激情av| 最近中文字幕高清免费大全6| 日韩国内少妇激情av| 久久精品久久精品一区二区三区| 不卡视频在线观看欧美| 一区二区三区四区激情视频| 国产精品一区二区三区四区免费观看| 欧美激情在线99| 亚洲第一区二区三区不卡| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕| 亚洲美女搞黄在线观看| 色综合站精品国产| 岛国毛片在线播放| 五月天丁香电影| 精品人妻视频免费看| 精品人妻一区二区三区麻豆| 国产男人的电影天堂91| 亚洲欧美日韩东京热| 日韩亚洲欧美综合| 中文欧美无线码| 国产高清三级在线| 天美传媒精品一区二区| 大话2 男鬼变身卡| 欧美成人精品欧美一级黄| 高清日韩中文字幕在线| 国产午夜精品论理片| 啦啦啦中文免费视频观看日本| 三级国产精品欧美在线观看| 18禁动态无遮挡网站| 亚洲内射少妇av| 色综合色国产| 狂野欧美白嫩少妇大欣赏| 婷婷色综合大香蕉| 久久这里有精品视频免费| 午夜福利视频精品| 亚洲av中文av极速乱| 国产乱人偷精品视频| 我要看日韩黄色一级片| 欧美日韩一区二区视频在线观看视频在线 | 一个人免费在线观看电影| 欧美xxⅹ黑人| 中文字幕亚洲精品专区| 欧美另类一区| 美女脱内裤让男人舔精品视频| 天堂俺去俺来也www色官网 | 久久精品国产亚洲av天美| 男女国产视频网站| 九九在线视频观看精品| 亚洲成人一二三区av| 99久久精品热视频| 日日摸夜夜添夜夜添av毛片| 国语对白做爰xxxⅹ性视频网站| 国产精品福利在线免费观看| 亚洲在线观看片| 69人妻影院| 两个人的视频大全免费| 高清毛片免费看| 久久久精品94久久精品| 国产精品一区二区在线观看99 | 国产精品伦人一区二区| 九色成人免费人妻av| 欧美不卡视频在线免费观看| 国产中年淑女户外野战色| 久久久成人免费电影| 91狼人影院| 五月天丁香电影| 精品国产三级普通话版| 国产在线男女| 我要看日韩黄色一级片| 国产成人精品久久久久久| 亚洲av电影不卡..在线观看| 91久久精品国产一区二区三区| 亚洲成色77777| 成年免费大片在线观看| 午夜激情久久久久久久| 内射极品少妇av片p| 人妻少妇偷人精品九色| 久久这里只有精品中国| 色5月婷婷丁香| 国产老妇伦熟女老妇高清| 亚洲美女视频黄频| 人人妻人人看人人澡| 美女cb高潮喷水在线观看| 亚洲熟女精品中文字幕| 免费人成在线观看视频色| 亚洲婷婷狠狠爱综合网| 亚洲自拍偷在线| 真实男女啪啪啪动态图| 亚洲最大成人av| 国产精品.久久久| 国产在视频线精品| 亚洲av中文av极速乱| 国产爱豆传媒在线观看| 国产成人91sexporn| 国产色婷婷99| 亚洲最大成人av| 伦理电影大哥的女人| 大香蕉97超碰在线| 人妻夜夜爽99麻豆av| 色视频www国产| 亚洲欧美成人综合另类久久久| 晚上一个人看的免费电影| 免费电影在线观看免费观看| 国产午夜精品一二区理论片| av在线播放精品|