• <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

    日本vs欧美在线观看视频| av国产久精品久网站免费入址| 久久精品亚洲熟妇少妇任你| 高清在线视频一区二区三区| 深夜精品福利| 国产精品偷伦视频观看了| 日韩中文字幕视频在线看片| 国产精品嫩草影院av在线观看| 美女视频免费永久观看网站| 91精品三级在线观看| 丝袜人妻中文字幕| 多毛熟女@视频| 日韩欧美精品免费久久| 久久久久久久国产电影| 欧美激情 高清一区二区三区| 男女下面插进去视频免费观看| 一级毛片 在线播放| 亚洲成人手机| 新久久久久国产一级毛片| 日本黄色日本黄色录像| 欧美精品一区二区免费开放| 校园人妻丝袜中文字幕| 欧美成人精品欧美一级黄| 国产不卡av网站在线观看| 亚洲精品aⅴ在线观看| 国产精品人妻久久久影院| av电影中文网址| 久久久久久久大尺度免费视频| 一区二区日韩欧美中文字幕| 国产成人精品无人区| 少妇的丰满在线观看| 天堂8中文在线网| 丁香六月欧美| 成人18禁高潮啪啪吃奶动态图| 国产精品av久久久久免费| 国产精品偷伦视频观看了| 亚洲精品乱久久久久久| 男女高潮啪啪啪动态图| 日本av手机在线免费观看| 欧美日韩国产mv在线观看视频| 亚洲欧美中文字幕日韩二区| 丰满饥渴人妻一区二区三| 国产成人系列免费观看| 亚洲精品一二三| 久久人妻熟女aⅴ| 精品一区二区三区四区五区乱码 | 在线天堂最新版资源| 一级毛片黄色毛片免费观看视频| 欧美日韩亚洲高清精品| 人人妻人人添人人爽欧美一区卜| av在线观看视频网站免费| 亚洲欧美成人精品一区二区| 国产av精品麻豆| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一区蜜桃| 亚洲伊人色综图| 精品亚洲乱码少妇综合久久| 亚洲成人国产一区在线观看 | 免费少妇av软件| 国产精品蜜桃在线观看| 老汉色∧v一级毛片| 国产成人免费观看mmmm| 精品一区二区免费观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲美女视频黄频| 精品久久久久久电影网| 日本av免费视频播放| 国产成人免费观看mmmm| 久久性视频一级片| 欧美日韩视频精品一区| av福利片在线| 免费不卡黄色视频| 各种免费的搞黄视频| 777米奇影视久久| 亚洲精品日本国产第一区| 亚洲熟女毛片儿| 国产一区有黄有色的免费视频| 男女床上黄色一级片免费看| 久久久欧美国产精品| 99久久人妻综合| 一区二区三区激情视频| 欧美精品亚洲一区二区| 精品亚洲成国产av| 国产xxxxx性猛交| 看免费成人av毛片| 欧美xxⅹ黑人| 日本av免费视频播放| 菩萨蛮人人尽说江南好唐韦庄| 欧美少妇被猛烈插入视频| 国产视频首页在线观看| 国产精品 欧美亚洲| 毛片一级片免费看久久久久| 自线自在国产av| 国产一区二区三区综合在线观看| 狠狠婷婷综合久久久久久88av| 亚洲激情五月婷婷啪啪| 成人国产av品久久久| 国产欧美日韩综合在线一区二区| 国产精品秋霞免费鲁丝片| 国产激情久久老熟女| 黄片播放在线免费| 亚洲欧美一区二区三区久久| 亚洲熟女毛片儿| 丝袜喷水一区| 赤兔流量卡办理| 香蕉丝袜av| 观看av在线不卡| 国产高清国产精品国产三级| 中文字幕制服av| 99精国产麻豆久久婷婷| 多毛熟女@视频| 天堂俺去俺来也www色官网| 国产探花极品一区二区| 亚洲成国产人片在线观看| 热99国产精品久久久久久7| 欧美黑人精品巨大| 这个男人来自地球电影免费观看 | 成人免费观看视频高清| 国产男人的电影天堂91| 一区二区三区激情视频| 精品国产乱码久久久久久男人| 精品午夜福利在线看| 国产精品一区二区精品视频观看| 精品人妻熟女毛片av久久网站| 中文字幕人妻熟女乱码| 在现免费观看毛片| 精品少妇黑人巨大在线播放| av福利片在线| 下体分泌物呈黄色| 亚洲av日韩精品久久久久久密 | 国产片内射在线| 麻豆乱淫一区二区| 亚洲人成77777在线视频| 日韩大片免费观看网站| 在线免费观看不下载黄p国产| 免费在线观看完整版高清| 免费在线观看黄色视频的| 精品国产国语对白av| 国产xxxxx性猛交| 18在线观看网站| 观看美女的网站| av线在线观看网站| 在线看a的网站| 伦理电影免费视频| 国产在线一区二区三区精| 天堂中文最新版在线下载| 国产精品三级大全| 国产精品久久久久久久久免| 欧美激情高清一区二区三区 | 捣出白浆h1v1| 自拍欧美九色日韩亚洲蝌蚪91| 97精品久久久久久久久久精品| 久久久久精品人妻al黑| 精品久久蜜臀av无| 丝袜在线中文字幕| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 国产无遮挡羞羞视频在线观看| 亚洲av日韩精品久久久久久密 | 观看美女的网站| 女人高潮潮喷娇喘18禁视频| 两性夫妻黄色片| 色播在线永久视频| 国产成人午夜福利电影在线观看| 最近最新中文字幕免费大全7| 男人爽女人下面视频在线观看| 久久av网站| 亚洲av男天堂| 国产一级毛片在线| 亚洲一码二码三码区别大吗| 丁香六月欧美| 免费看不卡的av| 欧美人与性动交α欧美精品济南到| 欧美乱码精品一区二区三区| 亚洲第一青青草原| 各种免费的搞黄视频| 青春草国产在线视频| 久久女婷五月综合色啪小说| 两个人看的免费小视频| 日本欧美视频一区| 一区二区av电影网| 国产亚洲av高清不卡| 国产一区二区激情短视频 | 精品一区二区三区四区五区乱码 | 18禁观看日本| 久久99一区二区三区| 9色porny在线观看| 欧美在线黄色| 午夜老司机福利片| 日韩中文字幕欧美一区二区 | 桃花免费在线播放| 老汉色∧v一级毛片| 黄频高清免费视频| 国产免费视频播放在线视频| 秋霞伦理黄片| 99热网站在线观看| 亚洲三区欧美一区| 人成视频在线观看免费观看| 大片免费播放器 马上看| 美女中出高潮动态图| 欧美人与性动交α欧美精品济南到| 中文字幕精品免费在线观看视频| www.av在线官网国产| 国产一区二区 视频在线| av在线app专区| 日韩av在线免费看完整版不卡| 在线观看免费高清a一片| 国产又爽黄色视频| 日本wwww免费看| 亚洲欧美中文字幕日韩二区| 国产熟女欧美一区二区| 久久97久久精品| 免费少妇av软件| 别揉我奶头~嗯~啊~动态视频 | 丰满少妇做爰视频| 欧美激情极品国产一区二区三区| 亚洲欧洲日产国产| 男人舔女人的私密视频| 国产黄色免费在线视频| 中文乱码字字幕精品一区二区三区| 国产探花极品一区二区| av天堂久久9| videos熟女内射| 丝袜人妻中文字幕| 男人爽女人下面视频在线观看| 国产在线免费精品| 国产亚洲欧美精品永久| 精品人妻一区二区三区麻豆| 99香蕉大伊视频| av女优亚洲男人天堂| 女性生殖器流出的白浆| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| av福利片在线| 久久久久久免费高清国产稀缺| 叶爱在线成人免费视频播放| av国产精品久久久久影院| 一区福利在线观看| 少妇人妻 视频| 亚洲五月色婷婷综合| 波野结衣二区三区在线| 免费观看性生交大片5| 男人添女人高潮全过程视频| 久久久久精品久久久久真实原创| 天堂俺去俺来也www色官网| 51午夜福利影视在线观看| 精品一区二区三区av网在线观看 | 七月丁香在线播放| 亚洲专区中文字幕在线 | 亚洲国产成人一精品久久久| 国产精品国产三级专区第一集| 女人爽到高潮嗷嗷叫在线视频| 黄片小视频在线播放| 国产97色在线日韩免费| 日韩制服丝袜自拍偷拍| 又大又爽又粗| 国产精品嫩草影院av在线观看| av有码第一页| 久久久久久久精品精品| 免费看不卡的av| 久久久久国产一级毛片高清牌| 亚洲国产成人一精品久久久| 视频在线观看一区二区三区| a级毛片黄视频| 免费高清在线观看日韩| 中文字幕av电影在线播放| 国产不卡av网站在线观看| 可以免费在线观看a视频的电影网站 | 一二三四在线观看免费中文在| videosex国产| 国产av国产精品国产| 两个人免费观看高清视频| 色婷婷av一区二区三区视频| 在线亚洲精品国产二区图片欧美| 亚洲一区二区三区欧美精品| 高清欧美精品videossex| 欧美人与善性xxx| 午夜久久久在线观看| 亚洲国产精品国产精品| 日韩伦理黄色片| 亚洲精品成人av观看孕妇| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美一区二区三区久久| 色94色欧美一区二区| 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 成人三级做爰电影| 久久国产精品男人的天堂亚洲| 成人国产麻豆网| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| av电影中文网址| 久久鲁丝午夜福利片| 久久久久精品国产欧美久久久 | 国产老妇伦熟女老妇高清| 久久精品人人爽人人爽视色| 国产av国产精品国产| 黄片小视频在线播放| 国产精品三级大全| 精品人妻在线不人妻| 中国三级夫妇交换| 无遮挡黄片免费观看| 国产麻豆69| 国产无遮挡羞羞视频在线观看| 日韩av在线免费看完整版不卡| 免费日韩欧美在线观看| 你懂的网址亚洲精品在线观看| 91精品伊人久久大香线蕉| 91精品三级在线观看| 久久久久人妻精品一区果冻| 国产精品久久久人人做人人爽| 色综合欧美亚洲国产小说| 国产乱人偷精品视频| av免费观看日本| 99热全是精品| 只有这里有精品99| 最近手机中文字幕大全| 国产精品香港三级国产av潘金莲 | 九色亚洲精品在线播放| 亚洲美女搞黄在线观看| a级毛片在线看网站| 综合色丁香网| 亚洲图色成人| a级毛片在线看网站| 一区二区三区乱码不卡18| 亚洲视频免费观看视频| 少妇的丰满在线观看| 在线亚洲精品国产二区图片欧美| 亚洲精品aⅴ在线观看| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 国产熟女欧美一区二区| svipshipincom国产片| 蜜桃国产av成人99| 男女床上黄色一级片免费看| 日韩不卡一区二区三区视频在线| 午夜福利影视在线免费观看| 亚洲第一av免费看| 欧美精品人与动牲交sv欧美| 亚洲,一卡二卡三卡| 免费黄频网站在线观看国产| 免费观看av网站的网址| 国产高清国产精品国产三级| av线在线观看网站| 高清黄色对白视频在线免费看| 国产日韩一区二区三区精品不卡| 亚洲国产成人一精品久久久| 国产免费福利视频在线观看| 欧美久久黑人一区二区| 一区福利在线观看| 免费久久久久久久精品成人欧美视频| 午夜免费观看性视频| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 色婷婷av一区二区三区视频| tube8黄色片| 国产成人欧美在线观看 | 熟女av电影| 青草久久国产| 国产亚洲av高清不卡| 亚洲精品视频女| 日韩 亚洲 欧美在线| 亚洲av日韩在线播放| 丁香六月欧美| 黑丝袜美女国产一区| 99久久综合免费| 国产高清不卡午夜福利| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线| 在现免费观看毛片| 日韩av不卡免费在线播放| 一边摸一边抽搐一进一出视频| 在线看a的网站| 国产片特级美女逼逼视频| 满18在线观看网站| 国产精品免费视频内射| 青草久久国产| 狠狠精品人妻久久久久久综合| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲美女搞黄在线观看| 又大又黄又爽视频免费| 99久久人妻综合| 免费不卡黄色视频| 国产xxxxx性猛交| 欧美日韩亚洲国产一区二区在线观看 | 婷婷色av中文字幕| 69精品国产乱码久久久| 老熟女久久久| 精品午夜福利在线看| 国产精品香港三级国产av潘金莲 | 精品少妇一区二区三区视频日本电影 | 伦理电影免费视频| 国产精品二区激情视频| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 亚洲精品久久午夜乱码| 国产精品.久久久| 好男人视频免费观看在线| 美女脱内裤让男人舔精品视频| 又黄又粗又硬又大视频| 亚洲av欧美aⅴ国产| 免费黄色在线免费观看| 亚洲天堂av无毛| 男人舔女人的私密视频| 一本色道久久久久久精品综合| 亚洲成人免费av在线播放| 亚洲成国产人片在线观看| 欧美97在线视频| 熟女少妇亚洲综合色aaa.| 亚洲精品久久成人aⅴ小说| 国产男女内射视频| 午夜福利,免费看| 亚洲色图综合在线观看| 欧美最新免费一区二区三区| 在现免费观看毛片| 18在线观看网站| 亚洲av综合色区一区| 99热国产这里只有精品6| 免费黄网站久久成人精品| 性少妇av在线| 巨乳人妻的诱惑在线观看| 国产成人啪精品午夜网站| 久久久欧美国产精品| 国产成人av激情在线播放| www.av在线官网国产| av国产久精品久网站免费入址| 男女之事视频高清在线观看 | 亚洲精品国产av蜜桃| 少妇人妻 视频| 午夜福利视频在线观看免费| 国产亚洲午夜精品一区二区久久| 久久av网站| 丁香六月欧美| 男女之事视频高清在线观看 | 狠狠精品人妻久久久久久综合| 国产成人91sexporn| 在线观看人妻少妇| 免费久久久久久久精品成人欧美视频| 亚洲男人天堂网一区| 欧美国产精品一级二级三级| 国产色婷婷99| 亚洲精品视频女| 久久久精品国产亚洲av高清涩受| 丰满饥渴人妻一区二区三| 欧美精品一区二区大全| 校园人妻丝袜中文字幕| 婷婷色综合www| 亚洲av在线观看美女高潮| 99香蕉大伊视频| 国产无遮挡羞羞视频在线观看| 国产精品偷伦视频观看了| 免费看av在线观看网站| 汤姆久久久久久久影院中文字幕| 熟女av电影| 国产精品秋霞免费鲁丝片| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 最近的中文字幕免费完整| 中文天堂在线官网| 亚洲综合精品二区| 久久久久久久久久久免费av| 综合色丁香网| 欧美激情 高清一区二区三区| 欧美国产精品一级二级三级| 国产又爽黄色视频| 一本—道久久a久久精品蜜桃钙片| 男人舔女人的私密视频| av福利片在线| 国产精品久久久久久久久免| 成人黄色视频免费在线看| 亚洲成人一二三区av| 国产一区二区三区av在线| 天堂俺去俺来也www色官网| 免费高清在线观看视频在线观看| 伊人久久大香线蕉亚洲五| xxx大片免费视频| 国产精品久久久人人做人人爽| 久久久久久久精品精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品成人av观看孕妇| 欧美日韩福利视频一区二区| 亚洲成人免费av在线播放| 伦理电影免费视频| videosex国产| 丝袜脚勾引网站| 在线观看免费日韩欧美大片| 亚洲少妇的诱惑av| 欧美另类一区| 亚洲av日韩精品久久久久久密 | 99九九在线精品视频| 亚洲av欧美aⅴ国产| 婷婷成人精品国产| 国产成人免费观看mmmm| kizo精华| 卡戴珊不雅视频在线播放| 久久99一区二区三区| 另类精品久久| 各种免费的搞黄视频| 熟女av电影| 99久久99久久久精品蜜桃| 国产精品久久久久久人妻精品电影 | 欧美日韩亚洲综合一区二区三区_| 丝袜脚勾引网站| 最黄视频免费看| 国产免费又黄又爽又色| 日韩成人av中文字幕在线观看| 叶爱在线成人免费视频播放| 亚洲欧美激情在线| 七月丁香在线播放| 国产伦理片在线播放av一区| 久久精品国产综合久久久| 久久天躁狠狠躁夜夜2o2o | 黑丝袜美女国产一区| 久热爱精品视频在线9| 日本av免费视频播放| 丝袜喷水一区| 我的亚洲天堂| 成年美女黄网站色视频大全免费| 免费黄色在线免费观看| 极品人妻少妇av视频| 秋霞伦理黄片| 91精品伊人久久大香线蕉| 欧美精品人与动牲交sv欧美| 成年av动漫网址| 一级片'在线观看视频| 国产一区二区三区综合在线观看| 色婷婷久久久亚洲欧美| 美女主播在线视频| 久久婷婷青草| 男男h啪啪无遮挡| 精品福利永久在线观看| 国产成人欧美| 大码成人一级视频| 成年人免费黄色播放视频| 国产老妇伦熟女老妇高清| 亚洲一卡2卡3卡4卡5卡精品中文| 国产淫语在线视频| av国产精品久久久久影院| 欧美在线黄色| 久久 成人 亚洲| 母亲3免费完整高清在线观看| 18禁观看日本| av一本久久久久| 黄色一级大片看看| 老司机在亚洲福利影院| 欧美日韩福利视频一区二区| 久久精品aⅴ一区二区三区四区| 人人妻人人添人人爽欧美一区卜| 国产免费福利视频在线观看| 欧美 日韩 精品 国产| 中文字幕色久视频| 日韩电影二区| 激情五月婷婷亚洲| 亚洲成色77777| 麻豆乱淫一区二区| 一区二区av电影网| 天天影视国产精品| 亚洲欧洲精品一区二区精品久久久 | 看非洲黑人一级黄片| 精品一品国产午夜福利视频| 嫩草影院入口| 麻豆av在线久日| 精品亚洲成a人片在线观看| 人人澡人人妻人| 免费观看人在逋| 婷婷成人精品国产| 午夜久久久在线观看| 亚洲av电影在线进入| 国产精品国产三级国产专区5o| 成人国产av品久久久| 男女高潮啪啪啪动态图| 99精国产麻豆久久婷婷| 国产精品久久久久久精品古装| 国产97色在线日韩免费| 亚洲在久久综合| 五月天丁香电影| 少妇人妻 视频| 黄片小视频在线播放| 久久精品亚洲熟妇少妇任你| 国产欧美日韩综合在线一区二区| 久久毛片免费看一区二区三区| 午夜福利影视在线免费观看| 免费观看a级毛片全部| 午夜免费观看性视频| av国产久精品久网站免费入址| 成人手机av| 国产淫语在线视频| netflix在线观看网站| 黄片无遮挡物在线观看| 国产免费又黄又爽又色| av卡一久久| 国产精品国产三级专区第一集| 亚洲欧美色中文字幕在线| 日日摸夜夜添夜夜爱| 日韩人妻精品一区2区三区| 涩涩av久久男人的天堂| 久久青草综合色| 视频在线观看一区二区三区| 亚洲国产毛片av蜜桃av| av国产精品久久久久影院| 亚洲精品日韩在线中文字幕| 国产一区有黄有色的免费视频| 啦啦啦视频在线资源免费观看| 亚洲精品日韩在线中文字幕| 午夜影院在线不卡| 黄色毛片三级朝国网站| 超色免费av| 精品视频人人做人人爽| 丰满少妇做爰视频| 69精品国产乱码久久久| 热99久久久久精品小说推荐| 日本欧美国产在线视频| 日本wwww免费看| 精品人妻熟女毛片av久久网站| 一级a爱视频在线免费观看|