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

    Experimental Research on Mechanism of Hull Plates Curved Forming by a Clean Energy Source

    2021-12-31 07:52:42,-,,,,-,-,
    船舶力學(xué) 2021年12期

    ,-,,,,-,-,

    (Ship Intelligent Manufacturing and Ship Intelligence Integrated Laboratory,Jiangsu University of Science and Technology,Zhenjiang 212003,China)

    Abstract: Line heating is the main practical method for manufacturing curved surfaces in the outer hull plates of a ship. Normally, oxyacetylene flames are used as the heat source. However, performing the line heating method with oxyacetylene flames has challenges in most of shipyards, including large energy consumption, heavy pollution and low automation degree. In order to solve the aforementioned problem,it was intended to utilize the oxyhydrogen flame as the non-emission heat source.Many studies on the deformation mechanism and law of plates were performed through establishing mechanism models and finite element numerical simulation.Then experiments were carried out to evaluate the accuracy and availability of models and energy-saving effect.Different aspects,including the determination of processing parameters,verification of the temperature field and deformation field, and economical comparisons between oxyhydrogen and acetylene flames, were investigated in this regards. It was found that the results obtained from the finite element analysis had a good agreement with that of the experiment.Moreover,it was observed that the processing effect and efficiency of the oxyhydrogen flame met the requirements of curved surface deformation process of the plates.The present study shows that the overall expense of the oxyhydrogen flame is fairly less than that of the oxyacetylene flame.

    Key words:line heating;curved surface deformation;oxyhydrogen flame;clean and non-emission energy

    0 Introduction

    The outer hull surface of a ship is commonly composed of a complex and non-expandable three-dimensional shape with special curves[1]. The marine steel plates are processed into the designed curved shape by means of the dieless forming. At present, shipbuilding companies all over the world generally adopt the line heating method in this regards (Terms‘the process’and‘the method’refer hereafter to‘the ship’s outer surface forming process’and‘the line heating method’in the ship manufacturing process respectively)[2]. It is widely used in shipyards around the world due to its fast processing, flexible operation, no need for other equipment and suitability for the complex forming processing. It should be indicated that‘the method’belongs to the range of the thermal stress forming,whose heating process schematic diagram is shown in Fig.1

    Fig.1 Schematic diagram of plate forming by the line heating process

    However, in view of the long-term application and modern ship manufacturing requirements,the oxyacetylene flame has considerable drawbacks, including low heating efficiency, heavy environmental pollution,serious harm to workers’health and high operational costs.

    Many researchers have tried to replace oxyacetylene with other sources of heat.Siqueira[3]studied the laser forming of high-strength aluminum alloy plates, and investigated the influence of different powers, scanning speeds and spot diameters on the bending angle of a 1.6 mm thick highstrength aluminum alloy plate through experiments. Shi[4]studied the effect of multi-channel laser forming on the forming accuracy of sheet metal. However, the laser is of a high operational cost so that only low-power lasers are used to form thin plates. In other words, the laser cannot meet the power demand of thick plates for the hot stress forming. Therefore, the laser is not popular in shipbuilding[5].

    Jin[5]studied the effect of high frequency induction heating on diagonal deformation,and established the efficiency of high frequency induction heating system, and compared and analyzed the temperature of measurement and finite element analysis. Li Rui[6]processed 5083 aluminum alloy plates with a high frequency electromagnetic induction heat source, used ANSYS software to carry out finite element modeling and simulation, and studied the influence of the processing time of induction heat source and the pulse current amplitude of electromagnetic induction heat source on the forming effect of ship’s plate.But high frequency induction heating is still in the research stage of forming mechanism,it needs to be studied for a long time.

    With the development of technology,hydrogen as a clean energy,has attracted more and more attention of all walks of life. Oxyhydrogen flame was used as the heat source to heat ship’s plates and its forming mechanism was studied.It is of great significance to the development of ship’s plate processing industry.

    1 Mathematical model of the heating process for hull plates

    1.1 Heat source model of the gas flame

    The Gaussian distribution is used to model the heat flux and study the combustion of the oxyhydrogen gas[7]. The combustion equation of the oxyhydrogen gas is:

    Fig.2 shows the heat flux in accordance with the Gaussian distribution. The applied heat flux density of the oxyhydrogen flame to the surface of the plates varies along the radius.

    Considering the Gaussian distribution,q″ can be expressed as

    Fig.2 Heat flux model of the Gauss distribution

    1.2 Mathematical model of the temperature field

    The oxyhydrogen gas initially (i.e.t= 0) moves from the edge along the pre-designed heating line with a uniform speed ofvf.The temperature fieldT(x,y,z,t)of the plate can be expressed as

    wherea=λ/( )ρ·cp.In the present study,it is assumed that the plate is infinitely large and the heating time is long enough.Then the relative coordinate system on the center of the heat source is used as the center of the moving axis so that the temperature field can be treated as a steady state problem.Consideringx′=x,y′=y-vf t,z′=z,Eq.(4)can be simplified to the following form:

    Eq.(5)should satisfy the following three boundary conditions:

    (1)Initial condition:t= 0;

    (2)T(x,y,z,0 )=T0,whereT0= 23 ℃is the room temperature;(

    3)Input condition of the oxyhydrogen gas.

    It is assumed that the linear equation of the heating line isx= 0,z= 0. Then, the heat flux density of thet-moment heating surfaceq″can be expressed as

    Forced conditions of the heat convection for the air-cooling and the water-cooling are shown as belows:

    When the air-cooling is utilized, the problem of the heat convection involves the heat transfer between the plate and air.Surface coefficients of the plate should be selected accordingly:

    On the other hand,the membrane boiling heat transfer theory should be applied to the process of water-cooling when the plate temperature exceeds 200 ℃.

    1.3 Mathematical model of the deformation field

    In this section,the deformation kinematics of the plate is used to determine the residual deformations[8-9]. The plate is initially flat, while the target surface is the desired shape in accordance with the hull design.The following assumptions are made in the surface forming process of the plate:

    (1)The thickness of the plate in the deformation process is constant;(2)Surface shrinkage and bending deformation occur during the plate deformation process;(3)The shear deformation is negligible so that it is ignored in the present study.

    1.4 Finite element theory

    1.4.1 Finite element theory for the temperature field

    The finite element solution steps for the temperature field are as follows[10]:

    (a)Discretization:The temperature fieldT(x,y,z,t)is discretized within the plate as

    Eq.(9) can be combined with the differential equation and the corresponding boundary condition.Applying the finite element approach to the heat source term of the gas flame results in the following expression for the conducted heat

    When the temperature field is solved, the boundary heat transfer should also be considered.The expressionQsis the boundary heat loss in the form shown below:

    In practical and engineering applications, the correlation between the heat transfer coefficient and the temperature is expressed in a piecewise linear model. Moreover, studies show that material characteristics vary as the temperature changes. If{ }Teshows the temperature of a node at timet{ }Te,thermophysical properties of the node can be calculated as the following:

    1.4.2 Finite element theory of the deformation field

    The finite element solution for the deformation field of the plate in "the method" is shown as belows[11-12]:

    Similar to the finite element approach of the temperature field, the displacement field increment Δu(x,y,z,t)is discretized on the plate as below

    ‘The method’can be summarized in a three-dimensional,transient,and thermo-elastic-plastic problem. It should be indicated that an appropriate time step should be selected to simulate the deformation process.

    2 Experiment of oxyhydrogen flame

    The research on the oxyhydrogen flame is composed of determination of processing parameters,verification of the temperature field,and deformation field.

    2.1 Experimental device

    The experimental devices consist of such basic measuring tools as the oxyhydrogen machine,infrared thermal imager,water-cooling device,welding torch,tape measure and the vernier caliper.

    The required oxyhydrogen gas for the experiment is produced by the oxyhydrogen machine.Fig.3 shows the appearance of the oxyhydrogen machine and Fig.4 shows the appearance of the infrared thermal imager. This machine effectively electrolyzes the water and produces a remarkable amount of oxyhydrogen gas, and the gas flux rate can be read from the HMI. In the experiment, an infrared thermal camera (Ruice AI50) with a measuring accuracy of ±2 ℃is utilized to collect the temperature data of the plate during the heating process.

    Fig.3 Oxyhydrogen machine

    Fig.4 Infrared thermal imager

    2.2 Calculating parameters of the oxyhydrogen flame processing technology

    A key problem in performing the finite element analysis for the temperature and deformation fields of the plate is calculating the corresponding processing parameters. In the present study, the Gaussian distribution is used to model the heat source. There are two parameters in this model, the radiusr0of the heat source and heating efficiencyη,which are normally identified by other parameters such as the flux of hydrogen gasQ,the height of the nozzleHand the moving speed of the heat source[13]. Therefore,ηandr0must be initially determined to utilize the model in the actual test case.

    2.2.1 Calculating the processing parameters

    In Eq.(2), many parameters, such as the nozzle diameter, nozzle height, flame movement speed,plate thickness and the material properties[14-18],can affectηandr0.In‘the process’,the effect of other factors except the gas flux is negligible and can be ignored for the optimal value ofH.Therefore,the values ofηandr0mainly depend on the amount of gas flux.

    In order to verify the temperature field of the oxyhydrogen flame,the optimal values ofHandr0should be initially measured.Then the value ofηis modified to adjust the finite element model until the numerical results are consistent with the experimental ones.Finally,the obtained optimal value ofηshows the heating efficiency of the corresponding flux.

    2.2.2 Experimental scheme

    (1)Determining the optimal value ofH

    In order to initially obtain the correlation between the hydrogen fluxQH2,ηandr0, the optimal value ofHshould be determined. The same testing plate is processed for different gas flux rates in this regard.The surface temperatureT1of the plate is measured by the infrared measuring device in the same heating time. The optimal height for different flux rates can be obtained through comparison of the results.

    (2)Calculatingηandr0

    After determining the optimal value ofH, the value ofr0for different hydrogen gas fluxesQH2can be measured. Moreover, the value ofηis obtained by finite element analysis. Six groups of different gas fluxes are taken within the gas production ranges of the oxyhydrogen machine, which are combined with the actual process to give six sets of processing parameters.

    2.2.3 Experimental results

    Tab.1 shows the experimental data of optimal value ofHfor the oxyhydrogen flame.

    Tab.1 Optimum height

    The experiment is carried out by the above mentioned scheme for the optimal value ofH, and the value ofr0is measured.

    The value ofηis modified to adjust the finite element model until the numerical results are consistent with the experimental ones. Then the heating efficiencyηfor different hydrogen flux rates is obtained.Tab.2 shows the specific data.

    Tab.2 Experimental determination of η and r0

    Tab.2 shows the values ofQH2,ηandr0for different fluxes.The values ofηandr0are constant at 1 533(L/h)of the hydrogen flux.With the change of the speed,the finite element numeration temperatures of each group are compared with the measuring temperature. The value ofηis verified through comparison.Tab.3 illustrates the experimental results.

    Tab.3 Experimental verification

    Tab.3 shows that there exist errors in the measurement and numerical analysis, and the error range between the experimental data and the computed data is acceptable. Therefore, it is concluded that the obtained heating efficiency by this experimental scheme is effective and reliable.

    2.3 Result verification of the oxyhydrogen flame deformation field

    2.3.1 Experimental scheme

    The verification steps of the deformation field results are as the following:

    (1)Mark four sets of measuring points on the upper and lower surfaces around the heating line as illustrated in Fig.5;

    (2)Measure the three-dimensional coordinate data of the measuring points before and after deformation of the plates;

    (3) Employ Eqs.(21)-(23) to calculate ΔLandθin accordance with the deformation principle.If the numerical simulation data approach to the experimental data, it is proved that the numerical simulation is valid.

    Fig.5 Marked ship points on the measuring plateFig.6 Schematic principle of the deformation

    Fig.6 shows schematic principle of the deformation and the principle of the numeration method is as follows:

    2.3.2 Measurement and analysis of the deformation field data of the oxyhydrogen flame

    The data of the deformation field are mainly obtained by calculating the displacement of measuring points[19].The data before and after point deformation of four groups are measured accordingly. According toL0,LTandLB, values of the line and angular deformations corresponding to each set are calculated.Then they are averaged to obtainLandθ.

    Tab.4 shows the comparison of the obtained numerical results through ANSYS with the ones of the experiment.

    Tab.4 Comparison of plate deformation results

    The results indicate that the corresponding errors of ΔLandθbetween the ANSYS simulation and the experimental results are 4% and 14% respectively[20-23]. The measurement process and the indirect process of the measured values produce different degrees of errors. It is concluded that results of the deformation experiment are reliable within the allowable range[4,6].

    2.4 Discussion on temperature field results

    According to the established finite element model for the temperature field of the plate,the numerical model of the ANSYS software results in the temperature field distribution.

    Contours of abovementioned temperatures show the range and trend of the temperature variation in the plate and the movement of the heat source in the heating process[24]. Fig.7(a) shows that during the heating phase of‘the method’, the high temperature region is concentrated near both ends of the heating line.Moreover,Fig.7(b)indicates that at the end of the heating process,the heat input spreads from the region near the heating line (high temperature region) to the surrounding region (cryogenic region). Furthermore, Fig.7(c) shows that when the plate is cooled, the heat continues to spread while the temperature of the plate continues to reduce.Eventually,the plate temperature approaches to the steady state value.

    Fig.7 Temperature distribution at 2 500(L/h)flux for different moments

    In order to analyze the temperature field distribution during‘the process’,the temperature change of the feature points on the heating line is studied and three feature points on the heating line, including the start, middle and end points, are selected. Fig.8 shows the three feature points. It is observed that the maximum temperatures in the three points are different due to the asymmetry of the material between the two ends and the middle point. It is found that the maximum temperatures of the plate surface at the start, end and middle points are about 670 ℃,800 ℃and 840 ℃respectively.

    Fig.8 Temperature variation trend of three feature points

    2.5 Discussion on the results of deformation field

    According to the finite element model for the deformation field of the plate, the displacement distribution for different moments is generated in‘the process’by ANSYS simulation.Figs.9(a)-(d)illustrate the displacement contours for the start, middle, end and cooling phases of‘the process’.Moreover,Fig.10 shows the overall deformation effect.

    Fig.9 Displacement distribution at 2 500(L/h)flux for different moments

    Fig.10 Overall deformation result

    Fig.11 Displacement variation trend of three feature points

    Fig.11 shows that thez-displacement change history of feature points at the plate is different.The obtained results in different positions of the displacement and the bending angle are not the same, which verifies the existence of the‘uneven deformation’phenomenon[12-13]. Uneven deformation brings challenges to the precise control of‘the process’.Moreover,Figs.9-10 indicate that the deformation effect of the plate with single heating line is small.

    3 Comparison and analysis of two flames

    In order to further compare the differences between the two flames,the temperature and deformation field of the two flames were established by finite element simulation, and the surface maximum temperature was made similar by adjusting heating speed at the same flow rate. Then the differences of heating width, heating depth, angular deformation and line deformation were compared and analyzed. Furthermore, in order to verify the energy-saving effect, the same plates were processed with similar deformation requirements by using the oxyhydrogen and oxyacetylene flames.Then the energy consumption was computed and the results were compared with each other to verify the energy-saving effect.

    3.1 Comparison of the two flames

    3.1.1 Temperature field contrast

    In previous studies, the surface temperature of a hull plate should be controlled between 450 ℃and 850 ℃[25]. If the temperature is higher than 850 ℃, the properties of the metallographic structure for the hull plate will change irreversibly, affecting the performance of the hull plate. If the temperature is lower than 450 ℃,the steel plate will not be deformed significantly.Therefore,in the actual processing, the maximum temperature of the heating surface for a hull plate is often required to be close to 850 ℃, and the maximum temperature of the back of the hull plate can also reach 450 ℃.

    Fig.12 shows the surface temperature field of oxyhydrogen and oxyacetylene gas flame. Fig.13 shows the heating depth of the two flames. Tab.5 shows the effective surface heating area of oxyhydrogen flame heat source is smaller than the oxyacetylene flame heat source,but the flame penetration of oxyhydrogen flame is stronger and the heating depth is deeper.

    Fig.12 Comparison of surface temperature fields

    Fig.13 Comparison of heating depths

    Tab.5 Comparison of two flame effects

    3.1.2 Deformation field contrast

    As can be seen from the data in Tab.6, under the same maximum temperature and flow rate,the deformation generated by the heating of oxyhydrogen flame is slightly smaller than that of oxyacetylene flame. The average linear deformation of oxyhydrogen flame is 87% of oxyacetylene flame,and the average angular deformation is 91%.In combination with the comparison of temperature field above[26], the calorific value and the radius of oxyhydrogen flame heat source are both smaller than those of oxyacetylene flame,so it is necessary to heat plates at a smaller heating speed.Therefore the‘penetrability’of oxyhydrogen flame is stronger.This enables ship’s plates to achieve better heating depth and makes up for the lack of heating width.

    Tab.6 Comparison of two flame deformation effects

    3.2 Energy-saving verification of the oxyhydrogen flame

    3.2.1 Experimental scheme

    According to the finite element simulation case, the oxyacetylene and oxyhydrogen flames are used to process two plates with the same specification and material and the same deformation requirements. The energy consumption of two plates is compared mainly considering the heat source and use cost.

    3.2.2 Energy consumption comparison

    The plates are curved after‘the process’through the actual experiment. Fig.14 illustrates the final processing effect of the oxyacetylene flame while the final processing effect of the oxyhydrogen flame is shown in Fig.15. The comparison of the two flame processing effects is shown in Fig.16.The comparative results indicate that the final deformation effect has a good consistency with each other.Then the energy consumption is computed from the viewpoint of raw material and operational costs.

    Fig.14 Plate deformation by oxyacetylene flame

    Fig.15 Plate deformation by oxyhydrogen flame

    Raw material cost numeration:

    According to the relevant statistics, each cubic meters of the acetylene consumes 21.8 kW of electricity.The production of the same amount of hydrogen requires 4.8 kW of electricity. In other words, the direct use of the hydrogen gas rather than using the acetylene, reduces the electricity costs by up to 22%.

    Use-cost numeration:

    According to the usage statement, two bottles of acetylene and one bottle of oxygen are consumed in the oxyacetylene processing experiment. The time consumption of the mixed oxyhydrogen gas processing experiment is 6 h. Tab.7 shows the numeration of costs.

    Fig.16 Comparison of plate deformations of two schemes

    Tab.7 Comparison of consumption costs

    According to the relevant literature, the oxyhydrogen flame is not only zero emission but also has obvious energy-saving effects. With the same processing effect, the cost of the oxyhydrogen flame is 66.6%lower than that of the oxy-acetylene flame on average.

    Energy consumption comparison shows that the oxyhydrogen machine can fully meet the process requirements and can also save energy costs.It is unnecessary to change the bottle when using the oxyhydrogen flame, which can save the bottle changing time and improve the work efficiency.Moreover, water is used as the raw material by electrolysis of water. After combustion of the oxyhydrogen flame,the oxyhydrogen gas reverts to the water,so no waste gas is generated and there is no pollution to the environment.So it is environment-friendly and safe.

    4 Concluding remarks

    This paper presents the use of the oxyhydrogen flame as the heat source for the first time in‘the process’.The experiments and ANSYS simulation show that it is feasible to use the oxyhydrogen flame as the heat source for‘the process’.Moreover,compared with the conventional oxyacetylene flame,the oxyhydrogen flame is more efficient and produces no carbon dioxide emissions.Furthermore,experimental data show that the operational cost of the oxyhydrogen flame is 66.6% lower than that of the acetylene flame. It is believed that the efficiency of the oxyhydrogen flame is only higher than the oxyacetylene flame if it is applied on a large scale in industrial applications.The research provides support for the subsequent automated processing of‘the process’.

    国产成人精品在线电影| 日本-黄色视频高清免费观看| 精品福利永久在线观看| 看非洲黑人一级黄片| 国产精品 欧美亚洲| videossex国产| 亚洲欧美精品综合一区二区三区 | 99国产精品免费福利视频| 国产淫语在线视频| 日韩电影二区| 精品视频人人做人人爽| 人妻 亚洲 视频| 免费在线观看完整版高清| 久久久久国产网址| 国产免费现黄频在线看| 国产一区二区激情短视频 | 五月伊人婷婷丁香| 国产麻豆69| 国产男人的电影天堂91| 在线观看国产h片| av女优亚洲男人天堂| 亚洲第一av免费看| 最近最新中文字幕大全免费视频 | 两性夫妻黄色片| av一本久久久久| 美女国产高潮福利片在线看| 久久毛片免费看一区二区三区| 最近最新中文字幕大全免费视频 | 少妇人妻久久综合中文| 高清欧美精品videossex| 精品亚洲乱码少妇综合久久| 国产精品av久久久久免费| 亚洲国产av影院在线观看| 美女xxoo啪啪120秒动态图| 免费黄色在线免费观看| 国产精品二区激情视频| 最黄视频免费看| 亚洲精品久久午夜乱码| 自线自在国产av| 亚洲国产日韩一区二区| 精品久久久久久电影网| 热re99久久国产66热| 天天躁夜夜躁狠狠躁躁| 国产精品av久久久久免费| 久久久久久伊人网av| 亚洲国产欧美网| 香蕉精品网在线| 国产精品二区激情视频| av免费观看日本| 午夜免费鲁丝| 国产在线一区二区三区精| 亚洲精品美女久久久久99蜜臀 | 成人国语在线视频| 亚洲欧美成人综合另类久久久| 黑人欧美特级aaaaaa片| 最近2019中文字幕mv第一页| 九九爱精品视频在线观看| 亚洲精品久久成人aⅴ小说| 成年人免费黄色播放视频| 男女免费视频国产| 美女国产视频在线观看| 亚洲,一卡二卡三卡| 男女边吃奶边做爰视频| 自拍欧美九色日韩亚洲蝌蚪91| 可以免费在线观看a视频的电影网站 | 国产黄色视频一区二区在线观看| 久久精品夜色国产| 三上悠亚av全集在线观看| 性色av一级| 欧美日韩亚洲高清精品| 中文字幕亚洲精品专区| 亚洲欧美日韩另类电影网站| 2021少妇久久久久久久久久久| 国产av国产精品国产| 成人18禁高潮啪啪吃奶动态图| 欧美国产精品一级二级三级| 18禁裸乳无遮挡动漫免费视频| 一区二区三区精品91| 2022亚洲国产成人精品| 亚洲一区二区三区欧美精品| 日韩av不卡免费在线播放| 免费女性裸体啪啪无遮挡网站| 欧美+日韩+精品| 精品视频人人做人人爽| 亚洲欧美色中文字幕在线| 久久久精品区二区三区| 精品卡一卡二卡四卡免费| 国产成人免费无遮挡视频| 亚洲精品久久午夜乱码| av国产精品久久久久影院| 各种免费的搞黄视频| 99久久人妻综合| 亚洲精品国产一区二区精华液| 日韩熟女老妇一区二区性免费视频| 飞空精品影院首页| 国产精品 国内视频| 在线精品无人区一区二区三| 考比视频在线观看| 久久久久久久国产电影| 久久久久久久国产电影| 少妇人妻精品综合一区二区| 男男h啪啪无遮挡| 久久人人爽av亚洲精品天堂| 亚洲国产欧美网| av电影中文网址| 人人妻人人添人人爽欧美一区卜| 我的亚洲天堂| 日日摸夜夜添夜夜爱| 亚洲国产欧美网| 国产老妇伦熟女老妇高清| 国产老妇伦熟女老妇高清| 一区二区三区乱码不卡18| 2018国产大陆天天弄谢| 欧美中文综合在线视频| 国产一区二区 视频在线| 最近中文字幕高清免费大全6| 亚洲av.av天堂| 1024视频免费在线观看| 中文精品一卡2卡3卡4更新| 有码 亚洲区| 国产成人av激情在线播放| 国产成人免费无遮挡视频| 久久国内精品自在自线图片| 国产精品免费视频内射| 美国免费a级毛片| av国产久精品久网站免费入址| 午夜福利,免费看| 九草在线视频观看| 一级毛片 在线播放| 国产成人aa在线观看| 日本av手机在线免费观看| 久久久国产一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲精品视频女| 男的添女的下面高潮视频| 伊人亚洲综合成人网| 美女大奶头黄色视频| 国产福利在线免费观看视频| 国产熟女欧美一区二区| 亚洲成人av在线免费| 十八禁网站网址无遮挡| 看免费成人av毛片| 欧美黄色片欧美黄色片| 亚洲成av片中文字幕在线观看 | 不卡视频在线观看欧美| 又大又黄又爽视频免费| 日韩中字成人| 亚洲国产精品一区二区三区在线| 日韩熟女老妇一区二区性免费视频| 久久精品国产综合久久久| 国产成人精品无人区| 久久精品国产亚洲av涩爱| 十八禁高潮呻吟视频| 亚洲图色成人| 国产成人精品一,二区| 男女午夜视频在线观看| 狂野欧美激情性bbbbbb| 自拍欧美九色日韩亚洲蝌蚪91| 少妇人妻久久综合中文| 男人操女人黄网站| 国产 一区精品| 搡女人真爽免费视频火全软件| 看免费av毛片| 亚洲美女搞黄在线观看| 亚洲第一青青草原| 久久久久久久久久久久大奶| 精品一区二区免费观看| 国产精品久久久久成人av| 国产野战对白在线观看| 人成视频在线观看免费观看| 亚洲国产精品一区三区| 久久久久久人妻| 久久久久久久久久久免费av| 十分钟在线观看高清视频www| 一级毛片我不卡| 国产日韩欧美亚洲二区| 天天操日日干夜夜撸| av网站在线播放免费| 亚洲精品日本国产第一区| 国产免费又黄又爽又色| 国产成人精品无人区| 国产精品三级大全| 999精品在线视频| 99re6热这里在线精品视频| 国产 一区精品| av在线播放精品| 男人舔女人的私密视频| 精品国产一区二区三区久久久樱花| 国产97色在线日韩免费| 国产成人aa在线观看| 欧美成人精品欧美一级黄| 国产一区二区三区综合在线观看| 一级毛片 在线播放| 男人爽女人下面视频在线观看| 亚洲第一av免费看| 国产极品天堂在线| 啦啦啦在线免费观看视频4| 国产精品久久久av美女十八| 亚洲精品自拍成人| 免费观看性生交大片5| 成年av动漫网址| 波野结衣二区三区在线| 亚洲内射少妇av| 最近的中文字幕免费完整| 国产精品国产av在线观看| 1024视频免费在线观看| 亚洲av福利一区| 老司机影院成人| 午夜免费观看性视频| 一级,二级,三级黄色视频| 少妇的丰满在线观看| 久久久久视频综合| 人人妻人人澡人人看| 在线观看三级黄色| 欧美少妇被猛烈插入视频| 又粗又硬又长又爽又黄的视频| 日韩不卡一区二区三区视频在线| 亚洲成国产人片在线观看| 国产高清国产精品国产三级| 日日爽夜夜爽网站| 国产成人aa在线观看| 午夜福利,免费看| 欧美日韩视频精品一区| 午夜日韩欧美国产| 久热这里只有精品99| 国产精品偷伦视频观看了| 亚洲 欧美一区二区三区| 国产乱人偷精品视频| 国产成人精品久久二区二区91 | 91午夜精品亚洲一区二区三区| 国产日韩一区二区三区精品不卡| 如何舔出高潮| 美女脱内裤让男人舔精品视频| 国产亚洲欧美精品永久| 日日啪夜夜爽| 免费高清在线观看视频在线观看| 丰满少妇做爰视频| 国产在线免费精品| 亚洲精品国产一区二区精华液| 精品少妇一区二区三区视频日本电影 | 美女国产视频在线观看| 一二三四在线观看免费中文在| 国产成人精品在线电影| 欧美中文综合在线视频| 久久久久久人妻| 天堂中文最新版在线下载| 在线观看美女被高潮喷水网站| 丝瓜视频免费看黄片| 日本-黄色视频高清免费观看| 色哟哟·www| 国产爽快片一区二区三区| 黑丝袜美女国产一区| 久久久国产一区二区| 妹子高潮喷水视频| 午夜福利在线观看免费完整高清在| 日韩精品有码人妻一区| 亚洲 欧美一区二区三区| 日韩制服丝袜自拍偷拍| 尾随美女入室| 99热国产这里只有精品6| 春色校园在线视频观看| 国产精品不卡视频一区二区| 久久久久久伊人网av| 日韩成人av中文字幕在线观看| 天天躁夜夜躁狠狠躁躁| 国产黄色免费在线视频| 色视频在线一区二区三区| 免费少妇av软件| 两个人免费观看高清视频| 美女大奶头黄色视频| 亚洲精品久久午夜乱码| 亚洲国产看品久久| 日韩精品有码人妻一区| 男人爽女人下面视频在线观看| 亚洲在久久综合| 可以免费在线观看a视频的电影网站 | 日日爽夜夜爽网站| 国产又爽黄色视频| 亚洲美女黄色视频免费看| 王馨瑶露胸无遮挡在线观看| 亚洲精华国产精华液的使用体验| av片东京热男人的天堂| 在线观看www视频免费| 国产男女超爽视频在线观看| 精品一区在线观看国产| 丝袜喷水一区| 五月伊人婷婷丁香| 97人妻天天添夜夜摸| 黄色配什么色好看| 91久久精品国产一区二区三区| 午夜精品国产一区二区电影| 91午夜精品亚洲一区二区三区| 丝袜美足系列| 国产又色又爽无遮挡免| 大码成人一级视频| 免费黄网站久久成人精品| 母亲3免费完整高清在线观看 | 97在线人人人人妻| 欧美精品一区二区免费开放| 中文字幕精品免费在线观看视频| 亚洲久久久国产精品| 桃花免费在线播放| 热re99久久精品国产66热6| 老司机影院毛片| 91久久精品国产一区二区三区| www.自偷自拍.com| 两性夫妻黄色片| 亚洲欧美一区二区三区久久| 久久精品aⅴ一区二区三区四区 | 国产成人免费无遮挡视频| h视频一区二区三区| 最黄视频免费看| a级毛片黄视频| www.av在线官网国产| 国产男女内射视频| 女人高潮潮喷娇喘18禁视频| 精品人妻偷拍中文字幕| 午夜91福利影院| 另类亚洲欧美激情| 色吧在线观看| 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 少妇人妻精品综合一区二区| 国产白丝娇喘喷水9色精品| 欧美少妇被猛烈插入视频| 99久久人妻综合| 2022亚洲国产成人精品| 老司机影院成人| 99国产精品免费福利视频| 性少妇av在线| 日韩中文字幕欧美一区二区 | 国产精品久久久久久久久免| 久久精品国产亚洲av高清一级| 色视频在线一区二区三区| 精品一区二区三卡| 七月丁香在线播放| 美女视频免费永久观看网站| 国产成人精品无人区| 9191精品国产免费久久| 亚洲精品美女久久久久99蜜臀 | 日韩一本色道免费dvd| 咕卡用的链子| 久久精品国产亚洲av涩爱| 黑人猛操日本美女一级片| 久久久a久久爽久久v久久| 日本色播在线视频| 秋霞在线观看毛片| 国产成人精品久久二区二区91 | 一区二区三区四区激情视频| 考比视频在线观看| 成人国语在线视频| 国产爽快片一区二区三区| 大香蕉久久成人网| 久久久久网色| 国产男女超爽视频在线观看| 校园人妻丝袜中文字幕| 日韩一区二区视频免费看| 亚洲一级一片aⅴ在线观看| 18禁动态无遮挡网站| 老鸭窝网址在线观看| 亚洲国产最新在线播放| 国产欧美亚洲国产| xxxhd国产人妻xxx| 国产成人一区二区在线| 丰满乱子伦码专区| 天天躁狠狠躁夜夜躁狠狠躁| 成人国产麻豆网| 成人二区视频| 青草久久国产| 免费看av在线观看网站| 男女边摸边吃奶| 亚洲婷婷狠狠爱综合网| 美女午夜性视频免费| 男女下面插进去视频免费观看| 精品久久蜜臀av无| 国产精品国产三级国产专区5o| 午夜福利一区二区在线看| 中文字幕人妻熟女乱码| 亚洲综合精品二区| 美女午夜性视频免费| 免费高清在线观看视频在线观看| 亚洲情色 制服丝袜| 国产片特级美女逼逼视频| 日韩欧美一区视频在线观看| 日韩一本色道免费dvd| 久久99热这里只频精品6学生| 亚洲伊人色综图| 天堂中文最新版在线下载| 精品国产一区二区久久| 99久久精品国产国产毛片| 久久毛片免费看一区二区三区| 欧美日韩一级在线毛片| 最近手机中文字幕大全| 伊人久久国产一区二区| 免费女性裸体啪啪无遮挡网站| 国产高清国产精品国产三级| 亚洲成人一二三区av| 欧美日韩综合久久久久久| 综合色丁香网| 十八禁网站网址无遮挡| 黄片播放在线免费| 国产伦理片在线播放av一区| 亚洲精品av麻豆狂野| 亚洲男人天堂网一区| 国产片内射在线| 精品一区二区三卡| 亚洲综合色网址| 一个人免费看片子| 日日撸夜夜添| 自线自在国产av| 久久久久精品久久久久真实原创| 青春草国产在线视频| 天天操日日干夜夜撸| 超碰成人久久| 肉色欧美久久久久久久蜜桃| 久久午夜福利片| 精品少妇黑人巨大在线播放| 国产精品一区二区在线不卡| 秋霞伦理黄片| 国产激情久久老熟女| 久久久a久久爽久久v久久| 午夜久久久在线观看| 99久久人妻综合| 日韩精品有码人妻一区| 精品99又大又爽又粗少妇毛片| 制服诱惑二区| 免费看av在线观看网站| 91午夜精品亚洲一区二区三区| 超碰成人久久| 叶爱在线成人免费视频播放| 国产成人a∨麻豆精品| 国产一区二区激情短视频 | 日韩中文字幕视频在线看片| 免费黄频网站在线观看国产| av免费在线看不卡| 97在线人人人人妻| 极品人妻少妇av视频| 精品99又大又爽又粗少妇毛片| 久久精品国产自在天天线| 亚洲国产成人一精品久久久| 69精品国产乱码久久久| 久久人人97超碰香蕉20202| 一本—道久久a久久精品蜜桃钙片| 精品国产一区二区三区四区第35| 91精品三级在线观看| 国产激情久久老熟女| 国产精品香港三级国产av潘金莲 | 亚洲,欧美,日韩| 久久综合国产亚洲精品| 天天影视国产精品| 国产男人的电影天堂91| 亚洲国产看品久久| 亚洲精品国产一区二区精华液| 成年动漫av网址| 国产成人精品婷婷| 国产亚洲精品第一综合不卡| 欧美精品国产亚洲| 欧美变态另类bdsm刘玥| 久久ye,这里只有精品| 国产精品不卡视频一区二区| 久久久久久久久免费视频了| 伦精品一区二区三区| 午夜免费鲁丝| 国产成人精品福利久久| 91精品国产国语对白视频| 大香蕉久久网| 少妇猛男粗大的猛烈进出视频| 国产一区二区激情短视频 | 精品99又大又爽又粗少妇毛片| 免费播放大片免费观看视频在线观看| 免费黄网站久久成人精品| 一区二区av电影网| 男女午夜视频在线观看| 各种免费的搞黄视频| 18禁动态无遮挡网站| 成人国产av品久久久| av免费在线看不卡| 丁香六月天网| 国产男人的电影天堂91| 日韩一区二区视频免费看| √禁漫天堂资源中文www| 一区二区日韩欧美中文字幕| 丰满迷人的少妇在线观看| 亚洲在久久综合| 久久97久久精品| 色94色欧美一区二区| 一区二区三区乱码不卡18| 99国产综合亚洲精品| 美国免费a级毛片| 免费观看av网站的网址| 国产免费视频播放在线视频| 亚洲成色77777| 亚洲综合精品二区| 亚洲一码二码三码区别大吗| 少妇 在线观看| 久久精品亚洲av国产电影网| 三上悠亚av全集在线观看| 国产一区二区在线观看av| 男女午夜视频在线观看| 大码成人一级视频| 在线观看免费视频网站a站| 一级毛片电影观看| 国产伦理片在线播放av一区| 国产成人av激情在线播放| 精品人妻在线不人妻| 99热网站在线观看| 久久精品国产a三级三级三级| 精品一区二区免费观看| 欧美精品高潮呻吟av久久| 中文精品一卡2卡3卡4更新| 久久久a久久爽久久v久久| 精品午夜福利在线看| 美女视频免费永久观看网站| 人人妻人人澡人人看| 日本91视频免费播放| 亚洲欧美成人综合另类久久久| 成人二区视频| 少妇人妻 视频| 涩涩av久久男人的天堂| 国产亚洲欧美精品永久| 国产成人午夜福利电影在线观看| 日韩成人av中文字幕在线观看| 黄色一级大片看看| 久久精品久久精品一区二区三区| 青春草视频在线免费观看| 久久精品国产自在天天线| 亚洲国产精品成人久久小说| 久久鲁丝午夜福利片| 亚洲国产av影院在线观看| 性色av一级| 黄频高清免费视频| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 国产一区亚洲一区在线观看| 亚洲精品在线美女| 亚洲欧美精品自产自拍| 色播在线永久视频| 午夜福利乱码中文字幕| 桃花免费在线播放| 精品99又大又爽又粗少妇毛片| 亚洲一区二区三区欧美精品| 嫩草影院入口| 午夜免费观看性视频| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 国产精品无大码| 97人妻天天添夜夜摸| 男女免费视频国产| 久久99热这里只频精品6学生| 在线观看一区二区三区激情| 日本av免费视频播放| 亚洲一区中文字幕在线| 成人毛片a级毛片在线播放| a级毛片黄视频| 人成视频在线观看免费观看| 国产精品嫩草影院av在线观看| 天堂中文最新版在线下载| 午夜91福利影院| 香蕉丝袜av| 亚洲欧美清纯卡通| 亚洲人成77777在线视频| 国产综合精华液| 人妻 亚洲 视频| 大码成人一级视频| 亚洲伊人色综图| 少妇被粗大猛烈的视频| 亚洲av综合色区一区| 这个男人来自地球电影免费观看 | 日韩中字成人| 亚洲成人一二三区av| 国产精品偷伦视频观看了| 在现免费观看毛片| 亚洲中文av在线| 亚洲av电影在线观看一区二区三区| √禁漫天堂资源中文www| 国产一区有黄有色的免费视频| 日韩制服骚丝袜av| 成年av动漫网址| 人人妻人人澡人人看| 国产爽快片一区二区三区| 99久久中文字幕三级久久日本| 亚洲精品久久成人aⅴ小说| 久久久久久久久免费视频了| videosex国产| 老司机影院成人| 国产精品成人在线| 亚洲伊人色综图| 欧美另类一区| av有码第一页| 亚洲欧美中文字幕日韩二区| 如日韩欧美国产精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 啦啦啦在线观看免费高清www| 亚洲,欧美精品.| 丝袜美足系列| 亚洲欧美精品自产自拍| 亚洲精品视频女| 一本—道久久a久久精品蜜桃钙片| 18禁动态无遮挡网站| 99久久中文字幕三级久久日本| 国产日韩欧美在线精品| 综合色丁香网| 高清黄色对白视频在线免费看| 久久影院123| 香蕉丝袜av| 在线观看www视频免费| 精品午夜福利在线看| 久久久久精品人妻al黑| 春色校园在线视频观看| 青草久久国产| 国产精品麻豆人妻色哟哟久久| 91在线精品国自产拍蜜月| 精品人妻偷拍中文字幕| 国产深夜福利视频在线观看| 久久影院123|