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

    大型航天器桁架式主承力結(jié)構(gòu)構(gòu)型拓?fù)鋬?yōu)化研究

    2014-12-28 05:45:52
    航天器工程 2014年2期
    關(guān)鍵詞:力點(diǎn)貯箱桿件

    (中國(guó)空間技術(shù)研究院,北京 100094)

    1 引言

    以桁架作為主承力結(jié)構(gòu)的大型航天器,如大型桁架式衛(wèi)星等,其構(gòu)型設(shè)計(jì)所受的約束條件頗多,運(yùn)載火箭、推進(jìn)方案、熱控方案和有效載荷類型等對(duì)構(gòu)型設(shè)計(jì)均有不同程度的影響。此外,還要考慮傳力、平臺(tái)設(shè)備之間的電連接、傳熱,以及裝配和維護(hù)等要求。在構(gòu)型設(shè)計(jì)初期,首先是滿足傳力要求。合理的傳力路徑能用較少的結(jié)構(gòu)材料獲得較高的結(jié)構(gòu)剛度和強(qiáng)度。

    隨著航天技術(shù)的發(fā)展,滿足現(xiàn)代航天器日趨復(fù)雜的結(jié)構(gòu)本身對(duì)結(jié)構(gòu)性能越來越高的要求,必須應(yīng)用結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[1]。國(guó)外某些航空航天公司已將結(jié)構(gòu)優(yōu)化作為結(jié)構(gòu)設(shè)計(jì)中的規(guī)范環(huán)節(jié)[2];在國(guó)內(nèi),優(yōu)化設(shè)計(jì)方法也開始在航天器結(jié)構(gòu)設(shè)計(jì)中得到初步應(yīng)用[3-4]?,F(xiàn)階段,結(jié)構(gòu)參數(shù)優(yōu)化的應(yīng)用研究較多,而結(jié)構(gòu)拓?fù)鋬?yōu)化的應(yīng)用研究相對(duì)較少。

    本文以大型航天器平臺(tái)桁架式承力結(jié)構(gòu)的研發(fā)為背景,在其概念設(shè)計(jì)階段,根據(jù)平臺(tái)設(shè)計(jì)要求和外部約束,應(yīng)用連續(xù)體-離散體兩級(jí)拓?fù)鋬?yōu)化策略得到了桁架式承力結(jié)構(gòu)的構(gòu)型形式,所得結(jié)果可為詳細(xì)設(shè)計(jì)階段提供依據(jù)。與尺寸優(yōu)化和截面形狀優(yōu)化相比,桁架結(jié)構(gòu)構(gòu)型拓?fù)鋬?yōu)化能夠在更大范圍內(nèi)改善結(jié)構(gòu)特性,進(jìn)而使所得承力結(jié)構(gòu)質(zhì)量更小,性能更優(yōu)。本文提出的優(yōu)化策略也可用于其他桁架式承力結(jié)構(gòu)的設(shè)計(jì)工作。

    2 連續(xù)體-離散體兩級(jí)拓?fù)鋬?yōu)化策略

    桁架式結(jié)構(gòu)強(qiáng)度高,通用性、組合性好,形式簡(jiǎn)單,傳力路線合理,易于實(shí)現(xiàn)模塊化且可擴(kuò)展,因此特別適于用作某些大型航天器的主承力結(jié)構(gòu)。桁架結(jié)構(gòu)能夠傳遞集中載荷,沿著載荷的主要傳力路徑布置桿件,能夠顯著減小結(jié)構(gòu)質(zhì)量,提高結(jié)構(gòu)性能。

    結(jié)構(gòu)拓?fù)鋬?yōu)化是一種尋找結(jié)構(gòu)傳力路徑的有效方法。根據(jù)優(yōu)化對(duì)象的結(jié)構(gòu)類型的不同,結(jié)構(gòu)拓?fù)鋬?yōu)化可分為連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化和離散體結(jié)構(gòu)拓?fù)鋬?yōu)化,在此前的研究工作中,兩者一般不同時(shí)應(yīng)用。前者用于從二維或三維設(shè)計(jì)區(qū)域中去除低效能材料,得到的剩余材料的分布情況可指導(dǎo)詳細(xì)設(shè)計(jì)[5];后者多采用基結(jié)構(gòu)法,優(yōu)化過程從大量初始桿件(即基結(jié)構(gòu))中刪除低效能桿件,得到所需的桁架結(jié)構(gòu),但其最初基結(jié)構(gòu)的主要節(jié)點(diǎn)和桿件連接關(guān)系則是人為給定的。

    為將連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化和離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)合起來進(jìn)行桁架結(jié)構(gòu)的設(shè)計(jì),提出連續(xù)體-離散體兩級(jí)拓?fù)鋬?yōu)化策略(見圖1)。首先,進(jìn)行連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化。根據(jù)桁架結(jié)構(gòu)總體設(shè)計(jì)要求,給出連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化的可行設(shè)計(jì)區(qū)域,并設(shè)定目標(biāo)函數(shù)和約束條件,形成優(yōu)化問題。利用工具軟件(如HyperWorks/OptiStruct)求解計(jì)算可得到類似桁架結(jié)構(gòu)的優(yōu)化結(jié)果,在此基礎(chǔ)上進(jìn)行歸納提煉,得到由桿件組成的標(biāo)準(zhǔn)桁架結(jié)構(gòu)形式。然后,進(jìn)行離散體結(jié)構(gòu)拓?fù)鋬?yōu)化。以連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化得到的桁架結(jié)構(gòu)為基礎(chǔ),根據(jù)經(jīng)驗(yàn)人為地添加節(jié)點(diǎn)和桿件,增加其他可能的傳力路徑,得到其拓展基結(jié)構(gòu),再按照離散體結(jié)構(gòu)拓?fù)鋬?yōu)化數(shù)學(xué)模型進(jìn)行優(yōu)化,即可得到最終的桁架式結(jié)構(gòu)構(gòu)型。如有多個(gè)方案,還可對(duì)比選優(yōu)。

    圖1 桁架結(jié)構(gòu)優(yōu)化設(shè)計(jì)流程Fig.1 Flow chart of optimal design for truss structure

    3 大型航天器桁架式承力結(jié)構(gòu)拓?fù)鋬?yōu)化

    3.1 問題描述

    設(shè)定大型航天器采用模塊化設(shè)計(jì),其主體結(jié)構(gòu)由載荷艙和平臺(tái)兩部分組成,其中平臺(tái)部分采用桁架式的主承力結(jié)構(gòu),用于連接平臺(tái)底部航天器和火箭的對(duì)接環(huán)和頂部載荷艙,見圖2。以對(duì)接環(huán)底面圓心為坐標(biāo)原點(diǎn),各坐標(biāo)軸方向如圖2所示。

    圖2 航天器主體結(jié)構(gòu)示意圖Fig.2 Sketch map for main body of spacecraft

    航天器的總質(zhì)量為7000kg時(shí),其剛度滿足一階橫向頻率不小于8 Hz、一階縱向頻率不小于30Hz的要求;航天器平臺(tái)應(yīng)滿足所采用運(yùn)載火箭整流罩的包絡(luò)要求;載荷艙要求平臺(tái)頂端應(yīng)有不少于5個(gè)的傳力點(diǎn);在保持航天器結(jié)構(gòu)構(gòu)型不變的前提下,結(jié)構(gòu)平臺(tái)的承載能力可擴(kuò)展;航天器和火箭的對(duì)接環(huán)直徑不小于2000 mm。在上述要求的約束下,對(duì)航天器平臺(tái)部分桁架式承力結(jié)構(gòu)的構(gòu)型形式進(jìn)行設(shè)計(jì)。

    3.2 連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化

    3.2.1 邊界條件和載荷條件

    受運(yùn)載火箭包絡(luò)尺寸的限制,并考慮平臺(tái)安裝太陽(yáng)翼、天線等設(shè)備和散熱的要求,平臺(tái)基本尺寸確定為2880mm×2880mm。推進(jìn)劑貯箱直接安裝在對(duì)接環(huán)上,四貯箱并聯(lián)安裝。根據(jù)推進(jìn)劑預(yù)算可確定單個(gè)貯箱容積,設(shè)計(jì)為直徑大于1000 mm、高度大于1300mm 的球頭貯箱。對(duì)接環(huán)直徑按要求暫取為2000mm。

    根據(jù)運(yùn)載方規(guī)定的航天器設(shè)計(jì)載荷條件,并考慮1.25的安全系數(shù),確定平臺(tái)結(jié)構(gòu)X和Y兩個(gè)橫向加速度載荷為3gn,縱向(Z向)加速度載荷為8gn。

    需要說明的是,除幾何包絡(luò)要求、邊界約束條件、載荷加載形式和初步確定的承力點(diǎn)位置外,連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果對(duì)設(shè)計(jì)區(qū)域結(jié)構(gòu)的材料特性、集中質(zhì)量點(diǎn)質(zhì)量的絕對(duì)大小以及載荷的絕對(duì)數(shù)值都是不敏感的。不過,各方向載荷的相對(duì)大小變化時(shí),其合矢量方向也會(huì)變化,對(duì)優(yōu)化結(jié)果會(huì)產(chǎn)生一定的影響。

    3.2.2 優(yōu)化方法及優(yōu)化問題數(shù)學(xué)模型

    連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化應(yīng)用HyperWorks軟件,采用變密度法進(jìn)行優(yōu)化[6],應(yīng)用的是固體各向同性懲罰微結(jié)構(gòu)模型(Solid Isotropic Microstructures with Penalization,SIMP)[7]。該材料模型將結(jié)構(gòu)設(shè)計(jì)域離散為n個(gè)有限單元,每個(gè)單元的彈性模量與其相對(duì)密度的關(guān)系可表示為

    式中:xi為單元i的相對(duì)密度;E0為實(shí)體材料的彈性模量;Ei為單元i材料的彈性模量;p表示懲罰因子,在OptiStruct軟件中,與p有關(guān)的優(yōu)化參數(shù)DISCRETE取值為(p-1)[8]。

    圖3給出了p取不同值時(shí)相對(duì)彈性模量(Ei/E0)與相對(duì)密度xi的關(guān)系曲線。

    圖3 相對(duì)彈性模量與相對(duì)密度的關(guān)系Fig.3 Relative stiffness vs relative density

    假定彈性體在受力作用的過程中始終保持平衡,因而沒有動(dòng)能的改變,且彈性體的非機(jī)械能也沒有變化,于是,外力勢(shì)能的減少(也就是外力所做的功)就完全轉(zhuǎn)變?yōu)樾巫儎?shì)能(又稱為應(yīng)變能)[9]。在這種情況下,外力功和應(yīng)變能在數(shù)值上是一致的,可稱該值為靜態(tài)應(yīng)變能,作為反映結(jié)構(gòu)靜剛度的一個(gè)標(biāo)量響應(yīng)。那么,優(yōu)化問題可描述為:受加速度載荷時(shí),設(shè)計(jì)區(qū)域體積分?jǐn)?shù)和航天器基頻約束下的結(jié)構(gòu)靜態(tài)應(yīng)變能最小,其設(shè)計(jì)變量為結(jié)構(gòu)有限元模型中設(shè)計(jì)區(qū)域內(nèi)每個(gè)單元的相對(duì)密度值。連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化問題的數(shù)學(xué)模型如下。

    在該模型中,X=(x1,x2,…,xn),為拓?fù)鋬?yōu)化設(shè)計(jì)變量向量,xi為第i個(gè)單元的相對(duì)密度;C(X)為反映結(jié)構(gòu)靜態(tài)應(yīng)變能的量;Ui為單元i的位移向量;Ki為單元i的剛度矩陣;V(X)為設(shè)計(jì)區(qū)域結(jié)構(gòu)總體積;fV為設(shè)計(jì)區(qū)域體積分?jǐn)?shù),vi為單元i的體積;fmin和分別為結(jié)構(gòu)基頻及其下限;δ為一極小正數(shù),用于避免剛度矩陣奇異。

    3.2.3 結(jié)構(gòu)有限元建模

    根據(jù)已確定的邊界條件,平臺(tái)結(jié)構(gòu)占據(jù)的幾何空間為一立方體,XY平面尺寸已知,高度則由貯箱高度確定。因?yàn)橥七M(jìn)劑貯箱不能與其他結(jié)構(gòu)發(fā)生干涉,所以立方體中還要預(yù)留出貯箱空間。氧化劑貯箱與燃料貯箱分別對(duì)角放置,使用梁?jiǎn)卧M貯箱結(jié)構(gòu),單元自身的結(jié)構(gòu)質(zhì)量用來模擬未填充推進(jìn)劑貯箱的結(jié)構(gòu)質(zhì)量,推進(jìn)劑質(zhì)量則通過非結(jié)構(gòu)質(zhì)量進(jìn)行模擬。貯箱下端直接連接到對(duì)接環(huán)上,其上端則通過拉桿與平臺(tái)結(jié)構(gòu)相連。平臺(tái)頂端的傳力點(diǎn)也是平臺(tái)結(jié)構(gòu)的主要承載點(diǎn),因此貯箱拉桿也是連接到這些傳力點(diǎn)上。有限元模型中使用集中質(zhì)量點(diǎn)來模擬載荷艙的質(zhì)量,考慮到設(shè)計(jì)要求中規(guī)定平臺(tái)頂端應(yīng)有不少于5個(gè)的傳力點(diǎn),提出五傳力點(diǎn)和七傳力點(diǎn)兩種質(zhì)量分布方案。

    基于以上考慮,本文構(gòu)建了如圖4所示的有限元模型。其中:五傳力點(diǎn)模型共有43 193個(gè)節(jié)點(diǎn),36 430個(gè)單元;七傳力點(diǎn)模型共有43 193個(gè)節(jié)點(diǎn),36 434個(gè)單元。

    圖4 設(shè)計(jì)區(qū)域的結(jié)構(gòu)有限元模型Fig.4 Finite element model for design domain

    3.2.4 優(yōu)化結(jié)果

    基于已建立的有限元模型和優(yōu)化的數(shù)學(xué)模型,通過不斷調(diào)整優(yōu)化數(shù)學(xué)模型中的體積分?jǐn)?shù)和結(jié)構(gòu)基頻約束進(jìn)行多次試算,并在已有優(yōu)化結(jié)果上進(jìn)行二次優(yōu)化,最終得到五傳力點(diǎn)方案和七傳力點(diǎn)方案的連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果。兩個(gè)方案的優(yōu)化結(jié)果傳力路徑清晰,易于提取出對(duì)應(yīng)的桁架形式,可作為后續(xù)離散體結(jié)構(gòu)拓?fù)鋬?yōu)化的輸入條件。連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果及對(duì)應(yīng)桁架形式如圖5和圖6所示。

    圖5 五傳力點(diǎn)方案連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果及對(duì)應(yīng)的桁架形式Fig.5 Optimization result of continuum structure for fivejuncture scheme and its corresponding truss configuration

    圖6 七傳力點(diǎn)方案連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果及其對(duì)應(yīng)的桁架形式Fig.6 Optimization result of continuum structure for sevenjuncture scheme and its corresponding truss configuration

    3.3 離散體結(jié)構(gòu)拓?fù)鋬?yōu)化

    通過連續(xù)體拓?fù)鋬?yōu)化,已經(jīng)得到了桁架結(jié)構(gòu)主要的節(jié)點(diǎn)信息和桿件分布形式,這時(shí)可對(duì)桁架結(jié)構(gòu)進(jìn)行進(jìn)一步優(yōu)化,增加必要的節(jié)點(diǎn)和桿件,再通過離散體結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)使桁架結(jié)構(gòu)的構(gòu)型更加合理。

    3.3.1 邊界條件和載荷條件

    離散體結(jié)構(gòu)拓?fù)鋬?yōu)化時(shí)的邊界條件和載荷條件,與連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化時(shí)基本相同,只是設(shè)計(jì)空間由三維實(shí)體結(jié)構(gòu)換成了空間桁架。

    3.3.2 優(yōu)化方法及優(yōu)化問題數(shù)學(xué)模型

    求解離散體結(jié)構(gòu)拓?fù)鋬?yōu)化問題時(shí)常用的是基結(jié)構(gòu)法。在設(shè)計(jì)空間內(nèi)人為地布置足夠多的節(jié)點(diǎn),并將每一個(gè)節(jié)點(diǎn)與其他所有節(jié)點(diǎn)用桿件連接起來形成基結(jié)構(gòu)。在此基礎(chǔ)上進(jìn)行截面尺寸優(yōu)化,并將截面過小的桿件刪除,從而實(shí)現(xiàn)結(jié)構(gòu)的拓?fù)鋬?yōu)化。這就是基結(jié)構(gòu)法的基本思想。

    離散體結(jié)構(gòu)拓?fù)鋬?yōu)化的設(shè)計(jì)變量為設(shè)計(jì)桿元的橫截面面積,非設(shè)計(jì)桿元在優(yōu)化過程中截面尺寸保持不變。優(yōu)化目標(biāo)為過載工況下的結(jié)構(gòu)靜態(tài)應(yīng)變能最?。患s束條件為優(yōu)化后桿件總體積不超過給定值(一般設(shè)置為初始桿件總體積的1/3),結(jié)構(gòu)一階頻率大于某給定值。離散體結(jié)構(gòu)拓?fù)鋬?yōu)化的數(shù)學(xué)模型如下。

    在該模型中,X′=(A1,A2,…,Am),為拓?fù)鋬?yōu)化設(shè)計(jì)中各桿件的橫截面面積,Ai為第i根桿的橫截面面積,m為桿件數(shù);C(X′)為結(jié)構(gòu)靜態(tài)應(yīng)變能;F為節(jié)點(diǎn)等效載荷向量;U為節(jié)點(diǎn)位移向量;V(X′)為結(jié)構(gòu)體積;li為第i根桿的長(zhǎng)度;ˉV為體積約束上限;K為結(jié)構(gòu)總剛度矩陣。

    該模型所描述的優(yōu)化問題中存在平衡方程約束,因此是一個(gè)非凸優(yōu)化問題,難以找到全局最優(yōu)解。為便于求解,可將優(yōu)化模型進(jìn)行轉(zhuǎn)化。應(yīng)用半定規(guī)劃理論,通過一定的數(shù)學(xué)處理將該非凸模型轉(zhuǎn)化為凸性的半定規(guī)劃模型,即可通過內(nèi)點(diǎn)法進(jìn)行求解[10]。需要說明的是,模型中的平衡方程F=KU必然滿足,但仍將其列出;靜態(tài)應(yīng)變能在此處選用外力功的等價(jià)表達(dá)式,這樣的處理使得優(yōu)化模型可以轉(zhuǎn)化為半定規(guī)劃的標(biāo)準(zhǔn)形式,便于求解。

    3.3.3 基結(jié)構(gòu)的拓展及有限元建模

    在連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化得到的桁架結(jié)構(gòu)基礎(chǔ)上進(jìn)行拓展,增加必要的節(jié)點(diǎn)和桿件,添加可能的傳力路徑,形成桁架的拓展基結(jié)構(gòu),如圖7所示。使用桁架式結(jié)構(gòu)拓?fù)湓O(shè)計(jì)優(yōu)化系統(tǒng)(TOpology Design Optimization System of Structure with Truss,TODOSST)[11]進(jìn)行離散體結(jié)構(gòu)拓?fù)鋬?yōu)化。與連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化相比,離散體結(jié)構(gòu)拓?fù)鋬?yōu)化有限元模型以桿結(jié)構(gòu)為主,為各桿件賦予材料和尺寸信息;為貯箱設(shè)置適當(dāng)?shù)姆墙Y(jié)構(gòu)質(zhì)量來模擬推進(jìn)劑質(zhì)量;對(duì)結(jié)構(gòu)相應(yīng)節(jié)點(diǎn)施加固定支撐約束,并設(shè)置過載工況。各傳力點(diǎn)的集中質(zhì)量可在TODOSST 的人機(jī)交互界面中添加。在界面中按照優(yōu)化模型完成設(shè)置后,即可在MATLAB軟件環(huán)境中運(yùn)行TODOSST,進(jìn)行離散體結(jié)構(gòu)拓?fù)鋬?yōu)化。

    圖7 五傳力點(diǎn)和七傳力點(diǎn)方案桁架的拓展基結(jié)構(gòu)Fig.7 Expanded base structures for both five-juncture and seven-juncture schemes

    3.3.4 優(yōu)化結(jié)果

    以過載工況下的結(jié)構(gòu)靜態(tài)應(yīng)變能最小為優(yōu)化目標(biāo),優(yōu)化后桿件體積約束設(shè)置為結(jié)構(gòu)初始總體積的1/3,并分別取不同的一階頻率約束,所得的優(yōu)化結(jié)果分別如圖8和圖9所示。圖中均去除了截面積小于結(jié)構(gòu)中最大桿件截面積1%的桿件;灰色部分為非設(shè)計(jì)桿元,黑色部分為設(shè)計(jì)桿元,桿件顯示的粗細(xì)僅表示其橫截面積的相對(duì)大小,不代表實(shí)際數(shù)值。

    從圖8和圖9可見,對(duì)五傳力點(diǎn)方案,當(dāng)一階頻率約束從12Hz開始不斷提高至22 Hz時(shí),所得的離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果的拓?fù)湫问绞峭耆恢碌?,且?dāng)外圍桿件截面積增大時(shí),結(jié)構(gòu)基頻也相應(yīng)提高。一階頻率約束大于24 Hz時(shí),未能得到合理的優(yōu)化結(jié)果。對(duì)于七傳力點(diǎn)方案,當(dāng)一階頻率約束從10Hz開始不斷提高至24 Hz時(shí),所得的離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果的拓?fù)湫问绞峭耆恢碌?。同樣的,?dāng)外圍桿件截面積增大時(shí),結(jié)構(gòu)基頻也相應(yīng)提高,一階頻率約束大于27Hz時(shí),未能得到合理的優(yōu)化結(jié)果。

    在以上優(yōu)化結(jié)果中,刪除了某些次元(截面尺寸相對(duì)較小的桿件)之后,有可能產(chǎn)生機(jī)構(gòu),此時(shí)要人為合并相應(yīng)節(jié)點(diǎn)以形成結(jié)構(gòu)。圖10給出了以上離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果所得出的最優(yōu)拓?fù)湫问剑渲衅邆髁c(diǎn)方案中暫時(shí)保留了外側(cè)小橫桿和“X”形交叉桿,可以認(rèn)為它們模擬了外側(cè)箱板的支撐作用,不屬于桁架結(jié)構(gòu)的組成部分。

    圖8 五傳力點(diǎn)方案中不同一階頻率約束時(shí)的離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果Fig.8 Topology optimization results of discrete structure for five-juncture scheme under various first frequency constraints

    圖9 七傳力點(diǎn)方案中不同一階頻率約束時(shí)的離散體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果Fig.9 Topology optimization results of discrete structure for seven-juncture scheme under various first frequency constraints

    圖10 五傳力點(diǎn)和七傳力點(diǎn)方案經(jīng)優(yōu)化后的最終拓?fù)湫问紽ig.10 Final topology after optimization for both five-juncture and seven-juncture schemes

    3.4 優(yōu)化結(jié)構(gòu)分析比較

    兩種方案的優(yōu)化設(shè)計(jì)結(jié)果均滿足設(shè)計(jì)要求中對(duì)于運(yùn)載火箭包絡(luò)尺寸、對(duì)接環(huán)直徑和與載荷艙接口數(shù)量的要求。圖8和圖9的優(yōu)化結(jié)果顯示,經(jīng)過后期的詳細(xì)設(shè)計(jì),所得的桁架構(gòu)型形式能滿足航天器基頻要求,并具有一定的拓展能力。

    對(duì)比圖8和圖9可發(fā)現(xiàn),在不改變桁架結(jié)構(gòu)構(gòu)型的前提下,僅通過修改模型中各桿件的橫截面尺寸,七傳力點(diǎn)方案在給定的條件下能夠達(dá)到比五傳力點(diǎn)方案更高的結(jié)構(gòu)基頻,因此具有更大的設(shè)計(jì)和優(yōu)化空間。在最初確定貯箱分布方式時(shí),有意將貯箱沿X方向的間隔拉大,以在中間位置預(yù)留更多的設(shè)備安裝空間。從五傳力點(diǎn)方案的優(yōu)化結(jié)果可以看到,預(yù)留空間部分已被承力結(jié)構(gòu)占據(jù),該處可用空間被承力桿件分割得較為零碎,不利于設(shè)備的安裝;若考慮加裝隔板來增加結(jié)構(gòu)的剛度和承載能力,由桿件的位置可以看出,隔板分布將會(huì)非常密集,且有多塊隔板為斜板(不與XOY平面垂直),即使結(jié)構(gòu)剛度和承載能力有所提高,也會(huì)增加較多的結(jié)構(gòu)質(zhì)量,浪費(fèi)更多的可用空間。另外,七傳力點(diǎn)方案中每個(gè)貯箱頂端的3個(gè)拉桿位置分布較均勻,拉桿長(zhǎng)度也較為合理;五傳力點(diǎn)方案中拉桿分布相對(duì)來說不是很均勻,其中一根拉桿的長(zhǎng)度明顯比七傳力點(diǎn)方案中的長(zhǎng),這樣既增加結(jié)構(gòu)質(zhì)量,又使傳力路徑變長(zhǎng),因而不是較為合理的設(shè)計(jì)方案。

    根據(jù)以上幾點(diǎn)分析,七傳力點(diǎn)方案的優(yōu)化結(jié)果具有比較明顯的優(yōu)勢(shì),建議以七傳力點(diǎn)方案的最優(yōu)拓?fù)湫问綖榛A(chǔ)開展后續(xù)設(shè)計(jì)工作。

    4 結(jié)構(gòu)參數(shù)優(yōu)化的討論

    在結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)階段,優(yōu)化對(duì)象僅是桁架式主承力結(jié)構(gòu)的構(gòu)型形式,即主要節(jié)點(diǎn)的位置和節(jié)點(diǎn)間的桿件連接情況,后續(xù)工作中必須進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計(jì),其必要性主要體現(xiàn)在:

    (1)從能量角度講,結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)時(shí)的載荷艙結(jié)構(gòu)僅由幾個(gè)集中質(zhì)量點(diǎn)來模擬,未考慮實(shí)際載荷艙的彈性變形對(duì)基礎(chǔ)振動(dòng)能量的吸收、儲(chǔ)存作用,也就是說,平臺(tái)桁架結(jié)構(gòu)要自行吸收全部振動(dòng)能量,設(shè)計(jì)偏保守。

    (2)在實(shí)際發(fā)射過程中,一般不會(huì)出現(xiàn)3個(gè)方向同時(shí)受最大載荷的過載工況,進(jìn)行結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)時(shí)應(yīng)用這種嚴(yán)酷的載荷條件也會(huì)使設(shè)計(jì)偏保守。

    (3)結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果對(duì)結(jié)構(gòu)的材料特性、加速度載荷的絕對(duì)數(shù)值等都是不敏感的,設(shè)計(jì)要求中的某些技術(shù)指標(biāo)在該階段起不到具體的約束作用。只有在航天器結(jié)構(gòu)進(jìn)行較細(xì)致建模時(shí),航天器模型的力學(xué)特性才能達(dá)到與技術(shù)指標(biāo)進(jìn)行對(duì)比的程度,有必要進(jìn)一步進(jìn)行優(yōu)化。

    (4)經(jīng)過初步的詳細(xì)設(shè)計(jì),主承力結(jié)構(gòu)中各桿件的橫截面形狀和材料特性均逐步明確,這時(shí)進(jìn)行結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計(jì),能夠針對(duì)明確的設(shè)計(jì)變量進(jìn)行優(yōu)化,所得優(yōu)化結(jié)果可直接用于指導(dǎo)結(jié)構(gòu)設(shè)計(jì),從而獲得更好的結(jié)構(gòu)特性。

    5 結(jié)束語(yǔ)

    目前,連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化和離散體結(jié)構(gòu)拓?fù)鋬?yōu)化在航天器工程中的應(yīng)用正在逐步普及,但大多只是單獨(dú)應(yīng)用其中一類拓?fù)鋬?yōu)化方法。本文提出針對(duì)桁架結(jié)構(gòu)的連續(xù)體-離散體兩級(jí)拓?fù)鋬?yōu)化策略,成功應(yīng)用于大型航天器桁架式承力結(jié)構(gòu)的優(yōu)化設(shè)計(jì),得到了基本滿足設(shè)計(jì)要求的桁架構(gòu)型形式,可為后續(xù)詳細(xì)設(shè)計(jì)提供參考,同時(shí)也驗(yàn)證了此種優(yōu)化策略在解決特定的桁架結(jié)構(gòu)優(yōu)化設(shè)計(jì)問題時(shí)的有效性。結(jié)構(gòu)優(yōu)化方法雖然能夠指導(dǎo)結(jié)構(gòu)設(shè)計(jì),但在優(yōu)化可行設(shè)計(jì)區(qū)域確定、連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化結(jié)果的提煉歸納、離散體結(jié)構(gòu)桁架基結(jié)構(gòu)的拓展等方面,工程經(jīng)驗(yàn)仍是必不可少的。拓?fù)鋬?yōu)化設(shè)計(jì)僅僅是得到了桁架結(jié)構(gòu)的構(gòu)型形式,在此基礎(chǔ)上要進(jìn)行進(jìn)一步的結(jié)構(gòu)參數(shù)優(yōu)化,才能確定桁架中桿件的最優(yōu)截面尺寸,達(dá)到減小結(jié)構(gòu)質(zhì)量、提高結(jié)構(gòu)性能的目的。

    (References)

    [1]陳烈民.航天器結(jié)構(gòu)與機(jī)構(gòu)[M].北京:中國(guó)科學(xué)技術(shù)出版社,2005:92-102 Chen Liemin.Spacecraft structure and mechanism[M].Beijing:China Science and Technology Press,2005:92-102(in Chinese)

    [2]澳汰爾工程軟件(上海)有限公司.澳汰爾優(yōu)化設(shè)計(jì)在航空部件設(shè)計(jì)中的應(yīng)用[J].CAD/CAM 與制造業(yè)信息化,2004(9):61-63 Altair Engineering Inc(Shanghai).Application of Altair Optimization on aerial components[J].CAD/CAM Yu Zhizao Xinxihua,2004(9):61-63(in Chinese)

    [3]周志成,曲廣吉,黃海.某衛(wèi)星平臺(tái)多結(jié)構(gòu)工況下的優(yōu)化設(shè)計(jì)[J].北京航空航天大學(xué)學(xué)報(bào),2009,35(7):821-823 Zhou Zhicheng,Qu Guangji,Huang Hai.Design optimization of satellite platform considering multiple structural cases[J].Journal of Beijing University of Aeronautics and Astronautics,2009,35(7):821-823(in Chinese)

    [4]盛聰,曾福明,濮海玲.應(yīng)用OptiStruct軟件的太陽(yáng)翼基板結(jié)構(gòu)優(yōu)化[J].航天器工程,2011,20(6):63-68 Sheng Cong,Zeng Fuming,Pu Hailing.Optimization of solar panel structures by using OptiStruct[J].Spacecraft Engineering,2011,20(6):63-68(in Chinese)

    [5]杜江華,陳昌亞,顧亦磊.拓?fù)鋬?yōu)化方法在某衛(wèi)星下艙隔板結(jié)構(gòu)設(shè)計(jì)中的應(yīng)用[C]//第一屆結(jié)構(gòu)及多學(xué)科優(yōu)化工程應(yīng)用與理論研討會(huì)論文集.北京:中國(guó)力學(xué)學(xué)會(huì),2009:1-5 Du Jianghua,Chen Changya,Gu Yilei.Application of topology optimization method to optimal design of clapboard in a satellite lower-cabin[C]//Proceedings of the 1st Structure and MDO Engineering Application and Theory Conference.Beijing:Chinese Society of Theoretical and Applied Mechanics,2009:1-5(in Chinese)

    [6]張勝蘭,鄭冬黎,郝琪,等.基于HyperWorks的結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[M].北京:機(jī)械工業(yè)出版社,2008 Zhang Shenglan,Zheng Dongli,Hao Qi,et al.Structural optimization technology based on HyperWorks[M].Beijing:China Machine Press,2008(in Chinese)

    [7]Eschenauer H A,Olhoff N.Topology optimization of continuum structures:a review[J].Appl Mech Rev,2001,54(4):331-390

    [8]Altair Engineering Inc.Altair HyperWorks version 10.0 help[EB/OL].[2013-04-16].http://www.altair.com.cn

    [9]徐芝綸.彈性力學(xué)(上冊(cè))[M].第4版.北京:高等教育出版社,2006:257-260 Xu Zhilun.Elastic mechanics(volume 1)[M].4th ed.Beijing:Higher Education Press,2006:257-260(in Chinese)

    [10]李東澤,于登云,馬興瑞.基頻約束下的桁架結(jié)構(gòu)半定規(guī)劃法拓?fù)鋬?yōu)化[J].工程力學(xué),2011,28(2):182-185 Li Dongze,Yu Dengyun,Ma Xingrui.Truss topology optimization with fundamental frequency constraints via semidefinite programming[J].Engineering Mechanics,2011,28(2):182-185(in Chinese)

    [11]周志成,曲廣吉.通信衛(wèi)星總體設(shè)計(jì)和動(dòng)力學(xué)分析[M].北京:中國(guó)科學(xué)技術(shù)出版社,2012:232-235 Zhou Zhicheng,Qu Guangji.System design and dynamics analysis of communication satellites[M].Beijing:China Science and Technology Press,2012:232-235(in Chinese)

    猜你喜歡
    力點(diǎn)貯箱桿件
    運(yùn)載火箭貯箱補(bǔ)償器結(jié)構(gòu)剛度的試驗(yàn)研究
    基于臨時(shí)支撐結(jié)構(gòu)的桿件初彎曲對(duì)其軸壓性能的影響
    四川建筑(2021年1期)2021-03-31 01:01:46
    塔式起重機(jī)拼裝式超長(zhǎng)附著桿設(shè)計(jì)與應(yīng)用
    基于Surface Evolver的推進(jìn)劑貯箱氣液界面分析
    貯箱爆炸碎片初始速度及影響因素
    貯箱輕量化設(shè)計(jì)幾何參數(shù)優(yōu)化方法
    解析重心找不著力點(diǎn)無(wú)須找
    少林與太極(2016年9期)2016-12-01 00:49:47
    KD379:便攜折疊式衣架
    某網(wǎng)架桿件彎曲的原因分析及處理
    論教師職業(yè)道德建設(shè)目標(biāo)設(shè)置的基點(diǎn)力點(diǎn)和視域
    亚洲性久久影院| 亚洲国产色片| 丝袜喷水一区| 七月丁香在线播放| 白带黄色成豆腐渣| av卡一久久| 美女脱内裤让男人舔精品视频| 午夜福利成人在线免费观看| 亚洲欧洲国产日韩| 国产真实乱freesex| 日韩人妻高清精品专区| 午夜久久久久精精品| 亚洲熟妇中文字幕五十中出| 一级毛片aaaaaa免费看小| 亚洲av熟女| 又粗又爽又猛毛片免费看| 在线播放国产精品三级| 黄片无遮挡物在线观看| 国产大屁股一区二区在线视频| 视频中文字幕在线观看| 欧美97在线视频| 我要搜黄色片| 亚洲精品乱码久久久v下载方式| 欧美一区二区精品小视频在线| 国产成人精品婷婷| 亚洲国产色片| 亚洲国产精品国产精品| 国产精品久久久久久精品电影| 亚洲无线观看免费| 国产黄片视频在线免费观看| 欧美三级亚洲精品| 直男gayav资源| 一级毛片久久久久久久久女| 国产精品1区2区在线观看.| 51国产日韩欧美| 亚洲精品aⅴ在线观看| 国产欧美日韩精品一区二区| 欧美性猛交黑人性爽| 欧美三级亚洲精品| 免费观看精品视频网站| 美女xxoo啪啪120秒动态图| 熟女电影av网| 男人狂女人下面高潮的视频| 日本av手机在线免费观看| 亚洲色图av天堂| 中文字幕制服av| 青春草亚洲视频在线观看| 国产精品电影一区二区三区| 99国产精品一区二区蜜桃av| 高清午夜精品一区二区三区| 成人av在线播放网站| 色噜噜av男人的天堂激情| 国产精品伦人一区二区| 久久精品影院6| 夜夜看夜夜爽夜夜摸| 男人狂女人下面高潮的视频| 国产一区二区亚洲精品在线观看| 一区二区三区高清视频在线| 99久久人妻综合| a级毛片免费高清观看在线播放| 永久网站在线| 久久久久久久久久久丰满| 亚洲国产精品成人综合色| 欧美又色又爽又黄视频| 亚洲人成网站高清观看| 国产精品女同一区二区软件| 美女脱内裤让男人舔精品视频| 色尼玛亚洲综合影院| 久久韩国三级中文字幕| 人体艺术视频欧美日本| 日韩一区二区三区影片| 97在线视频观看| 在线天堂最新版资源| 国产精品蜜桃在线观看| 视频中文字幕在线观看| 成人毛片60女人毛片免费| 久久久久久久亚洲中文字幕| 一级黄片播放器| 美女xxoo啪啪120秒动态图| 在线天堂最新版资源| 久久精品久久久久久久性| 国产精品一区二区三区四区久久| 国产熟女欧美一区二区| 九九爱精品视频在线观看| 久久国内精品自在自线图片| 五月伊人婷婷丁香| 亚洲人成网站在线观看播放| 国产精品熟女久久久久浪| 亚洲av成人av| 人体艺术视频欧美日本| 成年av动漫网址| 国产精品久久视频播放| 波多野结衣高清无吗| 亚洲国产日韩欧美精品在线观看| 99热全是精品| 免费看美女性在线毛片视频| 日韩av不卡免费在线播放| 老女人水多毛片| 最近的中文字幕免费完整| 亚洲av男天堂| 欧美精品国产亚洲| 中文字幕精品亚洲无线码一区| 日本wwww免费看| 免费无遮挡裸体视频| 日韩 亚洲 欧美在线| 亚洲最大成人av| 欧美又色又爽又黄视频| 美女被艹到高潮喷水动态| 亚洲精品国产av成人精品| 亚洲av福利一区| 91aial.com中文字幕在线观看| 久久久久久久久久成人| 女人被狂操c到高潮| 美女内射精品一级片tv| 村上凉子中文字幕在线| 成人一区二区视频在线观看| 男人和女人高潮做爰伦理| 精品午夜福利在线看| 看黄色毛片网站| 日韩亚洲欧美综合| 免费av毛片视频| 亚洲,欧美,日韩| 最近视频中文字幕2019在线8| 在线观看av片永久免费下载| 久久久精品大字幕| 少妇猛男粗大的猛烈进出视频 | 亚洲色图av天堂| 啦啦啦啦在线视频资源| av又黄又爽大尺度在线免费看 | www日本黄色视频网| 久久久精品大字幕| 亚洲国产精品成人综合色| 欧美又色又爽又黄视频| 精品久久久久久久久久久久久| 一本一本综合久久| 欧美性感艳星| 男女啪啪激烈高潮av片| 国产探花在线观看一区二区| 少妇人妻精品综合一区二区| 天天躁日日操中文字幕| 草草在线视频免费看| 九草在线视频观看| 亚洲av不卡在线观看| 好男人视频免费观看在线| 一边亲一边摸免费视频| 午夜日本视频在线| .国产精品久久| 久久久久久久亚洲中文字幕| 国产精品美女特级片免费视频播放器| 日韩强制内射视频| 亚洲激情五月婷婷啪啪| 欧美不卡视频在线免费观看| 真实男女啪啪啪动态图| 一夜夜www| 99久久中文字幕三级久久日本| 日本黄大片高清| 夜夜爽夜夜爽视频| 嘟嘟电影网在线观看| 人妻少妇偷人精品九色| 亚洲av福利一区| 寂寞人妻少妇视频99o| 青春草视频在线免费观看| 亚洲人与动物交配视频| 美女xxoo啪啪120秒动态图| 日韩一本色道免费dvd| 看免费成人av毛片| 精品久久久久久久人妻蜜臀av| 午夜精品国产一区二区电影 | 少妇熟女欧美另类| 久久99热这里只频精品6学生 | 狂野欧美激情性xxxx在线观看| 亚洲欧洲国产日韩| 久久久久久九九精品二区国产| 小蜜桃在线观看免费完整版高清| 一级av片app| 色5月婷婷丁香| 久久99热6这里只有精品| 亚洲av一区综合| 日本爱情动作片www.在线观看| 乱人视频在线观看| 久久久久久久久久黄片| 亚洲怡红院男人天堂| 亚洲av中文av极速乱| 最近最新中文字幕大全电影3| 久久精品久久久久久久性| 青青草视频在线视频观看| 国产黄a三级三级三级人| 国产不卡一卡二| 草草在线视频免费看| 国产精品一及| 欧美一区二区国产精品久久精品| 91在线精品国自产拍蜜月| 国产精品熟女久久久久浪| 色视频www国产| 国产精品麻豆人妻色哟哟久久 | 久久精品人妻少妇| 搡女人真爽免费视频火全软件| 国产精品人妻久久久影院| 综合色丁香网| 午夜福利视频1000在线观看| 少妇熟女欧美另类| 三级国产精品片| 久久久国产成人精品二区| 精品久久久久久成人av| 日日摸夜夜添夜夜添av毛片| 午夜福利高清视频| 在线观看一区二区三区| 成人特级av手机在线观看| 色播亚洲综合网| 国产精品蜜桃在线观看| 在线播放无遮挡| 国产黄片美女视频| 只有这里有精品99| 国产男人的电影天堂91| 九九热线精品视视频播放| 久久精品久久精品一区二区三区| 亚洲av免费高清在线观看| 韩国av在线不卡| 亚洲人成网站在线观看播放| 久久久久国产网址| 欧美三级亚洲精品| 97超视频在线观看视频| 高清av免费在线| 在线观看美女被高潮喷水网站| 久久久色成人| 国产在视频线在精品| 欧美日韩国产亚洲二区| 午夜福利在线观看免费完整高清在| 蜜臀久久99精品久久宅男| 久久婷婷人人爽人人干人人爱| 狂野欧美激情性xxxx在线观看| videossex国产| 欧美色视频一区免费| 成人二区视频| 日本一本二区三区精品| 欧美激情国产日韩精品一区| 国产午夜福利久久久久久| 国产精品国产三级专区第一集| 亚洲怡红院男人天堂| 亚洲av电影在线观看一区二区三区 | 丝袜喷水一区| a级一级毛片免费在线观看| 一个人看视频在线观看www免费| 伦理电影大哥的女人| 久久精品91蜜桃| 汤姆久久久久久久影院中文字幕 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | a级一级毛片免费在线观看| 国产成人午夜福利电影在线观看| 高清视频免费观看一区二区 | 国产精品av视频在线免费观看| 国产亚洲最大av| 精品久久久久久久末码| 韩国av在线不卡| 级片在线观看| 久久亚洲国产成人精品v| 女人久久www免费人成看片 | 国产精品蜜桃在线观看| 女的被弄到高潮叫床怎么办| 夜夜爽夜夜爽视频| 亚洲在久久综合| 欧美不卡视频在线免费观看| 日韩欧美精品v在线| 麻豆国产97在线/欧美| 欧美性猛交黑人性爽| 好男人视频免费观看在线| 如何舔出高潮| 色5月婷婷丁香| 寂寞人妻少妇视频99o| 欧美日韩在线观看h| 欧美一区二区国产精品久久精品| 免费搜索国产男女视频| 极品教师在线视频| 99久久成人亚洲精品观看| 亚洲国产日韩欧美精品在线观看| 日本-黄色视频高清免费观看| 在线观看av片永久免费下载| 久久精品熟女亚洲av麻豆精品 | 欧美激情国产日韩精品一区| 久久久久久久久中文| av福利片在线观看| 国产精品一区二区在线观看99 | 一级黄色大片毛片| 汤姆久久久久久久影院中文字幕 | 色综合亚洲欧美另类图片| 国产亚洲午夜精品一区二区久久 | 日本wwww免费看| 变态另类丝袜制服| 一区二区三区四区激情视频| 成人av在线播放网站| 欧美丝袜亚洲另类| 国语自产精品视频在线第100页| 国产精品综合久久久久久久免费| 午夜福利高清视频| 午夜a级毛片| 亚洲欧洲日产国产| 亚洲精品影视一区二区三区av| 国产美女午夜福利| 久久久a久久爽久久v久久| 国产白丝娇喘喷水9色精品| 国内精品美女久久久久久| 国产成人91sexporn| 好男人视频免费观看在线| 成年免费大片在线观看| 久久久精品大字幕| 国内精品一区二区在线观看| 中文字幕av成人在线电影| 中文精品一卡2卡3卡4更新| 欧美性猛交╳xxx乱大交人| 久久久久网色| 亚洲天堂国产精品一区在线| 国产又色又爽无遮挡免| 成人亚洲欧美一区二区av| 日韩一本色道免费dvd| 在线观看一区二区三区| 欧美精品国产亚洲| 国产黄片视频在线免费观看| 成人无遮挡网站| 成人性生交大片免费视频hd| 免费看av在线观看网站| 中文字幕av在线有码专区| av在线蜜桃| 久久精品国产自在天天线| 免费大片18禁| 乱系列少妇在线播放| 又爽又黄a免费视频| 高清视频免费观看一区二区 | 搡女人真爽免费视频火全软件| 久久人人爽人人爽人人片va| 国产亚洲91精品色在线| 卡戴珊不雅视频在线播放| 少妇人妻精品综合一区二区| 观看免费一级毛片| 中文亚洲av片在线观看爽| 日韩精品青青久久久久久| 亚洲精品,欧美精品| 美女黄网站色视频| 日韩av在线大香蕉| 久久久久网色| 国产久久久一区二区三区| 日韩欧美在线乱码| 亚洲第一区二区三区不卡| 美女大奶头视频| 一级毛片电影观看 | 精品久久久久久久久久久久久| 看黄色毛片网站| 国产高清不卡午夜福利| 国产三级中文精品| 国产一区亚洲一区在线观看| 亚洲va在线va天堂va国产| 日韩在线高清观看一区二区三区| 久久精品影院6| 国产在视频线在精品| 精品人妻视频免费看| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影| 性插视频无遮挡在线免费观看| 欧美激情国产日韩精品一区| 国产精品精品国产色婷婷| 日韩视频在线欧美| 亚洲av男天堂| 国产乱人偷精品视频| 亚洲精品成人久久久久久| 久久精品国产亚洲网站| 国产乱人视频| 国产美女午夜福利| 久久热精品热| 水蜜桃什么品种好| 午夜福利在线观看免费完整高清在| 六月丁香七月| 国产精品国产高清国产av| 欧美性猛交╳xxx乱大交人| 秋霞伦理黄片| 色哟哟·www| 亚洲人成网站高清观看| 亚洲精品乱久久久久久| 又爽又黄a免费视频| 最新中文字幕久久久久| 午夜激情福利司机影院| 国产免费一级a男人的天堂| 久久韩国三级中文字幕| 少妇高潮的动态图| 精品免费久久久久久久清纯| 精品少妇黑人巨大在线播放 | 人人妻人人澡人人爽人人夜夜 | 国产精品熟女久久久久浪| 色网站视频免费| 91精品伊人久久大香线蕉| 91精品一卡2卡3卡4卡| 国产成人91sexporn| 国产精品蜜桃在线观看| 美女大奶头视频| 国产亚洲一区二区精品| 午夜福利高清视频| 搡老妇女老女人老熟妇| 水蜜桃什么品种好| 成年版毛片免费区| 久久久精品94久久精品| 热99re8久久精品国产| 黄色欧美视频在线观看| 欧美不卡视频在线免费观看| 99久久精品热视频| 最近中文字幕高清免费大全6| 在线播放国产精品三级| 色噜噜av男人的天堂激情| 一级黄片播放器| 午夜精品一区二区三区免费看| 国产精品综合久久久久久久免费| 淫秽高清视频在线观看| 乱系列少妇在线播放| 色综合色国产| 亚洲av成人av| 狂野欧美白嫩少妇大欣赏| av黄色大香蕉| 国产精品一区二区三区四区免费观看| 日韩欧美三级三区| 好男人视频免费观看在线| 99热这里只有精品一区| 黄色日韩在线| av播播在线观看一区| 校园人妻丝袜中文字幕| 免费黄网站久久成人精品| 一个人看视频在线观看www免费| 国产三级在线视频| 亚洲va在线va天堂va国产| 精品无人区乱码1区二区| 看黄色毛片网站| 看十八女毛片水多多多| 免费看美女性在线毛片视频| 亚洲自拍偷在线| 舔av片在线| 国产乱来视频区| 少妇的逼水好多| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 免费在线观看成人毛片| 中文字幕人妻熟人妻熟丝袜美| 国产一区二区三区av在线| 69av精品久久久久久| 国产在视频线精品| 欧美+日韩+精品| 婷婷色麻豆天堂久久 | 国产一级毛片在线| 国产精品人妻久久久久久| 亚洲精品日韩在线中文字幕| 超碰av人人做人人爽久久| av卡一久久| 亚洲精华国产精华液的使用体验| 国产精品国产三级专区第一集| 免费大片18禁| 男女那种视频在线观看| 九九在线视频观看精品| 最近最新中文字幕免费大全7| 国产亚洲91精品色在线| 少妇裸体淫交视频免费看高清| 最近的中文字幕免费完整| 九草在线视频观看| av又黄又爽大尺度在线免费看 | 国产精品蜜桃在线观看| 国产熟女欧美一区二区| АⅤ资源中文在线天堂| 人妻少妇偷人精品九色| 亚洲在线观看片| 村上凉子中文字幕在线| 久久99热这里只频精品6学生 | 日本黄大片高清| 婷婷色综合大香蕉| 日韩欧美精品免费久久| 又黄又爽又刺激的免费视频.| 日韩欧美 国产精品| 美女国产视频在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲成人中文字幕在线播放| 18禁在线无遮挡免费观看视频| 老女人水多毛片| 日本-黄色视频高清免费观看| 亚洲精品国产成人久久av| 午夜激情欧美在线| 看片在线看免费视频| 午夜福利在线在线| 男女视频在线观看网站免费| 日韩欧美国产在线观看| 成人毛片a级毛片在线播放| 免费看美女性在线毛片视频| 一级av片app| 亚洲人成网站在线观看播放| 日韩精品有码人妻一区| 综合色丁香网| 亚洲美女搞黄在线观看| 国内精品美女久久久久久| 中文精品一卡2卡3卡4更新| av女优亚洲男人天堂| 亚洲欧美日韩东京热| 91午夜精品亚洲一区二区三区| 色哟哟·www| 国产69精品久久久久777片| 久久久久久久久久久丰满| 午夜日本视频在线| 免费观看性生交大片5| 99久久人妻综合| 免费观看在线日韩| 菩萨蛮人人尽说江南好唐韦庄 | 一卡2卡三卡四卡精品乱码亚洲| 亚洲国产高清在线一区二区三| 有码 亚洲区| 99久久成人亚洲精品观看| 人体艺术视频欧美日本| a级毛片免费高清观看在线播放| 天堂影院成人在线观看| 男的添女的下面高潮视频| 91久久精品国产一区二区成人| 伦精品一区二区三区| 亚洲欧美精品专区久久| 国语对白做爰xxxⅹ性视频网站| 国产伦精品一区二区三区四那| 日韩制服骚丝袜av| 久久久久九九精品影院| 国产精品乱码一区二三区的特点| 亚州av有码| 最近2019中文字幕mv第一页| 久久久久免费精品人妻一区二区| 免费电影在线观看免费观看| 精品人妻熟女av久视频| 免费看光身美女| 午夜日本视频在线| 一本一本综合久久| 精品久久久久久久末码| 亚洲av成人精品一区久久| 国产av在哪里看| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 日韩制服骚丝袜av| 女人被狂操c到高潮| 特级一级黄色大片| 我要搜黄色片| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 久久国内精品自在自线图片| h日本视频在线播放| 日韩成人av中文字幕在线观看| 日韩欧美在线乱码| 男女啪啪激烈高潮av片| 亚洲欧美清纯卡通| 欧美97在线视频| 亚洲av福利一区| 国产精品蜜桃在线观看| 欧美日本亚洲视频在线播放| 最近的中文字幕免费完整| 波野结衣二区三区在线| 日本爱情动作片www.在线观看| 啦啦啦观看免费观看视频高清| 七月丁香在线播放| 超碰97精品在线观看| 精品一区二区三区人妻视频| 日本av手机在线免费观看| 免费人成在线观看视频色| 91aial.com中文字幕在线观看| 在线观看美女被高潮喷水网站| 日韩在线高清观看一区二区三区| 亚洲最大成人av| 国产伦精品一区二区三区四那| 亚洲国产精品专区欧美| 成人鲁丝片一二三区免费| 夫妻性生交免费视频一级片| 最后的刺客免费高清国语| 国产成人一区二区在线| 激情 狠狠 欧美| 久久人人爽人人片av| 日本五十路高清| 久久精品综合一区二区三区| av在线老鸭窝| 国内揄拍国产精品人妻在线| 中文字幕精品亚洲无线码一区| 两性午夜刺激爽爽歪歪视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲天堂国产精品一区在线| 高清日韩中文字幕在线| 高清在线视频一区二区三区 | 2021天堂中文幕一二区在线观| 亚洲人与动物交配视频| 亚洲欧美成人精品一区二区| 日本午夜av视频| 晚上一个人看的免费电影| 亚洲国产高清在线一区二区三| 亚洲aⅴ乱码一区二区在线播放| 国产高清三级在线| 搡女人真爽免费视频火全软件| 日本黄色片子视频| 午夜福利高清视频| 能在线免费看毛片的网站| 国产黄片美女视频| 一夜夜www| 国产综合懂色| av免费在线看不卡| 国产精品国产三级国产专区5o | 最近中文字幕高清免费大全6| 乱码一卡2卡4卡精品| 亚洲国产精品专区欧美| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 国产黄片视频在线免费观看| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 久99久视频精品免费| 国产高清三级在线| 日韩中字成人| 久久精品国产99精品国产亚洲性色| 日韩,欧美,国产一区二区三区 | 国产成人精品一,二区| 国产欧美另类精品又又久久亚洲欧美| 又爽又黄a免费视频| 波多野结衣高清无吗| 欧美色视频一区免费| 网址你懂的国产日韩在线| 性插视频无遮挡在线免费观看|