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

    剛性變彎度機翼設(shè)計及低速風(fēng)洞試驗研究

    2021-11-17 12:07:04謝長川
    振動與沖擊 2021年21期
    關(guān)鍵詞:彎度迎角升力

    冒 森, 楊 超, 謝長川, 孟 楊

    (1.成都飛機工業(yè)(集團)有限責(zé)任公司技術(shù)中心,成都 510100;2.北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院,北京 100191)

    從萊特兄弟1903年第一架飛機的試飛成功開始,優(yōu)化飛機的氣動性能,實現(xiàn)更為高效、機動及安全的飛行一直是飛機研制的重點。未來飛行器朝飛行任務(wù)多樣化、飛行條件多變以及實時調(diào)節(jié)飛機的飛行性能等方向發(fā)展,傳統(tǒng)機翼已無法始終保持飛行過程中所需的最佳氣動外形。變彎度機翼無論從基礎(chǔ)的空氣動力學(xué)方面還是工程實用角度都具有極大潛力[1-2],采用變彎度機翼自適應(yīng)技術(shù)的飛行器能依據(jù)飛行環(huán)境和任務(wù)的變化實時調(diào)整機翼彎度來提高飛行效率,增強機動性,以適應(yīng)復(fù)雜多變的任務(wù)環(huán)境,被認為是未來實現(xiàn)航空技術(shù)新突破的主要研究方向之一。

    從以往的研究可以看出,變彎度機翼的研究方法主要有仿真計算和試驗研究兩方面。仿真計算具有成本低、速度快、精度高等優(yōu)點,在初始設(shè)計階段和后期結(jié)構(gòu)優(yōu)化方面發(fā)揮了重要作用。Fujiwara等[3]以可變后緣的遠程寬體跨音速客機為研究對象,采用CFD/CSD流固耦合計算并優(yōu)化了機翼后緣氣動外形和翼盒結(jié)構(gòu),使得巡航狀態(tài)下全機減少了4.7%的燃油消耗。聶雪媛等[4]通過相似的方法研究柔性飛行器大變形時的靜氣動彈性問題。李哲等[5]利用幾何拓撲的方法對柔性變彎度機翼結(jié)構(gòu)進行優(yōu)化設(shè)計。Mao等[6]基于幾何精確渦格法對采用壓電片驅(qū)動的變后緣二維翼段進行分析,很好地描述了壓電變彎度機翼的氣動非線性特性。倪迎鴿等[7]則是通過有理函數(shù)擬合非定常氣動力研究了具有鉸鏈間隙的折疊變形機翼的非線性氣動彈性時域響應(yīng)。但變彎度機翼結(jié)構(gòu)普遍難以建立精確的數(shù)學(xué)模型,并且變彎度機翼實際變形能力有限,這導(dǎo)致仿真計算結(jié)果和實際結(jié)果差距比較大。因此,仿真計算不能作為判斷變彎度機翼設(shè)計的唯一標準,需要和試驗分析相結(jié)合。

    試驗研究作為變形機翼主要的研究手段一直也是變彎度機翼的研究熱點。試驗主要有地面試驗,風(fēng)洞試驗和飛行試驗三種。地面試驗主要目的是驗證機翼的變彎度功能。Campanile等[8]基于柔順機構(gòu)提出了帶肋變形結(jié)構(gòu),用碳纖維和金屬鉸鏈制造帶肋結(jié)構(gòu),并在試驗中驗證了變形能力。Icardi等[9]設(shè)計并制作了形狀機翼合金(SMA)驅(qū)動的柔性蒙皮變體機翼模型,機翼結(jié)構(gòu)由夾心盒段、柔性蒙皮和柔性翼肋組成,能夠?qū)崿F(xiàn)機翼整體彎曲和局部表面變形。王曉宏[10]設(shè)計了一種利用形狀記憶合金驅(qū)動的可連續(xù)變彎度機翼機構(gòu),將SMA絲布置于波紋狀主動變形結(jié)構(gòu)的谷底位置,通過形狀記憶合金變形驅(qū)動整個結(jié)構(gòu)變形。這些地面試驗直觀地反映所設(shè)計機翼的變形能力,但是沒有引入氣動載荷,其真實變形能力需要后續(xù)的風(fēng)洞試驗進一步驗證。飛行試驗是驗證設(shè)計可行性的直接方案。NASA將可連續(xù)變彎度機翼機構(gòu)與機翼固定盒段連接成為一個整體,機翼表面覆蓋可變形蒙皮,將其運用到F-111飛機上,并進行了多次飛行驗證試驗,均取得了圓滿成功[11-12]。Probst等[13]采用MFC(Macro Fiber Composite)作為驅(qū)動器,在空中對展長為0.5 m的無人機進行彎曲控制。飛行試驗是變彎度機翼進行實際運用的最終手段,但是由于重復(fù)性差、準備周期長,并不適合于初步設(shè)計階段。風(fēng)洞試驗可以準確測量包括氣動特性和驅(qū)動功率,實際變形位移等變彎度機翼關(guān)注的物理量,具有重復(fù)性好、測控精度高等優(yōu)點,在試驗研究中承擔(dān)承前啟后的作用,因此絕大多數(shù)變彎度機翼的實驗研究都集中在這一階段。Woods等[14-15]提出了一種魚骨主動變彎度機翼結(jié)構(gòu)。風(fēng)洞性能測試表明,這種結(jié)構(gòu)能夠提供很大的向上向下偏轉(zhuǎn)變形的能力,并且對驅(qū)動能力要求較低。2015年,歐盟FP7項目開展了基于柔性結(jié)構(gòu)技術(shù)的機翼翼尖變形裝置的風(fēng)洞驗證[16]。Jenett等[17]公布一款模塊機翼,該機翼可實現(xiàn)翼展方向的連續(xù)扭轉(zhuǎn)變形。采用柔性臂實現(xiàn)翼尖扭轉(zhuǎn)驅(qū)動,并在風(fēng)洞試驗中證明其相對傳統(tǒng)機翼具有更加優(yōu)越的滾轉(zhuǎn)效率。劉逸峰等[18]設(shè)計了一種基于SMA驅(qū)動器厚度可變的機翼結(jié)構(gòu),并進行了地面加載和風(fēng)洞試驗,驗證了設(shè)計方案的可行性??梢钥闯鲈谘芯砍跗?,風(fēng)洞試驗可以快速準確的提供一些初始數(shù)據(jù),為變彎度機翼初步設(shè)計提供重要的參考。另一方面,不同形式的變彎度機翼在風(fēng)洞試驗需求上并沒有顯著區(qū)別,因此開發(fā)可靠準確的變彎度機翼風(fēng)洞試驗平臺,也是研究變彎度機翼的重要研究方向。

    對于變彎度機翼的設(shè)計方案,國內(nèi)外都開展了廣泛而深入的研究。從現(xiàn)有的研究成果來看,變彎度機翼設(shè)計可以大致分為基于傳統(tǒng)材料的機翼設(shè)計和基于新一代智能材料的機翼設(shè)計兩大類?,F(xiàn)階段在全尺寸飛機上得到飛行驗證的只有美國 FlexSys Inc.的無縫后緣變形方案[19-20]和歐盟Saristu的自適應(yīng)變形后緣方案[21-23]兩種,而這兩種方案都是基于傳統(tǒng)材料和常規(guī)驅(qū)動的變彎度機翼設(shè)計。側(cè)面反映出當(dāng)前階段智能變形機翼從技術(shù)成熟度和結(jié)構(gòu)可靠性角度還是有所欠缺的。因此,本文設(shè)計了一種基于傳統(tǒng)材料的剛性分段后緣變彎度機翼結(jié)構(gòu)(簡稱剛性變彎度機翼),以此作為風(fēng)洞試驗的研究平臺。該機翼結(jié)構(gòu)共分為四個翼盒,各個翼盒之間通過轉(zhuǎn)軸進行前后連接,并通過與轉(zhuǎn)軸相連的舵機實現(xiàn)機翼結(jié)構(gòu)彎度的變化。

    本文分為三部分。第1章介紹剛性變彎度機翼的詳細設(shè)計,通過氣動力計算軟件XFOIL對變彎度機翼的不同構(gòu)型進行數(shù)值仿真計算,研究其氣動性能。第2章設(shè)計變彎度機翼風(fēng)洞試驗平臺,搭建測試采集系統(tǒng),對剛性變彎度機翼進行風(fēng)洞試驗。第3章分析并研究仿真結(jié)果和風(fēng)洞測試結(jié)果的變化趨勢和差別。最后總結(jié)并評估剛性變彎度機翼的實際變形能力和風(fēng)洞測試平臺的性能。

    1 剛性變彎度機翼的設(shè)計與分析

    1.1 結(jié)構(gòu)設(shè)計

    本文提出的剛性分段后緣變彎度機翼結(jié)構(gòu)由四段剛性翼盒組成。區(qū)別于柔性變彎度機翼,本文提出的變彎度機翼變彎度功能不涉及材料變形,而是通過舵機驅(qū)動翼盒之間的轉(zhuǎn)軸實現(xiàn),結(jié)構(gòu)如圖1所示。主要設(shè)計包括以下部分:①變彎度機翼結(jié)構(gòu)段,包括機翼前緣段、機翼中段、機翼中后段以及機翼后緣段;②轉(zhuǎn)動機構(gòu),包括舵機、轉(zhuǎn)軸、聯(lián)軸器等;③角度測試機構(gòu),包括角度傳感器、可拆卸支架和轉(zhuǎn)動輪。各個機翼結(jié)構(gòu)段均是由縱墻、翼肋、蒙板等組成的翼盒結(jié)構(gòu),主梁安置在機翼前緣結(jié)構(gòu)段中。機翼結(jié)構(gòu)段之間通過部分翼肋與轉(zhuǎn)軸相連,通過舵機帶動與結(jié)構(gòu)段相連的轉(zhuǎn)軸實現(xiàn)機翼彎度變化。角度傳感器通過傳動輪與舵機相連,根據(jù)傳感器轉(zhuǎn)動角度和傳動比間接計算出舵機的轉(zhuǎn)動角度。機翼參考翼型是NACA0015,考慮到舵機的尺寸和安裝形式,以及后續(xù)改進方案的需求,將從前緣到后緣4段翼盒占總弦長的比例設(shè)置為40%,17.5%,17.5%,25%。

    (a)

    (c)

    1-機翼前緣結(jié)構(gòu)段;2-機翼中部結(jié)構(gòu)段;3-機翼中后部結(jié)構(gòu)段;4-機翼后緣結(jié)構(gòu)段;5-縱墻;6-翼肋;7-蒙板;8-轉(zhuǎn)軸;9-舵機;10-角度傳感器;11-可拆卸支架;12-傳動輪;13-普通翼肋;14-加強翼肋;15-舵機翼肋;16-蒙皮。

    圖1 剛性分段后緣變彎度結(jié)構(gòu)設(shè)計圖與實物圖

    Fig.1 Structural design drawing and physical drawing of the rigid variable camber wing

    為了保證剛性變彎度機翼為沒有縫隙的封閉機翼。剛性翼盒設(shè)計為通長結(jié)構(gòu),翼盒表面布置蒙板,只在兩個翼盒連接處留有細長的縫隙。翼盒之間的柔性蒙皮選擇有機硅精密薄膜,這種薄膜彈性模量比較小并且不容易斷裂,經(jīng)預(yù)拉伸后進行粘貼到翼盒連接縫隙上。

    給定機翼期望彎度,控制裝置輸出相應(yīng)的電信號給舵機在舵機的驅(qū)動下,轉(zhuǎn)軸帶動與之相連的加強翼肋轉(zhuǎn)動,從而使得與之相連的機翼結(jié)構(gòu)段實現(xiàn)轉(zhuǎn)動功能。同時,舵機帶動與之相連的角度傳感器轉(zhuǎn)動,角度傳感器讀取舵機的轉(zhuǎn)動角度并將其轉(zhuǎn)化為相應(yīng)電信號傳輸?shù)娇刂蒲b置,從而構(gòu)成一套完整的角度測試系統(tǒng)。

    1.2 氣動分析

    為了確定本文提出的剛性變彎度機翼的設(shè)計方案,使用XFOIL氣動軟件進行了初步分析,并比較剛性變彎度機翼,標準NACA 0015翼型以及傳統(tǒng)帶后緣控制面機翼的氣動性能。XFOIL代碼是基于勢流理論,加入黏性邊界層求解器來預(yù)測表面摩擦阻力和流動分離,提供了比無黏性代碼更完整的阻力預(yù)測,且可在一定范圍內(nèi)考慮氣流分離特性。XFOIL已被證明在計算二維翼型時與高精度的CFD計算結(jié)果具有可信度[24]。它的計算成本低,易于集成到MATLAB軟件中,因此為模型設(shè)計中的中等精度氣動力分析提供了一個理想的解決方案。

    XFOIL的輸入是翼型坐標和外部氣動條件。為了準確得到機翼主動變形后的翼型坐標,本文采用翼型參數(shù)化方法,將剛性變彎度機翼翼盒的相對轉(zhuǎn)角作為參數(shù),生成機翼變形后的翼型坐標。具體方法為:首先選取翼盒表面上相當(dāng)數(shù)量的坐標點為控制點,根據(jù)相對轉(zhuǎn)角對控制點進行坐標旋轉(zhuǎn)變換,生成新的控制點。并在翼盒關(guān)節(jié)處對控制點進行三次樣條插值處理,最后重新插值生成新的翼型坐標。另外選取M=0.2,Re=106的工況作為外部氣動力條件。這是無人機典型工況點,馬赫數(shù)足夠低,可以避免顯著的壓縮性影響,雷諾數(shù)足夠大,可以避免低雷諾數(shù)時阻力的高度非線性增長。

    本文中將前緣翼盒、第二段、第三段及后緣翼盒之間的相對轉(zhuǎn)角分別記為轉(zhuǎn)角1、轉(zhuǎn)角2、轉(zhuǎn)角3。每個翼盒向下偏轉(zhuǎn)4°作為剛性變彎度機翼的典型氣動外形,按轉(zhuǎn)角順序記做變彎度構(gòu)型4-4-4。為了進行有效的比較,取弦長相同,且傳統(tǒng)帶后緣舵面翼型的最大升力系數(shù)將和變彎度構(gòu)型4-4-4相同作為基準狀態(tài)。通過分析計算,選擇初始翼型為NACA0015的傳統(tǒng)舵面翼型為研究對象,舵面長度設(shè)定為25%弦長,舵面向下偏轉(zhuǎn)14.5°,達到兩種翼型具有相同的最大升力系數(shù)。分析時假設(shè)傳統(tǒng)帶后緣控制面翼型是一個沒有縫隙的封閉舵面,圖中及后文簡記為傳統(tǒng)舵面翼型;無控制面偏轉(zhuǎn)的NACA0015翼型簡記為初始翼型,各翼型如圖2所示。

    圖2 三種翼型的翼剖面Fig.2 Wing profiles of the three airfoil

    圖3(a)顯示后緣向下偏轉(zhuǎn)后,傳統(tǒng)舵面翼型和構(gòu)型4-4-4升力系數(shù)顯著大于初始翼型。剛性變彎度機翼迎角為11°時,升力系數(shù)達到最大1.44,之后發(fā)生失速,升力系數(shù)減少。傳統(tǒng)舵面翼型擁有相同的最大升力系數(shù),迎角為12°時到達最大點。圖3(b)顯示當(dāng)迎角大于10°時,阻力系數(shù)顯著增加。表明臨近失速時,通過增加迎角已經(jīng)不能提升翼型升力,反而會額外增加阻力??紤]到XFOIL無法在大迎角失速狀態(tài)下得到準確的阻力系數(shù),此時的數(shù)據(jù)僅做為趨勢參考。

    圖4顯示的是迎角與力矩系數(shù)的關(guān)系,負值代表是低頭,正值代表抬頭。迎角為0°時,初始翼型的力矩系數(shù)為0,而傳統(tǒng)舵面翼型和變彎度構(gòu)型4-4-4有很大低頭力矩系數(shù)。迎角大于-3.5°時,隨著迎角的增加,傳統(tǒng)舵面翼型和構(gòu)型4-4-4的低頭力矩系數(shù)減少,并且傳統(tǒng)舵面翼型低頭力矩系數(shù)大于構(gòu)型4-4-4的。當(dāng)迎角小于-3.5°,傳統(tǒng)舵面翼型的低頭力矩系數(shù)要小于構(gòu)型4-4-4的力矩系數(shù)。綜合升力系數(shù)與力矩系數(shù)的結(jié)果可以看出,對實際飛機而言,在正常飛行時,在同一迎角下,相對傳統(tǒng)舵面翼型,變彎度翼型能夠提供更大的升力而產(chǎn)生更小的低頭力矩,這是變彎度翼型在實際應(yīng)用中的優(yōu)勢之一。

    (a) 迎角與升力系數(shù)關(guān)系

    (b) 迎角與阻力系數(shù)關(guān)系圖3 二維翼型氣動系數(shù)與迎角的關(guān)系Fig.3 Relationship of the 2D airfoils between aerodynamiccoefficient and angle of attack

    圖4 二維翼型力矩系數(shù)與迎角的關(guān)系

    圖5為典型迎角下各翼型的壓力系數(shù)(Cp)分布,橫坐標為各翼型的弦線弧長。圖5(a)顯示傳統(tǒng)舵面翼型和剛性變彎曲構(gòu)型4-4-4的上下表面壓力差均大于初始翼型,所以升力系數(shù)更高。構(gòu)型4-4-4的Cp分布曲線整體上比較光滑,只在翼盒轉(zhuǎn)角處有細微的突變。而傳統(tǒng)舵面翼型在后緣部分則出現(xiàn)明顯的氣動壓力上升,變化幅度劇烈。不光滑的Cp分布也是阻力增大的原因,變彎度機翼可以通過柔性蒙皮技術(shù)進一步提高壓力分布的光滑度從而減少阻力,對比傳統(tǒng)舵面翼型具有氣動優(yōu)勢。圖5(b)可以看出,迎角為-4°時傳統(tǒng)舵面翼型的下表面壓力系數(shù)對比迎角為-3.5°時顯著減少,所以引起了氣動力矩的突變。同等情況下,剛性變彎度機翼Cp分布整體變化不大,前緣升力減小。這與圖4中氣動力矩分布有一致性,因此變彎度構(gòu)型4-4-4對比傳統(tǒng)舵面翼型力矩曲線隨迎角變化更平穩(wěn)。

    (a) 5°迎角,翼型壓力系數(shù)分布

    (b) -3.5°與-4°迎角,翼型壓力系數(shù)分布圖5 典型迎角下各翼型的壓力系數(shù)分布

    圖6表示三種翼型的升力系數(shù)與阻力系數(shù)關(guān)系圖6(a)為翼型極曲線,可以看出,當(dāng)升力系數(shù)小于1.165時,變彎度構(gòu)型4-4-4的阻力系數(shù)要明顯小于傳統(tǒng)舵面翼型,甚至小于初始翼型。當(dāng)升力系數(shù)大于1.165時,傳統(tǒng)舵面翼型和構(gòu)型4-4-4有相當(dāng)?shù)淖枇ο禂?shù)。此外,當(dāng)升力系數(shù)大于1.165時,剛性變彎度機翼的阻力系數(shù)顯著上升而此時升力系數(shù)增加緩慢,這對應(yīng)了機翼大迎角時的狀態(tài)。圖6(b)顯示了升阻比與升力系數(shù)的關(guān)系,該圖進一步顯示了用變彎度構(gòu)型4-4-4相對于初始翼型和傳統(tǒng)舵面翼型的優(yōu)越性。構(gòu)型4-4-4最大升阻比相對于初始機翼提高了34%,相對于傳統(tǒng)舵面翼型提高了72%。當(dāng)升力系數(shù)相同且都小于1.165時,構(gòu)型4-4-4對比其他兩種翼型都有較高的升阻比。這一點可以從圖6(c)上得到驗證,當(dāng)迎角小于5°時,構(gòu)型4-4-4的升阻比大于其他兩種翼型,此時構(gòu)型4-4-4的升力系數(shù)為1.165。值得注意的是,剛性變彎度翼型已顯示出較好的氣動性能,且此時未針對迎角狀態(tài)優(yōu)選配置各段偏角,而且對傳統(tǒng)舵面翼型也未考慮縫隙的作用。

    為進一步研究變彎度機翼的氣動性能,選取變彎度構(gòu)型4-4-4,3-3-2,4-0-2,1-0-3和0-0-2進行計算,得到升阻力與迎角的關(guān)系如圖7所示,實線表示了幾種構(gòu)型升阻比外包絡(luò)線。由圖7可知,不同變彎度構(gòu)型的最大升阻比對應(yīng)的迎角不盡相同。翼型轉(zhuǎn)角越大,最大升阻比對應(yīng)的迎角越小。另外還注意到,變彎度構(gòu)型越接近初始翼型,升阻比隨迎角變化越平緩;變彎度翼型中弦線彎曲程度越大,升阻比隨迎角變化越劇烈。圖7中各構(gòu)型的升阻比包絡(luò)線證明,可以通過調(diào)整各段轉(zhuǎn)角使得剛性變彎度機翼在較大的迎角范圍內(nèi)都能保持較高升阻比,這是傳統(tǒng)機翼不能做到的,滿足變彎度機翼始終保持最優(yōu)氣動外形的設(shè)計初衷??紤]到二維翼型和實際機翼的差別,剛性變彎度機翼的實際氣動性能可能會低于計算預(yù)測值,但是三種構(gòu)型的氣動系數(shù)均通過統(tǒng)一分析工具得到,其相對關(guān)系具有明顯的參考價值。

    (a) 翼型極曲線

    從上面的分析結(jié)果可以得到,本文設(shè)計的剛性變彎度機翼對比傳統(tǒng)舵面翼型具有更好的氣動性能,同等升力下,阻力更小,氣動力矩隨迎角變化更平穩(wěn),并且能在一個較大升阻比下擴展迎角范圍。提供的翼型參數(shù)化方法可以正確反映翼型變化,結(jié)合XFOIL氣動力計算軟件,可以很好的分析不同翼型的氣動力系數(shù),計算耗費小,計算速度快,作為變彎度機翼初步階段的氣動分析工具是可行的。

    圖7 剛性變彎度機翼不同構(gòu)型迎角與升阻比關(guān)系

    2 風(fēng)洞試驗設(shè)計

    針對本文研究需求設(shè)計了一種低成本的變彎度機翼風(fēng)洞測試專用平臺,可以便捷地改變機翼對來流地迎角,并對機翼在各種迎角來流下的升力、阻力、俯仰力矩進行測試,為變彎度機翼的氣動力辨識,控制系統(tǒng)建模,主動氣動彈性控制等研究提供初始數(shù)據(jù)。作為變彎度機翼研究的初始階段,本文重點測試剛性變彎度機翼不同構(gòu)型下的氣動力,對比仿真計算結(jié)果,評估其在氣動力作用下的真實變形能力。

    2.1 試驗?zāi)P兔枋?/h3>

    試驗?zāi)P蜑榫匦螜C翼,弦長360 mm,展長700 mm,如圖8所示。翼盒結(jié)構(gòu)為木質(zhì)材料,總質(zhì)量510 g。機翼主梁布置在弦長22.2%處,主梁為一通長的鋁管,壁厚2 mm,外徑20 mm,長度為800 mm。剛性變彎度機翼總共布置了三個舵機,兩個舵機跟翼盒轉(zhuǎn)軸同步直接驅(qū)動翼盒偏轉(zhuǎn),后緣翼盒舵機通過連桿結(jié)構(gòu)驅(qū)動最后一段翼盒轉(zhuǎn)動。翼盒中另外布置三個角度傳感器,通過傳動輪與轉(zhuǎn)軸相連。所有舵機都選用Futaba BLS173SVI型航模舵機,其驅(qū)動力可以達到6.8 kg/cm,并且重量和體積都相對較小。柔性蒙皮材料為100 μm厚PDMS有機硅精密薄膜。

    迎角調(diào)節(jié)盤可在±15°內(nèi)較為精確調(diào)節(jié)機翼迎角,調(diào)節(jié)精度由轉(zhuǎn)盤上的角度刻度盤確定,約為±0.1°。機翼軸套基座和迎角調(diào)節(jié)轉(zhuǎn)盤用螺栓螺母連接,其套筒上有插銷孔,用于固定剛性變彎度機翼的主軸,并防止其相對轉(zhuǎn)動。

    角鐵、測力天平安裝在風(fēng)洞框架后,風(fēng)洞試驗前調(diào)節(jié)迎角時,首先固定住軸套與迎角調(diào)節(jié)轉(zhuǎn)盤,隨后按照對應(yīng)的角度刻盤轉(zhuǎn)動迎角調(diào)節(jié)轉(zhuǎn)盤使測試翼段與來流達到一定迎角,再緊固相應(yīng)連接,實現(xiàn)迎角調(diào)節(jié)轉(zhuǎn)盤與測力天平一體固定。

    圖8 風(fēng)洞試驗支持系統(tǒng)Fig.8 Wind tunnel test support system

    2.2 數(shù)據(jù)采集與試驗監(jiān)控系統(tǒng)

    風(fēng)洞試驗相關(guān)設(shè)備主要包括穩(wěn)壓電源、低通濾波器、數(shù)據(jù)采集卡、工控機和監(jiān)控機等,整個系統(tǒng)示意圖見圖9(a),試驗實物圖見圖9(b)。穩(wěn)壓電源的作用是給測力天平供電。低通濾波器的作用是將天平平臺的采集信號中的高頻噪聲濾掉。數(shù)據(jù)采集卡以1 000 Hz的采樣頻率采集傳感器感受到的模擬電壓信號并將其轉(zhuǎn)換成數(shù)字信號,再將數(shù)字信號輸入到工控機,通過工控機計算產(chǎn)生控制指令傳給可變彎度機翼上的舵機,控制可變彎度機翼的翼型彎度。監(jiān)控機的作用是實時監(jiān)控角度數(shù)據(jù),保證試驗?zāi)P偷陌踩W罱K,數(shù)據(jù)采集系統(tǒng)可以將數(shù)據(jù)采集卡得到的試驗數(shù)據(jù)傳遞給上位機,通過Labview編寫的圖形界面直觀地展示數(shù)據(jù)。

    (a) 系統(tǒng)示意圖

    (b) 系統(tǒng)實物圖圖9 風(fēng)洞試驗數(shù)據(jù)采集與監(jiān)控系統(tǒng)Fig.9 Acquisition and monitoring system of wind tunnel test

    本次風(fēng)洞試驗主要通過六自由度測力天平監(jiān)測模型根部測力數(shù)據(jù),見圖10,具體參數(shù)見表1,滿足本文試驗的測試需求。

    圖10 風(fēng)洞試驗天平Fig.10 Balance of the wind tunnel test

    表1 測力天平參數(shù)Tab.1 Parameters of the force balance

    在整個試驗過程中,機翼的主軸與天平幾何中心重合,機翼翼型弦線與天平X方向中線重合,保證天平獲得的力的大小能代表機翼機體坐標系的受力。再由坐標系轉(zhuǎn)化獲得機翼的氣動力特性,對于本次試驗來說,具體轉(zhuǎn)換關(guān)系如下

    升力系數(shù):

    (1)

    阻力系數(shù):

    (2)

    力矩系數(shù):

    (3)

    式中:Fx為天平X方向測力數(shù)據(jù);Fy為天平Y(jié)方向測力數(shù)據(jù);Mz為天平Z方向的力矩數(shù)據(jù);l為主梁和氣動中心的距離;α為機翼迎角;ρ為試驗空氣密度;v為試驗段來流風(fēng)速;S為機翼參考面積。

    3 試驗結(jié)果分析與討論

    開啟風(fēng)洞,風(fēng)速設(shè)定為10 m/s,給舵機輸入控制信號改變機翼的三個轉(zhuǎn)角,通過監(jiān)控器監(jiān)控并記錄翼盒實際轉(zhuǎn)角,同時記錄天平測力數(shù)據(jù)。迎角選取0°、2°、4°三個角度進行變化,轉(zhuǎn)角選擇為0°、2°、4°、6°四個角度進行變化。從迎角0°,構(gòu)型0-0-0開始測量。首先固定迎角,然后按順序分別變動4次轉(zhuǎn)角3、轉(zhuǎn)角2和轉(zhuǎn)角1,每次測量僅獨立變動一個角度,共需測量3×4×4×4=192個狀態(tài)。運用式(1)~(3),將天平測力數(shù)據(jù)換算為剛性變形機翼的氣動力系數(shù),記為測力試驗值。通過迎角轉(zhuǎn)盤和各翼盒相對實測角度,按1.2節(jié)分析方法計算得到對應(yīng)構(gòu)型的氣動力系數(shù),記為實測角度解算值。實測角度和指令角度按照試驗順序排列如圖11所示。同時,按1.2節(jié)分析方法計算得到舵機輸入角度對應(yīng)構(gòu)型的氣動力系數(shù),記為仿真計算值。

    (a) 轉(zhuǎn)角1

    圖11為轉(zhuǎn)角實測角度和指令角度按照測試順序排列。整體上實測角度和指令角度符合較好,平均相對誤差為-4.32%,但是由于存在結(jié)構(gòu)間隙和摩擦等不利因素影響,翼盒實際轉(zhuǎn)角小于指令轉(zhuǎn)角。又觀察到機翼迎角越大,實測角度和指令角度差距越大,在轉(zhuǎn)角1和轉(zhuǎn)角3上體現(xiàn)得比較明顯。這是由于迎角越大,機翼承受的氣動力越大,航模舵機在大載荷下輸出角度不穩(wěn)定引起的。這是后續(xù)研究中需要改進的地方,這也說明了一個問題,開環(huán)形式的變彎度機翼很難保證精確的角度控制,即使是剛性變彎度機翼這種可靠性比較高的變彎度結(jié)構(gòu),需要通過有效的反饋控制才能到達預(yù)期的效果。

    圖12(a)為剛性變彎度機翼各構(gòu)型在不同迎角下的升力系數(shù),圖12(b)為0°迎角下各構(gòu)型升力系數(shù),圖12(c)為圖12(b)黑色虛線框中放大圖,顯示的是轉(zhuǎn)角1為2°時各構(gòu)型升力系數(shù)。從圖12(a)可以看出各構(gòu)型的升力系數(shù)隨著迎角增加而增大。三種結(jié)果隨轉(zhuǎn)角變化的趨勢是一致的,實測角度解算值略小于仿真計算結(jié)果,平均相對誤差為-5.81%,測力試驗值則與仿真計算值吻合度略低,測力試驗值和實測角度解算值的平均相對誤差為-8.92%。實測角度解算值小于仿真計算結(jié)果是因為轉(zhuǎn)角實測角度要小于指令角度。測力試驗值總體偏低,其差異來源于兩點,一方面XFOIL是基于面元法的氣動力計算軟件,二維翼型升力系數(shù)仿真計算結(jié)果會高于實際;另一方面由于三維效應(yīng)的存在,二維翼型仿真結(jié)果也會比實際機翼模型的氣動力系數(shù)高。

    圖12(b)分為四部分,每部分有16個變彎度機翼構(gòu)型,表示不同轉(zhuǎn)角1角度下的各構(gòu)型升力系數(shù)??梢钥闯?,變彎度機翼構(gòu)型升力系數(shù)隨著轉(zhuǎn)角1角度增加而增大。圖12(c)橫坐標為具體構(gòu)型。對比構(gòu)型2-0-0,2-0-2,2-0-4,2-0-6的升力系數(shù)變化,可以得出保持轉(zhuǎn)角1和轉(zhuǎn)角2不變,增加轉(zhuǎn)角3角度,變彎度構(gòu)型升力系數(shù)增加。當(dāng)構(gòu)型2-0-6變動到2-2-0構(gòu)型時,翼型的升力系數(shù)下降。這是由于變彎度機翼的轉(zhuǎn)角3減少了6°,而轉(zhuǎn)角2只增加了2°,轉(zhuǎn)角3角度變化對升力下降的影響超過了轉(zhuǎn)角2。這也是圖12(c)中升力系數(shù)連線出現(xiàn)起伏的原因,類似的情況還出現(xiàn)在構(gòu)型2-2-6到構(gòu)型2-4-0等。此外,從構(gòu)型2-0-0,2-2-0,2-4-0,2-6-0的升力系數(shù)變化,可以得出轉(zhuǎn)角1和轉(zhuǎn)角3不變,增加轉(zhuǎn)角2的角度,構(gòu)型升力系數(shù)也會增加。綜上,對于本文中的變彎度機翼,增加轉(zhuǎn)角角度,都會增加升力系數(shù)。原因是增加翼盒轉(zhuǎn)角,相當(dāng)于增加了局部迎角,使得翼型整體環(huán)量上升,從而引起升力增加。

    (a)

    為了研究各轉(zhuǎn)角變化對變彎度機翼構(gòu)型升力系數(shù)的具體影響,計算得到各轉(zhuǎn)角平均單位角度升力系數(shù)增量與轉(zhuǎn)角的關(guān)系,如圖13所示。從圖中看出,仿真計算結(jié)果和風(fēng)洞試驗值都表明三個轉(zhuǎn)角的平均單位角度升力系數(shù)增量都隨著轉(zhuǎn)角的增加而減小。測力試驗值的平均單位角度升力系數(shù)增量要小于實測角度解算值且小于仿真計算值,這和圖12(a)中測力試驗值升力系數(shù)要小于實測角度解算值和小于仿真計算值是一致的。相同轉(zhuǎn)角下,轉(zhuǎn)角1的單位角度升力系數(shù)增量大于轉(zhuǎn)角2大于轉(zhuǎn)角3。這是由于轉(zhuǎn)角1轉(zhuǎn)動,后續(xù)所有翼盒都保持相對轉(zhuǎn)角整體在偏轉(zhuǎn),而轉(zhuǎn)角3變化只相當(dāng)于后緣翼盒發(fā)生偏轉(zhuǎn)。不同轉(zhuǎn)角變化對應(yīng)不同的機翼面積偏轉(zhuǎn),所以變化轉(zhuǎn)角1帶來的機翼升力系數(shù)變化最大。因此對于剛性變彎度機翼,為了提高機翼升力增加轉(zhuǎn)角1角度收益最好。

    圖14為剛性變彎度機翼各構(gòu)型在不同迎角下的阻力系數(shù)。迎角和翼盒轉(zhuǎn)角越大,阻力系數(shù)越大,這與升力系數(shù)的變化趨勢一致。構(gòu)型6-6-6是各迎角下的阻力系數(shù)最大的構(gòu)型。各構(gòu)型的阻力系數(shù)的實測角度解算值和仿真計算值的平均相對誤差為-7.52%,測力試驗值和實測角度解算值的平均相對誤差為-13.22%,都要大于升力系數(shù)的誤差。阻力系數(shù)的角度推算值與仿真計算值誤差較大的主要原因是,XFOIL計算阻力系數(shù)時對剛性變彎度機翼構(gòu)型轉(zhuǎn)角誤差更為敏感。考慮到機翼的三維效應(yīng),并且阻力系數(shù)測力試驗值能正確反映轉(zhuǎn)角變化對于機翼阻力的影響,其與仿真計算值誤差是可以接受的。

    圖13 平均單位角度升力系數(shù)增量與轉(zhuǎn)角關(guān)系

    圖14 剛性變彎度機翼各構(gòu)型阻力系數(shù)

    圖15為各變形構(gòu)型下力矩系數(shù)試驗值和仿真計算值的比較。從中可以看出剛性變彎度機翼翼盒轉(zhuǎn)角增加,產(chǎn)生低頭力矩,隨著迎角增加,低頭力矩減少。最大的低頭力矩系數(shù)是0°迎角下構(gòu)型0-6-6產(chǎn)生的。變形機翼構(gòu)型的力矩系數(shù)的實測角度解算值和仿真計算值的平均相對誤差為7.28%,測力試驗值和實測角度解算值的平均相對誤差為6.77%,試驗值跟仿真計算結(jié)果吻合得更好。計算得到各轉(zhuǎn)角平均單位角度氣動力矩系數(shù)增量與轉(zhuǎn)角的關(guān)系,如圖16所示。

    圖15 剛性變彎度機翼各構(gòu)型力矩系數(shù)Fig.15 Moment coefficients of the configurations of the rigidvariable camber wing

    圖16顯示單位角度低頭力矩系數(shù)增量隨著轉(zhuǎn)角的增加而減少。轉(zhuǎn)角較小時,轉(zhuǎn)角1提供的單位角度低頭力矩系數(shù)增量最大,但是隨著轉(zhuǎn)角增加迅速減少。轉(zhuǎn)角較大時,轉(zhuǎn)角3的單位角度低頭力矩系數(shù)增量最大,而且隨角度變化不大。轉(zhuǎn)角2單位角度低頭力矩系數(shù)增量整體上小于轉(zhuǎn)角3。這意味想要通過改變剛性變彎度機翼的氣動力矩來調(diào)整飛機姿態(tài),最好的辦法是調(diào)整轉(zhuǎn)角3的角度。

    圖17為不同迎角下剛性變彎度機翼各構(gòu)型的升阻比。其中仿真計算值和實測角度推算值符合的很好,且都大于測力試驗值。這主要是因為前兩者是二維翼型的升阻比,忽略了機翼的三維效應(yīng),所以要大于測力試驗值得到的升阻比。實測角度推算結(jié)果顯示,0°迎角下4-4-2構(gòu)型升阻比最大為102.4,2°迎角下2-2-4構(gòu)型升阻比最大為85.32,4°迎角下構(gòu)型1-2-0升阻比最大為77.55。隨著迎角增加,變彎度機翼構(gòu)型最大升阻比減少;迎角越大,最大升阻比對應(yīng)的翼型轉(zhuǎn)角越小。這和本文第二小節(jié)圖7得到的結(jié)論是一致的。測力試驗值表明,0°迎角下4-4-2構(gòu)型升阻比最大為69.29,2°迎角下2-2-2構(gòu)型升阻比最大為68.94,4°迎角下構(gòu)型1-4-0升阻比最大為49.71。雖然測力試驗值得到的不同迎角下最大升阻比構(gòu)型與實測角度推算值有所不同,但是相關(guān)構(gòu)型轉(zhuǎn)角差別不大并且其變化趨勢一致,進一步驗證了本文提供的仿真計算方法的合理性。同時看到當(dāng)構(gòu)型轉(zhuǎn)角等于6°時升阻比都比較小。因此剛性變彎度機翼要保持較高的升阻比,進行轉(zhuǎn)角控制時盡量避免單一轉(zhuǎn)角角度過大的情況,迎角增加時需要減少轉(zhuǎn)角角度。

    圖16 平均單位角度力矩系數(shù)增量與轉(zhuǎn)角關(guān)系Fig.16 Relationship between average increment of momentcoefficients per unit angle and rotation angle

    圖17 剛性變彎度機翼各構(gòu)型升阻比

    結(jié)合上述分析,我們可以得到:① 本文設(shè)計的變形機翼風(fēng)洞試驗平臺滿足剛性后緣分段變彎度機翼氣動性能測試要求。對于升力系數(shù),實測角度解算值和仿真計算值氣動系數(shù)的誤差為-5.81%,測力試驗值和實測角度解算值誤差為-8.92%;對于阻力系數(shù),實測角度解算值和仿真計算值的平均相對誤差為-7.52%,測力試驗值和實測角度解算值的平均相對誤差為-13.22%;對于力矩系數(shù),實測角度解算值和仿真計算值的平均相對誤差為7.28%,測力試驗值和實測角度解算值的平均相對誤差為6.77%。總體來說,氣動力系數(shù)的測力試驗值的結(jié)果與仿真計算結(jié)果符合較好,能正確反映剛性變彎度機翼的轉(zhuǎn)角變化對氣動力的影響以及不同迎角下升阻比最大的變彎度機翼構(gòu)型。② 對于本文提供的剛性變彎度機翼,翼盒實際轉(zhuǎn)角和指令轉(zhuǎn)角平均相對誤差在5%以內(nèi),保證了結(jié)構(gòu)變形功能。對于其飛行控制而言,增加轉(zhuǎn)角1能顯著提高升力,增加轉(zhuǎn)角3則能顯著提高低頭力矩,避免單一轉(zhuǎn)角過大并且在迎角增加的情況下減少轉(zhuǎn)角可以保持高升阻比。

    4 結(jié) 論

    本文以一種剛性變彎度機翼為研究對象,將剛性變彎度機翼參數(shù)化后,結(jié)合XFOIL氣動力計算軟件,對變彎度機翼的不同構(gòu)型進行數(shù)值仿真計算。在此基礎(chǔ)上,設(shè)計變彎度機翼風(fēng)洞試驗平臺,搭建測試采集系統(tǒng),對剛性變彎度機翼進行低速風(fēng)洞試驗研究。分析仿真計算結(jié)果和風(fēng)洞試驗結(jié)果后,得出結(jié)論:

    (1) 本文提出的剛性變彎度機翼作為變彎度機翼的一種設(shè)計方案,結(jié)構(gòu)簡單可靠,變形效果顯著。仿真計算和風(fēng)洞試驗表明均表明其在增升減阻,保證高升阻比和提供機翼配平力矩的優(yōu)勢明顯,并提出了對應(yīng)的控制策略。但是實際控制中存在轉(zhuǎn)角變形不到位的情況,需要對其進行角度反饋控制。

    (2) 本文設(shè)計的剛性變彎度機翼風(fēng)洞試驗系統(tǒng)可以正確測量和記錄不同迎角和構(gòu)型下剛性變彎度機翼的氣動力。通過本次多工況的風(fēng)洞試驗,其結(jié)構(gòu)可靠性和易操作性都得到了驗證,為仿真計算提供了數(shù)據(jù)支撐。作為變彎度機翼風(fēng)洞試驗平臺可以在此基礎(chǔ)上進一步開展變形機翼的動力學(xué)試驗。

    (3) 將剛性變彎度機翼參數(shù)化后,結(jié)合XFOIL氣動力計算軟件,可以快速的計算不同剛性變彎度機翼構(gòu)型的氣動系數(shù)。其作為剛性變彎度機翼的氣動分析工具是高效可行的,能正確反映轉(zhuǎn)角變化對剛性變彎度機翼氣動力的影響,準確性得到了風(fēng)洞數(shù)據(jù)驗證。在后續(xù)研究中可以作為分析工具開展剛性變彎度機翼的氣動外形優(yōu)化工作。

    猜你喜歡
    彎度迎角升力
    高速列車車頂–升力翼組合體氣動特性
    連續(xù)變彎度翼型動態(tài)氣動特性數(shù)值模擬
    連續(xù)變迎角試驗數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    無人機升力測試裝置設(shè)計及誤差因素分析
    翼型彎度變化自適應(yīng)規(guī)律研究
    基于自適應(yīng)偽譜法的升力式飛行器火星進入段快速軌跡優(yōu)化
    機翼后緣連續(xù)變彎度對客機氣動特性影響
    升力式再入飛行器體襟翼姿態(tài)控制方法
    失速保護系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    獨自面對生銹的彎度
    久久久精品欧美日韩精品| 日本与韩国留学比较| 草草在线视频免费看| 亚洲av电影在线进入| 亚洲最大成人手机在线| 亚洲av成人精品一区久久| 国产精品一区二区三区四区免费观看 | 18禁裸乳无遮挡免费网站照片| 欧美中文日本在线观看视频| 亚洲最大成人中文| 亚洲成人精品中文字幕电影| 欧美性感艳星| 极品教师在线视频| 老熟妇乱子伦视频在线观看| www日本黄色视频网| 免费大片18禁| 老司机福利观看| 好看av亚洲va欧美ⅴa在| 亚洲 国产 在线| 久久精品影院6| 欧美乱妇无乱码| 观看免费一级毛片| 悠悠久久av| 日韩欧美精品v在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲不卡免费看| 亚洲不卡免费看| 国产 一区 欧美 日韩| 午夜精品久久久久久毛片777| 99热这里只有是精品在线观看 | 黄色日韩在线| 草草在线视频免费看| а√天堂www在线а√下载| 桃色一区二区三区在线观看| 午夜日韩欧美国产| 亚洲国产精品sss在线观看| 噜噜噜噜噜久久久久久91| 精品人妻熟女av久视频| 丁香六月欧美| 亚洲精品一区av在线观看| 成人特级av手机在线观看| 亚洲人成电影免费在线| 日日摸夜夜添夜夜添av毛片 | 亚洲最大成人av| 国产成人aa在线观看| 99热精品在线国产| 亚洲真实伦在线观看| 露出奶头的视频| 在线看三级毛片| a级一级毛片免费在线观看| 亚洲成人久久性| 天堂√8在线中文| 99在线视频只有这里精品首页| 十八禁国产超污无遮挡网站| 国产精品亚洲一级av第二区| www.999成人在线观看| 亚洲av免费高清在线观看| 国产在线精品亚洲第一网站| 午夜福利高清视频| av在线天堂中文字幕| xxxwww97欧美| 99久久九九国产精品国产免费| 中文字幕人成人乱码亚洲影| 18禁裸乳无遮挡免费网站照片| 级片在线观看| 大型黄色视频在线免费观看| 制服丝袜大香蕉在线| 成人鲁丝片一二三区免费| 国产伦一二天堂av在线观看| 国产精品av视频在线免费观看| 91久久精品国产一区二区成人| 日韩欧美在线二视频| 国产精品98久久久久久宅男小说| 国产成年人精品一区二区| 亚洲成av人片免费观看| 人人妻人人看人人澡| 亚洲av成人不卡在线观看播放网| 最近中文字幕高清免费大全6 | 国产大屁股一区二区在线视频| 欧美日韩综合久久久久久 | 久久精品影院6| 嫩草影院新地址| 在现免费观看毛片| 午夜老司机福利剧场| 亚洲av电影不卡..在线观看| 亚洲天堂国产精品一区在线| 亚洲国产日韩欧美精品在线观看| 在线天堂最新版资源| 国内揄拍国产精品人妻在线| 国产不卡一卡二| 噜噜噜噜噜久久久久久91| 久久久国产成人免费| 欧美成人性av电影在线观看| 两人在一起打扑克的视频| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品久久久com| 色综合欧美亚洲国产小说| 99久久精品一区二区三区| 在线十欧美十亚洲十日本专区| 国产午夜福利久久久久久| 一区福利在线观看| 国产成年人精品一区二区| 最近视频中文字幕2019在线8| 99久久无色码亚洲精品果冻| 国产精品嫩草影院av在线观看 | 丝袜美腿在线中文| 麻豆国产av国片精品| 女人被狂操c到高潮| 啦啦啦观看免费观看视频高清| 人妻制服诱惑在线中文字幕| 一进一出抽搐gif免费好疼| 99久久精品热视频| 免费观看人在逋| 国产淫片久久久久久久久 | 999久久久精品免费观看国产| 国产一区二区在线观看日韩| 午夜精品在线福利| 一个人看视频在线观看www免费| 3wmmmm亚洲av在线观看| 久久久久国内视频| 丁香六月欧美| 日本 av在线| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址| 少妇人妻一区二区三区视频| 久久久久九九精品影院| 免费一级毛片在线播放高清视频| 日韩中字成人| 麻豆成人午夜福利视频| 在线十欧美十亚洲十日本专区| 精品一区二区免费观看| 亚洲最大成人手机在线| 一本综合久久免费| 久久久久亚洲av毛片大全| 久久精品综合一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 精品人妻偷拍中文字幕| 免费人成视频x8x8入口观看| 窝窝影院91人妻| 午夜福利18| 欧美在线一区亚洲| 亚洲av中文字字幕乱码综合| 久久精品久久久久久噜噜老黄 | 麻豆国产av国片精品| 麻豆国产av国片精品| 亚洲avbb在线观看| 国产亚洲精品av在线| 亚洲国产精品sss在线观看| 性插视频无遮挡在线免费观看| 最近最新中文字幕大全电影3| 国产乱人视频| 亚洲五月天丁香| 淫妇啪啪啪对白视频| 亚洲精品亚洲一区二区| 国模一区二区三区四区视频| 国产亚洲精品av在线| 亚洲片人在线观看| 男人和女人高潮做爰伦理| 99久久99久久久精品蜜桃| 亚洲成人中文字幕在线播放| 97热精品久久久久久| 成人一区二区视频在线观看| 精华霜和精华液先用哪个| 日本一二三区视频观看| 国产免费男女视频| 国产黄a三级三级三级人| 国产精品亚洲av一区麻豆| 日本成人三级电影网站| 极品教师在线免费播放| 尤物成人国产欧美一区二区三区| 国产av在哪里看| 性插视频无遮挡在线免费观看| 免费观看精品视频网站| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品久久久com| 欧美日本视频| netflix在线观看网站| 一区福利在线观看| 国产日本99.免费观看| 久久久久国产精品人妻aⅴ院| 色吧在线观看| 97超级碰碰碰精品色视频在线观看| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 亚洲不卡免费看| 成人av一区二区三区在线看| 性插视频无遮挡在线免费观看| 窝窝影院91人妻| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区精品| 免费av不卡在线播放| 床上黄色一级片| 特大巨黑吊av在线直播| 日韩人妻高清精品专区| 性色avwww在线观看| 小蜜桃在线观看免费完整版高清| 亚洲激情在线av| 中文字幕精品亚洲无线码一区| 欧美激情久久久久久爽电影| 琪琪午夜伦伦电影理论片6080| 人妻久久中文字幕网| 身体一侧抽搐| 岛国在线免费视频观看| 在线国产一区二区在线| 精品午夜福利在线看| 老司机午夜十八禁免费视频| 小蜜桃在线观看免费完整版高清| 国产精品一区二区性色av| 中文字幕熟女人妻在线| 无人区码免费观看不卡| 国产精品人妻久久久久久| 日韩成人在线观看一区二区三区| 五月玫瑰六月丁香| 久久久久九九精品影院| 高清在线国产一区| av福利片在线观看| 男女做爰动态图高潮gif福利片| 丝袜美腿在线中文| 人妻久久中文字幕网| 亚洲色图av天堂| 九九久久精品国产亚洲av麻豆| 欧美在线一区亚洲| 色尼玛亚洲综合影院| 一级黄色大片毛片| 99久久精品热视频| 亚洲精品久久国产高清桃花| 真实男女啪啪啪动态图| 欧美日韩黄片免| 国产精品三级大全| 首页视频小说图片口味搜索| 欧美不卡视频在线免费观看| 久久人人爽人人爽人人片va | 国产精品综合久久久久久久免费| 午夜免费激情av| 成人永久免费在线观看视频| 亚洲无线在线观看| 怎么达到女性高潮| 欧美午夜高清在线| 欧美极品一区二区三区四区| 2021天堂中文幕一二区在线观| 嫁个100分男人电影在线观看| 真人做人爱边吃奶动态| 神马国产精品三级电影在线观看| 日韩欧美一区二区三区在线观看| 日本精品一区二区三区蜜桃| 免费大片18禁| 热99re8久久精品国产| 国产v大片淫在线免费观看| 国产高清三级在线| 亚洲成人免费电影在线观看| 国产精品免费一区二区三区在线| 国产精华一区二区三区| 免费观看精品视频网站| 亚洲 国产 在线| 偷拍熟女少妇极品色| 极品教师在线视频| 国产精品1区2区在线观看.| 色吧在线观看| 好男人在线观看高清免费视频| 日日夜夜操网爽| 精品人妻偷拍中文字幕| 在线观看美女被高潮喷水网站 | 亚洲av二区三区四区| 国产精品一及| 亚洲av电影不卡..在线观看| 美女大奶头视频| 日韩 亚洲 欧美在线| 国产亚洲精品久久久久久毛片| 性色av乱码一区二区三区2| 亚洲aⅴ乱码一区二区在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲乱码一区二区免费版| 精品欧美国产一区二区三| 欧美不卡视频在线免费观看| 岛国在线免费视频观看| 91在线观看av| 一个人免费在线观看的高清视频| 国产精品野战在线观看| 99在线人妻在线中文字幕| 久久精品国产亚洲av天美| 色在线成人网| 少妇人妻精品综合一区二区 | 女人十人毛片免费观看3o分钟| 国产精品野战在线观看| 变态另类丝袜制服| 成人美女网站在线观看视频| 久久精品人妻少妇| 国产一区二区三区在线臀色熟女| 99riav亚洲国产免费| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| 精品久久久久久久末码| 亚洲经典国产精华液单 | 波野结衣二区三区在线| 嫩草影院新地址| 午夜免费男女啪啪视频观看 | 国产主播在线观看一区二区| 婷婷丁香在线五月| 精品人妻视频免费看| 一进一出抽搐动态| 91字幕亚洲| 嫁个100分男人电影在线观看| 伦理电影大哥的女人| 极品教师在线免费播放| 国产熟女xx| 女生性感内裤真人,穿戴方法视频| 中文字幕人妻熟人妻熟丝袜美| 美女免费视频网站| 久久精品国产99精品国产亚洲性色| 日韩欧美 国产精品| 禁无遮挡网站| 在线天堂最新版资源| 91久久精品国产一区二区成人| 国产成+人综合+亚洲专区| 五月玫瑰六月丁香| 一本精品99久久精品77| 国产精品免费一区二区三区在线| 国产av一区在线观看免费| 中文字幕人成人乱码亚洲影| 极品教师在线视频| 丁香欧美五月| 精品日产1卡2卡| 久久国产乱子伦精品免费另类| 久久久久久国产a免费观看| 亚洲男人的天堂狠狠| 色在线成人网| aaaaa片日本免费| 91字幕亚洲| 国产人妻一区二区三区在| 不卡一级毛片| 午夜福利视频1000在线观看| 亚洲国产高清在线一区二区三| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av第一区精品v没综合| 免费观看精品视频网站| 床上黄色一级片| 十八禁国产超污无遮挡网站| 五月玫瑰六月丁香| 99热精品在线国产| 少妇的逼好多水| 亚洲av成人精品一区久久| 高清日韩中文字幕在线| 在线播放国产精品三级| 国产私拍福利视频在线观看| 91麻豆av在线| 精品人妻视频免费看| 日本黄色片子视频| av中文乱码字幕在线| 欧美日韩福利视频一区二区| 91字幕亚洲| 国产真实乱freesex| 国产白丝娇喘喷水9色精品| 国产亚洲欧美在线一区二区| 三级毛片av免费| 免费在线观看日本一区| 丁香欧美五月| 国产亚洲av嫩草精品影院| 99久久99久久久精品蜜桃| 深夜a级毛片| 给我免费播放毛片高清在线观看| 直男gayav资源| 成人特级黄色片久久久久久久| 午夜福利在线在线| 精华霜和精华液先用哪个| 国产成人av教育| av在线蜜桃| 中国美女看黄片| 人妻制服诱惑在线中文字幕| a级毛片免费高清观看在线播放| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| www.www免费av| 欧美+日韩+精品| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 婷婷色综合大香蕉| 真人一进一出gif抽搐免费| 在线十欧美十亚洲十日本专区| 一二三四社区在线视频社区8| 国产成人欧美在线观看| 国产一区二区三区视频了| 91在线精品国自产拍蜜月| av黄色大香蕉| 小说图片视频综合网站| 精品久久国产蜜桃| 男人狂女人下面高潮的视频| 不卡一级毛片| 亚洲av第一区精品v没综合| 成人一区二区视频在线观看| 少妇被粗大猛烈的视频| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡| 综合色av麻豆| 丰满人妻熟妇乱又伦精品不卡| av在线观看视频网站免费| 亚洲五月婷婷丁香| 一个人看视频在线观看www免费| 免费高清视频大片| 极品教师在线视频| 91字幕亚洲| 国产亚洲av嫩草精品影院| 黄色配什么色好看| 色哟哟·www| 国产在线精品亚洲第一网站| 久久人人爽人人爽人人片va | 不卡一级毛片| 国产美女午夜福利| 久久国产乱子伦精品免费另类| 欧美中文日本在线观看视频| 欧美bdsm另类| 国产av一区在线观看免费| 亚洲,欧美精品.| 看免费av毛片| 乱人视频在线观看| 老司机深夜福利视频在线观看| 亚洲av一区综合| 757午夜福利合集在线观看| 日韩av在线大香蕉| 熟女电影av网| 99精品久久久久人妻精品| 91字幕亚洲| 美女高潮的动态| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 亚州av有码| 精品国内亚洲2022精品成人| 90打野战视频偷拍视频| 搞女人的毛片| 久久亚洲真实| 身体一侧抽搐| 久久99热6这里只有精品| 精品人妻视频免费看| 亚洲男人的天堂狠狠| 国产91精品成人一区二区三区| 精品久久久久久久久久免费视频| 久久久久久久午夜电影| 少妇熟女aⅴ在线视频| 99国产精品一区二区蜜桃av| 成年女人看的毛片在线观看| 熟女人妻精品中文字幕| 亚洲精品粉嫩美女一区| 精品久久久久久久久久免费视频| 久久久久久久午夜电影| 天天一区二区日本电影三级| 最后的刺客免费高清国语| 亚洲中文字幕一区二区三区有码在线看| 熟女人妻精品中文字幕| 亚洲自拍偷在线| 欧美日韩亚洲国产一区二区在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| 精品久久久久久久久亚洲 | 欧美潮喷喷水| 俄罗斯特黄特色一大片| 免费无遮挡裸体视频| 国产亚洲精品久久久com| 91午夜精品亚洲一区二区三区 | 国产精品亚洲美女久久久| 欧美激情在线99| 一区二区三区激情视频| 亚洲人与动物交配视频| 3wmmmm亚洲av在线观看| 亚洲在线自拍视频| or卡值多少钱| 国产精品98久久久久久宅男小说| 麻豆成人午夜福利视频| 色在线成人网| 国产aⅴ精品一区二区三区波| 国产国拍精品亚洲av在线观看| 麻豆一二三区av精品| 天美传媒精品一区二区| 一个人看的www免费观看视频| 给我免费播放毛片高清在线观看| 精品人妻1区二区| 亚洲avbb在线观看| 亚洲18禁久久av| 国产亚洲精品av在线| 啪啪无遮挡十八禁网站| 永久网站在线| 首页视频小说图片口味搜索| 色综合亚洲欧美另类图片| 窝窝影院91人妻| 国内毛片毛片毛片毛片毛片| 欧美区成人在线视频| 可以在线观看的亚洲视频| 国产野战对白在线观看| 国产av一区在线观看免费| АⅤ资源中文在线天堂| 国产精品乱码一区二三区的特点| 给我免费播放毛片高清在线观看| av天堂在线播放| 欧洲精品卡2卡3卡4卡5卡区| 欧美绝顶高潮抽搐喷水| 国产精品久久视频播放| 人妻夜夜爽99麻豆av| 国产成人av教育| 久久精品久久久久久噜噜老黄 | 91久久精品电影网| 久久精品国产99精品国产亚洲性色| 欧美午夜高清在线| 欧美一级a爱片免费观看看| 综合色av麻豆| 变态另类成人亚洲欧美熟女| 精品人妻熟女av久视频| 最新中文字幕久久久久| 怎么达到女性高潮| 一个人免费在线观看电影| 久久香蕉精品热| 国产成人av教育| 精品国产三级普通话版| 99久久精品一区二区三区| 天美传媒精品一区二区| 精品国产亚洲在线| 99热这里只有是精品在线观看 | 丁香六月欧美| 深爱激情五月婷婷| 国产白丝娇喘喷水9色精品| 成人av在线播放网站| 国产真实伦视频高清在线观看 | 免费在线观看影片大全网站| 国产单亲对白刺激| 男女床上黄色一级片免费看| 老司机深夜福利视频在线观看| 国产成人影院久久av| 欧美又色又爽又黄视频| 中文字幕人妻熟人妻熟丝袜美| 色综合站精品国产| av福利片在线观看| 热99re8久久精品国产| 久久这里只有精品中国| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看| 12—13女人毛片做爰片一| 国产成人欧美在线观看| www.999成人在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲五月婷婷丁香| 深夜a级毛片| 亚洲欧美日韩高清在线视频| 欧美日韩乱码在线| 麻豆久久精品国产亚洲av| 人妻久久中文字幕网| www.999成人在线观看| 人妻制服诱惑在线中文字幕| 精品欧美国产一区二区三| av国产免费在线观看| 精品国内亚洲2022精品成人| 亚洲国产精品成人综合色| 国产高清视频在线播放一区| 一卡2卡三卡四卡精品乱码亚洲| 一二三四社区在线视频社区8| 757午夜福利合集在线观看| 久久久色成人| 少妇丰满av| 波多野结衣巨乳人妻| 性插视频无遮挡在线免费观看| 久久欧美精品欧美久久欧美| 又爽又黄a免费视频| 偷拍熟女少妇极品色| 国产一级毛片七仙女欲春2| 久久久国产成人精品二区| 亚洲18禁久久av| 又黄又爽又免费观看的视频| 欧美色视频一区免费| 久久精品夜夜夜夜夜久久蜜豆| 看黄色毛片网站| 97热精品久久久久久| 婷婷精品国产亚洲av在线| 观看免费一级毛片| 一个人免费在线观看电影| 每晚都被弄得嗷嗷叫到高潮| 日本三级黄在线观看| 嫩草影院精品99| 日本黄色视频三级网站网址| 色综合欧美亚洲国产小说| 久久精品综合一区二区三区| .国产精品久久| 美女高潮的动态| 免费高清视频大片| 婷婷亚洲欧美| 老司机午夜十八禁免费视频| 88av欧美| 亚洲最大成人中文| 久99久视频精品免费| 免费av毛片视频| 欧美日韩国产亚洲二区| 禁无遮挡网站| 男女下面进入的视频免费午夜| 天堂av国产一区二区熟女人妻| 色精品久久人妻99蜜桃| 国产午夜精品久久久久久一区二区三区 | 欧美性猛交黑人性爽| 亚洲中文字幕日韩| 亚洲av电影不卡..在线观看| 久久久久久久午夜电影| 日韩中文字幕欧美一区二区| 婷婷亚洲欧美| 国产中年淑女户外野战色| 亚洲精品色激情综合| 757午夜福利合集在线观看| 午夜激情欧美在线| 国产欧美日韩一区二区三| 国产伦一二天堂av在线观看| 简卡轻食公司| 亚洲av.av天堂| 亚洲一区二区三区不卡视频| 国产淫片久久久久久久久 | 免费看日本二区| 国产国拍精品亚洲av在线观看| 变态另类成人亚洲欧美熟女| 国产精品伦人一区二区| 成人特级av手机在线观看| 亚洲第一区二区三区不卡|