• <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*

    欧美 亚洲 国产 日韩一| 日本猛色少妇xxxxx猛交久久| 啦啦啦 在线观看视频| 一区二区日韩欧美中文字幕| 国产成人系列免费观看| 亚洲av在线观看美女高潮| 亚洲成人免费av在线播放| 嫩草影视91久久| www.自偷自拍.com| 91成人精品电影| 精品一品国产午夜福利视频| 日本欧美视频一区| 久久狼人影院| 国产免费又黄又爽又色| 日韩 欧美 亚洲 中文字幕| 国产野战对白在线观看| 日韩视频在线欧美| 考比视频在线观看| 热99久久久久精品小说推荐| 国产av一区二区精品久久| 久久久久精品人妻al黑| 亚洲av日韩在线播放| 亚洲国产毛片av蜜桃av| 美女视频免费永久观看网站| 国产片特级美女逼逼视频| 精品卡一卡二卡四卡免费| 黄色怎么调成土黄色| 多毛熟女@视频| 欧美另类一区| 日本av手机在线免费观看| 国产片内射在线| 亚洲精品国产一区二区精华液| 久久99热这里只频精品6学生| 高潮久久久久久久久久久不卡| 免费一级毛片在线播放高清视频 | 免费不卡黄色视频| 久久九九热精品免费| 99国产精品一区二区蜜桃av | 国产成人系列免费观看| www.999成人在线观看| 国语对白做爰xxxⅹ性视频网站| 婷婷色综合www| 一级a爱视频在线免费观看| 午夜视频精品福利| a级片在线免费高清观看视频| 国产精品一区二区在线观看99| 亚洲国产欧美一区二区综合| 亚洲一区二区三区欧美精品| e午夜精品久久久久久久| 青青草视频在线视频观看| 不卡av一区二区三区| 超碰成人久久| 丁香六月欧美| 国产免费一区二区三区四区乱码| 看免费av毛片| 一本久久精品| 亚洲三区欧美一区| 97精品久久久久久久久久精品| 久久人人97超碰香蕉20202| 精品亚洲乱码少妇综合久久| 久久人人97超碰香蕉20202| 亚洲三区欧美一区| 一边摸一边做爽爽视频免费| 日日摸夜夜添夜夜爱| 亚洲欧洲精品一区二区精品久久久| 交换朋友夫妻互换小说| 欧美在线黄色| 久久人人97超碰香蕉20202| 国产精品免费视频内射| 一本大道久久a久久精品| 国产一区亚洲一区在线观看| 日韩制服骚丝袜av| 麻豆av在线久日| 人成视频在线观看免费观看| 亚洲熟女毛片儿| 精品少妇内射三级| 亚洲情色 制服丝袜| 一本—道久久a久久精品蜜桃钙片| av网站免费在线观看视频| 欧美日韩成人在线一区二区| 久久中文字幕一级| 午夜福利,免费看| 国产又色又爽无遮挡免| 美女中出高潮动态图| 五月天丁香电影| 免费日韩欧美在线观看| 午夜福利影视在线免费观看| 亚洲欧美清纯卡通| 少妇猛男粗大的猛烈进出视频| 亚洲,一卡二卡三卡| 国产成人免费无遮挡视频| 精品高清国产在线一区| 男人舔女人的私密视频| 免费人妻精品一区二区三区视频| 国产欧美日韩精品亚洲av| 99国产精品免费福利视频| 久久久久久久大尺度免费视频| 久热爱精品视频在线9| 一级黄片播放器| 亚洲精品国产区一区二| 亚洲av综合色区一区| 中文字幕制服av| 午夜免费成人在线视频| 亚洲精品在线美女| 999精品在线视频| avwww免费| 五月开心婷婷网| 亚洲av日韩在线播放| 欧美精品人与动牲交sv欧美| 男人操女人黄网站| 2018国产大陆天天弄谢| 日日爽夜夜爽网站| 久久久久网色| 纯流量卡能插随身wifi吗| 91麻豆精品激情在线观看国产 | 久久精品国产亚洲av高清一级| 老鸭窝网址在线观看| 亚洲av电影在线进入| 久久久久精品人妻al黑| 亚洲 国产 在线| 国产欧美日韩一区二区三区在线| 久久人人爽人人片av| 婷婷色av中文字幕| 男人舔女人的私密视频| 热99久久久久精品小说推荐| 久久av网站| 精品免费久久久久久久清纯 | 国产一卡二卡三卡精品| 另类精品久久| 晚上一个人看的免费电影| 人人妻,人人澡人人爽秒播 | 国产主播在线观看一区二区 | xxx大片免费视频| 久久ye,这里只有精品| 大香蕉久久成人网| 久久99一区二区三区| 国产免费福利视频在线观看| 欧美日韩成人在线一区二区| 一区二区av电影网| 日韩一卡2卡3卡4卡2021年| 少妇的丰满在线观看| 啦啦啦在线观看免费高清www| 一级片'在线观看视频| 夫妻午夜视频| 一边亲一边摸免费视频| 熟女少妇亚洲综合色aaa.| 各种免费的搞黄视频| 亚洲精品国产色婷婷电影| 伦理电影免费视频| 亚洲精品久久午夜乱码| av又黄又爽大尺度在线免费看| 成人国产av品久久久| 国产成人精品久久二区二区91| 国产淫语在线视频| 亚洲精品久久成人aⅴ小说| 欧美变态另类bdsm刘玥| 五月开心婷婷网| 在线观看免费视频网站a站| 捣出白浆h1v1| 久久女婷五月综合色啪小说| cao死你这个sao货| 欧美日韩福利视频一区二区| 各种免费的搞黄视频| 免费在线观看黄色视频的| 国产成人系列免费观看| 亚洲成av片中文字幕在线观看| 51午夜福利影视在线观看| 国产欧美日韩综合在线一区二区| 2018国产大陆天天弄谢| 国产成人啪精品午夜网站| 大型av网站在线播放| 久久精品成人免费网站| 国产亚洲欧美在线一区二区| 亚洲一码二码三码区别大吗| 成年av动漫网址| 欧美 亚洲 国产 日韩一| 狂野欧美激情性bbbbbb| 久久久国产欧美日韩av| 久久毛片免费看一区二区三区| 99精国产麻豆久久婷婷| 欧美精品一区二区免费开放| av一本久久久久| 免费观看a级毛片全部| 女警被强在线播放| 久久综合国产亚洲精品| 免费在线观看视频国产中文字幕亚洲 | 黄色a级毛片大全视频| 曰老女人黄片| 日韩制服丝袜自拍偷拍| 午夜免费观看性视频| 看十八女毛片水多多多| 香蕉国产在线看| 一本—道久久a久久精品蜜桃钙片| 国产精品一区二区在线不卡| 国语对白做爰xxxⅹ性视频网站| 999久久久国产精品视频| h视频一区二区三区| 欧美人与善性xxx| 国产一区二区三区综合在线观看| 国产黄色免费在线视频| 久久久久国产一级毛片高清牌| 亚洲精品久久午夜乱码| 亚洲精品日本国产第一区| 97在线人人人人妻| 亚洲人成网站在线观看播放| 交换朋友夫妻互换小说| 99国产综合亚洲精品| 69精品国产乱码久久久| 婷婷色综合大香蕉| 脱女人内裤的视频| 日韩中文字幕视频在线看片| 国产免费视频播放在线视频| 十八禁高潮呻吟视频| 看十八女毛片水多多多| 777米奇影视久久| 一级黄片播放器| 视频区欧美日本亚洲| 欧美成狂野欧美在线观看| 中文字幕av电影在线播放| 丝瓜视频免费看黄片| 大话2 男鬼变身卡| 伦理电影免费视频| 韩国精品一区二区三区| 一级毛片黄色毛片免费观看视频| 18禁裸乳无遮挡动漫免费视频| 好男人视频免费观看在线| 亚洲精品中文字幕在线视频| 亚洲,欧美,日韩| 视频区图区小说| 中文字幕精品免费在线观看视频| 国产日韩欧美在线精品| 纵有疾风起免费观看全集完整版| 18禁裸乳无遮挡动漫免费视频| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区三区在线| 美女脱内裤让男人舔精品视频| 久久这里只有精品19| www日本在线高清视频| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 久久国产精品人妻蜜桃| 婷婷成人精品国产| 人妻一区二区av| 电影成人av| 欧美另类一区| 久久精品亚洲熟妇少妇任你| 黄片小视频在线播放| 欧美日韩成人在线一区二区| 国产成人a∨麻豆精品| 国产亚洲av高清不卡| 深夜精品福利| 日本av免费视频播放| av网站在线播放免费| 欧美日韩黄片免| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲中文av在线| av又黄又爽大尺度在线免费看| 久久国产精品影院| 国产又色又爽无遮挡免| 国产麻豆69| av国产久精品久网站免费入址| 国产主播在线观看一区二区 | 亚洲欧美精品综合一区二区三区| 欧美激情极品国产一区二区三区| 两人在一起打扑克的视频| 9191精品国产免费久久| 成人午夜精彩视频在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲成人手机| 欧美黑人精品巨大| 纵有疾风起免费观看全集完整版| 91精品三级在线观看| 菩萨蛮人人尽说江南好唐韦庄| 久久 成人 亚洲| 午夜免费成人在线视频| 国产福利在线免费观看视频| 大片电影免费在线观看免费| 手机成人av网站| 国产不卡av网站在线观看| 99热国产这里只有精品6| av不卡在线播放| 久久久精品免费免费高清| 日本一区二区免费在线视频| 亚洲欧美清纯卡通| 18禁国产床啪视频网站| 在线av久久热| 国产精品99久久99久久久不卡| 欧美黄色片欧美黄色片| 国产熟女午夜一区二区三区| 我的亚洲天堂| 伊人久久大香线蕉亚洲五| 波野结衣二区三区在线| 精品一区在线观看国产| 亚洲精品一区蜜桃| 男女下面插进去视频免费观看| 欧美xxⅹ黑人| 日韩人妻精品一区2区三区| 啦啦啦在线免费观看视频4| 一个人免费看片子| 欧美日韩视频精品一区| 亚洲精品国产av成人精品| tube8黄色片| 国产亚洲欧美在线一区二区| 久久天堂一区二区三区四区| 国产亚洲精品久久久久5区| 亚洲成国产人片在线观看| 久久国产精品影院| 久久久久精品人妻al黑| 91麻豆av在线| 日韩中文字幕欧美一区二区 | 亚洲伊人久久精品综合| 欧美国产精品va在线观看不卡| 中文字幕人妻丝袜一区二区| 亚洲色图 男人天堂 中文字幕| 美国免费a级毛片| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| av不卡在线播放| 国产精品久久久人人做人人爽| 亚洲熟女精品中文字幕| 丰满少妇做爰视频| 国产91精品成人一区二区三区 | 国产精品久久久人人做人人爽| 精品一区二区三区四区五区乱码 | 国产欧美亚洲国产| 成人国产av品久久久| 天天添夜夜摸| 十分钟在线观看高清视频www| 欧美日韩av久久| tube8黄色片| 天堂中文最新版在线下载| 在线观看国产h片| 亚洲av成人不卡在线观看播放网 | 91麻豆精品激情在线观看国产 | 欧美变态另类bdsm刘玥| av网站免费在线观看视频| 人人澡人人妻人| av在线app专区| 最新的欧美精品一区二区| 伦理电影免费视频| av天堂在线播放| 1024视频免费在线观看| 两个人看的免费小视频| 亚洲专区中文字幕在线| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 777久久人妻少妇嫩草av网站| www.999成人在线观看| 老司机午夜十八禁免费视频| 亚洲国产最新在线播放| av天堂在线播放| 少妇裸体淫交视频免费看高清 | bbb黄色大片| 午夜福利视频精品| 精品一区二区三卡| 99精品久久久久人妻精品| 女人爽到高潮嗷嗷叫在线视频| 国产精品国产av在线观看| 国产国语露脸激情在线看| 免费日韩欧美在线观看| 91精品三级在线观看| 91麻豆av在线| 国产欧美日韩一区二区三区在线| 人成视频在线观看免费观看| av国产精品久久久久影院| 亚洲第一av免费看| 中文字幕人妻熟女乱码| 女性被躁到高潮视频| 欧美日本中文国产一区发布| 啦啦啦在线观看免费高清www| 中文字幕亚洲精品专区| 满18在线观看网站| 欧美久久黑人一区二区| 亚洲av成人不卡在线观看播放网 | 成人亚洲欧美一区二区av| 国产亚洲欧美精品永久| 三上悠亚av全集在线观看| 久久精品成人免费网站| 国产1区2区3区精品| 国产主播在线观看一区二区 | 一二三四在线观看免费中文在| 久久久欧美国产精品| 少妇人妻久久综合中文| 一级a爱视频在线免费观看| 亚洲中文av在线| 免费av中文字幕在线| 精品人妻1区二区| 中文欧美无线码| 老司机亚洲免费影院| 捣出白浆h1v1| 欧美中文综合在线视频| 18禁观看日本| 色精品久久人妻99蜜桃| 人人澡人人妻人| 亚洲国产看品久久| 大香蕉久久成人网| a级毛片黄视频| 看免费成人av毛片| 欧美中文综合在线视频| 老熟女久久久| 日韩人妻精品一区2区三区| 美国免费a级毛片| 在线天堂中文资源库| www.熟女人妻精品国产| 91字幕亚洲| 丝袜人妻中文字幕| 又黄又粗又硬又大视频| 岛国毛片在线播放| 黑人猛操日本美女一级片| 国产亚洲av片在线观看秒播厂| 只有这里有精品99| 久久精品亚洲熟妇少妇任你| 亚洲第一青青草原| av网站在线播放免费| 嫩草影视91久久| 一区二区三区精品91| 欧美性长视频在线观看| 视频区欧美日本亚洲| 丰满迷人的少妇在线观看| tube8黄色片| 在线观看免费午夜福利视频| 亚洲欧美一区二区三区国产| 欧美激情极品国产一区二区三区| 男女国产视频网站| 亚洲欧美色中文字幕在线| 日本午夜av视频| 九色亚洲精品在线播放| 久久中文字幕一级| 中文字幕av电影在线播放| a级片在线免费高清观看视频| 少妇的丰满在线观看| 飞空精品影院首页| 婷婷色麻豆天堂久久| 男人爽女人下面视频在线观看| 久久精品国产亚洲av涩爱| av一本久久久久| 男人添女人高潮全过程视频| 后天国语完整版免费观看| 国产一区二区激情短视频 | 18在线观看网站| 欧美性长视频在线观看| 蜜桃在线观看..| 亚洲精品国产一区二区精华液| av在线app专区| 日韩伦理黄色片| 另类精品久久| 亚洲精品日韩在线中文字幕| 女性被躁到高潮视频| 成人亚洲欧美一区二区av| av在线播放精品| 亚洲国产精品999| 久久久国产一区二区| 涩涩av久久男人的天堂| 亚洲人成电影观看| 亚洲,一卡二卡三卡| 国产又色又爽无遮挡免| www.熟女人妻精品国产| 欧美黑人精品巨大| 欧美中文综合在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利免费观看在线| 亚洲男人天堂网一区| 多毛熟女@视频| 男人添女人高潮全过程视频| 国产成人欧美在线观看 | 久久精品熟女亚洲av麻豆精品| 精品国产一区二区久久| 久久久久久免费高清国产稀缺| 免费看十八禁软件| 日本欧美国产在线视频| 叶爱在线成人免费视频播放| 色播在线永久视频| 看免费成人av毛片| 黄色怎么调成土黄色| 久久亚洲国产成人精品v| 日韩av免费高清视频| 国产真人三级小视频在线观看| 成人国语在线视频| 在线观看人妻少妇| 国产成人av教育| 亚洲精品国产av蜜桃| 国产精品久久久久久人妻精品电影 | 美女国产高潮福利片在线看| 精品国产乱码久久久久久男人| netflix在线观看网站| 国产在线免费精品| av在线播放精品| 成年人免费黄色播放视频| 青青草视频在线视频观看| 尾随美女入室| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 国产一区亚洲一区在线观看| 中文字幕制服av| av福利片在线| 免费在线观看影片大全网站 | 日日摸夜夜添夜夜爱| 欧美激情极品国产一区二区三区| 日本av手机在线免费观看| 免费看不卡的av| 男女国产视频网站| 午夜免费男女啪啪视频观看| 国产激情久久老熟女| 日本欧美视频一区| 亚洲欧洲国产日韩| www.自偷自拍.com| 国产欧美日韩综合在线一区二区| 国语对白做爰xxxⅹ性视频网站| 蜜桃在线观看..| 国产成人免费观看mmmm| 亚洲三区欧美一区| 高清视频免费观看一区二区| 久久久久久亚洲精品国产蜜桃av| 黄网站色视频无遮挡免费观看| 国产免费现黄频在线看| 黑人猛操日本美女一级片| 亚洲专区中文字幕在线| 男女边摸边吃奶| 精品一区二区三区av网在线观看 | 精品一区二区三卡| 女性生殖器流出的白浆| 午夜久久久在线观看| 成人国产一区最新在线观看 | 久久精品国产综合久久久| 国产又色又爽无遮挡免| 免费不卡黄色视频| 国产无遮挡羞羞视频在线观看| 老司机影院毛片| av天堂久久9| 丝袜喷水一区| 午夜av观看不卡| 男人爽女人下面视频在线观看| 日韩,欧美,国产一区二区三区| 大片电影免费在线观看免费| 亚洲国产成人一精品久久久| 永久免费av网站大全| 91精品伊人久久大香线蕉| 赤兔流量卡办理| e午夜精品久久久久久久| 久久人人爽人人片av| 香蕉丝袜av| 国产av国产精品国产| 美女中出高潮动态图| 亚洲国产av新网站| 欧美亚洲日本最大视频资源| 后天国语完整版免费观看| 狠狠精品人妻久久久久久综合| 精品国产一区二区久久| 涩涩av久久男人的天堂| 在线av久久热| 精品人妻一区二区三区麻豆| av片东京热男人的天堂| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9 | 亚洲熟女毛片儿| 成人手机av| 18禁观看日本| 七月丁香在线播放| 新久久久久国产一级毛片| 精品欧美一区二区三区在线| 男女无遮挡免费网站观看| av一本久久久久| www.999成人在线观看| 我的亚洲天堂| 亚洲av日韩精品久久久久久密 | 涩涩av久久男人的天堂| 大码成人一级视频| 美女中出高潮动态图| 国产片内射在线| 91精品三级在线观看| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 亚洲国产精品成人久久小说| 又大又爽又粗| 黄色 视频免费看| 日本一区二区免费在线视频| 新久久久久国产一级毛片| 精品久久蜜臀av无| 涩涩av久久男人的天堂| 女人久久www免费人成看片| 少妇猛男粗大的猛烈进出视频| 国产精品av久久久久免费| 一区二区三区乱码不卡18| 永久免费av网站大全| 麻豆乱淫一区二区| 一区在线观看完整版| 婷婷成人精品国产| 日韩制服骚丝袜av| 欧美97在线视频| 青春草亚洲视频在线观看| 亚洲成av片中文字幕在线观看| 免费一级毛片在线播放高清视频 | 日韩av免费高清视频| 熟女av电影| 少妇猛男粗大的猛烈进出视频| 视频区欧美日本亚洲| 不卡av一区二区三区| 亚洲精品自拍成人| 每晚都被弄得嗷嗷叫到高潮| 国产极品粉嫩免费观看在线| 丝瓜视频免费看黄片| 啦啦啦在线观看免费高清www| 亚洲成人国产一区在线观看 | 女性生殖器流出的白浆| 99九九在线精品视频| 男的添女的下面高潮视频| 国产无遮挡羞羞视频在线观看| 老司机午夜十八禁免费视频| 国产淫语在线视频| 午夜免费鲁丝| 国产av精品麻豆| 亚洲欧美一区二区三区国产| 成人手机av| 涩涩av久久男人的天堂|