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

    面向復(fù)雜構(gòu)型飛機(jī)的非定常氣動(dòng)力建模與辨識(shí)

    2016-11-14 00:57:49楊文卜忱眭建軍尚祖銘
    航空學(xué)報(bào) 2016年8期
    關(guān)鍵詞:時(shí)間常數(shù)氣動(dòng)力迎角

    楊文, 卜忱, 眭建軍, 尚祖銘

    中航工業(yè)空氣動(dòng)力研究院 氣動(dòng)發(fā)展部, 哈爾濱 150001

    ?

    面向復(fù)雜構(gòu)型飛機(jī)的非定常氣動(dòng)力建模與辨識(shí)

    楊文, 卜忱*, 眭建軍, 尚祖銘

    中航工業(yè)空氣動(dòng)力研究院 氣動(dòng)發(fā)展部, 哈爾濱150001

    不論是現(xiàn)代高機(jī)動(dòng)隱身戰(zhàn)斗機(jī)的設(shè)計(jì)需求還是常規(guī)布局飛機(jī)的飛行動(dòng)力學(xué)分析,深入研究大迎角飛行時(shí)的非線性非定常氣動(dòng)力模型都極其重要。基于縱向運(yùn)動(dòng)小振幅及大振幅強(qiáng)迫振蕩試驗(yàn)數(shù)據(jù),分析了常規(guī)穩(wěn)定導(dǎo)數(shù)模型的準(zhǔn)確性,并從導(dǎo)數(shù)模型出發(fā)發(fā)展了簡(jiǎn)化渦流和分離流時(shí)間遲滯效應(yīng)的非定常氣動(dòng)力線性模型和非線性模型,最后應(yīng)用風(fēng)洞典型機(jī)動(dòng)歷程模擬試驗(yàn)驗(yàn)證了模型的有效性。結(jié)果表明:對(duì)于復(fù)雜構(gòu)型高機(jī)動(dòng)飛機(jī)模型,發(fā)展并改進(jìn)的非線性微分方程模型可以準(zhǔn)確預(yù)測(cè)飛機(jī)不同機(jī)動(dòng)下的非定常氣動(dòng)力特性,具有較強(qiáng)的工程可行性。

    大迎角; 非定常氣動(dòng)力; 時(shí)間遲滯; 微分方程; 典型機(jī)動(dòng)

    在大氣飛行力學(xué)中,確定并描述作用在飛機(jī)上的氣動(dòng)力及力矩是一項(xiàng)非常重要的任務(wù),而將飛行力學(xué)與其他力學(xué)分支區(qū)分開來的正是其中的氣動(dòng)力部分。嚴(yán)格地來講,氣動(dòng)力及力矩都是飛行狀態(tài)變量的泛函。大量的試驗(yàn)結(jié)果表明[1-3],他們不僅僅與這些飛行狀態(tài)變量的瞬時(shí)值有關(guān),還與這些變量在運(yùn)動(dòng)過程中的整個(gè)時(shí)間歷程相關(guān)。實(shí)際應(yīng)用中,僅采用這些變量及其導(dǎo)數(shù)的函數(shù)表征氣動(dòng)力及力矩,并拓展至關(guān)于這些變量的泰勒展開式。然而,如果將來要配備有多重控制器并具備大迎角區(qū)域、高機(jī)動(dòng)性的戰(zhàn)斗機(jī),那么在飛行動(dòng)力學(xué)分析及控制律設(shè)計(jì)中,就要求足夠準(zhǔn)確的非定常氣動(dòng)力模型。即使對(duì)于正常布局的飛機(jī),在其大迎角飛行條件下準(zhǔn)確預(yù)測(cè)其非定常氣動(dòng)力也是非常必要的。

    相對(duì)于數(shù)值模擬方法,基于風(fēng)洞試驗(yàn)數(shù)據(jù)建立非定常氣動(dòng)力模型仍然是目前非定常氣動(dòng)力建模研究的主要手段。目前大迎角非定常氣動(dòng)力建模方法中國(guó)內(nèi)外研究最多的主要有兩大類:一類是基于飛機(jī)表面流動(dòng)機(jī)理確定的非定常氣動(dòng)力和運(yùn)動(dòng)狀態(tài)變量的數(shù)學(xué)關(guān)系表達(dá)式,如穩(wěn)定導(dǎo)數(shù)模型[4]及其擴(kuò)展的多項(xiàng)式模型[5]、階躍響應(yīng)函數(shù)模型[6-7]及其簡(jiǎn)化模型[8]、狀態(tài)空間模型[9-11]和非線性微分方程模型[12-14]等;另一類是避開復(fù)雜的物理機(jī)理直接采用諸如模糊邏輯[15-16]、神經(jīng)網(wǎng)絡(luò)[17]及模糊神經(jīng)[18]等純數(shù)學(xué)的方法進(jìn)行非線性代數(shù)擬合建模。

    然而,目前在實(shí)際的飛行動(dòng)力學(xué)分析及飛行控制工程應(yīng)用中仍然采用穩(wěn)定導(dǎo)數(shù)模型,而上述關(guān)于氣動(dòng)力建模的研究還主要是集中在方法上,很少有人評(píng)價(jià)傳統(tǒng)導(dǎo)數(shù)模型表達(dá)式的準(zhǔn)確性,更談不上在此基礎(chǔ)上發(fā)展大迎角非定常氣動(dòng)力模型。其次,許多建模方法的介紹中,僅僅涉及到諸如大后掠三角翼等布局的氣動(dòng)力建模研究,而且很少應(yīng)用機(jī)動(dòng)歷程試驗(yàn)驗(yàn)證模型的有效性。實(shí)際的研究發(fā)現(xiàn),在工程中小展弦比復(fù)雜構(gòu)型飛機(jī)的氣動(dòng)力建模通常需要對(duì)模型結(jié)構(gòu)進(jìn)行修正并對(duì)辨識(shí)算法進(jìn)行改進(jìn)。本文[12]基于Goman提出的微分方程模型基本結(jié)構(gòu)形式,從分析常規(guī)穩(wěn)定導(dǎo)數(shù)模型的準(zhǔn)確性出發(fā),發(fā)展并改進(jìn)了同時(shí)適用于某小展弦比復(fù)雜構(gòu)型動(dòng)態(tài)標(biāo)模小振幅和大振幅運(yùn)動(dòng)的非定常氣動(dòng)力微分方程模型,最后應(yīng)用在FL-8低速風(fēng)洞中完成的尾沖機(jī)動(dòng)[19]模擬試驗(yàn)驗(yàn)證非線性微分方程模型的有效性和適用性。

    1 穩(wěn)定導(dǎo)數(shù)模型

    盡管動(dòng)導(dǎo)數(shù)與頻率相關(guān)性很強(qiáng),但在實(shí)際應(yīng)用中,大迎角區(qū)域的氣動(dòng)力模型仍然沿用式(1)所表達(dá)的常規(guī)導(dǎo)數(shù)模型:

    Ciδ(α,ω)δ

    (1)

    模型中導(dǎo)數(shù)與頻率的相關(guān)性使其與時(shí)間域的氣動(dòng)力建模任務(wù)不一致。在實(shí)際的穩(wěn)定性分析及飛行仿真中,式(1)中的頻率相關(guān)性往往被忽略。結(jié)果,減縮頻率ω隨意選取可能就會(huì)導(dǎo)致氣動(dòng)力導(dǎo)數(shù)模型具有極大地不確定性。圖1給出了小振幅俯仰振蕩法向力系數(shù)CN遲滯環(huán)和導(dǎo)數(shù)模型的對(duì)比圖,其中f表征頻率??梢钥闯鲈谑儆歉浇鼌^(qū)域,模型預(yù)測(cè)和試驗(yàn)結(jié)果存在較大差別。

    圖1 導(dǎo)數(shù)模型與試驗(yàn)數(shù)據(jù)的比較(f=1.0 Hz)Fig.1 Comparison between derivatives model and test data (f=1.0 Hz)

    2 線性微分方程模型

    2.1模型表達(dá)式推導(dǎo)

    針對(duì)第1節(jié)中導(dǎo)數(shù)模型遇到的問題,基于機(jī)翼表面流動(dòng)機(jī)理[11],Goman等提出了一種單自由度俯仰小幅振蕩線性微分方程模型[12],具體表達(dá)式為

    對(duì)于小振幅強(qiáng)迫振蕩試驗(yàn),迎角的變化規(guī)律為

    Ci(t)=Ci0+Ciααssin ω t+Ciqαsω cos ω t

    (4)

    式中:Ci0為氣動(dòng)力均值;Ci α和Ci q分別為所謂的“同相”及“異相”氣動(dòng)力導(dǎo)數(shù),并且在大迎角時(shí)與頻率、振幅相關(guān)。

    將式(3)代入氣動(dòng)力模型式(2)中,經(jīng)拉普拉斯變換、線性化、拉普拉斯反變換及歸并同類項(xiàng)后,與式(4)中同相及正交動(dòng)導(dǎo)數(shù)比較,整理可得

    (5)

    式中:ΔCiα,vb=Ciα,st-Ciα,att,從式(5)中消去非線性項(xiàng),可以得到同相導(dǎo)數(shù)與異相導(dǎo)數(shù)之間的線性關(guān)系式為

    Ciq=Ciq,att-τ(Ciα-Ciα,att)

    (6)

    式(5)還可以簡(jiǎn)化為

    (7)

    式中:

    (8)

    為待估計(jì)參數(shù),而

    (9)

    為已知量。

    從式(6)中能得到一個(gè)重要的結(jié)論,即以同相導(dǎo)數(shù)為自變量,異相導(dǎo)數(shù)為因變量,理論上在不同振蕩頻率下應(yīng)呈線性關(guān)系,而斜率就是這個(gè)迎角對(duì)應(yīng)的反映遲滯特性的時(shí)間常量,圖2給出了某復(fù)雜構(gòu)型飛機(jī)小振幅強(qiáng)迫振蕩的動(dòng)導(dǎo)數(shù)相圖,從圖中可見法向力系數(shù)和俯仰力矩系數(shù)對(duì)應(yīng)的同相和異相導(dǎo)數(shù)呈現(xiàn)出明顯的線性關(guān)系,特別是在大迎角下這種特征更加突出。

    圖2 不同中心迎角、不同頻率下的同相及異相導(dǎo)數(shù)Fig.2 In-phase and out-phase derivatives at different angles of attack and frequencies

    2.2線性模型參數(shù)辨識(shí)

    上述線性微分方程模型,在每個(gè)不同中心迎角處,特征時(shí)間常數(shù)τ采用線性回歸方法辨識(shí),而其他參數(shù)采用牛頓-拉夫遜算法最小化如式(10)所示的判據(jù)函數(shù)即可獲得未知參數(shù)b、c和d的估計(jì)結(jié)果。

    (10)

    式中:σ1為氣動(dòng)導(dǎo)數(shù)Ciα試驗(yàn)結(jié)果的相對(duì)誤差;σ2為氣動(dòng)導(dǎo)數(shù)Ciq試驗(yàn)結(jié)果的相對(duì)誤差。

    圖3給出了特征時(shí)間常數(shù)的辨識(shí)結(jié)果及參數(shù)估計(jì)的置信區(qū)間。圖4給出了所有參數(shù)辨識(shí)完后解微分方程得到的線性模型預(yù)測(cè)結(jié)果和小振幅振蕩氣動(dòng)力遲滯環(huán)的比較,在圖4(a)中CN,test為法向力系數(shù)小振幅遲滯環(huán)試驗(yàn)數(shù)據(jù),CN,model 1為法向力系數(shù)線性模型預(yù)測(cè)結(jié)果,CN,dyn為小振幅遲滯環(huán)中的非定常貢獻(xiàn)量,ΔCN為渦破裂引起的法向力系數(shù)差量,而CN,st為靜態(tài)試驗(yàn)數(shù)據(jù),CN,model 2為法向力系數(shù)導(dǎo)數(shù)模型結(jié)果,圖4(b)中的俯仰力矩元相關(guān)變量與法向力元類似。顯然,線性微分方程模型明顯改善了與大迎角小振幅試驗(yàn)氣動(dòng)遲滯環(huán)的擬合效果。

    圖3 特征時(shí)間常數(shù)隨迎角的變化曲線Fig.3 Characteristic time scale varying with angle of attack

    圖4 線性模型與試驗(yàn)數(shù)據(jù)的比較Fig.4 Comparison between linear model and test data

    3 非線性微分方程模型

    3.1模型結(jié)構(gòu)修正及辨識(shí)

    隨著振幅的增大,第2節(jié)中的線性數(shù)學(xué)模型可能失去其準(zhǔn)度,為了使微分方程模型能反映大振幅的非定常氣動(dòng)力,就需要對(duì)式(2)做適當(dāng)?shù)匦拚?,Goman的修正方程為

    (11)

    式中:n為多項(xiàng)式函數(shù)的階次,取n=3;方程中k1(α)與特征時(shí)間常數(shù)相關(guān),且當(dāng)k2(α)=k3(α)=0時(shí),k1(α)=1/τ(α),從本質(zhì)上講,上述修正是因?yàn)橛行卣鲿r(shí)間常數(shù)τeff不僅與迎角有關(guān),還與運(yùn)動(dòng)的時(shí)間歷程相關(guān),因而大振幅非定常氣動(dòng)力建模中需要辨識(shí)新的有效特征時(shí)間常數(shù)τeff。然而,不同于尖前緣大后掠三角翼[11],在復(fù)雜構(gòu)型布局的飛機(jī)模型氣動(dòng)力建模研究中發(fā)現(xiàn),k2(α)和k3(α)參數(shù)對(duì)大振幅非定常氣動(dòng)力模型辨識(shí)精度影響并不大,針對(duì)實(shí)際情況,本文對(duì)模型進(jìn)行如下形式的修正:

    (12)

    式中:Ci,st(α)為靜態(tài)氣動(dòng)力系數(shù);Ci,uns為非定常非線性氣動(dòng)力增量;τ1為反映流動(dòng)拓?fù)浣Y(jié)構(gòu)遲滯效應(yīng)的特征時(shí)間常數(shù),且與迎角有關(guān),如對(duì)于前機(jī)身帶邊條翼的戰(zhàn)斗機(jī)布局在較小迎角處流場(chǎng)始終是脫體渦結(jié)構(gòu),并沒有出現(xiàn)渦的破裂,因而特征時(shí)間常數(shù)通常較小,而在失速迎角附近,由于脫體渦的破裂與再附相對(duì)于運(yùn)動(dòng)本身會(huì)出現(xiàn)明顯遲滯,因而此時(shí)的特征時(shí)間常數(shù)將會(huì)顯著增大。

    (13)

    式中:τ1(α)、Aij為待辨識(shí)參數(shù)。實(shí)際上相關(guān)文獻(xiàn)[20]研究表明,基于運(yùn)動(dòng)狀態(tài)變量泛函的泰勒展開式的微分方程表達(dá)式是階躍響應(yīng)函數(shù)的一種簡(jiǎn)化形式,因而上述模型結(jié)構(gòu)的修正有其理論依據(jù)。當(dāng)模型的主體結(jié)構(gòu)確定后,接下來要做的是確定式(13)中右邊泰勒展開式的階次,對(duì)不同階次n計(jì)算模型的殘差平方和SEE,SEE最小的模型即視為最優(yōu)模型階次。研究發(fā)現(xiàn),隨著n的增大,SEE逐漸變小,但當(dāng)n達(dá)到一定數(shù)值時(shí),SEE變化的速率顯著減小,而此時(shí)如果繼續(xù)增大階次n反而會(huì)極大地增加模型辨識(shí)的復(fù)雜程度,綜合考慮后本文取n=6。

    3.2非線性微分方程模型參數(shù)估計(jì)

    數(shù)學(xué)模型結(jié)構(gòu)確定后,問題就成了根據(jù)辨識(shí)準(zhǔn)則和試驗(yàn)數(shù)據(jù)求取模型中的待定參數(shù),即參數(shù)估計(jì)問題,這是系統(tǒng)辨識(shí)定量研究的核心。而參數(shù)估計(jì)包括辨識(shí)準(zhǔn)則和優(yōu)化算法兩部分。本文選用最大似然法準(zhǔn)則作為辨識(shí)準(zhǔn)則,選用牛頓-拉夫遜算法作為優(yōu)化算法。實(shí)際上本文的氣動(dòng)力建模問題中只有一個(gè)觀測(cè)量且不考慮過程噪聲,因而最大似然準(zhǔn)則和最小平方差和準(zhǔn)則相同。

    實(shí)踐中發(fā)現(xiàn),牛頓-拉夫遜法雖然具有不低于二階的收斂速度,但對(duì)于非線性強(qiáng)的氣動(dòng)力模型容易發(fā)散,因?yàn)樵撍惴ㄒ竽繕?biāo)函數(shù)的Hesse矩陣G(θ)在每個(gè)迭代點(diǎn)θk處是正定的,否則難以保證牛頓方向下降,為了克服這一缺陷,本文引進(jìn)阻尼因子μk=0.5,使得矩陣Ak=G(θk)+μkI正定。

    3.3算例分析

    現(xiàn)在采用上述非線性微分方程模型和參數(shù)辨識(shí)算法來處理FL-8低速風(fēng)洞中完成的某復(fù)雜構(gòu)型飛機(jī)大幅諧波振蕩試驗(yàn)數(shù)據(jù),以檢驗(yàn)上述修正模型結(jié)構(gòu)及辨識(shí)算法的可行性。選取振幅為40°,頻率分別為0.2、0.4、0.6和0.8Hz的試驗(yàn)數(shù)據(jù)及振幅為30°,頻率分別為0.2、0.4和0.8Hz的試驗(yàn)數(shù)據(jù)作為訓(xùn)練樣本,而選取振幅為30°,頻率0.6Hz的振蕩試驗(yàn)數(shù)據(jù)及振幅為20°,頻率為0.2、0.4、0.6和0.8Hz的試驗(yàn)數(shù)據(jù)作為檢驗(yàn)樣本。特征時(shí)間常數(shù)隨迎角變化,在0°~80° 迎角范圍內(nèi)將τ1(α)等間隔離散41個(gè)點(diǎn),在0°~25° 迎角范圍內(nèi),其初值取0.2~2.0,由前所述,在小迎角時(shí)翼面流場(chǎng)結(jié)構(gòu)不變,流場(chǎng)內(nèi)部渦系拓?fù)浣Y(jié)構(gòu)調(diào)整的滯后小,因而時(shí)間尺度小。而在25° 迎角后,渦破裂的出現(xiàn)將使渦系拓?fù)浣Y(jié)構(gòu)調(diào)整的滯后顯著增大,其初值取線性模型中τ(α)的辨識(shí)結(jié)果,其他參數(shù)初值全取1。另外,由流動(dòng)的物理特性可知該參數(shù)值不小于零,如果參數(shù)的辨識(shí)結(jié)果出現(xiàn)負(fù)值,即使收斂也將失去其本身的物理意義,實(shí)踐發(fā)現(xiàn),在30° 迎角后,特征時(shí)間常數(shù)通常是一個(gè)較大的正數(shù),而在小迎角時(shí)可能出現(xiàn)負(fù)值,且參數(shù)出現(xiàn)負(fù)值時(shí),雖然模型的精度較高,但其預(yù)測(cè)能力一般較差,為了避免此時(shí)該參數(shù)出現(xiàn)負(fù)值,在迭代的過程中人為地限制其值大于等于零,從而保證收斂后辨識(shí)結(jié)果不至于喪失其物理意義。圖5給出了法向力系數(shù)參數(shù)辨識(shí)迭代收斂曲線。圖中SEE表征殘差平方和,而RR表征多層相關(guān)系數(shù)。

    圖5 迭代收斂曲線Fig.5 Iterative convergence curves

    圖6和圖7給出了模型輸出和訓(xùn)練樣本的比較,其中Test表征大幅振蕩試驗(yàn)數(shù)據(jù),而Model表征非線性微分方程模型計(jì)算數(shù)據(jù)。從圖中可以發(fā)現(xiàn),經(jīng)過修正后的非線性氣動(dòng)力模型與大幅振蕩俯仰運(yùn)動(dòng)的試驗(yàn)數(shù)據(jù)吻合得很好。

    圖6 非定常模型與試驗(yàn)數(shù)據(jù)的比較(f=0.2 Hz)Fig.6 Comparison between unsteady model and test data (f=0.2 Hz)

    圖7 非定常模型與試驗(yàn)數(shù)據(jù)的比較(f=0.6 Hz)Fig.7 Comparison between unsteady model and test data (f=0.6 Hz)

    圖8給出了模型預(yù)測(cè)和檢驗(yàn)樣本的比較,雖然模型的預(yù)測(cè)精度比建模精度略差,但依然能比較準(zhǔn)確地預(yù)測(cè)。這就說明上述的模型結(jié)構(gòu)辨識(shí)與參數(shù)估計(jì)足夠準(zhǔn)確。

    圖8 非定常模型與試驗(yàn)數(shù)據(jù)30° 振幅值比較(f=0.6 Hz)Fig.8 Comparison between unsteady model and test data at 30° of amplitude (f=0.6 Hz)

    3.4典型機(jī)動(dòng)試驗(yàn)驗(yàn)證

    為了進(jìn)一步檢驗(yàn)及評(píng)估上述非線性微分方程模型的有效性,這里將應(yīng)用尾沖機(jī)動(dòng)風(fēng)洞模擬試驗(yàn)來驗(yàn)證。圖9描述了尾沖機(jī)動(dòng)的時(shí)間歷程。圖10 和圖11給出了非線性微分方程模型預(yù)測(cè)和機(jī)動(dòng)歷程試驗(yàn)數(shù)據(jù)的比較,其中t表征機(jī)動(dòng)歷程的時(shí)間??梢婎A(yù)測(cè)結(jié)果和試驗(yàn)數(shù)據(jù)的一致性較好。

    鑒于上述訓(xùn)練樣本、檢驗(yàn)樣本及典型機(jī)動(dòng)試驗(yàn)數(shù)據(jù)均與模型預(yù)測(cè)擬合較好,這就說明針對(duì)某復(fù)雜構(gòu)型飛機(jī)布局的工程應(yīng)用研究,本文建立的改進(jìn)非線性微分方程模型能夠較準(zhǔn)確描述其大迎角非定常氣動(dòng)力特性。

    圖9 尾沖機(jī)動(dòng)Fig.9 Tail slide maneuvers

    圖10 非定常模型與尾沖機(jī)動(dòng)試驗(yàn)數(shù)據(jù)的比較(t=5 s)Fig.10 Comparison between unsteady model and test data in tail slide maneuvers (t=5 s)

    圖11 非定常模型與尾沖機(jī)動(dòng)試驗(yàn)數(shù)據(jù)比較(t=10 s)Fig.11 Comparison between unsteady model and test data in tail slide maneuvers (t=10 s)

    4 結(jié) 論

    1) 對(duì)于小振幅運(yùn)動(dòng)氣動(dòng)力遲滯環(huán),微分方程模型的線性化表達(dá)式具有較高的預(yù)測(cè)準(zhǔn)度,因而它為預(yù)測(cè)風(fēng)洞試驗(yàn)難以實(shí)現(xiàn)的高頻率動(dòng)導(dǎo)數(shù)提供了一種途徑。

    2) 雖然相較于諸如粒子群等智能方法,牛頓類辨識(shí)算法程序設(shè)計(jì)相對(duì)復(fù)雜,但它無需設(shè)置參數(shù)范圍且通過適當(dāng)?shù)男拚ǔ>哂辛己玫氖諗啃裕颐坎接忻鞔_的數(shù)學(xué)意義,適合面向工程的氣動(dòng)力建模研究。

    3) 本文發(fā)展的非線性微分方程模型結(jié)構(gòu)形式簡(jiǎn)單、物理意義清晰、工程針對(duì)性強(qiáng),基本能夠準(zhǔn)確反映復(fù)雜構(gòu)型飛機(jī)大迎角非定常氣動(dòng)力特性。

    [1]JARRAH M A. Low speed wind tunnel investigation of flow about delta wings oscillating in pitch to very high angle of attack: AIAA-1989-0295[R]. Reston: AIAA, 1989.

    [2]BRANDON J M, SHAH G H. Unsteady aerodynamic characteristics of fighter model undergoing large-amplitude pitching motions at high angles of attack: AIAA-1990-0309[R]. Reston: AIAA,1990.

    [3]ROBERT C N. The unsteady aerodynamics of slender wings and aircraft undergoing large amplitude maneuvers[J]. Progress in Aerospace Sciences, 2003, 39(3): 185-248.

    [4]ABZUG M J, LARRAB E. Airplane stability and control[M]. Cambridge: Cambridge University Press, 1977: 165.

    [5]LIN G F, LAN C E, BRANDON J M. A generalized dynamic aerodynamic coefficient model for flight dynamics applications: AIAA-1997-3643[R]. Reston: AIAA, 1997.

    [6]TOBAK M, SCHIFF L. On the formulation of the aerodynamic characteristics in aircraft dynamics: NASA TRR-456[R]. Washington, D.C.: NASA, 1976.

    [7]GHOREYSHI M, JIRASEK A, POST M L. A computational investigation into the use of response functions for a aerodynamic loads modeling: AIAA-2010-3518[R]. Reston: AIAA, 2010.

    [8]KLEIN V, MURPHY P C. Estimation of aircraft nonlinear unsteady parameters from wind tunnel data: NASA TM-1998-208[R]. Washington, D.C.: NASA, 1998.

    [9]GOMAN M G, KHRABROV A N. State-space representation of aerodynamic characteristics of an aircraft at high angels of attack[J]. Journal of Aircraft, 1994, 31(5): 1109-1115.

    [10]YI G F, FREDERICK H, LUTZE. Identification of unsteady aerodynamic modeling at high angle of attack: AIAA-1996-3407-CP[R]. Reston: AIAA, 1996.

    [11]楊文, 尚祖銘, 卜忱. 一種基于升力面分解思想的橫向狀態(tài)空間模型[J]. 飛行力學(xué), 2014, 32(5): 463-467.

    YANG W, SHANG Z M, BU C. An lateral-directional state-space model based the idea of splitting the aircraft into several lifting surfaces or panels[J]. Flight Dynamics, 2014, 32(5): 463-467 (in Chinese).

    [12]ABRAMOV M, GOMAN M G. Aircraft dynamics at high incidence flight with account of unsteady aerodynamic effects: AIAA-2004-5274[R]. Reston: AIAA,2004.

    [13]ABRAMOV N B, GOMAN M G, KHRABROV A N, et al. Simple wings unsteady aerodynamics at high angles of attack: Experimental and modeling results: AIAA-1999-4013[R]. Reston: AIAA,1999.

    [14]龔正, 沈宏良. 非定常氣動(dòng)力非線性微分方程建模方法[J]. 航空學(xué)報(bào), 2010, 32(1): 83-90.

    GONG Z, SHEN H L. Unsteady aerodynamic modeling method using nonlinear differential equations[J]. Acta Aeronautica et Astronautica Sinica, 2010, 32(1): 83-90 (in Chinese).

    [15]WANG Z J, LAN C E. Fuzzy logic modeling of lateral-directional unsteady aerodynamics: AIAA-1999-4012[R]. Reston: AIAA, 1999.

    [16]ZHONG J, WANG C, EDWARD L. Fuzzy logic modeling of nonlinear unsteady aerodynamics: AIAA-1998-4351[R]. Reston: AIAA,1998.

    [17]WANG Q, HE K F. Unsteady aerodynamics modeling for flight dynamics application[J]. Acta Mechanica Sinica, 2012, 28(1): 14-23.

    [18]史志偉, 明曉. 基于模糊聚類的模糊神經(jīng)網(wǎng)絡(luò)在非定常氣動(dòng)力建模中的應(yīng)用[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2005, 23(1): 21-24.

    SHI Z W, MING X. The application of FNN in unsteady aerodynamics modeling based on fuzzy clustering[J]. Acta Aerodynamic Sinica, 2005, 23(1): 21-24 (in Chinese).

    [19]ALCORN C W. The X-31 aircraft advances in aircraft agility and performance[J]. Progress Aerospace Science, 1996, 15(32): 377-413.

    [20]汪清, 蔡金獅.飛機(jī)大迎角非定常氣動(dòng)力建模與辨識(shí)[J]. 航空學(xué)報(bào), 1996, 17(4): 391-398.

    WANG Q, CAI J S. Unsteady aerodynamic modeling and identification of airplane at high angles of attack[J]. Acta Aeronautica et Astronautica Sinica, 1996, 17(4): 391-398 (in Chinese).

    楊文男, 碩士, 工程師。主要研究方向: 風(fēng)洞動(dòng)態(tài)試驗(yàn)技術(shù)及非定常氣動(dòng)力建模。

    Tel: 0451-87570255

    E-mail: yangwen19860804@163.com

    卜忱男, 碩士, 研究員。主要研究方向: 風(fēng)洞動(dòng)態(tài)試驗(yàn)技術(shù)。

    Tel: 0451-87571476

    E-mail: buchen.1975@126.com

    Unsteady aerodynamic modeling and identification fora complicated aircraft configurations

    YANG Wen, BU Chen*, SUI Jianjun, SHANG Zuming

    Department of Aerodynamics Development, AVIC Aerodynamics Research Institute, Harbin150001, China

    Adequate modeling of nonlinear and unsteady aerodynamics at high angle of attack flight is important for the design of future fighters with high maneuverability and stealth as well as for the improved prediction of normal aircraft configuration’s dynamics. The limitations for conventional aerodynamic derivatives model based on longitudinal small amplitude experimental date and large amplitude experimental date at high angle of attack had been analyzed. The dynamic linear and nonlinear aerodynamic model approximating the vertical and separated flow time lag effects is considered along with the conventional aerodynamic model. Using the model the unsteady aerodynamics of one aircraft in typical maneuver simulation tests is predicted. It is suggested that the structural modification of nonlinear differential equation model proposed in this paper is valid in different maneuvers for complicated aircraft configurations, which bears proof on the practicality of the flight dynamics analysis.

    large angle of attack; unsteady aerodynamic; time lag; differential equation; typical maneuver

    2016-01-25; Revised: 2016-02-15; Accepted: 2016-03-14; Published online: 2016-04-0611:39

    . Tel.: 0451-87571476E-mail: buchen.1975@126.com

    2016-01-25; 退修日期: 2016-02-15; 錄用日期: 2016-03-14;

    時(shí)間: 2016-04-0611:39

    www.cnki.net/kcms/detail/11.1929.V.20160406.1139.004.html

    .Tel.: 0451-87571476E-mail: buchen.1975@126.com

    10.7527/S1000-6893.2016.0094

    V212.1

    A

    1000-6893(2016)08-2464-08

    引用格式: 楊文,卜忱, 眭建軍, 等. 面向復(fù)雜構(gòu)型飛機(jī)的非定常氣動(dòng)力建模與辨識(shí) [J]. 航空學(xué)報(bào), 2016, 37(8): 2464-2471. YANG W, BU C, SUI J J, et al. Unsteady aerodynamic modeling and identification for a complicated aircraft configuration [J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2464-2471.

    http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

    URL: www.cnki.net/kcms/detail/11.1929.V.20160406.1139.004.html

    猜你喜歡
    時(shí)間常數(shù)氣動(dòng)力迎角
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    飛行載荷外部氣動(dòng)力的二次規(guī)劃等效映射方法
    熱電偶時(shí)間常數(shù)檢測(cè)分揀系統(tǒng)設(shè)計(jì)
    側(cè)風(fēng)對(duì)拍動(dòng)翅氣動(dòng)力的影響
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    瞬變電磁視時(shí)間常數(shù)tau成像分析與應(yīng)用研究
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    風(fēng)力機(jī)氣動(dòng)力不對(duì)稱故障建模與仿真
    關(guān)于一階電路時(shí)間常數(shù)求法的討論
    伺服電動(dòng)機(jī)的時(shí)間常數(shù)問題
    老汉色∧v一级毛片| 国语自产精品视频在线第100页| 国产精品美女特级片免费视频播放器 | 日本三级黄在线观看| 亚洲国产精品成人综合色| 欧美激情高清一区二区三区| 一二三四社区在线视频社区8| 午夜福利,免费看| 91精品三级在线观看| 国产xxxxx性猛交| 波多野结衣av一区二区av| 欧美一区二区精品小视频在线| ponron亚洲| 久久久久久人人人人人| 好男人电影高清在线观看| 动漫黄色视频在线观看| 激情视频va一区二区三区| 精品国产超薄肉色丝袜足j| 老司机在亚洲福利影院| 手机成人av网站| 亚洲av五月六月丁香网| 亚洲最大成人中文| 亚洲一区二区三区不卡视频| 搡老熟女国产l中国老女人| 99久久99久久久精品蜜桃| 欧美最黄视频在线播放免费| 日韩欧美国产在线观看| √禁漫天堂资源中文www| 日本vs欧美在线观看视频| 女警被强在线播放| av福利片在线| tocl精华| 91国产中文字幕| 啪啪无遮挡十八禁网站| 久久人妻熟女aⅴ| 99国产综合亚洲精品| 亚洲性夜色夜夜综合| 欧美黑人欧美精品刺激| 国产私拍福利视频在线观看| 亚洲国产毛片av蜜桃av| 国产男靠女视频免费网站| 亚洲精品一区av在线观看| 人人妻,人人澡人人爽秒播| 午夜福利18| 色在线成人网| 波多野结衣av一区二区av| 国产又爽黄色视频| 亚洲自拍偷在线| av欧美777| 男女下面进入的视频免费午夜 | 午夜久久久久精精品| 色综合欧美亚洲国产小说| 欧美中文综合在线视频| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久成人av| 午夜久久久在线观看| 熟女少妇亚洲综合色aaa.| 一卡2卡三卡四卡精品乱码亚洲| 美女扒开内裤让男人捅视频| 亚洲成a人片在线一区二区| 性少妇av在线| 成人国语在线视频| 人妻久久中文字幕网| 午夜成年电影在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 午夜久久久在线观看| 在线播放国产精品三级| 国内久久婷婷六月综合欲色啪| 欧美国产日韩亚洲一区| 老司机福利观看| 性欧美人与动物交配| 成熟少妇高潮喷水视频| 成人三级做爰电影| 欧美另类亚洲清纯唯美| 亚洲国产精品久久男人天堂| 欧美中文日本在线观看视频| 久久久久国产一级毛片高清牌| 99久久国产精品久久久| 国产极品粉嫩免费观看在线| 国产精品秋霞免费鲁丝片| 久久香蕉激情| 在线观看免费午夜福利视频| 久久人妻熟女aⅴ| 大陆偷拍与自拍| 男女之事视频高清在线观看| 亚洲av第一区精品v没综合| 叶爱在线成人免费视频播放| 国产av精品麻豆| 国产成人精品在线电影| 成在线人永久免费视频| 两个人视频免费观看高清| 日韩中文字幕欧美一区二区| 神马国产精品三级电影在线观看 | 伦理电影免费视频| 我的亚洲天堂| 啦啦啦观看免费观看视频高清 | 国产高清有码在线观看视频 | 乱人伦中国视频| 嫩草影院精品99| 欧美色欧美亚洲另类二区 | 黄片大片在线免费观看| 成人精品一区二区免费| 亚洲精品一区av在线观看| 免费女性裸体啪啪无遮挡网站| 欧美色欧美亚洲另类二区 | 给我免费播放毛片高清在线观看| 夜夜躁狠狠躁天天躁| 亚洲专区中文字幕在线| 免费观看人在逋| 亚洲欧美精品综合久久99| 欧美精品啪啪一区二区三区| 日本 av在线| 一本久久中文字幕| 日韩大尺度精品在线看网址 | 久久精品国产99精品国产亚洲性色 | 亚洲国产欧美日韩在线播放| 国产av精品麻豆| 午夜福利成人在线免费观看| 妹子高潮喷水视频| 欧美日韩亚洲综合一区二区三区_| 久久精品国产99精品国产亚洲性色 | 美女午夜性视频免费| 老鸭窝网址在线观看| 国产日韩一区二区三区精品不卡| 日本免费一区二区三区高清不卡 | 侵犯人妻中文字幕一二三四区| 在线免费观看的www视频| av福利片在线| 精品福利观看| 国产一区二区在线av高清观看| 精品无人区乱码1区二区| 国产精品乱码一区二三区的特点 | 中文字幕av电影在线播放| 国产精品久久久久久亚洲av鲁大| 高清黄色对白视频在线免费看| 国产高清videossex| 日本黄色视频三级网站网址| 欧美日韩亚洲综合一区二区三区_| 亚洲av电影不卡..在线观看| 12—13女人毛片做爰片一| 露出奶头的视频| 狂野欧美激情性xxxx| 欧美日韩精品网址| 91成年电影在线观看| 婷婷丁香在线五月| 午夜两性在线视频| videosex国产| 伦理电影免费视频| 午夜福利18| 中亚洲国语对白在线视频| 我的亚洲天堂| 亚洲美女黄片视频| 亚洲精华国产精华精| 人人澡人人妻人| 欧美+亚洲+日韩+国产| 免费看美女性在线毛片视频| 久久国产精品人妻蜜桃| av在线播放免费不卡| 视频在线观看一区二区三区| 国产99久久九九免费精品| 国产av在哪里看| 男女下面插进去视频免费观看| 最新美女视频免费是黄的| 好男人电影高清在线观看| 满18在线观看网站| 男男h啪啪无遮挡| 国产一区在线观看成人免费| 操出白浆在线播放| 亚洲一区中文字幕在线| 国产精品日韩av在线免费观看 | 最好的美女福利视频网| 青草久久国产| 搞女人的毛片| 99在线人妻在线中文字幕| 欧美黄色淫秽网站| 12—13女人毛片做爰片一| 性色av乱码一区二区三区2| 久久伊人香网站| 少妇粗大呻吟视频| 电影成人av| 丁香六月欧美| 国产精品久久久久久精品电影 | 在线国产一区二区在线| 国产精品,欧美在线| 国产国语露脸激情在线看| 久久久久国产精品人妻aⅴ院| 亚洲第一青青草原| 亚洲av电影不卡..在线观看| www.自偷自拍.com| 高清毛片免费观看视频网站| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| 国产精品精品国产色婷婷| 黄色丝袜av网址大全| 久久人人97超碰香蕉20202| 两个人视频免费观看高清| 国内精品久久久久精免费| 97人妻精品一区二区三区麻豆 | 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色 | 国产av一区二区精品久久| 久久久国产成人免费| 不卡av一区二区三区| 午夜福利,免费看| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 18禁美女被吸乳视频| 91老司机精品| 国产又爽黄色视频| 性少妇av在线| 国产欧美日韩一区二区三区在线| 嫁个100分男人电影在线观看| 9191精品国产免费久久| 在线免费观看的www视频| 国产精品亚洲av一区麻豆| 亚洲国产高清在线一区二区三 | 亚洲精品久久国产高清桃花| 天天一区二区日本电影三级 | 99国产精品一区二区三区| 亚洲成av片中文字幕在线观看| 夜夜夜夜夜久久久久| 欧美成人午夜精品| 一个人观看的视频www高清免费观看 | 欧美+亚洲+日韩+国产| 电影成人av| 午夜日韩欧美国产| 他把我摸到了高潮在线观看| 精品午夜福利视频在线观看一区| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 午夜福利影视在线免费观看| 国产一区二区激情短视频| 午夜福利成人在线免费观看| 欧美日本亚洲视频在线播放| 在线观看日韩欧美| 999久久久精品免费观看国产| 亚洲精品中文字幕一二三四区| 一级毛片女人18水好多| av视频免费观看在线观看| 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 满18在线观看网站| 国产1区2区3区精品| 男女做爰动态图高潮gif福利片 | 纯流量卡能插随身wifi吗| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频 | 脱女人内裤的视频| 久久亚洲精品不卡| 欧美黄色淫秽网站| 91字幕亚洲| 18禁裸乳无遮挡免费网站照片 | 亚洲成人免费电影在线观看| 亚洲黑人精品在线| 精品国产国语对白av| 亚洲成人国产一区在线观看| svipshipincom国产片| 免费看a级黄色片| 国产亚洲av嫩草精品影院| 欧美国产精品va在线观看不卡| 国产午夜精品久久久久久| 久久草成人影院| 老熟妇仑乱视频hdxx| 这个男人来自地球电影免费观看| 又黄又粗又硬又大视频| 国产三级黄色录像| 国内精品久久久久精免费| 国产高清激情床上av| 国产麻豆69| 国产99白浆流出| 69av精品久久久久久| 中文字幕高清在线视频| 国产欧美日韩精品亚洲av| 久久久久精品国产欧美久久久| 亚洲精品一区av在线观看| 国产熟女午夜一区二区三区| 亚洲欧美激情在线| 这个男人来自地球电影免费观看| 国产精品99久久99久久久不卡| 俄罗斯特黄特色一大片| 国产亚洲精品av在线| 精品乱码久久久久久99久播| 成人18禁在线播放| 日韩免费av在线播放| 久久人妻熟女aⅴ| 黄色毛片三级朝国网站| 亚洲美女黄片视频| 视频在线观看一区二区三区| 老司机在亚洲福利影院| 色播亚洲综合网| 咕卡用的链子| 免费高清在线观看日韩| 久久 成人 亚洲| 国产一卡二卡三卡精品| 日本三级黄在线观看| 国产精品久久久av美女十八| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区蜜桃| 国产av又大| 变态另类成人亚洲欧美熟女 | 亚洲欧美日韩另类电影网站| 日本 欧美在线| 欧美乱色亚洲激情| 免费高清在线观看日韩| 9191精品国产免费久久| 国产精品爽爽va在线观看网站 | 超碰成人久久| 午夜亚洲福利在线播放| 9191精品国产免费久久| 黄色视频,在线免费观看| 亚洲色图综合在线观看| 久久精品影院6| 最新美女视频免费是黄的| 亚洲片人在线观看| 女人被狂操c到高潮| 99国产极品粉嫩在线观看| 99在线视频只有这里精品首页| 少妇 在线观看| 日本a在线网址| 久久精品亚洲精品国产色婷小说| 亚洲九九香蕉| 色播在线永久视频| 午夜福利免费观看在线| 波多野结衣巨乳人妻| 成人av一区二区三区在线看| 国产亚洲精品综合一区在线观看 | 精品电影一区二区在线| 日本一区二区免费在线视频| 18禁观看日本| 又紧又爽又黄一区二区| 18美女黄网站色大片免费观看| 亚洲最大成人中文| 日韩av在线大香蕉| 精品福利观看| 欧美成人一区二区免费高清观看 | 欧美乱妇无乱码| 久久国产精品影院| 久久精品91无色码中文字幕| 亚洲欧美精品综合一区二区三区| 亚洲国产欧美一区二区综合| 一级a爱视频在线免费观看| 亚洲欧美一区二区三区黑人| 岛国在线观看网站| 99国产精品一区二区蜜桃av| 97人妻天天添夜夜摸| 两个人视频免费观看高清| 欧美 亚洲 国产 日韩一| 国产av精品麻豆| 91在线观看av| 国产伦人伦偷精品视频| 国产熟女xx| 久久中文字幕人妻熟女| 咕卡用的链子| 精品免费久久久久久久清纯| 757午夜福利合集在线观看| 久久精品国产亚洲av高清一级| 无限看片的www在线观看| 大陆偷拍与自拍| 精品国产乱子伦一区二区三区| 国产一卡二卡三卡精品| 久久久久久久久免费视频了| 老汉色av国产亚洲站长工具| 亚洲 欧美一区二区三区| www.999成人在线观看| 黄色a级毛片大全视频| 宅男免费午夜| 免费观看精品视频网站| 亚洲精品粉嫩美女一区| 黑人操中国人逼视频| 久久久久久久午夜电影| 亚洲国产日韩欧美精品在线观看 | 国产亚洲欧美98| 日韩欧美国产一区二区入口| 两人在一起打扑克的视频| 99国产极品粉嫩在线观看| 免费不卡黄色视频| 国产在线观看jvid| 夜夜夜夜夜久久久久| 日日夜夜操网爽| 成人国产一区最新在线观看| 大香蕉久久成人网| 91在线观看av| АⅤ资源中文在线天堂| 国产片内射在线| 中文字幕人成人乱码亚洲影| 国产国语露脸激情在线看| 日本黄色视频三级网站网址| 夜夜爽天天搞| 中文字幕av电影在线播放| 国产精品久久久人人做人人爽| 无人区码免费观看不卡| 欧美 亚洲 国产 日韩一| 国产黄a三级三级三级人| 成人亚洲精品av一区二区| 久久久久久免费高清国产稀缺| 日本免费a在线| 精品第一国产精品| 日韩成人在线观看一区二区三区| 18美女黄网站色大片免费观看| 免费一级毛片在线播放高清视频 | 性色av乱码一区二区三区2| 欧美一级a爱片免费观看看 | 精品高清国产在线一区| 久久天堂一区二区三区四区| 亚洲免费av在线视频| 免费在线观看影片大全网站| 亚洲aⅴ乱码一区二区在线播放 | 免费女性裸体啪啪无遮挡网站| 欧美一级a爱片免费观看看 | 黑人操中国人逼视频| 91精品三级在线观看| 午夜成年电影在线免费观看| 国产又爽黄色视频| 色尼玛亚洲综合影院| 宅男免费午夜| 窝窝影院91人妻| 国产精品爽爽va在线观看网站 | 亚洲中文字幕日韩| 一区在线观看完整版| 伊人久久大香线蕉亚洲五| 别揉我奶头~嗯~啊~动态视频| 亚洲av成人不卡在线观看播放网| e午夜精品久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 亚洲男人天堂网一区| 夜夜躁狠狠躁天天躁| 国产亚洲欧美在线一区二区| 亚洲专区字幕在线| av网站免费在线观看视频| 久热爱精品视频在线9| 桃红色精品国产亚洲av| ponron亚洲| 97超级碰碰碰精品色视频在线观看| 两性夫妻黄色片| 夜夜看夜夜爽夜夜摸| 后天国语完整版免费观看| 欧美一区二区精品小视频在线| 极品人妻少妇av视频| 国产精品久久久久久亚洲av鲁大| 如日韩欧美国产精品一区二区三区| 老熟妇乱子伦视频在线观看| 长腿黑丝高跟| 久久久国产成人精品二区| 成人av一区二区三区在线看| 麻豆av在线久日| 亚洲一区二区三区不卡视频| 一二三四社区在线视频社区8| 最近最新中文字幕大全电影3 | 一本综合久久免费| 在线免费观看的www视频| av福利片在线| 亚洲自偷自拍图片 自拍| 黑人欧美特级aaaaaa片| av欧美777| 久久精品国产99精品国产亚洲性色 | 人人妻人人澡人人看| 亚洲国产欧美一区二区综合| 又黄又爽又免费观看的视频| 亚洲精品中文字幕一二三四区| 在线十欧美十亚洲十日本专区| 久久人妻福利社区极品人妻图片| 丝袜美足系列| 一个人免费在线观看的高清视频| 成人欧美大片| 岛国视频午夜一区免费看| xxx96com| 国产欧美日韩一区二区三| 亚洲国产中文字幕在线视频| 日本欧美视频一区| av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 性欧美人与动物交配| 日韩 欧美 亚洲 中文字幕| 国产成人啪精品午夜网站| 欧美黑人精品巨大| 亚洲av电影不卡..在线观看| 两个人免费观看高清视频| 在线播放国产精品三级| 国产欧美日韩一区二区三| 美女扒开内裤让男人捅视频| 久久婷婷人人爽人人干人人爱 | 亚洲av成人不卡在线观看播放网| 69精品国产乱码久久久| 国产精品亚洲一级av第二区| 69av精品久久久久久| 老司机在亚洲福利影院| 国产精品 国内视频| 女同久久另类99精品国产91| 在线观看免费日韩欧美大片| 久久香蕉精品热| 男女下面插进去视频免费观看| 91老司机精品| 黄色a级毛片大全视频| 日韩成人在线观看一区二区三区| 制服诱惑二区| 国产视频一区二区在线看| 18禁美女被吸乳视频| 久久精品国产综合久久久| 国产黄a三级三级三级人| 校园春色视频在线观看| 国产精品久久久av美女十八| 丰满的人妻完整版| 久久午夜综合久久蜜桃| 黄片小视频在线播放| 一级毛片女人18水好多| 亚洲av日韩精品久久久久久密| 黑人操中国人逼视频| 香蕉丝袜av| 欧美激情久久久久久爽电影 | 精品第一国产精品| 免费高清视频大片| 丁香六月欧美| 国产精品自产拍在线观看55亚洲| 国产1区2区3区精品| 丝袜美腿诱惑在线| 老司机深夜福利视频在线观看| 久久精品国产亚洲av高清一级| 不卡av一区二区三区| 日日干狠狠操夜夜爽| 一边摸一边抽搐一进一出视频| 男女下面插进去视频免费观看| 亚洲专区国产一区二区| 中国美女看黄片| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| netflix在线观看网站| а√天堂www在线а√下载| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区激情短视频| 国产亚洲精品综合一区在线观看 | 日韩成人在线观看一区二区三区| 亚洲国产中文字幕在线视频| 一进一出好大好爽视频| 久久久久精品国产欧美久久久| 99国产精品免费福利视频| 可以在线观看毛片的网站| 亚洲中文av在线| 成人永久免费在线观看视频| 久久精品国产亚洲av高清一级| 亚洲国产精品合色在线| 日韩大码丰满熟妇| 丝袜在线中文字幕| 老司机午夜十八禁免费视频| 久久国产精品男人的天堂亚洲| www.熟女人妻精品国产| 亚洲精品美女久久久久99蜜臀| 麻豆一二三区av精品| 丝袜美腿诱惑在线| 激情在线观看视频在线高清| 国产成人影院久久av| 国产午夜福利久久久久久| 三级毛片av免费| 波多野结衣av一区二区av| 法律面前人人平等表现在哪些方面| 1024香蕉在线观看| 一级毛片女人18水好多| 大陆偷拍与自拍| 色老头精品视频在线观看| 电影成人av| 麻豆久久精品国产亚洲av| 一进一出抽搐动态| 精品久久久久久,| 天堂影院成人在线观看| 99国产精品免费福利视频| 亚洲国产日韩欧美精品在线观看 | 精品国产乱码久久久久久男人| 曰老女人黄片| 禁无遮挡网站| 午夜精品国产一区二区电影| 国内精品久久久久久久电影| 女人高潮潮喷娇喘18禁视频| 可以在线观看的亚洲视频| 999精品在线视频| 国产精品秋霞免费鲁丝片| 伊人久久大香线蕉亚洲五| 在线观看日韩欧美| 国产欧美日韩一区二区三区在线| 亚洲成国产人片在线观看| 久久婷婷人人爽人人干人人爱 | 成人免费观看视频高清| 搞女人的毛片| 禁无遮挡网站| 亚洲伊人色综图| 一卡2卡三卡四卡精品乱码亚洲| 他把我摸到了高潮在线观看| 亚洲欧美日韩高清在线视频| 午夜福利免费观看在线| 岛国在线观看网站| 午夜福利视频1000在线观看 | 国产麻豆成人av免费视频| 91成人精品电影| 国产精品久久电影中文字幕| 99热只有精品国产| netflix在线观看网站| 免费女性裸体啪啪无遮挡网站| 精品福利观看| 如日韩欧美国产精品一区二区三区| 好男人电影高清在线观看| 99久久精品国产亚洲精品| 欧美av亚洲av综合av国产av| 免费久久久久久久精品成人欧美视频| 宅男免费午夜| 精品久久久久久久人妻蜜臀av | 欧美日韩亚洲国产一区二区在线观看| 男女做爰动态图高潮gif福利片 | 午夜精品国产一区二区电影| 久久久久久亚洲精品国产蜜桃av| 黄色女人牲交| 中文字幕久久专区| 国产精品av久久久久免费| 激情在线观看视频在线高清| 亚洲色图 男人天堂 中文字幕|