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

    Low-Cycle Fatigue Crack-Propagation Behavior of Ship Cracked Plate Considering the Accumulative Plasticity under Variable Amplitude Loading

    2018-03-28 12:30:23DENGJunlinYANGPingMALiQIANWei
    船舶力學(xué) 2018年3期

    DENG Jun-lin,YANG Ping,MA Li,QIAN Wei

    (1a.Key Laboratory of High Performance Ship Technology(Wuhan University of Technology),Ministry of Education,Wuhan 430063,China;1b.Wuhan University of Technology,Wuhan 430063,China;2.Department of Naval Architecture&Ocean Engineering,Qinzhou University,Guangxi 535000,China)

    0 Introduction

    For many fatigue-critical parts of structures,vehicles,ship and machines,the fatigue crack propagation under service conditions generally involves random or variable amplitude,rather than constant amplitude loading conditions[1].The significant retardation or acceleration effect in fatigue crack propagation may occur as a result of the tensile/compression overload variations[2-3],where the variation of load size and load order has an important influence on the fa-tigue crack growth rate.In addition,the research results have revealed that the most overall fractures of ship hulls tend to be the coupling interaction result of low-cycle fatigue and incremental accumulative plastic failure[4-5].The actual ship hull structure under service loading may suffer from low-cycle fatigue damage,when the accumulative plastic deformation of the ship hull structures exceeds a certain value,it will lead to the initiation and expansion of the low-cycle fatigue crack of the local region of the ship hull structures,and finally until the structural failure.Therefore,it is of great practical significance to analyze the low-cycle fatigue crack growth life of ship hull under variable amplitude loading,and an accurate prediction of fatigue life requires an adequate evaluation of these load interaction effects.

    The effects of single peak tensile overloads have been reported in many investigations[6-10]simply because this type of loading can lead to significant load interaction effects.As one of the most common forms of variable amplitude loading,the low-cycle fatigue crack propagation under tensile overload has been extensively studied[6].The significant retardation effect in fatigue crack propagation may occur under the tensile overload variations.Many factors influence the retardation effect,where the most significant one is the overload ratio ROL=σOL/σmax.In addition,several mechanisms have been proposed to explain crack growth retardation,which include models based on residual stress[7],crack closure[8],crack tip blunting[9]and strain hardening[10].In spite of some controversy,the effect of residual plastic deformation,which leads to compressive stresses before the crack tip and raises the crack opening load on subsequent crack growth(crack closure),has been identified as the most important aspect in explaining the variation of the characteristic features of post-overload transients[7-10].

    Compared with the studies on the retardation effect of single peak tensile overload in fatigue crack propagation,the research on the acceleration effects of single peak compressive overload in fatigue crack propagation is relatively much fewer.Some researchers[11]pointed out that the compressive overload will cause the acceleration effects in fatigue crack propagation;Skorupa[12]pointed out that the compressive overload effect on fatigue crack propagation depended on the material and load parameters.In addition,based on the critical distance theory,Zheng[13]has proposed a numerical model based on the theory of critical distances to simulate the fatigue crack propagation of cast aluminum alloys under tensile/compressive overload.The accumulated plastic energy density at the critical distance point ahead of the crack tip is used as a measure of the fatigue damage.

    The low-cycle fatigue crack propagation is mainly in the plastic zone,where the stress field at the crack tip is quite complex.So the classical predicted model of fatigue crack propagation based on linear elastic mechanics is difficult to apply.Additionally,Yao and Kang,et al[14-15]have pointed out that the influence of accumulative incremental plasticity should be considered in the study of the low-cycle fatigue life of ship structures.The accumulative plastic deformation which will produce additional damage and reduce the low-cycle fatigue life of the hull structure is produced at the crack tip.Therefore,the low-cycle fatigue crack growth life analysis method of considering the accumulative damage at crack tip is more realistic to evaluate the overall fracture load bearing capacity of the ship hull plate members.Accordingly,the accumulative plastic deformation at crack tip decreases structure material stiffness on the crack tip plastic zone,then exacerbates the ductility loss of structure material in the crack tip plastic zone and causes the fatigue damage of the crack tip by large scale yielding under low-cycle fatigue loading.

    Therefore,it is necessary to analyze the low-cycle fatigue crack propagation of cracked plate under the variable amplitude loading based on the accumulative plasticity in the vicinity of crack tip.As mentioned above,this study based on the cyclic stress-strain field of crack tip,continues to use the prediction model proposed in previous study[1]to evaluate the low-cycle fatigue crack growth of cracked plate under the tensile/compressive overload of the variable amplitude low-cyclic fatigue loading,considering the effect of the accumulative plastic deformation of crack tip.

    1 Prediction model of low-cycle fatigue crack growth

    In this study,the low-cycle fatigue crack growth behavior of cracked plate under variable amplitude loading is studied based on the proposed model in the previous study[1].The maximum plastic zone radius and the cyclic plastic zone radius at crack tip under variable amplitude loading are shown in Eq.(1);If the maximum stress intensity factor (Kmax)and the stress intensity factor()K at the crack tip are known,the cycle plastic zone radius and the maximum plastic zone radius can be expressed as follows[16]:

    Fig.1 shows that the plate element structure in the monotonic plastic zone at the crack tip is composed of a series of micro fatigue unit δρ.The amount of the fatigue damage for micro fatigue units δρ in the monotonic plastic zone at the crack tip caused by each cycle load is about 1/Ni.After n time cycles,the total amount of the fatigue damage crack tip for micro fatigue units δρ in the cycle plastic zone at the crack tip is up toWhen the accumulative fatigue damage D for micro fatigue units δρ in the monotonic plastic zone at the crack tip reaches the critical value of D0,the crack tip extends a plastic zone size.So the accumulative fatigue damage for micro fatigue units δρ in the monotonic plastic zone is calculated by each fatigue unit under variable amplitudes loading.

    Fig.1 The element model in front of crack tip

    The current crack length aican be considered as a function of the initial crack size a0,the numbers of damaged fatigue units naand the width of the residual ligament at the crack tip of the current unit rl,as shown in Fig.1.

    At the beginning of the low-cycle fatigue crack growth,the width of the residual ligament at the crack tip of the current unit rl=δρ,and the number of damaged fatigue units in the plastic zone near the crack tip na=0,So the current crack size is equal to the initial crack size ai=a0;at the same time,all the crack damage variables of fatigue unit at the initial time were zero.

    According to Fig.1,when the accumulative damage of the ith fatigue units in the plastic zone at the crack tip is calculated,the distance rifrom the fatigue unit to the current crack tip under the variable amplitude low-cycle fatigue loading can be expressed as follows:

    After obtained the distance expression ri,the Eq.(3)is substituted into the previous study[1]to solve the distribution of the stress-strain field at the crack tip and the amount of fatigue damage under the ith fatigue cycle load.In the Eq.(3),the width of the residual ligament at the crack tip of the current fatigue unit rl can be obtained through the linear interpolation of the damage parameters of the front and rear boundary of the ith fatigue unit in the monotonic plastic zone at the crack tip,which is expressed as below:

    The Eqs.(2)and(3)are substituted into the previous study[1],which have obtained the stress-strain field and the fatigue damage of the crack tip under constant fatigue load,to solve respectively the distribution of the stress-strain field at the crack tip and the amount of fatigue damage under variable amplitude fatigue loading.

    The maximum plastic zone radius Eq.(1)under variable amplitude fatigue load is substitut-ed into the proposed model in the author’s previous study[1],which can be used to obtain the low-cycle fatigue crack propagation rate under variable amplitude low-cycle fatigue loading as follows:

    The Eq.(5)is the predicted model for the low-cycle fatigue crack growth of cracked plate under variable amplitude low-cycle fatigue loading,which considered the effect of accumulative plastic deformation at crack tip.By using the proposed prediction model Eq.(5),the lowcycle fatigue crack growth rate of hull cracked plate is analyzed respectively under the single peak tensile/compressive overload of variable amplitude loading in this study.

    2 Experimental investigations

    2.1 Experimental setup

    The fatigue crack growth experiments were made on 7 mm thick hull notch cracked plate specimens with other dimensions shown in Fig.2.The specimen and the test setup are shown in Fig.3.The steel Q235,is a low carbon steel which is widely used in hull structures,was employed to the tests,where the uniaxial tensile stress-strain curve of Q235 steel is obtained by tensile test as shown in Fig.4.And the basic material mechanics properties of Q235 steel are obtained by tensile test as shown in Tab.1.Chemical composition(in%wt)of this material is:C 0.18,Si 0.43,Mn 1.4,P 0.02,S 0.014 and V 0.13.The crack planes in the specimens were orientated perpendicular to the rolling direction of the plates.The fatigue crack growth experiments were performed in air and at room temperature using a computer controlled servohydraulic test machine,MTS322 250 kN.The crack length was simultaneously measured using strain extensometer with a 10 mm gauge length,±1 mm range,and 0.01%extensometer strain control accuracy whose operating temperature ranges from_80_C to 200_C,as shown in Fig.3.A uniaxial controlled load with maximum applied load was used to control the total stress range and a tensile-compression loading with a triangular waveform was used to ensure that the strain rate remained constant in a loop beginning with the tensile load.

    Fig.2 A notch cracked plate specimen(notch diameter is 4.8 mm,thickness is 7mm)

    A 0.25 Hz frequency was chosen based on other low-cycle fatigue test results available in the literature.The specimens were tested in cyclic loading,considered the influence of the compressive stress.The tests were conducted under constant amplitude loading with different stress ratios and mean stress to study the low-cycle fatigue crack growth.Each test was performed to failure and a mean of three measurements under each condition was taken.

    The real-time information such as time,load,crack tip opening displacement and extensometer strain were recorded.After the fatigue failure,fractographs and microscopic changes were examined by using LEO Electron Microscopy 1530 VP field emission scanning electron microscope(FESEM)to analyze the failure mechanisms.An example of the fracture process and fracture morphologies of low-cycle fatigue crack propagation is shown in Fig.5.

    Fig.3 Test setup diagram of notch cracked plate specimen

    Fig.4 Tensile stress strain curve of Q235 steel

    Tab.1 The mechanics properties of Q235

    First by sticking to the crack tip strain we measured to get accumulative plastic deformation at the crack tip and the corresponding relationship between the deformation rate and the number of cycles under different stress ratio and mean stress and provide validation for the verification of the model;then by the scale line prefabricated on the initial crack tip and combined with distribution in the sample on the high times magnifying glass,we can observe and record a crack moving cycles of a tick spacing and obtain low-cycle fatigue crack growth rate curve.

    By connecting the projector with the magnifying glass,it is clear that the position of the current crack tip can be obtained clearly and timely.When the fatigue load is stuck in the cyclic vale,the horizontal coordinate of the crack tip is recorded,and the crack length corresponding to the cycle can be obtained by subtracting the horizontal coordinate from the prior determination of the fatigue crack initiation and crack.The crack growth rate can be obtained by the increment of the crack length divided by the number of load cycles experienced by two adjacent crack measurements:

    where ai+1and aiare measured value of two adjacent crack length.Crack propagation velocity calculated by Eq.(6)is mean velocity of extension from aito ai+1.ΔKeffalso uses the mean crack lengthof two measurements;specific calculation expressions can refer to Ref.[13].

    The effective stress intensity factor range of the crack tip can be obtained by the following equation[13]:

    2.2 Fracture morphology analysis

    Fig.5 shows the fracture process and fracture morphology of low-cycle fatigue crack growth for specimen.It can be seen the obvious crack closure phenomenon in Fig.5(a),which is the crack morphology before overloading;the opening displacement at the crack tip increases obviously in Fig.5(b)after the overload,and the bifurcation can be emerged at crack tip,which caused two new crack tip.It is the morphology of a crack tip at unloading stage after applying an overload peak in Fig.5(c).It can not be seen the obvious crack closure in Fig.5(c),when the loading is unloaded to the minimum value;meanwhile,the acceleration effect can be produced at crack tip with the expansion of the low-cycle fatigue crack.It is the morphology of a crack tip expanded a short distance after overloading in Fig.5(d).Contrasted to Fig.5(a),the crack closure size caused by the tensile overload peak is larger than that by the constant loading in Fig.5(d).Moreover,the low-cycle fatigue crack growth continues to expand in the original expansion direction after the segregation occurred.The fracture morphology in the fatigue test showed transgranular fracture in Fig.5(e).It can be found that the accumulative plastic deformation was remarkably at the crack tip,and then the larger damage was produced in the process of cyclic loading,which led to the decrease of the ability to resist deformation for the material.With the increasing of stress ratio,the larger accumulative plastic deformation will lead to earlier failure.When the accumulative plastic deformation in front of crack tip is up to the critical value,the crack tip is gradually moved downward from a vertical scale to the next one with performance for the crack propagation.

    Fig.5 Fracture process and fracture morphology of low-cycle fatigue crack growth

    Therefore,it is necessary to consider the effect of accumulative incremental plastic damage at the Gauss integral point of the low-cycle fatigue crack tip in the study of the low-cycle fatigue crack growth rate of the cracked plate under asymmetrical stress cyclic loading.

    2.3 Experimental result and discussion

    According to the conventional mechanical properties and the material parameters of Q235 steel obtained from tests,the low-cycle fatigue crack growth rate of cracked plate under different low-cycle fatigue load is predicted by using the prediction model proposed in this study,and the analyzed results are finally compared with the experiment as follows:

    2.3.1 Low-cycle fatigue crack growth rate under single peak tensile overloads

    After the single peak tensile overload introduced into the constant amplitude low-cycle fatigue load,for the different stress ratios R,the mean stresses and tensile overload ratios OLR,the results of the comparison between the predicted value of the low-cycle fatigue crack growth rate used by the proposed model in this study and the experimental data are shown in Fig.6.

    It has been shown in Fig.6 that the retardation effect of the low-cycle fatigue crack propagation occurs under the tensile overload.The larger the tensile overload ratio OLR is,the more obvious in the retardation effect of the low-cycle fatigue crack propagation of the cracked plate is with the constant stress ratio and mean stress,and the relatively lower of the minimum low-cycle fatigue crack growth rate in retardation stage is.The length of the retardation period of the low-cycle fatigue crack propagation is also longer.Alike,the larger the stress ratio R is,the more obvious in the retardation effect of the low-cycle fatigue crack propagation of the cracked plate is with the same tensile overload ratio OLR and mean stress,and the relatively lower of the minimum low-cycle fatigue crack growth rate in retardation stage also is.The mean stress has a similar influence on the retardation effect of the low-cycle fatigue crack growth of the ship with the tensile overload ratio.

    Fig.6 Comparison results between the prediction and the experiment of low-cycle fatigue crack growth rate for different tensile overloads

    The mechanism of the retardation effect on the low-cycle fatigue crack growth rate of cracked plate with tensile overload is explained as follows:According to the calculation results of residual stress with the tensile overload,the accumulation progressive plastic deformation of crack tip significantly increased,which caused the large residual compressive stress in the vicinity of crack tip.The residual compressive stress decreases the effective stress intensity factor range of the crack tip,which reduces the low-cycle fatigue crack growth rate,and then produces the retardation effect.When the crack tip gradually exceeds out the forward plastic zone of crack tip which induced by the single peak tensile overload,the residual stress at crack tip gradually disappear,so the low-cycle fatigue crack growth rate gradually recovers the level of the constant amplitude loading.According to the experiment and the finite element analysis,Jiang[17]pointed out that the increase of the fatigue crack opening stress after the tensile overload peak led to the decrease of the effective stress intensity factor range,which produced the retardation effect of the low-cycle fatigue crack growth.

    2.3.2 Low-cycle fatigue crack growth under single peak compressive overload

    After the single peak compressive overload introduced into the constant amplitude lowcycle fatigue load,for the different cases of the stress ratio R,the mean stress and compressive overload ratios ULR,the results of the comparison between the predicted value of the lowcycle fatigue crack growth rate used by the proposed model in this study and the experimental data are shown in Fig.7.

    Fig.7 Comparison results between the prediction and the experiment of low-cycle fatigue crack growth rate for different compressive overloads

    It has been shown in Fig.7 that the acceleration effect of the low-cycle fatigue crack propagation occurs under the compressive overload.The larger the compressive overload ratio ULR is,the more obvious in the acceleration effect of the low-cycle fatigue crack propagation of the cracked plate is with the constant stress ratio and mean stress,and the relatively lower of the minimum low-cycle fatigue crack growth rate in retardation stage is.The length of the acceleration stage of the low-cycle fatigue crack propagation is also longer.Alike,the larger the stress ratio R is,the more obvious in the acceleration effect of the low-cycle fatigue crack propagation of the cracked plate is with the same compressive overload ratio ULR and mean stress,and the relatively lower of the minimum low-cycle fatigue crack growth rate in retardation stage also is.The mean stress has a similar influence on the acceleration effect of the lowcycle fatigue crack growth of the ship with the compressive overload ratio.

    The mechanism of the acceleration effect on the low-cycle fatigue crack growth rate of cracked plate with compressive overload is explained as follows:The residual stress of the crack tip is offset the compressive stress caused by the tensile load,which produced the acceleration effect of the fatigue crack growth under the single peak compressive overload.Meanwhile,the significantly increase of the accumulative plastic deformation of the crack tip in the plastic zone produced a larger reverse plastic yield region,which promoted the loss of ductility in the vicinity of the crack tip.A large surplus accumulated plastic damage had been produced in the vicinity of crack tip under compressive overload,which promoted the expansion of low-cycle fatigue crack.Additionally,the large surplus accumulative plastic damage of the crack tip promotes the low-cycle fatigue crack propagation before the crack tip extended out the reverse plastic yield zone after single peak compression overload;Then,with the crack tip gradually extending out the reverse plastic yield zone,the low-cycle fatigue propagation rate is gradually restored to the level of steady expansion before the compression overload.

    It can be seen from the above figures that the prediction values for three different influence factors under single peak tensile/compressive overload agree well with the experiment.The prediction model proposed in this study considering the increasing accumulative plasticity of crack tip can better reflect the low-cycle fatigue crack growth rate of hull cracked plate under single peak tensile/compressive overload,which has the feasibility and higher prediction accuracy.

    3 The low-cycle fatigue crack growth life prediction

    According to the prediction model Eq.(5),the calculated model of the low-cycle fatigue crack growth life of cracked plate under variable amplitude cyclic loading with a single peak is expressed as follows:

    where a0is the critical size of low-cycle fatigue crack initiation,and the critical crack initiation size a0=0.5 mm is selected in this study;afis the fracture failure size in fatigue crack propagation test of cracked plate.

    The low-cycle fatigue crack growth life of hull cracked plate under the tension/compression loading with a single peak is analyzed by proposed model Eq.(8)in this study.Confined to the space,the predicted values of the low-cycle fatigue crack growth life with different tension/compression overload ratios are only displayed to compare with the experimental results.

    As shown in Fig.8,the bigger the tensile overload ratio is,the more remarkable the increase of the low-cycle fatigue crack growth life caused by the tensile overload peak with the same mean stress and stress ratio;Similarly,the larger the compressive overload ratio is,the more obvious the decrease of the low-cycle fatigue crack growth life caused by the compressive overload peak.Finally,the predictions in this study are in good agreement with the experimental results.

    Fig.8 The comparison results of low-cycle fatigue crack growth life between the prediction andtest under different variable amplitude loadings

    The predicted model Eq.(8)can be well used to evaluate the low-cycle fatigue crack propagation life under the tensile/compressive overload,based on the proposed model of the lowcycle fatigue crack growth rate in this study,which have a good accuracy.Confined to the space and the same prediction relationship,the predicted values of the low-cycle fatigue crack growth life with different cases,such as different mean stresses and stress ratios,and so on,are not displayed repeatly in this study.

    4 Conclusions

    The low-cycle fatigue crack propagation behavior of ship cracked plate is studied under variable amplitude loading considering the effect of accumulative plasticity in this study.The lowcycle fatigue crack growth behavior of ship cracked plate can well evaluated by the predicted model proposed in this study under the tensile/compressive overload with a single peak.

    Additionally,the obvious retardation effect emerged at the crack tip of cracked plate with the Q235 steel after the tensile overload peak,and then the low-cycle fatigue crack propagation gradually recovered the level of the constant amplitude loading.The larger the tensile overload ratio is,the more remarkable the retardation effect is;however,the stress ratio has an opposite influence on the retardation effect.The acceleration effect will be produced under the compressive overload,but the affect range of the acceleration effect is much lower than the retardation effect.

    [1]Deng J L.Research on low-cycle fatigue strength of ship hull plate based on accumulated plastic damage[D].Wuhan:Wuhan University of Technology,2016.(in Chinese)

    [2]Ding Zhenyu.A study on intrinsic mechanism of crack propagation behavior subjected overload[D].Hangzhou:Zhejiang University of Technology,2013.(in Chinese)

    [3]Tvergaard V.Effect of underloads or overloads in fatigue crack growth by crack-tip blunting[J].Engineering Fracture Mechanics,2006,73:869-79.

    [4]Mansour A E,et al.An experimental investigation of ship hull ultimate strength[J].Trans.SNAME,1968,98:411-439.

    [5]Murray J M.Structural development of tankers[J].European Shipbuilding,1953,3(5):69-75.

    [6]Skorupa M.Load interaction effects during fatigue crack growth under variable amplitude loading-a literature review[J].Part II:Qualitative Interpretation.Fatigue Fract.Eng.Mater.Struct.,1999,22:905-26.

    [7]Shijve J,Broek D.The result of a test program based on a gust spectrum with variable amplitude loading[J].Aircraft Engng,1962,34:314-316.

    [8]Meggiolaro M A,de Castro J T P.On the dominant role of crack closure on fatigue crack growth modeling[J].International Journal of Fatigue,2003,25(9-11):843-854.

    [9]Christensen R H.Fatigue crack,fatigue damage and their detection[M].Metal Fatigue.New York:MacGraw-Hill,1959.

    [10]Jones R E.Fatigue crack growth retardation after single-cycle peak overload in Ti-6Al-4V titanium alloy[J].Engng.Fract.Mech.,1973,5:585-604.

    [11]Silva F S.The importance of compressive stresses on fatigue crack propagation rate[J].International Journal of Fatigue,2005,27(10-12):1441-1452.

    [12]Skorupa M.Load interaction effects during fatigue crack growth under variable amplitude loading-a literature review.Part I:Empirical trends[J].Fatigue Fract.Engng.Mater Struct.,1998,21:985-1006.

    [13]Zheng X,Cui H,Su X.Numerical modeling of fatigue crack propagation based on the theory of critical distances[J].Engineering Fracture Mechanics,2013,114:151-165.

    [14]Yao T,Nikolov P I.Buckling/plastic collapse of plates under cyclic loading[J].J Soc.Naval Arch Japan,1990(168):449-462.

    [15]Kang G Z,Gao Q,Cai L X,et al.Experimental study on the uniaxial and nonproportionally multiaxial of SS304 stainless steel at room and high temperatures[J].Unclear Engineering and Design,2002,216:13-26.

    [16]Creager M,Paris P C.Elastic field equations for blunt cracks with reference to stress corrosion cracking[J].Int.J Fract.Mech.,1967,3:247-52.

    [17]Jiang Y Y,Feng M L,Ding F.A reexamination of plastivity-induced crack closure in fatigue crack propagation[J].International Journal of Plasticity,2005,21(9):1720-1740.

    搡老妇女老女人老熟妇| 免费不卡的大黄色大毛片视频在线观看 | 一级毛片我不卡| 日本欧美国产在线视频| 国产亚洲欧美98| 精品国内亚洲2022精品成人| 小说图片视频综合网站| 18禁黄网站禁片免费观看直播| 久久久精品94久久精品| 久久久色成人| 在线播放国产精品三级| 黄色一级大片看看| 亚洲精品一区av在线观看| 亚洲内射少妇av| 亚洲最大成人手机在线| 中文字幕免费在线视频6| 国产午夜精品论理片| 男女啪啪激烈高潮av片| 久久精品影院6| 美女cb高潮喷水在线观看| 国产91av在线免费观看| 久久精品国产亚洲av涩爱 | 床上黄色一级片| 亚洲最大成人手机在线| 中文在线观看免费www的网站| 国产精品国产高清国产av| 亚洲乱码一区二区免费版| 国产一区二区激情短视频| 亚洲av免费高清在线观看| 国产在视频线在精品| 白带黄色成豆腐渣| 美女免费视频网站| 免费看美女性在线毛片视频| 婷婷精品国产亚洲av在线| 国产一区亚洲一区在线观看| 亚洲欧美清纯卡通| 国产久久久一区二区三区| 看免费成人av毛片| 无遮挡黄片免费观看| 精品久久国产蜜桃| 99久国产av精品国产电影| 国产大屁股一区二区在线视频| 国产黄片美女视频| 一a级毛片在线观看| 精品免费久久久久久久清纯| 99热这里只有是精品在线观看| 看片在线看免费视频| 国产成人精品久久久久久| 精品一区二区免费观看| 国产精品久久视频播放| 欧美极品一区二区三区四区| 欧美绝顶高潮抽搐喷水| 深夜a级毛片| 超碰av人人做人人爽久久| 日本黄色片子视频| 欧美人与善性xxx| 人人妻人人澡欧美一区二区| 乱人视频在线观看| 热99re8久久精品国产| 成人精品一区二区免费| 男女啪啪激烈高潮av片| 淫秽高清视频在线观看| 最近2019中文字幕mv第一页| 国产在线男女| 国产精品综合久久久久久久免费| 哪里可以看免费的av片| 美女cb高潮喷水在线观看| 久久久a久久爽久久v久久| 成年女人看的毛片在线观看| 99热只有精品国产| 精品久久久久久久末码| 亚洲精品乱码久久久v下载方式| 不卡视频在线观看欧美| 国产午夜精品久久久久久一区二区三区 | 欧美色视频一区免费| 久久久精品大字幕| 国产亚洲精品久久久久久毛片| 国产成人91sexporn| 成人无遮挡网站| 日韩欧美国产在线观看| 国产视频一区二区在线看| 久久中文看片网| 99riav亚洲国产免费| 十八禁国产超污无遮挡网站| 成年女人看的毛片在线观看| 免费高清视频大片| 亚洲国产色片| 国产激情偷乱视频一区二区| 午夜福利高清视频| 久久精品综合一区二区三区| 久99久视频精品免费| 午夜视频国产福利| 最好的美女福利视频网| 观看美女的网站| 成年女人毛片免费观看观看9| 国产精品国产高清国产av| 搡女人真爽免费视频火全软件 | 免费av观看视频| 嫩草影视91久久| 成人漫画全彩无遮挡| 91精品国产九色| 国产精品国产高清国产av| 国产蜜桃级精品一区二区三区| 日韩高清综合在线| 成人综合一区亚洲| 91午夜精品亚洲一区二区三区| 亚洲最大成人av| 国产精品伦人一区二区| 大香蕉久久网| 亚洲综合色惰| 99久久精品国产国产毛片| 99热精品在线国产| 久久欧美精品欧美久久欧美| 成人漫画全彩无遮挡| 如何舔出高潮| 国产白丝娇喘喷水9色精品| 国产成人91sexporn| 亚洲四区av| av专区在线播放| 此物有八面人人有两片| 日韩成人av中文字幕在线观看 | 欧美精品国产亚洲| 成人性生交大片免费视频hd| 97超碰精品成人国产| 村上凉子中文字幕在线| 国产高清有码在线观看视频| 欧美极品一区二区三区四区| av天堂中文字幕网| 日日撸夜夜添| 亚洲欧美日韩高清专用| 久久久久久久久中文| 精品人妻视频免费看| 亚洲高清免费不卡视频| 久久久a久久爽久久v久久| 久久久久久九九精品二区国产| 一级毛片电影观看 | 最近视频中文字幕2019在线8| 国产亚洲欧美98| 内射极品少妇av片p| 久久这里只有精品中国| 最近手机中文字幕大全| 真人做人爱边吃奶动态| 国产一区二区三区av在线 | 午夜精品国产一区二区电影 | 久久精品91蜜桃| 亚洲熟妇中文字幕五十中出| 99热全是精品| 国产探花极品一区二区| 日韩制服骚丝袜av| 99久国产av精品国产电影| 国产激情偷乱视频一区二区| 国产精品无大码| 天天躁日日操中文字幕| 午夜久久久久精精品| 亚洲欧美精品自产自拍| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 国产精品乱码一区二三区的特点| 久久久成人免费电影| 99久久精品国产国产毛片| 中文字幕av在线有码专区| 亚洲国产精品sss在线观看| 亚洲精品久久国产高清桃花| 精品久久久噜噜| 午夜福利在线在线| 女人十人毛片免费观看3o分钟| 十八禁国产超污无遮挡网站| 久久欧美精品欧美久久欧美| 国产成人a∨麻豆精品| 国产精品一区二区免费欧美| 久久久午夜欧美精品| 国产精品精品国产色婷婷| 校园春色视频在线观看| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 性欧美人与动物交配| 波野结衣二区三区在线| 国产一区二区在线av高清观看| 秋霞在线观看毛片| 亚洲国产精品合色在线| 亚洲一级一片aⅴ在线观看| 国产人妻一区二区三区在| 久久精品国产亚洲av涩爱 | 日日摸夜夜添夜夜爱| 国产在视频线在精品| 国产爱豆传媒在线观看| 精品免费久久久久久久清纯| 国产又黄又爽又无遮挡在线| 91午夜精品亚洲一区二区三区| 女人十人毛片免费观看3o分钟| 一本精品99久久精品77| 91久久精品电影网| 日日干狠狠操夜夜爽| av在线老鸭窝| 日产精品乱码卡一卡2卡三| 99久久无色码亚洲精品果冻| 可以在线观看毛片的网站| 美女大奶头视频| 久久久久国内视频| 看片在线看免费视频| eeuss影院久久| 1024手机看黄色片| 欧美+日韩+精品| 丰满的人妻完整版| 国产淫片久久久久久久久| 伦精品一区二区三区| 丰满的人妻完整版| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 国产成人福利小说| 人妻久久中文字幕网| 99久久久亚洲精品蜜臀av| 校园人妻丝袜中文字幕| 乱码一卡2卡4卡精品| 婷婷精品国产亚洲av在线| 国产精品久久久久久av不卡| 美女被艹到高潮喷水动态| 精品久久久久久久久亚洲| 久久久色成人| 又黄又爽又免费观看的视频| 亚洲精品粉嫩美女一区| 欧美成人精品欧美一级黄| 午夜福利在线观看免费完整高清在 | 日本免费a在线| 国产精品一区二区免费欧美| 自拍偷自拍亚洲精品老妇| 成人无遮挡网站| 赤兔流量卡办理| 热99在线观看视频| 熟女人妻精品中文字幕| 国产黄片美女视频| a级毛片a级免费在线| 亚洲图色成人| 久久热精品热| 久久韩国三级中文字幕| 中文字幕av在线有码专区| 日本 av在线| 亚洲成人久久爱视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品爽爽va在线观看网站| 国产乱人偷精品视频| 97碰自拍视频| 一区二区三区免费毛片| 国产亚洲精品综合一区在线观看| 国产蜜桃级精品一区二区三区| 国产美女午夜福利| 亚洲av二区三区四区| 日韩强制内射视频| 国产精品永久免费网站| 此物有八面人人有两片| 美女cb高潮喷水在线观看| 亚洲,欧美,日韩| 国产一区亚洲一区在线观看| 中出人妻视频一区二区| 亚洲成人久久性| 欧美成人a在线观看| 国产精品不卡视频一区二区| 99热只有精品国产| 蜜桃亚洲精品一区二区三区| 人妻夜夜爽99麻豆av| 亚洲av第一区精品v没综合| 人妻丰满熟妇av一区二区三区| www日本黄色视频网| 在线免费十八禁| 亚洲成人av在线免费| 激情 狠狠 欧美| 观看免费一级毛片| 搡老岳熟女国产| 成人av在线播放网站| 国语自产精品视频在线第100页| 久久久久久大精品| 精品一区二区三区av网在线观看| 国内揄拍国产精品人妻在线| 韩国av在线不卡| 三级国产精品欧美在线观看| 久久韩国三级中文字幕| 亚洲人成网站在线观看播放| 国产成年人精品一区二区| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 亚洲精品国产av成人精品 | 插阴视频在线观看视频| 熟妇人妻久久中文字幕3abv| 国产国拍精品亚洲av在线观看| 又爽又黄a免费视频| 亚洲自拍偷在线| 麻豆成人午夜福利视频| 国产高清不卡午夜福利| 欧美日韩乱码在线| 国产精品野战在线观看| 久久久久久伊人网av| 99久国产av精品| 男人舔奶头视频| 亚洲欧美日韩高清专用| h日本视频在线播放| 听说在线观看完整版免费高清| 国内精品久久久久精免费| 精品一区二区免费观看| 晚上一个人看的免费电影| 91午夜精品亚洲一区二区三区| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 精品久久国产蜜桃| av在线老鸭窝| 久久99热这里只有精品18| 国产精品久久久久久久久免| 精品人妻偷拍中文字幕| 最后的刺客免费高清国语| 久久精品国产亚洲av涩爱 | 中文字幕精品亚洲无线码一区| 中文亚洲av片在线观看爽| 成人午夜高清在线视频| 国产伦在线观看视频一区| 国产精品日韩av在线免费观看| 久久久久久久亚洲中文字幕| 婷婷六月久久综合丁香| 亚洲精品456在线播放app| 国产极品精品免费视频能看的| 成年女人永久免费观看视频| 久久中文看片网| 亚洲成人av在线免费| 欧美bdsm另类| 男人的好看免费观看在线视频| 国产真实乱freesex| 国产又黄又爽又无遮挡在线| av福利片在线观看| 亚洲精品日韩av片在线观看| 亚洲av一区综合| 99久国产av精品| 你懂的网址亚洲精品在线观看 | 亚洲国产精品国产精品| 日本黄大片高清| 中文字幕熟女人妻在线| 午夜福利成人在线免费观看| 在线天堂最新版资源| 国产av不卡久久| 美女高潮的动态| 变态另类成人亚洲欧美熟女| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 中国国产av一级| 中文在线观看免费www的网站| 精品熟女少妇av免费看| 一区二区三区免费毛片| 又黄又爽又免费观看的视频| 一级毛片久久久久久久久女| 中国美白少妇内射xxxbb| 中文字幕熟女人妻在线| 国产欧美日韩精品一区二区| 国产高清视频在线观看网站| 亚洲在线观看片| 国内揄拍国产精品人妻在线| 干丝袜人妻中文字幕| 久久久午夜欧美精品| 最近2019中文字幕mv第一页| 国产精品综合久久久久久久免费| 国产片特级美女逼逼视频| 日本黄大片高清| 亚洲最大成人手机在线| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 亚洲成人av在线免费| 六月丁香七月| 国产白丝娇喘喷水9色精品| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放 | 波多野结衣高清作品| 久久久久久久亚洲中文字幕| 波多野结衣高清作品| 天堂网av新在线| 精品午夜福利在线看| 精品久久久久久久久久免费视频| 深夜精品福利| 99热这里只有是精品50| 亚州av有码| 中文在线观看免费www的网站| 婷婷六月久久综合丁香| 日本在线视频免费播放| 99久久精品热视频| 中文在线观看免费www的网站| 欧美日韩在线观看h| 色视频www国产| 国产午夜精品久久久久久一区二区三区 | 精品国内亚洲2022精品成人| 床上黄色一级片| 亚洲人成网站高清观看| 搡老岳熟女国产| 欧美成人一区二区免费高清观看| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久com| 伦理电影大哥的女人| 日本五十路高清| 亚洲国产精品成人久久小说 | 婷婷精品国产亚洲av| 美女xxoo啪啪120秒动态图| 国产久久久一区二区三区| 亚洲成人精品中文字幕电影| 性色avwww在线观看| 人妻制服诱惑在线中文字幕| 亚洲精品粉嫩美女一区| 亚洲av不卡在线观看| 欧美又色又爽又黄视频| 亚洲国产精品成人久久小说 | 亚洲国产高清在线一区二区三| 我要看日韩黄色一级片| 国产大屁股一区二区在线视频| 欧美日韩在线观看h| 日韩亚洲欧美综合| 内射极品少妇av片p| 久久精品影院6| 亚洲不卡免费看| 日日啪夜夜撸| 麻豆国产av国片精品| 免费无遮挡裸体视频| 亚洲精品粉嫩美女一区| 亚洲人成网站在线播| 日韩av不卡免费在线播放| 欧美在线一区亚洲| 午夜日韩欧美国产| 国产精品,欧美在线| 国产亚洲欧美98| 三级国产精品欧美在线观看| 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 97人妻精品一区二区三区麻豆| 好男人在线观看高清免费视频| 白带黄色成豆腐渣| 成人漫画全彩无遮挡| 大又大粗又爽又黄少妇毛片口| avwww免费| 国产日本99.免费观看| 精品一区二区三区视频在线观看免费| av视频在线观看入口| 免费av观看视频| 亚洲av美国av| 亚洲久久久久久中文字幕| 国产乱人视频| 亚洲最大成人手机在线| 精品少妇黑人巨大在线播放 | 三级男女做爰猛烈吃奶摸视频| 天美传媒精品一区二区| 亚洲丝袜综合中文字幕| 久久这里只有精品中国| 啦啦啦啦在线视频资源| 老熟妇仑乱视频hdxx| 99久久无色码亚洲精品果冻| 国产精品电影一区二区三区| 久久久成人免费电影| 国产精品久久久久久久久免| 成人av一区二区三区在线看| 日韩欧美精品免费久久| 欧美人与善性xxx| 亚洲在线观看片| 成年女人看的毛片在线观看| 天美传媒精品一区二区| 啦啦啦观看免费观看视频高清| 成人国产麻豆网| 国产精品一二三区在线看| 亚洲第一区二区三区不卡| 女人十人毛片免费观看3o分钟| 女的被弄到高潮叫床怎么办| 色哟哟哟哟哟哟| 国产精品久久视频播放| 天天一区二区日本电影三级| 国产男人的电影天堂91| 高清毛片免费看| 亚洲美女黄片视频| 国产美女午夜福利| av.在线天堂| 免费搜索国产男女视频| 亚洲精华国产精华液的使用体验 | 热99在线观看视频| 中文字幕精品亚洲无线码一区| 亚洲欧美精品自产自拍| 国产精品一区二区性色av| av在线播放精品| 亚洲,欧美,日韩| 亚洲无线在线观看| 观看免费一级毛片| 亚洲欧美成人综合另类久久久 | 久久中文看片网| 精品日产1卡2卡| 亚洲成人精品中文字幕电影| 熟女人妻精品中文字幕| 精品福利观看| 免费观看人在逋| 少妇的逼水好多| 亚洲自偷自拍三级| 欧美国产日韩亚洲一区| 少妇人妻一区二区三区视频| 成人性生交大片免费视频hd| 亚洲电影在线观看av| 欧美成人一区二区免费高清观看| 内射极品少妇av片p| 99九九线精品视频在线观看视频| 天堂网av新在线| 少妇被粗大猛烈的视频| 欧美色欧美亚洲另类二区| 12—13女人毛片做爰片一| 最近手机中文字幕大全| 国产白丝娇喘喷水9色精品| 我要搜黄色片| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 啦啦啦观看免费观看视频高清| avwww免费| 国产精品伦人一区二区| 日韩人妻高清精品专区| 国产淫片久久久久久久久| 亚洲欧美精品自产自拍| 国产成人影院久久av| 亚洲经典国产精华液单| 精品免费久久久久久久清纯| 免费av观看视频| 丰满乱子伦码专区| 亚洲三级黄色毛片| 十八禁网站免费在线| 老女人水多毛片| 我的女老师完整版在线观看| 免费黄网站久久成人精品| 日本与韩国留学比较| 免费在线观看成人毛片| 中文字幕av成人在线电影| 久久久色成人| 成年女人看的毛片在线观看| 春色校园在线视频观看| 免费av不卡在线播放| 欧美zozozo另类| 日韩国内少妇激情av| 黄色视频,在线免费观看| 少妇的逼水好多| 国国产精品蜜臀av免费| 美女被艹到高潮喷水动态| 亚洲av免费在线观看| 国内久久婷婷六月综合欲色啪| 午夜影院日韩av| 麻豆乱淫一区二区| 男人舔奶头视频| 亚洲欧美日韩高清专用| 国产一级毛片七仙女欲春2| 日本一二三区视频观看| 在线天堂最新版资源| 最后的刺客免费高清国语| 看片在线看免费视频| 美女免费视频网站| 晚上一个人看的免费电影| 麻豆久久精品国产亚洲av| 欧美zozozo另类| 伊人久久精品亚洲午夜| 女人十人毛片免费观看3o分钟| 成人高潮视频无遮挡免费网站| 嫩草影院新地址| 欧美激情久久久久久爽电影| av在线播放精品| 国产精品一区二区三区四区免费观看 | 亚洲四区av| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品久久国产高清桃花| 成人亚洲欧美一区二区av| 久久午夜福利片| 三级男女做爰猛烈吃奶摸视频| 俺也久久电影网| 亚洲中文字幕一区二区三区有码在线看| 国产亚洲精品久久久com| 在线天堂最新版资源| 午夜福利视频1000在线观看| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 国产视频一区二区在线看| 午夜免费激情av| 国内精品美女久久久久久| 一边摸一边抽搐一进一小说| 国产高清不卡午夜福利| 97热精品久久久久久| 女同久久另类99精品国产91| 国产精品久久久久久av不卡| 免费av不卡在线播放| 亚洲精品久久国产高清桃花| 精品熟女少妇av免费看| 麻豆成人午夜福利视频| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 亚洲av中文av极速乱| 性欧美人与动物交配| eeuss影院久久| 国产91av在线免费观看| 国产成人91sexporn| 99热这里只有是精品50| 女生性感内裤真人,穿戴方法视频| 免费一级毛片在线播放高清视频| or卡值多少钱| 久久鲁丝午夜福利片| 免费搜索国产男女视频| 熟女电影av网| 国产片特级美女逼逼视频| 午夜福利在线在线| 老司机午夜福利在线观看视频| 国产一区二区在线观看日韩| 中文在线观看免费www的网站| 国产成年人精品一区二区| 亚洲欧美日韩高清在线视频| 亚洲国产欧美人成| 搡老熟女国产l中国老女人| 一区二区三区四区激情视频 | 一个人看的www免费观看视频| 岛国在线免费视频观看| av天堂中文字幕网| 高清毛片免费观看视频网站| 12—13女人毛片做爰片一| 97超级碰碰碰精品色视频在线观看| 五月伊人婷婷丁香| 99久久成人亚洲精品观看| 日本精品一区二区三区蜜桃| 免费观看精品视频网站|