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

    運(yùn)載火箭與航天器界面載荷設(shè)計(jì)優(yōu)化方法

    2023-04-15 13:04:28曾耀祥
    宇航學(xué)報(bào) 2023年3期
    關(guān)鍵詞:基頻質(zhì)心航天器

    何 巍,曾耀祥,劉 暉

    (北京宇航系統(tǒng)工程研究所,北京 100076)

    0 引 言

    運(yùn)載火箭在跨聲速飛行時(shí),受到脈動(dòng)壓力等外激勵(lì)影響[1],火箭及搭載的航天器會(huì)產(chǎn)生明顯的動(dòng)態(tài)加速度和內(nèi)力響應(yīng)。動(dòng)態(tài)內(nèi)力是航天器主體結(jié)構(gòu)、火箭適配器、器箭連接裝置等航天器內(nèi)分支結(jié)構(gòu)設(shè)計(jì)的重要輸入條件,而準(zhǔn)確合理地給出動(dòng)態(tài)內(nèi)力設(shè)計(jì)值能有效降低非必要結(jié)構(gòu)重量、有效提升有效載荷質(zhì)量,對(duì)提升器箭結(jié)構(gòu)效率和整體設(shè)計(jì)性能有著十分重要的意義。

    國(guó)內(nèi)外火箭用戶手冊(cè)中以航天器質(zhì)心等效加速度的形式規(guī)定器箭界面載荷條件,根據(jù)航天器質(zhì)量、質(zhì)心高度以及質(zhì)心加速度位置則可計(jì)算得到界面力及界面力矩。目前大多數(shù)火箭用戶手冊(cè)中,該條件對(duì)于任意航天器均取相同值[2-4],根據(jù)動(dòng)力學(xué)理論,航天器質(zhì)量、質(zhì)心及基頻等質(zhì)量與頻率特性必然會(huì)對(duì)動(dòng)態(tài)內(nèi)力響應(yīng)產(chǎn)生影響,而當(dāng)前的載荷設(shè)計(jì)條件未考慮航天器特性的影響。近些年,部分火箭的用戶手冊(cè)開(kāi)始依據(jù)不同的航天器參數(shù)區(qū)間給出界面載荷條件,宇宙神5火箭和“質(zhì)子”火箭針對(duì)不同的航天器質(zhì)量范圍(按照是否大于3175 kg進(jìn)行劃分)給出了不同的橫向質(zhì)心等效過(guò)載;同樣,獵鷹9火箭按照航天器質(zhì)量是否大于1814 kg給出不同的橫向質(zhì)心等效過(guò)載,不同區(qū)間對(duì)應(yīng)的橫向質(zhì)心等效過(guò)載最大差異可達(dá)100%,但用戶手冊(cè)中未對(duì)劃分依據(jù)進(jìn)行說(shuō)明[5-7]。由于動(dòng)態(tài)內(nèi)力隨自身特性參數(shù)變化規(guī)律的復(fù)雜性,需要開(kāi)展深入研究。

    器箭組合體動(dòng)力學(xué)模型復(fù)雜,直接通過(guò)理論分析研究質(zhì)心等效加速度的變化趨勢(shì)十分困難,考慮到器箭組合體屬于長(zhǎng)細(xì)比很大的薄壁結(jié)構(gòu),符合梁模型動(dòng)力學(xué)基本假設(shè),因此可基于梁動(dòng)力學(xué)理論分析航天器分支特性參數(shù)對(duì)質(zhì)心等效加速度的影響規(guī)律[8-10]。對(duì)于連續(xù)梁,其固有頻率的影響因素有截面形狀、約條件、材料及分段桿長(zhǎng)等[11-14],李順才等[15]、黃志杰等[16]針對(duì)不同邊界條件下的連續(xù)梁振動(dòng)進(jìn)行了研究,提出了對(duì)應(yīng)的微分方程形式,對(duì)固有頻率以及振型等振動(dòng)特性進(jìn)行了理論求解,而器箭組合體實(shí)際為多分支梁結(jié)構(gòu),其動(dòng)態(tài)響應(yīng)的理論解還需進(jìn)一步研究。

    簡(jiǎn)化的梁動(dòng)力學(xué)模型可用于定性分析質(zhì)心等效加速度與不同參數(shù)間的關(guān)聯(lián),但當(dāng)前對(duì)質(zhì)心等效加速度較為準(zhǔn)確的計(jì)算還依賴于器箭動(dòng)力學(xué)耦合分析,該方法經(jīng)過(guò)逐漸發(fā)展,目前已形成成熟工具軟件[17-20]。郭峰等[21]研究了有效載荷基頻變化對(duì)火箭點(diǎn)火起飛時(shí)刻動(dòng)響應(yīng)的影響,結(jié)果表明基頻降低會(huì)顯著增加軸向過(guò)載,但未對(duì)橫向過(guò)載變化進(jìn)行研究,未涉及其他特性參數(shù)。

    本文基于分支梁動(dòng)力學(xué)理論,構(gòu)建了器箭組合體的簡(jiǎn)化模型,對(duì)橫向激勵(lì)作用下的橫向質(zhì)心等效加速度進(jìn)行了理論分析,基于理論解析式,辨識(shí)了主要影響因素,并研究了橫向質(zhì)心等效加速度與航天器質(zhì)量、質(zhì)心高度以及基頻等因素的相關(guān)性以及變化規(guī)律。同時(shí),采用耦合分析方法,系統(tǒng)獲取了跨聲速階段不同航天器參數(shù)組合下的橫向質(zhì)心等效加速度,結(jié)合理論分析規(guī)律,對(duì)橫向質(zhì)心等效加速度分布開(kāi)展進(jìn)一步分類統(tǒng)計(jì)分析,提出了一種可行的器箭界面動(dòng)態(tài)內(nèi)力精細(xì)化設(shè)計(jì)思路及方法。

    1 器箭界面動(dòng)態(tài)內(nèi)力影響因素分析

    火箭跨聲速飛行時(shí),激勵(lì)一般作用在整流罩前錐與柱段連接點(diǎn)位置。因此,可通過(guò)施加橫向單頻激勵(lì),基于簡(jiǎn)化的分支梁動(dòng)力學(xué)理論,推導(dǎo)橫向質(zhì)心等效加速度求解理論方程組,基于理論方程組對(duì)質(zhì)心等效加速度隨航天器質(zhì)量、質(zhì)心高度以及基頻等參數(shù)的變化規(guī)律進(jìn)行系統(tǒng)研究。

    1.1 器箭組合體分支梁模型

    器箭組合體結(jié)構(gòu)可以劃分為三個(gè)部分,分別為箭體結(jié)構(gòu)、整流罩以及航天器(有效載荷)??紤]到器箭組合體的長(zhǎng)細(xì)比較大,可以近似用梁模型動(dòng)力學(xué)基本理論進(jìn)行分析,梁模型如圖1。

    圖1 器箭組合體簡(jiǎn)化梁模型Fig.1 Simplified beam model of the combination of rocket and spacecraft

    其中1~3均為等直梁,1代表火箭箭體(不包含整流罩),2代表整流罩,3代表航天器,虛線框表征器箭界面,三段梁在該處具有共同的節(jié)點(diǎn)。在2號(hào)梁的右端點(diǎn)即整流罩頭部施加頻率為為ωF、幅值為F的單頻橫向外激勵(lì)模擬火箭在跨聲速飛行階段的橫向脈動(dòng)壓力。

    構(gòu)建器箭組合體簡(jiǎn)化模型后,基于梁振動(dòng)理論,即可對(duì)航天器質(zhì)心等效加速度ac表達(dá)形式進(jìn)行理論推導(dǎo)。

    1.2 航天器質(zhì)心等效加速度理論求解

    火箭飛行過(guò)程處于自由狀態(tài),因此首先考慮等直梁的自由振動(dòng)方程:

    (1)

    式中:yi表征i號(hào)梁的橫向位移;x代表i號(hào)梁上各點(diǎn)坐標(biāo);Mi,Li分別為i號(hào)梁的質(zhì)量及長(zhǎng)度;Hi代表i號(hào)梁的彎曲剛度(EJ)i。

    對(duì)式(1)頻響形式解進(jìn)行求解,即梁振動(dòng)頻率與外激勵(lì)頻率一致,也即認(rèn)為梁的橫向位移滿足如下形式:

    yi(x,t)=yi(x)exp(jωFt)

    (2)

    將式(2)中的橫向位移形式代入式(1)中,則可以求解得到頻響形式解為

    (3)

    對(duì)圖1中模擬器箭組合體的分支梁模型,取器箭連接面處為坐標(biāo)原點(diǎn),考慮到簡(jiǎn)化模型為線性系統(tǒng),取單位外激勵(lì)對(duì)該系統(tǒng)開(kāi)展分析。

    根據(jù)結(jié)構(gòu)力學(xué)理論,等直梁轉(zhuǎn)角、彎矩以及剪力可由橫向位移曲線確定,即:

    (4)

    對(duì)于器箭組合體分支梁模型,火箭(1號(hào)梁)及航天器(3號(hào)梁)為自由邊界,即彎矩與剪力為0;而整流罩由于頭部施加橫向外激勵(lì),因此邊界處彎矩為0,而剪力則等于施加的橫向外激勵(lì)。則動(dòng)力學(xué)方程邊界條件如下:

    (5)

    考慮到器箭界面處位移及轉(zhuǎn)角需要滿足變形協(xié)調(diào)條件,界面處可建立以下方程:

    (6)

    另,考慮界面處平衡條件,即連接面兩側(cè)彎矩及剪力需互相平衡,可建立以下方程:

    (7)

    以上方程構(gòu)成了器箭組合體系統(tǒng)三段等直梁的各4個(gè)非齊次線性方程組,將橫向位移代入方程,對(duì)以上方程組進(jìn)行求解,可得到器箭界面彎矩:

    (8)

    考慮橫向質(zhì)心等效加速度計(jì)算公式,記航天器質(zhì)量為m,質(zhì)心高度為h,橫向質(zhì)心加速度可通過(guò)界面彎矩?fù)Q算得到:

    (9)

    將航天器簡(jiǎn)化為等直梁,M3=m,L3=2h,最終可得到橫向質(zhì)心等效加速度表達(dá)式為

    (10)

    λ1~λ8為與航天器特性(質(zhì)量及質(zhì)心高度)無(wú)關(guān)的參數(shù),具體形式如下:

    (11)

    其中p0~p4與航天器基頻相關(guān),稱為頻率參數(shù),而ξ1~ξ8只與箭體結(jié)構(gòu)與外激勵(lì)頻率相關(guān),稱為箭體參數(shù),兩組參數(shù)的具體形式分別見(jiàn)式(12)和(13)。

    (12)

    (13)

    1.3 航天器質(zhì)心等效加速度變化規(guī)律研究

    根據(jù)動(dòng)力學(xué)理論分析得到的航天器橫向質(zhì)心等效加速度表達(dá)式(10),橫向質(zhì)心等效加速度與有效載荷、箭體結(jié)構(gòu)形式以及外激勵(lì)特性均有關(guān)系,考慮到對(duì)于特定型號(hào)的火箭,箭體結(jié)構(gòu)形式相對(duì)固定,同時(shí)外激勵(lì)形式經(jīng)過(guò)多次飛行驗(yàn)證,因此重點(diǎn)研究不同頻率外激勵(lì)作用下航天器質(zhì)量、質(zhì)心高度以及基頻對(duì)橫向質(zhì)心等效加速度的影響。本文選定某型整流罩、火箭及航天器等效參數(shù)開(kāi)展研究,具體參數(shù)見(jiàn)表1。

    表1 整流罩及火箭參數(shù)Table 1 Fairing and rocket parameters

    由于質(zhì)心等效加速度的變化趨勢(shì)十分復(fù)雜,以隨質(zhì)量的變化為例,不同外激勵(lì)和器箭組合體參數(shù)下,質(zhì)心等效加速度可能呈現(xiàn)完全不同的變化趨勢(shì)。因此本文定義了趨勢(shì)指數(shù)以更好地描述質(zhì)心等效加速度隨不同參數(shù)的整體變化規(guī)律。以下對(duì)趨勢(shì)指數(shù)進(jìn)行簡(jiǎn)要介紹。

    1.3.1航天器質(zhì)量對(duì)質(zhì)心等效加速度的影響

    為了更直觀地研究航天器質(zhì)量對(duì)橫向質(zhì)心等效加速度的影響,在選定若干航天器基頻及質(zhì)心高度參數(shù)組合下分析橫向質(zhì)心等效加速度隨質(zhì)量的變化規(guī)律,基頻分別為30 rad/s, 90 rad/s, 150 rad/s,對(duì)應(yīng)質(zhì)心高度分別為1 m, 3 m, 5 m。

    圖2~圖4為航天器不同質(zhì)心高度以及基頻組合下質(zhì)心等效加速度變化趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)量變化的二維分布圖。

    圖2 趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)量分布(ω=30 rad/s)Fig.2 Distribution of trend index on frequency of external excitation and mass (ω=30 rad/s)

    根據(jù)分析結(jié)果,對(duì)于各參數(shù)組合,盡管不同的外激勵(lì)頻率作用下,橫向質(zhì)心等效加速度隨航天器質(zhì)量的變化規(guī)律存在差異,但橫向質(zhì)心等效加速度關(guān)于外激勵(lì)頻率ωF與航天器質(zhì)量m的分布圖中均以白色區(qū)域占主要部分,占比均在65%以上,即橫向質(zhì)心等效加速度與航天器的質(zhì)量有明顯的相關(guān)性,隨著航天器質(zhì)量增加其有顯著的下降趨勢(shì),并且隨著質(zhì)心高度的增加,橫向質(zhì)心等效加速度隨航天器質(zhì)量而降低的趨勢(shì)逐漸增強(qiáng)。

    圖3 趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)量分布(ω=90 rad/s)Fig.3 Distribution of trend index on frequency of external excitation and mass (ω=90 rad/s)

    對(duì)比相同質(zhì)心高度、不同基頻的分布圖,白色及灰色區(qū)域占比變化很小,說(shuō)明基頻對(duì)橫向質(zhì)心等效加速度隨航天器質(zhì)量的變化趨勢(shì)影響并不明顯。

    1.3.2航天器質(zhì)心高度對(duì)質(zhì)心等效加速度的影響

    采用同樣方法,在給定的航天器基頻以及質(zhì)量組合下分析質(zhì)心高度對(duì)橫向質(zhì)心等效加速度的影響,基頻分別為30 rad/s, 90 rad/s及150 rad/s,對(duì)應(yīng)質(zhì)量分別為5000 kg, 15000 kg, 25000 kg。

    繪制橫向質(zhì)心等效加速度變化趨勢(shì)指數(shù)隨外激勵(lì)頻率以及航天器質(zhì)心高度變化分布圖,如圖5~圖7所示。可以看出,不同航天器質(zhì)量以及基頻組合下,橫向質(zhì)心等效加速度隨質(zhì)心高度的變化規(guī)律并不一致,總體呈現(xiàn)出隨著航天器質(zhì)量增加,橫向質(zhì)心等效加速度隨質(zhì)心高度變化減小趨勢(shì)占比增加,且不同的基頻下,該變化趨勢(shì)差別較大。ω=30 rad/s及ω=90 rad/s時(shí),不同質(zhì)量分布圖中,白色灰色區(qū)域占比變化相對(duì)較小;ω=150 rad/s時(shí),隨著質(zhì)量增加,白色區(qū)域占比增加十分明顯,即橫向質(zhì)心等效加速度隨質(zhì)心高度增加而減小趨勢(shì)顯著增強(qiáng)。

    圖5 趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)心高度分布(ω=30 rad/s)Fig.5 Distribution of trend index on frequency of external excitation and centroid location (ω=30 rad/s)

    圖6 趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)心高度分布(ω=90 rad/s)Fig.6 Distribution of trend index on frequency of external excitation and centroid location (ω=90 rad/s)

    圖7 趨勢(shì)指數(shù)隨外激勵(lì)頻率及質(zhì)心高度分布(ω=150 rad/s)Fig.7 Distribution of trend index on frequency of external excitation and centroid location (ω=150 rad/s)

    1.3.3航天器基頻對(duì)質(zhì)心等效加速度的影響

    采用同樣方法,在給定的航天器質(zhì)量及質(zhì)心高度組合參數(shù)下分析基頻對(duì)橫向質(zhì)心等效加速度的影響,質(zhì)心高度分別為1 m,3 m,5 m,對(duì)應(yīng)質(zhì)量分別為5000 kg,15000 kg,25000 kg。

    繪制橫向質(zhì)心等效加速度變化趨勢(shì)指數(shù)隨外激勵(lì)頻率和基頻的分布圖,如圖8~圖10所示。從橫向質(zhì)心等效加速度隨基頻的變化趨勢(shì)中可以看出,不同質(zhì)量及質(zhì)心高度的航天器,其分布圖中白色灰色區(qū)域占比沒(méi)有明顯變化規(guī)律,可認(rèn)為單個(gè)基頻參數(shù)對(duì)橫向質(zhì)心等效加速度的影響較小,兩者相關(guān)性較弱。

    圖8 趨勢(shì)指數(shù)隨外激勵(lì)頻率及基頻分布(h=1 m)Fig.8 Distribution of trend index on frequency of external excitation and fundamental frequency of spacecraft (h=1 m)

    圖9 趨勢(shì)指數(shù)隨外激勵(lì)頻率及基頻分布(h=3 m)Fig.9 Distribution of trend index on frequency of external excitation and fundamental frequency of spacecraft (h=3 m)

    圖10 趨勢(shì)指數(shù)隨外激勵(lì)頻率及基頻分布(h=5 m)Fig.10 Distribution of trend index on frequency of external excitation and fundamental frequency of spacecraft (h=5 m)

    2 器箭界面動(dòng)態(tài)內(nèi)力仿真及精細(xì)化設(shè)計(jì)

    前面基于分支梁動(dòng)力學(xué)理論,建立了質(zhì)心等效加速度的解析表達(dá)式,并基于解析表達(dá)式對(duì)質(zhì)心等效加速度隨主要參數(shù)變化趨勢(shì)進(jìn)行了研究,闡明了質(zhì)心等效加速度與不同參數(shù)之間的相關(guān)性。

    考慮理論解析時(shí)對(duì)器箭組合體和激勵(lì)形式進(jìn)行了簡(jiǎn)化,理論解僅用于分析橫向質(zhì)心等效加速度隨不同參數(shù)的變化趨勢(shì),無(wú)法得到較為準(zhǔn)確的橫向質(zhì)心等效加速度分布用于器箭精細(xì)化設(shè)計(jì),因此本章以某型火箭為例,建立器箭動(dòng)力學(xué)有限元模型。通過(guò)設(shè)計(jì)不同航天器的質(zhì)量、質(zhì)心位置以及基頻等特性參數(shù)組,耦合打靶仿真分析得到跨聲速飛行工況下真實(shí)的質(zhì)心等效加速度分布,并結(jié)合理論分析的規(guī)律,對(duì)質(zhì)心等效加速度數(shù)據(jù)開(kāi)展分類統(tǒng)計(jì)分析?;诮y(tǒng)計(jì)分析結(jié)果,完成器箭界面質(zhì)心等效加速度條件的精細(xì)化設(shè)計(jì)。

    2.1 橫向質(zhì)心等效加速度分布分析

    根據(jù)航天器合理參數(shù)范圍,確定參數(shù)區(qū)間見(jiàn)表2,在給定區(qū)間內(nèi)隨機(jī)生成航天器質(zhì)量、質(zhì)心高度以及基頻參數(shù)組合,并建立航天器有限元模型,其與火箭有限元模型組合后開(kāi)展器箭動(dòng)力學(xué)耦合分析,通過(guò)打靶仿真得到各參數(shù)組合下的質(zhì)心等效加速度。

    表2 航天器參數(shù)范圍Table 2 Range of spacecraft parameters

    根據(jù)理論分析規(guī)律,橫向質(zhì)心等效加速度與航天器質(zhì)量有明顯的相關(guān)性,同時(shí)質(zhì)心高度會(huì)對(duì)變化趨勢(shì)產(chǎn)生一定影響,航天器基頻對(duì)質(zhì)心等效加速度的影響具有較大的不確定性,因此可按照質(zhì)心高度對(duì)仿真結(jié)果進(jìn)行劃分,在每個(gè)區(qū)間內(nèi)對(duì)質(zhì)心等效加速度隨航天器質(zhì)量的變化規(guī)律進(jìn)行分析,獲得不同質(zhì)心高度航天器對(duì)應(yīng)的橫向質(zhì)心等效加速度隨有效載荷質(zhì)量的分布曲線如圖11所示。

    圖11 不同質(zhì)心高度航天器質(zhì)心等效加速度隨質(zhì)量分布Fig.11 Distribution of equivalence acceleration of spacecraft with different centroid locations

    根據(jù)上述分布規(guī)律,可以再按照航天器質(zhì)量將不同高度區(qū)間的橫向質(zhì)心等效加速度的變化曲線劃分為三個(gè)區(qū)間,其中,當(dāng)航天器質(zhì)量<1000 kg時(shí),橫向質(zhì)心等效加速度隨著航天器質(zhì)量的增加迅速降低;而當(dāng)航天器質(zhì)量1000 kg≤m<4000 kg時(shí),隨著航天器質(zhì)量的增加,橫向質(zhì)心加速度也會(huì)有明顯的降低,但是其下降速率要明顯低于m<1000 kg時(shí);當(dāng)航天器質(zhì)量m≥4000 kg時(shí),橫向質(zhì)心等效加速度隨質(zhì)量的下降曲線較平緩,該區(qū)間航天器質(zhì)量已不是影響橫向質(zhì)心等效加速度的主要因素,質(zhì)心高度對(duì)橫向質(zhì)心等效加速度的影響增強(qiáng)。

    2.2 器箭界面動(dòng)態(tài)內(nèi)力精細(xì)化設(shè)計(jì)

    根據(jù)上節(jié)的分析結(jié)論,針對(duì)任意火箭構(gòu)型,可開(kāi)展器箭界面動(dòng)態(tài)內(nèi)力的精細(xì)化設(shè)計(jì)。仍以上節(jié)中的火箭為對(duì)象,對(duì)器箭界面動(dòng)態(tài)內(nèi)力精細(xì)化設(shè)計(jì)過(guò)程進(jìn)行說(shuō)明。

    首先,按照質(zhì)心高度區(qū)間,對(duì)橫向質(zhì)心等效加速度精細(xì)化設(shè)計(jì),選取1.0 m≤h<2.0 m,2.0 m≤h<3.0 m以及3.0 m≤h<4.0 m三個(gè)區(qū)間進(jìn)行分類統(tǒng)計(jì),在各區(qū)間內(nèi)繪制橫向質(zhì)心等效加速度關(guān)于航天器質(zhì)量的分布圖,如圖12~圖14所示。

    然后,在不同質(zhì)心高度類別再按照質(zhì)量區(qū)間進(jìn)行分檔,可按照各區(qū)間打靶仿真得到的橫向質(zhì)心等效加速度最大值(當(dāng)區(qū)間子樣較多時(shí)取99%概率統(tǒng)計(jì)值),考慮動(dòng)力學(xué)分析不確定性系數(shù)(一般取1.25),分別確定各區(qū)間的橫向質(zhì)心等效加速度設(shè)計(jì)值。圖12~圖14中陰影柱狀分布即為不同質(zhì)量范圍的橫向質(zhì)心加速度設(shè)計(jì)值(工程應(yīng)用需確保子樣數(shù)足夠),而該型火箭跨聲速工況原設(shè)計(jì)值為1.5g。可以看出,當(dāng)衛(wèi)星質(zhì)量m≥2000 kg時(shí),且質(zhì)心高度h≥3.0 m時(shí),設(shè)計(jì)值減小至約0.7g,優(yōu)化效果達(dá)到50%以上。

    圖12 質(zhì)心等效加速度分布及精細(xì)化設(shè)計(jì)(1.0 m≤h<2.0 m)Fig.12 Distribution and fine design of equivalence acceleration of spacecraft (1.0 m≤h<2.0 m)

    圖14 質(zhì)心等效加速度分布及精細(xì)化設(shè)計(jì)(3.0 m≤h<4.0 m)Fig.14 Distribution and fine design of equivalence acceleration of spacecraft (3.0 m≤h<4.0 m)

    當(dāng)然,為了確保精細(xì)化設(shè)計(jì)條件的正確性,通常在航天器和火箭設(shè)計(jì)完成后,需根據(jù)真實(shí)狀態(tài)再次開(kāi)展動(dòng)力學(xué)耦合分析,驗(yàn)證精細(xì)化設(shè)計(jì)條件的覆蓋性,避免出現(xiàn)復(fù)雜因素導(dǎo)致的小概率不確定性。

    圖13 質(zhì)心等效加速度分布及精細(xì)化設(shè)計(jì)(2.0 m≤h<3.0 m)Fig.13 Distribution and fine design of equivalence acceleration of spacecraft (2.0 m≤h<3.0 m)

    3 結(jié) 論

    本文基于分支梁動(dòng)力學(xué)理論,推導(dǎo)了跨聲速飛行時(shí)器箭界面動(dòng)態(tài)內(nèi)力(橫向質(zhì)心等效加速度形式)的解析表達(dá)式,辨識(shí)出外激勵(lì)頻率、航天器質(zhì)量、質(zhì)心高度、基頻等主要影響參數(shù),并對(duì)橫向質(zhì)心等效加速度隨航天器質(zhì)量、質(zhì)心高度、基頻的變化趨勢(shì)開(kāi)展了系統(tǒng)研究。結(jié)合理論研究規(guī)律,采用器箭動(dòng)力學(xué)耦合仿真方法,通過(guò)隨機(jī)生成航天器質(zhì)量、質(zhì)心高度、基頻,得到了不同質(zhì)心高度、質(zhì)量下橫向質(zhì)心等效加速度的分布規(guī)律,并基于加速度分布的統(tǒng)計(jì)分析結(jié)果,制定了一種器箭界面動(dòng)態(tài)內(nèi)力的精細(xì)化設(shè)計(jì)方法,研究主要得到以下結(jié)論:

    1) 根據(jù)理論解析推導(dǎo),橫向質(zhì)心等效加速度與外激勵(lì)頻率、航天器質(zhì)量、質(zhì)心高度、基頻等參數(shù)均存在聯(lián)系。

    2) 橫向質(zhì)心等效加速度與航天器質(zhì)量間有明顯的相關(guān)性,隨著航天器質(zhì)量增加有明顯的下降趨勢(shì);質(zhì)心高度對(duì)橫向質(zhì)心等效加速度的影響規(guī)律在不同條件下存在較大差異,但整體上隨著質(zhì)量增加,橫向質(zhì)心等效加速度隨著質(zhì)心高度增加而下降的趨勢(shì)增強(qiáng),而基頻對(duì)橫向質(zhì)心等效加速度的影響并不明顯。

    3) 針對(duì)任意火箭,可通過(guò)開(kāi)展器箭耦合分析打靶仿真,并按照質(zhì)心高度及質(zhì)量區(qū)間進(jìn)行分類統(tǒng)計(jì),對(duì)器箭界面動(dòng)態(tài)內(nèi)力完成精細(xì)化設(shè)計(jì)。該精細(xì)化方法已在長(zhǎng)征五號(hào)B和長(zhǎng)征七號(hào)甲等新一代運(yùn)載火箭上成功應(yīng)用,顯著降低了器箭設(shè)計(jì)載荷。

    本文聚焦剛度分布較均勻的航天器界面載荷設(shè)計(jì)規(guī)律和優(yōu)化方法研究,對(duì)于剛度分布存在較大差異的航天器,還需增加振型等影響參數(shù)。此外,未來(lái)將結(jié)合理論和仿真方法,進(jìn)一步系統(tǒng)研究外激勵(lì)及器箭多個(gè)參數(shù)的耦合作用效應(yīng),補(bǔ)充完善器箭載荷設(shè)計(jì)優(yōu)化方法。

    猜你喜歡
    基頻質(zhì)心航天器
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    重型半掛汽車質(zhì)量與質(zhì)心位置估計(jì)
    語(yǔ)音同一認(rèn)定中音段長(zhǎng)度對(duì)基頻分析的影響
    基于時(shí)域的基頻感知語(yǔ)音分離方法?
    基于GNSS測(cè)量的天宮二號(hào)質(zhì)心確定
    橋面鋪裝層對(duì)中小跨徑橋梁基頻影響分析
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    45000kHz基頻晶體濾波器
    電子制作(2017年20期)2017-04-26 06:57:35
    色综合站精品国产| 国产av不卡久久| 国产中年淑女户外野战色| 国内精品久久久久精免费| 日日干狠狠操夜夜爽| 国产亚洲精品久久久久久毛片| 男女之事视频高清在线观看| 97碰自拍视频| 日日摸夜夜添夜夜添av毛片 | 欧美另类亚洲清纯唯美| 欧美zozozo另类| 国内揄拍国产精品人妻在线| 国产av麻豆久久久久久久| 桃红色精品国产亚洲av| 久久精品久久久久久噜噜老黄 | 亚洲av不卡在线观看| 免费观看人在逋| 亚洲真实伦在线观看| 黄色一级大片看看| 精品99又大又爽又粗少妇毛片 | 成年女人毛片免费观看观看9| 99九九线精品视频在线观看视频| 久久人人精品亚洲av| 搞女人的毛片| 国产69精品久久久久777片| 亚洲性夜色夜夜综合| 亚洲一级一片aⅴ在线观看| 日本-黄色视频高清免费观看| 淫妇啪啪啪对白视频| 午夜福利在线观看吧| 国产男人的电影天堂91| 99国产极品粉嫩在线观看| 小蜜桃在线观看免费完整版高清| 熟女人妻精品中文字幕| 日韩欧美免费精品| 搡老岳熟女国产| 黄色日韩在线| 欧美高清性xxxxhd video| 久久午夜福利片| 亚洲18禁久久av| 精品人妻偷拍中文字幕| 成人永久免费在线观看视频| 美女高潮的动态| 免费人成视频x8x8入口观看| 久久亚洲精品不卡| 日本-黄色视频高清免费观看| 啦啦啦观看免费观看视频高清| 亚洲欧美精品综合久久99| 韩国av在线不卡| 亚洲av熟女| 亚洲18禁久久av| 亚洲成人精品中文字幕电影| 婷婷精品国产亚洲av| 国产高清三级在线| 1024手机看黄色片| a级毛片a级免费在线| 久久亚洲精品不卡| 色哟哟哟哟哟哟| 国产麻豆成人av免费视频| 亚洲av.av天堂| 国产精品乱码一区二三区的特点| 我要看日韩黄色一级片| .国产精品久久| 成人鲁丝片一二三区免费| 美女高潮喷水抽搐中文字幕| 欧美日韩综合久久久久久 | 欧美丝袜亚洲另类 | 亚洲精品456在线播放app | 亚洲男人的天堂狠狠| 91久久精品国产一区二区三区| 三级国产精品欧美在线观看| 日本一本二区三区精品| 最后的刺客免费高清国语| 国产精品久久久久久精品电影| 国产亚洲精品综合一区在线观看| 日本黄色片子视频| 听说在线观看完整版免费高清| 日本黄大片高清| 真实男女啪啪啪动态图| 色播亚洲综合网| 草草在线视频免费看| 日日夜夜操网爽| 国产精品自产拍在线观看55亚洲| 亚洲在线观看片| 性色avwww在线观看| а√天堂www在线а√下载| 国产av麻豆久久久久久久| 在线天堂最新版资源| 直男gayav资源| 国产精品一区www在线观看 | 国内久久婷婷六月综合欲色啪| 久久精品国产亚洲网站| 精品久久久久久久末码| 嫩草影院精品99| 色吧在线观看| 91av网一区二区| 亚洲成人久久性| 亚洲av五月六月丁香网| 91在线观看av| 亚洲第一电影网av| 免费在线观看成人毛片| 禁无遮挡网站| 欧美zozozo另类| 亚洲最大成人av| 久久人人爽人人爽人人片va| 波多野结衣高清无吗| 亚洲人成网站在线播放欧美日韩| 久久久久久久久中文| 少妇的逼水好多| 可以在线观看毛片的网站| av黄色大香蕉| 国产精品久久久久久久久免| 国产欧美日韩一区二区精品| 啪啪无遮挡十八禁网站| 亚洲av第一区精品v没综合| 精品午夜福利在线看| 午夜福利欧美成人| 国产精品综合久久久久久久免费| 少妇熟女aⅴ在线视频| 少妇的逼水好多| av专区在线播放| 九色国产91popny在线| 99热这里只有精品一区| 久久久久性生活片| 在线免费十八禁| 久久久久久国产a免费观看| 国产成年人精品一区二区| 村上凉子中文字幕在线| 九九爱精品视频在线观看| 男女视频在线观看网站免费| 亚洲av.av天堂| 亚洲国产精品sss在线观看| 22中文网久久字幕| 国产极品精品免费视频能看的| 99久久九九国产精品国产免费| 看黄色毛片网站| 国产又黄又爽又无遮挡在线| 男人狂女人下面高潮的视频| 日韩亚洲欧美综合| 免费大片18禁| 丰满乱子伦码专区| 免费人成视频x8x8入口观看| 久久热精品热| 99视频精品全部免费 在线| 精品福利观看| av在线亚洲专区| 免费在线观看日本一区| 亚洲在线自拍视频| 国产av一区在线观看免费| 一级黄片播放器| 日本 欧美在线| 91久久精品国产一区二区成人| 欧美日韩乱码在线| 色哟哟·www| 一进一出好大好爽视频| 日韩中字成人| 国产一级毛片七仙女欲春2| 女同久久另类99精品国产91| 内射极品少妇av片p| 日韩欧美在线乱码| 成人精品一区二区免费| 日本在线视频免费播放| 亚洲av电影不卡..在线观看| 中文字幕精品亚洲无线码一区| 人妻少妇偷人精品九色| 成人无遮挡网站| 嫩草影院新地址| 免费看美女性在线毛片视频| 亚洲精品粉嫩美女一区| 日韩强制内射视频| 国产一区二区在线av高清观看| 成人国产麻豆网| 男人舔奶头视频| 三级国产精品欧美在线观看| 看十八女毛片水多多多| 亚洲,欧美,日韩| 在线国产一区二区在线| 国产色爽女视频免费观看| 国产精品久久久久久久久免| 国产精品美女特级片免费视频播放器| 美女cb高潮喷水在线观看| 久久久久久久精品吃奶| 亚洲综合色惰| 热99在线观看视频| 欧美又色又爽又黄视频| 两个人的视频大全免费| 国内少妇人妻偷人精品xxx网站| 欧美日韩综合久久久久久 | 亚洲,欧美,日韩| 亚洲性久久影院| 精品国内亚洲2022精品成人| 老师上课跳d突然被开到最大视频| 内地一区二区视频在线| 成熟少妇高潮喷水视频| 亚洲欧美激情综合另类| 99riav亚洲国产免费| 91久久精品国产一区二区成人| 国产伦人伦偷精品视频| 国产蜜桃级精品一区二区三区| 午夜福利18| 日本黄大片高清| 伦理电影大哥的女人| 欧美激情国产日韩精品一区| 国产精品99久久久久久久久| 国产黄片美女视频| 中文亚洲av片在线观看爽| 成人国产综合亚洲| 免费看光身美女| 婷婷亚洲欧美| 精品不卡国产一区二区三区| 一个人看的www免费观看视频| 欧美最新免费一区二区三区| 99视频精品全部免费 在线| 国产一级毛片七仙女欲春2| 亚洲一区二区三区色噜噜| 国产亚洲精品久久久久久毛片| 国产精品1区2区在线观看.| 97碰自拍视频| 搡老岳熟女国产| 搡女人真爽免费视频火全软件 | av专区在线播放| 国产久久久一区二区三区| 亚洲在线观看片| 赤兔流量卡办理| 亚洲,欧美,日韩| 伦理电影大哥的女人| 国语自产精品视频在线第100页| 国产高清不卡午夜福利| 婷婷精品国产亚洲av| 国产高清不卡午夜福利| 在线观看av片永久免费下载| 国产高潮美女av| 午夜影院日韩av| 神马国产精品三级电影在线观看| 免费在线观看影片大全网站| 国产乱人视频| 欧美日本亚洲视频在线播放| 美女免费视频网站| 精品一区二区免费观看| 亚洲精品成人久久久久久| 免费电影在线观看免费观看| 女生性感内裤真人,穿戴方法视频| 18禁在线播放成人免费| 免费搜索国产男女视频| 男人和女人高潮做爰伦理| 三级毛片av免费| 干丝袜人妻中文字幕| 久久精品国产亚洲av天美| 免费黄网站久久成人精品| 亚洲图色成人| 直男gayav资源| 国产色婷婷99| 欧美日韩瑟瑟在线播放| 亚洲av二区三区四区| 人妻夜夜爽99麻豆av| 国产三级中文精品| 久久国产乱子免费精品| 久久久成人免费电影| 欧美一区二区国产精品久久精品| 免费看美女性在线毛片视频| 国产成人影院久久av| 狂野欧美白嫩少妇大欣赏| 内地一区二区视频在线| 蜜桃亚洲精品一区二区三区| 亚洲欧美精品综合久久99| 乱码一卡2卡4卡精品| 日本一本二区三区精品| 午夜精品一区二区三区免费看| 欧美人与善性xxx| 亚洲专区国产一区二区| 麻豆久久精品国产亚洲av| 99国产极品粉嫩在线观看| 婷婷六月久久综合丁香| 88av欧美| 91在线精品国自产拍蜜月| 麻豆国产av国片精品| 热99re8久久精品国产| 亚洲美女搞黄在线观看 | 乱系列少妇在线播放| 久久草成人影院| 成人美女网站在线观看视频| 色5月婷婷丁香| 国产精品一区二区性色av| 成人国产麻豆网| 久久久久久国产a免费观看| 男女之事视频高清在线观看| 99久久无色码亚洲精品果冻| 日韩高清综合在线| xxxwww97欧美| 久久久久久久亚洲中文字幕| 久久久久久久久中文| 搞女人的毛片| 国内少妇人妻偷人精品xxx网站| 久久精品国产清高在天天线| 我要搜黄色片| 美女cb高潮喷水在线观看| 婷婷丁香在线五月| 亚洲欧美日韩高清在线视频| 在线观看舔阴道视频| 欧美三级亚洲精品| 日本五十路高清| 国产精品福利在线免费观看| 我的女老师完整版在线观看| 亚洲狠狠婷婷综合久久图片| www.www免费av| 干丝袜人妻中文字幕| 成人一区二区视频在线观看| 伦理电影大哥的女人| 国内少妇人妻偷人精品xxx网站| www.色视频.com| or卡值多少钱| 午夜福利视频1000在线观看| 亚洲人成伊人成综合网2020| 国产久久久一区二区三区| 成人永久免费在线观看视频| 国产高清视频在线播放一区| 五月玫瑰六月丁香| 一个人观看的视频www高清免费观看| 熟女人妻精品中文字幕| 日本三级黄在线观看| 黄色配什么色好看| 精品人妻1区二区| 亚洲av美国av| 天天躁日日操中文字幕| 久久久国产成人精品二区| 可以在线观看毛片的网站| 久久精品人妻少妇| 69av精品久久久久久| 欧美高清性xxxxhd video| 搡老熟女国产l中国老女人| 老熟妇仑乱视频hdxx| 悠悠久久av| 欧美绝顶高潮抽搐喷水| 老师上课跳d突然被开到最大视频| 久久久久久伊人网av| 精品午夜福利视频在线观看一区| 国产精品,欧美在线| 久久人人爽人人爽人人片va| 性色avwww在线观看| 久久精品国产99精品国产亚洲性色| 日日摸夜夜添夜夜添av毛片 | 校园人妻丝袜中文字幕| 美女 人体艺术 gogo| 国产黄色小视频在线观看| 99精品久久久久人妻精品| 啦啦啦韩国在线观看视频| 国产精品亚洲一级av第二区| 日韩中字成人| 成人午夜高清在线视频| 性欧美人与动物交配| 91狼人影院| 深夜a级毛片| a级毛片免费高清观看在线播放| 给我免费播放毛片高清在线观看| 久久人妻av系列| 九九在线视频观看精品| 国产免费av片在线观看野外av| 嫩草影院入口| 欧洲精品卡2卡3卡4卡5卡区| 免费看光身美女| 亚洲精品日韩av片在线观看| 亚洲无线在线观看| 国产成人aa在线观看| 成人一区二区视频在线观看| 在线观看av片永久免费下载| 综合色av麻豆| 夜夜夜夜夜久久久久| 免费av观看视频| 午夜福利18| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| .国产精品久久| 国产私拍福利视频在线观看| 俺也久久电影网| 成人综合一区亚洲| 精品久久久久久,| 99九九线精品视频在线观看视频| 一级黄片播放器| 久久久久久久久中文| 久久欧美精品欧美久久欧美| 琪琪午夜伦伦电影理论片6080| 午夜爱爱视频在线播放| 午夜免费男女啪啪视频观看 | 欧美绝顶高潮抽搐喷水| 免费看a级黄色片| 麻豆国产97在线/欧美| 免费观看的影片在线观看| 一本精品99久久精品77| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久精品电影| 午夜老司机福利剧场| 亚洲欧美精品综合久久99| 在线播放无遮挡| 欧美日韩综合久久久久久 | av在线天堂中文字幕| 精品无人区乱码1区二区| 亚洲午夜理论影院| 久久久久久九九精品二区国产| 精品乱码久久久久久99久播| 国产高清有码在线观看视频| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 日韩精品有码人妻一区| 国产精品久久久久久av不卡| 国产精品伦人一区二区| 色5月婷婷丁香| 国产一区二区三区av在线 | 老司机深夜福利视频在线观看| avwww免费| 啪啪无遮挡十八禁网站| 黄色丝袜av网址大全| 日韩中文字幕欧美一区二区| 波多野结衣高清作品| 国产亚洲精品久久久com| 免费一级毛片在线播放高清视频| 午夜免费成人在线视频| 天美传媒精品一区二区| 日本黄色片子视频| 黄色日韩在线| 国产91精品成人一区二区三区| 一本久久中文字幕| av专区在线播放| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 内地一区二区视频在线| 日韩中文字幕欧美一区二区| 成人特级av手机在线观看| 桃色一区二区三区在线观看| 女同久久另类99精品国产91| 国产欧美日韩精品亚洲av| 久久久久性生活片| 变态另类成人亚洲欧美熟女| 深夜精品福利| 国产单亲对白刺激| 日本成人三级电影网站| 成年女人永久免费观看视频| 久久草成人影院| 不卡一级毛片| 听说在线观看完整版免费高清| 校园春色视频在线观看| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看| 干丝袜人妻中文字幕| 成年女人看的毛片在线观看| 少妇裸体淫交视频免费看高清| 十八禁国产超污无遮挡网站| 亚洲美女黄片视频| 午夜免费激情av| 在线播放无遮挡| 久久久成人免费电影| 乱系列少妇在线播放| 婷婷亚洲欧美| 十八禁网站免费在线| 女的被弄到高潮叫床怎么办 | 天天躁日日操中文字幕| 窝窝影院91人妻| 日本与韩国留学比较| 在现免费观看毛片| 婷婷亚洲欧美| 免费高清视频大片| 男人狂女人下面高潮的视频| av在线蜜桃| 日日干狠狠操夜夜爽| 国内久久婷婷六月综合欲色啪| 久久精品国产鲁丝片午夜精品 | 免费搜索国产男女视频| 日韩大尺度精品在线看网址| 91av网一区二区| 春色校园在线视频观看| 欧美+日韩+精品| 国产精品一区二区三区四区久久| 亚洲av电影不卡..在线观看| 国产乱人视频| 天天躁日日操中文字幕| 国产精品野战在线观看| 国产精品av视频在线免费观看| 午夜激情欧美在线| 熟女电影av网| 我要搜黄色片| 国内精品久久久久精免费| 五月伊人婷婷丁香| 99精品久久久久人妻精品| 免费无遮挡裸体视频| 国产一区二区三区在线臀色熟女| 国产久久久一区二区三区| 日韩欧美精品v在线| 久久午夜亚洲精品久久| 欧美激情久久久久久爽电影| 久久久久久久午夜电影| 国产在视频线在精品| 97人妻精品一区二区三区麻豆| 好男人在线观看高清免费视频| 91在线观看av| 久久精品影院6| 亚洲图色成人| 亚洲真实伦在线观看| 精品一区二区三区视频在线| 欧美绝顶高潮抽搐喷水| 琪琪午夜伦伦电影理论片6080| 级片在线观看| 国产一区二区在线观看日韩| av女优亚洲男人天堂| 亚洲av二区三区四区| 综合色av麻豆| 欧美一级a爱片免费观看看| 男人舔女人下体高潮全视频| 少妇人妻精品综合一区二区 | 一级a爱片免费观看的视频| 最近在线观看免费完整版| 99久久九九国产精品国产免费| 乱人视频在线观看| 免费人成在线观看视频色| 久久精品国产清高在天天线| 91精品国产九色| 波多野结衣高清无吗| 真人一进一出gif抽搐免费| 国语自产精品视频在线第100页| 我要看日韩黄色一级片| 成人欧美大片| 久久久久久久久久黄片| 久久久久久久久久久丰满 | 麻豆国产97在线/欧美| 国产在视频线在精品| 一a级毛片在线观看| 免费看av在线观看网站| 午夜福利高清视频| 毛片女人毛片| 黄色日韩在线| 一区二区三区四区激情视频 | 国产成人福利小说| 精品久久久久久久久久久久久| 精品久久久久久,| 成人性生交大片免费视频hd| 精品不卡国产一区二区三区| 久久6这里有精品| 特级一级黄色大片| 深爱激情五月婷婷| 久久久久久久久中文| 国产av不卡久久| 淫妇啪啪啪对白视频| 91麻豆精品激情在线观看国产| 欧美成人免费av一区二区三区| 亚洲成人精品中文字幕电影| 久久精品国产亚洲网站| 99久久久亚洲精品蜜臀av| 99久久精品国产国产毛片| 中文亚洲av片在线观看爽| 69av精品久久久久久| 又黄又爽又免费观看的视频| 深爱激情五月婷婷| 国产精品亚洲美女久久久| 久久久精品欧美日韩精品| 18+在线观看网站| 色哟哟·www| 久久久久久久亚洲中文字幕| 精品一区二区免费观看| 九九爱精品视频在线观看| 午夜福利欧美成人| 国产69精品久久久久777片| 99国产极品粉嫩在线观看| 少妇人妻精品综合一区二区 | 精品久久久久久久久av| 十八禁国产超污无遮挡网站| 高清日韩中文字幕在线| 亚洲avbb在线观看| 久久人人精品亚洲av| 成年免费大片在线观看| 此物有八面人人有两片| 国产中年淑女户外野战色| 男女啪啪激烈高潮av片| 人妻少妇偷人精品九色| 国产免费一级a男人的天堂| 真人一进一出gif抽搐免费| 精品一区二区三区视频在线观看免费| 一本久久中文字幕| 夜夜看夜夜爽夜夜摸| 俺也久久电影网| 动漫黄色视频在线观看| 一本一本综合久久| 国内少妇人妻偷人精品xxx网站| 无人区码免费观看不卡| 两个人的视频大全免费| 国产欧美日韩一区二区精品| 色哟哟哟哟哟哟| 男插女下体视频免费在线播放| 亚洲人成网站在线播放欧美日韩| 亚洲精品久久国产高清桃花| 一区二区三区激情视频| 美女cb高潮喷水在线观看| 成人特级黄色片久久久久久久| 亚洲va在线va天堂va国产| 又紧又爽又黄一区二区| 日韩高清综合在线| 成年女人看的毛片在线观看| 久9热在线精品视频| 22中文网久久字幕| 亚洲国产色片| 亚洲 国产 在线| 91麻豆av在线| 日本 av在线| 尤物成人国产欧美一区二区三区| 日本精品一区二区三区蜜桃| 国产精品电影一区二区三区| 丰满乱子伦码专区| 韩国av一区二区三区四区| 极品教师在线免费播放| 久久久久免费精品人妻一区二区| 国产伦人伦偷精品视频| 国内精品一区二区在线观看| 岛国在线免费视频观看| 99久久九九国产精品国产免费| 亚洲五月天丁香| 日韩中文字幕欧美一区二区| 最近在线观看免费完整版| 久久精品人妻少妇|