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

    基于速度分裂法的翼型陣風(fēng)響應(yīng)及減緩數(shù)值模擬

    2023-05-09 08:42:50郭同慶紀(jì)哲翰陸志良
    關(guān)鍵詞:陣風(fēng)升力黏性

    高 坤,郭同慶,紀(jì)哲翰,周 迪,陸志良

    (南京航空航天大學(xué) 非定??諝鈩?dòng)力學(xué)與流動(dòng)控制工業(yè)和信息化部重點(diǎn)實(shí)驗(yàn)室,南京 210016)

    0 引 言

    陣風(fēng)載荷響應(yīng)作為一種氣動(dòng)彈性動(dòng)力響應(yīng)問題,是指飛行器在大氣中飛行時(shí)遇到突然的陣風(fēng)擾動(dòng),這一擾動(dòng)氣流將改變飛機(jī)氣動(dòng)升力面的有效迎角,引起升力等的突然改變,使得飛行器產(chǎn)生涉及彈性變形的陣風(fēng)載荷響應(yīng)。當(dāng)飛行器處于陣風(fēng)環(huán)境中時(shí),飛行器結(jié)構(gòu)將承載較大的陣風(fēng)載荷,這將導(dǎo)致飛行器操縱困難且容易發(fā)生結(jié)構(gòu)疲勞損壞。隨著航空科學(xué)技術(shù)的發(fā)展,大型飛行器的結(jié)構(gòu)設(shè)計(jì)朝著低重量和大柔性方向發(fā)展,這導(dǎo)致陣風(fēng)對大型飛行器的飛行影響更為復(fù)雜[1]。因此,對于現(xiàn)代飛行器設(shè)計(jì),其陣風(fēng)響應(yīng)分析的非定常氣動(dòng)模型由傳統(tǒng)線化理論逐漸轉(zhuǎn)向更為準(zhǔn)確的CFD技術(shù),其中陣風(fēng)效應(yīng)的嵌入主要分為遠(yuǎn)場邊界條件法和給定流場陣風(fēng)速度法兩類。

    將陣風(fēng)加入CFD程序最直接的方法是修改計(jì)算域遠(yuǎn)場邊界條件,稱為遠(yuǎn)場邊界條件(far-field boundary condition, FBC)法[2]。該方法為防止因計(jì)算域網(wǎng)格分辨率不足引起數(shù)值耗散導(dǎo)致陣風(fēng)響應(yīng)計(jì)算精度下降,要求計(jì)算域網(wǎng)格從物體表面到遠(yuǎn)場邊界均需要非常精細(xì),這將導(dǎo)致高昂的計(jì)算成本。為此,Tang和Bader[3]提出根據(jù)陣風(fēng)在計(jì)算域中所處位置的變化對計(jì)算網(wǎng)格進(jìn)行局部細(xì)化,該方法采用高階多項(xiàng)式插值來重新分布網(wǎng)格點(diǎn)以獲得陣風(fēng)區(qū)域更新后的精細(xì)網(wǎng)格[4]。

    給定速度法根據(jù)每一時(shí)刻陣風(fēng)場在計(jì)算域中的空間位置來計(jì)算陣風(fēng)速度,因此不存在遠(yuǎn)場邊界條件法中因遠(yuǎn)離物面區(qū)域網(wǎng)格分辨率不足而引起的數(shù)值耗散問題,被廣泛應(yīng)用于陣風(fēng)響應(yīng)數(shù)值計(jì)算。給定速度法包括網(wǎng)格速度法(field velocity method, FVM)和速度分裂法(split velocity method, SVM)。其中,F(xiàn)VM由Parameswaran和Bader提出[5],該方法借用動(dòng)網(wǎng)格系統(tǒng)流動(dòng)控制Euler/N-S方程,將陣風(fēng)速度等效為流場當(dāng)?shù)鼐W(wǎng)格速度來模擬陣風(fēng)運(yùn)動(dòng),但實(shí)際并不會(huì)引起網(wǎng)格運(yùn)動(dòng)。他們采用FVM模擬翼型穿越階躍陣風(fēng)時(shí)的氣動(dòng)響應(yīng),計(jì)算結(jié)果與理論分析結(jié)果十分接近。然而,F(xiàn)VM并非通過嚴(yán)格理論推導(dǎo)得來[6],也沒有考慮翼型對陣風(fēng)的影響。顧寧等[7]指出當(dāng)FVM應(yīng)用于N-S方程時(shí)需要加入黏性通量修正。相比而言,由Wales等[8]提出的SVM則是以非定常Euler/N-S方程為基礎(chǔ),將陣風(fēng)條件下的流場速度分解為陣風(fēng)速度與無陣風(fēng)作用下的流場背景速度兩部分,嚴(yán)格理論推導(dǎo)出SVM陣風(fēng)模擬控制方程。與FVM相比,除陣風(fēng)速度外,SVM還包括了陣風(fēng)對黏性通量的影響以及與陣風(fēng)時(shí)空導(dǎo)數(shù)相關(guān)的源項(xiàng)。此方法已被證明在較短陣風(fēng)尺度、跨聲速等條件下的陣風(fēng)響應(yīng)模擬中能夠獲得更為準(zhǔn)確的計(jì)算結(jié)果[8]。

    目前,國內(nèi)還尚未開展基于SVM的飛行器陣風(fēng)響應(yīng)模擬和減緩研究。本文針對One-Minus-Cosine離散陣風(fēng),首先基于動(dòng)態(tài)網(wǎng)格系統(tǒng)非定常N-S方程,理論推導(dǎo)SVM陣風(fēng)響應(yīng)模擬的控制方程。然后耦合翼型結(jié)構(gòu)運(yùn)動(dòng)方程,發(fā)展基于SVM的彈性翼型陣風(fēng)響應(yīng)CFD/CSD時(shí)域耦合算法。進(jìn)而采用控制翼型旋轉(zhuǎn),發(fā)展彈性翼型陣風(fēng)載荷減緩的數(shù)值模擬方法。最后采用上述方法,開展剛性及彈性翼型陣風(fēng)響應(yīng)和減緩分析。

    1 陣風(fēng)模型

    本文針對One-Minus-Cosine型離散陣風(fēng)開展翼型陣風(fēng)響應(yīng)分析。飛行器以速度V水平飛行,突然受到速度為wg的垂向陣風(fēng)作用。wg為陣風(fēng)速度Vg在直角坐標(biāo)系下的垂向分量,定義如下[9]:

    其中:Lg為陣風(fēng)尺度;wg=Vtan(Δα)為陣風(fēng)速度幅值;Δα為陣風(fēng)速度幅值等效迎角。

    2 陣風(fēng)響應(yīng)CFD模擬方法

    2.1 網(wǎng)格速度法(FVM)

    本節(jié)簡要給出FVM[10-11]的控制方程,便于與SVM方法進(jìn)行比較。動(dòng)態(tài)網(wǎng)格系統(tǒng)下非定常N-S方程為:

    其中:? 、為背景流場速度的速度分量;xt、zt為網(wǎng)格速度Vt的分量;ρ為流體密度;E為單位質(zhì)量流體總能;p為流體壓強(qiáng);γ為比熱比。

    根據(jù)相對運(yùn)動(dòng)原理,如果計(jì)算域網(wǎng)格的運(yùn)動(dòng)速度為Vt,則等價(jià)于在計(jì)算域靜止情況下受到速度為-Vt的氣流作用[12]。因此,基于動(dòng)態(tài)網(wǎng)格系統(tǒng)非定常N-S方程(2),F(xiàn)VM將陣風(fēng)速度表示為相反方向的網(wǎng)格速度來模擬陣風(fēng)作用[13],即:

    需要注意的是,F(xiàn)VM只是利用了網(wǎng)格速度的概念,在實(shí)施過程中并不引起實(shí)際的網(wǎng)格運(yùn)動(dòng)。盡管FVM是目前陣風(fēng)響應(yīng)CFD分析的主要方法,但該方法并未得到嚴(yán)格的理論闡明。

    2.2 速度分裂法(SVM)

    SVM將陣風(fēng)條件下的流場速度分解為陣風(fēng)速度與流場背景速度兩部分,即

    其中,u、w為流場總速度V的速度分量。為了開展考慮氣動(dòng)彈性的陣風(fēng)響應(yīng)分析,接下來將現(xiàn)有的基于固定網(wǎng)格系統(tǒng)的SVM方法[14]拓展至動(dòng)態(tài)網(wǎng)格系統(tǒng)。包含陣風(fēng)作用的動(dòng)態(tài)網(wǎng)格系統(tǒng)下非定常N-S方程可寫為:

    聯(lián)合式(5)~式(7),并利用背景流場的連續(xù)方程和動(dòng)量方程,可得SVM的控制方程為:

    其中黏性應(yīng)力張量的表達(dá)式為:

    其中,μ為黏性系數(shù)。式(8)中附加源項(xiàng)為:

    若采用Euler方程,黏性通量為0,則附加源項(xiàng)為:

    Sm、Se分別定義如下:

    不難看出,相比于SVM,傳統(tǒng)FVM的動(dòng)量方程和能量方程缺少了上述附加源項(xiàng)。該源項(xiàng)與陣風(fēng)場的時(shí)間和空間導(dǎo)數(shù)相關(guān),反映了背景流場對陣風(fēng)場的影響。對于小尺度陣風(fēng)或存在激波情形,附加源項(xiàng)的作用將會(huì)凸顯。

    本文離散陣風(fēng)形式由式(1)給定,因此源項(xiàng)中陣風(fēng)相關(guān)的導(dǎo)數(shù)可以解析求解。采用有限體積法、雙時(shí)間推進(jìn)法[15]求解包含陣風(fēng)作用的動(dòng)態(tài)網(wǎng)格系統(tǒng)下的非定常N-S方程(8),采用并行算法提高計(jì)算效率;基于RANS模擬湍流,湍流模型選用SA一方程模型[16]。遠(yuǎn)場采用無反射邊界條件,物面采用黏性無滑移邊界條件,即本文后續(xù)彈性翼型陣風(fēng)響應(yīng)及載荷減緩計(jì)算均采用無限插值法(TFI)[17]生成動(dòng)態(tài)網(wǎng)格。

    3 彈性翼型陣風(fēng)響應(yīng)CFD/CSD時(shí)域模擬方法

    3.1 三自由度結(jié)構(gòu)運(yùn)動(dòng)方程

    圖1給出了一個(gè)具有沉浮、俯仰和控制面偏轉(zhuǎn)的三自由度典型二元翼段力學(xué)模型。圖中三個(gè)自由度為:剛心E的沉浮位移h,向下為正;翼段繞剛心的俯仰角位移α,抬頭為正;控制面偏轉(zhuǎn)角 β ,向下偏轉(zhuǎn)為正。彈性軸E距離翼弦中點(diǎn)的距離為ahb,ah為基于半弦長b的無量綱距離,彈性軸在翼弦中點(diǎn)之后時(shí)ah>0。Cβb表示控制面鉸鏈軸距翼型弦線中點(diǎn)的距離,其氣動(dòng)彈性運(yùn)動(dòng)方程為[18-21]:

    圖1 三自由度典型翼段力學(xué)模型[21]Fig.1 Sketch of an airfoil with three degrees of freedom[21]

    其中:m為翼段質(zhì)量;Iα、Iβ分別為翼段繞剛心和控制面繞鉸鏈軸的轉(zhuǎn)動(dòng)慣量;Sα、Sβ分別為翼段和控制面的質(zhì)量靜矩;Th、Tα和Tβ分別為翼段的沉浮、俯仰和控制面轉(zhuǎn)動(dòng)阻尼系數(shù);Kh、Kα和Kβ分別為翼段的沉浮、俯仰和控制面轉(zhuǎn)動(dòng)剛度;L、Mα和Hβ分別為升力、俯仰力矩和鉸鏈力矩;CL、Cm和CH分別為升力系數(shù)、俯仰力矩系數(shù)和鉸鏈力矩系數(shù); ρ∞為自由來流密度;c為翼型弦長。

    式(16)的無量綱形式為:

    為了便于在時(shí)域中求解結(jié)構(gòu)運(yùn)動(dòng)方程(18),將其轉(zhuǎn)換為狀態(tài)空間形式:

    其中,I為單位矩陣。

    本文采用隱式歐拉法求解方程(20):

    得到n+ 1時(shí)刻的狀態(tài)變量:

    3.2 CFD/CSD時(shí)域耦合計(jì)算流程

    給定計(jì)算狀態(tài),首先不考慮陣風(fēng)作用,CFD計(jì)算獲得該狀態(tài)下的定常流場,然后基于定常流場開展陣風(fēng)響應(yīng)時(shí)域計(jì)算。基于定常流場,引入陣風(fēng)場,在時(shí)域內(nèi)耦合求解氣動(dòng)方程和結(jié)構(gòu)運(yùn)動(dòng)方程,從而求解出陣風(fēng)作用下翼型氣動(dòng)響應(yīng)和結(jié)構(gòu)位移響應(yīng)。

    4 基于翼型俯仰運(yùn)動(dòng)的陣風(fēng)減緩

    類似旋轉(zhuǎn)翼尖,本文嘗試通過控制翼型俯仰運(yùn)動(dòng)來進(jìn)行陣風(fēng)減緩。將翼型的俯仰控制角度αc與其位移響應(yīng)、速度響應(yīng)和加速度響應(yīng)聯(lián)系起來,

    其中,G1、G2、G3、G4、G5、G6為控制系統(tǒng)增益。

    則 式 (18) 右 側(cè) 的 σa應(yīng) 改 寫 為 σ =σa+σc, 其 中為控制力向量。

    5 算例分析

    5.1 剛性翼型陣風(fēng)響應(yīng)計(jì)算

    分別采用FVM和SVM方法,開展NACA0012剛性翼型的One-Minus-Cosine陣風(fēng)響應(yīng)計(jì)算分析。計(jì)算條件為:Ma= 0.3,迎角 0°,Re= 7.0 × 106;力矩參考點(diǎn)位于0.4c;定義無量綱時(shí)間S=2V∞t/c。采用C型結(jié)構(gòu)網(wǎng)格,陣風(fēng)從遠(yuǎn)場邊界外進(jìn)入計(jì)算域。根據(jù)文獻(xiàn)[8]中不同網(wǎng)格數(shù)對比計(jì)算結(jié)果,翼型周向和法向網(wǎng)格數(shù)取為257 × 129。翼型弦長c= 7.5m,遠(yuǎn)場邊界位于15倍弦長處,第一層網(wǎng)格高度為1×10-3c,如圖2所示。本節(jié)除黏性對陣風(fēng)響應(yīng)影響分析采用N-S方程計(jì)算之外,其余算例均采用Euler方程。

    圖2 NACA0012翼型計(jì)算網(wǎng)格Fig.2 Computational mesh for NACA0012

    5.1.1 時(shí)間步長的確定

    陣風(fēng)尺度Lg= 1c,等效迎角 Δα= 2.0°。分別取時(shí)間步長 ΔS為 0.1、0.01和 0.001,計(jì)算結(jié)果如圖3所示。對于FVM和SVM, ΔS= 0.1時(shí)的升力系數(shù)最大值和力矩系數(shù)最大值均大于 ΔS= 0.01或 ΔS= 0.001時(shí)的最大值,而 ΔS= 0.01與 ΔS= 0.001的結(jié)果相吻合,因此本文時(shí)間步長取 ΔS= 0.01。

    圖3 時(shí)間步長對FVM和SVM陣風(fēng)響應(yīng)的影響Fig.3 Effects of time step on gust response obtained by FVM and SVM

    圖4比較了陣風(fēng)尺度為1c時(shí)本文計(jì)算得到的升力系數(shù)陣風(fēng)響應(yīng)與文獻(xiàn)[8]的結(jié)果。對于小尺度陣風(fēng),陣風(fēng)梯度增大、作用域更為集中,在陣風(fēng)離開翼型后緣的時(shí)間段內(nèi),升力系數(shù)產(chǎn)生小幅波動(dòng);除波動(dòng)幅值存在一定差異外,本文FVM和SVM兩種方法計(jì)算得到的升力系數(shù)峰值及其出現(xiàn)時(shí)刻均與文獻(xiàn)結(jié)果吻合良好。

    圖4 FVM和SVM氣動(dòng)力陣風(fēng)響應(yīng)與文獻(xiàn)結(jié)果對比Fig.4 Comparisons of aerodynamic responses between present and ref.[8]data

    進(jìn)一步研究表明,陣風(fēng)離開翼型后緣時(shí)升力系數(shù)小幅波動(dòng)的幅值與翼型后緣厚度的處理方式直接相關(guān)。圖5比較了尖后緣和鈍后緣NACA0012翼型計(jì)算得到的升力系數(shù)陣風(fēng)響應(yīng),鈍后緣翼型后緣厚度為1.26‰c。相比于尖后緣,陣風(fēng)離開鈍后緣翼型時(shí)升力系數(shù)響應(yīng)的波動(dòng)幅值更大,峰值略微降低。

    5.1.2 陣風(fēng)尺度的影響

    分別采用FVM、SVM計(jì)算陣風(fēng)尺度 1c、5c和10c,等效迎角 Δα= 2.0°的陣風(fēng)響應(yīng),計(jì)算結(jié)果如圖6所示。在陣風(fēng)尺度5c和10c時(shí),F(xiàn)VM和SVM的計(jì)算結(jié)果幾乎一致。當(dāng)陣風(fēng)尺度減小到1c時(shí),陣風(fēng)場的時(shí)空導(dǎo)數(shù)增大,SVM方法中源項(xiàng)的作用凸顯出來,導(dǎo)致FVM和SVM的計(jì)算結(jié)果出現(xiàn)差異,SVM計(jì)算得到的氣動(dòng)參數(shù)陣風(fēng)響應(yīng)峰值略高于FVM。SVM方法的推導(dǎo)過程表明FVM只是SVM忽略源項(xiàng)的一種近似方法,因此理論上SVM的計(jì)算結(jié)果更為準(zhǔn)確。隨著陣風(fēng)尺度的增加,陣風(fēng)場的時(shí)空導(dǎo)數(shù)減小,陣風(fēng)離開翼型后緣時(shí)氣動(dòng)系數(shù)的小幅波動(dòng)現(xiàn)象消失。

    圖6 陣風(fēng)尺度對FVM和SVM的陣風(fēng)響應(yīng)影響Fig.6 Gust length effect on gust responses obtained by FVM and SVM

    5.1.3 黏性對陣風(fēng)響應(yīng)的影響

    針對陣風(fēng)尺度1c和10c、等效迎角 Δα= 2.0°,采用SVM方法開展Euler、N-S方程的陣風(fēng)響應(yīng)對比計(jì)算,其中N-S方程計(jì)算網(wǎng)格的第一層網(wǎng)格高度為1×10-5c,其他網(wǎng)格參數(shù)與Euler方程計(jì)算網(wǎng)格相同。計(jì)算結(jié)果對比如圖7所示,黏性效應(yīng)導(dǎo)致陣風(fēng)響應(yīng)的升力系數(shù)峰值降低,陣風(fēng)尺度越大,黏性對陣風(fēng)場的影響越大,升力減小量越明顯;黏性作用下,1c陣風(fēng)離開翼型后緣時(shí)的升力系數(shù)波動(dòng)明顯減弱。

    圖7 黏性對陣風(fēng)響應(yīng)的影響Fig.7 Viscosity effects on gust response

    5.2 彈性翼型陣風(fēng)響應(yīng)計(jì)算

    采用非定常Euler方程計(jì)算彈性NACA0012翼型的陣風(fēng)響應(yīng)。翼型結(jié)構(gòu)模型參數(shù)[22]由表1給出。

    表1 NACA0012翼型結(jié)構(gòu)參數(shù)Table 1 NACA0012 parameters

    5.2.1 陣風(fēng)響應(yīng)CFD/CSD耦合算法驗(yàn)證

    選取文獻(xiàn)[22]的計(jì)算參數(shù),Ma= 0.3,V*= 0.2,翼型弦長c= 1.0 m。翼型結(jié)構(gòu)位移的初始狀態(tài) ξ (0) = 0,α(0)= 0°,陣風(fēng)尺度Lg= 6.25 m,等效迎角 Δα= 0.57°??紤]到文獻(xiàn)計(jì)算方法為FVM,本算例也采用FVM。計(jì)算結(jié)果對比如圖8所示。本文計(jì)算結(jié)果與文獻(xiàn)[22]數(shù)據(jù)吻合,驗(yàn)證了本文彈性翼型陣風(fēng)響應(yīng)CFD/CSD耦合算法的可靠性。

    圖8 彈性翼型陣風(fēng)響應(yīng)與文獻(xiàn)結(jié)果對比Fig.8 Comparison of gust responses of an elastic airfoil between present and reference results

    5.2.2 彈性對翼型陣風(fēng)響應(yīng)的影響

    采用SVM計(jì)算陣風(fēng)尺度1c和10c、等效迎角 Δα=2.0°的陣風(fēng)響應(yīng),V*= 0.2,翼型結(jié)構(gòu)初始狀態(tài) ξ (0) =0,α(0)= 0°。圖9為剛性翼型和彈性翼型的氣動(dòng)系數(shù)陣風(fēng)響應(yīng)對比。圖10為彈性翼型的結(jié)構(gòu)響應(yīng)結(jié)果。在One-Minus-Cosine陣風(fēng)穿越彈性翼型期間,翼型的氣動(dòng)參數(shù)和俯仰角均先增大然后減小至極小值,沉浮位移則一直增大至峰值;10c陣風(fēng)穿越彈性翼型時(shí),升力和俯仰角近乎同時(shí)達(dá)到峰值,導(dǎo)致升力峰值較剛性翼型增大;1c陣風(fēng)穿越彈性翼型時(shí),俯仰角峰值遠(yuǎn)滯后于升力峰值,彈性對升力峰值幾乎沒有影響。

    圖9 剛性和彈性翼型氣動(dòng)力陣風(fēng)響應(yīng)對比Fig.9 Comparison of aerodynamic gust response between rigid and elastic airfoils

    圖10 彈性翼型結(jié)構(gòu)響應(yīng)Fig.10 Structural response of an elastic airfoil

    5.3 彈性翼型陣風(fēng)載荷減緩

    以NACA64A010翼型為研究對象,翼型弦長c=1.0 m,其他網(wǎng)格參數(shù)與第5.1節(jié)Euler方程計(jì)算網(wǎng)格相同,如圖11所示。翼型結(jié)構(gòu)參數(shù)[23]如表2所示,取Ma= 0.3、V*= 0.337 4。采用SVM進(jìn)行翼型陣風(fēng)載荷減緩響應(yīng)計(jì)算。陣風(fēng)尺度Lg= 10c,等效迎角 Δα=2.0°。翼型結(jié)構(gòu)位移 的 初 始狀態(tài) ξ (0) = 0、α(0) = 0°、β(0)= 0°。本文以減緩翼型在通過陣風(fēng)時(shí)產(chǎn)生的最大升力系數(shù)為目標(biāo)。

    表2 NACA64A010翼型結(jié)構(gòu)參數(shù)Table 2 NACA64A010 parameters

    圖11 NACA64A010翼型計(jì)算網(wǎng)格Fig.11 Computational mesh for NACA64A010

    彈性翼型陣風(fēng)響應(yīng)分析表明,在陣風(fēng)穿越翼型期間,升力系數(shù)先增大后減小至極小值,而翼型沉浮位移則一直增大至極值,因此本文選取沉浮速度值作為控制輸入量。圖12、圖13為增益G3分別取5.0和10.0時(shí)閉環(huán)控制系統(tǒng)與開環(huán)系統(tǒng)陣風(fēng)響應(yīng)的結(jié)果對比。采用G3閉環(huán)控制情況下,沉浮運(yùn)動(dòng)峰值顯著降低并迅速收斂;在升力系數(shù)峰值之前,由于翼型沉浮速度為負(fù)值,G3等效于給彈性翼型增加低頭俯仰控制力矩,導(dǎo)致俯仰角、升力和力矩系數(shù)峰值明顯降低,G3= 10時(shí)升力系數(shù)峰值降低49.758%;當(dāng)陣風(fēng)離開翼型后,翼型沉浮速度為正值,俯仰控制力矩使翼型俯仰角逐漸增大至正極大值點(diǎn),并高于無控制情況,導(dǎo)致此時(shí)氣動(dòng)系數(shù)響應(yīng)極值也均略高于無控制情況,并接近或高于陣風(fēng)穿越翼型時(shí)的峰值。

    圖12 開環(huán)系統(tǒng)和閉環(huán)控制系統(tǒng)翼型氣動(dòng)力陣風(fēng)響應(yīng)對比Fig.12 Comparison of aerodynamic gust responses between open- and closed-loop control systems

    圖13 開環(huán)系統(tǒng)和閉環(huán)控制系統(tǒng)翼型結(jié)構(gòu)響應(yīng)對比Fig.13 Compsrison of structural response between open- and closed-loop control systems

    接下來采用G2= -5.0、G3= 10.0同時(shí)控制翼型的沉浮和俯仰運(yùn)動(dòng)。圖14、圖15分別為G2、G3組合的閉環(huán)控制陣風(fēng)響應(yīng)結(jié)果和僅G3的響應(yīng)結(jié)果對比。增加俯仰控制后,俯仰角位移幅值顯著減小,導(dǎo)致陣風(fēng)作用期間的氣動(dòng)參數(shù)峰值高于無俯仰控制情況,但陣風(fēng)離開翼型后的氣動(dòng)參數(shù)極值則顯著降低。

    圖14 G2 = -5.0、G3 = 10.0和僅有G3 = 10.0氣動(dòng)力響應(yīng)對比Fig.14 Comparison of aerodynamic response between G2 = -5.0, G3 = 10.0 and only G3 = 10.0

    圖15 G2 = -5.0、G3 = 10.0和僅有G3 = 10.0結(jié)構(gòu)響應(yīng)對比Fig.15 Comparison of structural response between G2 = -5.0, G3 = 10.0 and only G3 = 10.0

    6 結(jié) 論

    本文基于動(dòng)態(tài)網(wǎng)格系統(tǒng)下非定常N-S方程,理論推導(dǎo)出SVM陣風(fēng)模擬的控制方程,發(fā)展出基于SVM的彈性翼型陣風(fēng)響應(yīng)CFD/CSD時(shí)域耦合算法和基于翼型俯仰控制的陣風(fēng)減緩模擬方法。算例分析了NACA0012剛性、彈性翼型的One-Minus-Cosine陣風(fēng)響應(yīng),計(jì)算結(jié)果與文獻(xiàn)數(shù)據(jù)吻合;分析了陣風(fēng)尺度、流體黏性和結(jié)構(gòu)彈性對NACA0012翼型陣風(fēng)響應(yīng)的影響;實(shí)現(xiàn)了基于俯仰控制的NACA64A010彈性翼型陣風(fēng)減緩。主要結(jié)論如下:

    1) 理論推導(dǎo)表明,傳統(tǒng)FVM是忽略SVM中陣風(fēng)相關(guān)附加源項(xiàng)的一種近似方法,對于大尺度陣風(fēng),兩種方法的計(jì)算結(jié)果一致,源項(xiàng)的作用隨著陣風(fēng)尺度的減小而凸顯,理論上SVM更為準(zhǔn)確;

    2) 黏性效應(yīng)導(dǎo)致NACA0012翼型升力系數(shù)陣風(fēng)響應(yīng)峰值降低,陣風(fēng)尺度越大,升力減小越明顯,黏性同時(shí)減弱了1c陣風(fēng)離開翼型后緣時(shí)的升力系數(shù)波動(dòng);

    3) 結(jié)構(gòu)彈性引起NACA0012翼型俯仰角陣風(fēng)響應(yīng)峰值增加,陣風(fēng)尺度越大,升力系數(shù)增加越顯著;

    4) 俯仰控制能夠有效減緩NACA64A010彈性翼型的陣風(fēng)載荷,沉浮速度控制輸入量G3= 10.0時(shí)升力系數(shù)陣風(fēng)響應(yīng)峰值減緩49.758%,結(jié)合G2、G3則能夠同時(shí)減緩沉浮和俯仰響應(yīng)。

    猜你喜歡
    陣風(fēng)升力黏性
    高速列車車頂–升力翼組合體氣動(dòng)特性
    陣風(fēng)戰(zhàn)斗機(jī)
    軍事文摘(2021年19期)2021-10-10 13:29:36
    法國陣風(fēng)戰(zhàn)斗機(jī)
    軍事文摘(2021年17期)2021-09-24 01:11:46
    無人機(jī)升力測試裝置設(shè)計(jì)及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進(jìn)入段快速軌跡優(yōu)化
    富硒產(chǎn)業(yè)需要強(qiáng)化“黏性”——安康能否玩轉(zhuǎn)“硒+”
    如何運(yùn)用播音主持技巧增強(qiáng)受眾黏性
    傳媒評論(2019年4期)2019-07-13 05:49:28
    陣風(fēng)勁吹
    航空世界(2018年12期)2018-07-16 08:34:50
    玩油灰黏性物成網(wǎng)紅
    基層農(nóng)行提高客戶黏性淺析
    欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 日韩成人av中文字幕在线观看| 午夜福利在线观看免费完整高清在| 亚洲成av片中文字幕在线观看 | 久久精品久久久久久噜噜老黄| 国产精品麻豆人妻色哟哟久久| 久久综合国产亚洲精品| 国产精品免费大片| av免费观看日本| 国产免费现黄频在线看| 国产精品免费视频内射| 久久久精品94久久精品| 色婷婷av一区二区三区视频| 视频在线观看一区二区三区| 欧美日韩综合久久久久久| 日日啪夜夜爽| 下体分泌物呈黄色| 中国国产av一级| 久久综合国产亚洲精品| 日韩av免费高清视频| 日本色播在线视频| 边亲边吃奶的免费视频| 丝袜喷水一区| 99久久人妻综合| 99久久精品国产国产毛片| 久久韩国三级中文字幕| 欧美少妇被猛烈插入视频| 欧美+日韩+精品| 欧美在线黄色| 黑人巨大精品欧美一区二区蜜桃| 国产黄色视频一区二区在线观看| 91精品伊人久久大香线蕉| 久久这里只有精品19| 日本av免费视频播放| 中文字幕人妻丝袜制服| 亚洲精品久久成人aⅴ小说| 99热全是精品| 麻豆乱淫一区二区| 有码 亚洲区| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美在线一区| 国产97色在线日韩免费| 男女免费视频国产| 欧美97在线视频| 亚洲精品成人av观看孕妇| 国精品久久久久久国模美| 久久久久久久久免费视频了| 精品少妇久久久久久888优播| 中国国产av一级| 制服人妻中文乱码| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美精品自产自拍| 亚洲美女搞黄在线观看| 天堂俺去俺来也www色官网| 精品少妇久久久久久888优播| 久久精品国产自在天天线| 亚洲综合色惰| 精品国产乱码久久久久久男人| 欧美人与善性xxx| 精品视频人人做人人爽| 另类精品久久| 国产 精品1| 在线 av 中文字幕| 最近中文字幕高清免费大全6| 国产激情久久老熟女| 精品少妇内射三级| 狂野欧美激情性bbbbbb| 在现免费观看毛片| 久久精品久久精品一区二区三区| 2021少妇久久久久久久久久久| 日本av手机在线免费观看| 午夜激情久久久久久久| 777米奇影视久久| 9色porny在线观看| 99九九在线精品视频| 曰老女人黄片| 亚洲精品第二区| 夜夜骑夜夜射夜夜干| 亚洲欧美一区二区三区久久| 久久韩国三级中文字幕| 91精品国产国语对白视频| 日本av手机在线免费观看| 国产精品三级大全| 亚洲精品久久午夜乱码| 丝袜美足系列| 老司机影院成人| 亚洲欧美精品综合一区二区三区 | 久久久久久人人人人人| 一级黄片播放器| 亚洲欧美精品综合一区二区三区 | 18禁动态无遮挡网站| 我要看黄色一级片免费的| 国产老妇伦熟女老妇高清| 欧美成人午夜精品| 欧美日韩精品网址| 综合色丁香网| 亚洲精品国产色婷婷电影| 久久青草综合色| 中文字幕最新亚洲高清| 嫩草影院入口| 国产综合精华液| 高清av免费在线| 亚洲av国产av综合av卡| 永久免费av网站大全| 大话2 男鬼变身卡| 香蕉国产在线看| 五月伊人婷婷丁香| 亚洲美女视频黄频| 欧美中文综合在线视频| 成人国产av品久久久| 成人漫画全彩无遮挡| 下体分泌物呈黄色| 蜜桃在线观看..| 只有这里有精品99| 校园人妻丝袜中文字幕| 亚洲欧美色中文字幕在线| 成人毛片60女人毛片免费| 亚洲综合色网址| 午夜激情av网站| 亚洲精品一二三| 国产熟女午夜一区二区三区| 2021少妇久久久久久久久久久| 亚洲第一av免费看| 99国产综合亚洲精品| 高清欧美精品videossex| 日韩一卡2卡3卡4卡2021年| 亚洲av.av天堂| 大片免费播放器 马上看| 亚洲欧美精品综合一区二区三区 | 亚洲精品一二三| 老鸭窝网址在线观看| 日韩制服骚丝袜av| 黄色怎么调成土黄色| a 毛片基地| 天天影视国产精品| 欧美日韩一区二区视频在线观看视频在线| 日韩中字成人| 2018国产大陆天天弄谢| 最近2019中文字幕mv第一页| 人妻 亚洲 视频| 亚洲久久久国产精品| 纵有疾风起免费观看全集完整版| 成人黄色视频免费在线看| 丝袜人妻中文字幕| 永久免费av网站大全| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| 国产精品久久久久久av不卡| 麻豆精品久久久久久蜜桃| 少妇精品久久久久久久| 久久亚洲国产成人精品v| 午夜久久久在线观看| av又黄又爽大尺度在线免费看| av在线播放精品| 人妻系列 视频| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 久久久精品94久久精品| 99热全是精品| 亚洲久久久国产精品| 日本色播在线视频| 综合色丁香网| 秋霞在线观看毛片| 涩涩av久久男人的天堂| 热99久久久久精品小说推荐| 老汉色av国产亚洲站长工具| 黄色配什么色好看| 亚洲av欧美aⅴ国产| 久久久久久久久久久久大奶| 亚洲天堂av无毛| 日韩精品有码人妻一区| 日韩中文字幕视频在线看片| 卡戴珊不雅视频在线播放| 在线观看www视频免费| 久久婷婷青草| 日本-黄色视频高清免费观看| 日本wwww免费看| 久久97久久精品| 久久久久久免费高清国产稀缺| 另类亚洲欧美激情| 18在线观看网站| 91精品国产国语对白视频| 欧美av亚洲av综合av国产av | 成人国语在线视频| 人妻系列 视频| 看免费成人av毛片| 多毛熟女@视频| 亚洲精品国产av蜜桃| 最近最新中文字幕大全免费视频 | 成人国产av品久久久| 男女午夜视频在线观看| 亚洲熟女精品中文字幕| 丝袜脚勾引网站| 黄色一级大片看看| av天堂久久9| 日日啪夜夜爽| 免费人妻精品一区二区三区视频| 在线 av 中文字幕| 国产精品.久久久| 大片电影免费在线观看免费| 18禁动态无遮挡网站| 午夜福利在线观看免费完整高清在| av天堂久久9| 黑人巨大精品欧美一区二区蜜桃| 成人手机av| av卡一久久| 久久久精品免费免费高清| 免费播放大片免费观看视频在线观看| 黄网站色视频无遮挡免费观看| 夫妻性生交免费视频一级片| 色视频在线一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美成人综合另类久久久| 黄色视频在线播放观看不卡| 国产成人精品在线电影| 久久久久久久久免费视频了| 久久久久久久亚洲中文字幕| 少妇人妻 视频| 在线免费观看不下载黄p国产| 黄色配什么色好看| 美女xxoo啪啪120秒动态图| 欧美av亚洲av综合av国产av | 色视频在线一区二区三区| 一级爰片在线观看| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美视频二区| 日本-黄色视频高清免费观看| 亚洲情色 制服丝袜| tube8黄色片| 欧美日韩视频精品一区| 亚洲精品第二区| 叶爱在线成人免费视频播放| 国产成人av激情在线播放| 免费少妇av软件| 哪个播放器可以免费观看大片| 不卡av一区二区三区| 国产一区二区三区综合在线观看| 男女国产视频网站| videos熟女内射| 欧美日韩视频精品一区| 一边摸一边做爽爽视频免费| 亚洲av欧美aⅴ国产| 好男人视频免费观看在线| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 中文字幕精品免费在线观看视频| 成人免费观看视频高清| 在线观看美女被高潮喷水网站| 国产精品麻豆人妻色哟哟久久| 观看av在线不卡| 午夜激情av网站| 人体艺术视频欧美日本| 久久 成人 亚洲| 久久久欧美国产精品| 日韩大片免费观看网站| 国产在线一区二区三区精| 国产精品偷伦视频观看了| www日本在线高清视频| av有码第一页| 精品少妇久久久久久888优播| 成年女人在线观看亚洲视频| 日韩大片免费观看网站| 国产乱人偷精品视频| 亚洲精品,欧美精品| 国产精品熟女久久久久浪| 日韩熟女老妇一区二区性免费视频| 亚洲精品国产色婷婷电影| 国产精品久久久av美女十八| 亚洲av男天堂| www.自偷自拍.com| 日韩中字成人| 香蕉精品网在线| 日韩av在线免费看完整版不卡| 一区二区三区乱码不卡18| 99国产精品免费福利视频| 国产精品欧美亚洲77777| 女人高潮潮喷娇喘18禁视频| 在线观看免费日韩欧美大片| www.熟女人妻精品国产| 久久韩国三级中文字幕| 久久精品夜色国产| 国产亚洲精品第一综合不卡| 中文字幕另类日韩欧美亚洲嫩草| 免费观看a级毛片全部| 日韩大片免费观看网站| 97在线视频观看| 色婷婷久久久亚洲欧美| 美女视频免费永久观看网站| 国产xxxxx性猛交| 国产成人一区二区在线| 国产av精品麻豆| 亚洲欧美中文字幕日韩二区| 99热全是精品| 日韩欧美一区视频在线观看| 国产av精品麻豆| 自拍欧美九色日韩亚洲蝌蚪91| 老司机影院成人| 日韩电影二区| 精品人妻偷拍中文字幕| 美女福利国产在线| 国产精品麻豆人妻色哟哟久久| 这个男人来自地球电影免费观看 | 99香蕉大伊视频| 国产成人av激情在线播放| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 日本av免费视频播放| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美| 日韩熟女老妇一区二区性免费视频| 成人手机av| 久久影院123| 久久精品国产综合久久久| 欧美日韩亚洲高清精品| 亚洲精品日韩在线中文字幕| 免费黄色在线免费观看| 国产女主播在线喷水免费视频网站| 国产一区二区 视频在线| 老熟女久久久| 免费观看a级毛片全部| 久久久精品区二区三区| 国产黄频视频在线观看| 国产无遮挡羞羞视频在线观看| 日韩一本色道免费dvd| 亚洲欧美清纯卡通| 亚洲av电影在线进入| 日日爽夜夜爽网站| 十八禁网站网址无遮挡| 久久综合国产亚洲精品| 尾随美女入室| 国产精品女同一区二区软件| 久久精品国产亚洲av涩爱| 成人影院久久| 18+在线观看网站| 高清av免费在线| 国产xxxxx性猛交| 色吧在线观看| 日本欧美视频一区| 欧美精品国产亚洲| 免费少妇av软件| 少妇的丰满在线观看| www.自偷自拍.com| av片东京热男人的天堂| 熟妇人妻不卡中文字幕| 成人国产麻豆网| 狠狠精品人妻久久久久久综合| 街头女战士在线观看网站| 一级黄片播放器| 欧美精品人与动牲交sv欧美| 高清欧美精品videossex| 国产av码专区亚洲av| 亚洲婷婷狠狠爱综合网| 久久国产精品大桥未久av| 亚洲婷婷狠狠爱综合网| 色网站视频免费| 男男h啪啪无遮挡| 免费大片黄手机在线观看| 观看美女的网站| 国产毛片在线视频| 午夜免费鲁丝| 美女高潮到喷水免费观看| 国产精品一区二区在线不卡| 免费在线观看完整版高清| 欧美精品av麻豆av| 热99久久久久精品小说推荐| 久久午夜综合久久蜜桃| 亚洲,欧美,日韩| 国产片内射在线| 久久久精品免费免费高清| 国产精品欧美亚洲77777| 一区福利在线观看| 激情五月婷婷亚洲| 欧美激情 高清一区二区三区| 国产精品三级大全| 亚洲第一青青草原| 国产高清国产精品国产三级| 亚洲精品久久成人aⅴ小说| 国产亚洲午夜精品一区二区久久| 另类亚洲欧美激情| 精品一区在线观看国产| 国产精品不卡视频一区二区| 久久影院123| 欧美老熟妇乱子伦牲交| 日本av手机在线免费观看| 99国产精品免费福利视频| 精品少妇久久久久久888优播| 免费在线观看完整版高清| 国产亚洲最大av| 亚洲精华国产精华液的使用体验| 久久久久人妻精品一区果冻| 亚洲,欧美,日韩| 午夜免费鲁丝| 黄色配什么色好看| 丰满少妇做爰视频| 中文字幕人妻丝袜制服| 男人爽女人下面视频在线观看| 国产成人精品无人区| 欧美国产精品va在线观看不卡| 黄频高清免费视频| 亚洲欧美精品综合一区二区三区 | 国产免费一区二区三区四区乱码| 人人妻人人澡人人爽人人夜夜| 一二三四中文在线观看免费高清| videosex国产| 国产在视频线精品| 亚洲欧美一区二区三区黑人 | 1024视频免费在线观看| 欧美成人午夜精品| 成人黄色视频免费在线看| 午夜精品国产一区二区电影| 又大又黄又爽视频免费| av在线app专区| 久久久久视频综合| 久久这里有精品视频免费| 亚洲欧洲精品一区二区精品久久久 | 午夜老司机福利剧场| 99国产综合亚洲精品| 久久午夜福利片| 精品国产露脸久久av麻豆| 一级毛片我不卡| 久久久久久久久久人人人人人人| 免费在线观看视频国产中文字幕亚洲 | 久久久久久久精品精品| videosex国产| 欧美人与性动交α欧美软件| 麻豆av在线久日| 亚洲精品成人av观看孕妇| 国产精品一国产av| 天天躁夜夜躁狠狠躁躁| 国产精品亚洲av一区麻豆 | 高清视频免费观看一区二区| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花| 夜夜骑夜夜射夜夜干| 高清黄色对白视频在线免费看| 久久久久国产精品人妻一区二区| 香蕉国产在线看| 日韩制服丝袜自拍偷拍| 亚洲成国产人片在线观看| 男人添女人高潮全过程视频| 免费av中文字幕在线| 18禁国产床啪视频网站| 国产av国产精品国产| 最近中文字幕2019免费版| 婷婷色综合www| 在线看a的网站| 日韩 亚洲 欧美在线| 欧美精品人与动牲交sv欧美| 丰满饥渴人妻一区二区三| 久久精品亚洲av国产电影网| 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 青春草视频在线免费观看| 男女国产视频网站| 日本黄色日本黄色录像| 日韩欧美精品免费久久| 狠狠婷婷综合久久久久久88av| 日韩 亚洲 欧美在线| 亚洲中文av在线| 久久影院123| 少妇的逼水好多| av一本久久久久| 美女中出高潮动态图| 十分钟在线观看高清视频www| 久久久久久人妻| 欧美精品一区二区免费开放| a 毛片基地| 亚洲一码二码三码区别大吗| 好男人视频免费观看在线| 看免费成人av毛片| 亚洲av电影在线观看一区二区三区| 国产日韩欧美在线精品| 国产97色在线日韩免费| 999久久久国产精品视频| 国产 一区精品| 青春草国产在线视频| 亚洲国产精品一区二区三区在线| 青青草视频在线视频观看| 考比视频在线观看| 国产成人91sexporn| 韩国精品一区二区三区| av在线观看视频网站免费| 女性被躁到高潮视频| 久久人妻熟女aⅴ| 一区二区三区乱码不卡18| 超碰97精品在线观看| 又大又黄又爽视频免费| 国产又爽黄色视频| 色视频在线一区二区三区| 免费高清在线观看日韩| 黄色视频在线播放观看不卡| 两性夫妻黄色片| 久久精品国产亚洲av天美| 欧美精品国产亚洲| 捣出白浆h1v1| 国产精品久久久久成人av| 国产女主播在线喷水免费视频网站| 久久久久久免费高清国产稀缺| 亚洲av成人精品一二三区| 亚洲美女搞黄在线观看| 亚洲av成人精品一二三区| tube8黄色片| www.熟女人妻精品国产| 亚洲,欧美精品.| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 日韩免费高清中文字幕av| 日韩一本色道免费dvd| 18禁观看日本| 自拍欧美九色日韩亚洲蝌蚪91| 欧美人与善性xxx| 一级毛片我不卡| 亚洲精品久久久久久婷婷小说| 久久精品国产a三级三级三级| 亚洲内射少妇av| 免费大片黄手机在线观看| 国语对白做爰xxxⅹ性视频网站| 亚洲av国产av综合av卡| 九九爱精品视频在线观看| 欧美中文综合在线视频| 黄色配什么色好看| 9热在线视频观看99| 一区二区三区精品91| 精品亚洲成a人片在线观看| 男女国产视频网站| 久久久久久久久久人人人人人人| 少妇人妻精品综合一区二区| 极品人妻少妇av视频| 亚洲图色成人| 欧美中文综合在线视频| 欧美激情极品国产一区二区三区| 老司机亚洲免费影院| 国产一区二区三区av在线| 十八禁高潮呻吟视频| 黑人猛操日本美女一级片| 人妻 亚洲 视频| 91精品国产国语对白视频| 边亲边吃奶的免费视频| 久久久久久人人人人人| 国产亚洲午夜精品一区二区久久| 国产一区二区激情短视频 | 男女边摸边吃奶| 国产精品久久久久久久久免| 国产一区亚洲一区在线观看| 久久精品国产亚洲av涩爱| 99久久精品国产国产毛片| 街头女战士在线观看网站| 高清在线视频一区二区三区| 中国国产av一级| 国产白丝娇喘喷水9色精品| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 一区二区日韩欧美中文字幕| 美女福利国产在线| 男女边摸边吃奶| 乱人伦中国视频| 久久影院123| 国产有黄有色有爽视频| 婷婷成人精品国产| 国产男女内射视频| 国产成人精品福利久久| 少妇 在线观看| 日本wwww免费看| a级片在线免费高清观看视频| 午夜免费鲁丝| av网站免费在线观看视频| 亚洲精品中文字幕在线视频| 国产一区亚洲一区在线观看| www.av在线官网国产| 久久精品夜色国产| 久久久精品国产亚洲av高清涩受| 1024香蕉在线观看| 大片免费播放器 马上看| 这个男人来自地球电影免费观看 | 亚洲国产日韩一区二区| videossex国产| 国产黄频视频在线观看| 亚洲欧美色中文字幕在线| 亚洲,一卡二卡三卡| 久久精品国产亚洲av高清一级| 如何舔出高潮| 男女国产视频网站| 99国产综合亚洲精品| 老汉色∧v一级毛片| 桃花免费在线播放| 国产精品二区激情视频| 一本大道久久a久久精品| 日韩中字成人| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 纵有疾风起免费观看全集完整版| 午夜精品国产一区二区电影| 国产人伦9x9x在线观看 | 午夜激情av网站| 少妇被粗大猛烈的视频| 久久婷婷青草| 黄色毛片三级朝国网站| 2022亚洲国产成人精品| 亚洲精品一区蜜桃| 校园人妻丝袜中文字幕| 在线看a的网站| av线在线观看网站| 国产一区二区 视频在线| 纵有疾风起免费观看全集完整版| 亚洲精品国产av蜜桃| 亚洲av电影在线观看一区二区三区| 熟女电影av网| 国产欧美日韩一区二区三区在线| 午夜激情久久久久久久| 久久久久国产网址| 秋霞在线观看毛片| 国产精品人妻久久久影院|