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

    水面艦船航行性能MDO系統(tǒng)頂層設(shè)計(jì)研究

    2012-06-22 05:07:36程素斌韋喜忠李勝忠
    船舶力學(xué) 2012年11期
    關(guān)鍵詞:艦船航行水面

    趙 峰 ,程素斌 ,楊 磊 ,韋喜忠 ,李勝忠

    (1中國船舶科學(xué)研究中心,江蘇無錫 214082;2同濟(jì)大學(xué),上海 20092)

    水面艦船航行性能MDO系統(tǒng)頂層設(shè)計(jì)研究

    趙 峰1,程素斌2,楊 磊1,韋喜忠1,李勝忠1

    (1中國船舶科學(xué)研究中心,江蘇無錫 214082;2同濟(jì)大學(xué),上海 20092)

    水面艦船航行性能設(shè)計(jì)是一個(gè)典型的多學(xué)科設(shè)計(jì)優(yōu)化問題,適宜于采用基于系統(tǒng)工程思想的MDO創(chuàng)新設(shè)計(jì)理念和方法。文章首先分析了水面艦船航行性能設(shè)計(jì)領(lǐng)域開展MDO應(yīng)用研究的重要意義,之后闡述了MDO概念內(nèi)涵;制定了總體研究思路,確定了系統(tǒng)設(shè)計(jì)目標(biāo)、技術(shù)指標(biāo)及頂層設(shè)計(jì)技術(shù)方案,剖析了需要首先解決的若干關(guān)鍵技術(shù),并給出關(guān)鍵技術(shù)的攻關(guān)建議。

    水面艦船;航行性能;MDO;頂層設(shè)計(jì)

    1 引 言

    航行性能是水面艦船的基本性能,它直接影響著艦船的作戰(zhàn)效能和海上生存能力。因此,航行性能設(shè)計(jì)一直被作為水面艦船總體設(shè)計(jì)的核心基礎(chǔ)內(nèi)容。

    航行性能設(shè)計(jì)是一個(gè)系統(tǒng)性很強(qiáng)的復(fù)雜過程,它涉及多個(gè)學(xué)科,這些學(xué)科高度關(guān)聯(lián)、相互影響,使得直接進(jìn)行系統(tǒng)綜合設(shè)計(jì)相當(dāng)困難。傳統(tǒng)上將其分解為快速性、耐波性和操縱性三大子學(xué)科,采用序列式設(shè)計(jì)方法,設(shè)計(jì)優(yōu)化在各子學(xué)科內(nèi)部進(jìn)行,各子學(xué)科之間的交互作用通常用一組經(jīng)驗(yàn)性的“靜態(tài)常數(shù)”來描述,最后將各學(xué)科的設(shè)計(jì)結(jié)果裝配起來,通過設(shè)計(jì)師的綜合協(xié)調(diào)得到最終方案。這種傳統(tǒng)的設(shè)計(jì)模式中,設(shè)計(jì)師試圖利用先驗(yàn)的、經(jīng)驗(yàn)性的“靜態(tài)常數(shù)”來協(xié)調(diào)學(xué)科間的動(dòng)態(tài)交互作用,造成學(xué)科設(shè)計(jì)環(huán)境的“失真”;經(jīng)“裝配”和再協(xié)調(diào)使得各子學(xué)科的性能被進(jìn)一步“中和”,勢必大幅降低設(shè)計(jì)方案的整體性能。基于系統(tǒng)工程概念的多學(xué)科設(shè)計(jì)優(yōu)化(Multidisciplinary Design Optimization,MDO)技術(shù)以求解多系統(tǒng)的耦合協(xié)同為核心,利用有效的MDO策略來重新規(guī)劃、組織復(fù)雜系統(tǒng)的設(shè)計(jì)過程,有效克服傳統(tǒng)設(shè)計(jì)模式的缺陷,能夠獲得系統(tǒng)整體最優(yōu)解。開展水面艦船航行性能MDO應(yīng)用研究對提高艦船總體設(shè)計(jì)水平和質(zhì)量,縮短設(shè)計(jì)周期、降低研制成本和設(shè)計(jì)風(fēng)險(xiǎn)有著重要意義。

    本研究以水面艦船為主要應(yīng)用對象,以加強(qiáng)基礎(chǔ)、提高能力為著眼點(diǎn),按照系統(tǒng)工程思想,提出水面艦船航行性能MDO系統(tǒng)頂層策劃思路,確定了系統(tǒng)設(shè)計(jì)目標(biāo)、技術(shù)指標(biāo)及頂層設(shè)計(jì)技術(shù)方案,分析了需要首先解決的若干關(guān)鍵技術(shù),并給出了關(guān)鍵技術(shù)的攻關(guān)建議。

    2 MDO應(yīng)用研究現(xiàn)狀

    MDO作為一門單獨(dú)的研究領(lǐng)域,始于20世紀(jì)80年代后期。1982年Sobieski在研究大型結(jié)構(gòu)優(yōu)化問題時(shí),首次提出了MDO設(shè)想[1],并在其后發(fā)表的一系列論文中對MDO問題做了進(jìn)一步的闡述[2-3]。MDO概念一經(jīng)提出便引起學(xué)術(shù)界的極大關(guān)注,并率先在航空航天飛行器設(shè)計(jì)領(lǐng)域得到成功應(yīng)用。美國航空航天局(NASA)在“塞式噴管的MDO研究計(jì)劃”中,曾采用傳統(tǒng)串行設(shè)計(jì)方法和MDO方法對某塞式噴管發(fā)動(dòng)機(jī)進(jìn)行對比設(shè)計(jì),結(jié)果表明,后者獲得的最佳方案其最小起始推重比較前者降低4%;在其高速民用飛機(jī)(HSCT)研究項(xiàng)目中,同時(shí)考慮了271個(gè)設(shè)計(jì)變量、31 868個(gè)約束。

    MDO技術(shù)對復(fù)雜工程系統(tǒng)的設(shè)計(jì)能力掀起了MDO研究的熱潮,美國、歐洲、俄羅斯、日本、韓國等紛紛啟動(dòng)各自的MDO研究項(xiàng)目,MDO的應(yīng)用領(lǐng)域已從航空航天擴(kuò)展到汽車、通訊、機(jī)械、醫(yī)療、建筑、運(yùn)輸?shù)阮I(lǐng)域[4]。

    在艦船設(shè)計(jì)領(lǐng)域,美國新一代CVN-78“福特級”航母、DDG-1000“朱姆沃特”級驅(qū)逐艦和SSN-774“弗吉尼亞”級攻擊型核潛艇開始應(yīng)用頂層設(shè)計(jì)思想,其設(shè)計(jì)理念與MDO思想已非常接近;美國在新型驅(qū)逐艦DD-21的概念設(shè)計(jì)階段開始應(yīng)用MDO技術(shù)進(jìn)行總體優(yōu)化設(shè)計(jì);Hart[5]將MDO與PSO優(yōu)化算法相結(jié)合,應(yīng)用于艦船概念設(shè)計(jì)優(yōu)化和復(fù)雜工程系統(tǒng)的成本預(yù)算上,收到了不錯(cuò)的效果。

    MDO技術(shù)的深入研究和應(yīng)用,推動(dòng)了iSIGHT、ModelCenter、Optimus等商業(yè)軟件系統(tǒng)的發(fā)展,加速了有針對性的專業(yè)MDO集成平臺的開發(fā),如美國在HSTC項(xiàng)目中,針對飛行器總體設(shè)計(jì)優(yōu)化開發(fā)了FIDO、CJOpt等MDO軟件框架,俄羅斯科學(xué)研究院開發(fā)了“多學(xué)科設(shè)計(jì)優(yōu)化環(huán)境(Indirect Optimization on the basis of Self-Organization,IOSO)”[6]等。

    國內(nèi)對MDO技術(shù)及應(yīng)用研究的進(jìn)展情況進(jìn)行了積極的跟蹤研究[7],在MDO理論,如MDO算法、MDO建模、分解策略等方面開展了一些探索性的研究[8-10],并嘗試在一些簡單的小型系統(tǒng)設(shè)計(jì)中應(yīng)用,如陳小前將基于響應(yīng)面的MDO算法應(yīng)用于飛行器的概念設(shè)計(jì)中[11];陳琪鋒將協(xié)同優(yōu)化算法應(yīng)用于導(dǎo)彈總體參數(shù)化設(shè)計(jì)和衛(wèi)星系統(tǒng)的設(shè)計(jì)中[12]。在艦船設(shè)計(jì)領(lǐng)域,MDO的應(yīng)用研究也在逐步受到重視,并在潛器、水面艦船、海洋平臺、魚雷等方面開展了一些探索性研究工作,如劉蔚等提出大深度載人潛水器多學(xué)科設(shè)計(jì)優(yōu)化的初步思考[13],并建立了基于BLH的MDO模型[14];趙敏等結(jié)合BLISS方法和CO方法的特點(diǎn),提出BLISCO算法,并嘗試用于大深度載人潛水器的設(shè)計(jì)優(yōu)化[15];潘彬彬等分析了船舶概念設(shè)計(jì)階段子系統(tǒng)的劃分,總結(jié)了船舶設(shè)計(jì)中應(yīng)用多學(xué)科設(shè)計(jì)優(yōu)化方法應(yīng)包含的研究內(nèi)容[16],并提出了一種船舶概念設(shè)計(jì)階段的MDO模型[17];趙敏[18]、姜哲[19]將MDO方法應(yīng)用于桁架式Spar平臺的方案設(shè)計(jì)。趙加鵬提出基于并行子空間的魚雷外形設(shè)計(jì)MDO模型[20];馮百威開展了艦船多學(xué)科綜合優(yōu)化設(shè)計(jì)計(jì)算環(huán)境研究[21]。從發(fā)表的文獻(xiàn)來看,國內(nèi)對MDO技術(shù)及應(yīng)用研究相當(dāng)活躍,尤以航空、航天領(lǐng)域最為突出;艦船設(shè)計(jì)領(lǐng)域的研究雖相對滯后,但已處在持續(xù)升溫狀態(tài);但國內(nèi)對MDO基礎(chǔ)理論及相關(guān)技術(shù)的研究缺乏系統(tǒng)性和長效性,以簡單模型研究為主,鮮有實(shí)例驗(yàn)證,與工程應(yīng)用的要求尚有一定差距。以工程應(yīng)用為目標(biāo),以水面艦船航行性能為對象的MDO應(yīng)用系統(tǒng)性研究,在國內(nèi)尚屬空白。

    MDO作為一種正在崛起的新興工程學(xué)科,將成為優(yōu)化設(shè)計(jì)的大趨勢。MDO的創(chuàng)新設(shè)計(jì)優(yōu)化理念和方法契合現(xiàn)代水面艦船總體設(shè)計(jì)對于技術(shù)創(chuàng)新和跨越式發(fā)展的迫切需求,也響應(yīng)船舶工業(yè)發(fā)展對于提高設(shè)計(jì)能力、縮短研制周期的需求。

    3 MDO概念內(nèi)涵

    3.1 MDO定義及數(shù)學(xué)模型

    目前,在學(xué)術(shù)界,對MDO的定義尚未完全統(tǒng)一,美國航天航空學(xué)會(huì)(AIAA)的MDO技術(shù)委員會(huì)從不同的角度給出了三種定義,在此援引被認(rèn)為相對嚴(yán)謹(jǐn)、引用次數(shù)較多的一種來定義MDO問題。

    MDO定義:MDO是一種通過充分探索和利用系統(tǒng)中相互作用的協(xié)同機(jī)制來設(shè)計(jì)復(fù)雜系統(tǒng)和子系統(tǒng)的方法論[22]。

    MDO方法的基本指導(dǎo)思想是利用MDO策略來組織和管理復(fù)雜系統(tǒng)及其子系統(tǒng)的優(yōu)化設(shè)計(jì)過程,通過分解和協(xié)調(diào)來實(shí)現(xiàn)多系統(tǒng)的協(xié)同優(yōu)化,MDO技術(shù)實(shí)質(zhì)上是一種設(shè)計(jì)方法學(xué)。在MDO系統(tǒng)中,每一個(gè)子系統(tǒng)可用圖1來描述。

    圖1中:X、Xi:設(shè)計(jì)變量,其中X為系統(tǒng)設(shè)計(jì)變量,Xi為學(xué)科i的局部設(shè)計(jì)變量;Yji:耦合狀態(tài)變量,由學(xué)科j輸出到學(xué)科i的量;(j=1,2,…,n);Yi:學(xué)科i的輸出,它是 X、Xi、Yji的函數(shù),可表達(dá)為: Yi=Y(jié)( X , Xi,Yji)(j= 1 ,2,…,n );hi、gi: 約束條件,h 表示等式約束,g表示不等式約束。

    則多學(xué)科設(shè)計(jì)優(yōu)化問題可表述為:

    圖1 MDO系統(tǒng)的系統(tǒng)模型Fig.1 Subsystem model of the MDO system

    其中:f( f1, f2,…,fn):目標(biāo)函數(shù); hi(X ,Xi,Yi):學(xué)科 i等式約束;gi(X ,Xi,Yi):學(xué)科 i不等式約束;Ei():學(xué)科i一致性約束。

    從MDO問題的數(shù)學(xué)描述及其子系統(tǒng)模型可以看出,由于MDO過程中要協(xié)調(diào)處理學(xué)科間的動(dòng)態(tài)交互作用,使得MDO問題的復(fù)雜度、計(jì)算量將顯著大于幾個(gè)學(xué)科設(shè)計(jì)過程的簡單疊加,如何降低系統(tǒng)的復(fù)雜度、減少計(jì)算量是MDO技術(shù)及應(yīng)用需要重點(diǎn)研究的內(nèi)容。

    MDO作為一種新興的工程學(xué)科,其理論和方法仍在不斷發(fā)展,根據(jù)MDO技術(shù)委員會(huì)的分類方法,MDO的主要研究內(nèi)容如表1所列[4]。

    表1 MDO主要研究內(nèi)容Tab.1 Primary research content of the MDO

    3.2 MDO與其他優(yōu)化過程的區(qū)別

    近年來,國內(nèi)MDO應(yīng)用研究相當(dāng)活躍,見諸報(bào)道的文獻(xiàn)很多,但是對MDO的理解差異較大,MDO與涉及多個(gè)學(xué)科的優(yōu)化及多目標(biāo)優(yōu)化概念有相互混淆的趨勢。在此,本文以舵的綜合優(yōu)化為例來闡述本文所述的MDO概念。

    如圖2所示,設(shè)計(jì)目標(biāo)是使舵葉升阻比和結(jié)構(gòu)強(qiáng)度達(dá)到最大,選擇舵葉高度、寬度、剖面幾何形狀及舵葉結(jié)構(gòu)厚度作為設(shè)計(jì)變量。作用于舵葉上的水動(dòng)力載荷將引起舵葉結(jié)構(gòu)變形,而結(jié)構(gòu)變形又反過來影響舵葉的水動(dòng)力載荷分布。顯然這是一個(gè)涉及兩個(gè)學(xué)科的耦合系統(tǒng),但卻有兩種解決策略:若在優(yōu)化過程中忽略兩個(gè)學(xué)科之間的影響(如圖3左圖所示),則它僅是一個(gè)涉及多個(gè)學(xué)科的多目標(biāo)過程;而只有在優(yōu)化過程中考慮了兩個(gè)學(xué)科之間的動(dòng)態(tài)交互作用,它才是一個(gè)MDO過程(圖3右圖所示)。

    圖2 舵葉幾何形狀Fig.2 Shape of rudder

    圖3 優(yōu)化過程(左:多目標(biāo)過程;右:MDO過程)Fig.3 Optimization process(left:Multi-objective process;right:MDO process)

    由此可以看出,判斷一個(gè)設(shè)計(jì)優(yōu)化過程是否是MDO過程不是看它是否涉及到不同的學(xué)科領(lǐng)域,而是要看它是否應(yīng)用了MDO方法對學(xué)科之間的耦合關(guān)系進(jìn)行了協(xié)調(diào)處理。

    4 水面艦船航行性能MDO系統(tǒng)頂層設(shè)計(jì)

    4.1 設(shè)計(jì)思路

    本研究的功能需求目標(biāo)是構(gòu)建面向工程應(yīng)用、以快速性、耐波性和操縱性為基本學(xué)科內(nèi)涵的水面艦船航行性能MDO集成設(shè)計(jì)平臺,為艦船總體性能優(yōu)化設(shè)計(jì)以及新船型開發(fā)提供先進(jìn)的技術(shù)手段和基礎(chǔ)能力。由于航行性能各學(xué)科領(lǐng)域的設(shè)計(jì)、分析及預(yù)報(bào)技術(shù)在不斷發(fā)展,MDO技術(shù)本身也在不斷發(fā)展,艦船對航行性能的要求會(huì)隨著艦船裝備的提升而不斷發(fā)展,因此水面艦船航行性能MDO系統(tǒng)是一個(gè)不斷發(fā)展逐步完善的過程。水面艦船航行性能MDO系統(tǒng)頂層設(shè)計(jì)思路如圖4所示,其中需求分析、研究現(xiàn)狀及存在的問題在上文已有詳細(xì)描述;研究的目標(biāo)在于應(yīng)用MDO策略重新規(guī)劃、組織航行性能設(shè)計(jì)優(yōu)化過程,建立科學(xué)的、基于系統(tǒng)工程思想的學(xué)科耦合協(xié)調(diào)機(jī)制,實(shí)現(xiàn)單學(xué)科設(shè)計(jì)向多學(xué)科綜合優(yōu)化設(shè)計(jì)的歷史轉(zhuǎn)變。

    圖4 水面艦船航行性能MDO系統(tǒng)頂層設(shè)計(jì)思路Fig.4 General road map for top-layer design on the MDO system of ship navigational performance

    4.2 頂層設(shè)計(jì)技術(shù)方案

    依據(jù)上文需求分析情況及設(shè)計(jì)目標(biāo)的要求,水面艦船航行性能頂層設(shè)計(jì)技術(shù)方案如圖5所示,整個(gè)研究規(guī)劃分為四個(gè)部分:第一部分為綜合航行性能評價(jià)指標(biāo)體系研究及設(shè)計(jì)問題定義,實(shí)際上是從工程應(yīng)用的角度解析總體設(shè)計(jì)對航行性能的需求,以及如何在MDO系統(tǒng)中描述或者說定義這些需求,要求能夠體現(xiàn)航行性能設(shè)計(jì)一般性要求,同時(shí)也能定對具體對象的特殊性能設(shè)計(jì)要求;第二部分為MDO技術(shù)及應(yīng)用研究,主要包括MDO基本理論研究、MDO算法研究、MDO算法應(yīng)用研究、MDO過程中的近似技術(shù)應(yīng)用研究以及MDO建模技術(shù)研究,這是本方法的核心及區(qū)別于其他設(shè)計(jì)優(yōu)化過程的關(guān)鍵所在;第三部分為學(xué)科設(shè)計(jì)優(yōu)化共性技術(shù)研究,主要包括性能分析評估技術(shù)、船型幾何表達(dá)及重構(gòu)技術(shù)研究、設(shè)計(jì)空間搜索策略 (優(yōu)化算法)技術(shù),這三個(gè)部分可以組成一個(gè)相對獨(dú)立的子學(xué)科優(yōu)化設(shè)計(jì)過程,在MDO策略的組織下,實(shí)現(xiàn)各學(xué)科之間的協(xié)同優(yōu)化;第四部分是綜合集成及集成平臺的應(yīng)用與驗(yàn)證,主要是將上述研究成果以知識化的形式凝練起來,形成可應(yīng)用的手段。主要研究內(nèi)容分述如下:

    圖5 水面艦船航行性能MDO系統(tǒng)技術(shù)方案Fig.5 Technical scheme of the MDO system of ship navigational performance

    4.3 關(guān)鍵技術(shù)分析及攻關(guān)策略

    4.3.1 關(guān)鍵技術(shù)分析

    水面艦船航行性能MDO系統(tǒng)研究涉及MDO技術(shù)、水面艦船航行性能分析評估技術(shù)、船型設(shè)計(jì)技術(shù)、優(yōu)化技術(shù)、計(jì)算機(jī)技術(shù)等眾多學(xué)科領(lǐng)域,研究范圍廣,需要突破的關(guān)鍵技術(shù)多。根據(jù)國內(nèi)外在水面艦船航行性能領(lǐng)域MDO應(yīng)用技術(shù)研究現(xiàn)狀分析,當(dāng)前需要率先解決的關(guān)鍵技術(shù)可總結(jié)如下:

    (1)水面艦船綜合航行性能評價(jià)體系構(gòu)建;

    (2)水面艦船航行性能MDO建模技術(shù);

    (3)水面艦船航行性能MDO優(yōu)化算法應(yīng)用技術(shù);

    (4)近似技術(shù);

    (5)靈敏度分析技術(shù);

    (6)航行性能分析評估技術(shù);

    (7)船體幾何參數(shù)化表達(dá)及重構(gòu)技術(shù);

    (8)水面艦船航行性能MDO計(jì)算環(huán)境架構(gòu)。

    4.3.2 關(guān)鍵技術(shù)攻關(guān)策略

    以下對4.3.2節(jié)中關(guān)鍵技術(shù)的難點(diǎn)進(jìn)行分析,從實(shí)用性和可行性的角度提出攻關(guān)建議。

    (1)水面艦船綜合航行性能評價(jià)體系構(gòu)建

    水面艦船綜合航行性能評價(jià)體系構(gòu)建的困難在于艦船任務(wù)的多樣性、作業(yè)環(huán)境的復(fù)雜性使得不同的艦船對航行性能的要求各不相同,構(gòu)建統(tǒng)一的評價(jià)體系相當(dāng)困難;其次航行性能評價(jià)指標(biāo)眾多,且目標(biāo)不一、標(biāo)準(zhǔn)不一,評價(jià)方法之間還可能存在相互沖突的情況,從中提煉能夠反映綜合航行性能的評價(jià)指標(biāo)族相當(dāng)困難。

    構(gòu)建水面艦船航行性能評價(jià)體系應(yīng)該與艦船總體設(shè)計(jì)統(tǒng)一考慮,參考現(xiàn)行艦船總體設(shè)計(jì)過程中采用的性能評價(jià)方法,采取分步走的策略,按基本指標(biāo)族和較完整的指標(biāo)體系兩個(gè)步驟來實(shí)施??紤]到指標(biāo)分配也是一個(gè)需要不斷權(quán)衡、不斷優(yōu)化的過程,也可以應(yīng)用MDO技術(shù)來求解。

    (2)水面艦船航行性能MDO建模技術(shù)

    MDO建模是多學(xué)科設(shè)計(jì)優(yōu)化過程的基礎(chǔ),水面艦船航行性能MDO建模的難點(diǎn)在于對航行性能各子系統(tǒng)之間耦合關(guān)系的簡化和抽象及其數(shù)學(xué)模型的構(gòu)建,系統(tǒng)級和子系統(tǒng)級設(shè)計(jì)目標(biāo)、設(shè)計(jì)變量的選擇,同時(shí)還要考慮問題求解的復(fù)雜程度問題。根據(jù)水面艦船航行性能研究、設(shè)計(jì)組織體系及性能設(shè)計(jì)的特點(diǎn),可采用如圖6所示的建模策略。

    (3)水面艦船航行性能MDO算法應(yīng)用技術(shù)研究

    MDO算法是MDO技術(shù)及應(yīng)用最核心的部分,是MDO研究最活躍的領(lǐng)域。水面艦船航行性能MDO系統(tǒng)設(shè)計(jì)的重點(diǎn)在于對MDO算法的應(yīng)用和發(fā)展,其難點(diǎn)在于采用何種MDO策略協(xié)調(diào)各學(xué)科之間的耦合關(guān)系才是最有效的,以及如何提高M(jìn)DO算法的可靠性及收斂速度,如何增強(qiáng)系統(tǒng)自治性及降低信息傳遞量,以滿足工程設(shè)計(jì)的要求。

    圖6 水面艦船航行性能MDO建模策略Fig.6 Modeling strategy of the MDO system of ship navigational performance

    按照系統(tǒng)分解的層次關(guān)系,MDO算法可分為兩類:單級MDO算法和多級MDO算法。常見的單級優(yōu)化算法包括標(biāo)準(zhǔn)單級優(yōu)化算法[23]、基于全局敏度方程(Global Sensitivity Equations,GSE)的單級優(yōu)化算法[3]和同時(shí)設(shè)計(jì)與分析單級優(yōu)化算法[24]。這類方法適合于求解設(shè)計(jì)變量少、學(xué)科不多和計(jì)算量不大的簡單系統(tǒng),對于航行性能這樣的復(fù)雜系統(tǒng)將難以適用。

    多級優(yōu)化算法主要包括并行子空間優(yōu)化算法(Concurrent Subspace Optimization,CSSO)[25]、協(xié)同優(yōu)化算法 (Collaborative Optimization, CO)[17]和兩級系統(tǒng)綜合優(yōu)化算法 (Bi-Level Integrated System Synthesis,BLISS)[23]。多級優(yōu)化算法一般由系統(tǒng)級協(xié)調(diào)處理各學(xué)科之間的耦合關(guān)系,而各個(gè)學(xué)科可控制本學(xué)科的局部設(shè)計(jì)變量來進(jìn)行本學(xué)科的優(yōu)化設(shè)計(jì),其中協(xié)同優(yōu)化算法由于結(jié)構(gòu)簡單、學(xué)科自治性高及學(xué)科間的數(shù)據(jù)傳輸量小而備受關(guān)注。

    對于水面艦船航行性能領(lǐng)域內(nèi)的MDO算法應(yīng)用研究建議在掌握上述主要多級優(yōu)化算法(及其改進(jìn)形式)基本理論和基本方法的基礎(chǔ)上,借鑒航空、航天等領(lǐng)域MDO應(yīng)用研究的經(jīng)驗(yàn)和取得的成果,先模仿后改進(jìn)的技術(shù)路線來進(jìn)行。

    (4) 近似技術(shù)

    MDO過程涉及的學(xué)科多,學(xué)科分析的計(jì)算量大;對學(xué)科間交互作用的協(xié)調(diào)使得MDO過程的計(jì)算量急劇增加;另外,由于學(xué)科分析可能存在數(shù)值噪聲,若不對這些學(xué)科分析進(jìn)行必要的處理,將嚴(yán)重影響MDO過程的速度、質(zhì)量,降低工程應(yīng)用價(jià)值。近似技術(shù)作為學(xué)科分析工具和優(yōu)化算法之間的接口,可以大幅降低MDO過程的計(jì)算量,并具有平滑數(shù)值噪聲的作用,平衡MDO過程的精度和時(shí)間費(fèi)用。

    近似技術(shù)的本質(zhì)是通過構(gòu)造近似函數(shù),將復(fù)雜的學(xué)科分析從優(yōu)化過程中分離出來,而將便于計(jì)算的近似函數(shù)耦合到優(yōu)化過程中,進(jìn)行系列優(yōu)化,多次迭代循環(huán)后得到實(shí)際問題的近似最優(yōu)解,其主要難點(diǎn)在于對精度和時(shí)間的權(quán)衡,以及對近似誤差的修正。

    當(dāng)前,近似技術(shù)的研究已相當(dāng)成熟,工程系統(tǒng)中常用的近似方法有:響應(yīng)面模型(Response Surface Method)、變逼真度模型(Variable-fidelity Model)、Kriging 模型和 RBF 模型(Radical Basis Function)等,在本系統(tǒng)中可根據(jù)學(xué)科分析模型及工程精度需求,通過不同近似方法的應(yīng)用比較,選擇合適的方法并發(fā)展相適應(yīng)的誤差修正方法。

    (5)靈敏度分析技術(shù)

    MDO系統(tǒng)中的靈敏度分析是指對系統(tǒng)性能因設(shè)計(jì)變量或參數(shù)的變化而顯示出來的敏感程度的分析,其目的在于確定系統(tǒng)設(shè)計(jì)變量或參數(shù)對目標(biāo)函數(shù)或約束函數(shù)影響的大小,以及確定各子系統(tǒng)之間的耦合強(qiáng)度,最終用于指導(dǎo)設(shè)計(jì)與搜索方向、輔助決策。靈敏度分析技術(shù)是解決MDO問題的關(guān)鍵技術(shù)之一,其主要難點(diǎn)在于MDO過程中學(xué)科之間的函數(shù)關(guān)系模型復(fù)雜、維數(shù)多,且往往不能寫出顯式的函數(shù)表達(dá)式,使得靈敏度計(jì)算量大,效率低、求解困難。

    靈敏度分析又可分為學(xué)科靈敏度分析和多學(xué)科靈敏度分析,前者是后者的基礎(chǔ)。單學(xué)科靈敏度計(jì)算方法的發(fā)展已較為成熟,有符號微分法、有限差分法、復(fù)變量法、自動(dòng)微分法和解析方法等等;多學(xué)科靈敏度分析原則上可以采用學(xué)科靈敏度分析中已經(jīng)成熟的技術(shù),但在實(shí)際中,由于系統(tǒng)分析的維數(shù)太多,多學(xué)科靈敏度分析所需要的數(shù)據(jù)遠(yuǎn)比學(xué)科分析的更復(fù)雜,將學(xué)科靈敏度分析技術(shù)簡單的擴(kuò)展到多學(xué)科靈敏度分析并不現(xiàn)實(shí)。多學(xué)科靈敏度分析方法的應(yīng)用與MDO系統(tǒng)的分解策略,即MDO系統(tǒng)的類型相關(guān),如最優(yōu)靈敏度分析[26](Optimization Sensitivity Analysis,OSA)可用于層次系統(tǒng),全局靈敏度方程[27](Global Sensitivity Equation,GSE)和滯后耦合伴隨[28](Lagged-Coupled Adjoint,LCA)的靈敏度分析方法可用于耦合系統(tǒng)靈敏度分析。

    對于水面艦船航行性能MDO系統(tǒng)而言,因?qū)W科之間存在的耦合關(guān)系,宜選擇適合于耦合系統(tǒng)的靈敏度分析方法,如全局靈敏度或滯后耦合伴隨的靈敏度分析方法。

    (6)航行性能分析預(yù)報(bào)技術(shù);

    性能分析預(yù)報(bào)技術(shù)是進(jìn)行優(yōu)化設(shè)計(jì)的基礎(chǔ),性能分析預(yù)報(bào)的精度、可靠度直接影響著設(shè)計(jì)的質(zhì)量和水平。

    按照性能分析所依據(jù)的模型,可將性能分析方法分為兩大類:第一類是基于船型參數(shù)(如長寬比、寬度吃水比、方形系數(shù)和浮心坐標(biāo)等)的經(jīng)驗(yàn)方法,這類方法基于長期統(tǒng)計(jì)結(jié)果,簡單、快捷和穩(wěn)定,但適用范圍有限、局限性強(qiáng);第二類是基于船體幾何模型的數(shù)值方法,這類方法具有精細(xì)、客觀和適用范圍廣的特點(diǎn),并且能夠體現(xiàn)船型局部(或附體)變化對性能的影響,但這類方法建模困難、計(jì)算時(shí)間長,目前只對少數(shù)幾個(gè)性能指標(biāo)(如阻力)有一定的精度。這兩種方法互為補(bǔ)充,應(yīng)同時(shí)發(fā)展,其中第一類方法應(yīng)根據(jù)船型技術(shù)的發(fā)展,不斷修正、完善已有方法,并發(fā)展新的方法;第二類方法重點(diǎn)是提高分析的效率和可靠性,另一方面要不斷拓展其分析能力。

    值得注意的是,當(dāng)前快速性、耐波性和操縱性三個(gè)學(xué)科領(lǐng)域的性能分析預(yù)報(bào)技術(shù)的發(fā)展水平不一,其中快速性性能預(yù)報(bào)技術(shù)發(fā)展較快,而操縱性的發(fā)展卻相對滯后,建議采取邊研究邊應(yīng)用的策略,首先將成熟的、低精度的分析/預(yù)報(bào)技術(shù)納入系統(tǒng),并逐步由高效、高精度分析工具來替代。

    (7)船體幾何參數(shù)化表達(dá)及重構(gòu)技術(shù);

    船體幾何參數(shù)化表達(dá)及幾何重構(gòu)是航行性能設(shè)計(jì)優(yōu)化技術(shù)的共性技術(shù),它是聯(lián)系優(yōu)化算法與性能分析/預(yù)報(bào)模塊的橋梁,是航行性能優(yōu)化過程的關(guān)鍵環(huán)節(jié)。其難點(diǎn)在于對船體幾何曲面的準(zhǔn)確表達(dá),以及如何應(yīng)用盡可能少的參數(shù)來控制整體及局部的變形,以獲得盡可能大范圍的可行域,并保證船體曲面的光順性。

    與性能預(yù)報(bào)技術(shù)相適應(yīng),船體幾何參數(shù)化表達(dá)及重構(gòu)方法應(yīng)向兩個(gè)方向發(fā)展:一種是以若干船型參數(shù)為控制變量,實(shí)現(xiàn)對船體線型的幾何變換,如經(jīng)典的Lackenby變換方法、吃水函數(shù)法等;一種是借助于船體線型設(shè)計(jì)軟件或CAD軟件,以三維曲面或三向投影的二維曲線來直接表達(dá)船體幾何形狀,并用一組空間坐標(biāo)來控制船體幾何形狀的變形 (全船或局部),如疊加調(diào)和方法(Morphing Approach)、Bezier補(bǔ)丁方法(Bezier Patch Approach)、自由變形方法(Free-Form Deformation approach)和基于CAD變換方法等。

    (8)水面艦船航行性能MDO計(jì)算環(huán)境架構(gòu)

    建立便于MDO算法實(shí)現(xiàn)、提供適用的優(yōu)化算法庫、支持不同的水面艦船線型設(shè)計(jì)軟件、工具及性能分析工具的MDO計(jì)算框架是構(gòu)建水面艦船航行性能MDO集成設(shè)計(jì)平臺的關(guān)鍵技術(shù)之一。

    MDO計(jì)算環(huán)境不是功能模塊的簡單集成,而是旨在提供一個(gè)能夠描述和構(gòu)造MDO問題的通用平臺,它應(yīng)該對MDO技術(shù)的重要組成部分如近似模型、靈敏度分析數(shù)據(jù)的存儲、并行計(jì)算和性能分析數(shù)據(jù)的存儲等有良好的支撐能力。

    在當(dāng)前沒有太多經(jīng)驗(yàn)基礎(chǔ)的情況下,可參考成熟的MDO計(jì)算框架軟件,如iSIGHT、ModelCenter;還可以參考國內(nèi)艦船水動(dòng)力性能設(shè)計(jì)集成平臺開發(fā)經(jīng)驗(yàn),如中國船舶科學(xué)研究中心開發(fā)的SHIDS艦船綜合性能設(shè)計(jì)優(yōu)化系統(tǒng)、“先進(jìn)螺旋槳CAD軟件”等。

    5 結(jié) 語

    水面艦船航行性能MDO系統(tǒng)是一個(gè)極其復(fù)雜的系統(tǒng),不僅涉及到MDO技術(shù)、水面艦船水動(dòng)力性能分析預(yù)報(bào)技術(shù),還涉及到CAD技術(shù)、軟件工程技術(shù)和優(yōu)化技術(shù)等技術(shù)領(lǐng)域,其發(fā)展及應(yīng)用需要長期不懈的努力,在實(shí)際研究過程中可以采取邊開發(fā)邊應(yīng)用的模式。盡管在此MDO系統(tǒng)所應(yīng)用的性能分析工具、搜索策略等與傳統(tǒng)的優(yōu)化系統(tǒng)基本相同,但由于引入基于系統(tǒng)工程概念的MDO技術(shù),整個(gè)設(shè)計(jì)優(yōu)化過程在MDO算法的系統(tǒng)協(xié)調(diào)下開展,可以預(yù)見,綜合航行性能的提升幅度將是驚人的。

    航行性能設(shè)計(jì)與艦船總體設(shè)計(jì)技術(shù)是密切相關(guān)的,若視總體為一大系統(tǒng),則航行性能設(shè)計(jì)僅是大系統(tǒng)中的一個(gè)分系統(tǒng),故航行性能設(shè)計(jì)必須服從總體設(shè)計(jì)約束,滿足總體設(shè)計(jì)的需求,因此航行性能設(shè)計(jì)技術(shù)必須與總體設(shè)計(jì)技術(shù)統(tǒng)一考慮,所開發(fā)的系統(tǒng)才具有工程實(shí)用價(jià)值。因此,水面艦船航行性能MDO系統(tǒng)及集成平臺開發(fā),必須與艦船總體設(shè)計(jì)統(tǒng)一考慮,設(shè)計(jì)目標(biāo)、約束條件應(yīng)滿足總體設(shè)計(jì)要求。

    水面艦船航行性能MDO模型及集成平臺的研發(fā),將推動(dòng)航行性能設(shè)計(jì)模式的革命性轉(zhuǎn)變,將帶來綜合航行性能的大幅提升,推動(dòng)水面艦船水動(dòng)力構(gòu)型能力的跨越式發(fā)展。

    [1]Sobieszczanski-Sobieski J.A linear decomposition method for optimization problems-blueprint for development[R].NASA Technical Memorandum 83248,1982.

    [2]Sobieszczanski-Sobieski J.Optimization by decomposition:A step from hierarchic to non-hierachic systemps[C]//.2nd NASA/Air Force Symposium on Recent Advances in Multidisciplinary Analysis and Optimization,Hampton,VA,Sept.1988.NASA-TM-101494,NASA-CP-3031,Part1,1988.

    [3]Sobieszczanski-Sobieski J.The sensitivity of complex,internally coupled systems[J].AIAA Journal,1990,28(1):153-160.

    [4]王振東,陳小前,羅文彩,張為華.飛行器多學(xué)科設(shè)計(jì)優(yōu)化理論與應(yīng)用研究[M].北京:國防工業(yè)出版社,2006.

    [5]Christopher Gregory Hart.Multidisciplinary design optimization of complex engineering system for cost assessment under uncertainty[D].USA:Naval Architecture and Marine Engineering in The University of Michiga,2010.

    [6]Egorov I N,Kretinin G V,Leshchenko I A,et al.IOSO optimization toolkit-novel software to create better design[C].Proceedings of 8th AIAA/USAF/NASA/ISSMO Symposium on Multidisciplinary Analysis and Optimization,2002.

    [7]趙 敏,崔維成.多學(xué)科設(shè)計(jì)優(yōu)化研究應(yīng)用現(xiàn)狀綜述[J].中國造船,2007,48(3):63-72.

    [8]李世海.多學(xué)科設(shè)計(jì)優(yōu)化(MDO)算法研究[D].北京:中國地質(zhì)大學(xué),2009.

    [9]蘇子健.多學(xué)科設(shè)計(jì)優(yōu)化的分解、協(xié)調(diào)及不確定性研究[D].武漢:華中科技大學(xué),2008.

    [10]Gou Peng,Liu Wei,Cui Weicheng.A comparison of approximation methods for Multidisciplinary Design Optimization of ship structures[J].Journal of Ship Mechanics,2007,11(6):913-923.

    [11]陳小前.飛行器總體優(yōu)化設(shè)計(jì)理論與應(yīng)用研究[D].長沙:國防科技大學(xué),2001.

    [12]陳琪鋒.飛行器分布式協(xié)同進(jìn)化多學(xué)科設(shè)計(jì)優(yōu)化方法研究[D].長沙:國防科技大學(xué),2003.

    [13]Liu Wei,Cui Weicheng.Multidisciplinary Design Optimization(MDO):A promising tool for the design of HOV[J].Journal of Ship Mechanics,2004,8(6):95-112.

    [14]劉 蔚,操安喜,茍 鵬,崔維成.基于BLH框架的大深度載人潛水器總體性能的多學(xué)科設(shè)計(jì)優(yōu)化[J].船舶力學(xué),2008,12(1):110-117.

    [15]趙 敏,崔維成.BLISCO方法在載人潛水器設(shè)計(jì)中的應(yīng)用[J].船舶力學(xué),2009,13(2):259-268.

    [16]Pan Binbin,Cui Weicheng.Multidisciplinary Design Optimization methods for ship design[J].Journal of Ship Mechanics,2008,12(6):914-930.

    [17]Pan Binbin,Cui Weicheng,Leng Wenhao.Multidisciplinary Design Optimization of surface vessels[J].Journal of Ship Mechanics,2009.6.13(3):378-387.

    [18]Jiang Zhe,Cui Weicheng,Huang Xiaoping.Preliminary study of Multidisciplinary Design Optimization in truss spar concept design[J].Journal of Ship Mechanics,2009,13(3):444-457.

    [19]姜 哲,崔維成,黃曉平.基于響應(yīng)面的可變復(fù)雜度方法在桁架式Spar平臺方案設(shè)計(jì)中的應(yīng)用[J].船舶力學(xué),2010,14(7):771-781.

    [20]趙加鵬,余鋒等.魚雷外形多學(xué)科設(shè)計(jì)優(yōu)化方法研究[J].魚雷技術(shù),2009,17(5):10-14.

    [21]馮佰威,劉祖源,常海超.艦船多學(xué)科綜合優(yōu)化設(shè)計(jì)計(jì)算環(huán)境研究[J].艦船,2009(2):5-8.

    [22]Giesing Joseph P,Barthelemy J M.A summary of industry MDO application and needs[C].An AIAA White Paper,7th AIAA/USAF/NASA/ISSMO Symposium on Multidisciplinary Analysis and Optimization,1998.

    [23]Balling R J,Wilkison C A.Execution of multidisciplinary design optimization approaches on common test problems[J].AIAA Journal,1997,35(1):178-186.

    [24]Kroo J,Altus S,Braun R,et al.Multidisciplinary optimization methods for aircraft preliminary design[J].AIAA-94-4325,1994.

    [25]Sobieszczanski-Sobieski J.Application of global sensitivity equations in multidisciplinary aircraft synthesis[J].Aircraft,1990,27(12):1002-1010.

    [26]Sobieszczanski-Sobieski J,Barthelemy J F,Riley K M.Sensitivity of optimum solutions to problem parameters[C]//.Proceedings of the AIAA/ASME/ASCE/AHS 22nd Structures,Structural Dynamics and Materials Conference.Atlanta,Ga.,April 1981:184-205.

    [27]Sobieszczanski-Sobieski J.Sensitivity of complex,internally coupled systems[J].AIAA Journal,1990,28(1):153-162.

    [28]Joaquim R R,Martins A.A coupled-adjoint method for high-fidelity aero-structural optimization[D].USA:Stanford University,2002.

    Top-down design research for the navigational performance MDO system of naval surface combatant

    ZHAO Feng1,CHENG Su-bin2,YANG Lei1,WEI Xi-zhong1,LI Sheng-zhong1

    (1 China Ship Scientific Research Center,Wuxi 214082,China;2 Tongji University,Shanghai 200092,China)

    Navigational performance design of naval surface combatant is a typical multidisciplinary Design Optimization (MDO)problem.Therefore,it is suitable to adopt the innovational design idea and method which is based on system engineering idea to do,it no other than MDO method.This paper analyses the important meaning of the MDO applied research in the field of navigational performance design of the naval surface combatant,and introduces the concept and research content of the MDO method.The top-down design scheme of the navigational performance design system of the naval surface combatant is presented.The system design objection,evaluating indicators as well as the main content were all demonstrated.Finally,many crucial technologies which should be solved firstly were analyzed,and technical approaches were also suggested.

    naval surface combatant;navigational performance;MDO;top-down design

    U662

    A

    1007-7294(2012)11-1257-10

    2012-06-01

    趙 峰(1964-),男,中國船舶科學(xué)研究中心研究員;程素斌(1976-),女,同濟(jì)大學(xué)講師。

    猜你喜歡
    艦船航行水面
    艦船通信中的噪聲消除研究
    到慧骃國的航行
    艦船測風(fēng)傳感器安裝位置數(shù)值仿真
    水黽是怎樣浮在水面的
    創(chuàng)造足以亂真的水面反光
    小舟在河上航行
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    爭奪水面光伏
    能源(2016年3期)2016-12-01 05:11:02
    艦船腐蝕預(yù)防與控制系統(tǒng)工程
    一塊水面
    滇池(2014年5期)2014-05-29 07:30:25
    欧美一级毛片孕妇| 免费观看av网站的网址| 99国产综合亚洲精品| 久久国产精品男人的天堂亚洲| 国产一区二区在线观看av| 亚洲熟女毛片儿| 老司机午夜福利在线观看视频 | 国产精品98久久久久久宅男小说| 男女之事视频高清在线观看| 国产精品久久久久久精品古装| 午夜91福利影院| av超薄肉色丝袜交足视频| 麻豆乱淫一区二区| 日本wwww免费看| 国产免费现黄频在线看| 国产人伦9x9x在线观看| 国产91精品成人一区二区三区 | 后天国语完整版免费观看| 日本撒尿小便嘘嘘汇集6| 女性被躁到高潮视频| 丰满饥渴人妻一区二区三| 日韩成人在线观看一区二区三区| 精品国产一区二区三区久久久樱花| 1024视频免费在线观看| 亚洲成国产人片在线观看| 久久影院123| 99九九在线精品视频| 国产视频一区二区在线看| 美女视频免费永久观看网站| 亚洲精品久久成人aⅴ小说| 99re6热这里在线精品视频| 亚洲性夜色夜夜综合| 久久久久久久久久久久大奶| 国产精品欧美亚洲77777| 亚洲欧美日韩高清在线视频 | 最近最新中文字幕大全免费视频| 亚洲黑人精品在线| 桃花免费在线播放| 啦啦啦中文免费视频观看日本| 黄色怎么调成土黄色| av电影中文网址| 欧美成人午夜精品| 男女无遮挡免费网站观看| 日本黄色日本黄色录像| 成人亚洲精品一区在线观看| 老熟妇乱子伦视频在线观看| 91字幕亚洲| 青草久久国产| 免费看a级黄色片| 久久久久久亚洲精品国产蜜桃av| 国产一区二区三区在线臀色熟女 | 久久狼人影院| 成年人午夜在线观看视频| 最新在线观看一区二区三区| 国产三级黄色录像| 久久久水蜜桃国产精品网| 亚洲av第一区精品v没综合| 国产av又大| 亚洲国产欧美在线一区| 王馨瑶露胸无遮挡在线观看| 一二三四在线观看免费中文在| av免费在线观看网站| 久久久久久久国产电影| 深夜精品福利| 男女之事视频高清在线观看| 多毛熟女@视频| 搡老熟女国产l中国老女人| 中文亚洲av片在线观看爽 | 日本黄色日本黄色录像| 十八禁人妻一区二区| 国产成人免费观看mmmm| 国产精品电影一区二区三区 | 建设人人有责人人尽责人人享有的| 国产精品影院久久| 国产在线免费精品| 制服诱惑二区| 国产在线精品亚洲第一网站| 每晚都被弄得嗷嗷叫到高潮| 一本一本久久a久久精品综合妖精| 国产一区二区激情短视频| 一区在线观看完整版| 欧美成人免费av一区二区三区 | 国产男女内射视频| 丰满少妇做爰视频| 亚洲国产看品久久| 黄频高清免费视频| 国产97色在线日韩免费| 9191精品国产免费久久| 少妇粗大呻吟视频| videosex国产| 亚洲人成伊人成综合网2020| 精品少妇内射三级| 午夜视频精品福利| 肉色欧美久久久久久久蜜桃| 两人在一起打扑克的视频| 中文字幕精品免费在线观看视频| 久久精品国产a三级三级三级| 天堂动漫精品| 精品熟女少妇八av免费久了| 国产亚洲精品一区二区www | 久久精品aⅴ一区二区三区四区| 成人免费观看视频高清| 18禁黄网站禁片午夜丰满| 99国产精品99久久久久| 成人av一区二区三区在线看| 国产精品九九99| 欧美日韩成人在线一区二区| 黄色a级毛片大全视频| 色婷婷久久久亚洲欧美| 新久久久久国产一级毛片| 在线天堂中文资源库| 免费人妻精品一区二区三区视频| 国产日韩一区二区三区精品不卡| 建设人人有责人人尽责人人享有的| 国产一区二区三区综合在线观看| 国产欧美日韩综合在线一区二区| 久久99热这里只频精品6学生| 男人舔女人的私密视频| 久久久久久久精品吃奶| 亚洲欧洲精品一区二区精品久久久| 亚洲熟女毛片儿| 日韩三级视频一区二区三区| 亚洲精品自拍成人| 精品卡一卡二卡四卡免费| 天天躁狠狠躁夜夜躁狠狠躁| 久久九九热精品免费| 日本av手机在线免费观看| 精品一区二区三区视频在线观看免费 | 国产成人影院久久av| 国产精品久久电影中文字幕 | 亚洲精品乱久久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 久久ye,这里只有精品| 久久久久久亚洲精品国产蜜桃av| 国产97色在线日韩免费| 日韩人妻精品一区2区三区| 久久亚洲真实| 夜夜爽天天搞| 欧美日韩亚洲综合一区二区三区_| 交换朋友夫妻互换小说| 久久这里只有精品19| 久久午夜亚洲精品久久| 日韩一卡2卡3卡4卡2021年| 成人黄色视频免费在线看| 日本av免费视频播放| 熟女少妇亚洲综合色aaa.| 亚洲色图综合在线观看| 成人特级黄色片久久久久久久 | 国产成+人综合+亚洲专区| 动漫黄色视频在线观看| av线在线观看网站| 成年女人毛片免费观看观看9 | 又紧又爽又黄一区二区| 两人在一起打扑克的视频| 午夜成年电影在线免费观看| 色综合婷婷激情| 水蜜桃什么品种好| 亚洲 国产 在线| 丁香欧美五月| 人妻 亚洲 视频| 19禁男女啪啪无遮挡网站| 女人久久www免费人成看片| 大香蕉久久成人网| 午夜久久久在线观看| 欧美日韩国产mv在线观看视频| 国产精品久久久久成人av| 侵犯人妻中文字幕一二三四区| 国产在线观看jvid| 亚洲av美国av| 丝袜美腿诱惑在线| 午夜精品久久久久久毛片777| 亚洲天堂av无毛| 国产亚洲欧美精品永久| 一本久久精品| 亚洲精品在线美女| 国产精品1区2区在线观看. | a级片在线免费高清观看视频| 国产黄频视频在线观看| 法律面前人人平等表现在哪些方面| 欧美久久黑人一区二区| 国产aⅴ精品一区二区三区波| 91精品国产国语对白视频| 真人做人爱边吃奶动态| 欧美日韩av久久| 美女视频免费永久观看网站| 精品少妇一区二区三区视频日本电影| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕在线视频| 纯流量卡能插随身wifi吗| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品美女久久久久99蜜臀| 欧美黑人欧美精品刺激| 久久精品国产亚洲av高清一级| 国产精品偷伦视频观看了| 精品国产乱子伦一区二区三区| 乱人伦中国视频| 黑丝袜美女国产一区| 日韩大片免费观看网站| 少妇 在线观看| 久久久欧美国产精品| 狠狠狠狠99中文字幕| 熟女少妇亚洲综合色aaa.| 黄色视频在线播放观看不卡| 99re6热这里在线精品视频| 波多野结衣一区麻豆| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 中文字幕人妻丝袜制服| 国产高清国产精品国产三级| 亚洲五月色婷婷综合| 亚洲第一青青草原| 国产精品美女特级片免费视频播放器 | 如日韩欧美国产精品一区二区三区| 亚洲精品国产一区二区精华液| av免费在线观看网站| 91麻豆精品激情在线观看国产 | 午夜福利乱码中文字幕| 青草久久国产| 成人永久免费在线观看视频 | 法律面前人人平等表现在哪些方面| 麻豆成人av在线观看| 国产熟女午夜一区二区三区| www.熟女人妻精品国产| 免费看a级黄色片| 51午夜福利影视在线观看| 91老司机精品| 午夜福利视频在线观看免费| 亚洲av日韩在线播放| 亚洲avbb在线观看| 欧美 日韩 精品 国产| kizo精华| 久久久精品国产亚洲av高清涩受| 日日爽夜夜爽网站| 免费在线观看完整版高清| 啪啪无遮挡十八禁网站| 免费久久久久久久精品成人欧美视频| 亚洲精品一卡2卡三卡4卡5卡| 18禁观看日本| 欧美老熟妇乱子伦牲交| 久久av网站| 一本久久精品| 一区二区三区激情视频| 99久久人妻综合| 成年动漫av网址| 日韩熟女老妇一区二区性免费视频| 交换朋友夫妻互换小说| 波多野结衣一区麻豆| 午夜成年电影在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 国产在视频线精品| 色综合婷婷激情| 久久久精品区二区三区| 久热爱精品视频在线9| 无遮挡黄片免费观看| 在线永久观看黄色视频| 99在线人妻在线中文字幕 | 黄色视频不卡| 国产精品一区二区在线观看99| 最近最新中文字幕大全免费视频| 欧美日韩一级在线毛片| 伦理电影免费视频| 国产成+人综合+亚洲专区| 久久狼人影院| 中亚洲国语对白在线视频| 两人在一起打扑克的视频| 怎么达到女性高潮| 99精品在免费线老司机午夜| 无限看片的www在线观看| 伊人久久大香线蕉亚洲五| 国产野战对白在线观看| 亚洲精品粉嫩美女一区| 国产色视频综合| 久久这里只有精品19| netflix在线观看网站| 乱人伦中国视频| 国产一区二区三区综合在线观看| 男女边摸边吃奶| 国产区一区二久久| 脱女人内裤的视频| a级毛片黄视频| 韩国精品一区二区三区| 午夜福利一区二区在线看| 黄色视频,在线免费观看| 国产成人影院久久av| 成人三级做爰电影| 国产免费视频播放在线视频| 欧美激情久久久久久爽电影 | 少妇 在线观看| 免费女性裸体啪啪无遮挡网站| 国产免费视频播放在线视频| 国产在线免费精品| 久久久久久久久免费视频了| 欧美精品高潮呻吟av久久| 欧美亚洲 丝袜 人妻 在线| 窝窝影院91人妻| tube8黄色片| 日韩人妻精品一区2区三区| 法律面前人人平等表现在哪些方面| 黄色 视频免费看| 大陆偷拍与自拍| 在线播放国产精品三级| 人人妻,人人澡人人爽秒播| 老鸭窝网址在线观看| 免费一级毛片在线播放高清视频 | 十八禁网站网址无遮挡| 欧美日韩亚洲国产一区二区在线观看 | 最黄视频免费看| 不卡av一区二区三区| 热99久久久久精品小说推荐| 欧美日韩成人在线一区二区| 在线观看舔阴道视频| 精品人妻在线不人妻| 久9热在线精品视频| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 亚洲avbb在线观看| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 丝袜美足系列| 国产成人欧美在线观看 | 丁香六月欧美| 日本欧美视频一区| 亚洲精品自拍成人| 午夜福利影视在线免费观看| 十八禁网站网址无遮挡| 午夜激情久久久久久久| 老熟妇仑乱视频hdxx| 大型av网站在线播放| 我的亚洲天堂| 日韩 欧美 亚洲 中文字幕| 亚洲一区二区三区欧美精品| 亚洲av成人一区二区三| 人妻 亚洲 视频| 日本撒尿小便嘘嘘汇集6| 亚洲熟女毛片儿| 免费观看人在逋| 亚洲精品久久午夜乱码| 亚洲成av片中文字幕在线观看| 精品午夜福利视频在线观看一区 | 国产成人系列免费观看| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 亚洲第一青青草原| 婷婷丁香在线五月| 亚洲全国av大片| 最近最新免费中文字幕在线| 国产成人精品久久二区二区免费| 国产精品国产av在线观看| 日本撒尿小便嘘嘘汇集6| 日韩中文字幕欧美一区二区| 男女下面插进去视频免费观看| 男女高潮啪啪啪动态图| av免费在线观看网站| 涩涩av久久男人的天堂| 1024香蕉在线观看| 中国美女看黄片| 国产成人免费观看mmmm| 欧美日韩国产mv在线观看视频| 黄色 视频免费看| 亚洲成av片中文字幕在线观看| 国产精品偷伦视频观看了| 91大片在线观看| 久久久久久人人人人人| 国产男女超爽视频在线观看| 老司机午夜福利在线观看视频 | avwww免费| 一进一出抽搐动态| 久久精品国产综合久久久| 国产精品国产av在线观看| 免费人妻精品一区二区三区视频| 久久99热这里只频精品6学生| 欧美+亚洲+日韩+国产| 久久中文字幕人妻熟女| 在线观看免费高清a一片| 欧美在线黄色| 免费久久久久久久精品成人欧美视频| 丝袜喷水一区| 国产麻豆69| 在线观看免费视频日本深夜| 青草久久国产| 丝袜美足系列| 国产深夜福利视频在线观看| 一级a爱视频在线免费观看| 热99久久久久精品小说推荐| 日韩一卡2卡3卡4卡2021年| 变态另类成人亚洲欧美熟女 | 免费在线观看完整版高清| 欧美一级毛片孕妇| 热re99久久精品国产66热6| 亚洲中文av在线| 国产成人欧美在线观看 | 成人精品一区二区免费| 国产一区二区三区视频了| 视频区欧美日本亚洲| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 国产精品麻豆人妻色哟哟久久| 男女之事视频高清在线观看| 一级片'在线观看视频| 午夜91福利影院| 成年动漫av网址| 国产精品偷伦视频观看了| xxxhd国产人妻xxx| av天堂在线播放| 亚洲av日韩在线播放| 搡老乐熟女国产| 欧美日韩亚洲高清精品| 人人妻人人澡人人爽人人夜夜| 后天国语完整版免费观看| av线在线观看网站| 亚洲精品美女久久久久99蜜臀| 男女高潮啪啪啪动态图| avwww免费| 一本色道久久久久久精品综合| 欧美 日韩 精品 国产| 这个男人来自地球电影免费观看| 老司机深夜福利视频在线观看| 久久中文字幕人妻熟女| netflix在线观看网站| 老汉色av国产亚洲站长工具| 中文字幕最新亚洲高清| 99re在线观看精品视频| 精品欧美一区二区三区在线| 极品少妇高潮喷水抽搐| 亚洲人成电影免费在线| 丁香欧美五月| bbb黄色大片| 国产无遮挡羞羞视频在线观看| 色精品久久人妻99蜜桃| 国产精品一区二区精品视频观看| 亚洲国产精品一区二区三区在线| 99国产精品99久久久久| 91精品三级在线观看| 亚洲精品久久午夜乱码| 深夜精品福利| 久久国产精品人妻蜜桃| 亚洲成人国产一区在线观看| 久久久久精品国产欧美久久久| 另类亚洲欧美激情| 久久久精品94久久精品| 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看 | 丝袜人妻中文字幕| 黄色成人免费大全| 亚洲av欧美aⅴ国产| 欧美精品啪啪一区二区三区| 精品免费久久久久久久清纯 | 俄罗斯特黄特色一大片| 亚洲综合色网址| 人妻 亚洲 视频| 亚洲avbb在线观看| 亚洲国产看品久久| 亚洲avbb在线观看| 国产高清激情床上av| av片东京热男人的天堂| 老司机靠b影院| 成人手机av| 日本黄色视频三级网站网址 | 欧美日韩亚洲综合一区二区三区_| 国精品久久久久久国模美| 亚洲第一青青草原| 最近最新中文字幕大全电影3 | 91九色精品人成在线观看| videosex国产| 亚洲午夜理论影院| 日韩视频一区二区在线观看| 脱女人内裤的视频| 免费观看人在逋| 老司机深夜福利视频在线观看| 丝袜美腿诱惑在线| 黄色a级毛片大全视频| 高清视频免费观看一区二区| 午夜免费鲁丝| 美女高潮喷水抽搐中文字幕| 亚洲中文日韩欧美视频| 成人免费观看视频高清| 王馨瑶露胸无遮挡在线观看| 大陆偷拍与自拍| 热re99久久国产66热| 亚洲中文av在线| 99国产精品一区二区三区| 久久中文字幕人妻熟女| 一区二区三区乱码不卡18| 叶爱在线成人免费视频播放| 国产精品久久久久成人av| 精品国产乱码久久久久久男人| 黄色怎么调成土黄色| 人人妻人人爽人人添夜夜欢视频| 国产精品一区二区在线不卡| a级片在线免费高清观看视频| 制服人妻中文乱码| 国产91精品成人一区二区三区 | 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 蜜桃国产av成人99| 男人操女人黄网站| 中文亚洲av片在线观看爽 | videosex国产| 国产高清videossex| 日韩成人在线观看一区二区三区| 精品欧美一区二区三区在线| tube8黄色片| 少妇的丰满在线观看| 国产亚洲精品第一综合不卡| 成年人黄色毛片网站| 下体分泌物呈黄色| 成人精品一区二区免费| 少妇 在线观看| 亚洲人成77777在线视频| 国产高清videossex| 久久国产精品男人的天堂亚洲| 精品午夜福利视频在线观看一区 | 91成人精品电影| 亚洲人成伊人成综合网2020| videos熟女内射| 欧美黄色淫秽网站| 国产亚洲精品第一综合不卡| 性少妇av在线| 精品亚洲成a人片在线观看| 久久久久视频综合| 欧美黑人精品巨大| 多毛熟女@视频| 美女高潮到喷水免费观看| 一级毛片女人18水好多| 日本a在线网址| 97人妻天天添夜夜摸| 亚洲av国产av综合av卡| 妹子高潮喷水视频| 欧美成狂野欧美在线观看| 日本黄色视频三级网站网址 | 日韩一卡2卡3卡4卡2021年| 日韩成人在线观看一区二区三区| 亚洲情色 制服丝袜| 成年人免费黄色播放视频| 午夜日韩欧美国产| 国产欧美日韩精品亚洲av| 不卡av一区二区三区| av超薄肉色丝袜交足视频| 成人特级黄色片久久久久久久 | 国产成人精品久久二区二区免费| 俄罗斯特黄特色一大片| 在线观看舔阴道视频| 欧美日韩成人在线一区二区| 飞空精品影院首页| 欧美另类亚洲清纯唯美| 日韩一卡2卡3卡4卡2021年| 国产免费视频播放在线视频| 成人免费观看视频高清| 中文字幕精品免费在线观看视频| 欧美在线一区亚洲| 一级黄色大片毛片| 国产黄色免费在线视频| 天堂俺去俺来也www色官网| 中文欧美无线码| 欧美亚洲日本最大视频资源| 精品熟女少妇八av免费久了| 亚洲av成人不卡在线观看播放网| 无人区码免费观看不卡 | 午夜精品久久久久久毛片777| 久久久久久久国产电影| 色婷婷av一区二区三区视频| 国产在线观看jvid| 欧美另类亚洲清纯唯美| 国产精品国产av在线观看| 青青草视频在线视频观看| 欧美激情极品国产一区二区三区| 亚洲午夜理论影院| 亚洲av日韩在线播放| 变态另类成人亚洲欧美熟女 | 丰满迷人的少妇在线观看| 黄色视频不卡| 午夜福利影视在线免费观看| 欧美日韩视频精品一区| 在线观看免费视频网站a站| 真人做人爱边吃奶动态| 最新的欧美精品一区二区| 精品亚洲成a人片在线观看| 我的亚洲天堂| av网站免费在线观看视频| 丝瓜视频免费看黄片| 丰满迷人的少妇在线观看| 亚洲欧美一区二区三区黑人| 老熟妇仑乱视频hdxx| 色尼玛亚洲综合影院| 欧美日韩亚洲综合一区二区三区_| 亚洲精品中文字幕在线视频| 夜夜夜夜夜久久久久| 国产欧美日韩精品亚洲av| 午夜激情av网站| 国产成人免费观看mmmm| 久久久久精品人妻al黑| 九色亚洲精品在线播放| 国产区一区二久久| 欧美日本中文国产一区发布| 国产精品98久久久久久宅男小说| e午夜精品久久久久久久| 久久久久精品人妻al黑| 国产精品免费视频内射| e午夜精品久久久久久久| 老汉色av国产亚洲站长工具| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久亚洲精品国产蜜桃av| 精品人妻熟女毛片av久久网站| 9热在线视频观看99| 操美女的视频在线观看| 捣出白浆h1v1| 又紧又爽又黄一区二区| 老司机午夜福利在线观看视频 | 淫妇啪啪啪对白视频| 啪啪无遮挡十八禁网站| 丁香六月天网| 看免费av毛片| 丝袜人妻中文字幕| 国产真人三级小视频在线观看| 亚洲国产欧美网|