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

    Vibration response analysis of a lathe spindle by using the ANSYS finite element method

    2015-11-03 07:08:16ZexinZHOU
    機床與液壓 2015年3期

    Ze-xin ZHOU

    (School of Energy and Power Engineering, Wuhan University of Technology, Wuhan 430063, China)

    ?

    Vibration response analysis of a lathe spindle by using the ANSYS finite element method

    Ze-xin ZHOU*

    (School of Energy and Power Engineering, Wuhan University of Technology, Wuhan 430063, China)

    Making an analysis for vibration modal and frequency response of the lathe spindle, respectively by using finite element method based on ANSYS and experiment of CA6140 type lathe in machining, and the calculation results are compared and analyzed, which verified the accuracy of ANSYS method. Numerical simulation and experimental results show that: Spindle in the first order and fifth order are prone to resonance, but did not reach resonance, the low order natural frequency have more effect than the high order natural frequency of the spindle vibration; by the experiments can conclude that the maximum vibration of the main shaft in the working state is mainly concentrated in the vicinity of its two ends, therefore, the improved bearing is an important way to reduce the vibration of the main shaft and ensure the machining accuracy, and the research results can provide a theoretical reference for the structural optimization design of the lathe.

    ANSYS FEM, Lathe spindle, Modal analysis, Harmonic response analysis, Resonance

    1 Introduction

    With the rapid development of mechanical processing in high speed, high precision and high flexibility, and in addition to the requirements for machine tool have features that light quality, low cost, easy to use and also has a good process performance, and calls for machine tools has become more and more high processing performance. We all know that the dynamic characteristic is the most important performance index of the evaluation of machine tools [1]. But most design of ordinary lathe follows the usual concept, the theory basis of machine tool design deficiencies [2]. In the component part of the whole machine tool, the spindle is an important part, which drives the rotation of the chuck and work piece, and in the process of turning processing, often have a strong vibration between work piece and tool, damage and interfere with the normal cutting processing, which is a very harmful phenomenon. When the lathe vibrates, the surface quality of work piece is deteriorated, and it has obvious surface marks, which increasing the roughness of the work piece, at this time must be reduced to the working efficiency of the lathe greatly. When the lathe Strong vibrates, can cause collapse edge phenomenon, and make the machining process cannot continue. Due to the vibration, will make the lathe and the tool wear and to shorten the lathe and the cutting tool service life. At the same time vibration is also accompanied by the noise, making harms to workers health of body and mind, and making work environment deterioration. Lathe vibrations can be classified into three kinds, which are free vibration, forced vibration and the Self system vibration, according to estimates that of vibration is respectively 5%, 30%, 65% [3-4]. Since the free vibration, which motor driving occupies a large part of the main shaft vibration, so the dynamic analysis of the main shaft is necessary. In this paper, we making an analysis for vibration modal and harmonic response of the lathe spindle based on ANSYS finite element method of CA6140 type lathe in machining and its technical parameters as shown in Table 1, through the research of the natural frequency, vibration mode, harmonic response and other dynamics research to provide a reference for optimization design of this type machine tool.

    Table 1 Main technical parameters of CA6140 lathe

    Amaximumturningdiameter/mm210Spindlecenterofbedflatguidetodis-tance/mm205Mainmotorpower/kW7.5Motorpowerquickly/kW250Spindlespeed/(r·min-1)0-3000Spindleaperture/mm48

    2 Theoretical analysis

    2.1 The foundation theory of modal analysis

    Any deformable object has the inherent frequency and the vibration mode, and they are the important parameters in the structure design under the dynamic load. The main purpose of the modal analysis is to determine the vibration characteristics of the structure, i.e., the natural frequency, the vibration mode and the mode participation coefficient of vibration (the degree of vibration in a particular direction), importantly they are the frequency feature points of the structure that produce resonance. Calculating modal analysis based on the geometric features of the mechanical structure and material properties such as the original parameters, form the system of discrete mathematical model with finite element method-mass matrix and stiffness matrix, and then by solving eigenvalue problem to determine the modal parameters of the system [5]. Undamped modal analysis is a classical eigenvalue problem, the dynamics equations of motion can be described as:

    (1)

    Where:[M] is mass matrix;[K] is stiffness matrix.

    The free vibration of Structural is harmonic vibration, the displacement of sine function can be described as:

    (2)

    Taking (2) into (1) :

    (3)

    Modalanalysisactuallyissolvingtheeigenvalueandeigenvector,andthroughthemodalanalysiscandeterminethemodelorthevibrationofthestructurecharacteristics,andcanalsogetthemodelofnaturalfrequencies,modeshapesandmodalparticipationfactor(inaparticulardirectionvibrationmodetheextenttowhichinvolvedinthevibration).

    2.2Thefoundationtheoryofharmonicresponse

    Anycontinuousperiodicloadwillgenerateacontinuousperiodicresponseinthesystem,i.e.,theharmonicresponse,harmonicresponseanalysisisalsoknownasthefrequencyresponse,whichismainlyanalysisthesteadystateresponseoflinearstructureinunderharmonicload(withtimeregularchangeinharmonicload).Themainpurposeoftheanalysisistogetthestructureresponsecurveinseveralfrequencies,undertheharmonicload,andfoundthepeakresponseofthestructure,sothestressresponseofthestructureisfurtherobserved[6].Theharmonicresponseanalysisislinearwhichonlycalculatestheresponseofthestructuretotheforcedvibration,withoutconsideringthetransientvibrationatthestartoftheexcitation.Harmonicresponseanalysisofthemainaimistoensurethatthedesignofthestructurecanwithstandavarietyofdifferentfrequencyharmonicload,detectingstructuralresonantresponse,andtoavoidresonancewhennecessary.

    Theharmonicresponseanalysisisthevibrationresponseofthestructurewhentheexternalforceformistheharmonicfunction(i.e.thefunctionformissinorcos),theexpressionis[7]

    (4)

    Where:

    Where:

    Fmax—Load amplitude;ψ—Phase angle of load function;F1=Fmaxcosψ;F2=Fmaxsinψ;

    Dmax—Displacement amplitude;D1=Dmaxcosψ;D2=Dmaxsinψ;

    Taking (5)and (6)into(3)can get the equations of motion for harmonic analysis.

    (7)

    To solve the harmonic motion equation (7), there are two methods, one is a complete method, the other is a modal superposition method.

    In ANSYS, using the complete method is a default method, using complete structural matrix, and allows the asymmetric matrix.

    In the complete method, the matrix equation is directly solved in the nodal coordinate system, except for the use of the complex number, which is similar to the linear static analysis.

    Order:

    (8)

    There is

    (9)

    3 The ANSYS finite element analysis

    3.1 Modeling and material properties

    In the environment of Solidworks, according to the similarity theory, the model of the lathe spindle is simplified as shown in Fig.1 (a), main shaft length is 870 mm, and the length is divided into eight parts, the diameter is 80-195 mm, central hole of spindle is 48 mm, Sectional view of the model is shown in Fig.1 (b).

    Fig.1 3 d model of Main shaft

    According to the actual lathe spindle material model entity set to 45 steel, the performance of the material parameters are shown in Table 2.

    Figures are the density, Poisson’s ratio and Young’s modulus of 45 steel with temperature variation, we can find that it is not affected by temperature.

    Table 2 Material parameters of 45 steel

    DensityPoisson’sratioYoung’smodulusShearmodulusBulkmodulus0.785×1040.32×10117.692×10101.667×1011

    Fig.2 Density with the temperature variation

    Fig.3 Poisson’s ratio with the temperature variation

    Fig.4 Young’s elastic modulus with temperature variation

    3.2 The modal analysis

    The material properties of the model are set up, in the meshing, mesh with the method is insert sizes, the size of the grid size is 0.005 m, a total of 515381 nodes, the cell number is 357554. To differentiate good model grid, as shown in Fig.5, we examined the grid quality which are good, and there was no big distortion phenomenon.

    Fig.5 Meshing figure

    For order selection, the study found that the vibration system with more than 1 000 degrees of freedom, to extract the first three order vibration mode, its precision can achieve 90% above, if extract of the first ten order modes, its accuracy can reach above 99% [8], so for the convenience of calculation, for some higher order modes can be neglected, mainly analyzes its low-order modes. In the process of calculation, extraction of the first six order vibration mode were analyzed, and modal analysis after the first six order natural frequency as shown in Table 3.

    Table 3 Natural frequency of the first six order

    Ordertime123Frequency/Hz562.83562.851384.4Ordertime456Frequency/Hz1384.517382460.4

    The first order frequency of main shaft is 562.83 Hz, and vibration mode is shown in Fig.6, as can be seen from its vibration mode: the torsional vibration of main shaft occurred along theXaxis, the biggest variable occurs in the central spindle. The second order frequency of main shaft is 562.85 Hz,and vibration mode is shown in Fig.7, as can be seen from its vibration mode: main shaft torsional vibration occurred along theYaxis, the biggest variable also occurs in the central spindle, and the main shaft of first and second order frequency is closed to. Third order frequency is 1 384.4 Hz, vibration mode is shown in Fig.8, as can be seen from its vibration mode: spindle happened along theXaxis bending vibration, the largest occurs in the central spindle shape variables on the right side. The forth order frequency is 1 384.5 Hz, and vibration mode is shown in Fig.9, as can be seen from its vibration mode: happened to the main shaft bending vibration along theYaxis, the biggest right form variables also occurs in the central spindle, and the main shaft of third-order and fourth-order frequency is closed to. The fifth order frequency is 1 738 Hz, vibration mode as shown in Fig.10, as can be seen from its vibration mode: spindle happened along theZaxis transverse compression vibration, maximize maximum shape variables in the main shaft diameter. The sixth order frequency is 2 460.4 Hz, and vibration mode is shown in Fig.11, as can be seen from its vibration mode: spindle happened along theXaxis bending vibration, the biggest variable occurs in the central spindle right. In the first six order modal, the maximum deformation is the sixth order.

    Fig.6 The first-order modal

    Fig.7 The second-order modal

    Fig.8 The third-order modal

    Fig.9 The forth-order modal

    Fig.10 The fifth-order modal

    Fig.11 The sixth-order modal

    3.3 The harmonic response analysis

    On the basis of the modal analysis and the harmonic response analysis was carried out on the spindle work condition by using ANSYS [9]. Because the spindle of the highest frequency is 2 460.4 Hz, here set the harmonic response frequency range is 0-1 640 Hz, harmonic response analysis need to impose a sine vibration force, here in the middle of the main force set 250 N, and the displacement response harmonic response curves analysis, as shown in Fig.12 and the normal stress response curve as shown in Fig.13, from curves we can see that harmonic response in the first order and second order is maximum, and deformation phase is 180 degrees, in the third and fourth order have the minimum deformation. So in the first and second order spindle are prone to resonance, and the low order natural frequency is higher than the natural frequency of the vibration of the spindle and the influence is large, the lower the frequency, the more violent vibration, the impact on the structure is large.

    4 Experiments and analysis

    4.1 Experiment systems

    The modal experiment system is composed of 4 parts, which are the object, the excitation system, the data acquisition system and the data analysis system. In this experiment, the main shaft quality and the size are small, and the vibration device in the same number of magnitude, so the experimental use of non-stationary excitation system excitation device, which is exciting-hammer composed of hammer and sensor, and the input signal generated by excited vibration is pulse signal. Data acquisition system includes sensor, charge amplifier and signal acquisition system, in which the acceleration sensor will be collected signal by the charge amplifier, and the force on the force sensor to the signal collected together with the multi-channel data acquisition instrument for processing. Data acquisition instrument transfer data to a PC, using modal testing analysis software to estimate the frequency response function, to identify the measured object modal parameters and the modal change of structure. The modal test equipment of the spindle is shown in Table 4.

    Fig.12 Displacement of the frequency response curve

    Fig.13 Under the action of normal stress response curve

    Table 4 Modal test equipment of the spindle

    EquipmentTypeExciting-hammerB&K8206-002PiezoelectricaccelerationsensorB&K4326A001ChargeamplifierB&KZG0429Multi-channeldataacquisitioninstrumentLMSSC305-UTPTestanalysissoftwareLMSTest.labRev6v

    4.2 Experiment analysis

    Modal experiment is the first to choose the correct object support, this experiment adopted free support and used with elastic rubber tube spindle hanging up, making the spindle in the free state approximation, with most of the degrees of freedom. In the process of the experiment, the hammer is used to knock on the main shaft, and a sensor is used for reference, and by moving the other sensors to measure the signal of each measurement point. In the experiment, the 3 measuring points are arranged in the vicinity of the front and back ends of the spindle and other 3 measuring points is arranged in the middle of the spindle, and a total of 9 measuring points are arranged in the total spindle. In the analysis of experimental results, the first response in two exciting point all the frequency response function of overall goodness of fit, which can show a spindle inherent resonance frequency more accurately, the modal of the spindle measured by experiment is shown in Table 5.

    Table 5 Modal parameters of experiment

    OrderFrequency/HzDampingratio/%1556.81.052559.41.1231380.51.2841380.91.2951710.61.3262400.61.45

    Through the experiment that spindle modal as shown in Table 5, and Table 3 is calculated by the ANSYS finite element method, both error analysis as shown in Fig.14. Through the contrast analysis diagram can be seen the error is very small, and to verify the accuracy of the finite element method.

    Because the vibration of the spindle changed with the motor variation speed, which belongs to the inherent characteristics, and has a smaller association with grinding state, therefore, the vibration experiment of the spindle with the variation rotating speed is also studied. In this paper, the measuring experiment is carried out with 360 r/min and 3000 r/min of two rotating speed, and the time domain signal of the maximum vibration displacement of the bearing area at both ends of the front and rear ends are shown in Fig.15 and Fig.16, and the frequency domain analysis of the bearing are shown in Fig.17, Fig.18.

    Fig.14 Error analysis chart

    Fig.15 Time domain signal at 360 r/min

    Fig.16 Time domain signal at 3 000 r/min

    Fig.17 Frequency domain analysis at 360 r/min

    Fig.18 Frequency domain analysis at 3 000 r/min

    As the accuracy and rigidity of the main shaft has not been unified the international standards, so the experimental reference to the German manufacturers of well-known manufacturers of GMN precision acceptance criteria as a reference. For the vibration of grinding spindle head with cutting tool and using ceramic ball bearing, the test standard of GMN company is less than 2 μm, but the main shaft used in this experiment did not install the tool, and the use of the ball bearing is still greater than 2 μm, from its spectrum also can be seen in the first and fifth order near the vibration. From the experimental data, the domestic production of precision spindle with foreign advanced level, there are still large compared to the gap, the main reason might be the spindle bearing, the manufacturing and assembly process and dynamic balance technology in three aspects.

    5 Conclusions

    This paper made an analysis for vibration modal and frequency response of the lathe spindle, by using finite element method based on ANSYS and the experiment respectively, through the research can get the following conclusion.

    1) Through the modal analysis of the main shaft is obtained the 6 order natural frequency, the natural frequency of the first order and the second order is closed to, and the vibration is mainly torsional. The natural frequency of the third order and fourth order is also closed to, and the vibration is mainly bending.

    2) Through the harmonic response analysis of the spindle can be obtained easily in one order and two order at the occurrence of resonance, and the low order natural frequency influence than the oscillation of high order natural frequencies of the spindle, the lower the frequency, the vibration is more intensive, which have large influence on the structure, So the work speed of the spindle system as far as possible to avoid first and second order critical speed, to improve reliability and service life of the spindle.

    3) The experimental results show that the maximum vibration of the main shaft in the working state is mainly concentrated in the area near the ends of the bearing. So the improvement of bearing is an important way to reduce the vibration of the main spindle and ensure the machining accuracy.

    Acknowledgement

    This work is financially supported by Independent Innovation Research Fund of Wuhan University of Technology(No. 2014-ND-B1-09).

    [1]FENG chengguo, CAO jujiang, ZHANG lei. Research on dynamic character of milling machining center worktable based on ANSYSWorkbench[J]. Hydromechatronics Engineering, 2013(4):175-177.

    [2]YANG yongliang. The structural optimization of Lathe-Bed based on finite element [D]. Dalian: Dalian University of Technology, 2006.

    [3]YANG hao, QIN ping, LI bolin, et al. Modle analysis of floor type boring and milling machine tool based on ANSYS[J]. Hydromechatronics Engineering, 2010(05):91-93.

    [4]LI yaofang, LI duanneng, LIU zhi, et al. Modal analysis of cross platform on using hydrostatic guide-way hard turning machine tool[J]. Hydromechatronics Engineering, 2012(12): 21-24.

    [5]YANG xiaowen. Modal Analysis of Vibration of Rolling Washing Machine[D].Shanghai: Shanghai Jiaotong University.

    [6]Yang mingya, Yang yingjie. Dynamic performance analysis of high speed motorized spindle box[J]. Digital design, 2009(5): 68-70.

    [7]Lv jianfa. Research on dynamic characteristics of NC boring and milling machine[D]. Chengdu: Xinan Jiaotong University, 2010.

    [8]MAN jia, XU yanshen, ZHANG xueling. Modal analysis of vertical lathe beam based on ANSYS workbench[J]. Combination machine tools and automatic processing technology, 2006(6): 45-49.

    [9]YAO Xiaobo, LIU deping, GAO Jianshe, et al. Dynamic analysis of machining center based on ANSYS Workbench[J]. Hydromechatronics Engineering, 2013(8):164-167.

    (Continued on 69 page)

    25 November 2014; revised 21 March 2015;

    Ze-xin ZHOU, E-mail: 351268909@qq.com

    10.3969/j.issn.1001-3881.2015.18.011 Document code: A

    TB533+.1

    accepted 2 July 2015

    Hydromechatronics Engineering

    http://jdy.qks.cqut.edu.cn

    E-mail: jdygcyw@126.com

    精品卡一卡二卡四卡免费| 在线观看www视频免费| 免费少妇av软件| 成人亚洲精品一区在线观看| 最黄视频免费看| 国产精品欧美亚洲77777| 欧美 亚洲 国产 日韩一| 久久久久久久国产电影| 精品少妇黑人巨大在线播放| 亚洲天堂av无毛| 国产成人啪精品午夜网站| 中文字幕精品免费在线观看视频| 国产日韩欧美亚洲二区| 色94色欧美一区二区| 亚洲色图 男人天堂 中文字幕| 黑人猛操日本美女一级片| 婷婷色综合www| 国产深夜福利视频在线观看| 精品国产一区二区三区久久久樱花| 大香蕉久久成人网| 亚洲欧美日韩高清在线视频 | 国产麻豆69| 校园人妻丝袜中文字幕| 晚上一个人看的免费电影| 欧美精品高潮呻吟av久久| 乱人伦中国视频| 日韩中文字幕视频在线看片| 狂野欧美激情性xxxx| 熟女少妇亚洲综合色aaa.| 色婷婷av一区二区三区视频| 国产极品粉嫩免费观看在线| 精品国产乱码久久久久久小说| 一级片免费观看大全| 亚洲精品日韩在线中文字幕| 欧美成狂野欧美在线观看| 婷婷色综合大香蕉| 大片电影免费在线观看免费| 操出白浆在线播放| 久久这里只有精品19| 精品久久蜜臀av无| 国产高清国产精品国产三级| 国产日韩一区二区三区精品不卡| 亚洲欧美日韩高清在线视频 | 亚洲欧美成人综合另类久久久| 黄色怎么调成土黄色| 亚洲一区中文字幕在线| 熟女少妇亚洲综合色aaa.| 黑丝袜美女国产一区| 少妇被粗大的猛进出69影院| 久热这里只有精品99| 国产伦理片在线播放av一区| 中文字幕人妻丝袜一区二区| 深夜精品福利| 成年人免费黄色播放视频| 看十八女毛片水多多多| 肉色欧美久久久久久久蜜桃| 最近最新中文字幕大全免费视频 | 高潮久久久久久久久久久不卡| 久久国产亚洲av麻豆专区| 婷婷色综合www| 亚洲精品成人av观看孕妇| 男人舔女人的私密视频| 久久狼人影院| 交换朋友夫妻互换小说| 青春草亚洲视频在线观看| 一级黄片播放器| 欧美变态另类bdsm刘玥| 丝袜人妻中文字幕| 一级毛片女人18水好多 | 日韩一卡2卡3卡4卡2021年| 久久国产亚洲av麻豆专区| 嫩草影视91久久| 亚洲成色77777| 精品国产一区二区久久| 亚洲国产欧美网| 男女午夜视频在线观看| 国产女主播在线喷水免费视频网站| 国产男人的电影天堂91| 侵犯人妻中文字幕一二三四区| 国产一卡二卡三卡精品| 日本一区二区免费在线视频| 亚洲精品国产av成人精品| 91成人精品电影| 免费观看av网站的网址| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 精品一品国产午夜福利视频| 亚洲美女黄色视频免费看| 人人妻人人添人人爽欧美一区卜| 丝瓜视频免费看黄片| 国产高清不卡午夜福利| 国产男女内射视频| 精品国产超薄肉色丝袜足j| 91精品三级在线观看| 日本午夜av视频| 精品高清国产在线一区| 亚洲国产毛片av蜜桃av| 这个男人来自地球电影免费观看| 一区二区三区四区激情视频| 成人国产一区最新在线观看 | 亚洲欧美中文字幕日韩二区| 热99久久久久精品小说推荐| 操出白浆在线播放| 欧美性长视频在线观看| 久热这里只有精品99| 欧美亚洲 丝袜 人妻 在线| 久久热在线av| 亚洲精品一区蜜桃| 久久人妻熟女aⅴ| 免费久久久久久久精品成人欧美视频| 婷婷色av中文字幕| 视频区欧美日本亚洲| 少妇人妻久久综合中文| 女人爽到高潮嗷嗷叫在线视频| av在线老鸭窝| 久久99精品国语久久久| 国产99久久九九免费精品| 亚洲精品国产av蜜桃| 久久久国产精品麻豆| 最新在线观看一区二区三区 | 午夜福利免费观看在线| 一级片'在线观看视频| 亚洲专区国产一区二区| 国产在线免费精品| 久久久久国产一级毛片高清牌| 女性被躁到高潮视频| 亚洲视频免费观看视频| 国产日韩欧美亚洲二区| 亚洲视频免费观看视频| 日本wwww免费看| 中文字幕精品免费在线观看视频| 亚洲,欧美精品.| 69精品国产乱码久久久| 国产精品久久久久成人av| 成人午夜精彩视频在线观看| 在线观看免费午夜福利视频| 热99国产精品久久久久久7| 久热这里只有精品99| 成年人黄色毛片网站| 亚洲中文日韩欧美视频| 在线观看免费视频网站a站| 十八禁人妻一区二区| 久久99热这里只频精品6学生| 亚洲 国产 在线| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 9热在线视频观看99| 熟女av电影| 亚洲人成电影免费在线| 99国产精品一区二区三区| 考比视频在线观看| 男人操女人黄网站| 免费在线观看视频国产中文字幕亚洲 | 曰老女人黄片| 五月开心婷婷网| 一本—道久久a久久精品蜜桃钙片| 一边摸一边抽搐一进一出视频| 国产伦理片在线播放av一区| 在线av久久热| 亚洲精品成人av观看孕妇| 成年人免费黄色播放视频| 十分钟在线观看高清视频www| 天天躁狠狠躁夜夜躁狠狠躁| 成人手机av| 热re99久久国产66热| 巨乳人妻的诱惑在线观看| 国产一区亚洲一区在线观看| 国产精品免费视频内射| 波野结衣二区三区在线| 婷婷丁香在线五月| 中文字幕制服av| 免费女性裸体啪啪无遮挡网站| 一级片'在线观看视频| 久久久久久免费高清国产稀缺| 夫妻性生交免费视频一级片| 两个人看的免费小视频| 极品少妇高潮喷水抽搐| 亚洲五月婷婷丁香| 超碰成人久久| 岛国毛片在线播放| 国产精品三级大全| 精品人妻1区二区| cao死你这个sao货| 午夜福利视频精品| kizo精华| 爱豆传媒免费全集在线观看| 亚洲精品国产区一区二| 亚洲av日韩在线播放| 成人18禁高潮啪啪吃奶动态图| 伊人久久大香线蕉亚洲五| 国产成人精品久久二区二区91| 亚洲国产精品一区三区| 日韩制服丝袜自拍偷拍| 90打野战视频偷拍视频| 一级,二级,三级黄色视频| 久久性视频一级片| 国产97色在线日韩免费| 91老司机精品| 亚洲精品一卡2卡三卡4卡5卡 | 免费一级毛片在线播放高清视频 | videosex国产| 女警被强在线播放| 欧美人与性动交α欧美精品济南到| 另类亚洲欧美激情| 国产欧美亚洲国产| 在线观看免费日韩欧美大片| 亚洲美女黄色视频免费看| 午夜激情久久久久久久| 国产视频首页在线观看| 黄色视频在线播放观看不卡| 日本欧美视频一区| 国产av精品麻豆| 久9热在线精品视频| 亚洲免费av在线视频| 亚洲自偷自拍图片 自拍| 乱人伦中国视频| 亚洲欧美一区二区三区黑人| 精品一区在线观看国产| 国产精品一二三区在线看| 一区二区日韩欧美中文字幕| 最近最新中文字幕大全免费视频 | 久久人人爽人人片av| 成年av动漫网址| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 国精品久久久久久国模美| 97在线人人人人妻| 国产有黄有色有爽视频| 亚洲综合色网址| 在现免费观看毛片| 精品人妻在线不人妻| 欧美成狂野欧美在线观看| 午夜日韩欧美国产| 亚洲少妇的诱惑av| 成人免费观看视频高清| 七月丁香在线播放| 国产黄色视频一区二区在线观看| 免费看不卡的av| 亚洲人成电影观看| 啦啦啦在线免费观看视频4| 久久人人97超碰香蕉20202| 国产精品av久久久久免费| 美女主播在线视频| 大片免费播放器 马上看| 波多野结衣av一区二区av| 国产精品av久久久久免费| 一级毛片 在线播放| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 国产一区二区三区av在线| 欧美黄色淫秽网站| 手机成人av网站| 免费看不卡的av| 午夜老司机福利片| 国产成人a∨麻豆精品| 免费黄频网站在线观看国产| 欧美日韩成人在线一区二区| 亚洲精品国产一区二区精华液| 不卡av一区二区三区| 亚洲国产成人一精品久久久| 久久亚洲精品不卡| 叶爱在线成人免费视频播放| 亚洲av成人精品一二三区| 亚洲精品av麻豆狂野| 中文字幕制服av| 看免费成人av毛片| 国产精品国产三级专区第一集| 久久久精品94久久精品| 精品福利观看| 人妻 亚洲 视频| 日韩人妻精品一区2区三区| 亚洲成国产人片在线观看| 欧美精品啪啪一区二区三区 | 午夜福利影视在线免费观看| 在线天堂中文资源库| 多毛熟女@视频| 亚洲少妇的诱惑av| 精品国产一区二区久久| 精品国产超薄肉色丝袜足j| 国产成人精品在线电影| 搡老乐熟女国产| 每晚都被弄得嗷嗷叫到高潮| 欧美少妇被猛烈插入视频| 亚洲专区国产一区二区| 热99久久久久精品小说推荐| 国产成人影院久久av| 2018国产大陆天天弄谢| 亚洲欧美精品综合一区二区三区| 韩国精品一区二区三区| 成人黄色视频免费在线看| 男人操女人黄网站| 啦啦啦中文免费视频观看日本| 韩国精品一区二区三区| 一级毛片黄色毛片免费观看视频| 国产成人欧美在线观看 | 国产在视频线精品| 日本欧美视频一区| 一级片免费观看大全| 国产片内射在线| 女人被躁到高潮嗷嗷叫费观| 美女脱内裤让男人舔精品视频| 成人影院久久| 在线观看免费视频网站a站| 精品国产乱码久久久久久小说| 国产亚洲一区二区精品| 国产成人影院久久av| 国产精品久久久久久精品电影小说| 大陆偷拍与自拍| 国精品久久久久久国模美| 亚洲第一青青草原| 国产日韩一区二区三区精品不卡| 亚洲精品一二三| 美女扒开内裤让男人捅视频| 国产精品香港三级国产av潘金莲 | 国产成人系列免费观看| 七月丁香在线播放| 91成人精品电影| 少妇人妻 视频| 超碰97精品在线观看| 视频在线观看一区二区三区| 国产精品免费视频内射| 欧美精品一区二区大全| 欧美人与性动交α欧美软件| 精品免费久久久久久久清纯 | 51午夜福利影视在线观看| netflix在线观看网站| 黄色一级大片看看| 高清不卡的av网站| 亚洲av美国av| 国产一区二区 视频在线| av视频免费观看在线观看| 又粗又硬又长又爽又黄的视频| 久久久国产欧美日韩av| 在线av久久热| 国产在线观看jvid| 国产男人的电影天堂91| 久久狼人影院| 男女边吃奶边做爰视频| 美女中出高潮动态图| 成年女人毛片免费观看观看9 | 亚洲人成电影观看| 一区二区三区激情视频| 香蕉国产在线看| 一边摸一边做爽爽视频免费| 在线观看免费日韩欧美大片| 婷婷色av中文字幕| 婷婷色综合大香蕉| 国产爽快片一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 精品少妇黑人巨大在线播放| 亚洲 欧美一区二区三区| 欧美人与善性xxx| 欧美大码av| 国产又色又爽无遮挡免| 免费av中文字幕在线| 日韩 欧美 亚洲 中文字幕| 久久久久国产一级毛片高清牌| 无遮挡黄片免费观看| 99久久精品国产亚洲精品| 人人妻,人人澡人人爽秒播 | 晚上一个人看的免费电影| 欧美日韩视频精品一区| 丝袜在线中文字幕| 久久精品久久久久久噜噜老黄| 在线亚洲精品国产二区图片欧美| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 操出白浆在线播放| 亚洲精品一区蜜桃| 一本色道久久久久久精品综合| 国产亚洲av片在线观看秒播厂| 少妇人妻久久综合中文| 天天躁狠狠躁夜夜躁狠狠躁| 一级片'在线观看视频| 日韩 欧美 亚洲 中文字幕| 午夜激情久久久久久久| 国产av一区二区精品久久| 免费av中文字幕在线| av视频免费观看在线观看| 精品一区在线观看国产| 深夜精品福利| 亚洲精品av麻豆狂野| 国产三级黄色录像| 国产免费一区二区三区四区乱码| 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 在线看a的网站| 亚洲欧美中文字幕日韩二区| 这个男人来自地球电影免费观看| 热re99久久精品国产66热6| 国产精品偷伦视频观看了| 久久久久久免费高清国产稀缺| 制服诱惑二区| 亚洲 国产 在线| 精品一区二区三区四区五区乱码 | 香蕉国产在线看| 国产一级毛片在线| 日本wwww免费看| 新久久久久国产一级毛片| 国产成人91sexporn| 午夜免费成人在线视频| 亚洲精品乱久久久久久| 在线观看免费视频网站a站| 亚洲综合色网址| 午夜免费观看性视频| 多毛熟女@视频| 国产不卡av网站在线观看| 制服人妻中文乱码| 看免费成人av毛片| 久久精品国产综合久久久| 免费观看av网站的网址| 九色亚洲精品在线播放| 亚洲人成电影观看| 美女视频免费永久观看网站| 国产一区二区在线观看av| 精品亚洲成国产av| 亚洲精品国产一区二区精华液| av视频免费观看在线观看| 国产日韩一区二区三区精品不卡| 啦啦啦 在线观看视频| 精品少妇内射三级| 你懂的网址亚洲精品在线观看| 欧美 亚洲 国产 日韩一| 交换朋友夫妻互换小说| 一区二区三区激情视频| 日本一区二区免费在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品欧美亚洲77777| 蜜桃在线观看..| 高潮久久久久久久久久久不卡| 国产精品九九99| 日日夜夜操网爽| videosex国产| 色婷婷久久久亚洲欧美| 最新的欧美精品一区二区| 天堂俺去俺来也www色官网| 久9热在线精品视频| av线在线观看网站| 精品久久蜜臀av无| 91麻豆精品激情在线观看国产 | 亚洲精品一卡2卡三卡4卡5卡 | 日本一区二区免费在线视频| 午夜福利视频在线观看免费| av国产精品久久久久影院| 看免费成人av毛片| 最近手机中文字幕大全| cao死你这个sao货| 婷婷色综合www| 国产成人精品在线电影| 亚洲av电影在线观看一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产成人影院久久av| netflix在线观看网站| 亚洲,欧美,日韩| 国产成人一区二区三区免费视频网站 | 一边亲一边摸免费视频| 中文字幕制服av| 在线观看免费日韩欧美大片| 欧美精品亚洲一区二区| 亚洲国产精品一区二区三区在线| 好男人视频免费观看在线| 国产欧美亚洲国产| 久久综合国产亚洲精品| 一区福利在线观看| 国产成人欧美在线观看 | 高清视频免费观看一区二区| 丰满人妻熟妇乱又伦精品不卡| 日韩视频在线欧美| 欧美黄色淫秽网站| 久久精品国产a三级三级三级| 日本色播在线视频| 国产亚洲一区二区精品| 免费看十八禁软件| 狂野欧美激情性bbbbbb| 黄频高清免费视频| 欧美成人午夜精品| 日本wwww免费看| 美女午夜性视频免费| 制服人妻中文乱码| 午夜福利在线免费观看网站| 午夜激情av网站| 免费观看a级毛片全部| 十八禁高潮呻吟视频| 国产成人啪精品午夜网站| 亚洲久久久国产精品| 99久久综合免费| 国产女主播在线喷水免费视频网站| 一区二区三区激情视频| 国产成人91sexporn| 国产主播在线观看一区二区 | 香蕉丝袜av| 狂野欧美激情性bbbbbb| 操出白浆在线播放| 国产亚洲精品第一综合不卡| 丰满人妻熟妇乱又伦精品不卡| 成年人免费黄色播放视频| 欧美人与性动交α欧美软件| 欧美亚洲 丝袜 人妻 在线| www.自偷自拍.com| 国产精品人妻久久久影院| 另类亚洲欧美激情| 亚洲成人手机| 久久久亚洲精品成人影院| 久久久久久人人人人人| 亚洲国产欧美在线一区| 人妻一区二区av| 天天躁日日躁夜夜躁夜夜| 久久久亚洲精品成人影院| 三上悠亚av全集在线观看| 亚洲国产欧美在线一区| 亚洲精品美女久久久久99蜜臀 | 免费看十八禁软件| 丝瓜视频免费看黄片| www.精华液| 黄色一级大片看看| 国产精品久久久久久精品电影小说| 男人操女人黄网站| 国产又爽黄色视频| 久久99热这里只频精品6学生| 少妇裸体淫交视频免费看高清 | 韩国高清视频一区二区三区| 两人在一起打扑克的视频| 国产又色又爽无遮挡免| av国产久精品久网站免费入址| 国产精品 国内视频| 老汉色∧v一级毛片| 亚洲色图 男人天堂 中文字幕| 亚洲国产中文字幕在线视频| 成年人午夜在线观看视频| 亚洲av日韩在线播放| 2021少妇久久久久久久久久久| 男女午夜视频在线观看| 免费观看人在逋| 国产免费一区二区三区四区乱码| 成人午夜精彩视频在线观看| 国产男女内射视频| 久热这里只有精品99| 伦理电影免费视频| 亚洲国产欧美网| 国产真人三级小视频在线观看| 午夜激情久久久久久久| 人妻 亚洲 视频| 成人影院久久| 欧美人与善性xxx| 亚洲色图综合在线观看| 国产欧美日韩精品亚洲av| 亚洲欧美日韩高清在线视频 | 在线观看人妻少妇| 国产麻豆69| 男女午夜视频在线观看| 国产成人欧美在线观看 | 性高湖久久久久久久久免费观看| 免费一级毛片在线播放高清视频 | 大片电影免费在线观看免费| h视频一区二区三区| 欧美少妇被猛烈插入视频| 91麻豆精品激情在线观看国产 | 50天的宝宝边吃奶边哭怎么回事| 精品免费久久久久久久清纯 | 一级毛片 在线播放| 人妻一区二区av| 热99国产精品久久久久久7| 婷婷色综合大香蕉| 久久亚洲精品不卡| 色播在线永久视频| 大香蕉久久成人网| 免费少妇av软件| 日日摸夜夜添夜夜爱| 亚洲黑人精品在线| 久久人人97超碰香蕉20202| 伦理电影免费视频| 你懂的网址亚洲精品在线观看| 国产精品久久久久久精品古装| 国产精品人妻久久久影院| 啦啦啦 在线观看视频| 午夜激情av网站| 亚洲伊人久久精品综合| 1024香蕉在线观看| 高潮久久久久久久久久久不卡| 女人久久www免费人成看片| 色94色欧美一区二区| 亚洲精品久久午夜乱码| 18禁国产床啪视频网站| 久久国产亚洲av麻豆专区| 国产伦人伦偷精品视频| 色婷婷久久久亚洲欧美| 亚洲欧美精品自产自拍| 国产老妇伦熟女老妇高清| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的| 亚洲熟女精品中文字幕| 国产在线一区二区三区精| 18禁裸乳无遮挡动漫免费视频| 男女无遮挡免费网站观看| 搡老乐熟女国产| 欧美激情极品国产一区二区三区| 一级毛片我不卡| 亚洲黑人精品在线| 性色av一级| 国产精品秋霞免费鲁丝片| 在线观看一区二区三区激情| svipshipincom国产片| 91精品三级在线观看| 亚洲欧美一区二区三区久久| 高清黄色对白视频在线免费看| 欧美成人午夜精品| 国产极品粉嫩免费观看在线| 欧美av亚洲av综合av国产av| 国产精品三级大全| 高清视频免费观看一区二区| 国产1区2区3区精品| 男女下面插进去视频免费观看| 成年av动漫网址| 久久综合国产亚洲精品| 欧美精品人与动牲交sv欧美|