CN112949104A - 一种协作机器人实时模态分析方法 - Google Patents
一种协作机器人实时模态分析方法 Download PDFInfo
- Publication number
- CN112949104A CN112949104A CN201911169700.2A CN201911169700A CN112949104A CN 112949104 A CN112949104 A CN 112949104A CN 201911169700 A CN201911169700 A CN 201911169700A CN 112949104 A CN112949104 A CN 112949104A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- matrix
- module
- robot
- real
- 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
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Numerical Control (AREA)
Abstract
Description
技术领域
本发明涉及协作机器人领域,具体地说是一种协作机器人实时模态分析方法。
背景技术
协作机器人是指能够在指定的协作区域内与人进行直接交互的机器人,其具有人机融合、安全易用、灵敏精准及灵活通用的特征。轻质、高负载自重比的设计理念使机器人引入了大量的柔性因素,导致机器人在运动过程中不可避免的产生振动,影响了机器人的动态性能和运动精度。机器人模态分析获得的固有频率、振幅及振型等动态特性是机器人结构优化、轨迹规划和振动抑制的重要参数,考虑结构形状的多位姿动态性能结构优化对机器人动态特性求解的实时性提出了要求,轨迹规划以及振动抑制对机器人的固有频率、振幅及振型等动态特性的求解精度及实时性要求较高,所以高精度实时的模态分析方法是实现协作机器人动态性能提高及振动抑制的重要理论基础,但关于机器人实时模态分析方法的研究目前还未见开展。
以往的模态分析方法主要有解析法、有限元法以及振动试验法。解析法是通过建立机器人的刚柔耦合或柔性动力学模型获得机器人***的无阻尼自由振动方程,进而得到机器人的固有频率和模态。解析法能够满足机器人动态性能结构优化和振动抑制对实时性的要求,因而广泛应用于并联、串联机器人的动态性能结构优化和振动抑制中,但该方法对机器人模型进行了大量的简化,将机器人关节简化为线性扭转弹簧,将机器人结构件简化为刚体或多段梁结构,虽然该方法降低了计算量,但关节及连杆的简化降低了计算精度,为提高人机协作的安全性和舒适性,协作机器人连杆通常采用光顺曲面设计,而将连杆过度简化等效为多段梁结构会造成建模精度降低。
高精度的模态分析方法包括有限元方法和振动实验法。有限元方法是基于有限元理论对机器人进行有限元建模、分析,从而获得机器人的固有频率及振型,振动实验法是利用频率响应函数(FRF)测量装置通过对机械结构进行激励分析获得***的固有频率、阻尼和振型等模态参数。有限元法可用于机器人设计阶段的动力学分析及结构优化等,在物理样机搭建完成后可通过振动实验法进行验证和有限元模型修正。有限元分析法和振动实验法能够直观真实的反映机械结构的动态性能,但均存在实时性差、工作量大的缺点,每一次分析结果只对所计算的位姿和边界条件有效,一旦位姿或边界条件改变,则需要重新建模或布置试验进行计算。采用上述方法分析机器人在工作空间内的动态特性,计算量巨大,难以实现实时模态分析。
发明内容
本发明的目的在于提供一种协作机器人实时模态分析方法,具有计算量小、效率高、精度高的优点,能够在保证精度的同时大幅提高实时性。
本发明的目的是通过以下技术方案来实现的:
一种协作机器人实时模态分析方法,包括如下步骤:
步骤一:将机器人拆分成多个独立的模块,并将各个模块等效为超单元;
步骤四:建立机械臂的无阻尼自由振动方程,并将步骤三中获得的整机刚度矩阵K和质量刚度矩阵M代入方程实现协作机器人固有频率及振型的实时求解,为协作机器人实时模态分析提供依据。
步骤一中,将机器人整机模型拆分为多个相互独立的连杆模块和关节模块,并将每个模块等效为具有两个主节点的超单元,各模块对应的超单元依次串联组成机械臂的等效有限元模型。
步骤二中,先通过三维建模软件对模块三维几何模型进行简化,再通过有限元分析软件完成各模块有限元的建模、参数设定以及单元矩阵的提取,并使用四面体单元划分有限元模型,然后采用有限元子结构法提取各个超单元在自身坐标系下的的刚度矩阵和质量矩阵
步骤三中:
3.1、先建立机器人各模块自身坐标系,并推导各模块运动关系;
上述式(1)中,表示第i个模块输入坐标系到第i-1个模块输入坐标系的齐次变换矩阵,表示第i个模块输入坐标系到第i-1个模块输入坐标系的旋转矩阵,[xi,yi,zi]分别表示第i个模块输入坐标系原点在第i-1个模块输入坐标系X,Y,Z三个方向的位置;
3.3、推导第i(i>1)个超单元自身坐标系到整机坐标系的变换矩阵,超单元自身坐标系到整机坐标系(整机坐标系与基座输入坐标系重合)的映射通过各超单元变换矩阵连乘获得:
3.4、根据3.3中超单元自身坐标系到整机坐标系的映射关系,获得第i个超单元到整机坐标系下的转换矩阵Si:
3.6、通过分块叠加原则形成整机刚度矩阵K和质量刚度矩阵M:
步骤3.1中,在各模块的输入法兰定义输入坐标系{Oi1},输出法兰定义输出坐标系{Oi2}用于描述各模块自身及相互之间的几何特征及运动特征。
步骤3.3中,机器人整机模型拆分为多个相互独立的连杆模块和关节模块,当第i个超单元为关节模块时,T为常数矩阵,当第i(i>1)个超单元为连杆模块时,T为关于关节运动变量的函数矩阵。
步骤四中,首先建立协作机器人无阻尼自由振动方程,方程可以表示为:
再将步骤三中获得整机刚度矩阵K和质量刚度矩阵M代入式(8),实现协作机器人固有频率及振型的实时求解。
步骤四中,由于自由振动分解为一系列简谐振动的叠加,因此所述式(8)的解设为:
q=Φejωt (9);
上述式(9)中:ω为简谐振动圆频率;Φ=[φ1,φ2…φn]为节点振幅列向量;φi(i=1,2…n)为自由度i方向上的振幅。
将式(9)代入式(8)并同时消去因子ejωt,可得:
(K-ω2M)×Φ=0 (10);
在进行边界条件设定后,将步骤三中获得的K、M代入公式(10)求解整机结构固有圆频率和模态振型,其中固有圆频率ωi(rad/s)与固有频率fi(Hz)的关系为:
ωi=2πfi (11)。
本发明的优点与积极效果为:
1、本发明具有计算量小、效率高、精度高的优点,能够全面考虑机器人各部件的弹性变形,弥补了以往分析方法的不足,相比于现有技术中的有限元法、振动实验法等,本发明方法实现了协作机器人固有频率及振型的实时求解,且能够在大幅提高实时性同时保证计算精度。本发明针对SHIR5型机器人进行了实验与仿真,结果表明,依据本发明的计算结果与仿真,和现场实验结果误差小于6%,处理周期约为5ms。
2、本发明采用有限元与解析法相结合的方法求解协作机器人动态性能等,为机器人结构优化、轨迹规划以及振动抑制提供了重要理论基础,并且可推广到其它串联、并联或串并联机器人***,为机器人***的设计提供了有效的分析理论和仿真手段。
附图说明
图1为本发明的方法流程示意图,
图2为本发明涉及的协作机器人简化模型示意图,
图3为本发明涉及的机器人关节模块示意图,
图4为本发明涉及的机器人连杆模块示意图,
图5为模块等效模型示意图,
图6为超单元模型示意图,
图7为提取单元刚度和质量矩阵的关节模块有限元模型示意图。
具体实施方式
下面结合附图对本发明作进一步详述。
如图1所示,本发明包括以下步骤:
步骤一:将机器人拆分成独立的连杆模块和关节模块,并将各个模块等效为超单元。
本步骤具体为:如图2所示,基于模块化设计理念的协作机器人,由于各关节和连杆模块的功能、机械接口、电气接口等相互独立,所以机器人模型可以看作由若干相互独立的连杆模块和关节模块串联连接而成。在Solidworks三维建模软件中,机器人整机三维模型可以拆分为若干相互独立的连杆和关节三维模型,其中关节模块如图3所示,连杆模块如图4所示,每个模块都可以看作一个超单元,并如图5所示,可以采用结构质量点单元来模拟模块之间的约束关系,质量点单元分别位于模块的输入端和输出端中心,和端面上的节点通过多点约束(MPC)连接,如图6所示,将结构质量点单元定义为主节点,所以每一个模块都可以简化成具有两个主节点的超单元,其中位于输入端的主节点编号为1,位于输出端的主节点编号为2,主节点自由度对应于超单元自由度。诸多模块对应的超单元依次串联即可组成如图2所示的机械臂等效有限元模型,如图2所示的七自由度机械臂可以等效为由15个超单元、16个主节点串联组成的有限元模型,其中E1~E15表示超单元,n1~n16表示整机有限元模型中的主节点。ki,mi分别表示单个超单元的刚度矩阵和质量矩阵,i表示超单元编号。
步骤二:采用Solidworks三维建模软件对模块三维几何模型进行简化,对模块刚度影响可以忽略的结构如倒角、圆角、螺钉孔等进行简化,并使用ANSYS有限元分析软件完成各模块建模、参数设定及单元矩阵提取,获得超单元在自身坐标系下的刚度矩阵和质量矩阵
使用ANSYS有限元分析软件完成各模块有限元的建模、参数设定以及单元矩阵的提取时,使用四面体单元划分有限元模型,结构件之间的机械接口简化为刚性连接,对模块刚度影响可以忽略的元器件及结构件均简化为结构质量点,如制动器、电机等元件,与结构件刚性连接;支撑元件如电机轴承、十字交叉滚珠轴承等简化为具有质量的六维弹簧单元;传动元器件如谐波减速器等简化为具有质量的六维弹簧单元;使用MPC将模块输入和输出端面上的各节点与结构质量点单元刚性连接,结构质量点定义为超单元主节点。以关节模型为例,如图7所示为用于提取单元刚度和质量矩阵的关节模块有限元模型。
本步骤过程具体为:
3.1、首先建立机器人各模块自身坐标系,采用Matlab数值分析软件推导各模块运动关系,在各模块的输入法兰定义输入坐标系{Oi1},输出法兰定义输出坐标系{Oi2}用于描述各模块自身及相互之间的几何特征及运动特征。
上述式(1)中,表示第i个模块输入坐标系到第i-1个模块输入坐标系的齐次变换矩阵,表示第i个模块输入坐标系到第i-1个模块输入坐标系的旋转矩阵,[xi,yi,zi]分别表示第i个模块输入坐标系原点在第i-1个模块输入坐标系X,Y,Z三个方向的位置。
3.3、推导第i(i>1)个超单元自身坐标系到整机坐标系{O0}的变换矩阵。如图2所示的模型中,各模块输入坐标系到前一模块输入坐标系的映射分别为:
当第i个超单元为关节模块时,T为常数矩阵。当第i(i>1)个超单元为连杆模块时,T为关于关节运动变量的函数矩阵,超单元自身坐标系到整机坐标系(整机坐标系与基座输入坐标系重合)的映射可以通过各超单元变换矩阵连乘获得:
3.4、根据3.3中超单元自身坐标系到整机坐标系的映射关系,获得第i个超单元到整机坐标系下的转换矩阵Si:
上述式(3)的推导可见参考文献《王新敏.ANSYS结构动力分析与应用[M].北京:人民交通出版社,2014.》。
上述式(4)和(5)中,Si称为转换矩阵,对于3D超单元,其为12×12的对角正交矩阵。
3.6、通过分块叠加原则形成整机刚度矩阵K和质量刚度矩阵M。
整机刚度矩阵K和质量矩阵M是按照单元位移编号和***整***移编号的关系进行叠加组装的。
依据图2所示机器人模型的编号叠加获得整机刚度和质量矩阵为:
本步骤通过Matlab软件实现。
步骤四:通过建立机械臂的无阻尼自由振动方程,实现协作机器人固有频率及振型的实时求解,为协作机器人实时模态分析提供依据。
本步骤具体如下:
4.1、建立协作机器人无阻尼自由振动方程,方程可以表示为:
上述式(8)根据参考文献“程丽,刘玉旺,骆海涛,等.165kg焊接机器人有限元模态分析[J].机械设计与制造,2012(1):147-149.”获得。
4.2、采用Matlab数值分析软件求解无阻尼自由振动方程。
由于协作机器人通过基座安装在工作台上,所以整机有限元模型的第一个主节点应施加固定约束,即主节点n1的节点位移为0,整机刚度矩阵和质量矩阵降阶为90×90矩阵,由于自由振动可分解为一系列简谐振动的叠加,因此式(8)的解可设为:
q=Φejωt (9);
上述式(9)中:ω为简谐振动圆频率;Φ=[φ1,φ2…φn]为节点振幅列向量;φi(i=1,2…n)为自由度i方向上的振幅。
将式(9)代入式(8)并同时消去因子ejωt,可得
(K-ω2M)×Φ=0 (10);
在进行边界条件设定后,将步骤三中的K、M代入公式(10)求解整机结构固有圆频率和模态振型,其中固有圆频率ωi(rad/s)与固有频率fi(Hz)的关系为:
ωi=2πfi (11)。
本步骤基于MATLAB数值分析软件并采用2.3GHz Core i5-6200U处理器求解上述方程处理周期约为3.9ms,相比于有限元法、振动实验法等,该方法能够大幅提高实时性,实现了协作机器人固有频率及振型的实时求解,为机器人的实时模态分析提供依据。
Claims (8)
2.根据权利要求1所述的协作机器人实时模态分析方法,其特征在于:步骤一中,将机器人整机模型拆分为多个相互独立的连杆模块和关节模块,并将每个模块等效为具有两个主节点的超单元,各模块对应的超单元依次串联组成机械臂的等效有限元模型。
4.根据权利要求1所述的协作机器人实时模态分析方法,其特征在于:步骤三中:
3.1、先建立机器人各模块自身坐标系,并推导各模块运动关系;
上述式(1)中,表示第i个模块输入坐标系到第i-1个模块输入坐标系的齐次变换矩阵,表示第i个模块输入坐标系到第i-1个模块输入坐标系的旋转矩阵,[xi,yi,zi]分别表示第i个模块输入坐标系原点在第i-1个模块输入坐标系X,Y,Z三个方向的位置;
3.3、推导第i(i>1)个超单元自身坐标系到整机坐标系的变换矩阵,超单元自身坐标系到整机坐标系(整机坐标系与基座输入坐标系重合)的映射通过各超单元变换矩阵连乘获得:
3.4、根据3.3中超单元自身坐标系到整机坐标系的映射关系,获得第i个超单元到整机坐标系下的转换矩阵Si:
3.6、通过分块叠加原则形成整机刚度矩阵K和质量刚度矩阵M:
5.根据权利要求4所述的协作机器人实时模态分析方法,其特征在于:步骤3.1中,在各模块的输入法兰定义输入坐标系{Oi1},输出法兰定义输出坐标系{Oi2}用于描述各模块自身及相互之间的几何特征及运动特征。
6.根据权利要求4所述的协作机器人实时模态分析方法,其特征在于:步骤3.3中,机器人整机模型拆分为多个相互独立的连杆模块和关节模块,当第i个超单元为关节模块时,T为常数矩阵,当第i(i>1)个超单元为连杆模块时,T为关于关节运动变量的函数矩阵。
8.根据权利要求7所述的协作机器人实时模态分析方法,其特征在于:步骤四中,由于自由振动分解为一系列简谐振动的叠加,因此所述式(8)的解设为:
q=Φejωt (9);
上述式(9)中:ω为简谐振动圆频率;Φ=[φ1,φ2…φn]为节点振幅列向量;φi(i=1,2…n)为自由度i方向上的振幅。
将式(9)代入式(8)并同时消去因子ejωt,可得:
(K-ω2M)·Φ=0 (10);
在进行边界条件设定后,将步骤三中获得的K、M代入公式(10)求解整机结构固有圆频率和模态振型,其中固有圆频率ωi(rad/s)与固有频率fi(Hz)的关系为:
ωi=2πfi (11)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911169700.2A CN112949104B (zh) | 2019-11-26 | 2019-11-26 | 一种协作机器人实时模态分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911169700.2A CN112949104B (zh) | 2019-11-26 | 2019-11-26 | 一种协作机器人实时模态分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112949104A true CN112949104A (zh) | 2021-06-11 |
CN112949104B CN112949104B (zh) | 2023-07-25 |
Family
ID=76224952
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911169700.2A Active CN112949104B (zh) | 2019-11-26 | 2019-11-26 | 一种协作机器人实时模态分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112949104B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113418591A (zh) * | 2021-07-29 | 2021-09-21 | 哈尔滨理工大学 | 五轴机床工作台轴向振动双位非接触检测装置及预测方法 |
CN113954079A (zh) * | 2021-11-23 | 2022-01-21 | 北京邮电大学 | 一种同质模块化机器人通用数学表征方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266632A (zh) * | 2008-04-18 | 2008-09-17 | 清华大学 | 一种面向产品级仿真的零部件模型降阶方法 |
CN102637228A (zh) * | 2012-04-05 | 2012-08-15 | 天津大学 | 串联机器人结构设计优化方法 |
CN104802167A (zh) * | 2015-04-17 | 2015-07-29 | 东南大学 | 基于MAPLESIM的Delta机器人快速建模方法 |
CN106295028A (zh) * | 2016-08-16 | 2017-01-04 | 西北工业大学 | 一种局部结构动力学建模方法及装置 |
CN109255188A (zh) * | 2018-09-14 | 2019-01-22 | 哈尔滨工业大学 | 基于有限元的六轴工业机器人动态性能优化方法 |
CN109522622A (zh) * | 2018-10-31 | 2019-03-26 | 中国运载火箭技术研究院 | 一种多自由度太阳电池阵在轨载荷工况确定方法及*** |
CN109992856A (zh) * | 2019-03-19 | 2019-07-09 | 东北大学 | 一种适用于不同涂层厚度的硬涂层叶盘有限元建模方法 |
CN110188417A (zh) * | 2019-05-14 | 2019-08-30 | 东南大学 | 一种基于多级超单元的复杂结构有限元模型修正方法 |
CN110414071A (zh) * | 2019-07-03 | 2019-11-05 | 上海航天设备制造总厂有限公司 | 一种重载机器人机构的参数同步优化方法 |
-
2019
- 2019-11-26 CN CN201911169700.2A patent/CN112949104B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266632A (zh) * | 2008-04-18 | 2008-09-17 | 清华大学 | 一种面向产品级仿真的零部件模型降阶方法 |
CN102637228A (zh) * | 2012-04-05 | 2012-08-15 | 天津大学 | 串联机器人结构设计优化方法 |
CN104802167A (zh) * | 2015-04-17 | 2015-07-29 | 东南大学 | 基于MAPLESIM的Delta机器人快速建模方法 |
CN106295028A (zh) * | 2016-08-16 | 2017-01-04 | 西北工业大学 | 一种局部结构动力学建模方法及装置 |
CN109255188A (zh) * | 2018-09-14 | 2019-01-22 | 哈尔滨工业大学 | 基于有限元的六轴工业机器人动态性能优化方法 |
CN109522622A (zh) * | 2018-10-31 | 2019-03-26 | 中国运载火箭技术研究院 | 一种多自由度太阳电池阵在轨载荷工况确定方法及*** |
CN109992856A (zh) * | 2019-03-19 | 2019-07-09 | 东北大学 | 一种适用于不同涂层厚度的硬涂层叶盘有限元建模方法 |
CN110188417A (zh) * | 2019-05-14 | 2019-08-30 | 东南大学 | 一种基于多级超单元的复杂结构有限元模型修正方法 |
CN110414071A (zh) * | 2019-07-03 | 2019-11-05 | 上海航天设备制造总厂有限公司 | 一种重载机器人机构的参数同步优化方法 |
Non-Patent Citations (1)
Title |
---|
李志荣;许海;孙春华;: "兰杰文振子的结构动力学有限元模型", 苏州市职业大学学报, no. 01 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113418591A (zh) * | 2021-07-29 | 2021-09-21 | 哈尔滨理工大学 | 五轴机床工作台轴向振动双位非接触检测装置及预测方法 |
CN113418591B (zh) * | 2021-07-29 | 2024-01-26 | 哈尔滨理工大学 | 五轴机床工作台轴向振动双位非接触检测装置及预测方法 |
CN113954079A (zh) * | 2021-11-23 | 2022-01-21 | 北京邮电大学 | 一种同质模块化机器人通用数学表征方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112949104B (zh) | 2023-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dong et al. | A screw theory-based semi-analytical approach for elastodynamics of the tricept robot | |
CN108015763A (zh) | 一种抗噪声干扰的冗余度机械臂路径规划方法 | |
CN112949103A (zh) | 一种协作机器人非线性刚度建模方法 | |
CN112949104B (zh) | 一种协作机器人实时模态分析方法 | |
Ma et al. | Static and dynamic performance evaluation of a 3-DOF spindle head using CAD–CAE integration methodology | |
CN110826132B (zh) | 一种结构分散振动控制***设计方法 | |
CN105938329A (zh) | 一种数控机床动态特性仿真及其与数控***集成的方法 | |
CN112861388B (zh) | 一种基于正交设计的协作机器人全域结构优化设计方法 | |
Zhang et al. | Mathematical modeling and dynamic characteristic analysis of a novel parallel tracking mechanism for inter-satellite link antenna | |
Liao et al. | The design of flexible robotic manipulators with optimal arm geometries fabricated from composite laminates with optimal material properties | |
CN103034759B (zh) | 基于参数模块化加工中心刚度计算的机床配置方法 | |
CN115407680A (zh) | 一种并行化的高精度时域气动弹性控制耦合仿真方法 | |
Lamberson et al. | Aeroelastic simulations with modal and finite-element structural solvers using CREATE-AV/Kestrel v5 | |
CN116127669A (zh) | 一种基于超参数单元的管路***结构评估方法 | |
CN110334459A (zh) | 一种输电塔塔线体系快速精细化建模***及方法 | |
CN111027125B (zh) | 一种蜂窝等效单元及基于该蜂窝等效单元的参数计算方法 | |
Xiaodong et al. | The high efficient dynamics modeling method for modular manipulator based on Space Operator Algebra | |
CN111238755A (zh) | 一种多层框架子结构远程协同试验方法及*** | |
CN113910211B (zh) | 考虑自重的协作机器***性变形建模方法及终端 | |
Ma et al. | Dynamic modeling and analysis of the 3-P RS power head based on the screw theory and rigid multipoint constraints | |
Chen et al. | Research on Mechanical Environment Prediction Method for Structure Considering Random Uncertain Elastic Joint | |
CN117610380B (zh) | 基于正交实验联合拓扑优化的机器人轻量化设计方法 | |
CN116062198B (zh) | 面向超大型航天结构的虚实融合地面试验***及方法 | |
CN114492000A (zh) | 基于Modelica的足式机器人建模***、方法、设备及介质 | |
CN117332636A (zh) | 一种基于apdl命令流的复杂结构有限元模态分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |