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

    Design and research of magnetically levitated testbed with composite superconductor bearing for micro thrust measurement

    2021-10-31 08:15:02FawziDERKAOUIZhaoxinLIU劉朝鑫WenjiangYANG楊文將YuQIN秦瑀KunlongWU吳坤隆PengZHAO趙鵬JuzhuangYAN閆炬壯JunxueREN任軍學(xué)andHaibinTANG湯海濱
    Plasma Science and Technology 2021年10期
    關(guān)鍵詞:趙鵬海濱

    Fawzi DERKAOUI, Zhaoxin LIU (劉朝鑫), Wenjiang YANG (楊文將),Yu QIN (秦瑀), Kunlong WU (吳坤隆), Peng ZHAO (趙鵬), Juzhuang YAN(閆炬壯), Junxue REN (任軍學(xué)) and Haibin TANG (湯海濱)

    1 School of Astronautics, Beihang University, Beijing 100191, People’s Republic of China

    2 School of Power Engineering,Liuzhou Railway Vocational Technical College,Liuzhou 545616,People’s Republic of China

    3 School of Space and Environment, Beihang University, Beijing 100191, People’s Republic of China

    Abstract A high temperature superconducting (HTS) magnetically levitated testbed has been developed for the steady thrust measurement of miniature ion electrospray thruster.The structure of the testbed mainly consists of an HTS composite bearing, a magnetic shielding plate, an active electromagnetic brake and a laser displacement sensor.The steady thrust is described as a function of the equilibrium angle displacement of the floating frame.Furthermore, the mechanical behaviors of HTS composite bearing were studied via finite element simulation and experiments, which include the load capacity, levitation stiffness and background noise.The results show that the thrust testbed can keep in low noise and have a load capacity up to 4 kg.According to the ignition testing of the electrospray thruster, the thrust force of 25.2μN(yùn) was measured by the testbed,which is close to the design value of miniature ion electrospray thruster.

    Keywords: micro-force measurement, maglev testbed, composite HTS bearing

    1.Introduction

    With the rapid development of micro-thrusters, more potential space applications such as micro-satellites formations are activated.Among them, the task of gravitational wave detection requires micro-thrusters to compensate for the space interference force by providing a very precise thrust,so the thrust performance has to be evaluated very accurately to prove the design of the precise micro-thrusters[1–3].In the past time,a lot of researches on micro-thrust measurement have been carried out and have made great progress.Up to now, the thrust in the sub-micro newton range can be measured by the mechanical pivot stand[4].

    Hathaway from Canada designed a sensitive vacuum microthrust test stand,which introduced a pair of jewel-pivot bearings to measure the thrust of 1–100μN(yùn) and the resolution is 20 nN.This type of pivot bearing can sustain significant load with extremely low rotation friction[4].Chakrabortyet alhas designed a thrust test stand which can measure the thrust up to 100μN(yùn),and the resolution is 10 nN.The device used the homodyne thrust measurement technique to effectively isolate the interference of facility vibration on thrust reading,and can be used for different types of micro thrusters [5].Soni and Roy form university of Florida designed a nano-Newton resolution thrust stand [6], and the logarithmic decrement method was proposed to calibrate the test stand.The uncertainty of calibration process can be limited to 5%by reducing the number of error sources.The signal analysis technology was used to analyze the noise source and its corresponding frequency band, and the noise correction sensor was introduced to further improve the resolution of the test stand.However, there are also some problems such as inability to effectively isolate noise, internal force interference and low load capacity in the mechanical pivot.To solve these problems, a micro-thrust testbed based on magnetic levitation with high temperature superconducting(HTS)composite bearings is shown in this paper.

    The phenomenon of HTS magnetic levitation was first discovered in 1987.Nowadays,HTS magnetic bearing is one of the most prospective applications of bulk superconductors in stable levitation with considerable stiffness and self-stability due to its unique flux pinning property.HTS magnetic bearing has great potential application of emerging technologies such as magnetic levitation transportation [7] and flywheel energy storage systems (FESS) [8–11].Boeing has developed an HTS FESS with a speed of up to 22500 rpm and 5 kWh, which uses an axial magnetic bearing structure.The radial bearing structures were developed by Germany’s Nexans Superconducting Company and Siemens [12].Furthermore, the German ATZ company studied the levitation and stiffness characteristics of HTS bearings [13, 14].

    To improve the levitation characteristics of micro-thrust testbed, an axial-radial composite HTS bearing was introduced in this paper, which can realize the self-stabilizing levitation of testbed without complex external control system and rotation loss influence.

    2.Measuring principle and thruster

    The measurement principle of the testbed based on the composite HTS bearing is similar to that of the traditional mechanical pivot stand.The thrust is calculated by the measured displacement or deformation of the floating frame.The oscillation of floating frame follows second-order equation:

    where,θis the rotation angle,Iis the floating frame’s inertia moment,cis the damping coefficient of frame,kthe is restoring moment coefficient,RTthe is equivalent thrust arm length andFTis the thrust force.The schematic diagram of the thrust measurement is shown in figure 1.

    According to vibration mechanics, if ignition frequency of the thruster is largely smaller than the natural frequency of floating frame, then the formula (1) can be simplified as:

    This means that the thrust to be measured is mainly calculated from the restoring moment coefficientkand the amount of change in angle.The floating frame deviates from initial position under the action of thrust until the restoring force of the floating frame is equal to the small thrust and the frame oscillates near the new equilibrium point.The restoring moment coefficientkin formula(2)can be obtained through calibration,or obtained from the free oscillation curve of floating frame.The natural angular frequencyω0of the floating frame satisfies:

    According to formula(3),the natural periodTof the floating frame can be obtained as:

    The restoring moment coefficientkis thus obtained as:

    The rotation angle of floating frame can be obtained by the laser displacement sensor.The natural period is calculated according to the oscillation curve of floating frame.Finally, the thrust force can be calculated according to the formulas (2)–(5).

    In this micro-thrust measurement, the miniature ion electrospray thruster developed by Beihang University is evaluated.The effective emission current reach 200μA and the thrust is about 20μN(yùn) under the accelerating voltage of 2500 V.The thruster is shown in figure 2.

    3.Design of thrust testbed

    The structure of HTS magnetically levitated testbed is shown in figure 3.The whole propulsion system, including thrusters and power units,is installed on the floating frame to avoid the interference between power units and the HTS composite bearing to suppress noise of the testbed, which is an improvement on the previous version [15].The thruster is ignited when the floating frame is in stable levitation,and the rotation movement of floating frame is measured and analyzed to calculate the thrust.

    3.1.Testbed structure

    The challenge of using the deflection angle of the floating platform to calculate the thrust is to reduce the hysteresis in the floating frame.The sources of hysteresis mainly include restoring moments and different resistances such as rotational friction.Feasible measures to reduce hysteresis include the more reasonable structure design, precise counterweight and the horizontal adjustment of floating frame [16].

    The restoring moment is mainly composed of gravity restoring moment and electromagnetic torque.The error caused by the fabrication and installation makes the floating frame have an inherent gravitational potential well.The nonuniform distribution of surrounding magnetic field produces an electromagnetic torque.The restoring moment coefficientkdetermined by the inherent restoring moment can be calculated by the passive calibration method.

    In addition,the performance evaluation of thruster for the gravitational wave detection mission requires that the thrust testbed has very low noise.The LISA Pathfinder spacecraft jointly developed by ESA and NASA requires the background noise of thrust testbed is less than 0.1 μN(yùn)/Hz0.5within 1 mHz–4 Hz [1].

    The measurement testbed designed in this work is mainly optimized on the structural design.The testbed adopts innovatively axial-radial composite HTS bearing to improve the load capacity and reduce noise of the testbed.The introduction of composite bearing improves the radial stiffness of rotation floating frame and can better absorb radial noise.More importantly, compared with electromagnetic levitation,HTS bearing has the outstanding advantage of self-stabilizing levitation due to its unique magnetic flux pinning characteristics, which can eliminate the complex closed-loop control system and effectively avoid the introduction of control noise.HTS bearing also has better thermal noise characteristics.

    The magnetic shielding plate in figure 3 is placed on the upper surface of floating frame,which can greatly weaken the vertical magnetic field strength, and effectively suppress the electromagnetic interference on the thrust measurement.

    Low-friction levitation means that oscillation of floating frame will last for a long time, so an active electromagnetic brake is used to slow down the oscillation.The electromagnetic brake is usually composed of a controllable magnetic field source and conductive plate, as shown in figure 3.The principle of electromagnetic brake is that the electromagnet produces magnetic field and generates eddy current in conductive plate, and then eddy current in the magnetic field will produce Lorentz force to resist oscillation.

    3.2.Simulation and design

    The design parameters of HTS composite bearing are calculated by electromagnetic finite simulations, which also study the magnetic flux density and leakage of the magnetic bearing at the same time.Table 1 shows the specific parameters of composite bearing structure, and the structure model of composite bearing is shown in figure 4.

    Table 1.Simulation structure parameters of composite bearing.

    Table 2.The data of five repeated experiments.

    Both of the axial and radial bearings are composed of two permanent magnet (PM) and an iron ring, the iron ring is introduced to increase magnetic field density at iron positions.In order to obtain larger levitation force and make the bearing capacity stronger, the thicknesses of PMs are designed as large as possible in the limited space.Furthermore, the authors added an axial PM container to make the diameters of axial PMs match the width of HTS bulk,which can reduce the magnetic flux leakage caused by the coupling of axial and radial PMs magnetic field according to the simulation results.

    Figure 5 shows the simulation results of the magnetic flux density.An axial two-dimensional section line along the surface of superconducting stator is drawn to calculate the magnetic flux density, and it is shown that the maximum magnetic flux density is 0.16 T.The two-dimensional crosssections at different gap heights along the radial direction are also drawn and the authors obtain that when the gap height is between 3 mm and 7 mm, the corresponding magnetic flux density is from 0.3 T to 0.1 T.

    Figure 1.Schematic diagram of micro thrust measurement.

    Figure 2.Ionic liquid electrospray thruster.

    Figure 3.Magnetically levitated testbed with HTS composite bearing.

    Figure 4.Simulation structure model of composite bearing.

    Figure 5.Simulation results of composite bearing.

    Figure 6.Composite bearing (left: HTS stator; right: composite PM rotor).

    According to parameters analysis obtained from the simulation, the composite bearing is designed and manufactured.The composite bearing mainly includes an arcshaped superconducting bulk,a copper fixing frame,axial and radial PMs and axial PMs container.The physical diagram is shown in figure 6.

    Me: How s the salmon2?Server: Fantastic!Me: Does it come with rice?Server: Absolutely!Would a good and a yes have been sufficient? Undeniably!At Starbucks, the smallest coffee you can order is a Tall

    3.3.Motion detector

    The measurement method used in this work is to use an Archimedes measuring block to convert the angular displacement change of the floating frame into the distance change between the measuring block surface and the laserdisplacement sensor, as shown in figure 7.The method can not only accurately measure the rotation angleθof the floating frame, but also obtain the angular velocity and angular acceleration by the displacement differentiation,without causing any interference to the motion state of floating frame.More importantly,the sensor can also identify the offset of the floating frame in the radial direction.The accuracy of the laser measuring system based on the Archimedes measuring block is 0.03°.

    Deflection angle of the floating frame has a linear relationship with the displacement, the measurement equation is designed:

    It can be seen that the surface of measuring block should be designed as Archimedes spiral, as shown in the figure 8.

    The minimum radius of Archimedes measuring block is 15 mm, and the maximum is 21 mm.The radius increases or decreases linearly with the change of the angle.The angle of Archimedes measuring block changes by 180°, while the linear displacement recognized by the laser displacement sensor changes linearly by 6 mm at the same time, as shown in figure 8(a).

    The resolution of laser displacement sensor is 1 μm,which means that the resolution of angular displacement is

    4.Experiments and analysis

    4.1.Levitation force and stiffness testing

    The levitation force and stiffness determine the loading characteristics of composite HTS bearing.The load capacity can be evaluated by the levitation force-displacement curve.Stiffness is used to study the loading stability of HTS bearing and shows the ability that effectively resists the interference to the floating frame.The stiffness can be obtained by differentiating the levitation force-displacement curve,

    In general, larger axial bearing capacity means smaller radial stiffness, which requires the authors to make a balance between the bearing capacity and stability, in order to make the testbed have better radial stiffness and low rotating noise.

    In the experiment, the HTS bulk was cooled down to 77 K at a distance of 8 mm below the surface of PM rotor to capture magnetic flux and get a large pinning force to obtain radial stiffness, which is also named field cooling (FC) for bulk superconductors.The PM rotor is moved vertically downward at 1 mm s?1until the gap between the axial PM and the HTS bulk is 1 mm and there is no pause for PM rotor moving vertically at 1 mm s?1upward to its initial position.Finally, the levitation force displacement curve of PM rotor can be obtained.Figure 10 shows the experimental process.

    4.2.Micro-thrust measurement

    The thruster ignition test was carried out in the vacuum condition including a vacuum sub-system and cooling structures.The vacuum sub-system consists of a non-magnetic vacuum chamber and vacuum pump.The minimum vacuum is lower than 10?4Pa.

    In order to achieve rapid and uniform cooling for the superconducting bulks, the double-channel conduction-cooling structure shown in figure 11 is designed.The structure includes a cold head of which the minimum temperature can reach 4.2 K, copper heat conduction structure and alumina ceramic cover plate.The heat of the superconductor is transferred from the upper and lower directions to the cold head.This cooling structure greatly improves the cooling speed and effectively guarantees the temperature uniformity of superconducting bulks.

    The authors assessed the background noise level of floating frame based on composite superconducting bearing before passive calibration and ignition measurement.Figure 12 shows the oscillation state of the floating frame and power spectrum density (PSD) of the background noise.

    Figure 12(a) shows that the floating frame has low zero drift and stable oscillation state, suggesting that the floating frame has a stable oscillation angle around the vertical axis with no radial offset.The data is further processed to obtain the PSD of floating frame, so as to obtain the noise level of different frequency bands, as shown in figure 12(b).Compared with a single axial bearing, the composite bearing designed in this article shows better low-frequency noise characteristics by improving the radial stiffness, as shown in figure 13.

    Figure 7.The laser measurement system.

    Figure 8.(a) Archimedes spiral diagram, and (b) Archimedes block physical diagram.

    Figure 9.Two-dimensional measuring device.

    Figure 10.Experimental process.

    Figure 11.Cooling structures and vacuum sub-system.

    Figure 12.(a) Oscillation state of floating frame, and (b) power spectrum density of the background noise.

    Figure 13.(a) Oscillation state of floating frame based on different bearing, and (b) PSD of single axial superconducting bearing.

    Figure 14.Free oscillation curve of the floating frame.

    Figure 15.Thruster ignition.

    Figure 16.The relationship between levitation force and air gap.

    Figure 17.Thruster ignition test.

    The PSD of floating frame oscillation based on two different superconducting bearings can be compared in figures 12(b) and 13(b).It can be found that between 1 mHz and 10 Hz, as shown in the blue box, the maximum noise value of floating frame based on the superconducting composite bearing is 1 mm/Hz0.5,while the value is 10 mm/Hz0.5for single axial superconducting bearing.The average noise of superconducting composite bearing is also smaller than that of the single axial bearing above the level of 1 μm/Hz0.5.

    Before ignition, the curve in the free oscillation state needs to be captured to calculate the natural periodTand use formula (5) to calculate the restoring torque coefficient.The free oscillation curve is shown in figure 14.

    It can be seen from the figure 14 that the curve is sinusoidal, and the amplitude is basically unchanged, indicating that the floating frame experiences extremely low damping and noise in the testing conditions.The curve fitting equation is:

    The authors carried out five repetitive calibration tests to obtain five oscillation periods, and calculated the average value via equation (10)

    According to equations (9) and (10), the average free oscillation period of floating frame can be obtained as ?T=30.206 s.

    The electromagnetic damper is turned on to weak the oscillation amplitude of the floating frame.The floating frame still has a small oscillation near the equilibrium point after 30 mins.The authors turn off the electromagnetic damper and continue to ignite, while the testbed works in open-loop mode.The displacement change of floating frame is collected via the motion detector.The thruster ignition picture is shown in figure 15.

    4.3.Results and discussion

    In order to study the influence of FC height on the levitation force, the experiments measured and analyzed the levitation forces when the FC height was 5.5 mm, 13.5 mm, and 26.5 mm (zero FC).The results are shown in figure 16.

    It can be seen from the figure 16 that the vertical levitation force curves in the three cases generally have the same trend, which shows the higher FC height brings the larger levitation force.In the horizontal direction, the test results of the radial force show the opposite trend to the vertical levitation force.The results show that the higher FC height causes the smaller radial force.The reason for this phenomenon is that the radial force is mainly the flux pinning force, and the larger axial cooling height will cause a bigger levitation force and but a smaller pinning force.

    At the same time, the stiffness characteristics of HTS bearing with FC heights of 5.5 mm,13.5 mm and 26.5 mm are analyzed via the levitation force curves.It can be concluded that the FC height has a great influence on the levitation stiffness of HTS bearing,and the larger FC height leads to the smaller stiffness for both vertical and horizontal directions.

    The floating frame can remain stable under the field cold height of 8 mm when the mass is 4 kg,and the corresponding levitation air gap is 2 mm, which can provide about 8 N mm?1stiffness for the floating frame.

    After the floating frame is under stable levitation and the oscillation amplitude is suppressed by the damper,the authors turn off the electromagnetic damper and ignite after 200 s,and the thruster continues to ignite for 400 s.It can be clearly seen from the figure 17 that the oscillation equilibrium point of floating frame is changed after ignition.According to the calculation, the change angleθof equilibrium point is 0.75°.The inertia momentIof the floating frame is 0.00445 kg m2,and the equivalent thrust radius is 10 cm.By substituting the above parameters into equations (5) and (2), we can get the thrust was about 25.2μN(yùn).

    The resolution of this measurement testbed mainly depends on the accuracy of laser displacement sensor according to the low zero drift state of floating frame.The resolution of laser displacement sensor is 1 μm, and the corresponding angular deflection resolution is 0.03°.While the elasticity coefficient ˉkin this article is 33.602 μN(yùn)/°, so the resolution of the testbed is: 33.602 μN(yùn)/°×0.03°=1.008 μN(yùn).

    The authors analyzed the data of five repeated experiments in order to obtain the error of the measurement testbed.The data is shown in table 2.

    The standard deviation of the sample with respect to the restoring moment coefficientkcan be determined by equation (11)for a nysatisfies the 3σ principle.

    The value of 3Sk is 2.37 μN(yùn)/°,which can conclude that Thus,it is measured from the free oscillation experiment that the value of=33.602±2.37 μN(yùn)/°, corresponding thrust is 25.2±1.78 μN(yùn).

    5.Conclusion

    This paper introduces a micro-thrust measuring method based on the magnetic levitation thrust testbed with the key design of an HTS composite bearing, and the new testbed was designed to measure the steady thrust of miniature ion electrospray thruster.An electromagnetic finite simulation was executed to give an optimized design for the HTS composite bearing, the magnetic flux leakage and the electromagnetic interference of the bearing was reduced largely.A 2D measurement device was used to quantitatively analyze the axial and radial levitation forces for the HTS composite bearing.The results show that the HTS bearing has the levitation force larger than 40 N, suggesting that the thrust testbed using the bearing has a load capacity up to 4 kg.Furthermore, the free oscillation curve indicates that the testbed has low noise by introducing HTS composite bearing to increase the radial stiffness.Finally, the floating frame has an obvious response to the ignition of ion electrospray thruster, and the analyzed result is about 25.2μN(yùn) with an error of 1.78μN(yùn).The magnetic levitation thrust testbed also shows a low noise background because of the small radial offset of floating frame during the measurement.

    猜你喜歡
    趙鵬海濱
    《平行四邊形》拓展精練
    夏日海濱
    Structural Reliability Analysis Based on Support Vector Machine and Dual Neural Network Direct Integration Method
    海濱書簡(jiǎn)
    “同分異構(gòu)體”知識(shí)解談
    尋找失蹤的少女
    海濱1
    尋找失蹤的少女(上)
    社團(tuán)少年
    社團(tuán)少年
    国内精品久久久久久久电影| 午夜福利视频1000在线观看| 五月玫瑰六月丁香| 国产私拍福利视频在线观看| 精品久久久久久,| 亚洲专区中文字幕在线| 动漫黄色视频在线观看| h日本视频在线播放| 国产色婷婷99| 三级国产精品欧美在线观看| 88av欧美| 能在线免费观看的黄片| 国产精品人妻久久久久久| 一区二区三区激情视频| 美女黄网站色视频| 亚洲精品久久国产高清桃花| 日韩欧美国产在线观看| 在线天堂最新版资源| 国产精品av视频在线免费观看| 最近视频中文字幕2019在线8| 日韩 亚洲 欧美在线| 啦啦啦观看免费观看视频高清| 日韩中文字幕欧美一区二区| 亚洲av免费在线观看| 久久久久国内视频| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 日本黄色视频三级网站网址| 日本色播在线视频| 男女视频在线观看网站免费| 国产av一区在线观看免费| 亚洲av五月六月丁香网| 美女免费视频网站| 性插视频无遮挡在线免费观看| 特大巨黑吊av在线直播| 少妇高潮的动态图| 男女做爰动态图高潮gif福利片| 久久6这里有精品| 成人鲁丝片一二三区免费| 亚洲国产日韩欧美精品在线观看| 91av网一区二区| 听说在线观看完整版免费高清| 少妇人妻精品综合一区二区 | 91麻豆精品激情在线观看国产| 欧美在线一区亚洲| 精品久久久久久久久久久久久| 亚洲av二区三区四区| 欧美日韩精品成人综合77777| 淫秽高清视频在线观看| 夜夜夜夜夜久久久久| av在线亚洲专区| 一区二区三区高清视频在线| 精品久久国产蜜桃| 亚洲成a人片在线一区二区| netflix在线观看网站| 男女啪啪激烈高潮av片| 亚洲性夜色夜夜综合| 国产精品一区www在线观看 | 久久久久久久精品吃奶| 蜜桃久久精品国产亚洲av| 欧美激情久久久久久爽电影| 国产伦精品一区二区三区视频9| 啪啪无遮挡十八禁网站| 中文在线观看免费www的网站| av天堂中文字幕网| 国产精品一及| 国产精品一区www在线观看 | 亚洲精品久久国产高清桃花| 午夜久久久久精精品| 欧美最新免费一区二区三区| 国产精品,欧美在线| 一区二区三区免费毛片| 日韩欧美国产在线观看| 国产免费一级a男人的天堂| 国产男人的电影天堂91| 亚洲欧美日韩高清专用| 免费观看精品视频网站| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 久久99热6这里只有精品| 欧美+日韩+精品| 少妇人妻一区二区三区视频| 欧美黑人巨大hd| 久久草成人影院| 亚洲最大成人中文| 久久亚洲精品不卡| 亚洲最大成人手机在线| 亚洲成av人片在线播放无| 亚洲美女黄片视频| 亚洲精华国产精华精| 精品久久久久久久久久久久久| 全区人妻精品视频| 成人特级av手机在线观看| 亚洲人成伊人成综合网2020| 国内精品一区二区在线观看| 国产主播在线观看一区二区| 午夜免费激情av| 淫妇啪啪啪对白视频| 国产女主播在线喷水免费视频网站 | 亚洲真实伦在线观看| 欧美性感艳星| 亚洲avbb在线观看| 亚洲中文字幕日韩| 亚洲精品国产成人久久av| 97超级碰碰碰精品色视频在线观看| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久精品电影| 一本久久中文字幕| 亚洲精华国产精华液的使用体验 | 日韩人妻高清精品专区| 听说在线观看完整版免费高清| 国产精品伦人一区二区| 夜夜看夜夜爽夜夜摸| 亚洲熟妇熟女久久| 日本黄色片子视频| 欧美日韩瑟瑟在线播放| 国产成年人精品一区二区| 国产真实乱freesex| 午夜视频国产福利| 少妇人妻一区二区三区视频| 毛片女人毛片| 国产成人一区二区在线| 久久久久久伊人网av| 亚洲欧美日韩高清专用| 日韩国内少妇激情av| 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 波多野结衣巨乳人妻| 午夜影院日韩av| 两个人视频免费观看高清| 直男gayav资源| 日本成人三级电影网站| 一a级毛片在线观看| 最近视频中文字幕2019在线8| 两个人视频免费观看高清| 蜜桃久久精品国产亚洲av| 一区二区三区激情视频| 超碰av人人做人人爽久久| 亚洲人成网站在线播| 在线观看66精品国产| 深夜精品福利| 一级av片app| 欧美最新免费一区二区三区| 国产精品爽爽va在线观看网站| 免费看av在线观看网站| 国产极品精品免费视频能看的| 欧美一区二区国产精品久久精品| 91麻豆av在线| 国内精品美女久久久久久| 91久久精品国产一区二区三区| 国产精品久久久久久久久免| 身体一侧抽搐| 成人午夜高清在线视频| 免费观看精品视频网站| 最新中文字幕久久久久| 亚洲人与动物交配视频| 日韩国内少妇激情av| 亚洲综合色惰| 国产亚洲欧美98| 免费无遮挡裸体视频| 赤兔流量卡办理| 国产精品无大码| 日本 欧美在线| 可以在线观看毛片的网站| 久久人人精品亚洲av| 性插视频无遮挡在线免费观看| 成人综合一区亚洲| avwww免费| 国产精品嫩草影院av在线观看 | 九色成人免费人妻av| 欧美3d第一页| 国产亚洲精品久久久久久毛片| 琪琪午夜伦伦电影理论片6080| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼| 亚洲av一区综合| 国产在视频线在精品| 日韩一本色道免费dvd| 两个人视频免费观看高清| 久久久久久久亚洲中文字幕| 色在线成人网| 别揉我奶头~嗯~啊~动态视频| 国产伦一二天堂av在线观看| 国产亚洲91精品色在线| 婷婷丁香在线五月| 国产精品一区二区性色av| 国内毛片毛片毛片毛片毛片| 中文字幕人妻熟人妻熟丝袜美| 22中文网久久字幕| 中亚洲国语对白在线视频| 日韩一区二区视频免费看| 少妇被粗大猛烈的视频| 三级毛片av免费| 日韩精品中文字幕看吧| 亚洲无线在线观看| 99久久成人亚洲精品观看| 欧美+亚洲+日韩+国产| 免费人成视频x8x8入口观看| 亚洲成av人片在线播放无| 久久久成人免费电影| 国产成人aa在线观看| 伦理电影大哥的女人| 亚洲国产欧洲综合997久久,| 国产精品久久电影中文字幕| 亚洲精品一区av在线观看| 国产精品久久视频播放| 熟妇人妻久久中文字幕3abv| 国国产精品蜜臀av免费| 一本精品99久久精品77| 97人妻精品一区二区三区麻豆| 亚洲最大成人手机在线| 免费看光身美女| 亚洲精品色激情综合| 女人被狂操c到高潮| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 精品人妻一区二区三区麻豆 | 国产主播在线观看一区二区| 性欧美人与动物交配| 免费观看在线日韩| 亚洲国产精品成人综合色| 可以在线观看毛片的网站| 日韩精品有码人妻一区| 国产美女午夜福利| 欧美日韩黄片免| 亚洲七黄色美女视频| 联通29元200g的流量卡| avwww免费| 无人区码免费观看不卡| 亚洲中文字幕一区二区三区有码在线看| 超碰av人人做人人爽久久| 精品一区二区三区av网在线观看| 婷婷六月久久综合丁香| 成年女人毛片免费观看观看9| 午夜福利在线观看免费完整高清在 | 国产一区二区三区av在线 | 久久亚洲真实| 国产精品乱码一区二三区的特点| 亚洲无线在线观看| 99精品在免费线老司机午夜| www.色视频.com| 久久久久性生活片| 国产高清视频在线播放一区| 一区二区三区免费毛片| 少妇高潮的动态图| 一夜夜www| 国产精品电影一区二区三区| 88av欧美| 老师上课跳d突然被开到最大视频| 2021天堂中文幕一二区在线观| 国模一区二区三区四区视频| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 国产欧美日韩精品一区二区| 97超视频在线观看视频| 日本在线视频免费播放| www日本黄色视频网| 狂野欧美白嫩少妇大欣赏| 国产视频一区二区在线看| 欧美性猛交╳xxx乱大交人| 久久久精品大字幕| 九九久久精品国产亚洲av麻豆| 一级黄片播放器| 国产午夜福利久久久久久| 丝袜美腿在线中文| 熟妇人妻久久中文字幕3abv| 舔av片在线| 黄色欧美视频在线观看| 久久精品国产亚洲av香蕉五月| 日本撒尿小便嘘嘘汇集6| 午夜福利在线观看免费完整高清在 | 97人妻精品一区二区三区麻豆| 国产成人aa在线观看| 亚洲人成网站在线播放欧美日韩| 免费搜索国产男女视频| 日韩,欧美,国产一区二区三区 | 伦理电影大哥的女人| 我要看日韩黄色一级片| 麻豆一二三区av精品| 高清日韩中文字幕在线| 精品久久久久久久久久久久久| 亚洲国产精品sss在线观看| 好男人在线观看高清免费视频| 国产精品98久久久久久宅男小说| 免费人成视频x8x8入口观看| 亚洲国产精品合色在线| 天美传媒精品一区二区| 久久久久久久久久成人| 免费大片18禁| 亚洲性久久影院| 成年女人永久免费观看视频| 一本精品99久久精品77| 久久久久久伊人网av| 成熟少妇高潮喷水视频| 日韩欧美精品v在线| 精品人妻熟女av久视频| 久久久久久九九精品二区国产| 嫁个100分男人电影在线观看| 日韩,欧美,国产一区二区三区 | 国产视频一区二区在线看| 国产高清有码在线观看视频| 内射极品少妇av片p| 综合色av麻豆| 亚洲中文日韩欧美视频| 色哟哟·www| 国产aⅴ精品一区二区三区波| 欧美bdsm另类| 中出人妻视频一区二区| 一进一出抽搐gif免费好疼| 亚洲五月天丁香| 特级一级黄色大片| 国产黄色小视频在线观看| 日韩一区二区视频免费看| 欧美黑人欧美精品刺激| 男人狂女人下面高潮的视频| 久久人人精品亚洲av| 午夜激情福利司机影院| 国产乱人视频| 亚洲最大成人手机在线| 午夜福利在线观看免费完整高清在 | 搡老岳熟女国产| 两人在一起打扑克的视频| 精品人妻偷拍中文字幕| 成人av在线播放网站| 午夜福利在线观看免费完整高清在 | 久久国产乱子免费精品| 国产精品久久久久久久久免| 3wmmmm亚洲av在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人精品中文字幕电影| 最新在线观看一区二区三区| 精品人妻视频免费看| 午夜激情欧美在线| 亚洲精品亚洲一区二区| 国产亚洲精品久久久com| 精品不卡国产一区二区三区| 成人无遮挡网站| 国产成人一区二区在线| 18禁在线播放成人免费| 国产伦在线观看视频一区| 日韩强制内射视频| 亚洲成av人片在线播放无| 九九热线精品视视频播放| 美女被艹到高潮喷水动态| 亚洲熟妇中文字幕五十中出| 日日夜夜操网爽| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲在线自拍视频| 天堂√8在线中文| 国内精品久久久久精免费| 深夜精品福利| 国产精品电影一区二区三区| 久久久久久久午夜电影| 欧美高清成人免费视频www| 久久久久久久午夜电影| 国产精品不卡视频一区二区| 少妇裸体淫交视频免费看高清| 91狼人影院| 国产精品自产拍在线观看55亚洲| 久久久久免费精品人妻一区二区| 亚洲av一区综合| 伦精品一区二区三区| 最近视频中文字幕2019在线8| 国产亚洲精品综合一区在线观看| 俄罗斯特黄特色一大片| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 欧美日韩国产亚洲二区| 嫩草影院新地址| 国产亚洲精品av在线| 日韩强制内射视频| 熟妇人妻久久中文字幕3abv| 日韩av在线大香蕉| 一级黄色大片毛片| 午夜福利视频1000在线观看| 成人av一区二区三区在线看| x7x7x7水蜜桃| 午夜影院日韩av| 又爽又黄无遮挡网站| 国产 一区精品| 久久久国产成人精品二区| 亚洲国产精品sss在线观看| 国产精品久久久久久久久免| 搡女人真爽免费视频火全软件 | 国产av麻豆久久久久久久| 九色成人免费人妻av| 99久久精品一区二区三区| 亚洲精品粉嫩美女一区| 欧美日本视频| 天美传媒精品一区二区| 国语自产精品视频在线第100页| 毛片一级片免费看久久久久 | 国产黄色小视频在线观看| 不卡视频在线观看欧美| 国产精品久久久久久av不卡| 全区人妻精品视频| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 可以在线观看的亚洲视频| 成人毛片a级毛片在线播放| 亚洲四区av| 精品一区二区三区人妻视频| 国产成人aa在线观看| 国产精品久久久久久亚洲av鲁大| 性欧美人与动物交配| 国产 一区精品| 极品教师在线视频| 一进一出好大好爽视频| 欧美一级a爱片免费观看看| 国产男人的电影天堂91| 日韩欧美在线乱码| 欧美成人性av电影在线观看| 女同久久另类99精品国产91| 久久久色成人| 人人妻,人人澡人人爽秒播| 夜夜夜夜夜久久久久| 亚洲成人中文字幕在线播放| 精品久久久久久久久久久久久| 亚州av有码| 国产欧美日韩精品一区二区| 国内毛片毛片毛片毛片毛片| 国产黄片美女视频| 禁无遮挡网站| 在线观看av片永久免费下载| 久久精品国产99精品国产亚洲性色| 欧美日韩精品成人综合77777| 日韩欧美在线二视频| 久久中文看片网| 成人综合一区亚洲| 极品教师在线免费播放| 黄色配什么色好看| 春色校园在线视频观看| 欧美一区二区精品小视频在线| 无人区码免费观看不卡| 精华霜和精华液先用哪个| 九九爱精品视频在线观看| 国产91精品成人一区二区三区| 中文字幕久久专区| 国产精品精品国产色婷婷| 欧美日韩瑟瑟在线播放| 狂野欧美激情性xxxx在线观看| 窝窝影院91人妻| av在线天堂中文字幕| 又黄又爽又免费观看的视频| 高清毛片免费观看视频网站| 久久精品人妻少妇| 精品日产1卡2卡| 国产高清激情床上av| 国产综合懂色| 淫秽高清视频在线观看| 国产精品一区二区三区四区免费观看 | 男插女下体视频免费在线播放| 在线天堂最新版资源| 国产人妻一区二区三区在| 2021天堂中文幕一二区在线观| 校园春色视频在线观看| 在线免费观看的www视频| 精品国内亚洲2022精品成人| 在线观看美女被高潮喷水网站| 亚洲色图av天堂| 国产亚洲精品久久久com| 成人国产一区最新在线观看| 免费看av在线观看网站| 亚洲综合色惰| 国产在线男女| 丝袜美腿在线中文| 日本黄大片高清| 久久精品国产亚洲网站| 男女下面进入的视频免费午夜| 国产高清三级在线| 综合色av麻豆| h日本视频在线播放| 一本一本综合久久| av在线观看视频网站免费| 九九爱精品视频在线观看| 亚洲国产欧洲综合997久久,| 欧美日韩国产亚洲二区| 在线观看午夜福利视频| 大型黄色视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品国产三级国产av玫瑰| 午夜视频国产福利| 亚洲黑人精品在线| 悠悠久久av| 国产单亲对白刺激| 美女黄网站色视频| 欧美日本视频| 久久久久国内视频| 波多野结衣高清无吗| 欧美成人免费av一区二区三区| 91麻豆精品激情在线观看国产| 午夜视频国产福利| 女人十人毛片免费观看3o分钟| 日本成人三级电影网站| 乱系列少妇在线播放| 国产 一区精品| 久久国产精品人妻蜜桃| 极品教师在线免费播放| 欧美区成人在线视频| 三级毛片av免费| 日韩精品有码人妻一区| 成人欧美大片| 狠狠狠狠99中文字幕| 永久网站在线| 久久久国产成人精品二区| 热99在线观看视频| 真人做人爱边吃奶动态| 欧美黑人巨大hd| 日本免费一区二区三区高清不卡| 精品99又大又爽又粗少妇毛片 | 内射极品少妇av片p| 18禁裸乳无遮挡免费网站照片| 日韩精品青青久久久久久| 久久久久国内视频| 日日啪夜夜撸| 黄色丝袜av网址大全| 久久久久久久午夜电影| 国产精品98久久久久久宅男小说| 51国产日韩欧美| 男人舔女人下体高潮全视频| av专区在线播放| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看影片大全网站| 精品日产1卡2卡| 日本 av在线| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 免费在线观看日本一区| 男女边吃奶边做爰视频| 小蜜桃在线观看免费完整版高清| 听说在线观看完整版免费高清| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区亚洲| 内射极品少妇av片p| 性插视频无遮挡在线免费观看| 人妻夜夜爽99麻豆av| 国产v大片淫在线免费观看| 久久久久久久久久黄片| 内射极品少妇av片p| 午夜精品在线福利| 国产真实伦视频高清在线观看 | 深夜精品福利| 伊人久久精品亚洲午夜| 国产淫片久久久久久久久| 亚洲精品成人久久久久久| 啦啦啦观看免费观看视频高清| 狠狠狠狠99中文字幕| 99久久久亚洲精品蜜臀av| 中出人妻视频一区二区| 国产黄a三级三级三级人| 成人国产综合亚洲| 日韩 亚洲 欧美在线| 黄色女人牲交| 日日摸夜夜添夜夜添av毛片 | 人人妻,人人澡人人爽秒播| av黄色大香蕉| 亚洲狠狠婷婷综合久久图片| 国产精品久久久久久久电影| 性欧美人与动物交配| 亚洲熟妇中文字幕五十中出| 国产精品国产高清国产av| 精品人妻熟女av久视频| www.色视频.com| av在线亚洲专区| 国产精品嫩草影院av在线观看 | 国产亚洲欧美98| 欧美黑人欧美精品刺激| 日韩精品中文字幕看吧| 三级毛片av免费| 国产精品一及| 一个人看视频在线观看www免费| 精品一区二区三区视频在线| 国产高清激情床上av| 国产亚洲精品综合一区在线观看| 日本一二三区视频观看| 成人三级黄色视频| 99精品在免费线老司机午夜| 色精品久久人妻99蜜桃| 午夜福利18| 联通29元200g的流量卡| 最近最新免费中文字幕在线| 国产三级在线视频| 亚洲三级黄色毛片| 亚洲精华国产精华精| 亚洲美女黄片视频| 免费在线观看日本一区| 亚洲精华国产精华精| 亚洲乱码一区二区免费版| 午夜精品一区二区三区免费看| 久久6这里有精品| 中文字幕av成人在线电影| 亚洲天堂国产精品一区在线| 精品免费久久久久久久清纯| 亚洲专区中文字幕在线| 国产又黄又爽又无遮挡在线| 嫩草影院新地址| 又粗又爽又猛毛片免费看| 国产又黄又爽又无遮挡在线| 国产高清有码在线观看视频| 国产三级中文精品| 久久人妻av系列| 久久久久久久精品吃奶| 精品人妻1区二区| 中文字幕久久专区| 一区二区三区免费毛片| 成人三级黄色视频| 波多野结衣巨乳人妻| 99久久精品国产国产毛片| 午夜a级毛片| 亚洲性久久影院| 国产精品一区www在线观看 | 久久久久久久久中文| 久久九九热精品免费|