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

    Influence of stealth aircraft dynamic RCS peak on radar detection probability

    2023-04-22 02:04:52XiaoqiangLUJunHUANGYacongWULeiSONG
    CHINESE JOURNAL OF AERONAUTICS 2023年3期

    Xiaoqiang LU, Jun HUANG, Yacong WU, Lei SONG

    School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China

    KEYWORDS Maneuver;Penetration;RCS;RCS fluctuation statistical model;RCS peak exposure time;RCS peak width;Stealth;Trajectory planning

    Abstract For modern stealth aircraft,it is important to analyze the influence of Radar Cross Section(RCS)peak exposure on penetration for guiding stealth design and penetration trajectory planning,which needs to reflect the RCS statistical uncertainty and the RCS difference with the change of incident angle.Based on the RCS characteristics of typical stealth aircraft,this paper established a simplified RCS dynamic fluctuation statistical model with the parameters log mean and log standard deviation.According to the detection probability algorithm in radar signal processing field,this paper built the algorithm of radar detection probability based on the RCS dynamic fluctuation statistical model.The analysis of examples concluded that the key to successful penetration is to shorten the RCS peak exposure time, which can be reduced by decreasing the RCS peak width or increasing velocity.Based on the conclusion,this paper proposed the method of turning maneuvering to reduce RCS peak exposure time dramatically.

    1.Introduction

    In the stealth design of modern aircraft,radar echoes are often concentrated in several small non-head angular regions to ensure the stealth performance during frontal penetration.However,in the real penetration process,it is difficult to maintain the frontal penetration, which will inevitably make the peak angle area toward the radar.Besides, the developing metasurface technology could control the angular distribution of electromagnetic wave energy,1,2and might change the angular position of Radar Cross Section (RCS) peak.Therefore, it is necessary to analyze the influence of RCS peak exposure on penetration, in order to guide the stealth design and penetration trajectory planning.

    To analyze the penetration, radar detection probability should be used as an index to evaluate trajectory.3As a result,the problem is transformed into the research on the influence of stealth aircraft dynamic RCS peak on radar detection probability.

    Since the RCS of stealth aircraft changes drastically with the azimuth angle,4–7and the incidence of radar waves in dynamic flight is a random event, the process of the analysis needs to reflect the RCS statistical uncertainty and the variation of radar incident angle with time during the dynamic flight.

    Firstly, to reflect the RCS statistical uncertainty, researchers usually divide the aircraft into multiple angular regions in the circumferential direction, and use fluctuation statistical model8–13with different parameters in each angular region to carry out the analysis.Chen14proposed a 30-degree-angularregion fluctuation statistical model, and developed a set of detection probability algorithm suitable for various specific aircraft targets.Li et al.15proposed target circumferential RCS scattering model.They used this model to calculate the probability of radar discovery.This kind of method is mainly used to describe the stealth performance of the aircraft itself,and to measure the penetration capability of the aircraft by calculating the radar discovery probability in a frontal penetration scenario(facing the radar,fast approaching).Although the model parameters are different in each angular domain,they do not change continuously with the azimuth angle.The focus is on the description of the aircraft’s own capabilities.That is, it pays attention to radar stealth performance rather than trajectory planning.However,considering the trajectory planning factors,the radar incident angle of the aircraft is constantly rapidly changing under a specific penetration threat and flight trajectory.Therefore, when conducting the simulation of the radar detection process in a given scenario,it is necessary to use a statistical model of RCS fluctuations which can reflect the difference in radar incident angle.

    Secondly, to reflect the variation of radar incident angle with time during the dynamic flight,researchers establish algorithm and calculate the radar detection probability by different typical maneuver penetration methods.Yuan et al.16used the RCS deterministic model(the method in which the RCS is constant on a specific radar incident azimuth)to express the target characteristics of the aircraft, and analyzed the beneficial effects of the serpentine maneuver on the penetration.However,for the RCS deterministic model,the calculation interval of radar incident angle is in a dilemma.If the interval is not small enough,the drastic RCS changes caused by small differences in azimuth angles may cause large errors;if the interval is small enough, the calculation time will be unacceptably long.Zhang et al.17and others superimposed jitter information on the RCS deterministic model to characterize the statistical characteristics of the echo, but its radar detection probability is calculated based on Swerling I model, which is different from the typical echo statistical characteristics of stealth aircraft.The methods above reflect the time-varying radar incident angle, mainly to obtain the calculation results of typical scenarios, and select a better penetration trajectory for stealth aircraft.However,the real trajectory planning process is highly nonlinear, and the optimization efficiency can be accelerated by adding appropriate constraints.Therefore, it is necessary to explore the influence of the change of parameters on the penetration effect, analyze the causes of the influence, and use it to constrain the trajectory planning,which can also provide some guidance for the stealth design.

    It can be seen that we need to combine the two perspectives above.Firstly, an RCS fluctuation statistical model should be established under the time-varying radar incident azimuth to analyze the change of radar echo during the dynamic flight,and the process of radar detecting aircraft should be described accurately.Secondly, a variety of penetration examples in given scenarios need to be calculated.Based on the results, it is necessary to analyze the factors that affect the radar detection probability in penetration, and find the design methods and maneuver measures that are beneficial to penetration.

    In this paper,based on the RCS data of a typical stealth aircraft, we introduce a simplified RCS dynamic fluctuation statistical model.For this model, we define statistical model parameters to characterize the dynamic stealth performance of the aircraft, including log means, log standard deviation,and peak width.Referencing the detection probability algorithm in radar signal processing, we construct an algorithm of the radar detection probability based on the simplified RCS dynamic fluctuation statistical model.Using this method,we explore the influence of changes in aircraft dynamic stealth performance parameters on the radar detection probability in different penetration trajectories, and obtain the influence of aircraft RCS peak exposure and maneuver on the penetration effect.Meantime, we make suggestions for maneuver method that can be used for trajectory planning.

    2.Method

    2.1.Simplified RCS dynamic fluctuation statistical model and basic assumptions

    2.1.1.Simplified RCS dynamic fluctuation statistical model

    According to the analysis of the angular domain RCS fluctuation model,18,19we choose the lognormal model to build a simplified model of stealth aircraft dynamic target characteristics.

    The RCS lognormal distribution model could be expressed as

    where Λ represents the radar incident azimuth, and σRCSrepresents the RCS value.The parameter μ(Λ) represents the azimuth-varying log mean of the RCS, and the parameter s(Λ) represents the azimuth-varying log standard deviation of the RCS.

    According to the RCS data of typical stealth aircraft,20in the non-RCS-peak angular region, the RCS fluctuates with a low level, corresponding to a small parameter μ; in the RCSpeak-center angular region, the RCS fluctuates with a high level, corresponding to a large parameter μ.At the edge of the RCS-peak angular region (the high value-low value junction),the RCS value fluctuates significantly,so the corresponding log standard deviation parameter s is relatively large.

    In order to make the simplified RCS dynamic fluctuation statistical model be factual, we referenced the RCS data of a flying wing, and the schematic diagram is shown in Fig.1.

    Based on the information above,a simplified RCS dynamic fluctuation statistical model of stealth aircraft is established.The peak value of the RCS model is assumed to be 10 m2.4 peaks symmetrically distribute in the circumferential direction of RCS, and the azimuth of the peak center is ±35° and±145°.The valley value of the RCS model is assumed to be 0.01 m2.

    The parameters of the simplified RCS dynamic fluctuation statistical model are shown in Fig.2.

    Fig.1 Typical RCS data of a flying wing.

    Fig.2 Parameters of fluctuation model.

    As shown in Fig.2, the peak value μpeakand valley value μvalleyof the azimuth-varying log mean are

    According to the previous analysis,the azimuth-varying log standard deviation of the model is defined as:svalleyis 0.5 at the valley, sfootis 1.01 at the junction of peak and valley, speakis 0.01 at the peak, the largest log standard deviation is at the junction of peak and valley,and the smallest one is at the center of peak.

    The definition of peak width is shown in Fig.2.The peak angular region of the model is 34.4°–35.6°, and the RCS Peak Width (RPW) is 1.2°.

    2.1.2.Basic assumptions

    To analyze the factors influencing penetration, this paper is based on the following basic assumptions:

    (1) The simplified RCS dynamic fluctuation statistical model is used to reflect the RCS fluctuation of aircraft.

    (2) According to the typical radar detection level, the radar model is defined as being able to detect 1 m2target at 240 km(i.e.20 dB signal-to-noise ratio).The pulse accumulation number is 2;the correlation coefficient is 1;the scanning period is 12 s; the false alarm probability pfis 1×10-6.

    (3) The calculation of radar detection probability is based on the single-detection-found assumption, that is, the aircraft is found at any point in the trajectory, which means the failure of the penetration trajectory.

    (4) The influence of roll angle on target characteristics is not considered when the aircraft maneuvers.

    (5) The turning radius is restricted by 2G load factor when the aircraft maneuvers.

    2.2.Algorithm of radar detection probability

    The establishment process of the method is as follows:

    (1) Based on the mathematical and physical model of radar detection,21,22the probability of radar’s detection of a specific attitude aircraft at a specific spatial position is obtained.

    (2) In combination with the real flight situation, the radar incident angle is calculated during the aircraft motion simulation process.During this process,the RCS fluctuation characteristics in the angle domain of the aircraft are converted into the fluctuation characteristics in the time domain, which can be substituted into the detection.The probability of the aircraft being detected by a radar’s complete scan can be obtained with this model.

    (3) The detection probabilities of multiple scans are superimposed to obtain the probability of detection of the aircraft during the entire process of flying over a specific trajectory.

    2.2.1.Algorithm of radar detection probability in single position of aircraft

    The most widely used in radar signal detection is the threshold detection based on the Neyman–Pearson criterion.Once a radar receives echoes, the receiving and detecting system first performs matched filtering on a single pulse signal in the intermediate frequency strip, and then N pulses are weighted accumulating.The cumulative output is compared with a specific threshold voltage.If the output envelope exceeds the threshold, the target is considered to exist, and vice versa.In this case, the detection probability of the target can be shown by

    where pdis the probability of detection,r is the integrator output normalized to the average per pulse noise power, pr(r ) is the probability density function of r, and rbis the detection threshold, which can be obtained from the false alarm probability pf.The false alarm probability refers to the probability that the noise exceeds the threshold level and is mistaken for the signal.Lowering the threshold will increase the false alarm probability,while increasing the threshold will lose some weak signals.The Neyman–Pearson criterion, which maximizes the probability of detection under the constraint of false alarm probability, can help operators or electronic systems to select the appropriate threshold level rb, that is,

    In order to obtain the probability density function pr(r ),according to Kanter’s correlation function model,21do the Laplace transform to pr(r )

    In summary, the radar detection probability can be calculated with the RCS fluctuation.

    2.2.2.Algorithm of single-scan radar detection probability in time domain

    In the real penetration, during the process of radar scanning,the distance and azimuth of the aircraft relative to the radar will change greatly.Therefore, according to the specific flight trajectory and radar scanning time, the angle-domain RCS fluctuation characteristics should be transformed into the time-domain RCS fluctuation characteristics, that is, putting the change of azimuth with time into the probability density function in angle domain, and obtaining the probability density function in time domain.

    where μ(Λ (t ))and s(Λ (t ))represent the parameters μ and s of fluctuation statistical model corresponding to a specific azimuth Λ(t ).

    In the real calculation, because the time of scanning to the aircraft is uncertain in a single scan, the probability density function of a single scan is obtained by integrating the probability density function in time domain.

    where tsand teare the start and end time of a single scan respectively.

    Significantly, due to the numerical algorithm, in order to prevent the RCS fluctuation characteristics from changing too much in each infinitesimal segment and affecting the calculation results, it is necessary to adjust the start and end time interval of each infinitesimal segment in the integration, in order to make the azimuth change within the segment less than 0.1° and ensure the accuracy of numerical calculation.

    Substituting time-domain RCS fluctuation characteristics into algorithm of radar detection probability in single position of aircraft (presented in Section 2.2.1), we have

    Since there are multiple scans in the whole trajectory, the detection probability of the whole trajectory (pd_route) can be obtained by accumulating the detection probability of the multi-segment single scans in the trajectory,

    3.Penetration evaluation under different influencing factors

    3.1.Influence of RCS peaks and valleys on the probability of radar detection

    In this section,the top-view schematic diagram of the penetration scene is represented by the Cartesian coordinate system XoY, in kilometer.

    For the mission,

    (1) Defining a typical mission, the aircraft needs to penetrate from (0, -100) to (0,100).

    (2) Defining the Straight Trajectory:a straight-line segment from (0, -100) to (0,100).

    For the threat environment:

    (1) Defining the Radar Distance(R): the distance between radar position to (0,0).

    (2) Defining the Valley Situation: the radar is located at(R,0), which means that the peak value is not exposed to the radar during the whole trajectory.

    (3) Defining the Peak Situation:the radar is located at(Rsin(35°), Rcos(35°)),that is, when the aircraft flies over the point (0,0), its maximum peak angle (35°) is exposed to the radar.

    Fig.3 shows the schematic diagram of two situations when the R equals 500 km.

    In this part,the RPW of the aircraft is 1.6,and the velocity along the trajectory is 800 km/h.

    By calculating the radar detection probability in the two threat environments, we can know the impact of RCS peak exposure on penetration.

    As shown in Fig.4, we assume that the Radar Detection Range (RDR) is radar distance (R) when the radar detection probability is 90%.

    In the Valley Situation and the Peak Situation,the RDR is about 213 km and 586 km respectively,and the detection probability when R equals 300 km is about 4.497% and 99.87%respectively.Thus when R is between 300 km and 600 km,the influence of peak exposure on radar detection is decisive.It can be seen that the peak exposure will sharply increase the radar detection probability and the RDR.For the 200-km-long trajectory, the peak exposure of this trajectory leads to the increase of the RDR by nearly 400 km.

    3.2.Influence of RCS peak width on radar detection probability

    The typical characteristic of aircraft RCS peak is its peak width, that is, the RPW.To study the influence of the RPW on radar detection probability, the dynamic stealth performance of aircraft can be evaluated by using the RPW.

    As shown in Fig.5,the trajectory is a straight-line segment from (0, -100) to (0,100).The velocity along the track is 800 km/h.The threat environment is Peak Situation.

    In the calculation, the RCS Peak Width (RPW) of the aircraft is taken as 0.4°, 0.8°, 1.2° and 1.6°, and the radar detection probability of targets with different RPW on the same trajectory is calculated respectively.

    In Fig.6,it can be seen that when R is greater than 250 km,the influence of RPW on detection probability begins to show.When R is greater than 300 km,the effect of RPW on detection probability is pronounced.In the penetration process of a specific situation (e.g.R = 500 km), 10 radar scans when the aircraft is nearby the location(0,0)are drawn out for analysis.Through the detection probability of the 10 radar scans, how RPW affects the detection probability can be analyzed.

    Fig.3 Schematic diagram of two situations.

    Fig.4 Detection probability comparison between Peak Situation and Valley Situation.

    Fig.5 Schematic diagram of penetration scene.

    Fig.6 Probability of detection with different RPW.

    As shown in Fig.7, the curve of ordered-scan detection probability shows a peak shape.When the RPW is 0.4°, there are 2 scans which can detect the aircraft probabilistically, and the total radar detection probability is 46.76%;when the RPW is 1.6°, there are 6 scans which can detect the aircraft probabilistically,and the total radar detection probability is 96.42%.

    Fig.7 Probability of detection by ordered scan with different RPW (R = 500 km).

    At this distance, the influence of peak exposure on radar detection is decisive (shown in Section 3.1).Combined with Fig.7,the number of scans which can detect the aircraft probabilistically is related to the RCS Peak Exposure Time(RPET).The RPW is positively correlated with radar detection probability, and its essence is that RPW changes the RPET.

    Therefore, on the one hand, from the perspective of evaluation, when other parameters are the same, the smaller the RPW is, the better the dynamic stealth performance is; on the other hand, from the perspective of design, stealth design can be made with minimum RPW as the goal.However, in the trajectory planning of the specified aircraft, the RPW of the aircraft cannot be changed.What can be changed at this time is the RCS peak exposure time in the trajectory.

    3.3.Influence of RCS peak exposure time on radar detection probability

    The trajectory for calculation is also straight trajectory, same as Fig.5 in Section 3.2.The threat environment is Peak Situation.The RPW of aircraft target is 0.4°.The radar distance R equals 500 km.10 radar scans when the aircraft is nearby the location (0,0) are drawn out for analysis.Through the detection probability of 10 radar scans, how RPET affects the detection probability can be analyzed.

    As shown in Fig.8, v represents the velocity of aircraft.With the increase of v, the RPET decreases, and the probability of radar detection also decreases from 96.58%at 200 km/h to 46.76% at 800 km/h.This shows that reducing the RPET can directly reduce the probability of radar detection.In fact,reducing RPW is a way to reduce RPET for stealth design.In the case of a certain peak height,the detection probability at a single point is also certain.The shorter the time is, the smaller the total detection probability is.Therefore, when aircraft has to expose its RCS peak to the radar, the target of trajectory planning is to reduce RPET.

    3.4.Influence of turning maneuver on radar detection probability

    3.4.1.Definition and principle of turning maneuver

    In a real trajectory, the velocity of aircraft has an upper limit.Therefore, if we want to reduce the RPET and the probability of detection,we need to make the RCS peaks quickly‘‘sweep”the radar.

    Fig.8 Probability of detection by ordered scan at different velocities (R = 500 km).

    Therefore, we select a typical case (the RPW is 0.4°, the velocity is 800 km/h, Peak Situation) for maneuver planning:

    When the radar incident angle is 34°,the aircraft turns right till facing the radar;when reaching the y=0 line,the aircraft makes a turning maneuver,so that the trajectory of the second half is the mirror image of the first half.

    Fig.9 shows the schematic diagram of 2 trajectories when R is 300 km.

    The maneuver changes the Straight Trajectory when radar incident angle is 1°less than 35°(34°).Thus,we define that the maneuver angle field is 35±1°,the Turning Angle(TA)is 2°,and the trajectory is ± 1° Avoidance Maneuver Trajectory.

    Under different R, the radar detection probabilities of Straight Trajectory and Maneuver Trajectory are calculated respectively, and the results are shown in Fig.10.

    When the R is more than 250 km,the effect of maneuvering turn begins to appear.When the R is more than 300 km, the effect of maneuvering turn is obvious.When the R is 300 km, the radar detection probability of Maneuver Trajectory is only 24.84%, which is far less than 73.84% of Straight Trajectory.The reason is that when the aircraft maneuvers,the RCS peak value quickly passes the radar, which is equivalent to greatly increasing the velocity of the peak direction relative to the radar, reducing the RPET, thus reducing the radar detection probability.

    Fig.9 Schematic diagram of Straight Trajectory and Maneuver Trajectory.

    Fig.10 Probability of detection in different trajectory.

    In order to show the essence of turning maneuver, Fig.11 and Fig.12 show the change of radar incident azimuth with time during turning maneuver.

    As shown in Fig.11, it can be seen that both curves cross the peak centerline (azimuth = 35°), but when crossing the peak centerline,the slope of curve in Maneuver Trajectory Situation is much greater than that in Straight Trajectory Situation.

    The specific RPET can be seen from Fig.12.For the Straight Trajectory Situation, the RCS peak exposure is from 441.74 s to 458.18 s,and the RPET is 16.44 s.For the Maneuvering Trajectory Situation, the peak RCS exposure is from 457.21 s to 457.30 s, and the RPET is 0.09 s, which is two orders of magnitude lower than that of Straight Trajectory Situation.Therefore,the turning maneuver can greatly reduce the RPET, thus reducing radar detection probability.

    3.4.2.Influence of velocity on turning maneuver effect

    Under different R,the radar detection probabilities of different trajectory and velocity are calculated respectively, and the results are shown in Fig.13.

    As shown in Fig.13,turning maneuver has a good effect on both high-velocity and low-velocity targets.For the same mission,when R is greater than about 300 km,the radar detection probability of low-velocity Maneuvering Trajectory is less than that of high-velocity Straight Trajectory, which indicates that the effect of turning maneuver is more obvious on the middle- and long-distance trajectory.It is consistent with the conclusion in Section 3.4.1.

    Fig.11 Change of azimuth with time during turning maneuver(R = 500 km, overall).

    Fig.12 Change of azimuth with time during turning maneuver(R = 500 km, partially magnified).

    Fig.13 Probability of detection with different trajectory and velocity.

    3.4.3.Influence of RPW and TA on turning maneuver effect

    Under different R,the radar detection probabilities of different trajectory and RPW are calculated respectively,and the results are shown in Fig.14.

    Fig.14 Probability of detection with different trajectory and RPW.

    Fig.15 Probability of detection on different maneuver trajectory.

    Fig.16 Probability of detection with different RPW and TA.

    As shown in Fig.14, for a fixed maneuver mode (±1°Turning Maneuver Trajectory), the change of the RPW of the aircraft has little effect on the radar detection probability after maneuver.Therefore, we can infer that the maneuver planning only needs to make the maneuver angle domain cover the dangerous angle domain.

    In order to verify this conjecture,it is necessary to calculate the radar detection probabilities of different RPW and TA.The schematic diagram of different maneuver is shown in Fig.15, and the results are shown in Fig.16.

    As shown in Fig.16, the turning maneuver could have a very good effect, as long as the TA is greater than the RPW.Compared to Fig.14, when the TA is less than the RPW,the turning maneuver works, but could not reach the effect above.

    4.Conclusions

    In this paper, a simplified RCS dynamic fluctuation statistical model is constructed for stealth aircraft,and a radar detection probability algorithm based on the model is proposed.The dynamic stealth performance of the model is characterized by the azimuth-varying log mean, the azimuth-varying log standard deviation, and the peak width.By calculating different examples, the influence on the radar detection probability with different parameters is analyzed, and the following conclusions are obtained:

    (1) The key to successful penetration is to keep the RCS peak unexposed or the exposure time short.The RCS Peak Exposure Time (RPET) is positively correlated with the radar detection probability.By decreasing the RCS Peak Width (RPW) or increasing velocity, the RPET can be reduced, so as to reduce the radar detection probability to a certain extent.

    (2) The RPET can be greatly reduced by turning maneuver defined in this paper,so as to reduce the radar detection probability.Compared with increasing velocity, the results show that the turning maneuver has more obvious effects on the middle- and long-distance trajectory.Meanwhile, the turning maneuver could have a very good effect as long as the TA is greater than the RPW.Therefore, from the aspect of stealth aircraft design, there is no need to excessively pursue too small RPW and too large velocity.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    亚洲精品国产区一区二| 伊人久久大香线蕉亚洲五| 国产无遮挡羞羞视频在线观看| 日韩制服丝袜自拍偷拍| 久久国产精品男人的天堂亚洲| 国产一卡二卡三卡精品| 成年人免费黄色播放视频| 国产免费现黄频在线看| 韩国精品一区二区三区| h视频一区二区三区| 亚洲,欧美,日韩| 人体艺术视频欧美日本| 国产熟女欧美一区二区| 日本午夜av视频| 精品久久蜜臀av无| 美女福利国产在线| 建设人人有责人人尽责人人享有的| 伊人久久大香线蕉亚洲五| 日韩电影二区| 国产一区二区三区综合在线观看| 天堂中文最新版在线下载| 视频在线观看一区二区三区| 婷婷色综合大香蕉| 一级毛片电影观看| 国产男女内射视频| 夫妻午夜视频| 精品卡一卡二卡四卡免费| 亚洲精品第二区| 午夜精品国产一区二区电影| 国产欧美日韩一区二区三 | 亚洲第一青青草原| 丝袜美腿诱惑在线| 亚洲精品乱久久久久久| 日日摸夜夜添夜夜爱| 99久久人妻综合| 在线精品无人区一区二区三| 黄频高清免费视频| 亚洲国产av影院在线观看| 午夜激情久久久久久久| 天天躁日日躁夜夜躁夜夜| 欧美日韩精品网址| 男人舔女人的私密视频| 免费在线观看完整版高清| 色网站视频免费| 超色免费av| 亚洲色图 男人天堂 中文字幕| 日韩中文字幕视频在线看片| 亚洲中文日韩欧美视频| 国产精品久久久久久精品电影小说| 免费在线观看影片大全网站 | 香蕉丝袜av| 国产黄色免费在线视频| 91麻豆av在线| av国产久精品久网站免费入址| 日韩一区二区三区影片| avwww免费| 91精品三级在线观看| 大型av网站在线播放| 别揉我奶头~嗯~啊~动态视频 | av在线播放精品| 免费在线观看视频国产中文字幕亚洲 | 久久中文字幕一级| 国产精品一国产av| 免费久久久久久久精品成人欧美视频| 无限看片的www在线观看| 欧美日韩黄片免| 日本五十路高清| 少妇被粗大的猛进出69影院| 纵有疾风起免费观看全集完整版| 久久久久久亚洲精品国产蜜桃av| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 蜜桃国产av成人99| 午夜视频精品福利| 久久精品国产a三级三级三级| 成年av动漫网址| 欧美 日韩 精品 国产| 久久精品久久久久久久性| 老司机深夜福利视频在线观看 | 国产精品国产av在线观看| 老熟女久久久| 女人高潮潮喷娇喘18禁视频| 亚洲成人免费av在线播放| 午夜福利视频在线观看免费| 丝袜人妻中文字幕| 亚洲欧美精品自产自拍| 国产亚洲欧美精品永久| 免费看av在线观看网站| 欧美激情极品国产一区二区三区| 纵有疾风起免费观看全集完整版| 国产成人91sexporn| 久久久久久久久免费视频了| 亚洲一区中文字幕在线| 精品人妻1区二区| 欧美97在线视频| 国产在线免费精品| 这个男人来自地球电影免费观看| 精品少妇久久久久久888优播| 亚洲精品国产av蜜桃| videosex国产| 777米奇影视久久| 亚洲一区中文字幕在线| 一二三四社区在线视频社区8| 最黄视频免费看| 亚洲中文日韩欧美视频| 丝袜喷水一区| 男女免费视频国产| 精品第一国产精品| 黄色 视频免费看| 五月开心婷婷网| 国产麻豆69| 午夜福利影视在线免费观看| 两个人免费观看高清视频| 最新的欧美精品一区二区| 日本av手机在线免费观看| 久久久久久久大尺度免费视频| videosex国产| 国产男女内射视频| 韩国高清视频一区二区三区| 国产免费现黄频在线看| 熟女少妇亚洲综合色aaa.| 女性被躁到高潮视频| 亚洲精品久久成人aⅴ小说| 亚洲av国产av综合av卡| 亚洲成人手机| 乱人伦中国视频| 99精国产麻豆久久婷婷| 狂野欧美激情性xxxx| 热99国产精品久久久久久7| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| 日日爽夜夜爽网站| 午夜福利乱码中文字幕| 久久精品久久久久久噜噜老黄| 男男h啪啪无遮挡| 亚洲欧美色中文字幕在线| 交换朋友夫妻互换小说| 久久99热这里只频精品6学生| 亚洲精品国产av蜜桃| 国产熟女欧美一区二区| 91老司机精品| 精品亚洲成a人片在线观看| 久久精品国产亚洲av高清一级| 激情五月婷婷亚洲| www.自偷自拍.com| 男女之事视频高清在线观看 | a 毛片基地| a级片在线免费高清观看视频| 国产主播在线观看一区二区 | 好男人电影高清在线观看| 大陆偷拍与自拍| 亚洲精品国产一区二区精华液| 亚洲欧美激情在线| 国产成人精品在线电影| 欧美日韩亚洲国产一区二区在线观看 | 欧美av亚洲av综合av国产av| 91国产中文字幕| 啦啦啦视频在线资源免费观看| av在线app专区| cao死你这个sao货| 人人澡人人妻人| 伊人久久大香线蕉亚洲五| 大话2 男鬼变身卡| 久久青草综合色| 女人高潮潮喷娇喘18禁视频| 成人手机av| 中文字幕高清在线视频| 精品视频人人做人人爽| 欧美日韩福利视频一区二区| 亚洲欧美中文字幕日韩二区| 国产一区二区激情短视频 | 在现免费观看毛片| 国产一区二区激情短视频 | 亚洲一区中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 午夜福利乱码中文字幕| 18禁黄网站禁片午夜丰满| 咕卡用的链子| 久久av网站| 日本一区二区免费在线视频| 欧美国产精品va在线观看不卡| 午夜福利一区二区在线看| 国产精品偷伦视频观看了| 中文字幕高清在线视频| 欧美激情 高清一区二区三区| 视频区图区小说| 91九色精品人成在线观看| 久久国产精品影院| 午夜免费男女啪啪视频观看| 亚洲精品美女久久久久99蜜臀 | 婷婷色综合www| 国产免费又黄又爽又色| 精品福利永久在线观看| 国产日韩一区二区三区精品不卡| 久久久久久久大尺度免费视频| 免费高清在线观看日韩| 国产成人一区二区三区免费视频网站 | 久久精品久久精品一区二区三区| 欧美日韩亚洲高清精品| 色94色欧美一区二区| 欧美黄色淫秽网站| 五月天丁香电影| 青春草视频在线免费观看| 亚洲欧美精品自产自拍| 男女边吃奶边做爰视频| 五月天丁香电影| 国产精品香港三级国产av潘金莲 | 美国免费a级毛片| 国产又爽黄色视频| 国产成人精品久久二区二区91| 欧美另类一区| 又紧又爽又黄一区二区| 赤兔流量卡办理| 成人黄色视频免费在线看| 精品高清国产在线一区| 欧美精品高潮呻吟av久久| 少妇裸体淫交视频免费看高清 | 亚洲成国产人片在线观看| 国产免费视频播放在线视频| 亚洲图色成人| 人妻一区二区av| 久久精品久久久久久噜噜老黄| 老司机深夜福利视频在线观看 | 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 精品高清国产在线一区| 搡老岳熟女国产| 精品免费久久久久久久清纯 | 啦啦啦中文免费视频观看日本| 少妇 在线观看| 久久久国产欧美日韩av| 色网站视频免费| 亚洲精品乱久久久久久| 中文字幕高清在线视频| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 久久天躁狠狠躁夜夜2o2o | 韩国精品一区二区三区| 女人久久www免费人成看片| 18禁黄网站禁片午夜丰满| 久久精品久久精品一区二区三区| 妹子高潮喷水视频| 久久av网站| 久久久久网色| 多毛熟女@视频| 老司机靠b影院| 久久久国产精品麻豆| 国产精品欧美亚洲77777| 久久人人97超碰香蕉20202| 欧美精品人与动牲交sv欧美| 久久精品亚洲熟妇少妇任你| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| av天堂在线播放| 婷婷色麻豆天堂久久| 无限看片的www在线观看| 成年人免费黄色播放视频| svipshipincom国产片| 极品少妇高潮喷水抽搐| 精品熟女少妇八av免费久了| 亚洲欧美日韩另类电影网站| 欧美精品人与动牲交sv欧美| 精品国产国语对白av| 熟女av电影| 亚洲国产欧美在线一区| www.自偷自拍.com| 在线观看免费午夜福利视频| 欧美亚洲日本最大视频资源| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 老汉色av国产亚洲站长工具| 国产人伦9x9x在线观看| 久久亚洲精品不卡| 91麻豆av在线| 国产成人系列免费观看| 国产97色在线日韩免费| 成人国产一区最新在线观看 | 亚洲成人免费电影在线观看 | 国产精品久久久久久精品古装| 久久国产精品大桥未久av| 亚洲欧美日韩高清在线视频 | 亚洲欧美色中文字幕在线| 又大又黄又爽视频免费| 亚洲精品成人av观看孕妇| 你懂的网址亚洲精品在线观看| xxxhd国产人妻xxx| 欧美黄色片欧美黄色片| 激情五月婷婷亚洲| 久久久精品免费免费高清| 亚洲成国产人片在线观看| 国产爽快片一区二区三区| 最新的欧美精品一区二区| 五月天丁香电影| 成在线人永久免费视频| 丁香六月欧美| 欧美人与善性xxx| √禁漫天堂资源中文www| 亚洲少妇的诱惑av| 亚洲国产av新网站| 午夜视频精品福利| 最黄视频免费看| 啦啦啦啦在线视频资源| 两性夫妻黄色片| 久久久久国产精品人妻一区二区| 精品福利永久在线观看| 丝袜脚勾引网站| 叶爱在线成人免费视频播放| 9191精品国产免费久久| 在线天堂中文资源库| 免费久久久久久久精品成人欧美视频| 美女中出高潮动态图| 极品人妻少妇av视频| 国产高清视频在线播放一区 | www.精华液| 精品人妻1区二区| bbb黄色大片| 免费在线观看影片大全网站 | 热99国产精品久久久久久7| 亚洲av电影在线进入| 久久久久久久精品精品| 国产一卡二卡三卡精品| 国产日韩一区二区三区精品不卡| 黄片小视频在线播放| 国产男女内射视频| 久久精品aⅴ一区二区三区四区| 国产精品人妻久久久影院| 新久久久久国产一级毛片| 女警被强在线播放| 日韩制服丝袜自拍偷拍| 精品一区二区三区四区五区乱码 | 亚洲精品自拍成人| 97人妻天天添夜夜摸| 久久99一区二区三区| 高清欧美精品videossex| 色94色欧美一区二区| 亚洲精品av麻豆狂野| 十八禁网站网址无遮挡| 看免费av毛片| 美女高潮到喷水免费观看| 欧美精品av麻豆av| 久久天堂一区二区三区四区| 精品一区在线观看国产| 久久亚洲国产成人精品v| 国产精品欧美亚洲77777| 国产淫语在线视频| 成年女人毛片免费观看观看9 | 免费在线观看黄色视频的| 美女视频免费永久观看网站| 别揉我奶头~嗯~啊~动态视频 | 在线天堂中文资源库| 色视频在线一区二区三区| 亚洲av美国av| 青青草视频在线视频观看| 亚洲视频免费观看视频| 久久午夜综合久久蜜桃| 1024香蕉在线观看| 成人18禁高潮啪啪吃奶动态图| 国产精品国产av在线观看| 欧美另类一区| 国产av一区二区精品久久| 日韩中文字幕视频在线看片| 男女无遮挡免费网站观看| 欧美日韩亚洲国产一区二区在线观看 | av有码第一页| 亚洲av片天天在线观看| 一级毛片黄色毛片免费观看视频| 美国免费a级毛片| 国产熟女欧美一区二区| 欧美国产精品va在线观看不卡| 亚洲欧美精品自产自拍| 久久久久网色| 免费看十八禁软件| 日韩一区二区三区影片| 后天国语完整版免费观看| 建设人人有责人人尽责人人享有的| 免费看十八禁软件| 99热全是精品| 男人操女人黄网站| 成人国产一区最新在线观看 | 欧美97在线视频| 成人三级做爰电影| 日韩精品免费视频一区二区三区| 男女边吃奶边做爰视频| 欧美激情 高清一区二区三区| 久久国产精品人妻蜜桃| 人人妻人人爽人人添夜夜欢视频| 久久久精品94久久精品| av视频免费观看在线观看| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院| 人人澡人人妻人| 99国产精品免费福利视频| 一级毛片女人18水好多 | 久久99热这里只频精品6学生| 国产在线视频一区二区| 午夜日韩欧美国产| 欧美xxⅹ黑人| 老司机亚洲免费影院| 国产福利在线免费观看视频| 手机成人av网站| 国产精品久久久久久人妻精品电影 | 满18在线观看网站| 老司机在亚洲福利影院| 免费看十八禁软件| 免费在线观看影片大全网站 | 在线观看www视频免费| 男人操女人黄网站| 人人妻,人人澡人人爽秒播 | 精品人妻1区二区| 国产精品亚洲av一区麻豆| 精品亚洲成a人片在线观看| 中文字幕色久视频| 国产淫语在线视频| 精品第一国产精品| 日韩熟女老妇一区二区性免费视频| netflix在线观看网站| 成人免费观看视频高清| 精品久久久精品久久久| 欧美在线一区亚洲| 亚洲欧美日韩另类电影网站| 日本午夜av视频| 国产深夜福利视频在线观看| 精品欧美一区二区三区在线| 国产野战对白在线观看| 高清av免费在线| 久久毛片免费看一区二区三区| 在线观看国产h片| 男女床上黄色一级片免费看| 大片免费播放器 马上看| 成在线人永久免费视频| 亚洲欧美中文字幕日韩二区| 国产色视频综合| 18禁观看日本| 操美女的视频在线观看| 男女午夜视频在线观看| 久久久国产一区二区| 亚洲国产看品久久| 9色porny在线观看| 美女高潮到喷水免费观看| 一级毛片我不卡| 成年美女黄网站色视频大全免费| 又紧又爽又黄一区二区| 男人操女人黄网站| 国产野战对白在线观看| 免费av中文字幕在线| 亚洲第一青青草原| 日日摸夜夜添夜夜爱| 精品国产一区二区三区久久久樱花| 欧美av亚洲av综合av国产av| 中文字幕人妻熟女乱码| 波多野结衣av一区二区av| 色婷婷av一区二区三区视频| 在线观看人妻少妇| 激情五月婷婷亚洲| www.av在线官网国产| 久久 成人 亚洲| 国产有黄有色有爽视频| 一级黄色大片毛片| 后天国语完整版免费观看| 男人舔女人的私密视频| 老司机午夜十八禁免费视频| 国产日韩欧美亚洲二区| 国产精品人妻久久久影院| 午夜免费男女啪啪视频观看| 啦啦啦视频在线资源免费观看| 大陆偷拍与自拍| 欧美日韩av久久| 大香蕉久久成人网| 嫩草影视91久久| 各种免费的搞黄视频| 欧美97在线视频| 亚洲国产精品一区三区| 亚洲精品美女久久av网站| 美国免费a级毛片| 亚洲av欧美aⅴ国产| 黄网站色视频无遮挡免费观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲九九香蕉| 欧美日韩av久久| 后天国语完整版免费观看| 乱人伦中国视频| 熟女少妇亚洲综合色aaa.| 久久亚洲国产成人精品v| 午夜福利视频精品| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 午夜福利在线免费观看网站| 日本av手机在线免费观看| 国产成人一区二区在线| 久久影院123| 亚洲视频免费观看视频| 欧美另类一区| 久久毛片免费看一区二区三区| svipshipincom国产片| 日韩电影二区| 欧美国产精品一级二级三级| 国产成人系列免费观看| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 国产女主播在线喷水免费视频网站| 一本大道久久a久久精品| 国产成人精品在线电影| 国产成人免费无遮挡视频| 一二三四社区在线视频社区8| 国产日韩一区二区三区精品不卡| 1024香蕉在线观看| 男女国产视频网站| 日韩视频在线欧美| 大片电影免费在线观看免费| 丝袜喷水一区| 国产精品亚洲av一区麻豆| 欧美精品高潮呻吟av久久| 首页视频小说图片口味搜索 | 十分钟在线观看高清视频www| 中文字幕制服av| 欧美日韩福利视频一区二区| 嫩草影视91久久| 国产精品秋霞免费鲁丝片| 丝瓜视频免费看黄片| 免费在线观看黄色视频的| 国产欧美日韩一区二区三区在线| 在线观看免费午夜福利视频| 狠狠精品人妻久久久久久综合| 日韩大片免费观看网站| kizo精华| 丝袜美足系列| 国产女主播在线喷水免费视频网站| 久久久久久免费高清国产稀缺| 日韩大片免费观看网站| 久久亚洲国产成人精品v| 欧美日韩国产mv在线观看视频| 精品福利永久在线观看| 18禁观看日本| 欧美中文综合在线视频| 岛国毛片在线播放| 一边摸一边抽搐一进一出视频| 美女主播在线视频| 国产精品久久久av美女十八| 亚洲情色 制服丝袜| 老司机影院毛片| 国产精品免费视频内射| 国产成人av教育| 精品一品国产午夜福利视频| 久久99热这里只频精品6学生| 满18在线观看网站| 欧美成人精品欧美一级黄| 久久久国产一区二区| 欧美日本中文国产一区发布| 9热在线视频观看99| 99精品久久久久人妻精品| 日韩精品免费视频一区二区三区| 国精品久久久久久国模美| 久久久久久久久久久久大奶| 成年av动漫网址| 亚洲国产欧美在线一区| 久久久久精品人妻al黑| 老司机影院成人| 欧美成狂野欧美在线观看| 午夜福利视频在线观看免费| 久久久久久久久免费视频了| 黄频高清免费视频| 久久久久久久大尺度免费视频| 高清av免费在线| videosex国产| 午夜两性在线视频| 日韩大片免费观看网站| 老汉色av国产亚洲站长工具| 中文乱码字字幕精品一区二区三区| 免费高清在线观看视频在线观看| 国产亚洲精品第一综合不卡| 高潮久久久久久久久久久不卡| 你懂的网址亚洲精品在线观看| 国产一级毛片在线| 久久久国产一区二区| 久久精品人人爽人人爽视色| 一级黄色大片毛片| 日韩大片免费观看网站| 在线精品无人区一区二区三| 手机成人av网站| 少妇人妻 视频| 亚洲成人国产一区在线观看 | 久久性视频一级片| 91成人精品电影| 9色porny在线观看| 永久免费av网站大全| 女人被躁到高潮嗷嗷叫费观| 一级毛片电影观看| 你懂的网址亚洲精品在线观看| 欧美日韩黄片免| 91麻豆av在线| 亚洲成色77777| 97人妻天天添夜夜摸| 两人在一起打扑克的视频| 精品人妻1区二区| av天堂久久9| 大话2 男鬼变身卡| 国产深夜福利视频在线观看| 在线观看国产h片| 在线看a的网站| 婷婷成人精品国产| 美女扒开内裤让男人捅视频| 九草在线视频观看| 日本欧美国产在线视频| 亚洲欧美成人综合另类久久久| 人妻 亚洲 视频| 精品久久久精品久久久| 美女主播在线视频| 欧美变态另类bdsm刘玥| 女人爽到高潮嗷嗷叫在线视频| 汤姆久久久久久久影院中文字幕| 91字幕亚洲| 午夜激情久久久久久久| videosex国产| 妹子高潮喷水视频| svipshipincom国产片|