• <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| 观看美女的网站| 卡戴珊不雅视频在线播放| 亚洲va在线va天堂va国产| 久久久久网色| 特大巨黑吊av在线直播| 老女人水多毛片| 99热精品在线国产| 少妇熟女aⅴ在线视频| av福利片在线观看| 日本黄色视频三级网站网址| 亚洲欧美日韩无卡精品| 日韩制服骚丝袜av| 淫秽高清视频在线观看| 亚洲精品国产av成人精品| 免费看美女性在线毛片视频| 男插女下体视频免费在线播放| 成人av在线播放网站| 国产黄色视频一区二区在线观看 | 亚洲欧美日韩高清专用| 久久精品国产99精品国产亚洲性色| 内地一区二区视频在线| 久久久久久久久中文| 国产精品国产高清国产av| 精品国产一区二区三区久久久樱花 | ponron亚洲| 汤姆久久久久久久影院中文字幕 | 国产综合懂色| 国产成人午夜福利电影在线观看| 国产精品.久久久| 国产极品精品免费视频能看的| 婷婷色麻豆天堂久久 | a级毛色黄片| 国产精品一二三区在线看| 欧美极品一区二区三区四区| 可以在线观看毛片的网站| 色5月婷婷丁香| 成人欧美大片| 国产黄片美女视频| 三级经典国产精品| 国语自产精品视频在线第100页| 中国美白少妇内射xxxbb| 白带黄色成豆腐渣| 亚洲精品乱码久久久久久按摩| 日韩精品有码人妻一区| 国产高清国产精品国产三级 | 久久久久久九九精品二区国产| 免费看a级黄色片| 久久亚洲精品不卡| 中国国产av一级| 能在线免费看毛片的网站| 国产精品美女特级片免费视频播放器| 亚洲图色成人| 看片在线看免费视频| 久久久久久久久久成人| 国产成人精品久久久久久| 久久99精品国语久久久| 亚洲精品乱久久久久久| 一夜夜www| 亚洲综合精品二区| 国产一级毛片在线| 日韩一区二区视频免费看| 欧美激情久久久久久爽电影| 国产免费一级a男人的天堂| 久久久亚洲精品成人影院| 亚洲美女视频黄频| 青青草视频在线视频观看| 久久国内精品自在自线图片| 99久久人妻综合| 国产又色又爽无遮挡免| 国产高清不卡午夜福利| 纵有疾风起免费观看全集完整版 | 亚洲av男天堂| 观看美女的网站| 亚洲最大成人手机在线| 亚洲国产精品久久男人天堂| av在线老鸭窝| 日本熟妇午夜| 国产爱豆传媒在线观看| eeuss影院久久| 精品午夜福利在线看| 精品免费久久久久久久清纯| 人妻系列 视频| 国产在线一区二区三区精 | 内射极品少妇av片p| 在线观看66精品国产| 婷婷色麻豆天堂久久 | 国产精品久久久久久久电影| av播播在线观看一区| 亚洲av电影在线观看一区二区三区 | 国产精品人妻久久久影院| 中文乱码字字幕精品一区二区三区 | 亚洲欧美成人精品一区二区| 日本午夜av视频| 亚洲av熟女| 国产精品久久电影中文字幕| 日本-黄色视频高清免费观看| 五月玫瑰六月丁香| 最近的中文字幕免费完整| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 日本av手机在线免费观看| 国产精品蜜桃在线观看| 全区人妻精品视频| 久久久久久大精品| 免费看日本二区| 岛国毛片在线播放| 国产欧美另类精品又又久久亚洲欧美| 建设人人有责人人尽责人人享有的 | 汤姆久久久久久久影院中文字幕 | 春色校园在线视频观看| 啦啦啦韩国在线观看视频| 91午夜精品亚洲一区二区三区| 舔av片在线| 婷婷色av中文字幕| 免费观看在线日韩| 欧美性猛交╳xxx乱大交人| 欧美3d第一页| 夫妻性生交免费视频一级片| 成人特级av手机在线观看| 女人被狂操c到高潮| 国内精品一区二区在线观看| 国产欧美另类精品又又久久亚洲欧美| 日本黄色片子视频| 波野结衣二区三区在线| 国产一区亚洲一区在线观看| 国产69精品久久久久777片| 伊人久久精品亚洲午夜| 爱豆传媒免费全集在线观看| 精品无人区乱码1区二区| 可以在线观看毛片的网站| 看十八女毛片水多多多| 永久免费av网站大全| 三级国产精品片| 国产精品麻豆人妻色哟哟久久 | 国产淫语在线视频| 三级男女做爰猛烈吃奶摸视频| 国产一级毛片七仙女欲春2| 哪个播放器可以免费观看大片| 亚洲自偷自拍三级| 一级毛片电影观看 | 成人美女网站在线观看视频| 亚洲中文字幕日韩| 日韩av不卡免费在线播放| 国国产精品蜜臀av免费| 久久精品久久精品一区二区三区| 蜜桃亚洲精品一区二区三区| 99久久人妻综合| 久久久久国产网址| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 免费一级毛片在线播放高清视频| 日本免费在线观看一区| 亚洲精品日韩av片在线观看| 亚洲18禁久久av| 亚洲五月天丁香| 欧美成人精品欧美一级黄| 午夜福利成人在线免费观看| 九九在线视频观看精品| 婷婷色综合大香蕉| 国产一区二区亚洲精品在线观看| av卡一久久| 中文在线观看免费www的网站| 搞女人的毛片| 欧美精品一区二区大全| 97超视频在线观看视频| 国产亚洲91精品色在线| av国产久精品久网站免费入址| 99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 亚洲av二区三区四区| 国产乱人视频| 日本猛色少妇xxxxx猛交久久| 男人的好看免费观看在线视频| 高清毛片免费看| 特级一级黄色大片| 岛国毛片在线播放| 日本熟妇午夜| 亚洲精品,欧美精品| 亚洲av熟女| 国产黄色视频一区二区在线观看 | 久久久亚洲精品成人影院| av播播在线观看一区| 99在线人妻在线中文字幕| 校园人妻丝袜中文字幕| 能在线免费看毛片的网站| 又爽又黄无遮挡网站| 日韩一区二区视频免费看| 人体艺术视频欧美日本| 日本与韩国留学比较| 亚洲av成人av| 日本一二三区视频观看| 国产精品一二三区在线看| 国产午夜精品一二区理论片| 国产白丝娇喘喷水9色精品| 国产一区二区三区av在线| 91久久精品国产一区二区三区| 99久国产av精品国产电影| 高清在线视频一区二区三区 | 久久精品影院6| 男女视频在线观看网站免费| 国产综合懂色| av国产免费在线观看| 最后的刺客免费高清国语| 国产精品野战在线观看| 久久99热这里只有精品18| 久久精品综合一区二区三区| 美女cb高潮喷水在线观看| 亚洲成人久久爱视频| 成人av在线播放网站| 精品久久久久久久人妻蜜臀av| 只有这里有精品99| av视频在线观看入口| 国产精品久久电影中文字幕| 伊人久久精品亚洲午夜| 免费播放大片免费观看视频在线观看 | 欧美zozozo另类| 国产av一区在线观看免费| 午夜福利网站1000一区二区三区| 色哟哟·www| 人妻系列 视频| 国产精品1区2区在线观看.| 我的女老师完整版在线观看| 男人和女人高潮做爰伦理| 一级毛片我不卡| 久久精品熟女亚洲av麻豆精品 | 成人高潮视频无遮挡免费网站| 人人妻人人澡欧美一区二区| 久久久久久久午夜电影| 免费看av在线观看网站| 国产成人freesex在线| www.色视频.com| 国内精品宾馆在线| 精品久久久久久久人妻蜜臀av| 久久99精品国语久久久| 最近中文字幕高清免费大全6| 欧美成人免费av一区二区三区| 国产毛片a区久久久久| 免费观看a级毛片全部| 天美传媒精品一区二区| 久热久热在线精品观看| 久久精品国产亚洲av涩爱| 91av网一区二区| 91在线精品国自产拍蜜月| 国产中年淑女户外野战色| 看片在线看免费视频| 日本五十路高清| 午夜亚洲福利在线播放| 婷婷色av中文字幕| 九草在线视频观看| 久久99热这里只有精品18| 一边摸一边抽搐一进一小说| 尤物成人国产欧美一区二区三区| 美女脱内裤让男人舔精品视频| 亚州av有码| 国产精品嫩草影院av在线观看| 成人美女网站在线观看视频| 一区二区三区高清视频在线| 欧美性猛交黑人性爽| 国产黄色视频一区二区在线观看 | 免费搜索国产男女视频| 国产精品麻豆人妻色哟哟久久 | 日日撸夜夜添| 人妻少妇偷人精品九色| 麻豆成人av视频| 亚洲av.av天堂| 免费黄色在线免费观看| 熟妇人妻久久中文字幕3abv| 日韩中字成人| 国产精品久久视频播放| 国产精品人妻久久久久久| 亚洲美女视频黄频| 边亲边吃奶的免费视频| 狂野欧美激情性xxxx在线观看| 久久精品国产99精品国产亚洲性色| 亚洲成人av在线免费| 淫秽高清视频在线观看| 久久精品国产亚洲av天美| 日韩欧美 国产精品| 中国国产av一级| 日韩欧美国产在线观看| 亚洲真实伦在线观看| 久久精品久久精品一区二区三区| 久久人妻av系列| 日韩欧美 国产精品| 一个人免费在线观看电影| 内地一区二区视频在线| 黄色配什么色好看| 国产真实伦视频高清在线观看| 免费观看性生交大片5| 午夜福利成人在线免费观看| 高清午夜精品一区二区三区| 长腿黑丝高跟| 非洲黑人性xxxx精品又粗又长| 国内精品一区二区在线观看| 黑人高潮一二区| 亚洲人与动物交配视频| 不卡视频在线观看欧美| 亚洲国产精品专区欧美| 成人性生交大片免费视频hd| 特级一级黄色大片| 我的女老师完整版在线观看| 亚洲人与动物交配视频| 国产毛片a区久久久久| 色哟哟·www| 亚洲在线自拍视频| 亚洲一级一片aⅴ在线观看| 99热这里只有精品一区| 成人国产麻豆网| 女人被狂操c到高潮| 日日啪夜夜撸| 国产视频内射| 亚洲乱码一区二区免费版| 简卡轻食公司| 日韩一本色道免费dvd| 国产精华一区二区三区| 永久网站在线| 2021少妇久久久久久久久久久| 黄色日韩在线| av卡一久久| 岛国毛片在线播放| 国产精品久久久久久av不卡| 白带黄色成豆腐渣| av卡一久久| 色综合色国产| 亚洲国产精品成人综合色| 伦精品一区二区三区| 噜噜噜噜噜久久久久久91| 99热精品在线国产| 久久久久久国产a免费观看| 久久久精品大字幕| 亚洲av免费在线观看| 午夜久久久久精精品| 亚洲怡红院男人天堂| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 99热全是精品| 蜜臀久久99精品久久宅男| 非洲黑人性xxxx精品又粗又长| 日本免费一区二区三区高清不卡| 久久精品夜色国产| 男人舔奶头视频| 国产真实乱freesex| 亚洲第一区二区三区不卡| 一本久久精品| 身体一侧抽搐| 欧美色视频一区免费| 国产精品一区二区性色av| 最近手机中文字幕大全| 日本黄大片高清| 国产精品熟女久久久久浪| 精品国产露脸久久av麻豆 | 99在线视频只有这里精品首页| 天堂中文最新版在线下载 | 国产国拍精品亚洲av在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲精品国产成人久久av| 精品久久久久久久人妻蜜臀av| 有码 亚洲区| 亚洲三级黄色毛片| 国产成人福利小说| 午夜视频国产福利| 日韩强制内射视频| 特级一级黄色大片| 亚洲国产精品成人综合色| 99久久九九国产精品国产免费| 日韩av在线大香蕉| 精品不卡国产一区二区三区| 亚洲精品国产av成人精品| 日本免费在线观看一区| 亚洲国产色片| 91午夜精品亚洲一区二区三区| 欧美xxxx性猛交bbbb| 中文字幕av成人在线电影| 热99re8久久精品国产| 2022亚洲国产成人精品| 欧美丝袜亚洲另类| 国产亚洲5aaaaa淫片| 免费在线观看成人毛片| 亚洲欧美精品专区久久| 欧美极品一区二区三区四区| 亚洲综合色惰| 中文字幕久久专区| 久久久久久国产a免费观看| 赤兔流量卡办理| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 欧美性感艳星| 伦理电影大哥的女人| 美女内射精品一级片tv| 免费观看在线日韩| 亚洲欧美一区二区三区国产| 久久99蜜桃精品久久| 嘟嘟电影网在线观看| 色噜噜av男人的天堂激情| 亚洲精品乱码久久久v下载方式| 色综合亚洲欧美另类图片| 在线播放国产精品三级| 爱豆传媒免费全集在线观看| 看片在线看免费视频| 激情 狠狠 欧美| 欧美+日韩+精品| 亚洲欧美日韩卡通动漫| 人体艺术视频欧美日本| 亚洲国产精品成人久久小说| 国产av在哪里看| 91精品国产九色| 久久亚洲精品不卡| 国产成人a区在线观看| 国产成人福利小说| 国产伦在线观看视频一区| 国产高清有码在线观看视频| 久久久久久久久久久免费av| 插阴视频在线观看视频| 纵有疾风起免费观看全集完整版 | 成年av动漫网址| 美女国产视频在线观看| 成人一区二区视频在线观看| 国产高清国产精品国产三级 | 禁无遮挡网站| 1000部很黄的大片| 桃色一区二区三区在线观看| 久久精品夜夜夜夜夜久久蜜豆| av线在线观看网站| 五月玫瑰六月丁香| 欧美高清成人免费视频www| 午夜亚洲福利在线播放| 国产高清不卡午夜福利| 在线免费十八禁| 国产精品国产三级专区第一集| 亚洲精品456在线播放app| 久久精品久久久久久噜噜老黄 | or卡值多少钱| 伊人久久精品亚洲午夜| 亚洲va在线va天堂va国产| 亚洲第一区二区三区不卡| 国产美女午夜福利| 国产极品天堂在线| 特大巨黑吊av在线直播| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久人妻蜜臀av| 成人漫画全彩无遮挡| 日本熟妇午夜| 日本一二三区视频观看| 日韩人妻高清精品专区| 亚洲人与动物交配视频| 五月伊人婷婷丁香| 国产一区有黄有色的免费视频 | 国产高清有码在线观看视频| 高清午夜精品一区二区三区| 热99re8久久精品国产| 国产 一区精品| 免费黄网站久久成人精品| 我的女老师完整版在线观看| 赤兔流量卡办理| 69人妻影院| 欧美最新免费一区二区三区| 日本一二三区视频观看| 2022亚洲国产成人精品| 亚洲精品乱码久久久v下载方式| 国产精品人妻久久久影院| 天堂av国产一区二区熟女人妻| 99九九线精品视频在线观看视频| 国产精品1区2区在线观看.| 好男人在线观看高清免费视频| 蜜臀久久99精品久久宅男| 成人亚洲精品av一区二区| 久久久久久伊人网av| 国产精品野战在线观看| 亚洲欧美日韩卡通动漫| 国内精品美女久久久久久| 99热这里只有是精品在线观看| 免费播放大片免费观看视频在线观看 | 日本黄色视频三级网站网址| 秋霞在线观看毛片| 久久精品久久久久久噜噜老黄 | 我要搜黄色片| 人体艺术视频欧美日本| 热99re8久久精品国产| 色噜噜av男人的天堂激情| 国产大屁股一区二区在线视频| 少妇的逼好多水| 国产午夜精品论理片| 男女视频在线观看网站免费| 人妻夜夜爽99麻豆av| 免费看av在线观看网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆成人午夜福利视频| 日本色播在线视频| 久久久精品大字幕| 视频中文字幕在线观看| 欧美三级亚洲精品| 欧美3d第一页| 97超视频在线观看视频| 国产成人一区二区在线| 91在线精品国自产拍蜜月| 天堂中文最新版在线下载 | 两个人视频免费观看高清| 久久精品91蜜桃| 久久99热6这里只有精品| 国产精品久久久久久精品电影| 麻豆一二三区av精品| 国产免费视频播放在线视频 | 国产一区二区在线av高清观看| 伊人久久精品亚洲午夜| 男人和女人高潮做爰伦理| 亚洲成人精品中文字幕电影| 国产高清三级在线| 免费观看人在逋| 你懂的网址亚洲精品在线观看 | av在线观看视频网站免费| 国产精品久久电影中文字幕| 色哟哟·www| 国产精品国产高清国产av| 久久久国产成人免费| 日本午夜av视频| 黄色欧美视频在线观看| 免费观看性生交大片5| 国产免费福利视频在线观看| 韩国高清视频一区二区三区| 色综合亚洲欧美另类图片| 人人妻人人看人人澡| 99久久无色码亚洲精品果冻| 亚洲av日韩在线播放| 91精品国产九色| 国产亚洲91精品色在线| 黄色配什么色好看| 久99久视频精品免费| 久久精品夜夜夜夜夜久久蜜豆| 国产欧美日韩精品一区二区| 亚洲av熟女| 中文天堂在线官网| 亚洲美女搞黄在线观看| 有码 亚洲区| h日本视频在线播放| 乱系列少妇在线播放| 精品国产三级普通话版| 中文字幕av在线有码专区| 国产精品久久久久久av不卡| 午夜日本视频在线| 国产在线一区二区三区精 | 欧美潮喷喷水| 亚洲av熟女| 午夜福利高清视频| 99久久人妻综合| 舔av片在线| 你懂的网址亚洲精品在线观看 | 在线免费观看的www视频| 91精品一卡2卡3卡4卡| 亚洲国产最新在线播放| 午夜福利在线在线| 中文字幕人妻熟人妻熟丝袜美| 国产精品精品国产色婷婷| av.在线天堂| 国产69精品久久久久777片| 国产精品人妻久久久影院| 亚洲最大成人中文| 麻豆成人午夜福利视频| 久久亚洲精品不卡| 欧美日本视频| 春色校园在线视频观看| 精华霜和精华液先用哪个| 欧美97在线视频| 99久久无色码亚洲精品果冻| 午夜福利在线在线| 久久久久久久亚洲中文字幕| 日韩成人av中文字幕在线观看| 1000部很黄的大片| 别揉我奶头 嗯啊视频| 中国国产av一级| 国产午夜福利久久久久久| 国产精品伦人一区二区| 久久精品影院6| 一个人看的www免费观看视频| 免费在线观看成人毛片| 国产成年人精品一区二区| 黄色一级大片看看| 国产黄色视频一区二区在线观看 | 18+在线观看网站| 日本一二三区视频观看| eeuss影院久久| 看黄色毛片网站| 免费观看的影片在线观看| 日本色播在线视频| a级一级毛片免费在线观看| 美女内射精品一级片tv| 亚洲欧美精品专区久久| 国产精品电影一区二区三区| 91精品伊人久久大香线蕉| 午夜日本视频在线| 99热这里只有是精品50| 岛国在线免费视频观看| 国产乱人视频| av在线老鸭窝| 日韩av在线大香蕉| 人人妻人人澡欧美一区二区| 亚洲三级黄色毛片| 欧美xxxx黑人xx丫x性爽| 精品人妻视频免费看| 免费人成在线观看视频色| 能在线免费看毛片的网站| 哪个播放器可以免费观看大片| 精品久久久久久久人妻蜜臀av| 亚洲精品国产av成人精品| 色噜噜av男人的天堂激情| 国产亚洲精品av在线| 麻豆成人午夜福利视频| 亚洲av男天堂| 黄色一级大片看看| 韩国高清视频一区二区三区|