• <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
    氮末端腦鈉肽前體與急性心肌梗死預后的相關性
    97碰自拍视频| 成年女人永久免费观看视频| 中文资源天堂在线| 非洲黑人性xxxx精品又粗又长| 亚洲av.av天堂| 国产 一区 欧美 日韩| 亚洲国产欧美人成| 亚洲一区二区三区色噜噜| 国产久久久一区二区三区| 高清毛片免费观看视频网站| 亚洲综合色惰| 九九在线视频观看精品| 男女那种视频在线观看| 高清日韩中文字幕在线| 午夜福利在线观看免费完整高清在 | 啦啦啦观看免费观看视频高清| 欧美在线一区亚洲| 亚洲专区国产一区二区| 性欧美人与动物交配| 一个人免费在线观看的高清视频| 欧美精品啪啪一区二区三区| 啪啪无遮挡十八禁网站| 在线观看免费视频日本深夜| 国产亚洲欧美98| 欧美日韩中文字幕国产精品一区二区三区| 99久久久亚洲精品蜜臀av| 夜夜躁狠狠躁天天躁| 成年人黄色毛片网站| 久久人人爽人人爽人人片va | 色哟哟哟哟哟哟| 日本成人三级电影网站| 在线播放无遮挡| 日韩 亚洲 欧美在线| 亚洲av一区综合| 午夜免费成人在线视频| 真人做人爱边吃奶动态| 国产麻豆成人av免费视频| 久久6这里有精品| 高清日韩中文字幕在线| 精品一区二区三区视频在线观看免费| 男女那种视频在线观看| 很黄的视频免费| 91久久精品国产一区二区成人| 精品乱码久久久久久99久播| 在线观看av片永久免费下载| 伊人久久精品亚洲午夜| 国产精品一区二区免费欧美| 欧美色视频一区免费| 我要看日韩黄色一级片| 精品国产三级普通话版| 亚洲国产高清在线一区二区三| 欧美最黄视频在线播放免费| 免费在线观看成人毛片| 亚洲精品影视一区二区三区av| 丰满人妻一区二区三区视频av| 日韩欧美 国产精品| 欧美日韩国产亚洲二区| 别揉我奶头 嗯啊视频| 毛片女人毛片| 欧美性猛交黑人性爽| 久久精品久久久久久噜噜老黄 | 少妇丰满av| 成人av一区二区三区在线看| 亚洲三级黄色毛片| 久久精品国产99精品国产亚洲性色| 国产亚洲精品久久久com| 久久久精品欧美日韩精品| 丰满的人妻完整版| 波野结衣二区三区在线| 91麻豆av在线| 久久久久九九精品影院| 色噜噜av男人的天堂激情| 亚洲午夜理论影院| 一本一本综合久久| 亚洲激情在线av| 亚洲avbb在线观看| 久久精品国产亚洲av香蕉五月| 国产老妇女一区| netflix在线观看网站| 老女人水多毛片| 亚洲一区高清亚洲精品| 亚洲国产精品成人综合色| 亚洲狠狠婷婷综合久久图片| 美女大奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 一进一出抽搐gif免费好疼| 在线国产一区二区在线| 精品99又大又爽又粗少妇毛片 | 男女下面进入的视频免费午夜| 国产午夜福利久久久久久| 18+在线观看网站| 欧美高清性xxxxhd video| 麻豆成人午夜福利视频| 亚洲av第一区精品v没综合| 精品久久久久久久久亚洲 | 国产高清三级在线| 亚洲成a人片在线一区二区| 一边摸一边抽搐一进一小说| 九九热线精品视视频播放| 男女做爰动态图高潮gif福利片| 国产精品综合久久久久久久免费| 搞女人的毛片| 亚洲中文日韩欧美视频| 国产又黄又爽又无遮挡在线| 亚洲精品456在线播放app | 亚洲av一区综合| 欧美黄色片欧美黄色片| 少妇被粗大猛烈的视频| 亚洲av成人不卡在线观看播放网| 成人高潮视频无遮挡免费网站| 欧美最黄视频在线播放免费| 国产久久久一区二区三区| 色视频www国产| 久久久久九九精品影院| 日日干狠狠操夜夜爽| 国产一区二区亚洲精品在线观看| 搡女人真爽免费视频火全软件 | 神马国产精品三级电影在线观看| 小说图片视频综合网站| 日本熟妇午夜| 午夜免费男女啪啪视频观看 | 亚洲精品在线美女| 色5月婷婷丁香| 亚洲在线观看片| 亚洲久久久久久中文字幕| 18禁裸乳无遮挡免费网站照片| 久久6这里有精品| 国产精品精品国产色婷婷| 久久久久久国产a免费观看| www.999成人在线观看| 搡女人真爽免费视频火全软件 | 色综合婷婷激情| 99热精品在线国产| 精品人妻一区二区三区麻豆 | 欧美zozozo另类| 哪里可以看免费的av片| 亚洲美女视频黄频| 亚洲美女黄片视频| 男人舔奶头视频| 欧美三级亚洲精品| 国产精品98久久久久久宅男小说| 亚洲一区二区三区不卡视频| 美女被艹到高潮喷水动态| 国产成人a区在线观看| 精品午夜福利视频在线观看一区| 国产麻豆成人av免费视频| 亚洲成av人片免费观看| 日日干狠狠操夜夜爽| 国产一区二区亚洲精品在线观看| 精品不卡国产一区二区三区| 91麻豆av在线| 久久国产乱子免费精品| 国产爱豆传媒在线观看| 黄色丝袜av网址大全| 成年女人毛片免费观看观看9| x7x7x7水蜜桃| 国产日本99.免费观看| 亚洲人成电影免费在线| 毛片女人毛片| 午夜精品在线福利| 男人和女人高潮做爰伦理| 一级毛片久久久久久久久女| 日韩欧美国产一区二区入口| 亚洲专区中文字幕在线| 91字幕亚洲| 久久热精品热| 国产精品综合久久久久久久免费| 人妻夜夜爽99麻豆av| 麻豆久久精品国产亚洲av| 精品一区二区三区人妻视频| 精品一区二区三区人妻视频| 午夜福利在线观看免费完整高清在 | 国产亚洲欧美在线一区二区| 一二三四社区在线视频社区8| 国产私拍福利视频在线观看| 久久亚洲精品不卡| 亚洲av美国av| 久久久久久久午夜电影| 国产一区二区亚洲精品在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲精品粉嫩美女一区| 毛片一级片免费看久久久久 | 在线观看66精品国产| 亚洲专区国产一区二区| 国产老妇女一区| 亚洲成人免费电影在线观看| 国产免费一级a男人的天堂| 成人毛片a级毛片在线播放| 色哟哟哟哟哟哟| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 精品人妻偷拍中文字幕| 日韩精品青青久久久久久| 直男gayav资源| 午夜精品久久久久久毛片777| 精品欧美国产一区二区三| 午夜精品久久久久久毛片777| 亚洲人成伊人成综合网2020| 男人狂女人下面高潮的视频| 欧美另类亚洲清纯唯美| 丰满人妻熟妇乱又伦精品不卡| 日本在线视频免费播放| 国产综合懂色| 成人美女网站在线观看视频| 可以在线观看毛片的网站| 国产精品人妻久久久久久| 特级一级黄色大片| netflix在线观看网站| 色综合婷婷激情| 国产伦一二天堂av在线观看| 国产精品不卡视频一区二区 | 91久久精品电影网| 亚洲 欧美 日韩 在线 免费| 国产黄片美女视频| h日本视频在线播放| 欧美绝顶高潮抽搐喷水| 欧美激情国产日韩精品一区| 午夜免费男女啪啪视频观看 | 日韩欧美 国产精品| 亚洲专区国产一区二区| 少妇熟女aⅴ在线视频| 免费在线观看日本一区| 国产精品永久免费网站| 老司机午夜福利在线观看视频| 午夜精品久久久久久毛片777| 制服丝袜大香蕉在线| 国产精品,欧美在线| 3wmmmm亚洲av在线观看| 又紧又爽又黄一区二区| 久久久久免费精品人妻一区二区| 高清毛片免费观看视频网站| 一本一本综合久久| 成年人黄色毛片网站| 禁无遮挡网站| 桃红色精品国产亚洲av| 国产人妻一区二区三区在| 亚洲人成电影免费在线| 夜夜夜夜夜久久久久| 国产视频内射| 欧美日韩亚洲国产一区二区在线观看| 在线播放国产精品三级| 在线观看午夜福利视频| 免费搜索国产男女视频| 国产亚洲av嫩草精品影院| 国产成人aa在线观看| 国产日本99.免费观看| 搡女人真爽免费视频火全软件 | 搡老熟女国产l中国老女人| 窝窝影院91人妻| 国产欧美日韩精品一区二区| 好看av亚洲va欧美ⅴa在| 欧美激情国产日韩精品一区| 老司机深夜福利视频在线观看| 久久久久久久久中文| 脱女人内裤的视频| 久久久精品大字幕| 成人永久免费在线观看视频| 在线观看66精品国产| 成人国产一区最新在线观看| 色噜噜av男人的天堂激情| 天堂动漫精品| 乱人视频在线观看| 天天躁日日操中文字幕| 午夜精品在线福利| 中出人妻视频一区二区| 国产探花在线观看一区二区| 成人永久免费在线观看视频| 亚洲 欧美 日韩 在线 免费| 久久国产乱子伦精品免费另类| 波多野结衣高清作品| 久久久久久久精品吃奶| 国产精品伦人一区二区| 亚洲av成人精品一区久久| 国产探花极品一区二区| 每晚都被弄得嗷嗷叫到高潮| 亚洲美女视频黄频| 深爱激情五月婷婷| 男人和女人高潮做爰伦理| 欧美日韩亚洲国产一区二区在线观看| 国产精品野战在线观看| 国产精品98久久久久久宅男小说| 午夜免费成人在线视频| 国产成年人精品一区二区| 一区二区三区四区激情视频 | 51午夜福利影视在线观看| 一级作爱视频免费观看| 国产久久久一区二区三区| 国产午夜精品论理片| 1000部很黄的大片| avwww免费| 国内毛片毛片毛片毛片毛片| 久99久视频精品免费| 别揉我奶头 嗯啊视频| 91九色精品人成在线观看| 成人特级黄色片久久久久久久| 人人妻人人澡欧美一区二区| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 性插视频无遮挡在线免费观看| 男人和女人高潮做爰伦理| 国产视频内射| 99久久成人亚洲精品观看| 一个人看视频在线观看www免费| 亚洲中文字幕一区二区三区有码在线看| 丁香欧美五月| 男女床上黄色一级片免费看| 亚洲国产欧美人成| 久久天躁狠狠躁夜夜2o2o| 狂野欧美白嫩少妇大欣赏| 成年免费大片在线观看| 婷婷丁香在线五月| av专区在线播放| 国产午夜精品论理片| 国产一区二区三区视频了| 亚洲成人精品中文字幕电影| 精品99又大又爽又粗少妇毛片 | 99国产极品粉嫩在线观看| 999久久久精品免费观看国产| 国产三级中文精品| 久久精品91蜜桃| 自拍偷自拍亚洲精品老妇| 午夜a级毛片| 成人特级av手机在线观看| 亚洲av熟女| 亚洲成av人片免费观看| 久久久国产成人精品二区| 免费在线观看亚洲国产| 亚洲avbb在线观看| 美女被艹到高潮喷水动态| 激情在线观看视频在线高清| 国产男靠女视频免费网站| 他把我摸到了高潮在线观看| 国产av不卡久久| 日日干狠狠操夜夜爽| 性色av乱码一区二区三区2| 在线观看美女被高潮喷水网站 | 俺也久久电影网| 亚洲国产色片| 国产综合懂色| 在线观看av片永久免费下载| 国产精品综合久久久久久久免费| 男女下面进入的视频免费午夜| 如何舔出高潮| 别揉我奶头 嗯啊视频| 日韩国内少妇激情av| 九九在线视频观看精品| 国产蜜桃级精品一区二区三区| 有码 亚洲区| 亚洲av第一区精品v没综合| 日本一本二区三区精品| 亚洲av成人精品一区久久| 看黄色毛片网站| 老熟妇乱子伦视频在线观看| 亚洲精品在线美女| 五月伊人婷婷丁香| 色吧在线观看| 国产一区二区三区在线臀色熟女| 成年人黄色毛片网站| 99久国产av精品| 看免费av毛片| 国产精品1区2区在线观看.| 日韩免费av在线播放| 五月伊人婷婷丁香| 亚洲国产欧美人成| 久久性视频一级片| 欧美中文日本在线观看视频| 2021天堂中文幕一二区在线观| 亚洲激情在线av| 91在线精品国自产拍蜜月| 18禁黄网站禁片午夜丰满| 亚洲国产精品久久男人天堂| 亚洲在线观看片| 国语自产精品视频在线第100页| 两个人的视频大全免费| 少妇人妻一区二区三区视频| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美国产在线观看| 白带黄色成豆腐渣| 美女高潮喷水抽搐中文字幕| 久久午夜亚洲精品久久| 亚洲人成电影免费在线| 亚洲在线自拍视频| .国产精品久久| 色综合亚洲欧美另类图片| 国产69精品久久久久777片| 国产av在哪里看| 亚洲自偷自拍三级| 俺也久久电影网| 色综合婷婷激情| 一本综合久久免费| 两个人的视频大全免费| 最新中文字幕久久久久| 波多野结衣高清作品| 欧美乱妇无乱码| 国产精品久久久久久久久免 | 宅男免费午夜| 免费观看人在逋| 国产美女午夜福利| 亚洲男人的天堂狠狠| 国产伦精品一区二区三区四那| 免费在线观看影片大全网站| 一进一出好大好爽视频| 国产国拍精品亚洲av在线观看| 成人午夜高清在线视频| 日本撒尿小便嘘嘘汇集6| 成人欧美大片| 国产黄片美女视频| 男人舔奶头视频| 亚洲18禁久久av| 免费大片18禁| 最新在线观看一区二区三区| 精品国内亚洲2022精品成人| 激情在线观看视频在线高清| 久久热精品热| 最近中文字幕高清免费大全6 | 日本与韩国留学比较| 在线免费观看的www视频| 亚洲成人精品中文字幕电影| 国产精品一及| 成人特级av手机在线观看| 观看美女的网站| 成人午夜高清在线视频| 91久久精品国产一区二区成人| 啪啪无遮挡十八禁网站| 亚洲人成伊人成综合网2020| 波多野结衣巨乳人妻| 97碰自拍视频| 亚州av有码| 嫩草影院入口| 99久国产av精品| 韩国av一区二区三区四区| 在线天堂最新版资源| 国产欧美日韩一区二区三| 欧美zozozo另类| 国内精品久久久久久久电影| 欧美潮喷喷水| 日本成人三级电影网站| 乱码一卡2卡4卡精品| 亚洲在线观看片| 久久香蕉精品热| 亚洲乱码一区二区免费版| 亚洲精品粉嫩美女一区| 亚洲国产精品成人综合色| 欧美色欧美亚洲另类二区| 午夜日韩欧美国产| 久久精品人妻少妇| 一区二区三区四区激情视频 | 在线观看一区二区三区| 我要搜黄色片| xxxwww97欧美| 日日夜夜操网爽| 成年女人看的毛片在线观看| 免费看a级黄色片| 一区二区三区四区激情视频 | 3wmmmm亚洲av在线观看| 国产精品三级大全| 非洲黑人性xxxx精品又粗又长| 怎么达到女性高潮| 亚洲成人久久性| 琪琪午夜伦伦电影理论片6080| 久99久视频精品免费| 欧美一级a爱片免费观看看| 亚洲第一区二区三区不卡| 99在线人妻在线中文字幕| 18禁裸乳无遮挡免费网站照片| 欧美国产日韩亚洲一区| 国产精品亚洲美女久久久| 国产三级中文精品| 久久久久久久午夜电影| 精品一区二区三区视频在线| 99在线视频只有这里精品首页| 日韩有码中文字幕| 久久久久性生活片| 欧美高清性xxxxhd video| 日韩欧美在线二视频| 淫妇啪啪啪对白视频| 伦理电影大哥的女人| 少妇高潮的动态图| 九九在线视频观看精品| 好男人电影高清在线观看| 国产真实伦视频高清在线观看 | 欧美成人免费av一区二区三区| 中出人妻视频一区二区| 看片在线看免费视频| 亚洲最大成人中文| 国产精品亚洲一级av第二区| 日韩欧美 国产精品| 人人妻,人人澡人人爽秒播| 久久精品国产自在天天线| 91麻豆精品激情在线观看国产| 精品一区二区三区视频在线| 真实男女啪啪啪动态图| av天堂在线播放| 午夜福利免费观看在线| 国产成人aa在线观看| 老司机深夜福利视频在线观看| 日本与韩国留学比较| 舔av片在线| 别揉我奶头 嗯啊视频| 婷婷精品国产亚洲av| 美女大奶头视频| 亚洲精品久久国产高清桃花| 久久久精品大字幕| 国产精品一区二区性色av| 超碰av人人做人人爽久久| 成人鲁丝片一二三区免费| 黄色一级大片看看| 精品国产三级普通话版| 日韩欧美精品免费久久 | 少妇的逼水好多| 国产乱人视频| 深爱激情五月婷婷| 国产精品伦人一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 欧美又色又爽又黄视频| 俄罗斯特黄特色一大片| 国产久久久一区二区三区| 精品福利观看| 亚洲激情在线av| 性欧美人与动物交配| 亚洲国产精品合色在线| 一a级毛片在线观看| 午夜福利在线观看吧| 久久久久免费精品人妻一区二区| 欧美在线黄色| 亚洲欧美清纯卡通| 老司机深夜福利视频在线观看| 搡女人真爽免费视频火全软件 | 亚洲精品色激情综合| 高清毛片免费观看视频网站| 精品国内亚洲2022精品成人| 人人妻,人人澡人人爽秒播| 日本a在线网址| 精品久久久久久久久久久久久| 人妻夜夜爽99麻豆av| 亚洲aⅴ乱码一区二区在线播放| 高清在线国产一区| 免费在线观看影片大全网站| 国产国拍精品亚洲av在线观看| 色综合欧美亚洲国产小说| 国产精品亚洲一级av第二区| 国产一区二区亚洲精品在线观看| 内地一区二区视频在线| 久久香蕉精品热| 国产欧美日韩精品亚洲av| 色综合欧美亚洲国产小说| 亚洲人成网站在线播放欧美日韩| 美女高潮喷水抽搐中文字幕| 天天躁日日操中文字幕| www.熟女人妻精品国产| 成人av一区二区三区在线看| 男女之事视频高清在线观看| 美女被艹到高潮喷水动态| 每晚都被弄得嗷嗷叫到高潮| 久久久精品欧美日韩精品| 欧美中文日本在线观看视频| av福利片在线观看| 变态另类丝袜制服| x7x7x7水蜜桃| 小说图片视频综合网站| www.999成人在线观看| 极品教师在线视频| 免费无遮挡裸体视频| 97碰自拍视频| 久久人人爽人人爽人人片va | 欧美精品国产亚洲| 狠狠狠狠99中文字幕| 美女高潮的动态| 欧美激情国产日韩精品一区| 老熟妇仑乱视频hdxx| 日韩欧美三级三区| 久久九九热精品免费| 国产色爽女视频免费观看| 国产高清视频在线播放一区| 丰满人妻一区二区三区视频av| 最近最新免费中文字幕在线| 男人狂女人下面高潮的视频| 婷婷亚洲欧美| 成人毛片a级毛片在线播放| 欧美三级亚洲精品| 全区人妻精品视频| 99国产极品粉嫩在线观看| 亚洲国产欧洲综合997久久,| 精品一区二区三区av网在线观看| 一级a爱片免费观看的视频| 男人的好看免费观看在线视频| 国语自产精品视频在线第100页| 亚洲av电影在线进入| 亚洲av中文字字幕乱码综合| 在线观看一区二区三区| 精品一区二区三区视频在线| 日本黄大片高清| 久久亚洲真实| 欧美精品啪啪一区二区三区| 99热只有精品国产| 一级黄片播放器| 国产免费一级a男人的天堂| 久久久久久久精品吃奶| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区久久| 精品午夜福利视频在线观看一区| 午夜福利在线在线| 真人一进一出gif抽搐免费| 亚洲中文字幕一区二区三区有码在线看| 毛片一级片免费看久久久久 | 亚州av有码| 欧美zozozo另类| 啦啦啦观看免费观看视频高清| 精品一区二区免费观看| 51国产日韩欧美| 国产亚洲精品久久久久久毛片| 男女之事视频高清在线观看| 熟妇人妻久久中文字幕3abv| 色哟哟·www| 日本一本二区三区精品|