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

    A novel variable topology design for a multi-flexible ejection mechanism

    2020-05-23 07:09:28HoLiuJunZhouZongliMoShiweiZhng
    Defence Technology 2020年2期

    Ho Liu , Jun Zhou ,*, Zong-li Mo , Shi-wei Zhng

    a School of Astronautics, Northwestern Polytechnical University, Xi'an, 710072, China

    b School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China

    c Luoyang Opto-Electro Technology Development Center, 471009, China

    Keywords:Variable topology Gap dynamics Multi-flexible body Flexible control Ejection mechanism

    ABSTRACT The airborne missile launch mechanism often subjects to significant deformations induced by the large ejection force during high-speed actuations in missile ejection process,leading to a substantial deviation of separation parameters from designed values that threats safety of the carrier. This study proposes a novel variable topology design for launch mechanism, achieved via a Prismatic-Revolute-Revolute pair(PR-R) motion formed by the structural gap with a specific direction. It enables launch mechanism variability during missile ejection process and optimizes the ejection force given by the front and back ejection arms, and greatly optimizes the separation parameters during missile ejection. The kinetics simulation analysis is conducted under working conditions of the original ejection mechanism and the novel mechanism with variable topology design, respectively. The results show that the novel variable topology design is more befitting for the launch process in terms of system safety and controllability,effectively improving the separation posture, restraining the flexible effect of the mechanism, and fulfilling the effectiveness of the design value of multi-rigid body.

    1. Introduction

    The embedded ejection technology of air-to-air missile is to provide the concealed shooting effect for the advanced fighter.Before the missile launching, the carrier first opens the weapon cabin door, and then the launching mechanism ejects the missile with a certain downward velocity and angular velocity, and the missile finally fires away after falling for a certain distance.During the ejection process, the safety of the missile separation is mainly associated with the separation speed and separation angular velocity of the missile, given by the launch mechanism.

    In the last years,some researches on the launch safety of missile have been carried out, and some good results have been achieved,including the influence of the pneumatic flow field [1-4], the ejection dynamics of separation [5-7], and the ground test verification techniques [8], etc. Merrick et al. [1] studied the free-fall trajectory and attitude of the embedded missile under the action of airflow, and pointed out that the missile attitude changed dramatically under the influence of complex airflow, leading to hidden dangers to the carrier safety. By utilizing the overlapping grid technology,Zhang et al.[3]conducted a numerical simulation comparison between the external hanging delivery and the embedded launching,and obtained the falling rules of the external hanging delivery and the embedded launching body under the conditions of subsonic velocity and supersonic velocity. It was pointed out that the initial downward angular velocity of the missile was crucial for the purpose of the launching safety of the embedded missile. Liu et al. [4] analyzed the impact of overload centrifugal force generated by large maneuvering of aircraft on the ejection separation parameters of air-to-air missile based on Lagrange dynamic equation. Yan et al. [6] utilized the improved delayed separation vortex-simulation method combined with the dynamic grid technology to carry out unsteady calculation of the three-dimensional flow field in the process of internal ejection separation, and obtained the influence law of different inlet Mach number, projectile launching velocity, initial angular velocity and other parameters on the trajectory of projectile body.Based on the rigid-flexible coupling multi-body dynamics model, Wang [7]studied the influence of structural parameters of the air-to-air missile embedded catapult launching system on its catapult separation parameters, and pointed out that the initial angular velocity of the missile was “rise” due to the flexibility of the mechanism during the catapult launching.

    In order to ensure the launch safety, the embedded missile needs a downward velocity and downward angular velocity relative to the carrier. However, due to the large ejection force and the strictly lightweight design, the components of the ejection mechanism suffer from serious flexible deformation,resulting in critical deviation of the angular velocity of the missile from the designed value, even leading to the “head-up” problem which threats the safety of the carrier. Currently, there are two common methods to solve the flexible deformation problem for high-speed motion mechanism:

    1) Increase the stiffness of the mechanism components. It can reduce the flexibility deformation effect of the mechanism. But it also significantly increases the weight and the volume of the ejection mechanism, which is conflict with the lightweight design requirement.

    2) Passive control method using damping layer [9,10] and active control method based on piezoelectric actuation [11-13]. The airborne missile catapult launcher requires a stringent, lightweight and compact design with large ejection force and short actuation time (less than 0.15 s). The passive control method of damping layer and the active control method of piezoelectric actuator are difficult to be implemented due to the limited control effect and high engineering cost.

    Variable topology design is a new method [14-17] used for optimizing the mechanism function.Currently,the optimization of the dynamic behavior of the flexible mechanism is rarely reported.This paper proposes a novel variable topology design to optimize the airborne missile ejection separation parameters by using the variable kinematic pair to realize the variable topological design.This method focuses on the suppression of the flexibility effect of the entire ejection mechanism instead of individual component.The variable topology design can also be used for the dynamic performance optimization of other parallel high-speed flexible actuators.

    2. Dynamics model of the missile ejection based on a fixed topological launcher

    The dynamics model of the missile ejection process is based on the combination of the finite element method and the multi-body dynamic Lagrange equation. Considering the complexity of the mechanism model, the mechanical components that flexible processed in finite element software(ANSYS)are assembled according to the actual motion constraint pair.

    Firstly,the modal analysis of each flexible component is carried out via the finite element software(ANSYS),and the modal neutral file is generated, which is required by the multi-body dynamics software (ADAMS) for multi-flexible body dynamics modeling. After the components being assembled according to the actual motion constraint pairs,the ejection power is set up.Finally,the virtual prototype simulation model of the ejection mechanism of airborne missile is established.

    The detailed flexible process of components can be divided into four steps:

    1) Three-dimensional (3D) modeling

    Components with simple shapes can be directly modeled in ANSYS.But for those of complex shapes,the modeling can be built in software Unigraphics NX(UG)or other professional 3D modeling software, and then imported into ANSYS. This paper adopts the second method.

    2) Material attribute input

    Input material properties, including material density, elastic modulus and Poisson's ratio.

    3) Grid generation

    For 3D solid model,the Brick 8 node 185 mesh in ANSYS can be selected,with the mesh side length smaller than 3 mm.For modal calculation,default Block Lanczos method is applied,which is often used in solid grid,with many advantages such as wide application,high precision and fast calculation speed.

    4) Output files to ADAMS

    The file output for ADAMS needs to set the order of the modal,which can be generally taken based on engineering experience.The more orders, the greater the computational burden of the model.Here,the first 15-30 orders of the modal are selected.

    2.1. Mechanism configuration of the fixed topological launcher

    The airborne missile launch mechanism is based on the configuration of the US AMREL mechanism,which has a double yshaped mechanism, with ejection power derived from the hydraulic piston rod which can only move horizontally, as shown in Fig.1.The hydraulic power drives the piston of the horizontal piston rod,pushing the horizontal piston rod(denoted as no.4 in Fig.1)to slide forward horizontally.The front head of the missile is located at the left end of the missile shown in Fig.1.

    Each y-shaped mechanism constitutes two isosceles triangles for converting the horizontal thrust given by the hydraulic piston rod into a vertical downward thrust of the push arm, thereby pushing the missile to move downwards at a high speed. After reaching a certain distance, the missile leaves the launcher at a certain speed and attitude. The y-shaped mechanism at the front end of the launch device is designed to be longer than the y-shaped mechanism at the rear end.The launch mechanism has an assembly height of 109 mm in initial contraction state and an ejection stroke of 350 mm.Ejection stroke represents that the missile is separated from the launcher structure, after 350 mm downward of the mass center of the missile from the initial position.

    2.2. Lagrange dynamic equation for multi-flexible mechanism

    Fig.1. The configuration of the AMELT launcher in its initial stage: 1. Structure of the ejection mechanism, 2. Front support arm, 3. Front push arm, 4. Horizontal hydraulic piston rod, 5. Rear support arm, 6. Hydraulic actuated cylinder, 7. Rear push arm, 8.Missile.

    Fig. 2. Kinematic relations in the flexible body coordinate system.

    The multi-flexible dynamics is described by using the coordinate system as shown in Fig. 2, which consists of the inertial coordinate system(er)and the floating coordinate system(eb).The former does not change with time,and the latter is to describe the movement and deformation of the flexible body,which is changed with the deformation of the flexible body.

    In the coordinate system er, the motion of any pointPon the flexible body can be decomposed into a rigid translation and a rigid rotation relative to the system origin, or a deformation motion relative to the floating coordinate system.The positional motion ofPrelative toGcan be expressed as:

    where r for the vector of pointPin the inertial coordinate system;rBis the vector of the origin of the floating coordinate system in the inertial coordinate system; A for the transformation matrix; spis the vector of pointPin floating coordinate system when the flexible body is not deformed; upfor the deformation vector of pointP,calculated by the modal coordinates:

    where φPfor the assumed deformation modal matrix at pointPthat satisfies the ritz basis vector requirement, q is the generalized coordinate of deformation.

    1) Kinetic energy T

    The position,direction and mode of any nodePbefore and after deformation in the generalized coordinates of the flexible body are defined as:

    The kinetic energyTof the flexible body is:

    The mass matrix M(ξ)is the 3×3 dimension matrix,which can be expressed as:

    wheret,r,andmare the translational,rotational and modal degrees of freedom,respectively.

    2) Potential energy W

    Potential energyWincludes two parts:gravity potential energy and elastic potential energy, which can be expressed by Eq. (8):

    where K is the stiffness matrix of structural components, corresponding to the modal coordinates of q.

    3) Dynamic equation of Multi-flexible bodies

    Lagrange equation for flexible bodies is described by:

    among them: Ψ is constraint equation; Q is generalized force projected in defined generalized coordinates ξ; λ is Lagrange multiplier corresponding to the constraint equation;Lis the Lagrange term, defined asL=T-W. Then the final differential equation is:

    3. Dynamics simulation analysis of ejection process of a fixed topological launch mechanism

    The simulation and analysis are performed on missile speed,the missile angular velocity and the vertical thrust, under the conditions of multi-rigid ejection mechanism and the multi-flexible ejection mechanism. As the missile ejection device is a highprecision mechanism, the assembly clearance of each joint is not considered in this study.

    3.1. Main simulation parameters

    Table 1 shows the main simulation input parameters of the airborne missile ejection dynamics model. The front support arm,the front push arm,the rear support arm and the rear push arm are made of aluminum alloy, with an elastic modulus ofE=70×109and the Poisson's ratio is 0.3.The missile is treated as a rigid body.The material of the piston rod is high-strength alloy steel, with an elastic modulus ofE=210×109, and the Poisson's ratio is 0.3.

    3.2. Simulation results

    The dynamic simulation analysis of the ejection process is carried out, under the assumption of the multi-rigidity ejection mechanism and the flexible state of the mechanism,respectively.In the flexible simulation model,the front support arm,the front push arm, the rear support arm, the rear push arm, and the horizontal piston rod are flexibly processed in finite element software.

    Fig. 3 shows the ejection power generated by the horizontal hydraulic cylinder. The first 10 ms is the ejection power establishment phase, with a peak value of 58367 N at 10 ms. After that, the ejection power gradually decreases. The missile and the launcher are separated at 110 ms, with a power value of 26788 N.

    Table 1 Main simulation parameters.

    Fig. 3. Horizontal hydraulic ejection force.

    Fig. 4 shows the separation speeds during the missile ejection process. It indicates that, under multi-rigid body conditions (black line in Fig. 4), the speed between the missile and the launching device is 7.7 m/s at the time of separation;the blue dotted line is the simulation result of the flexible mechanism, with a separation speed of 7.85 m/s. As the difference between the two results is small, the flexibility of the mechanism has little effect on the separation speed.

    Fig. 4. Separation speed simulation during missile ejection.

    Fig. 5. Angular velocity of the missile during ejection process.

    Fig. 5 shows the angular velocity of the missile during ejection process.Among them,the simulation result of the multi-rigid body demonstrates that the angular velocity at separation time is -23 deg/s,while the value is+3 deg/s at the time of ejection separation for the flexible mechanism.It can be seen that the flexibility of the mechanism has a great influence on the separation angular velocity. For the air-to-air missile, the ejection angular velocity range is about -20 ± 10 deg/s.

    For air-to-air missile ejection separation of angular velocity, to ensure that a missile can rapidly away from the aircraft after separation, missiles should have certain downward angular velocity.Normally, this downward angular velocity generally should be greater than -15 deg/s. Therefore, it seems that a downward angular velocity of-23 deg/s could well meet the requirements of the multi-rigid design. However, Fig. 5 shows that the angular velocity of the flexible mechanism is+3 deg/s,which does not fit the range of -20 ± 10 deg/s. It is already in the “head-up” angular velocity state and may fly towards the carrier, giving a serious threat to launch safety.

    The analysis of the vertical force on the missile induced by the push arms can better explain the reason of “head up” problem of the missile during ejection process.Fig.6 is the downward force on the missile induced by the front and rear push arms,under multirigid and multi-flexible simulation conditions, respectively.

    As shown in Fig. 6, the force from the front push arm is in a continuously increasing state,and its value is 14057 N at the time of ejection separation; while the rear push arm force is in a continuous increase state in the early stage but sharply reduced in the later period, with a force of 14631 N at the time of ejection separation.

    Fig. 6. The downward force on the missile by the front and rear push arms under the multi-rigid and multi-flexible body conditions.

    Simulation results demonstrate that after enabling the mechanism flexibility,the vertical force results of the front and rear push arms on the missile shows significant difference with that of the multi-rigid dynamics model, especially at the end of the ejection,where the force of the front push arm is decrease rapidly instead of continuously increasing. The force is -9732 N at the time of separation; In contrast, for the rear push arm, its value continues increasing instead of decreasing at the end,reached 34072 N at the time of ejection separation.

    Therefore,the huge force changes of the missile by the front and rear push arms could result in seriously deviating from the ideal design value of the angular velocity of the missile at the moment of ejection separation, even the occurrence of “head-up”, which seriously threatens the safety of missile launch. Consequently, to resolve this problem, it is necessary to reduce the influence of thrust force change caused by the front push arm and the rear push arm.Specifically,in the last stage of the ejection process, the force acted by the front push arm should be increased,while the force of the rear push arm needs to be decreased.

    4. Dynamics simulation analysis of ejection process using variable topological design

    4.1. Variable topological design of ejection mechanism

    To prevent the angular velocity of the airborne missile from deviating from the designed value, it is necessary to suppress the flexibility effect of the entire ejection mechanism.Considering that the ejection time is only 0.11 s and the ejection force reaches up to several tons, the existing flexible suppression methods of various mechanisms have great limitations.

    Fig.7 is the variable topology design of ejection mechanism.At the hinge of the shaft hole of the front main arm and the horizontal piston rod, a clearance having a specific direction and size is designed in the rotary hole, which is different from the random machining error generated during manufacturing. The clearance designed in this research is a horizontal U-shaped gap located behind the rotating shaft. In the initial state, the rotating shaft is located at the front of the U-shaped gap.

    In the initial stage of ejection,due to the existence of horizontal U-shaped gap,the piston rod cannot apply force on the front push arm. The ejection power of the piston rod can only act on the missile via the rear push arm. During the ejection process, the rotating shaft in U-shaped gap slips backward, reaching the rear end of the U-shaped gap after the mechanism stretching a certain distance. The horizontal piston rod then can generate horizontal force on the rotating shaft, further pushing the front push arm which acts on the missile in vertical downward direction. The mechanism continues to operate until the missile is separated from aircraft.

    Fig. 7. Variable topology design of ejection mechanism.

    During the ejection process, the mechanism performs a topological transformation when the rotating shaft reaches the rear end of the U-shaped gap. The front pushing arm is changed from the rotating and the sliding to the relative rotation with respect to the piston rod, that is, the movement of the Prismatic-Revolute (PR)pair is changed to a Revolute(R)pair.Therefore,the U-shaped gap topology design here is referred to as the Prismatic-Revolute-Revolute (PR-R) kinematic pair. After the topology of the mechanism is changed, the ejection power is provided by both the front and rear push arms instead of only rear arm.

    4.2. Dynamic contact force model of the PR-R kinematic pair

    In the case of the variable topology design, the rotation pair adopted in the fixed topology mechanism is not fit as the motion constraint of the front push arm,the horizontal piston rod and the rotating shaft.Instead,the contact force between the rotating shaft and the horizontal piston rod is used to judge whether the rotating shaft and the horizontal piston rod have the interaction force.While considering the practical design of the structure, the front push arm and the rotating shaft continue to adopt the rotation constraint.

    The model of normal contact forceFof PR-R motion constraint considers both the material stiffness at the contact location and the damping effect in impact event.In this case,the classical Lankarani-Nikravesh dynamics model [18] is used:

    Among them:Kfor the contact stiffness of the rotating shaft and the horizontal piston rod at the gap,5×107 N/mm is taken in this paper.Cfor the damping coefficient of the contact,δ for embedding depth,nfor the index coefficient, take 2.2 in this paper. The relationship between multiple groups of embedded depth δ and contact forceFcan be calculated according to the structural design parameters at the hinge of the launching mechanism by using the finite element method, and thenKandnin Eq. (11) can be determined by combining results of the least square method and engineering experience.

    whereeis the material recovery coefficient of the horizontal plug rod structure, take 0.8 in this paper; Pc5-Antarctica_JASTP_keywords_V2 is the initial relative velocity when the rotating shaft is in contact with the horizontal piston rod.

    4.3. Simulation results and analysis

    The simulation was conducted on the forces of the front and rear push arm, and the ejection separation angular velocity of the variable topological ejection mechanism after using the PR-R gap pair.The PR-R gap in the simulation model is set to 2 mm.

    Fig.8 is a normal contact force at the PR-R gap pair.It shows that the rotation axis of the PR-R gap pair begins to contact with the piston rod at 6.7 ms, and the contact force is 18700 N. At this moment,the topological form of the ejection mechanism changes.Before that, the piston rod has no force applied on the front push arm. While after that, the piston rod has a force on the front push arm, and the force gradually increases, reaching a peak of 24542 N at 55 ms, and then the contact force gradually decreases and approaches zero at 110 ms.

    Fig. 8. Contact force at the PR-R gap.

    Fig. 9 shows the missile displacement during ejection process.At 6.7 ms, the missile is in the initial position of the ejection compared with the entire ejection stroke of 350 mm, with only 0.4 mm downward displacement. Therefore, the topology switching of the variable topology mechanism completes in the initial startup phase.

    The purpose of the clearance design is to change the acting force on the missile by the front and rear push arm at the last stage of the ejection process, further to avoid the “head up” problem of the missile. Fig. 10 shows the influence of clearance design on the vertical force acted by the front and rear push arms on the missile during the ejection process.

    Fig. 9. Missile displacement during ejection process.

    Fig.10. Effect of clearance design on ejection force of the front and the rear push arm.

    Fig.11. Separation angular velocity simulations for different gap sizes.

    Before 6.7 ms,the vertical force of the variable topology ejection mechanism to the missile is 0; after about 50 ms, the variable topology design is adopted, and the force by the front push arm is higher than that of the fixed topological mechanism,while the rear push arm has a smaller force on the missile than that of the fixed topology mechanism. Therefore,the variable topology mechanism can generate a larger pitch down torque on missile in the late stage of the ejection, preventing the missile from “heading up” event.

    Fig.11 shows the separation angular velocities of missile with different sizes of PR-R gap pairs. It is simulated under the gaps of 0 mm, 0.5 mm,1 mm,1.5 mm, and 2 mm, respectively.

    Results show that the larger the gap, the higher the angular velocities.With the gap size of 2 mm,the missile achieves the same separation angular velocity of -23 deg/s as in the condition of multiple rigid bodies at the time of separation. It can be seen that the PR-R gap pair can greatly adjust the ejection separation angular velocity of the airborne missile ejection mechanism to ensure the safety of the launch.

    5. Discussion

    In order to facilitate readers to reproduce the simulation process and results,the mechanism in this paper is a simplified mechanism of the original mechanism,and the members of the mechanism are given in simple shape, namely rectangular members with holes,without considering the mutual structural interference of the members during the assembly of the mechanism and during the movement process. The fabrication and subsequent verification of the experimental mechanism require cutting and grooving of the interferential parts of the mechanism during assembly and movement. However, the shape of such mechanism components is relatively complex,which is not suitable for readers to simulate and reproduce. Therefore, there is no experimental data of physical prototype which is completely consistent with the virtual prototype state in this paper.However,this design has been verified by a real object prototype and has a good effect of inhibiting the flexible effect of the mechanism. After optimization design, the ejection separation parameters consistent with the rigid hypothesis can be achieved.

    Readers can manufacture a physical prototype for experimental verification based on the design parameters provided in the paper,and the design parameters of the reader's physical prototype does not need to be completely consistent with those in the paper. The relevant design data of ejection mechanism given in this study is sufficient for readers to simulate and reproduce the results.

    6. Conclusion

    Due to the flexibility of the various components of the airborne missile ejection mechanism, the missile ejection separation parameters seriously deviate from the ideal design value of the multirigid body,which threatens the launch safety.The variable topology flexible mechanism design strategy proposed in this paper aims to solve this problem. The main conclusions are:

    1) The flexibility of the ejection mechanism has little effect on the separation speed of the missile, but it can greatly influence the separation angular velocity and threaten the launch safety;

    2) Adopting the variable topology design via clearance design that can effectively adjust the ejection force on the missile produced by the front and rear push arms, thereby realizing a great adjustment of the separation angular velocity of the missile.

    3) Different ejection separation angular velocities can be obtained via different gap sizes.A separation angular velocity of-23 deg/s can be achieved when the gap size is 2 mm.

    4) The dynamic design of the variable topology multi-flexible mechanism is to optimize the dynamic performance of the mechanism from the overall perspective of the flexible mechanism, and is not limited to each individual component of the mechanism.This design is easy to be implemented and effective in engineering,and can be applied to similar high-speed flexible mechanisms.

    Acknowledgement

    I would like to thank one of my funding support, which is Aerospace Science Foundation (NO. 2017ZC12008).

    国产精品嫩草影院av在线观看| 精品久久久久久久人妻蜜臀av| 日韩高清综合在线| 亚洲国产精品合色在线| 久久久国产成人免费| 欧美成人a在线观看| 日本欧美国产在线视频| 欧美色视频一区免费| 亚洲欧美中文字幕日韩二区| 国产黄色视频一区二区在线观看 | 国产精品永久免费网站| av在线老鸭窝| 亚洲三级黄色毛片| 97超碰精品成人国产| 99久久无色码亚洲精品果冻| 青青草视频在线视频观看| 日日干狠狠操夜夜爽| 国产成人aa在线观看| 国产成人a区在线观看| 三级男女做爰猛烈吃奶摸视频| 欧美日本亚洲视频在线播放| 日本爱情动作片www.在线观看| 久久99热6这里只有精品| 少妇的逼水好多| 精品人妻偷拍中文字幕| 一边亲一边摸免费视频| 在线免费观看不下载黄p国产| 99热这里只有是精品50| 亚洲电影在线观看av| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 欧美97在线视频| av.在线天堂| 免费av毛片视频| 免费搜索国产男女视频| 精品久久久久久久久久久久久| 久久久久久久久久久免费av| 国产成人免费观看mmmm| 亚洲精品影视一区二区三区av| 如何舔出高潮| 亚洲经典国产精华液单| 免费不卡的大黄色大毛片视频在线观看 | 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| 啦啦啦韩国在线观看视频| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 乱系列少妇在线播放| 国产探花在线观看一区二区| 内射极品少妇av片p| 亚洲性久久影院| 日本免费a在线| av播播在线观看一区| 国产极品天堂在线| 人人妻人人澡欧美一区二区| 久久精品国产自在天天线| 日韩欧美三级三区| 亚洲最大成人av| 精品久久久久久久末码| 亚洲欧美中文字幕日韩二区| 久久久国产成人精品二区| 亚洲在线自拍视频| 少妇熟女aⅴ在线视频| 成人毛片a级毛片在线播放| 少妇丰满av| 丰满乱子伦码专区| 精品人妻熟女av久视频| 啦啦啦观看免费观看视频高清| 色网站视频免费| 中文欧美无线码| 国产精品三级大全| 91精品伊人久久大香线蕉| 国产真实乱freesex| 啦啦啦啦在线视频资源| 中国国产av一级| 精品欧美国产一区二区三| 精品人妻一区二区三区麻豆| 一级爰片在线观看| 国产免费又黄又爽又色| 国产 一区 欧美 日韩| 97热精品久久久久久| 一夜夜www| 国产在视频线精品| 免费观看的影片在线观看| 成人性生交大片免费视频hd| 99久久无色码亚洲精品果冻| 午夜视频国产福利| 国产精品一区二区在线观看99 | 亚洲av不卡在线观看| 色尼玛亚洲综合影院| 亚洲人成网站在线观看播放| 赤兔流量卡办理| 秋霞伦理黄片| 黄色配什么色好看| 成人毛片60女人毛片免费| 欧美日本亚洲视频在线播放| 天堂√8在线中文| 国产黄色小视频在线观看| 国产精品熟女久久久久浪| 久久草成人影院| 欧美潮喷喷水| 国产高清不卡午夜福利| 人人妻人人澡欧美一区二区| 亚洲内射少妇av| 一区二区三区四区激情视频| 久99久视频精品免费| 久久精品国产亚洲av天美| 亚洲人成网站在线观看播放| 国产午夜福利久久久久久| 亚洲精品亚洲一区二区| 国产av码专区亚洲av| 亚洲美女搞黄在线观看| 国产精品熟女久久久久浪| 三级毛片av免费| 99久久人妻综合| 在线播放无遮挡| 中文亚洲av片在线观看爽| 人人妻人人澡欧美一区二区| 亚洲精品色激情综合| 成人午夜高清在线视频| 最近2019中文字幕mv第一页| 中文在线观看免费www的网站| 两性午夜刺激爽爽歪歪视频在线观看| ponron亚洲| 久久婷婷人人爽人人干人人爱| 日本午夜av视频| 亚洲精品影视一区二区三区av| 日本一本二区三区精品| 亚洲自拍偷在线| 亚洲欧美日韩卡通动漫| 青春草视频在线免费观看| 中文天堂在线官网| 亚州av有码| 91aial.com中文字幕在线观看| 黄片无遮挡物在线观看| 日本与韩国留学比较| 有码 亚洲区| 毛片一级片免费看久久久久| 国产精品久久久久久久久免| 久久精品熟女亚洲av麻豆精品 | 欧美不卡视频在线免费观看| 国产毛片a区久久久久| 亚洲精华国产精华液的使用体验| 午夜日本视频在线| 精品午夜福利在线看| 卡戴珊不雅视频在线播放| 九九爱精品视频在线观看| 精品酒店卫生间| 亚洲精品亚洲一区二区| 日本猛色少妇xxxxx猛交久久| 日本欧美国产在线视频| 国产精品伦人一区二区| 亚洲精华国产精华液的使用体验| av在线观看视频网站免费| 少妇的逼好多水| 插阴视频在线观看视频| 亚洲美女搞黄在线观看| 亚洲国产精品合色在线| 少妇丰满av| 亚洲av电影在线观看一区二区三区 | 一级爰片在线观看| 国产午夜福利久久久久久| 最近最新中文字幕免费大全7| 日韩av在线大香蕉| 建设人人有责人人尽责人人享有的 | 最近2019中文字幕mv第一页| 国产色爽女视频免费观看| 91狼人影院| 高清午夜精品一区二区三区| 色噜噜av男人的天堂激情| 18禁在线无遮挡免费观看视频| av黄色大香蕉| 精品国内亚洲2022精品成人| 亚洲欧美日韩东京热| 身体一侧抽搐| 男女边吃奶边做爰视频| 亚洲av免费高清在线观看| 精品国内亚洲2022精品成人| 亚洲激情五月婷婷啪啪| 七月丁香在线播放| 亚洲五月天丁香| 国产av一区在线观看免费| 国产精品一区二区在线观看99 | 日本黄大片高清| 国产精品久久久久久精品电影| 亚洲av电影在线观看一区二区三区 | 亚洲av福利一区| 麻豆精品久久久久久蜜桃| 国产精品.久久久| 亚洲国产欧洲综合997久久,| a级毛色黄片| 亚洲av二区三区四区| 天堂√8在线中文| 晚上一个人看的免费电影| 国内精品美女久久久久久| 我要搜黄色片| 免费一级毛片在线播放高清视频| 久久久成人免费电影| 国产亚洲精品久久久com| 亚洲av成人精品一区久久| 十八禁国产超污无遮挡网站| 丰满乱子伦码专区| 一边摸一边抽搐一进一小说| 久久人妻av系列| 中文字幕av成人在线电影| 国产成人a∨麻豆精品| 国产激情偷乱视频一区二区| 色吧在线观看| 18禁动态无遮挡网站| 午夜福利在线观看免费完整高清在| 激情 狠狠 欧美| av线在线观看网站| 国产极品天堂在线| 国产真实乱freesex| 大香蕉久久网| 免费大片18禁| 黄色欧美视频在线观看| 亚洲怡红院男人天堂| 韩国av在线不卡| 99视频精品全部免费 在线| 亚洲真实伦在线观看| 亚洲婷婷狠狠爱综合网| 亚洲精华国产精华液的使用体验| АⅤ资源中文在线天堂| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 日韩亚洲欧美综合| 国产av码专区亚洲av| 视频中文字幕在线观看| 少妇的逼水好多| 99久久精品热视频| 国产精品蜜桃在线观看| 日本五十路高清| 国产在视频线精品| 国产精品1区2区在线观看.| 中文字幕熟女人妻在线| 欧美97在线视频| 久久精品久久久久久噜噜老黄 | 最近中文字幕2019免费版| 中文字幕熟女人妻在线| 秋霞伦理黄片| 亚洲熟妇中文字幕五十中出| 精品久久久久久久久久久久久| 久久99蜜桃精品久久| 99久国产av精品国产电影| 伊人久久精品亚洲午夜| 99在线人妻在线中文字幕| 99热这里只有精品一区| 黄片无遮挡物在线观看| 精品无人区乱码1区二区| 欧美日韩一区二区视频在线观看视频在线 | 欧美激情久久久久久爽电影| 亚洲av电影在线观看一区二区三区 | 国语对白做爰xxxⅹ性视频网站| 麻豆成人av视频| 少妇的逼好多水| 国产精品永久免费网站| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 日韩精品青青久久久久久| 最后的刺客免费高清国语| 国产伦一二天堂av在线观看| 精品久久久久久久久av| 一级爰片在线观看| 久久99热这里只有精品18| 国产69精品久久久久777片| 午夜福利在线观看免费完整高清在| 欧美激情在线99| 日韩欧美国产在线观看| 午夜免费男女啪啪视频观看| 青春草亚洲视频在线观看| 亚洲av免费高清在线观看| 国产免费又黄又爽又色| 在线观看66精品国产| 亚洲精品乱久久久久久| 内地一区二区视频在线| 国产av在哪里看| 日韩亚洲欧美综合| 三级男女做爰猛烈吃奶摸视频| 天美传媒精品一区二区| 美女被艹到高潮喷水动态| 久久亚洲精品不卡| 精品人妻熟女av久视频| 日日摸夜夜添夜夜添av毛片| 久久精品国产亚洲av天美| 最近中文字幕2019免费版| 丰满人妻一区二区三区视频av| 乱码一卡2卡4卡精品| 晚上一个人看的免费电影| 亚洲中文字幕一区二区三区有码在线看| 国产黄片视频在线免费观看| 天堂影院成人在线观看| 亚洲精品,欧美精品| 97热精品久久久久久| 国产欧美日韩精品一区二区| av在线观看视频网站免费| 免费搜索国产男女视频| 久久久精品大字幕| 只有这里有精品99| 国产免费视频播放在线视频 | 日本三级黄在线观看| 1000部很黄的大片| 欧美性感艳星| 天堂影院成人在线观看| 久久热精品热| av在线播放精品| 久久久精品大字幕| 国产三级在线视频| 男女视频在线观看网站免费| 成人综合一区亚洲| 欧美日本亚洲视频在线播放| 好男人视频免费观看在线| 大又大粗又爽又黄少妇毛片口| 亚洲真实伦在线观看| 久久午夜福利片| 欧美变态另类bdsm刘玥| www日本黄色视频网| 亚洲va在线va天堂va国产| 亚洲成av人片在线播放无| 久久亚洲精品不卡| 亚洲四区av| 可以在线观看毛片的网站| 久久99蜜桃精品久久| 午夜老司机福利剧场| 夜夜爽夜夜爽视频| av播播在线观看一区| 插阴视频在线观看视频| 在线播放无遮挡| 亚洲综合色惰| 99热全是精品| 欧美最新免费一区二区三区| a级毛片免费高清观看在线播放| 国产v大片淫在线免费观看| 熟女电影av网| 中国国产av一级| 国产在视频线在精品| 又粗又爽又猛毛片免费看| 日韩av在线大香蕉| 国产伦精品一区二区三区四那| av免费观看日本| 日本欧美国产在线视频| 亚洲在久久综合| 黑人高潮一二区| 嘟嘟电影网在线观看| 天堂中文最新版在线下载 | 欧美变态另类bdsm刘玥| 搡女人真爽免费视频火全软件| av女优亚洲男人天堂| 国产成人freesex在线| 欧美精品国产亚洲| 国产v大片淫在线免费观看| 寂寞人妻少妇视频99o| 免费看美女性在线毛片视频| 午夜久久久久精精品| 久久久午夜欧美精品| 日韩欧美精品免费久久| kizo精华| 亚洲精品乱久久久久久| 亚洲精品456在线播放app| 国产成人精品久久久久久| 免费大片18禁| 天堂网av新在线| 免费观看a级毛片全部| 99在线视频只有这里精品首页| www.av在线官网国产| 欧美成人a在线观看| av专区在线播放| 日日撸夜夜添| 女人被狂操c到高潮| 日韩欧美 国产精品| 两个人视频免费观看高清| 欧美变态另类bdsm刘玥| 国产成人a区在线观看| 成人毛片60女人毛片免费| 舔av片在线| 久久精品国产亚洲av天美| ponron亚洲| 亚洲va在线va天堂va国产| 少妇熟女aⅴ在线视频| 女人十人毛片免费观看3o分钟| 国内精品一区二区在线观看| 91精品伊人久久大香线蕉| 精品酒店卫生间| av卡一久久| 亚洲成人av在线免费| 日韩高清综合在线| 国产免费一级a男人的天堂| 成人性生交大片免费视频hd| 国产亚洲5aaaaa淫片| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| 久久精品国产鲁丝片午夜精品| 亚洲综合色惰| 夜夜爽夜夜爽视频| 综合色av麻豆| 亚洲在久久综合| 99国产精品一区二区蜜桃av| 老司机影院毛片| 男人的好看免费观看在线视频| 国产精品久久久久久久电影| 美女cb高潮喷水在线观看| 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 中文天堂在线官网| 麻豆av噜噜一区二区三区| 国产午夜精品久久久久久一区二区三区| 婷婷色av中文字幕| 精品久久久久久电影网 | 最新中文字幕久久久久| 日韩欧美国产在线观看| 一边亲一边摸免费视频| 中国美白少妇内射xxxbb| 一夜夜www| h日本视频在线播放| 国产在视频线精品| 国产乱人视频| 久久精品久久久久久久性| 国产真实伦视频高清在线观看| 午夜福利成人在线免费观看| 国产精品,欧美在线| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| 免费观看精品视频网站| 婷婷六月久久综合丁香| 亚洲四区av| 熟女人妻精品中文字幕| 国产午夜精品一二区理论片| eeuss影院久久| 97热精品久久久久久| 亚洲性久久影院| 精品久久久久久电影网 | 久久99热这里只有精品18| 日本一本二区三区精品| 成人特级av手机在线观看| 日韩成人伦理影院| 日韩成人av中文字幕在线观看| 免费无遮挡裸体视频| 成人无遮挡网站| 在线a可以看的网站| 好男人在线观看高清免费视频| 亚洲精品亚洲一区二区| 欧美成人精品欧美一级黄| 日韩精品有码人妻一区| 久久久a久久爽久久v久久| 色哟哟·www| 伦精品一区二区三区| 69av精品久久久久久| 你懂的网址亚洲精品在线观看 | 3wmmmm亚洲av在线观看| 久久久a久久爽久久v久久| 亚洲国产欧美人成| 成人高潮视频无遮挡免费网站| videossex国产| 大话2 男鬼变身卡| 免费黄色在线免费观看| 黄片wwwwww| 老女人水多毛片| 精品久久国产蜜桃| 伊人久久精品亚洲午夜| 日本一二三区视频观看| 又粗又爽又猛毛片免费看| 国产一级毛片在线| 日韩欧美三级三区| 国产精品不卡视频一区二区| 91久久精品国产一区二区成人| 日本免费在线观看一区| 久久久久性生活片| 99久久精品热视频| a级一级毛片免费在线观看| a级一级毛片免费在线观看| 亚洲自偷自拍三级| 97超碰精品成人国产| 女人十人毛片免费观看3o分钟| 亚洲电影在线观看av| 99热6这里只有精品| 中文欧美无线码| 免费看av在线观看网站| 国产亚洲5aaaaa淫片| 欧美一区二区亚洲| 春色校园在线视频观看| 久久鲁丝午夜福利片| 午夜亚洲福利在线播放| 波多野结衣巨乳人妻| av又黄又爽大尺度在线免费看 | 亚洲在线自拍视频| 久久精品国产自在天天线| 麻豆精品久久久久久蜜桃| 午夜日本视频在线| 亚洲成人久久爱视频| 三级男女做爰猛烈吃奶摸视频| 女人十人毛片免费观看3o分钟| 久久久久久久久久成人| 中文字幕熟女人妻在线| 国产又黄又爽又无遮挡在线| 成人国产麻豆网| 国产精品蜜桃在线观看| 国产亚洲5aaaaa淫片| 韩国av在线不卡| 国产极品天堂在线| 韩国av在线不卡| 三级国产精品欧美在线观看| 午夜福利在线观看吧| 欧美激情国产日韩精品一区| 校园人妻丝袜中文字幕| 小蜜桃在线观看免费完整版高清| 久久99热这里只有精品18| 亚洲精品色激情综合| 日本三级黄在线观看| 国内揄拍国产精品人妻在线| 最近2019中文字幕mv第一页| av在线蜜桃| 亚洲精品国产成人久久av| 亚洲国产精品成人久久小说| 全区人妻精品视频| 好男人视频免费观看在线| 国产精品一区二区在线观看99 | 免费观看a级毛片全部| 内地一区二区视频在线| 欧美一区二区国产精品久久精品| 日本熟妇午夜| 国产淫片久久久久久久久| 亚洲精品成人久久久久久| 午夜精品一区二区三区免费看| 美女xxoo啪啪120秒动态图| 精品国产露脸久久av麻豆 | 久久精品国产亚洲av涩爱| 美女cb高潮喷水在线观看| av播播在线观看一区| 久久久a久久爽久久v久久| 成人漫画全彩无遮挡| 在线观看66精品国产| 在线免费观看不下载黄p国产| 草草在线视频免费看| 赤兔流量卡办理| 18禁在线无遮挡免费观看视频| 日韩精品青青久久久久久| 亚洲精品一区蜜桃| 综合色av麻豆| 熟女电影av网| 超碰97精品在线观看| 五月玫瑰六月丁香| 国内精品一区二区在线观看| 能在线免费观看的黄片| 日本wwww免费看| 欧美变态另类bdsm刘玥| 成人无遮挡网站| 人人妻人人澡欧美一区二区| 日本一本二区三区精品| 最近中文字幕高清免费大全6| 麻豆久久精品国产亚洲av| 亚洲av不卡在线观看| 色哟哟·www| 人人妻人人澡人人爽人人夜夜 | 欧美高清性xxxxhd video| 丰满乱子伦码专区| 人体艺术视频欧美日本| 超碰97精品在线观看| 中文精品一卡2卡3卡4更新| 亚洲性久久影院| 在现免费观看毛片| 久久精品国产自在天天线| 成人午夜精彩视频在线观看| 好男人视频免费观看在线| 亚洲精品久久久久久婷婷小说 | 亚洲av成人精品一二三区| 国产又色又爽无遮挡免| 日韩欧美精品免费久久| 亚洲欧美日韩无卡精品| 国产精品电影一区二区三区| 亚洲五月天丁香| 亚洲精品,欧美精品| 别揉我奶头 嗯啊视频| 亚洲av免费在线观看| 欧美高清成人免费视频www| 91在线精品国自产拍蜜月| 国产亚洲精品久久久com| 色5月婷婷丁香| 欧美激情在线99| 最近中文字幕2019免费版| 国产真实乱freesex| 国产精品福利在线免费观看| 亚洲人与动物交配视频| 午夜日本视频在线| 中文字幕制服av| 精品一区二区免费观看| 国产免费又黄又爽又色| 日韩强制内射视频| 两性午夜刺激爽爽歪歪视频在线观看| 一个人观看的视频www高清免费观看| 日本免费在线观看一区| 非洲黑人性xxxx精品又粗又长| 亚洲成人中文字幕在线播放| 久久精品国产鲁丝片午夜精品| 长腿黑丝高跟| 免费av观看视频| 狠狠狠狠99中文字幕| 亚洲精品aⅴ在线观看| 18禁动态无遮挡网站| 美女黄网站色视频| 国产精品野战在线观看| av又黄又爽大尺度在线免费看 | 51国产日韩欧美| 亚洲精品日韩在线中文字幕| 啦啦啦啦在线视频资源| 看非洲黑人一级黄片| 男女那种视频在线观看| a级毛片免费高清观看在线播放| 精品少妇黑人巨大在线播放 | 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 欧美bdsm另类| 大又大粗又爽又黄少妇毛片口| 我要搜黄色片| 日韩,欧美,国产一区二区三区 | 97超碰精品成人国产|