• <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
    氮末端腦鈉肽前體與急性心肌梗死預后的相關性
    精品高清国产在线一区| 美女免费视频网站| 欧美性猛交╳xxx乱大交人| 成人手机av| 精品久久蜜臀av无| 久久 成人 亚洲| 一边摸一边做爽爽视频免费| 久久精品夜夜夜夜夜久久蜜豆 | 日韩大码丰满熟妇| 在线播放国产精品三级| 一边摸一边做爽爽视频免费| 日韩欧美三级三区| 好男人在线观看高清免费视频| 亚洲欧美精品综合久久99| netflix在线观看网站| 亚洲一区二区三区色噜噜| 日韩成人在线观看一区二区三区| 国产亚洲精品久久久久5区| 99热6这里只有精品| 成人一区二区视频在线观看| 久久久久久九九精品二区国产 | 亚洲国产欧美网| 亚洲aⅴ乱码一区二区在线播放 | 免费高清视频大片| 午夜精品在线福利| 一级毛片女人18水好多| 久热爱精品视频在线9| 99re在线观看精品视频| 夜夜爽天天搞| 久久伊人香网站| 欧美不卡视频在线免费观看 | 亚洲人成网站高清观看| 色哟哟哟哟哟哟| 久久中文看片网| 午夜福利在线观看吧| 亚洲精品国产一区二区精华液| 欧美绝顶高潮抽搐喷水| av国产免费在线观看| 欧美日韩乱码在线| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 午夜免费观看网址| 久久九九热精品免费| 午夜久久久久精精品| 婷婷精品国产亚洲av| 亚洲成人中文字幕在线播放| 色精品久久人妻99蜜桃| 亚洲在线自拍视频| 国产爱豆传媒在线观看 | 亚洲va日本ⅴa欧美va伊人久久| av天堂在线播放| 国产亚洲欧美在线一区二区| 午夜精品久久久久久毛片777| 午夜精品一区二区三区免费看| 久久国产精品影院| 男女做爰动态图高潮gif福利片| 成人午夜高清在线视频| 校园春色视频在线观看| 国产精品一区二区精品视频观看| 色综合亚洲欧美另类图片| 国产精品久久久久久人妻精品电影| 亚洲第一欧美日韩一区二区三区| 国语自产精品视频在线第100页| 五月伊人婷婷丁香| 亚洲aⅴ乱码一区二区在线播放 | 亚洲性夜色夜夜综合| 无人区码免费观看不卡| 特级一级黄色大片| 97人妻精品一区二区三区麻豆| 一个人免费在线观看电影 | 亚洲一区中文字幕在线| 1024香蕉在线观看| 嫁个100分男人电影在线观看| 国产高清videossex| 一级黄色大片毛片| 九色成人免费人妻av| 午夜精品久久久久久毛片777| 一级毛片精品| 波多野结衣高清无吗| 亚洲人成77777在线视频| 精品无人区乱码1区二区| 99久久99久久久精品蜜桃| 日本五十路高清| 一个人免费在线观看的高清视频| 色综合欧美亚洲国产小说| 成人高潮视频无遮挡免费网站| 国产亚洲精品第一综合不卡| 免费人成视频x8x8入口观看| 亚洲一卡2卡3卡4卡5卡精品中文| 高潮久久久久久久久久久不卡| 手机成人av网站| 日韩欧美 国产精品| 亚洲av成人av| 欧美色欧美亚洲另类二区| 91九色精品人成在线观看| 亚洲欧美日韩无卡精品| 午夜福利在线观看吧| 九色国产91popny在线| 动漫黄色视频在线观看| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 777久久人妻少妇嫩草av网站| 国产av一区二区精品久久| 1024视频免费在线观看| 天堂av国产一区二区熟女人妻 | 亚洲色图 男人天堂 中文字幕| 啦啦啦观看免费观看视频高清| 久久久久久人人人人人| 欧美成人性av电影在线观看| 精品第一国产精品| 精品人妻1区二区| 久久人妻福利社区极品人妻图片| 可以免费在线观看a视频的电影网站| 婷婷亚洲欧美| 长腿黑丝高跟| 99热6这里只有精品| 亚洲专区字幕在线| 色噜噜av男人的天堂激情| 狠狠狠狠99中文字幕| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 精品国产乱码久久久久久男人| 精品午夜福利视频在线观看一区| 夜夜看夜夜爽夜夜摸| 哪里可以看免费的av片| 精品一区二区三区四区五区乱码| 成人高潮视频无遮挡免费网站| 婷婷六月久久综合丁香| 一级作爱视频免费观看| 日本熟妇午夜| 黄片小视频在线播放| 午夜精品久久久久久毛片777| 亚洲av成人一区二区三| 亚洲片人在线观看| www.www免费av| 免费在线观看日本一区| av国产免费在线观看| www.www免费av| 国产精品电影一区二区三区| 五月伊人婷婷丁香| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品亚洲精品国产色婷小说| 1024香蕉在线观看| 精品一区二区三区四区五区乱码| 麻豆成人av在线观看| 五月玫瑰六月丁香| 久久这里只有精品19| 又黄又粗又硬又大视频| 男女床上黄色一级片免费看| 国产v大片淫在线免费观看| 亚洲专区国产一区二区| 欧美性猛交黑人性爽| 精品久久蜜臀av无| 少妇粗大呻吟视频| 国产亚洲精品久久久久5区| 在线播放国产精品三级| 2021天堂中文幕一二区在线观| 夜夜爽天天搞| 欧美高清成人免费视频www| 欧美黄色片欧美黄色片| 又爽又黄无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 天堂动漫精品| 精品日产1卡2卡| 国产视频一区二区在线看| 天堂√8在线中文| 亚洲一区高清亚洲精品| 国产午夜福利久久久久久| 天天添夜夜摸| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| 精品久久久久久久末码| 真人做人爱边吃奶动态| 又粗又爽又猛毛片免费看| 两个人免费观看高清视频| 午夜免费激情av| 麻豆一二三区av精品| 五月伊人婷婷丁香| 无人区码免费观看不卡| 一二三四在线观看免费中文在| 老汉色av国产亚洲站长工具| xxxwww97欧美| 亚洲成人免费电影在线观看| 极品教师在线免费播放| 精品一区二区三区av网在线观看| 国产野战对白在线观看| 久久午夜亚洲精品久久| 国内精品一区二区在线观看| 国产精品日韩av在线免费观看| 特大巨黑吊av在线直播| 91字幕亚洲| 亚洲一码二码三码区别大吗| 99久久99久久久精品蜜桃| 国产蜜桃级精品一区二区三区| 久久香蕉激情| 国产视频一区二区在线看| or卡值多少钱| 国产在线精品亚洲第一网站| 99re在线观看精品视频| 亚洲成人久久爱视频| 欧美日韩瑟瑟在线播放| 亚洲精品美女久久av网站| 久久久久久亚洲精品国产蜜桃av| 国产一级毛片七仙女欲春2| 久久人妻福利社区极品人妻图片| 成人一区二区视频在线观看| 欧美午夜高清在线| 免费在线观看成人毛片| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 在线国产一区二区在线| 欧美黑人欧美精品刺激| 特大巨黑吊av在线直播| 神马国产精品三级电影在线观看 | 午夜视频精品福利| 中国美女看黄片| 日韩欧美在线乱码| 国产精品免费一区二区三区在线| 亚洲avbb在线观看| 欧美日本亚洲视频在线播放| 成人手机av| 日本熟妇午夜| 亚洲精品国产精品久久久不卡| 欧美成人免费av一区二区三区| 悠悠久久av| 两个人视频免费观看高清| 国内精品久久久久精免费| 精华霜和精华液先用哪个| 麻豆国产av国片精品| 又黄又爽又免费观看的视频| 亚洲无线在线观看| 国产探花在线观看一区二区| 1024香蕉在线观看| 可以免费在线观看a视频的电影网站| 国产一区二区三区视频了| 亚洲精品中文字幕一二三四区| 啦啦啦韩国在线观看视频| 中文字幕人成人乱码亚洲影| 国产视频内射| 亚洲精品中文字幕在线视频| 99国产综合亚洲精品| 精品无人区乱码1区二区| 老司机在亚洲福利影院| 久久香蕉国产精品| 在线看三级毛片| 黄色成人免费大全| 亚洲欧美一区二区三区黑人| 搡老熟女国产l中国老女人| 亚洲国产欧洲综合997久久,| 国产片内射在线| 亚洲真实伦在线观看| 午夜精品在线福利| 中文亚洲av片在线观看爽| 黄频高清免费视频| 亚洲国产日韩欧美精品在线观看 | avwww免费| 色哟哟哟哟哟哟| 俄罗斯特黄特色一大片| 99精品在免费线老司机午夜| 免费电影在线观看免费观看| 国产久久久一区二区三区| 人人妻,人人澡人人爽秒播| 黄片大片在线免费观看| 欧美精品亚洲一区二区| 一夜夜www| 亚洲精品av麻豆狂野| 免费av毛片视频| 久久草成人影院| 久久久久久九九精品二区国产 | 国产成人影院久久av| 中文字幕av在线有码专区| 一区福利在线观看| www国产在线视频色| 91字幕亚洲| 免费看日本二区| 又紧又爽又黄一区二区| 亚洲国产欧美一区二区综合| 国产伦在线观看视频一区| 国产精品免费一区二区三区在线| 在线播放国产精品三级| 久久精品aⅴ一区二区三区四区| 精品一区二区三区四区五区乱码| 91字幕亚洲| 亚洲成av人片免费观看| 在线观看免费视频日本深夜| 男女视频在线观看网站免费 | av超薄肉色丝袜交足视频| 国产亚洲精品av在线| 波多野结衣高清作品| 日本三级黄在线观看| 少妇人妻一区二区三区视频| 亚洲欧美激情综合另类| 国产高清视频在线播放一区| 久久久国产成人免费| 女同久久另类99精品国产91| 99久久无色码亚洲精品果冻| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕久久专区| 五月玫瑰六月丁香| 精品第一国产精品| 91麻豆精品激情在线观看国产| 亚洲第一欧美日韩一区二区三区| 在线视频色国产色| e午夜精品久久久久久久| ponron亚洲| 高清毛片免费观看视频网站| 又紧又爽又黄一区二区| 亚洲精品久久国产高清桃花| 变态另类成人亚洲欧美熟女| 夜夜爽天天搞| 手机成人av网站| 国产伦人伦偷精品视频| 欧美日韩国产亚洲二区| 狠狠狠狠99中文字幕| 成人精品一区二区免费| videosex国产| 看黄色毛片网站| 亚洲国产欧洲综合997久久,| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| 中文在线观看免费www的网站 | 中文字幕久久专区| 免费高清视频大片| 国产精品电影一区二区三区| 亚洲精品国产精品久久久不卡| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 天堂动漫精品| 亚洲 欧美一区二区三区| 国产精品免费一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| av在线播放免费不卡| 欧美丝袜亚洲另类 | 最近视频中文字幕2019在线8| 国内精品久久久久精免费| 在线免费观看的www视频| 国产高清videossex| 久久精品aⅴ一区二区三区四区| 成熟少妇高潮喷水视频| 亚洲精品久久成人aⅴ小说| 国产单亲对白刺激| av国产免费在线观看| 少妇的丰满在线观看| 99国产极品粉嫩在线观看| 操出白浆在线播放| 国产三级中文精品| 桃色一区二区三区在线观看| 国产野战对白在线观看| 久久草成人影院| 嫩草影视91久久| 99精品久久久久人妻精品| 不卡一级毛片| 久久香蕉激情| 嫩草影视91久久| 黄色片一级片一级黄色片| 国内少妇人妻偷人精品xxx网站 | 国产亚洲精品综合一区在线观看 | 亚洲一区中文字幕在线| 国产精品 国内视频| 成人国产综合亚洲| 中亚洲国语对白在线视频| 手机成人av网站| 日韩av在线大香蕉| 老司机在亚洲福利影院| www.www免费av| 精品不卡国产一区二区三区| 国产亚洲av嫩草精品影院| 精品国产亚洲在线| 成人av在线播放网站| 12—13女人毛片做爰片一| 午夜精品久久久久久毛片777| 99久久久亚洲精品蜜臀av| 99精品久久久久人妻精品| 一进一出抽搐gif免费好疼| 好男人在线观看高清免费视频| 亚洲一区二区三区色噜噜| 亚洲欧美日韩高清专用| 亚洲 国产 在线| 中文字幕精品亚洲无线码一区| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 久久久久亚洲av毛片大全| x7x7x7水蜜桃| 午夜影院日韩av| 两人在一起打扑克的视频| 男女午夜视频在线观看| 精品日产1卡2卡| 国内揄拍国产精品人妻在线| 欧美国产日韩亚洲一区| 丰满人妻一区二区三区视频av | 国产av又大| 999久久久精品免费观看国产| 欧美乱妇无乱码| 亚洲中文字幕一区二区三区有码在线看 | 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 听说在线观看完整版免费高清| 韩国av一区二区三区四区| 久久久久久久久久黄片| 久久伊人香网站| 草草在线视频免费看| 精品日产1卡2卡| 亚洲av成人不卡在线观看播放网| 欧美色视频一区免费| 日韩欧美 国产精品| 99久久国产精品久久久| 国产野战对白在线观看| 亚洲男人的天堂狠狠| 男女做爰动态图高潮gif福利片| 国产片内射在线| 亚洲av美国av| 一级毛片精品| 制服丝袜大香蕉在线| 亚洲最大成人中文| 日韩高清综合在线| 12—13女人毛片做爰片一| 国产成人一区二区三区免费视频网站| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 五月玫瑰六月丁香| 国产亚洲精品综合一区在线观看 | 国产亚洲av高清不卡| 国产精华一区二区三区| 成人国产一区最新在线观看| e午夜精品久久久久久久| av国产免费在线观看| 久久伊人香网站| 岛国视频午夜一区免费看| 成人欧美大片| 午夜亚洲福利在线播放| 国产精品1区2区在线观看.| 18禁裸乳无遮挡免费网站照片| 一级毛片女人18水好多| 国产一区二区激情短视频| 亚洲av片天天在线观看| 一夜夜www| 手机成人av网站| 俺也久久电影网| 国产高清视频在线播放一区| 国产69精品久久久久777片 | 精品久久久久久久末码| 悠悠久久av| 久久久久性生活片| 国产不卡一卡二| 国产人伦9x9x在线观看| 99在线人妻在线中文字幕| netflix在线观看网站| 亚洲精品中文字幕一二三四区| 国产精品免费视频内射| av福利片在线| 国产精品久久久av美女十八| 久久午夜综合久久蜜桃| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 岛国在线免费视频观看| 精品国产超薄肉色丝袜足j| 狂野欧美白嫩少妇大欣赏| 欧美av亚洲av综合av国产av| 色尼玛亚洲综合影院| 深夜精品福利| 欧美高清成人免费视频www| 岛国在线免费视频观看| 精品一区二区三区视频在线观看免费| 非洲黑人性xxxx精品又粗又长| www日本黄色视频网| 观看免费一级毛片| 国产精品爽爽va在线观看网站| 超碰成人久久| 又大又爽又粗| 欧美+亚洲+日韩+国产| 国产v大片淫在线免费观看| 91麻豆精品激情在线观看国产| 国产亚洲欧美98| 国产av一区二区精品久久| 中文字幕久久专区| 性色av乱码一区二区三区2| 久热爱精品视频在线9| 欧美人与性动交α欧美精品济南到| 国产三级在线视频| 欧美国产日韩亚洲一区| 国产三级黄色录像| 岛国在线观看网站| 好男人在线观看高清免费视频| 国产三级在线视频| 欧美大码av| 91成年电影在线观看| 亚洲最大成人中文| 一本一本综合久久| 欧美av亚洲av综合av国产av| 久久热在线av| 久久久精品欧美日韩精品| 国产乱人伦免费视频| 99在线视频只有这里精品首页| 人人妻人人澡欧美一区二区| 成人18禁在线播放| 国产精品日韩av在线免费观看| 亚洲最大成人中文| 1024手机看黄色片| 18禁黄网站禁片免费观看直播| 在线播放国产精品三级| 九九热线精品视视频播放| 白带黄色成豆腐渣| 在线观看日韩欧美| 男人舔女人下体高潮全视频| 国产成人av激情在线播放| 我要搜黄色片| 亚洲成人免费电影在线观看| 97超级碰碰碰精品色视频在线观看| 中文字幕人成人乱码亚洲影| 午夜福利18| 欧美色欧美亚洲另类二区| www.自偷自拍.com| 精品无人区乱码1区二区| 97碰自拍视频| 亚洲欧美激情综合另类| 国产男靠女视频免费网站| 国产日本99.免费观看| 国产av在哪里看| 国产区一区二久久| 日韩欧美精品v在线| 特大巨黑吊av在线直播| 他把我摸到了高潮在线观看| 正在播放国产对白刺激| 亚洲18禁久久av| 国产1区2区3区精品| 精品久久久久久久毛片微露脸| 一本综合久久免费| 18禁美女被吸乳视频| 国产黄片美女视频| 18禁裸乳无遮挡免费网站照片| 妹子高潮喷水视频| 在线观看日韩欧美| 女人高潮潮喷娇喘18禁视频| 中出人妻视频一区二区| 一a级毛片在线观看| 日韩欧美国产在线观看| 亚洲熟妇中文字幕五十中出| 免费一级毛片在线播放高清视频| 露出奶头的视频| 日本熟妇午夜| cao死你这个sao货| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月婷婷丁香| 精品电影一区二区在线| 久久久精品国产亚洲av高清涩受| 亚洲成人精品中文字幕电影| 午夜福利成人在线免费观看| 国产精品一及| 国产精品电影一区二区三区| 久久久精品欧美日韩精品| 欧美日韩国产亚洲二区| 岛国在线免费视频观看| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 国产精品1区2区在线观看.| 十八禁人妻一区二区| 亚洲中文字幕一区二区三区有码在线看 | 嫩草影院精品99| 成人高潮视频无遮挡免费网站| 国产精品 国内视频| 一区二区三区激情视频| 国产真实乱freesex| 天天躁夜夜躁狠狠躁躁| 国内精品久久久久久久电影| 成人三级黄色视频| 国产精品永久免费网站| 久久久久久久久久黄片| 丁香六月欧美| 免费看美女性在线毛片视频| 午夜激情av网站| 一夜夜www| 一边摸一边做爽爽视频免费| 首页视频小说图片口味搜索| www.999成人在线观看| 欧美不卡视频在线免费观看 | 久久精品aⅴ一区二区三区四区| 日韩欧美在线乱码| 91国产中文字幕| 日韩 欧美 亚洲 中文字幕| 免费一级毛片在线播放高清视频| 欧美成人一区二区免费高清观看 | 日本a在线网址| 99久久久亚洲精品蜜臀av| 日本三级黄在线观看| 极品教师在线免费播放| 人妻夜夜爽99麻豆av| 美女扒开内裤让男人捅视频| 亚洲狠狠婷婷综合久久图片| 国产99白浆流出| 亚洲国产中文字幕在线视频| 久久久精品欧美日韩精品| 国产99白浆流出| 国产激情久久老熟女| 听说在线观看完整版免费高清| 身体一侧抽搐| 国内久久婷婷六月综合欲色啪| 国产精品av久久久久免费| 国产99白浆流出| 日韩欧美在线二视频| 叶爱在线成人免费视频播放| 欧美日韩黄片免| 亚洲国产中文字幕在线视频| 黄色视频,在线免费观看| 久久精品91蜜桃| 亚洲男人的天堂狠狠| 久久中文字幕人妻熟女| 日韩欧美 国产精品| 一本精品99久久精品77| av视频在线观看入口| 国产三级中文精品| 国产成人av激情在线播放| 中亚洲国语对白在线视频| 国产精品久久久av美女十八|