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

    亚洲男人天堂网一区| 一本久久精品| 在线观看人妻少妇| 最近中文字幕2019免费版| 欧美精品一区二区免费开放| 人人澡人人妻人| 黄色视频在线播放观看不卡| 日韩视频在线欧美| 最近最新中文字幕大全免费视频 | 久久99精品国语久久久| 国产片内射在线| 老司机影院成人| 女人高潮潮喷娇喘18禁视频| 夫妻性生交免费视频一级片| 免费在线观看黄色视频的| 亚洲情色 制服丝袜| a 毛片基地| 国产免费又黄又爽又色| 免费高清在线观看日韩| 日本vs欧美在线观看视频| 免费观看无遮挡的男女| 老汉色∧v一级毛片| 另类亚洲欧美激情| 欧美亚洲 丝袜 人妻 在线| 亚洲内射少妇av| 久久国内精品自在自线图片| 国产野战对白在线观看| 狠狠婷婷综合久久久久久88av| 国产精品秋霞免费鲁丝片| 久久久久久久精品精品| 三上悠亚av全集在线观看| 深夜精品福利| 黄频高清免费视频| 国产一区二区三区av在线| 国产精品一国产av| 免费看av在线观看网站| 自线自在国产av| 亚洲视频免费观看视频| 国产毛片在线视频| 叶爱在线成人免费视频播放| 超色免费av| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品94久久精品| 国语对白做爰xxxⅹ性视频网站| 麻豆精品久久久久久蜜桃| 国产乱人偷精品视频| 日韩电影二区| 国产毛片在线视频| 狠狠婷婷综合久久久久久88av| 老鸭窝网址在线观看| 在线观看免费日韩欧美大片| 国产一区二区三区av在线| 亚洲av.av天堂| 97在线视频观看| 久久精品久久精品一区二区三区| av电影中文网址| 2021少妇久久久久久久久久久| 秋霞在线观看毛片| 韩国精品一区二区三区| 国产又爽黄色视频| 精品少妇黑人巨大在线播放| 亚洲精品日韩在线中文字幕| 久久ye,这里只有精品| 免费黄色在线免费观看| 亚洲成av片中文字幕在线观看 | 国产人伦9x9x在线观看 | 中文字幕人妻熟女乱码| av电影中文网址| av在线老鸭窝| 亚洲国产精品一区二区三区在线| www.熟女人妻精品国产| 蜜桃国产av成人99| 国产精品av久久久久免费| 色婷婷av一区二区三区视频| 国产老妇伦熟女老妇高清| 国产日韩欧美亚洲二区| 男女免费视频国产| 亚洲精品国产色婷婷电影| 建设人人有责人人尽责人人享有的| 亚洲视频免费观看视频| 精品国产乱码久久久久久小说| 国产免费视频播放在线视频| 欧美av亚洲av综合av国产av | 中文精品一卡2卡3卡4更新| 欧美国产精品va在线观看不卡| 国产日韩欧美在线精品| freevideosex欧美| av国产精品久久久久影院| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 亚洲天堂av无毛| 亚洲色图综合在线观看| 国产亚洲精品第一综合不卡| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 韩国高清视频一区二区三区| 伦理电影大哥的女人| 成人国产麻豆网| 久久97久久精品| 狠狠精品人妻久久久久久综合| 国产白丝娇喘喷水9色精品| 精品一区二区三区四区五区乱码 | 一本色道久久久久久精品综合| 99久久综合免费| 日韩中文字幕欧美一区二区 | 亚洲av免费高清在线观看| 亚洲精品中文字幕在线视频| 国产精品一区二区在线观看99| 亚洲伊人久久精品综合| 亚洲国产色片| 欧美人与善性xxx| 日本色播在线视频| 亚洲av免费高清在线观看| 精品一品国产午夜福利视频| 亚洲精品视频女| 久久人人爽av亚洲精品天堂| 国产在线一区二区三区精| 丝袜美腿诱惑在线| 激情五月婷婷亚洲| 中文字幕制服av| 亚洲av电影在线进入| 久久久久国产一级毛片高清牌| 亚洲四区av| 日韩欧美一区视频在线观看| 哪个播放器可以免费观看大片| 国产av国产精品国产| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 激情五月婷婷亚洲| av线在线观看网站| 国产成人精品久久二区二区91 | 久久国产亚洲av麻豆专区| 精品亚洲成a人片在线观看| 亚洲第一区二区三区不卡| 欧美日韩亚洲国产一区二区在线观看 | 久久人妻熟女aⅴ| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 极品人妻少妇av视频| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 日韩一卡2卡3卡4卡2021年| 中文字幕精品免费在线观看视频| 欧美+日韩+精品| 又大又黄又爽视频免费| 亚洲国产色片| 90打野战视频偷拍视频| 黑人猛操日本美女一级片| 蜜桃国产av成人99| 最近最新中文字幕大全免费视频 | 青草久久国产| 国产片内射在线| 精品久久久精品久久久| 热re99久久精品国产66热6| 久久狼人影院| 99国产综合亚洲精品| 最新中文字幕久久久久| 国产一区二区 视频在线| 美女国产高潮福利片在线看| 亚洲国产精品国产精品| 黄色一级大片看看| 亚洲精品自拍成人| 精品午夜福利在线看| 男男h啪啪无遮挡| 久久久久久久大尺度免费视频| 高清视频免费观看一区二区| 久久99蜜桃精品久久| 欧美少妇被猛烈插入视频| 飞空精品影院首页| 人成视频在线观看免费观看| 久久 成人 亚洲| 波多野结衣av一区二区av| 精品人妻偷拍中文字幕| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| 不卡av一区二区三区| kizo精华| 一本—道久久a久久精品蜜桃钙片| 国产精品99久久99久久久不卡 | 又大又黄又爽视频免费| 欧美成人午夜精品| 看十八女毛片水多多多| 免费av中文字幕在线| 18禁观看日本| 亚洲国产精品一区二区三区在线| 免费观看无遮挡的男女| 亚洲精华国产精华液的使用体验| 一级毛片 在线播放| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 成年女人毛片免费观看观看9 | 91成人精品电影| 91午夜精品亚洲一区二区三区| 色视频在线一区二区三区| 日本91视频免费播放| 国产精品国产三级专区第一集| 99久久综合免费| 日韩欧美一区视频在线观看| 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 免费女性裸体啪啪无遮挡网站| 成人影院久久| 成年av动漫网址| 国产av码专区亚洲av| 欧美精品av麻豆av| 丝瓜视频免费看黄片| 人成视频在线观看免费观看| videossex国产| 免费看不卡的av| 视频区图区小说| 中文乱码字字幕精品一区二区三区| 亚洲欧美精品综合一区二区三区 | 欧美精品国产亚洲| 女人被躁到高潮嗷嗷叫费观| 国产极品天堂在线| 国产精品嫩草影院av在线观看| 久久av网站| 国产成人免费观看mmmm| 亚洲精品国产av蜜桃| 午夜福利网站1000一区二区三区| 热re99久久国产66热| 一区二区日韩欧美中文字幕| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 电影成人av| 国产成人免费无遮挡视频| 精品亚洲乱码少妇综合久久| 街头女战士在线观看网站| 国产在视频线精品| av有码第一页| 夫妻午夜视频| 免费少妇av软件| 一区在线观看完整版| 下体分泌物呈黄色| 人人妻人人澡人人爽人人夜夜| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 在线天堂中文资源库| 又粗又硬又长又爽又黄的视频| 美女视频免费永久观看网站| 精品人妻偷拍中文字幕| 久久久久网色| 香蕉丝袜av| 久久久久久久久免费视频了| videos熟女内射| 只有这里有精品99| 久久久久精品性色| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 欧美老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 9191精品国产免费久久| 日本黄色日本黄色录像| 搡老乐熟女国产| 亚洲欧美色中文字幕在线| 欧美日韩一级在线毛片| 在线观看www视频免费| 国产成人精品久久二区二区91 | 久久久久网色| av国产精品久久久久影院| 午夜福利一区二区在线看| 亚洲美女视频黄频| 另类精品久久| 欧美亚洲 丝袜 人妻 在线| 在现免费观看毛片| 中文字幕人妻丝袜制服| 毛片一级片免费看久久久久| 91精品三级在线观看| 又粗又硬又长又爽又黄的视频| 久久久久久久国产电影| 日本午夜av视频| 亚洲欧美一区二区三区国产| 亚洲欧美一区二区三区黑人 | 男女下面插进去视频免费观看| 啦啦啦在线免费观看视频4| 久久久欧美国产精品| 精品久久蜜臀av无| 国产日韩欧美亚洲二区| 日韩制服骚丝袜av| 国产精品久久久久久精品古装| 久久久久久久久久久久大奶| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 少妇猛男粗大的猛烈进出视频| 亚洲国产精品999| 亚洲内射少妇av| 91成人精品电影| 26uuu在线亚洲综合色| 欧美成人午夜精品| 麻豆av在线久日| 欧美国产精品va在线观看不卡| 亚洲成色77777| 日韩成人av中文字幕在线观看| av视频免费观看在线观看| 成年女人毛片免费观看观看9 | 在线观看免费高清a一片| 国产不卡av网站在线观看| 最近最新中文字幕免费大全7| 国产亚洲一区二区精品| 国产黄色免费在线视频| 在线免费观看不下载黄p国产| 中文字幕人妻丝袜制服| 午夜日韩欧美国产| 国产精品.久久久| 天堂8中文在线网| 免费看不卡的av| 亚洲内射少妇av| 欧美av亚洲av综合av国产av | 欧美日韩国产mv在线观看视频| 日韩,欧美,国产一区二区三区| av天堂久久9| 秋霞伦理黄片| 男女边吃奶边做爰视频| 两个人免费观看高清视频| 亚洲伊人久久精品综合| 熟妇人妻不卡中文字幕| 亚洲欧美成人综合另类久久久| 久久精品久久久久久久性| 亚洲av中文av极速乱| 亚洲欧美中文字幕日韩二区| 成人毛片a级毛片在线播放| 超碰成人久久| 大话2 男鬼变身卡| 男人爽女人下面视频在线观看| 最近中文字幕2019免费版| 国精品久久久久久国模美| 精品人妻熟女毛片av久久网站| 男女国产视频网站| 女人高潮潮喷娇喘18禁视频| 欧美精品一区二区大全| 91在线精品国自产拍蜜月| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产无遮挡羞羞视频在线观看| 成年美女黄网站色视频大全免费| 青青草视频在线视频观看| 久久久久久久亚洲中文字幕| 欧美成人午夜免费资源| 老汉色av国产亚洲站长工具| 啦啦啦在线观看免费高清www| 亚洲国产欧美日韩在线播放| 男的添女的下面高潮视频| 精品卡一卡二卡四卡免费| 国产精品国产av在线观看| 女人精品久久久久毛片| 久久久久国产精品人妻一区二区| 欧美+日韩+精品| 最近最新中文字幕免费大全7| 国产精品久久久久久精品古装| 国产成人精品久久久久久| 毛片一级片免费看久久久久| 免费少妇av软件| 最新的欧美精品一区二区| 春色校园在线视频观看| 亚洲男人天堂网一区| 日本av免费视频播放| 日本欧美国产在线视频| 亚洲,欧美,日韩| 哪个播放器可以免费观看大片| 激情视频va一区二区三区| 大香蕉久久成人网| 亚洲成人一二三区av| 亚洲美女黄色视频免费看| 丰满饥渴人妻一区二区三| 青草久久国产| 18禁动态无遮挡网站| 哪个播放器可以免费观看大片| 婷婷成人精品国产| 18在线观看网站| 男人爽女人下面视频在线观看| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 精品久久久久久电影网| 久久精品夜色国产| 国产在线免费精品| 国产亚洲最大av| 久久久a久久爽久久v久久| 国产 一区精品| 亚洲欧美一区二区三区国产| 久久精品人人爽人人爽视色| 人人澡人人妻人| 精品人妻在线不人妻| 菩萨蛮人人尽说江南好唐韦庄| 午夜激情av网站| 在线天堂最新版资源| 最近中文字幕高清免费大全6| 久久久欧美国产精品| 爱豆传媒免费全集在线观看| 久久久久久久国产电影| 国产黄频视频在线观看| 免费日韩欧美在线观看| 久久女婷五月综合色啪小说| 97在线视频观看| 久久久久久久大尺度免费视频| 国产女主播在线喷水免费视频网站| 国产1区2区3区精品| 黄色配什么色好看| 好男人视频免费观看在线| 最近中文字幕2019免费版| 国产精品久久久久成人av| 久久精品国产亚洲av高清一级| 2021少妇久久久久久久久久久| 男女免费视频国产| 欧美日韩精品网址| 久久精品亚洲av国产电影网| 纯流量卡能插随身wifi吗| 国产麻豆69| 日本91视频免费播放| 大码成人一级视频| 国产成人午夜福利电影在线观看| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 欧美日韩国产mv在线观看视频| 岛国毛片在线播放| 国产人伦9x9x在线观看 | 一级毛片 在线播放| av在线app专区| av又黄又爽大尺度在线免费看| 日韩精品免费视频一区二区三区| av天堂久久9| 精品少妇黑人巨大在线播放| 精品卡一卡二卡四卡免费| 中文字幕人妻丝袜制服| 国产精品久久久久久精品电影小说| 中文天堂在线官网| 超碰97精品在线观看| 欧美日韩一级在线毛片| 免费黄频网站在线观看国产| 男女边吃奶边做爰视频| 免费在线观看完整版高清| 啦啦啦中文免费视频观看日本| 成人亚洲精品一区在线观看| 亚洲av免费高清在线观看| 中文天堂在线官网| 日韩中文字幕欧美一区二区 | 男的添女的下面高潮视频| 久久人人爽人人片av| 日韩av免费高清视频| 免费播放大片免费观看视频在线观看| 青春草国产在线视频| 男女啪啪激烈高潮av片| 九色亚洲精品在线播放| 制服丝袜香蕉在线| 一区二区三区乱码不卡18| 亚洲第一青青草原| 久久久国产一区二区| 国产亚洲精品第一综合不卡| 国产精品无大码| 国产精品蜜桃在线观看| av片东京热男人的天堂| 男人添女人高潮全过程视频| 9色porny在线观看| 性高湖久久久久久久久免费观看| videossex国产| 国产成人精品无人区| 最近2019中文字幕mv第一页| 日韩中字成人| 国产一区有黄有色的免费视频| 一级毛片电影观看| 岛国毛片在线播放| 99久久精品国产国产毛片| 亚洲激情五月婷婷啪啪| 丝袜在线中文字幕| 一本色道久久久久久精品综合| 精品国产一区二区三区四区第35| 美女xxoo啪啪120秒动态图| 女性生殖器流出的白浆| www日本在线高清视频| 最新的欧美精品一区二区| 亚洲精品国产av蜜桃| 免费不卡的大黄色大毛片视频在线观看| 人妻 亚洲 视频| 最近中文字幕高清免费大全6| 国产成人精品在线电影| 人妻少妇偷人精品九色| 五月伊人婷婷丁香| 午夜福利视频在线观看免费| 国产精品 欧美亚洲| 久久亚洲国产成人精品v| 毛片一级片免费看久久久久| 边亲边吃奶的免费视频| 婷婷色av中文字幕| 美女国产视频在线观看| 一区二区日韩欧美中文字幕| av国产精品久久久久影院| 日本免费在线观看一区| 丝袜脚勾引网站| 如日韩欧美国产精品一区二区三区| 香蕉精品网在线| 香蕉国产在线看| 热re99久久精品国产66热6| 永久免费av网站大全| 一本大道久久a久久精品| 天堂中文最新版在线下载| 国产精品久久久久久精品电影小说| 中文天堂在线官网| 久久人人爽人人片av| 91国产中文字幕| 久热久热在线精品观看| 成人手机av| 久久人人爽人人片av| 电影成人av| 水蜜桃什么品种好| 国产女主播在线喷水免费视频网站| 国产精品久久久久成人av| 久久人妻熟女aⅴ| 国产精品久久久久成人av| 中文字幕人妻丝袜制服| 一本大道久久a久久精品| 欧美日韩综合久久久久久| 丰满少妇做爰视频| 国产av精品麻豆| 人妻人人澡人人爽人人| 国产精品av久久久久免费| 日韩欧美精品免费久久| 黑丝袜美女国产一区| 亚洲国产欧美网| 亚洲成人av在线免费| 欧美日韩国产mv在线观看视频| 99国产综合亚洲精品| 亚洲综合精品二区| 人体艺术视频欧美日本| av卡一久久| 精品国产一区二区三区久久久樱花| 国语对白做爰xxxⅹ性视频网站| 中文字幕另类日韩欧美亚洲嫩草| 高清欧美精品videossex| 最近最新中文字幕大全免费视频 | 中文精品一卡2卡3卡4更新| 高清不卡的av网站| 国产男女超爽视频在线观看| 国产精品免费视频内射| 国产精品秋霞免费鲁丝片| 80岁老熟妇乱子伦牲交| 王馨瑶露胸无遮挡在线观看| 最新的欧美精品一区二区| 亚洲伊人色综图| 久久久久人妻精品一区果冻| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 一本—道久久a久久精品蜜桃钙片| 亚洲第一青青草原| 咕卡用的链子| 亚洲精品,欧美精品| 在线观看人妻少妇| 亚洲情色 制服丝袜| 天堂中文最新版在线下载| 蜜桃国产av成人99| 熟女少妇亚洲综合色aaa.| 美女高潮到喷水免费观看| 午夜激情av网站| 亚洲国产av新网站| 熟女av电影| 亚洲综合精品二区| 国产精品一国产av| 青春草视频在线免费观看| 国产精品国产av在线观看| 亚洲,一卡二卡三卡| 欧美精品一区二区大全| 成年av动漫网址| 国产视频首页在线观看| 国产熟女欧美一区二区| 丝袜在线中文字幕| 国产激情久久老熟女| 久久久久久久久久久免费av| 欧美成人午夜精品| 各种免费的搞黄视频| 日本黄色日本黄色录像| 国产精品99久久99久久久不卡 | 日韩成人av中文字幕在线观看| 满18在线观看网站| 另类精品久久| 日韩av免费高清视频| 久久精品人人爽人人爽视色| 日韩一区二区三区影片| 日韩中文字幕视频在线看片| 美女高潮到喷水免费观看| 最近的中文字幕免费完整| 色吧在线观看| 91精品三级在线观看| 下体分泌物呈黄色| 亚洲天堂av无毛| 一区福利在线观看| 亚洲国产最新在线播放| 久久免费观看电影| 日本91视频免费播放| 午夜免费鲁丝| 久久久久久久久久久久大奶| 亚洲图色成人| 美女xxoo啪啪120秒动态图| 中文字幕人妻丝袜制服| 女的被弄到高潮叫床怎么办| www.熟女人妻精品国产| 香蕉丝袜av| 亚洲精品国产av蜜桃| 亚洲欧洲精品一区二区精品久久久 | 国产又爽黄色视频| 少妇猛男粗大的猛烈进出视频| 欧美亚洲 丝袜 人妻 在线| 国产又爽黄色视频| 99久久综合免费| 超色免费av| 久久狼人影院| 午夜日韩欧美国产| 制服诱惑二区| 赤兔流量卡办理| 中文字幕精品免费在线观看视频| 国产欧美日韩综合在线一区二区| 亚洲精华国产精华液的使用体验| www日本在线高清视频|