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

    電動飛機(jī)飛行品質(zhì)快速評估與驗(yàn)證

    2023-10-29 01:49:46楊鳳田周文雅
    計算機(jī)仿真 2023年9期
    關(guān)鍵詞:升降舵風(fēng)洞試驗(yàn)迎角

    張 磊,楊鳳田,周文雅,黃 俊

    (1. 遼寧通用航空研究院,遼寧 沈陽 110136;2. 沈陽航空航天大學(xué)遼寧省通用航空重點(diǎn)實(shí)驗(yàn)室,遼寧 沈陽110136;3. 沈陽飛機(jī)設(shè)計研究所,遼寧 沈陽110135;4. 大連理工大學(xué)航空航天學(xué)院 遼寧 大連116024;5. 北京航空航天大學(xué)航空科學(xué)與工程學(xué)院 北京100083)

    1 引言

    隨著中國采取更加有力的政策和措施,減少二氧化碳排放,逐步實(shí)現(xiàn)碳中和,新能源飛機(jī)的發(fā)展進(jìn)入了快車道[1,2]。區(qū)別于采用傳統(tǒng)燃油發(fā)動機(jī)的飛機(jī),電動飛機(jī)采用電機(jī)提供推力,電池提供能源,具有綠色環(huán)保的優(yōu)點(diǎn)[3-5]。在飛行器初步設(shè)計階段,通常需要快速計算出飛行器氣動參數(shù),進(jìn)而對飛行器穩(wěn)定性和操縱性進(jìn)行快速評估[6],并為飛行控制律設(shè)計提供基礎(chǔ)數(shù)據(jù)。因此,各種飛行器氣動計算和仿真軟件被開發(fā)出來。其中較為常用的氣動力估算軟件包括DATCOM軟件[7-9]。盡管飛行器氣動參數(shù)計算和仿真軟件功能越來越強(qiáng)大,計算能力越來越強(qiáng),但飛行器風(fēng)洞試驗(yàn)仍然是不可或缺的重要環(huán)節(jié)[10],尤其是在需要精確氣動參數(shù)數(shù)據(jù)的詳細(xì)設(shè)計階段[11]。盡管采用氣動計算仿真軟件和風(fēng)洞試驗(yàn)方法,可以得到較為準(zhǔn)確的飛機(jī)氣動參數(shù)和飛行性能,但飛行試驗(yàn)仍然是不可或缺的驗(yàn)證手段[12-14]。

    本文以沈陽航空航天大學(xué)和遼寧通用航空研究院研制的RX1E-A電動飛機(jī)為基準(zhǔn)模型,采用DATCOM氣動估算軟件計算其氣動參數(shù)和飛行品質(zhì)。再采用RX1E-A飛機(jī)縮比驗(yàn)證機(jī)在中航工業(yè)氣動院FL-8風(fēng)洞[15]所做風(fēng)洞試驗(yàn)氣動數(shù)據(jù),計算了該型飛機(jī)飛行品質(zhì),進(jìn)而驗(yàn)證DATCOM計算的氣動參數(shù)和飛行品質(zhì)的準(zhǔn)確程度。最后對RX1E-A飛機(jī)進(jìn)行飛行試驗(yàn),采集飛機(jī)舵偏產(chǎn)生的姿態(tài)角變化曲線,并根據(jù)飛行員評價,以此驗(yàn)證DATCOM計算的操穩(wěn)特性的準(zhǔn)確程度。

    RX1E-A電動飛機(jī)為上單翼常規(guī)布局螺旋槳電動力飛機(jī),RX1E-A電動飛機(jī)的基本參數(shù)如表1所示。

    表1 RX1E-A電動飛機(jī)基本參數(shù)

    2 電動飛機(jī)運(yùn)動建模

    本文著重推導(dǎo)繞質(zhì)心轉(zhuǎn)動動力學(xué)方程[16]。在地面坐標(biāo)系中,有下式

    (1)

    式中:∑M為飛機(jī)所受外合力矩,L為飛機(jī)動量矩。

    則有下式

    (2)

    將動量矩L寫成下式

    (3)

    式中:J為飛機(jī)慣性張量,Jx、Jy和Jz為飛機(jī)對機(jī)體坐標(biāo)系各軸的轉(zhuǎn)動慣量,Jxy、Jyx、Jxz、Jzx、Jyz和Jzy為飛機(jī)對機(jī)體坐標(biāo)系各軸的慣量積。

    本文飛機(jī)建模時,沒有按照常規(guī)的飛機(jī)為面對稱飛行器的假設(shè),而是把飛機(jī)出現(xiàn)機(jī)體不對稱、質(zhì)量分布不均勻的情況考慮進(jìn)來。假設(shè)飛機(jī)慣性張量J在某一段時間內(nèi)為時不變的常量,(2)式第一項(xiàng)可以可寫成

    (4)

    且,(2)式第二項(xiàng)可以可寫成

    Ω×L=S(Ω)JΩ

    (5)

    其中:S(Ω)定義如下

    (6)

    式中:Ω定義見(7)式。

    Ω=pib+qjb+rkb

    (7)

    式中:p為滾轉(zhuǎn)角速度,q為俯仰角速度,r為偏航角速度。

    將外合力矩∑M在機(jī)體坐標(biāo)系中分解,有

    (8)

    將(4)式、(5)式代入(1)式,可得

    (9)

    對(9)式兩邊慣性張量J求逆,則可得到在機(jī)體坐標(biāo)系中建立的基于不對稱飛機(jī)慣性張量的矩陣形式繞質(zhì)心轉(zhuǎn)動動力學(xué)方程

    (10)

    3 風(fēng)洞試驗(yàn)對DATCOM計算氣動參數(shù)驗(yàn)證

    DATCOM氣動力工程估算軟件是為快速計算飛行器氣動參數(shù)而由美國開發(fā)的,對于常規(guī)構(gòu)型飛行器均能得到比較理想的氣動數(shù)據(jù)。除此以外,DATCOM程序還可以計算高超聲速飛行器氣動舵效參數(shù),程序計算時的限制條件是迎角需在0°到20°范圍內(nèi),馬赫數(shù)大于5;DATCOM程序還可以計算橫向射流控制效果(Transverse-jet Control Effectiveness),程序計算時要求馬赫數(shù)范圍為2馬赫到20馬赫。

    將DATCOM軟件計算的RX1E-A飛機(jī)氣動參數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行對比,升力系數(shù)隨迎角變化曲線如圖1所示。

    圖1 升力系數(shù)隨迎角變化曲線

    DATCOM阻力系數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)對比如圖2所示。

    圖2 DATCOM阻力系數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)對比圖

    俯仰力矩系數(shù)隨迎角變化曲線如圖3所示。

    圖3 俯仰力矩系數(shù)隨迎角變化曲線

    滾轉(zhuǎn)力矩系數(shù)隨側(cè)滑角變化曲線如圖4所示。

    圖4 滾轉(zhuǎn)力矩系數(shù)隨側(cè)滑角變化曲線

    副翼產(chǎn)生的滾轉(zhuǎn)力矩系數(shù)隨側(cè)滑角變化曲線如圖5所示。

    圖5 副翼產(chǎn)生的滾轉(zhuǎn)力矩系數(shù)隨側(cè)滑角變化曲線

    升降舵產(chǎn)生的俯仰力矩系數(shù)隨迎角變化曲線如圖6所示。

    圖6 升降舵產(chǎn)生的俯仰力矩系數(shù)隨迎角變化曲線

    從圖1可知:DATCOM計算的升力系數(shù)在迎角處于-4°到1°范圍內(nèi)與風(fēng)洞試驗(yàn)數(shù)據(jù)相近,在迎角處于1°到15°范圍內(nèi)比風(fēng)洞試驗(yàn)數(shù)據(jù)小,在迎角為9°時兩者差值最大,DATCOM升力系數(shù)為風(fēng)洞數(shù)據(jù)的77%。

    從圖2可知,DATCOM修正后的阻力系數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)接近,尤其是0度到8度迎角范圍內(nèi)幾乎重合。

    從圖3可知:DATCOM俯仰力矩系數(shù)在迎角處于-4°到0°范圍時相差不大,在0°到16°范圍內(nèi)下降速率大約是風(fēng)洞試驗(yàn)數(shù)據(jù)下降速率的一半,雖然此時DATCOM升力系數(shù)比風(fēng)洞試驗(yàn)小一些,但筆者認(rèn)為主要原因還是由于DATCOM升力作用點(diǎn)隨迎角增加而后移的速度較小。

    從圖4可知:DATCOM滾轉(zhuǎn)力矩系數(shù)隨側(cè)滑角變化速率為風(fēng)洞試驗(yàn)數(shù)據(jù)1.5倍。

    從圖5可知:DATCOM計算得到的副翼偏轉(zhuǎn)產(chǎn)生的滾轉(zhuǎn)力矩系數(shù)約為風(fēng)洞數(shù)據(jù)的77%。

    從圖6可知:DATCOM計算得到的升降舵偏轉(zhuǎn)產(chǎn)生的俯仰力矩系數(shù)增量約為風(fēng)洞數(shù)據(jù)的1.3倍。

    4 風(fēng)洞試驗(yàn)數(shù)據(jù)對DATCOM計算飛行品質(zhì)驗(yàn)證

    飛行品質(zhì)對保證飛機(jī)飛行任務(wù)性能和飛行安全起到關(guān)鍵作用。其中飛機(jī)本體飛行品質(zhì)是飛機(jī)操縱和飛行控制的基礎(chǔ)。一架飛機(jī)在設(shè)計時需要關(guān)注飛行員操縱起來的感受,是否安全、有效和好飛。

    GJB185-86《有人駕駛飛機(jī)(固定翼)飛行品質(zhì)》對我國軍用有人駕駛飛機(jī)在空中和地面的飛行和操縱品質(zhì)提出了具體要求,用于保證完成任務(wù)和飛行安全。

    電動飛機(jī)本體飛行品質(zhì)基本要求和油動飛機(jī)基本一致,但電動飛機(jī)相比于油動飛機(jī)有如下特點(diǎn):

    1)電動飛機(jī)使用動力電池作為能源,在電動飛機(jī)整個飛行階段,動力電池重量幾乎不發(fā)生改變,因此電動飛機(jī)在整個飛行階段,重量重心和轉(zhuǎn)動慣量都不發(fā)生改變。

    2)電動飛機(jī)依靠電機(jī)提供推力,電機(jī)在運(yùn)行時不產(chǎn)生二氧化碳等有害氣體,并且噪聲很低,因此電動飛機(jī)具有對環(huán)境友好的特點(diǎn)。

    3)電動飛機(jī)電機(jī)工作效率幾乎不隨飛行高度變化,飛行高度適應(yīng)性較好。

    4)由于現(xiàn)階段動力電池能量密度沒有燃油高,因此電動飛機(jī)通常設(shè)計成大展弦比,因此電動飛機(jī)的穩(wěn)定性通常較好,但往往會犧牲一部分操縱性。

    由于RX1E-A電動飛機(jī)重量不超過4500kg,最大法向過載不超過4.5,因此本文按照輕小類(QX)飛機(jī),針對航向階段(B種)整理其飛行品質(zhì)要求,對應(yīng)RX1E-A電動飛機(jī)的飛行品質(zhì)標(biāo)準(zhǔn)如表2所示。

    表2 RX1E-A電動飛機(jī)選用飛行品質(zhì)標(biāo)準(zhǔn)

    基于DATCOM計算氣動數(shù)據(jù)進(jìn)行RX1E-A飛機(jī)飛行品質(zhì)分析,并用基于風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行的RX1E-A飛機(jī)飛行品質(zhì)分析結(jié)果對其進(jìn)行驗(yàn)證,具體如表3所示。

    表3 基于 DATCOM與風(fēng)洞試驗(yàn)飛行品質(zhì)對比

    從表3可知,基于DATCOM計算數(shù)據(jù)的RX1E-A電動飛機(jī)的飛行品質(zhì),除了螺旋模態(tài)滿足標(biāo)準(zhǔn)3要求,其它模態(tài)都滿足標(biāo)準(zhǔn)1要求?;陲L(fēng)洞試驗(yàn)數(shù)據(jù)的RX1E-A電動飛機(jī)進(jìn)行飛行品質(zhì),除了長周期模態(tài)滿足標(biāo)準(zhǔn)2要求,其它模態(tài)都滿足標(biāo)準(zhǔn)1要求?;贒ATCOM計算得到的長周期模態(tài)阻尼比ζsp是基于風(fēng)洞試驗(yàn)得到的數(shù)據(jù)403%,基于DATCOM得到的螺旋模態(tài)倍幅時間T2S是基于風(fēng)洞試驗(yàn)得到的數(shù)據(jù)的11%,其它數(shù)據(jù)都近似相等。

    RX1E-A電動飛機(jī)滾轉(zhuǎn)操縱性能也是按照輕小類(QX)飛機(jī)航向階段(B種)要求進(jìn)行判定,如表4所示。

    表4 RX1E-A電動飛機(jī)滾轉(zhuǎn)操縱性能要求

    基于風(fēng)洞試驗(yàn)數(shù)據(jù),RX1E-A飛機(jī)副翼給定-1°階躍輸入和-1°舵偏角持續(xù)1s的短時輸入對滾轉(zhuǎn)角作用曲線如圖7所示。其中階躍響應(yīng)為一條斜線,滾轉(zhuǎn)角隨時間線性增加,增加速率為2°/s,短時舵偏響應(yīng)在開始到1s時間段線性增加,增加速率與階躍響應(yīng)相同,1s后響應(yīng)曲線近似為一條水平線,但略有上升,俯仰角振蕩周期約為20s。根據(jù)階躍響應(yīng)可知,副翼偏轉(zhuǎn)20°時,滾轉(zhuǎn)角從30°變化到-30°需要1.5s時間。滿足輕小類(QX)飛機(jī)航向階段(B種)滾轉(zhuǎn)操縱性能要求中60度滾轉(zhuǎn)操縱性能標(biāo)準(zhǔn)1。

    圖7 基于風(fēng)洞試驗(yàn)數(shù)據(jù)RX1E-A飛機(jī)-1°副翼偏角的階躍響應(yīng)和擾動響應(yīng)

    基于DATCOM數(shù)據(jù),RX1E-A飛機(jī)副翼給定-1°階躍輸入和-1°舵偏角持續(xù)1s的短時輸入對滾轉(zhuǎn)角作用曲線如圖8所示。其中階躍響應(yīng)在0s到5s時間段近似為一條斜線,滾轉(zhuǎn)角隨時間線性增加,增加速率為1.1°/s,短時舵偏響應(yīng)在開始到1s時間段線性增加,增加速率與階躍響應(yīng)相同,1s后響應(yīng)曲線近似為一條水平線,但略有下降,俯仰角振蕩周期約為20s。根據(jù)階躍響應(yīng)可知,副翼偏轉(zhuǎn)20°時,滾轉(zhuǎn)角從30°變化到-30°需要2.7s時間。滿足輕小類(QX)飛機(jī)航向階段(B種)滾轉(zhuǎn)操縱性能要求中60度滾轉(zhuǎn)操縱性能標(biāo)準(zhǔn)3。

    圖8 基于DATCOM RX1E-A飛機(jī)-1°副翼偏角的階躍響應(yīng)和擾動響應(yīng)

    將基于DATCOM數(shù)據(jù)和風(fēng)洞試驗(yàn)數(shù)據(jù)得到的滾轉(zhuǎn)操縱性能參數(shù)進(jìn)行對比驗(yàn)證,如表5所示??芍贒ATCOM數(shù)據(jù)的60°滾轉(zhuǎn)操縱時間是基于風(fēng)洞試驗(yàn)數(shù)據(jù)得到的參數(shù)的180%。

    表5 基于 DATCOM與風(fēng)洞試驗(yàn)滾轉(zhuǎn)操縱參數(shù)對比

    基于DATCOM計算氣動數(shù)據(jù),RX1E-A飛機(jī)升降舵給定-5°階躍輸入和-5°舵偏角持續(xù)1s的擾動輸入對俯仰角作用曲線如圖9所示。其中俯仰角階躍響應(yīng)經(jīng)過300s仍振蕩衰減,最終收斂至15°附近。迎角階躍響應(yīng)在開始后的1s快速變化,然后隨著俯仰角的振蕩而輕微振蕩。俯仰角擾動響應(yīng)經(jīng)過200s振蕩衰減,最終收斂至0°俯仰角。迎角擾動響應(yīng)在開始后的2s快速振蕩后收斂至0°附近,后期隨著俯仰角的振蕩而輕微振蕩。

    圖9 基于DATCOM計算數(shù)據(jù)RX1E-A飛機(jī)-1°升降舵偏角的階躍響應(yīng)和擾動響應(yīng)

    基于風(fēng)洞試驗(yàn)數(shù)據(jù),RX1E-A飛機(jī)升降舵給定-5°階躍輸入和-5°舵偏角持續(xù)1s的擾動輸入對俯仰角作用曲線如圖10所示。其中俯仰角階躍響應(yīng)經(jīng)過500s后仍然在振蕩衰減,經(jīng)過約1000s最終收斂至5°附近,迎角階躍響應(yīng)在開始后的1s快速變化,然后隨著俯仰角的振蕩而輕微振蕩。俯仰角擾動響應(yīng)經(jīng)過500s仍然在振蕩衰減,經(jīng)過約700s最終收斂至0°俯仰角,迎角擾動響應(yīng)在開始后的2.2s快速振蕩后收斂至0°附近,后期隨著俯仰角的振蕩而輕微振蕩。可知本機(jī)俯仰穩(wěn)定性收斂穩(wěn)定,滿足正常飛行操縱要求。

    圖10 基于風(fēng)洞試驗(yàn)數(shù)據(jù)RX1E-A飛機(jī)-1°升降舵偏角階躍響應(yīng)和擾動響應(yīng)

    從圖9到圖10可以看出,DATCOM計算數(shù)據(jù)的升降舵偏角階躍輸入產(chǎn)生的俯仰角響應(yīng)在經(jīng)過振蕩穩(wěn)定后是風(fēng)洞試驗(yàn)數(shù)據(jù)得到的俯仰角的3倍。

    5 飛行試驗(yàn)對DATCOM計算飛行品質(zhì)驗(yàn)證

    RX1E-A飛機(jī)進(jìn)行飛行試驗(yàn)照片如圖11所示。

    圖11 RX1E-A電動飛機(jī)飛行試驗(yàn)照片

    試飛時,在RX1E-A飛機(jī)安裝飛控采集系統(tǒng),采集實(shí)際飛行中升降舵偏角、飛機(jī)俯仰角和飛機(jī)迎角,其中飛行試驗(yàn)中升降舵偏轉(zhuǎn)作用產(chǎn)生的飛機(jī)俯仰角和迎角變化的曲線如圖12所示。

    圖12 RX1E-A飛機(jī)飛行試驗(yàn)升降舵偏轉(zhuǎn)產(chǎn)生的飛機(jī)俯仰角和迎角變化曲線

    圖12中,橫軸為飛機(jī)飛行時間,從本次飛行試驗(yàn)計時的第1150s開始,給升降舵一個不超過-5°偏角的偏轉(zhuǎn),可以看到飛機(jī)迎角從不到10°迅速上升到接近20°,經(jīng)過快速的振蕩后逐漸平穩(wěn),在升降舵偏角回中時,飛機(jī)迎角也迅速穩(wěn)定在10°附近。飛機(jī)俯仰角也從10°附近上升到接近20°,然后緩慢振蕩,在升降舵在1165s回中時,俯仰角仍然在緩慢振蕩,在1180s時俯仰角仍然沒有停止振蕩,只是振蕩收斂,幅度越來越小,振蕩周期約為17s。本次飛行試驗(yàn)中升降舵偏轉(zhuǎn)產(chǎn)生的迎角變化與圖10基于風(fēng)洞試驗(yàn)數(shù)據(jù)和圖9中基于DATCOM數(shù)據(jù)的MATLAB仿真的迎角變化趨勢一致,迎角變化屬于縱向短周期運(yùn)動,由于具有較大的靜穩(wěn)定恢復(fù)力矩,快速振蕩后馬上(幾秒以內(nèi))就能穩(wěn)定;飛行試驗(yàn)中升降舵偏轉(zhuǎn)產(chǎn)生的俯仰角與圖10基于風(fēng)洞試驗(yàn)數(shù)據(jù)和圖9中基于DATCOM數(shù)據(jù)的MATLAB仿真的俯仰角變化趨勢一致,俯仰角變化屬于縱向長周期運(yùn)動,振蕩雖收斂但緩慢,在升降舵回中后仍然需要較長時間(數(shù)十秒到幾百秒)恢復(fù)到穩(wěn)態(tài),這個過程可以看成飛機(jī)勢能與動能相互轉(zhuǎn)化過程(假設(shè)飛機(jī)推力與阻力相等條件下),在這個過程中迎角幾乎不發(fā)生變化,主要是飛機(jī)速度和高度的緩慢變化,由于飛機(jī)質(zhì)量較大,慣性較大,而起恢復(fù)和阻尼作用的氣動力相對較小,因此恢復(fù)過程非常緩慢。

    在RX1E-A飛機(jī)進(jìn)行飛行驗(yàn)證中,根據(jù)飛行員反饋,該機(jī)在給定短時副翼偏轉(zhuǎn)并回中后,滾轉(zhuǎn)角先增加,然后穩(wěn)定在一個角度,并緩慢增加,與圖7基于風(fēng)洞試驗(yàn)數(shù)據(jù)的MATLAB仿真的變化趨勢一致,但不同于圖8中基于DATCOM數(shù)據(jù)的MATLAB仿真曲線。

    6 結(jié)論

    本文采用DATCOM氣動計算軟件,快速計算了RX1E-A飛機(jī)氣動參數(shù)和飛行品質(zhì),并采用風(fēng)洞試驗(yàn)數(shù)據(jù)和飛行試驗(yàn)數(shù)據(jù)對其進(jìn)行驗(yàn)證,得到以下結(jié)論:

    1)使用DATCOM計算飛機(jī)氣動參數(shù)與風(fēng)洞試驗(yàn)數(shù)據(jù)相比具有一定偏差,文中所列舉各氣動參數(shù)偏差在50%范圍內(nèi),但基本滿足初步設(shè)計階段的需求。

    2)基于DATCOM計算數(shù)據(jù)得到的飛機(jī)飛行品質(zhì)與基于風(fēng)洞試驗(yàn)數(shù)據(jù)得到飛行品質(zhì),在長周期模態(tài)ζp和螺旋模態(tài)T2S有較大偏差,在其它模態(tài)偏差很小,滾轉(zhuǎn)操縱參數(shù)偏差可達(dá)80%。

    3)由DATCOM計算仿真得到的升降舵階躍響應(yīng)和擾動響應(yīng)與飛行試驗(yàn)采集實(shí)際升降舵偏對俯仰角和迎角影響趨勢相近。

    4)在飛機(jī)初步設(shè)計階段,DATCOM軟件具有快速高效的評估飛機(jī)氣動參數(shù)和飛行品質(zhì)的優(yōu)點(diǎn),盡量部分參數(shù)與風(fēng)洞試驗(yàn)和飛行試驗(yàn)所得數(shù)據(jù)有一定偏差,但DATCOM快速估算方法仍有一定價值。

    猜你喜歡
    升降舵風(fēng)洞試驗(yàn)迎角
    無人機(jī)升降舵位移異常分析與改進(jìn)設(shè)計
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    淺談空客A320飛機(jī)升降舵伺服控制邏輯與作動原理
    智富時代(2018年11期)2018-01-15 09:52:06
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    滾轉(zhuǎn)機(jī)動載荷減緩風(fēng)洞試驗(yàn)
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    高速鐵路接觸線覆冰后氣動力特性的風(fēng)洞試驗(yàn)研究
    自拍偷自拍亚洲精品老妇| 中文欧美无线码| 成人亚洲精品一区在线观看 | 久久99精品国语久久久| 伦精品一区二区三区| 男人添女人高潮全过程视频| 91精品国产九色| 欧美xxxx性猛交bbbb| 成年人午夜在线观看视频| 午夜福利视频精品| 亚洲欧洲日产国产| av免费观看日本| 欧美精品人与动牲交sv欧美| 国产av码专区亚洲av| 日韩,欧美,国产一区二区三区| 另类亚洲欧美激情| 一级毛片aaaaaa免费看小| 尤物成人国产欧美一区二区三区| 18禁动态无遮挡网站| 99久久人妻综合| 中文乱码字字幕精品一区二区三区| 欧美zozozo另类| 少妇裸体淫交视频免费看高清| 有码 亚洲区| 黑人高潮一二区| 九九久久精品国产亚洲av麻豆| 亚洲成人中文字幕在线播放| 久久精品人妻少妇| 亚洲欧美一区二区三区黑人 | 国产精品不卡视频一区二区| 一级av片app| 干丝袜人妻中文字幕| 日韩 亚洲 欧美在线| 亚洲精品日本国产第一区| 国产免费视频播放在线视频| 精品少妇久久久久久888优播| 久久久久久国产a免费观看| 亚州av有码| 色吧在线观看| 亚洲av男天堂| 亚洲欧洲国产日韩| 黄片wwwwww| 少妇人妻一区二区三区视频| 一区二区三区精品91| 亚洲精品日韩在线中文字幕| 永久免费av网站大全| 免费观看av网站的网址| 大陆偷拍与自拍| 国产白丝娇喘喷水9色精品| 欧美激情久久久久久爽电影| 色5月婷婷丁香| 国产精品人妻久久久久久| 97超碰精品成人国产| av在线亚洲专区| 日韩成人伦理影院| 亚洲性久久影院| www.av在线官网国产| 99久久九九国产精品国产免费| 免费大片黄手机在线观看| 内射极品少妇av片p| 国产伦理片在线播放av一区| 亚洲欧美成人精品一区二区| 日韩成人伦理影院| 麻豆成人av视频| 尾随美女入室| 中国美白少妇内射xxxbb| 国产精品不卡视频一区二区| 中国美白少妇内射xxxbb| 国产精品不卡视频一区二区| 国产成人一区二区在线| 国产毛片在线视频| 国产毛片在线视频| 国产精品99久久久久久久久| 日韩国内少妇激情av| 伊人久久精品亚洲午夜| 久久精品国产亚洲av涩爱| 免费大片黄手机在线观看| 看十八女毛片水多多多| 五月开心婷婷网| 国产白丝娇喘喷水9色精品| 五月开心婷婷网| 人人妻人人澡人人爽人人夜夜| 精品酒店卫生间| 伦理电影大哥的女人| 一个人看的www免费观看视频| 夫妻性生交免费视频一级片| h日本视频在线播放| 三级经典国产精品| 老司机影院毛片| 日本-黄色视频高清免费观看| 麻豆久久精品国产亚洲av| 男女啪啪激烈高潮av片| 2021天堂中文幕一二区在线观| 一区二区三区精品91| 少妇被粗大猛烈的视频| 九九久久精品国产亚洲av麻豆| 亚洲一级一片aⅴ在线观看| 国产国拍精品亚洲av在线观看| 亚洲内射少妇av| 国产黄a三级三级三级人| 一级a做视频免费观看| 国产国拍精品亚洲av在线观看| 最近最新中文字幕免费大全7| 国产精品一及| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜日本视频在线| 国产黄频视频在线观看| av国产免费在线观看| 欧美日韩精品成人综合77777| 婷婷色综合大香蕉| 免费高清在线观看视频在线观看| 哪个播放器可以免费观看大片| 最新中文字幕久久久久| 久久久成人免费电影| 18禁裸乳无遮挡动漫免费视频 | 18禁动态无遮挡网站| 蜜桃久久精品国产亚洲av| av在线app专区| 国产免费又黄又爽又色| 成人一区二区视频在线观看| 少妇的逼水好多| 91狼人影院| 久久久久国产网址| 99精国产麻豆久久婷婷| 夜夜看夜夜爽夜夜摸| 两个人的视频大全免费| 日韩一本色道免费dvd| 久久99精品国语久久久| 国产一区二区在线观看日韩| 大又大粗又爽又黄少妇毛片口| kizo精华| 国产毛片在线视频| www.av在线官网国产| 在线a可以看的网站| 国内精品美女久久久久久| 国内少妇人妻偷人精品xxx网站| 日韩人妻高清精品专区| 国产一区二区三区av在线| 亚洲欧美精品专区久久| 3wmmmm亚洲av在线观看| 草草在线视频免费看| 成人国产av品久久久| 各种免费的搞黄视频| 亚洲精品色激情综合| 午夜激情久久久久久久| 一级a做视频免费观看| 亚洲最大成人av| 黄色怎么调成土黄色| 亚洲,一卡二卡三卡| 能在线免费看毛片的网站| 亚洲精品日本国产第一区| 国产69精品久久久久777片| 午夜精品国产一区二区电影 | 人人妻人人澡人人爽人人夜夜| 欧美三级亚洲精品| 日韩人妻高清精品专区| av国产免费在线观看| 免费不卡的大黄色大毛片视频在线观看| 日韩av免费高清视频| 97人妻精品一区二区三区麻豆| 哪个播放器可以免费观看大片| 亚洲av成人精品一二三区| 日本午夜av视频| 国产爱豆传媒在线观看| 国产男女超爽视频在线观看| 白带黄色成豆腐渣| 全区人妻精品视频| 人妻夜夜爽99麻豆av| 91午夜精品亚洲一区二区三区| 自拍偷自拍亚洲精品老妇| 日日撸夜夜添| 男女国产视频网站| 黄色怎么调成土黄色| 高清欧美精品videossex| 久久久久性生活片| 欧美极品一区二区三区四区| 最近最新中文字幕大全电影3| a级一级毛片免费在线观看| 亚洲,欧美,日韩| 成年版毛片免费区| 欧美日韩一区二区视频在线观看视频在线 | 又黄又爽又刺激的免费视频.| 一级毛片aaaaaa免费看小| 国产精品不卡视频一区二区| www.av在线官网国产| 国产熟女欧美一区二区| 交换朋友夫妻互换小说| 国产极品天堂在线| 美女主播在线视频| 国产成人91sexporn| 自拍欧美九色日韩亚洲蝌蚪91 | 中文在线观看免费www的网站| 在线观看一区二区三区激情| 夜夜爽夜夜爽视频| 国产探花极品一区二区| 亚洲经典国产精华液单| 两个人的视频大全免费| 黄片无遮挡物在线观看| 免费av毛片视频| 少妇的逼好多水| 中国国产av一级| 校园人妻丝袜中文字幕| 91精品伊人久久大香线蕉| 欧美xxxx黑人xx丫x性爽| 老司机影院成人| 99热全是精品| 午夜精品国产一区二区电影 | 少妇高潮的动态图| 天堂中文最新版在线下载 | 日本wwww免费看| 久久久国产一区二区| 激情 狠狠 欧美| 色播亚洲综合网| 乱系列少妇在线播放| 国产欧美日韩一区二区三区在线 | 国产91av在线免费观看| 夜夜爽夜夜爽视频| 91精品一卡2卡3卡4卡| 一个人看的www免费观看视频| 男的添女的下面高潮视频| 街头女战士在线观看网站| 久久99热这里只有精品18| 中文精品一卡2卡3卡4更新| 久久久精品免费免费高清| 一个人看视频在线观看www免费| 亚洲电影在线观看av| 国产极品天堂在线| 国产成人91sexporn| 欧美性感艳星| 男人和女人高潮做爰伦理| 建设人人有责人人尽责人人享有的 | 在线观看av片永久免费下载| 精品一区二区免费观看| 亚洲精品国产成人久久av| 一级毛片我不卡| a级毛色黄片| 精品一区二区三卡| 久久精品国产亚洲av涩爱| 午夜福利网站1000一区二区三区| 国产精品无大码| 免费av毛片视频| 日韩强制内射视频| 国精品久久久久久国模美| 最近的中文字幕免费完整| 欧美成人精品欧美一级黄| 成人美女网站在线观看视频| 黄色一级大片看看| 交换朋友夫妻互换小说| 3wmmmm亚洲av在线观看| xxx大片免费视频| 国产免费视频播放在线视频| 简卡轻食公司| 97在线人人人人妻| av在线亚洲专区| 国产熟女欧美一区二区| 精品久久久久久久久亚洲| 五月玫瑰六月丁香| 最近中文字幕高清免费大全6| 亚洲国产色片| 国产久久久一区二区三区| 在线播放无遮挡| 在线天堂最新版资源| 精华霜和精华液先用哪个| xxx大片免费视频| 91精品伊人久久大香线蕉| 国产中年淑女户外野战色| 老女人水多毛片| 国产一区二区三区综合在线观看 | 性插视频无遮挡在线免费观看| 最近2019中文字幕mv第一页| 亚洲精品乱码久久久v下载方式| 国产精品国产av在线观看| 亚洲成人av在线免费| 日本色播在线视频| 九色成人免费人妻av| 一区二区三区免费毛片| 美女视频免费永久观看网站| 精品视频人人做人人爽| 最近中文字幕2019免费版| 亚洲成人精品中文字幕电影| 少妇的逼好多水| 国产免费视频播放在线视频| 人妻系列 视频| 国产大屁股一区二区在线视频| 久久久久国产网址| 久久精品人妻少妇| 国产一区亚洲一区在线观看| 久久精品久久精品一区二区三区| 男人爽女人下面视频在线观看| 高清毛片免费看| 国产成人freesex在线| 永久免费av网站大全| www.色视频.com| 欧美zozozo另类| 建设人人有责人人尽责人人享有的 | 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美在线一区| 亚洲av一区综合| 午夜激情福利司机影院| 一级爰片在线观看| 亚洲av二区三区四区| 少妇人妻一区二区三区视频| 久久久久久国产a免费观看| 插阴视频在线观看视频| 91精品国产九色| 国产毛片在线视频| 久久久久久九九精品二区国产| 爱豆传媒免费全集在线观看| 成人无遮挡网站| 嘟嘟电影网在线观看| 国产又色又爽无遮挡免| 制服丝袜香蕉在线| 国产精品av视频在线免费观看| 成人漫画全彩无遮挡| 午夜爱爱视频在线播放| 最近最新中文字幕大全电影3| 国产毛片a区久久久久| 亚洲真实伦在线观看| 亚洲国产欧美人成| 日本一二三区视频观看| 久久久久久久久久久丰满| 在线播放无遮挡| 日韩一本色道免费dvd| 97超视频在线观看视频| 久久午夜福利片| 在线亚洲精品国产二区图片欧美 | 性色av一级| 新久久久久国产一级毛片| 亚洲av二区三区四区| 黄片无遮挡物在线观看| 99九九线精品视频在线观看视频| 国产精品熟女久久久久浪| 国产成人a区在线观看| 欧美日韩综合久久久久久| 一区二区av电影网| 观看美女的网站| 秋霞在线观看毛片| kizo精华| 久久99蜜桃精品久久| 国产精品爽爽va在线观看网站| 男人爽女人下面视频在线观看| 夜夜爽夜夜爽视频| 永久免费av网站大全| 爱豆传媒免费全集在线观看| 日韩免费高清中文字幕av| 免费黄网站久久成人精品| 天天一区二区日本电影三级| 国产精品成人在线| 亚洲av电影在线观看一区二区三区 | 亚洲精品中文字幕在线视频 | 国产免费一级a男人的天堂| 爱豆传媒免费全集在线观看| 亚洲婷婷狠狠爱综合网| 国产久久久一区二区三区| 免费av不卡在线播放| 99久久精品热视频| 一个人观看的视频www高清免费观看| 亚洲精品国产av蜜桃| 白带黄色成豆腐渣| av国产免费在线观看| 观看美女的网站| 国产欧美日韩一区二区三区在线 | 亚洲国产精品999| 国产高潮美女av| 又粗又硬又长又爽又黄的视频| 在线 av 中文字幕| 男女那种视频在线观看| 亚洲av日韩在线播放| 又粗又硬又长又爽又黄的视频| 亚洲精品日本国产第一区| 国产日韩欧美在线精品| 午夜福利网站1000一区二区三区| 2018国产大陆天天弄谢| av女优亚洲男人天堂| 午夜视频国产福利| 汤姆久久久久久久影院中文字幕| 精品久久国产蜜桃| 99热6这里只有精品| 97超碰精品成人国产| 亚洲自拍偷在线| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 三级经典国产精品| 国产黄频视频在线观看| 五月玫瑰六月丁香| 久久国内精品自在自线图片| 亚洲在线观看片| 老司机影院毛片| 边亲边吃奶的免费视频| 自拍欧美九色日韩亚洲蝌蚪91 | 尤物成人国产欧美一区二区三区| .国产精品久久| 99久久中文字幕三级久久日本| 男人和女人高潮做爰伦理| 联通29元200g的流量卡| 91午夜精品亚洲一区二区三区| 噜噜噜噜噜久久久久久91| 日韩一区二区三区影片| 久久久久久久久久成人| 日本午夜av视频| 在线观看免费高清a一片| 欧美高清性xxxxhd video| 国产老妇女一区| 五月玫瑰六月丁香| 亚洲国产色片| 亚洲精华国产精华液的使用体验| 精品人妻视频免费看| 69av精品久久久久久| 亚洲图色成人| 国产精品久久久久久久电影| 国产亚洲91精品色在线| 亚洲欧美日韩另类电影网站 | 性色avwww在线观看| 国产成人精品一,二区| 国产乱人视频| 美女xxoo啪啪120秒动态图| 中文字幕制服av| av福利片在线观看| 国产欧美另类精品又又久久亚洲欧美| 成年人午夜在线观看视频| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 欧美xxxx黑人xx丫x性爽| tube8黄色片| 国产黄色免费在线视频| 久久综合国产亚洲精品| 在线亚洲精品国产二区图片欧美 | 亚洲丝袜综合中文字幕| 2022亚洲国产成人精品| 亚洲欧美日韩卡通动漫| 久久久久久久大尺度免费视频| av福利片在线观看| 亚洲内射少妇av| 成年免费大片在线观看| 网址你懂的国产日韩在线| 国产精品国产三级国产av玫瑰| 成年人午夜在线观看视频| 91狼人影院| www.色视频.com| 一本一本综合久久| 国产伦精品一区二区三区四那| 一区二区三区精品91| 在线观看免费高清a一片| tube8黄色片| 男女边摸边吃奶| 免费大片18禁| 香蕉精品网在线| 国产精品一区二区三区四区免费观看| 亚洲av男天堂| 99热6这里只有精品| 亚洲国产日韩一区二区| 成人国产av品久久久| 亚洲av男天堂| 久久99热6这里只有精品| 免费观看av网站的网址| 在线观看国产h片| 女的被弄到高潮叫床怎么办| 三级男女做爰猛烈吃奶摸视频| 七月丁香在线播放| 亚洲成人久久爱视频| 人妻一区二区av| 免费播放大片免费观看视频在线观看| 伊人久久精品亚洲午夜| 日韩欧美精品v在线| 大片免费播放器 马上看| 丝袜喷水一区| 自拍偷自拍亚洲精品老妇| 免费观看的影片在线观看| 91狼人影院| 午夜激情久久久久久久| 内地一区二区视频在线| 纵有疾风起免费观看全集完整版| 国产综合精华液| 美女被艹到高潮喷水动态| 草草在线视频免费看| 人妻系列 视频| 欧美日韩综合久久久久久| 黄色怎么调成土黄色| 久久久久久久亚洲中文字幕| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 国产精品.久久久| av一本久久久久| 精品人妻一区二区三区麻豆| 只有这里有精品99| 99久国产av精品国产电影| 国产综合懂色| 纵有疾风起免费观看全集完整版| 亚洲国产精品成人综合色| 欧美老熟妇乱子伦牲交| 国产成年人精品一区二区| 精品99又大又爽又粗少妇毛片| 久久精品人妻少妇| 婷婷色综合www| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 亚洲人成网站高清观看| 九九久久精品国产亚洲av麻豆| 国产精品伦人一区二区| 日韩欧美 国产精品| 国产在视频线精品| 国产欧美亚洲国产| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 精品熟女少妇av免费看| 各种免费的搞黄视频| 精品一区在线观看国产| 免费看av在线观看网站| 国产亚洲av嫩草精品影院| 99热这里只有是精品50| 黄色配什么色好看| 一级片'在线观看视频| 国产淫片久久久久久久久| 国产淫语在线视频| 亚洲精品日韩在线中文字幕| 综合色丁香网| 成人二区视频| www.av在线官网国产| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 亚洲,欧美,日韩| 日韩伦理黄色片| 毛片女人毛片| 色综合色国产| 免费观看av网站的网址| 国产亚洲91精品色在线| 成年av动漫网址| 国产成人精品福利久久| 日韩制服骚丝袜av| 丰满乱子伦码专区| 一级爰片在线观看| 人妻夜夜爽99麻豆av| 激情 狠狠 欧美| 欧美极品一区二区三区四区| 2022亚洲国产成人精品| 少妇裸体淫交视频免费看高清| 特大巨黑吊av在线直播| 国产精品.久久久| 久久精品久久久久久噜噜老黄| 69av精品久久久久久| 免费黄频网站在线观看国产| 欧美激情国产日韩精品一区| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说 | 午夜视频国产福利| 欧美bdsm另类| 国产黄色视频一区二区在线观看| 少妇的逼水好多| 亚洲自拍偷在线| 久久99热这里只频精品6学生| 国产人妻一区二区三区在| 国产精品不卡视频一区二区| 久久久久精品久久久久真实原创| 最后的刺客免费高清国语| 美女被艹到高潮喷水动态| 菩萨蛮人人尽说江南好唐韦庄| 在现免费观看毛片| 夜夜爽夜夜爽视频| 国产黄色免费在线视频| 亚洲av电影在线观看一区二区三区 | 国产男人的电影天堂91| 久久久色成人| 久久久欧美国产精品| 亚洲自偷自拍三级| 亚洲成人久久爱视频| 亚洲欧美成人精品一区二区| 又大又黄又爽视频免费| 亚洲va在线va天堂va国产| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级专区第一集| 校园人妻丝袜中文字幕| 亚洲欧美日韩另类电影网站 | 国产一区二区三区av在线| 少妇猛男粗大的猛烈进出视频 | 欧美激情国产日韩精品一区| 亚洲成人精品中文字幕电影| 一区二区三区乱码不卡18| 午夜精品国产一区二区电影 | 午夜视频国产福利| 深夜a级毛片| a级毛色黄片| 黄色配什么色好看| 777米奇影视久久| 秋霞在线观看毛片| 亚洲成色77777| 听说在线观看完整版免费高清| 日日啪夜夜爽| 国产精品无大码| 欧美性感艳星| 亚洲av中文字字幕乱码综合| 亚洲精品日韩av片在线观看| 亚洲精品影视一区二区三区av| 日本一本二区三区精品| 国产高潮美女av| 舔av片在线| 男人狂女人下面高潮的视频| 久久女婷五月综合色啪小说 | 色网站视频免费| 婷婷色综合大香蕉| 国产黄片视频在线免费观看| 99久久精品一区二区三区| 色哟哟·www| 99热这里只有精品一区| 在线精品无人区一区二区三 | 国产免费又黄又爽又色| 97人妻精品一区二区三区麻豆| 日韩一区二区三区影片| 国产熟女欧美一区二区| 97在线视频观看| 内地一区二区视频在线| 成人国产麻豆网| 人妻少妇偷人精品九色| 久久久久久久久久人人人人人人| 观看美女的网站| 99热这里只有精品一区|