CN110188417A - 一种基于多级超单元的复杂结构有限元模型修正方法 - Google Patents

一种基于多级超单元的复杂结构有限元模型修正方法 Download PDF

Info

Publication number
CN110188417A
CN110188417A CN201910397001.7A CN201910397001A CN110188417A CN 110188417 A CN110188417 A CN 110188417A CN 201910397001 A CN201910397001 A CN 201910397001A CN 110188417 A CN110188417 A CN 110188417A
Authority
CN
China
Prior art keywords
matrix
hyperelement
finite element
multistage
boundary
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
CN201910397001.7A
Other languages
English (en)
Other versions
CN110188417B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201910397001.7A priority Critical patent/CN110188417B/zh
Publication of CN110188417A publication Critical patent/CN110188417A/zh
Application granted granted Critical
Publication of CN110188417B publication Critical patent/CN110188417B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • 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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • 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

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于多级超单元的复杂结构有限元模型修正方法，包括：建立初始有限元模型，将初始有限元模型划分成多级超单元有限元模型；对多级超单元模型采用模态综合法进行缩聚，将超单元缩聚成质量矩阵和刚度矩阵，将所得的质量矩阵和刚度矩阵装配到残余结构，对残余结构进行模态分析得到模态振型和固有频率；测量实际结构的试验模态，采用摄动法构造试验固有频率和模态分析得到的固有频率的残差，建立目标函数，通过优化反问题求解，实现模型修正。本发明在初始有限元模型的基础上，先进行多级超单元建模分析，选择合适结构作为残余结构，对残余结构进行模态分析；对残余结构进行模型修正，有效地提高了有限元模型修正的效率。

Description

一种基于多级超单元的复杂结构有限元模型修正方法
技术领域
本发明涉及有限元模型修正方法，尤其涉及一种基于多级超单元的复杂结构有限元模型修正方法。
背景技术
对于传动系统这类的工程结构，常常需要快速准确地分析整个结构的动态特性，因此迫切需要提高计算精度与效率，而数学建模分析是结构动力分析中最关键的一环。由于工程结构的日益大型复杂化和建模过程中的一些因素，如单元类型、材料属性和几何形状等，往往存在一定的误差，影响分析结果的准确性。因此，建立复杂结构的精细化的有限元模型并进行准确的计算分析，这本身就是一件非常困难的工作。
因此，亟待解决上述问题。
发明内容
发明目的：针对大型复杂结构模型规模较大导致修正时间过长、运算困难的问题，本发明提供一种基于多级超单元的复杂结构有限元模型修正方法，在多级超单元建模的基础上进行模型修正，有效地提高了有限元模型修正的效率。
技术方案：为实现以上目的，本发明公开了一种基于多级超单元的复杂结构有限元模型修正方法，包括以下步骤：
(1)建立初始有限元模型，将初始有限元模型划分成多级超单元有限元模型；
(2)对多级超单元模型采用模态综合法进行缩聚，将超单元缩聚成质量矩阵和刚度矩阵，将所得的质量矩阵和刚度矩阵装配到残余结构，对残余结构进行模态分析得到模态振型和固有频率；
(3)测量实际结构的试验模态，采用摄动法构造试验固有频率和步骤(2)中模态分析得到的固有频率的残差，建立目标函数，通过优化反问题求解，实现基于多级超单元的复杂结构有限元模型修正。
其中，所述步骤(1)中具体划分方法包括如下步骤：
(1.1)建立初始有限元模型；
(1.2)考虑到研究对象本身的材料分布特性和连接属性，对初始有限元模型进行分级，得到各级超单元和残余结构，具体分级方法包括：首先将初始有限元模型中具有装配和分层关系的结构划分为不同级的超单元，其次，将同一级超单元中具有不同材料属性的结构划分为同一级下的不同超单元，最后将变形量可能大于最小变形量50％以上的结构作为残余结构。
优选的，所述步骤(2)中模态综合法的具体步骤以下步骤：
(2.1)超单元无阻尼运动方程为
式中
其中M表示超单元的刚度矩阵，K表示超单元的质量矩阵，表示超单元的加速度，u表示位移，f表示外力，下标ii、ij(ji)和jj分别表示内部自由度、边界和内部耦合自由度、约束边界自由度；
以模态坐标变换矩阵将物理坐标变换到缩减的模态坐标p下
式中
其中分别表示模态坐标下的质量矩阵和刚度矩阵，表示模态坐标二次求导，Φc0表示全部界面自由度的约束模态，φbL表示约束界面主模态，表示约束界面特征向量；Ij表示单位矩阵；
(2.2)根据超单元的界面位移协调和力平衡条件，将多级超单元和残余结构进行模态综合，建立整体结构的无阻尼自由振动方程
式中
其中，表示整体结构的刚度矩阵，表示整体结构的质量矩阵，Mres是残余结构的刚度矩阵，Kres是残余结构的质量矩阵；q表示广义坐标，表示广义坐标二次求导，I表示单位矩阵，和分别表示超单元A刚度矩阵中的内部矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A中刚度矩阵的约束边界矩阵、超单元B中刚度矩阵的内部矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的约束边界矩阵， 和分别表示超单元A质量矩阵中的内部矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的约束边界矩阵、超单元B质量矩阵中的内部矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的约束边界矩阵；
对整体结构的无阻尼自由振动方程进行求解，即可得到整体结构的模态振型和固有频率。
再者，所述步骤(3)中具体方法包括以下步骤：
(3.1)若fe、fa分别表示试验特征值向量和分析特征值向量，求解特征值问题，进一步将结构模型修正问题转化为优化问题，其修正目标函数为
式中：ε＝fe-fa(p)是试验和分析得到的特征值残差向量，Wε是对称加权矩阵，p和分别表示待修正参数的上、下边界值；
(3.2)有限元模型修正过程中，通过调整结构设计参数p迭代修正目标函数公式，直至满足迭代收敛条件。
进一步，所述步骤(3.2)中迭代收敛条件为
有益效果：与现有技术相比，本发明具有以下显著优点：本发明基于多级超单元的复杂结构有限元模型修正方法考虑了大型复杂结构模型修正的困难，在初始有限元模型的基础上，先进行多级超单元建模分析，选择合适结构作为残余结构，对残余结构进行模态分析；只需要对残余结构进行模型修正，有效地提高了有限元模型修正的效率。
附图说明
图1为本发明的流程示意图；
图2为本发明实施例中的行星轮系和支撑轴的精细化有限元模型；
图3为本发明实施例中的行星轮系和支撑轴多级超单元划分示意图；
图4为本发明基于多级超单元的行星轮系和支撑轴有限元模型修正的残余结构频率误差收敛曲线；
图5为本发明基于多级超单元的行星轮系和支撑轴有限元模型修正的残余结构参数误差收敛曲线。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示，本发明一种基于多级超单元的复杂结构有限元模型修正方法，包括以下步骤：
1、一种基于多级超单元的复杂结构有限元模型修正方法，其特征在于，包括以下步骤：
(1)建立初始有限元模型，将初始有限元模型划分成多级超单元有限元模型；
其中具体划分方法包括如下步骤：
(1.1)建立初始有限元模型；
(1.2)考虑到研究对象本身的材料分布特性和连接属性，对初始有限元模型进行分级，得到各级超单元和残余结构，具体分级方法包括：首先将初始有限元模型中具有装配和分层关系的结构划分为不同级的超单元，其次，将同一级超单元中具有不同材料属性的结构划分为同一级下的不同超单元，最后将变形量可能大于最小变形量50％以上的结构作为残余结构；
(2)对多级超单元模型采用模态综合法进行缩聚，将超单元缩聚成质量矩阵和刚度矩阵，将所得的质量矩阵和刚度矩阵装配到残余结构，对残余结构进行模态分析得到模态振型和固有频率；
其中模态综合法的具体步骤以下步骤：
(2.1)超单元无阻尼运动方程为
式中
其中M表示超单元的刚度矩阵，K表示超单元的质量矩阵，表示超单元的加速度，u表示位移，f表示外力，下标ii、ij(ji)和jj分别表示内部自由度、边界和内部耦合自由度、约束边界自由度；
以模态坐标变换矩阵将物理坐标变换到缩减的模态坐标p下
式中
其中分别表示模态坐标下的质量矩阵和刚度矩阵，表示模态坐标二次求导，Φc0表示全部界面自由度的约束模态，φbL表示约束界面主模态，表示约束界面特征向量；Ij表示单位矩阵；
(2.2)根据超单元的界面位移协调和力平衡条件，将多级超单元和残余结构进行模态综合，建立整体结构的无阻尼自由振动方程
式中
其中，表示整体结构的刚度矩阵，表示整体结构的质量矩阵，Mres是残余结构的刚度矩阵，Kres是残余结构的质量矩阵；q表示广义坐标，表示广义坐标二次求导，I表示单位矩阵，和分别表示超单元A刚度矩阵中的内部矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A中刚度矩阵的约束边界矩阵、超单元B中刚度矩阵的内部矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的约束边界矩阵， 和分别表示超单元A质量矩阵中的内部矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的约束边界矩阵、超单元B质量矩阵中的内部矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的约束边界矩阵；
对整体结构的无阻尼自由振动方程进行求解，即可得到整体结构的模态振型和固有频率；
(3)测量实际结构的试验模态，采用摄动法构造试验固有频率和步骤(2)中模态分析得到的固有频率的残差，建立目标函数，通过优化反问题求解，实现基于多级超单元的复杂结构有限元模型修正；
其中具体方法包括以下步骤：
(3.1)若fe、fa分别表示试验特征值向量和分析特征值向量，求解特征值问题，进一步将结构模型修正问题转化为优化问题，其修正目标函数为
式中：ε＝fe-fa(p)是试验和分析得到的特征值残差向量，Wε是对称加权矩阵，p和分别表示待修正参数的上、下边界值；
(3.2)有限元模型修正过程中，通过调整结构设计参数p迭代修正目标函数公式，直至满足迭代收敛条件，其中迭代收敛条件为
实施例
一种基于多级超单元的复杂结构有限元模型修正方法，包括以下步骤：
(1)建立初始有限元模型，将初始有限元模型划分成多级超单元有限元模型；具体划分步骤包括：
(1.1)以行星轮系和支撑轴为研究对象，包括一个太阳轮、四个行星轮、一个内齿圈和一个支撑轴；行星轮系的几何参数如表1所示，其整体结构为钢材料，弹性模量E＝190000MPa，泊松比μ＝0.3，密度ρ＝7500kg/m3；在商用有限元软件Hypermesh采用六面体单元进行建模，图2所示为该结构的有限元模型；共建立57492个节点，45485个单元；定义全模型为包含所有计算自由度的有限元模型；
表1行星轮系几何参数
(1.2)考虑齿轮之间的耦合，选取太阳轮、行星轮和内齿圈的齿作为残余结构，各个齿轮的轮盘和支撑轴作为超单元，将支撑轴作为一级超单元，所有齿轮的轮盘作为二级超单元；
(2)对多级超单元模型采用模态综合法进行缩聚，将超单元缩聚成质量矩阵和刚度矩阵，将所得的质量矩阵和刚度矩阵装配到残余结构，对残余结构进行模态分析得到模态振型和固有频率；
其中模态综合法的具体步骤以下步骤：
(2.1)超单元无阻尼运动方程为
式中
其中M表示超单元的刚度矩阵，K表示超单元的质量矩阵，表示超单元的加速度，u表示位移，f表示外力，下标ii、ij(ji)和jj分别表示内部自由度、边界和内部耦合自由度、约束边界自由度；
以模态坐标变换矩阵将物理坐标变换到缩减的模态坐标p下
式中
其中分别表示模态坐标下的质量矩阵和刚度矩阵，表示模态坐标二次求导，Φc0表示全部界面自由度的约束模态，φbL表示约束界面主模态，表示约束界面特征向量；Ij表示单位矩阵；
(2.2)根据超单元的界面位移协调和力平衡条件，将多级超单元和残余结构进行模态综合，建立整体结构的无阻尼自由振动方程
式中
其中，表示整体结构的刚度矩阵，表示整体结构的质量矩阵，Mres是残余结构的刚度矩阵，Kres是残余结构的质量矩阵；q表示广义坐标，表示广义坐标二次求导，I表示单位矩阵，和分别表示超单元A刚度矩阵中的内部矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A中刚度矩阵的约束边界矩阵、超单元B中刚度矩阵的内部矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的约束边界矩阵， 和分别表示超单元A质量矩阵中的内部矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的约束边界矩阵、超单元B质量矩阵中的内部矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的约束边界矩阵；
对整体结构的无阻尼自由振动方程进行求解，即可得到整体结构的模态振型和固有频率；
(3)测量实际结构的试验模态，构造试验固有频率和步骤(2)中模态分析得到的固有频率的残差，建立目标函数，通过优化反问题求解，实现基于多级超单元的复杂结构有限元模型修正；
其中具体方法包括以下步骤：
(3.1)若fe、fa分别表示试验特征值向量和分析特征值向量，求解特征值问题，进一步将结构模型修正问题转化为优化问题，其修正目标函数为
式中：ε＝fe-fa(p)是试验和分析得到的特征值残差向量，Wε是对称加权矩阵，p和分别表示待修正参数的上、下边界值；
(3.2)有限元模型修正过程中，通过调整结构设计参数p迭代修正目标函数公式，直至满足迭代收敛条件，其中迭代收敛条件为
基于多级超单元的复杂结构有限元模型修正方法的验证，结果如下：为比较有限元模型和多级超单元模型的分析精度，对两种模型进行自由模态分析，如表2和表3所示；从表2和表3可以看出，有限元模型分析得到的结构模态与多级超单元分析所得残余结构的模态振型一致；定义频率分析误差ef＝100％×(fw-fs)/fw，其中w和s分别表示有限元模型和多级超单元模型；分析表的频率误差可知，该复杂模型的最大误差出现在第5阶，为0.00253％，其余4阶均小于0.0025％，计算偏差较小。
表2有限元模型和多级超单元模型对应模态振型比较
表3有限元模型和多级超单元模型固有频率误差分析
基于行星轮系和支撑轴的多级超单元有限元模型修正方法，采用摄动法，构造误差模型。假设残余结构的弹性模量E和泊松比μ分别有20％和-20％的初始误差，即E0＝228GPa和μ0＝0.24，选取结构第七阶到第十一阶特征频率，结构前六阶模态为刚体模态，以E和μ为待修正参数进行修正。从图4和图5可以看出，随着迭代次数的增加，目标函数逐渐收敛，而各阶固有频率和结构参数逐渐趋向于修正值。从表2列举的数据可以看出，行星齿轮系结构模型修正后的固有频率与目标频率误差很小，很好地说明了方法的有效性。
表4行星轮系和支撑轴结构计算频率修正前后误差分析

Claims (5)

1.一种基于多级超单元的复杂结构有限元模型修正方法，其特征在于，包括以下步骤：
(1)建立初始有限元模型，将初始有限元模型划分成多级超单元有限元模型；
(2)对多级超单元模型采用模态综合法进行缩聚，将超单元缩聚成质量矩阵和刚度矩阵，将所得的质量矩阵和刚度矩阵装配到残余结构，对残余结构进行模态分析得到模态振型和固有频率；
(3)测量实际结构的试验模态，采用摄动法构造试验固有频率和步骤(2)中模态分析得到的固有频率的残差，建立目标函数，通过优化反问题求解，实现基于多级超单元的复杂结构有限元模型修正。
2.根据权利要求1所述的基于多级超单元的复杂结构有限元模型修正方法，其特征在于，所述步骤(1)中具体划分方法包括如下步骤：
(1.1)建立初始有限元模型；
(1.2)考虑到研究对象本身的材料分布特性和连接属性，对初始有限元模型进行分级，得到各级超单元和残余结构，具体分级方法包括：首先将初始有限元模型中具有装配和分层关系的结构划分为不同级的超单元，其次，将同一级超单元中具有不同材料属性的结构划分为同一级下的不同超单元，最后将变形量可能大于最小变形量50％以上的结构作为残余结构。
3.根据权利要求1所述的基于多级超单元的复杂结构有限元模型修正方法，其特征在于，所述步骤(2)中模态综合法的具体步骤以下步骤：
(2.1)超单元无阻尼运动方程为
式中
其中M表示超单元的刚度矩阵，K表示超单元的质量矩阵，表示超单元的加速度，u表示位移，f表示外力，下标ii、ij(ji)和jj分别表示内部自由度、边界和内部耦合自由度、约束边界自由度；
以模态坐标变换矩阵将物理坐标变换到缩减的模态坐标p下
式中
其中分别表示模态坐标下的质量矩阵和刚度矩阵，表示模态坐标二次求导，Φc0表示全部界面自由度的约束模态，φbL表示约束界面主模态，表示约束界面特征向量；Ij表示单位矩阵；
(2.2)根据超单元的界面位移协调和力平衡条件，将多级超单元和残余结构进行模态综合，建立整体结构的无阻尼自由振动方程
式中
其中，表示整体结构的刚度矩阵，表示整体结构的质量矩阵，Mres是残余结构的刚度矩阵，Kres是残余结构的质量矩阵；q表示广义坐标，表示广义坐标二次求导，I表示单位矩阵，和分别表示超单元A刚度矩阵中的内部矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A刚度矩阵中的边界和内部耦合矩阵、超单元A中刚度矩阵的约束边界矩阵、超单元B中刚度矩阵的内部矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的边界和内部耦合矩阵、超单元B刚度矩阵中的约束边界矩阵， 和分别表示超单元A质量矩阵中的内部矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的边界和内部耦合矩阵、超单元A质量矩阵中的约束边界矩阵、超单元B质量矩阵中的内部矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的边界和内部耦合矩阵、超单元B质量矩阵中的约束边界矩阵；
对整体结构的无阻尼自由振动方程进行求解，即可得到整体结构的模态振型和固有频率。
4.根据权利要求1所述的基于多级超单元的复杂结构有限元模型修正方法，其特征在于，所述步骤(3)中具体方法包括以下步骤：
(3.1)若fe、fa分别表示试验特征值向量和分析特征值向量，求解特征值问题，进一步将结构模型修正问题转化为优化问题，其修正目标函数为
式中：ε＝fe-fa(p)是试验和分析得到的特征值残差向量，Wε是对称加权矩阵，p和分别表示待修正参数的上、下边界值；
(3.2)有限元模型修正过程中，通过调整结构设计参数p迭代修正目标函数公式，直至满足迭代收敛条件。
5.根据权利要求4所述的基于多级超单元的复杂结构有限元模型修正方法，其特征在于，所述步骤(3.2)中迭代收敛条件为
CN201910397001.7A 2019-05-14 2019-05-14 一种基于多级超单元的复杂结构有限元模型修正方法 Active CN110188417B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910397001.7A CN110188417B (zh) 2019-05-14 2019-05-14 一种基于多级超单元的复杂结构有限元模型修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910397001.7A CN110188417B (zh) 2019-05-14 2019-05-14 一种基于多级超单元的复杂结构有限元模型修正方法

Publications (2)

Publication Number Publication Date
CN110188417A true CN110188417A (zh) 2019-08-30
CN110188417B CN110188417B (zh) 2023-03-24

Family

ID=67716190

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910397001.7A Active CN110188417B (zh) 2019-05-14 2019-05-14 一种基于多级超单元的复杂结构有限元模型修正方法

Country Status (1)

Country Link
CN (1) CN110188417B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112507585A (zh) * 2020-11-30 2021-03-16 中南大学 基于emd和模型缩聚的动力响应重构方法及系统、存储介质
CN112580239A (zh) * 2020-12-15 2021-03-30 中南大学 一种子结构响应重构方法及系统、存储介质
CN112781717A (zh) * 2020-12-31 2021-05-11 长江勘测规划设计研究有限责任公司 一种复杂建筑结构自振特性测量方法
CN112818580A (zh) * 2021-02-07 2021-05-18 上海机电工程研究所 基于扩充模态矩阵的间隙结构动力学模型降阶方法及系统
CN112949104A (zh) * 2019-11-26 2021-06-11 中国科学院沈阳自动化研究所 一种协作机器人实时模态分析方法
CN113111547A (zh) * 2021-03-15 2021-07-13 南京航空航天大学 基于缩减基的频域有限元模型修正方法
CN114398808A (zh) * 2021-12-10 2022-04-26 航天科工火箭技术有限公司 螺栓连接结合面快速模型修正方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091033A (zh) * 2014-07-25 2014-10-08 哈尔滨工业大学 基于超单元结合虚拟变形法的桥梁静力有限元模型修正方法
CN109241559A (zh) * 2018-08-01 2019-01-18 东南大学 一种基于子结构的复合材料弹性参数识别方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091033A (zh) * 2014-07-25 2014-10-08 哈尔滨工业大学 基于超单元结合虚拟变形法的桥梁静力有限元模型修正方法
CN109241559A (zh) * 2018-08-01 2019-01-18 东南大学 一种基于子结构的复合材料弹性参数识别方法

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949104A (zh) * 2019-11-26 2021-06-11 中国科学院沈阳自动化研究所 一种协作机器人实时模态分析方法
CN112949104B (zh) * 2019-11-26 2023-07-25 中国科学院沈阳自动化研究所 一种协作机器人实时模态分析方法
CN112507585A (zh) * 2020-11-30 2021-03-16 中南大学 基于emd和模型缩聚的动力响应重构方法及系统、存储介质
CN112507585B (zh) * 2020-11-30 2022-05-24 中南大学 基于emd和模型缩聚的动力响应重构方法及系统、存储介质
CN112580239A (zh) * 2020-12-15 2021-03-30 中南大学 一种子结构响应重构方法及系统、存储介质
CN112580239B (zh) * 2020-12-15 2022-07-05 中南大学 一种子结构响应重构方法及系统、存储介质
CN112781717A (zh) * 2020-12-31 2021-05-11 长江勘测规划设计研究有限责任公司 一种复杂建筑结构自振特性测量方法
CN112818580A (zh) * 2021-02-07 2021-05-18 上海机电工程研究所 基于扩充模态矩阵的间隙结构动力学模型降阶方法及系统
CN112818580B (zh) * 2021-02-07 2022-08-16 上海机电工程研究所 基于扩充模态矩阵的间隙结构动力学模型降阶方法及系统
CN113111547A (zh) * 2021-03-15 2021-07-13 南京航空航天大学 基于缩减基的频域有限元模型修正方法
CN114398808A (zh) * 2021-12-10 2022-04-26 航天科工火箭技术有限公司 螺栓连接结合面快速模型修正方法
CN114398808B (zh) * 2021-12-10 2024-03-29 航天科工火箭技术有限公司 螺栓连接结合面快速模型修正方法

Also Published As

Publication number Publication date
CN110188417B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
CN110188417A (zh) 一种基于多级超单元的复杂结构有限元模型修正方法
CN108959686A (zh) 一种基于灵敏度分析的有限元模型修正方法
Chittick et al. An asymmetric suboptimization approach to aerostructural optimization
CN107391891A (zh) 一种基于模型融合方法的大展弦比机翼优化设计方法
CN104679956A (zh) 考虑动态特性的高速压力机底座可靠性稳健设计方法
CN108763658A (zh) 基于等几何方法的组合薄壁结构固有频率设计方法
CN105930562A (zh) 一种非概率条件下的结构性能优化设计方法
CN105678015B (zh) 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法
CN106529035A (zh) 一种基于三维分形理论的栓接结构动态特性分析方法
CN107016173B (zh) 一种高速压力机底座动态特性的可靠性设计方法
CN103995914A (zh) 基于动态特性分析的磨齿机立柱结构优化设计的方法
Zhang et al. A construction method of digital twin model for contact characteristics of assembly interface
CN110516365A (zh) 一种测试螺栓连接刚度的方法
Han et al. Numerical improvement for the mechanical performance of bikes based on an intelligent PSO-ABC algorithm and WSN technology
CN110399675A (zh) 一种基于遗传算法的电梯层门多目标优化设计方法
CN109255141A (zh) 一种汽车车身正向概念设计截面形状优化方法
CN107657073A (zh) 一种三明治复合材料结构的参数识别方法
CN114996995A (zh) 一种超材料隔振单元性能预报方法和系统
CN110334440A (zh) 基于神经网络的转子有限元模型优化设计方法
CN113283143A (zh) 一种基于动静力数据的叠合梁有限元模型修正方法
CN106354944A (zh) 一种基于振动模态的磨机基础边界支承刚度的识别方法
CN103870614B (zh) 一种结构概率优化设计方法
CN111737908A (zh) 一种基于动载荷静力等效的蒙皮桁条结构快速动态优化设计方法
Weule et al. Structural Optimization of Machine Tools including the static and dynamic Workspace Behavior
CN115455754A (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