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

    Crack Evaluation of an Ultra-High-Pressure Chamber Used for Deep-Sea Environment Simulation

    2021-11-03 14:02:22JiangsuProvincialKeyLaboratoryofAdvancedManufactureandProcessforMarineMechanicalEquipmentJiangsuUniversityofScienceandTechnologyZhenjiang2200ChinaShanghaiEngineeringResearchCenterofMarineRenewableEnergyCollegeofEngineering
    船舶力學(xué) 2021年10期

    ,,,,,4(.Jiangsu Provincial Key Laboratory of Advanced Manufacture and Process for Marine Mechanical Equipment,Jiangsu University of Science and Technology,Zhenjiang 2200,China;2.Shanghai Engineering Research Center of Marine Renewable Energy,College of Engineering Science and Technology,Shanghai Ocean University,Shanghai 2006,China;.Shanghai Engineering Research Center of Hadal Science and Technology,College of Marine Sciences,Shanghai Ocean University,Shanghai 2006,China;4.Westlake University,Hangzhou 0024,China)

    Abstract:Ultra-high-pressure(UHP)thick chamber for simulating the underwater pressure environment bears high cyclic load and is prone to cracks caused by fatigue.Fatigue crack growth is an important factor affecting their fracture.The objective of this paper is to analyze the behavior of semi-elliptical crack growth in an aged UHP chamber and to estimate the safety margin of the chamber.For this purpose,crack growth properties of the material 20MnMoNb steel were experimentally investigated first considering triangular and trapezoidal load history respectively.Two series of results were compared to examine the sensitivity of the material on holding time.A three-dimensional finite element(FE)model was used for fatigue crack growth calculation based on the unified crack growth rate model and verified by comparing groups of numerical and experimental results on CT specimen.A series of FE models for cracks with different initial size,aspect ratio and inclination level are presented and residual life of the chamber is obtained based on the allowable criterion depth of surface cracks on the inner wall of the chamber.The analysis results will be used as a reference to evaluate the reliability of the UHP chamber.

    Key words:ultra-high-pressure chamber;crack growth;pressure vessel;deep-sea environment simulation;residual life

    0 Introduction

    Deep sea exploration attracts human’s attention in recent years,which has led to an upsurge in research and development of deep-sea equipment.Pressure tests on the structure strength and sealing performance are required for components of deep-sea equipment before they are applied in deep-sea to ensure safety.These pressure tests will be carried out in ultra-high-pressure(UHP)thick chambers for simulating the underwater pressure environment with its pressure-resistant structure,as shown in Fig.1 for illustration of a component pressured in the UHP chamber.

    Fig.1 Sketch of UHP chamber for deep-sea environment simulation

    Research has been done on the design and analysis of a variety of UHP chamber systems with the functions of pressure resistance,water-tightness,detecting leaks,video transmission etc.[1-5].The core problem of pressured hull design is to examine the ability of the hull to resist deformation and damage under the action of load[6-8],that is,to guarantee its strength,stiffness and stability,and to ensure the safety.UHP chambers used for deep-sea environment simulation with the form of a typical cylindrical shell are of normal temperature chambers.And the design internal medium pressure may be higher than hundred mega-pascals and much greater than the external pressure.The failure form of a general internal pressure hull such as underwater pipelines includes leakage,rupture and excessive deformation or even destruction[9-10].Strength failure criteria with fatigue and fracture failure included will be adopted to avoid leakage or structural loss stability.

    However,micro-cracks as shown in Fig.2 are sometimes buried below the surface of shell which can expand under severe stress.Strict flaw detection should be performed before use.Undetected defects may grow to new macro-cracks especially near the surface due to the interaction of more rigorous stress processes,which may propagate along the wall thickness direction.Fatigue crack growth properties must be taken into account.

    Fig.2 Sketch of the distributed microcracks inside the material

    Most of studies on crack growth in pressure vessels are carried out in pipe and nuclear industries.Margolin and Karzov[11]applied the fatigue crack propagation model based on the analysis of material behavior in the local zones near the crack tip on nuclear pressure vessels made of 2CrMoNiV.Price and Ibrahim[12]studied crack growth in aluminum cylinders with cracks initiating at the neck shoulder of nuclear vessels.Meaningful results are obtained for damage tolerance analysis.Tan et al[13]presented creep constraint and fracture parameter for axial semi-elliptical surface cracks with high aspect ratio in pressurized pipes made of 25Cr2NiMo1V steel,in which creep life assessment method was provided.Various studies can be found in other literatures[14-17].The researchers of these studies have grasped some methods and laws of crack propagation in pressure vessels from different levels,but they have not paid enough attention to ultra-high pressure thickwalled vessels.

    The UHP chambers equipped in the authors’laboratory have been in use for more than five years.In order to ensure the safety of the pressure chambers,it is necessary to evaluate their performance in a comprehensive way.At present,it is difficult to detect the inner walls of the UHP chambers totally.However,according to past experience,there will inevitably be some severe deformation and surface defects in the inner walls after cycles of high pressure.These surface defects can be conservatively treated as surface cracks.The evaluation of the cyclic capacity of the cracked bodies can be specified.20MnMoNb steel is a material which is usually used for construction of such chambers.But the crack growth data are not enough.In this paper,the crack growth properties of the material under two types of loading conditions are experimentally studied to examine the sensitivity of the material on holding time.Then a numerical study on a certain UHP chamber with semielliptical cracks in the internal side is established to evaluate its safety against fracture and fatigue failure under acceptable crack size criterion.

    1 Problem description

    1.1 Structure and material

    The analyzed UHP chamber has an inner diameter of 450 mm and an outer diameter of 890 mm,which was designed for simulating the full ocean depth environment with a pressure of 115 MPa.Its configuration is shown in Fig.3.It is a cylindrical thick-walled structure with a flat cover.

    Fig.3 Configuration of UHP chamber

    The chamber is made of 20MnMoNb steel.The chemical composition ratio of the material is listed in Tab.1.Four tensile samples tested to obtain the mechanical properties of the material are listed in Tab.2,showing a satisfactory plasticity and high strength.

    Tab.1 Chemical composition of material 20MnMoNb steel(weight%)

    Tab.2 Material mechanical properties

    1.2 Loading history

    In the course of service,deep-sea equipment will bear the pressure of sea water,and the pressure will change with the process of dive-operation-float for many times,that is to say,it will bear the effect of cyclic load.Both the static strength and fatigue strength of the equipment need to meet the requirements.In the process of strength checking and fatigue checking of these components by means of a UHP chamber,the chamber itself is also affected by cyclic load.In addition,it is sometimes necessary to consider whether the retaining load during the operation in deep sea has an impact on the component.Fig.4 records the pressure of the chamber analyzed in the present study during a certain period of time.Retaining load exists in some loading cycles for several minutes or even hours.Generally,in fatigue analysis,the treatment of similar trapezoidal loads can be approximated by linear superposition method.Crack propagates under the joint action of triangular cyclic load and retaining load.The crack growth may be different under the two load sequences respectively,as shown in Fig.5(a)and Fig.5(b).But the extent of difference also depends on the properties of the material itself and the duration of magnitude of the sustaining load.Research on the dwell fatigue properties of several candidate high strength alloy materials used in deep sea equipment were conducted by Wang et al[18-19].Dwell fatigue properties of the present material for UHP chambers 20MnMoNb steel have not been reported yet.Fatigue life analysis based on fracture mechanics needs considering the effect of sustaining load on crack growth performance of the material which is presented in the next section.

    Fig.4 Loading history of the ultra-high-pressure chamber(with an inner diameter of 450 mm)

    Fig.5 Sketch of two types of simplified loading histories

    1.3 Crack growth properties by experiments

    The crack propagation rate test was carried out on MTS810 fatigue testing machine.The standard compact tensile specimens(hereinafter referred to as CT specimens)withW=60 mm,B=15 mm anda0=12 mm are used to prefabricate cracks of 1-2 mm.The test was carried out in room temperature,using sine wave loading mode with a loading frequency of 10 Hz and the maximum load is 20 kN with a stress ratio of 0.The time of holding load determines the cost of test.To save the test cost,the reference time of retaining load in the present study was set to 2 minutes.

    The crack growth rate da/dNis obtained from thea-Ncurve.The magnitude of the stress intensity factorΔKis calculated as follows:

    whereais the crack length,Wis the sample width,Bis the sample thickness,andΔFis the load magnitude.

    The experimental results of twoa-Ncurves under normal fatigue loading(triangular load history)and dwell fatigue loading(trapezoidal load history)are shown in Fig.6.The crack growth rate curves(da/dN~ΔKcurves)are shown in Fig.7.In the stable growth stage,the two curves almost coincide,and the crack growth is not affected by the sustaining load.

    Fig.6 Crack length versus cycles(a-N curve)under dwell fatigue loading(triangular load history)and normal fatigue loading(trapezoidal load history)

    Fig.7 Crack growth rate curves(da/dN~ΔK curve)under normal fatigue loading(triangularload history)and dwell fatigue loading(trapezoidal load history)

    A unified crack growth rate model(UCGRM)proposed by Cui et al(2013)will be used for crack growth calculation,which can be expressed by the following equations,

    whereAis a material-and environmentally-sensitive constant of dimensions(MPa)-2;mis a constant representing the slope of the corresponding fatigue crack growth rate curve;nis the index indicating the unstable fracture;KIcis the plane strain fracture toughness of the material;KCfis the fracture toughness of the material under fatigue loading;reis an empirical material constant of the inherent flaw length of the order of 1 mm;ais the modified crack length which is equal toreplus the actual crack length;σmaxis the maximum applied stress;σminis the minimum applied stress;Y(a)is a geometrical factor;Y(re)is a geometrical factor whenais equal tore;Ris the stress ratio(=σmin/σmax);ΔKthis the threshold value of stress intensity factor range;athis the threshold value corresponding to ΔKth;ΔKeffis the effective range of the stress intensity factor;ΔKeffthis the effective range of the stress intensity factor at the threshold level;Kopis the stress intensity factor at the opening level;α′is the crack tip stress/strain constraint ratio,which is 1 for the plane stress state and 1/(1-2ν)for the plane strain state.The effect ofnis significant only in the unstable propagation region;a constant value of 6 is recommended for a quick and simple engineering analysis.It is generally assumed thatKCf=KIC.

    Fig.8 is the crack growth rate curve simulated by UCGRM.All the model parameters are listed with the curve.The model will be used for calculating the crack length in 3D simulation.

    Fig.8 Crack growth rate curve simulated by UCGRM with model parameters

    2 Theoretical and numerical analysis

    2.1 Failure criteria

    The failure form of a general internal pressure chamber includes leakage,rupture and excessive deformation or even destruction.Failure criteria can be divided into strength failure criteria and stiffness failure criteria.Strength failure criteria includes elastic failure,plastic failure,elasticplastic failure,fatigue failure,fracture failure and corrosion failure etc.and stiffness failure criteria typically refer to those of leakage and structural loss stability.Sometimes,an acceptable critical dimension is artificially defined to limit the extent of defect expansion before ensuring that the structure does not break.The plasticity and toughness of this 20MnMoNb steel are very strong,especially the fracture toughness reaching more than 120 MPa·m0.5.In the design stage,the deformation of the structure will be checked and the excessive deformation damage of the cylinder will not occur when the structure is used properly.In the stage of crack propagation,high fracture toughness prevents the occurrence of unstable propagation.Therefore,in the service life calculation based on fracture mechanics,the failure criterion proposed in this paper is the criterion of acceptable crack sizes.

    2.2 Modelling approach

    2.2.1 Simulation validation

    A general finite element software Abaqus and a special analytical software for three-dimensional crack propagation Franc3D were combined to simulate the fatigue crack growth of structure for the present study.The accuracy of the stress intensity factors at the crack front determines the accuracy of the simulation results.Therefore,the calculation results of stress intensity factor at crack tip of standard CT specimens were compared with those from Franc3D software firstly.

    Before the finite element calculation,the material properties of CT samples were defined according to the experimental results of 20MnMoNb as listed in Tab.1 and the Poisson’s ratio was 0.3.The surfaces of the two cylindrical holes contacted with the pins were coupled to their respective centers and loaded with a concentrated force of 20 kN at the two points of the coupling.The loads were directed toward both sides of the specimen and the boundary conditions were applied on the central line of the surface farther away from the circular hole.Fig.9 shows the comparison between thea-Ncurves by test and simulation.The calculated results of stress intensity factor at crack tip of CT specimens obtained by Franc3D simulation were compared with those obtained by formula,as shown in Fig.10.According to the comparison of the curves in the figure,it can be seen that the calculated stress intensity factors of Franc3D software are in good agreement with those from the formula except for some deviated points from unstable calculation.Then the model can be used for numerical simulation of fatigue crack growth with satisfactory accuracy.

    Fig.9 Comparison between the a-N curves by test and simulation

    Fig.10 Comparison between theΔK-a curves from test and simulation

    2.2.2 3D simulation of surface cracks analysis

    For a cylindrical internal pressure vessel,the stress on the inner wall of the vessel is the largest,and the internal defects of the material are most likely to occur there.It is assumed that a semielliptical macro-surface crack occurs on the inner wall of the chamber as shown in Fig.11 and the semi-elliptical crack is also a simplified form of the general surface cracks.The crack will be calculated through the sub-model of the cylindrical structure after mesh refinement in Franc3D software.aandcare respectively the short and long half-axis lengths of a semi-elliptical shape,i.e.crack depth and surface half-length.The initial size of the crack,c0=5 mm anda0=3 mm,is presumed for initial analysis.When the maximum value of the cyclic internal pressure is 110 MPa,the variation of crack dimensions with the number of loading cycles can be calculated,as shown in Fig.12.When the number of cycles is 150 000,the crack depth grows to 65 mm and the length reaches 95 mm.The thickness of the cylinder is 220 mm,that is to say,after 150 000 times,the assumed macroscopic crack depth has reached 30% of the thickness of the chamber.

    Fig.11 FE global model and sub-model with a surface crack

    Fig.12 Crack length versus cycles when c0=5 mm and a0=3 mm respectively

    If one tenth of the wall thickness is taken as the acceptable crack depth,the loading cycles are about 110 000 times.The ultimate strength test frequency by using the present pressure cylinder is about 400 times per year,or about 12 000 times in 30 years.If we consider cyclic loading tests of some special components,the frequency of loading cycles can be increased by 200 times a year according to the current frequency records,then it is about 18 000 cycles in thirty years.Under the current criterion that the crack depth is less than one tenth of the wall thickness,the life requirement is satisfied.And according to the analysis,the service life of the pressure cylinder in the design pressure state is about six times the safety margin.

    2.2.3 Effect of initial aspect ratio and crack inclination on crack growth properties

    Fig.14 Effects of initial aspect ratio on crack growth in length direction when a0 is set to 3 mm

    Crack characteristics are one of the factors affecting life.The influences of aspect ratio of an initial crack and angle of crack relative to the center axis of the pressure cylinder on life will be analyzed in this section.The acceptable crack size is still 1/10 of the cylinder wall thickness.Figs.13-14 show the crack growth curves under different initial aspect ratios.Fig.10 presents the analysis of the crack growth at different initial depths under the same initial crack half-lengthc0of 5 mm while Fig.13 shows the crack growth with different initial lengths at the same initial crack deptha0of 3 mm.When the initial crack length was 5 mm,the crack growth at initial depths of 2 mm,3 mm,4 mm and 5 mm was investigated.The minimum residual life was 96 668 cycles.When the initial crack depth was fixed to 3 mm,the crack growth was investigated when the initial length was 3 mm,4 mm,5 mm and 6 mm respectively.The minimum residual life was 102 566 cycles.The remaining life can meet the application requirement.

    Fig.13 Effects of initial aspect ratio on crack growth in depth direction when c0 is set to 5 mm

    The effect of the crack inclination angle to the cylinder axis was also studied under the condition ofa0=3 mm andc0=5 mm.The configuration of the surface crack inclination in the FE model is shown in Fig.15.The inclination angles of cracks are 0°,5°,10°,20°and 30°respectively.Results in Fig.16 show that the angle of inclination retards crack growth.Then conservative calculation can be done by not taking inclination into account.

    Fig.15 Configuration of surface crack inclination in FE model

    Fig.16 Effects of inclination on crack growth in depth direction when a0=3 mm and c0=5 mm

    3 Conclusions

    UHP chambers used to simulate deep-sea environment need to bear a pressure of more than 110 MPa.After several years of use,the UHP chambers need to be checked in a comprehensive way.The residual life of a certain aged UHP chamber under the condition of cracking was studied.One tenth of the chamber’s wall thickness was defined as the allowable depth of surface cracks on the inner wall.Based on this criterion,the residual life of the chamber was analyzed under the assumption that there were macroscopic inner surface cracks.

    The UHP chambers sometimes withstand trapezoidal load with holding time at peak stress in service.Therefore,in this paper,the sensitivity of material to topographic load was studied experimentally at first.The sustaining time of 2 min at peak stress was set for each cycle.The experimental results were compared with those under triangular fatigue load,which proves that the material is not sensitive to holding time.The material parameters obtained from conventional fatigue crack growth test can be used to evaluate the crack growth in UHP chambers.

    Franc3D was used to calculate the crack growth,and its accuracy was verified by comparing the computed results of CT test with the experimental results.The crack growth with different initial sizes and depth to length ratios is analyzed,and the effect of the crack inclination angle to the cylinder axis was studied.Results show that the angle of inclination retards crack growth,and that,under the current criterion,the crack depth is less than one tenth of the wall thickness,and the life requirement is satisfied.According to the analysis,the service life of the pressure cylinder in the design pressure state is about six times the safety margin.The analysis results can be used as a reference for evaluation of the reliability of UHP chambers.

    亚洲国产成人一精品久久久| 国产精品.久久久| 美女脱内裤让男人舔精品视频| 亚洲天堂av无毛| 自拍欧美九色日韩亚洲蝌蚪91 | 久久97久久精品| 在现免费观看毛片| 少妇人妻久久综合中文| 欧美日韩国产mv在线观看视频 | 一本色道久久久久久精品综合| 色视频在线一区二区三区| 男人爽女人下面视频在线观看| 熟妇人妻不卡中文字幕| 岛国毛片在线播放| 日本午夜av视频| 久久久久人妻精品一区果冻| 乱系列少妇在线播放| videos熟女内射| 成人免费观看视频高清| 欧美日韩视频精品一区| 日本爱情动作片www.在线观看| 51国产日韩欧美| 国产 一区 欧美 日韩| 国产欧美亚洲国产| 最黄视频免费看| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品国产色婷婷电影| 日本欧美国产在线视频| 国产日韩欧美在线精品| 午夜福利在线观看免费完整高清在| 日本wwww免费看| 欧美日韩视频高清一区二区三区二| 久久精品夜色国产| 国产精品国产av在线观看| 我的老师免费观看完整版| 最近最新中文字幕免费大全7| 日日摸夜夜添夜夜爱| 欧美3d第一页| 欧美激情极品国产一区二区三区 | 高清在线视频一区二区三区| 国产精品久久久久久久电影| 日韩中文字幕视频在线看片 | 男人舔奶头视频| 国产男女内射视频| 秋霞在线观看毛片| 国产午夜精品一二区理论片| 自拍偷自拍亚洲精品老妇| a级毛片免费高清观看在线播放| 国产片特级美女逼逼视频| 在线观看一区二区三区激情| 午夜福利在线观看免费完整高清在| 久久精品久久久久久噜噜老黄| 欧美精品国产亚洲| 国产亚洲欧美精品永久| 色5月婷婷丁香| 欧美日韩综合久久久久久| 欧美精品一区二区大全| 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 日韩欧美精品免费久久| 男人爽女人下面视频在线观看| 国产一级毛片在线| 99热这里只有是精品50| 国产精品不卡视频一区二区| 少妇高潮的动态图| 男女免费视频国产| 免费高清在线观看视频在线观看| 欧美日韩在线观看h| 黄片无遮挡物在线观看| 3wmmmm亚洲av在线观看| 丰满人妻一区二区三区视频av| 亚洲精品国产色婷婷电影| 国产精品久久久久久av不卡| 亚洲国产精品成人久久小说| 99热这里只有是精品在线观看| 精品亚洲成a人片在线观看 | 日韩制服骚丝袜av| av国产久精品久网站免费入址| 老熟女久久久| 波野结衣二区三区在线| 国产亚洲91精品色在线| 亚洲欧美一区二区三区国产| 亚洲性久久影院| 欧美精品一区二区免费开放| 好男人视频免费观看在线| 99久久人妻综合| 老司机影院毛片| 高清日韩中文字幕在线| 国产免费视频播放在线视频| 成年免费大片在线观看| 另类亚洲欧美激情| 亚洲av国产av综合av卡| 一级片'在线观看视频| 18禁裸乳无遮挡免费网站照片| 色视频www国产| 大又大粗又爽又黄少妇毛片口| 成人免费观看视频高清| 国产真实伦视频高清在线观看| 激情五月婷婷亚洲| 日本爱情动作片www.在线观看| 国产中年淑女户外野战色| 久久毛片免费看一区二区三区| 2022亚洲国产成人精品| 久久青草综合色| 亚洲av不卡在线观看| 少妇高潮的动态图| 国产成人精品婷婷| 亚洲成人一二三区av| 看免费成人av毛片| 亚洲精品乱久久久久久| 狂野欧美激情性bbbbbb| 久久久久国产网址| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 久久久欧美国产精品| 国产精品偷伦视频观看了| 国产高清有码在线观看视频| 又黄又爽又刺激的免费视频.| 在线看a的网站| 有码 亚洲区| 99热全是精品| 黄色日韩在线| 国产男人的电影天堂91| 亚洲欧美日韩另类电影网站 | 又粗又硬又长又爽又黄的视频| 亚洲人与动物交配视频| 看十八女毛片水多多多| 国产探花极品一区二区| 91aial.com中文字幕在线观看| 久久99热这里只频精品6学生| 一级av片app| 国产精品.久久久| 只有这里有精品99| 一级二级三级毛片免费看| 亚洲国产成人一精品久久久| 人体艺术视频欧美日本| 免费大片18禁| 一级爰片在线观看| 午夜日本视频在线| 毛片一级片免费看久久久久| 日韩欧美一区视频在线观看 | 欧美性感艳星| 久久精品国产亚洲av天美| 日本免费在线观看一区| 国产高清国产精品国产三级 | 精品午夜福利在线看| 欧美xxxx性猛交bbbb| 亚洲精品一二三| 日本色播在线视频| 国产精品熟女久久久久浪| 国产精品国产av在线观看| av免费观看日本| 最后的刺客免费高清国语| 日韩视频在线欧美| 成人黄色视频免费在线看| 好男人视频免费观看在线| 精品一区二区三卡| 日本av免费视频播放| 嫩草影院入口| 亚洲精华国产精华液的使用体验| 最黄视频免费看| 我要看黄色一级片免费的| 欧美变态另类bdsm刘玥| 最黄视频免费看| 国产亚洲精品久久久com| 国产伦精品一区二区三区四那| 日本av免费视频播放| 欧美日韩国产mv在线观看视频 | 亚洲欧美日韩卡通动漫| 秋霞伦理黄片| 哪个播放器可以免费观看大片| 国模一区二区三区四区视频| 久久久精品免费免费高清| 观看av在线不卡| 国语对白做爰xxxⅹ性视频网站| av专区在线播放| 成人黄色视频免费在线看| 一本久久精品| 成人高潮视频无遮挡免费网站| 国产精品国产三级国产av玫瑰| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| 亚洲精品国产成人久久av| 在线亚洲精品国产二区图片欧美 | 国产高清不卡午夜福利| 丰满乱子伦码专区| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品亚洲一区二区| 91精品国产国语对白视频| 一个人免费看片子| 观看av在线不卡| 热99国产精品久久久久久7| 中国美白少妇内射xxxbb| 草草在线视频免费看| 成年人午夜在线观看视频| 国产成人一区二区在线| www.色视频.com| 午夜激情福利司机影院| 久久久久国产网址| 久久影院123| 新久久久久国产一级毛片| 亚洲激情五月婷婷啪啪| 高清欧美精品videossex| 熟女av电影| 久久99蜜桃精品久久| 国产亚洲最大av| 久久精品国产亚洲av天美| 身体一侧抽搐| 色视频在线一区二区三区| 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 深夜a级毛片| 成人18禁高潮啪啪吃奶动态图 | videossex国产| a级一级毛片免费在线观看| 少妇被粗大猛烈的视频| 少妇 在线观看| 波野结衣二区三区在线| 综合色丁香网| 人妻 亚洲 视频| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 人人妻人人澡人人爽人人夜夜| 成人特级av手机在线观看| 成年av动漫网址| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 欧美丝袜亚洲另类| 大话2 男鬼变身卡| 美女内射精品一级片tv| 日韩在线高清观看一区二区三区| 中文精品一卡2卡3卡4更新| 欧美成人精品欧美一级黄| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 人妻 亚洲 视频| 精品一区二区免费观看| 久热久热在线精品观看| 边亲边吃奶的免费视频| 最近手机中文字幕大全| 亚洲色图av天堂| 黑人猛操日本美女一级片| 久久精品国产鲁丝片午夜精品| 蜜臀久久99精品久久宅男| 亚洲va在线va天堂va国产| 亚洲无线观看免费| 免费久久久久久久精品成人欧美视频 | 晚上一个人看的免费电影| 纵有疾风起免费观看全集完整版| 亚洲va在线va天堂va国产| 麻豆成人av视频| 大片电影免费在线观看免费| 精品亚洲成a人片在线观看 | 精品一品国产午夜福利视频| 最近最新中文字幕免费大全7| 国产精品伦人一区二区| 亚洲真实伦在线观看| 国内精品宾馆在线| 大陆偷拍与自拍| 黄色怎么调成土黄色| 一级爰片在线观看| 七月丁香在线播放| 特大巨黑吊av在线直播| 亚洲欧美日韩无卡精品| 日韩 亚洲 欧美在线| 精品国产露脸久久av麻豆| 观看美女的网站| 久久婷婷青草| 99热网站在线观看| 99热这里只有精品一区| 国产 一区精品| 人妻少妇偷人精品九色| 成年美女黄网站色视频大全免费 | 在线观看国产h片| 观看免费一级毛片| 日本vs欧美在线观看视频 | 亚洲欧美成人精品一区二区| 日韩av免费高清视频| 丝瓜视频免费看黄片| 免费高清在线观看视频在线观看| 欧美xxxx性猛交bbbb| 肉色欧美久久久久久久蜜桃| 一区二区三区免费毛片| 网址你懂的国产日韩在线| 欧美一区二区亚洲| 国产午夜精品一二区理论片| 多毛熟女@视频| tube8黄色片| 美女主播在线视频| 嘟嘟电影网在线观看| 国产乱人视频| a 毛片基地| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 黄色一级大片看看| 嫩草影院新地址| a级一级毛片免费在线观看| 18禁裸乳无遮挡动漫免费视频| 欧美成人精品欧美一级黄| 婷婷色麻豆天堂久久| 2021少妇久久久久久久久久久| 网址你懂的国产日韩在线| 黄色怎么调成土黄色| 午夜免费男女啪啪视频观看| 五月玫瑰六月丁香| 国产一区有黄有色的免费视频| 最新中文字幕久久久久| 18禁动态无遮挡网站| 只有这里有精品99| 久久99热6这里只有精品| 边亲边吃奶的免费视频| 亚洲国产色片| 寂寞人妻少妇视频99o| 99久久综合免费| 人妻制服诱惑在线中文字幕| 18禁在线播放成人免费| 国产毛片在线视频| 久久影院123| 人妻 亚洲 视频| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产高清在线一区二区三| 男男h啪啪无遮挡| 精品少妇黑人巨大在线播放| 欧美性感艳星| 欧美人与善性xxx| 激情 狠狠 欧美| 人妻 亚洲 视频| 国产成人精品婷婷| 免费人成在线观看视频色| 黑人高潮一二区| av黄色大香蕉| 国产精品一及| av在线老鸭窝| 中文精品一卡2卡3卡4更新| 国产大屁股一区二区在线视频| 黄片wwwwww| 嘟嘟电影网在线观看| 国产午夜精品久久久久久一区二区三区| 熟女人妻精品中文字幕| 久久久久久九九精品二区国产| 五月天丁香电影| 免费观看在线日韩| 国产精品嫩草影院av在线观看| 夫妻性生交免费视频一级片| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲av涩爱| 十分钟在线观看高清视频www | 直男gayav资源| 国产一区二区在线观看日韩| 在线亚洲精品国产二区图片欧美 | 在线观看免费日韩欧美大片 | 下体分泌物呈黄色| 男女啪啪激烈高潮av片| 一级黄片播放器| 久久精品人妻少妇| 免费在线观看成人毛片| 国产精品成人在线| 干丝袜人妻中文字幕| 国产黄频视频在线观看| 国产v大片淫在线免费观看| 日韩av在线免费看完整版不卡| 好男人视频免费观看在线| 国产毛片在线视频| 18禁裸乳无遮挡动漫免费视频| 干丝袜人妻中文字幕| 日本-黄色视频高清免费观看| 国产综合精华液| 精品一区二区三卡| 欧美日韩视频精品一区| 精品一区二区三区视频在线| 在线精品无人区一区二区三 | 又爽又黄a免费视频| 男女边摸边吃奶| 少妇猛男粗大的猛烈进出视频| 欧美97在线视频| 精品人妻视频免费看| 91久久精品国产一区二区成人| 日韩欧美一区视频在线观看 | 久久久久久伊人网av| 在线免费观看不下载黄p国产| 国产成人精品福利久久| 亚洲综合色惰| 国产乱来视频区| 国产av精品麻豆| 超碰av人人做人人爽久久| 成人一区二区视频在线观看| 亚洲精品国产成人久久av| 97超碰精品成人国产| 欧美性感艳星| 黄色一级大片看看| 最后的刺客免费高清国语| 国产精品久久久久久久电影| 成年人午夜在线观看视频| av网站免费在线观看视频| 亚洲欧美一区二区三区国产| 国产精品成人在线| 午夜精品国产一区二区电影| 国产成人91sexporn| 亚洲图色成人| 最近最新中文字幕免费大全7| 亚洲国产毛片av蜜桃av| 国产精品三级大全| 国产毛片在线视频| 国产免费视频播放在线视频| 夜夜看夜夜爽夜夜摸| 大香蕉97超碰在线| 日韩欧美 国产精品| 狂野欧美白嫩少妇大欣赏| 午夜福利高清视频| 国产69精品久久久久777片| 国产黄频视频在线观看| 久久久久久久国产电影| 国产精品一及| 秋霞伦理黄片| 天堂8中文在线网| 久久久久国产网址| 国语对白做爰xxxⅹ性视频网站| 永久网站在线| 在线免费观看不下载黄p国产| 国产黄频视频在线观看| 熟妇人妻不卡中文字幕| 一本色道久久久久久精品综合| 欧美97在线视频| 丰满乱子伦码专区| 久久精品国产亚洲av涩爱| 久久韩国三级中文字幕| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久精品电影小说 | 久久 成人 亚洲| 久久久久久久久久人人人人人人| 五月天丁香电影| 亚洲一区二区三区欧美精品| 欧美 日韩 精品 国产| 亚洲欧洲国产日韩| 免费少妇av软件| 男女边吃奶边做爰视频| 久久av网站| 免费观看性生交大片5| 女人久久www免费人成看片| 能在线免费看毛片的网站| 亚洲av日韩在线播放| 51国产日韩欧美| 91aial.com中文字幕在线观看| 中文字幕免费在线视频6| 亚洲国产精品专区欧美| 免费人成在线观看视频色| 你懂的网址亚洲精品在线观看| 国产av精品麻豆| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 中文天堂在线官网| 午夜老司机福利剧场| 伊人久久精品亚洲午夜| 美女xxoo啪啪120秒动态图| 少妇的逼水好多| 各种免费的搞黄视频| 亚洲美女搞黄在线观看| 国产午夜精品久久久久久一区二区三区| 国产精品嫩草影院av在线观看| 看免费成人av毛片| 狂野欧美激情性xxxx在线观看| 亚洲第一av免费看| 日本与韩国留学比较| 亚洲国产色片| 国产男女超爽视频在线观看| 国产精品精品国产色婷婷| 成人美女网站在线观看视频| 26uuu在线亚洲综合色| 久久久久性生活片| 成人毛片60女人毛片免费| 性色avwww在线观看| 看十八女毛片水多多多| 综合色丁香网| 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| videos熟女内射| 国产亚洲午夜精品一区二区久久| 伦精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 日韩欧美 国产精品| 干丝袜人妻中文字幕| 国产成人91sexporn| 成年女人在线观看亚洲视频| 最近中文字幕高清免费大全6| 老司机影院毛片| 国产 精品1| 日本午夜av视频| 少妇精品久久久久久久| 成人18禁高潮啪啪吃奶动态图 | 97精品久久久久久久久久精品| 观看av在线不卡| 久久久久久久亚洲中文字幕| 久久综合国产亚洲精品| 国产有黄有色有爽视频| 秋霞在线观看毛片| 国产伦精品一区二区三区视频9| 大又大粗又爽又黄少妇毛片口| 一级二级三级毛片免费看| 老司机影院毛片| 最近2019中文字幕mv第一页| 视频区图区小说| 大又大粗又爽又黄少妇毛片口| 综合色丁香网| 久久人妻熟女aⅴ| 国产精品一二三区在线看| 99久久精品热视频| 欧美变态另类bdsm刘玥| 99久久综合免费| 久久热精品热| 内地一区二区视频在线| 久久久成人免费电影| 午夜视频国产福利| 嫩草影院入口| 男女免费视频国产| 乱系列少妇在线播放| 中文字幕av成人在线电影| 中文字幕制服av| xxx大片免费视频| 青春草视频在线免费观看| 国产成人精品婷婷| 亚洲,欧美,日韩| 高清欧美精品videossex| 日本与韩国留学比较| 国产大屁股一区二区在线视频| 在线观看国产h片| 欧美日韩综合久久久久久| 午夜激情福利司机影院| 国产成人91sexporn| 久久精品熟女亚洲av麻豆精品| 亚洲内射少妇av| 亚洲人成网站在线观看播放| 永久免费av网站大全| 国产高潮美女av| 国产 精品1| 在线观看免费日韩欧美大片 | 久久久久精品久久久久真实原创| 男人狂女人下面高潮的视频| 观看av在线不卡| 精品亚洲乱码少妇综合久久| 国产欧美另类精品又又久久亚洲欧美| 在线观看免费视频网站a站| 国产精品欧美亚洲77777| av视频免费观看在线观看| 九色成人免费人妻av| 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 国产深夜福利视频在线观看| 久久99热这里只有精品18| 精品亚洲成a人片在线观看 | 又黄又爽又刺激的免费视频.| 久久久久视频综合| 免费少妇av软件| 免费不卡的大黄色大毛片视频在线观看| 一个人免费看片子| 一边亲一边摸免费视频| 赤兔流量卡办理| 欧美一区二区亚洲| 精品酒店卫生间| 精品久久久精品久久久| 久久人妻熟女aⅴ| 国产精品一区www在线观看| 亚洲精品自拍成人| 男人狂女人下面高潮的视频| 成人国产麻豆网| 最近最新中文字幕大全电影3| 国产亚洲欧美精品永久| 国产精品人妻久久久久久| 午夜免费鲁丝| 国产精品.久久久| 国产亚洲精品久久久com| 日韩av不卡免费在线播放| 成年免费大片在线观看| 最新中文字幕久久久久| 丝袜喷水一区| 亚洲国产欧美人成| 尾随美女入室| 免费观看无遮挡的男女| 国产有黄有色有爽视频| 简卡轻食公司| 黄片无遮挡物在线观看| 精品一区二区三卡| 99久久综合免费| 国产精品福利在线免费观看| 欧美97在线视频| 午夜福利视频精品| 久久99精品国语久久久| 亚洲第一av免费看| 少妇猛男粗大的猛烈进出视频| 婷婷色麻豆天堂久久| 国产一区有黄有色的免费视频| 国精品久久久久久国模美| 日韩人妻高清精品专区| 国产亚洲av片在线观看秒播厂| 欧美一区二区亚洲| 国产女主播在线喷水免费视频网站| 欧美成人午夜免费资源| 欧美一区二区亚洲| 日韩欧美 国产精品| 精品亚洲乱码少妇综合久久| 视频中文字幕在线观看| 久久热精品热| tube8黄色片| 国产免费视频播放在线视频| 国产精品国产三级专区第一集| 日本欧美视频一区| 卡戴珊不雅视频在线播放| 亚洲精品自拍成人| 一级爰片在线观看| 男男h啪啪无遮挡| 国产午夜精品久久久久久一区二区三区| 联通29元200g的流量卡| 亚洲国产最新在线播放| 少妇人妻一区二区三区视频|