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

    低速風(fēng)洞雙自由度大幅振蕩試驗(yàn)技術(shù)

    2016-11-14 00:58:18劉春明趙志軍卜忱王建鋒牟偉強(qiáng)
    航空學(xué)報(bào) 2016年8期
    關(guān)鍵詞:氣動(dòng)力迎角風(fēng)洞

    劉春明, 趙志軍, 卜忱, 王建鋒, 牟偉強(qiáng)

    1.中國航空工業(yè)空氣動(dòng)力研究院, 哈爾濱 150001 2.中國空氣動(dòng)力研究與發(fā)展中心, 綿陽 621000

    ?

    低速風(fēng)洞雙自由度大幅振蕩試驗(yàn)技術(shù)

    劉春明1,*, 趙志軍2, 卜忱1, 王建鋒1, 牟偉強(qiáng)1

    1.中國航空工業(yè)空氣動(dòng)力研究院, 哈爾濱150001 2.中國空氣動(dòng)力研究與發(fā)展中心, 綿陽621000

    為滿足國內(nèi)戰(zhàn)斗機(jī)型號(hào)研制對(duì)大迎角過失速機(jī)動(dòng)過程中非定常氣動(dòng)力問題研究的需求,中國航空工業(yè)空氣動(dòng)力研究院在FL-8風(fēng)洞開發(fā)了一套基于電液伺服馬達(dá)和電機(jī)耦合驅(qū)動(dòng)的雙自由度大幅振蕩試驗(yàn)技術(shù)。簡單介紹了該系統(tǒng)的基本組成和試驗(yàn)原理,給出了FL-8風(fēng)洞動(dòng)態(tài)標(biāo)準(zhǔn)模型雙自由度大幅振蕩試驗(yàn)的典型結(jié)果,分析表明:該系統(tǒng)運(yùn)行性能穩(wěn)定,試驗(yàn)數(shù)據(jù)可靠,重復(fù)性精度較高,可以有效應(yīng)用于飛行器雙自由度耦合運(yùn)動(dòng)氣動(dòng)特性研究。

    大迎角; 大振幅; 雙自由度; 非定常氣動(dòng)力; 風(fēng)洞試驗(yàn)

    高機(jī)動(dòng)性、高敏捷性和過失速機(jī)動(dòng)能力已經(jīng)成為現(xiàn)代及下一代戰(zhàn)斗機(jī)的戰(zhàn)術(shù)性能和作戰(zhàn)效能的重要指標(biāo)。過失速機(jī)動(dòng)可有效規(guī)避敵方攻擊,具有很高的戰(zhàn)術(shù)價(jià)值。一些先進(jìn)戰(zhàn)斗機(jī)已經(jīng)具有大迎角、過失速可控飛行能力,能夠完成各種高難度的機(jī)動(dòng)飛行,如蘇-27的眼鏡蛇機(jī)動(dòng)、X-31與F-22的Herbst機(jī)動(dòng)等[1-4]。復(fù)雜的機(jī)動(dòng)飛行過程中,飛機(jī)的飛行包線已經(jīng)擴(kuò)大到大迎角、高速率和大振幅區(qū)域;在大迎角、過失速飛行時(shí)產(chǎn)生的三維分離流動(dòng)、旋渦的相互干擾、渦的破裂、流動(dòng)滯后及流動(dòng)再附等使得在該區(qū)域內(nèi)氣動(dòng)力呈現(xiàn)出高度的非線性特性和非定常特性[5-9]。精確地獲知非線性、非定常氣動(dòng)力載荷對(duì)于飛機(jī)大迎角、過失速區(qū)域的氣動(dòng)力特性的研究、飛行力學(xué)特性的分析和飛行控制系統(tǒng)的設(shè)計(jì)極為重要。

    為了開展大迎角非定常氣動(dòng)力問題研究,世界各航空發(fā)達(dá)國家都開發(fā)了相應(yīng)的風(fēng)洞動(dòng)態(tài)試驗(yàn)技術(shù),用來進(jìn)行飛機(jī)各種運(yùn)動(dòng)模態(tài)的模擬及機(jī)動(dòng)飛行時(shí)間歷程的模擬,開展針對(duì)飛行器機(jī)動(dòng)性、敏捷性,特別是大迎角、過失速機(jī)動(dòng)和可控飛行能力的研究,為新機(jī)研制提供了有效的動(dòng)態(tài)試驗(yàn)手段[10]。試驗(yàn)技術(shù)從一些比較簡單的模型姿態(tài)變化和流速變化,如等速俯仰、俯仰角正弦振蕩、加減流速或流速正弦波動(dòng)等,到一些比較復(fù)雜的動(dòng)態(tài)試驗(yàn)裝置和技術(shù),如多軸試驗(yàn)系統(tǒng)(Multi-Axis Test Rig)、模型定位系統(tǒng)(Model Positioning Mechanism)、振蕩圓錐運(yùn)動(dòng)試驗(yàn)技術(shù)、機(jī)翼搖滾試驗(yàn)技術(shù)、氣動(dòng)力模型風(fēng)洞驗(yàn)證技術(shù)、大迎角大振幅多自由度耦合試驗(yàn)技術(shù)以及旋轉(zhuǎn)與振蕩的耦合運(yùn)動(dòng)試驗(yàn)技術(shù)[11-18]等。

    本文介紹了中國航空工業(yè)空氣動(dòng)力研究院研發(fā)的一種新型的低速風(fēng)洞雙自由度大幅振蕩試驗(yàn)技術(shù),包括系統(tǒng)組成、試驗(yàn)原理及動(dòng)態(tài)標(biāo)準(zhǔn)模型驗(yàn)證試驗(yàn)典型結(jié)果分析。

    1 雙自由度大幅振蕩試驗(yàn)系統(tǒng)

    FL-8風(fēng)洞雙自由度大幅振蕩試驗(yàn)系統(tǒng)可驅(qū)動(dòng)模型進(jìn)行諧波振蕩或任意給定軌跡的運(yùn)動(dòng),包括大幅振蕩試驗(yàn)、動(dòng)導(dǎo)數(shù)試驗(yàn)和典型機(jī)動(dòng)歷程模擬試驗(yàn)等。系統(tǒng)適用于1 m量級(jí)模型,試驗(yàn)風(fēng)速范圍為20~45 m/s。主要運(yùn)動(dòng)指標(biāo)見表1。

    表1雙自由度大幅振蕩試驗(yàn)系統(tǒng)運(yùn)動(dòng)指標(biāo)

    Table 1Motion indexs of double degree-of-freedom large amplitude oscillation test system

    ModelMotionindexInnerframefrequencyresponse/(Hz,(°))2,45Outerframefrequencyresponse/(Hz,(°))1,45Innerframeangularvelocity/((°)·s-1)560Outerframeangularvelocity/((°)·s-1)280Innerframeangularacceleration/((°)·s-2)3500Outerframeangularacceleration/((°)·s-2)1700

    1.1系統(tǒng)組成

    FL-8風(fēng)洞雙自由度大幅振蕩試驗(yàn)系統(tǒng)由運(yùn)動(dòng)系統(tǒng)和測(cè)控系統(tǒng)組成。

    雙自由度大幅振蕩試驗(yàn)系統(tǒng)如圖1所示。機(jī)械系統(tǒng)整體結(jié)構(gòu)采用內(nèi)、外框結(jié)構(gòu),內(nèi)框?yàn)闈L轉(zhuǎn)軸,外框?yàn)楦┭鲚S/偏航軸的框架結(jié)構(gòu),兩框均為單支撐形式,外框液壓馬達(dá)固定在支撐架上。內(nèi)框由電機(jī)驅(qū)動(dòng)滾轉(zhuǎn)支桿帶動(dòng)模型運(yùn)動(dòng)。當(dāng)模型采用尾撐形式時(shí),內(nèi)框機(jī)構(gòu)帶動(dòng)模型繞模型縱軸滾轉(zhuǎn)運(yùn)動(dòng)。外框由液壓馬達(dá)驅(qū)動(dòng)碳纖維彎刀帶動(dòng)模型運(yùn)動(dòng)。當(dāng)模型正裝時(shí)(如圖1所示狀態(tài)),外框機(jī)構(gòu)帶動(dòng)模型繞模型橫軸俯仰運(yùn)動(dòng)。當(dāng)模型側(cè)裝時(shí)(圖1中模型繞機(jī)體縱軸旋轉(zhuǎn)90°),外框機(jī)構(gòu)帶動(dòng)模型繞模型豎軸偏航運(yùn)動(dòng)。在設(shè)定好的運(yùn)動(dòng)軌跡下,同步驅(qū)動(dòng)內(nèi)、外框,可實(shí)現(xiàn)模型繞體軸的雙自由度運(yùn)動(dòng),包括俯仰與滾轉(zhuǎn)耦合運(yùn)動(dòng)、偏航與滾轉(zhuǎn)耦合運(yùn)動(dòng)。

    測(cè)控系統(tǒng)包括運(yùn)動(dòng)控制子系統(tǒng)和輔助信號(hào)子系統(tǒng),測(cè)控系統(tǒng)的結(jié)構(gòu)組成如圖2所示??刂葡到y(tǒng)由外框俯仰(或偏航)運(yùn)動(dòng)控制和內(nèi)框滾轉(zhuǎn)運(yùn)動(dòng)控制兩部分組成,分別模擬飛機(jī)在飛行過程中的滾動(dòng)、俯仰(或偏航)運(yùn)動(dòng)姿態(tài)角的變化。

    圖1 雙自由度大幅振蕩試驗(yàn)系統(tǒng)示意圖Fig.1 Schematic diagram of double degree-of-freedom large amplitude oscillation test system

    圖2 測(cè)控系統(tǒng)的組成Fig.2 Components of measurement and control system

    1.2系統(tǒng)雙軸同步控制原理與策略

    滾轉(zhuǎn)軸和俯仰軸(或偏航軸)分別閉環(huán)控制,實(shí)現(xiàn)模型雙自由度運(yùn)動(dòng)的模擬。該控制系統(tǒng)具有較高的響應(yīng)速度和定時(shí)精度,能夠定時(shí)觸發(fā)采集系統(tǒng)采集天平數(shù)據(jù),并且同步輸出運(yùn)動(dòng)位移信號(hào), 圖3給出了雙軸同步運(yùn)動(dòng)控制策略框圖。

    圖3 雙軸同步運(yùn)動(dòng)控制策略Fig.3 Control strategy for dual axis synchronizingmotion

    為保證兩軸同步運(yùn)動(dòng),采用實(shí)時(shí)系統(tǒng)同步驅(qū)動(dòng)液壓軸及電機(jī)軸,二者循環(huán)周期均為0.4 ms,在循環(huán)周期內(nèi)二者同時(shí)獨(dú)立讀取編碼器信息,計(jì)算實(shí)時(shí)位置輸入及輸出;將液壓軸的位置輸出作為電機(jī)軸的位置輸入,即電機(jī)軸跟蹤液壓軸運(yùn)動(dòng)。通過對(duì)采集回來的兩軸運(yùn)動(dòng)位置信息進(jìn)行曲線擬合,得到兩軸運(yùn)動(dòng)相位差值在0.03° 左右。

    1.3大幅振蕩試驗(yàn)數(shù)據(jù)采集與處理方法

    為了提高測(cè)量精度,每一個(gè)試驗(yàn)狀態(tài)的天平測(cè)量值都取若干周期運(yùn)動(dòng)的平均值。一般每個(gè)試驗(yàn)點(diǎn)采集10個(gè)振蕩周期的數(shù)據(jù),每一個(gè)周期等時(shí)間間隔采集240個(gè)數(shù)據(jù),取10個(gè)周期對(duì)應(yīng)點(diǎn)的平均值作為一個(gè)試驗(yàn)狀態(tài)的數(shù)據(jù)。吹風(fēng)試驗(yàn)前需在無風(fēng)狀態(tài)下測(cè)取模型的重力、慣性力和力矩,將吹風(fēng)試驗(yàn)所得的原始數(shù)據(jù)減去無風(fēng)時(shí)的數(shù)據(jù),然后轉(zhuǎn)換成氣動(dòng)系數(shù)。試驗(yàn)數(shù)據(jù)處理過程中進(jìn)行了落差系數(shù)修正、天平彈性角修正和兩心距修正。

    為了有效消除試驗(yàn)過程中設(shè)備振動(dòng)以及電磁干擾帶來的高頻噪聲,數(shù)據(jù)處理過程中采用了傅里葉變換數(shù)字濾波方法對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行消噪處理。將一個(gè)周期中氣動(dòng)力和力矩系數(shù)變換成試驗(yàn)振蕩頻率的傅里葉級(jí)數(shù),取6倍試驗(yàn)頻率為截止頻率,即前6項(xiàng)級(jí)數(shù)之和為最終試驗(yàn)結(jié)果。以偏航與滾轉(zhuǎn)耦合振蕩試驗(yàn)中的偏航力矩系數(shù)Cn為例,圖4給出了運(yùn)動(dòng)頻率f為0.2、0.4、0.6、0.8 Hz下偏航力矩的頻譜特性曲線。圖中:不同顏色圓圈代表不同頻率的6倍值所在位置??梢钥闯?,6倍運(yùn)動(dòng)頻率之后偏航力矩已經(jīng)衰減至很小,說明能量主要集中在模型運(yùn)動(dòng)頻率的前6階整數(shù)倍上。因此,將濾波截止頻率取為模型運(yùn)動(dòng)頻率的6倍是基本合理的[19-20]。

    圖4 偏航力矩頻譜分析Fig.4 Spectrum analysis of yawing moment

    1.4機(jī)構(gòu)角度與模型姿態(tài)角轉(zhuǎn)換關(guān)系

    1) 設(shè)轉(zhuǎn)盤角為θ1,規(guī)定從風(fēng)洞上方俯視時(shí)逆時(shí)針方向?yàn)檎?/p>

    2) 設(shè)液壓控制外框碳纖維彎刀擺角為φ1,規(guī)定下擺為正。

    3) 設(shè)伺服電機(jī)控制內(nèi)框滾轉(zhuǎn)支桿滾轉(zhuǎn)角為φ1,規(guī)定迎著來流方向,順時(shí)針為正。

    4) 模型迎角為α,側(cè)滑角為β。

    模型機(jī)體軸坐標(biāo)系定義見圖5。

    圖5 機(jī)體軸坐標(biāo)系示意圖Fig.5 Schematic diagram of body axis coordinate system

    以俯仰/滾轉(zhuǎn)耦合運(yùn)動(dòng)為例,如圖1所示,此時(shí)各機(jī)構(gòu)角為零,模型正裝時(shí),模型機(jī)體軸坐標(biāo)系與洞體軸坐標(biāo)系重合。通過將速度矢量在機(jī)體軸坐標(biāo)系上分解,可推導(dǎo)出俯仰/滾轉(zhuǎn)耦合運(yùn)動(dòng)的機(jī)構(gòu)角轉(zhuǎn)換關(guān)系為

    (1)

    sin β= -sin θ1cos φ1-cos θ1sin φ1sin φ1

    (2)

    令φ1=φ1+90°,代入式(1)、式(2)便可得到偏航/滾轉(zhuǎn)耦合運(yùn)動(dòng)的機(jī)構(gòu)角轉(zhuǎn)換關(guān)系為

    (3)

    sin β= -sin θ1sin φ1+cos θ1sin φ1cos φ1

    (4)

    驗(yàn)證試驗(yàn)過程中,各軸運(yùn)動(dòng)均采用正弦運(yùn)動(dòng)。

    2 驗(yàn)證試驗(yàn)方案

    采用動(dòng)態(tài)標(biāo)準(zhǔn)模型在FL-8風(fēng)洞對(duì)雙自由度大幅振蕩試驗(yàn)系統(tǒng)進(jìn)行驗(yàn)證試驗(yàn)。

    2.1試驗(yàn)?zāi)P?/p>

    試驗(yàn)?zāi)P蜑閯?dòng)態(tài)標(biāo)準(zhǔn)模型B4-05,模型結(jié)構(gòu)見圖6。該模型氣動(dòng)外形為四代機(jī)典型布局,采用翼身融合體、雙發(fā)雙垂尾,綜合優(yōu)化曲面外形,截尖菱形上單翼,V形傾斜雙垂尾,全動(dòng)平尾,具有較好的大迎角氣動(dòng)特性。模型全長為1.182 5 m,展長為0.847 5 m,機(jī)翼參考面積為0.304 688 m2,平均氣動(dòng)弦長為0.426 9 m。模型質(zhì)量為8 kg。

    圖6 動(dòng)態(tài)標(biāo)準(zhǔn)模型三維效果圖Fig.6 3D picture of standard dynamic model

    2.2試驗(yàn)內(nèi)容

    采用B4-05動(dòng)態(tài)標(biāo)準(zhǔn)模型在FL-8風(fēng)洞中進(jìn)行了雙自由度大幅振蕩重復(fù)性試驗(yàn)(見圖7),包括俯仰/滾轉(zhuǎn)耦合振蕩、偏航/滾轉(zhuǎn)耦合振蕩兩種運(yùn)動(dòng)的重復(fù)性試驗(yàn)以及振幅、頻率影響試驗(yàn)。試驗(yàn)條件為:試驗(yàn)風(fēng)速為25 m/s;迎角范圍為-10°~90°;側(cè)滑角范圍為-40°~40°;頻率范圍為0.2~0.8 Hz;振幅為20°和40°。

    圖7 動(dòng)態(tài)標(biāo)準(zhǔn)模型試驗(yàn)照片F(xiàn)ig.7 Photo of standard dynamic model test

    3 試驗(yàn)結(jié)果及分析

    3.1偏航/滾轉(zhuǎn)耦合振蕩試驗(yàn)

    試驗(yàn)中進(jìn)行了7次重復(fù)性試驗(yàn)(平衡角θ0=30°,振幅A=40°,頻率f=0.4 Hz,兩種運(yùn)動(dòng)的相位差Δψ為0° 和180°)。

    偏航/滾轉(zhuǎn)耦合運(yùn)動(dòng)的相位差定義為:相位差為0° 時(shí),偏航軸處于負(fù)的最大偏航角處,滾轉(zhuǎn)軸處于正的最大滾轉(zhuǎn)角處,即模型機(jī)頭向右指向運(yùn)動(dòng)的同時(shí)向左滾轉(zhuǎn);相位差為180° 時(shí),偏航軸處于負(fù)的最大偏航角處,滾轉(zhuǎn)軸處于負(fù)的最大滾轉(zhuǎn)角處,即模型機(jī)頭向右指向運(yùn)動(dòng)的同時(shí)向右滾轉(zhuǎn)。

    Yawingangle/(°)σClσCn-360.000500.00034-300.000460.00033-240.000320.00026-180.000170.00029-120.000140.00036-60.000170.0004100.000190.0004160.000210.00038120.000220.00034180.000220.00030240.000230.00025300.000280.00023360.000350.00031

    在偏航/滾轉(zhuǎn)耦合振蕩試驗(yàn)中,偏航、滾轉(zhuǎn)振蕩振幅均為40°,相位差有0° 和180° 兩種,平衡角θ0有10°、30°及50° 三種。由于篇幅所限,下文僅對(duì)平衡角θ0=30°、頻率f=0.4 Hz試驗(yàn)條件下的結(jié)果進(jìn)行分析。

    圖8為偏航/滾轉(zhuǎn)耦合振蕩Cl-φ、Cn-φ頻率影響。可見,在30° 平衡角下,從大滾轉(zhuǎn)角向小滾轉(zhuǎn)角運(yùn)動(dòng)過程中,遲滯現(xiàn)象不明顯,從小滾轉(zhuǎn)角向大滾轉(zhuǎn)角運(yùn)動(dòng)過程中,遲滯現(xiàn)象比較明顯,初步分析這種現(xiàn)象與運(yùn)動(dòng)過程中渦的再附過程有關(guān)。整個(gè)運(yùn)動(dòng)周期形成的耦合運(yùn)動(dòng)遲滯環(huán)面積隨著頻率的增大而增加,整體表現(xiàn)出的非定常效應(yīng)增強(qiáng)。

    圖8 偏航/滾轉(zhuǎn)耦合振蕩Cl- φ、Cn- φ頻率影響(Δψ =0°)Fig.8 Effect of frequency on Cl- φ, Cn- φ in yaw/rollcoupled oscillation test (Δψ=0°)

    圖9 偏航/滾轉(zhuǎn)耦合振蕩Cl- φ、Cn- φ相位影響Fig.9 Effect of phase on Cl- φ, Cn- φ in yaw/roll coupled oscillation test

    圖9為偏航/滾轉(zhuǎn)耦合振蕩Cl-φ、Cn-φ相位影響。可見,相位差對(duì)動(dòng)態(tài)遲滯環(huán)的形態(tài)和面積影響明顯。相比0° 相位差試驗(yàn)結(jié)果,180° 相位差的結(jié)果符號(hào)相反,量值增大,遲滯環(huán)方向相反。滾轉(zhuǎn)力矩元的遲滯環(huán)變得更飽滿,偏航力矩元的遲滯環(huán)則變窄。相位差給橫航向氣動(dòng)特性帶來的這些差異,主要是因?yàn)樵趦煞N運(yùn)動(dòng)過程中,模型側(cè)滑角不同導(dǎo)致的。

    3.2俯仰/滾轉(zhuǎn)耦合振蕩試驗(yàn)

    試驗(yàn)中進(jìn)行了7次重復(fù)性試驗(yàn)(平衡角θ0=40°,振幅A=40°,頻率f=0.4 Hz,兩種運(yùn)動(dòng)的相位差Δψ為0° 和180°)。

    俯仰/滾轉(zhuǎn)耦合運(yùn)動(dòng)相位差定義為:相位差為0° 時(shí),俯仰軸處于最低點(diǎn),滾轉(zhuǎn)軸處于正的最大滾轉(zhuǎn)角處,即模型上仰運(yùn)動(dòng)的同時(shí)向左滾轉(zhuǎn);相位差為180° 時(shí),俯仰軸處于最低點(diǎn),滾轉(zhuǎn)軸處于負(fù)的最大滾轉(zhuǎn)角處,即模型上仰運(yùn)動(dòng)的同時(shí)向右滾轉(zhuǎn)。

    Rollingangle/(°)σCNσCmσClσCn-360.004260.002090.000460.00066-300.005310.002400.000480.00049-240.005170.001600.000410.00020-180.004120.001170.000360.00051-120.003120.002060.000450.00079-60.002910.003040.000650.0009300.003560.003550.000840.0008960.004190.003440.000950.00071120.004300.002770.000940.00049180.003830.001800.000810.00039240.002890.000720.000530.00043300.002220.000410.000180.00039360.002290.000470.000130.00015

    在俯仰/滾轉(zhuǎn)耦合振蕩試驗(yàn)中,俯仰、滾轉(zhuǎn)振蕩振幅均為40°,相位差有0° 和180° 兩種。下文給出頻率及相位差的影響,并進(jìn)行簡要分析。

    圖10為俯仰/滾轉(zhuǎn)耦合振蕩CN-θ、Cm-θ、Cl-φ和Cn-φ頻率影響。圖中:θ為俯仰角??梢钥闯觯瑢?duì)于法向力元,上仰過程中,32° 迎角后,非定常氣動(dòng)力增量受頻率影響明顯增大,而在下俯過程中,60° 迎角以下,非定常氣動(dòng)力增量受頻率影響才顯現(xiàn)。顯然在上仰過程中,流場結(jié)構(gòu)從附著渦流向完全渦破裂狀態(tài)過渡時(shí),頻率的變化強(qiáng)烈地影響到氣動(dòng)力的非線性變化,而下俯過程中,同樣是渦破裂流動(dòng)向附著渦流過渡時(shí)氣動(dòng)力非線性變化更顯著,這一現(xiàn)象本質(zhì)上就是源于流動(dòng)結(jié)構(gòu)的轉(zhuǎn)換效應(yīng)。

    對(duì)于俯仰力矩元,隨著頻率的增大,遲滯環(huán)面積增大,但“8”字環(huán)幾乎在同一點(diǎn)相交,即阻尼特性區(qū)域沒有隨頻率發(fā)生變化。

    對(duì)于滾轉(zhuǎn)力矩元和偏航力矩元,由于其遲滯環(huán)面積受迎角變化顯著,大迎角時(shí),模型處于負(fù)滾轉(zhuǎn)位置,所以對(duì)于正負(fù)滾轉(zhuǎn)位置,其力矩曲線嚴(yán)重不對(duì)稱,但其環(huán)的面積仍符合隨頻率增加而增大的規(guī)律。

    圖10 俯仰/滾轉(zhuǎn)耦合振蕩CN- θ、Cm- θ、Cl- φ和Cn- φ頻率影響(Δψ =0°)Fig.10 Effect of frequency on CN- θ, Cm- θ, Cl- φ and Cn- φ in pitch/roll coupled oscillation test (Δψ =0°)

    圖11為俯仰/滾轉(zhuǎn)耦合振蕩CN-θ、Cm-θ、Cl-φ和Cn-φ相位影響。可以看出,同一頻率下,不同相位差的俯仰/滾轉(zhuǎn)耦合振蕩,其縱向氣動(dòng)力特性變化不明顯,主要是因?yàn)橄鄬?duì)于機(jī)身對(duì)稱面來說,兩種耦合運(yùn)動(dòng)是對(duì)稱的,繞機(jī)身流動(dòng)對(duì)飛機(jī)的縱向氣動(dòng)特性影響基本是相同的;由于側(cè)滑角相反,不同相位差的俯仰/滾轉(zhuǎn)耦合振蕩橫航向氣動(dòng)力矩系數(shù)剛好相反,但阻尼發(fā)散特性不變[5]。

    圖11 俯仰/滾轉(zhuǎn)耦合振蕩CN- θ、Cm- θ、Cl- φ和Cn- φ相位影響Fig.11 Effect of phase on CN- θ, Cm- θ, Cl- φ and Cn- φ in pitch/roll coupled oscillation test

    4 結(jié) 論

    1) 中國航空工業(yè)空氣動(dòng)力研究院研制的低速風(fēng)洞雙自由度大幅振蕩試驗(yàn)系統(tǒng)采用電機(jī)和液壓伺服控制方式與實(shí)時(shí)操作系統(tǒng),保證了系統(tǒng)控制與數(shù)據(jù)采集的精度;采用液壓驅(qū)動(dòng)方式,保證了系統(tǒng)運(yùn)動(dòng)穩(wěn)定,結(jié)構(gòu)間隙大幅減小,試驗(yàn)重復(fù)性精度較高。

    2) 驗(yàn)證試驗(yàn)結(jié)果顯示,模型進(jìn)行耦合振蕩時(shí),由于模型繞流流態(tài)變化復(fù)雜,氣動(dòng)力呈現(xiàn)明顯的非定常特性與非線性。隨振蕩頻率增大,各元?dú)鈩?dòng)力系數(shù)或氣動(dòng)力矩系數(shù)的遲滯環(huán)面積增加;不同元“8”字環(huán)的交點(diǎn)位置變化特性也有所不同。兩軸運(yùn)動(dòng)相位差改變180°對(duì)耦合運(yùn)動(dòng)試驗(yàn)結(jié)果影響明顯,對(duì)于偏航/滾轉(zhuǎn)耦合振蕩來說,遲滯環(huán)的形狀、面積和方向均有較大變化;對(duì)于俯仰/滾轉(zhuǎn)耦合振蕩來說,橫航向力矩系數(shù)曲線關(guān)于原點(diǎn)對(duì)稱,法向力系數(shù)和俯仰力矩系數(shù)沒有變化。

    5 展 望

    雙自由度大幅振蕩試驗(yàn)技術(shù)為獲取飛行器非定常氣動(dòng)力特性提供了試驗(yàn)手段,在此基礎(chǔ)上,可將模型表面測(cè)壓技術(shù)與之結(jié)合起來,用以獲取雙自由度大幅運(yùn)動(dòng)過程中模型表面非定常壓力分布,為開展飛行器表面復(fù)雜流動(dòng)機(jī)理研究和飛行器結(jié)構(gòu)設(shè)計(jì)提供依據(jù)。

    [1]BRIDGES D H. The asymmetric vortex wake problem-asking the right question: AIAA-2006-3553[R]. Reston: AIAA, 2006.

    [2]WILLIAM J, GILLARD P E. AFRL F-22 dynamic wind tunnel test results: AIAA-1999-4015[R]. Reston: AIAA, 1999.

    [3]BRANDON J M, SHAN 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.

    [4]MICHAEL C, FRANCIS S. X-31 enhanced fighter maneuver ability demonstrator: Flight test achievements: ICAS-94-7.2.2[R]. [S.l.:s.n.], 1994.

    [5]李周復(fù). 風(fēng)洞特種試驗(yàn)技術(shù)[M]. 北京: 航空工業(yè)出版社, 2010: 257-259.

    LI Z F. Special wind tunnel test[M]. Beijing: Aviation Industry Press, 2010: 257-259 (in Chinese).

    [6]ERICSSON L E. Sources of high alpha vortex asymmetry at zero sideslip[J]. Journal of Aircraft, 1992, 29(6): 1086-1090.

    [7]楊勐, 黃達(dá). 大迎角耦合運(yùn)動(dòng)非定常空氣動(dòng)力特性[J]. 實(shí)驗(yàn)流體力學(xué), 2011, 25(6): 19-22.

    YANG M, HUANG D. The unsteady aerodynamic characteristics of coupled motion at high angle of attack[J]. Journal of Experiments in Fluid Mechanics, 2011, 25(6): 19-22 (in Chinese).

    [8]黃達(dá), 談佳楨. 飛行器大振幅運(yùn)動(dòng)實(shí)驗(yàn)與氣動(dòng)力建模[J]. 南京航空航天大學(xué)學(xué)報(bào), 2010, 42(4): 407-410.

    HUANG D, TAN J Z. Experiment about aircraft motion in large amplitude and aerodynamic modeling[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2010, 42(4): 407-410 (in Chinese).

    [9]黃達(dá). 飛行器大幅振蕩運(yùn)動(dòng)非定??諝鈩?dòng)力特性研究[D]. 南京: 南京航空航天大學(xué), 2007: 1-3.

    HUANG D. Unsteady aerodynamic characteristics for the aircraft oscilation in large amplitude[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2007: 1-3 (in Chinese).

    [10]史志偉, 吳根興, 黃達(dá). 基于大振幅諧波運(yùn)動(dòng)的非定常氣動(dòng)模型風(fēng)洞實(shí)驗(yàn)驗(yàn)證[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2010, 28(6): 650-654.

    SHI Z W, WU G X, HUANG D. The validation of unsteady aerodynamic model based on the high amplitude harmonic rolling motion test[J]. Acta Aerodynamica Sinica, 2010, 28(6): 650-654 (in Chinese).

    [11]黃達(dá), 李志強(qiáng), 史志偉, 等. 飛機(jī)大振幅非定常滾轉(zhuǎn)運(yùn)動(dòng)的非線性穩(wěn)定性分析[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2000, 18(4): 401-406.

    HUANG D, LI Z Q, SHI Z W, et al. A nonlinear stability analysis for an aircraft rolling in large amplitude motion[J]. Acta Aerodynamica Sinica, 2000, 18(4): 401-406 (in Chinese).

    [12]高劍軍, 卜忱, 杜希奇. 旋轉(zhuǎn)流場下模型大幅俯仰運(yùn)動(dòng)實(shí)驗(yàn)研究[J]. 實(shí)驗(yàn)力學(xué), 2010, 25(5): 207-211.

    GAO J J, BU C, DU X Q. Experimental investigation on aerodynamic characteristics of model in large amplitude pitching movement in a rotational flow field[J]. Journal of Experimental Mechanics, 2010, 25(5): 207-211 (in Chinese).

    [13]楊勐, 黃達(dá), 吳根興. 基于風(fēng)洞試驗(yàn)的非定常空氣動(dòng)力模型識(shí)別[J]. 南京航空航天大學(xué)學(xué)報(bào), 2011, 43(2): 149-153.

    YANG M, HUANG D, WU G X. Identification of unsteady aerodynamic model based on wind tunnel test[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2011, 43(2): 149-153 (in Chinese).

    [14]黃達(dá), 吳根興. 三角翼俯仰滾轉(zhuǎn)耦合運(yùn)動(dòng)氣動(dòng)特性研究[J]. 航空學(xué)報(bào), 1999, 20(6): 485-488.

    HUANG D, WU G X. Investigation of unsteady aerodynamic characteristics for a delta wing oscillating in large amplitude pitching roll motion[J]. Acta Aeronautica et Astronautica Sinica, 1999, 20(6): 485-488 (in Chinese).

    [15]黃達(dá), 吳根興. 飛機(jī)偏航-滾轉(zhuǎn)耦合運(yùn)動(dòng)非定??諝鈩?dòng)力實(shí)驗(yàn)[J]. 南京航空航天大學(xué)學(xué)報(bào), 2005, 37(4): 408-411.

    HUANG D, WU G X. Experiment on fighter oscillating in large amplitude yaw-roll motion[J]. Journal of Nanjing University of Aeronautics and Astronautics, 2005, 37(4): 408-411 (in Chinese).

    [16]黃達(dá), 吳根興. 飛機(jī)大振幅滾轉(zhuǎn)運(yùn)動(dòng)動(dòng)態(tài)氣動(dòng)特性實(shí)驗(yàn)研究[J]. 空氣動(dòng)力學(xué)學(xué)報(bào), 2005, 23(2): 135-138.

    HUANG D, WU G X. Investigation of the unsteady aerodynamic characteristics for a fighter rolling in large amplitude[J]. Acta Aerodynamica Sinica, 2005, 23(2): 135-138 (in Chinese).

    [17]李志強(qiáng), 黃達(dá), 史志偉, 等. 俯仰-滾轉(zhuǎn)耦合兩自由度大振幅非定常實(shí)驗(yàn)技術(shù)[J]. 南京航空航天大學(xué)學(xué)報(bào), 1999, 31(2): 121-126.

    LI Z Q, HUANG D, SHI Z W, et al. Test technology of unsteady aerodynamic characteristic for a model oscillating in large amplitude pitching rolling motion[J]. Journal of Nanjing University of Aeronautics and Astronautics, 1999, 31(2): 121-126 (in Chinese).

    [18]謝志江, 孫小勇, 孫海生, 等. 低速風(fēng)洞動(dòng)態(tài)試驗(yàn)的高速并聯(lián)機(jī)構(gòu)設(shè)計(jì)及動(dòng)力學(xué)分析[J]. 航空學(xué)報(bào), 2013, 34(3): 487-494.

    XIE Z J, SUN X Y, SUN H S, et al. Mechanism design and dynamics analysis of high speed parallel robot for dynamic test in low speed wind tunnel[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(3): 487-494 (in Chinese).

    [19]蔣永, 孫海生, 沈志宏, 等. 非定常大振幅振蕩試驗(yàn)數(shù)據(jù)處理研究[J]. 實(shí)驗(yàn)流體力學(xué), 2015, 29(1): 97-102.

    JIANG Y, SUN H S, SHEN Z H, et al. Research on the data processing of unsteady large-amplitude-oscilltion experiment[J]. Jourmal of Experiments in Fluid Mechanics, 2015, 29(1): 97-102 (in Chinese).

    [20]徐永長, 鄭世華. 動(dòng)態(tài)失速風(fēng)洞實(shí)驗(yàn)數(shù)據(jù)處理中頻譜分析與數(shù)字濾波[J]. 流體力學(xué)實(shí)驗(yàn)與測(cè)量, 2000, 14(3): 79-82.

    XU Y C, ZHENG S H. Spectrum analysis and digital filtering to experimental data of dynamic stall test in wind tunnel[J]. Experiments and Measurements in Fluid Mechanics, 2000, 14(3): 79-82 (in Chinese).

    劉春明男, 碩士, 高級(jí)工程師。主要研究方向: 風(fēng)洞動(dòng)態(tài)試驗(yàn)技術(shù)。

    Tel.: 0451-87570255

    E-mail: lcmqdy@163.com

    趙志軍男, 碩士, 高級(jí)工程師。主要研究方向: 空氣動(dòng)力學(xué)試驗(yàn)。

    Tel.: 010-62654599

    E-mail: 1943223664@qq.com

    Double degree-of-freedom large amplitude oscillation test technology in low speed wind tunnel

    LIU Chunming1, *, ZHAO Zhijun2, BU Chen1, WANG Jianfeng1, MU Weiqiang1

    1. AVIC Aerodynamics Research Institute, Harbin150001, China 2. China Aerodynamics Research and Development Center, Mianyang621000, China

    AVIC Aerodynamics Research Institute developed a double degree-of-freedom large amplitude oscillation testing technique based on driving technology coupled hydraulic motor and servo motor in FL-8 wind tunnel in order to study the aircraft unsteady aerodynamics during post maneuver at high angle-of-attack. This paper presents the basic components and test principle, as well as the FL-8 wind tunnel test results of the standard model test. It is indicated that the system has high stability, reliability and repeatability precision. The technology can be effectively applied to studying and predicting the aerodynamic characteristics of double degree-of-freedom coupled motion.

    high angle-of-attack; large amplitude oscillation; double degree-of-freedom; unsteady aerodynamics; wind tunnel test

    2016-01-25; Revised: 2016-02-15; Accepted: 2016-05-23; Published online: 2016-06-1214:16

    Aeronautical Science Foundation of China (20143203002)

    . Tel.: 0451-87570255E-mail: lcmqdy@163.com

    2016-01-25; 退修日期: 2016-02-15; 錄用日期: 2016-05-23;

    時(shí)間: 2016-06-1214:16

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

    航空科學(xué)基金 (20143203002)

    .Tel.: 0451-87570255E-mail: lcmqdy@163.com

    10.7527/S1000-6893.2016.0154

    V212.1

    A

    1000-6893(2016)08-2417-09

    引用格式: 劉春明, 趙志軍, 卜忱, 等. 低速風(fēng)洞雙自由度大幅振蕩試驗(yàn)技術(shù)[J]. 航空學(xué)報(bào), 2016, 37(8): 2417-2425. LIU C M, ZHAO Z J, BU C, et al. Double degree-of-freedom large amplitude oscillation test technology in low speed wind tunnel[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2417-2425.

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

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

    猜你喜歡
    氣動(dòng)力迎角風(fēng)洞
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    飛行載荷外部氣動(dòng)力的二次規(guī)劃等效映射方法
    斑頭雁進(jìn)風(fēng)洞
    黃風(fēng)洞貂鼠精
    基于NI cRIO平臺(tái)的脈沖燃燒風(fēng)洞控制系統(tǒng)設(shè)計(jì)
    側(cè)風(fēng)對(duì)拍動(dòng)翅氣動(dòng)力的影響
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    風(fēng)力機(jī)氣動(dòng)力不對(duì)稱故障建模與仿真
    飛行器風(fēng)洞模型的快速制造技術(shù)
    日本一区二区免费在线视频| 国产亚洲午夜精品一区二区久久| 日韩人妻精品一区2区三区| 亚洲,欧美精品.| 无遮挡黄片免费观看| 精品亚洲成a人片在线观看| 天堂俺去俺来也www色官网| 欧美人与性动交α欧美软件| 久久精品亚洲熟妇少妇任你| 亚洲精品国产av成人精品| 精品一区在线观看国产| 日本猛色少妇xxxxx猛交久久| 视频区图区小说| 波多野结衣一区麻豆| 久久性视频一级片| 色精品久久人妻99蜜桃| 人成视频在线观看免费观看| av线在线观看网站| 亚洲欧美日韩另类电影网站| a级毛片黄视频| 男人爽女人下面视频在线观看| 亚洲欧美一区二区三区久久| 热re99久久国产66热| 亚洲欧洲日产国产| 久久久精品94久久精品| 亚洲欧洲精品一区二区精品久久久 | 国产1区2区3区精品| 欧美日韩视频高清一区二区三区二| 一级毛片电影观看| 久久午夜综合久久蜜桃| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费日韩欧美大片| a级毛片在线看网站| 国产免费福利视频在线观看| 叶爱在线成人免费视频播放| 大陆偷拍与自拍| 国产av码专区亚洲av| 久久久国产精品麻豆| 亚洲成av片中文字幕在线观看| 欧美黄色片欧美黄色片| 亚洲国产精品999| 国产在线视频一区二区| 亚洲国产日韩一区二区| 人妻一区二区av| 午夜免费男女啪啪视频观看| 精品一品国产午夜福利视频| 欧美97在线视频| 亚洲国产毛片av蜜桃av| 午夜福利视频在线观看免费| 亚洲少妇的诱惑av| 涩涩av久久男人的天堂| 少妇精品久久久久久久| 天天影视国产精品| 人妻 亚洲 视频| 久久婷婷青草| 日韩欧美一区视频在线观看| 欧美97在线视频| 巨乳人妻的诱惑在线观看| 国产又爽黄色视频| 赤兔流量卡办理| 日韩精品免费视频一区二区三区| 下体分泌物呈黄色| 美女国产高潮福利片在线看| av不卡在线播放| 777米奇影视久久| 一区二区av电影网| 日韩av不卡免费在线播放| 大话2 男鬼变身卡| 制服丝袜香蕉在线| 美女大奶头黄色视频| 大话2 男鬼变身卡| 天天躁夜夜躁狠狠久久av| 免费高清在线观看视频在线观看| 亚洲av中文av极速乱| 亚洲精品日韩在线中文字幕| 亚洲国产欧美一区二区综合| 欧美在线一区亚洲| 国产黄频视频在线观看| 亚洲美女黄色视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| bbb黄色大片| 亚洲国产av影院在线观看| 大片电影免费在线观看免费| 亚洲欧洲精品一区二区精品久久久 | 欧美日韩福利视频一区二区| 日日撸夜夜添| 18禁裸乳无遮挡动漫免费视频| 少妇的丰满在线观看| 国产精品亚洲av一区麻豆 | 亚洲图色成人| 亚洲成人免费av在线播放| 免费黄频网站在线观看国产| 久久狼人影院| 午夜免费男女啪啪视频观看| 国产成人一区二区在线| 美女高潮到喷水免费观看| 欧美最新免费一区二区三区| 满18在线观看网站| 亚洲精品久久久久久婷婷小说| 如日韩欧美国产精品一区二区三区| 99热国产这里只有精品6| 视频区图区小说| 日韩一卡2卡3卡4卡2021年| 大码成人一级视频| av一本久久久久| 综合色丁香网| 久久韩国三级中文字幕| 亚洲国产最新在线播放| 国产精品一区二区精品视频观看| 丝袜在线中文字幕| 深夜精品福利| 赤兔流量卡办理| 欧美日本中文国产一区发布| 国产伦人伦偷精品视频| 国产女主播在线喷水免费视频网站| 这个男人来自地球电影免费观看 | a 毛片基地| 激情视频va一区二区三区| 九九爱精品视频在线观看| 丁香六月天网| 赤兔流量卡办理| 岛国毛片在线播放| 中文字幕人妻丝袜制服| 亚洲精品日韩在线中文字幕| 操出白浆在线播放| 国产亚洲欧美精品永久| 日韩中文字幕视频在线看片| svipshipincom国产片| 色视频在线一区二区三区| 超碰成人久久| 超碰97精品在线观看| 十八禁高潮呻吟视频| 男女无遮挡免费网站观看| 欧美人与善性xxx| 欧美日韩一区二区视频在线观看视频在线| 1024视频免费在线观看| 久久久精品国产亚洲av高清涩受| 国产亚洲欧美精品永久| 青青草视频在线视频观看| 亚洲国产最新在线播放| 黄网站色视频无遮挡免费观看| 熟妇人妻不卡中文字幕| 国产又爽黄色视频| 黄色视频不卡| 51午夜福利影视在线观看| 成人手机av| 国产精品人妻久久久影院| 青春草亚洲视频在线观看| 夜夜骑夜夜射夜夜干| 一级黄片播放器| www.自偷自拍.com| 欧美黑人精品巨大| 国产成人精品福利久久| 亚洲在久久综合| 秋霞伦理黄片| 伦理电影大哥的女人| 日韩av在线免费看完整版不卡| 欧美亚洲 丝袜 人妻 在线| 亚洲美女视频黄频| 午夜久久久在线观看| 亚洲色图综合在线观看| 五月天丁香电影| 啦啦啦啦在线视频资源| 男女午夜视频在线观看| 久久精品aⅴ一区二区三区四区| 黄片小视频在线播放| 亚洲av日韩精品久久久久久密 | 精品一区二区三卡| 亚洲第一青青草原| 你懂的网址亚洲精品在线观看| 汤姆久久久久久久影院中文字幕| 久久久精品国产亚洲av高清涩受| 男女国产视频网站| 嫩草影院入口| 侵犯人妻中文字幕一二三四区| 在线观看人妻少妇| a级片在线免费高清观看视频| 欧美日韩一区二区视频在线观看视频在线| 免费高清在线观看视频在线观看| 免费不卡黄色视频| 精品第一国产精品| 亚洲自偷自拍图片 自拍| 热re99久久精品国产66热6| 午夜福利视频在线观看免费| 亚洲精品美女久久久久99蜜臀 | 国产av精品麻豆| 一区二区三区精品91| 777米奇影视久久| 在线观看免费高清a一片| 又黄又粗又硬又大视频| 天堂俺去俺来也www色官网| 国产又爽黄色视频| 亚洲成av片中文字幕在线观看| 大话2 男鬼变身卡| 免费在线观看黄色视频的| 别揉我奶头~嗯~啊~动态视频 | 亚洲图色成人| 亚洲精品美女久久久久99蜜臀 | 一本色道久久久久久精品综合| 亚洲国产av新网站| 男女国产视频网站| 一本大道久久a久久精品| 欧美日韩亚洲高清精品| 免费少妇av软件| av有码第一页| 亚洲av国产av综合av卡| 成年女人毛片免费观看观看9 | 精品国产超薄肉色丝袜足j| 青春草视频在线免费观看| 国产一区二区激情短视频 | 精品亚洲乱码少妇综合久久| 人妻一区二区av| 久久天躁狠狠躁夜夜2o2o | 777米奇影视久久| 91精品伊人久久大香线蕉| 国产欧美日韩一区二区三区在线| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 亚洲国产av新网站| 久久久久久久精品精品| 一二三四在线观看免费中文在| 在线观看一区二区三区激情| 亚洲在久久综合| 超色免费av| 精品午夜福利在线看| 色播在线永久视频| 婷婷成人精品国产| 考比视频在线观看| 国产伦人伦偷精品视频| 香蕉国产在线看| av一本久久久久| 日本欧美国产在线视频| 男人舔女人的私密视频| 亚洲国产精品999| 人成视频在线观看免费观看| 国产成人免费观看mmmm| 最近中文字幕高清免费大全6| 国产日韩欧美视频二区| 国产xxxxx性猛交| 亚洲精品美女久久av网站| 免费高清在线观看视频在线观看| 熟女少妇亚洲综合色aaa.| 99久久综合免费| 日韩一区二区三区影片| 男女无遮挡免费网站观看| 黑丝袜美女国产一区| 99久久精品国产亚洲精品| 国产高清国产精品国产三级| 欧美国产精品va在线观看不卡| 亚洲少妇的诱惑av| 看免费成人av毛片| 中文欧美无线码| 一本色道久久久久久精品综合| 王馨瑶露胸无遮挡在线观看| 欧美日韩亚洲高清精品| 亚洲,欧美,日韩| 国产一区二区三区av在线| 国产亚洲最大av| 又黄又粗又硬又大视频| 十八禁高潮呻吟视频| 亚洲国产日韩一区二区| 久久天躁狠狠躁夜夜2o2o | 精品国产一区二区久久| 人人妻人人添人人爽欧美一区卜| 精品少妇一区二区三区视频日本电影 | 国产黄色视频一区二区在线观看| 日本91视频免费播放| 五月开心婷婷网| 日本色播在线视频| 亚洲精品国产色婷婷电影| 久久国产亚洲av麻豆专区| 纯流量卡能插随身wifi吗| 熟女av电影| 国产黄色免费在线视频| 精品一区在线观看国产| 狂野欧美激情性bbbbbb| 日韩一区二区视频免费看| 满18在线观看网站| 国产一区二区在线观看av| 精品少妇内射三级| 久久人人97超碰香蕉20202| 亚洲国产日韩一区二区| 99精品久久久久人妻精品| 黄色怎么调成土黄色| 晚上一个人看的免费电影| 宅男免费午夜| 欧美在线一区亚洲| 久久久久视频综合| 欧美日韩国产mv在线观看视频| 亚洲精品国产av蜜桃| 国产精品.久久久| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕高清在线视频| www.自偷自拍.com| 99精国产麻豆久久婷婷| 高清av免费在线| 在线观看www视频免费| kizo精华| 汤姆久久久久久久影院中文字幕| 国产精品 欧美亚洲| 亚洲欧美精品综合一区二区三区| 精品视频人人做人人爽| 美女大奶头黄色视频| 亚洲精品国产一区二区精华液| 日韩 亚洲 欧美在线| 成年动漫av网址| 麻豆av在线久日| 制服人妻中文乱码| 一级a爱视频在线免费观看| kizo精华| 美女主播在线视频| 在线观看免费高清a一片| 亚洲国产欧美日韩在线播放| 啦啦啦 在线观看视频| 亚洲欧洲国产日韩| 亚洲国产最新在线播放| 国产精品久久久av美女十八| 51午夜福利影视在线观看| 一级片'在线观看视频| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 国产精品 国内视频| 男女高潮啪啪啪动态图| 最近2019中文字幕mv第一页| 亚洲视频免费观看视频| 欧美亚洲 丝袜 人妻 在线| 国产精品一区二区在线不卡| 激情视频va一区二区三区| 国产一区二区三区综合在线观看| 国产一区二区三区综合在线观看| 两性夫妻黄色片| 亚洲精品av麻豆狂野| 国产在线一区二区三区精| 免费女性裸体啪啪无遮挡网站| 一边亲一边摸免费视频| 狂野欧美激情性bbbbbb| 欧美人与善性xxx| 久久久久精品久久久久真实原创| 卡戴珊不雅视频在线播放| 男人舔女人的私密视频| 成人亚洲精品一区在线观看| 亚洲av男天堂| 国产片特级美女逼逼视频| 久久毛片免费看一区二区三区| 久久精品亚洲av国产电影网| 丁香六月欧美| 午夜福利视频在线观看免费| 天天操日日干夜夜撸| √禁漫天堂资源中文www| 女人爽到高潮嗷嗷叫在线视频| 国产精品香港三级国产av潘金莲 | 久久99一区二区三区| 国产日韩一区二区三区精品不卡| 国产成人欧美| 久久精品人人爽人人爽视色| 少妇人妻 视频| 亚洲精品av麻豆狂野| 99久久人妻综合| 免费观看av网站的网址| 黑人猛操日本美女一级片| 国产亚洲一区二区精品| 啦啦啦 在线观看视频| 日日撸夜夜添| 狂野欧美激情性bbbbbb| 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻熟女乱码| 欧美日韩国产mv在线观看视频| 亚洲第一av免费看| 老鸭窝网址在线观看| 9色porny在线观看| 亚洲美女搞黄在线观看| 精品免费久久久久久久清纯 | 在线观看免费视频网站a站| 免费高清在线观看视频在线观看| 啦啦啦视频在线资源免费观看| 国产亚洲精品第一综合不卡| 色精品久久人妻99蜜桃| 人人妻人人添人人爽欧美一区卜| 美女视频免费永久观看网站| 香蕉丝袜av| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 19禁男女啪啪无遮挡网站| 亚洲av男天堂| 久久久久国产精品人妻一区二区| 宅男免费午夜| bbb黄色大片| 2021少妇久久久久久久久久久| 国产精品香港三级国产av潘金莲 | 汤姆久久久久久久影院中文字幕| av国产精品久久久久影院| tube8黄色片| 在线观看免费日韩欧美大片| 精品卡一卡二卡四卡免费| 午夜福利一区二区在线看| 99香蕉大伊视频| 丝袜脚勾引网站| 久久久久精品人妻al黑| 最新在线观看一区二区三区 | 青春草亚洲视频在线观看| www.自偷自拍.com| 成人毛片60女人毛片免费| 久久久久久久精品精品| 无遮挡黄片免费观看| 亚洲av中文av极速乱| 国产99久久九九免费精品| 国产日韩欧美在线精品| 国产黄频视频在线观看| 久久精品国产综合久久久| 日韩欧美一区视频在线观看| 日韩精品有码人妻一区| 国产片特级美女逼逼视频| 在线精品无人区一区二区三| 18禁裸乳无遮挡动漫免费视频| 精品少妇一区二区三区视频日本电影 | 国产精品国产av在线观看| 少妇被粗大猛烈的视频| 丰满少妇做爰视频| 午夜福利一区二区在线看| 操出白浆在线播放| 下体分泌物呈黄色| 国产高清不卡午夜福利| 超色免费av| 如日韩欧美国产精品一区二区三区| 国产男女内射视频| 亚洲国产欧美日韩在线播放| 满18在线观看网站| 久久国产亚洲av麻豆专区| 亚洲国产看品久久| 久久精品久久久久久噜噜老黄| 一边摸一边抽搐一进一出视频| 国产不卡av网站在线观看| 十分钟在线观看高清视频www| 超碰97精品在线观看| 精品国产一区二区三区久久久樱花| 欧美精品av麻豆av| 国产一区有黄有色的免费视频| 午夜av观看不卡| 美国免费a级毛片| 人人妻人人爽人人添夜夜欢视频| xxx大片免费视频| 丝袜在线中文字幕| 中文天堂在线官网| 精品少妇内射三级| 久久av网站| 国产亚洲欧美精品永久| 久久久久久人妻| 这个男人来自地球电影免费观看 | 青春草亚洲视频在线观看| 精品一区在线观看国产| 欧美成人精品欧美一级黄| 天堂中文最新版在线下载| 精品一区二区三区av网在线观看 | 婷婷色av中文字幕| 不卡视频在线观看欧美| 欧美中文综合在线视频| 国产野战对白在线观看| 日本午夜av视频| 最近的中文字幕免费完整| 亚洲第一区二区三区不卡| 91精品国产国语对白视频| 80岁老熟妇乱子伦牲交| 亚洲成国产人片在线观看| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 国产精品熟女久久久久浪| 黄色一级大片看看| 久久久久精品国产欧美久久久 | 日韩av不卡免费在线播放| 午夜福利视频精品| 男女床上黄色一级片免费看| 在线观看三级黄色| 飞空精品影院首页| 免费人妻精品一区二区三区视频| 老司机影院成人| 性高湖久久久久久久久免费观看| 国产一区二区 视频在线| 日韩一区二区视频免费看| 成人国产av品久久久| 啦啦啦在线免费观看视频4| 亚洲欧美成人精品一区二区| 亚洲国产欧美在线一区| 男人舔女人的私密视频| 一区二区三区精品91| 九色亚洲精品在线播放| 久久人妻熟女aⅴ| 国语对白做爰xxxⅹ性视频网站| 黄色一级大片看看| 国产精品一区二区在线不卡| 欧美在线黄色| 操美女的视频在线观看| 国产精品蜜桃在线观看| 一区二区三区精品91| 欧美精品亚洲一区二区| 成人毛片60女人毛片免费| 王馨瑶露胸无遮挡在线观看| 岛国毛片在线播放| 综合色丁香网| 国产伦理片在线播放av一区| 免费av中文字幕在线| 美女中出高潮动态图| 免费不卡黄色视频| 国产无遮挡羞羞视频在线观看| 视频在线观看一区二区三区| 美女国产高潮福利片在线看| 人体艺术视频欧美日本| 久久国产亚洲av麻豆专区| 亚洲成人国产一区在线观看 | 2018国产大陆天天弄谢| 男女午夜视频在线观看| 日本一区二区免费在线视频| 老汉色∧v一级毛片| 亚洲一区二区三区欧美精品| 精品人妻熟女毛片av久久网站| 亚洲免费av在线视频| 国产精品蜜桃在线观看| 久久精品亚洲av国产电影网| 婷婷色综合www| 久久国产精品大桥未久av| 最近中文字幕高清免费大全6| 免费黄网站久久成人精品| 国产av码专区亚洲av| 国产精品一区二区在线不卡| h视频一区二区三区| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| 成人免费观看视频高清| 精品一区在线观看国产| 国产精品香港三级国产av潘金莲 | 伦理电影大哥的女人| 亚洲精品av麻豆狂野| 老汉色av国产亚洲站长工具| 亚洲欧洲精品一区二区精品久久久 | 悠悠久久av| 超色免费av| 久久久精品区二区三区| 这个男人来自地球电影免费观看 | 人成视频在线观看免费观看| www日本在线高清视频| 久久人妻熟女aⅴ| 亚洲自偷自拍图片 自拍| 国产精品av久久久久免费| 咕卡用的链子| 五月开心婷婷网| 中国国产av一级| 久久天堂一区二区三区四区| 国产精品免费大片| 国产xxxxx性猛交| 色综合欧美亚洲国产小说| 精品国产一区二区三区四区第35| 狂野欧美激情性xxxx| 国产乱人偷精品视频| 国产国语露脸激情在线看| 自线自在国产av| 久久这里只有精品19| 久久精品国产a三级三级三级| 亚洲欧美成人综合另类久久久| 久久久久久人人人人人| 亚洲一卡2卡3卡4卡5卡精品中文| 麻豆精品久久久久久蜜桃| 汤姆久久久久久久影院中文字幕| 韩国av在线不卡| 十八禁人妻一区二区| 青春草国产在线视频| 日韩人妻精品一区2区三区| 国产 精品1| 亚洲精品自拍成人| 中文乱码字字幕精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 天天躁日日躁夜夜躁夜夜| 午夜福利在线免费观看网站| 午夜91福利影院| 成人亚洲欧美一区二区av| 美国免费a级毛片| 成人亚洲欧美一区二区av| av电影中文网址| 国产熟女欧美一区二区| 午夜福利影视在线免费观看| 精品国产露脸久久av麻豆| 国产免费视频播放在线视频| 九草在线视频观看| 午夜免费鲁丝| 黄色视频不卡| 久久韩国三级中文字幕| 国产男女超爽视频在线观看| 毛片一级片免费看久久久久| 一本一本久久a久久精品综合妖精| 国产成人av激情在线播放| 日本黄色日本黄色录像| av国产精品久久久久影院| 哪个播放器可以免费观看大片| 婷婷成人精品国产| av福利片在线| 99久久人妻综合| 一本—道久久a久久精品蜜桃钙片| 亚洲少妇的诱惑av| 色综合欧美亚洲国产小说| 精品人妻在线不人妻| 国产高清国产精品国产三级| 国产精品二区激情视频| 最近中文字幕高清免费大全6| 岛国毛片在线播放| 久久精品国产综合久久久| 欧美日韩一区二区视频在线观看视频在线| 日本欧美国产在线视频| 日韩av不卡免费在线播放| 麻豆精品久久久久久蜜桃| 国产熟女午夜一区二区三区| 天天躁夜夜躁狠狠躁躁| 亚洲免费av在线视频| 亚洲第一av免费看| 中文精品一卡2卡3卡4更新| 成人午夜精彩视频在线观看| 两个人看的免费小视频| 在线天堂中文资源库|