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

    Preliminary Simulation and Manufacture of a 40-element PAF Prototype for FAST

    2023-09-03 01:36:30JunWangChengjinJinBoPengandYanZhu

    Jun Wang ,Chengjin Jin ,Bo Peng ,and Yan Zhu

    1 National Astronomical Observatories,Chinese Academy of Sciences,Beijing 100101,China;wangjun@nao.cas.cn

    2 University of Chinese Academy of Sciences,Beijing 100049,China

    Abstract This paper presents the preliminary simulation results of the prototype 40-element phased array feed (PAF) array for the Five-hundred-meter Aperture Spherical radio Telescope (FAST) and the measurement of a manufactured single Vivaldi element.The simulation results include the characteristic of a standalone Vivaldi element.The fullwave simulation of the explicit PAF array is carried out,from which the performance of the embedded elements and the mutual coupling among the elements could be accurately inferred.The performance of the combination of the PAF and the main reflector of FAST is calculated.The results imply that the PAF combined with the FAST antenna could operate over a much wider bandwidth from 650 to 1450 MHz compared with the standalone element.The antenna temperature remains almost the same for zenith angles from 0 to 40°.From the simulation results,the PAF shows good potential to improve the performance of FAST.A measurement of a manufactured Vivaldi element is also presented.

    Key words: instrumentation: detectors– methods: observational– telescopes

    1.Introduction

    The Five-hundred-meter Aperture Spherical radio Telescope(FAST) is the largest and the most sensitive single dish radio telescope in the world (Nan 2006).Following the innovative adoption of an active main reflector and a lightweight focus suspension mechanism,a traditional point feed can be utilized on FAST.The neutral shape of the main reflector of FAST is spherical,with a radius of 300 m and an aperture of 500 m in diameter.When doing observations,part of the main reflector is deformed into a paraboloid.The diameter of the aperture of the actively formed paraboloid is about 300 m.The metal surface outside of the paraboloid is still in its neutral spherical position.This outer spherical metal surface blocks thermal emission from the ground effectively.As the zenith angle becomes larger,the edge of the paraboloid moves to the edge of the main reflector.When the zenith angle is larger than 26°.4,the paraboloid starts to be truncated by the edge of the main reflector.This results in a gain reduction and an increased system noise temperature.FAST can be seen as a prime-focus radio telescope with a zenith angle less than 40°.The gain starts to reduce when the zenith angle is larger than 26°.4.The precision of the fitting to the paraboloid by the active main reflector,combined with the accuracy of positioning the feed horns,dictates the high-frequency limit of FAST.

    The current receiver installed on FAST is a 19-horn receiver(Dunning et al.2019).There are 19 identical horns closely packed in a hexagonal pattern,with each horn working independently.Compared with the horn at the center,the offcentered horns present reduced gain due to the coma effect.The reduction in gain becomes larger as the off-center distance increases.There are gaps in the sky coverage corresponding to the horns and tessellation mode is usually adopted to fill the gaps by multiple pointing on the sky.

    A phased array feed (PAF) utilizes closely packed small feeding elements on the focal plane and beam-forming is done by a weighted summation of the voltage data received by the feeding elements.PAFs can provide continuous sky coverage by forming 3 dB overlapped beams within the field of view(FoV)and the sensitivity of the formed beams are close or even better than that obtained by using a traditional horn receiver.PAF technology for radio astronomy has been intensively investigated and developed during the last two decades with remarkable achievements.The ambient-temperature PAFs include APERTIF(APERture Tile In Focus)at the Westerbork Synthesis Radio Telescope(van Cappellen et al.2022)and the chequerboard PAF at the Australian Square Kilometre Array Pathfinder radio telescope (Chippendale &Schinckel 2011).Research on cryogenically cooled PAFs include the Focal L-band Array for the Green Bank Telescope (FLAG) (Pingel et al.2021),a prototype for Arecibo (Cortes-Medellin et al.2015),the Advanced L-Band Phased Array Camera for Arecibo (Cortes-Medellin et al.2016) and the Rocket PAF for the Parkes telescope (Dunning et al.2016).

    In this paper,we present the design of a prototype 40-element PAF for FAST,a simulation of the performance of this PAF combined with FAST’s main reflector and the measurement of a manufactured single feeding element.A brief introduction to the PAF and the simulation of the feeding array are presented in the next section.Simulation of the performance of the PAF on FAST is presented in the third section.The measurement of a single manufactured element is presented in the fourth section.Concluding remarks are given in the last section.

    2.A Prototype CryoPAF Designed for FAST

    A 40-element dual-polarization cryogenically cooled PAF for FAST is currently under development.This prototype PAF array consists of 101 identical elements,49 for one linear polarization and 52 for another;the two polarizations are orthogonal to each other.The elements are arranged in a squared grid,as shown in Figure 1.Each element is an allmetal Vivaldi feed,as shown in Figure 3.The interconnection of the elements is done by a slide-in-groove interlocking mechanism as shown in Figure 2.Vivaldi antennae are widely used for wideband operation.The all-metal design helps to achieve low ohmic loss and the inter-locked metal ridges together with the ground plane provide mechanical strength to support the array.If the f/D ratio of a paraboloid is given and the PAF is put at the prime focus,the size of the element is normally chosen to ensure an appreciable gain after proper beam forming.In the case of FAST,the f/D of the actively formed paraboloid of the main reflector is 0.4611.The size(same as the distance between the adjacent elements)is chosen as 125 mm,which is about 0.6 times the wavelength of the radio wave at 1.45 GHz.

    Figure 1.A square array of Vivaldi elements.There are 49 horizontal and 52 vertical elements,which correspond to two orthogonal linear polarizations.The distance between adjacent elements is 125 mm.

    Figure 2.An all-metal Vivaldi feeding element.The slide-in-groove interlocking mechanism is clearly seen at the joining position of the ridges.

    Figure 3.The 3D model of an all-metal Vivaldi feeding element.The labels illustrate the various parts of the element.The element is excited by the coaxial port at its bottom.

    The 3D model of the Vivaldi element is shown in Figure 3.The only excitation is at the coaxial port at the bottom.The thickness and the profile of the ridge,width of the slot at various parts,the size of the resonant cavity and the position of the coaxial probe are adjusted to optimize the return loss of the coaxial port for the frequency range of 1.05–1.45 GHz.The ground plane is a square plate with a size of 125 mm.

    The simulated return loss of a standalone single Vivaldi element is shown in Figure 4.The S11 is below ?20 dB from 1.05 to 1.55 GHz,and below ?15 dB from 1 to 1.65 GHz.The far-field pattern of this isolated vivaldi element is shown in Figure 5.The edge taper of this pattern w.r.t.the paraboloid part of FAST is only roughly ?5 dB,which is certainly not good for a single feed receiver.The high edge taper of this small Vivaldi element is good for PAF,as can be inferred from the simulation results presented in Section 3.

    Figure 4.Simulated S11 of a standalone Vivaldi feeding element.The S11 is below ?20 dB from 1.05 to 1.55 GHz,and below ?15 dB from 1.0 to 1.65 GHz.The orange dashed line indicates a return loss of ?20 dB.

    Figure 5.Far-field pattern of an isolated single Vivaldi element.The taper at 56°.4 is roughly from 4 to 5 dB.

    Imagine a plane wave coming from direction of the optical axis of the paraboloid,it will be reflected by the paraboloid and form an Airy dot on the focal plane.In order to collect the reflected energy on the focal plane,the Vivaldi elements are interconnected to form an array to cover the Airy dot as much and as efficiently as possible.Since the elements are closely packed together,mutual coupling among them is inevitable.The return loss,far-field pattern of any Vivaldi elements in the array will then be different from that of a standalone element.

    A full-wave electromagnetic (EM)simulation of the explicit array as shown in Figure 1 is done using the HFSS software package.Figure 6 shows the coupling level between the central and some of the other elements.At all frequencies,the level of mutual coupling becomes weaker as the distance between the elements increases.The S12,which is the level of the coupling between adjacent elements,decreases as the frequency increases.However,for element pairs with larger distances,the level of mutual coupling becomes larger as the frequency increases.For frequencies higher than 1.5 GHz,the level of mutual coupling does not change much as the distance increases.This may imply that the mutual coupling at higher frequencies is in different modes due to the size of the element relative to wavelength,and this may partly explain the gain reduction of the PAF combined with the main reflector of FAST at higher frequencies,as mentioned in Section 3.1.A PAF with smaller-sized elements would be required if operation at higher frequencies is desired.

    Figure 6.Simulated mutual coupling level between the central element and several other elements,as indicated in the left figure.The coupling normally becomes weaker as the distance between the elements gets larger as expected.

    Figure 7 shows the simulated S11 of all the embedded elements in the array,together with the S11 of a standalone element.The return loss is a bit worse than a standalone element from 1.05 to 1.45 GHz,but is still below ?10 dB.One remarkable phenomenon is that the S11 for frequencies below 1 GHz becomes better because of mutual coupling.This may enable the PAF receiver to operate over a wider frequency range than a standalone element.The simulated gain of the PAF combined with the main reflector of FAST presented in Section 3.1 also illustrates the wideband feature of the PAF.

    Figure 7.Simulated S11 of all the embedded elements in the PAF array and a standalone element(red dashed line).The S11 of embedded elements are below ?10 dB for frequency range of 1.0–1.6 GHz,but lower down to ?15 dB at 650 MHz.This is due to mutual coupling effect among the elements and enables the PAF to operate over a wider frequency range than a standalone element.

    The far-field patterns of the embedded elements are also different compared with a standalone element due to mutual coupling among them.Depending on the position within the array,the far-field patterns of the elements are not the same.Figure 8 shows the far-field patterns of two elements in the PAF array.The left figure is for the element on the edge.Compared with the far-field pattern of the element at the center of the PAF,the far-field pattern of the element at the edge has a somewhat wider beamwidth and a lower peak gain.The element at the edge will only be affected by the elements on one side of it and thus will be affected less than the central element which is surrounded by elements around it.

    Figure 8.Simulated far-field pattern of two elements in the PAF array.The left figure is for one element at the edge of the array,while the right is for the central element.The far-field pattern of the central element shows a narrower beam and a higher on-axis gain.This is due to a stronger coupling effect for the central element.

    This mutual coupling effect will affect the noise matching between the feeding element and the low noise amplifier(LNA),since the noise of the LNA that goes into the feed will partly be coupled into other feeding elements and will go into the beam-forming calculation eventually.Compared with a standalone feed receiver,this additional noise would require co-design of the LNA with the PAF array with a beam-forming algorithm.This effect is usually referred to as active impedance matching.

    3.Simulation of the PAF Combined with the Main Reflector of FAST

    A simulation of the efficiency (gain) and system noise temperature of the PAF combined with the main reflector of FAST is presented in this section.The main reflector is modeled by discrete points separated by approximately 5 m.The PAF array is put at the focal plane,approximately 162.12 m below the center of the neutral spherical surface.

    Considering an incoming plane wave,the response of each element is calculated following the reciprocity principle(Warnick et al.2018).First,the far-field pattern of one element of the PAF combined with the main reflector of FAST is calculated using the GRASP software package.The embedded far-field pattern of the element is taken as the input.The complex voltage of the far-field pattern at a certain direction is then taken as the response of this element,assuming that an incoming plane wave is coming from that particular direction.The effective area of this element combined with the main reflector can be estimated by using Equation (1).

    In Equation(1),the Gain corresponds to the far-field pattern at a certain direction of the embedded element combined with the main reflector of FAST.

    By repeating the above procedure for all the elements,the complex voltage received by each embedded element can be obtained.Notice that the mutual coupling effect among the feeding element is already taken into account because the embedded far-field pattern of each element is used.

    Each feeding element is followed by an LNA,a banddefining filter and an RF circuit.The amplified complex voltages from all the elements are then weighted and summed to form a beam.The weighting factors will affect the gain and system noise temperature of the formed beam.

    3.1.Gain or Effective Collecting Area

    Assuming that the noise of the LNAs following each element are uncorrelated and the mutual coupling among the feeding elements are ignored,the conjugate field matching (CFM)beam-forming mode would provide the highest signal-to-noise ratio (SNR).Though the mutual coupling among the feeding element mentioned in Section 2 would induce the active impedance matching between the feeding elements and the LNAs,this will normally make the receiver noise a bit higher.Thus,the CFM beam-forming mode may not be exact for the highest SNR.The CFM beam-forming will still provide an idea of how much gain could be achieved.In the practical design of a PAF,the mutual coupling among the feeding elements and the beam-forming algorithm may need to be taken into account when designing the matching between the feeding element and the LNA.In the following simulation,the CFM beam-forming mode will be adopted.

    The efficiency relative to a 300 m aperture of the CFM formed beam when pointing to the zenith of one linear polarization at several frequencies is shown in Figure 9.The efficiency of the other orthogonal linear polarization presents similar efficiency.The efficiency is about 0.63 at 1400 MHz and decreases to 0.53 at 950 MHz.However,the gain remains greater than 0.52 down to 650 MHz.At all frequencies,the gain remains almost constant up to an off-axis angle of 0°.06.The efficiency starts to decrease to 0.49 at 1500 MHz and 0.38 at 1650 MHz.

    Figure 9.Simulated efficiency relative to a 300 m aperture at several RF frequencies of one polarization when the telescope is pointing to the zenith.The abscissa is the off-axis angle,which corresponds to the off-axis direction of incoming plane wave.The efficiency for all frequencies remains unchanged up to an off-axis angle of 0.06 and starts to decrease.This is mainly due to the limited size of this PAF prototype.

    3.2.System Noise Temperature

    In the CFM beam-forming mode,the voltages that come from each feeding element are weighted by the conjugate of the complex voltages received by each feeding element if a point source is put at a certain direction within the FoV.Under this weighting scheme,the PAF array itself could be thought of as a feed,and the far-field pattern of this feed could be obtained by the weighted sum of the far-field pattern of each element,taking into account of the different position of each feeding element on the focal plane.

    Figure 10 shows two slices of the far-field patterns of the PAF under CFM beam forming at 1400 MHz when pointing at zenith angles of 0°and 40°.The slices are cut through the plane formed by the optical axis of the paraboloid at zenith angles of 0° and 40°.Clearly,the beams have a sharp cutoff at the edge of the paraboloidal part of the main reflector.This is mainly because the PAF has covered most of the Airy dot on the focal plane,while the traditional feed normally only receives the main lobe of the Airy dot.The antenna temperature is then obtained by integrating this pattern over the whole space,weighted by the radiation temperature of the environment.The sky and the reflector is set to be 5 K at L band,and the ground and mountain is set to be 300 K.The T-antenna is roughly 5–7 K when the antenna is pointing at the zenith.Assuming the receiver temperature to be 10 K,the Tsyswould then be 15–17 K.

    Figure 10.Far-field pattern of the PAF under CFM beam-forming mode.

    In Figure 10,the left part of the blue curve(40°zenith angle)shows similar sharp cutoff compared with the orange curve(0°zenith angle),but shrinks inwardly by a certain amount.This is because the paraboloid is truncated by the ring girder of the main reflector.Compared with a traditional horn feed,this reduces the spillover dramatically.The calculated T-antenna at 40° zenith angle is almost the same as that when the zenith angle is 0°.Thus,the system noise temperature changes very little for zenith angles from 0°to 40°.In the case of a traditional horn feed,the Tsysmay go up to more than 40 K at a zenith angle of 40° and the gain may decrease by nearly 30%.If the PAF prototype presented in this paper is used,the gain may decrease by 20% at a zenith angle of 40°,but the Tsysalmost remains the same as that at a zenith angle of 0°.Thus,the adoption of the PAF may help to obtain higher sensitivity at large zenith angles for FAST.

    3.3.Far-field Pattern of the Formed Beam

    Under a certain weighting scheme,the PAF combined with the reflector will have different response to the incoming plane waves from various directions,which naturally forms a far-field pattern of the combination of the PAF and the main reflector.Under a CFM weighting scheme,the weighting factor to each element happens to be the conjugate of the received complex voltage for an incoming plane wave from a certain direction.The response of all the feeding elements of the PAF is calculated following the method mentioned in the beginning of this section.Applying the weighting factors to the response of all the feeding elements for incoming plane waves from different directions,a far-field pattern of the combination of the PAF and the main reflector of the antenna would then be obtained.

    Figure 11 shows the far-field pattern of the PAF combined with the main reflector of FAST.The overlap of the adjacent beams is approximately at 3 dB,thus enabling continuous sky coverage,which could not be achieved by using multihorn feed array.The far-field patterns of small off-axis angles remains to be symmetric for the near-side lobes.This is mainly because the PAF would cover most of the Airy dot for small off-axis angles,thus the PAF will not suffer from the coma effect that the traditional horn feed often encounters.

    Figure 11.Far-field pattern of the PAF combined with the main reflector of FAST.The adjacent beams overlap at approximately 3 dB to allow for a continuous sky coverage.

    By a finer adjustment of the phase and amplitude of the weighting factor,more closely packed beams could be formed simultaneously,as shown in Figure 12.When observing a compact source with a radio telescope with larger aperture and when the antenna pointing is not perfect,these simultaneous closely packed multiple beams could be adopted to emulate a mechanical scan to calibrate the pointing of the antenna.This electrical pointing mode could make the pointing calibration much faster and could readily be realized by using a PAF receiver.

    Figure 12.The simultaneously formed closely packed multiple beams.The peak gain of these beams are nearly the same.These beams could be used to realize electrical pointing for a large aperture antenna.

    4.Measurement of a Single Element and the Manufacture of the Array

    One single element is manufactured according to the simulation model.The feed is made of bare aluminium and no special coating is used.Figure 13 shows a photo of the manufactured single element.The inner conductor of the N-type connector is silver plated with a thickness of three times the skin depth and coated with a thin layer of gold with 0.05 um thickness to avoid oxidation.

    Figure 13.The manufactured single element with a ground plane and an N-type connector.The element is fed by the coaxial port from the bottom.

    The PAF array is made by interconnecting the individual Vivaldi elements to form a square array.The interconnection is realized by a sliding-in-groove mechanism,as shown in Figure 14.The S-parameter is affected by the addition of the ridges that belongs to other Vivaldi elements in the array,as can be seen from Figure 15.The mutual coupling among the closely packed elements improves the performance of the Vivaldi elements at lower frequencies,thus allowing the PAF array to operate over a wider frequency range compared with a standalone single Vivaldi element.

    Figure 14.The ridges of the Vivaldi elements (left).The single feeding element with the ridges belongs to adjacent installed elements (right).

    Figure 15.The measured S11 of the single element(brown curve)and with the added ridges that belongs to adjacent elements(blue curve).The red dashed line is for the simulated standalone single element.

    5.Conclusion

    A 40-element dual-polarization prototype PAF feeding array is introduced in this paper.EM simulation results manifest the mutual coupling among the closely packed feeding elements.The simulated performance of the combination of the PAF receiver and FAST’s main reflector indicates an efficiency of 62% at a zenith angle of 0° and a system noise temperature of 15–17 K.At zenith angles higher than 26°.4,although the actively formed paraboloid is truncated by the edge of the main reflector,the spillover is still small due to the fact that the PAF covers both the main and side lobes of the airy dot on the focal plane.Compared with the current multi-horn 19 beam receiver,the far-field patterns of PAF combined with FAST’s main reflector show better symmetry for off-axis beams and the PAF may enlarge the frequency coverage.The PAF may enable continuous sky coverage within the FoV.A manufactured single Vivaldi element confirms the simulation result and demonstrates the mutual coupling effect by adding several ridges that belong to adjacent elements.

    Acknowledgments

    This work is supported by the National Key R&D Program of China under grant No.2018YFA0404703,and the Open Project Program of the CAS Key Laboratory of FAST,NAOC,Chinese Academy of Sciences.We acknowledge the CASMPG LEGACY funding “Low Frequency Gravitational Wave Astronomy and Gravitational Physics in Space”.The authors thank Prof.Sander Weinreb for valuable discussions.

    小说图片视频综合网站| 999久久久精品免费观看国产| 国产精品久久久久久亚洲av鲁大| 老汉色av国产亚洲站长工具| 国产亚洲精品久久久com| 午夜成年电影在线免费观看| 欧美绝顶高潮抽搐喷水| 欧美av亚洲av综合av国产av| 麻豆国产97在线/欧美| 国模一区二区三区四区视频 | 国产精品久久久久久久电影 | 色播亚洲综合网| 两人在一起打扑克的视频| 国产av麻豆久久久久久久| 国产精品一区二区精品视频观看| 欧美av亚洲av综合av国产av| av视频在线观看入口| 亚洲人成伊人成综合网2020| 国产伦一二天堂av在线观看| 午夜福利18| 露出奶头的视频| 麻豆成人午夜福利视频| 亚洲av中文字字幕乱码综合| 日韩欧美国产一区二区入口| 久久久国产精品麻豆| 别揉我奶头~嗯~啊~动态视频| 99热这里只有是精品50| 午夜福利高清视频| 国产一区二区在线av高清观看| 午夜精品在线福利| 天堂av国产一区二区熟女人妻| 国产极品精品免费视频能看的| 国产欧美日韩一区二区三| 久久亚洲真实| 观看免费一级毛片| 一区二区三区高清视频在线| 国产午夜精品久久久久久| 精品欧美国产一区二区三| 伦理电影免费视频| 亚洲熟女毛片儿| 香蕉av资源在线| 麻豆av在线久日| svipshipincom国产片| 免费在线观看视频国产中文字幕亚洲| 久久精品aⅴ一区二区三区四区| 日日干狠狠操夜夜爽| www国产在线视频色| 国产伦一二天堂av在线观看| 国内揄拍国产精品人妻在线| 日韩高清综合在线| 天堂影院成人在线观看| 国产精品亚洲一级av第二区| 老汉色av国产亚洲站长工具| av女优亚洲男人天堂 | 在线观看美女被高潮喷水网站 | 午夜亚洲福利在线播放| 精品日产1卡2卡| 成人国产综合亚洲| 丰满人妻一区二区三区视频av | 亚洲成a人片在线一区二区| 91在线观看av| 国产av在哪里看| 99国产综合亚洲精品| 成年免费大片在线观看| 亚洲欧美日韩东京热| 亚洲成人精品中文字幕电影| a级毛片在线看网站| 高清在线国产一区| 我的老师免费观看完整版| 日本精品一区二区三区蜜桃| 亚洲欧美日韩高清专用| 成人精品一区二区免费| 国产真人三级小视频在线观看| 最近最新中文字幕大全免费视频| 日韩欧美国产在线观看| 欧美三级亚洲精品| 哪里可以看免费的av片| 中文字幕精品亚洲无线码一区| 99国产精品99久久久久| 午夜福利在线在线| 免费搜索国产男女视频| 淫妇啪啪啪对白视频| 午夜视频精品福利| 一夜夜www| 色综合欧美亚洲国产小说| 欧美不卡视频在线免费观看| 亚洲,欧美精品.| 毛片女人毛片| 丁香六月欧美| 三级国产精品欧美在线观看 | 欧美中文综合在线视频| 欧美高清成人免费视频www| 国产高清视频在线播放一区| 欧美日韩乱码在线| 一本精品99久久精品77| 一进一出抽搐gif免费好疼| 麻豆国产97在线/欧美| 热99re8久久精品国产| 久久久久久久久中文| 免费观看的影片在线观看| 日本一本二区三区精品| 好男人电影高清在线观看| 成在线人永久免费视频| 成人18禁在线播放| 国产爱豆传媒在线观看| 国产高清有码在线观看视频| 国产高清三级在线| 97碰自拍视频| 久久亚洲真实| 美女被艹到高潮喷水动态| 91字幕亚洲| 叶爱在线成人免费视频播放| 国产又黄又爽又无遮挡在线| 精品人妻1区二区| 一个人免费在线观看电影 | 99在线人妻在线中文字幕| 亚洲欧美日韩高清专用| 窝窝影院91人妻| 此物有八面人人有两片| 国产精品久久久久久人妻精品电影| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看| 国产精品九九99| 国产视频内射| 黑人巨大精品欧美一区二区mp4| 久久久久久久久久黄片| 夜夜夜夜夜久久久久| www日本在线高清视频| 国产精品av视频在线免费观看| www.www免费av| 中文字幕久久专区| www日本黄色视频网| 欧美xxxx黑人xx丫x性爽| 99热这里只有是精品50| 男人舔女人的私密视频| 久久久久久久精品吃奶| 黑人欧美特级aaaaaa片| 99国产极品粉嫩在线观看| 亚洲七黄色美女视频| 久久久久精品国产欧美久久久| 免费在线观看影片大全网站| 黄频高清免费视频| 99热这里只有是精品50| 少妇人妻一区二区三区视频| a级毛片在线看网站| 精品国产乱子伦一区二区三区| 又粗又爽又猛毛片免费看| 日本免费a在线| 成人18禁在线播放| 男人舔女人的私密视频| 亚洲精品中文字幕一二三四区| 国产精品 欧美亚洲| 国产极品精品免费视频能看的| 久久久久久大精品| 美女扒开内裤让男人捅视频| 色哟哟哟哟哟哟| 极品教师在线免费播放| 高清在线国产一区| 波多野结衣高清无吗| 亚洲精品中文字幕一二三四区| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av在线| 国产成人精品久久二区二区91| 观看免费一级毛片| 欧美日本亚洲视频在线播放| 一本一本综合久久| 亚洲自偷自拍图片 自拍| 九九久久精品国产亚洲av麻豆 | 亚洲精品一卡2卡三卡4卡5卡| 黄片小视频在线播放| 非洲黑人性xxxx精品又粗又长| 日韩精品中文字幕看吧| 九九热线精品视视频播放| 99久久综合精品五月天人人| 午夜福利在线观看免费完整高清在 | 亚洲在线观看片| 久久久久亚洲av毛片大全| 91久久精品国产一区二区成人 | 搡老熟女国产l中国老女人| 久久婷婷人人爽人人干人人爱| 麻豆成人午夜福利视频| 一个人看的www免费观看视频| 亚洲精品美女久久久久99蜜臀| 成年免费大片在线观看| 国产99白浆流出| 国产精品av久久久久免费| 亚洲第一电影网av| 村上凉子中文字幕在线| 在线观看日韩欧美| 国产精品国产高清国产av| 免费高清视频大片| 在线播放国产精品三级| 久久久国产欧美日韩av| 可以在线观看的亚洲视频| 国产成+人综合+亚洲专区| 色综合亚洲欧美另类图片| 国产成人一区二区三区免费视频网站| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 日本黄色片子视频| 成人特级黄色片久久久久久久| 最近视频中文字幕2019在线8| 蜜桃久久精品国产亚洲av| 亚洲成人免费电影在线观看| 亚洲av美国av| 国内毛片毛片毛片毛片毛片| 中文字幕人妻丝袜一区二区| 日韩av在线大香蕉| 成人亚洲精品av一区二区| 午夜福利在线观看免费完整高清在 | 噜噜噜噜噜久久久久久91| 无人区码免费观看不卡| 黑人操中国人逼视频| 啦啦啦免费观看视频1| 久久久久九九精品影院| 啦啦啦韩国在线观看视频| 无遮挡黄片免费观看| 日本 欧美在线| 国产高清三级在线| 熟女电影av网| 国产一区二区在线av高清观看| 在线看三级毛片| 婷婷亚洲欧美| 在线永久观看黄色视频| 国产亚洲av嫩草精品影院| 久久午夜综合久久蜜桃| 亚洲专区国产一区二区| 1024香蕉在线观看| 好男人在线观看高清免费视频| 少妇的逼水好多| 精品久久久久久,| 亚洲七黄色美女视频| 一卡2卡三卡四卡精品乱码亚洲| 手机成人av网站| 国产午夜精品久久久久久| 法律面前人人平等表现在哪些方面| 国产乱人视频| 麻豆国产av国片精品| 丁香欧美五月| 欧美日韩一级在线毛片| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件 | 成人特级av手机在线观看| 在线播放国产精品三级| 村上凉子中文字幕在线| 1024手机看黄色片| 亚洲国产高清在线一区二区三| 亚洲欧美一区二区三区黑人| 天堂av国产一区二区熟女人妻| 国产蜜桃级精品一区二区三区| 成年免费大片在线观看| 亚洲专区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 曰老女人黄片| 狠狠狠狠99中文字幕| 啦啦啦观看免费观看视频高清| 日本 av在线| 美女 人体艺术 gogo| 久久香蕉精品热| 久久久精品欧美日韩精品| 久久午夜综合久久蜜桃| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美精品v在线| av女优亚洲男人天堂 | 三级毛片av免费| 国产精品一区二区精品视频观看| 俄罗斯特黄特色一大片| 日日夜夜操网爽| 国产精品久久久久久精品电影| 国产精品亚洲av一区麻豆| ponron亚洲| 亚洲精品久久国产高清桃花| 网址你懂的国产日韩在线| 久久中文字幕一级| 国产亚洲精品av在线| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 精品久久久久久,| bbb黄色大片| 久久精品91蜜桃| cao死你这个sao货| 欧美日本视频| 露出奶头的视频| 看免费av毛片| 波多野结衣高清作品| 高清在线国产一区| 欧美在线黄色| 老司机福利观看| 最好的美女福利视频网| 18禁黄网站禁片午夜丰满| 午夜福利在线观看免费完整高清在 | 国产伦精品一区二区三区视频9 | 精品国产乱子伦一区二区三区| 九九久久精品国产亚洲av麻豆 | 99国产极品粉嫩在线观看| 99在线视频只有这里精品首页| 麻豆成人av在线观看| 亚洲av熟女| 色综合婷婷激情| 亚洲熟妇熟女久久| 高潮久久久久久久久久久不卡| 一进一出好大好爽视频| 国产av一区在线观看免费| 国产不卡一卡二| 色视频www国产| 嫩草影视91久久| 久9热在线精品视频| 九九在线视频观看精品| 啪啪无遮挡十八禁网站| 欧美不卡视频在线免费观看| 啪啪无遮挡十八禁网站| 一个人看的www免费观看视频| 一夜夜www| 中文资源天堂在线| 叶爱在线成人免费视频播放| 一本综合久久免费| 精品久久久久久久毛片微露脸| 欧美一级a爱片免费观看看| 欧美午夜高清在线| 中文字幕最新亚洲高清| 国产精品一及| 国产 一区 欧美 日韩| 国内揄拍国产精品人妻在线| 欧洲精品卡2卡3卡4卡5卡区| 久久婷婷人人爽人人干人人爱| 老司机福利观看| 日韩成人在线观看一区二区三区| 欧美中文日本在线观看视频| 91av网一区二区| 国产亚洲欧美98| 婷婷精品国产亚洲av在线| 可以在线观看的亚洲视频| 麻豆久久精品国产亚洲av| 国产亚洲精品综合一区在线观看| 中文字幕人成人乱码亚洲影| 久久久国产成人免费| 久久久久精品国产欧美久久久| 中文字幕久久专区| 少妇丰满av| av片东京热男人的天堂| 亚洲av成人一区二区三| 中文字幕人成人乱码亚洲影| 亚洲五月天丁香| 91九色精品人成在线观看| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 日韩国内少妇激情av| 亚洲欧美日韩无卡精品| 91久久精品国产一区二区成人 | 国内揄拍国产精品人妻在线| 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 久久婷婷人人爽人人干人人爱| 一本一本综合久久| 久久久久国产精品人妻aⅴ院| 国内揄拍国产精品人妻在线| 18禁黄网站禁片午夜丰满| 黄色成人免费大全| 久久久久久国产a免费观看| 国产成人啪精品午夜网站| 精品国产超薄肉色丝袜足j| 叶爱在线成人免费视频播放| 国产又黄又爽又无遮挡在线| 日本与韩国留学比较| 18禁裸乳无遮挡免费网站照片| 搡老熟女国产l中国老女人| 在线免费观看的www视频| 久久草成人影院| 999久久久精品免费观看国产| 国产伦人伦偷精品视频| 久久热在线av| 亚洲国产精品999在线| 久久午夜亚洲精品久久| 超碰成人久久| 叶爱在线成人免费视频播放| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 精品福利观看| 黄片大片在线免费观看| 日韩欧美 国产精品| 操出白浆在线播放| 人人妻人人澡欧美一区二区| 少妇的丰满在线观看| 日韩欧美精品v在线| 久久久久国产一级毛片高清牌| 亚洲国产精品999在线| 亚洲欧美精品综合久久99| 日韩欧美精品v在线| 精品久久蜜臀av无| 99久久精品热视频| 淫妇啪啪啪对白视频| 变态另类丝袜制服| 亚洲中文日韩欧美视频| 国产伦人伦偷精品视频| h日本视频在线播放| 欧美成人一区二区免费高清观看 | 国产麻豆成人av免费视频| 国产精品久久久久久亚洲av鲁大| 两人在一起打扑克的视频| 在线视频色国产色| 熟女少妇亚洲综合色aaa.| 亚洲成人精品中文字幕电影| 又爽又黄无遮挡网站| 18禁黄网站禁片免费观看直播| 精品午夜福利视频在线观看一区| 国产精品久久久久久精品电影| 最新美女视频免费是黄的| 美女高潮喷水抽搐中文字幕| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 国内毛片毛片毛片毛片毛片| 免费无遮挡裸体视频| 国模一区二区三区四区视频 | e午夜精品久久久久久久| 无遮挡黄片免费观看| 老司机福利观看| 99久久精品热视频| 一级毛片女人18水好多| 亚洲狠狠婷婷综合久久图片| 少妇熟女aⅴ在线视频| cao死你这个sao货| 欧美三级亚洲精品| 午夜激情福利司机影院| 啦啦啦免费观看视频1| 色尼玛亚洲综合影院| 国产精品免费一区二区三区在线| 国产精品电影一区二区三区| 国内毛片毛片毛片毛片毛片| 婷婷精品国产亚洲av| 男人的好看免费观看在线视频| 不卡一级毛片| 免费在线观看亚洲国产| 五月伊人婷婷丁香| 观看免费一级毛片| 九九热线精品视视频播放| 亚洲自偷自拍图片 自拍| 欧美黑人巨大hd| 亚洲电影在线观看av| 熟女人妻精品中文字幕| 成人国产一区最新在线观看| 在线免费观看不下载黄p国产 | 久久久久亚洲av毛片大全| 午夜精品久久久久久毛片777| 黄色成人免费大全| 观看美女的网站| 欧美高清成人免费视频www| 欧美一区二区国产精品久久精品| 人妻丰满熟妇av一区二区三区| 国产一区二区在线观看日韩 | 国产精品av视频在线免费观看| 午夜激情福利司机影院| 一区二区三区激情视频| 亚洲狠狠婷婷综合久久图片| 成人欧美大片| 99热只有精品国产| 日韩欧美国产在线观看| av欧美777| 色综合欧美亚洲国产小说| 国产野战对白在线观看| 日韩有码中文字幕| 日本黄大片高清| 男人的好看免费观看在线视频| 国产精品99久久99久久久不卡| а√天堂www在线а√下载| 在线观看舔阴道视频| 亚洲无线在线观看| 成人欧美大片| 老汉色av国产亚洲站长工具| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 精品日产1卡2卡| 国产精品久久久av美女十八| 成年免费大片在线观看| 日韩欧美国产一区二区入口| 熟女人妻精品中文字幕| 2021天堂中文幕一二区在线观| 久久久久国产一级毛片高清牌| 午夜免费观看网址| 成人精品一区二区免费| 免费人成视频x8x8入口观看| 精品国产超薄肉色丝袜足j| 精品一区二区三区四区五区乱码| 久久亚洲精品不卡| 一级毛片女人18水好多| 国产黄色小视频在线观看| 在线观看免费视频日本深夜| 国产成人av教育| 国产又色又爽无遮挡免费看| 人妻丰满熟妇av一区二区三区| 国产精品av久久久久免费| 亚洲av五月六月丁香网| 成人av在线播放网站| 丁香欧美五月| 午夜精品一区二区三区免费看| 叶爱在线成人免费视频播放| 小说图片视频综合网站| 国产美女午夜福利| 精品国产三级普通话版| 日日摸夜夜添夜夜添小说| 日日干狠狠操夜夜爽| 国产极品精品免费视频能看的| av片东京热男人的天堂| 国产私拍福利视频在线观看| 亚洲午夜理论影院| www.www免费av| 国产精品影院久久| 99久久成人亚洲精品观看| 亚洲中文av在线| 成人国产综合亚洲| av女优亚洲男人天堂 | 久久热在线av| 香蕉丝袜av| 夜夜看夜夜爽夜夜摸| av片东京热男人的天堂| 免费电影在线观看免费观看| 中文资源天堂在线| 男人的好看免费观看在线视频| 精品人妻1区二区| av黄色大香蕉| av在线天堂中文字幕| 日本成人三级电影网站| 在线看三级毛片| 免费av毛片视频| 99久久99久久久精品蜜桃| 午夜福利在线观看吧| 国产三级中文精品| 在线播放国产精品三级| 舔av片在线| 亚洲人成网站在线播放欧美日韩| 美女 人体艺术 gogo| 婷婷六月久久综合丁香| 日韩欧美在线乱码| 99riav亚洲国产免费| 国产男靠女视频免费网站| www.熟女人妻精品国产| 成在线人永久免费视频| 亚洲 欧美 日韩 在线 免费| 少妇的丰满在线观看| 88av欧美| 不卡av一区二区三区| 亚洲色图 男人天堂 中文字幕| 一级a爱片免费观看的视频| 日韩大尺度精品在线看网址| av视频在线观看入口| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| www.999成人在线观看| 老司机深夜福利视频在线观看| 男人和女人高潮做爰伦理| 一级a爱片免费观看的视频| 久久久精品大字幕| 久久久久久久精品吃奶| 又粗又爽又猛毛片免费看| www.熟女人妻精品国产| 99国产综合亚洲精品| 成年女人看的毛片在线观看| 99久久综合精品五月天人人| 国产一区二区在线av高清观看| h日本视频在线播放| 精品一区二区三区四区五区乱码| 宅男免费午夜| 两性午夜刺激爽爽歪歪视频在线观看| 男人和女人高潮做爰伦理| 国产一区二区在线观看日韩 | 成人国产综合亚洲| 99久久国产精品久久久| 1024香蕉在线观看| 超碰成人久久| 精品熟女少妇八av免费久了| 精华霜和精华液先用哪个| 午夜精品一区二区三区免费看| 国产成人精品无人区| 老熟妇仑乱视频hdxx| 国产精品 欧美亚洲| 亚洲狠狠婷婷综合久久图片| 精品不卡国产一区二区三区| 性欧美人与动物交配| www.999成人在线观看| 免费无遮挡裸体视频| 老司机午夜十八禁免费视频| 国产精品一及| 在线十欧美十亚洲十日本专区| 可以在线观看的亚洲视频| 国产欧美日韩一区二区三| 特级一级黄色大片| 精品一区二区三区视频在线 | 十八禁人妻一区二区| 97人妻精品一区二区三区麻豆| 日日摸夜夜添夜夜添小说| 亚洲精品久久国产高清桃花| 国产精品久久久av美女十八| 亚洲av中文字字幕乱码综合| 国产成人影院久久av| 美女免费视频网站| 窝窝影院91人妻| 日日干狠狠操夜夜爽| 此物有八面人人有两片| 久久人妻av系列| 日韩欧美免费精品| 久久久久久人人人人人| 亚洲欧洲精品一区二区精品久久久| 性色av乱码一区二区三区2| 日韩免费av在线播放| 国产av一区在线观看免费| 国产黄色小视频在线观看| 大型黄色视频在线免费观看| 最新美女视频免费是黄的| 亚洲欧美日韩卡通动漫| 亚洲精品在线观看二区| 色综合站精品国产| 国产99白浆流出| 色尼玛亚洲综合影院| 男女之事视频高清在线观看|