• <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)研究
    久久久久久人妻| 亚洲精品第二区| 免费观看无遮挡的男女| 爱豆传媒免费全集在线观看| 91午夜精品亚洲一区二区三区| h视频一区二区三区| 亚洲国产成人一精品久久久| 亚洲精品日本国产第一区| 亚洲婷婷狠狠爱综合网| 欧美另类一区| 美女国产视频在线观看| 亚洲精品,欧美精品| videosex国产| 欧美老熟妇乱子伦牲交| 亚洲欧美色中文字幕在线| 日韩一区二区三区影片| 在现免费观看毛片| 国产 一区精品| 国产成人精品在线电影| 18+在线观看网站| 久久久久久人人人人人| a 毛片基地| 一级毛片我不卡| 国精品久久久久久国模美| 亚洲精品国产色婷婷电影| 国产成人免费观看mmmm| 国产精品无大码| 欧美精品一区二区免费开放| 日韩一本色道免费dvd| 高清黄色对白视频在线免费看| 大香蕉久久网| 人人妻人人爽人人添夜夜欢视频| 国产av一区二区精品久久| 亚洲国产看品久久| 久久久国产欧美日韩av| 人妻系列 视频| 国产精品久久久久久精品电影小说| 9191精品国产免费久久| av一本久久久久| 精品人妻偷拍中文字幕| 亚洲,欧美精品.| 看免费av毛片| 精品国产一区二区久久| 久久久精品免费免费高清| 久久久久久人人人人人| 大码成人一级视频| 免费高清在线观看视频在线观看| 波野结衣二区三区在线| 亚洲中文av在线| 亚洲中文av在线| 国产伦理片在线播放av一区| 亚洲精华国产精华液的使用体验| 乱码一卡2卡4卡精品| 久久久久久久大尺度免费视频| 97精品久久久久久久久久精品| 春色校园在线视频观看| 黄色毛片三级朝国网站| 免费观看无遮挡的男女| 精品少妇久久久久久888优播| 久久韩国三级中文字幕| 9色porny在线观看| 99久久中文字幕三级久久日本| 精品国产一区二区久久| 亚洲欧洲国产日韩| av福利片在线| 99热国产这里只有精品6| 欧美亚洲 丝袜 人妻 在线| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区黑人 | 久久久久久久大尺度免费视频| 黄色 视频免费看| 婷婷成人精品国产| 天堂俺去俺来也www色官网| 久久精品夜色国产| 曰老女人黄片| 一级黄片播放器| 在线免费观看不下载黄p国产| 久久久久精品性色| a 毛片基地| 少妇熟女欧美另类| 亚洲精品456在线播放app| 久久精品久久久久久久性| 久久人人爽人人片av| 久久99热6这里只有精品| 亚洲精品第二区| 黄色视频在线播放观看不卡| 久久国内精品自在自线图片| 又黄又粗又硬又大视频| 性色avwww在线观看| 一二三四中文在线观看免费高清| 久久综合国产亚洲精品| 在线观看美女被高潮喷水网站| av在线播放精品| 日日摸夜夜添夜夜爱| 王馨瑶露胸无遮挡在线观看| 国产亚洲精品久久久com| 国产亚洲精品久久久com| av女优亚洲男人天堂| 男人爽女人下面视频在线观看| 午夜激情av网站| 一级黄片播放器| 婷婷色av中文字幕| 欧美日韩亚洲高清精品| 纯流量卡能插随身wifi吗| 精品亚洲乱码少妇综合久久| 欧美精品高潮呻吟av久久| 日韩视频在线欧美| 午夜福利乱码中文字幕| 1024视频免费在线观看| 视频在线观看一区二区三区| 久久99一区二区三区| 1024视频免费在线观看| 欧美日韩成人在线一区二区| av播播在线观看一区| av福利片在线| 极品人妻少妇av视频| 亚洲av.av天堂| 精品久久蜜臀av无| 色婷婷av一区二区三区视频| 桃花免费在线播放| 亚洲综合色网址| 国产片内射在线| xxxhd国产人妻xxx| 王馨瑶露胸无遮挡在线观看| 99久久精品国产国产毛片| 美女国产高潮福利片在线看| 日韩 亚洲 欧美在线| 亚洲精品456在线播放app| 十八禁高潮呻吟视频| 国产xxxxx性猛交| 日韩在线高清观看一区二区三区| 91精品三级在线观看| 国产精品一区二区在线观看99| 在线天堂最新版资源| 母亲3免费完整高清在线观看 | 在线观看免费视频网站a站| 亚洲人成网站在线观看播放| 人人妻人人爽人人添夜夜欢视频| 五月开心婷婷网| 国产精品久久久久久久电影| 这个男人来自地球电影免费观看 | 女的被弄到高潮叫床怎么办| 极品少妇高潮喷水抽搐| 久久精品国产鲁丝片午夜精品| 亚洲第一av免费看| 精品人妻在线不人妻| 在线精品无人区一区二区三| 精品少妇黑人巨大在线播放| 亚洲第一区二区三区不卡| a级毛片在线看网站| 最黄视频免费看| 啦啦啦在线观看免费高清www| 在线观看国产h片| 久久韩国三级中文字幕| 九色亚洲精品在线播放| 黄网站色视频无遮挡免费观看| 高清毛片免费看| 久久精品熟女亚洲av麻豆精品| 岛国毛片在线播放| 岛国毛片在线播放| 久久女婷五月综合色啪小说| 九草在线视频观看| 欧美最新免费一区二区三区| 熟女电影av网| 亚洲av国产av综合av卡| av在线观看视频网站免费| 亚洲国产av影院在线观看| 最近最新中文字幕免费大全7| 国产男女内射视频| 亚洲欧美日韩卡通动漫| 日本欧美国产在线视频| 人妻系列 视频| av线在线观看网站| 97在线人人人人妻| 亚洲三级黄色毛片| 国产毛片在线视频| 90打野战视频偷拍视频| 中文字幕另类日韩欧美亚洲嫩草| 婷婷色麻豆天堂久久| 中文字幕人妻熟女乱码| 亚洲av.av天堂| 中国国产av一级| 性色avwww在线观看| 国产精品一二三区在线看| 日日啪夜夜爽| 久久久欧美国产精品| 亚洲av电影在线进入| 美女大奶头黄色视频| 日韩中字成人| 久久人人爽人人爽人人片va| 国产综合精华液| 人妻人人澡人人爽人人| 欧美xxxx性猛交bbbb| 18禁观看日本| 亚洲国产看品久久| 九草在线视频观看| 亚洲欧美一区二区三区黑人 | 五月开心婷婷网| 亚洲一码二码三码区别大吗| 午夜av观看不卡| 亚洲av.av天堂| 久久久久久久精品精品| 2021少妇久久久久久久久久久| 久久久精品区二区三区| 精品一区二区三区视频在线| 亚洲,一卡二卡三卡| a级毛片黄视频| 亚洲第一区二区三区不卡| 久久国产精品大桥未久av| 91国产中文字幕| 国产亚洲欧美精品永久| 黄色怎么调成土黄色| 最近手机中文字幕大全| 超碰97精品在线观看| 国产色婷婷99| 女人久久www免费人成看片| 18在线观看网站| 亚洲av.av天堂| 欧美日韩视频高清一区二区三区二| 精品卡一卡二卡四卡免费| 久久久精品区二区三区| 免费黄网站久久成人精品| 久热这里只有精品99| 插逼视频在线观看| 五月伊人婷婷丁香| 18禁裸乳无遮挡动漫免费视频| 精品一区在线观看国产| 午夜免费鲁丝| 免费日韩欧美在线观看| 高清毛片免费看| 国产av一区二区精品久久| 大话2 男鬼变身卡| 99久久中文字幕三级久久日本| av国产精品久久久久影院| 国产亚洲精品久久久com| 中文字幕亚洲精品专区| 精品国产一区二区三区四区第35| 黄色一级大片看看| 婷婷色麻豆天堂久久| 在线观看三级黄色| av卡一久久| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 大陆偷拍与自拍| 成人手机av| 久久久精品94久久精品| 人妻人人澡人人爽人人| 亚洲精品国产av蜜桃| 国产伦理片在线播放av一区| 中文天堂在线官网| 国产精品国产三级国产av玫瑰| 久久精品久久精品一区二区三区| 黄色视频在线播放观看不卡| 乱码一卡2卡4卡精品| 在线天堂最新版资源| 菩萨蛮人人尽说江南好唐韦庄| 黑人欧美特级aaaaaa片| 少妇 在线观看| h视频一区二区三区| 人人妻人人澡人人爽人人夜夜| 欧美人与性动交α欧美精品济南到 | 一本大道久久a久久精品| 日韩一区二区三区影片| 久久婷婷青草| 2022亚洲国产成人精品| 亚洲精品色激情综合| 嫩草影院入口| 国产精品蜜桃在线观看| 在线观看免费视频网站a站| 国产av码专区亚洲av| 日本wwww免费看| 久久精品久久精品一区二区三区| 丰满饥渴人妻一区二区三| 成人综合一区亚洲| 日本-黄色视频高清免费观看| 精品国产一区二区久久| 午夜福利视频精品| 女人被躁到高潮嗷嗷叫费观| 欧美+日韩+精品| 丝袜人妻中文字幕| 大片电影免费在线观看免费| 2021少妇久久久久久久久久久| av黄色大香蕉| videosex国产| 免费播放大片免费观看视频在线观看| 日本av手机在线免费观看| 国产乱来视频区| 日韩欧美一区视频在线观看| 亚洲成人av在线免费| 亚洲欧洲国产日韩| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 国产视频首页在线观看| a级毛色黄片| 亚洲精华国产精华液的使用体验| 久久久久久久久久久久大奶| 亚洲内射少妇av| 水蜜桃什么品种好| av视频免费观看在线观看| 精品人妻偷拍中文字幕| 国产激情久久老熟女| 国产精品久久久久久精品电影小说| 91成人精品电影| 日本黄大片高清| 日韩av免费高清视频| 男人爽女人下面视频在线观看| 男人添女人高潮全过程视频| 亚洲欧美一区二区三区黑人 | av福利片在线| 好男人视频免费观看在线| 精品亚洲乱码少妇综合久久| 大陆偷拍与自拍| 我的女老师完整版在线观看| 丰满少妇做爰视频| 亚洲欧美色中文字幕在线| 亚洲一码二码三码区别大吗| 中文精品一卡2卡3卡4更新| 美女国产高潮福利片在线看| 日韩一区二区视频免费看| 日韩av不卡免费在线播放| 色哟哟·www| 秋霞在线观看毛片| 国产伦理片在线播放av一区| 亚洲熟女精品中文字幕| 啦啦啦视频在线资源免费观看| 国产精品蜜桃在线观看| 国产乱来视频区| 欧美日本中文国产一区发布| 一区二区av电影网| 亚洲综合精品二区| 久久久a久久爽久久v久久| 久久99一区二区三区| 高清av免费在线| 亚洲五月色婷婷综合| 久久久久久久久久成人| 少妇被粗大猛烈的视频| 一级毛片我不卡| 精品亚洲成国产av| 亚洲国产精品一区二区三区在线| 色哟哟·www| 麻豆精品久久久久久蜜桃| 不卡视频在线观看欧美| 国产一区二区在线观看av| 一区二区三区乱码不卡18| 成人综合一区亚洲| 最新的欧美精品一区二区| 免费观看性生交大片5| 99久久人妻综合| 美女脱内裤让男人舔精品视频| 少妇的丰满在线观看| 在现免费观看毛片| 欧美老熟妇乱子伦牲交| 一区在线观看完整版| 色视频在线一区二区三区| 女人被躁到高潮嗷嗷叫费观| 国产一区二区三区综合在线观看 | 午夜激情av网站| 午夜免费观看性视频| 日本av免费视频播放| av国产久精品久网站免费入址| 久久久久久久亚洲中文字幕| 午夜日本视频在线| 校园人妻丝袜中文字幕| 蜜桃国产av成人99| 不卡视频在线观看欧美| 桃花免费在线播放| 中文乱码字字幕精品一区二区三区| 97超碰精品成人国产| 又大又黄又爽视频免费| 亚洲精品国产av成人精品| h视频一区二区三区| 亚洲精品,欧美精品| av卡一久久| 亚洲三级黄色毛片| 超色免费av| 免费观看性生交大片5| 国产深夜福利视频在线观看| 最近最新中文字幕大全免费视频 | 综合色丁香网| 久久人人97超碰香蕉20202| 日本黄大片高清| 成人漫画全彩无遮挡| 久久精品国产鲁丝片午夜精品| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 国产日韩欧美亚洲二区| 国产黄色免费在线视频| 国产成人a∨麻豆精品| 日韩视频在线欧美| 国产精品无大码| 亚洲欧美清纯卡通| a级片在线免费高清观看视频| 两性夫妻黄色片 | 在线免费观看不下载黄p国产| 色94色欧美一区二区| av天堂久久9| 欧美精品亚洲一区二区| 久久久久久人人人人人| 看非洲黑人一级黄片| 国产一区二区激情短视频 | 一级毛片 在线播放| 美女内射精品一级片tv| 国产日韩欧美亚洲二区| 国产精品国产三级专区第一集| 99久久综合免费| 免费人成在线观看视频色| av福利片在线| 女性被躁到高潮视频| 亚洲在久久综合| 精品久久蜜臀av无| 啦啦啦啦在线视频资源| 国产精品.久久久| 日韩三级伦理在线观看| 欧美xxxx性猛交bbbb| 日韩,欧美,国产一区二区三区| 精品久久久久久电影网| 99热6这里只有精品| 欧美日韩综合久久久久久| 久久ye,这里只有精品| 97在线人人人人妻| 男女啪啪激烈高潮av片| 亚洲高清免费不卡视频| 日韩欧美一区视频在线观看| 视频区图区小说| a 毛片基地| 少妇高潮的动态图| 777米奇影视久久| 久久久精品区二区三区| 国产成人精品久久久久久| 欧美日韩视频高清一区二区三区二| 久久久a久久爽久久v久久| 午夜av观看不卡| 中文字幕精品免费在线观看视频 | 99久久中文字幕三级久久日本| 麻豆乱淫一区二区| 美女脱内裤让男人舔精品视频| 国产日韩一区二区三区精品不卡| 亚洲国产毛片av蜜桃av| av在线观看视频网站免费| 一区二区三区乱码不卡18| 深夜精品福利| 大码成人一级视频| 女的被弄到高潮叫床怎么办| 国产福利在线免费观看视频| 热re99久久国产66热| 宅男免费午夜| 成人国产av品久久久| 在线观看www视频免费| 亚洲欧美一区二区三区黑人 | 久久人人爽av亚洲精品天堂| 亚洲国产精品国产精品| 欧美亚洲 丝袜 人妻 在线| 亚洲国产最新在线播放| 国产成人一区二区在线| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久久久免| 午夜免费男女啪啪视频观看| 久久久久久久久久久久大奶| 午夜福利视频精品| kizo精华| 国产男人的电影天堂91| 国产精品久久久久久久久免| 国产视频首页在线观看| 欧美日韩av久久| 人妻一区二区av| 国产亚洲欧美精品永久| 老熟女久久久| 日本av免费视频播放| 水蜜桃什么品种好| 美女脱内裤让男人舔精品视频| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美色中文字幕在线| 亚洲美女黄色视频免费看| 国产精品一区www在线观看| 国产有黄有色有爽视频| 亚洲 欧美一区二区三区| 成人午夜精彩视频在线观看| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放| 视频在线观看一区二区三区| 午夜91福利影院| 久久鲁丝午夜福利片| 午夜激情久久久久久久| 亚洲少妇的诱惑av| 99热国产这里只有精品6| 一边摸一边做爽爽视频免费| 你懂的网址亚洲精品在线观看| 一区二区三区四区激情视频| 熟女av电影| 亚洲内射少妇av| 国产高清三级在线| 国产精品秋霞免费鲁丝片| 另类精品久久| 成人毛片60女人毛片免费| av国产久精品久网站免费入址| 国产毛片在线视频| 黑人猛操日本美女一级片| 亚洲国产毛片av蜜桃av| 亚洲欧洲国产日韩| 天堂俺去俺来也www色官网| 亚洲av日韩在线播放| 一本久久精品| 老司机影院毛片| 国产精品久久久久久av不卡| 久久毛片免费看一区二区三区| 亚洲综合色惰| 熟妇人妻不卡中文字幕| av播播在线观看一区| 九色亚洲精品在线播放| 亚洲欧美精品自产自拍| 美女内射精品一级片tv| 久久精品国产亚洲av涩爱| 亚洲欧洲精品一区二区精品久久久 | 成人无遮挡网站| 欧美老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 国产一级毛片在线| av播播在线观看一区| 国产 精品1| 精品国产露脸久久av麻豆| 卡戴珊不雅视频在线播放| 久久久久精品人妻al黑| 看免费av毛片| 亚洲性久久影院| 少妇的逼好多水| 亚洲第一av免费看| 99热网站在线观看| 久久久久久久大尺度免费视频| 日产精品乱码卡一卡2卡三| 亚洲精品第二区| 亚洲av欧美aⅴ国产| 国产在视频线精品| 亚洲中文av在线| 在线观看美女被高潮喷水网站| 免费看不卡的av| 国产亚洲午夜精品一区二区久久| 老女人水多毛片| 久久久久精品久久久久真实原创| 如何舔出高潮| 久久99蜜桃精品久久| 热re99久久精品国产66热6| 久久久久久久久久久久大奶| 国产精品熟女久久久久浪| 最近最新中文字幕大全免费视频 | 欧美人与性动交α欧美软件 | 黄色怎么调成土黄色| 久久久欧美国产精品| 在线观看三级黄色| 亚洲欧美色中文字幕在线| 国产一区二区在线观看av| 自线自在国产av| 免费在线观看黄色视频的| 亚洲国产成人一精品久久久| 日本午夜av视频| 青春草国产在线视频| 国产精品国产三级专区第一集| 麻豆乱淫一区二区| 国产黄色免费在线视频| 97精品久久久久久久久久精品| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 久久国产精品大桥未久av| 亚洲欧美日韩卡通动漫| 日韩,欧美,国产一区二区三区| 黄色怎么调成土黄色| 国产在线视频一区二区| 啦啦啦视频在线资源免费观看| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 少妇高潮的动态图| 一二三四中文在线观看免费高清| 久久久久国产网址| 国产av国产精品国产| 国产精品.久久久| 国产成人a∨麻豆精品| 日韩熟女老妇一区二区性免费视频| 国精品久久久久久国模美| 免费大片黄手机在线观看| 国产成人a∨麻豆精品| 国产探花极品一区二区| 国产综合精华液| 精品少妇黑人巨大在线播放| 最近中文字幕高清免费大全6| 欧美日本中文国产一区发布| 黄色毛片三级朝国网站| 伦理电影大哥的女人| 午夜福利视频精品| 国产精品一区二区在线观看99| 午夜福利影视在线免费观看| av.在线天堂| 国产精品嫩草影院av在线观看| 成人毛片a级毛片在线播放| 欧美xxxx性猛交bbbb| 久久99精品国语久久久| 黄色 视频免费看| 26uuu在线亚洲综合色| 亚洲伊人色综图| 各种免费的搞黄视频| 久久这里有精品视频免费| 男女国产视频网站| 久久精品国产亚洲av涩爱| 久久韩国三级中文字幕| 国产av精品麻豆| av不卡在线播放| 国产永久视频网站| 建设人人有责人人尽责人人享有的| 免费观看a级毛片全部| 久久精品国产鲁丝片午夜精品| 各种免费的搞黄视频| 久久精品aⅴ一区二区三区四区 | 91aial.com中文字幕在线观看|