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

    基于多目標(biāo)模擬退火算法的微管網(wǎng)拓?fù)鋬?yōu)化

    2017-06-19 19:15:46黃金勇
    華東交通大學(xué)學(xué)報 2017年3期
    關(guān)鍵詞:份數(shù)沿程微管

    李 鵬,劉 遠(yuǎn),黃金勇

    (華東交通大學(xué)機(jī)電與車輛工程學(xué)院,江西 南昌330013)

    基于多目標(biāo)模擬退火算法的微管網(wǎng)拓?fù)鋬?yōu)化

    李 鵬,劉 遠(yuǎn),黃金勇

    (華東交通大學(xué)機(jī)電與車輛工程學(xué)院,江西 南昌330013)

    針對損傷自修復(fù)中微膠囊(或空心玻纖)載體在結(jié)構(gòu)上相互獨立、互不貫通等不足,研究基于多目標(biāo)模擬退火算法(multiobject simulated annealing,MOSA),對微管網(wǎng)的拓?fù)浣Y(jié)構(gòu)進(jìn)行優(yōu)化,旨在設(shè)計一種具有低體積份數(shù)和高流動效率的,具有相互貫通結(jié)構(gòu)特點的微管網(wǎng)載體。首先建立微管網(wǎng)數(shù)學(xué)模型;其次采用哈迪克羅斯迭代法設(shè)計了多目標(biāo)優(yōu)化函數(shù),并結(jié)合非劣分層和MOSA算法對微管網(wǎng)的拓?fù)浣Y(jié)構(gòu)進(jìn)行研究;最后與非劣分層遺傳算法(non-dominated sort genetic algorithm II,NSGA-II)進(jìn)行對比研究,驗證了MOSA的優(yōu)化性能。研究結(jié)果表明:微管網(wǎng)的體積份數(shù)和沿程水頭損失呈負(fù)相關(guān)性;與NSGA-II相比,MOSA的收斂時長略長,但MOSA解集的支配度及分布均勻程度均優(yōu)于NSGA-II。

    微管網(wǎng);多目標(biāo)優(yōu)化;模擬退火;非劣分層;哈迪克羅斯迭代

    復(fù)合材料因質(zhì)輕、耐磨及易加工等特點,已在航空航天、汽車、電子和建筑等領(lǐng)域得到了廣泛應(yīng)用。然而,作為一種由多種異質(zhì)/異形組元材料復(fù)合而成的新材料,其力學(xué)性能分散性較大,致使材料在生產(chǎn)和服役中易產(chǎn)生各類損傷,尤其是微裂紋的產(chǎn)生和擴(kuò)展將引起材料力學(xué)性能的劣化,導(dǎo)致構(gòu)件過早失效。目前,復(fù)合材料的損傷修復(fù)問題已經(jīng)成為結(jié)構(gòu)失效與安全服役科學(xué)中的關(guān)鍵問題之一[1]。

    上世紀(jì)80年代,美國軍方首先提出了復(fù)合材料自修復(fù)的概念,其基本原理是將內(nèi)置有修復(fù)劑的載體埋入材料內(nèi)部,當(dāng)材料受損產(chǎn)生裂紋時誘發(fā)載體破裂,釋放修復(fù)劑,實現(xiàn)材料的自修復(fù)[2]。此后,美國伊利諾伊大學(xué)Carolyn Dry研究了基于空心玻纖載體的自修復(fù)體系[3];White、Sottos和Moore等對基于微膠囊載體的自修復(fù)體系開展了大量研究,并論證了該方法的可行性[4-5];Bond等對自修復(fù)載體的性能進(jìn)行了研究[6]。在國內(nèi),哈爾濱工業(yè)大學(xué)的王榮國[7]和中山大學(xué)的章明秋[8-9]等先后也開展了基于微膠囊(或空心玻纖)載體的自修復(fù)研究;此外,本文作者也就空心玻纖載體與材料的兼容性等內(nèi)容開展了研究[10]。

    然而,目前自修復(fù)研究中普遍采用微膠囊和空心玻纖作為修復(fù)劑載體,由于載體彼此間互不貫通,限制了材料的損傷自修復(fù)效果(如,多次修復(fù),修復(fù)劑更換等)。因此,研究旨在設(shè)計一種具有相互貫通結(jié)構(gòu)特點的微管網(wǎng)載體,為保證載體與材料的兼容性和修復(fù)劑的流動效率,研究基于哈迪克羅斯迭代設(shè)計了多目標(biāo)函數(shù),并結(jié)合非劣分層和多目標(biāo)模擬退火(multi-object simulated annealing,MOSA)算法,對微管網(wǎng)的拓?fù)浣Y(jié)構(gòu)進(jìn)行了優(yōu)化,獲得了理想的結(jié)果。

    1 微管網(wǎng)優(yōu)化模型及目標(biāo)函數(shù)

    1.1 優(yōu)化模型

    研究以3×3正交微管網(wǎng)拓?fù)浣Y(jié)構(gòu)為對象,該管網(wǎng)由長度l為20 mm的24根管道構(gòu)成9個閉環(huán)回路,各管道的內(nèi)徑取值為決策空間,管道編號如圖1所示,其中復(fù)合材料的幾何尺寸為100 mm×100 mm×4 mm,S/T分別表示修復(fù)劑入/出口[11]。

    圖1 正交微管網(wǎng)模型Fig.1 Orthogonal network model

    1.2 多目標(biāo)優(yōu)化函數(shù)

    作為損傷自修復(fù)載體,埋入材料內(nèi)部的微管網(wǎng)要在不損壞材料初始性能的前提下,保證管道中的修復(fù)劑具有較小的沿程水頭損失(即具有良好的流動效率)。因此,微管網(wǎng)的拓?fù)鋬?yōu)化屬于多目標(biāo)優(yōu)化問題,即

    式中:φv為體積份數(shù);φf為沿程水頭損失;D為管道內(nèi)徑的解向量;dk為第k根管道的內(nèi)徑。

    1)體積份數(shù)φv。微管網(wǎng)占復(fù)合材料的體積份數(shù)φv可表示為

    式中:k為管道編號;l為管道長度;dk為第k根管道的內(nèi)徑;p為管道壁厚;V為復(fù)合材料體積。

    2)沿程水頭損失φf。修復(fù)劑流經(jīng)管道的沿程水頭損失表示了驅(qū)動修復(fù)劑流動所需能量的大小,為保證修復(fù)劑的流動效率,研究以全部管道的沿程水頭損失之和作為目標(biāo)函數(shù):

    式中:U為阻力因數(shù);n和m為修復(fù)劑材料系數(shù);Qk為第k根管道的流量,在忽略各環(huán)路之間相互影響的前提下,采用哈迪克羅斯(Hardy Cross)迭代方法獲取各管道的流量Qk[12]。

    1.3 非支配解

    微管網(wǎng)拓?fù)鋬?yōu)化的多目標(biāo)函數(shù) (即體積份數(shù)φv和沿程水頭損失φf)屬于非一致性優(yōu)化目標(biāo),即隨著體積份數(shù)下降,沿程水頭損失將呈現(xiàn)上升趨勢。因此,期望獲得解空間的最優(yōu)解是困難的,研究采用非劣分層方法,獲取對多目標(biāo)函數(shù)均有效的非支配解,即對于一個解集中任意解D1,該解集不存在解D2可同時滿足① φf(D2)≤φf(D1),且 φv(D2)≤φv(D1);② φf(D2)≤φf(D1)或 φv(D2)≤φv(D1),則該解集即為非支配解集,其中的解D1為非劣解。

    圖2 MOSA算法流程圖Fig.2 The flow chart of MOSA

    2 多目標(biāo)模擬退火算法

    多目標(biāo)模擬退火源于對熱力學(xué)中退火過程的模擬,是基于模擬退火 (simulated annealing,SA)算法,針對多目標(biāo)優(yōu)化的一種隨機(jī)尋優(yōu)算法。MOSA在搜索過程中結(jié)合逐漸趁于零且隨時間變化的概率突跳,實現(xiàn)了全局尋優(yōu)。MOSA算法包括初始解的生成、鄰域解的構(gòu)造和新解的產(chǎn)生等過程,如圖2所示。

    2.1 初始解集的生成

    以 24根管道的內(nèi)徑構(gòu)成解向量Di,j={dki,j},j∈(1,2,…,800),k=1,2,…,24,其中i為迭代次數(shù),j為當(dāng)前解中解的編號,k為管道編號。研究中采用隨機(jī)的方式產(chǎn)生規(guī)模為800的初始解集D0。

    2.2 鄰域解集的構(gòu)造

    2.3 新解集的產(chǎn)生及接受

    新解集Di+1通過對鄰域解中支配解與非支配解進(jìn)行不同處理而獲得,由以下兩部分組成:

    (2)第二部分。對于鄰域解中非支配解根據(jù)式(4)計算其內(nèi)能增量△E:

    式中:T和T′分別為本次和上次迭代的溫度值;λ為衰減因子。當(dāng)內(nèi)部循環(huán)次數(shù)達(dá)到鏈長L時,執(zhí)行降溫T=λT′。優(yōu)化中,溫度初值T0=240,衰減因子λ=0.99,鏈長L=10。

    2.4 終止條件

    優(yōu)化迭代的終止條件主要由退火溫度來控制,當(dāng)溫度達(dá)到精度要求時迭代終止。

    圖3 鄰域解構(gòu)造示意圖Fig.3 Schematic diagram of constructing neighborhood

    3 優(yōu)化結(jié)果及分析

    3.1 拓?fù)浣Y(jié)構(gòu)優(yōu)化結(jié)果

    通過MOSA算法對自修復(fù)微管網(wǎng)載體的拓?fù)浣Y(jié)構(gòu)進(jìn)行優(yōu)化,并獲得了其迭代過程及最終非支配解,如圖4和圖5所示。

    圖4 MOSA迭代過程Fig.4 The iterative process of MOSA

    圖5 MOSA和NSGA-II非支配解Fig.5 Non-dominated solution of MOSA and NSGA-II

    圖4中由下到上分別為MOSA第100,200,400,600,800,1 000次和1 233次迭代中的非支配解,結(jié)果表明:

    1)最終收斂迭代次數(shù)為1 233;

    2)隨著迭代次數(shù)增加,非劣解個數(shù)由64增至190;

    3)各次迭代所獲得的非支配解,其體積份數(shù)與沿程水頭損失均呈負(fù)相關(guān)性。

    圖5給出了第1 233次迭代的第一層190個非支配解,取非支配解(a),(b)和(c)進(jìn)行分析,其對應(yīng)的管網(wǎng)拓?fù)浣Y(jié)構(gòu)如圖6所示。

    圖6(a),(b)和(c)分別對應(yīng)于圖5中(a),(b)和(c)處的非支配解,線條的粗細(xì)代表管徑的大小,管道下方和左方為對應(yīng)管道的直徑,單位為mm。體積份數(shù)依次為2.73%,2.88%和3.25%的3個非支配解,其相應(yīng)的沿程水頭損失分別為2.810 6×10-6,1.601 0×10-6,7.769 0×10-6m。管網(wǎng)中內(nèi)徑較大的管道一般先從入口方向開始,并逐步向出口方向延伸,形成的大直徑通道對管道中的流體起到一定的疏導(dǎo)作用,降低水頭損失。

    圖6 拓?fù)浣Y(jié)構(gòu)優(yōu)化結(jié)果Fig.6 Topology structure of optimization resuits

    3.2 MOSA算法性能分析

    為驗證MOSA的優(yōu)化性能,在相同條件下,將多目標(biāo)退火算法與非劣分層遺傳算法(non-dominated sorting genetic algorithm II,NSGA-II)進(jìn)行了對比研究[13-14],主要包括:① 兩種算法之間非支配解的支配關(guān)系,即解的優(yōu)劣;②非支配解分布的均勻程度;③ 非支配解的收斂性。

    3.2.1 非支配解的支配關(guān)系

    不同解集之間的支配關(guān)系是評價解集優(yōu)劣的重要標(biāo)準(zhǔn),首先采用兩種算法分別獲得了非支配解(如圖6所示),其次,對兩個解集進(jìn)行非劣分層分析,并通過分析非支配解中,MOSA和NSGA-II解集所占的比例(PSA和PGA),實現(xiàn)兩個解集整體的優(yōu)劣程度的判斷,對比結(jié)果表明,MOSA和NSGA-II的非支配解共有252個,其中MOSA解集190個均為非支配解,NSGA-II解集62個為非支配解,占比分別為:PSA=75.4%,PGA=24.6%,顯然,采用MOSA算法獲得的解集顯著優(yōu)于后者。

    3.2.2 非支配解分布的均勻程度

    非支配解分布的均勻程度表征了解集的多樣性(即非劣解的選擇空間)。采用平均距離對非劣解的均勻程度進(jìn)行評價,計算公式如下:

    式中:φl為第l個目標(biāo)函數(shù);n為非支配解的規(guī)模。平均距離越短,表明非劣解分布越均勻。結(jié)果表明:① 以φv進(jìn)行評價,MOSA和NSGA-II的分布均勻程度為0.001 2和0.001 3,MOSA提高了8.3%;②以φf進(jìn)行評價,MOSA和NSGA-II的分布均勻程度分別為3.102 7×10-9和 5.188 3×10-9,MOSA提高了65.3%。

    3.2.3 非支配解的收斂性

    圖7中由下到上分別為NSGA-II第1 000,2 000,4 000,6 000,8 000,1 000次和12 000次迭代中獲得的非支配解,圖中X、Y軸分別為目標(biāo)函數(shù)φv、φf,Z軸為迭代次數(shù)。由于MOSA存在著內(nèi)部循環(huán),其時間復(fù)雜高于NSGA-II,因此NSGA-II的迭代次數(shù)設(shè)置為12 000。非支配解個數(shù)均達(dá)到解集規(guī)模設(shè)定值800。MOSA和NSGAII的收斂時長分別為19 985 s和18 674 s,NSGA-II略優(yōu)于MOSA。

    圖7 NSGA-II迭代過程Fig.7 The iterative process of NSGA-II

    4 結(jié)束語

    針對損傷自修復(fù)中微管網(wǎng)載體的拓?fù)浣Y(jié)構(gòu)優(yōu)化問題,研究基于哈迪克羅斯迭代法,設(shè)計了優(yōu)化目標(biāo)函數(shù),采用MOSA算法并結(jié)合非劣分層方法,研究了一種多目標(biāo)微管網(wǎng)拓?fù)浣Y(jié)構(gòu)優(yōu)化方法。

    1)采用哈迪克羅斯迭代法,修正了微管網(wǎng)中各管道的流量,獲得了表征微管網(wǎng)流動效率的沿程水頭損失,并以其作為目標(biāo)函數(shù)之一;

    2)結(jié)合MOSA和非劣分層方法,對微管網(wǎng)的拓?fù)浣Y(jié)構(gòu)進(jìn)行了優(yōu)化,經(jīng)1 233次迭代后,共獲得190個非支配解,結(jié)果表明:體積份數(shù)與沿程水頭損失呈現(xiàn)負(fù)相關(guān),且管網(wǎng)中內(nèi)徑較大的管道主要沿入口方向,并逐步向出口方向延伸;

    3)通過與NSGA-II的優(yōu)化結(jié)果進(jìn)行對比,結(jié)果表明:在NSGA-II的收斂時長略優(yōu)于MOSA的情況下,MOSA的優(yōu)化結(jié)果支配度75.4%優(yōu)于NSGA-II的24.6%;以φv和φf進(jìn)行解集分布均勻程度評價,MOSA分別提高了8.3%和65.3%。

    [1]祁恒治,趙蘊慧,朱孔營,等.自修復(fù)聚合物材料的研究進(jìn)展[J].化學(xué)進(jìn)展,2011(12):2560-2567.

    [2]王曉崗,張星,李原芃,等.自修復(fù)功能防腐涂膜研究進(jìn)展[J].功能材料,2012,43(19):2584-2587.

    [3]DRY C M.Self-repairing of composites[J].Proceeding of SPIE-the International Society for Optical Engineering,2003,5055:376-379.

    [4]OLUGEBEFOLA S C,ARAGON A M,HANSEN C J,et al.Polymer microvascular network composites[J].Composite Materials,2010,44(22):2587-2603.

    [5]PATRICK J F,SOTTOS N R,WHITE S R.Microvascular based self-healing polymeric foam[J].Polymer,2012,53(19):4231-4240.

    [6]KNIPPRATH C,MC COMBE G P,TRASK R S,et al.Predicting self-healing strength recovery using a multi-objective genetic algorithm[J].Composites Science&Technology,2012,72(6):752-759.

    [7]李海燕,王榮國,劉文博,等.微膠囊自修復(fù)聚合物材料[J].宇航材料工藝,2010,40(1):1-4.

    [8]章明秋,容敏智.結(jié)構(gòu)用自修復(fù)型高分子材料的制備[J].高分子學(xué)報,2012(11):1183-1199.

    [9]汪海平,容敏智,章明秋.微膠囊填充型自修復(fù)聚合物及其復(fù)合材料[J].化學(xué)進(jìn)展,2010,22(12):2397-2407.

    [10]李鵬.智能復(fù)合結(jié)構(gòu)損傷光修復(fù)與監(jiān)測的相關(guān)技術(shù)研究[D].南京:南京航空航天大學(xué),2009.

    [11]王憲杰,張洵安.基于宏觀性能的微觀多孔材料拓?fù)鋬?yōu)化[J].功能材料,2014,45(18):18078-18082.

    [12]丁祖榮,單雪雄,姜楫.流體力學(xué)[M].下冊.北京:高等教學(xué)出版社,2003:14-19.

    [13]曾毅,朱旭生,廖國勇.一種基于鄰域空間的混合粒子群優(yōu)化算法[J].華東交通大學(xué)學(xué)報,2013,30(3):44-49.

    [14]傅軍棟,鄒歡,康水華.PSO-SVM算法在智能建筑環(huán)境監(jiān)控系統(tǒng)中的應(yīng)用[J].華東交通大學(xué)學(xué)報,2016,33(1):121-127.

    Topology Optimization of Micro-vascular Network Based on Multi-objective Simulated Annealing Algorithms

    Li Peng,Liu Yuan,Huang Jinyong
    (School of Mechanotronics&Vehicle Engineering,East China Jiaotong University,Nanchang 330013,China)

    In the self-healing research,the microcapsule(or the hollow glass giber)is served as the carrier of self-healing agent.However,this carrier has some deficiencies such as being independent and uncorrelated with each other.In this paper,a novel carrier of micro-vascular network based on multi-objective simulated annealing(MOSA)algorithms was studied and topology of micro-vascular network was optimized to design a carrier with low volume fraction,high flow efficiency and inter-connective structure.Firstly,the mathematical model was established.Secondly,the multi-objective optimization function was designed by iteration of Hardy Cross.And then the topology of micro-vascular by MOSA algorithms combined with non-dominated sorting was studied.Finally,the superiority of MOSA was verified compared with non-dominated sort genetic algorithm (NSGA-II)algorithms.The results show that volume fraction and flow efficiency of micro-vascular network are negatively correlated.The convergence time of NSGA-II is slightly better than NSGA-II,but the degree of dominant and uniformity of MOSA are significantly better than that of NSGA-II.

    micro-vascular network;multi-objective optimization;simulated annealing;non-dominated sorting; iteration of Hardy Cross

    TB381

    A

    1005-0523(2017)03-0125-06

    (責(zé)任編輯 劉棉玲)

    2016-12-16

    國家自然科學(xué)基金項目(51365012);研究生創(chuàng)新專項資金項目(YC2015-S237)

    李鵬(1976—),男,副教授,博士,研究方向為智能結(jié)構(gòu),優(yōu)化算法。

    猜你喜歡
    份數(shù)沿程微管
    首張人類細(xì)胞微管形成高清圖繪出
    簡單和可控的NiO/ZnO孔微管的制備及對痕量H2S氣體的增強(qiáng)傳感
    不同微納米曝氣滴灌入口壓力下迷宮流道沿程微氣泡行為特征
    如何利用題組訓(xùn)練提高分?jǐn)?shù)“量”與“率”的區(qū)分度
    對提單及保單出具份數(shù)的思考
    中國外匯(2020年10期)2020-11-25 22:44:49
    典型生活垃圾爐排焚燒鍋爐沿程受熱面飛灰理化特性分析
    基于井下長管線沿程阻力損失的計算研究
    液壓與氣動(2020年5期)2020-05-22 03:34:40
    “份數(shù)法”的妙用
    胸腔微管引流并注入尿激酶治療結(jié)核性胸膜炎
    條條道路通羅馬
    欧美 日韩 精品 国产| 亚洲av欧美aⅴ国产| 欧美人与性动交α欧美精品济南到 | 哪个播放器可以免费观看大片| 亚洲精品国产一区二区精华液| 国产av国产精品国产| 曰老女人黄片| 十八禁高潮呻吟视频| 国产伦理片在线播放av一区| 日韩人妻精品一区2区三区| 色婷婷av一区二区三区视频| 我的亚洲天堂| 超碰97精品在线观看| 亚洲精品久久成人aⅴ小说| 精品一品国产午夜福利视频| 国产白丝娇喘喷水9色精品| 亚洲激情五月婷婷啪啪| 欧美国产精品一级二级三级| 97在线视频观看| 啦啦啦中文免费视频观看日本| 校园人妻丝袜中文字幕| 一二三四在线观看免费中文在| 亚洲精品一区蜜桃| 久久国产精品大桥未久av| 两性夫妻黄色片| 亚洲国产欧美在线一区| 大片免费播放器 马上看| 麻豆精品久久久久久蜜桃| 免费少妇av软件| 国产高清不卡午夜福利| 精品一区二区免费观看| 一二三四中文在线观看免费高清| 中文字幕亚洲精品专区| 亚洲欧美成人精品一区二区| 交换朋友夫妻互换小说| 有码 亚洲区| 夜夜骑夜夜射夜夜干| 黄色 视频免费看| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 爱豆传媒免费全集在线观看| 亚洲四区av| 久久精品国产亚洲av天美| 老汉色av国产亚洲站长工具| 日本午夜av视频| 国产日韩一区二区三区精品不卡| 青春草国产在线视频| 国产淫语在线视频| 国产精品 国内视频| 各种免费的搞黄视频| 热re99久久国产66热| 欧美亚洲日本最大视频资源| 看免费成人av毛片| 日韩精品有码人妻一区| 大片免费播放器 马上看| 欧美日韩亚洲高清精品| 久久女婷五月综合色啪小说| av片东京热男人的天堂| 夜夜骑夜夜射夜夜干| 高清在线视频一区二区三区| 亚洲人成77777在线视频| av免费观看日本| 黑人欧美特级aaaaaa片| 亚洲成人一二三区av| 日韩在线高清观看一区二区三区| 亚洲美女黄色视频免费看| 欧美另类一区| 亚洲三区欧美一区| 久久精品国产亚洲av涩爱| 国产精品偷伦视频观看了| 国产精品 欧美亚洲| www日本在线高清视频| 91成人精品电影| 久久女婷五月综合色啪小说| 久久精品国产自在天天线| 免费观看性生交大片5| 天堂俺去俺来也www色官网| 18+在线观看网站| 国产免费福利视频在线观看| 国产精品 欧美亚洲| 欧美日韩av久久| 青春草视频在线免费观看| 中文字幕人妻丝袜制服| 高清av免费在线| 热99国产精品久久久久久7| 天天躁夜夜躁狠狠躁躁| 人妻系列 视频| 一区二区三区四区激情视频| 久久这里只有精品19| 欧美日韩综合久久久久久| 亚洲av在线观看美女高潮| 99久国产av精品国产电影| 黑人猛操日本美女一级片| 亚洲一区二区三区欧美精品| 男女无遮挡免费网站观看| 午夜激情av网站| 丰满少妇做爰视频| 欧美中文综合在线视频| √禁漫天堂资源中文www| 精品一区在线观看国产| 最近中文字幕高清免费大全6| 一二三四中文在线观看免费高清| 国产精品99久久99久久久不卡 | 丝袜脚勾引网站| 一二三四在线观看免费中文在| 天堂俺去俺来也www色官网| 丝袜喷水一区| 欧美精品一区二区大全| 在线观看国产h片| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 免费人妻精品一区二区三区视频| 国产一区二区 视频在线| 韩国av在线不卡| 国产又色又爽无遮挡免| 五月伊人婷婷丁香| 美国免费a级毛片| 涩涩av久久男人的天堂| 满18在线观看网站| 国产视频首页在线观看| 国产成人欧美| 久久久久久久久久人人人人人人| 久久久久视频综合| 欧美精品亚洲一区二区| 99re6热这里在线精品视频| 一区二区三区精品91| 成人18禁高潮啪啪吃奶动态图| 亚洲精品成人av观看孕妇| 国产午夜精品一二区理论片| 青草久久国产| 91精品伊人久久大香线蕉| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区| 青草久久国产| 卡戴珊不雅视频在线播放| 久久久久久人妻| 国产白丝娇喘喷水9色精品| 亚洲欧美成人综合另类久久久| 啦啦啦视频在线资源免费观看| 三上悠亚av全集在线观看| 黄色配什么色好看| 日韩大片免费观看网站| 下体分泌物呈黄色| 在线精品无人区一区二区三| 日韩一卡2卡3卡4卡2021年| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 日韩一区二区三区影片| 中文字幕人妻丝袜一区二区 | 亚洲精品国产色婷婷电影| 麻豆精品久久久久久蜜桃| 国产精品二区激情视频| 欧美人与性动交α欧美精品济南到 | 韩国av在线不卡| 亚洲,欧美,日韩| 亚洲天堂av无毛| 男男h啪啪无遮挡| 在线亚洲精品国产二区图片欧美| 美女脱内裤让男人舔精品视频| av又黄又爽大尺度在线免费看| 色婷婷久久久亚洲欧美| 午夜福利,免费看| 美女主播在线视频| 高清欧美精品videossex| 男人添女人高潮全过程视频| 国产 一区精品| 免费观看无遮挡的男女| 日本欧美视频一区| 成人黄色视频免费在线看| 日本爱情动作片www.在线观看| 在线亚洲精品国产二区图片欧美| 亚洲av日韩在线播放| 夫妻午夜视频| 国产爽快片一区二区三区| 精品国产一区二区三区久久久樱花| av片东京热男人的天堂| 亚洲精品成人av观看孕妇| 亚洲男人天堂网一区| 亚洲色图 男人天堂 中文字幕| 久久99热这里只频精品6学生| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av高清一级| 亚洲精品国产一区二区精华液| 久久久久国产精品人妻一区二区| 久久鲁丝午夜福利片| 免费播放大片免费观看视频在线观看| 99久国产av精品国产电影| 精品一品国产午夜福利视频| 香蕉国产在线看| 亚洲,一卡二卡三卡| 久久人人97超碰香蕉20202| 人人妻人人澡人人爽人人夜夜| 国产一区二区在线观看av| 青春草亚洲视频在线观看| 亚洲精品在线美女| av.在线天堂| 亚洲精品国产色婷婷电影| 国产精品久久久久久久久免| 午夜福利在线免费观看网站| 丰满乱子伦码专区| 伊人久久国产一区二区| 青草久久国产| 亚洲精品日本国产第一区| 熟女av电影| 日韩中字成人| 一二三四在线观看免费中文在| 九草在线视频观看| 欧美中文综合在线视频| 欧美 日韩 精品 国产| 在线天堂最新版资源| 亚洲精品乱久久久久久| 国产1区2区3区精品| av国产久精品久网站免费入址| 在线观看www视频免费| 91成人精品电影| 欧美亚洲日本最大视频资源| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人欧美| 欧美精品高潮呻吟av久久| 人妻一区二区av| 一级片免费观看大全| 日韩大片免费观看网站| 热re99久久国产66热| 精品99又大又爽又粗少妇毛片| 国产精品国产三级专区第一集| 久久久久精品久久久久真实原创| 人人妻人人添人人爽欧美一区卜| 黄网站色视频无遮挡免费观看| 国产淫语在线视频| 亚洲美女黄色视频免费看| 日韩免费高清中文字幕av| kizo精华| 亚洲综合色惰| 欧美 日韩 精品 国产| 久久精品国产亚洲av高清一级| 欧美日韩视频高清一区二区三区二| 亚洲少妇的诱惑av| 如日韩欧美国产精品一区二区三区| 日日爽夜夜爽网站| 亚洲国产精品999| 99热网站在线观看| 久久久久久免费高清国产稀缺| 香蕉丝袜av| 青青草视频在线视频观看| av在线老鸭窝| 国产极品粉嫩免费观看在线| 天堂俺去俺来也www色官网| 美女脱内裤让男人舔精品视频| 精品第一国产精品| 久久狼人影院| 日日撸夜夜添| 999久久久国产精品视频| 在线看a的网站| 久久99热这里只频精品6学生| 青春草视频在线免费观看| 看非洲黑人一级黄片| 人妻少妇偷人精品九色| 色94色欧美一区二区| 亚洲精品国产色婷婷电影| 涩涩av久久男人的天堂| 婷婷色综合大香蕉| 一级黄片播放器| 亚洲成色77777| www.熟女人妻精品国产| 涩涩av久久男人的天堂| 免费播放大片免费观看视频在线观看| 亚洲精品在线美女| 麻豆av在线久日| 咕卡用的链子| 国产精品.久久久| 欧美精品av麻豆av| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 亚洲精品成人av观看孕妇| 你懂的网址亚洲精品在线观看| 久久久久久久久免费视频了| 国产精品国产av在线观看| 日本vs欧美在线观看视频| 丰满饥渴人妻一区二区三| 国产成人精品久久久久久| 久久青草综合色| 国产精品香港三级国产av潘金莲 | 亚洲一码二码三码区别大吗| 日韩欧美精品免费久久| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 秋霞伦理黄片| 久久久久国产一级毛片高清牌| 国产爽快片一区二区三区| 丝袜人妻中文字幕| 欧美日韩精品网址| 男女免费视频国产| 欧美人与性动交α欧美软件| 国产精品.久久久| 在线观看国产h片| 欧美97在线视频| 精品一区二区三区四区五区乱码 | 亚洲国产精品一区三区| 赤兔流量卡办理| 亚洲一级一片aⅴ在线观看| 超碰97精品在线观看| 永久免费av网站大全| 男男h啪啪无遮挡| 亚洲国产精品999| 一边亲一边摸免费视频| tube8黄色片| 国产国语露脸激情在线看| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| 777久久人妻少妇嫩草av网站| 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 国产精品不卡视频一区二区| 男女国产视频网站| 久久久精品免费免费高清| av免费在线看不卡| 久久精品国产鲁丝片午夜精品| 亚洲熟女精品中文字幕| 超碰97精品在线观看| 婷婷成人精品国产| 久久免费观看电影| 下体分泌物呈黄色| 国产福利在线免费观看视频| 热99国产精品久久久久久7| 夫妻性生交免费视频一级片| 久久国产亚洲av麻豆专区| 熟女av电影| 少妇的丰满在线观看| 老熟女久久久| 精品一区二区三区四区五区乱码 | 久久热在线av| 热99久久久久精品小说推荐| 男女高潮啪啪啪动态图| 两个人免费观看高清视频| xxx大片免费视频| 国产野战对白在线观看| 少妇人妻久久综合中文| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 国产无遮挡羞羞视频在线观看| 男女啪啪激烈高潮av片| 日韩 亚洲 欧美在线| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| 亚洲av日韩在线播放| 最近手机中文字幕大全| 人体艺术视频欧美日本| 咕卡用的链子| 精品国产乱码久久久久久男人| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 免费观看无遮挡的男女| 成人免费观看视频高清| 午夜日韩欧美国产| 人人妻人人澡人人爽人人夜夜| 免费观看无遮挡的男女| 五月开心婷婷网| 精品少妇久久久久久888优播| 黄片小视频在线播放| 亚洲精品国产色婷婷电影| 午夜免费观看性视频| 大码成人一级视频| 国产一区有黄有色的免费视频| 国产成人免费观看mmmm| 波多野结衣av一区二区av| 女人精品久久久久毛片| 国产精品麻豆人妻色哟哟久久| 极品人妻少妇av视频| 亚洲国产成人一精品久久久| 精品人妻偷拍中文字幕| 日日爽夜夜爽网站| 国产xxxxx性猛交| 一级,二级,三级黄色视频| 免费观看性生交大片5| 老女人水多毛片| 母亲3免费完整高清在线观看 | 一二三四在线观看免费中文在| 亚洲国产精品999| 女人精品久久久久毛片| 午夜激情av网站| 日韩熟女老妇一区二区性免费视频| 中文字幕最新亚洲高清| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 国产男女内射视频| 在线观看www视频免费| 黄色配什么色好看| 蜜桃在线观看..| 日韩精品有码人妻一区| 97精品久久久久久久久久精品| 久久久国产精品麻豆| 国产成人av激情在线播放| 人妻一区二区av| 久久人人爽av亚洲精品天堂| 精品少妇黑人巨大在线播放| 在线观看美女被高潮喷水网站| 国产成人一区二区在线| xxxhd国产人妻xxx| 美女高潮到喷水免费观看| 亚洲,一卡二卡三卡| www.av在线官网国产| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 91精品伊人久久大香线蕉| 国产免费现黄频在线看| 日韩免费高清中文字幕av| 黄频高清免费视频| 欧美国产精品va在线观看不卡| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 成人毛片60女人毛片免费| 国产成人免费无遮挡视频| 日韩一区二区视频免费看| 国产精品免费大片| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图| 伦精品一区二区三区| 91在线精品国自产拍蜜月| 黄色一级大片看看| 午夜91福利影院| 国产有黄有色有爽视频| 国产精品久久久久久精品电影小说| 一级片免费观看大全| 韩国av在线不卡| 91aial.com中文字幕在线观看| 晚上一个人看的免费电影| 又大又黄又爽视频免费| 国产野战对白在线观看| 亚洲av福利一区| 精品人妻在线不人妻| av天堂久久9| 国产精品二区激情视频| 热re99久久国产66热| 国产av码专区亚洲av| 老司机亚洲免费影院| 黄频高清免费视频| 最近的中文字幕免费完整| av在线观看视频网站免费| 国产成人a∨麻豆精品| 精品少妇黑人巨大在线播放| 亚洲婷婷狠狠爱综合网| 久久精品国产亚洲av高清一级| 国产 精品1| 亚洲av中文av极速乱| 性色av一级| 日韩中字成人| 在线观看人妻少妇| 国产精品免费大片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费av中文字幕在线| 久久精品aⅴ一区二区三区四区 | 91成人精品电影| 国产男人的电影天堂91| 一区二区日韩欧美中文字幕| 国产人伦9x9x在线观看 | 亚洲国产最新在线播放| 在线看a的网站| 精品少妇一区二区三区视频日本电影 | 国产精品av久久久久免费| 免费日韩欧美在线观看| 日韩av在线免费看完整版不卡| 久久精品久久久久久久性| 久久精品久久久久久噜噜老黄| 国产欧美日韩一区二区三区在线| 亚洲综合色惰| 久热这里只有精品99| 久久久久人妻精品一区果冻| 大香蕉久久网| 久久国内精品自在自线图片| 亚洲av国产av综合av卡| 在线亚洲精品国产二区图片欧美| 亚洲成国产人片在线观看| 久热这里只有精品99| 天堂8中文在线网| 国产一区二区激情短视频 | 亚洲精品aⅴ在线观看| 久久久久国产精品人妻一区二区| 亚洲色图 男人天堂 中文字幕| 韩国精品一区二区三区| 啦啦啦在线观看免费高清www| 另类亚洲欧美激情| 高清不卡的av网站| 你懂的网址亚洲精品在线观看| 国产成人aa在线观看| 国产精品一国产av| 亚洲国产精品国产精品| 日本猛色少妇xxxxx猛交久久| 欧美黄色片欧美黄色片| 午夜日韩欧美国产| 亚洲av欧美aⅴ国产| 妹子高潮喷水视频| 日本vs欧美在线观看视频| 成年动漫av网址| 在线观看美女被高潮喷水网站| 色94色欧美一区二区| 老司机影院成人| 大香蕉久久成人网| 丝瓜视频免费看黄片| 久久久国产一区二区| 高清视频免费观看一区二区| 天天躁日日躁夜夜躁夜夜| 亚洲av欧美aⅴ国产| 蜜桃国产av成人99| 18禁裸乳无遮挡动漫免费视频| 80岁老熟妇乱子伦牲交| 黄色毛片三级朝国网站| 考比视频在线观看| 国产1区2区3区精品| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 欧美bdsm另类| 久久狼人影院| 九色亚洲精品在线播放| 日韩中字成人| 亚洲三级黄色毛片| 最近中文字幕高清免费大全6| 肉色欧美久久久久久久蜜桃| 国产日韩欧美视频二区| 亚洲美女视频黄频| 日韩欧美一区视频在线观看| 久久精品久久精品一区二区三区| 亚洲欧美一区二区三区黑人 | 国产欧美日韩综合在线一区二区| 欧美在线黄色| 日本av手机在线免费观看| 肉色欧美久久久久久久蜜桃| 丰满饥渴人妻一区二区三| 亚洲av男天堂| 亚洲精华国产精华液的使用体验| 亚洲成人手机| 1024视频免费在线观看| www.自偷自拍.com| 欧美另类一区| 亚洲成色77777| 九色亚洲精品在线播放| 亚洲av成人精品一二三区| h视频一区二区三区| 久久99热这里只频精品6学生| 中文字幕人妻丝袜制服| 激情视频va一区二区三区| 七月丁香在线播放| 欧美激情 高清一区二区三区| 久久午夜福利片| 亚洲av国产av综合av卡| 97在线人人人人妻| 久久影院123| 叶爱在线成人免费视频播放| 97精品久久久久久久久久精品| 亚洲欧洲国产日韩| 在线观看一区二区三区激情| 观看美女的网站| 亚洲第一区二区三区不卡| 亚洲av综合色区一区| av又黄又爽大尺度在线免费看| 99re6热这里在线精品视频| 亚洲综合精品二区| 黄片无遮挡物在线观看| 高清在线视频一区二区三区| 在线观看国产h片| 在线天堂最新版资源| 在线观看www视频免费| 爱豆传媒免费全集在线观看| 国产探花极品一区二区| 热99国产精品久久久久久7| 国产熟女午夜一区二区三区| 精品国产一区二区三区四区第35| 亚洲 欧美一区二区三区| 边亲边吃奶的免费视频| 男人舔女人的私密视频| 亚洲欧洲日产国产| 人人妻人人澡人人爽人人夜夜| 亚洲五月色婷婷综合| 亚洲av男天堂| 国产毛片在线视频| 99久国产av精品国产电影| 欧美日韩精品成人综合77777| h视频一区二区三区| 99久久综合免费| 亚洲中文av在线| 免费黄色在线免费观看| 一级毛片 在线播放| 深夜精品福利| 热99国产精品久久久久久7| 99香蕉大伊视频| 秋霞伦理黄片| av国产久精品久网站免费入址| 免费人妻精品一区二区三区视频| 亚洲欧美色中文字幕在线| 天美传媒精品一区二区| 成人毛片a级毛片在线播放| 欧美成人午夜免费资源| 国产精品一区二区在线不卡| 99国产综合亚洲精品| 国产亚洲一区二区精品| 欧美精品国产亚洲| 岛国毛片在线播放| 狂野欧美激情性bbbbbb| 性少妇av在线| 亚洲国产精品999| 在线 av 中文字幕| 国产极品粉嫩免费观看在线| 热99久久久久精品小说推荐| 最近手机中文字幕大全| 欧美日韩精品成人综合77777| 午夜福利一区二区在线看| 亚洲欧美一区二区三区久久| 国产精品无大码| 亚洲精品久久成人aⅴ小说| 99九九在线精品视频| 日本av免费视频播放| 秋霞伦理黄片| 天天躁日日躁夜夜躁夜夜| 如日韩欧美国产精品一区二区三区|