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

    A New Method for Fast Estimation of Collision Dissipated Energy and Its Application in Ship-Platform Collision Accident

    2018-01-04 08:25:05GUJiayangDENGBinglinTAOYanwuWANGPuZHOUJia
    船舶力學(xué) 2017年12期
    關(guān)鍵詞:海洋工程碩士研究院

    GU Jia-yang,DENG Bing-lin,TAO Yan-wu,WANG Pu,ZHOU Jia

    (1.Jiangsu University of Science and Technology,Zhenjiang 212003,China;2.Marine Design&Research Institute of China,Shanghai 200011,China)

    A New Method for Fast Estimation of Collision Dissipated Energy and Its Application in Ship-Platform Collision Accident

    GU Jia-yang1,DENG Bing-lin1,TAO Yan-wu1,WANG Pu2,ZHOU Jia2

    (1.Jiangsu University of Science and Technology,Zhenjiang 212003,China;2.Marine Design&Research Institute of China,Shanghai 200011,China)

    In the course of operation,the offshore platform has the potential risk of collision with ships.Once an accident occurs,it may cause serious damages to the structure of the platform and affect the safety of operation.Therefore,more attention should be paid to the crashworthiness of the offshore platform in the design stage.In this paper,collisions between an ocean supply ship and a semi-submersible drilling platform are studied.Based on the collision external dynamic theory and DNV specification calculation method,taking into account the influence of yawing motion response to platforms,a new dissipated energy estimation method,named DNV`R method,is proposed.Using the numerical simulation results of nonlinear finite element of ship-platform collision for reference,this analytic method is compared with the specification method,Zhang’s method and Liu’s method.The comparison shows that the DNV`R method fills the gap of the limitation of the specification method that is conservative and only applicable for central impact.The DNV`R method is simple and convenient.And its calculation results are close to those of Zhang’s method or those of Liu’s method,which meets the accuracy requirements of fast estimation of ship-platform collision dissipated energy.

    ship-platform collision;dissipated energy;analytical method;external dynamics of collision

    0 Introduction

    The offshore platform has the potential risk of collision with ships in the course of operation.The resultant structural damage is likely to bring about economic losses,and may even cause serious consequences to people’s life and property and ecological environment,so it is necessary to consider the collision resistance of the platform against ship impacts in the earliest design stages of an platform structure.In terms of energy,the nature of the collision is theprocess of energy conversion:the strain energy generated by the platform structural damage derives from the energy dissipated by the impact ship,in other words,the magnitude of the strain energy due to a collision accident reflects the severity of the structural damage.Therefore,it is important to estimate the energy dissipated by the collision vessel accurately and quickly as it can provide important information for the crashworthiness design of the platform structure.

    For collisions,it has strong nonlinear characteristic,and the collision process is accompanied by the large deformation of the structure and the coupling reaction of the fluid around the two sides.It is difficult to consider the above factors at the same time in the actual analysis.Wang’s[1]research shows that the collision process of ships can be split into collision external dynamics and internal dynamics,respectively.The former mainly analyzes the movement of the hull and the energy loss after the collision,and the latter analyzes the deformation resistance,deformation energy dissipation and structural damage of the collisional structure based on the theory of elastoplastic theory.

    In the study of external dynamics,Minorsky[2]pioneered the method of estimating the kinetic energy dissipation in the ship collision.Petersen[3]analyzed the external mechanism of ship collision by time domain simulation.The hydrodynamic force of ship was calculated by slicing theory.It was assumed that the main body of the ship was rigid body and the deformation caused by collision mainly occurred on the contact surface,and the contact surface was simplified as a nonlinear spring.Woison[4]analyzed the energy loss of the ship collision and assumed that the collided ship was completely plastic,the hull surface was completely rough,and no relative slippage occurred on the collision contact surface.Besides,Zhang[5]developed an algorithm for the external mechanism of plane motion of ship-ship collision,which solved the energy dissipation and impulse change value in each direction by integrating the contact force and relative displacement in each direction.Based on Stronge’s work[6]that developed an advanced solution for 3D impacts,Liu[7]extended the method of Zhang’s[5]to three dimensions,developed a new formulation for the analysis of the impact mechanics of ship collision and applied the formula to the ship-iceberg impact field.Besides,Hu et al[8]analyzed the external dynamics of the ship-ship collision model,and proposed a method for determining the coefficient of conservative restitution.

    In the internal dynamics study,McDermott et al[9]established a set of low-energy collision analysis methods for oil tankers.The method is based on a series of basic assumptions that include bending and buckling of stiffened plates,film stretching of outer panels and decks,and failure of frames.In addition,Hairs[10],Gao et al[11]and Liu[12-13]studied the damage deformation and failure of the components in the collision,and established the corresponding calculation model.

    In recent years,some researchers have coupled the external dynamics and the internal dynamics collision analysis[14-15].Compared to the analysis of two individual processes,it includes more aspects and provides a relatively comprehensive revelation of the basic rules of the collision.However,the cost of coupling analysis is high,and independent external analysis and internal mechanism analysis are still efficient and the analytical results obtained are still within a certain permit range of analytical methods.

    For the ship-platform collision,due to the relatively smaller probability of ship-bridge collision and ship-ship collision and the relatively low attention,it is studied mainly on the research theory of ship-ship collision and ship-bridge collision at present.The method of calculating the dissipated energy after the ship impact platform is divided into analytic method and numerical method,i.e.,finite element method.Based on the general theory of rigid body collision mechanics,the mathematical equation deduced by analytical method is simple,easy to apply and widely applicable.And the numerical method requires accurate mechanical models,including meshing,boundary condition and external force processing.Although the result is more accurate,it requires a lot of modeling and calculation cost input,and only applies to a specific mechanical model,which means the applicability is limited.

    In this paper,collisions between an ocean supply ship and a semi-submersible drilling platform are studied.Based on the external dynamic theory of collision and DNV specification[16],the dynamic analysis of the ship-platform eccentric impact scenario is carried out.Considering the influence of yawing motion response to a platform,a new dissipative energy estimation method is developed,called DNV`R analytical method.Furthermore,a calculation example of collision scene with different impact velocity and different impact angle is carried out.Referring to the numerical simulation results of nonlinear finite element of ship-platform collision the result of DNV`R analysis method is compared to those of the standard method,Zhang’s[5]analysis method and Liu’s[7]analysis method.

    1 Calculation of dissipated energy by DNV specification method

    Based on the characteristics of energy conversion during ship-platform collision,DNV specification[16]proposed the energy dissipation experience of floating platform.When calculating the impact energy dissipation of the collision vessel,it is generally believed that the kinetic energy dissipated by the impact vessel is all transformed into the structural strain energy,and the collision platform reaches the common velocity.With the energy conservation equation and the momentum conservation equation of the collision system,the dissipated energy calculation formula can be obtained as follows:

    where Edisiptis the dissipated energy of the striking ship,ms,miare the mass of the ship and platform,as,aiare the added mass of the ship and platform,and vs,viare the velocity before collision of the ship and platform.

    In the operation process on the platform,it is generally considered that the speed of the platform before the collision is vi=0,then Eq.(1)can be reduced to:

    2 The dissipated energy estimation method considering the influence of yawing motion response to platform

    Fig.1 A collision scenario with an offset between a striking ship and a semi-submersible drilling platform

    With the semi-submersible drilling platform for the collision object,the platform column is the impact site.Assuming that the ship-platform collision is entirely plastic collision,and the struck platform is unconstrained before the collision,considering that the speed of the platform before the collision is vi0=0 and the yawing motion response to platform,we develop a new analytical method for estimating dissipated energy.A collision scenario with an offset between a striking ship and a semi-submersible platform is shown in Fig.1.

    Assume the mass of a striking ship is Ms,the ship sails at a forward speed of vs0,the added mass coefficient is ms,and the velocity at the end of the collision is vs1.The mass of the struck platform is Mi,the added mass coefficient is mi,the velocity in the COG of the platform at the end of the collision is vi1,and the final angular velocity is wi1.The vertical distance between the direction of the impact and the COG of the platform is R,the moment of inertia of semi-submersible drilling platform around its COG including added inertia is Ji,and the radius of gyration of semi-submersible drilling platform is ρi,i.e.,As the collision duration is short,and the collision force is very large,so all other external forces are neglected,the momentum of the collision system is conserved during the collision.This leads to the two following equations:

    At the critical state of the end of the collision,assuming that the striking ship and the platform reach a common speed,that is:

    According to Eq.(3),Eqs.(4)and(5),we can solve the speed at the end of the platform impact:

    Similarly,the speed of the striking ship at the end of the collision can be solved:

    On the basis of the conservation of energy,the dissipation energy of the impact ship at the end of collision is:

    By substituting Eq.(5),Eq.(7)and Eq.(8)into Eq.(9),we can obtain:

    Note that the expression of the moment of inertia of platform,i.e.,with the def-inition of p and q,we can obtaincan be expressed as:

    It can be seen from the Eq.(11)that the energy dissipated by the impact vessel can be quickly estimated if the ratio of the impact ship and the platform,p,the vertical distance between the impact direction and the center of gravity of the platform,R,and the radius of the platform ρiare known.

    In order to facilitate the comparison with the following examples,the method of estimating dissipative energy of Eq.(2)is denoted as DNV method,and Eq.(11),the new formulation since the eccentricity R is taken into account,is named DNV`R method.It can be seen from the two formulas that the formulas of the two methods have certain similarity,i.e.,when the eccentricity is 0,the Eq.(11)can be reduced to Eq.(2).

    3 Principal dimensions and the FE model of platform and OSV

    Due to the lack of the calculation examples of the dissipated energy of the ship platform collision and the relevant experimental data,in this paper,the striking ship is chosen as the supply vessel,and the struck platform is a semi-submersible drilling platform.The nonlinear finite element analysis(NLFEA)method is adopted based on MSC.DYTRAN nonlinear analysis software,so as to realize the virtual collision test between OSV and selected platform,and the calculated results are used as a reference object of the estimation results of each analytic method.

    3.1 Principal dimensions of platform and OSV

    The semi-submersible drilling platform is 3 000 m of operating depth,whose principal dimensions are as shown in Tab.1.The selected striking ship is an ocean supply ship(OSV)that is the main service vessel of the platform.The principal dimensions of the OSV are presented in Tab.2.

    Tab.1 Principal dimensions of semi-submersible drilling platform

    Tab.2 Principal dimensions of OSV

    3.2 Collision scenarios between OSV and platform

    There are three conditions of ship collision platform,i.e.,bow impact,stern impact and side impact.According to the data statistics of IACS about collision accidents between supply ships and platforms in the past,we can find that the probability of stern impact is 70%,theside impact 20%,and the bow impact very unlikely.Considering the characteristics of the selected OSV whose accommodation space is located at the bow and with a big deck for cargo laid out at the stern,the berthing way of OSV is mainly stern toward the target platform.From an integrative consideration,it is reasonable and necessary to choose the stern impact as the main type of collision.When OSV is berthing stern to target platform,the possible impact position of the supply ship on the work platform is the straight wall or the transition surface of the column,as shown in Fig.2.

    3.3 The FE model of platform and OSV

    The FE model of semi-submersible drilling platform and OSV is shown in Fig.3.The element size of the collision zone is generally 160 mm.For large offshore structures in full-scale collision simulations,the mesh size can accurately simulate the damage and deformation of structures and capture the strain energy absorption values of the plastic deformation structure[17].At the meantime,for the structures that far from the collision zone,larger elements are adopted not only to accurately simulate the impact of the inertia of the OSV and a semi-submersible drilling platform in the process of the collision,but to reduce modeling and calculation costs.

    Fig.3 The FE model of platform and OSV

    3.4 Main parameter setting in numerical simulation of ship-platform collision

    To be as close as possible to the actual collision scene,all of the structures of the striking ship and platform are made of the elastic-plastic material,which is ElasPlas(DAMTEP24)in the DYTRAN code.The OSV uses low-carbon steel material parameters,whose parameters of strain rate are obtained by D=40.4/s and q=5;correspondingly,high strength steel material parameter is adopted in the platform and the parameters of strain rate are obtained by D=4 000/s and q=5[18].Moreover,the failure strain of these two materials are set to equal to εf=0.15.

    The collision process of ship-platform involves the dynamic contact problem.The key technique of solving is to deal with the contact algorithms of different structure interface.Penalty-based contact algorithms are used to simulate the contact between the striking ship and the column of the platform and the internal of either the ship or the platform.Besides,the static and dynamic friction coefficients of 0.3 are used in all of the contact formulations[19].Furthermore,when the platform is in operation,it is a typical floating platform[16].Therefore,in numerical simulation,the struck platform is free of any boundary conditions,so that it is in a fully free state,in order to make numerical simulation closer to the actual collision scenario.

    4 Numerical examples

    In this paper,when the analytical method is analyzed,Zhang’s analytic method and Liu’s analytic method in the ship collision dissipation energy estimation are applied to the ship-platform collision scene.The results are compared with those of the DNV`R analytic method and the new formulation proposed in this paper.In other words,DNV,DNV`R,Zhang,Liu`2D and Liu`3D dissipative energy analysis are used to analyze the cases in the following numerical examples.On using the Liu`2D method,it is assumed that the vertical position of the collision point is 0,i.e.,Zsc=Zic′=0.In addition,taking into account the complexity of the computation process on using Zhang’s method and Liu’s methods,we will work out programs to realize the dissipation energy of ship-platform collision scenarios respectively.Finally,the results of nonlinear finite element numerical simulation of ship-platform collision are compared as a reference with those of the analytical methods.

    On calculating by the analytic method,the coefficient of restitution is assumed to e=0,i.e.,that the ship-platform collision is entirely plastic collision.And taking into account the interaction effects between the striking ship and the struck platform and the surrounding water,the additional water mass method is adopted,and the added mass coefficients obtained are shown in Tab.3.The radius of gyration of them can be determined by empirical formula,as presented in Tab.4.In the table,Liis the length of the lower hull of the platform,Biis the width of the outer side of the two lower hulls,and Lsand Bsare the length and the width of the supply ship respectively.

    Tab.3 Hydrodynamic coefficients

    4.1 The supply vessel collision with the semi-submersible platform along the Y′direction

    The supply vessel has the same impact position,and the impact velocity is chosen from 0.5 m/s to 10 m/s.At the end of the collision,the total strain energy absorbed by the structure is the calculation result of Nonlinear Finite Element Analysis(NLFEA)when the impact force is decreased to 0.Meanwhile,according to the parameters such as mass ratio p,eccentricity R and impact velocity of the ship and the platform,the dissipated energy is estimated by different analytic methods.The results are shown in Tab.5.

    Tab.4 The radius of gyration

    Tab.5 The dissipated energy of different impact velocity at the same impact position

    From Tab.5,the coefficient of energy dissipation ηdisipt,and relative error of each calculation method can be obtained,ηdisiptis defined as follows:

    where Edisiptis the dissipated energy of OSV,Ekis the initial kinetic energy of the OSV.

    And the relative error of analytic method,Error,is defined as:

    where EAnalyticalMethodis the dissipated energy estimated by the analytic method,ENLFEAis the strain energy calculated by NLFEA.

    Fig.4 is a graph showing the coefficient of energy dissipation of each calculation result obtained by Eq.(12),and Fig.5 is a graph showing the relative errors of various analytical methods with the impact velocity based on Eq.(13).

    Fig.4 The coefficient curve of energy dissipation and impact velocity relationship

    Fig.5 The relative errors of various analytical methods for impact velocity

    From Tab.5 and Fig.4,it is clear that the dissipation energy results calculated by each analytical method are larger than those obtained by NLFEA.This is mainly due to the fact that on using NLFEA,the partial of the initial kinetic energy of the striking ship is transformed into the hourglass energy of the shell element and the friction energy,so the strain energy dissipation transformed by kinetic energy of the ship is relatively small.

    It can be seen from Fig.5,compared to NLFEA strain energy,that DNV results have greater errors,and the relative error is enlarging gradually with the increase of the impact velocity.Specifically,the error is increased from about 14%when the impact velocity is 0.5 m/s to about 22%when the velocity is 10 m/s,thus,the dissipative energy results of DNV estimation are relatively conservative.In addition,the main difference between Zhang and Liu`2D is the coordinate system when the position parameter is input.When the vertical height of the impact point is neglected,the calculated results are consistent with NLFEA results given by MSC.DYTRAN.The results show that it is feasible to apply Zhang’s method and Liu’s method to estimate the dissipation energy of the ship-platform collision.More importantly,the calculation results of DNV`R method developed in this paper are very close to those of NLFEA,and to those calculated by Zhang’s method and Liu’s method.In particular,for the same impact position,the smaller the impact speed is,the more accurate the estimated result of energy is,and the maximum relative error is only 8.2%when the impact speed is 10 m/s.

    4.2 The supply vessel impacts the semi-submersible platform with different impact angles

    In the actual collision accident,the impact direction of the supply ship impact platform has certain randomness due to human factors and sea conditions and many other uncertainties,so it is necessary to consider the collision scenarios of the ship’s stern at different angles against the semi-submersible column.The analysis shows that due to the particularity of the outer surface of the platform column,each impact point is perpendicular to the longitudinal direction of the supply ship in the tangential direction of the collision plane,i.e.,the waterline angles α,are 90°.In order to analyze the impact of different impact angles θ,on the calculation results,the 4 m/s impact velocity and the same height over the COG of the platform are analyzed.A schematic diagram of the supply vessel against the platform column with different impact angles,for simplicity the longitudinal direction of the supply vessel is indicated by the longitudinal profile of the striking ship,is shown in Fig.6.

    With the collision parameters of different impact angles,the dissipation energy of the supply ship can be obtained by the analytic methods(the initial impact kinetic energy of the supply vessel is 44MJ),and the calculation results of energy dissipation corresponding to impact angles are shown in Tab.6.

    Tab.6 The calculation results of dissipation energy of supply vessel

    According to Eq.(12)and Eq.(13),it can be calculated for each calculation of the dissipation energy coefficient curve and the relative error with the impact velocity curve respectively,as shown in Fig.7 and Fig.8.

    Fig.7 The coefficient curve of energy dissipation and impact angle relationship

    Fig.8 The relative errors of various analytical methods for impact angle

    It can be seen from Fig.7 that the dissipation energy coefficient of each method is decreasing as the impact angle increases at the same impact velocity,indicating that the supply ship dissipative energy is declining,and the ship-platform motion response is increasing with the impact angle increasing.Besides,as we can see from Fig.8,the results of analytical methods have some deviations compared with numerical simulation.This kind of error is mainly caused by the existence of hourglass and other damping in numerical simulation,which is unavoidable.And we can find that the results obtained by Liu`2D and Liu`3D are very similar,indicating that the height of the collision is not significant for the roll motion of the platform.Finally,it is important to know that by comparing the relative error of each analytic method,it can be seen that the error of DNV`R method is basically close to that of Zhang’s method and Liu’s method.

    4.3 The supply vessel impacts the semi-submersible platform column transition at different angles

    Fig.9 The supply vessel impacts the platform in the transition direction with different impact angles

    The supply of the ship stern may also impact the platform column transition,supply ship stern impact point corresponding to the platform transition for the cylindrical straight wall,the corresponding water angle is 90°,the collision scene is shown in Fig.9.

    Similar to Section 4.2,take different impact angle,4 m/s impact speed and the same impact height for calculation and analysis respectively.With inputting the collision parameters of each impact angle,the dissipation energy of the supply ship can be obtained by each analytic method,see Tab.7.

    Tab.7 The dissipated energy of OSV impacting the platform along the column transition

    According to the Eq.(12)and Eq.(13),it can also be calculated for each calculation result of the dissipation energy coefficient curve and the relative error with the impact velocity curve respectively,as shown in Fig.10 and Fig.11.

    Fig.10 The coefficient curve of energy dissipation and impact angle relationship

    Fig.11 The relative errors of various analytical methods for impact angle

    It can be seen from Fig.10 that when the impact angle is about 45°,the dissipation energy coefficients obtained by each method are larger,indicating that the energy dissipation of the supply vessel is more and the ship-platform movement response is relatively small.When the impact angle is larger than 45°,the greater the vertical distance of the ship impact direction and the center of gravity of the platform,the greater the motion response to the platform after the collision,and correspondingly,the dissipation coefficient decreases in Fig.10.From the relative errors of various analytical methods for impact angle in Fig.11,and we can also find that the results obtained by Zhang’s method,Liu`2D and Liu`3D are very similar,indicating that the height of the collision is not significant for the roll motion of the platform.On the other hand,we can find that the relative error of Zhang’s method,Liu`2D and Liu`3D is relatively stable,basically between 12%and 15%.Last but not least,the DNV`R method,the maximum error is about 14%when the impact angle is about 45°;when the impact angle is 0°or 90°,DNV`R’s estimation results are very close to the results of NLFEA,and the minimum error is only about 4%.

    5 Conclusions

    Based on DNV specification calculation method,a new analytic method for estimating the dissipated energy of OSV impacting the platform was developed,called DNV`R analytical method.Moreover,we applied Zhang’s method and Liu’s method in the ship collision dissipation energy estimation to the ship-platform collision scene,and worked out programs to realize the dissipation energy of ship platform collision scenarios respectively.Based on the numerical simulation results of NLFEA of ship-platform collision,the results of each analytical method were compared and analyzed.The conclusions are drawn as follows:

    (1)The dissipative energy estimation method proposed by DNV,the specification method,is relatively conservative,and only applicable for central impact.For collision with an offset,the calculation error of the DNV method is more obvious,so the range of application has a larger limitation.Compared with the NLFEA strain energy results,it can be found that the DNV`R method proposed in this paper is suitable for the collision with an offset,and especially suitable for the collision scene of the ocean structure with small response to rolling and pitching motion.

    (2)From the calculation results of the dissipated energy of the striking ship with different impact velocities and different impact angles,it can be found that:the relative errors of various analytical methods are below 15%,except for the DNV dissipation energy estimation method.When estimating the dissipated energy of the striking ship in such collision scenarios,we can use the DNV`R method proposed in this paper or the program based on Zhang’s method and Liu’s method to make a fast estimation.

    (3)In general,the accuracy of the DNV`R method proposed in this paper is close to that of Zhang’s method and Liu’s method,which meets the accuracy requirements of fast estimation of ship-platform collision dissipation energy.As DNV`R’s analytic form is simple and the process of calculating is convenient,it can be used to estimate the dissipated energy of given collision scenarios quickly,so as to provide some reference for the design of the collision structure of marine structures such as semi-submersible drilling platforms.

    [1]Wang Z,Gu Y.Motion lag of struck ship in collision[J].Shipbuilding of China,2001,42(2):56-62.(in Chinese)

    [2]Minorsky V U.An analysis of ship collision to protection of nuclear powered plant[J].J Ship Rsearch,1959(1):1-4.

    [3]Petersen M J.Dynamics of ship collisions[J].Ocean Engineering,1982,9(4):295-329.

    [4]Woisin G.Instantaneous loss of energy with unsymmetrical ship collisions[J].Schiff&Hafen,1988,40(1):50-55.

    [5]Zhang S.The mechanics of ship collisions[D].Denmark:Technical University of Denmark,1999.

    [6]Stronge W J.Impact mechanics[M].Cambridge:Cambridge University Press,2004.

    [7]Liu Z,Amdahl J.A new formulation of the impact mechanics of ship collisions and its application to a ship-iceberg collision[J].Marine Structures,2010,23(3):360-384.

    [8]Liu J,Hu Z.3D analytical method for the external dynamics of ship collisions and investigation of the coefficient of restitution[J].Chinese Journal of Ship Research,2017,12(2):84-91.(in Chinese)

    [9]Mcdermott J,Kline R,Jones E,et al.Tanker structural analysis for minor collisions[J].Deformation,1974.

    [10]Sabril Haris,Jorgen Amdahl.An analytical model to assess a ship side during a collision[J].Ships and Offshore Structures,2012,7(4):431-448.

    [11]Gao Z,Hu Z,Wang G.A simplified analytical method for prediction of anti-collision capability of FPSO side structures[J].Engineering Mechanics,2014,31(s1):155-160.(in Chinese)

    [12]Liu B,Villavicencio R,Soares C G.Simplified analytical method to evaluate tanker side panels during minor collision incidents[J].International Journal of Impact Engineering,2015,78(78):20-33.

    [13]Liu B,Soares C G.Analytical method to determine the crushing behavior of girders with stiffened web[J].International Journal of Impact Engineering,2016,93:49-61.

    [14]Liu K,Wang Z,Zhang Y,et al.Collision behavior of structural analysis in ship collisions based on full-coupling technology[J].Journal of Ship Mechanics,2015,19(5):574-581.(in Chinese)

    [15]Yu Z,Amdahl J,Storheim M.A new approach for coupling external dynamics and internal mechanics in ship collisions[J].Marine Structures,2016,45:110-132.

    [16]DNV-RP-C204.Design against accidental loads[S].2010.

    [17]K?rgesaar M,Romanoff J.Influence of mesh size,stress triaxiality and damage induced softening on ductile fracture of large-scale shell structures[J].Marine Structures,2014,38(5):1-17.

    [18]DNV-RP-208.Determination of structural capacity by non-linear FE analysis methods[S].2013.

    [19]Hogstr?m P,Ringsberg J W.Assessment of the crashworthiness of a selection of innovative ship structures[J].Ocean Engineering,2013,59(59):58-72.

    一種新的碰撞耗散能快速估算方法及其在船—平臺碰撞事故中的應(yīng)用研究

    谷家揚1,鄧炳林1,陶延武1,王 璞2,周 佳2

    (1.江蘇科技大學(xué) 海洋裝備研究院,江蘇 鎮(zhèn)江212003;2.中國船舶及海洋工程設(shè)計研究院,上海 200011)

    海洋平臺在作業(yè)過程中存在遭受船舶碰撞的潛在危險,一旦事故發(fā)生將可能導(dǎo)致結(jié)構(gòu)嚴重損傷從而影響平臺的作業(yè)安全,因此在平臺結(jié)構(gòu)設(shè)計階段需要考慮其抗撞性能。文章以海洋供應(yīng)船撞擊半潛式鉆井平臺為研究對象,基于碰撞外部動力學(xué)理論和DNV規(guī)范計算方法,計及平臺艏搖運動響應(yīng)的影響,提出了一種新的耗散能估算方法-DNV`R解析法。以非線性有限元船-平臺碰撞數(shù)值仿真結(jié)果為參考,將該解析法與規(guī)范法、Zhang解析法及Liu解析法計算結(jié)果進行對比,研究得出:文中所提出的DNV`R解析法彌補了規(guī)范法較為保守且僅適用于對心碰撞場景的不足,計算過程簡單、便捷,且計算精度與Zhang解析法和Liu解析法基本接近,能夠滿足船-平臺碰撞耗散能快速估算的精度要求。

    船—平臺碰撞;耗散能;解析法;碰撞外部動力學(xué)

    U661.1 U674.38+1

    A

    谷家揚(1979-),男,博士,江蘇科技大學(xué)副教授,通訊作者,E-mail:gujayang@126.com;

    鄧炳林(1990-),男,江蘇科技大學(xué)碩士研究生;

    陶延武(1990-),男,碩士,江蘇科技大學(xué)海洋裝備研究院,工程師;

    王 璞(1977-),男,碩士,中國船舶及海洋工程設(shè)計研究院,研究員;

    周 佳(1981-),男,博士,中國船舶及海洋工程設(shè)計研究院,高級工程師。

    U661.1 U674.38+1

    A

    10.3969/j.issn.1007-7294.2017.12.011

    date:2017-09-21

    Supported by Special Project for innovation for‘Seventh generation ultra-deepwater drilling platform(ship)’of high-tech ship of Ministry of Industry and Information Technology;Special Project for innovation for‘Research and development of deepwater semi-submersible support platform’ of hightech ship Ministry of Industry and Information Technology;the National Natural Science Foundation of China(51779109);Natural Science Foundation of Jiangsu,China(BK20171306);the State Key Laboratory of Ocean Engineering of Opening Foundation(1407)

    Biography:GU Jia-yang(1979-),male,Ph.D.,associate professor,corresponding author,E-mail:gujiayang@126.com;DENG Bing-lin(1990-),male,master student.

    1007-7294(2017)12-1564-14

    猜你喜歡
    海洋工程碩士研究院
    昆明理工大學(xué)工商管理碩士(MBA)簡介
    北京食品科學(xué)研究院
    肉類研究(2022年5期)2022-06-16 05:53:24
    工程技術(shù)研究院簡介
    從心所欲不逾矩——為中國戲曲研究院成立70周年作
    戲曲研究(2021年3期)2021-06-05 07:06:46
    不是我!是他搗亂!
    如何寫好碩士博士學(xué)位論文
    海洋工程專家 劉培林
    《海洋工程》第二屆理事會
    海洋工程(2015年1期)2015-10-28 01:36:21
    海洋工程學(xué)會第四屆理事會
    海洋工程(2015年1期)2015-10-28 01:29:14
    北斗RDSS在海洋工程數(shù)據(jù)傳輸中的應(yīng)用
    国产精品国产三级专区第一集| 一区二区日韩欧美中文字幕| 国语对白做爰xxxⅹ性视频网站| www.精华液| 国产一区二区三区av在线| 叶爱在线成人免费视频播放| 亚洲九九香蕉| 精品国产一区二区三区四区第35| 精品少妇内射三级| 麻豆av在线久日| 午夜免费男女啪啪视频观看| 亚洲第一青青草原| 97人妻天天添夜夜摸| av一本久久久久| 国产黄频视频在线观看| 精品熟女少妇八av免费久了| 韩国精品一区二区三区| 啦啦啦视频在线资源免费观看| 亚洲人成电影免费在线| 香蕉国产在线看| 国产黄色免费在线视频| 日韩大码丰满熟妇| 性色av乱码一区二区三区2| 美女高潮到喷水免费观看| 赤兔流量卡办理| 免费黄频网站在线观看国产| 免费在线观看影片大全网站 | 丁香六月欧美| 国产精品一区二区在线不卡| 日本猛色少妇xxxxx猛交久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲熟女精品中文字幕| 日韩大码丰满熟妇| 亚洲欧美精品自产自拍| 国产男女超爽视频在线观看| 老汉色av国产亚洲站长工具| 久久久欧美国产精品| 国产精品久久久av美女十八| 纯流量卡能插随身wifi吗| 午夜两性在线视频| a级毛片黄视频| 最近手机中文字幕大全| 少妇猛男粗大的猛烈进出视频| www.熟女人妻精品国产| 爱豆传媒免费全集在线观看| 99香蕉大伊视频| 欧美xxⅹ黑人| 久久免费观看电影| 国产成人精品无人区| 日韩一本色道免费dvd| 久久久亚洲精品成人影院| 亚洲精品美女久久av网站| 一区福利在线观看| 亚洲国产成人一精品久久久| 97精品久久久久久久久久精品| 在线亚洲精品国产二区图片欧美| 韩国精品一区二区三区| 国产视频一区二区在线看| 欧美性长视频在线观看| 亚洲精品国产av成人精品| 中国国产av一级| 国产一区二区在线观看av| 只有这里有精品99| 精品福利永久在线观看| 夫妻性生交免费视频一级片| 91麻豆精品激情在线观看国产 | 香蕉国产在线看| 日韩av在线免费看完整版不卡| 久久人妻福利社区极品人妻图片 | 中文字幕最新亚洲高清| av又黄又爽大尺度在线免费看| 亚洲国产av新网站| 亚洲一卡2卡3卡4卡5卡精品中文| 超色免费av| 天天躁夜夜躁狠狠久久av| 性色av乱码一区二区三区2| 精品久久久久久久毛片微露脸 | 亚洲国产精品一区二区三区在线| 捣出白浆h1v1| 精品视频人人做人人爽| 久久久久精品国产欧美久久久 | 国产淫语在线视频| tube8黄色片| 国产爽快片一区二区三区| 黄色 视频免费看| 激情视频va一区二区三区| 纯流量卡能插随身wifi吗| 国语对白做爰xxxⅹ性视频网站| av天堂久久9| 18禁裸乳无遮挡动漫免费视频| 欧美日韩成人在线一区二区| 天天添夜夜摸| 国产欧美日韩一区二区三 | 国产精品 欧美亚洲| 美女主播在线视频| 欧美日韩成人在线一区二区| 国产精品国产三级专区第一集| 国产黄频视频在线观看| 日本一区二区免费在线视频| 欧美日韩视频精品一区| 如日韩欧美国产精品一区二区三区| 女警被强在线播放| 亚洲欧洲国产日韩| 18禁黄网站禁片午夜丰满| 亚洲人成电影免费在线| 又大又黄又爽视频免费| 午夜老司机福利片| 99久久精品国产亚洲精品| 少妇人妻久久综合中文| 亚洲国产av新网站| 一本—道久久a久久精品蜜桃钙片| 狂野欧美激情性xxxx| 777米奇影视久久| 夫妻午夜视频| 国产片内射在线| 中国国产av一级| 高潮久久久久久久久久久不卡| 久久青草综合色| 美女福利国产在线| 美国免费a级毛片| 亚洲国产精品成人久久小说| 下体分泌物呈黄色| 97精品久久久久久久久久精品| 亚洲精品第二区| 久久久精品区二区三区| 精品卡一卡二卡四卡免费| 亚洲色图综合在线观看| 男女国产视频网站| 丰满迷人的少妇在线观看| 午夜日韩欧美国产| 精品少妇一区二区三区视频日本电影| 欧美成人午夜精品| 久久久国产欧美日韩av| 久久久欧美国产精品| 啦啦啦 在线观看视频| 男女免费视频国产| 免费女性裸体啪啪无遮挡网站| 精品免费久久久久久久清纯 | 在线观看人妻少妇| 精品人妻一区二区三区麻豆| h视频一区二区三区| 久久女婷五月综合色啪小说| 美女高潮到喷水免费观看| 亚洲精品美女久久av网站| 女人久久www免费人成看片| 国产免费现黄频在线看| 90打野战视频偷拍视频| videosex国产| 久久国产精品大桥未久av| 叶爱在线成人免费视频播放| 满18在线观看网站| 色精品久久人妻99蜜桃| 亚洲精品自拍成人| 一区二区三区精品91| 观看av在线不卡| 亚洲av成人精品一二三区| 七月丁香在线播放| 大码成人一级视频| 欧美少妇被猛烈插入视频| 国产精品 欧美亚洲| 一级片'在线观看视频| 国产精品久久久久久人妻精品电影 | 最近中文字幕2019免费版| 男女之事视频高清在线观看 | 久久久久久亚洲精品国产蜜桃av| 国产麻豆69| 日本午夜av视频| 一级片免费观看大全| 色网站视频免费| 19禁男女啪啪无遮挡网站| 天天添夜夜摸| 亚洲欧美精品自产自拍| 久久久亚洲精品成人影院| 成人手机av| av福利片在线| 成人国语在线视频| 欧美少妇被猛烈插入视频| 黄色a级毛片大全视频| 大片免费播放器 马上看| 高清欧美精品videossex| 亚洲专区中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 一级黄色大片毛片| 亚洲熟女毛片儿| 性欧美人与动物交配| 久久人妻福利社区极品人妻图片| 久热爱精品视频在线9| 老司机深夜福利视频在线观看| 一级黄色大片毛片| 国产黄片美女视频| 黑人操中国人逼视频| 精品少妇一区二区三区视频日本电影| 国产精品亚洲美女久久久| 在线观看66精品国产| 国产一级毛片七仙女欲春2 | 亚洲熟妇熟女久久| 久久精品国产亚洲av高清一级| 可以在线观看的亚洲视频| 国产片内射在线| 男女那种视频在线观看| 国产精品电影一区二区三区| 久久久久久久久中文| 免费在线观看影片大全网站| 免费在线观看完整版高清| 免费看美女性在线毛片视频| 午夜精品在线福利| 久久伊人香网站| 亚洲专区字幕在线| 波多野结衣av一区二区av| 88av欧美| 一夜夜www| 久久九九热精品免费| 91大片在线观看| 欧美日韩瑟瑟在线播放| 女人高潮潮喷娇喘18禁视频| 99精品欧美一区二区三区四区| 两个人视频免费观看高清| 久久久国产成人免费| 免费在线观看亚洲国产| 校园春色视频在线观看| 黑人操中国人逼视频| 亚洲全国av大片| 欧美性猛交黑人性爽| 久久久国产成人免费| 欧美一级a爱片免费观看看 | 成人国产综合亚洲| 美女午夜性视频免费| 亚洲av日韩精品久久久久久密| 亚洲七黄色美女视频| 亚洲精品久久国产高清桃花| 香蕉国产在线看| 亚洲国产中文字幕在线视频| 黄片播放在线免费| 亚洲三区欧美一区| 国产91精品成人一区二区三区| 中国美女看黄片| 在线av久久热| 91在线观看av| 久久精品91蜜桃| 两个人免费观看高清视频| 制服丝袜大香蕉在线| 日本一区二区免费在线视频| 看片在线看免费视频| 777久久人妻少妇嫩草av网站| 久久亚洲精品不卡| 俺也久久电影网| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| 欧美不卡视频在线免费观看 | av在线播放免费不卡| 琪琪午夜伦伦电影理论片6080| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 日韩国内少妇激情av| 999精品在线视频| 国产精品久久久人人做人人爽| or卡值多少钱| av福利片在线| 精品国产国语对白av| 国产野战对白在线观看| 精品国产乱码久久久久久男人| 非洲黑人性xxxx精品又粗又长| 搡老熟女国产l中国老女人| 精品电影一区二区在线| 欧美日韩中文字幕国产精品一区二区三区| 国产真实乱freesex| 伦理电影免费视频| 亚洲欧美精品综合一区二区三区| 97人妻精品一区二区三区麻豆 | www.熟女人妻精品国产| 欧美日韩黄片免| 香蕉av资源在线| 好男人在线观看高清免费视频 | 黄色女人牲交| 久久中文字幕一级| 国产成人欧美在线观看| 亚洲国产欧美网| 亚洲国产毛片av蜜桃av| 91成年电影在线观看| 亚洲欧美精品综合久久99| 精品不卡国产一区二区三区| 亚洲成人国产一区在线观看| 成年人黄色毛片网站| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 韩国av一区二区三区四区| 999精品在线视频| 在线天堂中文资源库| 国产国语露脸激情在线看| 国产片内射在线| 亚洲 国产 在线| 国产黄色小视频在线观看| 日韩欧美在线二视频| 精品卡一卡二卡四卡免费| 久久婷婷成人综合色麻豆| av在线天堂中文字幕| 亚洲精品国产一区二区精华液| 亚洲电影在线观看av| 淫秽高清视频在线观看| 制服诱惑二区| 精品久久久久久久人妻蜜臀av| av天堂在线播放| 国产亚洲精品久久久久5区| 高清在线国产一区| ponron亚洲| av在线天堂中文字幕| 午夜激情福利司机影院| 两性午夜刺激爽爽歪歪视频在线观看 | 精品久久久久久久久久久久久 | 欧美在线黄色| 91av网站免费观看| 亚洲成a人片在线一区二区| 夜夜看夜夜爽夜夜摸| 国产午夜精品久久久久久| 国内精品久久久久久久电影| 精品国产国语对白av| 免费av毛片视频| 好男人在线观看高清免费视频 | 国产精品av久久久久免费| 伦理电影免费视频| 悠悠久久av| АⅤ资源中文在线天堂| 久久久精品国产亚洲av高清涩受| 国产精品亚洲美女久久久| 长腿黑丝高跟| 久久久久国产一级毛片高清牌| 日本熟妇午夜| 日韩欧美一区视频在线观看| 一边摸一边抽搐一进一小说| 热re99久久国产66热| 日韩精品免费视频一区二区三区| 久久九九热精品免费| 热re99久久国产66热| 国产高清激情床上av| 国产精品美女特级片免费视频播放器 | 在线观看www视频免费| 成在线人永久免费视频| 精品久久蜜臀av无| 国产高清视频在线播放一区| 亚洲精品粉嫩美女一区| 最新在线观看一区二区三区| 熟妇人妻久久中文字幕3abv| www日本黄色视频网| 午夜精品在线福利| www日本在线高清视频| 国产91精品成人一区二区三区| 亚洲国产欧美一区二区综合| 两个人看的免费小视频| 亚洲精品一区av在线观看| 免费搜索国产男女视频| 国内精品久久久久精免费| 19禁男女啪啪无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 麻豆一二三区av精品| 亚洲精品久久成人aⅴ小说| 久99久视频精品免费| 欧美日韩瑟瑟在线播放| 亚洲av成人av| 国产成人欧美在线观看| 日本免费一区二区三区高清不卡| 久久久久久久精品吃奶| 怎么达到女性高潮| 男女做爰动态图高潮gif福利片| 午夜精品久久久久久毛片777| 一本综合久久免费| 亚洲av美国av| 大香蕉久久成人网| 日韩精品中文字幕看吧| 亚洲片人在线观看| bbb黄色大片| 777久久人妻少妇嫩草av网站| 99国产综合亚洲精品| 国产高清videossex| 麻豆久久精品国产亚洲av| 黑丝袜美女国产一区| 亚洲国产中文字幕在线视频| 亚洲国产欧美日韩在线播放| 露出奶头的视频| 校园春色视频在线观看| 熟妇人妻久久中文字幕3abv| АⅤ资源中文在线天堂| 中文字幕人妻熟女乱码| 国产成人av激情在线播放| 国产av在哪里看| 国产精品久久久久久人妻精品电影| 国产99久久九九免费精品| 国产精品,欧美在线| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 日日干狠狠操夜夜爽| 国产激情欧美一区二区| 亚洲片人在线观看| 桃色一区二区三区在线观看| 亚洲国产精品合色在线| 国产一卡二卡三卡精品| 99精品久久久久人妻精品| 亚洲三区欧美一区| 免费人成视频x8x8入口观看| 国产亚洲精品综合一区在线观看 | 日本精品一区二区三区蜜桃| 香蕉国产在线看| 日本熟妇午夜| av有码第一页| 久久午夜综合久久蜜桃| 人人妻人人澡欧美一区二区| 日本免费a在线| 在线播放国产精品三级| 久久久久久九九精品二区国产 | 欧美色欧美亚洲另类二区| 欧美一区二区精品小视频在线| www.www免费av| 亚洲午夜精品一区,二区,三区| 视频区欧美日本亚洲| 中出人妻视频一区二区| 国产精品影院久久| 免费人成视频x8x8入口观看| 久久精品亚洲精品国产色婷小说| 成人亚洲精品一区在线观看| 夜夜躁狠狠躁天天躁| 午夜a级毛片| 婷婷丁香在线五月| ponron亚洲| 午夜福利在线观看吧| 99国产精品99久久久久| 两个人免费观看高清视频| av片东京热男人的天堂| 免费在线观看影片大全网站| 黄片大片在线免费观看| av免费在线观看网站| 精品免费久久久久久久清纯| 久久久久久久久免费视频了| 国产精品香港三级国产av潘金莲| 老司机福利观看| 久久精品人妻少妇| 欧美色视频一区免费| 久热爱精品视频在线9| 最近最新免费中文字幕在线| 欧美国产精品va在线观看不卡| 亚洲男人天堂网一区| 热99re8久久精品国产| 亚洲全国av大片| 欧美人与性动交α欧美精品济南到| 精品一区二区三区视频在线观看免费| 国产精品 国内视频| 岛国在线观看网站| 男女之事视频高清在线观看| 男女那种视频在线观看| 中文字幕人妻丝袜一区二区| 久久久国产成人免费| 亚洲欧美一区二区三区黑人| 中文资源天堂在线| 久久久久久国产a免费观看| 啦啦啦观看免费观看视频高清| 国产高清激情床上av| bbb黄色大片| 亚洲人成网站在线播放欧美日韩| 国产午夜精品久久久久久| 亚洲 欧美一区二区三区| 亚洲av美国av| 日本一本二区三区精品| 黄色女人牲交| 校园春色视频在线观看| www.自偷自拍.com| 麻豆成人午夜福利视频| 欧美日韩黄片免| 一本一本综合久久| 亚洲三区欧美一区| а√天堂www在线а√下载| 黄色毛片三级朝国网站| 88av欧美| 两个人看的免费小视频| 精品乱码久久久久久99久播| 观看免费一级毛片| 18禁国产床啪视频网站| √禁漫天堂资源中文www| 男女做爰动态图高潮gif福利片| 亚洲成a人片在线一区二区| 国产亚洲精品第一综合不卡| 好男人在线观看高清免费视频 | 欧美色视频一区免费| 国产午夜福利久久久久久| 成年人黄色毛片网站| 妹子高潮喷水视频| 欧美成人一区二区免费高清观看 | 俺也久久电影网| 一区二区三区激情视频| 亚洲精品久久国产高清桃花| 午夜两性在线视频| 亚洲五月天丁香| 亚洲 欧美一区二区三区| 国产一卡二卡三卡精品| x7x7x7水蜜桃| 国产精品久久久久久亚洲av鲁大| 18禁黄网站禁片午夜丰满| 精品国产超薄肉色丝袜足j| 老司机靠b影院| 啪啪无遮挡十八禁网站| 国产成+人综合+亚洲专区| 亚洲国产欧洲综合997久久, | 一进一出好大好爽视频| 国产精品永久免费网站| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 国产精品久久久人人做人人爽| www.熟女人妻精品国产| 亚洲片人在线观看| 久久精品国产99精品国产亚洲性色| 18禁黄网站禁片免费观看直播| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 国产av又大| 不卡一级毛片| 国产国语露脸激情在线看| 午夜福利在线在线| 久久草成人影院| 日韩 欧美 亚洲 中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲第一av免费看| 精品午夜福利视频在线观看一区| 成人三级黄色视频| 国产单亲对白刺激| 国产野战对白在线观看| 亚洲一区高清亚洲精品| av中文乱码字幕在线| 美女 人体艺术 gogo| 午夜福利在线在线| 18禁黄网站禁片午夜丰满| av天堂在线播放| 午夜久久久在线观看| 亚洲国产精品999在线| 最新在线观看一区二区三区| 婷婷精品国产亚洲av在线| 狂野欧美激情性xxxx| 一区福利在线观看| 日本在线视频免费播放| 亚洲精品久久成人aⅴ小说| 欧美另类亚洲清纯唯美| 午夜a级毛片| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久av网站| 亚洲在线自拍视频| 中文字幕人妻丝袜一区二区| 50天的宝宝边吃奶边哭怎么回事| 亚洲av日韩精品久久久久久密| 欧美 亚洲 国产 日韩一| 俄罗斯特黄特色一大片| 精品高清国产在线一区| 美女高潮到喷水免费观看| 久久精品亚洲精品国产色婷小说| 亚洲国产欧美一区二区综合| 精品少妇一区二区三区视频日本电影| 十八禁人妻一区二区| 超碰成人久久| 欧美zozozo另类| 亚洲精品在线美女| 亚洲人成电影免费在线| 欧美亚洲日本最大视频资源| 1024手机看黄色片| 亚洲全国av大片| 91麻豆av在线| 91麻豆精品激情在线观看国产| 中文在线观看免费www的网站 | 国产欧美日韩一区二区三| 国产激情偷乱视频一区二区| 黄片小视频在线播放| 国产aⅴ精品一区二区三区波| 色综合亚洲欧美另类图片| 久久久久精品国产欧美久久久| 日韩三级视频一区二区三区| 看免费av毛片| 19禁男女啪啪无遮挡网站| 国产激情偷乱视频一区二区| 国产一区二区三区在线臀色熟女| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美av亚洲av综合av国产av| 国产高清videossex| 黄色视频不卡| 欧美乱色亚洲激情| 日本 欧美在线| 国产精品综合久久久久久久免费| 夜夜夜夜夜久久久久| 亚洲国产看品久久| 黄色丝袜av网址大全| 亚洲精品美女久久久久99蜜臀| 欧美乱色亚洲激情| 久久精品成人免费网站| 777久久人妻少妇嫩草av网站| 99久久无色码亚洲精品果冻| 国产午夜福利久久久久久| 国产欧美日韩一区二区精品| 大型黄色视频在线免费观看| xxxwww97欧美| 国产主播在线观看一区二区| 我的亚洲天堂| 久久中文字幕一级| 亚洲五月色婷婷综合| 99在线视频只有这里精品首页| 亚洲真实伦在线观看| 91麻豆av在线| 亚洲av成人av| 怎么达到女性高潮| 一二三四社区在线视频社区8| 国产熟女xx| 亚洲精华国产精华精| 日本免费一区二区三区高清不卡| 人妻久久中文字幕网| 日本精品一区二区三区蜜桃| 又黄又粗又硬又大视频| 一本久久中文字幕| 成人国产一区最新在线观看| 校园春色视频在线观看|