• <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在线观看| 欧美日韩视频精品一区| 久久鲁丝午夜福利片| 亚洲精品一二三| 久久久欧美国产精品| 免费少妇av软件| 亚洲欧美日韩另类电影网站| 在现免费观看毛片| 欧美 亚洲 国产 日韩一| 久久精品国产综合久久久| 桃花免费在线播放| 久久精品国产亚洲av高清一级| 国产伦理片在线播放av一区| 飞空精品影院首页| 不卡av一区二区三区| 亚洲精品国产av成人精品| 观看美女的网站| 精品少妇内射三级| 2021少妇久久久久久久久久久| 欧美 亚洲 国产 日韩一| 国产精品女同一区二区软件| 午夜激情久久久久久久| 亚洲国产欧美在线一区| 欧美日韩一级在线毛片| 操出白浆在线播放| bbb黄色大片| 久久性视频一级片| 这个男人来自地球电影免费观看 | 日本欧美国产在线视频| 97精品久久久久久久久久精品| 亚洲 欧美一区二区三区| 色视频在线一区二区三区| 一区二区三区四区激情视频| 丰满少妇做爰视频| 久久人人爽人人片av| 九九爱精品视频在线观看| 亚洲图色成人| 国产又爽黄色视频| 亚洲第一青青草原| 亚洲,欧美精品.| 久久av网站| 国产有黄有色有爽视频| 叶爱在线成人免费视频播放| 美女高潮到喷水免费观看| 日本黄色日本黄色录像| 久久亚洲国产成人精品v| 秋霞伦理黄片| 精品一区二区三卡| 天天影视国产精品| 久久狼人影院| 男人添女人高潮全过程视频| 亚洲视频免费观看视频| 高清视频免费观看一区二区| 久久久久精品人妻al黑| 日韩熟女老妇一区二区性免费视频| 亚洲七黄色美女视频| 久久99一区二区三区| 日韩av在线免费看完整版不卡| 亚洲欧美成人综合另类久久久| 女性生殖器流出的白浆| 在线精品无人区一区二区三| 欧美黑人欧美精品刺激| 精品视频人人做人人爽| 美女午夜性视频免费| 午夜福利乱码中文字幕| 宅男免费午夜| 国产一区二区激情短视频 | 国产av国产精品国产| 精品少妇久久久久久888优播| 国产成人免费观看mmmm| 免费不卡黄色视频| 久久久国产一区二区| av天堂久久9| 啦啦啦中文免费视频观看日本| 观看美女的网站| 热re99久久国产66热| 丰满少妇做爰视频| 成年人午夜在线观看视频| 秋霞在线观看毛片| 国产亚洲av高清不卡| 久久久精品区二区三区| 亚洲美女搞黄在线观看| 亚洲欧美激情在线| 国产 精品1| 纯流量卡能插随身wifi吗| tube8黄色片| 男人添女人高潮全过程视频| 午夜福利免费观看在线| videosex国产| 欧美黑人精品巨大| 日韩中文字幕欧美一区二区 | 大香蕉久久网| 亚洲欧美激情在线| 日韩中文字幕欧美一区二区 | 曰老女人黄片| 国产精品免费视频内射| 一区二区三区激情视频| 久久久久国产精品人妻一区二区| 丰满少妇做爰视频| 午夜福利视频精品| 丝袜脚勾引网站| 国产人伦9x9x在线观看| 亚洲成人国产一区在线观看 | 亚洲精品久久午夜乱码| 国产免费一区二区三区四区乱码| 麻豆乱淫一区二区| 日韩一区二区三区影片| 免费黄色在线免费观看| 亚洲成人一二三区av| 在线观看免费视频网站a站| 九草在线视频观看| 青草久久国产| 欧美日韩成人在线一区二区| 日本av免费视频播放| 欧美久久黑人一区二区| 国产精品久久久久久精品电影小说| 国产亚洲一区二区精品| 无遮挡黄片免费观看| 国产亚洲一区二区精品| 高清欧美精品videossex| 国产 精品1| 在线观看人妻少妇| 男女午夜视频在线观看| 精品人妻熟女毛片av久久网站| 九草在线视频观看| 欧美在线一区亚洲| 免费高清在线观看视频在线观看| 精品国产一区二区久久| 精品少妇黑人巨大在线播放| 男人添女人高潮全过程视频| 国产爽快片一区二区三区| 久久久久精品久久久久真实原创| 国产激情久久老熟女| 成人18禁高潮啪啪吃奶动态图| 视频区图区小说| 亚洲免费av在线视频| 综合色丁香网| 国产精品免费大片| 麻豆av在线久日| 日韩av在线免费看完整版不卡| 国产精品无大码| 精品午夜福利在线看| 看免费av毛片| 国产精品免费视频内射| 久久人人97超碰香蕉20202| 老司机在亚洲福利影院| 欧美日韩成人在线一区二区| 久久精品人人爽人人爽视色| 九草在线视频观看| 精品国产露脸久久av麻豆| 免费观看av网站的网址| 日韩大码丰满熟妇| 国产麻豆69| 久久久久久久久免费视频了| 最近手机中文字幕大全| 超碰97精品在线观看| av在线app专区| 午夜免费男女啪啪视频观看| 波野结衣二区三区在线| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av国产av综合av卡| 亚洲av成人精品一二三区| 亚洲国产av新网站| 国产精品免费大片| 精品少妇内射三级| 少妇 在线观看| 国产深夜福利视频在线观看| 亚洲成人国产一区在线观看 | 丁香六月天网| 秋霞伦理黄片| 国产av码专区亚洲av| 国产精品无大码| av网站在线播放免费| 日日摸夜夜添夜夜爱| 下体分泌物呈黄色| 久久久久久人妻| 人人妻人人爽人人添夜夜欢视频| 99国产综合亚洲精品| 久久久亚洲精品成人影院| 99香蕉大伊视频| 免费久久久久久久精品成人欧美视频| 丰满少妇做爰视频| 国产亚洲一区二区精品| 九九爱精品视频在线观看| 精品亚洲乱码少妇综合久久| 午夜影院在线不卡| 不卡视频在线观看欧美| 十八禁人妻一区二区| 欧美国产精品一级二级三级| 不卡av一区二区三区| 高清在线视频一区二区三区| 精品少妇黑人巨大在线播放| 一级片免费观看大全| 精品国产一区二区久久| 乱人伦中国视频| 亚洲av美国av| 电影成人av| 人妻久久中文字幕网| 精品国产乱子伦一区二区三区| 欧美激情久久久久久爽电影 | 久久久久久国产a免费观看| 日本a在线网址| 国内精品久久久久久久电影| 啦啦啦观看免费观看视频高清 | 母亲3免费完整高清在线观看| 免费不卡黄色视频| 免费久久久久久久精品成人欧美视频| 一本综合久久免费| 大型黄色视频在线免费观看| 国产欧美日韩精品亚洲av| 十八禁人妻一区二区| 人妻久久中文字幕网| 国产伦一二天堂av在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久精品国产清高在天天线| 色播亚洲综合网| 国产区一区二久久| 黑人操中国人逼视频| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 老鸭窝网址在线观看| 99久久综合精品五月天人人| 欧美 亚洲 国产 日韩一| 欧美日本亚洲视频在线播放| 一个人免费在线观看的高清视频| 国产熟女午夜一区二区三区| 91老司机精品| 变态另类成人亚洲欧美熟女 | 国产亚洲欧美精品永久| 国产日韩一区二区三区精品不卡| 9热在线视频观看99| 精品人妻1区二区| 久久香蕉国产精品| 成年人黄色毛片网站| 日本在线视频免费播放| 中文字幕av电影在线播放| 国产99白浆流出| 黄色女人牲交| 亚洲黑人精品在线| 欧美丝袜亚洲另类 | 在线观看日韩欧美| 久久草成人影院| 咕卡用的链子| 999精品在线视频| 国产黄a三级三级三级人| 国产熟女午夜一区二区三区| 国产亚洲欧美精品永久| 亚洲va日本ⅴa欧美va伊人久久| 村上凉子中文字幕在线| 国产伦人伦偷精品视频| 操美女的视频在线观看| 一级毛片高清免费大全| 国产精品久久久av美女十八| 国产av精品麻豆| 欧美黄色淫秽网站| av电影中文网址| 涩涩av久久男人的天堂| 大陆偷拍与自拍| 99在线视频只有这里精品首页| 黄片播放在线免费| 久久精品91蜜桃| 国产欧美日韩精品亚洲av| 最近最新免费中文字幕在线| 热99re8久久精品国产| 成人欧美大片| 国产精品1区2区在线观看.| 久久国产精品男人的天堂亚洲| 欧美激情高清一区二区三区| 免费在线观看亚洲国产| 国产av一区在线观看免费| 精品不卡国产一区二区三区| 久久精品aⅴ一区二区三区四区| 亚洲一码二码三码区别大吗| 最新美女视频免费是黄的| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 一级,二级,三级黄色视频| 久久伊人香网站| 色哟哟哟哟哟哟| 女性生殖器流出的白浆| 日本撒尿小便嘘嘘汇集6| 国产在线精品亚洲第一网站| 9热在线视频观看99| 精品国产美女av久久久久小说| 一级,二级,三级黄色视频| 黄色丝袜av网址大全| 变态另类成人亚洲欧美熟女 | 熟女少妇亚洲综合色aaa.| 免费在线观看视频国产中文字幕亚洲| 搞女人的毛片| 亚洲专区中文字幕在线| 99热只有精品国产| 久久九九热精品免费| 美女 人体艺术 gogo| 91老司机精品| 亚洲熟女毛片儿| 亚洲少妇的诱惑av| 欧美日韩亚洲综合一区二区三区_| 欧美日韩福利视频一区二区| 9色porny在线观看| 日韩中文字幕欧美一区二区| 十八禁网站免费在线| 国产av精品麻豆| 国产av又大| 久久热在线av| 日韩免费av在线播放| 黄色视频不卡| 国产精品秋霞免费鲁丝片| 国产免费男女视频| 亚洲自拍偷在线| av超薄肉色丝袜交足视频| 久久久久久久久免费视频了| 亚洲一区中文字幕在线| 国产精品免费视频内射| 91精品三级在线观看| 日韩有码中文字幕| 国产精品野战在线观看| 亚洲精品在线观看二区| 51午夜福利影视在线观看| 色尼玛亚洲综合影院| 久久香蕉精品热| 99精品欧美一区二区三区四区| 制服诱惑二区| 黄色毛片三级朝国网站| av天堂在线播放| 人成视频在线观看免费观看| 欧美一级毛片孕妇| 久久精品91蜜桃| 美女高潮喷水抽搐中文字幕| 后天国语完整版免费观看| 真人做人爱边吃奶动态| 久久国产乱子伦精品免费另类| 国产精品av久久久久免费| 亚洲男人的天堂狠狠| 变态另类成人亚洲欧美熟女 | svipshipincom国产片| 欧美精品啪啪一区二区三区| 丰满的人妻完整版| 在线观看日韩欧美| 国产精品亚洲美女久久久| 国产不卡一卡二| 午夜精品在线福利| 亚洲av成人一区二区三| 亚洲美女黄片视频| 久久精品aⅴ一区二区三区四区| 国产亚洲av高清不卡| 一个人观看的视频www高清免费观看 | 国产精品一区二区在线不卡| 一进一出抽搐动态| 琪琪午夜伦伦电影理论片6080| 可以在线观看毛片的网站| 变态另类成人亚洲欧美熟女 | 可以免费在线观看a视频的电影网站| 少妇熟女aⅴ在线视频| 欧美黄色片欧美黄色片| 精品无人区乱码1区二区| 一二三四在线观看免费中文在| 欧美日本中文国产一区发布| 免费久久久久久久精品成人欧美视频| 一级黄色大片毛片| 久久天堂一区二区三区四区| 美女 人体艺术 gogo| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 国产欧美日韩综合在线一区二区| 成人精品一区二区免费| 一区二区日韩欧美中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一区二区三区不卡视频| 国产精品久久久人人做人人爽| 麻豆国产av国片精品| 美国免费a级毛片| 国产欧美日韩一区二区三| 午夜精品久久久久久毛片777| 亚洲最大成人中文| 十八禁网站免费在线| 成年女人毛片免费观看观看9| 亚洲中文字幕一区二区三区有码在线看 | 免费少妇av软件| 69精品国产乱码久久久| 欧美日韩瑟瑟在线播放| 午夜福利,免费看| 最新在线观看一区二区三区| 一二三四在线观看免费中文在| 久久久久久大精品| 伦理电影免费视频| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 女人爽到高潮嗷嗷叫在线视频| 欧美午夜高清在线| 麻豆国产av国片精品| 欧美日韩瑟瑟在线播放| 免费高清视频大片| 夜夜夜夜夜久久久久| 手机成人av网站| 国产一区二区在线av高清观看| 人妻久久中文字幕网| 动漫黄色视频在线观看| 久久久精品欧美日韩精品| 国产精品一区二区免费欧美| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av高清一级| 久久国产精品人妻蜜桃| 免费少妇av软件| 99在线视频只有这里精品首页| 精品福利观看| 国产高清激情床上av| 中亚洲国语对白在线视频| 可以在线观看毛片的网站| 国产欧美日韩一区二区三区在线| 国产一卡二卡三卡精品| 三级毛片av免费| 亚洲中文字幕日韩| 女性生殖器流出的白浆| 精品久久蜜臀av无| 侵犯人妻中文字幕一二三四区| 大陆偷拍与自拍| 免费在线观看日本一区| 在线观看免费视频网站a站| 久久香蕉激情| 亚洲一区二区三区不卡视频| АⅤ资源中文在线天堂| 法律面前人人平等表现在哪些方面| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 久久久久久亚洲精品国产蜜桃av| 999精品在线视频| 午夜福利在线观看吧| 欧美乱色亚洲激情| 啦啦啦 在线观看视频| 高清在线国产一区| 乱人伦中国视频| 成人手机av| 亚洲午夜理论影院| av在线播放免费不卡| 国产精品爽爽va在线观看网站 | 别揉我奶头~嗯~啊~动态视频| 色婷婷久久久亚洲欧美| 免费在线观看完整版高清| 成人三级做爰电影| 亚洲成人久久性| 亚洲国产欧美一区二区综合| 无限看片的www在线观看| 天堂动漫精品| 国产麻豆69| 国产成人欧美| 久久精品国产亚洲av高清一级| 黄色a级毛片大全视频| 久久天堂一区二区三区四区| e午夜精品久久久久久久| 久久婷婷成人综合色麻豆| 亚洲国产高清在线一区二区三 | 免费一级毛片在线播放高清视频 | 国产97色在线日韩免费| 亚洲中文字幕一区二区三区有码在线看 | 午夜精品国产一区二区电影| 国产成人精品久久二区二区91| 啦啦啦韩国在线观看视频| 久久久国产欧美日韩av| 国产精品爽爽va在线观看网站 | 亚洲av日韩精品久久久久久密| 人人澡人人妻人| 久久婷婷人人爽人人干人人爱 | 村上凉子中文字幕在线| 热99re8久久精品国产| 久久婷婷成人综合色麻豆| 欧美成人一区二区免费高清观看 | 国产99白浆流出| 禁无遮挡网站| 亚洲电影在线观看av| 久久精品91无色码中文字幕| 免费无遮挡裸体视频| 久热这里只有精品99| 亚洲无线在线观看| 精品久久久久久成人av| 美女免费视频网站| 三级毛片av免费| 乱人伦中国视频| av片东京热男人的天堂| 午夜精品国产一区二区电影| 禁无遮挡网站| 亚洲一区二区三区色噜噜| 母亲3免费完整高清在线观看| 亚洲欧美日韩另类电影网站| 午夜福利在线观看吧| 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 精品不卡国产一区二区三区| 男女下面插进去视频免费观看| 搞女人的毛片| 高清黄色对白视频在线免费看| 亚洲va日本ⅴa欧美va伊人久久| 日日干狠狠操夜夜爽| 岛国在线观看网站| 一个人观看的视频www高清免费观看 | 亚洲国产看品久久| 99热只有精品国产| 国产精品一区二区在线不卡| 少妇被粗大的猛进出69影院| 国产精品久久电影中文字幕| 巨乳人妻的诱惑在线观看| 久久久久久亚洲精品国产蜜桃av| 乱人伦中国视频| x7x7x7水蜜桃| 久久人人97超碰香蕉20202| 国产国语露脸激情在线看| 999久久久国产精品视频| 国产熟女午夜一区二区三区| 欧美日韩黄片免| 男人的好看免费观看在线视频 | 国产成人精品久久二区二区91| 亚洲最大成人中文| 国产伦人伦偷精品视频| 国产亚洲精品av在线| 国产精品爽爽va在线观看网站 | 老司机午夜福利在线观看视频| 无遮挡黄片免费观看| 性欧美人与动物交配| 色在线成人网| 亚洲一区高清亚洲精品| 欧美日本亚洲视频在线播放| 麻豆av在线久日| 久久 成人 亚洲| 国产一卡二卡三卡精品| 午夜福利成人在线免费观看| 少妇粗大呻吟视频| 亚洲国产欧美一区二区综合| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 757午夜福利合集在线观看| 国产伦人伦偷精品视频| 好看av亚洲va欧美ⅴa在| 国产av一区二区精品久久| 亚洲国产精品久久男人天堂| 男女下面进入的视频免费午夜 | 99国产综合亚洲精品| 91在线观看av| 久久精品人人爽人人爽视色| 麻豆一二三区av精品| 亚洲狠狠婷婷综合久久图片| 99国产精品99久久久久| 国产精品爽爽va在线观看网站 | 精品一区二区三区四区五区乱码| 韩国精品一区二区三区| 老司机午夜十八禁免费视频| 色在线成人网| 欧美日韩亚洲国产一区二区在线观看| 男人舔女人下体高潮全视频| 欧美日韩亚洲国产一区二区在线观看| 黑人巨大精品欧美一区二区蜜桃| 国内精品久久久久久久电影| av天堂在线播放| 性欧美人与动物交配| 日本欧美视频一区| 亚洲国产中文字幕在线视频| 男女做爰动态图高潮gif福利片 | 女性被躁到高潮视频| 两个人视频免费观看高清| 老司机深夜福利视频在线观看| 天天添夜夜摸| 一本大道久久a久久精品| 日日干狠狠操夜夜爽| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利影视在线免费观看| av在线天堂中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久精品吃奶| 国产精品免费视频内射| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 操出白浆在线播放| 精品久久久久久久毛片微露脸| 免费观看人在逋| 又黄又爽又免费观看的视频| 亚洲国产欧美一区二区综合| 亚洲精品在线美女| 国产色视频综合| 免费在线观看影片大全网站| 给我免费播放毛片高清在线观看| 久久香蕉精品热| 日韩av在线大香蕉| 黄色成人免费大全| 高清黄色对白视频在线免费看| 99久久99久久久精品蜜桃| 少妇裸体淫交视频免费看高清 | 熟妇人妻久久中文字幕3abv| 制服丝袜大香蕉在线| 熟妇人妻久久中文字幕3abv| cao死你这个sao货| 在线观看免费视频日本深夜| 亚洲五月婷婷丁香| 动漫黄色视频在线观看| 午夜福利视频1000在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 日韩国内少妇激情av| 色播亚洲综合网| 欧美绝顶高潮抽搐喷水| 嫩草影院精品99| 好看av亚洲va欧美ⅴa在| 日日干狠狠操夜夜爽| 曰老女人黄片| 久久精品成人免费网站| 中文亚洲av片在线观看爽| tocl精华| 国产三级在线视频| 免费久久久久久久精品成人欧美视频| 欧美老熟妇乱子伦牲交| 国产精品一区二区免费欧美| 成人免费观看视频高清|