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

    Simulation and Analysis of Co-phasing Errors of the Segmented Primary Mirror Tiled by Petal-shaped Segments in LOT*

    2022-12-12 08:24:30SHENShidongCUIXiangqunZHANGYong
    天文學報 2022年6期

    SHEN Shi-dong CUI Xiang-qun ZHANG Yong

    (1 National Astronomical Observatories/Nanjing Institute of Astronomical Optics & Technology,Chinese Academy of Sciences,Nanjing 210042)(2 CAS Key Laboratory of Astronomical Optics & Technology,Nanjing Institute of Astronomical Optics & Technology,Nanjing 210042)(3 University of Chinese Academy of Sciences,Beijing 100049)

    ABSTRACT In order to fulfill the requirement for astronomical observation in the opticalinfrared band using large universal aperture optical telescopes,China is planning to construct a 12-meter large aperture optical-infrared telescope(LOT).LOT is of particular significance to Chinese astronomy and would be co-phased in the near-infrared waveband.The primary mirror of LOT will be constructed from either hexagonal segments or petal-shaped segments.Analysis of co-phasing errors using a segmented primary mirror tiled by petal-shaped segments for LOT has been presented.The proposed primary mirror of LOT has 72 petal-shaped segments and a 4.5-meter circular mirror in the center.Detailed theoretical derivation of co-phasing errors and simulation about out-plane errors separately and comprehensively have been provided.SR(Strehl Ratio)of LOT is greater than 0.8 when the RMS(Root Mean Square)value of normally distributed tip-tilt errors is less than 0.016′′or the RMS value of normally distributed piston error is less than 42.5 nm at the co-phasing wavelength 1 μm.The simulation results are of great importance for the segmented primary mirror with active optics in LOT.

    Key words telescopes,methods:analytical,methods:data analysis,techniques:high angular resolution

    1 Introduction

    With the rapid development of astronomical science and technology,and the urgent need for large telescopes with higher light collecting capacity and resolution,many large aperture opticalinfrared telescopes have been built or are under construction both in China and abroad.The main approaches for the design of large aperture telescopes include:(1)the thin mirror active optics approach used in the VLT(Very Large Telescope)of the European Southern Observatory[1];(2)the segmented mirror active optics approach used in the KECK(Keck telescope)of the United States[2];and(3)the combined thin mirror and segmented mirror active optics approach pioneered by LAMOST(Large Sky Area Multi-Object Fiber Spectroscopic Telescope)of China[3–5].However,due to mirror size limitations,the current aperture limit of a single mirror is 8.4 m,which was successfully constructed by the Mirror Laboratory at University of Arizona.In order to construct larger or giant aperture telescopes,segmented mirror active optics has become the approach of choice taking into consideration factors such as technical difficulty,cost,and construction time.

    To fulfill the requirements for astronomical observation in the optical-infrared band with large universal aperture optical telescopes,CFGT(Chinese Future Giant Telescope)was proposed in 2000,using a 30-meter primary mirror design with fan-shaped segments[6–11].Japan subsequently proposed JELT(Japan Extremely Large Telescope),the primary mirror of which would also have been 30 m in size[12].Japan has since built a middle aperture telescope,namely the Seimei telescope.The Seimei telescope is a 3.8-meter opticalinfrared telescope,the primary mirror of which consists of 18 petal-shaped segments[13–14].

    Chinese astronomers proposed LOT(12-meter large aperture optical-infrared telescope)in 2016 and have considered the use of a segmented primary mirror[15–17].However,there are two proposed segmentation schemes:one type is a segmented primary mirror tiled into 84 hexagonal segments which is commonly used in segmented mirror telescopes,and the other is an innovative design using a 4.5-meter circular center mirror and 72 petal-shaped segments.Fig.1 depicts two types of segmentation schemes.

    Fig.1 Two types of segmentation schemes

    The segmented primary mirror tiled into petal-shaped segments has numerous advantages compared to the segmented primary mirror tiled into hexagonal segments.For example,the primary mirror can have a circular profile,and the type of segments are consistent in every ring,regardless of physical size or profile.The number of petalshaped segments required is also fewer than would be required using hexagonal segments for the same size of primary mirror,even though the area of an individual hexagonal segment and a petal-shaped segment is almost the same[18–20].This paper will therefore focus on developing an innovative scheme as described in Table 1,in which the size parameters and quantity of each segment are introduced in detail.We can also consider adopting all fanshaped segments for LOT.However,changing the inner 4.5-meter diameter to fan-shaped segments will not affect the research described in this paper from a practical point of view.Furthermore,a perfect co-phasing wavefront in visible imaging can be obtained by using this 4.5-meter mirror;thin mirror active optics technology can also be used in the central mirror to achieve a combination of thin mirror active optics and fan-shaped segmented mirror active optics if necessary.

    Table 1 The primary mirror of LOT tiled with fan-shaped segments

    The development of LOT is of great significance to Chinese astronomy,and the segmented mirror active optics approach is the preferred choice for LOT taking into consideration co-focus in the visible band and co-phase in the infrared band[21–23].Theoretical derivation and simulation of co-phasing errors are presented in this paper when the primary mirror of LOT reaches co-phase at 1 μm[24–27].

    Error analysis and simulation are the prerequisites and necessary conditions for error allocation during the feasibility evaluation for any large aperture segmented telescope.Chanan et al.[28]got SR(Strehl Ratio)for KECK as a function of the RMS(Root Mean Square)segment co-phasing error.Besides,Troy et al.[29]analyzed the influence of gaps,secondary mirror obstruction and reflectivity of mirrors for TMT(Thirty Meter Telescope).In addition,Yaitskova et al.[30–32]derived the analytical expression for formulas describing highly segmented telescopes,established an optical calculation model of PSF(Point Spread Function)based on FFT(Fast Fourier Transform),and simulated the influence of piston and tip-tilt on SR of the telescope,providing data support for the construction of ELT(European Extremely Large Telescope).

    In order to evaluate our proposed scheme and study the sensitivity of co-phasing errors,in this study we undertake theoretical derivation and simulation analysis of co-phasing errors in LOT in which the primary mirror is tiled by petal-shaped segments.The results provide a reference and basis for the selection and construction of the primary mirror shape scheme of the LOT.

    2 Analysis of co-phasing errors

    Fig.2 introduces the global coordinate system of the primary mirror and the local coordinate system of the segment.TheXCYplane is the projection plane of the primary mirror in Fig.2,and theZ-axis is the direction of the optical axis perpendicular to theXCYplane.PointCis the origin of the global coordinate systemXY Z.The origin point,o,of the local coordinate system(xyz)of each segment is located at the center of circumcircle of the segment.They-direction is the radial direction of the primary mirror and thez-direction is the normal direction of a segment,wherexis perpendicular toyandz.

    Fig.2 XCY plane in the primary mirror global coordinate system and yoz plane in the segment local coordinate system

    Wis assumed to be the wavefront of the primary mirror,andΔWis the wavefront error of the primary mirror.ωiis the wavefront of theith segment,and(Δω)iis the wavefront error of theith segment.The formula to describe the approximate relationship between the Strehl Ratio and wavefront is as follows[33]:

    whereσwfis the RMS value of wavefront error(in wavelength units,λ)in Eq.(1).

    In addition,σwfcan be expressed as[34]:

    where

    Note that in Eq.(3),Arepresents the integral region of the primary mirror.

    In Eq.(4)and Eq.(5),Nrepresents the number of segments,and dzirepresents the integral region of the segment.

    The radius of the circumcircle of a petalshaped segment isa,the center of which is the center of its circumcircle.If the polar coordinate system(ρ,θ)is considered with regard to the local coordinate system whereρ=r/a,the surface of segment in the local coordinates can be expressed as a expansion[35]:

    αmnis the local coordinate expansion coefficient,bothmandnare integers andm≥n≥0 andm-n=even in Eq.(6).

    The radius of curvature of the primary mirror,Rc,is-38400 mm.prepresents an off-axis distance from the center of a segment to the center of the primary mirror,andp=(wherexoandyoare the coordinates of pointoin Fig.2).The conic constant of the primary,k,has a value of-0.9837843.α20is the coefficient of the first term in Eq.(6),represents defocus.The coefficient of the second term in Eq.(6)that represents astigmatism,α22,iska2p2/4.The coefficient of the third term in Eq.(6),representing coma,α31,iska3p/2.In this study,we primarily focus on the effects of astigmatism as this dominates in terms of wavefront aberrations[36].

    2.1 The influence of piston

    If there is some random piston errorαiwhich obeys normal distributionNof theith segment,the standard deviationσαis the RMS value of piston error,we can obtain the wavefront of theith segment:

    The wavefront of the segmented primary mirror can also be obtained:

    wherezirepresents the integral region of theith segment and segirepresents area of theith segment.By substituting this result in Eq.(2),it can be concluded that:

    At this time,the RMS value of piston error is 0.0375λwhen SR is 0.8 orσwfis 0.075λ.The RMS value of piston should be less than 37.5 nm when the co-phasing wavelength is 1 μm.

    2.2 The influence of tip-tilt

    We use the notationai,γi,βito represent the radius of circumcircle of theith petal-shaped segment,the amount of tilt in thex-direction of theith petal-shaped segment(inλunits),and the amount of tilt in they-direction of theith petal-shaped segment(inλunits)respectively.Ifγi~Nandβi~N,the standard deviationσγis the RMS value of tip and the standard deviationσβis the RMS value of tilt,the wavefront of theith petal-shaped segment can be obtained:

    Assuming that the standard deviation of the normal distribution of tip-tilt is the same,we can conclude that:

    The RMS value of tip-tilt errorσis 0.065λwhen SR is 0.8 orσwfis 0.075λ.The RMS value of tiptilt error should be less than 0.016′′(σwf/ai)when the co-phasing wavelength is 1 μm.The theoretical co-phasing errors of the two segmented mirror schemes using hexagonal segments[37]and petalshaped segments are compared in Table 2.

    Table 2 Comparison of RMS value ofco-phasing errors between two segment shapes

    2.3 The influence of in-plane errors

    Decenter inxaxis can be ignored among the three in-plane errors,based on the rotational symmetry of the quadric surface.Here,the radial error decenter inyaxis and clocking error are mainly analyzed.We assume that the radial error of radial distancepof theith fan-shaped segment can be expressed as(Δp)i(inλunits),and,σΔpis the RMS value of radial translation error.The radial distance of fanshaped segments in the first ring isp1,with a value of 2.919352 m.The radial distance of fan-shaped segments in the second ring isp2,4.160594 m and the radial distance of fan-shaped segments in the third ring isp3,5.404607 m,then it can be calculated as follows:

    wherepiis the radial distance of theith segment.

    In this case,theσwfvalue is 0.075λwhen SR is 0.8.The RMS value of radial translation errorσΔpis 6162.8λ.The RMS value of radial translation error of each fan-shaped segment should be less than 6.162 mm when the co-phasing wavelength is 1 μm.

    Assuming that the clocking error of theith fan-shaped segment can be expressed as(Δθ)i(radian units),and(Δθ)i~N,the wavefront error of theith fan-shaped segment can be obtained:

    The house on Beechwood Avenue had a big Front yard with shade trees and fruite trees. In the middle of the yard was a water pump where the four little girls pumped water for cooking, cleaning and watering the garden. On one side of the yard, Grandma and Grandpa planted tomatoes, beans, squash2(), cucumbers, peppers and strawberries to feed their growing family. They planted roses geraniums lilacs() and irises3 on the other side of the yard,around the statue of the Blessed Mother.

    In this case,theσwfvalue is 0.075λwhen SR is 0.8.The RMS value of radial translation errorσΔθis 2.93′when the co-phasing wavelength is 1 μm.

    3 Simulation results of cophasing errors

    As described above,a number of approximate calculations are required in the theoretical derivation of the co-phasing errors of the primary mirror which need to be implemented in the simulation software.The LOT model was established to simulate real co-phasing errors,which were used for the final reference and comparative verification of the simulation analysis.Nasmyth focal system of LOT will be built as a priority in order to meet the requirements of domestic astronomers for fine detail observation.The parameters are shown in Table 3.

    Table 3 Main parameters of Nasmyth focal system of LOT(Focal ratio F number/14.4)

    Fig.3 is a ring of petal-shaped segments in the aperture.The inner radius and outer radius of a segment are given byr,Rrespectively and the space between adjacent segments is 2din Fig.3.A1 and B1 represent two vertices of a segment respectively.

    Fig.3 Diagram of calculating the aperture of petal-shaped segment

    The four vertices of a segment are calculated,and a program was written to generate the segment aperture,which would then be imported into ZEMAX software for the next stage of the simulation.

    The fan-shaped segments in each ring can be solved by calculating values for one of the fanshaped segments in a given position,and then rotated to obtain the positions of the other segments in the same ring.For example,the coordinates A1 and B1 of a segment in the second ring are calculated according to

    and

    The segment aperture program is used to generate the corresponding aperture,which would be imported into ZEMAX software for the next simulation.

    The loss of light-gathering energy is 0.41% which is very small compared to the light loss caused by the secondary mirror and the central hole of the M4 mirror in LOT,assuming the space between the fan-shaped segments is 3 mm.It is further calculated that the energy loss values are 1.3% and 2.6% when the gap between fan-shaped segments is 10 mm and 20 mm,respectively.The SR value of LOT is 1 in both cases,and it can therefore be concluded that the gap between the segments has no effect on the image quality except for the shadow effect on the final diffraction pattern.

    Each petal-shaped segment has six degrees of freedom and different reflective surface properties,and all segments are independent in the nonsequence mode of ZEMAX software.This method can simulate the actual segmented mirror system as shown in Fig.4.Since three in-plane degrees of freedom of a segment,such as radial displacement and rotation about the center of a segment,can be guaranteed by the machining accuracy,the influence of the three out-of-plane degrees of freedom will be mainly simulated.

    Fig.4 3D layout of LOT

    3.1 Comprehensive simulation analysis of tip-tilt

    The Monte Carlo method is used to analyze the sensitivity of tip-tilt errors.The tilt angles of 72 segments in ZEMAX software obey the same normal distribution.The tolerance width is±2σ(whereσis the standard deviation),where the two ends of the width are equal to the maximum and minimum values of the tolerance operands,and RMS of tilt angles isσif wavefront error RMS is used as the error evaluation standard.If the tilt angles of all segments are uniformly distributed,the tolerance width is±ν,and both ends of the width are equal to the maximum and minimum values of the tolerance operands,the RMS value of tilt angles is given byν/

    Firstly,a comprehensive simulation analysis of tip-tilt is carried out.It is assumed that the two types of errors obey the same normal distribution with a tolerance range of±0.032′′.The simulation results are shown in Table 4,and the wavefront of the primary mirror is better than 0.0755λ(λis 1 μm)for 90% of the processes.The calculated SR of telescope will be larger than 0.8 when the tip-tilt RMS value is less than 0.016′′,which is consistent with the theoretically derived value.A simulation result is shown in Fig.5,in which the RMS value of tip-tilt is 0.015′′,the RMS value of the wavefront of LOT is 0.068λ(λis 1 μm)and the SR value of LOT is 0.84.

    Table 4 Simulation results of the influence of tip-tilt on SR(subject to normal distribution)

    Fig.5 Wavefront of LOT(RMS value of tip-tilt is 0.015′′)

    In addition,it is assumed that tip-tilt is uniformly distributed at the same time and the tolerance range is±0.025′′.The simulation results are as follows in Table 5.

    Table 5 Simulation results of the influence of tip-tilt on SR(subject to uniform distribution)

    The simulation results show that the wavefront of the telescope in 90% of the processes is better than 0.0758λ,indicating that the SR value of the telescope is greater than 0.8 when the RMS value of the uniformly distributed tip-tilt errors is less than 0.0145′′.One of the simulation results is shown in Fig.6,where the RMS value of tip-tilt is 0.013′′,the RMS value of the wavefront of LOT is 0.0588λ(co-phasing wavelengthλis 1 μm)and the SR of LOT is 0.847.

    Fig.6 Wavefront of LOT in a simulation of uniformly distributed tip-tilt(RMS value of tip-tilt is 0.013′′)

    3.2 Simulation analysis of piston

    Monte Carlo analysis for the piston of all segments is undertaken,as the piston of each segment can be assumed to follow a normal distribution.The average value of piston is 0,and the tolerance range is±85 nm.90% of processes will be better than 0.075λ(λis 1 μm),which is taken as the evaluation criterion,and simulation results are shown in Table 6.The RMS value of the wavefront of the telescope in 90% of the processes is less than 0.072λ,indicating that the SR of telescope is greater than 0.8 when RMS of piston is less than 42.5 nm.One of the simulation results is shown in Fig.7,in which the RMS value of piston is 30.5 nm,the RMS value of the wavefront of LOT is 0.0596λ(co-phasing wavelengthλis 1 μm)and the SR of LOT is 0.87.

    Fig.7 Wavefront of LOT(RMS value of piston is 30.5 nm)

    Table 6 Simulation results of the influence of piston on SR(subject to normal distribution)

    In addition,the simulation results indicate the RMS value of the wavefront in 90% of processes in Monte Carlo analysis is superior to 0.075λ,assuming that piston of segments are uniformly distributed with a tolerance range of±71 nm.Simulation results are shown in Table 7.The SR of the telescope is greater than 0.8 when the RMS value of uniformly distributed piston is less than 41 nm.One of the results of the simulation is shown in Fig.8,in which the RMS value of piston is 36.1 nm,the RMS value of the wavefront of LOT is 0.07λ(cophasing wavelengthλis 1 μm)and the SR of LOT is 0.822.

    Fig.8 Wavefront of LOT(RMS value of piston is 36.1 nm)

    3.3 Comprehensive simulation analysis of tip-tilt and piston

    Table 8 shows theoretical values of tip-tilt and piston that are normally distributed with the same weighting.Monte Carlo analysis of the out-ofplane errors is presented below,aiming for 90% of the simulation processes to be better than 0.075λ(λis 1 μm).The piston of all segments is normally distributed and has a range of±56 nm and an average value of 0.In addition,tip-tilt values of all segments follow a normal distribution,with a range of±0.024′′and an average value of 0.The simulation results are shown in Table 9.

    Table 7 Simulation results of the influence of piston on SR(subject to uniform distribution)

    Table 8 Theoretical values of tip-tilt and piston with equal weighting effect

    Table 9 Simulation results of out-of-plane errors(subject to normal distribution,equal weight)

    The simulation results show that the RMS value of the wavefront for 90% of processes is better than 0.07235λ.Therefore,the SR of the telescope is greater than 0.8 when the RMS value of tip-tilt with a normal distribution is less than 0.012′′and the RMS value of piston is less than 28 nm.The RMS value of the wavefront of the telescope is 0.06125λand SR=0.867,given one of the integrated simulation processes,seen in Fig.9,where the RMS value of tip-tilt is 0.009′′and the RMS value of piston is 23.4 nm.

    Fig.9 Wavefront of LOT(RMS value of tip-tilt is 0.009′′and RMS value of piston is 23.4 nm).

    The requirement for the tip-tilt error of segments can be relaxed and the RMS value of piston can be reduced to 20 nm or 15 nm in practical engineering.The comprehensive simulation results are shown in Table 10 and Table 11.

    RMS of the wavefront of LOT is better than 0.074λduring 90% of the processes in Table 10,indicating that the SR of the telescope is better than 0.8 when the RMS value of piston is less than 20 nm and the RMS value of tip-tilt is less than 0.014′′.The RMS of the wavefront of LOT is better than 0.0746λduring 90% of the processes in Table 11,indicating that the SR of the telescope is greater than 0.8 when the RMS of piston is less than 15 nm and the RMS of tip-tilt is less than 0.015′′.

    Table 10 Monte Carlo simulation results of out-of-plane errors(RMS value of piston is 20 nm and RMS value of tip-tilt is 0.014′′)

    Table 11 Monte Carlo simulation results of out-of-plane errors(RMS value of piston is 15 nm and RMS value of tip-tilt is 0.015′′)

    4 Conclusion

    The Chinese optical infrared astronomy community has strong potential to develop into an international-class community;however,it is challenging to grow this potential without a large telescope like LOT as a platform.Therefore,there is a strong need of a general-purpose optical telescope for China,which could be equipped with different scientific instruments.

    We study tiling schemes using fan-shaped segments,and consider the influence of various tiling errors from the segmented primary mirror on image quality,especially the Strehl Ratio,and complete the theoretical derivation for the whole segmented mirror active optics simulation procedure.The segment errors,including three out-plane segment errors,are simulated using both numerical and synthetic approaches,with a lot of simulation results successfully achieved.The SR of LOT is greater than 0.8 when the RMS value of normally distributed tip-tilt errors is less than 0.016′′or the RMS value of normally distributed piston error is less than 42.5 nm with the co-phasing wavelength 1 μm.The comprehensive simulation analysis of the three out-of-plane errors shows that the RMS value of tip-tilt should be less than 0.012′′and the RMS value of piston should be less than 28 nm when tip-tilt and piston are normally distributed with the same weighting.In addition,the RMS value of normally distributed tip-tilt should be less than 0.014′′when the RMS value of normally distributed piston is 20 nm.The RMS value of normally distributed tip-tilt should be less than 0.015′′when the RMS value of normally distributed piston is 15 nm.

    The fan-shaped segmented mirror scheme is similar to the hexagonal segmented mirror scheme in terms of co-phasing errors,but the fan-shaped segmented mirror scheme is a better choice based on the processing cost and time.The results can be used as a reference for LOT and will provide LOT with the right to choose the scheme.Our results have important implications for the future development of large diameter optical-infrared telescopes.

    国产精品.久久久| 亚洲一码二码三码区别大吗| 岛国毛片在线播放| videosex国产| 欧美日韩综合久久久久久| 欧美日本中文国产一区发布| 亚洲成av片中文字幕在线观看 | 青春草亚洲视频在线观看| 久久午夜综合久久蜜桃| 青青草视频在线视频观看| 日本av手机在线免费观看| 中文精品一卡2卡3卡4更新| 午夜久久久在线观看| av天堂久久9| 久久久久人妻精品一区果冻| 成人影院久久| videossex国产| 春色校园在线视频观看| 久久午夜福利片| 亚洲国产精品成人久久小说| 国产片特级美女逼逼视频| 美国免费a级毛片| 亚洲精品av麻豆狂野| 久久精品国产综合久久久 | 免费高清在线观看日韩| 永久网站在线| 大码成人一级视频| 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 肉色欧美久久久久久久蜜桃| 色吧在线观看| 在现免费观看毛片| 午夜免费观看性视频| 在线看a的网站| 亚洲精品乱久久久久久| 亚洲av免费高清在线观看| 三上悠亚av全集在线观看| 欧美激情国产日韩精品一区| av福利片在线| 精品卡一卡二卡四卡免费| 自线自在国产av| 黄色一级大片看看| 色婷婷av一区二区三区视频| 午夜91福利影院| 亚洲国产成人一精品久久久| 色94色欧美一区二区| 久久精品久久久久久噜噜老黄| 久久久久久久精品精品| 韩国高清视频一区二区三区| 日韩在线高清观看一区二区三区| 欧美日韩成人在线一区二区| 综合色丁香网| 国产精品蜜桃在线观看| 亚洲av电影在线进入| 国产麻豆69| 午夜福利视频精品| 乱码一卡2卡4卡精品| 色网站视频免费| 七月丁香在线播放| 在线 av 中文字幕| 热99久久久久精品小说推荐| av播播在线观看一区| 亚洲av电影在线观看一区二区三区| 最近2019中文字幕mv第一页| 人人妻人人爽人人添夜夜欢视频| 婷婷色综合大香蕉| 亚洲婷婷狠狠爱综合网| videos熟女内射| 久久女婷五月综合色啪小说| 亚洲av男天堂| 国产老妇伦熟女老妇高清| 国产精品久久久av美女十八| 国产精品 国内视频| 国产亚洲午夜精品一区二区久久| 亚洲欧美一区二区三区国产| 99国产精品免费福利视频| 国产精品一国产av| 久久午夜福利片| av在线app专区| 一级毛片黄色毛片免费观看视频| 啦啦啦在线观看免费高清www| 一区二区日韩欧美中文字幕 | 制服人妻中文乱码| 亚洲五月色婷婷综合| 亚洲丝袜综合中文字幕| 2022亚洲国产成人精品| 国产精品女同一区二区软件| 妹子高潮喷水视频| 欧美激情极品国产一区二区三区 | 日日啪夜夜爽| 成年动漫av网址| 欧美xxⅹ黑人| 如日韩欧美国产精品一区二区三区| 男女边吃奶边做爰视频| 熟女人妻精品中文字幕| a级片在线免费高清观看视频| 2022亚洲国产成人精品| 免费观看无遮挡的男女| 精品亚洲乱码少妇综合久久| 国产精品.久久久| 日本vs欧美在线观看视频| 9191精品国产免费久久| 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕| 久久韩国三级中文字幕| 在线观看人妻少妇| 国产福利在线免费观看视频| 两个人看的免费小视频| 黄色 视频免费看| videosex国产| 国产亚洲午夜精品一区二区久久| 晚上一个人看的免费电影| 国产成人精品福利久久| 各种免费的搞黄视频| 一级毛片我不卡| 在线观看一区二区三区激情| 成人毛片a级毛片在线播放| 欧美另类一区| 天天躁夜夜躁狠狠躁躁| 亚洲欧美一区二区三区黑人 | 国产老妇伦熟女老妇高清| 80岁老熟妇乱子伦牲交| 日韩伦理黄色片| 精品卡一卡二卡四卡免费| 亚洲国产毛片av蜜桃av| 午夜福利网站1000一区二区三区| 亚洲欧洲日产国产| 久久久久久人人人人人| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧洲精品一区二区精品久久久 | 久久毛片免费看一区二区三区| 国产黄色免费在线视频| 国产1区2区3区精品| 精品人妻一区二区三区麻豆| 丝袜人妻中文字幕| 国产成人a∨麻豆精品| 乱人伦中国视频| 免费观看在线日韩| 搡老乐熟女国产| 丝袜在线中文字幕| 99国产综合亚洲精品| 五月玫瑰六月丁香| 人妻系列 视频| 欧美日本中文国产一区发布| 99久久中文字幕三级久久日本| av女优亚洲男人天堂| 亚洲丝袜综合中文字幕| 美女大奶头黄色视频| 日本午夜av视频| 欧美精品国产亚洲| 亚洲高清免费不卡视频| 欧美成人午夜精品| 欧美精品高潮呻吟av久久| 丝袜喷水一区| 免费少妇av软件| 日韩,欧美,国产一区二区三区| 交换朋友夫妻互换小说| 永久网站在线| 免费观看a级毛片全部| 亚洲美女视频黄频| 亚洲av福利一区| 9191精品国产免费久久| 晚上一个人看的免费电影| a级片在线免费高清观看视频| 美女大奶头黄色视频| 成人国语在线视频| 18禁在线无遮挡免费观看视频| 18禁观看日本| 久久久久久久久久成人| 欧美精品一区二区免费开放| av天堂久久9| 亚洲国产精品成人久久小说| 在线观看人妻少妇| 欧美激情 高清一区二区三区| 国产女主播在线喷水免费视频网站| 色婷婷av一区二区三区视频| xxxhd国产人妻xxx| 97人妻天天添夜夜摸| 亚洲国产毛片av蜜桃av| 国产男女超爽视频在线观看| av免费观看日本| 久久热在线av| 老司机影院成人| 亚洲少妇的诱惑av| 性高湖久久久久久久久免费观看| 国产一区二区三区av在线| 亚洲欧洲日产国产| 亚洲欧美中文字幕日韩二区| 老熟女久久久| 久久精品人人爽人人爽视色| 欧美bdsm另类| 午夜免费鲁丝| 久久精品久久久久久久性| 色5月婷婷丁香| 亚洲精品久久久久久婷婷小说| 在线观看三级黄色| 18禁在线无遮挡免费观看视频| 国产精品蜜桃在线观看| 亚洲性久久影院| 最新的欧美精品一区二区| 亚洲欧洲日产国产| 99热6这里只有精品| av视频免费观看在线观看| 亚洲激情五月婷婷啪啪| av在线观看视频网站免费| 母亲3免费完整高清在线观看 | 精品亚洲乱码少妇综合久久| 国产有黄有色有爽视频| www.色视频.com| 一二三四中文在线观看免费高清| 亚洲精品中文字幕在线视频| 国产欧美亚洲国产| 国产极品粉嫩免费观看在线| 日日摸夜夜添夜夜爱| 黑人巨大精品欧美一区二区蜜桃 | av不卡在线播放| 欧美97在线视频| 亚洲精品视频女| av一本久久久久| 日韩 亚洲 欧美在线| 夜夜爽夜夜爽视频| 一级片免费观看大全| 亚洲一级一片aⅴ在线观看| 午夜老司机福利剧场| 亚洲精品一区蜜桃| 国产男女超爽视频在线观看| 成人国语在线视频| 国产精品一区www在线观看| 97超碰精品成人国产| 成人毛片a级毛片在线播放| 国产男女内射视频| 久久国内精品自在自线图片| 日本欧美视频一区| 久久99一区二区三区| 少妇的逼好多水| 午夜激情av网站| 精品视频人人做人人爽| 成年美女黄网站色视频大全免费| 99久久中文字幕三级久久日本| 人妻人人澡人人爽人人| 新久久久久国产一级毛片| 亚洲中文av在线| 精品久久久久久电影网| 王馨瑶露胸无遮挡在线观看| 十八禁高潮呻吟视频| 日本免费在线观看一区| 国产高清国产精品国产三级| av在线app专区| 亚洲精品美女久久久久99蜜臀 | 免费高清在线观看视频在线观看| 免费av中文字幕在线| 美女国产高潮福利片在线看| 国产xxxxx性猛交| 中文天堂在线官网| 国产精品三级大全| 久久精品夜色国产| 高清视频免费观看一区二区| 中文字幕人妻熟女乱码| 亚洲高清免费不卡视频| 91成人精品电影| 久久99一区二区三区| 多毛熟女@视频| a级毛色黄片| 人妻一区二区av| 一区在线观看完整版| 欧美日韩国产mv在线观看视频| 久久久久久久亚洲中文字幕| 亚洲精品乱码久久久久久按摩| 伦理电影免费视频| 性色av一级| 99国产精品免费福利视频| 哪个播放器可以免费观看大片| 中文字幕av电影在线播放| av线在线观看网站| 国产成人精品久久久久久| 国产精品欧美亚洲77777| 欧美成人午夜精品| 亚洲在久久综合| 亚洲欧洲国产日韩| 一本久久精品| 久久这里有精品视频免费| 国产黄频视频在线观看| 国精品久久久久久国模美| 乱人伦中国视频| 肉色欧美久久久久久久蜜桃| 日韩三级伦理在线观看| 久久精品夜色国产| 人人妻人人爽人人添夜夜欢视频| 免费观看a级毛片全部| 成人无遮挡网站| 久久狼人影院| 男人操女人黄网站| 日韩欧美精品免费久久| 国产日韩欧美在线精品| 日韩av在线免费看完整版不卡| 免费观看无遮挡的男女| 亚洲国产毛片av蜜桃av| 午夜福利在线观看免费完整高清在| 又粗又硬又长又爽又黄的视频| 在线 av 中文字幕| 少妇的逼水好多| 如日韩欧美国产精品一区二区三区| 毛片一级片免费看久久久久| 夜夜爽夜夜爽视频| 国产黄频视频在线观看| 久久韩国三级中文字幕| 老女人水多毛片| 亚洲国产精品一区二区三区在线| 久久综合国产亚洲精品| 免费黄色在线免费观看| 成年av动漫网址| 国产男女超爽视频在线观看| 日日撸夜夜添| 黄色视频在线播放观看不卡| 亚洲,欧美,日韩| 亚洲国产av影院在线观看| 99久久精品国产国产毛片| 成人综合一区亚洲| 国产一区二区三区av在线| 最后的刺客免费高清国语| 成年动漫av网址| 精品人妻偷拍中文字幕| 男女下面插进去视频免费观看 | 乱人伦中国视频| 肉色欧美久久久久久久蜜桃| 一级a做视频免费观看| 久久久久久人妻| 男女免费视频国产| 亚洲欧美成人综合另类久久久| 久久韩国三级中文字幕| 秋霞伦理黄片| 久久久久精品久久久久真实原创| 老司机影院成人| 国产不卡av网站在线观看| 亚洲av成人精品一二三区| 亚洲四区av| kizo精华| 国产麻豆69| 少妇人妻精品综合一区二区| 在线精品无人区一区二区三| 又粗又硬又长又爽又黄的视频| 如何舔出高潮| 免费观看av网站的网址| 一级毛片电影观看| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 黑人欧美特级aaaaaa片| 精品亚洲乱码少妇综合久久| 国产成人精品无人区| 国产成人av激情在线播放| 国产日韩欧美亚洲二区| 日韩一本色道免费dvd| 卡戴珊不雅视频在线播放| 日韩视频在线欧美| 波多野结衣一区麻豆| 欧美国产精品va在线观看不卡| 伦精品一区二区三区| 国产成人免费观看mmmm| 久久人人爽人人片av| 色94色欧美一区二区| 多毛熟女@视频| 内地一区二区视频在线| 天天操日日干夜夜撸| 中文字幕最新亚洲高清| 亚洲av在线观看美女高潮| 亚洲欧美一区二区三区黑人 | 日本av手机在线免费观看| 久久国产精品男人的天堂亚洲 | 女人被躁到高潮嗷嗷叫费观| 欧美xxⅹ黑人| 久久久久久久精品精品| 日韩视频在线欧美| 国产一区二区在线观看日韩| 欧美日韩av久久| 国产免费现黄频在线看| 亚洲欧洲精品一区二区精品久久久 | 精品99又大又爽又粗少妇毛片| 男女免费视频国产| 卡戴珊不雅视频在线播放| 国产不卡av网站在线观看| 精品午夜福利在线看| av在线老鸭窝| 免费久久久久久久精品成人欧美视频 | 国产不卡av网站在线观看| 一本色道久久久久久精品综合| 大片电影免费在线观看免费| 午夜福利视频精品| 欧美日本中文国产一区发布| 高清av免费在线| 在线观看人妻少妇| freevideosex欧美| 毛片一级片免费看久久久久| 亚洲图色成人| 国产成人精品在线电影| 亚洲av成人精品一二三区| 国产1区2区3区精品| 国产爽快片一区二区三区| 考比视频在线观看| 母亲3免费完整高清在线观看 | 老司机亚洲免费影院| 丝袜喷水一区| 制服丝袜香蕉在线| 国产精品秋霞免费鲁丝片| 欧美性感艳星| 国产精品一二三区在线看| 久久av网站| 亚洲美女黄色视频免费看| 丰满乱子伦码专区| 国产精品一区www在线观看| 成年动漫av网址| 插逼视频在线观看| 日韩大片免费观看网站| 久久午夜福利片| 男女午夜视频在线观看 | 国产男人的电影天堂91| 夜夜骑夜夜射夜夜干| videos熟女内射| 日本-黄色视频高清免费观看| 大片电影免费在线观看免费| 成年人午夜在线观看视频| 亚洲精品久久成人aⅴ小说| 国产福利在线免费观看视频| 久久久久久久久久久免费av| 国产精品偷伦视频观看了| 一本大道久久a久久精品| 精品人妻熟女毛片av久久网站| 99热网站在线观看| 久久99热6这里只有精品| 国产有黄有色有爽视频| 精品久久国产蜜桃| 亚洲精品一区蜜桃| 一区在线观看完整版| 一二三四在线观看免费中文在 | 国产精品国产三级专区第一集| 91精品国产国语对白视频| 在线观看免费视频网站a站| 成人国产av品久久久| 在线观看一区二区三区激情| 观看av在线不卡| 侵犯人妻中文字幕一二三四区| 亚洲国产色片| 免费人妻精品一区二区三区视频| 久久久久久人人人人人| 成人国产av品久久久| 99热6这里只有精品| av播播在线观看一区| 欧美日韩亚洲高清精品| 亚洲av.av天堂| 夜夜爽夜夜爽视频| 成人18禁高潮啪啪吃奶动态图| 欧美人与善性xxx| 久久国产亚洲av麻豆专区| 人体艺术视频欧美日本| 男人爽女人下面视频在线观看| 亚洲av国产av综合av卡| 丝袜在线中文字幕| av免费在线看不卡| 午夜福利乱码中文字幕| 天堂8中文在线网| 看十八女毛片水多多多| 老司机亚洲免费影院| 亚洲精品自拍成人| 丝袜在线中文字幕| 日韩大片免费观看网站| 亚洲人成网站在线观看播放| 精品国产一区二区久久| av线在线观看网站| 久久人人爽人人爽人人片va| 国产毛片在线视频| 99久久综合免费| 免费观看av网站的网址| 少妇熟女欧美另类| 看非洲黑人一级黄片| 国产亚洲精品第一综合不卡 | 欧美激情国产日韩精品一区| 99九九在线精品视频| 亚洲av日韩在线播放| 侵犯人妻中文字幕一二三四区| 哪个播放器可以免费观看大片| 亚洲伊人色综图| 久久人人爽人人爽人人片va| 人人妻人人爽人人添夜夜欢视频| videossex国产| 男女下面插进去视频免费观看 | 国产成人91sexporn| 精品视频人人做人人爽| 青青草视频在线视频观看| 国产一区二区激情短视频 | 亚洲国产精品国产精品| 免费久久久久久久精品成人欧美视频 | 少妇被粗大的猛进出69影院 | 成人国产av品久久久| 日韩av不卡免费在线播放| 天堂8中文在线网| 亚洲少妇的诱惑av| 国产精品 国内视频| 亚洲精品久久午夜乱码| 涩涩av久久男人的天堂| 精品人妻偷拍中文字幕| 看免费av毛片| 日本91视频免费播放| 天天操日日干夜夜撸| xxxhd国产人妻xxx| 成人综合一区亚洲| 日本午夜av视频| 老司机影院成人| 99re6热这里在线精品视频| 777米奇影视久久| 久久国产精品男人的天堂亚洲 | 中文字幕制服av| 五月玫瑰六月丁香| 亚洲,欧美,日韩| 午夜福利影视在线免费观看| 丝袜人妻中文字幕| 97在线人人人人妻| 国产高清国产精品国产三级| 久久女婷五月综合色啪小说| 亚洲欧美一区二区三区国产| 国产国语露脸激情在线看| 男女无遮挡免费网站观看| 午夜免费鲁丝| 制服丝袜香蕉在线| 婷婷成人精品国产| 日韩中文字幕视频在线看片| 欧美+日韩+精品| 人妻人人澡人人爽人人| 免费人妻精品一区二区三区视频| 精品亚洲成a人片在线观看| 精品一区在线观看国产| 国产成人精品在线电影| 少妇人妻精品综合一区二区| 一区二区av电影网| 另类亚洲欧美激情| 丝袜美足系列| 日韩av不卡免费在线播放| 内地一区二区视频在线| 一级毛片电影观看| 国产男女超爽视频在线观看| 一级黄片播放器| 国产成人91sexporn| 国产精品人妻久久久久久| 9191精品国产免费久久| 日产精品乱码卡一卡2卡三| 精品人妻一区二区三区麻豆| 久久热在线av| 99视频精品全部免费 在线| 麻豆精品久久久久久蜜桃| 极品人妻少妇av视频| kizo精华| 婷婷色综合www| 如日韩欧美国产精品一区二区三区| 国产一区有黄有色的免费视频| 黄片播放在线免费| 精品熟女少妇av免费看| 亚洲精品中文字幕在线视频| xxx大片免费视频| 亚洲精品美女久久久久99蜜臀 | 午夜福利网站1000一区二区三区| 曰老女人黄片| 日韩av在线免费看完整版不卡| 日韩一区二区三区影片| 在线看a的网站| 国产精品99久久99久久久不卡 | 久久久久久久国产电影| 少妇被粗大猛烈的视频| 男人舔女人的私密视频| 一区二区av电影网| 国产亚洲精品第一综合不卡 | 在线观看免费视频网站a站| 久久国产亚洲av麻豆专区| 久久免费观看电影| 老熟女久久久| 欧美精品av麻豆av| 免费观看av网站的网址| 晚上一个人看的免费电影| 高清av免费在线| 国产免费现黄频在线看| 日韩一区二区三区影片| 日韩中文字幕视频在线看片| www.av在线官网国产| 97人妻天天添夜夜摸| 精品一区二区免费观看| 宅男免费午夜| 老司机亚洲免费影院| 免费播放大片免费观看视频在线观看| 另类亚洲欧美激情| 夜夜骑夜夜射夜夜干| 婷婷色综合大香蕉| 久热这里只有精品99| 最后的刺客免费高清国语| 欧美精品一区二区免费开放| 男女无遮挡免费网站观看| 亚洲高清免费不卡视频| 成人手机av| 夜夜爽夜夜爽视频| 亚洲成人手机| 欧美亚洲 丝袜 人妻 在线| 免费少妇av软件| 美女福利国产在线| 久久久久精品性色| 亚洲欧美成人综合另类久久久| 免费看不卡的av| 亚洲欧美清纯卡通| 日韩中字成人| 18+在线观看网站| 日本欧美国产在线视频| 欧美bdsm另类| 18禁国产床啪视频网站| 国产精品熟女久久久久浪| 亚洲av成人精品一二三区| 国产伦理片在线播放av一区| 成年人午夜在线观看视频| 少妇人妻久久综合中文| 不卡视频在线观看欧美| 久久久久国产精品人妻一区二区|