CN109885968A - 一种复杂运动的运动描述建模方法及系统 - Google Patents
一种复杂运动的运动描述建模方法及系统 Download PDFInfo
- Publication number
- CN109885968A CN109885968A CN201910187086.6A CN201910187086A CN109885968A CN 109885968 A CN109885968 A CN 109885968A CN 201910187086 A CN201910187086 A CN 201910187086A CN 109885968 A CN109885968 A CN 109885968A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- sports immunology
- movement
- immunology
- sports
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开一种复杂运动的运动描述建模方法及系统，所述方法包括：获取目标的实际运动；对所述实际运动进行分解，得到多级运动；根据各级运动获取各级运动的参考系和运动描述方程；根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；根据所述两组坐标系建立装配关系；根据所述装配关系增加运动描述；将所述装配关系和所述运动描述合成，得到最终的运动模型，用于CFD数值仿真得到仿真对象的流场分布及气动作用力。本发明中的上述方法使用多级坐标系变换，实现复合形式的运动，或可对观察到的运动进行级数拟合，来描述复杂运动。
Description
技术领域
本发明涉及流体动力学领域，特别是涉及一种复杂运动的运动描述建模方法及系统。
背景技术
计算流体动力学，(Computational Fluid Dynamics，CFD)是近代流体力学，数值数学和计算机科学结合的产物，是一门具有强大生命力的交叉科学。它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题，在模拟中往往伴随着物体运动，且难以仅通过参考系变换转变为静止状态，如直升机桨叶旋转等，此时就需要添加合适的运动描述来实现这一运动。
针对上述运动问题，现有方案并未很好解决。绝大部分软件仅支持如单方向平移，旋转。其主要问题在于实际物体运动复杂多变，难以用公式进行充分描述。于是，部分商业CFD软件直接将该问题交由用户自行处理，如Fluent中依赖UDF，由用户直接编写代码实现运动，这虽然满足需求，但也给软件的使用带来了极大的不便。
由用户自行编写代码实现运动，这除了需要用户能清楚每个运动细节外，还要求用户具有较强的编程能力。而这些，都给CFD模拟带来了而外的工作，且针对设计分析及优化中细微的调整或改变需要重新编写代码，不利于快速改动。于是，建立一套适用于CFD模拟的运动描述，使得便于用户实现复杂运动建模及编写相应配置文件。
发明内容
本发明的目的是提供一种复杂运动的运动描述建模方法及系统，来模拟复杂运动。
为实现上述目的，本发明提供了如下方案：
一种复杂运动的运动描述建模方法，所述方法包括：
获取目标的实际运动；
对所述实际运动进行分解，得到多级运动；
根据各级运动获取各级运动的参考系和运动描述方程；
根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
根据所述两组坐标系建立装配关系；
根据所述装配关系增加运动描述；
将所述装配关系和所述运动描述合成，得到最终的运动模型。
可选的，所述根据所述两组坐标系建立装配关系具体包括：
固定参考点，所述参考点用b表示；
通过所述参考点b建立第一坐标系O1和第二坐标系O2，且所述第一坐标系和所述第二坐标系中的基向量在所述第一坐标系和所述第二坐标系中的表示为第一基向量第二基向量
所述参考点b与所述第一基向量以及所述第二基向量在实际空间位置中重合，构成两组新坐标系，记为Ob1和Ob2； 至此，装配得到的坐标系间的坐标变换已经确定。
可选的，所述根据所述装配关系增加运动描述具体包括：
定义坐标系Ob2，设坐标系Ob2由0时刻与Ob1重合变成了时刻t下的O'b2；
定义函数向量为时刻t时坐标系O'b2的原点在坐标Ob1中表示；
定义函数矩阵为t时刻坐标系O'b2的基向量组在坐标系Ob2的基向量组下的表示，且满足|Rsi(t),Rsj(t),Rsk(t)|＝1，|Ris(t),Rjs(t),Rks(t)|＝1，其中s＝i,j,k，初始位置重合；
建立运动后的坐标系O'b1，b'1＝b1+T(t)；
可选的，所述将所述装配关系和所述运动描述合成，得到最终的运动模型具体包括：
将运动坐标系下的点P2变换成固定坐标系下的点P'1；
记
P'1＝A(t)×(p1+b)+c(t)。
本发明还另外提供一种复杂运动的运动描述建模系统，所述系统包括：
实际运动获取模块，用于获取目标的实际运动；
分解模块，用于对所述实际运动进行分解，得到多级运动；
参考系和描述方程获取模块，用于根据各级运动获取各级运动的参考系和运动描述方程；
坐标系建立模块，用于根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
装配关系建立模块，用于根据所述两组坐标系建立装配关系；
运动描述建立模块，用于根据所述装配关系增加运动描述；
运动模型确定模块，用于将所述装配关系和所述运动描述合成，得到最终的运动模型；
运动结果接口模块，用于将最终的运动模型的运动结果导入至CFD数值仿真中使用，以便计算得到流场分布及气动作用力。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明中的上述方法适用于CFD数值仿真中复杂刚体运动的嵌套形式，使用多级坐标系变换，实现复合形式的运动，或可对观察到的运动进行级数拟合，来描述复杂运动；且该运动描述在课题组关于直升机CFD数值仿真中被大量应用，验证了有效性；使用多级坐标变换，将复杂运动转换为简单基本运动，描述简单且精确，并结合解释型描述函数，指定运动方式更简单；支持对任意运动进行所需精度的级数拟合，表达准确，仅需依次计算各级变换间的数值解，且最终表达式使用矩阵运算实现，便于CPU/GPU矢量运算。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例复杂运动的运动描述建模方法流程图；
图2为本发明实施例复杂运动的运动描述建模系统结构示意图；
图3为本发明实施例CFD计算具有周期变距旋翼运动效果图；
图4为本发明实施例CFD计算具有主动襟翼控制旋翼运动效果图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种复杂运动的运动描述建模方法及系统，来模拟复杂运动。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例复杂运动的运动描述建模方法流程图，如图1所示，所述方法包括：
步骤101：获取目标的实际运动；
步骤102：对所述实际运动进行分解，得到多级运动；
步骤103：根据各级运动获取各级运动的参考系和运动描述方程；
步骤104：根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
步骤105：根据所述两组坐标系建立装配关系；
步骤106：根据所述装配关系增加运动描述；
步骤107：将所述装配关系和所述运动描述合成，得到最终的运动模型。
具体的，步骤101中，实际运动是指CFD数值模拟的分析对象的运动，例如普通螺旋桨叶运动。
具体的，步骤102和步骤103中，针对CFD数值仿真中的运动进行分析，将实际运动分解为尽可能简单的多级运动，并获得各级运动下参考系和相对运动描述方程，对于不便直接给出运动方程的实验观测结果，可采用函数拟合形式。
具体的，步骤105中，各级运动坐标系之间的关系，每两两间运动，以其中一个为参照，另一个对其存在相对运动，建立两组坐标系，确定坐标系间的变化，记为建立装配关系，在装配中需要先固定参考点，所述参考点用b表示；在空间中取任意点，在第一坐标系O1下的表示为在第二坐标系O2下的表示为再固定参考点b的位置，最后固定参考点b的参考方向，即通过所述参考点b建立第一坐标系O1和第二坐标系O2，且所述第一坐标系和所述第二坐标系中的基向量在所述第一坐标系和所述第二坐标系中的表示为第一基向量第二基向量
所述参考点b与所述第一基向量组以及所述第二基向量组在实际空间位置中重合，构成两组新坐标系，记为Ob1和Ob2，那么空间中任意点，分别用坐标系Ob1和坐标系Ob2的表示应相同，于是有其中P1和P2分别为点P在坐标系O1和O2中的表示，即已知可求出P1为：上式中， 至此，装配得到的坐标系间的坐标变换已经确定。
具体的，步骤106中，在装配完成的基础上，可增加运动描述，运动描述为运动坐标系对参照系之间，坐标系变换矩阵值随时间变化的函数表示，设坐标系Ob2由0时刻与Ob1重合变成了时刻t下的O'b2
定义函数向量为时刻t时坐标系O'b2的原点在坐标Ob1中表示；
定义函数矩阵为t时刻坐标系O'b2的基向量组在坐标系Ob2的基向量组下的表示，且满足|Rsi(t),Rsj(t),Rsk(t)|＝1，|Ris(t),Rjs(t),Rks(t)|＝1，其中s＝i,j,k，即初始位置重合，函数矩阵中行或列构成基向量组；
再建立运动后的坐标系O'b1，b'1＝b1+T(t)；
同理，仍有满足
运动后的坐标系仍然满足P2与坐标系Ob2相对位置关系不变，则有
P'2＝P2，于是解得：
具体的，步骤107中，将运动坐标系下的点P2变换成固定坐标系下的点P'1；
记
P'1＝A(t)×(p1+b)+c(t)。
将上述得到的运动参数及关系编写成运动配置文件，包含如下结构：
运动方程：上述步骤建模步骤中运动描述方程的公式书写表示；
坐标系表示：
参考坐标系：表示该坐标系依赖的坐标系，缺省值为参考地坐标系；
参考点：在参考系中选取的装配点(可视情况任取，便于设置本地点即可)，记为Pref；
本地点：本地坐标系中和参考点实际位置重合的点的本地坐标，记为Plocal；
参考轴：在参考系中选取的基准轴系(可视情况任取，便于设置本地轴即可)，记为Axisref；
本地轴：本地坐标系和参考点实际位置重合的点，记为Axislocal；
运动方程索引：根据索引从运动方程中检索，包含旋转运动方程T(t)和平移运动方程x(t)两部分；
装配关联：指定模型关联到哪个运动坐标系上；
使用程序读取并解释运动配置文件，并由用户输入所需的计算时刻；
对各模型，根据其坐标系表示计算运动后的模型相对坐标点，有原始坐标Pold得到新坐标的计算公式为：
Pnew＝(T(t)*Axislocal)-1*Axislocal*(Pold-Plocal)+Pref+x(t)；
若坐标系存在参考坐标系，将上步骤得到的Pnew作为下次执行的Pold跳转到上一步；
输出最终运动后的模型坐标点并结束程序。
图2为本发明实施例复杂运动的运动描述建模系统结构示意图，如图2所示，所述系统包括：
实际运动获取模块201，用于获取目标的实际运动；
分解模块202，用于对所述实际运动进行分解，得到多级运动；
参考系和描述方程获取模块203，用于根据各级运动获取各级运动的参考系和运动描述方程；
坐标系建立模块204，用于根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
装配关系建立模块205，用于根据所述两组坐标系建立装配关系；
运动描述建立模块206，用于根据所述装配关系增加运动描述；
运动模型确定模块207，用于将所述装配关系和所述运动描述合成，得到最终的运动模型；
运动结果接口模块208，用于将最终的运动模型的运动结果导入至CFD数值仿真中使用，以便计算得到流场分布及气动作用力。
图3、图4分别为本发明实施例运动效果图，按照图2所述流程，将首先将运动进行分解，其次构建配置文件，再次编写程序，并将配置文件作为输入，最后将输出结果用于CFD数值仿真。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。
Claims (5)
1.一种复杂运动的运动描述建模方法，其特征在于，所述方法包括：
获取目标的实际运动；
对所述实际运动进行分解，得到多级运动；
根据各级运动获取各级运动的参考系和运动描述方程；
根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
根据所述两组坐标系建立装配关系；
根据所述装配关系增加运动描述；
将所述装配关系和所述运动描述合成，得到最终的运动模型。
2.根据权利要求1所述的复杂运动的运动描述建模方法，其特征在于，所述根据所述两组坐标系建立装配关系具体包括：
固定参考点，所述参考点用b表示；
通过所述参考点b建立第一坐标系O1和第二坐标系O2，且所述第一坐标系和所述第二坐标系中的基向量在所述第一坐标系和所述第二坐标系中的表示为第一基向量第二基向量
所述参考点b与所述第一基向量以及所述第二基向量在实际空间位置中重合，构成两组新坐标系，记为Ob1和Ob2； 至此，装配得到的坐标系间的坐标变换已经确定。
3.根据权利要求1所述的复杂运动的运动描述建模方法，其特征在于，所述根据所述装配关系增加运动描述具体包括：
定义坐标系Ob2，设坐标系Ob2由0时刻与Ob1重合变成了时刻t下的O'b2；
定义函数向量为时刻t时坐标系O'b2的原点在坐标Ob1中表示；
定义函数矩阵为t时刻坐标系O'b2的基向量组在坐标系Ob2的基向量组下的表示，且满足|Rsi(t),Rsj(t),Rsk(t)|＝1，|Ris(t),Rjs(t),Rks(t)|＝1，其中s＝i,j,k，初始位置重合；
建立运动后的坐标系O'b1，b'1＝b1+T(t)；
4.根据权利要求1所述的复杂运动的运动描述建模方法，其特征在于，所述将所述装配关系和所述运动描述合成，得到最终的运动模型具体包括：
将运动坐标系下的点P2变换成固定坐标系下的点P'1；
记
P'1＝A(t)×(p1+b)+c(t)。
5.一种复杂运动的运动描述建模系统，其特征在于，所述系统包括：
实际运动获取模块，用于获取目标的实际运动；
分解模块，用于对所述实际运动进行分解，得到多级运动；
参考系和描述方程获取模块，用于根据各级运动获取各级运动的参考系和运动描述方程；
坐标系建立模块，用于根据所述参考系和所述描述方程，每两两间运动，以其中一个运动为参考建立两组坐标系；
装配关系建立模块，用于根据所述两组坐标系建立装配关系；
运动描述建立模块，用于根据所述装配关系增加运动描述；
运动模型确定模块，用于将所述装配关系和所述运动描述合成，得到最终的运动模型；
运动结果接口模块，用于将最终的运动模型的运动结果导入至CFD数值仿真中使用，以便计算得到流场分布及气动作用力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910187086.6A CN109885968B (zh) | 2019-03-13 | 2019-03-13 | 一种复杂运动的运动描述建模方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910187086.6A CN109885968B (zh) | 2019-03-13 | 2019-03-13 | 一种复杂运动的运动描述建模方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109885968A true CN109885968A (zh) | 2019-06-14 |
| CN109885968B CN109885968B (zh) | 2020-11-20 |
Family
ID=66931946
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910187086.6A Active CN109885968B (zh) | 2019-03-13 | 2019-03-13 | 一种复杂运动的运动描述建模方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109885968B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050262933A1 (en) * | 2004-05-26 | 2005-12-01 | Fang Zhu | Tire uniformity through compensation between radial run out and stiffness variation |
| CN102161153A (zh) * | 2011-02-28 | 2011-08-24 | 西安交通大学 | 一种用于自动化装配的模块化柔性六自由度并联冗余驱动调姿机构及其调整方法 |
| CN106055901A (zh) * | 2016-06-02 | 2016-10-26 | 西北工业大学 | 一种自由漂浮空间机器人捕获翻滚目标的时机确定方法 |
| CN106485035A (zh) * | 2016-11-23 | 2017-03-08 | 长春理工大学 | 基于随机有限元分析的红外诱饵弹气动特性建模方法 |
| CN108920789A (zh) * | 2018-06-20 | 2018-11-30 | 北京控制工程研究所 | 一种柔性航天器多级复合控制动力学建模方法 |
| CN108917753A (zh) * | 2018-04-08 | 2018-11-30 | 中国人民解放军63920部队 | 基于从运动恢复结构的飞行器位置确定方法 |
| CN109033600A (zh) * | 2018-07-18 | 2018-12-18 | 安徽大学 | 一种单侧人型杆缠绕过程有限元仿真分析方法 |
-
2019
- 2019-03-13 CN CN201910187086.6A patent/CN109885968B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050262933A1 (en) * | 2004-05-26 | 2005-12-01 | Fang Zhu | Tire uniformity through compensation between radial run out and stiffness variation |
| CN102161153A (zh) * | 2011-02-28 | 2011-08-24 | 西安交通大学 | 一种用于自动化装配的模块化柔性六自由度并联冗余驱动调姿机构及其调整方法 |
| CN106055901A (zh) * | 2016-06-02 | 2016-10-26 | 西北工业大学 | 一种自由漂浮空间机器人捕获翻滚目标的时机确定方法 |
| CN106485035A (zh) * | 2016-11-23 | 2017-03-08 | 长春理工大学 | 基于随机有限元分析的红外诱饵弹气动特性建模方法 |
| CN108917753A (zh) * | 2018-04-08 | 2018-11-30 | 中国人民解放军63920部队 | 基于从运动恢复结构的飞行器位置确定方法 |
| CN108920789A (zh) * | 2018-06-20 | 2018-11-30 | 北京控制工程研究所 | 一种柔性航天器多级复合控制动力学建模方法 |
| CN109033600A (zh) * | 2018-07-18 | 2018-12-18 | 安徽大学 | 一种单侧人型杆缠绕过程有限元仿真分析方法 |
Non-Patent Citations (3)
| Title |
|---|
| 叶靓 等: ""共轴式双旋翼悬停流场和气动力的CFD计算"", 《空气动力学学报》 * |
| 胡志远 等: ""基于CFD方法的主动襟翼控制旋翼翼型涡特性研究"", 《南京航空航天大学学报》 * |
| 雷静桃 等: ""基于单元划分和连接关系矩阵的模块化并联机器人运动学建模"", 《高技术通讯》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109885968B (zh) | 2020-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Real-time finite element structural analysis in augmented reality | |
| Terze et al. | Lie-group integration method for constrained multibody systems in state space | |
| CN108015763A (zh) | 一种抗噪声干扰的冗余度机械臂路径规划方法 | |
| CN104200052B (zh) | 液压驱动六足机器人仿真系统及仿真方法 | |
| Terze et al. | Singularity-free time integration of rotational quaternions using non-redundant ordinary differential equations | |
| EP1966687A2 (en) | Digital effects analysis in modeling environments | |
| CN104123411B (zh) | 一种汽车动力总成系统转动惯量合成的简洁推导方法 | |
| Blajer | On the determination of joint reactions in multibody mechanisms | |
| Antonya et al. | Design evaluation and modification of mechanical systems in virtual environments | |
| Hoffmann et al. | On optical data-guided optimal control simulations of human motion | |
| CN106202735A (zh) | 一种局部非线性地基土‑结构相互作用子结构试验方法 | |
| CN111554355B (zh) | 一种基于非冯诺依曼架构的分子动力学计算方法 | |
| CN109885968A (zh) | 一种复杂运动的运动描述建模方法及系统 | |
| JP4612989B2 (ja) | メカトロニクスシステムのシミュレーション方法 | |
| Duan et al. | Dynamic modeling and experimental research on position-dependent behavior of twin ball screw feed system | |
| CN115407680A (zh) | 一种并行化的高精度时域气动弹性控制耦合仿真方法 | |
| Wei et al. | Design and control of a multi-axis servo motion chair system based on a microcontroller | |
| Hu et al. | Inverse kinematics model’s parameter simulation for stewart platform design of driving simulator | |
| Ge et al. | Blending isogeometric and Lagrangian elements in three-dimensional analysis | |
| Cammarata et al. | A system-based reduction method for spatial deformable multibody systems using global flexible modes | |
| Otani | Software tools for dynamic and kinematic modeling of human motion | |
| Hoshino | Interpolation and extrapolation of repeated motions obtained with magnetic motion capture | |
| Yang et al. | Research on Longitudinal Control and Visual Simulation System for Civil Aircraft Based on Simulink/FlightGear | |
| Shi et al. | Local Extrema Refinement Based Tensor Product Model Transformation Controller Design with Vary Input Methods | |
| Binante et al. | NOSA-ITACA 1.1 documentation |
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 |