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

    Compact voice coil deformable mirror with high wavefront fitting precision

    2023-12-01 05:50:14HULifaJIANGLvHUQiliXUXingyuHUANGYangWUJingjingYULin
    中國光學(xué) 2023年6期

    HU Li-fa,JIANG Lv,HU Qi-li,XU Xing-yu,HUANG Yang ,WU Jing-jing,YU Lin

    (1.School of Science, Jiangnan University, Wuxi 214122, China;2.Key Laboratory of Electro-Optical Countermeasures Test &Evaluation Technology, Luoyang 471003, China;3.Jiangsu Provincial Research Center of Light Industrial Opto-electronic Engineering and Technology, Wuxi 214122, China)

    Abstract: To meet the requirements of wavefront distortion correction for miniaturized adaptive optics systems,a Deformable Mirror (DM) using micro voice coil actuators was designed based on systematic theoretical analysis.The structural parameters of the micro voice coil actuator were optimized by electromagnetic theory and the finite element method.The DM was optimized with respect to thermal deformation,resonance frequency,coupling coefficient and other parameters.Finally,wavefront fitting and residual calculation were completed according to the influence function.The optimized 69-element Voice Coil Deformable Mirror (VCDM) has a large phase stroke,good thermal stability,and a large first resonance of 2 220 Hz.The RMS of the fitting residuals of the VCDM for the first 35 Zernike modes with a PV value of 1 μm are all below 30 nm.For complex random aberrations,the compact VCDM can reduce the wavefront RMS to less than 10%.Compared with a traditional VCDMs,the results of our compact VCDM indicate that it has a higher wavefront fitting precision.The compact VCDM with high performance and low cost has good potential applications in human retinal or airborne imaging systems.

    Key words: adaptive optics;deformable mirror;voice coil actuator;multiparameter analysis

    1 Introduction

    In adaptive optical imaging systems,the Deformable Mirror (DM) is usually used to correct aberrations induced by atmospheric turbulence[1].Different kinds of DMs have been investigated,such as the piezoelectric (PZT) DM[2],magnetic fluid DM[3],electrostatic DM[4],and voice coil DM (VCDM)[5].Compared with others,the VCDM has many merits,such as linear responses,a large phase stroke,and low voltage,which is suitable for the next generation of DMs with research spanning back to the 1990 s[6].

    In 1993,P.Salinari demonstrated a VCDM as a deformable secondary mirror of the Large Binoculars Telescope (LBT).After that,a prototype of P36 was manufactured and demonstrated based on the principle.The MMT336 DM was fabricated in 2000 and installed on the Multi-mirror Telescope(MMT) in 2002[7].MMT was the first ground-based optical telescope equipped with a VCDM as the adaptive secondary mirror,followed by the Large Binoculars Telescope,the Magellan Telescope,the Very Large Telescope,and others[8].However,the diameter of a Voice Coil Actuator (VCA) in these DMs ranges between 20 mm and 30 mm,which leads to a large VCDM with an aperture of about 1 000 mm.Large VCAs also require cooling systems.In addition,its coupling coefficient is generally larger than 60%,which is much larger than that of PZT-based DM.For a compact VCDM,a relatively small coupling coefficient is valuable to obtain high wavefront correction precision.Finally,it should be noted that it is expensive to use custom VCDMs with complex structures.In applications,such as retinal imaging,biological microscopy,and optical communications,a compact and relatively inexpensive DM with high electro-optical performance is necessary.

    Several kinds of compact DMs have been investigated,such as the membrane DM[9],micro-electro-mechanical system (MEMS) DM[10],and magnetic DM[11].The membrane DM is composed of a flexible silicon membrane supported by an underlying array of electrostatic parallel plate actuators.Usually,its phase stoke is not large,and its driving voltage is relatively high at about several tens of volts[12].In addition,a dedicated electronic device is needed to linearize its actuation.Although the MEMS DM has a large number of actuators,its phase stroke is only about several micrometers[13],limiting its applications as a wavefront corrector in adaptive optics systems.The magnetic DM uses small planar coils on silicon.Therefore,to generate large enough force to push or pull in the thin mirror of a magnetic DM,a relatively large current is necessary,which leads to high power dissipation.VCDMs have the merits of low cost and linear response compared to PZT DM.However,standard VCDMs have two obvious drawbacks: (1) large coil diameters with a high cost and relatively low bandwidth,(2) high coupling coefficient and a large current that cannot be used in a compact VCDM.Multiple compact VCDMs from ALPAO have a pitch of 0.8 mm,1.5 mm,2.5 mm,and 5 mm,and the number of actuators ranges from 69 to 3 228.They have a non-linearity below 3%,and a hysteresis as low as 2%.However,these VCDMs use the MEMS manufacturing process resulting in high costs.In this paper,an alternative compact VCDM is presented by deriving theoretical equations that could guide a reasonable design.The designed micro VCA has a small diameter of about 2 mm,which is only about a tenth of a standard one and allows a compact VCDM with a high density of actuators.Different factors such as phase stroke,resonant frequency,and power dissipation were analyzed.The VCDM's structure has been simplified and optimized to realize a low coupling coefficient and good electro-optical performances.In fact,these parameters are influenced by the interactions among multi-physical fields,such as electromagnetic,thermal,stress and strain fields.For example,increasing the current in a voice coil can produce a large electromagnetic force and increase its phase stroke.However,it can also lead to a large thermal dissipation,and even distort the VCDM’s thin mirror.Therefore,theoretical mathematical equations were derived to find reasonable requirements from applications.Further,numerical simulations based on the finite element method were used to optimize the compact VCDM.

    In the paper,the compact VCDM with high electro-optical performance is theoretically designed and demonstrated.In section 2,a theoretical analysis is given according to their typical applications.In section 3,its structure is designed,and its micro VCA is optimized based on electromagnetic finite element theory.In section 4,the temporal and spatial properties of the compact VCDM using optimized micro VCA are discussed.Finally,conclusions are given in section 5.

    2 Theory of the compact VCDM

    There are many factors that affect the electrooptical performance of a compact VCDM,such as phase stroke,the electromagnetic force of the micro VCA,the power dissipation and efficiency of the micro VCA,the response time of the VCDM and the coupling coefficient.A DM with only one or two of them with excellent values as reported in the references could not really lead to a high electro-optical performance for a compact VCDM[5].Therefore,it is necessary to systematically analyze the requirements from the applications according to electromagnetic theory.

    2.1 Requirement of phase stroke and frequency for the VCDM

    The phase stroke of the VCDM is determined by the stroke of its micro VCA.The phase stroke requirement of a DM could be different according to the media that the light beam propagates through,such as atmospheric turbulence or the human eye.For airborne/spaceborne adaptive optics applications,according to Noll’s theory,if the first J Zernike modes of the distortion have been corrected,the mean square residual error can be calculated as follows:

    where〈?2〉is the phase variance andajis the Zernike polynomial coefficient[14].WhenJis 1,Eq.(1) can be approximated as

    whereDis the telescope aperture andr0is the atmospheric turbulence coherent length generally ranging from 5 cm to 20 cm forλ=550 nm[15].For anr0of 5 cm andDof 8.4 m,Δ1is equal to 5 268 rad2according to Eq.(2).The wavefront that the DM has to correct is about 5 times larger than the residual[16].In addition,considering the reflection doubling effect,the required mechanical stroke of a VCA for atmospheric turbulence correction could be calculated as

    Unlike the aberrations of atmospheric turbulence,those in the human eye are difficult to represent using a theoretical model.Therefore,we take the approximate requirements for the phase stroke from experimental statistical results reported by other groups.Doble N[17]measured human eye aberrations on two different populations at the University of Rochester and Indiana University using a Shack-Hartmann wavefront sensor.The maximum PV of the wavefront is up to 53 μm for a pupil diameter ranging from 2 mm to 7.5 mm.Zhao J L[18]performed an aberration analysis on 332 healthy eyes and 344 diseased eyes.For a wavelength of 550 nm and a pupil diameter of 6 mm,the PV value of the total-order aberrations was (24.1±10.9) μm in the normal group and (19.0±10.05) μm in the abnormal group[18].Therefore,the mechanical phase stroke of a micro VCA should be at least 26.5 μm for human eye aberration correction.

    According to the discussions above,retinal imaging requires a larger phase stroke compared to astronomical observation.Considering the static aberration in the optical path,the compact VCDM needs a margin of phase modulation.The phase stroke of the compact VCDM should be at least 27 μm.

    The bandwidth of the VCDM should be high enough to meet the frequency of atmospheric turbulence aberrations[19].The Greenwood frequency is usually used to characterize atmospheric turbulence and it could be approximated as:

    wherevis the wind speed,andr0is the atmospheric turbulence coherence length.And the Greenwood frequency in an observatory station is necessarily larger than 50 Hz.

    2.2 Key parameters of micro VCA

    The electromagnetic force generated by the micro VCA directly determines the phase stroke of the VCDM[20].A micro VCA as shown in Fig.1.includes mainly two parts: one is the permanent magnet on the top and the other is the voice coil.When a current is passed through the voice coil,the magnet connected to the thin mirror of the DM is movable,which pushes or pulls the local thin mirror and generates deformation.

    Fig.1 Structure of micro VCA

    In Fig.1,rmis the magnet’s radius,hmis its height,dc-inis the coil’s inner diameter,dc-outis the coil’s outer diameter,hcis the coil’s height,andhgis the height of the air gap.When the permanent magnet in the VCA is in the magnetic field generated by the current-carrying coil,the electromagnetic force could be calculated as

    whereVis the volume of a permanent magnet,andBris the remanence of the permanent magnet,His field vector,Idlis the current element,anderis unit vector along a position vectorr.Eq.(5) connects the force to some important parameters such as the material and direction of a magnet,the current,the geometric size of the coil,and the gaphg.The finite element method is used to simulate the electric-magnetic field and calculate the force.

    To quantify analyze the performance of the micro VCA,efficiencyεis used as the evaluation factor[21].It could be calculated as

    where the unit ofεisFis the electromagnetic force,Iis the current through the coil,andRis the resistance of the coil.The larger the efficiency of the VCA,the better it is.

    3 Design and optimization of the compact VCDM

    The sectional view of the compact VCDM is shown in Fig.2(a).The structure of compact VCDM mainly includes ten components.The pupil diameter in the front shell is 17 mm.Fig.2(b) shows the arrangement of VCAs.Sixty-nine VCAs are arranged in a square array under the thin mirror.A single micro VCA contains a permanent magnet and a voice coil.One end of the strut is attached to a micro VCA,while the other is attached to the back of the thin mirror.The strut transmits the force of an actuator to the thin mirror and deforms the local surface of the thin mirror.The shaft of a strut is mounted with a spring in the spring-fixing plate.This structure provides the micro VCA with axial stiffness,which overcomes the problem where an actuator’s mover is directly suspended on the mirror in traditional VCDMs[22].The thin mirror is formed by spinning and depositing polyimide liquid on a highly smooth substrate.The assembly process of the DM can be found in reference [16].

    Fig.2 Structure of VCDM.(a) Sectional view;(b) arrangement of VCAs

    To realize the compact VCDM,the pitch of actuators must be small enough.The geometrical shape of the VCA array as shown in Fig.2 is square with a pitch of 2.5 mm.To ensure that there is an installation gap between the actuators,the diameter of a VCA should not be greater than 2.2 mm.At the same time,the mover of the VCA should be as small as possible to get a faster response,so the initial height of the magnet is set to 0.1 mm.The geometrical model of the VCAs is built and simulated with FEM.The initial geometrical sizes for the micro VCA are defined as follows: (1) a permanent magnet with a radius of 1 mm and a height of 0.1 mm;(2) a coil with an inner diameter of 0.2 mm,an outer diameter of 2.2 mm,and a height of 1 mm.To meet the phase stroke requirements calculated in Section 2.1 and reserve redundancy,the initial size of the air gap is set to 50 μm.The permanent magnet is made of NdFe35.The parameters that need to be optimized are listed in Tab.1.

    Tab.1 Parameters for micro VCAs to be optimized

    In our previous work[21],for axial and radial magnetization directions of the permanent magnet shown in Fig.3(a) and 3(c),the outputted force generated by radial magnetization must be larger than that of the axial magnetization.However,it is difficult to realize radial magnetization for the permanent magnet with so small a size.In addition,parallel magnetization as shown in Fig.3(b) is easy to achieve.

    Fig.3 Magnetization directions of the permanent magnet.(a) Axial magnetization,(b) parallel magnetization,(c) radial magnetization

    Fig.4 (color online) shows the comparison of electromagnetic force in three magnetization directions.From Fig.4 It can be seen that the effect of parallel magnetization is very good,and the electromagnetic force under the same current is close to that of radial magnetization.

    Fig.4 Electromagnetic force as a function of current for the permanent magnet at different magnetization directions

    Fig.5 (color online) shows the optimization results of the structural parameters under parallel magnetization and radial magnetization where parallel magnetization is represented as (P) radial magnetization as (R).

    Fig.5 Force and efficiency as a function of the VCA’s structural parameters: (a) permanent magnet radius;(b) permanent magnet height;(c) coil inner diameter;(d) coil outer diameter;(e) coil height;(f) air gap

    Fig.5(a) and Fig.5(b) show the electromagnetic force and efficiencyεas a function of magnet radius and height,respectively.The trade-off between efficiency and size is what determines the parameters.The larger the magnet radius and magnet height,the larger the electromagnetic force and efficiency.Their differences are the slopes of the electromagnetic force and efficiency curves.It should be noted that the permanent magnet is the determinant of the VCA,and increasing the permanent magnet’s volume will reduce the response speed of the VCA.Therefore,the optimized radius of the magnet is selected as 1 mm and its optimized height is selected as 0.2 mm.The results in Fig.5(c) indicate that as the inner diameter of the coil increases,both the electromagnetic force and the efficiency decrease.The optimized inner diameter of the coil is 0.2 mm and according to Fig.5(d),the optimized outer diameter of the coil is 2 mm.The change of force and efficiency in Fig.5(e) are slightly different from the others.Although the electromagnetic force increases with the height of the coil,its efficiency first increases and then decreases.This is because the wire being far away from the permanent magnet means that it contributes less to the electromagnetic force.Considering the efficiency of the actuator,the optimized coil height is 0.4 mm.Fig.5(f) shows that the larger the air gap,the smaller the electromagnetic force and efficiency.Therefore,the optimized air gap of the micro VCA is 50 μm.

    The optimized VCA has a diameter of 2 mm and a height of 0.65 mm.It is only one-tenth the size of the actuators for traditional VCDMs.Fig.6(color online) shows electromagnetic force and efficiency as a function of current before and after optimization.The obtained results indicate that the electromagnetic force and efficiency increased substantially after optimization.The output electromagnetic force of the VCA with optimized structural parameters increased by a factor of 10%.The results confirm that the efficiency also notably increased.The efficiency of the VCA before optimization is 0.016and it is up to 0.032after optimization.

    Fig.6 Electromagnetic force and efficiency as a function of current (Circle line,before optimization;pentagram line,after optimization)

    4 Discussions

    4.1 Aberration of the thin mirror due to the thermal effect

    To investigate the effect of temperature increases from the current in VCAs on the thin mirror,the model of a VCA is analyzed using Maxwell and steady-state thermal modules in Ansys software.The procedure of thermal analysis is shown in Fig.7(color online).

    Fig.7 The thermal analysis process of the VCDM

    In compact DMs,active cooling elements such as air cooling and water cooling are generally not used.The thermal convection of all parts in the VCDM is through the air.In the ANSYS steadystate thermal module,the air convection heat transfer coefficient is set to 5 W/(m2K)[23].According to the mechanical model of the VCDM,as shown in Fig.2,the internal heat generation of the voice coil is introduced by the ohmic loss calculated by the Maxwell module.After that,the thin mirror temperature obtained by the steady-state thermal module is transferred to the static structural module to calculate the thermal deformation of the thin mirror.

    The mirror material is CP1 Polyimide from NeXolve[24].O-rings,springs,and struts are made of 316 Stainless Steel.The adhesive between the struts and the mirror is Epoxy.The material of the permanent magnet is NdFe35,and the material of the coil is copper.To accelerate heat dissipation,the shell of the DM,spring fixing plate,and baseplate are all made of aluminum alloy.The parameters of the adopted materials are listed in Tab.2.

    Tab.2 Material parameters uesd in thermal analysis

    Fig.8(a) (color online) shows the temperature of the VCDM’s thin mirror as a function of currents.The maximum temperature of the mirror surface increases along with the current.When the current is 0.1 A corresponding to a current density of 10 A/mm2,the temperature rise of the VCDM is about 1 degree Celsius.Fig.8(b) (color online)shows the thermal deformation of the thin mirror due to non-uniform distribution of rising temperatures.It indicates that the PV of the thin mirror deformation increases along with the current.When the current density is 6 A/mm2corresponding to 0.06 A,the PV is 12.36 nm,which is about 0.022λforλ=550 nm.Therefore,the maximum control current of a single VCA should not be larger than 0.06 A for good thermal stability.At this time,the temperature and the deformation of the thin mirror are shown in the inserted figures of Fig.8(a) and Fig.8(b).

    Fig.8 Temperature and thermal deformation of the VCDM’s thin mirror as a function of current.(a) The temperature of the thin mirror as a function of current.The inset shows the temperature chart of the thin mirror when the current is 0.1 A.(b) The thermal deformation chart of the mirror surface as a function of current.The inset shows the deformation chart of the thin mirror when the current is 0.06 A

    4.2 First resonance frequency

    The first resonance frequency of the VCDM is mainly related to the mirror and the actuator’s stiffness.When the mirror material is selected,the mirror’s stiffness is determined by its thickness,and the stiffness of the VCA is dominated by the spring.As shown in Fig.9,the first resonance frequency is positively correlated with the spring’s stiffness for a given mirror thickness.However,for unchanging spring stiffnesses,the first resonance frequency increases and then decreases with an increasing mirror thickness,which shows that there is an extreme point in the relationship.The frame rate of the Adaptive Optics (AO) system should be 10 to 30 times higher than the bandwidth of the closed-loop system.Therefore,a mirror thickness ranging from 20 μm to 50 μm is preferable.

    Fig.9 Relationship between first resonance frequency of the DM and mirror thickness under different spring's stiffnesses

    4.3 Wavefront fitting precision

    The wavefront fitting precision of the VCDM is important for AO wavefront correction.The influence function of the VCDM could also be expressed as[25]:

    whereZ(r) is the thin mirror deformation at a distance ofr,kis the deformation of the mirror at the position of the central actuator,ωis the coupling coefficient that indicating the influence of the central actuator on neighboring actuators,the actuator spacingd0is 2.5 mm,andαis the Gaussian index.

    The coupling coefficientω,a key indicator of the DM,is mainly influenced by the thickness of the mirror and the stiffness of the actuator.There is an optimal range of coupling coefficients for a given DM;too large or too small will affect the performance of the system[26].As shown in Fig.10,the coupling coefficient increases with mirror thickness.In addition,the larger the stiffness,the smaller the coupling coefficient.

    Fig.10 Relationship between coupling coefficient of the DM and mirror thickness under different spring's stiffnesses

    The wavefront fitting errors as a function of the coupling coefficient are shown in Fig.11 for four Zernike modes Z1,Z3,Z6,and Z10 with a PV value of 8 μm.The results indicate that with an increase in the coupling coefficient,errors decrease rapidly at first and then increase slowly.The optimal coupling coefficient for the compact VCDM is about 25%,which is much smaller than that of a standard VCDM with a value larger than 60%.Generally,a relatively low coupling coefficient is preferable for high wavefront correction precision.A VCDM model named S1 with a mirror thickness of 50 μm and an actuator stiffness of 90 N/m was built.Its influence function is shown in Fig.12.Sixty-nine VCAs are arranged in a square array under the thin mirror.The pitch of S1 is 2.5 mm,the coupling coefficient is 24.9% and the first resonance frequency is 2 220 Hz.

    Fig.11 Correction ability varying with the coupling coefficient

    Fig.12 The influence function of S1

    The MATLAB codes were programmed to complete the wavefront fitting and residual calculation according to the influence functions.It should be pointed out that,to simplify the calculation model,the error caused by thermal deformation and response delay is not considered.In the fitting experiment,S1 is actuated to generate Zernike shapes with a PV value of 1 μm.Fig.13 shows the wavefront Root Mean Square (RMS) of the first 35 Zernike modes and the residual wavefront RMS after fitting.The results indicate that the fitting effect of S1 for the first 35 Zernike modes is very good and the RMS is always less than 30 nm.Compared with the large 121-element VCDM designed by Zhang Z G,the wavefront fitting precision of S1 is higher by 15%[27].Fig.14 (color online) shows the wavefront fitting results of S1 for a complex random aberration.The complex random aberration consists of 90 Zernike modes and its PV value is 1.647 μm.It can be seen that the fitting wavefront generated by S1 is very similar to the original wavefront.The VCDM S1 reduces the RMS of the complex random aberration from 0.272 2 μm to 0.023 9 μm.The above results show that S1 has high wavefront fitting precision.

    Fig.13 Wavefront fitting results of S1 for a single Zernike mode

    Fig.14 The wavefront fitting results of S1 for complex random aberrations.From left to right,there is the original wavefront,the fitting wavefront and the residual wavefront

    5 Conclusion

    The VCDM has many merits,such as being hysteresis-free,requiring low voltage,having a large phase stroke and so on.However,there are many limitations and factors in VCDMs that will affect their performance.In the paper,to design a compact VCDM with improved properties,an analytical expression based on electromagnetic theory is systematically built,which includes its parameters such as phase stroke,force,efficiency,and frequency.They guide the design and reasonable optimization of the compact VCDM.The thin mirror of the VCDM is made of polyimide and is coated with a thin layer of evaporated aluminum in order to make it optically reflective.The struts of micro VCAs are installed in the spring fixing plate,whose springs provide micro VCAs with axial stiffness.This structure makes it possible for the VCDM to have a higher operating bandwidth.Based on multiparameter analysis and finite element analysis using ANSYS,a well-optimized VCDM is obtained with a micro VCA of 2 mm in diameter and an efficiency of 0.032.The maximum current of a micro VCA is 0.06 A,which gives the VCDM excellent thermal stability.The temperature difference of the thin mirror is less than 0.4 °C at the maximum current of 0.06 A,and the thermal deformation of the thin mirror is only 12.36 nm.The first resonance frequency of the 69-element compact VCDM S1 is 2 220 Hz.High wavefront fitting precision with a relatively low coupling coefficient of about 25% is also demonstrated in this paper.The fitting results of Zernike aberrations show that the wavefront fitting precision of the compact VCDM increased by 15% compared with traditional VCDM.In addition,S1 can reduce the wavefront RMS of complex random aberrations to less than 10%.The above results indicate that the compact VCDM can satisfy the re-quirements of miniaturized adaptive optics systems.The design and optimization methods are also valu-able for the design of other kinds of DMs.Ourdesign decreases the development cost of VCDMsand obtains a compact version with high electro-op-tical performance.

    满18在线观看网站| 国产精品国产三级国产专区5o| 黄色视频,在线免费观看| 亚洲熟女毛片儿| 日本猛色少妇xxxxx猛交久久| 久久性视频一级片| 免费少妇av软件| 制服诱惑二区| 免费观看人在逋| 91精品三级在线观看| 欧美激情久久久久久爽电影 | a在线观看视频网站| 在线观看人妻少妇| 欧美精品亚洲一区二区| www.熟女人妻精品国产| 精品福利永久在线观看| 国产真人三级小视频在线观看| 天堂8中文在线网| 国产成人免费观看mmmm| 亚洲七黄色美女视频| 黄色 视频免费看| 精品视频人人做人人爽| 久热爱精品视频在线9| 亚洲免费av在线视频| 一区福利在线观看| 亚洲欧美精品综合一区二区三区| 大片免费播放器 马上看| 日本wwww免费看| 亚洲欧美精品自产自拍| 一区二区日韩欧美中文字幕| 日日爽夜夜爽网站| 日本vs欧美在线观看视频| 午夜福利视频在线观看免费| 亚洲精品乱久久久久久| 在线看a的网站| 女性被躁到高潮视频| 嫩草影视91久久| 欧美人与性动交α欧美精品济南到| 蜜桃在线观看..| av在线app专区| 91大片在线观看| 人成视频在线观看免费观看| 亚洲国产精品成人久久小说| 婷婷丁香在线五月| 国产精品 国内视频| 99久久99久久久精品蜜桃| 国产亚洲精品久久久久5区| 国产淫语在线视频| 欧美亚洲日本最大视频资源| 少妇 在线观看| 少妇猛男粗大的猛烈进出视频| e午夜精品久久久久久久| 亚洲精品国产区一区二| 国产一区二区三区av在线| 欧美日韩黄片免| 啦啦啦视频在线资源免费观看| 后天国语完整版免费观看| 最黄视频免费看| 另类精品久久| 欧美国产精品va在线观看不卡| 首页视频小说图片口味搜索| 肉色欧美久久久久久久蜜桃| 免费少妇av软件| 男女无遮挡免费网站观看| 午夜福利一区二区在线看| 亚洲 欧美一区二区三区| e午夜精品久久久久久久| 国产色视频综合| 亚洲人成电影观看| 女人久久www免费人成看片| 国产深夜福利视频在线观看| 日韩人妻精品一区2区三区| 下体分泌物呈黄色| 人人妻人人澡人人看| 日韩欧美一区二区三区在线观看 | www.精华液| 中文精品一卡2卡3卡4更新| 嫁个100分男人电影在线观看| 老熟妇乱子伦视频在线观看 | 国产视频一区二区在线看| 国产精品一二三区在线看| 99国产极品粉嫩在线观看| 日韩欧美一区二区三区在线观看 | 一二三四社区在线视频社区8| 考比视频在线观看| 日韩视频在线欧美| 巨乳人妻的诱惑在线观看| 精品人妻一区二区三区麻豆| 水蜜桃什么品种好| 91精品国产国语对白视频| 成人18禁高潮啪啪吃奶动态图| 性少妇av在线| 可以免费在线观看a视频的电影网站| 亚洲免费av在线视频| 国产精品久久久久久精品电影小说| 国产精品一区二区精品视频观看| 建设人人有责人人尽责人人享有的| 在线亚洲精品国产二区图片欧美| av天堂在线播放| av线在线观看网站| 一边摸一边做爽爽视频免费| 老司机在亚洲福利影院| 免费日韩欧美在线观看| 精品久久久精品久久久| 欧美日韩国产mv在线观看视频| 一本色道久久久久久精品综合| 日韩大片免费观看网站| 亚洲熟女精品中文字幕| 亚洲欧美成人综合另类久久久| 欧美成狂野欧美在线观看| tocl精华| 黄片大片在线免费观看| 性高湖久久久久久久久免费观看| 免费在线观看影片大全网站| 国产精品久久久久成人av| 99国产精品99久久久久| 伦理电影免费视频| 黑人猛操日本美女一级片| 色婷婷av一区二区三区视频| 欧美精品一区二区免费开放| 国产精品偷伦视频观看了| 美女主播在线视频| 美女大奶头黄色视频| 欧美亚洲 丝袜 人妻 在线| 久久久精品国产亚洲av高清涩受| 黄片大片在线免费观看| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品久久久久5区| 满18在线观看网站| 色94色欧美一区二区| 黄色视频,在线免费观看| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区在线观看99| 国产一区有黄有色的免费视频| 久热这里只有精品99| 两个人免费观看高清视频| 免费少妇av软件| www.自偷自拍.com| 999久久久精品免费观看国产| 欧美人与性动交α欧美软件| 一级毛片女人18水好多| 一级毛片女人18水好多| 中文精品一卡2卡3卡4更新| 亚洲成国产人片在线观看| av在线播放精品| 欧美人与性动交α欧美软件| 夜夜骑夜夜射夜夜干| 一区二区三区四区激情视频| 亚洲av成人一区二区三| 熟女少妇亚洲综合色aaa.| 成年女人毛片免费观看观看9 | 午夜免费成人在线视频| 欧美激情 高清一区二区三区| 男女下面插进去视频免费观看| 亚洲精品日韩在线中文字幕| 操出白浆在线播放| 亚洲国产欧美一区二区综合| 午夜福利在线观看吧| 亚洲美女黄色视频免费看| 午夜福利免费观看在线| 99精品欧美一区二区三区四区| 久久久精品免费免费高清| 亚洲国产毛片av蜜桃av| 久久精品人人爽人人爽视色| 黄片小视频在线播放| 黄片小视频在线播放| 美女高潮喷水抽搐中文字幕| 91成年电影在线观看| 日韩三级视频一区二区三区| 亚洲av电影在线进入| 51午夜福利影视在线观看| 亚洲美女黄色视频免费看| 精品国产乱码久久久久久男人| 19禁男女啪啪无遮挡网站| 岛国毛片在线播放| 麻豆乱淫一区二区| 久久人人爽人人片av| 视频区欧美日本亚洲| 无遮挡黄片免费观看| 狠狠精品人妻久久久久久综合| 热99国产精品久久久久久7| 淫妇啪啪啪对白视频 | 在线观看免费高清a一片| 黑人巨大精品欧美一区二区mp4| 手机成人av网站| 手机成人av网站| 水蜜桃什么品种好| 国产欧美日韩一区二区三区在线| av网站在线播放免费| 日韩制服骚丝袜av| 久久久精品94久久精品| 精品少妇内射三级| 在线十欧美十亚洲十日本专区| 精品一区在线观看国产| 亚洲av电影在线观看一区二区三区| 国产成人精品久久二区二区免费| 亚洲精品国产av成人精品| 国产亚洲午夜精品一区二区久久| 国产av又大| 女性生殖器流出的白浆| 亚洲 国产 在线| 欧美日韩精品网址| 热re99久久国产66热| 国产精品一区二区免费欧美 | 午夜老司机福利片| 国产av又大| 国产一区二区在线观看av| 久久久国产一区二区| 国产精品欧美亚洲77777| 午夜久久久在线观看| 欧美另类亚洲清纯唯美| 18禁裸乳无遮挡动漫免费视频| 亚洲欧洲日产国产| 高潮久久久久久久久久久不卡| 性色av一级| 美女视频免费永久观看网站| 亚洲成人手机| 黑人欧美特级aaaaaa片| 99久久99久久久精品蜜桃| 国产高清国产精品国产三级| 天堂8中文在线网| 久久99热这里只频精品6学生| 色老头精品视频在线观看| 天天影视国产精品| 考比视频在线观看| 18禁国产床啪视频网站| 国产亚洲精品一区二区www | 精品国产乱码久久久久久男人| 亚洲激情五月婷婷啪啪| 一本色道久久久久久精品综合| 精品一区二区三区四区五区乱码| 国产1区2区3区精品| 韩国高清视频一区二区三区| 亚洲男人天堂网一区| 51午夜福利影视在线观看| 精品久久久久久电影网| 成人影院久久| 中文字幕精品免费在线观看视频| 国产免费一区二区三区四区乱码| 王馨瑶露胸无遮挡在线观看| 动漫黄色视频在线观看| 久久国产精品男人的天堂亚洲| 人人妻人人澡人人爽人人夜夜| 久久人妻福利社区极品人妻图片| 欧美人与性动交α欧美精品济南到| 精品高清国产在线一区| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三 | 黄色片一级片一级黄色片| 丰满迷人的少妇在线观看| 精品少妇一区二区三区视频日本电影| 一区福利在线观看| 少妇裸体淫交视频免费看高清 | 18在线观看网站| 久久这里只有精品19| 青春草亚洲视频在线观看| 午夜福利在线免费观看网站| 91成人精品电影| 久久久久国产精品人妻一区二区| 韩国精品一区二区三区| 国产一区二区激情短视频 | 日韩视频一区二区在线观看| 叶爱在线成人免费视频播放| 丝袜人妻中文字幕| 欧美精品一区二区免费开放| 精品欧美一区二区三区在线| 久久午夜综合久久蜜桃| 嫁个100分男人电影在线观看| 中文欧美无线码| 叶爱在线成人免费视频播放| 国产精品av久久久久免费| 电影成人av| 国产精品亚洲av一区麻豆| 国产高清videossex| 精品久久蜜臀av无| 亚洲精品国产色婷婷电影| 国产av精品麻豆| 又大又爽又粗| 黄色片一级片一级黄色片| 老司机午夜福利在线观看视频 | 亚洲 欧美一区二区三区| 一级毛片女人18水好多| 久久精品国产综合久久久| 一区二区三区四区激情视频| 精品福利观看| 久久久久久久国产电影| 国产伦人伦偷精品视频| 久久国产亚洲av麻豆专区| 亚洲精品在线美女| 亚洲成av片中文字幕在线观看| 亚洲第一青青草原| 亚洲精华国产精华精| cao死你这个sao货| 麻豆国产av国片精品| 欧美亚洲 丝袜 人妻 在线| 亚洲国产日韩一区二区| 亚洲久久久国产精品| 欧美黑人精品巨大| 久久久久精品人妻al黑| 亚洲国产看品久久| www.自偷自拍.com| 欧美老熟妇乱子伦牲交| 亚洲国产欧美一区二区综合| svipshipincom国产片| 美女福利国产在线| 777久久人妻少妇嫩草av网站| 亚洲专区国产一区二区| 亚洲欧美色中文字幕在线| 在线观看人妻少妇| 蜜桃在线观看..| 在线观看人妻少妇| 午夜福利,免费看| 在线看a的网站| 国产成人影院久久av| 桃红色精品国产亚洲av| 99久久国产精品久久久| 老司机福利观看| 男男h啪啪无遮挡| 亚洲成人免费电影在线观看| tube8黄色片| 十八禁网站网址无遮挡| 欧美日韩中文字幕国产精品一区二区三区 | 波多野结衣一区麻豆| 久久精品久久久久久噜噜老黄| 欧美另类亚洲清纯唯美| 午夜91福利影院| 久久久久国产精品人妻一区二区| 黄色a级毛片大全视频| 亚洲av日韩在线播放| 黑人猛操日本美女一级片| 精品少妇内射三级| 国产精品久久久久久精品电影小说| 午夜影院在线不卡| 精品一品国产午夜福利视频| 久久精品国产亚洲av香蕉五月 | 亚洲伊人色综图| av在线app专区| 蜜桃在线观看..| 男女边摸边吃奶| 丰满人妻熟妇乱又伦精品不卡| 亚洲 国产 在线| 久久99一区二区三区| 午夜福利视频精品| 免费高清在线观看日韩| 免费在线观看日本一区| 国产在线观看jvid| 国产麻豆69| 超碰97精品在线观看| 男女免费视频国产| 女人精品久久久久毛片| 国产av又大| 亚洲av电影在线观看一区二区三区| 一本综合久久免费| 手机成人av网站| 亚洲欧美激情在线| 黄片小视频在线播放| 日韩熟女老妇一区二区性免费视频| 男人操女人黄网站| 叶爱在线成人免费视频播放| 国产黄频视频在线观看| 熟女少妇亚洲综合色aaa.| 大片免费播放器 马上看| 女人高潮潮喷娇喘18禁视频| 久久免费观看电影| 欧美国产精品va在线观看不卡| 美女午夜性视频免费| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 免费一级毛片在线播放高清视频 | 国产成人av教育| 日本猛色少妇xxxxx猛交久久| 亚洲欧美日韩高清在线视频 | 国产激情久久老熟女| 老司机午夜福利在线观看视频 | 精品视频人人做人人爽| 最近中文字幕2019免费版| 免费观看a级毛片全部| 精品国内亚洲2022精品成人 | 菩萨蛮人人尽说江南好唐韦庄| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 国产精品香港三级国产av潘金莲| 人成视频在线观看免费观看| 亚洲免费av在线视频| 精品高清国产在线一区| 99久久综合免费| 亚洲五月色婷婷综合| 精品人妻在线不人妻| 亚洲精品日韩在线中文字幕| 欧美久久黑人一区二区| 婷婷色av中文字幕| 国产免费福利视频在线观看| 欧美中文综合在线视频| 黄网站色视频无遮挡免费观看| 高清视频免费观看一区二区| 少妇人妻久久综合中文| 午夜激情av网站| 午夜福利,免费看| 一区在线观看完整版| 精品国内亚洲2022精品成人 | 另类亚洲欧美激情| 欧美人与性动交α欧美精品济南到| 亚洲自偷自拍图片 自拍| 另类亚洲欧美激情| 国产97色在线日韩免费| 悠悠久久av| 又大又爽又粗| 中文字幕精品免费在线观看视频| 两个人免费观看高清视频| 91国产中文字幕| 美国免费a级毛片| 国产精品香港三级国产av潘金莲| 一区二区三区精品91| 成年女人毛片免费观看观看9 | 一区二区av电影网| 嫩草影视91久久| 欧美人与性动交α欧美精品济南到| 一区二区三区激情视频| 久久免费观看电影| 9热在线视频观看99| 国产亚洲欧美在线一区二区| www.精华液| 久久久精品区二区三区| 老司机深夜福利视频在线观看 | 男人操女人黄网站| 免费一级毛片在线播放高清视频 | 国产av又大| 美女高潮到喷水免费观看| 一进一出抽搐动态| 国产成人a∨麻豆精品| videos熟女内射| 欧美精品人与动牲交sv欧美| 美国免费a级毛片| 国产在线免费精品| 亚洲国产中文字幕在线视频| 80岁老熟妇乱子伦牲交| 夫妻午夜视频| av天堂在线播放| 免费看十八禁软件| 中文字幕av电影在线播放| 久久久精品区二区三区| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 精品国产乱子伦一区二区三区 | 在线观看一区二区三区激情| av超薄肉色丝袜交足视频| 国产高清videossex| 久久久久久人人人人人| videosex国产| 最新在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲天堂av无毛| 国产精品国产三级国产专区5o| 久9热在线精品视频| 老司机靠b影院| 日本av免费视频播放| 亚洲欧美清纯卡通| 18禁国产床啪视频网站| 久久久精品国产亚洲av高清涩受| 老司机深夜福利视频在线观看 | 丰满饥渴人妻一区二区三| 欧美黑人精品巨大| 亚洲国产欧美一区二区综合| 精品第一国产精品| 欧美成人午夜精品| 日韩,欧美,国产一区二区三区| 午夜精品久久久久久毛片777| 亚洲精品第二区| 亚洲精品日韩在线中文字幕| 动漫黄色视频在线观看| 99精品久久久久人妻精品| 不卡一级毛片| 少妇人妻久久综合中文| 亚洲精品在线美女| 久久99一区二区三区| 黄色a级毛片大全视频| 国产成人欧美| 中文字幕另类日韩欧美亚洲嫩草| 亚洲免费av在线视频| 欧美性长视频在线观看| 国产不卡av网站在线观看| 亚洲第一欧美日韩一区二区三区 | 亚洲av美国av| 别揉我奶头~嗯~啊~动态视频 | 男女床上黄色一级片免费看| 国产黄频视频在线观看| 一区二区三区乱码不卡18| 国产免费av片在线观看野外av| 国产麻豆69| 欧美亚洲日本最大视频资源| 视频区图区小说| 97在线人人人人妻| 色综合欧美亚洲国产小说| 亚洲国产成人一精品久久久| 久久九九热精品免费| 午夜视频精品福利| 黄色片一级片一级黄色片| 精品人妻1区二区| 亚洲欧洲日产国产| 亚洲av美国av| 国产免费视频播放在线视频| 国产在线视频一区二区| 淫妇啪啪啪对白视频 | 三级毛片av免费| 亚洲少妇的诱惑av| 国产精品亚洲av一区麻豆| 宅男免费午夜| 久久99热这里只频精品6学生| 成人免费观看视频高清| 99国产精品99久久久久| 国产男女超爽视频在线观看| 我要看黄色一级片免费的| 性色av乱码一区二区三区2| 老司机亚洲免费影院| 如日韩欧美国产精品一区二区三区| 男人舔女人的私密视频| 在线av久久热| 女人被躁到高潮嗷嗷叫费观| 中亚洲国语对白在线视频| 又黄又粗又硬又大视频| 久久毛片免费看一区二区三区| 亚洲国产av影院在线观看| 色婷婷久久久亚洲欧美| 男女床上黄色一级片免费看| 桃花免费在线播放| 脱女人内裤的视频| 国产精品久久久av美女十八| 麻豆国产av国片精品| 一级黄色大片毛片| 最近最新中文字幕大全免费视频| av免费在线观看网站| 日韩人妻精品一区2区三区| 女人爽到高潮嗷嗷叫在线视频| 欧美激情 高清一区二区三区| 2018国产大陆天天弄谢| 老司机靠b影院| 中文字幕另类日韩欧美亚洲嫩草| 国产成人av激情在线播放| 亚洲国产av影院在线观看| xxxhd国产人妻xxx| 免费观看a级毛片全部| 又黄又粗又硬又大视频| 黄色视频不卡| 色综合欧美亚洲国产小说| 天天躁夜夜躁狠狠躁躁| cao死你这个sao货| 国产1区2区3区精品| 深夜精品福利| 久久久久久久精品精品| 国产又色又爽无遮挡免| 国产男人的电影天堂91| 老熟妇仑乱视频hdxx| 久久久久久亚洲精品国产蜜桃av| 欧美午夜高清在线| 咕卡用的链子| 亚洲五月婷婷丁香| 涩涩av久久男人的天堂| 淫妇啪啪啪对白视频 | 免费高清在线观看视频在线观看| 丰满少妇做爰视频| 中文字幕另类日韩欧美亚洲嫩草| 老司机影院成人| 90打野战视频偷拍视频| 国产在视频线精品| 51午夜福利影视在线观看| 精品高清国产在线一区| 国产精品.久久久| 精品高清国产在线一区| 高清视频免费观看一区二区| 国产1区2区3区精品| 在线天堂中文资源库| 窝窝影院91人妻| 人妻 亚洲 视频| 伦理电影免费视频| 热99久久久久精品小说推荐| 黄片播放在线免费| 大香蕉久久网| 不卡av一区二区三区| 麻豆国产av国片精品| 青春草视频在线免费观看| 美女福利国产在线| 超色免费av| 久久国产精品男人的天堂亚洲| 丁香六月欧美| 色精品久久人妻99蜜桃| 91成人精品电影| xxxhd国产人妻xxx| 亚洲国产欧美一区二区综合| 悠悠久久av| 国产97色在线日韩免费| 日韩视频在线欧美| 亚洲欧洲日产国产| 亚洲第一青青草原| 国产精品 欧美亚洲| 亚洲国产看品久久| 中文字幕色久视频| 亚洲国产看品久久| 国产激情久久老熟女| 18禁国产床啪视频网站| 人人澡人人妻人| 久久久水蜜桃国产精品网| 50天的宝宝边吃奶边哭怎么回事| a级片在线免费高清观看视频| 一级毛片女人18水好多| 美女脱内裤让男人舔精品视频| 99九九在线精品视频| 少妇裸体淫交视频免费看高清 | 日韩人妻精品一区2区三区| 精品少妇一区二区三区视频日本电影| 伊人亚洲综合成人网| 亚洲精品一卡2卡三卡4卡5卡 | 欧美中文综合在线视频| 美女高潮喷水抽搐中文字幕| 脱女人内裤的视频| 色94色欧美一区二区|