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

    船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化算法

    2020-11-19 01:35:30劉冰潔畢曉君
    關(guān)鍵詞:約束條件設(shè)計(jì)方案種群

    劉冰潔,畢曉君

    (1.哈爾濱工程大學(xué) 信息與通信工程學(xué)院,黑龍江 哈爾濱 150001;2.中央民族大學(xué) 信息工程學(xué)院,北京 100081)

    隨著造船和航運(yùn)業(yè)的發(fā)展,船舶向著高速化、大型化、專用化發(fā)展,航道上船舶航運(yùn)密集程度增加,船舶操縱與控制的難度加大,發(fā)生碰撞、觸底等海難事故的概率及由此帶來(lái)的危害增大,提高船舶固有操縱性日漸顯示出其重要性,成為船型設(shè)計(jì)中至關(guān)重要的一環(huán)[1-2]。船舶操縱性優(yōu)化是一個(gè)多變量、多約束的復(fù)雜多目標(biāo)優(yōu)化問(wèn)題。傳統(tǒng)的優(yōu)化過(guò)程是定性分析比較和人工不斷湊試的過(guò)程,這種半經(jīng)驗(yàn)半理論的優(yōu)化方法通常帶有一定的盲目性,費(fèi)時(shí)費(fèi)力,很難為設(shè)計(jì)者提供多目標(biāo)優(yōu)化的設(shè)計(jì)方案。

    隨著人工智能的發(fā)展,群智能進(jìn)化算法逐漸應(yīng)用到船型優(yōu)化設(shè)計(jì)中,取得了一定的成果[3],但目前對(duì)船舶操縱性優(yōu)化的研究仍舊相對(duì)較少,其中多數(shù)研究都是將船舶操縱性優(yōu)化的多個(gè)目標(biāo)轉(zhuǎn)化成單目標(biāo)后進(jìn)行求解[4-6],這些方法由于事先設(shè)定了偏好信息,所求得的解受設(shè)計(jì)者先驗(yàn)知識(shí)的影響,不能覆蓋所有的解集,極易遺漏更好的解。還有學(xué)者研究利用多目標(biāo)進(jìn)化算法對(duì)此進(jìn)行求解,文獻(xiàn)[7]通過(guò)一種Pareto遺傳算法得到均勻分布的Pareto解集,但約束條件僅為設(shè)計(jì)變量的上下限,對(duì)船舶操縱性優(yōu)化數(shù)學(xué)模型的描述不夠準(zhǔn)確;文獻(xiàn)[8] 將NSGA-II和MOPSO應(yīng)用到船舶水動(dòng)力性能綜合優(yōu)化當(dāng)中,尋求適用于船舶水動(dòng)力性能優(yōu)化的多目標(biāo)進(jìn)化算法及其參數(shù)配置,但考慮的約束條件僅為設(shè)計(jì)變量的上下限,并未涉及約束處理技術(shù);文獻(xiàn)[9]采用解析目標(biāo)分流法,實(shí)現(xiàn)快速性、操縱性等5個(gè)學(xué)科組成的船舶總體設(shè)計(jì)優(yōu)化問(wèn)題,較好地解決了船舶總體設(shè)計(jì)學(xué)科間復(fù)雜耦合所帶來(lái)的計(jì)算困難,但其采用的約束處理技術(shù)為罰函數(shù)法,其中的參數(shù)較難設(shè)置,而且影響收斂精度;文獻(xiàn)[10]通過(guò)在基于分解的多目標(biāo)進(jìn)化算法(multi-objective evolutionary algorithm based on decomposition, MOEA/D)中引入自適應(yīng)約束處理技術(shù)和目標(biāo)值歸一化機(jī)制,獲得船舶水動(dòng)力性能優(yōu)化的Pareto解集,但其可行度參數(shù)僅隨種群中可行解比例變化,沒(méi)有考慮到進(jìn)化不同階段對(duì)不可行解的重視程度不同,且其約束處理技術(shù)僅在個(gè)體選擇階段,無(wú)法有效利用優(yōu)秀的不可行解,影響了優(yōu)化的精度;不能兼顧良好的多樣性和局部搜索能力,無(wú)法應(yīng)對(duì)多約束條件導(dǎo)致的復(fù)雜空間拓?fù)浣Y(jié)構(gòu),影響了優(yōu)化的精度。

    針對(duì)上述問(wèn)題,本文建立了船舶操縱性優(yōu)化問(wèn)題的數(shù)學(xué)模型,并提出一種船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化的方法:采用MOEA/D作為進(jìn)化算法框架,將差分進(jìn)化(differential evolution, DE)策略引入到MOEA/D中,并充分利用優(yōu)秀的不可行解對(duì)DE策略進(jìn)行改進(jìn),增強(qiáng)各搜索方向間的信息交互,提高船舶操縱性設(shè)計(jì)方案的收斂性能;其次,充分利用優(yōu)秀的不可行解,設(shè)計(jì)新的個(gè)體比較準(zhǔn)則,提高種群的分布性,使得船舶操縱性設(shè)計(jì)方案集在具有較好收斂性的同時(shí)獲得較好的分布性,給設(shè)計(jì)者提供多個(gè)不同需求下的優(yōu)化結(jié)果,設(shè)計(jì)者可以根據(jù)對(duì)操縱性衡準(zhǔn)指標(biāo)的不同需求自行選擇設(shè)計(jì)方案。

    1 船舶操縱性優(yōu)化的數(shù)學(xué)模型

    船舶操縱性優(yōu)化是一個(gè)多變量、多約束的復(fù)雜多目標(biāo)優(yōu)化問(wèn)題,不失公平性,本文與文獻(xiàn)[5]一致,以一艘載重量35 000 t原油船為例,討論船舶的操縱性優(yōu)化問(wèn)題。

    1.1 設(shè)計(jì)變量

    選取目前船舶操縱性優(yōu)化中常用的6個(gè)船型參數(shù)作為設(shè)計(jì)變量:船長(zhǎng)L、船寬B、吃水T、方形系數(shù)CB、舵的展高h(yuǎn)R和舵的弦長(zhǎng)bR[5]。

    1.2 目標(biāo)函數(shù)

    我國(guó)船舶標(biāo)準(zhǔn)化委員會(huì)操縱性分組通過(guò)的《海洋船舶操縱性標(biāo)準(zhǔn)》中,規(guī)定主要從船舶的回轉(zhuǎn)性、直線穩(wěn)定性和航向機(jī)動(dòng)性3個(gè)性能指標(biāo)來(lái)衡量船舶操縱性能[11]。由于航向機(jī)動(dòng)性目前在國(guó)內(nèi)外沒(méi)有統(tǒng)一的評(píng)價(jià)標(biāo)準(zhǔn),而且需要較多實(shí)驗(yàn)資料來(lái)表述,所以,與現(xiàn)有文獻(xiàn)一致,在本優(yōu)化模型中只選取相對(duì)回轉(zhuǎn)直徑D′和直線穩(wěn)定性指數(shù)C′作為衡量指標(biāo)[5]。

    1)直線穩(wěn)定性指數(shù)為:

    (1)

    2)相對(duì)回轉(zhuǎn)直徑D′的計(jì)算采用Lyster和Knights基于實(shí)船試驗(yàn)數(shù)據(jù)提出的相對(duì)回轉(zhuǎn)直徑回歸公式:

    D′=4.19-203CB/|δ|+47.4Trim/L-13B/L+

    194/δ-3.82(hR×bR)/(L×T)+7.79Ab/(L×T)

    (2)

    式中:δ為舵角;Trim為船舶傾斜;Ab為艏部浸濕面積。

    上述2個(gè)目標(biāo)性能指標(biāo)中,C′>0表示船舶在水平面內(nèi)的運(yùn)動(dòng)具有直線穩(wěn)定性,C′<0表示不足有直線穩(wěn)定性。直線穩(wěn)定性指數(shù)C′越大,表明船舶穩(wěn)定性越好,相對(duì)回轉(zhuǎn)直徑D′越小,表示船舶回轉(zhuǎn)性越好。

    1.3 約束條件

    根據(jù)文獻(xiàn)[5],船舶操作性優(yōu)化的11個(gè)約束條件為:

    (3)

    式中:Δ=CBρLBd是排水量;ρ=1 025 kg/m3為海水密度。

    1.4 船舶操縱性優(yōu)化數(shù)學(xué)模型

    綜上所述,船舶操縱性優(yōu)化數(shù)學(xué)模型描述為:

    minF(X)=[f1(X)f2(X)]T=[-C′D′]T

    s.tgj(X)>0,j=1,2,…,11

    (4)

    式中:X=[LBTCBhRbR]T為決策向量,F(xiàn)(X)為目標(biāo)向量,由于直線穩(wěn)定性指數(shù)C′越大越好,此處加上負(fù)號(hào)使其轉(zhuǎn)化為最小化問(wèn)題,gj(X)為式(3)所示的11個(gè)不等式約束條件。如果決策向量X滿足所有約束條件,則稱其為可行解;如果不完全滿足式(3)中的約束條件,則稱其為不可行解。根據(jù)這些約束條件計(jì)算個(gè)體的約束違反度函數(shù),如式(5)所示,其值稱為個(gè)體的約束違反度。

    (5)

    式中:當(dāng)gj(X)≤0,〈gj(X)〉返回其絕對(duì)值;當(dāng)gj(X)>0時(shí),〈gj(X)〉返回0。因此,當(dāng)X為可行解時(shí),G(X)=0;當(dāng)X為不可行解時(shí),G(X)>0。

    2 船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化算法

    本文建立的船舶操縱性優(yōu)化模型包含6個(gè)設(shè)計(jì)變量、2個(gè)目標(biāo)函數(shù),11個(gè)不等式約束條件,是一個(gè)多變量、多約束的非線性約束多目標(biāo)優(yōu)化問(wèn)題。為解決此問(wèn)題,本文提出一種船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化算法。采用MOEA/D作為多目標(biāo)進(jìn)化算法的框架,引入DE進(jìn)化策略,并充分利用優(yōu)秀的不可行解對(duì)DE變異算子進(jìn)行改進(jìn),增強(qiáng)各搜索方向間的信息交互,提高船舶操縱性優(yōu)化算法的收斂性能;其次,考慮到進(jìn)化不同時(shí)期不可行解的重要程度,設(shè)計(jì)新的個(gè)體選擇準(zhǔn)則,提高種群的分布性,使得船舶操縱性優(yōu)化的方案集在具有較好收斂性的同時(shí)獲得較好的分布性,給設(shè)計(jì)者提供多個(gè)不同需求下的優(yōu)化結(jié)果,設(shè)計(jì)者可以根據(jù)對(duì)C′和D′的不同需求自行選擇設(shè)計(jì)方案。

    2.1 基于分解的多目標(biāo)進(jìn)化算法框架

    MOEA/D[12]能夠產(chǎn)生均勻分布的Pareto最優(yōu)解集,更有利于向設(shè)計(jì)者提供均勻的、適合不同需求的設(shè)計(jì)方案。因此,本文采用MOEA/D進(jìn)化算法框架求解船舶操縱性優(yōu)化問(wèn)題。MOEA/D通過(guò)設(shè)定N個(gè)均勻分布的權(quán)重向量將一個(gè)m目標(biāo)問(wèn)題分解成m個(gè)單目標(biāo)子問(wèn)題,從而獲得N個(gè)均勻分布的Pareto最優(yōu)解。Tchebycheff 分解策略[12]是目前MOEA/D算法中最常用的一種分解策略,且該分解策略的收斂精度受Pareto前沿面形狀的影響不大,適合船舶操縱性優(yōu)化問(wèn)題Pareto前沿未知的情況。因此,本文采用Tchebycheff分解策略,對(duì)于船舶操縱性優(yōu)化問(wèn)題,其數(shù)學(xué)表達(dá)式為:

    (6)

    2.2 改進(jìn)差分變異策略

    DE是目前在種群多樣性和全局搜索能力方面效果較好的一種進(jìn)化算法[13],因此,本文將DE變異、交叉策略引入到MOEA/D框架中,用以求解船舶操縱性優(yōu)化問(wèn)題。考慮到原始的DE變異操作在處理設(shè)計(jì)變量維數(shù)較高、約束條件較多的約束優(yōu)化問(wèn)題時(shí),存在著對(duì)可行域的探索能力不強(qiáng),種群易陷入局部最優(yōu)的缺點(diǎn)[14],而船舶操縱性優(yōu)化問(wèn)題恰好是一個(gè)多變量、多約束的非線性約束多目標(biāo)優(yōu)化問(wèn)題。因此,需對(duì)DE變異操作進(jìn)行改進(jìn),以增強(qiáng)種群的搜索能力,從而提高船舶操縱性優(yōu)化問(wèn)題的求解精度。

    船舶操縱性優(yōu)化問(wèn)題由于存在多個(gè)約束條件,其可行域的拓?fù)浣Y(jié)構(gòu)變得極為復(fù)雜,種群能否跨越孤立的可行域、探索可行域邊界關(guān)系到船舶操縱性優(yōu)化問(wèn)題求解精度的高低。研究表明,進(jìn)化過(guò)程中有效利用優(yōu)秀不可行解所攜帶的有效信息,有助于種群跨越孤立的可行域,增強(qiáng)可行域邊界的探索能力,提高種群的多樣性和收斂性[15]。因此,本文利用優(yōu)秀不可行解對(duì)DE變異策略進(jìn)行改進(jìn),改進(jìn)的差分變異操作為:

    Vi=Xi+F1×(Xn-Xr1)+F2×(Xr3-Xr2)

    (7)

    式中:F1、F2為縮放因子;Xn為從當(dāng)代優(yōu)秀不可行解集A中隨機(jī)選取的不可行解;Xr1、Xr2、Xr3為Xi對(duì)應(yīng)的鄰域中隨機(jī)選擇的可行解;r1~r3為1~N上互不相等的正整數(shù)。

    在進(jìn)化的前期,種群中可能只包含不可行解,在這種情況下,應(yīng)使種群迅速探索可行域,找到更多的可行解。此時(shí),側(cè)重優(yōu)秀不可行解可以擴(kuò)大搜索的范圍,提高種群的多樣性,避免種群陷入局部最優(yōu),從而有助于搜索到更多的可行解。因此,縮放因子F1應(yīng)該較大,而F2應(yīng)該較小。在進(jìn)化的后期,種群中存在較多的可行解,為引導(dǎo)種群向Pareto前沿收斂,變異操作應(yīng)側(cè)重可行解,加強(qiáng)對(duì)Pareto最優(yōu)解附近區(qū)域的開(kāi)發(fā)。因此,縮放因子F1應(yīng)該較小,而F2應(yīng)該較大??s放因子及進(jìn)化前期和進(jìn)化后期的劃分根據(jù)實(shí)際問(wèn)題來(lái)確定,這里通過(guò)后邊的實(shí)驗(yàn)來(lái)確定。

    2.3 提出個(gè)體選擇策略

    雖然優(yōu)秀的不可行解參與到種群的進(jìn)化中,有助于擴(kuò)大探索范圍,提高種群的收斂性和多樣性。但是參與進(jìn)化的不可行解過(guò)多或者不夠“優(yōu)秀”,反而會(huì)影響種群收斂。因此,如何選擇可行解與不可行解,從而在保證種群多樣性的同時(shí)又能促使種群收斂到Pareto最優(yōu)前沿顯得尤為重要。

    為了更有效地利用優(yōu)秀的不可行解,兼顧種群的多樣性和局部搜索能力,應(yīng)對(duì)船舶操縱性優(yōu)化問(wèn)題中多約束條件導(dǎo)致的復(fù)雜空間拓?fù)浣Y(jié)構(gòu),提高優(yōu)化的精度,本文提出一種新的個(gè)體選擇準(zhǔn)則,進(jìn)化不同階段的選擇準(zhǔn)則對(duì)可行解和不可行解的側(cè)重不同。首先,引入約束違反度最大擴(kuò)增度ε(t)的概念,對(duì)個(gè)體需滿足的約束條件進(jìn)行松弛,讓約束違反度小于ε(t)的不可行解有機(jī)會(huì)參與到進(jìn)化中。接著,提出如式(8)所示的個(gè)體選擇準(zhǔn)則:如果2個(gè)個(gè)體的約束違反度均小于等于約束違反度最大擴(kuò)增度ε(t),或者2個(gè)個(gè)體的約束違反度相等,那么選擇Tchebycheff函數(shù)值小的個(gè)體,否則選擇約束違反度值小的個(gè)體。最后,提出如式(9)所示的“優(yōu)秀”不可行解選擇準(zhǔn)則:如果2個(gè)個(gè)體的約束違反度均大于約束違反度最大擴(kuò)增度ε(t),那么將Tchebycheff函數(shù)值小且約束違反度值大的優(yōu)秀不可行個(gè)體存儲(chǔ)在外部歸檔集A中。

    Ui,t代替Xi,t?

    (8)

    根據(jù)不同進(jìn)化階段,可行解與不可行解的重要性不同,本文引入的約束違反度最大擴(kuò)增度ε(t)是一個(gè)隨迭代次數(shù)動(dòng)態(tài)變化的值,如式(10)所示。

    Ui,t存儲(chǔ)到A?

    gte(Ui,t)G(Xi,t),

    (9)

    如果G(Ui,t),G(Xi,t)>ε(t)

    (10)

    式中:t為進(jìn)化迭代次數(shù);Gmax為最大進(jìn)化迭代次數(shù)。ε(0)為初始值,其設(shè)置方式如式(11)所示。

    (11)

    式中:N為種群規(guī)模;G(Xi)為初始種群個(gè)體的約束違反度。

    本文通過(guò)式(10)和式(11)動(dòng)態(tài)調(diào)節(jié)ε(t),逐漸加強(qiáng)對(duì)約束違反度擴(kuò)增程度的限制,在進(jìn)化初始,ε(t)為種群約束違反度的均值,種群中若只包含不可行解,或者不可行解距離可行域較遠(yuǎn),則ε(t)很大,可以讓更多的不可行解參與到進(jìn)化中,擴(kuò)大了搜索范圍,提高了種群多樣性。同時(shí),在進(jìn)化前期,ε(t)隨進(jìn)化迭代的進(jìn)行逐漸減小,從而不斷縮小滿足約束條件的區(qū)域,使種群逐漸收縮到可行域內(nèi)進(jìn)行搜索。在種群進(jìn)化的后期,約束違反度最大擴(kuò)增度ε(t)為0,避免過(guò)多的不可行解參與進(jìn)化,從而保障算法能夠收斂到可行的Pareto前沿,使求得的解都是滿足約束條件的。

    2.4 船舶操縱性優(yōu)化實(shí)現(xiàn)流程

    基于上面的2個(gè)創(chuàng)新點(diǎn),本文提出一種船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化算法,采用MOEA/D框架,有效利用優(yōu)秀不可行解對(duì)DE變異操作進(jìn)行改進(jìn),并提出一種個(gè)體選擇策略,其實(shí)現(xiàn)步驟如下所示,算法流程如圖1所示。

    圖1 船舶操縱性優(yōu)化流程Fig.1 Flow chart of ship maneuverability optimization

    1)初始化階段。

    ①設(shè)置實(shí)驗(yàn)參數(shù):種群大小N,權(quán)重鄰域大小T,最大進(jìn)化迭代次數(shù)Gmax,縮放因子F1、F2,交叉因子CR,確定船舶設(shè)計(jì)中變量的上限、下限;

    ②生成N個(gè)均勻分布的權(quán)重向量λ1,λ2,…,λN;

    ③根據(jù)任意2個(gè)權(quán)重向量之間的歐氏距離,求得距離每一個(gè)權(quán)重向量λi最近的T個(gè)權(quán)重向量的索引集合—鄰域B(i)={i1,i2,…,iT};

    ④隨機(jī)生成初始種群{X1,X2,…,XN};

    2)進(jìn)化階段。

    ①種群開(kāi)始進(jìn)化,令迭代次數(shù)t=1。

    ②從每一個(gè)B(i) 中隨機(jī)選取3個(gè)個(gè)體,從不可行解集A中隨機(jī)選取一個(gè)個(gè)體,與Xi按照式(7)進(jìn)行DE變異操作和交叉操作生成試驗(yàn)個(gè)體U;

    ③計(jì)算新生個(gè)體的目標(biāo)函數(shù)值。

    3)更新階段。

    ②個(gè)體比較階段:利用2.3節(jié)方法,若個(gè)體U與Xi比較,滿足式(8),則令Xi=U;若滿足式(9),將個(gè)體U與存儲(chǔ)在優(yōu)秀不可行解集A中。

    4)判斷終止條件。若t=Gmax,則算法停止并將種群中的Pareto最優(yōu)解作為結(jié)果輸出,否則,t=t+1,返回到2)。

    5)輸出。輸出所求得的Pareto解集及對(duì)應(yīng)的最優(yōu)解集合。

    3 仿真實(shí)驗(yàn)與結(jié)果分析

    為了驗(yàn)證本文方法的有效性和先進(jìn)性,將其與現(xiàn)有的求解船舶操縱性優(yōu)化效果較好的算法:自適應(yīng)加權(quán)算法[5]、Pareto遺傳算法[7]和基于MOEA/D的算法[10]進(jìn)行對(duì)比實(shí)驗(yàn)。所有實(shí)驗(yàn)在硬件配置為Intel(R) Core(TM)i7-7500U CPU@2.70 GHz 2.90 GHz、8G內(nèi)存,win10 64位操作系統(tǒng)的計(jì)算機(jī)上進(jìn)行,程序采用MATLAB R2016b編寫(xiě)。

    3.1 實(shí)驗(yàn)參數(shù)設(shè)置

    設(shè)置本文算的種群規(guī)模N=100,鄰域大小T=20,最大迭代次數(shù)Gmax=2 500,為保證實(shí)驗(yàn)結(jié)果的公平性,3種對(duì)比算法的種群規(guī)模N與最大迭代次數(shù)Gmax與本文算法相等。交叉率CR=0.9。不考慮縱傾,艏部浸濕面積Ab=15 m2,舵角δ=35°。設(shè)計(jì)變量的上下限如表1所示。

    表 1 決策變量的范圍Table 1 The scope of decision variables

    式(7)中縮放因子F1、F2需要人為設(shè)定,下面通過(guò)對(duì)其取不同的值對(duì)算法性能的影響來(lái)確定其取值。通過(guò)反復(fù)實(shí)驗(yàn),發(fā)現(xiàn)F1、F2在進(jìn)化前期取值對(duì)算法性能的影響,會(huì)隨進(jìn)化后期其取值的變化而變化,而F1、F2在進(jìn)化后期取值對(duì)算法性能的影響,也會(huì)隨著其在進(jìn)化前期的變化而變化。因此,下面通過(guò)正交試驗(yàn)設(shè)計(jì)確定進(jìn)化前、后期的縮放因子取值。固定進(jìn)化后期的縮放因子,并保證F1≥F2的前提下,F(xiàn)1、F2以0.3為基準(zhǔn),0.1為步長(zhǎng),獨(dú)立運(yùn)行28次實(shí)驗(yàn),選取Pareto解集不受其他解集支配的一組縮放因子(F1=0.8,F(xiàn)2=0.4)為最后的參數(shù)。其次,確定進(jìn)化后期的縮放因子:固定進(jìn)化前期的縮放因子,保證F1≤F2的前提下,F(xiàn)1、F2以0.3為基準(zhǔn),0.1為步長(zhǎng),獨(dú)立運(yùn)行28次實(shí)驗(yàn),選取Pareto解集不受其他解集支配的一組縮放因子(F1=0.4,F(xiàn)2=0.8)為最后的參數(shù)。

    3.2 對(duì)比結(jié)果分析

    為驗(yàn)證本文算法的有效性和先進(jìn)性,將其與3個(gè)目前求解船舶操縱性優(yōu)化效果較優(yōu)的算法:自適應(yīng)加權(quán)算法[5]、Pareto遺傳算法[7]和基于MOEA/D的算法[10]進(jìn)行對(duì)比實(shí)驗(yàn),3種對(duì)比算法與本文算法的實(shí)驗(yàn)對(duì)比結(jié)果分別如表2~4所示。需要注意的是,本文對(duì)Pareto最優(yōu)解集不進(jìn)行設(shè)計(jì)方案的決策,設(shè)計(jì)者可以根據(jù)對(duì)C′和D′的不同需求自行選擇設(shè)計(jì)方案。

    表2為本文算法和文獻(xiàn)[5]求得設(shè)計(jì)方案的對(duì)比。文獻(xiàn)[5]通過(guò)自適應(yīng)加權(quán)系數(shù),將2目標(biāo)優(yōu)化問(wèn)題轉(zhuǎn)化成單目標(biāo)優(yōu)化進(jìn)行尋優(yōu),表2中列出的文獻(xiàn)[5]方案是根據(jù)其給出的屬性權(quán)重值得到的最終最優(yōu)解,本文給出更多的設(shè)計(jì)方案,而且對(duì)文獻(xiàn)[5]給出的效果最好的2個(gè)方案(方案1和方案2),都能從本文方法求得的方案集中找到優(yōu)于它們的方案,如表2所示。對(duì)于文獻(xiàn)[5]給出的方案1,可從本文方案集中找到較其具有更大的直線穩(wěn)定性和更小的相對(duì)回轉(zhuǎn)直徑的方案,如方案1和方案2;同樣的對(duì)于文獻(xiàn)[5]給出的方案2,可從本文方案集中找到較其具有更大的直線穩(wěn)定性和更小的相對(duì)回轉(zhuǎn)直徑的方案,如方案3和方案4。由此可見(jiàn),本文求得的船舶操縱性優(yōu)化設(shè)計(jì)方案優(yōu)于文獻(xiàn)[5]自適應(yīng)加權(quán)方法所求得的設(shè)計(jì)方案。

    表 2 本文算法與文獻(xiàn)[5]方法結(jié)果對(duì)比Table 2 The comparison between the results of this paper and reference [5]

    表3為本文算法與文獻(xiàn)[7]求得設(shè)計(jì)方案的對(duì)比。文獻(xiàn)[7]沒(méi)有應(yīng)用約束處理技術(shù),為了比較公平性,將用文獻(xiàn)[7]方法求得的解集去除掉不滿足約束條件的解再與本文算法進(jìn)行比較。比較結(jié)果發(fā)現(xiàn)文獻(xiàn)[7]求得的方案均受本文算法的支配,而本文算法求得的所有解均不受文獻(xiàn)[7]的支配。為了比較的公平性,隨機(jī)從文獻(xiàn)[7]中選擇3個(gè)方案。對(duì)于文獻(xiàn)[7]求得的方案1,能從本文求得的方案集中找到方案1優(yōu)于它,即同時(shí)具有更大的直線穩(wěn)定性和更小的相對(duì)回轉(zhuǎn)直徑;而對(duì)于文獻(xiàn)[7]的方案2,可以從本文求得的方案集中找到2個(gè)方案(方案2和方案3)優(yōu)于它;對(duì)于文獻(xiàn)[7]中的方案3,可以從本文方案集中找到2個(gè)方案(方案4和方案5)優(yōu)于它。由此可見(jiàn),本文提出的船舶操縱性優(yōu)化設(shè)計(jì)方案優(yōu)于文獻(xiàn)[7]。

    表3 本文算法與文獻(xiàn)[7]方法對(duì)比Table 3 The comparison between the results of this paper and reference [7]

    表 4 本文算法與文獻(xiàn)[10]方法對(duì)比Table 4 The comparison between the results of this paper and reference [10]

    表4為本文算法與文獻(xiàn)[10]求得設(shè)計(jì)方案的對(duì)比。實(shí)驗(yàn)結(jié)果表明,對(duì)于文獻(xiàn)[10]求得的設(shè)計(jì)方案中有28個(gè)能在本文求得的方案中找到優(yōu)于他們的方案,而本文求得的所有方案都不差于文獻(xiàn)[10]方案。表4為28個(gè)方案中隨機(jī)選取的3個(gè)實(shí)例,從表中可以看出,本文求得的方案相比于文獻(xiàn)[10]具有更好的直線穩(wěn)定性和回轉(zhuǎn)性。由此可見(jiàn),本文提出的船舶操縱性設(shè)計(jì)方案優(yōu)于文獻(xiàn)[10]的。

    圖2為4種算法求得的Pareto最優(yōu)前沿,從圖2中可以看出,按照文獻(xiàn)[7]求得的方案都受到本文方案的支配,而本文求得的所有方案都不受文獻(xiàn)[7]方案的支配。文獻(xiàn)[10]求得的方案部分被本文方案中找到支配,而本文求得的所有方案都不差于文獻(xiàn)[10]方案,而且本文能求得更多的解,解集的分布性更好。

    圖2 4種算法Pareto最優(yōu)前沿對(duì)比Fig.2 Comparison of the Pareto optimal front obtained by the 4 algorithms

    為了比較本文算法與3個(gè)對(duì)比算法的優(yōu)化效率,表5給出了4種算法獨(dú)立運(yùn)行10次的平均時(shí)間。從表中可以看出文獻(xiàn)[5]、文獻(xiàn)[10]和本文算法由于不需要Pareto非支配比較,因而具有更快的優(yōu)化效率,但文獻(xiàn)[5]需要依靠專家經(jīng)驗(yàn),而文獻(xiàn)[10]對(duì)優(yōu)秀不可行的利用不夠充分,Pareto解集的收斂性沒(méi)有本文算法高。

    綜上所述,本文提出的船舶操縱性優(yōu)化方法求得的Pareto最優(yōu)解集收斂性更好,能夠使船舶具有更大的直線穩(wěn)定性、更小的相對(duì)回轉(zhuǎn)直徑;而且可以給設(shè)計(jì)人員提供分布更加均勻的設(shè)計(jì)方案,使設(shè)計(jì)者可以根據(jù)自己的偏好進(jìn)行選擇。

    4 結(jié)論

    1)本文提出一種船舶操縱性優(yōu)化的約束多目標(biāo)進(jìn)化算法,在優(yōu)化過(guò)程發(fā)現(xiàn)船舶的直線穩(wěn)定性和船舶的回轉(zhuǎn)性2個(gè)目標(biāo)是相互沖突的,多目標(biāo)進(jìn)化算法可以為設(shè)計(jì)者提供多種設(shè)計(jì)方案,供其根據(jù)不同的需求進(jìn)行選擇。

    2)本文算法充分利用優(yōu)秀不可行解,較對(duì)比算法的收斂性和分布性更好,可以給設(shè)計(jì)人員提供更多、更好的選擇。

    猜你喜歡
    約束條件設(shè)計(jì)方案種群
    邢氏水蕨成功繁衍并建立種群 等
    山西省發(fā)現(xiàn)刺五加種群分布
    基于可持續(xù)理念舊建筑改造設(shè)計(jì)方案探討
    基于一種改進(jìn)AZSVPWM的滿調(diào)制度死區(qū)約束條件分析
    數(shù)據(jù)中心ECC設(shè)計(jì)方案研究
    A literature review of research exploring the experiences of overseas nurses in the United Kingdom (2002–2017)
    線性規(guī)劃的八大妙用
    高壓電力系統(tǒng)規(guī)劃設(shè)計(jì)方案探討
    電子制作(2016年21期)2016-05-17 03:53:23
    某輕卡線束設(shè)計(jì)方案
    崗更湖鯉魚(yú)的種群特征
    成年人黄色毛片网站| 狂野欧美白嫩少妇大欣赏| 久久精品国产99精品国产亚洲性色| 国产精品 国内视频| 亚洲五月婷婷丁香| 免费在线观看亚洲国产| 成人18禁在线播放| 亚洲午夜理论影院| 精品电影一区二区在线| 十八禁网站免费在线| 天美传媒精品一区二区| 俺也久久电影网| 国产精品亚洲av一区麻豆| 国产探花在线观看一区二区| 亚洲中文字幕日韩| 在线国产一区二区在线| 在线看三级毛片| 午夜精品一区二区三区免费看| 尤物成人国产欧美一区二区三区| netflix在线观看网站| 伊人久久精品亚洲午夜| 欧美又色又爽又黄视频| 国产黄a三级三级三级人| 亚洲精品影视一区二区三区av| 99在线视频只有这里精品首页| 少妇的丰满在线观看| 啦啦啦韩国在线观看视频| 国产精品99久久99久久久不卡| 一个人看视频在线观看www免费 | 日本 av在线| 琪琪午夜伦伦电影理论片6080| 亚洲内射少妇av| 久久天躁狠狠躁夜夜2o2o| 两人在一起打扑克的视频| 黑人欧美特级aaaaaa片| 老熟妇乱子伦视频在线观看| 好男人电影高清在线观看| 99精品欧美一区二区三区四区| 俄罗斯特黄特色一大片| 欧美成狂野欧美在线观看| 老司机午夜福利在线观看视频| 人妻丰满熟妇av一区二区三区| 国产精品98久久久久久宅男小说| 国产三级在线视频| 午夜亚洲福利在线播放| 精品午夜福利视频在线观看一区| 91九色精品人成在线观看| 91久久精品电影网| 琪琪午夜伦伦电影理论片6080| 国产精品99久久久久久久久| 国产成人系列免费观看| 久久久久久久精品吃奶| 婷婷丁香在线五月| 三级男女做爰猛烈吃奶摸视频| 国产午夜精品久久久久久一区二区三区 | 国产一区二区三区视频了| 国产精品免费一区二区三区在线| 亚洲av第一区精品v没综合| 国产97色在线日韩免费| 免费看光身美女| 国产真实伦视频高清在线观看 | 亚洲国产精品合色在线| 看片在线看免费视频| 午夜两性在线视频| 一级a爱片免费观看的视频| 国产精品日韩av在线免费观看| 午夜亚洲福利在线播放| 内地一区二区视频在线| 国产精品98久久久久久宅男小说| 丁香欧美五月| 国产精品久久视频播放| 日本一本二区三区精品| 午夜两性在线视频| 久久中文看片网| 18禁裸乳无遮挡免费网站照片| 欧美色欧美亚洲另类二区| 真人做人爱边吃奶动态| 嫩草影视91久久| 欧美一区二区国产精品久久精品| 99精品久久久久人妻精品| 男女那种视频在线观看| 午夜福利在线观看免费完整高清在 | 嫩草影视91久久| 亚洲精品国产精品久久久不卡| 最近最新中文字幕大全免费视频| 中出人妻视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲乱码一区二区免费版| 床上黄色一级片| 麻豆成人午夜福利视频| 看黄色毛片网站| 在线观看一区二区三区| av视频在线观看入口| 国产v大片淫在线免费观看| 内射极品少妇av片p| 国产成人系列免费观看| 亚洲av一区综合| xxxwww97欧美| 亚洲av五月六月丁香网| 国产不卡一卡二| 99久久无色码亚洲精品果冻| 亚洲人成网站高清观看| 久久精品人妻少妇| 国产黄色小视频在线观看| 国产精品久久久久久久久免 | 亚洲最大成人中文| 九九久久精品国产亚洲av麻豆| 亚洲av一区综合| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 国产精品久久久人人做人人爽| 亚洲性夜色夜夜综合| aaaaa片日本免费| 亚洲一区二区三区不卡视频| 精品国产美女av久久久久小说| 久久亚洲真实| 天堂网av新在线| 久久国产精品人妻蜜桃| 成人国产综合亚洲| 久久精品国产清高在天天线| av女优亚洲男人天堂| 免费一级毛片在线播放高清视频| 国产男靠女视频免费网站| 性色av乱码一区二区三区2| 99久国产av精品| 国产精品98久久久久久宅男小说| 国产精品一区二区三区四区久久| 热99re8久久精品国产| 国产伦一二天堂av在线观看| 黄色女人牲交| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 人人妻,人人澡人人爽秒播| 一个人免费在线观看电影| 亚洲avbb在线观看| 国产视频一区二区在线看| 久久久久久久久大av| 极品教师在线免费播放| 母亲3免费完整高清在线观看| 一级作爱视频免费观看| 国产探花在线观看一区二区| 搡老岳熟女国产| 天天躁日日操中文字幕| 免费在线观看日本一区| 3wmmmm亚洲av在线观看| 国产淫片久久久久久久久 | 18+在线观看网站| 人人妻人人澡欧美一区二区| 日韩高清综合在线| 小蜜桃在线观看免费完整版高清| 亚洲国产精品999在线| 夜夜爽天天搞| 午夜视频国产福利| 亚洲片人在线观看| 国产成人aa在线观看| 国产亚洲精品av在线| 五月伊人婷婷丁香| 精品欧美国产一区二区三| 色综合婷婷激情| 少妇人妻一区二区三区视频| 在线观看午夜福利视频| netflix在线观看网站| 亚洲成人中文字幕在线播放| 午夜福利视频1000在线观看| 久久精品国产亚洲av涩爱 | 亚洲国产精品合色在线| 99久久精品国产亚洲精品| 欧美日韩瑟瑟在线播放| 亚洲av免费在线观看| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 午夜福利18| 中亚洲国语对白在线视频| 免费人成在线观看视频色| 久久精品91蜜桃| 亚洲美女黄片视频| 欧美乱妇无乱码| 最新美女视频免费是黄的| 亚洲中文字幕一区二区三区有码在线看| 不卡一级毛片| 色吧在线观看| 国产精品亚洲av一区麻豆| 在线国产一区二区在线| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美在线一区二区| 操出白浆在线播放| 在线观看66精品国产| 无遮挡黄片免费观看| 最后的刺客免费高清国语| 亚洲avbb在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美人成| 国产精品嫩草影院av在线观看 | 一本综合久久免费| 搡老岳熟女国产| 69人妻影院| 日本一二三区视频观看| 啦啦啦韩国在线观看视频| 成人av一区二区三区在线看| 亚洲av第一区精品v没综合| 精品日产1卡2卡| 亚洲午夜理论影院| 亚洲精品美女久久久久99蜜臀| 男女视频在线观看网站免费| 国产aⅴ精品一区二区三区波| 久久精品亚洲精品国产色婷小说| 我的老师免费观看完整版| 在线免费观看的www视频| 色精品久久人妻99蜜桃| 欧美日本视频| 午夜福利18| 久久午夜亚洲精品久久| 天堂av国产一区二区熟女人妻| 又黄又爽又免费观看的视频| 丰满的人妻完整版| 久久久久九九精品影院| 欧美精品啪啪一区二区三区| 久久久久久久久大av| 精品欧美国产一区二区三| 天堂网av新在线| 18禁黄网站禁片免费观看直播| 日韩欧美精品免费久久 | 日本一本二区三区精品| 搡老岳熟女国产| 亚洲一区二区三区色噜噜| 国产在线精品亚洲第一网站| 精品久久久久久久人妻蜜臀av| 午夜精品在线福利| 日韩亚洲欧美综合| 久9热在线精品视频| 18禁裸乳无遮挡免费网站照片| 亚洲内射少妇av| 亚洲成av人片在线播放无| 亚洲av中文字字幕乱码综合| 夜夜躁狠狠躁天天躁| 午夜精品一区二区三区免费看| 亚洲人与动物交配视频| 真人一进一出gif抽搐免费| 天天躁日日操中文字幕| 国产精品女同一区二区软件 | 黄色成人免费大全| 欧美+日韩+精品| 国产免费男女视频| 老司机深夜福利视频在线观看| 男女那种视频在线观看| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久| 窝窝影院91人妻| 男女那种视频在线观看| www.www免费av| 免费观看的影片在线观看| 亚洲美女视频黄频| a级毛片a级免费在线| 国产精品免费一区二区三区在线| 日韩欧美国产一区二区入口| 最好的美女福利视频网| 亚洲成人精品中文字幕电影| 欧美bdsm另类| 精品人妻1区二区| 黄片小视频在线播放| 少妇的逼好多水| 亚洲成人中文字幕在线播放| 亚洲av免费在线观看| 黄色女人牲交| 国产精品久久久久久精品电影| www国产在线视频色| 国产97色在线日韩免费| 校园春色视频在线观看| 亚洲人成电影免费在线| www.熟女人妻精品国产| 日本一本二区三区精品| 亚洲成人免费电影在线观看| 在线a可以看的网站| 国产高潮美女av| 综合色av麻豆| 国产私拍福利视频在线观看| 欧美黑人巨大hd| 麻豆国产av国片精品| 精品一区二区三区视频在线 | 99精品在免费线老司机午夜| 久久国产精品人妻蜜桃| 99热精品在线国产| 一级毛片高清免费大全| 国产高潮美女av| 久久草成人影院| 成人精品一区二区免费| 久久精品国产亚洲av涩爱 | 禁无遮挡网站| 日韩亚洲欧美综合| 欧美日韩福利视频一区二区| 国产在视频线在精品| 岛国在线观看网站| 成人特级黄色片久久久久久久| 国产伦精品一区二区三区四那| 亚洲无线观看免费| 日本熟妇午夜| 1000部很黄的大片| 此物有八面人人有两片| 精品一区二区三区视频在线观看免费| 国产免费男女视频| 在线a可以看的网站| 此物有八面人人有两片| 最新中文字幕久久久久| 国产精品日韩av在线免费观看| 免费高清视频大片| 国产精品精品国产色婷婷| 乱人视频在线观看| 男女视频在线观看网站免费| 中文字幕人妻丝袜一区二区| 五月伊人婷婷丁香| 久久精品国产清高在天天线| 精品久久久久久久毛片微露脸| 神马国产精品三级电影在线观看| av女优亚洲男人天堂| 国产视频内射| 精品人妻偷拍中文字幕| 99久久综合精品五月天人人| 搡老熟女国产l中国老女人| 99久久九九国产精品国产免费| 亚洲一区二区三区不卡视频| 精品不卡国产一区二区三区| 人妻夜夜爽99麻豆av| 久久精品国产亚洲av涩爱 | 色在线成人网| 欧美bdsm另类| 免费高清视频大片| 亚洲精品影视一区二区三区av| 久久欧美精品欧美久久欧美| 丰满乱子伦码专区| а√天堂www在线а√下载| 性色av乱码一区二区三区2| av片东京热男人的天堂| 色尼玛亚洲综合影院| 国产精品亚洲一级av第二区| 97碰自拍视频| 精品国产美女av久久久久小说| 亚洲国产精品sss在线观看| 日韩欧美一区二区三区在线观看| 12—13女人毛片做爰片一| 久久人人精品亚洲av| 在线观看一区二区三区| 成人欧美大片| 在线观看66精品国产| 999久久久精品免费观看国产| 两性午夜刺激爽爽歪歪视频在线观看| 日韩精品青青久久久久久| 1000部很黄的大片| 国产精品香港三级国产av潘金莲| 亚洲精品亚洲一区二区| 日韩欧美国产一区二区入口| 国产一区在线观看成人免费| 夜夜爽天天搞| 女警被强在线播放| 熟女少妇亚洲综合色aaa.| 看黄色毛片网站| 在线观看66精品国产| 欧美区成人在线视频| 国产野战对白在线观看| 中文字幕人妻丝袜一区二区| 成人午夜高清在线视频| 国产精华一区二区三区| 免费在线观看亚洲国产| 亚洲人成网站高清观看| 午夜免费男女啪啪视频观看 | 男女床上黄色一级片免费看| 日本熟妇午夜| 国产在线精品亚洲第一网站| 国产av一区在线观看免费| 国产探花在线观看一区二区| 国产亚洲精品av在线| 国产免费男女视频| 亚洲国产高清在线一区二区三| 欧美高清成人免费视频www| 中文字幕人妻熟人妻熟丝袜美 | 久久香蕉国产精品| 日韩av在线大香蕉| 夜夜夜夜夜久久久久| 熟女少妇亚洲综合色aaa.| 中文字幕av在线有码专区| av天堂中文字幕网| 成人三级黄色视频| 麻豆国产97在线/欧美| 久久亚洲真实| 一进一出抽搐动态| 国产久久久一区二区三区| 成人高潮视频无遮挡免费网站| 日本在线视频免费播放| 18+在线观看网站| 国产精品亚洲一级av第二区| 亚洲人成网站在线播放欧美日韩| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美在线乱码| 丰满的人妻完整版| 午夜福利成人在线免费观看| svipshipincom国产片| 亚洲久久久久久中文字幕| 午夜日韩欧美国产| 免费无遮挡裸体视频| 国内精品美女久久久久久| 我的老师免费观看完整版| 一区二区三区高清视频在线| 午夜免费成人在线视频| 久久精品国产清高在天天线| 日韩成人在线观看一区二区三区| 免费人成在线观看视频色| 久久久久国产精品人妻aⅴ院| 婷婷亚洲欧美| 成人特级黄色片久久久久久久| 久久久久久九九精品二区国产| 国产精品亚洲一级av第二区| 国产真实伦视频高清在线观看 | 国产69精品久久久久777片| 内地一区二区视频在线| 日韩欧美国产在线观看| 小蜜桃在线观看免费完整版高清| 黄色丝袜av网址大全| 国产色婷婷99| 香蕉久久夜色| 中文在线观看免费www的网站| 久久性视频一级片| 久久精品91无色码中文字幕| 性色av乱码一区二区三区2| 国产麻豆成人av免费视频| 丝袜美腿在线中文| 啦啦啦观看免费观看视频高清| 99久久99久久久精品蜜桃| 欧美日韩国产亚洲二区| 在线观看免费视频日本深夜| 搡女人真爽免费视频火全软件 | 久久精品人妻少妇| 18禁黄网站禁片免费观看直播| 日韩 欧美 亚洲 中文字幕| 亚洲精品在线美女| 俄罗斯特黄特色一大片| 性欧美人与动物交配| 999久久久精品免费观看国产| 嫩草影视91久久| av中文乱码字幕在线| 99热只有精品国产| 久久草成人影院| 麻豆一二三区av精品| 国产v大片淫在线免费观看| 美女被艹到高潮喷水动态| 色哟哟哟哟哟哟| 亚洲av免费在线观看| 最近在线观看免费完整版| 啦啦啦韩国在线观看视频| 成人国产一区最新在线观看| 美女大奶头视频| 美女cb高潮喷水在线观看| 国产高清有码在线观看视频| 99在线视频只有这里精品首页| 亚洲五月天丁香| 欧美日韩中文字幕国产精品一区二区三区| 亚洲最大成人中文| 精华霜和精华液先用哪个| 性色av乱码一区二区三区2| 美女高潮的动态| 精品免费久久久久久久清纯| aaaaa片日本免费| 在线十欧美十亚洲十日本专区| 在线观看66精品国产| 男人舔女人下体高潮全视频| 欧美在线一区亚洲| 成年女人毛片免费观看观看9| 听说在线观看完整版免费高清| 桃红色精品国产亚洲av| 国内久久婷婷六月综合欲色啪| 女警被强在线播放| 18+在线观看网站| 亚洲精品一区av在线观看| 国产黄色小视频在线观看| 国产成人aa在线观看| 国产亚洲精品一区二区www| 亚洲自拍偷在线| 成人av在线播放网站| 精品人妻一区二区三区麻豆 | 高清毛片免费观看视频网站| av黄色大香蕉| 在线a可以看的网站| 亚洲av美国av| av视频在线观看入口| 黄色片一级片一级黄色片| 少妇的逼水好多| 少妇人妻一区二区三区视频| 草草在线视频免费看| 国产黄a三级三级三级人| 亚洲精华国产精华精| 久99久视频精品免费| 成人永久免费在线观看视频| 久久精品国产自在天天线| 中国美女看黄片| 国内精品久久久久精免费| 日本与韩国留学比较| 舔av片在线| 老司机福利观看| av福利片在线观看| 美女黄网站色视频| 在线观看美女被高潮喷水网站 | 日韩欧美精品v在线| 天堂动漫精品| 一个人免费在线观看的高清视频| 欧美性猛交╳xxx乱大交人| 日本撒尿小便嘘嘘汇集6| 国产免费男女视频| 成人国产一区最新在线观看| 噜噜噜噜噜久久久久久91| 国产成人av教育| 18+在线观看网站| 亚洲美女黄片视频| 他把我摸到了高潮在线观看| 欧美色视频一区免费| 好看av亚洲va欧美ⅴa在| 最近视频中文字幕2019在线8| 午夜精品久久久久久毛片777| 国产成人av激情在线播放| 又粗又爽又猛毛片免费看| 久久久久久久精品吃奶| 免费观看精品视频网站| 一个人看视频在线观看www免费 | 亚洲av电影在线进入| 久久久久性生活片| 18禁国产床啪视频网站| 日本与韩国留学比较| 操出白浆在线播放| 日本撒尿小便嘘嘘汇集6| 欧美日韩国产亚洲二区| www日本黄色视频网| 真人一进一出gif抽搐免费| 天堂av国产一区二区熟女人妻| 国产精品乱码一区二三区的特点| 成年女人看的毛片在线观看| 日韩av在线大香蕉| xxx96com| 精品不卡国产一区二区三区| 麻豆国产97在线/欧美| aaaaa片日本免费| 三级男女做爰猛烈吃奶摸视频| 村上凉子中文字幕在线| 亚洲七黄色美女视频| 中文字幕人妻熟人妻熟丝袜美 | 欧美日本亚洲视频在线播放| 久久天躁狠狠躁夜夜2o2o| 少妇的逼水好多| 国内揄拍国产精品人妻在线| 久久婷婷人人爽人人干人人爱| 亚洲av第一区精品v没综合| 99热6这里只有精品| 91麻豆av在线| 国产极品精品免费视频能看的| 97碰自拍视频| 男女视频在线观看网站免费| 国产精品亚洲av一区麻豆| 精品电影一区二区在线| 夜夜爽天天搞| 一个人看的www免费观看视频| 久久国产精品影院| 黄片小视频在线播放| 欧美3d第一页| 精品不卡国产一区二区三区| 国产一区二区在线av高清观看| 综合色av麻豆| 日韩欧美在线二视频| 免费在线观看影片大全网站| 国产精品精品国产色婷婷| 人妻久久中文字幕网| 亚洲成人免费电影在线观看| 18禁在线播放成人免费| 国产av不卡久久| 国产97色在线日韩免费| 国产极品精品免费视频能看的| 亚洲成人久久性| 婷婷六月久久综合丁香| 1024手机看黄色片| 国语自产精品视频在线第100页| 内地一区二区视频在线| 网址你懂的国产日韩在线| 长腿黑丝高跟| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 亚洲国产欧洲综合997久久,| 97超级碰碰碰精品色视频在线观看| 色综合欧美亚洲国产小说| 欧美丝袜亚洲另类 | 啦啦啦免费观看视频1| 乱人视频在线观看| 观看免费一级毛片| 亚洲乱码一区二区免费版| 十八禁人妻一区二区| 熟妇人妻久久中文字幕3abv| 国产欧美日韩一区二区三| 亚洲av电影不卡..在线观看| 精品一区二区三区人妻视频| 长腿黑丝高跟| 国产 一区 欧美 日韩| 精品乱码久久久久久99久播| 少妇人妻精品综合一区二区 | 欧美成狂野欧美在线观看| 日本免费a在线| 中文字幕人妻熟人妻熟丝袜美 | 亚洲18禁久久av| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 国产伦人伦偷精品视频| 精品久久久久久久久久免费视频| 国产精品久久久久久久电影 | 中文字幕高清在线视频| 国产高清激情床上av| 一级黄片播放器| 精品国产超薄肉色丝袜足j| 在线观看免费午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 国产一区二区亚洲精品在线观看| 两人在一起打扑克的视频| 亚洲成人免费电影在线观看|