CN113297730A - 基于自适应模态的柔性多体系统动态响应计算方法和系统 - Google Patents

基于自适应模态的柔性多体系统动态响应计算方法和系统 Download PDF

Info

Publication number
CN113297730A
CN113297730A CN202110523265.XA CN202110523265A CN113297730A CN 113297730 A CN113297730 A CN 113297730A CN 202110523265 A CN202110523265 A CN 202110523265A CN 113297730 A CN113297730 A CN 113297730A
Authority
CN
China
Prior art keywords
body system
flexible multi
flexible
dynamic response
representing
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
Application number
CN202110523265.XA
Other languages
English (en)
Other versions
CN113297730B (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202110523265.XA priority Critical patent/CN113297730B/zh
Publication of CN113297730A publication Critical patent/CN113297730A/zh
Application granted granted Critical
Publication of CN113297730B publication Critical patent/CN113297730B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开一种基于自适应模态的柔性多体系统动态响应计算方法，包括以下步骤：S1：将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；S2：设计柔性多体系统的运动方程的抽样矩阵；S3：求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；S4：利用贪婪高斯‑牛顿算法求解l1范数优化问题，得到柔性多体系统动态响应。本发明定义了可以用模态形状稀疏表示的弹性坐标，导出了具有少量运动方程约束的宽范数坐标下的l1范数优化问题，作为每个时间步长的动态响应分析。提出了GGN算法，有效地解决了这种优化问题，使模态坐标能够自适应地选择各种工作模式。

Description

基于自适应模态的柔性多体系统动态响应计算方法和系统
技术领域
本发明涉及柔性多体系统技术领域，更具体地，涉及一种基于自适应模态的柔性多体系统动态响应计算方法和系统。
背景技术
柔性多体系统(FMS)是由刚性和变形构件以某种方式连接而成的系统。利用有限元法将柔性体离散为有限自由度。对于复杂的组件，利用有限元法将产生数千个自由度。这将给柔性多体动态响应分析的计算效率带来负担。因此，需要采用模型降阶技术来降低自由度。模态截断是一种常用的约简技术。它通过选取低频模态形状将弹性坐标转换为模态坐标。应用此方法后，柔性部件的自由度将显著降低。然而，模态截断主要缺点是不能直接考虑执行器的控制系统或位置，对许多低频模态集对结构和控制系统之间的相互作用没有贡献。
公开日为2020年05月15日，公开号为CN111159636A的中国专利公开了一种基于绝对节点坐标描述的柔性多体系统动力学半解析灵敏度分析方法，首先，基于绝对节点坐标方法，建立柔性多体系统的质量矩阵、刚度矩阵和广义力列阵；其次，建立柔性多体系统的动力学方程和优化目标函数；再次，基于直接微分法或伴随变量法，建立柔性多体系统动力学的半解析灵敏度计算公式；最后，求解柔性多体系统动力学微分代数方程，获得灵敏度计算结果。该专利的计算方法效率较低。
发明内容
本发明的首要目的是提供一种基于自适应模态的柔性多体系统动态响应计算方法，提高柔性多体动态响应分析的计算效率且突出适应不同的工况。
本发明的次要目的是提供一种基于自适应模态的柔性多体系统动态响应计算系统。
为解决上述技术问题，本发明的技术方案如下：
一种基于自适应模态的柔性多体系统动态响应计算方法，包括以下步骤：
S1：将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
S2：设计柔性多体系统的运动方程的抽样矩阵；
S3：求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
S4：利用贪婪高斯-牛顿算法求解l1范数优化问题，得到柔性多体系统动态响应。
优选地，步骤S1中将柔性多体系统的模态坐标定义为正交振型上的稀疏系数，具体为：
设模态坐标是稀疏的或近似稀疏的，柔性多体系统中有N个柔性体，则广义坐标q表示如下：
Figure BDA0003064799310000021
式中，i表示第i个柔性体柔性体，
Figure BDA0003064799310000022
Figure BDA00030647993100000211
为柔性体i的弹性坐标的向量，RiT为柔性体i相对于惯性系的位置向量，θiT为柔性体i相对于惯性系的角位移，I是单位矩阵，Bi是由柔性体i的全模态形状组成，
Figure BDA0003064799310000023
Figure BDA0003064799310000024
其中
Figure BDA0003064799310000025
表示柔性多体系统的弹性坐标，
Figure BDA0003064799310000026
表示柔性多体系统的模态坐标；
结合拉格朗日乘子向量λ，可得：
Figure BDA0003064799310000027
式中，
Figure BDA0003064799310000028
表示信号，
Figure BDA0003064799310000029
表示基，
Figure BDA00030647993100000210
表示系数向量。
优选地，步骤S2中柔性多体系统的运动方程的抽样矩阵，具体为：
Figure BDA0003064799310000031
式中，
Figure BDA0003064799310000032
表示相对于柔性体的位置和旋转的质量矩阵，
Figure BDA0003064799310000033
是耦合刚性运动和变形的质量矩阵，
Figure BDA0003064799310000034
表示单个柔性体与弹性坐标有关的质量矩阵，
Figure BDA0003064799310000035
表示与单个柔性体的弹性坐标有关的对称刚度矩阵，
Figure BDA0003064799310000036
和
Figure BDA0003064799310000037
分别是相对于刚性和弹性坐标的外力，
Figure BDA0003064799310000038
和
Figure BDA0003064799310000039
分别是关于刚性和弹性坐标的二次速度矢量，其中Gr和Gf表示与刚性和弹性坐标相关的g(q)的Jacobin矩阵，Ωi表示剩余方程的指数，I(Ωi,:)表示提取单元矩阵的Ωi行。
优选地，步骤S2中柔性多体系统的运动方程，简化写为：
Figure BDA00030647993100000310
式中，(·)Ω表示柔性体ith采样后的矩阵或向量，
Figure BDA00030647993100000311
表示柔性多体系统的质量矩阵，
Figure BDA00030647993100000312
表示柔性多体系统的刚度矩阵，
Figure BDA00030647993100000313
表示柔性多体系统的约束雅可比矩阵转置形式，
Figure BDA00030647993100000314
表示柔性多体系统的广义外力，
Figure BDA00030647993100000315
表示柔性多体系统的二次速度矢量。
优选地，步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
优选地，步骤S3中每个时间步长下柔性多体系统的动态响应，非线性映射算子定义为：
Figure BDA00030647993100000316
式中，Φ∈Rm×n是一个测量矩阵，Φ(·)∈Rm×l表示非线性映射算子，
Figure BDA0003064799310000041
表示t时刻的弹性坐标和模态坐标，h表示步长，M(pt)表示t时刻质量矩阵，(·)Ωt表示t时刻柔性多体系统采样后的矩阵或向量，f(pt)表示t时刻作用于柔性体的广义外力，GT(pt)表示t时刻约束雅可比矩阵的转置，
Figure BDA0003064799310000042
表示t时刻拉格朗日乘子，g(pt)表示t时刻约束方程，且
Figure BDA0003064799310000043
a表示系数向量。
优选地，步骤S3中每个时间步长动态响应的l1范数优化问题，定义为：
Figure BDA00030647993100000413
约束
Figure BDA0003064799310000044
式中，
Figure BDA00030647993100000414
表示t时刻柔性多体系统的模态坐标，||·||1表示1范数，s表示柔性多体系统模态坐标的稀疏度，则
Figure BDA00030647993100000415
的稀疏性可表示为:
Figure BDA0003064799310000045
式中，si表示柔性体i模态坐标的稀疏度。
优选地，步骤S4中利用GGN算法求解l1范数优化问题，具体为：
S4.1：读取每个柔性体的有限元信息，计算每个柔性体的惯性形状积分模态，设置初始参数，所述惯性形状积分模态包括mff和Kff；
S4.2：求解
Figure BDA0003064799310000046
得到初始时刻的加速度
Figure BDA0003064799310000047
和拉格朗日乘子λ0，式中，M表示柔性多体系统的质量矩阵，G表示柔性多体系统的约束雅可比矩阵，
Figure BDA0003064799310000048
表示柔性多体系统的加速度，λ表示柔性多体系统的拉格朗日乘子，
Figure BDA0003064799310000049
表示作用于柔性多体系统的广义外力，
Figure BDA00030647993100000410
表示与加速度有关的向量；
S4.3：对每个柔性体进行测量，得到其待定运动方程Φ(c)；
S4.4：利用GGN算法求解
Figure BDA00030647993100000411
得到p1,
Figure BDA00030647993100000412
表示在t时刻的坐标，速度，加速度和拉格朗日乘子；
S4.5：判断t是否大于Time，若是，则结束计算，若不是，令t＝t+h，返回步骤S4.3。
优选地，步骤S4.1中所述初始参数包括p0,
Figure BDA0003064799310000051
Time,h，其中p0表示初始时刻的坐标，
Figure BDA0003064799310000052
表示初始时刻的速度，Time表示模拟时间，h表示时间步长。
一种基于自适应模态的柔性多体系统动态响应计算系统，所述系统基于上述所述的基于自适应模态的柔性多体系统动态响应计算方法，包括：
模态模块，所述模态模块将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
抽样模块，所述抽样模块设计柔性多体系统的运动方程的抽样矩阵；
范数优化模块，所述范数优化模块求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
GGN求解模块，所述GGN求解模块利用GGN算法求解l1范数优化问题，得到柔性多体系统动态响应。
与现有技术相比，本发明技术方案的有益效果是：
本发明定义了可以用模态形状稀疏表示的弹性坐标，导出了具有少量运动方程约束的宽范数坐标下的l1范数优化问题，作为每个时间步长的动态响应分析。提出了GGN算法，有效地解决了这种优化问题，使模态坐标能够自适应地选择各种工作模式。
附图说明
图1为本发明的方法流程示意图。
图2为实施例1中GGN算法流程图。
图3为实施例1中曲柄滑块机构示意图。
图4为实施例1中曲柄滑块机构点C的x位移图。
图5为使用实施例1所述的方法与FOM之间的相对误差。
图6为使用实施例1所述的方法计算模态坐标的值和选择示意图。
图7为实施例2的系统示意图。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；
对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
本实施例提供一种基于自适应模态的柔性多体系统动态响应计算方法，如图1，包括以下步骤：
S1：将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
S2：设计柔性多体系统的运动方程的抽样矩阵；
S3：求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
S4：利用贪婪高斯-牛顿算法求解l1范数优化问题，得到柔性多体系统动态响应。
步骤S1中将柔性多体系统的模态坐标定义为正交振型上的稀疏系数，具体为：
设模态坐标是稀疏的或近似稀疏的，柔性多体系统中有N个柔性体，则广义坐标q表示如下：
Figure BDA0003064799310000061
式中，i表示第i个柔性体柔性体，
Figure BDA0003064799310000062
为柔性体i的弹性坐标的向量，RiT为柔性体i相对于惯性系的位置向量，θiT为柔性体i相对于惯性系的角位移，I是单位矩阵，Bi是由柔性体i的全模态形状组成，
Figure BDA0003064799310000063
Figure BDA0003064799310000064
其中
Figure BDA0003064799310000065
表示柔性多体系统的弹性坐标，
Figure BDA0003064799310000066
表示柔性多体系统的模态坐标；
结合拉格朗日乘子向量λ，可得：
Figure BDA0003064799310000071
式中，
Figure BDA0003064799310000072
表示信号，
Figure BDA0003064799310000073
表示基，
Figure BDA0003064799310000074
表示系数向量。
步骤S2中柔性多体系统的运动方程的抽样矩阵，具体为：
Figure BDA0003064799310000075
式中，
Figure BDA0003064799310000076
表示相对于柔性体的位置和旋转的质量矩阵，
Figure BDA0003064799310000077
是耦合刚性运动和变形的质量矩阵，
Figure BDA0003064799310000078
表示单个柔性体与弹性坐标有关的质量矩阵，
Figure BDA0003064799310000079
表示与单个柔性体的弹性坐标有关的对称刚度矩阵，
Figure BDA00030647993100000710
和
Figure BDA00030647993100000711
分别是相对于刚性和弹性坐标的外力，
Figure BDA00030647993100000712
和
Figure BDA00030647993100000713
分别是关于刚性和弹性坐标的二次速度矢量，其中Gr和Gf表示与刚性和弹性坐标相关的g(q)的Jacobin矩阵，Ωi表示剩余方程的指数，I(Ωi,:)表示提取单元矩阵的Ωi行。
步骤S2中柔性多体系统的运动方程，简化写为：
Figure BDA00030647993100000714
式中，(·)Ω表示柔性体ith采样后的矩阵或向量，
Figure BDA00030647993100000715
表示柔性多体系统的质量矩阵，
Figure BDA00030647993100000716
表示柔性多体系统的刚度矩阵，
Figure BDA00030647993100000717
表示柔性多体系统的约束雅可比矩阵转置形式，
Figure BDA00030647993100000718
表示柔性多体系统的广义外力，
Figure BDA00030647993100000719
表示柔性多体系统的二次速度矢量。
步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
步骤S3中每个时间步长下柔性多体系统的动态响应，非线性映射算子定义为：
Figure BDA0003064799310000081
式中，Φ∈Rm×n是一个测量矩阵，Φ(·)∈Rm×l表示非线性映射算子，
Figure BDA0003064799310000082
表示t时刻的弹性坐标和模态坐标，h表示步长，M(pt)表示t时刻质量矩阵，(·)Ωt表示t时刻柔性多体系统采样后的矩阵或向量，f(pt)表示t时刻作用于柔性体的广义外力，GT(pt)表示t时刻约束雅可比矩阵的转置，
Figure BDA0003064799310000083
表示t时刻拉格朗日乘子，g(pt)表示t时刻约束方程，且
Figure BDA0003064799310000084
a表示系数向量。
步骤S3中每个时间步长动态响应的l1范数优化问题，定义为：
Figure BDA00030647993100000812
约束
Figure BDA0003064799310000085
式中，
Figure BDA00030647993100000813
表示t时刻柔性多体系统的模态坐标，||·||1表示1范数，s表示柔性多体系统模态坐标的稀疏度，则
Figure BDA00030647993100000814
的稀疏性可表示为:
Figure BDA0003064799310000086
式中，si表示柔性体i模态坐标的稀疏度。
步骤S4中利用GGN算法求解l1范数优化问题，具体为：
S4.1：读取每个柔性体的有限元信息，计算每个柔性体的惯性形状积分模态，设置初始参数，所述惯性形状积分模态包括mff和Kff；
S4.2：求解
Figure BDA0003064799310000087
得到初始时刻的加速度
Figure BDA0003064799310000088
和拉格朗日乘子λ0，式中，M表示柔性多体系统的质量矩阵，G表示柔性多体系统的约束雅可比矩阵，
Figure BDA0003064799310000089
表示柔性多体系统的加速度，λ表示柔性多体系统的拉格朗日乘子，
Figure BDA00030647993100000810
表示作用于柔性多体系统的广义外力，
Figure BDA00030647993100000811
表示与加速度有关的向量；
S4.3：对每个柔性体进行测量，得到其待定运动方程Φ(c)；
S4.4：利用GGN算法求解
Figure BDA0003064799310000098
得到p1,
Figure BDA0003064799310000091
表示在t时刻的坐标，速度，加速度和拉格朗日乘子，所述利用GGN算法求解的具体过程如图2所示，在算法中，符号Sk表示支持集，包括刚性坐标和拉格朗日乘子的索引，和
Figure BDA0003064799310000092
Figure BDA0003064799310000099
表示提取
Figure BDA0003064799310000093
的tth列，
Figure BDA0003064799310000094
由
Figure BDA0003064799310000095
列组成，对应于支持集Sn。tol和ε表示公差。运算符(·)+表示伪逆；
S4.5：判断t是否大于Time，若是，则结束计算，若不是，令t＝t+h，返回步骤S4.3。
步骤S4.1中所述初始参数包括p0,
Figure BDA0003064799310000096
Time,h，其中p0表示初始时刻的坐标，
Figure BDA0003064799310000097
表示初始时刻的速度，Time表示模拟时间，h表示时间步长。
为了进一步说明本发明的可行性和有效性，选用了一个曲柄滑块系统，具体如图3所示。在本例中，曲柄和滑块被设置为刚性，连杆被视为柔性。连杆机构由3088个低阶三角形单元组成，自由度为6060。连杆的厚度、材料密度、杨氏模量和泊松比分别设置为0.01m、2700kg/m3、70Gpa和0.3。链接的参考条件设置为固定。曲柄和滑块的z轴质量和惯性矩分别设置为0.475kg，2.8169E-4kg.m2和0.290kg，2.31889E-4kg.m2。总仿真时间和时间步长分别设置为0.2s和1E-4s，联动稀疏度设置为120，采样次数设置为400。并将扭矩M＝15Nm施加到曲柄上。在初始配置中，曲柄和连杆沿x方向放置。
分别计算FOM算法和GGN算法。C点的x位移如图4所示。FOM算法的结果和GGN的算法高度接近。如图5所示，FOM算法与GGN算法之间的相对误差。其相对误差的大小为1E-4，这也意味着该方法的结果是准确的。如图6所示，所GGN算法自适应地选择了120个模态坐标的最大绝对值。由计算结果可知，FOM算法的计算时间为37725秒，而GGN算法的计算时间仅为18786秒，使用该方法计算时间几乎减少了50％。
实施例2
本实施例提供一种基于自适应模态的柔性多体系统动态响应计算系统，如图7，所述系统基于实施例1所述的基于自适应模态的柔性多体系统动态响应计算方法，包括：
模态模块，所述模态模块将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
抽样模块，所述抽样模块设计柔性多体系统的运动方程的抽样矩阵；
范数优化模块，所述范数优化模块求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
GGN求解模块，所述GGN求解模块利用GGN算法求解l1范数优化问题，得到柔性多体系统动态响应。
相同或相似的标号对应相同或相似的部件；
附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制；
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (10)

1.一种基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，包括以下步骤：
S1：将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
S2：设计柔性多体系统的运动方程的抽样矩阵；
S3：求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
S4：利用贪婪高斯-牛顿算法求解l1范数优化问题，得到柔性多体系统动态响应。
2.根据权利要求1所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S1中将柔性多体系统的模态坐标定义为正交振型上的稀疏系数，具体为：
设模态坐标是稀疏的或近似稀疏的，柔性多体系统中有N个柔性体，则广义坐标q表示如下：
Figure FDA0003064799300000011
式中，i表示第i个柔性体柔性体，
Figure FDA0003064799300000012
Figure FDA0003064799300000013
为柔性体i的弹性坐标的向量，RiT为柔性体i相对于惯性系的位置向量，θiT为柔性体i相对于惯性系的角位移，I是单位矩阵，Bi是由柔性体i的全模态形状组成，
Figure FDA0003064799300000014
Figure FDA0003064799300000015
其中
Figure FDA0003064799300000016
表示柔性多体系统的弹性坐标，
Figure FDA0003064799300000017
表示柔性多体系统的模态坐标；
结合拉格朗日乘子向量λ，可得：
Figure FDA0003064799300000021
式中，
Figure FDA0003064799300000022
表示信号，
Figure FDA0003064799300000023
表示基，
Figure FDA0003064799300000024
表示系数向量。
3.根据权利要求2所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S2中柔性多体系统的运动方程的抽样矩阵，具体为：
Figure FDA0003064799300000025
式中，
Figure FDA0003064799300000026
表示相对于柔性体的位置和旋转的质量矩阵，
Figure FDA0003064799300000027
是耦合刚性运动和变形的质量矩阵，
Figure FDA0003064799300000028
表示单个柔性体与弹性坐标有关的质量矩阵，
Figure FDA0003064799300000029
表示与单个柔性体的弹性坐标有关的对称刚度矩阵，
Figure FDA00030647993000000210
和
Figure FDA00030647993000000211
分别是相对于刚性和弹性坐标的外力，
Figure FDA00030647993000000212
和
Figure FDA00030647993000000213
分别是关于刚性和弹性坐标的二次速度矢量，其中Gr和Gf表示与刚性和弹性坐标相关的g(q)的Jacobin矩阵，Ωi表示剩余方程的指数，I(Ωi,:)表示提取单元矩阵的Ωi行。
4.根据权利要求3所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S2中柔性多体系统的运动方程，简化写为：
Figure FDA00030647993000000214
式中，(·)Ω表示柔性体ith采样后的矩阵或向量，
Figure FDA00030647993000000215
表示柔性多体系统的质量矩阵，
Figure FDA00030647993000000216
表示柔性多体系统的刚度矩阵，
Figure FDA00030647993000000217
表示柔性多体系统的约束雅可比矩阵转置形式，
Figure FDA00030647993000000218
表示柔性多体系统的广义外力，
Figure FDA00030647993000000219
表示柔性多体系统的二次速度矢量。
5.根据权利要求4所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
6.根据权利要求5所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S3中每个时间步长下柔性多体系统的动态响应，非线性映射算子定义为：
Figure FDA0003064799300000031
式中，Φ∈Rm×n是一个测量矩阵，Φ(·)∈Rm×l表示非线性映射算子，
Figure FDA0003064799300000032
表示t时刻的弹性坐标和模态坐标，h表示步长，M(pt)表示t时刻质量矩阵，
Figure FDA0003064799300000033
表示t时刻柔性多体系统采样后的矩阵或向量，f(pt)表示t时刻作用于柔性体的广义外力，GT(pt)表示t时刻约束雅可比矩阵的转置，
Figure FDA0003064799300000034
表示t时刻拉格朗日乘子，g(pt)表示t时刻约束方程，且
Figure FDA0003064799300000035
a表示系数向量。
7.根据权利要求6所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S3中每个时间步长动态响应的l1范数优化问题，定义为：
Figure FDA0003064799300000036
式中，
Figure FDA0003064799300000038
表示t时刻柔性多体系统的模态坐标，||·||1表示1范数，s表示柔性多体系统模态坐标的稀疏度，则
Figure FDA0003064799300000039
的稀疏性可表示为:
Figure FDA0003064799300000037
式中，si表示柔性体i模态坐标的稀疏度。
8.根据权利要求7所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S4中利用GGN算法求解l1范数优化问题，具体为：
S4.1：读取每个柔性体的有限元信息，计算每个柔性体的惯性形状积分模态，设置初始参数，所述惯性形状积分模态包括mff和Kff；
S4.2：求解
Figure FDA0003064799300000041
得到初始时刻的加速度
Figure FDA0003064799300000042
和拉格朗日乘子λ0，式中，M表示柔性多体系统的质量矩阵，G表示柔性多体系统的约束雅可比矩阵，
Figure FDA0003064799300000043
表示柔性多体系统的加速度，λ表示柔性多体系统的拉格朗日乘子，
Figure FDA0003064799300000044
表示作用于柔性多体系统的广义外力，
Figure FDA0003064799300000045
表示与加速度有关的向量；
S4.3：对每个柔性体进行测量，得到其待定运动方程Φ(c)；
S4.4：利用GGN算法求解
Figure FDA0003064799300000046
得到p1,
Figure FDA0003064799300000047
表示在t时刻的坐标，速度，加速度和拉格朗日乘子；
S4.5：判断t是否大于Time，若是，则结束计算，若不是，令t＝t+h，返回步骤S4.3。
9.根据权利要求8所述的基于自适应模态的柔性多体系统动态响应计算方法，其特征在于，步骤S4.1中所述初始参数包括p0,
Figure FDA0003064799300000048
Time,h，其中p0表示初始时刻的坐标，
Figure FDA0003064799300000049
表示初始时刻的速度，Time表示模拟时间，h表示时间步长。
10.一种基于自适应模态的柔性多体系统动态响应计算系统，其特征在于，所述系统基于所述权利要求1至9任一项所述的基于自适应模态的柔性多体系统动态响应计算方法，包括：
模态模块，所述模态模块将柔性多体系统的模态坐标定义为正交振型上的稀疏系数；
抽样模块，所述抽样模块设计柔性多体系统的运动方程的抽样矩阵；
范数优化模块，所述范数优化模块求解每个时间步长下柔性多体系统的动态响应，得到每个时间步长动态响应的l1范数优化问题；
GGN求解模块，所述GGN求解模块利用GGN算法求解l1范数优化问题，得到柔性多体系统动态响应。
CN202110523265.XA 2021-05-13 2021-05-13 基于自适应模态的柔性多体系统动态响应计算方法和系统 Active CN113297730B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110523265.XA CN113297730B (zh) 2021-05-13 2021-05-13 基于自适应模态的柔性多体系统动态响应计算方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110523265.XA CN113297730B (zh) 2021-05-13 2021-05-13 基于自适应模态的柔性多体系统动态响应计算方法和系统

Publications (2)

Publication Number Publication Date
CN113297730A true CN113297730A (zh) 2021-08-24
CN113297730B CN113297730B (zh) 2023-04-11

Family

ID=77321843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110523265.XA Active CN113297730B (zh) 2021-05-13 2021-05-13 基于自适应模态的柔性多体系统动态响应计算方法和系统

Country Status (1)

Country Link
CN (1) CN113297730B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114357651A (zh) * 2022-01-04 2022-04-15 北京理工大学 一种高效的不确定性柔性多体系统动力学响应预测方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060139347A1 (en) * 2004-12-27 2006-06-29 Choi Min G Method and system of real-time graphical simulation of large rotational deformation and manipulation using modal warping
CN104915498A (zh) * 2015-06-08 2015-09-16 广东工业大学 基于模型识别与等效简化的高速平台运动参数自整定方法
US20160350462A1 (en) * 2014-06-10 2016-12-01 Guangdong University Of Technology Method of planning asymmetric variable acceleration based on non-linear finite element dynamic response simulation
CN107122515A (zh) * 2017-03-17 2017-09-01 北京航空航天大学 基于绝对节点坐标法的绳系运输系统的动力学分析方法
CN107220421A (zh) * 2017-05-18 2017-09-29 北京理工大学 一种空间复杂柔性结构多体系统动力学建模与计算方法
CN107545126A (zh) * 2017-09-28 2018-01-05 大连理工大学 一种基于多体系统滑移绳索单元的聚合式张拉整体结构动力响应分析方法
US20180307786A1 (en) * 2017-04-23 2018-10-25 Cmlabs Simulations Inc. System and method for executing a simulation of a constrained multi-body system
CN109543264A (zh) * 2018-11-12 2019-03-29 天津理工大学 一种基于多维度重构校正的柔性多体机器人建模与求解方法
US20190286786A1 (en) * 2018-03-16 2019-09-19 Autodesk, Inc. Efficient sensitivity analysis for generative parametric design of dynamic mechanical assemblies
CN111159636A (zh) * 2019-12-04 2020-05-15 大连理工大学 基于绝对节点坐标描述的柔性多体系统动力学半解析灵敏度分析方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060139347A1 (en) * 2004-12-27 2006-06-29 Choi Min G Method and system of real-time graphical simulation of large rotational deformation and manipulation using modal warping
US20160350462A1 (en) * 2014-06-10 2016-12-01 Guangdong University Of Technology Method of planning asymmetric variable acceleration based on non-linear finite element dynamic response simulation
CN104915498A (zh) * 2015-06-08 2015-09-16 广东工业大学 基于模型识别与等效简化的高速平台运动参数自整定方法
US20170124249A1 (en) * 2015-06-08 2017-05-04 Guangdong University Of Technology High-speed platform motion parameter self-tuning method based on model identification and equivalent simplification
CN107122515A (zh) * 2017-03-17 2017-09-01 北京航空航天大学 基于绝对节点坐标法的绳系运输系统的动力学分析方法
US20180307786A1 (en) * 2017-04-23 2018-10-25 Cmlabs Simulations Inc. System and method for executing a simulation of a constrained multi-body system
CN107220421A (zh) * 2017-05-18 2017-09-29 北京理工大学 一种空间复杂柔性结构多体系统动力学建模与计算方法
CN107545126A (zh) * 2017-09-28 2018-01-05 大连理工大学 一种基于多体系统滑移绳索单元的聚合式张拉整体结构动力响应分析方法
US20190286786A1 (en) * 2018-03-16 2019-09-19 Autodesk, Inc. Efficient sensitivity analysis for generative parametric design of dynamic mechanical assemblies
CN109543264A (zh) * 2018-11-12 2019-03-29 天津理工大学 一种基于多维度重构校正的柔性多体机器人建模与求解方法
CN111159636A (zh) * 2019-12-04 2020-05-15 大连理工大学 基于绝对节点坐标描述的柔性多体系统动力学半解析灵敏度分析方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
JIANG, QS ET AL.: "Characteristic Studies of Coupled Rigid and Flexible Multi-body Dynamic for Quick Moving Transfer Robot with Big Inertia", 《PROCEEDINGS OF 2013 INTERNATIONAL FORUM ON SPECIAL EQUIPMENTS AND ENGINEERING MECHANICS》 *
WANG BO-XING ET AL.: "Research on template-based vehicle multi-body dynamics modeling", 《COMPUTER ENGINEERING AND DESIGN》 *
冯诚: "快速柔性多体动力学分析及其在机械臂优化中的应用", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
刘延柱 等: "《多体系统动力学》", 30 December 2014, 高等教育出版社 *
刘树青等: "基于零空间法的柔性多体系统动力学计算及仿真", 《东南大学学报(自然科学版)》 *
孙世基 等: "内燃机曲柄连杆机构运动动力学分析的柔性多体建模方法", 《设计研究》 *
洪嘉振: "《计算多体系统动力学》", 30 December 1999, 北京高等教育出版社 *
葛翔赟 等: "铰接连接结构动态响应与强度计算方法", 《机械设计与制造》 *
黄鹏程 等: "基于注意力机-多任务网络的液压系统多状态监测识别方法研究", 《机电工程》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114357651A (zh) * 2022-01-04 2022-04-15 北京理工大学 一种高效的不确定性柔性多体系统动力学响应预测方法

Also Published As

Publication number Publication date
CN113297730B (zh) 2023-04-11

Similar Documents

Publication Publication Date Title
CN113111430B (zh) 基于非线性气动力降阶的弹性飞机飞行动力学建模方法
Long et al. A meshless local Petrov-Galerkin method for solving the bending problem of a thin plate
CN108986220B (zh) 一种加速有限元求解物体网格模型弹性变形的方法
CN107589671B (zh) 一种基于事件驱动的卫星姿态控制方法
Svendsen et al. Resonant vibration control of rotating beams
CN107633132B (zh) 空间分布动态载荷的等效识别方法
CN106354954B (zh) 一种基于叠层基函数的三维力学模态仿真模拟方法
Ruiz et al. Eigenvector sensitivity when tracking modes with repeated eigenvalues
CN111159636A (zh) 基于绝对节点坐标描述的柔性多体系统动力学半解析灵敏度分析方法
CN110457785B (zh) 一种用于结构大变形响应的物质点法的物质信息映射方法
CN110287505B (zh) 飞行器稳定性分析方法
CN113297730A (zh) 基于自适应模态的柔性多体系统动态响应计算方法和系统
CN115157250A (zh) 一种七自由度机械臂动力学参数辨识方法
CN109902418B (zh) 一种计算eacld中心刚体-悬臂梁模型动力学响应的仿真方法
CN110147571B (zh) 一种组件结构的拓扑优化方法及装置
CN103235515A (zh) 一种利用零运动避免单框架控制力矩陀螺群框架轴转速死区的方法
CN112859948B (zh) 基于分布扰动的柔性杆振动抑制方法及系统
JP4584905B2 (ja) 車両用シミュレーションモデル作成装置
CN113505405A (zh) 等效荷载获取方法、基于等效荷载的拓扑优化方法及系统
CN116611295A (zh) 一种海上风力机叶片动力响应计算方法及系统
JPH0921720A (ja) 構造振動解析方法
CN115533911A (zh) 基于块脉冲函数的柔性机械臂分数阶建模与补偿方法
Trier et al. Design sensitivities by the adjoint variable method in nonlinear structural dynamics
Bicchi et al. Optimal design of dynamic multi-axis force/torque sensor
CN113836621A (zh) 一种基于无迹变换的结构状态/参数/载荷联合识别方法

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