• <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ú)的種群特征
    亚洲电影在线观看av| 又黄又粗又硬又大视频| www.熟女人妻精品国产| 又大又爽又粗| 美女扒开内裤让男人捅视频| 日本三级黄在线观看| 午夜久久久久精精品| √禁漫天堂资源中文www| 天堂av国产一区二区熟女人妻 | 国产精品一区二区精品视频观看| 少妇熟女aⅴ在线视频| 成年版毛片免费区| 欧美日韩黄片免| 欧美日本亚洲视频在线播放| 亚洲国产精品sss在线观看| 草草在线视频免费看| 中文字幕人妻丝袜一区二区| 天堂√8在线中文| a级毛片a级免费在线| 免费观看人在逋| 免费人成视频x8x8入口观看| 男插女下体视频免费在线播放| 久99久视频精品免费| 欧美色欧美亚洲另类二区| 国内久久婷婷六月综合欲色啪| 中文资源天堂在线| 亚洲人成网站在线播放欧美日韩| 国产亚洲av高清不卡| 亚洲专区中文字幕在线| 99久久精品热视频| 88av欧美| 亚洲七黄色美女视频| 99riav亚洲国产免费| 可以在线观看的亚洲视频| 久久午夜亚洲精品久久| 99久久国产精品久久久| 国产精品美女特级片免费视频播放器 | 欧美精品啪啪一区二区三区| 亚洲人成网站高清观看| 人人妻人人看人人澡| 色老头精品视频在线观看| 亚洲精品中文字幕在线视频| 日韩欧美精品v在线| 国产精品久久久久久精品电影| 欧美一区二区国产精品久久精品 | 欧美国产日韩亚洲一区| 欧美最黄视频在线播放免费| 搞女人的毛片| 777久久人妻少妇嫩草av网站| 美女大奶头视频| 香蕉丝袜av| 国产日本99.免费观看| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 露出奶头的视频| 亚洲欧美激情综合另类| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 18美女黄网站色大片免费观看| 亚洲男人的天堂狠狠| 岛国在线观看网站| 免费av毛片视频| 99riav亚洲国产免费| 成人18禁高潮啪啪吃奶动态图| 在线a可以看的网站| 久久这里只有精品中国| 99riav亚洲国产免费| 黄色丝袜av网址大全| av福利片在线观看| 国产精品久久久人人做人人爽| 久久亚洲真实| 中文字幕久久专区| 亚洲精品在线观看二区| 午夜亚洲福利在线播放| 午夜成年电影在线免费观看| 日韩中文字幕欧美一区二区| 国产三级黄色录像| 最新在线观看一区二区三区| av福利片在线观看| 亚洲精品国产精品久久久不卡| 国产欧美日韩一区二区精品| 亚洲av成人一区二区三| 欧美在线黄色| 俺也久久电影网| 国产aⅴ精品一区二区三区波| 欧美成人午夜精品| 国产三级黄色录像| 老司机福利观看| 亚洲国产精品合色在线| 宅男免费午夜| tocl精华| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 亚洲av第一区精品v没综合| 欧美日韩亚洲国产一区二区在线观看| 亚洲中文字幕一区二区三区有码在线看 | 深夜精品福利| 国产成人一区二区三区免费视频网站| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 久久精品91无色码中文字幕| 久久久国产成人免费| 中国美女看黄片| 欧美zozozo另类| 99riav亚洲国产免费| 午夜福利视频1000在线观看| 欧美中文日本在线观看视频| 一级毛片女人18水好多| 又紧又爽又黄一区二区| 国产野战对白在线观看| 亚洲精品美女久久av网站| a级毛片在线看网站| 欧美日韩福利视频一区二区| 日韩 欧美 亚洲 中文字幕| 国产三级黄色录像| 成人手机av| 国产激情久久老熟女| 国产高清videossex| 久久人妻av系列| 亚洲中文字幕日韩| 亚洲av电影在线进入| 桃色一区二区三区在线观看| 欧美成人午夜精品| 久久久久久大精品| 欧美国产日韩亚洲一区| 成人国产一区最新在线观看| 欧美zozozo另类| 在线播放国产精品三级| 国产爱豆传媒在线观看 | 亚洲av五月六月丁香网| 九色成人免费人妻av| 99久久精品国产亚洲精品| 757午夜福利合集在线观看| 免费av毛片视频| 欧美人与性动交α欧美精品济南到| 黄色毛片三级朝国网站| 亚洲成人精品中文字幕电影| 久久九九热精品免费| 一区二区三区国产精品乱码| av在线播放免费不卡| 欧美成人一区二区免费高清观看 | 黄色 视频免费看| 久久久久久久午夜电影| 成人三级做爰电影| 99久久久亚洲精品蜜臀av| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 国产av在哪里看| 亚洲欧美一区二区三区黑人| 九色国产91popny在线| 久久久久久九九精品二区国产 | 一边摸一边做爽爽视频免费| 国产精品av久久久久免费| 亚洲国产精品久久男人天堂| 亚洲精品中文字幕在线视频| 亚洲av熟女| 成人18禁在线播放| 最好的美女福利视频网| 丝袜人妻中文字幕| 国产成+人综合+亚洲专区| 国产不卡一卡二| 俺也久久电影网| 亚洲中文日韩欧美视频| 久久亚洲精品不卡| 老司机午夜十八禁免费视频| 一区二区三区激情视频| 我要搜黄色片| 脱女人内裤的视频| 亚洲国产看品久久| 18禁黄网站禁片午夜丰满| www日本在线高清视频| 丰满人妻一区二区三区视频av | 国产精品日韩av在线免费观看| 亚洲第一欧美日韩一区二区三区| 亚洲精品中文字幕在线视频| 啪啪无遮挡十八禁网站| 免费一级毛片在线播放高清视频| 美女午夜性视频免费| 十八禁人妻一区二区| 亚洲精品av麻豆狂野| 巨乳人妻的诱惑在线观看| 久久性视频一级片| 亚洲午夜精品一区,二区,三区| 日韩欧美三级三区| 午夜福利成人在线免费观看| 国产爱豆传媒在线观看 | 999久久久精品免费观看国产| 欧美性猛交╳xxx乱大交人| 蜜桃久久精品国产亚洲av| 妹子高潮喷水视频| 听说在线观看完整版免费高清| 法律面前人人平等表现在哪些方面| 欧美日韩黄片免| 我要搜黄色片| 日日摸夜夜添夜夜添小说| 午夜福利欧美成人| 久久久精品国产亚洲av高清涩受| 亚洲精品在线美女| 一个人观看的视频www高清免费观看 | 观看免费一级毛片| 美女大奶头视频| 给我免费播放毛片高清在线观看| 欧美性长视频在线观看| 91成年电影在线观看| 亚洲精品国产一区二区精华液| 99久久久亚洲精品蜜臀av| 亚洲av熟女| 国产激情欧美一区二区| 1024香蕉在线观看| 精品久久蜜臀av无| 好男人电影高清在线观看| 久久久久久久精品吃奶| 成人一区二区视频在线观看| av欧美777| 黄色片一级片一级黄色片| 欧美精品亚洲一区二区| 亚洲av电影在线进入| 国产人伦9x9x在线观看| 国产在线观看jvid| 免费在线观看成人毛片| 99久久久亚洲精品蜜臀av| 欧美成人性av电影在线观看| 男男h啪啪无遮挡| 一进一出好大好爽视频| 亚洲av电影不卡..在线观看| 国产精品影院久久| 真人一进一出gif抽搐免费| 香蕉久久夜色| 亚洲精品粉嫩美女一区| 日韩欧美国产在线观看| av在线播放免费不卡| 香蕉国产在线看| 亚洲男人的天堂狠狠| 成人高潮视频无遮挡免费网站| 亚洲专区中文字幕在线| 久久中文看片网| 男女之事视频高清在线观看| 欧美黄色淫秽网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av电影在线进入| 免费在线观看成人毛片| 亚洲专区字幕在线| 美女免费视频网站| 国产成人精品无人区| 亚洲精品国产一区二区精华液| 亚洲欧美精品综合一区二区三区| 他把我摸到了高潮在线观看| 后天国语完整版免费观看| 真人做人爱边吃奶动态| 亚洲美女黄片视频| av超薄肉色丝袜交足视频| 人人妻人人看人人澡| 丰满人妻一区二区三区视频av | 国产成人精品无人区| 99国产精品99久久久久| 一二三四社区在线视频社区8| 精品欧美一区二区三区在线| 久久久久久人人人人人| 99国产精品一区二区三区| 深夜精品福利| 色av中文字幕| 午夜影院日韩av| 麻豆成人午夜福利视频| 欧美成人午夜精品| 色精品久久人妻99蜜桃| 成人国产综合亚洲| 日本a在线网址| 中文字幕高清在线视频| 中文字幕人妻丝袜一区二区| www.精华液| 日韩欧美国产一区二区入口| 久久性视频一级片| 麻豆成人av在线观看| 久久人妻福利社区极品人妻图片| 日韩大尺度精品在线看网址| 香蕉av资源在线| 欧美3d第一页| 久久欧美精品欧美久久欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品中文字幕在线视频| 亚洲成av人片免费观看| 亚洲中文字幕日韩| a在线观看视频网站| 日韩精品中文字幕看吧| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 美女高潮喷水抽搐中文字幕| 久久久久精品国产欧美久久久| 美女免费视频网站| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| 亚洲人与动物交配视频| 两个人的视频大全免费| 真人一进一出gif抽搐免费| 一级黄色大片毛片| 免费搜索国产男女视频| 免费在线观看视频国产中文字幕亚洲| 国产精品自产拍在线观看55亚洲| 日韩欧美在线乱码| 一个人免费在线观看的高清视频| 啦啦啦韩国在线观看视频| 宅男免费午夜| 最近最新中文字幕大全电影3| 美女午夜性视频免费| av片东京热男人的天堂| av中文乱码字幕在线| 精品欧美国产一区二区三| 三级男女做爰猛烈吃奶摸视频| 国产黄色小视频在线观看| 国产1区2区3区精品| 深夜精品福利| 黄色a级毛片大全视频| 一进一出好大好爽视频| 中文字幕最新亚洲高清| 国产精品久久久人人做人人爽| 特大巨黑吊av在线直播| 日韩 欧美 亚洲 中文字幕| 国产97色在线日韩免费| 亚洲乱码一区二区免费版| 热99re8久久精品国产| 妹子高潮喷水视频| 精品欧美国产一区二区三| 一边摸一边做爽爽视频免费| 高清毛片免费观看视频网站| 夜夜看夜夜爽夜夜摸| 男女视频在线观看网站免费 | 久久九九热精品免费| 小说图片视频综合网站| 俺也久久电影网| 亚洲精品一卡2卡三卡4卡5卡| 一个人免费在线观看电影 | 亚洲乱码一区二区免费版| 热99re8久久精品国产| 午夜福利高清视频| 无人区码免费观看不卡| 亚洲专区中文字幕在线| 美女 人体艺术 gogo| 欧美一区二区精品小视频在线| 制服诱惑二区| 精品国产亚洲在线| 一本久久中文字幕| 国产一级毛片七仙女欲春2| 俺也久久电影网| 又紧又爽又黄一区二区| 黄片大片在线免费观看| 夜夜看夜夜爽夜夜摸| 99riav亚洲国产免费| 变态另类成人亚洲欧美熟女| 18美女黄网站色大片免费观看| 国语自产精品视频在线第100页| 黄频高清免费视频| 亚洲国产精品合色在线| 母亲3免费完整高清在线观看| 久久精品亚洲精品国产色婷小说| 国内毛片毛片毛片毛片毛片| 亚洲va日本ⅴa欧美va伊人久久| 午夜亚洲福利在线播放| 超碰成人久久| 亚洲欧美一区二区三区黑人| 白带黄色成豆腐渣| 啪啪无遮挡十八禁网站| 小说图片视频综合网站| 国模一区二区三区四区视频 | 久久这里只有精品中国| 日韩欧美国产在线观看| 深夜精品福利| 成人欧美大片| 国产亚洲欧美在线一区二区| 免费在线观看日本一区| 亚洲无线在线观看| 一本综合久久免费| av视频在线观看入口| 久久久久国产精品人妻aⅴ院| а√天堂www在线а√下载| 亚洲一区二区三区不卡视频| 在线观看美女被高潮喷水网站 | 床上黄色一级片| 国产精品电影一区二区三区| a级毛片a级免费在线| 亚洲av成人一区二区三| 日本免费a在线| 欧美日韩国产亚洲二区| 欧美成人性av电影在线观看| av免费在线观看网站| 中国美女看黄片| 欧美成人午夜精品| 99在线视频只有这里精品首页| 少妇粗大呻吟视频| 狂野欧美白嫩少妇大欣赏| 久久久久久九九精品二区国产 | avwww免费| 国产精品 欧美亚洲| 免费av毛片视频| 久久国产精品人妻蜜桃| 中文字幕人妻丝袜一区二区| 99国产精品99久久久久| 91av网站免费观看| 亚洲国产精品合色在线| 又大又爽又粗| 12—13女人毛片做爰片一| 可以在线观看的亚洲视频| 国产三级中文精品| www国产在线视频色| 国产精品久久久av美女十八| 99热这里只有精品一区 | 欧美黄色片欧美黄色片| 中出人妻视频一区二区| 午夜福利高清视频| 日本a在线网址| 蜜桃久久精品国产亚洲av| 悠悠久久av| 99久久无色码亚洲精品果冻| 99久久99久久久精品蜜桃| 波多野结衣高清作品| 日本在线视频免费播放| 少妇人妻一区二区三区视频| 国产欧美日韩一区二区精品| 国产乱人伦免费视频| 日韩欧美免费精品| 男男h啪啪无遮挡| 美女高潮喷水抽搐中文字幕| 亚洲av电影不卡..在线观看| 久久久久亚洲av毛片大全| 给我免费播放毛片高清在线观看| 曰老女人黄片| 亚洲狠狠婷婷综合久久图片| 91大片在线观看| 男女下面进入的视频免费午夜| 两人在一起打扑克的视频| 国内少妇人妻偷人精品xxx网站 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲全国av大片| 亚洲人成伊人成综合网2020| 国产一区二区在线观看日韩 | 欧美成人一区二区免费高清观看 | 变态另类丝袜制服| 欧美乱码精品一区二区三区| 免费在线观看影片大全网站| 无人区码免费观看不卡| 亚洲自拍偷在线| 国产精品久久电影中文字幕| 黄片小视频在线播放| 国产精品电影一区二区三区| 91av网站免费观看| 三级男女做爰猛烈吃奶摸视频| 久久午夜亚洲精品久久| 波多野结衣高清无吗| 国产69精品久久久久777片 | 国产成人精品久久二区二区91| 日韩精品免费视频一区二区三区| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 国产亚洲精品综合一区在线观看 | 在线观看免费视频日本深夜| 欧美日韩福利视频一区二区| 成人三级黄色视频| 午夜久久久久精精品| 99在线人妻在线中文字幕| 男男h啪啪无遮挡| 一夜夜www| 国产精品99久久99久久久不卡| 欧美乱妇无乱码| 18禁观看日本| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 日本一本二区三区精品| 天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 国产av又大| av中文乱码字幕在线| 亚洲avbb在线观看| 免费高清视频大片| 91九色精品人成在线观看| 啪啪无遮挡十八禁网站| 每晚都被弄得嗷嗷叫到高潮| 桃色一区二区三区在线观看| 丰满的人妻完整版| 99精品在免费线老司机午夜| 性色av乱码一区二区三区2| 欧美日韩精品网址| 午夜精品久久久久久毛片777| 亚洲精品在线观看二区| 日本在线视频免费播放| 久久久久久国产a免费观看| 日韩国内少妇激情av| 最近最新免费中文字幕在线| 老司机靠b影院| 首页视频小说图片口味搜索| 免费电影在线观看免费观看| 一区二区三区国产精品乱码| 99久久无色码亚洲精品果冻| 午夜视频精品福利| 老司机福利观看| 99精品欧美一区二区三区四区| 亚洲av成人一区二区三| 色精品久久人妻99蜜桃| 久久精品人妻少妇| 欧美一级毛片孕妇| 亚洲,欧美精品.| 国产亚洲精品综合一区在线观看 | 午夜视频精品福利| 色综合欧美亚洲国产小说| 美女高潮喷水抽搐中文字幕| 欧美日韩国产亚洲二区| 俺也久久电影网| 亚洲精品在线观看二区| 国产片内射在线| 亚洲熟妇熟女久久| 又黄又爽又免费观看的视频| 好男人在线观看高清免费视频| 欧美日本亚洲视频在线播放| 国产又黄又爽又无遮挡在线| 国产精品精品国产色婷婷| 久久国产精品影院| 国内精品久久久久久久电影| 他把我摸到了高潮在线观看| 99久久99久久久精品蜜桃| 色精品久久人妻99蜜桃| 亚洲av成人不卡在线观看播放网| 亚洲精品国产一区二区精华液| 亚洲精品美女久久av网站| 欧美性长视频在线观看| 1024手机看黄色片| 99re在线观看精品视频| 99精品在免费线老司机午夜| 成人午夜高清在线视频| 国产一区二区三区视频了| 又爽又黄无遮挡网站| 1024香蕉在线观看| netflix在线观看网站| 一级片免费观看大全| 三级国产精品欧美在线观看 | 免费在线观看视频国产中文字幕亚洲| 十八禁人妻一区二区| 婷婷精品国产亚洲av在线| 伦理电影免费视频| 黄色视频,在线免费观看| 国产视频一区二区在线看| 非洲黑人性xxxx精品又粗又长| 男人舔女人下体高潮全视频| 免费一级毛片在线播放高清视频| 亚洲一码二码三码区别大吗| 午夜福利高清视频| 欧美国产日韩亚洲一区| 亚洲狠狠婷婷综合久久图片| 国产亚洲精品久久久久久毛片| 亚洲成人久久爱视频| 国产高清激情床上av| 精品国产乱子伦一区二区三区| 日日干狠狠操夜夜爽| 成人av一区二区三区在线看| 国产精品久久电影中文字幕| 国产精品九九99| 天天一区二区日本电影三级| 在线观看午夜福利视频| 久久人人精品亚洲av| 亚洲人成网站高清观看| 亚洲七黄色美女视频| 久久久久久人人人人人| 亚洲精品国产一区二区精华液| av有码第一页| 一本大道久久a久久精品| 最近最新免费中文字幕在线| 免费在线观看影片大全网站| 法律面前人人平等表现在哪些方面| videosex国产| 三级男女做爰猛烈吃奶摸视频| av片东京热男人的天堂| 久久中文看片网| 亚洲一区二区三区不卡视频| 真人一进一出gif抽搐免费| 一a级毛片在线观看| 一级毛片高清免费大全| 色播亚洲综合网| 亚洲专区中文字幕在线| 国产高清videossex| 精品国产乱码久久久久久男人| 老汉色∧v一级毛片| 人妻久久中文字幕网| 久久久国产精品麻豆| 亚洲人成网站在线播放欧美日韩| 免费看日本二区| 国产精品九九99| 黄色成人免费大全| 最新在线观看一区二区三区| 国产成年人精品一区二区| 天天一区二区日本电影三级| 18禁黄网站禁片午夜丰满| 久久精品国产清高在天天线| 国产伦在线观看视频一区| 久久久国产成人免费| 久久婷婷人人爽人人干人人爱| 一级毛片高清免费大全| 99riav亚洲国产免费| 精品久久久久久,| 可以在线观看毛片的网站| 中国美女看黄片| 久久人人精品亚洲av| 成人一区二区视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 日本三级黄在线观看| 舔av片在线| 我的老师免费观看完整版| 特大巨黑吊av在线直播| 九九热线精品视视频播放| 亚洲专区国产一区二区| 又爽又黄无遮挡网站| 黄色视频不卡| 国产精品亚洲一级av第二区| 亚洲欧美精品综合久久99| 亚洲,欧美精品.| 激情在线观看视频在线高清| 99久久国产精品久久久| 88av欧美| 亚洲欧美精品综合一区二区三区|