• <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)用
    国产一区有黄有色的免费视频| 日本wwww免费看| 18禁国产床啪视频网站| 新久久久久国产一级毛片| 免费观看性生交大片5| 一级毛片 在线播放| 日本色播在线视频| 国产男人的电影天堂91| 有码 亚洲区| 制服人妻中文乱码| 亚洲四区av| 最近2019中文字幕mv第一页| 国产野战对白在线观看| 欧美最新免费一区二区三区| 观看av在线不卡| 成人黄色视频免费在线看| 国产人伦9x9x在线观看 | av不卡在线播放| 丝瓜视频免费看黄片| 一区二区日韩欧美中文字幕| 国产一区二区激情短视频 | 亚洲精品av麻豆狂野| 99久国产av精品国产电影| 蜜桃国产av成人99| 91成人精品电影| 国产毛片在线视频| 免费观看a级毛片全部| av又黄又爽大尺度在线免费看| 在线观看美女被高潮喷水网站| 曰老女人黄片| 国产精品久久久av美女十八| 成人免费观看视频高清| 男女边摸边吃奶| 菩萨蛮人人尽说江南好唐韦庄| 777米奇影视久久| 一区二区三区乱码不卡18| 免费黄色在线免费观看| 国产色婷婷99| 精品国产乱码久久久久久男人| 日韩一本色道免费dvd| 免费久久久久久久精品成人欧美视频| 91精品伊人久久大香线蕉| 久久久久国产精品人妻一区二区| a级片在线免费高清观看视频| 亚洲第一区二区三区不卡| 波多野结衣一区麻豆| 日韩成人av中文字幕在线观看| 五月天丁香电影| 精品99又大又爽又粗少妇毛片| 亚洲在久久综合| 天天躁夜夜躁狠狠久久av| 亚洲伊人久久精品综合| 久久ye,这里只有精品| 亚洲五月色婷婷综合| 十八禁高潮呻吟视频| 黄片小视频在线播放| 欧美xxⅹ黑人| 99热网站在线观看| 亚洲四区av| 男人舔女人的私密视频| 十八禁网站网址无遮挡| 天天躁夜夜躁狠狠久久av| 欧美日韩亚洲国产一区二区在线观看 | 免费观看av网站的网址| 欧美bdsm另类| 在线观看国产h片| 久久精品夜色国产| 精品亚洲乱码少妇综合久久| 日本-黄色视频高清免费观看| 99香蕉大伊视频| 又大又黄又爽视频免费| 夜夜骑夜夜射夜夜干| 国产精品久久久久久久久免| 午夜免费男女啪啪视频观看| 中文欧美无线码| 中文精品一卡2卡3卡4更新| 免费观看无遮挡的男女| 亚洲综合精品二区| 久久精品国产鲁丝片午夜精品| 人妻一区二区av| 国产在线视频一区二区| 精品久久久久久电影网| 热re99久久国产66热| 9热在线视频观看99| 丝袜喷水一区| 美女脱内裤让男人舔精品视频| 一区在线观看完整版| xxxhd国产人妻xxx| 免费看av在线观看网站| 考比视频在线观看| av免费观看日本| 久久久精品94久久精品| 最近2019中文字幕mv第一页| 老汉色av国产亚洲站长工具| 男人操女人黄网站| 久久人人爽av亚洲精品天堂| 亚洲精品一区蜜桃| 人人妻人人澡人人爽人人夜夜| 如日韩欧美国产精品一区二区三区| 久久人妻熟女aⅴ| 啦啦啦中文免费视频观看日本| 亚洲欧洲国产日韩| 国产男人的电影天堂91| 超碰成人久久| 9191精品国产免费久久| 夜夜骑夜夜射夜夜干| 天堂俺去俺来也www色官网| 国产1区2区3区精品| 精品少妇久久久久久888优播| 王馨瑶露胸无遮挡在线观看| 午夜福利视频精品| 99久久精品国产国产毛片| 9191精品国产免费久久| 夜夜骑夜夜射夜夜干| 精品少妇一区二区三区视频日本电影 | 国产精品欧美亚洲77777| 亚洲四区av| 国产亚洲精品第一综合不卡| 亚洲中文av在线| 桃花免费在线播放| 最近手机中文字幕大全| 欧美av亚洲av综合av国产av | 男女啪啪激烈高潮av片| 国产精品免费大片| 热99国产精品久久久久久7| 久热这里只有精品99| 欧美xxⅹ黑人| 777米奇影视久久| 少妇人妻精品综合一区二区| 国产毛片在线视频| 久久久久精品久久久久真实原创| a级毛片黄视频| 久久精品国产a三级三级三级| 亚洲国产日韩一区二区| 这个男人来自地球电影免费观看 | 两个人看的免费小视频| 建设人人有责人人尽责人人享有的| 欧美人与善性xxx| 国产黄频视频在线观看| 校园人妻丝袜中文字幕| 成人二区视频| 成年人免费黄色播放视频| 欧美精品一区二区免费开放| 久久青草综合色| 激情视频va一区二区三区| 美女高潮到喷水免费观看| 亚洲天堂av无毛| 国产精品女同一区二区软件| 蜜桃国产av成人99| 满18在线观看网站| 夫妻午夜视频| 亚洲欧美一区二区三区久久| 国产毛片在线视频| 国产国语露脸激情在线看| 久久韩国三级中文字幕| 高清在线视频一区二区三区| 成年美女黄网站色视频大全免费| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美网| 在线观看www视频免费| 欧美黄色片欧美黄色片| 少妇人妻久久综合中文| 超碰成人久久| 久久久久久久亚洲中文字幕| 日本午夜av视频| 在线观看美女被高潮喷水网站| 在线观看www视频免费| 久久精品熟女亚洲av麻豆精品| 国产精品偷伦视频观看了| 91久久精品国产一区二区三区| 最近2019中文字幕mv第一页| 国产男女超爽视频在线观看| 亚洲国产欧美在线一区| 亚洲欧美中文字幕日韩二区| 在线天堂中文资源库| 看非洲黑人一级黄片| 91在线精品国自产拍蜜月| 伊人久久大香线蕉亚洲五| 岛国毛片在线播放| 99久久中文字幕三级久久日本| 啦啦啦中文免费视频观看日本| 丝袜喷水一区| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 久久久久久人妻| 欧美国产精品一级二级三级| 丝袜脚勾引网站| 一区福利在线观看| 亚洲精品aⅴ在线观看| 啦啦啦在线观看免费高清www| av不卡在线播放| 亚洲国产看品久久| 夫妻性生交免费视频一级片| 国产97色在线日韩免费| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 纵有疾风起免费观看全集完整版| 亚洲精品乱久久久久久| 亚洲精品第二区| 1024香蕉在线观看| 美女大奶头黄色视频| 免费播放大片免费观看视频在线观看| a级毛片黄视频| 人妻 亚洲 视频| 如日韩欧美国产精品一区二区三区| 精品少妇久久久久久888优播| 久久久久久人人人人人| 女性生殖器流出的白浆| 欧美中文综合在线视频| 少妇熟女欧美另类| 国产成人一区二区在线| www.精华液| 国产日韩欧美视频二区| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久 | av又黄又爽大尺度在线免费看| 老司机影院成人| 少妇被粗大的猛进出69影院| 69精品国产乱码久久久| 99久久精品国产国产毛片| 欧美黄色片欧美黄色片| 看非洲黑人一级黄片| 午夜91福利影院| 国产又爽黄色视频| 亚洲精品久久久久久婷婷小说| 亚洲国产精品一区三区| 久久精品夜色国产| 美女高潮到喷水免费观看| 亚洲一码二码三码区别大吗| √禁漫天堂资源中文www| 麻豆精品久久久久久蜜桃| 不卡视频在线观看欧美| 国产精品亚洲av一区麻豆 | 激情五月婷婷亚洲| 久久免费观看电影| 亚洲av电影在线进入| 欧美精品一区二区大全| 久久久久精品久久久久真实原创| 亚洲熟女精品中文字幕| 日日撸夜夜添| av一本久久久久| 亚洲欧美一区二区三区国产| 国产亚洲精品第一综合不卡| 在线精品无人区一区二区三| 美女大奶头黄色视频| 亚洲内射少妇av| 99热全是精品| 久久久国产精品麻豆| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 成人亚洲欧美一区二区av| 国产高清国产精品国产三级| 天天躁夜夜躁狠狠躁躁| 一级片免费观看大全| 亚洲婷婷狠狠爱综合网| 在线观看免费视频网站a站| 国产极品天堂在线| 国产成人av激情在线播放| 国产精品香港三级国产av潘金莲 | 黑人巨大精品欧美一区二区蜜桃| 母亲3免费完整高清在线观看 | 建设人人有责人人尽责人人享有的| 伦理电影大哥的女人| 成人国语在线视频| 久久久精品区二区三区| 人成视频在线观看免费观看| 天堂8中文在线网| 婷婷成人精品国产| 电影成人av| 人人妻人人添人人爽欧美一区卜| 丝袜喷水一区| 美女国产视频在线观看| 日韩中文字幕欧美一区二区 | 91在线精品国自产拍蜜月| 国产亚洲一区二区精品| 十八禁高潮呻吟视频| 黄色一级大片看看| 久久99精品国语久久久| 免费久久久久久久精品成人欧美视频| 又黄又粗又硬又大视频| 久久亚洲国产成人精品v| 大码成人一级视频| 婷婷成人精品国产| 成人毛片60女人毛片免费| 午夜福利视频精品| 久久毛片免费看一区二区三区| 大香蕉久久成人网| 黄色视频在线播放观看不卡| 性高湖久久久久久久久免费观看| 97在线视频观看| 波多野结衣av一区二区av| 亚洲一区二区三区欧美精品| 日韩 亚洲 欧美在线| 免费女性裸体啪啪无遮挡网站| 少妇猛男粗大的猛烈进出视频| 黄片小视频在线播放| 十八禁网站网址无遮挡| 久久ye,这里只有精品| 国产又爽黄色视频| 精品久久久久久电影网| 亚洲伊人色综图| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡 | av免费观看日本| 精品少妇一区二区三区视频日本电影 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品秋霞免费鲁丝片| 国产精品国产三级国产专区5o| 国产福利在线免费观看视频| 久久国内精品自在自线图片| 最近中文字幕2019免费版| 亚洲,一卡二卡三卡| 国产毛片在线视频| 久久久久久久久免费视频了| 高清视频免费观看一区二区| 青春草亚洲视频在线观看| 免费少妇av软件| 免费观看a级毛片全部| 国产在线一区二区三区精| 亚洲熟女精品中文字幕| 欧美精品亚洲一区二区| 久久精品人人爽人人爽视色| 啦啦啦视频在线资源免费观看| 一个人免费看片子| 日本色播在线视频| 国产一级毛片在线| 交换朋友夫妻互换小说| 美女国产高潮福利片在线看| 国产白丝娇喘喷水9色精品| 又粗又硬又长又爽又黄的视频| 亚洲少妇的诱惑av| 中文欧美无线码| 久久国内精品自在自线图片| 大陆偷拍与自拍| 亚洲精品国产色婷婷电影| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 国产成人aa在线观看| 波野结衣二区三区在线| 高清av免费在线| 亚洲精品在线美女| 精品人妻偷拍中文字幕| 一级毛片我不卡| 成人亚洲精品一区在线观看| 国产一区二区 视频在线| 亚洲图色成人| 性少妇av在线| 另类精品久久| 国产成人精品久久二区二区91 | 天天躁夜夜躁狠狠躁躁| 男的添女的下面高潮视频| 99久久综合免费| 午夜福利一区二区在线看| 国产亚洲欧美精品永久| 亚洲av在线观看美女高潮| 电影成人av| 人妻少妇偷人精品九色| 伦精品一区二区三区| 777米奇影视久久| 亚洲国产欧美日韩在线播放| 成人毛片60女人毛片免费| 天堂8中文在线网| 国产乱人偷精品视频| 人人妻人人澡人人看| 一边摸一边做爽爽视频免费| 午夜福利,免费看| 国产片内射在线| 亚洲少妇的诱惑av| 欧美日韩视频精品一区| 人妻系列 视频| 黄片小视频在线播放| 精品一区二区三区四区五区乱码 | 久久精品国产a三级三级三级| 最新中文字幕久久久久| 婷婷色av中文字幕| 最近最新中文字幕大全免费视频 | 国产极品天堂在线| 在线 av 中文字幕| 国产av码专区亚洲av| 国产成人a∨麻豆精品| 免费在线观看完整版高清| av又黄又爽大尺度在线免费看| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 国产精品久久久久久精品电影小说| 国产成人精品在线电影| 国产综合精华液| h视频一区二区三区| 婷婷色麻豆天堂久久| 中文字幕最新亚洲高清| 久久精品国产鲁丝片午夜精品| 青草久久国产| 中文字幕制服av| 欧美日韩成人在线一区二区| 九色亚洲精品在线播放| 亚洲激情五月婷婷啪啪| 色婷婷久久久亚洲欧美| 欧美激情高清一区二区三区 | 久久精品久久久久久久性| 成人毛片60女人毛片免费| 精品国产国语对白av| 美女脱内裤让男人舔精品视频| 亚洲 欧美一区二区三区| 18+在线观看网站| av免费在线看不卡| 中文字幕制服av| 亚洲精品,欧美精品| 久久国内精品自在自线图片| 99精国产麻豆久久婷婷| 亚洲欧洲精品一区二区精品久久久 | 午夜久久久在线观看| 女人高潮潮喷娇喘18禁视频| 少妇被粗大猛烈的视频| 一二三四在线观看免费中文在| 男女边摸边吃奶| 一本大道久久a久久精品| 免费大片黄手机在线观看| 精品少妇内射三级| 国产成人精品福利久久| 最近中文字幕高清免费大全6| 欧美激情极品国产一区二区三区| 丰满迷人的少妇在线观看| 成年动漫av网址| 91久久精品国产一区二区三区| 日韩电影二区| 在线观看三级黄色| 久久毛片免费看一区二区三区| 久久久国产一区二区| 免费黄频网站在线观看国产| 熟女少妇亚洲综合色aaa.| 国产亚洲精品第一综合不卡| 精品一区二区免费观看| 精品亚洲乱码少妇综合久久| 日韩制服骚丝袜av| 丝瓜视频免费看黄片| 日日撸夜夜添| 亚洲少妇的诱惑av| 精品一区二区免费观看| 国产爽快片一区二区三区| 黄片无遮挡物在线观看| 秋霞在线观看毛片| 看非洲黑人一级黄片| 老汉色∧v一级毛片| 国产精品国产三级专区第一集| 在线免费观看不下载黄p国产| 久久99一区二区三区| 亚洲婷婷狠狠爱综合网| 成人免费观看视频高清| 亚洲综合精品二区| 亚洲成色77777| 久久人人爽av亚洲精品天堂| 赤兔流量卡办理| 高清av免费在线| 日韩视频在线欧美| 99国产综合亚洲精品| av一本久久久久| 日本爱情动作片www.在线观看| 日本vs欧美在线观看视频| 国产一区二区 视频在线| 丰满饥渴人妻一区二区三| 人妻一区二区av| av福利片在线| 99久久精品国产国产毛片| 精品国产一区二区三区久久久樱花| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 99热国产这里只有精品6| 中文字幕制服av| 国产熟女午夜一区二区三区| 麻豆乱淫一区二区| 中文字幕人妻熟女乱码| 纯流量卡能插随身wifi吗| 亚洲天堂av无毛| 建设人人有责人人尽责人人享有的| av福利片在线| 狠狠精品人妻久久久久久综合| 欧美日韩精品成人综合77777| 久久久久国产精品人妻一区二区| a级毛片在线看网站| 女性被躁到高潮视频| 国产福利在线免费观看视频| 免费日韩欧美在线观看| 亚洲av日韩在线播放| 久久99热这里只频精品6学生| 1024香蕉在线观看| 制服人妻中文乱码| 亚洲,欧美精品.| 久久韩国三级中文字幕| 夫妻性生交免费视频一级片| 亚洲国产精品成人久久小说| 日本色播在线视频| 两个人免费观看高清视频| 国产精品免费大片| 热99久久久久精品小说推荐| 在线观看美女被高潮喷水网站| 黑丝袜美女国产一区| 国产精品99久久99久久久不卡 | 最近最新中文字幕免费大全7| 午夜免费观看性视频| 亚洲综合色网址| 欧美精品高潮呻吟av久久| 日韩大片免费观看网站| 国产乱人偷精品视频| 国产国语露脸激情在线看| 999精品在线视频| 亚洲精品,欧美精品| 人妻系列 视频| 亚洲人成电影观看| 亚洲欧美日韩另类电影网站| 婷婷色麻豆天堂久久| 99热全是精品| 极品人妻少妇av视频| 精品人妻在线不人妻| 啦啦啦在线免费观看视频4| 永久网站在线| 蜜桃国产av成人99| av又黄又爽大尺度在线免费看| 人妻一区二区av| 国语对白做爰xxxⅹ性视频网站| 视频在线观看一区二区三区| 男人舔女人的私密视频| 亚洲欧洲日产国产| 久久久久国产网址| 五月天丁香电影| 一区在线观看完整版| 国产成人精品一,二区| 成人漫画全彩无遮挡| 亚洲成人av在线免费| 日韩 亚洲 欧美在线| 2022亚洲国产成人精品| 日韩制服丝袜自拍偷拍| 日本猛色少妇xxxxx猛交久久| 九草在线视频观看| a 毛片基地| 免费在线观看黄色视频的| 成年动漫av网址| 国精品久久久久久国模美| 边亲边吃奶的免费视频| 欧美日韩亚洲高清精品| 国产激情久久老熟女| 岛国毛片在线播放| 成年人午夜在线观看视频| 欧美日韩一级在线毛片| 日韩一区二区视频免费看| 亚洲色图综合在线观看| 欧美精品一区二区大全| www.自偷自拍.com| 亚洲精品国产色婷婷电影| 成人亚洲欧美一区二区av| 国产乱来视频区| 色婷婷av一区二区三区视频| 宅男免费午夜| 又黄又粗又硬又大视频| 国产 一区精品| 久久久久久久亚洲中文字幕| 久久精品久久精品一区二区三区| 亚洲一码二码三码区别大吗| 狠狠精品人妻久久久久久综合| 色哟哟·www| 精品酒店卫生间| 777久久人妻少妇嫩草av网站| 一区二区av电影网| 亚洲欧美成人精品一区二区| 美国免费a级毛片| 欧美精品一区二区免费开放| 日韩视频在线欧美| 美女福利国产在线| 1024视频免费在线观看| 精品亚洲成国产av| 久久久久久久精品精品| 国产女主播在线喷水免费视频网站| 三上悠亚av全集在线观看| 国产精品 欧美亚洲| 国产高清不卡午夜福利| 久久久精品区二区三区| 啦啦啦在线观看免费高清www| 亚洲欧美成人精品一区二区| 中文字幕av电影在线播放| 欧美最新免费一区二区三区| 亚洲国产欧美网| 丝袜美腿诱惑在线| 在现免费观看毛片| 日韩人妻精品一区2区三区| 人妻系列 视频| 欧美bdsm另类| 少妇精品久久久久久久| 一边亲一边摸免费视频| 免费看不卡的av| 天堂8中文在线网| 黄频高清免费视频| 国产精品蜜桃在线观看| 99久久中文字幕三级久久日本| 久久精品熟女亚洲av麻豆精品| 黑人欧美特级aaaaaa片| 美女高潮到喷水免费观看| 一本久久精品| 亚洲av中文av极速乱| 久久久a久久爽久久v久久| 亚洲色图综合在线观看| 午夜日本视频在线| 久久久久久久亚洲中文字幕| 国产成人免费无遮挡视频| 汤姆久久久久久久影院中文字幕| 在线观看一区二区三区激情| 一二三四中文在线观看免费高清| 男女边摸边吃奶| 欧美精品高潮呻吟av久久| 国产1区2区3区精品| 免费高清在线观看视频在线观看| 久久久久久久国产电影| 国产精品人妻久久久影院| 午夜免费男女啪啪视频观看| 国产av国产精品国产| 香蕉丝袜av| 日韩一区二区三区影片|