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

    Ship Manoeuvring Performance and Construction of its Metamodels

    2015-05-02 09:29:14LIDongqinPhilipWILSONJIANGZhiyong
    船舶力學(xué) 2015年12期
    關(guān)鍵詞:南安普頓船型工程學(xué)院

    LI Dong-qin,Philip A.WILSON,JIANG Zhi-yong

    (1 School of Naval Architecture&Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang, China,212003;2 Fluid Structure Interactions Research Group,Faculty of Engineering and the Environment, University of Southampton,Southampton,UK,SO17 1BJ)

    Ship Manoeuvring Performance and Construction of its Metamodels

    LI Dong-qin1,Philip A.WILSON2,JIANG Zhi-yong1

    (1 School of Naval Architecture&Ocean Engineering,Jiangsu University of Science and Technology,Zhenjiang, China,212003;2 Fluid Structure Interactions Research Group,Faculty of Engineering and the Environment, University of Southampton,Southampton,UK,SO17 1BJ)

    To improve the efficiency of ship manoeuvring calculation technique,NAPA-based calculations and Support Vector Regression(SVR)were combined to predict the manoeuvring performance of the Offshore Support Vessels(OSV),and enough ship information about the offshore support vessel were gathered;a series of 30 similar hull bodies was adopted by the Latin Hypercube Design which were employed to explore the design space and to sample data for covering the design space.The ship hull series were generated from the affine and displacement transformation and some adjustments in NAPA according to the ship design variables,thus creating the calculation model for the hull geometry expression.For each of the ship models,5 different manoeuvrability criteria were calculated,which were the advance,tactical diameter,transfer,10°/10°first overshoot angle and 20° /20°first overshoot angle.To improve the efficiency of manoeuvring calculation,the Single-parameter Lagrangian Support Vector Regression(SPL-SVR)was adopted and trained to establish the metamodels and predict the manoeuvring performance and this new algorithm was first proposed by the author and combined with Laplace loss function,which has only one parameter to control the errors and adds b2/2 to the item of confidence interval at the same time,For the OSV case,the manoeuvrability criteria were predicted with the SPL-SVR and compared with the NAPA-based calculation results with manoeuvring manager,the Artificial Neural Network results and classical SVR results. The results agreed well with each other.Instead of requiring the evaluation of expensive simulation codes,the metamodels of ship manoeuvring performance were suitable for the practical application in ship preliminary design stage and all the numerical results show the effectiveness and practicability of the new approximation algorithms.

    metamodel;Support Vector Machine;ship manoeuvring

    0 Introduction

    In real world ship design problems,there is often more than one area of element in the overall design;that is,there may be different disciplines that contribute to the ship design[1-2],such as structures,economics,resistance,manoeuvring and seakeeping performance,and so on. How to evaluate these ship performances efficiently is still a research hotspot[3-4].Now manoeuvring performance is considered in the preliminary ship design stage because ship manoeuvrability is directly related to navigation safety and economy.Ship manoeuvrability is the ability for a ship to keep or change motion state under the control actions,i.e.,to keep the straight-ahead course with constant speed,or to change the speed,the course and or the position of the ship,according to the intention of the helmsman.

    Each ship will has its own manoeuvring characteristics.The position of the pivot point will vary performance,while performance itself can be affected by numerous factors.The development of new ship design technology is dependent upon a cooperative,multidisciplinary design approach.To reduce the computational cost of implementing computer-based simulations and analyses in ship design,a variety of metamodeling techniques have been developed[5-7],for instance Response Surface Model(RSM),kriging and ANN.Metamodel is a key element of the Multidisciplinary Design Optimization(MDO).In this paper,a new simple and effective algorithm of Support Vector Machines is proposed and used to establish the alternative metamodels of ship manoeuvring to provide the approximation of a design space.

    1 Mathematical model of SPL-SVR

    Accordingly,it is really important to improve the accuracy and robustness performance of metamodeling techniques especially when the sample size becomes small,limited and scarce. The Support Vector Machines(SVM)[8-10]aims at the limited samples and has a good generalization performance as well as global optimal extremum which have been proved by many researchers[11].In this paper,we will use a new support vector regression algorithm SPL-SVR which is proposed by authors in Ref.[12]to construct the implicit metamodels of ship manoeuvring performance as shown in red frame of Fig.1.The detailed description of this algorithm and its applications can be found in the author’s previous work[12].Here,we will recall the mathematical theory of this algorithm for the readers’convenience and it is named Single-parameter Lagrangian Support Vector Regression(SPL-SVR).The formula is listed as follows:

    The above optimization problem can be stated as standard formulized quadratic programming.The dual problem can then be expressed as the following standard quadratic programming form.

    Thus,the estimation function is calculated as follows:

    With this simpler algorithm,we can obtain the black box which describes the complicated mapping relation without knowing the connection between the dependent variables and independent variables.Therefore,this new algorithm is well suited for the construction of ship manoeuvring approximation model in the preliminary and early-stage design.

    2 Establishment of metamodels of ship manoeuvring performance

    In this work,the Single-parameter Lagrangian Support Vector Regression is used to establish the metamodels of ship manoeuvring performance,which can reduce the huge calculation cost of the simulation codes and satisfy the growing needs of computation in the multidisciplinary ship design optimization.

    2.1 Distribution of ship samples

    Before constructing the metamodels of ship manoeuvring performance,we need to gatherplenty of ship information to decide the scope of design variables and select the training ship data set.At the same time,plenty of data about the offshore support vessels were gathered from the shipping companies and design institutions.The distributions of main principal characteristics are showed in Fig.1,in which the red points represent the 30 training ship data from DOE method.

    2.2 Design of experiments

    Latin Hypercube Designs[13]is adopted here as design of experiments(DOE)to sample out the points of the ship design variables.The Latin Hypercube method chooses points to maximize distance between design points,but with a constraint which maintains the even spacing between factor levels.

    As is well known,there are many variables which affect the ship manoeuvring performance.We chose the length between perpendiculars,breadth,depth,design draught,the longitudinal centre of buoyancy,ship velocity and diameter of propeller as the design variables, and these seven parameters can show the geometrical characteristics of ship hull.Also,the standard Model-based calibration toolbox from commercial software Matlab is used to choose 30 testing data set with Latin Hypercube Design and the space distribution of these ship types is showed in Fig.2.Once the fixed parameters are established,the design variables are chosen and we just select 30 testing ships data here listing in the Tab.1.

    Tab.1 Design variables of selected 15 training ship data

    Continue Tab.1

    Fig.2 The space distribution of ship training data set

    2.3 Establishment of metamodels of ship manoeuvring performance

    The ship hulls of these thirty ship samples were generated from the parent ship by thesoftware NAPA within several steps.Firstly,construct the hull of the parent ship in NAPA. Secondly,transform the parent frame area distribution into the needed one according to the difference of the longitudinal coordinates of buoyancy between the parent ship and the design ship.Thirdly,combine the affine and displacement transformation into the design values Lpp, B,D,T under the task TRA in NAPA.Finally,make some adjustments of the ship hull in local place from the personal experience.The ship hull of one training ship is shown in Fig.3.

    Fig.3 Transverse section and 3D lay-out of ship hull

    Before establishing the metamodels of manoeuvring performance in Multidisciplinary Ship Design Optimization,the ship manoeuvring performance of OSV will be calculated at first. To develop a practical approximation models of ship manoeuvring performance in the hydrodynamic-based multidisciplinary design optimization of offshore support vessel at the early design stage,a practical calculation tool,based on the MMG(Ship Manoeuvring Mathematical Model Group)Model called Manoeuvring Manager from the commercial software NAPA,is used to compute the manoeuvring criteria.

    The coordinate systems used in the manoeuvring calculations are as follows:for all input data of manoeuvring devices and hydrostatics,etc,the normal NAPA coordinates are used with X axis starting from the aft perpendicular and the positive Y coordinate to port for the default right handed coordinates;maneuvering derivatives are defined in the manoeuvring coordinate system fixed to the midship;all results for location,velocities,etc,are presented in the manoeuvring coordinate system fixed to the centre of gravity.These coordinate systems and rotations in NAPA are shown in Fig.4.The simulation of turning circle manoeuvre is shown in Fig.5; the simulations of Zigzag tests for different turning angles are shown in Fig.6 and Fig.7.

    Fig.4 System of coordinates in NAPA

    Fig.5 The simulation of turning circle manoeuvre

    The International Maritime Organization(IMO)proposed criteria for parameters derived from the standard manoeuvres.These criteria,described in IMO Resolution A.751(18)(1993)[14],are commonly used to judge the manoeuvring characteristics of a vessel.Here we choose advance,tactical diameter,transfer,10°/10°first overshoot angle and 20°/20°first overshoot angle as manoeuvring criteria to evaluate the performance for the offshore support vessel.These calculated manoeuvring criteria of 30 selected ship types are listed in Tab.2.

    Fig.6 The simulation of Zigzag test 10°/10°

    Fig.7 The simulation of Zigzag test 20°/20°

    Tab.2 Calculated manoeuvring criteria of selected 30 training ship data

    Here,the first 20 ship types created with DOE method are selected as training data set and the last 10 ship types as test data set.The chosen variables are the length between perpendiculars,breadth,depth,design draught,the longitudinal centre of buoyancy,ship velocity and diameter of propeller as the design variables,and the ship advance,tactical diameter,transfer,10°/10°first overshoot angle and 20°/20°first overshoot angle as each output variable.The calculation results were compared with Manoeuvring Manager,ANN and classic SVR which were shown as Fig.8.Here the calculation results for advance and 10°/10°first overshoot angle are just listed in Tabs.3 and 4 for instance.The Relative Error(RE)and Mean relative error(MRE)are applied as performance indexes:

    Tab.3 Calculation results with Relative Error(RE)for advance

    Tab.4 Calculation results with Relative Error(RE)for 10°/10°first overshoot angle

    Fig.8 Approximation results of manoeuvring criteria for ship type 20 to 30

    Considering the above circumstances,the proposed SPL-SVR algorithm shows a good approximation and prediction performance.The comparison of relative errors is listed in Tab.5. The maximum MRE within the five manoeuvring criteria comparing to the result of Manoeuvring Manager for ANN is 7.77%and the minimum MRE is 5.97%;the maximum MRE for SVR is 4.89%and the minimum MRE is 3.41%;the maximum MRE for SPL-SVR is 3.15%and the minimum MRE is 1.69%.It can be seen that the results are acceptable and agree well with each other.Obviously,if the training ships data set,the kernel parameters and the calculation method for manoeuvring are chosen properly,we can use these metamodels to predict the ship manoeuvring performance in the preliminary ship design stage.

    Tab.5 The errors comparison for the five manoeuvring criteria

    Using these calculated values of manoeuvring criteria,it is possible to establish metamodels of manoeuvring performance of offshore support vessels in the process of multidisciplinary design optimization.The programs are written in Matlab which are easy to integrate into the Optimus software and the metamodels of ship manoeuvring performance are constructed basing on the proposed Support Vector Regression algorithm.Without running expensive model tests or time consuming CFD calculations,the benchmarking methodology presented in this paper can give an insight of the relationship between input design variables and output responses,facilitate the integration of discipline dependent analysis codes and also reduce the time consumption of computer simulation and saving time in the preliminary ship design stage.

    3 Conclusions

    The assessment of ship manoeuvring performance is studied in this paper including advance,tactical diameter,transfer,10°/10°first overshoot angle and 20°/20°first overshoot angle.The metamodels for manoeuvring performance established by a new SPL-SVR algorithm in conjunction with LHS are employed in place of expensive simulation and analysis codes and these metamodels can be used to evaluate the ship manoeuvring performance efficiently at preliminary design stages of offshore support vessel.Without using computationally expensive methods such as CFD or model tests,the main advantage of this methodology is that it can provide detailed and realistic operational profiles of ship designs at an early stage of the design process.

    As part of the future work,metamodels of ship resistance,seakeeping and manoeuvring can be combined in the framework of Multidisciplinary Design Optimization to improve its convergence efficiency.Multidisciplinary and multiobjective optimization design problems widely exist in the field of ship design,development of effective optimization framework for multidisciplinary and multi-objective problems can be also considered as a future research direc-tion.Additionally,model tests and high-precision simulations of manoeuvring performance will be undertaken in case to enhance the calculation precision of those metamodels for offshore support vessels.

    Acknowledgments

    The authors wish to thank The Natural Science Foundation of Jiangsu Province of China (Grant No.BK2012696)and the Project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)for their financial support.The author also would like to thank Prof.Philip A.Wilson for his constructive and insightful suggestions and comments.

    [1]Campana Emilio Fortunato,Fasano Giovanni,Peri Daniele.Penalty function approaches for ship multidisciplinary design optimisation(MDO)[J].European Journal of Industrial Engineering,2012,6(6):765-784.

    [2]Cui Hao,Osman T,Philip S.Learning-based ship design optimization approach[J].CAD Computer Aided Design,2012, 44(3):186-195.

    [3]Woodward M D.Evaluation of inter-facility uncertainty for ship manoeuvring performance prediction[J].Ocean Engineering,2014,88:598-606.

    [4]Broglia R,Dubbioso G,Durante D,Di Mascio A.Simulation of turning circle by CFD:Analysis of different propeller models and their effect on manoeuvring prediction[J].Applied Ocean Research,2012,39:1-10.

    [5]He Jim,Hannapel Shari,Vlahopoulos Nickolas.Multidisciplinary design optimization of ship hull forms using metamodels [C].Proceedings of the ASME Design Engineering Technical Conference,2011,5:847-856.

    [6]Xiao Mi,Gao Liang,Shao Xinyu,Qiu Haobo,Jiang Ping.A generalised collaborative optimisation method and its combination with kriging metamodels for engineering design[J].Journal of Engineering Design,2012,23(5):379-399.

    [7]Hesham G,Chu Xuezheng,Gao Liang,Li Peigen.An approach combined response surface method and particle swarm optimization to ship multidisciplinary design and optimization[C].IEEE International Conference on Industrial Engineering and Engineering Management,2009:1810-1814.

    [8]He Wei,Zou Zaojian.A Lease Square Support Vector Machine Metamodel for ship performance in uncertainty quantification study[J].Journal of Ship Mechanics,2013,17(6):604-615.

    [9]Siamak M,Tillmann F,Suykens Johan A K.Approximate solutions to ordinary differential equations using least squares support vector machines[J].IEEE Transactions on Neural Networks and Learning Systems,2012,23(9):1356-1367.

    [10]Vapnik V N.Universal learning technology:support vector machines[J].NEC Journal of Advanced Technology,2005,2 (2):137-144.

    [11]Liu Changde,Zhang Hua,Han Yang,Shi Chao.Black-box modeling and prediction of ship maneuverability based on least square support vector machine[J].Journal of Ship Mechanics,2013,17(8):872-877.

    [12]Li Dongqin,Guan Yifeng,Wang Qingfeng,Chen Zhitong.Support vector regression-based multidisciplinary design optimization for ship design[C]//Proceedings of the 31th International Conference on Ocean,Offshore and Arctic Engineering,2012,June 10-15.Rio de Janeiro,Brazil,2012.

    [13]Fabian Fuerle,Johann Sienz.Decomposed surrogate based optimization of carbon-fiber bicycle frames using Optimum Latin Hypercubes for constraineddesign spaces[J].Computers&Structures,2013,119:48-59.

    [14]IMO Resolution A.751(18)[M].Interim Standards for Ship Manoeuvrability,1993.

    船舶操縱性能及其代理模型構(gòu)建

    李冬琴1,Philip A.WILSON2,蔣志勇1
    (1江蘇科技大學(xué)船舶與海洋工程學(xué)院,江蘇鎮(zhèn)江,212003;2南安普頓大學(xué)工程與環(huán)境學(xué)院,流體結(jié)構(gòu)耦合研究組,英國(guó)南安普頓,SO17 1BJ)

    為了提高船舶操縱性能計(jì)算效率,文章結(jié)合基于NAPA計(jì)算仿真結(jié)果和支持向量回歸方法預(yù)報(bào)海洋平臺(tái)支援船的操縱性能,在收集足夠相關(guān)船型信息前提下,采用能夠合理探索設(shè)計(jì)空間和抽樣的拉丁超立方方法獲取了30條樣本船型數(shù)據(jù)。通過NAPA仿射變化、位移轉(zhuǎn)換及根據(jù)船型設(shè)計(jì)變量進(jìn)行局部調(diào)整從而生成系列船型以表達(dá)船體幾何形狀。針對(duì)每個(gè)船型,分別計(jì)算了5項(xiàng)操縱性衡準(zhǔn)指標(biāo):進(jìn)距、戰(zhàn)術(shù)直徑、橫距、10°舵角第一超越角和20°舵角第一超越角。為提高船舶操縱性能的計(jì)算效率,文中利用作者早先新提出的一種單參數(shù)Lagrangian支持向量回歸算法來訓(xùn)練并構(gòu)建代理模型以預(yù)報(bào)船舶操縱性能,該算法整合了Laplace損失函數(shù),僅采用單參數(shù)控制計(jì)算誤差并于置信區(qū)間中增加了b2/2項(xiàng)。以海洋平臺(tái)支援船為例,采用SPL-SVR算法預(yù)報(bào)船舶操縱性能,并與基于NAPA計(jì)算仿真結(jié)果、人工神經(jīng)網(wǎng)絡(luò)、經(jīng)典支持向量回歸算法進(jìn)行對(duì)比。不需要昂貴的仿真代碼計(jì)算,文中采用SPL-SVR算法建立的船舶操縱性能響應(yīng)面模型比較適合船型初步設(shè)計(jì)的工程實(shí)際應(yīng)用,并具有較好的效率和適用性。

    代理模型;支持向量機(jī);船舶操縱性

    U661.33

    :A

    李冬琴(1979-),女,江蘇科技大學(xué)船舶與海洋工程學(xué)院副教授;Philip A.WILSON(1954-),男,英國(guó)南安普頓大學(xué)流體結(jié)構(gòu)耦合研究小組教授;

    1007-7294(2015)12-1463-12

    U661.33

    :A

    10.3969/j.issn.1007-7294.2015.12.003

    蔣志勇(1956-),男,江蘇科技大學(xué)船舶與海洋工程學(xué)院教授。

    Received date:2015-05-22

    Foundation item:Supported by the National Natural Science Foundation of China(Grant No.51509114); the Natural Science Foundation of Jiangsu Province of China(Grant No.BK2012696)

    Biography:LI Dong-qin(1979-),female,associate professor,E-mail:mandy_ldq@163.com; Philip A.WILSON(1954-),male,professor.

    猜你喜歡
    南安普頓船型工程學(xué)院
    福建工程學(xué)院
    福建工程學(xué)院
    Discuss Learner-related Factors and How These Might Influence the Process and Outcomes of Classroom Language Learning
    福建工程學(xué)院
    福建工程學(xué)院
    中國(guó)商人收購南安普頓
    中國(guó)古代著名的海船船型
    ——福船
    中國(guó)古代著名的海船船型—廣船
    Anti-ageing effects of a new Dimethylaminoethanol-based formulation on DGalactose induced skin ageing model of rat
    Rapid Prototyping Technology of Tissue Engineering Scaffold
    科技視界(2014年21期)2014-08-21 02:38:06
    videos熟女内射| 美女国产高潮福利片在线看| 美女主播在线视频| 免费高清在线观看日韩| 在线观看免费视频网站a站| 国产无遮挡羞羞视频在线观看| 一二三四中文在线观看免费高清| 黑人高潮一二区| 午夜福利在线观看免费完整高清在| 久久久久久久国产电影| 久久韩国三级中文字幕| 亚洲av.av天堂| 伊人亚洲综合成人网| 久久亚洲国产成人精品v| 欧美成人精品欧美一级黄| 国产日韩欧美视频二区| 亚洲四区av| 美女内射精品一级片tv| 丰满迷人的少妇在线观看| 亚洲伊人久久精品综合| 人成视频在线观看免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产不卡av网站在线观看| 国产乱人偷精品视频| 国产亚洲一区二区精品| 久久国内精品自在自线图片| 久久热在线av| 国产一区亚洲一区在线观看| 久久99精品国语久久久| 欧美xxⅹ黑人| 国产av码专区亚洲av| 九草在线视频观看| 亚洲欧美成人精品一区二区| 精品一区二区三卡| 黑人巨大精品欧美一区二区蜜桃 | 大陆偷拍与自拍| 国产精品不卡视频一区二区| 久久久国产欧美日韩av| 高清毛片免费看| 国产精品久久久久久精品古装| 久久精品夜色国产| 亚洲精品国产av成人精品| 天美传媒精品一区二区| av天堂久久9| 亚洲国产欧美在线一区| 亚洲欧洲国产日韩| 两性夫妻黄色片 | 哪个播放器可以免费观看大片| 日韩人妻精品一区2区三区| 香蕉丝袜av| 伦理电影大哥的女人| 亚洲成av片中文字幕在线观看 | 狂野欧美激情性bbbbbb| 中文乱码字字幕精品一区二区三区| 国产在视频线精品| 久久午夜福利片| 王馨瑶露胸无遮挡在线观看| 亚洲 欧美一区二区三区| 欧美性感艳星| 国产亚洲精品第一综合不卡 | 成人免费观看视频高清| 日韩一区二区视频免费看| 少妇精品久久久久久久| 精品99又大又爽又粗少妇毛片| 欧美成人午夜精品| 亚洲欧美一区二区三区国产| 五月天丁香电影| 卡戴珊不雅视频在线播放| 精品第一国产精品| 久久综合国产亚洲精品| 中国国产av一级| 精品久久国产蜜桃| 搡老乐熟女国产| 国产视频首页在线观看| 国产高清国产精品国产三级| 男女边吃奶边做爰视频| 夫妻午夜视频| 亚洲精品国产色婷婷电影| www.熟女人妻精品国产 | 久久久久久久亚洲中文字幕| 日本91视频免费播放| 亚洲少妇的诱惑av| 自线自在国产av| 黄色视频在线播放观看不卡| 婷婷色综合www| 亚洲国产欧美在线一区| 又粗又硬又长又爽又黄的视频| 亚洲欧美成人综合另类久久久| av免费观看日本| 99香蕉大伊视频| 欧美人与性动交α欧美精品济南到 | 亚洲天堂av无毛| 黄片无遮挡物在线观看| 久久99精品国语久久久| 观看av在线不卡| 99香蕉大伊视频| 99热国产这里只有精品6| 亚洲国产av新网站| 少妇人妻久久综合中文| 少妇的丰满在线观看| 女人精品久久久久毛片| 精品一区二区三区四区五区乱码 | 赤兔流量卡办理| 国产午夜精品一二区理论片| 秋霞伦理黄片| 18+在线观看网站| 高清视频免费观看一区二区| 人妻一区二区av| 亚洲在久久综合| 亚洲精品美女久久av网站| 十八禁网站网址无遮挡| 亚洲,欧美,日韩| 黄色配什么色好看| 黄色配什么色好看| 免费av不卡在线播放| 大陆偷拍与自拍| 人体艺术视频欧美日本| 中文字幕另类日韩欧美亚洲嫩草| 国产女主播在线喷水免费视频网站| 寂寞人妻少妇视频99o| a级毛片在线看网站| 男女边摸边吃奶| 香蕉国产在线看| 久久久久精品久久久久真实原创| 久久久久精品久久久久真实原创| 七月丁香在线播放| 99热网站在线观看| 一本色道久久久久久精品综合| 午夜久久久在线观看| 狂野欧美激情性xxxx在线观看| 在线观看免费高清a一片| 日韩欧美精品免费久久| 美国免费a级毛片| 亚洲国产色片| 蜜桃国产av成人99| 国产乱人偷精品视频| 丁香六月天网| 在线观看免费视频网站a站| 热99久久久久精品小说推荐| 欧美变态另类bdsm刘玥| 亚洲成av片中文字幕在线观看 | 高清毛片免费看| 国产成人精品一,二区| 国产不卡av网站在线观看| 国产毛片在线视频| 久久久国产精品麻豆| 成年美女黄网站色视频大全免费| 一级片免费观看大全| 最近中文字幕2019免费版| 色视频在线一区二区三区| 最近的中文字幕免费完整| 午夜精品国产一区二区电影| 尾随美女入室| 国产在线一区二区三区精| 亚洲美女搞黄在线观看| 亚洲国产毛片av蜜桃av| 精品一品国产午夜福利视频| 少妇被粗大的猛进出69影院 | 九九在线视频观看精品| 制服丝袜香蕉在线| 日韩av不卡免费在线播放| 久久99蜜桃精品久久| 观看av在线不卡| 99热网站在线观看| 久久久久视频综合| 男女边吃奶边做爰视频| 亚洲情色 制服丝袜| 久久午夜福利片| 欧美日韩av久久| 国产精品久久久久久精品古装| 内地一区二区视频在线| 精品一区二区三卡| 亚洲性久久影院| 99国产综合亚洲精品| 超色免费av| 久久久久视频综合| 国产精品欧美亚洲77777| 丰满少妇做爰视频| 久久精品国产综合久久久 | 在线精品无人区一区二区三| 人妻一区二区av| 免费av不卡在线播放| 午夜福利网站1000一区二区三区| 两性夫妻黄色片 | 国产高清国产精品国产三级| 啦啦啦视频在线资源免费观看| 女性被躁到高潮视频| 亚洲精品国产av蜜桃| 国产日韩欧美在线精品| 日韩成人伦理影院| 制服丝袜香蕉在线| 看非洲黑人一级黄片| 又粗又硬又长又爽又黄的视频| 天堂8中文在线网| 亚洲精品国产av成人精品| 麻豆乱淫一区二区| 国产日韩欧美亚洲二区| 久久亚洲国产成人精品v| 全区人妻精品视频| 国产精品无大码| 日韩三级伦理在线观看| 人人妻人人添人人爽欧美一区卜| 大香蕉97超碰在线| 女的被弄到高潮叫床怎么办| 亚洲图色成人| 亚洲精品一二三| 大陆偷拍与自拍| 2022亚洲国产成人精品| 极品人妻少妇av视频| 色哟哟·www| 王馨瑶露胸无遮挡在线观看| 韩国av在线不卡| 欧美精品高潮呻吟av久久| 伦精品一区二区三区| 国产成人精品久久久久久| 精品国产乱码久久久久久小说| 又粗又硬又长又爽又黄的视频| 亚洲精品av麻豆狂野| 9热在线视频观看99| av线在线观看网站| 91精品三级在线观看| 一级,二级,三级黄色视频| 中文精品一卡2卡3卡4更新| 人人澡人人妻人| 国产精品一区www在线观看| 日日啪夜夜爽| 大香蕉久久网| 18在线观看网站| 色吧在线观看| 国产精品欧美亚洲77777| 王馨瑶露胸无遮挡在线观看| 免费看av在线观看网站| 夜夜爽夜夜爽视频| 日韩熟女老妇一区二区性免费视频| 永久免费av网站大全| 国产麻豆69| 亚洲 欧美一区二区三区| 中文乱码字字幕精品一区二区三区| 伊人久久国产一区二区| 国产毛片在线视频| 视频在线观看一区二区三区| 我的女老师完整版在线观看| 日韩成人av中文字幕在线观看| 曰老女人黄片| 亚洲欧美色中文字幕在线| 女人久久www免费人成看片| 久久99精品国语久久久| 精品久久国产蜜桃| 欧美 亚洲 国产 日韩一| 日韩制服丝袜自拍偷拍| 国产伦理片在线播放av一区| 国产男女超爽视频在线观看| 91国产中文字幕| 色哟哟·www| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久久久免| 秋霞伦理黄片| 精品熟女少妇av免费看| 亚洲美女视频黄频| 青青草视频在线视频观看| 九色亚洲精品在线播放| 亚洲图色成人| 成人免费观看视频高清| 最近中文字幕2019免费版| 成年美女黄网站色视频大全免费| 欧美人与性动交α欧美精品济南到 | 久久国产亚洲av麻豆专区| 99热6这里只有精品| 天天影视国产精品| 老司机影院毛片| 亚洲欧洲国产日韩| 国产成人av激情在线播放| 久久99热6这里只有精品| 午夜福利乱码中文字幕| av.在线天堂| 色婷婷av一区二区三区视频| 日日爽夜夜爽网站| 国产无遮挡羞羞视频在线观看| 成人综合一区亚洲| 亚洲人成77777在线视频| 免费播放大片免费观看视频在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲综合精品二区| 亚洲图色成人| 亚洲精品美女久久av网站| 亚洲精品久久久久久婷婷小说| 亚洲国产欧美在线一区| 高清欧美精品videossex| 高清在线视频一区二区三区| 尾随美女入室| 午夜日本视频在线| 日本av手机在线免费观看| 国产精品久久久久久久久免| 国产成人午夜福利电影在线观看| 亚洲精品aⅴ在线观看| 国产成人免费观看mmmm| 一二三四中文在线观看免费高清| 香蕉精品网在线| 另类亚洲欧美激情| 高清黄色对白视频在线免费看| 日韩大片免费观看网站| 一级毛片我不卡| 十八禁高潮呻吟视频| 岛国毛片在线播放| 99热国产这里只有精品6| 亚洲精品av麻豆狂野| 国产一级毛片在线| 国产精品人妻久久久久久| 91午夜精品亚洲一区二区三区| 国产男女内射视频| 国产精品一区www在线观看| 免费少妇av软件| 熟女电影av网| 国产国语露脸激情在线看| 综合色丁香网| 亚洲精品中文字幕在线视频| 在线精品无人区一区二区三| av免费观看日本| av在线播放精品| 美女视频免费永久观看网站| 在线精品无人区一区二区三| 欧美日韩视频精品一区| 亚洲情色 制服丝袜| 亚洲国产欧美日韩在线播放| 在线精品无人区一区二区三| 精品国产乱码久久久久久小说| 久久久久久久精品精品| 蜜桃国产av成人99| 国产成人精品福利久久| 99香蕉大伊视频| 少妇猛男粗大的猛烈进出视频| 久久人人97超碰香蕉20202| 国产白丝娇喘喷水9色精品| 国产精品一国产av| 亚洲内射少妇av| 少妇人妻精品综合一区二区| 欧美老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| 美女福利国产在线| 国产爽快片一区二区三区| 免费大片黄手机在线观看| 日韩不卡一区二区三区视频在线| 黄色毛片三级朝国网站| 飞空精品影院首页| 亚洲经典国产精华液单| 亚洲成国产人片在线观看| 国产亚洲精品久久久com| 国产精品麻豆人妻色哟哟久久| 久久99一区二区三区| 黑人高潮一二区| 美女国产高潮福利片在线看| 久久久精品免费免费高清| 1024视频免费在线观看| 亚洲欧洲精品一区二区精品久久久 | 人妻人人澡人人爽人人| 乱码一卡2卡4卡精品| 日本黄大片高清| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区国产| 国产又色又爽无遮挡免| 国产一区二区三区av在线| av在线播放精品| 国产日韩欧美亚洲二区| 色网站视频免费| xxx大片免费视频| 久久精品夜色国产| 精品少妇黑人巨大在线播放| 又大又黄又爽视频免费| 精品国产国语对白av| 又黄又爽又刺激的免费视频.| 国产福利在线免费观看视频| 精品一区在线观看国产| 免费黄色在线免费观看| 日韩成人伦理影院| 日韩伦理黄色片| 亚洲婷婷狠狠爱综合网| 精品亚洲成a人片在线观看| 亚洲精品aⅴ在线观看| 搡老乐熟女国产| 婷婷色av中文字幕| 免费高清在线观看日韩| 又粗又硬又长又爽又黄的视频| 国产精品99久久99久久久不卡 | 国产精品三级大全| 久久精品夜色国产| 久久久久久久久久成人| 七月丁香在线播放| 水蜜桃什么品种好| 精品少妇黑人巨大在线播放| 看非洲黑人一级黄片| 一级,二级,三级黄色视频| 人人澡人人妻人| 欧美成人午夜免费资源| 免费观看av网站的网址| 久久人人爽人人爽人人片va| 欧美精品亚洲一区二区| 高清视频免费观看一区二区| 久久久久精品久久久久真实原创| 欧美人与性动交α欧美软件 | 两性夫妻黄色片 | 欧美丝袜亚洲另类| 成年人免费黄色播放视频| 欧美人与性动交α欧美精品济南到 | 成人午夜精彩视频在线观看| 在线观看免费日韩欧美大片| 看免费av毛片| 免费黄频网站在线观看国产| 22中文网久久字幕| av国产精品久久久久影院| 午夜福利视频在线观看免费| 男的添女的下面高潮视频| 人人妻人人添人人爽欧美一区卜| 97在线人人人人妻| 欧美精品一区二区大全| 视频中文字幕在线观看| 亚洲国产精品国产精品| 国产成人精品久久久久久| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 成人毛片a级毛片在线播放| 黑人欧美特级aaaaaa片| 亚洲av在线观看美女高潮| 永久网站在线| 在线观看国产h片| 精品人妻一区二区三区麻豆| av网站免费在线观看视频| 中文字幕av电影在线播放| 九九爱精品视频在线观看| 成人毛片a级毛片在线播放| 国产1区2区3区精品| 国产成人精品一,二区| 欧美人与性动交α欧美精品济南到 | 一区二区日韩欧美中文字幕 | 精品少妇内射三级| 亚洲国产精品国产精品| 母亲3免费完整高清在线观看 | 黄片无遮挡物在线观看| 性色avwww在线观看| 青春草视频在线免费观看| 免费观看a级毛片全部| 日日爽夜夜爽网站| 性色av一级| av.在线天堂| 亚洲国产欧美在线一区| 亚洲国产色片| 大香蕉久久网| 午夜久久久在线观看| 九九爱精品视频在线观看| 蜜臀久久99精品久久宅男| av不卡在线播放| 黄色一级大片看看| 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡 | 人妻人人澡人人爽人人| 国产精品不卡视频一区二区| 成人国语在线视频| 晚上一个人看的免费电影| 亚洲av欧美aⅴ国产| 91精品伊人久久大香线蕉| 久久久欧美国产精品| 又粗又硬又长又爽又黄的视频| 久热久热在线精品观看| 男女免费视频国产| 国产成人免费无遮挡视频| 热99久久久久精品小说推荐| 一级毛片 在线播放| 亚洲欧美日韩另类电影网站| 国产又爽黄色视频| 大话2 男鬼变身卡| 国精品久久久久久国模美| av视频免费观看在线观看| 国产精品国产三级国产av玫瑰| 在线亚洲精品国产二区图片欧美| 男女边摸边吃奶| 国产精品嫩草影院av在线观看| www.av在线官网国产| 狂野欧美激情性xxxx在线观看| 老女人水多毛片| 国产精品成人在线| a级毛片黄视频| 成人亚洲欧美一区二区av| 色视频在线一区二区三区| 热99久久久久精品小说推荐| 中文精品一卡2卡3卡4更新| 精品亚洲成国产av| 免费在线观看黄色视频的| 日韩一区二区三区影片| 在线看a的网站| 欧美日韩国产mv在线观看视频| 久久亚洲国产成人精品v| 国产男女内射视频| 97在线人人人人妻| 免费在线观看黄色视频的| 激情视频va一区二区三区| 欧美xxⅹ黑人| 99九九在线精品视频| 久久鲁丝午夜福利片| 久久影院123| 少妇被粗大的猛进出69影院 | 日韩欧美精品免费久久| 国产探花极品一区二区| 欧美国产精品一级二级三级| 又黄又爽又刺激的免费视频.| 国产精品秋霞免费鲁丝片| 丝袜在线中文字幕| 亚洲精品,欧美精品| 国产黄色视频一区二区在线观看| 99国产综合亚洲精品| 久久青草综合色| 色5月婷婷丁香| 国产片内射在线| 精品一区二区三卡| av免费在线看不卡| 久久精品国产亚洲av天美| 亚洲一区二区三区欧美精品| 久久久久网色| 国产av精品麻豆| 国产成人免费无遮挡视频| 久久精品国产亚洲av涩爱| 中国国产av一级| 国产黄频视频在线观看| av不卡在线播放| 肉色欧美久久久久久久蜜桃| 最后的刺客免费高清国语| 日韩熟女老妇一区二区性免费视频| 在线观看美女被高潮喷水网站| 国产伦理片在线播放av一区| 久久久精品94久久精品| 免费少妇av软件| 中国三级夫妇交换| av在线app专区| 亚洲国产精品专区欧美| 女人被躁到高潮嗷嗷叫费观| 男人操女人黄网站| 热re99久久精品国产66热6| 最新的欧美精品一区二区| 亚洲人与动物交配视频| 人成视频在线观看免费观看| 欧美人与性动交α欧美软件 | 久久青草综合色| 精品少妇黑人巨大在线播放| 久久热在线av| 免费看不卡的av| www.熟女人妻精品国产 | 亚洲av福利一区| 国产日韩欧美亚洲二区| 久久精品久久久久久噜噜老黄| 欧美 日韩 精品 国产| 91aial.com中文字幕在线观看| 免费黄色在线免费观看| 男女边摸边吃奶| 亚洲五月色婷婷综合| freevideosex欧美| a级毛色黄片| 91国产中文字幕| 妹子高潮喷水视频| 午夜av观看不卡| 亚洲国产欧美日韩在线播放| 哪个播放器可以免费观看大片| 婷婷成人精品国产| 哪个播放器可以免费观看大片| 亚洲精品456在线播放app| 在线观看免费日韩欧美大片| 激情视频va一区二区三区| 中文欧美无线码| 久久精品久久久久久久性| 欧美日韩精品成人综合77777| 成人二区视频| 中文字幕人妻熟女乱码| 欧美日韩亚洲高清精品| av福利片在线| 男人舔女人的私密视频| 国产熟女午夜一区二区三区| 晚上一个人看的免费电影| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 成人手机av| 久热久热在线精品观看| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 观看av在线不卡| 国产爽快片一区二区三区| 日韩制服骚丝袜av| 亚洲伊人色综图| 亚洲国产精品专区欧美| 日本91视频免费播放| 日韩中文字幕视频在线看片| 国产淫语在线视频| 男女国产视频网站| 高清视频免费观看一区二区| 综合色丁香网| 日韩,欧美,国产一区二区三区| 亚洲精品久久午夜乱码| 五月开心婷婷网| av电影中文网址| 国产福利在线免费观看视频| 一二三四在线观看免费中文在 | 丰满乱子伦码专区| 99re6热这里在线精品视频| 久久国产亚洲av麻豆专区| 欧美日韩综合久久久久久| 免费日韩欧美在线观看| 91午夜精品亚洲一区二区三区| 天堂俺去俺来也www色官网| 成人黄色视频免费在线看| 欧美成人午夜免费资源| 亚洲成国产人片在线观看| 精品一品国产午夜福利视频| 中文字幕最新亚洲高清| 久久久久精品人妻al黑| av卡一久久| 九九在线视频观看精品| www日本在线高清视频|