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

    Effect of Radial Resistance Gap on the Pressure Drop of a Compact Annular-Radial-Orifice Flow Magnetorheological Valve

    2019-01-17 01:11:14GuoliangHuJiaweiZhangMingkeLiaoandRuqiDing

    Guoliang Hu, Jiawei Zhang, Mingke Liao and Ruqi Ding

    (Key Laboratory of Conveyance and Equipment, Ministry of Education, East China Jiaotong University, Nanchang 330013, China)

    Abstract: A compact annular-radial-orifice flow magnetorheological (MR) valve was developed to investigate the effects of radial resistance gap on pressure drop. The fluid flow paths of this proposed MR valve consist of a single annular flow channel, a single radial flow channel and an orifice flow channel through structure design. The finite element modelling and simulation analysis of the MR valve was carried out using ANSYS/Emag software to investigate the changes of the magnetic flux density and yield stress along the fluid flow paths under the four different radial resistance gaps. Moreover, the experimental tests were also conducted to evaluate the pressure drop, showing that the proposed MR valve has significantly improved its pressure drop at 0.5 mm width of the radial resistance gap when the annular resistance gap is fixed at 1 mm.

    Key words: magnetorheological (MR) valve; annular-radial-orifice flow; pressure drop; radial resistance gap; finite element analysis

    As a kind of smart material, magnetorheological (MR) fluid exhibits unusual characteristics in that its rheological properties can be continuously and reversibly changed within milliseconds by applying or removing a magnetic field[1]. This interesting property has inspired the design of a large variety of MR fluid-based devices in various engineering applications, such as MR shock absorbers and dampers, MR engine mounts, MR brakes and clutch systems. Since the working principles of most MR fluid-based devices are based on the manipulation of fluid flow rate in the resistance gaps, the key performance of MR fluid-based devices is determined by the performance of MR valve[2-4]. Therefore, a successful improvement of MR valve performance can significantly impact the development of other MR fluid-based devices.

    In general, a typical MR valve consists of a valve body, a valve spool, an excitation coil and a fluid flow path. The detailed configuration of each MR valve could be different depending on the specific design, but the basic working principle is usually similar. The magnetic field of the excitation coil is applied to the MR fluid in the fluid flow path to make its rheological properties change. This leads to changes in the fluid flow resistance and therefore changes in the pressure drop of the valve that allows the MR fluid flow to be slowed or even stopped.

    In recent years the research of MR valves aims to achieve large pressure drops or fast responses by either designing a novel valve structure or optimizing an existing valve structure. Among them, there are three typical fluid flow paths to configure the geometrical arrangements of the MR valves: the annular flow path, the radial flow path, and both the annular and radial flow paths.

    The annular type MR valve is the most common type of valve. It has already been altered into several design varieties such as the inner coil design[5], the outer coil design[6-7]and the multiple-coil design[8]. Normally, the annular gap is characterized by the MR fluid flow channel longitudinal to the axial fluid flow, where it is usually parallel to the inlet and outlet of the valve. So, elongating the annular gap is the only way to increase the effective region as it will increase the pressure drop. Nevertheless, the elongation of the flow channel will enlarge the overall dimension of the valve as well, which might not be preferable in some valve designs.

    Radial type MR valves have also been introduced as an alternative to the annular type MR valves with the capability of increasing the pressure drop capacity by providing modular stages[9-11]. It is characterized by the way MR fluid flow path is transversal in the valve; that is, the radial gap is placed perpendicularly to the inlet and outlet of the valve. The radial type MR valve can help generate a large pressure drop with a larger width and a shorter longitudinal flow channel length. For example, Aydar et al. developed a large-scale modular MRF bypass damper with a two-stage disk type bypass MR valve which can provide a pressure drop of over 9.6 MPa at 5A of activation electric input current[9]. Liao et al. developed a novel MR valve with multi-stage radial fluid flow paths, and an MR damper was also prototyped and integrated with this bypass MR valve. The simulation and experimental results showed that this new MR damper is superior to the traditional MR damper with annular damping channel[10].

    The latest advance for MR valves is a combination of annular and radial type valves, which has shown promising performance by effectively multiplying the number of effective gaps. For example, Wang and Ai proposed an MR valve simultaneously possessing annular and radial fluid flow resistance channels with the assumption that the magnetic flux densities at the annular and radial fluid flow gaps are identical[12-13]. The results showed that the radial fluid flow gaps in the MR valve were more efficient in providing a larger, controllable range than those with annular fluid flow gaps. Imaduddin and Mazlan further developed a compact MR valve with multiple annular and radial gaps[14-16]. The simulation and experimental results revealed that this novel design could increase the achievable pressure drop of the valve by more than 2.5 MPa.

    As mentioned above, the performance of the valve with the combination of radial and annular fluid flow paths was the best, followed by the valve with the single radial fluid flow path, and then the worst was the valve with the single annular fluid flow path. This was based on the assumption that the valves were constrained in a specific volumetric dimension and the annular and radial resistance gaps or both of them were fixed. The thickness of the fluid flow resistance gap is an important factor that affects the performance of the valves, especially the pressure drop. Because the magnetic resistance in the fluid flow gap is one thousand times that of the magnetic materials used in the valve spool and body, a large gap thickness will change the magnetic circuit of the MR valve and decrease the magnetic flux density, which in turn will decrease the MR effect. If the gap thickness is too small, it will be difficult for the MR fluid to pass through the gaps even without the magnetic field being applied to the gaps due to the oil film effect. In a word, the feasible gap thickness should be determined when an MR valve is designed. In general, the reasonable gap thickness ranges between 0.5 mm and 1.5 mm[12].

    In this paper, a compact MR valve with annular-radial-orifice flow paths was developed. First, the fluid flow paths were composed of a single annular flow channel, a single radial flow channel and an orifice flow channel through compact structural design. This valve structure not only guarantees a larger pressure drop but also reduces the geometry of the proposed MR valve to some extent. Then the effect of the radial resistance gap of the proposed MR valve on the pressure drop was investigated through simulation analysis and experimental verification. Finally, the effects of the radial resistance gaps on the pressure drop using experimental methodology is being discussed.

    1 Structure Design and Model of the Annular-Radial-Orifice Flow MR Valve

    1.1 Structural analysis of the proposed MR valve

    The structure of the annular-radial-orifice flow MR valve is shown in Fig.1. In the MR valve, the right and left end covers are made from stainless steel materials, and they have six threaded holes as a holder for the screw and one fluid flow channel port acting as the inlet/outlet port of the MR valve. The flow guiding plate, which is made from stainless steel material, acts as a connecting passage between the inlet port and the annular flow channel, and it can allow the fluid to be in a state of laminated flow so as to avoid turbulence. The excitation coil was wound onto the groove formed by the fixed plate and the valve spool and protruded out through a hole in the right end cover. The assembly between the fixed plate and the valve spool was dependent on the screw joints at the left end of the valve spool. The screw was used to fix the magnetic plate, the washer and the fixed plate together. The radial resistance gap was guaranteed by the thickness of the washer. The radial resistance gap can be changed by replacing the different sizes of the washers. Here, four different sizes of washers were manufactured to investigate the effects of radial resistance gaps on the pressure drop in the following experimental tests.

    1—left end cover; 2—flowguiding plate; 3—screw; 4.magnetic plate; 5—washer; 6—fixed plate; 7—valve body; 8—excitation coil; 9—valve spool; 10—right end cover Fig.1 Schematic diagram of the annular-radial-orifice flow MR valve

    The MR fluid flow paths in the proposed MR valve were comprised of a single annular flow path, a radial flow path and an orifice flow path. In detail, the annular flow path was formed between the outer circumference of the magnetic plate and the inner circumference of the valve body, and the annular resistance gap was fixed at 1.0 mm[12]. The radial fluid flow path was looped between the right surface of the magnetic plate and the left surface of the valve spool. There was a through hole with a diameter of 4 mm in the middle of the valve spool, and it formed the orifice flow path of the MR valve. In addition, the proposed MR valve has an outer diameter size of 62 mm and an overall length of 80 mm.

    When a direct current is applied to the excitation coil, the closed loop magnetic circuits will be generated in the valve body, the annular flow resistance gap, the magnetic plate, the radial resistance gap and the valve spool. The intense magnetic fields will generate in the annular and radial resistance gaps that are perpendicular to the gaps. Finally, a pressure drop between the inlet port and the outlet port of the MR valve will be produced because of the shear stress in both of the resistance gaps. Therefore by adjusting the direct current, this pressure drop can be controlled.

    1.2 Magnetic circuit of the proposed MR valve

    Fig.2 shows the simplified magnetic circuit of the proposed MR valve. According to the principle of the continuity of magnetic flux that is used to determine the magnetic flux density throughout its conduit, the magnetic flux can be considered as

    ΦMR,r=ΦMR,a=Φsteel=Φ

    (1)

    whereΦis the magnetic flux of the circuit,ΦMR,aandΦMR,rare the magnetic flux of the MR fluid in the annular fluid flow path and radial fluid flow path, respectively, andΦsteelis the magnetic flux of the valve body, valve spool and magnetic plate in the primary magnetic path.

    Fig.2 Simplified magnetic circuit of the proposed MR valve

    The magnetic circuit is analyzed using the magnetic Kirchhoff’s law as

    (2)

    whereNcis the number of turns of the excitation coil, andIis the current applied to the excitation coil;Hiis the magnetic field intensity in theith link of the circuit, andliis the overall effective length of that link.

    On the other hand, the magnetic flux conservation rule of the circuit is given by

    Φ=∮cBdS=BiSi

    (3)

    whereBiandSiare the magnetic flux density and cross-sectional area of theith link, respectively.

    For the proposed MR valve that is shown in Fig.1 and Fig.2, the effective length of theith link of the circuit can be defined as

    (4)

    The cross-sectional area of theith link in the magnetic flux path can be calculated by

    (5)

    According to electromagnetic theory, the magnetic flux density,B, can be approximately expressed in terms of the magnetic field intensityHas

    Bi=μ0μiHi

    (6)

    whereμ0=4π×10-7T·mA-1is the magnetic permeability of free space, andμiis the relative magnetic permeability of magnetic materials.

    Rirepresents the magnetic resistance of theith link in the magnetic flux path, and the equation ofRicab be expressed as

    (7)

    So, the equation of magnetic Kirchhoff’s law can be rewritten as

    (8)

    The magnetic flux density of each part of the MR valve can be expressed in the following formand respectively upper bounded by the magnetization property of the used magnetic materials.

    (9)

    whereBjsatis the saturated magnetic flux density of the corresponding material in thejth link.

    According to Eqs.(7) (9), the magnetic flux density of the MR fluid in the annular fluid flow path and radial fluid flow path can be obtained, respectively. They are expressed as

    (10)

    (11)

    1.3 Control model of the MR valve

    In general, the performance of an MR valve can be determined by its capability to create and regulate pressure drop between the inlet port and the outlet port; although, the other capabilities of MR valve such as time constant, linearity and hysteresis are also important.

    Ignoring the pressure drop induced in the left and right end covers and the flow guiding plate, it can be seen that the flow paths consist of one annular gap, one radial gap, and one orifice gap together in the MR valve from Fig.2. So the pressure drop Δpof the MR valve can be defined as

    Δp=Δpa+Δpr+Δpo

    (12)

    where Δpa, Δpr, Δporepresent the pressure drop through the annular fluid flow path, radial fluid flow path and the orifice flow channel, respectively. The pressure drop Δpain the annular resistance gap is expressed as[12]

    (13)

    where Δpa,ηand Δpa,τare the pressure drop from the viscous properties of the fluid and the pressure drop from the field dependent yield stress of the fluid in the annular resistance gap, respectively.ηis the viscosity without the magnetic field,qis the system flow rate,gais the thickness of the annular resistance gap,Ris the radius of the valve, andthandLaare the thickness of the valve body and magnetic plate, respectively.τy,ais the change of yield stress in the annular resistance gap that responds to an applied magnetic field, andcis the coefficient which depends on the flow velocity profile and ranges from a minimum value of 1 to a maximum value of 3.

    The pressure drop Δprin the radial resistance gap is defined as[12]

    (14)

    where Δpr,ηand Δpr,τare the pressure drop from the viscous properties of the fluid and the pressure drop from the field dependent yield stress of the fluid in the radial resistance gap, respectively.τy,ris the change of yield stress in the radial resistance gap that responds to an applied magnetic field,gris the thickness of the radial resistance gap,R0is the radius of the central hole of the valve spool,Wcis the thickness of the groove of the valve spool,Rcis the radius of the valve spool at the left end without screw thread, andRdis the radius of the valve spool at the left end with screw thread.

    The pressure drop Δpoin the orifice resistance gap is defined as

    (15)

    whereLis the length of the valve spool.

    By substituting Eqs.(13)-(15) into Eq.(12), the pressure drop Δpof the proposed MR valve is obtained. In order to make it clear whether pressure drop is controllable or not, the range change is defined as the ratio of active pressure drop Δpτand passive pressure drop Δpη.

    (16)

    2 Modeling and Simulation of the Annular-Radial-Orifice Flow MR Valve

    The MR fluid with the type MRF-J01T provided by the Chongqing Instrument Material Research Institute in China was used in the following simulations and experiments[11].

    Fig.4 Magnetic flux density of the defined path

    Fig.3 shows the axisymmetric two-dimensional finite element model of the proposed MR valve using ANSYS/EMAG software. The magnetic plate and the valve spool and body were made from No.10 steel, which has a saturated magnetic flux density of 2.4 T when the applied magnetic field intensity exceeds 320 kA/m. The excitation coil is made of copper with a relative permeability of 1. The annular and radial resistance gaps are full of MR fluid whose permeability is defined by the B-H curve of MRF-J01T. In the simulation, the annular resistance gap was set to 1.0 mm, while the radial resistance gap was set to 1.5 mm, 1.0 mm, 0.8 mm and 0.5 mm, respectively. The current density in the excitation coil was set at 2.51 A/mm2when the applied current to the excitation coil equaled 2.0 A. The diameter of the excitation coil was 0.6 mm, and the number of the turns was 400. In the figure, the entity model was meshed using a quadrangular element with a total number of elements of 1 629 and a total number of nodes of 5 036.

    In order to investigate the distribution of the magnetic flux density in the fluid flow paths clearly, the pathP1along the annular fluid flow path and the pathP2along the radial fluid flow path were defined separately. The dimensions of the path position along the fluid flow paths are also listed, as shown in Fig.3.

    Fig.3 Two dimensional finite element model of the MR valve

    Fig.4 shows the magnetic flux density along the defined paths in the annular fluid flow path and in the radial fluid flow path, respectively. It can be seen that the magnetic flux density, both in the annular fluid flow path and in the radial fluid flow path, distributes evenly though it decreases sharply at the end of the defined path. In addition, the magnetic flux density at the radial fluid flow path is larger than that at the annular fluid flow path under the same radial resistance gap. Moreover, the less the radial resistance gap is, the bigger the magnetic flux density is, both in the radial resistance gap and in the annular resistance gap. The average magnetic flux density in the annular fluid flow path is 0.737 T, and the average magnetic flux density in the radial fluid flow path is 0.837 T, when the current applied to the excitation coil is 2.0 A, the radial resistance gap is 0.5 mm, and the annular resistance gap is 1.0 mm.

    Fig.5 shows the distribution of the magnetic flux density along the fluid flow paths under the different radial resistance gaps when the applied current is 2.0 A. It can be seen that the magnetic flux density increased both in the annular and radial resistance gaps when the radial resistance gap decreased from 1.5 mm to 0.5 mm, and the magnetic flux density in the radial fluid flow path was bigger than that in the annular fluid flow path. However, the deviation of the magnetic flux density in the radial and annular resistance gap is not big, which shows the high utilization of the magnetic flux density through the compact design for the proposed MR valve.

    Fig.5 Magnetic flux density along the fluid flow paths under the different radial resistance gaps

    Fig.7 Simulations of the magnetic flux density under different applied currents

    Fig.6 shows the changes of yield stress along the fluid flow paths under four different radial resistance gaps. Here,currents of 0.8 A and 2.0 A were applied to the excitation coil, respectively. Observing the figure, the trend of the yield stress variation is shown to be similar to the trend of the magnetic flux density variation as depicted in Fig.5 when the applied current is 0.8 A. However, the values of the yield stress remain stable when the applied current exceeds 2.0 A; the reason is that the change of yield stress is saturated under this applied current.

    Fig.6 Changes of yield stress along the fluid flow path under the different radial resistance gaps

    Fig.7a shows the relationships between the magnetic flux densities at the radial resistance gap under different applied currents with four different radial resistance gaps, and Fig.7b shows the relationships between the magnetic flux densities at the annular resistance gap. Both figures show that the magnetic flux density increased as the applied current increased when the resistance gap was fixed but decreased as the radial resistance gap increased. This occurred because the magnetic resistance increased as the resistance gap became wider, which in turn led to a decrease in the magnetic flux density.

    Fig.8 is the estimation of the pressure drop under four types of radial resistance gaps. The pressure drop increased with the increase of the applied current, and the pressure drop of the MR valve increased with the decreasing of the radial resistance gaps. Moreover, the less applied current was needed for the saturation of the pressure drop when the radial resistance gap changed from 1.5 mm to 0.5 mm.

    Fig.8 Effect of radial resistance gap size on the pressure drop

    3 Performance Analysis of the Annular-Radial-Orifice Flow MR Valve

    3.1 Prototype of the MR valve

    In order to make full use of magnetic flux density in the resistance gap, the magnetic flux density in the radial flow channel was designed to be roughly equal to that in the annular flow channels, so the reasonable resistance length ofLrandLacan also be obtained.

    From Eq.(3), it can be deduced as

    BMR,rSMR,r=BMR,aSMR,a

    (17)

    As mentioned above,BMR,ris roughly equal toBMR,a, soSMR,ris also roughly equal toSMR,a. According to Eq.(5), the relationship ofLrandLacan be determined. Here,Lris 11 mm andLais 4.2 mm.

    Fig.9 shows the prototype of the proposed MR valve with annular-radial-orifice fluid flow paths.

    3.2 Experimental test rig setup

    In order to validate the valve performance of the proposed MR valve, an experimental test rig was built up and shown in Fig.10. A motor-driven fixed gear pump was used as a power unit. Pressure transducers (a) and (b) were used to measure the inlet pressure and the outlet pressure of the MR valve, respectively. A relief valve (a) was used as a safety valve to protect the hydraulic system, and a relief valve (b) was adopted to simulate the load cases in the hydraulic system. A DC power (a) was used to supply power to the two pressure transducers, and the DC power (b) was applied to supply power to the excitation coil of the proposed MR valve. A data acquisition board was used to capture the pressures. A host computer was used to monitor the relevant test parameters of the hydraulic system in real time.

    Fig.9 Photograph of the manufactured MR valve

    Fig.10 Experimental test rig of the proposed MR valve

    3.3 Effect of different radial resistance gap sizes on the pressure drop

    In order to investigate the effects of radial resistance gaps on the pressure drop of the proposed MR valve, four different widths of washer were manufacturedas shown in Fig.11. In other words, four different radial resistance gaps from 0.5 mm to 1.5 mm can be obtained by changing the washers shown in Fig.11.

    Fig.11 Four sizes of washers used in the experiment

    Fig.12 shows the experimental pressure of the MR valve under applied currents with four different radial resistance gaps.The inlet port pressure and the pressure drop increased as the applied currents increased under the four different radial resistance gaps, while the outlet port pressure remained at a value of 500 kPa. The reason

    is that the pressure drop equals a subtraction of the inlet port pressure and outlet port pressure, while the outlet port pressure was held at a certain value by adjusting the outside valve knob of the relief valve (b) in Fig.10. However, the outlet port pressure fluctuated a little because the gears in the fixed gear pump were worn out due to long operating time and the abrasive nature of the MR fluid.

    Fig.12 Pressure change under applied currents with four different radial resistance gaps

    Fig.13 shows the experimental pressure drop of the MR valve under applied currents with four different radial resistance gaps. Itcan be seen that the pressure drop increased as the applied current increased. In addition, the pressure drop also increased as the radial resistance gap decreased from 1.5 mm to 0.5 mm. The maximum pressure drop is approximately 2 650 kPa at the applied current of 1.2 A when the radial resistance gap is 0.5 mm.

    Fig.13 Pressure drop under applied currents with four different radial resistance gaps

    Fig.14 shows the comparison of experimental pressure drop and simulated pressure drop under applied currents with radial gap sizes of 0.5 mm, 0.8 mm, 1.0 mm and 1.5 mm, respectively. As shown in the figure, the simulated pressure drop is lower than the experimental pressure drop at first. The possible reason may be that the pressure drop from the viscous properties of the fluid makes dominated contributions when the current is small. At the same time, the thickening effect is not considered in the magnetic field simulation based on the Bingham model, which causes the simulated viscosity to be less than the experimental viscosity. As the applied current increases, the pressure drop from the field dependent yield stress of the fluid plays a leading role; the simulated pressure drop is bigger than that of experimental pressure drop. The possible reason may be that the enhanced yield stress phenomenon occurs when several single chain structures of magnetic particles join together and form a thick column, which makes the saturated yield stress value in the simulation bigger than the experimental conditions.

    Fig.14 Pressure drop comparsion of simulation result and experimental result

    The range changes under four different radial resistance gaps can be obtained by Eq.(16) and the experimental results. It can be deduced that the range changesare 7.31, 5.21, 5.03 and 4.20 when the radial gap size was changed from 0.5 mm to 1.5 mm. The smaller the radial gap size is, the bigger the range changes are. Moreover, the range changes increase sharply when the radial gap size changes from 0.8 mm to the 0.5 mm, which shows a good pressure regulating ability.

    3.4 Effect of different load cases on the pressure drop of the MR valve

    The relief valve (b) was applied in the experimental test rig to simulate the loading conditions in the hydraulic system, and two typical load cases were selected by adjusting the outside valve knob in this study. The load case 1 was denoted by adjusting the outside valve knob clockwise once at the initial state, and the load case 2 was denoted by adjusting the outside valve knob clockwise twice at the initial state. Fig.15 shows the variation of pressure drop and load cases under different applied currents for the four radial resistance gaps. Here the pressure drop remained at a stable value when the current was increased from 0 A to 1.8 A under the two different load cases, which means the load cases did not influence the pressure drop. It is also seen that the pressure drop increased as the applied current increased, specifically.

    Fig.15 Variation of pressure drop under two typical load cases

    The pressure drop is not affected by load cases; this is an advantage because it makes the proposed MR valve useful as an active control element, For example, it can be used as a bypass valve for controlling the damper, which constitutes a bypass type MR valve controlled damper. It can also effectively expand the adjustable range of the damping force so that the MR damper can be used in the vibration damping and seismic application under different damping conditions.

    4 Conclusions

    A compact annular-radial-orifice flow MR valve was developed and prototyped. The fluid flow paths of the proposed MR valve consisted of a single annular flow channel, a single radial flow channel and an orifice flow channel through structural design. In order to investigate the change of the magnetic flux density and yield stress along the fluid flow paths under the four different radial resistance gaps, magnetic circuit design and finite element analysis were carried out. The simulation results show that this proposed MR valve can provide a larger pressure drop of 3 342 kPa at a 0.5 mm width of the radial resistance gap and a 1.0 mm width of annular resistance gap. Meanwhile, the volume space of MR valve has an outer diameter of 62 mm and overall length of 80 mm.

    The experimental pressure drop under the four different radial resistance gap was also carried out on the test rig. The results show that the pressure drop increases as the applied current increases, and the pressure drop also increases as the radial resistance gap decreases from 1.5 mm to 0.5 mm. The maximum pressure drop is approximately 2 650 kPa at the applied current of 1.2 A when the radial resistance gap is 0.5 mm, which is well accorded with the simulation results.

    This proposed MR valve has significantly improved its efficiency through compact design and changeable radial resistance gap. This is beneficial to control the MR damper as a bypass control valve for its good pressure regulating capability.

    91精品国产九色| 身体一侧抽搐| 国产精品秋霞免费鲁丝片| 噜噜噜噜噜久久久久久91| 99久久精品国产国产毛片| 别揉我奶头 嗯啊视频| 国产伦理片在线播放av一区| 婷婷色麻豆天堂久久| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 国产欧美日韩一区二区三区在线 | 国产成人免费观看mmmm| 91久久精品国产一区二区成人| 啦啦啦啦在线视频资源| 亚洲性久久影院| 亚洲精品乱码久久久久久按摩| 亚洲精品色激情综合| 18禁在线无遮挡免费观看视频| 韩国高清视频一区二区三区| 91狼人影院| 日韩一本色道免费dvd| 18禁裸乳无遮挡动漫免费视频 | 成人毛片a级毛片在线播放| 18禁裸乳无遮挡动漫免费视频 | 日韩av在线免费看完整版不卡| 国产成年人精品一区二区| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 黄色一级大片看看| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 国产欧美日韩一区二区三区在线 | 国内精品宾馆在线| 观看免费一级毛片| 国产免费福利视频在线观看| 中文精品一卡2卡3卡4更新| 国产亚洲精品久久久com| 老女人水多毛片| .国产精品久久| 亚洲精品,欧美精品| 成人国产av品久久久| 国产精品久久久久久精品古装| 午夜免费观看性视频| 欧美xxⅹ黑人| 国产成人福利小说| 亚洲精品第二区| 亚洲人成网站在线观看播放| 国产亚洲5aaaaa淫片| 一级毛片电影观看| 三级国产精品片| 人人妻人人看人人澡| 男人舔奶头视频| a级毛片免费高清观看在线播放| 丝袜脚勾引网站| 精品一区二区三卡| 亚洲精品久久久久久婷婷小说| 久久久久久久午夜电影| 国产一区亚洲一区在线观看| 亚洲欧美清纯卡通| 国产淫语在线视频| 成人毛片60女人毛片免费| 国产亚洲91精品色在线| 国产中年淑女户外野战色| 国产午夜精品一二区理论片| 亚洲精品,欧美精品| 欧美性感艳星| 欧美极品一区二区三区四区| 国产精品久久久久久久电影| 国产亚洲91精品色在线| 亚洲欧美一区二区三区国产| 高清视频免费观看一区二区| 3wmmmm亚洲av在线观看| av在线播放精品| 成人二区视频| 网址你懂的国产日韩在线| 国产精品秋霞免费鲁丝片| 一区二区三区免费毛片| 国产日韩欧美在线精品| 久久精品人妻少妇| 日本av手机在线免费观看| 亚洲精品中文字幕在线视频 | 人妻系列 视频| 中国美白少妇内射xxxbb| 久久99蜜桃精品久久| 在线观看美女被高潮喷水网站| 欧美3d第一页| 人人妻人人澡人人爽人人夜夜| 少妇人妻 视频| 好男人视频免费观看在线| 国产日韩欧美亚洲二区| 1000部很黄的大片| 亚洲av电影在线观看一区二区三区 | 白带黄色成豆腐渣| 国产欧美另类精品又又久久亚洲欧美| 水蜜桃什么品种好| 欧美+日韩+精品| av国产免费在线观看| 日日摸夜夜添夜夜爱| 亚洲不卡免费看| 国产av码专区亚洲av| 日韩三级伦理在线观看| 国产黄色视频一区二区在线观看| 精品少妇黑人巨大在线播放| 久久久久久久久大av| 春色校园在线视频观看| 人人妻人人看人人澡| 国产一区二区三区av在线| 久久精品久久精品一区二区三区| 直男gayav资源| 国产av码专区亚洲av| 午夜激情久久久久久久| 青春草国产在线视频| 伦精品一区二区三区| 少妇丰满av| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 97人妻精品一区二区三区麻豆| 欧美zozozo另类| 久久精品国产鲁丝片午夜精品| 亚洲第一区二区三区不卡| 国产国拍精品亚洲av在线观看| 51国产日韩欧美| 色5月婷婷丁香| 日本午夜av视频| 国产成人a∨麻豆精品| 极品少妇高潮喷水抽搐| 日韩一本色道免费dvd| 超碰av人人做人人爽久久| 特大巨黑吊av在线直播| 免费观看av网站的网址| 舔av片在线| 极品教师在线视频| 久久久久久久久久人人人人人人| 又粗又硬又长又爽又黄的视频| 久久精品夜色国产| 黄色一级大片看看| 2018国产大陆天天弄谢| 国产精品伦人一区二区| 老司机影院毛片| 国产亚洲91精品色在线| 欧美成人精品欧美一级黄| 免费av不卡在线播放| 爱豆传媒免费全集在线观看| 麻豆精品久久久久久蜜桃| 亚洲最大成人中文| 国产成年人精品一区二区| 免费黄色在线免费观看| 国产成人精品福利久久| 又粗又硬又长又爽又黄的视频| 三级国产精品片| 看十八女毛片水多多多| 男女边吃奶边做爰视频| 91精品国产九色| 99精国产麻豆久久婷婷| 中国国产av一级| 99热这里只有精品一区| 三级国产精品欧美在线观看| 天美传媒精品一区二区| 99热这里只有是精品50| 亚洲自拍偷在线| 国模一区二区三区四区视频| 熟女av电影| 99热国产这里只有精品6| 亚洲四区av| 熟女av电影| 国产精品av视频在线免费观看| 国产精品嫩草影院av在线观看| 可以在线观看毛片的网站| 97人妻精品一区二区三区麻豆| 18禁裸乳无遮挡免费网站照片| 日韩电影二区| 亚洲欧美成人综合另类久久久| 国产大屁股一区二区在线视频| 久久久久久久久久人人人人人人| 日韩强制内射视频| 伊人久久国产一区二区| 亚洲精品日韩在线中文字幕| 国产成人免费无遮挡视频| 日本一本二区三区精品| 国产精品国产三级专区第一集| 一个人观看的视频www高清免费观看| 国产免费又黄又爽又色| 国产永久视频网站| 国产一区亚洲一区在线观看| 搡老乐熟女国产| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久久久免费av| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 欧美成人一区二区免费高清观看| 婷婷色av中文字幕| 国产又色又爽无遮挡免| 久久久久久九九精品二区国产| 免费不卡的大黄色大毛片视频在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲精品乱码久久久v下载方式| 国产欧美亚洲国产| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 看黄色毛片网站| 97精品久久久久久久久久精品| 人人妻人人看人人澡| 久久精品久久久久久噜噜老黄| 一本一本综合久久| 成人高潮视频无遮挡免费网站| 国产成年人精品一区二区| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 欧美zozozo另类| 日韩视频在线欧美| 一区二区av电影网| 天堂俺去俺来也www色官网| 老司机影院毛片| 亚洲精品国产成人久久av| 成人亚洲精品av一区二区| 成人毛片60女人毛片免费| av国产精品久久久久影院| 在线观看一区二区三区激情| 在线观看一区二区三区| 亚洲国产精品成人久久小说| 成人免费观看视频高清| 久久久精品欧美日韩精品| 日日啪夜夜爽| 久久99蜜桃精品久久| 男插女下体视频免费在线播放| 熟女av电影| 一本色道久久久久久精品综合| 久久久久国产精品人妻一区二区| 午夜福利在线观看免费完整高清在| 久久综合国产亚洲精品| 美女被艹到高潮喷水动态| 日韩精品有码人妻一区| 国产精品久久久久久精品古装| 国内精品美女久久久久久| 国精品久久久久久国模美| 日韩一区二区视频免费看| 视频区图区小说| 少妇 在线观看| 欧美日韩综合久久久久久| 久久综合国产亚洲精品| 少妇的逼水好多| www.色视频.com| 久久精品久久久久久久性| 最近中文字幕2019免费版| 成人漫画全彩无遮挡| 又黄又爽又刺激的免费视频.| 成人黄色视频免费在线看| 在线免费十八禁| 亚洲成色77777| 综合色av麻豆| 国产一区有黄有色的免费视频| 深夜a级毛片| 丝袜脚勾引网站| 久久久久久久久久人人人人人人| 免费观看性生交大片5| 青春草视频在线免费观看| 久久精品国产自在天天线| 国产精品99久久99久久久不卡 | 国产成人aa在线观看| 波野结衣二区三区在线| av国产精品久久久久影院| 久久久精品免费免费高清| 亚洲综合色惰| 直男gayav资源| 欧美日韩在线观看h| 亚洲国产色片| 永久网站在线| 女人十人毛片免费观看3o分钟| 国产成人午夜福利电影在线观看| 中国三级夫妇交换| 直男gayav资源| 国产午夜精品久久久久久一区二区三区| 成年人午夜在线观看视频| 白带黄色成豆腐渣| 国精品久久久久久国模美| kizo精华| 亚洲精品国产色婷婷电影| 国产成人a区在线观看| 少妇被粗大猛烈的视频| 干丝袜人妻中文字幕| 国产高潮美女av| 好男人视频免费观看在线| 免费人成在线观看视频色| 只有这里有精品99| 国产综合精华液| 麻豆乱淫一区二区| 欧美激情久久久久久爽电影| 精品一区二区免费观看| 午夜福利视频1000在线观看| 成人午夜精彩视频在线观看| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 性色avwww在线观看| 高清视频免费观看一区二区| 欧美另类一区| 成人欧美大片| 色哟哟·www| 国产综合懂色| 亚洲av男天堂| 亚洲精品一区蜜桃| 在线观看一区二区三区| 久久久久久久国产电影| 婷婷色综合大香蕉| 欧美日韩亚洲高清精品| 97人妻精品一区二区三区麻豆| 久久人人爽人人片av| av在线蜜桃| 国语对白做爰xxxⅹ性视频网站| 大香蕉久久网| 久久韩国三级中文字幕| 五月开心婷婷网| 国产精品女同一区二区软件| 91狼人影院| 亚洲精品国产成人久久av| 国产淫语在线视频| 国产黄a三级三级三级人| 男女那种视频在线观看| 亚洲精品色激情综合| 又爽又黄a免费视频| 日韩伦理黄色片| 韩国高清视频一区二区三区| 国产精品成人在线| 国产精品一二三区在线看| 男的添女的下面高潮视频| 国产精品嫩草影院av在线观看| 插逼视频在线观看| 亚洲综合精品二区| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 成人高潮视频无遮挡免费网站| 国产在视频线精品| 99久久九九国产精品国产免费| 黄色怎么调成土黄色| 亚洲精品一二三| 夫妻午夜视频| 亚洲人成网站在线观看播放| 嫩草影院精品99| 男人添女人高潮全过程视频| 777米奇影视久久| 老师上课跳d突然被开到最大视频| 国产精品秋霞免费鲁丝片| 精品99又大又爽又粗少妇毛片| 国产午夜精品久久久久久一区二区三区| 国产中年淑女户外野战色| 天天躁夜夜躁狠狠久久av| 成年人午夜在线观看视频| 日本黄大片高清| 毛片女人毛片| av在线蜜桃| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 日韩亚洲欧美综合| 午夜福利高清视频| 97精品久久久久久久久久精品| 午夜福利高清视频| 99热全是精品| 国产精品无大码| 天堂俺去俺来也www色官网| 精品酒店卫生间| 欧美激情在线99| 成人高潮视频无遮挡免费网站| 黄色怎么调成土黄色| 欧美+日韩+精品| 综合色丁香网| 久久鲁丝午夜福利片| av在线播放精品| 岛国毛片在线播放| 精品久久久久久久末码| 熟女电影av网| 男人爽女人下面视频在线观看| av在线蜜桃| 在线免费十八禁| 亚洲av欧美aⅴ国产| 亚洲国产色片| freevideosex欧美| 国产精品久久久久久精品古装| 深夜a级毛片| 香蕉精品网在线| 日本午夜av视频| 中文字幕制服av| 99热国产这里只有精品6| 亚洲精品国产av蜜桃| eeuss影院久久| 欧美日韩亚洲高清精品| 久久99热这里只有精品18| av线在线观看网站| 免费电影在线观看免费观看| 欧美一级a爱片免费观看看| 在线观看三级黄色| 人妻 亚洲 视频| 激情五月婷婷亚洲| 特级一级黄色大片| 毛片一级片免费看久久久久| 插逼视频在线观看| 男人爽女人下面视频在线观看| 亚洲av欧美aⅴ国产| 狂野欧美白嫩少妇大欣赏| 午夜日本视频在线| 超碰97精品在线观看| 欧美日韩视频精品一区| 免费黄网站久久成人精品| 免费播放大片免费观看视频在线观看| 伊人久久国产一区二区| 亚洲av福利一区| 丝瓜视频免费看黄片| 日日撸夜夜添| 男女无遮挡免费网站观看| 国国产精品蜜臀av免费| 亚洲国产精品999| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 午夜福利视频精品| 你懂的网址亚洲精品在线观看| 男插女下体视频免费在线播放| 三级国产精品欧美在线观看| 秋霞在线观看毛片| 哪个播放器可以免费观看大片| 国产v大片淫在线免费观看| 少妇丰满av| av黄色大香蕉| 五月玫瑰六月丁香| 国产精品国产av在线观看| 欧美性猛交╳xxx乱大交人| 久久久久久久久久成人| 亚洲婷婷狠狠爱综合网| 我的老师免费观看完整版| 熟女人妻精品中文字幕| 国产精品偷伦视频观看了| 如何舔出高潮| 日韩大片免费观看网站| 精品视频人人做人人爽| 欧美一区二区亚洲| h日本视频在线播放| 一区二区三区免费毛片| 永久免费av网站大全| 国产中年淑女户外野战色| 亚洲欧美日韩无卡精品| 蜜桃久久精品国产亚洲av| 成人国产av品久久久| 啦啦啦在线观看免费高清www| 午夜免费观看性视频| 丰满乱子伦码专区| 国产伦在线观看视频一区| 国产 一区 欧美 日韩| 久久99热这里只有精品18| 一级爰片在线观看| 男女边吃奶边做爰视频| 赤兔流量卡办理| 观看美女的网站| 岛国毛片在线播放| 久久国内精品自在自线图片| 日韩一本色道免费dvd| 亚洲欧美中文字幕日韩二区| 亚洲真实伦在线观看| 18禁动态无遮挡网站| 日日撸夜夜添| 另类亚洲欧美激情| 国产成人免费无遮挡视频| 一级毛片黄色毛片免费观看视频| 97精品久久久久久久久久精品| 亚洲精品日韩av片在线观看| 97在线人人人人妻| 欧美另类一区| 国产精品一区二区三区四区免费观看| 午夜福利视频1000在线观看| tube8黄色片| 777米奇影视久久| 日韩在线高清观看一区二区三区| 在线观看免费高清a一片| 夜夜看夜夜爽夜夜摸| 亚洲欧美一区二区三区国产| 美女主播在线视频| 国产高清三级在线| 青春草亚洲视频在线观看| 大香蕉久久网| 亚洲精品日本国产第一区| 亚洲精品,欧美精品| 日本三级黄在线观看| 久久久国产一区二区| 国产毛片在线视频| 能在线免费看毛片的网站| 欧美精品人与动牲交sv欧美| 91精品国产九色| 最近最新中文字幕免费大全7| 人妻 亚洲 视频| 69人妻影院| 欧美日韩综合久久久久久| 亚洲成人久久爱视频| 91aial.com中文字幕在线观看| 别揉我奶头 嗯啊视频| 高清av免费在线| 在线播放无遮挡| 青春草亚洲视频在线观看| 尾随美女入室| 麻豆国产97在线/欧美| 成人美女网站在线观看视频| 免费播放大片免费观看视频在线观看| 如何舔出高潮| 成人黄色视频免费在线看| 干丝袜人妻中文字幕| av播播在线观看一区| 成年免费大片在线观看| 国产一区二区三区av在线| 国产爽快片一区二区三区| 少妇人妻久久综合中文| 在线观看一区二区三区激情| 国产一区二区亚洲精品在线观看| 日日摸夜夜添夜夜添av毛片| 精品少妇久久久久久888优播| 精品酒店卫生间| 亚洲综合精品二区| 交换朋友夫妻互换小说| 青春草国产在线视频| 日韩中字成人| 国产大屁股一区二区在线视频| 久久久国产一区二区| 欧美日本视频| 日日摸夜夜添夜夜爱| 高清午夜精品一区二区三区| 国产精品久久久久久精品电影小说 | 国产av码专区亚洲av| 久久久久精品性色| 免费观看在线日韩| 国模一区二区三区四区视频| 99久久九九国产精品国产免费| 性插视频无遮挡在线免费观看| 在线a可以看的网站| 97精品久久久久久久久久精品| 少妇人妻 视频| 久久久精品免费免费高清| 亚洲国产精品成人综合色| av网站免费在线观看视频| 国产欧美日韩一区二区三区在线 | 一个人看视频在线观看www免费| 男插女下体视频免费在线播放| 久久99蜜桃精品久久| 国产亚洲av片在线观看秒播厂| 亚洲av电影在线观看一区二区三区 | 亚洲av一区综合| 天堂中文最新版在线下载 | 亚洲精品自拍成人| 久久精品国产a三级三级三级| 成人特级av手机在线观看| 久久女婷五月综合色啪小说 | 1000部很黄的大片| 成人亚洲精品av一区二区| 亚洲内射少妇av| 欧美激情久久久久久爽电影| 亚洲成人久久爱视频| 九九久久精品国产亚洲av麻豆| 少妇的逼水好多| videossex国产| 免费看a级黄色片| 最近最新中文字幕大全电影3| 超碰av人人做人人爽久久| 国产午夜福利久久久久久| 美女高潮的动态| 亚洲av福利一区| 最后的刺客免费高清国语| 国产黄频视频在线观看| tube8黄色片| 嫩草影院新地址| 国产精品一区二区性色av| 极品教师在线视频| 国产亚洲午夜精品一区二区久久 | 黄色视频在线播放观看不卡| 99热国产这里只有精品6| 人妻 亚洲 视频| 欧美日韩国产mv在线观看视频 | av在线播放精品| 国内少妇人妻偷人精品xxx网站| 久久精品国产亚洲av天美| 综合色av麻豆| av.在线天堂| 男人狂女人下面高潮的视频| 国产成人福利小说| 老司机影院毛片| 成年免费大片在线观看| 国产成人精品婷婷| 国产白丝娇喘喷水9色精品| 最近的中文字幕免费完整| 日日撸夜夜添| 国产探花极品一区二区| 在线免费十八禁| www.色视频.com| 特级一级黄色大片| 国产人妻一区二区三区在| 成人国产麻豆网| 国产伦理片在线播放av一区| 好男人在线观看高清免费视频| 精品国产三级普通话版| 在线播放无遮挡| 色婷婷久久久亚洲欧美| 日韩不卡一区二区三区视频在线| 欧美老熟妇乱子伦牲交| 国产老妇伦熟女老妇高清| 最新中文字幕久久久久| 亚洲人成网站在线观看播放| 美女cb高潮喷水在线观看| 日韩中字成人| kizo精华| 亚洲精品日本国产第一区| 亚洲美女视频黄频| 尾随美女入室| 久久久久精品性色| 久久久久久伊人网av| 国产精品99久久久久久久久| av女优亚洲男人天堂| 久久久成人免费电影| 欧美人与善性xxx| 国产精品99久久99久久久不卡 | 搞女人的毛片| 国产 精品1| 久久久久久伊人网av| 看免费成人av毛片|