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

    Fabrication and cold test of prototype of spatially periodic radio frequency quadrupole focusing linac

    2021-03-03 10:08:46PeiYanYuBinZhangFengFengWangChenXingLiGuoZhenSunZhiJunWangLuBeiLiuChenZhangYuanYuanHeHuShanXu
    Nuclear Science and Techniques 2021年1期

    Pei-Yan Yu? Bin Zhang ? Feng-Feng Wang? Chen-Xing Li ?Guo-Zhen Sun ? Zhi-Jun Wang ? Lu-Bei Liu ? Chen-Zhang Yuan ?Yuan He? Hu-Shan Xu

    Abstract A 325 MHz aluminum prototype of a spatially periodic RF quadrupole focusing linac was developed at the Institute of Modern Physics, Chinese Academy of Sciences, as a promising candidate for the front end of a high-current linac. It consists of an alternating series of crossbar H-type drift tubes and RF quadrupole sections.Owing to its special geometry,cavity fabrication is a major hurdle for its engineering development and application. In this paper,we report the detailed mechanical design of this structure and describe its fabrication process, including machining,assembly,and inspection.The field distribution was measured by the bead-pull technique.The results show that the field errors of both the accelerating and focusing fields are within an acceptable range. A tuning scheme for this new structure is proposed and verified. The cold test process and results are presented in detail. The development of this prototype provides valuable guidance for the application of the spatially periodic RF quadrupole structure.

    Keywords Spatially periodic RF quadrupole focusing linac · Mechanical structure design · Bead-pull measurement

    1 Introduction

    Accelerating structures with electric field focusing are routinely used in the front end of high-current proton or heavy-ion linacs. For example, the spatially uniform RF quadrupole (RFQ) focusing structure has become the optimal choice for use after the ion source because it integrates beam acceleration, bunching, and transverse focusing in one compact cavity [1, 2]. However, the accelerating efficiency of the RFQ structure decreases as the beam energy increases because the longitudinal electric field is provided only by vane modulation. Therefore, the RFQ structure is generally used for beam energies ranging from tens of keV to several MeV.

    To extend the application range of the RFQ focusing principle, several more efficient accelerating structures have been developed as alternatives to the higher-energy section of an RFQ linac and its subsequent accelerating structures; these structures are obtained by combining the RFQs with accelerating gaps formed between drift tubes.An early structure proposed by Teplyakov[3]supplements the accelerating gap with four conducting fingers arranged in a quadrupole geometry, which generate a focusing field in the accelerating gap. Structures based on this principle have been employed for many years in a linac injector for the proton synchrotron at the Institute for High Energy Physics [4]. A comparable scheme is the RF-focused drift tube(RFD)linac structure,in which the RFQ focusing field is placed inside a drift tube by splitting the drift tube into two independent electrodes that support two fingers each[5]. As a modified version of the RFD linac, the RF-focused interdigital linac structure adopts a combination of the interdigital H-mode structure and the RFQ focusing structure to reduce the cavity size and improve the shunt impedance [6].

    In recent years, new structures with spatially periodic RFQ focusing have been developed, which are characterized by the presence of separate accelerating and focusing zones. This concept was first suggested in the hybrid RFQ for the acceleration of low-charge-state heavy-ion beams[7, 8]. In addition, a structure suitable for all ion species was proposed; it consists of an alternating series of crossbar H-type drift tubes and four-vane-type RFQ sections.A detailed beam dynamics study showed that this structure can provide a high energy gain rate and high shunt impedance with moderate transverse emittance growth [9]. The electrode geometry is somewhat simpler than the focusing finger geometry; it avoids an overly high surface electric field and reduces the difficulty of RF tuning. In particular,this structure can be built using the conventional four-vanetype RFQ cavity,which makes it easier to cool and affords the potential for continuous-wave operation.

    The spatially periodic RFQ structure is a promising choice for the design of a compact and efficient front end of a high-current linac. Consequently, the research and development of this structure was begun at the Institute of Modern Physics, Chinese Academy of Sciences, as a prestudy of the High Intensity Heavy-Ion Accelerator Facility(HIAF) project. The HIAF is being designed to provide intense primary and radioactive ion beams for nuclear physics,atomic physics,and application research[10].The Ion-Linac used as the injector in the HIAF is a superconducting (SC) heavy-ion accelerator; it consists of a lowenergy beam transport line, an RFQ, a medium-energy beam transport line, and an SC section [11, 12]. The spatially periodic RFQ structure is intended to serve as an intermediate structure between the RFQ and SC cavities.It can be used with a shorter RFQ and improve the input energy of the SC section. As a key step in the study, an aluminum prototype was designed and fabricated to investigate the feasibility of the mechanical design and process technology for the new structure and also to determine the field distribution and its tunability. In this paper,the mechanical design,RF measurement,and tuning of the prototype are presented in detail.

    2 RF structure

    The RF structure of the prototype cavity is based on that of a conventional four-vane-type RFQ in which the tips of the vanes are cut off,and the remaining parts make up four girders. Drift tubes connected to two stems are installed crosswise on the girders to form the drift tube linac (DTL)sections. The RFQ sections are built by mounting the electrode blocks on the girders at the designed positions.Figure 1a shows the configuration of the alternating DTL and RFQ sections of the prototype cavity.

    The cavity is operated in an H210-like mode, and no dipole modes exist near the operating mode because the drift tubes cause the opposite girders to be electrically shorted [9, 13]. The longitudinal electric field is plotted in Fig.1b, which shows that charged particle acceleration will take place in both the gaps between the conventional drift tubes and the gaps between the drift tubes and RFQs.The RFQs provide transverse focusing for the beam and can also provide acceleration if the electrodes are modulated. The main parameters of the RF model are summarized in Table 1.The frequency of the prototype cavity was set to 325 MHz to keep its size relatively small.The cavity has 14 drift tubes(including two in the endplates) and two unmodulated RFQs.Twenty slug tuners with a diameter of 40 mm were arranged in four quadrants to study the tunability of the cavity. The cells were set to be identical in length to obtain a naturally flat field distribution and to simplify fabrication.

    3 Mechanical structure design and fabrication

    3.1 Mechanical design

    As shown in Fig.2a,the main structure of the prototype cavity is a box resembling that employed for the RFQ of the C-ADS Injector II [14]. The external structure of the cavity is divided into four cavity walls and two endplates.The four girders are integrated with the four corresponding cavity walls.Inside the cavity,the drift tubes and RFQs are distributed alternately on top of the girders.Each drift tube is connected with girders on either side by two collinear stems, and the stems of adjacent drift tubes are perpendicular to each other.

    The cavity is an overconstrained structure owing to the configuration of the stems inside the cavity; that is, one or more degrees of freedom (DOFs) of the components are repeatedly constrained by different fixtures. In the mechanical design process,a drift tube connected to a stem at each end is referred to as a drift tube component to avoid confusion with individual drift tubes. To install one drift tube component,the stems on each end have to be inserted into the installation holes on each symmetrical girder, as shown on the left side of Fig.2b. The DOFs of the drift tube component, except for the translational DOF in the axial direction of the stem, are overconstrained by the two installation holes. In turn, each cavity wall is subject to redundant constraints from the drift tube components,as it is fully constrained by the other cavity walls and endplates.The RFQ section is formed by inserting electrode blocks into the installation holes of the girders. Consequently, thesingle-electrode block is underconstrained because the translational DOF in the installation direction is free. The fabrication of the overconstrained structure requires highprecision machining because insufficient machining precision tends to cause assembly failure.

    Fig.1 (Color online) a RF structure of the prototype,b distribution of the accelerating field along the central axis of the cavity with 1 J of stored energy

    Table 1 Main parameters of the RF model

    One of the major differences between the proposed prototype cavity and the conventional DTL cavity is that the stems are positioned on the girders instead of on the outer cavity walls, which means that the drift tube component cannot be operated manually in certain assembly steps. During the assembly of the third or fourth cavity wall,the six installation holes on the girder must be fitted with the six drift tube components mounted on the opposite girder. Even if all the critical dimensions are guaranteed to be within tolerance, there is still a high risk of assembly failure without manual operation.

    Owing to the mechanical structure of the prototype cavity, threaded connections are added to the drift tube components and the electrode blocks,as shown on the right side of Fig.2b. A groove is designed on the cavity wall to enable the use of short, stable bolts. The remaining free DOF of the electrode block is constrained by the bolt. For the drift tube components, the constraint on the translational DOF in the axial direction of the stem is transferred from the bottom of the installation holes to the groove on the cavity wall.In addition,the overconstraint on the cavity wall in this direction is removed. Furthermore, this design makes it possible to manually control the drift tube components during the assembly process, thereby reducing the risk of assembly failure.

    Fig.2 (Color online)a Mechanical model of the prototype cavity, b illustration of the positioning of the drift tube component without threaded connection (left) and with threaded connection(right)

    3.2 Fabrication

    In the processing of the prototype cavity, several methods are used to control the machining precision of the parts. The concentricity of the stems in the drift tube component is a crucial parameter and strongly affects the overall assembly precision of the cavity. Therefore, a unique processing method is proposed to fabricate the drift tube components, as shown in Fig.3a. First, the external cylindrical surface of the drift tube is machined to its design size, and a radial hole for the connection with the stem and an axial hole as the center reference for subsequent processing are drilled. Next, the stem is penetrated through the drift tube using cold assembly technology.Finally, the internal cylindrical surface of the drift tube is machined to its design size; during this process, the long stem is split into two parts. The magnitude of the interference is set to 0.08 mm to ensure that the connection between the drift tube and the stem is reliable. Coordinate measuring machine results show that the concentricity of the stems in the drift tube component can reach φ0.05 mm.

    The machining of the aluminum alloy is likely to generate intense local stress and deformation, which would cause the workpiece to be scrapped for being out of tolerance. To reduce the machining deformation, the milling of the cavity walls and endplates proceeded as follows:rough machining, semi-finishing, and finishing. To eliminate the residual stress, heat treatment was applied after rough machining, and natural aging was performed after the other process steps.

    3.3 Assembly and inspection

    The assembly of the prototype cavity was performed under the guidance of the in-process inspection with a coordinate measuring arm [15]. Before the assembly process, the mating dimensions of the four cavity walls were inspected individually. The results showed that only the length of the cavity wall is out of tolerance owing to thermal shrinkage, which is relatively large (- 0.3 to- 0.4 mm), but would not hinder the assembly process.Then the drift tube components were mounted on the corresponding girder, the concentricity and spacing dimensions of which were measured. The measurement results can provide a preliminary check on the fitting accuracy and verify the feasibility of assembly.

    In the assembly process, the bottom and left cavity walls, each of which was equipped with six drift tubecomponents, were initially assembled, and the right cavity wall was installed (see Fig.3b). When the girder of the right cavity wall approached the drift tube components mounted on the opposite girder, a thread rod was used to pull the drift tube components out to align with both girders.After all six components were mounted on the two girders, the right cavity wall was slowly pushed into the design position, and the three assembled cavity walls formed a U-type subassembly. The top cavity wall was assembled by the same process. Figure 3c shows the assembled cavity without the front endplate.

    Fig.3 (Color online)Fabrication and assembly process: a processing of drift tube component, b assembly process of right cavity wall,c assembled cavity without the front endplate

    The fabrication precision of the cavity determines the concentricity of the drift tubes, which is a critical parameter for efficient beam transmission and is hard to measure for the finished cavity,because the probe of the coordinate measuring arm cannot reach the drift tubes in the middle of the cavity. Therefore, the concentricity of the drift tubes was initially measured on the U-type subassembly,and the results are shown in Fig.4. Note that the measured results differ slightly from those of the complete assembly owing to the absence of the top cavity wall. In the U-type subassembly, the measured positions of the vertical drift tube components (drift tubes 1, 3, 6, 8, 9, and 11) are slightly lower than the design positions because these components are bolted at the bottom and unconstrained at the top.Conversely, the positions of the horizontal drift tube components (drift tubes 2, 4, 5, 7, 10, and 12) are slightly higher than the design positions because the left and right cavity walls are not precisely parallel but have a small open angle under the constrained condition.

    Fig.4 Measurement results of the drift tube(DT)positions using the central axis as a reference

    A scheme for measuring the concentricity of the drift tubes in the finished cavity using a laser tracker was proposed. A rod with a magnet attached to the front end was designed that can fix a ferromagnetic spherically mounted retroreflector (SMR) by the attractive magnetic force. By using this rod, the SMR could reach the surfaces of all the drift tubes to perform the measurement with the laser tracker. The measurement results are shown in Fig.4.Taking the central axis defined by the cavity’s internal surfaces as a reference,the concentricity of the twelve drift tubes is φ0.22 mm; that is, the maximum deviation from the central axis is 0.11 mm. This result will be taken into account in the error study of the beam dynamics design,and the technical process will be optimized further to improve the precision.

    4 RF measurement

    4.1 Experimental setup

    In our early measurements, we found many errors caused by the noise floor of the VNA at the bottom of the field distribution curve.This phenomenon can be explained as follows. First, the relatively low Q value of the aluminum cavity results in a wider bandwidth,and more noise is collected during the measurement.In addition,according to the relationship between the electric field amplitude and the frequency shift, the noise will give rise to a larger deviation in the low-field areas, for example, in the drift tubes [18].

    To improve the measurement accuracy, the noise floor of the VNA was reduced by combining the averaging and the intermediate-frequency bandwidth; these two noise reduction features of the VNA are commonly used to address the measurement issue mentioned above [19].However, the use of both features prolongs the measurement time,during which environmental parameters such as temperature and humidity might change,introducing errors to the measurements [20]. Therefore, preliminary tests were performed to determine an appropriate setting of the VNA at which both the noise floor power and the measurement time are acceptable. Several dielectric spheres with different diameters were used as the perturbation bodies in these tests. It was found that the relative amplitude of the errors decreased with increasing sphere diameter, because a larger sphere would produce a larger frequency shift. However, this finding does not mean that the perturbation body should be as large as possible. The reason will be explained in Sect. 4.2.1. The perturbation body chosen for the field distribution measurement is a dielectric sphere with a diameter of 6 mm, and the maximum frequency shift it produces is approximately 7 kHz.

    Fig.5 (Color online) Photograph of the bead-pull measurement apparatus

    4.2 Results of RF measurement

    The frequency, the Q value, and the distribution of the accelerating and focusing fields of the prototype cavity were measured with all the tuners adjusted to the design insertion depth of 20 mm. The measured frequency was 326.72 MHz, which is 1.72 MHz higher than the design value. This deviation is attributed mainly to the errors in the frequency solution and fabrication errors. Thermal shrinkage makes the major contribution to the fabrication errors because the cavity was not processed in a constant temperature environment.The measured Q value was 6061 when no RF seals were used,which is 65%of the simulated value.

    4.2.1 Distribution of the accelerating field

    Figure 6a compares the measured and simulated field distributions along the central axis of the cavity, which show good agreement overall. The field amplitudes in the gaps between the internal drift tubes are flat. The normalized amplitude reaches 0.6 in the gaps at both ends of the cavity, and it is approximately 0.5 in the gaps between the drift tubes and the RFQs.To quantify the uniformity of the overall accelerating field, the amplitude deviation in the nine gaps between the internal drift tubes is defined as the field unevenness, which is 1.3%.

    In addition, inconsistencies appear between the measured and simulated results in the two RFQ sections.According to the simulated results, the amplitude should decrease to zero in the middle of the RFQ section,whereas the measured value at this position is 0.22. To explain this deviation, the simulated distribution of the transverse field at this position is investigated, as shown in Fig.6b. The amplitude along the angle bisector of the quadrants is zero at the center, but it increases linearly with increasing distance from the center. The Slater perturbation theorem assumes that the unperturbed field at the measurement position is homogeneous [17]. In our experiment, the measured field amplitude corresponds to the average amplitude over the 6 mm sphere, which apparently is not zero. The field amplitude at this position measured using spheres of different diameters is plotted using the preliminary test data in Fig.6c. The amplitude has a linear relationship with the sphere diameter, and the fitting line passes through the origin, confirming that the above explanation is correct.

    4.2.2 Distribution of the focusing field

    During the measurement of the focusing field, the distance between the trajectory of the bead and the cavity axis was set to 5 mm, because the inner diameter of the drift tube is 20 mm.The field distribution along the entire cavity was measured, and the field curves corresponding to the two RFQs were extracted from the measured curves of the entire cavity. The measurement results are summarized in Fig.6d and compared with the simulation results. The amplitudes are normalized by the maximum value over the entire trajectory of the bead, and the normalization coefficient is essentially the same as that of the accelerating field. Owing to the inhomogeneity of the nearby field, the measured amplitudes in the four quadrants are all higher than the simulated values.

    For a conventional RFQ, the frequency of the operating quadrupole mode is close to that of the dipole modes, and thus, accidental degeneracy between them is possible. The amplitudes of the dipole modes have to be minimized during the tuning process to produce the best approximation to a pure quadrupole mode [21]. For the spatially periodic RFQ structure, no dipole mode exists near the operating quadrupole mode because of the stem array.The amplitude of the quadrupole mode Uqis defined as the average of the measured amplitudes U1, U2, U3, and U4in the four quadrants. To evaluate the uniformity of the quadrupole field, the relative error of the quadrupole field is defined as the shape similarity error between the measured and simulated curves of Uq. The similarity of the curves is compared using the metrology software of the coordinate measurement machine. The relative errors of the quadrupole fields of the two RFQ sections are within ± 0.4% and ± 1.0%, respectively.

    Fig.6 (Color online) RF measurement results: a field distribution along the central axis of the cavity (normalized by the maximum value), b simulated field distribution along the angle bisector of the quadrants in the middle plane of the RFQ when the stored energy of the cavity is 1 J, c amplitude at the center of the first RFQ section measured using spheres of different diameters, d measured field distribution in the four quadrants and the measured and simulated field distribution of the quadrupole field (top: first section, bottom:second section)

    Distinct differences in the field amplitudes of the four quadrants are shown in Fig.6d.Specifically,U1and U2are relatively low, whereas U3and U4are relatively high. The differences are attributed mainly to the effects of gravity,which makes the bead sag from the desired position during the measurement [22]. Fortunately, the differences have little effect on Uqbecause the deviations of the four quadrants’ results caused by gravity cancel each other out in the calculation process.

    5 Tuning experiment

    5.1 Tuning algorithm

    Although the measurement results show that the field errors of both the accelerating and focusing fields are within a small and acceptable range, a tuning experiment aimed at improving the field uniformity was also performed to investigate the tunability of the cavity. A tuning algorithm based on a response matrix is employed for this new structure, which is built assuming a linear relationship between each tuner setting and the field amplitude at any given longitudinal position[23].It can be expressed by the following formula:

    where Vdiis the desired field amplitude at the ith measurement position and V0iis the actual field amplitude obtained from the measurement. Tdjand T0jrepresent the unknown desired setting and the actual setting of the jth tuner, respectively. Each ?Vi/?Tjvalue is the measured derivative of the ith amplitude with respect to the jth tuner,which describes the sensitivity of the jth tuner to adjust the ith amplitude. The formula above can be rearranged in matrix form as follows:

    V =M·T,

    Fig.7 (Color online) a Peak amplitude distribution of all the accelerating gaps, b frequency responses of the tuners at different longitudinal positions

    where V and T correspond to the variations of the field amplitudes and tuner settings, respectively. In addition,Mij=?Vi/?Tjis the response matrix. The dimensions of the response matrix are determined by the number of measurement positions and tuners. For this prototype, the measurement positions corresponding to the accelerating field and focusing field should both be included in the matrix.However,the desired values of the quadrupole field amplitudes are difficult to predict accurately owing to the systematic errors of the measurement, as mentioned in Sect. 4.2.1.Considering that the defined unevenness of the accelerating field is sufficient to characterize the uniformity of the overall field distribution, only the measurement positions corresponding to the accelerating field are selected,and the response matrix M is defined as a 9 × 20 non-square matrix.

    To establish the response matrix, 20 slug tuners were individually inserted 5 mm into the cavity.After each tuner was inserted,the accelerating field was measured to obtain the field variations at the measurement positions. Singular value decomposition can be used to invert the response matrix [24].

    5.2 Tuning results

    5.2.1 Field tuning

    The normalized peak amplitudes in all the accelerating gaps after the first tuning iteration are plotted in Fig.7a with those of the simulated and measured results before tuning. The field unevenness defined by the amplitudes in the nine gaps between the internal drift tubes is still 1.3%.However,the field unevenness is within 0.8% if the eighth gap is not considered.Furthermore,the peak amplitudes in the gaps at both ends of the cavity and the gaps between the drift tubes and the RFQs indicate that a more uniform field distribution along the entire cavity was obtained after the first tuning iteration. After the second tuning iteration, the peak amplitude in the eighth gap remained almost unchanged according to the measurement results. An examination of the response matrix revealed that the quantified tuning capability of each tuner on the eighth gap is one order of magnitude smaller than that of the tuners on the other gaps. After the first tuning iteration, the relative errors of the quadrupole fields of the two RFQ sections were within ± 0.7% and ± 0.8%, respectively.

    5.2.2 Frequency tuning

    The frequency of the cavity was tuned after field tuning by inserting all the tuners to the same depth to maintain the field distribution. The tuners at the same longitudinal position have the same effect on the cavity frequency;therefore, only the five tuners in the first quadrant were used to obtain the frequency response curves.As shown in Fig.7b,the cavity frequency and the insertion depth of the tuner seem to have a linear relationship when the insertion depth is greater than 6 mm. When the insertion depth is less than 6 mm, the linear relationship is invalid because the head of the tuner was machined with a fillet of 6 mm.For tuners at different longitudinal positions, the tuning coefficients of the tuners at both ends of the cavity are smaller than those of the tuners in the middle. The cavity frequency was tuned to 324.979 MHz by pulling all the tuners out by 11.02 mm. Then, the field distribution was confirmed without significant changes.

    6 Summary

    An aluminum prototype of the spatially periodic RFQ structure was designed, built, and tested to determine the accelerating and focusing field distributions, tuning performance, and mechanical structure and process technology of the cavity. Mechanically, the cavity is an overconstrained structure owing to the array of stems;thus,it requires high-precision machining and is difficult to assemble. The use of threaded connections between the drift tube components and the cavity walls mitigates the issue of overconstraint and enables manual operation of the drift tube components during the assembly process,thereby reducing the risk of assembly failure. A method of processing the drift tube components was proposed to ensure precise processing. After the cavity was assembled, the concentricity of the drift tubes measured by a laser tracker was φ0.22 mm. The RF measurement results of the prototype cavity showed good agreement with the simulation overall.The unevenness of the accelerating field was 1.3%,and the relative errors of the quadrupole fields of the two RFQ sections were within ± 0.4% and ± 1.0%, respectively.A more uniform field distribution was obtained after a tuning experiment.

    The successful fabrication and testing of the aluminum prototype provided valuable experience for the development of the spatially periodic RFQ structure. Next, the thermal performance of the cavity will be analyzed,and the structure and fabrication technology for a full-power prototype with vacuum seals and cooling channels will be studied.

    AcknowledgementsThe authors would like to thank Lanzhou Ruiyuan Machinery Equipment Co. Ltd. for collaboration on the mechanical design and manufacture of the cavity.Thanks are also due to Dr. Shichun Huang for language revision.

    Author contributionsZhi-Jun Wang and Yuan He conceived the idea of the study. The design of mechanical structure and process scheme was performed by Bin Zhang, Feng-Feng Wang and Pei-Yan Yu. Guo-Zhen Sun implemented the inspection of the cavity.The cold test was performed by Pei-Yan Yu and Chen-Xing Li. The remaining authors provided valuable suggestions and support for this study.The initial draft of the manuscript was written by Pei-Yan Yan,and all authors contributed to the writing and revisions of the manuscript.

    两个人看的免费小视频| 国产男女内射视频| 正在播放国产对白刺激| a级毛片黄视频| 999久久久精品免费观看国产| 老司机靠b影院| 啦啦啦免费观看视频1| 亚洲精品国产精品久久久不卡| 免费看十八禁软件| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 欧美人与性动交α欧美软件| 97在线人人人人妻| 美女脱内裤让男人舔精品视频| 久久性视频一级片| 操美女的视频在线观看| 午夜福利视频精品| 欧美+亚洲+日韩+国产| 中文字幕人妻丝袜一区二区| 日本av免费视频播放| 男女国产视频网站| 亚洲中文字幕日韩| 精品乱码久久久久久99久播| 久久精品亚洲av国产电影网| 久久综合国产亚洲精品| 侵犯人妻中文字幕一二三四区| 999久久久精品免费观看国产| 9191精品国产免费久久| 国产亚洲精品久久久久5区| 五月开心婷婷网| 亚洲久久久国产精品| 日韩三级视频一区二区三区| 视频区欧美日本亚洲| 天堂中文最新版在线下载| 性色av乱码一区二区三区2| a级毛片在线看网站| 青草久久国产| 99久久国产精品久久久| 亚洲免费av在线视频| 大码成人一级视频| 亚洲精品国产色婷婷电影| 免费一级毛片在线播放高清视频 | 久久久久视频综合| 久久精品亚洲熟妇少妇任你| 天天躁狠狠躁夜夜躁狠狠躁| 天天添夜夜摸| 久久女婷五月综合色啪小说| 一本久久精品| 免费av中文字幕在线| 爱豆传媒免费全集在线观看| 免费观看a级毛片全部| 操出白浆在线播放| 亚洲成国产人片在线观看| 亚洲欧洲日产国产| 欧美久久黑人一区二区| 美女福利国产在线| 欧美性长视频在线观看| 侵犯人妻中文字幕一二三四区| 999久久久国产精品视频| 国产精品久久久av美女十八| 精品卡一卡二卡四卡免费| 99九九在线精品视频| 欧美人与性动交α欧美软件| 老司机影院成人| av不卡在线播放| 精品国产乱子伦一区二区三区 | 999久久久精品免费观看国产| 国产精品影院久久| 国产xxxxx性猛交| 91麻豆精品激情在线观看国产 | 又大又爽又粗| 18禁黄网站禁片午夜丰满| 亚洲国产欧美在线一区| 丰满迷人的少妇在线观看| 欧美日韩国产mv在线观看视频| 亚洲五月色婷婷综合| 国产深夜福利视频在线观看| 国产亚洲av片在线观看秒播厂| 99久久人妻综合| 日韩制服丝袜自拍偷拍| 国产欧美日韩一区二区三区在线| 国产亚洲av片在线观看秒播厂| av欧美777| 黄色视频不卡| 久久人妻福利社区极品人妻图片| 麻豆av在线久日| 亚洲国产av影院在线观看| 国产av又大| 日韩制服骚丝袜av| 国产人伦9x9x在线观看| 人妻人人澡人人爽人人| 超碰97精品在线观看| 12—13女人毛片做爰片一| 免费日韩欧美在线观看| a级毛片黄视频| 国产成人欧美| 亚洲成人免费av在线播放| a级片在线免费高清观看视频| 人人妻人人澡人人爽人人夜夜| av超薄肉色丝袜交足视频| 日本wwww免费看| 亚洲精品国产色婷婷电影| 黑人操中国人逼视频| 亚洲中文字幕日韩| 一个人免费看片子| 在线天堂中文资源库| 久久精品成人免费网站| 欧美一级毛片孕妇| 黄片播放在线免费| 久久久久精品国产欧美久久久 | 一边摸一边抽搐一进一出视频| 日韩制服丝袜自拍偷拍| 他把我摸到了高潮在线观看 | 亚洲国产欧美在线一区| 一区在线观看完整版| 午夜久久久在线观看| 精品卡一卡二卡四卡免费| 韩国高清视频一区二区三区| 无遮挡黄片免费观看| 黄片小视频在线播放| 他把我摸到了高潮在线观看 | 一本一本久久a久久精品综合妖精| 精品国产乱码久久久久久男人| 免费看十八禁软件| 国产区一区二久久| 免费不卡黄色视频| 成人影院久久| 超碰成人久久| 黑人巨大精品欧美一区二区蜜桃| 欧美黄色片欧美黄色片| 在线看a的网站| 久久久久国内视频| a级片在线免费高清观看视频| 亚洲av欧美aⅴ国产| 亚洲综合色网址| 99国产精品一区二区三区| 日韩制服丝袜自拍偷拍| 成年人免费黄色播放视频| 纵有疾风起免费观看全集完整版| 亚洲免费av在线视频| 日韩免费高清中文字幕av| 亚洲精品一区蜜桃| 91av网站免费观看| 两人在一起打扑克的视频| 午夜福利视频精品| 国产亚洲av高清不卡| 久久精品aⅴ一区二区三区四区| 亚洲欧美一区二区三区黑人| 久久久久久亚洲精品国产蜜桃av| 99国产精品99久久久久| 韩国精品一区二区三区| 久久中文看片网| 桃红色精品国产亚洲av| 黑人猛操日本美女一级片| 欧美日韩av久久| 午夜福利一区二区在线看| 最近最新免费中文字幕在线| 色综合欧美亚洲国产小说| 两性午夜刺激爽爽歪歪视频在线观看 | av电影中文网址| 激情视频va一区二区三区| 国产淫语在线视频| 最近中文字幕2019免费版| 18禁黄网站禁片午夜丰满| 久久国产精品人妻蜜桃| 国产精品一区二区精品视频观看| 国产一区有黄有色的免费视频| 美女高潮到喷水免费观看| 国产精品偷伦视频观看了| 午夜成年电影在线免费观看| 夜夜骑夜夜射夜夜干| 涩涩av久久男人的天堂| 90打野战视频偷拍视频| 一区二区av电影网| 国产成人免费无遮挡视频| 老汉色∧v一级毛片| 一区二区日韩欧美中文字幕| 国产高清视频在线播放一区 | 黄网站色视频无遮挡免费观看| 在线观看人妻少妇| 老熟女久久久| 亚洲国产成人一精品久久久| 久久久久国产精品人妻一区二区| 一级a爱视频在线免费观看| 午夜91福利影院| 一区二区三区乱码不卡18| 人妻人人澡人人爽人人| 欧美人与性动交α欧美软件| 老司机影院成人| 亚洲精品第二区| 性色av一级| 午夜精品国产一区二区电影| 精品人妻一区二区三区麻豆| 午夜福利一区二区在线看| 欧美黄色片欧美黄色片| 黑丝袜美女国产一区| cao死你这个sao货| 亚洲激情五月婷婷啪啪| 91精品国产国语对白视频| avwww免费| 美女主播在线视频| 亚洲精品第二区| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 少妇精品久久久久久久| 亚洲人成电影免费在线| 欧美成人午夜精品| 欧美在线一区亚洲| 天天躁日日躁夜夜躁夜夜| 在线永久观看黄色视频| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 男女下面插进去视频免费观看| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 亚洲五月婷婷丁香| 成年人午夜在线观看视频| 国产欧美亚洲国产| 亚洲欧洲精品一区二区精品久久久| 亚洲专区中文字幕在线| 两性夫妻黄色片| 不卡av一区二区三区| 老汉色av国产亚洲站长工具| 男女国产视频网站| 国产真人三级小视频在线观看| 99香蕉大伊视频| 97人妻天天添夜夜摸| 午夜福利免费观看在线| 高潮久久久久久久久久久不卡| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 久久久久网色| 五月开心婷婷网| 久久这里只有精品19| 免费观看人在逋| 亚洲黑人精品在线| 可以免费在线观看a视频的电影网站| 日日爽夜夜爽网站| 欧美日韩av久久| 男人添女人高潮全过程视频| 最近最新中文字幕大全免费视频| 久久久国产一区二区| 亚洲欧美精品综合一区二区三区| 国产av一区二区精品久久| 亚洲精品久久久久久婷婷小说| 久久综合国产亚洲精品| 女人被躁到高潮嗷嗷叫费观| 欧美激情极品国产一区二区三区| 女警被强在线播放| 免费一级毛片在线播放高清视频 | 亚洲 欧美一区二区三区| 在线亚洲精品国产二区图片欧美| 成年人免费黄色播放视频| 这个男人来自地球电影免费观看| 亚洲欧洲日产国产| 一本大道久久a久久精品| 亚洲第一av免费看| 国产精品久久久久久人妻精品电影 | 老鸭窝网址在线观看| 久久久久精品国产欧美久久久 | 亚洲国产精品999| 国产色视频综合| 丁香六月欧美| 亚洲精品一区蜜桃| 精品国产乱码久久久久久男人| 午夜成年电影在线免费观看| 老司机靠b影院| 俄罗斯特黄特色一大片| 午夜成年电影在线免费观看| 精品少妇一区二区三区视频日本电影| 欧美97在线视频| 国产精品一区二区在线观看99| 国产成人av教育| 青青草视频在线视频观看| 亚洲精品国产区一区二| 亚洲av日韩精品久久久久久密| 成人国产av品久久久| 国产99久久九九免费精品| 啦啦啦视频在线资源免费观看| 水蜜桃什么品种好| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 大片免费播放器 马上看| 国产亚洲av片在线观看秒播厂| 伊人亚洲综合成人网| 91老司机精品| www.自偷自拍.com| 久久国产精品大桥未久av| 亚洲国产av新网站| 又黄又粗又硬又大视频| 午夜视频精品福利| 久久av网站| 中文欧美无线码| 亚洲精品av麻豆狂野| 99国产极品粉嫩在线观看| 满18在线观看网站| 国产片内射在线| 蜜桃国产av成人99| 色94色欧美一区二区| 成年女人毛片免费观看观看9 | 成在线人永久免费视频| 亚洲国产欧美在线一区| 久久人人爽人人片av| 亚洲精品粉嫩美女一区| 菩萨蛮人人尽说江南好唐韦庄| a级毛片黄视频| 蜜桃在线观看..| 一区二区av电影网| 99国产精品一区二区三区| 波多野结衣av一区二区av| 国产成人精品久久二区二区91| 一级片'在线观看视频| 午夜两性在线视频| 亚洲精品国产精品久久久不卡| 日本vs欧美在线观看视频| 亚洲va日本ⅴa欧美va伊人久久 | 久久天躁狠狠躁夜夜2o2o| 脱女人内裤的视频| av在线老鸭窝| 丝袜人妻中文字幕| 精品国产乱子伦一区二区三区 | 久久中文看片网| 国产欧美日韩精品亚洲av| 亚洲一码二码三码区别大吗| 亚洲精品久久久久久婷婷小说| 91精品伊人久久大香线蕉| 亚洲av美国av| 一本色道久久久久久精品综合| 狠狠精品人妻久久久久久综合| 国产免费福利视频在线观看| 亚洲专区字幕在线| 国产亚洲欧美在线一区二区| 黄片播放在线免费| 午夜福利一区二区在线看| 亚洲精品乱久久久久久| 精品亚洲乱码少妇综合久久| 亚洲国产av新网站| 黑人巨大精品欧美一区二区蜜桃| 69精品国产乱码久久久| 婷婷色av中文字幕| 麻豆国产av国片精品| 久热爱精品视频在线9| 国产精品九九99| 黑人欧美特级aaaaaa片| 国产又色又爽无遮挡免| 国产片内射在线| 午夜影院在线不卡| 狂野欧美激情性xxxx| 99国产极品粉嫩在线观看| 最黄视频免费看| 欧美成人午夜精品| 久久久久国内视频| 又紧又爽又黄一区二区| 最近最新中文字幕大全免费视频| 超碰成人久久| 男女之事视频高清在线观看| 又黄又粗又硬又大视频| 日韩欧美免费精品| 高清欧美精品videossex| 久久久国产精品麻豆| 高潮久久久久久久久久久不卡| 国产欧美日韩一区二区三 | 久久久久久久国产电影| 亚洲天堂av无毛| 亚洲熟女毛片儿| 亚洲国产欧美一区二区综合| 777久久人妻少妇嫩草av网站| 亚洲一码二码三码区别大吗| 日韩 亚洲 欧美在线| 欧美精品一区二区大全| 老熟妇乱子伦视频在线观看 | 少妇猛男粗大的猛烈进出视频| 亚洲精品美女久久av网站| 丰满迷人的少妇在线观看| 他把我摸到了高潮在线观看 | 一边摸一边做爽爽视频免费| 亚洲人成电影观看| 亚洲熟女毛片儿| 久久精品亚洲av国产电影网| 亚洲精品在线美女| 最黄视频免费看| 久久久久久人人人人人| 日本一区二区免费在线视频| 国产精品影院久久| 国产精品亚洲av一区麻豆| 精品亚洲成国产av| 少妇的丰满在线观看| 一个人免费看片子| 老司机福利观看| 巨乳人妻的诱惑在线观看| 久久久久视频综合| 日韩精品免费视频一区二区三区| 久久久久精品国产欧美久久久 | 久久久久国产一级毛片高清牌| 无遮挡黄片免费观看| 日韩大片免费观看网站| 天天影视国产精品| 国产亚洲午夜精品一区二区久久| 国产免费av片在线观看野外av| 高清在线国产一区| 在线精品无人区一区二区三| 真人做人爱边吃奶动态| 国产精品自产拍在线观看55亚洲 | 免费在线观看日本一区| 久久久国产精品麻豆| 精品国产国语对白av| 国产男女超爽视频在线观看| 无限看片的www在线观看| 国产免费福利视频在线观看| 亚洲专区国产一区二区| 亚洲五月婷婷丁香| 18禁观看日本| 性少妇av在线| 亚洲国产毛片av蜜桃av| 2018国产大陆天天弄谢| 日韩欧美一区二区三区在线观看 | 中亚洲国语对白在线视频| 国产精品一区二区精品视频观看| 桃红色精品国产亚洲av| 丝袜在线中文字幕| 日本vs欧美在线观看视频| 97在线人人人人妻| 色综合欧美亚洲国产小说| 亚洲国产看品久久| 搡老乐熟女国产| 亚洲av美国av| 男男h啪啪无遮挡| 亚洲人成77777在线视频| 丰满少妇做爰视频| 在线看a的网站| 777久久人妻少妇嫩草av网站| 汤姆久久久久久久影院中文字幕| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 久久久久久久大尺度免费视频| 国产一区二区在线观看av| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 久久av网站| 亚洲 欧美一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 大型av网站在线播放| 欧美人与性动交α欧美软件| 精品福利永久在线观看| 亚洲欧美激情在线| 欧美人与性动交α欧美精品济南到| 亚洲全国av大片| 夜夜夜夜夜久久久久| 午夜影院在线不卡| 精品一区在线观看国产| 国产一区二区三区综合在线观看| 9191精品国产免费久久| 大型av网站在线播放| 欧美成人午夜精品| 女警被强在线播放| 巨乳人妻的诱惑在线观看| 国产片内射在线| 国产在线视频一区二区| 丰满少妇做爰视频| 国产精品 国内视频| 777久久人妻少妇嫩草av网站| 黄色视频,在线免费观看| 91大片在线观看| 精品一区二区三卡| 18禁裸乳无遮挡动漫免费视频| 国产淫语在线视频| 青青草视频在线视频观看| 久久久久国产一级毛片高清牌| 欧美变态另类bdsm刘玥| 啦啦啦中文免费视频观看日本| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 欧美亚洲 丝袜 人妻 在线| 欧美 日韩 精品 国产| 国产麻豆69| 男人添女人高潮全过程视频| 淫妇啪啪啪对白视频 | 黑人欧美特级aaaaaa片| 国产男女超爽视频在线观看| 欧美日韩精品网址| 久久精品亚洲av国产电影网| 精品人妻1区二区| 国产成人啪精品午夜网站| 99精品欧美一区二区三区四区| 一级毛片精品| 成年美女黄网站色视频大全免费| 又紧又爽又黄一区二区| 亚洲欧美成人综合另类久久久| 美女午夜性视频免费| 91麻豆av在线| 人人妻,人人澡人人爽秒播| 国产国语露脸激情在线看| 亚洲精品中文字幕在线视频| 国产av又大| 亚洲精品中文字幕在线视频| 国产精品亚洲av一区麻豆| 国产精品久久久人人做人人爽| 搡老乐熟女国产| 人人妻人人澡人人爽人人夜夜| 女性被躁到高潮视频| 亚洲人成电影观看| 精品人妻1区二区| 国产精品欧美亚洲77777| 精品少妇黑人巨大在线播放| 一区二区三区乱码不卡18| 午夜福利在线观看吧| 亚洲国产成人一精品久久久| 久久久久国产一级毛片高清牌| 国产伦理片在线播放av一区| 天天躁日日躁夜夜躁夜夜| 久久免费观看电影| 老熟妇仑乱视频hdxx| 亚洲成国产人片在线观看| 亚洲视频免费观看视频| 99久久综合免费| 交换朋友夫妻互换小说| 在线观看www视频免费| 一区二区av电影网| 日韩视频在线欧美| 黄色a级毛片大全视频| 丝袜美腿诱惑在线| 天堂俺去俺来也www色官网| 国产一区二区在线观看av| 午夜福利乱码中文字幕| 色视频在线一区二区三区| 一个人免费在线观看的高清视频 | 免费在线观看完整版高清| 亚洲国产精品一区三区| 亚洲av日韩精品久久久久久密| 欧美激情久久久久久爽电影 | 人人妻人人澡人人看| 精品一区二区三卡| 亚洲精品av麻豆狂野| 亚洲欧美一区二区三区黑人| 欧美日韩黄片免| 一本色道久久久久久精品综合| 免费在线观看完整版高清| 日韩免费高清中文字幕av| 亚洲av国产av综合av卡| 黑人巨大精品欧美一区二区蜜桃| av一本久久久久| 人人妻人人澡人人爽人人夜夜| 久热爱精品视频在线9| 亚洲国产av影院在线观看| 久久毛片免费看一区二区三区| 亚洲av成人一区二区三| 这个男人来自地球电影免费观看| 99久久人妻综合| 一级毛片电影观看| 亚洲性夜色夜夜综合| 亚洲精品第二区| 久久久精品国产亚洲av高清涩受| 久久天堂一区二区三区四区| 王馨瑶露胸无遮挡在线观看| 色精品久久人妻99蜜桃| 久久精品国产a三级三级三级| svipshipincom国产片| 女人久久www免费人成看片| 丝袜脚勾引网站| 一区二区av电影网| 亚洲成人国产一区在线观看| 91九色精品人成在线观看| a在线观看视频网站| 国产野战对白在线观看| 久久热在线av| 51午夜福利影视在线观看| 亚洲精品中文字幕在线视频| 亚洲avbb在线观看| 国产精品一区二区在线观看99| 一级a爱视频在线免费观看| 日本vs欧美在线观看视频| 免费在线观看完整版高清| 不卡av一区二区三区| 欧美人与性动交α欧美精品济南到| 男女无遮挡免费网站观看| 国产1区2区3区精品| 电影成人av| 新久久久久国产一级毛片| 亚洲av电影在线观看一区二区三区| 久久人妻熟女aⅴ| 欧美亚洲日本最大视频资源| 午夜久久久在线观看| 成人亚洲精品一区在线观看| 国产亚洲av高清不卡| 首页视频小说图片口味搜索| 中文字幕人妻丝袜一区二区| 欧美精品啪啪一区二区三区 | 天堂俺去俺来也www色官网| 91麻豆av在线| 91成年电影在线观看| 亚洲中文字幕日韩| av网站在线播放免费| 女人久久www免费人成看片| 99热全是精品| 精品欧美一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区 | 一本久久精品| 午夜精品久久久久久毛片777| av欧美777| 亚洲国产欧美网| 欧美成狂野欧美在线观看| 久久狼人影院| 国产免费一区二区三区四区乱码| 国产黄频视频在线观看| 日韩电影二区| 精品国产国语对白av| 欧美日韩视频精品一区| 久久狼人影院| 欧美乱码精品一区二区三区| 成年人黄色毛片网站| 久久精品aⅴ一区二区三区四区| 王馨瑶露胸无遮挡在线观看| 亚洲色图 男人天堂 中文字幕| 国产av一区二区精品久久| 国产在视频线精品|