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

    Novel method for operating characteristics test of a hydraulic pump

    2015-04-24 05:30:24WANGXiangzhou王向周LUYe陸葉ZHENGShuhua鄭戍華
    關(guān)鍵詞:清陽(yáng)耳病肝氣

    WANG Xiang-zhou(王向周), LU Ye(陸葉), ZHENG Shu-hua(鄭戍華)

    (School of Automation, Beijing Institute of Technology, Beijing 100081, China)

    ?

    Novel method for operating characteristics test of a hydraulic pump

    WANG Xiang-zhou(王向周), LU Ye(陸葉), ZHENG Shu-hua(鄭戍華)

    (School of Automation, Beijing Institute of Technology, Beijing 100081, China)

    The traditional testing method of a hydraulic pump requires a large amount of test data from a variety of pump working conditions. The test is usually time-consuming and energy-consuming. And the accurate characteristic curves of the pump were hardly obtained due to a limited amount and discreteness of the test data. In order to simplify the test procedure and improve the test accuracy, a novel method for measuring hydraulic pump operating characteristic based on multi-element nonlinear regression (NLMR) modeling is proposed in this paper. The main idea of this modeling method is establishing a mathematical model to predict the performance parameters of the hydraulic pump, only a small amount of test data is needed. Consequently, the pump operating characteristics in any working conditions are obtained. And the test results of the pump are easily characterized in the graphs, charts, tables and so on. The evaluations of the model are carried out and discussed in this paper. The result shows that the test error of the novel method can be controlled to be about 0.1%. Compared with the traditional test method, the proposed method reduces greatly the test time and the random error of the test data, and improves the efficiency and accuracy of the pump test.

    hydraulic pump; test; model; work characteristics

    The traditional testing method demands that the hydraulic pump should be tested under a variety of pump working conditions. A lot of time and energy were consumed during the test procedure, and the test data were mixed with the random error. In recent years, more and more emerging technologies have been introduced to the operating characteristics test of pumps. For example, H Yang measured the flow of a centrifugal pump with 3-D particle image velocimetry (PIV), the maximum test error was about 3.14%[1]. ProMinent Dosiertechnik Gmbh adapted a computer to control and parameterize the individual components in a new pump test bench, all the input parameters were monitored, analyzed, and documented by an operating data collection system. The new pump test bench improved the test efficiency and the measured value has an accuracy of 0.1%[2]. Tao increased the measurement accuracy and automated the measurement operation in a nonlinear feedback with a proportional-integral (PI) controller[3]. Although these new technologies were helpful in the work characteristics test of pump to some extent, they can’t overcome the defects existing in the traditional testing method completely[4]. According to the traditional testing method, the test data was usually processed with simple numerical methods (such as an interpolation method) to present characteristic curves which could describe work characteristics of the pump. For example, Warren E Wilson proposed the method of calculating the performance curves with the test data, including the flow curve, the torque curve, the power curve, the efficiency curve, the equivalent efficiency curve and so on[5]. Shaharm Derekhshan developed the performance curves with an interpolation method and predicted the best efficiency point of a pump[6]. Kochevsky developed a method of computational fluid dynamics (CFD) to predict performance parameters of centrifugal pumps treating multiphase mixture[7]. The simple numerical methods couldn’t eliminate the random error mixed with the test data, sometimes the test error might be even amplified. So the accurate characteristic curves of the tested pump were hardly obtained with traditional testing methods. And due to a limited amount and discreteness of the test data, only a few characteristic curves could be acquired with the traditional testing method. Hence the work characteristics under various working conditions couldn’t be obtained by users. In addition, traditional testing method was difficult to meet the requirements of the performance test on hydraulic pumps. It required a large amount of experimental data under a variety of pump working conditions and complex experiments. In the test process, not only a lot of time and energy were consumed, but also test data was mixed with some measurement errors. Moreover, the representation of the test results was monotonous.

    The goal of this paper is to discuss a novel method for performance test of pump based onthe multi-element nonlinear regression (NLMR) theory. The test points can be reduced significantly while the efficiency and accuracy of the test are both improved. And the accurate work characteristics in any working condition can be obtained with just a small amount of test data, and the work characteristics can be characterized as the graphs, charts and tables and so on.

    1 Acquisition and analysis of the test data

    The modeling object should be determined before the mathematical model is established. In this paper, a work characteristics test on the hydraulic pump was conducted and a small amount of test data is acquired. By analyzing the test data, the appropriate modeling object can be determined.

    1.1 Test procedure and analysis of the test data

    The pump is tested according to the National Standard GB/T 17491—2011[8]. And the complete test system diagram is shown in Fig.1.

    Fig.1 Experiment system

    As shown in Fig.1, the pump is working in an open circuit and driven by an inverter motor. The relief valve is used to load for the pump. The torque and the rotational speed are measured by the torque gauge. The flows are measured by the volumetric flow sensor and the pressures are measured by the pressure sensors. The entire test data are collected in the computer through the DAQ device. The major components used to measure and collect data are shown in Tab.1.

    In the test, the pump of axial piston A7V80DR manufactured by REXROTH is Used. The inlet pressure of the pump is set at 0 MPa and the oil temperature is controlled in the range from 42.0 ℃ to 44.9 ℃. The test procedure is listed as follows:

    Keeping the pump rotational speed at 400 r/min and increasing the load of the pump gradually, the parameters are measured when the pressure is at 3.5, 7.0, 12.0, 16.0, 20.0, 25.0, 30.0 MPa, respectively.

    Increasing the pump rotational speed gradually, the above test procedure is repeated when the pump rotational speed is fixed at 400, 600, 825, 1 050, 1 275, 1 500, 1 600 r/min.

    Tab.1 List of the major components of the experiment system

    After measuring all the parameters which can be tested directly, other parameters such as the power and the efficiency can be calculated based on the test data. Finally 49 sets of test data are acquired in Tab.2.

    Tab.2 Experimental data from pump A7V80DR

    The rotational speed and the torque are pump’s inputs while the flows and pressures are the outputs of the pump. By analyzing the test datain in Tab.1, it can be found that the pump driving torque is determined by the differential pressure between the inlet and outlet of the pump and influenced by the pump rotational speed. Similarly, the output flow of the pump is determined by the rotational speed and influenced by the differential pressure.

    In theory, the input torque of the pump is proportional to the output pressure of the pump. But the actual values of input torque of the pump are greater than the ones in theory calculation due to the existence of coulomb friction between the bearings and other parts inside the pump. The difference between the theoretical torque and the actual torque is called the torque loss of the pump. The loss of torque is usually evaluated by mechanical efficiencyηt[9]

    (1)

    whereTtis the torque in theory;Tis the torque in actual.

    The output flow of the pump is proportional to the input rotational speed of the pump. But the actual flow is less than the one in theory due to the leakage of oil. The flow loss of the pump is usually evaluated by volumetric efficiencyηvas follows

    (2)

    whereQtis the flow in theory;Qis the flow in actual.

    In fact, the volumetric efficiency and the mechanical efficiency are the most important indicators in evaluating the performance of the pump[10]. Meanwhile, the experiments show that the volumetric efficiency and the mechanical efficiency are both evaluated by the loss characteristics of the pump. So the losing characteristics of the pump are selected for modeling the pump in this paper.

    1.2 Losing characteristics analysis

    The losing characteristics of the pump can fully reflect pump’s work characteristics under any working condition. And loss characteristics of the pump are more suitable for the modeling method than the actual torques and flows of the pump due to their smaller values. Therefore, the torque and flow loss characteristics of the pump are adopted to build the NLMR models.

    According to the 49 sets of test data, the torque loss and flow loss in test points are acquired as shown in Fig.2 and Fig.3. The scatter points represent the torque loss and the flow loss calculated by test points. The figures show that the torque loss and flow loss are the nonlinear and multi-parameter coupling systems. Hence a nonlinear function which can describe the losing characteristics should be determined first.

    Fig.2 The torque loss

    Fig.3 The flow loss

    2 NLMR models description

    2.1 NLMR theory

    The multi-element nonlinear regression (NLMR) is a simple and effective method in data processing. Specifically, it is a kind of method to find a nonlinear model describing the relationship between the dependent variable and a set of independent variables[11]. The mathematical model based on NLMR can extract a valuable common law from great amounts of engineering on-site data and the results can rapidly provide suggestions back to engineering practice[12]. Therefore, it has been more and more widely used in complex engineering problems which are nonlinear and multi-parameter coupling. Li et al. employed the model to reasonably eliminate the transmission congestion and optimize the electric power system[13]. Zhang et al. developed a nonlinear regression model which could successfully predict the toxicity of alkyphenol compounds[14]. However, the multi-element nonlinear regression NLMR theory is still not introduced into the performance work characteristics test on hydraulic pump.

    2.2 Establishment of NLMR models

    It can be known from section 2 that a nonlinear function is suitable for building the loss models. First, a class of functions are defined as

    φ=span{φ0,φ1,φ2,…,φn}

    (3)

    {φ0,φ1,φ2,…,φn} are the continuous and linear independent functions.

    Weierstrass approximation theorem states: iffis a continuous real-valued function on a closed interval and if anyε>0 is given, then there exists a polynomialPon the closed interval such that

    |f-P|<ε

    (4)

    In other words, every continuous function defined on a closed interval can be uniformly approximated by a polynomial function to any degree of accuracy. Consequently, polynomial functions can be selected from classφ, the structure is

    (5)

    whereaiis the undetermined coefficient.

    Considering that the mathematical model must be able to describe the relationship between the loss and the rotational speed or the differential pressure, a mathematical model which can describe the torque loss and flow loss is acquired after comparative researches and repeated experiments.

    The losing characteristics of pumps can be described in the NLMR models as follows:

    (6)

    where ΔPis the differential pressure;nis the rotational speed;mis the order of model;akis the undetermined coefficient;is the loss of torque or flow.

    2.3 Calculation of model coefficient

    The methods to calculate the model coefficients can be roughly categorized into two classes: the interpolation and the approximation. The interpolation method requires that the function curve must pass through all the original data points regardless of the random errors included in the test data. So the approximation method is preferred than the interpolation method. The least square method is a typical case of the approximation method.

    The third-order model is selected as an example to show how to calculate the coefficientakwith the least square method. The structure of the third-order model is as follows:

    Δ=a1+a2ΔP+a3n+a4ΔP2+a5ΔPn+a6n2+a7ΔP3+a8ΔP2n+a9ΔPn2+a10n3

    (7)

    Define

    (8)

    the third-order model can be linearized as follows:

    Δ=a1+a2x2+a3x3+a4x4+a5x5+a6x6+a7x7+a8x8+a9x9+a10x10

    (9)

    the sum of error squaresvis defined as:

    (10)

    where Δyiis the torque loss or the flow loss calculated by test data and Δiis the one calculated by the loss models.

    According to the least square theory,akis calculated based on the principle:

    v=[v]min

    (11)

    or

    (12)

    Therefore we have

    (13)

    Then the equations are simplified as follows:

    (14)

    The matrix representation is

    XTXA=XTΔy

    (15)

    Ais the coefficient matrix:

    A=[a1,a2,a3…a10]T

    (16)

    Δyis observation vector:

    與其他臟腑的關(guān)系:中醫(yī)又有“心開(kāi)竅于耳”之說(shuō),所以耳屬心腎二臟之竅,但以腎為主,以心為客。耳與肝膽相連,肝氣通于耳,肝氣調(diào)達(dá),則聽(tīng)力聰敏。肝膽主升發(fā),喜條達(dá),若肝膽失調(diào),膽經(jīng)有熱,易上逆于耳而為病。脾主運(yùn)化而升清,脾氣健旺,氣血充沛,清陽(yáng)之氣上奉耳,則耳的功能正常;若脾失健運(yùn),氣血不足,則耳失所養(yǎng)而失聰;若濕邪困脾,清陽(yáng)不升,濁陰不降,則蒙蔽耳竅而為病。耳與肺亦有一定關(guān)系,在臨床上,耳病初起,往往出現(xiàn)邪氣在表的肺經(jīng)癥狀。耳通過(guò)經(jīng)脈與臟腑和全身廣泛地聯(lián)系,因此有人將耳竅分區(qū),分別隸屬于人體各部,并以此作為耳穴診斷疾病和治療疾病的依據(jù)。

    Δy=[Δy1,Δy2,Δy3…Δy10]T

    (17)

    Xis observation matrix:

    (18)

    thenAis calculated by

    A=(XTX)-1XTΔy

    (19)

    In the above equations, the observation matrixXcan be calculated with the rotational speed and the differential pressure. The observation vector Δywhich represents the torque loss or the flow loss can also be calculated with the test data. Consequently the coefficient matrixAcan be obtained and the complete loss models are established.

    3 Evaluation of the loss models

    3.1 The evaluation criterion

    In this section, the further research is conducted by combining with the test data to discuss the accuracy of the models with different orders. The mean-square-root error (RMSE,Rfor short), the average relative error (ARE,Efor short) and the maximum relative error (MRE,Mfor short) are introduced for evaluation purpose.

    Denotedias the absolute error andεias the relative error between the parameters calculated by test data and model respectively.nrepresents the number of test points.

    RMSE is an important parameter for evaluating the accuracy of regression models. It’s a meaningful measure of dispersion of data.

    (20)

    ARE is the arithmetic average of the relative errors calculated above. It’s defined as

    (21)

    M={εi}max

    (22)

    3.2 Evaluation of the loss models with different orders

    Combined with the test data, RMSE, ARE and MRE are calculated to evaluate the accuracy of the loss models with different orders. The results are shown in Tab.3.

    Tab.3 gives evaluation results of the loss models with different orders. It can be seen that the variation of MRE is much larger than the variation of RMSE and ARE. It demonstrates that the modeling method is able to eliminate the random errors included in the test data and improve the test accuracy.

    Tab.3 Results of evaluation

    Based on Tab.3, a histogram is acquired as shown in Fig.4.

    Fig.4 Results of evaluation

    It can be seen from Tab.3 and Fig.4 that the outputs of the fourth-order and the fifth-order models are close to the original data. Both of their average relative errors are controlled at about 0.1%, less than the third-order model. The mean-square-root error of the fourth-order model is about 0.12% and the maximum relative error is controlled under 0.3%. It demonstrates that the fourth-order model has a good fitting effect in all the test points and meets the requirements of engineering test.

    The volumetric efficiency, the mechanical efficiency and the overall efficiency in any working condition can be predicted with the loss models. To evaluate the loss models with different orders more intuitively, 10 sets of test data are randomly selected. And the relative errors between the predicted efficiency and the true one calculated with the test data are gained as follows.

    It can be seen from Figs.5-7 that the fourth-order and fifth-order models are more accurate than the third-order model. The largest relative errors of the fourth-order and fifth-order models are both controlled under 0.3% while the largest relative error of third-order model is over 0.7%. Considering that the fifth-order model needs more test data to calculate the coefficients, the fourth-order polynomial function is selected for modeling the torque loss and the flow loss of pump. The loss model can be represented by

    y=a1+a2Δp+a3n+a4Δp2+a5Δpn+a6n2+a7Δp3+a8Δp2n+a9Δpn2+a10n3+a11Δp4+a12Δp3n+a13Δp2n2+a14Δpn3+a15n4

    (22)

    Fig.5 Relative error of volumetric efficiency calculated by loss models with different orders

    Fig.6 Relative error of mechanical efficiency calculated by loss models with different orders

    The complete loss models are established and can be used to predict the work characteristics of the pump. The fitting effects of the loss models are shown in Figs.8-9.

    The scatter pointsin Figs.8-9 are original data. The surfaces are drawn according to the loss models. It can be seen that about 98% of the scatter points are located near the surface. Therefore, the fourth-order polynomial function is suitable for modeling the loss of torque and flow in the performance test on pump.

    Fig.7 Relative error of overall efficiency calculated by loss models with different orders

    Fig.8 The loss of torque

    Fig.9 The loss of flow

    4 Application of the modeling method in pump test

    The test results can only be represented with the incomplete performance curves or the characteristic curves with the traditional testing method of pump. Thanks to the modeling method, all the performance parameters in any working condition can be obtained by the loss models. And test results can be easily characterized as the graphs, charts, tables, even 3D graphs.

    As shown in Figs.10-12, the test results are characterized as the 3D graphs on LabVIEW. The graphs contain the torque surface, the input power surface and the overall efficiency surface which can reflect the pump’s work characteristics at various rotational speeds and pressures. And the pump work characteristics in any working condition could be easily obtained by users.

    Fig.10 The torque surface

    Fig.11 The input power surface

    Fig.12 The overall efficiency surface

    Figs.13-15 are the contour charts reflecting pump’s performance, including the equivalent torque curves, the equivalent power curves and the equivalent efficiency curves. With the contour charts, pumps’ work characteristics can be acquired quite easily. The traditional test method only gave several inaccurate curves. Thanks to the modeling method proposed in this paper, the contour charts can be derived from the 3D graphs directly and users can obtain the contour curves at different values as many as they want. Besides that, the modeling method can eliminate the random errors mixed with the test data by surface-fitting. Consequently, the contour curves drawn with the modeling method are much more accurate than ones by the traditional method.

    Fig.13 Equivalent torque curve

    Fig.15 Equivalent efficiency curve

    5 Conclusions

    The modeling method for operating performance test of pump is presented. The loss models for predicting the torque loss and the flow loss are discussed and evaluated. Compared with the traditional test method, the modeling method has some advantages listed in the following.

    First, the number of test points has been reduced to half of the conventional method. A lot of time and energy are saved. The testing efficiency is improved greatly.

    Second, the modeling method can suppress the random errors of the raw test data by surface-fitting with the least square method. And the test accuracy is improved greatly.

    Finally, the work characteristics of the pump in any working condition can be obtained easily. The test results of pump can be characterized as the graphs, charts, tables, even 3D graphs.

    [1] Yang H, Xu H R, Liu C. Flow measurements in a model centrifugal pump by 3-D PIV[J]. IOP Conference Series: Earth and Environmental Science, 2012, 15(6): 562-570.

    [2] ProMinent. API 675 test bed for ProMInent process pumps[J]. World Pumps, 2011, 9: 10.

    [3] Tao Jianfeng, Liu Chengliang, Gu Jianjiang. Improvements to the measurement of electrically controlled hydraulic pumps’ flow /pressure characteristics[J]. Measurement Science and Technology, 2011, 22(6): 354-364.

    [4] Li Zhongshan, Zhang Jintao, Zhang Yaojuan, et al. The development of intelligent measurement and control system about the test and experimental platform of hydraulic pump performance[C]∥ MEC, 2011:576-579.

    [5] Warren Wilson. Positive-displacement pumps and fluid motors[M]. New York: Pitman Publishing Corporation Press, 1950:152-180.

    [6] Shahram Derekhshan, Ahmad Nourbakhsh. Experimental study of characteristic curves of centrifugal pumps working as turbines in different specific speeds[J]. Experimental Thermal and Fluid Science, 2008, 32: 800-807.

    [7] Alexey N Kochevsky. Capabilities of numerical simulation of multiphase flows in centrifugal pumps using modern CFD software[EB/OL]. [2005-09-22].http:∥arxiv:physics/0509193.

    [8] GB/T 17491—2011, Hydraulic fluid power-positive displacement pumps, motors and integral transmissions—Methods of testing and presenting basic steady state performance[S]. Beijing: Chinese Standards Press, 2011.

    [9] Choi S B, Yoo J K, Cho M S, et al. Position control of a cylinder system using a piezoactuator-driven pump[J]. Mechatronics, 2005, 15: 239-249.

    [10] Wu H W, Lee C B. Influence of a relief valve on the performance of a pump/inverter controlled hydraulic motor system[J]. Mechatronics, 1996, 6:1-19.

    [11] Yagiz S. Utilizing rock mass properties for predicting TBM performance in hard rock condition[J]. Tunnelling and Underground Space Technology, 2008, 23: 326-39.

    [12] Zhang Qian, Kang Yilan, Cai Zongxi, et al. A nonlinear multiple regression model for load identification of shield machine combined with mechanical analysis[C]∥FSKD, 2011:1212-1216.

    [13] Li Chunyan, Yu Jihui, Xie Kaigui, et al. Optimal model of congestion management based on non-linear regression[J]. Transactions of China University of Mining & Technology, 2008, 2: 861-863.

    [14] Zhang Hui, Kong Fanxiang, Yu Yang, et al. Combined estrogenic effects of alkyphenolic compounds[J]. Environmental Chemistry, 2009,3:329-33.

    (Edited by Wang Yuxia)

    10.15918/j.jbit1004-0579.201524.0111

    TG 806 Document code: A Article ID: 1004- 0579(2015)01- 0074- 09

    Received 2013- 05- 15

    Supported by the National Science & Technology Pillar Program During the 12th Five-Year Plan Period (2014BAF08B06); the Research Foundation of Beijing Institute of Technology (20110601061)

    E-mail: zhengshuhua@bit.edu.cn

    猜你喜歡
    清陽(yáng)耳病肝氣
    逍遙散加減治療肝氣郁滯型月經(jīng)不調(diào)的臨床效果
    春季養(yǎng)肝食何“味”
    汽車發(fā)動(dòng)機(jī)機(jī)械故障的非接觸式檢測(cè)技術(shù)研究
    耳病治鼻,事出有因
    中老年保健(2021年4期)2021-08-22 07:07:22
    福晉與王爺?shù)膼?ài)情輕喜劇
    女士(2017年8期)2017-08-08 18:39:14
    68例肝氣郁結(jié)型慢性膽囊炎腹腔鏡術(shù)后中醫(yī)護(hù)理的分析研究
    豬藍(lán)耳病的防控
    廣東飼料(2016年2期)2016-12-01 03:43:08
    一例高致病性豬藍(lán)耳病的診治
    高致病性藍(lán)耳病的診斷和治療
    柴胡疏肝散加減治療肝氣犯胃型胃脘痛的療效觀察
    69av精品久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品成人综合色| 国产麻豆成人av免费视频| 国产精品永久免费网站| 欧美国产精品va在线观看不卡| 欧美乱码精品一区二区三区| 伦理电影免费视频| 久久久久国产一级毛片高清牌| 国产极品粉嫩免费观看在线| 免费看十八禁软件| 亚洲精华国产精华精| 99久久精品国产亚洲精品| 啦啦啦韩国在线观看视频| 在线观看免费视频日本深夜| 黄色女人牲交| 老鸭窝网址在线观看| 欧美中文日本在线观看视频| 欧美性长视频在线观看| 男人舔女人下体高潮全视频| av视频在线观看入口| e午夜精品久久久久久久| 一区二区三区国产精品乱码| 99国产极品粉嫩在线观看| 欧美激情极品国产一区二区三区| 国产成人一区二区三区免费视频网站| 男男h啪啪无遮挡| 麻豆久久精品国产亚洲av| 久久中文看片网| 老鸭窝网址在线观看| 99久久综合精品五月天人人| 色婷婷久久久亚洲欧美| av欧美777| 国产欧美日韩一区二区三| 波多野结衣av一区二区av| 成人国产一区最新在线观看| 国产精品综合久久久久久久免费| 女警被强在线播放| 日韩大码丰满熟妇| 香蕉丝袜av| 丝袜美腿诱惑在线| 中出人妻视频一区二区| 白带黄色成豆腐渣| 欧美午夜高清在线| 久久九九热精品免费| 亚洲欧洲精品一区二区精品久久久| 岛国在线观看网站| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 十分钟在线观看高清视频www| 一区二区三区国产精品乱码| 亚洲专区中文字幕在线| 长腿黑丝高跟| 免费在线观看亚洲国产| 国产精品乱码一区二三区的特点| 亚洲国产精品sss在线观看| 少妇裸体淫交视频免费看高清 | 18禁裸乳无遮挡免费网站照片 | 91成人精品电影| 国产精品爽爽va在线观看网站 | 熟女电影av网| 久久精品影院6| 一本久久中文字幕| 黄色片一级片一级黄色片| 国产成人欧美在线观看| 国产成人欧美在线观看| 亚洲天堂国产精品一区在线| 国产精品电影一区二区三区| 久久青草综合色| 亚洲一区二区三区不卡视频| 色哟哟哟哟哟哟| 日本精品一区二区三区蜜桃| 色播亚洲综合网| 99在线视频只有这里精品首页| 久久热在线av| 神马国产精品三级电影在线观看 | 男人的好看免费观看在线视频 | 妹子高潮喷水视频| 琪琪午夜伦伦电影理论片6080| 欧美人与性动交α欧美精品济南到| 伦理电影免费视频| 中文字幕av电影在线播放| 成人欧美大片| 黑人操中国人逼视频| 老司机在亚洲福利影院| 日本黄色视频三级网站网址| 日韩欧美国产一区二区入口| 日韩欧美一区视频在线观看| 国产亚洲欧美在线一区二区| 国产精品久久久久久人妻精品电影| 无人区码免费观看不卡| 在线十欧美十亚洲十日本专区| 一级a爱片免费观看的视频| cao死你这个sao货| 男女那种视频在线观看| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 欧美日韩亚洲国产一区二区在线观看| 国产成+人综合+亚洲专区| 天天躁狠狠躁夜夜躁狠狠躁| av视频在线观看入口| √禁漫天堂资源中文www| 亚洲专区国产一区二区| 亚洲av美国av| 很黄的视频免费| 国产伦在线观看视频一区| 国产成人系列免费观看| a在线观看视频网站| 精品一区二区三区av网在线观看| 日韩欧美 国产精品| 女性生殖器流出的白浆| 性欧美人与动物交配| 国产免费男女视频| 丁香欧美五月| 久久久久久久久久黄片| 18禁国产床啪视频网站| 国产亚洲欧美在线一区二区| 日日干狠狠操夜夜爽| 十八禁人妻一区二区| 国产国语露脸激情在线看| 桃色一区二区三区在线观看| 精品一区二区三区av网在线观看| 国产精品国产高清国产av| 午夜福利在线在线| www.999成人在线观看| 精品欧美国产一区二区三| 99国产精品99久久久久| av欧美777| 一本久久中文字幕| 可以在线观看毛片的网站| 日韩中文字幕欧美一区二区| 国产精品亚洲av一区麻豆| 国产日本99.免费观看| 99久久综合精品五月天人人| 精品久久久久久久久久免费视频| 日本熟妇午夜| 少妇粗大呻吟视频| 亚洲 国产 在线| 日韩三级视频一区二区三区| 这个男人来自地球电影免费观看| 免费在线观看亚洲国产| 国产亚洲精品av在线| 久久久久久九九精品二区国产 | 亚洲全国av大片| 国产一卡二卡三卡精品| 欧美成人性av电影在线观看| 91国产中文字幕| 亚洲色图 男人天堂 中文字幕| 国产97色在线日韩免费| 亚洲国产精品久久男人天堂| 国产伦在线观看视频一区| 免费av毛片视频| 91成年电影在线观看| www国产在线视频色| 一级a爱视频在线免费观看| 深夜精品福利| 久久久久国产一级毛片高清牌| 美女国产高潮福利片在线看| 97碰自拍视频| 国产高清videossex| 国产精品永久免费网站| 中文字幕精品免费在线观看视频| 午夜视频精品福利| 亚洲av熟女| 最近最新中文字幕大全电影3 | 国产午夜精品久久久久久| 国产视频内射| 亚洲国产高清在线一区二区三 | 精品久久久久久久久久免费视频| 国产精品久久久av美女十八| 亚洲五月天丁香| 国产成人av激情在线播放| 国产av在哪里看| 母亲3免费完整高清在线观看| 九色国产91popny在线| 91字幕亚洲| 性色av乱码一区二区三区2| 丰满人妻熟妇乱又伦精品不卡| 国产高清视频在线播放一区| 99久久无色码亚洲精品果冻| 亚洲av成人一区二区三| 久久天堂一区二区三区四区| 国产亚洲欧美98| 亚洲精品一卡2卡三卡4卡5卡| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 男女之事视频高清在线观看| 精品国内亚洲2022精品成人| 精品国产一区二区三区四区第35| 国产精品美女特级片免费视频播放器 | av视频在线观看入口| 国产精品一区二区免费欧美| 国产不卡一卡二| 国产欧美日韩一区二区精品| 日韩国内少妇激情av| 午夜福利18| 90打野战视频偷拍视频| 别揉我奶头~嗯~啊~动态视频| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 男女视频在线观看网站免费 | 亚洲成人久久性| 日本免费a在线| 天堂影院成人在线观看| 精品人妻1区二区| 丝袜人妻中文字幕| 免费在线观看视频国产中文字幕亚洲| 啦啦啦 在线观看视频| 亚洲欧美精品综合久久99| 国产一级毛片七仙女欲春2 | 免费在线观看影片大全网站| 国产精品 国内视频| 中文字幕av电影在线播放| 国产精品香港三级国产av潘金莲| 97超级碰碰碰精品色视频在线观看| 少妇 在线观看| 成人亚洲精品一区在线观看| 91麻豆av在线| 午夜视频精品福利| 午夜福利在线观看吧| 久久久久久久久中文| 国产蜜桃级精品一区二区三区| 热99re8久久精品国产| 国产精品一区二区精品视频观看| 国产精品日韩av在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美亚洲日本最大视频资源| 久久久久久久精品吃奶| 国产伦在线观看视频一区| 深夜精品福利| www.熟女人妻精品国产| 亚洲国产精品999在线| 大香蕉久久成人网| 亚洲最大成人中文| www.999成人在线观看| 成人亚洲精品av一区二区| 日韩免费av在线播放| 最近最新免费中文字幕在线| 亚洲最大成人中文| 美女高潮到喷水免费观看| www.熟女人妻精品国产| av在线播放免费不卡| 国产极品粉嫩免费观看在线| 亚洲精品国产区一区二| 午夜福利在线观看吧| 女人爽到高潮嗷嗷叫在线视频| 亚洲av美国av| 久久久久精品国产欧美久久久| 国产av在哪里看| 91国产中文字幕| 757午夜福利合集在线观看| 美女高潮到喷水免费观看| 国产片内射在线| 91成人精品电影| 好男人在线观看高清免费视频 | 欧美色视频一区免费| 成年版毛片免费区| 极品教师在线免费播放| 国产激情偷乱视频一区二区| 亚洲男人天堂网一区| 亚洲av电影不卡..在线观看| 亚洲欧美日韩无卡精品| 日本撒尿小便嘘嘘汇集6| 亚洲国产高清在线一区二区三 | 人人妻人人看人人澡| 亚洲,欧美精品.| 老司机福利观看| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频| 久久青草综合色| 黄色女人牲交| 欧美中文综合在线视频| 国产精品久久久av美女十八| x7x7x7水蜜桃| ponron亚洲| 在线观看午夜福利视频| 欧美不卡视频在线免费观看 | 非洲黑人性xxxx精品又粗又长| 精品少妇一区二区三区视频日本电影| 男女视频在线观看网站免费 | 国产野战对白在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩中文字幕欧美一区二区| 精品久久久久久久末码| 老熟妇乱子伦视频在线观看| 搡老妇女老女人老熟妇| 久久久久免费精品人妻一区二区 | 精品久久久久久久久久久久久 | 午夜老司机福利片| 老汉色∧v一级毛片| 麻豆av在线久日| 国内精品久久久久精免费| 自线自在国产av| 成人午夜高清在线视频 | 黄色视频,在线免费观看| 一二三四社区在线视频社区8| 日本三级黄在线观看| 久久久国产欧美日韩av| 国产亚洲欧美98| 久久久久久九九精品二区国产 | 日韩欧美在线二视频| 国产一级毛片七仙女欲春2 | 两个人看的免费小视频| 久久国产精品影院| 国内揄拍国产精品人妻在线 | 级片在线观看| 亚洲国产欧美日韩在线播放| 叶爱在线成人免费视频播放| 国产亚洲欧美精品永久| 在线观看免费视频日本深夜| 国产三级黄色录像| 欧美丝袜亚洲另类 | 老汉色av国产亚洲站长工具| 欧美日韩福利视频一区二区| 精品卡一卡二卡四卡免费| 国产真人三级小视频在线观看| 亚洲五月婷婷丁香| 女同久久另类99精品国产91| 日韩大码丰满熟妇| 午夜两性在线视频| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 香蕉丝袜av| 亚洲av片天天在线观看| 亚洲无线在线观看| 人人妻人人看人人澡| 啦啦啦 在线观看视频| 黄片播放在线免费| 免费在线观看完整版高清| 国产成人啪精品午夜网站| 精华霜和精华液先用哪个| xxxwww97欧美| 999久久久国产精品视频| 给我免费播放毛片高清在线观看| 国产一级毛片七仙女欲春2 | 99国产极品粉嫩在线观看| 在线观看免费午夜福利视频| 午夜福利高清视频| 无人区码免费观看不卡| 亚洲精品在线观看二区| 国产区一区二久久| 2021天堂中文幕一二区在线观 | 美女午夜性视频免费| 夜夜夜夜夜久久久久| 日韩欧美国产一区二区入口| 99久久精品国产亚洲精品| 久久欧美精品欧美久久欧美| 非洲黑人性xxxx精品又粗又长| 老熟妇仑乱视频hdxx| 久久久国产成人免费| 亚洲精品中文字幕在线视频| 99精品久久久久人妻精品| 黑人巨大精品欧美一区二区mp4| 国产一区二区激情短视频| 国产激情久久老熟女| 美女午夜性视频免费| 中出人妻视频一区二区| 中文字幕人成人乱码亚洲影| 亚洲av电影不卡..在线观看| 亚洲第一av免费看| 国产亚洲av嫩草精品影院| 日本免费a在线| 欧美人与性动交α欧美精品济南到| 色播在线永久视频| 午夜a级毛片| 非洲黑人性xxxx精品又粗又长| 色婷婷久久久亚洲欧美| 长腿黑丝高跟| xxxwww97欧美| 日本黄色视频三级网站网址| 国产一区二区三区在线臀色熟女| 天堂影院成人在线观看| 精品人妻1区二区| 嫁个100分男人电影在线观看| 中文字幕高清在线视频| 成人三级黄色视频| 亚洲五月天丁香| 真人做人爱边吃奶动态| 久久欧美精品欧美久久欧美| 欧美zozozo另类| 真人做人爱边吃奶动态| 他把我摸到了高潮在线观看| 人人妻人人澡人人看| 国产精品国产高清国产av| 国产成年人精品一区二区| 国产精品亚洲av一区麻豆| 少妇粗大呻吟视频| 久久精品成人免费网站| 中文资源天堂在线| 18禁国产床啪视频网站| aaaaa片日本免费| 国产精品一区二区精品视频观看| 久久久水蜜桃国产精品网| 成年免费大片在线观看| 又紧又爽又黄一区二区| 免费高清在线观看日韩| www.熟女人妻精品国产| 午夜影院日韩av| 午夜福利视频1000在线观看| www日本在线高清视频| 亚洲国产欧洲综合997久久, | 亚洲国产精品合色在线| 免费av毛片视频| 久久午夜亚洲精品久久| 一级毛片精品| 精品人妻1区二区| 黄色视频,在线免费观看| 老司机午夜福利在线观看视频| 精品免费久久久久久久清纯| 91大片在线观看| 亚洲中文字幕日韩| 美女高潮到喷水免费观看| 国产日本99.免费观看| 不卡av一区二区三区| avwww免费| 大香蕉久久成人网| 91在线观看av| 亚洲av美国av| 亚洲成国产人片在线观看| 日日夜夜操网爽| 熟女电影av网| 在线观看一区二区三区| 露出奶头的视频| 欧美乱妇无乱码| 国产精品永久免费网站| av电影中文网址| 老汉色∧v一级毛片| 精品国产一区二区三区四区第35| 热re99久久国产66热| 亚洲午夜精品一区,二区,三区| 亚洲无线在线观看| 色哟哟哟哟哟哟| 国产精品一区二区三区四区久久 | 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 国内久久婷婷六月综合欲色啪| 国产精品亚洲美女久久久| 亚洲 国产 在线| 国产欧美日韩精品亚洲av| 免费看十八禁软件| 美女扒开内裤让男人捅视频| 亚洲一区二区三区不卡视频| 亚洲专区中文字幕在线| 中文字幕人妻丝袜一区二区| 精品久久久久久久末码| 精品午夜福利视频在线观看一区| 国产精品久久久久久亚洲av鲁大| 国产亚洲av嫩草精品影院| 国产乱人伦免费视频| 一夜夜www| 中文资源天堂在线| 国产又爽黄色视频| 亚洲免费av在线视频| 精品国产乱子伦一区二区三区| 又紧又爽又黄一区二区| 一区福利在线观看| 最近最新中文字幕大全免费视频| 两人在一起打扑克的视频| 亚洲中文日韩欧美视频| 18禁观看日本| 老汉色av国产亚洲站长工具| 午夜精品在线福利| 亚洲在线自拍视频| 日韩中文字幕欧美一区二区| 国产蜜桃级精品一区二区三区| 国产精品久久久久久人妻精品电影| 精品久久蜜臀av无| 免费在线观看视频国产中文字幕亚洲| 国产男靠女视频免费网站| 久久久精品欧美日韩精品| 三级毛片av免费| 亚洲中文日韩欧美视频| 亚洲av成人不卡在线观看播放网| 国产精品日韩av在线免费观看| 亚洲成人久久爱视频| 黄色毛片三级朝国网站| 18禁观看日本| 国产精品爽爽va在线观看网站 | 久久久久免费精品人妻一区二区 | 国产成人精品久久二区二区91| 特大巨黑吊av在线直播 | 十八禁网站免费在线| 国产色视频综合| 我的亚洲天堂| 国产爱豆传媒在线观看 | 国产一区在线观看成人免费| 亚洲欧美一区二区三区黑人| 91麻豆av在线| 婷婷精品国产亚洲av| 国产私拍福利视频在线观看| 禁无遮挡网站| 可以在线观看毛片的网站| 成人国产一区最新在线观看| 国产一区二区三区在线臀色熟女| 97超级碰碰碰精品色视频在线观看| 国产亚洲精品一区二区www| 在线视频色国产色| 免费搜索国产男女视频| 精品久久久久久久久久久久久 | 免费看a级黄色片| 亚洲精品中文字幕一二三四区| 国产激情久久老熟女| 精品不卡国产一区二区三区| 深夜精品福利| 国产成人av教育| 99国产精品一区二区三区| 日本 av在线| 91国产中文字幕| 真人做人爱边吃奶动态| 午夜福利18| 欧美成人一区二区免费高清观看 | 宅男免费午夜| 亚洲天堂国产精品一区在线| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o| 高清在线国产一区| 久久婷婷成人综合色麻豆| a级毛片a级免费在线| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 国产三级黄色录像| www国产在线视频色| 女同久久另类99精品国产91| 国产黄a三级三级三级人| 最新在线观看一区二区三区| 婷婷亚洲欧美| 免费高清视频大片| 亚洲免费av在线视频| 亚洲精品一卡2卡三卡4卡5卡| АⅤ资源中文在线天堂| 国产区一区二久久| 变态另类成人亚洲欧美熟女| 亚洲精品粉嫩美女一区| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜精品一区,二区,三区| 88av欧美| 性欧美人与动物交配| 国产视频一区二区在线看| videosex国产| 久久草成人影院| 999久久久精品免费观看国产| 成人亚洲精品一区在线观看| 老司机福利观看| 每晚都被弄得嗷嗷叫到高潮| 男女之事视频高清在线观看| 一本综合久久免费| 日本a在线网址| 色老头精品视频在线观看| 一二三四社区在线视频社区8| 国产成人欧美在线观看| 日本黄色视频三级网站网址| 国产亚洲av嫩草精品影院| 一进一出抽搐gif免费好疼| 男人舔女人的私密视频| 亚洲人成电影免费在线| 神马国产精品三级电影在线观看 | 亚洲熟妇中文字幕五十中出| av片东京热男人的天堂| 国产亚洲精品久久久久久毛片| 国产一区二区三区在线臀色熟女| 精品一区二区三区av网在线观看| 亚洲久久久国产精品| 一级毛片高清免费大全| 丝袜美腿诱惑在线| 女性被躁到高潮视频| 亚洲欧美一区二区三区黑人| 国产成+人综合+亚洲专区| 国产精品自产拍在线观看55亚洲| 脱女人内裤的视频| 亚洲国产欧洲综合997久久, | 18禁黄网站禁片免费观看直播| 两人在一起打扑克的视频| 国产真人三级小视频在线观看| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 满18在线观看网站| 国产亚洲精品综合一区在线观看 | 欧美+亚洲+日韩+国产| 中国美女看黄片| 国产伦一二天堂av在线观看| 波多野结衣高清作品| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 欧美黑人巨大hd| 国产熟女xx| 免费高清视频大片| 国产精品 欧美亚洲| 国产三级在线视频| 麻豆成人午夜福利视频| 国产精品98久久久久久宅男小说| 国产精品自产拍在线观看55亚洲| 国产精品久久电影中文字幕| 国产精品影院久久| 1024手机看黄色片| 久久久久久九九精品二区国产 | 丝袜人妻中文字幕| 欧美国产精品va在线观看不卡| 黄频高清免费视频| 一本综合久久免费| 欧美日韩亚洲综合一区二区三区_| 色综合婷婷激情| 免费高清在线观看日韩| 国产亚洲精品综合一区在线观看 | 日韩欧美免费精品| 熟妇人妻久久中文字幕3abv| 免费女性裸体啪啪无遮挡网站| 可以在线观看毛片的网站| 亚洲av美国av| 亚洲国产精品久久男人天堂| 美女大奶头视频| 99热6这里只有精品| 成人特级黄色片久久久久久久| 免费人成视频x8x8入口观看| 亚洲一区二区三区不卡视频|