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

    Fatigue Life Prediction under Service Load Considering Strengthening Effect of Loads below Fatigue Limit

    2014-02-07 12:44:58ZHAOLihuiZHENGSonglinandFENGJinzhi

    ZHAO Lihui ,ZHENG Songlin *,and FENG Jinzhi

    1 College of Mechanical Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China

    2 College of Automobile Engineering,Shanghai University of Engineering Science,Shanghai 201620,China

    3 Machinery Industry Key Laboratory for Mechanical Strength &Reliability Evaluation of Auto Chassis Components,University of Shanghai for Science and Technology,Shanghai 200093,China

    1 Introduction*

    Accurately predicting the fatigue life according to the service loading histories is essential to the lightweight and reliability design of mechanical products.Traditionally,this was achieved based on the number of load cycles counted from loading history and load-life relations(S-N curve)from constant amplitude load(CAL)tests by the Palmgren-Miner Rule(PMR)D=ni/Ni,and when D=1 failure occurs[1].Here D is the cumulative damage,niis the mumber of cycles of load with an amplitude Si,and Niis the fatigue life under this amplitude.However,numerous studies have shown that the cumulative damage D calculated from experimental tests under variable amplitude loads often deviate from 1.SCHUTZ[2]pointed that damage sum D typically ranges from 0.3 to 3,but can reach as much as 0.1 to 10.Studies carried out by WIRSCHING,et al[3]and LEE,et al[4]also found that D may range from 0.15 to 1.06 under different loadings.This means PMR would result in poor predictions,cumulative damage less than one will lead to unsafe(non-conservative)prediction,while larger than one to conservative results[5].

    With regard to the situation that D<1,it is believed that the fatigue damage of load cycles below the fatigue limit were ignored,which was induced by the reduction of fatigue strength caused by previous amplitudes above the fatigue limit[6].So,many alternatives were proposed to cover this deficiency.Three major alternatives are:(i)modified miner’s rule(MMR),which extends the S-N curve to lower amplitudes with a slope k′=k or k′=2k-11[7];(ii)relative miner rule(RMR),which replaces the damage sum till failure D=1 with a relative smaller value D=q(q<1)[8];(iii)combination of MMR and RMR,which both extends the S-N curve and replaces the damage value,such as the design recommendation of welding structures from International Institute of Welding(IIW)[9].STANZEL,et al[10],tested AISI C1020 steel under random loadings and good agreement was found for measured and calculated results according to MMR.SCHUTZ[8]studied the results reported in literature and suggested that a damage value ranging from 0.4 to 0.7 will lead to accurate predictions for RMR.SONSINO,et al[11],conducted fatigue tests of four industrial sectors and found experimental results were all comparable to IIW Recommendations.

    Although the application of these methods effectively improved the accuracy of fatigue life prediction of this situation,they were not suitable for the situation that D>1,under which predicted lives were conservative compared with test lives.LIN,et al[12],calculated the fatigue life of motorcycle handlebar under service loadings using PMR and results show that the ratio of predicted life to test life was between 0.37 and 0.66.FRICKE,et al[13],studied the accuracy of three stress approaches for the fatigue assessment of welded structures,and pointed that predicted lives were shorter than test life except for one occasion.XU,et al[14],predicted the fatigue life of an automotive rear axle,and found that the ratio of life prediction to experiment was 0.42,0.45 and 0.59 for different location.Under these situations,the aforementioned alternatives would result in more conservative predictions.Due to a conservative prediction would lead to the safe side of a design;it did not attractive enough consideration for a long time.But with the lightweight design becomes more of a concern,it is necessary to improve the accuracy of fatigue life prediction,because too conservative a predicted life will lead to a heavy structure design because of a higher safe factor.

    The situation that D>1 was caused by an excessive calculation of damage of load cycles[6].This is because the aforementioned methods were established on the baiss that load cycles,no matter above or below the fatigue limit,are damaging.Actually,not all of the load cycles were damaging.Many studies show that the fatigue strength of the materials with strain aging could be improved by gradually increasing the applied load from a level below the fatigue limit[15–16].KAWAMOTO,et al[17]and IKAI,et al[18],investigated the strengthening effect of pre-stressing on fatigue strength using low carbon steel,and concluded that the fatigue strength was increased and gradually saturated with the increase of number of load cycles.WANG,et al[19],studied the influence of different load amplitudes on fatigue damage by means of comparative fatigue tests,and found that fatigue life will reach valley point when omitting the amplitudes around 20% of the biggest load in load spectrum.Refs.[20–22]study the strengthening effects of loads using three carbon steels with different ultimate tensile strength(UTS).The results show that the load amplitudes with strengthening effect are between 65% and 95% of the fatigue limit(low-amplitude-load strengthening,LALS),and the fatigue strength and fatigue life of the research objects are increased obviously,by 3%–22% and 80%–130%,separately.An increase in fatigue life will lead to a decrease of calculated damage,thus the accuracy of fatigue life prediction could be effectively improved using this concept.

    The purpose of this study is to propose a practical method for fatigue life prediction under service loading histories by taking into account the strengthening effect of load cycles below the fatigue limit.First,the phenomenon and mechanism of LALS were briefly reviewed.Then the mechanics model was derived and the method was established.Finally,validations were made using two examples,half shaft and torsion beam rear axle of automobile.Also fatigue lives are predicted with MMR(k′=k and k′=2k-1)for a further comparison.

    2 Fatigue Life Prediction Based on LALS

    2.1 Brief review of LALS

    Refs.[20–22]study LALS by fatigue tests of three vehicle components made of carbon steel,which are front axle(No.20 steel),half shaft(40Cr)and spur gear(20MnCr5).In the experiments,specimens were divided into two groups.For one group,fatigue test of each specimen was conducted under a load cycle above the fatigue limit(proving load,Ld),and fatigue life N0was obtained.For the other group,each specimen was first loaded by a load below the fatigue limit(strengthening load,Ls)for a certain number of cycles(strengthening number,ns),then fatigue test was conducted under the same proving load and the new fatigue life Nswas obtained.Results showed that the fatigue performances of all the three components could be significantly improved.The fatigue life under proving load after strengthened varies with Lsand ns,and can reach as much as 1.81,2.30,2.09 times of the original life for front axle,half shaft and spur gear.The relationship of strengthened life Nfvs.strengthening load Lsand strengthening number of cycles nswas shown in Fig.1.

    Using scanning electron microscopy(SEM),transmission electron microscopy(TEM),Zheng studied the microscopic mechanisms of LALS and conclusions were summarized as follows[23].

    i)Results from SEM show that the pores in the material became smaller,and secondary cracks were lesser and thinner after strengthened,see in Fig.2(a)and Fig.2(b)[23].This means the material become more uniform and compact,postponing the fatigue crack initiation and propagation,and thus resulting in an increased fatigue strength and fatigue life.

    ii)TEM results indicate that after strengthened by loads below the fatigue limit,dislocations within grains gradually moved to the grain boundaries.The accumulation of dislocation along grain boundaries forms a dislocation wall which restricts further movement of dislocations.The strong strain energy gathered in the wall increases the friction among lattices and raises the local yield stress,and finally leading to an improvement in fatigue strength and fatigue life.

    Fig.2.Microscopic mechanism of LALS

    2.2 Mechanics model of LALS

    The increase of fatigue life means an improvement in fatigue strength.However,the results can only verify that the fatigue strength was improved under a certain load level.It's necessary to know the strengthening effect under different load levels for the application of fatigue life prediction.According to Basquin relation[1],the S-N curve for high cycle fatigue(HCF)could be represented as

    where S—Stress amplitude,

    m—Inverse slope of S-N curve,

    N—Fatigue life,

    C—Fatigue constant.

    If use σ–1to present the fatigue limit,and subscript 0 and 1 to identify the original and new S-N curve after strengthened by a load below fatigue limit,the strengthening effect can be regard as fatigue strength at life N1and Ncwere increased from S0to S1and from σ–1,0to σ–1,1separately,or the fatigue life at load level S1was increased from N0to N1,as shown in Fig.3.

    From Fig.3 and Eq.(1),the following results can be obtained:

    i.e.,

    From Eq.(4),it can be easily seen that the previous S-N curve AC was parallel to the new one DE,which means the improvement of fatigue strength at different load levels are equal.So the new S-N curve,and eventually the new life can be easily calculated if knowing either of the two parameters S1or σ–1,1combining Eqs.(3)and(4).

    Taking the percentage increment of fatigue strength ΔS%to represent the strengthening effect,the new fatigue life after strengthened could be easily derived from Eq.(4),and can be expressed as

    2.3 Fatigue life prediction

    Fig.1 shows that the fatigue life under proving load will first increase then decrease with the strengthening number of cycles,and there is a best strengthening number(Nb),at which the strengthening effect reaches the maximum.This indicates the strengthening effect of loads below fatigue limit accumulates with the strengthening number nsbefore it reaches Nb,then gradually decreases to zero when nsreaches a certain value.This could be interpreted as that the loads below fatigue limit are of strengthening effect when ns<Nb,and of damaging effect when ns>Nb.

    Therefore,the process for fatigue assessment considering LALS under service loading history could be divided into two stages by the best strengthening number(Nb):stage I and stage II.During stage I,small load cycles are of strengthening effect and damage is calculated by the load cycles above the fatigue limit using a new fatigue life.In stage II,damage is calculated from both the loads above the fatigue limit and those with strengthening effect.While those smaller than strengthening loads can be ignored because they are of no effect on the fatigue strength.Combining with Miner’s Rule[1],the method for fatigue assessment of components under VAL taking into account LALS could be expressed as follows:

    where D1—Damage value in stage I,

    D2—Damage value in stage II,

    bI—Number of repetitions of load block in stage I,

    bII—Number of repetitions of load block in stage II,

    ld—Number of load levels with damaging effect,

    ls—Number of load levels with strengthening effect,

    ni—Number of cycles of the ith load level,

    Ni(j)—Fatigue life under the ith load level after strengthened by the previous j-1 load blocks.

    When D=1,failure occurs.The total life in the form of repetitions of loading block is:Nblock=bI+bII.

    Therefore,the procedure for fatigue life calculation under service loading blocks based on LALS could be performed as follows:

    i)Count the number of cycles niof each load level Siin the load history,and determine the number of levels of strengthening load lsand damaging load ld;

    ii)Calculate the number of repetitions of loading block bIto finish stage I based on lsand Nb,and calculate the fatigue life at the jth loading blockaccording to the original S-N curve and Eq.(7);

    iii)Calculate the fatigue damage during stage I D1according to Eq.(6a),and the remaining damage during stage II with D2=1-D1according to Eq.(6c);

    iv)Calculate the number of repetitions of loading block in Stage II using both the damaging and strengthening load levels according to Eq.(6b)by

    v)Calculate the total life in the form of loading block:Nblock=bI+bII.

    3 Validations

    Two examples,half shaft and rear axle,are used to validate the proposed method for the fatigue life prediction under service loading.For a further comparison,the fatigue lives of the two examples are also predicted with PMR and MMR.

    3.1 Half shaft

    The half shaft(Fig.4)was made of 40Cr.The mechanical properties of the material were as follows:ultimate tensile strength 980 MPa,yield strength 785 MPa.Service loading histories of half shaft were recorded on proving ground,and nominal stresses were calculated.The load spectrum with a sequence length Ls=21 262 was obtained after rain-flow cycle counting[24],see in Fig.5.In order to simplify calculation,an eight level block load was developed for fatigue assessment based on the rule of equivalent damage[25].The load levels and corresponding cycles were listed in Table 1.

    Fig.4.Half shaft of automobile

    Fig.5.Load spectrum of half shaft

    Table 1.Summary of block load of half shaft

    The S-N curve with Ps=50% was determined by conducting torsional fatigue tests under constant amplitude loads at a stress ratio of R=-1,and could be expressed as

    From Eq.(8),the fatigue limit τ–1(at N=2×106)of half shaft was 100.87 MPa,therefore the loads with strengthening effect range from 65.57 to 95.83 MPa(0.65–0.95τ–1).So,load levels 1 to 4 were above fatigue limit;level 5(0.84τ–1)was of strengthening effect;levels 6 to 8 were out of strengthening region and can be neglected for fatigue life prediction.

    For the mechanical characteristics of the material of half shaft were similar to that of transmission gear reported in Ref.[21],the strengthening effect of gear was taken for the fatigue life prediction of half shaft,of which the fatigue strength was increased by 4.2% after strengthened under a load at 0.84σ–1for 300 000 cycles.So the number of repetitions of load blocks to finish stage I was 208.Considering the strengthening effect to strengthening number as linear,the increment of fatigue strength of each load block was 0.020 2%,and the new fatigue life can be calculated according to Eq.(7).Then the fatigue life can be calculated according to Eq.(6):

    So the total number of repetitions of loading block to failure is 1114.

    Torsional fatigue tests of half shaft were conducted under random sequence of the eight-level block loading,Fig.6.Failure criterion was defined as the breakthrough of half shaft.Six specimens from the same batch were tested,and the fatigue lives were 1115,1289,930,1168,1395 and 1050 repetitions.The average value of experimental lives and predicted lives according to PMR,MMR and LALS are listed in Table 2.

    Fig.6.Test set-up of half shaft

    Table 2.Comparisons of predicted and test lives of half shaft

    3.2 Torsion beam rear axle

    The torsion beam rear axle was made of TL1114 and comprised of eight components welded together.The mechanical properties of the material were as follows:ultimate tensile strength 500 MPa,yield strength 360 MPa.Fatigue failures occurred at the welding tip of torsion beam and reinforcing part,see in Fig.7.Service loading histories were recorded on proving ground using strain gauge and stress histories were calculated.Loading spectrum with a sequence length Ls=5951 was obtained after rain-flow counting,and an eight level loading block was derived for simplification of calculation,see in Fig.8.The block load levels and corresponding cycles were listed in Table 3.

    Fig.7.Torsion beam rear axle

    Fig.8.Load spectrum of rear axle

    Table 3.Block loading of rear axle from service load on proving ground

    The S-N curve with Ps=50% was determined by conducting fatigue tests under torsional loads with constant amplitude at a stress ratio of R=-1,and can be expressed as

    From Eq.(9),the fatigue limit σ–1(at N=2×106)was 82.48 MPa.Load with strengthening effect was between 53.62 and 78.36 MPa.So load level 1 to 5 were above fatigue limit;level 6(0.87σ–1)was of strengthening effect;level 7 to 8 can be omitted for fatigue life prediction.

    Taking the strengthening effect of rear axle the same as front axle studied in Ref.[20],for they are of similar mechanical characteristics.According to the curve surface,the maximum increment of fatigue strength under a load at 0.87σ–1was 9.8%,and the best strengthening number Nbwas 200 000 cycles.So the number of repetitions of loading block to finish process I was 454,and the increase in fatigue strength of each block was 0.021 6%.Then the fatigue life prediction of rear axle under service loading history could be expressed as follows:

    So the total number of repetitions of loading block was 501.

    Fatigue tests were carried out on MTS 6DOF road load simulator,as shown in Fig.9.Six specimens from the same batch were tested,and the fatigue lives were 486,510,530,476,519 and 539 repetitions of loading block.The results of predicted and test lives were summarized in Table 4.

    Fig.9.Test set-up of rear axle

    Table 4.Comparisons of fatigue lives and test life of rear axle

    4 Discussion

    The accuracy of the proposed model could be evaluated by comparing the calculated lives with experimental lives.All values for predicted and experimental lives of half shaft and rear axle are summarized and correlations between them are presented in Fig.10.

    For the two examples,predicted lives based on different methods are all conservative compared with experimental lives,and LALS is of the highest accuracy.The other three methods are of poor accuracy.Though PMR is of high accuracy for assessment of half shaft,it fails for rear axle.Comparisons and discussions are made as follows.

    (1)Comparison of LALS with PMR.The fatigue life calculated with PMR uses only the loads above fatigue limit,while those below fatigue limit was ignored.However,according to LALS,loads are of strengthening effect before reaching the best strengthening number of cycles Nband of damaging effect after it.In the application of half shaft,the strengthening load(level 5)first caused an increase in fatigue strength in stage I(208 repetitions of loading block),then causes damage in stage II(906 repetitions).The strengthening and damaging counteracted each other,resulting in an accurate prediction.With regard to the rear axle,the strengthening effect of load level 6 continues for 454 repetitions of loading block,leaving only 47 repetitions in stage II.The damage effect caused by strengthening is much smaller than strengthening effect,resulting in a net increase of fatigue strength.PMR neglects the increase of fatigue strength,leading to a much more conservative prediction.

    (2)Comparison of LALS with MMR.The fatigue lives of half shaft and rear axle predicted according to PMR are carried out by extending the S-N curve with a slope k′=k and k′=2k-1 below the fatigue limit.Due to damage effect of loads below fatigue limit was taken into account,fatigue lives predicted based on MMR are more conservative than PMR,not to mention with LALS.It should be noted that,with regard to the rear axle,the fatigue lives predicted according to MMR using an extending slope k′=k and k′=2k-1 are close;while to the half shaft,they are of great difference.That’s because the damage caused by loads below fatigue in the spectrum of rear axle takes up a small proportion,about 1.7%(k′=2k-1)and 3.8%(k′=k),while in the spectrum of half shaft,it takes up a large proportion,to 10.7%(k′=2k-1)and 27.6%(k′=k).

    5 Conclusions

    (1)A method for the fatigue life prediction under service load was introduced considering LALS based on PMR.In this method,the fatigue process was divided into two stages by the best strengthening number of cycles Nbaccording to the characteristics of low-amplitude-load strengthening(LALS).

    (2)The LALS method was validated by fatigue tests of two vehicle components,half shaft and torsion beam rear axle,under service loadings.The percentage errors of the predicted life with proposed method to test lives of the two examples were-3.78% and-1.76% separately.

    (3)Traditional methods,such as PMR and MMR will lead to more conservative predictions for both the half shaft and rear axle compared with LALS.

    (4)Using the new method,lightweight design can be further achieved in the design stage,because a lower safety factor can be selected due to higher accuracy of fatigue life prediction.

    [1]LEE Y L,PAN J,HATHAWAY R,et al.Fatigue testing and analysis:theory and practice[M].Boston:Butterworth-Heinemann,2005.

    [2]SCHUTZ W.A history of fatigue[J].Engineering Fracture Mechanics,1996,54(2):263–300.

    [3]WIRSCHING P H,PAEZ T L,ORTIZ K.Random vibrations:theory and practice[M].New York:Courier Dover Publications,2006.

    [4]LEE Y L,LU M W,SEGAR R C,et al.Reliability-based cumulative fatigue damage assessment in crack initiation[J].International Journal of Materials and Product Technology,1999,14(1):1–16.

    [5]FRICKE W.Fatigue analysis of welded joints:state of development[J].Marine Structures,2003,16(3):185–200.

    [6]SCHIJVE J.Fatigue of structures and materials[M].New York:Springer,2008.

    [7]HAIBACH E.Modified linear damage accumulation hypothesis accounting for a decreasing fatigue strength during increasing fatigue damage[R].Fraunhofer Institute for Structural Durability-LBF,Darmstadt,1970.

    [8]SCHUTZ W.The prediction of fatigue life in the crack initiation and propagation stages-a state of the art survey[J].Engineering Fracture Mechanics,1979,11(2):405–421.

    [9]HOBBACHER A.Recommendations for fatigue design of welded joints and components[M].New York:Welding Research Council,2009.

    [10]STANZL S E,TSCHEGG E K,MAYER H.Lifetime measurements for random loading in the very high cycle fatigue range[J].International Journal of Fatigue,1986,8(4):195–200.

    [11]SONSINO C M,FRICKE W,DE BRUYNE F,et al.Notch stress concepts for the fatigue assessment of welded joints–background and applications[J].International Journal of Fatigue,2012,34(1):2–16.

    [12]LIN K Y,HWANG J R,CHANG J M.Accelerated durability assessment of motorcycle components in real-time simulation testing[J].Proceedings of the Institution of Mechanical Engineers,Part D:Journal of Automobile Engineering,2010,224(2):245–259.

    [13]FRICKE W,KAHL A.Comparison of different structural stress approaches for fatigue assessment of welded ship structures[J].Marine Structures,2005,18(7):473–488.

    [14]XU Gang,ZHOU Hong,CHEN Donghua.Virtual test rig-based study on fatigue life prediction[J].Journal of Tong Ji University(Natural Science),2009,37(1):97–100.(in Chinese)

    [15]SINCLAIR G M.An investigation of the coaxing effect in fatigue of metals[R].Urbana:Illinois University,1952.

    [16]LAWNICKI J J,CUTLER V C.Effect of under stressing and coaxing on the fatigue limit of a transverse butt-welded joint[J].Journal of Testing and Evaluation,1981,9(1):39–43.

    [17]KAWAMOTO M,NAKAGAWA T,FUJIWARA M.Fatigue strength and hardness of mild steel subjected to the repeated pre-understressing[J].Journal of the Japan Society for Testing Materials,1961,10(91):294–303.(in Japanese)

    [18]IKAI Y,NAKAJIMA M,NAKAGAWA T.Strain aging and fatigue limit on high purity iron[J].Journal of the Japan Society for Testing Materials,1982,31(346):644–649.(in Japanese)

    [19]WANG Z,CHEN Z W.Influence of small load cycle omission on fatigue damage accumulation[C]//Proceedings of the 7th International Fatigue Congress(Fatigue’99),Beijing,China,June 8–12,1999.Beijing:Higher Education Press,1999,2:1113–1118.

    [20]ZHENG Songlin.Studying the effect of low amplitude loading on fatigue life of truck front axle[J].Journal of Mechanical Strength,2002,24(4):547–549.(in Chinese)

    [21]LU Xi,ZHENG Songlin.Strengthening and damaging under low-amplitude loads below the fatigue limit[J].International Journal of Fatigue,2009,31:341–345.

    [22]ZHENG Songlin,YUAN Feng,FENG Jinzhi,et al.Dual-effect experiment study on strengthening and damaging under low amplitude load[J].Journal of Mechanical Strength,2009,31(6):1007–1011.(in Chinese)

    [23]ZHENG Songlin,LU Xi.Microscopic mechanism of strengthening under low-amplitude loads below the fatigue limit[J].Journal of Materials Engineering and Performance,2012,21(7):1526–1533.

    [24]DOWNING S D,SOCIE D F.Simple rainflow counting algorithms[J].International Journal of Fatigue,1982,4(1):31–40.

    [25]SONSINO C M.Fatigue testing under variable amplitude loading[J].International Journal of Fatigue,2007,29(6):1080–1089.

    黄色成人免费大全| 色综合亚洲欧美另类图片| 真人一进一出gif抽搐免费| 国产精品99久久久久久久久| 久久久久久久精品吃奶| 美女大奶头视频| 午夜视频国产福利| 亚洲精品久久国产高清桃花| 18禁黄网站禁片免费观看直播| 脱女人内裤的视频| 成人亚洲精品av一区二区| 欧美性猛交╳xxx乱大交人| 色哟哟哟哟哟哟| 久久人妻av系列| 国产视频内射| 99国产综合亚洲精品| 久久欧美精品欧美久久欧美| 淫妇啪啪啪对白视频| 女人十人毛片免费观看3o分钟| 少妇丰满av| 亚洲电影在线观看av| 久久久精品大字幕| 一级作爱视频免费观看| 香蕉av资源在线| 久久国产精品影院| 真人做人爱边吃奶动态| 国产黄色小视频在线观看| 成人三级黄色视频| 午夜两性在线视频| 亚洲av成人精品一区久久| 久久久国产成人免费| 一夜夜www| 久久99热这里只有精品18| 在线a可以看的网站| 成人欧美大片| 99久久精品国产亚洲精品| 日本免费a在线| 久久99热这里只有精品18| 欧美日韩黄片免| 免费一级毛片在线播放高清视频| 色综合欧美亚洲国产小说| 女人十人毛片免费观看3o分钟| 亚洲国产精品sss在线观看| 香蕉久久夜色| 亚洲国产精品久久男人天堂| 亚洲aⅴ乱码一区二区在线播放| 亚洲专区中文字幕在线| 免费观看精品视频网站| 97超视频在线观看视频| 国产乱人视频| 国产成+人综合+亚洲专区| 亚洲成a人片在线一区二区| 日韩欧美免费精品| 99久久精品热视频| 精品不卡国产一区二区三区| 啦啦啦韩国在线观看视频| 在线观看免费午夜福利视频| 欧美成人a在线观看| 欧美日韩中文字幕国产精品一区二区三区| 人人妻人人澡欧美一区二区| 99久久九九国产精品国产免费| 久久人妻av系列| bbb黄色大片| 亚洲精华国产精华精| 久久精品国产亚洲av香蕉五月| 亚洲电影在线观看av| 中文资源天堂在线| 色综合欧美亚洲国产小说| 久久久色成人| 欧美性感艳星| 丰满人妻一区二区三区视频av | 欧洲精品卡2卡3卡4卡5卡区| 色播亚洲综合网| 亚洲激情在线av| 啦啦啦观看免费观看视频高清| 久久6这里有精品| 美女黄网站色视频| av国产免费在线观看| 欧美日韩黄片免| 9191精品国产免费久久| 黄色丝袜av网址大全| 特大巨黑吊av在线直播| 69av精品久久久久久| 美女黄网站色视频| 欧美日本亚洲视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 国产精品 欧美亚洲| 国产美女午夜福利| 99riav亚洲国产免费| 亚洲性夜色夜夜综合| 色吧在线观看| 村上凉子中文字幕在线| 三级男女做爰猛烈吃奶摸视频| 亚洲成av人片在线播放无| 国产在视频线在精品| 99久久精品国产亚洲精品| 真人一进一出gif抽搐免费| 女人十人毛片免费观看3o分钟| 深夜精品福利| 大型黄色视频在线免费观看| 久久久久国产精品人妻aⅴ院| 美女高潮的动态| 久久亚洲真实| 日本与韩国留学比较| 欧美三级亚洲精品| 亚洲一区二区三区不卡视频| 亚洲不卡免费看| 老司机午夜福利在线观看视频| 国产精品久久久久久人妻精品电影| 欧美黑人欧美精品刺激| 国产精品av视频在线免费观看| 国产高清视频在线播放一区| 免费av不卡在线播放| 国产午夜精品论理片| 亚洲avbb在线观看| 一本一本综合久久| 在线观看美女被高潮喷水网站 | 日本a在线网址| 国产精品,欧美在线| 亚洲乱码一区二区免费版| 精品国产超薄肉色丝袜足j| 久久精品亚洲精品国产色婷小说| 国产探花在线观看一区二区| 久久久久久久亚洲中文字幕 | 国产黄a三级三级三级人| 啦啦啦免费观看视频1| 亚洲专区中文字幕在线| 在线观看免费视频日本深夜| 伊人久久精品亚洲午夜| 免费人成在线观看视频色| e午夜精品久久久久久久| 久久久久九九精品影院| 狠狠狠狠99中文字幕| 久久久久久久午夜电影| 男女做爰动态图高潮gif福利片| 精品国产亚洲在线| 亚洲专区国产一区二区| 久久精品国产综合久久久| 亚洲在线自拍视频| 国产久久久一区二区三区| 亚洲一区高清亚洲精品| 婷婷丁香在线五月| 99精品久久久久人妻精品| 母亲3免费完整高清在线观看| 久久久久久久久久黄片| 男人舔女人下体高潮全视频| 午夜福利在线在线| 老熟妇仑乱视频hdxx| 亚洲av二区三区四区| 亚洲无线观看免费| 美女被艹到高潮喷水动态| 欧洲精品卡2卡3卡4卡5卡区| 少妇人妻精品综合一区二区 | 国产真人三级小视频在线观看| 国产av在哪里看| 亚洲成人免费电影在线观看| 婷婷亚洲欧美| 久久精品国产综合久久久| 制服丝袜大香蕉在线| 色综合亚洲欧美另类图片| 少妇的逼好多水| 亚洲av免费在线观看| 国产精品久久久久久精品电影| 美女免费视频网站| 亚洲 国产 在线| 国产一级毛片七仙女欲春2| 小说图片视频综合网站| 亚洲无线观看免费| 白带黄色成豆腐渣| 一区二区三区激情视频| 99在线视频只有这里精品首页| 18禁裸乳无遮挡免费网站照片| 国产视频内射| 中出人妻视频一区二区| 国产高清videossex| 国产探花在线观看一区二区| 在线观看舔阴道视频| 国产极品精品免费视频能看的| 两个人的视频大全免费| 身体一侧抽搐| 国产一级毛片七仙女欲春2| 波野结衣二区三区在线 | 午夜久久久久精精品| 国产又黄又爽又无遮挡在线| 欧美中文日本在线观看视频| 国产黄色小视频在线观看| 亚洲精品粉嫩美女一区| 久久香蕉国产精品| 久久欧美精品欧美久久欧美| 丁香六月欧美| 欧美3d第一页| 午夜福利成人在线免费观看| 成年免费大片在线观看| 久久草成人影院| 久久久国产成人精品二区| 91麻豆av在线| av片东京热男人的天堂| 久99久视频精品免费| 最近在线观看免费完整版| 午夜精品久久久久久毛片777| 亚洲精品影视一区二区三区av| 最近视频中文字幕2019在线8| 嫩草影院精品99| 长腿黑丝高跟| 成人特级av手机在线观看| 一进一出好大好爽视频| 精品国产三级普通话版| 亚洲精品日韩av片在线观看 | 又紧又爽又黄一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 搡老熟女国产l中国老女人| 少妇的逼水好多| 国产av麻豆久久久久久久| av专区在线播放| 国产av麻豆久久久久久久| 成人精品一区二区免费| 免费观看精品视频网站| 黑人欧美特级aaaaaa片| x7x7x7水蜜桃| 蜜桃亚洲精品一区二区三区| 国产精品 欧美亚洲| 成人精品一区二区免费| 激情在线观看视频在线高清| 久久精品91蜜桃| avwww免费| 国产欧美日韩精品亚洲av| 国产黄色小视频在线观看| 久9热在线精品视频| 性欧美人与动物交配| 午夜视频国产福利| 日韩国内少妇激情av| 欧美区成人在线视频| 午夜免费成人在线视频| 国产精品久久久久久精品电影| АⅤ资源中文在线天堂| 国产亚洲精品久久久com| 老鸭窝网址在线观看| 99riav亚洲国产免费| 亚洲精品影视一区二区三区av| 日韩高清综合在线| 麻豆国产av国片精品| bbb黄色大片| xxxwww97欧美| 国产精品久久电影中文字幕| 真实男女啪啪啪动态图| 国产色婷婷99| 中文字幕高清在线视频| 日韩国内少妇激情av| 国产精品一区二区三区四区免费观看 | 精品久久久久久久人妻蜜臀av| 少妇的逼水好多| 又黄又爽又免费观看的视频| 国语自产精品视频在线第100页| 日韩成人在线观看一区二区三区| 国产乱人视频| 精华霜和精华液先用哪个| 久久精品夜夜夜夜夜久久蜜豆| 美女黄网站色视频| 亚洲精品色激情综合| 天堂√8在线中文| 偷拍熟女少妇极品色| 成年女人永久免费观看视频| 欧美一区二区亚洲| 狠狠狠狠99中文字幕| 俺也久久电影网| 久久久国产成人免费| 两人在一起打扑克的视频| 亚洲人成网站在线播| 精品人妻一区二区三区麻豆 | 欧美极品一区二区三区四区| 女人被狂操c到高潮| 天堂影院成人在线观看| 中文字幕久久专区| 国产亚洲精品av在线| 久久久久国产精品人妻aⅴ院| 在线天堂最新版资源| 久久久久久人人人人人| 在线观看av片永久免费下载| 我的老师免费观看完整版| 亚洲精品在线观看二区| 97人妻精品一区二区三区麻豆| 亚洲成人精品中文字幕电影| 午夜a级毛片| 在线免费观看的www视频| 亚洲真实伦在线观看| 亚洲精品色激情综合| 亚洲,欧美精品.| 露出奶头的视频| 丝袜美腿在线中文| 亚洲成av人片免费观看| 国产69精品久久久久777片| 欧美一级毛片孕妇| 国产成人欧美在线观看| 国产aⅴ精品一区二区三区波| 啦啦啦韩国在线观看视频| 久久精品亚洲精品国产色婷小说| 黄色成人免费大全| 悠悠久久av| 精品国产美女av久久久久小说| 免费高清视频大片| 日本a在线网址| 最后的刺客免费高清国语| 日韩有码中文字幕| 国产v大片淫在线免费观看| 两个人看的免费小视频| 人人妻人人看人人澡| 黄色成人免费大全| 国产爱豆传媒在线观看| 一a级毛片在线观看| 亚洲熟妇中文字幕五十中出| 制服丝袜大香蕉在线| 给我免费播放毛片高清在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲人成电影免费在线| 老司机福利观看| 久久国产精品人妻蜜桃| 精品久久久久久久人妻蜜臀av| 亚洲av五月六月丁香网| 欧美日韩黄片免| 亚洲欧美日韩无卡精品| 色吧在线观看| 国产精品女同一区二区软件 | 国产伦人伦偷精品视频| 午夜福利视频1000在线观看| 美女高潮的动态| 真实男女啪啪啪动态图| 久久国产精品人妻蜜桃| 波多野结衣巨乳人妻| 国语自产精品视频在线第100页| 国产高潮美女av| 久久中文看片网| 最近最新中文字幕大全免费视频| 18禁黄网站禁片午夜丰满| 两个人看的免费小视频| 欧美日韩瑟瑟在线播放| 精品国产超薄肉色丝袜足j| 欧美日本亚洲视频在线播放| 国产蜜桃级精品一区二区三区| a级毛片a级免费在线| 成人午夜高清在线视频| 国产精品美女特级片免费视频播放器| 久久久久精品国产欧美久久久| 成人欧美大片| 麻豆成人av在线观看| av欧美777| 国产精品,欧美在线| 欧美成人a在线观看| 色视频www国产| 亚洲精品在线美女| 久久久久久久精品吃奶| 搡老岳熟女国产| 久久九九热精品免费| 伊人久久精品亚洲午夜| 天堂√8在线中文| www国产在线视频色| 久久久久久九九精品二区国产| 最近最新中文字幕大全免费视频| 久久精品国产综合久久久| 日本精品一区二区三区蜜桃| 欧美日韩亚洲国产一区二区在线观看| 日韩中文字幕欧美一区二区| 美女高潮喷水抽搐中文字幕| 色精品久久人妻99蜜桃| 成人高潮视频无遮挡免费网站| 超碰av人人做人人爽久久 | 免费看日本二区| 成人av一区二区三区在线看| 免费人成视频x8x8入口观看| xxx96com| 午夜福利高清视频| 亚洲国产精品成人综合色| 少妇的逼水好多| 日本精品一区二区三区蜜桃| 首页视频小说图片口味搜索| 国产午夜精品论理片| 精品国产三级普通话版| av视频在线观看入口| 小蜜桃在线观看免费完整版高清| 十八禁网站免费在线| 欧美日韩亚洲国产一区二区在线观看| 身体一侧抽搐| 18禁裸乳无遮挡免费网站照片| 在线国产一区二区在线| 18禁黄网站禁片免费观看直播| 性色avwww在线观看| 亚洲午夜理论影院| 中文字幕av成人在线电影| 在线天堂最新版资源| 亚洲熟妇中文字幕五十中出| 国产淫片久久久久久久久 | 女生性感内裤真人,穿戴方法视频| 丝袜美腿在线中文| 国产v大片淫在线免费观看| 国语自产精品视频在线第100页| 免费电影在线观看免费观看| 一进一出抽搐动态| 波野结衣二区三区在线 | 啪啪无遮挡十八禁网站| 中文字幕人妻熟人妻熟丝袜美 | av视频在线观看入口| 成人18禁在线播放| 亚洲国产高清在线一区二区三| 久久精品影院6| 免费在线观看亚洲国产| 在线十欧美十亚洲十日本专区| www.色视频.com| 好男人在线观看高清免费视频| 色综合婷婷激情| 日本 av在线| 美女高潮的动态| 久久精品91蜜桃| 久久久久久国产a免费观看| 欧美3d第一页| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出| 国产成人欧美在线观看| 亚洲av五月六月丁香网| 亚洲熟妇中文字幕五十中出| 中文字幕精品亚洲无线码一区| www日本在线高清视频| 久久精品国产清高在天天线| 成人无遮挡网站| 欧美最黄视频在线播放免费| 最好的美女福利视频网| 狠狠狠狠99中文字幕| 欧美大码av| 伊人久久大香线蕉亚洲五| 91在线精品国自产拍蜜月 | eeuss影院久久| 欧美成狂野欧美在线观看| 成人精品一区二区免费| 色在线成人网| 男插女下体视频免费在线播放| 国产伦精品一区二区三区四那| 搡老岳熟女国产| 国产亚洲精品久久久com| 亚洲中文日韩欧美视频| 国产主播在线观看一区二区| 欧美色视频一区免费| 国产精品亚洲美女久久久| 国产欧美日韩一区二区三| 操出白浆在线播放| 国产精品98久久久久久宅男小说| 九色成人免费人妻av| 美女高潮的动态| 国产极品精品免费视频能看的| 中亚洲国语对白在线视频| 少妇高潮的动态图| 岛国在线免费视频观看| 国产免费男女视频| 噜噜噜噜噜久久久久久91| 国产色婷婷99| 床上黄色一级片| 成人欧美大片| 国产国拍精品亚洲av在线观看 | 黄色日韩在线| 男女之事视频高清在线观看| 久久久久久久久中文| 亚洲,欧美精品.| www.熟女人妻精品国产| 国产av在哪里看| 亚洲色图av天堂| 国产亚洲精品久久久久久毛片| 欧美性感艳星| 99久久精品一区二区三区| 国内揄拍国产精品人妻在线| 欧美日韩一级在线毛片| 国产黄片美女视频| 天天添夜夜摸| 天堂网av新在线| 亚洲内射少妇av| 国产精品久久久久久久久免 | 国产综合懂色| 久久国产乱子伦精品免费另类| 国产精品嫩草影院av在线观看 | 日韩欧美一区二区三区在线观看| 亚洲成a人片在线一区二区| 亚洲 欧美 日韩 在线 免费| av视频在线观看入口| 可以在线观看毛片的网站| 日韩欧美一区二区三区在线观看| 在线看三级毛片| 色哟哟哟哟哟哟| 精品久久久久久久末码| 国产免费男女视频| 琪琪午夜伦伦电影理论片6080| 日韩欧美国产在线观看| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 精品国产亚洲在线| 亚洲欧美一区二区三区黑人| 免费一级毛片在线播放高清视频| 美女高潮喷水抽搐中文字幕| 亚洲久久久久久中文字幕| 99国产极品粉嫩在线观看| 亚洲成av人片免费观看| 老司机在亚洲福利影院| 香蕉丝袜av| 国产av麻豆久久久久久久| 美女黄网站色视频| 精品久久久久久成人av| 午夜影院日韩av| 在线观看免费午夜福利视频| 午夜两性在线视频| 亚洲国产精品成人综合色| 国产亚洲精品一区二区www| 亚洲精品粉嫩美女一区| svipshipincom国产片| 亚洲精品日韩av片在线观看 | 免费人成在线观看视频色| 18禁国产床啪视频网站| 亚洲人成网站在线播| 亚洲精品亚洲一区二区| 午夜精品在线福利| 少妇高潮的动态图| 亚洲av一区综合| 在线国产一区二区在线| 美女高潮喷水抽搐中文字幕| 日本在线视频免费播放| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 欧美日韩黄片免| 亚洲精华国产精华精| 欧美日韩瑟瑟在线播放| 极品教师在线免费播放| 真人做人爱边吃奶动态| 亚洲av免费在线观看| 蜜桃久久精品国产亚洲av| 男女午夜视频在线观看| 国产成人av教育| 亚洲成人久久性| 毛片女人毛片| 又黄又粗又硬又大视频| 亚洲精品国产精品久久久不卡| 国产成人a区在线观看| 久久精品国产清高在天天线| 蜜桃久久精品国产亚洲av| 少妇人妻一区二区三区视频| av片东京热男人的天堂| 亚洲 国产 在线| 在线观看日韩欧美| 999久久久精品免费观看国产| www.www免费av| av天堂中文字幕网| 亚洲aⅴ乱码一区二区在线播放| 一本综合久久免费| 久久久精品欧美日韩精品| 窝窝影院91人妻| 少妇的逼好多水| 国内精品久久久久精免费| 亚洲av免费在线观看| 观看免费一级毛片| 亚洲激情在线av| 国产在视频线在精品| 黄片小视频在线播放| 久久久精品大字幕| 麻豆国产97在线/欧美| 少妇的逼水好多| 一夜夜www| 国产亚洲精品久久久久久毛片| 成人av在线播放网站| 好男人电影高清在线观看| e午夜精品久久久久久久| 热99在线观看视频| 最新美女视频免费是黄的| www.999成人在线观看| 人人妻人人看人人澡| 午夜亚洲福利在线播放| 小说图片视频综合网站| 99精品久久久久人妻精品| 亚洲美女黄片视频| 国产精品 欧美亚洲| 久久精品影院6| 国产亚洲精品一区二区www| 欧美在线一区亚洲| 亚洲国产高清在线一区二区三| 精品久久久久久久末码| 亚洲成人免费电影在线观看| 亚洲人与动物交配视频| 国产三级在线视频| 午夜福利成人在线免费观看| 亚洲av免费在线观看| 脱女人内裤的视频| 久久99热这里只有精品18| 中文在线观看免费www的网站| 国产私拍福利视频在线观看| 成人欧美大片| 九九热线精品视视频播放| 亚洲无线观看免费| 男女视频在线观看网站免费| 悠悠久久av| 少妇的丰满在线观看| 色噜噜av男人的天堂激情| av黄色大香蕉| 国产私拍福利视频在线观看| 小蜜桃在线观看免费完整版高清| 日本a在线网址| 一个人看视频在线观看www免费 | 内射极品少妇av片p| 色精品久久人妻99蜜桃| 一进一出抽搐gif免费好疼| 国产视频内射| 色噜噜av男人的天堂激情| 99热这里只有是精品50| 午夜福利免费观看在线| 成人永久免费在线观看视频| a在线观看视频网站| 神马国产精品三级电影在线观看| 性欧美人与动物交配| 中文字幕人成人乱码亚洲影| 欧美成人免费av一区二区三区| 美女黄网站色视频| 午夜精品一区二区三区免费看| 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩东京热|