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

    Dynamic Evaluation Strategies for Multiple Aircrafts Formation Using Collision and Matching Probabilities

    2021-04-13 10:48:04HongboZhaoYongmingWenSentangWuandJiaDeng
    IEEE/CAA Journal of Automatica Sinica 2021年4期
    關(guān)鍵詞:堅持原則資源庫原則

    Hongbo Zhao, Yongming Wen, Sentang Wu, and Jia Deng

    Abstract—Configuration evaluation is a key technology to be considered in the design of multiple aircrafts formation (MAF)configurations with high dynamic properties in engineering applications. This paper deduces the relationship between relative velocity, dynamic safety distance and dynamic adjacent distance of formation members, then divides the formation states into collision-state and matching-state. Meanwhile, probability models are constructed based on the binary normal distribution of relative distance and relative velocity. Moreover, configuration evaluation strategies are studied by quantitatively analyzing the denseness and the basic capabilities according to the MAF collision-state probability and the MAF matching-state probability, respectively. The scale of MAF is grouped into 5 levels, and previous lattice-type structures are extended into four degrees by taking the relative velocities into account to instruct the configuration design under complex task conditions. Finally,hardware-in-loop (HIL) simulation and outfield flight test results are presented to verify the feasibility of these evaluation strategies.

    I. INTRODUCTION

    AT present, there are many applications for multiple aircrafts formation (MAF) systems. For example, MAF can improve comprehensive combat effectiveness significantly by carrying out tasks cooperatively. The configuration evaluation of MAF plays an important role in cooperative task execution [1]–[5]. Many efforts have been devoted towards the development of control methods for MAF, such as the leader-follower method [6]–[10], the virtual leader/virtual structure method [11], [12], the behavior-based method[13]–[15], consensus control [16]–[20] and neural network control [21], [22]. Towards realizing an actual MAF system,numerous interesting results about reducing the communication burden are proposed in [23]–[26]. However, few studies have paid attention to the configuration evaluation of the MAF system.

    Many studies have used qualitative descriptions, such as“l(fā)arge-scale”, “complex formation”, “close formation”, and“multilayer formation”, to represent the tightness and scale of the formation configuration. For example, a new formation framework of “l(fā)arge-scale” intelligent autonomous vehicles was developed in [27], which can realize “complex formations” while reducing data exchange. A relative positioning solution for unmanned aerial vehicle (UAV)rendezvous was proposed in [28] and the “close formation”flight had been verified in multiple flight tests. A multilayer formation control problem was defined in [29], in which the agents can receive and transmit information among the layers.Previous studies only described the configuration denseness qualitatively rather than quantitatively. In order to give efficient guidance for the tightness and scale design of MAF,we classify the denseness of MAF configuration quantitatively based on the MAF collision-state probability model.

    The basic capability of the configuration, i.e., the formation forming and holding ability, is constrained by the MAF members’ maneuverability, the sensor measurement and control accuracy, the network service quality, the task environment, and the denseness and scale of configuration in engineering applications. There was almost no quantitative analysis of them. For example, the formation control problem of multiple quadrotors tasked to achieve aggressive trajectory tracking with prescribed formation patterns was of concern in[30]. In [31], the authors analyzed the formation’s stability and solved the decentralized multivehicle formation control problem. The authors in [32] used Dudek’s taxonomy to evaluate the performance of the autonomous UAV group cooperating in mission execution. A taxonomy was presented to classify the multi-agent systems according to communication, computation and other capabilities in [33]. Most studies provided the qualitative evaluation of the configuration’s basic capabilities, but it is difficult for designers to optimize the formation’s properties to maximize the efficiency-cost ratio in different conditions. Therefore, we evaluate the basic capability of MAF configuration quantitatively according to the MAF matching-state probability model.

    Moreover, MAF’s lattice-type structure is an important indicator that can reflect basic capabilities. The more rigid the lattice-type structure, the stronger the basic capabilities. Two lattice-type structures called α-lattice and quasi α-lattice were presented in [34], which only require the distance between formation members to be equal to or approximately equal to the fixed pre-set distance. However, the relative distance is always fluctuating near the pre-set distance; the MAF members are not necessarily isomorphic, thus the pre-set distances between two adjacent members are not always the same. In addition, the α-lattice, quasi α-lattice, and quasi σ-adjacent configuration designed in [1] were all defined according to the relative distance between two members in the formation but the velocity maintaining capability of MAF was not taken into consideration, thereby leading to a lack of completeness. Therefore, the other two lattice-type structures (the σ-adjacent configuration and the quasi σ-adjacent configuration) are complemented in this paper by considering relative velocity to reflect the dynamic properties of the formation configuration.

    Interesting results on real-time experimental evaluation of formation control methods have also been reported. Both numerical simulations and field tests with a flock of five quadcopter drones were presented in [35], showing that the proposed flocking algorithm performs well in both the convergence rate and its ability to track the desired path. A real-time experiment was implemented using a team of ground mobile robots in [36], and the experimental formation tracking error and the trajectory illustrated the effectiveness of the proposed methods. In [19], four iRobot Create robots were used in the experiment, and the results showed that the robots’linear and angular velocities are reduced to zero while achieving the desired formation and heading consensus. Most real-time experiments are based on low-speed vehicles such as ground robots and quadcopter drones, and evaluation strategies are based on static criticism. To the best of our knowledge, there are limited experimental results on the evaluation of formation design and containment of cooperative aircraft with high velocity and high dynamic. In order to ensure reliability and safety of high dynamics multiagent systems in the outfield work, the development of dynamic configuration evaluation strategy is of critical necessity.

    In a recent paper [37], a formation generation method was proposed, and an actual outfield experiment was carried out with six UAVs. In addition, outfield flight experiments were carried out to verify the effectiveness of the proposed control method in [38] and [39]. In [38], the result curves of two aircrafts’ flight paths, velocity, and relative distance showed that the formation collisions and the chain effect are avoided effectively. In this paper, the results of [37]–[39] are significantly extended by considering configuration evaluation strategies for a team of cooperative fixed-wing unmanned aircraft.

    Motivated by the above discussion, this study focuses on quantitative configuration evaluation strategies for MAF systems based on collision probability and matching probability. The main contributions and innovations are summarized as follows.

    1) A novel evaluation method for highly dynamic MAF systems is discussed in this paper, which considers both the relative distance and relative velocity between formation members. In order to quantitatively describe the evaluation indicators, the MAF matching-state probability model and the collision-state probability model are constructed based on Reynolds three criterions [40].

    2) The denseness (including the tightness and the scale) of the MAF configuration is one of the evaluation objects. It is evaluated quantitatively and explicitly according to the MAF collision-state probability for the first time in this paper. The MAF tightness and scale are classified into 3 and 5 levels,respectively.

    3) The basic capabilities (the ability to form and hold a specified lattice-type structure) of the MAF configuration is another evaluation subject. It is evaluated quantitatively and explicitly according to the MAF matching-state probability for the first time in this paper. A more complete lattice-type structure classification is presented in this paper. Accordingly,four levels of MAF basic capabilities are graded.

    4) Most studies involving MAF focuses on the theoretical analysis of the control approaches, whereas the evaluation technology proposed in this paper is motivated by engineering applications and has been verified by both hardware-in-loop(HIL) simulations and outfield flight tests.

    The rest of this paper is organized as follows. In Section II,the background, basic knowledge, symbol definitions, and assumptions of the MAF are given. In Section III, the division of MAF collision-state and matching-state are presented,along with their probability models. Configuration evaluation strategies based on the MAF collision-state probability model and MAF matching state probability model are given in Sections IV and V, respectively. In Section VI, HIL simulation and outfield experiment results are presented to verify the feasibility of these evaluation methods. Finally,concluding remarks are given in Section VII.

    II. PRELIMINARIES

    A. Background and Basic Knowledge

    more complicated and closer to engineering applications than that of MAF in 3D space. Therefore, based on the above analysis, the configuration evaluation of the flat MAF in a 2D horizontal plane will be studied in this paper.

    B. The Distribution of the Relative Vector

    Lemma 1 (Wu [1]): Suppose that the distance dij(t) is an ergodic stochastic process, according to the central limit theorem, under the combined effects of various indistinctive random factors, and dij(t) approximately obeys the normal distribution

    C. Assumptions [38]

    1) Members of MAF can obtain information about local absolute positions and velocities through the navigation system in real-time.

    Fig.1. Probability density surface of Dij(t).

    2) The communication network of MAF only needs to be connected in the proximity groups through the one-hop or multi-hop pattern.

    III. DYNAMIC MAF STATES PROBABILITY MODELING

    Based on the Reynolds three criterions [40] (speed matching criterion, mutual attraction criterion, and mutual exclusion criterion), the matching-state and the collision-state probability models are constructed in this section to indicate the dynamic mutual exclusion and speed matching processes,respectively. With these probability models, we can quantitatively describe the evaluation indicators in the following sections.

    In this section, the dynamic safety distance and dynamic adjacent distance are first introduced to divide the state of two flight members. After that, the states of two members are extended to the states of the dynamic proximity groups and then coming to the states of MAF.

    A. Dynamic Safety Distance

    When the relative velocity is 0, the static safety distance is dsij0, which is the sum of the sensor error and control error,the maneuverability, the network induced noises, and the task environment disturbances. When the relative distance is less than the static safety distance, there will be a risk of collision,and the aircraft need to take maneuver immediately to avoid collision [1].

    When there is a nonzero relative velocity, the static safety distance turns to be the dynamic safety distance dsij(t). If the relative velocity is greater than 0, i.e., the two formation members fly far away from each other, the safe distance is dsij0; if the relative velocity is less than 0, i.e., the two formation members fly towards each other, extra distance is required for acceleration. Therefore, the dynamic safe distance will increase correspondingly. The relationship between dsij(t)and relative velocity is as follows.

    According to the static safety distance dsij0, if

    then a collision occurs between xiand xj.

    When Vlsij≤Vij(t)<0, and there is a relative acceleration amij>0 between xiand xj, Vij(t) is made to accelerate to 0 as quickly as possible to avoid collision. amijis the maximum relative acceleration, which is determined by the maneuverability of the aircraft. The collision condition between xiand xjis

    dij(t)?Vij2(t)/(2amij)

    where Vlsijis the minimum safe relative velocity, which is determined by the static safety distance, the maximum relative acceleration, and the static adjacent distance. The static adjacent distance daij0is the maximum distance required for communication between xiand xj, and satisfies daij0≥dsij0.When dij(t)

    Considering the situation that the two flight members fly towards each other with the adjacent distance as the relative distance, they should make full use of their maneuverability to decelerate, then by substituting dij(t)=daij0into (8), the minimum relative velocity they can tolerant is

    Due to the limitation of daij0, when Vij(t)

    According to (7), (8), and (10), we get the relationship of dsij(t) and Vij(t)

    In summary, the collision condition between xiand xjis

    The relationship curve of the dynamic safety distance is shown as in Fig.2, i.e., the dividing curve of the green region and red region.

    Fig.2. Relationship of relative distance and relative velocity.

    B. Dynamic Adjacent Distance

    Allow the static adjacent distance to be daij0. Taking relative velocity into consideration, the static adjacent distance turns into the dynamic adjacent distance. If it satisfies

    xiand xjdo not need to communicate with each other, and they are in a state of mutual isolation. When 00 between xiand xjmaking Vij(t) decelerate to zero, according to (13),the condition under which xiand xjdo not need to communicate with is given by

    where Vlaijis the maximum adjacent relative velocity. By substituting dij(t)=dsij0into (14), we can get

    When Vij(t)>Vlaijand dij(t) is larger than the static safe distance dsij0, then xiand xjdo not need to communicate, i.e.,

    According to (13), (14), and (16), we can get the relationship of the adjacent distance daij(t) and relative velocity Vij(t)

    In summary, when

    xiand xjdo not need to communicate with each other, and they are in a state of mutual isolation. Thus, the network flow pressure can be reduced, and the MAF can use the limited network resources among the members in the dynamic proximity groups to improve the network service quality,which is conducive to the safety, stability, and efficiency of MAF.

    The relationship of relative distance and relative velocity is illustrated in Fig.2, where ?ijdenotes the relationship plane.The dividing curve of the green region and the red region is the safety distance in (11). The dividing curve of the green region and the blue region is the adjacent distance in (17).

    C. Two Formation Member’s States Division

    Definition 1 (Collision-State): In the relative distance and relative velocity relationship plane shown in Fig.2, the area constrained by (12), i.e., the red region Ψij

    is defined as the collision-area. xiand xjare in the collisionstate when the relative vector Dij(t) is in the collision-area Φij.

    According to (6), (11), and (19) we can get collision probability p(Φij) of xiand xjas follows:

    where P{·} denotes the probability. τ and ω are the integral symbols of dij(t) and Vij(t), respectively.

    Definition 2 (Isolation-State): In the distance and relative velocity plane ?ijshown in Fig.2, the plane constrained by(18), i.e., the blue area Ψij

    Ψij={(dij(t),Vij(t)):dij(t)≥daij(t)} (21)

    is defined as the isolation-area. xiand xjare in the isolationstate when the relative vector Dij(t) is in the isolation-area Ψi j.

    Definition 3 (Matching-State [2]): In the distance and relative velocity plane ?ijshown in Fig.2, after removing the collision-area and the isolation-area, the remaining plane, i.e.,the green area Πij

    在開發(fā)及利用課程資源的過程中不能夠任性妄為,而堅持原則能夠?qū)θ藗兊男袨檫M(jìn)行一定的規(guī)范,指導(dǎo)人們正確的行為,在開發(fā)以及利用課程資源的過程中也需要一定的原則來進(jìn)行規(guī)范。在收集到的資料中大多沒有形成規(guī)矩沒有順序,不能夠直接簡單、盲目地投放到資源庫中或直接使用。所以,為了能夠充分地、合理地利用收集到的資料,應(yīng)該進(jìn)行篩選。在篩選的過程中應(yīng)該遵循以下原則:個性化原則、發(fā)展性原則、生活性原則、適應(yīng)性原則、科學(xué)性原則。

    is defined as the matching-area. When the relative vector Dij(t) is in the matching-area, i.e., Dij(t)∈Πij, we say that xiand xjare in the matching-state, which means that xiand xjneed to communicate with each other and there will be no collision between them.

    Based on (6), (11), (17), and (22), we get the matching-state probability p(Πij) between xiand xjas follows:

    By combining Figs. 1 and 2, we can get the probability division of the three states as shown in Fig.3.

    Fig.3. Probability division of states between formation members.

    From Fig.3, we can observe that the probabilities of collision-state, matching-state, and isolation-state of flight members are the volumes formed by collision-area, matchingarea, and isolation-area on the probability density surface of Dij(t) in Fig.1, respectively.

    D. MAF States Division

    Definition 4 (Dynamic Adjacency Degree [2]): The dynamic adjacency degree naiof xiis the number of adjacent neighbors that are in the matching-state with xi. The dynamic proximity group of xiis

    where the group members satisfy

    Fig.10. MAF initial configurations and turning configurations of different dynamic adjacent degree. (a) line; (b) diamond; (c) two-line; (d) hexagon; (e) line turning; (f) diamond turning; (g) two-line turning; (h) hexagon turning.

    basic capability of MAF. Moreover, the quasi σ-adjacent is the most widely used one in practical application and is more suitable to describe the lattice-type of MAF.

    VI. HIL EXPERIMENT AND OUTFIELD FLIGHT TESTS

    A. System Composition

    The system overview of one HIL flight member is shown in Fig.9(a), “NET” denotes network port communication,“COM” is serial communication, the green arrows denote the data link of formation flight control, and the red arrows denote the data link of member flight control. The HIL experiment system formed by seven aircraft is shown in Fig.9(b).

    It should be noted that each member of MAF consists of the formation control module, the flight control module, the flight simulation module, the virtual communication module, and the flight monitoring module, but not including the formation monitoring module.

    B. Comparative Experiment and Result Analysis of Different Dynamic Adjacent Degree

    With the above HIL system, we design the line configuration (Fig.10(a)), the diamond configuration (Fig.10(b)), the two-line configuration (Fig.10(c)), and the hexagon configuration (Fig.10(d)), where the red dot and the white dot represent the central member and the adjacent members in the dynamic proximity group, respectively, and the numbers before and after the colon indicate the member number and flight altitude, respectively.

    In order to verify the influence of different dynamic adjacent degrees on the denseness and the basic capability of MAF configuration, we observe the collision probability and the configuration holding capability of the seven aircraft at it turns as shown in Figs. 10(e)–10(h).

    As shown in Fig.10(a), the dynamic adjacent degree of A4(A represents aircraft) is na4=2 and its adjacent members are A3 and A5; in Fig.10(b), the dynamic adjacent degree of A4 is na4=3 and its adjacent members are A0, A3, and A5; in Fig.10(c), the dynamic adjacent degree of A6 is na6=4 and its adjacent members are A1, A3, A4, and A5; in Fig.10(d),the dynamic adjacent degree of A6 is na6=6 and its adjacent members are A0, A1, A2, A3, A4, and A5.

    We can observe from Figs. 10(e)–10(h) that: i) when MAF is turning, the configuration on the inner side of the corner is squeezed, and the collision risks of the dynamic proximity group increase. The greater the dynamic adjacent degree, the smaller the safe space that can be used to avoid a collision, the greater the probability of dynamic proximity group collisionstate and the denser the MAF configuration. In other words,the hexagon configuration is the densest while having higher risks of being in the collision state, while the line configuration is the loosest and has more maneuver space to avoid collision; ii) when MAF is in turning, the smaller the dynamic adjacent degree of dynamic proximity group is, the better the configuration holds and the better the basic capability of MAF is. That is to say, the hexagon configuration owns the worse configuration holding capacity while the line configuration holds the best.

    The dynamic adjacent group’s collision probability and matching probability are calculated based on the results in Section III, the denseness and the basic capability classification are given based on the results in Tables I and III.Specific comparison data are shown in Table IV.

    Fig.11. MAF initial hexagon configurations and turning configurations of different expected distances. (a) μij=50.0 m ; (b) μij=60.0 m ; (c) μi j=74.0 m;(d) μ ij=90.0 m ; (e) turning of μ ij=50.0 m ; (f) turning of μ ij=60.0 m ; (g) turning of μ ij=74.0 m ; (h) turning of μ ij=90.0 m.

    TABLE IV COMPARISON DATA OF DIFFERENT DYNAMIC ADJACENT DEGREE

    In conclusion, if the task requires a relatively complex configuration with a higher adjacent degree like the hexagon configuration, the platform should be equipped with higher maneuverability to adjust the relative distance and velocity with higher frequency to avoid a collision. Due to the frequent adjustment in this case, the commander needs to be more tolerant with the bad MAF configuration holding ability.

    C. Comparative Experiment and Result Analysis of Different Expected Distances

    In order to verify the influence of different expected distances on the denseness and the basic capability of MAF configuration, we designed the hexagon configuration for MAF with different expected distances under the instruction of Result 1. In this paper, we chose to present four representative distances: μij=50.0 m (Fig.11(a), λdij0=0.0),μij=60.0 m (Fig.11(b), λdij0=2.0 ), μij=74.0 m (Fig.11(c),λdij0=4.8) , and μij=90.0 m (Fig.11(d), λdij0=8.0). We observe the collision probability and the configuration holding degree of the dynamic proximity group at the turn of MAF shown in Figs. 11(e)–(h).

    We can observe from Fig.11 that: i) when MAF is turning,the configuration on the inner side of the corner is squeezed,and the collision risk of the dynamic proximity group increases. The smaller the expected distance (or the static distance safe factor), the smaller the safe space that can be used to avoid collision, the greater the probability of dynamic proximity group collision-state and the denser the MAF configuration; ii) when MAF is turning, the closer to the expected distance, the better the configuration holds and the basic capability of MAF configuration improves. We got the expected distance 74 m after many tests which corresponding to a good basic capability of hexagon configuration, and the matching-state probability is relatively high compared with others.

    The collision probabilities and matching probabilities of the hexagon configuration with different expected distances are calculated based on the results in Section III, and their denseness and basic capability are graded based on the results in Tables I and III. Specific comparison data are shown in Table V.

    TABLE V COMPARISON DATA OF DIFFERENT EXPECTED DISTANCE

    The results reach a reasonable evaluation of MAF denseness and basic capability. It also arrives at the same conclusion as Result 4 that μijshould take the intermediate value of safety distance and adjacent distance. To be more specific, when μijis less than the safety distance, the collision probability between two formation members will increase; when μijis larger than adjacent distance, the isolation probability will increase.

    D. Outfield Flight Test

    In order to guide the engineering applications better, we carried out MAF outfield flight tests after verifying the reasonability and safety of the designed configuration and expected relative distance in the HIL experiment. The denseness was evaluated quantitatively by the proposed method.

    Fig.12. Outfield flight test equipment. (a) Six UAVs; (b) The loads in the cabin of the UAV.

    Fig.13. MAF flight. (a) MAF flocking (before turning); (b) MAF flocking (turning).

    The experiment is composed of six unmanned aircraft as shown in Fig.12(a). One of them serves as the aerial photographic UAV to capture the configuration. The other five of them are equipped with the power module, the autopilot module, the formation cooperative guidance module,the formation communication module, the GPS, etc.(Fig.12(b)).

    The expected flight height is 200.0 m, the expected velocity is 25 m/s, the static safety distance is 50 m and the adjacent distance is 100 m. We set the height of the aerial photographic UAV to be 400 meters to capture the configuration.

    The specific comparison data of these MAF flight experiments are shown in Table VI.

    TABLE VI DATA FOR MAF FLIGHT EXPERIMENTS

    Note that we have applied collision avoidance algorithms to the formation control in outfield flight tests to keep the distance between two flight members always larger than the dynamic safe distance to ensure the safety of aircrafts.

    E. Experiment Conclusion

    In summary, we verify the influences of the dynamic adjacent degree naiand the expected distance μij(or the static distance safe factor λdij0) on the denseness and the basic capability of MAF configuration respectively by the HIL simulation experiment and outfield flight test, and the experiment results are the same as the conclusions in Sections IV and V. Quantitative evaluation of formation collision-state and matching-state provides quantitative evaluation indicators for engineering practices. The feasibility of these two evaluation methods of MAF configuration is thus verified.

    VII. CONCLUSION

    In this paper, with the aim to evaluate high dynamic MAF configuration, the dynamic safety distance and dynamic adjacent distance which are related with relative velocity are introduced to divide the MAF states. Based on the division,the collision-state and matching-state probability models are built. Then the two evaluation methods of MAF configuration are studied: a) The denseness of MAF configuration,including the tightness (dense, tight, and loose configuration)and the scale (massive, large, medium, small, and adjacent configuration), are evaluated quantitatively by the probability of MAF collision-state. The influences of the static safety distance, the maximum relative acceleration, the distance standard deviation, the relative velocity standard deviation,the correlation coefficient of distance and relative velocity, the expected distance, and the expected relative velocity on the collision-state probability are analyzed. The methods to reduce the probability of MAF collision-state to strengthen the formation safety are given according to the evaluation. The MAF tightness and scale are graded based on the collisionstate probability; b) The basic capability of MAF configuration is evaluated quantitatively by the probability of MAF matching-state. In this process, the influences of the same factors in a) on the matching-state probability are analyzed, and the methods to increase the probability of MAF matching-state to strengthen the basic capability of MAF are given. In addition, based on the relative distance and the relative velocity, two lattice-type structures ( σ-adjacent configuration and the quasi σ-adjacent configuration) are compensated for the previous studies to grade the configuration basic capability, so that the lattice-type structure of configuration becomes more complete. Finally, hardwarein-the-loop simulation and outfield flight test results are presented to verify the rationality of these evaluation methods.

    猜你喜歡
    堅持原則資源庫原則
    健身氣功開放課程資源庫建設(shè)研究
    貴州●石斛種質(zhì)資源庫
    堅守原則,逐浪前行
    女漢子的新標(biāo)準(zhǔn)
    知識窗(2016年12期)2016-12-13 17:57:58
    高中歷史信息化教育資源庫應(yīng)用探索
    福建基礎(chǔ)教育教學(xué)資源庫建設(shè)研究——以福建基礎(chǔ)教育網(wǎng)資源庫為例
    無罪推定原則的理解與完善
    堅持原則 凈化聯(lián)壇——潘濤會長捍衛(wèi)《聯(lián)律通則》二三事
    對聯(lián)(2011年20期)2011-09-19 06:24:36
    惹人喜愛的原則(二)
    知識窗(2010年9期)2010-05-14 09:07:50
    惹人喜愛的原則
    知識窗(2010年8期)2010-05-14 09:07:49
    久久久精品大字幕| 亚洲真实伦在线观看| 国产精品久久久人人做人人爽| 亚洲性夜色夜夜综合| 亚洲成av人片免费观看| 国产又黄又爽又无遮挡在线| 亚洲国产高清在线一区二区三| 最新中文字幕久久久久 | 国产精品久久电影中文字幕| 亚洲最大成人中文| 巨乳人妻的诱惑在线观看| 国产精品日韩av在线免费观看| 这个男人来自地球电影免费观看| 观看免费一级毛片| 成人特级黄色片久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 免费观看人在逋| 日本撒尿小便嘘嘘汇集6| 日日摸夜夜添夜夜添小说| 一级黄色大片毛片| 99热只有精品国产| 亚洲五月婷婷丁香| 亚洲天堂国产精品一区在线| 精品国产超薄肉色丝袜足j| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 亚洲国产精品久久男人天堂| 国产精品久久久久久久电影 | 国模一区二区三区四区视频 | 日韩三级视频一区二区三区| 在线免费观看不下载黄p国产 | 亚洲欧美一区二区三区黑人| 婷婷丁香在线五月| 国产精品99久久99久久久不卡| 国产精品久久久av美女十八| 国产精品久久久久久人妻精品电影| 欧美中文日本在线观看视频| 国产毛片a区久久久久| 日韩成人在线观看一区二区三区| 香蕉丝袜av| av片东京热男人的天堂| 熟妇人妻久久中文字幕3abv| 国产又色又爽无遮挡免费看| 亚洲黑人精品在线| 欧美极品一区二区三区四区| 超碰成人久久| 搡老岳熟女国产| 国产伦人伦偷精品视频| 国产成人啪精品午夜网站| 757午夜福利合集在线观看| 日韩成人在线观看一区二区三区| 日本 欧美在线| 欧美中文综合在线视频| 欧美不卡视频在线免费观看| 国产91精品成人一区二区三区| 中国美女看黄片| 在线播放国产精品三级| 亚洲中文字幕一区二区三区有码在线看 | 最近在线观看免费完整版| 九色国产91popny在线| 黑人巨大精品欧美一区二区mp4| 久久99热这里只有精品18| 91久久精品国产一区二区成人 | 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 免费av毛片视频| 麻豆av在线久日| 国产精品永久免费网站| 在线a可以看的网站| 国产精品久久久久久精品电影| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 亚洲精品在线美女| 哪里可以看免费的av片| 午夜激情福利司机影院| 免费看光身美女| 18禁黄网站禁片免费观看直播| 一夜夜www| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| av天堂在线播放| 舔av片在线| 国产日本99.免费观看| 搡老熟女国产l中国老女人| 午夜福利18| 91在线精品国自产拍蜜月 | 亚洲av五月六月丁香网| 国产精品免费一区二区三区在线| 可以在线观看的亚洲视频| 丁香六月欧美| 久久久久亚洲av毛片大全| 免费看美女性在线毛片视频| 99久久综合精品五月天人人| av天堂中文字幕网| 长腿黑丝高跟| 很黄的视频免费| 在线观看舔阴道视频| 亚洲人成网站高清观看| 久久精品国产综合久久久| 亚洲午夜精品一区,二区,三区| 不卡一级毛片| 久久婷婷人人爽人人干人人爱| 长腿黑丝高跟| www.www免费av| 老司机午夜十八禁免费视频| 久久久精品大字幕| 变态另类成人亚洲欧美熟女| 欧美三级亚洲精品| 最新中文字幕久久久久 | 亚洲人成伊人成综合网2020| 亚洲第一电影网av| 亚洲欧洲精品一区二区精品久久久| 日日夜夜操网爽| 特级一级黄色大片| 午夜视频精品福利| 欧美午夜高清在线| 精品免费久久久久久久清纯| 亚洲人与动物交配视频| 欧美大码av| 亚洲av成人一区二区三| 激情在线观看视频在线高清| 亚洲成人久久性| 神马国产精品三级电影在线观看| 亚洲国产欧洲综合997久久,| 午夜福利在线观看吧| 在线观看免费午夜福利视频| 精品久久久久久久人妻蜜臀av| www.自偷自拍.com| 亚洲国产精品成人综合色| 国内久久婷婷六月综合欲色啪| 夜夜看夜夜爽夜夜摸| 国产高清视频在线观看网站| 桃色一区二区三区在线观看| 国产精品 国内视频| 午夜成年电影在线免费观看| 不卡一级毛片| 亚洲精品中文字幕一二三四区| 成人亚洲精品av一区二区| 一个人看的www免费观看视频| 他把我摸到了高潮在线观看| 女警被强在线播放| 欧美激情久久久久久爽电影| 在线免费观看的www视频| 亚洲人成电影免费在线| 亚洲最大成人中文| 99精品欧美一区二区三区四区| 麻豆成人午夜福利视频| 亚洲成av人片在线播放无| 美女大奶头视频| 国产精品 国内视频| 久久久久久久久免费视频了| 午夜福利18| 亚洲国产欧美一区二区综合| 精品久久久久久久人妻蜜臀av| 亚洲自拍偷在线| 日本 欧美在线| 99热这里只有是精品50| av天堂在线播放| 久久久久久久精品吃奶| 天堂av国产一区二区熟女人妻| 村上凉子中文字幕在线| 丝袜人妻中文字幕| 91在线观看av| 国产私拍福利视频在线观看| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 中文字幕人妻丝袜一区二区| 狂野欧美激情性xxxx| 国语自产精品视频在线第100页| 99热只有精品国产| 亚洲第一欧美日韩一区二区三区| 亚洲国产高清在线一区二区三| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 一级黄色大片毛片| 偷拍熟女少妇极品色| 欧美中文综合在线视频| 18禁观看日本| 午夜福利在线观看吧| 成人欧美大片| 久久久成人免费电影| 毛片女人毛片| 国产毛片a区久久久久| 国产精品免费一区二区三区在线| 欧美丝袜亚洲另类 | 日本黄色视频三级网站网址| 午夜免费观看网址| 嫩草影院入口| 国产精品 欧美亚洲| 国产亚洲精品一区二区www| 国产一级毛片七仙女欲春2| 日韩 欧美 亚洲 中文字幕| 久久久久久大精品| 99热6这里只有精品| 91老司机精品| 1000部很黄的大片| 99精品欧美一区二区三区四区| 一夜夜www| 母亲3免费完整高清在线观看| 欧美zozozo另类| 一区福利在线观看| 深夜精品福利| 欧美性猛交╳xxx乱大交人| 欧美中文日本在线观看视频| 少妇熟女aⅴ在线视频| 丁香欧美五月| 亚洲国产精品成人综合色| 国产高清videossex| 日韩精品中文字幕看吧| 久久精品夜夜夜夜夜久久蜜豆| 国产精品 国内视频| 久久人妻av系列| 亚洲精品在线美女| 99在线人妻在线中文字幕| 99久久国产精品久久久| 国产激情偷乱视频一区二区| 欧美一区二区精品小视频在线| 亚洲18禁久久av| 久久九九热精品免费| 欧美一级a爱片免费观看看| av福利片在线观看| 亚洲国产精品999在线| 最近视频中文字幕2019在线8| 99热6这里只有精品| 欧美高清成人免费视频www| 国产精品影院久久| 嫩草影院精品99| 2021天堂中文幕一二区在线观| 综合色av麻豆| 午夜激情福利司机影院| 亚洲av成人一区二区三| 免费看a级黄色片| 欧美性猛交╳xxx乱大交人| 最新在线观看一区二区三区| 特级一级黄色大片| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品合色在线| 亚洲精品美女久久久久99蜜臀| 欧美成人免费av一区二区三区| 1000部很黄的大片| 十八禁人妻一区二区| 好看av亚洲va欧美ⅴa在| 国产成年人精品一区二区| 99久久综合精品五月天人人| 人人妻人人澡欧美一区二区| 1024手机看黄色片| 国产精品久久久久久久电影 | bbb黄色大片| 老司机在亚洲福利影院| 18禁美女被吸乳视频| 脱女人内裤的视频| 国产免费av片在线观看野外av| 色视频www国产| 国内精品美女久久久久久| 岛国在线免费视频观看| 特级一级黄色大片| 久久香蕉国产精品| 高潮久久久久久久久久久不卡| 亚洲国产看品久久| 成年免费大片在线观看| 精品电影一区二区在线| 高清毛片免费观看视频网站| 又黄又爽又免费观看的视频| 黄频高清免费视频| 成人精品一区二区免费| 国产91精品成人一区二区三区| 狂野欧美激情性xxxx| 精品电影一区二区在线| 午夜免费激情av| 欧美激情久久久久久爽电影| 性色av乱码一区二区三区2| 99精品在免费线老司机午夜| 97超视频在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一区av在线观看| 成年人黄色毛片网站| 成人国产综合亚洲| 观看美女的网站| 窝窝影院91人妻| 一进一出好大好爽视频| 午夜成年电影在线免费观看| 在线观看一区二区三区| 日本撒尿小便嘘嘘汇集6| 1000部很黄的大片| 亚洲色图 男人天堂 中文字幕| 亚洲av五月六月丁香网| 国产成人精品久久二区二区91| 真人做人爱边吃奶动态| 女生性感内裤真人,穿戴方法视频| 一本精品99久久精品77| 国产伦精品一区二区三区视频9 | 亚洲人成网站高清观看| 丰满人妻一区二区三区视频av | 亚洲 欧美一区二区三区| 99久久精品热视频| 欧美另类亚洲清纯唯美| 色噜噜av男人的天堂激情| 97超级碰碰碰精品色视频在线观看| 麻豆久久精品国产亚洲av| a级毛片在线看网站| 老司机深夜福利视频在线观看| 欧美一区二区国产精品久久精品| 一级毛片精品| 久久香蕉国产精品| 精品国产乱子伦一区二区三区| 九九久久精品国产亚洲av麻豆 | 日本五十路高清| 成人av在线播放网站| 国产亚洲av高清不卡| 亚洲黑人精品在线| 白带黄色成豆腐渣| 91久久精品国产一区二区成人 | 日日夜夜操网爽| 亚洲精品一区av在线观看| 噜噜噜噜噜久久久久久91| 成人鲁丝片一二三区免费| 免费av不卡在线播放| 九色成人免费人妻av| 99久国产av精品| 我的老师免费观看完整版| 亚洲七黄色美女视频| 亚洲av成人精品一区久久| 日韩欧美精品v在线| 午夜福利视频1000在线观看| 久久精品国产综合久久久| 1024香蕉在线观看| x7x7x7水蜜桃| 国产黄片美女视频| 成年女人永久免费观看视频| 亚洲男人的天堂狠狠| 国产爱豆传媒在线观看| 国产一区二区三区在线臀色熟女| svipshipincom国产片| 一进一出抽搐gif免费好疼| 男女视频在线观看网站免费| 色综合欧美亚洲国产小说| 亚洲国产高清在线一区二区三| 亚洲无线在线观看| 日本黄色片子视频| 757午夜福利合集在线观看| 免费av毛片视频| 国内精品久久久久久久电影| 国产黄片美女视频| 女生性感内裤真人,穿戴方法视频| 成年女人毛片免费观看观看9| 国产视频一区二区在线看| or卡值多少钱| 高清在线国产一区| 男人舔奶头视频| 中亚洲国语对白在线视频| 小蜜桃在线观看免费完整版高清| 国产高清激情床上av| 国产淫片久久久久久久久 | 亚洲成av人片免费观看| 999精品在线视频| 中国美女看黄片| 婷婷精品国产亚洲av| 亚洲国产欧美一区二区综合| 99视频精品全部免费 在线 | 免费看a级黄色片| 在线看三级毛片| 久久欧美精品欧美久久欧美| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| av片东京热男人的天堂| 91av网一区二区| 成年女人毛片免费观看观看9| 美女大奶头视频| 成年女人毛片免费观看观看9| 国产成人系列免费观看| 99热这里只有精品一区 | 色视频www国产| 午夜影院日韩av| 在线看三级毛片| 最好的美女福利视频网| 最新在线观看一区二区三区| 99re在线观看精品视频| 亚洲精品乱码久久久v下载方式 | 禁无遮挡网站| 麻豆久久精品国产亚洲av| 国产淫片久久久久久久久 | 精品久久久久久久人妻蜜臀av| 国产亚洲欧美在线一区二区| 叶爱在线成人免费视频播放| 亚洲国产欧美网| 欧美日韩中文字幕国产精品一区二区三区| 亚洲成av人片在线播放无| or卡值多少钱| 日韩国内少妇激情av| 91在线精品国自产拍蜜月 | 日本 av在线| svipshipincom国产片| 在线永久观看黄色视频| 久久热在线av| 国产一区二区激情短视频| 国产激情偷乱视频一区二区| 国产精品av久久久久免费| 国产伦在线观看视频一区| 亚洲九九香蕉| 又黄又爽又免费观看的视频| 在线观看午夜福利视频| 黄色 视频免费看| 欧美日韩一级在线毛片| 日本五十路高清| 一区二区三区高清视频在线| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 欧美+亚洲+日韩+国产| 精品一区二区三区四区五区乱码| 精品人妻1区二区| 在线看三级毛片| 亚洲成人精品中文字幕电影| 欧美极品一区二区三区四区| 免费电影在线观看免费观看| 啦啦啦韩国在线观看视频| 19禁男女啪啪无遮挡网站| 亚洲av中文字字幕乱码综合| 欧美日韩乱码在线| 亚洲欧美精品综合一区二区三区| 欧美乱妇无乱码| 久久久久久九九精品二区国产| 99热6这里只有精品| 色在线成人网| 日本三级黄在线观看| 91麻豆精品激情在线观看国产| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 成人精品一区二区免费| 51午夜福利影视在线观看| 亚洲午夜精品一区,二区,三区| 18禁国产床啪视频网站| 国产 一区 欧美 日韩| АⅤ资源中文在线天堂| 人人妻人人看人人澡| 少妇裸体淫交视频免费看高清| 在线十欧美十亚洲十日本专区| 国产精品久久久久久人妻精品电影| 小说图片视频综合网站| 欧美日韩一级在线毛片| 无遮挡黄片免费观看| 亚洲男人的天堂狠狠| 1024香蕉在线观看| 国产伦精品一区二区三区四那| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| 最近最新免费中文字幕在线| 两个人的视频大全免费| 日本黄色片子视频| 久久婷婷人人爽人人干人人爱| 亚洲av成人精品一区久久| 色综合站精品国产| 亚洲va日本ⅴa欧美va伊人久久| 免费人成视频x8x8入口观看| 国产99白浆流出| 久久这里只有精品19| 一级a爱片免费观看的视频| 国内久久婷婷六月综合欲色啪| 国产成年人精品一区二区| 麻豆一二三区av精品| 性色av乱码一区二区三区2| 黄色日韩在线| 1024手机看黄色片| 午夜福利18| 国产精品 欧美亚洲| 国产精品一区二区三区四区久久| 日韩欧美在线乱码| 啪啪无遮挡十八禁网站| 午夜福利在线观看免费完整高清在 | 一区福利在线观看| 精品国产乱子伦一区二区三区| 亚洲狠狠婷婷综合久久图片| 精品福利观看| 亚洲欧美日韩高清专用| 午夜影院日韩av| 国产久久久一区二区三区| 日本精品一区二区三区蜜桃| 亚洲成av人片在线播放无| 99国产精品一区二区蜜桃av| 亚洲一区二区三区色噜噜| 亚洲人成伊人成综合网2020| 国产精品av视频在线免费观看| 怎么达到女性高潮| 欧美日韩乱码在线| 日韩中文字幕欧美一区二区| 女警被强在线播放| 一区福利在线观看| 美女被艹到高潮喷水动态| 午夜亚洲福利在线播放| 夜夜看夜夜爽夜夜摸| 在线免费观看不下载黄p国产 | 亚洲欧美日韩高清在线视频| 香蕉丝袜av| 精品不卡国产一区二区三区| 日韩有码中文字幕| 美女免费视频网站| 真实男女啪啪啪动态图| 亚洲av电影在线进入| 日韩精品中文字幕看吧| 亚洲欧美日韩高清专用| 亚洲国产精品999在线| 最新美女视频免费是黄的| 露出奶头的视频| 色老头精品视频在线观看| 久久久色成人| 97人妻精品一区二区三区麻豆| 一本精品99久久精品77| 精品久久久久久久毛片微露脸| 日韩精品中文字幕看吧| 别揉我奶头~嗯~啊~动态视频| 美女大奶头视频| 黑人欧美特级aaaaaa片| 精品国产乱码久久久久久男人| 久久国产精品影院| 久久这里只有精品19| 草草在线视频免费看| 窝窝影院91人妻| 成人午夜高清在线视频| 12—13女人毛片做爰片一| 久久久国产成人精品二区| 欧美午夜高清在线| 91在线观看av| 十八禁人妻一区二区| 人人妻,人人澡人人爽秒播| 国产精品永久免费网站| 淫秽高清视频在线观看| 国产精品一区二区精品视频观看| 99久久国产精品久久久| 黄色女人牲交| 小蜜桃在线观看免费完整版高清| 日本 欧美在线| 嫩草影视91久久| 亚洲国产中文字幕在线视频| 国产成人福利小说| 亚洲色图av天堂| 女警被强在线播放| 男女做爰动态图高潮gif福利片| 丁香六月欧美| 亚洲五月天丁香| 国产乱人视频| 国产久久久一区二区三区| 欧美一区二区国产精品久久精品| 黄色日韩在线| 在线看三级毛片| 18禁国产床啪视频网站| 国产精品亚洲美女久久久| 禁无遮挡网站| 五月伊人婷婷丁香| 可以在线观看毛片的网站| 亚洲黑人精品在线| www.自偷自拍.com| 中文字幕最新亚洲高清| 天天添夜夜摸| 久久精品91蜜桃| 中出人妻视频一区二区| 国产成人福利小说| 网址你懂的国产日韩在线| 午夜激情福利司机影院| 91麻豆精品激情在线观看国产| 香蕉久久夜色| 免费观看的影片在线观看| 亚洲欧美日韩无卡精品| 黑人欧美特级aaaaaa片| 日日干狠狠操夜夜爽| 国产爱豆传媒在线观看| 国产精品永久免费网站| 国产精品av久久久久免费| 国产午夜福利久久久久久| 亚洲欧美精品综合久久99| 国产高清视频在线观看网站| 日韩精品中文字幕看吧| 国产成人一区二区三区免费视频网站| 久久久成人免费电影| 亚洲一区二区三区不卡视频| 美女cb高潮喷水在线观看 | 国产成人av激情在线播放| 最近最新中文字幕大全免费视频| 999久久久国产精品视频| 国产精品永久免费网站| 亚洲九九香蕉| 亚洲成人久久性| 国产久久久一区二区三区| 成在线人永久免费视频| 日本黄大片高清| 久久精品国产99精品国产亚洲性色| 国产高清视频在线观看网站| 激情在线观看视频在线高清| 老司机午夜十八禁免费视频| 日韩欧美三级三区| 又黄又爽又免费观看的视频| 哪里可以看免费的av片| h日本视频在线播放| 97超级碰碰碰精品色视频在线观看| 日本精品一区二区三区蜜桃| 手机成人av网站| 欧美色欧美亚洲另类二区| 国产高清视频在线观看网站| 亚洲av五月六月丁香网| 欧美色欧美亚洲另类二区| 国产高清视频在线观看网站| 网址你懂的国产日韩在线| 成年免费大片在线观看| 变态另类丝袜制服| av福利片在线观看| 国产成人欧美在线观看| 国产亚洲欧美在线一区二区| 亚洲中文字幕日韩| 老司机福利观看| 天天添夜夜摸| 亚洲成人久久性| 中文字幕最新亚洲高清| 91在线精品国自产拍蜜月 | 在线观看美女被高潮喷水网站| 我的老师免费观看完整版|