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

    An Improved Prediction Method for Structural Response Analysis of Offshore Wind Turbines

    2023-12-19 08:13:32-
    船舶力學(xué) 2023年12期

    -

    (a.State Key Laboratory of Ocean Engineering;b.School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

    Abstract:In order to overcome the shortcomings of the traditional Rosenblatt-IFORM (Inverse-First-Order Reliability Method) environmental contour method, the use of the PCA (Principal Component Analysis)-IFORM contour method is proposed for generating more accurate environmental contours.The environmental contours at a chosen offshore site obtained by using the PCA-IFORM contour method are compared with those obtained by using the Rosenblatt-IFORM environmental contour method,and the accuracy and effectiveness of the PCA-IFORM contour method are convincingly validated.The 50-year extreme structural dynamic responses of a monopile-supported 5 MW offshore wind turbine installed at this chosen offshore site are calculated based on the highest nonlinear sea state from the PCA-IFORM contour method (or the traditional Rosenblatt-IFORM environmental contour method).The calculation results are systematically analyzed and compared,and the necessity and importance of using more realistic environmental contours (such as those generated using the PCA-IFORM contour method)are finally highlighted.

    Key words:novel environmental contour method;principal component analysis;inverse-first-order reliability method;extreme structural dynamic response;offshore wind turbine

    0 Introduction

    In the ship classification rules for designing and building offshore wind turbines,the characteristic responses(extreme responses)are defined regarding a specific return period(e.g.a 50-year design service life).Accurate prediction of these characteristic response values requires correct identification of the long-term extreme sea states.In recent years, the environmental contour method has been advocated as a rational approach to identify an extreme sea state for predicting the long-term characteristic responses of an offshore wind turbine.The predicted environmental contours should be accurate enough to ensure the safety and economic viability of the offshore wind turbine whose dynamic responses are calculated based on the highest sea states along these contours during the design process.

    In the traditional Rosenblatt-IFORM environmental contour method[1-2],Rosenblatt transformation relates the environmental variables (e.g.the significant wave heightHSand the spectral peak periodTp) in the physical space with independent standard Gaussian variables (u1,u2) in the socalled standard Gaussian space.This Rosenblatt-IFORM environmental contour method was used by Saranyasoontorn and Manuel[3]to assess the accuracy of offshore wind turbine reliability-based design loads.The extreme response variables selected for the reliability-based design studies in Saranyasoontorn and Manuel[3]include the 50-year design shear forces and overturning moments at the sea bed of a bottom fixed offshore wind turbine.Karmakar et al[4]analyzed the reliability-based design loads based on the Rosenblatt-IFORM environmental contour method to estimate the longterm extreme loads for a 5 MW offshore wind turbine.In the environmental contour method of Karmakar et al[4],the identification of the environmental contours is carried out for a given return period(e.g.50-year)to predict the maximum response for the 5 MW offshore wind turbine.Horn and Winterstein[5]performed the extreme response estimation of a 10 MW offshore wind turbine mounted on a monopile using the Rosenblatt-IFORM environmental contour method.In their study, they created two-dimensional 50-year environmental contours regarding the significant wave heightHSand the spectral peak periodTp.Velarde et al[6]conducted the probabilistic analysis of a 10 MW monopile-supported offshore wind turbine under extreme response.In their study, they derived 50-yearHS-Tpcontours based on in-situ wave observations at the North Sea using the Rosenblatt-IFORM environmental contour method.

    However, the 50-year environmental contours created by using the Rosenblatt-IFORM method often fail to cover the measured data taken over a relatively short (8-20 year)time period.Therefore,the extreme structural dynamic responses predicted based on the 50-year Rosenblatt-IFORM environmental contours as reported in the aforementioned offshore wind turbine studies can hardly be deemed reliable.The shortcomings of the Rosenblatt-IFORM environmental contours obviously demonstrate that a new approach for generating more realistic and accurate environmental contours is imperative in order to predict more accurate extreme structural dynamic responses of offshore wind turbines.

    In this study, we will try to improve the traditional Rosenblatt-IFORM environmental contour method through principal component analysis (PCA) to create a set of uncorrelated variables that can be used to generate better-fitting environmental contours.In the following we will call this proposed approach the PCA -IFORM environmental contour method.The environmental contours at a chosen offshore site obtained by using the PCA-IFORM contour method will be compared with those obtained by using the Rosenblatt-IFORM environmental contour method, and the accuracy and effectiveness of the PCA-IFORM contour method will be convincingly validated.Next, we will calculate the 50-year extreme structural dynamic responses of a monopile-supported 5 MW offshore wind turbine installed at this chosen offshore site based on the highest nonlinear sea state from the PCA-IFORM contour method (or the traditional Rosenblatt-IFORM environmental contour method).Our calculation results will be systematically analyzed and compared, and the necessity and importance of using more realistic environmental contours (such as those generated using the PCA-IFORM contour method)will finally be highlighted.

    This paper begins in Chapter 1 by specifying the theories and formulas of the traditional approaches and the proposed novel approach for generating environmental contours.It continues in Chapter 2 to elucidate the theories for simulating nonlinear random waves after the highest sea state is identified along the generated environmental contour.In Chapter 3 the environmental contours at a chosen offshore site obtained by using the PCA-IFORM contour method will be compared with those obtained by using the Rosenblatt-IFORM environmental contour method.In Chapter 4, the extreme structural dynamic responses of a 5 MW offshore wind turbine based on the highest nonlinear sea state from the PCA-IFORM contour method (or the traditional Rosenblatt-IFORM environmental contour method) will be calculated and compared, with concluding remarks finally summarized in Chapter 5.

    1 Environmental contour methods

    1.1 Rosenblatt transformation

    For designing offshore structures, the environmental contour method is utilized for choosing the appropriate short-term design sea states leading to load and response extremes corresponding to a prescribed return period, which otherwise has to be obtained from a full long-term analysis.The environmental contour method can be applied for an offshore site if the joint distribution modelFHSTp(h,t)of the significant wave heightHSand the spectral peak periodTpis available.Environmental contour lines corresponding to a constant annual exceedance probability can be obtained by transforming the joint model to a space consisting of independent,standard normal variables (u1,u2)by first using the Rosenblatt transformation and then utilizing the Inverse First-Order Reliability Method(IFORM).The mathematical equations for the Rosenblatt transformation[7]are written as

    In Eqs.(1)-(2)Φ( )is the standard normal distribution function.In the standard normal space,the environmental contour line corresponding to an annual exceedance probability ofqwill be circles with radiusr=Φ-1( )1 -q/8760 where 8760 is the number of 1-hsea states per year[8].Theqprobability environmental contour lines can subsequently be obtained by transforming the aforementioned circles back to the physical parameter space using the following equations:

    In Eqs.(3)-(4),θis the equally-distanced angle within the range of [0,2π].In the following we call the afore-mentioned procedures for obtaining theq-probability environmental contour lines the Rosenblatt-IFORM environmental contour method.

    1.2 IFORM environmental contour method based on Principal Component Analysis(PCA)

    In some circumstances,the 50-year environmental contours obtained using the Rosenblatt-IFORM method fail to cover the observations taken over a relatively short (10-20 year) period of record.To overcome this shortcoming,the proposed principal component analysis can be utilized to create a set of uncorrelated variables that can be used to generate better-fitting environmental contours.

    Given a collection of points in two,three,or higher dimensional space,a‘best fitting’line can be defined as one that minimizes the average squared distance from a point to the line.The next best-fitting line can be similarly chosen from directions perpendicular to the first.Repeating this process yields an orthogonal basis in which different individual dimensions of the data are uncorrelated.These basis vectors are called principal components.Principal component analysis is defined as an orthogonal linear transformation that transforms the data to a new coordinate system such that the greatest variance by some scalar projection of the data comes to lie on the first coordinate(called the first principal component),the second greatest variance on the second coordinate,and so on.

    The principal component analysis is used to remove the correlation between the significant wave heightHSand the spectral peak periodTpfor a measured wave data set at an offshore site.The application of the principal component analysis generates two new variablesC1(the first principal component) andC2(the second principal component), andC1is the component that has the greatest variance.The application of the principal component analysis to the originalHSandTpdata set yields a coefficient matrix as shown in Eq.(5).This coefficient matrix defines a linear combination that allows for rotation into the principal component space[9].

    The components of the coefficient matrix in Eq.(5)have the following relations:v1,1= -v2,2andv2,1=v1,2.The values for each component can then be obtained by applying the coefficient matrix in Eq.(5)for each data point(HSi,Tp)as follows[9]:

    whereC1iis theith rotated first component corresponding to anith original significant wave height data pointHSiand anith original energy period data pointTei,C2iis theith rotated second component corresponding to anith original significant wave height data pointHSiand anith original energy period data pointTei,sis a shifting distance.By utilizing the following two equations, we can transform the point (C1i,C2i) in the principal component space back to the corresponding point in the original input space:

    After rotating the input variables into the principal component space,we can fit the cumulative distribution function ofC1with an inverse Gaussian distribution:

    whereΦ( )is the standard Gaussian distribution function, andμandλare the location and scale of the inverse Gaussian distribution respectively.In Eq.(12)xstands forC1.In probability theory,the inverse Gaussian distribution (also known as the Wald distribution)is a two-parameter family of continuous probability distributions.The inverse Gaussian distribution has several properties analogous to a Gaussian distribution.The name can be misleading: it is an‘inverse’only in that, while the Gaussian describes a Brownian motion's level at a fixed time, the inverse Gaussian describes the distribution of the time of a Brownian motion with positive drift takes to reach a fixed positive level.Its cumulant generating function (logarithm of the characteristic function)is the inverse of the cumulant generating function of a Gaussian random variable.

    The values ofC2are sorted according to their correspondingC1values and split into bins.The parameters of theC2distributions are fit as functions of the representative value ofC1for each bin in order to represent the dependency betweenC1andC2that has not yet been captured.Based on the shape of the cumulative distribution functions for each bin ofC2, a Gaussian distribution expressed as follows is chosen to fit each CDF ofC2:

    whereerfis the error function,and the meanμand the standard deviationσare the parameters defining the Gaussian distribution.In the following we will use a simple linear approximation to fitμas a function ofC1.Furthermore, we will use a constrained quadratic approximation to fitσas shown in Eqs.(15)-(17):

    As soon as the sun rose the King came, and when he perceived the gold he was astonished and delighted, but his heart only lusted more than ever after the precious metal

    In the following we call the procedures in this section for obtaining theq-probability environmental contour lines the PCA-IFORM environmental contour method.Robertson et al[10]performed the holistic marine energy resource assessments from a wave and offshore wind perspective of metocean conditions.In their research they generated accurate environmental contours using the WDRT toolbox[11]based on the PCA-IFORM environmental contour method.Li and Zhang[12]conducted long-term fatigue damage assessment for a floating offshore wind turbine under realistic environmental conditions.They pointed out that building an accurate environmental model that could accurately describe the multivariate dependence among environmental parameters is essential to ensure structural lifetime safety and reliability,and they particularly mentioned that the PCA-IFORM environmental contour statistical model could be used to achieve this objective.Yang et al[13]studied the characteristics and variability of the nearshore wave resource on the U.S.West Coast.They pointed out that environmental contours are widely used in the maritime industry for design,including the design of wave energy converters.Taking advantage of the long-term hindcast data set,environmental contours were computed by them at the four study locations using the PCA-IFORM environmental contour method as described by Eckert-Gallup et al[9].They found that the wave heights corresponding to the 100-year event(1%yearly probably of occurrence)decrease with latitude.Coe et al[14]performed a full long-term design response analysis of a wave energy converter.Environmental characterization and selection of sea states for this analysis at the intended deployment site were performed using principle-components analysis.Sea states for full long-term analyses in this study were obtained using a modified inverse first-order reliability (I-FORM) method (i.e., PCA-IFORM environmental contour method).This process was conducted using the tools provided in the WEC Design Response Toolbox(WDRT)[11].Finally,Coe et al[15]carried out a classical survey of WEC reliability, survival and design practices.In their paper they pointed out that‘An improved version of the I-FORM process is presented in Ref.[9].Here, a principle component analysis (PCA) approach is applied to better utilize the physical relationship between wave period and height.’.

    2 Theories for simulating nonlinear random waves

    In a classical paper[16]regarding estimating the extreme structural dynamic responses of offshore structures based upon the environmental contour method, the authors pointed out that‘The worst sea state turns out to be the highest sea state along the contour.’.After the highest sea state is identified along the PCA- IFORM (or Rosenblatt-IFORM) environmental contour, we can simulate nonlinear wave surface elevations based on the data of this highest sea state according to the following theories:

    wheredis the water depth.

    3 Environmental contours calculation examples based on the measured data at National Data Buoy Center(NDBC)Station 42012

    Fig.2 shows the 50-year environmental contours generated by the PCA-IFORM method and the Rosenblatt-IFORM method for NDBC.The red‘points’in Fig.2 show the measured data for the NDBC station 42012.The solid blue line in Fig.2 represents the 50-year environmental contour generated using the Rosenblatt-IFORM environmental contour method.We can

    Fig.1 Buoy of the NDBC Station 42012

    In this chapter we will present our environmental contours calculation examples based on the measured data at NDBC Station 42012.This buoy is installed at Orange Beach-44 nautical miles southeast of Mobile, AL, U.S.A.The water depth at this location is 25.9 m.Fig.1 shows the buoy of the NDBC Station 42012.In the following we will use the PCA- IFORM method and the Rosenblatt-IFORM method to generate environmental contours based on the 86089 hourly observations ofHSandTptaken from Jan.1, 2009 to Dec.31, 2019 at NDBC Station 42012.clearly find that this 50-year solid blue contour does not include all of the measured NDBC 42012 data of about 11 years.Furthermore, this solid blue contour does not follow the shape of trends of the measured NDBC 42012 data either.This solid blue contour fails for coverage in the area where theHSvalues are very high.These calculation results demonstrate that the obtained Rosenblatt-IFORM environmental contour does not have good predictivity.

    Fig.2 50-year environmental contours created by the PCA-IFORM method and the Rosenblatt-IFORM method for NDBC 42012

    The solid green curve in Fig.2 shows the 50-year environmental contour generated by using the PCA-IFORM method.We can clearly find that this solid green curve fit quite well with both higher and smaller values ofHS.We can also clearly notice that this solid green curve does include most of the measured NDBC 42012 data, missing only three data points.Furthermore, this solid green curve also follows quite well the shape of trends of the measured NDBC 42012 data.This solid green curve is very successful for coverage in the area where theHSvalues are very high.These calculation results indicate that the environmental contour generated using the PCA-IFORM method does have good predictivity.

    4 Extreme structural dynamic responses of a 5 MW offshore wind turbine

    4.1 Chosen 5 MW offshore wind turbine

    In the following we will calculate the extreme structural dynamic responses of the NREL (National Renewable Energy Laboratory) 5 MW monopilesupported offshore wind turbine (as shown in Fig.3)based on the theoretical methods as illustrated in Chapter 1 and Chapter 2.This chosen offshore wind turbine system is composed of a 5 MW wind turbine, a tower and a monopile foundation[17].The NREL 5MW baseline offshore wind turbine has three blades.The rotor diameter is 126 m and the hub diameter is 3 m.The rotor mass is 110 000 kg, and the nacelle mass is 240 000 kg.The NREL 5 MW fixed-bottom monopile wind turbine is sited at a water depth of 25.9 m.At the seabed mudline, the tower attaches to the monopile foundation of the same dimensions (with an outer diameter of 6 m and a wall thickness of 6 cm).The foundation is embedded 36 m into the seabed.The hub height of the wind turbine is 84.1 m from the mean sea water line.The draft of this offshore wind turbine is 25.9 m.These detailed properties and dimensions of this offshore wind turbine system is also contained in Re.[17].

    Fig.3 NREL 5 MW monopile-supported offshore wind turbine system

    4.2 Calculation examples

    From Fig.2 (50-year environmental contours created by the PCA-IFORM method and the Rosenblatt-IFORM method for NDBC 42012) we can find that in the highest sea state (50 years of extreme sea conditions) along the PCA-IFORM environmental contour,HS=6.47 m andTp=10.16 s.We can also find that in the highest sea state (50 years of extreme sea conditions)along the Rosenblatt-IFORM environmental contour,HS=4.03 m andTp=9.87 s.In the following calculation examples,the aforementioned NREL 5 MW offshore wind turbine is assumed to be installed at a site with an IEC (International Electrotechnical Commission) Class I-B wind regime.The turbulence wind random field over the turbine rotor plane is simulated using the computer code TurbSim with a wind speed of 1.6 m/s at the hub height.

    The structural dynamic responses of the aforementioned NREL 5 MW offshore wind turbine have been calculated by implementing time domain simulations in FAST.Hydrodynamic loads acting on the offshore wind turbine support structure are computed by applying Morison's equation.The nonlinear random ocean waves have been simulated using the theories in Chapter 2 based on a JONSWAP wave spectrum with a significant wave height of 6.47 m and a spectral peak period of 10.16 s.For this highest sea state along the PCA-IFORM environmental contour in Fig.2, we have imple-mented six 10 min stochastic time domain simulations.Tab.1 shows the first and second random seeds of incident waves during these six simulations.Fig.4 shows one simulated 10 min time series of the fore-aft shear forces at the mudline of the offshore wind turbine.In our study we have joined the aforementioned six 10 min time series into one and tried to forecast the 50-year extreme shear force based on this joint one series by using the peaks-over-threshold method[18-20].

    Tab.1 1st and 2nd random seeds of incident waves during 6 simulations(PCA-IFORM case)

    Fig.4 One simulated 10 min time series of the fore-aft shear forces at the mudline of the offshore wind turbine

    For the joint 60 min time series of the foreaft shear forces at the mudline of the offshore wind turbine we have obtained an optimal threshold level of 2 334 500 N as shown by the green vertical line in the mean excess plot (i.e.,Fig.5).We have extracted 18 local extremes from the aforementioned time series above the threshold level of 2 334 500 N and fitted a Generalized Pareto distribution as shown by the solid green line in Fig.6 using the maximum likelihood method.In Fig.6,the blue straight lines are the empirical distribution (histogram) obtained from these extracted 18 local extremes.We can find that the fitted Generalized Pareto distribution fits the empirical distribution (histogram) very well.The two red dashed curves in Fig.6 are the upper and lower 95%confidence interval lines respectively.

    In this study we want to predict the 50-year extreme structural dynamic responses of the NREL 5 MW monopile-supported wind turbine.That is to say, our aimed exceedance probability target is 2×10-2.As pointed out in the theory of the environmental contour method in Ref.[8]:‘The fractile level,p, will depend on the aimed exceedance probability target,q, …… For most practical systems,p=0.90 seems reasonable forq=10-2, whilep=0.95 may be more adequate forq=10-4.’.Therefore,it is reasonable for us to choosep=0.85 because our aimed exceedance probability target is 2×10-2.Consequently, based on the aforementioned fitted Generalized Pareto distribution we finally obtain that the 0.85 fractile level load (fore-aft shear force at the mudline) is 2.922 2E+06 N, and this number is summarized in Tab.3.Other data in Tab.3 have been obtained in a similar manner.

    Fig.5 Mean excess plot for the fore-aft shear forces at the mudline of the offshore wind turbine

    Fig.6 Fitted generalized Pareto distribution as shown by the solid green line

    Tab.2 1st and 2nd random seeds of incident waves during 6 simulations(Rosenblatt-IFORM case)

    Here we should further explain the calculation processes for obtaining the data in Column 3 of Tab.3, i.e., the calculation results based on the extreme sea state (HS=4.03 m andTp=9.87 s) obtained using the Rosenblatt-IFORM environmental contour method.For this highest sea state along the environmental contour based on the Rosenblatt-IFORM contour method, we have also carried out six 10 min stochastic time domain simulations.Tab.2 shows the first and second random seeds of incident waves during these six simulations.Subsequently, using the similar procedures for obtaining the fore-aft shear force result list in Column 2 of Tab.3 we have finally obtained the data as listed in Column 3 of Tab.3.

    Tab.3 50-year extreme response values based on 50-year extreme sea states obtained using different methods

    By carefully studying the calculation results listed in Tab.3 we can find that the 50-year extreme fore-aft shear force value based on the 50-year extreme sea state obtained using the PCA-IFORM contour method is 27.6% larger than the corresponding value obtained based on the Rosenblatt-IFORM contour method.The 50-year extreme fore-aft bending moment at the mudline based on the extreme sea state obtained using the PCA-IFORM contour method is 10.6% larger than the corresponding value obtained based on the Rosenblatt-IFORM contour method.Because we have already demonstrated in Chapter 3 that the PCA-IFORM contour method can generate more accurate and realistic environmental contours than the Rosenblatt-IFORM contour method,the 50-year extreme responses(fore-aft shear force and bending moment at the mudline)predicted based on the 50-year extreme sea state obtained using the PCA-IFORM contour method are undoubtfully more accurate.This means that an offshore wind turbine design based on the extreme dynamic response values predicted using the highest sea state data along a PCA-IFORM contour will be safer.On the contrary, structural failures will likely occur for an offshore wind turbine design when the input extreme sea state is forecasted using the Rosenblatt-IFORM environmental contour method.

    5 Concluding remarks

    In order to overcome the shortcomings of the traditional Rosenblatt-IFORM environmental contour method, the use of the PCA-IFORM contour method for generating more accurate environmental contours was proposed.Based on the measured data at a chosen offshore site, it was demonstrated that the PCA-IFORM contour method can generate more accurate and realistic environmental contours than the Rosenblatt-IFORM environmental contour method.The 50-year extreme structural dynamic responses of a monopile-supported 5 MW offshore wind turbine installed at this chosen offshore site were calculated based on the PCA-IFORM contour method and the traditional Rosenblatt-IFORM environmental contour method.It is concluded that an offshore wind turbine design based on the extreme dynamic response values predicted by using the highest sea state data along a PCA-IFORM contour will be safer while an offshore wind turbine design will likely result in structural failures when the input extreme sea state is forecasted using the Rosenblatt-IFORM environmental contour method.In summary, this study highlights the necessity and importance of using more realistic environmental contours (such as those generated using the PCA-IFORM contour method)in prediction of the extreme structural dynamic responses of offshore wind turbines.

    国产精品免费一区二区三区在线| 国产精品人妻久久久久久| av黄色大香蕉| 国产精品1区2区在线观看.| 久久国产乱子免费精品| 天堂√8在线中文| 狠狠狠狠99中文字幕| 久9热在线精品视频| 色尼玛亚洲综合影院| 日本熟妇午夜| 一a级毛片在线观看| 成人特级av手机在线观看| 久久久色成人| 免费观看精品视频网站| 夜夜看夜夜爽夜夜摸| 久久婷婷人人爽人人干人人爱| 免费大片18禁| 九色成人免费人妻av| 国产爱豆传媒在线观看| 欧美日韩综合久久久久久 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲欧美日韩卡通动漫| 麻豆国产97在线/欧美| 国产av不卡久久| 深夜精品福利| 欧美性猛交黑人性爽| 美女xxoo啪啪120秒动态图 | 免费人成在线观看视频色| netflix在线观看网站| 亚洲一区二区三区不卡视频| 欧美绝顶高潮抽搐喷水| 精品不卡国产一区二区三区| 久久精品人妻少妇| 深夜a级毛片| 色综合亚洲欧美另类图片| 99在线人妻在线中文字幕| 成人特级av手机在线观看| 日本a在线网址| 男人的好看免费观看在线视频| 18禁黄网站禁片午夜丰满| 性欧美人与动物交配| 中文字幕人妻熟人妻熟丝袜美| 日日摸夜夜添夜夜添小说| 国产成+人综合+亚洲专区| 亚洲av二区三区四区| 欧美+亚洲+日韩+国产| 最新中文字幕久久久久| 午夜a级毛片| 久久人人爽人人爽人人片va | 国产精品亚洲av一区麻豆| 欧美中文日本在线观看视频| 欧美另类亚洲清纯唯美| 人妻丰满熟妇av一区二区三区| a级一级毛片免费在线观看| 熟妇人妻久久中文字幕3abv| 99视频精品全部免费 在线| 大型黄色视频在线免费观看| 亚洲乱码一区二区免费版| 亚洲激情在线av| 丰满人妻一区二区三区视频av| 丰满人妻一区二区三区视频av| 男人舔女人下体高潮全视频| 国产爱豆传媒在线观看| 人妻制服诱惑在线中文字幕| 国内少妇人妻偷人精品xxx网站| 一a级毛片在线观看| 久久久久久大精品| 18禁裸乳无遮挡免费网站照片| 麻豆av噜噜一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 能在线免费观看的黄片| 久久九九热精品免费| 桃色一区二区三区在线观看| 老司机午夜十八禁免费视频| 欧美乱妇无乱码| 久久久国产成人精品二区| 中文在线观看免费www的网站| 中文字幕久久专区| av黄色大香蕉| 久久伊人香网站| 免费观看精品视频网站| 亚洲一区高清亚洲精品| 毛片女人毛片| 久久午夜福利片| 岛国在线免费视频观看| 午夜亚洲福利在线播放| 国产野战对白在线观看| 人人妻人人澡欧美一区二区| 国产91精品成人一区二区三区| 亚洲成人中文字幕在线播放| 精品人妻视频免费看| 久久欧美精品欧美久久欧美| 麻豆久久精品国产亚洲av| 一进一出抽搐gif免费好疼| 免费av不卡在线播放| 麻豆成人av在线观看| 精品福利观看| 18+在线观看网站| 窝窝影院91人妻| 免费人成在线观看视频色| 看片在线看免费视频| 乱码一卡2卡4卡精品| 欧美高清性xxxxhd video| 亚洲精品456在线播放app | 成人午夜高清在线视频| 亚洲av日韩精品久久久久久密| 精品免费久久久久久久清纯| 国产精品久久久久久人妻精品电影| 午夜激情欧美在线| 亚洲人成网站高清观看| 国产精品国产高清国产av| 亚洲电影在线观看av| 在线观看一区二区三区| 深爱激情五月婷婷| 最近最新免费中文字幕在线| 成人性生交大片免费视频hd| 国产精品久久久久久久久免 | 18+在线观看网站| 亚洲最大成人中文| 又粗又爽又猛毛片免费看| 亚洲自拍偷在线| 一a级毛片在线观看| 国产精品一及| 精品久久国产蜜桃| 麻豆国产av国片精品| 如何舔出高潮| 人妻丰满熟妇av一区二区三区| 国产av麻豆久久久久久久| 看十八女毛片水多多多| 欧美区成人在线视频| 不卡一级毛片| 国产极品精品免费视频能看的| 亚洲精品亚洲一区二区| 免费电影在线观看免费观看| 免费观看人在逋| bbb黄色大片| 99热这里只有是精品50| .国产精品久久| 成人国产一区最新在线观看| 亚洲自偷自拍三级| 一区二区三区免费毛片| 欧美另类亚洲清纯唯美| 国产探花极品一区二区| 国产男靠女视频免费网站| 久久久久久国产a免费观看| av女优亚洲男人天堂| 啦啦啦韩国在线观看视频| 久久香蕉精品热| 久久精品国产清高在天天线| 在线看三级毛片| 黄色一级大片看看| 久久久久久久午夜电影| 黄色丝袜av网址大全| 色哟哟哟哟哟哟| 久久久久久大精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人久久性| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| 波多野结衣高清无吗| 丝袜美腿在线中文| 天堂av国产一区二区熟女人妻| 亚洲成人精品中文字幕电影| 免费观看的影片在线观看| 国产精品人妻久久久久久| avwww免费| 欧美成狂野欧美在线观看| 午夜日韩欧美国产| 欧美午夜高清在线| 在现免费观看毛片| 欧美精品啪啪一区二区三区| 国产av不卡久久| 久久久国产成人精品二区| 变态另类成人亚洲欧美熟女| 噜噜噜噜噜久久久久久91| 人妻久久中文字幕网| 国产精品亚洲美女久久久| 国产精品综合久久久久久久免费| 国产高清三级在线| 伊人久久精品亚洲午夜| 搡老岳熟女国产| 中文字幕熟女人妻在线| 成人性生交大片免费视频hd| 国产野战对白在线观看| 亚洲av成人av| 99riav亚洲国产免费| 国产精品久久电影中文字幕| 久久久久久久久久成人| 亚洲乱码一区二区免费版| or卡值多少钱| 非洲黑人性xxxx精品又粗又长| 中文字幕高清在线视频| 精品人妻偷拍中文字幕| 床上黄色一级片| bbb黄色大片| netflix在线观看网站| 99国产极品粉嫩在线观看| 淫秽高清视频在线观看| 免费在线观看成人毛片| 18+在线观看网站| 99国产极品粉嫩在线观看| 中文字幕熟女人妻在线| 在线看三级毛片| 搡老熟女国产l中国老女人| 亚洲av成人不卡在线观看播放网| or卡值多少钱| 日本黄色片子视频| 一个人看视频在线观看www免费| 在线观看一区二区三区| 日韩av在线大香蕉| 色哟哟·www| 国产极品精品免费视频能看的| 亚洲一区二区三区色噜噜| 2021天堂中文幕一二区在线观| 噜噜噜噜噜久久久久久91| 国产日本99.免费观看| 在线播放无遮挡| 久久人人精品亚洲av| 午夜亚洲福利在线播放| 蜜桃亚洲精品一区二区三区| 久久久久久久久久成人| 搡老岳熟女国产| 欧美zozozo另类| 亚洲第一电影网av| 亚洲人成网站在线播放欧美日韩| 国产精品免费一区二区三区在线| 男人的好看免费观看在线视频| 看免费av毛片| 久久热精品热| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 悠悠久久av| 亚洲欧美日韩卡通动漫| 看黄色毛片网站| 国产蜜桃级精品一区二区三区| 午夜激情欧美在线| 亚洲成人久久性| 欧美黑人巨大hd| 国产亚洲精品综合一区在线观看| 床上黄色一级片| 色尼玛亚洲综合影院| av天堂在线播放| a在线观看视频网站| 在线播放无遮挡| 九九久久精品国产亚洲av麻豆| 制服丝袜大香蕉在线| 一本久久中文字幕| 最近中文字幕高清免费大全6 | 中文字幕人妻熟人妻熟丝袜美| 精品国内亚洲2022精品成人| 看免费av毛片| 欧美+亚洲+日韩+国产| 我要看日韩黄色一级片| 欧美色视频一区免费| 一本精品99久久精品77| 精品久久久久久久人妻蜜臀av| 在线天堂最新版资源| 少妇人妻一区二区三区视频| 亚洲无线在线观看| 一区二区三区免费毛片| 久久精品国产亚洲av天美| 成人高潮视频无遮挡免费网站| 欧美+亚洲+日韩+国产| 五月玫瑰六月丁香| 高清毛片免费观看视频网站| 亚洲国产精品久久男人天堂| 国产午夜精品论理片| 最好的美女福利视频网| 免费高清视频大片| 18美女黄网站色大片免费观看| 国产成人啪精品午夜网站| 网址你懂的国产日韩在线| 在线观看av片永久免费下载| 黄色女人牲交| 精品久久国产蜜桃| 一本久久中文字幕| 中文字幕人妻熟人妻熟丝袜美| 熟女电影av网| 国产精品日韩av在线免费观看| bbb黄色大片| 看片在线看免费视频| 亚洲精品在线观看二区| 男女之事视频高清在线观看| 别揉我奶头~嗯~啊~动态视频| 乱码一卡2卡4卡精品| 国产黄片美女视频| 亚洲18禁久久av| 俺也久久电影网| 久久草成人影院| 中文字幕av在线有码专区| 999久久久精品免费观看国产| 婷婷六月久久综合丁香| 18禁在线播放成人免费| 99riav亚洲国产免费| 日韩欧美精品免费久久 | 91午夜精品亚洲一区二区三区 | 丁香欧美五月| 亚州av有码| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美激情综合另类| 国产 一区 欧美 日韩| 国产毛片a区久久久久| 日韩中文字幕欧美一区二区| eeuss影院久久| 成年女人永久免费观看视频| 亚洲最大成人手机在线| 日本三级黄在线观看| 久久精品国产亚洲av涩爱 | 久久精品国产亚洲av天美| 在线国产一区二区在线| 免费在线观看影片大全网站| 美女高潮喷水抽搐中文字幕| 亚洲精品久久国产高清桃花| 免费黄网站久久成人精品 | 亚洲熟妇熟女久久| 久久久国产成人免费| 国产又黄又爽又无遮挡在线| 一个人看的www免费观看视频| 97碰自拍视频| 制服丝袜大香蕉在线| 欧美国产日韩亚洲一区| 身体一侧抽搐| or卡值多少钱| 精品日产1卡2卡| 中文亚洲av片在线观看爽| 亚洲自偷自拍三级| 桃色一区二区三区在线观看| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 尤物成人国产欧美一区二区三区| 欧美一区二区精品小视频在线| 91久久精品国产一区二区成人| 内射极品少妇av片p| a级毛片a级免费在线| 丁香六月欧美| 久久久久久国产a免费观看| 国产探花在线观看一区二区| 亚洲精品粉嫩美女一区| 麻豆一二三区av精品| 麻豆国产97在线/欧美| 亚洲av成人精品一区久久| 亚洲综合色惰| 久久国产精品人妻蜜桃| 亚洲,欧美精品.| 丰满人妻熟妇乱又伦精品不卡| 中文字幕人成人乱码亚洲影| 内地一区二区视频在线| 搡女人真爽免费视频火全软件 | а√天堂www在线а√下载| 一本综合久久免费| 俄罗斯特黄特色一大片| 99热只有精品国产| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 欧美潮喷喷水| 色综合站精品国产| 日韩欧美在线乱码| 我要搜黄色片| 少妇高潮的动态图| 亚洲国产精品成人综合色| 国产精品乱码一区二三区的特点| 日本免费a在线| 国产熟女xx| 亚洲最大成人av| 99久久99久久久精品蜜桃| 国产精品爽爽va在线观看网站| 看黄色毛片网站| 露出奶头的视频| 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 国产探花在线观看一区二区| 欧美一区二区国产精品久久精品| 韩国av一区二区三区四区| 99国产综合亚洲精品| 在线播放无遮挡| 国产精品精品国产色婷婷| 久久久久免费精品人妻一区二区| 美女xxoo啪啪120秒动态图 | 国产精品嫩草影院av在线观看 | 亚洲国产精品sss在线观看| 首页视频小说图片口味搜索| 搡老岳熟女国产| 精品日产1卡2卡| 久久精品久久久久久噜噜老黄 | av福利片在线观看| 久久性视频一级片| 每晚都被弄得嗷嗷叫到高潮| bbb黄色大片| 波多野结衣高清无吗| 国产日本99.免费观看| 日本与韩国留学比较| 男人狂女人下面高潮的视频| 日本 av在线| 别揉我奶头 嗯啊视频| 欧美午夜高清在线| 高潮久久久久久久久久久不卡| 国产精品一区二区三区四区久久| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 91在线精品国自产拍蜜月| 人妻制服诱惑在线中文字幕| 色精品久久人妻99蜜桃| 欧美日韩黄片免| 好男人在线观看高清免费视频| 国产精品综合久久久久久久免费| xxxwww97欧美| 五月伊人婷婷丁香| 黄色丝袜av网址大全| 亚州av有码| 看免费av毛片| av在线天堂中文字幕| 亚洲一区二区三区不卡视频| 啪啪无遮挡十八禁网站| 99视频精品全部免费 在线| 午夜精品在线福利| 波多野结衣高清无吗| 美女高潮的动态| 精品久久久久久久久亚洲 | 国产亚洲精品久久久久久毛片| 午夜福利视频1000在线观看| 老熟妇仑乱视频hdxx| 热99re8久久精品国产| 亚洲乱码一区二区免费版| 欧美成人免费av一区二区三区| 一区福利在线观看| 国产成人a区在线观看| 精品久久国产蜜桃| 国产精品99久久久久久久久| 国产精品久久视频播放| 麻豆久久精品国产亚洲av| 热99re8久久精品国产| 午夜福利18| 校园春色视频在线观看| 亚洲,欧美精品.| 国产av麻豆久久久久久久| 亚洲性夜色夜夜综合| 国内精品一区二区在线观看| 嫁个100分男人电影在线观看| 12—13女人毛片做爰片一| 国内精品久久久久久久电影| 精品人妻视频免费看| 内地一区二区视频在线| 国产高清视频在线观看网站| 亚洲国产日韩欧美精品在线观看| 亚洲国产精品成人综合色| 熟女电影av网| 51国产日韩欧美| 少妇人妻一区二区三区视频| 日韩欧美国产在线观看| 琪琪午夜伦伦电影理论片6080| 欧美乱色亚洲激情| 老鸭窝网址在线观看| 午夜福利在线观看免费完整高清在 | 校园春色视频在线观看| 国产成年人精品一区二区| 少妇的逼水好多| 少妇裸体淫交视频免费看高清| 欧美激情久久久久久爽电影| 久久久久九九精品影院| 成人美女网站在线观看视频| 亚洲国产精品999在线| 色尼玛亚洲综合影院| 香蕉av资源在线| 人妻丰满熟妇av一区二区三区| 亚洲狠狠婷婷综合久久图片| 又黄又爽又免费观看的视频| 久久午夜亚洲精品久久| av天堂在线播放| 精品久久久久久成人av| 亚洲自偷自拍三级| 真人做人爱边吃奶动态| 69av精品久久久久久| 亚洲第一区二区三区不卡| 精品无人区乱码1区二区| 亚洲精品成人久久久久久| 国产又黄又爽又无遮挡在线| 噜噜噜噜噜久久久久久91| 又黄又爽又免费观看的视频| 男人和女人高潮做爰伦理| 在线观看免费视频日本深夜| 最近最新免费中文字幕在线| 欧美不卡视频在线免费观看| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 五月玫瑰六月丁香| 国产av一区在线观看免费| 美女免费视频网站| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区 | av在线老鸭窝| 国产91精品成人一区二区三区| 久久久久精品国产欧美久久久| 色综合亚洲欧美另类图片| 直男gayav资源| 久久婷婷人人爽人人干人人爱| 国产成人啪精品午夜网站| 国内揄拍国产精品人妻在线| 午夜福利成人在线免费观看| 少妇人妻一区二区三区视频| 老司机午夜福利在线观看视频| 亚洲成人精品中文字幕电影| 一级av片app| 成人三级黄色视频| 欧美中文日本在线观看视频| av国产免费在线观看| 久久草成人影院| 国内揄拍国产精品人妻在线| 国产探花在线观看一区二区| 一区二区三区四区激情视频 | 最新在线观看一区二区三区| 国产精品久久久久久精品电影| 看片在线看免费视频| 69av精品久久久久久| 我的老师免费观看完整版| 男人狂女人下面高潮的视频| 亚洲自偷自拍三级| 久久精品国产亚洲av涩爱 | 亚洲自拍偷在线| 国产 一区 欧美 日韩| 乱人视频在线观看| 亚洲成人免费电影在线观看| 夜夜夜夜夜久久久久| 亚洲av.av天堂| 亚洲国产日韩欧美精品在线观看| 99riav亚洲国产免费| 精品人妻1区二区| 色尼玛亚洲综合影院| 亚洲无线在线观看| 在线观看66精品国产| 成人午夜高清在线视频| 国产白丝娇喘喷水9色精品| 老鸭窝网址在线观看| 在线观看免费视频日本深夜| 国产精品嫩草影院av在线观看 | 波野结衣二区三区在线| 美女cb高潮喷水在线观看| 我的老师免费观看完整版| 99热6这里只有精品| 国产三级在线视频| 极品教师在线视频| 色综合亚洲欧美另类图片| 男女那种视频在线观看| 99在线人妻在线中文字幕| 国产在线男女| 一级毛片久久久久久久久女| 亚洲av五月六月丁香网| 中国美女看黄片| 麻豆久久精品国产亚洲av| 久久久久性生活片| 亚洲国产精品sss在线观看| 老司机午夜福利在线观看视频| 亚洲乱码一区二区免费版| 最新在线观看一区二区三区| 精品人妻1区二区| 又爽又黄无遮挡网站| 国产一区二区在线观看日韩| 国内久久婷婷六月综合欲色啪| 他把我摸到了高潮在线观看| 国产精品一区二区免费欧美| 99在线人妻在线中文字幕| 国产精品人妻久久久久久| 一个人免费在线观看电影| 国产白丝娇喘喷水9色精品| 欧美成人性av电影在线观看| 亚洲av二区三区四区| 色精品久久人妻99蜜桃| 亚洲人成网站在线播放欧美日韩| 亚洲av免费高清在线观看| 成人毛片a级毛片在线播放| 日本撒尿小便嘘嘘汇集6| 久久国产乱子免费精品| 亚洲avbb在线观看| 91麻豆精品激情在线观看国产| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 欧美绝顶高潮抽搐喷水| 欧美成人a在线观看| 午夜福利在线观看免费完整高清在 | 久久天躁狠狠躁夜夜2o2o| 欧美日韩黄片免| 日本撒尿小便嘘嘘汇集6| 少妇丰满av| 老女人水多毛片| 丰满乱子伦码专区| 婷婷亚洲欧美| 一级黄片播放器| 久久中文看片网| 成人高潮视频无遮挡免费网站| 日韩有码中文字幕| 亚洲激情在线av| 国产精品99久久久久久久久| 99精品久久久久人妻精品| 国产主播在线观看一区二区| 在线观看舔阴道视频| 日韩亚洲欧美综合| 草草在线视频免费看| 欧美日韩国产亚洲二区| 国产真实乱freesex| 亚洲综合色惰| 伊人久久精品亚洲午夜| 国产麻豆成人av免费视频| 美女黄网站色视频| 欧美激情久久久久久爽电影| 免费在线观看亚洲国产| 亚洲熟妇中文字幕五十中出| 亚洲精品一卡2卡三卡4卡5卡| 欧美黄色淫秽网站| 色哟哟·www| 国产精品久久久久久久电影| 欧美日韩亚洲国产一区二区在线观看| 综合色av麻豆| 精品国产亚洲在线| 内地一区二区视频在线| 国产精品女同一区二区软件 |