• <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-09-06 10:53:32LimingDIWanliSUNYingSHIWenfengWANG
    機床與液壓 2014年1期

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

    1School 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

    ?

    Cam profile optimization design of variable cycle reciprocating piston engine*

    Li-ming DI? 1, Wan-li SUN2, Ying SHI3, Wen-feng WANG1

    1School 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

    Abstract: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 four-stroke or two-stroke mode. In order to validate the rationality of the cam profile design, the multi-body 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.

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

    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 four-stroke 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 four-stroke 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 single-cylinder 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 multi-body dynamics simulation model of the corresponding valve gear was established and analyzed by SolidWorks, results from Solidworks multi-body dynamics simulation were used to verify the accuracy of the Excite_TD dynamics simulation.

    2.Intake cam profile design of four-stroke 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:

    cqxq+crxr+csxs)

    (1)

    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

    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 dynamics performance 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 SolidWorks is 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 four-stroke 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

    FullnessRatioMax.Jerk(mm·rad-3)Max.Contactforce/(N·mm-2)LengthofPositiveAccelera-tion/(°)Min.Radius/mmκ0.545189.45662511.74.3

    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 multi-body 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 closing 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 four-stroke 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 for intake cam of two-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

    FullnessRatioMax.Jerk/(mm·rad-3)Max.Contactforce/(N·mm-2)LengthofPositiveAccelera-tion/(°)Min.Radius/mmκ0.5293212.25752510.24.4

    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

    FullnessRatioMax.Jerk/(mm·rad-3)Max.Contactforce/(N·mm-2)LengthofPositiveAccelera-tion/(°)Min.Radius/mmκ0.5326241.56892510.53.8

    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 two-stroke 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. SolidWorks 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.

    References

    [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.

    (Continued from 45 page)

    DOI:10.3969/j.issn.1001-3881.2014.06.003

    Received: 2013-10-26

    *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

    av国产精品久久久久影院| 成人三级做爰电影| 色精品久久人妻99蜜桃| 亚洲精华国产精华精| 热99re8久久精品国产| √禁漫天堂资源中文www| 在线观看免费高清a一片| 建设人人有责人人尽责人人享有的| 精品视频人人做人人爽| 亚洲国产精品一区二区三区在线| 久久久久久亚洲精品国产蜜桃av| 国产成人精品久久二区二区免费| 在线亚洲精品国产二区图片欧美| 99国产极品粉嫩在线观看| 亚洲 国产 在线| 伊人久久大香线蕉亚洲五| 狠狠精品人妻久久久久久综合| 不卡一级毛片| 国产欧美日韩一区二区三| cao死你这个sao货| 可以免费在线观看a视频的电影网站| 亚洲成人手机| 999精品在线视频| 国产精品免费一区二区三区在线 | 叶爱在线成人免费视频播放| 三级毛片av免费| 精品少妇久久久久久888优播| 亚洲欧美激情在线| 美女午夜性视频免费| 97在线人人人人妻| 国产亚洲欧美精品永久| 大香蕉久久成人网| av超薄肉色丝袜交足视频| 交换朋友夫妻互换小说| 午夜福利在线免费观看网站| 国产午夜精品久久久久久| 久久久水蜜桃国产精品网| 精品久久久久久久毛片微露脸| 欧美精品人与动牲交sv欧美| 国产精品免费大片| 曰老女人黄片| 国产又爽黄色视频| 亚洲九九香蕉| aaaaa片日本免费| 少妇猛男粗大的猛烈进出视频| 12—13女人毛片做爰片一| 深夜精品福利| 欧美精品一区二区免费开放| 欧美久久黑人一区二区| 亚洲熟女毛片儿| 国产不卡av网站在线观看| 最近最新免费中文字幕在线| 男女高潮啪啪啪动态图| 精品少妇一区二区三区视频日本电影| 国产免费av片在线观看野外av| 亚洲精华国产精华精| 色尼玛亚洲综合影院| 97人妻天天添夜夜摸| 久久精品熟女亚洲av麻豆精品| 中文字幕人妻丝袜制服| 成年女人毛片免费观看观看9 | 国精品久久久久久国模美| 久久婷婷成人综合色麻豆| 精品少妇黑人巨大在线播放| 欧美成狂野欧美在线观看| 人人妻人人添人人爽欧美一区卜| 在线观看www视频免费| 久久久久久久精品吃奶| 嫁个100分男人电影在线观看| 国产高清videossex| 一级黄色大片毛片| 久久天堂一区二区三区四区| 日韩大码丰满熟妇| 久久人人爽av亚洲精品天堂| 欧美日韩黄片免| www.999成人在线观看| 制服诱惑二区| 精品国内亚洲2022精品成人 | 国产亚洲精品第一综合不卡| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 一本综合久久免费| 在线观看一区二区三区激情| 国产亚洲精品久久久久5区| 香蕉国产在线看| 久久久国产成人免费| 精品国产一区二区久久| 亚洲第一欧美日韩一区二区三区 | 国产精品久久久久久精品古装| tube8黄色片| 亚洲精品国产精品久久久不卡| 欧美日韩福利视频一区二区| 亚洲国产欧美网| 97在线人人人人妻| 国产福利在线免费观看视频| 色尼玛亚洲综合影院| 午夜激情av网站| 久久久久久久国产电影| 亚洲成人免费av在线播放| 欧美变态另类bdsm刘玥| 国产在线视频一区二区| 俄罗斯特黄特色一大片| 在线观看免费日韩欧美大片| 黄片播放在线免费| 18禁观看日本| 精品国产一区二区三区久久久樱花| 视频区图区小说| 亚洲午夜理论影院| 日本欧美视频一区| 久久精品91无色码中文字幕| 久久人人爽av亚洲精品天堂| 亚洲专区中文字幕在线| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 国产日韩欧美视频二区| 日本黄色日本黄色录像| av福利片在线| 搡老熟女国产l中国老女人| 国产在线观看jvid| 久久天躁狠狠躁夜夜2o2o| 久久人妻av系列| 老汉色∧v一级毛片| 午夜福利免费观看在线| 老司机午夜福利在线观看视频 | 国产av精品麻豆| 首页视频小说图片口味搜索| 水蜜桃什么品种好| 99re6热这里在线精品视频| 人人妻人人澡人人看| 美国免费a级毛片| 在线观看免费视频网站a站| 免费看十八禁软件| 精品国产乱子伦一区二区三区| 飞空精品影院首页| 叶爱在线成人免费视频播放| svipshipincom国产片| 丰满饥渴人妻一区二区三| 国产在线免费精品| 亚洲精品乱久久久久久| 一二三四在线观看免费中文在| 99国产精品免费福利视频| 女警被强在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 一级毛片女人18水好多| 午夜免费成人在线视频| 交换朋友夫妻互换小说| 午夜免费鲁丝| 女人高潮潮喷娇喘18禁视频| 久久久久久久精品吃奶| 色视频在线一区二区三区| 国产在线一区二区三区精| 亚洲七黄色美女视频| 国产精品一区二区在线不卡| 午夜精品国产一区二区电影| 少妇裸体淫交视频免费看高清 | 久久久久国内视频| 日日摸夜夜添夜夜添小说| tube8黄色片| 视频在线观看一区二区三区| 亚洲avbb在线观看| 国产精品一区二区免费欧美| 亚洲午夜理论影院| 手机成人av网站| 好男人电影高清在线观看| 欧美变态另类bdsm刘玥| 国产在线免费精品| 久久精品成人免费网站| 日韩中文字幕视频在线看片| 黑人操中国人逼视频| 精品人妻1区二区| cao死你这个sao货| 亚洲精品久久午夜乱码| 色在线成人网| 18在线观看网站| 成人精品一区二区免费| av欧美777| 欧美人与性动交α欧美精品济南到| 亚洲欧美一区二区三区久久| 一进一出好大好爽视频| 久久久精品国产亚洲av高清涩受| 亚洲午夜精品一区,二区,三区| 欧美老熟妇乱子伦牲交| 精品国产一区二区久久| 九色亚洲精品在线播放| av有码第一页| 一级片'在线观看视频| 欧美日韩黄片免| 欧美激情高清一区二区三区| 最新美女视频免费是黄的| 男人舔女人的私密视频| 欧美在线一区亚洲| 久久青草综合色| 国产成人系列免费观看| 久久九九热精品免费| 美国免费a级毛片| 久久中文看片网| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久精品久久久| 巨乳人妻的诱惑在线观看| 十八禁人妻一区二区| 欧美日韩国产mv在线观看视频| 久热这里只有精品99| 一进一出抽搐动态| 亚洲色图av天堂| 国产午夜精品久久久久久| 成人三级做爰电影| 免费在线观看影片大全网站| 免费观看av网站的网址| 久久人妻av系列| 丁香六月天网| 成人手机av| 搡老熟女国产l中国老女人| 久久天堂一区二区三区四区| 桃花免费在线播放| 在线观看www视频免费| 狠狠婷婷综合久久久久久88av| 国产精品一区二区免费欧美| 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| 亚洲一区中文字幕在线| 人人妻人人澡人人爽人人夜夜| 色尼玛亚洲综合影院| 亚洲伊人色综图| 一边摸一边抽搐一进一小说 | 99国产极品粉嫩在线观看| 久久久久国内视频| 亚洲专区字幕在线| 亚洲成av片中文字幕在线观看| 十分钟在线观看高清视频www| 丰满饥渴人妻一区二区三| 国产伦人伦偷精品视频| 日本av手机在线免费观看| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美国产精品va在线观看不卡| 高清视频免费观看一区二区| 亚洲avbb在线观看| 黑人巨大精品欧美一区二区蜜桃| 久久精品成人免费网站| 超碰成人久久| 在线天堂中文资源库| √禁漫天堂资源中文www| 99国产极品粉嫩在线观看| 黄片大片在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看www视频免费| 一夜夜www| 99精品在免费线老司机午夜| 热99久久久久精品小说推荐| 中文字幕色久视频| 91字幕亚洲| a在线观看视频网站| 高清欧美精品videossex| 欧美老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 国产一区二区 视频在线| 久久国产精品影院| av视频免费观看在线观看| 亚洲天堂av无毛| 窝窝影院91人妻| 日本五十路高清| 亚洲人成77777在线视频| 欧美精品一区二区免费开放| 91大片在线观看| 欧美 亚洲 国产 日韩一| 亚洲欧美色中文字幕在线| av福利片在线| 天堂俺去俺来也www色官网| 亚洲精品中文字幕一二三四区 | 久久国产精品大桥未久av| 久久国产精品人妻蜜桃| 久久国产精品男人的天堂亚洲| 成人黄色视频免费在线看| 精品福利观看| 国产有黄有色有爽视频| 成人精品一区二区免费| 亚洲精品美女久久久久99蜜臀| 高清av免费在线| 老熟女久久久| 亚洲天堂av无毛| 久热爱精品视频在线9| 久久国产精品大桥未久av| 免费在线观看日本一区| 热99久久久久精品小说推荐| 精品午夜福利视频在线观看一区 | 国产精品偷伦视频观看了| 免费日韩欧美在线观看| 亚洲精品国产精品久久久不卡| 久久ye,这里只有精品| 亚洲成人免费av在线播放| 别揉我奶头~嗯~啊~动态视频| 国产精品免费视频内射| 国产免费福利视频在线观看| 亚洲国产欧美网| 捣出白浆h1v1| 青青草视频在线视频观看| 考比视频在线观看| 高清av免费在线| 亚洲成人免费av在线播放| 国产精品美女特级片免费视频播放器 | 男女之事视频高清在线观看| 午夜日韩欧美国产| 亚洲精品粉嫩美女一区| 黄色a级毛片大全视频| 日韩欧美一区二区三区在线观看 | 亚洲精品国产色婷婷电影| 日本欧美视频一区| 亚洲精品乱久久久久久| 涩涩av久久男人的天堂| 99久久国产精品久久久| 精品国产超薄肉色丝袜足j| 精品福利永久在线观看| 99国产精品一区二区三区| 老司机影院毛片| 免费在线观看日本一区| 天堂动漫精品| 国产一区二区 视频在线| 久久久精品免费免费高清| 狂野欧美激情性xxxx| 老熟女久久久| 黄色片一级片一级黄色片| 亚洲精品在线美女| 狠狠狠狠99中文字幕| 91老司机精品| 波多野结衣av一区二区av| 最近最新免费中文字幕在线| av不卡在线播放| 91成年电影在线观看| 精品国内亚洲2022精品成人 | 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 国产伦理片在线播放av一区| 国产高清国产精品国产三级| 国产av一区二区精品久久| 大码成人一级视频| 十八禁高潮呻吟视频| 18禁黄网站禁片午夜丰满| 黄色 视频免费看| 欧美日韩成人在线一区二区| 亚洲精品成人av观看孕妇| 欧美在线一区亚洲| 老熟妇乱子伦视频在线观看| 69av精品久久久久久 | 两性夫妻黄色片| 亚洲三区欧美一区| 嫁个100分男人电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品中文字幕一二三四区 | 99riav亚洲国产免费| 欧美老熟妇乱子伦牲交| 亚洲人成伊人成综合网2020| 少妇 在线观看| 韩国精品一区二区三区| 久久亚洲真实| 9191精品国产免费久久| 欧美黑人欧美精品刺激| 成人国语在线视频| 老司机靠b影院| 欧美在线一区亚洲| 精品国产国语对白av| h视频一区二区三区| 亚洲成a人片在线一区二区| 亚洲国产成人一精品久久久| 男女免费视频国产| 日本av手机在线免费观看| 精品亚洲成国产av| 久久久久久久精品吃奶| av一本久久久久| 精品人妻1区二区| 国产福利在线免费观看视频| 丰满人妻熟妇乱又伦精品不卡| 一个人免费看片子| 18禁黄网站禁片午夜丰满| 亚洲精品成人av观看孕妇| 51午夜福利影视在线观看| 久久99热这里只频精品6学生| 免费黄频网站在线观看国产| 桃红色精品国产亚洲av| 国产成人精品久久二区二区免费| 男女无遮挡免费网站观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美国产精品va在线观看不卡| 老汉色av国产亚洲站长工具| 色综合欧美亚洲国产小说| 成人国产一区最新在线观看| 1024视频免费在线观看| 欧美日韩国产mv在线观看视频| 国产不卡av网站在线观看| 日本撒尿小便嘘嘘汇集6| 午夜福利视频精品| 午夜福利乱码中文字幕| 精品熟女少妇八av免费久了| 9热在线视频观看99| av有码第一页| 久久久久精品国产欧美久久久| 午夜免费成人在线视频| 亚洲天堂av无毛| 国产高清视频在线播放一区| 免费在线观看黄色视频的| 999久久久国产精品视频| 99热国产这里只有精品6| 丝瓜视频免费看黄片| 女警被强在线播放| 成人av一区二区三区在线看| 黄色片一级片一级黄色片| 精品人妻1区二区| 亚洲精品乱久久久久久| 久久亚洲精品不卡| 搡老乐熟女国产| 宅男免费午夜| 黄色怎么调成土黄色| 国产三级黄色录像| 亚洲国产欧美网| 欧美激情高清一区二区三区| 精品福利观看| 两性夫妻黄色片| 亚洲国产成人一精品久久久| 日本wwww免费看| 极品少妇高潮喷水抽搐| 亚洲国产欧美在线一区| 久久久精品免费免费高清| 久久午夜综合久久蜜桃| 国产欧美亚洲国产| 欧美乱码精品一区二区三区| 最近最新中文字幕大全免费视频| 一二三四在线观看免费中文在| 日本a在线网址| 一级毛片精品| 精品国产一区二区三区四区第35| 久9热在线精品视频| 精品卡一卡二卡四卡免费| 精品一品国产午夜福利视频| 国产成人av教育| 国产精品亚洲一级av第二区| 国产一区二区三区综合在线观看| 久久久久国产一级毛片高清牌| 在线观看66精品国产| 国产精品av久久久久免费| xxxhd国产人妻xxx| 亚洲,欧美精品.| 国产男女内射视频| 天天躁狠狠躁夜夜躁狠狠躁| 操出白浆在线播放| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 国产淫语在线视频| 亚洲自偷自拍图片 自拍| 日韩三级视频一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 精品国产乱码久久久久久男人| a级毛片黄视频| 超色免费av| 中文字幕制服av| 日本撒尿小便嘘嘘汇集6| 少妇精品久久久久久久| 麻豆av在线久日| 久久国产精品男人的天堂亚洲| 欧美激情 高清一区二区三区| 国产成人影院久久av| 久久香蕉激情| 国产免费视频播放在线视频| 精品亚洲乱码少妇综合久久| 在线看a的网站| 人人妻人人添人人爽欧美一区卜| 国产深夜福利视频在线观看| 亚洲精品中文字幕一二三四区 | 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 亚洲欧美色中文字幕在线| 国产一区二区三区在线臀色熟女 | 飞空精品影院首页| 欧美日韩亚洲综合一区二区三区_| www.精华液| 国产高清视频在线播放一区| 国产精品av久久久久免费| 亚洲avbb在线观看| 另类精品久久| 亚洲九九香蕉| 久久人人爽av亚洲精品天堂| 涩涩av久久男人的天堂| 老鸭窝网址在线观看| 一区二区三区国产精品乱码| 国产又爽黄色视频| tocl精华| 欧美国产精品va在线观看不卡| 国产伦理片在线播放av一区| 久久久久久亚洲精品国产蜜桃av| 激情视频va一区二区三区| 色94色欧美一区二区| 亚洲欧洲日产国产| 成人国产一区最新在线观看| 人成视频在线观看免费观看| 日本av手机在线免费观看| 丁香欧美五月| 波多野结衣一区麻豆| 黄色视频在线播放观看不卡| 欧美激情极品国产一区二区三区| 午夜激情久久久久久久| 国产精品久久久久久精品古装| 在线看a的网站| 丝袜美足系列| 他把我摸到了高潮在线观看 | 热re99久久精品国产66热6| 免费高清在线观看日韩| 免费久久久久久久精品成人欧美视频| 大型黄色视频在线免费观看| 99riav亚洲国产免费| 女人爽到高潮嗷嗷叫在线视频| 欧美老熟妇乱子伦牲交| 色尼玛亚洲综合影院| 日韩精品免费视频一区二区三区| 婷婷丁香在线五月| 人人妻,人人澡人人爽秒播| 一级,二级,三级黄色视频| 大片免费播放器 马上看| 国产不卡av网站在线观看| 日韩大片免费观看网站| 精品一区二区三区视频在线观看免费 | 欧美乱码精品一区二区三区| 国产野战对白在线观看| 欧美激情高清一区二区三区| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 国产午夜精品久久久久久| 夫妻午夜视频| 精品一区二区三卡| 在线av久久热| 午夜视频精品福利| 国产成人影院久久av| 美女国产高潮福利片在线看| 老司机影院毛片| 啦啦啦中文免费视频观看日本| 日韩欧美一区二区三区在线观看 | 精品国产一区二区三区久久久樱花| 亚洲精品久久午夜乱码| 国产免费av片在线观看野外av| 中文字幕色久视频| 久久久久久久国产电影| a在线观看视频网站| 国产精品熟女久久久久浪| 大香蕉久久网| 91老司机精品| 热99久久久久精品小说推荐| 国产单亲对白刺激| 国产在线一区二区三区精| 欧美乱妇无乱码| 侵犯人妻中文字幕一二三四区| 男女免费视频国产| av电影中文网址| 国产亚洲欧美在线一区二区| 成人18禁高潮啪啪吃奶动态图| 久久国产精品男人的天堂亚洲| 日本撒尿小便嘘嘘汇集6| 12—13女人毛片做爰片一| 1024视频免费在线观看| 我要看黄色一级片免费的| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看一区二区三区激情| 免费不卡黄色视频| 久久久久久人人人人人| 一级,二级,三级黄色视频| www日本在线高清视频| 久久中文字幕一级| 国产1区2区3区精品| www.999成人在线观看| 国产日韩一区二区三区精品不卡| 日韩视频在线欧美| 桃花免费在线播放| 亚洲伊人色综图| 操出白浆在线播放| 亚洲精品中文字幕在线视频| 亚洲人成电影观看| 水蜜桃什么品种好| 精品福利永久在线观看| 精品福利观看| 高清黄色对白视频在线免费看| 丁香六月欧美| 十八禁网站免费在线| 天堂8中文在线网| 欧美黑人欧美精品刺激| 伦理电影免费视频| 久久久国产成人免费| 国产视频一区二区在线看| 欧美在线一区亚洲| 久久久久网色| 亚洲国产欧美一区二区综合| 欧美精品啪啪一区二区三区| 黄频高清免费视频| 午夜91福利影院| 国内毛片毛片毛片毛片毛片| 在线 av 中文字幕| 久久精品国产亚洲av高清一级| 亚洲中文日韩欧美视频| 国产精品九九99| 脱女人内裤的视频| 99国产精品免费福利视频| 久久人妻av系列| 男女床上黄色一级片免费看| 搡老乐熟女国产| 国产欧美日韩一区二区三| 久久久精品94久久精品| 精品一品国产午夜福利视频| 婷婷丁香在线五月| 欧美黑人精品巨大| 免费av中文字幕在线| 蜜桃在线观看..| 中文字幕最新亚洲高清| 国产精品久久久人人做人人爽| 午夜久久久在线观看| 国产精品久久久人人做人人爽| 欧美精品av麻豆av| 在线永久观看黄色视频| av天堂久久9| 中文字幕人妻丝袜一区二区| 久久久久精品国产欧美久久久|