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

    Study on the Connector Loads Characteristics of a Very Large Floating Structure in the Impact Accident

    2016-05-16 02:41:54GUJiayangZHUYueXIEYulinGUANGYifengQIEnrongLUYe
    船舶力學(xué) 2016年12期
    關(guān)鍵詞:浮體碩士江蘇

    GU Jia-yang,ZHU Yue,XIE Yu-lin,GUANG Yi-feng,QI En-rong,LU Ye

    (1.School of Naval Architecture and Marine Engineering,Jiangsu University of Science and Technology, Zhenjiang 212003,China;2.China Ship Scientific Research Center,Wuxi 214082,China)

    Study on the Connector Loads Characteristics of a Very Large Floating Structure in the Impact Accident

    GU Jia-yang1,ZHU Yue1,XIE Yu-lin1,GUANG Yi-feng1,QI En-rong2,LU Ye2

    (1.School of Naval Architecture and Marine Engineering,Jiangsu University of Science and Technology, Zhenjiang 212003,China;2.China Ship Scientific Research Center,Wuxi 214082,China)

    Accurate assessment of the connector loads characteristics has an important engineering significance for the hydrodynamic design and structural safety of the very large floating structure (VLFS).In this paper,the combined effects of ship impact,wind,wave and current load are considered to study a semi-submersible VLFS consisting of five modules.The AQWA software is used to carry out the numerical calculation for loads characteristics of connectors.The ship impact load is simplified as a rectangular pulse,and the steady wind,steady flow and JONSWAP spectrum are used to discuss the influence of the ship impact angle,the impact position,the environmental load and other factors on the connector loads characteristics.The studied results can provide some basis for the design and strength calculation of the VLFS.

    VLFS;connector;impact load;environmental load;motion response

    0 Introduction

    In recent years,with the rapid development of the exploiting and utilization of marine resources,the design schemes for different types of VLFS have been proposed.Very large floating structures(VLFSs)are generally connected by flexible connectors that allow the relative motion of modules.In the harsh marine environment,VLFSs suffer the combined effects of wind, wave and current load,which cause large relative deformation between modules,and the stress caused by the deformation is borne by the connector.At present,there is no mature product can be used as the connectors between the very large floating modules.The connectors remain in the research phase,so it has an important engineering significance to carry out the dynamic load forecasting and structure design of connectors.

    Wang et al[1]adopted a time sequence analysis method to study the extreme motion re-sponse and connector extreme load of a VLFS which consists of three semi-submersible modules by considering the effect of connector stiffness and using the Airy wave theory and Morrison equation.Qi et al[2]designed a connector model which is composed of a ball hinge and a spring sliding block device,and conducted experiments to study the connector stiffness characteristic of shallow draft VLFS,as well as floating body motion response under dynamic load. Based on the three dimensional potential theory,Liu et al[3]adopted rigid module flexible connector(RMFC)model to carry out a short term forecast for the hydrodynamic performance and connector load of VLFS consisting of 5-modules.This paper also compared the response rules of connectors under different size of floating body structure,connector stiffness and two typical draft(300 m and 50 m).By an extended hydroelasticity theory taking into account hinge rigid modes,Fu et al[4]studied the effects of connector and module stiffness on the hydroelastic response of a two-module interconnected structure with flexible modules and connectors.The characteristics of the bending moment distribution are also discussed.Based on three dimensional hydroelastic theory,Fu et al[5]carried out numerical calculations for a VLFS model in regular waves with varied bending stiffness and wave frequency.The effect of bending stiffness on hydroelasticity response of VLFS was discussed in detail.The numerical results of hydroelastic behaviors included vertical response amplitude,generalized displacement and bending moment.Lai et al[6]used the numerical simulation method to calculate the motion response of VLFS consisting of five semi-submersible modules connected by hinge type connector.Based on hydroelasticity theory considering the second-order fluid loads induced by the coupling of first-order wave potentials,Chen et al[7]calculated the bending moments and the vertical displacements of a floating plate in multidirectional irregular waves.Gu et al[8]proposed a time domain method to calculate the connector load of the VLFS in waves,and adopted RMFC model to study the motion characteristics of connector load under impact load.Considering a ribbon bridge subjected to a moving load,Wang et al[9]analyzed the hydroelastic response of the floating bridge and the fatigue behavior of the connectors by employing the local stress-strain approach.The result showed that fatigue life of the connector sharply goes down with the increase of the passing speed of the moving load.Basing on the potential flow theory,Yu et al[10-11]analyzed the effect of sea conditions,connector stiffness and interaction of modules on the connector load by self-fabricated software HYCOFELS.Miguel et al[12]introduced the development history,functional characteristics,structural form and connector form of VLFS,Mega-Float,MOB(Mobile offshore Base),and PSP(Pneumatically Stabilized Platform).Riggs et al[13]studied the effect of connector stiffness on the motion response of VLFS consisting of five modules by adopting RMFC,and discussed the connector load and the motion characteristics of VLFS under different wave frequency and wave direction.Riggs et al[14]simplified MOB as RMFC model and investigated the influence of the connection scheme and the damping of the connector on the motion of VLFS.Riggs et al[15]used two-dimensional and three-dimensional hydroelasticity methods to predict the responses of a 5-module and 16-module VLFS in regular waves.Then,the effects of the interaction between the modules on the connector load and the motion of module were analyzed.Riggs et al[16]studied the motion characteristics of twotypical structures of 5-modules MOB.

    The dynamic motion characteristics of the VLFS subjected to impact load are studied by other related researchers,too.Watanabe et al[17]used the finite element program to analyze the transient response of VLFS under impact loading caused by plane landing.Hisayoshi et al[18]studied the motion response of the VLFS under the combined action of plane takeoff,landing and wave motion,and calculated the drag loads acting on the aircraft caused by the deformation of the VLFS.Kashiwagiet al[19]calculated the elastic transient deformation of a pontoontype VLFS during plane takeoff and landing.

    1 VLFS and mooring system

    This paper studies a VLFS consisting of five-semi-submersible modules with a flexible hinged connector linking the head and tail of modules.The single module is shown in Fig.1. The length is 300 m,the width is 100 m,the depth is 27 m,and there are 5 floats,the float is assigned laterally along the width direction of platform,the longitudinal spacing of the floats center is 30 m.There are two columns 16 m high and 18 m in diameter between the float and the deck,the transverse spacing of the column center is 60 m.There are two cross-brace 30 m long and 3 m in diameter designed to support structure between the floats,and the transverse distance between the central axis of cross-braces is 60 m.

    Fig.1 Single semi-submersible model

    Fig.2 Schematic diagram of VLFS mooring scheme

    The VLFS is assumed to work in the 40 m uniform water depth of the lagoon by using catenary mooring system.Schematic diagram of VLFS mooring scheme is shown in Fig.2,the length of the mooring cable is 510 m,and the mooring radius is 500 m.The upper part of the mooring cable is connected with the floating body,and the length of unstudded cable is 60 m. The diameter cable is 146 mm,breaking tension is 1.405×107N,the initial tension of the mooring system is provided by the mooring line.The weight and stiffness of the mooring system are increased by using the form of the three chains in the middle segment.The diameter of chain is 252 mm,the length of chain is 170 m,and the breaking tension is 3.088×107N.The bottom section connected with the seabed anchor is adopted by steel cable,the diameter of cable is 150 mm,the breaking tension of cable is 1.452×107N,the length of cable is 280 m,which provide enough lying area to prevent the up-pull of anchor.

    2 Ship impact load,ocean environment and calculation condition

    2.1 Impact condition and transient impact load

    VLFSs are mainly used in the coastal areas.So the impact accident may occur when theships pass or berthed.According to the arrangement of 5-modules floating body and connectors,the position of the ship impact is assumed to be located in the column.As both ends of 3 modules in the middle are connected with connectors,the modules in the left and right sides only have one end subject to the role of the connector which suffer more impact of ship impact.At the same time,due to the symmetry of the structure,there is only one floating body in the left or right sides is considered in the impact condition.The top view of left-side single module is shown in Fig.3.The modules are labeled as M1-M5.The arrow indicates the ship impact angle,and the height position of the impact position is assumed at the waterline.

    Fig.3 Schematic diagram of the connection type

    Fig.4 The impact load time history curve

    Ship impact load is selected by referencing impact load formula of Norway bridge load code and modified Worthing formula.By assuming that a bulk carrier with a load of 20 000 T at speed 12 kns hits column,the calculation result with Norway bridge load code is 98.2 MN, and the calculation result with modified Worthing formula is 104.7 MN.In this paper,the ship impact load is selected as 100 MN.The load model is simplified as a rectangular load with a duration of 3 s.The impact load time history curve is shown in Fig.4.

    2.2 Environmental load

    The environmental load of VLFS is mainly influenced by wind,current and wave.In order to facilitate the study,the velocity of constant wind is chosen as 20 m/s,the velocity of uniform inflow is selected as 0.5 m/s.Wind load and flow load are calculated by NUMECA Fine Marine software.When the flow direction angles are 0°,45°or 90°,the load of wind and flow affected on a single floating body module is shown in Tab.1.The wave is chosen as JONSWAP spectrum,the significant wave height is 4.5 m,spectra peak period is 9.63 s.

    2.3 Calculation condition

    In this paper,the combined effects of ship impact,wind,wave and current load are considered to study the load characteristics of VLFS.The wind,wave and current are assumed in the same direction as ship impact to obtain the maximum value of the connector load.M1-M5 columns are impacted by ship respectively.The ship impact angles of 0°,45°,90°and 135°are considered to act on M1,and the ship impact angles of 45°,90°and 135°are considered toact on M2-M5.There are 16 kinds of working condition.Hydrodynamic software AQWA is used to simulate various impact conditions and environmental loads.The total computing time is 2 000 s, time step is 0.2 s.The Schematic diagram of hydrodynamic model of VLFS in AQWA is shown in Fig.5.The four connectors are labeled as C1-C4,and the global coordinate system of the whole model and the local coordinate system of the connectors are given,too.

    Tab.1 The load of wind and flow in a single floating body module

    Fig.5 The Schematic diagram of hydrodynamic model of VLFS in AQWA

    3 Calculation results and analysis

    3.1 The load of connection point

    The statistic value of load on the connector in the three directions when the ship hits the column M1 by four typical angles are given in Fig.6.As shown in figure,when the ship impact angle is 90°,the connector C1 which is near the impact point suffer the largest load in the Y direction.The value of the largest load is 3.49×107N.The maximum loads in the X direction and Z direction appear at the connector C1 with the ship impact angle 0°,which are 7.09×107N and 5.13×107N respectively.As the difference of the ship impact angles,there is a difference in the directions of the maximum stress on the connector.When the ship impact angle is 90°,the load in Y direction reach the maximum.Similarly,the load in X direction is large at the 0°ship impact angle.When the ship impact angle is close to 0°or 180°,the load in X direction becomes larger.The maximum load in Z direction also appears at connector C1 under ship impact angle 0°,which is far greater than the loads of other connectors and the loads under other ship impact angle.

    Fig.6 The statistic value of load on the connector in the three directions

    The loads statistics of each connector when different columns are hit are shown in Fig.7. As can be seen from Fig.7(b),when the ship hits the column M5 at 90°,the maximum loads of all connectors in the Y direction appear and the largest loads is 4.98×107N.As shown in the Fig.7(a)and Fig.7(c),the load of each connector in X direction and Z direction dose not change with the different impact columns.The reason is that the effect of ship impact on load of each connector in Z or X direction is much less than that of the environment.The loads in X direction and Z direction are mainly produced by environmental load.

    Fig.7 The loads statistics of each connector under different working conditions

    The load time history curve of each connector in X direction when the ship impacts the column M1 with angle 0°is shown in Fig.8(a).The load time history curve of each connectorin Y direction when the ship impacts the column M5 with angle 90°is shown in Fig.8(b).From the Fig.8,it is obvious that the Fxand Fyof each connector have a clear oscillation in the 1 000 s,the maximum load oscillation amplitude occurs at C1.As shown in Fig.8(c),the time history curve of Fzat each connector has no obvious change,indicating that the effect of ship impact on load of each connector in Z direction is small.

    Fig.8 The load time history curve of each connector in three directions

    3.2 The bending moment characteristic of hinge point

    Each floating body module is connected through a flexible hinged connector,and each connector can be rotated around the Y axis.In this paper,the calculated value of Myis very small,and the maximum amplitude is only 138N·m,So it is unnecessary to consider the influence of bending moment My.From Fig.9(a),it is found that the amplitude of Mxis larger when the ship impact column M1 with angle 90o,and the maximum amplitude appears at the connector C2 which is 1.67×109N·m.The bending moments Mxof two connectors in the middle are larger than that of two connectors at sides.Fig.9(b)shows the statistics value of Mzof each connector when the column M1 is impacted.When the ship impact angle is 90°,the maximum amplitude appears at the connector C2,which is 2.03×1010N·m.

    As shown in Fig.10(a),the impact point is farther away from the connector,the value of the bending moment Mxin each connector is larger.The distance between the impact points is greater,the lever arm is larger,the value of bending moment is larger,which is consistent with common sense.The maximum amplitude of bending moment Mxappears at the connector C2when the ship impact column M1,which is 1.67×109N·m.From Fig.10(b),it is found that the maximum amplitude Mzoccurs at connector C2 when the column M1 is impacted,and the value is 2.03×1010N·m.

    Fig.10 The statistics of bending moments when ship impact angle is 90°

    As shown in Fig.11,the maximum value of Mxappears at 1 015 s.Before the impact,the fluctuation of bending moment at the connector caused by the environmental load is small. During the impact,the bending moment has an obvious oscillation,but the oscillation of Mxis weaker than that of Mz.After 1 040 s,the oscillation of Mxbegins to decay.In the 1 000 s of ship impact time,the bending moment Mzhas a larger oscillation,the duration time is 3 s.The bending moment Mzat connector C2 has the greatest oscillation,but the bending moment Mzat connector C4 has the smallest oscillation.

    Fig.11 The time history curve of bending moments when the ship hits the column M1 by 90°

    4 The influence of environmental load on the load of the connectors

    When the ship impact different columns with different angles,the time history curve of load at each connector is shown in Fig.12.The influence of environmental load on the amplitude of Fxat each connector is negligible.The load time history curve of each connector without environmental loads is almost coincide with that under environmental loads.The change trend of time history curve of Fyis influenced by environmental load.The environmental loadhas a certain contribution to the amplitude of Fz,which causes an obvious oscillation on the time history curve of Fz.

    Fig.12 The time history curve of load of each connector with or without environmental load

    Fig.13 The time history curves of Mxand Mzat connector C2

    Due to the hinge type is adopted for the connectors,the value of Mycan be ignored.When the ship impacts the column M1 at 90°,the time history curves of Mxand Mzat connector C2 are shown in Fig.13.It is found that the oscillation of Mxwith environmental load is more severe than the oscillation of Mxwithout environmental load.But the environmental load has small effect on Mz.

    5 Conclusions

    In this paper,the AQWA software was used to carry out the numerical calculation for the connector loads characteristics of VLFS under ship impact.The ship impact angle,the impact position and the environmental load are considered to calculate and analyze the load of each connector.The conclusions are summarized as follows:

    (1)The bending moment along the vertical direction at each connector is the greatest;the ship impact angle and the impact position have a great influence on the connector.The largest bending moment appears when the impact point is far from the connector and the ship impact angle is 90°.

    (2)The impact position has a great influence on the load of each connector.The influence of load on connectors in the middle is greater than connectors at sides.

    (3)The bending moments along the vertical direction and the transverse loads of different connectors are less affected by the environmental load.But the environmental load has a greater impact on vertical load,longitudinal load and bending moment along the longitudinal direction.

    [1]Wang Pu,Yu Lan,Li Runpei,A time sequence analysis method for predicting connector loads in a semi-submersible Very Large Floating Structure[J].Ocean Engineering,2002,20(3):9-13.

    [2]Qi Enrong,Liu Chao,Xia Jinsong,Lu Ye,Li Zhiwei,Yue Yalin.Experimental study of functional simulation for flexible connectors of very large floating structure[J].Journal of Ship Mechanics,2015,19(10):1245-1254.

    [3]Liu Chao,Qi Enrong,Lu Ye.Dynamic response of connectors of Very Large Floating Structures under shallow draft[J]. Journal of Ship Mechanics,2015,19(10):1245-1254.

    [4]Wei Wei,Fu Shixiao,Fei Guo,Liang Yuanhua.The effect of bending stiffness on hydroelastic response ofV LFS[C]//Proceedings of the ASME 33rd International Conference on Ocean,Offshore and Arctic Engineering,San Francisco,June 8-13,2014.California,USA,2014-23710.

    [5]Fu Shixiao,Moan Torgeir,Chen Xujun,Cui Weicheng.Hydroelastic analysis of flexible floating interconnected structures [J].Ocean Engineering,2007,34:1516-1531.

    [6]Lai Zhimeng,Cheng Xiaoming,Tian Chao,Zhang Kai.Effect of the connector stiffness on the responses of multi-moduled Very Large Floating Structures in waves[C].Proceedings of 13th National Conference on water dynamics and the 26th National Symposium on water dynamics,2014:1484-1494.

    [7]Chen Xujun,Moan Torgeir,Fu Shixiao.Extreme response of very large floating structure considering second-order hydroelastic effects in multidirectional irregular waves[J].Journal of Offshore Mechanics and Arctic Engineering,2010,132(4): 1-11.

    [8]Gu Jiayang,Wu Jie,Qi Enrong,Guan Yifeng,Yuan Yu-bo.Time domain calculation of connector loads of a Very Large Floating Structure[J].Journal of Marine Science and Application,2015,14(2):183-188.

    [9]Wang Cong,Fu Shixiao,Cui Weicheng.Hydroelasticity based fatigue assessment of the connector for a ribbon bridge subjected to a moving load[J].Marine Structures,2008,22(2):246-260.

    [10]Yu Lan,Li Runpei,Shu Zhi.Dynamic characteristics of mobile offshore base connectors[J].Journal of Shanghai Jiaotong University,2003,37(8):1159-1163.

    [11]Yu Lan,Ding Wei,Li Runpei.Effect of the multiple modules interaction on MOB connector loads[J].Ocean Engineering, 2004,22(1):25-31.

    [12]Miguel L P,Gregorio I,Luis C.A review of Very Large Floating Structures(VLFS)for coastal and offshore uses[J].Ocean Engineering,2015,109:677-690.

    [13]Riggs H R,Ertekin R C,Mills T R J.Impact of connector stiffness on the response of a multi-module mobile offshore base[C]//Proceedings of the International Offshore and Polar Engineering Conference,May 24-29,1998.Montreal,Can, 1998:200-207.

    [14]Riggs H R,Ertekin R C.Response characteristics of serially connected semisubmersibles[J].Journal of Ship Research,1999,43(4):229-240.

    [15]Riggs H R,Ertekin R C.Approximate methods for dynamic response of multi-module floating structures[J].Marine Structure,1993,6(2-3):117-141.

    [16]Riggs H R,Ertekin R C,Mills T R J.Comparative study of RMFC and FEA models for the wave-induced response of a MOB[J].Marine Structure,2000,13(4-5):217-232.

    [17]Watanabe Eiichi,Utsunomiya Tomoaki,Tanigaki Shinkichi.Transient response analysis of a Very Large Floating Structure by Finite Element Method[J].Structural Engineering/Earthquake Engineering,1998,15(2):155-163.

    [18]Hisayoshi Endo.The behavior of a VLFS and an airplane during takeoff/landing run in wave condition[J].Marine Structures,2000,13(4-5):477-491.

    [19]Kashiwagi M,Higashimachi N.Numerical simulations of transient responses of VLFS during landing and take-off of an airplane[C]//Proceedings of the International symposium of Ocean Space Utilization Technology,January 28-31,2003. Tokyo,Japan,2003:83-97.

    碰撞事故時(shí)超大型浮體模塊連接處荷載特性研究

    谷家揚(yáng)1,朱玥1,謝玉林1,管義鋒1,祁恩榮2,陸曄2

    (1.江蘇科技大學(xué)船舶與海洋工程學(xué)院,江蘇鎮(zhèn)江212003;2.中國(guó)船舶科學(xué)研究中心,江蘇無(wú)錫214082)

    精確評(píng)估模塊連接處荷載特性對(duì)超大型浮體的水動(dòng)力設(shè)計(jì)及結(jié)構(gòu)安全性研究具有重要的工程意義。文章以5模塊半潛式超大型浮體為研究對(duì)象,考慮碰撞荷載和風(fēng)、浪、流環(huán)境荷載的聯(lián)合作用并采用AQWA軟件對(duì)超大型浮體模塊連接點(diǎn)處的受力特性開(kāi)展了數(shù)值計(jì)算。船舶碰撞荷載簡(jiǎn)化為矩形脈沖,采用定常風(fēng)、定常流以及JONSWAP波譜,探討了船舶碰撞角度、撞擊位置、環(huán)境荷載等因素對(duì)模塊連接點(diǎn)受力特性的影響。該文的研究成果可為事故荷載作用時(shí)超大型浮體連接器結(jié)構(gòu)設(shè)計(jì)和強(qiáng)度計(jì)算提供一定依據(jù)。

    超大型浮體;連接器;碰撞荷載;環(huán)境荷載;運(yùn)動(dòng)響應(yīng)

    U661.43

    A

    谷家揚(yáng)(1979-),男,博士,江蘇科技大學(xué)副教授,通訊作者,E-mail:gujayang@126.com;朱玥(1991-),男,江蘇科技大學(xué)碩士研究生;謝玉林(1994-),男,江蘇科技大學(xué)碩士研究生;管義鋒(1966-),男,江蘇科技大學(xué)教授;祁恩榮(1965-),男,博士,中國(guó)船舶科學(xué)研究中心研究員;陸曄(1983-),男,碩士,中國(guó)船舶科學(xué)研究中心工程師。

    U661.43 < class="emphasis_bold">Document code:A

    A

    10.3969/j.issn.1007-7294.2016.12.007

    1007-7294(2016)12-1573-11

    Received date:2016-08-29

    Foundation item:Supported by National Key Basic Research and Development Program(973)(No.2013CB036104);the National Natural Science Foundation of China(51309123);the Open Foundation of State Key Laboratory of Ocean Engineering(1407)and‘Qing Lan Project’of Colleges and Universities in Jiangsu Province, Academic Program Development of Jiangsu Higher Education Institutions(PAPD),the Collaborative Innovation Center Funded Projects in Jiangsu University(High Technology Ship Category)

    Biography:GU Jia-yang(1979-),male,associate professor,corresponding author,E-mail:gujiayang@126.com; QI En-rong(1965-),male,Ph.D.,researcher.

    猜你喜歡
    浮體碩士江蘇
    浮體結(jié)構(gòu)沉浮過(guò)程周圍水流特性研究
    昆明理工大學(xué)工商管理碩士(MBA)簡(jiǎn)介
    數(shù)讀江蘇
    物探船硬浮體陣列自擴(kuò)變量分析與應(yīng)用
    超大型浮體結(jié)構(gòu)碰撞損傷研究
    數(shù)獨(dú)江蘇
    數(shù)讀江蘇
    江蘇
    有限流動(dòng)水域浮體受力及側(cè)傾研究
    如何寫好碩士博士學(xué)位論文
    亚洲图色成人| 欧美bdsm另类| 国产精品 国内视频| 国产精品久久久久久久电影| 免费高清在线观看日韩| 人体艺术视频欧美日本| 欧美xxⅹ黑人| 美女脱内裤让男人舔精品视频| 99国产精品免费福利视频| 人妻 亚洲 视频| 国产日韩欧美亚洲二区| 亚洲精品中文字幕在线视频| 熟女人妻精品中文字幕| 老司机影院成人| 精品久久久精品久久久| 国产免费视频播放在线视频| 免费观看无遮挡的男女| 免费黄色在线免费观看| 日韩一区二区视频免费看| 亚洲成色77777| 国产精品成人在线| 亚洲精品国产色婷婷电影| 毛片一级片免费看久久久久| 五月天丁香电影| 韩国av在线不卡| 蜜桃国产av成人99| 精品99又大又爽又粗少妇毛片| 97精品久久久久久久久久精品| 久久精品国产综合久久久 | 一二三四在线观看免费中文在 | 国产精品麻豆人妻色哟哟久久| av黄色大香蕉| 18禁动态无遮挡网站| 黄色 视频免费看| 久久精品夜色国产| 秋霞在线观看毛片| 中文字幕制服av| 18禁国产床啪视频网站| 欧美日韩成人在线一区二区| 国产高清国产精品国产三级| 99热6这里只有精品| 亚洲欧美色中文字幕在线| 日韩熟女老妇一区二区性免费视频| 在线观看免费日韩欧美大片| 观看美女的网站| 97人妻天天添夜夜摸| av在线app专区| 日韩,欧美,国产一区二区三区| 亚洲一区二区三区欧美精品| 午夜久久久在线观看| 日韩三级伦理在线观看| 性色avwww在线观看| 亚洲国产成人一精品久久久| 在线精品无人区一区二区三| 波多野结衣一区麻豆| 欧美另类一区| 麻豆乱淫一区二区| 精品国产一区二区三区四区第35| 国产成人精品无人区| 久久久久久伊人网av| 久久97久久精品| 国产精品一二三区在线看| 在线观看免费视频网站a站| 国产在线免费精品| 草草在线视频免费看| 中国国产av一级| 最近的中文字幕免费完整| 日本欧美视频一区| 欧美精品人与动牲交sv欧美| 中文精品一卡2卡3卡4更新| 一本大道久久a久久精品| 内地一区二区视频在线| 国产在线视频一区二区| 少妇人妻 视频| 日本av手机在线免费观看| 国产乱来视频区| 午夜免费观看性视频| 国产成人精品久久久久久| 色婷婷av一区二区三区视频| av国产久精品久网站免费入址| 免费黄网站久久成人精品| 国产在线视频一区二区| 亚洲国产av新网站| 2021少妇久久久久久久久久久| 久久精品国产鲁丝片午夜精品| 欧美激情极品国产一区二区三区 | 国产在线一区二区三区精| 久久久久精品人妻al黑| 亚洲av成人精品一二三区| 久久韩国三级中文字幕| 成年av动漫网址| 亚洲国产日韩一区二区| 久久这里只有精品19| 免费观看无遮挡的男女| 中文字幕制服av| 成年人免费黄色播放视频| 欧美日韩视频精品一区| 99久久综合免费| 精品99又大又爽又粗少妇毛片| 免费大片黄手机在线观看| 亚洲色图 男人天堂 中文字幕 | 亚洲国产看品久久| 人成视频在线观看免费观看| av天堂久久9| 97精品久久久久久久久久精品| av视频免费观看在线观看| 91午夜精品亚洲一区二区三区| 免费av中文字幕在线| 我要看黄色一级片免费的| 涩涩av久久男人的天堂| 22中文网久久字幕| 久久午夜综合久久蜜桃| 欧美成人精品欧美一级黄| 深夜精品福利| 精品一区二区免费观看| 色婷婷久久久亚洲欧美| 18+在线观看网站| 亚洲一区二区三区欧美精品| 啦啦啦中文免费视频观看日本| 999精品在线视频| 欧美+日韩+精品| 欧美人与性动交α欧美精品济南到 | 香蕉国产在线看| 在线观看美女被高潮喷水网站| 妹子高潮喷水视频| 日本黄色日本黄色录像| 亚洲欧美成人综合另类久久久| 亚洲高清免费不卡视频| 国产色婷婷99| 国产黄色免费在线视频| 少妇猛男粗大的猛烈进出视频| 男人舔女人的私密视频| 日本-黄色视频高清免费观看| 亚洲综合色惰| 色婷婷久久久亚洲欧美| 亚洲,欧美精品.| freevideosex欧美| 日本爱情动作片www.在线观看| 晚上一个人看的免费电影| 国产一区二区三区综合在线观看 | 丁香六月天网| 亚洲一码二码三码区别大吗| 亚洲欧洲日产国产| 超碰97精品在线观看| 欧美+日韩+精品| 一区在线观看完整版| 91在线精品国自产拍蜜月| 婷婷成人精品国产| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 亚洲精品成人av观看孕妇| 久久亚洲国产成人精品v| 免费人妻精品一区二区三区视频| 亚洲欧洲日产国产| 一二三四在线观看免费中文在 | 在线观看三级黄色| 免费看av在线观看网站| 国产av精品麻豆| 日韩中文字幕视频在线看片| 久久精品夜色国产| 欧美日本中文国产一区发布| 18禁裸乳无遮挡动漫免费视频| 在线 av 中文字幕| 亚洲性久久影院| 少妇被粗大猛烈的视频| a 毛片基地| 深夜精品福利| 国产一区亚洲一区在线观看| 久久久精品区二区三区| 你懂的网址亚洲精品在线观看| 日韩三级伦理在线观看| 蜜桃国产av成人99| 国产成人精品无人区| 欧美精品av麻豆av| 亚洲伊人色综图| 大话2 男鬼变身卡| 久久97久久精品| 91成人精品电影| 天天躁夜夜躁狠狠久久av| 欧美性感艳星| 精品亚洲乱码少妇综合久久| 日韩一本色道免费dvd| 精品第一国产精品| 免费黄网站久久成人精品| 90打野战视频偷拍视频| 免费久久久久久久精品成人欧美视频 | 免费少妇av软件| 日韩av不卡免费在线播放| 国产成人av激情在线播放| 丝瓜视频免费看黄片| 最近最新中文字幕大全免费视频 | 国产午夜精品一二区理论片| 亚洲美女视频黄频| 亚洲欧美一区二区三区国产| 女人久久www免费人成看片| 久久久久久伊人网av| 亚洲欧美一区二区三区国产| 女性生殖器流出的白浆| 国产男人的电影天堂91| 麻豆精品久久久久久蜜桃| 亚洲久久久国产精品| 又大又黄又爽视频免费| 麻豆乱淫一区二区| 大陆偷拍与自拍| 免费久久久久久久精品成人欧美视频 | 欧美激情 高清一区二区三区| 亚洲av成人精品一二三区| 成年人午夜在线观看视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区av电影网| a级毛色黄片| 国产白丝娇喘喷水9色精品| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影| 亚洲精品乱码久久久久久按摩| av国产精品久久久久影院| 日本91视频免费播放| 亚洲一区二区三区欧美精品| 高清欧美精品videossex| 亚洲,欧美精品.| 久久久久久伊人网av| 亚洲五月色婷婷综合| 亚洲人成77777在线视频| av又黄又爽大尺度在线免费看| 一边摸一边做爽爽视频免费| 国产白丝娇喘喷水9色精品| 男女边吃奶边做爰视频| 十八禁高潮呻吟视频| 日韩av在线免费看完整版不卡| 在线观看美女被高潮喷水网站| 午夜福利视频在线观看免费| 日本午夜av视频| 色婷婷久久久亚洲欧美| 亚洲人成77777在线视频| 国产片内射在线| 九九爱精品视频在线观看| av网站免费在线观看视频| 一二三四中文在线观看免费高清| 久久婷婷青草| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的| 黑丝袜美女国产一区| 精品一区在线观看国产| 亚洲欧美色中文字幕在线| 制服诱惑二区| 少妇熟女欧美另类| 在线天堂最新版资源| 成人漫画全彩无遮挡| 91成人精品电影| 欧美精品一区二区免费开放| 国产精品久久久久久精品电影小说| 看免费成人av毛片| 成人影院久久| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 亚洲成人一二三区av| 大话2 男鬼变身卡| 美女福利国产在线| 捣出白浆h1v1| 久久久久人妻精品一区果冻| 免费少妇av软件| 久久精品国产自在天天线| 狂野欧美激情性bbbbbb| 亚洲av日韩在线播放| av在线播放精品| 午夜免费观看性视频| www.色视频.com| 欧美日韩av久久| 伦理电影大哥的女人| 国产一区亚洲一区在线观看| 天天躁夜夜躁狠狠躁躁| 午夜福利网站1000一区二区三区| 中文字幕最新亚洲高清| 精品一区二区三区视频在线| 久久女婷五月综合色啪小说| 日日爽夜夜爽网站| 精品少妇黑人巨大在线播放| 中文精品一卡2卡3卡4更新| 国产黄色视频一区二区在线观看| 国产免费一级a男人的天堂| 精品午夜福利在线看| 久久99热6这里只有精品| 日本-黄色视频高清免费观看| 国国产精品蜜臀av免费| 国产色婷婷99| 中文欧美无线码| 精品久久久精品久久久| www.av在线官网国产| 七月丁香在线播放| 伊人久久国产一区二区| 欧美xxⅹ黑人| 又黄又粗又硬又大视频| 国产成人精品福利久久| 啦啦啦在线观看免费高清www| 内地一区二区视频在线| 天天操日日干夜夜撸| 啦啦啦中文免费视频观看日本| 赤兔流量卡办理| av免费观看日本| 一二三四在线观看免费中文在 | 精品一区二区三区四区五区乱码 | 久久久久久人妻| 午夜日本视频在线| 欧美人与性动交α欧美精品济南到 | 成人毛片a级毛片在线播放| 亚洲婷婷狠狠爱综合网| 久久久国产欧美日韩av| 女人精品久久久久毛片| 熟妇人妻不卡中文字幕| 中文字幕人妻丝袜制服| 综合色丁香网| 国产高清三级在线| 一区在线观看完整版| 少妇熟女欧美另类| 免费高清在线观看日韩| 亚洲色图综合在线观看| 亚洲,欧美,日韩| 草草在线视频免费看| 999精品在线视频| 1024视频免费在线观看| 99九九在线精品视频| www.av在线官网国产| 国产男人的电影天堂91| 日韩视频在线欧美| 亚洲国产欧美在线一区| 久久久久久人妻| 免费观看a级毛片全部| 国产精品久久久久成人av| 欧美最新免费一区二区三区| 黑人高潮一二区| 精品亚洲乱码少妇综合久久| 国国产精品蜜臀av免费| 精品一品国产午夜福利视频| 啦啦啦在线观看免费高清www| 日韩一区二区三区影片| 高清黄色对白视频在线免费看| 嫩草影院入口| 晚上一个人看的免费电影| 99国产综合亚洲精品| 男女边吃奶边做爰视频| 男人爽女人下面视频在线观看| 国产精品不卡视频一区二区| 在线观看免费日韩欧美大片| 最近手机中文字幕大全| 日韩成人伦理影院| 国产精品一国产av| 中文字幕精品免费在线观看视频 | 欧美日韩视频精品一区| av国产久精品久网站免费入址| 熟妇人妻不卡中文字幕| 亚洲精品乱码久久久久久按摩| 精品久久蜜臀av无| 18禁动态无遮挡网站| 狠狠婷婷综合久久久久久88av| 亚洲综合色惰| 丰满饥渴人妻一区二区三| 尾随美女入室| 亚洲国产精品专区欧美| 亚洲高清免费不卡视频| tube8黄色片| freevideosex欧美| 久久精品久久精品一区二区三区| 精品国产一区二区三区四区第35| 久久精品久久久久久噜噜老黄| 久久久久网色| 国产女主播在线喷水免费视频网站| 校园人妻丝袜中文字幕| 精品酒店卫生间| 日韩成人av中文字幕在线观看| 久久国产精品男人的天堂亚洲 | 欧美 日韩 精品 国产| 国产一区二区三区av在线| 麻豆精品久久久久久蜜桃| 亚洲第一区二区三区不卡| 春色校园在线视频观看| 免费大片黄手机在线观看| av在线播放精品| 不卡视频在线观看欧美| 夫妻午夜视频| 一区二区日韩欧美中文字幕 | 免费看av在线观看网站| 伊人亚洲综合成人网| 国产精品久久久久久精品电影小说| 欧美激情 高清一区二区三区| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 老女人水多毛片| 18在线观看网站| 国产乱来视频区| 91国产中文字幕| 一本大道久久a久久精品| 免费观看在线日韩| 日日摸夜夜添夜夜爱| 女性生殖器流出的白浆| 亚洲高清免费不卡视频| 人人妻人人爽人人添夜夜欢视频| 久久97久久精品| 国产一区有黄有色的免费视频| 国产精品国产三级国产av玫瑰| 国产精品久久久av美女十八| 日本wwww免费看| 香蕉精品网在线| 日韩av不卡免费在线播放| 少妇熟女欧美另类| 欧美日韩视频精品一区| 少妇熟女欧美另类| 免费观看a级毛片全部| 狠狠精品人妻久久久久久综合| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 婷婷色麻豆天堂久久| 亚洲五月色婷婷综合| 精品酒店卫生间| 美女内射精品一级片tv| 中文字幕另类日韩欧美亚洲嫩草| 大片免费播放器 马上看| 99视频精品全部免费 在线| 人人妻人人添人人爽欧美一区卜| 欧美+日韩+精品| 亚洲精品av麻豆狂野| 成年女人在线观看亚洲视频| 精品一区在线观看国产| tube8黄色片| 久久精品久久久久久噜噜老黄| 久久精品国产a三级三级三级| 亚洲色图 男人天堂 中文字幕 | 久久久久久久久久成人| 欧美日韩视频高清一区二区三区二| 日韩,欧美,国产一区二区三区| 伦理电影免费视频| 在线天堂中文资源库| 日韩av免费高清视频| 久久久欧美国产精品| 欧美日韩视频高清一区二区三区二| 两个人免费观看高清视频| 成人漫画全彩无遮挡| 午夜影院在线不卡| 精品一区二区三区四区五区乱码 | 欧美激情国产日韩精品一区| 国产成人a∨麻豆精品| 在线天堂最新版资源| 一二三四中文在线观看免费高清| 中文字幕另类日韩欧美亚洲嫩草| 卡戴珊不雅视频在线播放| 尾随美女入室| 国产成人免费无遮挡视频| 伦精品一区二区三区| 91精品国产国语对白视频| 伦理电影大哥的女人| 人人妻人人添人人爽欧美一区卜| 波野结衣二区三区在线| 中文字幕人妻丝袜制服| 久久97久久精品| 亚洲人成77777在线视频| 日本爱情动作片www.在线观看| 久久女婷五月综合色啪小说| 黄片播放在线免费| 亚洲成人av在线免费| 国产精品久久久久成人av| 成年人免费黄色播放视频| 丰满少妇做爰视频| a级毛片在线看网站| av国产久精品久网站免费入址| 国产精品.久久久| 久久影院123| 青春草视频在线免费观看| 桃花免费在线播放| 自线自在国产av| 日产精品乱码卡一卡2卡三| 下体分泌物呈黄色| 日本色播在线视频| 精品人妻一区二区三区麻豆| 婷婷色综合大香蕉| 久久久a久久爽久久v久久| 97在线视频观看| 久久精品国产鲁丝片午夜精品| 国产成人精品婷婷| 欧美人与善性xxx| 91久久精品国产一区二区三区| 欧美人与善性xxx| 精品国产国语对白av| 伦理电影大哥的女人| 2021少妇久久久久久久久久久| 视频在线观看一区二区三区| 最近的中文字幕免费完整| 日韩中字成人| 日本91视频免费播放| 日日撸夜夜添| 免费看不卡的av| 少妇猛男粗大的猛烈进出视频| 黄色一级大片看看| 免费日韩欧美在线观看| 国产精品成人在线| 超碰97精品在线观看| 色婷婷av一区二区三区视频| 人妻 亚洲 视频| 多毛熟女@视频| 久久鲁丝午夜福利片| 一级毛片我不卡| 亚洲人成77777在线视频| a 毛片基地| 一个人免费看片子| 国产高清国产精品国产三级| 亚洲人成77777在线视频| 熟女人妻精品中文字幕| 国产极品粉嫩免费观看在线| 欧美精品高潮呻吟av久久| 中国三级夫妇交换| 国产精品人妻久久久久久| 91精品伊人久久大香线蕉| 久久亚洲国产成人精品v| 精品卡一卡二卡四卡免费| 永久免费av网站大全| 日本欧美国产在线视频| av播播在线观看一区| 国产在线一区二区三区精| 最新中文字幕久久久久| h视频一区二区三区| 青青草视频在线视频观看| 日本91视频免费播放| 久久久久视频综合| 日本欧美视频一区| 亚洲精品久久午夜乱码| 黄色一级大片看看| 亚洲国产av影院在线观看| 天天躁夜夜躁狠狠躁躁| 久久av网站| 亚洲成人av在线免费| 欧美日韩视频高清一区二区三区二| 高清在线视频一区二区三区| 一区在线观看完整版| 在线 av 中文字幕| 母亲3免费完整高清在线观看 | 国产精品久久久久久精品电影小说| 亚洲性久久影院| 丰满乱子伦码专区| 国产在线一区二区三区精| 草草在线视频免费看| 亚洲精品成人av观看孕妇| 亚洲精品日韩在线中文字幕| 国产免费视频播放在线视频| 五月开心婷婷网| 国产成人精品婷婷| 亚洲精品久久成人aⅴ小说| 有码 亚洲区| 国产高清三级在线| xxx大片免费视频| 男人舔女人的私密视频| 亚洲av欧美aⅴ国产| 王馨瑶露胸无遮挡在线观看| 黑丝袜美女国产一区| 久久99热这里只频精品6学生| 在线精品无人区一区二区三| 99久久中文字幕三级久久日本| 在线观看人妻少妇| 香蕉国产在线看| 午夜精品国产一区二区电影| 亚洲伊人久久精品综合| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到 | 欧美另类一区| 国产男女超爽视频在线观看| 亚洲高清免费不卡视频| 超色免费av| 丁香六月天网| 午夜激情久久久久久久| 亚洲久久久国产精品| 国产精品嫩草影院av在线观看| 日本午夜av视频| 久久人人97超碰香蕉20202| av一本久久久久| 国产av码专区亚洲av| 亚洲欧美清纯卡通| 亚洲美女黄色视频免费看| 欧美97在线视频| 2018国产大陆天天弄谢| 国产成人精品在线电影| 在线免费观看不下载黄p国产| 国产精品无大码| 国产免费又黄又爽又色| 卡戴珊不雅视频在线播放| 夫妻性生交免费视频一级片| 天堂俺去俺来也www色官网| 各种免费的搞黄视频| 涩涩av久久男人的天堂| 少妇的逼好多水| 免费av不卡在线播放| 国产日韩欧美视频二区| 亚洲色图 男人天堂 中文字幕 | 亚洲在久久综合| 成年人午夜在线观看视频| 搡女人真爽免费视频火全软件| 久久久精品免费免费高清| 午夜免费鲁丝| 成人毛片a级毛片在线播放| 国产在线一区二区三区精| 国产男人的电影天堂91| 国产免费现黄频在线看| 国产亚洲一区二区精品| 日韩三级伦理在线观看| 18+在线观看网站| 亚洲av免费高清在线观看| 人人妻人人澡人人看| 少妇精品久久久久久久| 啦啦啦中文免费视频观看日本| 久久免费观看电影| 91成人精品电影| 欧美另类一区| 精品熟女少妇av免费看| 中文字幕制服av| 亚洲成国产人片在线观看| 国产一区二区三区综合在线观看 | 尾随美女入室| a级毛片黄视频|