• <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ù)
    一区二区三区精品91| 欧美日韩一级在线毛片| 欧美激情 高清一区二区三区| 丝袜美足系列| 麻豆国产av国片精品| 搡老岳熟女国产| 两性午夜刺激爽爽歪歪视频在线观看 | 国产片内射在线| a级片在线免费高清观看视频| 窝窝影院91人妻| 美女大奶头黄色视频| 女人被躁到高潮嗷嗷叫费观| 国产日韩一区二区三区精品不卡| 国产男人的电影天堂91| 深夜精品福利| 人妻人人澡人人爽人人| 免费观看人在逋| 精品亚洲成国产av| 狠狠精品人妻久久久久久综合| 91成年电影在线观看| 亚洲五月色婷婷综合| 久久热在线av| 欧美人与性动交α欧美精品济南到| 在线观看一区二区三区激情| 欧美人与性动交α欧美软件| 国产精品成人在线| 99精品欧美一区二区三区四区| 少妇 在线观看| 高清在线国产一区| a在线观看视频网站| 亚洲视频免费观看视频| 久久久久国产精品人妻一区二区| 国产精品秋霞免费鲁丝片| 91精品三级在线观看| 一区二区三区四区激情视频| 国产精品av久久久久免费| 日韩 亚洲 欧美在线| 天堂8中文在线网| 欧美日韩成人在线一区二区| 午夜影院在线不卡| 夜夜骑夜夜射夜夜干| 波多野结衣一区麻豆| 久久中文看片网| 欧美日韩黄片免| 黄色怎么调成土黄色| 国产欧美日韩综合在线一区二区| 男女床上黄色一级片免费看| 男人添女人高潮全过程视频| 国产深夜福利视频在线观看| 免费久久久久久久精品成人欧美视频| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美一区二区综合| 亚洲精品日韩在线中文字幕| 亚洲av欧美aⅴ国产| 人成视频在线观看免费观看| 国产在线观看jvid| 中文字幕精品免费在线观看视频| 亚洲av国产av综合av卡| kizo精华| 老司机靠b影院| 91大片在线观看| 叶爱在线成人免费视频播放| 老司机在亚洲福利影院| 夜夜夜夜夜久久久久| 在线天堂中文资源库| 黑人巨大精品欧美一区二区蜜桃| 91精品伊人久久大香线蕉| av超薄肉色丝袜交足视频| 日韩 欧美 亚洲 中文字幕| 日日爽夜夜爽网站| 操美女的视频在线观看| 黑丝袜美女国产一区| 亚洲精品久久成人aⅴ小说| 一级a爱视频在线免费观看| 9色porny在线观看| bbb黄色大片| 中国美女看黄片| 日韩熟女老妇一区二区性免费视频| 久久久国产欧美日韩av| 美女视频免费永久观看网站| 午夜福利影视在线免费观看| 久久影院123| 老司机午夜福利在线观看视频 | 久久午夜综合久久蜜桃| 成年人黄色毛片网站| 男女床上黄色一级片免费看| 女人被躁到高潮嗷嗷叫费观| 五月天丁香电影| 欧美久久黑人一区二区| 亚洲国产中文字幕在线视频| 99国产精品一区二区蜜桃av | 中文字幕另类日韩欧美亚洲嫩草| 久热这里只有精品99| 黄色视频在线播放观看不卡| 一级毛片女人18水好多| 亚洲欧洲日产国产| 夜夜夜夜夜久久久久| 国产黄色免费在线视频| 国产免费现黄频在线看| 免费少妇av软件| 少妇粗大呻吟视频| 欧美日韩精品网址| 国产成人免费无遮挡视频| 99国产精品99久久久久| 国产成人av激情在线播放| 亚洲国产中文字幕在线视频| 亚洲成国产人片在线观看| 女警被强在线播放| 精品国产乱码久久久久久小说| 欧美日韩视频精品一区| 国产精品久久久久久人妻精品电影 | 青春草视频在线免费观看| 9热在线视频观看99| 老司机亚洲免费影院| 91九色精品人成在线观看| 丁香六月欧美| 高潮久久久久久久久久久不卡| 可以免费在线观看a视频的电影网站| 性少妇av在线| 99久久国产精品久久久| 国产一区二区三区av在线| 老鸭窝网址在线观看| 国产成人精品久久二区二区91| 亚洲av美国av| 午夜福利,免费看| 在线永久观看黄色视频| 国产精品影院久久| 婷婷色av中文字幕| 777久久人妻少妇嫩草av网站| 精品国产国语对白av| 久久人人爽av亚洲精品天堂| 欧美另类一区| 亚洲精品国产色婷婷电影| 老司机福利观看| 久久久久久亚洲精品国产蜜桃av| 久久精品久久久久久噜噜老黄| 如日韩欧美国产精品一区二区三区| 日本av免费视频播放| 精品久久久久久久毛片微露脸 | 在线精品无人区一区二区三| 午夜影院在线不卡| 国产真人三级小视频在线观看| 精品人妻一区二区三区麻豆| 久久综合国产亚洲精品| 99久久国产精品久久久| 老司机午夜十八禁免费视频| 亚洲免费av在线视频| 久9热在线精品视频| 欧美日韩视频精品一区| 丝袜在线中文字幕| 亚洲五月色婷婷综合| 黄色视频在线播放观看不卡| 在线观看舔阴道视频| 91精品国产国语对白视频| 国产精品九九99| 久久九九热精品免费| 亚洲伊人色综图| a 毛片基地| 在线十欧美十亚洲十日本专区| 久久精品国产综合久久久| 91精品国产国语对白视频| 9191精品国产免费久久| 亚洲欧美日韩高清在线视频 | 人成视频在线观看免费观看| 亚洲av欧美aⅴ国产| 99国产精品一区二区蜜桃av | 亚洲av日韩精品久久久久久密| 在线观看人妻少妇| 亚洲国产成人一精品久久久| 少妇精品久久久久久久| 国产在线观看jvid| 搡老岳熟女国产| 日韩视频在线欧美| 亚洲伊人色综图| 97在线人人人人妻| 国产精品麻豆人妻色哟哟久久| 美女中出高潮动态图| 蜜桃在线观看..| 成年人免费黄色播放视频| 视频在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 成人国语在线视频| 夫妻午夜视频| 日本av手机在线免费观看| 啦啦啦视频在线资源免费观看| 中文字幕av电影在线播放| 国产av国产精品国产| 亚洲色图综合在线观看| av在线播放精品| 国产99久久九九免费精品| 亚洲精品久久午夜乱码| 亚洲中文日韩欧美视频| 欧美性长视频在线观看| 亚洲午夜精品一区,二区,三区| 曰老女人黄片| 啦啦啦啦在线视频资源| 久久国产精品影院| 欧美性长视频在线观看| 老司机在亚洲福利影院| 国产精品熟女久久久久浪| 黄网站色视频无遮挡免费观看| 国产精品久久久久久精品古装| 女人高潮潮喷娇喘18禁视频| 1024香蕉在线观看| 两个人免费观看高清视频| 亚洲久久久国产精品| av欧美777| videos熟女内射| 国产亚洲一区二区精品| 曰老女人黄片| 天堂俺去俺来也www色官网| 每晚都被弄得嗷嗷叫到高潮| 18禁观看日本| 国产精品久久久久久精品古装| 午夜精品国产一区二区电影| 国产日韩欧美亚洲二区| 久久久精品免费免费高清| 亚洲第一青青草原| 18禁国产床啪视频网站| www.av在线官网国产| 午夜福利免费观看在线| 男女之事视频高清在线观看| 中国美女看黄片| 香蕉国产在线看| 69精品国产乱码久久久| 日本一区二区免费在线视频| 亚洲国产av新网站| 国产精品亚洲av一区麻豆| 三级毛片av免费| 精品熟女少妇八av免费久了| 日韩制服骚丝袜av| 在线av久久热| 亚洲精品乱久久久久久| 国产成人欧美在线观看 | 色婷婷久久久亚洲欧美| www日本在线高清视频| 777久久人妻少妇嫩草av网站| 日韩一区二区三区影片| 午夜免费成人在线视频| 99国产精品一区二区三区| 国产免费一区二区三区四区乱码| 在线 av 中文字幕| 女人被躁到高潮嗷嗷叫费观| av欧美777| 一个人免费在线观看的高清视频 | 18在线观看网站| 成年人黄色毛片网站| 国产成人免费观看mmmm| 精品国产乱子伦一区二区三区 | 国产一区二区三区综合在线观看| 亚洲精品美女久久av网站| 国产精品1区2区在线观看. | 亚洲国产精品成人久久小说| 一级毛片电影观看| av国产精品久久久久影院| 动漫黄色视频在线观看| 在线观看免费视频网站a站| 最近最新免费中文字幕在线| 欧美国产精品一级二级三级| 久久人人爽人人片av| 深夜精品福利| 视频区图区小说| 99九九在线精品视频| 免费在线观看影片大全网站| 国产xxxxx性猛交| 日本猛色少妇xxxxx猛交久久| 亚洲 欧美一区二区三区| 成人av一区二区三区在线看 | 无遮挡黄片免费观看| 日本91视频免费播放| 精品视频人人做人人爽| 国产免费现黄频在线看| 日韩大码丰满熟妇| 国产欧美日韩一区二区精品| 成人国语在线视频| 免费在线观看日本一区| 午夜激情av网站| 亚洲熟女毛片儿| 国产成人系列免费观看| 99精品欧美一区二区三区四区| 天堂中文最新版在线下载| 国产成人av教育| 久久久精品94久久精品| 老司机影院成人| 精品亚洲成国产av| 99国产精品免费福利视频| 高清黄色对白视频在线免费看| 一区二区av电影网| 久久99热这里只频精品6学生| 咕卡用的链子| 日韩一区二区三区影片| 欧美精品高潮呻吟av久久| 日韩三级视频一区二区三区| 一区二区三区精品91| 人妻久久中文字幕网| 国精品久久久久久国模美| a 毛片基地| 天天影视国产精品| 国产三级黄色录像| 丝袜人妻中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品一区三区| 欧美激情久久久久久爽电影 | 女人久久www免费人成看片| 亚洲成国产人片在线观看| 色播在线永久视频| 在线亚洲精品国产二区图片欧美| 欧美午夜高清在线| 天天操日日干夜夜撸| av片东京热男人的天堂| 精品一区二区三区av网在线观看 | 少妇的丰满在线观看| 色视频在线一区二区三区| 久久久久久免费高清国产稀缺| 80岁老熟妇乱子伦牲交| 亚洲av日韩在线播放| 在线观看免费午夜福利视频| 欧美久久黑人一区二区| 免费日韩欧美在线观看| 飞空精品影院首页| 捣出白浆h1v1| 欧美日韩视频精品一区| 淫妇啪啪啪对白视频 | 精品亚洲成a人片在线观看| av有码第一页| 成人黄色视频免费在线看| 欧美精品av麻豆av| a级毛片黄视频| 老汉色av国产亚洲站长工具| 别揉我奶头~嗯~啊~动态视频 | 在线天堂中文资源库| 各种免费的搞黄视频| 妹子高潮喷水视频| 国产免费av片在线观看野外av| 51午夜福利影视在线观看| 国产精品偷伦视频观看了| 脱女人内裤的视频| 一个人免费看片子| 国内毛片毛片毛片毛片毛片| 99香蕉大伊视频| 久久久久久久精品精品| 建设人人有责人人尽责人人享有的| 另类亚洲欧美激情| 日本a在线网址| 亚洲性夜色夜夜综合| 国产欧美日韩综合在线一区二区| 女性被躁到高潮视频| 免费少妇av软件| 黄频高清免费视频| 777米奇影视久久| 精品久久久精品久久久| 午夜免费成人在线视频| 黑丝袜美女国产一区| 黄色片一级片一级黄色片| 精品国内亚洲2022精品成人 | 久久久欧美国产精品| 男女边摸边吃奶| 久久久久精品国产欧美久久久 | 精品国产一区二区三区四区第35| 高清黄色对白视频在线免费看| 中文字幕色久视频| 亚洲自偷自拍图片 自拍| 在线观看一区二区三区激情| 少妇裸体淫交视频免费看高清 | 曰老女人黄片| 涩涩av久久男人的天堂| 国产成人精品在线电影| 国产国语露脸激情在线看| 日日夜夜操网爽| 亚洲国产成人一精品久久久| 亚洲精品国产精品久久久不卡| 在线观看免费午夜福利视频| 免费在线观看日本一区| 亚洲国产成人一精品久久久| 99国产精品免费福利视频| 欧美黄色淫秽网站| 国产精品国产三级国产专区5o| 精品人妻熟女毛片av久久网站| 中文欧美无线码| 久久青草综合色| av天堂在线播放| 精品免费久久久久久久清纯 | 99久久人妻综合| 女性生殖器流出的白浆| 亚洲免费av在线视频| 9热在线视频观看99| 国产精品久久久人人做人人爽| 国产精品成人在线| 久久99热这里只频精品6学生| 亚洲全国av大片| 大陆偷拍与自拍| 日韩制服骚丝袜av| 精品久久久精品久久久| 天天影视国产精品| 久久免费观看电影| 欧美大码av| 老汉色∧v一级毛片| 啦啦啦 在线观看视频| 两人在一起打扑克的视频| 亚洲欧美一区二区三区久久| tube8黄色片| 色精品久久人妻99蜜桃| 黄色a级毛片大全视频| 自线自在国产av| 欧美xxⅹ黑人| 欧美中文综合在线视频| 丰满饥渴人妻一区二区三| 成年人免费黄色播放视频| e午夜精品久久久久久久| 国产精品久久久久久精品电影小说| 在线永久观看黄色视频| 国产av又大| 亚洲av片天天在线观看| 久久亚洲精品不卡| 国产黄频视频在线观看| 精品少妇内射三级| 男人添女人高潮全过程视频| 久久久久国内视频| 久久国产亚洲av麻豆专区| 黑人巨大精品欧美一区二区蜜桃| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 国产一级毛片在线| 黑丝袜美女国产一区| 美女扒开内裤让男人捅视频| 婷婷成人精品国产| 欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦 在线观看视频| 老司机午夜十八禁免费视频| 成年人午夜在线观看视频| 久久精品国产亚洲av香蕉五月 | 99精国产麻豆久久婷婷| 曰老女人黄片| 亚洲一卡2卡3卡4卡5卡精品中文| 18禁观看日本| 久久精品人人爽人人爽视色| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品久久二区二区免费| 美女午夜性视频免费| 久久九九热精品免费| 一区二区三区乱码不卡18| av视频免费观看在线观看| 欧美日韩视频精品一区| 我要看黄色一级片免费的| 国产又色又爽无遮挡免| 一进一出抽搐动态| 热99国产精品久久久久久7| 亚洲男人天堂网一区| 成人国产av品久久久| 人人妻人人澡人人爽人人夜夜| 精品国产一区二区三区久久久樱花| 久久久精品94久久精品| 婷婷色av中文字幕| 狠狠婷婷综合久久久久久88av| 制服人妻中文乱码| 久久精品国产a三级三级三级| 男女边摸边吃奶| 男男h啪啪无遮挡| av欧美777| 久久精品久久久久久噜噜老黄| 黄色视频在线播放观看不卡| 国产精品九九99| √禁漫天堂资源中文www| 超碰97精品在线观看| 黑丝袜美女国产一区| 男女国产视频网站| 少妇粗大呻吟视频| 一本一本久久a久久精品综合妖精| 一区在线观看完整版| 人人澡人人妻人| 一边摸一边做爽爽视频免费| 免费高清在线观看日韩| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 精品福利永久在线观看| 亚洲成人手机| 欧美午夜高清在线| 成人国语在线视频| 亚洲精品国产区一区二| 大型av网站在线播放| 啦啦啦视频在线资源免费观看| 久热这里只有精品99| 国产国语露脸激情在线看| av免费在线观看网站| 国产免费福利视频在线观看| 精品久久久精品久久久| 久久中文看片网| 狂野欧美激情性xxxx| 999久久久精品免费观看国产| 亚洲欧美精品自产自拍| 日日爽夜夜爽网站| 精品人妻在线不人妻| 国产成人av教育| 激情视频va一区二区三区| 亚洲欧美一区二区三区黑人| 日日摸夜夜添夜夜添小说| 国产淫语在线视频| 亚洲伊人久久精品综合| 12—13女人毛片做爰片一| 久久99一区二区三区| 久久国产亚洲av麻豆专区| 少妇被粗大的猛进出69影院| 日韩一卡2卡3卡4卡2021年| 动漫黄色视频在线观看| 最近最新中文字幕大全免费视频| 50天的宝宝边吃奶边哭怎么回事| 美女中出高潮动态图| 亚洲国产av新网站| 丰满饥渴人妻一区二区三| 51午夜福利影视在线观看| 青春草亚洲视频在线观看| 男人添女人高潮全过程视频| 精品一区二区三区四区五区乱码| 久久人人爽人人片av| 亚洲精品中文字幕一二三四区 | 中国国产av一级| 亚洲第一av免费看| 视频区图区小说| 999久久久精品免费观看国产| 母亲3免费完整高清在线观看| 一区二区三区乱码不卡18| 99re6热这里在线精品视频| 欧美少妇被猛烈插入视频| av天堂在线播放| 久9热在线精品视频| 男人爽女人下面视频在线观看| 成人免费观看视频高清| 久久影院123| 午夜影院在线不卡| 日本一区二区免费在线视频| 天堂俺去俺来也www色官网| 国产精品久久久久久人妻精品电影 | 亚洲七黄色美女视频| 99国产精品99久久久久| 国产91精品成人一区二区三区 | 国产成人免费观看mmmm| 欧美精品啪啪一区二区三区 | 亚洲黑人精品在线| 国产成+人综合+亚洲专区| 国产一区二区激情短视频 | 精品乱码久久久久久99久播| 人妻久久中文字幕网| 免费黄频网站在线观看国产| 视频区图区小说| 国产欧美日韩一区二区三区在线| 欧美日韩成人在线一区二区| 免费在线观看视频国产中文字幕亚洲 | 电影成人av| 欧美性长视频在线观看| 在线天堂中文资源库| 一级片'在线观看视频| 亚洲欧美激情在线| 叶爱在线成人免费视频播放| 亚洲色图综合在线观看| 欧美午夜高清在线| av网站在线播放免费| 黑人巨大精品欧美一区二区mp4| 欧美激情 高清一区二区三区| 熟女少妇亚洲综合色aaa.| 秋霞在线观看毛片| 久久久水蜜桃国产精品网| 制服人妻中文乱码| 午夜福利一区二区在线看| 免费久久久久久久精品成人欧美视频| 亚洲av成人不卡在线观看播放网 | 搡老熟女国产l中国老女人| 啦啦啦 在线观看视频| 黄色怎么调成土黄色| 国产淫语在线视频| 极品人妻少妇av视频| 老熟女久久久| 精品久久久久久电影网| 999久久久国产精品视频| 波多野结衣一区麻豆| 97人妻天天添夜夜摸| 亚洲黑人精品在线| 亚洲专区字幕在线| 国产亚洲精品久久久久5区| 夜夜骑夜夜射夜夜干| 免费观看人在逋| 黑人操中国人逼视频| 国产男女超爽视频在线观看| 啦啦啦中文免费视频观看日本| 欧美老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 国产一卡二卡三卡精品| 精品第一国产精品| 91麻豆精品激情在线观看国产 | 国产97色在线日韩免费| 女人高潮潮喷娇喘18禁视频| 亚洲国产av新网站| 国产欧美日韩综合在线一区二区| 999精品在线视频| 亚洲天堂av无毛| 丰满饥渴人妻一区二区三| 免费高清在线观看视频在线观看| 少妇 在线观看| 久久精品成人免费网站| 日韩一区二区三区影片| 国产成人影院久久av| 日韩中文字幕欧美一区二区| 大香蕉久久成人网| 黄色 视频免费看| 亚洲人成电影免费在线| 91字幕亚洲| 精品一区二区三区四区五区乱码| 久久精品国产a三级三级三级| 大片电影免费在线观看免费| 男人舔女人的私密视频| bbb黄色大片| 亚洲av成人不卡在线观看播放网 | 性色av一级| 18禁黄网站禁片午夜丰满| 女人高潮潮喷娇喘18禁视频| 爱豆传媒免费全集在线观看|