CN113297730A - 基于自适应模态的柔性多体系统动态响应计算方法和系统 - Google Patents
基于自适应模态的柔性多体系统动态响应计算方法和系统 Download PDFInfo
- 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
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表示如下:
式中,i表示第i个柔性体柔性体, 为柔性体i的弹性坐标的向量,RiT为柔性体i相对于惯性系的位置向量,θiT为柔性体i相对于惯性系的角位移,I是单位矩阵,Bi是由柔性体i的全模态形状组成, 其中表示柔性多体系统的弹性坐标,表示柔性多体系统的模态坐标;
结合拉格朗日乘子向量λ,可得:
优选地,步骤S2中柔性多体系统的运动方程的抽样矩阵,具体为:
式中,表示相对于柔性体的位置和旋转的质量矩阵,是耦合刚性运动和变形的质量矩阵,表示单个柔性体与弹性坐标有关的质量矩阵,表示与单个柔性体的弹性坐标有关的对称刚度矩阵,和分别是相对于刚性和弹性坐标的外力,和分别是关于刚性和弹性坐标的二次速度矢量,其中Gr和Gf表示与刚性和弹性坐标相关的g(q)的Jacobin矩阵,Ωi表示剩余方程的指数,I(Ωi,:)表示提取单元矩阵的Ωi行。
优选地,步骤S2中柔性多体系统的运动方程,简化写为:
式中,(·)Ω表示柔性体ith采样后的矩阵或向量,表示柔性多体系统的质量矩阵,表示柔性多体系统的刚度矩阵,表示柔性多体系统的约束雅可比矩阵转置形式,表示柔性多体系统的广义外力,表示柔性多体系统的二次速度矢量。
优选地,步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
优选地,步骤S3中每个时间步长下柔性多体系统的动态响应,非线性映射算子定义为:
式中,Φ∈Rm×n是一个测量矩阵,Φ(·)∈Rm×l表示非线性映射算子,表示t时刻的弹性坐标和模态坐标,h表示步长,M(pt)表示t时刻质量矩阵,(·)Ωt表示t时刻柔性多体系统采样后的矩阵或向量,f(pt)表示t时刻作用于柔性体的广义外力,GT(pt)表示t时刻约束雅可比矩阵的转置,表示t时刻拉格朗日乘子,g(pt)表示t时刻约束方程,且a表示系数向量。
优选地,步骤S3中每个时间步长动态响应的l1范数优化问题,定义为:
式中,si表示柔性体i模态坐标的稀疏度。
优选地,步骤S4中利用GGN算法求解l1范数优化问题,具体为:
S4.1:读取每个柔性体的有限元信息,计算每个柔性体的惯性形状积分模态,设置初始参数,所述惯性形状积分模态包括mff和Kff;
S4.2:求解得到初始时刻的加速度和拉格朗日乘子λ0,式中,M表示柔性多体系统的质量矩阵,G表示柔性多体系统的约束雅可比矩阵,表示柔性多体系统的加速度,λ表示柔性多体系统的拉格朗日乘子,表示作用于柔性多体系统的广义外力,表示与加速度有关的向量;
S4.3:对每个柔性体进行测量,得到其待定运动方程Φ(c);
S4.5:判断t是否大于Time,若是,则结束计算,若不是,令t=t+h,返回步骤S4.3。
一种基于自适应模态的柔性多体系统动态响应计算系统,所述系统基于上述所述的基于自适应模态的柔性多体系统动态响应计算方法,包括:
模态模块,所述模态模块将柔性多体系统的模态坐标定义为正交振型上的稀疏系数;
抽样模块,所述抽样模块设计柔性多体系统的运动方程的抽样矩阵;
范数优化模块,所述范数优化模块求解每个时间步长下柔性多体系统的动态响应,得到每个时间步长动态响应的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表示如下:
式中,i表示第i个柔性体柔性体,为柔性体i的弹性坐标的向量,RiT为柔性体i相对于惯性系的位置向量,θiT为柔性体i相对于惯性系的角位移,I是单位矩阵,Bi是由柔性体i的全模态形状组成, 其中表示柔性多体系统的弹性坐标,表示柔性多体系统的模态坐标;
结合拉格朗日乘子向量λ,可得:
步骤S2中柔性多体系统的运动方程的抽样矩阵,具体为:
式中,表示相对于柔性体的位置和旋转的质量矩阵,是耦合刚性运动和变形的质量矩阵,表示单个柔性体与弹性坐标有关的质量矩阵,表示与单个柔性体的弹性坐标有关的对称刚度矩阵,和分别是相对于刚性和弹性坐标的外力,和分别是关于刚性和弹性坐标的二次速度矢量,其中Gr和Gf表示与刚性和弹性坐标相关的g(q)的Jacobin矩阵,Ωi表示剩余方程的指数,I(Ωi,:)表示提取单元矩阵的Ωi行。
步骤S2中柔性多体系统的运动方程,简化写为:
式中,(·)Ω表示柔性体ith采样后的矩阵或向量,表示柔性多体系统的质量矩阵,表示柔性多体系统的刚度矩阵,表示柔性多体系统的约束雅可比矩阵转置形式,表示柔性多体系统的广义外力,表示柔性多体系统的二次速度矢量。
步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
步骤S3中每个时间步长下柔性多体系统的动态响应,非线性映射算子定义为:
式中,Φ∈Rm×n是一个测量矩阵,Φ(·)∈Rm×l表示非线性映射算子,表示t时刻的弹性坐标和模态坐标,h表示步长,M(pt)表示t时刻质量矩阵,(·)Ωt表示t时刻柔性多体系统采样后的矩阵或向量,f(pt)表示t时刻作用于柔性体的广义外力,GT(pt)表示t时刻约束雅可比矩阵的转置,表示t时刻拉格朗日乘子,g(pt)表示t时刻约束方程,且a表示系数向量。
步骤S3中每个时间步长动态响应的l1范数优化问题,定义为:
式中,si表示柔性体i模态坐标的稀疏度。
步骤S4中利用GGN算法求解l1范数优化问题,具体为:
S4.1:读取每个柔性体的有限元信息,计算每个柔性体的惯性形状积分模态,设置初始参数,所述惯性形状积分模态包括mff和Kff;
S4.2:求解得到初始时刻的加速度和拉格朗日乘子λ0,式中,M表示柔性多体系统的质量矩阵,G表示柔性多体系统的约束雅可比矩阵,表示柔性多体系统的加速度,λ表示柔性多体系统的拉格朗日乘子,表示作用于柔性多体系统的广义外力,表示与加速度有关的向量;
S4.3:对每个柔性体进行测量,得到其待定运动方程Φ(c);
S4.4:利用GGN算法求解得到p1,表示在t时刻的坐标,速度,加速度和拉格朗日乘子,所述利用GGN算法求解的具体过程如图2所示,在算法中,符号Sk表示支持集,包括刚性坐标和拉格朗日乘子的索引,和 表示提取的tth列,由列组成,对应于支持集Sn。tol和ε表示公差。运算符(·)+表示伪逆;
S4.5:判断t是否大于Time,若是,则结束计算,若不是,令t=t+h,返回步骤S4.3。
为了进一步说明本发明的可行性和有效性,选用了一个曲柄滑块系统,具体如图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范数优化问题,得到柔性多体系统动态响应。
5.根据权利要求4所述的基于自适应模态的柔性多体系统动态响应计算方法,其特征在于,步骤S3中每个时间步长下柔性多体系统的动态响应采用一阶后向欧拉法。
8.根据权利要求7所述的基于自适应模态的柔性多体系统动态响应计算方法,其特征在于,步骤S4中利用GGN算法求解l1范数优化问题,具体为:
S4.1:读取每个柔性体的有限元信息,计算每个柔性体的惯性形状积分模态,设置初始参数,所述惯性形状积分模态包括mff和Kff;
S4.2:求解得到初始时刻的加速度和拉格朗日乘子λ0,式中,M表示柔性多体系统的质量矩阵,G表示柔性多体系统的约束雅可比矩阵,表示柔性多体系统的加速度,λ表示柔性多体系统的拉格朗日乘子,表示作用于柔性多体系统的广义外力,表示与加速度有关的向量;
S4.3:对每个柔性体进行测量,得到其待定运动方程Φ(c);
S4.5:判断t是否大于Time,若是,则结束计算,若不是,令t=t+h,返回步骤S4.3。
10.一种基于自适应模态的柔性多体系统动态响应计算系统,其特征在于,所述系统基于所述权利要求1至9任一项所述的基于自适应模态的柔性多体系统动态响应计算方法,包括:
模态模块,所述模态模块将柔性多体系统的模态坐标定义为正交振型上的稀疏系数;
抽样模块,所述抽样模块设计柔性多体系统的运动方程的抽样矩阵;
范数优化模块,所述范数优化模块求解每个时间步长下柔性多体系统的动态响应,得到每个时间步长动态响应的l1范数优化问题;
GGN求解模块,所述GGN求解模块利用GGN算法求解l1范数优化问题,得到柔性多体系统动态响应。
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114357651A (zh) * | 2022-01-04 | 2022-04-15 | 北京理工大学 | 一种高效的不确定性柔性多体系统动力学响应预测方法 |
Citations (10)
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 | 大连理工大学 | 基于绝对节点坐标描述的柔性多体系统动力学半解析灵敏度分析方法 |
-
2021
- 2021-05-13 CN CN202110523265.XA patent/CN113297730B/zh active Active
Patent Citations (11)
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)
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)
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 |