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

    Energy saving by using asymmetric aftbodies for merchant ships-design methodology, numerical simulation and validation*

    2016-12-26 06:51:50JieDANGHaoCHEN陳灝

    Jie DANG, Hao CHEN (陳灝)

    1. Maritime Research Institute Netherlands (MARIN), Wageningen, The Netherlands, E-mail: J.Dang@marin.nl

    2. Guangzhou Shipyard International Company Ltd., Guangzhou 510382, China

    Energy saving by using asymmetric aftbodies for merchant ships-design methodology, numerical simulation and validation*

    Jie DANG1, Hao CHEN (陳灝)2

    1. Maritime Research Institute Netherlands (MARIN), Wageningen, The Netherlands, E-mail: J.Dang@marin.nl

    2. Guangzhou Shipyard International Company Ltd., Guangzhou 510382, China

    The methodology and procedures are discussed on designing merchant ships to achieve fully-integrated and optimized hull-propulsion systems by using asymmetric aftbodies. Computational fluid dynamics (CFD) has been used to evaluate the powering performance through massive calculations with automatic deformation algorisms for the hull forms and the propeller blades. Comparative model tests of the designs to the optimized symmetric hull forms have been carried out to verify the efficiency gain. More than 6% improvement on the propulsive efficiency of an oil tanker has been measured during the model tests. Dedicated sea-trials show good agreement with the predicted performance from the test results.

    energy saving, hull form, asymmetric aftbody, model tests, sea trials

    Introduction

    Ship design for its powering performance has experienced continuous improvement on both hull forms and propulsors. Hull forms are, in principle, designed to accommodate certain types of propulsors, so that a single screw ship defers from a multiple-screw ship and a ship with azimuthing thrusters or pods differs from a ship with water jet systems, as the examples. Propulsors are tailored for the wake fields generated by the hulls, as so-called the “wake-adapted” propulsor designs[1]. For merchant ships, although cares have been taken in the hull form designs to achieve better wake fields at the propellers for a better cavitation performance and to lead more viscous wake into the propeller discs in order to increase hull efficiency, a real ‘propeller-adapted’ hull form design has not yet been fully developed.

    For decades, MARIN has been dedicated at developing sophisticated CFD-aided optimization procedures and tools to design both the hull forms and propeller geometry in order to achieve the fully-integra-ted designs with maximized propulsive efficiency and comfort for ships. It started from the optimization method for the resistance and the wake quality, for the shaft power based on series propellers[2], for the shaft power with actual design propellers by using RANSBEM coupling technique and gradually developed into“fully-integrated” optimization procedures recently, where the hull form and the propeller geometry can be designed simultaneously[3].

    The above mentioned techniques have been proven in practice to improve the powering performance of the single screw ships significantly, by both model tests and sea trials.

    When a ship with single screw and symmetric hull form is already optimized, the remaining energy losses behind the ship are then mainly the asymmetric slip-stream behind the propeller which rotates in one single direction, and the non-uniform axial flow which has typically a higher speed behind the down-running blades and lower speed behind the up-running blades. The energy losses and the principle of the ship propulsion has been already fully established in Ref.[4] and further elaborated in Ref.[5] with practical considerations. Recent studies have shown that both the transverse kinetic energy loss due to the swirl of the flow and also the axial kinetic energy loss due to the nonuniform axial velocity in the far field, behind the propeller and the rudder, are the two major energy losses that can be recovered by using the Energy SavingDevices (ESDs)[6,7].

    As already made clear by studies, ESDs often surfer from scale effects and fouling in full scale. Losing the stator blades of the ESDs, in some cases, and damaging the hull structure with cracks often make the ship owners worried about applying ESDs to their ships. Some ESDs have been cut off after a period of the ships’ life in service due to practical reasons[8-10].

    The merit of applying an asymmetric aftbody for a single-screw ship lies in the fact that the swirl and the non-uniform flow of the ship’s wake can be regulated in the same way (or even a better way) as by applying ESDs to improve the propeller-hull interactions, but without adding “additional appendages”onto the hull. This idea is hence reviving in recent years.

    Despite of the maneuvering problem of “Carlotti”, studied and tested at MARIN in the early 1950’s that confirmed 15% shaft power savings by applying an asymmetric aftbody[11,12], some container ships with asymmetric aftbodies have been built in the 1980’s[13]and further studies went on into late 1980’s and early 1990’s[14-17]. It could be due to the difficulties in designing and judging the asymmetric aftbodies by using only model tests in the 1980’s, and it could also be due to the difficulties in the ship’s arrangement and construction in full scale, asymmetric aftbodies were not prevailing. In the last decades however, the rapid growth of practical applications of the CFD in hull form and propeller blade designs opened a new era for applying asymmetric aftbodies for ships where a fullyintegrated hull-propeller design could be achieved to minimize the energy losses[18-21].

    The present paper will elaborate the design methodology and the optimization procedures, with a focus on the validation of the asymmetric ship designs by using both model tests and also full scale trial data. For the ship used as an example and discussed in this paper, the CFD evaluations of the propulsive performance of the asymmetric aftbody, the resistance and propulsion tests, the cavitation observations and the hull pressure fluctuation measurements, the course stability and maneuvering tests, and the sea trials have all been carried out with a final design propeller at MARIN. More than 6% improvement on the total propulsive efficiency for this ship with an asymmetric aftbody, with respect to an optimized symmetric ship, has been found at the design draught by model tests. The sea trial data show also a good agreement with the prediction from the model test results.

    1. Calculations

    1.1 Fully-integrated design and massive CFD calculations

    The RANS code-PARNASSOS, a code developed and frequently applied by MARIN, was used to calculate the flow along the hull. It predicts the steady turbulent flow around the ship hulls and solves the discretised Reynolds-averaged Navier-Stokes (RANS) equations for steady incompressible flow. Turbulence models are implemented in the code. Structured multi-block grids are used with a finite-difference discretization of second and third-order schemes for the various terms.

    This RANS code is coupled with a code based on Boundary Element Method (BEM) for potential flow simulation of propellers with blade sheet cavitationcalled PROCAL, developed under the framework of the Cooperative Research Ships (CRS). The coupling is done iteratively between the propeller loads and the volume forces for the RANS, until a self-propulsion balance has been found at a given shaft power or a given speed. This technique is called RANS-BEM coupling.

    Fig.1(a) Planes defined in front of the propeller for the extrapolation of the effective wake to the propeller plane. Plane No. 0 is at propeller plane, Plane Nos. 1, 2 and 3 are located at a distance x/ R=0.20, 0.25 and 0.30 in front of the propeller, respectively[9]

    Fig.1(b) The Pareto front of a 78 000 DWT bulk carrier, as an example[8]

    Fig.2 Comparisons of the total velocity field behind an operating propeller in model and full scales, for ships with symmetric and asymmetric aftbodies[8]

    The merit of this technique is that the induced velocity of the propeller can be calculated from the BEM code and subtracted from the total velocity field of the RANS calculations so that the effective wake field can be obtained, as shown in Fig.1(a). An effective wake field is essentially important for the propeller design.

    With thousands of hull form and propeller geometry combinations, a huge matrix of variations can be prepared for the analysis. CPU and memory requirements of the PARNASSOS are quite modest compared to most other commercial codes. The BEM calculations consumes negligible amount of CPU time, even if with cavitation simulations. However, for the hundreds of thousands of variations of hull forms and propellers, parallel computations are needed. At MARIN, the massive computations for all variants can be done at the computer cluster with 4 000 computational cores or distributed over a large number of idle PCs at night, using the Condor tools. Typically, more than thousands of variants can be calculated over one night.

    To judge the performance of the variants, various objective functions can be defined. For the power, either the resistance (the effective power) or the propeller shaft power of the ship can be used. For the comfort, either the pressure fluctuations/induced excitations or the uniformity of the nominal/effective wake (change of angle of attacks to the section profile of the propeller blades) can be used, where a wake objective function (WOF) is often defined as,

    The best compromise of the efficiency and comfort can be found at a Pareto front. As an example shown in Fig.1(b), the relative changes on the shaft power and the quality of the effective wake field are used. Based on the Pareto front, an optimized hull form and propeller geometry can be determined. For a fullyintegrated full-propeller design, the Pareto front is often based on the shaft power against the vertical excitation force of the cavitating propeller.

    More details of the optimization using automatic deformation of the hull and propeller geometry, the massive calculations and the numerical settings can be found in Refs.[8], [9] and [22].

    1.2 Wake field and kinetic energy

    Although the shaft power can be calculated directly by the RANS-BEM coupling technique and assessed for the ship’s powering performance, the kinetic energy losses in the ship’s far wake field, where the static pressure is completely recovered, provide indirect but easy assessment on the propulsive efficiency gain and are better illustrations of the velocity fields. The complete theory can be found in Ref.[4], however the kinetic energy losses in the wake can be simplified into Eqs.(2) to (4) below,

    where the total kinetic energy losses are split into the axial part Kaxand the transverse part Ktr.

    It might be not a bad idea to make an assessment of the kinetic energy losses just behind the propeller, as shown by the “behind wake disc” in Fig.10 of Ref.[8], instead of the real far field where numericaldissipation may become too strong. Although the static pressure just behind the propeller is not completely recovered, it is assumed that the difference between the symmetric hull and the asymmetric hull may not differ too much. The kinetic energy difference at that location may provide important information on the energy saving by applying an asymmetric aftbody in a relative sense.

    Table 1 Comparison of kinetic energy behind an operating propeller at self-propulsion condition[8]

    On the “behind wake disc” with a propeller in operation at the self-propulsion condition, the velocity wake fields are compared to each other in Fig.2, between the symmetric and the asymmetric aftbodies for both the model scale and the full scale.

    Fig.3(a) The large scale ship model with asymmetric aftbody with the rudder and the final design propeller

    Fig.3(b) The measured nominal wake field at the propeller disc in model scale[9]

    It can be seen from the wake field that: compared to the symmetric aftbody, the axial velocity field becomes more uniform for the asymmetric aftbody, especially the reduction of the peak on the starboard side where the propeller blades rotate downwards; on the portside, the tangential velocity is reduced significantly by the asymmetric aftbody. When using Eqs.(2) to (4), the kinetic energy losses can be quantified in a relative sense in Table 1.

    Fig.4(a) Shaft power savings and shaft rotational rate changes due to asymmetric aftbody (at design draught)

    Fig.4(b) Screen shots of the high speed underwater video’s of the cavitation observation at design draught and speed of the ship[8]

    Fig.5 MV Front Cheetah, the asymmetric aftbody in building and the sea trial results compared to the model test results

    It is clearly seen that the total reductions of the kinetic energy losses is about 15.5% in model scale and about 29.0% in full scale by the asymmetric aftbody, in the domain studied. The reduction of the kinetic energy in full scale doubles almost that in model scale. This may indicate more energy saving in full scale than in model scale by asymmetric aftbodies, as already experienced in full scale trials by Ref.[12].

    2. Model tests

    Ship models for both powering and maneuvering tests of the performance have been manufactured with the asymmetric aftbodies, as shown in Figs.3(a) and 3(b), in order to verify both efficiency gain and course stability of the ship for which people worried. The measured nominal wake filed is shown also in this figure where the pre-swirl against the rotation of the propeller and the reduced axial flow on the port side are clearly seen.

    It is rather similar to the test results of other asymmetric ships developed in the 1980’s[11,12]that slightly lower resistance has been measured in various draughts for the present ship with an asymmetric aftbody. Although it cannot be generalized that an asymmetric aftbody can reduce the ships resistance, both the kinetic energy analyses and the resistance tests show that an asymmetric aftbody has limited effect on the ship’s resistance if designed properly.

    By assuming that the amount of viscous wake, which goes through the propeller disc, is the same for the symmetric and asymmetric aftbodies, the same wake scaling can be used in the extrapolation of the self-propulsion test results. The speed-power-RPM relations are obtained and plotted in Fig.4(a) on the lefthand side for the design draught compared to the ship with the symmetric aftbody. A 6.6% shaft power saving has been found by the asymmetric aftbody. In addition, stable sheet cavitations have been observed and low pressure fluctuations on the hull surface have been measured, see the photos in Fig.4(b) on the righthand side.

    3. Sea trials

    A series of oil tankers discussed above have beenbuilt at Guangzhou Shipyard International Company Ltd. (GSI). Dedicated sea trials have been carried out by MARIN’s Trials and Monitoring Department, including speed-power trials and also manoeuvring trials.

    As the first ship of the series, the speed power trials have been performed on board MV Front Cheetah in January (see Fig.5), 2016 in the South China Sea, south of Guangzhou. The analysis of the measurements and the corrections to arrive at the specified contract conditions are done according to the STAIMO speed trial analysis protocols. The speed trial results are analysed using STAIMO V2.2, following ISO (2015).

    Fig.6 Comparison of the maneuvering trial results to the model test results

    The speed trials have been conducted in the ballast, design and scantling loading conditions. The ship speed, course, heading and position were measured by a Hemisphere twin DGPS. The DGPS antenna was installed on top of the bridge. The delivered propulsion power is determined by the power measurement system (PMS) on the shaft. This system measures with a strain gauge the deflection of the shaft, which is a measure for the torque. An optical sensor records the RPM of the shaft. During the trials, the water depth was around 70 m to 100 m and weather condition was moderate with wave height around 2.2 m to 2.3 m.

    The trial results were corrected for wind, wave and current. Water depth and displacement were carefully monitored and corrected too. The trial results, corrected for the ideal trial condition for the design draught, are shown on the right-hand side in Fig.5. The deviation between the trial results and the prediction from the model tests is less than 2%.

    Manoeuvring tests were also tried at the same time period. The results are compared to the model test results in Fig.6 and correlated very well with the predictions which confirmed the course stability of the ship.

    4. Conclusion

    The design methodology and design tools for a fully-integrated ship’s hull form and propeller geometry have been developed at MARIN in the last decades, with emphasis on single screw ships in order to maximize the total propulsive efficiency and comfort gains for merchant ships. The asymmetric aftbodies have been applied and shown promising results for the integration of hull and propeller designs, without adding any appendages such as ESDs. For the ship used as an example in this paper, more than 6% improvement on the total propulsive efficiency has been achieved and measured by model tests which were further verified by full scale trials with good cavitation performance and low vibration levels, alongside the ship’s good course stability.

    Acknowledgements

    The authors are grateful for the valuable support from Guangzhou Shipbuilding International Company Ltd. (GSI) for the present study. Thanks are also extended to Harry Willemsen, Luis Rueda, Auke van der Ploeg and Geerhard Janse of MARIN for lines and CFD calculations.

    [1] Oosterveld M. W. C. Wake adapted ducted propellers [D]. Doctoral Thesis, Delft, The Netherlands: Delft Technical University, 1970.

    [2] Van Der Ploeg A., Raven H. C. CFD-based optimization for minimal power and wake field quality [C]. Proceedings of 11th International Symposium on Practical Design of Ships and Other Floating Structures (PRADS). Rio de Janeiro, Brazil, 2010.

    [3] Van Der Ploeg A., Foeth E.-J. Optimization of a chemical tanker and propeller with CFD [C]. Proceedings of ECCOMAS International Conference on Computational Methods in Marine Engineering V. Hamburg, Germany, 2013.

    [4] Wald Q. Performance of a propeller in a wake and the interaction of propeller and hull [J]. Journal of Ship Research, 1965, 9: 1-8, 36.

    [5] Dyne G. The principles of propulsion optimization [J]. Transactions RINA, 1995, 137: 189-208.

    [6] Dang J., Chen H., Dong G. et al. An exploratory study on the working principles of energy saving devices (ESDs) [C]. Greenship’2011. Wuxi, China, 2011.

    [7] Dang J., Dong G., Chen H. An exploratory study on the working principles of energy saving devices (ESDs)-PIV, CFD investigations and ESD design guidelines [C]. The 31st International Conference on Ocean, Offshore and Arctic Engineering (OMAE 2012). 2012, paper No. OMAE2012-83053.

    [8] Dang J., Chen H., Rueda L. et al. Integrated design of asymmetric aftbody and propeller for an aframax tanker to maximize energy efficiency [C]. SMP’15. Austin, Texas, USA, 2015.

    [9] Rueda L., Dang J. Integrated design of asymmetric aftbody and propeller to maximize energy efficiency [C]. Proceedings of Energy Efficient Ships of the Royal Institution of Naval Architects (RINA Conference). Rotterdam, The Netherlands, 2015.

    [10] Abramowski T., ?elazny K., Szelangiewicz T. Numerical analysis of effect of asymmetric stern of ship on its screw propeller efficiency [J]. Polish Marine Research, 2010, 15(4): 13-16.

    [11] Collatz G. Fuel saving by asymmetric afterbodies for single screw vessels [J]. Ocean Engineering, 1983, 10(4): 213-226.

    [12] Collatz G. Non-symmetrical aft-bodies [C]. Proceedings of Workshop on Developments in Hull Form Design. Wageningen, The Netherland, 1985.

    [13] Stierman E. J., Osborne M. G. A review of the performance and application of the asymmetric stern [R]. Shell International Marine, SIM Report No. MR 13/86, 1986.

    [14] Stierman E. J. The design of an energy saving, wake adapted duct [C]. Proceedings of International Symposium on Practical Design of Ship and Other Floating Structures (PRADS’87). Trondheim, Norway, 1987.

    [15] Stierman E. J. Energy saving with an asymmetric stern and wake adapted ducts [C]. Proceedings of International Symposium on Ship Resistance and Powering Performance. Shanghai, China, 1986, 124-132.

    [16] Piskorz-Nalecki J. W. The influence of the asymmetric afterbody and the horizontal fins on the rudder blade on the propulsive efficiency [C]. Proceedings of Workshop on Developments in Hull Form Design. Wageningen, The Netherlands, 1985.

    [17] Blaurock J. An appraisal of unconventional aftbody configurations and propulsion devices [J]. Marine Technology, 1990, 27(6): 325-336.

    [18] Vaz G., Bosschers J. Modelling three dimensional sheet cavitation on marine propellers using a boundary element method [C]. Sixth International Symposium on Cavitation CAV2006. Wageningen, The Netherlands, 2006.

    [19] Bosschers J., Vaz G., Starke B. et al. Computational analysis of propeller sheet cavitation and propeller-ship interaction [C]. RINA conference “MARINE CFD2008”. Southampton, UK, 2008.

    [20] Starke B., Bosschers J. Analysis of scale effects in ship powering performance using a hybrid RANS-BEM approach [C]. 29th Symposium on Naval Hydrodynamics. Gothenburg, Sweden, 2012.

    [21] Rijpkema D., Starke B., Bosschers J. Numerical simulation of propeller-hull interaction and determination of the effective wake field using a hybrid RANS-BEM approach [C]. SMP’13. Launceston, Tasmania, Australia, 2015.

    [22] Veldhuis C., Gornicz T. and Scholcz T. Ship optimization using viscous flow computations in combination with generic shape variations and design of experiments [C]. Proceedings of PRADS 2016. Copenhagen, Denmark, 2016.

    (Received June 20, 2016, Revised August 13, 2016)

    * Biography:Jie DANG, Male, Ph. D., Senior Project Manager

    Hao CHEN,

    E-mail:chenhao@chinagsi.com

    麻豆成人av在线观看| 最近最新免费中文字幕在线| 午夜精品久久久久久毛片777| 国产精品久久电影中文字幕| 777久久人妻少妇嫩草av网站| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 国产av又大| 亚洲国产精品久久男人天堂| 婷婷精品国产亚洲av| 国产野战对白在线观看| 亚洲精品国产精品久久久不卡| 色精品久久人妻99蜜桃| 亚洲avbb在线观看| 老司机在亚洲福利影院| 久久国产精品男人的天堂亚洲| or卡值多少钱| 99久久国产精品久久久| 色综合婷婷激情| 久久久久久人人人人人| 国产一区二区在线av高清观看| 久久久久九九精品影院| 日本一本二区三区精品| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产欧美一区二区综合| 国产精品 欧美亚洲| 欧美国产精品va在线观看不卡| 波多野结衣高清无吗| 深夜精品福利| 国产真实乱freesex| 久久久久久久久免费视频了| 老汉色∧v一级毛片| 国产一卡二卡三卡精品| 午夜福利视频1000在线观看| 免费观看精品视频网站| 午夜福利高清视频| 久久精品影院6| 国产成人影院久久av| 亚洲国产中文字幕在线视频| 亚洲va日本ⅴa欧美va伊人久久| 国产野战对白在线观看| 日韩av在线大香蕉| 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 国产欧美日韩精品亚洲av| 给我免费播放毛片高清在线观看| 高潮久久久久久久久久久不卡| 熟女电影av网| 亚洲精品美女久久久久99蜜臀| 欧美在线一区亚洲| 欧洲精品卡2卡3卡4卡5卡区| 女同久久另类99精品国产91| 日本精品一区二区三区蜜桃| 亚洲人成电影免费在线| 91在线观看av| 亚洲精品一区av在线观看| 亚洲无线在线观看| 亚洲一区高清亚洲精品| 色老头精品视频在线观看| 美女午夜性视频免费| 免费女性裸体啪啪无遮挡网站| 桃色一区二区三区在线观看| 色婷婷久久久亚洲欧美| 久久99热这里只有精品18| 99久久无色码亚洲精品果冻| 观看免费一级毛片| 男人舔女人的私密视频| 日韩欧美免费精品| 亚洲国产精品久久男人天堂| 色哟哟哟哟哟哟| 日本精品一区二区三区蜜桃| 99热6这里只有精品| 亚洲av中文字字幕乱码综合 | 首页视频小说图片口味搜索| 曰老女人黄片| 国产91精品成人一区二区三区| 免费看十八禁软件| 久久青草综合色| 亚洲国产欧美日韩在线播放| 亚洲欧洲精品一区二区精品久久久| 国产精品一区二区精品视频观看| 日韩欧美一区二区三区在线观看| 国产精品久久视频播放| 日韩有码中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 夜夜躁狠狠躁天天躁| 成人欧美大片| 午夜激情av网站| ponron亚洲| 久久婷婷人人爽人人干人人爱| 视频区欧美日本亚洲| 中亚洲国语对白在线视频| 美女 人体艺术 gogo| 高清在线国产一区| 9191精品国产免费久久| 一区二区三区激情视频| 国产三级在线视频| 国产精品一区二区三区四区久久 | 亚洲av中文字字幕乱码综合 | 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 国产高清激情床上av| 免费看十八禁软件| 久久久久久久久中文| 神马国产精品三级电影在线观看 | 亚洲成a人片在线一区二区| 久久久久免费精品人妻一区二区 | 国产精品美女特级片免费视频播放器 | 两性夫妻黄色片| 夜夜躁狠狠躁天天躁| 日韩免费av在线播放| 90打野战视频偷拍视频| 精品欧美一区二区三区在线| 侵犯人妻中文字幕一二三四区| 国产91精品成人一区二区三区| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 亚洲欧美日韩无卡精品| 欧美成人一区二区免费高清观看 | 精品日产1卡2卡| 国产乱人伦免费视频| 国产欧美日韩一区二区三| 两性夫妻黄色片| 真人做人爱边吃奶动态| 国产真人三级小视频在线观看| 中亚洲国语对白在线视频| 人人澡人人妻人| 午夜福利高清视频| 免费搜索国产男女视频| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 久久中文字幕人妻熟女| 一本精品99久久精品77| 午夜视频精品福利| 国产成人精品无人区| 国产av一区在线观看免费| 久久精品91蜜桃| 亚洲精品国产精品久久久不卡| 桃色一区二区三区在线观看| 亚洲中文字幕日韩| 久久人妻福利社区极品人妻图片| 制服人妻中文乱码| 国产av在哪里看| 俄罗斯特黄特色一大片| 久久九九热精品免费| 欧美黄色淫秽网站| 国产99久久九九免费精品| 亚洲成人久久性| 亚洲精品美女久久久久99蜜臀| 俺也久久电影网| 特大巨黑吊av在线直播 | 亚洲片人在线观看| 久久久久九九精品影院| 国产成+人综合+亚洲专区| 午夜a级毛片| 国产色视频综合| 在线国产一区二区在线| 1024香蕉在线观看| 日韩一卡2卡3卡4卡2021年| 日韩国内少妇激情av| 久久精品国产亚洲av香蕉五月| 后天国语完整版免费观看| 久久 成人 亚洲| 色在线成人网| 99在线视频只有这里精品首页| 人人妻人人看人人澡| 国产一区在线观看成人免费| 久久久国产成人免费| 亚洲午夜理论影院| 欧美午夜高清在线| 国产97色在线日韩免费| 麻豆久久精品国产亚洲av| 天天添夜夜摸| 久久久久精品国产欧美久久久| 国产高清视频在线播放一区| 免费人成视频x8x8入口观看| 麻豆一二三区av精品| 国产在线观看jvid| 精品熟女少妇八av免费久了| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利高清视频| 久久久久久久精品吃奶| 色播在线永久视频| 国产免费av片在线观看野外av| 国产真实乱freesex| 不卡av一区二区三区| 在线免费观看的www视频| 美女 人体艺术 gogo| 免费无遮挡裸体视频| 天堂动漫精品| www国产在线视频色| 黄网站色视频无遮挡免费观看| 午夜福利免费观看在线| 一本精品99久久精品77| 亚洲国产毛片av蜜桃av| 午夜精品久久久久久毛片777| 不卡av一区二区三区| 亚洲精品色激情综合| 日韩一卡2卡3卡4卡2021年| 国产精品1区2区在线观看.| 国产精品,欧美在线| 久久久精品国产亚洲av高清涩受| 国产精品98久久久久久宅男小说| 成在线人永久免费视频| 怎么达到女性高潮| 精华霜和精华液先用哪个| 一区二区三区精品91| 9191精品国产免费久久| 免费高清视频大片| 国产精品 欧美亚洲| 91字幕亚洲| 亚洲国产精品成人综合色| 欧美最黄视频在线播放免费| 久久99热这里只有精品18| 69av精品久久久久久| 日本一本二区三区精品| 神马国产精品三级电影在线观看 | 日本三级黄在线观看| 亚洲一区二区三区色噜噜| 午夜免费观看网址| 免费高清在线观看日韩| 欧美乱妇无乱码| 欧美久久黑人一区二区| 啦啦啦免费观看视频1| 麻豆成人午夜福利视频| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 成年免费大片在线观看| 亚洲五月天丁香| 国产极品粉嫩免费观看在线| 久久午夜综合久久蜜桃| 一本大道久久a久久精品| 99久久久亚洲精品蜜臀av| 免费电影在线观看免费观看| 精品欧美一区二区三区在线| 人人妻,人人澡人人爽秒播| 男女视频在线观看网站免费 | 亚洲专区字幕在线| 国产欧美日韩精品亚洲av| 国产亚洲欧美精品永久| 村上凉子中文字幕在线| 在线观看日韩欧美| 白带黄色成豆腐渣| 国内毛片毛片毛片毛片毛片| 日日夜夜操网爽| 国产精品香港三级国产av潘金莲| 国产区一区二久久| 999精品在线视频| 国产高清有码在线观看视频 | 超碰成人久久| 国产精品自产拍在线观看55亚洲| 免费观看精品视频网站| 亚洲第一青青草原| 法律面前人人平等表现在哪些方面| 日本免费一区二区三区高清不卡| 香蕉国产在线看| 日韩精品中文字幕看吧| 午夜亚洲福利在线播放| 中文资源天堂在线| 国产午夜精品久久久久久| 亚洲国产中文字幕在线视频| 精华霜和精华液先用哪个| 日本一区二区免费在线视频| 丝袜在线中文字幕| 禁无遮挡网站| 午夜两性在线视频| 久久精品成人免费网站| 激情在线观看视频在线高清| 国产欧美日韩精品亚洲av| 在线看三级毛片| 热re99久久国产66热| 国产片内射在线| 99久久久亚洲精品蜜臀av| 日日干狠狠操夜夜爽| 久久久精品国产亚洲av高清涩受| 又紧又爽又黄一区二区| 亚洲av日韩精品久久久久久密| 在线观看日韩欧美| 中文字幕精品免费在线观看视频| 91麻豆精品激情在线观看国产| cao死你这个sao货| 久久午夜综合久久蜜桃| 日韩精品青青久久久久久| 黄频高清免费视频| 午夜久久久久精精品| 日本熟妇午夜| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 国产成人影院久久av| 热re99久久国产66热| 一区二区日韩欧美中文字幕| 久久久国产成人免费| 欧美国产精品va在线观看不卡| 国产精品,欧美在线| 18禁黄网站禁片午夜丰满| 757午夜福利合集在线观看| 亚洲激情在线av| 久久香蕉激情| 久久精品亚洲精品国产色婷小说| 丰满的人妻完整版| 欧美乱妇无乱码| 好看av亚洲va欧美ⅴa在| 婷婷六月久久综合丁香| 国产av又大| 一本大道久久a久久精品| 91麻豆精品激情在线观看国产| 18禁黄网站禁片午夜丰满| 国产精品国产高清国产av| 亚洲精品色激情综合| 精品久久蜜臀av无| 99久久久亚洲精品蜜臀av| 国产成人精品久久二区二区91| 精品一区二区三区av网在线观看| 999久久久国产精品视频| 黄频高清免费视频| 国产aⅴ精品一区二区三区波| 国产99白浆流出| 精品久久久久久久毛片微露脸| 亚洲成国产人片在线观看| 亚洲欧美精品综合一区二区三区| 精品国产美女av久久久久小说| 久9热在线精品视频| 成人欧美大片| 久久人人精品亚洲av| 日本免费a在线| 国产精品永久免费网站| 999久久久精品免费观看国产| 久久久久久国产a免费观看| 久久草成人影院| 亚洲欧美精品综合一区二区三区| 深夜精品福利| e午夜精品久久久久久久| 久久香蕉精品热| 国产av不卡久久| 亚洲一区高清亚洲精品| 国产乱人伦免费视频| 长腿黑丝高跟| 成人国语在线视频| 欧美不卡视频在线免费观看 | 日韩免费av在线播放| 午夜亚洲福利在线播放| 亚洲天堂国产精品一区在线| 999精品在线视频| 免费观看人在逋| 国产高清videossex| 欧美三级亚洲精品| 90打野战视频偷拍视频| 欧美日韩亚洲国产一区二区在线观看| 一级a爱视频在线免费观看| 亚洲精品色激情综合| 久久久久久亚洲精品国产蜜桃av| 婷婷亚洲欧美| 亚洲男人天堂网一区| 桃色一区二区三区在线观看| 久热爱精品视频在线9| 国产91精品成人一区二区三区| 两性夫妻黄色片| 亚洲av成人一区二区三| 美女午夜性视频免费| 欧美中文综合在线视频| netflix在线观看网站| 精品国产一区二区三区四区第35| 日韩av在线大香蕉| 国产人伦9x9x在线观看| 国产午夜福利久久久久久| 亚洲九九香蕉| 国产人伦9x9x在线观看| 男人舔女人的私密视频| 亚洲七黄色美女视频| 一本精品99久久精品77| 欧美+亚洲+日韩+国产| 欧美黑人巨大hd| 波多野结衣巨乳人妻| 一区福利在线观看| 精品欧美国产一区二区三| 中文字幕精品免费在线观看视频| 免费av毛片视频| 热99re8久久精品国产| 女性被躁到高潮视频| 长腿黑丝高跟| 久久久久国产一级毛片高清牌| 亚洲中文字幕一区二区三区有码在线看 | 99久久国产精品久久久| 老熟妇乱子伦视频在线观看| 亚洲国产欧洲综合997久久, | 亚洲精品一区av在线观看| 欧美成人一区二区免费高清观看 | 国产精品 国内视频| 一级毛片精品| av在线播放免费不卡| 精品熟女少妇八av免费久了| 国产精品亚洲美女久久久| 亚洲av美国av| 性欧美人与动物交配| 欧美日韩黄片免| 制服丝袜大香蕉在线| 两个人看的免费小视频| 丁香欧美五月| 欧美成狂野欧美在线观看| 亚洲午夜理论影院| 1024香蕉在线观看| 国产97色在线日韩免费| 亚洲国产高清在线一区二区三 | 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 日韩视频一区二区在线观看| www.熟女人妻精品国产| 此物有八面人人有两片| 18美女黄网站色大片免费观看| 一级毛片高清免费大全| 色综合亚洲欧美另类图片| 日韩一卡2卡3卡4卡2021年| 日韩中文字幕欧美一区二区| 天天添夜夜摸| 国产成人av激情在线播放| 精品久久久久久,| 国产av又大| 久热爱精品视频在线9| 国产精品精品国产色婷婷| 最新美女视频免费是黄的| 国产三级在线视频| 国产亚洲av嫩草精品影院| 日本熟妇午夜| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| 18禁观看日本| 中亚洲国语对白在线视频| 精品卡一卡二卡四卡免费| 精品乱码久久久久久99久播| 国产视频内射| 色哟哟哟哟哟哟| a级毛片在线看网站| 搡老熟女国产l中国老女人| 中文字幕人妻丝袜一区二区| 精品欧美一区二区三区在线| 91字幕亚洲| 久久九九热精品免费| 国产精品乱码一区二三区的特点| 久久精品影院6| 亚洲国产精品久久男人天堂| 97碰自拍视频| 亚洲五月婷婷丁香| 欧美性长视频在线观看| 啦啦啦 在线观看视频| 日本在线视频免费播放| 国产男靠女视频免费网站| 成人手机av| 日日夜夜操网爽| 1024香蕉在线观看| 老司机福利观看| 亚洲自拍偷在线| 国产麻豆成人av免费视频| 亚洲熟女毛片儿| 人成视频在线观看免费观看| 亚洲国产日韩欧美精品在线观看 | 亚洲av第一区精品v没综合| 国产1区2区3区精品| 亚洲男人天堂网一区| 色尼玛亚洲综合影院| 搡老岳熟女国产| 国产亚洲欧美98| 欧美日韩黄片免| 一本久久中文字幕| 精品福利观看| 51午夜福利影视在线观看| 中文资源天堂在线| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久久久久久久 | 日韩欧美免费精品| 在线国产一区二区在线| 亚洲熟妇中文字幕五十中出| 国产精品久久电影中文字幕| 亚洲国产精品999在线| 在线观看免费午夜福利视频| 亚洲 欧美 日韩 在线 免费| 色播亚洲综合网| 老司机深夜福利视频在线观看| 色综合站精品国产| 我的亚洲天堂| 在线免费观看的www视频| www.熟女人妻精品国产| 久久国产乱子伦精品免费另类| 成年人黄色毛片网站| 国产1区2区3区精品| 久久精品国产亚洲av香蕉五月| 97超级碰碰碰精品色视频在线观看| 免费无遮挡裸体视频| 国产亚洲欧美在线一区二区| 欧美一级a爱片免费观看看 | 国产av不卡久久| 黄片播放在线免费| 深夜精品福利| 久久人妻av系列| 一区福利在线观看| 国产单亲对白刺激| 亚洲真实伦在线观看| 午夜激情福利司机影院| 亚洲电影在线观看av| 18美女黄网站色大片免费观看| 精品国产一区二区三区四区第35| x7x7x7水蜜桃| 韩国精品一区二区三区| 最近在线观看免费完整版| 午夜两性在线视频| 麻豆av在线久日| 午夜激情av网站| 中亚洲国语对白在线视频| 又紧又爽又黄一区二区| 日本 欧美在线| 亚洲国产精品合色在线| 搡老妇女老女人老熟妇| 无限看片的www在线观看| av在线天堂中文字幕| 国产亚洲精品av在线| 精品久久久久久久毛片微露脸| 日日摸夜夜添夜夜添小说| АⅤ资源中文在线天堂| 啦啦啦观看免费观看视频高清| 欧美激情 高清一区二区三区| 最新美女视频免费是黄的| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 久久婷婷人人爽人人干人人爱| 久久伊人香网站| 婷婷六月久久综合丁香| 日本三级黄在线观看| 他把我摸到了高潮在线观看| www.精华液| 国产精品久久久久久精品电影 | 午夜a级毛片| 法律面前人人平等表现在哪些方面| 日韩欧美一区视频在线观看| 日本在线视频免费播放| 午夜福利一区二区在线看| 丰满人妻熟妇乱又伦精品不卡| 国产主播在线观看一区二区| 美女扒开内裤让男人捅视频| 日本 av在线| 午夜成年电影在线免费观看| 久久精品夜夜夜夜夜久久蜜豆 | 精品欧美国产一区二区三| 一区二区三区激情视频| 男人的好看免费观看在线视频 | 97碰自拍视频| 日本撒尿小便嘘嘘汇集6| 亚洲专区中文字幕在线| 三级毛片av免费| 婷婷精品国产亚洲av在线| av片东京热男人的天堂| 无遮挡黄片免费观看| 91麻豆av在线| 欧美黑人精品巨大| 日韩欧美国产在线观看| 久久久久久久久久黄片| 人妻丰满熟妇av一区二区三区| АⅤ资源中文在线天堂| 中文字幕高清在线视频| 俄罗斯特黄特色一大片| 超碰成人久久| 99国产精品一区二区蜜桃av| 成人国产综合亚洲| 最近最新免费中文字幕在线| av免费在线观看网站| 十八禁人妻一区二区| 亚洲三区欧美一区| 欧美色视频一区免费| 亚洲av成人不卡在线观看播放网| 国产区一区二久久| 美女免费视频网站| svipshipincom国产片| 国产精品久久久久久精品电影 | 在线免费观看的www视频| 亚洲人成网站高清观看| 十八禁网站免费在线| 看片在线看免费视频| 国产一卡二卡三卡精品| 50天的宝宝边吃奶边哭怎么回事| 国产三级黄色录像| 搡老熟女国产l中国老女人| 狠狠狠狠99中文字幕| 侵犯人妻中文字幕一二三四区| 成人手机av| 亚洲专区中文字幕在线| 中文字幕精品亚洲无线码一区 | 日韩欧美免费精品| 色老头精品视频在线观看| 一边摸一边抽搐一进一小说| 久久香蕉精品热| av福利片在线| 午夜福利18| 欧美中文综合在线视频| 国产精品,欧美在线| 黄色毛片三级朝国网站| 免费高清在线观看日韩| 欧美激情久久久久久爽电影| 18禁裸乳无遮挡免费网站照片 | 变态另类丝袜制服| 中文资源天堂在线| 免费搜索国产男女视频| 亚洲三区欧美一区| 三级毛片av免费| 国产又爽黄色视频| 日韩欧美国产在线观看| 99在线人妻在线中文字幕| 人人妻人人看人人澡| 少妇裸体淫交视频免费看高清 | 国产成人av激情在线播放| 亚洲无线在线观看| 欧美 亚洲 国产 日韩一| 亚洲精华国产精华精| 精品国产亚洲在线| 亚洲欧美激情综合另类|