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

    Multi-body dynamic simulation research on the influence of piston engine crankshaft thrust bearings abrasion①

    2016-12-22 05:46:08YaoZiyun么子云JiangZhinongLiXinJiaYangZhangJinjie
    High Technology Letters 2016年4期

    Yao Ziyun (么子云), Jiang Zhinong, Li Xin, Jia Yang, Zhang Jinjie②

    (*Diagnosis and Self-Recovery Engineering Research Center, Beijing University of Chemical Technology, Beijing 100029, P.R.China)(**Safety Detection and Intelligent System R&D of Beijing Rail Transit Engineering Research Center,Beijing Rail and Transit Design & Research Institute CO., LTD, Beijing 100089, P.R.China)

    ?

    Multi-body dynamic simulation research on the influence of piston engine crankshaft thrust bearings abrasion①

    Yao Ziyun (么子云)*, Jiang Zhinong*, Li Xin**, Jia Yang*, Zhang Jinjie②*

    (*Diagnosis and Self-Recovery Engineering Research Center, Beijing University of Chemical Technology, Beijing 100029, P.R.China)(**Safety Detection and Intelligent System R&D of Beijing Rail Transit Engineering Research Center,Beijing Rail and Transit Design & Research Institute CO., LTD, Beijing 100089, P.R.China)

    Fault diagnosis studying on piston engine, crankshaft and gearbox is focused in this paper. The thrust bearing abrasion caused by axial movement of the crankshaft will affect the force of timing gears and oil pump gears, which will result in the fracture of gears, abnormal ignition, connecting rod cracking and collision of cylinder. Simulation based on CREO software is done to build three-dimensional models of crankshaft and gears of a WP10 diesel engine. The models are imported into ADAMS to complete multi-body dynamics simulations. The force analysis of gears in different kinds of axial movements is finished and variations rules of gear dynamic load is obtained. The presented results show that the crankshaft axial movement can cause overload and vibration on gears. Combined with the realistic case data, the fault feature through simulation research is validated and early warming parameters of gear fault are proposed.

    crankshaft, piston engine, axial movement, gear, multi-body simulation

    0 Introduction

    The piston engines, represented by diesel engine, gasoline engine and gas engine, are widely applied in industries such as ships, vehicles, petrochemical and engineering machinery. Taking ship building industry as an example, driving diesel engines and generator diesel engines take almost 5%~ 8% and 1%~ 3% of total construction costs. The piston engine works based on crankshaft, connecting rods and piston, and the forces acting on the crankshaft of the engine are complex because of the burst pressure of the cylinder.

    Considering the diesel engine as the research focus, to date studies are tended to concentrate on optimal design of the diesel engine to improve performance and operating life. Some researchers have focused on the parts structure, seeking optimization of design of various parts such as crankshaft, camshaft and cooling system[1,2]; other academics have studied the dynamic performance focusing on: vibration, noise and fluidity of operation[3,4]; some pay attention to parts friction and want to improve the friction performance[5,6]; some academics want to improve the fuel economy and reduce diesel engine exhaust pollution and use the computational fluid dynamics simulation technology to optimize the performance of working processes such as combustion, fluid flow and energy transfer[7-10].

    Diesel engine fault monitoring systems can complete data collection, analysis, processing of vibration signals, thermodynamic parameters and lubricating oil characteristics signals. Faults can be detected through signal process and analysis. By extracting fault condition thermal parameters, the system can diagnose the ignition and leakage faults[11,12]. But the approaches based on thermal parameters analysis are limited and have got poor effect for machinery faults. Oil analysis is presented by detecting metal particles content in lubricating oil to indirectly judge the abrasion of metal parts, which is given priority for the detection of ferromagnetic material. Vibration analysis methods use different sensors to collect vibration signals of different working conditions, and process the data with the diesel engine transient speed signals in time and frequency domains[13-15]. However, a diesel engine fault diagnosis system has not been used widely as that of reciprocating compressor[16,17].

    Abrasion of crankshaft thrust bearings results in axial movement, that will cause some serious faults, such as stopping of the engine, serious wear of main bearings, and transmission gears breaking. This research focuses on the influence of piston engine crankshaft axial movement with multi-body dynamic simulation. The forces of gears with the axial movement of crankshaft are simulated by ADAMS software. Some rules of force change with different axial movements of the crankshaft are got which proves that crankshaft axial movement could cause overload and impact on gears. Based on actual data, the fault characteristics are obtained and fault early warning parameters based on vibration signals are proposed.

    1 Force analysis of gears

    Letters in the formulas of this study are defined in the following:

    NomenclatureTTorque(N·mm)FtCircumferentialforce(N)T1Crankshaftgeartorque(N·mm)FaAxialforce(N)nRotatespeed(r/min)Fa1Axialforceincludingtheexternalforceactingonthecrankshaft(N)dPitchdiameter(mm)FrRadialforce(N)zNumberofgearteethFnResultantforce(N)mnNormalmodule(mm)FpContactforce(N)CEquivalentdampingcoefficientFn1Resultantforceincludingthecontactforce(N)kStiffnesscoefficient(N/mm)FdExternalforceactingonthecrankshaft(N)bGeartoothwidth(mm)RCompositecurveradius(mm)ΔvSpeeddifferencesbeforeandafterthecollision(mm/s)R1Equivalentradiusofdrivinggear(mm)eNonlineardampingpowerexponentR2Equivalentradiusofdrivengear(mm)uCollisionforceexponentECompositeelasticmodulus(MPa)mMassofthecrankshaft(kg)E1Elasticmodulusofdrivinggear(MPa) E2Elasticmodulusofdrivengear(MPa) GreeklettersαPressureangle(°)εRestitutioncoefficientβHelicalangle(°)μPoisson’sratioδPenetrationcoefficient(mm)μ1Poisson’sratioofdrivinggearδNormaldeformationvelocity(mm/s)μ2Poisson’sratioofdrivengear

    Fig.1 shows a WP10 type diesel engine as the research object. The crankshaft and gear transmission system of the diesel engine is shown in Fig.2.

    Fig.1 The diesel engine experiment platform

    Fig.2 The internal mechanisms of the crankshaft and gear transmission system

    Action forces between two gears are equal but in opposite direction when they mesh and work in a stable condition normally. T1is defined as the torque of the crankshaft gear. The load of diesel engine crankshaft gear in simulation is about 120000N·mm, which is about 10% of output torque of diesel engine.

    T=T1=120000N·mm

    (1)

    Circumferential force calculation formula:

    (2)

    Axial force calculation formula is

    Fa=Fttanβ=2587.5×tan14°=645.1N

    (3)

    Radial force calculation formula is

    (4)

    In conclusion, the resultant force calculation formula acting on the gear is as follows:

    (5)

    When the crankshaft thrust bearing is worn, if there is an axial force f acting on the crankshaft, axial displacement may appear on the crankshaft thereby causing the crankshaft gear to generate axial displacement along with the crankshaft. There will be a reactive force equal to axial force Fdin the position of the gear contact in order to maintain a balance of the crankshaft.

    The wear of the thrust bearing influencing gears meshing state was simulated by changing the axial displacement of the crankshaft, and there is a positive correlation between the axial displacement and the wear of the thrust bearing. Therefore, in the research a dynamic simulation method is used to analyze the effect of thrust bearing wear on the gears meshing force. In order to avoid the problems such as large calculation scale caused by the direct application of the axial contact constraint to the crankshaft and the errors caused by the setting of contact pair parameters, an equivalent method is applied to simulate the wear of the bearing through the axial displacement constraint to the crankshaft.

    Setting the crankshaft displacement function as f(x), the crankshaft gear axial force calculation formula is

    (6)

    In normal meshing state f(x)=0, and Fa1=Fa. If crankshaft does not move or reach a new steady state with axial movement, the gear stress will remain unchanged. However, the stress will change because of the shorter contact lines of the helical gears in the new steady state. If the axial acceleration always exists, the gear axial force will change with the law affected by crankshaft axial movement. Meanwhile, when crankshaft movement happens, the gear cannot properly connect with the intermediate gear, causing collision and extrusion. The calculation formula is

    (7)

    where k is Hertz stiffness coefficient affected by material properties and interface curvature radius, u is collision force exponent, usually equal to 1.5 for steel, δ is normal penetration coefficient and C is equivalent damping coefficient and δ is the normal deformation velocity. With the theory of Hertz:

    (8)

    E is elastic modulus, reflecting non-linearity degree of the material.

    According to the research in Ref.[18], equivalent damping coefficient C is connected with restitution coefficient ε and stiffness coefficient k. The relationship can be expressed as

    (9)

    Because the input power, output power and rotate speed remain constant, the circumferential force and radial force also remain constant. So the resultant force acting on the gear is

    (10)

    Therefore, when the crankshaft thrust bearing is worn, the crankshaft will drive the axial displacement of the gears, and the connection between gears will be out of the normal state. This will cause the emergence of shock, extrusion and other phenomena, which will result in a sharp increase in the force on the gears. Meanwhile, the two gears are in bad meshing relations, the vibration will be caused by moment and torque, which can be detected by vibration transducers installed on the surface of the body of the engine.

    2 Simulation and result analysis

    2.1 Models of the crankshaft and gears

    Firstly, parts entity models are set up and assembled in CREO software; the models are then transferred into ADAMS software to finish the simulation condition setting. The modeling processes are as follows.

    (1) In CREO, the tooth profile is established according to the gear geometry size and the relationship of different parameters. Translation and rotation commands based on feature copying are used to create the transverse tooth profile in both sides. The helical gear entity models are established by scan, mix and array commands. The size parameters of the gears are shown in Table 1.

    Table 1 Size parameters of the gears

    (2) The crankshaft model is established by rotation and extrusion commands. Size parameters of the crankshaft are shown in Table 2.

    Table 2 Parameters of crankshaft

    (3) In assembly mode, the crankshaft is set to the default location, the crankshaft gear is assembled whose axis coincides with the crankshaft axis, an intermediate gear is assembled to coincide with the end face of the crankshaft gear and set the center distance to the theoretical calculation value. The correct meshing relation of the intermediate gear and the crankshaft gear should be assured.

    (4) There are lots of mode formats imported to ADAMS software. Parasolid file has the best compatibility among them. Therefore, the model assembled in CREO as x_t format is saved, then it is imported to ADAMS.

    (5) Set the component parts material property as “Steel”.

    (6) Set up the connectors, motions and contacts. Set a cylindrical joint between the crankshaft and the ground, a fix joint between the crankshaft gear and the crankshaft. Then set a revolute joint between the intermediate gear and the ground, a contact is set between the two gears. The motions include the crankshaft speed and axial displacement. Meanwhile, setting an opposite torque on the intermediate gear as the load.

    (7) Choose a solver and set simulation time, step size and steps. Then making iterative solution.

    In order to reduce calculation and save time, the models will be simplified. There are no crankshaft oil holes, chamfering and keyway. Since the main analysis object is the tooth stress of gears, the structure of wheel has also been simplified. The established models in software is shown in Fig.3.

    Fig.3 The models in software

    2.2 The contact parameters of calculation

    In ADAMS the main parameters that are required are stiffness coefficient k, collision force exponent n, damping coefficient C and penetration coefficient δ. According to the Hertz contact theory:

    (11)

    (12)

    where R1, R2are the equivalent radius of contact bodies in contact respectively. The reference radius is set as the equivalent radius. R1=46.1mm, R2=76.85mm.

    (13)

    All material of gears is 45# steel, so, μ1=μ2=0.29, E1=E2=207GPa; k=8.09×105N/mm. The other parameters can be set based on experience at u=1.5, C=10, δ=0.05mm.

    2.3 Simulation condition

    Crankshaft speed is 1500rpm, the load is 200Nm, the simulation time is set to 5s and the sampling frequency is 1000Hz. The displacements are 0mm, 0.3mm and 0.5mm. The operating condition of unidirectional movement of the crankshaft is simulated. The crankshaft has moved to a certain position on one side of the gearbox. Crankshaft unidirectional movement is very common in reality, and can lead to bush-burning of the thrust bearing, the failure of gears and wear of main bearings.

    2.4 Reliability of the simulation results

    The stress on the gear surface is analyzed. The simulation model is a rigid body model, and there is a certain deviation between the calculation result and the simulation result. The reference value of the simulation result can be judged through the result that whether the effective value of the force in the non displacement group is close to the theoretical value. The above calculation parameters are the same as the simulation setting parameters and the theoretical calculation and in normal working conditions, the stress of gear is 2837.8N, and the value is close to the effective value of non displacement simulation result of 2722N. Therefore, it can be considered that the simulation result is in the error tolerance range.

    2.5 Analysis of the simulation results of unidirectional motion

    Through the simulation, the contact force on the transmission gears under 0mm, 0.3mm and 0.5mm axial displacements with unidirectional movement conditions are obtained, shown in Fig.4, Fig.5 and Fig.6.

    Fig.4 The contact force on the transmission gears with 0mm

    Fig.5 The contact force on the transmission gears with 0.3mm

    Fig.6 The contact force on the transmission gears with 0.5mm

    Comparing the contact forces on the gear of unidirectional displacement with those of normal working condition, the force will increase as the displacement increases, as shown in Fig.7. The result can be observed in Table 3. Taking the Fourier transform of the force signals of the gear, the frequency spectrums are obtained in Fig.8 to Fig.10. It can be observed that the components of the crankshaft gear rotating frequency and each order frequency multiplication increase as the displacement increases. The main frequency components are shown in Table 4.

    Fig.7 Gear load contrast diagram

    ThecontactforceMaximumvalueAveragevalueEffectivevalueNormal16173N767N2722NUnidirectional0.3mm34893N2795N6921NUnidirectional0.5mm63190N6463N13523N

    Fig.8 Frequency spectrum of contact force data with 0mm

    Fig.9 Frequency spectrum of contact force data with 0.3mm

    The forces of the gear are decomposed into the circumferential, radial and axial directions as shown in Table 5. This highlights that all three reverse forces increase, with the radial force becoming the largest. The forces of the three directions of the gear are not in accordance with the above equations, therefore, the size and the direction of the gear forces are changed. After unidirectional displacement has reached a new steady state, gear stress will increase mainly due to the transmission ratio of the driven wheel and the driving wheel remains unchanged in simulation process. The axial movement of the crankshaft causes the squeeze of thetwo gears, which will lead to a great collision based on the collision force theory in the second section of this paper. This is the primary cause of rapidly increased gear stress. In reality, the driving wheel speed and the driven wheel speed in the movement process cannot maintain the same proportion. The speed fluctuation and vibration will appear until the gear reaches a new stable operation state. After reaching the steady state, the force of the gear will not increase sharply. However, the process of vibration and slippage have an impact between teeth and teeth. The higher the speed, the bigger the impact. This may lead to fracture on the gears.

    Fig.10 Frequency spectrum of contact force data with 0.5mm

    DisplacementGearforceunderthemajorfrequency25Hz50Hz75Hz100Hz450HzNormal759N641N583N570N416NUnidirectional0.3mm2378N1151N865N995N1782NUnidirectional0.5mm4988N1398N1216N2096N4416N

    Table 5 The force on the gear in unidirectional movement

    3 Validation study

    Since it is hard to simulate engine crankshaft thrust bearing wear and fault, in addition it is also difficult to research the influence of crankshaft motion on engine gears through experiment, research is made on the fault diagnosis method of gears after the crankshaft thrust bearing wear of a real case.

    The engine is a V-shaped-16-cylinders diesel engine, the rated speed is 1500rpm and rated output is more than 3000kW. Online monitoring system comprising following transducers: accelerometer, proximity probe and key-phase transducer is installed on the engine. These transducers record the faults in data during different periods. There are sixteen acceleration transducers installed on the cylinder, four on the crankcase and one on the gearbox. It can achieve online monitoring of the vibration of cylinders, crankcase and gearbox. The transducers installation diagram is shown in Fig.11. During running, the crankshaft axial displacement occurs because of the wear of crankshaft axial thrust bearing, which triggers fracture failure of gear in the gearbox. The engine inspection and maintenance results are shown in Fig.12. The online monitoring system captures the signals of the fault.

    Fig.11 The transducers installation diagram of the diesel engine

    Fig.12 The broken gear of the diesel engine

    Fig.13 and Fig.14 are the vibration data of the crankcase and gearbox of the engine monitored by the transducers. It can be observed from the historical trends shown in Fig.13, that at approximately 4:00 am on August 17, 2014, at V4 point of the crankcase vibration and G1 point of gearbox axial vibration are increasing rapidly at the same time. As can be seen from the figure, the maximum vibration RMS values are respectively over 50mm/s and 900mm/s.

    Fig.13 Vibration trend of diesel engine crankcase

    Fig.14 Vibration trend of diesel engine gearbox

    Fig.15 and Fig.16 are the vibration frequency spectrograms of gearbox vibration in common and fault conditions. As can be seen from the spectrograms, in the fault, the measured vibration peak value of 25Hz, 75Hz and 150Hz are significantly increased. These are similar to the simulation results of the above theoretical models.

    Fig.15 Vibration frequency spectrogram of diesel engine gearbox in normal

    Fig.16 Vibration frequency spectrogram of diesel engine gearbox in fault

    4 Conclusion

    This study focuses on the regularity of the force on the gear in a piston engine with crankshaft axial movement. The simulation results prove that the axial displacement of the engine crankshaft will seriously influence the safety and stable operation of the engine; in addition it causes increased stress on the gear. This could lead to shock vibration and fracture of gear even misfire and fracture of connecting rod. The conclusions are as follows:

    (1) From the graph of gear contact force and time, it can be obtained that the unidirectional motion of the crankshaft will cause increased load on the gears. The load is proportional to the displacement.

    As shown in Table 3, the effective value of load under 0.3mm displacement increases to 2.542 times of that in normal, and the effective value of load under 0.5mm displacement increases to 4.968 times of that in normal, which will lead to overload and breakage of the gear teeth.

    (2) From the spectrum graphs of gear contact force, it can be found that with the displacement increasing, different frequencies of the contact force are increasing accordingly.

    (3) Through the analysis of actual engine fault data, theoretical simulation results and the actual fault signal feature has a good consistency. In the fault of the wear on crank thrust bearing, the gears of the engine will be broken and the vibrations of engine crankcase and gearbox will increase significantly. Analyzing the vibration signal frequency components, vibration values of rotating frequency and frequency multiplication increase significantly.

    Therefore, the actual operation of the engine can be monitored by online monitory system through measuring vibrations of crankcase and gearbox. Based on vibration feature extraction and analysis of time domain and frequency domain signals, it can effectively capture the characteristics of the gear faults under crankshaft axial movement condition.

    [ 1] Sun J, Fu Y, Deng M, et al. Study on design optimization of a crankshaft bearing considering crankshaft strength for an internal combustion engine. Proceedings of the Institution of Mechanical Engineers Part J Journal of Engineering Tribology, 2009, 223(6): 919-924

    [ 2] Fontanesi S, Giacopini M. Multiphase CFD-CHT optimization of the cooling jacket and FEM analysis of the engine head of a V6 diesel engine. Applied Thermal Engineering, 2013, 52: 293-303

    [ 3] Sellerbeck P, Nettelbeck C, Heinrichs R, et al. Improving diesel sound quality on engine level and vehicle level-a holistic approach. SAE paper, No.2007-01-2372

    [ 4] Li W, GuF S, Andrew D B, et al. A study of the noise from diesel engines using the independent component analysis. Mechanical Systems and Signal Processing, 2001, 15: 1165-1184

    [ 5] Tamminen J, Sandstrom C E, Andersson P. Influence of load on the tribological conditions in piston ring and cylinder liner contacts in a medium-speed diesel engine. Tribology International, 2006, 39(12): 1643-1652

    [ 6] Jahansson S, Frennfelt C, Killinger A, et al. Frictional evaluation of thermally sprayed coatings applied on the cylinder liner of a heavy duty diesel engine: pilot tribometer analysis and full scale enginetest. Wear, 2011, 273: 82-92

    [ 7] Senecal P K, Reitz R D. Simultaneous reduction of engine emissions and fuel consumption using genetic algorithms and multi-dimensional spray and combustion modeling. SAE paper, No. 2000-01-1890

    [ 8] Wickman D D, Senecal P K, Reitz RD. Diesel engine combustion chamber geometry optimization using genetic algorithms and multidimensional spray and combustion modelling. SAE paper, No. 2001-01-0547

    [ 9] Shrivastav R, Hessel R, Reitz R D. CFD optimization of DI diesel engine performance and emissions using variable intake valve actuation with boost pressure, EGR and multiple injections. SAE paper, No. 2002-01-0959

    [10] Abdelaal M M, Hegab A H. Combustion and emission characteristics of a natural gas-fuelled diesel engine with EGR. Energy Conversion and Management, 2012, 64: 301-312

    [11] Wang X J, Cai Y P, Lin X Z. Diesel engine PT pump fault diagnosis based on the characteristics of its fuel pressure. Ieri Procedia, 2014, 7:84-89

    [12] Tamura M, Saito H, Murata Y, et al. Misfire detection on internal combustion engines using exhaust gas temperature with low sampling rate. Applied Thermal Engineering, 2011, 31:4125-4131

    [13] Li Y, Tse P W, Yang X, et al. EMD-based fault diagnosis for abnormal clearance between contacting components in a diesel engine. Mechanical Systems & Signal Processing, 2010, 24:193-210

    [14] Geveci M, Andrew W O, Franchek M A. An investigation of crankshaft oscillations for cylinder health diagnostics. Mechanical Systems & Signal Processing, 2005, 19:1107-1134

    [15] Liu S, Du R, Yang S. Fault diagnosis for diesel engines by wavelet packet analysis of vibration signal measured on cylinder head. Journal of Vibration Engineering, 2000, 13(4): 577-584 (In Chinese)

    [16] Jin M, Jiang Z N. A new feature reduction method and its application in the reciprocating engine fault diagnosis. Wseas Transactions on Signal Processing, 2012, 8(2): 54-63

    [17] Jiang Z N, Zhang J J, Jin M, et al. An expert system based on multi-source signal integration for reciprocating compressor. WSEAS Transactions on Systems, 2013, 12(5): 266-279

    [18] Lankarani H M, Nikracesh P E. Continuous contact force model for impact analysis in multi body systems. Nonlinear Dynamics, 1994, 5: 193-207

    Yao Ziyun, born in 1978, is studying for his Ph.D degree in Diagnosis and Self-recovery Engineering Research Center of Beijing University of Chemical Technology. His research interests focus on piston engine fault diagnosis and expert system.

    10.3772/j.issn.1006-6748.2016.04.012

    ① Supported by the National Basic Research Program of China (No. 2012CB026000), the Beijing Municipal Science and Technology Commission Program (No. Z131100004113004) and the Fundamental Research Funds for the Central Universities (JD1506).

    ② To whom correspondence should be addressed. E-mail: zjj87427@163.com Received on May 31, 2016*

    久久久久久久久久人人人人人人| 国产爽快片一区二区三区| 久久久久人妻精品一区果冻| 欧美精品一区二区大全| 妹子高潮喷水视频| 久久热在线av| 欧美日韩精品成人综合77777| 国产免费现黄频在线看| 欧美精品高潮呻吟av久久| 激情视频va一区二区三区| 热re99久久国产66热| 99久久中文字幕三级久久日本| 亚洲成人av在线免费| 欧美精品高潮呻吟av久久| 少妇熟女欧美另类| 在线观看美女被高潮喷水网站| 日产精品乱码卡一卡2卡三| 久久久久精品人妻al黑| 国产日韩欧美亚洲二区| 中文乱码字字幕精品一区二区三区| 亚洲国产精品一区三区| 91aial.com中文字幕在线观看| 美女高潮到喷水免费观看| 香蕉精品网在线| av卡一久久| 黑人猛操日本美女一级片| 男人操女人黄网站| 午夜福利在线免费观看网站| 久久精品国产鲁丝片午夜精品| 国产在视频线精品| 少妇人妻 视频| 日本黄色日本黄色录像| 免费看av在线观看网站| 我的亚洲天堂| 亚洲综合色惰| 青草久久国产| 欧美xxⅹ黑人| 99re6热这里在线精品视频| tube8黄色片| 免费黄色在线免费观看| 午夜免费鲁丝| 搡女人真爽免费视频火全软件| 永久网站在线| 宅男免费午夜| 女人高潮潮喷娇喘18禁视频| 热99国产精品久久久久久7| 多毛熟女@视频| 国产精品 欧美亚洲| 欧美日韩成人在线一区二区| 深夜精品福利| 亚洲欧美色中文字幕在线| 中文字幕精品免费在线观看视频| 国产精品三级大全| 菩萨蛮人人尽说江南好唐韦庄| 人体艺术视频欧美日本| 日韩电影二区| 一本久久精品| 日韩精品有码人妻一区| 欧美日韩综合久久久久久| 国产深夜福利视频在线观看| 久久韩国三级中文字幕| 亚洲伊人久久精品综合| 91精品三级在线观看| 韩国高清视频一区二区三区| 日韩免费高清中文字幕av| 久久综合国产亚洲精品| 哪个播放器可以免费观看大片| 少妇 在线观看| 久久久久视频综合| 国产精品不卡视频一区二区| 国产精品嫩草影院av在线观看| 国产亚洲最大av| 2022亚洲国产成人精品| 久久久久精品人妻al黑| 高清在线视频一区二区三区| 亚洲国产色片| 国产亚洲午夜精品一区二区久久| 亚洲第一区二区三区不卡| 青春草视频在线免费观看| 国产精品蜜桃在线观看| 一级黄片播放器| 亚洲国产精品成人久久小说| 久久影院123| 一本—道久久a久久精品蜜桃钙片| 黄片无遮挡物在线观看| 国产男女超爽视频在线观看| 日韩精品免费视频一区二区三区| 久久久精品区二区三区| 超碰成人久久| 黄色毛片三级朝国网站| av国产精品久久久久影院| 伦理电影免费视频| 成人国产麻豆网| 中文欧美无线码| 久久99精品国语久久久| 天天操日日干夜夜撸| 久久国产精品大桥未久av| 五月开心婷婷网| 久久精品aⅴ一区二区三区四区 | 女的被弄到高潮叫床怎么办| 久久午夜福利片| 国产精品免费大片| 在线观看三级黄色| freevideosex欧美| 精品国产露脸久久av麻豆| av有码第一页| 午夜福利,免费看| 少妇人妻精品综合一区二区| 免费日韩欧美在线观看| 人人妻人人添人人爽欧美一区卜| 免费观看av网站的网址| 国产成人一区二区在线| 日韩视频在线欧美| 欧美 日韩 精品 国产| 亚洲国产最新在线播放| 丰满乱子伦码专区| av福利片在线| 欧美xxⅹ黑人| 日韩一区二区视频免费看| 亚洲国产精品国产精品| 美女国产高潮福利片在线看| 久久韩国三级中文字幕| 综合色丁香网| 制服人妻中文乱码| 国产精品人妻久久久影院| 久久精品久久久久久噜噜老黄| 在线天堂最新版资源| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级黄片播放器| 一本大道久久a久久精品| 亚洲精品久久午夜乱码| 国产欧美日韩综合在线一区二区| 国产男女超爽视频在线观看| 欧美老熟妇乱子伦牲交| 国产日韩欧美在线精品| 国产老妇伦熟女老妇高清| 七月丁香在线播放| 亚洲欧美精品自产自拍| 80岁老熟妇乱子伦牲交| 国产无遮挡羞羞视频在线观看| 国产免费视频播放在线视频| av又黄又爽大尺度在线免费看| 永久免费av网站大全| 久久久a久久爽久久v久久| 少妇熟女欧美另类| 午夜福利视频精品| 少妇的丰满在线观看| 午夜福利网站1000一区二区三区| 午夜福利在线观看免费完整高清在| av有码第一页| 日本欧美视频一区| 日日啪夜夜爽| 免费看av在线观看网站| 在线观看www视频免费| 久久久国产欧美日韩av| 少妇猛男粗大的猛烈进出视频| 亚洲综合色网址| 欧美日韩成人在线一区二区| 亚洲精华国产精华液的使用体验| 久久久久精品人妻al黑| 一级片'在线观看视频| 国产精品二区激情视频| 国产精品欧美亚洲77777| 母亲3免费完整高清在线观看 | 成人毛片60女人毛片免费| 一区二区日韩欧美中文字幕| 国产亚洲精品第一综合不卡| 国产精品国产三级国产专区5o| 国产成人a∨麻豆精品| 国产综合精华液| 亚洲欧洲精品一区二区精品久久久 | 国产一区亚洲一区在线观看| 免费观看av网站的网址| 一区福利在线观看| 国产深夜福利视频在线观看| 飞空精品影院首页| 久久人人97超碰香蕉20202| 香蕉精品网在线| 纯流量卡能插随身wifi吗| 中文字幕亚洲精品专区| 亚洲精品在线美女| 国产精品国产三级国产专区5o| 国产精品一二三区在线看| 精品少妇一区二区三区视频日本电影 | 国产深夜福利视频在线观看| 人妻系列 视频| 一本久久精品| 天天躁日日躁夜夜躁夜夜| 婷婷色综合www| 国产精品嫩草影院av在线观看| 高清欧美精品videossex| 99久久精品国产国产毛片| 天天躁日日躁夜夜躁夜夜| 99久久精品国产国产毛片| 欧美 亚洲 国产 日韩一| 女性被躁到高潮视频| 国产精品熟女久久久久浪| 中文字幕精品免费在线观看视频| 毛片一级片免费看久久久久| 久久精品国产综合久久久| 亚洲四区av| 侵犯人妻中文字幕一二三四区| 免费日韩欧美在线观看| 日韩三级伦理在线观看| 看非洲黑人一级黄片| 丝袜人妻中文字幕| 国产白丝娇喘喷水9色精品| 两个人免费观看高清视频| 美女视频免费永久观看网站| 欧美精品亚洲一区二区| 国产精品久久久久久av不卡| 中文字幕av电影在线播放| 丰满迷人的少妇在线观看| 亚洲欧洲精品一区二区精品久久久 | 久久久精品94久久精品| 一本—道久久a久久精品蜜桃钙片| 日韩av在线免费看完整版不卡| 久久99蜜桃精品久久| 一级片免费观看大全| 成人亚洲欧美一区二区av| 人体艺术视频欧美日本| 99国产综合亚洲精品| 欧美激情高清一区二区三区 | 亚洲国产看品久久| 在线观看一区二区三区激情| 国产亚洲av片在线观看秒播厂| 乱人伦中国视频| av线在线观看网站| 在线 av 中文字幕| 欧美日韩av久久| 建设人人有责人人尽责人人享有的| 国产亚洲最大av| 天堂8中文在线网| 午夜老司机福利剧场| 午夜福利影视在线免费观看| 欧美日韩国产mv在线观看视频| 丝袜人妻中文字幕| 午夜日本视频在线| 久久免费观看电影| 国产精品偷伦视频观看了| 国产国语露脸激情在线看| 高清av免费在线| 国产亚洲午夜精品一区二区久久| 91久久精品国产一区二区三区| 毛片一级片免费看久久久久| 国产在视频线精品| 国产片特级美女逼逼视频| 狂野欧美激情性bbbbbb| 999久久久国产精品视频| 欧美av亚洲av综合av国产av | 老司机影院成人| 日韩三级伦理在线观看| 欧美精品人与动牲交sv欧美| 老汉色av国产亚洲站长工具| 国产爽快片一区二区三区| 一区二区三区乱码不卡18| av又黄又爽大尺度在线免费看| 这个男人来自地球电影免费观看 | 国产综合精华液| 桃花免费在线播放| 久久99一区二区三区| 国产成人午夜福利电影在线观看| 一区二区三区激情视频| 一区二区三区乱码不卡18| 亚洲精品国产一区二区精华液| 国产黄频视频在线观看| 国产不卡av网站在线观看| 亚洲精品国产av蜜桃| 国产精品嫩草影院av在线观看| 亚洲欧洲精品一区二区精品久久久 | 日日撸夜夜添| av网站在线播放免费| 国产视频首页在线观看| 日日爽夜夜爽网站| 欧美人与性动交α欧美软件| 色婷婷久久久亚洲欧美| 亚洲欧美中文字幕日韩二区| 97精品久久久久久久久久精品| 男人爽女人下面视频在线观看| 9热在线视频观看99| 欧美成人午夜精品| 在线观看www视频免费| 久久精品国产鲁丝片午夜精品| 日韩制服丝袜自拍偷拍| 97在线人人人人妻| 伊人久久国产一区二区| 久久人人97超碰香蕉20202| 18禁动态无遮挡网站| 精品一品国产午夜福利视频| 精品人妻熟女毛片av久久网站| 久久女婷五月综合色啪小说| 男女边摸边吃奶| 久久99蜜桃精品久久| 亚洲国产av新网站| 可以免费在线观看a视频的电影网站 | 少妇被粗大的猛进出69影院| 亚洲国产日韩一区二区| 国产成人精品久久二区二区91 | 啦啦啦视频在线资源免费观看| 亚洲国产精品成人久久小说| 久久久久久久大尺度免费视频| 伊人久久国产一区二区| av视频免费观看在线观看| 久久久国产精品麻豆| freevideosex欧美| 一区二区av电影网| 女人精品久久久久毛片| 中文字幕最新亚洲高清| 国产av国产精品国产| 欧美日韩精品网址| www.自偷自拍.com| 在线观看免费视频网站a站| 国产免费现黄频在线看| 日韩一本色道免费dvd| 男人操女人黄网站| 最近最新中文字幕免费大全7| 亚洲欧洲精品一区二区精品久久久 | 欧美亚洲日本最大视频资源| 在线观看人妻少妇| 亚洲欧美成人精品一区二区| 欧美xxⅹ黑人| av免费观看日本| 一级毛片 在线播放| 人人澡人人妻人| 日本av手机在线免费观看| 免费在线观看视频国产中文字幕亚洲 | 欧美另类一区| 欧美中文综合在线视频| 人人妻人人添人人爽欧美一区卜| 久久99热这里只频精品6学生| 一级毛片我不卡| 国产精品嫩草影院av在线观看| 亚洲四区av| 亚洲欧美中文字幕日韩二区| 成年人午夜在线观看视频| 久久久久久免费高清国产稀缺| 巨乳人妻的诱惑在线观看| 久久精品国产亚洲av涩爱| 人成视频在线观看免费观看| 又大又黄又爽视频免费| av天堂久久9| av国产久精品久网站免费入址| 国产精品久久久久久精品电影小说| 国产欧美日韩一区二区三区在线| 国产淫语在线视频| 黄片小视频在线播放| 亚洲av电影在线进入| 日韩伦理黄色片| 综合色丁香网| 久久人人爽人人片av| 久久精品aⅴ一区二区三区四区 | 国产一区二区三区av在线| 亚洲综合色网址| 人妻人人澡人人爽人人| 久久久久久久大尺度免费视频| 亚洲精华国产精华液的使用体验| 电影成人av| 亚洲一级一片aⅴ在线观看| 王馨瑶露胸无遮挡在线观看| 在线观看www视频免费| 午夜免费男女啪啪视频观看| 国产精品 国内视频| 亚洲综合色网址| 我要看黄色一级片免费的| 国产乱人偷精品视频| 国产在线视频一区二区| 亚洲精品久久久久久婷婷小说| 七月丁香在线播放| 国语对白做爰xxxⅹ性视频网站| 老汉色∧v一级毛片| 熟女av电影| 精品福利永久在线观看| 各种免费的搞黄视频| 亚洲婷婷狠狠爱综合网| 狠狠精品人妻久久久久久综合| 日韩,欧美,国产一区二区三区| 青草久久国产| 制服人妻中文乱码| 国产黄色视频一区二区在线观看| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人 | 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 久久久久久久久久久久大奶| 日本av免费视频播放| 亚洲国产精品999| 男女国产视频网站| 如何舔出高潮| av有码第一页| 毛片一级片免费看久久久久| 一级毛片我不卡| 亚洲美女搞黄在线观看| 伦精品一区二区三区| 久久人人爽人人片av| 精品一区二区三卡| 你懂的网址亚洲精品在线观看| 七月丁香在线播放| 我要看黄色一级片免费的| 亚洲成人一二三区av| 妹子高潮喷水视频| 亚洲精品久久成人aⅴ小说| 欧美日本中文国产一区发布| 久久精品亚洲av国产电影网| 各种免费的搞黄视频| av电影中文网址| 久久久a久久爽久久v久久| 日本-黄色视频高清免费观看| 国产精品无大码| 亚洲精品乱久久久久久| 国产av国产精品国产| 少妇熟女欧美另类| 黄片小视频在线播放| 91国产中文字幕| 亚洲av电影在线观看一区二区三区| 亚洲五月色婷婷综合| 一区二区日韩欧美中文字幕| 欧美精品人与动牲交sv欧美| 国产亚洲午夜精品一区二区久久| 婷婷成人精品国产| 日韩,欧美,国产一区二区三区| av免费观看日本| av国产精品久久久久影院| 韩国精品一区二区三区| 久久精品国产亚洲av高清一级| 搡老乐熟女国产| 中文字幕亚洲精品专区| 免费观看av网站的网址| www.熟女人妻精品国产| 午夜精品国产一区二区电影| 人体艺术视频欧美日本| 多毛熟女@视频| 999精品在线视频| 成人亚洲精品一区在线观看| av又黄又爽大尺度在线免费看| 在线免费观看不下载黄p国产| 亚洲视频免费观看视频| 日韩精品有码人妻一区| av女优亚洲男人天堂| 熟女av电影| 国产黄频视频在线观看| 99久久综合免费| 国产淫语在线视频| 国产亚洲一区二区精品| 一边亲一边摸免费视频| 久久精品国产自在天天线| 波多野结衣一区麻豆| 免费黄频网站在线观看国产| 91精品三级在线观看| 亚洲国产av新网站| 国产精品嫩草影院av在线观看| 秋霞伦理黄片| 亚洲av欧美aⅴ国产| 97在线视频观看| 啦啦啦啦在线视频资源| 日韩大片免费观看网站| 一边亲一边摸免费视频| 欧美日韩国产mv在线观看视频| 精品少妇一区二区三区视频日本电影 | 久久精品夜色国产| 国产免费一区二区三区四区乱码| 侵犯人妻中文字幕一二三四区| 成人二区视频| 蜜桃在线观看..| 国产精品av久久久久免费| 国产精品一国产av| 乱人伦中国视频| 日本欧美国产在线视频| 日韩电影二区| 午夜av观看不卡| 热99国产精品久久久久久7| 欧美+日韩+精品| 又粗又硬又长又爽又黄的视频| 18禁动态无遮挡网站| 国产成人精品久久二区二区91 | 国产精品久久久久久精品古装| 欧美黄色片欧美黄色片| 免费播放大片免费观看视频在线观看| 欧美日韩精品网址| 久久久久久久久久久免费av| 成人影院久久| 男女午夜视频在线观看| 中文欧美无线码| 在线 av 中文字幕| 久久久久精品人妻al黑| 色哟哟·www| 国产日韩欧美亚洲二区| 久久99精品国语久久久| 久久精品国产综合久久久| 考比视频在线观看| 一级毛片 在线播放| 999精品在线视频| 肉色欧美久久久久久久蜜桃| 校园人妻丝袜中文字幕| 在线观看免费视频网站a站| 日韩精品有码人妻一区| 亚洲三区欧美一区| 国产成人精品无人区| 大片免费播放器 马上看| 1024香蕉在线观看| 各种免费的搞黄视频| 黄色 视频免费看| 老鸭窝网址在线观看| 男男h啪啪无遮挡| 1024香蕉在线观看| 国产精品亚洲av一区麻豆 | 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 国产一区亚洲一区在线观看| 婷婷成人精品国产| 成年美女黄网站色视频大全免费| 日韩成人av中文字幕在线观看| 精品一区二区三卡| 又大又黄又爽视频免费| 少妇熟女欧美另类| 黄频高清免费视频| 国产一区二区 视频在线| 免费黄频网站在线观看国产| 日韩中字成人| 久久久久网色| 叶爱在线成人免费视频播放| 亚洲情色 制服丝袜| 欧美av亚洲av综合av国产av | 亚洲国产精品一区二区三区在线| 观看美女的网站| 日韩视频在线欧美| 日韩制服骚丝袜av| 婷婷色av中文字幕| 在线看a的网站| 国产在线视频一区二区| 国产麻豆69| 久久久久国产精品人妻一区二区| 91精品三级在线观看| 黑人猛操日本美女一级片| 一级爰片在线观看| 亚洲欧美成人综合另类久久久| 性色av一级| 欧美激情高清一区二区三区 | 性少妇av在线| 麻豆乱淫一区二区| 亚洲精品美女久久av网站| 各种免费的搞黄视频| 久久久久久久国产电影| 美女大奶头黄色视频| 免费在线观看视频国产中文字幕亚洲 | 日本av手机在线免费观看| 欧美av亚洲av综合av国产av | 久久99蜜桃精品久久| 中文乱码字字幕精品一区二区三区| 亚洲av电影在线进入| 成人影院久久| 国产野战对白在线观看| 欧美日韩一级在线毛片| 日韩中文字幕欧美一区二区 | 精品少妇久久久久久888优播| 亚洲成人一二三区av| 男人添女人高潮全过程视频| 色网站视频免费| av在线观看视频网站免费| 国产国语露脸激情在线看| 校园人妻丝袜中文字幕| 国产日韩欧美视频二区| 欧美精品高潮呻吟av久久| 免费观看在线日韩| 一区二区三区乱码不卡18| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 国产日韩欧美在线精品| 亚洲精品乱久久久久久| 精品亚洲乱码少妇综合久久| 又大又黄又爽视频免费| 国产精品一二三区在线看| 美女大奶头黄色视频| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费| 久久精品国产综合久久久| 久久精品人人爽人人爽视色| 久久青草综合色| 伦精品一区二区三区| 久久影院123| 久久亚洲国产成人精品v| 中文字幕最新亚洲高清| 日韩一区二区视频免费看| 成人国语在线视频| 极品少妇高潮喷水抽搐| 免费播放大片免费观看视频在线观看| 中文字幕亚洲精品专区| 欧美bdsm另类| 三上悠亚av全集在线观看| av福利片在线| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人看| 日韩制服丝袜自拍偷拍| 久久这里只有精品19| 精品人妻一区二区三区麻豆| 色婷婷av一区二区三区视频| 久久人人97超碰香蕉20202| 日韩伦理黄色片| 天天躁夜夜躁狠狠躁躁| 亚洲欧美成人综合另类久久久| 国产欧美日韩一区二区三区在线| 五月伊人婷婷丁香| 18禁观看日本| 又大又黄又爽视频免费| 丝瓜视频免费看黄片| 99热网站在线观看| 亚洲国产精品国产精品| 亚洲综合色惰| av在线老鸭窝| 一级片'在线观看视频| h视频一区二区三区| 国产片内射在线| 国产精品国产av在线观看| 91久久精品国产一区二区三区| 巨乳人妻的诱惑在线观看| 2018国产大陆天天弄谢|