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

    俯仰角對鴨式布局飛機搖滾運動的影響與機理分析

    2017-04-28 03:46:24魏德宸史志偉昂海松付軍泉朱佳晨
    空氣動力學學報 2017年2期
    關鍵詞:主翼迎角前體

    魏德宸, 史志偉, 昂海松, 付軍泉, 朱佳晨

    (南京航空航天大學 航空宇航學院, 江蘇 南京 210016)

    ?

    俯仰角對鴨式布局飛機搖滾運動的影響與機理分析

    魏德宸, 史志偉*, 昂海松, 付軍泉, 朱佳晨

    (南京航空航天大學 航空宇航學院, 江蘇 南京 210016)

    為了研究鴨式布局飛機的搖滾運動,設計了一種包括鴨翼、脊型前體、邊條翼、主翼和垂尾的模型,進行滾轉(zhuǎn)自由度釋放、靜態(tài)測力、動導數(shù)試驗和煙線流場顯示多種技術手段相結(jié)合的風洞試驗。首先為了得到搖滾從發(fā)生到消失的全過程,進行了俯仰角范圍為12° ~52°、間隔1°測量的滾轉(zhuǎn)自由度釋放試驗;結(jié)果表明對應于不同的俯仰角,鴨翼布局飛機的搖滾運動也具有不同的性質(zhì),其平衡點和運動形式均發(fā)生變化。接著靜態(tài)測力和動導數(shù)試驗證實:在翼體結(jié)構的多渦系影響下,搖滾可在零度平衡點和非零平衡點位置處發(fā)生;且運動可為極限環(huán)和非極限環(huán)形式。最后通過流場顯示,分析了在不同迎角時鴨式布局飛機形成翼體搖滾的可能的流動機理。

    鴨翼;翼體搖滾;非零平衡點;非極限環(huán);流場顯示

    0 引 言

    隨著人們對非線性渦升力認識的不斷深化和先進控制技術的應用,現(xiàn)代飛行器的可控迎角不斷增大。然而大迎角飛行時,常常會誘發(fā)非指令的自激運動,搖滾正是其中一種嚴重的動態(tài)表現(xiàn)形式。這種以繞體軸滾轉(zhuǎn)振蕩為主的自激運動不僅縮小了飛行包線,限制了戰(zhàn)斗機的機動性和敏捷性, 而且嚴重影響了安定性和操縱性, 甚至引起結(jié)構強度問題導致飛機解體。

    Ericsson根據(jù)飛行器產(chǎn)生搖滾的機理對搖滾現(xiàn)象分為三類[1-4]: (1) 細長體搖滾; (2) 常規(guī)機翼搖滾; (3) 翼身組合體搖滾。翼身組合體搖滾是最接近真實飛行器搖滾的類型,然而由于問題的復雜性等原因,翼體搖滾相應研究開展的較少。隸屬于北約組織(NATO)的航空研究與發(fā)展顧問小組(AGARD)發(fā)起的一項針對解決飛行器高機動飛行的項目中,翼體搖滾的研究和預測作為其中重要組成部分,開展了針對編號為“WG16A”翼身組合體模型的風洞試驗研究[5]。Guglieri則通過前體形狀分別為圓柱和圓錐不同構型的80° 后掠三角翼的風洞試驗,研究前體渦的改變對搖滾特性的影響[6]。作為NASA和美國海軍聯(lián)合牽頭的翼突然失速(AWS)計劃的一部分[7-9],自由搖滾(FTR)和測力試驗被應用在多種真實飛行器縮比模型的風洞試驗中。其中Owens等研究了前后緣襟翼角度變化、渦流發(fā)生器等對F-35戰(zhàn)斗機搖滾的影響[10]。Green等研究了F/A-18戰(zhàn)斗機的機翼厚度、彎度等參數(shù)對搖滾運動的影響[11-12]。類似的AV-8B和F-16C的縮比模型也通過風洞試驗來用于驗證和預測實際飛行中可能出現(xiàn)的搖滾現(xiàn)象[13]。近年國內(nèi)多家單位也開展了翼體搖滾的研究工作, 張涵信等針對翼體搖滾運動中的靜穩(wěn)定性和動穩(wěn)定性進行了大量研究和分析,提出了靜態(tài)和動態(tài)穩(wěn)定性分析判據(jù)[14-16]。馬寶峰、榮臻等重點分析了前體渦在細長前體結(jié)構的翼身組合體搖滾中的誘導作用[17-18]。 張杰等針對一種先進作戰(zhàn)飛機進行了跨聲速的自由滾轉(zhuǎn)試驗,并對橫滾運動進行了量化的評估[19]。孫海生和姜裕標研究了一種戰(zhàn)斗機模型的搖滾風洞試驗,結(jié)果表明單獨前體和翼身組合體結(jié)構均可產(chǎn)生搖滾;且翼身組合體在不同迎角范圍,具有不同的搖滾特征[20]。王兵等開展了尖拱旋成體和小后掠梯形機翼組合體在人工轉(zhuǎn)捩條件下的搖滾研究,擴展了低速風洞中進行搖滾試驗的雷諾數(shù)范圍[21]。但目前國內(nèi)發(fā)表的翼體搖滾文獻多基于細長旋成體和機翼的簡單結(jié)構模型,國外公開文獻也缺乏針對復雜結(jié)構飛行器特別是鴨式布局飛機所引發(fā)的翼體搖滾的流場研究和機理分析。

    1 試驗設備及模型

    1.1 試驗風洞

    試驗是在南京航空航天大學開口回流式風洞中

    進行的,風洞實驗段長×寬×高:1.7 m×1.5 m×1.0 m,最大風速35 m/s,最小穩(wěn)定風速3 m/s。風洞流場品質(zhì)良好,俯仰方向氣流偏角|α|≤0.5°,偏航方向氣流偏角|β|≤0.5°,湍流度ε≤0.07%。

    1.2 試驗模型

    試驗所采用的為鴨式布局飛機模型,全長0.66 m,主翼展長0.5 m。前體為脊型截面;鴨翼具有-20° 的固定偏角(前緣向下),在主翼迎角較大時,鴨翼的實際迎角較小;鴨翼與主翼前緣后掠角均為Λ=50°。具體參數(shù)如圖1所示。

    圖1 試驗模型Fig.1 Test model

    1.3 試驗設備

    試驗系統(tǒng)臺為圖2所示的五自由度并聯(lián)機構,可實現(xiàn)模型俯仰、偏航和滾轉(zhuǎn)三個方向的運動,以及這三個方向上的組合運動[22]。本次試驗在此機構上進行了靜態(tài)測力與動導數(shù)測量,所用數(shù)據(jù)采集卡的分辨率為16位, 最大采樣速率250 kS/s。將天平支桿更換為自由搖滾裝置可進行自由搖滾試驗[23]。自由搖滾裝置如圖3所示,滾轉(zhuǎn)角測量精度為0.1°。

    圖2 五自由度并聯(lián)機構Fig.2 The 5 DOF dynamic test mechanism

    圖3 自由搖滾裝置[23]Fig.3 Free-to-roll test sting[23]

    2 滾轉(zhuǎn)自由度釋放與穩(wěn)定性分析

    2.1 翼身組合體自由度釋放

    搖滾本質(zhì)上是多自由度耦合的運動,但其運動形式為以繞體軸滾轉(zhuǎn)振蕩為主的自維持動態(tài)運動,滾轉(zhuǎn)是最主要的自由度,有文獻[24-25]指出,一定情況下,可采用單滾轉(zhuǎn)自由度假設來研究飛行器的動態(tài)搖滾特性。為分析鴨式布局飛機翼體搖滾的特性,進行了θ=12°~52° 大俯仰角范圍的自由搖滾試驗,包括了搖滾從發(fā)生到消失的全過程。圖4給出了鴨式布局飛機不同俯仰角θ范圍內(nèi)典型的滾轉(zhuǎn)角φ時間歷程,試驗風速為V=25 m/s。

    從搖滾運動的形式分析:如圖4(a)所示,θ=12° 時滾轉(zhuǎn)運動呈現(xiàn)微小抖動。原因可能為小迎角下翼體表面以附著流動為主,搖滾尚未發(fā)生,非均勻來流中小擾動的作用使得滾轉(zhuǎn)角發(fā)生變化[26]。圖4(b~d)所示,θ=23°、31° 和35° 時搖滾運動頻率和振幅存在明顯波動,伴有間歇、猝發(fā)等特征??赡艿脑驗榍绑w渦、鴨翼前緣渦、邊條渦和主翼前緣渦的強弱變化和相互干擾[13],使得搖滾運動不是規(guī)則的等幅振蕩(不同于常見的細長三角翼搖滾),而是非極限環(huán)運動,其機理將在后續(xù)的流場測量中進行分析。在圖4(e)中,θ=45° 時搖滾呈現(xiàn)中小迎角范圍內(nèi)未出現(xiàn)的準極限環(huán)運動形式,這時的搖滾運動可能類似于細長體搖滾。此時對搖滾起主要作用的可能是脊型前體產(chǎn)生的前體渦,其他渦系基本完全破裂。這時的主機翼也許只是提供滾轉(zhuǎn)力矩所需力的作用面;或去除主機翼,機身相當于細長體模型,搖滾特征參數(shù)會發(fā)生變化,但極限環(huán)形式的搖滾現(xiàn)象仍可能發(fā)生[20,27]。而如圖4(f)所示,θ=52° 時滾轉(zhuǎn)運動再次呈現(xiàn)無規(guī)則地微小抖動,搖滾運動消失。此時的運動狀態(tài)和θ=12° 相似,但機理完全不同,可能的原因為此迎角下包括前體渦在內(nèi)的多個渦系已完全破裂,翼體表面被湍流覆蓋,不能產(chǎn)生較大的氣動力矩。搖滾運動形式隨俯仰角變化的機理分析具體見第3節(jié)。

    (a) θ=12°

    (c) θ=31°

    (e) θ=45°

    (b) θ=23°

    (d) θ=35°

    (f) θ=52°

    圖4 鴨式布局飛機不同俯仰角下的典型搖滾時間歷程
    Fig.4 Typical Roll angle time history for the canard-configuration aircraft model at different pitch angles

    從平衡點位置分析,此鴨式布局飛機隨俯仰角變化,翼體搖滾從發(fā)生到消失的過程中滾轉(zhuǎn)中心角也隨之改變的。除了常規(guī)的φmean≈0°,在滾轉(zhuǎn)一側(cè)也具有多個非零平衡點。如圖4所示θ=12° 時,初始時平衡點為靜態(tài)平衡點,中心角φmean≈0°;θ=23° 時,模型處于翼體搖滾狀態(tài)(證明見2.3節(jié)),此時的搖滾平衡點為φmean≈-12°;隨著俯仰角的增加,如圖4(c)所示θ=31° 時,模型仍處于翼體搖滾狀態(tài),但此時的平衡點回到φmean≈0°;θ=35° 時, 鴨式布局飛機圍繞更大的滾轉(zhuǎn)中心角φmean≈26° 做搖滾運動;俯仰角度增加,如圖4(e)所示θ=45° 時,搖滾平衡點再次回到φmean≈0°;俯仰角繼續(xù)增加,如圖4(f)所示θ=52° 時,搖滾現(xiàn)象逐漸消失,0° 搖滾平衡點變?yōu)?° 靜態(tài)平衡點。

    2.2 滾轉(zhuǎn)中心角隨俯仰角變化規(guī)律分析

    為了更好的分析搖滾發(fā)展規(guī)律,對自由搖滾試驗得到的滾轉(zhuǎn)角度序列進行數(shù)值分析。處理公式如下所示[28],其中φmean為滾轉(zhuǎn)中心角;φi為時間歷程下的滾轉(zhuǎn)角。

    如圖4(b)、(d)所示搖滾運動可在滾轉(zhuǎn)對稱兩側(cè)發(fā)生,為便于分析搖滾發(fā)展規(guī)律,將中心角取絕對值得到平衡點位置的絕對值隨俯仰角變化曲線,如圖5所示。

    圖5 鴨式布局飛機滾轉(zhuǎn)中心角隨俯仰角變化曲線Fig.5 Mean roll angles change with increasing pitch angle for canard-configuration aircraft

    如圖5所示,θ=12°~18° 時,初始時φmean≈0°為靜態(tài)平衡點;在θ=19°~26° 范圍,隨著俯仰角的增加,搖滾在這一階段產(chǎn)生,平衡點位置由φmean≈0° 逐漸變動到φmean≈16°;θ=27°~32° 時,搖滾平衡點再次回到φmean≈0°;θ=33°~39° 時,搖滾平衡點又一次位于非零位置;θ=40°~50° 范圍內(nèi)搖滾平衡點再次回到零度位置;θ≥51° 時,搖滾消失,模型在滾轉(zhuǎn)角零度位置附近穩(wěn)定。通過θ變化范圍40°、間隔Δ=1° 的滾轉(zhuǎn)自由度釋放試驗,可知在不同的俯仰角范圍,翼體搖滾的影響渦系可能也不同,具體分析見第3節(jié)。

    2.3 翼體搖滾的穩(wěn)定性分析

    翼體搖滾為單自由度滾轉(zhuǎn)運動時,其滾轉(zhuǎn)運動公式可寫為:

    試驗風速仍為25m/s,通過兩側(cè)滾轉(zhuǎn)位置的靜態(tài)測力和動導數(shù)試驗,獲得了圖4中非極限環(huán)和極限環(huán)兩種搖滾運動形式典型俯仰角時的滾轉(zhuǎn)力矩系數(shù)Cl和滾轉(zhuǎn)阻尼導數(shù)Clp,見圖6和圖7。

    圖6 不同俯仰角的滾轉(zhuǎn)力矩系數(shù)ClFig.6 Static roll moment coefficients at different pitch angles

    圖7 不同俯仰角的滾轉(zhuǎn)阻尼導數(shù)ClpFig.7 Dynamic damping derivatives at different pitch angles

    首先對不同俯仰角時,φ=0° 的滾轉(zhuǎn)力矩系數(shù)相對滾轉(zhuǎn)角的斜率進行總結(jié)。滾轉(zhuǎn)角φ=0°附近滾轉(zhuǎn)力矩系數(shù)Cl均過零點,力矩是平衡的。如圖6所示,θ=23° 時斜率為正,說明模型在φ=0° 是靜態(tài)不穩(wěn)定的,搖滾不能在此中心角度發(fā)生;θ=31° 時斜率為負,說明模型在φ=0° 靜態(tài)穩(wěn)定,該滾轉(zhuǎn)角為翼體搖滾的一個靜態(tài)平衡點或者搖滾的平衡點;同理θ=35° 時零度滾轉(zhuǎn)角的Clp斜率再次變?yōu)檎?;?45° 時其斜率又為負值。Clp曲線斜率在φ=0° 時的正負變化和圖4的零度平衡點和非零平衡點的變化是一致的,理論分析結(jié)果和試驗結(jié)果重合。再觀測圖6中滾轉(zhuǎn)力矩系數(shù)為零時對應的非零滾轉(zhuǎn)角,θ=23° 且φ≈±12°時斜率為負,θ=35° 且φ≈±26° 時斜率也為負。最后結(jié)合圖7動導數(shù)結(jié)果,當靜態(tài)力矩斜率為負且動導數(shù)為正值,對應滾轉(zhuǎn)角才為搖滾平衡點,可知φ≈±12°、0°、±26°、0° 依次為θ=23°、31°、35°、45° 時的搖滾平衡點。

    3 流場顯示與機理分析

    鴨式布局飛機表面的渦系分布如圖8所示。其中脊形前體渦與鴨翼前緣渦同向旋轉(zhuǎn)且位置較近,在鴨翼翼面上繞合在一起,形成前體與鴨翼前緣卷繞渦;邊條渦和主翼前緣渦繞合在一起,形成了邊條與主翼前緣卷繞渦。前體與鴨翼前緣卷繞渦位于內(nèi)側(cè),且沿流向發(fā)展距離較長,故在主翼面的位置較高;邊條與主翼前緣卷繞渦受到同向旋轉(zhuǎn)的前體與鴨翼前緣卷繞渦的外推作用,位于主翼外側(cè),且沿流向發(fā)展距離短,渦心位置也較低[31]。

    圖8 鴨式布局飛機的渦系分布Fig.8 Vortex structures over canard-configuration aircraft

    為了研究不同俯仰角時,鴨式布局飛機翼體搖滾具有不同性質(zhì)的原因,應用煙線和激光片光技術在如圖1截面S所示的主翼1/2弦線位置,φ=0° 時進行流場觀測。

    如圖9(a)所示,距離鴨式布局飛機中軸線距離較近的渦為脊型前體和鴨翼前緣形成的前體與鴨翼前緣卷繞渦,在圖示中將其命名為“Vortex A”; 距離軸線更遠的渦則為邊條和主翼前緣形成的邊條與主翼前緣卷繞渦,在圖示中將其命名為“Vortex B”。

    在θ=12° 時,前緣后掠角為50° 的下偏鴨翼和主翼產(chǎn)生的前緣渦強度較弱;而脊型前體和邊條均可產(chǎn)生穩(wěn)定的渦流;故此時翼身組合體兩側(cè)的渦流分布較為對稱,飛行器處于靜態(tài)平衡狀態(tài)。θ=23° 時,脊型前體渦仍然穩(wěn)定[32-33],而此時鴨翼的實際迎角為23°-20°=3°,所產(chǎn)生的前緣渦強度弱,故前體與鴨翼前緣卷繞渦如圖9(b)所示仍較為對稱;主翼前緣渦變得不對稱,使得兩側(cè)的邊條與主翼前緣卷繞渦的對稱性遭到破壞,如圖9(b)所示機身左側(cè)遠離中軸線的渦位置較右側(cè)高,φ≈0° 時模型不再靜態(tài)穩(wěn)定。鴨式布局模型在如圖4(b)所示的φ≈-12° 中心角產(chǎn)生翼體搖滾的主要因素,可能為較弱的前體與鴨翼前緣卷繞渦和不對稱性的邊條與主翼前緣卷繞渦的共同影響。

    (a) θ=12°

    (b) θ=23°

    (c) θ=31°

    (d) θ=35°

    (e) θ=45°

    (f) θ=52°

    θ=31° 時,如圖9(c)所示前體渦進一步增強使得前體與鴨翼前緣卷繞渦也較為穩(wěn)定;而此時由于主翼前緣渦的進一步失衡使得邊條與主翼前緣卷繞渦的不對稱性增加。在此迎角下?lián)u滾發(fā)生在如圖4(c)所示φ≈0° 的原因,可能為更強且較為對稱的前體與鴨翼前緣卷繞渦的約束[33]和不對稱的邊條與主翼前緣卷繞渦的擾動。

    θ=35° 時,由于此時鴨翼前緣渦的不對稱性增加使得前體與鴨翼前緣卷繞渦也較為不對稱;同時兩側(cè)的邊條與主翼前緣卷繞渦如圖9(d)所示差異也較大。此時搖滾發(fā)生在如圖4(d)所示φ≈26° 的原因,可能為均失去對稱性的前體與鴨翼前緣卷繞渦和邊條與主翼前緣卷繞渦共同所致。類似的,在Gresham等的研究工作中,發(fā)現(xiàn)與鴨翼和主翼后掠角度相同的50° 后掠三角翼具有較大的非零平衡點[34];且比較可知,對于本次試驗中采用的鴨式布局模型, 搖滾可在俯仰角更大時發(fā)生,可能的原因為脊型前體[32-33],以及邊條翼帶來的穩(wěn)定性作用[35-36]。

    θ=45° 時,此時的搖滾和細長體搖滾類似,如圖4(e)所示,不同之前俯仰角范圍下的非極限環(huán),運動已表現(xiàn)為準極限環(huán)形式,且搖滾中心角為φ≈0°;可能的原因為如圖9(e)所示,此時對翼體搖滾起主要作用的是脊型前體產(chǎn)生的前體渦,其他的渦系均已完全破裂,而機身相當于細長體模型,主機翼也許只是提供滾轉(zhuǎn)力距的作用面,類似的情況在Ross等在只保留機身的搖滾試驗中也有發(fā)生[37]。θ=52° 時,隨著俯仰角的繼續(xù)增大,如圖9(f)所示脊型前體渦也已完全破裂,全機渦的破裂使得搖滾現(xiàn)象消失,θ=45° 時的零度搖滾平衡點變?yōu)槿鐖D4(f)所示的零度靜態(tài)平衡點。

    4 結(jié) 論

    1) 隨俯仰角變化,由于多渦系的影響,鴨式布局飛機搖滾具有不同的平衡點和運動形式,可為非零小角度平衡點的非極限環(huán)運動、零度平衡點非極限環(huán)運動、非零大角度平衡點的非極限環(huán)運動和零度平衡點準極限環(huán)運動。

    2) 穩(wěn)定性分析表明,對于存在多渦結(jié)構的鴨式布局飛機,通過地面的測力和動導數(shù)試驗可以判斷和預測其搖滾運動。

    3)θ=19°~26°范圍,鴨式布局飛機的不對稱性的邊條與主翼前緣卷繞渦可能是促發(fā)搖滾的主要渦系,前體與鴨翼前緣卷繞渦具有一定穩(wěn)定作用,但渦強較弱,共同的作用使得搖滾在較小中心滾轉(zhuǎn)角發(fā)生;θ=27°~32° 時,前體渦強度的增加使得搖滾中心回到零度位置;θ=33°~39° 時,前體與鴨翼前緣卷繞渦不對稱性的增加使得搖滾平衡點位于較大的非零位置;θ=40°~50° 范圍內(nèi)其他渦均已破裂,不對稱的前體渦可能是引起搖滾的主要原因。

    4)隨著俯仰角的改變,一方面前體或主翼等結(jié)構所產(chǎn)生渦的相對強弱變化,形成了不同的搖滾促發(fā)流場;另一方面前緣渦、前體渦、邊條渦等渦系對于翼體搖滾或誘導或穩(wěn)定的不同作用,共同施加于翼身組合體。這兩方面因素的影響使得鴨式布局飛機在不同俯仰角具有不同性質(zhì)的搖滾運動。

    [1]Ericsson L E. Slender wing rock revisited[J]. Journal of Aircraft, 1993, 30(3): 352-356.

    [2]Ericsson L E. Wing rock analysis of slender delta wings, Review and extension[J]. Journal of Aircraft, 1995, 32(6): 1221-1226.

    [3]Ericsson L E. Wing rock generated by forebody vortices[J]. Journal of Aircraft, 1989, 26(2): 110-116.

    [4]Ericsson L E. Various sources of wing rock[J]. Journal of Aircraft, 1990, 27(6): 488-494.

    [5]North Atlantic Treaty Organization. Cooperative programme on dynamic wind tunnel experiments for manoevring aircraft. AGARD AR-305[R]. Brussels: AGARD, 1996.

    [6]Guglieri G. A comprehensive analysis of wing rock dynamics for slender delta wing configurations[J]. Nonlinear Dyn, 2012, 69: 1559-1575.

    [7]Hall R M, Woodson S H. Introduction to the abrupt wing stall program[J]. Journal of Aircraft, 2004, 41(3): 425-435.

    [8]Hall R M, Woodson S H, Chambers J R. Overview of the abrupt wing stall program[J]. Progress in Aerospace Sciences, 2004, 40(7): 417-452.

    [9]Hall R M. Introduction: Abrupt wing stall program, part 2[J]. Journal of Aircraft, 2005, 42(3): 577-577.

    [10]Owens D B, McConnell J K, Brandon J M, et al. Transonic free-to-roll analysis of the F-35 (joint strike fighter) aircraft[J]. Journal of Aircraft, 2006, 43(3): 608-615.

    [11]Green B E, Ott J D. FA-18C to E wing morphing study for the abrupt-wing-stall program. AIAA 2003-0925[R]. Reston: AIAA, 2003.

    [12]Green B E, Ott J D. FA-18C to E wing morphing study for the abrupt-wing-stall program[J]. Journal of Aircraft, 2005, 42(3): 617-626.

    [13]Owens D B, Capone F J, Hall R M, et al. Free-to-roll analysis of abrupt wing stall on military aircraft at transonic speeds. AIAA 2003-0750[R]. Reston: AIAA, 2003.

    [14]Zhang Hanxin, Liu Wei, Xie Yufei, et al. On the rocking motion and its dynamic evolution of a sweep delta wing[J]. Acta Aerodynamica Sinica, 2006, 24(1): 5-9. (in Chinese)張涵信, 劉偉, 謝昱飛, 等. 后掠三角翼的搖滾及其動態(tài)演化問題[J]. 空氣動力學學報, 2006, 24(1): 5-9.

    [15]Ye Youda, Zhao Zhongliang, Tian Hao, et al. The stability analysis of rolling motion of hypersonic vehicles and its validations[J]. Science China, 2014, 57(12): 2194-2204.

    [16]Liu Wei, Zhang Hanxin. Nonlinear dynamic analysis and simulation of wing rock for a slender wing[C]//Progress in Computation Fluid Dynamics-Proceeding of 12thChina Computational Fluid Dynamics Conference. Beijing: Aerodynamics Society of China, 2004: 370-376. (in Chinese)劉偉, 張涵信. 細長機翼搖滾的非線性動力學分析及數(shù)值模擬[C]//計算流體力學會議論文集. 北京: 中國空氣動力學會, 2004: 370-376.

    [17]Ma Baofeng, Deng Xueying, Rong Zhen, et al. The self-excited rolling oscillations induced by fore-body vortices[J]. Aerospace Science and Technology, 2015, 47: 299-313.

    [18]Rong Zhen, Deng Xueying, Wang Bing, et al. Experimental investigation on flow characteristics during double limit cycle oscillation induced by forebody vortex[J]. Journal of Experiments in Fluid Mechanics, 2010, 24(3): 20-24. (in Chinese)榮臻, 鄧學鎣, 王兵, 等. 前體渦誘導雙極限環(huán)搖滾流動特性的實驗研究[J]. 實驗流體力學, 2010, 24(3): 20-24.

    [19]Zhang Jie, Cai Yi, Wu Jiali, et al. Research on a free-to-roll transonic test capability[J]. Acta Aerodynamica Sinica, 2016, 34(5): 611-616. (in Chinese)張杰, 才義, 吳佳莉, 等. 跨聲速自由滾轉(zhuǎn)試驗技術研究[J]. 空氣動力學學報, 2016, 34(5): 611-616.

    [20]Sun Haisheng. Investigation on wing rock in low speed wind tunnel for a fighter configuration[J]. Experiments and Measurements in Fluid Mechanics, 2000, 14(4): 32-35. (in Chinese) 孫海生. 飛機機翼搖滾低速風洞實驗研究[J]. 流體力學實驗與測量, 2000, 14(4): 32-35.

    [21]Wang Bing, Deng Xueying, Ma Baofeng, et al. An artificial transition technique applied to wing rock induced by forebody vortex[J]. Acta Aerodynamica Sinica, 2010, 28(5): 525-535. (in Chinese) 王兵, 鄧學鎣, 馬寶峰, 等. 前體渦誘導機翼搖滾的人工轉(zhuǎn)捩技術研究[J]. 空氣動力學學報, 2010, 28(5): 525-535.

    [22]Shi Zhiwei, Li Luyang, Cheng Keming, et al. Wind-tunnel high attack angle dynamic testing device based on hybrid mechanism: China, CN101929915B[P]. 2011-11-30. (in Chinese)http://www.pss-system.gov.cn/sipopublicsearch/portal/uiIndex.shtml史志偉, 李鷺揚, 程克明, 等. 基于混聯(lián)機構的風洞大迎角動態(tài)實驗裝置: 中國, CN101929915B[P]. 2011-11-30. http://www.pss-system.gov.cn/sipopublicsearch/portal/uiIndex.shtml

    [23]Geng Xi, Shi Zhiwei, Cheng Keming, et al. A new hybrid mechanism for dynamic wind tunnel test of high maneuverable air vehicle[J]. Proc IMechE, Part G: J Aerospace Engineering, (Online) DOI: 10. 1177/0954410015620448

    [24]Nelson R C, Pelletier A. The unsteady aerodynamics of slender wings and aircraft undergoing large amplitude maneuvers[J]. Progress in Aerospace Sciences, 2003, 39(2): 185-248.

    [25]Katz J. Wing/vortex interactions and wing rock[J]. Progress in Aerospace Sciences, 1999, 35: 727-750.

    [26]Wei Longkun, Ma Baofeng. Wing rock induced by a hemisphere-cylinder forebody[J]. Journal of Aircraft, 2014, 51(2): 606-613.

    [27]Liu Wei, Yang Xiaoliang, Zhang Hanxin, et al. A review on investigations of wing rock problems under high angles of attack[J]. Advances in Mechanics, 2008, 38(2): 214-218. (in Chinese)劉偉, 楊小亮, 張涵信, 等. 大迎角運動時的機翼搖滾問題研究綜述[J]. 力學進展, 2008, 38(2): 214-218.

    [28]Tao Yang, Zhao Zhongliang, Yang Haiyong. Investigation on wing rock of wing-body configuration at high speed wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2011, 25(6): 45-48. (in Chinese) 陶洋, 趙忠良, 楊海泳. 翼身組合體搖滾特性高速試驗研究[J]. 實驗流體力學, 2011, 25(6): 45-48.

    [29]Takashi M, Shigeru Y, Yoshiaki N. The effect of leading-edge profile of self-induced oscillation of 45-degree delta wings. AIAA 2000-4004[R]. Reston: AIAA, 2000.

    [30]Takashi M, Shigeru Y, Yoshiaki N. Data-based modeling of non-slender wing rock. AIAA 2003-738[R]. Reston: AIAA, 2003.

    [31]Wei Dechen, Shi Zhiwei, Geng Xi, et al. Wind tunnel test for wing-body rock of canard-configuration aircraft[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(10), 3003-3010. (in Chinese)魏德宸, 史志偉, 耿璽, 等. 鴨式布局飛行器的翼體搖滾特性風洞試驗[J]. 航空學報, 2016, 37(10): 3003-3010.

    [32]Brandon J M, Nguyen L T. Experimental study of effects of forebody geometry on high angle of attack static and dynamic stability. AIAA-86-0331[R]. Reston: AIAA, 1986.

    [33]Brandon J M, Nguyen L T. Experimental study of effects of forebody geometry on high angle of attack stability[J]. Journal of Aircraft, 1988, 25(7): 591-597.

    [34]Vortex dynamics of free-to-roll slender and nonslender delta wings[J]. Journal of Aircraft, 2010, 47(1): 292-302.

    [35]Grismer, Deborah S, Robert C N. Aerodynamics of a double delta wing undergoing pitching motions with and without sideslip. AIAA-94-1855[R]. Reston: AIAA, 1994.

    [36]Grismer, Deborah S. An experimental study of the steady and unsteady aerodynamics of double delta wings for zero and non-zero sideslip[D]. Canada: The University of Notre Dame, 1995.

    [37]Ross A J, Nguyen L T. Some observations regarding wing-rock oscillations at high angles of attack. AIAA-88-4371[R]. Reston: AIAA, 1988.

    Effects of pitch angle on wing-body rock of canard-configuration aircraft

    Wei Dechen, Shi Zhiwei*, Ang Haisong, Fu Junquan, Zhu Jiachen

    (CollegeofAerospaceEngineering,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China)

    To study wing-body rock of canard-configuration aircrafts, a model for wind tunnel test was designed having a chine forebody, canard wings, strake wings,main wings, and a vertical fin. This model was tested by different methods including free rolling motion, static force balance test, dynamic derivative test, and smoke wire technique. First, the complete rock process was studied by free-to-roll experiment in the range of pitch angelθ=12°~52° andΔ=1°. The results show that, rolling motions of the canard-configuration aircraft have several different mean roll angles and motion forms vary with different pitch angles. Second, the results of the static force balance test and dynamic derivative test confirm that the rolling motion can occur around zero or nonzero equilibrium roll angles with limit-cycle forms or nolimit cycle ones due to the influence of multi-vortex structures. Finally, the flow visualization tests show that the main vortices generate the wing-body rock at different stages. The possible mechanism for the wing-body rock was also analyzed for the canard-configuration aircraft.

    canard wing; wing-body rock; nonzero equilibrium roll angle; non-limit cycle; flow visualization

    0258-1825(2017)02-0198-07

    2016-11-13;

    2016-12-28

    江蘇高校優(yōu)勢學科建設工程資助項目

    魏德宸(1986-),男,江蘇徐州人,博士研究生,研究方向:飛行器大迎角氣動特性,非定??諝鈩恿W. E-mail: weimoving@163.com

    史志偉*(1972-),男,教授. E-mail: szwam@nuaa.edu.cn

    魏德宸, 史志偉, 昂海松, 等. 俯仰角對鴨式布局飛機搖滾運動的影響與機理分析[J]. 空氣動力學學報, 2017, 35(2): 198-204.

    10.7638/kqdlxxb-2016.0162 Wei D C, Shi Z W, Ang H S, et al. Effects of pitch angle on wing-body rock of canard-configuration aircraft[J]. Acta Aerodynamica Sinica, 2017, 35(2): 198-204.

    V211.7

    A doi: 10.7638/kqdlxxb-2016.0162

    猜你喜歡
    主翼迎角前體
    多段翼低雷諾數(shù)繞流渦-邊界層相互干擾
    航空學報(2023年12期)2023-07-28 10:44:12
    三黃雞
    某型民機低速巡航構型平尾抖振特性風洞試驗研究
    N-末端腦鈉肽前體與糖尿病及糖尿病相關并發(fā)癥呈負相關
    連續(xù)變迎角試驗數(shù)據(jù)自適應分段擬合濾波方法
    探究鴨式布局模型飛機
    航空模型(2017年3期)2017-07-28 21:22:01
    N-端腦鈉肽前體測定在高血壓疾病中的應用研究
    失速保護系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    茶葉香氣前體物研究進展
    茶葉通訊(2014年2期)2014-02-27 07:55:40
    氮末端腦鈉肽前體與急性心肌梗死預后的相關性
    久久久久免费精品人妻一区二区| 91av网一区二区| 日本av手机在线免费观看| 亚洲中文字幕日韩| 亚洲国产日韩欧美精品在线观看| 色视频www国产| 91aial.com中文字幕在线观看| 国产精品永久免费网站| av视频在线观看入口| 人人妻人人澡人人爽人人夜夜 | 精品久久国产蜜桃| 波多野结衣高清无吗| 中文欧美无线码| 汤姆久久久久久久影院中文字幕 | 精品国产一区二区三区久久久樱花 | 少妇裸体淫交视频免费看高清| 国产精品精品国产色婷婷| 日韩三级伦理在线观看| eeuss影院久久| 哪个播放器可以免费观看大片| 自拍偷自拍亚洲精品老妇| 日韩欧美精品免费久久| 一级毛片电影观看 | 精华霜和精华液先用哪个| 色网站视频免费| 亚洲在线自拍视频| 国产乱来视频区| 婷婷六月久久综合丁香| av在线蜜桃| 中文字幕av在线有码专区| 内地一区二区视频在线| 日韩亚洲欧美综合| 国产精品三级大全| 久久99精品国语久久久| 国产毛片a区久久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品国产av成人精品| 亚洲精品乱久久久久久| 国产老妇女一区| av在线老鸭窝| 国产精品一区www在线观看| 成年av动漫网址| 大话2 男鬼变身卡| 国产又黄又爽又无遮挡在线| 国产麻豆成人av免费视频| 免费大片18禁| 精品国产三级普通话版| 国产成人一区二区在线| 亚洲精品乱码久久久v下载方式| 亚洲四区av| 日日干狠狠操夜夜爽| 丰满人妻一区二区三区视频av| 国产免费一级a男人的天堂| 国产精品爽爽va在线观看网站| 丰满人妻一区二区三区视频av| 婷婷色麻豆天堂久久 | 男人舔女人下体高潮全视频| 日本黄大片高清| 国产精品乱码一区二三区的特点| 在线播放国产精品三级| 99热全是精品| 国产黄色小视频在线观看| 国产成人91sexporn| 午夜福利在线观看免费完整高清在| 日本熟妇午夜| 天堂网av新在线| 成人二区视频| 最近手机中文字幕大全| 欧美性猛交╳xxx乱大交人| 91aial.com中文字幕在线观看| 欧美激情国产日韩精品一区| 夜夜爽夜夜爽视频| 十八禁国产超污无遮挡网站| 国产综合懂色| 国产色爽女视频免费观看| 最近最新中文字幕大全电影3| 99久久人妻综合| 老女人水多毛片| 久久人妻av系列| 天堂网av新在线| 亚洲欧美一区二区三区国产| 又粗又硬又长又爽又黄的视频| 天堂中文最新版在线下载 | 99热网站在线观看| h日本视频在线播放| 亚洲精品自拍成人| 精品国产一区二区三区久久久樱花 | 麻豆久久精品国产亚洲av| 国产在线男女| 亚洲人成网站高清观看| 久热久热在线精品观看| 日本欧美国产在线视频| 好男人视频免费观看在线| 国产黄色小视频在线观看| 欧美性猛交╳xxx乱大交人| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩精品成人综合77777| 中文字幕av成人在线电影| 国产极品精品免费视频能看的| 看十八女毛片水多多多| 精品少妇黑人巨大在线播放 | 国产 一区精品| 男人舔奶头视频| 亚洲精品aⅴ在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 99热这里只有是精品50| 久久99热6这里只有精品| 春色校园在线视频观看| 国产免费视频播放在线视频 | 有码 亚洲区| 国产精品一及| 精品一区二区免费观看| 亚洲精品色激情综合| 久久精品夜夜夜夜夜久久蜜豆| 少妇人妻一区二区三区视频| 神马国产精品三级电影在线观看| 色网站视频免费| 网址你懂的国产日韩在线| 国产成人freesex在线| 97超碰精品成人国产| 国产激情偷乱视频一区二区| av线在线观看网站| 亚洲人与动物交配视频| 大香蕉97超碰在线| 毛片女人毛片| 18禁裸乳无遮挡免费网站照片| 国产高清有码在线观看视频| 丰满乱子伦码专区| АⅤ资源中文在线天堂| 欧美精品国产亚洲| 国产高清视频在线观看网站| 免费播放大片免费观看视频在线观看 | 中国美白少妇内射xxxbb| 免费一级毛片在线播放高清视频| 久久精品国产鲁丝片午夜精品| 国产91av在线免费观看| 久久精品国产亚洲av天美| 晚上一个人看的免费电影| 又爽又黄无遮挡网站| 国产精品一二三区在线看| 久久精品国产亚洲网站| 亚洲四区av| 少妇猛男粗大的猛烈进出视频 | videos熟女内射| 99热6这里只有精品| 插阴视频在线观看视频| 国产免费视频播放在线视频 | 久久精品国产自在天天线| 村上凉子中文字幕在线| 国产精品福利在线免费观看| 成人毛片60女人毛片免费| 国产白丝娇喘喷水9色精品| 波多野结衣巨乳人妻| 国产精品电影一区二区三区| 国产探花在线观看一区二区| 国产亚洲av片在线观看秒播厂 | 一级毛片电影观看 | 插阴视频在线观看视频| 国产成人a∨麻豆精品| 美女cb高潮喷水在线观看| 老司机影院成人| 人妻少妇偷人精品九色| 久99久视频精品免费| 日韩亚洲欧美综合| 亚洲五月天丁香| 性插视频无遮挡在线免费观看| 久久精品久久精品一区二区三区| 视频中文字幕在线观看| 精品久久久噜噜| 亚洲精品456在线播放app| 国产亚洲一区二区精品| 免费黄色在线免费观看| 简卡轻食公司| 久99久视频精品免费| 精品人妻熟女av久视频| 日韩 亚洲 欧美在线| 免费在线观看成人毛片| 99热这里只有是精品在线观看| 三级毛片av免费| 麻豆成人午夜福利视频| 男人的好看免费观看在线视频| 午夜福利视频1000在线观看| 一本久久精品| 一夜夜www| 国产精品久久久久久av不卡| 午夜福利在线观看吧| 能在线免费看毛片的网站| 一区二区三区乱码不卡18| 免费观看的影片在线观看| 能在线免费观看的黄片| 青春草国产在线视频| 嘟嘟电影网在线观看| 建设人人有责人人尽责人人享有的 | 啦啦啦观看免费观看视频高清| 中国美白少妇内射xxxbb| 国产白丝娇喘喷水9色精品| 高清午夜精品一区二区三区| 亚洲精品亚洲一区二区| 亚洲成人中文字幕在线播放| 国产精品一及| 在线观看美女被高潮喷水网站| 91在线精品国自产拍蜜月| 午夜免费激情av| 国产麻豆成人av免费视频| 午夜精品国产一区二区电影 | 精品久久久久久久人妻蜜臀av| 日韩成人伦理影院| 波野结衣二区三区在线| 人体艺术视频欧美日本| 国产精品一区二区性色av| 日韩av不卡免费在线播放| 汤姆久久久久久久影院中文字幕 | 天堂av国产一区二区熟女人妻| 亚洲在久久综合| 精品熟女少妇av免费看| 免费看a级黄色片| 国产精品久久久久久精品电影| 免费看美女性在线毛片视频| 久久韩国三级中文字幕| 99在线视频只有这里精品首页| 欧美3d第一页| 亚洲天堂国产精品一区在线| 天堂网av新在线| 国产精品人妻久久久久久| 国产精品国产高清国产av| 亚洲精品aⅴ在线观看| 最新中文字幕久久久久| 亚洲av成人精品一二三区| av又黄又爽大尺度在线免费看 | 久久久国产成人免费| 欧美成人精品欧美一级黄| 国产大屁股一区二区在线视频| 国产午夜精品论理片| 伦理电影大哥的女人| av在线蜜桃| 日韩欧美精品v在线| 亚洲熟妇中文字幕五十中出| 精品免费久久久久久久清纯| 亚洲国产最新在线播放| 国产精品熟女久久久久浪| 人妻系列 视频| 国产乱人偷精品视频| 国产伦精品一区二区三区四那| 久久精品国产亚洲av涩爱| 一个人看的www免费观看视频| 欧美一区二区精品小视频在线| 婷婷色综合大香蕉| 亚洲一级一片aⅴ在线观看| 嫩草影院精品99| 一区二区三区乱码不卡18| a级毛色黄片| 日韩欧美国产在线观看| 日韩在线高清观看一区二区三区| 两个人的视频大全免费| 日韩中字成人| 亚洲精品国产成人久久av| 国产在视频线精品| 久久精品久久久久久久性| 国产精品麻豆人妻色哟哟久久 | 欧美变态另类bdsm刘玥| 国产毛片a区久久久久| 波多野结衣高清无吗| 色综合亚洲欧美另类图片| 少妇猛男粗大的猛烈进出视频 | 在线观看av片永久免费下载| 男女下面进入的视频免费午夜| 亚洲av中文av极速乱| 国产精品一区二区三区四区久久| 日本爱情动作片www.在线观看| 国产精品一区www在线观看| av线在线观看网站| 欧美日韩国产亚洲二区| 三级国产精品欧美在线观看| 2021天堂中文幕一二区在线观| 日韩欧美三级三区| 亚洲aⅴ乱码一区二区在线播放| 国产美女午夜福利| 最近最新中文字幕大全电影3| 国产成人freesex在线| 免费观看a级毛片全部| 狠狠狠狠99中文字幕| 韩国av在线不卡| 黄色欧美视频在线观看| 伦理电影大哥的女人| 亚洲av中文av极速乱| 精品久久久久久久久av| 少妇猛男粗大的猛烈进出视频 | 22中文网久久字幕| 日韩大片免费观看网站 | 最近2019中文字幕mv第一页| 国产淫语在线视频| 日本-黄色视频高清免费观看| av免费观看日本| 日韩强制内射视频| 亚洲av免费高清在线观看| 欧美成人精品欧美一级黄| 麻豆一二三区av精品| 三级经典国产精品| 免费看a级黄色片| 亚洲av福利一区| 好男人视频免费观看在线| 99久国产av精品| 免费av观看视频| 亚洲av日韩在线播放| 亚洲欧美成人综合另类久久久 | 久久久a久久爽久久v久久| 亚洲最大成人中文| 国产极品天堂在线| 高清av免费在线| 毛片女人毛片| 久久久久久久国产电影| 人人妻人人看人人澡| 人妻夜夜爽99麻豆av| 五月玫瑰六月丁香| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 中文字幕av在线有码专区| 精品一区二区三区人妻视频| 色噜噜av男人的天堂激情| 美女xxoo啪啪120秒动态图| 国产淫片久久久久久久久| 国产大屁股一区二区在线视频| 成年版毛片免费区| 2021天堂中文幕一二区在线观| 国国产精品蜜臀av免费| 国产片特级美女逼逼视频| 国产高清三级在线| 国产精品久久久久久久电影| 久久久久久伊人网av| 欧美zozozo另类| 亚洲久久久久久中文字幕| 白带黄色成豆腐渣| 3wmmmm亚洲av在线观看| 69av精品久久久久久| 高清毛片免费看| 一级黄色大片毛片| 免费电影在线观看免费观看| 欧美成人一区二区免费高清观看| 99热这里只有精品一区| 国产成人freesex在线| 最新中文字幕久久久久| 欧美成人a在线观看| 直男gayav资源| 网址你懂的国产日韩在线| 亚洲av成人精品一二三区| 国产精品久久视频播放| 精品久久久久久久久久久久久| 夫妻性生交免费视频一级片| 日韩制服骚丝袜av| 亚洲婷婷狠狠爱综合网| 国产在视频线在精品| 一级黄色大片毛片| 国产精品久久久久久精品电影| 亚洲五月天丁香| 蜜桃久久精品国产亚洲av| 免费av不卡在线播放| 国产成人免费观看mmmm| 一区二区三区乱码不卡18| 国产综合懂色| 精品熟女少妇av免费看| 国产男人的电影天堂91| 国产黄色视频一区二区在线观看 | 在线观看一区二区三区| 日韩av在线大香蕉| 亚洲久久久久久中文字幕| 国产精品.久久久| 国产不卡一卡二| 久久韩国三级中文字幕| 蜜桃久久精品国产亚洲av| h日本视频在线播放| 久久久久久久久大av| 日韩高清综合在线| 国产一区有黄有色的免费视频 | 男人狂女人下面高潮的视频| 国产伦在线观看视频一区| 国产真实伦视频高清在线观看| 久久久精品欧美日韩精品| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 91aial.com中文字幕在线观看| 国产黄色小视频在线观看| 女的被弄到高潮叫床怎么办| 亚洲在久久综合| 色噜噜av男人的天堂激情| 日韩av在线大香蕉| 久久久久久久久久久丰满| 熟女电影av网| 亚洲国产欧洲综合997久久,| 久久精品夜色国产| 欧美+日韩+精品| 我的女老师完整版在线观看| 久久久久性生活片| 国产免费男女视频| 少妇被粗大猛烈的视频| 免费搜索国产男女视频| 级片在线观看| 亚洲av日韩在线播放| 成人特级av手机在线观看| 国产乱来视频区| 成人高潮视频无遮挡免费网站| 热99在线观看视频| 毛片女人毛片| 赤兔流量卡办理| 看非洲黑人一级黄片| 国语对白做爰xxxⅹ性视频网站| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 国产黄片视频在线免费观看| 欧美成人午夜免费资源| 亚洲第一区二区三区不卡| 综合色av麻豆| 99热这里只有精品一区| 免费观看在线日韩| 色综合站精品国产| 99久国产av精品| 你懂的网址亚洲精品在线观看 | 高清在线视频一区二区三区 | 一级黄色大片毛片| 99久久精品热视频| 偷拍熟女少妇极品色| 蜜桃亚洲精品一区二区三区| 最近中文字幕2019免费版| 亚洲内射少妇av| 成人二区视频| 99视频精品全部免费 在线| 精品久久国产蜜桃| 亚洲精品456在线播放app| 人体艺术视频欧美日本| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 久久草成人影院| av.在线天堂| 狂野欧美激情性xxxx在线观看| 午夜精品国产一区二区电影 | 日韩三级伦理在线观看| 久久人人爽人人片av| 欧美最新免费一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲欧美精品综合久久99| 日本黄色视频三级网站网址| 欧美性猛交╳xxx乱大交人| 久久精品国产鲁丝片午夜精品| 22中文网久久字幕| 69av精品久久久久久| 亚洲aⅴ乱码一区二区在线播放| 天堂网av新在线| 91aial.com中文字幕在线观看| 久久国内精品自在自线图片| 色5月婷婷丁香| 99久久中文字幕三级久久日本| a级毛色黄片| av.在线天堂| 午夜亚洲福利在线播放| 国产高清三级在线| 亚洲成人精品中文字幕电影| 男人和女人高潮做爰伦理| 精品少妇黑人巨大在线播放 | 尤物成人国产欧美一区二区三区| 久久欧美精品欧美久久欧美| 一边摸一边抽搐一进一小说| 国产伦一二天堂av在线观看| 爱豆传媒免费全集在线观看| 成人鲁丝片一二三区免费| 成年女人看的毛片在线观看| 一个人免费在线观看电影| 欧美一级a爱片免费观看看| 日本午夜av视频| 麻豆成人午夜福利视频| 黄色配什么色好看| 亚洲国产欧美人成| 久久久久久久久久久丰满| 国产av在哪里看| 一本一本综合久久| 国产毛片a区久久久久| 国国产精品蜜臀av免费| 级片在线观看| 久久国产乱子免费精品| 99在线视频只有这里精品首页| 色综合色国产| 国产亚洲最大av| 乱码一卡2卡4卡精品| 精品人妻视频免费看| 亚洲欧美中文字幕日韩二区| 女人久久www免费人成看片 | 联通29元200g的流量卡| 亚洲怡红院男人天堂| 欧美3d第一页| 欧美精品一区二区大全| 亚洲精品影视一区二区三区av| 在线播放无遮挡| 亚洲18禁久久av| 亚洲人成网站在线播| 婷婷色麻豆天堂久久 | 一个人看的www免费观看视频| 久久久久久大精品| 丝袜喷水一区| 看非洲黑人一级黄片| 国产中年淑女户外野战色| 在线观看66精品国产| 69人妻影院| 亚洲精华国产精华液的使用体验| 美女黄网站色视频| 老司机福利观看| 日韩一本色道免费dvd| 亚洲国产最新在线播放| 久久久久国产网址| 国产一区二区亚洲精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 日本午夜av视频| 青春草国产在线视频| 永久免费av网站大全| 亚洲国产最新在线播放| 天天一区二区日本电影三级| 综合色av麻豆| 一级毛片aaaaaa免费看小| 国产精品久久久久久精品电影| 国产一级毛片在线| 亚洲真实伦在线观看| 乱系列少妇在线播放| 日本-黄色视频高清免费观看| 毛片一级片免费看久久久久| 亚洲欧美日韩东京热| 乱人视频在线观看| 秋霞伦理黄片| 国产精品一区二区三区四区久久| 亚洲国产精品久久男人天堂| 深夜a级毛片| 99久国产av精品国产电影| 少妇高潮的动态图| 舔av片在线| 免费观看精品视频网站| 久久婷婷人人爽人人干人人爱| 国产69精品久久久久777片| 日韩欧美 国产精品| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩东京热| 中文字幕熟女人妻在线| 国产成人a区在线观看| 午夜a级毛片| 婷婷色av中文字幕| 午夜日本视频在线| 男人和女人高潮做爰伦理| .国产精品久久| 久久精品久久精品一区二区三区| kizo精华| 久久精品综合一区二区三区| 国产免费又黄又爽又色| 国产成人91sexporn| 午夜福利成人在线免费观看| 亚洲欧美精品专区久久| 最近的中文字幕免费完整| 大话2 男鬼变身卡| 亚洲精品国产av成人精品| 亚洲成色77777| 亚洲欧美日韩卡通动漫| 高清视频免费观看一区二区 | 全区人妻精品视频| 久久综合国产亚洲精品| 欧美高清性xxxxhd video| 亚洲精品亚洲一区二区| 午夜爱爱视频在线播放| 国产午夜福利久久久久久| 性色avwww在线观看| 波野结衣二区三区在线| 日韩av在线大香蕉| 老司机影院毛片| 六月丁香七月| 一个人看的www免费观看视频| 91在线精品国自产拍蜜月| 国产成人一区二区在线| av在线天堂中文字幕| 国产成人午夜福利电影在线观看| 精品久久久久久久久av| 成人综合一区亚洲| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久久丰满| 草草在线视频免费看| 国产91av在线免费观看| 3wmmmm亚洲av在线观看| 如何舔出高潮| 日韩中字成人| 成人亚洲精品av一区二区| 五月伊人婷婷丁香| 超碰97精品在线观看| 国产午夜精品久久久久久一区二区三区| 麻豆乱淫一区二区| 又粗又爽又猛毛片免费看| 最近2019中文字幕mv第一页| 午夜福利高清视频| 乱系列少妇在线播放| 欧美一区二区精品小视频在线| 在线免费十八禁| 日韩欧美 国产精品| 午夜爱爱视频在线播放| 午夜福利高清视频| 精品人妻熟女av久视频| 日本免费一区二区三区高清不卡| 人妻少妇偷人精品九色| 天天躁夜夜躁狠狠久久av| 国产成人aa在线观看| 91午夜精品亚洲一区二区三区| 日韩 亚洲 欧美在线| 久久久亚洲精品成人影院| 久久精品夜色国产| 汤姆久久久久久久影院中文字幕 | 天堂√8在线中文| 久久精品熟女亚洲av麻豆精品 | 亚洲国产高清在线一区二区三| 午夜免费激情av| 永久网站在线| 国产精品一区www在线观看| 久久热精品热| 精品无人区乱码1区二区|