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

    基于自由變形的自主水下航行器殼體優(yōu)化設(shè)計方法研究

    2022-09-21 03:38:10張代雨王志東呂海寧楊建民
    船舶力學(xué) 2022年9期
    關(guān)鍵詞:型線殼體控制點

    張 貝,張代雨,王志東,呂海寧,楊建民

    (1.上海交通大學(xué)海洋工程國家重點實驗室,上海 200240;2.上海交通大學(xué)三亞崖州灣深??萍佳芯吭?,海南三亞 572000;3.江蘇科技大學(xué)船舶與海洋工程學(xué)院,江蘇鎮(zhèn)江 212003)

    0 引 言

    隨著人類社會的發(fā)展,陸地上的資源將被消耗殆盡,海洋在資源、環(huán)境空間等方面擁有得天獨厚的優(yōu)勢,已經(jīng)成為人類發(fā)展的重要基地,但在新世紀(jì)里如何合理開發(fā)、利用、保護(hù)海洋仍面臨巨大挑戰(zhàn)。海洋環(huán)境探測及資源開發(fā)利用,離不開先進(jìn)的海洋觀測技術(shù)及探測技術(shù),自主水下航行器(autonomous underwater vehicle,AUV)將人工智能、探測識別、信息融合、智能控制、系統(tǒng)集成等多方面技術(shù)集中應(yīng)用于同一水下載體,通過自主決策、控制完成復(fù)雜海洋環(huán)境中預(yù)定任務(wù),能夠高效率地開展各種作業(yè)及探測任務(wù),具有觀測范圍大、效率高、機(jī)動性能好等優(yōu)點[1-4],已在大范圍海洋探測和資源開發(fā)領(lǐng)域發(fā)揮出越來越重要的作用。

    AUV 殼體不僅為內(nèi)部設(shè)備提供空間、密封、承受外壓,還與外部流體相互作用,產(chǎn)生航行阻力、顛覆力矩及流體噪聲等[5-6],因此,AUV 殼體外形的優(yōu)劣對其整體性能具有重要影響。AUV 殼體一般要求體積豐滿系數(shù)越大越好,以提供更大的內(nèi)部容積來安裝探測設(shè)備;頭部的轉(zhuǎn)捩點越靠后越好,以降低湍流及轉(zhuǎn)捩區(qū)流動噪聲對制導(dǎo)系統(tǒng)的嚴(yán)重干擾;航行阻力越小越好,以降低整個航行過程的能量需求;最大減壓系數(shù)越小越好,避免航行過程中的空化現(xiàn)象等。為了更好地達(dá)到上述要求,目前,優(yōu)化設(shè)計方法被廣泛應(yīng)用于AUV的殼體設(shè)計中[7-10],且發(fā)揮出越來越重要的作用。

    Alvarez 等[11]使用Nystrom 型線對一型工作在近自由表面(浮潛深度)的觀測型AUV 殼體進(jìn)行幾何參數(shù)化,以殼體總體積為約束條件,波浪阻力為目標(biāo)函數(shù),采用模擬退火算法對AUV 殼體進(jìn)行優(yōu)化設(shè)計,最大減少了25%的航行阻力;Joung 等[12]采用Myring 型線對AUV 殼體進(jìn)行幾何參數(shù)化,用計算流體力學(xué)方法對AUV 總阻力進(jìn)行精確計算,并實現(xiàn)了網(wǎng)格的自動化快速生成,引入Ansys中的實驗設(shè)計優(yōu)化方法進(jìn)行了AUV 殼體的優(yōu)化設(shè)計,取得了滿意的優(yōu)化結(jié)果;Alam 等[13]基于NSGA-II(non-dominated sorting genetic algorithm)和IDEA(infeasibility driven evolutionary algorithm)兩種流行的優(yōu)化算法,提出了一種針對多種AUV 殼體外形的優(yōu)化框架,該框架使用Nystrom 型線對AUV 殼體進(jìn)行幾何參數(shù)化,并引入CATIA 和Matlab 軟件進(jìn)行AUV 殼體的三維建模和性能計算;Gao 等[6]采用Myring 型線對AUV 殼體進(jìn)行幾何參數(shù)化,利用ISIGHT、ICEM 和Fluent 等軟件建立了AUV 殼體的優(yōu)化設(shè)計流程,并對MIGA(multi-island genetic algorithm)和PSO(particle swarm optimization)兩種優(yōu)化方法的結(jié)果進(jìn)行了對比,發(fā)現(xiàn)MIGA 可以得到更好的殼體外形;Hou 等[14]、Hang 等[15]在AUV 殼體的優(yōu)化設(shè)計過程中,考慮了不確定參數(shù)對航行阻力的影響,進(jìn)行了AUV殼體的可靠性優(yōu)化設(shè)計,結(jié)果表明,設(shè)計的AUV殼體在海洋環(huán)境干擾的情況下具有更高的可靠性。此外,Chen 等[16-17]、Zhang 等[18]和Wang 等[19]還使用格蘭維爾型線對AUV 殼體進(jìn)行了幾何建模,并在此基礎(chǔ)上進(jìn)行了多學(xué)科優(yōu)化設(shè)計,均取得了很好的設(shè)計效果。

    綜上,雖然眾多學(xué)者針對AUV 殼體進(jìn)行了一系列優(yōu)化設(shè)計研究,并取得了豐碩的研究成果,但均采用傳統(tǒng)的Nystrom 型線、Myring型線、卡洛斯型線及格蘭維爾型線等進(jìn)行幾何參數(shù)化建模,頭部曲線段或尾部曲線段的形狀可調(diào)參數(shù)局限于1~2 個,導(dǎo)致優(yōu)化過程中控制頭部曲線段和尾部曲線段的優(yōu)化變量較少,未能充分探索AUV 殼體外形的可變空間,從而得到最優(yōu)的殼體外形。針對這一現(xiàn)狀,本文提出了一種基于自由變形(free-form deformation,F(xiàn)FD)的AUV 殼體優(yōu)化設(shè)計方法,在優(yōu)化過程中實現(xiàn)頭部曲線段和尾部曲線段的自由變形,充分挖掘AUV 殼體的最優(yōu)特性。自由變形優(yōu)化設(shè)計方法主要基于B樣條技術(shù),建立AUV殼體的自由變形控制體,并在此基礎(chǔ)上提出一種全局參數(shù)化方法對FFD控制體進(jìn)行間接變形操縱,在實現(xiàn)自由變形的前提下,大大減少參數(shù)化變量個數(shù),提高優(yōu)化效率。最后,為了驗證所提方法的有效性,分別使用本文方法和傳統(tǒng)方法對一型AUV 殼體進(jìn)行優(yōu)化設(shè)計,并對比分析了所得的優(yōu)化結(jié)果。

    1 基于傳統(tǒng)參數(shù)化的殼體優(yōu)化設(shè)計方法

    1.1 傳統(tǒng)的殼體參數(shù)化方法

    AUV殼體外形多選擇回轉(zhuǎn)體結(jié)構(gòu),原因在于回轉(zhuǎn)體便于加工,成本低;具有較好的流線型,流體動力性能優(yōu)良;在容積一定的情況下表面積最小,可減少航行阻力與殼體用料等。此外,回轉(zhuǎn)體由母線繞旋轉(zhuǎn)軸旋轉(zhuǎn)而成,故回轉(zhuǎn)體的幾何形狀可由母線的型線表示,描述簡單。對于AUV 殼體,依據(jù)其幾何特點,傳統(tǒng)參數(shù)化方法一般將其母線型線分為三部分:頭部曲線段AB、平行中段BC及尾部曲線段CD,如圖1所示。

    圖1 AUV殼體外形示意圖Fig.1 Profile of the AUV hull shape

    圖1中,Lh為頭部曲線段長度,Lc為平行中段長度,Lt為尾部曲線段長度,d為平行中段直徑,Rh(x)為頭部曲線段的型線,Rt(x)為尾部曲線段的型線。其中,Rh(x)和Rt(x)通常采用參數(shù)化型線描述,如Nystrom 型線、Myring 型線、卡洛斯型線、格蘭維爾型線等。下面詳細(xì)介紹Nystrom、Myring 和格蘭維爾三種常用的傳統(tǒng)殼體參數(shù)化方法。

    (1)Nystrom 型線[20-21],其由Nystrom 在1868 年提出,使用1/4 橢圓加一段拋物線來描述AUV 殼體的型線,又稱水滴型。頭部和尾部的型線具體表述為

    式中:nf和na是可調(diào)參數(shù),通過改變nf和na的值,可改變頭部和尾部曲線段的形狀;此外,nf和na有一定的可行域,通常取1<nf<3,1<na<3。

    (2)Myring型線[22-23],其頭部采用類Nystrom 型線,尾部為Myring提出的尖尾型線方程,著名的航行器Remus采用該型線,具體的型線表達(dá)式如式(2)~(3)所示。

    式中:n和θ為可調(diào)參數(shù),且θ表示尾部曲線末端與回轉(zhuǎn)體轉(zhuǎn)軸間的夾角;n的可行域與Nystrom類似,取1<n<3;θ影響尾部曲線段去流特性,取7°<θ<16°。

    (3)格蘭維爾型線[24-25],其主要利用多項式構(gòu)造型線,在型線端點、連接點、曲率等處添加約束條件,構(gòu)造無拐點、無極值方程。頭部和尾部的具體型線方程如式(4)~(5)所示:

    式中,qh1和qh2為頭部曲線段的可調(diào)參數(shù),qt1和qt2為尾部曲線段的可調(diào)參數(shù)。另外需注意的一點是,并不是所有的可調(diào)參數(shù)值均可構(gòu)造適用的型線,為了得到無拐點、無極值的型線,公式(4)和公式(5)的可調(diào)參數(shù)可行域取0<qh1<3,0<qh2<12,0<qt1<2,0<qt2<20。

    1.2 傳統(tǒng)的AUV殼體優(yōu)化流程

    對于AUV 殼體優(yōu)化設(shè)計,一般要求:(1)殼體型線具有較大的內(nèi)部轉(zhuǎn)載空間,具體可利用豐滿度系數(shù)ψ提出;(2)內(nèi)部能源搭載量有限,為了航行更遠(yuǎn)距離,航行阻力Cx越小越好;(3)頭部曲線段轉(zhuǎn)捩點的位置越靠后越好,以降低湍流即轉(zhuǎn)捩區(qū)流動噪聲對制導(dǎo)系統(tǒng)的嚴(yán)重干擾,但考慮到轉(zhuǎn)捩點計算復(fù)雜、難以準(zhǔn)確確定,直接作為設(shè)計參數(shù)是不實用的,從定性上講,轉(zhuǎn)捩點發(fā)生在逆壓區(qū),即最小壓力系數(shù)Cpmin的位置之后,故常使用最小壓力系數(shù)點位置xCpmin進(jìn)行替代;(4)頭部曲線段不產(chǎn)生空化,且不發(fā)生流動分離,其中,流動分離的產(chǎn)生,與逆壓梯度dCp/dx具有重要關(guān)系,一般來說,逆壓梯度越小,越不易發(fā)生流動分離,因此壓力分布越平坦越好。綜上,AUV殼體的總體優(yōu)化設(shè)計準(zhǔn)則是豐滿度系數(shù)越大越好,阻力越小越好,最大減壓系數(shù)位置越后越好,壓力分布曲線越平坦越好。

    對于傳統(tǒng)的AUV 殼體優(yōu)化設(shè)計流程,主要利用Nystrom、Myring 或格蘭維爾等型線對殼體幾何外形進(jìn)行參數(shù)化建模,取相應(yīng)的可調(diào)參數(shù)作為優(yōu)化設(shè)計變量,在優(yōu)化過程中改變AUV 殼體的外形;選取豐滿度系數(shù)ψ、航行阻力Cx、最小壓力系數(shù)點位置xCpmin、逆壓梯度dCp/dx等性能參數(shù)中的一個或多個作為目標(biāo)函數(shù),并且滿足根據(jù)任務(wù)指標(biāo)設(shè)定的約束條件gi(ψ,Cx,xCpmin,dCp/dx)≥0;使用計算流體力學(xué)(computational fluid dynamics,CFD)方法[26]進(jìn)行航行阻力、殼體表面壓力分布等流體動力參數(shù)的計算;采用序列二次規(guī)劃算法(sequential quadratic programming,SQP)[27]、粒子群算法(particle swarm optimization,PSO)[28]及NSGAII 算法(non-dominated sorting genetic algorithm)[29]等控制整個優(yōu)化設(shè)計的運(yùn)行,得到最優(yōu)的設(shè)計方案。詳細(xì)的AUV殼體外形優(yōu)化流程如圖2所示。

    圖2 傳統(tǒng)的AUV殼體外形優(yōu)化流程Fig.2 Traditional optimization process of AUV hull shape

    2 基于自由變形的殼體優(yōu)化設(shè)計方法

    對于AUV殼體的傳統(tǒng)優(yōu)化方法,雖然Nystrom型線、Myring型線、卡洛斯型線及格蘭維爾型線具有良好的流線性能,但頭部曲線段和尾部曲線段僅僅具有1 或2 個可調(diào)參數(shù),在一定程度上限制了AUV殼體的外形變化范圍,減少了優(yōu)化過程中的探索空間,不利于得到最優(yōu)的AUV殼體外形。FFD方法既可保持初始外形的流線性能,又具備優(yōu)良的局部變形能力[30-31],本文引入FFD 方法進(jìn)行殼體外形的參數(shù)化,提出了一種基于自由變形的AUV殼體優(yōu)化設(shè)計方法,最大程度地提高AUV的殼體性能。首先,基于B樣條方法,提出了AUV 殼體外形的FFD 控制體建立方法,實現(xiàn)對AUV 殼體外形的自由變形;然后,考慮到AUV 殼體的回轉(zhuǎn)體特性,基于軸變形方法,提出了一種全局參數(shù)化方法對FFD 控制體進(jìn)行間接變形操縱,大大減少了參數(shù)化變量個數(shù);最后,給出了基于FFD參數(shù)化的豐滿度計算方法,建立了基于自由變形方法的AUV殼體外形優(yōu)化設(shè)計框架。下面將分節(jié)進(jìn)行詳細(xì)闡述。

    2.1 基于B樣條的自由變形方法

    FFD 方法為Sederberg 和Parry[32]在1986 年提出,該方法主要借鑒彈性物體受力變形的思想,假設(shè)將研究對象放置于FFD 控制體內(nèi),給控制體施加變形后,則控制體內(nèi)部的空間會發(fā)生形變,進(jìn)而嵌入其中的研究對象也會根據(jù)某種映射關(guān)系發(fā)生變化。

    針對AUV 殼體優(yōu)化設(shè)計來說,既可以將具有良好流線性能的初始幾何外形嵌入到FFD 控制體中進(jìn)行自由變形,也可以將由初始幾何外形進(jìn)行網(wǎng)格劃分得到的表面網(wǎng)格嵌入到FFD 控制體中進(jìn)行自由變形。對于前一種嵌入方法,在優(yōu)化過程中,需要進(jìn)行AUV 殼體幾何重新生成、網(wǎng)格劃分、流體參數(shù)計算共三個步驟以得到所需的目標(biāo)函數(shù)和約束值;對于后一種嵌入方法,由于直接對表面網(wǎng)格進(jìn)行變形,在優(yōu)化過程中不需對AUV 殼體進(jìn)行幾何重新生成,節(jié)省了優(yōu)化過程步驟。又考慮到通過表面網(wǎng)格可以逆向生成2 階連續(xù)的幾何外形,本文選擇在FFD 控制體中嵌入AUV 殼體的表面網(wǎng)格以實現(xiàn)優(yōu)化過程中的AUV外形變化,如圖3所示。

    圖3 嵌入AUV殼體表面網(wǎng)格的FFD控制體Fig.3 FFD control volume embedded in the surface mesh of AUV hull

    圖中,紅色點為FFD 控制體的控制點,AUV 殼體的表面網(wǎng)格嵌入到FFD 控制體內(nèi),通過改變控制點的位置,移動FFD 控制點可間接實現(xiàn)AUV 殼體表面網(wǎng)格的自由變形控制。此處,還需指出的是對于AUV 的回轉(zhuǎn)殼體外形,為了便于控制變形,建立的FFD 控制體一般為橫截面為正方形的長方體,并且長方體與回轉(zhuǎn)殼體的中心軸共軸。

    下面詳細(xì)描述采用FFD 方法實現(xiàn)AUV 殼體表面網(wǎng)格自由變形的方法。首先,需要在表面網(wǎng)格周圍建立FFD 控制體。Bernstein 多項式、B 樣條函數(shù)等[33]均可作為基函數(shù)來創(chuàng)建FFD 控制體,但考慮到Bernstein多項式的局部變形性能不佳,采用B樣條函數(shù)創(chuàng)建FFD控制體。

    曲線節(jié)點矢量U上的p階B樣條基函數(shù)Ni,p(u)可表述如下:

    由B樣條基函數(shù)創(chuàng)建的控制體方程可表述成如下形式:

    式中:Pi,j,k為FFD 控制體上的控制點在笛卡爾坐標(biāo)系中的坐標(biāo)值;Ni,p(u)、Nj,p(v)和Nk,p(w)分別為參數(shù)化空間(u,v,w)上的p階B樣條基函數(shù);K為笛卡爾坐標(biāo)系中的坐標(biāo)值。

    然后,計算待變形的各個網(wǎng)格點坐標(biāo)值X在FFD 控制體中的局部坐標(biāo)(u,v,w),其中局部坐標(biāo)的三個分量均位于[0,1]之間,通過逆向求解公式(8)確定。需注意的是,每個網(wǎng)格點的局部坐標(biāo)只需計算一次,不需要在每次CFD網(wǎng)格變形的過程中重復(fù)計算。

    接著,改變FFD 控制點的位置,產(chǎn)生位移量ΔPi,j,k,則任一個局部坐標(biāo)為(u,v,w)的網(wǎng)格點的位移ΔX可由公式(9)確定:

    最后,可得到變形后表面網(wǎng)格上每一個網(wǎng)格點的位置X?(u,v,w)為

    2.2 基于軸變形的全局參數(shù)化方法

    盡管FFD參數(shù)化方法具有不依賴于外形拓?fù)浣Y(jié)構(gòu)、可保持初始外形連續(xù)光滑、變形簡單靈活等優(yōu)點,但若把FFD 控制體的每一個控制點均作為優(yōu)化變量,則會導(dǎo)致優(yōu)化設(shè)計變量規(guī)模過大,進(jìn)而造成優(yōu)化設(shè)計計算量大、效率低、收斂慢等問題。此外,對于AUV 殼體外形,其一般要求設(shè)計為回轉(zhuǎn)體結(jié)構(gòu),若采用每一個FFD控制點作為優(yōu)化變量,還需要添加大量的幾何約束以保持優(yōu)化后的殼體外形保持為回轉(zhuǎn)體結(jié)構(gòu)。因此,在使用FFD 方法進(jìn)行AUV 殼體變形時,如何既能充分利用FFD 參數(shù)化方法的優(yōu)點,又能合理控制優(yōu)化設(shè)計變量的規(guī)模,是一個亟需解決的問題。

    本文提出一種基于軸變形的全局FFD參數(shù)化方法,在FFD控制體內(nèi)引入一條參考軸,并將FFD的所有控制點關(guān)聯(lián)到該參考軸上,通過設(shè)置全局變量對參考軸進(jìn)行變形控制,進(jìn)而實現(xiàn)對FFD控制點的分組集中控制,減少優(yōu)化設(shè)計變量的規(guī)模。該方法詳細(xì)闡述如下:

    首先,使用如下的p階B 樣條曲線定義一條參考軸,且該軸為正方形橫截面的FFD 控制體的對稱中心軸,如圖4中的綠色軸線。

    圖4 FFD控制點到參考軸的映射圖Fig.4 Mapping diagram of FFD control points to the reference axis

    式中,Qi(i= 1,…,N)是曲線的N個控制點,其對應(yīng)于曲線上N個點的坐標(biāo)值;Ni,p(s)是曲線節(jié)點矢量S上的p階B樣條基函數(shù)。

    然后,將FFD 控制點映射到參考軸上。映射的主要思路為:對于每個FFD 控制點Pi,j,k,分別求出參考軸上與其距離最近的點,則該點為相應(yīng)FFD 控制點在參考軸上的映射點,如圖4 中的綠色點,并計算出其在參考軸上的參數(shù)化坐標(biāo)值為s?i,j,k。該映射方法保證每個FFD 控制點與相應(yīng)映射點的直線與參考軸在映射點處的切線具有垂直關(guān)系。求出映射點后,計算每個FFD 控制點到映射點的向量T→i,j,k=Pi,j,k-C(s?i,j,k) ,如圖4中的藍(lán)色虛線,以備后續(xù)步驟使用。

    接著,在參考軸的N個控制點Qi處,給定在這些點處的回轉(zhuǎn)半徑比例因子Ri,且基于公式(11)中定義的p階B 樣條曲線,對參考軸上每一參數(shù)化坐標(biāo)處的回轉(zhuǎn)半徑比例因子進(jìn)行擬合,具體公式如式(12)所示:

    式中,僅采用Ri(i= 1,…,N)代替原曲線的N個控制點,其他參數(shù)和原曲線相同,以保證在曲線控制點Qi的參數(shù)化坐標(biāo)si處,求得的回轉(zhuǎn)半徑比例因子R(si)與Ri相等。

    再者,根據(jù)每一個FFD 控制點映射在參考軸上的參數(shù)化坐標(biāo)值s?i,j,k,和公式(12)求出的對應(yīng)的回轉(zhuǎn)半徑比例因子值R(s?i,j,k),進(jìn)而通過式(13)計算形狀改變后的FFD控制點P?i,j,k:

    最后,根據(jù)FFD 控制點的新坐標(biāo)P?i,j,k,計算出位移量ΔPi,j,k,進(jìn)而通過公式(9)~(10)求得嵌入FFD控制體中的表面網(wǎng)格的新坐標(biāo)X?。

    綜上,在上述整個過程中,主要通過給定不同的回轉(zhuǎn)半徑比例因子Ri(i=1,…,N)控制FFD控制體的變形,進(jìn)而控制AUV 表面網(wǎng)格的變形。此處,回轉(zhuǎn)半徑比例因子Ri(i=1,…,N)即為全局參數(shù)化變量,其變量個數(shù)可根據(jù)優(yōu)化精度需要進(jìn)行給定。

    2.3 基于自由變形的AUV殼體優(yōu)化框架

    基于2.2 節(jié)提出的AUV 殼體全局參數(shù)化方法,本文提出了一種基于自由變形的AUV 殼體外形優(yōu)化設(shè)計框架。該框架選取公式(12)的回轉(zhuǎn)半徑比例因子Ri(i=1,…,N)作為優(yōu)化設(shè)計變量,優(yōu)化器通過改變Ri的值,實現(xiàn)對FFD 控制點的位移控制,進(jìn)而實現(xiàn)AUV 殼體表面網(wǎng)格的自由變形;選取豐滿度系數(shù)ψ、航行阻力Cx、最小壓力系數(shù)點位置xCpmin、逆壓梯度dCp/dx中的多個性能參數(shù)構(gòu)建目標(biāo)函數(shù),建立目標(biāo)優(yōu)化模型,以平衡各個性能指標(biāo),在設(shè)計變量的取值范圍內(nèi)搜索最佳設(shè)計點,使設(shè)計目標(biāo)的整體性能達(dá)到最優(yōu);使用CFD方法對變形后的AUV殼體表面網(wǎng)格進(jìn)行計算,得到所需的流體動力參數(shù)。

    但采用所提優(yōu)化框架進(jìn)行AUV殼體優(yōu)化時,其主要通過將AUV殼體表面網(wǎng)格嵌入到FFD控制體中以實現(xiàn)優(yōu)化過程中的AUV外形變化,不生成殼體的幾何外形,因此,基于幾何外形的傳統(tǒng)AUV殼體豐滿度計算方法不再適用。為此,本文提出了一種基于網(wǎng)格點的AUV 殼體豐滿度計算方法,在每一優(yōu)化迭代步中,根據(jù)變形后的AUV殼體表面網(wǎng)格對AUV殼體豐滿度進(jìn)行計算。該方法詳述如下:

    首先,將AUV 殼體x軸方向的總長度L均分為n等份,每段長度采用x坐標(biāo)區(qū)間表示為Uk=[kL/n,(k+1)L/n](k=0,…,n-1);對所有網(wǎng)格點(xi,yi,zi)(i=1,…,m)進(jìn)行分類,若網(wǎng)格點的x軸坐標(biāo)值xi位于區(qū)間Uk,則將其放置于集合Ak={(xi,yi,zi)|xi∈Uk,i=1,…,m}中。

    然后,針對每個集合A(kk=1,…,n-1)按照式(14)求解各個集合中所包含網(wǎng)格點的回轉(zhuǎn)半徑的平均值:

    式中,mk表示包含在幾何Ak的網(wǎng)格點個數(shù)。

    最后,按照式(15)計算AUV殼體的總體積和豐滿度為

    為了更為清晰、直觀地描述上述針對AUV殼體的自由變形優(yōu)化框架,圖5給出了該優(yōu)化框架的結(jié)構(gòu)圖。

    圖5 針對AUV殼體外形的自由變形優(yōu)化框架Fig.5 Free-form deformation optimization framework for AUV hull shape

    3 優(yōu)化結(jié)果與討論

    在本章中,以長度為5 m、半徑為0.533 m 的AUV 型殼體為例,分別采用傳統(tǒng)優(yōu)化設(shè)計方法和本文優(yōu)化設(shè)計方法對殼體外形進(jìn)行優(yōu)化設(shè)計,并對比分析所得的優(yōu)化結(jié)果,驗證所提方法的有效性。

    3.1 優(yōu)化問題描述

    首先,確定該AUV 殼體的參數(shù)化方法和優(yōu)化變量。對于傳統(tǒng)優(yōu)化設(shè)計方法,使用Nystrom、Myring和Granville 曲線對AUV 殼體形狀進(jìn)行參數(shù)化,其詳細(xì)幾何參數(shù)如表1 所示。考慮到在優(yōu)化之前已根據(jù)任務(wù)確定了圖1中AB、BC及CD段的軸向長度,選擇頭部曲線段Rh(x)和尾部曲線段Rt(x)中的可調(diào)參數(shù)作為優(yōu)化變量。

    表1 AUV殼體的主要幾何尺寸Tab.1 Main geometry parameters of the AUV hull

    對于本文提出的優(yōu)化設(shè)計方法,以表1所描述的AUV殼體作為初始幾何外形進(jìn)行表面網(wǎng)格劃分,并在表面網(wǎng)格外部設(shè)置長方體形FFD 控制體,具體如圖6所示。圖中,考慮到AUV 殼體平行中段(BC段)的長度和直徑在優(yōu)化過程中不需要更改,故僅在頭部曲線段(AB)和尾部曲線段(CD)布置FFD 控制點和回轉(zhuǎn)半徑比例因子進(jìn)行變形控制。對于頭部曲線段,如圖6 所示,共布置9 個回轉(zhuǎn)半徑比例因子,但為了保持頭部首端面和頭部曲線段、頭部尾端面和平行中段連接處的二階連續(xù)光滑,取中間的5個回轉(zhuǎn)半徑比例因子Ri(i=1,…,5)定義為頭部曲線段的全局優(yōu)化變量,以控制頭部曲線段在優(yōu)化過程中的幾何變形。對于尾部曲線段,如圖6所示,共設(shè)置6個回轉(zhuǎn)半徑比例因子,同樣為了保持連接處的二階連續(xù)光滑,取中間的2個回轉(zhuǎn)半徑比例因子Ri(i=6,7)定義為尾部曲線段的全局優(yōu)化變量,以控制尾部曲線段在優(yōu)化過程中的幾何變形。

    圖6 AUV殼體表面網(wǎng)格的FFD控制體Fig.6 FFD control volume outside the surface mesh of the AUV hull

    然后,確定AUV 殼體的優(yōu)化目標(biāo)函數(shù)。本文主要考慮將流動噪聲對制導(dǎo)系統(tǒng)的干擾降為最低且提升內(nèi)部裝載空間,故選用最小壓力系數(shù)點位置xCpmin和殼體豐滿度系數(shù)ψ兩個性能參數(shù)作為目標(biāo)函數(shù),且設(shè)定兩個性能參數(shù)同等重要,分配各個性能參數(shù)的權(quán)重為0.5,可得到AUV 殼體的目標(biāo)函數(shù)為0.5xCpmin+0.5ψ。

    最后,確定AUV 殼體的優(yōu)化約束。為了滿足設(shè)備的負(fù)載和運(yùn)行要求,要求AUV 殼體在軸向長度0.2 m、0.3 m、0.5 m 和0.75 m 處的回轉(zhuǎn)半徑必須不小于預(yù)設(shè)值。此外,優(yōu)化變量的可行區(qū)域也作為優(yōu)化約束被引入。

    綜上,表2詳細(xì)描述了上述四種優(yōu)化方法對應(yīng)的AUV殼體優(yōu)化問題。

    表2 不同優(yōu)化方法對應(yīng)的AUV殼體優(yōu)化問題Tab.2 Optimization problems of the AUV hull under different optimization methods

    3.2 優(yōu)化結(jié)果和比較

    對于表2 中的優(yōu)化問題,為了進(jìn)行全局優(yōu)化,選擇非支配排序遺傳算法(NSGA-II)作為四個優(yōu)化問題的優(yōu)化器,以控制整個優(yōu)化設(shè)計的運(yùn)行,并且設(shè)置相同的優(yōu)化參數(shù)和收斂條件。此外,為了減少優(yōu)化過程中的計算量,采用面元法計算目標(biāo)函數(shù)中包含的最小壓力系數(shù)點位置。

    表3列出了不同優(yōu)化方法獲得的最優(yōu)值。結(jié)果表明,相較于三種傳統(tǒng)的優(yōu)化設(shè)計方法,本文所提優(yōu)化方法在滿足所有約束的同時,取得了最大目標(biāo)值。產(chǎn)生此結(jié)果的原因是:基于自由變形的殼體參數(shù)化方法具有更大的優(yōu)化探索空間,可獲得更優(yōu)的殼體外形。

    表3 不同優(yōu)化方法得到的最優(yōu)值Tab.3 Optimal values obtained by different optimization design methods

    4 結(jié) 語

    為了更大程度地提高AUV 的殼體性能,本文提出了一種基于自由變形的AUV 殼體優(yōu)化設(shè)計方法。該方法基于B 樣條方法,建立了AUV 殼體的FFD 控制體,實現(xiàn)了AUV 殼體頭部曲線段和尾部曲線段的自由變形;并針對FFD控制體帶來的優(yōu)化設(shè)計變量規(guī)模過大等問題,提出一種全局參數(shù)化方法對FFD 控制體進(jìn)行間接變形操縱,在實現(xiàn)自由變形的前提下,大大減少參數(shù)化變量個數(shù),提高了優(yōu)化效率。此外,本文還給出了基于FFD 參數(shù)化的殼體豐滿度計算方法,搭建了基于自由變形的AUV 殼體優(yōu)化設(shè)計框架。

    最后,以一型5 m 長、0.533 m 寬的AUV 殼體為例,分別采用傳統(tǒng)優(yōu)化設(shè)計方法和本文提出的優(yōu)化設(shè)計方法對該殼體進(jìn)行了優(yōu)化設(shè)計,結(jié)果顯示相對于傳統(tǒng)優(yōu)化設(shè)計方法,本文優(yōu)化設(shè)計方法可以得到更好的優(yōu)化結(jié)果,更大地提升AUV殼體的性能,進(jìn)而證明了該方法的有效性,展示了其在AUV設(shè)計中的應(yīng)用前景。

    猜你喜歡
    型線殼體控制點
    減速頂殼體潤滑裝置的研制與構(gòu)想
    高次曲線組合型線渦旋盤性能研究*
    汽車變速箱殼體零件自動化生產(chǎn)線
    型線絞合導(dǎo)體ZC-YJLHV22-103×630鋁合金電纜的設(shè)計和生產(chǎn)
    電線電纜(2018年3期)2018-06-29 07:41:00
    NFFD控制點分布對氣動外形優(yōu)化的影響
    基于風(fēng)險管理下的項目建設(shè)內(nèi)部控制點思考
    油泵殼體的沖壓工藝及模具設(shè)計
    TH型干燥筒殼體的焊接
    焊接(2016年2期)2016-02-27 13:01:22
    相似材料模型中控制點像點坐標(biāo)定位研究
    變截面復(fù)雜渦旋型線的加工幾何與力學(xué)仿真
    国产av不卡久久| x7x7x7水蜜桃| 男女那种视频在线观看| 久久草成人影院| 成人亚洲精品av一区二区| 夜夜躁狠狠躁天天躁| 99久久精品一区二区三区| 国产亚洲精品一区二区www| www日本在线高清视频| 三级男女做爰猛烈吃奶摸视频| 欧美一区二区精品小视频在线| 亚洲精品日韩av片在线观看 | 99久久精品热视频| www.999成人在线观看| 国产真人三级小视频在线观看| av欧美777| 啦啦啦免费观看视频1| 欧美黑人欧美精品刺激| 国产黄色小视频在线观看| 成年人黄色毛片网站| 成年人黄色毛片网站| 精品日产1卡2卡| 啦啦啦观看免费观看视频高清| netflix在线观看网站| 精品不卡国产一区二区三区| 婷婷精品国产亚洲av在线| 亚洲精品成人久久久久久| 国产国拍精品亚洲av在线观看 | 中文字幕高清在线视频| 啦啦啦韩国在线观看视频| 真人做人爱边吃奶动态| 波多野结衣巨乳人妻| 久久精品国产亚洲av涩爱 | 91久久精品电影网| 亚洲成人中文字幕在线播放| 国产一区二区亚洲精品在线观看| 日韩av在线大香蕉| 国产淫片久久久久久久久 | 国产三级黄色录像| 欧美日韩乱码在线| x7x7x7水蜜桃| 99热精品在线国产| 成人高潮视频无遮挡免费网站| 免费在线观看影片大全网站| avwww免费| 女生性感内裤真人,穿戴方法视频| 日韩av在线大香蕉| 少妇的丰满在线观看| 人妻丰满熟妇av一区二区三区| 欧美三级亚洲精品| 18禁黄网站禁片免费观看直播| 久久久久久久久大av| 19禁男女啪啪无遮挡网站| 男女床上黄色一级片免费看| 97人妻精品一区二区三区麻豆| 母亲3免费完整高清在线观看| 色噜噜av男人的天堂激情| 免费av毛片视频| 国产黄a三级三级三级人| 九九热线精品视视频播放| 亚洲aⅴ乱码一区二区在线播放| 91在线精品国自产拍蜜月 | 国产黄片美女视频| 动漫黄色视频在线观看| 亚洲自拍偷在线| 久久精品影院6| 51午夜福利影视在线观看| 麻豆国产av国片精品| 成人亚洲精品av一区二区| 亚洲国产欧美网| 亚洲中文字幕日韩| 一a级毛片在线观看| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| 午夜激情福利司机影院| 99热这里只有精品一区| 国产激情偷乱视频一区二区| 国产免费男女视频| 亚洲中文字幕一区二区三区有码在线看| 午夜福利18| www.色视频.com| 国产精品爽爽va在线观看网站| 亚洲精品影视一区二区三区av| 51国产日韩欧美| 特大巨黑吊av在线直播| 久久国产精品人妻蜜桃| 国产精品一区二区免费欧美| www国产在线视频色| 午夜久久久久精精品| 在线观看舔阴道视频| h日本视频在线播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品一区av在线观看| 3wmmmm亚洲av在线观看| 国产成年人精品一区二区| 国产精品三级大全| 久久香蕉国产精品| 国产一区二区在线av高清观看| 麻豆国产97在线/欧美| 嫩草影院精品99| 欧美成狂野欧美在线观看| 成年免费大片在线观看| 精品一区二区三区av网在线观看| 在线观看舔阴道视频| 一个人看的www免费观看视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品在线美女| 精品不卡国产一区二区三区| 在线观看免费视频日本深夜| 99久久九九国产精品国产免费| av在线蜜桃| 国内精品美女久久久久久| 国产三级在线视频| 少妇丰满av| 亚洲精品粉嫩美女一区| 深夜精品福利| 国产国拍精品亚洲av在线观看 | 在线观看午夜福利视频| 国产高清视频在线播放一区| 日韩欧美 国产精品| 国产一区在线观看成人免费| 亚洲精品粉嫩美女一区| 欧美xxxx黑人xx丫x性爽| 亚洲人成电影免费在线| 少妇人妻精品综合一区二区 | 少妇丰满av| 床上黄色一级片| 亚洲aⅴ乱码一区二区在线播放| 黄片大片在线免费观看| 女人被狂操c到高潮| 看免费av毛片| 欧美日韩瑟瑟在线播放| 熟女电影av网| 国产三级在线视频| 久久久久久久午夜电影| 欧美色视频一区免费| 国产精品久久电影中文字幕| 好看av亚洲va欧美ⅴa在| 国产一级毛片七仙女欲春2| 国产主播在线观看一区二区| 精品久久久久久久毛片微露脸| 亚洲欧美一区二区三区黑人| 亚洲国产色片| 国产成人欧美在线观看| 亚洲国产欧美网| 国内少妇人妻偷人精品xxx网站| www.色视频.com| 久久亚洲真实| 特大巨黑吊av在线直播| 免费观看精品视频网站| 亚洲欧美激情综合另类| 国产视频内射| 成人精品一区二区免费| 一区二区三区高清视频在线| 国产伦精品一区二区三区视频9 | 女生性感内裤真人,穿戴方法视频| www日本黄色视频网| 在线观看一区二区三区| 露出奶头的视频| 老鸭窝网址在线观看| 日韩有码中文字幕| 国产欧美日韩精品一区二区| 国产精品久久久久久久电影 | 美女cb高潮喷水在线观看| 精品99又大又爽又粗少妇毛片 | 国产麻豆成人av免费视频| 国产又黄又爽又无遮挡在线| 老司机在亚洲福利影院| 国产不卡一卡二| 免费观看人在逋| 欧美色欧美亚洲另类二区| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品久久久com| 精品99又大又爽又粗少妇毛片 | 亚洲欧美日韩高清专用| 亚洲人成网站在线播| 中文字幕久久专区| 首页视频小说图片口味搜索| 91九色精品人成在线观看| 国产男靠女视频免费网站| 叶爱在线成人免费视频播放| 成人av在线播放网站| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久亚洲av鲁大| 午夜福利免费观看在线| 久久久久性生活片| 啦啦啦韩国在线观看视频| 人妻夜夜爽99麻豆av| 他把我摸到了高潮在线观看| 久久精品91无色码中文字幕| 精品久久久久久成人av| 99久久九九国产精品国产免费| 久久国产乱子伦精品免费另类| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 美女被艹到高潮喷水动态| 国产亚洲精品久久久久久毛片| 成人三级黄色视频| 国产精品精品国产色婷婷| 亚洲成a人片在线一区二区| 免费av观看视频| 亚洲av成人不卡在线观看播放网| 麻豆国产97在线/欧美| 欧美日韩福利视频一区二区| 亚洲av成人精品一区久久| 精品一区二区三区视频在线观看免费| 18禁国产床啪视频网站| 99riav亚洲国产免费| 亚洲第一电影网av| 在线观看66精品国产| 国产主播在线观看一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲不卡免费看| 手机成人av网站| 最新中文字幕久久久久| 亚洲黑人精品在线| 蜜桃久久精品国产亚洲av| 国产精品1区2区在线观看.| 丰满的人妻完整版| 亚洲五月婷婷丁香| 日韩大尺度精品在线看网址| 免费在线观看影片大全网站| 国产三级黄色录像| 免费在线观看成人毛片| 在线观看舔阴道视频| 在线视频色国产色| 一区二区三区免费毛片| 亚洲人成网站在线播| 欧美精品啪啪一区二区三区| 午夜福利欧美成人| 精品久久久久久久末码| 欧美大码av| 亚洲一区二区三区色噜噜| 亚洲内射少妇av| 国产真实乱freesex| 免费在线观看日本一区| 91在线观看av| 成人特级黄色片久久久久久久| 国产高清视频在线播放一区| 男女做爰动态图高潮gif福利片| 国产私拍福利视频在线观看| 精品久久久久久久末码| 99久久成人亚洲精品观看| 九九在线视频观看精品| 最新中文字幕久久久久| 日韩av在线大香蕉| 亚洲国产精品久久男人天堂| 在线观看免费视频日本深夜| 亚洲,欧美精品.| 成年免费大片在线观看| 日韩免费av在线播放| 国产精品 国内视频| 真人做人爱边吃奶动态| 亚洲熟妇中文字幕五十中出| 日本在线视频免费播放| 日韩欧美精品v在线| 日韩亚洲欧美综合| 久久精品91蜜桃| 51午夜福利影视在线观看| 免费一级毛片在线播放高清视频| 在线免费观看的www视频| 国产成人欧美在线观看| 国产爱豆传媒在线观看| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 亚洲五月天丁香| 91在线观看av| 19禁男女啪啪无遮挡网站| 国产在视频线在精品| 日韩欧美国产一区二区入口| av女优亚洲男人天堂| 亚洲va日本ⅴa欧美va伊人久久| 宅男免费午夜| 黄色片一级片一级黄色片| 亚洲欧美日韩高清专用| 性色avwww在线观看| 国产久久久一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 久久这里只有精品中国| 欧美成人一区二区免费高清观看| 国产高清三级在线| 亚洲专区中文字幕在线| 午夜亚洲福利在线播放| 国产久久久一区二区三区| 动漫黄色视频在线观看| 麻豆国产97在线/欧美| 日本撒尿小便嘘嘘汇集6| 少妇高潮的动态图| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 99久久99久久久精品蜜桃| 啦啦啦免费观看视频1| 99riav亚洲国产免费| 美女被艹到高潮喷水动态| 免费看a级黄色片| 夜夜爽天天搞| 蜜桃久久精品国产亚洲av| 国产午夜福利久久久久久| 国产亚洲欧美在线一区二区| 两人在一起打扑克的视频| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 国产精华一区二区三区| 国产精品电影一区二区三区| 精品久久久久久,| 欧美成人a在线观看| 99热精品在线国产| 国产精品影院久久| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 国产免费一级a男人的天堂| 欧美性猛交╳xxx乱大交人| 欧美极品一区二区三区四区| av天堂中文字幕网| 老汉色av国产亚洲站长工具| 色视频www国产| 三级男女做爰猛烈吃奶摸视频| 19禁男女啪啪无遮挡网站| 又爽又黄无遮挡网站| 国产精品98久久久久久宅男小说| 久久香蕉国产精品| 欧美一级毛片孕妇| 麻豆久久精品国产亚洲av| a在线观看视频网站| 亚洲色图av天堂| 精品免费久久久久久久清纯| 成人高潮视频无遮挡免费网站| 国产精品女同一区二区软件 | 美女 人体艺术 gogo| 成人av一区二区三区在线看| 1024手机看黄色片| 不卡一级毛片| 亚洲成av人片免费观看| 亚洲一区二区三区色噜噜| 黄片大片在线免费观看| 国产99白浆流出| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久 | 国产伦精品一区二区三区视频9 | 欧美在线黄色| 99久国产av精品| 国产三级在线视频| 老司机午夜福利在线观看视频| 村上凉子中文字幕在线| 欧美黄色淫秽网站| 18禁黄网站禁片免费观看直播| 首页视频小说图片口味搜索| 欧美日本亚洲视频在线播放| 国产成人欧美在线观看| 国产一区二区在线观看日韩 | 日本五十路高清| 一进一出抽搐动态| 久久欧美精品欧美久久欧美| 黄色片一级片一级黄色片| 国内少妇人妻偷人精品xxx网站| 最近在线观看免费完整版| 身体一侧抽搐| 欧美乱妇无乱码| 亚洲电影在线观看av| 美女cb高潮喷水在线观看| 国产精品久久久久久久电影 | 国产av在哪里看| 少妇人妻精品综合一区二区 | 99久久精品热视频| 观看免费一级毛片| 天堂√8在线中文| 欧美最黄视频在线播放免费| 国产一区二区亚洲精品在线观看| 亚洲成av人片在线播放无| 噜噜噜噜噜久久久久久91| 午夜视频国产福利| 舔av片在线| 亚洲无线观看免费| 观看美女的网站| 91在线观看av| 亚洲欧美一区二区三区黑人| 99久久九九国产精品国产免费| 免费人成视频x8x8入口观看| 欧美区成人在线视频| 男插女下体视频免费在线播放| 久久精品国产综合久久久| 午夜福利免费观看在线| 特级一级黄色大片| 日韩av在线大香蕉| 欧美区成人在线视频| 男插女下体视频免费在线播放| 日本黄色视频三级网站网址| 亚洲人成电影免费在线| 国产男靠女视频免费网站| 欧美黄色淫秽网站| 最新中文字幕久久久久| 国产又黄又爽又无遮挡在线| 天堂√8在线中文| 免费一级毛片在线播放高清视频| 极品教师在线免费播放| 757午夜福利合集在线观看| 免费搜索国产男女视频| 亚洲 欧美 日韩 在线 免费| 在线观看免费午夜福利视频| 国产精品影院久久| 特大巨黑吊av在线直播| www.色视频.com| 国产在视频线在精品| 在线免费观看的www视频| 午夜精品一区二区三区免费看| 亚洲内射少妇av| 小蜜桃在线观看免费完整版高清| 黄色视频,在线免费观看| 国产伦在线观看视频一区| 天美传媒精品一区二区| 亚洲精品一卡2卡三卡4卡5卡| 99riav亚洲国产免费| 国产99白浆流出| 在线观看免费午夜福利视频| 午夜影院日韩av| 久久久国产成人精品二区| 午夜福利欧美成人| 在线播放国产精品三级| 高清日韩中文字幕在线| 日本黄色片子视频| 欧美日本亚洲视频在线播放| 色av中文字幕| 亚洲,欧美精品.| xxx96com| 久久香蕉国产精品| 又黄又粗又硬又大视频| 国产精品久久久久久人妻精品电影| 国语自产精品视频在线第100页| 日韩欧美精品免费久久 | 国产国拍精品亚洲av在线观看 | 亚洲精品美女久久久久99蜜臀| 亚洲在线观看片| 国产精品一及| 可以在线观看毛片的网站| 久久久国产精品麻豆| 窝窝影院91人妻| 亚洲性夜色夜夜综合| 757午夜福利合集在线观看| 一区二区三区免费毛片| 午夜福利在线观看免费完整高清在 | 欧美激情在线99| 在线观看66精品国产| tocl精华| 国产精品久久久久久亚洲av鲁大| 成人亚洲精品av一区二区| 欧美午夜高清在线| 看片在线看免费视频| 国产亚洲精品综合一区在线观看| 啦啦啦免费观看视频1| 三级毛片av免费| 每晚都被弄得嗷嗷叫到高潮| 久久久久九九精品影院| 手机成人av网站| 岛国在线免费视频观看| 在线播放无遮挡| 亚洲内射少妇av| 国产成人a区在线观看| 久久精品国产清高在天天线| 欧美极品一区二区三区四区| 国产午夜精品久久久久久一区二区三区 | 欧美性猛交黑人性爽| 日韩大尺度精品在线看网址| 美女高潮喷水抽搐中文字幕| 国产在视频线在精品| 婷婷六月久久综合丁香| 国产单亲对白刺激| 男女午夜视频在线观看| 亚洲人与动物交配视频| 老汉色∧v一级毛片| 男插女下体视频免费在线播放| 日韩免费av在线播放| 国产单亲对白刺激| 黄片小视频在线播放| 亚洲中文日韩欧美视频| 亚洲黑人精品在线| 两个人的视频大全免费| 国产97色在线日韩免费| 男女做爰动态图高潮gif福利片| 高潮久久久久久久久久久不卡| 乱人视频在线观看| 在线观看66精品国产| 欧美日韩福利视频一区二区| 五月伊人婷婷丁香| 夜夜躁狠狠躁天天躁| av在线天堂中文字幕| 黄色成人免费大全| 可以在线观看的亚洲视频| 美女免费视频网站| 亚洲在线观看片| 深爱激情五月婷婷| av福利片在线观看| 51国产日韩欧美| 色噜噜av男人的天堂激情| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱| 成人亚洲精品av一区二区| 久久这里只有精品中国| 亚洲熟妇中文字幕五十中出| 欧洲精品卡2卡3卡4卡5卡区| 99久久精品热视频| 久久久久久久亚洲中文字幕 | 女人高潮潮喷娇喘18禁视频| 在线看三级毛片| 啦啦啦免费观看视频1| 法律面前人人平等表现在哪些方面| eeuss影院久久| 欧美另类亚洲清纯唯美| 五月伊人婷婷丁香| 校园春色视频在线观看| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 国产成年人精品一区二区| 亚洲av中文字字幕乱码综合| 欧美一区二区亚洲| 老司机在亚洲福利影院| 日本一二三区视频观看| 两个人的视频大全免费| xxxwww97欧美| 亚洲美女黄片视频| 内地一区二区视频在线| 色吧在线观看| 男女那种视频在线观看| 亚洲熟妇熟女久久| 天堂av国产一区二区熟女人妻| 欧美日韩中文字幕国产精品一区二区三区| 天堂av国产一区二区熟女人妻| 免费高清视频大片| 国产精品电影一区二区三区| 有码 亚洲区| 午夜日韩欧美国产| 欧美性感艳星| 国产探花极品一区二区| 亚洲av第一区精品v没综合| 亚洲国产精品合色在线| 一个人观看的视频www高清免费观看| 国产黄片美女视频| 国产精品久久视频播放| 国产精品久久久久久久电影 | av在线蜜桃| 久99久视频精品免费| 亚洲精品在线观看二区| 欧美黄色淫秽网站| eeuss影院久久| www.熟女人妻精品国产| 一进一出好大好爽视频| 亚洲欧美日韩高清在线视频| АⅤ资源中文在线天堂| 天堂av国产一区二区熟女人妻| 欧美性感艳星| 国产一区二区三区视频了| 首页视频小说图片口味搜索| 国产高清视频在线观看网站| 欧美激情久久久久久爽电影| 国产成人aa在线观看| 九九久久精品国产亚洲av麻豆| 亚洲av中文字字幕乱码综合| 免费看十八禁软件| 精品久久久久久,| 免费在线观看日本一区| 99久久无色码亚洲精品果冻| 成人国产一区最新在线观看| 国内久久婷婷六月综合欲色啪| 九九热线精品视视频播放| 国产精品亚洲一级av第二区| 国产97色在线日韩免费| 老汉色av国产亚洲站长工具| 成人无遮挡网站| 欧美+亚洲+日韩+国产| www.熟女人妻精品国产| 亚洲熟妇中文字幕五十中出| 亚洲成av人片在线播放无| 我要搜黄色片| 成人欧美大片| 日韩欧美 国产精品| 亚洲成人久久爱视频| 无遮挡黄片免费观看| 国产精品久久电影中文字幕| 深夜精品福利| 免费av观看视频| 淫秽高清视频在线观看| 亚洲国产欧美网| 露出奶头的视频| 91在线观看av| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频| 欧美午夜高清在线| 日韩大尺度精品在线看网址| 亚洲性夜色夜夜综合| 欧美绝顶高潮抽搐喷水| 母亲3免费完整高清在线观看| 亚洲电影在线观看av| 亚洲成人久久爱视频| 国产男靠女视频免费网站| 欧美一级毛片孕妇| 久久久久久久久大av| 看免费av毛片| 欧美一区二区亚洲| 亚洲av二区三区四区| 久久久色成人| 国产三级在线视频| 99热这里只有精品一区| 深夜精品福利| 欧美乱码精品一区二区三区| 免费搜索国产男女视频| 国内精品美女久久久久久| av片东京热男人的天堂| 久久这里只有精品中国| 一级毛片女人18水好多| 乱人视频在线观看| h日本视频在线播放| 麻豆国产av国片精品| 无遮挡黄片免费观看| 级片在线观看| 亚洲av日韩精品久久久久久密| 亚洲国产精品合色在线|