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

    深空探測器模塊化結(jié)構(gòu)動力學(xué)研究

    2016-10-20 03:41:04于登云張玉花褚英志李昊王建煒杜冬
    深空探測學(xué)報 2016年3期
    關(guān)鍵詞:子結(jié)構(gòu)航天器探測器

    于登云,張玉花,褚英志,李昊,王建煒,杜冬

    (1. 中國航天科技集團公司,北京 100048;2. 上海航天技術(shù)研究院,上海 201109;3. 上海衛(wèi)星工程研究所,上海 200240)

    深空探測器模塊化結(jié)構(gòu)動力學(xué)研究

    于登云1,張玉花2,褚英志2,李昊3,王建煒3,杜冬3

    (1. 中國航天科技集團公司,北京 100048;2. 上海航天技術(shù)研究院,上海 201109;3. 上海衛(wèi)星工程研究所,上海 200240)

    隨著深空探測器技術(shù)的不斷發(fā)展,深空探測器設(shè)計逐漸轉(zhuǎn)向模塊化和標(biāo)準(zhǔn)化。串聯(lián)結(jié)構(gòu)是模塊化深空探測器常見并較為合理的結(jié)構(gòu)形式,而如何將運載對探測器整體模態(tài)頻率的指標(biāo)要求合理轉(zhuǎn)化為對探測器各個模塊結(jié)構(gòu)的設(shè)計要求是當(dāng)前工作中的一項挑戰(zhàn)。根據(jù)子結(jié)構(gòu)綜合模態(tài)法并基于固定邊界法和懸臂梁理論,提出了探測器子結(jié)構(gòu)模塊動力學(xué)模態(tài)指標(biāo)的分配方法。采用3種子結(jié)構(gòu)模塊極限處理方法,用全解析方式推導(dǎo)了子結(jié)構(gòu)模塊的基頻表達(dá)式,并分析了子結(jié)構(gòu)模塊的動力學(xué)特性與探測器整體結(jié)構(gòu)動力學(xué)特性之間的關(guān)系。提出的子結(jié)構(gòu)模塊的動力學(xué)指標(biāo)分解方法操作簡單,精度較高,可為串聯(lián)式模塊化深空探測器的設(shè)計提供關(guān)鍵性的指導(dǎo)。

    深空探測器;模塊化;子結(jié)構(gòu);動力學(xué);指標(biāo)分解

    引用格式:于登云,張玉花,褚英志,等. 深空探測器模塊化結(jié)構(gòu)動力學(xué)研究[J]. 深空探測學(xué)報,2016,3(3):268-274.

    Reference format: Yu D Y,Zhang Y H,Chu Y Z,et al. Structure dynamics analysis of modularized deep space detector[J]. Journal of Deep Space Exploration,2016,3(3):268-274.

    0 引 言

    隨著航天技術(shù)的不斷發(fā)展,航天器結(jié)構(gòu)設(shè)計逐漸轉(zhuǎn)向模塊化、標(biāo)準(zhǔn)化。在深空探測領(lǐng)域,執(zhí)行深空探測任務(wù)往往需要同時具備軌道轉(zhuǎn)移功能、軌道進入功能和科學(xué)試驗功能,所以通常設(shè)置不同獨立模塊完成相應(yīng)功能。受運載包絡(luò)空間的限制,模塊化串聯(lián)結(jié)構(gòu)成為深空探測航天器的典型形式。例如1989年美國發(fā)射的“麥哲倫號”金星探測器,俄羅斯與中國聯(lián)合研制的“福布斯”與“螢火1號”火星探測器組合體,以及擬于2022年發(fā)射的金星氣候探測器都是典型的模塊化串聯(lián)結(jié)構(gòu),如圖1所示。

    航天器結(jié)構(gòu)動力學(xué)的基本任務(wù)是分析研究航天器結(jié)構(gòu)的固有特性(包括模態(tài)頻率和模態(tài)振型),以及在外載激勵下的動力學(xué)響應(yīng),包括結(jié)構(gòu)模態(tài)分析、動靜強度計算、剛度分析、響應(yīng)分析與結(jié)構(gòu)優(yōu)化等。航天器結(jié)構(gòu)動力學(xué)分析的目的是:1)為航天器結(jié)構(gòu)與機構(gòu)優(yōu)化設(shè)計與試驗提供預(yù)示和設(shè)計依據(jù)。2)為運載火箭提供星箭模態(tài)綜合與試驗的全星結(jié)構(gòu)動力學(xué)特性。3)為全星與結(jié)構(gòu)振動有關(guān)的各類耦合動力學(xué)分析提供附件級乃至全星級結(jié)構(gòu)模態(tài)參數(shù)等[2]。

    航天對結(jié)構(gòu)設(shè)計各方面的要求都極為苛刻,其中模態(tài)參數(shù)等是航天器設(shè)計極為關(guān)鍵的參數(shù),這些參數(shù)的精確性與可靠性決定了航天器控制系統(tǒng)是否能正常工作,這些參數(shù)的錯誤會造成控制的錯誤,使航天器因無法正常飛行而失敗。因此航天器結(jié)構(gòu)動力學(xué)的研究具有重要意義。在工程應(yīng)用中,對于復(fù)雜的航天器結(jié)構(gòu),通??梢苑纸獬扇舾蓚€比較簡單的子結(jié)構(gòu)。復(fù)雜結(jié)構(gòu)的子結(jié)構(gòu)分析法通過對各個子結(jié)構(gòu)的動態(tài)特性進行綜合,并利用子結(jié)構(gòu)間的邊界條件,最終得到整個結(jié)構(gòu)的動態(tài)特性。

    圖1 模塊化串聯(lián)結(jié)構(gòu)深空探測器Fig.1 Modularized series-connected deep space detectors

    1 航天器結(jié)構(gòu)動力學(xué)子結(jié)構(gòu)模態(tài)分析方法

    自20世紀(jì)60年代初Hurty(1960)[3]和Gladwell(1964)[4]奠定模態(tài)綜合技術(shù)以來,子結(jié)構(gòu)方法[5-7]已廣泛應(yīng)用于航天航空和各種大型工程領(lǐng)域,這是一種復(fù)雜結(jié)構(gòu)建模與分析的有效方法。采用這種方法通過模態(tài)坐標(biāo)變換可以把結(jié)構(gòu)動力學(xué)問題化為縮聚自由度的問題,從而大大簡化計算,提高分析效率。動態(tài)子結(jié)構(gòu)方法[8-9]以瑞利-李茲法(Rayleigh-Ritz)為理論基礎(chǔ),并在20世紀(jì)60年代由Hurty和Gladwell等人提出模態(tài)綜合技術(shù)后蓬勃發(fā)展,經(jīng)Craig和Bampton等學(xué)者改進完善后,目前已廣泛用于航空航天等工程領(lǐng)域。

    動態(tài)子結(jié)構(gòu)方法的基本思想是按照工程的觀點或結(jié)構(gòu)的幾何輪廓,在遵循某些原則要求的基礎(chǔ)上,把完整的大型復(fù)雜結(jié)構(gòu)人為地抽象成若干個子結(jié)構(gòu),首先對各子結(jié)構(gòu)進行凝聚得到其主要模態(tài)信息,然后再綜合得到總體結(jié)構(gòu)的動態(tài)特性。動態(tài)子結(jié)構(gòu)法在大大降低結(jié)構(gòu)自由度的基礎(chǔ)上,保留了結(jié)構(gòu)的動態(tài)特性信息。一方面,能夠以較少的時間和資源成本解決大型結(jié)構(gòu)系統(tǒng)的動力學(xué)計算問題;另一方面,在模態(tài)坐標(biāo)下給出有限元模型的動態(tài)特性信息,能夠屏蔽結(jié)構(gòu)設(shè)計參數(shù)等技術(shù)秘密,便于企業(yè)間和國際間合作的開展。

    動態(tài)子結(jié)構(gòu)模態(tài)綜合法按照對界面的處理方法,分為自由界面法、固定界面法和混合和綜合界面法。自由界面法將子結(jié)構(gòu)的界面自由度完全放松,利用子結(jié)構(gòu)的模態(tài)來構(gòu)造變換,通過界面協(xié)調(diào)條件來得到整體結(jié)構(gòu)的縮聚系統(tǒng)。自由界面法由于綜合方程不包含界面自由度,所以自產(chǎn)生之初一直備受關(guān)注。Hou(1969)[10]和Goldman(1969)[11]于1969年分別提出了自由界面法。由于子結(jié)構(gòu)的低階主模態(tài)只對總體結(jié)構(gòu)的低階模態(tài)有較大影響,因此Hou和Goldman只采用低階主模態(tài)為基礎(chǔ)來描述子結(jié)構(gòu)的運動。初始的基于低階模態(tài)的子結(jié)構(gòu)法由于未考慮高階模態(tài),計算精度較差。在自由界面法的后續(xù)發(fā)展中,Macneal(1971)[12]、Rubin(1975)等[13]考慮了高階模態(tài)的影響,大大提高了自由界面法的計算精度。固定界面法通過在子結(jié)構(gòu)的連接面上施加固定約束,以此結(jié)構(gòu)的主模態(tài)和約束模態(tài)構(gòu)造整體結(jié)構(gòu)的綜合模態(tài)。Hurty(1960)[14]首先提出了模態(tài)坐標(biāo)和模態(tài)綜合法的概念,奠定了固定界面法的基礎(chǔ)。Craig和Bampton對Hurty的方法進行改進,認(rèn)為應(yīng)將復(fù)雜結(jié)構(gòu)的模態(tài)分為主模態(tài)和約束模態(tài)兩部分。改進的方法擁有綜合精度高和編制程序方便的優(yōu)點,因此被廣泛使用[15]。Craig和Chang(1977)[16]針對固定界面法自由度較多的缺點,提出了3種自由度減縮方法,分別是主模態(tài)坐標(biāo)減縮、Guyan減縮和Ritz減縮?;旌辖缑娣ń橛谧杂山缑娣ê凸潭ń缑娣ㄖg,兼具兩種方法的特點,可施加約束或解除約束來生成彼此獨立的子結(jié)構(gòu),并分別予以綜合。Benfield和Hruda(1971)[15]提出綜合界面法。該方法按照某一界面固定而對應(yīng)的另一界面自由的原則,將一個復(fù)雜結(jié)構(gòu)分離成若干主從結(jié)構(gòu)和從子結(jié)構(gòu)。

    隨著航天器研制過程中部門間、企業(yè)間乃至國際間合作的增多,航天器各部件有限元模型的交互日益頻繁。但考慮到保護其自身技術(shù)的需要,企業(yè)常常不愿將反映大量結(jié)構(gòu)設(shè)計信息的有限元模型直接提供給合作方。迫切需要將有限元物理模型進行適當(dāng)技術(shù)處理,在保留模型力學(xué)特性的情況下,屏蔽具體的設(shè)計信息。動態(tài)子結(jié)構(gòu)分析方法的出現(xiàn),能較好地解決上述問題。

    2 基于子結(jié)構(gòu)模態(tài)法的航天器模塊模態(tài)指標(biāo)設(shè)計

    隨著航天領(lǐng)域內(nèi)的合作與分工日益復(fù)雜化,模塊化航天器各個模塊又往往分屬于不同研制單位,如果能夠?qū)⒄髂B(tài)要求合理分配成各個模塊的結(jié)構(gòu)設(shè)計指標(biāo),則對于任務(wù)的分解和總體集成極富應(yīng)用價值。因此,模塊化航天器的設(shè)計過程中,需首先確定不同模塊的各自模態(tài)頻率與航天整體結(jié)構(gòu)模態(tài)頻率之間的定量化關(guān)系。在運載限定整器模態(tài)頻率下限的情況下,合理分配各自的設(shè)計指標(biāo)對于各個模塊的結(jié)構(gòu)設(shè)計具有重要意義。

    懸臂梁結(jié)構(gòu)與串聯(lián)模塊化探測器的發(fā)射狀態(tài)極為相似,本研究基于懸臂梁假設(shè),探索串聯(lián)結(jié)構(gòu)基頻的快速估算方法,可為總體指標(biāo)分解與集成提供有益的技術(shù)支持。目前,航天領(lǐng)域內(nèi)對于深空探測器結(jié)構(gòu)動力學(xué)設(shè)計的慣用做法是由各個研制單位提供各自負(fù)責(zé)模塊的有限元模型,經(jīng)總體復(fù)核整器模態(tài)特性后作為實物研制的參考依據(jù)。這種設(shè)計方法的缺點是,以有限元模型作為參考依據(jù)可操作性較差,尤其難于通過各個模塊單獨的力學(xué)試驗考核設(shè)計的合理性,研制任務(wù)書中無法提煉出關(guān)鍵性的技術(shù)指標(biāo),實際上也就造成了盲目設(shè)計,一旦出現(xiàn)整器模態(tài)頻率不達(dá)標(biāo),往往不能明確劃分各方的責(zé)任界限。

    2.1串聯(lián)結(jié)構(gòu)理論建模

    針對串聯(lián)型結(jié)構(gòu),為簡化分析,得到更為普適的解析解,將串聯(lián)結(jié)構(gòu)的基礎(chǔ)模塊和附加模塊假設(shè)為串聯(lián)的兩段懸臂梁,如圖2所示。

    圖2中,梁Ⅰ代表基礎(chǔ)模塊,梁Ⅱ代表附加模塊。兩段梁的線密度分別為ρ1和ρ2,抗彎剛度分別為E1I1和E2I2,長度分別為a和b。懸臂梁總長度為l,根據(jù)圖中所示應(yīng)有

    圖2 懸臂梁模型示意圖Fig.2 Sketch map of cantilever model

    按照經(jīng)典分析方法,可假定懸臂梁基頻振型與懸臂梁受慣性過載產(chǎn)生的靜變形曲線完全一致。設(shè)懸臂梁的末端位移為y0。不難計算懸臂梁的靜變形曲線為(推導(dǎo)過程較為復(fù)雜,限于篇幅,此處僅列出結(jié)果)

    其中

    根據(jù)簡單梁假設(shè),懸臂梁最大彈性勢能為

    將式(2)~(5)代入式(6)、(7),并利用瑞利原理,令

    式(9)、(10)代表了一般情況下串聯(lián)結(jié)構(gòu)基頻的解析解。為檢驗解析解的正確性,取一種最簡單的情況,即兩端懸臂梁的抗彎剛度和線密度完全相同,此時

    將式(15)代入式(11)、(12)和式(9),最終得到

    或者,令梁Ⅰ的抗彎剛度E1I1無窮大,則懸臂梁退化為梁Ⅱ,此時

    將式(17)代入式(9),可得

    長度為L,抗彎剛度為EI,線密度為ρ的懸臂梁基頻之精確解為

    式(16)、(18)和式(19)的誤差在0.43%之內(nèi),可見上述利用假設(shè)振型并通過瑞利原理計算基頻的近似方法精度是相當(dāng)高的。

    2.2模型研究

    從力學(xué)角度作定性考慮,圖2所示串聯(lián)式懸臂結(jié)構(gòu)的一彎模態(tài)中,梁Ⅰ作為基礎(chǔ)模塊大體上是以剛度單元參與作用的,梁Ⅱ作為附加模塊大體上以質(zhì)量單元參與作用。為便于比較,將模型作3種極限處理。

    將式(21)代入式(11)、(12)和式(10),最終得到

    此時梁Ⅰ(基礎(chǔ)模塊)僅具有剛度特性(既有線剛度,又有角剛度)。

    將式(24)代入式(11)、(12)和式(10),最終得到

    此時梁Ⅱ(附加模塊)僅具有質(zhì)量特性(包括質(zhì)量、質(zhì)心以及慣量3個特性)。

    將式(28)代入式(11)、(12)和式(10),最終得到

    此時梁Ⅰ(基礎(chǔ)模塊)僅具有剛度特性,梁Ⅱ(附加模塊)僅具有質(zhì)量特性。

    2.3比較與分析

    為顯示真實情況(見式(10))與3種極限處理方法(見式(22)、(25)、(29))之間的差異,取4個頻率解析式中無量綱部分,即

    并考慮基礎(chǔ)模塊(梁Ⅰ)與附加模塊(梁Ⅱ)剛度、質(zhì)量密度的差異(即不同的α、β),列于圖3中。橫軸為,即基礎(chǔ)模塊高度所占的比例??紤]到工程中可能的情況,分別列舉了α取10、1、0.1和β取10、1、0.1的組合,總計9種狀態(tài)。

    對圖3分析如下。

    1)圖3(a)~(h)顯示,僅考慮基礎(chǔ)模塊(梁Ⅰ,下同)的剛度特性,忽略質(zhì)量特性,同時僅考慮附加模塊(梁Ⅱ,下同)的質(zhì)量特性,忽略剛度特性,這種處理方法與真實情況誤差最大,不推薦作為串聯(lián)結(jié)構(gòu)基頻估算方法。不過圖3(i)表明,當(dāng)且僅當(dāng)基礎(chǔ)模塊剛度較大,同時附加模塊質(zhì)量較大時,該方法計算誤差較小。還應(yīng)當(dāng)看到,這種估算方法得到的無量綱曲線具有平坦的“谷底”,相對而言,該區(qū)域與真實情況接近——在對模塊更詳細(xì)的模態(tài)信息無從知曉的情況下,不妨將“谷底”區(qū)域作為初步的近似計算手段。

    2)僅考慮基礎(chǔ)模塊的剛度特性,忽略其質(zhì)量特性,如果基礎(chǔ)模塊高度占整體結(jié)構(gòu)高度比例不太大,該方法與真實情況誤差較小,見圖3(d)~(f)。如果基礎(chǔ)模塊剛度較附加模塊大很多,則即便基礎(chǔ)模塊高度比例偏大,該方法也可能有效,見圖3(a)~(c)。

    圖3 串聯(lián)式結(jié)構(gòu)模態(tài)頻率計算方法比較Fig.3 Comparison of model frequency prediction using different method

    3)僅考慮附加模塊的質(zhì)量特性,忽略其剛度特性,如果附加模塊高度占整體結(jié)構(gòu)高度比例不太大,該方法與真實情況誤差較小,見圖3(a)、(b)、(d)、(e)、(g)、(h)。如果基礎(chǔ)模塊質(zhì)量較附加模塊大,那么即便附加模塊高度比例偏大,該方法也可能有效,見圖3(g)~(i)。

    4)從圖3中總結(jié)得到

    也就是說,可以將考慮基礎(chǔ)模塊的剛度特性并忽略其質(zhì)量特性的極限情況、考慮附加模塊的質(zhì)量特性并忽略其剛度特性的極限情況作為串聯(lián)結(jié)構(gòu)模態(tài)頻率的簡單估算方法。真實頻率就是兩種極限處理方法的最小下界,其誤差不超過5%。

    上述分析的第4條可作為串聯(lián)式結(jié)構(gòu)模態(tài)頻率總體指標(biāo)分解的重要依據(jù),具體操作方法如下。

    首先,在對運載提出的模態(tài)頻率指標(biāo)分解之前,應(yīng)該估算附加模塊的質(zhì)量特性(質(zhì)量、質(zhì)心和慣量)以及整星外形基本尺寸,利用單自由度彈簧振子原理,根據(jù)質(zhì)量特性和運載模態(tài)頻率要求估算基礎(chǔ)模塊與附加模塊連接界面的剛度矩陣,此剛度矩陣為基礎(chǔ)模塊整體的靜剛度參數(shù),包含線剛度、角剛度以及耦合剛度。估算的附加模塊質(zhì)量特性與基礎(chǔ)模塊剛度矩陣是對研制方的輸入內(nèi)容;然后,對研制方提出要求。

    對基礎(chǔ)模塊結(jié)構(gòu)研制方的要求如下。

    1)基礎(chǔ)模塊在頂部連接了質(zhì)量特性與附加模塊完全相同的剛體之后,基頻應(yīng)高于運載要求。

    2)基礎(chǔ)模塊結(jié)構(gòu)與附加模塊連接界面處的靜剛度矩陣元素絕對值必須大于估算的靜剛度矩陣。

    對附加模塊結(jié)構(gòu)研制方的要求如下。

    1)附加模塊在底部連接了剛度矩陣與估算剛度矩陣相同的彈性單元(無質(zhì)量)之后,基頻應(yīng)高于運載要求。

    2)附加模塊的質(zhì)量、質(zhì)心、慣量等數(shù)據(jù)與估算質(zhì)量特性之間的誤差不得大于5%。

    如果基礎(chǔ)模塊結(jié)構(gòu)研制方能夠滿足總體要求,則Λ3滿足基頻要求。如果附加模塊研制方能夠滿足總體要求,則Λ2也滿足基頻要求。根據(jù)式(34),Λ必然也滿足運載要求,其誤差在5%以內(nèi)。為保證指標(biāo)可靠性,總體部門在實際下達(dá)頻率指標(biāo)要求時可適當(dāng)留有一定余量。

    針對具有多個模塊的串聯(lián)結(jié)構(gòu),該方法也可采用遞歸方式予以推廣。以三模塊串聯(lián)(即上、中、下模塊串聯(lián))為例,可先定義上部模塊為附加模塊,中、下部模塊整體作為基礎(chǔ)模塊,完成雙串聯(lián)結(jié)構(gòu)頻率指標(biāo)分配后,將得到的基礎(chǔ)模塊頻率指標(biāo)作為中下部模塊的總體指標(biāo),進行再分配,得到中下部兩模塊各自的頻率指標(biāo)。必須指出,多次遞歸式的誤差可能將增加。

    3 結(jié) 論

    本文對模塊化航天器設(shè)計中的子結(jié)構(gòu)動力學(xué)分析方法進行研究,并基于子結(jié)構(gòu)模態(tài)綜合法,提出了一種對兩模塊串聯(lián)航天器結(jié)構(gòu)各個模塊指標(biāo)進行分配的設(shè)計方法,即利用兩種極限情況擬合真實結(jié)構(gòu)的模態(tài)頻率,其誤差在5%以內(nèi)。該方法可以為總體分配探測器兩個模塊結(jié)構(gòu)設(shè)計指標(biāo),或者分配星體結(jié)構(gòu)與大型有效載荷結(jié)構(gòu)設(shè)計指標(biāo)提供關(guān)鍵性支持。

    [1]張偉,方寶東,成玫,等. 空間飛行器分離式構(gòu)型設(shè)計[J]. 上海航天,2013,30(1):1-7. Zhang W,F(xiàn)ang B D,Cheng M,et al. Design of separation configuration for spacecraft[J]. Aerospace Shanghai,2013,30(1):1-7.

    [2]黃文虎,曹登慶,韓增堯. 航天器動力學(xué)與控制的研究進展與展望[J]. 力學(xué)進展,42(4):68-76. Huang W H,Cao D Q,Han Z Y. Advances and trends in dynamics and control of spacecrafts[J]. Advances in Mechanics,42(4):68-76.

    [3]Hurty W C. Vibrations of structural systems by component mode synthesis[J]. Journal of the Engineering Mechanics Division,ASCE,1960,86(4):51-70.

    [4]Gladwell G M L. Branch mode analysis of vibrating systems[J]. Journal of Sound and Vibration,1964,1(1):41-59.

    [5]王文亮,杜作潤. 結(jié)構(gòu)振動與動態(tài)子結(jié)構(gòu)方法[M]. 上海:復(fù)旦大學(xué)出版社,1985. Wang W L,Du Z R. Structural vibration and dynamic substructure method[M]. Shanghai:Fudan University Press,1985.

    [6]殷學(xué)剛,陳淮,塞開林. 結(jié)構(gòu)振動分析的子結(jié)構(gòu)方法[M]. 北京:中國鐵道出版社,1991. Yin X G,Chen H,Sai K L. Substructure method for structural vibration analysis[M]. Beijing:China Railway Publishing House,1991.

    [7]鄭兆昌. 復(fù)雜結(jié)構(gòu)振動研究的模態(tài)綜合技術(shù)[J]. 振動與沖擊,1982(1):28-36. Zheng Z C. The study on vibration of complex structure systems by component mode synthesis techniques[J]. Journal of Vibration and Shock,1982 (1):28-36.

    [8]宋景濤,方明霞. 模態(tài)綜合法在ANSYS中的應(yīng)用[J]. 計算機輔助工程,2007,16(3):145-148. Song J T,F(xiàn)ang M X. Application of component mode synthesis in ansys[J]. Computer Aided Engineering,2007,16(3):145-148.

    [9]王永巖.動態(tài)子結(jié)構(gòu)方法理論及應(yīng)用[M]. 北京:科學(xué)出版社,1999. Wang Y Y. Theory and application of dynamic substructure method[M]. Beijing:Science Press,1999.

    [10]Hou S N. Review of modal synthesis techniques and a newapproach[J]. Shock Vib. bull ,1969,40(4):25-39.

    [11]Goldman R L. Vibration analysis by dynamic partitioning[J]. AIAA J,1969,7(6):1152-1154.

    [12]MacNeal R H. A hybrid method of component mode synthesis[J]. Computers Structures,1971,1(4):581-601.

    [13]Rubin S. Improved component-mode representation for structure dynamic analysis[J]. AIAA J,1975,13(8):995-1006.

    [14]Hurty W C. Vibration of structure by component mode synthesis[J]. J. Engr. Div.,ASCE,1960,86:51-59.

    [15]杜飛平,譚永華,陳建華,等. 航天器子結(jié)構(gòu)模態(tài)綜合法研究現(xiàn)狀及進展[J]. 力學(xué)進展,2010,36(3):39-62. Du F P,Tan Y H,Chen J H,et al. Status and progress on substructural modal synthesis techniques for space vehicls[J] Advances in Mechanics,2010,36(3):39-62.

    [16]Craig R R,Chang C J. On the use of attachment modes in substructure coupling for dynamic analysis[C]. AIAA/ASME 18th Struct. Dyn. and Materials. Conf.,SanDiego,CA:1977.

    通信地址:北京市海淀區(qū)阜成路16號航天科技大廈(100048)

    電話:(010)68370292

    E-mail:yudyun@sina.com

    Structure Dynamics Analysis of Modularized Deep Space Detector

    YU Dengyun1,ZHANG Yuhua2,CHU Yingzhi2,LI Hao3,WANG Jianwei3,DU Dong3
    (1. China Aerospace Science And Technology Corporation,Beijing 100048,China;2. Shanghai Academy Of Spaceflight Technology,Shanghai 201109,China;3. Shanghai Institute Of Satellite Engineering,Shanghai 200240,China)

    Modularized and standardized concepts of deeps pace detectors are increasing along with the development of deep space exploration technique. Series connection between different sub-modules of deep space detectors are reasonable and popular,and it is a challenge to propose a reasonable dynamic response index for a single sub-module according to the dynamic response index of the overall series-connected structure. Based on the substructure modal synthesis method and the cantilever beam theory, It proposes a method for assigning design index of dynamic response for sub-modules of series-connected structures. The equation of fundamental frequencies of sub-modules are derived analytically based on three consumptions,and the relation between the dynamics of sub-modules and that of overall series-connected structure are discussed. The proposed dynamics index assignment method shows high accuracy and can provide critical instructions for the design of series-connected modularized deep space detectors.

    deep space detector;modularize;sub-structure;dynamics;index assignment

    V476

    A

    2095-7777(2016)03-0268-7

    10.15982/j.issn.2095-7777.2016.03.011

    于登云(1961- ),男,研究員,博士生導(dǎo)師,中國航天科技集團公司科技委副主任,中國探月工程副總設(shè)計師,主要研究方向:航天器動力學(xué)與總體設(shè)計技術(shù)。

    [責(zé)任編輯:楊曉燕]

    2016-04-20;

    2016-05-21

    猜你喜歡
    子結(jié)構(gòu)航天器探測器
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    完全對換網(wǎng)絡(luò)的結(jié)構(gòu)連通度和子結(jié)構(gòu)連通度
    第二章 探測器有反應(yīng)
    EN菌的引力波探測器
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    第二章 探測器有反應(yīng)
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    鋼框架腹板雙角鋼連接梁柱子結(jié)構(gòu)抗倒塌性能分析
    基于子結(jié)構(gòu)的柴油機曲軸有限元建模方法研究
    在线天堂中文资源库| 成人三级做爰电影| 九色国产91popny在线| 中文字幕高清在线视频| 午夜福利高清视频| 琪琪午夜伦伦电影理论片6080| 成人午夜高清在线视频 | 丁香欧美五月| or卡值多少钱| 在线av久久热| 国产精品乱码一区二三区的特点| 一区福利在线观看| 成人午夜高清在线视频 | 国产成人系列免费观看| 国产亚洲精品久久久久5区| 在线观看日韩欧美| 窝窝影院91人妻| 动漫黄色视频在线观看| 亚洲成人国产一区在线观看| 免费看十八禁软件| 亚洲第一欧美日韩一区二区三区| 美国免费a级毛片| 999精品在线视频| 亚洲av美国av| 大香蕉久久成人网| 亚洲第一电影网av| 一级a爱视频在线免费观看| 成人三级黄色视频| 日日干狠狠操夜夜爽| 国产视频内射| 久久精品亚洲精品国产色婷小说| 国产野战对白在线观看| 国产精品av久久久久免费| 亚洲人成电影免费在线| 精品高清国产在线一区| 久久久久亚洲av毛片大全| 啦啦啦 在线观看视频| 男人的好看免费观看在线视频 | 午夜激情av网站| 亚洲国产精品sss在线观看| 淫秽高清视频在线观看| 国产区一区二久久| 不卡av一区二区三区| 午夜免费成人在线视频| 成人亚洲精品一区在线观看| av片东京热男人的天堂| 欧美 亚洲 国产 日韩一| 男人舔奶头视频| 国产精品日韩av在线免费观看| 欧美乱妇无乱码| 成人18禁高潮啪啪吃奶动态图| 精品福利观看| 免费电影在线观看免费观看| 好男人在线观看高清免费视频 | videosex国产| 三级毛片av免费| 国产激情欧美一区二区| 男人舔奶头视频| 成人亚洲精品一区在线观看| 国产亚洲欧美98| 在线国产一区二区在线| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 桃红色精品国产亚洲av| 黄色丝袜av网址大全| av免费在线观看网站| 日韩精品免费视频一区二区三区| 中国美女看黄片| 欧美国产日韩亚洲一区| 在线观看www视频免费| 日韩欧美一区二区三区在线观看| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 欧美一区二区精品小视频在线| 一二三四社区在线视频社区8| 嫁个100分男人电影在线观看| 日韩欧美国产在线观看| 国产99白浆流出| 丰满的人妻完整版| 成年人黄色毛片网站| 久久国产精品人妻蜜桃| 久久精品91蜜桃| 精品乱码久久久久久99久播| 别揉我奶头~嗯~啊~动态视频| 精品国内亚洲2022精品成人| 国产av在哪里看| 亚洲九九香蕉| 亚洲人成伊人成综合网2020| 超碰成人久久| 国产蜜桃级精品一区二区三区| 一级片免费观看大全| 欧美另类亚洲清纯唯美| 91大片在线观看| 激情在线观看视频在线高清| 88av欧美| 精品一区二区三区四区五区乱码| 女性被躁到高潮视频| 国产又爽黄色视频| 亚洲精品久久国产高清桃花| 免费观看精品视频网站| 婷婷精品国产亚洲av在线| 成在线人永久免费视频| 99久久综合精品五月天人人| 欧美成人午夜精品| 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 一边摸一边抽搐一进一小说| 啦啦啦观看免费观看视频高清| 老汉色∧v一级毛片| 亚洲久久久国产精品| 精品欧美国产一区二区三| 日韩大尺度精品在线看网址| 人人妻,人人澡人人爽秒播| 久久中文字幕人妻熟女| 欧美成人免费av一区二区三区| 757午夜福利合集在线观看| 日韩 欧美 亚洲 中文字幕| 99国产精品一区二区蜜桃av| 黄片播放在线免费| 自线自在国产av| 成人亚洲精品一区在线观看| 老司机福利观看| 激情在线观看视频在线高清| 99精品欧美一区二区三区四区| 日本一区二区免费在线视频| 日韩国内少妇激情av| 国产精品香港三级国产av潘金莲| 久久性视频一级片| 真人做人爱边吃奶动态| 日韩成人在线观看一区二区三区| 久久精品91无色码中文字幕| 日日爽夜夜爽网站| 欧美av亚洲av综合av国产av| 日本a在线网址| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久久久久久久 | 深夜精品福利| 男女视频在线观看网站免费 | 日韩大尺度精品在线看网址| ponron亚洲| 1024香蕉在线观看| 欧美zozozo另类| 中文字幕精品亚洲无线码一区 | 一进一出抽搐gif免费好疼| 国产三级在线视频| 国产亚洲欧美在线一区二区| 国产成人精品久久二区二区91| 性欧美人与动物交配| 国产熟女xx| 欧美av亚洲av综合av国产av| 国产高清有码在线观看视频 | 精品熟女少妇八av免费久了| 精品第一国产精品| 曰老女人黄片| 亚洲片人在线观看| 18美女黄网站色大片免费观看| 精品久久蜜臀av无| 免费在线观看视频国产中文字幕亚洲| 最好的美女福利视频网| 99国产综合亚洲精品| 午夜老司机福利片| 日日夜夜操网爽| 亚洲熟女毛片儿| 久久国产乱子伦精品免费另类| 久久久国产精品麻豆| 一边摸一边抽搐一进一小说| 免费看日本二区| 亚洲av成人一区二区三| 999久久久精品免费观看国产| 午夜福利一区二区在线看| 国内毛片毛片毛片毛片毛片| 最近在线观看免费完整版| 国产区一区二久久| 人妻丰满熟妇av一区二区三区| 亚洲最大成人中文| 日本免费一区二区三区高清不卡| 亚洲激情在线av| 国产成人精品久久二区二区91| 国产区一区二久久| 亚洲精品av麻豆狂野| 亚洲中文字幕一区二区三区有码在线看 | 18禁国产床啪视频网站| 俄罗斯特黄特色一大片| 热99re8久久精品国产| 黄片播放在线免费| 欧美日韩瑟瑟在线播放| 丰满的人妻完整版| 日韩一卡2卡3卡4卡2021年| 欧美日本视频| 1024香蕉在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲精品国产区一区二| 狠狠狠狠99中文字幕| 最近最新中文字幕大全电影3 | 一本久久中文字幕| 国产熟女xx| 国产高清激情床上av| 久久性视频一级片| 日本熟妇午夜| 天堂√8在线中文| 国产区一区二久久| 一二三四在线观看免费中文在| 一二三四社区在线视频社区8| 欧美黄色淫秽网站| 久热爱精品视频在线9| 免费高清视频大片| 亚洲avbb在线观看| 此物有八面人人有两片| 1024视频免费在线观看| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 免费在线观看黄色视频的| 久久久久久久久中文| 国产精品久久久av美女十八| 麻豆成人午夜福利视频| 岛国在线观看网站| 国产91精品成人一区二区三区| 在线观看日韩欧美| 啦啦啦观看免费观看视频高清| 国产成人欧美| 中文在线观看免费www的网站 | 丰满的人妻完整版| 欧美日韩中文字幕国产精品一区二区三区| 亚洲五月婷婷丁香| 国产激情偷乱视频一区二区| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| 亚洲国产中文字幕在线视频| 熟女电影av网| АⅤ资源中文在线天堂| 亚洲aⅴ乱码一区二区在线播放 | 嫁个100分男人电影在线观看| 99热这里只有精品一区 | 欧美黑人欧美精品刺激| 又大又爽又粗| 成人18禁在线播放| 黄色成人免费大全| 婷婷六月久久综合丁香| 1024视频免费在线观看| 正在播放国产对白刺激| 黑人巨大精品欧美一区二区mp4| 18禁国产床啪视频网站| 欧美乱妇无乱码| 怎么达到女性高潮| 人人妻人人澡欧美一区二区| 91国产中文字幕| 看片在线看免费视频| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区精品| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 熟女少妇亚洲综合色aaa.| 99久久久亚洲精品蜜臀av| 欧美激情高清一区二区三区| 成人欧美大片| 亚洲激情在线av| 欧美色欧美亚洲另类二区| 国产精品亚洲一级av第二区| 亚洲欧美一区二区三区黑人| 一区二区三区国产精品乱码| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 午夜影院日韩av| 国产色视频综合| 国产免费男女视频| 久久久国产精品麻豆| 亚洲国产精品久久男人天堂| 国产在线观看jvid| 亚洲成a人片在线一区二区| 夜夜看夜夜爽夜夜摸| 一夜夜www| 欧美最黄视频在线播放免费| 国产精品免费一区二区三区在线| 精品欧美一区二区三区在线| 亚洲美女黄片视频| 99热6这里只有精品| 黑人巨大精品欧美一区二区mp4| 久久中文字幕一级| 亚洲人成网站高清观看| 午夜免费观看网址| 国产精品 欧美亚洲| 国产三级在线视频| 99热6这里只有精品| 日本黄色视频三级网站网址| 久久天躁狠狠躁夜夜2o2o| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 国产人伦9x9x在线观看| 男女下面进入的视频免费午夜 | 一区二区三区激情视频| 久久久久久久精品吃奶| 国产熟女xx| 一二三四社区在线视频社区8| 色av中文字幕| 亚洲精品一区av在线观看| 欧美乱妇无乱码| 欧美国产精品va在线观看不卡| 国产黄a三级三级三级人| 久久伊人香网站| 亚洲精品一区av在线观看| 一本久久中文字幕| 国产精品98久久久久久宅男小说| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 欧美国产日韩亚洲一区| 欧美av亚洲av综合av国产av| 夜夜夜夜夜久久久久| 久久婷婷成人综合色麻豆| 亚洲五月天丁香| 天天一区二区日本电影三级| 国产伦人伦偷精品视频| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美网| 99在线人妻在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 国产又色又爽无遮挡免费看| 精品欧美一区二区三区在线| 精品久久久久久,| 可以在线观看的亚洲视频| 在线观看免费视频日本深夜| 国产精品 国内视频| 免费在线观看视频国产中文字幕亚洲| 免费看美女性在线毛片视频| 亚洲av成人一区二区三| 久久久久久久午夜电影| 亚洲精品av麻豆狂野| 国产私拍福利视频在线观看| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 首页视频小说图片口味搜索| 好男人电影高清在线观看| 久久九九热精品免费| 曰老女人黄片| 99久久无色码亚洲精品果冻| 日韩精品中文字幕看吧| 午夜亚洲福利在线播放| 久久精品亚洲精品国产色婷小说| 久久人妻福利社区极品人妻图片| 热re99久久国产66热| 十八禁网站免费在线| 亚洲天堂国产精品一区在线| 亚洲成人国产一区在线观看| 久久国产精品人妻蜜桃| 亚洲av美国av| 99精品在免费线老司机午夜| 亚洲熟妇熟女久久| 国产精品野战在线观看| 黄色视频不卡| 亚洲欧美日韩高清在线视频| 一二三四社区在线视频社区8| 91大片在线观看| 国产精品久久久人人做人人爽| 色av中文字幕| 搡老岳熟女国产| 欧美黄色淫秽网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产成人精品二区| 国产精品影院久久| 亚洲成a人片在线一区二区| 午夜免费成人在线视频| 国产区一区二久久| 亚洲国产精品sss在线观看| 国产亚洲欧美98| 日本五十路高清| 无人区码免费观看不卡| 午夜福利在线在线| 国产视频内射| 国产亚洲精品第一综合不卡| 女人被狂操c到高潮| 在线视频色国产色| 久久精品国产99精品国产亚洲性色| 国产亚洲精品第一综合不卡| 少妇的丰满在线观看| 国产野战对白在线观看| 久久草成人影院| 无遮挡黄片免费观看| 成年免费大片在线观看| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 他把我摸到了高潮在线观看| 久久热在线av| 中文字幕人妻丝袜一区二区| 50天的宝宝边吃奶边哭怎么回事| 国产高清有码在线观看视频 | 国产99白浆流出| 两性夫妻黄色片| 亚洲电影在线观看av| 99久久无色码亚洲精品果冻| 少妇 在线观看| 国产黄色小视频在线观看| 国产伦一二天堂av在线观看| 满18在线观看网站| av免费在线观看网站| 又大又爽又粗| 999精品在线视频| 久9热在线精品视频| 久久婷婷成人综合色麻豆| 真人一进一出gif抽搐免费| 国产一级毛片七仙女欲春2 | 巨乳人妻的诱惑在线观看| 精品午夜福利视频在线观看一区| 午夜福利欧美成人| 国产又黄又爽又无遮挡在线| 99精品欧美一区二区三区四区| 50天的宝宝边吃奶边哭怎么回事| 国产三级黄色录像| 亚洲国产日韩欧美精品在线观看 | 久久香蕉精品热| 黄片大片在线免费观看| 免费在线观看影片大全网站| 午夜激情福利司机影院| 日韩成人在线观看一区二区三区| 亚洲精品国产一区二区精华液| 日韩欧美一区视频在线观看| 国产高清videossex| 91麻豆精品激情在线观看国产| 午夜免费鲁丝| 午夜福利免费观看在线| 亚洲国产欧美网| 欧美日韩亚洲综合一区二区三区_| 国产成+人综合+亚洲专区| 亚洲精品中文字幕一二三四区| 在线十欧美十亚洲十日本专区| 搡老妇女老女人老熟妇| 激情在线观看视频在线高清| 午夜免费激情av| 国产99白浆流出| 中文字幕人妻熟女乱码| 神马国产精品三级电影在线观看 | 久久精品亚洲精品国产色婷小说| cao死你这个sao货| 亚洲熟女毛片儿| 真人一进一出gif抽搐免费| 国产爱豆传媒在线观看 | 亚洲午夜理论影院| 51午夜福利影视在线观看| 别揉我奶头~嗯~啊~动态视频| 精品国产亚洲在线| 国产国语露脸激情在线看| 国产av在哪里看| 精品久久蜜臀av无| 国产精品久久久人人做人人爽| 欧美久久黑人一区二区| 久久久久久国产a免费观看| 9191精品国产免费久久| а√天堂www在线а√下载| 激情在线观看视频在线高清| 又大又爽又粗| 免费在线观看影片大全网站| 国产日本99.免费观看| 亚洲精品一区av在线观看| 精品午夜福利视频在线观看一区| 亚洲激情在线av| 中文字幕久久专区| 制服丝袜大香蕉在线| 日韩精品中文字幕看吧| 一区二区三区精品91| 国产高清视频在线播放一区| 久久久久久亚洲精品国产蜜桃av| 欧美黑人精品巨大| 欧美一级a爱片免费观看看 | 国产成人欧美在线观看| 久久精品aⅴ一区二区三区四区| 91大片在线观看| 欧美激情极品国产一区二区三区| 欧美性长视频在线观看| 99riav亚洲国产免费| 丝袜在线中文字幕| 国产区一区二久久| 男女做爰动态图高潮gif福利片| 亚洲在线自拍视频| 亚洲无线在线观看| 女性被躁到高潮视频| www.www免费av| 久久久久久人人人人人| 欧美丝袜亚洲另类 | 久久精品人妻少妇| 国产视频内射| avwww免费| 精品国产乱子伦一区二区三区| 日韩视频一区二区在线观看| 国产精品一区二区精品视频观看| 午夜精品久久久久久毛片777| 波多野结衣av一区二区av| 欧美国产日韩亚洲一区| 国产亚洲欧美在线一区二区| 婷婷亚洲欧美| 别揉我奶头~嗯~啊~动态视频| 欧美av亚洲av综合av国产av| 亚洲人成网站高清观看| 国内毛片毛片毛片毛片毛片| 亚洲人成网站在线播放欧美日韩| 老司机午夜十八禁免费视频| 在线永久观看黄色视频| 人人妻人人看人人澡| 久久中文字幕一级| 亚洲中文字幕日韩| 非洲黑人性xxxx精品又粗又长| 99re在线观看精品视频| 看免费av毛片| av视频在线观看入口| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久国产a免费观看| 日本熟妇午夜| 欧美色视频一区免费| 成人三级黄色视频| 欧美一级a爱片免费观看看 | 国产激情久久老熟女| 亚洲熟妇中文字幕五十中出| 搞女人的毛片| 真人一进一出gif抽搐免费| 在线视频色国产色| 欧美一级a爱片免费观看看 | 亚洲精品国产一区二区精华液| 亚洲av成人av| 日韩国内少妇激情av| 成年人黄色毛片网站| 精品久久久久久,| 国产精品一区二区免费欧美| 国产伦人伦偷精品视频| 999精品在线视频| 老司机午夜福利在线观看视频| 欧美黑人欧美精品刺激| 欧美性长视频在线观看| 欧美人与性动交α欧美精品济南到| 欧美av亚洲av综合av国产av| 我的亚洲天堂| 亚洲av第一区精品v没综合| 亚洲中文字幕一区二区三区有码在线看 | 丝袜美腿诱惑在线| 在线观看免费视频日本深夜| 亚洲成a人片在线一区二区| 亚洲成人国产一区在线观看| 国产伦在线观看视频一区| 在线观看免费日韩欧美大片| 亚洲欧美日韩无卡精品| 91老司机精品| 两个人视频免费观看高清| 免费在线观看影片大全网站| 校园春色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 黄色成人免费大全| 日本撒尿小便嘘嘘汇集6| 国产精品,欧美在线| 精品卡一卡二卡四卡免费| 日韩高清综合在线| 村上凉子中文字幕在线| 高潮久久久久久久久久久不卡| 亚洲精品中文字幕一二三四区| 成人av一区二区三区在线看| 99热这里只有精品一区 | 国产野战对白在线观看| 国产成人av激情在线播放| 精品免费久久久久久久清纯| 一进一出抽搐gif免费好疼| www.精华液| 亚洲无线在线观看| 午夜两性在线视频| 老司机午夜福利在线观看视频| 亚洲精品美女久久av网站| 欧美激情久久久久久爽电影| 三级毛片av免费| 久久99热这里只有精品18| 亚洲av片天天在线观看| 国产成人影院久久av| 亚洲国产欧美日韩在线播放| 免费看日本二区| 18禁国产床啪视频网站| 亚洲欧美日韩高清在线视频| 91成年电影在线观看| 久久国产精品人妻蜜桃| 19禁男女啪啪无遮挡网站| 少妇熟女aⅴ在线视频| 国产极品粉嫩免费观看在线| 国产欧美日韩一区二区三| 欧美日本视频| 亚洲七黄色美女视频| 午夜影院日韩av| 欧美zozozo另类| 国产一级毛片七仙女欲春2 | 亚洲真实伦在线观看| 国产区一区二久久| 久久精品影院6| 国产精品亚洲美女久久久| 国产成人精品无人区| 淫妇啪啪啪对白视频| 亚洲中文av在线| 亚洲成人免费电影在线观看| 亚洲欧美精品综合久久99| 少妇被粗大的猛进出69影院| 国产成人av激情在线播放| 国产一区二区激情短视频| 国产成人av教育| 午夜a级毛片| 国产熟女xx| 麻豆成人午夜福利视频| 亚洲精品在线美女| 欧美中文综合在线视频| 亚洲人成77777在线视频| 亚洲国产日韩欧美精品在线观看 | 又黄又爽又免费观看的视频| 免费看十八禁软件| 一级a爱片免费观看的视频| 亚洲中文日韩欧美视频| 亚洲九九香蕉| 国产爱豆传媒在线观看 | 亚洲 欧美一区二区三区| 母亲3免费完整高清在线观看| 又黄又爽又免费观看的视频| 日本撒尿小便嘘嘘汇集6| 久久久水蜜桃国产精品网|