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

    基于響應(yīng)面和遺傳算法的尾座式無人機(jī)結(jié)構(gòu)參數(shù)優(yōu)化

    2019-06-04 01:10:46劉文帥姚小敏李超群張夢(mèng)飛淡煦珈韓文霆
    關(guān)鍵詞:后掠角尾座小翼

    劉文帥 姚小敏 李超群 張夢(mèng)飛 淡煦珈 韓文霆

    (1.西北農(nóng)林科技大學(xué)機(jī)械與電子工程學(xué)院, 陜西楊凌 712100; 2.西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室, 陜西楊凌 712100; 3.南京禾譜航空科技有限公司, 南京 211300; 4.西北農(nóng)林科技大學(xué)水土保持研究所, 陜西楊凌 712100)

    0 引言

    尾座式無人機(jī)是一種尾座坐地式垂直起降飛行器,融合了多旋翼的垂直起降特性和固定翼的高效巡航特性,有效解決了多旋翼航程短和飛行效率低,固定翼無人機(jī)起降難、部署時(shí)間長(zhǎng)和輔助保障設(shè)備多等問題,可以滿足日益復(fù)雜的農(nóng)業(yè)低空遙感作業(yè)要求[1-3]。相較于其他垂直起降無人機(jī),尾座式無人機(jī)無需額外的轉(zhuǎn)動(dòng)控制部件,具有結(jié)構(gòu)簡(jiǎn)單、質(zhì)量輕、整機(jī)自然穩(wěn)定性好和操作簡(jiǎn)便等優(yōu)點(diǎn),近些年得到廣泛關(guān)注。

    國(guó)內(nèi)外學(xué)者常用計(jì)算流體動(dòng)力學(xué)(Computational fluid dynamics,CFD)方法進(jìn)行尾座式無人機(jī)的設(shè)計(jì)和氣動(dòng)優(yōu)化,通過CFD數(shù)值模擬方法能夠獲得無人機(jī)的升力、阻力和升阻比等氣動(dòng)特性參數(shù),通過氣動(dòng)參數(shù)的分析,在機(jī)翼、螺旋槳、小翼和副翼等部件的結(jié)構(gòu)參數(shù)優(yōu)化、氣動(dòng)布局優(yōu)化和動(dòng)力學(xué)建模等方面已取得了一定的研究成果[4-17]。由于CFD計(jì)算方法只能獲得離散的數(shù)據(jù)點(diǎn),很難進(jìn)行全局范圍內(nèi)最優(yōu)結(jié)構(gòu)參數(shù)的求解,因此基于全局最優(yōu)點(diǎn)求解的響應(yīng)面和多目標(biāo)遺傳算法(Multi-objective genetic algorithm,MOGA)受到越來越多的關(guān)注,并逐漸應(yīng)用到結(jié)構(gòu)設(shè)計(jì)中。在復(fù)合材料無人機(jī)的結(jié)構(gòu)強(qiáng)度改進(jìn)、機(jī)翼的減重、動(dòng)力系統(tǒng)飛行效率的提高、換熱器結(jié)構(gòu)參數(shù)優(yōu)化和加工機(jī)床的結(jié)構(gòu)設(shè)計(jì)等方面取得了很多研究成果[18-25]。

    本文將響應(yīng)面和多目標(biāo)遺傳算法應(yīng)用于尾座式無人機(jī)的優(yōu)化設(shè)計(jì),構(gòu)建無人機(jī)的翼展長(zhǎng)、后掠角、小翼高和小翼厚4個(gè)結(jié)構(gòu)參數(shù)與升阻比和阻力的響應(yīng)面模型,以升阻比最大、阻力最小為優(yōu)化目標(biāo),采用多目標(biāo)遺傳算法求解最優(yōu)結(jié)構(gòu)參數(shù)。

    1 材料與方法

    1.1 尾座式無人機(jī)結(jié)構(gòu)與參數(shù)

    1.1.1結(jié)構(gòu)布局

    設(shè)計(jì)的尾座式無人機(jī)采用飛翼雙動(dòng)力前拉布局形式,如圖1所示,主要由左右對(duì)稱的機(jī)翼、小翼、電機(jī)座、正反轉(zhuǎn)電機(jī)和螺旋槳組成。碳桿貫穿機(jī)翼、電機(jī)座和小翼,起到機(jī)身加固和安裝定位的作用。機(jī)翼通過翼根膠粘的方式組成機(jī)身,小翼和電機(jī)座通過膠粘的方式固定在機(jī)身上,并且保持小翼和電機(jī)座軸線與機(jī)身軸線平行。電機(jī)和螺旋槳通過螺紋固定在電機(jī)座上。

    圖1 尾座式無人機(jī)結(jié)構(gòu)示意圖Fig.1 Schematic of tail-sitter UAV1.機(jī)翼 2.小翼 3.舵面 4.電機(jī)座 5.電機(jī) 6.螺旋槳

    尾座式無人機(jī)有多旋翼起降和固定翼巡航兩種飛行模式。多旋翼起降時(shí),通過正反轉(zhuǎn)螺旋槳為無人機(jī)提供垂直上升的動(dòng)力;并通過兩個(gè)反轉(zhuǎn)差速電機(jī)的轉(zhuǎn)速進(jìn)行左右的姿態(tài)控制;通過舵面的俯仰運(yùn)動(dòng)進(jìn)行無人機(jī)的俯仰姿態(tài)調(diào)節(jié)。當(dāng)飛機(jī)達(dá)到一定高度時(shí),切換為固定翼巡航模式,電機(jī)工作模式轉(zhuǎn)換為反轉(zhuǎn)同速電機(jī),為固定翼巡航狀態(tài)提供前拉動(dòng)力;舵面的工作模式轉(zhuǎn)換為固定翼混控模式,進(jìn)行固定翼的俯仰和盤旋姿態(tài)的控制。

    1.1.2結(jié)構(gòu)參數(shù)范圍

    尾座式無人機(jī)采用飛翼布局形式,根據(jù)文獻(xiàn)[4-5]可知,翼展長(zhǎng)、后掠角、小翼高和小翼厚是影響尾座式無人機(jī)升阻比和阻力的主要因素,如圖2所示。

    圖2 尾座式無人機(jī)的結(jié)構(gòu)參數(shù)Fig.2 Structural parameters of tail-sitter UAV1.機(jī)翼 2.小翼

    由文獻(xiàn)知現(xiàn)有的小型固定翼無人機(jī)的翼展長(zhǎng)b的范圍為900~1 200 mm,后掠角Λw為15°~60°,小翼高lv為0.05b~0.3b,小翼厚h為3~7 mm,確定尾座式無人機(jī)結(jié)構(gòu)參數(shù)范圍為:翼展長(zhǎng)900~1 300 mm,后掠角為15°~35°,小翼高25~45 mm,小翼厚3~7 mm。本文以翼展長(zhǎng)1 000 mm,后掠角20°、小翼高30 mm和小翼厚6 mm為原始樣機(jī)進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化。

    1.1.3結(jié)構(gòu)參數(shù)優(yōu)化流程

    尾座式無人機(jī)結(jié)構(gòu)參數(shù)優(yōu)化流程如圖3所示。

    圖3 結(jié)構(gòu)參數(shù)優(yōu)化流程圖Fig.3 Flow chart of structural parameter optimization

    本文對(duì)尾座式無人機(jī)的翼展長(zhǎng)、后掠角、小翼高和小翼厚4個(gè)結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化,通過CCD進(jìn)行4個(gè)結(jié)構(gòu)參數(shù)的組合設(shè)計(jì)得到25個(gè)組合樣本點(diǎn)和3個(gè)對(duì)照樣本點(diǎn),在CFD中模擬計(jì)算得到各樣本點(diǎn)的升阻比和阻力,并通過風(fēng)洞試驗(yàn)驗(yàn)證CFD數(shù)值模擬方法的可靠性。利用多項(xiàng)式回歸的方法分別建立了升阻比和阻力與結(jié)構(gòu)參數(shù)的響應(yīng)面模型,并由對(duì)照樣本點(diǎn)的CFD模擬值進(jìn)行驗(yàn)證,以升阻比最大和阻力最小為優(yōu)化目標(biāo),利用多目標(biāo)遺傳算法進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化得到最優(yōu)結(jié)構(gòu)參數(shù)。

    1.2 無人機(jī)外流場(chǎng)數(shù)值模擬方法

    尾座式無人機(jī)的結(jié)構(gòu)優(yōu)化是反復(fù)權(quán)衡續(xù)航時(shí)間和姿態(tài)轉(zhuǎn)換穩(wěn)定性的過程,其中較大的升阻比能夠提高無人機(jī)的氣動(dòng)效率,延長(zhǎng)續(xù)航時(shí)間;較小的阻力能夠降低無人機(jī)的側(cè)向滑移,提高轉(zhuǎn)換過程的穩(wěn)定性。本文采用CFD數(shù)值模擬的方法對(duì)尾座式無人機(jī)的升阻比和阻力進(jìn)行模擬計(jì)算。由于螺旋和電機(jī)布置在機(jī)翼的氣動(dòng)弦長(zhǎng)上,螺旋槳在機(jī)翼上作用力可視為0°攻角下的氣動(dòng)力,相對(duì)較小可以忽略。并且由于無人機(jī)為左右對(duì)稱布局,利用CATIA軟件對(duì)機(jī)身左側(cè)建立三維模型,并去掉螺旋槳、電機(jī)和電機(jī)座等結(jié)構(gòu)。將三維模型導(dǎo)入ANSYS的Geometry模塊構(gòu)建外流場(chǎng)的三維實(shí)體模型,利用ICEM對(duì)外流場(chǎng)進(jìn)行非結(jié)構(gòu)體網(wǎng)格劃分,而且為了更好地捕捉無人機(jī)近壁面外流場(chǎng)區(qū)域的流動(dòng)細(xì)節(jié),對(duì)無人機(jī)周圍網(wǎng)格進(jìn)行局部加密處理,如圖4a所示,最終生成的外流場(chǎng)網(wǎng)格數(shù)量為230萬個(gè)節(jié)點(diǎn),無人機(jī)近壁面處的網(wǎng)格如圖4b所示。

    圖4 無人機(jī)外流場(chǎng)網(wǎng)格Fig.4 Grid of UAV outflow

    尾座式無人機(jī)的姿態(tài)轉(zhuǎn)換過程是攻角逐漸增大的過程,當(dāng)攻角超過18°時(shí)進(jìn)入失速狀態(tài),為非定常流動(dòng)模式,現(xiàn)有的CFD軟件無法進(jìn)行精準(zhǔn)模擬,本文模擬不同結(jié)構(gòu)的無人機(jī)在攻角為4°~12°時(shí)的升阻比和阻力。無人機(jī)的飛行速度為12 m/s、雷諾數(shù)為50 000時(shí)為湍流,選用SSTk-ω模型,入口速度為12 m/s,速度夾角為4°~12°,出口壓力為標(biāo)準(zhǔn)大氣壓,計(jì)算區(qū)域的邊界假定為無滑移條件,用標(biāo)準(zhǔn)的邊界方程進(jìn)行近壁面處理,用SIMPLE算法進(jìn)行方程求解。

    1.3 組合樣本點(diǎn)設(shè)計(jì)及模型構(gòu)建方法

    1.3.1組合樣本點(diǎn)設(shè)計(jì)

    CCD設(shè)計(jì)是在兩因子全水平試驗(yàn)的基礎(chǔ)上增加了重復(fù)使用的中心點(diǎn)和軸向點(diǎn),使得每個(gè)因子有5個(gè)變量水平(如圖5所示),以適應(yīng)二次模型的擬合。利用CCD對(duì)翼展長(zhǎng)、后掠角、小翼高和小翼厚4個(gè)結(jié)構(gòu)參數(shù)進(jìn)行設(shè)計(jì)得到16個(gè)立方點(diǎn)、8個(gè)軸向點(diǎn)和1個(gè)中心點(diǎn),總計(jì)25個(gè)組合樣本點(diǎn)。在組合樣本點(diǎn)之外隨機(jī)抽取3個(gè)對(duì)照樣本點(diǎn),用于驗(yàn)證響應(yīng)面模型的通用性。

    圖5 CCD樣本點(diǎn)設(shè)計(jì)Fig.5 Design of CCD sample points

    1.3.2模型構(gòu)建方法

    利用CFD計(jì)算25個(gè)組合樣本點(diǎn)的升阻比和阻力,并采用二次多項(xiàng)式擬合的方法分別構(gòu)建升阻比和阻力與4個(gè)結(jié)構(gòu)參數(shù)的響應(yīng)面模型。響應(yīng)面模型表示為

    yi=Cti+ε

    (1)

    其中

    C=[C1C2C3C4]

    式中yi——第i個(gè)點(diǎn)的響應(yīng)面模型估算值

    ti——第i個(gè)點(diǎn)的回歸系數(shù)

    C——回歸參數(shù)矩陣

    ε——相對(duì)誤差

    1.4 多目標(biāo)遺傳算法

    尾座式無人機(jī)以固定翼巡航模式進(jìn)行目標(biāo)區(qū)域的遙感圖像采集,在設(shè)計(jì)過程中應(yīng)保證無人機(jī)獲得較大的升阻比。姿態(tài)轉(zhuǎn)換過程中俯仰角從0°逐漸增加到90°,在該過程中為保證姿態(tài)轉(zhuǎn)換過程的平穩(wěn)性,無人機(jī)的阻力應(yīng)盡量小,則無人機(jī)的多目標(biāo)遺傳算法公式表示為

    maxK(x)=[K1(x)K2(x) …KN(x)]

    (2)

    minD(x)=[D1(x)D2(x) …DN(x)]

    (3)

    其中

    x=(x1,x2,x3,x4)

    式中x——結(jié)構(gòu)參數(shù)向量

    xd——結(jié)構(gòu)參數(shù)

    MOGA的求解過程是一個(gè)種群構(gòu)建、個(gè)體求解和個(gè)體尋優(yōu)的循環(huán)過程,直至達(dá)到設(shè)定的目標(biāo)函數(shù)。利用1.3.2節(jié)構(gòu)建的響應(yīng)面模型進(jìn)行個(gè)體求解,并通過約束條件進(jìn)行個(gè)體篩選,并重新構(gòu)建種群,循環(huán)迭代直至尋找到升阻比最大且阻力最小的最優(yōu)個(gè)體,即為最優(yōu)結(jié)構(gòu)參數(shù)。

    1.5 數(shù)值模擬方法驗(yàn)證

    1.5.1試驗(yàn)設(shè)備

    基于升阻比和阻力的CFD模擬值建立響應(yīng)面模型,采用風(fēng)洞試驗(yàn)進(jìn)行CFD數(shù)值模擬方法的驗(yàn)證。在25個(gè)CCD組合樣本中隨機(jī)選擇一個(gè)樣本點(diǎn)參數(shù)為:翼展長(zhǎng)1 000 mm、后掠角27°、小翼高40 mm和小翼厚5 mm,并分別測(cè)定樣本在不同風(fēng)速下的升阻比和阻力。由于風(fēng)洞的風(fēng)場(chǎng)直徑(800 mm)限制,根據(jù)風(fēng)洞試驗(yàn)的相似準(zhǔn)則,將樣本按照1∶0.6的比例進(jìn)行縮放,得到風(fēng)洞模型尺寸為翼展長(zhǎng)600 mm、后掠角27°、小翼高24 mm和小翼厚3 mm,采用3D打印技術(shù)進(jìn)行風(fēng)洞模型的加工,如圖6所示。

    圖6 風(fēng)洞試驗(yàn)樣機(jī)Fig.6 Wind tunnel experimental prototype

    1.5.2試驗(yàn)條件與方案

    風(fēng)洞試驗(yàn)系統(tǒng)包括無人機(jī)模型、飛機(jī)支架、六分量應(yīng)變天平、風(fēng)機(jī)、變頻器和數(shù)據(jù)采集控制系統(tǒng)等,如圖7所示。無人機(jī)模型加工過程中,為保證機(jī)翼的流線分布形式和表面光潔,采用展向掃描遞推的打印方式進(jìn)行模型的加工,并對(duì)凹槽和縫隙進(jìn)行蒙皮處理。為降低支架對(duì)無人機(jī)尾部氣流的干擾,采用天平內(nèi)嵌機(jī)身和支桿延伸的安裝方式進(jìn)行樣機(jī)的固定。為保證模型在風(fēng)場(chǎng)中數(shù)據(jù)的準(zhǔn)確性,試驗(yàn)過程中無人機(jī)保持在風(fēng)場(chǎng)的中心位置(離地1.4 m)。試驗(yàn)過程中變頻器控制風(fēng)機(jī)轉(zhuǎn)速待風(fēng)速平穩(wěn)后,由六分量應(yīng)變天平進(jìn)行升力和阻力的采集,并由數(shù)據(jù)分析系統(tǒng)對(duì)采集到的1 000組數(shù)據(jù)進(jìn)行處理,最終得到無人機(jī)的升阻比和阻力。

    圖7 風(fēng)洞試驗(yàn)系統(tǒng)Fig.7 Wind tunnel experiment system1.天平支架 2.風(fēng)機(jī) 3.試驗(yàn)樣機(jī) 4.六分量應(yīng)變天平

    風(fēng)洞試驗(yàn)在西北工業(yè)大學(xué)翼型葉柵空氣動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室內(nèi)進(jìn)行,試驗(yàn)中測(cè)定無人機(jī)俯仰角為8°,通過變頻器將風(fēng)速穩(wěn)定在12~20 m/s之間,間隔2 m/s,測(cè)定樣機(jī)在不同風(fēng)速下的升阻比和阻力,并將CFD的模擬值與風(fēng)洞試驗(yàn)測(cè)量值進(jìn)行對(duì)比。

    1.5.3數(shù)值模擬結(jié)果驗(yàn)證

    無人機(jī)工況的改變,實(shí)際上是雷諾數(shù)的變化。依據(jù)雷諾數(shù)相似理論,幾何相似的前提下,相同的雷諾數(shù)具有相似的氣動(dòng)特性[26]。繪制樣本和風(fēng)洞模型在不同雷諾數(shù)下的升阻比和阻力散點(diǎn)圖,并與風(fēng)洞數(shù)據(jù)進(jìn)行對(duì)比,如圖8所示。

    圖8 風(fēng)洞試驗(yàn)驗(yàn)證結(jié)果Fig.8 Wind tunnel test results

    風(fēng)洞模型與樣本的升阻比相對(duì)誤差小于2.3%,阻力相對(duì)誤差小于2.7%,表明基于等比縮放的方法進(jìn)行數(shù)值模擬具有較高的可信度;風(fēng)洞模型的仿真數(shù)據(jù)與試驗(yàn)的升阻比相對(duì)誤差小于7.3%,阻力相對(duì)誤差小于5.8%;樣本的仿真數(shù)據(jù)與試驗(yàn)的升阻比相對(duì)誤差小于8.0%,阻力的平均相對(duì)誤差為5.7%,在允許的誤差范圍內(nèi)。風(fēng)洞模型與樣本的誤差原因是縮放帶來的原理誤差;風(fēng)洞模型仿真數(shù)據(jù)與試驗(yàn)數(shù)據(jù)的誤差原因是風(fēng)洞模型采用三維打印加工,模型的表面粗糙度會(huì)影響機(jī)翼近壁面的流場(chǎng)分布,風(fēng)洞模型阻力較小,試驗(yàn)過程中模型的抖動(dòng)和信號(hào)干擾等因素均會(huì)影響阻力的測(cè)量精度;樣本仿真數(shù)據(jù)與試驗(yàn)的誤差原因是由縮放原理和試驗(yàn)環(huán)境共同導(dǎo)致,但試驗(yàn)結(jié)果誤差小于8%,在允許范圍內(nèi),表明CFD數(shù)值模擬方法可靠,可以用于尾座式無人機(jī)的外流場(chǎng)模擬。

    2 結(jié)果與分析

    2.1 結(jié)構(gòu)參數(shù)對(duì)升阻比和阻力的影響

    2.1.1結(jié)構(gòu)參數(shù)對(duì)升阻比的影響

    無人機(jī)的升阻比主要與翼展長(zhǎng)b、后掠角Λw、小翼厚h和小翼高lv等結(jié)構(gòu)參數(shù)有關(guān),建立各結(jié)構(gòu)參數(shù)在不同攻角α下的升阻比變化曲線,如圖9~12所示。

    由圖9可知,翼展長(zhǎng)從900 mm增加到1 300 mm的過程中無人機(jī)的升阻比呈線性增加。因?yàn)闊o人機(jī)的升阻比與展弦比有關(guān),展弦比越大升阻比越大。當(dāng)其它結(jié)構(gòu)參數(shù)確定時(shí),展弦比與翼展長(zhǎng)呈線性關(guān)系,即翼展長(zhǎng)與升阻比線性相關(guān)。

    圖10 后掠角與升阻比的關(guān)系曲線Fig.10 Relationship curves between sweep angle and lift drag ratio

    圖11 小翼高與升阻比的關(guān)系曲線Fig.11 Relationship curves between wing height and lift drag ratio

    在攻角不同的情況下,升阻比隨后掠角和小翼高的變化如圖10、11所示。在相同的攻角下升阻比隨后掠角的增大而略有減小,隨小翼高的增加而略有增大。當(dāng)攻角為12°時(shí),后掠角為35°的升阻比比后掠角為15°的降低了1%;小翼高為25 mm的升阻比比小翼高為45 mm的升阻比下降了1%,表明后掠角和小翼高對(duì)升阻比的影響較小。因?yàn)樯璞戎饕c展弦比相關(guān),展弦比主要與翼展長(zhǎng)和氣動(dòng)弦長(zhǎng)相關(guān),在梯形翼中隨著后掠角和小翼高的增加氣動(dòng)弦長(zhǎng)的變化較小,因此后掠角對(duì)升阻比K影響較小。

    圖12 小翼厚與升阻比的關(guān)系曲線Fig.12 Relationship curves between wing thickness and lift drag ratio

    由圖12可知,在小攻角4°~8°時(shí),升阻比隨小翼厚的增加而逐漸減小,并且4°攻角時(shí)升阻比的下降速率最大;在大攻角10°~12°時(shí),隨著小翼厚的增加,升阻比先減小后增大。因?yàn)樵谛」ソ菚r(shí),小翼主要用于降低翼尖擾流,減小誘導(dǎo)阻力提高升阻比,并且隨著小翼厚的增加,小翼表面的擾流情況加劇,反而會(huì)增加誘導(dǎo)阻力降低升阻比;在大攻角的情況下,小翼表面會(huì)出現(xiàn)渦的脫落,進(jìn)而影響機(jī)翼表面的氣流分布。

    2.1.2結(jié)構(gòu)參數(shù)對(duì)阻力的影響

    在不同翼展長(zhǎng)的情況下阻力隨攻角的變化曲線如圖13所示。計(jì)算結(jié)果表明,在相同的攻角下,阻力隨翼展長(zhǎng)的增加而增大。此外,在其它幾何參數(shù)確定時(shí),隨著翼展長(zhǎng)的增加,阻力的增加速率變大。因?yàn)闄C(jī)翼上下表面的壓強(qiáng)差產(chǎn)生了升力,并且由于氣流壓差的作用在翼尖附近會(huì)產(chǎn)生一個(gè)由下向上繞的誘導(dǎo)阻力,隨著翼展長(zhǎng)的增加,誘導(dǎo)阻力增大,因此增加了無人機(jī)的阻力。

    圖13 攻角與阻力的變化曲線Fig.13 Relationship curves between attack angle and drag

    圖14顯示了在不同的小翼厚時(shí),阻力隨后掠角的變化關(guān)系。結(jié)果表明隨著后掠角的增大,阻力先增大后減小。在指定小翼厚的情況下存在阻力的最大值。因?yàn)樵诤舐咏菫?5°~20°的變化過程中,機(jī)翼的后緣逐漸由后掠變?yōu)樗?,由于后緣氣流附著面積增加,導(dǎo)致阻力增加;在后掠角在20°~35°的變化過程中,隨著前緣后掠角的增加導(dǎo)致機(jī)翼壁面附近氣流的提前分離,無人機(jī)的阻力逐漸減小。此外,在較小的小翼厚時(shí),阻力的增加速率較慢。它可以解釋為當(dāng)小翼厚小于4 mm時(shí),小翼可以視作為薄片用于降低翼尖擾流,進(jìn)而減小誘導(dǎo)阻力;當(dāng)小翼厚大于4 mm時(shí),小翼可以視作翼展的延伸,并且在小翼厚方向上會(huì)生成新的誘導(dǎo)阻力,并且隨著小翼厚的增加,誘導(dǎo)阻力變大。結(jié)果表明,當(dāng)后掠角增加時(shí),阻力先增大后減小;當(dāng)小翼厚增加時(shí),阻力逐漸增加。

    圖14 后掠角與阻力的變化曲線Fig.14 Relationship curves between sweep angle and drag

    2.2 響應(yīng)面模型建立與精度驗(yàn)證

    2.2.1響應(yīng)面模型建立

    利用Design-Expert軟件,采用逐步回歸的方法構(gòu)建了翼展長(zhǎng)、后掠角、小翼高和小翼厚4個(gè)結(jié)構(gòu)參數(shù)與升阻比和阻力的二次響應(yīng)面方程

    (4)

    (5)

    根據(jù)響應(yīng)面方程對(duì)結(jié)構(gòu)參數(shù)進(jìn)行了敏感性分析,靈敏度越高,表明該結(jié)構(gòu)參數(shù)對(duì)性能的影響越顯著,結(jié)果如圖15所示。

    圖15 結(jié)構(gòu)參數(shù)敏感性分析結(jié)果Fig.15 Sensitivity of structural parameters

    在阻力方面翼展長(zhǎng)和小翼厚為正效應(yīng),后掠角和小翼高為負(fù)效應(yīng),并且在靈敏度上b>Λw≈h>lv。在升阻比方面,翼展長(zhǎng)和小翼高為正效應(yīng),后掠角和小翼厚為負(fù)效應(yīng);并且在靈敏度上b>h>Λw>lv。

    2.2.2響應(yīng)面模型精度評(píng)價(jià)及驗(yàn)證

    利用構(gòu)建的升阻比和阻力響應(yīng)面模型進(jìn)行了MOGA中個(gè)體的求解,因此需要對(duì)響應(yīng)面模型的精度和通用性進(jìn)行驗(yàn)證。用3個(gè)對(duì)照樣本點(diǎn)的CFD模擬值作為真值,并與響應(yīng)面預(yù)測(cè)值進(jìn)行對(duì)比,結(jié)果如表1所示。

    阻力的相對(duì)誤差隨著翼展長(zhǎng)的增加逐漸減小,升阻比的相對(duì)誤差隨著翼展長(zhǎng)的增加先減小后增大,3個(gè)對(duì)照點(diǎn)的相對(duì)誤差小于3%,表明響應(yīng)面模型具有較高的精度和良好的通用性,可用于尾座式無人機(jī)升阻比和阻力的預(yù)測(cè)。

    2.3 最優(yōu)結(jié)構(gòu)參數(shù)

    以升阻比取得最大值、阻力取得最小值為目標(biāo)函數(shù),利用MOGA在全局范圍內(nèi)尋找最優(yōu)結(jié)構(gòu)參數(shù)。采用篩選法構(gòu)建初始種群,利用響應(yīng)面模型進(jìn)行升阻比和阻力的計(jì)算,如圖16所示,得到了3組最優(yōu)的結(jié)構(gòu)參數(shù)如表2所示。其中最優(yōu)點(diǎn)1的升阻比提高了16.3%,阻力降低了4.1%。3個(gè)優(yōu)化模型的升阻比平均提高了14.2%,阻力降低了4.8%,均優(yōu)于原始樣機(jī)。對(duì)比3架樣機(jī)的面積可知,3號(hào)樣機(jī)機(jī)翼面積最小,在姿態(tài)轉(zhuǎn)換過程中能減小側(cè)向偏移量,因此本文選定最優(yōu)點(diǎn)3為最終模型。則尾座式無人機(jī)的結(jié)構(gòu)參數(shù)為:翼展長(zhǎng)1 123 mm、后掠角34°、小翼高39 mm、小翼厚3 mm。

    表1 響應(yīng)面模型精度驗(yàn)證Tab.1 Accuracy verification of response surface model

    圖16 多目標(biāo)遺傳算法尋優(yōu)結(jié)果Fig.16 MOGA optimization result

    參數(shù)原始參數(shù)值最優(yōu)點(diǎn)1最優(yōu)點(diǎn)2最優(yōu)點(diǎn)3翼展長(zhǎng)/mm10001193.91153.41123.3后掠角/(°)2034.5834.4034.37小翼高/mm3040.0839.7939.26小翼厚/mm63.013.223.21面積/m20.4740.5550.5390.526升阻比12.17114.15613.87113.680升阻比優(yōu)化比率/%16.314.012.4阻力/N0.6580.6310.6250.623阻力優(yōu)化比率/%-4.1-5.0-5.3

    3 結(jié)論

    (1)利用CFD對(duì)25組樣本點(diǎn)的升阻比和阻力進(jìn)行了數(shù)值模擬,樣本模擬數(shù)據(jù)、風(fēng)洞模型模擬數(shù)據(jù)和試驗(yàn)數(shù)據(jù)的對(duì)比分析表明,樣本與風(fēng)洞模型的相對(duì)誤差小于2.7%,風(fēng)洞模型與試驗(yàn)的相對(duì)誤差小于7.3%,樣本與試驗(yàn)的相對(duì)誤差小于8.0%,說明數(shù)值模擬方法可靠,可以用于尾座式無人機(jī)的外流場(chǎng)模擬。

    (2)構(gòu)建了結(jié)構(gòu)參數(shù)與升阻比和阻力的響應(yīng)面模型,模型相對(duì)誤差小于3%;基于響應(yīng)面模型的敏感性分析結(jié)果表明,在阻力方面翼展長(zhǎng)和小翼厚為正效應(yīng),后掠角和小翼高為負(fù)效應(yīng),其中翼展長(zhǎng)為主要影響因素;在升阻比方面,翼展長(zhǎng)和小翼高為正效應(yīng),后掠角和小翼厚為負(fù)效應(yīng),其中翼展長(zhǎng)和小翼厚為主要影響因素。

    (3)升阻比隨著翼展長(zhǎng)和小翼高的增加而增大,后掠角和小翼高對(duì)升阻比的影響較小,當(dāng)攻角為4°~8°時(shí),升阻比隨小翼厚的增加而減小,當(dāng)攻角為10°~12°時(shí),升阻比隨小翼厚的增加而增大;阻力隨著翼展長(zhǎng)和小翼厚的增加而增大,隨小翼高的增加而減小,隨后掠角的增加先增大后減小。

    (4)利用多目標(biāo)遺傳算法確定了尾座式無人機(jī)的最優(yōu)結(jié)構(gòu)參數(shù)為:翼展長(zhǎng)1 123 mm、后掠角34°、小翼高39 mm、小翼厚3 mm。與原始樣機(jī)相比升阻比提高了12.4%,阻力降低了5.3%。

    猜你喜歡
    后掠角尾座小翼
    漢字獵人(一)
    不同后掠角大展弦比復(fù)合材料機(jī)翼氣動(dòng)特性
    我家養(yǎng)了一只紙精靈(二)
    我家養(yǎng)了一只紙精靈(四)
    對(duì)故障尾座開裂原因的分析意見
    我是霸王龍
    新型重型臥式車床尾座鉆孔裝置設(shè)計(jì)
    基于CFD的最優(yōu)變后掠規(guī)律研究
    不同后掠角三角翼的靜態(tài)地面效應(yīng)數(shù)值模擬
    中頻感應(yīng)設(shè)備在槍管尾座熱處理中的應(yīng)用
    亚洲精品在线美女| 欧美黑人欧美精品刺激| 国产在视频线精品| 国产精品一区二区在线不卡| 99热网站在线观看| 肉色欧美久久久久久久蜜桃| 国产午夜精品一二区理论片| 国产av一区二区精品久久| 亚洲国产精品国产精品| 中文字幕亚洲精品专区| 如何舔出高潮| 麻豆乱淫一区二区| 中文字幕精品免费在线观看视频| 久久ye,这里只有精品| 成人免费观看视频高清| 国产又色又爽无遮挡免| 伊人久久国产一区二区| 水蜜桃什么品种好| 中文字幕精品免费在线观看视频| 亚洲国产欧美网| 最近中文字幕高清免费大全6| 欧美人与性动交α欧美精品济南到| 欧美日韩一区二区视频在线观看视频在线| 午夜福利乱码中文字幕| www.自偷自拍.com| av电影中文网址| 啦啦啦在线观看免费高清www| 在线天堂最新版资源| 国产野战对白在线观看| 建设人人有责人人尽责人人享有的| 丝袜美足系列| 熟女少妇亚洲综合色aaa.| 看非洲黑人一级黄片| 国产探花极品一区二区| 国产成人午夜福利电影在线观看| 午夜影院在线不卡| 一级a爱视频在线免费观看| 黄色毛片三级朝国网站| 欧美另类一区| 日韩大码丰满熟妇| 国产高清不卡午夜福利| 日韩 欧美 亚洲 中文字幕| 精品国产超薄肉色丝袜足j| 亚洲国产欧美网| 久久影院123| 国产精品秋霞免费鲁丝片| 母亲3免费完整高清在线观看| 久久亚洲国产成人精品v| 在线亚洲精品国产二区图片欧美| 久久精品亚洲熟妇少妇任你| 国产亚洲最大av| 精品一区二区免费观看| 交换朋友夫妻互换小说| 女人久久www免费人成看片| av视频免费观看在线观看| 九九爱精品视频在线观看| 国产老妇伦熟女老妇高清| 啦啦啦中文免费视频观看日本| 各种免费的搞黄视频| 国产精品成人在线| 亚洲欧美一区二区三区黑人| 人妻一区二区av| 一个人免费看片子| 三上悠亚av全集在线观看| 欧美人与性动交α欧美软件| 中国国产av一级| xxx大片免费视频| 又大又黄又爽视频免费| 日韩av在线免费看完整版不卡| 美女高潮到喷水免费观看| 久久婷婷青草| av线在线观看网站| 9热在线视频观看99| 少妇猛男粗大的猛烈进出视频| 日本猛色少妇xxxxx猛交久久| 别揉我奶头~嗯~啊~动态视频 | 精品国产露脸久久av麻豆| 国产av一区二区精品久久| 国产日韩欧美亚洲二区| 无遮挡黄片免费观看| 免费看av在线观看网站| 建设人人有责人人尽责人人享有的| 天堂8中文在线网| xxxhd国产人妻xxx| 日韩av免费高清视频| 国产成人啪精品午夜网站| 日本猛色少妇xxxxx猛交久久| 美女扒开内裤让男人捅视频| 爱豆传媒免费全集在线观看| 少妇的丰满在线观看| 亚洲av男天堂| 久久综合国产亚洲精品| 丝袜脚勾引网站| 一区二区日韩欧美中文字幕| 亚洲色图 男人天堂 中文字幕| 日本爱情动作片www.在线观看| 国产一级毛片在线| 国产精品av久久久久免费| 大码成人一级视频| 亚洲精品在线美女| 人体艺术视频欧美日本| 中文精品一卡2卡3卡4更新| 黄色视频在线播放观看不卡| 看免费成人av毛片| 日韩制服骚丝袜av| 又黄又粗又硬又大视频| 精品视频人人做人人爽| 成人手机av| 中文字幕av电影在线播放| 99久久人妻综合| 丝袜美腿诱惑在线| 男女国产视频网站| 纯流量卡能插随身wifi吗| 女人久久www免费人成看片| 女人精品久久久久毛片| 一级a爱视频在线免费观看| 人人妻人人爽人人添夜夜欢视频| 十分钟在线观看高清视频www| www日本在线高清视频| 午夜影院在线不卡| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 免费少妇av软件| 久久久久久久久久久久大奶| 久久 成人 亚洲| av线在线观看网站| 视频在线观看一区二区三区| 一二三四中文在线观看免费高清| 搡老岳熟女国产| 人成视频在线观看免费观看| 国产精品女同一区二区软件| 丝袜喷水一区| av在线观看视频网站免费| 午夜免费鲁丝| av在线老鸭窝| 一级片免费观看大全| 久久ye,这里只有精品| 亚洲av综合色区一区| 女人久久www免费人成看片| 国产激情久久老熟女| 午夜福利影视在线免费观看| 日韩一区二区视频免费看| 久久精品亚洲av国产电影网| 一级毛片我不卡| 欧美黑人精品巨大| 免费在线观看黄色视频的| 亚洲国产最新在线播放| 另类亚洲欧美激情| 一区二区av电影网| 免费黄色在线免费观看| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| 九九爱精品视频在线观看| 亚洲第一区二区三区不卡| svipshipincom国产片| 久久久精品区二区三区| 天天添夜夜摸| 久久影院123| 成年av动漫网址| 在线天堂最新版资源| 国产精品 欧美亚洲| 啦啦啦啦在线视频资源| 久久久久久久久久久免费av| 丰满乱子伦码专区| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 嫩草影院入口| 亚洲视频免费观看视频| 国产精品亚洲av一区麻豆 | 国产精品偷伦视频观看了| 热re99久久国产66热| 国产精品成人在线| 亚洲久久久国产精品| 制服丝袜香蕉在线| 观看av在线不卡| 最黄视频免费看| 一区二区av电影网| 久久午夜综合久久蜜桃| 午夜福利视频精品| 欧美日韩精品网址| 欧美成人精品欧美一级黄| 国产亚洲精品第一综合不卡| 国产淫语在线视频| 亚洲国产日韩一区二区| 国产成人精品久久久久久| 欧美日韩福利视频一区二区| 成人漫画全彩无遮挡| 侵犯人妻中文字幕一二三四区| 久久毛片免费看一区二区三区| 亚洲精品自拍成人| 日韩中文字幕视频在线看片| 国产片特级美女逼逼视频| 亚洲少妇的诱惑av| 国产成人午夜福利电影在线观看| 99re6热这里在线精品视频| 黄片小视频在线播放| 老司机影院毛片| 母亲3免费完整高清在线观看| av福利片在线| 国产亚洲最大av| 最近手机中文字幕大全| 亚洲av综合色区一区| 久久狼人影院| 九草在线视频观看| svipshipincom国产片| av电影中文网址| 人人妻,人人澡人人爽秒播 | 女人被躁到高潮嗷嗷叫费观| 欧美在线黄色| 精品国产露脸久久av麻豆| 成人国产av品久久久| 两个人看的免费小视频| 国产免费一区二区三区四区乱码| 超碰97精品在线观看| 亚洲国产看品久久| 免费观看人在逋| 看非洲黑人一级黄片| 搡老乐熟女国产| 亚洲精品日本国产第一区| 国产精品久久久久久精品古装| 91精品三级在线观看| 天天躁夜夜躁狠狠久久av| 精品一区二区三区四区五区乱码 | 亚洲av福利一区| 午夜激情av网站| 九草在线视频观看| 国产一区二区 视频在线| 天天躁狠狠躁夜夜躁狠狠躁| 99热国产这里只有精品6| 777久久人妻少妇嫩草av网站| 黄色一级大片看看| 青春草亚洲视频在线观看| 精品亚洲成a人片在线观看| 国产精品成人在线| 国产 精品1| 国产精品一区二区在线观看99| 黑人猛操日本美女一级片| 极品少妇高潮喷水抽搐| 国产黄频视频在线观看| 最近手机中文字幕大全| 十八禁网站网址无遮挡| 国产亚洲午夜精品一区二区久久| 午夜影院在线不卡| 精品国产乱码久久久久久小说| 一级a爱视频在线免费观看| 黑丝袜美女国产一区| 久久综合国产亚洲精品| 久久天堂一区二区三区四区| 色综合欧美亚洲国产小说| 晚上一个人看的免费电影| 亚洲国产欧美一区二区综合| 9191精品国产免费久久| 捣出白浆h1v1| 天堂8中文在线网| 国产淫语在线视频| 色婷婷av一区二区三区视频| 熟妇人妻不卡中文字幕| 这个男人来自地球电影免费观看 | 男女边摸边吃奶| 王馨瑶露胸无遮挡在线观看| 亚洲中文av在线| www.熟女人妻精品国产| 日韩视频在线欧美| 99热网站在线观看| 亚洲av在线观看美女高潮| 国产精品国产三级专区第一集| 久久精品亚洲熟妇少妇任你| 欧美激情高清一区二区三区 | 欧美激情高清一区二区三区 | 在线看a的网站| 啦啦啦在线免费观看视频4| 人妻 亚洲 视频| 亚洲精品乱久久久久久| 99久久人妻综合| 日本av免费视频播放| 亚洲国产精品成人久久小说| av.在线天堂| 啦啦啦在线免费观看视频4| kizo精华| 成人影院久久| 久久久精品免费免费高清| 七月丁香在线播放| 亚洲美女搞黄在线观看| 午夜av观看不卡| 大片免费播放器 马上看| 999精品在线视频| 热re99久久国产66热| 建设人人有责人人尽责人人享有的| 捣出白浆h1v1| 国产97色在线日韩免费| 亚洲一区中文字幕在线| 人人妻人人爽人人添夜夜欢视频| 一区在线观看完整版| 国产av国产精品国产| 天堂俺去俺来也www色官网| 十八禁网站网址无遮挡| 亚洲天堂av无毛| 黄色怎么调成土黄色| 亚洲精品国产区一区二| 99re6热这里在线精品视频| 国产极品天堂在线| 制服人妻中文乱码| 老司机靠b影院| 亚洲精品美女久久av网站| 日本色播在线视频| 啦啦啦在线免费观看视频4| 99久久人妻综合| 成年美女黄网站色视频大全免费| 一级片'在线观看视频| 国产在线一区二区三区精| 最黄视频免费看| 在线观看www视频免费| 久久天躁狠狠躁夜夜2o2o | 韩国高清视频一区二区三区| 在线观看免费日韩欧美大片| 丰满饥渴人妻一区二区三| 天天操日日干夜夜撸| 亚洲四区av| 中文字幕人妻丝袜一区二区 | 亚洲人成电影观看| 国产97色在线日韩免费| 午夜激情av网站| 亚洲在久久综合| 亚洲欧美成人综合另类久久久| 久久影院123| 啦啦啦啦在线视频资源| 亚洲精品美女久久av网站| 少妇被粗大猛烈的视频| 欧美乱码精品一区二区三区| 极品少妇高潮喷水抽搐| 午夜福利免费观看在线| e午夜精品久久久久久久| av女优亚洲男人天堂| 精品视频人人做人人爽| 亚洲国产欧美一区二区综合| 一级毛片电影观看| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 国产成人欧美在线观看 | 成年人免费黄色播放视频| 亚洲国产中文字幕在线视频| 欧美在线黄色| 欧美中文综合在线视频| 亚洲人成网站在线观看播放| 欧美另类一区| 亚洲成人国产一区在线观看 | 亚洲综合精品二区| 90打野战视频偷拍视频| 中文字幕高清在线视频| 人体艺术视频欧美日本| 亚洲自偷自拍图片 自拍| 夫妻性生交免费视频一级片| 美国免费a级毛片| 精品国产一区二区三区久久久樱花| 91精品国产国语对白视频| 国产探花极品一区二区| 国产成人免费无遮挡视频| 伦理电影大哥的女人| 久久 成人 亚洲| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产精品嫩草影院av在线观看| 亚洲欧美一区二区三区黑人| 国产 精品1| 久久久亚洲精品成人影院| 九色亚洲精品在线播放| 1024视频免费在线观看| 国产日韩欧美视频二区| 日本色播在线视频| 嫩草影视91久久| 精品国产乱码久久久久久小说| 女的被弄到高潮叫床怎么办| 国产国语露脸激情在线看| 免费黄网站久久成人精品| 人人妻人人澡人人看| 久久久久久久国产电影| 精品国产一区二区久久| 一个人免费看片子| 日韩av免费高清视频| 精品国产一区二区三区四区第35| 国产男女内射视频| 少妇被粗大的猛进出69影院| 国产精品女同一区二区软件| 九草在线视频观看| 丰满少妇做爰视频| 成人毛片60女人毛片免费| 五月天丁香电影| 国产成人精品无人区| 国产精品99久久99久久久不卡 | 熟妇人妻不卡中文字幕| 亚洲免费av在线视频| av在线app专区| 亚洲av男天堂| 日韩电影二区| 青春草亚洲视频在线观看| 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 日韩,欧美,国产一区二区三区| 久久 成人 亚洲| 亚洲人成77777在线视频| 国产在视频线精品| 久久久久久久久久久久大奶| 激情视频va一区二区三区| 亚洲熟女精品中文字幕| 亚洲精品一区蜜桃| 午夜激情久久久久久久| 少妇被粗大的猛进出69影院| 黄色怎么调成土黄色| 黄色一级大片看看| 日韩免费高清中文字幕av| 免费观看av网站的网址| 国产精品免费视频内射| 青草久久国产| 亚洲精品国产av蜜桃| 少妇人妻 视频| 国产乱人偷精品视频| 久久久久视频综合| 日韩一区二区视频免费看| 久久久久人妻精品一区果冻| 九草在线视频观看| 亚洲欧美清纯卡通| a级毛片在线看网站| 亚洲精品成人av观看孕妇| 如日韩欧美国产精品一区二区三区| 久久久久久久国产电影| 各种免费的搞黄视频| avwww免费| www.熟女人妻精品国产| 亚洲国产日韩一区二区| 一本久久精品| 一级,二级,三级黄色视频| 欧美成人午夜精品| 如日韩欧美国产精品一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲欧美成人综合另类久久久| 久久久国产欧美日韩av| 日本猛色少妇xxxxx猛交久久| 亚洲欧美成人精品一区二区| 久久国产亚洲av麻豆专区| 老司机影院毛片| www日本在线高清视频| 久久久久久久久久久久大奶| 蜜桃国产av成人99| 精品久久久久久电影网| 国产又爽黄色视频| 成人影院久久| 久久久国产精品麻豆| 国产精品久久久久成人av| 久久久久久久久久久久大奶| 欧美亚洲日本最大视频资源| 日韩制服骚丝袜av| 日韩一区二区视频免费看| 青青草视频在线视频观看| 久久精品人人爽人人爽视色| 丰满少妇做爰视频| 天天添夜夜摸| 午夜日韩欧美国产| 日韩,欧美,国产一区二区三区| 国产精品三级大全| 狠狠精品人妻久久久久久综合| 亚洲人成77777在线视频| 免费女性裸体啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 热99久久久久精品小说推荐| 日本猛色少妇xxxxx猛交久久| 国产午夜精品一二区理论片| 久久青草综合色| 波野结衣二区三区在线| 9热在线视频观看99| 街头女战士在线观看网站| 亚洲国产欧美日韩在线播放| 18禁裸乳无遮挡动漫免费视频| 赤兔流量卡办理| 欧美日韩视频高清一区二区三区二| 大香蕉久久网| 高清欧美精品videossex| 亚洲国产欧美一区二区综合| 中文字幕制服av| 桃花免费在线播放| 亚洲视频免费观看视频| 丝袜美足系列| a级毛片在线看网站| 如何舔出高潮| 一本色道久久久久久精品综合| 久久精品aⅴ一区二区三区四区| 午夜免费鲁丝| 丝袜美腿诱惑在线| 一边摸一边做爽爽视频免费| 亚洲国产av影院在线观看| 巨乳人妻的诱惑在线观看| 亚洲国产精品一区三区| 90打野战视频偷拍视频| 九色亚洲精品在线播放| 又大又黄又爽视频免费| 在线观看免费视频网站a站| 看免费av毛片| 满18在线观看网站| 悠悠久久av| 精品久久蜜臀av无| 亚洲久久久国产精品| 美国免费a级毛片| 又大又爽又粗| 黑人巨大精品欧美一区二区蜜桃| 多毛熟女@视频| 美女高潮到喷水免费观看| 国产精品久久久久成人av| av在线观看视频网站免费| 自线自在国产av| 亚洲国产中文字幕在线视频| 大香蕉久久网| 777久久人妻少妇嫩草av网站| 欧美日韩视频精品一区| 日本欧美视频一区| 人人澡人人妻人| 亚洲成人av在线免费| 777久久人妻少妇嫩草av网站| 中文字幕人妻熟女乱码| 男女无遮挡免费网站观看| 啦啦啦中文免费视频观看日本| 亚洲欧美精品自产自拍| 黄色 视频免费看| 国产精品亚洲av一区麻豆 | 中文字幕高清在线视频| 日韩av不卡免费在线播放| 免费黄网站久久成人精品| 午夜久久久在线观看| 国产一区二区三区综合在线观看| 免费在线观看视频国产中文字幕亚洲 | 国产精品熟女久久久久浪| 久久性视频一级片| 国产日韩欧美亚洲二区| 国产一区二区 视频在线| 亚洲国产日韩一区二区| 亚洲一区中文字幕在线| 99热国产这里只有精品6| 久久婷婷青草| 热re99久久国产66热| 亚洲一区二区三区欧美精品| 国产深夜福利视频在线观看| 免费av中文字幕在线| 国产精品一区二区在线观看99| 久久久久国产精品人妻一区二区| 日韩一卡2卡3卡4卡2021年| 午夜免费男女啪啪视频观看| 亚洲欧美一区二区三区久久| 伊人亚洲综合成人网| 久久人人97超碰香蕉20202| 男女高潮啪啪啪动态图| 最近中文字幕2019免费版| 一区二区三区激情视频| 亚洲欧美成人综合另类久久久| av在线app专区| 又粗又硬又长又爽又黄的视频| 99re6热这里在线精品视频| 我要看黄色一级片免费的| 美女大奶头黄色视频| 亚洲五月色婷婷综合| 99精品久久久久人妻精品| 午夜免费观看性视频| 国产成人免费无遮挡视频| 男女边摸边吃奶| 极品少妇高潮喷水抽搐| 天堂8中文在线网| 男女之事视频高清在线观看 | 新久久久久国产一级毛片| 黑人欧美特级aaaaaa片| 制服丝袜香蕉在线| 999久久久国产精品视频| 男女国产视频网站| 精品国产露脸久久av麻豆| 亚洲精品国产区一区二| 国产男女内射视频| 国产高清不卡午夜福利| 热re99久久精品国产66热6| 色综合欧美亚洲国产小说| 五月开心婷婷网| 日韩 亚洲 欧美在线| 国产成人免费观看mmmm| 精品人妻熟女毛片av久久网站| 亚洲精品日本国产第一区| 亚洲av成人精品一二三区| 亚洲人成网站在线观看播放| 亚洲精品国产av成人精品| 婷婷色综合大香蕉| h视频一区二区三区| 欧美日韩国产mv在线观看视频| 黄片无遮挡物在线观看| 国产免费又黄又爽又色| 丁香六月欧美| 亚洲人成电影观看| 日韩视频在线欧美| 99久久综合免费| 国产精品国产三级专区第一集| 天天影视国产精品| 精品少妇内射三级| 91aial.com中文字幕在线观看| 日韩精品免费视频一区二区三区| 好男人视频免费观看在线| 久久性视频一级片| h视频一区二区三区| 精品久久久精品久久久| 国产精品无大码| 久久精品aⅴ一区二区三区四区| 一级毛片我不卡| 美女主播在线视频| 国产不卡av网站在线观看| 一本大道久久a久久精品| 国产熟女午夜一区二区三区| 午夜免费男女啪啪视频观看| 99久国产av精品国产电影| 久久精品久久久久久噜噜老黄| 熟女av电影| 9热在线视频观看99| 成人国产av品久久久|