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

    基于自適應模擬退火法的UUV艏部線型優(yōu)化設計

    2011-05-27 08:52:40張宇文王亞東袁緒龍
    水下無人系統(tǒng)學報 2011年2期
    關鍵詞:聲功率線型模擬退火

    裴 譞, 張宇文, 王亞東, 袁緒龍

    ?

    基于自適應模擬退火法的UUV艏部線型優(yōu)化設計

    裴 譞, 張宇文, 王亞東, 袁緒龍

    (西北工業(yè)大學 航海學院, 陜西 西安, 710072)

    基于多目標空間的全局搜索理論, 將無人水下航行器(UUV)艏部豐滿度, 流體動力與流噪聲作為集成判定指標的綜合一體化UUV艏部線型優(yōu)化設計方法。建立了針對UUV艏部外形設計的一體化設計優(yōu)化模型, 結合經典粘流理論、CFD算法實現(xiàn)了集成優(yōu)化設計, 針對UUV外形優(yōu)化設計要求, 提出了復雜外形參數(shù)化和設計變量的選取原則, 并對具體算例進行外形綜合仿真優(yōu)化設計。結果表明, 基于自適應模擬退火法能夠有效實現(xiàn)對UUV期望目標的一體化優(yōu)化設計, 在相互制約的多種目標函數(shù)所組成的目標域中獲得最優(yōu)解, 優(yōu)化后UUV艏部流體動力和聲學性能都有較大提高。

    無人水下航行器(UUV); 多目標優(yōu)化; 模擬退火法; 艏部線型; 集成優(yōu)化設計

    0 引言

    無人水下航行器(unmanned undersea vehicle, UUV)是未來海洋開發(fā)和海洋軍事領域應用不可缺少的工具[1], 具有極為廣闊的發(fā)展前景。其中外形設計是總體綜合優(yōu)化設計的重要組成部分, 直接關系到UUV諸多戰(zhàn)術技術性能, 也是研究艇體結構水動力外載荷和水動力噪聲的依據和基礎[2]。UUV艏部線型設計一般基于雙參數(shù)方程[3], 其計算域理論上趨于無窮大, 且解空間具有明顯的多峰性、離散性和耦合性, 目標函數(shù)存在很強的相互制約關系, 往往為解空間的強函數(shù)[4]。因此,對優(yōu)化技術的全局搜索路徑與決策方向提出了更高的要求。文獻[5]和[6]指出, 現(xiàn)有的設計優(yōu)化方法一般采用線性搜索的方法, 將阻力或噪聲作為單一目標函數(shù)或將多目標函數(shù)離散線性化, 依照目標優(yōu)先級作為單一函數(shù)優(yōu)化, 以期獲得設計參數(shù)之間的協(xié)調綜合, 但是線型搜索法所獲得的解空間區(qū)域的全局化程度較差, 其搜索方向往往趨向于局部化, 從而造成優(yōu)化結果不理想。本文基于自適應模擬退火算法, 將UUV艏部的邊界層轉捩、阻力系數(shù)、豐滿度與表面聲功率級作為優(yōu)化目標函數(shù)集進行一體化設計, 從而獲得滿足該目標凸空間域的UUV外形期望最優(yōu)解。

    1 優(yōu)化目標模型

    一般意義下, 多個目標之間往往是相互沖突的, 即在目標投影凸空間內, 不存在一個最優(yōu)解使所有目標同時達到最優(yōu)。一個目標性能的改善往往以其他一個或多個目標性能的降低為代價[6]。要在滿足給定約束條件的前提下, 從設計變量的取值范圍內搜索最佳設計點, 使多個設計目標決定的設計對象整體性能達到最優(yōu), 就要適當選擇多目標的優(yōu)化方法。

    多目標優(yōu)化問題可描述為[7]: 尋找一組設計變量=(1,2,…,x)T, 使max:[f()]=1, …,,約束條件為

    式中:,,,分別為設計變量、目標函數(shù)、等式約束和不等式約束的個數(shù);為優(yōu)化問題的設計變量組成的向量。

    模擬退火算法來源于固體退火原理, 將固體加溫至充分高, 再讓其徐徐冷卻, 加溫時, 固體內部粒子隨溫度升高變?yōu)闊o序狀, 內能增大, 而冷卻時粒子漸趨有序, 在每個溫度都達到平衡態(tài), 最后在常溫時達到基態(tài), 內能減為最小。假設在狀態(tài)old時, 系統(tǒng)受到某種擾動而使其狀態(tài)變?yōu)閚ew, 對應的系統(tǒng)能量也從(old)變成(new), 系統(tǒng)由狀態(tài)old變?yōu)闋顟B(tài)new的接受概率為[8]。模擬退火算法與初始值無關, 算法求得的解與初始解狀態(tài)(是算法迭代的起點)無關; 模擬退火算法具有漸近收斂性, 已在理論上被證明是一種以概率收斂于全局最優(yōu)解的全局優(yōu)化算法。

    2 優(yōu)化設計系統(tǒng)

    UUV艏部線型設計準則的提出采用具體設計具體分析的方法, 從總體角度將阻力、轉捩點、豐滿度和噪聲作為判定函數(shù), 明確所設計的艏部線型的主要指標要求, 依據解決主要矛盾的原則提出?;谧赃m應模擬退火算法, 利用并行融合優(yōu)化的設計思路, 建立了包括邊界層轉捩點、UUV艏部流體動力計算, UUV艏部表面聲功率級計算與外形參數(shù)一體化融合優(yōu)化設計的綜合優(yōu)化設計系統(tǒng)。UUV艏部線型優(yōu)化設計的目標函數(shù)由3個目標參數(shù)與1個約束參數(shù)組成[9], 分別為艏部阻力系數(shù)C(目標參數(shù))、艏部邊界層轉捩點x(目標參數(shù))、艏部表面噪聲源表面聲功率級S(目標參數(shù))和艏部線型豐滿度(約束參數(shù))。

    系統(tǒng)流程構架示意如圖1所示。

    圖1 優(yōu)化流程框圖

    2.1 艏部線型設計模型

    式中:k1為數(shù)學線型在=0處的曲率半徑;k2為數(shù)學線型在=1處的曲率變化率。

    其與物理線型對應量0和K1的關系為

    函數(shù)()與K1()分別為0與k1的影響函數(shù), 反映了0與k1對線型的影響程度[10]。

    2.2 艏部邊界層計算模型

    UUV艏部表面軸對稱邊界層動量積分方程

    式中:為卡門形狀因子;為動量厚度;C為表面摩擦應力系數(shù)。

    軸對稱邊界層的算法分為層流邊界層和湍流邊界層2種情況。層流邊界層采用Holstein- Bohlen單參數(shù)法進行計算[11], 引入無量綱參數(shù)

    將式(5)代入式(4), 可得

    由式(7)可獲得艏部邊界層局部動量雷諾數(shù)關系表達式

    式中:tr為當?shù)貏恿亢穸壤字Z數(shù), 其表達式

    Re為層流轉捩點處子午線弧長雷諾數(shù),其表達式

    在向下游推進計算中, 當tr首次大于或等于式(8)與右端項時, 表明該點已經轉捩或者正好是轉捩點。

    2.3 艏部流噪聲計算模型

    聲場的控制方程采用的FW-H 聲學模型, 該模型是以Lighthill聲類比理論為基礎, 結合Kirchhof公式和FW-H公式, 其控制方程如下。

    Lighthill密度波動方程[12]

    K-FWH控制方程

    式中:T為Lighthill應力張量, 其表達式

    P為壓應力張量, 其表達式

    ¢為遠場聲壓, 其表達式

    ¢=-0(14)

    ()和()分別為狄拉克函數(shù)和海維賽德函數(shù);為壁面函數(shù);uu分別為來流速度在x方向和垂直壁面方向的分量;vv分別為運動表面響應的速度分量;0為遠場聲速;n為單位法向矢量由固體邊界指向流場。

    式(12)為N-S方程與連續(xù)性方程的重新組合, 可以化為具有2個面源項(單、偶極子)和1個體源項(四極子)的非齊次波動方程形式。由于FW-H方程基于連續(xù)性方程, 因此其分母在非線性區(qū)域仍然滿足波動方程。

    3 優(yōu)化設計算例

    本文的優(yōu)化設計程序通過自適應模擬退火算法實現(xiàn)了低阻低噪一體化優(yōu)化設計。這里選擇總體部分參數(shù)參考SUBOFF數(shù)據, 生成艏部曲線段與部分中段圓柱段, 總長為3.5 m, 頭部曲線段長度1.5 m, 直徑為0.51 m, 速度為6 kn, 設計參數(shù)k1和k2的初始值分別為2。根據工程經驗和結構設計方面的要求, 確定了優(yōu)化的約束條件k1≤3,k2≤10, 以保證在圓弧段與直線段銜接處不會出現(xiàn)明顯的凸起, 且保證一定的豐滿度。優(yōu)化方向是圓弧段與直線段的組合, 設計變量的搜索區(qū)間如圖2所示。

    圖2 設計變量計算域

    本文選取在定常層/湍流轉換情況下的流體動力特性作為目標函數(shù), 其中第1個目標函數(shù)取為在設計狀態(tài)阻力系數(shù)C最小, 第2個目標函數(shù)為在設計狀態(tài)轉捩點位置x盡可能靠后布置, 第3個目標函數(shù)為在設計狀態(tài)下表面噪聲源白噪聲表面聲功率級S最低。約束條件為在滿足外形約束的條件下, 具有一定程度的豐滿度。因此, 該優(yōu)化問題可以描述為, 在優(yōu)化設計點尋找一組設計變量集(k1,k2)[13]。則UUV艏部一體化優(yōu)化準則為

    將第2個目標函數(shù)作最小化處理, 則3個目標函數(shù)的目標值可以獲得[14]

    其中,W為目標函數(shù)權重系數(shù)。

    約束條件

    由式(17)可知, 目標優(yōu)化的搜索方向是沿著使減小的方向實現(xiàn)退火處理, 但由于單個目標值不同, 其對的貢獻亦不同[15]。

    根據式(17)可知, 權重系數(shù)決定了優(yōu)化目標值的大小以及決策方向的轉移, 可以根據不同要求采用適當?shù)臋嘀叵禂?shù)以改變解空間的搜索方向。本文采用3種權重比例工況如表1所示。

    表1 不同權重系數(shù)對比計算工況

    Table 1 Comparison of different weight coefficient calcu- lation conditions

    從優(yōu)化結果可看出, 隨著k1的降低與k2的提高,C,xS均較初始值有所改善。圖3給出了優(yōu)化工況1與優(yōu)化前的UUV艏部外形線型與表面壓力系數(shù)分布對比曲線。由圖3可知, 優(yōu)化后的表面壓力系數(shù)最小值明顯小于優(yōu)化前, 且壓力系數(shù)變化趨于平緩。圖4給出了優(yōu)化前后UUV表面的聲功率級分布,由圖4可知,優(yōu)化后圓弧段聲功率級峰值和均值均小于初始值, 且圓弧段與圓柱段銜接處初始值出現(xiàn)突變, 且峰值較大, 優(yōu)化后銜接處聲功率級變化非常平緩, 有效改善了UUV艏部噪聲源聲功率級。可見, 基于自適應退火模擬法對UUV艏部線型的優(yōu)化是較為成功的, 能夠同時兼顧對C,xS等多目標的優(yōu)化, 可作為UUV外形一體化設計的有效工具。

    圖3 UUV艏部線型及壓力系數(shù)分布對比

    圖4 UUV艏部表面噪聲源對比

    表2給出了優(yōu)化結果和優(yōu)化前初始外形的目標函數(shù)值。由表2可知, 不同權重系數(shù)所構成的凸多邊形對優(yōu)化結果具有較為明顯的影響, 顯著改變了參數(shù)的全域搜索方向與優(yōu)化結果。突出相關權重系數(shù)時, 其對應的優(yōu)化目標亦有相應的重點優(yōu)化提高, 但同時帶來的結果是導致其余相應對象優(yōu)化指標一定程度的下降。但從優(yōu)化結果分析, 各種工況目標參數(shù)指標均較原有線型有明顯的提高。因此, 在對航行器艏部線型進行優(yōu)化設計時, 應根據實際需要選取恰當?shù)膬?yōu)化權重凸多邊形, 從而實現(xiàn)艏部線型最終設計符合期望設計要求。

    表2 參數(shù)優(yōu)化對比

    4 結束語

    本文針對UUV線型變化對流體動力系數(shù)的優(yōu)化問題, 提出了在初始設計過程中將多種目標函數(shù)作為一個整體進行綜合設計的優(yōu)化方法, 能夠從一開始就對相互制約的各種流體動力參數(shù)進行系統(tǒng)化配置, 在目標函數(shù)中綜合考慮了流體動力與聲學兩個雙學科之間的耦合關系, 對魚雷外形進行了優(yōu)化設計, 從而得到最優(yōu)設計。本文基于格蘭維爾線型對UUV艏部線型進行了優(yōu)化仿真計算, 從流體動力系數(shù)/流噪聲特性優(yōu)化前后仿真數(shù)據對比結果表明, 運用自適應模擬退火算法對UUV艏部線型進行優(yōu)化是可行的, 在此基礎上分析了不同優(yōu)化權重凸多邊形條件對優(yōu)化結果的影響。

    [1] 王曉武, 林志民, 崔立軍. 無人潛水器及其動力系統(tǒng)技術發(fā)展現(xiàn)狀及趨勢分析[J]. 艦船科學技術, 2009, 31(8): 31- 34. Wang Xiao-wu, Lin Zhi-min, Chui Li-jun.Analysis of Tech- nology Status and Development Trend for Unmanned Under- water Vehicle and Its Propulsion System[J]. Ship Science and Technology, 2009, 31(8): 31-34.

    [2] 張懷新, 潘雨村. CFD在潛艇外形方案比較中的應用[J]. 船舶力學, 2006, 10(4): 1-8. Zhang Huai-xin, Pan Yu-cun.Application CFD to Compare Submarine Hull Forms[J].Journal of Ship Mechanics, 2006, 10(4): 1-8.

    [3] 張宇文. 魚雷外形設計[M]. 西安:西北工業(yè)大學出版社, 1998.

    [4] 余德海, 宋保維, 楊世興, 等. 魚雷外形多目標多學科綜合優(yōu)化設計方法研究[J]. 兵工學報, 2008, 3(29): 337-340. Yu De-hai, Song Bao-wei, Yang Shi-xing, et al. Investigation of Integrated Multi-disciplinary and Multi-objective Optimi- zation of the Torpedo Shape Design Method[J]. ACTA Armamentirii, 2008, 3(29): 337-340.

    [5] 胡坤, 王樹宗, 徐亦凡. 基于免疫遺傳算法的潛艇水動力系數(shù)優(yōu)化研究[J]. 兵工學報, 2009, 29(12): 1532-1536. Hu Kun, Wang Shu-zong, Xu Yi-fan. Optimization Research on Hydrodynamic Coefficients of Submarine Based on Imm- une Genetic Algorithm[J]. ACTA Aamentirii, 2009, 29(12): 1532-1536.

    [6] 成林, 張文明. 基于模擬退火法的礦用車輛卸料系統(tǒng)協(xié)同優(yōu)化[J]. 中北大學學報, 2009, 3(20): 206-212. Cheng Lin, Zhang Wen-ming, Collaborative Optimization of Mining Truck¢s Hoist System Based on Simulated Annealing Algorithm[J]. Journal of North University of China (Natural Science Edition), 2009, 3(20): 206-212.

    [7] Yukish M, Bennett L, Simpson T W. Requirements on MDO Imposed by the Undersea Vehicle Conceptual Design Probl- em [J]. AIAA 2000-4816: 1-8.

    [8] Barr R A.A Review and Comparison of Ship Maneuvering Methods[J].SNAME Transaction, 1993, 1(10): 609-635.

    [9] Sen D. A Study on Sensitivity of Maneuverability Perfor- mance on the Hydrodynamic Coefficients for Submerged Bo- dies[J].Journal of Ship Research, 2000, 44(3): 186-196.

    [10] 黃景泉, 張宇文. 魚雷流體力學[M]. 西安: 西北工業(yè)大學, 1988.

    [11] Ramajeyathilagam K, Vendhan C P. Bhujanga R V. Non- linear Transient Dynamic Response of Rectangular Plates un- der Shock Loading[J]. Impact Engineering, 2000, 24(10): 999 -1015.

    [12] 盧云濤. 全附體潛艇的流場和流噪聲的數(shù)值研究[D]. 上海: 上海交通大學, 2008. Lu Yun-tao.Numerical Simulation of the Flow-Field and Flow-Noise of Fully Appendage Submarine[D].Shanghai: Shanghai JiaoTong University, 2008.

    [13] Namgoong H, Crossley W A, Lyrintzis A S.Aerodynamic O- ptimization of a Morphing Airfoil Using Energy as an Objective[J]. AIAA 2007, 15(9): 2113-2124.

    [14] 蘇偉, 高正紅, 夏露. 隱身性能約束的多目標氣動外形優(yōu)化設計[J]. 空氣動力學學報, 2006, 1(24): 137-140. Su Wei, Gao Zheng-hong, Xia Lu.Multiobjective Optimiza- tion Design of Aerodynamic Configuration Constrained by Stealth Performance[J],ACTA Aerodynamica Sinica, 2006, 1(24): 137-140.

    [15] Roth B D, Crossley W A.Application of Optimization Tech- niques in the Conceptual Design of Morphing Aircraft[C]// AIAA’S 3rdAnnual Aviation Technology Technology, Inte- gration, and Operations(ATIO) Tech, 2003, 1-11.

    UUV Bow Profile Optimization Design Based on Adaptive Simulated Annealing Algorithm

    PEI Xuan, ZHANG Yu-wen, WANG Ya-dong, YUAN Xu-long

    (College of Marine Engineering, Northwestern Polytechnical University, Xi′an 710072, China)

    Based on the theory of global searching in the multi-objective domain, this paper presents a method of optimizing an unmanned underwater vehicle (UUV) bow profile by adopting the simulated annealing algorithm. This method takes UUV bow fullness, hydrodynamics and flow noise as the goal function. An integrated UUV bow shape optimal design model is established by utilizing the classical viscous flow theory and the computational fluid dynamics (CFD) algorithm. The principle of choosing the complex shape parameters and design variables is described according to the optimization of UUV shape. The optimization design of an UUV bow profile is simulated, and the result shows that with the present method, we can effectively achieve the desired goals of UUV optimization design, and get the optimal solution from the domain composed of various objective functions which are restricted each other. The method greatly improves the acoustic performance and the hydrodynamic characteristic of UUV bow.

    unmanned underwater vehicle (UUV); multi-objective optimization; simulated annealing algorithm; bow profile; integrated optimization design

    TJ630.2

    A

    1673-1948(2011)02-0086-05

    2010-06-01;

    2011-08-01.

    裴譞(1983-), 男, 在讀博士, 研究方向為水下航行器總體設計與彈道計算.

    (責任編輯: 陳 曦)

    猜你喜歡
    聲功率線型模擬退火
    模擬退火遺傳算法在機械臂路徑規(guī)劃中的應用
    測控技術(2018年3期)2018-11-25 09:45:08
    整體道床軌道扣件剛度對鋼軌聲功率特性的影響
    高等級公路幾何線型優(yōu)化設計分析與評價
    江西建材(2018年4期)2018-04-10 12:37:28
    基于模糊自適應模擬退火遺傳算法的配電網故障定位
    核安全1級設備線型支承分析方法研究
    SOA結合模擬退火算法優(yōu)化電容器配置研究
    電源技術(2015年5期)2015-08-22 11:18:24
    自由風扇聲功率級測量方法與測量不確定度
    基于遺傳-模擬退火算法的城市軌道交通快慢車停站方案
    一種新的壓縮機噪聲聲功率的室內測量方法
    一種非均勻線型的互連線能量分布模型
    計算物理(2014年1期)2014-03-11 17:00:58
    男女边摸边吃奶| 纵有疾风起免费观看全集完整版| 一区二区日韩欧美中文字幕 | 人人妻人人爽人人添夜夜欢视频| 黑丝袜美女国产一区| 你懂的网址亚洲精品在线观看| 男的添女的下面高潮视频| 亚洲精品久久久久久婷婷小说| 日韩中字成人| 久久久久久久久大av| 一级片'在线观看视频| 黄片播放在线免费| 春色校园在线视频观看| 狂野欧美激情性bbbbbb| 亚洲不卡免费看| 欧美97在线视频| 日韩在线高清观看一区二区三区| 午夜福利在线观看免费完整高清在| 18在线观看网站| freevideosex欧美| 亚洲精品久久久久久婷婷小说| 亚洲激情五月婷婷啪啪| 一本色道久久久久久精品综合| 久久久精品94久久精品| 99久久精品国产国产毛片| 麻豆精品久久久久久蜜桃| av在线app专区| 国产av一区二区精品久久| 伦理电影大哥的女人| 交换朋友夫妻互换小说| 菩萨蛮人人尽说江南好唐韦庄| 欧美国产精品一级二级三级| 99久久精品国产国产毛片| 亚洲成人一二三区av| 免费少妇av软件| 亚洲在久久综合| 免费观看在线日韩| 视频区图区小说| 91久久精品国产一区二区三区| av卡一久久| 国内精品宾馆在线| 国产日韩欧美亚洲二区| 久久久久国产网址| 久久久午夜欧美精品| 亚洲综合色惰| 亚洲欧洲日产国产| 丝袜美足系列| 中文字幕亚洲精品专区| 国产毛片在线视频| 人成视频在线观看免费观看| 少妇人妻 视频| 免费av中文字幕在线| 美女福利国产在线| 欧美亚洲日本最大视频资源| 国产成人午夜福利电影在线观看| 久久久久国产网址| 久久久久久久久久久久大奶| 七月丁香在线播放| 免费观看无遮挡的男女| 青春草国产在线视频| 极品少妇高潮喷水抽搐| 久久久国产精品麻豆| 满18在线观看网站| 一区在线观看完整版| 日韩三级伦理在线观看| 久久韩国三级中文字幕| 欧美 日韩 精品 国产| 在线观看一区二区三区激情| 少妇高潮的动态图| 国产精品无大码| 国产成人91sexporn| 国产精品一区二区在线不卡| 久久久久网色| 中文字幕人妻熟人妻熟丝袜美| 日日摸夜夜添夜夜添av毛片| 国产永久视频网站| 啦啦啦在线观看免费高清www| 亚洲色图综合在线观看| 妹子高潮喷水视频| 国产成人午夜福利电影在线观看| 欧美精品人与动牲交sv欧美| 大片电影免费在线观看免费| 校园人妻丝袜中文字幕| 亚洲av福利一区| 国产毛片在线视频| 在线精品无人区一区二区三| 久久精品熟女亚洲av麻豆精品| 人妻少妇偷人精品九色| 国产乱来视频区| 99久久人妻综合| 99热6这里只有精品| 国产黄色视频一区二区在线观看| 午夜激情av网站| 韩国av在线不卡| 日韩一区二区三区影片| 性色av一级| 日本黄色片子视频| 97精品久久久久久久久久精品| 免费观看性生交大片5| 中文字幕最新亚洲高清| 26uuu在线亚洲综合色| av免费观看日本| 亚洲美女视频黄频| 成人国语在线视频| 国产男女内射视频| 黑人猛操日本美女一级片| 18禁在线播放成人免费| 久久青草综合色| 最近手机中文字幕大全| 夜夜骑夜夜射夜夜干| 日韩亚洲欧美综合| 一本一本综合久久| 国产免费一区二区三区四区乱码| av电影中文网址| 国产无遮挡羞羞视频在线观看| h视频一区二区三区| 久久国产精品大桥未久av| 如日韩欧美国产精品一区二区三区 | 国产精品人妻久久久久久| 国产精品一区二区在线观看99| 中文字幕人妻熟人妻熟丝袜美| 青青草视频在线视频观看| 视频中文字幕在线观看| 在线观看免费高清a一片| 精品少妇久久久久久888优播| 国产高清有码在线观看视频| 久久久国产一区二区| 国产一区二区在线观看日韩| 国产精品无大码| 18+在线观看网站| 亚洲丝袜综合中文字幕| 午夜影院在线不卡| 国产欧美日韩一区二区三区在线 | 两个人免费观看高清视频| 黑人巨大精品欧美一区二区蜜桃 | 午夜av观看不卡| 亚洲怡红院男人天堂| 桃花免费在线播放| 国产精品人妻久久久影院| 高清av免费在线| 王馨瑶露胸无遮挡在线观看| 欧美 日韩 精品 国产| 人妻一区二区av| 亚洲精品,欧美精品| 亚洲av电影在线观看一区二区三区| 女性被躁到高潮视频| 亚洲人成网站在线播| 国产不卡av网站在线观看| 最近中文字幕2019免费版| 亚洲综合色网址| 免费黄色在线免费观看| 色94色欧美一区二区| 各种免费的搞黄视频| 男男h啪啪无遮挡| 亚洲精品国产av蜜桃| 成年人免费黄色播放视频| 免费人妻精品一区二区三区视频| 国产黄片视频在线免费观看| 18禁动态无遮挡网站| 精品酒店卫生间| 国产熟女午夜一区二区三区 | 国产69精品久久久久777片| av视频免费观看在线观看| kizo精华| 国产黄色免费在线视频| 激情五月婷婷亚洲| 最后的刺客免费高清国语| 欧美成人午夜免费资源| 插逼视频在线观看| 欧美人与性动交α欧美精品济南到 | 91精品三级在线观看| 国产熟女午夜一区二区三区 | 亚洲精品久久成人aⅴ小说 | 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| 国产日韩欧美亚洲二区| 人妻夜夜爽99麻豆av| 97超视频在线观看视频| 亚洲高清免费不卡视频| 久久久久久久大尺度免费视频| 国产毛片在线视频| 美女主播在线视频| 91久久精品国产一区二区成人| 超碰97精品在线观看| 蜜桃在线观看..| a级毛片在线看网站| 九色成人免费人妻av| 成年人午夜在线观看视频| 午夜影院在线不卡| 久久99一区二区三区| 精品国产乱码久久久久久小说| 综合色丁香网| 纯流量卡能插随身wifi吗| 99re6热这里在线精品视频| 日韩av在线免费看完整版不卡| 免费久久久久久久精品成人欧美视频 | 久久婷婷青草| 久久热精品热| 狂野欧美激情性xxxx在线观看| 免费黄色在线免费观看| 亚洲精品美女久久av网站| 亚洲精品日韩av片在线观看| 国产免费视频播放在线视频| 欧美成人精品欧美一级黄| 在线精品无人区一区二区三| 男人操女人黄网站| 青春草亚洲视频在线观看| av女优亚洲男人天堂| 亚洲精品久久成人aⅴ小说 | 成人免费观看视频高清| 少妇人妻久久综合中文| 一区二区三区免费毛片| 高清在线视频一区二区三区| 我的老师免费观看完整版| 婷婷色综合大香蕉| 91精品伊人久久大香线蕉| 亚洲不卡免费看| 人人妻人人澡人人爽人人夜夜| 亚洲欧美成人综合另类久久久| 国产成人免费无遮挡视频| 18禁观看日本| 九九在线视频观看精品| 国产精品成人在线| 啦啦啦在线观看免费高清www| 2021少妇久久久久久久久久久| 精品亚洲乱码少妇综合久久| 夜夜看夜夜爽夜夜摸| 精品久久蜜臀av无| 亚洲,一卡二卡三卡| 亚洲av不卡在线观看| 99re6热这里在线精品视频| 最黄视频免费看| 久久精品久久精品一区二区三区| 99久久精品一区二区三区| 一区二区av电影网| 美女主播在线视频| 日韩欧美精品免费久久| 久久久久国产网址| 美女xxoo啪啪120秒动态图| 丁香六月天网| 久久精品国产自在天天线| 丝袜脚勾引网站| av有码第一页| 免费av中文字幕在线| 亚洲少妇的诱惑av| 五月伊人婷婷丁香| 啦啦啦啦在线视频资源| 精品人妻偷拍中文字幕| 亚洲欧美清纯卡通| 免费av不卡在线播放| 9色porny在线观看| 日韩欧美一区视频在线观看| 777米奇影视久久| 99热这里只有精品一区| 日韩人妻高清精品专区| 五月天丁香电影| 国产免费现黄频在线看| 国产黄片视频在线免费观看| 中国美白少妇内射xxxbb| 国产精品国产三级国产av玫瑰| 婷婷成人精品国产| 国产极品天堂在线| 一本一本综合久久| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久av不卡| 伊人久久精品亚洲午夜| 国产综合精华液| 在线观看免费视频网站a站| 国产成人91sexporn| 成人二区视频| 欧美日韩成人在线一区二区| 97精品久久久久久久久久精品| 免费观看在线日韩| 能在线免费看毛片的网站| 人人妻人人爽人人添夜夜欢视频| 国产高清国产精品国产三级| 久久精品熟女亚洲av麻豆精品| 美女国产视频在线观看| 女的被弄到高潮叫床怎么办| 免费人妻精品一区二区三区视频| 99热6这里只有精品| 日本黄大片高清| 毛片一级片免费看久久久久| 麻豆成人av视频| 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 日本-黄色视频高清免费观看| 国产黄片视频在线免费观看| 国产毛片在线视频| 亚洲av电影在线观看一区二区三区| 青春草视频在线免费观看| 日本wwww免费看| 午夜久久久在线观看| 丁香六月天网| 一级毛片 在线播放| 亚洲国产成人一精品久久久| 五月伊人婷婷丁香| 欧美 日韩 精品 国产| 这个男人来自地球电影免费观看 | 人妻夜夜爽99麻豆av| 精品久久久精品久久久| 久久国产精品大桥未久av| 免费观看性生交大片5| a级片在线免费高清观看视频| 极品人妻少妇av视频| 嘟嘟电影网在线观看| 国产精品国产三级专区第一集| 国产精品嫩草影院av在线观看| 国产成人精品在线电影| 在线观看免费高清a一片| 99re6热这里在线精品视频| 日韩强制内射视频| 国产av一区二区精品久久| 中文精品一卡2卡3卡4更新| 久久精品国产鲁丝片午夜精品| 国产 一区精品| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 欧美日韩国产mv在线观看视频| 久久鲁丝午夜福利片| 国产精品一区二区在线不卡| 久久亚洲国产成人精品v| 一级片'在线观看视频| av在线老鸭窝| 爱豆传媒免费全集在线观看| 丝袜喷水一区| 亚洲av男天堂| 乱人伦中国视频| 99热国产这里只有精品6| 国产69精品久久久久777片| 18在线观看网站| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 亚洲综合精品二区| 国产男女超爽视频在线观看| 亚洲高清免费不卡视频| 99视频精品全部免费 在线| 18禁在线无遮挡免费观看视频| 久久久久国产网址| 卡戴珊不雅视频在线播放| 中国美白少妇内射xxxbb| 狂野欧美激情性xxxx在线观看| 午夜视频国产福利| 国产片内射在线| 秋霞伦理黄片| 免费av中文字幕在线| av专区在线播放| 最近最新中文字幕免费大全7| 久久精品国产亚洲av天美| 最黄视频免费看| 欧美成人精品欧美一级黄| 亚洲熟女精品中文字幕| 亚洲成人一二三区av| 国产精品无大码| 另类精品久久| 亚洲天堂av无毛| 国产亚洲一区二区精品| 日本av免费视频播放| 国产成人aa在线观看| 中文天堂在线官网| 在线免费观看不下载黄p国产| 美女大奶头黄色视频| 久久久久网色| 日韩av在线免费看完整版不卡| 人妻少妇偷人精品九色| 亚洲精品自拍成人| 国产亚洲一区二区精品| 观看av在线不卡| 亚洲国产成人一精品久久久| 亚洲国产av影院在线观看| 乱人伦中国视频| 美女国产视频在线观看| 天天影视国产精品| 夫妻午夜视频| 久久久精品免费免费高清| 久久久精品94久久精品| 亚洲国产精品999| 久久综合国产亚洲精品| 国产成人精品福利久久| 大片免费播放器 马上看| 成人免费观看视频高清| 一本久久精品| 久久99一区二区三区| av女优亚洲男人天堂| 伊人亚洲综合成人网| 国产黄色视频一区二区在线观看| 视频在线观看一区二区三区| 久久精品国产亚洲网站| 69精品国产乱码久久久| 熟妇人妻不卡中文字幕| 一本—道久久a久久精品蜜桃钙片| 国精品久久久久久国模美| 国产视频首页在线观看| 国产欧美日韩一区二区三区在线 | 国产高清有码在线观看视频| 国产又色又爽无遮挡免| 少妇的逼好多水| 国产一区有黄有色的免费视频| 亚洲国产av影院在线观看| 性色av一级| 一级,二级,三级黄色视频| 免费观看无遮挡的男女| 欧美老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频 | 国产黄频视频在线观看| 国产精品久久久久久精品电影小说| 最近中文字幕2019免费版| av天堂久久9| 少妇被粗大的猛进出69影院 | 99国产综合亚洲精品| 黄色怎么调成土黄色| 亚洲三级黄色毛片| 久久精品国产自在天天线| 97精品久久久久久久久久精品| av.在线天堂| 日日摸夜夜添夜夜爱| 亚洲精华国产精华液的使用体验| 99热网站在线观看| 国产精品久久久久久精品古装| 国产免费一级a男人的天堂| 国精品久久久久久国模美| 另类亚洲欧美激情| 亚洲精品色激情综合| 亚洲精品国产av蜜桃| 欧美日韩成人在线一区二区| 热99国产精品久久久久久7| 亚洲经典国产精华液单| 丰满迷人的少妇在线观看| 夜夜爽夜夜爽视频| www.色视频.com| 国产无遮挡羞羞视频在线观看| 91午夜精品亚洲一区二区三区| 一本大道久久a久久精品| 看免费成人av毛片| 男女国产视频网站| 免费少妇av软件| 国产精品不卡视频一区二区| 99热这里只有精品一区| 2022亚洲国产成人精品| 久久国内精品自在自线图片| 人妻制服诱惑在线中文字幕| 99视频精品全部免费 在线| 日日爽夜夜爽网站| 熟女av电影| 一边亲一边摸免费视频| 中文字幕人妻熟人妻熟丝袜美| 少妇丰满av| 晚上一个人看的免费电影| 街头女战士在线观看网站| 午夜激情av网站| 亚洲第一区二区三区不卡| 一区二区三区免费毛片| 国产精品一区二区三区四区免费观看| 制服人妻中文乱码| 一级毛片aaaaaa免费看小| 草草在线视频免费看| 中文字幕制服av| 亚洲精品日韩av片在线观看| 丰满饥渴人妻一区二区三| 天天操日日干夜夜撸| 一本大道久久a久久精品| 黑丝袜美女国产一区| 久久久久网色| 欧美日韩视频精品一区| 十八禁高潮呻吟视频| 又粗又硬又长又爽又黄的视频| 亚洲丝袜综合中文字幕| 啦啦啦中文免费视频观看日本| videosex国产| 成年人午夜在线观看视频| 亚洲精品日韩在线中文字幕| 纯流量卡能插随身wifi吗| 爱豆传媒免费全集在线观看| 免费人成在线观看视频色| 久久人人爽av亚洲精品天堂| 久久女婷五月综合色啪小说| 国产 精品1| 国产一区二区三区av在线| 久久久久久久久久久丰满| 午夜激情久久久久久久| 精品少妇黑人巨大在线播放| 春色校园在线视频观看| 免费人妻精品一区二区三区视频| 久久精品国产a三级三级三级| 一级片'在线观看视频| 麻豆成人av视频| 亚洲国产最新在线播放| 精品人妻偷拍中文字幕| 狠狠精品人妻久久久久久综合| 亚洲精品色激情综合| 熟女av电影| 久久午夜福利片| 日韩人妻高清精品专区| 久久久久久人妻| 五月玫瑰六月丁香| 尾随美女入室| 一本久久精品| 亚洲精品,欧美精品| 国产在线一区二区三区精| 一级a做视频免费观看| 十八禁网站网址无遮挡| 99久国产av精品国产电影| 免费高清在线观看视频在线观看| 亚洲少妇的诱惑av| 狂野欧美激情性xxxx在线观看| 久久久国产精品麻豆| a 毛片基地| 日韩伦理黄色片| 亚洲三级黄色毛片| 91在线精品国自产拍蜜月| 99国产综合亚洲精品| 国产精品熟女久久久久浪| 日韩一本色道免费dvd| 亚洲人成网站在线播| 少妇的逼好多水| 成年人免费黄色播放视频| 一本一本综合久久| 日本色播在线视频| 亚洲欧美色中文字幕在线| 亚洲情色 制服丝袜| 国产男女超爽视频在线观看| 韩国av在线不卡| 国产成人91sexporn| 国产不卡av网站在线观看| 一区二区三区免费毛片| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| 搡老乐熟女国产| 波野结衣二区三区在线| 狂野欧美激情性bbbbbb| 日本wwww免费看| 午夜免费男女啪啪视频观看| 国模一区二区三区四区视频| 欧美日韩成人在线一区二区| 久久久国产欧美日韩av| 视频中文字幕在线观看| 亚洲av日韩在线播放| 特大巨黑吊av在线直播| 一级a做视频免费观看| 女性生殖器流出的白浆| 久久久久国产精品人妻一区二区| 国产一区二区在线观看av| 蜜桃在线观看..| 国语对白做爰xxxⅹ性视频网站| 免费人妻精品一区二区三区视频| 国产高清三级在线| 日本爱情动作片www.在线观看| 国产免费福利视频在线观看| 亚洲国产成人一精品久久久| √禁漫天堂资源中文www| 性高湖久久久久久久久免费观看| 狂野欧美白嫩少妇大欣赏| 一区二区av电影网| 日韩在线高清观看一区二区三区| 一本大道久久a久久精品| 亚洲色图 男人天堂 中文字幕 | 久久精品国产a三级三级三级| 日韩强制内射视频| 在线天堂最新版资源| 婷婷色av中文字幕| 国产男人的电影天堂91| 亚洲中文av在线| 欧美xxxx性猛交bbbb| 日本wwww免费看| 熟女av电影| 亚洲四区av| 观看av在线不卡| 精品人妻在线不人妻| av卡一久久| 国产免费福利视频在线观看| 国产免费一区二区三区四区乱码| 亚洲精品国产av蜜桃| 精品一区二区免费观看| av在线app专区| 在线观看三级黄色| 少妇被粗大的猛进出69影院 | 少妇人妻 视频| 国产淫语在线视频| 久久综合国产亚洲精品| 欧美最新免费一区二区三区| a级毛片黄视频| 国产精品成人在线| 中文字幕人妻熟人妻熟丝袜美| 欧美国产精品一级二级三级| 在线精品无人区一区二区三| 99久久综合免费| 在线看a的网站| 麻豆精品久久久久久蜜桃| 欧美人与性动交α欧美精品济南到 | 久久久亚洲精品成人影院| 久久精品国产亚洲av天美| 美女xxoo啪啪120秒动态图| 国产免费一级a男人的天堂| 亚洲不卡免费看| 黄片播放在线免费| 成人毛片a级毛片在线播放| 日日啪夜夜爽| 观看美女的网站| 久久热精品热| 丝袜喷水一区| 在线观看免费高清a一片| 一区二区三区精品91| 久久精品熟女亚洲av麻豆精品| 国模一区二区三区四区视频| 日本色播在线视频| 考比视频在线观看| 国产亚洲精品久久久com| 3wmmmm亚洲av在线观看| 久久久精品免费免费高清| 高清欧美精品videossex| 视频区图区小说| 人人妻人人爽人人添夜夜欢视频| a 毛片基地|