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

    Dynamic simulation and safety evaluation of high-speed trains meeting in open air

    2016-09-06 11:38:17SongyanLiZhijunZhengJilinYuChunqiangQian
    Acta Mechanica Sinica 2016年2期

    Songyan Li·Zhijun Zheng·Jilin Yu·Chunqiang Qian

    ?

    Dynamic simulation and safety evaluation of high-speed trains meeting in open air

    Songyan Li1·Zhijun Zheng1·Jilin Yu1·Chunqiang Qian1

    ?The Chinese Society of Theoretical and Applied Mechanics;Institute of Mechanics,Chinese Academy of Sciences and Springer-Verlag Berlin Heidelberg 2015

    AbstractDynamic responses of a carriage under excitation with the German high-speed low-interference track spectrum together with the air pressure pulse generated as high-speed trains passing each other are investigated with a multi-body dynamics method.The variations of degrees of freedom(DOFs:horizontal movement,roll angle,and yaw angle),thelateralwheel-railforce,thederailmentcoefficient and the rate of wheel load reduction with time when two carriages meet in open air are obtained and compared with the results of a single train travelling at specifie speeds.Results show that the rate of wheel load reduction increases with the increase of train speed and meets some safety standard at a certain speed,but exceeding the value of the rate of wheel load reduction does not necessarily mean derailment. The evaluation standard of the rate of wheel load reduction is somewhat conservative and may be loosened.The pressure pulse has significan effects on the train DOFs,and the evaluations of these safety indexes are strongly suggested in practice.Thepressurepulsehasalimitedeffectonthederailment coefficien and the lateral wheel-rail force,and,thus, their further evaluations may be not necessary.

    KeywordsHigh-speedtrain·Pressurepulse·Derailment· Dynamics simulation·Safety standard

    1 Introduction

    Aerodynamics problems and structural safety involved in modern high-speed trains have gained much attention[1]. When two trains pass each other,an air pressure pulse generates and sweeps through car bodies,which is a type of transientload.Atypicalfeatureofapressurepulseisthatthe value of pressure fluctuatio alternates positively and negatively[2],which is caused by the fore-/after-body of one train passing the other train.The spatial interval between the positive-and the negative-peak pressure values is constant and the peak value of the air pressure pulse is proportional to the square of the relative speed of trains[3].For example, Xiong and Liang[4]measured the pressure pulse generated in the meeting of China Railways CRH2 Electric Multiple Unit(EMU)trains and presented a fittinfunction of the pressure amplitude versus meeting speeds.In addition, Tian and He[5]simulated the three-dimensional flw fiel with a flui dynamics method to obtain the pressure pulse. Furthermore,Cui and Zhang[6]employed flui dynamics and train multibody system dynamics methods,analyzed the dynamicsresponseofatraintravelingunderacrosswind,and evaluated the safety indexes under a crosswind environment. Also,Lai et al.[7]indicated that strong transient loads generated by high-speed trains meeting have significan effects on horizontal movement,roll angle,and yaw angle of trains. When subjected to a pressure pulse,a train is more likely to have snaking motion,overturn,and derail.Lateral wheelrail force[8],derailment coefficien[9],and rate of wheel load reduction[10]are key indexes to evaluate the safety of running vehicle.A track irregularity has a significan effect on the running safety of trains.An investigation by Choi et al.[11]showed the irregularities of cross level,vertical profile and gauge to understand their influenc on the safetyindexes of trains.Wheel-rail unsymmetrical contact caused bywheelprofil wearofhigh-speedtrainsmayleadtolateral instability and a large derailment coefficien[12].In another study,Bocciolone et al.[13]studied the crosswind action on rail vehicles.Meanwhile,Li et al.[14]investigated the characteristics of a pressure wave induced by high-speed trainsmeetinginopenair.AnewmethodfornumericalsimulationwasfoundbyZhaoandSun[15]fortwotrainspassing each other at the same speed.In another study,Liu et al.[16]studied the influenc factors and developed a regularity for an air pressure pulse as trains pass by each other.Also,Qian et al.[17]simulated a dynamic response for side windows of high-speed trains subjected to crossing air pressure pulse.

    Amulti-bodydynamicsmethodisemployedhereintoanalyzethedynamicresponseofamiddlecarofChinaRailways CRH3EMU,whosegeometricandkineticparametersarethe same as the actual vehicle.The air pressure measured from high-speed trains on the Zhengzhou-Xi’an railway passenger dedicated line by China Academy of Railway Sciences is adopted,and excitation with the German high-speed lowinterference track spectrum is applied to the running train model.The vibrations of degrees of freedom(DOFs:horizontal movement,roll angle,and yaw angle),the lateral wheel-rail force,the derailment coefficient and the rate of wheel load reduction of two trains meeting in open air are comparedwiththeresultsofasingletraintravelingatcertain speeds.The influence of air pressure applying to traveling trains on all safety indexes are obtained.

    2 Numerical model

    2.1Multi-body dynamics model

    A China Railways CRH3 EMU has eight carriages.The firs and the last carriages are cabs,and the other carriages are passenger areas.One of the passenger carriages is simplifie andshowninFig.1.Itsdynamicresponsesarecarriedoutby using a multi-body dynamics simulation method in the commercial software SIMPACK.A simplifie model of bogies is shown in Fig.2.The model is established as follows:

    Fig.1 Model of a car body

    Fig.2 Models of bogies

    (1)The models including one car body,two bogies,four wheelsets,and eight rotating arms are assumed to be rigid.Each rotating arm has one DOF,and any of the other free rigid bodies has three translational DOFs and three rotational DOFs.Thus,there are in total 42 independent DOFs and eight constrained DOFs.

    (2)The total mass and the rotational inertias of a carriage, and the parameters of bogies are consistent with those of an actual carriage.

    (3)Some tiny structures,such as door handles and power boxes,are not included in the car body model.

    (4)The air pressure generated by the firs cab is considered, and the peak value of the air pressure is taken to be constant in the vertical direction of a carriage.

    2.2Nonlinear details

    The nonlinear parts of the model mainly include the wheel railcontactnonlineargeometryrelationship,wheelrailcreep, free transverse momentum of the wheel set,and nonlinear suspension system.In this study,the wheel tread of the train is S1002GH and the rail outline is Rail_60H [18].The wheel-rail relationship is shown in Fig.3.The nonlinear factors of the suspension system mainly include the vertical dampers of primary suspension and secondary suspension,thelateraldamperofsecondarysuspension,antisnake movement damper,and the lateral stoppingblock.The force-displacementcurveofalateralstoppingblockisshown in Fig.4.

    2.3Simplificatio and loading of pressure pulse

    Consideratypicalpressurepulsegeneratedbytrainspassing each other in open air applyied to a passenger carriage.Two trainstravelatthesamespeed,denotedas V.Whentheypass each other,the pressure pulse relative to a train moves with speedv,whichistwicethemeetingspeedoftrains,writtenasAnexampleofthepressurepulsemeasuredisgiveninFig.5. The waveform of the central compartment point was measured,when the CRHA380 trains passed each other in open air in the line of Zhengzhou-Xi’an in September,2010. When the heads of the trains meet each other,the amplitude of the pressure wave generated is much higher than that when the tails of the trains meet each other.Here,only the waveform of the pressure pulse when the heads of the trains meet each other is taken into consideration.The simplifie waveform when two trains pass each other at the same speed of V=350 km/h is shown in Fig.6,and the wavelength of the pressure pulse does not vary with the change of the meeting speed.

    Fig.3 Spatial wheel/rail contact geometry relationship

    Fig.4 The force-displacement curve of a lateral stop

    Fig.5 The pressure pulse was measured from the middle point of the passenger carriage by China Academy of Railway Sciences when two CRHA380trainspassedeachotherinopenairontheZhengzhou-Xi’an railway at 350 km/h in September,2010

    Fig.6 Simplifie pressure pulse of CRH3 when two trains pass each other at 350 km/h

    Applyingthepressurepulsesweepingthroughthecarriage to the centroid of the train(the pressure pulse in the vertical direction is uniformly distributed)and using the interaction method,we get the force perpendicular to the side of the carriage Fy,the shaking moment Mz,and the overturning moment Mx.The expressions are given by

    where L is the length of a single carriage,h the effective height of train,D the height difference between the geometric center and the centroid of carriage,P(x)the spatialdistributionfunctionofthepressurepulse,andxcthelocation of the centroid of carriage.The force and moments applying to the carriage when two trains pass each other in open air at 350 km/h are shown in Fig.7.According to the rule that the amplitude of the pressure pulse is in proportion to the square of the meeting speed,we have

    Fig.7 The force and moments of CRH3 when two trains pass each other at 350 km/h

    where p is the absolute value of the difference between the positive and the negative amplitude of the pressure pulse, whose unit is kPa,and V the meeting speed of trains,whose unit is km/h.Furthermore,we can get the force,Fy,and the two moments,Mzand Mx,at different meeting speeds.

    2.4Working conditions

    The rails of railway are not in the ideal straight and fla situation.They have several types of irregularity,including alignmentirregularity,crosslevelirregularity,verticalprofil irregularity,and gauge irregularity.These geometric irregularities make the wheel have much interaction with the rail, cause rail surface abrasion and track geometry changes,and make the rail have random and imbalanced conditions[11]. The faster a train travels,the larger frequency range of the track irregularity effects.This increases the disturbance and decreases the stability of traveling.The international standard of the German high-speed low-disturbance spectrum [11]is applied in the present simulation.For the comparability of results,the track spectrum is used at different traveling speeds of trains and different track irregularities.

    Theworkingconditionsforthemulti-bodydynamicssimulation in this study are:

    (1)Transient pressure pulses are generated under the condition of trains meeting at constant speeds of 250,300, 350,400,450,and 500 km/h.

    Fig.8 Irregularities of the German high-speed low-disturbance spectrum.a Vertical profil irregularity and alignment irregularity.b Cross level irregularity(The data are from the commercial software SIMPACK.)

    (2)The track excitation is the German high-speed lowdisturbance spectrum[19],including vertical profile cross level and alignment irregularities,see Fig.8.

    3 Results and discussion

    3.1Dynamic response of train DOFs

    When a high-speed train moves on irregular tracks,it has small vibration responses.The maximum amplitude that the

    Fig.9 Comparisonofvibrationsbetweentwotrainspassingeachother and a single train traveling at 400 km/h.a Horizontal movement.b Roll angle.c Yaw angle

    vibration responses can reach is an important condition for measuring the operation’s stability.When two trains pass each other,the swings of the trains may be much evident. The variations of horizontal movement,roll angle,and yaw angle when two trains pass each other in open air at 400 km/h,and when a train travels alone are shown in Fig.9. Vibrations in all directions are significantl increased with the increase of time.The maximum values of vibrations are given in Table 1.It is found that the pressure wave when two trains pass each other in open air has significan effects on the vibrations of train DOFs.The horizontal movement and the yaw angle of the meeting trains vary more strongly than those of a single traveling train.The severe vibrations may increase the likelihood of capsizing and serpentine locomotion of trains,and decrease the train operation stability and the passenger comfort.

    Table 1 The maximums of vibration when two trains pass each other and when a single train travels at 400 km/h

    3.2Lateral wheel-rail force

    The lateral wheel-rail force is the transverse component of theforcebetween trackand wheel,whosemaximum value is used to judge whether the train will lead to gauge widening or produce some severe deformation on the line.The firs and the second limit values of the lateral wheel-rail force are given by[8]

    where H is the maximum of the lateral wheel-rail force with theunitofkN,and Pst1and Pst2thestaticloadsoftheleft-and right-side wheels,respectively,with the unit of kN.In general,the lateral wheel-rail force requires meeting a tougher second limit.In the other words,it has certain standards for a safety margin.According to the railway vehiclesspecificatio for evaluation,the dynamic performance,and accreditation test[8],the weight of the locomotive axle load studiedinthisstudyis116.35kN,andthesecondstandardin Eq.(4)gives H≤41.47 kN.The lateral wheel-rail forces of the four wheelsets versus train speeds are shown in Fig.10. It shows that the lateral wheel-rail force of each wheelset is increased with the increase of the meeting speed.When the meeting speed is under 500 km/h,all the lateral wheel-rail forces are less than the critical value 41.47 kN.The lateral wheel-railforceisaffectedbyboththerailrandomexcitation and the air pressure pulse.The lateral wheel-rail force at 250 km/h could be slightly greater than that at 300 km/h in some wheelsets.

    3.3Derailment coefficien

    Derailment means that the geometric relationship between wheel and rail is broken.This possibility of derailment is measured by the derailment coefficient which is define as the ratio of wheel lateral(Qw)to vertical(Pw)forces

    Fig.10 The lateral wheel-rail force versus train speed.a Wheelset 1. b Wheelset 2.c Wheelset 3.d Wheelset 4

    Table 2 Evaluation criteria for the derailment coefficien[8]

    at the flange-railheacontact surface.The Chinese derailment coefficien standard is shown in Table 2.According to the code for design of a high-speed railway(for trial implementation)promulgated in 2009,the value of the derailment coefficien is Qw/Pw≤0.8.

    The vehicle model used in this study has four axles. Each axle has two wheels,and there are eight wheels in total.The calculation of the derailment coefficien is for two working conditions:a single train travels on tracks and two trains pass each other on tracks.The train speeds for the two working conditions are 250,300,350,400,450, and 500 km/h.The derailment coefficientfor the eight wheels at 400 km/h are shown in Table 3.The derailment coefficien of each wheel when the trains pass each other is slightly larger than that of a single traveling train. The derailment coefficien of the fourth right wheel is larger than that of others.The derailment coefficient of the fourth right wheel at different traveling speeds are shown in Fig.11.At a low traveling speed(<400 km/h),the pressure pulse does not have obvious effects on the derailment coefficient At a high traveling speed(>400 km/h),the pressure pulse makes the derailment coefficien obviously increase.For example,at 500 km/h,the derailment coeffi cient increases from 0.248 to 0.289,16.54%higher.It is found that all the derailment coefficient are far lower than 0.6 and,thus,the levels of evaluation are all in the excellent range.

    3.4Rate of wheel load reduction

    The rate of wheel load reduction is another important index to evaluate the safety of a running train.It is define as the ratio of the load reduction of the wheel vertical force?ΔPw to the average of the left and the right wheel vertical forces Pw,i.e.,?ΔPw/Pw.Therateofwheelloadreductionplaysa supporting and complementary role to the derailment coeffi cient.ThewheelloadreductionevaluationstandardsofChina are shown in Table 4.The wheel load reductions at 400 km/h are presented in Table 5.The rates of wheel load reduction for the four wheelsets when traveling at different speeds are showninFig.12.Farbelow400km/h,theratesofwheelload reductionwhenatraintravelsaloneortwotrainspassingeach other in open air are less than 0.65,which is the first-l vel standard of China;above 400 km/h,the values exceed 0.65 and are close to 1 at 500 km/h.

    Table 3 Derailment coefficient of two trains passing each other and a single train traveling at 400 km/h

    Fig.11 The derailment coefficien of wheelset 4 right versus train speed

    Table 4 Evaluation criteria for the rate of wheel load reduction[9]

    Table 5 The rate of wheel load reduction of two trains passing each other and a single train traveling at 400 km/h

    There is a special explanation needed for the contradiction that the rate of wheel load reduction and the derailment coefficien sometimes does appear in the EMU test.Sometimes the derailment coefficien is not too large,but the rate of wheel load reduction could exceed the standard.In the EMU speed test of the“Pioneer”and“China Star”,even if the rate of wheel load reduction standard is relaxed to 0.8, there are some situations where the standard is exceeded. Some European scholars[3]think that the dynamic rate of wheelloadreductionisnotnecessary.SomeNorthAmerican scholars[9]relax the standard to 0.9.The rate of wheel load reduction in the dynamic conditions have a relationship with both magnitude and direction of lateral wheel-rail force.If the lateral wheel-rail force is not conducive to derailment, the derailment coefficien will not be too large and a train will not derail.Japanese scholars[9]even allow the rate of wheel load reduction larger than 1.0,but the duration can not exceed 0.015 s.Therefore,the derailment coefficien is a main standard while the rate of wheel load reduction is an auxiliary standard.If the derailment coefficien and the rate of wheel load reduction both exceed the maximum values allowed,the rate of wheel load reduction is the main factor causing the derailment.If the derailment coefficien is not large but the rate of wheel load reduction is exceeded,we should analyze and examine the cause.

    From the comparison in Fig.11,the rate of wheel load reduction when trains pass each other in open air is larger than that when a single train is traveling.Especially at a high speed(>400 km/h),the increase of the rate of wheel load reduction is much more obvious,for which one should be pay much attention.

    4 Conclusions

    A carriage model under excitation with the German highspeed low-interference track spectrum and action with a pressure pulse is proposed with a multi-body dynamics simulation method in the commercial software SIMPACK.The carriage model is applied with the pressure pulse generated by trains passing each other at different speeds to simulate various conditions.It is focused on understanding the effects of the pressure pulse on the safety of trains meeting in open air.The investigations of the vibrations of DOFs(horizontal movement,roll angle,and yaw angle),the lateral wheel-rail force,the derailment coefficient and the rate of wheel load reduction are carried out.

    (1)When the train speed is not less than 400 km/h,all safety indexes calculated meet the requirements for the two cases of a single train traveling and two trains meeting in open air.When the train speed is higher than 450 km/h, therateofwheel loadreductionexceeds 0.65.According to the safety standard for the rate of wheel load reduction in Europe and United States,the derailment coefficien playsthemostimportantroleinderailment,whiletherate of wheel load reduction is a complementary role to the derailment,so exceeding the value of the rate of wheel load reduction does not necessarily mean derailment.It is suggested that the evaluation standard in China is too conservative,which may be loosened to 1 with limiting the exceeding duration.

    Fig.12 Therateofwheelloadreductionversustrainspeed.aWheelset 1 b Wheelset 2 c Wheelset 3 d Wheelset 4

    (2)The variations of train DOFs(horizontal movement,roll angle,and yaw angle),the lateral wheel-rail force,the derailment coefficient and the rate of wheel load reduction with time when trains pass each other in open air are larger than those when a single train is traveling,which means that,the pressure generated by trains meeting,significantlaffects the dynamic response of high-speed trains.The horizontal movement and the yaw angle apparently change much.Thus,these two indexes shouldbeevaluatedinpractice.Theevaluationoftheroll angle is also necessary.

    (3)The results show that the pressure pulse has a limited effectonthederailmentcoefficien andthelateralwheelrailforce.Eveniftwotrainspasseachotherinopenairat ahighspeedabove400km/h,thevaluesofthederailment coefficient andthelateralwheel-railforcesaremuchless than the safety standard values.The further evaluations of the derailment coefficien and the lateral wheel-rail force may be not necessary in practice.

    AcknowledgmentsThe financia supports from the National Basic Research Program of China(973 program,Grant 2011CB711100),the National Natural Science Foundation of China(Project No.11372307)and the Chinese Academy of Sciences(Grant KJCX2-EW-L03)are acknowledged.

    References

    1.Baker,C.:The flw around high speed trains.J.Wind Eng.Ind. Aerodyn.98,277-298(2010)

    2.Komatsu,N.,Yamada,F.:The Reduction of the Train Draft Pressure in Passing by Each Other.Proceedings of World Congrence Railway Research,Tokyo(1999)

    3.Raghunathan,R.S.,Kim,H.D.,Setoguchi,T.:Aerodynamics of high-speed railway train.Prog.Aerosp.Sci.38,469-514(2002)

    4.Xiong,X.H.,Liang,X.F.:Analysisofairpressurepulsesinmeeting ofCRH2EMUtrains.TiedaoXuebao/J.ChinaRailw.Soc.31,15-20(2009).(in Chinese)

    5.Tian,H.Q.,He,D.X.:3-DNumericalcalculationoftheairpressure pulse from two trains passing by each other.Tiedao Xuebao/J. China Railw.Soc.23,18-22(2001).(in Chinese)

    6.Cui,T.,Zhang,W.H.:A new studying method for safety of a train with attitude changing in side wind.Zhendong yu Chongji/J.Vib. Shock 30,143-146(2011).(in Chinese)

    7.Lai,J.,Zeng,X.H.,Sun,Z.X.,et al.:The dynamic responses of high-speed trains passing tunnel,passing by each other on open trackandinsidetunnel.Proc.Twenty-FirstNatl.Conf.Struct.Eng. II,222-225(2012)

    8.GB/T 5599-1985.Railway vehicles-Specificatiofor evaluation the dynamic performance and accreditation test.(1985)

    9.Zeng,Y.Q.,Wang,W.D.,Shu,X.G.,et al.:Dynamic limits for derailment safety evaluation.J.China Acad.Railw.Sci.20,70-77 (1999)

    10.Zhang,H.,Yang,G.Z.:Derailment and wheel unloading of passenger car bogies.J.Rolling Stock 43,10-16(2005).(in Chinese)

    11.Choi,I.Y.,Um,J.H.,Lee,J.S.,et al.:The influenc of track irregularities on the running behavior of high-speed trains.Proc.Inst. Mech.Eng.Part F J.Rail Rapid Transit 227,94-102(2013)

    12.Liu,W.,Ma,S.H.,Luo,S.H.,et al.:Wheel/Rail unsymmetrical contact of a high-speed train due to wheel wear.Zhendong yu Chongji/J.Vib.Shock 32,107-107(2013).(in Chinese)

    13.Bocciolone,M.,Cheli,F.,Corradi,R.,et al.:Crosswind action on rail vehicles:wind tunnel experimental analyses.J.Wind Eng.Ind. Aerodyn.96,584-610(2008)

    14.Li,X.J.,Liu,C.H.,Yang,G.W.:Characteristics of pressure wave induced by high speed trains meeting in open air.J.Sci.Technol. Eng.12,3689-3694(2012).(in Chinese)

    15.Zhao,X.,Sun,Z.:A new method for numerical simulation of two trains passing by each other at the same speed.J.Hydrodyn.22, 697-702(2010)

    16.Liu,J.,Li,R.X.,Qi,Z.D.:The influencfactors and developing regularity of air pressure pulse as trains passing through each other.Third International Conference on Transportation Engineering(ICTE),2707-2713.Chengdu,China 2011

    17.Qian,C.Q.,Zheng,Z.J.,Yu,J.L.,et al.:Dynamics response of side windows of high-speed trains subjected to crossing air pressure pulse.J.Mech.Eng.49,30-36(2013).(in Chinese)

    18.Railway Industrial Standard TB/T 449-2003:Profile of wheel flang tread for locomotives and cars.Ministry of railway of the People’s Republic of China(2003).(in Chinese)

    19.Zhao,Y.,Zhang,Y.H.,Lin,J.H.:Summaryonthepseudo-excitation method for vehicle random vibration PSD analysis.Appl.Math. Mech.34,107-117(2013).(in Chinese)

    24 September 2014/Revised:28 December 2014/Accepted:3 May 2015/Published online:15 August 2015

    ?Zhijun Zheng

    zjzheng@ustc.edu.cn

    1CAS Key Laboratory of Mechanical Behavior and Design of Materials,University of Science and Technology of China, Hefei 230026,China

    天天一区二区日本电影三级| 国产一区二区三区在线臀色熟女| 男人舔女人下体高潮全视频| 亚洲精品久久成人aⅴ小说| 最好的美女福利视频网| 男女午夜视频在线观看| 欧美在线黄色| 久久亚洲精品不卡| 国内揄拍国产精品人妻在线 | 香蕉av资源在线| 亚洲片人在线观看| 草草在线视频免费看| 国产一卡二卡三卡精品| 久久亚洲真实| 亚洲欧美精品综合久久99| 在线国产一区二区在线| 夜夜看夜夜爽夜夜摸| 国内揄拍国产精品人妻在线 | 久久青草综合色| 国产精品乱码一区二三区的特点| 在线十欧美十亚洲十日本专区| 18禁国产床啪视频网站| 亚洲成人免费电影在线观看| 黑人欧美特级aaaaaa片| 欧美色视频一区免费| 色综合站精品国产| 日韩精品青青久久久久久| 欧美精品啪啪一区二区三区| 正在播放国产对白刺激| 给我免费播放毛片高清在线观看| 视频区欧美日本亚洲| 亚洲欧美日韩无卡精品| 午夜两性在线视频| 国产精品亚洲美女久久久| 黄色丝袜av网址大全| 精品日产1卡2卡| 欧美一级a爱片免费观看看 | 国产又爽黄色视频| 老鸭窝网址在线观看| 成年免费大片在线观看| 天天一区二区日本电影三级| 午夜福利一区二区在线看| 夜夜爽天天搞| www.熟女人妻精品国产| 精品不卡国产一区二区三区| 天堂动漫精品| 亚洲人成伊人成综合网2020| 又黄又粗又硬又大视频| 亚洲精品久久成人aⅴ小说| 啪啪无遮挡十八禁网站| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 国产高清激情床上av| 夜夜爽天天搞| 免费观看精品视频网站| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| 欧美成人免费av一区二区三区| 不卡av一区二区三区| 97碰自拍视频| 波多野结衣高清作品| 亚洲黑人精品在线| 午夜激情福利司机影院| 婷婷精品国产亚洲av在线| 国产欧美日韩一区二区三| 听说在线观看完整版免费高清| 久久香蕉激情| 正在播放国产对白刺激| 99久久99久久久精品蜜桃| 757午夜福利合集在线观看| 国内少妇人妻偷人精品xxx网站 | √禁漫天堂资源中文www| 男男h啪啪无遮挡| 一边摸一边做爽爽视频免费| 精品人妻1区二区| 亚洲精品国产精品久久久不卡| 亚洲精品中文字幕一二三四区| 日本一区二区免费在线视频| 看免费av毛片| 亚洲精品在线美女| 午夜福利一区二区在线看| cao死你这个sao货| 国产视频内射| 欧美一级a爱片免费观看看 | 又黄又粗又硬又大视频| 久久精品91蜜桃| 91麻豆精品激情在线观看国产| 国产一区二区激情短视频| 黄频高清免费视频| 日韩有码中文字幕| 欧美色欧美亚洲另类二区| 午夜久久久在线观看| 91老司机精品| 亚洲自偷自拍图片 自拍| 欧美黄色淫秽网站| 可以在线观看的亚洲视频| 欧美黄色淫秽网站| 国产亚洲精品av在线| 自线自在国产av| 一a级毛片在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟妇熟女久久| 美女高潮到喷水免费观看| 欧美日韩福利视频一区二区| 中文字幕精品亚洲无线码一区 | 午夜福利成人在线免费观看| 黑人操中国人逼视频| ponron亚洲| 一进一出好大好爽视频| 欧美性猛交╳xxx乱大交人| 一级作爱视频免费观看| 黄色a级毛片大全视频| 国产精品98久久久久久宅男小说| av欧美777| 亚洲精品中文字幕一二三四区| 国产真实乱freesex| 日韩欧美国产在线观看| 岛国在线观看网站| 亚洲国产毛片av蜜桃av| 黄色视频,在线免费观看| 99久久国产精品久久久| 精品久久久久久成人av| 一级毛片女人18水好多| 白带黄色成豆腐渣| 日韩欧美国产一区二区入口| 国产私拍福利视频在线观看| 日日干狠狠操夜夜爽| 国产成人av教育| 高潮久久久久久久久久久不卡| 国产黄片美女视频| 波多野结衣av一区二区av| 国产午夜精品久久久久久| 国产精品久久久久久人妻精品电影| 夜夜夜夜夜久久久久| 久久青草综合色| 日本黄色视频三级网站网址| 一区二区日韩欧美中文字幕| 亚洲精品中文字幕一二三四区| 禁无遮挡网站| 亚洲男人的天堂狠狠| 欧美乱码精品一区二区三区| 精品不卡国产一区二区三区| 免费在线观看日本一区| 久久伊人香网站| 成年版毛片免费区| 国产爱豆传媒在线观看 | 日韩免费av在线播放| videosex国产| 好男人在线观看高清免费视频 | 久久欧美精品欧美久久欧美| 美女大奶头视频| 中文字幕另类日韩欧美亚洲嫩草| 一夜夜www| av在线播放免费不卡| 亚洲av片天天在线观看| 国产精品1区2区在线观看.| 国产男靠女视频免费网站| 男人舔女人下体高潮全视频| 免费搜索国产男女视频| 1024视频免费在线观看| 美女免费视频网站| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| 国产v大片淫在线免费观看| 久久久久免费精品人妻一区二区 | 亚洲国产精品999在线| 视频区欧美日本亚洲| 一本一本综合久久| 欧美大码av| 黄色成人免费大全| 巨乳人妻的诱惑在线观看| 午夜亚洲福利在线播放| 国产区一区二久久| 一本综合久久免费| 91成人精品电影| 成在线人永久免费视频| 国产区一区二久久| 亚洲国产欧美网| 国产爱豆传媒在线观看 | 免费高清在线观看日韩| 男人操女人黄网站| 精品第一国产精品| 亚洲专区国产一区二区| 天堂动漫精品| 18禁黄网站禁片午夜丰满| 国产成人精品久久二区二区免费| 午夜福利在线在线| 黄频高清免费视频| 国产亚洲欧美精品永久| 国产精品98久久久久久宅男小说| 美女大奶头视频| 亚洲九九香蕉| 精品一区二区三区四区五区乱码| 啦啦啦观看免费观看视频高清| 欧美成人午夜精品| 男女做爰动态图高潮gif福利片| 欧美成人性av电影在线观看| 少妇的丰满在线观看| 久久精品国产综合久久久| 好男人电影高清在线观看| 中文字幕人妻丝袜一区二区| 香蕉国产在线看| xxx96com| 日本黄色视频三级网站网址| 岛国视频午夜一区免费看| 大香蕉久久成人网| 高潮久久久久久久久久久不卡| 此物有八面人人有两片| 久久九九热精品免费| 国产成+人综合+亚洲专区| 香蕉av资源在线| 1024视频免费在线观看| 一本大道久久a久久精品| 岛国视频午夜一区免费看| 性欧美人与动物交配| 久久精品91蜜桃| 免费女性裸体啪啪无遮挡网站| www.熟女人妻精品国产| 成人三级黄色视频| 99久久国产精品久久久| 精品国产亚洲在线| 日本免费a在线| 国产又黄又爽又无遮挡在线| 欧美在线黄色| av片东京热男人的天堂| 国产成人欧美| 69av精品久久久久久| 日日爽夜夜爽网站| 亚洲成人国产一区在线观看| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| cao死你这个sao货| 国产免费男女视频| 90打野战视频偷拍视频| 一级毛片女人18水好多| 亚洲成人国产一区在线观看| 欧美中文综合在线视频| 无遮挡黄片免费观看| 成人手机av| 在线十欧美十亚洲十日本专区| 国产欧美日韩精品亚洲av| 欧美成人免费av一区二区三区| 日本免费一区二区三区高清不卡| 18禁裸乳无遮挡免费网站照片 | 日韩三级视频一区二区三区| 日韩欧美国产在线观看| 亚洲精品久久成人aⅴ小说| 99热6这里只有精品| 在线观看免费视频日本深夜| 国产高清有码在线观看视频 | 搞女人的毛片| or卡值多少钱| 国产精品电影一区二区三区| 日韩欧美三级三区| 久久久久精品国产欧美久久久| 国产男靠女视频免费网站| 午夜福利视频1000在线观看| 又黄又粗又硬又大视频| 精华霜和精华液先用哪个| 99国产精品一区二区蜜桃av| 久久中文看片网| 久久精品影院6| 最近在线观看免费完整版| 人妻丰满熟妇av一区二区三区| 麻豆久久精品国产亚洲av| videosex国产| 国产单亲对白刺激| 老熟妇乱子伦视频在线观看| 1024香蕉在线观看| 最近在线观看免费完整版| 午夜免费鲁丝| 久久久久久免费高清国产稀缺| 99在线人妻在线中文字幕| xxxwww97欧美| 女性被躁到高潮视频| 日韩欧美一区视频在线观看| 欧美zozozo另类| 国产aⅴ精品一区二区三区波| 欧美亚洲日本最大视频资源| 国产一区二区三区在线臀色熟女| 亚洲avbb在线观看| 国产精品久久久av美女十八| 免费在线观看成人毛片| 精品电影一区二区在线| 一本一本综合久久| 亚洲精品在线美女| 久久久久国内视频| 国产精品99久久99久久久不卡| 免费看美女性在线毛片视频| 亚洲av中文字字幕乱码综合 | 午夜成年电影在线免费观看| 男女下面进入的视频免费午夜 | 日韩免费av在线播放| 大香蕉久久成人网| 一级黄色大片毛片| 免费搜索国产男女视频| 天天一区二区日本电影三级| 男女之事视频高清在线观看| 欧美 亚洲 国产 日韩一| 国产一区二区在线av高清观看| 正在播放国产对白刺激| 丰满的人妻完整版| svipshipincom国产片| 国产av在哪里看| 麻豆成人午夜福利视频| 亚洲成人久久性| 丝袜人妻中文字幕| 手机成人av网站| 91av网站免费观看| 欧美绝顶高潮抽搐喷水| 少妇裸体淫交视频免费看高清 | 久热爱精品视频在线9| 满18在线观看网站| 欧美日韩瑟瑟在线播放| 制服丝袜大香蕉在线| 亚洲国产欧美网| 国产亚洲欧美在线一区二区| 国产激情偷乱视频一区二区| av有码第一页| 一级毛片高清免费大全| 女生性感内裤真人,穿戴方法视频| 黄片播放在线免费| 亚洲熟妇中文字幕五十中出| 亚洲在线自拍视频| 国产av不卡久久| 精品久久久久久久久久免费视频| 国产精品免费一区二区三区在线| 日韩欧美国产在线观看| 亚洲第一青青草原| 成人18禁高潮啪啪吃奶动态图| 婷婷六月久久综合丁香| 首页视频小说图片口味搜索| 亚洲精品一卡2卡三卡4卡5卡| 女性被躁到高潮视频| 国产精品免费视频内射| 国产99白浆流出| 久久中文字幕人妻熟女| 欧美日韩亚洲国产一区二区在线观看| 夜夜看夜夜爽夜夜摸| 亚洲性夜色夜夜综合| 国产人伦9x9x在线观看| 男人舔奶头视频| 日韩高清综合在线| 欧美乱色亚洲激情| 免费av毛片视频| 成年人黄色毛片网站| 日日夜夜操网爽| 美女免费视频网站| 国产熟女午夜一区二区三区| 老司机深夜福利视频在线观看| av在线播放免费不卡| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 婷婷精品国产亚洲av在线| 久久亚洲真实| 看黄色毛片网站| 欧美在线黄色| 中文字幕人妻熟女乱码| 国语自产精品视频在线第100页| 欧美日韩中文字幕国产精品一区二区三区| 久久草成人影院| 麻豆成人av在线观看| 久久中文字幕一级| 国产又黄又爽又无遮挡在线| 女人高潮潮喷娇喘18禁视频| 久久午夜亚洲精品久久| 精品卡一卡二卡四卡免费| 岛国视频午夜一区免费看| 成人精品一区二区免费| 国产又爽黄色视频| 给我免费播放毛片高清在线观看| 欧美绝顶高潮抽搐喷水| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| www日本黄色视频网| 免费在线观看成人毛片| 亚洲人成网站在线播放欧美日韩| 久久中文字幕一级| 啦啦啦韩国在线观看视频| 啦啦啦免费观看视频1| 国产麻豆成人av免费视频| 色老头精品视频在线观看| 亚洲熟妇熟女久久| 99热这里只有精品一区 | 一区福利在线观看| 无人区码免费观看不卡| 欧美成人午夜精品| 男女之事视频高清在线观看| 欧美最黄视频在线播放免费| 草草在线视频免费看| 无人区码免费观看不卡| 亚洲成人免费电影在线观看| 午夜亚洲福利在线播放| 一级毛片精品| 国产一卡二卡三卡精品| 精品久久蜜臀av无| 国产在线观看jvid| 欧美日韩亚洲国产一区二区在线观看| 一进一出好大好爽视频| 一级毛片女人18水好多| or卡值多少钱| 亚洲最大成人中文| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播放欧美日韩| 99在线视频只有这里精品首页| 日本免费一区二区三区高清不卡| e午夜精品久久久久久久| 91成人精品电影| 精品欧美一区二区三区在线| 黄色女人牲交| 高潮久久久久久久久久久不卡| 亚洲精品粉嫩美女一区| 成人精品一区二区免费| 18禁美女被吸乳视频| 757午夜福利合集在线观看| 欧美黑人巨大hd| 女同久久另类99精品国产91| 午夜精品久久久久久毛片777| 18禁观看日本| 在线av久久热| 久久99热这里只有精品18| 精品久久久久久久久久久久久 | 久久九九热精品免费| 黑人巨大精品欧美一区二区mp4| 国产精品98久久久久久宅男小说| 久久久水蜜桃国产精品网| 亚洲男人的天堂狠狠| 波多野结衣高清无吗| 国产色视频综合| 婷婷精品国产亚洲av| 黄频高清免费视频| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 国产精品久久视频播放| 精品日产1卡2卡| 91字幕亚洲| 免费搜索国产男女视频| 久久热在线av| 午夜福利欧美成人| 99久久精品国产亚洲精品| 免费在线观看影片大全网站| 母亲3免费完整高清在线观看| 欧美性长视频在线观看| 国产精品精品国产色婷婷| 99精品在免费线老司机午夜| 免费电影在线观看免费观看| 午夜久久久在线观看| 在线视频色国产色| 国产精品永久免费网站| 日韩三级视频一区二区三区| 又黄又粗又硬又大视频| www.www免费av| 久久国产精品影院| 国产熟女午夜一区二区三区| 村上凉子中文字幕在线| 精品久久久久久久久久久久久 | 丁香欧美五月| 窝窝影院91人妻| 午夜视频精品福利| 69av精品久久久久久| 精品国产乱子伦一区二区三区| 午夜福利在线在线| 欧美日韩一级在线毛片| 精品一区二区三区视频在线观看免费| 禁无遮挡网站| 亚洲熟女毛片儿| 日韩大码丰满熟妇| 亚洲欧美精品综合一区二区三区| 久久久久久久午夜电影| 亚洲第一电影网av| 国产精品久久久av美女十八| 日韩有码中文字幕| 91在线观看av| 99精品欧美一区二区三区四区| 特大巨黑吊av在线直播 | av在线天堂中文字幕| 久久 成人 亚洲| 19禁男女啪啪无遮挡网站| 人人澡人人妻人| 国产私拍福利视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 色在线成人网| e午夜精品久久久久久久| 国产三级在线视频| 国产97色在线日韩免费| 亚洲美女黄片视频| 美女大奶头视频| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| 中国美女看黄片| 变态另类成人亚洲欧美熟女| 人成视频在线观看免费观看| 免费无遮挡裸体视频| a级毛片在线看网站| 变态另类成人亚洲欧美熟女| 一区二区三区精品91| 国产一级毛片七仙女欲春2 | 热re99久久国产66热| 琪琪午夜伦伦电影理论片6080| 国产激情偷乱视频一区二区| 亚洲精品国产精品久久久不卡| 亚洲精品av麻豆狂野| 亚洲熟妇中文字幕五十中出| 男女那种视频在线观看| 国产精品久久久久久精品电影 | 在线观看一区二区三区| xxx96com| 少妇裸体淫交视频免费看高清 | 一区二区三区高清视频在线| 老司机福利观看| 国产精品一区二区三区四区久久 | 国产三级黄色录像| 白带黄色成豆腐渣| 欧美日韩一级在线毛片| 国产真人三级小视频在线观看| 两个人看的免费小视频| 久久中文字幕人妻熟女| 国产亚洲精品久久久久久毛片| bbb黄色大片| 男女下面进入的视频免费午夜 | 久久精品国产清高在天天线| 在线国产一区二区在线| 最近最新免费中文字幕在线| 欧美+亚洲+日韩+国产| 色综合欧美亚洲国产小说| 99精品欧美一区二区三区四区| 香蕉久久夜色| 男女视频在线观看网站免费 | 美女午夜性视频免费| 啦啦啦观看免费观看视频高清| 日本a在线网址| 国产99白浆流出| 可以免费在线观看a视频的电影网站| 欧美激情 高清一区二区三区| 在线国产一区二区在线| 丝袜人妻中文字幕| 2021天堂中文幕一二区在线观 | 听说在线观看完整版免费高清| 一级a爱视频在线免费观看| 国产国语露脸激情在线看| 久热这里只有精品99| 国产成人精品久久二区二区免费| 亚洲精品久久国产高清桃花| 人人妻人人澡欧美一区二区| 日本一本二区三区精品| 精品国产乱码久久久久久男人| 国产精品亚洲一级av第二区| aaaaa片日本免费| 亚洲第一青青草原| 中文资源天堂在线| 嫁个100分男人电影在线观看| 亚洲自拍偷在线| 亚洲天堂国产精品一区在线| 俄罗斯特黄特色一大片| 欧美日韩黄片免| 成熟少妇高潮喷水视频| 女性生殖器流出的白浆| 精品少妇一区二区三区视频日本电影| e午夜精品久久久久久久| 美女高潮喷水抽搐中文字幕| 欧美 亚洲 国产 日韩一| 久9热在线精品视频| av欧美777| 色综合婷婷激情| 一本久久中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 一级毛片女人18水好多| 母亲3免费完整高清在线观看| 日韩精品青青久久久久久| 午夜久久久久精精品| 中文字幕另类日韩欧美亚洲嫩草| 麻豆久久精品国产亚洲av| 国产激情久久老熟女| 久久久久久久午夜电影| 久久香蕉精品热| 非洲黑人性xxxx精品又粗又长| 亚洲成人国产一区在线观看| 99国产精品一区二区三区| av福利片在线| 啦啦啦免费观看视频1| 日日爽夜夜爽网站| 久久久久国内视频| 美女大奶头视频| 丰满人妻熟妇乱又伦精品不卡| 免费高清视频大片| 色av中文字幕| 亚洲国产欧美网| www.999成人在线观看| 成年人黄色毛片网站| 亚洲人成网站在线播放欧美日韩| 中文字幕精品亚洲无线码一区 | 亚洲aⅴ乱码一区二区在线播放 | 伊人久久大香线蕉亚洲五| 在线看三级毛片| 此物有八面人人有两片| 成人亚洲精品av一区二区| 亚洲精品久久成人aⅴ小说| 黄片播放在线免费| 亚洲av电影在线进入| 成人国产综合亚洲| 给我免费播放毛片高清在线观看| 欧美又色又爽又黄视频| 国产精品久久电影中文字幕| 国产精品一区二区三区四区久久 | 久久久久亚洲av毛片大全| 狂野欧美激情性xxxx| 一进一出好大好爽视频| 真人做人爱边吃奶动态| АⅤ资源中文在线天堂| 国产片内射在线| 亚洲精品色激情综合| 色播在线永久视频| 国产伦人伦偷精品视频| 成人精品一区二区免费| 国产精品自产拍在线观看55亚洲|