CN106529066A - Dynamic modelling method for connection process of indexing mechanism capture connection device - Google Patents

Dynamic modelling method for connection process of indexing mechanism capture connection device Download PDF

Info

Publication number
CN106529066A
CN106529066A CN201611033381.9A CN201611033381A CN106529066A CN 106529066 A CN106529066 A CN 106529066A CN 201611033381 A CN201611033381 A CN 201611033381A CN 106529066 A CN106529066 A CN 106529066A
Authority
CN
China
Prior art keywords
indexing
base
capture
head
centerdot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201611033381.9A
Other languages
Chinese (zh)
Inventor
王寅
陈金宝
王治易
耿海峰
王小涛
王道波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201611033381.9A priority Critical patent/CN106529066A/en
Publication of CN106529066A publication Critical patent/CN106529066A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种转位机构捕获连接装置连接过程动力学建模方法,属于空间站在轨服务技术领域。首先,确定转位机构捕获连接装置的结构和参数;确定连接过程中捕获头与转位基座之间相互接触滑移的运动约束;其次,确定捕获头与转位基座之间的接触点在转位机构捕获连接装置上的相对滑移速度;再次,根据捕获头与转位基座连接过程的约束动力学方程;最后,解析约束动力学方程。将转位机构捕获连接装置连接过程复杂的运动模态,转化成了空间任意运动的直线与直线之间以及直线与支架之间相互接触滑动的运动形式,极大的降低了转位机构捕获连接装置连接过程动力学分析的难度。

The invention discloses a dynamic modeling method for the connection process of a capture connection device of an indexing mechanism, and belongs to the technical field of space station orbit services. First, determine the structure and parameters of the capture connection device of the indexing mechanism; determine the motion constraints of the mutual contact and slippage between the capture head and the index base during the connection process; secondly, determine the contact point between the capture head and the index base The relative sliding velocity on the connecting device is captured in the indexing mechanism; again, according to the constraint dynamic equation of the connection process between the capture head and the index base; finally, the constraint dynamic equation is analyzed. The complex motion mode of the connecting process captured by the indexing mechanism is transformed into a movement form of contact and sliding between the straight line and the straight line and between the straight line and the support in any space, which greatly reduces the capture and connection of the indexing mechanism. Difficulty in kinetic analysis of device connection process.

Description

转位机构捕获连接装置连接过程动力学建模方法Dynamic modeling method for the connection process of the indexing mechanism capture connection device

技术领域technical field

本发明涉及一种动力学建模方法,具体讲是一种转位机构捕获连接装置连接过程动力学建模方法,属于空间站在轨服务技术领域。The invention relates to a dynamic modeling method, in particular to a dynamic modeling method for the connection process of an indexing mechanism capturing a connecting device, and belongs to the technical field of space station orbit services.

背景技术Background technique

空间站转位机构捕获连接装置能够用于大型航天器与空间站对接后在轨运行所需的组装、对接和分离等操作活动,是载人航天的重要组成部分之一。转位机构由旋转臂和旋转基座组成。旋转臂安装于待转舱体上,转位基座安装于空间站节点舱上,两者能实现对接、锁紧、解锁与复位。转臂捕获头与转位基座连接过程的动力学与运动学特性复杂,当导向瓣发生接触后,存在空间上的运动约束。常规的基于牛顿力学的建模方法不能直接应用于含有滑动约束的运动过程建模,制约了转位机构捕获连接装置的研究发展。The capture and connection device of the space station indexing mechanism can be used for the assembly, docking and separation required for the on-orbit operation of large spacecraft after docking with the space station, and is one of the important components of manned spaceflight. The indexing mechanism consists of a rotating arm and a rotating base. The rotating arm is installed on the body to be turned, and the indexing base is installed on the node module of the space station. The two can realize docking, locking, unlocking and reset. The dynamics and kinematics characteristics of the connection process between the capture head of the rotating arm and the indexing base are complex, and there are spatial motion constraints when the guide flap comes into contact. Conventional modeling methods based on Newtonian mechanics cannot be directly applied to the modeling of motion processes with sliding constraints, which restricts the research and development of indexing mechanism capture connection devices.

发明内容Contents of the invention

本发明所要解决的技术问题在于克服现有技术缺陷,提供一种能够准确描述转位机构捕获连接装置在连接过程中各组成部分之间相互作用下的动力学关系的建模方法。The technical problem to be solved by the present invention is to overcome the defects of the prior art and provide a modeling method capable of accurately describing the dynamic relationship of the indexing mechanism capturing the interaction between the various components of the connection device during the connection process.

为了解决上述技术问题,本发明提供的转位机构捕获连接装置连接过程动力学建模方法,包括以下步骤:In order to solve the above-mentioned technical problems, the dynamic modeling method for the connection process of the indexing mechanism capture connection device provided by the present invention includes the following steps:

步骤1:确定转位机构捕获连接装置的结构和参数,所述转位机构捕获连接装置包括捕获头、捕获杆和转位基座,所述捕获头与捕获杆相连接,所述捕获头与转位基座可相对移动;所述捕获头端面的导向装置与转位基座均为圆环形结构,沿外周向均匀间隔设有3个向外展开的梯形导向瓣,捕获头与转位基座之间各导向瓣的位置一一对应;Step 1: Determine the structure and parameters of the capture and connection device of the indexing mechanism. The capture and connection device of the indexing mechanism includes a capture head, a capture rod and an index base, the capture head is connected with the capture rod, and the capture head is connected with the capture rod. The indexing base can be moved relatively; the guide device on the end face of the capturing head and the indexing base are circular structures, and three trapezoidal guide lobes are evenly spaced along the outer circumference, and the capturing head and the indexing The positions of the guide petals between the bases are in one-to-one correspondence;

步骤2:确定连接过程中捕获头与转位基座之间相互接触滑移的运动约束;Step 2: Determine the motion constraints of the mutual contact slip between the capture head and the indexing base during the connection process;

步骤3:确定捕获头与转位基座之间的接触点在转位机构捕获连接装置上的相对滑移速度;Step 3: Determine the relative sliding speed of the contact point between the capturing head and the indexing base on the capturing connection device of the indexing mechanism;

步骤4:根据步骤2和3确定捕获头与转位基座连接过程的约束动力学方程;Step 4: According to steps 2 and 3, determine the constraint dynamic equation of the connection process between the capture head and the indexing base;

步骤5:解析约束动力学方程。Step 5: Solve the constrained dynamic equations.

本发明中,所述步骤2为:设捕获头与转位基座接触时,第i个接触点的绝对速度分别为:In the present invention, the step 2 is: when the capture head is in contact with the indexing base, the absolute speeds of the i-th contact point are respectively:

其中,表示捕获头的第i个接触点的绝对速度,表示转位基座的第i个接触点的绝对速度,ω1表示捕获头导向瓣相对于捕获头连接面坐标系的旋转角速度,ω2表示转位基座导向瓣在转位基座连接面坐标系内的旋转角速度,ω3表示由捕获头连接面与转位基座连接面所构成的局部坐标系相对于惯性系的旋转角速度,捕获头连接面中心到第i个接触点的位置矢量,转位基座连接面中心到第i个接触点的位置矢量;表示捕获头与转位基座导向瓣在接触点1,3的绝对速度,表示捕获头与转位基座导向瓣在接触点2,4的绝对速度。in, represents the absolute velocity of the i-th contact point of the capture head, Indicates the absolute velocity of the i-th contact point of the indexing base, ω 1 indicates the rotational angular velocity of the guide lobe of the capture head relative to the coordinate system of the interface of the capture head, ω 2 indicates that the guide lobe of the index base is on the interface of the index base The angular velocity of rotation in the coordinate system, ω 3 represents the angular velocity of rotation of the local coordinate system formed by the connecting surface of the capture head and the connecting surface of the indexing base relative to the inertial system, The position vector from the center of the connecting surface of the capture head to the i-th contact point, The position vector from the center of the connecting surface of the indexing base to the i-th contact point; Indicates the absolute velocity of the capture head and the guide lobe of the indexing base at the contact point 1,3, Indicates the absolute velocity of the capture head at the point of contact with the guide lobe of the indexing base2,4.

本发明中,所述步骤3为:当捕获头边缘g7,q7均为顺时针方向时,In the present invention, the step 3 is: when the edges g 7 and q 7 of the capture head are both clockwise,

其中,分别表示与捕获头和转位基座相关的变量,其中下标数字表示接触的边缘;in, denote the variables associated with the capture head and the indexing base, respectively, where the subscript numbers denote the edges of contact;

因此,第i个接触滑移点速度约束方程为:Therefore, the velocity constraint equation of the i-th contact slip point is:

其中,表示捕获头连接面坐标系相对于捕获头连接面与转位基座坐连接面构成的局部坐标的旋转角速度,表示捕获头与转位基座在第2个接触点的绝对速度,表示转位基座导向瓣在转位基座连接面坐标系内的旋转角速度;in, Indicates the rotational angular velocity of the capture head connection surface coordinate system relative to the local coordinates formed by the capture head connection surface and the index base seat connection surface, Indicates the absolute velocity of the capture head and the indexing base at the second contact point, and indicates the rotational angular velocity of the guide lobe of the indexing base in the coordinate system of the connecting surface of the indexing base;

捕获头的导向瓣与基座第i个接触点的加速度约束为:The acceleration constraint of the i-th contact point between the guide flap of the capture head and the base is:

本发明中,所述步骤4为:In the present invention, the step 4 is:

其中,m1表示转臂的质量,m2表示转位基座所在机构的质量,m3表示转臂捕获头的质量,m4表示转位基座的质量,AI1表示惯性系向捕获头连接面坐标系转换的旋转矩阵,A1I表示捕获头连接面坐标系向惯性系转换的旋转矩阵,m表示接触点的个数,第i个接触点的约束力,表示转臂在惯性系下的等效张量矩阵、表示捕获头与转位基座在接触点1处的加速度、捕获头导向瓣与转位基座导向瓣在接触点1,3接触处的加速度、AI3表示惯性系向局部坐标系变换的旋转矩,为捕获头与转位基座之间的作用力、A1I捕获头连接面坐标系向惯性系变换的旋转矩阵、为捕获头与转位基座之间的作用力矩、A3I表示局部坐标系向惯性系变换的旋转矩阵,表示捕获头与转位基座在接触点2处的加速度,A24表示转位基座连接面坐标系向惯性系变换的转换矩阵、A42表示惯性系变换向基座连接面坐标系的转换矩阵,A2I同A24,Fctrl为作用在转位基座上的主动控制力,pi为转位基座与捕获头在e3坐标系内描述的位置矢量,ai为转位基座与捕获头之间接触点在e4坐标系内秒速的矢量位置,I3转位基座在e3坐标系内描述的等效惯性张量,I1为转位基座在e1坐标系内描述的等效惯性张量。Among them, m 1 represents the mass of the rotating arm, m 2 represents the mass of the mechanism where the indexing base is located, m 3 represents the mass of the capturing head of the rotating arm, m 4 represents the mass of the indexing base, and A I1 represents the inertial system towards the capturing head The rotation matrix of the connection surface coordinate system transformation, A 1I represents the rotation matrix of the capture head connection surface coordinate system to the inertial system transformation, m represents the number of contact points, the binding force at the i-th contact point, represents the equivalent tensor matrix of the arm in the inertial system, Indicates the acceleration of the capture head and the indexing base at the contact point 1, The acceleration of the guide lobe of the capture head and the guide lobe of the indexing base at the contact points 1 and 3, A I3 represents the rotational moment transformed from the inertial system to the local coordinate system, is the force between the capture head and the indexing base, the rotation matrix of A 1I transformation from the coordinate system of the connection surface of the capture head to the inertial system, is the acting moment between the capture head and the indexing base, A 3I represents the rotation matrix transformed from the local coordinate system to the inertial system, Indicates the acceleration of the capture head and the indexing base at the contact point 2, A 24 indicates the transformation matrix from the coordinate system of the joint surface of the indexing base to the inertial system, and A 42 indicates the transformation from the inertial system transformation to the coordinate system of the joint surface of the base Matrix, A 2I is the same as A 24 , F ctrl is the active control force acting on the indexing base, p i is the position vector described by the indexing base and the capture head in the e 3 coordinate system, a i is the indexing The contact point between the base and the capture head is the vector position of the second speed in the e 4 coordinate system, I 3 is the equivalent inertia tensor described by the indexing base in the e 3 coordinate system, and I 1 is the indexing base in the e 1 The equivalent inertial tensor described in the coordinate system.

5、根据权利要求1至4任一项所述的转位机构捕获连接装置连接过程动力学建模方法,其特征在于所述步骤5为:5. The dynamic modeling method for the connection process of the indexing mechanism capture connection device according to any one of claims 1 to 4, characterized in that the step 5 is:

步骤51:将步骤4捕获头与转位基座的约束动力学方程设为矩阵:Step 51: Set the constraint dynamic equation of the capture head and the indexing base in step 4 as a matrix:

其中,M表示捕获头与转位基座连接后所形成整体的质量、AT约束方程系数的转置矩阵;Among them, M represents the mass of the whole formed by the connection between the capture head and the indexing base, and the transpose matrix of the coefficients of the A T constraint equation;

捕获头与转位基座之间的接触滑移的运动约束,Kinematic constraints for the contact slip between the capture head and the indexing base,

式中,为整个系统位形坐标一阶倒数列阵,F表示转位基座与捕获头之间的约束力,f=(f1,f2…fm)为捕获头与转位基座的接触约束力,A为约束方程系数矩阵;In the formula, is the first-order reciprocal array of the configuration coordinates of the entire system, F represents the constraint force between the indexing base and the capture head, f=(f 1 ,f 2 …f m ) is the contact constraint between the capture head and the indexing base Force, A is the constraint equation coefficient matrix;

步骤52:构造约束动力学方程的A矩阵的正交补矩阵:Step 52: Construct the orthogonal complement matrix of the A matrix constraining the dynamic equation:

G=ATA (14)G=A T A (14)

式中,AT为A的转置矩阵,由式(9)可得G为n×n方阵,其秩为m;In the formula, AT is the transposed matrix of A, from formula (9), it can be obtained that G is an n×n square matrix, and its rank is m;

由|λΕ-G|=0得到G的特征值,利用该特征值构造矩阵Get the eigenvalue of G by |λΕ-G|=0, use the eigenvalue to construct the matrix

式中,Ε表示矩阵G的特征值,C’和C分别表示由非零特征值所对应特征向量和零特征值所对应特征向量所构成的矩阵,矩阵L满秩In the formula, E represents the eigenvalue of the matrix G, C' and C represent the matrix composed of the eigenvectors corresponding to the non-zero eigenvalues and the eigenvectors corresponding to the zero eigenvalues, respectively, and the matrix L is full rank

δX=LδZ (16)δX=LδZ (16)

其中,δX和δZ之间存在一个单值一一对应关系,由此in, There is a single-value one-to-one correspondence between δX and δZ, thus

因为C矩阵的每一列分别是由矩阵G的特征值的特征向量组成,因此Since each column of the C matrix is composed of the eigenvectors of the eigenvalues of the matrix G respectively, so

GC=0 (18)GC=0 (18)

代入(14)式Substitute into (14) formula

由此可以得到From this you can get

AC=0 (20)AC = 0 (20)

矩阵C即为矩阵A的正交补矩阵,因此方程(19)中的δZ”取任意值均成立,有Matrix C is the orthogonal complement matrix of matrix A, so any value of δZ” in Equation (19) is valid, and we have

把(16)式可写成Formula (16) can be written as

δX=CδZ” (22)δX=CδZ" (22)

把(22)式代入(12)式可得Substitute (22) into (12) to get

结合δZ”的任意性可得Combined with the arbitrariness of δZ", we get

约束动力学的正交补方程组为The orthogonal complementary equations for constrained dynamics are

本发明将转位机构捕获连接装置的动力学模型简化为由三个个体组成的多体系统,即捕获头、捕获杆和转位基座。在建立动力学模型时,考虑捕获头与转位基座之间的接触滑移约束的复杂性和不确定性,将捕获头与转位基座分开考虑,将接触点的约束力作为外力分别作用与捕获头和转位基座连接面上。此外,不考虑导向瓣在校正过程中的弹性变形,并利用虚功原理(Jourdain-Bertrand原理),对转位机构捕获连接装置的连接过程进行动力学建模。The invention simplifies the dynamic model of the capture connection device of the index mechanism into a multi-body system composed of three individuals, that is, the capture head, the capture rod and the index base. When establishing the dynamic model, considering the complexity and uncertainty of the contact slip constraint between the capture head and the indexing base, the capture head and the indexing base are considered separately, and the constraint force at the contact point is taken as an external force, respectively It acts on the connecting surface of the capture head and the indexing base. In addition, without considering the elastic deformation of the guide flap during the correction process, and using the principle of virtual work (Jourdain-Bertrand principle), the dynamic modeling of the connection process of the indexing mechanism capturing the connecting device is carried out.

本发明的有益效果在于:(1)、将转位机构捕获连接装置连接过程复杂的运动模态,转化成了空间任意运动的直线与直线之间以及直线与支架之间相互接触滑动的运动形式,极大的降低了转位机构捕获连接装置连接过程动力学分析的难度;(2)、通过分析上述两种运动模态的运动约束,给出了其相互接触点的相对滑移速度的解析表达式,解决了转位机构捕获连接装置在连接过程动力学特性的解析描述。The beneficial effects of the present invention are as follows: (1), the complex motion mode of the connection process captured by the indexing mechanism is transformed into a movement form of mutual contact and sliding between straight lines and between straight lines and supports in any space. , which greatly reduces the difficulty of dynamic analysis of the connection process of the indexing mechanism to capture the connection device; (2), through the analysis of the motion constraints of the above two motion modes, the analysis of the relative slip velocity of their mutual contact points is given The expression solves the analytical description of the indexing mechanism capturing the dynamic characteristics of the connection device during the connection process.

附图说明Description of drawings

图1为转位机构捕获连接装置结构示意图;Fig. 1 is a schematic diagram of the structure of the capture and connection device of the indexing mechanism;

图2为捕获头端面导向装置结构示意图;(a)为俯视图、(b)为侧视图;Fig. 2 is a schematic diagram of the structure of the end surface guide device of the capture head; (a) is a top view, (b) is a side view;

图3为转位机构基座结构示意图;(a)为俯视图、b)为侧视图。Fig. 3 is a schematic diagram of the base structure of the indexing mechanism; (a) is a top view, and b) is a side view.

具体实施方式detailed description

下面结合附图对本发明作进一步详细说明。The present invention will be described in further detail below in conjunction with the accompanying drawings.

如图1所述,本发明转位机构捕获连接装置连接过程动力学建模方法中采用的转位机构捕获连接装置捕获头12、捕获杆11和转位基座14。捕获头12与捕获杆11相连接,捕获头12的前端面安装导向装置13,捕获头12与转位基座14可相对移动,并进行接触。导向装置13和转位基座14的外周安装向外展开的导向瓣14。As shown in FIG. 1 , the indexing mechanism used in the dynamic modeling method for the connection process of the indexing mechanism capturing connecting device of the present invention captures the capturing head 12 , the capturing rod 11 and the indexing base 14 of the connecting device. The capture head 12 is connected with the capture rod 11, the front end of the capture head 12 is equipped with a guide device 13, and the capture head 12 and the indexing base 14 can move relatively and make contact. The outer periphery of the guide device 13 and the indexing base 14 is equipped with guide petals 14 that expand outwards.

如图2所示,本发明捕获头的导向装置为圆环形结构,捕获头的导向装置沿外周向均匀间隔设有3个向外展开的导向瓣,导向瓣为梯形结构。捕获头导向装置坐标系e1,图中D为导向装置内环直径,Rk为导向装置外环半径,β为导向瓣相对于导向装置端面的安装角,H为导向瓣高度,L为导向瓣长度。捕获头导向装置的导向瓣边缘编号依次为1,2,3,…6,端面边缘的编号7。图中,x1、y1和z1表示捕获头连接面坐标系的方向。As shown in Fig. 2, the guide device of the capture head of the present invention is a circular structure, and the guide device of the capture head is provided with three outwardly expanding guide petals evenly spaced along the outer circumference, and the guide petals are trapezoidal structures. Coordinate system e 1 of the guide device of the capture head, D in the figure is the diameter of the inner ring of the guide device, R k is the radius of the outer ring of the guide device, β is the installation angle of the guide flap relative to the end face of the guide device, H is the height of the guide flap, and L is the guide petal length. The numbering of the guide flap edge of the capture head guide device is 1, 2, 3, ... 6 in sequence, and the number of the edge of the end face is 7. In the figure, x1, y1 and z1 represent the direction of the coordinate system of the connection surface of the capture head.

如图3所示,本发明转位基座与捕获头的导向装置的结构一致,其为坐标系e2(即表示转位基座连接面坐标系)。转位基座上的3个导向瓣与捕获头导向装置上3个导向瓣的位置一一相对应。转位基座上的导向瓣边缘编号依次为1,2,3,…6,端面边缘的编号7。图中,x2,y2和z2表示转位基座连接面坐标系。As shown in FIG. 3 , the structure of the indexing base of the present invention is consistent with the guide device of the capture head, which is the coordinate system e 2 (that is, the coordinate system of the connecting surface of the indexing base). The positions of the three guide petals on the indexing base correspond to the positions of the three guide petals on the guide device of the capture head. The numbering of the edge of the guide flap on the transposition base is 1, 2, 3, ... 6, and the number of the edge of the end face is 7. In the figure, x2, y2 and z2 represent the coordinate system of the connecting surface of the indexing base.

本发明转位机构捕获连接装置连接过程动力学建模方法的具体步骤如下:The specific steps of the dynamic modeling method for the connection process of the indexing mechanism capture connection device of the present invention are as follows:

步骤1:确定转位机构捕获连接装置的结构和结构参数,具体结合及附图1-3。Step 1: Determine the structure and structural parameters of the indexing mechanism capture connection device, specifically combined with Figures 1-3.

步骤2:确定连接过程中捕获头与转位基座之间相互接触滑移的运动约束:当转位机构捕获连接装置的捕获头与转位基座在捕获连接过程中,捕获头与转位基座锁定并回收后,由于导向瓣的结构特点使得接触点向减小捕获头和转位基座之间相对位置和姿态偏差的方向滑动。因此,建立这种接触约束条件下的动态约束是建立转位机构捕获连接装置捕获连接过程动力学模型的首要前提。Step 2: Determine the motion constraints of the mutual contact and slippage between the capture head and the index base during the connection process: when the capture head and the index base of the connection device are captured by the index mechanism during the capture connection process, the capture head and the index After the base is locked and retrieved, due to the structural characteristics of the guide flap, the contact point slides in the direction of reducing the relative position and attitude deviation between the capture head and the indexing base. Therefore, the establishment of dynamic constraints under such contact constraints is the first prerequisite for establishing the dynamic model of the capture connection process of the indexing mechanism capture connection device.

根据转位机构捕获连接装置的结构特点,将捕获头与转位基座导向瓣与导向瓣之间的接触约束,简化为空间两条运动直线的相互接触滑移。假设捕获头与转位基座接触时,第i个接触点的绝对速度分别为:According to the structural characteristics of the capture and connection device of the indexing mechanism, the contact constraint between the capture head and the guide lobe of the index base is simplified as the mutual contact and sliding of two moving straight lines in space. Assuming that when the capture head is in contact with the indexing base, the absolute velocities of the i-th contact point are:

其中,分别表示捕获头与转位基座的第i个接触点的绝对速度,ω1表示捕获头导向瓣相对于捕获头连接面坐标系的旋转角速度,ω2表示转位基座导向瓣在转位基座连接面坐标系内的旋转角速度,ω3表示由捕获头连接面与转位基座连接面所构成的局部坐标系相对于惯性系的旋转角速度,捕获头连接面中心到第i个接触点的位置矢量,转位基座连接面中心到第i个接触点的位置矢量;表示捕获头与转位基座导向瓣在接触点1,3的绝对速度,表示捕获头与转位基座导向瓣在接触点2,4的绝对速度。in, and Respectively represent the absolute velocity of the i-th contact point between the capture head and the indexing base, ω 1 represents the rotational angular velocity of the guide lobe of the capture head relative to the coordinate system of the connection surface of the capture head, and ω 2 represents the guide lobe of the index base when indexing The angular velocity of rotation in the coordinate system of the base connecting surface, ω 3 represents the angular velocity of rotation of the local coordinate system formed by the connecting surface of the capture head and the connecting surface of the indexing base relative to the inertial system, The position vector from the center of the connecting surface of the capture head to the i-th contact point, The position vector from the center of the connecting surface of the indexing base to the i-th contact point; Indicates the absolute velocity of the capture head and the guide lobe of the indexing base at the contact point 1,3, Indicates the absolute velocity of the capture head at the point of contact with the guide lobe of the indexing base2,4.

步骤3:确定约束接触点在转位机构捕获连接装置上的相对滑移速度:Step 3: Determine the relative slip velocity of the constrained contact point on the indexing mechanism capture linkage:

当捕获头边缘g7,q7均为顺时针方向时,When the edges g 7 and q 7 of the capture head are both clockwise,

其中,分别表示与捕获头和基座相关的变量,其中下标数字表示接触的边缘。in, Denote the variables associated with the capture head and base, respectively, where the subscript numbers indicate the edge of contact.

因此,第i个接触滑移点速度约束方程为:Therefore, the velocity constraint equation of the i-th contact slip point is:

其中,表示捕获头连接面坐标系相对于捕获头连接面与转位基座坐连接面构成的局部坐标的旋转角速度,表示捕获头与转位基座在第2个接触点的绝对速度;in, Indicates the rotational angular velocity of the capture head connection surface coordinate system relative to the local coordinates formed by the capture head connection surface and the index base seat connection surface, Indicates the absolute speed of the capture head and the indexing base at the second contact point;

导向瓣与转位基座第i个接触点的加速度约束为:The acceleration constraint of the i-th contact point between the guide flap and the indexing base is:

步骤4:确定转位机构捕获连接装置连接过程动力学方程Step 4: Determine the kinetic equation of the connection process of the indexing mechanism capture connection device

其中,m1表示转臂的质量,m2表示转位基座所在机构的质量,m3表示转臂捕获头的质量,m4表示转位基座的质量,AI1表示惯性系向捕获头连接面坐标系转换的旋转矩阵,A1I表示捕获头连接面坐标系向惯性系转换的旋转矩阵,m(无下标时)表示接触点的个数,第i个接触点的约束力,表示转臂在惯性系下的等效张量矩阵、表示捕获头与转位基座在接触点1处的加速度、捕获头导向瓣与转位基座导向瓣在接触点1,3接触处的加速度、AI3表示惯性系向局部坐标系变换的旋转矩,为捕获头与转位基座之间的作用力、A1I捕获头连接面坐标系向惯性系变换的旋转矩阵、第i个接触点的约束力,为捕获头与转位基座之间的作用力矩、A3I表示局部坐标系向惯性系变换的旋转矩阵,表示捕获头与转位基座在接触点2处的加速度,A24表示转位基座连接面坐标系向惯性系变换的转换矩阵、A42表示惯性系变换向基座连接面坐标系的转换矩阵,A2I同A24,Fctrl为作用在转位基座上的主动控制力,pi为转位基座与捕获头在e3坐标系(即捕获头连接面与转位基座连接面所构成的坐标系)内描述的位置矢量;ai为转位基座与捕获头之间接触点在e4坐标系(即惯性坐标系)内秒速的矢量位置,I3转位基座在e3坐标系(即捕获头连接面与转位基座连接面所构成的坐标系)内描述的等效惯性张量,I1为转位基座在e1坐标系(即捕获头连接面坐标系)内描述的等效惯性张量,如图1所示。Among them, m 1 represents the mass of the rotating arm, m 2 represents the mass of the mechanism where the indexing base is located, m 3 represents the mass of the capturing head of the rotating arm, m 4 represents the mass of the indexing base, and A I1 represents the inertial system towards the capturing head The rotation matrix of the connection surface coordinate system transformation, A 1I represents the rotation matrix of the capture head connection surface coordinate system to the inertial system transformation, m (without subscript) represents the number of contact points, the binding force at the i-th contact point, represents the equivalent tensor matrix of the arm in the inertial system, Indicates the acceleration of the capture head and the indexing base at the contact point 1, The acceleration of the guide lobe of the capture head and the guide lobe of the indexing base at the contact points 1 and 3, A I3 represents the rotational moment transformed from the inertial system to the local coordinate system, is the force between the capture head and the indexing base, the rotation matrix of A 1I transformation from the coordinate system of the connection surface of the capture head to the inertial system, The binding force of the i-th contact point, is the acting moment between the capture head and the indexing base, A 3I represents the rotation matrix transformed from the local coordinate system to the inertial system, Indicates the acceleration of the capture head and the indexing base at the contact point 2, A 24 indicates the transformation matrix from the coordinate system of the joint surface of the indexing base to the inertial system, and A 42 indicates the transformation from the inertial system transformation to the coordinate system of the joint surface of the base Matrix, A 2I is the same as A 24 , F ctrl is the active control force acting on the transposition base, p i is the e 3 coordinate system between the transposition base and the capture head (that is, the connecting surface of the capture head is connected to the transposition base The position vector described in the coordinate system formed by the surface); a i is the vector position of the second speed of the contact point between the indexing base and the capture head in the e 4 coordinate system (that is, the inertial coordinate system), and I 3 the indexing base The equivalent inertia tensor described in the e 3 coordinate system (that is, the coordinate system formed by the connecting surface of the capture head and the connecting surface of the indexing base), I 1 is the index of the indexing base in the e 1 coordinate system (that is, the connection surface of the capturing head The equivalent inertia tensor described in the surface coordinate system) is shown in Fig. 1.

步骤5:确定约束动力学方程的解析形式Step 5: Determine the Analytical Form of the Constrained Kinetic Equations

将转位机构的捕获头与转位基座的动力学方程写成矩阵形式Write the dynamic equations of the capturing head and the indexing base of the indexing mechanism in matrix form

其中,M表示捕获头与转位基座连接后所形成整体的质量、AT约束方程系数的转置矩阵;Among them, M represents the mass of the whole formed by the connection between the capture head and the indexing base, and the transpose matrix of the coefficients of the A T constraint equation;

捕获头与转位基座之间的接触滑移运动约束,Contact-slip motion constraints between capture head and index base,

上式中,为整个系统位形坐标一阶倒数列阵,F表示转位基座与捕获头之间的约束力,f=(f1,f2…fm)为捕获头与转位基座的接触约束力,A为约束方程系数矩阵。由于捕获头与转位基座在捕获连接过程中,捕获头与转位基座之间相互运动的空间性和空间作业的过程复杂,直接对式(10)和(11)的约束方程组进行求解。本发明运用D’Alembert正交补原理来处理变维约束问题,目的是将原变维动力学方程组转换为固定的n为位形空间微分方程组进行求解。根据虚功原理对X的虚位位移需满足:In the above formula, is the first-order reciprocal array of the configuration coordinates of the entire system, F represents the constraint force between the indexing base and the capture head, f=(f 1 ,f 2 …f m ) is the contact constraint between the capture head and the indexing base Force, A is the constraint equation coefficient matrix. Since the capture head and the indexing base are in the process of capturing and connecting, the mutual movement between the capturing head and the indexing base is complex and the process of the space operation is complicated, the constraint equations of equations (10) and (11) are directly calculated solve. The invention uses the D'Alembert orthogonal complement principle to deal with the variable-dimensional constraint problem, and aims to convert the original variable-dimensional dynamic equations into fixed n-configuration space differential equations for solution. According to the principle of virtual work, the virtual position displacement of X needs to satisfy:

约束动力学方程需要满足The constrained dynamic equations need to satisfy

AδX=0 (13)AδX=0 (13)

考虑到系统受到动态运动约束,所以X并不独立,为此需要通过约束方程的A矩阵的正交补矩阵进行求解,首先约束动力学方程A矩阵的正交补矩阵Considering that the system is constrained by dynamic motion, X is not independent. Therefore, it is necessary to solve it through the orthogonal complement matrix of the A matrix of the constraint equation. First, the orthogonal complement matrix of the A matrix of the constraint equation is constrained

G=ATA (14)G=A T A (14)

其中AT为A的转置矩阵,有由式(9)可得G为n×n方阵,其秩为m。因此,由|λΕ-G|=0可以得到G的特征值,利用特征值构造矩阵Among them, AT is the transposed matrix of A, and G can be obtained from formula (9) as n×n square matrix, and its rank is m. Therefore, the eigenvalue of G can be obtained by |λΕ-G|=0, and the matrix is constructed using the eigenvalue

其中,Ε表示矩阵G的特征值,C’和C分别表示由非零特征值所对应特征向量和零特征值所对应特征向量所构成的矩阵。因此矩阵L必定满秩Among them, E represents the eigenvalue of the matrix G, and C' and C represent the matrix formed by the eigenvectors corresponding to the non-zero eigenvalues and the eigenvectors corresponding to the zero eigenvalues, respectively. So the matrix L must be full rank

δX=LδZ (16)δX=LδZ (16)

其中,δZ,δX是构造的矩阵,仅有数学含义,无物理意义,δX和δZ之间存在一个单值一一对应关系,由此in, δZ and δX are constructed matrices, which have only mathematical meaning and no physical meaning. There is a single-value one-to-one correspondence between δX and δZ, thus

因为C矩阵的每一列分别是由矩阵G的特征值的特征向量组成,因此Since each column of the C matrix is composed of the eigenvectors of the eigenvalues of the matrix G respectively, so

GC=0 (18)GC=0 (18)

代入(14)式Substitute into (14) formula

由此可以得到From this you can get

AC=0 (20)AC = 0 (20)

矩阵C即为矩阵A的正交补矩阵,因此方程(19)中的δZ”取任意值均成立,有Matrix C is the orthogonal complement matrix of matrix A, so any value of δZ” in Equation (19) is valid, and we have

把(16)式可写成Formula (16) can be written as

δX=CδZ” (22)δX=CδZ" (22)

把上式代入(12)式可得Substituting the above formula into (12) can get

由于δZ”的任意性,可得Due to the arbitrariness of δZ", we can get

因此,约束动力学的正交补方程组Therefore, the system of orthogonal complement equations for constrained dynamics

式(25)即为在n维位形空间描述的约束动力学基本方程,而约束力f已不出现在方程中了。Equation (25) is the basic equation of constraint dynamics described in n-dimensional configuration space, and the constraint force f does not appear in the equation.

本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。There are many specific application approaches of the present invention, and the above description is only a preferred embodiment of the present invention. It should be pointed out that for those of ordinary skill in the art, some improvements can also be made without departing from the principles of the present invention. Improvements should also be regarded as the protection scope of the present invention.

Claims (5)

1. A dynamic modeling method for a connection process of a capture connection device of an indexing mechanism is characterized by comprising the following steps:
step 1: determining the structure and parameters of an indexing mechanism capturing connection device, wherein the indexing mechanism capturing connection device comprises a capturing head, a capturing rod and an indexing base, the capturing head is connected with the capturing rod, and the capturing head and the indexing base can move relatively; the guide device on the end face of the capturing head and the transposition base are both in circular structures, 3 trapezoidal guide flaps which are unfolded outwards are uniformly arranged at intervals along the outer periphery, and the positions of the guide flaps between the capturing head and the transposition base correspond to each other one by one;
step 2: determining the motion constraint of mutual contact and slippage between the capture head and the indexing base in the connection process;
and step 3: determining a relative sliding speed of a contact point between the capture head and the indexing base on the indexing mechanism capture connection device;
and 4, step 4: determining a constraint kinetic equation of the connection process of the capture head and the transposition base according to the steps 2 and 3;
and 5: and (6) resolving a constraint kinetic equation.
2. The method for modeling the dynamics of an indexing mechanism capture link assembly attachment process according to claim 1, wherein step 2 is:
when the capture head is in contact with the indexing base, the absolute speed of the ith contact point is respectively as follows:
R → i p = R → 1 + ω → 1 × ( R → 13 + p → i ) + R → 13 + ω → 3 × p → i + p → i R → i a = R → 1 + ω → 2 × ( a → i + R → 24 ) + a → i - - - ( 1 )
wherein,representing the absolute velocity of the ith contact point of the capture head,indicating the absolute speed, ω, of the ith contact point of the indexing carriage1Representing the angular velocity, ω, of rotation of the guide lobe of the capture head relative to the coordinate system of the connection face of the capture head2Representing the angular velocity, omega, of the guide lobe of the indexing base in the coordinate system of the surface of the indexing base connection3Represents the rotation angular velocity of a local coordinate system formed by the connecting surface of the capturing head and the connecting surface of the indexing base relative to the inertial system,capturing the position vector of the center of the head connection surface to the ith contact point,a position vector from the center of the connecting surface of the indexing base to the ith contact point;indicating the absolute velocity of the capture head and indexing base guide lobe at contact points 1,3,indicating the absolute velocity of the capture head at the contact points 2,4 with the indexing base guide lobes.
3. The method for modeling the dynamics of an indexing mechanism capture link assembly attachment process according to claim 2, wherein step 3 is:
when capturing the head edge g7,q7When the direction of the water is clockwise,
( R &RightArrow; i p - R &RightArrow; i a ) &CenterDot; ( g &RightArrow; i &times; q &RightArrow; i ) < 0 , ( i = 1 , 3 , 5 , 7 ... 12 ) ( R &RightArrow; i p - R &RightArrow; i a ) &CenterDot; ( g &RightArrow; i &times; q &RightArrow; i ) > 0 , ( i = 2 , 4 , 6 , 8 ... 18 ) - - - ( 2 )
wherein,representing variables associated with the capture head and the index base, respectively, where subscript numbers represent the edges of the contacts;
therefore, the ith contact slip point velocity constraint equation is:
wherein,representing the angular velocity of rotation of the capture head joint plane coordinate system relative to the local coordinate formed by the capture head joint plane and the indexing base joint plane,the absolute speed of the capture head and the indexing base at the 2 nd contact point is shown, and the rotation angular speed of the indexing base guide lobe in the indexing base connection surface coordinate system is shown;
the acceleration constraint of the guide lobe of the capture head to the ith contact point of the base is:
4. the method for modeling the dynamics of an indexing mechanism capture link assembly attachment process according to claim 3, wherein step 4 is:
( m 1 + m 2 ) R &RightArrow; &CenterDot;&CenterDot; 1 - m 3 A I 1 R ~ 13 &omega; &RightArrow; &CenterDot; 1 + m 3 A I 1 R ~ 13 = - &Sigma; i = 1 m F &RightArrow; i - m 3 A I 1 &omega; ~ 1 &omega; ~ 1 R ~ 13 - 2 m 3 A I 1 &omega; ~ 1 R ~ 13 + A I 1 F C t r l
( m 2 + m 4 ) R &RightArrow; &CenterDot;&CenterDot; 2 = &Sigma; i m F &RightArrow; i - - - ( 8 )
wherein m is1Representing the mass of the arm, m2Denotes the mass of the mechanism in which the index base is located, m3Representing the mass of the pivoted-arm catch head, m4Denotes the mass of the index pedestal, AI1Rotation matrix representing the transformation of the inertial system into the plane of coordinates of the attachment surface of the capture head, A1IA rotation matrix representing the conversion of the coordinate system of the connection surface of the capturing head to the inertial system, m represents the number of contact points,the restraining force of the ith point of contact,an equivalent tensor matrix representing the rotating arm under the inertial system,Representing the acceleration of the capture head at the contact point 1 with the indexing base,Acceleration of the Capture head guide lobe in contact with the indexing base guide lobe at contact points 1,3, AI3Representing the rotational moment of the transformation of the inertial system into the local coordinate system,for the force between the catch head and the indexing base, A1IA rotation matrix for transforming the coordinate system of the connection surface of the capture head to the inertial system,For the moment of action between the catch head and the indexing base, A3IA rotation matrix representing a transformation of the local coordinate system into the inertial system,representing the acceleration of the capture head at the contact point 2 with the indexing base, A24Transformation matrix representing transformation from coordinate system of connection surface of indexing base to inertial system, A42A transformation matrix representing the transformation of the inertial system into the coordinate system of the base connection surface, A2ISame as A24,FctrlFor active control forces acting on the indexing table, piFor indexing the base and the capture head at e3Position vector described in a coordinate system, aiFor the contact point between the indexing base and the capture head at e4Vector position in the coordinate system at second speed, I3The indexing base is at3Equivalent inertia tensor, I, described in a coordinate system1For indexing the base at e1An equivalent inertia tensor described in a coordinate system.
5. The method for modeling dynamics of an indexing mechanism capture link assembly connection process according to any of claims 1-4, wherein the step 5 is:
step 51: and 4, setting a constraint kinetic equation of the capture head and the indexing base in the step 4 as a matrix:
M X &RightArrow; &CenterDot;&CenterDot; = F - A T f - - - ( 10 )
wherein M represents the mass of the capture head as a whole when it is connected to the indexing base, ATA transposed matrix of constraint equation coefficients;
the motion constraint of the contact slip between the capture head and the index base,
A X &RightArrow; &CenterDot;&CenterDot; = 0 A X &CenterDot;&CenterDot; = - A &CenterDot; X &CenterDot; - - - ( 11 )
in the formula,for the whole system position coordinate first order reciprocal array, F represents the restraining force between the indexing base and the capture head, and F ═ F (F ═ F)1,f2…fm) A is a constraint equation coefficient matrix for the contact constraint force of the capture head and the transposition base;
step 52: constructing an orthogonal complement matrix of an A matrix of a constraint kinetic equation:
G=ATA (14)
in the formula, ATIs a transposed matrix of A, G is a square matrix of n × n obtained from formula (9), and the rank is m;
obtaining the eigenvalue of G from | Lambda E-G | ═ 0, and constructing matrix by using the eigenvalue
In the formula, e represents the eigenvalue of matrix G, C' and C represent the matrix formed by eigenvector corresponding to non-zero eigenvalue and eigenvector corresponding to zero eigenvalue, respectively, and matrix L full rank
X=LZ (16)
Wherein,there is a single value one-to-one correspondence between X and Z, whereby
A L &delta; Z = 0 AC &prime; &delta;Z &prime; + ACZ &prime; &prime; = 0 - - - ( 17 )
Since each column of the C matrix is composed of the eigenvectors of the eigenvalues of the matrix G, respectively
GC=0 (18)
Into formula (14)
A T A C = 0 , C T A T A C C = ( A C ) T A C = 0 - - - ( 19 )
Thereby can obtain
AC=0 (20)
The matrix C is an orthogonal complement of the matrix A, so that Z' in equation (19) holds for any value, some
AC &prime; &delta;Z &prime; &prime; = 0 &CenterDot; &CenterDot; &CenterDot; det ( AC &prime; ) &NotEqual; 0 &CenterDot; &CenterDot; &CenterDot; &delta;Z &prime; = 0 - - - ( 21 )
Can write formula (16) into
X=CZ” (22)
Substituting the formula (22) into the formula (12) to obtain
&delta;Z &prime; &prime; T C T ( F - A T f - m X &CenterDot;&CenterDot; ) = 0 - - - ( 23 )
Optional availability of binding Z ″)
C T ( F - A T f - m X &CenterDot;&CenterDot; ) = C T ( F - M X &CenterDot;&CenterDot; ) = 0 - - - ( 24 )
The orthogonal complement equation system of constraint dynamics is
C T M A X &CenterDot;&CenterDot; = C T F - A &CenterDot; X &CenterDot; - - - ( 25 ) .
CN201611033381.9A 2016-11-16 2016-11-16 Dynamic modelling method for connection process of indexing mechanism capture connection device Pending CN106529066A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611033381.9A CN106529066A (en) 2016-11-16 2016-11-16 Dynamic modelling method for connection process of indexing mechanism capture connection device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611033381.9A CN106529066A (en) 2016-11-16 2016-11-16 Dynamic modelling method for connection process of indexing mechanism capture connection device

Publications (1)

Publication Number Publication Date
CN106529066A true CN106529066A (en) 2017-03-22

Family

ID=58357654

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611033381.9A Pending CN106529066A (en) 2016-11-16 2016-11-16 Dynamic modelling method for connection process of indexing mechanism capture connection device

Country Status (1)

Country Link
CN (1) CN106529066A (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809276A (en) * 2015-04-14 2015-07-29 长安大学 Multi-finger robot dynamic analytical model and modeling method thereof

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104809276A (en) * 2015-04-14 2015-07-29 长安大学 Multi-finger robot dynamic analytical model and modeling method thereof

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
王兴贵 等: "周边式对接机构航天器的对接动力学仿真(第1部分)-运动学约束方程", 《系统仿真学报》 *
王兴贵 等: "周边式对接机构航天器的对接动力学仿真(第2部分)-约束动力学", 《系统仿真学报》 *
王翔 等: "对接过程中捕获锁质量和锁合面间隙影响分析", 《航天器工程》 *

Similar Documents

Publication Publication Date Title
Yoshida et al. Zero reaction maneuver: flight validation with ETS-VII space robot and extension to kinematically redundant arm
Ma et al. On-orbit identification of inertia properties of spacecraft using a robotic arm
CN101733749B (en) Multidomain uniform modeling and emulation system of space robot
CN101891018B (en) Single frame control moment gyro control method based on moment output capability optimization
CN104570742B (en) Feedforward PID (proportion, integration and differentiation) control based rapid high-precision relative pointing control method of noncoplanar rendezvous orbit
CN102759358B (en) Relative posture dynamics modeling method based on dead satellite surface reference points
CN106864776B (en) A kind of method and system of the capture target satellite based on butt joint ring
CN103955611B (en) Method for establishing universal mechanical model of tethered space robot
CN106289641A (en) Spacecraft centroid position and rotary inertia parametric joint discrimination method
CN104635740B (en) A kind of deep space probe autonomous attitude maneuver control method
CN105930627A (en) Free-floating space manipulator modeling method under condition of considering spacial flexible deformation of arm lever
CN106055901A (en) Method for determining timing of capturing tumbling target by free floating space robot
CN108469737A (en) A kind of dynamic control method and system of space non-cooperative target navigation capture
CN106249616A (en) One server mechanical arm dynamic modeling method and system in-orbit
Yu et al. A comparative study on motion characteristics of three two-degree-of-freedom pointing mechanisms
Pilz et al. Robust controller design for formation flight of quad-rotor helicopters
CN110641738B (en) Trajectory tracking control method of space five-degree-of-freedom free flying mechanical arm
Rosyid et al. A novel walking parallel robot for on-structure three-axis machining of large structures
Estrada et al. Force and moment constraints of a curved surface gripper and wrist for assistive free flyers
CN106529066A (en) Dynamic modelling method for connection process of indexing mechanism capture connection device
CN107807656B (en) Integrated modeling method for double-star formation
CN103438892B (en) A kind of astronomical autonomous orbit determination algorithm based on EKF of improvement
CN113060309A (en) A parameter identification method of space target under double-arm capture based on dynamic model
Shen et al. Translational and rotational maneuvers of an underactuated space robot using prismatic actuators
CN109975838A (en) A kind of world equivalent experiment method for the closed loop rope system fleet system that spins

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170322