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

    仿生飛行器非定常氣動優(yōu)化設計研究進展與挑戰(zhàn)

    2018-03-09 07:54:56肖天航羅東明鄭祥明昂海松吉愛紅南京航空航天大學江蘇南京210016
    空氣動力學學報 2018年1期
    關鍵詞:外形氣動飛行器

    肖天航, 羅東明, 鄭祥明, 昂海松, 吉愛紅(南京航空航天大學, 江蘇 南京 210016)

    0 引 言

    自然界的鳥、昆蟲等飛行生物具有高超的飛行技巧,翅膀的大幅撲動是其奧秘所在,是其克服尺寸小、速度低引起的低雷諾數(shù)空氣動力問題的重要方式。受自然界飛行生物啟發(fā),仿生撲翼在微型飛行器上獲得成功并重新引起學界的廣泛關注。通過實驗和數(shù)值模擬研究,國際、國內(nèi)學者在飛行原理、氣動機理、飛行力學以及與之相關的實驗和數(shù)值模擬方法等方面開展了大量的工作,揭示了多種撲翼高升力氣動機理,增進了對撲翼飛行原理的了解和掌握,為仿生飛行器的設計提供了原理依據(jù),相關研究進展和成果可參見Shyy[1]、孫茂[2]、Sane[3]、Wang[4]、楊文青等[5]的綜述文獻和著作。基于對撲翼空氣動力學原理的認識,近年來,多種形式的仿生飛行器相繼被研發(fā)問世。但航空界普遍認為,當前仿生飛行器的進一步發(fā)展和應用遇到了新的技術瓶頸。

    一方面,現(xiàn)階段仿生飛行器的設計仍是對飛行生物的簡單模仿,效率低,性能遠不及自然界飛行生物,在空氣動力方面亟待進一步精細設計優(yōu)化。受非定常氣動設計優(yōu)化理論與方法不足的制約,現(xiàn)有仿生飛行器的設計主要沿用簡單的估參數(shù)-預分析的單向設計思路;通過實驗或數(shù)值模擬,研究撲動運動參數(shù)或形狀外形參數(shù)對氣動的影響,以此獲得氣動特性隨設計變量的變化規(guī)律,并據(jù)此來確定仿生撲翼的設計參數(shù)(或大概的設計范圍),由于設計參數(shù)眾多且設計范圍較大,這種逐個分析參數(shù)的方法需要的計算量巨大;更嚴重的問題是,該方法常常割裂各個設計參數(shù)之間的內(nèi)在聯(lián)系,并不能實現(xiàn)真正的最佳設計。

    另一方面,仿生撲翼流場具有強烈的非定常效應、渦流動復雜、粘性效應強,且除氣動外形之外,翼的復雜大幅撲動運動也是其重要的氣動措施,兩者耦合使得設計變量眾多,屬大參數(shù)設計問題。撲翼氣動問題的特殊性給氣動建模和設計優(yōu)化理論與方法帶來了極大的挑戰(zhàn),常規(guī)定常氣動設計方法已不能適應非定常仿生飛行器這類新問題。仿生撲翼這一新型飛行器技術的發(fā)展迫切需要基礎的設計優(yōu)化理論方法有大的革新和新的突破。

    基于仿生飛行器的發(fā)展需要,仿生撲翼的非定常氣動設計優(yōu)化問題逐漸受到重視,近年,國內(nèi)外學者針對撲翼氣動優(yōu)化設計問題,從氣動建模、設計優(yōu)化理論方法和應用研究等角度開展工作,構建了不同保真度的優(yōu)化設計方法和模型用于撲翼氣動的優(yōu)化設計,取得了較大的進展;但鑒于撲翼非定常流動問題的復雜性和特殊性,仿生撲翼非定常氣動優(yōu)化設計仍面臨很大的困難和挑戰(zhàn)。本文擬從氣動建模、優(yōu)化方法策略等方面,對當前仿生飛行器的非定常氣動優(yōu)化設計研究取得的進展、面臨的困難和挑戰(zhàn)進行分析、評述和總結,以與同行交流。為方便起見,本文按照氣動模型的保真度和優(yōu)化方法策略的不同,將現(xiàn)有仿生撲翼氣動優(yōu)化設計的研究工作歸納為基于傳統(tǒng)方法和基于現(xiàn)代新型方法兩個大類進行論述。

    1 運用傳統(tǒng)方法的仿生撲翼氣動優(yōu)化設計

    這一類的仿生撲翼氣動優(yōu)化設計研究,按照氣動模型的保真度從低到高歸納為基于葉素理論氣動建模、基于快速數(shù)值方法氣動建模、基于非定常N-S方程氣動建模和實驗測量幾種。

    1.1 基于葉素理論氣動建模的仿生撲翼氣動優(yōu)化設計

    (1)

    (2)

    圖1 基于葉素理論的撲翼氣動建模示意圖[8]Fig.1 Schematics of aerodynamic model of flapping wing based on blade element theory[8]

    Berman和Wang等[7]即基于這類準定常模型,運用遺傳全局優(yōu)化和Nelder-Mead單純形的梯度局部優(yōu)化相結合的組合算法,對昆蟲撲翼懸停狀態(tài)的11個撲動運動參數(shù)進行優(yōu)化設計。以周期平均升力等于重力為約束、以平均功率消耗(氣動力+慣性力)最小為目標,先用遺傳算法進行大范圍尋優(yōu),后用梯度優(yōu)化算法進行局部優(yōu)化,最終優(yōu)化設計得到的撲翼運動軌跡與對自然界生物的觀測結果相似。Kurdi等[9]運用同樣的氣動模型,采用梯度優(yōu)化方法對昆蟲翼懸停狀態(tài)的30個撲動運動變量進行優(yōu)化,撲動運動規(guī)律采用樣條函數(shù)擬合使其光滑連續(xù),優(yōu)化設計模型甚至考慮了彈性影響。 Ke 和 Zhang等[10]同樣基于準定常模型,運用遺傳算法和梯度優(yōu)化方法,以撲翼的平面幾何外形和撲動運動參數(shù)為設計變量,以功耗最小為目標進行優(yōu)化設計,分析認為撲翼外形參數(shù)和運動參數(shù)之間通過展弦比和雷諾數(shù)建立有較強的耦合關系而共同影響著目標函數(shù)。Yan和Taha等[11]對懸停狀態(tài)的撲翼運動參數(shù)進行優(yōu)化設計以使功耗最小,優(yōu)化過程采用了幾種基于葉素理論建立的精度不一的氣動模型,包括上述準定常模型和Haithen等[12]建立的考慮前緣渦等非定常效應的全非定常氣動模型;研究認為準定常模型較高估計了升力系數(shù),使得最小功耗值比全非定常模型優(yōu)化結果更小,論文通過綜合比較建議撲翼的非定常氣動優(yōu)化設計應該采用較高精度的非定常氣動模型。

    Gogulapati等[13-14]完成了懸停狀態(tài)柔性撲翼的優(yōu)化設計,氣動模型采用基于葉素理論的結合勢流求解器的非定常氣動模型[15],即葉素模型的微元氣動力通過求解二維速勢方程得到,速勢方程的求解假定分離從前緣開始,渦強度由強化的前緣滯止條件和后緣庫塔條件決定,同時計及前緣渦高升力;優(yōu)化設計也包括了非線性結構有限元模型;由于氣動模型和結構有限元模型計算量相對較大,作者通過Kriging插值建立時均推力和時均功率的代理模型;設計變量為撲動和俯仰運動參數(shù)以及柔性翼的剛度共13個設計變量;優(yōu)化算法是全局優(yōu)化和梯度優(yōu)化的組合算法;優(yōu)化結果表明在一定撲動幅度和頻率下,撲動與俯仰運動的相位差對效率影響很大,柔性翼的效率峰值出現(xiàn)于較小俯仰運動角度情況,在若干設計變量區(qū)間,柔性翼相比剛性翼并不占優(yōu)。

    基于葉素理論氣動模型的優(yōu)化設計可快速地在全局范圍內(nèi)尋得最優(yōu)解,但氣動模型本身的問題決定了優(yōu)化設計的精度不會很高;原因是,葉素理論模型從原理上就不考慮展向流動和翼尖渦,而實際上由于撲翼展弦比一般不大,且強烈的展向流動是撲翼延遲失速的重要機理;因此,從設計精度角度來說,這類優(yōu)化設計只適合初步概念設計階段大量參數(shù)的粗略分析和篩選。

    1.2 基于快速數(shù)值方法的仿生撲翼氣動優(yōu)化設計

    撲翼非定常氣動的快速數(shù)值方法主要是非定常渦格法(Unsteady Vortex-Lattice Method, UVLM)。UVLM在運動的升力面上進行渦格劃分,并在每個渦格上布置馬蹄渦和控制點,通過在控制點處滿足升力面法向無穿透的邊界條件來進行求解。相對葉素理論氣動模型,該方法具有較高的精度,同時計算量又遠小于N-S方程的求解,在撲翼非定常氣動力計算分析上得到較多應用,同時也在不斷改進發(fā)展中[16-17]。

    目前有少數(shù)學者運用UVLM進行氣動建模用于仿生撲翼的氣動優(yōu)化設計。Stanford等[18]將UVLM氣動模型嵌入梯度優(yōu)化器,在升力、阻力約束條件下以氣動功率最小為目標,對前飛狀態(tài)的主動變形仿鳥撲翼的弦向扭轉和展向彎曲進行優(yōu)化,與Kurdi[9]類似,運動規(guī)律采用樣條函數(shù)擬合以保證其光滑連續(xù);優(yōu)化結果表明,在低縮減頻率下,展向彎曲變形對增加時均推力、提高效率起主導作用,在較高縮減頻率下,弦向扭轉變形則在減少功耗上更為重要,該研究也表明增加設計變量數(shù)量可獲得更優(yōu)的設計和更高的計算效率。 Ghommen等[19]基于非定常渦格法結合確定全局優(yōu)化方法,建立基于三角函數(shù)和樣條函數(shù)的兩種撲翼運動規(guī)律的參數(shù)化模型,對主動變形仿鳥撲翼的運動規(guī)律進行優(yōu)化設計,設計參數(shù)僅4-8個;結果表明樣條函數(shù)的參數(shù)化建模由于三角函數(shù)模型,翼變形對效率有益,且扭轉和彎曲中的二階模態(tài)對效率提升有明顯作用。之后,Ghommen等[19]又對撲翼外形(包括展弦比、彎度、前緣后緣彎曲形狀等)進行基于全局優(yōu)化和梯度優(yōu)化結合的優(yōu)化設計,表明合理配置展弦比、彎度和前后緣外形,對改善尾渦結構(如圖2所示)、提高推進效率有積極意義,其結果可為仿鳥撲翼的工程設計提供有價值的參考。Stewart等[21]將UVLM和平板有限元模型結合,運用局部優(yōu)化方法,對柔性撲翼的外形和結構進行了多目標的優(yōu)化設計,結果顯示增強結構彈性變形,有利于提高推力但也增加了功率消耗。賀紅林等[22]運用非定常渦格法建立氣動模型,結合模式搜索算法,以升力和推力最大化為目標,對撲動和俯仰運動相位差、撲動頻率和柔性扭轉角進行優(yōu)化設計,對三者的配匹給出了一定的建議。

    圖2 優(yōu)化設計展弦比及前后緣外形對撲翼尾渦的形狀和強度影響[19]Fig.2 Optimal shape on the wake patterns and vorticity strength of a bird-like flapping wing[19]

    總的來說,基于ULVM氣動模型的撲翼優(yōu)化設計,計算量適中,用于仿中大型鳥類撲翼這類較弱或中等強度非定常氣動設計,也能取得較高的精度;但用于非定常效應更嚴重的仿昆蟲撲翼的優(yōu)化設計,則還需ULVM方法的進一步發(fā)展和完善。

    1.3 基于非定常N-S方程和實驗測量的仿生撲翼氣動優(yōu)化設計

    對非定常效應強烈、渦流動復雜、粘性效應強的仿生撲翼流場,基于非定常N-S方程的氣動模型能獲得相對較高的保真度,但N-S方程求解復雜,動態(tài)網(wǎng)格處理繁瑣,計算量大,若按照傳統(tǒng)思路,采用全局優(yōu)化或常規(guī)梯度優(yōu)化算法,對撲翼尤其三維復雜撲翼的氣動優(yōu)化設計將由于超大的計算量而難以實現(xiàn)。目前只有極少數(shù)工作進行了嘗試,但也僅限于二維或簡單三維情況,且設計變量數(shù)目非常有限。

    Tuncer等[23]基于前向模式差分的梯度優(yōu)化算法,建立二維N-S方程的優(yōu)化設計模型,對在自由來流中做俯仰和沉浮運動的翼型進行推力和推進效率最大化的優(yōu)化設計;揭示了撲動運動與流場結構、進而與氣動效率的因果關系,表明能產(chǎn)生前緣渦的撲動運動可增加升力但降低了推進效率,而不產(chǎn)生前緣渦的撲動運動則有利于提高推進效率。Soueid等[24]通過建立基于敏感度方程和二維不可壓N-S方程氣動模型的優(yōu)化設計模型,完成了NACA0012翼型撲動運動規(guī)律(沉浮、俯仰和相位差)的設計。Culbreth等[25]運用有限差分的梯度優(yōu)化方法和三維N-S方程氣動模型,分析了矩形撲翼的4個展向扭轉參數(shù)對推進效率的敏感度。

    最近也有個別學者嘗試通過實驗測量或基于有限的實驗數(shù)據(jù),利用相關的優(yōu)化算法,以撲翼氣動性能參數(shù)為目標進行撲動運動參數(shù)的優(yōu)化,Milano[26]、Chaudhuri[27]、Nan[28]等的工作是其中的主要代表。但基于實驗測量的優(yōu)化設計,成本高 、難度大,不易掌控,應用于工程實際問題存在比較大的困難。

    1.4 傳統(tǒng)優(yōu)化設計面臨的問題

    雖然仿生飛行器的發(fā)展備受關注,但目前對仿生撲翼的非定常氣動設計優(yōu)化方法仍缺乏系統(tǒng)的研究?,F(xiàn)有研究大多數(shù)都沿用傳統(tǒng)的優(yōu)化設計理論,利用遺傳算法等全局尋優(yōu)技術或利用常規(guī)的基于梯度的局部尋優(yōu)方法,構建仿生撲翼的優(yōu)化設計模型,這一思路在初步設計或少量設計參數(shù)問題上可發(fā)揮一定的作用,但用于非定常仿生飛行器會面臨很大的問題和挑戰(zhàn)。正如前文所述,仿生撲翼流動具有強非定常、低雷諾數(shù)特點,渦流動復雜,粘性效應強,其氣動的計算分析需采用基于N-S方程的高保真度數(shù)值方法,計算量大;另一方面,仿生撲翼的設計變量,包括氣動外形參數(shù)和撲動運動參數(shù),數(shù)量可達數(shù)十甚至上百個。傳統(tǒng)優(yōu)化設計方法的計算量與設計變量數(shù)量相關,因此,用于仿生飛行器問題,要么以犧牲精度為代價,氣動分析采用較低精度的模型,如準定常氣動模型、渦格法等;或者氣動模型采用高階模型,但嚴格限制設計變量數(shù)量,即使采用代理模型技術構建響應面,其設計參數(shù)數(shù)量也相當有限??梢哉f,基于傳統(tǒng)優(yōu)化設計理論的思路很難實現(xiàn)有效的高精度的優(yōu)化設計,對仿生飛行器問題,要另辟蹊徑,需在適于大參數(shù)非定常氣動問題的優(yōu)化設計理論和方法研究上尋求突破。

    2 仿生撲翼氣動的非定常伴隨優(yōu)化設計

    近年來,航空界新發(fā)展一種基于控制理論的伴隨優(yōu)化方法,代表著當前國際先進的設計優(yōu)化理論和發(fā)展趨勢。該方法將設計問題當作最優(yōu)控制問題,通過引入并求解一個伴隨方程系統(tǒng)實現(xiàn)設計目標對設計變量的靈敏度分析,一輪優(yōu)化,只需求解一次流動控制方法和一次伴隨方程,其計算量與設計變量數(shù)量無關,成為解決大參數(shù)氣動設計優(yōu)化問題技術瓶頸的重要措施和發(fā)展方向。目前,對伴隨優(yōu)化理論和方法的研究已經(jīng)在多個領域的氣動外形設計優(yōu)化上(如葉輪機械[29]、翼型和三維機翼[30]、氣動噪聲[31]等)取得了很多成果,國內(nèi)的北京大學劉峰團隊[32-33]、西安交通大學豐鎮(zhèn)平團隊[34-35]、西北工業(yè)大學[36-38]、上海交大[31]以及北京理工大學[29,39]、中國空氣動力研究與發(fā)展中心[40-41]、蘭州理工[42]和南京航空航天大學[43-45]等,在該領域的研究都相當出色。但需要指出的是,上述研究的重點主要針對定常流動的氣動外形設計優(yōu)化。

    隨著非定常流動及其設計優(yōu)化問題日益受到重視,為其發(fā)展有效的優(yōu)化設計理論與方法也迫在眉睫。對非定常問題,時均性能參數(shù)通常是受關心的量,理論上來說,非定常流場的優(yōu)化設計問題,也可以看成是受非定常流體控制方程和動網(wǎng)格方程約束的時均目標函數(shù)的最小值優(yōu)化問題,同樣可以運用最優(yōu)控制理論求解。對某段物理時間,通過建立與時均目標函數(shù)和非定常流動控制方程求解相容的非定常伴隨方程,在前向積分求解出流場后,再逆時間積分求出各物理時間步的伴隨變量,從而實現(xiàn)該段時間內(nèi)時均目標函數(shù)對設計變量的敏感度分析,這一思路為解決大參數(shù)非定常設計優(yōu)化問題提供了可能,在仿生撲翼非定常氣動問題上也具有可行性。

    對仿生飛行器而言,時均氣動參數(shù)是衡量性能的重要指標。我們將仿生飛行器的非定常氣動優(yōu)化設計問題,當成是受非定常流體控制方程和動態(tài)網(wǎng)格方程約束的時均目標函數(shù)最小值優(yōu)化問題,即:

    (3)

    基于最優(yōu)控制理論,可找到適當?shù)腝和D使得受約束的目標函數(shù)值最小,達到解決式(3)最小值優(yōu)化問題的目的。通過拉格朗日乘子Λf,Λg將式(3)中的約束引入到目標函數(shù)中,形成新的拉格朗日函數(shù)L,即:

    (4)

    目前已有少數(shù)學者開展了這方面的研究。如Nielsen[46-47]、Mavriplis[48]、Alonso[49]等在非定常離散伴隨和非定常連續(xù)伴隨優(yōu)化理論及求解方法上開展了開拓性的研究工作,并在旋翼、動態(tài)失速翼型等的氣動外形設計上進行了初步應用研究;北京理工大學雷娟棉[50-51]發(fā)展了非定常離散伴隨方法并在二維跨聲速俯仰翼型和低雷諾數(shù)翼型的外形優(yōu)化上進行了驗證;上海應用數(shù)學力學研究所[52]研究了二維沉浮翼型的連續(xù)伴隨優(yōu)化方法;清華大學[53]初步探索了基于非定常流動控制方程瞬態(tài)解和非定常伴隨方程瞬態(tài)解的動態(tài)演化伴隨方法,以達到提高定常氣動外形優(yōu)化效率的目的,其思路和優(yōu)勢在二維翼型外形設計算例上得到驗證。

    圖3 非定常氣動離散伴隨優(yōu)化設計思路和流程Fig.3 Flow chart of unsteady discrete adjoint optimization for flapping wings

    對仿生撲翼的非定常氣動問題,伴隨優(yōu)化方法也表現(xiàn)出良好的應用前景,近兩年逐漸受到關注。Lee和Liou[54-55]對伴隨優(yōu)化在仿生撲翼問題上的應用進行了初步嘗試,分別研究了二維撲動翼型的外形及簡單運動軌跡的優(yōu)化設計,由于伴隨方法的計算量與設計變量數(shù)目無關,該方法才有可能將高精度N-S方程的氣動模型用于較多設計變量的非定常氣動優(yōu)化。Jones等[56]則將應用對象擴展到三維情況,分別對仿昆蟲撲翼的撲動運動參數(shù)、幾何外形參數(shù)進行了優(yōu)化設計,并嘗試了將運動參數(shù)和外形參數(shù)聯(lián)立,構成統(tǒng)一的設計空間進行設計,其結果也表明兩者對優(yōu)化目標而言存在較強的非線性耦合關系。Xu等[57]針對撲翼問題的特殊性,發(fā)展了一種新的伴隨方法中動態(tài)邊界的處理辦法,并將其應用到剛性、柔性平板撲翼的優(yōu)化設計,表明柔性對推力特性和推進效率都有幫助;對三維矩形撲翼懸停狀態(tài)的運動參數(shù)和柔性變形參數(shù)的優(yōu)化設計結果證實,非定常伴隨優(yōu)化方法能找到最佳的運動參數(shù)和柔性變形,使得前緣渦更大,下洗流動更強,柔性變形能更好地捕捉尾渦,而這些特征是低階氣動模型無法做到的,顯示了伴隨優(yōu)化在仿生飛行器上的良好愿景。就目前而言,仿生撲翼的非定常離散伴隨優(yōu)化設計的研究不多,但已有工作發(fā)揮了良好的推動作用。

    3 仿生飛行器非定常氣動優(yōu)化設計的問題與展望

    仿生撲翼流場具有強烈的非定常效應,流動機理復雜多樣,雷諾數(shù)低、粘性效應強;同時,氣動外形和翼的復雜大幅撲動運動都是其重要的氣動措施,兩者相互耦合,設計變量數(shù)目眾多;仿生撲翼的復雜性和特殊性給仿生飛行器的氣動優(yōu)化設計帶來很大的困難和挑戰(zhàn)。

    針對仿生飛行器發(fā)展的需要,現(xiàn)階段已經(jīng)從不同途徑,或基于傳統(tǒng)方法,或發(fā)展新的技術,發(fā)展了若干不同精度的氣動模型或氣動評估方法,運用不同優(yōu)化設計策略進行了仿生飛行器氣動的優(yōu)化設計。這些氣動模型或方法各有優(yōu)勢,從實際工程來說,不必片面強調高精度而忽略其他。對于仿生飛行器的具體設計,從概念設計到詳細設計,不同階段需要不同精度的氣動優(yōu)化模型和策略,發(fā)展并應用分層次、分階段的多保真度仿生撲翼非定常氣動優(yōu)化設計方法和模型,將是仿生飛行器領域值得考慮的問題。

    對仿生飛行器非定常氣動的優(yōu)化設計,未來以下幾個方面的問題需要重點關注和尋求解決。

    1) 仿生撲翼的低階氣動模型可在仿生飛行器的概念設計階段發(fā)揮重要作用,未來應結合撲翼非定常氣動機理、相關空氣動力學原理和理論,發(fā)展和完善計算模型,充分計及前緣渦、延遲失速、旋轉效應等高升機制和展向流動、翼尖渦等在氣動力產(chǎn)生方面的作用或影響,形成一套效率好、精度高的仿生飛行器非定常氣動快速優(yōu)化設計的技術方案。

    2) 基于傳統(tǒng)方法的優(yōu)化設計因為模型精度與計算量的矛盾而不能滿足仿生飛行器這一復雜問題的精細設計需要,伴隨優(yōu)化是解決這一問題的有效途徑,但非定常伴隨優(yōu)化理論和方法的研究剛剛起步,在仿生飛行器上應用還不成熟,在非定常伴隨方程構造、高效求解、初邊值條件(空間的和時間的)以及與仿生飛行器特殊問題相關的適定性、目標泛函選取原則、設計準則等方面,其理論和方法的發(fā)展完善還需要學者持續(xù)深入研究;此外,仿生撲翼的大幅撲動運動常常要求流場解算基于動態(tài)嵌套網(wǎng)格進行,在動態(tài)嵌套網(wǎng)格上,非定常伴隨方程的構造求解理論與方法及其擴展遷移也是仿生飛行器非定常伴隨優(yōu)化研究需要解決的問題。

    3) 就仿生撲翼氣動優(yōu)化設計本身來說,現(xiàn)有研究,無論是基于傳統(tǒng)優(yōu)化方法還是應用新型伴隨優(yōu)化方法的,大多數(shù)將氣動外形和撲動運動參數(shù)分開,單一地就某一方面進行設計,較少將運動或軌跡參數(shù)與外形參數(shù)耦合構建一個統(tǒng)一的設計空間;而對仿生飛行器而言,外形與撲動運動的耦合又十分重要,如何將復雜非連續(xù)撲動運動(如翼翻轉、往返非對稱等)參數(shù)化并與外形參數(shù)聯(lián)合,構造連續(xù)光順的利于優(yōu)化設計收斂的設計空間,是仿生飛行器非定常氣動優(yōu)化設計值得探討的問題。

    4) 對仿生撲翼的優(yōu)化設計,目前多數(shù)研究工作也主要集中在懸停狀態(tài),前飛狀態(tài)、機動飛行狀態(tài)甚至陣風狀態(tài)下的撲翼氣動優(yōu)化設計也應該是下一步要解決的問題;此外,翼的柔性是影響氣動效率的重要因素,翼的結構柔性也應是其氣動優(yōu)化設計的變量,柔性撲翼的氣動建模及優(yōu)化設計、甚至柔性撲翼氣動-結構多學科耦合的優(yōu)化設計,都將是仿生飛行器領域的重要課題。

    總的來說,仿生飛行器非定??諝鈩恿Φ膬?yōu)化設計仍是一個具有挑戰(zhàn)性的課題,需要從基礎理論到工程技術研究的全面推動和發(fā)展。

    4 結 論

    仿生撲翼是一種新型的飛行器技術,由于其撲翼流場強烈的非定常特征和低雷諾數(shù)效應、外形與撲動運動耦合的大量設計變量,仿生飛行器的氣動優(yōu)化設計存在很大的難度,是一個非常復雜的不同于定常氣動問題的研究課題。由于缺乏成熟有效的設計方法和工具,目前仿生撲翼的優(yōu)化設計大多數(shù)按照傳統(tǒng)的優(yōu)化設計思路,存在精度低或計算量大的缺陷,較難取得良好的優(yōu)化設計效果。伴隨優(yōu)化方法是仿生飛行器非定常氣動高精度高效率優(yōu)化設計的有效途徑,但該方法仍有待完善發(fā)展。對仿生撲翼氣動的優(yōu)化設計,需耦合外形參數(shù)和撲翼運動參數(shù),建立聯(lián)合的設計空間;前飛狀態(tài)、機動飛行狀態(tài)或陣風狀態(tài)下的撲翼氣動優(yōu)化設計也是未來工作關注的重點。對仿生飛行器設計來說,采用不同保真度模型,分層次、分階段地進行氣動優(yōu)化設計是從概念設計走向詳細設計的有效策略。

    [1]Shyy W, Aono H, Kang C K, et al. An introduction to flapping wing aerodynamics[M]. Cambridge University Press, 2013

    [2]孫茂. 昆蟲飛行的空氣動力學[J]. 力學進展, 2015, 45: 1-28

    [3]Sane S P. The aerodynamics of insect flight[J]. Journal of Experimental Biology, 2003, 206: 4191-4208

    [4]Wang Z J. Dissecting insect flight[J]. Annual Review of Fluid Mechanics, 2005, 37: 183-210

    [5]楊文青, 宋筆鋒, 宋文萍. 仿生微型撲翼飛行器中的空氣動力學問題研究進展與挑戰(zhàn)[J]. 實驗流體力學, 2015, 29: 1-10

    [6]Sane S P, Dickinson M H. The aerodynamic effects of wing rotation and a revised quasi-steady model of flapping flight[J]. The Journal of Experimental Biology, 2002, 205: 1087-1096

    [7]Berman G J, Wang Z J. Energy-minimizing kinematics in hovering insect flight[J]. Journal of Fluid Mechanics, 2007, 582: 153

    [8]Whitney J P, Wood R J. Conceptual design of flapping-wing micro air vehicles[J]. Bioinspiration & Biomimetics, 2012, 7: 36001

    [9]Kurdi M, Stanford B, Beran P. Optimal kinematics of hovering insect flight for minimum mechanical power[C]//48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Reston, Virigina: American Institute of Aeronautics and Astronautics, 2010

    [10]Ke X, Zhang W, Cai X, et al. Wing geometry and kinematic parameters optimization of flapping wing hovering flight for minimum energy[J]. Aerospace Science and Technology, 2017, 64: 192-203

    [11]Yan Z, Taha H E, Hajj M R. Effects of aerodynamic modeling on the optimal wing kinematics for hovering MAVs[J]. Aerospace Science and Technology, 2015, 45: 39-49

    [12]Haithem E Tahaa, Muhammad R Hajja, Philip S Beranb. State-space representation of the unsteady aerodynamics of flapping flight[J]. Aerospace Science and Technology, 2014, 34: 1-11

    [13]Gogulapati A, Friedmann P P, Martins J R R A. Optimization of flexible flapping-wing kinematics in hover[J]. AIAA Journal, 2014, 52: 2342-2354

    [14]Gogulapati A, Friedmann P P, Martins J R R. Optimization of the kinematics of a flapping wing MAV in hover for enhanced performance[C]//54th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Boston, Massachusetts: 2013

    [15]Gogulapati A, Friedmann P P, Kheng E, et al. Approximate aeroelastic modeling of flapping wings in hover[J]. AIAA Journal, 2013, 51: 567-583

    [16]Gardiner J, Razak N A, Dimitriadis G, et al. Simulation of bird wing flapping using the unsteady vortex lattice method[C]//IFASD 2013-International Forum on Aeroelasticity and Structural Dynamics, 24-27 June, Bristol, 2013: 1-14

    [17]Nguyen A T, Kim J K, Han J S, et al. Extended unsteady vortex-lattice method for insect flapping wings[J]. Journal of Aircraft, 2016, 53(6): 1709-1718

    [18]Stanford B K, Beran P S. Analytical sensitivity analysis of an unsteady vortex-lattice method for flapping-wing optimization[J]. Journal of Aircraft, 2010, 47: 647-662

    [19]Ghommem M, Collier N, Niemi A H, et al. On the shape optimization of flapping wings and their performance analysis[J]. Aerospace Science and Technology, 2014, 32: 274-292

    [20]Ghommem M, Hajj M R, Mook D T, et al. Global optimization of actively morphing flapping wings[J]. Journal of Fluids and Structures, 2012, 33: 210-228

    [21]Stewart E C, Patil M J, Canfield R A, et al. Aeroelastic shape optimization of a flapping wing[J]. Journal of Aircraft, 2016, 53: 636-650

    [22]賀紅林, 周翔, 朱保利, 等. 基于非定常渦格法的撲翼飛行器氣動特性優(yōu)化[J]. 空氣動力學學報, 2012, 30(1): 74-79

    [23]Tuncer I, Kaya M. Optimization of flapping airfoils for maximum thrust and propulsive efficiency[J]. AIAA Journal, 2005, 43: 2329-2336

    [24]Soueid H, Guglielmini L, Airiau C, et al. Optimization of the motion of a flapping airfoil using sensitivity functions[J]. Computers and Fluids, 2009, 38: 861-874

    [25]Culbreth M, Allaneau Y, Jameson A. High-fidelity optimization of flapping airfoils and wings[C]//29th AIAA Applied Aerodynamics Conference, 2011

    [26]Milano M, Gharib M. Uncovering the physics of flapping flat plates with artificial evolution[J]. Journal of Fluid Mechanics, 2005, 534: 403-409

    [27]Chaudhuri A, Haftka R T, Ifju P G, et al. Experimental flapping wing optimization and uncertainty quantification using limited samples[J]. Structural and Multidisciplinary Optimization, 2014, 51: 957-970

    [29]季路成, 李偉偉, 伊衛(wèi)林. 伴隨方法用于葉輪機優(yōu)化設計的回顧與展望[J]. 航空發(fā)動機, 2011, 37(5): 53-57+62

    [30]Osusky L, Buckley H, Reist T, et al. Drag minimization based on the navier-stokes equations using a Newton-Krylov approach[J]. AIAA Journal, 2015, 53: 1555-1577

    [31]邱昇. 基于伴隨方法的渦扇發(fā)動機涵道的氣動噪聲優(yōu)化設計研究[D]. 上海交通大學博士學位論文, 2013

    [32]Luo J, Zhou C, Liu F. Multipoint design optimization of a transonic compressor blade by using an adjoint method[J]. Journal of Turbomachinery, 2013, 136: 51005-51005-10

    [33]Luo J, Xiong J, Liu F. Aerodynamic design optimization by using a continuous adjoint method[J]. Science China-Physics, Mechanics & Astronomy, 2014, 57: 1363-1375

    [34]張朝磊, 盧娟, 豐鎮(zhèn)平. 基于離散伴隨方法的透平葉柵反設計[J]. 工程熱物理學報, 2012, 33(4): 583-586

    [35]盧娟, 張朝磊, 豐鎮(zhèn)平. 離散伴隨氣動優(yōu)化設計方法的紊流擴展[J]. 西安交通大學學報, 2013, 47(1): 37-42

    [36]王龍, 宋文萍, 楊旭東. 一種基于伴隨方程迎角不變的多目標優(yōu)化減阻方法[J]. 西北工業(yè)大學學報, 2012, 30(1): 68-72

    [37]白俊強, 陳頌, 華俊. 基于伴隨方程和自由變形技術的跨聲速機翼氣動設計方法研究[J]. 空氣動力學學報, 2014, 32(6): 820-826+833

    [38]左英桃, 蘇偉, 高正紅. 基于離散共軛方法的高超聲速導彈氣動外形優(yōu)化設計[J]. 計算力學學報, 2012, 29(2): 284-289

    [39]李偉偉, 季路成, 伊衛(wèi)林. 基于伴隨方法的多級葉輪機三維葉片優(yōu)化設計[J]. 工程熱物理學報, 2014, 35(11): 2164-2167

    [40]李彬, 鄧有奇, 唐靜, 等. 基于三維非結構混合網(wǎng)格的離散伴隨優(yōu)化方法[J]. 航空學報, 2014, 35(3): 674-686

    [41]吳文華, 范召林, 陳德華, 等. 基于伴隨算子的大飛機氣動布局精細優(yōu)化設計[J]. 空氣動力學學報, 2012, 30(6): 719-724+760

    [42]張人會, 劉宜, 郭苗, 等. 基于伴隨方法的離心泵葉輪優(yōu)化設計[J]. 排灌機械工程學報, 2014, 32(11): 943-947+954

    [43]徐兆可, 夏健, 高宜勝. 基于三維非結構網(wǎng)格的連續(xù)伴隨優(yōu)化方法[J]. 南京航空航天大學學報, 2015, 47(1): 145-152

    [44]高宜勝, 伍貽兆, 夏健. 基于非結構網(wǎng)格離散型伴隨方法的翼型優(yōu)化[J]. 空氣動力學學報, 2013, 31(2): 244-265

    [45]Tang Z, Dong J. Couplings in multi-criterion aerodynamic optimization problems using adjoint methods and game strategies[J]. Chinese Journal of Aeronautics, 2009, 22: 1-8

    [46]Nielsen E J, Diskin B, Yamaleev N K. Discrete adjoint-based

    design optimization of unsteady turbulent flows on dynamic unstructured grids[J]. AIAA Journal, 2010, 48: 1195-1206

    [47]Nielsen E J, Diskin B, Yamaleev N K. Discrete adjoint-based design optimization of unsteady turbulent flows on dynamic unstructured grids[J]. AIAA Journal, 2013, 51: 1355-1373

    [48]Mani K, Mavriplis D. Geometry optimization in three-dimensional unsteady flow problems using the discrete adjoint[C]//51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Grapevine (Dallas/Ft. Worth Region), Texas: 2013

    [49]Economon T D, Palacios F, Alonso J J. An unsteady continuous adjoint approach for aerodynamic design on dynamic meshes[J]. AIAA Journal, 2015, 53: 2437-2453

    [50]何建東, 雷娟棉. 基于伴隨方法的動態(tài)非定常氣動外形優(yōu)化設計[J]. 北京理工大學學報, 2015, 35(2): 127-132

    [51]Lei J, He J. Adjoint-based aerodynamic shape optimization for low reynolds number airfoils[J]. Journal of Fluids Engineering, Transactions of the ASME, 2016, 138: 21401-21401-6

    [52]Fang L, Li X. Design optimization of unsteady airfoils with continuous adjoint method[J]. Applied Mathematics and Mechanics, 2015, 36: 1329-1336

    [53]高瑩瑩, 何樹沈, 孟育. 非定常動態(tài)演化伴隨優(yōu)化設計方法[J]. 物理學報, 2012, 61(20): 200206-1-200206-11

    [54]Lee B J, Padulo M, Liou M, et al. Non-sinusoidal trajectory optimization of flapping airfoil using unsteady adjoint approach[C]//49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, Florida: 2011

    [55]Lee B J, Liou M S. Unsteady adjoint approach for design optimization of flapping airfoils[J]. AIAA Journal, 2012, 50: 2460-2475

    [56]Jones M, Yamaleev N K. Adjoint-based optimization of flapping-wing flows[J]. AIAA Journal, 2015, 53: 934-947

    [57]Xu M, Wei M. Using adjoint-based optimization to study kinematics and deformation of flapping wings[J]. Journal of Fluid Mechanics, 2016, 799: 56-99.

    猜你喜歡
    外形氣動飛行器
    中寰氣動執(zhí)行機構
    高超聲速飛行器
    比外形,都不同
    基于NACA0030的波紋狀翼型氣動特性探索
    基于反饋線性化的RLV氣動控制一體化設計
    測控技術(2018年9期)2018-11-25 07:44:24
    復雜飛行器的容錯控制
    電子制作(2018年2期)2018-04-18 07:13:25
    論袁牧之“外形的演技”
    足趾移植再造手指術后外形的整形
    神秘的飛行器
    KJH101-127型氣動司控道岔的改造
    又大又爽又粗| 日韩免费av在线播放| 美女被艹到高潮喷水动态| 别揉我奶头~嗯~啊~动态视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美日韩无卡精品| 黄色 视频免费看| 少妇的逼水好多| 99久国产av精品| 少妇裸体淫交视频免费看高清| 美女大奶头视频| 国产伦一二天堂av在线观看| 又爽又黄无遮挡网站| 一本一本综合久久| 成在线人永久免费视频| 亚洲 国产 在线| 国产精品1区2区在线观看.| 国产成人精品久久二区二区免费| 国产精品 欧美亚洲| aaaaa片日本免费| 色精品久久人妻99蜜桃| 亚洲精品粉嫩美女一区| 99久久精品热视频| 午夜福利在线在线| 成年女人看的毛片在线观看| 亚洲人成网站在线播放欧美日韩| 久久久久久久久中文| 中文字幕久久专区| 最新美女视频免费是黄的| 免费在线观看成人毛片| 天天添夜夜摸| 婷婷六月久久综合丁香| 男人舔女人的私密视频| av天堂中文字幕网| 小蜜桃在线观看免费完整版高清| 夜夜躁狠狠躁天天躁| 日韩欧美国产在线观看| 很黄的视频免费| cao死你这个sao货| 国产av不卡久久| 亚洲,欧美精品.| 九九久久精品国产亚洲av麻豆 | 亚洲乱码一区二区免费版| 色哟哟哟哟哟哟| 日韩国内少妇激情av| 小蜜桃在线观看免费完整版高清| 国产探花在线观看一区二区| 亚洲av电影不卡..在线观看| 精品无人区乱码1区二区| 亚洲色图 男人天堂 中文字幕| 母亲3免费完整高清在线观看| 一个人免费在线观看电影 | 国产综合懂色| 狂野欧美激情性xxxx| 日本一二三区视频观看| 少妇裸体淫交视频免费看高清| 九九热线精品视视频播放| 国产欧美日韩一区二区三| 日韩人妻高清精品专区| 99视频精品全部免费 在线 | 久久精品综合一区二区三区| 国产成年人精品一区二区| 精品国产美女av久久久久小说| 午夜a级毛片| 黄色女人牲交| ponron亚洲| 男人舔女人的私密视频| 女同久久另类99精品国产91| 老司机在亚洲福利影院| 亚洲自偷自拍图片 自拍| 中文字幕高清在线视频| 精品国产乱子伦一区二区三区| 国产av一区在线观看免费| 成熟少妇高潮喷水视频| 亚洲欧洲精品一区二区精品久久久| 午夜免费激情av| 岛国在线免费视频观看| 国产精华一区二区三区| 色播亚洲综合网| 国产精品,欧美在线| 日韩欧美精品v在线| 久久久国产欧美日韩av| 一区二区三区高清视频在线| 久久香蕉精品热| 日韩欧美 国产精品| 小说图片视频综合网站| 亚洲一区二区三区不卡视频| 中文字幕久久专区| 国产精品一区二区免费欧美| 久久热在线av| 淫妇啪啪啪对白视频| 成人三级做爰电影| 老司机福利观看| 成年女人毛片免费观看观看9| 在线免费观看不下载黄p国产 | 性色av乱码一区二区三区2| 欧美性猛交╳xxx乱大交人| 久久香蕉国产精品| 十八禁网站免费在线| 精品乱码久久久久久99久播| 宅男免费午夜| 欧美日本亚洲视频在线播放| 国产av在哪里看| 国产精品影院久久| 亚洲av成人不卡在线观看播放网| 身体一侧抽搐| 最好的美女福利视频网| 人人妻人人澡欧美一区二区| 2021天堂中文幕一二区在线观| avwww免费| 国产一区二区在线观看日韩 | 99久久精品一区二区三区| 国产av麻豆久久久久久久| 精品一区二区三区av网在线观看| 91麻豆av在线| 国产精品亚洲av一区麻豆| 国产精品亚洲av一区麻豆| 欧美极品一区二区三区四区| 午夜福利成人在线免费观看| 久久伊人香网站| 成人三级做爰电影| 国产高清videossex| 午夜精品一区二区三区免费看| 脱女人内裤的视频| 国产欧美日韩精品一区二区| 亚洲国产欧美网| 亚洲中文日韩欧美视频| 国产激情偷乱视频一区二区| 黑人欧美特级aaaaaa片| 99热这里只有是精品50| 久久久久久久久免费视频了| 久久人妻av系列| 精品久久久久久久久久免费视频| 国产成人av激情在线播放| 久久久久久久久中文| 国产精品一区二区三区四区免费观看 | 亚洲中文日韩欧美视频| 99久久精品一区二区三区| 丁香六月欧美| 国产亚洲精品一区二区www| 成年版毛片免费区| 不卡av一区二区三区| 搞女人的毛片| 色老头精品视频在线观看| 亚洲成av人片在线播放无| 欧美乱码精品一区二区三区| 在线十欧美十亚洲十日本专区| 9191精品国产免费久久| 国产亚洲精品综合一区在线观看| 成人永久免费在线观看视频| 网址你懂的国产日韩在线| 亚洲欧美精品综合一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲专区中文字幕在线| 国产三级中文精品| 亚洲国产精品sss在线观看| 精品乱码久久久久久99久播| 天堂动漫精品| 国产午夜精品久久久久久| 51午夜福利影视在线观看| 两性夫妻黄色片| 国内精品久久久久久久电影| 精品日产1卡2卡| 久久午夜综合久久蜜桃| 可以在线观看毛片的网站| 久久精品亚洲精品国产色婷小说| 国产精品乱码一区二三区的特点| 一二三四在线观看免费中文在| 老熟妇仑乱视频hdxx| 波多野结衣高清无吗| 久久久久亚洲av毛片大全| 日韩精品青青久久久久久| 男女下面进入的视频免费午夜| 亚洲色图 男人天堂 中文字幕| 国内精品久久久久久久电影| 国产精品久久电影中文字幕| 久久国产精品人妻蜜桃| 日日摸夜夜添夜夜添小说| 搞女人的毛片| 亚洲成人免费电影在线观看| ponron亚洲| 亚洲激情在线av| 国产精品日韩av在线免费观看| 国产成人精品久久二区二区91| 久久九九热精品免费| 99久久精品国产亚洲精品| 国内久久婷婷六月综合欲色啪| 久久人妻av系列| 免费一级毛片在线播放高清视频| 91九色精品人成在线观看| 搡老熟女国产l中国老女人| 免费无遮挡裸体视频| 动漫黄色视频在线观看| 日韩人妻高清精品专区| av国产免费在线观看| 色吧在线观看| 长腿黑丝高跟| 最近最新中文字幕大全免费视频| 亚洲色图 男人天堂 中文字幕| 夜夜夜夜夜久久久久| 欧美中文综合在线视频| 一级毛片高清免费大全| 亚洲午夜精品一区,二区,三区| 免费av不卡在线播放| av在线蜜桃| 97人妻精品一区二区三区麻豆| 亚洲乱码一区二区免费版| 草草在线视频免费看| 国产精品九九99| 岛国在线免费视频观看| 国产欧美日韩一区二区精品| 国产单亲对白刺激| 欧美另类亚洲清纯唯美| 国产1区2区3区精品| 日韩精品中文字幕看吧| 亚洲中文字幕一区二区三区有码在线看 | 国产三级在线视频| 久久中文看片网| 亚洲九九香蕉| www.www免费av| 日本撒尿小便嘘嘘汇集6| 蜜桃久久精品国产亚洲av| 国产av一区在线观看免费| 小蜜桃在线观看免费完整版高清| 少妇熟女aⅴ在线视频| 日本成人三级电影网站| 最好的美女福利视频网| 波多野结衣高清无吗| 午夜免费观看网址| 国产精品日韩av在线免费观看| 男女视频在线观看网站免费| 一进一出抽搐动态| 天堂网av新在线| 日韩精品中文字幕看吧| 又粗又爽又猛毛片免费看| 免费电影在线观看免费观看| 日韩欧美免费精品| 亚洲成人久久爱视频| 亚洲av电影不卡..在线观看| 国产探花在线观看一区二区| 在线观看66精品国产| 丁香欧美五月| 国产成人福利小说| a在线观看视频网站| av片东京热男人的天堂| 亚洲avbb在线观看| 丁香六月欧美| 国产精品1区2区在线观看.| 人人妻人人澡欧美一区二区| 国产精品一区二区三区四区久久| 18美女黄网站色大片免费观看| 亚洲人成伊人成综合网2020| 国产人伦9x9x在线观看| 久久国产精品人妻蜜桃| 国产伦精品一区二区三区四那| 精品国内亚洲2022精品成人| 免费一级毛片在线播放高清视频| 国产极品精品免费视频能看的| 在线观看一区二区三区| 听说在线观看完整版免费高清| 日本 欧美在线| 欧美3d第一页| 99热精品在线国产| 婷婷亚洲欧美| 18禁黄网站禁片免费观看直播| 久久久水蜜桃国产精品网| 精品国产乱码久久久久久男人| 两性夫妻黄色片| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩国产亚洲二区| 熟女电影av网| 欧美在线黄色| 久久久久久大精品| 久久久国产成人免费| 日本一本二区三区精品| 免费看光身美女| 日韩欧美国产一区二区入口| 日本精品一区二区三区蜜桃| 高清毛片免费观看视频网站| 99久久精品国产亚洲精品| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区色噜噜| 色综合亚洲欧美另类图片| 91久久精品国产一区二区成人 | 久久久久久大精品| 俺也久久电影网| 日本在线视频免费播放| 小蜜桃在线观看免费完整版高清| 亚洲一区二区三区色噜噜| 久久99热这里只有精品18| 国产成人精品久久二区二区91| 国产97色在线日韩免费| 99国产精品99久久久久| 国产高清视频在线播放一区| 18禁国产床啪视频网站| 日韩av在线大香蕉| 国产激情久久老熟女| 亚洲自偷自拍图片 自拍| 伊人久久大香线蕉亚洲五| 1024香蕉在线观看| 91在线精品国自产拍蜜月 | 国产真实乱freesex| 久久久久精品国产欧美久久久| www日本在线高清视频| 巨乳人妻的诱惑在线观看| 国产伦人伦偷精品视频| 国产日本99.免费观看| 亚洲无线观看免费| 又粗又爽又猛毛片免费看| 老汉色av国产亚洲站长工具| 欧美成狂野欧美在线观看| 天堂影院成人在线观看| 婷婷亚洲欧美| 午夜免费成人在线视频| 国产av不卡久久| 好男人在线观看高清免费视频| 久久久久精品国产欧美久久久| 欧美日韩国产亚洲二区| 18禁观看日本| 床上黄色一级片| 久久久国产成人免费| 亚洲av成人精品一区久久| 亚洲中文字幕日韩| 欧美午夜高清在线| 久久亚洲真实| 亚洲国产中文字幕在线视频| 国产精品国产高清国产av| avwww免费| 久久亚洲真实| 欧美大码av| 亚洲成人久久性| 久久久水蜜桃国产精品网| av女优亚洲男人天堂 | 草草在线视频免费看| 国产精品久久久久久人妻精品电影| 人妻久久中文字幕网| h日本视频在线播放| 久久99热这里只有精品18| 性色avwww在线观看| 国内精品久久久久精免费| 两个人的视频大全免费| 一二三四社区在线视频社区8| 久久久国产成人免费| 看片在线看免费视频| ponron亚洲| 亚洲av电影在线进入| 国产精品亚洲av一区麻豆| 老司机深夜福利视频在线观看| 一区福利在线观看| cao死你这个sao货| 亚洲欧洲精品一区二区精品久久久| 精品不卡国产一区二区三区| 丰满人妻一区二区三区视频av | 日本精品一区二区三区蜜桃| 可以在线观看的亚洲视频| 国产成人精品久久二区二区免费| 婷婷精品国产亚洲av在线| av在线天堂中文字幕| 色综合欧美亚洲国产小说| 国语自产精品视频在线第100页| 亚洲激情在线av| 欧美日韩福利视频一区二区| 特级一级黄色大片| 热99在线观看视频| 久久性视频一级片| 午夜精品一区二区三区免费看| 97碰自拍视频| 两个人视频免费观看高清| 亚洲第一电影网av| 国产极品精品免费视频能看的| 久久久国产精品麻豆| 狂野欧美激情性xxxx| h日本视频在线播放| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 久久精品亚洲精品国产色婷小说| 色在线成人网| 中亚洲国语对白在线视频| 免费在线观看日本一区| 亚洲乱码一区二区免费版| 久久天躁狠狠躁夜夜2o2o| 国产av麻豆久久久久久久| 精华霜和精华液先用哪个| 午夜免费激情av| 成人精品一区二区免费| 国产一区在线观看成人免费| 成年人黄色毛片网站| 91九色精品人成在线观看| 国产不卡一卡二| 久久午夜亚洲精品久久| 国产免费男女视频| 亚洲,欧美精品.| 国产成人av教育| 国内精品一区二区在线观看| 岛国视频午夜一区免费看| 久久久久性生活片| 免费无遮挡裸体视频| 亚洲国产欧洲综合997久久,| 国产视频内射| 亚洲精品一卡2卡三卡4卡5卡| 十八禁网站免费在线| x7x7x7水蜜桃| 黄色成人免费大全| 丰满人妻一区二区三区视频av | 日本精品一区二区三区蜜桃| x7x7x7水蜜桃| 激情在线观看视频在线高清| 亚洲国产欧洲综合997久久,| 精品无人区乱码1区二区| 岛国在线观看网站| 日本精品一区二区三区蜜桃| 女人被狂操c到高潮| 日本成人三级电影网站| 久久久久久久午夜电影| 午夜福利成人在线免费观看| 后天国语完整版免费观看| 亚洲 欧美一区二区三区| 久久人人精品亚洲av| 欧美xxxx黑人xx丫x性爽| 一级毛片高清免费大全| 热99在线观看视频| 成人av在线播放网站| 欧美日韩瑟瑟在线播放| 五月玫瑰六月丁香| 亚洲成a人片在线一区二区| av福利片在线观看| 99在线人妻在线中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜福利久久久久久| av在线天堂中文字幕| 黄色丝袜av网址大全| 午夜福利在线在线| 亚洲五月天丁香| 免费一级毛片在线播放高清视频| 国产精品久久久久久亚洲av鲁大| 久久精品91无色码中文字幕| 在线国产一区二区在线| 热99在线观看视频| 国产美女午夜福利| 成人国产综合亚洲| 亚洲av日韩精品久久久久久密| 亚洲欧洲精品一区二区精品久久久| 神马国产精品三级电影在线观看| 99riav亚洲国产免费| 在线观看舔阴道视频| 久久久久久久久久黄片| a级毛片a级免费在线| 99久久精品国产亚洲精品| 一个人看的www免费观看视频| 视频区欧美日本亚洲| 亚洲欧洲精品一区二区精品久久久| 国产男靠女视频免费网站| 黄频高清免费视频| 中文字幕久久专区| 久久精品夜夜夜夜夜久久蜜豆| 一进一出抽搐动态| 在线观看66精品国产| 99国产综合亚洲精品| 日韩大尺度精品在线看网址| 99re在线观看精品视频| 无人区码免费观看不卡| 免费无遮挡裸体视频| 亚洲熟妇熟女久久| 亚洲精品乱码久久久v下载方式 | 亚洲一区二区三区色噜噜| 91在线精品国自产拍蜜月 | 怎么达到女性高潮| 99re在线观看精品视频| 欧美日韩福利视频一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 久久中文看片网| 女人被狂操c到高潮| 好男人电影高清在线观看| 99久久无色码亚洲精品果冻| 国产精品一区二区免费欧美| 国产真人三级小视频在线观看| 99久久国产精品久久久| 亚洲中文字幕日韩| 99国产综合亚洲精品| 日韩欧美国产在线观看| 久久国产乱子伦精品免费另类| a在线观看视频网站| 老司机深夜福利视频在线观看| 亚洲av美国av| 99久久综合精品五月天人人| 嫩草影视91久久| 这个男人来自地球电影免费观看| 老司机深夜福利视频在线观看| 亚洲国产欧美人成| 久久午夜综合久久蜜桃| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 精品不卡国产一区二区三区| 精品免费久久久久久久清纯| 久久久久久久久中文| 国产高清视频在线观看网站| 法律面前人人平等表现在哪些方面| 午夜免费成人在线视频| tocl精华| 亚洲人成网站在线播放欧美日韩| 在线观看美女被高潮喷水网站 | 久久香蕉精品热| 狠狠狠狠99中文字幕| 亚洲av成人一区二区三| 久久久久久久久久黄片| 熟女少妇亚洲综合色aaa.| 午夜福利在线在线| 亚洲色图 男人天堂 中文字幕| h日本视频在线播放| 日韩欧美免费精品| 看黄色毛片网站| АⅤ资源中文在线天堂| 国产高清三级在线| 久久精品国产亚洲av香蕉五月| 性欧美人与动物交配| 日韩欧美国产一区二区入口| 男人和女人高潮做爰伦理| 欧美日韩综合久久久久久 | 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 国产精品影院久久| 午夜精品久久久久久毛片777| 欧美日韩综合久久久久久 | 精品国产美女av久久久久小说| www国产在线视频色| 亚洲精品一区av在线观看| 午夜精品久久久久久毛片777| 亚洲 欧美 日韩 在线 免费| 亚洲国产欧美一区二区综合| 国产精品国产高清国产av| 国内毛片毛片毛片毛片毛片| 国产精品亚洲美女久久久| 99国产综合亚洲精品| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| www.熟女人妻精品国产| 草草在线视频免费看| 国产精品久久电影中文字幕| 一区福利在线观看| 成人亚洲精品av一区二区| 欧美av亚洲av综合av国产av| 色视频www国产| 欧美+亚洲+日韩+国产| tocl精华| 三级国产精品欧美在线观看 | 亚洲专区国产一区二区| 久久亚洲真实| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久精品吃奶| av天堂中文字幕网| 两性夫妻黄色片| 身体一侧抽搐| 午夜福利视频1000在线观看| 麻豆成人av在线观看| 国产精品 国内视频| 亚洲国产精品999在线| 欧美色欧美亚洲另类二区| 日韩欧美国产一区二区入口| 黄色片一级片一级黄色片| 男人舔奶头视频| 最新在线观看一区二区三区| 日本与韩国留学比较| 精品无人区乱码1区二区| 日韩欧美 国产精品| 无人区码免费观看不卡| 人人妻人人澡欧美一区二区| 免费在线观看成人毛片| 久久久久亚洲av毛片大全| 很黄的视频免费| 在线视频色国产色| 免费看日本二区| 国产一区在线观看成人免费| 色综合欧美亚洲国产小说| 一区福利在线观看| 国内久久婷婷六月综合欲色啪| 在线永久观看黄色视频| 久久午夜综合久久蜜桃| 国产91精品成人一区二区三区| 一本精品99久久精品77| 亚洲av免费在线观看| 男女那种视频在线观看| 国产欧美日韩精品亚洲av| 成人一区二区视频在线观看| 成年女人永久免费观看视频| 欧美乱码精品一区二区三区| 一个人免费在线观看电影 | 男人舔女人下体高潮全视频| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 一本精品99久久精品77| 五月玫瑰六月丁香| 99精品在免费线老司机午夜| 国产高清视频在线观看网站| 可以在线观看毛片的网站| 亚洲无线在线观看| 99久久99久久久精品蜜桃| 免费观看人在逋| 亚洲最大成人中文| a级毛片在线看网站| 亚洲国产高清在线一区二区三| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三| 麻豆成人午夜福利视频| 日韩人妻高清精品专区| 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av在线| 又爽又黄无遮挡网站| 色噜噜av男人的天堂激情| 成人18禁在线播放| 99热只有精品国产| 真人做人爱边吃奶动态| 极品教师在线免费播放| 老熟妇乱子伦视频在线观看| 免费看十八禁软件| 国产高潮美女av|