CN104808512A - 一种航天器多级驱动刚柔耦合响应的获取方法 - Google Patents
一种航天器多级驱动刚柔耦合响应的获取方法 Download PDFInfo
- Publication number
- CN104808512A CN104808512A CN201510095003.2A CN201510095003A CN104808512A CN 104808512 A CN104808512 A CN 104808512A CN 201510095003 A CN201510095003 A CN 201510095003A CN 104808512 A CN104808512 A CN 104808512A
- Authority
- CN
- China
- Prior art keywords
- omega
- overbar
- parts
- joint
- matrix
- 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.)
- Granted
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
一种航天器多级驱动刚柔耦合响应的获取方法,步骤为:(1)将航天器视为由本体、部件和关节三类结构组成;(2)分别建立本体、部件和关节的动能和势能;(3)根据连接条件,建立整个航天器的动能和势能;(4)利用步骤(3)得到的航天器动能和势能,采用Lagrange方程建立航天器的动力学方程;(5)求解步骤(4)得到的动力学方程,获得多级驱动过程中航天器的刚柔耦合动力学响应。本发明方法适用于具有柔性多体运动特征的复杂航天器,获取这类复杂多体结构航天器运动过程精确动力学响应,对控制***设计和***级仿真验证具有重要的意义。
Description
技术领域
本发明属于航天器动力学领域,涉及一种针对柔性多体类航天器的动力学响应获取方法。
背景技术
随着载人航天和通信技术的发展,航天器正逐渐由单轴驱动太阳翼形式向着多体方向发展,例如载人空间机械臂、双轴驱动太阳翼、多自由度驱动天线等。这种刚体和柔性体的混合***,在各体之间存在大范围的转动驱动运动的同时,还叠加有柔性体的弹性振动。这种复杂的刚柔耦合动力学效应,以及与控制***的控制-结构交叉耦合,使得这类复杂航天器在动力学特性分析,特别是GNC控制方案验证时对柔性多体***的动力学模型及其解算需求十分迫切。对于这种刚-柔混合多体***,基于刚体假设的动力学模型已无法描述其复杂的动力学性态,必须采取柔性多体动力学建模方法,同时考虑部件大范围运动和部件本身的变形。
20世纪70年代中期,有关柔性多体***动力学的理论研究工作实际上已经展开。到目前为止,柔性多体***动力学的研究虽然取得了不少成果,但是还没有达到多刚体***动力学的研究水平,其主要原因是对物体大范围运动与弹性变形耦合问题的认识与处理方法上遇到困难。而目前的商业软件都是封装形式,只能进行响应计算,而无法给出解析式的动力学方程及其中间计算过程,导致无法提供控制***设计所需的解析式模型,以及无法开展动力学与控制的联合仿真。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种航天器多级驱动刚柔耦合响应的获取方法,通过给出一套通用化的柔性多体动力学解析式模型,并对模型求解获取航天器的动力学响应,解决此类航天器复杂的机构运动与弹性振动的耦合动力学解算问题,可以为刚柔耦合多体结构航天器的动力学建模、刚柔耦合响应计算、控制***设计,与控制***的联合仿真提供技术支持。
本发明的技术解决方案是:一种航天器多级驱动刚柔耦合响应的获取方法,包括如下步骤:
(1)将航天器视为由本体、部件和关节三类结构组成,其中本体和部件之间,以及部件和部件之间通过关节连接,或者固定连接,部件为刚体或者柔性体,关节为带驱动功能的柔性体,包括电机和减速机构;
(2)获取本体的动能Tb,
其中,vb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的速度,在ObXbYbZb下描述,ωb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的角速度,在ObXbYbZb下描述,Mb是本体的质量,Pb是本体相对于本体参考系ObXbYbZb的静矩,Ib是本体相对于本体参考系ObXbYbZb的转动惯量;
(3)获取部件的动能Tai和势能Vai,
其中,ηai是部件i的正则模态集的模态坐标,ξai,1是部件i的约束模态集相对于部件i-1连接位置的界面节点物理坐标,ξai,2是部件i的约束模态集相对于部件i+1连接位置的界面节点物理坐标,和中的各矩阵是质量阵和刚度阵根据ηai、ξai,1和ξai,2划分后的分块矩阵;
(4)获取关节的动能Tri和势能Vti,
其中,vti是关节i电机转子参考系OtiXtiYtiZti相对于惯性系OiXiYiZi的速度,ωti是关节i电机转子参考系OtiXtiYtiZti相对于惯性坐标系OiXiYiZi的角速度,Mti是关节i电机转子质量,Pti是关节i电机转子相对于电机转子参考系OtiXtiYtiZti的静矩,Iti是关节i电机转子相对于关节i电机转子参考系OtiXtiYtiZti的转动惯量,ξti是关节i减速机构与电机转子连接点相对于关节i电机转子参考系OtiXtiYtiZti的弹性变形位移,Kti是减速机构的刚度阵,kti,tx,kti,ty,kti,tz是3个平动刚度值,kti,rx,kti,ry,kti,rz是3个转动刚度值;
(5)根据连接条件,结合步骤(2)~(4)的结果,建立整个航天器的动能T和势能V,
其中Ms是航天器的质量,Is是航天器相对于本体参考系ObXbYbZb的转动惯量,Ps是航天器相对于本体参考系ObXbYbZb的静矩, ω0=[vb ωb]T是本体平动速度和转动角速度组成的向量, 是各部件或各关节的单轴转动角速度,q=[q1 q2 ... qn]T是各部件或各关节的单轴转动角,为包含各部件的正则模态坐标、界面点坐标、关节弹性变形位移的速度项,I0(q)、I1(q)、I2(q)、I10(q)、I20(q)、I21(q)是广义质量阵I关于的分块矩阵;
所述的连接条件为:
其中,ra1是部件1与本体固定连接时部件1与本体连接点在ObXbYbZb下的位置矢量,是ra1的斜对称矩阵,是部件1与本体通过关节1连接时关节1参考系Ot1Xt1Yt1Zt1相对于本体参考系ObXbYbZb的角速度,在Ot1Xt1Yt1Zt1下描述,rt1是部件1与关节1连接点在Ot1Xt1Yt1Zt1下的位置矢量,是rt1的斜对称矩阵,Cb→t1是ObXbYbZb到Ot1Xt1Yt1Zt1的坐标变换矩阵,Ct1→a1是Ot1Xt1Yt1Zt1到部件1参考系Oa1Xa1Ya1Za1的坐标变换矩阵,是关节i参考系OtiXtiYtiZti相对于部件i-1参考系Oai-1Xai-1Yai-1Zai-1的角速度,在OtiXtiYtiZti下描述,rai是部件i与部件i-1的连接点在Oai-1Xai-1Yai-1Zai-1下的位置矢量,是rai的斜对称矩阵,rti是部件i与关节i连接点在OtiXtiYtiZti下的位置矢量,是rti的斜对称矩阵,Cai-1→ti是Oai-1Xai-1Yai-1Zai-1到OtiXtiYtiZti的坐标变换矩阵,Cti→ai是OtiXtiYtiZti到OaiXaiYaiZai的坐标变换矩阵,是ξai-1,1的平动分量,是ξai-1,1的转动分量;
(6)利用步骤(5)得到的航天器的动能和势能,采用Lagrange方程建立航天器的动力学方程,
T=[Fb Tb τ1 τ2 ... τn]T
其中是ωb的斜对称矩阵,是vb的斜对称矩阵, 是广义质量阵I关于ω的分块矩阵,ξ是阻尼系数,Ω是振动频率矩阵,是I对时间t的导数,是I对qi的偏导数,是qi对时间t的导数,Fb,Tb是本体作用力和力矩,τ1,τ2,...τn是各部件电子转子驱动力矩;
(7)求解步骤(6)得到的动力学方程,获得航天器多级驱动过程中的刚柔耦合动力学响应。
本发明与现有技术相比的优点在于:本发明采用Lagrange方程的模态综合-混合坐标动力学建模方法。首先,将航天器划分为本体、部件和关节,分别建立其动能和势能;其次,根据本体、部件和关节之间的连接关系,建立整个航天器的动能和势能;最后,利用Lagrange方程,建立航天器动力学方程,积分求解获得航天器多级驱动刚柔耦合响应。本发明方法提供了用于控制***设计的解析式动力学模型,并且能够精确的给出刚柔耦合响应结果,实现航天器多级驱动动力学与控制联合仿真。
附图说明
图1为本发明树形拓扑结构示意图;
图2为本发明柔性关节建模示意图;
图3为本发明方法的流程框图;
图4为本发明柔性体连接示意图。
具体实施方式
本发明所研究的多级驱动航天器是指航天器除本体外,带有多轴驱动附件,如多轴驱动天线、两轴太阳翼、空间机械臂等,这类航天器具有多体运动特征,同时部件和关节存在弹性变形。本文所指的关节为带驱动功能的机构,其主要由电机和减速机构组成。根据图1,本发明将航天器***划分为本体和附件(如太阳翼、天线、机械臂等);附件由若干个部件(如机械臂的臂杆)和关节组成,本体和部件,以及部件之间可以通过关节连接,或固定连接;部件可以考虑为刚体,或柔性体。本发明中,首先对部件、关节建模,其次给出本体、部件和关节之间的连接条件,再次推导出***动能和势能,最后得到***动力学方程并求解,得到航天器的耦合响应。
整个发明的实施流程如图3所示。具体实施方式详述如下:
(1)部件柔性建模
针对图3中的第i个部件,下面建立其动力学模型,其振动方程可以写为
其中m为质量矩阵,k为刚度矩阵,f为包含有与第i-1个结构和第i+1个结构连接处的内力,x是节点响应。
当部件i与前后部件连接时,需要考虑其连接界面点,如图4中的ξai,1和ξai,2。因此,根据内部坐标和界面坐标,将方程(1)划分为
其中xi为内部节点响应,xj为界面节点响应,fj为耦合界面连接内力,mii,mij,mjj为质量矩阵m按照内、外部节点划分后的分块矩阵,kii,kij,kjj为刚度矩阵k按照内、外部节点划分后的分块矩阵。
对于固定界面有xj=0,由方程(2)可获得子结构自由振动方程为
于是,可求得子结构第n阶正则模态φn,满足
其中ωn为固定界面第n阶频率。
假设φii=[φ1,φ2,…,φn],则子结构的正则模态集为
对式(2),其静力方程为
由方程(7)可以得到
令得子结构的约束模态集
其中E是单位矩阵。
利用正则模态集Φn和约束模态集Φc组成子结构模态矩阵
Φ=[Φn Φc] (10)
将部件i物理坐标按照内外节点变换到模态坐标下有
x=Φn,aiηai+Φc,ai,1ξai,1+Φc,ai,2ξai,2 (11)
其中,ηai是部件i的正则模态集的模态坐标,ξai,1是部件i的约束模态集相对于部件i-1连接位置的界面节点物理坐标,ξai,2是部件i的约束模态集相对于部件i+1连接位置的界面节点物理坐标,Φn,ai是部件i的正则模态集,Φc,ai,1是部件i关于界面点ξai,1的约束模态集,Φc,ai,2是部件i关于界面点ξai,2的约束模态集。
那么,根据式(11),将式(2)划分为
其中和是模型(2)经过变换式(11)减缩后的质量阵和刚度阵,式(13)中的各矩阵是质量阵和刚度阵根据ηai、ξai,1和ξai,2划分后的分块矩阵。
将界面点ξai,2变换为与界面点ξai,1的相对运动关系
ξai,2=Tξai,1+ξai,2 (14)
其中rai是部件i与部件i-1的连接位置矢量,是rai的斜对称矩阵。
那么,利用变换和分别左乘和右乘方程(12)
变换后,最终获得柔性部件动力学方程(12)的系数矩阵
当部件作为刚体考虑时,方程(12)退化为
(2)关节柔性建模
对于由电机和减速机构组成的驱动机构(本发明中称为关节),本发明将其模化为电机转子+弹簧系,其中电机转子为刚体,描述关节的质量特性,弹簧系描述减速机构的柔性,如图2所示,弹簧系由3个线弹簧和3个卷簧组成。在***建模中,将电机转子按动能考虑为
其中,vti是关节i电机转子参考系OtiXtiYtiZti相对于惯性系OiXiYiZi的速度,ωti是电机转子参考系OtiXtiYtiZti相对于惯性坐标系OiXiYiZi的角速度。Mti是电机转子质量,Pti是电机转子相对于其参考系OtiXtiYtiZti的静矩,Iti是电机转子相对于其参考系OtiXtiYtiZti的转动惯量。电机转子参考系OtiXtiYtiZti的原点为电机转子连接端面的形心,三轴指向根据具体建模的需求进行定义。
在***建模中,将电机转子按势能考虑为
其中,ξti是图1中的弹簧系与驱动部件连接点相对于电机转子参考系的弹性变形位移,是弹簧系的刚度阵,kti,tx,kti,ty,kti,tz是3个平动刚度值,kti,rx,kti,ry,kti,rz是3个转动刚度值。
式(20)和式(21)为最终获得的柔性关节模型,当不考虑关节的柔性时,关节模型退化为式(20)。
(3)建立整个航天器的动能和势能
本发明所考虑的航天器力学模型如图1所示的树形拓扑结构。其中,OiXiYiZi是惯性参考系,ObXbYbZb是本体参考系,OaiXaiYaiZai是部件i参考系,OtiXtiYtiZti是关节i参考系。本体参考系ObXbYbZb原点位于本体质心,三轴指向根据具体建模的需求进行定义。部件参考系OaiXaiYaiZai原点位于部件与前一个部件的连接点,三轴指向根据具体建模的需求进行定义。
采用Lagrange方程建立航天器动力学方程,下面首先建立各个体的动能和势能表达式,建模过程中不考虑重力影响。
本体动能表达式为
其中,vb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的速度,在ObXbYbZb下描述,ωb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的角速度,在ObXbYbZb下描述。Mb是本体质量,Pb是本体相对于本体参考系ObXbYbZb的静矩,Ib是本体相对于本体参考系ObXbYbZb的转动惯量。
类似地,根据式(17)和式(18),各部件动能和势能表达式为
其中界面点ξai,1代表部件i的刚体运动,界面点ξai,2代表部件i与部件i+1连接点相对于部件i参考系OaiXaiYaiZai的振动物理位移,代表部件i的正则模态坐标。
当考虑关节柔性时,其动能和势能表达式如式(20)和式(21)所示。
下面考虑本体、各部件和关节之间的连接条件。根据考虑关节柔性与否,本体和部件1的连接关系分为如下两种情况:
1)不考虑关节柔性
其中,ωa1是部件1参考系Oa1Xa1Ya1Za1相对于本体参考系ObXbYbZb的角速度,在Oa1Xa1Ya1Za1下描述,ra1是部件1与本体连接点在ObXbYbZb下的位置矢量,是ra1的斜对称矩阵,Cb→a1是ObXbYbZb到Oa1Xa1Ya1Za1的坐标变换矩阵。
2)考虑关节柔性
其中,是关节1参考系Ot1Xt1Yt1Zt1相对于本体参考系ObXbYbZb的角速度,在Ot1Xt1Yt1Zt1下描述,rt1是部件1与关节1连接点在Ot1Xt1Yt1Zt1下的位置矢量,是rt1的斜对称矩阵,Cb→t1是ObXbYbZb到Ot1Xt1Yt1Zt1的坐标变换矩阵,Ct1→a1是Ot1Xt1Yt1Zt1到Oa1Xa1Ya1Za1的坐标变换矩阵。
类似地,根据考虑关节柔性与否,部件i和部件i-1的连接关系分为如下两种情况:
1)不考虑关节柔性
其中,是ξai-1,1的平动分量,是ξai-1,1的转动分量,ωai是部件i参考系OaiXaiYaiZai相对于部件i-1参考系Oai-1Xai-1Yai-1Zai-1的角速度,在OaiXaiYaiZai下描述,rai是部件i与部件i-1的连接点在Oai-1Xai-1Yai-1Zai-1下的位置矢量,是rai的斜对称矩阵,Cai-1→ai是Oai-1Xai-1Yai-1Zai-1到OaiXaiYaiZai的坐标变换矩阵。
2)考虑关节柔性
其中,是关节i参考系OtiXtiYtiZti相对于部件i-1参考系Oai-1Xai-1Yai-1Zai-1的角速度,在OtiXtiYtiZti下描述,rti是部件i与关节i连接点在OtiXtiYtiZti下的位置矢量,是rti的斜对称矩阵,Cai-1→ti是Oai-1Xai-1Yai-1Zai-1到OtiXtiYtiZti的坐标变换矩阵,Cti→ai是OtiXtiYtiZti到OaiXaiYaiZai的坐标变换矩阵。
最后,根据连接条件(26)~(31),得到***的动能和势能表达式
T=Tb+Ta1+Ta2+...+Tan (32)
V=Va1+Va2+...+Van (33)
(4)建立航天器动力学方程
下面根据(3)推导的***动能方程和势能方程,利用Lagrange方程,推导***动力学方程。将***动能T和***势能V描述如下
其中,***坐标 ω0=[vb ωb]T是本体平动速度和转动角速度, 是各部件或各关节的单轴转动角速度,q=[q1 q2 ... qn]T是各部件或各关节的单轴转动角,为包含各部件的正则模态坐标、界面点坐标、关节弹性变形位移的速度项。Ωi是第i阶振动频率矩阵。I0(q)、I1(q)、I2(q)、I10(q)、I20(q)、I21(q)是***广义质量阵I关于的分块矩阵。
广义坐标和伪坐标形式的Lagrange方程如下
其中Qk为对应全部主动力的广义力。
其中,Fb为本体所受外力,Tb为本体所受外力矩,是ωb的斜对称矩阵。
利用式(36)和式(37),分别对式(34)和(35)中的本体伪坐标和广义坐标求解,得到***动力学方程
T=[Fb Tb τ1 τ2 ... τn]T (40)
其中是ωb的斜对称矩阵,是vb的斜对称矩阵, ξ是阻尼系数,Ω是振动频率矩阵,是I对时间t的导数,是I对qi的偏导数,是q对时间t的导数,τ1,τ2,...τn是各部件驱动力矩。
综上所述,按照步骤(1)~(4)建立起来的方程(38)为最终获得的航天器多级驱动刚柔耦合动力学模型,对方程(38)进行数值积分运算,即可获得多级驱动过程中的刚柔耦合动力学响应。对方程(38)和控制***进行编程实现,闭环求解即可实现航天器的动力学与控制联合仿真。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (1)
1.一种航天器多级驱动刚柔耦合响应的获取方法,其特征在于包括如下步骤:
(1)将航天器视为由本体、部件和关节三类结构组成,其中本体和部件之间,以及部件和部件之间通过关节连接,或者固定连接,部件为刚体或者柔性体,关节为带驱动功能的柔性体,包括电机和减速机构;
(2)获取本体的动能Tb,
其中,vb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的速度,在ObXbYbZb下描述,ωb是本体参考系ObXbYbZb相对于惯性坐标系OiXiYiZi的角速度,在ObXbYbZb下描述,Mb是本体的质量,Pb是本体相对于本体参考系ObXbYbZb的静矩,Ib是本体相对于本体参考系ObXbYbZb的转动惯量;
(3)获取部件的动能Tai和势能Vai,
其中,ηai是部件i的正则模态集的模态坐标,ξai,1是部件i的约束模态集相对于部件i-1连接位置的界面节点物理坐标,ξai,2是部件i的约束模态集相对于部件i+1连接位置的界面节点物理坐标,和中的各矩阵是质量阵和刚度阵根据ηai、ξai,1和ξai,2划分后的分块矩阵;
(4)获取关节的动能Tri和势能Vti,
其中,vti是关节i电机转子参考系OtiXtiYtiZti相对于惯性系OiXiYiZi的速度,ωti是关节i电机转子参考系OtiXtiYtiZti相对于惯性坐标系OiXiYiZi的角速度,Mti是关节i电机转子质量,Pti是关节i电机转子相对于电机转子参考系OtiXtiYtiZti的静矩,Iti是关节i电机转子相对于关节i电机转子参考系OtiXtiYtiZti的转动惯量,ξti是关节i减速机构与电机转子连接点相对于关节i电机转子参考系OtiXtiYtiZti的弹性变形位移,是减速机构的刚度阵,kti,tx,kti,ty,kti,tz是3个平动刚度值,kti,rx,kti,ry,kti,rz是3个转动刚度值;
(5)根据连接条件,结合步骤(2)~(4)的结果,建立整个航天器的动能T和势能V,
其中Ms是航天器的质量,Is是航天器相对于本体参考系ObXbYbZb的转动惯量,Ps是航天器相对于本体参考系ObXbYbZb的静矩,ω0=[vb ωb]T是本体平动速度和转动角速度组成的向量,是各部件或各关节的单轴转动角速度,q=[q1 q2 ... qn]T是各部件或各关节的单轴转动角,为包含各部件的正则模态坐标、界面点坐标、关节弹性变形位移的速度项,I0(q)、I1(q)、I2(q)、I10(q)、I20(q)、I21(q)是广义质量阵I关于的分块矩阵;
所述的连接条件为:
其中,ra1是部件1与本体固定连接时部件1与本体连接点在ObXbYbZb下的位置矢量,是ra1的斜对称矩阵,是部件1与本体通过关节1连接时关节1参考系Ot1Xt1Yt1Zt1相对于本体参考系ObXbYbZb的角速度,在Ot1Xt1Yt1Zt1下描述,rt1是部件1与关节1连接点在Ot1Xt1Yt1Zt1下的位置矢量,是rt1的斜对称矩阵,Cb→t1是ObXbYbZb到Ot1Xt1Yt1Zt1的坐标变换矩阵,Ct1→a1是Ot1Xt1Yt1Zt1到部件1参考系Oa1Xa1Ya1Za1的坐标变换矩阵,是关节i参考系OtiXtiYtiZti相对于部件i-1参考系Oai-1Xai-1Yai-1Zai-1的角速度,在OtiXtiYtiZti下描述,rai是部件i与部件i-1的连接点在Oai-1Xai-1Yai-1Zai-1下的位置矢量,是rai的斜对称矩阵,rti是部件i与关节i连接点在OtiXtiYtiZti下的位置矢量,是rti的斜对称矩阵,Cai-1→ti是Oai-1Xai-1Yai-1Zai-1到OtiXtiYtiZti的坐标变换矩阵,Cti→ai是OtiXtiYtiZti到OaiXaiYaiZai的坐标变换矩阵,是ξai-1,1的平动分量,是ξai-1,1的转动分量;
(6)利用步骤(5)得到的航天器的动能和势能,采用Lagrange方程建立航天器的动力学方程,
T=[Fb Tb τ1 τ2 ...τn]T
其中是ωb的斜对称矩阵,是vb的斜对称矩阵,是广义质量阵I关于ω的分块矩阵,ξ是阻尼系数,Ω是振动频率矩阵,是I对时间t的导数,是I对qi的偏导数,是qi对时间t的导数,Fb,Tb是本体作用力和力矩,τ1,τ2,...τn是各部件电子转子驱动力矩;
(7)求解步骤(6)得到的动力学方程,获得航天器多级驱动过程中的刚柔耦合动力学响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510095003.2A CN104808512B (zh) | 2015-03-03 | 2015-03-03 | 一种航天器多级驱动刚柔耦合响应的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510095003.2A CN104808512B (zh) | 2015-03-03 | 2015-03-03 | 一种航天器多级驱动刚柔耦合响应的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104808512A true CN104808512A (zh) | 2015-07-29 |
CN104808512B CN104808512B (zh) | 2017-05-31 |
Family
ID=53693447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510095003.2A Active CN104808512B (zh) | 2015-03-03 | 2015-03-03 | 一种航天器多级驱动刚柔耦合响应的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104808512B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631167A (zh) * | 2016-03-03 | 2016-06-01 | 北京空间飞行器总体设计部 | 一种航天器热致振动动力学响应评估方法 |
CN105843074A (zh) * | 2016-03-28 | 2016-08-10 | 北京空间飞行器总体设计部 | 一种获取天线在轨振动影响的动力学建模方法 |
CN106707211A (zh) * | 2016-11-30 | 2017-05-24 | 上海卫星工程研究所 | 针对全极化微波辐射计角速度波动量的高精度测试方法 |
CN106777587A (zh) * | 2016-11-30 | 2017-05-31 | 上海卫星工程研究所 | 针对空间环境下挠性附件振动仿真的测量方法 |
CN107186548A (zh) * | 2017-06-08 | 2017-09-22 | 大连理工大学 | 一种五轴数控机床回转轴几何误差检测方法 |
CN107436605A (zh) * | 2017-07-26 | 2017-12-05 | 西北工业大学 | 考虑多体耦合的水下机器人控制方法 |
CN108388699A (zh) * | 2018-01-29 | 2018-08-10 | 南京理工大学 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
CN111814375A (zh) * | 2020-07-08 | 2020-10-23 | 中国电子科技集团公司第五十四研究所 | 基于机电一体化联合仿真的天线伺服***优化设计方法 |
US11400590B2 (en) | 2017-05-10 | 2022-08-02 | Moog Bv | Optimal control of coupled admittance controllers |
CN117921745A (zh) * | 2024-03-25 | 2024-04-26 | 中国科学院长春光学精密机械与物理研究所 | 多方向运动转换的时变刚度基座*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020116078A1 (en) * | 2001-02-20 | 2002-08-22 | Best System Inc. | Control method of open type motion simulation system |
CN101526801A (zh) * | 2009-04-17 | 2009-09-09 | 天津理工大学 | 履带式移动机械手的广义***统一动力学建模方法 |
CN103412491A (zh) * | 2013-08-27 | 2013-11-27 | 北京理工大学 | 一种挠性航天器特征轴姿态机动指数时变滑模控制方法 |
-
2015
- 2015-03-03 CN CN201510095003.2A patent/CN104808512B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020116078A1 (en) * | 2001-02-20 | 2002-08-22 | Best System Inc. | Control method of open type motion simulation system |
CN101526801A (zh) * | 2009-04-17 | 2009-09-09 | 天津理工大学 | 履带式移动机械手的广义***统一动力学建模方法 |
CN103412491A (zh) * | 2013-08-27 | 2013-11-27 | 北京理工大学 | 一种挠性航天器特征轴姿态机动指数时变滑模控制方法 |
Non-Patent Citations (3)
Title |
---|
史纪鑫等: "复合柔性结构航天器动力学降阶模型研究", 《200全国结构动力学学术研讨会》 * |
曲广吉: "柔性航天器动力学分析", 《航天器工程》 * |
缪炳祺等: "柔性航天器动力学模型的模态综合-混合坐标法建模研究", 《航天器工程》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631167A (zh) * | 2016-03-03 | 2016-06-01 | 北京空间飞行器总体设计部 | 一种航天器热致振动动力学响应评估方法 |
CN105631167B (zh) * | 2016-03-03 | 2018-08-14 | 北京空间飞行器总体设计部 | 一种航天器热致振动动力学响应评估方法 |
CN105843074B (zh) * | 2016-03-28 | 2018-08-14 | 北京空间飞行器总体设计部 | 一种获取天线在轨振动影响的动力学建模方法 |
CN105843074A (zh) * | 2016-03-28 | 2016-08-10 | 北京空间飞行器总体设计部 | 一种获取天线在轨振动影响的动力学建模方法 |
CN106707211A (zh) * | 2016-11-30 | 2017-05-24 | 上海卫星工程研究所 | 针对全极化微波辐射计角速度波动量的高精度测试方法 |
CN106777587A (zh) * | 2016-11-30 | 2017-05-31 | 上海卫星工程研究所 | 针对空间环境下挠性附件振动仿真的测量方法 |
CN106707211B (zh) * | 2016-11-30 | 2019-08-02 | 上海卫星工程研究所 | 针对全极化微波辐射计角速度波动量的高精度测试方法 |
US11400590B2 (en) | 2017-05-10 | 2022-08-02 | Moog Bv | Optimal control of coupled admittance controllers |
CN107186548A (zh) * | 2017-06-08 | 2017-09-22 | 大连理工大学 | 一种五轴数控机床回转轴几何误差检测方法 |
CN107186548B (zh) * | 2017-06-08 | 2018-12-18 | 大连理工大学 | 一种五轴数控机床回转轴几何误差检测方法 |
CN107436605B (zh) * | 2017-07-26 | 2020-11-27 | 西北工业大学 | 考虑多体耦合的水下机器人控制方法 |
CN107436605A (zh) * | 2017-07-26 | 2017-12-05 | 西北工业大学 | 考虑多体耦合的水下机器人控制方法 |
CN108388699A (zh) * | 2018-01-29 | 2018-08-10 | 南京理工大学 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
CN108388699B (zh) * | 2018-01-29 | 2022-02-18 | 南京理工大学 | 中心刚体-fgm楔形梁***末端动力学响应计算方法 |
CN111814375A (zh) * | 2020-07-08 | 2020-10-23 | 中国电子科技集团公司第五十四研究所 | 基于机电一体化联合仿真的天线伺服***优化设计方法 |
CN117921745A (zh) * | 2024-03-25 | 2024-04-26 | 中国科学院长春光学精密机械与物理研究所 | 多方向运动转换的时变刚度基座*** |
CN117921745B (zh) * | 2024-03-25 | 2024-05-24 | 中国科学院长春光学精密机械与物理研究所 | 多方向运动转换的时变刚度基座*** |
Also Published As
Publication number | Publication date |
---|---|
CN104808512B (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104808512B (zh) | 一种航天器多级驱动刚柔耦合响应的获取方法 | |
CN104723340B (zh) | 基于连接和阻尼配置的柔性关节机械臂的阻抗控制方法 | |
CN108714896B (zh) | 一种基于柔体动力学模型的机器人柔性力矩前馈控制方法 | |
CN103213129B (zh) | 一种空间机械臂位置力混合控制方法 | |
CN102207988B (zh) | 一种多自由度机械臂高效动力学建模方法 | |
CN104723341A (zh) | 基于连接和阻尼配置的柔性关节机械臂的位置控制方法 | |
CN107263466B (zh) | 空间机器人基于二次规划问题的基座无扰控制方法 | |
CN103869704B (zh) | 基于扩展雅克比矩阵的空间机器人星臂协调控制方法 | |
CN103076807B (zh) | 一种欠驱动柔性航天器姿态稳定的控制方法 | |
CN105160051A (zh) | 基于多体分析试验的桁架天线反射器展开动力学建模方法 | |
CN105479459A (zh) | 机器人零力控制方法以及*** | |
CN104656447A (zh) | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 | |
CN109522622B (zh) | 一种多自由度太阳电池阵在轨载荷工况确定方法及*** | |
Hu et al. | Maneuver and vibration control of flexible manipulators using variable-speed control moment gyros | |
CN103399986A (zh) | 基于微分几何的空间机械臂建模方法 | |
CN104267614A (zh) | 一种无人机实时仿真***及其开发方法 | |
CN102495550A (zh) | 并联机器人的正、逆动力学响应分析与控制方法 | |
CN109657282B (zh) | 一种基于拉格朗日动力学的h型运动平台建模方法 | |
Rong et al. | New efficient method for dynamic modeling and simulation of flexible multibody systems moving in plane | |
CN103577244A (zh) | 负载模拟器的速度同步控制方法和*** | |
CN105912007A (zh) | 空间机械臂抗干扰姿态稳定的微分几何非线性控制方法 | |
CN106249616A (zh) | 一种在轨服务机械臂动力学建模方法和*** | |
CN104742151A (zh) | 模块化双自由度球形关节及蛇形机器人及运动控制方法 | |
CN107894775B (zh) | 一种欠驱动无人水下航行器轨迹生成与控制方法 | |
CN104793497A (zh) | 基于多体***离散时间传递矩阵法的机器人动力学建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |