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

    Cam profile optimization design of variable cycle reciprocating piston engine*

    2014-03-09 02:05:38LimingDIWanliSUNYingSHIWenfengWANG
    機床與液壓 2014年6期
    關鍵詞:活塞式配氣型線

    Li-ming DI,Wan-li SUN,Ying SHI,Wen-feng WANG,,,;

    2Shenyang Jiashi Judicial Identification,Shenyang 110023,China;

    3Research Institute of Highway Ministry of Transport,Beijing 100088,China

    Cam profile optimization design of variable cycle reciprocating piston engine*

    Li-ming DI?1,Wan-li SUN2,Ying SHI3,Wen-feng WANG11School of Vehicle and Energy Sources,Yanshan University,Qinhuangdao 066004,China;

    2Shenyang Jiashi Judicial Identification,Shenyang 110023,China;

    3Research Institute of Highway Ministry of Transport,Beijing 100088,China

    The kinematics and dynamics simulation model for valve gear of HONDA 125 mL single-cylinder four-stroke gasoline engine was studied with the software of Excite_TD,by analysing and comparing the performance indexes of the cam profile design methods between polydynamic cam(POL)and individual section acceleration cam(ISAC).ISAC was chosen to optimize the intake and exhaust cam profile for the variable cycle reciprocating piston engine(VCRPE)in fourstroke or two-stroke mode.In order to validate the rationality of the cam profile design,the multibody dynamics simulation model of the corresponding valve gear was established and analyzed by SolidWorks.Through comparison of the dynamics simulation results between Excite_TD and SolidWorks,it shows that the results from Excite_TD are consistent with that from SolidWorks,which verifies the accuracy of cam profile optimization as well.

    Reciprocating piston engine, Variable cycle, Valve gear, Cam profile,Optimization design

    *Project supported by Hebei Education Department in China Natural and Science Research of Hebei Education Department in China(No.2009160),Doctor Key Foundation of Yanshan University(2013925)and Corporation Unit Exploitation Project of Yanshan University(2013101)

    ?Li-ming DI,Associate professor.E-mail:diliming@ysu.edu.cn

    1.Introduction

    Variable cycle reciprocating piston engine (VCRPE)combines the two-stroke and four-stroke mode together.For improving engine load factor,VCRPE can switch between the two-stroke and fourstroke mode timely according to engine operating condition,which provides a new method for engine miniaturization and lightweight design[1-4].Valve train directly affects the engine performance,and good valve gear can ensure engine’s smooth operation, less vibration and noise,and good life cycle.Cam is an important part of the valve gear,and its profile has a direct impact on valve motion and valve timing.Reasonable organization with the valve timing will ensure the gas exchange thoroughly and reduce the exchange losses.

    VCRPE was realized by transformation of fourstroke engine with existing air supply system,fuel system,cooling system,lubrication system and exhaust purification system[3-4].In the VCRPE working process,four-stroke mode and two-stroke mode can be realized by controlling the intake and exhaust valve opening and closing time.Thus,reasonable design of the cam profile has an important influence on the performance of the VCRPE.

    Cam design includes ramp section and working section.Ramp section has three design methods,which are rectangle function,cosine function and trapezoid function,and includes three parameters that are ramp height at the valve,total length of ramp and velocity at section end.Rationalization with the design methods and parameters ensure the smooth connection between ramp and working section,and avoid the speed,acceleration and the contact force exceeding the allowable limits.Working section also has three design methods,which are constant acceleration cam(STAC),polydynamic cam(POL)and individual section acceleration cam(ISAC),and each method has its application field.STAC is used for large and medium-sized low-speed engines,and POL is mostly applied for the tappet valve gear.STAC,which can freely mix arbitrary function and deposit flexibly,is applied to any kind of valve gear,and it’s the preferred method in cam optimization[5-10].

    The research object is HONDA 125 mL singlecylinder four-stroke gasoline engine in this paper.According to the design requirements of VCRPE valve gear,the professional software Excite_TD of AVL Company was used to establish the kinematics and dynamics simulation model.Then POL and STAC method were adopted to design the cam profile.By analyzing and comparing the performance index of the cam profile,the best design method was finally chosen to optimize the intake and exhaust cam profile for VCRPE in four-stroke or two-stroke mode.In the case of lacking in experimental conditions,the multibody dynamics simulation model of the corresponding valve gear was established and analyzed by Solid-Works,results from Solidworks multi-body dynamics simulation were used to verify the accuracy of the Excite_TD dynamics simulation.

    2.Intake cam profile design offourstroke mode

    The various parameters required for POL design method were reasonable selected by the orthogonal test in reference[11],adopted POL symmetric method and designed the intake and exhaust cam profile of four-stroke mode of VCRPE,the results basically meet the design requirements.However,the negative radius of curvature appears in the kinematics evaluation,which increases the difficulty of cam manufacture.If the minimum radius of curvature is smaller than the radius of cam follower,the follower couldn’t contact with the cam profile.On the basis of kinematics and dynamic simulation model in reference[11],the cam profiles were optimized by Excite_TD in this paper.

    2.1.Establishment of simulation models

    Cam design speed was at 4 000 r/min which is the maximum torque crankshaft speed of four-stroke mode,and the valve timing when valve lift at 1 mm is shown in Table 1.In ramp section design,cosine ramp was selected.Ramp height aims to eliminate the valve clearance and elastic deformation of the valve gear,in this paper it was set at 0.15 mm.Velocity at section end not only affects the shape of ramp,but also influences characteristic of the working section[12],this paper chose 250 mm/deg.

    Table 1.Valve timing for valve lift at 1 mm

    The Excite_TD kinematics and dynamics simulation models of four-stroke mode valve gear were shown in Figure 1 and Figure 2.

    Figure 1.Kinematics simulation model of valve gear

    Figure 2.Dynamics simulation model of valve gear

    To compare the advantages and disadvantages of different cam design methods,two cam profile design plans are set up.Plan 1 was based on the asymmetric design,the opening section used POL method while the closing section used the STAC method.Plan 2 was also based on the asymmetric design,the opening section used STAC method and the closing section used STAC method too.

    2.2.POL and STAC cam profile design

    1)POL cam profile design

    The mathematical expression for working section of POL cam profile design method[13]is:

    whereHmaxis the cam maximum lift;hris the ramp height;x=1-φ/Φ,and Φ is the half main cam length which is 65 deg in this paper;φ is the cam angle,the starting point of cam angle is defined as φ=0;constantc4is mainly used to control the shape of negative acceleration and its generally recommended value is 0.1,0.2 and 0.4.This paper chose 0.1;constantscp,cq,cr,csare decided by the boundary conditions,positive integer exponentsp,q,r,sare greater than or equal 8 and at the same timep<q<r<s;cam profile curve junction need to meet three conditions:the lift curve in junction is continuous,the velocity curve is consecutive and the velocity is zero,and the acceleration curve is continuous and reaches the negative maximum number.Exponent parameters chosep=12,q=24,r=32,s=46 in this paper.

    2)STAC cam profile design

    STAC method has three functions to define the cam profile which includes linear function,polynomial and sine function.STAC connects different sections of cam profile by circular arc to improve the stability of the valve gear,and ensure that the fullness ratio of the cam profile meets the design requirements,it’s a kind of flexible and effective cam profile design method[14].In cam profile design process,the design curve type is mainly determined according to the speed,acceleration,jerk,torque and complexity of the curve matched[15-16].

    The closing section of plan 1,opening and closing section of plan 2 were all adopted of STAC method,and their parameters are shown in Table 2 and Table 3.

    Table 2.Closing section of plan 1

    Table 3.Parameters of plan 2

    2.3.Kinematics and dynamics simulation analysis

    1)Kinematics simulation analysis

    Kinematics evaluation rules include valve lift,velocity,acceleration,radius of the curvature,fullness ratio of the valve lift,the maximum jerk and cam contact force.Cam design should ensure that all the parameters meet design requirements,and then carrying preliminary kinematics analysis on the valve gear.

    Comparing Figure 3 and Figure 4 in kinematics analysis,the valve lift and valve timing both meet the design requirements in two different plans.Curves of valve lift and speed are smooth and continuous,which indicate that the valve gear operates Smoothly.The peak valve acceleration is in a reasonable range.The acceleration curve is smooth and continuous and there is no large mutation,indicating that the valve gear operates well in this speed.

    Figure 3.Valve lift,velocity and acceleration of plan 1

    Figure 4.Valve lift,velocity and acceleration of plan 2

    Other kinematics parameters are shown in Table 4,and all of the parameters meet the design requirements.In cam design,fullness ratio is used to evaluate the time-area value.Bigger of the fullness ratio,better of the charging efficiency,and generally it approximately takes about 0.55.The maximum jerk is usually used to evaluate the stability of valve gear,and it is generally defined less than or equal to 1 000 mm/rad3.The materials of cam and follower in this valve gear are cold cast iron and steel,the maximum permissible contact force between these two materials is 1 250 N/mm2.Coefficient κ is used to describe the relationship between length of positive acceleration section and the vibration.in general,κ is required to greater than 1.2~1.3[12],if κ does not meet the design requirements,it easily lead to resonance.Length of positive acceleration section directly influences the minimum radius of curvature,however,there is no specification on the length of positive acceleration section.As shown in Figure 5,radius of curvature in plan 1 appears to be negative at the opening section,the reason is that,valve timing at 1 mm valve lift is provided which means cam profile length is fixed,and this problem can not be solved by increasing the cam length.STAC method which used in plan 2 can subjective control the cam profile length of each section and avoid the negative radius of curvature.

    Table 4.The kinematic parameters

    2)Dynamics Simulation Analysis

    Comparing the dynamic simulation results of the two plans in the engine’s maximum power speed of 7 000 r/min,the maximum valve lift both is 5.4 mm which is consistent with the design value(Figure 6).Considering of the elastic deformation and impact in dynamic model,the velocity curve has slight fluctuations without large mutations(Figure 7).Both of the valve velocity curves are consistent.The valve-seating velocity is less than 0.7 m/s,meeting the design requirements.The fluctuation of acceleration curve at the opening section of plan 1 is bigger than that of plan 2(Figure 8).In the comparison of cam contact force,the contact force of plan 1 has appeared 0 N in 7 000 r/min crankshaft speed,indicating that the valve gear has fly off phenomenon.However,there is no such phenomenon in plan 2.The comparison of contact force between plan 1 and plan 2 is shown in Figure 9.

    Figure 5.Radius of curvature for plan 1

    Figure 6.Comparison of valve lift

    Figure 7.Comparison of valve velocity

    Figure 8.Comparison of valve acceleration

    Figure 9.Comparison of the cam contact force

    According to the above analysis,STAC design method can effectively solve the negative curvature radius and obviously improves the dynamic performance in high speed.Therefore,STAC method should be adopted in opening and closing section design of cam profile.

    2.4.Comparison on dynamics simulation results of Excite_TD and SolidWorks

    Limited by the experimental conditions,the dynamic test of valve gear cannot be carried out.For testing the accuracy of dynamic simulation results,SolidWorks was used to establish the multi-body dynamics simulation model.Cam profile needed by SolidWorks was exported from Excite_TD,and dynamicsperformance was calculated in different speed.The dynamics simulation results of Excite_TD and SolidWorks were analyzed separately in the valve lift,valve speed,acceleration and cam contact force.Consistency and tendency of the curve would be compared to evaluate the dynamic performance of the valve gear.

    These two kinds of software use different simulation calculation methods,different simulation environment and different processing method with the model related components,so they are independent.Due to lack of experiment condition,using general multi-body dynamics software SolidWorks analysis results to verify the accuracy of the professional software Excite_TD simulation results has certain feasibility.

    Take the comparison of analysis results in the maximum torque crankshaft speed of 4 000 r/min as example,the various dynamics performances would be elaborated separately,as shown in Figures 10~13.

    Figure 10.Valve lift in 4 000 r/min

    Figure 11.Valve velocity in 4 000 r/min

    Figure 12.Valve acceleration in 4 000 r/min

    Figure 13.Cam contact force in 4 000 r/min

    From above figures we can know that,in 4 000 r/min crankshaft speed,Excite_TD’s valve lift is consistent with SolidWorks’s,and the curves are smooth and continuous without bounce,showing that the valve gear is well operated.So does the velocity curve.In the comparison of acceleration,the fluctuation obtained from Excite_TD is bigger than that from SolidWorks.Because Excite_TD is one-dimensional dynamic simulation software,and SolidWorksis three-dimensional dynamic simulation software,they are different in the algorithm,environment and accuracy of the simulation which lead to the inconformity in simulation results.However,the tendency of the two software simulation results is still consistent with each other.

    3.Exhaust cam profile design of fourstroke mode

    According to the experience of intake cam design,the exhaust cam of four-stroke mode adopts the same method with the intake cam,meaning that both the opening and closing section use STAC method.The concrete parameters are shown in Table 5.

    Unlike the intake cam,cylinder pressure and exhaust back pressure should be taken into account in exhaust cam design.Data of cylinder pressure and exhaust back pressure are obtained by performance simulation with software of BOOST.

    Table 5.Parameters for exhaust CAM

    3.1.Kinematics and dynamics evaluation of exhaust cam

    The simulation calculation process is consistent with the intake cam,thus there is no need to repeat.The simulation speed still chose the maximum torque crankshaft speed of~4 000 r/min,and import the cam profile,cylinder pressure and exhaust back pressure into the model to simulate the kinematics and dynamics performance.

    Kinematics evaluation results are shown in Figure 14 and Table 6.All of the parameters meet the design requirements and the valve gear operates stably in 4 000 r/min.

    Figure 14.Valve lift,velocity and acceleration for four-stroke exhaust cam

    Table 6.The kinematic parameters

    Dynamic evaluation results are shown in Table 7,and all of the parameters meet the dynamic design requirements in this speed.

    Table 7.The dynamics parameters

    3.2.Dynamics validation of exhaust cam

    SolidWorks is also used to simulate the multibody dynamics performance in exhaust cam design.The comparison of simulation results between Excite_ TD and SolidWorks show that,the valve lift curve coincides with each other in all speed,the curve is smooth and continuous without significant inverse and the maximum valve lift is 5.23 mm,complied with the design requirements.The valve velocity enhances with the speed increasing gradually,but the simulation results are kept in the reasonable range,the valve seating velocity is less than 0.7 m/s.With the increase of speed,valve acceleration is gradually increasing and the fluctuations at the opening and clos-ing section are still exist in Excite_TD simulation results,which in line with the actual situation.Because of the diversity in material setting and simulation environment,and the contact force is different too.SolidWorks responses slower than Excite_TD in low-speed,but the maximum force is consistent with that of the Excite_TD.As the speed increasing,the response condition of SolidWorks contact force is improved and the tendency is consistent with the tendency of Excite_TD simulation results.The contact force is less than six times of spring preload force and without 0 N contact force situation within the work section.Dynamics validation between Excite_TD and SolidWorks indicates that,the higher speed,the more they consistence with each other.

    4.Cam profile design for two-stroke mode

    The optimization design for four-stroke intake and exhaust cam indicates that STAC method has better superiority,and it provides a effective approach to the two-stroke cam design.

    When crankshaft speed is 4 000 r/min in fourstroke mode,the camshaft speed is 2 000 r/min.So the two-stroke cam is designed using Excite_TD when the speed ratio from crankshaft to camshaft is 1∶1,the design speed is 2 000 r/min crankshaft speed,right now the camshaft speed is 2 000 r/min too when the simulation environment of valve gear is consistent with four-stroke design environment.At the same time,related simulation data show that the maximum torque of two-stroke mode of VCRPE appears in 2 000 r/min crankshaft speed.Therefore,2 000 r/ min crankshaft speed is selected to design speed in two-stroke cam design.The design method also chose STAC.Parameters of the ramp section select the ramp height and velocity at section end,and shape of the ramp is cosine function.

    4.1.Intake cam profile design

    Design parameters of opening and closing section forintake cam oftwo-stroke mode are shown in Table 8.

    Kinematics evaluation results are shown in Figure 15 and Table 9.All of the parameters meet the design requirements and the valve gear operates stably in 2 000 r/min crankshaft speed.Dynamic evaluation results are shown in Table 10,and all of the parameters meet the dynamic design requirements in this speed.

    Table 8.Parameters for intake cam

    Figure 15.Valve lift,velocity and acceleration for two-stroke intake cam

    Table 9.The kinematic parameters

    Table 10.The dynamics parameters

    4.2.Exhaust cam profile design

    Design parameters of opening and closing section for exhaust cam of two-stroke cam are shown in Table 11.

    Kinematics evaluation results are shown in Figure 16 and Table 12.All of the parameters meet the design requirements and the valve gear operates stably in the crankshaft speed of 2 000 r/min.Dynamic evaluation results are shown in Table 13,and all of the parameters meet the dynamic design requirements in this speed.

    Table 11.Parameters for exhaust cam

    Figure 16.Valve lift,velocity and acceleration for two-stroke exhaust cam

    Table 12.The kinematic parameters

    Table 13.The dynamics parameters

    5.Conclusion

    Through analysis and comparison of POL and STAC method for VCRPE four-stroke cam profile design,it shows that the STAC cam design method has the advantages in cam optimization design,and it can match the cam design requirement flexibly.STAC method is used respectively to design the intake and exhaust cam profile for VCRPE four-stroke and twostroke mode,the optimized cam profile through analysis of the kinematics and dynamics simulation results meet all design requirements and valve gear operated stable,which reach the expected purpose.Solid-Works is used to simulate the multi-body dynamics performance for corresponding VCRPE valve gear,and the dynamics simulation results of Excite_TD and SolidWorks are consistent with each other,which verifies the rationality of cam optimization design.

    [1] Ricardo.Engine uses two and four stroke[J].Professional Engineering,2008(9):55.

    [2] Baccile A,Ceccarini D,Cheng L,et al.An innovative control system for a 2/4 stroke switchable engine[J].SAE paper 07011199,2007.

    [3] Di L M.Study on multiple working mode of reciprocating piston engine[D].Haerbin:Northeast Forestry University,2006.

    [4] Di L M,Yang C J,Zhao Y S.Mathematical Model for Closed-form Solution of Parameters of Cylinder Exhaust Port in Variable-stroke Engine[J].Advanced Materials Research,2010,97:2724-2727.

    [5] Feng R H.Optimized improvement design of the engine valve gear[D].Changsha:Hunan University,2009.

    [6] David J W.Optional design of high-speed valve gear systems[C].SAE paper 94250,1994.

    [7] Cardona A,Geradin M.Kinematical and dynamic analysis of mechanical with cams[J].Computer Methods in Applied Mechanics and Engineering 1993,103:115-134.

    [8] Adametal M.Application of numerical simulation for the analysis of the dynamic behavior of valve gear system[J].International Journal of Vehicle Design 1990(1): 281-29l.

    [9] Roland E.Advanced optimization techniques in valve gear design[C].SAE paper 982004,1998.

    [10]Ademar R F.Procedure to modification of a camshaft profile[C].SAE paper 942407,1994.

    [11]Li Z D.Research on the valve timing technology of variable stroke engine[D].Qinhuangdao:Yanshan University,CHINA,2010.

    [12]Shang H J.Design and calculation of engine valve gear[D].Shanghai:Fudan University Press,1988.

    [13]Sun J.Optimization design of polynomial dynamic cam for 475 diesel engine[J].Small Internal Combustion Engine and Motorcycle,2003,32:20-21.

    [14]Zhang G C.Dynamics analysis and cam profile optimization design of the diesel engine valve gear[D].Tianjin: Tianjin University,2009.

    [15]Ye H F,Hao Z Y,Guo L,et al.Design of unsymmetrical cam profiles for a high speed diesel engine with side mounted camshaft[J].Chinese Internal Combustion Engine Engineering,2010,31:25-29.

    [16]Li H Z,Yuan Z C,Le J B,et al.Development of valve cam design[J].Chinese Internal Combustion Engine Engineering,1989,10:32-37.

    可變循環(huán)往復活塞式發(fā)動機凸輪型線優(yōu)化設計*

    邸立明?1,孫萬利2,師 穎3,王文峰1

    1燕山大學車輛與能源學院,河北秦皇島066004;

    2沈陽佳實司法鑒定所,沈陽 110023;

    3交通運輸部公路科學研究院,北京 100088)

    利用Excite_TD軟件構建了本田125 mL單缸四沖程汽油機配氣機構運動學與動力學模型,通過分析多項動力加速度函數(POL)與分段加速度函數(ISAC)凸輪型線設計方法的性能指標,選擇ISAC優(yōu)化設計可變循環(huán)往復活塞式發(fā)動機(VCRPE)二沖程與四沖程工作模式的進、排氣凸輪型線。為了驗證凸輪型線設計的合理性,利用SolidWorks構建并分析了相應配氣機構的多體動力學模型,對比Excite_TD和SolidWorks模型的動力學分析結果表明,兩者分析數據吻合較好,驗證了凸輪型線優(yōu)化設計的準確性。

    往復活塞式發(fā)動機;可變循環(huán);配氣機構;凸輪型線;優(yōu)化設計

    TK423.41

    10.3969/j.issn.1001-3881.2014.06.003

    2013-10-26

    猜你喜歡
    活塞式配氣型線
    高次曲線組合型線渦旋盤性能研究*
    機械制造(2020年8期)2020-09-30 06:32:24
    型線絞合導體ZC-YJLHV22-103×630鋁合金電纜的設計和生產
    電線電纜(2018年3期)2018-06-29 07:41:00
    基于GT-POWER的配氣相位對柴油機性能影響的研究
    差動活塞式燃氣自增壓系統(tǒng)參數設計方法
    無閥配氣鑿巖機利用沖擊‘瞬?!档蜎_擊背壓的探索思考
    某柴油機配氣相位分析及裝配工藝改進
    某柴油機配氣機構異常磨損分析及改進
    萊康明航空活塞式發(fā)動機氣門卡阻故障淺析
    往復活塞式壓縮機節(jié)能降耗技術探究
    水輪機活塞式主軸密封的改進
    久久久国产成人精品二区| 中文资源天堂在线| 亚洲一卡2卡3卡4卡5卡精品中文| 色播亚洲综合网| 日日干狠狠操夜夜爽| 中文字幕人妻丝袜一区二区| 少妇被粗大的猛进出69影院| 18禁观看日本| 亚洲全国av大片| 悠悠久久av| 日韩欧美国产一区二区入口| 国产av在哪里看| 亚洲 国产 在线| 国产精品综合久久久久久久免费| 亚洲电影在线观看av| 女同久久另类99精品国产91| 黑人巨大精品欧美一区二区mp4| 亚洲 欧美 日韩 在线 免费| 国产99白浆流出| 三级国产精品欧美在线观看 | 中文字幕久久专区| 久久久精品大字幕| 久久热在线av| 久久久久久久久免费视频了| 午夜福利欧美成人| 亚洲熟妇熟女久久| 亚洲国产精品sss在线观看| 人人妻,人人澡人人爽秒播| 曰老女人黄片| 日本熟妇午夜| 欧美绝顶高潮抽搐喷水| 午夜日韩欧美国产| 一区二区三区国产精品乱码| 午夜精品一区二区三区免费看| 亚洲av电影不卡..在线观看| 亚洲男人的天堂狠狠| 国产精品永久免费网站| 久久热在线av| 黄色毛片三级朝国网站| 两个人看的免费小视频| 欧美中文综合在线视频| 精品欧美一区二区三区在线| 亚洲一区高清亚洲精品| 国产精品美女特级片免费视频播放器 | 久久天堂一区二区三区四区| 99re在线观看精品视频| 日韩成人在线观看一区二区三区| 午夜福利在线在线| 人人妻人人澡欧美一区二区| 亚洲 国产 在线| 人妻夜夜爽99麻豆av| 女同久久另类99精品国产91| 亚洲精品一卡2卡三卡4卡5卡| 亚洲午夜精品一区,二区,三区| 成人高潮视频无遮挡免费网站| 欧美精品亚洲一区二区| 在线十欧美十亚洲十日本专区| 国产精品国产高清国产av| 婷婷亚洲欧美| 全区人妻精品视频| 在线观看免费视频日本深夜| 欧美精品啪啪一区二区三区| 99国产精品一区二区三区| 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 国产免费男女视频| 曰老女人黄片| 亚洲精品在线美女| 亚洲熟女毛片儿| 动漫黄色视频在线观看| 国产熟女午夜一区二区三区| 亚洲熟妇熟女久久| 亚洲av熟女| 国产三级黄色录像| 国产在线观看jvid| 美女大奶头视频| 久久99热这里只有精品18| 精华霜和精华液先用哪个| 国产在线观看jvid| 午夜免费观看网址| 国产精品一区二区三区四区久久| 岛国视频午夜一区免费看| 久久久久久大精品| 一个人观看的视频www高清免费观看 | 精品第一国产精品| 日韩大尺度精品在线看网址| 免费在线观看成人毛片| 这个男人来自地球电影免费观看| 久久久久久久精品吃奶| 色综合站精品国产| 日本五十路高清| 久久久久亚洲av毛片大全| 日韩免费av在线播放| 欧美av亚洲av综合av国产av| 亚洲乱码一区二区免费版| av欧美777| avwww免费| 美女高潮喷水抽搐中文字幕| 亚洲国产精品sss在线观看| 日韩精品免费视频一区二区三区| 日韩 欧美 亚洲 中文字幕| 国产成人系列免费观看| 欧美日本视频| 天天一区二区日本电影三级| 亚洲成人久久性| 亚洲成人精品中文字幕电影| 色噜噜av男人的天堂激情| 欧美黑人巨大hd| 亚洲黑人精品在线| 午夜精品在线福利| 老司机在亚洲福利影院| 亚洲一区中文字幕在线| 国产精品久久久久久久电影 | 午夜影院日韩av| 日韩 欧美 亚洲 中文字幕| 日韩精品免费视频一区二区三区| 久久久久久久久久黄片| 天堂√8在线中文| 性色av乱码一区二区三区2| 香蕉国产在线看| 日韩中文字幕欧美一区二区| 最新美女视频免费是黄的| 久久久久久久久免费视频了| netflix在线观看网站| 国产免费男女视频| 熟女电影av网| av福利片在线观看| 久久久久久人人人人人| 成人三级黄色视频| 婷婷六月久久综合丁香| 男女下面进入的视频免费午夜| 欧美一区二区国产精品久久精品 | av超薄肉色丝袜交足视频| 香蕉国产在线看| 成年女人毛片免费观看观看9| 色av中文字幕| 午夜精品一区二区三区免费看| 亚洲一区二区三区不卡视频| 一边摸一边做爽爽视频免费| 国产成人系列免费观看| 欧美黑人欧美精品刺激| 久久香蕉国产精品| 久久人妻福利社区极品人妻图片| 国产精品久久视频播放| 亚洲人成电影免费在线| 蜜桃久久精品国产亚洲av| 一本精品99久久精品77| 久久精品亚洲精品国产色婷小说| 国产精品一区二区精品视频观看| 亚洲人成网站在线播放欧美日韩| 日本五十路高清| 一区福利在线观看| 国产精品av久久久久免费| 高清毛片免费观看视频网站| 日本撒尿小便嘘嘘汇集6| www.精华液| xxx96com| 亚洲国产精品合色在线| 亚洲成人久久性| 少妇被粗大的猛进出69影院| 久久久国产精品麻豆| 久久亚洲真实| 欧美色欧美亚洲另类二区| 色尼玛亚洲综合影院| 免费观看精品视频网站| 成人国产综合亚洲| 久久久久久九九精品二区国产 | 久9热在线精品视频| 国产精品亚洲av一区麻豆| 夜夜躁狠狠躁天天躁| 一边摸一边做爽爽视频免费| 久久精品国产清高在天天线| 黄色视频,在线免费观看| 两个人视频免费观看高清| xxxwww97欧美| 黑人操中国人逼视频| 黑人操中国人逼视频| 亚洲人成网站在线播放欧美日韩| 日韩成人在线观看一区二区三区| 久久精品影院6| 母亲3免费完整高清在线观看| 国产麻豆成人av免费视频| 国产v大片淫在线免费观看| 国产成人av激情在线播放| 国产视频一区二区在线看| 午夜福利欧美成人| av有码第一页| 午夜福利视频1000在线观看| 成人亚洲精品av一区二区| 精品高清国产在线一区| 成人永久免费在线观看视频| 国产99白浆流出| 啦啦啦免费观看视频1| 国产精品免费一区二区三区在线| 欧美日韩乱码在线| 99热6这里只有精品| 黄色毛片三级朝国网站| 亚洲精品中文字幕在线视频| 免费在线观看完整版高清| 18禁美女被吸乳视频| 国产成人aa在线观看| 国产精品一区二区免费欧美| 久久久久免费精品人妻一区二区| 国产精品综合久久久久久久免费| 国产私拍福利视频在线观看| а√天堂www在线а√下载| 国产真实乱freesex| 露出奶头的视频| 精品免费久久久久久久清纯| 1024手机看黄色片| 国产精品亚洲av一区麻豆| 丁香六月欧美| 熟女电影av网| 一进一出好大好爽视频| 搡老妇女老女人老熟妇| 亚洲欧美激情综合另类| 亚洲av片天天在线观看| 一边摸一边做爽爽视频免费| 亚洲国产精品sss在线观看| 亚洲av美国av| 久久亚洲真实| 亚洲第一欧美日韩一区二区三区| 久久久久国产一级毛片高清牌| 成人一区二区视频在线观看| 黑人巨大精品欧美一区二区mp4| 巨乳人妻的诱惑在线观看| 欧美性猛交黑人性爽| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 黄色a级毛片大全视频| 精品少妇一区二区三区视频日本电影| 国产成人系列免费观看| 嫩草影院精品99| 欧美乱妇无乱码| 国产av一区在线观看免费| 天天一区二区日本电影三级| 男人舔奶头视频| 亚洲熟妇熟女久久| 一进一出抽搐动态| 精品第一国产精品| 三级男女做爰猛烈吃奶摸视频| 亚洲美女黄片视频| 在线观看美女被高潮喷水网站 | 亚洲七黄色美女视频| 成人国语在线视频| 丁香六月欧美| 一区二区三区高清视频在线| 丰满人妻一区二区三区视频av | 亚洲av片天天在线观看| 一边摸一边做爽爽视频免费| 精品久久久久久,| 麻豆国产97在线/欧美 | 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 国产视频一区二区在线看| 在线观看一区二区三区| 久久久久久久精品吃奶| 精品久久久久久成人av| 久久久久国产精品人妻aⅴ院| 国产真人三级小视频在线观看| 99国产综合亚洲精品| 精品一区二区三区四区五区乱码| 怎么达到女性高潮| 国产精品综合久久久久久久免费| www.999成人在线观看| 最近视频中文字幕2019在线8| 美女 人体艺术 gogo| 国产精品野战在线观看| 老熟妇仑乱视频hdxx| 精品乱码久久久久久99久播| 国产精品久久久久久人妻精品电影| 99久久精品国产亚洲精品| 亚洲成av人片在线播放无| 人妻夜夜爽99麻豆av| 亚洲熟妇熟女久久| 亚洲五月婷婷丁香| 欧美中文综合在线视频| 成人特级黄色片久久久久久久| 中文字幕人妻丝袜一区二区| 99热6这里只有精品| 真人做人爱边吃奶动态| 午夜免费成人在线视频| 天天添夜夜摸| 亚洲欧美激情综合另类| 9191精品国产免费久久| 2021天堂中文幕一二区在线观| videosex国产| 久9热在线精品视频| 亚洲av成人精品一区久久| 老司机深夜福利视频在线观看| www日本黄色视频网| 亚洲狠狠婷婷综合久久图片| 亚洲五月婷婷丁香| 久久久久精品国产欧美久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 国产野战对白在线观看| 黄色视频,在线免费观看| 成年免费大片在线观看| 久久久久免费精品人妻一区二区| 女人高潮潮喷娇喘18禁视频| 丝袜美腿诱惑在线| 18禁裸乳无遮挡免费网站照片| 校园春色视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 女人高潮潮喷娇喘18禁视频| 嫩草影院精品99| x7x7x7水蜜桃| netflix在线观看网站| 久久久久亚洲av毛片大全| 久久久久久国产a免费观看| 在线观看日韩欧美| 国产aⅴ精品一区二区三区波| 久热爱精品视频在线9| 国产欧美日韩一区二区三| 最近最新中文字幕大全免费视频| 色尼玛亚洲综合影院| 又紧又爽又黄一区二区| 国内揄拍国产精品人妻在线| 琪琪午夜伦伦电影理论片6080| 日本撒尿小便嘘嘘汇集6| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久久电影 | 久久精品成人免费网站| 伦理电影免费视频| a级毛片在线看网站| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 婷婷精品国产亚洲av在线| 亚洲av第一区精品v没综合| 久久精品国产99精品国产亚洲性色| 亚洲va日本ⅴa欧美va伊人久久| 国产不卡一卡二| www.www免费av| 级片在线观看| 精品福利观看| 日韩成人在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 日韩欧美 国产精品| 亚洲精品国产精品久久久不卡| 97超级碰碰碰精品色视频在线观看| 欧美性长视频在线观看| 美女大奶头视频| 亚洲一码二码三码区别大吗| 欧美一级毛片孕妇| 亚洲色图av天堂| 一区福利在线观看| 99国产综合亚洲精品| 欧美人与性动交α欧美精品济南到| 可以在线观看毛片的网站| 欧美又色又爽又黄视频| 久久久国产精品麻豆| 国产69精品久久久久777片 | 人妻丰满熟妇av一区二区三区| 亚洲成人久久性| 欧美性猛交╳xxx乱大交人| 精品欧美一区二区三区在线| 手机成人av网站| 国产91精品成人一区二区三区| 日本一二三区视频观看| 精品一区二区三区av网在线观看| 国产爱豆传媒在线观看 | 日韩大尺度精品在线看网址| 999精品在线视频| 国产精品九九99| 国产单亲对白刺激| 精品不卡国产一区二区三区| avwww免费| 免费在线观看完整版高清| 一本一本综合久久| 日日夜夜操网爽| 好男人电影高清在线观看| 18禁国产床啪视频网站| a在线观看视频网站| 精品久久久久久成人av| 国产伦在线观看视频一区| 久久精品国产亚洲av香蕉五月| 黄片大片在线免费观看| 国产亚洲精品av在线| 日韩欧美一区二区三区在线观看| 日韩有码中文字幕| 午夜精品一区二区三区免费看| 国产爱豆传媒在线观看 | 久久久久久久久免费视频了| 亚洲全国av大片| 两个人视频免费观看高清| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 露出奶头的视频| 免费看日本二区| 国产黄片美女视频| 美女黄网站色视频| 三级毛片av免费| 搡老熟女国产l中国老女人| 日韩成人在线观看一区二区三区| 久久精品成人免费网站| 国产激情偷乱视频一区二区| 亚洲乱码一区二区免费版| 1024视频免费在线观看| 国内精品久久久久精免费| 神马国产精品三级电影在线观看 | 亚洲中文日韩欧美视频| 久久人人精品亚洲av| 一级片免费观看大全| 伦理电影免费视频| ponron亚洲| 国产成人精品久久二区二区91| 国产aⅴ精品一区二区三区波| 亚洲 国产 在线| 一级片免费观看大全| 可以在线观看毛片的网站| avwww免费| 毛片女人毛片| 国产免费男女视频| 亚洲精品美女久久av网站| 国产v大片淫在线免费观看| 妹子高潮喷水视频| 国产黄色小视频在线观看| 久久午夜综合久久蜜桃| 午夜视频精品福利| 欧美成人性av电影在线观看| 亚洲av成人精品一区久久| 国产三级在线视频| 日本一二三区视频观看| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 国产视频内射| 成人特级黄色片久久久久久久| 淫秽高清视频在线观看| 亚洲第一电影网av| 精品人妻1区二区| 高清毛片免费观看视频网站| 9191精品国产免费久久| 亚洲成人国产一区在线观看| 很黄的视频免费| 午夜两性在线视频| 午夜福利免费观看在线| 国产激情久久老熟女| 我的老师免费观看完整版| av视频在线观看入口| 啦啦啦观看免费观看视频高清| 美女免费视频网站| 真人一进一出gif抽搐免费| 欧美性长视频在线观看| 日本黄大片高清| avwww免费| 亚洲人成电影免费在线| 免费无遮挡裸体视频| 一个人免费在线观看的高清视频| 国内久久婷婷六月综合欲色啪| av中文乱码字幕在线| 亚洲五月婷婷丁香| 久久中文字幕一级| 成人亚洲精品av一区二区| 大型黄色视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 在线国产一区二区在线| 少妇被粗大的猛进出69影院| 两人在一起打扑克的视频| 国产精品一区二区三区四区免费观看 | 免费在线观看影片大全网站| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3| 人人妻人人澡欧美一区二区| 给我免费播放毛片高清在线观看| 欧美乱码精品一区二区三区| 88av欧美| 可以免费在线观看a视频的电影网站| 熟女电影av网| 精品免费久久久久久久清纯| 男男h啪啪无遮挡| 1024香蕉在线观看| 成人亚洲精品av一区二区| 久久久久久九九精品二区国产 | 欧美性长视频在线观看| 亚洲国产欧洲综合997久久,| 亚洲欧美精品综合久久99| 老汉色∧v一级毛片| 久久精品国产亚洲av高清一级| 丁香六月欧美| 美女黄网站色视频| 国产69精品久久久久777片 | 身体一侧抽搐| 麻豆国产97在线/欧美 | 免费无遮挡裸体视频| 嫩草影院精品99| 一级毛片女人18水好多| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| www.自偷自拍.com| 母亲3免费完整高清在线观看| 国产aⅴ精品一区二区三区波| 妹子高潮喷水视频| 国产av不卡久久| 美女 人体艺术 gogo| 国产熟女xx| 每晚都被弄得嗷嗷叫到高潮| 一进一出好大好爽视频| 免费观看人在逋| 国产69精品久久久久777片 | 精品久久久久久久人妻蜜臀av| 91麻豆精品激情在线观看国产| 国产单亲对白刺激| 18禁观看日本| 禁无遮挡网站| 黄色成人免费大全| 69av精品久久久久久| 日本熟妇午夜| 叶爱在线成人免费视频播放| 精品国内亚洲2022精品成人| 国产欧美日韩一区二区三| 久久久国产成人精品二区| 日韩免费av在线播放| 俺也久久电影网| 亚洲精品一卡2卡三卡4卡5卡| 叶爱在线成人免费视频播放| 国产激情偷乱视频一区二区| 九色成人免费人妻av| 久久热在线av| 国产蜜桃级精品一区二区三区| 国产精品亚洲av一区麻豆| 日韩免费av在线播放| 给我免费播放毛片高清在线观看| 成年版毛片免费区| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 听说在线观看完整版免费高清| 看黄色毛片网站| 黄色视频不卡| 精品国产乱子伦一区二区三区| 黄片小视频在线播放| 正在播放国产对白刺激| 色综合亚洲欧美另类图片| 国产单亲对白刺激| ponron亚洲| 亚洲欧美日韩高清专用| 亚洲精品一卡2卡三卡4卡5卡| 国产高清有码在线观看视频 | 免费在线观看影片大全网站| 久久久久九九精品影院| 在线a可以看的网站| 母亲3免费完整高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 最近最新免费中文字幕在线| 日韩国内少妇激情av| 听说在线观看完整版免费高清| 看片在线看免费视频| 美女高潮喷水抽搐中文字幕| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 狂野欧美激情性xxxx| 国产在线观看jvid| 99久久精品国产亚洲精品| 神马国产精品三级电影在线观看 | 欧美日韩精品网址| av视频在线观看入口| 欧美激情久久久久久爽电影| 最近最新中文字幕大全免费视频| 欧美精品啪啪一区二区三区| 一夜夜www| 国产伦在线观看视频一区| 国产精品精品国产色婷婷| 中文资源天堂在线| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 女人爽到高潮嗷嗷叫在线视频| 香蕉丝袜av| 亚洲最大成人中文| 欧美zozozo另类| 成年版毛片免费区| 精品久久久久久久末码| 黑人巨大精品欧美一区二区mp4| 久久精品成人免费网站| 在线观看免费视频日本深夜| 午夜免费激情av| 欧美三级亚洲精品| 亚洲欧美日韩高清专用| 一进一出抽搐gif免费好疼| 一级作爱视频免费观看| 悠悠久久av| 午夜激情福利司机影院| 人人妻人人澡欧美一区二区| 欧美3d第一页| 一a级毛片在线观看| 婷婷丁香在线五月| 国内揄拍国产精品人妻在线| 日本在线视频免费播放| 国产精品久久久人人做人人爽| www.精华液| 久久精品国产亚洲av香蕉五月| 蜜桃久久精品国产亚洲av| 91麻豆精品激情在线观看国产| 久久国产精品人妻蜜桃| 国产精品一区二区三区四区久久| 亚洲av成人精品一区久久| 亚洲精品美女久久av网站| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| 韩国av一区二区三区四区| 香蕉丝袜av| 日韩精品中文字幕看吧| 欧美在线一区亚洲| 国产亚洲精品久久久久久毛片| 最近最新免费中文字幕在线| 日本黄色视频三级网站网址| 精品欧美国产一区二区三| 久久久久久九九精品二区国产 | 少妇粗大呻吟视频| aaaaa片日本免费| 日韩欧美精品v在线| 欧美国产日韩亚洲一区| 在线免费观看的www视频| 国产精品电影一区二区三区| 国内毛片毛片毛片毛片毛片| 国产精品一区二区精品视频观看|