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

    Design and dynamic analysis of a scissors hoop-rib truss deployable antenna mechanism

    2023-07-04 08:02:24BoHnXingkunLiJinSunYundouXuJintoYoYongshengZho
    Defence Technology 2023年6期

    Bo Hn , Xingkun Li , Jin Sun , Yundou Xu , Jinto Yo , Yongsheng Zho ,*

    a Parallel Robot andMechatronic System Laboratoryof HebeiProvince,Yanshan University,Qinhuangdao,066004,China

    b KeyLaboratory ofAdvanced Forging& StampingTechnologyand Science,Ministryof Education ofChina,Yanshan University, Qinhuangdao, 066004,China

    Keywords:Scissors mechanism Deployable antenna Screw theory Kinematic analysis Dynamic analysis

    ABSTRACT

    1. Introduction

    Artificial satellites officially became a high-end equipment around the middle of the 20 centuries,with the Soviet Union 6 for launch of Sputnik 1. Fast forward to now, there are thousands of artificial satellites up orbiting earth. From the perspective of national defence security,the satellite navigation system provides the country with three independent functional guarantees of navigation,positioning and communication[1—4].As an important part of artificial satellite, deployable antenna mechanism mainly include planar antenna mechanisms, rib antenna mechanisms, ring truss antenna mechanisms,membrane antennas mechanisms,and solid surface antennas mechanisms [5—9].

    In defence technology, in order to make the electromagnetic wave of radar reconnaissance satellite have lower frequency and longer wavelength,and to make the artificial satellite have a larger service range, the aperture of the deployable antenna mechanism needs to be larger. Jia [10] attempted to improve the synthesis method of metamorphic mechanisms, enabling the synthesis of parallel mechanisms. Wang [11] proposed a petal-inspired space deployable-foldable mechanism for space applications combining features common to both the flower blooming process and deployable-foldable mechanisms.Based on screw theory,Yang and Shi [12,13] synthesized rectangular pyramidal overconstrained deployable units and linear foldable over-constrained deployable unit, respectively. Han, Dai and Meng [14—16] synthesized deployable mechanisms for ring truss antennas, respectively. Liu[17] proposed a new configuration that can realize twodimensional (2D) planar deployment and designed a novel largescale 2D deployable planar antenna mechanism. Intended for large mesh antennas, Sun [18] introduced a novel deployable antenna comprising a double-ring deployable truss and a cable net reflector.Li[19]proposed a double-level guyed membrane antenna to improve the stiffness of a large-scale tri-prism deployable mast using a collapsible tubular mast. In Ref. [20], a parametric design approach to the petal-type solid surface deployable reflector was presented. Wu [21] proposed the concept of a single-layer deployable truss structure driven by elastic components, which can be applicable to small satellites,and from its stowed state,the structure is self-deployable to a planar regular hexagon configuration. All the above studies have developed typical deployable antenna mechanisms, however, their stiffness will inevitably reduce when it comes to the formation of large aperture deployable antennas,which will also reduce the accuracy of the reflective cable mesh.

    Whether the deployable antenna mechanism can be successfully deployed determines the success or failure of the artificial satellite launch mission. To ensure that the deployable antenna mechanism will be successfully deployed in orbit,it is necessary to analyze the kinematic characteristics of the antenna mechanism.In Ref. [22], a method based on the equivalent concept of first removing a link and then restoring it was proposed for the DOF analysis of multiloop coupled deployable tetrahedral mechanisms.In Ref. [23], based on the aspects of flexibility and controllability through modularity and simple actuation requirements, a kinematic approach applied to a lightweight linkage structure was presented. Wei [24] studied geometry and kinematic of a planarspherical overconstrained mechanisms via closed-form equations.In Refs.[25,26],the possibility of actively controlling the stiffness of tensegrity plate-like structures throughout the changes in the selfstress level and the support conditions were explored. In another study, a general dynamic method for the deployment analysis of cable networks was proposed [27]. In Ref. [28], a large-scale deployable ring truss was introduced, which is equipped with a complete rope-driven driving method and supplements a cable net system to form a complete space antenna. Shi [29] studied the optimization of the antenna mechanism synthesis and obtained the antenna configuration with optimal stiffness and mass synthesis with a genetic algorithm. Nevertheless, the above antenna mechanisms have more than one DOFs,which will lead to a consequence of excessive using of synchronous joints,it means that the weight of the deployable antenna mechanism will become larger, and the deployment reliability of the mechanism will be reduced as the extreme temperature conditions in space may cause synchronous joints to stick.

    Artificial satellites provide three functional guarantees for navigation, positioning and communication. To achieve the above functions, the dynamic characteristics of antenna mechanisms during the deployment process are also of particular importance.Wu and Zhang [30,31] proposed a dynamic method for constructing an origami dynamic model.In Refs.[32—35],a dynamic models of frame structures with clearance joints were developed.Refs. [36,37] investigated the frequency characteristics of deployable mechanisms under thermal load and in the case where the antenna was extended and locked. Siqueira [38] developed a flexible actuator finite element model which can be applied for the modeling of spatial mechanisms existing in several industrial applications. You [39] studied the deployment dynamic of a space telescope via absolute node coordinate method. Zhang [40] introduced a simplified method, which can reduce the threedimensional (3D) cell structure of the ring truss into a 2D “flattened”one.Zhao[41]presented a novel computational method for the in-situ finite element modeling of inflatable membrane structures. Their method is based on geometrical shape measurements using photogrammetry and can be used for further structural analysis.

    The basic scissors mechanism unit is a planar 5R mechanism,which is symmetrical about the center of the rotational joint and has only one DOF,and it has a wide range of applications in security defence field [42—45]. A generic structure for a planar remote center-of-motion (RCM) mechanism with dual scissors-like mechanisms has been presented in Ref.[46].A centralized-driven flexible continuous robot based on a multiple scissors unit mechanism has been proposed in Ref. [47]. In that study, the kinematic and dynamic of the robot were analyzed, and its workspace and deformation performance were verified experimentally.

    As introduced above,when the aperture of deployable antenna mechanism increases, the structural stiffness will decrease, which will reduce the accuracy of the reflective cable mesh. Besides, the antenna mechanism with more than one DOF will greatly increase the weight due to the excessive using of synchronous joints, and will decrease the deployment reliability of the mechanism as the extreme temperature conditions in space may cause synchronous joints to stick.In order to solve the two problems,in this paper,the position of the intermediate rotational joint of the basic scissors mechanism is changed to form an ULSM,and the original parabolic rib is replaced by the ULSM, all the ULSM s are connected to each other so that the aperture of the antenna can be constructed larger.Moreover, in order to improve the structural stiffness of the deployable antenna mechanism,the upper and lower layers of the antenna mechanism are connected by a 3R planar mechanism,the antenna mechanism forms multiple closed-loop mechanisms, so that the SHRTDAM can be obtained. When the two rods of the 3R planar mechanism are collinear, the antenna mechanism is in its fully deployed state,and the boundary singular configuration of the 3R planar mechanism prevents the antenna mechanism from deploying further.As the mechanism has only one DOF,there is no need for the synchronous joints, which can also ensure the deployment reliability. And the reflective cable mesh can be arranged at the both ends of each scissors rod,so that the SHRTDAM proposed in this paper also has the better reflective accuracy of the cable mesh.

    The remainder of this paper is organized as follows.In Section 2,the SHRTDAM is proposed and its configuration and geometric characteristics are analyzed, and the DOF of the SHRTDAM is calculated in Section 3,the result showed that the mechanism has only one DOF. In Sections 4 and 5, the kinematic and dynamic models of the SHRTDAM are established based on screw theory and Lagrange equation. In Section 6, the kinematic and dynamic characteristics are calculated and simulated through MATLAB and Adams software,a ground experiment prototype of 1.5-m diameter was designed and fabricated and a deployment test is conducted,which demonstrated the mobility and deployment performance of the whole mechanism.Finally,conclusions are presented in Section 7, wherein the present work is summarized.

    2. SHRTDAM configuration and geometric characteristics analysis

    As illustrated in Fig. 1, the SHRTDAM consists of scissors mechanisms, nodes, and 3R planar mechanisms. The latter constitute the outer ring truss of the antenna and each rib is composed of three sets of scissors mechanisms, which connects the inner and outer nodes, each type of scissors has the same length.

    To facilitate the description of the SHRTDAM,a unit mechanism is isolated to explain the characteristics of each component. The scissors mechanisms close to the inner node are called inner scissors rods, the scissors mechanisms close to the outer node are called outer scissors rods, the components connecting the inner and outer scissors rods are called middle scissors rods,and the 3R planar mechanisms on the outer ring truss are called chord rods.As it can be seen in Fig.2,the reflective cable mesh can be arranged at the both ends of each scissors rod. So that reflective cable mesh with different curvatures can be obtained by changing the lengths of the scissors rods.Besides,through changing the number of ULSM and scissors rods, the deployable mechanisms with various apertures,curvatures and stiffness can also be constructed.

    Fig.1. 3D SSRLDAM model.

    Fig. 2. SHRTDAM connection.

    The distance between the rotation axis and the inner node center,as well as that between the rotation axis and the outer node center are set as m. The lengths of the scissor rods at the different levels areL1,L2,L3,andL4,whereL1

    As it can be observed in Fig. 4, the vertical projection of the SHRTDAM is a regular polygon with diagonal lines. The projection of each unit mechanism is an isosceles triangle, and the angle between the two waists of the triangle is α.

    Fig. 4. Vertical projection of the SHRTDAM and its folding process.

    Based on the size of each component and the angular relationship between components, the relationship between α,θ, and φ is as follows:

    3. SHRTDAM DOF analysis

    The SHRTDAM can be decomposed into several triangularshaped units. The adjacent units are connected by an inner node,an outer node, and a rib comprising three sets of scissors mechanisms, all units share two inner nodes. Therefore, the SHRTDAM containsNunit mechanisms, which can be further divided into a closed-loop deployable unit mechanism andN-1 open-loop deployable unit mechanisms (Fig. 5).

    Fig. 5. Splitting unit of SHRTDAM.

    3.1. DOF analysis of the closed-loop deployable unit mechanism

    The SHRTDAM can be regarded as a combination of closed-and open-loop deployable unit mechanisms. The basis for the DOF analysis of the SHRTDAM is the DOF analysis of the closed-loop unit mechanism. To this end, in this subsection, a closed-loop unit mechanism is selected to perform the DOF analysis.

    Taking the centroid of the node P1as the originO,the direction of theX-axis is parallel to the rotation axis of the three central rotational joints of a rib, the direction of theZ-axis is vertical upward, and that of theY-axis is determined by the right-hand rule.On this basis, the Cartesian coordinate system of the closed-loop unit mechanism is established (Fig. 6).

    Fig. 6. Cartesian coordinate system of the closed-loop unit mechanism.

    According to screw and graph theories, the component is represented by a circle,that is,the nodes are denoted by capital letters,and the chord rods are denoted by the letter X plus numbers.Moreover, the rotational joints between the components are represented by lines, and the rotational line vectors are distinguished by different subscript numbers. This way, the screw constraint diagram of the closed-loop unit mechanism can be established and is illustrated in Fig. 7.

    Screw theory is used to describe the motion of joints.Taking as an example the rotational joint 5 which connects scissors X11and X12(Fig.7),the twist of the rotational joint 5 can be determined as follows:

    Similarly, the twists of the other rotational joints can be obtained via screw operations.Since the scissors mechanisms on each rib are the same, the twist on other ribs can be obtained quickly via coordinate transformation. Taking the rotational joints 5 and 16 as examples,it can be found that they have the same position on each rib and the angle between the two ribs is 30?, thus, the rotational axis direction of the rotational joint 16 is

    Similarly,the spatial position and the twist of the rotational joint 16 can be obtained easily, consequently, all twists of the closed-loop unit mechanism can be obtained.

    As shown in Fig. 7, in loop I, the node P1is taken as the fixed platform, the node P2is taken as the moving platform, the following equation can be obtained based on the screw theory:

    The arrows pointing in the counterclockwise direction as positive, as a result, the above equation can be written as

    Similarly, for the seven closed-loops (I—VII) in Fig. 7, the corresponding equations can be established as the screw constraint equations of the closed-loop unit mechanisms.

    where ωidenotes the scalar angular velocity of the rotational jointiin the closed-loop unit mechanism,and0represents a null matrix vector with 6 rows and 1 column.

    As for the two redundant 3R planar mechanisms on the outer ring truss, the two outer chord rods perform the opposite movement around the common rotational joint, and their steering is opposite. The relationship of the angular velocity between the outer chord rods is

    Eqs. (6) and (7) can be written in a matrix form as

    where M represents the screw constraint matrix of the closed-loop unit mechanism, N represents a matrix composed of the scalar angular velocities,and0′represents a null matrix with 42 rows and 1 column.

    The matrix M can be expressed by 28 column vectors as

    Therefore,all 28 columns of the matrix M can be expressed as

    The matrix N can be expressed as

    The screw constraint matrix M contains 42 rows and 28 columns. The DOFs of the closed-loop unit mechanism can be determined by solving the zero-space dimension of the screw constraint matrix M. By introducing all twists into the MATLAB software, the rank of the matrix M can be calculated as

    The number of columns of the matrix M is 28,thus,the DOFs of the closed-loop unit mechanism are 28—27=1,i.e.,the closed-loop unit mechanism has only one DOF of retractable motion.

    3.2. DOF analysis of the SHRTDAM

    According to above analysis,the SHRTDAM can be divided into a closed-loop unit mechanism and multiple open-loop unit mechanisms.When a closed-loop unit mechanism and an open-loop unit mechanism are combined together, a double-closed-loop mechanism is formed. The projection diagram of the double-closed-loop mechanism on the horizontal plane is two isosceles triangles with a bevel. The combined unit mechanism and its coordinate system are depicted in Fig. 8.

    Fig. 8. Double-closed-loop unit mechanism: (a) Closed-loop unit mechanism; (b) Open-loop unit mechanism.

    The mechanism on the left side of the combined unit mechanism is the closed-loop unit mechanism exhibited in Fig. 8, while the open-loop unit mechanism is on the right side.The coordinate system is the same as that established by the closed-loop unit mechanism, which is also the global coordinate system of the SHRTDAM.Similarly,based on screw and graph theories,the screw constraint diagram of the open-loop unit mechanism can be obtained (Fig. 9). The components where the open-loop unit mechanism overlaps with the closed-loop one and the corresponding joints are represented by the dotted circle. It can be observed that the shared components are two outer nodes and a rib.

    Fig. 9. Screw constraint diagram of the open-loop unit mechanism.

    According to Fig.9,there are seven closed loops(I—VII)in total,thus,seven screw constraint equations can be obtained.The screw constraint equations of the open-loop unit mechanism are

    where ωidenotes the scalar angular velocity of the rotational jointiin the open-loop unit mechanism, and0represents a null matrix vector with 6 rows and 1 column.

    Since there are also two redundant 3R planar mechanisms in the open-loop unit mechanism, i.e., similar to Eq. (7), the relationship of the angular velocity between the two outer chord rods is as follows:

    Eqs. (13) and (14) can be written in a matrix form as

    where M′represents the screw constraint matrix of the open-loop unit mechanism, N′represents a matrix composed of the scalar angular velocities,and0′′represents a null matrix with 42 rows and 1 column.

    The matrix M′can be expressed by 28 column vectors as

    Thus, all the 28 columns of the matrix M′can be expressed as

    The matrix N′can be expressed as

    By combining Eqs. (6) and (13), the screw constraint equations of the double-closed-loop unit mechanism can be obtained as follows:

    Eq. (19) can be written in a matrix form as

    where

    In Fig. 8, the rib connected by the nodes P1, P2, C, and D and its 11 rotational joints are shared by two unit mechanisms.In Eqs.(6)and(13),the scissors mechanism and the rotational joints shared by the two unit mechanisms are calculated repeatedly. Therefore, when calculating the DOFs of the double-closed-loop unit mechanism,the twists of the repeated calculations need to be removed,that is,$12to $22should be removed once. Consequently, the screw constraint matrix of the double-closed-loop unit mechanism can be expressed as

    The matrix P can also be expressed by column vectors as

    The column vectors in Eq. (23) can be expressed as follows:

    where042represents a null matrix with 42 rows and 1 column.

    The rank of the diagonal matrix P is equal to the sum of the ranks of each diagonal matrix(M and M′′),the rank of the matrix P can be obtained as

    By introducing all twists into the MATLAB software,the rank of the matrices M and M′′can be calculated as

    Thus, the DOFs of the double-closed-loop unit mechanism can be determined as

    whereFdenotes the number of DOFs of the double-closed-loop unit mechanism and t denotes the number of columns of the screw constraint matrix P.

    It is clear that the matrix M′′is a column full rank matrix,Eq.(27)can be written as follows:

    wheretMdenotes the number of columns of the screw constraint matrix M.

    According to Eq. (28), the DOF of the double-closed-loop unit mechanism still depends on the matrix M, i.e., the DOF of the double-closed-loop unit mechanism is the same as that of the closed-loop unit mechanism.

    Based on the previous analysis,it can be understood that,when the open-loop unit mechanism is added on the basis of the combined mechanism(Fig.8),the overall DOF is still the same as that of the closed-loop unit mechanism. Hence, the DOF of SHRTDAM is the same as that of the closed-loop unit mechanism,i.e.,it has only one DOF.

    4. SHRTDAM kinematic analysis

    4.1. SHRTDAM velocity analysis

    The SHRTDAM is a spatially symmetric mechanism,the closedloop mechanism unit can be isolated to perform velocity analysis of the whole SHRTDAM. According to the DOF analysis in Section 3,the closed-loop unit mechanism has only one DOF.Therefore,when an input such as ω1is added to the closed-loop unit mechanism,each angular velocity can be solved via Eq. (6). Based on the geometric parameters of the closed-loop unit mechanism and the screw constraint diagram,the six-dimensional(6D)screw velocity of each node and rod in the closed-loop unit mechanism in Fig. 7 can be obtained through screw algebraic operations.

    In this paper,the inner node P1is set as the fixed platform,and its motion remains static in the 3D space. The screw constraint diagram of the closed-loop unit mechanism is exhibited in Fig. 7,and the screw velocity of each component in Loop I can be expressed as

    where0is a null vector with 6 rows and 1 column that represents the 6D screw velocity of each component.

    In Loop II, the screw velocities of the X11and X12components have been obtained via Eq. (29). Using the obtained screw velocities, the screw velocities of each component in Loop II can be derived as follows:

    The screw velocity of each component in Loops III-VII can be obtained in a similar manner. According to screw velocity, the angular velocity of each component is

    where ω(Vi) represents the first three elements of the screw velocity.

    The centroid linear velocity of each component can be expressed as

    where v(Vi)is a vector that represents the last three elements of the screw velocity and riis a vector from the coordinate origin to the centroid position of componenti.

    Since the entire SHRTDAM is a fully centrosymmetric mechanism, the components in each closed-loop unit mechanism have exactly the same size. Moreover, the center of the inner node is centrosymmetric,if a self-coordinate system is established in each unit mechanism, the components with the same position in each unit mechanism will have the same velocity in the self-coordinate system.

    As it can be seen in Fig. 8, the coordinate systemO-X1Y1Z1is established at the same position of the right open-loop unit mechanism referring to the left closed-loop unit mechanism and the coordinate systemO-XYZ. Since the origins of the two coordinate systems are the same, the velocity of the node C in the coordinate systemO-XYZand that of the nodeEin the coordinate systemO-X1Y1Z1are the same.

    The corresponding relationship regarding the velocity of the other components in Fig. 9 is similar to that between the nodes C and E, in this way, the velocity relationship can be transmitted to the entire SHRTDAM.

    Assuming that the entire SHRTDAM has N closed-loop unit mechanismsintotal,thenNcoordinatesystemsO-XYZ~O-XN-1YN-1ZN-1can be established.O-XYZis set as the global coordinate system, the remaining coordinate systems are relative coordinate systems.As shown in Fig.10,the rotation angle between adjacent coordinate systems around theZ-axis is α.

    In the global coordinate systemO-XYZ, the velocity of each component can be expressed as follows:

    whereirepresents the number of components andjrepresents the number of coordinate systems.

    The rotation transformation matrix in Eq. (34) is defined as

    Through the above analysis,the angular and linear velocities at the centroid of each component in the entire SHRTDAM can be derived and expressed in the global coordinate system through coordinate transformation.

    4.2. SHRTDAM acceleration analysis

    The screw acceleration can be determined by deriving the generalized velocity of each component based on the screw derivative. The screw acceleration equation is as follows:

    where Aidenotes the screw acceleration of componentiin the SHRTDAM,and represents the 6D acceleration column vector of the coincidence point with the origin of the componentireference coordinate system, εirepresents the angular acceleration of the coincidence point with the origin of the componentireference coordinate system,and a is the linear acceleration at the centroid of componenti.

    According to Eq.(36),in the 6D screw acceleration of componenti, the first three elements are its angular velocity, while the last three elements are its linear acceleration minus its centripetal acceleration.

    According to the screw acceleration, the acceleration of the inner node P2in Loop I of the screw constraint diagram in Fig.7 can be derived via theLieoperation [ ] as follows:

    whereP1εP2denotes the angular acceleration of component P2relative to the component P1.The expressions of the two Lie[]are as follows:

    By combining Eqs. (37) and (38), the following equation can be obtained:

    where 0 is a null matrix with 6 rows and 1 column.

    The acceleration screw constraint equations of Loops II—VII in Fig.7 can be obtained in a similar manner.Given the input angular acceleration ε1, the screw acceleration of each component can be calculated using the acceleration screw constraint equations.

    When the screw acceleration of each component has been calculated, the angular and centroid linear accelerations of each component can be obtained based on Eq. (36) as follows:

    where ε(εi$i)denotes the origin of the screw acceleration,which is represented by the first three elements, and a(εi$i) denotes the dual part of the screw acceleration,which is represented by the last three elements.

    Similar to the velocity analysis, all components that have the same position in each unit mechanism have the same acceleration in their own coordinate system. This can be represented in the global coordinate system as follows:

    Through the above analysis, the angular and centroid line accelerations of each component in the SHRTDAM can be calculated in the global coordinate system.

    5. SHRTDAM dynamic analysis

    5.1. Lagrange equation of the SHRTDAM

    The SHRTDAM has one DOF,the mechanism has only one set of independent generalized coordinates. To simplify the calculations,the angle between chord rods is set to 2φ,and that between scissors at all levels is set to 2θ. The input angle φ is selected as the generalized coordinate, and the dynamic equation can be established according to Lagrange equation as follows:

    whereEV,EP, andQare the kinetic energy, potential energy, and generalized force of the SHRTDAM, respectively.

    In dynamic analysis, the clearance and friction of the moving joints are ignored,Qis the generalized driving force applied to the SHRTDAM.

    5.2. SHRTDAM kinetic energy

    5.2.1. Kinetic energy of the nodes

    In this paper, the upper inner node P1at the central position is fixed to facilitate the dynamic analysis. According to the above analysis, the rest of the nodes perform the translational motion.Assuming that the mass of the inner nodes ismH1and that of the outer nodes ismH2, the kinetic energy of the nodes can be expressed as follows:

    5.2.2. Kinetic energy of the scissors mechanisms

    In the SHRTDAM, there are 12 inner scissors rods connected to the fixed node, and their kinetic energy can be expressed as

    whereJ1represents the rotational inertia of the inner scissors rod connected to the fixed node andm1represents the mass of the inner scissors rod.

    The remaining scissors rods, which include 12 inner scissors rods, 24 middle scissors rods, and 24 outer scissors rods, perform plane motion.The kinetic energy expressions are given below.

    The kinetic energy of 12 inner scissors rods can be calculated as

    whereJ2represents the rotational inertia of the inner scissors rod andΟν1denotes the centroid velocity of the inner scissors rod in the global coordinate system.

    The kinetic energy of the 24 middle scissors rods can be calculated as

    whereJ3=m2(L2+L3)2/12 represents the rotational inertia of the middle scissors rod andΟν2denotes the centroid velocity of the middle scissors rod in the global coordinate system.

    The kinetic energy of the 24 outer scissors rods can be calculated as

    whereJ4=m3(L3+L4)2/12 represents the rotational inertia of the outer scissors rod andΟν3denotes the centroid velocity of the outer scissors rod in the global coordinate system.

    Therefore,the kinetic energy of the remaining scissors rods is

    5.2.3. Kinetic energy of the chord rods

    In the SHRTDAM,there are in total 48 chord rods,which perform the composite motion. Their kinetic energy can be calculated as follows:

    5.3. Potential energy and generalized driving forces of the SHRTDAM

    The Potential energy and generalized driving forces are both set as 0 in this paper.

    The velocity of each component can be determined based on the analysis in Section 4 and the dynamic equation of the SHRTDAM can be obtained by substituting Eqs. (43)—(50) into Eq. (42). Thus far, the dynamic model of the SHRTDAM has been established.

    6. SHRTDAM dynamic simulation

    6.1. SHRTDAM simulation model analysis

    A novel deployable antenna mechanism named SHRTDAM, has been designed via the variant scissors mechanism, and its kinematic and dynamic have been analyzed. Due to the error between the theoretical analysis and the feasibility condition of the antenna mechanism,the correctness of the theoretical analysis needs to be verified. The SHRTDAM proposed in this paper is composed of too many components. Thus, the analysis and calculation of the dynamic model would be very complex. As a result, computer simulation is an inevitable choice for analyzing the SHRTDAM.

    For the established model, it is necessary to verify that the simulation model is correct.The Adams/View module provides the function of model verification, which can check the constraints,DOFs, and drive of the SHRTDAM. After the constraints of the imported SHRTDAM model have been created in Adams/View, it is verified that the DOFs of the antenna mechanism are one.

    Another method that can be used to verify the DOFs of the SHRTDAM is motion simulation.In Section 3,it was found the DOFs of the mechanism are one, the mechanism needs only one drive.

    The driving angle function is added to the chord rod,as well as the simulation time and timestep size are set to simulate the deployment process of the mechanism. The simulation model of the mechanism is illustrated in Fig.11.

    Fig.11. SHRTDAM simulation model.

    After adding one drive, the SHRTDAM can move up throughAdams/motionsimulation and complete the deployment process from the folded configuration to the deployed one. This process is depicted in Fig.12.

    Fig.12. Simulation of the deployment process in the Adams software: (a) Folded state; (b) Deploying state; (c) Deployed state.

    The motion of the antenna mechanism is analyzed. When only one drive is added,the motion form of the antenna mechanism can be determined, indicating that the antenna mechanism has only one DOF. This proves that the previous theoretical analysis on the DOFs of the mechanism is completely correct.

    To further verify that the SHRTDAM has one DOF, a ground experiment prototype of 1.5-m diameter was fabricated and a deployment test is conducted. As shown in Fig.13, the main components include five parts:first-level scissors-like linkages,secondlevel scissors-like linkages, three types of different nodes, and 3R linkages.Considering the inevitable errors in the assembly process,a ball screw pair was selected to drive the mechanism. The deployment process is illustrated in Fig. 14. The deployment test proved that the mechanism has only one DOF and can be deployed smoothly and continuously.

    Fig.13. Prototype components.

    Fig.14. 1.5-m prototype: (a) Folded state; (b) Deploying state; (c) Deployed state.

    6.2. SHRTDAM dynamic simulation analysis

    In the previous sections,kinematic and dynamic analyses of the theoretical model of the SHRTDAM were conducted, and the theoretical values of each motion parameter were obtained. In order to verify the correctness and feasibility of the theoretical model,the dynamic simulation software Adams and MATLAB were used.Subsequently,the consistency between the results obtained by the two software was assessed.

    There are several components in the SHRTDAM, and each type of component is located at the same position in each closed-loop unit mechanism. To facilitate the description, in this section, the closed-loop unit mechanism whose own coordinate system coincides with the global coordinate system is selected for the analysis.Taking the closed-loop unit mechanism as an example(Fig.8),the global coordinate system is established at the centroid of the node P1,and the inner scissors rod X13,the middle scissors rod X24,the outer scissors rod X33,the chord rod T1,the inner node P2,and the outer node C are selected for simulation analysis and verification of the theoretical model. Theoretical value calculated by MATLAB software are represented by lines of different colors and simulation value of Adams software are represented by symbols of different shapes.

    Their displacement, angular velocity, centroid linear velocity,angular acceleration,and linear acceleration of the theoretical and simulation value are presented in Fig.15.

    Fig.15. Displacement, velocity, and acceleration of the different components: (a) Displacement of the different components; (b) Centroid linear velocity of the different components; (c) Angular velocity of the different components; (d) Centroid linear acceleration of the different components; (e) Angular acceleration of the different components.

    The curve of the input driving torque of the theoretical and simulation value was obtained as well, and is presented in Fig.16.

    Fig.16. Input driving torque of the SHRTDAM.

    By analyzing the motion curves of the inner node P2and the outer node C in Fig.15,it can be found that the angular velocity and acceleration of the inner node P2and the outer node C are both 0.This indicates that, when the SHRTDAM expands from its folded configuration to the deployed one,the spatial motion form of each inner node P2and outer node C in the entire process is translational motion. By analyzing the motion characteristic curves of each component in Fig.15, when only one drive is applied, the motion curve of each component is certain, indicating that the SHRTDAM has only one DOF. The theoretical results in Figs.15 and 16 were consistent with the simulation results,indicating that the relevant motion theory derived above is correct.

    7. Conclusions

    A new type of deployable antenna mechanism with high reliability and simple structure was proposed in this study, and the kinematic and dynamic characteristics of this mechanism were analyzed, a ground experiment prototype was fabricated and the deployment test was conducted. The experimental results demonstrated the mobility and deployment performance of the proposed mechanism. The highlights of this work are as follows:

    (1) An ULSM was obtained by changing the position of the central rotational joint of traditional scissors mechanism,based on ULSM, the SHRTDAM with a high reliability and high stiffness is proposed, which can be used as large aperture satellite antenna mechanism. Moreover, a ground experiment prototype of 1.5-m diameter was designed and fabricated and a deployment test is conducted, the deployment test result showed that the mechanism has only one DOF and it can be deployed smoothly and continuously.

    (2) The motion relationship between the antenna mechanism components was represented as the screw constraint diagram. Subsequently, the screw constraint equations were extracted from the screw constraint diagram, and the number of DOFs of the antenna mechanism was determined by analyzing and solving them. This method relies on matrix operations, which can be applied to solve problems concerning complex mechanism.

    (3) The kinematic analysis of SHRTDAM was conducted via the screw constraint diagram, the screw and Lie algebraic operations were conducted based on the position of each component in the screw constraint diagram, and through rotation transformation between the coordinate systems,the velocity and acceleration expressions of each component relative to the global coordinate system were deduced.

    (4) Lagrange equation was adopt to establish the dynamic model of SHRTDAM, the torque required to drive the antenna mechanism was simulated and verified by the Adams and MATLAB software,the results verified the correctness of the dynamic modeling and simulation analysis.

    In summary, the mechanism proposed in this paper had the characteristic of high reliability and single-DOF,which can be used as a large aperture antenna mechanism,it could also provide a good reference for the design and analysis of large aperture space deployable antennas.

    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.

    Acknowledgements

    This work was supported by the National Natural Science Foundation of China (Grant Nos. 52105035 and 52075467), the Natural Science Foundation of Hebei Province of China (Grant No.E2021203109), the State Key Laboratory of Robotics and Systems(HIT) (Grant No. SKLRS-2021-KF-15), and the Industrial Robot Control and Reliability Technology Innovation Center of Hebei Province (Grant No.JXKF2105).

    日本91视频免费播放| 91在线精品国自产拍蜜月| 亚洲av.av天堂| 国产一区有黄有色的免费视频| 亚洲情色 制服丝袜| 国产精品秋霞免费鲁丝片| 亚洲国产精品专区欧美| 亚洲激情五月婷婷啪啪| 80岁老熟妇乱子伦牲交| 久久国产精品男人的天堂亚洲 | 99久久精品一区二区三区| 久久久国产欧美日韩av| 乱系列少妇在线播放| 久久久久网色| 亚洲自偷自拍三级| 色吧在线观看| 黑丝袜美女国产一区| 久久精品国产亚洲av涩爱| 18+在线观看网站| 国产在视频线精品| 午夜福利影视在线免费观看| 免费久久久久久久精品成人欧美视频 | 777米奇影视久久| 国产 一区精品| 在线观看免费日韩欧美大片 | 蜜臀久久99精品久久宅男| 一个人看视频在线观看www免费| 80岁老熟妇乱子伦牲交| 丝袜喷水一区| 天堂俺去俺来也www色官网| 天堂俺去俺来也www色官网| 亚洲精品国产成人久久av| 久久人妻熟女aⅴ| 国产成人精品无人区| 尾随美女入室| 亚洲中文av在线| 一本一本综合久久| 国产高清有码在线观看视频| 亚洲精品乱久久久久久| 精品国产一区二区三区久久久樱花| 国产高清不卡午夜福利| 精品人妻一区二区三区麻豆| 亚洲国产精品成人久久小说| 日韩亚洲欧美综合| 人妻少妇偷人精品九色| 欧美丝袜亚洲另类| av卡一久久| 亚洲国产色片| 亚洲国产最新在线播放| 曰老女人黄片| 99热全是精品| 亚洲成人手机| 日本欧美国产在线视频| 久久韩国三级中文字幕| 日本欧美国产在线视频| 国模一区二区三区四区视频| 久久久久网色| 王馨瑶露胸无遮挡在线观看| 精品少妇久久久久久888优播| 男男h啪啪无遮挡| 色哟哟·www| 国产成人a∨麻豆精品| 高清欧美精品videossex| 亚洲欧美中文字幕日韩二区| 国产老妇伦熟女老妇高清| 久久久久国产精品人妻一区二区| 一级毛片电影观看| 伦理电影免费视频| 蜜臀久久99精品久久宅男| 桃花免费在线播放| 高清黄色对白视频在线免费看 | 久久 成人 亚洲| 国产伦精品一区二区三区视频9| 最黄视频免费看| 成人国产麻豆网| 国国产精品蜜臀av免费| 哪个播放器可以免费观看大片| 夜夜看夜夜爽夜夜摸| 人体艺术视频欧美日本| 国产精品不卡视频一区二区| 亚洲电影在线观看av| 欧美xxⅹ黑人| 伊人亚洲综合成人网| 国产av码专区亚洲av| a级毛片在线看网站| av有码第一页| av.在线天堂| 久久99一区二区三区| 亚洲av中文av极速乱| 99国产精品免费福利视频| 99国产精品免费福利视频| 涩涩av久久男人的天堂| 九九爱精品视频在线观看| 久久久a久久爽久久v久久| 国产淫语在线视频| 亚洲国产精品999| 国产日韩欧美在线精品| a 毛片基地| 热re99久久精品国产66热6| 高清不卡的av网站| 观看美女的网站| 一级二级三级毛片免费看| 王馨瑶露胸无遮挡在线观看| 99久久精品国产国产毛片| 免费在线观看成人毛片| 91精品国产国语对白视频| 大片免费播放器 马上看| 午夜福利影视在线免费观看| 欧美日韩国产mv在线观看视频| 纯流量卡能插随身wifi吗| 国产一区有黄有色的免费视频| 三上悠亚av全集在线观看 | 国内少妇人妻偷人精品xxx网站| 亚洲成人av在线免费| 99视频精品全部免费 在线| 大又大粗又爽又黄少妇毛片口| 精品酒店卫生间| 赤兔流量卡办理| 人人妻人人爽人人添夜夜欢视频 | 国产高清不卡午夜福利| 久久热精品热| 亚洲一区二区三区欧美精品| 亚洲伊人久久精品综合| 日本免费在线观看一区| 老司机影院毛片| 久久精品国产a三级三级三级| 欧美97在线视频| 亚洲欧美一区二区三区国产| 蜜桃久久精品国产亚洲av| 乱码一卡2卡4卡精品| 日韩一区二区视频免费看| 亚洲第一av免费看| 久久精品国产亚洲av涩爱| 久久久久久久久久久久大奶| 色94色欧美一区二区| 看免费成人av毛片| 美女脱内裤让男人舔精品视频| 九九在线视频观看精品| 天堂中文最新版在线下载| 美女主播在线视频| 26uuu在线亚洲综合色| 成人黄色视频免费在线看| 又粗又硬又长又爽又黄的视频| 国产一区二区三区av在线| 国产精品秋霞免费鲁丝片| 一级毛片电影观看| 一级毛片电影观看| xxx大片免费视频| 久久 成人 亚洲| 人妻人人澡人人爽人人| 成年人午夜在线观看视频| 国产一级毛片在线| 极品教师在线视频| 熟女av电影| 又爽又黄a免费视频| 精品人妻一区二区三区麻豆| 韩国av在线不卡| 在线观看免费视频网站a站| 伦精品一区二区三区| 亚洲丝袜综合中文字幕| a级片在线免费高清观看视频| 日韩一区二区视频免费看| 熟女电影av网| 三级经典国产精品| 国产极品天堂在线| 在线精品无人区一区二区三| 人妻 亚洲 视频| 久久人妻熟女aⅴ| 一区在线观看完整版| 搡女人真爽免费视频火全软件| 精品久久久久久久久av| 亚洲精品乱码久久久久久按摩| 少妇被粗大猛烈的视频| 观看av在线不卡| 我要看日韩黄色一级片| 亚洲国产精品一区三区| 精品少妇久久久久久888优播| 亚洲av.av天堂| 午夜福利视频精品| 亚洲,一卡二卡三卡| 天堂俺去俺来也www色官网| 亚洲欧美清纯卡通| 内地一区二区视频在线| 赤兔流量卡办理| 伦理电影大哥的女人| 男人舔奶头视频| 国产成人精品无人区| 国产精品国产三级国产av玫瑰| 欧美国产精品一级二级三级 | h视频一区二区三区| 亚洲,欧美,日韩| 成人美女网站在线观看视频| 看十八女毛片水多多多| 人人妻人人爽人人添夜夜欢视频 | 成人黄色视频免费在线看| 少妇的逼好多水| 欧美日韩视频高清一区二区三区二| av国产精品久久久久影院| 精品国产国语对白av| 欧美日韩av久久| √禁漫天堂资源中文www| 日韩欧美精品免费久久| 欧美日韩国产mv在线观看视频| 人妻系列 视频| 美女cb高潮喷水在线观看| 欧美日韩亚洲高清精品| 免费看日本二区| 久久久久人妻精品一区果冻| 久久久久精品性色| 曰老女人黄片| 最近2019中文字幕mv第一页| 一级av片app| 亚洲av成人精品一区久久| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品999| 亚洲精品亚洲一区二区| 国产熟女午夜一区二区三区 | 亚洲国产精品一区二区三区在线| av在线观看视频网站免费| 寂寞人妻少妇视频99o| 亚洲国产欧美在线一区| 亚洲国产精品一区三区| 久久久精品94久久精品| 人妻少妇偷人精品九色| 丝袜在线中文字幕| kizo精华| 又黄又爽又刺激的免费视频.| 欧美xxⅹ黑人| av网站免费在线观看视频| 久久精品国产自在天天线| 高清不卡的av网站| videossex国产| 亚洲va在线va天堂va国产| 人人澡人人妻人| 免费大片18禁| 最后的刺客免费高清国语| 国产精品一二三区在线看| 亚洲国产精品成人久久小说| 久久精品熟女亚洲av麻豆精品| 成人午夜精彩视频在线观看| 精品熟女少妇av免费看| 日韩三级伦理在线观看| 久久综合国产亚洲精品| 免费看日本二区| 如何舔出高潮| 99热这里只有是精品50| av.在线天堂| 久久久久久久久久久久大奶| 亚洲内射少妇av| 美女福利国产在线| 久久精品国产自在天天线| 视频中文字幕在线观看| 十分钟在线观看高清视频www | 超碰97精品在线观看| 亚洲精品456在线播放app| 免费看不卡的av| 成人亚洲欧美一区二区av| 少妇人妻一区二区三区视频| 亚洲精品国产色婷婷电影| 六月丁香七月| 色视频www国产| 亚洲av成人精品一区久久| 永久免费av网站大全| 国产一区二区三区av在线| 美女内射精品一级片tv| 午夜老司机福利剧场| 久久久a久久爽久久v久久| 九九在线视频观看精品| 九九久久精品国产亚洲av麻豆| 国产高清国产精品国产三级| a级毛片免费高清观看在线播放| 少妇人妻精品综合一区二区| 国产成人午夜福利电影在线观看| 伊人亚洲综合成人网| xxx大片免费视频| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 在线亚洲精品国产二区图片欧美 | av有码第一页| 国产欧美日韩综合在线一区二区 | 精品少妇久久久久久888优播| 亚洲av不卡在线观看| 三级经典国产精品| 久久亚洲国产成人精品v| 日本午夜av视频| 欧美国产精品一级二级三级 | 一区在线观看完整版| 伊人亚洲综合成人网| 成人影院久久| 青春草亚洲视频在线观看| 国产精品99久久99久久久不卡 | av播播在线观看一区| 在线观看美女被高潮喷水网站| 18禁在线播放成人免费| 老司机亚洲免费影院| 麻豆成人午夜福利视频| 亚洲精华国产精华液的使用体验| 国产永久视频网站| 亚洲精品日韩av片在线观看| 国产亚洲最大av| 国产 精品1| 日日摸夜夜添夜夜添av毛片| 又黄又爽又刺激的免费视频.| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 搡老乐熟女国产| 内地一区二区视频在线| 九草在线视频观看| 中文欧美无线码| 91精品国产国语对白视频| 久久久a久久爽久久v久久| 这个男人来自地球电影免费观看 | 中文欧美无线码| av国产久精品久网站免费入址| 国产日韩欧美视频二区| 国产精品一区二区三区四区免费观看| 交换朋友夫妻互换小说| 精品一区在线观看国产| 国产在视频线精品| 亚洲色图综合在线观看| 丁香六月天网| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 大码成人一级视频| 亚洲欧洲精品一区二区精品久久久 | 97在线视频观看| 久久综合国产亚洲精品| 精品99又大又爽又粗少妇毛片| 国产一区二区三区综合在线观看 | 涩涩av久久男人的天堂| av在线老鸭窝| 亚洲一区二区三区欧美精品| 日本vs欧美在线观看视频 | h日本视频在线播放| 日本色播在线视频| 亚洲精品乱久久久久久| 女性被躁到高潮视频| 欧美日韩国产mv在线观看视频| 丝瓜视频免费看黄片| 交换朋友夫妻互换小说| 男女啪啪激烈高潮av片| 9色porny在线观看| 高清在线视频一区二区三区| 久久久久国产网址| 日本vs欧美在线观看视频 | 亚洲经典国产精华液单| 国产精品国产三级国产av玫瑰| 中文乱码字字幕精品一区二区三区| 国产日韩欧美在线精品| 两个人的视频大全免费| 中文资源天堂在线| 在线观看一区二区三区激情| 国产熟女欧美一区二区| 国产黄频视频在线观看| 一二三四中文在线观看免费高清| 亚洲欧美日韩另类电影网站| 精品酒店卫生间| 亚洲经典国产精华液单| 日本vs欧美在线观看视频 | 啦啦啦啦在线视频资源| 精品国产露脸久久av麻豆| 草草在线视频免费看| 一级黄片播放器| 赤兔流量卡办理| 亚洲综合色惰| 国产黄色视频一区二区在线观看| 成人免费观看视频高清| √禁漫天堂资源中文www| 熟妇人妻不卡中文字幕| 日韩av在线免费看完整版不卡| 一二三四中文在线观看免费高清| 日韩欧美一区视频在线观看 | 国产69精品久久久久777片| 欧美区成人在线视频| 久久久久人妻精品一区果冻| 国产高清三级在线| 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 国产精品国产av在线观看| 又大又黄又爽视频免费| 狂野欧美激情性bbbbbb| 伦理电影大哥的女人| 色哟哟·www| 十八禁高潮呻吟视频 | 水蜜桃什么品种好| 伦理电影大哥的女人| 美女xxoo啪啪120秒动态图| 一级爰片在线观看| 亚洲真实伦在线观看| 美女cb高潮喷水在线观看| 久热这里只有精品99| 女人精品久久久久毛片| 三级国产精品片| 天堂中文最新版在线下载| 综合色丁香网| h视频一区二区三区| 一区二区av电影网| 国产精品国产三级国产专区5o| 国产精品久久久久久久电影| 国产综合精华液| 成人综合一区亚洲| 亚洲伊人久久精品综合| 亚洲精品第二区| 久久99一区二区三区| 97在线视频观看| av在线老鸭窝| 久久国产乱子免费精品| 99久久综合免费| 久热这里只有精品99| 精品久久久久久久久亚洲| 免费大片18禁| 国产老妇伦熟女老妇高清| 亚洲经典国产精华液单| av天堂久久9| 国产深夜福利视频在线观看| 亚洲真实伦在线观看| av卡一久久| 成人美女网站在线观看视频| 成人国产麻豆网| 亚洲精品国产av成人精品| 美女中出高潮动态图| a级毛片在线看网站| 最近手机中文字幕大全| av专区在线播放| 久久99精品国语久久久| 视频中文字幕在线观看| 成人午夜精彩视频在线观看| 美女主播在线视频| 国产精品免费大片| 欧美成人午夜免费资源| 国产在线免费精品| 97超视频在线观看视频| 久久久久网色| 成人免费观看视频高清| 国产国拍精品亚洲av在线观看| 久久久久国产精品人妻一区二区| 日韩在线高清观看一区二区三区| 久久久久精品久久久久真实原创| 午夜激情福利司机影院| 免费大片黄手机在线观看| 欧美日韩在线观看h| 国产成人a∨麻豆精品| 亚洲,欧美,日韩| 国产熟女午夜一区二区三区 | 亚州av有码| 欧美日韩av久久| 另类亚洲欧美激情| 日本vs欧美在线观看视频 | 日韩,欧美,国产一区二区三区| 色哟哟·www| 欧美丝袜亚洲另类| 少妇高潮的动态图| 又粗又硬又长又爽又黄的视频| 在线亚洲精品国产二区图片欧美 | 中文乱码字字幕精品一区二区三区| 国产一级毛片在线| 成人免费观看视频高清| 免费播放大片免费观看视频在线观看| 中文乱码字字幕精品一区二区三区| 国产色爽女视频免费观看| 国产熟女午夜一区二区三区 | 啦啦啦视频在线资源免费观看| 亚洲人成网站在线播| 亚洲精品aⅴ在线观看| 91精品国产国语对白视频| 91精品一卡2卡3卡4卡| 少妇人妻久久综合中文| 成年人免费黄色播放视频 | 美女cb高潮喷水在线观看| 亚洲人成网站在线观看播放| 免费久久久久久久精品成人欧美视频 | 极品教师在线视频| 大陆偷拍与自拍| 久久久久视频综合| 久久精品国产亚洲网站| av有码第一页| 校园人妻丝袜中文字幕| 中国国产av一级| 午夜久久久在线观看| 韩国av在线不卡| 夫妻性生交免费视频一级片| 青青草视频在线视频观看| 91精品国产九色| 成人亚洲精品一区在线观看| 亚洲国产欧美在线一区| 高清在线视频一区二区三区| 欧美精品一区二区大全| 国国产精品蜜臀av免费| 高清毛片免费看| 成人毛片60女人毛片免费| 我要看黄色一级片免费的| 精品人妻熟女毛片av久久网站| 在线观看三级黄色| 99热这里只有是精品在线观看| 亚洲欧美一区二区三区黑人 | 久久久久久久久久人人人人人人| 国产精品女同一区二区软件| 国产高清三级在线| 国产黄频视频在线观看| 久久久久网色| 天天操日日干夜夜撸| 一二三四中文在线观看免费高清| 久久久久久久亚洲中文字幕| 中文欧美无线码| 亚洲av综合色区一区| 久久99一区二区三区| 欧美日韩av久久| 波野结衣二区三区在线| 卡戴珊不雅视频在线播放| 在线播放无遮挡| 国产日韩一区二区三区精品不卡 | 亚洲高清免费不卡视频| 伊人久久精品亚洲午夜| 国产亚洲91精品色在线| 男人狂女人下面高潮的视频| 亚洲精品日韩在线中文字幕| 久久人人爽人人片av| 99精国产麻豆久久婷婷| 免费观看的影片在线观看| 亚洲第一区二区三区不卡| 亚洲成人手机| 国产精品一区二区在线观看99| 久久ye,这里只有精品| 91精品一卡2卡3卡4卡| 免费大片18禁| 丝袜喷水一区| 18禁裸乳无遮挡动漫免费视频| 亚洲av成人精品一区久久| 成人国产麻豆网| 亚洲精品乱码久久久v下载方式| 欧美日韩视频精品一区| 欧美成人精品欧美一级黄| 午夜日本视频在线| 国产乱来视频区| 日韩人妻高清精品专区| 女性被躁到高潮视频| 国产欧美亚洲国产| 伦理电影免费视频| 国产伦理片在线播放av一区| 欧美国产精品一级二级三级 | 日韩av不卡免费在线播放| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 国产高清国产精品国产三级| 色5月婷婷丁香| 自拍欧美九色日韩亚洲蝌蚪91 | 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 久久99一区二区三区| 99久久精品国产国产毛片| 成年美女黄网站色视频大全免费 | 偷拍熟女少妇极品色| 免费少妇av软件| 黑人猛操日本美女一级片| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 丝瓜视频免费看黄片| 我的女老师完整版在线观看| 亚洲不卡免费看| 亚洲一级一片aⅴ在线观看| 亚洲精品中文字幕在线视频 | 国产精品99久久99久久久不卡 | 97精品久久久久久久久久精品| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 高清午夜精品一区二区三区| videossex国产| av.在线天堂| 午夜91福利影院| 黄色欧美视频在线观看| 又大又黄又爽视频免费| 久久久国产欧美日韩av| 日本vs欧美在线观看视频 | 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 婷婷色麻豆天堂久久| 麻豆成人午夜福利视频| 日本-黄色视频高清免费观看| 下体分泌物呈黄色| 久久久久久久久久久久大奶| 国产无遮挡羞羞视频在线观看| 亚洲欧美精品专区久久| 亚洲欧美中文字幕日韩二区| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 久久热精品热| 男人狂女人下面高潮的视频| 黄片无遮挡物在线观看| 亚洲国产色片| 久久久久久久久久久久大奶| 国产无遮挡羞羞视频在线观看| 日本av免费视频播放| 97在线人人人人妻| 婷婷色综合大香蕉| 国产视频内射| 精品国产一区二区久久| 韩国高清视频一区二区三区| av福利片在线| 国产精品伦人一区二区| 日本wwww免费看| 精品亚洲成a人片在线观看| 又粗又硬又长又爽又黄的视频| 极品少妇高潮喷水抽搐| 多毛熟女@视频| a级一级毛片免费在线观看| 久久久久久久精品精品| 午夜福利,免费看| 男男h啪啪无遮挡| 大陆偷拍与自拍| 国产黄色视频一区二区在线观看| 在线观看av片永久免费下载| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 国产男女内射视频| 又黄又爽又刺激的免费视频.|