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

    Study of Vortex Induced Characteristics of Multi-columns with Low Mass Ratio

    2014-01-19 05:49:42GUJiayangYANGJianminXIAOLongfei
    船舶力學(xué) 2014年9期
    關(guān)鍵詞:約化統(tǒng)計法低質(zhì)量

    GU Jia-yang,YANG Jian-min,XIAO Long-fei

    (1.School of Naval Architecture and Marine Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003,China;2.State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

    Study of Vortex Induced Characteristics of Multi-columns with Low Mass Ratio

    GU Jia-yang1,YANG Jian-min2,XIAO Long-fei2

    (1.School of Naval Architecture and Marine Engineering,Jiangsu University of Science and Technology,Zhenjiang 212003,China;2.State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

    Vortex-induced motion(VIM)of an offshore platform with multi-columns is an area of significant focus in the field ocean engineering.VIM involves the complex fluid-structure interaction problem.In this paper,the Reynolds-averaged Navier-Stokes solver combined with the shear-stress transport k-ω turbulence model for the Navier-Stokes equations was used to simulate the VIM of multi-columns with low mass ratio for different flow-approaching angles and reduced velocities.The values of natural period and dimensionless damping coefficient for columns were taken from the model test in the State Key Laboratory of Ocean Engineering of Shanghai Jiaotong University.A computational grid was set up using GAMBIT software.The coupling fluid-structure interaction between the columns and the flow field was obtained by calculating the instantaneous lift and drag forces acting on the multi-columns caused by the external flow field.The differential motion equation was solved by fourth-order Runge-Kutta method,which was manually written into the user-defined functions,and dynamic mesh technology was then adopted to update flow field.The maximum and nominal statistics methods were used to analyze the amplitude of VIM.When the reduced velocities were small,the ‘jump’ phenomenon was observed in the amplitude curve for 0°flow-approaching angle regardless of whether motion in the stream-wise direction was limited.The maximum amplitude was 0.5D when motion in the stream-wise direction was not limited and was slightly smaller than the maximum amplitude of the limited motion in the stream-wise direction 0.57 D).The restriction of motion in the stream-wise direction and the selection of statistical method significantly affected the size of the ‘lock-in’region.Galloping at the low-frequency and multiple-frequency lock in phenomena was observed for the flow-approaching angles of 0°and 45°.Variation trends of equilibrium position with reduced velocity were very close for the three typical flow-approaching angles.When the reduced velocity was greater than 6.0,the maximum stream-wise amplitude heading was 22.5°,whereas the minimum stream-wise amplitude heading was 0°heading.The vortex shedding pattern of multi-columns showed the 2P mode.

    multi-columns;VIM;dynamic mesh;lock in;low mass ratio

    1 Introduction

    The study of numerical simulation and model test for vortex-induced motion(VIM)started with research on the VIM of a spring-constrained rigid circular cylinder,which is the most representative work by Willamson et al[1-4].Based on this study,the VIM of multi-columns has recently become an area of research focus.The effect on VIM of different configurations,such as serial,parallel,and multi-columns,is emphasized.The spacing ratio of columns is also considered[5-6].

    The primary underwater structures of tension-leg platform(TLP)and semi-submersible platform are similar.Both platforms commonly have four columns of the cross-section type,which primarily includes square and circular shapes.Research on the VIM of offshore platforms with multi-columns can be conducted using two primary methods.The first method is numerical simulation.Commercial software,such as Fluent and CFX,are used to conduct numerical simulation with a focus on the VIM of a single column or four columns.The other method is experimental study.

    Yalpaniyan et al employed the volume-of-fluid method to simulate flow around the TLP with a laminar and turbulent solver[7].The existence of upstream column vortices affects the flow field of downstream columns,which then causes the vortex intensity of upstream columns to be greater than that of the downstream columns.Magee et al discussed model tests for the VIM of a TLP near the southern tip of the Malaysian Peninsula[8].First,decay tests were conducted on the TLP to obtain the values of the natural period and dimensionless decay coefficient.Second,VIM tests primarily studied the effect of flow velocity,current heading,and dimensionless draft(ratio of draft to characteristic length of column)on the VIM of a TLP.Kim et al investigated the VIM of a TLP by employing AcuSolve and STAR-CCM+software,respectively[9].The principal platform dimensions were 57.3 m×57.4 m,and the column diameter was 12.8 m.However,the study only considerered one current heading and three different reduced velocities.Thus,reflecting the vortex-induced characteristics of TLP is difficult.Waals et al conducted a model test for the VIM of four different platforms[10].The platforms were as follows:deep-draft semi-submersible(four pontoons;mass ratio,0.84),deep-draft semi-submersible(two pontoons;mass ratio,0.83),deep-draft TLP(four pontoons;mass ratio,0.57),and deep-draft TLP(four pontoons;mass ratio,0.83).The study found that lower mass ratios,such as those for conventional TLPs,may result in larger sway response.Two pontoons can generate larger motion response at higher Ur(Ur>10).The yaw motion amplitude for the 0°tow direction linearly increased with reduced velocity,and the yaw motion for the 45°tow direction was small(<1°).

    Zhao et al used a finite element method to solve the two-dimensional Reynolds-averaged Navier-Stokes equations and shear-stress transport(SST)k-w turbulence equations[11].Simua-ations were conducted for flow-approaching angles of α=0°,15°,30°and 45°,and the reduced velocities ranged from 1 to 20 with an interval of 1.The results of this research can be referred to the study of the VIM of TLPs in two dimensions.Xu et al have investigated the vortex-induced vibrations of four cylinders in an in-line square arrangement with low Reynolds number using a finite element method[12].The cylinder motion was modeled using a mass-springdamper system.The nonlinear coupling interactions between the cylinder and the flow field were studied through a dynamic mesh model.The mean and fluctuating aerodynamic forces,response characteristics,and vortex shedding pattern in the wake for each cylinder were primarily analyzed with a spacing ratio L/D from 2.5 to 6.0.Multi-columns of TLP move as a whole,which is the primary difference from the VIM work of Xu,considering that their research method has a certain guidance function.

    Rijen et al performed a model test to analyze the effect of current and wave on the VIM of a semi-submersible platform[13].The results showed that the presence of waves delayed the onset of VIM,but did not affect the VIM response.Rijen et al also investigated the effects of hull appendages,mooring system,and external damping on VIM.Ng et al studied the characteristics of flow in the wake region of a vertical bluff cylinder in waves,currents,and combined wave-current flows using a model test[14].Hong et al analyzed the effect of waves on the VIM of a deep semi-submersible platform[15].The research indicated that the sway motion amplitude is evident as VIM occurred.The sway motion amplitude was not only related to currents nor waves.Under a strong current condition,a smaller wave generates greater motion amplitude.

    The VIM of a TLP is more complicated because of the low mass ratio of such platform when compared with a the vortex-induced vibration of a single rise,regardless of mesh optimization,movement zone setup,selection of turbulence model,determination of the relaxation factor,and so on.Recent research on VIM focused on the spar platform.The investigation of VIM on semi-submersible platform and TLP,as types of multi-columns platforms,was the primary focus.Based on the combined self-development program of commercial software,this paper focuses on the vortex-induced characteristics of a TLP.

    2 Numerical computation model

    2.1 Governing equations

    2.1.1 Governing equations of computational fluid dynamics

    The SST k-ω model serves as a turbulent flow model.The transport equation of the turbulence kinetic energy k and its rate of dissipation ω are given by

    2.1.2 Dynamical governing equation in dimensionless forms

    2.2 Computational model

    A multi-column structure is simplified into a spring-mass-damping system.Considering whether motion in the stream-wise direction is limited,two kinds of computational model are established for the 0°flow-approaching angle.For a flow-approaching angle of 22.5°or 45°,only the condition of unlimited motion in the stream-wise direction is considered.The rectangle computational region is 20 D×60 D(D=0.435 m)with 15 D upstream,45 D downstream,and 8 D move zone set up.In this model,mixed grids are used,the region near the column is meshed by an intensively structured mesh,and the remainder of the region adopts an unstruc tured mesh.Fig.1 shows the computational grid for different flow-approaching angles as well as the zoomed-in view of the grid near the column corner.

    The SST k-ω turbulence model and the unsteady,first-order implicit solver are employed for the numerical calculation.The coupled algorithm is used to solve the velocity and pressure coupling equations,whereas the second-order upwind is selected to determine the momentum,turbulence kinetic energy,and rate of dissipation.The flow direction is from left to right,with the left side as the velocity inlet and the right side as the pressure outlet.The relative pressure is set as 0,upper and lower free-slip boundaries are set as symmetrical,and the model surface is set as a wall.

    Fig.1 Computational grid of four columns

    Multi-columns move in the water as a rigid body,and the interaction between multicolumns and the flow field must be considered.The governing equation of fluid,velocity field,and pressure field is solved by Fluent software,through which the lift-and-drag force acting on columns can be acquired.User-defined function is adopted to acquire the force acting on columns,which will then be taken into the dynamical governing equation.Using the fourthorder Runge-Kutta program,the acceleration,velocity,and displacement of columns can be obtained.The column is then provided a new position through the DEFINE_CG_MOTION function,and the flow field parameters are updated.The new flow field consequently affects the movement of multi-columns,thus causing the cycle to continue.Dynamic mesh technology is adopted to realize the movement of columns and to update the flow field.The time step of this simulation is 0.05 s.The new position and force of columns in the current time step are based on the velocity and force from the previous time step.

    3 Natural period and dimensionless decay coefficient

    To determine the oscillation period of multi-columns,a decay model test must be performed.The model test was conducted in the State Key Laboratory of Ocean Engineering of Shanghai Jiaotong University.According to the decay test results in still water,natural periods of surge,sway,heave,roll,pitch,and yaw are 102.4,104.5,2.76,2.37,2.32 and 75.4 s,respectively.Dimensionless damping coefficient of surge,sway,heave,roll,pitch,and yaw are 0.086,0.029,0.007,0.038,0.038,and 0.17,respectively.

    Fig.2 Surge and sway decay curves obtained from model test

    4 Calculation results and analyses

    4.1 VIM response of multi-columns for 0°flow-approaching angle

    This paper adopted two kinds of statistical methods for amplitude statistics.The results are analyzed in a comparative manner,and the statistical methods are as follows:

    Figs.3 and 4 show the time history of VIM at reduced velocities of 3.0,6.0,8.0,11.0 and 14.0 for limited and unlimited motions in the stream-wise direction.The lift-and-drag forces induced by multi-column VIM determined the VIM variation trend in the stream-wise and transverse directions.The vortex-induced characteristics for unlimited motion in the streamwise direction are first analyzed.When the reduced velocity is 3.0,the time history curves of stream-wise and transverse displacement gradually stabilized and became regular after 300 s,which is close to the sine wave.The maximum displacement is approximately 0.17 D,and the minimum displacement is 0.12 D in the stream-wise direction,such that the amplitude is only 0.05 D.In the transverse direction,the amplitude is 0.45 D,which is greater than the amplitude in stream-wise direction.A phase difference exists between the transverse and streamwise time history curves.The stream-wise phase is larger than the transverse phase.With an increase in reduced velocity,the effect of the upstream column on the downstream column gradually strengthens,such that the regularity of the time-history curves for stream-wise and transverse displacement is gradually decreased.When stream-wise motion is limited,the variation trend of transverse amplitude using two different statistical methods remains consistent(Figs.5 and 6),but the size of amplitude differs,and a “jump” phenomenon occurs in the vari ation trend of multi-column amplitude at low reduced velocity.The maximum amplitude of 0.50 D appears at a reduced velocity of 3.5 for two kinds of statistical methods.If the reduced velocity exceeds 4.0,the value of maximum amplitude when adopting the maximum statistical method locks in at 0.47 D.When adopting the nominal response method,the maximum amplitude locks in at 0.31 D.Therefore,with the increase in reduced velocity,transverse amplitude becomes stable.The equilibrium position has a linear relationship with the increase of reduced velocity.At the reduced velocity of 15.0°for 0°heading,the distance that multi-columns deviate from the initial position is 2.87 D.

    Fig.3 Time history of stream-wise and transverse amplitudes at different reduced velocities for unlimited motion in the stream-wise direction

    Fig.4 Time history of stream-wise and transverse amplitudes at different reduced velocities for limited motion in the stream-wise direction

    Vortex-induced characteristics for unlimited motion in the stream-wise direction have been analyzed.According to the numerical simulation of a circular cylinder VIM,the unlimited motion condition can generate greater transverse amplitude.Comparing Figs.3 and 4,the transverse displacement curve is more regular under limited motion in the stream-wise direction,the maximum amplitude for maximum statistical method is 0.57 D,and the result is 0.56 D for the nominal statistical method,indicating that the results are close.When the reduced velocity exceeds 5.0,the value of transverse amplitude for unlimited motion condition is greater than that of the limited motion in stream-wise direction regardless of the statistical method employed,indicating that the result is consistent with the classical theory of the VIM of a circular cylinder.The computational results also verify the correctness and reliability of the fluidstructure interaction program.Nevertheless,results using the maximum statistical method are significantly greater than those using the nominal response method.Regardless of the statistical method employed and regardless of whether the motion in the stream-wise direction is limited,a significant amplitude ‘jump’ phenomenon occurs when the reduced velocity is small,especially when the reduced velocities vary from 2.5 to 3.0 and 3.0 to 3.5.A similar phenomenon was observed when investigating the VIM of a single circular cylinder at a high mass ratio.Given the space constraints,this paper does not provide more time-history curves.However,an apparent“clap” phenomenon is observed in the time-history curve of transverse amplitude when the reduced velocity exceeds 10.0.

    Fig.5 Variation trend of transverse amplitude with reduced velocity for 0°flow-approaching angle

    Fig.6 Time series of lift coefficient on each column and total for 0°heading,U*=3.5

    In this study,multi-columns move as a whole.While modeling using Gambit software,the four columns can be set as a whole wall,and the mechanical characteristics of the platform can be analyzed.However,the method is incapable of analyzing the mechanical characteristics of each column.To analyze the interference problem among multiple columns,the four columns will be set as four walls,but the 8 D move zone associated with multi-columns is set as a whole zone with superimposed hydrodynamic force to facilitate calculation.Fig.6 shows the mechanical characteristics of each column and the total lift force.The lift coefficient curve is no longer a simple regular curve at steady stage because of the coupling effect among the four columns,which is identified with jagged shape at the top and bottom of the lift coefficient curve.The space ratio of columns is 2.9.A greater space ratio causes the coupling effect among the columns to be weaker.The jagged shape of the lift coefficient for upstream and downstream columns is complemented,such that the total lift coefficient curve becomes stable.Comparing the amplitude of the lift coefficient,the lift coefficient amplitude of the downstream column is greater than that of the upstream column.At the same time,the lift coefficient amplitude of downstream column 2 is approximately equal to that of downstream column 3.This rule is also applied to upstream columns 1 and 4.

    Fig.7 FFT analysis of lift coefficient at different reduced velocities for 0°flow approach angle

    Fig.7 presents the fast Fourier transform(FFT)analysis of the lift coefficient at different velocities of 3.5,6.0,8.0,10.0,12.0,and 14.0 for 0°flow-approaching angle.Spectrum transform curves under different reduced velocities are relatively regular except for the reduced velocity of 6.0,which has two spectral peaks,whereas others only have one spectral peak.Thus,the energy of lift coefficient is concentrated under different reduced velocities,and vortex shedding frequency increased with the reduced velocity.Considering that the transverse natural frequency of multi-columns is 0.06 Hz,one of the spectral peaks at the reduced velocity of 6.0 corresponds to the natural frequency of multi-columns,whereas the other corresponds to actual vortex shedding frequency.When the vortex shedding frequency is equal to the natural frequency of multi-columns,the lock-in phenomenon appears.

    Fig.8 Equilibrium position of stream-wise amplitude at different velocities

    Fig.9 Normalized frequency at different reduced velocities for 0°

    Fig.8 shows changing trend of the equilibrium position with different reduced velocities.The equilibrium position increases gradually with incoming flow velocity.However,the slope of the equilibrium position curve for different flow-approaching angles differed;a higher reduced velocity results in a greater slope.At a reduced velocity of 15.0,the equilibrium positions are 2.87,6.38,and 8.10 D,corresponding to 0°,22.5°,and 45°headings,respectively.Under the same reduced velocity,the characteristic length of multi-columns change with the flow-approaching angle.The characteristic length of 45°heading is 1.41 times larger than the size of the 0°heading.The equilibrium position from the initial position is larger than higher mesh quality requirement.

    Fig.9 shows that spectral transform of the lift coefficient for 0°under different reduced velocities.The changing trends for limited and unlimited motion in stream-wise direction are identical.With an increase in reduced velocity,the frequency of the lift coefficient increases linearly.For limited motion condition, ‘lock-in’ phenomenon appears at a reduced velocity of 2.5 to 3.5.The author investigated galloping under the low-frequency phenomenon for a square column at a mass ratio of 2.6.Galloping under the low-frequency phenomenon occurred for limited motion in the stream-wise direction at a reduced velocity of 4.0 to 5.8.Notably,the step-jump phenomena appeared for normalized frequency when the reduced velocity was less than 11.0.The stream-wise motion is limited in Fig.9,and the size of vortex shedding frequency is the same as that of the adjacent reduced velocity.

    4.2 Research on the influence of flow-approaching angle on VIM

    The paper has researched the influence of different flow-approaching angle on the amplitude,region of lock in and equilibrium position in stream-wise direction.Time history of VIM has been given for 22.5°and 45°headings in Figs.10 and 11.The study found that synchronization phenomena evidently emerged in the time-history curve of VIM in the stream-wise and transverse directions at a reduced velocity of 3.0 for 22.5°flow-approaching angle.The peaks and valleys of the time-history curve appeared almost simultaneously.The underlying cause lies in the characteristics of the lift-and-drag coefficient.The transverse amplitude of VIM is 0.43 D,whereas the stream-wise amplitude is 0.39°.Thus,the transverse amplitude is slightly greater than the stream-wise amplitude.The equilibrium position of multi-columns is about 0.31 D.With an increase in reduced velocity,the time-history curve shows a random characteristic.When the reduced velocity further increased to 6.0,8.0 10.0,12.0,and 14.0,the space between transverse and stream-wise time history curves enlarged,multi-columns gradually vibrated farther from the initial position,and randomicity was enhanced.If the maximum statistical method is used,larger errors may occur,which may affect the judgment and understanding of the motion trend.Therefore,the nominal statistical method must be used from the randomicity angle.

    Fig.10 Time histories of stream-wise and transverse amplitudes for 22.5°(unlimited motion in the stream-wise direction)

    Fig.11 Time histories of stream-wise and transverse amplitudes for 45°(unlimited motion in the stream-wise direction)

    Fig.12 Stream-wise amplitude in related with the induced velocity

    Fig.13 Transverse amplitude relative to induced velocity

    Figs.12 and 13 show the changing trend of stream-wise and transverse amplitudes with reduced velocity using the maximum and nominal statistical methods,respectively.The variations of stream-wise amplitude along with reduced velocity are analyzed.The changing trends of the two different statistical results are coincident when comparing Figs.12(a)and(b).Flowapproaching angles significantly affect the transverse amplitude.The amplitude curve is relatively flat at a 0°approaching angle.The threshold value is approximately 4.0 for the maxi-mum statistical method,whereas the nominal statistical threshold value is 5.5.Regardless of the statistical method employed,the transverse amplitude will lock in a valve and will change in the vicinity of this value when the reduced velocity exceeds the threshold value.The value for the maximum statistical method is 0.28 D,which is larger than the nominal statistical value of 0.18 D.When the flow-approaching angle increases to 22.5°,the stream-wise amplitude increases linearly with reduced velocity ranging from 1.5 to 4.0 for both statistical methods.A buffer stage for reduced velocity exists between 4.0 and 6.0.At this stage,transverse amplitude is relatively stable.When the reduced velocity exceeds 6.0,the growth rate of streamwise amplitude significantly increases.When the flow-approaching angle is 45°,variation trend of stream-wise amplitude is similar to that of 22.5°heading but enters a sharp vibration stage faster than other cases.The present studies on the vortex-induced characteristics for the spar platform,semi-submersible platform,and TLP primarily consider the 0°flow approaching-angle condition.Given that the cross section of a spar platform is a circular cylinder,the effect of flow-approaching angle can be ignored.Semi-submersible platform and TLP are types of multi-column platforms,and regardless of whether the section of platform is a circular cylinder or a square cylinder,the flow-approaching angle significantly affects the amplitude of VIM.Researches on the VIM of a single column have shown that transverse amplitude is larger than the stream-wise amplitude.Although VIM is present in the other freedom,the motion response of heaving and other rotational freedom are one to two orders of magnitude smaller than transverse and stream-wise motion,which is also the reason why VIM research considers the two main freedoms.Using the limited computing resources and computing time for these two leading freedoms can also further explain the VIM mechanisms of emergence,evolution,and development.

    At the steady stage for 0°heading,the transverse amplitude is 0.45 D,and stream-wise amplitude is 0.28 D for the maximum statistical method.For the nominal statistical method,the transverse amplitude is 0.32 D,and the stream-wise amplitude is 0.18 D.Regardless of the statistical method employed,the transverse amplitude is 1.7 times larger than the streamwise amplitude,which similar to the conclusion of on the VIM of a single column for 0°heading.

    5 Conclusions

    The VIM of multi-columns has been numerically simulated based on computational fluid dynamics combined with the self-development program of fluid-structure interaction.The effect of limited or unlimited motion in stream-wise direction on VIM has been studied.GAMBIT software was utilized to mesh the model,and a movement zone is set up near the multicolumns,where structured grids near the wall are used.A program in C language was manually written into the user-defined functions and embedded into the solver.Dynamic mesh tech-nology was adopted to update the flow field.Simulations were conducted for flow-approaching angles of 0°,22.5°and 45°.The conclusions are as follows:

    (1)This paper has taken multi-columns with low mass ratio as the research object.The cross sections of the columns are square cylinders.Given the difference of the vortex-induced characteristics between circular cylinder and square cylinder as well as the interaction among the columns,the maximum and the nominal statistical methods were adopted to analyze the data.Regardless of whether stream-wise motion is limited,the “l(fā)eap” phenomenon of transverse amplitudes appeared at 0°heading at a reduced velocity of 2.5 to 3.5.When the stream-wise motion was limited,the maximum transverse amplitude was 0.57 D.For the unlimited streamwise motion condition,the maximum transverse amplitude was 0.53 D.The statistical results of VIM for the two different conditions are very close.When the reduced velocity exceeded 5.0,the magnitude of transverse amplitude for unlimited condition became larger than that of the limited condition.A “l(fā)ock-in” phenomenon occurred under both conditions,but the magnitude of“l(fā)ock-in” amplitude and “l(fā)ock-in” region differed.

    (2)Under the unlimited motion condition,the variation trend of the equilibrium position with reduced velocity was more regular,and the distance between the stream-wise equilibrium position and the initial position for 45°was larger than that of the 22.5°heading and 0°heading at the same reduced velocity.The distance of the 0°flow-approaching angle was the least.When the reduced velocity was lower than 11.0°for 0°,a multiple-frequency “l(fā)ock-in”phenomenon appeared under the unlimited motion condition.At a reduced velocity of 2.5 to 3.0,the phenomenon of galloping at low frequency occurred.

    (3)In terms of low reduced velocity of the three typical flow-approaching angles,the timehistory curve of VIM became more regular,and the phase synchronization phenomenon between lift coefficient curve and drag coefficient curve appeared.As the reduced velocity increased,the randomness of hydrodynamic coefficient curve evidently increased.For 0°flowapproaching angle,the variation trend of transverse and stream-wise amplitudes with reduced velocity changed more gently.For a 22.5°heading,the variation trend of stream-wise amplitude can be divided into three stages:linear increasing stage,buffer stage,and sharp linearly increasing stage.The variation trend of 45°is similar to that of 22.5°but increased more rapidly.Using the nominal statistical method at low reduced velocity,the transverse and streamwise amplitudes of 22.5°heading and 45°heading are very close to each other,with both being larger than the amplitude of 0°heading.

    (4)The vortex shedding pattern of multi-columns with low mass ratio showed the 2P mode.Flow-approaching angle,time and reduced velocity significantly affected the vortex shedding position and vortex shedding pattern.

    [1]Govardban R,Williamson C H K.Modes of vortex formation and frequency response of a freely vibrating cylinder[J].Journal of Fluid Mechanics,2000,420:85-130.

    [2]Jauvtis N,Williamson C H K.Vortex-induced vibration of a cylinder with two degrees of freedom[J].Journal of Fluids and Structures,2003,l7:l035-l042.

    [3]Jauvtis N,Williamson C H K.The effect of two degrees of freedom on vortex-induced vibration at low mass and damping[J].Journal of Fluids Mechanics,2004,509:23-62.

    [4]Govardhan R,Williamson C H K.Defining the modified Griffin plot in vortex-induced vibration revealing the effect of Reynolds number using controlled damping[J].Journal of Fluid Mechanics,2006,561:147-180.

    [5]Sumner D.Closely spaced circular cylinders in cross-flow and a universal wake number[J].ASME,Journal of Fluids Engineering,2004,126:245-249.

    [6]Sumner D,Richards M D,Akosile O O.Two staggered circular cylinders of equal diameter in cross-flow[J].Journal of Fluids and Structures,2005,20:55-276.

    [7]Yalpaniyan A,Goodarzi M.The free surface flow around a TLP[C].Proceedings of the ASME 2010 29th International Conference on Ocean,Offshore and Arctic Engineering,2010,OMAE,20021.

    [8]Allan Magee,Rizwan Sheikh,Kenneth Yeoh Hock Guan.Medel tests for VIM of multi-column floating platform[C].Proceedings of the ASME 2011 30th International Conference on Ocean,Offshore and Arctic Engineering,OMAE,2011-49151.

    [9]Jang-Whan Kim,Allan Magee,Kenneth Yeoh Hock Guan.CFD simulation of flow-induced motions of a multi-column floating platform[C].Proceedings of the 30th International Conference on Ocean,Offshore and Arctic Engineering,OMAE,2011,49437.

    [10]Waals O J,Bultema S.Flow induced motions of multi column floaters[C].Proceedings of the 26th International Conference on Offshore Mechanics and Arctic Engineering,OMAE,2007:669-678.

    [11]Huang Zhiyong.Pan Zhiyuan,Cui Weicheng.Numerical simulation of VIV of a circular cylinder with two degrees of freedom and low mass-ratio[J].Journal of Ship Mechanics,2007,11(1):1-9.

    [12]Zhao Ming,Cheng Liang.Numerical simulation of vortex-induced vibration of four circular cylinders in a square configuration[J].Journal of Fluids and Structures,2012,31:125-140.

    [13]Xu Feng,Ou Jinping,Xiao Yiqing.Numerical study on vortex induced vibrations of four cylinders in an in-Line square configuration[J].Computational Structural Engineering,2009:553-567.

    [14]Rijken O,Leverette S.Experimental study into vortex induced motion response of semi submersibles with square columns[C].Proceedings of the ASME 27th International Conference on Offshore Mechanics and Arctic Engineering,2008,OMAE,57396.

    [15]Jimmy Ng K T,Halkyard J E,Soon C Eng.Statistical characteristics of flow in the wake region of a vertical bluff cylinder in waves,currents and combined wave-current flows[C].OMAE,2010-21181.

    [16]Hong Yongpyo,Choi Yongho,Lee Junyoung,Kim Youngbok.Vortex-induced motion of a deep-draft semi-submersible in current and waves[C].Proceedings of the 18th International Society Offshore and Polar Engineering,ISOPE,2008:453-459.

    低質(zhì)量比多立柱渦激特性研究

    谷家揚1,楊建民2,肖龍飛2

    (1.江蘇科技大學(xué) 船舶與海洋工程學(xué)院,江蘇鎮(zhèn)江 212003;
    2.上海交通大學(xué) 海洋工程國家重點實驗室,上海200240)

    多立柱海洋平臺的渦激運動是非常復(fù)雜的流固耦合問題也是近年來海洋工程領(lǐng)域的研究熱點。文章引入雷諾平均法求解NS方程、結(jié)合SST k-ω湍流模型采用耦合算法對低質(zhì)量比多立柱在不同流向及流速下的渦激運動進行了數(shù)值研究,立柱運動固有周期及無因次阻尼系數(shù)取自上海交通大學(xué)海洋工程國家重點實驗室的相關(guān)實驗數(shù)據(jù)。采用分區(qū)法建立了計算網(wǎng)格,多立柱和流場的流固耦合通過計算外部流場作用于柱體的瞬時升力和拖曳力,將求解運動微分方程的四階Runge-Kutta代碼嵌入UDF求解器中,運用動網(wǎng)格技術(shù)實現(xiàn)流場更新并采用最大值統(tǒng)計法和標(biāo)稱統(tǒng)計法對渦激運動振幅幅值進行統(tǒng)計,研究發(fā)現(xiàn):0°流向較小約化速度時,無論限制流向是否,振幅都存在跳躍現(xiàn)象,不限制流向下橫向最大振幅為0.50 D,略小于限制流向的0.57 D,較大約化速度時,橫向振幅出現(xiàn)了鎖定現(xiàn)象,但是否限制流向和統(tǒng)計方法的選取對鎖定區(qū)間長短有所影響,0°及45°流向時候觀測到了低頻馳振及多頻短暫鎖定現(xiàn)象,三個典型流向下流向平衡位置隨流速的變化趨勢基本相近。不限制流向下約化速度大于6.0時,采用標(biāo)稱統(tǒng)計法得到的22.5°流向橫向振幅最大,0°流向最小,多立柱尾渦結(jié)構(gòu)形式表現(xiàn)為2P模式。

    多立柱;渦激運動;動網(wǎng)格;頻率鎖定;低質(zhì)量比

    P751

    A

    谷家揚(1979-),男,博士,江蘇科技大學(xué)船舶與海洋工程學(xué)院副教授;

    肖龍飛(1973-),男,博士,上海交通大學(xué)研究員。

    1007-7294(2014)09-1055-17

    date:2014-04-01

    Supported by National Natural Science Foundation of China(51309123,51279104);Jiangsu Province Natural Science Research projects in Colleges and Universities(13KJB570002);Jiangsu Province Key Laboratory of Advanced Design and Manufacturing Technology of Ship Open Research Fund Project(CJ1203);Academic Program Development of Jiangsu Higher Education Institutions(PAPD).

    Biography:GU Jia-yang(1979-),male,Ph.D.,associate professor of Jiangsu University of Science and Technology,E-mail:gujiayang@126.com;Yang Jian-min(1958-),male,professor/tutor of Shanghai Jiaotong University;Xiao Long-fei(1973-),male Ph.D.,researcher.

    book=1056,ebook=475

    P751

    A

    10.3969/j.issn.1007-7294.2014.09.005

    book=1070,ebook=489

    book=9,ebook=490

    楊建民(1958-),男,上海交通大學(xué)教授,博士生導(dǎo)師;

    猜你喜歡
    約化統(tǒng)計法低質(zhì)量
    雷人畫語
    美文(2023年5期)2023-03-26 03:15:00
    約化的(3+1)維Hirota方程的呼吸波解、lump解和半有理解
    穩(wěn)健統(tǒng)計法在實驗室能力驗證中的應(yīng)用
    常用統(tǒng)計法處理實驗室間比對結(jié)果的探討
    智能城市(2018年7期)2018-07-10 08:30:02
    低質(zhì)量的婚姻不如高質(zhì)量的單身,是真的嗎?(一)
    婦女生活(2017年5期)2017-05-16 21:04:56
    基于經(jīng)驗?zāi)B(tài)分解/高階統(tǒng)計法實現(xiàn)微機械陀螺降噪
    破解學(xué)前教育低質(zhì)量現(xiàn)象
    閾值隨機共振及其在低質(zhì)量濃度氣體檢測中的應(yīng)用
    M-強對稱環(huán)
    概率統(tǒng)計法在儲量估算中的應(yīng)用
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    国精品久久久久久国模美| 人人妻人人爽人人添夜夜欢视频| 国产一区亚洲一区在线观看| 99热全是精品| 2021少妇久久久久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 如日韩欧美国产精品一区二区三区| 国产精品二区激情视频| 天天影视国产精品| 老汉色∧v一级毛片| 韩国av在线不卡| 国产亚洲一区二区精品| 成年人午夜在线观看视频| 国产精品熟女久久久久浪| 最近手机中文字幕大全| 嫩草影院入口| 9191精品国产免费久久| 国产精品国产av在线观看| 亚洲一区中文字幕在线| 色播在线永久视频| 桃花免费在线播放| 日本一区二区免费在线视频| 亚洲色图 男人天堂 中文字幕| 伊人亚洲综合成人网| 肉色欧美久久久久久久蜜桃| 天天添夜夜摸| 久久久久久久久久久免费av| 精品少妇黑人巨大在线播放| 人人妻人人爽人人添夜夜欢视频| 日本av免费视频播放| 亚洲欧洲日产国产| 97精品久久久久久久久久精品| 亚洲国产精品国产精品| 99精国产麻豆久久婷婷| 菩萨蛮人人尽说江南好唐韦庄| 久久精品久久精品一区二区三区| 久久久久久免费高清国产稀缺| 一级爰片在线观看| 欧美激情高清一区二区三区 | av国产久精品久网站免费入址| 国产亚洲一区二区精品| 国产精品二区激情视频| 最新的欧美精品一区二区| 中文精品一卡2卡3卡4更新| 午夜日韩欧美国产| 91成人精品电影| 黄色一级大片看看| 成人国产麻豆网| 91老司机精品| 黄色 视频免费看| 久久亚洲国产成人精品v| 日日撸夜夜添| 夜夜骑夜夜射夜夜干| 韩国精品一区二区三区| 夫妻性生交免费视频一级片| 99热网站在线观看| 婷婷色综合www| 老熟女久久久| 精品国产露脸久久av麻豆| 国产欧美日韩综合在线一区二区| 精品一区二区三区av网在线观看 | av在线app专区| 成年人免费黄色播放视频| 欧美日韩成人在线一区二区| 国产男女内射视频| 国产精品熟女久久久久浪| 波多野结衣av一区二区av| 国产一卡二卡三卡精品 | 我的亚洲天堂| 菩萨蛮人人尽说江南好唐韦庄| 99热国产这里只有精品6| 男女边摸边吃奶| 大香蕉久久成人网| 国产麻豆69| 男人操女人黄网站| 操美女的视频在线观看| 女人精品久久久久毛片| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品亚洲av一区麻豆 | 精品视频人人做人人爽| 久久毛片免费看一区二区三区| 亚洲中文av在线| 久久国产精品大桥未久av| 免费观看a级毛片全部| 在线观看国产h片| 一二三四在线观看免费中文在| 青春草视频在线免费观看| 成年美女黄网站色视频大全免费| 色综合欧美亚洲国产小说| 久久 成人 亚洲| 老汉色∧v一级毛片| 99re6热这里在线精品视频| 国产精品国产三级专区第一集| 亚洲欧洲日产国产| 亚洲欧美激情在线| 深夜精品福利| av网站免费在线观看视频| 又大又黄又爽视频免费| 99久久综合免费| 久久久久久久大尺度免费视频| 91精品三级在线观看| 青草久久国产| av卡一久久| 久久久久久人妻| 亚洲av日韩精品久久久久久密 | 18禁裸乳无遮挡动漫免费视频| 看免费av毛片| 伊人久久国产一区二区| 大陆偷拍与自拍| 亚洲美女视频黄频| 国产在视频线精品| 成人午夜精彩视频在线观看| 男女边摸边吃奶| 久久ye,这里只有精品| 亚洲精品第二区| 免费看不卡的av| 只有这里有精品99| 日日摸夜夜添夜夜爱| 久久久久久人人人人人| 最近中文字幕2019免费版| 极品人妻少妇av视频| 国产免费一区二区三区四区乱码| 国产一区有黄有色的免费视频| 在线天堂最新版资源| 国产成人欧美在线观看 | 男女无遮挡免费网站观看| 天天操日日干夜夜撸| 国产亚洲欧美精品永久| 国产精品一区二区精品视频观看| 一本一本久久a久久精品综合妖精| 街头女战士在线观看网站| 两个人看的免费小视频| 美女视频免费永久观看网站| 欧美最新免费一区二区三区| 国产精品免费大片| 热re99久久精品国产66热6| 精品人妻熟女毛片av久久网站| 国产熟女午夜一区二区三区| 毛片一级片免费看久久久久| 天天影视国产精品| 欧美日韩成人在线一区二区| 久久韩国三级中文字幕| 涩涩av久久男人的天堂| 伊人亚洲综合成人网| 成人影院久久| 亚洲专区中文字幕在线 | 狂野欧美激情性bbbbbb| 深夜精品福利| 亚洲欧洲国产日韩| 热99国产精品久久久久久7| 久久婷婷青草| 国产一区二区三区av在线| 亚洲精品av麻豆狂野| www.自偷自拍.com| 中文天堂在线官网| 黑人欧美特级aaaaaa片| 亚洲av在线观看美女高潮| 欧美激情高清一区二区三区 | 国产探花极品一区二区| 成人国产麻豆网| 一区二区三区乱码不卡18| 日本欧美视频一区| 波多野结衣一区麻豆| 色播在线永久视频| 少妇被粗大的猛进出69影院| 操出白浆在线播放| 亚洲一区中文字幕在线| 十八禁网站网址无遮挡| 精品午夜福利在线看| 97人妻天天添夜夜摸| 宅男免费午夜| 亚洲av男天堂| 在线 av 中文字幕| 亚洲欧美成人精品一区二区| 别揉我奶头~嗯~啊~动态视频 | av国产精品久久久久影院| 久久久久久人妻| 你懂的网址亚洲精品在线观看| 高清av免费在线| 日韩一区二区三区影片| 成人亚洲欧美一区二区av| 成年女人毛片免费观看观看9 | 超碰97精品在线观看| 日韩熟女老妇一区二区性免费视频| 亚洲伊人久久精品综合| 国产精品一区二区精品视频观看| 岛国毛片在线播放| 日韩av免费高清视频| 精品国产超薄肉色丝袜足j| av有码第一页| 色婷婷av一区二区三区视频| 久久久久久久精品精品| 午夜91福利影院| 热re99久久国产66热| 99久久综合免费| 水蜜桃什么品种好| 999久久久国产精品视频| 亚洲人成网站在线观看播放| 捣出白浆h1v1| 亚洲欧美精品自产自拍| 看非洲黑人一级黄片| 伦理电影免费视频| 美女福利国产在线| 亚洲国产欧美日韩在线播放| 国产一区有黄有色的免费视频| 中文字幕人妻丝袜一区二区 | 国产片特级美女逼逼视频| 免费黄网站久久成人精品| 午夜免费观看性视频| 天天添夜夜摸| 久久午夜综合久久蜜桃| 在线亚洲精品国产二区图片欧美| 国产免费福利视频在线观看| 一区二区三区四区激情视频| 美女中出高潮动态图| 成人手机av| 精品卡一卡二卡四卡免费| 国产xxxxx性猛交| 一区二区三区精品91| 99久久99久久久精品蜜桃| 久久婷婷青草| 波多野结衣一区麻豆| 亚洲精品国产av成人精品| 两个人免费观看高清视频| 日本色播在线视频| 国产 精品1| 久久鲁丝午夜福利片| 日韩大片免费观看网站| 老汉色∧v一级毛片| 狠狠婷婷综合久久久久久88av| 男女高潮啪啪啪动态图| 男女午夜视频在线观看| 国产精品一国产av| www日本在线高清视频| 尾随美女入室| 亚洲av福利一区| 国产日韩欧美视频二区| 你懂的网址亚洲精品在线观看| 久久久久人妻精品一区果冻| 亚洲av在线观看美女高潮| 久久免费观看电影| 天天躁日日躁夜夜躁夜夜| 两个人免费观看高清视频| 日韩伦理黄色片| 五月天丁香电影| 国产探花极品一区二区| √禁漫天堂资源中文www| 久久久久久人妻| 伊人亚洲综合成人网| 黄色毛片三级朝国网站| 亚洲欧美色中文字幕在线| 亚洲精品,欧美精品| 王馨瑶露胸无遮挡在线观看| 国产精品二区激情视频| 国产亚洲午夜精品一区二区久久| av在线app专区| 超色免费av| av电影中文网址| 蜜桃在线观看..| 日韩制服骚丝袜av| 国产成人欧美| 日韩一本色道免费dvd| 久久韩国三级中文字幕| 深夜精品福利| 免费av中文字幕在线| 午夜精品国产一区二区电影| 这个男人来自地球电影免费观看 | 精品福利永久在线观看| 啦啦啦中文免费视频观看日本| 精品一区二区三卡| av福利片在线| 国产亚洲一区二区精品| 精品福利永久在线观看| 最近手机中文字幕大全| 精品少妇黑人巨大在线播放| 2018国产大陆天天弄谢| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 在现免费观看毛片| 成人亚洲欧美一区二区av| 亚洲国产精品一区三区| 青春草视频在线免费观看| 日本一区二区免费在线视频| 另类精品久久| 日韩人妻精品一区2区三区| 大片免费播放器 马上看| 国产精品香港三级国产av潘金莲 | 如日韩欧美国产精品一区二区三区| 国产精品女同一区二区软件| 伊人久久国产一区二区| 天天操日日干夜夜撸| 欧美少妇被猛烈插入视频| 精品少妇久久久久久888优播| 韩国高清视频一区二区三区| 日韩,欧美,国产一区二区三区| 亚洲天堂av无毛| 亚洲精品久久午夜乱码| 青春草国产在线视频| 黑人巨大精品欧美一区二区蜜桃| 老司机亚洲免费影院| 国产野战对白在线观看| 啦啦啦 在线观看视频| 日本av免费视频播放| 91老司机精品| 伊人久久大香线蕉亚洲五| 性高湖久久久久久久久免费观看| 国产伦人伦偷精品视频| 少妇被粗大猛烈的视频| 欧美乱码精品一区二区三区| 国产xxxxx性猛交| 婷婷色麻豆天堂久久| 人人妻人人添人人爽欧美一区卜| 婷婷成人精品国产| 国产在线免费精品| 国产欧美日韩综合在线一区二区| 午夜91福利影院| 人妻 亚洲 视频| 国产 一区精品| 色播在线永久视频| 亚洲av综合色区一区| 女性生殖器流出的白浆| 99精品久久久久人妻精品| 中国三级夫妇交换| 热re99久久国产66热| 日本91视频免费播放| 国产成人av激情在线播放| 麻豆精品久久久久久蜜桃| 精品人妻在线不人妻| 亚洲av电影在线进入| 最新在线观看一区二区三区 | 性高湖久久久久久久久免费观看| av福利片在线| 黑人欧美特级aaaaaa片| 侵犯人妻中文字幕一二三四区| 国产亚洲精品第一综合不卡| 欧美日本中文国产一区发布| 日本色播在线视频| 啦啦啦中文免费视频观看日本| 国产无遮挡羞羞视频在线观看| 亚洲第一av免费看| 久久免费观看电影| 中文字幕色久视频| 国产成人a∨麻豆精品| 新久久久久国产一级毛片| 成年美女黄网站色视频大全免费| 一级毛片电影观看| 麻豆精品久久久久久蜜桃| 亚洲欧美成人精品一区二区| 制服丝袜香蕉在线| 国产亚洲欧美精品永久| 国产精品嫩草影院av在线观看| 欧美变态另类bdsm刘玥| 亚洲美女视频黄频| a 毛片基地| 日韩电影二区| av.在线天堂| 一级片'在线观看视频| 亚洲熟女毛片儿| 国产毛片在线视频| 国产精品国产三级专区第一集| 国产日韩一区二区三区精品不卡| 久久久精品区二区三区| 黄片小视频在线播放| 国产成人免费观看mmmm| 国产日韩欧美视频二区| 欧美日本中文国产一区发布| 日韩伦理黄色片| 亚洲婷婷狠狠爱综合网| 久久人人97超碰香蕉20202| 亚洲精品一二三| 在线观看免费日韩欧美大片| 男女国产视频网站| 中文字幕亚洲精品专区| 捣出白浆h1v1| 视频区图区小说| 日本欧美国产在线视频| 成人午夜精彩视频在线观看| 又大又黄又爽视频免费| 最近的中文字幕免费完整| 色网站视频免费| 日韩大片免费观看网站| 国产黄色视频一区二区在线观看| a级片在线免费高清观看视频| 亚洲一卡2卡3卡4卡5卡精品中文| 我的亚洲天堂| 国产精品久久久av美女十八| videosex国产| 视频区图区小说| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 国产欧美日韩综合在线一区二区| 亚洲 欧美一区二区三区| 久久精品国产a三级三级三级| 亚洲七黄色美女视频| 青春草亚洲视频在线观看| 在线亚洲精品国产二区图片欧美| 成年美女黄网站色视频大全免费| 国产日韩一区二区三区精品不卡| 国产在线一区二区三区精| a级毛片在线看网站| 国产成人免费无遮挡视频| 久久久精品免费免费高清| 免费在线观看黄色视频的| 精品国产国语对白av| 亚洲精品久久午夜乱码| 婷婷色av中文字幕| 大香蕉久久网| 亚洲国产中文字幕在线视频| 欧美乱码精品一区二区三区| 欧美激情 高清一区二区三区| 少妇人妻久久综合中文| 最近手机中文字幕大全| 欧美黄色片欧美黄色片| 国产精品一国产av| 精品国产国语对白av| 亚洲国产成人一精品久久久| 久久国产精品大桥未久av| 免费久久久久久久精品成人欧美视频| 亚洲国产欧美日韩在线播放| 母亲3免费完整高清在线观看| 欧美精品av麻豆av| 精品少妇一区二区三区视频日本电影 | 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| xxxhd国产人妻xxx| 午夜影院在线不卡| 日韩不卡一区二区三区视频在线| av女优亚洲男人天堂| 精品少妇黑人巨大在线播放| 美女扒开内裤让男人捅视频| 国产爽快片一区二区三区| 欧美日韩亚洲综合一区二区三区_| 日韩中文字幕视频在线看片| 黄色 视频免费看| 51午夜福利影视在线观看| 亚洲欧洲精品一区二区精品久久久 | 色婷婷久久久亚洲欧美| 午夜福利免费观看在线| 亚洲精品久久午夜乱码| 欧美97在线视频| 操出白浆在线播放| 国产精品一区二区精品视频观看| 欧美人与善性xxx| 又粗又硬又长又爽又黄的视频| 激情视频va一区二区三区| 国产av码专区亚洲av| 丁香六月欧美| 欧美黑人欧美精品刺激| 一边亲一边摸免费视频| 高清视频免费观看一区二区| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 国产av码专区亚洲av| 国产精品一区二区在线不卡| 精品亚洲成国产av| 菩萨蛮人人尽说江南好唐韦庄| 日韩免费高清中文字幕av| 久热这里只有精品99| 日本一区二区免费在线视频| 欧美 亚洲 国产 日韩一| 久久天堂一区二区三区四区| 亚洲欧洲国产日韩| 久久久精品区二区三区| 亚洲欧美精品自产自拍| 亚洲av电影在线观看一区二区三区| 少妇精品久久久久久久| 日韩不卡一区二区三区视频在线| 777久久人妻少妇嫩草av网站| 久久女婷五月综合色啪小说| 观看美女的网站| 久久狼人影院| e午夜精品久久久久久久| 亚洲精品成人av观看孕妇| 精品国产乱码久久久久久男人| 一二三四中文在线观看免费高清| 国产淫语在线视频| 波多野结衣一区麻豆| 久久久久精品性色| 天堂中文最新版在线下载| 久久精品国产综合久久久| a级片在线免费高清观看视频| 中国三级夫妇交换| 国产深夜福利视频在线观看| 一区二区三区精品91| 老司机影院成人| 久久鲁丝午夜福利片| 欧美中文综合在线视频| 国产精品麻豆人妻色哟哟久久| 观看av在线不卡| 中文乱码字字幕精品一区二区三区| 性高湖久久久久久久久免费观看| 久久人人爽av亚洲精品天堂| 国产精品国产av在线观看| 国产亚洲午夜精品一区二区久久| 久久精品国产a三级三级三级| 日韩大码丰满熟妇| 在现免费观看毛片| 尾随美女入室| 国产欧美日韩综合在线一区二区| 免费黄色在线免费观看| 国产精品av久久久久免费| 观看美女的网站| 秋霞伦理黄片| 色吧在线观看| 男女午夜视频在线观看| 极品人妻少妇av视频| 国产精品蜜桃在线观看| 中文天堂在线官网| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 国产男女内射视频| 午夜影院在线不卡| 亚洲综合色网址| 亚洲欧美色中文字幕在线| 99久久精品国产亚洲精品| 18禁观看日本| 久久99精品国语久久久| 国产精品99久久99久久久不卡 | 日韩成人av中文字幕在线观看| 男人添女人高潮全过程视频| 国产成人免费观看mmmm| 欧美 日韩 精品 国产| 男女边摸边吃奶| 9热在线视频观看99| 中文字幕人妻丝袜制服| 亚洲自偷自拍图片 自拍| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 久久久欧美国产精品| 国产精品一区二区在线不卡| 少妇被粗大的猛进出69影院| 欧美中文综合在线视频| 亚洲精品久久成人aⅴ小说| 久久这里只有精品19| 亚洲专区中文字幕在线 | 妹子高潮喷水视频| 天天躁夜夜躁狠狠久久av| 免费在线观看视频国产中文字幕亚洲 | 色吧在线观看| 18禁动态无遮挡网站| 国产一区二区三区av在线| 啦啦啦啦在线视频资源| av在线观看视频网站免费| 国产亚洲av高清不卡| 天天操日日干夜夜撸| 热99国产精品久久久久久7| 国产熟女欧美一区二区| 欧美日韩国产mv在线观看视频| 久久女婷五月综合色啪小说| 男男h啪啪无遮挡| 国产免费现黄频在线看| 美女福利国产在线| 国产免费现黄频在线看| 美女福利国产在线| 欧美最新免费一区二区三区| 国产免费一区二区三区四区乱码| 亚洲av电影在线进入| 好男人视频免费观看在线| 久久久国产精品麻豆| 丝袜美足系列| 午夜福利在线免费观看网站| 一二三四中文在线观看免费高清| 国产欧美亚洲国产| 黄色怎么调成土黄色| 国产毛片在线视频| 日本av手机在线免费观看| 国产伦理片在线播放av一区| 一边亲一边摸免费视频| 亚洲国产日韩一区二区| 在线 av 中文字幕| 别揉我奶头~嗯~啊~动态视频 | 18禁动态无遮挡网站| 中文字幕色久视频| 亚洲国产精品999| 好男人视频免费观看在线| 免费观看人在逋| 99热全是精品| 桃花免费在线播放| 免费在线观看完整版高清| 美国免费a级毛片| 在线观看免费日韩欧美大片| 亚洲成人免费av在线播放| 成人影院久久| 亚洲人成77777在线视频| 国产精品99久久99久久久不卡 | 天天影视国产精品| 日本av手机在线免费观看| 制服丝袜香蕉在线| 丝袜在线中文字幕| 在线观看三级黄色| 成年人免费黄色播放视频| 日本欧美国产在线视频| 哪个播放器可以免费观看大片| 国产又爽黄色视频| 亚洲色图综合在线观看| 最近最新中文字幕大全免费视频 | 日本黄色日本黄色录像| 在线 av 中文字幕| 在线观看免费午夜福利视频| 久热这里只有精品99| 色播在线永久视频| 国产精品国产三级专区第一集| 午夜精品国产一区二区电影| 男女边吃奶边做爰视频| 久久精品国产亚洲av涩爱| 无限看片的www在线观看| 午夜av观看不卡| 在线观看一区二区三区激情| 国产亚洲午夜精品一区二区久久| 看免费成人av毛片| 97人妻天天添夜夜摸| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 我的亚洲天堂| 永久免费av网站大全|