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

    免费av不卡在线播放| 婷婷精品国产亚洲av| 1000部很黄的大片| 日韩欧美三级三区| 最近最新中文字幕大全免费视频| 国产欧美日韩一区二区精品| 亚洲av成人精品一区久久| 网址你懂的国产日韩在线| 色在线成人网| 欧美日韩乱码在线| 成年女人毛片免费观看观看9| 亚洲av成人精品一区久久| 欧美日韩综合久久久久久 | 国产精品野战在线观看| 欧美三级亚洲精品| 亚洲精品色激情综合| 午夜福利在线观看吧| 老司机福利观看| 国产精品99久久久久久久久| 美女大奶头视频| 亚洲性夜色夜夜综合| 黄色片一级片一级黄色片| 天堂网av新在线| 中文亚洲av片在线观看爽| 色综合亚洲欧美另类图片| 亚洲狠狠婷婷综合久久图片| 欧美3d第一页| 搞女人的毛片| 18禁黄网站禁片免费观看直播| 久久久久久大精品| 在线看三级毛片| 精品久久久久久久人妻蜜臀av| 色老头精品视频在线观看| 村上凉子中文字幕在线| 日本a在线网址| 欧美黄色淫秽网站| 欧美日韩瑟瑟在线播放| 搡女人真爽免费视频火全软件 | 一级毛片女人18水好多| 91字幕亚洲| 久久香蕉国产精品| 熟女人妻精品中文字幕| 18+在线观看网站| 又粗又爽又猛毛片免费看| 国产一区二区在线av高清观看| 一本久久中文字幕| av专区在线播放| 三级国产精品欧美在线观看| 亚洲熟妇中文字幕五十中出| 51国产日韩欧美| 窝窝影院91人妻| 精品一区二区三区人妻视频| 嫩草影院精品99| 很黄的视频免费| 中文字幕久久专区| 99国产极品粉嫩在线观看| 欧美性猛交黑人性爽| 亚洲成人久久爱视频| 中文字幕久久专区| 国产三级在线视频| 欧美成人性av电影在线观看| 蜜桃久久精品国产亚洲av| 少妇高潮的动态图| 97超视频在线观看视频| 嫁个100分男人电影在线观看| 中文字幕精品亚洲无线码一区| 桃红色精品国产亚洲av| 欧美日韩黄片免| 国产精品综合久久久久久久免费| 国产亚洲精品久久久久久毛片| 又紧又爽又黄一区二区| 中文字幕av成人在线电影| 国内精品美女久久久久久| 动漫黄色视频在线观看| 老汉色∧v一级毛片| 男女午夜视频在线观看| 精品国产美女av久久久久小说| 岛国视频午夜一区免费看| a级毛片a级免费在线| 久久精品91无色码中文字幕| 无人区码免费观看不卡| 亚洲人成网站在线播| 美女被艹到高潮喷水动态| 岛国在线观看网站| 国产精品综合久久久久久久免费| 最好的美女福利视频网| 免费av毛片视频| 亚洲av第一区精品v没综合| 亚洲五月天丁香| 精品久久久久久久久久久久久| 婷婷精品国产亚洲av| 国产91精品成人一区二区三区| 色播亚洲综合网| 亚洲国产精品久久男人天堂| 午夜福利在线观看免费完整高清在 | 成人亚洲精品av一区二区| www.色视频.com| 久久6这里有精品| 日本在线视频免费播放| 白带黄色成豆腐渣| 91麻豆精品激情在线观看国产| 国产三级黄色录像| 男人和女人高潮做爰伦理| 看免费av毛片| 国产爱豆传媒在线观看| 叶爱在线成人免费视频播放| 日本三级黄在线观看| 99视频精品全部免费 在线| 久久99热这里只有精品18| 国产蜜桃级精品一区二区三区| 最新中文字幕久久久久| 亚洲性夜色夜夜综合| 亚洲五月婷婷丁香| 狂野欧美白嫩少妇大欣赏| 中文字幕久久专区| 成人三级黄色视频| www.www免费av| 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 哪里可以看免费的av片| www国产在线视频色| 国语自产精品视频在线第100页| 亚洲成人中文字幕在线播放| 桃红色精品国产亚洲av| 可以在线观看毛片的网站| 日韩免费av在线播放| 国产成人福利小说| 欧美丝袜亚洲另类 | svipshipincom国产片| 小说图片视频综合网站| or卡值多少钱| 尤物成人国产欧美一区二区三区| 国产极品精品免费视频能看的| 日韩欧美一区二区三区在线观看| 嫁个100分男人电影在线观看| 真实男女啪啪啪动态图| 三级国产精品欧美在线观看| а√天堂www在线а√下载| 国产高清videossex| 亚洲一区二区三区色噜噜| 变态另类丝袜制服| 成年免费大片在线观看| 在线天堂最新版资源| 19禁男女啪啪无遮挡网站| 国产97色在线日韩免费| 国产亚洲精品久久久久久毛片| 中文字幕人妻丝袜一区二区| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费 | 中文字幕av在线有码专区| 99国产极品粉嫩在线观看| 成人午夜高清在线视频| 天堂av国产一区二区熟女人妻| 一区二区三区高清视频在线| 母亲3免费完整高清在线观看| 色尼玛亚洲综合影院| 999久久久精品免费观看国产| 亚洲内射少妇av| 国产69精品久久久久777片| 欧美日韩精品网址| 成年女人永久免费观看视频| 特大巨黑吊av在线直播| 欧美一区二区精品小视频在线| 国产精品日韩av在线免费观看| x7x7x7水蜜桃| 日本成人三级电影网站| 欧美性猛交╳xxx乱大交人| 亚洲avbb在线观看| 精品国产亚洲在线| 俄罗斯特黄特色一大片| 男女那种视频在线观看| 丁香欧美五月| 日本一二三区视频观看| 国产亚洲精品一区二区www| 少妇丰满av| 男插女下体视频免费在线播放| 免费在线观看日本一区| 精品国产三级普通话版| 国产视频内射| 久久香蕉精品热| 色综合站精品国产| 午夜久久久久精精品| 日韩免费av在线播放| 成人午夜高清在线视频| 人人妻人人看人人澡| 日本 av在线| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 99热只有精品国产| 国产v大片淫在线免费观看| 国产亚洲精品久久久久久毛片| 国产黄a三级三级三级人| 久久精品国产99精品国产亚洲性色| 一个人看视频在线观看www免费 | 久久精品亚洲精品国产色婷小说| 亚洲激情在线av| 国产麻豆成人av免费视频| 国产精品乱码一区二三区的特点| 在线免费观看的www视频| 国产麻豆成人av免费视频| 国产乱人伦免费视频| 老汉色av国产亚洲站长工具| 最新在线观看一区二区三区| 日韩欧美 国产精品| 狠狠狠狠99中文字幕| 中文在线观看免费www的网站| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 叶爱在线成人免费视频播放| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| 亚洲成a人片在线一区二区| 天天一区二区日本电影三级| 亚洲精品成人久久久久久| 男女视频在线观看网站免费| 中文在线观看免费www的网站| 精品不卡国产一区二区三区| 亚洲熟妇熟女久久| 无人区码免费观看不卡| 一级毛片高清免费大全| 国产男靠女视频免费网站| 69av精品久久久久久| 精品久久久久久久久久久久久| 欧美乱色亚洲激情| 国产三级在线视频| 啪啪无遮挡十八禁网站| 国内精品久久久久久久电影| 久久久久免费精品人妻一区二区| 麻豆成人av在线观看| 国产精品免费一区二区三区在线| 一夜夜www| 桃色一区二区三区在线观看| 午夜视频国产福利| 久久人人精品亚洲av| 中文字幕熟女人妻在线| 在线看三级毛片| 亚洲一区二区三区不卡视频| 欧美色视频一区免费| 丁香六月欧美| 波多野结衣高清无吗| 69人妻影院| 亚洲男人的天堂狠狠| 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 久久九九热精品免费| 可以在线观看毛片的网站| av天堂在线播放| 香蕉丝袜av| 成年女人毛片免费观看观看9| 亚洲午夜理论影院| 精品人妻1区二区| 超碰av人人做人人爽久久 | 成熟少妇高潮喷水视频| 欧美最黄视频在线播放免费| 日本黄色视频三级网站网址| 天美传媒精品一区二区| 国产av麻豆久久久久久久| 一级a爱片免费观看的视频| 性色avwww在线观看| 一区二区三区国产精品乱码| 他把我摸到了高潮在线观看| 亚洲精品456在线播放app | 少妇人妻精品综合一区二区 | 久久午夜亚洲精品久久| 老司机深夜福利视频在线观看| 国产主播在线观看一区二区| xxxwww97欧美| 午夜福利高清视频| 久久久久国内视频| 欧美日本视频| 波多野结衣高清作品| 性欧美人与动物交配| 国产亚洲精品久久久久久毛片| 有码 亚洲区| 亚洲avbb在线观看| 亚洲欧美日韩无卡精品| 在线观看免费午夜福利视频| 国产亚洲精品久久久久久毛片| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 国产伦人伦偷精品视频| 看片在线看免费视频| 午夜福利欧美成人| 国产亚洲欧美98| 成人特级黄色片久久久久久久| 国产一区二区亚洲精品在线观看| 99热这里只有精品一区| 夜夜看夜夜爽夜夜摸| 搡老妇女老女人老熟妇| 在线国产一区二区在线| 日本精品一区二区三区蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 国内精品美女久久久久久| 一夜夜www| 免费人成在线观看视频色| 亚洲精品成人久久久久久| 亚洲人成电影免费在线| 一区二区三区激情视频| 一个人观看的视频www高清免费观看| 亚洲精品乱码久久久v下载方式 | 亚洲av一区综合| 久久伊人香网站| 亚洲 欧美 日韩 在线 免费| 老汉色av国产亚洲站长工具| 午夜激情欧美在线| 波多野结衣巨乳人妻| 日韩欧美精品v在线| or卡值多少钱| 久久午夜亚洲精品久久| 啦啦啦韩国在线观看视频| 亚洲人成网站在线播放欧美日韩| 精品久久久久久,| 乱人视频在线观看| 亚洲一区高清亚洲精品| 在线十欧美十亚洲十日本专区| 2021天堂中文幕一二区在线观| 国产欧美日韩精品亚洲av| 高清毛片免费观看视频网站| 日本三级黄在线观看| 99久久九九国产精品国产免费| 国产精品久久久久久亚洲av鲁大| 90打野战视频偷拍视频| 国产真实乱freesex| 淫妇啪啪啪对白视频| 国产av在哪里看| 午夜福利视频1000在线观看| 亚洲av中文字字幕乱码综合| 久久精品夜夜夜夜夜久久蜜豆| 国产精品久久视频播放| 国产亚洲精品综合一区在线观看| 日本 av在线| 熟女人妻精品中文字幕| 91字幕亚洲| 麻豆成人午夜福利视频| 久久久久久九九精品二区国产| a在线观看视频网站| 免费av不卡在线播放| 亚洲黑人精品在线| 免费av不卡在线播放| 亚洲黑人精品在线| 波野结衣二区三区在线 | 99热这里只有是精品50| 十八禁网站免费在线| 国产精品一区二区免费欧美| 亚洲18禁久久av| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 午夜久久久久精精品| 一进一出好大好爽视频| av黄色大香蕉| 久久这里只有精品中国| 在线免费观看不下载黄p国产 | 又紧又爽又黄一区二区| 午夜免费激情av| 亚洲中文字幕日韩| 好男人在线观看高清免费视频| 国产亚洲精品久久久久久毛片| 久久久久久久精品吃奶| 精品久久久久久久人妻蜜臀av| 中文字幕人成人乱码亚洲影| 观看美女的网站| 国产精品久久久人人做人人爽| 国产精品一区二区三区四区久久| 高清毛片免费观看视频网站| netflix在线观看网站| 久久中文看片网| 禁无遮挡网站| 一本综合久久免费| 日韩欧美三级三区| 精品无人区乱码1区二区| 午夜福利在线观看免费完整高清在 | 手机成人av网站| 成人特级黄色片久久久久久久| 欧美一区二区精品小视频在线| 国产精品一及| 亚洲乱码一区二区免费版| 免费高清视频大片| 国产v大片淫在线免费观看| 精品国产超薄肉色丝袜足j| 免费观看人在逋| 久久国产精品人妻蜜桃| 日韩欧美国产一区二区入口| 成人无遮挡网站| 9191精品国产免费久久| 他把我摸到了高潮在线观看| 在线看三级毛片| 午夜福利在线观看免费完整高清在 | 美女高潮的动态| 色哟哟哟哟哟哟| 亚洲五月天丁香| 香蕉久久夜色| 国产成年人精品一区二区| 久久久国产成人精品二区| 身体一侧抽搐| 69av精品久久久久久| 91在线观看av| 有码 亚洲区| 中文在线观看免费www的网站| 神马国产精品三级电影在线观看| 国产精华一区二区三区| 日韩免费av在线播放| 99热6这里只有精品| 日本五十路高清| 狂野欧美白嫩少妇大欣赏| 999久久久精品免费观看国产| 美女高潮喷水抽搐中文字幕| a级一级毛片免费在线观看| 熟妇人妻久久中文字幕3abv| 日韩欧美一区二区三区在线观看| 夜夜爽天天搞| 怎么达到女性高潮| 国产蜜桃级精品一区二区三区| 日韩亚洲欧美综合| 亚洲欧美日韩东京热| 国产三级黄色录像| 国产在视频线在精品| 国产午夜精品久久久久久一区二区三区 | 婷婷亚洲欧美| 亚洲国产欧洲综合997久久,| 蜜桃久久精品国产亚洲av| 午夜免费成人在线视频| 高清毛片免费观看视频网站| 免费无遮挡裸体视频| 免费观看人在逋| 免费av观看视频| 国产成人福利小说| 男女那种视频在线观看| 3wmmmm亚洲av在线观看| 在线播放国产精品三级| 五月玫瑰六月丁香| netflix在线观看网站| 18禁美女被吸乳视频| 国产精品99久久99久久久不卡| 亚洲aⅴ乱码一区二区在线播放| 久久精品亚洲精品国产色婷小说| 两个人视频免费观看高清| av天堂中文字幕网| 成人av在线播放网站| 男女午夜视频在线观看| 欧美中文综合在线视频| 十八禁网站免费在线| 色综合婷婷激情| 国内毛片毛片毛片毛片毛片| www日本在线高清视频| 欧美黑人巨大hd| 亚洲 欧美 日韩 在线 免费| e午夜精品久久久久久久| 午夜福利成人在线免费观看| 好男人电影高清在线观看| 亚洲成av人片免费观看| 女同久久另类99精品国产91| 制服丝袜大香蕉在线| 中文字幕久久专区| 色噜噜av男人的天堂激情| 国产成人aa在线观看| 亚洲 欧美 日韩 在线 免费| 久久天躁狠狠躁夜夜2o2o| 免费看光身美女| 国产精品自产拍在线观看55亚洲| 精品99又大又爽又粗少妇毛片 | 亚洲精品影视一区二区三区av| 在线播放无遮挡| 亚洲精品国产精品久久久不卡| 色av中文字幕| 91av网一区二区| 桃色一区二区三区在线观看| 日本免费一区二区三区高清不卡| 国产精品久久视频播放| 精品福利观看| 久久6这里有精品| 在线十欧美十亚洲十日本专区| or卡值多少钱| 欧美日本视频| 90打野战视频偷拍视频| 免费一级毛片在线播放高清视频| 88av欧美| 国产成人a区在线观看| 高清日韩中文字幕在线| 内地一区二区视频在线| 日本 av在线| 久久伊人香网站| 国产精品av视频在线免费观看| 久久草成人影院| 成人国产一区最新在线观看| 在线观看免费午夜福利视频| 一区二区三区高清视频在线| www国产在线视频色| 宅男免费午夜| 欧美日本视频| 久久久久国内视频| 久久久久亚洲av毛片大全| 老司机午夜福利在线观看视频| 国产午夜精品久久久久久一区二区三区 | 欧美成人a在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 久久这里只有精品中国| 首页视频小说图片口味搜索| 日韩欧美 国产精品| 亚洲最大成人手机在线| 国产探花极品一区二区| 男女做爰动态图高潮gif福利片| 99riav亚洲国产免费| 可以在线观看毛片的网站| 欧美黑人欧美精品刺激| 久久久久性生活片| 成人特级黄色片久久久久久久| 噜噜噜噜噜久久久久久91| 真人做人爱边吃奶动态| 国产精品久久久人人做人人爽| 久9热在线精品视频| 欧美区成人在线视频| 欧美黑人巨大hd| 国产精品1区2区在线观看.| 俄罗斯特黄特色一大片| 好看av亚洲va欧美ⅴa在| 免费大片18禁| 特级一级黄色大片| 搡老妇女老女人老熟妇| 成人欧美大片| 亚洲av成人av| 成人欧美大片| 成年女人永久免费观看视频| 在线国产一区二区在线| 淫秽高清视频在线观看| 一本精品99久久精品77| 99久久精品热视频| 丰满人妻熟妇乱又伦精品不卡| 人人妻,人人澡人人爽秒播| 在线看三级毛片| 日韩欧美国产一区二区入口| 亚洲成a人片在线一区二区| 午夜影院日韩av| 亚洲av不卡在线观看| 波多野结衣高清作品| 亚洲欧美日韩东京热| 99在线视频只有这里精品首页| 最好的美女福利视频网| 黑人欧美特级aaaaaa片| 国产黄色小视频在线观看| 国产一级毛片七仙女欲春2| 美女高潮的动态| 亚洲美女视频黄频| 亚洲男人的天堂狠狠| 91麻豆精品激情在线观看国产| 成人亚洲精品av一区二区| 亚洲成人中文字幕在线播放| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 国内精品美女久久久久久| 亚洲黑人精品在线| 国产亚洲精品一区二区www| 免费高清视频大片| 淫妇啪啪啪对白视频| 国产精品日韩av在线免费观看| 亚洲欧美精品综合久久99| 国产野战对白在线观看| 在线观看日韩欧美| 老汉色∧v一级毛片| 成人av在线播放网站| 又粗又爽又猛毛片免费看| 欧美日韩一级在线毛片| 麻豆久久精品国产亚洲av| 欧美乱妇无乱码| 麻豆久久精品国产亚洲av| 日韩有码中文字幕| 午夜福利在线在线| 亚洲精品456在线播放app | 亚洲精品乱码久久久v下载方式 | 99热这里只有精品一区| 99久久精品热视频| 精品久久久久久久久久免费视频| 深夜精品福利| 成人精品一区二区免费| 午夜亚洲福利在线播放| 一个人看的www免费观看视频| 老鸭窝网址在线观看| 中文亚洲av片在线观看爽| 欧美日韩一级在线毛片| 日本与韩国留学比较| 蜜桃久久精品国产亚洲av| 88av欧美| 少妇的逼水好多| 狠狠狠狠99中文字幕| 女同久久另类99精品国产91| 美女 人体艺术 gogo| 中文字幕人成人乱码亚洲影| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 叶爱在线成人免费视频播放| 999久久久精品免费观看国产| 狂野欧美激情性xxxx| 麻豆国产97在线/欧美| 黄色成人免费大全| 国产视频内射| 亚洲av免费在线观看| 白带黄色成豆腐渣| e午夜精品久久久久久久| 国产欧美日韩精品亚洲av| 老司机在亚洲福利影院| e午夜精品久久久久久久| 天堂√8在线中文| 久久国产精品影院| 久久草成人影院| 久久精品影院6| 国语自产精品视频在线第100页| 热99re8久久精品国产| 少妇人妻一区二区三区视频| 国产激情欧美一区二区| 日本 欧美在线| 亚洲欧美日韩无卡精品| 精品久久久久久久人妻蜜臀av| 欧美黑人巨大hd| 亚洲成人精品中文字幕电影|