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

    Multi-physics Coupling of Hydraulic System

    2013-12-07 07:33:19ZHANGJianLUONianningJIANGJihai
    機床與液壓 2013年1期
    關(guān)鍵詞:氣穴哈爾濱工業(yè)大學(xué)工程學(xué)院

    ZHANG Jian, LUO Nianning, JIANG Jihai*

    1.School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150080, China;2.School of Automobile Engineering, Harbin Institute of Technology, Weihai 264209, China

    Multi-physicsCouplingofHydraulicSystem

    ZHANG Jian1, LUO Nianning2, JIANG Jihai1*

    1.SchoolofMechatronicsEngineering,HarbinInstituteofTechnology,Harbin150080,China;2.SchoolofAutomobileEngineering,HarbinInstituteofTechnology,Weihai264209,China

    Inthispaper,asummaryofmulti-physicscouplingofhydraulicsystemisgiven.Thepaperintroducesthetaxonomyofvariouscouplingrelationshipandthemainlyphysicsfieldsinthehydraulicsystem.Andthispaperintroducesasummaryofmulti-physicscouplingofhydraulicsystembothwithinChinaandabroad.Accordingtosomecommonproblemsofthermal,cavitation,pressurerippleandflowrippleonhydraulicsystem,thispaperanalyzeswhichfieldsinfluenceontheseproblems,andbrieflyintroducestheexistingmulti-physicscouplingphenomenaofproblems.Thispaperpredictstheresearchtrendoftheseproblems,andpointsoutthatthemathematicalmodelofmulti-physicscouplingneedstobeestablishedinthefutureresearchwork.studiesonthemainlyphysicsfieldsareinvolvedinhydraulicsystemwhentheseproblemsinfluencehydraulicsystem.Finally,thispaperputsforwardtheresearchmethodofmulti-physicscouplingofhydraulicsystem.

    hydraulicsystem,thermal,cavitation,pressureripple,flowripple,multi-physicscoupling

    1.Introduction

    Multi-physics coupling refers to the interaction of two or more than two fields’ formation of physical phenomena[1], which exists widely in nature and engineering application. Physics holds that field is interaction of substances[2]. However, the fields in hydraulic system cite the conception of physics, mainly refer to flow field, pressure field, temperature field, noise field, structure field, and so on. Multi-physics coupling of hydraulic system mainly refers to the relationship of interaction of different kinds of fields.

    In 1648, Pascal put forward the principle of hydrostatic transmission, but this principle was first applied to hydraulic crane by England until 1850[3]. After this, hydraulic transmission has been widely applied in various machinery equipments, because of these series of advantages. The advantages of hydraulic transmission include high power density, components flexible disposal, convenient control, good dynamic performance, and so on. With the development of science and technology, hydraulic system constantly develops towards high pressure, high velocity and high power which put forward higher requirements for structure, friction, vibration, noise, seal and thermal effect of hydraulic system and hydraulic components which constitute the hydraulic system. However, because the coupled relationship of interaction of several physical factors which are flow field, temperature field, noise field, pressure field and structure field of hydraulic system is complicated, it is difficult to make breakthrough on single performance of hydraulic component or hydraulic system for us, which seriously restricts the application range of hydraulic transmission technology.

    This paper mainly analyzes the properties of kinds of fields of hydraulic system and interaction relationship between these fields, and makes classification of multi-physics coupling from several perspectives. And this paper introduces a summary of multi-physics coupling of hydraulic system both within China and abroad. Through the theory analysis, this paper reveals the action characters and the interactional relationship of physical fields of hydraulic system. Moreover, this paper seeks the ways and the methods to solve the present problem of hydraulic system, from the view of improving the reliability and the integrated transmission efficiency of hydraulic system.

    2.Classification of multi-physics coupling and fields of hydraulic system

    There are many classification methods of multi-physics coupling. According to the region of coupling occurrence, multi-physics coupling is divided into boundary coupling and domain coupling. According to the interaction of coupling, multi-physics coupling is divided into two-way coupling and one-way coupling. According to the action pathway of coupling, multi-physics coupling is divided into direct coupling and indirect coupling. According to the form of coupling equation, multi-physics coupling is divided into differential coupling and algebra coupling. According to the occurrence mechanism of coupling, multi-physics coupling is divided into source coupling, fluid coupling, attribute coupling and geometry coupling[4]. According to the coupling relationship, multi-physics coupling is divided into one-to-one coupling, one-to-many coupling, many-to-one coupling and many-to-many coupling.

    Multi-physics coupling problem of hydraulic system mainly involves in the coupling of structure field, flow field, noise field and temperature field. There are kinds of coupling relationship between related fields. Moreover, with the development of the electronic technology, the trend of electronic control of hydraulic component is more and more obvious in recent years. So based on above several fields, electromagnetic field is also introduced into hydraulic system. And because of the electromagnetic field, the analysis of multi-physics coupling problem of hydraulic system becomes more and more complex.

    People mainly concern system efficiency, system controllability, system noise and reliability of hydraulic system in engineering application. The coupling relationship which has great influence on above factors mainly includes heat-fluid coupling, thermal-structure coupling, fluid-structure interaction, magnetism and solid coupling, electromagnetic coupling, fluid-pressure coupling, solid-pressure coupling, fluid-noise coupling, and so on.

    Heat-fluid coupling: temperature field couples with fluid field.

    Thermal-structure coupling: temperature field couples with structure field.

    Fluid-structure interaction: fluid field couples with structure field.

    Magnetism-solid coupling: magnetism field couples with structure field.

    Electromagnetic coupling: electric field couples with magnetic field.

    Fluid-pressure coupling: fluid field couples with pressure field.

    Solid-pressure coupling: structure field couples with pressure field.

    Fluid-noise coupling: fluid field couples with noise field.

    Temperature field is the most influential physical fields for the performance of hydraulic system. All other physical fields are more or less affected by temperature. Meanwhile, most of the final results of interaction of physical fields in hydraulic system are displayed through heating in the hydraulic system.

    The route loss, local pressure loss, throttle loss, and mechanical loss and volumetric loss of each element in hydraulic system are the main reasons of heating in the hydraulic system[5]. Except for driving load and outputting the useful work, other powers produced by hydraulic system are carried away by the oil, and the temperature of hydraulic system will increase in the form of heat during hydraulic system operation. The route pressure loss is caused by fluid field of oil. The mechanical loss is caused by structure field. Throttle loss, local pressure loss and volumetric loss are not only caused by fluid field but also by structure field. Conversely, the increase of temperature of oil influences directly on oil flow and properties of oil. Meanwhile, the increase of temperature of oil makes the element of hydraulic system to produce volume change, which will influence the tolerance clearance of system, and has an indirect effect on structure field of system. At the same, the system or structure change will directly influence on the fluid field and the pressure field of hydraulic system, and also has an indirect effect on temperature field and noise field.

    3.Review on research status of multi-physics coupling

    The researches on multi-physics coupling of hydraulic system within China and abroad mainly focus on fluid-solid coupling and thermal performance of system at present.

    3.1.Theresearchonthecharacteristicsoffluid-solidcoupling

    Fluid-solid coupling mainly studies the coupling relationship between fluid field and structure field.

    In 1968, Wood[6] considered pipeline as spring-mass system, and he researched on the vibration response of simple pipe system which was carried out under the periodic excitation and sudden valve closure, respectively. In 1977, Walker & Phillips[7] considered the stress waves in the pipe caused by the pressure wave of the pipeline fluid and the poisson coupling between pipelines. Walker & Phillips put forward four governing equations model, then used this model to calculate the reepense of fluid pressure at a valve to sudden valve closure in RPV system. Regetz[8] studied the pressure and flow rate fluctuation of a straight pipe in rocket fuel supply system. And the end of the straight pipe was free closed. Meanwhile, Regetz measured the vibration velocity of the free end of the pipeline through experimental method. Then the measurement results were compared with calculation results in frequency area. The comparison results prove that the influence pipes moving on motion characteristics of fluid is very notabale. In 1980, El1is[9] studied the coupling between valve and joint part of bifurcated pipe in a practical piping. Ellis considered the influence of the pressure wave of fluid, axial stress waves and moment wave in the pipeline. And Ellis used characteristic method to solve the motion equation of the coupling. In 1995, Lee et al[10] combined Paidoussis’s kinetic equation of pipeline and Wiggert’s kinetic equation of pipeline, and ignored the effect of poisson coupling. Lee et al obtained the first nonlinear model which describes the fluid-solid coupled motion. In 1998, WANG Zhongmin et al[11] studied the pressure ripple and dynamic response characteristics of pipeline of hydraulic system through establishing fully fluid-solid coupling model. ZHONG-MIN WANG et al drew the conclusion that fluid-solid coupling model is very important in studying the dynamic characteristics of flexible pipe system. In 2010, according to different types of link coupling, A. Ahmadi and A. Keramat[12] studied the fluid-solid coupling, and they extended the conception of link coupling. In 2011, R. Jaiman et al[13] studied the instable fluid-solid coupling phenomenon used in combined interface boundary condition method, and they proposed and applied a kind of weak coupling method to fluid-solid coupling problem so as to overcome potential instability in coupling interface. In 2012, A. Keramat et al[14] studied the viscoelastic fluid-solid coupling function of tube wall when water hammer happened, and they carried out theoretical analysis and numerical simulation on fluid-solid coupling and viscoelastic of axial movable straight pipeline.

    Research on fluid-solid coupling in China include: In 1989, ZHU Geqi et al[15] introduced an approach for the mathematical modeling of piping system with fluid-solid coupling. The approach which was introduced by Zhu Geqi et al combined complex transfer matrix of fluid and complex mode matrix of structure to establish the fluid-solid coupling model. In 1997, FEI Wenping et al[16] set up a fluid-solid coupling general mathematical model of pipeline system with complex boundary conditions from the two angles of that the free body was never taken and the free body was taken, and they made a theoretical investigation. In the same year, XU Mubing et al[17] from Huazhong University of Science and Technology studied the energy flow in fluid-filled cylindrical shells. In 1999, Jiao Zongxia, Hua Qing et al[18] carried out a modal analysis on the fluid-solid coupling vibration in transmission pipeline. Jiao Zongxia et al had successful carried out a node coupling simulation by using the method of twice coordinate transfers, and deduced a high precision fluid-solid coupling vibration model. In 2008, Tao Yuhua et al[19] carried out a fluid-solid coupling simulation of aircraft hydraulic system, and preliminarily studied the pulsation of hydraulic pipeline. In 2011, according to the clip problem of hydraulic spool valve, WANG Anlin et al[20] studied the fluid solid heat couple of slide valve structure. In 2012, ZHANG Huixian et al[21] studied the fluid-solid coupling of pipeline under the hydraulic exciting wave, and they set up a fluid-solid coupling dynamic model.

    3.2.Researchonthethermalcharacteristicsofhydraulicsystem

    Research on the thermal characteristics of hydraulic system mainly refers to research on the relationship of among temperature, pressure, flow rate and structure of hydraulic system.

    In 1996, J A sidders,D G Tilley,P J ChapPle et al[22] from the British university of Bath studied the thermal characteristics of open circuit hydraulic system, set up the thermodynamic model of pump, valve, oil tank, radiator, and so on, by using the first order differential equation which was established based on a series of control volume energy conservation equations. According to different conditions, J A sidders et al simulated the characteristics of system pressure, flow rate and temperature. In 1998, Dominique Legendre et al[23] studied the thermal and dynamics development of spherical bubble in superheated or supercooling fluid. In their research, they assumed that the development of bubble was controlled by heat conduction, and they carried out a direct numerical simulation on the growth and collapse of spherical steam bubble in the fluid. Meanwhile, they obtained the fluid force acted on bubble when there is mean flow. In 2000, the thermodynamic model based on aircraft hydraulic system was developed by Engelhardr J[24]. Then Engelhardr J set up a heat transfer mathematical model of fluid and solid, and combined heat generation of hydraulic oil with heat-transfer of contact solid.

    In 2003, C. Vortmann et al[25] set up a new cavitation model. The phase change process of memory alloy was extended to the phase change of water in this new cavitation model. And the thermodynamic description of cavitation phase change process of water was made by this new cavitation model, where involved density, pressure and temperature of water. In 2006, the thermodynamic model of gear pump based on the first law of thermodynamics was established by Eduardo Dalla Lana et al[26]. And this model considered the power loss when the hydraulic oil passes through the gear pump, and the heat transfer occurs between the hydraulic oil and the ambient air. Then Eduardo Dalla Lana et al obtained a new kind of gear pump efficiency formula. In 2008, the thermodynamic model of vehicle one-dimensional transient power steering hydraulic system based on the first law of thermodynamics was established by Timothv C Scott[27]. Timothv C Scott set up the thermodynamic model of ordinary hydraulic component, pipeline and radiator. Then Timothv C Scott considered the effect of the distance between radiators on cooling air distribution in thermal module by using one-dimensional method, and used system model to simulate system temperature characteristics. In 2010, Maria Grazia De Giorgi et al[28] studied the effect of temperature and heat transfer on cavitation in vapor-liquid conversion process when the cavitation occurs in fluid. Maria Grazia De Giorgi et al mainly studied the water and cryogenic fluid, and concluded under the condition of that the same cavitation number cavitation is increased with temperature increasing. In 2011, D. T. Frate et al[29] studied the flow and heat transfer in the hydraulic reservoir of thrust vector control system. They obtained that for the laminar flows, when the inlet temperature of the working fluid is maintained at a lowlevel, the buoyancy effects at both sides of the reservoir walls are not dominant in the flows. However, at a higher inlet temperature of the working fluid, this greatly changes both flow patterns and heat transfer performance of the working fluid and the cooling air around the reservoir due to the stronger buoyancy effects.

    In China: In 1994, LU Tinghai[30] studied the effect of temperature on main performance parameters of hydraulic pump, and set up a simple mathematical model for the relationship between some of performance parameters of hydraulic pump and temperature, then pointed out that only simultaneous studing the pressure and temperature characteristics of hydraulic pump, could we comparatively detailed assess the performance of pump. In 2001, Yang Yung-Kuang and Jeng Ming-Chang[31] from Taiwan studied the thermal effect of asymmetrical hydraulic servo cylinder. Yang Yung-Kuang and Jeng Ming-Chang analyzed the influence of eccentricity and misalignment factor between hydraulic bearing and piston rod using numerical method. In 2006, according to the thermodynamic properties of axial piston pump, LI Chenggong et al[32] set up a thermodynamics model of aviation piston pump, and carried out a detail heat transfer analysis for internal structure of piston pump. In 2008, LI Chenggong et al[33] put forward the basic function of establishing simple model of thermal-hydraulic elements, and introduced the association rule and method of thermal- hydraulic system model. With the above methods, it can fully automatic generate the thermal-hydraulic system. In 2009, based on energy conservation principle, LI Yonglin et al[34] deduced the calculation formula for temperature variation in control volume, and combined with the pressure flow characteristic of servo valve, they put forward the method of establishing thermodynamics model of servo valve. Meanwhile, LI Yonglin et al carried out some modelings and simulations for hydraulic system which includes a four-through slide-valve. In 2010, LI Yonglin et al[35] applied the control volume method to set up the thermodynamics model of hydraulic poppet valve. LI Yonglin et al considered the effect of the change of thermophysical characteristics of hydraulic oil on the pressure loss of valve in model, then carried out a simulation study on hydraulic system which includes a poppet valve. In 2011, XI Renguo et al[36] studied the thermal analysis of aircraft hydraulic system. They used steady analysis method to estimate the equilibrium temperature of hydraulic system, used transient analysis to predict the law of system temperature change, and introduced a method of neural network based on analyzing.

    4.Thermal field analysis of hydraulic system

    The thermal-fluid coupling and the thermal-structure coupling are eqaul to temperature field couples with fliud field and structure field.

    The hydraulic system heating is the concentrated expression of interaction of multi-physics fields of hydraulic system which includes fluid field, pressure field, structure field and temperature field. Setting the throttle loss caused by various hydraulic valves as an example, when the hydraulic oil flows through hydraulic valve, for the opening size of hydraulic valve compares with hydraulic pipeline has a big change, it would cause the sudden change of oil flow velocity, and then lead to the sudden change of oil pressure, which causes the power loss of hydraulic system and the system heating. Meanwhile, the sudden change of pressure causes structural vibration, and produces noise.

    4.1.Basicthermodynamicsmodelofhydraulic

    system

    To evaluate the impact relation of physics fields in hydraulic system which include thermal field with fluid field, pressure field and temperature field, this paper analyses a typical thermodynamics model of hydraulic system which on the law of conservation of energy.

    This paper chooses the control volume as shown in Fig.1 for one-dimensional flow fluid.

    Fig.1 Control volume

    The energy conservation equation of control volume for

    (1)

    If the kinetic and potential energies are neglected, the time rate of change of the energy can be expressed according to

    (2)

    where,uis the specific internal energy of control volume.

    The time derivative ofhcan be expressed as

    (3)

    where,αpis the cubical expansion coefficient,cpis the specific heat at constant pressure.cpis expressed as

    (4)

    The specific enthalpy is defined as

    h=u+pv

    (5)

    After introducing Eq. (3) and (5) into Eq. (2), it is obtained that

    (6)

    where,Vis the fluid volume.

    The continuity equation for one-dimensional flow is given as

    (7)

    Combining Eq. (6), (7) and (1), it is obtained that

    (8)

    (9)

    (10)

    Introducing Eq. (9), (10) into Eq. (8), it is given

    (11)

    Assuming that the average enthalpy within the control volume equates to the leaving enthalpy regardless of the inlet conditions, Eq. (11) can be expressed as

    (12)

    The change in specific enthalpy within the control volume is related to the change in pressure and temperature by

    (13)

    According to above thermodynamics model, it can be gotten that the thermal field of hydraulic system has a close relation with temperature field, pressure field and fluid field of hydraulic system. The change of temperature, pressure and fluid flowing situation can cause the change of thermal field in hydraulic system. The change of thermal field of hydraulic system is a kind of multi-physics coupling problem.

    This mathematical model is the classic one-dimensional flow thermodynamics model at present. More perfect multi-dimensional flow thermodynamics model will be described in the future research.

    4.2.Asummaryofthermalfieldcoupleswithotherphysicsfieldsphenomenoninhydraulicsystem

    The hydraulic system heating will cause physical change in hydraulic oil. The viscosity of hydraulic oil is decreased as the temperature increases. The viscosity decrease of hydraulic oil causes the leakage increase of hydraulic system, and results in the reduction of the volume efficiency of system. Meanwhile, the increase of oil temperature also affects the Reynolds number of oil. The temperature increases causes the increase of the Reynolds number of oil and the flow instability of oil increases, easy leads to turbulent, and induces greatly pressure ripple. The process of change in physical propertics of oil is caused by hydraulic system heating including thermal field, fluid field and pressure field. The coupling relationships in these fields including thermal field causes the change of fluid field, then the change of fluid field causes the pressure field ripple. Conversely, the changes of fluid field and pressure field will decrease the efficiency of hydraulic system, cause the power loss and the system heating, and make the continual system temperature rise.

    There is another impact of the change of thermal field of hydraulic system on physical properties of hydraulic oil. That is the bulk modulus of elasticity of hydraulic oil. The bulk modulus of elasticity of hydraulic oil is an important physical parameter for hydraulic oil. It has significant effect on position precision, power level, response time and stability of hydraulic system[37-38]. The influencing factors of bulk modulus of elasticity of hydraulic oil mainly include the pressure, temperature and gas content of hydraulic oil. The effect of thermal field on bulk modulus of elasticity of hydraulic oil mainly includes two respects. On the one hand, the temperature increase causes the decrease of oil viscosity, then reduces the bulk modulus of elasticity of hydraulic oil. On the other hand, the temperature increase causes the change of gas content of hydraulic oil, then impacts the bulk modulus of elasticity of hydraulic oil. The effects of the change of thermal field on gas content and bulk modulus of elasticity of hydraulic oil are expressed as follows:the gas solubility decreases with increasing the temperature under the same gaseous phase partial pressure, and which causes the increase of air bubble, then decreases the effective bulk modulus of elasticity of hydraulic oil. The decrease of effective bulk modulus of elasticity of hydraulic oil causes the decrease of the efficiency of hydraulic system. And the loss of power of hydraulic system will transfer into the quantity of heat. This causes the increase of the system heating. The increase of the system heating makes the temperature of hydraulic system continue increase, causes the dissolved gas continue to precipitate, and which impacts the flowability of hydraulic oil in hydraulic system. Meanwhile, the increase of the system heating causes the change of pressure of hydraulic oil along the way. This also shows the characteristics of coupling of multi-physics of hydraulic system.

    The thermal-structure coupling which caused by the change of thermal field displays that the temperature field forms the temperature difference, and leads to the expansion or reduction of the structure, which generates thermal stress[39]. The effect of temperature on element structure of hydraulic system is manifested as the element of hydraulic system expansion deformation caused by hydraulic system superheat, and it decreases the tolerance clearance between different parts of the kinematic pair, then causes the increase of frictional resistance, even the clamping stagnation, and eventually the hydraulic control element failure. Because of the structural deformation, it causes the increase of frictional resistance of kinematic pair, which will cause the power loss, then cause the increase of hydraulic system heating.

    Therefore, the change of thermal field of hydraulic system is not only manifested as increasing or decreasing the hydraulic system temperature, but can also produce the influence on fluid field, structure field and pressure field of hydraulic system, meanwhile, fluid field, structure field and pressure field react on temperature field. Hydraulic system heating phenomenon is a multi-physics coupling problem which involves temperature field, fluid field, structure field, pressure field, and so on.

    4.3.Theforecastofresearchonthermalproblemofhydraulicsystem

    For the research on thermal problem of hydraulic system, researchers should consider this problem from the view of multi-physics coupling. Firstly, researchers should make a thorough study on the heating mechanism of hydraulic system, optimize the existing research results, through comprehensive consideration of the coupling effect among multi-physics fields, find the coupling relationship among the multi-physics fields effecting on hydraulic system heating; Secondly, researchers should set up the mathematical model of thermal fields and other physics fields, and carry out simulation and experiment researches on establishing and optimizing the mathematical model, so as to achieve twice the result with half the effort for research on thermal of hydraulic system.

    5.The multi-physics coupling of cavitation of hydraulic system

    5.1.Briefintroductiononcapitationphenomenon

    Cavitation was defined earliest by Knapp et al[40]. The solution behavior of gases in hydraulic oil follows Henry law, that the volume of gases which dissolve in hydraulic oil is proportional to the absolute pressure of hydraulic oil under certain temperature, and increases with the increase of time. When the hydraulic system appears local low pressure, the volume of original micro-bubble dropping in hydraulic oil expands continuously, and these micro-bubbles polymerize each other, finally form certain volume bubbles and dissociate out. In this case, slight cavitation is formed. If the pressure continuously decreases to lower than the air apart pressure, the gases dissolving in hydraulic oil would separate out the hydraulic oil rapidly and form numerous bubbles. In this case, serious cavitation is formed. If the pressure continuously decreases to lower than the saturated vapor pressure of hydraulic oil, the hydraulic oil would vaporize and form numerous bubbles. In this case, intense cavition is formed. Cavitation often occurs the pressure of hydraulic oil lower than vapor pressure of hydraulic oil under current thermodynamics conditions[41]. The formation of low pressure usually relates to the pressure ripple of hydraulic system caused by flow ripple of hydraulic system.

    According to the research results of cavitation problem, it is one of the necessary conditions of hydraulic oil generating cavitation that there is micro gas nucleus in hydraulic oil. In order that the hydraulic oil does not generate cavitation, the free gas nucleus and the circle flowing liquid should satisfy the equilibrium condition. Besides satisfying the force equilbrium condition, the free gas nucleus and the circle flowing liquid need to meet the thermal equilibrium condition and the phase equilibrium condition. Namely when gas-liquid two-phase coexists in the equilibrium state, the second law of thermodynamics requires that both of the temperatue and the chemical potential must keep continuity at phase interface. The process of cavitation has intense thermodynamics imbalance[42]. It is showed that the generation and development of cavitation relate to multi-physics.

    When the cavitation occurs, the produced bubbles flow to the high pressure region with hydraulic oil, and the bubbles are compressed by the shock of ambient high pressure oil, the volume of bubble reduction rapidly then collapse, finally bubbles condense fluid. In this case, the space occupied by original bubbles forms a vacuum, then the nearby high pressure oil fills the space with a high flow speed. Because the collapse of bubble occurs in transient, it will cause a strong pressure shock, and both the pressure and the temperature increase drastically at the occurrence station of the bubble collapse. According to literature[43], the continuous bubble collapse will cause the increase of pressure up to 1 000 MPa. The pressure rised instantly spreads as pressure waves, and it makes the hydraulic system produce a strong vibration and noise. Partial pressure energy can be transformed into thermal energy by particle oscillation, and it makes the temperature of partial high pressure region be up to more than 1 000℃. The structure surface of hydraulic system at the bubbles accumulative station endures the hydraulic impact and the high temperature. The erosion and spalling appears, and this phenomenon is named cavitation erosion.

    People are not willing to see the phenomenon of cavitation in hydraulic system. The cavitation will cause the decrease of flow rate in hydraulic system, the increase of low pressure in pump, asymmetric load, vibration, noise, cavitation erosion, and so on[44]. It is the result of multi-physics coupling effect for the occurrence, development and collapse of cavitation. The cavitation relates to fluid field, pressure field, temperature field, noise field and structure field.

    5.2.Bubbledynamicequation

    The micro gas nucleus in hydraulic oil is one of the necessary preconditions of hydraulic oil generating cavitation. When free gas nucleus which contains gas and vapor, circle flowing liquid keep a dynamic balance, the force equilibrium equation is

    p=pv+pg-2S/R

    (14)

    where,pis the pressure of liquid aroud gas nucleus (Pa),pvis the vapor pressure (Pa),pgis the gas pressure inside the bubble;Sis surface tension of liquid,Ris gas nucleus radius under equilibrium state (m).

    According to Eq. (14), the condition of gas nucleus expansion is that the right side of the equation is larger than the left side of the equation, which means that the pressure inside gas nucleus is larger than the pressure of liquid aroud gas nucleus. When the balance is broken, the bubble dynamic equation is written as the Rayleigh-Plesset equation[45] as

    (15)

    where,Rbis the bubble radius (m),Sis the surface tension of liquid,Pis the pressure of liquid aroud bubble (Pa),vlis the kinematic viscosity of liquid;ρlis the liquid density (kg/m3),pbis the bubble surface pressure.pbis expressed as

    (16)

    whereTcis the temperature inside the bubble (K),T∞is the temperature of liquid (K),pvis the vapor pressure (Pa),pg0is the initial gas pressure;Ris the initial gas radius (m),Rbis the bubble radius (m).

    Assuming that the behavior of the gas in the bubble is polytropic,

    (17)

    whereγis the polytropic exponent,pgis the gas pressure at arbitrary time.

    According to Eq. (15)~(17), the occurrence and development of cavitation in hydraulic oil have intimate relationship with the pressure, temperature, viscosity and density of hydraulic oil. Meanwhile, the gas pressure decreases gradually with the development of bubble, but when the volume of bubble decreases rapidly, the pressure inside bubble will increase greatly. So it will produce very high pressure at the centre of collapse when the bubble collapses, then cause vibration and noise in hydraulic system.

    5.3.Asummaryofthephenomenonofmulti-physicscouplingincavitation

    After the occurrence of cavitation, the bubbles collapse. Besides producing vibration and noise, because the hydraulic oil around bubbles occupies the space,where formed bubbles collapse when bubbles collapse, it will cause the continuity of fluid to be broken, decrease the ability of through oil of oil suction pipe, then decrease the ability of through oil of hydraulic system.

    There is gas nucleus at the hydraulic system structure solid surface. The surface roughness of solid surface is a key factor for the influence on the size and quantity of gas nucleus. Besides, the solid surface irregularity has an effect on fluid flow, then produces the change of pressure. Commonly, the higher surface roughness is, the easier occurrence of cavitation will be[46]. The occurrence of cavitation is influenced by the structure surface quality. The structure surface quality has an effect on the fluid field and pressure field of hydraulic system. When cavitation occurs, there are the effects among structure field, fluid field and pressure field.

    The cavitation erosion is the cavitation effect on the hydraulic system structure surface. The occurrence of cavitation erosion not only damages the hydraulic system structure surface and decreases the service life of hydraulic element, but also produces the flow rate ripple and the pressure ripple in hydraulic pipe, then causes the power loss, leads to the system heating. So there are interactions among structure field, fluid field, pressure field and temperature field in cavitation erosion.

    5.4.Theforecastofresearchoncavitation

    According to the analysis of cavitation in hydraulic system, it can be known that the cavitation relates to several physics fields, but not an independent action by a certain field. There are multi-physics interactions among temperature field, pressure field, fluid field, structure field and noise field in cavitation. The occurrence and development of cavitation are typical multi-physics coupling problems.

    In the past study on cavitation, researchers often ignore the multi-physics coupling problem, and more studies of the cavitation problem is from a single angle, so research results mostly have limitation. So, according to the multi-physics coupling characteristic of cavitation, people study cavitation should apply multi-physics coupling ideological which contains thermal, pressure and fluid fields, and so on. It is need to inquiry into the change of thermal, pressure and oil flow condition effect on cavitation, and more deeply seek for the mechanism of occurrence, development and disappearance of cavitation. Specially, it is a trend that when the mathematical models of occurrence, development and disappearance of cavitation are built, researchers should consider more effect factors to model, and get more detailed mathematical models.

    6.Multi-physics coupling of hydraulic system pressure and flow rate ripple

    The hydraulic pump output pressure and the flow rate are not absolute stable in the hydraulic system practical work. Because of pump’s structural characteristic, the pump output flow rate changes periodically, and it is the inherent attribute of pump. The periodic flow rate of pump causes the pressure ripple of pump outlet and pipeline, and spreads to the whole system.

    Taking the axial piston pump as an example, as Fig.2 shows, the displacement of axial piston pumpVis

    (18)

    where,dis the piston diameter (m),lis the piston stroke (m),Ris the pitch circle diameter of piston centre (m),γis the theswashplate angle (°),Zis piston number.

    The flow rate of axial piston pump is

    (19)

    where,nBis the rotation speed of pump (r/min),ηwis the volume efficiency of pump.

    When the cylinder block rotatesφangle, the displacement of piston ball-end centre relative to cylinder block centre is

    s=R(1-cosφ)tanγ

    (20)

    where,φis the rotated angle of piston relative to vertical centre line (°).

    The piston axial movement is expressed as

    (21)

    So, the instantaneous flow rate of every pistonQ′ is expressed as

    (22)

    The instantaneous flow rate of pumpQis expressed as

    (23)

    According to equations (18)~(23), the instantaneous flow rate of piston pump is ripple and changes periodically. So with the periodical change of the flow rate, the pressure of pump changes periodically.

    Fig.2 The nomogram of instantaneous flow on piston-pump

    The hydraulic pump is the main noise source in hydraulic system. The noise of hydraulic pump can be classified into fluid-borne noise and structure-borne noise. The fluid-borne noise includes the flow rate ripple, the pressure shock and the cavitation noise. The same, taking the axial piston pump as an example, the key root cause of the fluid-borne noise is the flow rate ripple which then causes the pressure ripple[47]. And the fluid-borne noise will be transferred to the whole hydraulic system by hydraulic pipeline. Study on the generation mechanism of the fluid-borne noise of piston pump needs to begin with the flow rate ripple of piston pump, then study on pressure ripple, and they are complementing each other.

    Because of the flow rate ripple of piston pump, when there is no enough hydraulic oil in piston chamber, and which causes the pressure of piston chamber is lower than the oil vapor pressure, the piston pump will generate cavitation[48-49]. The occurrence of cavitation will bring a series of multi-physics coupling problems as noted in the previous section for hydraulic system.

    Due to the change of the flow velocity, the flow hydraulic oil in hydraulic system causes the change of pressure, and it sometimes forms a local low pressure region in hydraulic system. When the pressure of local low pressure region is lower than air apart pressure of gases which is dissolved into hydraulic oil, the dissolved gases will be isolated from hydraulic oil, then generates cavitation. It forms the complicated multi-physics coupling problem and produces a series of adverse effects on hydraulic system.

    When pump output fluid units which pipeline, hydraulic valves, hydraulic cylinder, hydraulic motors, if the pipeline or the element has great stiffness, great oscillating friction force and zero defect, the pressure ripple would be gradually damped and tend to stable the pressure flow; when the frequency of pressure oscillations of pulsating fluid is consistent or close to the inherent frequency of element or system, the resonance will occur, which makes the system can not work normally, even cause the element damage[50]. This embodies the coupling relationships among fluid field, pressure field and structure field in hydraulic system which flow rate ripple causes pressure ripple, then causes structure vibration.

    Due to the sudden closing of valve, the flow hydraulic oil in pipeline will form a high peak pressure in pipeline. The peak pressure will sometimes enhance with source pressure ripple or sutained oscillation in pipeline. It is the hydraulic shock[50].The hydraulic shock not only causes great vibration, but also can cause the damage of hydraulic system. The hydraulic shock also causes the hydraulic system power loss and the hydraulic system heating. The hydraulic shock relates to multi-physics coupling problem which contains fluid field, pressure field, structure field and temperature field.

    In a word, the pressure ripple and the flow rate ripple of hydraulic system have many-sided effects on hydraulic system working performance, it relates to multi-physics coupling problem which contains fluid field, pressure field, noise field, temperature field and structure field, and is not simple pressure field or fluid field problem.

    The effect pressure ripple and flow rate ripple of hydraulic system on hydraulic system working performance is the process of multi-physics coupling. When the pressure and the flow rate ripple of hydraulic system are studied, researchers apply the multi-physics coupling theory to study the generation mechanism of them, and build the detailed mathematical model of pressure and flow rate ripple couple with other physics fields, then find that the methods and the measures of decreasing pressure and flow rate ripple of hydraulic system is one of the research topics in pressure and flow rate ripple of hydraulic system.

    7.Research method of multi-physics coupling problem

    At present, little investigation has been made on multi-physics coupling characteristics of hydraulic system. Because the performance of hydraulic system is required to be increasingly improved, people should develop to study the multi-physics coupling characteristics of hydraulic system. The multi-physics coupling problem of hydraulic system is a complicate physics phenomenon, and it relates to multi-physics fields interaction. The multi-physics fields contain temperature field, pressure field, noise field, fluid field and structure field. The action mechanism between physics fields is complex. Based on the complicated of multi-physics coupling of hydraulic system, studying multi-physics coupling should concentrate on the following aspects:

    1) Set up mathematical model. According to thermal, cavitation, pressure ripple and flow rate ripple, these problems involved multi-physics coupling characteristics are theoretically analyzed. Researchers set up the interaction of parameters relationship between physics fields, establish and optimize the mathematical model of multi-physics coupling, then understand the interaction mechanism of physics fields.

    2) Numerical simulation. With the development of computer technology, numerical simulation has become a very practical research method. Researchers apply numerical simulation technology and combine with specific problem to carry out the simulation on multi-physics coupling characteristics of hydraulic system, verify the correctness of mathematical model, find the interaction characteristics between physics fields, and get the optimal solution of multi-physics characteristics through improving the simulation parameters.

    3) Experimental research. According to the numerical simulation results of multi-physics coupling characteristics of hydraulic system, researchers design and build corresponding experimental platform, carry out verification test on the numerical simulation, then get compellent research results.

    4) Design theory and design method of multi-physics coupling of hydraulic system. Researchers improve the design level of hydraulic system through theoretical analysis and summarizing experimental research, and find design theory and design method of hydraulic system suitable for the requirement of modern industry.

    8.Conclusions

    The working process of hydraulic system relates to fluid field, pressure field, noise field, temperature field, structure field, and so on. Through the analyses of thermal, cavitation, pressure ripple and flow rate ripple, this paper can obtain that the working process of hydraulic system not only relates to single physics field, related physics fields not isolated effect on working performance of hydraulic system, between fields interaction and mutual influence. It is a trend that people research on and design the hydraulic system in the view of multi-physics coupling.

    Based on that the multi-physics coupling characteristics of hydraulic system has important effect on performance of hydraulic system, this paper predicts the research direction on thermal, cavitation, pressure ripple and flow rate ripple. The research direction includes studying the interaction relationship between physics fields, setting up the mathematical model of coupling relationship between physics fields, carrying out the simulation and experimental verification on mathematical model, finally finding methods and measures to optimize the performance of hydraulic system, and improving the reliability and the integrated transmission efficiency of hydraulic system.

    [1] Carlos A,Felippa K C.Partitioned analysis of coupled mechanical systems[J].Computer Methods in Applied Mechanics and Engineering,2001,190(24/25):3247-3270.

    [2] ZHANG Qiuguang.Field Theory[M].Beijing: Geological Publishing House,1988.

    [3] Burrows C.Fluid power systems-some research issues[J].PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS PART C-JOURNAL OF MECHANICAL ENGINEERING SCIENCE,2000,214(1):203-220.

    [4] SONG Shaoyun.Modeling of Multiphysics Problem and Research of Coupling Relation[J].Journal of Wuhan Polytechnic University,2008,27(3):46-49.

    [5] FAN Xiaowang,DING Guolong,ZHANG Yong,et al.Discussion on Hydraulic System Heating[J].ENGINEERING & TEST,2011,51(1):65-67.

    [6] Wood D J.A Study of the Response of Coupled Liquid Flow-structural Systems Subjected to Periodic Disturbances[J].ASME Journal of Basic Engineering,1968,90:532-540.

    [7] Walker J S,Phillips J W.Pulse Propagation in Fluid-filled Tubes[J].Journal of Applied Mechanics,1977,44:31-35.

    [8] Regetz J D.An Experimental Determination of the Dynamic Response of a Long Hydraulic Line[Z].Washington: Natonal Aeronautics and Space Administration,Technical Note:D-576.

    [9] Ellis J.A Study of Pipe-liquid Interaction Following Pump-trip and Check-Valvel.JClosurein a Pumping Station.In Proceedings of the 3th Intenrational Conference on Pressure Surges[Z].BHRA,Canterbury,U.K.March.1980:203-220.

    [10] Lee V,Pak C H,Hong S C.The Dynamic of a Piping System with Internal Unsteady Flow[J].Jounral of Sound and Vibration,1995,182:297-311.

    [11] ZHONG-MIN WANG,VIBRATION S T.PRESSURE FLUCTUATION IN A FLEXIBLE HYDRAULIC POWER SYSTEM ON AN AIRCRAFT[J].Computers & Fluids,1998,27(1):1-9.

    [12] Ahmadi A,Keramat A.Investigation of fluid-structure interaction with various types of junction coupling[J].Journal of Fluids and Structures,2010,26(7/8):1123-1141.

    [13] Jaiman R,Geubelle P,Loth E,et al.Combined interface boundary condition method for unsteady fluid-structure interaction[J].Computer Methods in Applied Mechanics and Engineering,2011,200(1-4):27-39.

    [14] Keramat A,Tijsseling A S,Hou Q,et al.Fluid-structure interaction with pipe-wall viscoelasticity during water hammer[J].Journal of Fluids and Structures,2012,28:434-455.

    [15] ZHU Geqi,CAI Yigang,YANG Shichao.An Approach to Mathematical Modeling of Piping System With Fluid-Structure Interaction[J].Journal of Hydrodynamics,1989,4(1):6-12.

    [16] FEI Wenping,YANG Jiandong,SUN Lihua.Analysis Method of Complex Boundary Condition of Fluid structure Interaction in Piping Systems[J].J Wuhan Univ of Hydr & Elec Eng,1997,30(6):11-15.

    [17] XU Mubing,ZHANG Xiaoming,ZHANG Weiheng.The Vibrational Energy Flow in a Cyl indrical Shell Filled with Fluid[J].J.Huazhong Univ of Sci & Tech,1997,25(2):86-88.

    [18] JIAO Zongxia,HUA Qing,YU Kai.FREQUENCY DOMAIN ANALYSIS OF VIBRATIONS IN LIQUID FILLED PIPING SYSTEMS[J].ACTA AERONAUT ICA ET ASTRONAU TICA SINICA,1999,20(4):29-33.

    [19] TAO Yuhua,HUANG You,ZOU Tao.The Simulation of Fluid Structure Interaction and Research on Stress Fluction of Aircraft Hydraulic System[J].Machine Tool & Hydraulics,2008,36(10):161-162.

    [20] WANG Anlin,DONG Yaning,ZHOU Pengju,et al.Robust Design Method for the Seizure Problem of Hydraulic Slide Valve[J].JOURNAL OF SHANGHAI JIAOTONG UNIVERSITY,2011,45(11):1637-1642.

    [21] ZHANG Huixian,OU Ziming,WU Juan,et al.Dynamic Modeling and Experiments for Fluid Structure Interaction in Pipe under HydraulicShock Wave[J].JOURNAL OF XI’AN JIAOTONG UNIVERSITY,2012,46(3).

    [22] Sidders J A,Tilley D G,Chapple P J.Thermal-hydraulic performance prediction in fluid power systems[J].PROCEEDINGS OF THE INSTITUTION OF MECHANICAL ENGINEERS PART I-JOURNAL OF SYSTEMS AND CONTROL ENGINEERING,1996,210(4):231-242.

    [23] Legendre D,Borée J,Magnaudet J.Thermal and dynamic evolution of a spherical bubble moving steadily in a superheated or subcooled liquid[J].Physics of Fluids,1998,10(6):1256-1272.

    [24] Joerg Engelhardt.Thermal Simulation of an Aircraft Fluid Power System with Hydraulic-Electrical Power Conversion Units[C]//Proceedings of 1st FPIN-PhD Symp.Germany: FPIN,2000:435-448.

    [25] Vortmann C,Schnerr G H,Seelecke S.Thermodynamic modeling and simulation of cavitating nozzle flow[J].International Journal of Heat and Fluid Flow,2003,24(5):774-783.

    [26] Eduardo Dalla Lana.A New Evaluation Method for Hydraulic Gear Pump Efficiency through Temperature Measurements[C]//SAE 2006-01-3503.

    [27] Timothy C S.Thermal Modeling of Power Steering System Performance[C]//SAE 2008-01-1432.

    [28] De Giorgi M G,Bello D,Ficarella A.Analysis of Thermal Effects in a Cavitating Orifice Using Rayleigh Equation and Experiments[J].Journal of Engineering for Gas Turbines and Power,2010,132(9):92901-92910.

    [29] Frate D T.Flow and Heat Transfer in Hydraulic Reservoir of Thrust Vector Control System[J].JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER,2011,25(1):147-154.

    [30] LU Ting-hai.INFLUENCE OF TEMPERATURE ON MAIN PERFORMANCE PARAMETERS OF HYDRAULIC PUMP[J].CHINESE JOURNAL OF MECHANICAL ENGINEERING,1994(3):60-64.

    [31] YANG Y,JENG M.Analysis of thermal effects on the misaligned hydraulic servo cylinder[J].Tribology International,2001,34(2):95-106.

    [32] LI Cheng-gong,JIAO Zongxia.Thermal-hydraulic Modeling and Simulation of Piston Pump[J].Chinese Journal of Aeronautics,2006,19(4):354-358.

    [33] Chenggong L,Zongxia J.Calculation Method for Thermal-Hydraulic System Simulation[J].Journal of Heat Transfer,2008,130(8):84503-84505.

    [34] LI Yonglin,LI Baorui,SHEN Yanliang,et al.Thermal-hydraulic Modeling and Simulation of Hydraulic Servo Valve[J].Journal of System Simulation,2009,21(2):340-343.

    [35] LI Yonglin,CAO Keqiang,XU Haojun,et al.Thermal-hydraulic Modeling and Experimental Investigation of a Hydraulic Poppet Valve[J].Mechanical Science and Technology for Aerospace Engineering,2010,29(11):1521-1524.

    [36] XI Renguo,LIU Weiguo,CHEN Huanming,et al.Research on Thermal Analysis Method of Aircraft Hydraulic System[J].Machine Tool & Hydraulics,2011,39(23):41-44.

    [37] GEORGE H F,BARBER A.What is bulk modulus and when is it important[J].Hydraulic & Pneumatics,2007,(7):34-39.

    [38] MERRITT H E.Hydraulic control systems[M].New York: John Wiley and Sons,1967.

    [39] ZHANG Qiang.Numerical simulation method of multi-physics coupling for heat exchange tube and fluid[D].Daqing: Chemical Process Equipment of Northeast Petroleum University,2011.

    [40] HE Y,LIU Y.Experimental research into time-frequency characteristics of cavitation noise using wavelet scalogram[J].Applied Acoustics,2011,72(10):721-731.

    [41] Wei Y,Tseng C,Wang G.Turbulence and cavitation models for time-dependent turbulent cavitating flows[J].Acta Mechanica Sinica,2011,27(4):473.

    [42] HEguogeng,LUO Jun,|HUANG Suyi.Influence of Incipient Cavitation[J].JHuazhong Univof Sci & Tech,1999,27(1):67-69.

    [43] Singh R,Tiwari S,Mishra S.Cavitation Erosion in Hydraulic Turbine Components and Mitigation by Coatings: Current Status and Future Needs[J].Journal of Materials Engineering and Performance,2011:1-13.

    [44] Goncalvès E.Numerical study of unsteady turbulent cavitating flows[J].European Journal of Mechanics-B/Fluids,2011,30(1):26-40.

    [45] De Giorgi M G,Bello D,Ficarella A.Analysis of Thermal Effects in a Cavitating Orifice Using Rayleigh Equation and Experiments[J].Journal of Engineering for Gas Turbines and Power,2010,132(9):92901-92910.

    [46] GAO Qiusheng.Further Exploration of Liquid Cavitation Mechanism[J].JOURNAL OF HOHAI UNIVERSITY,1999,27(5):63-67

    [47] YANG Huayong,ZHANG Bin,XU Bing.Development of Axial Piston Pump/motor Technology[J].Chinese Journal of Mechanical Engineering,2008,44(10):1-8.

    [48] Wang S.The Analysis of Cavitation Problems in the Axial Piston Pump[J].Journal of Fluids Engineering,2010,132(7):74502-74506.

    [49] Seeniraj G K,Ivantysynova M.Impact of Valve Plate Design on Noise,Volumetric Efficiency and Control Effort in an Axial Piston Pump[J].ASME Conference Proceedings,2006,2006(47713):77-84.

    [50] QI Renjun.Mechanism Research of Pressure Ripple for Hydraulic Systems[J].JOURNAL OF TONGJI UNIVERSITY,2001,29(9):1017-1022.

    液壓系統(tǒng)多場耦合

    張 健1,羅念寧2,姜繼海1*

    1.哈爾濱工業(yè)大學(xué) 機電工程學(xué)院,哈爾濱 150080;2.哈爾濱工業(yè)大學(xué) 汽車工程學(xué)院,山東 威海 264209

    概述了液壓系統(tǒng)多物理場耦合問題,介紹了各種耦合關(guān)系的分類方法以及液壓系統(tǒng)中所涉及的主要物理場,并介紹了液壓系統(tǒng)多場耦合問題的國內(nèi)外研究現(xiàn)狀。針對液壓系統(tǒng)中常見的熱、氣穴、壓力與流量脈動問題,分析了這些常見問題主要受到了哪些物理場的影響,并簡要介紹了這些問題中存在的多物理場耦合現(xiàn)象,對這些問題的研究趨勢進行了預(yù)測,指出在日后的研究工作中應(yīng)建立多場耦合問題數(shù)學(xué)模型。最后提出了針對液壓系統(tǒng)多場耦合特性的研究方法。

    液壓系統(tǒng);熱;氣穴;壓力脈動;流量脈動;多場耦合

    TH137.1

    2012-08-25

    National Natural Science Key Foundation of China (50875054), Science Key Foundation of Zhejiang University (GZKF-2008003), Natural Science Key Foundation of Shandong(ZR2009FL002)*JIANG Jihai, Professor. E-mail: jjhlxw@hit.edu.cn

    10.3969/j.issn.1001-3881.2013.06.003

    猜你喜歡
    氣穴哈爾濱工業(yè)大學(xué)工程學(xué)院
    純水介質(zhì)雙向液壓鎖研究
    福建工程學(xué)院
    閥芯旋轉(zhuǎn)式高速開關(guān)閥的氣穴特性研究
    福建工程學(xué)院
    《哈爾濱工業(yè)大學(xué)學(xué)報》征稿簡則
    溫度效應(yīng)對液壓錐型節(jié)流閥內(nèi)氣穴形態(tài)的影響
    福建工程學(xué)院
    福建工程學(xué)院
    《哈爾濱工業(yè)大學(xué)學(xué)報》征稿簡則
    《哈爾濱工業(yè)大學(xué)學(xué)報》征稿簡則
    色94色欧美一区二区| 2021少妇久久久久久久久久久| 2022亚洲国产成人精品| 免费大片黄手机在线观看| 纯流量卡能插随身wifi吗| 亚洲国产最新在线播放| 免费av中文字幕在线| 波多野结衣一区麻豆| 中文字幕另类日韩欧美亚洲嫩草| 两个人免费观看高清视频| 999久久久国产精品视频| 亚洲精品国产av成人精品| 成人国语在线视频| av在线观看视频网站免费| 看十八女毛片水多多多| 热99久久久久精品小说推荐| 丝瓜视频免费看黄片| 九九爱精品视频在线观看| 亚洲四区av| 欧美日韩视频高清一区二区三区二| 精品国产一区二区久久| 国产黄频视频在线观看| 欧美激情高清一区二区三区 | 一边摸一边做爽爽视频免费| 热re99久久国产66热| 一级片'在线观看视频| 亚洲欧美成人综合另类久久久| 国产免费福利视频在线观看| 日本黄色日本黄色录像| 两个人免费观看高清视频| 亚洲美女搞黄在线观看| 一级黄片播放器| av电影中文网址| 久久97久久精品| 国产高清国产精品国产三级| 精品国产一区二区久久| 七月丁香在线播放| 精品一区二区三区四区五区乱码 | 国产免费福利视频在线观看| 老汉色av国产亚洲站长工具| 精品人妻在线不人妻| 一本大道久久a久久精品| 美女福利国产在线| 激情五月婷婷亚洲| 一边亲一边摸免费视频| 国产一区二区在线观看av| 91aial.com中文字幕在线观看| 亚洲一码二码三码区别大吗| 老汉色av国产亚洲站长工具| 熟女av电影| 人妻 亚洲 视频| 在线观看一区二区三区激情| 国产欧美日韩一区二区三区在线| 久久久精品区二区三区| 国产精品99久久99久久久不卡 | 国产精品人妻久久久影院| 成人黄色视频免费在线看| 天天操日日干夜夜撸| 黄色配什么色好看| 国产探花极品一区二区| 热99久久久久精品小说推荐| 国产日韩欧美在线精品| 精品人妻一区二区三区麻豆| 久久久精品区二区三区| 女性被躁到高潮视频| 一级毛片我不卡| 国产免费视频播放在线视频| 国产综合精华液| 男女国产视频网站| 亚洲av.av天堂| 亚洲天堂av无毛| 亚洲人成电影观看| 国产成人精品久久二区二区91 | 一区二区三区激情视频| 中文乱码字字幕精品一区二区三区| 日韩欧美精品免费久久| 精品卡一卡二卡四卡免费| 中文字幕av电影在线播放| 亚洲内射少妇av| 国产成人免费观看mmmm| 亚洲综合色惰| 欧美日韩一级在线毛片| 婷婷色麻豆天堂久久| 99精国产麻豆久久婷婷| 久久久久精品人妻al黑| kizo精华| 亚洲精华国产精华液的使用体验| 波多野结衣av一区二区av| 18禁动态无遮挡网站| 国产精品国产三级国产专区5o| 在线免费观看不下载黄p国产| 曰老女人黄片| 日韩免费高清中文字幕av| 国产免费视频播放在线视频| 久久综合国产亚洲精品| 多毛熟女@视频| 精品久久蜜臀av无| 色吧在线观看| 最近最新中文字幕免费大全7| 亚洲精品,欧美精品| 国产黄色免费在线视频| 日韩一本色道免费dvd| 久久免费观看电影| 最近最新中文字幕大全免费视频 | 久久久久久久精品精品| 国产成人精品婷婷| 久久久久精品久久久久真实原创| 91国产中文字幕| 麻豆av在线久日| 纵有疾风起免费观看全集完整版| 日本-黄色视频高清免费观看| 91精品三级在线观看| h视频一区二区三区| 欧美老熟妇乱子伦牲交| 老女人水多毛片| 午夜福利视频在线观看免费| 亚洲第一av免费看| 国产淫语在线视频| 久久这里有精品视频免费| 国产精品二区激情视频| 人人妻人人澡人人爽人人夜夜| 日本-黄色视频高清免费观看| 香蕉丝袜av| 国产av码专区亚洲av| 黄片无遮挡物在线观看| 观看av在线不卡| 欧美日韩一级在线毛片| 亚洲美女视频黄频| 亚洲国产av新网站| 国产97色在线日韩免费| 大码成人一级视频| 国产精品人妻久久久影院| 国产成人精品婷婷| 男人添女人高潮全过程视频| 妹子高潮喷水视频| 丝袜在线中文字幕| 999久久久国产精品视频| 亚洲精品一区蜜桃| 国产精品久久久av美女十八| 久久狼人影院| 天天影视国产精品| 在线观看免费视频网站a站| 一区二区日韩欧美中文字幕| 午夜影院在线不卡| 国产成人午夜福利电影在线观看| 国产一区有黄有色的免费视频| 伊人久久国产一区二区| 91成人精品电影| 免费少妇av软件| 久久久精品免费免费高清| 久久久久精品人妻al黑| 国产一区二区激情短视频 | 丁香六月天网| 十八禁高潮呻吟视频| 少妇被粗大的猛进出69影院| 欧美+日韩+精品| 美女高潮到喷水免费观看| 女人久久www免费人成看片| 欧美激情高清一区二区三区 | 考比视频在线观看| 国产毛片在线视频| 国产精品二区激情视频| 久久久久久久久久人人人人人人| 精品国产乱码久久久久久男人| 亚洲精品自拍成人| 国产熟女午夜一区二区三区| 精品亚洲成a人片在线观看| 男女边吃奶边做爰视频| 欧美中文综合在线视频| 在线天堂中文资源库| 十八禁网站网址无遮挡| 久久人妻熟女aⅴ| 国产女主播在线喷水免费视频网站| 丝瓜视频免费看黄片| 国产野战对白在线观看| 午夜影院在线不卡| 国产xxxxx性猛交| 高清视频免费观看一区二区| 国产不卡av网站在线观看| 色94色欧美一区二区| 国产精品无大码| 亚洲欧美成人综合另类久久久| a级片在线免费高清观看视频| 久久婷婷青草| 1024香蕉在线观看| 亚洲精品美女久久av网站| 欧美+日韩+精品| 欧美最新免费一区二区三区| 高清黄色对白视频在线免费看| 男女啪啪激烈高潮av片| 男女下面插进去视频免费观看| 欧美精品一区二区大全| 又大又黄又爽视频免费| 高清av免费在线| 少妇人妻 视频| www.精华液| 国产av国产精品国产| 在线 av 中文字幕| 丰满饥渴人妻一区二区三| 黄片小视频在线播放| 久久久久国产精品人妻一区二区| 热re99久久精品国产66热6| 男人舔女人的私密视频| 日日啪夜夜爽| 人妻 亚洲 视频| www.av在线官网国产| 久久精品国产亚洲av高清一级| 国产成人精品无人区| 国产成人精品无人区| 亚洲综合色惰| 国产精品三级大全| 久久99热这里只频精品6学生| av不卡在线播放| 精品少妇一区二区三区视频日本电影 | 亚洲国产精品一区三区| 精品久久久精品久久久| 侵犯人妻中文字幕一二三四区| 久久久久久久大尺度免费视频| 亚洲精品国产av蜜桃| 久久狼人影院| a级片在线免费高清观看视频| 国产在线免费精品| av.在线天堂| 街头女战士在线观看网站| 日本vs欧美在线观看视频| 亚洲av电影在线进入| 高清视频免费观看一区二区| 最近中文字幕2019免费版| 欧美精品一区二区大全| 少妇被粗大猛烈的视频| 亚洲综合色网址| xxx大片免费视频| 成人18禁高潮啪啪吃奶动态图| 好男人视频免费观看在线| 一级a爱视频在线免费观看| 高清在线视频一区二区三区| 免费在线观看完整版高清| 免费高清在线观看日韩| 免费人妻精品一区二区三区视频| 成人国产av品久久久| 一边摸一边做爽爽视频免费| 午夜福利视频在线观看免费| 国产精品二区激情视频| 精品国产一区二区三区久久久樱花| 18在线观看网站| 久久精品久久久久久噜噜老黄| 多毛熟女@视频| 久久久久久伊人网av| 中文字幕色久视频| 精品国产一区二区久久| 久久人人爽av亚洲精品天堂| 亚洲精品国产av成人精品| a级毛片黄视频| 久久人人97超碰香蕉20202| 伦精品一区二区三区| 中文天堂在线官网| 最近最新中文字幕免费大全7| 观看美女的网站| 国产成人aa在线观看| 男女午夜视频在线观看| 亚洲第一青青草原| 永久网站在线| 99香蕉大伊视频| 欧美少妇被猛烈插入视频| 一级毛片我不卡| 自线自在国产av| 日韩熟女老妇一区二区性免费视频| 亚洲激情五月婷婷啪啪| 免费播放大片免费观看视频在线观看| 夫妻午夜视频| 99久久综合免费| 婷婷色av中文字幕| 久久狼人影院| 欧美日韩综合久久久久久| 亚洲色图 男人天堂 中文字幕| 中文字幕色久视频| 国精品久久久久久国模美| 一边摸一边做爽爽视频免费| av在线老鸭窝| 欧美激情极品国产一区二区三区| 26uuu在线亚洲综合色| 咕卡用的链子| 欧美97在线视频| 久久精品亚洲av国产电影网| 毛片一级片免费看久久久久| 亚洲图色成人| 欧美日韩一区二区视频在线观看视频在线| 亚洲av福利一区| 最近手机中文字幕大全| 国产日韩欧美亚洲二区| 精品国产一区二区三区久久久樱花| 亚洲国产欧美在线一区| 国产精品熟女久久久久浪| 久久狼人影院| 亚洲av中文av极速乱| 成人影院久久| 久久久久久免费高清国产稀缺| 国产成人一区二区在线| 汤姆久久久久久久影院中文字幕| 综合色丁香网| 性色av一级| 久久精品aⅴ一区二区三区四区 | 国产福利在线免费观看视频| 国产片内射在线| 成人二区视频| 亚洲国产日韩一区二区| 最黄视频免费看| 欧美日韩国产mv在线观看视频| 亚洲av在线观看美女高潮| 国产成人精品福利久久| 久久久久视频综合| 日韩免费高清中文字幕av| 波多野结衣av一区二区av| 国产精品一区二区在线观看99| 波多野结衣一区麻豆| 亚洲欧洲精品一区二区精品久久久 | 国产又爽黄色视频| 黑丝袜美女国产一区| 99热国产这里只有精品6| 久久ye,这里只有精品| 五月天丁香电影| 91成人精品电影| 99国产综合亚洲精品| av福利片在线| 久久青草综合色| 久久久久久免费高清国产稀缺| 久久这里有精品视频免费| 最近最新中文字幕大全免费视频 | 国产乱来视频区| 秋霞伦理黄片| 亚洲精品aⅴ在线观看| 中文字幕色久视频| 丝袜美足系列| 在线亚洲精品国产二区图片欧美| 尾随美女入室| 大话2 男鬼变身卡| av免费在线看不卡| 久久精品国产综合久久久| 成年av动漫网址| 免费观看无遮挡的男女| 韩国av在线不卡| 久热这里只有精品99| 人成视频在线观看免费观看| 亚洲欧洲日产国产| a级毛片黄视频| tube8黄色片| 99国产综合亚洲精品| 99热网站在线观看| 精品少妇内射三级| 纯流量卡能插随身wifi吗| 午夜福利,免费看| 2022亚洲国产成人精品| 一级毛片电影观看| 少妇人妻精品综合一区二区| 日韩一区二区视频免费看| 国产成人精品久久久久久| 少妇的丰满在线观看| 久久精品国产综合久久久| 久久精品国产a三级三级三级| 国产一区二区激情短视频 | 丝袜人妻中文字幕| 多毛熟女@视频| 一二三四在线观看免费中文在| 精品一区二区三卡| 国产成人精品久久二区二区91 | 人妻一区二区av| 久久99精品国语久久久| 免费av中文字幕在线| 五月天丁香电影| 一区福利在线观看| 亚洲精品aⅴ在线观看| 国产成人精品婷婷| 国产片特级美女逼逼视频| 美女午夜性视频免费| 你懂的网址亚洲精品在线观看| 成人国产av品久久久| 国产亚洲精品第一综合不卡| 午夜福利在线免费观看网站| 一区二区av电影网| 王馨瑶露胸无遮挡在线观看| 精品人妻一区二区三区麻豆| 欧美激情高清一区二区三区 | 日本午夜av视频| 亚洲情色 制服丝袜| 亚洲精品在线美女| 高清视频免费观看一区二区| 两性夫妻黄色片| 国产日韩一区二区三区精品不卡| 精品国产超薄肉色丝袜足j| 日本黄色日本黄色录像| 99香蕉大伊视频| 毛片一级片免费看久久久久| 丝瓜视频免费看黄片| 美女脱内裤让男人舔精品视频| 亚洲成国产人片在线观看| 免费日韩欧美在线观看| 午夜免费观看性视频| 人成视频在线观看免费观看| 伦理电影免费视频| 另类精品久久| 国产白丝娇喘喷水9色精品| 日韩大片免费观看网站| 老司机亚洲免费影院| 午夜福利一区二区在线看| 久久99蜜桃精品久久| 久久99热这里只频精品6学生| 熟妇人妻不卡中文字幕| 欧美日韩一区二区视频在线观看视频在线| 国产一区亚洲一区在线观看| 九色亚洲精品在线播放| 亚洲人成网站在线观看播放| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| 午夜福利视频在线观看免费| 美女大奶头黄色视频| 观看美女的网站| 亚洲少妇的诱惑av| 亚洲av欧美aⅴ国产| 亚洲精品av麻豆狂野| videos熟女内射| 人成视频在线观看免费观看| 我的亚洲天堂| 久久精品aⅴ一区二区三区四区 | 久久久精品94久久精品| 18禁动态无遮挡网站| 国产极品天堂在线| 国产免费一区二区三区四区乱码| 国产精品香港三级国产av潘金莲 | 水蜜桃什么品种好| 久久鲁丝午夜福利片| 中文字幕色久视频| 人妻系列 视频| 精品人妻熟女毛片av久久网站| 大码成人一级视频| 亚洲av综合色区一区| 日本-黄色视频高清免费观看| 久久人人爽av亚洲精品天堂| 久久久a久久爽久久v久久| 婷婷色麻豆天堂久久| 咕卡用的链子| 日韩在线高清观看一区二区三区| 亚洲美女视频黄频| 欧美日韩av久久| 欧美日本中文国产一区发布| 亚洲国产欧美在线一区| 999精品在线视频| 国产成人精品久久二区二区91 | 91午夜精品亚洲一区二区三区| 一二三四中文在线观看免费高清| tube8黄色片| 日本色播在线视频| 久久国内精品自在自线图片| 精品一品国产午夜福利视频| 中文欧美无线码| 午夜老司机福利剧场| 亚洲成av片中文字幕在线观看 | 亚洲 欧美一区二区三区| 国产日韩欧美在线精品| 日本av免费视频播放| 日本欧美国产在线视频| 国产探花极品一区二区| av又黄又爽大尺度在线免费看| 亚洲精品中文字幕在线视频| 肉色欧美久久久久久久蜜桃| 久久人人爽av亚洲精品天堂| 麻豆乱淫一区二区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲久久久国产精品| 街头女战士在线观看网站| 亚洲欧美一区二区三区久久| 免费看不卡的av| 十分钟在线观看高清视频www| 97精品久久久久久久久久精品| 成人免费观看视频高清| 男女高潮啪啪啪动态图| 亚洲精品一区蜜桃| 亚洲精品日韩在线中文字幕| 中文字幕av电影在线播放| 青春草视频在线免费观看| 欧美日韩精品成人综合77777| 欧美 亚洲 国产 日韩一| 亚洲国产精品999| 日韩大片免费观看网站| 婷婷色av中文字幕| 成人毛片60女人毛片免费| 91精品伊人久久大香线蕉| 成人午夜精彩视频在线观看| 日韩中字成人| 一级黄片播放器| 天堂中文最新版在线下载| av卡一久久| 亚洲在久久综合| 欧美国产精品一级二级三级| 九九爱精品视频在线观看| 国产精品偷伦视频观看了| 性色avwww在线观看| 十八禁网站网址无遮挡| www.精华液| 国产乱人偷精品视频| 亚洲色图 男人天堂 中文字幕| 亚洲一区中文字幕在线| 成人免费观看视频高清| 伦理电影免费视频| 国产精品麻豆人妻色哟哟久久| 成人国产av品久久久| 少妇被粗大猛烈的视频| 久久久久久免费高清国产稀缺| av网站免费在线观看视频| 夜夜骑夜夜射夜夜干| 亚洲国产日韩一区二区| 久久午夜综合久久蜜桃| 97在线人人人人妻| 日韩伦理黄色片| 麻豆乱淫一区二区| 亚洲一码二码三码区别大吗| videossex国产| 午夜福利,免费看| 男女边吃奶边做爰视频| 精品人妻一区二区三区麻豆| 制服诱惑二区| 日韩一区二区三区影片| 1024视频免费在线观看| 午夜福利影视在线免费观看| 欧美激情极品国产一区二区三区| 丰满迷人的少妇在线观看| 日韩欧美一区视频在线观看| 亚洲av男天堂| 晚上一个人看的免费电影| 1024视频免费在线观看| 午夜老司机福利剧场| 在线天堂中文资源库| 欧美中文综合在线视频| 国产色婷婷99| 日本wwww免费看| 亚洲av免费高清在线观看| 美女高潮到喷水免费观看| www.精华液| 日韩欧美一区视频在线观看| 一本色道久久久久久精品综合| 少妇精品久久久久久久| 久久久精品国产亚洲av高清涩受| 成人二区视频| 91午夜精品亚洲一区二区三区| 春色校园在线视频观看| 国产av精品麻豆| 亚洲 欧美一区二区三区| 国产xxxxx性猛交| 国产一区二区在线观看av| 999精品在线视频| 超色免费av| 亚洲欧美精品综合一区二区三区 | 亚洲精品视频女| 欧美日韩国产mv在线观看视频| 欧美激情高清一区二区三区 | 99久久综合免费| 国产欧美亚洲国产| 看免费av毛片| 深夜精品福利| 91久久精品国产一区二区三区| 一区二区日韩欧美中文字幕| 各种免费的搞黄视频| 久久久a久久爽久久v久久| 国产97色在线日韩免费| 老女人水多毛片| 亚洲综合色惰| 国产在线一区二区三区精| 一区在线观看完整版| 一区二区三区激情视频| 久久99蜜桃精品久久| 久热久热在线精品观看| 国产精品成人在线| 国产乱来视频区| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 午夜老司机福利剧场| 免费高清在线观看日韩| 亚洲精品av麻豆狂野| 老汉色av国产亚洲站长工具| 亚洲综合色网址| 91成人精品电影| 日韩不卡一区二区三区视频在线| 亚洲国产精品成人久久小说| 成人亚洲精品一区在线观看| 天堂8中文在线网| 蜜桃在线观看..| 国产精品 欧美亚洲| 美女中出高潮动态图| 中文字幕色久视频| 涩涩av久久男人的天堂| 女人精品久久久久毛片| 精品人妻一区二区三区麻豆| 最黄视频免费看| 韩国av在线不卡| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产av新网站| 精品亚洲成a人片在线观看| 日韩三级伦理在线观看| 亚洲精品久久成人aⅴ小说| 一级毛片电影观看| 国产精品国产av在线观看| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| 久久av网站| 久久精品久久精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 一二三四中文在线观看免费高清| 女性生殖器流出的白浆| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 97精品久久久久久久久久精品| 91国产中文字幕| 丝袜在线中文字幕|