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

    Micro-Cantilever Electric Field Sensor Driven by Electrostatic Force

    2023-11-14 07:44:20ZhifeiHnJunHuLichengLiJinlingHe
    Engineering 2023年5期

    Zhifei Hn, Jun Hu, Licheng Li, Jinling He,*

    a State Key Lab of Power Systems, Department of Electrical Engineering, Tsinghua University, Beijing 100084, China

    bChina Southern Power Grid, Guangzhou 510623, China

    Keywords:Microsensors Electric sensing devices Microfabrication Piezoresistive devices

    ARTICLEINFO With the development of smart grids and the energy internet,large-scale monitoring of voltage and electric field data is required in all aspects of power systems, which requires the arrangement of various advanced sensors.Measurement of the electric field can replace traditional voltage transformers to realize the non-contact measurement of voltage,which reduces the insulation cost and the difficulty of operation and maintenance.Electric field measurement can also be applied in various other areas, such as equipment fault diagnosis, lightning warning, and electromagnetic environment measurement.Traditional electric field measurement devices, such as field mills, are bulky and costly, so they cannot be arranged flexibly on a large scale.In this paper, we present an electrostatically actuated microelectric field sensor (E-sensor) with a piezoresistive sensing structure.The presented E-sensor is fabricated into a four-cantilever structure using microfabrication technology.The cantilevers are displaced under the drive of the electrostatic force,and the generated strain is transformed into measurable signals through piezoresistive materials.The presented E-sensor has the advantages of small size, low cost, low power consumption, and easy mass production.Moreover, the E-sensor has a high signal-to-noise ratio,high resolution,and wide measuring range.The experimental results show that the E-sensor has a linear electric field measurement range from 1.1 to 1100.0 kV·m-1 with an alternating current(AC)resolution of up to 112 V·m-1·Hz-1/2 and a cut-off frequency of 496 Hz, making it suitable for most applications in smart grids and the energy internet.

    1.Introduction

    Large-scale information monitoring is one of the foundations for the realization of smart grids and the energy internet [1–3].Through the large-scale collection of voltage and current data, a control strategy for power systems can be optimized, and status monitoring and fault diagnosis of power systems and equipment can be realized [4].High voltage is traditionally measured by a voltage transformer, which is bulky and has poor performance.Moreover, the need for insulation increases the cost of measurement and the difficulty of operation and maintenance.Given the need for an increasing number of sensing nodes, such largevolume and high-cost voltage-measuring equipment can no longer meet the demand.According to the Gauss–Legendre algorithm,voltage can be reversed by measuring the electric field at several points [5,6].Therefore, with the use of an electric field sensor (Esensor), non-contact voltage measurement can be realized, which can reduce the cost of measurement and increase the quantity of sensing nodes [7,8].In addition to voltage measurement, Esensors can be used for the fault monitoring of electrical equipment [9–11].It is difficult to implement status monitoring and fault diagnosis of electrical equipment online.When insulation failure or partial discharge occurs in electrical equipment, the internal electric field of the equipment will change.By arranging micro-E-sensors and extracting fingerprint characteristics of the electric field signal, such faults can be monitored and diagnosed online.

    E-sensors have important applications in fields other than power systems as well, as shown in Fig.1 [12].Wearable Esensors can be used for near-electricity warnings for electric workers[13].E-sensors can also be used to monitor the electric field of the atmospheric, in order to realize lightning warning, weather forecasting, and the selection of space launch windows [14–17].Furthermore, measurement of the electrostatic field can be used in industrial production as a means of safety monitoring,especially in the petrochemical industry.In addition, E-sensors can be used for cabin monitoring in advanced electric vehicles and aircraft[18].Application in such areas requires E-sensors to have high resolution, a wide measurable range, and low power consumption, in addition to a small size and low cost.

    Fig.1.Applications of micro-E-sensors.E: electric field.

    At present,electric field measurement presents numerous difficulties.For one thing, large and high-permittivity measuring devices will cause electric field distortion.For another, a complex electromagnetic environment often results in high coupling interference.In a power system,the amplitude of electric field measurement can reach the order of megavolts per meter, and the measurement frequency varies from direct current (DC) to kilohertz.At the same time, in order to meet the needs of most application scenarios in a power grid, the resolution of electric field measurement must be within 200 V·m-1.Field mills are traditionally used in the measurement of electric fields.However,due to the large size and low accuracy,this approach can no longer meet current needs.Microsensors are small advanced sensors with a low cost and high performance that can easily be mass produced using microfabrication technology.The principles behind micro-Esensors include the electro-optical effect, charge induction, and the piezoelectric effect.The measurement of optical E-sensors is based on the change in the refractive index of electro-optic crystals under electric fields[19–23].Optical E-sensors have a high resolution and wide frequency bandwidth and are currently used in applications such as bus voltage measurement in substations.Although such E-sensors are small, they need to be connected to expensive and bulky light sources and receivers through optical fibers, which cannot be arranged flexibly in equipment or on a large scale.

    Electric fields can also be measured using the induced charges of a conductor in an electric field [24–27].Micro-electromechanical system (MEMS) E-sensors based on charge induction are composed of a vibrating shield electrode, a fixed electrode for inducing charge, and back-end circuits.Such sensors have the advantages of small size and low cost.However, they often introduce high noise interference,and it is difficult to improve their frequency response and resolution, due to the limitations of their sensing area.Micro-E-sensors based on the piezoelectric effect measure electric fields by sensing the deformation of piezoelectric materials under an electric field[28–31].Such E-sensors are small and can be easily integrated.Nevertheless,the fabrication technology for piezoelectric materials and issues with the sensor performance—such as resolution and the signal-to-noise ratio—limit the application of such sensors.

    In this paper, we present an electrostatically actuated micro-E-sensor with a piezoresistive sensing structure.The sensor is fabricated into a four-cantilever structure with grounded metal films arranged at the ends of the cantilevers.Driven by electrostatic force,the cantilevers bend up in an electric field.The strain generated by the bending of the cantilever causes a change in the resistance of the piezoresistive materials at the roots of the cantilevers.The electric field is measured according to the output signal of a Wheatstone bridge circuit with piezoresistive materials arranged on its arms.Compared with optical E-sensors, the presented Esensor has a small size of 8 mm×8 mm and can be mass produced at a low cost using silicon-based microfabrication technology.Compared with MEMS E-sensors, the presented E-sensor does not require additional shielded electrode drive sources.In addition,since the response frequency of the E-sensor is twice the frequency of the electric field,the coupling interference of the E-sensor can be effectively reduced.Moreover,because the harmonics in the power grid are mostly odd harmonics, the output signal at a double frequency makes the sensor more suitable for power systems.The experimental results show that the cut-off frequency of the Esensor is as high as 496 Hz,and the measurable electric field ranges from 1.1 to 1100.0 kV·m-1with an alternating current(AC)resolution of 112 V·m-1·Hz-1/2.The performance of the presented Esensor meets most of the application requirements for smart grids and the energy internet.

    2.Working principle and structural design

    When a metal film is placed in an electric field, the free charge carriers will redistribute, resulting in opposite charges at the surface of the metal film,as shown in Fig.2(a).This polarization compensates for the electric field inside the metal.In an external electric field, these generated charges are subject to an outward electrostatic field force,while the resultant force on the metal film remains at zero.If the lower surface of the metal film is grounded,the charges generated on the lower surface are transferred away,and the metal film is subject to an upward electrostatic force.When a metal film is prepared on the free end of a silicon cantilever, the cantilever will bend up under the drive of the metal film,as shown in Fig.2(b).The electrostatic force(Fes)on the metal film can be calculated as follows:

    Fig.2.Response of a silicon cantilever with metal film in an electric field.(a)Induced charges of the metal film;(b)force of the grounded metal film and displacement of the silicon cantilever.Fes: electrostatic force.

    where A is the area of the metal film,E is the electric field along the normal direction of the film, and ε is the dielectric constant of air.

    The structure of the presented E-sensor is shown in Fig.3.The sensor is designed as a four-cantilever structure, where the four cantilevers are symmetrically distributed.One end of each cantilever is fixed in the center of the silicon substrate column, while the other end of the cantilever can vibrate freely.A grounded metal film is prepared on the free end of each cantilever.Under the action of an electric field, the four cantilevers bend up with the same amplitudes, as shown in Fig.3(a), generating stress inside the silicon cantilevers.At the fixed end of the cantilevers, four zigzag piezoresistive areas are prepared using ion doping.The piezoresistive areas are connected by metal electrodes to form a Wheatstone bridge circuit, as shown in Fig.3(b).

    When the strain of the four silicon cantilevers changes, the energy band of the doped silicon also changes, moving the energy of the energy valley.The stress of the doped silicon satisfies the following:

    where σlis the stress along the length of the cantilever, σtis the stress perpendicular to the length of the cantilever, and ν is the Poisson ratio of the doped silicon.

    Under the silicon strain,the resistance of the four piezoresistive areas changes.Among the four piezoresistive areas, R1and R3are designed as a zigzag shape, which couples the strain along the length of the cantilever.R2and R4are designed as a straight strip shape, which couples the strain along the width of the cantilever.Thus, the resistance changes of the adjacent piezoresistive areas are reversed.The initial resistance values of the piezoresistive areas are designed to be the same to reduce zero drift.If we set the initial resistance of the four piezoresistive areas as R0and mark the resistivity changes in R1and R3as ΔRaand the resistivity changes in R2and R4as ΔRb, then:

    where πl(wèi)and πtare the piezoresistive coefficients in the longitudinal and transverse directions,respectively.In the designed E-sensor,the values of πl(wèi)and πtare +71.3 and -66.7, respectively.P(N, T) is the coefficient that depends on the semiconductor doping density N and the temperature T.

    When applying a constant voltage source Usat the input ports Us(+)and Us(-) of the Wheatstone bridge,according to Kirchhoff’s laws, the output of the Wheatstone bridge Uocan be calculated as follows:

    According to Eq.(4),under fixed doping conditions,the piezoresistive areas where the strain is the largest should be selected to increase the response of the E-sensor.

    Fig.3.Structure of the E-sensor.(a)Diagram of the E-sensor structure;(b)overhead view of the E-sensor structure and circuit.GND:ground;R1–R4:piezoresistive areas;Us:constant voltage source; Uo: output voltage.

    To verify the working principle of the presented E-sensor, we use the finite-element method (FEM) via Comsol Multiphysics to simulate the E-sensor response under different parameters and electric fields.To simplify the model, the surrounding fixed area of the E-sensor is omitted,while the fixed silicon column and four cantilevers are retained.The simulation uses the physics of electrostatics and solid mechanics and the multiphysics of electromechanical force.We set the cantilever to free boundary conditions.The boundary conditions of other areas are fixed.Fig.4 shows the response of the E-sensor when an electric field of 1 MV·m-1is applied along the z-axis.

    Figs.4(a) and (b) show the distribution of the cantilever displacement.Under the electrostatic force, the metal film at the end of the cantilever drives the cantilever to bend upward.The four cantilevers have the same displacement.It is worth noting that the force on the cantilever is always upward, regardless of whether the direction of the electric field is along the positive or negative z-axis.Figs.4(c) and (d) shows the distribution of the XX strain component on the cantilever.The YY strain component is symmetrical with the XX component.As shown in the figure,the largest strain component on the cantilever is located at the root of the cantilever.In order to maximize the sensor response,the ion-doped area is chosen to be close to the fixed area on the cantilever.

    According to simulations, the maximum strain of the E-sensor increases linearly with the electric field.When the electric field is applied in the reverse direction along the z-axis, the sensor response is the same as the electric field is applied along the zaxis.In addition, the response of the E-sensor is related to the length and thickness of the cantilever.Simulations can provide a theoretical reference for designing the optimal parameters and layout of the E-sensor.Increasing the maximum strain of the cantilever under the same electric field can increase the output of the Wheatstone bridge, thereby increasing the resolution of the E-sensor.The relationship between the maximum strain tensor XX(eXX)on the cantilever and the length and thickness of the cantilever is shown in Fig.5.

    According to the principle of leverage, under the same electrostatic force, increasing the length of the cantilever increases the torque.Therefore,under the same electric field,the longer the cantilever is,the greater the maximum strain on the cantilever will be,as shown in Fig.5(a).Within a reasonable range, increasing the cantilever length improves the response of the E-sensor.However,an excessive cantilever length both reduces the mechanical strength of the cantilever and makes it more difficult to prepare the device.Moreover, as the length of the cantilever increases,the resonant frequency of the cantilever drops, causing the cutoff frequency of the E-sensor to decrease.Considering various factors such as device fabrication, lifetime, and device performance,the length of the fabricated E-sensor cantilever was selected to be within the range of 1250–2250 μm.

    Similarly,reducing the thickness of the cantilever increases the maximum strain on the cantilever,thereby increasing the E-sensor response, as shown in Fig.5(b).However, a cantilever that is too thin will affect the fabrication difficulty and mechanical strength of the E-sensor.Thus, the cantilever thickness of our E-sensor was chosen to be 5 and 10 μm.

    Changing the metal film area affects the sensor response.According to Eq.(1), increasing the metal film area A increases the electrostatic force, thus increasing the sensor response.Therefore, the area of the metal film should be as large as possible.However, when the length of the metal film increases past a certain value, the increase in the sensor response becomes insignificant.Moreover, when the metal is close to the root of the cantilever, it may affect the fabrication of the piezoresistive material.Therefore, the metal film was fabricated at the end of the cantilever, and a length of 800 μm was chosen.Section 4 shows a performance comparison of E-sensors with different parameters.

    Fig.4.Simulated response of the E-sensor.(a) Distribution of the displacement; (b) displacement curve of the cantilever; (c) distribution of the strain tensor XX (eXX);(d) cantilever curve of the eXX.

    Fig.5.Relationship between the maximum eXX on the cantilever and the length and thickness of the cantilever.(a) Maximum eXX under different cantilever lengths l;(b) maximum eXX under different cantilever thicknesses Tc.

    3.Device fabrication

    Based on the theoretical analysis and simulation results of the designed E-sensor structure,the presented E-sensor was fabricated using silicon-based microfabrication technology.The fabrication process was based on 10 cm diameter silicon-on-insulator (SOI)wafers with a device layer thickness of 5 and 10 μm and a box thickness of 1 μm.The crystal orientation of the device layer was selected to be <100>and the resistivity of the device silicon was 6–8 Ω·cm.About 90 E-sensors can be produced on one 10 cm diameter wafer.

    A schematic illustration of the E-sensor fabrication process is provided in Fig.6.The wafer is patterned and fabricated through a combination of photolithography technology and other microfabrication technologies.The six photolithography masks used in the fabrication process are shown in Fig.7(a).After the SOI wafer is cleaned with a mixture of hot concentrated sulfuric acid and hydrogen peroxide, the piezoresistive area is patterned onto the device layer.In order to ion dope the piezoresistive area,the wafer is heated in an oxygen atmosphere to form a 20 nm silicon oxide layer on the surface.The oxide layer acts as a shielding layer for the ion doping.Next,boron ions are doped onto the device layer with a concentration of 5×1018cm-3,forming a PN junction with a depth of 1.5 μm.The piezoresistive areas are created within the designed shape using photolithography, and ohmic contact areas are provided to ensure good ohmic contact between the piezoresistive areas and the electrodes.The wafer is annealed at 1000 °C for 30 min in a nitrogen atmosphere to activate the doped areas.After the piezoresistive areas are fabricated,a silicon oxide insulating layer with a thickness of 400 nm is grown on the surface of the wafer in an oxygen atmosphere.The insulating layer is then patterned using inductively coupled plasma(ICP)etching technology to release the silicon oxide in the ohmic contact area and the cantilever area.Finally,a metal layer composed of a chromium(Cr)layer and a gold(Au)layer is prepared on the surface of the device layer by means of magnetron sputtering.The metal is patterned into a Wheatstone bridge shape using photolithography and wet etching.In this way,the processing of structure on the device layer surface is completed.

    Fig.6.Schematic illustration of the E-sensor fabrication.PCB:printed circuit board.

    Fig.7.The fabricated micro-E-sensor.(a) Masks for photolithography; (b) overall view of the 10 cm diameter wafer before dicing; (c) microscopic image of the Esensor structure; (d) diagram showing the E-sensor size.

    Processing in the depth direction of the wafer is then performed.First, the device layer is patterned using ICP etching to form the shape of the cantilevers.Then, the handle layer under the cantilever structure is etched away.Finally, the box layer is etched away to release the four cantilevers.During the cantilever release,the temperature and time of the processing must be accurately controlled to ensure that the cantilevers have the same initial state.The fabricated wafer before dicing is shown in Fig.7(b).The microstructure of the sensor chip is shown in Fig.7(c).

    Finally, the wafer is diced into 8 mm × 8 mm chips using laser cutting,as shown in Fig.7(d).The metal electrodes on the E-sensor are welded to the printed circuit board (PCB) with gold wire.As discussed in Section 4,the PCB was used in an electric field test station to further test the E-sensor performance.

    4.Device performance

    Fig.8.Electric field measurement setup.(a) Schematic of the electric field measurement setup;(b) schematic of the PCB; (c) back-end differential amplification circuit.PC:personal computer; Sync: synchronization signal; f0: frequency of generator signal; RG: external gain setting resistor; Vin: input signal; Vout: output signal; VREF: output reference signal; Vs: positive supply terminal.

    An electric field test station was built for the performance testing of the E-sensors,as shown in Fig.8(a).The test station can generate electric fields with different field strengths, frequencies, and waveforms.In the test station, an electric field is generated by applying a voltage to a pair of circular metal electrodes.The ratio of the electrode diameter to the distance should be high enough to ensure a uniform electric field between the electrodes.A column for fixing the PCB board is arranged beside the two electrodes.The PCB board is designed as a ‘‘T” shape,as shown in Fig.8(b).The Esensor is placed on one end of the PCB,and the other end of the PCB is fixed by screws.The input voltage signal is generated by a generator (DG1022, Rigol Technology Inc., China) and amplified by a high-voltage amplifier (TREK MODEL 10/10B-HS, Trek Inc., USA)with a fixed gain of 1000.The output signal of the sensor is differentially amplified by the back-end circuit based on AD623, as shown in Fig.8(c).The back-end circuit is arranged on the PCB.The amplified signal is used for frequency selection through a lock-in amplifier (MFLI, Zurich Instruments, Switzerland), based on a generator signal frequency twice as high as the synchronization signal frequency.In addition, a high-precision voltage source(B2912A, Keysight Technology Inc., USA) is used to provide a 6 V input signal to the Wheatstone bridge on the sensor chip.Based on a LabVIEW program, the output signal of the generator and the output of the lock-in amplifier are controlled and collected through a personal computer (PC).

    Based on the electric test station, we tested the prepared Esensors with different parameters under electric fields with different amplitudes and frequencies.The tests were carried out under stable environmental conditions with the air pressure maintained at 101 kPa, the temperature at (25 ± 1) °C, and the humidity at 50% ± 10%.For electric field measurement, important parameters include the dynamic range,linearity,cut-off frequency,and resolution.The fabricated E-sensors were divided into four groups,marked as S1, S2, S3, and S4.Each group consisted of 10 sensor chips.The cantilever lengths of the E-sensors were 1250,1750,and 2250 μm, respectively.The cantilever thickness of the E-sensors was set to 5 and 10 μm.Table 1 displays the performance of Esensors with different parameters.

    The four groups of E-sensors were tested under different electric field strengths at 50 Hz.The output voltages of the E-sensors are shown in Fig.9(a).At a fixed frequency, the amplitude response of the E-sensors exhibits high linearity.The output voltage curve is linearly fitted,and the coefficient of determination(r2)values of S1–S4 are 0.998,0.998,0.998,and 0.996,respectively.In line with the theoretical analysis and simulation results,increasing the length or reducing the thickness of the cantilever increased the response of the E-sensor.Combining the response curves of the E-sensors with the actual sources of noise and interference,the AC resolution of the E-sensors was calculated, as shown in Table 1.When the electric field strength is too large or too small, the response of the E-sensor no longer exhibits linearity.When the electric field is high,the response curve of the sensor bends downward, because the response of the E-sensors is too large, and the working mode is changed.Under a low electric field,the influence of noise increases significantly.In addition, the interference of the environment and temperature causes the sensor to exhibit zero drift, which is manifested as the intercept of the output voltage curve.Taking the linearity of the response curve r2>0.99 as the limiting condition, the linear measurable electric field range(i.e.,the dynamic range)of the E-sensors can be obtained,as shown in Table 1.

    Fig.9(b)shows the frequency response curves of E-sensors with different parameters.Below the resonant frequency of the device,the E-sensor response curves are relatively flat.Taking a 3 dB change as the range, the cut-off frequencies of S1–S4 are 496,474, 441, and 296 Hz, respectively.The difference in the resonant frequencies of the sensors is caused by the difference in the size parameters of the sensor structure.Consistent with the theoretical analysis,reducing the length or increasing the thickness of the cantilever increases the resonance frequency of the mechanical structure, thereby increasing the cut-off frequency of the E-sensor.However, the increase in cut-off frequency is accompanied by a decrease in sensor response.Therefore, the size parameters of the E-sensor must be designed according to actual application scenarios.Above the resonance frequency, the response of the Esensor decreases rapidly.When the frequency is too low, the response of the E-sensor gradually decreases as well.It is worth noting that the sensor also responds to the DC electric field, but the output voltage will gradually decrease over time.

    Fig.9(c) shows the frequency response curve of S2 under different electric field strengths.As the electric field strength increases, the cut-off frequency of the E-sensor remains basically unchanged.When the field strength is low, the signal-to-noise ratio of the sensor is greatly affected by noise.At a high field strength, the signal-to-noise ratio of the E-sensor output signal is high, but it is limited by the maximum output voltage of the lock-in amplifier.

    The experimental results show that the E-sensor has a high performance and high consistency under different electric fields.The AC resolution of the E-sensor can reach 112 V·m-1·Hz-1/2, the linear dynamic range can reach 1.1–1100.0 kV·m-1, and the cut-off frequency is 496 Hz.Under normal operating conditions,the power consumption of the sensor is on the order of microwatts.In addition, the double frequency response enables the output signal of the E-sensor to be easily separated from the coupling signal,which greatly improves the signal-to-noise ratio.Table 2 compares the performance and characteristics of the presented E-sensor and Esensors based on other principles.The presented E-sensor exhibits good performance in terms of measurement range, resolution,elimination of coupled noise, and power consumption, and can meet the needs of most electric field measurement scenarios.In the future, the cut-off frequency of the presented E-sensor needs to be further improved.

    5.Conclusions

    In this paper, we presented an electrostatically actuated Esensor with a piezoresistive cantilever sensing structure.The presented E-sensor can be used for high-resolution electric field measurement in power systems and electrical equipment.Through theoretical calculation and a finite-element simulation method,the working principle of the E-sensor was verified,and the parameters of the E-sensor were optimized.We designed the fabrication process based on microfabrication technology and carried out sensor fabrication using 10 cm diameter SOI wafers.The experimental results showed that the presented E-sensor has high linearity and resolution.The AC resolution of the E-sensor can reach 112 V·m-1·Hz-1/2, the linear dynamic range can reach 1.1–1100.0 kV·m-1, and the cut-off frequency is 496 Hz; thus, the sensor can meet the needs of most applications in power grids.The presented E-sensor also has the advantages of small size,low cost,easy mass production, and low power consumption.Moreover, the output signal of the E-sensor can be separated from the coupled noise by a double frequency lock using a lock-in amplifier, improving the signal-to-noise ratio of the output signal.The cut-off frequency of the presented E-sensor still requires further improvement in the future.In addition,by integrating micro-E-sensors with an energyacquisition module, communication module, and data-processing module, the construction of high-performance sensor nodes can be realized.This high-performance electric field measurement system will play an important role in system detection and sensor network architecture in smart grids and the energy internet.

    Table 1Results of electric field measurement for E-sensors with different parameters.

    Fig.9.Response curves of E-sensors under electric fields.(a) Experimental static output Uo versus electric field E of E-sensors with different parameters; (b) frequency (f)response curves of E-sensors with different parameters; (c) frequency response curves of S2 under different electric field strengths.

    Table 2Comparison of the characteristics of E-sensors based on different principles.

    Acknowledgments

    This work was supported in part by the National Natural Science Foundation of China (51720105004, 51921005, and 52125703).

    Compliance with ethics guidelines

    Zhifei Han,Jun Hu,Licheng Li,and Jinliang He declare that they have no conflict of interest or financial conflicts to disclose.

    日韩精品有码人妻一区| 国产成人精品婷婷| 国产精品一国产av| 2022亚洲国产成人精品| 亚洲美女黄色视频免费看| 亚洲精品久久久久久婷婷小说| 波多野结衣一区麻豆| 国产免费现黄频在线看| 男女无遮挡免费网站观看| 免费看不卡的av| 午夜福利在线免费观看网站| 国产av一区二区精品久久| 少妇的逼水好多| 久久久久久免费高清国产稀缺| 久久久久久免费高清国产稀缺| 久久精品国产亚洲av天美| 最近中文字幕高清免费大全6| 免费高清在线观看日韩| 男人添女人高潮全过程视频| 九草在线视频观看| 国产福利在线免费观看视频| 亚洲精品久久久久久婷婷小说| 日韩不卡一区二区三区视频在线| 一本—道久久a久久精品蜜桃钙片| 久久这里只有精品19| 成年女人毛片免费观看观看9 | 国产精品二区激情视频| videosex国产| 国产极品天堂在线| 久久久精品国产亚洲av高清涩受| 中文字幕色久视频| 久久久精品免费免费高清| 女人被躁到高潮嗷嗷叫费观| 精品人妻在线不人妻| 亚洲久久久国产精品| 免费高清在线观看日韩| 两性夫妻黄色片| 日韩制服骚丝袜av| 久久热在线av| 夜夜骑夜夜射夜夜干| 91精品国产国语对白视频| 99久久综合免费| 久久精品国产综合久久久| 王馨瑶露胸无遮挡在线观看| 两个人看的免费小视频| 男女高潮啪啪啪动态图| 久久久久国产精品人妻一区二区| 一边摸一边做爽爽视频免费| 丁香六月天网| 亚洲在久久综合| 精品国产一区二区三区四区第35| 欧美精品高潮呻吟av久久| 国产免费又黄又爽又色| 久久亚洲国产成人精品v| tube8黄色片| 天天影视国产精品| 国产日韩欧美在线精品| 美女高潮到喷水免费观看| 日本黄色日本黄色录像| 亚洲精品aⅴ在线观看| 欧美日韩综合久久久久久| 在线天堂中文资源库| 一级爰片在线观看| 国产又色又爽无遮挡免| 成人午夜精彩视频在线观看| 人人妻人人爽人人添夜夜欢视频| 人人妻人人爽人人添夜夜欢视频| 满18在线观看网站| 亚洲在久久综合| 色94色欧美一区二区| 亚洲经典国产精华液单| 婷婷色麻豆天堂久久| 精品国产乱码久久久久久小说| 性少妇av在线| 欧美精品国产亚洲| 午夜福利视频精品| 久久99热这里只频精品6学生| 国产精品免费大片| 可以免费在线观看a视频的电影网站 | 日日爽夜夜爽网站| 精品少妇一区二区三区视频日本电影 | av网站在线播放免费| 一级片'在线观看视频| 国产成人av激情在线播放| av在线播放精品| 国产成人a∨麻豆精品| 免费在线观看视频国产中文字幕亚洲 | 日韩av不卡免费在线播放| 亚洲av欧美aⅴ国产| 制服诱惑二区| av又黄又爽大尺度在线免费看| 一级片'在线观看视频| 国产精品人妻久久久影院| 日本wwww免费看| 国产 一区精品| 一本久久精品| 国产成人精品一,二区| 老汉色av国产亚洲站长工具| 久久这里有精品视频免费| 丝袜在线中文字幕| 自线自在国产av| 国产成人午夜福利电影在线观看| 制服丝袜香蕉在线| 汤姆久久久久久久影院中文字幕| 日韩中文字幕视频在线看片| 综合色丁香网| 大片电影免费在线观看免费| 国产熟女午夜一区二区三区| av国产精品久久久久影院| 99久国产av精品国产电影| 在线亚洲精品国产二区图片欧美| 亚洲欧美成人精品一区二区| 国产福利在线免费观看视频| 国产精品亚洲av一区麻豆 | 国产综合精华液| 啦啦啦视频在线资源免费观看| 青春草亚洲视频在线观看| 国产高清国产精品国产三级| 久久国内精品自在自线图片| 青春草国产在线视频| 日日啪夜夜爽| 黄色一级大片看看| 国产极品天堂在线| 免费av中文字幕在线| 热re99久久国产66热| 最近2019中文字幕mv第一页| 一本久久精品| 美女视频免费永久观看网站| 一级黄片播放器| 日韩中文字幕欧美一区二区 | 天美传媒精品一区二区| 成年人午夜在线观看视频| 999精品在线视频| 日本91视频免费播放| 亚洲三级黄色毛片| 大陆偷拍与自拍| 青草久久国产| 久久影院123| 一本久久精品| 女人高潮潮喷娇喘18禁视频| 国产一区二区三区综合在线观看| 中文字幕人妻丝袜制服| 欧美精品国产亚洲| 一区二区三区乱码不卡18| 春色校园在线视频观看| 亚洲图色成人| 在线天堂中文资源库| 亚洲激情五月婷婷啪啪| 国产精品人妻久久久影院| 亚洲精品国产色婷婷电影| 国产日韩一区二区三区精品不卡| 国产激情久久老熟女| 热re99久久国产66热| 婷婷色av中文字幕| 久久婷婷青草| 伦精品一区二区三区| 视频在线观看一区二区三区| 青春草亚洲视频在线观看| 亚洲国产色片| 夜夜骑夜夜射夜夜干| 黑人巨大精品欧美一区二区蜜桃| 亚洲第一区二区三区不卡| 国产国语露脸激情在线看| 久久人人爽av亚洲精品天堂| 建设人人有责人人尽责人人享有的| 国产日韩欧美视频二区| 国产成人精品一,二区| 香蕉丝袜av| 一级毛片黄色毛片免费观看视频| 少妇人妻 视频| 超碰成人久久| 美女午夜性视频免费| 中文字幕另类日韩欧美亚洲嫩草| 国产精品人妻久久久影院| 精品久久久精品久久久| av线在线观看网站| 不卡视频在线观看欧美| av福利片在线| 国产精品秋霞免费鲁丝片| 国产成人精品婷婷| 97人妻天天添夜夜摸| 亚洲精品一区蜜桃| 在线天堂最新版资源| 国产探花极品一区二区| 日韩av在线免费看完整版不卡| 中文字幕色久视频| 99久久精品国产国产毛片| 秋霞伦理黄片| 免费观看av网站的网址| 亚洲成国产人片在线观看| 亚洲人成77777在线视频| 黄片无遮挡物在线观看| 国产精品 欧美亚洲| 我要看黄色一级片免费的| 国产乱人偷精品视频| 国产又色又爽无遮挡免| 叶爱在线成人免费视频播放| 91国产中文字幕| 亚洲精品国产色婷婷电影| 少妇 在线观看| 亚洲三级黄色毛片| 免费高清在线观看视频在线观看| 久久久久久久久久久久大奶| 日韩中字成人| 啦啦啦在线免费观看视频4| 热99国产精品久久久久久7| 精品99又大又爽又粗少妇毛片| 好男人视频免费观看在线| 黄网站色视频无遮挡免费观看| 少妇被粗大的猛进出69影院| 精品亚洲乱码少妇综合久久| 男人操女人黄网站| 国产成人aa在线观看| 免费在线观看完整版高清| 色吧在线观看| 欧美日韩亚洲高清精品| 亚洲激情五月婷婷啪啪| 日本欧美国产在线视频| 国产精品.久久久| 午夜激情久久久久久久| 国产精品嫩草影院av在线观看| 久久国产精品男人的天堂亚洲| 国产av码专区亚洲av| 1024香蕉在线观看| 欧美精品亚洲一区二区| 成人影院久久| 精品少妇久久久久久888优播| 欧美人与性动交α欧美精品济南到 | 欧美精品亚洲一区二区| 麻豆乱淫一区二区| 色婷婷久久久亚洲欧美| 免费av中文字幕在线| 亚洲欧洲日产国产| 大片免费播放器 马上看| 国产精品99久久99久久久不卡 | 91午夜精品亚洲一区二区三区| 香蕉国产在线看| 中文字幕色久视频| 久久影院123| 黄色毛片三级朝国网站| 不卡av一区二区三区| 九色亚洲精品在线播放| 三上悠亚av全集在线观看| 国产精品香港三级国产av潘金莲 | 最近最新中文字幕免费大全7| 丰满乱子伦码专区| 欧美少妇被猛烈插入视频| 久久99热这里只频精品6学生| 午夜日韩欧美国产| 国产 精品1| 久久久久久久国产电影| 天天躁夜夜躁狠狠躁躁| 久久久欧美国产精品| 久久av网站| 国产免费福利视频在线观看| 亚洲美女黄色视频免费看| 久久这里只有精品19| 国产综合精华液| 日韩一区二区三区影片| 精品亚洲成a人片在线观看| 亚洲成人av在线免费| 2022亚洲国产成人精品| 高清av免费在线| 亚洲国产色片| 亚洲精品一区蜜桃| 亚洲精品久久久久久婷婷小说| 看免费av毛片| videos熟女内射| 久久精品国产亚洲av高清一级| 亚洲av成人精品一二三区| 一区二区三区乱码不卡18| 久久精品熟女亚洲av麻豆精品| 成人毛片a级毛片在线播放| 91久久精品国产一区二区三区| 麻豆乱淫一区二区| 1024视频免费在线观看| 成年人午夜在线观看视频| 午夜福利在线免费观看网站| 久久韩国三级中文字幕| 亚洲精品美女久久久久99蜜臀 | 一级,二级,三级黄色视频| 9191精品国产免费久久| 免费在线观看完整版高清| 亚洲精品国产一区二区精华液| 97在线人人人人妻| 99香蕉大伊视频| 狂野欧美激情性bbbbbb| 成人毛片60女人毛片免费| av卡一久久| 亚洲男人天堂网一区| 欧美国产精品一级二级三级| 日韩一区二区视频免费看| 精品一区二区免费观看| 男人舔女人的私密视频| 国产极品粉嫩免费观看在线| 在线亚洲精品国产二区图片欧美| 美女主播在线视频| 成人二区视频| 成人黄色视频免费在线看| 狂野欧美激情性bbbbbb| 久久精品国产亚洲av高清一级| videossex国产| 日韩,欧美,国产一区二区三区| av天堂久久9| 亚洲精品久久成人aⅴ小说| 女人精品久久久久毛片| 九九爱精品视频在线观看| 久久精品亚洲av国产电影网| 国产精品免费视频内射| 18禁观看日本| 国产成人一区二区在线| 永久网站在线| 啦啦啦在线免费观看视频4| 黄色视频在线播放观看不卡| 欧美av亚洲av综合av国产av | 天天躁狠狠躁夜夜躁狠狠躁| 电影成人av| av在线播放精品| 亚洲欧美色中文字幕在线| 亚洲国产av新网站| 久久人人97超碰香蕉20202| 久久精品亚洲av国产电影网| 另类精品久久| 青草久久国产| 国产一级毛片在线| 肉色欧美久久久久久久蜜桃| 9色porny在线观看| tube8黄色片| 中文字幕另类日韩欧美亚洲嫩草| 国产精品免费大片| 国产欧美日韩一区二区三区在线| 肉色欧美久久久久久久蜜桃| 青青草视频在线视频观看| 老女人水多毛片| 国产在线免费精品| 一本—道久久a久久精品蜜桃钙片| 亚洲图色成人| 波多野结衣av一区二区av| 91国产中文字幕| av电影中文网址| 80岁老熟妇乱子伦牲交| 美女国产视频在线观看| 日韩一区二区三区影片| 日韩一本色道免费dvd| 91精品伊人久久大香线蕉| 韩国精品一区二区三区| 电影成人av| 国产精品久久久久久av不卡| 精品国产一区二区三区四区第35| 亚洲av成人精品一二三区| 爱豆传媒免费全集在线观看| 精品少妇一区二区三区视频日本电影 | av在线老鸭窝| 中文天堂在线官网| 免费观看性生交大片5| 亚洲精品国产av蜜桃| 国产精品一区二区在线观看99| 久久久久精品久久久久真实原创| 国产精品免费视频内射| 丝袜美腿诱惑在线| 免费高清在线观看视频在线观看| 国产成人精品婷婷| 十八禁高潮呻吟视频| 国产极品天堂在线| 夜夜骑夜夜射夜夜干| 美女午夜性视频免费| 精品久久蜜臀av无| 纯流量卡能插随身wifi吗| 高清黄色对白视频在线免费看| 国产人伦9x9x在线观看 | 91午夜精品亚洲一区二区三区| 老司机亚洲免费影院| 亚洲久久久国产精品| 国产激情久久老熟女| 欧美 日韩 精品 国产| 亚洲欧美一区二区三区黑人 | 免费大片黄手机在线观看| 中文精品一卡2卡3卡4更新| 久久久久久人人人人人| freevideosex欧美| 一本大道久久a久久精品| 日韩一区二区视频免费看| 99热网站在线观看| 国产一级毛片在线| 国产黄色免费在线视频| 婷婷色综合www| 少妇精品久久久久久久| 国产黄频视频在线观看| 搡老乐熟女国产| 欧美精品亚洲一区二区| 熟女av电影| 桃花免费在线播放| 国产乱来视频区| 女性生殖器流出的白浆| 午夜久久久在线观看| 精品一品国产午夜福利视频| 中国三级夫妇交换| 中文字幕制服av| 伦精品一区二区三区| 极品人妻少妇av视频| 啦啦啦啦在线视频资源| 边亲边吃奶的免费视频| 9热在线视频观看99| 在线天堂中文资源库| 欧美另类一区| 天堂俺去俺来也www色官网| 99精国产麻豆久久婷婷| 9色porny在线观看| 亚洲精品久久成人aⅴ小说| 少妇人妻久久综合中文| 国产精品久久久久久久久免| 男女啪啪激烈高潮av片| 一边摸一边做爽爽视频免费| 欧美精品人与动牲交sv欧美| 一级毛片我不卡| 亚洲在久久综合| 99热网站在线观看| 欧美97在线视频| 成人国产麻豆网| 国产日韩欧美视频二区| 伊人亚洲综合成人网| 精品少妇内射三级| 韩国av在线不卡| 黄色配什么色好看| 国产激情久久老熟女| 亚洲国产毛片av蜜桃av| 久久精品国产自在天天线| 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| videos熟女内射| 一边亲一边摸免费视频| 亚洲视频免费观看视频| 国产精品一区二区在线不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄频高清免费视频| 久久久久久久精品精品| 新久久久久国产一级毛片| 国产精品国产三级专区第一集| 超碰成人久久| 最近中文字幕2019免费版| 日韩成人av中文字幕在线观看| 国产成人精品福利久久| av不卡在线播放| 视频在线观看一区二区三区| 中文字幕亚洲精品专区| 久久久久视频综合| 国产av一区二区精品久久| 人人澡人人妻人| 久久久国产精品麻豆| 性色avwww在线观看| 人成视频在线观看免费观看| 国产精品嫩草影院av在线观看| 亚洲精品日韩在线中文字幕| 亚洲国产成人一精品久久久| 国产一区亚洲一区在线观看| 亚洲精品一二三| 亚洲精品av麻豆狂野| 国产精品av久久久久免费| 宅男免费午夜| 边亲边吃奶的免费视频| 午夜福利视频精品| 晚上一个人看的免费电影| 亚洲第一区二区三区不卡| 亚洲欧美精品综合一区二区三区 | 搡老乐熟女国产| 国产在线视频一区二区| 国产成人免费观看mmmm| 曰老女人黄片| 成人手机av| 久热久热在线精品观看| 精品一区二区三区四区五区乱码 | 亚洲精品美女久久久久99蜜臀 | 国产极品粉嫩免费观看在线| 制服人妻中文乱码| 巨乳人妻的诱惑在线观看| av片东京热男人的天堂| 亚洲国产最新在线播放| 国产免费现黄频在线看| 精品卡一卡二卡四卡免费| 久久精品国产a三级三级三级| 观看av在线不卡| 欧美成人精品欧美一级黄| 日韩人妻精品一区2区三区| 中文字幕精品免费在线观看视频| 精品国产一区二区三区久久久樱花| 制服人妻中文乱码| 欧美黄色片欧美黄色片| 亚洲av男天堂| 少妇精品久久久久久久| av片东京热男人的天堂| 久久久a久久爽久久v久久| 亚洲欧洲国产日韩| 黑人巨大精品欧美一区二区蜜桃| 99久久中文字幕三级久久日本| 天天躁狠狠躁夜夜躁狠狠躁| 一本色道久久久久久精品综合| 国产成人欧美| 如日韩欧美国产精品一区二区三区| 亚洲精品第二区| 美女视频免费永久观看网站| 久久久久国产网址| 精品久久久久久电影网| 黄色毛片三级朝国网站| 日韩,欧美,国产一区二区三区| 看免费av毛片| 国产一区二区三区综合在线观看| 久久综合国产亚洲精品| 久久久精品免费免费高清| 日韩中文字幕视频在线看片| 精品亚洲乱码少妇综合久久| 国产欧美日韩综合在线一区二区| 香蕉精品网在线| 日韩av在线免费看完整版不卡| 水蜜桃什么品种好| 欧美日韩一区二区视频在线观看视频在线| 女性被躁到高潮视频| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 最近最新中文字幕大全免费视频 | 亚洲国产精品一区二区三区在线| 不卡av一区二区三区| 成人毛片a级毛片在线播放| 黄色视频在线播放观看不卡| 免费高清在线观看视频在线观看| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 超碰成人久久| 亚洲一码二码三码区别大吗| 国产黄色视频一区二区在线观看| 亚洲一码二码三码区别大吗| 国产精品久久久久久精品电影小说| 国产成人精品久久久久久| 久久午夜福利片| 久久久国产精品麻豆| 午夜福利视频精品| av在线老鸭窝| 精品午夜福利在线看| 国产精品麻豆人妻色哟哟久久| 日本wwww免费看| 久久精品国产a三级三级三级| av片东京热男人的天堂| 男女午夜视频在线观看| 黄片无遮挡物在线观看| 亚洲四区av| 黑丝袜美女国产一区| 啦啦啦视频在线资源免费观看| 在线精品无人区一区二区三| 精品少妇久久久久久888优播| 999久久久国产精品视频| 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| kizo精华| 一级毛片 在线播放| 国产亚洲精品第一综合不卡| 搡女人真爽免费视频火全软件| 久久精品国产鲁丝片午夜精品| 色吧在线观看| 一个人免费看片子| 免费不卡的大黄色大毛片视频在线观看| 青春草视频在线免费观看| 国产激情久久老熟女| 少妇的逼水好多| 一级毛片电影观看| 国产淫语在线视频| 极品人妻少妇av视频| 一个人免费看片子| 亚洲国产av新网站| 亚洲图色成人| 久久久久国产一级毛片高清牌| h视频一区二区三区| 亚洲精品国产色婷婷电影| h视频一区二区三区| 精品一区在线观看国产| 国产毛片在线视频| 国产成人av激情在线播放| 看免费成人av毛片| 久久久久久久久免费视频了| 日本黄色日本黄色录像| 丝袜美足系列| 国产高清不卡午夜福利| 99久国产av精品国产电影| 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 制服诱惑二区| 韩国精品一区二区三区| 亚洲一级一片aⅴ在线观看| 日韩 亚洲 欧美在线| 香蕉国产在线看| 97在线人人人人妻| 有码 亚洲区| 国产免费又黄又爽又色| 中国三级夫妇交换| 在线观看人妻少妇| 日韩大片免费观看网站| 亚洲第一av免费看| 色哟哟·www| 1024香蕉在线观看| 欧美日韩视频高清一区二区三区二| 久久av网站| 啦啦啦视频在线资源免费观看| 在线看a的网站| 三级国产精品片| 中文欧美无线码| 国产精品一二三区在线看| 电影成人av| 国产精品熟女久久久久浪| 精品国产一区二区三区久久久樱花| 中文字幕av电影在线播放| 亚洲男人天堂网一区| 久久午夜综合久久蜜桃| 亚洲男人天堂网一区| 男男h啪啪无遮挡| 在线看a的网站| 超碰97精品在线观看| 亚洲国产av影院在线观看|